Gene duplication drives genome expansion in a major lineage of Thaumarchaeota

An evolutionary view of the Fusarium core genome

Fusarium, a member of the Ascomycota fungi, encompasses several pathogenic species significant to plants and animals. Some phytopathogenic species have received special attention due to their negative economic impact on the agricultural industry around the world. Traditionally, identification and taxonomic analysis of Fusarium have relied on morphological and phenotypic features, including the fungal host, leading to taxonomic conflicts that have been solved using molecular systematic technologies. In this work, we applied a phylogenomic approach that allowed us to resolve the evolutionary history of the species complexes of the genus and present evidence that supports the F. ventricosum species complex as the most basal lineage of the genus. Additionally, we present evidence that proposes modifications to the previous hypothesis of the evolutionary history of the F. staphyleae, F. newnesense, F. nisikadoi, F. oxysporum, and F. fujikuroi species complexes. Evolutionary analysis showed that the genome GC content tends to be lower in more modern lineages, in both, the whole-genome and core-genome coding DNA sequences. In contrast, genome size gain and losses are present during the evolution of the genus. Interestingly, core genome duplication events positively correlate with genome size. Evolutionary and genome conservation analysis supports the F3 hypothesis of Fusarium as a more compact and conserved group in terms of genome conservation. By contrast, outside of the F3 hypothesis, the most basal clades only share 8.8% of its genomic sequences with the F3 clade.

Comparative genomics reveals the adaptation of ammonia-oxidising Thaumarchaeota to arid soils

Thaumarchaeota are predominant in oligotrophic habitats such as deserts and arid soils, but their adaptations to these arid conditions are not well understood. In this study, we assembled 23 Thaumarchaeota genomes from arid and semi-arid soils collected from the Inner Mongolia Steppe and the Qinghai-Tibet Plateau. Using a comparative genomics approach, integrated with 614 Thaumarchaeota genomes from public databases, we identified the traits and evolutionary forces that contribute to their adaptations to aridity. Our results showed that the newly assembled genomes represent an early diverging group within the lineage of ammonia-oxidising Thaumarchaeota. While the genomic functions previously identified in arid soil lineages were conserved across terrestrial, shallow-ocean and deep-ocean lineages, several traits likely contribute to Thaumarchaeota's adaptation to aridity. These include chlorite dismutase, arsenate reductase, V-type ATPase and genes dealing with oxidative stresses. The acquisition and loss of traits at the last common ancestor of arid soil lineages may have facilitated the specialisation of Thaumarchaeota in arid soils. Additionally, the acquisition of unique adaptive traits, such as a urea transporter, Ca2+ :H+ antiporter, mannosyl-3-phosphoglycerate synthase and phosphatase, DNA end-binding protein Ku and phage shock protein A, further distinguishes arid soil Thaumarchaeota. This study provides evidence for the adaptations of Thaumarchaeota to arid soil, enhancing our understanding of the nitrogen and carbon cycling driven by Thaumarchaeota in drylands.

Gene inversion led to the emergence of brackish archaeal heterotrophs in the aftermath of the Cryogenian Snowball Earth

Land-ocean interactions greatly impact the evolution of coastal life on earth. However, the ancient geological forces and genetic mechanisms that shaped evolutionary adaptations and allowed microorganisms to inhabit coastal brackish waters remain largely unexplored. In this study, we infer the evolutionary trajectory of the ubiquitous heterotrophic archaea Poseidoniales (Marine Group II archaea) presently occurring across global aquatic habitats. Our results show that their brackish subgroups had a single origination, dated to over 600 million years ago, through the inversion of the magnesium transport gene corA that conferred osmotic-stress tolerance. The subsequent loss and gain of corA were followed by genome-wide adjustment, characterized by a general two-step mode of selection in microbial speciation. The coastal family of Poseidoniales showed a rapid increase in the evolutionary rate during and in the aftermath of the Cryogenian Snowball Earth (∼700 million years ago), possibly in response to the enhanced phosphorus supply and the rise of algae. Our study highlights the close interplay between genetic changes and ecosystem evolution that boosted microbial diversification in the Neoproterozoic continental margins, where the Cambrian explosion of animals soon followed.

Metabolic potential of Nitrososphaera-associated clades

Soil ammonia-oxidizing archaea (AOA) play a crucial role in converting ammonia to nitrite, thereby mobilizing reactive nitrogen species into their soluble form, with a significant impact on nitrogen losses from terrestrial soils. Yet, our knowledge regarding their diversity and functions remains limited. In this study, we reconstructed 97 high-quality AOA metagenome-assembled genomes (MAGs) from 180 soil samples collected in Central Germany during 2014-2019 summers. These MAGs were affiliated with the order Nitrososphaerales and clustered into four family-level clades (NS-α/γ/δ/ε). Among these MAGs, 75 belonged to the most abundant but least understood δ-clade. Within the δ-clade, the amoA genes in three MAGs from neutral soils showed a 99.5% similarity to the fosmid clone 54d9, which has served as representative of the δ-clade for the past two decades since even today no cultivated representatives are available. Seventy-two MAGs constituted a distinct δ sub-clade, and their abundance and expression activity were more than twice that of other MAGs in slightly acidic soils. Unlike the less abundant clades (α, γ, and ε), the δ-MAGs possessed multiple highly expressed intracellular and extracellular carbohydrate-active enzymes responsible for carbohydrate binding (CBM32) and degradation (GH5), along with highly expressed genes involved in ammonia oxidation. Together, these results suggest metabolic versatility of uncultured soil AOA and a potential mixotrophic or chemolithoheterotrophic lifestyle among 54d9-like AOA.

Novel order-level lineage of ammonia-oxidizing archaea widespread in marine and terrestrial environments

Ammonia-oxidizing archaea (AOA) are among the most ubiquitous and abundant archaea on Earth, widely distributed in marine, terrestrial, and geothermal ecosystems. However, the genomic diversity, biogeography, and evolutionary process of AOA populations in subsurface environments are vastly understudied compared to those in marine and soil systems. Here, we report a novel AOA order Candidatus (Ca.) Nitrosomirales which forms a sister lineage to the thermophilic Ca. Nitrosocaldales. Metagenomic and 16S rRNA gene-read mapping demonstrates the abundant presence of Nitrosomirales AOA in various groundwater environments and their widespread distribution across a range of geothermal, terrestrial, and marine habitats. Terrestrial Nitrosomirales AOA show the genetic capacity of using formate as a source of reductant and using nitrate as an alternative electron acceptor. Nitrosomirales AOA appear to have acquired key metabolic genes and operons from other mesophilic populations via horizontal gene transfer, including genes encoding urease, nitrite reductase, and V-type ATPase. The additional metabolic versatility conferred by acquired functions may have facilitated their radiation into a variety of subsurface, marine, and soil environments. We also provide evidence that each of the four AOA orders spans both marine and terrestrial habitats, which suggests a more complex evolutionary history for major AOA lineages than previously proposed. Together, these findings establish a robust phylogenomic framework of AOA and provide new insights into the ecology and adaptation of this globally abundant functional guild.

Temperature, pH, and oxygen availability contributed to the functional differentiation of ancient Nitrososphaeria

Ammonia-oxidizing Nitrososphaeria are among the most abundant archaea on Earth and have profound impacts on the biogeochemical cycles of carbon and nitrogen. In contrast to these well-studied ammonia-oxidizing archaea (AOA), deep-branching non-AOA within this class remain poorly characterized because of a low number of genome representatives. Here, we reconstructed 128 Nitrososphaeria metagenome-assembled genomes from acid mine drainage and hot spring sediment metagenomes. Comparative genomics revealed that extant non-AOA are functionally diverse, with capacity for carbon fixation, carbon monoxide oxidation, methanogenesis, and respiratory pathways including oxygen, nitrate, sulfur, or sulfate, as potential terminal electron acceptors. Despite their diverse anaerobic pathways, evolutionary history inference suggested that the common ancestor of Nitrososphaeria was likely an aerobic thermophile. We further surmise that the functional differentiation of Nitrososphaeria was primarily shaped by oxygen, pH, and temperature, with the acquisition of pathways for carbon, nitrogen, and sulfur metabolism. Our study provides a more holistic and less biased understanding of the diversity, ecology, and deep evolution of the globally abundant Nitrososphaeria.

Metagenomic characterization of a novel non-ammonia-oxidizing Thaumarchaeota from hadal sediment

BACKGROUND

The hadal sediment, found at an ocean depth of more than 6000 m, is geographically isolated and under extremely high hydrostatic pressure, resulting in a unique ecosystem. Thaumarchaeota are ubiquitous marine microorganisms predominantly present in hadal environments. While there have been several studies on Thaumarchaeota there, most of them have primarily focused on ammonia-oxidizing archaea (AOA). However, systematic metagenomic research specifically targeting heterotrophic non-AOA Thaumarchaeota is lacking.

RESULTS

In this study, we explored the metagenomes of Challenger Deep hadal sediment, focusing on the Thaumarchaeota. Functional analysis of sequence reads revealed the potential contribution of Thaumarchaeota to recalcitrant dissolved organic matter degradation. Metagenome assembly binned one new group of hadal sediment-specific and ubiquitously distributed non-AOA Thaumarchaeota, named Group-3.unk. Pathway reconstruction of this new type of Thaumarchaeota also supports heterotrophic characteristics of Group-3.unk, along with ABC transporters for the uptake of amino acids and carbohydrates and catabolic utilization of these substrates. This new clade of Thaumarchaeota also contains aerobic oxidation of carbon monoxide-related genes. Complete glyoxylate cycle is a distinctive feature of this clade in supplying intermediates of anabolic pathways. The pan-genomic and metabolic analyses of metagenome-assembled genomes belonging to Group-3.unk Thaumarchaeota have highlighted distinctions, including the dihydroxy phthalate decarboxylase gene associated with the degradation of aromatic compounds and the absence of genes related to the synthesis of some types of vitamins compared to AOA. Notably, Group-3.unk shares a common feature with deep ocean AOA, characterized by their high hydrostatic pressure resistance, potentially associated with the presence of V-type ATP and di-myo-inositol phosphate syntheses-related genes. The enrichment of organic matter in hadal sediments might be attributed to the high recruitment of sequence reads of the Group-3.unk clade of heterotrophic Thaumarchaeota in the trench sediment. Evolutionary and genetic dynamic analyses suggest that Group-3 non-AOA consists of mesophilic Thaumarchaeota organisms. These results indicate a potential role in the transition from non-AOA to AOA Thaumarchaeota and from thermophilic to mesophilic Thaumarchaeota, shedding light on recent evolutionary pathways.

CONCLUSIONS

One novel clade of heterotrophic non-AOA Thaumarchaeota was identified through metagenome analysis of sediments from Challenger Deep. Our study provides insight into the ecology and genomic characteristics of the new sub-group of heterotrophic non-AOA Thaumarchaeota, thereby extending the knowledge of the evolution of Thaumarchaeota. Video Abstract.

Genomic insight into strategy, interaction and evolution of nitrifiers in metabolizing key labile-dissolved organic nitrogen in different environmental niches

Ammonia-oxidizing archaea (AOA) and bacteria (AOB), nitrite-oxidizing bacteria (NOB), and complete ammonia oxidizers (comammox) are responsible for nitrification in nature; however, some groups have been reported to utilize labile-dissolved organic nitrogen (LDON) for satisfying nitrogen demands. To understand the universality of their capacity of LDON metabolism, we collected 70 complete genomes of AOA, AOB, NOB, and comammox from typical environments for exploring their potentials in the metabolism of representative LDON (urea, polyamines, cyanate, taurine, glycine betaine, and methylamine). Genomic analyses showed that urea was the most popular LDON used by nitrifiers. Each group harbored unique urea transporter genes (AOA: dur3 and utp, AOB: utp, and NOB and comammox: urtABCDE and utp) accompanied by urease genes ureABC. The differentiation in the substrate affinity of these transporters implied the divergence of urea utilization efficiency in nitrifiers, potentially driving them into different niches. The cyanate transporter (cynABD and focA/nirC) and degradation (cynS) genes were detected mostly in NOB, indicating their preference for a wide range of nitrogen substrates to satisfy high nitrogen demands. The lack of genes involved in the metabolism of polyamines, taurine, glycine betaine, and methylamines in most of nitrifiers suggested that they were not able to serve as a source of ammonium, only if they were degraded or oxidized extracellularly as previously reported. The phylogenetic analyses assisted with comparisons of GC% and the Codon Adaptation Index between target genes and whole genomes of nitrifiers implied that urea metabolic genes dur3 and ureC in AOA evolved independently from bacteria during the transition from Thaumarchaeota to AOA, while utp in terrestrial AOA was acquired from bacteria via lateral gene transfer (LGT). Cyanate transporter genes cynS and focA/nirC detected only in a terrestrial AOA Candidadus Nitrsosphaera gargensis Ga9.2 could be gained synchronously with Nitrospira of NOB by an ancient LGT. Our results indicated that LDON utilization was a common feature in nitrifiers, but metabolic potentials were different among nitrifiers, possibly being intensely interacted with their niches, survival strategies, and evolutions.

Genomics of soil depth niche partitioning in the Thaumarchaeota family Gagatemarchaeaceae

Knowledge of deeply-rooted non-ammonia oxidising Thaumarchaeota lineages from terrestrial environments is scarce, despite their abundance in acidic soils. Here, 15 new deeply-rooted thaumarchaeotal genomes were assembled from acidic topsoils (0-15 cm) and subsoils (30-60 cm), corresponding to two genera of terrestrially prevalent Gagatemarchaeaceae (previously known as thaumarchaeotal Group I.1c) and to a novel genus of heterotrophic terrestrial Thaumarchaeota. Unlike previous predictions, metabolic annotations suggest Gagatemarchaeaceae perform aerobic respiration and use various organic carbon sources. Evolutionary divergence between topsoil and subsoil lineages happened early in Gagatemarchaeaceae history, with significant metabolic and genomic trait differences. Reconstruction of the evolutionary mechanisms showed that the genome expansion in topsoil Gagatemarchaeaceae resulted from extensive early lateral gene acquisition, followed by progressive gene duplication throughout evolutionary history. Ancestral trait reconstruction using the expanded genomic diversity also did not support the previous hypothesis of a thermophilic last common ancestor of the ammonia-oxidising archaea. Ultimately, this study provides a good model for studying mechanisms driving niche partitioning between spatially related ecosystems.

Nitrogen and phosphorous acquisition strategies drive coexistence patterns among archaeal lineages in soil

Soil represents the largest reservoir of Archaea on Earth. Present-day archaeal diversity in soils globally is dominated by members of the class Nitrososphaeria. The evolutionary radiation of this class is thought to reflect adaptations to a wide range of temperatures, pH, and other environmental conditions. However, the mechanisms that govern competition and coexistence among Nitrososphaeria lineages in soil remain poorly understood. Here we show that predominant soil Nitrososphaeria lineages compose a patchwork of gene inventory and expression profiles for ammonia, urea, and phosphate utilization. In contrast, carbon fixation, respiration, and ATP synthesis genes are conserved and expressed consistently among predominant phylotypes across 12 major evolutionary lineages commonly found in soil. In situ gene expression profiles closely resemble pure culture reference strains under optimal growth conditions. Together, these results reveal resource-based coexistence patterns among Nitrososphaeria lineages and suggest complementary ecophysiological niches associated with differential nutrient acquisition strategies among globally predominant archaeal lineages in soil.

Habitat Transition in the Evolution of Bacteria and Archaea

Related groups of microbes are widely distributed across Earth's habitats, implying numerous dispersal and adaptation events over evolutionary time. However, relatively little is known about the characteristics and mechanisms of these habitat transitions, particularly for populations that reside in animal microbiomes. Here, we review the literature concerning habitat transitions among a variety of bacterial and archaeal lineages, considering the frequency of migration events, potential environmental barriers, and mechanisms of adaptation to new physicochemical conditions, including the modification of protein inventories and other genomic characteristics. Cells dependent on microbial hosts, particularly bacteria from the Candidate Phyla Radiation, have undergone repeated habitat transitions from environmental sources into animal microbiomes. We compare their trajectories to those of both free-living cells-including the Melainabacteria, Elusimicrobia, and methanogenic archaea-and cellular endosymbionts and bacteriophages, which have made similar transitions. We conclude by highlighting major related topics that may be worthy of future study.

Distribution and abundance of tetraether lipid cyclization genes in terrestrial hot springs reflect pH

Many Archaea produce membrane-spanning lipids that enable life in extreme environments. These isoprenoid glycerol dibiphytanyl glycerol tetraethers (GDGTs) may contain up to eight cyclopentyl and one cyclohexyl ring, where higher degrees of cyclization are associated with more acidic, hotter or energy-limited conditions. Recently, the genes encoding GDGT ring synthases, grsAB, were identified in two Sulfolobaceae; however, the distribution and abundance of grs homologs across environments inhabited by these and related organisms remain a mystery. To address this, we examined the distribution of grs homologs in relation to environmental temperature and pH, from thermal springs across Earth, where sequences derive from metagenomes, metatranscriptomes, single-cell and cultivar genomes. The abundance of grs homologs shows a strong negative correlation to pH, but a weak positive correlation to temperature. Archaeal genomes and metagenome-assembled genomes (MAGs) that carry two or more grs copies are more abundant in low pH springs. We also find grs in 12 archaeal classes, with the most representatives in Thermoproteia, followed by MAGs of the uncultured Korarchaeia, Bathyarchaeia and Hadarchaeia, while several Nitrososphaeria encodes >3 copies. Our findings highlight the key role of grs-catalysed lipid cyclization in archaeal diversification across hot and acidic environments.

Replicated life-history patterns and subsurface origins of the bacterial sister phyla Nitrospirota and Nitrospinota

The phyla Nitrospirota and Nitrospinota have received significant research attention due to their unique nitrogen metabolisms important to biogeochemical and industrial processes. These phyla are common inhabitants of marine and terrestrial subsurface environments and contain members capable of diverse physiologies in addition to nitrite oxidation and complete ammonia oxidation. Here, we use phylogenomics and gene-based analysis with ancestral state reconstruction and gene-tree-species-tree reconciliation methods to investigate the life histories of these two phyla. We find that basal clades of both phyla primarily inhabit marine and terrestrial subsurface environments. The genomes of basal clades in both phyla appear smaller and more densely coded than the later-branching clades. The extant basal clades of both phyla share many traits inferred to be present in their respective common ancestors, including hydrogen, one-carbon, and sulfur-based metabolisms. Later-branching groups, namely the more frequently studied classes Nitrospiria and Nitrospinia, are both characterized by genome expansions driven by either de novo origination or laterally transferred genes that encode functions expanding their metabolic repertoire. These expansions include gene clusters that perform the unique nitrogen metabolisms that both phyla are most well known for. Our analyses support replicated evolutionary histories of these two bacterial phyla, with modern subsurface environments representing a genomic repository for the coding potential of ancestral metabolic traits.

Unexpected complexity of the ammonia monooxygenase in archaea

Ammonia oxidation, as the first step of nitrification, constitutes a critical process in the global nitrogen cycle. However, fundamental knowledge of its key enzyme, the copper-dependent ammonia monooxygenase, is lacking, in particular for the environmentally abundant ammonia-oxidizing archaea (AOA). Here the structure of the enzyme is investigated by blue-native gel electrophoresis and proteomics from native membrane complexes of two AOA. Besides the known AmoABC subunits and the earlier predicted AmoX, two new protein subunits, AmoY and AmoZ, were identified. They are unique to AOA, highly conserved and co-regulated, and their genes are linked to other AMO subunit genes in streamlined AOA genomes. Modeling and in-gel cross-link approaches support an overall protomer structure similar to the distantly related bacterial particulate methane monooxygenase but also reveals clear differences in extracellular domains of the enzyme. These data open avenues for further structure-function studies of this ecologically important nitrification complex.

Postglacial adaptations enabled colonization and quasi-clonal dispersal of ammonia-oxidizing archaea in modern European large lakes

Ammonia-oxidizing archaea (AOA) play a key role in the aquatic nitrogen cycle. Their genetic diversity is viewed as the outcome of evolutionary processes that shaped ancestral transition from terrestrial to marine habitats. However, current genome-wide insights into AOA evolution rarely consider brackish and freshwater representatives or provide their divergence timeline in lacustrine systems. An unbiased global assessment of lacustrine AOA diversity is critical for understanding their origins, dispersal mechanisms, and ecosystem roles. Here, we leveraged continental-scale metagenomics to document that AOA species diversity in freshwater systems is remarkably low compared to marine environments. We show that the uncultured freshwater AOA, "Candidatus Nitrosopumilus limneticus," is ubiquitous and genotypically static in various large European lakes where it evolved 13 million years ago. We find that extensive proteome remodeling was a key innovation for freshwater colonization of AOA. These findings reveal the genetic diversity and adaptive mechanisms of a keystone species that has survived clonally in lakes for millennia.

Fine-scale evaluation of two standard 16S rRNA gene amplicon primer pairs for analysis of total prokaryotes and archaeal nitrifiers in differently managed soils

The advance of high-throughput molecular biology tools allows in-depth profiling of microbial communities in soils, which possess a high diversity of prokaryotic microorganisms. Amplicon-based sequencing of 16S rRNA genes is the most common approach to studying the richness and composition of soil prokaryotes. To reliably detect different taxonomic lineages of microorganisms in a single soil sample, an adequate pipeline including DNA isolation, primer selection, PCR amplification, library preparation, DNA sequencing, and bioinformatic post-processing is required. Besides DNA sequencing quality and depth, the selection of PCR primers and PCR amplification reactions arguably have the largest influence on the results. This study tested the performance and potential bias of two primer pairs, i.e., 515F (Parada)-806R (Apprill) and 515F (Parada)-926R (Quince) in the standard pipelines of 16S rRNA gene Illumina amplicon sequencing protocol developed by the Earth Microbiome Project (EMP), against shotgun metagenome-based 16S rRNA gene reads. The evaluation was conducted using five differently managed soils. We observed a higher richness of soil total prokaryotes by using reverse primer 806R compared to 926R, contradicting to in silico evaluation results. Both primer pairs revealed various degrees of taxon-specific bias compared to metagenome-derived 16S rRNA gene reads. Nonetheless, we found consistent patterns of microbial community variation associated with different land uses, irrespective of primers used. Total microbial communities, as well as ammonia oxidizing archaea (AOA), the predominant ammonia oxidizers in these soils, shifted along with increased soil pH due to agricultural management. In the unmanaged low pH plot abundance of AOA was dominated by the acid-tolerant NS-Gamma clade, whereas limed agricultural plots were dominated by neutral-alkaliphilic NS-Delta/NS-Alpha clades. This study stresses how primer selection influences community composition and highlights the importance of primer selection for comparative and integrative studies, and that conclusions must be drawn with caution if data from different sequencing pipelines are to be compared.

Comparative Genomic Insights into the Evolution of Halobacteria-Associated "Candidatus Nanohaloarchaeota"

Members of the phylum "Candidatus Nanohaloarchaeota," a representative lineage within the DPANN superphylum, are characterized by their nanosized cells and symbiotic lifestyle with Halobacteria. However, the development of the symbiosis remains unclear. Here, we propose two novel families, "Candidatus Nanoanaerosalinaceae" and "Candidatus Nanohalalkaliarchaeaceae" in "Ca. Nanohaloarchaeota," represented by five dereplicated metagenome-assembled genomes obtained from hypersaline sediments or related enrichment cultures of soda-saline lakes. Phylogenetic analyses reveal that the two novel families are placed at the root of the family "Candidatus Nanosalinaceae," including the cultivated taxa. The two novel families prefer hypersaline sediments, and the acid shift of predicted proteomes indicates a "salt-in" strategy for hypersaline adaptation. They contain a lower proportion of putative horizontal gene transfers from Halobacteria than "Ca. Nanosalinaceae," suggesting a weaker association with Halobacteria. Functional prediction and historical events reconstruction disclose that they exhibit divergent potentials in carbohydrate and organic acid metabolism and environmental responses. Globally, comparative genomic analyses based on the new families enrich the taxonomic and functional diversity of "Ca. Nanohaloarchaeota" and provide insights into the evolutionary process of "Ca. Nanohaloarchaeota" and their symbiotic relationship with Halobacteria. IMPORTANCE The DPANN superphylum is a group of archaea widely distributed in various habitats. They generally have small cells and have a symbiotic lifestyle with other archaea. The archaeal symbiotic interaction is vital to understanding microbial communities. However, the formation and evolution of the symbiosis between the DPANN lineages and other diverse archaea remain unclear. Based on phylogeny, habitat distribution, hypersaline adaptation, host prediction, functional potentials, and historical events of "Ca. Nanohaloarchaeota," a representative phylum within the DPANN superphylum, we report two novel families representing intermediate stages, and we infer the evolutionary process of "Ca. Nanohaloarchaeota" and their Halobacteria-associated symbiosis. Altogether, this research helps in understanding the evolution of symbiosis in "Ca. Nanohaloarchaeota" and provides a model for the evolution of other DPANN lineages.

Comparative Genomic Insights into Chemoreceptor Diversity and Habitat Adaptation of Archaea

Archaea are capable of sensing and responding to environmental changes by several signal transduction systems with different mechanisms. Much attention is paid to model organisms with complex signaling networks to understand their composition and function, but general principles regarding how an archaeal species organizes its chemoreceptor diversity and habitat adaptation are poorly understood.

Recovery of Lutacidiplasmatales archaeal order genomes suggests convergent evolution in Thermoplasmatota

The Terrestrial Miscellaneous Euryarchaeota Group has been identified in various environments, and the single genome investigated thus far suggests that these archaea are anaerobic sulfite reducers. We assemble 35 new genomes from this group that, based on genome analysis, appear to possess aerobic and facultative anaerobic lifestyles and may oxidise rather than reduce sulfite. We propose naming this order (representing 16 genera) "Lutacidiplasmatales" due to their occurrence in various acidic environments and placement within the phylum Thermoplasmatota. Phylum-level analysis reveals that Thermoplasmatota evolution had been punctuated by several periods of high levels of novel gene family acquisition. Several essential metabolisms, such as aerobic respiration and acid tolerance, were likely acquired independently by divergent lineages through convergent evolution rather than inherited from a common ancestor. Ultimately, this study describes the terrestrially prevalent Lutacidiciplasmatales and highlights convergent evolution as an important driving force in the evolution of archaeal lineages.

Pangenome evolution in environmentally transmitted symbionts of deep-sea mussels is governed by vertical inheritance

Microbial pangenomes vary across species; their size and structure are determined by genetic diversity within the population and by gene loss and horizontal gene transfer (HGT). Many bacteria are associated with eukaryotic hosts where the host colonization dynamics may impact bacterial genome evolution. Host-associated lifestyle has been recognized as a barrier to HGT in parentally transmitted bacteria. However, pangenome evolution of environmentally acquired symbionts remains understudied, often due to limitations in symbiont cultivation. Using high-resolution metagenomics, here we study pangenome evolution of two co-occurring endosymbionts inhabiting Bathymodiolus brooksi mussels from a single cold seep. The symbionts, sulfur-oxidizing (SOX) and methane-oxidizing (MOX) gamma-proteobacteria, are environmentally acquired at an early developmental stage and individual mussels may harbor multiple strains of each symbiont species. We found differences in the accessory gene content of both symbionts across individual mussels, which are reflected by differences in symbiont strain composition. Compared to core genes, accessory genes are enriched in genome plasticity functions. We found no evidence for recent horizontal gene transfer between both symbionts. A comparison between the symbiont pangenomes revealed that the MOX population is less diverged and contains fewer accessory genes, supporting that the MOX association with B. brooksi is more recent in comparison to that of SOX. Our results show that the pangenomes of both symbionts evolved mainly by vertical inheritance. We conclude that genome evolution of environmentally transmitted symbionts that associate with individual hosts over their lifetime is affected by a narrow symbiosis where the frequency of HGT is constrained.​.

Effects of Strigolactone on Torreya grandis Gene Expression and Soil Microbial Community Structure Under Simulated Nitrogen Deposition

Nitrogen enters the terrestrial ecosystem through deposition. High nitrogen levels can affect physical and chemical properties of soil and inhibit normal growth and reproduction of forest plants. Nitrogen modulates the composition of soil microorganisms. Strigolactones inhibits plant branching, promotes root growth, nutrient absorption, and promotes arbuscular fungal mycelia branching. Plants are subjected to increasing atmospheric nitrogen deposition. Therefore, it is imperative to explore the relationship between strigolactone and nitrogen deposition of plants and abundance of soil microorganisms. In the present study, the effects of strigolactone on genetic responses and soil microorganisms of Torreya grandis, under simulated nitrogen deposition were explored using high-throughput sequencing techniques. T. grandis is a subtropical economic tree species in China. A total of 4,008 differentially expressed genes were identified in additional N deposition and GR24 treatment. These genes were associated with multiple GO terms and metabolic pathways. GO enrichment analysis showed that several DEGs were associated with enrichment of the transporter activity term. Both additional nitrogen deposition and GR24 treatment modulated the content of nutrient elements. The content of K reduced in leaves after additional N deposition treatment. The content of P increased in leaves after GR24 treatment. A total of 20 families and 29 DEGs associated with transporters were identified. These transporters may be regulated by transcription factors. A total of 1,402,819 clean reads and 1,778 amplicon sequence variants (ASVs) were generated through Bacterial 16S rRNA sequencing. Random forest classification revealed that Legionella, Lacunisphaera, Klebsiella, Bryobacter, and Janthinobacterium were significantly enriched in the soil in the additional N deposition group and the GR24 treatment group. Co-occurrence network analysis showed significant differences in composition of soil microbial community under different treatments. These results indicate a relationship between N deposition and strigolactones effect. The results provide new insights on the role of strigolactones in plants and composition of soil microorganisms under nitrogen deposition.

Phylogenetic divergence and adaptation of Nitrososphaeria across lake depths and freshwater ecosystems

Thaumarchaeota (now the class Nitrososphaeria in the phylum Thermoproteota in GTDB taxonomy) are abundant across marine and soil habitats; however, their genomic diversity and evolutionary history in freshwater environments remain elusive. Here, we reconstructed 17 high-quality metagenome-assembled genomes of Nitrososphaeria from a deep lake and two great rivers, and compared all available genomes between freshwater and marine habitats regarding their phylogenetic positions, relative abundance, and genomic content. We found that freshwater Nitrososphaeria were dominated by the family Nitrosopumilaceae and could be grouped into three distinct clades closely related to the genera Nitrosopumilus, Nitrosoarchaeum, and Nitrosotenuis. The Nitrosopumilus-like clade was exclusively from deep lakes, while the Nitrosoarchaeum-like clade was dominated by species from deep lakes and rivers, and the Nitrosotenuis-like clade was mainly from rivers, deep lakes, and estuaries. Interestingly, there was vertical niche separation between two clades in deep lakes, showing that the Nitrosopumilus-like species dominated shallow layers, whereas the relative abundance of the Nitrosoarchaeum-like clade increased toward deep waters. Phylogenetic clustering patterns in the Nitrosopumilaceae supported at least one freshwater-to-marine and two marine-to-freshwater transitions, the former of which refined the potential terrestrial-to-marine evolutionary path as previously proposed. The occurrence of the two marine-to-freshwater transitions were accompanied by horizontal transfer of the genes involved in nutrition regulation, osmoregulation, and cell motility during their colonization to freshwater habitats. Specifically, the Nitrosopumilus-like clade showed losses of genes encoding flagella assembly and ion transport, whereas the Nitrosoarchaeum-like clade had losses of intact genes involved in urea uptake and utilization and gains of genes encoding osmolarity-mediated mechanosensitive channels. Collectively, our results reveal for the first time the high genomic diversity of the class Nitrososphaeria across freshwater ecosystems and provide novel insights into their adaptive mechanisms and evolutionary histories.

Diverse ecophysiological adaptations of subsurface Thaumarchaeota in floodplain sediments revealed through genome-resolved metagenomics

The terrestrial subsurface microbiome contains vastly underexplored phylogenetic diversity and metabolic novelty, with critical implications for global biogeochemical cycling. Among the key microbial inhabitants of subsurface soils and sediments are Thaumarchaeota, an archaeal phylum that encompasses ammonia-oxidizing archaea (AOA) as well as non-ammonia-oxidizing basal lineages. Thaumarchaeal ecology in terrestrial systems has been extensively characterized, particularly in the case of AOA. However, there is little knowledge on the diversity and ecophysiology of Thaumarchaeota in deeper soils, as most lineages, particularly basal groups, remain uncultivated and underexplored. Here we use genome-resolved metagenomics to examine the phylogenetic and metabolic diversity of Thaumarchaeota along a 234 cm depth profile of hydrologically variable riparian floodplain sediments in the Wind River Basin near Riverton, Wyoming. Phylogenomic analysis of the metagenome-assembled genomes (MAGs) indicates a shift in AOA population structure from the dominance of the terrestrial Nitrososphaerales lineage in the well-drained top ~100 cm of the profile to the typically marine Nitrosopumilales in deeper, moister, more energy-limited sediment layers. We also describe two deeply rooting non-AOA MAGs with numerous unexpected metabolic features, including the reductive acetyl-CoA (Wood-Ljungdahl) pathway, tetrathionate respiration, a form III RuBisCO, and the potential for extracellular electron transfer. These MAGs also harbor tungsten-containing aldehyde:ferredoxin oxidoreductase, group 4f [NiFe]-hydrogenases and a canonical heme catalase, typically not found in Thaumarchaeota. Our results suggest that hydrological variables, particularly proximity to the water table, impart a strong control on the ecophysiology of Thaumarchaeota in alluvial sediments.

The adjustment of life history strategies drives the ecological adaptations of soil microbiota to aridity

Soil microbiota increase their fitness to local habitats by adjusting their life history strategies. Yet, how such adjustments drive their ecological adaptations in xeric grasslands remains elusive. In this study, shifts in the traits that potentially represent microbial life history strategies were studied along two aridity gradients with different climates using metagenomic and trait-based approaches. The results indicated that resource acquisition (e.g., higher activities of β-D-glucosidase and N-acetyl-beta-D-glucosidase, higher degradation rates of cellulose and chitin as well as genes involved in cell motilities, biodegradations, transportations and competitions) and growth yield (e.g., higher biomass and respiration) strategies were depleted with a higher aridity. However, cellular and high growth potential maintenance (e.g., higher metabolic quotients and genes related to DNA replication, transcription, translation, central carbon metabolisms and biosynthesis) and stress tolerance (e.g., genes involved in DNA damage repair, cation transporter, sporulation and osmolyte biosynthesis) strategies were enriched with a higher aridity. This implied that microbiota have lower growth yields but are probably well primed for rapid responses to pulse rainfalls in more arid soils, whereas those in less arid soils may have stronger resource acquisition and growth yield abilities. By integrating a large amount of evidence involved in taxonomic, metagenomic, genomic and biochemical properties, this study demonstrated that the ecological adaptations of soil microbiota to aridity made by adjusting and optimizing their life history strategies are universal in xeric grasslands and provided an underlying mechanistic understanding of soil microbial responses to climate changes.

Long-Term Adaptation of Acidophilic Archaeal Ammonia Oxidisers Following Different Soil Fertilisation Histories

Ammonia oxidising archaea (AOA) are ecologically important nitrifiers in acidic agricultural soils. Two AOA phylogenetic clades, belonging to order-level lineages of Nitrososphaerales (clade C11; also classified as NS-Gamma-2.3.2) and family-level lineage of Candidatus Nitrosotaleaceae (clade C14; NT-Alpha-1.1.1), usually dominate AOA population in low pH soils. This study aimed to investigate the effect of different fertilisation histories on community composition and activity of acidophilic AOA in soils. High-throughput sequencing of ammonia monooxygenase gene (amoA) was performed on six low pH agricultural plots originating from the same soil but amended with different types of fertilisers for over 20 years and nitrification rates in those soils were measured. In these fertilised acidic soils, nitrification was likely dominated by Nitrososphaerales AOA and ammonia-oxidising bacteria, while Ca. Nitrosotaleaceae AOA activity was non-significant. Within Nitrososphaerales AOA, community composition differed based on the fertilisation history, with Nitrososphaerales C11 only representing a low proportion of the community. This study revealed that long-term soil fertilisation selects for different acidophilic nitrifier communities, potentially through soil pH change or through direct effect of nitrogen, potassium and phosphorus. Comparative community composition among the differently fertilised soils also highlighted the existence of AOA phylotypes with different levels of stability to environmental changes, contributing to the understanding of high AOA diversity maintenance in terrestrial ecosystems.

Evolutionary History of RNA Modifications at N6-Adenosine Originating from the R-M System in Eukaryotes and Prokaryotes

Simple Summary

The m⁶A is the most abundant and well-studied modification of mRNA, and plays an important role in transcription and translation. It is known to be evolutionarily conserved machinery present in the last eukaryotic common ancestor (LECA). The writers and erasers responsible for adding or removing m⁶A belong to specific protein families, respectively, suggesting that these members are evolutionarily related. However, only some of these mRNA m⁶A modification-associated proteins have been studied from an evolutionary perspective, while there has been no comprehensive and systematic analysis of the distributions and evolutionary history of N6mA-associated proteins in the three kingdoms of life. In this study, we identified orthologues of all the reported N6mA-associated proteins in 88 organisms from three kingdoms of life and comprehensively reconstructed the evolutionary history of the RNA N6mA modification machinery. The results demonstrate that RNA N6mA-MTases are derived from at least two different types of prokaryotic DNA MTases (class α and β MTases). As the m⁶A reader, YTH proteins may be acquired by LECA from an individual prokaryotic YTH-domain protein that evolved from the N-terminals of an R-M system endonuclease. In addition, the origin of eukaryotic ALKBH family proteins is inferred to be driven by at least two occasions of independent HTG from the bacterial ALKB family.

Abstract

Methylation at the N6-position of adenosine (N6mA) on mRNA (m⁶A) is one of the most widespread, highly selective and dynamically regulated RNA modifications and plays an important role in transcription and translation. In the present study, a comprehensive analysis of phylogenetic relationships, conserved domain sequence characteristics and protein structure comparisons were employed to explore the distribution of RNA N6mA modification (m⁶A, m^6,6A, m⁶Am, m^{6, 6}Am and m⁶t⁶A)-associated proteins (writers, readers and erasers) in three kingdoms of life and reveal the evolutionary history of these modifications. These findings further confirmed that the restriction-modification (R-M) system is the origin of DNA and RNA N6mA modifications. Among them, the existing mRNA m⁶A modification system derived from the last eukaryotic common ancestor (LECA) is the evolutionary product of elements from the last universal common ancestor (LUCA) or driven by horizontal gene transfer (HGT) from bacterial elements. The subsequent massive gene gains and losses contribute to the development of unique and diverse functions in distinct species. Particularly, RNA methyltransferases (MTases) as the writer responsible for adding N6mA marks on mRNA and ncRNAs may have evolved from class α and β prokaryotic “orphan” MTases originating from the R-M system. The reader, YTH proteins that specifically recognize the m⁶A deposit, may be acquired by LECA from an individual prokaryotic YTH-domain protein that evolved from N-terminals of an R-M system endonuclease. The eraser, which emerged from the ALKB family (ALKBH5 and FTO) in eukaryotes, may be driven by independent HTG from bacterial ALKB proteins. The evolutionary history of RNA N6mA modifications was inferred in the present study, which will deepen our understanding of these modifications in different species.

Comparative Genomic Analyses of the Genus Nesterenkonia Unravels the Genomic Adaptation to Polar Extreme Environments

The members of the Nesterenkonia genus have been isolated from various habitats, like saline soil, salt lake, sponge-associated and the human gut, some of which are even located in polar areas. To identify their stress resistance mechanisms and draw a genomic profile across this genus, we isolated four Nesterenkonia strains from the lakes in the Tibetan Plateau, referred to as the third pole, and compared them with all other 30 high-quality Nesterenkonia genomes that are deposited in NCBI. The Heaps’ law model estimated that the pan-genome of this genus is open and the number of core, shell, cloud, and singleton genes were 993 (6.61%), 2782 (18.52%), 4117 (27.40%), and 7132 (47.47%), respectively. Phylogenomic and ANI/AAI analysis indicated that all genomes can be divided into three main clades, named NES-1, NES-2, and NES-3. The strains isolated from lakes in the Tibetan Plateau were clustered with four strains from different sources in the Antarctic and formed a subclade within NES-2, described as NES-AT. Genome features of this subclade, including GC (guanine + cytosine) content, tRNA number, carbon/nitrogen atoms per residue side chain (C/N-ARSC), and amino acid composition, in NES-AT individuals were significantly different from other strains, indicating genomic adaptation to cold, nutrient-limited, osmotic, and ultraviolet conditions in polar areas. Functional analysis revealed the enrichment of specific genes involved in bacteriorhodopsin synthesis, biofilm formation, and more diverse nutrient substance metabolism genes in the NES-AT clade, suggesting potential adaptation strategies for energy metabolism in polar environments. This study provides a comprehensive profile of the genomic features of the Nesterenkonia genus and reveals the possible mechanism for the survival of Nesterenkonia isolates in polar areas.

Progress and Challenges in Studying the Ecophysiology of Archaea

It has been less than two decades since the study of archaeal ecophysiology has become unshackled from the limitations of cultivation and amplicon sequencing through the advent of metagenomics. As a primer to the guide on producing archaeal genomes from metagenomes, we briefly summarize here how different meta'omics, imaging, and wet lab methods have contributed to progress in understanding the ecophysiology of Archaea. We then peer into the history of how our knowledge on two particularly important lineages was assembled: the anaerobic methane and alkane oxidizers, encountered primarily among Euryarchaeota, and the nanosized, mainly parasitic, members of the DPANN superphylum.

Rooting Species Trees Using Gene Tree-Species Tree Reconciliation

Interpreting phylogenetic trees requires a root, which provides the direction of evolution and polarizes ancestor-descendant relationships. But inferring the root using genetic data is difficult, particularly in cases where the closest available outgroup is only distantly related, which are common for microbes. In this chapter, we present a workflow for estimating rooted species trees and the evolutionary history of the gene families that evolve within them using probabilistic gene tree-species tree reconciliation. We illustrate the pipeline using a small dataset of prokaryotic genomes, for which the example scripts can be run using modest computer resources. We describe the rooting method used in this work in the context or other rooting strategies and discuss some of the limitations and opportunities presented by probabilistic gene tree-species tree reconciliation methods.

Reconstruction of State-Dependent Diversification: Integrating Phenotypic Traits into Molecular Phylogenies

The relative contribution of speciation and extinction into current diversity is certainly unknown, but mathematical frameworks that use genetic information have been developed to provide estimates of these processes. To that end, it is necessary to reconstruct molecular phylogenetic trees which summarize ancestor-descendant relationships as well as the timing of evolutionary events (i.e., rates). Nevertheless, diversification models show poor fit when assuming that single rate of speciation/extinction is constant over time and across lineages: species exhibit such a great variation in features that it is unlikely they give birth and die at the same pace. The state-dependent diversification framework (SSE) reconciles the species phenotypic variation with heterogeneous rates of diversification observed in a clade. This family of models allows testing contrasting hypotheses on mode of speciation, trait evolution, and its influence on speciation/extinction regimes. Although microbial species richness outnumbers diversity in plants and animals, diversification models are underused in microbiology. Here, we introduce microbiologists to models that estimate diversification rates and provide a detailed description of SSE models. Besides theoretical principles underlying the method, we also show how SSE analysis should be set up in R. We use pH evolution in Thaumarchaeota to explain its evolutionary dynamic in the light of SSE model. We hope this chapter spurs the study of trait evolution and evolutionary outcomes in microorganisms.

Leave no stone unturned: individually adapted xerotolerant Thaumarchaeota sheltered below the boulders of the Atacama Desert hyperarid core

BACKGROUND

The hyperarid core of the Atacama Desert is an extremely harsh environment thought to be colonized by only a few heterotrophic bacterial species. Current concepts for understanding this extreme ecosystem are mainly based on the diversity of these few species, yet a substantial area of the Atacama Desert hyperarid topsoil is covered by expansive boulder accumulations, whose underlying microbiomes have not been investigated so far. With the hypothesis that these sheltered soils harbor uniquely adapted microbiomes, we compared metagenomes and geochemistry between soils below and beside boulders across three distantly located boulder accumulations in the Atacama Desert hyperarid core.

RESULTS

Genome-resolved metagenomics of eleven samples revealed substantially different microbial communities in soils below and beside boulders, despite the presence of shared species. Archaea were found in significantly higher relative abundance below the boulders across all samples within distances of up to 205 km. These key taxa belong to a novel genus of ammonia-oxidizing Thaumarchaeota, Candidatus Nitrosodeserticola. We resolved eight mid-to-high quality genomes of this genus and used comparative genomics to analyze its pangenome and site-specific adaptations. Ca. Nitrosodeserticola genomes contain genes for ammonia oxidation, the 3-hydroxypropionate/4-hydroxybutyrate carbon fixation pathway, and acetate utilization indicating a chemolithoautotrophic and mixotrophic lifestyle. They also possess the capacity for tolerating extreme environmental conditions as highlighted by the presence of genes against oxidative stress and DNA damage. Site-specific adaptations of the genomes included the presence of additional genes for heavy metal transporters, multiple types of ATP synthases, and divergent genes for aquaporins.

CONCLUSION

We provide the first genomic characterization of hyperarid soil microbiomes below the boulders in the Atacama Desert, and report abundant and highly adapted Thaumarchaeaota with ammonia oxidation and carbon fixation potential. Ca. Nitrosodeserticola genomes provide the first metabolic and physiological insight into a thaumarchaeal lineage found in globally distributed terrestrial habitats characterized by various environmental stresses. We consequently expand not only the known genetic repertoire of Thaumarchaeota but also the diversity and microbiome functioning in hyperarid ecosystems. Video Abstract.

Gene Duplications Are At Least 50 Times Less Frequent than Gene Transfers in Prokaryotic Genomes

The contribution of gene duplications to the evolution of eukaryotic genomes is well studied. By contrast, studies of gene duplications in prokaryotes are scarce and generally limited to a handful of genes or careful analysis of a few prokaryotic lineages. Systematic broad-scale studies of prokaryotic genomes that sample available data are lacking, leaving gaps in our understanding of the contribution of gene duplications as a source of genetic novelty in the prokaryotic world. Here, we report conservative and robust estimates for the frequency of recent gene duplications within prokaryotic genomes relative to recent lateral gene transfer (LGT), as mechanisms to generate multiple copies of related sequences in the same genome. We obtain our estimates by focusing on evolutionarily recent events among 5,655 prokaryotic genomes, thereby avoiding vagaries of deep phylogenetic inference and confounding effects of ancient events and differential loss. We find that recent, genome-specific gene duplications are at least 50 times less frequent and probably 100 times less frequent than recent, genome-specific, gene acquisitions via LGT. The frequency of gene duplications varies across lineages and functional categories. The findings improve our understanding of genome evolution in prokaryotes and have far-reaching implications for evolutionary models that entail LGT to gene duplications ratio as a parameter.

Comparative genomic analysis reveals metabolic flexibility of Woesearchaeota

The archaeal phylum Woesearchaeota, within the DPANN superphylum, includes phylogenetically diverse microorganisms that inhabit various environments. Their biology is poorly understood due to the lack of cultured isolates. Here, we analyze datasets of Woesearchaeota 16S rRNA gene sequences and metagenome-assembled genomes to infer global distribution patterns, ecological preferences and metabolic capabilities. Phylogenomic analyses indicate that the phylum can be classified into ten subgroups, termed A-J. While a symbiotic lifestyle is predicted for most, some members of subgroup J might be host-independent. The genomes of several Woesearchaeota, including subgroup J, encode putative [FeFe] hydrogenases (known to be important for fermentation in other organisms), suggesting that these archaea might be anaerobic fermentative heterotrophs.

Patterns of Gene Content and Co-occurrence Constrain the Evolutionary Path toward Animal Association in Candidate Phyla Radiation Bacteria

Candidate Phyla Radiation (CPR) bacteria are small, likely episymbiotic organisms found across Earth's ecosystems. Despite their prevalence, the distribution of CPR lineages across habitats and the genomic signatures of transitions among these habitats remain unclear. Here, we expand the genome inventory for Absconditabacteria (SR1), Gracilibacteria, and Saccharibacteria (TM7), CPR bacteria known to occur in both animal-associated and environmental microbiomes, and investigate variation in gene content with habitat of origin. By overlaying phylogeny with habitat information, we show that bacteria from these three lineages have undergone multiple transitions from environmental habitats into animal microbiomes. Based on co-occurrence analyses of hundreds of metagenomes, we extend the prior suggestion that certain Saccharibacteria have broad bacterial host ranges and constrain possible host relationships for Absconditabacteria and Gracilibacteria. Full-proteome analyses show that animal-associated Saccharibacteria have smaller gene repertoires than their environmental counterparts and are enriched in numerous protein families, including those likely functioning in amino acid metabolism, phage defense, and detoxification of peroxide. In contrast, some freshwater Saccharibacteria encode a putative rhodopsin. For protein families exhibiting the clearest patterns of differential habitat distribution, we compared protein and species phylogenies to estimate the incidence of lateral gene transfer and genomic loss occurring over the species tree. These analyses suggest that habitat transitions were likely not accompanied by large transfer or loss events but rather were associated with continuous proteome remodeling. Thus, we speculate that CPR habitat transitions were driven largely by availability of suitable host taxa and were reinforced by acquisition and loss of some capacities. IMPORTANCE Studying the genetic differences between related microorganisms from different environment types can indicate factors associated with their movement among habitats. This is particularly interesting for bacteria from the Candidate Phyla Radiation because their minimal metabolic capabilities require associations with microbial hosts. We found that shifts of Absconditabacteria, Gracilibacteria, and Saccharibacteria between environmental ecosystems and mammalian mouths/guts probably did not involve major episodes of gene gain and loss; rather, gradual genomic change likely followed habitat migration. The results inform our understanding of how little-known microorganisms establish in the human microbiota where they may ultimately impact health.

Linking 16S rRNA Gene Classification to amoA Gene Taxonomy Reveals Environmental Distribution of Ammonia-Oxidizing Archaeal Clades in Peatland Soils

A highly resolved taxonomy for ammonia-oxidizing archaea (AOA) based on the alpha subunit of ammonia monooxygenase (amoA) was recently established, which uncovered novel environmental patterns of AOA, challenging previous generalizations. However, many microbiome studies target the 16S rRNA gene as a marker; thus, the usage of this novel taxonomy is currently limited. Here, we exploited the phylogenetic congruence of archaeal amoA and 16S rRNA genes to link 16S rRNA gene classification to the novel amoA taxonomy. We screened publicly available archaeal genomes and contigs for the co-occurring amoA and 16S rRNA genes and constructed a 16S rRNA gene database with the corresponding amoA clade taxonomy. Phylogenetic trees of both marker genes confirmed congruence, enabling the identification of clades. We validated this approach with 16S rRNA gene amplicon data from peatland soils. We succeeded in linking 16S rRNA gene amplicon sequence variants belonging to the class Nitrososphaeria to seven different AOA (amoA) clades, including two of the most frequently detected clades (Nitrososphaerales γ and δ clades) for which no pure culture is currently available. Water status significantly impacted the distribution of the AOA clades as well as the whole AOA community structure, which was correlated with pH, nitrate, and ammonium, consistent with previous clade predictions. Our study emphasizes the need to distinguish among AOA clades with distinct ecophysiologies and environmental preferences, for a better understanding of the ecology of the globally abundant AOA. IMPORTANCE The recently established phylogeny of amoA provides a finer resolution than previous studies, allowing clustering of AOA beyond the order level and thus revealing novel clades. While the 16S rRNA gene is mostly appreciated in microbiome studies, this novel phylogeny is in limited use. Here, we provide an alternative path to identifying AOA with this novel and highly resolved amoA taxonomy by using 16S rRNA gene sequencing data. We constructed a 16S rRNA gene database with the associated amoA clade taxonomy based on their phylogenetic congruence. With this database, we were able to assign 16S rRNA gene amplicons from peatland soils to different AOA clades, with a level of resolution provided previously only by amoA phylogeny. As 16S rRNA gene amplicon sequencing is still widely employed in microbiome studies, our database may have a broad application for interpreting the ecology of globally abundant AOA.

Predicting Lifestyle from Positive Selection Data and Genome Properties in Oomycetes

As evidenced in parasitism, host and niche shifts are a source of genomic and phenotypic diversification. Exemplary is a reduction in the core metabolism as parasites adapt to a particular host, while the accessory genome often maintains a high degree of diversification. However, selective pressures acting on the genome of organisms that have undergone recent lifestyle or host changes have not been fully investigated. Here, we developed a comparative genomics approach to study underlying adaptive trends in oomycetes, a eukaryotic phylum with a wide and diverse range of economically important plant and animal parasitic lifestyles. Our analysis reveals converging evolution on biological processes for oomycetes that have similar lifestyles. Moreover, we find that certain functions, in particular carbohydrate metabolism, transport, and signaling, are important for host and environmental adaptation in oomycetes. Given the high correlation between lifestyle and genome properties in our oomycete dataset, together with the known convergent evolution of fungal and oomycete genomes, we developed a model that predicts plant pathogenic lifestyles with high accuracy based on functional annotations. These insights into how selective pressures correlate with lifestyle may be crucial to better understand host/lifestyle shifts and their impact on the genome.

Protein Family Content Uncovers Lineage Relationships and Bacterial Pathway Maintenance Mechanisms in DPANN Archaea

DPANN are small-celled archaea that are generally predicted to be symbionts, and in some cases are known episymbionts of other archaea. As the monophyly of the DPANN remains uncertain, we hypothesized that proteome content could reveal relationships among DPANN lineages, constrain genetic overlap with bacteria, and illustrate how organisms with hybrid bacterial and archaeal protein sets might function. We tested this hypothesis using protein family content that was defined in part using 3,197 genomes including 569 newly reconstructed genomes. Protein family content clearly separates the final set of 390 DPANN genomes from other archaea, paralleling the separation of Candidate Phyla Radiation (CPR) bacteria from all other bacteria. This separation is partly driven by hypothetical proteins, some of which may be symbiosis-related. Pacearchaeota with the most limited predicted metabolic capacities have Form II/III and III-like Rubisco, suggesting metabolisms based on scavenged nucleotides. Intriguingly, the Pacearchaeota and Woesearchaeota with the smallest genomes also tend to encode large extracellular murein-like lytic transglycosylase domain proteins that may bind and degrade components of bacterial cell walls, indicating that some might be episymbionts of bacteria. The pathway for biosynthesis of bacterial isoprenoids is widespread in Woesearchaeota genomes and is encoded in proximity to genes involved in bacterial fatty acids synthesis. Surprisingly, in some DPANN genomes we identified a pathway for synthesis of queuosine, an unusual nucleotide in tRNAs of bacteria. Other bacterial systems are predicted to be involved in protein refolding. For example, many DPANN have the complete bacterial DnaK-DnaJ-GrpE system and many Woesearchaeota and Pacearchaeota possess bacterial group I chaperones. Thus, many DPANN appear to have mechanisms to ensure efficient protein folding of both archaeal and laterally acquired bacterial proteins.

Genome-Wide Analysis of nsLTP Gene Family and Identification of SiLTPs Contributing to High Oil Accumulation in Sesame (Sesamum indicum L.)

The biosynthesis and storage of lipids in oil crop seeds involve many gene families, such as nonspecific lipid-transfer proteins (nsLTPs). nsLTPs are cysteine-rich small basic proteins essential for plant development and survival. However, in sesame, information related to nsLTPs was limited. Thus, the objectives of this study were to identify the Sesamum indicum nsLTPs (SiLTPs) and reveal their potential role in oil accumulation in sesame seeds. Genome-wide analysis revealed 52 SiLTPs, nonrandomly distributed on 10 chromosomes in the sesame variety Zhongzhi 13. Following recent classification methods, the SiLTPs were divided into nine types, among which types I and XI were the dominants. We found that the SiLTPs could interact with several transcription factors, including APETALA2 (AP2), DNA binding with one finger (Dof), etc. Transcriptome analysis showed a tissue-specific expression of some SiLTP genes. By integrating the SiLTPs expression profiles and the weighted gene co-expression network analysis (WGCNA) results of two contrasting oil content sesame varieties, we identified SiLTPI.23 and SiLTPI.28 as the candidate genes for high oil content in sesame seeds. The presumed functions of the candidate gene were validated through overexpression of SiLTPI.23 in Arabidopsis thaliana. These findings expand our knowledge on nsLTPs in sesame and provide resources for functional studies and genetic improvement of oil content in sesame seeds.

Gene duplication drives genome expansion in a major lineage of Thaumarchaeota

Ammonia-oxidising archaea of the phylum Thaumarchaeota are important organisms in the nitrogen cycle, but the mechanisms driving their radiation into diverse ecosystems remain underexplored. Here, existing thaumarchaeotal genomes are complemented with 12 genomes belonging to the previously under-sampled Nitrososphaerales to investigate the impact of lateral gene transfer (LGT), gene duplication and loss across thaumarchaeotal evolution. We reveal a major role for gene duplication in driving genome expansion subsequent to early LGT. In particular, two large LGT events are identified into Nitrososphaerales and the fate of these gene families is highly lineage-specific, being lost in some descendant lineages, but undergoing extensive duplication in others, suggesting niche-specific roles. Notably, some genes involved in carbohydrate transport or coenzyme metabolism were duplicated, likely facilitating niche specialisation in soils and sediments. Overall, our results suggest that LGT followed by gene duplication drives Nitrososphaerales evolution, highlighting a previously under-appreciated mechanism of genome expansion in archaea.