Count: evolutionary analysis of phylogenetic profiles with parsimony and likelihood

Assembling bacterial puzzles: piecing together functions into microbial pathways

Functional metagenomics enables the study of unexplored bacterial diversity, gene families, and pathways essential to microbial communities. However, discovering biological insights with these data is impeded by the scarcity of quality annotations. Here, we use a co-occurrence-based analysis of predicted microbial protein functions to uncover pathways in genomic and metagenomic biological systems. Our approach, based on phylogenetic profiles, improves the identification of functional relationships, or participation in the same biochemical pathway, between enzymes over a comparable homology-based approach. We optimized the design of our profiles to identify potential pathways using minimal data, clustered functionally related enzyme pairs into multi-enzymatic pathways, and evaluated our predictions against reference pathways in the KEGG database. We then demonstrated a novel extension of this approach to predict inter-bacterial protein interactions amongst members of a marine microbiome. Most significantly, we show our method predicts emergent biochemical pathways between known and unknown functions. Thus, our work establishes a basis for identifying the potential functional capacities of the entire metagenome, capturing previously unknown and abstract functions into discrete putative pathways.

Host specificity and cophylogeny in the "animal-gut bacteria-phage" tripartite system

Cophylogeny has been identified between gut bacteria and their animal host and is highly relevant to host health, but little research has extended to gut bacteriophages. Here we use bee model to investigate host specificity and cophylogeny in the "animal-gut bacteria-phage" tripartite system. Through metagenomic sequencing upon different bee species, the gut phageome revealed a more variable composition than the gut bacteriome. Nevertheless, the bacteriome and the phageome showed a significant association of their dissimilarity matrices, indicating a reciprocal interaction between the two kinds of communities. Most of the gut phages were host generalist at the viral cluster level but host specialist at the viral OTU level. While the dominant gut bacteria Gilliamella and Snodgrassella exhibited matched phylogeny with bee hosts, most of their phages showed a diminished level of cophylogeny. The evolutionary rates of the bee, the gut bacteria and the gut phages showed a remarkably increasing trend, including synonymous and non-synonymous substitution and gene content variation. For all of the three codiversified tripartite members, however, their genes under positive selection and genes involving gain/loss during evolution simultaneously enriched the functions into metabolism of nutrients, therefore highlighting the tripartite coevolution that results in an enhanced ecological fitness for the whole holobiont.

Manesh MJ, Bing RG, Alexander BH, Kelly RM. Phenotype-driven assessment of the ancestral trajectory of sulfur biooxidation in the thermoacidophilic archaea Sulfolobaceae

Certain members of the family Sulfolobaceae represent the only archaea known to oxidize elemental sulfur, and their evolutionary history provides a framework to understand the development of chemolithotrophic growth by sulfur oxidation. Here, we evaluate the sulfur oxidation phenotype of Sulfolobaceae species and leverage comparative genomic and transcriptomic analysis to identify the key genes linked to sulfur oxidation. Metabolic engineering of the obligate heterotroph Sulfolobus acidocaldarius revealed that the known cytoplasmic components of sulfur oxidation alone are not sufficient to drive prolific sulfur oxidation. Imaging analysis showed that Sulfolobaceae species maintain proximity to the sulfur surface but do not necessarily contact the substrate directly. This indicates that a soluble form of sulfur must be transported to initiate cytoplasmic sulfur oxidation. Conservation patterns and transcriptomic response implicate an extracellular tetrathionate hydrolase and putative thiosulfate transporter in a newly proposed mechanism of sulfur acquisition in the Sulfolobaceae.IMPORTANCESulfur is one of the most abundant elements on earth (2.9% by mass), so it makes sense that the earliest biology found a way to use sulfur to create and sustain life. However, beyond evolutionary significance, sulfur and the molecules it comprises have important technological significance, not only in chemicals such as sulfuric acid and in pyritic ores containing critical metals but also as a waste product from oil and gas production. The thermoacidophilic Sulfolobaceae are unique among the archaea as sulfur oxidizers. The trajectory for how sulfur biooxidation arose and evolved can be traced using experimental and bioinformatic analyses of the available genomic data set. Such analysis can also inform the process by which extracellular sulfur is acquired and transported by thermoacidophilic archaea, a phenomenon that is critical to these microorganisms but has yet to be elucidated.

Dosage sensitivity shapes balanced expression and gene longevity of homoeologs after whole-genome duplications in angiosperms

Following whole-genome duplication (WGD), duplicate gene pairs (homoeologs) can evolve varying degrees of expression divergence. However, the determinants influencing these relative expression level differences (RFPKM) between homoeologs remain elusive. Here, we analyzed the RFPKM between homoeologs in three angiosperms, Nymphaea colorata, Nelumbo nucifera, and Acorus tatarinowii, all having undergone a single WGD since the origin of angiosperms. Our results show significant positive correlations in RFPKM of homoeologs among tissues within the same species, and among orthologs across these three species, indicating convergent expression balance/bias between homoeologous gene copies following independent WGDs. We linked RFPKM between homoeologs to gene attributes associated with dosage balance constraints, such as protein-protein interactions, lethal-phenotype scores in Arabidopsis (Arabidopsis thaliana) orthologs, domain numbers, and expression breadth. Notably, homoeologs with lower RFPKM often had more interactions and higher lethal-phenotype scores, indicating selective pressures favoring balanced expression. Also, homoeologs with lower RFPKM were more likely to be retained after WGDs in angiosperms. Within Nelumbo, greater RFPKM between homoeologs correlated with increased cis- and trans-regulatory differentiation between species, highlighting the ongoing escalation of gene expression divergence. We further found that expression degeneration in one copy of homoeologs is inclined towards nonfunctionalization. Our research highlights the importance of balanced expression, shaped by dosage balance constraints, in the evolutionary retention of homoeologs in plants.

A comparative genomic study of a hydrocarbon-degrading marine bacterial consortium

Ocean oil pollution has a large impact on the environment and the health of living organisms. Bioremediation cleaning strategies are promising eco-friendly alternatives for tackling this problem. Previously, we designed and reported a hydrocarbon (HC) degrading microbial consortium of four marine strains belonging to the species Alloalcanivorax xenomutans, Halopseudomonas aestusnigri, Paenarthrobacter sp., and Pseudomonas aeruginosa. However, the knowledge about the metabolic potential of this bacterial consortium for HC bioremediation is not yet well understood. Here, we analyzed the complete genomes of these marine bacterial strains accompanied by a phylogenetic reconstruction along with 138 bacterial strains. Synteny between complete genomes of the same species or genus, revealed high conservation among strains of the same species, covering over 91% of their genomic sequences. Functional predictions highlighted a high abundance of genes related to HC degradation, which may result in functional redundancy within the consortium; however, unique and complete gene clusters linked to aromatic degradation were found in the four genomes, suggesting substrate specialization. Pangenome gain and loss analysis of genes involved in HC degradation provided insights into the evolutionary history of these capabilities, shedding light on the acquisition and loss of relevant genes related to alkane and aromatic degradation. Our work, including comparative genomic analyses, identification of secondary metabolites, and prediction of HC-degrading genes, enhances our understanding of the functional diversity and ecological roles of these marine bacteria in crude oil-contaminated marine environments and contributes to the applied knowledge of bioremediation.

MAGs-centric crack: how long will, spore-positive Frankia and most Protofrankia, microsymbionts remain recalcitrant to axenic growth?

Nearly 50 years after the ground-breaking isolation of the primary Comptonia peregrina microsymbiont under axenic conditions, efforts to isolate a substantial number of Protofrankia and Frankia strains continue with enduring challenges and complexities. This study aimed to streamline genomic insights through comparative and predictive tools to extract traits crucial for isolating specific Frankia in axenic conditions. Pangenome analysis unveiled significant genetic diversity, suggesting untapped potential for cultivation strategies. Shared metabolic strategies in cellular components, central metabolic pathways, and resource acquisition traits offered promising avenues for cultivation. Ecological trait extraction indicated that most uncultured strains exhibit no apparent barriers to axenic growth. Despite ongoing challenges, potential caveats, and errors that could bias predictive analyses, this study provides a nuanced perspective. It highlights potential breakthroughs and guides refined cultivation strategies for these yet-uncultured strains. We advocate for tailored media formulations enriched with simple carbon sources in aerobic environments, with atmospheric nitrogen optionally sufficient to minimize contamination risks. Temperature adjustments should align with strain preferences-28-29°C for Frankia and 32-35°C for Protofrankia-while maintaining an alkaline pH. Given potential extended incubation periods (predicted doubling times ranging from 3.26 to 9.60 days, possibly up to 21.98 days), patience and rigorous contamination monitoring are crucial for optimizing cultivation conditions.

Cultivation of novel Atribacterota from oil well provides new insight into their diversity, ecology, and evolution in anoxic, carbon-rich environments

BACKGROUND

The Atribacterota are widely distributed in the subsurface biosphere. Recently, the first Atribacterota isolate was described and the number of Atribacterota genome sequences retrieved from environmental samples has increased significantly; however, their diversity, physiology, ecology, and evolution remain poorly understood.

RESULTS

We report the isolation of the second member of Atribacterota, Thermatribacter velox gen. nov., sp. nov., within a new family Thermatribacteraceae fam. nov., and the short-term laboratory cultivation of a member of the JS1 lineage, Phoenicimicrobium oleiphilum HX-OS.bin.34TS, both from a terrestrial oil reservoir. Physiological and metatranscriptomics analyses showed that Thermatribacter velox B11T and Phoenicimicrobium oleiphilum HX-OS.bin.34TS ferment sugars and n-alkanes, respectively, producing H2, CO2, and acetate as common products. Comparative genomics showed that all members of the Atribacterota lack a complete Wood-Ljungdahl Pathway (WLP), but that the Reductive Glycine Pathway (RGP) is widespread, indicating that the RGP, rather than WLP, is a central hub in Atribacterota metabolism. Ancestral character state reconstructions and phylogenetic analyses showed that key genes encoding the RGP (fdhA, fhs, folD, glyA, gcvT, gcvPAB, pdhD) and other central functions were gained independently in the two classes, Atribacteria (OP9) and Phoenicimicrobiia (JS1), after which they were inherited vertically; these genes included fumarate-adding enzymes (faeA; Phoenicimicrobiia only), the CODH/ACS complex (acsABCDE), and diverse hydrogenases (NiFe group 3b, 4b and FeFe group A3, C). Finally, we present genome-resolved community metabolic models showing the central roles of Atribacteria (OP9) and Phoenicimicrobiia (JS1) in acetate- and hydrocarbon-rich environments.

CONCLUSION

Our findings expand the knowledge of the diversity, physiology, ecology, and evolution of the phylum Atribacterota. This study is a starting point for promoting more incisive studies of their syntrophic biology and may guide the rational design of strategies to cultivate them in the laboratory. Video Abstract.

Convergent reductive evolution in bee-associated lactic acid bacteria

Distantly related organisms may evolve similar traits when exposed to similar environments or engaging in certain lifestyles. Several members of the Lactobacillaceae (LAB) family are frequently isolated from the floral niche, mostly from bees and flowers. In some floral LAB species (henceforth referred to as bee-associated), distinctive genomic (e.g., genome reduction) and phenotypic (e.g., preference for fructose over glucose or fructophily) features were recently documented. These features are found across distantly related species, raising the hypothesis that specific genomic and phenotypic traits evolved convergently during adaptation to the floral environment. To test this hypothesis, we examined representative genomes of 369 species of bee-associated and non-bee-associated LAB. Phylogenomic analysis unveiled seven independent ecological shifts towards the floral niche in LAB. In these bee-associated LAB, we observed pervasive, significant reductions of genome size, gene repertoire, and GC content. Using machine leaning, we could distinguish bee-associated from non-bee-associated species with 94% accuracy, based on the absence of genes involved in metabolism, osmotic stress, or DNA repair. Moreover, we found that the most important genes for the machine learning classifier were seemingly lost, independently, in multiple bee-associated lineages. One of these genes, adhE, encodes a bifunctional aldehyde-alcohol dehydrogenase associated with the evolution of fructophily, a rare phenotypic trait that was recently identified in many floral LAB species. These results suggest that the independent evolution of distinctive phenotypes in bee-associated LAB has been largely driven by independent loss of the same set of genes.

Importance

Several lactic acid bacteria (LAB) species are intimately associated with bees and exhibit unique biochemical properties with potential for food applications and honeybee health. Using a machine-learning based approach, our study shows that adaptation of LAB to the bee environment was accompanied by a distinctive genomic trajectory deeply shaped by gene loss. Several of these gene losses occurred independently in distantly related species and are linked to some of their unique biotechnologically relevant traits, such as the preference of fructose over glucose (fructophily). This study underscores the potential of machine learning in identifying fingerprints of adaptation and detecting instances of convergent evolution. Furthermore, it sheds light onto the genomic and phenotypic particularities of bee-associated bacteria, thereby deepening the understanding of their positive impact on honeybee health.

Meta-analysis of the Cetacea gut microbiome: Diversity, co-evolution, and interaction with the anthropogenic pathobiome

Despite their critical roles in marine ecosystems, only few studies have addressed the gut microbiome (GM) of cetaceans in a comprehensive way. Being long-living apex predators with a carnivorous diet but evolved from herbivorous ancestors, cetaceans are an ideal model for studying GM-host evolutionary drivers of symbiosis and represent a valuable proxy of overall marine ecosystem health. Here, we investigated the GM of eight different cetacean species, including both Odontocetes (toothed whales) and Mysticetes (baleen whales), by means of 16S rRNA-targeted amplicon sequencing. We collected faecal samples from free-ranging cetaceans circulating within the Pelagos Sanctuary (North-western Mediterranean Sea) and we also included publicly available cetacean gut microbiome sequences. Overall, we show a clear GM trajectory related to host phylogeny and taxonomy (i.e., phylosymbiosis), with remarkable GM variations which may reflect adaptations to different diets between baleen and toothed whales. While most samples were found to be infected by protozoan parasites of potential anthropic origin, we report that this phenomenon did not lead to severe GM dysbiosis. This study underlines the importance of both host phylogeny and diet in shaping the GM of cetaceans, highlighting the role of neutral processes as well as environmental factors in the establishment of this GM-host symbiosis. Furthermore, the presence of potentially human-derived protozoan parasites in faeces of free-ranging cetaceans emphasizes the importance of these animals as bioindicators of anthropic impact on marine ecosystems.

Integrated genomics provides insights into the evolution of the polyphosphate accumulation trait of Ca. Accumulibacter

Candidatus Accumulibacter, a prominent polyphosphate-accumulating organism (PAO) in wastewater treatment, plays a crucial role in enhanced biological phosphorus removal (EBPR). The genetic underpinnings of its polyphosphate accumulation capabilities, however, remain largely unknown. Here, we conducted a comprehensive genomic analysis of Ca. Accumulibacter-PAOs and their relatives within the Rhodocyclaceae family, identifying 124 core genes acquired via horizontal gene transfer (HGT) at its least common ancestor. Metatranscriptomic analysis of an enrichment culture of Ca. Accumulibacter revealed active transcription of 44 of these genes during an EBPR cycle, notably including the polyphosphate kinase 2 (PPK2) gene instead of the commonly recognized polyphosphate kinase 1 (PPK1) gene. Intriguingly, the phosphate regulon (Pho) genes showed minimal transcriptions, pointing to a distinctive fact of Pho dysregulation, where PhoU, the phosphate signaling complex protein, was not regulating the high-affinity phosphate transport (Pst) system, resulting in continuous phosphate uptake. To prevent phosphate toxicity, Ca. Accumulibacter utilized the laterally acquired PPK2 to condense phosphate into polyphosphate, resulting in the polyphosphate-accumulating feature. This study provides novel insights into the evolutionary emergence of the polyphosphate-accumulating trait in Ca. Accumulibacter, offering potential advancements in understanding the PAO phenotype in the EBPR process.

Resistance to freezing conditions of endemic Antarctic polychaetes is enhanced by cryoprotective proteins produced by their microbiome

The microbiome plays a key role in the health of all metazoans. Whether and how the microbiome favors the adaptation processes of organisms to extreme conditions, such as those of Antarctica, which are incompatible with most metazoans, is still unknown. We investigated the microbiome of three endemic and widespread species of Antarctic polychaetes: Leitoscoloplos geminus, Aphelochaeta palmeri, and Aglaophamus trissophyllus. We report here that these invertebrates contain a stable bacterial core dominated by Meiothermus and Anoxybacillus, equipped with a versatile genetic makeup and a unique portfolio of proteins useful for coping with extremely cold conditions as revealed by pangenomic and metaproteomic analyses. The close phylosymbiosis between Meiothermus and Anoxybacillus and these Antarctic polychaetes indicates a connection with their hosts that started in the past to support holobiont adaptation to the Antarctic Ocean. The wide suite of bacterial cryoprotective proteins found in Antarctic polychaetes may be useful for the development of nature-based biotechnological applications.

A non-methanogenic archaeon within the order Methanocellales

Serpentinization, a geochemical process found on modern and ancient Earth, provides an ultra-reducing environment that can support microbial methanogenesis and acetogenesis. Several groups of archaea, such as the order Methanocellales, are characterized by their ability to produce methane. Here, we generate metagenomic sequences from serpentinized springs in The Cedars, California, and construct a circularized metagenome-assembled genome of a Methanocellales archaeon, termed Met12, that lacks essential methanogenesis genes. The genome includes genes for an acetyl-CoA pathway, but lacks genes encoding methanogenesis enzymes such as methyl-coenzyme M reductase, heterodisulfide reductases and hydrogenases. In situ transcriptomic analyses reveal high expression of a multi-heme c-type cytochrome, and heterologous expression of this protein in a model bacterium demonstrates that it is capable of accepting electrons. Our results suggest that Met12, within the order Methanocellales, is not a methanogen but a CO2-reducing, electron-fueled acetogen without electron bifurcation.

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.

Pan-genome insights into adaptive evolution of bacterial symbionts in mixed host-microbe symbioses represented by human gut microbiota Bacteroides cellulosilyticus

Animal hosts harbor diverse assemblages of microbial symbionts that play crucial roles in the host's lifestyle. The link between microbial symbiosis and host development remains poorly understood. In particular, little is known about the adaptive evolution of gut bacteria in host-microbe symbioses. Recently, symbiotic relationships have been categorized as open, closed, or mixed, reflecting their modes of inter-host transmission and resulting in distinct genomic features. Members of the genus Bacteroides are the most abundant human gut microbiota and possess both probiotic and pathogenic potential, providing an excellent model for studying pan-genome evolution in symbiotic systems. Here, we determined the complete genome of an novel clinical strain PL2022, which was isolated from a blood sample and performed pan-genome analyses on a representative set of Bacteroides cellulosilyticus strains to quantify the influence of the symbiotic relationship on the evolutionary dynamics. B. cellulosilyticus exhibited correlated genomic features with both open and closed symbioses, suggesting a mixed symbiosis. An open pan-genome is characterized by abundant accessory gene families, potential horizontal gene transfer (HGT), and diverse mobile genetic elements (MGEs), indicating an innovative gene pool, mainly associated with genomic islands and plasmids. However, massive parallel gene loss, weak purifying selection, and accumulation of positively selected mutations were the main drivers of genome reduction in B. cellulosilyticus. Metagenomic read recruitment analyses showed that B. cellulosilyticus members are globally distributed and active in human gut habitats, in line with predominant vertical transmission in the human gut. However, existence and/or high abundance were also detected in non-intestinal tissues, other animal hosts, and non-host environments, indicating occasional horizontal transmission to new niches, thereby creating arenas for the acquisition of novel genes. This case study of adaptive evolution under a mixed host-microbe symbiosis advances our understanding of symbiotic pan-genome evolution. Our results highlight the complexity of genetic evolution in this unusual intestinal symbiont.

Evolution and adaptation of terrestrial plant-associated Plantibacter species into remote marine environments

Microbes are thought to be distributed and circulated around the world, but the connection between marine and terrestrial microbiomes remains largely unknown. We use Plantibacter, a representative genus associated with plants, as our research model to investigate the global distribution and adaptation of plant-related bacteria in plant-free environments, particularly in the remote Southern Ocean and the deep Atlantic Ocean. The marine isolates and their plant-associated relatives shared over 98% whole-genome average nucleotide identity (ANI), indicating recent divergence and ongoing speciation from plant-related niches to marine environments. Comparative genomics revealed that the marine strains acquired new genes via horizontal gene transfer from non-Plantibacter species and refined existing genes through positive selection to improve adaptation to new habitats. Meanwhile, marine strains retained the ability to interact with plants, such as modifying root system architecture and promoting germination. Furthermore, Plantibacter species were found to be widely distributed in marine environments, revealing an unrecognized phenomenon that plant-associated microbiomes have colonized the ocean, which could serve as a reservoir for plant growth-promoting microbes. This study demonstrates the presence of an active reservoir of terrestrial plant growth-promoting bacteria in remote marine systems and advances our understanding of the microbial connections between plant-associated and plant-free environments at the genome level.

Pan-evolutionary and regulatory genome architecture delineated by an integrated macro- and microsynteny approach

The forthcoming massive genome data generated by the Earth BioGenome Project will open up a new era of comparative genomics, for which genome synteny analysis provides an important framework. Profiling genome synteny represents an essential step in elucidating genome architecture, regulatory blocks/elements and their evolutionary history. Here we describe PanSyn, ( https://github.com/yhw320/PanSyn ), the most comprehensive and up-to-date genome synteny pipeline, providing step-by-step instructions and application examples to demonstrate its usage. PanSyn inherits both basic and advanced functions from existing popular tools, offering a user-friendly, highly customized approach for genome macrosynteny analysis and integrated pan-evolutionary and regulatory analysis of genome architecture, which are not yet available in public synteny software or tools. The advantages of PanSyn include: (i) advanced microsynteny analysis by functional profiling of microsynteny genes and associated regulatory elements; (ii) comprehensive macrosynteny analysis, including the inference of karyotype evolution from ancestors to extant species; and (iii) functional integration of microsynteny and macrosynteny for pan-evolutionary profiling of genome architecture and regulatory blocks, as well as integration with external functional genomics datasets from three- or four-dimensional genome and ENCODE projects. PanSyn requires basic knowledge of the Linux environment and Perl programming language and the ability to access a computer cluster, especially for large-scale genomic comparisons. Our protocol can be easily implemented by a competent graduate student or postdoc and takes several days to weeks to execute for dozens to hundreds of genomes. PanSyn provides yet the most comprehensive and powerful tool for integrated evolutionary and functional genomics.

Ancient eukaryotic protein interactions illuminate modern genetic traits and disorders

All eukaryotes share a common ancestor from roughly 1.5 - 1.8 billion years ago, a single-celled, swimming microbe known as LECA, the Last Eukaryotic Common Ancestor. Nearly half of the genes in modern eukaryotes were present in LECA, and many current genetic diseases and traits stem from these ancient molecular systems. To better understand these systems, we compared genes across modern organisms and identified a core set of 10,092 shared protein-coding gene families likely present in LECA, a quarter of which are uncharacterized. We then integrated >26,000 mass spectrometry proteomics analyses from 31 species to infer how these proteins interact in higher-order complexes. The resulting interactome describes the biochemical organization of LECA, revealing both known and new assemblies. We analyzed these ancient protein interactions to find new human gene-disease relationships for bone density and congenital birth defects, demonstrating the value of ancestral protein interactions for guiding functional genetics today.

The massive 340 megabase genome of Anisogramma anomala, a biotrophic ascomycete that causes eastern filbert blight of hazelnut

BACKGROUND

The ascomycete fungus Anisogramma anomala causes Eastern Filbert Blight (EFB) on hazelnut (Corylus spp.) trees. It is a minor disease on its native host, the American hazelnut (C. americana), but is highly destructive on the commercially important European hazelnut (C. avellana). In North America, EFB has historically limited commercial production of hazelnut to west of the Rocky Mountains. A. anomala is an obligately biotrophic fungus that has not been grown in continuous culture, rendering its study challenging. There is a 15-month latency before symptoms appear on infected hazelnut trees, and only a sexual reproductive stage has been observed. Here we report the sequencing, annotation, and characterization of its genome.

RESULTS

The genome of A. anomala was assembled into 108 scaffolds totaling 342,498,352 nt with a GC content of 34.46%. Scaffold N50 was 33.3 Mb and L50 was 5. Nineteen scaffolds with lengths over 1 Mb constituted 99% of the assembly. Telomere sequences were identified on both ends of two scaffolds and on one end of another 10 scaffolds. Flow cytometry estimated the genome size of A. anomala at 370 Mb. The genome exhibits two-speed evolution, with 93% of the assembly as AT-rich regions (32.9% GC) and the other 7% as GC-rich (57.1% GC). The AT-rich regions consist predominantly of repeats with low gene content, while 90% of predicted protein coding genes were identified in GC-rich regions. Copia-like retrotransposons accounted for more than half of the genome. Evidence of repeat-induced point mutation (RIP) was identified throughout the AT-rich regions, and two copies of the rid gene and one of dim-2, the key genes in the RIP mutation pathway, were identified in the genome. Consistent with its homothallic sexual reproduction cycle, both MAT1-1 and MAT1-2 idiomorphs were found. We identified a large suite of genes likely involved in pathogenicity, including 614 carbohydrate active enzymes, 762 secreted proteins and 165 effectors.

CONCLUSIONS

This study reveals the genomic structure, composition, and putative gene function of the important pathogen A. anomala. It provides insight into the molecular basis of the pathogen's life cycle and a solid foundation for studying EFB.

Phylogeny and Expansion of Serine/Threonine Kinases in Phagocytotic Bacteria in the Phylum Planctomycetota

The recently isolated bacterium "Candidatus Uabimicrobium amorphum" is the only known prokaryote that can engulf other bacterial cells. Its proteome contains a high fraction of proteins involved in signal transduction systems, which is a feature normally associated with multicellularity in eukaryotes. Here, we present a protein-based phylogeny which shows that "Ca. Uabimicrobium amorphum" represents an early diverging lineage that clusters with the Saltatorellus clade within the phylum Planctomycetota. A gene flux analysis indicated a gain of 126 protein families for signal transduction functions in "Ca. Uabimicrobium amorphum", of which 66 families contained eukaryotic-like Serine/Threonine kinases with Pkinase domains. In total, we predicted 525 functional Serine/Threonine kinases in "Ca. Uabimicrobium amorphum", which represent 8% of the proteome and is the highest fraction of Serine/Threonine kinases in a bacterial proteome. The majority of Serine/Threonine kinases in this species are membrane proteins and 30% contain long, tandem arrays of WD40 or TPR domains. The pKinase domain was predicted to be located in the cytoplasm, while the WD40 and TPR domains were predicted to be located in the periplasm. Such domain combinations were also identified in the Serine/Threonine kinases of other species in the Planctomycetota, although in much lower abundances. A phylogenetic analysis of the Serine/Threonine kinases in the Planctomycetota inferred from the Pkinase domain alone provided support for lineage-specific expansions of the Serine/Threonine kinases in "Ca. Uabimicrobium amorphum". The results imply that expansions of eukaryotic-like signal transduction systems are not restricted to multicellular organisms, but have occurred in parallel in prokaryotes with predatory lifestyles and phagocytotic-like behaviors.

Dollo Parsimony Overestimates Ancestral Gene Content Reconstructions

Ancestral reconstruction is a widely used technique that has been applied to understand the evolutionary history of gain and loss of gene families. Ancestral gene content can be reconstructed via different phylogenetic methods, but many current and previous studies employ Dollo parsimony. We hypothesize that Dollo parsimony is not appropriate for ancestral gene content reconstruction inferences based on sequence homology, as Dollo parsimony is derived from the assumption that a complex character cannot be regained. This premise does not accurately model molecular sequence evolution, in which false orthology can result from sequence convergence or lateral gene transfer. The aim of this study is to test Dollo parsimony's suitability for ancestral gene content reconstruction and to compare its inferences with a maximum likelihood-based approach that allows a gene family to be gained more than once within a tree. We first compared the performance of the two approaches on a series of artificial data sets each of 5,000 genes that were simulated according to a spectrum of evolutionary rates without gene gain or loss, so that inferred deviations from the true gene count would arise only from errors in orthology inference and ancestral reconstruction. Next, we reconstructed protein domain evolution on a phylogeny representing known eukaryotic diversity. We observed that Dollo parsimony produced numerous ancestral gene content overestimations, especially at nodes closer to the root of the tree. These observations led us to the conclusion that, confirming our hypothesis, Dollo parsimony is not an appropriate method for ancestral reconstruction studies based on sequence homology.

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.

Somatic genome architecture and molecular evolution are decoupled in "young" linage-specific gene families in ciliates

The evolution of lineage-specific gene families remains poorly studied across the eukaryotic tree of life, with most analyses focusing on the recent evolution of de novo genes in model species. Here we explore the origins of lineage-specific genes in ciliates, a ~1 billion year old clade of microeukaryotes that are defined by their division of somatic and germline functions into distinct nuclei. Previous analyses on conserved gene families have shown the effect of ciliates' unusual genome architecture on gene family evolution: extensive genome processing-the generation of thousands of gene-sized somatic chromosomes from canonical germline chromosomes-is associated with larger and more diverse gene families. To further study the relationship between ciliate genome architecture and gene family evolution, we analyzed lineage specific gene families from a set of 46 transcriptomes and 12 genomes representing x species from eight ciliate classes. We assess how the evolution lineage-specific gene families occurs among four groups of ciliates: extensive fragmenters with gene-size somatic chromosomes, non-extensive fragmenters with "large'' multi-gene somatic chromosomes, Heterotrichea with highly polyploid somatic genomes and Karyorelictea with 'paradiploid' somatic genomes. Our analyses demonstrate that: 1) most lineage-specific gene families are found at shallow taxonomic scales; 2) extensive genome processing (i.e., gene unscrambling) during development likely influences the size and number of young lineage-specific gene families; and 3) the influence of somatic genome architecture on molecular evolution is increasingly apparent in older gene families. Altogether, these data highlight the influences of genome architecture on the evolution of lineage-specific gene families in eukaryotes.

Unraveling the phylogenomic diversity of Methanomassiliicoccales and implications for mitigating ruminant methane emissions

BACKGROUND

Methanomassiliicoccales are a recently identified order of methanogens that are diverse across global environments particularly the gastrointestinal tracts of animals; however, their metabolic capacities are defined via a limited number of cultured strains.

RESULTS

Here, we profile and analyze 243 Methanomassiliicoccales genomes assembled from cultured representatives and uncultured metagenomes recovered from various biomes, including the gastrointestinal tracts of different animal species. Our analyses reveal the presence of numerous undefined genera and genetic variability in metabolic capabilities within Methanomassiliicoccales lineages, which is essential for adaptation to their ecological niches. In particular, gastrointestinal tract Methanomassiliicoccales demonstrate the presence of co-diversified members with their hosts over evolutionary timescales and likely originated in the natural environment. We highlight the presence of diverse clades of vitamin transporter BtuC proteins that distinguish Methanomassiliicoccales from other archaeal orders and likely provide a competitive advantage in efficiently handling B12. Furthermore, genome-centric metatranscriptomic analysis of ruminants with varying methane yields reveal elevated expression of select Methanomassiliicoccales genera in low methane animals and suggest that B12 exchanges could enable them to occupy ecological niches that possibly alter the direction of H2 utilization.

CONCLUSIONS

We provide a comprehensive and updated account of divergent Methanomassiliicoccales lineages, drawing from numerous uncultured genomes obtained from various habitats. We also highlight their unique metabolic capabilities involving B12, which could serve as promising targets for mitigating ruminant methane emissions by altering H2 flow.

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.

Novel evolutionary insights on the interactions of the Holosporales (Alphaproteobacteria) with eukaryotic hosts from comparative genomics

Holosporales are an alphaproteobacterial order engaging in obligate and complex associations with eukaryotes, in particular protists. The functional and evolutionary features of those interactions are still largely undisclosed. Here, we sequenced the genomes of two members of the species Bealeia paramacronuclearis (Holosporales, Holosporaceae) intracellularly associated with the ciliate protist Paramecium, which resulted in high correspondence. Consistent with the short-branched early-divergent phylogenetic position, Bealeia presents a larger functional repertoire than other Holosporaceae, comparable to those of other Holosporales families, particularly for energy metabolism and motility. Our analyses indicate that different Holosporales likely experienced at least partly autonomous genome reduction and adaptation to host interactions, for example regarding dependence on host biotin driven by multiple independent horizontal acquisitions of transporters. Among Alphaproteobacteria, this is reminiscent of the convergently evolved Rickettsiales, which however appear more diverse, possibly due to a probably more ancient origin. We identified in Bealeia and other Holosporales the plasmid-encoded putative genetic determinants of R-bodies, which may be involved in a killer trait towards symbiont-free hosts. While it is not clear whether these genes are ancestral or recently horizontally acquired, an intriguing and peculiar role of R-bodies is suggested in the evolution of the interactions of multiple Holosporales with their hosts.

Two species with a peculiar evolution within the genus Pectobacterium suggest adaptation to a new environmental niche

Historically, research on Soft Rot Pectobacteriacea (SRP) has focused on economically important crops and ornamentals and knowledge of these bacteria outside the plant context remains poorly investigated. Recently, two closely related species Pectobacterium aquaticum and Pectobacterium quasiaquaticum were isolated from water and have not been isolated from any plant yet. To identify the distinctive characteristics of these two species, we performed a comparative genomic analysis of 80 genomes representing 19 Pectobacterium species and performed an evolutionary reconstruction. Both water species underwent a reduction in genome size associated with a high pseudogene content. A high gene loss was predicted at the emergence of both species. Among the 199 gene families missing from both P. aquaticum and P. quasiaquaticum genomes but present in at least 80% of other Pectobacterium genomes, COG analysis identified many genes involved in nutrient transport systems. In addition, many type II secreted proteins were also missing in both species. Phenotypic analysis revealed that both species had reduced pectinolytic activity, a biofilm formation defect, were highly motile and had reduced virulence on several plants. These genomic and phenotypic data suggest that the ecological niche of P. aquaticum and P. quasiaquaticum may differ from that of other Pectobacterium species.

Emergence of putative energy parasites within Clostridia revealed by genome analysis of a novel endosymbiotic clade

The Clostridia is a dominant bacterial class in the guts of various animals and are considered to nutritionally contribute to the animal host. Here, we discovered clostridial endosymbionts of cellulolytic protists in termite guts, which have never been reported with evidence. We obtained (near-)complete genome sequences of three endosymbiotic Clostridia, each associated with a different parabasalid protist species with various infection rates: Trichonympha agilis, Pseudotrichonympha grassii, and Devescovina sp. All these protists are previously known to harbor permanently-associated, mutualistic Endomicrobia or Bacteroidales that supplement nitrogenous compounds. The genomes of the endosymbiotic Clostridia were small in size (1.0-1.3 Mbp) and exhibited signatures of an obligately-intracellular parasite, such as an extremely limited capability to synthesize amino acids, cofactors, and nucleotides and a disrupted glycolytic pathway with no known net ATP-generating system. Instead, the genomes encoded ATP/ADP translocase and, interestingly, regulatory proteins that are unique to eukaryotes in general and are possibly used to interfere with host cellular processes. These three genomes formed a clade with metagenome-assembled genomes (MAGs) derived from the guts of other animals, including human and ruminants, and the MAGs shared the characteristics of parasites. Gene flux analysis suggested that the acquisition of the ATP/ADP translocase gene in a common ancestor was probably key to the emergence of this parasitic clade. Taken together, we provide novel insights into the multilayered symbiotic system in the termite gut by adding the presence of parasitism and present an example of the emergence of putative energy parasites from a dominant gut bacterial clade.

Comparison of gene clustering criteria reveals intrinsic uncertainty in pangenome analyses

BACKGROUND

A key step for comparative genomics is to group open reading frames into functionally and evolutionarily meaningful gene clusters. Gene clustering is complicated by intraspecific duplications and horizontal gene transfers that are frequent in prokaryotes. In consequence, gene clustering methods must deal with a trade-off between identifying vertically transmitted representatives of multicopy gene families, which are recognizable by synteny conservation, and retrieving complete sets of species-level orthologs. We studied the implications of adopting homology, orthology, or synteny conservation as formal criteria for gene clustering by performing comparative analyses of 125 prokaryotic pangenomes.

RESULTS

Clustering criteria affect pangenome functional characterization, core genome inference, and reconstruction of ancestral gene content to different extents. Species-wise estimates of pangenome and core genome sizes change by the same factor when using different clustering criteria, allowing robust cross-species comparisons regardless of the clustering criterion. However, cross-species comparisons of genome plasticity and functional profiles are substantially affected by inconsistencies among clustering criteria. Such inconsistencies are driven not only by mobile genetic elements, but also by genes involved in defense, secondary metabolism, and other accessory functions. In some pangenome features, the variability attributed to methodological inconsistencies can even exceed the effect sizes of ecological and phylogenetic variables.

CONCLUSIONS

Choosing an appropriate criterion for gene clustering is critical to conduct unbiased pangenome analyses. We provide practical guidelines to choose the right method depending on the research goals and the quality of genome assemblies, and a benchmarking dataset to assess the robustness and reproducibility of future comparative studies.

Origins and Evolution of Novel Bacteroides in Captive Apes

Bacterial strains evolve in response to the gut environment of their hosts, with genomic changes that influence their interactions with hosts as well as with other members of the gut community. Great apes in captivity have acquired strains of Bacteroides xylanisolvens, which are common within gut microbiome of humans but not typically found other apes, thereby enabling characterization of strain evolution following colonization. Here, we isolate, sequence and reconstruct the history of gene gain and loss events in numerous captive-ape-associated strains since their divergence from their closest human-associated strains. We show that multiple captive-ape-associated B. xylanisolvens lineages have independently acquired gene complexes that encode functions related to host mucin metabolism. Our results support the finding of high genome fluidity in Bacteroides, in that several strains, in moving from humans to captive apes, have rapidly gained large genomic regions that augment metabolic properties not previously present in their relatives.

Comparative genomics of Mycoplasma feriruminatoris, a fast-growing pathogen of wild Caprinae

Mycoplasma feriruminatoris is a fast-growing Mycoplasma species isolated from wild Caprinae and first described in 2013. M. feriruminatoris isolates have been associated with arthritis, kerato conjunctivitis, pneumonia and septicemia, but were also recovered from apparently healthy animals. To better understand what defines this species, we performed a genomic survey on 14 strains collected from free-ranging or zoo-housed animals between 1987 and 2017, mostly in Europe. The average chromosome size of the M. feriruminatoris strains was 1,040±0,024 kbp, with 24 % G+C and 852±31 CDS. The core genome and pan-genome of the M. feriruminatoris species contained 628 and 1312 protein families, respectively. The M. feriruminatoris strains displayed a relatively closed pan-genome, with many features and putative virulence factors shared with species from the M. mycoides cluster, including the MIB-MIP Ig cleavage system, a repertoire of DUF285 surface proteins and a complete biosynthetic pathway for galactan. M. feriruminatoris genomes were found to be mostly syntenic, although repertoires of mobile genetic elements, including Mycoplasma Integrative and Conjugative Elements, insertion sequences, and a single plasmid varied. Phylogenetic- and gene content analyses confirmed that M. feriruminatoris was closer to the M. mycoides cluster than to the ruminant species M. yeatsii and M. putrefaciens. Ancestral genome reconstruction showed that the emergence of the M. feriruminatoris species was associated with the gain of 17 gene families, some of which encode defence enzymes and surface proteins, and the loss of 25 others, some of which are involved in sugar transport and metabolism. This comparative study suggests that the M. mycoides cluster could be extended to include M. feriruminatoris. We also find evidence that the specific organization and structure of the DnaA boxes around the oriC of M. feriruminatoris may contribute to drive the remarkable fast growth of this minimal bacterium.

Indications for a genetic basis for big bacteria and description of the giant cable bacterium Candidatus Electrothrix gigas sp. nov

Bacterial cells can vary greatly in size, from a few hundred nanometers to hundreds of micrometers in diameter. Filamentous cable bacteria also display substantial size differences, with filament diameters ranging from 0.4 to 8 µm. We analyzed the genomes of cable bacterium filaments from 11 coastal environments of which the resulting 23 new genomes represent 10 novel species-level clades of Candidatus Electrothrix and two clades that putatively represent novel genus-level diversity. Fluorescence in situ hybridization with a species-level probe showed that large-sized cable bacteria belong to a novel species with the proposed name Ca. Electrothrix gigas. Comparative genome analysis suggests genes that play a role in the construction or functioning of large cable bacteria cells: the genomes of Ca. Electrothrix gigas encode a novel actin-like protein as well as a species-specific gene cluster encoding four putative pilin proteins and a putative type II secretion platform protein, which are not present in other cable bacteria. The novel actin-like protein was also found in a number of other giant bacteria, suggesting there could be a genetic basis for large cell size. This actin-like protein (denoted big bacteria protein, Bbp) may have a function analogous to other actin proteins in cell structure or intracellular transport. We contend that Bbp may help overcome the challenges of diffusion limitation and/or morphological complexity presented by the large cells of Ca. Electrothrix gigas and other giant bacteria. IMPORTANCE In this study, we substantially expand the known diversity of marine cable bacteria and describe cable bacteria with a large diameter as a novel species with the proposed name Candidatus Electrothrix gigas. In the genomes of this species, we identified a gene that encodes a novel actin-like protein [denoted big bacteria protein (Bbp)]. The bbp gene was also found in a number of other giant bacteria, predominantly affiliated to Desulfobacterota and Gammaproteobacteria, indicating that there may be a genetic basis for large cell size. Thus far, mostly structural adaptations of giant bacteria, vacuoles, and other inclusions or organelles have been observed, which are employed to overcome nutrient diffusion limitation in their environment. In analogy to other actin proteins, Bbp could fulfill a structural role in the cell or potentially facilitate intracellular transport.

Insights into early evolutionary adaptations of the Akkermansia genus to the vertebrate gut

Akkermansia, a relevant mucin degrader from the vertebrate gut microbiota, is a member of the deeply branched Verrucomicrobiota, as well as the only known member of this phylum to be described as inhabitants of the gut. Only a few Akkermansia species have been officially described so far, although there is genomic evidence addressing the existence of more species-level variants for this genus. This niche specialization makes Akkermansia an interesting model for studying the evolution of microorganisms to their adaptation to the gastrointestinal tract environment, including which kind of functions were gained when the Akkermansia genus originated or how the evolutionary pressure functions over those genes. In order to gain more insight into Akkermansia adaptations to the gastrointestinal tract niche, we performed a phylogenomic analysis of 367 high-quality Akkermansia isolates and metagenome-assembled genomes, in addition to other members of Verrucomicrobiota. This work was focused on three aspects: the definition of Akkermansia genomic species clusters and the calculation and functional characterization of the pangenome for the most represented species; the evolutionary relationship between Akkermansia and their closest relatives from Verrucomicrobiota, defining the gene families which were gained or lost during the emergence of the last Akkermansia common ancestor (LAkkCA) and; the evaluation of the evolutionary pressure metrics for each relevant gene family of main Akkermansia species. This analysis found 25 Akkermansia genomic species clusters distributed in two main clades, divergent from their non-Akkermansia relatives. Pangenome analyses suggest that Akkermansia species have open pangenomes, and the gene gain/loss model indicates that genes associated with mucin degradation (both glycoside hydrolases and peptidases), (micro)aerobic metabolism, surface interaction, and adhesion were part of LAkkCA. Specifically, mucin degradation is a very ancestral innovation involved in the origin of Akkermansia. Horizontal gene transfer detection suggests that Akkermansia could receive genes mostly from unknown sources or from other Gram-negative gut bacteria. Evolutionary metrics suggest that Akkemansia species evolved differently, and even some conserved genes suffered different evolutionary pressures among clades. These results suggest a complex evolutionary landscape of the genus and indicate that mucin degradation could be an essential feature in Akkermansia evolution as a symbiotic species.

Acquisition of Type I methyltransferase via horizontal gene transfer increases the drug resistance of Aeromonas veronii

Aeromonas veronii is an opportunistic pathogen that affects both fish and mammals, including humans, leading to bacteraemia, sepsis, meningitis and even death. The increasing virulence and drug resistance of A. veronii are of significant concern and pose a severe risk to public safety. The Type I restriction-modification (RM) system, which functions as a bacterial defence mechanism, can influence gene expression through DNA methylation. However, little research has been conducted to explore its origin, evolutionary path, and relationship to virulence and drug resistance in A. veronii. In this study, we analysed the pan-genome of 233 A. veronii strains, and the results indicated that it was 'open', meaning that A. veronii has acquired additional genes from other species. This suggested that A. veronii had the potential to adapt and evolve rapidly, which might have contributed to its drug resistance. One Type I methyltransferase (MTase) and two complete Type I RM systems were identified, namely AveC4I, AveC4II and AveC4III in A. veronii strain C4, respectively. Notably, AveC4I was exclusive to A. veronii C4. Phylogenetic analysis revealed that AveC4I was derived from horizontal gene transfer from Thiocystis violascens and exchanged genes with the human pathogen Comamonas kerstersii. Single molecule real-time sequencing was applied to identify the motif methylated by AveC4I, which was unique and not recognized by any reported MTases in the REBASE database. We also annotated the functions and pathways of the genes containing the motif, revealing that AveC4I may control drug resistance in A. veronii C4. Our findings provide new insight on the mechanisms underlying drug resistance in pathogenic bacteria. By identifying the specific genes and pathways affected by AveC4I, this study may aid in the development of new therapeutic approaches to combat A. veronii infections.

Emergence of Raoultella ornithinolytica in human infections from different hospitals in Ecuador with OXA-48-producing resistance

Purpose

The purpose of this study was to highlight the clinical and molecular features of 13 Raoultella ornithinolytica strains isolated from clinical environments in Ecuador, and to perform comparative genomics with previously published genomes of Raoultella spp. As Raoultella is primarily found in environmental, clinical settings, we focused our work on identifying mechanisms of resistance that can provide this bacterium an advantage to establish and persist in hospital environments.

Methods

We analyzed 13 strains of Raoultella ornithinolytica isolated from patients with healthcare associated infections (HAI) in three hospitals in Quito and one in Santo Domingo de Los Tsáchilas, Ecuador, between November 2017 and April 2018. These isolates were subjected to phenotypic antimicrobial susceptibility testing, end-point polymerase chain reaction (PCR) to detect the presence of carbapenemases and whole-genome sequencing.

Results

Polymerase chain reaction revealed that seven isolates were positive isolates for blaOXA-48 and one for blaKPC-2 gene. Of the seven strains that presented the blaOXA-48 gene, six harbored it on an IncFII plasmid, one was inserted into the bacterial chromosome. The blaKPC gene was detected in an IncM2/IncR plasmid. From the bioinformatics analysis, nine genomes had the gene blaOXA-48, originating from Ecuador. Moreover, all R. ornithinolytica strains contained the ORN-1 gene, which confers resistance for β-lactams, such as penicillins and cephalosporins. Comparative genome analysis of the strains showed that the pangenome of R. ornithinolytica is considered an open pangenome, with 27.77% of core genes, which could be explained by the fact that the antibiotic resistance genes in the ancestral reconstruction are relatively new, suggesting that this genome is constantly incorporating new genes.

Conclusion

These results reveal the genome plasticity of R. ornithinolytica, particularly in acquiring antibiotic-resistance genes. The genomic surveillance and infectious control of these uncommon species are important since they may contribute to the burden of antimicrobial resistance and human health.

Reconstruction of the last bacterial common ancestor from 183 pangenomes reveals a versatile ancient core genome

BACKGROUND

Cumulative sequencing efforts have yielded enough genomes to construct pangenomes for dozens of bacterial species and elucidate intraspecies gene conservation. Given the diversity of organisms for which this is achievable, similar analyses for ancestral species are feasible through the integration of pangenomics and phylogenetics, promising deeper insights into the nature of ancient life.

RESULTS

We construct pangenomes for 183 bacterial species from 54,085 genomes and identify their core genomes using a novel statistical model to estimate genome-specific error rates and underlying gene frequencies. The core genomes are then integrated into a phylogenetic tree to reconstruct the core genome of the last bacterial common ancestor (LBCA), yielding three main results: First, the gene content of modern and ancestral core genomes are diverse at the level of individual genes but are similarly distributed by functional category and share several poorly characterized genes. Second, the LBCA core genome is distinct from any individual modern core genome but has many fundamental biological systems intact, especially those involving translation machinery and biosynthetic pathways to all major nucleotides and amino acids. Third, despite this metabolic versatility, the LBCA core genome likely requires additional non-core genes for viability, based on comparisons with the minimal organism, JCVI-Syn3A.

CONCLUSIONS

These results suggest that many cellular systems commonly conserved in modern bacteria were not just present in ancient bacteria but were nearly immutable with respect to short-term intraspecies variation. Extending this analysis to other domains of life will likely provide similar insights into more distant ancestral species.

Phylogenomic analysis of the Porphyromonas gingivalis - Porphyromonas gulae duo: approaches to the origin of periodontitis

Porphyromonas gingivalis is an oral human pathogen associated with the onset and progression of periodontitis, a chronic immune-inflammatory disease characterized by the destruction of the teeth-supporting tissue. P. gingivalis belongs to the genus Porphyromonas, which is characterized by being composed of Gram-negative, asaccharolytic, non-spore-forming, non-motile, obligatory anaerobic species, inhabiting niches such as the oral cavity, urogenital tract, gastrointestinal tract and infected wound from different mammals including humans. Among the Porphyromonas genus, P. gingivalis stands out for its specificity in colonizing the human oral cavity and its keystone pathogen role in periodontitis pathogenesis. To understand the evolutionary process behind P. gingivalis in the context of the Pophyoromonas genus, in this study, we performed a comparative genomics study with publicly available Porphyromonas genomes, focused on four main objectives: (A) to confirm the phylogenetic position of P. gingivalis in the Porphyromonas genus by phylogenomic analysis; (B) the definition and comparison of the pangenomes of P. gingivalis and its relative P. gulae; and (C) the evaluation of the gene family gain/loss events during the divergence of P. gingivalis and P. gulae; (D) the evaluation of the evolutionary pressure (represented by the calculation of Tajima-D values and dN/dS ratios) comparing gene families of P. gingivalis and P. gulae. Our analysis found 84 high-quality assemblies representing P. gingivalis and 14 P. gulae strains (from a total of 233 Porphyromonas genomes). Phylogenomic analysis confirmed that P. gingivalis and P. gulae are highly related lineages, close to P. loveana. Both organisms harbored open pangenomes, with a strong core-to-accessory ratio for housekeeping genes and a negative ratio for unknown function genes. Our analyses also characterized the gene set differentiating P. gulae from P. gingivalis, mainly associated with unknown functions. Relevant virulence factors, such as the FimA, Mfa1, and the hemagglutinins, are conserved in P. gulae, P. gingivalis, and P. loveana, suggesting that the origin of those factors occurred previous to the P. gulae - P. gingivalis divergence. These results suggest an unexpected evolutionary relationship between the P. gulae - P. gingivalis duo and P. loveana, showing more clues about the origin of the role of those organisms in periodontitis.

Evolutionary patterns of archaea predominant in acidic environment

BACKGROUND

Archaea of the order Thermoplasmatales are widely distributed in natural acidic areas and are amongst the most acidophilic prokaryotic organisms known so far. These organisms are difficult to culture, with currently only six genera validly published since the discovery of Thermoplasma acidophilum in 1970. Moreover, known great diversity of uncultured Thermoplasmatales represents microbial dark matter and underlines the necessity of efforts in cultivation and study of these archaea. Organisms from the order Thermoplasmatales affiliated with the so-called "alphabet-plasmas", and collectively dubbed "E-plasma", were the focus of this study. These archaea were found predominantly in the hyperacidic site PM4 of Parys Mountain, Wales, UK, making up to 58% of total metagenomic reads. However, these archaea escaped all cultivation attempts.

RESULTS

Their genome-based metabolism revealed its peptidolytic potential, in line with the physiology of the previously studied Thermoplasmatales isolates. Analyses of the genome and evolutionary history reconstruction have shown both the gain and loss of genes, that may have contributed to the success of the "E-plasma" in hyperacidic environment compared to their community neighbours. Notable genes among them are involved in the following molecular processes: signal transduction, stress response and glyoxylate shunt, as well as multiple copies of genes associated with various cellular functions; from energy production and conversion, replication, recombination, and repair, to cell wall/membrane/envelope biogenesis and archaella production. History events reconstruction shows that these genes, acquired by putative common ancestors, may determine the evolutionary and functional divergences of "E-plasma", which is much more developed than other representatives of the order Thermoplasmatales. In addition, the ancestral hereditary reconstruction strongly indicates the placement of Thermogymnomonas acidicola close to the root of the Thermoplasmatales.

CONCLUSIONS

This study has analysed the metagenome-assembled genome of "E-plasma", which denotes the basis of their predominance in Parys Mountain environmental microbiome, their global ubiquity, and points into the right direction of further cultivation attempts. The results suggest distinct evolutionary trajectories of organisms comprising the order Thermoplasmatales, which is important for the understanding of their evolution and lifestyle.

Eighty million years of rapid evolution of the primate Y chromosome

The Y chromosome usually plays a critical role in determining male sex and comprises sequence classes that have experienced unique evolutionary trajectories. Here we generated 19 new primate sex chromosome assemblies, analysed them with 10 existing assemblies and report rapid evolution of the Y chromosome across primates. The pseudoautosomal boundary has shifted at least six times during primate evolution, leading to the formation of a Simiiformes-specific evolutionary stratum and to the independent start of young strata in Catarrhini and Platyrrhini. Different primate lineages experienced different rates of gene loss and structural and chromatin change on their Y chromosomes. Selection on several Y-linked genes has contributed to the evolution of male developmental traits across the primates. Additionally, lineage-specific expansions of ampliconic regions have further increased the diversification of the structure and gene composition of the Y chromosome. Overall, our comprehensive analysis has broadened our knowledge of the evolution of the primate Y chromosome.

Pangenomic Study of Fusobacterium nucleatum Reveals the Distribution of Pathogenic Genes and Functional Clusters at the Subspecies and Strain Levels

Fusobacterium nucleatum is a prevalent periodontal pathogen and is associated with many systemic diseases. Our knowledge of the genomic characteristics and pathogenic effectors of different F. nucleatum strains is limited. In this study, we completed the whole genome assembly of the 4 F. nucleatum strains and carried out a comprehensive pangenomic study of 30 strains with their complete genome sequences. Phylogenetic analysis revealed that the F. nucleatum strains are mainly divided into 4 subspecies, while 1 of the sequenced strains was classified into a new subspecies. Gene composition analysis revealed that a total of 517 "core/soft-core genes" with housekeeping functions widely distributed in almost all the strains. Each subspecies had a unique gene cluster shared by strains within the subspecies. Analysis of the virulence factors revealed that many virulence factors were widely distributed across all the strains, with some present in multiple copies. Some virulence genes showed no consistent occurrence rule at the subspecies level and were specifically distributed in certain strains. The genomic islands mainly revealed strain-specific characteristics instead of subspecies level consistency, while CRISPR types and secondary metabolite biosynthetic gene clusters were identically distributed in F. nucleatum strains from the same subspecies. The variation in amino acid sites in the adhesion protein FadA did not affect the monomer and dimer 3D structures, but it may affect the binding surface and the stability of binding to host receptors. This study provides a basis for the pathogenic study of F. nucleatum at the subspecies and strain levels. IMPORTANCE We used F. nucleatum as an example to analyze the genomic characteristics of oral pathogens at the species, subspecies, and strain levels and elucidate the similarities and differences in functional genes and virulence factors among different subspecies/strains of the same oral pathogen. We believe that the unique biological characteristics of each subspecies/strain can be attributed to the differences in functional gene clusters or the presence/absence of certain virulence genes. This study showed that F. nucleatum strains from the same subspecies had similar functional gene compositions, CRISPR types, and secondary metabolite biosynthetic gene clusters, while pathogenic genes, such as virulence genes, antibiotic resistance genes, and GIs, had more strain level specificity. The findings of this study suggest that, for microbial pathogenicity studies, we should carefully consider the subspecies/strains being used, as different strains may vary greatly.

Draft genomes, phylogenomic reconstruction and comparative genome analysis of three Xenorhabdus strains isolated from soil-dwelling nematodes in Kenya

As a proven source of potent and selective antimicrobials, Xenorhabdus bacteria are important to an age plagued with difficult-to-treat microbial infections. Yet, only 27 species have been described to date. In this study, a novel Xenorhabdus species was discovered through genomic studies on three isolates from Kenyan soils. Soils in Western Kenya were surveyed for steinernematids and Steinernema isolates VH1 and BG5 were recovered from red volcanic loam soils from cultivated land in Vihiga and clay soils from riverine land in Bungoma respectively. From the two nematode isolates, Xenorhabdus sp. BG5 and Xenorhabdus sp. VH1 were isolated. The genomes of these two, plus that of X. griffiniae XN45 - this was previously isolated from Steinernema sp. scarpo that also originated from Kenyan soils - were sequenced and assembled. Nascent genome assemblies of the three isolates were of good quality with over 70 % of their proteome having known functions. These three isolates formed the X. griffiniae clade in a phylogenomic reconstruction of the genus. Their species were delineated using three overall genome relatedness indices: an unnamed species of the genus, Xenorhabdus sp. BG5, X. griffiniae VH1 and X. griffiniae XN45. A pangenome analysis of this clade revealed that over 70 % of species-specific genes encoded unknown functions. Transposases were linked to genomic islands in Xenorhabdus sp. BG5. Thus, overall genome-related indices sufficiently delineated species of two new Xenorhabdus isolates from Kenya, both of which were closely related to X. griffiniae . The functions encoded by most species-specific genes in the X. griffiniae clade remain unknown.

Obligate autotrophy at the thermodynamic limit of life in a new acetogenic bacterium

One of the important current issues of bioenergetics is the establishment of the thermodynamic limits of life. There is still no final understanding of what is the minimum value of the energy yield of a reaction that is sufficient to be used by an organism (the so-called "biological quantum of energy"). A reasonable model for determination of the minimal energy yield would be microorganisms capable of living on low-energy substrates, such as acetogenic prokaryotes. The most prominent metabolic feature of acetogens is autotrophic growth with molecular hydrogen and carbon dioxide as the substrates, which is hardly competitive in environments. Most probably, that is why only facultative autotrophic acetogens have been known so far. Here, we describe the first obligately autotrophic acetogenic bacterium Aceticella autotrophica gen. nov., sp. nov., strain 3443-3Ac^T. Phylogenetically, the new genus falls into a monophyletic group of heterotrophic bacteria of the genera Thermoanaerobacterium, Thermoanaerobacter, and Caldanaerobacter (hereinafter referred to as TTC group), where the sole acetogenic representative has so far been the facultatively autotrophic Thermoanaerobacter kivui. A. autotrophica and T. kivui both are acetogens employing energy-converting hydrogenase (Ech-acetogens) that are likely to have inherited the acetogenesis capacity vertically from common ancestor. However, their acetogenic machineries have undergone different adjustments by gene replacements due to horizontal gene transfers from different donors. Obligate autotrophy of A. autotrophica is associated with the lack of many sugar transport systems and carbohydrate catabolism enzymes that are present in other TTC group representatives, including T. kivui.

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.

Clustering pattern and evolution characteristic of microRNAs in grass carp (Ctenopharyngodon idella)

BACKGROUND

A considerable fraction of microRNAs (miRNAs) are highly conserved, and certain miRNAs correspond to genomic clusters. The clustering of miRNAs can be advantageous, possibly by allowing coordinated expression. However, little is known about the evolutionary forces responsible for the loss and acquisition of miRNA and miRNA clusters.

RESULTS

The results demonstrated that several novel miRNAs arose throughout grass carp evolution. Duplication and de novo production were critical strategies for miRNA cluster formation. Duplicates accounted for a smaller fraction of the expansion in the grass carp miRNA than de novo creation. Clustered miRNAs are more conserved and change slower, whereas unique miRNAs usually have high evolution rates and low expression levels. The expression level of miRNA expression in clusters is strongly correlated.

CONCLUSIONS

This study examines the genomic distribution, evolutionary background, and expression regulation of grass carp miRNAs. Our findings provide novel insights into the genesis and development of miRNA clusters in teleost.

Genomic characterization and molecular dating of the novel bacterium Permianibacter aggregans HW001T, which originated from Permian ground water

The Permian Basin is a unique ecosystem located in the southwest of the USA. An unanswered question is whether the bacteria in the Permian Basin adapted to the changing paleomarine environment and survived in the remnants of Permian groundwater. In our previous study, a novel bacterial strain, Permianibacter aggregans HW001T, was isolated from microalgae cultures incubated with Permian Basin waters, and was shown to originate from the Permian Ocean. In this study, strain HW001T was shown to be the representative strain of a novel family, classified as 'Permianibacteraceae'. The results of molecular dating suggested that the strain HW001T diverged ~ 447 million years ago (mya), which is the early Permian period (~ 250 mya). Genome analysis was used to access its potential energy utilization and biosynthesis capacity. A large number of transporters, carbohydrate-active enzymes and protein-degradation related genes have been annotated in the genome of strain HW001T. In addition, a series of important metabolic pathways, such as peptidoglycan biosynthesis, osmotic stress response system and multifunctional quorum sensing were annotated, which may confer the ability to adapt to various unfavorable environmental conditions. Finally, the evolutionary history of strain HW001T was reconstructed and the horizontal transfer of genes was predicted, indicating that the adaptation of P. aggregans to a changing marine environment depends on the evolution of their metabolic capabilities, especially in signal transmission. In conclusion, the results of this study provide genomic information for revealing the adaptive mechanism of strain HW001T to the changing ancient oceans.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00164-3.

Newly identified sex chromosomes in the Sphagnum (peat moss) genome alter carbon sequestration and ecosystem dynamics

Peatlands are crucial sinks for atmospheric carbon but are critically threatened due to warming climates. Sphagnum (peat moss) species are keystone members of peatland communities where they actively engineer hyperacidic conditions, which improves their competitive advantage and accelerates ecosystem-level carbon sequestration. To dissect the molecular and physiological sources of this unique biology, we generated chromosome-scale genomes of two Sphagnum species: S. divinum and S. angustifolium. Sphagnum genomes show no gene colinearity with any other reference genome to date, demonstrating that Sphagnum represents an unsampled lineage of land plant evolution. The genomes also revealed an average recombination rate an order of magnitude higher than vascular land plants and short putative U/V sex chromosomes. These newly described sex chromosomes interact with autosomal loci that significantly impact growth across diverse pH conditions. This discovery demonstrates that the ability of Sphagnum to sequester carbon in acidic peat bogs is mediated by interactions between sex, autosomes and environment.

The evolution of the gliotoxin biosynthetic gene cluster in Penicillium fungi

Fungi biosynthesize a diversity of secondary metabolites, small organic bioactive molecules that play diverse roles in fungal ecology. Fungal secondary metabolites are often encoded by physically clustered sets of genes known as biosynthetic gene clusters (BGCs). Fungi in the genus Penicillium produce diverse secondary metabolites that have been both useful (e.g., the antibiotic penicillin and the cholesterol-lowering drug mevastatin) and harmful (e.g., the mycotoxin patulin and the immunosuppressant gliotoxin) to human affairs. BGCs often also encode resistance genes that confer self-protection to the secondary metabolite-producing fungus. Some Penicillium species, such as Penicillium lilacinoechinulatum and Penicillium decumbens, are known to produce gliotoxin, a secondary metabolite with known immunosuppressant activity; however, an evolutionary characterization of the BGC responsible for gliotoxin biosynthesis among Penicillium species is lacking. Here, we examine the conservation of genes involved in gliotoxin biosynthesis and resistance in 35 Penicillium genomes from 23 species. We found homologous, less fragmented gliotoxin BGCs in 12 genomes, mostly fragmented remnants of the gliotoxin BGC in 21 genomes, whereas the remaining two Penicillium genomes lacked the gliotoxin BGC altogether. In contrast, we observed broad conservation of homologs of resistance genes that reside outside the BGC across Penicillium genomes. Evolutionary rate analysis revealed that BGCs with higher numbers of genes evolve slower than BGCs with few genes. Even though the gliotoxin BGC is fragmented to varying degrees in nearly all genomes examined, ancestral state reconstruction suggests that the ancestor of Penicillium species possessed the gliotoxin BGC. Our analyses suggest that genes that are part of BGCs can be retained in genomes long after the loss of secondary metabolite biosynthesis.

Gene gain facilitated endosymbiotic evolution of Chlamydiae

Chlamydiae is a bacterial phylum composed of obligate animal and protist endosymbionts. However, other members of the Planctomycetes-Verrucomicrobia-Chlamydiae superphylum are primarily free living. How Chlamydiae transitioned to an endosymbiotic lifestyle is still largely unresolved. Here we reconstructed Planctomycetes-Verrucomicrobia-Chlamydiae species relationships and modelled superphylum genome evolution. Gene content reconstruction from 11,996 gene families suggests a motile and facultatively anaerobic last common Chlamydiae ancestor that had already gained characteristic endosymbiont genes. Counter to expectations for genome streamlining in strict endosymbionts, we detected substantial gene gain within Chlamydiae. We found that divergence in energy metabolism and aerobiosis observed in extant lineages emerged later during chlamydial evolution. In particular, metabolic and aerobic genes characteristic of the more metabolically versatile protist-infecting chlamydiae were gained, such as respiratory chain complexes. Our results show that metabolic complexity can increase during endosymbiont evolution, adding an additional perspective for understanding symbiont evolutionary trajectories across the tree of life.

Chromosome-scale genome assembly of Glycyrrhiza uralensis revealed metabolic gene cluster centred specialized metabolites biosynthesis

A high-quality genome assembly is imperative to explore the evolutionary basis of characteristic attributes that define chemotype and provide essential resources for a molecular breeding strategy for enhanced production of medicinal metabolites. Here, using single-molecule high-fidelity (HiFi) sequencing reads, we report chromosome-scale genome assembly for Chinese licorice (Glycyrrhiza uralensis), a widely used herbal and natural medicine. The entire genome assembly was achieved in eight chromosomes, with contig and scaffold N50 as 36.02 and 60.2 Mb, respectively. With only 17 assembly gaps and half of the chromosomes having no or one assembly gap, the presented genome assembly is among the best plant genomes to date. Our results showed an advantage of using highly accurate long-read HiFi sequencing data for assembling a highly heterozygous genome including its complexed repeat content. Additionally, our analysis revealed that G. uralensis experienced a recent whole-genome duplication at approximately 59.02 million years ago post a gamma (γ) whole-genome triplication event, which contributed to its present chemotype features. The metabolic gene cluster analysis identified 355 gene clusters, which included the entire biosynthesis pathway of glycyrrhizin. The genome assembly and its annotations provide an essential resource for licorice improvement through molecular breeding and the discovery of valuable genes for engineering bioactive components and understanding the evolution of specialized metabolites biosynthesis.

The methanogen core and pangenome: conservation and variability across biology's growth temperature extremes

Temperature is a key variable in biological processes. However, a complete understanding of biological temperature adaptation is lacking, in part because of the unique constraints among different evolutionary lineages and physiological groups. Here we compared the genomes of cultivated psychrotolerant and thermotolerant methanogens, which are physiologically related and span growth temperatures from -2.5°C to 122°C. Despite being phylogenetically distributed amongst three phyla in the archaea, the genomic core of cultivated methanogens comprises about one-third of a given genome, while the genome fraction shared by any two organisms decreases with increasing phylogenetic distance between them. Increased methanogenic growth temperature is associated with reduced genome size, and thermotolerant organisms-which are distributed across the archaeal tree-have larger core genome fractions, suggesting that genome size is governed by temperature rather than phylogeny. Thermotolerant methanogens are enriched in metal and other transporters, and psychrotolerant methanogens are enriched in proteins related to structure and motility. Observed amino acid compositional differences between temperature groups include proteome charge, polarity and unfolding entropy. Our results suggest that in the methanogens, shared physiology maintains a large, conserved genomic core even across large phylogenetic distances and biology's temperature extremes.

Lactate-mediated neural plasticity genes emerged during the evolution of memory systems

The ability to record experiences and learning is present to different degrees in several species; however, the complexity and diversity of memory processes are cognitive function features that differentiate humans from other species. Lactate has recently been discovered to act as a signaling molecule for neuronal plasticity linked to long-term memory. Because lactate is not only an energy substrate for neurons but also a signaling molecule for plasticity (Magistretti and Allaman in Nat Rev Neurosci 19:235-249, 2018. https://doi.org/10.1038/nrn.2018.19 ), it is of particular interest to understand how and when memory-related genes and lactate-mediated neural plasticity (LMNP) genes emerged and evolved in humans. To understand the evolutionary origin and processes of memory and LMNP genes, we first collected information on genes related to memory and LMNP from the literature and then conducted a comparative analysis of these genes. We found that the memory and LMNP genes have different origins, suggesting that these genes may have become established gradually in evolutionarily and functional terms and not at the same time. We also found that memory and LMNP systems have a similar evolutionary history, having been formed with the gradual participation of newly emerging genes throughout their evolution. We propose that the function of LMNP as a signaling process may be evolutionarily associated with memory systems through an unidentified system that is linked by 13 common genes between memory and LMNP gene sets. This study provides evolutionary insight into the possible relationship between memory and the LMNP systems that deepens our understanding of the evolution of memory systems.