Structure of the endosomal Commander complex linked to Ritscher-Schinzel syndrome.

The Commander complex is required for endosomal recycling of diverse transmembrane cargos and is mutated in Ritscher-Schinzel syndrome. It comprises two sub-assemblies: Retriever composed of VPS35L, VPS26C, and VPS29; and the CCC complex which contains twelve subunits: COMMD1-COMMD10 and the coiled-coil domain-containing (CCDC) proteins CCDC22 and CCDC93. Combining X-ray crystallography, electron cryomicroscopy, and in silico predictions, we have assembled a complete structural model of Commander. Retriever is distantly related to the endosomal Retromer complex but has unique features preventing the shared VPS29 subunit from interacting with Retromer-associated factors. The COMMD proteins form a distinctive hetero-decameric ring stabilized by extensive interactions with CCDC22 and CCDC93. These adopt a coiled-coil structure that connects the CCC and Retriever assemblies and recruits a 16th subunit, DENND10, to form the complete Commander complex. The structure allows mapping of disease-causing mutations and reveals the molecular features required for the function of this evolutionarily conserved trafficking machinery.

Bacterial Vipp1 and PspA are members of the ancient ESCRT-III membrane-remodeling superfamily.

Membrane remodeling and repair are essential for all cells. Proteins that perform these functions include Vipp1/IM30 in photosynthetic plastids, PspA in bacteria, and ESCRT-III in eukaryotes. Here, using a combination of evolutionary and structural analyses, we show that these protein families are homologous and share a common ancient evolutionary origin that likely predates the last universal common ancestor. This homology is evident in cryo-electron microscopy structures of Vipp1 rings from the cyanobacterium Nostoc punctiforme presented over a range of symmetries. Each ring is assembled from rungs that stack and progressively tilt to form dome-shaped curvature. Assembly is facilitated by hinges in the Vipp1 monomer, similar to those in ESCRT-III proteins, which allow the formation of flexible polymers. Rings have an inner lumen that is able to bind and deform membranes. Collectively, these data suggest conserved mechanistic principles that underlie Vipp1, PspA, and ESCRT-III-dependent membrane remodeling across all domains of life.

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: Opisthokonta1,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 isolation3-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.

AleRax: a tool for gene and species tree co-estimation and reconciliation under a probabilistic model of gene duplication, transfer, and loss.

MOTIVATION: Genomes are a rich source of information on the pattern and process of evolution across biological scales. How best to make use of that information is an active area of research in phylogenetics. Ideally, phylogenetic methods should not only model substitutions along gene trees, which explain differences between homologous gene sequences, but also the processes that generate the gene trees themselves along a shared species tree. To conduct accurate inferences, one needs to account for uncertainty at both levels, that is, in gene trees estimated from inherently short sequences and in their diverse evolutionary histories along a shared species tree. RESULTS: We present AleRax, a software that can infer reconciled gene trees together with a shared species tree using a simple, yet powerful, probabilistic model of gene duplication, transfer, and loss. A key feature of AleRax is its ability to account for uncertainty in the gene tree and its reconciliation by using an efficient approximation to calculate the joint phylogenetic-reconciliation likelihood and sample reconciled gene trees accordingly. Simulations and analyses of empirical data show that AleRax is one order of magnitude faster than competing gene tree inference tools while attaining the same accuracy. It is consistently more robust than species tree inference methods such as SpeciesRax and ASTRAL-Pro 2 under gene tree uncertainty. Finally, AleRax can process multiple gene families in parallel thereby allowing users to compare competing phylogenetic hypotheses and estimate model parameters, such as duplication, transfer, and loss probabilities for genome-scale datasets with hundreds of taxa. AVAILABILITY AND IMPLEMENTATION: GNU GPL at https://github.com/BenoitMorel/AleRax and data are made available at https://cme.h-its.org/exelixis/material/alerax_data.tar.gz.

SNX27-Retromer directly binds ESCPE-1 to transfer cargo proteins during endosomal recycling.

Coat complexes coordinate cargo recognition through cargo adaptors with biogenesis of transport carriers during integral membrane protein trafficking. Here, we combine biochemical, structural, and cellular analyses to establish the mechanistic basis through which SNX27-Retromer, a major endosomal cargo adaptor, couples to the membrane remodeling endosomal SNX-BAR sorting complex for promoting exit 1 (ESCPE-1). In showing that the SNX27 FERM (4.1/ezrin/radixin/moesin) domain directly binds acidic-Asp-Leu-Phe (aDLF) motifs in the SNX1/SNX2 subunits of ESCPE-1, we propose a handover model where SNX27-Retromer captured cargo proteins are transferred into ESCPE-1 transport carriers to promote endosome-to-plasma membrane recycling. By revealing that assembly of the SNX27:Retromer:ESCPE-1 coat evolved in a stepwise manner during early metazoan evolution, likely reflecting the increasing complexity of endosome-to-plasma membrane recycling from the ancestral opisthokont to modern animals, we provide further evidence of the functional diversification of yeast pentameric Retromer in the recycling of hundreds of integral membrane proteins in metazoans.

Asteroid: a new algorithm to infer species trees from gene trees under high proportions of missing data.

MOTIVATION: Missing data and incomplete lineage sorting (ILS) are two major obstacles to accurate species tree inference. Gene tree summary methods such as ASTRAL and ASTRID have been developed to account for ILS. However, they can be severely affected by high levels of missing data. RESULTS: We present Asteroid, a novel algorithm that infers an unrooted species tree from a set of unrooted gene trees. We show on both empirical and simulated datasets that Asteroid is substantially more accurate than ASTRAL and ASTRID for very high proportions (>80%) of missing data. Asteroid is several orders of magnitude faster than ASTRAL for datasets that contain thousands of genes. It offers advanced features such as parallelization, support value computation and support for multi-copy and multifurcating gene trees. AVAILABILITY AND IMPLEMENTATION: Asteroid is freely available at https://github.com/BenoitMorel/Asteroid. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Signatures of hybridization in Trypanosoma brucei.

Genetic exchange among disease-causing micro-organisms can generate progeny that combine different pathogenic traits. Though sexual reproduction has been described in trypanosomes, its impact on the epidemiology of Human African Trypanosomiasis (HAT) remains controversial. However, human infective and non-human infective strains of Trypanosoma brucei circulate in the same transmission cycles in HAT endemic areas in subsaharan Africa, providing the opportunity for mating during the developmental cycle in the tsetse fly vector. Here we investigated inheritance among progeny from a laboratory cross of T. brucei and then applied these insights to genomic analysis of field-collected isolates to identify signatures of past genetic exchange. Genomes of two parental and four hybrid progeny clones with a range of DNA contents were assembled and analysed by k-mer and single nucleotide polymorphism (SNP) frequencies to determine heterozygosity and chromosomal inheritance. Variant surface glycoprotein (VSG) genes and kinetoplast (mitochondrial) DNA maxi- and minicircles were extracted from each genome to examine how each of these components was inherited in the hybrid progeny. The same bioinformatic approaches were applied to an additional 37 genomes representing the diversity of T. brucei in subsaharan Africa and T. evansi. SNP analysis provided evidence of crossover events affecting all 11 pairs of megabase chromosomes and demonstrated that polyploid hybrids were formed post-meiotically and not by fusion of the parental diploid cells. VSGs and kinetoplast DNA minicircles were inherited biparentally, with approximately equal numbers from each parent, whereas maxicircles were inherited uniparentally. Extrapolation of these findings to field isolates allowed us to distinguish clonal descent from hybridization by comparing maxicircle genotype to VSG and minicircle repertoires. Discordance between maxicircle genotype and VSG and minicircle repertoires indicated inter-lineage hybridization. Significantly, some of the hybridization events we identified involved human infective and non-human infective trypanosomes circulating in the same geographic areas.

Metagenome-assembled genome of the glacier alga Ancylonema yields insights into the evolution of streptophyte life on ice and land.

Contemporary glaciers are inhabited by streptophyte algae that balance photosynthesis and growth with tolerance of low temperature, desiccation and UV radiation. These same environmental challenges have been hypothesised as the driving force behind the evolution of land plants from streptophyte algal ancestors in the Cryogenian (720-635 million years ago). We sequenced, assembled and analysed the metagenome-assembled genome of the glacier alga Ancylonema nordenskiöldii to investigate its adaptations to life in ice, and whether this represents a vestige of Cryogenian exaptations. Phylogenetic analysis confirms the placement of glacier algae within the sister lineage to land plants, Zygnematophyceae. The metagenome-assembled genome is characterised by an expansion of genes involved in tolerance of high irradiance and UV light, while lineage-specific diversification is linked to the novel screening pigmentation of glacier algae. We found no support for the hypothesis of a common genomic basis for adaptations to ice and to land in streptophytes. Comparative genomics revealed that the reductive morphological evolution in the ancestor of Zygnematophyceae was accompanied by reductive genome evolution. This first genome-scale data for glacier algae suggests an Ancylonema-specific adaptation to the cryosphere, and sheds light on the genome evolution of land plants and Zygnematophyceae.

Recoding Amino Acids to a Reduced Alphabet may Increase or Decrease Phylogenetic Accuracy.

Common molecular phylogenetic characteristics such as long branches and compositional heterogeneity can be problematic for phylogenetic reconstruction when using amino acid data. Recoding alignments to reduced alphabets before phylogenetic analysis has often been used both to explore and potentially decrease the effect of such problems. We tested the effectiveness of this strategy on topological accuracy using simulated data on four-taxon trees. We simulated alignments in phylogenetically challenging ways to test the phylogenetic accuracy of analyses using various recoding strategies together with commonly used homogeneous models. We tested three recoding methods based on amino acid exchangeability, and another recoding method based on lowering the compositional heterogeneity among alignment sequences as measured by the Chi-squared statistic. Our simulation results show that on trees with long branches where sequences approach saturation, accuracy was not greatly affected by exchangeability-based recodings, but Chi-squared-based recoding decreased accuracy. We then simulated sequences with different kinds of compositional heterogeneity over the tree. Recoding often increased accuracy on such alignments. Exchangeability-based recoding was rarely worse than not recoding, and often considerably better. Recoding based on lowering the Chi-squared value improved accuracy in some cases but not in others, suggesting that low compositional heterogeneity by itself is not sufficient to increase accuracy in the analysis of these alignments. We also simulated alignments using site-specific amino acid profiles, making sequences that had compositional heterogeneity over alignment sites. Exchangeability-based recoding coupled with site-homogeneous models had poor accuracy for these data sets but Chi-squared-based recoding on these alignments increased accuracy. We then simulated data sets that were compositionally both site- and tree-heterogeneous, like many real data sets. The effect on the accuracy of recoding such doubly problematic data sets varied widely, depending on the type of compositional tree heterogeneity and on the recoding scheme. Interestingly, analysis of unrecoded compositionally heterogeneous alignments with the NDCH or CAT models was generally more accurate than homogeneous analysis, whether recoded or not. Overall, our results suggest that making trees for recoded amino acid data sets can be useful, but they need to be interpreted cautiously as part of a more comprehensive analysis. The use of better-fitting models like NDCH and CAT, which directly account for the patterns in the data, may offer a more promising long-term solution for analyzing empirical data. [Compositional heterogeneity; models of evolution; phylogenetic methods; recoding amino acid data sets.].

Recent expansion of metabolic versatility in Diplonema papillatum, the model species of a highly speciose group of marine eukaryotes.

BACKGROUND: Diplonemid flagellates are among the most abundant and species-rich of known marine microeukaryotes, colonizing all habitats, depths, and geographic regions of the world ocean. However, little is known about their genomes, biology, and ecological role. RESULTS: We present the first nuclear genome sequence from a diplonemid, the type species Diplonema papillatum. The ~ 280-Mb genome assembly contains about 32,000 protein-coding genes, likely co-transcribed in groups of up to 100. Gene clusters are separated by long repetitive regions that include numerous transposable elements, which also reside within introns. Analysis of gene-family evolution reveals that the last common diplonemid ancestor underwent considerable metabolic expansion. D. papillatum-specific gains of carbohydrate-degradation capability were apparently acquired via horizontal gene transfer. The predicted breakdown of polysaccharides including pectin and xylan is at odds with reports of peptides being the predominant carbon source of this organism. Secretome analysis together with feeding experiments suggest that D. papillatum is predatory, able to degrade cell walls of live microeukaryotes, macroalgae, and water plants, not only for protoplast feeding but also for metabolizing cell-wall carbohydrates as an energy source. The analysis of environmental barcode samples shows that D. papillatum is confined to temperate coastal waters, presumably acting in bioremediation of eutrophication. CONCLUSIONS: Nuclear genome information will allow systematic functional and cell-biology studies in D. papillatum. It will also serve as a reference for the highly diverse diplonemids and provide a point of comparison for studying gene complement evolution in the sister group of Kinetoplastida, including human-pathogenic taxa.

SpeciesRax: A Tool for Maximum Likelihood Species Tree Inference from Gene Family Trees under Duplication, Transfer, and Loss.

Species tree inference from gene family trees is becoming increasingly popular because it can account for discordance between the species tree and the corresponding gene family trees. In particular, methods that can account for multiple-copy gene families exhibit potential to leverage paralogy as informative signal. At present, there does not exist any widely adopted inference method for this purpose. Here, we present SpeciesRax, the first maximum likelihood method that can infer a rooted species tree from a set of gene family trees and can account for gene duplication, loss, and transfer events. By explicitly modeling events by which gene trees can depart from the species tree, SpeciesRax leverages the phylogenetic rooting signal in gene trees. SpeciesRax infers species tree branch lengths in units of expected substitutions per site and branch support values via paralogy-aware quartets extracted from the gene family trees. Using both empirical and simulated data sets we show that SpeciesRax is at least as accurate as the best competing methods while being one order of magnitude faster on large data sets at the same time. We used SpeciesRax to infer a biologically plausible rooted phylogeny of the vertebrates comprising 188 species from 31,612 gene families in 1 h using 40 cores. SpeciesRax is available under GNU GPL at https://github.com/BenoitMorel/GeneRax and on BioConda.

Relative Time Constraints Improve Molecular Dating.

Dating the tree of life is central to understanding the evolution of life on Earth. Molecular clocks calibrated with fossils represent the state of the art for inferring the ages of major groups. Yet, other information on the timing of species diversification can be used to date the tree of life. For example, horizontal gene transfer events and ancient coevolutionary interactions such as (endo)symbioses occur between contemporaneous species and thus can imply temporal relationships between two nodes in a phylogeny. Temporal constraints from these alternative sources can be particularly helpful when the geological record is sparse, for example, for microorganisms, which represent the majority of extant and extinct biodiversity. Here, we present a new method to combine fossil calibrations and relative age constraints to estimate chronograms. We provide an implementation of relative age constraints in RevBayes that can be combined in a modular manner with the wide range of molecular dating methods available in the software. We use both realistic simulations and empirical datasets of 40 Cyanobacteria and 62 Archaea to evaluate our method. We show that the combination of relative age constraints with fossil calibrations significantly improves the estimation of node ages. [Archaea, Bayesian analysis, cyanobacteria, dating, endosymbiosis, lateral gene transfer, MCMC, molecular clock, phylogenetic dating, relaxed molecular clock, revbayes, tree of life.].

Convergent Evolution of Hydrogenosomes from Mitochondria by Gene Transfer and Loss.

Hydrogenosomes are H2-producing mitochondrial homologs found in some anaerobic microbial eukaryotes that provide a rare intracellular niche for H2-utilizing endosymbiotic archaea. Among ciliates, anaerobic and aerobic lineages are interspersed, demonstrating that the switch to an anaerobic lifestyle with hydrogenosomes has occurred repeatedly and independently. To investigate the molecular details of this transition, we generated genomic and transcriptomic data sets from anaerobic ciliates representing three distinct lineages. Our data demonstrate that hydrogenosomes have evolved from ancestral mitochondria in each case and reveal different degrees of independent mitochondrial genome and proteome reductive evolution, including the first example of complete mitochondrial genome loss in ciliates. Intriguingly, the FeFe-hydrogenase used for generating H2 has a unique domain structure among eukaryotes and appears to have been present, potentially through a single lateral gene transfer from an unknown donor, in the common aerobic ancestor of all three lineages. The early acquisition and retention of FeFe-hydrogenase helps to explain the facility whereby mitochondrial function can be so radically modified within this diverse and ecologically important group of microbial eukaryotes.

Characterization and evolutionary origin of novel C2H2 zinc finger protein (ZNF648) required for both erythroid and megakaryocyte differentiation in humans.

Human ZNF648 is a novel poly C-terminal C2H2 zinc finger protein identified amongst the most dysregulated proteins in erythroid cells differentiated from iPSC. Its nuclear localisation and structure indicate it is likely a DNA-binding protein. Using a combination of ZNF648 overexpression in an iPSC line and primary adult erythroid cells, ZNF648 knockdown in primary adult erythroid cells and megakaryocytes, comparative proteomics and transcriptomics we show that ZNF648 is required for both erythroid and megakaryocyte differentiation. Orthologues of ZNF648 were detected across Mammals, Reptilia, Actinopterygii, in some Aves, Amphibia and Coelacanthiformes suggesting the gene originated in the common ancestor of Osteichthyes (Euteleostomi or bony fish). Conservation of the C-terminal zinc finger domain is higher, with some variation in zinc finger number but a core of at least six zinc fingers conserved across all groups, with the N-terminus recognisably similar within but not between major lineages. This suggests the N-terminus of ZNF648 evolves faster than the C-terminus, however this is not due to exon-shuffling as the entire coding region of ZNF648 is within a single exon. As for other such transcription factors, the N-terminus likely carries out regulatory functions, but showed no sequence similarity to any known domains. The greater functional constraint on the zinc finger domain suggests ZNF648 binds at least some similar regions of DNA in the different organisms. However, divergence of the N-terminal region may enable differential expression, allowing adaptation of function in the different organisms.

Integrative modeling of gene and genome evolution roots the archaeal tree of life.

A root for the archaeal tree is essential for reconstructing the metabolism and ecology of early cells and for testing hypotheses that propose that the eukaryotic nuclear lineage originated from within the Archaea; however, published studies based on outgroup rooting disagree regarding the position of the archaeal root. Here we constructed a consensus unrooted archaeal topology using protein concatenation and a multigene supertree method based on 3,242 single gene trees, and then rooted this tree using a recently developed model of genome evolution. This model uses evidence from gene duplications, horizontal transfers, and gene losses contained in 31,236 archaeal gene families to identify the most likely root for the tree. Our analyses support the monophyly of DPANN (Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanoarchaeota, Nanohaloarchaea), a recently discovered cosmopolitan and genetically diverse lineage, and, in contrast to previous work, place the tree root between DPANN and all other Archaea. The sister group to DPANN comprises the Euryarchaeota and the TACK Archaea, including Lokiarchaeum, which our analyses suggest are monophyletic sister lineages. Metabolic reconstructions on the rooted tree suggest that early Archaea were anaerobes that may have had the ability to reduce CO2 to acetate via the Wood-Ljungdahl pathway. In contrast to proposals suggesting that genome reduction has been the predominant mode of archaeal evolution, our analyses infer a relatively small-genomed archaeal ancestor that subsequently increased in complexity via gene duplication and horizontal gene transfer.

The Effect of Nonreversibility on Inferring Rooted Phylogenies.

Most phylogenetic models assume that the evolutionary process is stationary and reversible. In addition to being biologically improbable, these assumptions also impair inference by generating models under which the likelihood does not depend on the position of the root. Consequently, the root of the tree cannot be inferred as part of the analysis. Yet identifying the root position is a key component of phylogenetic inference because it provides a point of reference for polarizing ancestor-descendant relationships and therefore interpreting the tree. In this paper, we investigate the effect of relaxing the unrealistic reversibility assumption and allowing the position of the root to be another unknown. We propose two hierarchical models that are centered on a reversible model but perturbed to allow nonreversibility. The models differ in the degree of structure imposed on the perturbations. The analysis is performed in the Bayesian framework using Markov chain Monte Carlo methods for which software is provided. We illustrate the performance of the two nonreversible models in analyses of simulated data using two types of topological priors. We then apply the models to a real biological data set, the radiation of polyploid yeasts, for which there is robust biological opinion about the root position. Finally, we apply the models to a second biological alignment for which the rooted tree is controversial: the ribosomal tree of life. We compare the two nonreversible models and conclude that both are useful in inferring the position of the root from real biological data.

An archaeal origin of eukaryotes supports only two primary domains of life.

The discovery of the Archaea and the proposal of the three-domains 'universal' tree, based on ribosomal RNA and core genes mainly involved in protein translation, catalysed new ideas for cellular evolution and eukaryotic origins. However, accumulating evidence suggests that the three-domains tree may be incorrect: evolutionary trees made using newer methods place eukaryotic core genes within the Archaea, supporting hypotheses in which an archaeon participated in eukaryotic origins by founding the host lineage for the mitochondrial endosymbiont. These results provide support for only two primary domains of life--Archaea and Bacteria--because eukaryotes arose through partnership between them.

Molecular palaeontology illuminates the evolution of ecdysozoan vision.

Colour vision is known to have arisen only twice-once in Vertebrata and once within the Ecdysozoa, in Arthropoda. However, the evolutionary history of ecdysozoan vision is unclear. At the molecular level, visual pigments, composed of a chromophore and a protein belonging to the opsin family, have different spectral sensitivities and these mediate colour vision. At the morphological level, ecdysozoan vision is conveyed by eyes of variable levels of complexity; from the simple ocelli observed in the velvet worms (phylum Onychophora) to the marvellously complex eyes of insects, spiders, and crustaceans. Here, we explore the evolution of ecdysozoan vision at both the molecular and morphological level; combining analysis of a large-scale opsin dataset that includes previously unknown ecdysozoan opsins with morphological analyses of key Cambrian fossils with preserved eye structures. We found that while several non-arthropod ecdysozoan lineages have multiple opsins, arthropod multi-opsin vision evolved through a series of gene duplications that were fixed in a period of 35-71 million years (Ma) along the stem arthropod lineage. Our integrative study of the fossil and molecular record of vision indicates that fossils with more complex eyes were likely to have possessed a larger complement of opsin genes.

Clarifying the Relationships between Microsporidia and Cryptomycota.

Some protists with microsporidian-like cell biological characters, including Mitosporidium, Paramicrosporidium, and Nucleophaga, have SSU rRNA gene sequences that are much less divergent than canonical Microsporidia. We analysed the phylogenetic placement and environmental diversity of microsporidian-like lineages that group near the base of the fungal radiation and show that they group in a clade with metchnikovellids and canonical microsporidians, to the exclusion of the clade including Rozella, in line with what is currently known of their morphology and cell biology. These results show that the phylogenetic scope of Microsporidia has been greatly underestimated. We propose that much of the lineage diversity previously thought to be cryptomycotan/rozellid is actually microsporidian, offering new insights into the evolution of the highly specialized parasitism of canonical Microsporidia. This insight has important implications for our understanding of opisthokont evolution and ecology, and is important for accurate interpretation of environmental diversity. Our analyses also demonstrate that many opisthosporidian (aphelid+rozellid+microsporidian) SSU V4 OTUs from Neotropical forest soils group with the short-branching Microsporidia, consistent with the abundance of their protist and arthropod hosts in soils. This novel diversity of Microsporidia provides a unique opportunity to investigate the evolutionary origins of a highly specialized clade of major animal parasites.

Coupling of diversification and pH adaptation during the evolution of terrestrial Thaumarchaeota.

The Thaumarchaeota is an abundant and ubiquitous phylum of archaea that plays a major role in the global nitrogen cycle. Previous analyses of the ammonia monooxygenase gene amoA suggest that pH is an important driver of niche specialization in these organisms. Although the ecological distribution and ecophysiology of extant Thaumarchaeota have been studied extensively, the evolutionary rise of these prokaryotes to ecological dominance in many habitats remains poorly understood. To characterize processes leading to their diversification, we investigated coevolutionary relationships between amoA, a conserved marker gene for Thaumarchaeota, and soil characteristics, by using deep sequencing and comprehensive environmental data in Bayesian comparative phylogenetics. These analyses reveal a large and rapid increase in diversification rates during early thaumarchaeotal evolution; this finding was verified by independent analyses of 16S rRNA. Our findings suggest that the entire Thaumarchaeota diversification regime was strikingly coupled to pH adaptation but less clearly correlated with several other tested environmental factors. Interestingly, the early radiation event coincided with a period of pH adaptation that enabled the terrestrial Thaumarchaeota ancestor to initially move from neutral to more acidic and alkaline conditions. In contrast to classic evolutionary models, whereby niches become rapidly filled after adaptive radiation, global diversification rates have remained stably high in Thaumarchaeota during the past 400-700 million years, suggesting an ongoing high rate of niche formation or switching for these microbes. Our study highlights the enduring importance of environmental adaptation during thaumarchaeotal evolution and, to our knowledge, is the first to link evolutionary diversification to environmental adaptation in a prokaryotic phylum.