Large Evolutionary Distances Research Articles

Decoding the Functional Interactome of Non-Model Organisms with PHILHARMONIC.

Protein-protein interaction (PPI) networks are a fundamental resource for modeling cellular and molecular function, and a large and sophisticated toolbox has been developed to leverage their structure and topological organization to predict the functional roles of under-studied genes, proteins, and pathways. However, the overwhelming majority of experimentally-determined interactions from which such networks are constructed come from a small number of well-studied model organisms. Indeed, most species lack even a single experimentally-determined interaction in these databases, much less a network to enable the analysis of cellular function, and methods for computational PPI prediction are too noisy to apply directly. We introduce PHILHARMONIC, a novel computational approach that couples deep learning de novo network inference with robust unsupervised spectral clustering algorithms to uncover functional relationships and high-level organization in non-model organisms. Our clustering approach allows us to de-noise the predicted network, producing highly informative functional modules. We also develop a novel algorithm called ReCIPE, which aims to reconnect disconnected clusters, increasing functional enrichment and biological interpretability. We perform remote homology-based functional annotation by leveraging hmmscan and GODomainMiner to assign initial functions to proteins at large evolutionary distances. Our clusters enable us to newly assign functions to uncharacterized proteins through "function by association." We demonstrate the ability of PHILHARMONIC to recover clusters with significant functional coherence in the reef-building coral P. damicornis, its algal symbiont C. goreaui, and the well-annotated fruit fly D. melanogaster. We perform a deeper analysis of the P. damicornis network, where we show that PHILHARMONIC clusters correlate strongly with gene co-expression and investigate several clusters that participate in temperature regulation in the coral, including the first putative functional annotation of several previously uncharacterized proteins. Easy to run end-to-end and requiring only a sequenced proteome, PHILHARMONIC is an engine for biological hypothesis generation and discovery in non-model organisms. PHILHARMONIC is available at https://github.com/samsledje/philharmonic.

Pairwise sequence alignment at arbitrarily large evolutionary distance

Pairwise sequence alignment at arbitrarily large evolutionary distance

Abstract 5750: The mutational landscape of ultraconserved elements in human cancers

Abstract Background: Ultraconserved regions or elements (UCEs) are greater than 100-base pairs in length and are perfectly (100%) conserved across large evolutionary distances in the genomes of at least 3 of 5 placental mammals: human, cow, dog, rat, and mouse. Most UCEs are located in non-coding regions of the genome; however, some UCEs overlap with coding exons. In addition, some UCEs have been found to be transcriptionally active as long non-coding RNAs and are involved in a variety of cellular processes such as cellular proliferation. Although a few UCE sequence variations were reported to be associated with human diseases, the molecular functions of somatic mutations in UCEs remain largely unexplained in human cancers. In this study, we characterized the distribution of somatic UCE mutations throughout a spectrum of cancers and investigated the biological functions of cancer associated UCEs. Methods: We examined somatic UCE mutations in 2,449 cases of 22 cancer types using the PCAWG and ICGC platforms. To validate these mutation patterns, UCE sequencing was performed on the Illumina NovaSeq 6000 platform for MD Anderson patient cohorts. A custom AsCpf1 guide library was built to identify potential regulatory UCE functions in colorectal cancer models, and 3 guides were designed to target 2,247 UCEs. The proliferation ratio of knockout UCEs was analyzed. We then generated stable mutated UCE clones using the CRISPR AsCpf1 technology in colorectal cancer cells. We conducted RNA sequencing (RNA-seq) of mutated UCE cells. The expression of target proteins and genes were analyzed by qRT-PCR and western blotting. To address the effects of the mutated UCE_11311 tumor growth, DLD1 cells were injected in athymic nude mice. Results: We analyzed the WGS data on 2,449 cases of 22 cancer types and identified 24,039 somatic mutations in 10,090 (73.46%) UCEs. These were mostly located in non-coding DNAs, mainly in introns.Based on the RNA-seq data, one of the strongest impacts of mutated UCEs was on ARID1B. We confirmed that specific mutated UCE decreased ARID1B mRNA and protein levels. Our in vivo results demonstrated that UCE mutations enhanced the tumorigenicity of DLD1 xenograft tumors. As a summary, UCE_11311 has a transcriptional enhancer activity on the ARID1B gene and UCE_11311 mutations are actively participating in the tumorigenesis at least in part through regulating ARID1B expression. Discussion: We identified hundreds of unexplored tumorigenic UCEs that need to be further characterized functionally and clinically, and proved one of them to be a transcriptional enhancer of the tumor suppressor ARID1B. Collectively, these data support the concept that certain somatic UCE mutations are frequent and functional in cancer evolution, acting as driver mutations that can be used as new therapeutic targets; other somatic UCE mutations are cancer specific or patient specific, are markers of aggressiveness, and can be used to personalized therapy. Citation Format: Recep Bayraktar, Yitao Tang, Mihnea P. Dragomir, Linda Fabris, Giulio F. Draetta, Han Liang, George A. Calin. The mutational landscape of ultraconserved elements in human cancers. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 5750.

Glutamine Metabolism Underlies the Functional Similarity of T Cells between Nile Tilapia and Tetrapod.

As the lowest organisms possessing T cells, fish are instrumental for understanding T cell evolution and immune defense in early vertebrates. This study established in Nile tilapia models suggests that T cells play a critical role in resisting Edwardsiella piscicida infection via cytotoxicity and are essential for IgM+ B cell response. CD3 and CD28 monoclonal antibody crosslinking reveals that full activation of tilapia T cells requires the first and secondary signals, while Ca2+ -NFAT, MAPK/ERK, NF-κB, and mTORC1 pathways and IgM+ B cells collectively regulate T cell activation. Thus, despite the large evolutionary distance, tilapia and mammals such as mice and humans exhibit similar T cell functions. Furthermore, it is speculated that transcriptional networks and metabolic reprogramming, especially c-Myc-mediated glutamine metabolism triggered by mTORC1 and MAPK/ERK pathways, underlie the functional similarity of T cells between tilapia and mammals. Notably, tilapia, frogs, chickens, and mice utilize the same mechanisms to facilitate glutaminolysis-regulated T cell responses, and restoration of the glutaminolysis pathway using tilapia components rescues the immunodeficiency of human Jurkat T cells. Thus, this study provides a comprehensive picture of T cell immunity in tilapia, sheds novel perspectives for understanding T cell evolution, and offers potential avenues for intervening in human immunodeficiency.

Ancient origin and constrained evolution of the division and cell wall gene cluster in Bacteria.

The division and cell wall (dcw) gene cluster in Bacteria comprises 17 genes encoding key steps in peptidoglycan synthesis and cytokinesis. To understand the origin and evolution of this cluster, we analysed its presence in over 1,000 bacterial genomes. We show that the dcw gene cluster is strikingly conserved in both gene content and gene order across all Bacteria and has undergone only a few rearrangements in some phyla, potentially linked to cell envelope specificities, but not directly to cell shape. A large concatenation of the 12 most conserved dcw cluster genes produced a robust tree of Bacteria that is largely consistent with recent phylogenies based on frequently used markers. Moreover, evolutionary divergence analyses show that the dcw gene cluster offers advantages in defining high-rank taxonomic boundaries and indicate at least two main phyla in the Candidate Phyla Radiation (CPR) matching a sharp dichotomy in dcw gene cluster arrangement. Our results place the origin of the dcw gene cluster in the Last Bacterial Common Ancestor and show that it has evolved vertically for billions of years, similar to major cellular machineries such as the ribosome. The strong phylogenetic signal, combined with conserved genomic synteny at large evolutionary distances, makes the dcw gene cluster a robust alternative set of markers to resolve the ever-growing tree of Bacteria.

Multispecies transcriptomes reveal core fruit development genes.

During angiosperm evolution there have been repeated transitions from an ancestral dry fruit to a derived fleshy fruit, often with dramatic ecological and economic consequences. Following the transition to fleshy fruits, domestication may also dramatically alter the fruit phenotype via artificial selection. Although the morphologies of these fruits are well documented, relatively less is known about the molecular basis of these developmental and evolutionary shifts. We generated RNA-seq libraries from pericarp tissue of desert tobacco and both cultivated and wild tomato species at common developmental time points and combined this with corresponding, publicly available data from Arabidopsis and melon. With this broadly sampled dataset consisting of dry/fleshy fruits and wild/domesticated species, we applied novel bioinformatic methods to investigate conserved and divergent patterns of gene expression during fruit development and evolution. A small set of 121 orthologous "core" fruit development genes show a common pattern of expression across all five species. These include key players in developmental patterning such as orthologs of KNOLLE, PERIANTHIA, and ARGONAUTE7. GO term enrichment suggests that these genes function in basic cell division processes, cell wall biosynthesis, and developmental patterning. We furthermore uncovered a number of "accessory" genes with conserved expression patterns within but not among fruit types, and whose functional enrichment highlights the conspicuous differences between these phenotypic classes. We observe striking conservation of gene expression patterns despite large evolutionary distances, and dramatic phenotypic shifts, suggesting a conserved function for a small subset of core fruit development genes.

Structure of a mitochondrial ribosome with fragmented rRNA in complex with membrane-targeting elements

Mitoribosomes of green algae display a great structural divergence from their tracheophyte relatives, with fragmentation of both rRNA and proteins as a defining feature. Here, we report a 2.9 Å resolution structure of the mitoribosome from the alga Polytomella magna harbouring a reduced rRNA split into 13 fragments. We found that the rRNA contains a non-canonical reduced form of the 5S, as well as a permutation of the LSU domain I. The mt-5S rRNA is stabilised by mL40 that is also found in mitoribosomes lacking the 5S, which suggests an evolutionary pathway. Through comparison to other ribosomes with fragmented rRNAs, we observe that the pattern is shared across large evolutionary distances, and between cellular compartments, indicating an evolutionary convergence and supporting the concept of a primordial fragmented ribosome. On the protein level, eleven peripherally associated HEAT-repeat proteins are involved in the binding of 3′ rRNA termini, and the structure features a prominent pseudo-trimer of one of them (mL116). Finally, in the exit tunnel, mL128 constricts the tunnel width of the vestibular area, and mL105, a homolog of a membrane targeting component mediates contacts with an inner membrane bound insertase. Together, the structural analysis provides insight into the evolution of the ribosomal machinery in mitochondria.

Inferring mammalian tissue-specific regulatory conservation by predicting tissue-specific differences in open chromatin

BackgroundEvolutionary conservation is an invaluable tool for inferring functional significance in the genome, including regions that are crucial across many species and those that have undergone convergent evolution. Computational methods to test for sequence conservation are dominated by algorithms that examine the ability of one or more nucleotides to align across large evolutionary distances. While these nucleotide alignment-based approaches have proven powerful for protein-coding genes and some non-coding elements, they fail to capture conservation of many enhancers, distal regulatory elements that control spatial and temporal patterns of gene expression. The function of enhancers is governed by a complex, often tissue- and cell type-specific code that links combinations of transcription factor binding sites and other regulation-related sequence patterns to regulatory activity. Thus, function of orthologous enhancer regions can be conserved across large evolutionary distances, even when nucleotide turnover is high.ResultsWe present a new machine learning-based approach for evaluating enhancer conservation that leverages the combinatorial sequence code of enhancer activity rather than relying on the alignment of individual nucleotides. We first train a convolutional neural network model that can predict tissue-specific open chromatin, a proxy for enhancer activity, across mammals. Next, we apply that model to distinguish instances where the genome sequence would predict conserved function versus a loss of regulatory activity in that tissue. We present criteria for systematically evaluating model performance for this task and use them to demonstrate that our models accurately predict tissue-specific conservation and divergence in open chromatin between primate and rodent species, vastly out-performing leading nucleotide alignment-based approaches. We then apply our models to predict open chromatin at orthologs of brain and liver open chromatin regions across hundreds of mammals and find that brain enhancers associated with neuron activity have a stronger tendency than the general population to have predicted lineage-specific open chromatin.ConclusionThe framework presented here provides a mechanism to annotate tissue-specific regulatory function across hundreds of genomes and to study enhancer evolution using predicted regulatory differences rather than nucleotide-level conservation measurements.

Biogenesis and metabolic homeostasis of trypanosomatid glycosomes: New insights and new questions

Kinetoplastea and Diplonemea possess peroxisome-related organelles that, uniquely, contain most of the enzymes of the glycolytic pathway and are hence called glycosomes. Enzymes of several other core metabolic pathways have also been located in glycosomes, in addition to some characteristic peroxisomal systems such as pathways of lipid metabolism. A considerable amount of research has been performed on glycosomes of trypanosomes since their discovery four decades ago. Not only the role of the glycosomal enzyme systems in the overall cell metabolism appeared to be unique, but also the organelles display remarkable features regarding their biogenesis and structural properties. These features are similar to those of the well-studied peroxisomes of mammalian and plant cells and yeasts yet exhibit also differences reflecting the large evolutionary distance between these protists and the representatives of other major eukaryotic lineages. Despite all research performed, many questions remain about various properties and the biological roles of glycosomes and peroxisomes. Here, we review the current knowledge about glycosomes, often comparing it with information about peroxisomes. Furthermore, we highlight particularly many questions that remain about the biogenesis, and the heterogeneity in structure and content of these enigmatic organelles, and the properties of their boundary membrane.

Phylogenetic Tests of Models of Viral Transmission

The hunt for the immediate non-human host of SARS-CoV-2 has centered on bats of the genus Rhinolophus. We explored the phylogenetic predictions of two models of viral transmission, the Spillover Model and the Circulation Model and suggest that the Spillover Model can be eliminated. The Circulation Model suggests that viral transmission occurs among susceptible hosts irrespective of their phylogenetic relationships. Susceptibility could be mediated by the ACE2 gene (important for viral docking) and we constructed a phylogeny of this gene for 159 mammal species, finding a phylogenetic pattern consistent with established mammalian relationships. The tree indicates that viral transfer occurs over large evolutionary distances. Although lacking consensus, some studies identify a virus from a particular R. affinis individual (RaTG13) as being most closely related phylogenetically to human SARS-CoV-2. However, other R. affinis harbor viruses that are relatively unrelated to human viruses, and viruses found in this species exhibit sequence differences of up to 20%, suggesting multiple transfers over time. There is little correspondence between viral and host (bat) species limits or phylogenetic relationships. An ACE2 phylogeny for Rhinolophus followed species limits, unlike the pattern in the viral phylogeny indicating that phylogenetic similarity of ACE2 is not a predictor of viral transmission at the bat species level. The Circulation Model could be modified to apply to any individual of any species of Rhinolophus; more individuals and species must be examined.

Complete mitochondrial genome of Rhizosolenia setigera (Coscinodiscophyceae, Bacillariophyta)

Rhizosolenia is a species-rich genus with 144 described species, many of which are harmful algal species (HABs) with significant negative ecological impact. Despite their significance in primary production and their potential to induce HABs, genome data of these species remain extremely limited. In this study, the complete mitochondrial genome (mtDNA) of Rhizosolenia setigera Brightwell 1858 was determined for the first time, which also represented the first mtDNA of the order Rhizosoleniales. The circular mtDNA was 34,792 bp in length with GC content of 23.28%. It encoded 63 genes including 35 protein-coding genes (PCGs), 24 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes, and 2 conserved open reading frames (orfs). Phylogenetic analysis using concatenated PCGs revealed that R. setigera and Melosira undulate, which also belongs to the class Coscinodiscophyceae, clustered together as expected. However, comparison of these two mtDNAs revealed extensive genome rearrangement events, suggesting large evolutionary distance. The complete mtDNA of R. setigera will facilitate research on the phylogenetic relationship among Rhizosolenia species, which will in turn facilitate exploration of the evolutionary relationships in the class of Coscinodiscophyceae.

The cytochrome P450 (CYP) superfamily in cnidarians

The cytochrome P450 (CYP) superfamily is a diverse and important enzyme family, playing a central role in chemical defense and in synthesis and metabolism of major biological signaling molecules. The CYPomes of four cnidarian genomes (Hydra vulgaris, Acropora digitifera, Aurelia aurita, Nematostella vectensis) were annotated; phylogenetic analyses determined the evolutionary relationships amongst the sequences and with existing metazoan CYPs. 155 functional CYPs were identified and 90 fragments. Genes were from 24 new CYP families and several new subfamilies; genes were in 9 of the 12 established metazoan CYP clans. All species had large expansions of clan 2 diversity, with H. vulgaris having reduced diversity for both clan 3 and mitochondrial clan. We identified potential candidates for xenobiotic metabolism and steroidogenesis. That each genome contained multiple, novel CYP families may reflect the large evolutionary distance within the cnidarians, unique physiology in the cnidarian classes, and/or different ecology of the individual species.

Classification of Precursor MicroRNAs from Different Species Based on K-mer Distance Features

MicroRNAs (miRNAs) are short RNA sequences that are actively involved in gene regulation. These regulators on the post-transcriptional level have been discovered in virtually all eukaryotic organisms. Additionally, miRNAs seem to exist in viruses and might also be produced in microbial pathogens. Initially, transcribed RNA is cleaved by Drosha, producing precursor miRNAs. We have previously shown that it is possible to distinguish between microRNA precursors of different clades by representing the sequences in a k-mer feature space. The k-mer representation considers the frequency of a k-mer in the given sequence. We further hypothesized that the relationship between k-mers (e.g., distance between k-mers) could be useful for classification. Three different distance-based features were created, tested, and compared. The three feature sets were entitled inter k-mer distance, k-mer location distance, and k-mer first–last distance. Here, we show that classification performance above 80% (depending on the evolutionary distance) is possible with a combination of distance-based and regular k-mer features. With these novel features, classification at closer evolutionary distances is better than using k-mers alone. Combining the features leads to accurate classification for larger evolutionary distances. For example, categorizing Homo sapiens versus Brassicaceae leads to an accuracy of 93%. When considering average accuracy, the novel distance-based features lead to an overall increase in effectiveness. On the contrary, secondary-structure-based features did not lead to any effective separation among clades in this study. With this line of research, we support the differentiation between true and false miRNAs detected from next-generation sequencing data, provide an additional viewpoint for confirming miRNAs when the species of origin is known, and open up a new strategy for analyzing miRNA evolution.

One stem cell program to rule them all?

Many species of animals have stem cells that they maintain throughout their lives, which suggests that stem cells are an ancestral feature of all animals. From this, we take the viewpoint that cells with the biological properties of 'stemness'-self-renewal and multipotency-may share ancestral genetic circuitry. However, in practice is it very difficult to identify and compare stemness gene signatures across diverse animals and large evolutionary distances? First, it is critical to experimentally demonstrate self-renewal and potency. Second, genomic methods must be used to determine specific gene expression in stem cell types compared with non-stem cell types to determine stem cell gene enrichment. Third, gene homology must be mapped between diverse animals across large evolutionary distances. Finally, conserved genes that fulfill these criteria must be tested for role in stem cell function. It is our viewpoint that by comparing stem cell-specific gene signatures across evolution, ancestral programs of stemness can be uncovered, and ultimately, the dysregulation of stemness programs drives the state of cancer stem cells.

The ranging of amino acids substitution matrices of various types in accordance with the alignment accuracy criterion

BackgroundThe alignment of character sequences is important in bioinformatics. The quality of this procedure is determined by the substitution matrix and parameters of the insertion-deletion penalty function. These matrices are derived from sequence alignment and thus reflect the evolutionary process. Currently, in addition to evolutionary matrices, a large number of different background matrices have been obtained. To make an optimal choice of the substitution matrix and the penalty parameters, we conducted a numerical experiment using a representative sample of existing matrices of various types and origins.ResultsWe tested both the classical evolutionary matrix series (PAM, Blosum, VTML, Pfasum); structural alignment based matrices, contact energy matrix, and matrix based on the properties of the genetic code. This study presents results for two test set types: first, we simulated sequences that reflect the divergent evolution; second, we performed tests on Balibase sequences. In both cases, we obtained the dependences of the alignment quality (Accuracy, Confidence) on the evolutionary distance between sequences and the evolutionary distance to which the substitution matrices correspond. Optimization of a combination of matrices and the penalty parameters was carried out for local and global alignment on the values of penalty function parameters.Consequently, we found that the best alignment quality is achieved with matrices corresponding to the largest evolutionary distance. These matrices prove to be universal, i.e. suitable for aligning sequences separated by both large and small evolutionary distances. We analysed the correspondence of the correlation coefficients of matrices to the alignment quality. It was found that matrices showing high quality alignment have an above average correlation value, but the converse is not true.ConclusionsThis study showed that the best alignment quality is achieved with evolutionary matrices designed for long distances: Gonnet, VTML250, PAM250, MIQS, and Pfasum050. The same property is inherent in matrices not only of evolutionary origin, but also of another background corresponding to a large evolutionary distance. Therefore, matrices based on structural data show alignment quality close enough to its value for evolutionary matrices. This agrees with the idea that the spatial structure is more conservative than the protein sequence.

Computational analysis identifies druggable mutations in human rBAT mediated Cystinuria

Culex quinquefasciatus Cqm1 protein acts as the receptor for Lysinibacillus sphaericus mosquito-larvicidal binary (BinAB) toxin that is used worldwide for mosquito control. We found amino acid transporter protein, rBAT, as phylogenetically closest Cqm1 homolog in humans. The present study reveals large evolutionary distance between Cqm1 and rBAT, and rBAT ectodomain lacks the sequence motif which serves as binding-site for the BinAB toxin. Thus, BinAB toxin can be expected to remain safe for humans. rBAT (heavy subunit; SLC3A1) and catalytic b0,+AT (light subunit; SLC7A9), linked by single disulfide bond, mediate renal reabsorption of cystine and dibasic amino acids in Na+ independent manner. Mutations in rBAT cause type I Cystinuria disease which shows global prevalence, and rBAT can be thought as an important pharmacological target. However, 3D structures of rBAT and b0,+AT, the two components of b0,+ heteromeric amino acid transporter systems, are not available. We constructed a reliable homology model of rBAT using Cqm1 coordinates and that of transmembrane b0,+AT subunit using LAT1 coordinates. Mapping of pathogenic mutations onto rBAT ectodomain revealed their scattered distribution throughout the rBAT protein. Further, our computational simulations-based scoring of several known deleterious mutations of rBAT revealed that mutations those do not compromise the protein fold and stability, are localized on the same face of the molecule. These residues are expected to interact with the b0,+AT transporter. The present study thus identifies druggable sites on rBAT that could be targeted for the treatment of type I Cystinuria. Communicated by Ramaswamy H. Sarma

Functional Conservation of Divergent p63-Bound cis-Regulatory Elements.

The transcription factor p63 is an essential regulator of vertebrate ectoderm development, including epidermis, limbs, and craniofacial tissues. Here, we have investigated the evolutionary conservation of p63 binding sites (BSs) between zebrafish and human. First, we have analyzed sequence conservation of p63 BSs by comparing ChIP-seq data from human keratinocytes and zebrafish embryos, observing a very poor conservation. Next, we compared the gene regulatory network orchestrated by p63 in both species and found a high overlap between them, suggesting a high degree of functional conservation during evolution despite sequence divergence and the large evolutionary distance. Finally, we used transgenic reporter assays in zebrafish embryos to functionally validate a set of equivalent p63 BSs from zebrafish and human located close to genes involved in epidermal development. Reporter expression was driven by human and zebrafish BSs to many common tissues related to p63 expression domains. Therefore, we conclude that the gene regulatory network controlled by p63 is highly conserved across vertebrates despite the fact that p63-bound regulatory elements show high divergence.

Read-SpaM: assembly-free and alignment-free comparison of bacterial genomes with low sequencing coverage

BackgroundIn many fields of biomedical research, it is important to estimate phylogenetic distances between taxa based on low-coverage sequencing reads. Major applications are, for example, phylogeny reconstruction, species identification from small sequencing samples, or bacterial strain typing in medical diagnostics.ResultsWe adapted our previously developed software program Filtered Spaced-Word Matches (FSWM) for alignment-free phylogeny reconstruction to take unassembled reads as input; we call this implementation Read-SpaM.ConclusionsTest runs on simulated reads from semi-artificial and real-world bacterial genomes show that our approach can estimate phylogenetic distances with high accuracy, even for large evolutionary distances and for very low sequencing coverage.

A Non-amyloid Prion Particle that Activates a Heritable Gene Expression Program.

Spatiotemporal gene regulation is often driven by RNA-binding proteins that harbor long intrinsically disordered regions in addition to folded RNA-binding domains. We report that the disordered region of theevolutionarily ancient developmental regulator Vts1/Smaug drives self-assembly into gel-like condensates. These proteinaceous particles are not composed of amyloid, yet they are infectious, allowing them to act as a protein-based epigenetic element: a prion [SMAUG+]. In contrast to many amyloid prions, condensation of Vts1 enhances its function in mRNA decay, and its self-assembly properties are conserved over large evolutionary distances. Yeast cells harboring [SMAUG+] downregulate a coherent network of mRNAs and exhibit improved growth under nutrient limitation. Vts1 condensates formed from purified protein can transform naive cells to acquire [SMAUG+]. Our data establish that non-amyloid self-assembly of RNA-binding proteins can drive a form of epigenetics beyond the chromosome, instilling adaptive gene expression programs that are heritable over long biological timescales.

Developmental transcriptomes of the sea star, Patiria miniata, illuminate how gene expression changes with evolutionary distance

Understanding how changes in developmental gene expression alter morphogenesis is a fundamental problem in development and evolution. A promising approach to address this problem is to compare the developmental transcriptomes between related species. The echinoderm phylum consists of several model species that have significantly contributed to the understanding of gene regulation and evolution. Particularly, the regulatory networks of the sea star, Patiria miniata (P. miniata), have been extensively studied, however developmental transcriptomes for this species were lacking. Here we generated developmental transcriptomes of P. miniata and compared these with those of two sea urchins species. We demonstrate that the conservation of gene expression depends on gene function, cell type and evolutionary distance. With increasing evolutionary distance the interspecies correlations in gene expression decreases. The reduction is more severe in the correlations between morphologically equivalent stages (diagonal elements) than in the correlation between morphologically distinct stages (off-diagonal elements). This could reflect a decrease in the morphological constraints compared to other constraints that shape gene expression at large evolutionary divergence. Within this trend, the interspecies correlations of developmental control genes maintain their diagonality at large evolutionary distance, and peak at the onset of gastrulation, supporting the hourglass model of phylotypic stage conservation.