Modeling Site Heterogeneity with Posterior Mean Site Frequency Profiles Accelerates Accurate Phylogenomic Estimation

Gleaning Euglenozoa-specific DNA polymerases in public single-cell transcriptome data

Multiple genes encoding family A DNA polymerases (famA DNAPs), which are evolutionary relatives of DNA polymerase I (PolI) in bacteria and phages, have been found in eukaryotic genomes, and many of these proteins are used mainly in organelles. Among members of the phylum Euglenozoa, distinct types of famA DNAP, PolIA, PolIBCD+, POP, and eugPolA, have been found. It is intriguing how the suite of famA DNAPs had been established during the evolution of Euglenozoa, but the DNAP data have not been sampled from the taxa that sufficiently represent the diversity of this phylum. In particular, little sequence data were available for basal branching species in Euglenozoa until recently. Thanks to the single-cell transcriptome data from symbiontids and phagotrophic euglenids, we have an opportunity to cover the "hole" in the repertory of famA DNAPs in the deep branches in Euglenozoa. The current study identified 16 new famA DNAP sequences in the transcriptome data from 33 phagotrophic euglenids and two symbiontids, respectively. Based on the new famA DNAP sequences, the updated diversity and evolution of famA DNAPs in Euglenozoa are discussed.

Single-cell genomics reveals new rozellid lineages and supports their sister relationship to Microsporidia

The phylum Rozellomycota has been proposed for a group of early-branching holomycotan lineages representing obligate parasites and hyperparasites of zoosporic fungi, oomycotes or phytoplankton. Given their predominantly intracellular lifestyle, rozellids are typically known from environmental ribosomal DNA data, except for the well-studied Rozella species. To date, the phylogenetic relationship between rozellids and microsporidians (Microsporidia) is not fully understood and most reliable hypotheses are based on phylogenomic analyses that incorporate the only publicly available rozellid genome of Rozella allomycis. Here, we provide genomic data of three new rozellid lineages obtained by single-cell sequencing from environmental samples and show with a phylogenomic approach that rozellids form a monophyletic group that is sister to microsporidians, corroborating the previously proposed phylum Rozellomycota. Whereas no mitochondrial genes coding for the respiratory Complex I could be found, we discovered a gene coding for a nucleotide phosphate transporter in one of the three draft genomes. The scattered absence of Complex I genes and scattered presence of nucleotide transporter genes across diverse microsporidian and rozellid lineages suggest that these adaptations to a parasitic lifestyle, which reduce the parasite's capability to synthesize ATP but enables it to steal ATP from its host, evolved independently in microsporidians and rozellids.

A sequence-based evolutionary distance method for Phylogenetic analysis of highly divergent proteins

Because of the limited effectiveness of prevailing phylogenetic methods when applied to highly divergent protein sequences, the phylogenetic analysis problem remains challenging. Here, we propose a sequence-based evolutionary distance algorithm termed sequence distance (SD), which innovatively incorporates site-to-site correlation within protein sequences into the distance estimation. In protein superfamilies, SD can effectively distinguish evolutionary relationships both within and between protein families, producing phylogenetic trees that closely align with those based on structural information, even with sequence identity less than 20%. SD is highly correlated with the similarity of the protein structure, and can calculate evolutionary distances for thousands of protein pairs within seconds using a single CPU, which is significantly faster than most protein structure prediction methods that demand high computational resources and long run times. The development of SD will significantly advance phylogenetics, providing researchers with a more accurate and reliable tool for exploring evolutionary relationships.

DNA Sequences Are as Useful as Protein Sequences for Inferring Deep Phylogenies

Inference of deep phylogenies has almost exclusively used protein rather than DNA sequences based on the perception that protein sequences are less prone to homoplasy and saturation or to issues of compositional heterogeneity than DNA sequences. Here, we analyze a model of codon evolution under an idealized genetic code and demonstrate that those perceptions may be misconceptions. We conduct a simulation study to assess the utility of protein versus DNA sequences for inferring deep phylogenies, with protein-coding data generated under models of heterogeneous substitution processes across sites in the sequence and among lineages on the tree, and then analyzed using nucleotide, amino acid, and codon models. Analysis of DNA sequences under nucleotide-substitution models (possibly with the third codon positions excluded) recovered the correct tree at least as often as analysis of the corresponding protein sequences under modern amino acid models. We also applied the different data-analysis strategies to an empirical dataset to infer the metazoan phylogeny. Our results from both simulated and real data suggest that DNA sequences may be as useful as proteins for inferring deep phylogenies and should not be excluded from such analyses. Analysis of DNA data under nucleotide models has a major computational advantage over protein-data analysis, potentially making it feasible to use advanced models that account for among-site and among-lineage heterogeneity in the nucleotide-substitution process in inference of deep phylogenies.

Resolving Marine-Freshwater Transitions by Diatoms Through a Fog of Gene Tree Discordance

Despite the obstacles facing marine colonists, most lineages of aquatic organisms have colonized and diversified in freshwaters repeatedly. These transitions can trigger rapid morphological or physiological change and, on longer timescales, lead to increased rates of speciation and extinction. Diatoms are a lineage of ancestrally marine microalgae that have diversified throughout freshwater habitats worldwide. We generated a phylogenomic data set of genomes and transcriptomes for 59 diatom taxa to resolve freshwater transitions in one lineage, the Thalassiosirales. Although most parts of the species tree were consistently resolved with strong support, we had difficulties resolving a Paleocene radiation, which affected the placement of one freshwater lineage. This and other parts of the tree were characterized by high levels of gene tree discordance caused by incomplete lineage sorting and low phylogenetic signal. Despite differences in species trees inferred from concatenation versus summary methods and codons versus amino acids, traditional methods of ancestral state reconstruction supported six transitions into freshwaters, two of which led to subsequent species diversification. Evidence from gene trees, protein alignments, and diatom life history together suggest that habitat transitions were largely the product of homoplasy rather than hemiplasy, a condition where transitions occur on branches in gene trees not shared with the species tree. Nevertheless, we identified a set of putatively hemiplasious genes, many of which have been associated with shifts to low salinity, indicating that hemiplasy played a small but potentially important role in freshwater adaptation. Accounting for differences in evolutionary outcomes, in which some taxa became locked into freshwaters while others were able to return to the ocean or become salinity generalists, might help further distinguish different sources of adaptive mutation in freshwater diatoms.

The evolution of endoparasitism and complex life cycles in parasitic platyhelminths

Within flatworms, the vast majority of parasitism is innate to Neodermata, the most derived and diversified group of the phylum Platyhelminthes.1,2 The four major lineages of Neodermata maintain various combinations of life strategies.3 They include both externally (ecto-) and internally feeding (endo-) parasites. Some lineages complete their life cycles directly by infecting a single host, whereas others succeed only through serial infections of multiple hosts of various vertebrate and invertebrate groups. Food sources and modes of digestion add further combinatorial layers to the often incompletely understood mosaic of neodermatan life histories. Their evolutionary trajectories have remained molecularly unresolved because of conflicting evolutionary inferences and a lack of genomic data.4 Here, we generated transcriptomes for nine early branching neodermatan representatives and performed detailed phylogenomic analyses to address these critical gaps. Polyopisthocotylea, mostly hematophagous ectoparasites, form a group with the mostly hematophagous but endoparasitic trematodes (Trematoda), rather than sharing a common ancestor with Monopisthocotylea, ectoparasitic epithelial feeders. Phylogenetic placement of the highly specialized endoparasitic Cestoda alters depending on the model. Regardless of this uncertainty, this study brings an unconventional perspective on the evolution of platyhelminth parasitism, rejecting a common origin for the endoparasitic lifestyle intrinsic to cestodes and trematodes. Instead, our data indicate that complex life cycles and invasion of vertebrates' gut lumen, the hallmark features of these parasites, evolved independently within Neodermata. We propose the demise of the traditionally recognized class Monogenea and the promotion of its two subclasses to the class level as Monopisthocotyla new class and Polyopisthocotyla new class.

Phylogeny of Rhynchium and Its Related Genera (Hymenoptera: Eumeninae) Based on Universal Single-Copy Orthologs and Ultraconserved Elements

The subfamily Eumeninae is a large group of fierce predatory insects that prey mainly on the larvae of Lepidoptera pests. Because of the highly similar morphologies of the genus Rhynchium and its related genera in the subfamily, including Rhynchium Spinola, Allorhynchium van der Vecht, Anterhynchium de Saussure, Pararrhynchium de Saussure, it is essential to delineate their relationships. A previous phylogenetic analysis based on mitochondrial genomes suggested the inconsistent relationships of these genera under traditional classification based on morphological characters. In this study, we first used single-copy orthologs [USCO] and ultraconserved elements [UCE] extracted from 10 newly sequenced low-coverage whole genomes to resolve the phylogenetic relationships of the above genera. The newly sequenced genomes are 152.99 Mb to 211.49 Mb in size with high completeness (BUSCO complete: 91.5-95.6%) and G + C content (36.31-38.76%). Based on extracted 5811 USCOs and 2312 UCEs, the phylogenetic relationships of Rhynchium and its related genera were: ((Allorhynchium + Lissodynerus) + (Pararrhynchium + (Anterhynchium + (Dirhynchium + Rhynchium)))), which was consistent with the mitochondrial genome results. The results supported the genus Rhynchium as monophyletic, whereas Anterhynchium was recovered as paraphyletic, with Anterhynchium (Dirhynchium) as a sister to Rhynchium and hence deserving generic status; In addition, in the genus Pararrhynchium, P. septemfasciatus feanus and P. venkataramani were separated, not clustered on a branch, just as P. septemfasciatus feanus was not together with P. striatum based on mitochondrial genomes. Since Lissodynerus septemfasciatus, the type species of the genus Lissodynerus, was transferred to Pararrhynchium, it is considered that the genus Lissodynerus should be restituted as a valid genus, not a synonym of Pararrhynchium.

Evolutionary conservation of Trichomonas-mycoplasma symbiosis across the host species barrier

Introduction

The protozoan parasite Trichomonas vaginalis is the most common cellular sexually transmitted disease in humans, and the closely related species Trichomonas gallinae is an avian parasite of ecological and economic importance. Phylogenetic evidence suggests T. vaginalis arose during bird to human transmission of a T. gallinae-like ancestor. Trichomonas vaginalis shares a strong clinical association with the independent sexually transmitted pathogen Metamycoplasma (formerly Mycoplasma) hominis, and the uncultured bacterium "Candidatus Malacoplasma (formerly Mycoplasma) girerdii," with the latter association being an order of magnitude stronger. Both bacterial species have been shown to profoundly influence T. vaginalis growth, energy production and virulence-associated mechanisms.

Methods

Evidence for a novel Malacoplasma sp. was discovered by in vivo Illumina metatranscriptomics sequencing of the T. gallinae-infected pigeon mouth. We leveraged published 16S rDNA profiling data from digestive tract of 12 healthy and 24 T. gallinae-infected pigeons to investigate association between the novel Malacoplasma sp. and T. gallinae. We utilised Illumina metagenomics sequencing targeted to pigeon oral and crop samples infected with the novel Malacoplasma sp. to generate its full-length genome sequence. Sequence similarity network analysis was used to compare annotated proteins from the novel Malacoplasma sp. with a range of other related species.

Results

Here we present evidence for a novel Malacoplasma species, related to "Ca. M. girerdii," that is strongly associated with T. gallinae in the upper digestive tract of domestic pigeons. Analysis of the genome sequence revealed gene features apparently specific to a Trichomonas-symbiotic Malacoplasma lineage.

Discussion

These data support a model of long-term association between Trichomonas and Malacoplasma spp. that has been conserved across diversification of the Trichomonas lineage and the host species barrier from birds to human.

One high quality genome and two transcriptome datasets for new species of Mantamonas, a deep-branching eukaryote clade

Mantamonads were long considered to represent an "orphan" lineage in the tree of eukaryotes, likely branching near the most frequently assumed position for the root of eukaryotes. Recent phylogenomic analyses have placed them as part of the "CRuMs" supergroup, along with collodictyonids and rigifilids. This supergroup appears to branch at the base of Amorphea, making it of special importance for understanding the deep evolutionary history of eukaryotes. However, the lack of representative species and complete genomic data associated with them has hampered the investigation of their biology and evolution. Here, we isolated and described two new species of mantamonads, Mantamonas vickermani sp. nov. and Mantamonas sphyraenae sp. nov., for each of which we generated transcriptomic sequence data, as well as a high-quality genome for the latter. The estimated size of the M. sphyraenae genome is 25 Mb; our de novo assembly appears to be highly contiguous and complete with 9,416 predicted protein-coding genes. This near-chromosome-scale genome assembly is the first described for the CRuMs supergroup.

Early Divergence and Gene Exchange Highways in the Evolutionary History of Mesoaciditogales

The placement of a nonhyperthermophilic order Mesoaciditogales as the earliest branching clade within the Thermotogota phylum challenges the prevailing hypothesis that the last common ancestor of Thermotogota was a hyperthermophile. Yet, given the long branch leading to the only two Mesoaciditogales described to date, the phylogenetic position of the order may be due to the long branch attraction artifact. By testing various models and applying data recoding in phylogenetic reconstructions, we observed that early branching of Mesoaciditogales within Thermotogota is strongly supported by the conserved marker genes assumed to be vertically inherited. However, based on the taxonomic content of 1,181 gene families and a phylogenetic analysis of 721 gene family trees, we also found that a substantial number of Mesoaciditogales genes are more closely related to species from the order Petrotogales. These genes contribute to coenzyme transport and metabolism, fatty acid biosynthesis, genes known to respond to heat and cold stressors, and include many genes of unknown functions. The Petrotogales comprise moderately thermophilic and mesophilic species with similar temperature tolerances to that of Mesoaciditogales. Our findings hint at extensive horizontal gene transfer (HGT) between, or parallel independent gene gains by, the two ecologically similar lineages and suggest that the exchanged genes may be important for adaptation to comparable temperature niches.

A predicted CRISPR-mediated symbiosis between uncultivated archaea

CRISPR-Cas systems defend prokaryotic cells from invasive DNA of viruses, plasmids and other mobile genetic elements. Here, we show using metagenomics, metatranscriptomics and single-cell genomics that CRISPR systems of widespread, uncultivated archaea can also target chromosomal DNA of archaeal episymbionts of the DPANN superphylum. Using meta-omics datasets from Crystal Geyser and Horonobe Underground Research Laboratory, we find that CRISPR spacers of the hosts Candidatus Altiarchaeum crystalense and Ca. A. horonobense, respectively, match putative essential genes in their episymbionts' genomes of the genus Ca. Huberiarchaeum and that some of these spacers are expressed in situ. Metabolic interaction modelling also reveals complementation between host-episymbiont systems, on the basis of which we propose that episymbionts are either parasitic or mutualistic depending on the genotype of the host. By expanding our analysis to 7,012 archaeal genomes, we suggest that CRISPR-Cas targeting of genomes associated with symbiotic archaea evolved independently in various archaeal lineages.

Phylogenomics of weevils revisited: data curation and modelling compositional heterogeneity

Weevils represent one of the most prolific radiations of beetles and the most diverse group of herbivores on land. The phylogeny of weevils (Curculionoidea) has received extensive attention, and a largely satisfactory framework for their interfamilial relationships has been established. However, a recent phylogenomic study of Curculionoidea based on anchored hybrid enrichment (AHE) data yielded an abnormal placement for the family Belidae (strongly supported as sister to Nemonychidae + Anthribidae). Here we reanalyse the genome-scale AHE data for Curculionoidea using various models of molecular evolution and data filtering methods to mitigate anticipated systematic errors and reduce compositional heterogeneity. When analysed with the infinite mixture model CAT-GTR or using appropriately filtered datasets, Belidae are always recovered as sister to the clade (Attelabidae, (Caridae, (Brentidae, Curculionidae))), which is congruent with studies based on morphology and other sources of molecular data. Although the relationships of the 'higher Curculionidae' remain challenging to resolve, we provide a consistent and robust backbone phylogeny of weevils. Our extensive analyses emphasize the significance of data curation and modelling across-site compositional heterogeneity in phylogenomic studies.

Molecular phylogenetics of sessile Dolium sedentarium, a petalomonad euglenid

The euglenids are a species-rich group of flagellates with varying modes of nutrition that can be found in diverse habitats. Phagotrophic members of this group gave rise to phototrophs and hold the key to understanding the evolution of euglenids as a whole, including the evolution of complex morphological characters like the euglenid pellicle. Yet to understand the evolution of these characters, a comprehensive sampling of molecular data is needed to correlate morphological and molecular data, and to estimate a basic phylogenetic backbone of the group. While the availability of SSU rDNA and, more recently, multigene data from phagotrophic euglenids has improved, several "orphan" taxa remain without any molecular data whatsoever. Dolium sedentarium is one such taxon: It is a rarely-observed phagotrophic euglenid that inhabits tropical benthic environments and is one of few known sessile euglenids. Based on morphological characters, it has been thought of as part of the earliest branch of euglenids, the Petalomonadida. We report the first molecular sequencing data for Dolium using single-cell transcriptomics, adding another small piece in the puzzle of euglenid evolution. Both SSU rDNA and multigene phylogenies confirm it as a solitary branch within Petalomonadida.

Compositionally Constrained Sites Drive Long-Branch Attraction

Accurate phylogenies are fundamental to our understanding of the pattern and process of evolution. Yet, phylogenies at deep evolutionary timescales, with correspondingly long branches, have been fraught with controversy resulting from conflicting estimates from models with varying complexity and goodness of fit. Analyses of historical as well as current empirical datasets, such as alignments including Microsporidia, Nematoda, or Platyhelminthes, have demonstrated that inadequate modeling of across-site compositional heterogeneity, which is the result of biochemical constraints that lead to varying patterns of accepted amino acids along sequences, can lead to erroneous topologies that are strongly supported. Unfortunately, models that adequately account for across-site compositional heterogeneity remain computationally challenging or intractable for an increasing fraction of contemporary datasets. Here, we introduce "compositional constraint analysis," a method to investigate the effect of site-specific constraints on amino acid composition on phylogenetic inference. We show that more constrained sites with lower diversity and less constrained sites with higher diversity exhibit ostensibly conflicting signals under models ignoring across-site compositional heterogeneity that lead to long-branch attraction artifacts and demonstrate that more complex models accounting for across-site compositional heterogeneity can ameliorate this bias. We present CAT-posterior mean site frequencies (PMSF), a pipeline for diagnosing and resolving phylogenetic bias resulting from inadequate modeling of across-site compositional heterogeneity based on the CAT model. CAT-PMSF is robust against long-branch attraction in all alignments we have examined. We suggest using CAT-PMSF when convergence of the CAT model cannot be assured. We find evidence that compositionally constrained sites are driving long-branch attraction in two metazoan datasets and recover evidence for Porifera as the sister group to all other animals. [Animal phylogeny; cross-site heterogeneity; long-branch attraction; phylogenomics.].

Chromosome-Aware Phylogenomics of Assassin Bugs (Hemiptera: Reduvioidea) Elucidates Ancient Gene Conflict

Though the phylogenetic signal of loci on sex chromosomes can differ from those on autosomes, chromosomal-level genome assemblies for nonvertebrates are still relatively scarce and conservation of chromosomal gene content across deep phylogenetic scales has therefore remained largely unexplored. We here assemble a uniquely large and diverse set of samples (17 anchored hybrid enrichment, 24 RNA-seq, and 70 whole-genome sequencing samples of variable depth) for the medically important assassin bugs (Reduvioidea). We assess the performance of genes based on multiple features (e.g., nucleotide vs. amino acid, nuclear vs. mitochondrial, and autosomal vs. X chromosomal) and employ different methods (concatenation and coalescence analyses) to reconstruct the unresolved phylogeny of this diverse (∼7,000 spp.) and old (>180 Ma) group. Our results show that genes on the X chromosome are more likely to have discordant phylogenies than those on autosomes. We find that the X chromosome conflict is driven by high gene substitution rates that impact the accuracy of phylogenetic inference. However, gene tree clustering showed strong conflict even after discounting variable third codon positions. Alternative topologies were not particularly enriched for sex chromosome loci, but spread across the genome. We conclude that binning genes to autosomal or sex chromosomes may result in a more accurate picture of the complex evolutionary history of a clade.

Major Revisions in Pancrustacean Phylogeny and Evidence of Sensitivity to Taxon Sampling

The clade Pancrustacea, comprising crustaceans and hexapods, is the most diverse group of animals on earth, containing over 80% of animal species and half of animal biomass. It has been the subject of several recent phylogenomic analyses, yet relationships within Pancrustacea show a notable lack of stability. Here, the phylogeny is estimated with expanded taxon sampling, particularly of malacostracans. We show small changes in taxon sampling have large impacts on phylogenetic estimation. By analyzing identical orthologs between two slightly different taxon sets, we show that the differences in the resulting topologies are due primarily to the effects of taxon sampling on the phylogenetic reconstruction method. We compare trees resulting from our phylogenomic analyses with those from the literature to explore the large tree space of pancrustacean phylogenetic hypotheses and find that statistical topology tests reject the previously published trees in favor of the maximum likelihood trees produced here. Our results reject several clades including Caridoida, Eucarida, Multicrustacea, Vericrustacea, and Syncarida. Notably, we find Copepoda nested within Allotriocarida with high support and recover a novel relationship between decapods, euphausiids, and syncarids that we refer to as the Syneucarida. With denser taxon sampling, we find Stomatopoda sister to this latter clade, which we collectively name Stomatocarida, dividing Malacostraca into three clades: Leptostraca, Peracarida, and Stomatocarida. A new Bayesian divergence time estimation is conducted using 13 vetted fossils. We review our results in the context of other pancrustacean phylogenetic hypotheses and highlight 15 key taxa to sample in future studies.

The Fish Pathogen "Candidatus Clavichlamydia salmonicola"-A Missing Link in the Evolution of Chlamydial Pathogens of Humans

Chlamydiae like Chlamydia trachomatis and Chlamydia psittaci are well-known human and animal pathogens. Yet, the chlamydiae are a much larger group of evolutionary ancient obligate intracellular bacteria that includes predominantly symbionts of protists and diverse animals. This makes them ideal model organisms to study evolutionary transitions from symbionts in microbial eukaryotes to pathogens of humans. To this end, comparative genome analysis has served as an important tool. Genome sequence data for many chlamydial lineages are, however, still lacking, hampering our understanding of their evolutionary history. Here, we determined the first high-quality draft genome sequence of the fish pathogen "Candidatus Clavichlamydia salmonicola", representing a separate genus within the human and animal pathogenic Chlamydiaceae. The "Ca. Clavichlamydia salmonicola" genome harbors genes that so far have been exclusively found in Chlamydia species suggesting that basic mechanisms important for the interaction with chordate hosts have evolved stepwise in the history of chlamydiae. Thus, the genome sequence of "Ca. Clavichlamydia salmonicola" allows to constrain candidate genes to further understand the evolution of chlamydial virulence mechanisms required to infect mammals.

Confusion will be my epitaph: genome-scale discordance stifles phylogenetic resolution of Holothuroidea

Sea cucumbers (Holothuroidea) are a diverse clade of echinoderms found from intertidal waters to the bottom of the deepest oceanic trenches. Their reduced skeletons and limited number of phylogenetically informative traits have long obfuscated morphological classifications. Sanger-sequenced molecular datasets have also failed to constrain the position of major lineages. Noteworthy, topological uncertainty has hindered a resolution for Neoholothuriida, a highly diverse clade of Permo-Triassic age. We perform the first phylogenomic analysis of Holothuroidea, combining existing datasets with 13 novel transcriptomes. Using a highly curated dataset of 1100 orthologues, our efforts recapitulate previous results, struggling to resolve interrelationships among neoholothuriid clades. Three approaches to phylogenetic reconstruction (concatenation under both site-homogeneous and site-heterogeneous models, and coalescent-aware inference) result in alternative resolutions, all of which are recovered with strong support and across a range of datasets filtered for phylogenetic usefulness. We explore this intriguing result using gene-wise log-likelihood scores and attempt to correlate these with a large set of gene properties. While presenting novel ways of exploring and visualizing support for alternative trees, we are unable to discover significant predictors of topological preference, and our efforts fail to favour one topology. Neoholothuriid genomes seem to retain an amalgam of signals derived from multiple phylogenetic histories.

Subunit composition of mitochondrial dehydrogenase complexes in diplonemid flagellates

In eukaryotes, pyruvate, a key metabolite produced by glycolysis, is converted by a tripartite mitochondrial pyruvate dehydrogenase (PDH) complex to acetyl-coenzyme A, which is fed into the tricarboxylic acid cycle. Two additional enzyme complexes with analogous composition catalyze similar oxidative decarboxylation reactions albeit using different substrates, the branched-chain ketoacid dehydrogenase (BCKDH) complex and the 2-oxoglutarate dehydrogenase (OGDH) complex. Comparative transcriptome analyses of diplonemids, one of the most abundant and diverse groups of oceanic protists, indicate that the conventional E1, E2, and E3 subunits of the PDH complex are lacking. E1 was apparently replaced in the euglenozoan ancestor of diplonemids by an AceE protein of archaeal type, a substitution that we also document in dinoflagellates. Here we demonstrate that the mitochondrion of the model diplonemid Paradiplonema papillatum displays pyruvate and 2-oxoglutarate dehydrogenase activities. Protein mass spectrometry of mitochondria reveal that the AceE protein is as abundant as the E1 subunit of BCKDH. This corroborates the view that the AceE subunit is a functional component of the PDH complex. We hypothesize that by acquiring AceE, the diplonemid ancestor not only lost the eukaryotic-type E1, but also the E2 and E3 subunits of the PDH complex, which are present in other euglenozoans. We posit that the PDH activity in diplonemids seems to be carried out by a complex, in which the AceE protein partners with the E2 and E3 subunits from BCKDH and/or OGDH.

Phylogenomics of novel ploeotid taxa contribute to the backbone of the euglenid tree

Euglenids are a diverse group of flagellates that inhabit most environments and exhibit many different nutritional modes. The most prominent euglenids are phototrophs, but phagotrophs constitute the majority of phylogenetic diversity of euglenids. They are pivotal to our understanding of euglenid evolution, yet we are only starting to understand relationships amongst phagotrophs, with the backbone of the tree being most elusive. Ploeotids make up most of this backbone diversity-yet despite their morphological similarities, SSU rDNA analyses and multigene analyses show that they are non-monophyletic. As more ploeotid diversity is sampled, known taxa have coalesced into some subgroups (e.g. Alistosa), but the relationships amongst these are not always supported and some taxa remain unsampled for multigene phylogenetics. Here, we used light microscopy and single-cell transcriptomics to characterize five ploeotid euglenids and place them into a multigene phylogenetic framework. Our analyses place Decastava in Alistosa; while Hemiolia branches with Liburna, establishing the novel clade Karavia. We describe Hemiolia limna, a freshwater-dwelling species in an otherwise marine clade. Intriguingly, two undescribed ploeotids are found to occupy pivotal positions in the tree: Chelandium granulatum nov. gen. nov. sp. branches as sister to Olkasia, and Gaulosia striata nov. gen. nov. sp. remains an orphan taxon.

Parallel and convergent genomic changes underlie independent subterranean colonization across beetles

Adaptation to life in caves is often accompanied by dramatically convergent changes across distantly related taxa, epitomized by the loss or reduction of eyes and pigmentation. Nevertheless, the genomic underpinnings underlying cave-related phenotypes are largely unexplored from a macroevolutionary perspective. Here we investigate genome-wide gene evolutionary dynamics in three distantly related beetle tribes with at least six instances of independent colonization of subterranean habitats, inhabiting both aquatic and terrestrial underground systems. Our results indicate that remarkable gene repertoire changes mainly driven by gene family expansions occurred prior to underground colonization in the three tribes, suggesting that genomic exaptation may have facilitated a strict subterranean lifestyle parallelly across beetle lineages. The three tribes experienced both parallel and convergent changes in the evolutionary dynamics of their gene repertoires. These findings pave the way towards a deeper understanding of the evolution of the genomic toolkit in hypogean fauna.

Ancient Rapid Radiation Explains Most Conflicts Among Gene Trees and Well-Supported Phylogenomic Trees of Nostocalean Cyanobacteria

Prokaryotic genomes are often considered to be mosaics of genes that do not necessarily share the same evolutionary history due to widespread horizontal gene transfers (HGTs). Consequently, representing evolutionary relationships of prokaryotes as bifurcating trees has long been controversial. However, studies reporting conflicts among gene trees derived from phylogenomic data sets have shown that these conflicts can be the result of artifacts or evolutionary processes other than HGT, such as incomplete lineage sorting, low phylogenetic signal, and systematic errors due to substitution model misspecification. Here, we present the results of an extensive exploration of phylogenetic conflicts in the cyanobacterial order Nostocales, for which previous studies have inferred strongly supported conflicting relationships when using different concatenated phylogenomic data sets. We found that most of these conflicts are concentrated in deep clusters of short internodes of the Nostocales phylogeny, where the great majority of individual genes have low resolving power. We then inferred phylogenetic networks to detect HGT events while also accounting for incomplete lineage sorting. Our results indicate that most conflicts among gene trees are likely due to incomplete lineage sorting linked to an ancient rapid radiation, rather than to HGTs. Moreover, the short internodes of this radiation fit the expectations of the anomaly zone, i.e., a region of the tree parameter space where a species tree is discordant with its most likely gene tree. We demonstrated that concatenation of different sets of loci can recover up to 17 distinct and well-supported relationships within the putative anomaly zone of Nostocales, corresponding to the observed conflicts among well-supported trees based on concatenated data sets from previous studies. Our findings highlight the important role of rapid radiations as a potential cause of strongly conflicting phylogenetic relationships when using phylogenomic data sets of bacteria. We propose that polytomies may be the most appropriate phylogenetic representation of these rapid radiations that are part of anomaly zones, especially when all possible genomic markers have been considered to infer these phylogenies. [Anomaly zone; bacteria; horizontal gene transfer; incomplete lineage sorting; Nostocales; phylogenomic conflict; rapid radiation; Rhizonema.].

Inference and reconstruction of the heimdallarchaeial ancestry of eukaryotes

In the ongoing debates about eukaryogenesis-the series of evolutionary events leading to the emergence of the eukaryotic cell from prokaryotic ancestors-members of the Asgard archaea play a key part as the closest archaeal relatives of eukaryotes1. However, the nature and phylogenetic identity of the last common ancestor of Asgard archaea and eukaryotes remain unresolved2-4. Here we analyse distinct phylogenetic marker datasets of an expanded genomic sampling of Asgard archaea and evaluate competing evolutionary scenarios using state-of-the-art phylogenomic approaches. We find that eukaryotes are placed, with high confidence, as a well-nested clade within Asgard archaea and as a sister lineage to Hodarchaeales, a newly proposed order within Heimdallarchaeia. Using sophisticated gene tree and species tree reconciliation approaches, we show that analogous to the evolution of eukaryotic genomes, genome evolution in Asgard archaea involved significantly more gene duplication and fewer gene loss events compared with other archaea. Finally, we infer that the last common ancestor of Asgard archaea was probably a thermophilic chemolithotroph and that the lineage from which eukaryotes evolved adapted to mesophilic conditions and acquired the genetic potential to support a heterotrophic lifestyle. Our work provides key insights into the prokaryote-to-eukaryote transition and a platform for better understanding the emergence of cellular complexity in eukaryotic cells.

Evolutionary analysis of cellular reduction and anaerobicity in the hyper-prevalent gut microbe Blastocystis

Blastocystis is the most prevalent microbial eukaryote in the human and animal gut, yet its role as commensal or parasite is still under debate. Blastocystis has clearly undergone evolutionary adaptation to the gut environment and possesses minimal cellular compartmentalization, reduced anaerobic mitochondria, no flagella, and no reported peroxisomes. To address this poorly understood evolutionary transition, we have taken a multi-disciplinary approach to characterize Proteromonas lacertae, the closest canonical stramenopile relative of Blastocystis. Genomic data reveal an abundance of unique genes in P. lacertae but also reductive evolution of the genomic complement in Blastocystis. Comparative genomic analysis sheds light on flagellar evolution, including 37 new candidate components implicated with mastigonemes, the stramenopile morphological hallmark. The P. lacertae membrane-trafficking system (MTS) complement is only slightly more canonical than that of Blastocystis, but notably, we identified that both organisms encode the complete enigmatic endocytic TSET complex, a first for the entire stramenopile lineage. Investigation also details the modulation of mitochondrial composition and metabolism in both P. lacertae and Blastocystis. Unexpectedly, we identify in P. lacertae the most reduced peroxisome-derived organelle reported to date, which leads us to speculate on a mechanism of constraint guiding the dynamics of peroxisome-mitochondrion reductive evolution on the path to anaerobiosis. Overall, these analyses provide a launching point to investigate organellar evolution and reveal in detail the evolutionary path that Blastocystis has taken from a canonical flagellated protist to the hyper-divergent and hyper-prevalent animal and human gut microbe.

Mirusviruses link herpesviruses to giant viruses

DNA viruses have a major influence on the ecology and evolution of cellular organisms1-4, but their overall diversity and evolutionary trajectories remain elusive5. Here we carried out a phylogeny-guided genome-resolved metagenomic survey of the sunlit oceans and discovered plankton-infecting relatives of herpesviruses that form a putative new phylum dubbed Mirusviricota. The virion morphogenesis module of this large monophyletic clade is typical of viruses from the realm Duplodnaviria6, with multiple components strongly indicating a common ancestry with animal-infecting Herpesvirales. Yet, a substantial fraction of mirusvirus genes, including hallmark transcription machinery genes missing in herpesviruses, are closely related homologues of giant eukaryotic DNA viruses from another viral realm, Varidnaviria. These remarkable chimaeric attributes connecting Mirusviricota to herpesviruses and giant eukaryotic viruses are supported by more than 100 environmental mirusvirus genomes, including a near-complete contiguous genome of 432 kilobases. Moreover, mirusviruses are among the most abundant and active eukaryotic viruses characterized in the sunlit oceans, encoding a diverse array of functions used during the infection of microbial eukaryotes from pole to pole. The prevalence, functional activity, diversification and atypical chimaeric attributes of mirusviruses point to a lasting role of Mirusviricota in the ecology of marine ecosystems and in the evolution of eukaryotic DNA viruses.

New Mitogenomes of the Polypedilum Generic Complex (Diptera: Chironomidae): Characterization and Phylogenetic Implications

Mitochondrial genomics, as a useful marker for phylogenetics and systematics of organisms, are important for molecular biology studies. The phylogenetic relationships of the Polypedilum generic complex remains controversial, due to lack taxonomy and molecular information. In this study, we newly sequenced mitogenomes of 14 species of the Polypedilum generic complex. Coupled with three recently published sequences, we analyzed the nucleotide composition, sequence length, and evolutionary rate of this generic complex. The control region showed the highest AT content. The evolution rate of protein coding genes was as follows: ATP8 > ND6 > ND5 > ND3 > ND2 > ND4L > ND4 > COX1 > ND1 > CYTB > APT6 > COX2 > COX3. We reconstructed the phylogenetic relationships among the genera within the Polypedilum generic complex based on 19 mitochondrial genomes (seventeen ingroups and two outgroups), using Bayesian Inference (BI) and Maximum Likelihood (ML) methods for all databases. Phylogenetic analysis of 19 mitochondrial genomes demonstrated that the Endochironomus + Synendotendipes was sister to Phaenopsectra + Sergentia.

Phylogenetic placement of the protosteloid amoeba Microglomus paxillus identifies another case of sporocarpic fruiting in Discosea (Amoebozoa)

Protosteloid amoebae are a paraphyletic assemblage of amoeboid protists found exclusively in the eukaryotic assemblage Amoebozoa. These amoebae can facultatively form a dispersal structure known as a fruiting body, or more specifically, a sporocarp, from a single amoeboid cell. Sporocarps consist of one to a few spores atop a noncellular stalk. Protosteloid amoebae are known in two out of three well-established major assemblages of Amoebozoa. Amoebae with a protosteloid life cycle are known in the major Amoebozoa lineages Discosea and Evosea but not in Tubulinea. To date, only one genus, which is monotypic, lacks sequence data and, therefore, remains phylogenetically homeless. To further clarify the evolutionary milieu of sporocarpic fruiting we used single-cell transcriptomics to obtain data from individual sporocarps of isolates of the protosteloid amoeba Microglomus paxillus. Our phylogenomic analyses using 229 protein coding markers suggest that M. paxillus is a member of the Discosea lineage of Amoebozoa most closely related to Mycamoeba gemmipara. Due to the hypervariable nature of the SSU rRNA sequence we were unable to further resolve the phylogenetic position of M. paxillus in taxon rich datasets using only this marker. Regardless, our results widen the known distribution of sporocarpy in Discosea and stimulate the debate between a single or multiple origins of sporocarpic fruiting in Amoebozoa.

The Structure of Evolutionary Model Space for Proteins across the Tree of Life

The factors that determine the relative rates of amino acid substitution during protein evolution are complex and known to vary among taxa. We estimated relative exchangeabilities for pairs of amino acids from clades spread across the tree of life and assessed the historical signal in the distances among these clade-specific models. We separately trained these models on collections of arbitrarily selected protein alignments and on ribosomal protein alignments. In both cases, we found a clear separation between the models trained using multiple sequence alignments from bacterial clades and the models trained on archaeal and eukaryotic data. We assessed the predictive power of our novel clade-specific models of sequence evolution by asking whether fit to the models could be used to identify the source of multiple sequence alignments. Model fit was generally able to correctly classify protein alignments at the level of domain (bacterial versus archaeal), but the accuracy of classification at finer scales was much lower. The only exceptions to this were the relatively high classification accuracy for two archaeal lineages: Halobacteriaceae and Thermoprotei. Genomic GC content had a modest impact on relative exchangeabilities despite having a large impact on amino acid frequencies. Relative exchangeabilities involving aromatic residues exhibited the largest differences among models. There were a small number of exchangeabilities that exhibited large differences in comparisons among major clades and between generalized models and ribosomal protein models. Taken as a whole, these results reveal that a small number of relative exchangeabilities are responsible for much of the structure of the "model space" for protein sequence evolution. The clade-specific models we generated may be useful tools for protein phylogenetics, and the structure of evolutionary model space that they revealed has implications for phylogenomic inference across the tree of life.

Characterization of the Complete Mitochondrial Genome of a Flea Beetle Luperomorpha xanthodera (Coleoptera: Chrysomelidae: Galerucinae) and Phylogenetic Analysis

In this study, the mitochondrial genome of Luperomorpha xanthodera was assembled and annotated, which is a circular DNA molecule including 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, 2 ribosomal RNA genes (12S rRNA and 16S rRNA), and 1388 bp non-coding regions (A + T rich region), measuring 16,021 bp in length. The nucleotide composition of the mitochondrial genome is 41.3% adenine (A), 38.7% thymine (T), 8.4% guanine (G), and 11.6% cytosine (C). Most of the protein-coding genes presented a typical ATN start codon (ATA, ATT, ATC, ATG), except for ND1, which showed the start codon TTG. Three-quarters of the protein-coding genes showed the complete stop codon TAR (TAA, TAG), except the genes COI, COII, ND4, and ND5, which showed incomplete stop codons (T- or TA-). All the tRNA genes have the typical clover-leaf structure, except tRNASer1 (AGN), which has a missing dihydrouridine arm (DHU). The phylogenetic results determined by both maximum likelihood and Bayesian inference methods consistently supported the monophyly of the subfamily Galerucinae and revealed that the subtribe Luperina and genus Monolepta are polyphyletic groups. Meanwhile, the classification status of the genus Luperomorpha is controversial.

Chromosome-level genomes of multicellular algal sisters to land plants illuminate signaling network evolution

The filamentous and unicellular algae of the class Zygnematophyceae are the closest algal relatives of land plants. Inferring the properties of the last common ancestor shared by these algae and land plants allows us to identify decisive traits that enabled the conquest of land by plants. We sequenced four genomes of filamentous Zygnematophyceae (three strains of Zygnema circumcarinatum and one strain of Z. cylindricum) and generated chromosome-scale assemblies for all strains of the emerging model system Z. circumcarinatum. Comparative genomic analyses reveal expanded genes for signaling cascades, environmental response, and intracellular trafficking that we associate with multicellularity. Gene family analyses suggest that Zygnematophyceae share all the major enzymes with land plants for cell wall polysaccharide synthesis, degradation, and modifications; most of the enzymes for cell wall innovations, especially for polysaccharide backbone synthesis, were gained more than 700 million years ago. In Zygnematophyceae, these enzyme families expanded, forming co-expressed modules. Transcriptomic profiling of over 19 growth conditions combined with co-expression network analyses uncover cohorts of genes that unite environmental signaling with multicellular developmental programs. Our data shed light on a molecular chassis that balances environmental response and growth modulation across more than 600 million years of streptophyte evolution.

Micro and macroevolution of sea anemone venom phenotype

Venom is a complex trait with substantial inter- and intraspecific variability resulting from strong selective pressures acting on the expression of many toxic proteins. However, understanding the processes underlying toxin expression dynamics that determine the venom phenotype remains unresolved. By interspecific comparisons we reveal that toxin expression in sea anemones evolves rapidly and that in each species different toxin family dictates the venom phenotype by massive gene duplication events. In-depth analysis of the sea anemone, Nematostella vectensis, revealed striking variation of the dominant toxin (Nv1) diploid copy number across populations (1-24 copies) resulting from independent expansion/contraction events, which generate distinct haplotypes. Nv1 copy number correlates with expression at both the transcript and protein levels with one population having a near-complete loss of Nv1 production. Finally, we establish the dominant toxin hypothesis which incorporates observations in other venomous lineages that animals have convergently evolved a similar strategy in shaping their venom.

Building Phylogenies from Transcriptomic Data

Transcriptomic data (obtained from RNA sequencing) has become a very powerful source of information to reconstruct the evolutionary relationships among organisms. Although phylogenetic inference using transcriptomes retains the same core steps as when working with few molecular markers (viz., nucleic acid extraction and sequencing, sequence treatment, and tree inference), all of them show significant differences. First, the needed quantity and quality of the extracted RNA has to be very high. Although this may not represent a challenge when working with certain organisms, it may well be a headache with others, especially for those with small body sizes. Second, the tremendous increase in the quantity of sequences obtained requires a high computational power for both treating the sequences and inferring the subsequent phylogenies. This means that transcriptomic data can no longer be analyzed using personal computers nor local programs with a graphical interface. This, in turn, implies the requirement of an increased set of bioinformatic skills from the researchers. Finally, the genomic peculiarities of each group of organisms, such as the level of heterozygosity or the percentage of base composition, also need to be considered when inferring phylogenies using transcriptomic data.

Phylotranscriptomics interrogation uncovers a complex evolutionary history for the planarian genus Dugesia (Platyhelminthes, Tricladida) in the Western Mediterranean

The Mediterranean is one of the most biodiverse areas of the Paleartic region. Here, basing on large data sets of single copy orthologs obtained from transcriptomic data, we investigated the evolutionary history of the genus Dugesia in the Western Mediterranean area. The results corroborated that the complex paleogeological history of the region was an important driver of diversification for the genus, speciating as microplates and islands were forming. These processes led to the differentiation of three main biogeographic clades: Iberia-Apennines-Alps, Corsica-Sardinia, and Iberia-Africa. The internal relationships of these major clades were analysed with several representative samples per species. The use of large data sets regarding the number of loci and samples, as well as state-of-the-art phylogenomic inference methods allowed us to answer different unresolved questions about the evolution of particular groups, such as the diversification path of D. subtentaculata in the Iberian Peninsula and its colonization of Africa. Additionally, our results support the differentiation of D. benazzii in two lineages which could represent two species. Finally, we analysed here for the first time a comprehensive number of samples from several asexual Iberian populations whose assignment at the species level has been an enigma through the years. The phylogenies obtained with different inference methods showed a branching topology of asexual individuals at the base of sexual clades. We hypothesize that this unexpected topology is related to long-term asexuality. This work represents the first phylotranscriptomic analysis of Tricladida, laying the first stone of the genomic era in phylogenetic studies on this taxonomic group.

Reduced mitochondria provide an essential function for the cytosolic methionine cycle

The loss of mitochondria in oxymonad protists has been associated with the redirection of the essential Fe-S cluster assembly to the cytosol. Yet as our knowledge of diverse free-living protists broadens, the list of functions of their mitochondrial-related organelles (MROs) expands. We revealed another such function in the closest oxymonad relative, Paratrimastix pyriformis, after we solved the proteome of its MRO with high accuracy, using localization of organelle proteins by isotope tagging (LOPIT). The newly assigned enzymes connect to the glycine cleavage system (GCS) and produce folate derivatives with one-carbon units and formate. These are likely to be used by the cytosolic methionine cycle involved in S-adenosyl methionine recycling. The data provide consistency with the presence of the GCS in MROs of free-living species and its absence in most endobionts, which typically lose the methionine cycle and, in the case of oxymonads, the mitochondria.

Testing the systematic status of Homalictus and Rostrohalictus with weakened cross-vein groups within Halictini (Hymenoptera: Halictidae) using low-coverage whole-genome sequencing

The halictid genus Lasioglossum, as one of the most species-rich bee groups with persistently contentious subgeneric boundaries, is one of the most challenging bee groups from a systematic standpoint. An enduring question is the relationship of Lasioglossum and Homalictus, whether all halictine bees with weakened distal wing venation comprise one or multiple genera. Here, we analyzed the phylogenetic relationships among the subgroups within Lasioglossum s.l. based on thousands of single-copy orthologs and ultraconserved elements, which were extracted from 23 newly sequenced low-coverage whole genomes alongside a published genome (22 ingroups plus 2 outgroups). Both marker sets provided consistent results across maximum likelihood and coalescent-based species tree approaches. The phylogenetic and topology test results show that the Lasioglossum and Hemihalictus series are reciprocally monophyletic and Homalictus and Rostrohalictus are valid subgenera of Lasioglossum. Consequently, we lower Homalictus to subgenus status within Lasioglossum again, and we also raise Rostrohalictus to subgenus status from its prior synonymy with subgenus Hemihalictus. Lasioglossum przewalskyi is also transferred to the subgenus Hemihalictus. Ultimately, we redefine Lasioglossum to include all halictine bees with weakened distal wing venation.

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

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

Genomic diversity and biosynthetic capabilities of sponge-associated chlamydiae

Sponge microbiomes contribute to host health, nutrition, and defense through the production of secondary metabolites. Chlamydiae, a phylum of obligate intracellular bacteria ranging from animal pathogens to endosymbionts of microbial eukaryotes, are frequently found associated with sponges. However, sponge-associated chlamydial diversity has not yet been investigated at the genomic level and host interactions thus far remain unexplored. Here, we sequenced the microbiomes of three sponge species and found high, though variable, Chlamydiae relative abundances of up to 18.7% of bacteria. Using genome-resolved metagenomics 18 high-quality sponge-associated chlamydial genomes were reconstructed, covering four chlamydial families. Among these, Candidatus Sororchlamydiaceae shares a common ancestor with Chlamydiaceae animal pathogens, suggesting long-term co-evolution with animals. Based on gene content, sponge-associated chlamydiae resemble members from the same family more than sponge-associated chlamydiae of other families, and have greater metabolic versatility than known chlamydial animal pathogens. Sponge-associated chlamydiae are also enriched in genes for degrading diverse compounds found in sponges. Unexpectedly, we identified widespread genetic potential for secondary metabolite biosynthesis across Chlamydiae, which may represent an unexplored source of novel natural products. This finding suggests that Chlamydiae members may partake in defensive symbioses and that secondary metabolites play a wider role in mediating intracellular interactions. Furthermore, sponge-associated chlamydiae relatives were found in other marine invertebrates, pointing towards wider impacts of the Chlamydiae phylum on marine ecosystems.

Microbial predators form a new supergroup of eukaryotes

Molecular phylogenetics of microbial eukaryotes has reshaped the tree of life by establishing broad taxonomic divisions, termed supergroups, that supersede the traditional kingdoms of animals, fungi and plants, and encompass a much greater breadth of eukaryotic diversity¹. The vast majority of newly discovered species fall into a small number of known supergroups. Recently, however, a handful of species with no clear relationship to other supergroups have been described^2-4, raising questions about the nature and degree of undiscovered diversity, and exposing the limitations of strictly molecular-based exploration. Here we report ten previously undescribed strains of microbial predators isolated through culture that collectively form a diverse new supergroup of eukaryotes, termed Provora. The Provora supergroup is genetically, morphologically and behaviourally distinct from other eukaryotes, and comprises two divergent clades of predators-Nebulidia and Nibbleridia-that are superficially similar to each other, but differ fundamentally in ultrastructure, behaviour and gene content. These predators are globally distributed in marine and freshwater environments, but are numerically rare and have consequently been overlooked by molecular-diversity surveys. In the age of high-throughput analyses, investigation of eukaryotic diversity through culture remains indispensable for the discovery of rare but ecologically and evolutionarily important eukaryotes.

New Mitogenome Features of Philopotamidae (Insecta: Trichoptera) with Two New Species of Gunungiella

A total of 14 individuals of Philopotamidae, from China, were examined. Six species in four genera, including two new species of the genus Gunungiella, were recognized. Their COI barcode sequences were extracted, mitogenomes were sequenced, assembled and analyzed. All of these sequences were used to further reveal the phylogenetic relationships of the family Philopotamidae. In addition, two new species: Gunungiella wangi n. sp., Gunungiella flabellata n. sp. were described and illustrated.

Phylogenomics reveals deep relationships and diversification within phylactolaemate bryozoans

Bryozoans are mostly sessile colonial invertebrates that inhabit all kinds of aquatic ecosystems. Extant bryozoan species fall into two clades with one of them, Phylactolaemata, being the only exclusively freshwater clade. Phylogenetic relationships within the class Phylactolaemata have long been controversial owing to their limited distinguishable characteristics that reflect evolutionary relationships. Here, we present the first phylogenomic analysis of Phylactolaemata using transcriptomic data combined with dense taxon sampling of six families to better resolve the interrelationships and to estimate divergence time. Using maximum-likelihood and Bayesian inference approaches, we recovered a robust phylogeny for Phylactolaemata in which the interfamilial relationships are fully resolved. We show Stephanellidae is the sister taxon of all other phylactolaemates and confirm that Lophopodidae represents the second offshoot within the phylactolaemate tree. Plumatella fruticosa clearly falls outside Plumatellidae as previous investigations have suggested, and instead clusters with Pectinatellidae and Cristatellidae as the sister taxon of Fredericellidae. Our results demonstrate that cryptic speciation is very likely in F. sultana and in two species of Plumatella (P. repens and P. casmiana). Divergence time estimates show that Phylactolaemata appeared at the end of the Ediacaran and started to diverge in the Silurian, although confidence intervals were large for most nodes. The radiation of most extant phylactolaemate families occurred mainly in the Palaeogene and Neogene highlighting post-extinction diversification.

Genomic remnants of ancestral methanogenesis and hydrogenotrophy in Archaea drive anaerobic carbon cycling

Anaerobic methane metabolism is among the hallmarks of Archaea, originating very early in their evolution. Here, we show that the ancestor of methane metabolizers was an autotrophic CO₂-reducing hydrogenotrophic methanogen that possessed the two main complexes, methyl-CoM reductase (Mcr) and tetrahydromethanopterin-CoM methyltransferase (Mtr), the anaplerotic hydrogenases Eha and Ehb, and a set of other genes collectively called "methanogenesis markers" but could not oxidize alkanes. Overturning recent inferences, we demonstrate that methyl-dependent hydrogenotrophic methanogenesis has emerged multiple times independently, either due to a loss of Mtr while Mcr is inherited vertically or from an ancient lateral acquisition of Mcr. Even if Mcr is lost, Mtr, Eha, Ehb, and the markers can persist, resulting in mixotrophic metabolisms centered around the Wood-Ljungdahl pathway. Through their methanogenesis remnants, Thorarchaeia and two newly reconstructed order-level lineages in Archaeoglobi and Bathyarchaeia act as metabolically versatile players in carbon cycling of anoxic environments across the globe.

Divergent evolutionary trajectories of bryophytes and tracheophytes from a complex common ancestor of land plants

The origin of plants and their colonization of land fundamentally transformed the terrestrial environment. Here we elucidate the basis of this formative episode in Earth history through patterns of lineage, gene and genome evolution. We use new fossil calibrations, a relative clade age calibration (informed by horizontal gene transfer) and new phylogenomic methods for mapping gene family origins. Distinct rooting strategies resolve tracheophytes (vascular plants) and bryophytes (non-vascular plants) as monophyletic sister groups that diverged during the Cambrian, 515-494 million years ago. The embryophyte stem is characterized by a burst of gene innovation, while bryophytes subsequently experienced an equally dramatic episode of reductive genome evolution in which they lost genes associated with the elaboration of vasculature and the stomatal complex. Overall, our analyses reveal that extant tracheophytes and bryophytes are both highly derived from a more complex ancestral land plant. Understanding the origin of land plants requires tracing character evolution across a diversity of modern lineages.

A phylogenomically informed five-order system for the closest relatives of land plants

The evolution of streptophytes had a profound impact on life on Earth. They brought forth those photosynthetic eukaryotes that today dominate the macroscopic flora: the land plants (Embryophyta).¹ There is convincing evidence that the unicellular/filamentous Zygnematophyceae-and not the morphologically more elaborate Coleochaetophyceae or Charophyceae-are the closest algal relatives of land plants.^2-6 Despite the species richness (>4,000), wide distribution, and key evolutionary position of the zygnematophytes, their internal phylogeny remains largely unresolved.^7,8 There are also putative zygnematophytes with interesting body plan modifications (e.g., filamentous growth) whose phylogenetic affiliations remain unknown. Here, we studied a filamentous green alga (strain MZCH580) from an Austrian peat bog with central or parietal chloroplasts that lack discernible pyrenoids. It represents Mougeotiopsis calospora PALLA, an enigmatic alga that was described more than 120 years ago⁹ but never subjected to molecular analyses. We generated transcriptomic data of M. calospora strain MZCH580 and conducted comprehensive phylogenomic analyses (326 nuclear loci) for 46 taxonomically diverse zygnematophytes. Strain MZCH580 falls in a deep-branching zygnematophycean clade together with some unicellular species and thus represents a formerly unknown zygnematophycean lineage with filamentous growth. Our well-supported phylogenomic tree lets us propose a new five-order system for the Zygnematophyceae and provides evidence for at least five independent origins of true filamentous growth in the closest algal relatives of land plants. This phylogeny provides a robust and comprehensive framework for performing comparative analyses and inferring the evolution of cellular traits and body plans in the closest relatives of land plants.

An early origin of iron-sulfur cluster biosynthesis machineries before Earth oxygenation

Iron-sulfur (Fe-S) clusters are ubiquitous cofactors essential for life. It is largely thought that the emergence of oxygenic photosynthesis and progressive oxygenation of the atmosphere led to the origin of multiprotein machineries (ISC, NIF and SUF) assisting Fe-S cluster synthesis in the presence of oxidative stress and shortage of bioavailable iron. However, previous analyses have left unclear the origin and evolution of these systems. Here, we combine exhaustive homology searches with genomic context analysis and phylogeny to precisely identify Fe-S cluster biogenesis systems in over 10,000 archaeal and bacterial genomes. We highlight the existence of two additional and clearly distinct 'minimal' Fe-S cluster assembly machineries, MIS (minimal iron-sulfur) and SMS (SUF-like minimal system), which we infer in the last universal common ancestor (LUCA) and we experimentally validate SMS as a bona fide Fe-S cluster biogenesis system. These ancestral systems were kept in archaea whereas they went through stepwise complexification in bacteria to incorporate additional functions for higher Fe-S cluster synthesis efficiency leading to SUF, ISC and NIF. Horizontal gene transfers and losses then shaped the current distribution of these systems, driving ecological adaptations such as the emergence of aerobic lifestyles in archaea. Our results show that dedicated machineries were in place early in evolution to assist Fe-S cluster biogenesis and that their origin is not directly linked to Earth oxygenation.

A Practical Guide to Design and Assess a Phylogenomic Study

Over the last decade, molecular systematics has undergone a change of paradigm as high-throughput sequencing now makes it possible to reconstruct evolutionary relationships using genome-scale datasets. The advent of "big data" molecular phylogenetics provided a battery of new tools for biologists but simultaneously brought new methodological challenges. The increase in analytical complexity comes at the price of highly specific training in computational biology and molecular phylogenetics, resulting very often in a polarized accumulation of knowledge (technical on one side and biological on the other). Interpreting the robustness of genome-scale phylogenetic studies is not straightforward, particularly as new methodological developments have consistently shown that the general belief of "more genes, more robustness" often does not apply, and because there is a range of systematic errors that plague phylogenomic investigations. This is particularly problematic because phylogenomic studies are highly heterogeneous in their methodology, and best practices are often not clearly defined. The main aim of this article is to present what I consider as the ten most important points to take into consideration when planning a well-thought-out phylogenomic study and while evaluating the quality of published papers. The goal is to provide a practical step-by-step guide that can be easily followed by nonexperts and phylogenomic novices in order to assess the technical robustness of phylogenomic studies or improve the experimental design of a project.

Divergent genomic trajectories predate the origin of animals and fungi

Animals and fungi have radically distinct morphologies, yet both evolved within the same eukaryotic supergroup: Opisthokonta^1,2. Here we reconstructed the trajectory of genetic changes that accompanied the origin of Metazoa and Fungi since the divergence of Opisthokonta with a dataset that includes four novel genomes from crucial positions in the Opisthokonta phylogeny. We show that animals arose only after the accumulation of genes functionally important for their multicellularity, a tendency that began in the pre-metazoan ancestors and later accelerated in the metazoan root. By contrast, the pre-fungal ancestors experienced net losses of most functional categories, including those gained in the path to Metazoa. On a broad-scale functional level, fungal genomes contain a higher proportion of metabolic genes and diverged less from the last common ancestor of Opisthokonta than did the gene repertoires of Metazoa. Metazoa and Fungi also show differences regarding gene gain mechanisms. Gene fusions are more prevalent in Metazoa, whereas a larger fraction of gene gains were detected as horizontal gene transfers in Fungi and protists, in agreement with the long-standing idea that transfers would be less relevant in Metazoa due to germline isolation^3-5. Together, our results indicate that animals and fungi evolved under two contrasting trajectories of genetic change that predated the origin of both groups. The gradual establishment of two clearly differentiated genomic contexts thus set the stage for the emergence of Metazoa and Fungi.

Comparative Mitogenomic Analyses of Hydropsychidae Revealing the Novel Rearrangement of Protein-Coding Gene and tRNA (Trichoptera: Annulipalpia)

Gene rearrangement of the mitochondrial genome of insects, especially the rearrangement of protein-coding genes, has long been a hot topic for entomologists. Although mitochondrial gene rearrangement is common within Annulipalpia, protein-coding gene rearrangement is relatively rare. As the largest family in Annulipalpia, the available mitogenomes from Hydropsychidae Curtis, 1835 are scarce, and thus restrict our interpretation of the mitogenome characteristic. In this study, we obtained 19 novel mitogenomes of Hydropsychidae, of which the mitogenomes of the genus Arctopsyche are published for the first time. Coupled with published hydropsychid mitogenome, we analyzed the nucleotide composition evolutionary rates and gene rearrangements of the mitogenomes among subfamilies. As a result, we found two novel gene rearrangement patterns within Hydropsychidae, including rearrangement of protein-coding genes. Meanwhile, our results consider that the protein-coding gene arrangement of Potamyia can be interpreted by the tandem duplication/random loss (TDRL) model. In addition, the phylogenetic relationships within Hydropsychidae constructed by two strategies (Bayesian inference and maximum likelihood) strongly support the monophyly of Arctopscychinae, Diplectroninae, Hydropsychinae, and Macronematinae. Our study provides new insights into the mechanisms and patterns of mitogenome rearrangements in Hydropsychidae.

A novel structural maintenance of chromosomes (SMC)-related protein family specific to Archaea

The ATPases belonging to the structural maintenance of chromosomes (SMC) superfamily are involved in the maintenance of chromosome organization and dynamics, as well as DNA repair. The major proteins in this superfamily recognized to date are either conserved among the three domains of Life (i.e., SMC and Rad50) or specific to Bacteria (i.e., RecF, RecN, and MukB). In Archaea, no protein related to SMC (SMC-related protein) with a broad taxonomic distribution has been reported. Nevertheless, two SMC-related proteins, namely coalescin and Sph, have been identified in crenarchaea Sulfolobus spp. and the euryarchaeon Halobacterium salinarum, respectively, hinting that the diversity of SMC-related proteins has been overlooked in Archaea. In this study, we report a novel SMC-related protein that is distributed among broad archaeal lineages and termed “Archaea-specific SMC-related proteins” or “ASRPs.” We further demonstrate that the ASRP family encloses both coalescin and Sph but the two proteins represent only a tip of the diversity of this family.