FACIL: Fast and Accurate Genetic Code Inference and Logo

A cryptic plastid and a novel mitochondrial plasmid in Leucomyxa plasmidifera gen. and sp. nov. (Ochrophyta) push the frontiers of organellar biology

Complete plastid loss seems to be very rare among secondarily non-photosynthetic eukaryotes. Leukarachnion sp. PRA-24, an amoeboid colourless protist related to the photosynthetic algal class Synchromophyceae (Ochrophyta), is a candidate for such a case based on a previous investigation by transmission electron microscopy. Here, we characterize this organism in further detail and describe it as Leucomyxa plasmidifera gen. et sp. nov., additionally demonstrating it is the first known representative of a broader clade of non-photosynthetic ochrophytes. We recovered its complete plastid genome, exhibiting a reduced gene set similar to plastomes of other non-photosynthetic ochrophytes, yet being even more extreme in sequence divergence. Identification of components of the plastid protein import machinery in the L. plasmidifera transcriptome assembly corroborated that the organism possesses a cryptic plastid organelle. According to our bioinformatic reconstruction, the plastid contains a unique combination of biosynthetic pathways producing haem, a folate precursor and tocotrienols. As another twist to its organellar biology, L. plasmidifera turned out to contain an unusual long insertion in its mitogenome related to a newly discovered mitochondrial plasmid exhibiting unprecedented features in terms of its size and coding capacity. Combined, our work uncovered further striking outcomes of the evolutionary course of semiautonomous organelles in protists.

Dynamics and timing of diversification events of ciliated eukaryotes from a large phylogenomic perspective

Ciliophora, an exceptionally diverse lineage of unicellular eukaryotes, exhibits a remarkable range of species richness across classes in the ciliate Tree of Life. In this study, we have acquired transcriptome and genome data from 40 representative species in seven ciliate classes. Utilizing 247 genes and 105 taxa, we devised a comprehensive phylogenomic tree for Ciliophora, encompassing over 60 % of orders and constituting the most extensive dataset of ciliate species to date. We established a robust phylogenetic framework that encompasses ambiguous taxa and the major classes within the phylum. Our findings support the monophyly of each of two subphyla (Postciliodesmatophora and Intramacronucleata), along with three subclades (Protocruzia, CONTHREEP, and SAPML) nested within Intramacronucleata, and elucidate evolutionary positions among the major classes within the phylum. Drawing on the robust ciliate Tree of Life and three constraints, we estimated the radiation of Ciliophora around 1175 Ma during the middle of the Proterozoic Eon, and most of the ciliate classes diverged from their sister lineage during the latter half of this period. Additionally, based on the time-calibrated tree and species richness pattern, we investigated net diversification rates of Ciliophora and its classes. The global net diversification rate for Ciliophora was estimated at 0.004979 species/Ma. Heterogeneity in net diversification rates was evident at the class level, with faster rates observed in Oligohymenophorea and Spirotrichea than other classes within the subclades CONTHREEP and SAPML, respectively. Notably, our analysis suggests that variations in net diversification rates, rather than clade ages, appear to contribute to the differences in species richness in Ciliophora at the class level.

Giants' cooperation: a draft genome of the giant ciliate Muniziella cunhai suggests its ecological role in the capybara's digestive metabolism

Several hundred ciliate species live in animals' guts as a part of their microbiome. Among them, Muniziella cunhai (Trichostomatia, Pycnotrichidae), the largest described ciliate, is found exclusively associated with Hydrochoerus hydrochaeris (capybara), the largest known rodent reaching up to 90 kg. Here, we present the sequence, structural and functional annotation of this giant microeukaryote macronuclear genome and discuss its phylogenetic placement. The 85 Mb genome is highly AT rich (GC content 25.71 %) and encodes a total of 11 397 protein-coding genes, of which 2793 could have their functions predicted with automated functional assignments. Functional annotation showed that M. cunhai can digest recalcitrant structural carbohydrates, non-structural carbohydrates, and microbial cell walls, suggesting a role in diet metabolization and in microbial population control in the capybara's intestine. Moreover, the phylogenetic placement of M. cunhai provides insights on the origins of gigantism in the subclass Trichostomatia.

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.

MgCod: Gene Prediction in Phage Genomes with Multiple Genetic Codes

Massive sequencing of microbiomes has led to the discovery of a large number of phage genomes with intermittent stop codon recoding. We have developed a computational tool, MgCod, that identifies genomic regions (blocks) with distinct stop codon recoding simultaneously with the prediction of protein-coding regions. When MgCod was used to scan a large volume of human metagenomic contigs hundreds of viral contigs with intermittent stop codon recoding were revealed. Many of these contigs originated from genomes of known crAssphages. Further analyses had shown that intermittent recoding was associated with subtle patterns in the organization of protein-coding genes, such as 'single-coding' and 'dual-coding'. The dual-coding genes, clustered into blocks, could be translated by two alternative codes producing nearly identical proteins. It was observed that the dual-coded blocks were enriched with the early-stage phage genes, while the late-stage genes were residing in the single-coded blocks. MgCod can identify types of stop codon recoding in novel genomic sequences in parallel with gene prediction. It is available for download from https://github.com/gatech-genemark/MgCod.

Expression of Hairpin-Enriched Mitochondrial DNA in Two Hairworm Species (Nematomorpha)

Nematomorpha (hairworms) is a phylum of parasitic ecdysozoans, best known for infecting arthropods and guiding their hosts toward water, where the parasite can complete its life cycle. Over 350 species of nematomorphs have been described, yet molecular data for the group remain scarce. The few available mitochondrial genomes of nematomorphs are enriched with long inverted repeats, which are embedded in the coding sequences of their genes-a remarkably unusual feature exclusive to this phylum. Here, we obtain and annotate the repeats in the mitochondrial genome of another nematomorph species-Parachordodes pustulosus. Using genomic and transcriptomic libraries, we investigate the impact of inverted repeats on the read coverage of the mitochondrial genome. Pronounced drops in the read coverage coincide with regions containing long inverted repeats, denoting the 'blind spots' of short-fragment sequencing libraries. Phylogenetic inference with the novel data reveals multiple disagreements between the traditional system of Nematomorpha and molecular data, rendering several genera paraphyletic, including Parachordodes.

In silico characterization of cysteine-stabilized αβ defensins from neglected unicellular microeukaryotes

BACKGROUND

The emergence of multi-resistant pathogens have increased dramatically in recent years, becoming a major public-health concern. Among other promising antimicrobial molecules with potential to assist in this worldwide struggle, cysteine-stabilized αβ (CS-αβ) defensins are attracting attention due their efficacy, stability, and broad spectrum against viruses, bacteria, fungi, and protists, including many known human pathogens.

RESULTS

Here, 23 genomes of ciliated protists were screened and two CS-αβ defensins with a likely antifungal activity were identified and characterized, using bioinformatics, from a culturable freshwater species, Laurentiella sp. (LsAMP-1 and LsAMP-2). Although any potential cellular ligand could be predicted for LsAMP-2; evidences from structural, molecular dynamics, and docking analyses suggest that LsAMP-1 may form stably associations with phosphatidylinositol 4,5-bisphosphates (PIP2), a phospholipid found on many eukaryotic cells, which could, in turn, represent an anchorage mechanism within plasma membrane of targeted cells.

CONCLUSION

These data stress that more biotechnology-oriented studies should be conducted on neglected protists, such ciliates, which could become valuable sources of novel bioactive molecules for therapeutic uses.

Codetta: predicting the genetic code from nucleotide sequence

SUMMARY

Codetta is a Python program for predicting the genetic code table of an organism from nucleotide sequences. Codetta can analyze an arbitrary nucleotide sequence and needs no sequence annotation or taxonomic placement. The most likely amino acid decoding for each of the 64 codons is inferred from alignments of profile hidden Markov models of conserved proteins to the input sequence.

AVAILABILITY AND IMPLEMENTATION

Codetta 2.0 is implemented as a Python 3 program for MacOS and Linux and is available from http://eddylab.org/software/codetta/codetta2.tar.gz and at http://github.com/kshulgina/codetta.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genomic analysis reveals cryptic diversity in aphelids and sheds light on the emergence of Fungi

Over the past decade, molecular phylogenetics has reshaped our understanding of the fungal tree of life by unraveling a hitherto elusive diversity of the protistan relatives of Fungi. Aphelida constitutes one of these novel deep branches that precede the emergence of osmotrophic fungal lifestyle and hold particular significance as the pathogens of algae. Here, we obtain and analyze the genomes of aphelid species Amoeboaphelidium protococcarum and Amoeboaphelidium occidentale. Genomic data unmask the vast divergence between these species, hidden behind their morphological similarity, and reveal hybrid genomes with a complex evolutionary history in two strains of A. protococcarum. We confirm the proposed sister relationship between Aphelida and Fungi using phylogenomic analysis and chart the reduction of characteristic proteins involved in phagocytic activity in the evolution of Holomycota. Annotation of aphelid genomes demonstrates the retention of actin nucleation-promoting complexes associated with phagocytosis and amoeboid motility and also reveals a conspicuous expansion of receptor-like protein kinases, uncharacteristic of fungal lineages. We find that aphelids possess multiple carbohydrate-processing enzymes that are involved in fungal cell wall synthesis but do not display rich complements of algal cell-wall-processing enzymes, suggesting an independent origin of fungal plant-degrading capabilities. Aphelid genomes show that the emergence of Fungi from phagotrophic ancestors relied on a common cell wall synthetic machinery but required a different set of proteins for digestion and interaction with the environment.

Computational Tools for the Analysis of Uncultivated Phage Genomes

Over a century of bacteriophage research has uncovered a plethora of fundamental aspects of their biology, ecology, and evolution. Furthermore, the introduction of community-level studies through metagenomics has revealed unprecedented insights on the impact that phages have on a range of ecological and physiological processes. It was not until the introduction of viral metagenomics that we began to grasp the astonishing breadth of genetic diversity encompassed by phage genomes. Novel phage genomes have been reported from a diverse range of biomes at an increasing rate, which has prompted the development of computational tools that support the multilevel characterization of these novel phages based solely on their genome sequences. The impact of these technologies has been so large that, together with MAGs (Metagenomic Assembled Genomes), we now have UViGs (Uncultivated Viral Genomes), which are now officially recognized by the International Committee for the Taxonomy of Viruses (ICTV), and new taxonomic groups can now be created based exclusively on genomic sequence information. Even though the available tools have immensely contributed to our knowledge of phage diversity and ecology, the ongoing surge in software programs makes it challenging to keep up with them and the purpose each one is designed for. Therefore, in this review, we describe a comprehensive set of currently available computational tools designed for the characterization of phage genome sequences, focusing on five specific analyses: (i) assembly and identification of phage and prophage sequences, (ii) phage genome annotation, (iii) phage taxonomic classification, (iv) phage-host interaction analysis, and (v) phage microdiversity.

A computational screen for alternative genetic codes in over 250,000 genomes

The genetic code has been proposed to be a ‘frozen accident,’ but the discovery of alternative genetic codes over the past four decades has shown that it can evolve to some degree. Since most examples were found anecdotally, it is difficult to draw general conclusions about the evolutionary trajectories of codon reassignment and why some codons are affected more frequently. To fill in the diversity of genetic codes, we developed Codetta, a computational method to predict the amino acid decoding of each codon from nucleotide sequence data. We surveyed the genetic code usage of over 250,000 bacterial and archaeal genome sequences in GenBank and discovered five new reassignments of arginine codons (AGG, CGA, and CGG), representing the first sense codon changes in bacteria. In a clade of uncultivated Bacilli, the reassignment of AGG to become the dominant methionine codon likely evolved by a change in the amino acid charging of an arginine tRNA. The reassignments of CGA and/or CGG were found in genomes with low GC content, an evolutionary force that likely helped drive these codons to low frequency and enable their reassignment.

All life forms rely on a ‘code’ to translate their genetic information into proteins. This code relies on limited permutations of three nucleotides – the building blocks that form DNA and other types of genetic information. Each ‘triplet’ of nucleotides – or codon – encodes a specific amino acid, the basic component of proteins. Reading the sequence of codons in the right order will let the cell know which amino acid to assemble next on a growing protein. For instance, the codon CGG – formed of the nucleotides guanine (G) and cytosine (C) – codes for the amino acid arginine. From bacteria to humans, most life forms rely on the same genetic code. Yet certain organisms have evolved to use slightly different codes, where one or several codons have an altered meaning.

To better understand how alternative genetic codes have evolved, Shulgina and Eddy set out to find more organisms featuring these altered codons, creating a new software called Codetta that can analyze the genome of a microorganism and predict the genetic code it uses. Codetta was then used to sift through the genetic information of 250,000 microorganisms. This was made possible by the sequencing, in recent years, of the genomes of hundreds of thousands of bacteria and other microorganisms – including many never studied before.

These analyses revealed five groups of bacteria with alternative genetic codes, all of which had changes in the codons that code for arginine. Amongst these, four had genomes with a low proportion of guanine and cytosine nucleotides. This may have made some guanine and cytosine-rich arginine codons very rare in these organisms and, therefore, easier to be reassigned to encode another amino acid.

The work by Shulgina and Eddy demonstrates that Codetta is a new, useful tool that scientists can use to understand how genetic codes evolve. In addition, it can also help to ensure the accuracy of widely used protein databases, which assume which genetic code organisms use to predict protein sequences from their genomes.

New tyrosinases with putative action against contaminants of emerging concern

Tyrosinases (EC 1.14.18.1) are type-3 copper metalloenzymes with strong oxidative capacities and low allosteric selectivity to phenolic and non-phenolic aromatic compounds, which have been used as biosensors and biocatalysts to mitigate the impacts of environmental contaminants over aquatic ecosystems. However, the widespread use of these polyphenol oxidases is limited by elevated production costs and restricted knowledge on their spectrum of action. Here, six tyrosinase homologs were identified and characterized from the genomes of four widespread freshwater ciliates using bioinformatics. Next, we performed a virtual screening to calculate binding energies between 3D models of these homologs and ~ 1000 contaminants of emerging concern (CECs), as an indirect approach to identify likely and unlikely targets for tyrosinases. Many fine chemicals, pharmaceuticals, personal care products, illicit drugs, natural toxins, and pesticides exhibited strong binding energies to these new tyrosinases, suggesting the spectrum of targets of these enzymes might be considerably broader than previously thought. Many ciliates, including those carrying tyrosinase genes, are fast-growing unicellular microeukaryotes that can be efficiently cultured, at large scales, under in vitro conditions, suggesting these organisms should be regarded as potential low-cost sources of new environmental biotechnological molecules.

Evolution of the genetic code in the mitochondria of Labyrinthulea (Stramenopiles)

Mitochondrial translation often exhibits departures from the standard genetic code, but the full spectrum of these changes has certainly not yet been described and the molecular mechanisms behind the changes in codon meaning are rarely studied. Here we report a detailed analysis of the mitochondrial genetic code in the stramenopile group Labyrinthulea (Labyrinthulomycetes) and their relatives. In the genus Aplanochytrium, UAG is not a termination codon but encodes tyrosine, in contrast to the unaffected meaning of the UAA codon. This change is evolutionarily independent of the reassignment of both UAG and UAA as tyrosine codons recently reported from two uncultivated labyrinthuleans (S2 and S4), which we show are not thraustochytrids as proposed before, but represent the clade LAB14 previously recognised in environmental 18S rRNA gene surveys. We provide rigorous evidence that the UUA codon in the mitochondria of all labyrinthuleans serves as a termination codon instead of encoding leucine, and propose that a sense-to-stop reassignment has also affected the AGG and AGA codons in the LAB14 clade. The distribution of the different forms of sense-to-stop and stop-to-sense reassignments correlates with specific modifications of the mitochondrial release factor mtRF2a in different subsets of labyrinthuleans, and with the unprecedented loss of mtRF1a in Aplanochytrium and perhaps also in the LAB14 clade, pointing towards a possible mechanistic basis of the code changes observed. Curiously, we show that labyrinthulean mitochondria also exhibit a sense-to-sense codon reassignment, manifested as AUA encoding methionine instead of isoleucine. Furthermore, we show that this change evolved independently in the uncultivated stramenopile lineage MAST8b, together with the reassignment of the AGR codons from arginine to serine. Altogether, our study has uncovered novel variants of the mitochondrial genetic code and previously unknown modifications of the mitochondrial translation machinery, further enriching our understanding of the rules governing the evolution of one of the central molecular process in the cell.

Evolution and Unprecedented Variants of the Mitochondrial Genetic Code in a Lineage of Green Algae

Mitochondria of diverse eukaryotes have evolved various departures from the standard genetic code, but the breadth of possible modifications and their phylogenetic distribution are known only incompletely. Furthermore, it is possible that some codon reassignments in previously sequenced mitogenomes have been missed, resulting in inaccurate protein sequences in databases. Here we show, considering the distribution of codons at conserved amino acid positions in mitogenome-encoded proteins, that mitochondria of the green algal order Sphaeropleales exhibit a diversity of codon reassignments, including previously missed ones and some that are unprecedented in any translation system examined so far, necessitating redefinition of existing translation tables and creating at least seven new ones. We resolve a previous controversy concerning the meaning the UAG codon in Hydrodictyaceae, which beyond any doubt encodes alanine. We further demonstrate that AGG, sometimes together with AGA, encodes alanine instead of arginine in diverse sphaeroplealeans. Further newly detected changes include Arg-to-Met reassignment of the AGG codon and Arg-to-Leu reassignment of the CGG codon in particular species. Analysis of tRNAs specified by sphaeroplealean mitogenomes provides direct support for and molecular underpinning of the proposed reassignments. Furthermore, we point to unique mutations in the mitochondrial release factor mtRF1a that correlate with changes in the use of termination codons in Sphaeropleales, including the two independent stop-to-sense UAG reassignments, the reintroduction of UGA in some Scenedesmaceae, and the sense-to-stop reassignment of UCA widespread in the group. Codon disappearance seems to be the main drive of the dynamic evolution of the mitochondrial genetic code in Sphaeropleales.

Clades of huge phages from across Earth's ecosystems

Bacteriophages typically have small genomes1 and depend on their bacterial hosts for replication2. Here we sequenced DNA from diverse ecosystems and found hundreds of phage genomes with lengths of more than 200 kilobases (kb), including a genome of 735 kb, which is-to our knowledge-the largest phage genome to be described to date. Thirty-five genomes were manually curated to completion (circular and no gaps). Expanded genetic repertoires include diverse and previously undescribed CRISPR-Cas systems, transfer RNAs (tRNAs), tRNA synthetases, tRNA-modification enzymes, translation-initiation and elongation factors, and ribosomal proteins. The CRISPR-Cas systems of phages have the capacity to silence host transcription factors and translational genes, potentially as part of a larger interaction network that intercepts translation to redirect biosynthesis to phage-encoded functions. In addition, some phages may repurpose bacterial CRISPR-Cas systems to eliminate competing phages. We phylogenetically define the major clades of huge phages from human and other animal microbiomes, as well as from oceans, lakes, sediments, soils and the built environment. We conclude that the large gene inventories of huge phages reflect a conserved biological strategy, and that the phages are distributed across a broad bacterial host range and across Earth's ecosystems.

Conserved novel ORFs in the mitochondrial genome of the ctenophore Beroe forskalii

To date, five ctenophore species' mitochondrial genomes have been sequenced, and each contains open reading frames (ORFs) that if translated have no identifiable orthologs. ORFs with no identifiable orthologs are called unidentified reading frames (URFs). If truly protein-coding, ctenophore mitochondrial URFs represent a little understood path in early-diverging metazoan mitochondrial evolution and metabolism. We sequenced and annotated the mitochondrial genomes of three individuals of the beroid ctenophore Beroe forskalii and found that in addition to sharing the same canonical mitochondrial genes as other ctenophores, the B. forskalii mitochondrial genome contains two URFs. These URFs are conserved among the three individuals but not found in other sequenced species. We developed computational tools called pauvre and cuttlery to determine the likelihood that URFs are protein coding. There is evidence that the two URFs are under negative selection, and a novel Bayesian hypothesis test of trinucleotide frequency shows that the URFs are more similar to known coding genes than noncoding intergenic sequence. Protein structure and function prediction of all ctenophore URFs suggests that they all code for transmembrane transport proteins. These findings, along with the presence of URFs in other sequenced ctenophore mitochondrial genomes, suggest that ctenophores may have uncharacterized transmembrane proteins present in their mitochondria.

Coding palindromes in mitochondrial genes of Nematomorpha

Inverted repeats are common DNA elements, but they rarely overlap with protein-coding sequences due to the ensuing conflict with the structure and function of the encoded protein. We discovered numerous perfect inverted repeats of considerable length (up to 284 bp) embedded within the protein-coding genes in mitochondrial genomes of four Nematomorpha species. Strikingly, both arms of the inverted repeats encode conserved regions of the amino acid sequence. We confirmed enzymatic activity of the respiratory complex I encoded by inverted repeat-containing genes. The nucleotide composition of inverted repeats suggests strong selection at the amino acid level in these regions. We conclude that the inverted repeat-containing genes are transcribed and translated into functional proteins. The survey of available mitochondrial genomes reveals that several other organisms possess similar albeit shorter embedded repeats. Mitochondrial genomes of Nematomorpha demonstrate an extraordinary evolutionary compromise where protein function and stringent secondary structure elements within the coding regions are preserved simultaneously.

Rapid Genetic Code Evolution in Green Algal Mitochondrial Genomes

Genetic code deviations involving stop codons have been previously reported in mitochondrial genomes of several green plants (Viridiplantae), most notably chlorophyte algae (Chlorophyta). However, as changes in codon recognition from one amino acid to another are more difficult to infer, such changes might have gone unnoticed in particular lineages with high evolutionary rates that are otherwise prone to codon reassignments. To gain further insight into the evolution of the mitochondrial genetic code in green plants, we have conducted an in-depth study across mtDNAs from 51 green plants (32 chlorophytes and 19 streptophytes). Besides confirming known stop-to-sense reassignments, our study documents the first cases of sense-to-sense codon reassignments in Chlorophyta mtDNAs. In several Sphaeropleales, we report the decoding of AGG codons (normally arginine) as alanine, by tRNA(CCU) of various origins that carry the recognition signature for alanine tRNA synthetase. In Chromochloris, we identify tRNA variants decoding AGG as methionine and the synonymous codon CGG as leucine. Finally, we find strong evidence supporting the decoding of AUA codons (normally isoleucine) as methionine in Pycnococcus. Our results rely on a recently developed conceptual framework (CoreTracker) that predicts codon reassignments based on the disparity between DNA sequence (codons) and the derived protein sequence. These predictions are then validated by an evaluation of tRNA phylogeny, to identify the evolution of new tRNAs via gene duplication and loss, and structural modifications that lead to the assignment of new tRNA identities and a change in the genetic code.

A precedented nuclear genetic code with all three termination codons reassigned as sense codons in the syndinean Amoebophrya sp. ex Karlodinium veneficum

Amoebophrya is part of an enigmatic, diverse, and ubiquitous marine alveolate lineage known almost entirely from anonymous environmental sequencing. Two cultured Amoebophrya strains grown on core dinoflagellate hosts were used for transcriptome sequencing. BLASTx using different genetic codes suggests that Amoebophyra sp. ex Karlodinium veneficum uses the three typical stop codons (UAA, UAG, and UGA) to encode amino acids. When UAA and UAG are translated as glutamine about half of the alignments have better BLASTx scores, and when UGA is translated as tryptophan one fifth have better scores. However, the sole stop codon appears to be UGA based on conserved genes, suggesting contingent translation of UGA. Neither host sequences, nor sequences from the second strain, Amoebophrya sp. ex Akashiwo sanguinea had similar results in BLASTx searches. A genome survey of Amoebophyra sp. ex K. veneficum showed no evidence for transcript editing aside from mitochondrial transcripts. The dynein heavy chain (DHC) gene family was surveyed and of 14 transcripts only two did not use UAA, UAG, or UGA in a coding context. Overall the transcriptome displayed strong bias for A or U in third codon positions, while the tRNA genome survey showed bias against codons ending in U, particularly for amino acids with two codons ending in either C or U. Together these clues suggest contingent translation mechanisms in Amoebophyra sp. ex K. veneficum and a phylogenetically distinct instance of genetic code modification.

Megaphages infect Prevotella and variants are widespread in gut microbiomes

Bacteriophages (phages) dramatically shape microbial community composition, redistribute nutrients via host lysis and drive evolution through horizontal gene transfer. Despite their importance, much remains to be learned about phages in the human microbiome. We investigated the gut microbiomes of humans from Bangladesh and Tanzania, two African baboon social groups and Danish pigs; many of these microbiomes contain phages belonging to a clade with genomes >540 kilobases in length, the largest yet reported in the human microbiome and close to the maximum size ever reported for phages. We refer to these as Lak phages. CRISPR spacer targeting indicates that Lak phages infect bacteria of the genus Prevotella. We manually curated to completion 15 distinct Lak phage genomes recovered from metagenomes. The genomes display several interesting features, including use of an alternative genetic code, large intergenic regions that are highly expressed and up to 35 putative transfer RNAs, some of which contain enigmatic introns. Different individuals have distinct phage genotypes, and shifts in variant frequencies over consecutive sampling days reflect changes in the relative abundance of phage subpopulations. Recent homologous recombination has resulted in extensive genome admixture of nine baboon Lak phage populations. We infer that Lak phages are widespread in gut communities that contain the Prevotella species, and conclude that megaphages, with fascinating and underexplored biology, may be common but largely overlooked components of human and animal gut microbiomes.

CoreTracker: accurate codon reassignment prediction, applied to mitochondrial genomes

Motivation

Codon reassignments have been reported across all domains of life. With the increasing number of sequenced genomes, the development of systematic approaches for genetic code detection is essential for accurate downstream analyses. Three automated prediction tools exist so far: FACIL, GenDecoder and Bagheera; the last two respectively restricted to metazoan mitochondrial genomes and CUG reassignments in yeast nuclear genomes. These tools can only analyze a single genome at a time and are often not followed by a validation procedure, resulting in a high rate of false positives.

Results

We present CoreTracker, a new algorithm for the inference of sense-to-sense codon reassignments. CoreTracker identifies potential codon reassignments in a set of related genomes, then uses statistical evaluations and a random forest classifier to predict those that are the most likely to be correct. Predicted reassignments are then validated through a phylogeny-aware step that evaluates the impact of the new genetic code on the protein alignment. Handling simultaneously a set of genomes in a phylogenetic framework, allows tracing back the evolution of each reassignment, which provides information on its underlying mechanism. Applied to metazoan and yeast genomes, CoreTracker significantly outperforms existing methods on both precision and sensitivity.

Availability and implementation

CoreTracker is written in Python and available at https://github.com/UdeM-LBIT/CoreTracker.

Contact

mabrouk@iro.umontreal.ca.

Supplementary information

Supplementary data are available at Bioinformatics online.

Mitochondrial Genomes of Kinorhyncha: trnM Duplication and New Gene Orders within Animals

Many features of mitochondrial genomes of animals, such as patterns of gene arrangement, nucleotide content and substitution rate variation are extensively used in evolutionary and phylogenetic studies. Nearly 6,000 mitochondrial genomes of animals have already been sequenced, covering the majority of animal phyla. One of the groups that escaped mitogenome sequencing is phylum Kinorhyncha-an isolated taxon of microscopic worm-like ecdysozoans. The kinorhynchs are thought to be one of the early-branching lineages of Ecdysozoa, and their mitochondrial genomes may be important for resolving evolutionary relations between major animal taxa. Here we present the results of sequencing and analysis of mitochondrial genomes from two members of Kinorhyncha, Echinoderes svetlanae (Cyclorhagida) and Pycnophyes kielensis (Allomalorhagida). Their mitochondrial genomes are circular molecules approximately 15 Kbp in size. The kinorhynch mitochondrial gene sequences are highly divergent, which precludes accurate phylogenetic inference. The mitogenomes of both species encode a typical metazoan complement of 37 genes, which are all positioned on the major strand, but the gene order is distinct and unique among Ecdysozoa or animals as a whole. We predict four types of start codons for protein-coding genes in E. svetlanae and five in P. kielensis with a consensus DTD in single letter code. The mitochondrial genomes of E. svetlanae and P. kielensis encode duplicated methionine tRNA genes that display compensatory nucleotide substitutions. Two distant species of Kinorhyncha demonstrate similar patterns of gene arrangements in their mitogenomes. Both genomes have duplicated methionine tRNA genes; the duplication predates the divergence of two species. The kinorhynchs share a few features pertaining to gene order that align them with Priapulida. Gene order analysis reveals that gene arrangement specific of Priapulida may be ancestral for Scalidophora, Ecdysozoa, and even Protostomia.

Genetic Codes with No Dedicated Stop Codon: Context-Dependent Translation Termination

The prevailing view of the nuclear genetic code is that it is largely frozen and unambiguous. Flexibility in the nuclear genetic code has been demonstrated in ciliates that reassign standard stop codons to amino acids, resulting in seven variant genetic codes, including three previously undescribed ones reported here. Surprisingly, in two of these species, we find efficient translation of all 64 codons as standard amino acids and recognition of either one or all three stop codons. How, therefore, does the translation machinery interpret a “stop” codon? We provide evidence, based on ribosomal profiling and “stop” codon depletion shortly before coding sequence ends, that mRNA 3′ ends may contribute to distinguishing stop from sense in a context-dependent manner. We further propose that such context-dependent termination/readthrough suppression near transcript ends enables genetic code evolution.

•

Alternative nuclear genetic codes continue to be discovered in ciliates

•

Genetic codes with stops and all their codons encoding standard amino acids exist

•

Transcript ends may distinguish stop codons as such in ambiguous genetic codes

•

The ability to resolve genetic code ambiguity may enable genetic code evolution

Phylogenomic analysis of Candidatus 'Izimaplasma' species: free-living representatives from a Tenericutes clade found in methane seeps

Tenericutes are a unique class of bacteria that lack a cell wall and are typically parasites or commensals of eukaryotic hosts. Environmental 16S rDNA surveys have identified a number of tenericute clades in diverse environments, introducing the possibility that these Tenericutes may represent non-host-associated, free-living microorganisms. Metagenomic sequencing of deep-sea methane seep sediments resulted in the assembly of two genomes from a Tenericutes-affiliated clade currently known as 'NB1-n' (SILVA taxonomy) or 'RF3' (Greengenes taxonomy). Metabolic reconstruction revealed that, like cultured members of the Mollicutes, these 'NB1-n' representatives lack a tricarboxylic acid cycle and instead use anaerobic fermentation of simple sugars for substrate level phosphorylation. Notably, the genomes also contained a number of unique metabolic features including hydrogenases and a simplified electron transport chain containing an RNF complex, cytochrome bd oxidase and complex I. On the basis of the metabolic potential predicted from the annotated genomes, we devised an anaerobic enrichment media that stimulated the growth of these Tenericutes at 10 °C, resulting in a mixed culture where these organisms represented ~60% of the total cells by targeted fluorescence in situ hybridization (FISH). Visual identification by FISH confirmed these organisms were not directly associated with Eukaryotes and electron cryomicroscopy of cells in the enrichment culture confirmed an ultrastructure consistent with the defining phenotypic property of Tenericutes, with a single membrane and no cell wall. On the basis of their unique gene content, phylogenetic placement and ultrastructure, we propose these organisms represent a novel class within the Tenericutes, and suggest the names Candidatus 'Izimaplasma sp. HR1' and Candidatus 'Izimaplasma sp. HR2' for the two genome representatives.

Genome-based microbial ecology of anammox granules in a full-scale wastewater treatment system

Partial-nitritation anammox (PNA) is a novel wastewater treatment procedure for energy-efficient ammonium removal. Here we use genome-resolved metagenomics to build a genome-based ecological model of the microbial community in a full-scale PNA reactor. Sludge from the bioreactor examined here is used to seed reactors in wastewater treatment plants around the world; however, the role of most of its microbial community in ammonium removal remains unknown. Our analysis yielded 23 near-complete draft genomes that together represent the majority of the microbial community. We assign these genomes to distinct anaerobic and aerobic microbial communities. In the aerobic community, nitrifying organisms and heterotrophs predominate. In the anaerobic community, widespread potential for partial denitrification suggests a nitrite loop increases treatment efficiency. Of our genomes, 19 have no previously cultivated or sequenced close relatives and six belong to bacterial phyla without any cultivated members, including the most complete Omnitrophica (formerly OP3) genome to date.

ANaerobic AMMonium OXidation (ANAMMOX) combined with partial nitritation has been adopted for removal of ammonium from wastewater. Here, Speth et al. describe the bacterial metagenome of a partial-nitritation/anammox (PNA) reactor, and provide 23 draft genomes, 19 of which were previously uncharacterized/sequenced/cultivated.

Recoding of the stop codon UGA to glycine by a BD1-5/SN-2 bacterium and niche partitioning between Alpha- and Gammaproteobacteria in a tidal sediment microbial community naturally selected in a laboratory chemostat

Sandy coastal sediments are global hotspots for microbial mineralization of organic matter and denitrification. These sediments are characterized by advective porewater flow, tidal cycling and an active and complex microbial community. Metagenomic sequencing of microbial communities sampled from such sediments showed that potential sulfur oxidizing Gammaproteobacteria and members of the enigmatic BD1-5/SN-2 candidate phylum were abundant in situ (>10% and ~2% respectively). By mimicking the dynamic oxic/anoxic environmental conditions of the sediment in a laboratory chemostat, a simplified microbial community was selected from the more complex inoculum. Metagenomics, proteomics and fluorescence in situ hybridization showed that this simplified community contained both a potential sulfur oxidizing Gammaproteobacteria (at 24 ± 2% abundance) and a member of the BD1-5/SN-2 candidate phylum (at 7 ± 6% abundance). Despite the abundant supply of organic substrates to the chemostat, proteomic analysis suggested that the selected gammaproteobacterium grew partially autotrophically and performed hydrogen/formate oxidation. The enrichment of a member of the BD1-5/SN-2 candidate phylum enabled, for the first time, direct microscopic observation by fluorescent in situ hybridization and the experimental validation of the previously predicted translation of the stop codon UGA into glycine.

Reference-independent comparative metagenomics using cross-assembly: crAss

MOTIVATION

Metagenomes are often characterized by high levels of unknown sequences. Reads derived from known microorganisms can easily be identified and analyzed using fast homology search algorithms and a suitable reference database, but the unknown sequences are often ignored in further analyses, biasing conclusions. Nevertheless, it is possible to use more data in a comparative metagenomic analysis by creating a cross-assembly of all reads, i.e. a single assembly of reads from different samples. Comparative metagenomics studies the interrelationships between metagenomes from different samples. Using an assembly algorithm is a fast and intuitive way to link (partially) homologous reads without requiring a database of reference sequences.

RESULTS

Here, we introduce crAss, a novel bioinformatic tool that enables fast simple analysis of cross-assembly files, yielding distances between all metagenomic sample pairs and an insightful image displaying the similarities.

Data mining in the Life Sciences with Random Forest: a walk in the park or lost in the jungle?

In the Life Sciences 'omics' data is increasingly generated by different high-throughput technologies. Often only the integration of these data allows uncovering biological insights that can be experimentally validated or mechanistically modelled, i.e. sophisticated computational approaches are required to extract the complex non-linear trends present in omics data. Classification techniques allow training a model based on variables (e.g. SNPs in genetic association studies) to separate different classes (e.g. healthy subjects versus patients). Random Forest (RF) is a versatile classification algorithm suited for the analysis of these large data sets. In the Life Sciences, RF is popular because RF classification models have a high-prediction accuracy and provide information on importance of variables for classification. For omics data, variables or conditional relations between variables are typically important for a subset of samples of the same class. For example: within a class of cancer patients certain SNP combinations may be important for a subset of patients that have a specific subtype of cancer, but not important for a different subset of patients. These conditional relationships can in principle be uncovered from the data with RF as these are implicitly taken into account by the algorithm during the creation of the classification model. This review details some of the to the best of our knowledge rarely or never used RF properties that allow maximizing the biological insights that can be extracted from complex omics data sets using RF.

Evolution and diversification of the organellar release factor family

Translation termination is accomplished by proteins of the Class I release factor family (RF) that recognize stop codons and catalyze the ribosomal release of the newly synthesized peptide. Bacteria have two canonical RFs: RF1 recognizes UAA and UAG, RF2 recognizes UAA and UGA. Despite that these two release factor proteins are sufficient for de facto translation termination, the eukaryotic organellar RF protein family, which has evolved from bacterial release factors, has expanded considerably, comprising multiple subfamilies, most of which have not been functionally characterized or formally classified. Here, we integrate multiple sources of information to analyze the remarkable differentiation of the RF family among organelles. We document the origin, phylogenetic distribution and sequence structure features of the mitochondrial and plastidial release factors: mtRF1a, mtRF1, mtRF2a, mtRF2b, mtRF2c, ICT1, C12orf65, pRF1, and pRF2, and review published relevant experimental data. The canonical release factors (mtRF1a, mtRF2a, pRF1, and pRF2) and ICT1 are derived from bacterial ancestors, whereas the others have resulted from gene duplications of another release factor. These new RF family members have all lost one or more specific motifs relevant for bona fide release factor function but are mostly targeted to the same organelle as their ancestor. We also characterize the subset of canonical release factor proteins that bear nonclassical PxT/SPF tripeptide motifs and provide a molecular-model-based rationale for their retained ability to recognize stop codons. Finally, we analyze the coevolution of canonical RFs with the organellar genetic code. Although the RF presence in an organelle and its stop codon usage tend to coevolve, we find three taxa that encode an RF2 without using UGA stop codons, and one reverse scenario, where mamiellales green algae use UGA stop codons in their mitochondria without having a mitochondrial type RF2. For the latter, we put forward a “stop-codon reinvention” hypothesis that involves the retargeting of the plastid release factor to the mitochondrion.