Origins of genome-editing excisases as illuminated by the somatic genome of the ciliate Blepharisma.

Massive DNA excision occurs regularly in ciliates, ubiquitous microbial eukaryotes with somatic and germline nuclei in the same cell. Tens of thousands of internally eliminated sequences (IESs) scattered throughout the ciliate germline genome are deleted during the development of the streamlined somatic genome. The genus Blepharisma represents one of the two high-level ciliate clades (subphylum Postciliodesmatophora) and, unusually, has dual pathways of somatic nuclear and genome development. This makes it ideal for investigating the functioning and evolution of these processes. Here we report the somatic genome assembly of Blepharisma stoltei strain ATCC 30299 (41 Mbp), arranged as numerous telomere-capped minichromosomal isoforms. This genome encodes eight PiggyBac transposase homologs no longer harbored by transposons. All appear subject to purifying selection, but just one, the putative IES excisase, has a complete catalytic triad. We hypothesize that PiggyBac homologs were ancestral excisases that enabled the evolution of extensive natural genome editing.

MITE infestation accommodated by genome editing in the germline genome of the ciliate Blepharisma.

During their development following sexual conjugation, ciliates excise numerous internal eliminated sequences (IESs) from a copy of the germline genome to produce the functional somatic genome. Most IESs are thought to have originated from transposons, but the presumed homology is often obscured by sequence decay. To obtain more representative perspectives on the nature of IESs and ciliate genome editing, we assembled 40,000 IESs of Blepharisma stoltei, a species belonging to a lineage (Heterotrichea) that diverged early from those of the intensively studied model ciliate species. About a quarter of IESs were short (<115 bp), largely nonrepetitive, and with a pronounced ~10 bp periodicity in length; the remainder were longer (up to 7 kbp) and nonperiodic and contained abundant interspersed repeats. Contrary to the expectation from current models, the assembled Blepharisma germline genome encodes few transposases. Instead, its most abundant repeat (8,000 copies) is a Miniature Inverted-repeat Transposable Element (MITE), apparently a deletion derivative of a germline-limited Pogo-family transposon. We hypothesize that MITEs are an important source of IESs whose proliferation is eventually self-limiting and that rather than defending the germline genomes against mobile elements, transposase domestication actually facilitates the accumulation of junk DNA.

BleTIES: annotation of natural genome editing in ciliates using long read sequencing.

SUMMARY: Ciliates are single-celled eukaryotes that eliminate specific, interspersed DNA sequences (internally eliminated sequences, IESs) from their genomes during development. These are challenging to annotate and assemble because IES-containing sequences are typically much less abundant in the cell than those without, and IES sequences themselves often contain repetitive and low-complexity sequences. Long-read sequencing technologies from Pacific Biosciences and Oxford Nanopore have the potential to reconstruct longer IESs than has been possible with short reads but require a different assembly strategy. Here we present BleTIES, a software toolkit for detecting, assembling, and analyzing IESs using mapped long reads. AVAILABILITY AND IMPLEMENTATION: BleTIES is implemented in Python 3. Source code is available at https://github.com/Swart-lab/bleties (MIT license) and also distributed via Bioconda. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Determination of Abundant Metabolite Matrix Adducts Illuminates the Dark Metabolome of MALDI-Mass Spectrometry Imaging Datasets.

Spatial metabolomics using mass spectrometry imaging (MSI) is a powerful tool to map hundreds to thousands of metabolites in biological systems. One major challenge in MSI is the annotation of m/z values, which is substantially complicated by background ions introduced throughout the chemicals and equipment used during experimental procedures. Among many factors, the formation of adducts with sodium or potassium ions, or in case of matrix-assisted laser desorption ionization (MALDI)-MSI, the presence of abundant matrix clusters strongly increases total m/z peak counts. Currently, there is a limitation to identify the chemistry of the many unknown peaks to interpret their biological function. We took advantage of the co-localization of adducts with their parent ions and the accuracy of high mass resolution to estimate adduct abundance in 20 datasets from different vendors of mass spectrometers. Metabolites ranging from lipids to amines and amino acids form matrix adducts with the commonly used 2,5-dihydroxybenzoic acid (DHB) matrix like [M + (DHB-H2O) + H]+ and [M + DHB + Na]+. Current data analyses neglect those matrix adducts and overestimate total metabolite numbers, thereby expanding the number of unidentified peaks. Our study demonstrates that MALDI-MSI data are strongly influenced by adduct formation across different sample types and vendor platforms and reveals a major influence of so far unrecognized metabolite-matrix adducts on total peak counts (up to one third). We developed a software package, mass2adduct, for the community for an automated putative assignment and quantification of metabolite-matrix adducts enabling users to ultimately focus on the biologically relevant portion of the MSI data.

Desulfovibrio diazotrophicus sp. nov., a sulfate-reducing bacterium from the human gut capable of nitrogen fixation.

Sulfate-reducing bacteria (SRB) are widespread in human guts, yet their expansion has been linked to colonic diseases. We report the isolation, sequencing and physiological characterization of strain QI0027T , a novel SRB species belonging to the class Desulfovibrionia. Metagenomic sequencing of stool samples from 45 Chinese individuals, and comparison with 1690 Desulfovibrionaceae metagenome-assembled genomes recovered from humans of diverse geographic locations, revealed the presence of QI0027T in 22 further individuals. QI0027T encoded nitrogen fixation genes and based on the acetylene reduction assay, actively fixed nitrogen. Transcriptomics revealed that QI0027T overexpressed 42 genes in nitrogen-limiting conditions compared to cultures supplemented with ammonia, including genes encoding nitrogenases, a urea uptake system and the urease complex. Reanalyses of 835 public stool metatranscriptomes showed that nitrogenase genes from Desulfovibrio bacteria were expressed in six samples suggesting that nitrogen fixation might be active in the gut environment. Although frequently thought of as a nutrient-rich environment, nitrogen fixation can occur in the human gut. Animals are often nitrogen limited and have evolved diverse strategies to capture biologically active nitrogen, ranging from amino acid transporters to stable associations with beneficial microbes that provide fixed nitrogen. QI0027T is the first Desulfovibrio human isolate for which nitrogen fixation has been demonstrated, suggesting that some sulfate-reducing bacteria could also play a role in the availability of nitrogen in the gut.

phyloFlash: Rapid Small-Subunit rRNA Profiling and Targeted Assembly from Metagenomes.

The small-subunit rRNA (SSU rRNA) gene is the key marker in molecular ecology for all domains of life, but it is largely absent from metagenome-assembled genomes that often are the only resource available for environmental microbes. Here, we present phyloFlash, a pipeline to overcome this gap with rapid, SSU rRNA-centered taxonomic classification, targeted assembly, and graph-based binning of full metagenomic assemblies. We show that a cleanup of artifacts is pivotal even with a curated reference database. With such a filtered database, the general-purpose mapper BBmap extracts SSU rRNA reads five times faster than the rRNA-specialized tool SortMeRNA with similar sensitivity and higher selectivity on simulated metagenomes. Reference-based targeted assemblers yielded either highly fragmented assemblies or high levels of chimerism, so we employ the general-purpose genomic assembler SPAdes. Our optimized implementation is independent of reference database composition and has satisfactory levels of chimera formation. phyloFlash quickly processes Illumina (meta)genomic data, is straightforward to use, even as part of high-throughput quality control, and has user-friendly output reports. The software is available at https://github.com/HRGV/phyloFlash (GPL3 license) and is documented with an online manual.IMPORTANCE To track organisms across all domains of life, the SSU rRNA gene is the gold standard. Many environmental microbes are known only from high-throughput sequence data, but the SSU rRNA gene, the key to visualization by molecular probes and link to existing literature, is often missing from metagenome-assembled genomes (MAGs). The easy-to-use phyloFlash software suite tackles this gap with rapid, SSU rRNA-centered taxonomic classification, targeted assembly, and graph-based linking to MAGs. Starting from a cleaned reference database, phyloFlash profiles the taxonomic diversity and assembles the sorted SSU rRNA reads. The phyloFlash design is domain agnostic and covers eukaryotes, archaea, and bacteria alike. phyloFlash also provides utilities to visualize multisample comparisons and to integrate the recovered SSU rRNAs in a metagenomics workflow by linking them to MAGs using assembly graph parsing.

Sulfur-Oxidizing Symbionts without Canonical Genes for Autotrophic CO2 Fixation.

Since the discovery of symbioses between sulfur-oxidizing (thiotrophic) bacteria and invertebrates at hydrothermal vents over 40 years ago, it has been assumed that autotrophic fixation of CO2 by the symbionts drives these nutritional associations. In this study, we investigated "Candidatus Kentron," the clade of symbionts hosted by Kentrophoros, a diverse genus of ciliates which are found in marine coastal sediments around the world. Despite being the main food source for their hosts, Kentron bacteria lack the key canonical genes for any of the known pathways for autotrophic carbon fixation and have a carbon stable isotope fingerprint that is unlike other thiotrophic symbionts from similar habitats. Our genomic and transcriptomic analyses instead found metabolic features consistent with growth on organic carbon, especially organic and amino acids, for which they have abundant uptake transporters. All known thiotrophic symbionts have converged on using reduced sulfur to gain energy lithotrophically, but they are diverse in their carbon sources. Some clades are obligate autotrophs, while many are mixotrophs that can supplement autotrophic carbon fixation with heterotrophic capabilities similar to those in Kentron. Here we show that Kentron bacteria are the only thiotrophic symbionts that appear to be entirely heterotrophic, unlike all other thiotrophic symbionts studied to date, which possess either the Calvin-Benson-Bassham or the reverse tricarboxylic acid cycle for autotrophy.IMPORTANCE Many animals and protists depend on symbiotic sulfur-oxidizing bacteria as their main food source. These bacteria use energy from oxidizing inorganic sulfur compounds to make biomass autotrophically from CO2, serving as primary producers for their hosts. Here we describe a clade of nonautotrophic sulfur-oxidizing symbionts, "Candidatus Kentron," associated with marine ciliates. They lack genes for known autotrophic pathways and have a carbon stable isotope fingerprint heavier than other symbionts from similar habitats. Instead, they have the potential to oxidize sulfur to fuel the uptake of organic compounds for heterotrophic growth, a metabolic mode called chemolithoheterotrophy that is not found in other symbioses. Although several symbionts have heterotrophic features to supplement primary production, in Kentron they appear to supplant it entirely.

Chemosynthetic symbiont with a drastically reduced genome serves as primary energy storage in the marine flatworm Paracatenula.

Hosts of chemoautotrophic bacteria typically have much higher biomass than their symbionts and consume symbiont cells for nutrition. In contrast to this, chemoautotrophic Candidatus Riegeria symbionts in mouthless Paracatenula flatworms comprise up to half of the biomass of the consortium. Each species of Paracatenula harbors a specific Ca Riegeria, and the endosymbionts have been vertically transmitted for at least 500 million years. Such prolonged strict vertical transmission leads to streamlining of symbiont genomes, and the retained physiological capacities reveal the functions the symbionts provide to their hosts. Here, we studied a species of Paracatenula from Sant'Andrea, Elba, Italy, using genomics, gene expression, imaging analyses, as well as targeted and untargeted MS. We show that its symbiont, Ca R. santandreae has a drastically smaller genome (1.34 Mb) than the symbiont´s free-living relatives (4.29-4.97 Mb) but retains a versatile and energy-efficient metabolism. It encodes and expresses a complete intermediary carbon metabolism and enhanced carbon fixation through anaplerosis and accumulates massive intracellular inclusions such as sulfur, polyhydroxyalkanoates, and carbohydrates. Compared with symbiotic and free-living chemoautotrophs, Ca R. santandreae's versatility in energy storage is unparalleled in chemoautotrophs with such compact genomes. Transmission EM as well as host and symbiont expression data suggest that Ca R. santandreae largely provisions its host via outer-membrane vesicle secretion. With its high share of biomass in the symbiosis and large standing stocks of carbon and energy reserves, it has a unique role for bacterial symbionts-serving as the primary energy storage for its animal host.

Author Correction: Chemosynthetic symbionts of marine invertebrate animals are capable of nitrogen fixation.

In this Article, the completeness and number of contigs for draft genomes from two individuals of Laxus oneistus are incorrect in the main text, although the correct information is included in Table 1. The original and corrected versions of the relevant sentence are shown in the correction notice.

Specificity in diversity: single origin of a widespread ciliate-bacteria symbiosis.

Symbioses between eukaryotes and sulfur-oxidizing (thiotrophic) bacteria have convergently evolved multiple times. Although well described in at least eight classes of metazoan animals, almost nothing is known about the evolution of thiotrophic symbioses in microbial eukaryotes (protists). In this study, we characterized the symbioses between mouthless marine ciliates of the genus Kentrophoros, and their thiotrophic bacteria, using comparative sequence analysis and fluorescence in situ hybridization. Ciliate small-subunit rRNA sequences were obtained from 17 morphospecies collected in the Mediterranean and Caribbean, and symbiont sequences from 13 of these morphospecies. We discovered a new Kentrophoros morphotype where the symbiont-bearing surface is folded into pouch-like compartments, illustrating the variability of the basic body plan. Phylogenetic analyses revealed that all investigated Kentrophoros belonged to a single clade, despite the remarkable morphological diversity of these hosts. The symbionts were also monophyletic and belonged to a new clade within the Gammaproteobacteria, with no known cultured representatives. Each host morphospecies had a distinct symbiont phylotype, and statistical analyses revealed significant support for host-symbiont codiversification. Given that these symbioses were collected from two widely separated oceans, our results indicate that symbiotic associations in unicellular hosts can be highly specific and stable over long periods of evolutionary time.

Chemosynthetic symbionts of marine invertebrate animals are capable of nitrogen fixation.

Chemosynthetic symbioses are partnerships between invertebrate animals and chemosynthetic bacteria. The latter are the primary producers, providing most of the organic carbon needed for the animal host's nutrition. We sequenced genomes of the chemosynthetic symbionts from the lucinid bivalve Loripes lucinalis and the stilbonematid nematode Laxus oneistus. The symbionts of both host species encoded nitrogen fixation genes. This is remarkable as no marine chemosynthetic symbiont was previously known to be capable of nitrogen fixation. We detected nitrogenase expression by the symbionts of lucinid clams at the transcriptomic and proteomic level. Mean stable nitrogen isotope values of Loripes lucinalis were within the range expected for fixed atmospheric nitrogen, further suggesting active nitrogen fixation by the symbionts. The ability to fix nitrogen may be widespread among chemosynthetic symbioses in oligotrophic habitats, where nitrogen availability often limits primary productivity.

High rates of microbial dinitrogen fixation and sulfate reduction associated with the Mediterranean seagrass Posidonia oceanica.

Seagrass meadows of Posidonia oceanica represent hotspots of productivity in the oligotrophic Mediterranean Sea. The lack of dissolved inorganic nitrogen (DIN) in the seawater suggests that the N-demand of these meadows might be in part supported by microbial dinitrogen (N2) fixation. However, currently there are no direct N2 fixation measurements available for this productive marine macrophyte. Here we investigated N2 fixation activity associated with P. oceanica leaf, rhizome and root pieces. In 15N2 incubations, the roots exhibited highest rates of N2 fixation. The rates varied considerably between replicates, presumably due to a patchy microbial colonization of the roots. Additions of organic carbon compounds (acetate, glucose, sucrose or algal lysate) did not enhance the N2 fixation rates. Sulfate reduction rates measured alongside were also highest in root incubations. Correspondingly, sequences of the nifH gene (a marker gene for the iron protein of the N2-fixing enzyme nitrogenase) related to known sulfate-reducing bacteria were retrieved from P. oceanica roots. Other nifH sequences clustered with known heterotrophic diazotrophs previously identified in other marine macrophytes. In particular, many sequences obtained from P. oceanica roots were similar (>94%) to a saltmarsh rhizosphere-associated heterotrophic diazotroph, indicating that heterotrophic lifestyle might be common among marine macrophyte-associated diazotrophs.

gbtools: Interactive Visualization of Metagenome Bins in R.

Improvements in DNA sequencing technology have increased the amount and quality of sequences that can be obtained from metagenomic samples, making it practical to extract individual microbial genomes from metagenomic assemblies ("binning"). However, while many tools and methods exist for unsupervised binning with various statistical algorithms, there are few options for visualizing the results, even though visualization is vital to exploratory data analysis. We have developed gbtools, a software package that allows users to visualize metagenomic assemblies by plotting coverage (sequencing depth) and GC values of contigs, and also to annotate the plots with taxonomic information. Different sets of annotations, including taxonomic assignments from conserved marker genes or SSU rRNA genes, can be imported simultaneously; users can choose which annotations to plot. Bins can be manually defined from plots, or be imported from third-party binning tools and overlaid onto plots, such that results from different methods can be compared side-by-side. gbtools reports summary statistics of bins including marker gene completeness, and allows the user to add or subtract bins with each other. We illustrate some of the functions available in gbtools with two examples: the metagenome of Olavius algarvensis, a marine oligochaete worm that has up to five bacterial symbionts, and the metagenome of a synthetic mock community comprising 64 bacterial and archaeal strains. We show how instances of poor automated binning, sequencer GC% bias, and variation between samples can be quickly diagnosed by visualization, and demonstrate how the results from different binning tools can be combined and refined to yield manually curated bins with higher completeness. gbtools is open-source and written in R. The software package, documentation, and example data are available freely online at https://github.com/kbseah/genome-bin-tools.