Codetta: predicting the genetic code from nucleotide sequence

Single-cell transcriptomic analysis reveals genome evolution in predatory litostomatean ciliates

Many ciliated protists prey on other large microbial organisms, including other protists and microscopic metazoans. The ciliate class Litostomatea unites both predatory and endosymbiotic species. The evolution of predation ability in ciliates remains poorly understood, in part, due to a lack of genomic data. To fill this gap, we acquired the transcriptome profiles of six predatory litostomateans using single-cell sequencing technology and investigated their transcriptomic features. Our results show that: (1) in contrast to non-predatory ciliates, the predatory litostomateans have expanded gene families associated with transmembrane activity and reactive oxidative stress response pathways, potentially as a result of cellular behaviors such as fast contraction and extension; (2) the expansion of the calcium-activated BK potassium channel gene family, which hypothetically regulates cell contractility, is an ancient evolutionary event for the class Litostomatea, suggesting a rewired metabolism associated with the hunting behavior of predatory ciliates; and (3) three whole genome duplication (WGD) events have been detected in litostomateans, with genes associated with biosynthetic processes, transmembrane activity, and calcium-activated potassium channel activity being retained during the WGD events. In addition, we explored the evolutionary relationships among 17 ciliate species, including eight litostomateans, and provided a rich foundational dataset for future in-depth phylogenomic studies of Litostomatea. Our comprehensive analyses suggest that the rewired cellular metabolism via expanded gene families and WGD events might be the potential genetic basis for the predation ability of raptorial ciliates.

Single cell transcriptomics reveals UAR codon reassignment in Palmarella salina (Metopida, Armophorea) and confirms Armophorida belongs to APM clade

Anaerobes have emerged in several major lineages of ciliates, but the number of independent transitions to anaerobiosis among ciliates is unknown. The APM clade (Armophorea, Muranotrichea, Parablepharismea) represents the largest clade of obligate anaerobes among ciliates and contains free-living marine and freshwater representatives as well as gut endobionts of animals. The evolution of APM group has only recently started getting attention, and our knowledge on its phylogeny and genetics is still limited to a fraction of taxa. While ciliates portray a wide array of alternatives to the standard genetic code across numerous classes, the APM ciliates were considered to be the largest group using exclusively standard nuclear genetic code. In this study, we present a pan-ciliate phylogenomic analysis with emphasis on the APM clade, bringing the first phylogenomic analysis of the family Tropidoatractidae (Armophorea) and confirming the position of Armophorida within Armophorea. We include five newly sequenced single cell transcriptomes from marine, freshwater, and endobiotic APM ciliates - Palmarella salina, Anteclevelandella constricta, Nyctotherus sp., Caenomorpha medusula, and Thigmothrix strigosa. We report the first discovery of an alternative nuclear genetic code among APM ciliates, used by Palmarella salina (Tropidoatractidae, Armophorea), but not by its close relative, Tropidoatractus sp., and provide a comparative analysis of stop codon identity and frequency indicating the precedency to the UAG codon loss/reassignment over the UAA codon reassignment in the specific ancestor of Palmarella. Comparative genomic and proteomic studies of this group may help explain the constraints that underlie UAR stop-to-sense reassignment, the most frequent type of alternative nuclear genetic code, not only in ciliates, but eukaryotes in general.

The P10K database: a data portal for the protist 10 000 genomes project

Protists, a highly diverse group of microscopic eukaryotic organisms distinct from fungi, animals and plants, exert crucial roles within the earth's biosphere. However, the genomes of only a small fraction of known protist species have been published and made publicly accessible. To address this constraint, the Protist 10 000 Genomes Project (P10K) was initiated, implementing a specialized pipeline for single-cell genome/transcriptome assembly, decontamination and annotation of protists. The resultant P10K database (https://ngdc.cncb.ac.cn/p10k/) serves as a comprehensive platform, collating and disseminating genome sequences and annotations from diverse protist groups. Currently, the P10K database has incorporated 2959 genomes and transcriptomes, including 1101 newly sequenced datasets by P10K and 1858 publicly available datasets. Notably, it covers 45% of the protist orders, with a significant representation (53% coverage) of ciliates, featuring nearly a thousand genomes/transcriptomes. Intriguingly, analysis of the unique codon table usage among ciliates has revealed differences compared to the NCBI taxonomy system, suggesting a need to revise the codon tables used for these species. Collectively, the P10K database serves as a valuable repository of genetic resources for protist research and aims to expand its collection by incorporating more sequenced data and advanced analysis tools to benefit protist studies worldwide.

Multiple independent losses of sporulation and peptidoglycan in the Mycoplasmatales and related orders of the class Bacilli

Many peptidoglycan-deficient bacteria such as the Mycoplasmatales are known host-associated lineages, lacking the environmental resistance mechanisms and metabolic capabilities necessary for a free-living lifestyle. Several peptidoglycan-deficient and non-sporulating orders of interest are thought to be descended from Gram-positive sporulating Bacilli through reductive evolution. Here we annotate 2650 genomes belonging to the class Bacilli, according to the Genome Taxonomy Database, to predict the peptidoglycan and sporulation phenotypes of three novel orders, RFN20, RF39 and ML615J-28, known only through environmental sequence surveys. These lineages are interspersed between peptidoglycan-deficient non-sporulating orders including the Mycoplasmatales and Acholeplasmatales, and more typical Gram-positive orders such as the Erysipelotrichales and Staphylococcales. We use the extant genotypes to perform ancestral state reconstructions. The novel orders are predicted to have small genomes with minimal metabolic capabilities and to comprise a mix of peptidoglycan-deficient and/or non-sporulating species. In contrast to expectations based on cultured representatives, the order Erysipelotrichales lacks many of the genes involved in peptidoglycan and endospore formation. The reconstructed evolutionary history of these traits suggests multiple independent whole-genome reductions and loss of phenotype via intermediate transition states that continue into the present. We suggest that the evolutionary history of the reduced-genome lineages within the class Bacilli is one driven by multiple independent transitions to host-associated lifestyles, with the degree of reduction in environmental resistance and metabolic capabilities correlated with degree of host association.

The Genetic Code Assembles via Division and Fusion, Basic Cellular Events

Standard Genetic Code (SGC) evolution is quantitatively modeled in up to 2000 independent coding 'environments'. Environments host multiple codes that may fuse or divide, with division yielding identical descendants. Code division may be selected-sophisticated gene products could be required for an orderly separation that preserves the coding. Several unforeseen results emerge: more rapid evolution requires unselective code division rather than its selective form. Combining selective and unselective code division, with/without code fusion, with/without independent environmental coding tables, and with/without wobble defines 25 = 32 possible pathways for SGC evolution. These 32 possible histories are compared, specifically, for evolutionary speed and code accuracy. Pathways differ greatly, for example, by ≈300-fold in time to evolve SGC-like codes. Eight of thirty-two pathways employing code division evolve quickly. Four of these eight that combine fusion and division also unite speed and accuracy. The two most precise, swiftest paths; thus the most likely routes to the SGC are similar, differing only in fusion with independent environmental codes. Code division instead of fusion with unrelated codes implies that exterior codes can be dispensable. Instead, a single ancestral code that divides and fuses can initiate fully encoded peptide biosynthesis. Division and fusion create a 'crescendo of competent coding', facilitating the search for the SGC and also assisting the advent of otherwise uniformly disfavored wobble coding. Code fusion can unite multiple codon assignment mechanisms. However, via code division and fusion, an SGC can emerge from a single primary origin via familiar cellular events.

Identification of a non-canonical ciliate nuclear genetic code where UAA and UAG code for different amino acids

The genetic code is one of the most highly conserved features across life. Only a few lineages have deviated from the "universal" genetic code. Amongst the few variants of the genetic code reported to date, the codons UAA and UAG virtually always have the same translation, suggesting that their evolution is coupled. Here, we report the genome and transcriptome sequencing of a novel uncultured ciliate, belonging to the Oligohymenophorea class, where the translation of the UAA and UAG stop codons have changed to specify different amino acids. Genomic and transcriptomic analyses revealed that UAA has been reassigned to encode lysine, while UAG has been reassigned to encode glutamic acid. We identified multiple suppressor tRNA genes with anticodons complementary to the reassigned codons. We show that the retained UGA stop codon is enriched in the 3'UTR immediately downstream of the coding region of genes, suggesting that there is functional drive to maintain tandem stop codons. Using a phylogenomics approach, we reconstructed the ciliate phylogeny and mapped genetic code changes, highlighting the remarkable number of independent genetic code changes within the Ciliophora group of protists. According to our knowledge, this is the first report of a genetic code variant where UAA and UAG encode different amino acids.