MarkerScan: Separation and assembly of cobionts sequenced alongside target species in biodiversity genomics projects.

Contamination of public databases by mislabelled sequences has been highlighted for many years and the avalanche of novel sequencing data now being deposited has the potential to make databases difficult to use effectively. It is therefore crucial that sequencing projects and database curators perform pre-submission checks to remove obvious contamination and avoid propagating erroneous taxonomic relationships. However, it is important also to recognise that biological contamination of a target sample with unexpected species' DNA can also lead to the discovery of fascinating biological phenomena through the identification of environmental organisms or endosymbionts. Here, we present a novel, integrated method for detection and generation of high-quality genomes of all non-target genomes co-sequenced in eukaryotic genome sequencing projects. After performing taxonomic profiling of an assembly from the raw data, and leveraging the identity of small rRNA sequences discovered therein as markers, a targeted classification approach retrieves and assembles high-quality genomes. The genomes of these cobionts are then not only removed from the target species' genome but also available for further interrogation. Source code is available from https://github.com/CobiontID/MarkerScan. MarkerScan is written in Python and is deployed as a Docker container.

This article addresses a common issue in genetic research: the accidental mixing of genetic information from different species in public databases, often due to mislabelling or contamination. Interestingly, this 'contamination' can sometimes lead to exciting discoveries, like identifying DNA from unexpected species in a sample, revealing insights about organisms that live in the environment of the target organism. In our study, we developed a tool called MarkerScan for identifying these additional species found alongside the target species in eukaryotic genome sequencing projects. The method includes a way to sequence the whole genomes of the additional species. Our method involves sorting through the genetic data to identify certain small RNA sequences, which we then use as markers. These markers help to classify and assemble high-quality genomes from these additional species. This not only cleans up the main target species' genome data but also provides new, valuable genomes for further exploration.

Phylogenomic analysis of Wolbachia genomes from the Darwin Tree of Life biodiversity genomics project.

The Darwin Tree of Life (DToL) project aims to sequence all described terrestrial and aquatic eukaryotic species found in Britain and Ireland. Reference genome sequences are generated from single individuals for each target species. In addition to the target genome, sequenced samples often contain genetic material from microbiomes, endosymbionts, parasites, and other cobionts. Wolbachia endosymbiotic bacteria are found in a diversity of terrestrial arthropods and nematodes, with supergroups A and B the most common in insects. We identified and assembled 110 complete Wolbachia genomes from 93 host species spanning 92 families by filtering data from 368 insect species generated by the DToL project. From 15 infected species, we assembled more than one Wolbachia genome, including cases where individuals carried simultaneous supergroup A and B infections. Different insect orders had distinct patterns of infection, with Lepidopteran hosts mostly infected with supergroup B, while infections in Diptera and Hymenoptera were dominated by A-type Wolbachia. Other than these large-scale order-level associations, host and Wolbachia phylogenies revealed no (or very limited) cophylogeny. This points to the occurrence of frequent host switching events, including between insect orders, in the evolutionary history of the Wolbachia pandemic. While supergroup A and B genomes had distinct GC% and GC skew, and B genomes had a larger core gene set and tended to be longer, it was the abundance of copies of bacteriophage WO who was a strong determinant of Wolbachia genome size. Mining raw genome data generated for reference genome assemblies is a robust way of identifying and analysing cobiont genomes and giving greater ecological context for their hosts.

Mitotic recombination between homologous chromosomes drives genomic diversity in diatoms.

Diatoms, an evolutionarily successful group of microalgae, display high levels of intraspecific genetic variability in natural populations. However, the contribution of various mechanisms generating such diversity is unknown. Here we estimated the genetic micro-diversity within a natural diatom population and mapped the genomic changes arising within clonally propagated diatom cell cultures. Through quantification of haplotype diversity by next-generation sequencing and amplicon re-sequencing of selected loci, we documented a rapid accumulation of multiple haplotypes accompanied by the appearance of novel protein variants in cell cultures initiated from a single founder cell. Comparison of the genomic changes between mother and daughter cells revealed copy number variation and copy-neutral loss of heterozygosity leading to the fixation of alleles within individual daughter cells. The loss of heterozygosity can be accomplished by recombination between homologous chromosomes. To test this hypothesis, we established an endogenous readout system and estimated that the frequency of interhomolog mitotic recombination was under standard growth conditions 4.2 events per 100 cell divisions. This frequency is increased under environmental stress conditions, including treatment with hydrogen peroxide and cadmium. These data demonstrate that copy number variation and mitotic recombination between homologous chromosomes underlie clonal variability in diatom populations. We discuss the potential adaptive evolutionary benefits of the plastic response in the interhomolog mitotic recombination rate, and we propose that this may have contributed to the ecological success of diatoms.

Author Correction: The Seminavis robusta genome provides insights into the evolutionary adaptations of benthic diatoms.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Comprehensive and Functional Analysis of Horizontal Gene Transfer Events in Diatoms.

Diatoms are a diverse group of mainly photosynthetic algae, responsible for 20% of worldwide oxygen production, which can rapidly respond to favorable conditions and often outcompete other phytoplankton. We investigated the contribution of horizontal gene transfer (HGT) to its ecological success. A large-scale phylogeny-based prokaryotic HGT detection procedure across nine sequenced diatoms showed that 3-5% of their proteome has a horizontal origin and a large influx occurred at the ancestor of diatoms. More than 90% of HGT genes are expressed, and species-specific HGT genes in Phaeodactylum tricornutum undergo strong purifying selection. Genes derived from HGT are implicated in several processes including environmental sensing and expand the metabolic toolbox. Cobalamin (vitamin B12) is an essential cofactor for roughly half of the diatoms and is only produced by bacteria. Five consecutive genes involved in the final synthesis of the cobalamin biosynthetic pathway, which could function as scavenging and repair genes, were detected as HGT. The full suite of these genes was detected in the cold-adapted diatom Fragilariopsis cylindrus. This might give diatoms originating from the Southern Ocean, a region typically depleted in cobalamin, a competitive advantage. Overall, we show that HGT is a prevalent mechanism that is actively used in diatoms to expand its adaptive capabilities.

The Seminavis robusta genome provides insights into the evolutionary adaptations of benthic diatoms.

Benthic diatoms are the main primary producers in shallow freshwater and coastal environments, fulfilling important ecological functions such as nutrient cycling and sediment stabilization. However, little is known about their evolutionary adaptations to these highly structured but heterogeneous environments. Here, we report a reference genome for the marine biofilm-forming diatom Seminavis robusta, showing that gene family expansions are responsible for a quarter of all 36,254 protein-coding genes. Tandem duplications play a key role in extending the repertoire of specific gene functions, including light and oxygen sensing, which are probably central for its adaptation to benthic habitats. Genes differentially expressed during interactions with bacteria are strongly conserved in other benthic diatoms while many species-specific genes are strongly upregulated during sexual reproduction. Combined with re-sequencing data from 48 strains, our results offer insights into the genetic diversity and gene functions in benthic diatoms.

Distinctive Growth and Transcriptional Changes of the Diatom Seminavis robusta in Response to Quorum Sensing Related Compounds.

In aquatic habitats, diatoms are frequently found in association with Proteobacteria, many members of which employ cell-to-cell communication via N-acyl homoserine lactones (AHLs). It has been suggested that diatoms could distinguish between beneficial and algicidal bacteria in their surroundings by sensing AHLs. Although some microalgae can interfere with AHL signaling, e.g., by releasing AHL mimics or degrading them, molecular responses to AHLs in microalgae are still unclear. Therefore, we tested the effects of short-chained AHLs, i.e., N-hexanoyl homoserine lactone (C6-HSL), N-3-hydroxyhexanoyl homoserine lactone (OH-C6-HSL), and N-3-oxohexanoyl homoserine lactone (oxo-C6-HSL) and long-chained AHLs, i.e., N-tetradecanoyl homoserine lactone (C14-HSL), N-3-hydroxytetradecanoyl homoserine lactone (OH-C14-HSL), and N-3-oxotetradecanoyl homoserine lactone (oxo-C14-HSL), on growth of the benthic diatom Seminavis robusta. All tested short-chained AHLs did not affect diatom growth, while long-chained AHLs promoted (C14-HSL) or inhibited (OH-C14-HSL and oxo-C14-HSL) growth. To investigate the physiological effects of these long-chained AHLs in more detail, an RNA-seq experiment was performed during which S. robusta was treated with the growth-promoting C14-HSL and the growth-inhibiting oxo-C14-HSL. One tetramic acid was also tested (TA14), a structural rearrangement product of oxo-C14-HSL, which also induced growth inhibition in S. robusta. After 3 days of treatment, analysis revealed that 3,410 genes were differentially expressed in response to at least one of the compounds. In the treatment with the growth-promoting C14-HSL many genes involved in intracellular signaling were upregulated. On the other hand, exposure to growth-inhibiting oxo-C14-HSL and TA14 triggered a switch in lipid metabolism towards increased fatty acid degradation. In addition, oxo-C14-HSL led to downregulation of cell cycle genes, which is in agreement with the stagnation of cell growth in this treatment. Combined, our results indicate that bacterial signaling molecules with high structural similarity induce contrasting physiological responses in S. robusta.

Virus-host coexistence in phytoplankton through the genomic lens.

Virus-microbe interactions in the ocean are commonly described by "boom and bust" dynamics, whereby a numerically dominant microorganism is lysed and replaced by a virus-resistant one. Here, we isolated a microalga strain and its infective dsDNA virus whose dynamics are characterized instead by parallel growth of both the microalga and the virus. Experimental evolution of clonal lines revealed that this viral production originates from the lysis of a minority of virus-susceptible cells, which are regenerated from resistant cells. Whole-genome sequencing demonstrated that this resistant-susceptible switch involved a large deletion on one chromosome. Mathematical modeling explained how the switch maintains stable microalga-virus population dynamics consistent with their observed growth pattern. Comparative genomics confirmed an ancient origin of this "accordion" chromosome despite a lack of sequence conservation. Together, our results show how dynamic genomic rearrangements may account for a previously overlooked coexistence mechanism in microalgae-virus interactions.

A widespread alternative squalene epoxidase participates in eukaryote steroid biosynthesis.

Steroids are essential triterpenoid molecules that are present in all eukaryotes and modulate the fluidity and flexibility of cell membranes. Steroids also serve as signalling molecules that are crucial for growth, development and differentiation of multicellular organisms1-3. The steroid biosynthetic pathway is highly conserved and is key in eukaryote evolution4-7. The flavoprotein squalene epoxidase (SQE) catalyses the first oxygenation reaction in this pathway and is rate limiting. However, despite its conservation in animals, plants and fungi, several phylogenetically widely distributed eukaryote genomes lack an SQE-encoding gene7,8. Here, we discovered and characterized an alternative SQE (AltSQE) belonging to the fatty acid hydroxylase superfamily. AltSQE was identified through screening of a gene library of the diatom Phaeodactylum tricornutum in a SQE-deficient yeast. In accordance with its divergent protein structure and need for cofactors, we found that AltSQE is insensitive to the conventional SQE inhibitor terbinafine. AltSQE is present in many eukaryotic lineages but is mutually exclusive with SQE and shows a patchy distribution within monophyletic clades. Our discovery provides an alternative element for the conserved steroid biosynthesis pathway, raises questions about eukaryote metabolic evolution and opens routes to develop selective SQE inhibitors to control hazardous organisms.

Insights into the Evolution of Multicellularity from the Sea Lettuce Genome.

We report here the 98.5 Mbp haploid genome (12,924 protein coding genes) of Ulva mutabilis, a ubiquitous and iconic representative of the Ulvophyceae or green seaweeds. Ulva's rapid and abundant growth makes it a key contributor to coastal biogeochemical cycles; its role in marine sulfur cycles is particularly important because it produces high levels of dimethylsulfoniopropionate (DMSP), the main precursor of volatile dimethyl sulfide (DMS). Rapid growth makes Ulva attractive biomass feedstock but also increasingly a driver of nuisance "green tides." Ulvophytes are key to understanding the evolution of multicellularity in the green lineage, and Ulva morphogenesis is dependent on bacterial signals, making it an important species with which to study cross-kingdom communication. Our sequenced genome informs these aspects of ulvophyte cell biology, physiology, and ecology. Gene family expansions associated with multicellularity are distinct from those of freshwater algae. Candidate genes, including some that arose following horizontal gene transfer from chromalveolates, are present for the transport and metabolism of DMSP. The Ulva genome offers, therefore, new opportunities to understand coastal and marine ecosystems and the fundamental evolution of the green lineage.

Genome Analyses of the Microalga Picochlorum Provide Insights into the Evolution of Thermotolerance in the Green Lineage.

While the molecular events involved in cell responses to heat stress have been extensively studied, our understanding of the genetic basis of basal thermotolerance, and particularly its evolution within the green lineage, remains limited. Here, we present the 13.3-Mb haploid genome and transcriptomes of a halotolerant and thermotolerant unicellular green alga, Picochlorum costavermella (Trebouxiophyceae) to investigate the evolution of the genomic basis of thermotolerance. Differential gene expression at high and standard temperatures revealed that more of the gene families containing up-regulated genes at high temperature were recently evolved, and less originated at the ancestor of green plants. Inversely, there was an excess of ancient gene families containing transcriptionally repressed genes. Interestingly, there is a striking overlap between the thermotolerance and halotolerance transcriptional rewiring, as more than one-third of the gene families up-regulated at 35 °C were also up-regulated under variable salt concentrations in Picochlorum SE3. Moreover, phylogenetic analysis of the 9,304 protein coding genes revealed 26 genes of horizontally transferred origin in P. costavermella, of which five were differentially expressed at higher temperature. Altogether, these results provide new insights about how the genomic basis of adaptation to halo- and thermotolerance evolved in the green lineage.

PLAZA 4.0: an integrative resource for functional, evolutionary and comparative plant genomics.

PLAZA (https://bioinformatics.psb.ugent.be/plaza) is a plant-oriented online resource for comparative, evolutionary and functional genomics. The PLAZA platform consists of multiple independent instances focusing on different plant clades, while also providing access to a consistent set of reference species. Each PLAZA instance contains structural and functional gene annotations, gene family data and phylogenetic trees and detailed gene colinearity information. A user-friendly web interface makes the necessary tools and visualizations accessible, specific for each data type. Here we present PLAZA 4.0, the latest iteration of the PLAZA framework. This version consists of two new instances (Dicots 4.0 and Monocots 4.0) providing a large increase in newly available species, and offers access to updated and newly implemented tools and visualizations, helping users with the ever-increasing demands for complex and in-depth analyzes. The total number of species across both instances nearly doubles from 37 species in PLAZA 3.0 to 71 species in PLAZA 4.0, with a much broader coverage of crop species (e.g. wheat, palm oil) and species of evolutionary interest (e.g. spruce, Marchantia). The new PLAZA instances can also be accessed by a programming interface through a RESTful web service, thus allowing bioinformaticians to optimally leverage the power of the PLAZA platform.

Asgard archaea illuminate the origin of eukaryotic cellular complexity.

The origin and cellular complexity of eukaryotes represent a major enigma in biology. Current data support scenarios in which an archaeal host cell and an alphaproteobacterial (mitochondrial) endosymbiont merged together, resulting in the first eukaryotic cell. The host cell is related to Lokiarchaeota, an archaeal phylum with many eukaryotic features. The emergence of the structural complexity that characterizes eukaryotic cells remains unclear. Here we describe the 'Asgard' superphylum, a group of uncultivated archaea that, as well as Lokiarchaeota, includes Thor-, Odin- and Heimdallarchaeota. Asgard archaea affiliate with eukaryotes in phylogenomic analyses, and their genomes are enriched for proteins formerly considered specific to eukaryotes. Notably, thorarchaeal genomes encode several homologues of eukaryotic membrane-trafficking machinery components, including Sec23/24 and TRAPP domains. Furthermore, we identify thorarchaeal proteins with similar features to eukaryotic coat proteins involved in vesicle biogenesis. Our results expand the known repertoire of 'eukaryote-specific' proteins in Archaea, indicating that the archaeal host cell already contained many key components that govern eukaryotic cellular complexity.