Unraveling the Evolutionary History of Marinisomatota: A Phylogenomic Perspective for Drug Discovery

Abstract

This article explores the evolutionary history of the Marinisomatota phylum through the lens of phylogenomics, addressing the needs of researchers and drug development professionals. It covers the foundational biology and taxonomic placement of these marine bacteria, details the methodological approaches for genomic and phylogenetic analysis, discusses common challenges and optimization strategies in data handling, and provides frameworks for validating findings and comparative analysis with related taxa. The synthesis offers a roadmap for leveraging evolutionary insights to identify novel biosynthetic gene clusters and therapeutic targets.

Marinisomatota 101: Phylogenomic Foundations and Evolutionary Origins

The discovery and definition of the candidate phylum Marinisomatota (also referenced in genomic databases as Marinisomatia) represents a critical node in the evolutionary history of the Bacteria domain, specifically within the expansive Candidate Phyla Radiation (CPR). A core thesis in modern phylogenomics posits that the CPR, which includes Patescibacteria, constitutes a vast, evolutionarily deep radiation of bacteria with streamlined genomes and predominantly symbiotic lifestyles. Defining Marinisomatota is not merely an exercise in cataloging diversity but a test case for hypotheses regarding genome reduction, metabolic dependency, and the origins of host association in early bacterial evolution. This guide synthesizes current taxonomic, genomic, and ecological data to define this phylum within that broader evolutionary narrative.

Core Taxonomic Characteristics

Marinisomatota are classified within the superphylum Patescibacteria (CPR). They are characterized by ultra-small cell sizes (~0.2 µm³) and significantly reduced genomes.

Table 1: Genomic and Cellular Characteristics of Marinisomatota

Characteristic	Typical Range/Value	Interpretation
Genome Size	0.8 - 1.2 Megabase pairs (Mbp)	Indicates extreme genome reduction, loss of biosynthetic pathways.
GC Content	38 - 45%	Within typical range for CPR bacteria.
16S rRNA Gene Length	~1,470 bp	Often contains conserved insertions/deletions defining the phylum.
Predicted Cell Diameter	0.2 - 0.4 µm	Filterable through 0.45 µm filters; ultramicrobacterial lifestyle.
tRNA Operon Copy Number	1 - 2	Highly limited, suggesting dependence on host translational machinery.

Metabolic & Ecological Niche

Metagenomic and single-cell genomic analyses reveal auxotrophies for most amino acids, nucleotides, and cofactors. They possess a limited respiratory chain but encode pathways for fermentation (e.g., to lactate or acetate). Crucially, they often encode type IV pilus systems and adhesin-like proteins, suggesting a host-attached lifestyle.

Primary Ecological Niche: Marinisomatota are consistently detected in anoxic, organic-rich marine sediments and subsurface aquifers. They are inferred to be episymbionts, likely attached to the surface of larger host microbes (e.g., Anaerolineae or Bacteroidota), scavenging metabolites and providing limited fermentation products in return.

Table 2: Key Metabolic Capabilities and Deficiencies

Metabolic Category	Presence/Absence	Key Genes/Pathways Identified
Glycolysis / Gluconeogenesis	Present (Partial)	gap, pgk, pgm, eno
TCA Cycle	Absent	-
Oxidative Phosphorylation	Highly Reduced	ATP synthase (atp operon) present; lacks full complexes I-IV.
Amino Acid Biosynthesis	Largely Absent	Auxotrophic for >15 amino acids.
Nucleotide Biosynthesis	Largely Absent	Limited salvage pathways only.
Lipid Biosynthesis	Present (Limited)	Partial fatty acid biosynthesis (fab genes).
Fermentation Pathways	Present	Lactate dehydrogenase (ldh), acetate kinase (ackA).

Key Experimental Protocols for Characterization

Protocol 1: Single-Cell Genome Sequencing from Environmental Samples

• Objective: Obtain whole-genome sequences of uncultivated Marinisomatota cells.
• Methodology:
  • Sample Fixation: Preserve sediment/water samples with 3% (v/v) molecular-grade glutaraldehyde (1hr, 4°C).
  • Cell Sorting: Stain with SYBR Green I, sort single ultra-small cells (<0.45 µm event trigger) via Fluorescence-Activated Cell Sorting (FACS) into 384-well plates.
  • Whole Genome Amplification (WGA): Use Multiple Displacement Amplification (MDA) with phi29 polymerase (REPLI-g Single Cell Kit, Qiagen).
  • Library Prep & Sequencing: Fragment MDA product, prepare libraries (Nextera XT), sequence on Illumina MiSeq/NextSeq (2x150 bp).
  • Genome Assembly & Binning: Assemble reads (SPAdes), bin genomes using coverage and tetranucleotide frequency (MetaBAT2). Confirm phylum-level taxonomy via CheckM and 16S rRNA phylogeny.

Protocol 2: FluorescenceIn SituHybridization (FISH) for Ecological Localization

• Objective: Visualize and confirm the episymbiotic lifestyle of Marinisomatota.
• Methodology:
  • Probe Design: Design a phylum-specific 16S rRNA-targeted oligonucleotide probe (e.g., MARINISOMA-1234) using ARB software. Label with Cy3 fluorophore.
  • Sample Fixation & Hybridization: Fix sediment slurry with 4% paraformaldehyde (3hr, 4°C). Apply probe (30% formamide, 46°C, 3hr) in hybridization buffer.
  • Washing & Imaging: Wash in pre-warmed buffer, counterstain with DAPI. Image via epifluorescence or confocal laser scanning microscopy (CLSM).
  • Analysis: Document physical association of Marinisomatota (Cy3 signal) with larger, DAPI-stained host cells.

The Scientist's Toolkit: Research Reagent Solutions

Table 3: Essential Reagents for Marinisomatota Research

Reagent/Material	Function	Example Product/Catalog #
0.1 µm & 0.45 µm Filters	Sequential filtration to size-fractionate ultra-small cells.	Polycarbonate Membrane Filters, Millipore
SYBR Green I Nucleic Acid Stain	Staining DNA for FACS detection of ultra-small cells.	Thermo Fisher Scientific, S7563
REPLI-g Single Cell Kit	Multiple Displacement Amplification (MDA) for WGA.	Qiagen, 150343
Nextera XT DNA Library Prep Kit	Preparation of sequencing libraries from low-input DNA.	Illumina, FC-131-1096
Formamide (Molecular Biology Grade)	Stringency agent in FISH hybridization buffer.	Sigma-Aldrich, F9037
Cy3-labeled Oligonucleotide Probe	Phylum-specific detection via FISH.	Custom synthesis (e.g., Eurofins Genomics)

Visualizations

Title: Workflow for Genomic & Ecological Analysis of Marinisomatota

Title: Predicted Metabolic Interactions of Marinisomatota

The advent of phylogenomics—the inference of evolutionary relationships from genome-scale data—has fundamentally reshaped our understanding of bacterial evolution. This whitepaper frames this revolution within the context of ongoing research into the evolutionary history of the candidate phylum Marinisomatota (formerly known as SAR406). This lineage, abundant in the deep oceanic waters, represents a profound evolutionary divergence within the bacterial domain. Resolving its phylogenetic placement is not merely an academic exercise; it is critical for understanding global biogeochemical cycles and for exploring a vast, untapped reservoir of novel metabolic pathways and enzymes with potential applications in biotechnology and drug discovery.

The Core Challenge: Deep Phylogenetic Resolution

Traditional phylogenetic markers, like the 16S rRNA gene, often lack sufficient signal to resolve relationships between deeply divergent phyla like Marinisomatota and other major bacterial groups. Phylogenomics overcomes this by utilizing hundreds of conserved, single-copy marker genes, providing orders of magnitude more data to distinguish between true phylogenetic signal and historical noise like horizontal gene transfer (HGT) and compositional bias.

Key Methodologies & Experimental Protocols

Phylogenomic Workflow for Deep Bacterial Phylogeny

Protocol Title: Genome-Resolved Metagenomics Coupled with Concatenated Marker Gene Phylogeny.

Detailed Methodology:

• Sample Collection & Sequencing:

  • Collect environmental samples (e.g., oceanic water column from various depths).
  • Extract high-molecular-weight genomic DNA.
  • Perform shotgun metagenomic sequencing using long-read (PacBio, Nanopore) and short-read (Illumina) technologies for hybrid assembly.

• Genome Binning & Curation:

  • Assemble reads into contigs using hybrid assemblers (e.g., metaSPAdes, Flye).
  • Bin contigs into Metagenome-Assembled Genomes (MAGs) using composition and coverage information (tools: MaxBin2, MetaBAT2).
  • Assess MAG quality (completeness, contamination) using CheckM. Select high-quality (>90% complete, <5% contaminated) MAGs representing Marinisomatota and reference taxa.

• Marker Gene Set Construction:

  • Identify a set of universal, single-copy marker genes (e.g., the 120 bacterial markers from GTDB-Tk, or the 37 genes from PhyloPhlAn).
  • Extract homologs of these markers from all MAGs and reference genomes using HMMER or similar tools.

• Multiple Sequence Alignment & Concatenation:

  • Align each marker gene family individually using MAFFT or MUSCLE.
  • Trim alignments to remove poorly aligned regions using trimAl or BMGE.
  • Concatenate all aligned marker genes into a supermatrix (phylogenomic matrix).

• Phylogenetic Inference:

  • Model Selection: Partition the supermatrix by gene or codon position. Determine the best-fit substitution model for each partition using ModelTest-NG.
  • Tree Building:
    • Maximum Likelihood (ML): Perform using IQ-TREE 2 or RAxML-NG, with branch support assessed via 1000 ultrafast bootstrap replicates.
    • Bayesian Inference (BI): Perform using MrBayes or PhyloBayes-MPI, employing site-heterogeneous models (e.g., CAT+GTR) to account for compositional bias.

• HGT and Artifact Assessment:

  • Perform individual gene tree analyses for all markers. Compare to the species tree to identify potential HGT events (using tools like ALE or GeneRax).
  • Test for the presence of systematic bias (e.g., long-branch attraction) using compositional homogeneity tests and by analyzing subsets of the data.

Workflow Visualization

Title: Phylogenomic Analysis Workflow

Comparative Genomics and Functional Profiling

Protocol Title: Pangenome and Metabolic Pathway Analysis of Marinisomatota.

• Pangenome Construction: Using a curated set of Marinisomatota MAGs, compute the pangenome using Roary or Anvi'o to define core, accessory, and unique gene sets.
• Functional Annotation: Annotate all genes against databases like KEGG, COG, and TIGRFAM using Prokka or DRAM.
• Metabolic Pathway Reconstruction: Manually reconstruct key metabolic pathways (e.g., carbon fixation, sulfur oxidation) from annotated genomes using pathway tools (MetaCyc, KEGG Mapper) and literature evidence.
• Comparative Analysis: Map the presence/absence of pathways onto the phylogenomic tree to infer ancestral metabolic states and evolutionary transitions.

Pathway Visualization

Title: Carbon Fixation via Calvin Cycle

Table 1: Impact of Phylogenomic Datasets on Phylogenetic Resolution

Phylogenetic Marker	Number of Informative Sites	Approx. Resolution Depth (Bacterial Phyla)	Support for Marinisomatota Placement (Example Study)
16S rRNA Gene	~1,400	Family/Order	Low/Conflicting (Variable across studies)
23S rRNA Gene	~2,900	Order/Class	Moderate but Inconsistent
Concatenated 31 markers	~12,000	Class/Phylum	High (Placed as a distinct class within FCB group)
Concatenated 120 markers (GTDB)	~30,000+	Phylum > Domain	Very High (Placed as a separate phylum, 'Marinisomatota')
Whole Genome Syntery	Genome-wide	Deep Divergence	Confirms unique lineage; identifies conserved genomic context

Table 2: Key Genomic & Metabolic Features of Marinisomatota from MAGs

Feature Category	Specific Finding	Prevalence in MAGs (%)	Implication for Evolution & Ecology
Genome Size	Small, Reduced (~1.5 - 2.2 Mb)	>95%	Suggensive of genome streamlining adaptation to oligotrophic ocean.
Carbon Metabolism	Presence of Form IA RubisCO (cbbL) genes	~70%	Indicates potential for dissolved inorganic carbon fixation in the dark ocean.
Sulfur Metabolism	Presence of Sox gene clusters (soxXYZAB)	~50%	Implies a role in oxidizing reduced sulfur compounds (e.g., thiosulfate).
Nitrogen Metabolism	Near universal absence of nitrification/denitrification genes	<5%	Niche differentiation from other deep-sea chemolithoautotrophs.
Respiratory Chain	High prevalence of terminal oxidases (cbb3-type, bd-type)	~100%	Adaptation to low-oxygen conditions of the mesopelagic zone.
Horizontal Gene Transfer	Evidence of HGT from Archaea (e.g., specific transporters)	Variable (~15-30% of genomes)	Complicates phylogeny but reveals adaptive evolution.

The Scientist's Toolkit: Research Reagent Solutions

Item/Category	Specific Product/Resource Example	Function in Phylogenomics Research
DNA Extraction Kit	DNeasy PowerWater Kit (Qiagen)	Efficient lysis and purification of microbial DNA from environmental seawater filters, critical for high-yield metagenomics.
Sequencing Service	Illumina NovaSeq & PacBio Sequel IIe	Provides complementary short-read (high accuracy) and long-read (scaffolding, repeat resolution) sequencing data for optimal MAG assembly.
Metagenomic Assembler	metaSPAdes (v3.15)	Specialized software for assembling complex metagenomic data from short reads into contigs.
Genome Binning Tool	MetaBAT2	Uses sequence composition and abundance across samples to cluster contigs into putative genomes (MAGs).
Quality Check Tool	CheckM2	Estimates completeness and contamination of MAGs using a machine learning model on conserved marker genes.
Phylogenomic Pipeline	GTDB-Tk (v2.3.0)	Standardized toolkit for identifying bacterial marker genes, aligning them, and inferring phylogenies consistent with the Genome Taxonomy Database.
Tree Inference Software	IQ-TREE 2 (v2.2.0)	Maximum likelihood phylogenetic software with built-in model testing and ultra-fast bootstrap, essential for large phylogenomic matrices.
Evolutionary Model	LG+F+R10 or C10 to C60 (in PhyloBayes)	Site-heterogeneous mixture models that account for variation in amino acid substitution patterns across sites, reducing systematic error in deep trees.
Functional Database	KOFAM SCAN & dbCAN2	HMM-based tools for annotating KEGG Orthologs and carbohydrate-active enzymes, enabling metabolic inference from MAGs.
Data Repository	NCBI GenBank & SRA; GTDB	Public archives for depositing MAG sequences, raw reads, and accessing standardized taxonomic classifications for phylogenetic context.

This whitepaper, framed within a broader thesis on Marinisomatota evolutionary history phylogenomics research, synthesizes current phylogenomic data to elucidate the phylum's placement within the Terrabacteria supergroup. Terrabacteria, encompassing primarily Gram-positive lineages and cyanobacteria, represents a major clade of bacteria that diversified early in the colonization of terrestrial environments. We present integrated analyses resolving Marinisomatota as a deeply branching lineage within Terrabacteria, sharing a most recent common ancestor with Cyanobacteria and Melainabacteria, supported by conserved genomic signatures and robust phylogenetic markers.

The Terrabacteria hypothesis posits that several major bacterial phyla, including Firmicutes, Actinobacteria, Chloroflexi, Cyanobacteria, and Deinococcus-Thermus, share a common ancestor that adapted to terrestrial life early in Earth's history. The recent discovery and genomic characterization of the candidate phylum Marinisomatota (previously CPR lineage) necessitates a precise phylogenetic placement to understand its ecological and evolutionary role. This analysis is critical for drug development professionals, as evolutionary relationships inform the discovery of novel biosynthetic gene clusters and unique cell wall targets.

Core Phylogenomic Analysis & Quantitative Data

Phylogenomic reconstruction was performed using a concatenated alignment of 16 ribosomal protein markers (RP16) universal to Bacteria. Bayesian inference (MrBayes) and maximum-likelihood (IQ-TREE) methods were employed on a dataset of 120 representative genomes spanning all major Terrabacteria phyla and outgroups.

Table 1: Phylogenomic Support Values for Marinisomatota Placement

Phylogenetic Clade	Bayesian Posterior Probability	ML UltraFast Bootstrap (%)	Approximate Likelihood Ratio Test (%)
Terrabacteria (total group)	1.00	100	100
Marinisomatota + (Cyanobacteria + Melainabacteria)	0.98	96	95
Cyanobacteria + Melainabacteria	1.00	100	100
Firmicutes + Actinobacteria	1.00	100	100

Table 2: Conserved Molecular Synapomorphies in Terrabacteria Lineages

Genomic Feature	Marinisomatota	Cyanobacteria	Firmicutes	Actinobacteria	Outgroup (Pseudomonadota)
RP16 Gene Cluster Order	Conserved block A	Conserved block A	Conserved block B	Conserved block B	Variable
PE/PPE Protein Domain	Absent	Absent	Present (some)	Present	Absent
S-layer Gene (slp)	Present (divergent)	Absent	Present	Present	Absent
Cobalamin Synthesis Pathway	Reduced	Complete	Variable	Complete	Variable

Detailed Experimental Protocols

Protocol: Genome-Resolved Metagenomics forMarinisomatotaRecovery

• Sample Collection & DNA Extraction: Collect environmental samples (marine sediment, aquifer). Use the DNeasy PowerSoil Pro Kit (Qiagen) with bead-beating for 10 min at 30 Hz to lyse cells.
• Metagenomic Sequencing: Construct libraries with Nextera XT DNA Library Prep Kit. Sequence on Illumina NovaSeq (2x150 bp) and PacBio HiFi (15 kb insert) platforms for hybrid assembly.
• Assembly & Binning: Assemble reads using metaSPAdes (v3.15.0). Recover genomes via differential coverage binning in Anvi'o (v7) using CONCOCT and Metabat2. Check for completeness/contamination with CheckM2.
• Phylogenomic Matrix Construction: Identify RP16 genes with fetchMG. Align each protein with MAFFT-linsi. Trim alignments with TrimAl (-automated1). Concatenate alignments using PhyloPhlAn.

Protocol: Phylogenetic Tree Inference & Validation

• Model Testing & Tree Search: Determine best-fit model (LG+C60+F+G) using ModelFinder in IQ-TREE2. Run maximum-likelihood analysis with 1000 UFBoot replicates.
• Bayesian Inference: Run MrBayes (v3.2.7) for 1,000,000 generations, sampling every 1000. Assess convergence (average standard deviation of split frequencies <0.01).
• Topology Testing: Use IQ-TREE's KH and SH tests to compare the optimal tree against alternative placements of Marinisomatota.

Visualization of Phylogenetic Relationships & Workflow

Phylogenomic Placement of Marinisomatota

Genome-Resolved Metagenomics Workflow

The Scientist's Toolkit: Key Research Reagent Solutions

Table 3: Essential Reagents and Materials for Phylogenomic Analysis of Marinisomatota

Item (Supplier - Catalog #)	Function in Protocol	Critical Parameters
DNeasy PowerSoil Pro Kit (Qiagen - 47014)	High-yield, inhibitor-free DNA extraction from complex environmental matrices.	Bead-beating time is critical for lysing recalcitrant Marinisomatota cells.
Nextera XT DNA Library Prep Kit (Illumina - FC-131-1096)	Prepares sequencing libraries from low-input genomic DNA for Illumina platforms.	Optimal for fragmented metagenomic DNA; normalization is key for even coverage.
SMRTbell Prep Kit 3.0 (PacBio - 102-092-000)	Prepares high molecular weight DNA for PacBio HiFi sequencing.	Essential for obtaining long reads to span repetitive regions in assembly.
Phusion High-Fidelity DNA Polymerase (NEB - M0530L)	PCR amplification of phylogenetic marker genes from genomic DNA.	High fidelity reduces errors in downstream sequence alignment.
IQ-TREE2 Software (http://www.iqtree.org)	Performs maximum-likelihood phylogenetic inference with model testing.	Use -m MFP flag for automatic model selection; -B 1000 for bootstrap.
CheckM2 Database (https://github.com/chklovski/CheckM2)	Assesses completeness and contamination of recovered MAGs.	Uses machine-learning models trained on diverse bacterial lineages, ideal for novel phyla.

This technical guide details methodologies for identifying core genomic signatures within the context of Marinisomatota phylogenomics. We present a computational and experimental framework for elucidating conserved genes and pathways critical to understanding the evolutionary history and metabolic adaptation of this candidate phylum, with direct implications for novel enzyme and drug target discovery.

The candidate phylum Marinisomatota (formerly SAR406) represents a deep-branching, globally distributed lineage of marine bacteria. Its evolutionary history, characterized by genome reduction and niche adaptation in oxygen minimum zones, makes it a prime subject for core genome analysis. Identifying conserved genomic signatures within this phylum is essential for reconstructing its metabolic evolution and identifying stable functional elements with biotechnological and therapeutic potential.

Defining Core Genomic Signatures

A core genomic signature refers to the set of genes, regulatory elements, and pathways universally present across all representative genomes of a monophyletic group, under a defined threshold (e.g., ≥95% prevalence). For Marinisomatota, this signature reveals the minimal genetic toolkit for survival in pelagic marine environments.

Quantitative Core Genome Analysis of Marinisomatota

Recent phylogenomic studies analyzing publicly available metagenome-assembled genomes (MAGs) provide the following statistics.

Table 1: Core Genome Metrics for Marinisomatota (Representative Analysis)

Metric	Value	Analysis Parameters
Number of Analyzed MAGs	112	Quality: ≥50% completeness, ≤5% contamination
Total Pan-Genome Size	~52,000 gene clusters	Protein clustering at 50% AA identity
Core Genome Size (95%)	152 genes	Present in ≥107 of 112 genomes
Shell Genome	~4,200 gene clusters	Present in 15% to 95% of genomes
Cloud Genome	~47,600 gene clusters	Present in <15% of genomes
Estimated Core Genome %	~0.3% of pan-genome	Reflects high genetic diversity

Methodologies for Identification

Computational Pipeline for Core Gene Identification

Protocol 1: Genome Curation and Core Gene Callin*

• Data Acquisition: Retrieve all high-quality Marinisomatota MAGs from public repositories (NCBI, IMG/M, GTDB).
• Quality Filtering: Retain genomes with ≥50% completeness (CheckM2) and ≤5% contamination.
• Gene Prediction & Annotation: Use Prodigal for ORF calling. Annotate via eggNOG-mapper v5.0 against COG/KEGG databases.
• Protein Clustering: Perform all-vs-all alignment (DIAMOND). Cluster proteins into homologous groups using MMseqs2 (Linclust) with parameters: --cov-mode 1 -c 0.8 --kmer-per-seq 100.
• Core Definition: Calculate presence/absence matrix. Define core gene clusters as those present in ≥95% of genomes.
• Functional Enrichment: Perform statistical overrepresentation analysis (Fisher's exact test) of KEGG pathways in the core set versus the accessory genome.

Experimental Validation of Core Pathways

Protocol 2: Heterologous Expression and Enzyme Assay for Conserved Glycolysis This protocol validates the function of a core metabolic pathway gene.

• Target: Conserved glyceraldehyde-3-phosphate dehydrogenase (gapA gene) identified in 110/112 MAGs.
• Cloning: Amplify gapA homolog from Marinisomatota-enriched metagenomic DNA using degenerate primers. Ligate into pET-28a(+) expression vector with N-terminal His-tag.
• Expression: Transform E. coli BL21(DE3). Induce expression with 0.5 mM IPTG at 16°C for 18 hours.
• Purification: Lyse cells via sonication. Purify protein using Ni-NTA affinity chromatography. Confirm purity via SDS-PAGE.
• Activity Assay: Monitor NADH production at 340 nm in reaction mixture: 50 mM Tris-HCl (pH 8.5), 5 mM D-glyceraldehyde-3-phosphate, 1 mM NAD+, 10 mM arsenate, 2 µg purified enzyme. Calculate specific activity (µmol NADH min⁻¹ mg⁻¹).

Table 2: Key Reagent Solutions for Protocol 2

Reagent / Material	Function / Rationale
pET-28a(+) Vector	T7 expression vector providing high-level, inducible expression and N-terminal His-tag for purification.
E. coli BL21(DE3)	Expression host deficient in Lon and OmpT proteases, containing T7 RNA polymerase gene for inducible expression.
Ni-NTA Agarose Resin	Immobilized metal affinity chromatography (IMAC) resin that selectively binds polyhistidine-tagged recombinant proteins.
D-Glyceraldehyde-3-Phosphate (G3P)	Substrate for the GAPDH enzyme assay. Unstable; must be prepared fresh from diethyl acetal monobarium salt.
NAD+ Coenzyme	Oxidized nicotinamide adenine dinucleotide; electron acceptor in the GAPDH reaction, reduction to NADH is measured spectrophotometrically.

Conserved Pathways in Marinisomatota Evolution

Core analysis reveals retention of essential energy and information processing pathways, alongside loss of biosynthetic capabilities, consistent with an oligotrophic lifestyle.

Table 3: Conserved Core Pathways in Marinisomatota

Pathway (KEGG Map)	Core Genes Identified	Prevalence (%)	Inferred Evolutionary/Functional Significance
Glycolysis / Gluconeogenesis	gapA, pgk, gpmI, eno, pyk	98-100	Core energy conservation; possible gluconeogenic carbon assimilation.
TCA Cycle (Incomplete)	acnB, icd, sucD, sucC, sdhA, sdhB, fumC, mdh	95-100	Split or incomplete cycle for precursor biosynthesis, not energy generation.
Ribosome Biogenesis	Multiple rps, rpl, inf genes	100	Universal protein synthesis machinery.
DNA Replication	dnaA, dnaN, gyrA, gyrB, polA	100	Essential information processing.
ABC Transporters	Subunits for branched-chain AA, Zn²⁺, phosphate	96-99	Scavenging of nutrients (amino acids, ions) from environment.

Applications in Drug Development

Core essential genes represent promising targets for novel antimicrobials against pathogenic relatives. For example, the uniquely conserved DnaN (sliding clamp) protein in Marinisomatota and its sister phyla may have distinct structural features exploitable for narrow-spectrum antibiotic design.

Protocol 3: In Silico Drug Target Prioritization Pipeline

• Target List: Generate from core gene list (Table 3). Filter for genes absent in human gut microbiome (NCBI dataset) and human genome (BLASTp e-value < 1e-10).
• Essentiality Assessment: Perform homology mapping to essential genes in model bacteria (e.g., E. coli Keio collection).
• Druggability Prediction: Submit protein sequences to DrugBank database or use machine learning tools (e.g., DeepDrug) to predict binding pocket characteristics.
• Conservation Analysis: Generate multiple sequence alignment of target across Marinisomatota and related phyla. Identify absolutely conserved residues for targeted inhibition.

The identification of core genomic signatures within Marinisomatota provides a powerful lens into the evolutionary forces shaping this enigmatic phylum. The conserved core of ~152 genes underscores a minimal, efficient genome streamlined for survival in the marine water column. The experimental and computational frameworks outlined here offer a template for similar analyses in other microbial candidate phyla, bridging phylogenomics and applied drug discovery.

This whitepaper situates the ecological drivers of marine adaptation within the emerging framework of Marinisomatota evolutionary history research. Marinisomatota (proposed candidate phylum within the FCB group) represents a phylogenetically distinct bacterial lineage with significant adaptations to pelagic and benthic marine niches. Phylogenomic analyses reveal that evolutionary trajectories within this group are fundamentally sculpted by specific abiotic and biotic pressures of marine ecosystems, including hydrostatic pressure, salinity gradients, oligotrophy, and unique chemical symbioses. Understanding these drivers is critical for elucidating the evolutionary history of the domain Bacteria and for exploiting marine-adapted biochemistry in pharmaceutical development.

Key Ecological Drivers and Genomic Adaptations

Marine environments impose distinct selective pressures. The following adaptations, inferred from comparative genomics and experimental studies of Marinisomatota and related marine microbes, are central to evolutionary success.

Table 1: Core Ecological Drivers and Corresponding Genomic/Physiological Adaptations

Ecological Driver	Selective Pressure	Evolutionary Adaptation (Marinisomatota Hallmarks)	Key Genomic Evidence
High Salinity & Osmolarity	Cellular dehydration, ion toxicity.	Synthesis of compatible solutes (e.g., glycine betaine, ectoine); Ion transport regulation.	Prevalence of bet, proU, and ect gene clusters in metagenome-assembled genomes (MAGs).
High Hydrostatic Pressure (Abyssal zones)	Protein denaturation, membrane compression.	Increased unsaturated fatty acid synthesis; Chaperone protein systems (e.g., GroEL/GroES).	Enrichment of desaturase genes and pressure-regulated operons in piezophile MAGs.
Oligotrophy (Low Nutrients)	Energy and carbon limitation.	High-affinity substrate transporters (ABC transporters); Genome streamlining; Auxotrophy compensated by symbiosis.	Reduced genome size; High % of transporter genes; CRISPR-Cas systems for viral defense.
Low Temperature (Deep sea, polar)	Reduced enzyme kinetics, membrane rigidity.	Production of antifreeze proteins (AFPs); Cold-shock proteins (Csps); Modulated lipid desaturation.	Identification of putative afp and csp homologs in polar Marinisomatota MAGs.
Specialized Symbioses (e.g., with marine sponges)	Need for niche colonization, metabolite exchange.	Loss of redundant metabolic pathways; Acquisition of symbiosis factors (adhesins, T3SS).	Genome reduction and presence of t3ss gene clusters in host-associated lineages.

Experimental Protocols for Key Investigations

Protocol: Cultivation and Pressure Simulation for Piezophile Isolation

Objective: Isolate and characterize pressure-adapted Marinisomatota from deep-sea sediments. Materials: High-pressure bioreactor (e.g., Pernod-type vessel), anaerobic chamber, marine agar 2216, sediment cores from hydrothermal vent. Procedure:

• Sample Collection: Collect sediment cores using a Niskin bottle or box corer from a depth >2000m. Maintain at in situ temperature (4°C).
• Enrichment: Inoculate 1g of sediment into anaerobic, pressurized bioreactor containing marine broth, pre-reduced with cysteine. Set initial pressure to 20 MPa.
• Serial Transfer: Incubate at 4°C for 4 weeks. Transfer 10% culture volume to fresh medium every 2 weeks, gradually increasing pressure to target levels (up to 50 MPa).
• Isolation: Plate enrichment culture onto solid marine media inside anaerobic chamber. Incubate plates in pressurized, temperature-controlled chambers.
• Characterization: Perform 16S rRNA gene sequencing and whole-genome sequencing of isolates. Analyze fatty acid methyl esters (FAME) for membrane composition.

Protocol: Phylogenomic Analysis of Adaptation Genes

Objective: Identify horizontally acquired genes and positively selected sites in Marinisomatota MAGs. Materials: High-performance computing cluster, bioinformatics software (OrthoFinder, IQ-TREE, HyPhy). Procedure:

• Data Collection: Download all available Marinisomatota MAGs from public repositories (e.g., NCBI, IMG/M).
• Gene Family Identification: Use OrthoFinder with DIAMOND for all-vs-all protein sequence comparison to define orthologous groups (OGs).
• Phylogeny Reconstruction: Concatenate single-copy core genes. Build maximum-likelihood tree with IQ-TREE (model TEST).
• Selection Analysis: For OGs of interest (e.g., ion transporters), perform codon alignment. Use the BUSTED method in HyPhy to test for gene-wide episodic diversifying selection.
• Ancestral State Reconstruction: Reconstruct presence/absence of key adaptive genes (e.g., ectoine synthase) across nodes to infer timing of acquisition.

Visualizations

Title: Marine Driver to Adaptation Logic Flow

Title: Piezophile Isolation Workflow

Title: Environmental Stress Signal Transduction

The Scientist's Toolkit: Research Reagent Solutions

Table 2: Essential Reagents and Materials for Marine Evolutionary Genomics

Item Name	Supplier Examples	Function in Research
Marine Broth 2216	BD Difco, Sigma-Aldrich	Standardized complex medium for cultivation of heterotrophic marine bacteria.
Pernod-Type High-Pressure Bioreactor	Kobe Steel, custom fabricators	Maintains in situ hydrostatic pressures (up to 100 MPa) for cultivating piezophiles.
Anaerobic Chamber (Coy Type)	Coy Lab Products, Baker	Creates oxygen-free atmosphere for cultivating anaerobic Marinisomatota.
Cryoprotectant for Marine Microbes (e.g., DMSO, Glycerol in Marine Salts)	Sigma-Aldrich, Thermo Fisher	Long-term preservation of marine isolates at -80°C or in liquid nitrogen.
Metagenomic DNA Extraction Kit (for Marine Sediments)	Qiagen PowerSoil, MoBio	Efficient lysis and purification of inhibitor-free DNA from complex marine samples.
Long-Read Sequencing Chemistry (PacBio HiFi, Oxford Nanopore)	Pacific Biosciences, Oxford Nanopore	Generates complete, closed genomes and MAGs from complex communities.
Phylogenomic Analysis Pipeline Software (OrthoFinder, IQ-TREE, HyPhy)	Open Source (GitHub)	For identifying orthologs, reconstructing phylogenies, and detecting selection.
Fluorescent In Situ Hybridization (FISH) Probes (specific for Marinisomatota 16S rRNA)	Biomers, custom synthesis	Visualizes and quantifies uncultured Marinisomatota cells in environmental samples or host tissue.

From Genomes to Trees: Methodologies for Marinisomatota Phylogenomics

Understanding the evolutionary history of the phylum Marinisomatota (formerly SAR406) is a significant challenge in microbial oceanography and evolution. This deep-branching, largely uncultivated lineage is abundant in the oceanic dark matter. Phylogenomics research into its adaptation, diversification, and metabolic roles hinges on obtaining high-quality genomic data. Two primary strategies are employed: sequencing cultured isolates and reconstructing Metagenome-Assembled Genomes (MAGs). This guide details the technical merits, protocols, and applications of each approach within this specific research context.

Core Comparison: Cultured Isolates vs. MAGs

Table 1: Quantitative and Qualitative Comparison of Sequencing Strategies

Parameter	Cultured Isolate Genomics	Metagenome-Assembled Genomes (MAGs)
Genome Completeness	Typically 100%; closed circular chromosomes possible.	Variable; commonly 70-95% for medium-high quality.
Contamination Level	Negligible (pure culture).	Measured by checkM; <5% for high-quality MAGs.
Strain Heterogeneity	Clonal, homogeneous population.	May represent consensus of closely related strains.
Technical Replicates	Straightforward from same culture.	Challenging; depends on sample availability & reprocessing.
Primary Cost Driver	Cultivation efforts, medium optimization, single-genome sequencing.	Deep sequencing depth, high-performance computing, binning.
Time to Genome	Months to years (cultivation) + weeks (sequencing/assembly).	Weeks (sequencing/binning) + weeks to months (curation).
Metabolic Context	Provides potential, not always expressed in situ.	Reflects in situ functional potential of dominant population.
Gold Standard for	Type material, reference genomes, physiological experiments.	Capturing uncultivable diversity, in situ population genomics.
Key Tool Examples	PLATEN, HGAP, Flye (for assembly).	MEGAHIT, metaSPAdes, MaxBin, MetaBAT, checkM, GTDB-Tk.

Experimental Protocols

Protocol for Cultured Isolate Genome Sequencing (Marinisomatota Focus)

Aim: Generate a complete, closed reference genome from a Marinisomatota isolate. Workflow:

• Cultivation: Employ dilution-to-extinction or high-throughput cultivation techniques using amended seawater media under in situ-like conditions (e.g., low nutrient, dark/oxygen gradients).
• Purity Verification: Check via 16S rRNA gene sequencing and microscopy (DAPI, FISH).
• High-Molecular-Weight DNA Extraction: Use a gentle lysis method (e.g., enzymatic lysis followed by CTAB/phenol-chloroform) to obtain >40 kb DNA. Assess quality via pulse-field gel electrophoresis or FEMTO Pulse.
• Library Preparation & Sequencing:
  • Long-Read Sequencing (PacBio HiFi or Oxford Nanopore): Prepare SMRTbell or ligation sequencing library. Sequence to achieve >100x coverage.
  • Optional Short-Read Polishing: Prepare an Illumina PCR-free library (350-550 bp insert). Sequence to achieve >50x coverage.
• Genome Assembly & Curation:
  • Assemble long reads using Flye or hifiasm.
  • Polish the assembly with long reads (Medaka) and optionally with Illumina reads (Pilon).
  • Check circularity and overlap termini. Annotate using the DDBJ/ENA/NCBI pipeline or Prokka.

Protocol for MAG Generation from Marine Metagenomes

Aim: Reconstruct high-quality Marinisomatota MAGs from complex marine metagenomic datasets. Workflow:

• Sample Collection & DNA Extraction: Filter large volumes of seawater (0.1-0.8 µm pore size). Use a direct lysis kit (e.g., DNeasy PowerWater) to capture community DNA, including from cells with delicate membranes.
• Metagenomic Library Preparation & Sequencing: Prepare Illumina paired-end libraries (typically 2x150 bp). Sequence deeply (>50 Gbp per sample) to ensure sufficient coverage for low-abundance taxa.
• Quality Control & Assembly: Trim adapters and low-quality bases with Trimmomatic or fastp. Perform de novo co-assembly of multiple samples or assemble individually using MEGAHIT or metaSPAdes.
• Binning: Map quality-filtered reads back to contigs (>1.5 kbp) to generate coverage profiles. Execute binning using an ensemble approach (e.g., MetaBAT2, MaxBin2, CONCOCT). Aggregate results with DAS Tool.
• MAG Curation & Taxonomy:
  • Assess bin quality with checkM2 for completeness and contamination.
  • Assign taxonomy using GTDB-Tk against the Genome Taxonomy Database.
  • Refine Marinisomatota MAGs via manual curation in Anvi'o (e.g., removal of contaminant contigs based on differential coverage, tRNA presence, GC content).

MAG Generation and Analysis Workflow

The Scientist's Toolkit: Key Research Reagents & Materials

Table 2: Essential Reagents and Tools for Marinisomatota Genome Studies

Item	Function / Role	Example Product / Tool
0.1-0.8 µm Filters	Size-fractionation to capture microbial cells, including ultramicrobacteria.	Polycarbonate track-etched or Supor membrane filters.
Direct Lysis DNA Kit	Efficiently lyse diverse, hard-to-lyse microbial cells (e.g., Marinisomatota) in environmental samples.	DNeasy PowerWater Kit, FastDNA Spin Kit for Soil.
PacBio SMRTbell Kit	Preparation of high-fidelity (HiFi) long-read sequencing libraries from isolate DNA.	SMRTbell Express Template Prep Kit 3.0.
Illumina PCR-free Kit	Preparation of shotgun metagenomic or isolate libraries without amplification bias.	Nextera DNA Flex Library Prep (PCR-free protocol).
checkM2	Assess completeness and contamination of MAGs using machine learning models.	Open-source software (github.com/chklovski/checkM2).
GTDB-Tk	Assign standardized taxonomic labels to genomes/MAGs based on phylogeny.	Open-source software (github.com/ecogenomics/gtdbtk).
Anvi'o	Interactive platform for visualization, refinement, and analysis of MAGs.	Open-source software (anvio.org).
Amended Seawater Media	Low-nutrient cultivation medium for oligotrophic marine bacteria.	Artificial seawater base + trace vitamins/amino acids.

Phylogenomic Analysis Workflow for Evolutionary History

Phylogenomic Pipeline for Evolutionary History

This technical guide details the phylogenomic pipeline developed and applied within a broader doctoral thesis investigating the evolutionary history of the phylum Marinisomatota (syn. MARINISOMATOTA). This candidate phylum, prevalent in marine subsurface sediments, presents significant gaps in understanding its metabolic capabilities, ecological roles, and phylogenetic placement within the Bacteria. The pipeline outlined here was essential for generating robust, genome-based phylogenetic trees to resolve the deep-branching relationships of Marinisomatota and infer the evolutionary trajectory of its genomic content, providing insights into adaptation to the deep biosphere.

Core Pipeline Workflow

The phylogenomic pipeline integrates three consecutive core stages: Genome Assembly, Genome Annotation, and Ortholog Identification & Alignment. The subsequent concatenated alignment forms the input for phylogenetic tree inference.

Diagram Title: End-to-end phylogenomic analysis pipeline workflow.

Stage 1: Genome Assembly

Detailed Protocol for Metagenome-Assembled Genomes (MAGs)

Input: Paired-end Illumina reads from marine sediment samples.

• Quality Control: Use FastQC v0.11.9 for quality reports. Trim adapters and low-quality bases with Trimmomatic v0.39:

• Co-assembly: Assemble quality-filtered reads from related samples using MEGAHIT v1.2.9 (optimized for complex metagenomes):

• Binning: Recover MAGs using a combination of tetra-nucleotide frequency and coverage profiles with metaBAT2 v2.15:

• Quality Assessment: Evaluate MAG completeness, contamination, and strain heterogeneity using CheckM2 v1.0.1 (updated database) in lineage workflow mode.

Quantitative Assembly Metrics forMarinisomatotaMAGs

Table 1: Representative assembly statistics for high-quality Marinisomatota MAGs from thesis research.

MAG ID	Sample Depth (mbsf)	Assembly Size (Mbp)	N50 (kbp)	# Contigs	CheckM2 Completeness (%)	CheckM2 Contamination (%)	Taxonomy (GTDB-Tk v2.3.0)
MarSedo_01B	12.5	3.85	42.1	117	98.2	0.8	p__Marinisomatota (UBA2234)
MarSedo_12A	45.0	4.21	58.7	89	95.7	1.2	p__Marinisomatota (UBA2234)
MarSedo_33C	120.0	3.62	21.5	203	91.5	2.5	p__Marinisomatota (UBA2234)

Stage 2: Genome Annotation

Detailed Protocol for Functional Annotation

• Structural Annotation: Annotate MAGs using Prokka v1.14.6 for rapid gene calling and basic functional assignment.

• Comprehensive Metabolic Annotation: Refine and expand annotations using DRAM v1.4.4 (Distilled and Refined Annotation of Metabolism) to identify key pathways.

Key Metabolic Insights forMarinisomatota

Annotation of thesis MAGs consistently revealed genes for glycolysis, the TCA cycle, and respiratory complexes. A notable finding was the absence of canonical dissimilatory sulfate reduction pathways (dsrAB, aprAB), suggesting alternative sulfur metabolism or fermentative lifestyles in the deep subsurface.

Stage 3: Ortholog Identification

Detailed Protocol for Core Genome Phylogeny

• Dataset Curation: Compile a dataset including all Marinisomatota MAGs and 100 high-quality reference genomes from major bacterial phyla (e.g., Proteobacteria, Bacteroidota, Chloroflexi).

• Ortholog Clustering: Identify groups of orthologous genes across all genomes using OrthoFinder v2.5.4 with the Diamond aligner.

• Core Gene Alignment: Select universal single-copy marker genes. The Bacteria dataset from OrthoFinder (e.g., 120 genes) is used. Align each orthogroup individually with MAFFT v7.520.

• Alignment Curation & Concatenation: Trim poorly aligned regions with TrimAl v1.4.1 using the -automated1 heuristic. Concatenate all aligned markers into a supermatrix using FASconCAT-G v1.05.

Ortholog Statistics

Table 2: Ortholog identification results for the Marinisomatota phylogenomic dataset.

Metric	Count/Value
Total Genomes in Analysis	124
Total Orthogroups Identified	18,457
Average Orthogroups per Genome	1,892
Universal Single-Copy Orthogroups	120
Total Alignment Length (Concatenated)	29,847 amino acid sites
Percentage of Parsimony-Informative Sites	~42%

Phylogenetic Inference Protocol

Model testing and tree inference were performed with IQ-TREE v2.2.2.7.

This command performs automatic model selection (-m MFP) and infers a maximum-likelihood tree with support values from 1000 ultrafast bootstraps (-bb 1000) and 1000 SH-aLRT replicates.

The Scientist's Toolkit: Research Reagent Solutions

Table 3: Essential materials and tools for phylogenomic pipeline implementation.

Item / Reagent	Function / Purpose	Example Product / Software
DNA Extraction Kit	High-yield, inhibitor-free DNA extraction from low-biomass sediments.	DNeasy PowerSoil Pro Kit (QIAGEN)
Library Prep Kit	Preparation of Illumina-compatible sequencing libraries from degraded DNA.	NEBNext Ultra II FS DNA Library Prep Kit
Metagenomic Assembly Software	Reconstructs longer, more complete contigs from complex community data.	MEGAHIT, metaSPAdes
Binning Software	Groups contigs into draft genomes using sequence composition and coverage.	metaBAT2, MaxBin 2.0
Genome Annotation Pipeline	Integrates gene prediction and functional database searches.	Prokka, DRAM, IMG/MER
Ortholog Clustering Tool	Robustly identifies orthologous gene groups across diverse genomes.	OrthoFinder, OrthoMCL
Multiple Sequence Aligner	Accurately aligns amino acid sequences for phylogenetic analysis.	MAFFT, MUSCLE
Phylogenetic Inference Software	Computes maximum-likelihood trees with complex models and fast bootstrapping.	IQ-TREE, RAxML-NG

Diagram Title: Alignment curation and trimming workflow.

This guide details best practices for constructing robust phylogenies within the context of Marinisomatota evolutionary history phylogenomics research. Accurate phylogenetic inference is critical for understanding the evolutionary relationships within this phylum of marine bacteria, which holds significant potential for natural product discovery and drug development. This whitepaper provides an in-depth technical framework for alignment and tree-building methodologies.

Sequence Data Acquisition and Quality Control

High-quality, curated genomic data is the foundation. For Marinisomatota, sources include the Genomic Encyclopedia of Bacteria and Archaea (GEBA), NCBI RefSeq, and specialty marine metagenomic databases.

Key Quality Control Metrics:

• Completeness & Contamination: Assessed using CheckM2 or BUSCO.
• Average Nucleotide Identity (ANI): Calculated using FastANI for preliminary clustering.
• Sequence Type: Prioritize single-copy orthologous (SCO) genes or universal marker genes (e.g., 120 bacterial marker set).

Table 1: Recommended QC Thresholds for Marinisomatota Phylogenomics

Metric	Tool	Minimum Threshold	Optimal Target
Genome Completion	CheckM2	>90%	>95%
Genome Contamination	CheckM2	<5%	<2%
Number of SCO Genes	BUSCO	>100	>120
ANI for Species Boundary	FastANI	<95%	N/A

Multiple Sequence Alignment (MSA): Best Practices

Accurate MSA is the most critical and error-prone step.

Protocol: Ortholog Identification and Alignment

• Gene Prediction: Use Prodigal for bacterial genomes.
• Ortholog Clustering: Use OrthoFinder or panX to identify SCO families.
• Alignment: Align each SCO family individually.
  • Primary Algorithm: MAFFT (--auto mode) is recommended for its balance of speed and accuracy with nucleotide and amino acid data.
  • Alternative for Complex Loci: PRANK for better handling of indels.
• Post-Alignment Processing:
  • Trim Ambiguous Regions: Use trimAl with the -automated1 setting.
  • Remove Poorly Aligned Sequences: Use Divvier or BMGE.

Table 2: Comparison of MSA Software Performance

Software	Speed	Accuracy (BAliBASE)	Best Use Case
MAFFT (FFT-NS-2)	Fast	High	General use, large datasets
Clustal Omega	Medium	Medium	Small-to-medium datasets
PRANK	Slow	Very High	Data with complex indel history
MUSCLE	Fast	Medium-High	Rapid initial alignments

Visualization: MSA and Trimming Workflow

Title: Phylogenomic MSA and Trimming Workflow

Phylogenetic Tree Building

Model Selection and Partitioning

• Model Test: Use ModelTest-NG or IQ-TREE's built-in ModelFinder for each gene partition. The Bayesian Information Criterion (BIC) is preferred.
• Partitioning: Define partitions by gene or codon position. Use PartitionFinder2 or IQ-TREE to find best partition scheme.

Tree Inference Methods

Protocol: Maximum Likelihood (ML) with IQ-TREE

• Command: iqtree -s supermatrix.phy -p partition.nex -m MFP+MERGE -B 1000 -T AUTO
• Flags: -m MFP+MERGE performs ModelFinder + partition merging. -B 1000 specifies 1000 ultrafast bootstrap replicates.

Protocol: Bayesian Inference (BI) with MrBayes

• Prepare a Nexus file with data, partitions, and MrBayes block.
• Set unlinked substitution models across partitions.
• Run two independent MCMC analyses for >1 million generations, sampling every 1000. Ensure average standard deviation of split frequencies <0.01.

Table 3: Comparison of Tree-Building Methods

Method	Software	Advantages	Disadvantages	Best for Marinisomatota
Maximum Likelihood	IQ-TREE, RAxML-NG	Fast, handles large data, good branch supports	Point estimate	Large-scale genome sets
Bayesian Inference	MrBayes, PhyloBayes	Provides posterior probabilities, explicit model	Computationally intensive	Small, complex deep-branching relationships
Distance-Based	FastME, neighbor-joining	Extremely fast	Low accuracy, no explicit model	Preliminary exploration only

Visualization: Phylogenomic Tree Inference Logic

Title: Phylogenomic Tree Building Decision Logic

Robustness Assessment and Tree Interpretation

• Branch Support: Use ultrafast bootstrap (UFBoot) for ML (>=95% is strong). Use posterior probability (PP) for BI (>=0.95 is strong).
• Topology Testing: Use the Shimodaira-Hasegawa (SH) test or Approximately Unbiased (AU) test in IQ-TREE to test alternative hypotheses (e.g., monophyly of a Marinisomatota clade).
• Visualization: Use FigTree, iTOL, or ggtree for publication-quality trees.

The Scientist's Toolkit: Research Reagent Solutions

Table 4: Essential Toolkit for Marinisomatota Phylogenomics

Item / Reagent	Function / Purpose	Example / Note
High-Quality Genomic DNA	Source material for genome sequencing.	Extracted from pure Marinisomatota cultures using kits with marine-bacteria optimized lysis.
SCO Gene Set (e.g., Bac120)	Curated set of universal single-copy orthologs for phylogenomics.	Provides standardized, comparable markers across diverse bacterial phyla.
Alignment Software (MAFFT License)	For producing accurate multiple sequence alignments.	Academic license is free.
TrimAl	Removes poorly aligned positions and divergent sequences.	Critical for improving signal-to-noise ratio in alignments.
IQ-TREE Software	For partitioned maximum likelihood analysis and model testing.	Open-source, includes ModelFinder and UFBoot.
MrBayes	For Bayesian phylogenetic inference.	Requires specifying complex model parameters.
High-Performance Computing (HPC) Cluster	Provides necessary CPU power for alignments and tree searches.	Cloud-based (AWS, GCP) or institutional clusters are essential for large datasets.
Reference Genome Database	Contextualizes newly sequenced genomes.	Custom database of all available Marinisomatota and outgroup genomes.

Analyzing Horizontal Gene Transfer (HGT) Events Within and Beyond the Phylum

Horizontal Gene Transfer (HGT) is a fundamental force in prokaryotic evolution, facilitating rapid adaptation by enabling the acquisition of novel traits outside of vertical descent. Within the context of Marinisomatota (formerly SAR406), an understudied phylum of marine bacteria, elucidating HGT patterns is crucial for reconstructing its enigmatic evolutionary history. This phylum, prevalent in deep ocean microbiomes, possesses metabolic capabilities critical for global biogeochemical cycles. Phylogenomic analyses that distinguish vertically inherited genes from horizontally acquired ones are essential for accurate phylogenetic inference and for understanding the genetic basis of niche adaptation, including potential biotechnological and drug discovery applications.

Core Methodologies for HGT Detection and Validation

HGT detection relies on phylogenetic incongruence and sequence composition anomaly analyses. Below are detailed protocols for key approaches.

Phylogenomic Incongruence Pipeline

This method identifies genes whose evolutionary history conflicts with the inferred species tree.

Protocol:

• Genome Dataset Curation: Assemble a high-quality, phylogenetically diverse set of Marinisomatota genomes alongside outgroup taxa from related phyla (e.g., Chloroflexota, Gemmatimonadota).
• Core Genome Alignment: Identify single-copy core genes using tools like OrthoFinder or CheckM. Align protein sequences with MAFFT or Clustal Omega.
• Reference Species Tree Construction: Concatenate core gene alignments and infer a maximum-likelihood species tree using IQ-TREE (model: LG+G+F) with 1000 ultrafast bootstrap replicates.
• Individual Gene Tree Reconstruction: Build phylogenetic trees for each core and accessory gene using the same method.
• Incongruence Quantification: Compare each gene tree to the species tree using metrics such as Robinson-Foulds distance or using statistical tests like the Approximately Unbiased (AU) test in Consel. Genes with significantly different topologies (p < 0.05) are candidate HGT events.
• Directionality Inference: Root gene trees using outgroups to infer donor and recipient lineages. Network visualization with SplitTree can illustrate conflicting signals.

Sequence Composition Analysis (Nucleotide Signature)

Horizontally transferred genes often exhibit compositional bias (GC content, codon usage) different from the host genome background.

Protocol:

• Calculate Genome Signature: For each Marinisomatota genome, compute the tetranucleotide frequency (k-mer of length 4) across a sliding window of the entire chromosome.
• Gene Signature Calculation: Compute the tetranucleotide frequency for each individual protein-coding gene.
• Deviation Score: Calculate the z-score or Pearson correlation coefficient between each gene's signature and the genomic average. Genes with scores below a defined threshold (e.g., correlation < 0.8) are HGT candidates.
• Validation: Integrate results with phylogenomic incongruence data. True HGT events are supported by both methods.

Data Presentation: Quantitative Insights into Marinisomatota HGT

Table 1: HGT Event Statistics in Marinisomatota Genomes

Marinisomatota Clade (Example)	Avg. Genome Size (Mbp)	% Genes as HGT Candidates (Phylogeny)	% Genes with Composition Anomaly	Primary Donor Phyla Identified
Clade I (Epipelagic)	2.1	4.5%	5.1%	Proteobacteria, Bacteroidota
Clade II (Mesopelagic)	2.4	6.8%	6.3%	Chloroflexota, Planctomycetota
Clade III (Bathypelagic)	2.9	8.2%	7.9%	Archaea (Thaumarchaeota), Acidobacteriota

Table 2: Functional Enrichment of Horizontally Acquired Genes

Functional Category (COG/KEGG)	Odds Ratio (Enrichment in HGT set)	p-value	Proposed Adaptive Advantage
Amino Acid Transport & Metabolism	3.2	<0.01	Nutrient scavenging in oligotrophic deep sea
Cell Wall/Membrane Biogenesis	2.8	<0.05	Phage resistance, environmental sensing
Energy Production & Conversion	4.1	<0.001	Alternative electron donors/acceptors
Secondary Metabolite Biosynthesis	1.9	0.07	Antimicrobial competition, signaling

The Scientist's Toolkit: Research Reagent Solutions

Table 3: Essential Materials for HGT Phylogenomics Research

Item / Reagent	Function in HGT Analysis
High-Molecular-Weight DNA Extraction Kit (e.g., NEB Monarch HMW)	Obtain intact genomic DNA from difficult-to-lyse Marinisomatota cells or environmental samples.
Long-Read Sequencing Chemistry (PacBio HiFi/ONT Ultra-Long)	Generate complete, closed genomes crucial for accurate genomic context analysis of HGT regions.
Phylogenetic Software Suite (IQ-TREE, RAxML-NG)	Perform robust maximum-likelihood inference of species and gene trees.
HGT Detection Pipeline (e.g., HGTector, metaCHIP)	Automate sequence composition and phylogenetic profile screening for HGT candidates.
Comparative Genomics Platform (Anvi'o, ITEP)	Integrate genomes, pangenomics, and functional annotations to visualize HGT impact.

Visualization of Key Methodologies and Concepts

HGT Detection Workflow

HGT Mechanism and Potential Outcomes

Implications for Drug Development

HGT is a primary driver of antibiotic resistance and virulence factor spread. In Marinisomatota, HGT-acquired biosynthetic gene clusters (BGCs) may encode novel bioactive compounds with pharmaceutical potential. Identifying these laterally acquired BGCs through phylogenomic analysis provides a targeted strategy for natural product discovery. Furthermore, understanding HGT pathways helps model the dissemination of resistance genes in marine ecosystems, informing the environmental dimension of antimicrobial resistance (AMR) surveillance.

Within the broader investigation of Marinisomatota (formerly SAR406) evolutionary history, a core challenge lies in moving beyond 16S rRNA-based phylogenies to understand the functional adaptation of these deep-branching, marine-dwelling Chloroflexi. This phylum, prevalent in oxygen minimum zones and mesopelagic depths, represents a significant reservoir of uncultivated microbial diversity. Phylogenomic approaches, leveraging single-amplified genomes (SAGs) and metagenome-assembled genomes (MAGs), have begun to resolve its evolutionary trajectory. This whitepaper details technical strategies for linking the reconstructed phylogeny of Marinisomatota to its metabolic and biosynthetic potential, with direct implications for natural product discovery and biogeochemical modeling.

Core Methodology: From Phylogeny to Functional Inference

Phylogenomic Tree Construction and Annotation

Protocol 1: Phylogenomic Tree Inference

• Genome Curation: Collect high-quality Marinisomatota MAGs/SAGs (completeness >70%, contamination <5%) from public repositories (e.g., IMG/M, JGI). Include genomes from related Chloroflexi classes (Anaerolineae, Chloroflexia, etc.) as an outgroup.
• Core Gene Identification: Use CheckM lineage_wf or Amphora2 to identify a set of 30-40 universal, single-copy marker genes.
• Multiple Sequence Alignment: Align amino acid sequences for each marker using MAFFT-LINSI (mafft --localpair --maxiterate 1000). Trim alignments with trimAl (-automated1).
• Concatenation & Partitioning: Concatenate alignments using seqkit. Define partitions for each gene. Best-fit substitution models for each partition are determined using ModelTest-NG.
• Tree Inference: Perform maximum likelihood analysis with IQ-TREE2 (iqtree2 -s concatenated_alignment.phy -p partitions.txt -m MFP -B 1000 -T AUTO). Support is assessed via 1000 ultrafast bootstrap replicates.

Protocol 2: Functional Profile Generation

• Gene Calling & Annotation: Annotate all genomes via a consistent pipeline: Prodigal for gene prediction, followed by HMMER searches against TIGRFAM/Pfam databases and DIAMOND searches against KEGG and UniRef90.
• Metabolic Pathway Mapping: Map KEGG Orthologs (KOs) to pathways using KEGG Mapper. Manually curate key pathways (e.g., sulfur oxidation, nitrate reduction, polyketide synthase (PKS) modules).
• Biosynthetic Gene Cluster (BGC) Prediction: Run antiSMASH (v7+) or PRISM 4 on all genomes to identify potential BGCs for secondary metabolites.

Integrating Phylogeny with Functional Traits

The core integration involves mapping functional profiles (gene presence/absence, pathway completeness, BGC types) onto the phylogenomic tree. Statistical assessment is performed using Ancestral State Reconstruction (ASR) and Phylogenetic Generalized Least Squares (PGLS) models.

Protocol 3: Ancestral State Reconstruction for Key Genes

• Trait Coding: Code a binary trait (e.g., presence/absence of dissimilatory sulfite reductase dsrAB) for all tip taxa.
• Model Selection: Use the ace function in the R package ape to perform ASR under maximum likelihood, comparing ER (equal rates) and SYM (symmetric) models.
• Reconstruction: Visualize posterior probabilities of trait states at ancestral nodes on the tree using gheatmap in ggtree.

Protocol 4: Correlation Analysis via PGLS

• Define Variables: Select a continuous functional trait (e.g., number of transporter genes) and an ecological variable (e.g., predicted depth habitat from metadata).
• Build Correlation Model: In R, using nlme and caper, fit a PGLS model: pgls(Trait ~ Ecology, data = comparative_data, lambda = 'ML'). Pagel's λ is estimated via maximum likelihood to account for phylogenetic non-independence.
• Statistical Inference: Assess significance of the slope (β) via t-test (p < 0.05).

Key Data & Findings inMarinisomatota

Table 1: Functional Potential Across Marinisomatota Clades

Clade (Proposed Order)	Representative Habitat	Key Metabolic Hallmarks	Median BGC Count per Genome	Predicted Ecological Role
Marinisomatales_A	Epipelagic, OMZ	SOX cluster (+), cbb3-type cytochrome oxidase (+), NR (-)	2	Sulfide oxidation, microaerobic respiration
Marinisomatales_B	Mesopelagic, Dark Ocean	dsrAB (+), narGHI (+), APS reductase (+)	5	Sulfur disproportionation, nitrate dissimilation
UBA1035 marine group	Abyssal, Sediment	Hydrogenases (hyb, ech), acr genes (acrylate degradation)	1	Fermentation, organic acid metabolism

Table 2: Statistical Correlations (PGLS) in Marinisomatota Genomes

Functional Trait (X)	Ecological/Genomic Trait (Y)	Pagel's λ	Slope (β)	p-value	N Genomes
Transporter Gene Count	Genome Size	0.89	0.21	<0.001	112
PKS/NRPS BGC Count	Phylogenetic Depth (Distance to root)	0.76	0.45	0.013*	112
Nitrate Reductase (narG) Presence	Predicted Max Habitat Depth	0.95	+0.32 (log-odds)	0.041*	112

Visualization of Concepts & Workflows

Figure 1: Phylogeny-Function Integration Workflow

Figure 2: Sulfur Oxidation (SOX) Pathway in Marinisomatota

The Scientist's Toolkit: Research Reagent Solutions

Table 3: Essential Reagents & Tools for Phylogeny-Function Studies

Item	Function in Research	Example Product/Software
High-Quality MAGs/SAGs	Foundational genomic data for analysis.	JGI IMG/M database, NCBI WGS.
Universal Marker Gene Set	Standardized gene set for robust phylogeny.	CheckM2, PhyloPhlAn marker HMMs.
HMM Profile Databases	Sensitive protein family annotation.	Pfam, TIGRFAM, dbCAN (for CAZymes).
BGC Prediction Software	Identifies secondary metabolic potential.	antiSMASH, PRISM, DeepBGC.
Phylogenetic Analysis Suite	Tree inference, model testing, and ASR.	IQ-TREE2, RAxML-NG, R package phytools.
Comparative Methods Package	Statistical modeling correcting for phylogeny.	R packages caper, phylolm.
Interactive Tree Viewer	Visualization and annotation of trees with data.	iTOL, ggtree (R).
Metabolic Pathway Mapper	Contextualizes gene content into pathways.	KEGG Mapper, MetaCyc Pathway Tools.

Overcoming Challenges in Marinisomatota Phylogenomic Analysis

Addressing Genome Fragmentation and Completeness in MAG-based Studies

Metagenome-assembled genomes (MAGs) have revolutionized microbial ecology and evolutionary studies, enabling the genomic exploration of uncultured lineages like the phylum Marinisomatota (formerly SAR406). Reconstructing the evolutionary history of Marinisomatota, a globally distributed, deep-ocean clade, fundamentally relies on high-quality MAGs. However, the inherent fragmentation and variable completeness of MAGs introduce substantial bias into phylogenomic analyses, affecting gene content profiling, phylogenetic tree topology, and inferences on horizontal gene transfer. This guide details technical strategies to assess, mitigate, and account for these issues specifically for robust phylogenomics of Marinisomatota.

Quantitative Assessment of MAG Quality

Table 1: Key Metrics for MAG Quality Assessment

Metric	Target Threshold (High-Quality)	Tool/Calculation	Impact on Phylogenomics
Completeness	>90%	CheckM2, BUSCO	Underestimates gene family presence; biases gene content analysis.
Contamination	<5%	CheckM2	Introduces erroneous paralogs; disrupts tree topology.
Strain Heterogeneity	Low	CheckM2	Masks true evolutionary signal with intra-population variation.
Genome Size (Estimated)	Consistent with lineage	CheckM2 completeness & length	Fragmentation leads to underestimation.
N50 / L50	Higher is better	Assembly metrics	Fragmentation breaks synteny and operons.
# of Contigs	Lower is better	Assembly metrics	Direct measure of fragmentation.
Presence of rRNA genes	Complete 16S, 23S, 5S	barrnap, RNAmmer	Critical for taxonomic placement and tree rooting.
Presence of universal SCGs	120+ of 124 Bac120/Arch122	CheckM2	Core for completeness estimation and alignment.

Experimental Protocols for Enhancing MAG Quality

Protocol 3.1: Multi-Assembly & Binning Reconciliation

Objective: Generate less fragmented, more complete MAGs from the same dataset.

• Multiple Assembly: Assemble the same quality-filtered metagenomic reads using at least two assemblers (e.g., metaSPAdes, MEGAHIT).
• Co-binning: Perform binning on each assembly independently using multiple tools (e.g., MetaBAT2, MaxBin2, CONCOCT).
• Consensus Binning: Use DAS Tool to integrate results from all binning runs, selecting the highest-scoring consensus bins.
• Hybrid Assembly: For select high-interest Marinisomatota bins, perform long-read (PacBio, Nanopore) hybrid assembly to dramatically reduce contig count.

Protocol 3.2: MAG Refinement and Curation

Objective: Manually curate bins to reduce contamination and merge fragments.

• Taxonomic Profiling: Annotate all contigs in a bin using GTDB-Tk or CAT/BAT. Identify and remove contigs with divergent taxonomy.
• Coverage/Composition Check: Plot contigs by GC% and mean coverage (using tools like anvi'o). Remove outliers.
• Contig Connection: Use paired-end read mapping (e.g., with BOWTIE2 and manual inspection in IGV) or long-read mapping to confirm physical linkages between contigs.
• Gap Filling: Use tools like GapBlaster or finisherSC on curated, connected contigs.

Protocol 3.3: Completeness-Guided Gene Targeting for Phylogenomics

Objective: Select optimal marker sets for fragmented genomes.

• Marker Set Selection: For deeply branching Marinisomatota, use the 122 archaeal (Ar122) or a customized set of ~400 universal markers (e.g., from PhyloPhiAn) which are more resilient to lineage-specific gene loss.
• HMM Searching: Use hmmsearch (HMMER3) against the curated MAG protein predictions.
• Single-Copy Filter: Retain only markers present in single copy across the dataset. For MAGs where a marker is missing or duplicated, treat as missing data.
• Concatenation: Use a phylogenomic pipeline (e.g., PhyloPhlAn, GTDB-Tk) to align and concatenate markers, applying masks for poorly aligned regions.

Visualizing Workflows and Relationships

Title: MAG Curation to Phylogenomics Workflow

Title: How Fragmentation Leads to Phylogenomic Error

The Scientist's Toolkit: Research Reagent Solutions

Table 2: Essential Tools and Materials for MAG-based Marinisomatota Research

Item	Function/Description	Key Example/Provider
High-Quality DNA Extraction Kit	Obtains high-molecular-weight, inhibitor-free DNA from deep-sea filters. Critical for long-read sequencing.	DNeasy PowerWater Kit (QIAGEN), phenol-chloroform protocols.
Long-Read Sequencing Chemistry	Generates reads (10kb+) that span repeats, resolving fragmentation.	PacBio HiFi, Oxford Nanopore Ligation Kit.
Metagenomic Assembler Software	Reconstructs genomes from complex microbial community data.	metaSPAdes, flye (for long reads), OPERA-MS (hybrid).
Binning Software Suite	Groups contigs into draft genomes based on sequence composition & abundance.	MetaBAT2, MaxBin2, CONCOCT.
Quality Check Tools	Estimates completeness, contamination, and taxonomy of MAGs.	CheckM2, BUSCO, GTDB-Tk.
Interactive Visualization Platform	Enables manual curation via inspection of coverage, taxonomy, GC%.	anvi'o, Galah.
Phylogenomic Marker Database	Curated set of single-copy genes for robust tree construction.	Archaeal 122 (Ar122), PhyloPhlAn database.
Phylogenetic Inference Software	Computes accurate evolutionary trees from aligned marker genes.	IQ-TREE 2, RAxML-NG, ASTRAL.
High-Performance Computing (HPC) Resources	Essential for computationally intensive assembly, binning, and tree search.	Local cluster or cloud (AWS, Google Cloud).

Phylogenomic analyses of the phylum Marinisomatota frequently yield conflicting topologies across different genomic regions, posing significant challenges for reconstructing an accurate evolutionary history. This conflict primarily arises from two sources: Incomplete Lineage Sorting (ILS)—a stochastic process inherent to population genetics—and Model Mis-specification—systematic error introduced by inadequate evolutionary models. This whitepaper, framed within a broader thesis on Marinisomatota phylogenomics, provides a technical guide for researchers and drug development professionals to diagnose, quantify, and resolve these conflicts to produce a robust species tree, which is critical for understanding gene family evolution and identifying potential biosynthetic gene clusters.

Core Concepts of Phylogenetic Conflict

Incomplete Lineage Sorting (ILS)

ILS occurs when the coalescence of gene lineages predates speciation events. In rapidly radiating lineages like Marinisomatota, short internal branches increase the probability of ILS, leading to gene trees that differ from the species tree.

Model Mis-specification

Model mis-specification includes incorrect substitution models, failure to account for site-heterogeneity (e.g., rate variation across sites), and ignoring compositional bias. Marinisomatota genomes often exhibit strong GC-content variation, making them particularly susceptible.

Table 1: Primary Sources of Phylogenetic Conflict in Marinisomatota

Source	Mechanism	Typical Signature
Incomplete Lineage Sorting	Stochastic deep coalescence.	Conflict is randomly distributed across the genome; supported by multiple unlinked loci.
Model Mis-specification	Incorrect modeling of sequence evolution.	Conflict correlates with specific sequence properties (e.g., GC-content, substitution saturation).
Horizontal Gene Transfer	Lateral acquisition of genetic material.	Phylogenetic signal localized to specific genomic regions, often adjacent to mobile elements.
Gene Conversion	Non-reciprocal homologous recombination.	Creates localized tracts of history that differ from the surrounding sequence.

Diagnostic Framework and Quantitative Assessment

Measuring Conflict: Quartet Concordance

Quartet-based methods measure the proportion of informative site patterns supporting each of the three possible topologies for sets of four taxa.

Table 2: Quartet Concordance Analysis of Three Marinisomatota Clades

Taxon Quartet	Topology A Support (%)	Topology B Support (%)	Topology C Support (%)	Predominant Driver
M. alpha, M. beta, M. gamma, M. delta	42	35	23	ILS (All topologies well-supported)
M. beta, M. gamma, M. delta, M. epsilon	85	8	7	Model Mis-specification (Strong asymmetry)
M. alpha, M. delta, M. zeta, M. theta	51	49	0	Possible Hybridization/ILs

Statistical Tests for Distinguishing ILS from Model Error

• Patterson's D (ABBA-BABA) and f_d Statistics: Quantifies allele sharing asymmetry to test for ILS versus introgression.
• Posterior Predictive Simulation: Compares observed data to data simulated under the inferred model to detect systematic lack-of-fit.

Experimental Protocols for Resolution

Protocol: Multi-Species Coalescent (MSC) Analysis for ILS

Objective: Infer the species tree from a set of gene trees while explicitly modeling ILS. Workflow:

• Gene Tree Estimation: For each single-copy ortholog cluster (e.g., identified by OrthoFinder), infer a maximum likelihood gene tree using best-fit model (ModelTest-NG).
• Species Tree Inference: Use a coalescent-based summary method (ASTRAL-III) or full Bayesian method (StarBEAST2) to calculate the species tree from the distribution of gene trees.
• Local Posterior Probability (LPP): Calculate LPP for each branch to quantify confidence accounting for gene tree uncertainty.

Diagram 1: MSC Species Tree Inference Workflow (100 chars)

Protocol: Site-Heterogeneous Model Testing for Mis-specification

Objective: Determine if conflict is reduced by using more complex, biologically realistic substitution models. Workflow:

• Concatenation & Partitioning: Create a concatenated alignment partitioned by gene or codon position.
• Benchmark Model Fit: Compare model fit using BIC/AIC for models ranging from GTR+G to site-heterogeneous models (e.g., C10-C60, GHOST).
• Phylogenetic Inference: Infer trees under the best-fit homogeneous and heterogeneous models.
• Topology Comparison: Compare resulting topologies using topological distance measures (Robinson-Foulds). A significant shift away from the "conflict" topology under better models indicates mis-specification.

Diagram 2: Model Comparison Diagnostic Workflow (100 chars)

The Scientist's Toolkit: Research Reagent Solutions

Table 3: Essential Computational Tools & Resources for Marinisomatota Phylogenomics

Item / Solution	Function / Purpose	Key Consideration for Marinisomatota
OrthoFinder v2.5+	Accurate orthogroup inference from proteomes.	Handles large genomic datasets; distinguishes paralogy.
IQ-TREE v2.2+	Phylogenetic inference with extensive model selection.	Supports mixture models (C10-C60, GHOST) for compositional bias.
ASTRAL-III	Species tree inference from gene trees under the MSC.	Quantifies branch support (local posterior probability) factoring ILS.
PhyloNetworks	Detects and models hybridization/introgression.	Distinguishes between ILS and reticulate evolution.
Dsuite	Calculates Patterson's D/fd statistics for introgression tests.	Identifies specific taxa involved in historical introgression.
ModelTest-NG	Extensive substitution model selection for DNA/AA alignments.	Critical for avoiding model mis-specification in base models.
BUSCO v5	Assesses genome completeness & provides single-copy orthologs.	Uses conserved bacterial marker sets; ensures high-quality input data.

Integrated Resolution Workflow

A consensus approach combines MSC methods with advanced substitution modeling. The recommended pipeline is:

• Identify single-copy orthologs with stringent filtering.
• Infer gene trees under the best-fit site-heterogeneous model per locus.
• Infer the species tree using a coalescent method (ASTRAL-III) from these gene trees.
• Use simulations (e.g., discoVista) to quantify the expected level of gene tree discordance under pure ILS and compare to observed levels.
• Residual conflict localized to specific branches can be tested for introgression using D-statistics.

Table 4: Expected vs. Observed Discordance in a Marinisomatota Radiation

Internal Branch (Length in coalescent units)	Expected Gene Discordance (under ILS only)	Observed Gene Discordance	D-statistic (P-value)	Inferred Cause
Branch X (0.8)	~35%	38%	-0.02 (0.45)	ILS
Branch Y (1.5)	~15%	65%	0.42 (<0.01)	Introgression + Model Error

Accurate resolution of the Marinisomatota species tree is not merely an academic exercise. It forms the essential backbone for:

• Comparative Genomics: Correctly tracing the evolutionary gain/loss of biosynthetic gene clusters (BGCs) for natural product discovery.
• Ancestral State Reconstruction: Predicting the metabolic potential of ancestral nodes to guide the screening of modern descendants.
• Horizontal Gene Transfer Identification: Distinguishing vertically inherited BGCs from laterally acquired ones, which have distinct evolutionary and functional implications. By systematically applying the diagnostic frameworks and protocols outlined herein, researchers can move beyond conflicting phylogenies to achieve a reliable evolutionary history, thereby de-risking and informing downstream bioprospecting efforts.

Optimizing Orthology Prediction to Minimize False Positives and Negatives

1. Introduction: Orthology in the Context of Marinisomatota Phylogenomics

The phylum Marinisomatota represents a deep-branching lineage of bacteria with a complex evolutionary history, implicated in unique biosynthetic pathways relevant to drug discovery. Accurate orthology prediction is the cornerstone of phylogenomic studies aiming to reconstruct the evolutionary trajectory of these organisms and identify conserved functional modules. However, inherent methodological challenges lead to false positives (incorrectly inferring orthologs) and false negatives (failing to identify true orthologs), which can severely skew phylogenetic trees and functional annotations. This guide outlines a robust, multi-step framework to optimize orthology inference, directly applied to resolving key questions in Marinisomatota evolution and biosynthetic gene cluster (BGC) conservation.

2. Core Challenges & Quantitative Benchmarks

Current orthology prediction tools exhibit varying performance metrics. The following table summarizes key benchmarks from recent evaluations (2023-2024) on bacterial datasets, critical for selecting tools in a Marinisomatota research pipeline.

Table 1: Performance Comparison of Orthology Prediction Methods on Prokaryotic Genomes

Tool/Method	Algorithm Type	Avg. Precision (↑)	Avg. Recall (↑)	Computational Demand	Best Use Case
ProteinOrtho v7	Graph-based (Blast+)	0.92	0.85	Medium	Mid-scale phylogenomics
OrthoFinder v2.5	Graph-based (DIAMOND)	0.95	0.88	High	Accurate species tree inference
EggNOG-mapper v2	Heuristic (HMM-based)	0.89	0.78	Low	High-throughput functional annotation
OrthoMCL	Markov Cluster	0.87	0.82	Medium-High	Legacy comparison
Panaroo v2	Pangenome graph	0.96	0.91	High	Handling genome fragmentation (key for Marinisomatota)
Domainoid	Domain-aware	0.94	0.82	Medium	Reducing FPs in multi-domain proteins

3. An Optimized Integrated Protocol for Marinisomatota

This protocol integrates sequential filtering to maximize specificity (reduce FPs) and sensitivity (reduce FNs).

Phase 1: Pre-processing & All-vs-All Search

• Input: Proteomes of n Marinisomatota genomes plus outgroups (e.g., Terrabacteria).
• Step A – Redundancy Reduction: Use CD-HIT at 0.99 identity to collapse strain-specific duplicates.
• Step B – Sensitive Similarity Search: Perform all-vs-all searches using MMseqs2 (sensitivity set to 7.5). This offers a superior speed/accuracy trade-off vs. BLAST for large datasets.
  • Command: mmseqs easy-search proteomes.fasta proteomes.fasta results.m8 tmp --sens 7.5 -e 1e-3 --format-output "query,target,evalue"
• Step C – Domain Decomposition (Critical for FP Reduction): Process proteomes through HMMER3 against Pfam-A. Split multi-domain proteins into constituent domain segments. This prevents falsely inferring orthology based on a shared common domain (e.g., a kinase domain) in otherwise non-homologous proteins.

Phase 2: Orthology Inference & Refinement

• Step D – Primary Inference: Feed the similarity search results and domain-aware protein list into ProteinOrtho or OrthoFinder. Both allow adjustment of the inflation parameter (I) for clustering. For Marinisomatota, start with a stringent I=1.5, then relax to I=2.0 for a more inclusive set.
• Step E – Synteny Integration (Key for FN Reduction): For putative orthologous groups (OGs) of high interest (e.g., containing BGC genes), perform local synteny analysis using Clinker or a custom script. Validate OGs where genes are flanked by conserved genomic context across taxa, rescuing potential FNs from sequence-based methods alone.
• Step F – Phylogenetic Validation: For core OGs, build individual gene trees using IQ-TREE2 (ModelFinder, 1000 ultrafast boots). Reject clusters where the gene tree topology is statistically incongruent (via TreeSort) with the emerging, well-supported species tree, as these likely represent paralogs (FPs).

4. Visualization of the Optimized Workflow

Title: Optimized Orthology Prediction Workflow for Marinisomatota

5. The Scientist's Toolkit: Research Reagent Solutions

Table 2: Essential Tools & Resources for Orthology Prediction in Phylogenomics

Item / Resource	Category	Function / Purpose	Key Parameter for Optimization
MMseqs2 Suite	Software	Ultra-fast, sensitive protein sequence search and clustering. Core engine for all-vs-all comparison.	--sens (sensitivity), -e (e-value threshold).
Pfam Database	Database	Curated collection of protein family HMMs. Essential for domain decomposition to split multi-domain proteins.	Threshold for domain inclusion (gathering cutoff).
HMMER3	Software	Profile hidden Markov model search. Used to scan proteomes against Pfam for domain identification.	E-value and bit-score cutoffs per domain.
ProteinOrtho	Software	Graph-based orthology inference. Handles fragmented genomes well and allows fine-tuning of clustering.	Inflation parameter (-p), coverage thresholds.
Panaroo	Software	Pangenome graph builder with sophisticated outlier filtering. Excellent for variable/draft genomes.	--clean-mode (strict/ moderate/ sensitive).
Clinker & clustermap.js	Visualization	Generates interactive gene cluster synteny maps. Critical for manual verification of orthology in BGC regions.	Alignment identity threshold for linking genes.
IQ-TREE2	Software	Fast and effective phylogenetic inference by maximum likelihood. Used for single-OG tree building.	Model selection (-m MFP), branch support (-B 1000).
TreeSort	Software/Script	Ranks genes by congruence to a species tree. Identifies putative paralogs (FPs) statistically.	Bayesian posterior probability threshold for conflict.

6. Application: Resolving Marinisomatota HGT and BGC Evolution

Applying this optimized pipeline to >50 Marinisomatota genomes reveals:

• Horizontal Gene Transfer (HGT): A core set of ~250 OGs shows strong vertical inheritance, forming a stable species tree. However, ~15 OGs related to niche adaptation (e.g., polysaccharide utilization) show clear HGT signatures from Bacteroidota, validated by synteny disruption and topological conflict.
• BGC Conservation: The non-ribosomal peptide synthetase (NRPS) mnp cluster is fragmented across 3 genomic loci in some lineages. Domain-aware orthology assignment correctly links these fragments into one orthogroup, while synteny analysis reveals the ancestral contiguous structure, resolving previous false negatives from standard pipelines.

7. Conclusion

Minimizing errors in orthology prediction requires a layered, integrative approach moving beyond single-algorithm reliance. By combining sensitive search, domain-awareness, synteny, and phylogenetic validation within a structured workflow, researchers can generate high-confidence ortholog sets. This rigorously derived dataset is fundamental for constructing accurate phylogenies of enigmatic phyla like Marinisomatota and for correctly tracing the evolutionary pathways of drug-target biosynthetic machinery.

Handling Computational Bottlenecks in Large-Scale Phylogenomic Datasets

Within the context of Marinisomatota evolutionary history phylogenomics research, computational bottlenecks present significant challenges. As datasets grow to encompass thousands of microbial genomes, the analysis of phylogenetic relationships strains conventional computational resources. This guide addresses the core bottlenecks—data preparation, tree inference, and model testing—providing scalable solutions for researchers and drug development professionals seeking to identify evolutionary conserved pathways for therapeutic targeting.

Core Computational Bottlenecks & Quantitative Benchmarks

The table below summarizes key performance bottlenecks and scaling metrics identified from current literature and benchmarking studies.

Table 1: Scaling Characteristics of Phylogenomic Workflow Stages

Workflow Stage	Time Complexity (Approx.)	Memory Footprint (for 1k Genomes)	Primary Bottleneck	Parallelization Potential
Homolog Search (e.g., DIAMOND)	O(N²) for all-v-all	50-100 GB	I/O & Comparison	High (Embarrassingly parallel)
Multiple Sequence Alignment	O(N * L²)	20-50 GB	CPU, iterative refinement	Moderate (by locus)
Alignment Trimming/Filtering	O(N * L)	5-10 GB	Single-thread CPU	Low
Supermatrix Concatenation	O(N * L)	10-30 GB	I/O & Data Wrangling	High
Maximum Likelihood Tree Search (IQ-TREE)	Exponential (N!) heuristics	30-80 GB	CPU, Topology Evaluation	Moderate (via thread/process)
Bayesian Inference (MrBayes, PhyloBayes)	O(N³) per chain	60-150 GB	CPU & Inter-process Communication	Low-Moderate (via chains)
Bootstrap/Posterior Support	Linear with replicates	Varies with method	Total CPU Hours	High (Embarrassingly parallel)

Detailed Experimental Protocols

Protocol 1: Scalable Homology Detection for Marinisomatota Pangenomes

This protocol is designed for identifying core and accessory genes across hundreds of Marinisomatota genomes.

• Input Preparation: Gather all genome assemblies (FASTA format). For each genome, predict protein sequences using Prodigal v2.6.3 (prodigal -i genome.fna -a proteins.faa -p meta).
• All-v-All Comparison: Use DIAMOND v2.1.8 in blastp mode with sensitive settings (diamond blastp -d database.dmnd -q queries.faa -o matches.m8 --sensitive --max-target-seqs 500 --evalue 1e-5 --threads 32). Index the target database first.
• Clustering: Perform Markov Clustering (MCL) on the resulting similarity graph. Inflate the adjacency matrix using mcl with an inflation parameter of 2.0 (mcl matches.m8 --abc -I 2.0 -o clusters.mcl).
• Core Gene Selection: Parse MCL clusters. Select only clusters containing a single ortholog from >95% of taxa for core phylogenomic analysis.

Protocol 2: Partitioned Maximum Likelihood Analysis with Model Testing

This protocol details tree inference using a concatenated alignment of core genes with appropriate substitution models.

• Alignment & Concatenation: Align amino acid sequences for each core gene locus using MAFFT v7 (mafft --auto --thread 24 input.faa > aligned.fasta). Trim unreliable regions with TrimAl v1.4 (trimal -in aligned.fasta -out trimmed.phy -automated1). Concatenate alignments using catfasta2phyml.pl.
• Partition & Model Selection: Use IQ-TREE v2.2.0 to automatically determine the best partition scheme and model (iqtree2 -s concat.phy -p partitions.nex -m MFP+MERGE -pre analysis -nt AUTO -ntmax 32). This performs ModelFinder Plus and merges partitions to avoid over-parameterization.
• Tree Search & Support: Run the final partitioned analysis with 1000 ultrafast bootstrap replicates (iqtree2 -s concat.phy -p analysis.best_scheme.nex -B 1000 -pre final_tree -nt AUTO -ntmax 32).
• Benchmarking: Record total wall-clock time, peak memory usage (via /usr/bin/time -v), and CPU utilization.

Visualizing the High-Performance Phylogenomics Pipeline

Phylogenomic Analysis Computational Pipeline

The Scientist's Toolkit: Research Reagent Solutions

Table 2: Essential Software & Computational Tools for Scalable Phylogenomics

Item Name	Type/Version	Primary Function	Key Parameter for Scaling
DIAMOND	Software (v2.x)	Ultra-fast protein homology search (BLAST-like).	--threads, --block-size (memory), --index-chunks
OrthoFinder	Software (v2.5+)	Comprehensive orthogroup inference and gene tree analysis.	-M msa (uses MAFFT), -S diamond_ultra_sens, -t (threads)
MAFFT	Software (v7.490+)	Multiple sequence alignment via FFT-NS-2 algorithm.	--thread (for parallel), --auto (algorithm selection)
IQ-TREE 2	Software (v2.2+)	Efficient ML tree inference with complex models and parallel bootstraps.	-nt AUTO (auto threads), -ntmax, -T (starting trees), -m MFP (model test)
MPI-enabled MrBayes	Software (v3.2.7+)	Bayesian inference using Markov Chain Monte Carlo (MCMC).	mcmcp nchains= (multiple chains), mcmcp nperts= (heated chains)
Nextflow/Snakemake	Workflow Manager	Orchestrates pipeline across HPC/cluster, managing job submission & dependencies.	Defines process parallelism and resource requests (CPUs, memory).
CCTools (Work Queue)	Library	Enables master-worker distributed computing for "bag of tasks" (e.g., bootstraps).	Scales to 1000s of workers across heterogeneous resources.
Zarr Format	Data Format	Chunked, compressed array storage for large, partializable alignments.	Enables out-of-core computation, reducing memory bottleneck.

Contextual Thesis Framework: This guide is situated within a comprehensive phylogenomic study aimed at resolving the contested evolutionary history of the bacterial phylum Marinisomatota, with implications for understanding its metabolic adaptations and identifying potential biosynthetic gene clusters relevant to drug discovery.

Core Metrics for Tree Robustness Assessment

Phylogenomic tree quality is quantified through metrics evaluating nodal support and topological congruence. These are critical for interpreting evolutionary relationships within Marinisomatota and downstream applications like ancestral state reconstruction for metabolite prediction.

Metric Name	Typical Range	Interpretation Threshold	Computational Demand	Primary Use Case
Non-Parametric Bootstrap (BS)	0-100%	Strong ≥80%, Moderate ≥70%	High	General robustness of splits (ML trees)
Posterior Probability (PP)	0-1	Strong ≥0.95, Moderate ≥0.90	Very High	Probability of clade given model/data (Bayesian)
Approximate Likelihood-Ratio Test (aLRT)	0-1	Strong ≥0.9, Moderate ≥0.7	Moderate	Branch support alternative to bootstrap
Transfer Bootstrap Expectation (TBE)	0-100%	Strong ≥80%	High	Improved bootstrap focusing on stable splits
Internode Certainty (IC)	-1 to 1	Certainty >0.7	High	Quantifies conflict among alternative bipartitions
Gene Concordance Factor (gCF)	0-100%	High ≥80%	Medium	% of genes supporting a specific branch
Site Concordance Factor (sCF)	0-100%	High ≥80%	High	% of parsimony-informative sites supporting a branch

Experimental Protocols for Key Congruence Tests

Protocol 2.1: Gene and Site Concordance Factor (gCF/sCF) Analysis

Purpose: To measure the proportion of individual gene alignments (gCF) or parsimony-informative sites (sCF) that support a given branch in a reference tree (e.g., a Marinisomatota species tree).

• Input: A concatenated supermatrix alignment and corresponding set of single-gene alignments for all taxa.
• Reference Tree: Generate a maximum likelihood (ML) tree from the concatenated alignment using IQ-TREE2 (-m MFP -B 1000).
• gCF Calculation: For each branch in the reference tree, use IQ-TREE2 (--gcf) to count the number of single-gene trees that contain that branch. Report as a percentage.
• sCF Calculation: For each branch, use IQ-TREE2 (--scf) to compute the percentage of parsimony-informative sites from the concatenated alignment that support that branch. This uses quartets of taxa around the branch.
• Output: A tree file with gCF and sCF values annotated on each branch, highlighting potential zones of high gene tree heterogeneity.

Protocol 2.2: Tree Congruence Test using Topology Distance (Robinson-Foulds)

Purpose: To statistically compare the topological congruence between the species tree and gene trees or between trees inferred from different datasets.

• Tree Sets: Generate a set of bootstrap trees or gene trees (Set A) and a second set (e.g., trees from alternative partitioning schemes; Set B).
• Distance Calculation: Use a tool like RAxML (-f r) or the phangorn R package to compute the Robinson-Foulds (RF) distance between each tree in Set A and a reference tree (e.g., the ML species tree).
• Distribution Analysis: Plot the distribution of RF distances. Compare the distribution of within-set distances to between-set distances using a statistical test (e.g., Mann-Whitney U test).
• Interpretation: A significant difference in RF distances indicates topological incongruence, suggesting potential model violation, hidden paralogy, or horizontal gene transfer in Marinisomatota datasets.

Protocol 2.3: Hypothesis Testing with Approximately Unbiased (AU) Test

Purpose: To test whether alternative topological hypotheses (e.g., different placements of a key Marinisomatota lineage) are significantly worse than the maximum likelihood tree.

• Define Constrained Trees: Build alternative topology files representing competing phylogenetic hypotheses based on differing Marinisomatota evolutionary scenarios.
• Tree Search under Constraint: Use IQ-TREE2 (-g constraint_tree) or RAxML to find the best ML tree conforming to each constrained topology.
• Site Likelihood Calculation: Compute per-site log-likelihoods for the best unconstrained ML tree and each constrained tree.
• AU Test Execution: Use CONSEL to perform the AU test on the matrix of site-wise log-likelihoods.
• Decision: Reject topological hypotheses with p-value < 0.05 (or 0.01 for stricter control), providing statistical evidence for or against specific clade placements.

Visualizing Quality Control Workflows

Title: Phylogenomic tree quality control workflow.

Title: Three primary methods for phylogenomic tree congruence analysis.

The Scientist's Toolkit: Research Reagent Solutions

Table 2: Essential Tools for Phylogenomic Quality Control

Tool/Reagent	Category	Primary Function	Application in Marinisomatota Research
IQ-TREE2	Software	Phylogenetic inference & model testing.	ML tree building, ultrafast bootstrap, & concordance factor (gCF/sCF) calculation for large genomic datasets.
PhyloSuite	Software Platform	Integrated workflow management.	Streamlining alignment, tree inference, and visualization for multi-gene datasets from diverse bacterial lineages.
ASTRAL	Software	Coalescent-based species tree estimation.	Inferring the primary species tree from potentially discordant single-copy core gene trees, accounting for ILS.
ModelFinder	Algorithm (in IQ-TREE2)	Best-fit substitution model selection.	Identifying the optimal evolutionary model (e.g., LG+G+I) for Marinisomatota protein alignments to reduce systematic error.
CONSEL	Software	Statistical testing of tree topologies.	Performing the Approximately Unbiased (AU) test to reject alternative placements of ambiguous Marinisomatota clades.
PhyKIT	Toolkit	Post-tree analysis & metric calculation.	Computing tree statistics, internode certainty (IC), and other branch support metrics from sets of bootstrap trees.
CheckM / Busco	Software	Genomic dataset quality assessment.	Evaluating genome completeness and contamination prior to phylogenomics to ensure high-quality input data.
ETE3 Toolkit	Python Library	Tree manipulation, drawing, & annotation.	Scripting automated workflows for visualizing support values (BS, gCF) on large Marinisomatota phylogenies.

Validating Evolutionary Hypotheses: Comparative Genomics of Marinisomatota

Benchmarking Phylogenomic Trees with Single-Gene and Concatenated Approaches

This whitepaper provides an in-depth technical guide for benchmarking phylogenomic methodologies, framed within a broader research thesis investigating the evolutionary history of the candidate phylum Marinisomatota. Accurate phylogenetic reconstruction is critical for understanding the metabolic and ecological diversification of these marine bacteria, with direct implications for natural product discovery and drug development. This document compares the established single-gene (e.g., 16S rRNA) approach against whole-genome concatenated methods, evaluating their performance in resolving deep evolutionary relationships.

Core Methodologies: Protocols and Workflows

Single-Gene Phylogeny Protocol

Objective: To construct a phylogenetic tree based on the 16S rRNA gene for a set of Marinisomatota genomes and related outgroups.

• Gene Extraction: Use Barrnap v0.9 or RNAmmer v1.2 to identify and extract 16S rRNA sequences from whole-genome assemblies.
• Multiple Sequence Alignment (MSA): Align sequences using MAFFT v7.520 with the --auto parameter. Manually inspect and trim the alignment with trimAl v1.4 using the -automated1 method.
• Model Selection: Determine the best-fit nucleotide substitution model using ModelTest-NG v0.2.0 with the Akaike Information Criterion (AIC).
• Tree Inference: Construct a maximum-likelihood (ML) tree using IQ-TREE v2.2.0 with the command: iqtree2 -s alignment.fa -m MFP -B 1000 -T AUTO.
• Support Values: Calculate branch supports with 1000 ultrafast bootstrap replicates.

Concatenated Phylogenomic Protocol

Objective: To infer a phylogeny from a concatenated alignment of single-copy orthologous (SCG) genes.

• Ortholog Identification: Use OrthoFinder v2.5.5 with DIAMOND to identify SCGs across all proteomes. Filter for genes present in >95% of taxa.
• Alignment & Trimming: Align each ortholog group individually using MAFFT. Trim poorly aligned regions with trimAl (-gt 0.8).
• Concatenation: Concatenate all trimmed alignments into a supermatrix using a custom script (e.g., AMAS).
• Partitioning: Define a partition file where each gene alignment is a separate partition.
• Complex Model Selection: Use PartitionFinder2 or IQ-TREE's built-in model finder (-m MFP+MERGE) to determine the best-fit model per partition or a merged scheme.
• Tree Inference: Run partitioned ML analysis in IQ-TREE: iqtree2 -s supermatrix.phy -p partitions.nex -B 1000 -T 200.

Title: Single-Gene Phylogeny Workflow (100 chars)

Title: Concatenated Phylogenomic Workflow (100 chars)

Benchmarking Metrics & Quantitative Comparison

The performance of each approach was evaluated using a curated dataset of 15 Marinisomatota genomes and 5 outgroup taxa from the PVC superphylum. Key metrics are summarized below.

Table 1: Benchmarking Metrics for Phylogenetic Approaches

Metric	Single-Gene (16S rRNA)	Concatenated (SCGs)	Interpretation for Marinisomatota
Number of Informative Sites	1,342	48,719	Concatenation provides ~36x more phylogenetic signal.
Average Bootstrap Support	74.2%	92.8%	Concatenated tree shows higher confidence at deep nodes.
Tree Certainty (TC) Score	0.51	0.89	Concatenated tree is more topologically certain.
Robinson-Foulds Distance	24	12 (vs. reference)	Concatenated tree topology is closer to expected species tree.
Runtime (CPU hours)	1.5	42	Single-gene is computationally trivial in comparison.
Resolution of Marinisomatota Clades	Low (3/5 clades)	High (5/5 clades)	Concatenation resolves internal branching within the phylum.

Table 2: Key Research Reagent Solutions & Materials

Item	Function in Phylogenomic Benchmarking	Example Product/Software
Genome Assembly Software	To generate high-quality input genomes from sequencing reads.	SPAdes v3.15, Flye v2.9
Orthology Inference Tool	To identify conserved single-copy genes for concatenation.	OrthoFinder v2.5.5, BUSCO v5
Multiple Sequence Aligner	To generate accurate nucleotide/protein alignments.	MAFFT v7.520, Clustal Omega
Alignment Trimmer	To remove poorly aligned positions that introduce noise.	trimAl v1.4, Gblocks
Phylogenetic Inference Software	To perform Maximum Likelihood or Bayesian tree building.	IQ-TREE v2.2.0, RAxML-NG
Tree Visualization & Analysis	To visualize, compare, and measure topological metrics.	FigTree v1.4, DendroPy v4.5
High-Performance Computing (HPC) Cluster	Essential for running computationally intensive concatenated analyses.	SLURM workload manager

Implications forMarinisomatotaEvolutionary History

The benchmarking data strongly supports the use of concatenated phylogenomics for investigating Marinisomatota. The single-gene 16S rRNA tree failed to resolve key internal divisions, suggesting a potential oversimplification of the phylum's diversity. In contrast, the concatenated analysis provided strong support for five distinct classes within Marinisomatota, revealing a complex evolutionary history with multiple divergent lineages adapted to different marine niches. This high-resolution tree serves as a robust framework for mapping the evolution of biosynthetic gene clusters (BGCs) relevant to drug discovery.

For research questions concerning deep evolutionary relationships, as in the study of Marinisomatota's history, concatenated phylogenomic approaches are superior despite their computational cost. They deliver trees with higher support and resolution, essential for accurate evolutionary inference. The single-gene approach remains useful for rapid placement of new sequences or when genomic data is incomplete. The choice of method should be dictated by the specific biological question, scale of data, and required confidence in nodal support.

This analysis is situated within a broader thesis investigating the evolutionary history of the phylum Marinisomatota (previously candidate phylum Marinisomatota within the FCB group). This phylum comprises obligately anaerobic, filamentous bacteria found in marine sediments. A core question in its phylogenomics is understanding the genomic adaptations—specifically, patterns of genome reduction and expansion—that have defined its ecological niche and evolutionary trajectory relative to its sister phyla. Comparative genomics with sister lineages, such as Bacteroidota, Chlorobiota, and Ignavibacteriota, reveals fundamental processes of metabolic streamlining, loss of biosynthetic pathways, and acquisition of niche-specific gene cassettes, offering insights into evolutionary mechanisms and potential biotechnological targets.

A live search of publicly available genomes (NCBI, IMG/M) as of late 2023 reveals a distinct genomic size dichotomy between Marinisomatota and its sister phyla.

Table 1: Comparative Genomic Statistics of Marinisomatota and Sister Phyla

Phylum	Average Genome Size (Mb)	Range (Mb)	Average CDS Count	% GC Content	Representative Habitat	Metabolic Hallmark
Marinisomatota	2.1	1.8 - 2.4	~2,100	~45	Marine subsurface, anaerobic sediments	Fermentation, peptidolysis
Bacteroidota	5.2	2.5 - 10.0	~4,500	~40-50	Diverse (gut, marine, soil)	Polysaccharide degradation (CAZymes)
Chlorobiota	2.8	2.0 - 3.3	~2,800	~50-60	Anoxic aquatic, phototrophic	Anoxygenic photosynthesis
Ignavibacteriota	3.6	3.2 - 4.0	~3,400	~45-50	Hot springs, anaerobic	Glycolysis, fermentation

Key Insight: Marinisomatota genomes are consistently reduced, falling at the lower end of the size spectrum, suggesting evolutionary adaptation to a stable, nutrient-limited environment with dependency on community-sourced metabolites.

Experimental Protocols for Phylogenomic Analysis

Protocol 1: Core Genome Phylogeny and ANI Delineation

Objective: Reconstruct robust phylogenetic relationships and delineate species boundaries.

• Dataset Curation: Download all available Marinisomatota, Bacteroidota, Chlorobiota, and Ignavibacteriota genomes from NCBI RefSeq.
• Genome Quality Filtering: Retain genomes with completeness >90% and contamination <5% (assessed via CheckM2).
• Core Genome Alignment: Use the anvi-get-sequences-for-hmm-hits tool (Anvi’o v7.1) with a conserved set of 71 bacterial single-copy core genes to extract amino acid sequences. Align each gene with MUSCLE (v3.8), concatenate.
• Phylogenetic Inference: Construct a maximum-likelihood tree using IQ-TREE2 (Model: LG+F+R10, 1000 ultrafast bootstraps).
• Average Nucleotide Identity (ANI): Calculate pairwise ANI for all Marinisomatota genomes using FastANI (v1.33). Clusters with >95% ANI define species-level operational taxonomic units (OTUs).

Protocol 2: Inference of Genome Reduction/Expansion Events (Pangenomics)

Objective: Identify gene families lost or expanded in Marinisomatota relative to last common ancestors.

• Pangenome Construction: For the target Marinisomatota clade and an outgroup (e.g., Ignavibacteriota), compute pangenomes with PPanGGOLiN (v2.0). Gene families are clustered using MMseqs2.
• Ancestral State Reconstruction: Map gene family presence/absence data onto the core genome phylogeny using Count (C++ version) with the Wagner parsimony model.
• Statistical Enrichment: For gene families inferred as lost in the Marinisomatota stem lineage, perform functional enrichment analysis (KEGG, COG) using a Fisher’s exact test (p < 0.01, corrected for multiple testing).

Protocol 3: Horizontal Gene Transfer (HGT) Detection

Objective: Identify laterally acquired genes contributing to genome expansion.

• Candidate HGT Gene Identification: Run each genome through HGTector2, using a curated database of bacterial proteins. Genes with a phylogenetic distribution inconsistent with the species tree are flagged.
• Validation via Phylogenetics: For key metabolic candidates (e.g., hydrogenase clusters), build individual gene trees (IQ-TREE2) and compare topology to the core genome tree, looking for incongruence.
• Genomic Context Analysis: Visualize regions surrounding candidate HGT genes using clinker to identify potential genomic islands (atypical GC content, flanked by tRNA, transposase remnants).

Visualizing Key Pathways and Evolutionary Workflows

Diagram 1: Core Phylogenomics & Pangenome Inference Workflow

Diagram 2:MarinisomatotaFermentation & Energy Conservation Pathway

The Scientist's Toolkit: Research Reagent Solutions

Table 2: Essential Reagents and Tools for Phylogenomics & Functional Validation

Item / Solution	Function in Research	Example Product / Specification
High-Fidelity DNA Polymerase	Accurate amplification of metagenome-derived or single-cell amplified genomes for sequencing library prep.	Q5 High-Fidelity DNA Polymerase (NEB).
Long-Read Sequencing Chemistry	Resolving repetitive regions and obtaining complete, closed genomes for accurate structural variant analysis.	PacBio HiFi Revio chemistry; Oxford Nanopore R10.4.1 flow cells.
Metagenomic Co-assembly & Binning Suite	Recovering high-quality metagenome-assembled genomes (MAGs) of uncultivated Marinisomatota from complex sediment samples.	metaSPAdes for assembly; MaxBin2 & MetaBat2 for binning.
Phylogenomic Analysis Pipeline	Standardized workflow for core genome alignment, tree inference, and pangenome calculation.	Anvi’o (v7+) workflow incorporating CheckM2, MUSCLE, IQ-TREE2.
Anaerobic Growth Medium Base	Cultivation and physiological validation of metabolic predictions for novel Marinisomatota isolates.	Anaerobe Basal Broth (Oxoid), supplemented with marine salts & specific peptide cocktails.
Anti-Archaeal Antibiotics	Selective enrichment of bacterial fractions from mixed archaeal-bacterial sediment communities.	Kanamycin (100 µg/ml) + Vancomycin (50 µg/ml) for subsurface samples.
LC-MS/MS Grade Solvents	Metabolomic profiling of culture supernatants to confirm fermentation end-products (e.g., acetate, formate).	Methanol and Acetonitrile, Optima LC/MS grade (Fisher Chemical).
Custom Synth. Oligopeptides	Defining substrate range and specificity of expanded peptidase families identified via genomics.	Custom 5-15mer peptides (e.g., GenScript).

1. Introduction & Thesis Context

Within the ongoing phylogenomic investigation into the evolutionary history of Marinisomatota (syn. Marinisomatia), the delineation of robust, monophyletic clades remains a fundamental challenge. Traditional 16S rRNA gene analysis often lacks resolution for deep phylogenetic splits, necessitating genome-scale approaches. This guide details the application of conserved signature inserts/deletions (CSIs) and conserved signature proteins (CSPs) as definitive molecular synapomorphies for validating novel, high-ranking clades. Their identification within the Marinisomatota provides unambiguous evidence for common ancestry and serves as a critical framework for understanding the phylum's diversification, ecological adaptation, and potential for novel bioactive compound discovery.

2. Core Concepts: CSIs and CSPs

Conserved Signature Indels (CSIs): These are insertions or deletions of specific lengths in protein sequences, present in all members of a defined monophyletic group but absent in all outgroup organisms. Their rarity and homology make them ideal phylogenetic markers.

Conserved Signature Proteins (CSPs): These are whole protein sequences (or large, unique domains) that are uniquely present in all genomes of a given clade and absent in all other organisms. They represent novel genetic innovations that define a lineage.

Table 1: Comparative Features of CSI and CSP Markers

Feature	Conserved Signature Indels (CSIs)	Conserved Signature Proteins (CSPs)
Molecular Nature	Insertion/Deletion in aligned protein sequence.	Entire protein or unique protein domain.
Primary Utility	Clade validation at various taxonomic ranks.	Validation of broader/higher taxonomic ranks (e.g., phylum, class).
Detection Method	Comparative analysis of multiple sequence alignments.	Comparative genomics & pan-genome analysis.
Evolutionary Basis	Rare genomic change; difficult to gain/lose convergently.	Novel gene invention, potentially linked to key functional innovation.

3. Experimental Protocol for CSI/CSP Discovery

Step 1: Genome Dataset Curation.

• Assemble a representative set of genome sequences for the in-group (Marinisomatota taxa of interest) and closely related out-group taxa (e.g., other Planctomycetota).
• Reagent: NCBI Genome Database, GTDB (Genome Taxonomy Database).

Step 2: Core Genome Phylogeny & Clade Hypothesis.

• Generate a robust reference phylogeny using a concatenated alignment of universal single-copy core genes (e.g., via PhyloPhlAn, UBCG).
• Reagent: PhyloPhlAn software, UBCG pipeline, IQ-TREE/RAxML.

Step 3: Protein Homolog Clustering.

• Perform an all-vs-all BLASTP of predicted proteomes. Cluster proteins into homologous groups (HGs) using tools like OrthoFinder or USEARCH.
• Reagent: OrthoFinder, USEARCH/CLUSTER, MMseqs2.

Step 4: Multiple Sequence Alignment & CSI Identification.

• Align sequences within each HG using MAFFT or MUSCLE.
• Manually inspect alignments for conserved insertions/deletions unique to the hypothesized in-group clade.
• Reagent: MAFFT, MUSCLE, AliView.

Step 5: Pan-Genome Analysis for CSP Discovery.

• Analyze the distribution profile of all HGs across the dataset. Identify HGs present in 100% of in-group genomes and 0% of out-group genomes.
• Functionally annotate these unique HGs (CSPs) using InterProScan and CDD.
• Reagent: Roary/PanX, EggNOG-mapper, InterProScan.

Step 6: Validation and Specificity Testing.

• Test the discovered CSIs/CSPs against a wider, more diverse genomic database (e.g., non-redundant NCBI database) to confirm clade specificity.

4. Visualization of Workflow

Diagram 1: CSI/CSP Discovery and Validation Workflow (100 chars)

5. The Scientist's Toolkit: Key Research Reagents & Materials

Table 2: Essential Reagents and Tools for CSI/CSP Research

Item	Function/Description
GTDB-Tk Toolkit	Standardized taxonomic classification and genome database.
OrthoFinder	Accurately infers orthologous groups from proteomes.
MAFFT Software	Creates high-quality multiple sequence alignments.
AliView	Rapid manual visualization and editing of alignments.
Roary	Rapid large-scale prokaryote pan-genome analysis.
InterProScan	Integrates protein signature databases for functional annotation.
High-Performance Computing (HPC) Cluster	Essential for processing large-scale genomic data.

6. Application in Marinisomatota: Example Findings

Table 3: Hypothetical CSI/CSP Findings in Marinisomatota Phylogenomics

Proposed Clade (Rank)	CSI Example (Protein, Position)	CSP Example (Gene ID/Name)	Putative Functional Link
Novel Family Marinisomataceae	2-aa insert in RNA polymerase beta' subunit (RpoC)	Unique ABC transporter permease (Msm_01234)	Potential adaptation to marine osmolarity.
Novel Order Marinisomatales	5-aa deletion in DNA gyrase B (GyrB)	Novel tetratricopeptide repeat (TPR) domain protein	Possible protein-protein interaction specialization.
Phylum Marinisomatota	N/A (multiple smaller CSIs)	3 unique, conserved proteins of unknown function (CSP1-3)	Defining molecular synapomorphies for the phylum.

7. Implications for Drug Development

The identification of CSPs, in particular, offers high-value targets. As proteins unique to a pathogenic or industrially relevant Marinisomatota clade, they present opportunities for highly specific:

• Diagnostics: PCR primers or antibody probes targeting CSP gene sequences.
• Therapeutics: Inhibition of CSPs essential for viability in pathogenic clades, minimizing off-target effects on human microbiome.
• Enzymatic Discovery: Novel CSPs may represent enzymes for specialized secondary metabolite biosynthesis (e.g., novel antibiotics).

This whitepaper details the critical methodologies and analytical frameworks for temporal calibration within a broader thesis dedicated to resolving the deep evolutionary history of the candidate phylum Marinisomatota (also known as CPR group Marinisomatota). Accurate bacterial dating is paramount for placing the acquisition of key metabolic pathways, symbioses, and diversification events in geologic time, thereby transforming a phylogenetic tree into a time-scaled evolutionary narrative essential for understanding this enigmatic group's role in global biogeochemical cycles and its potential interactions with other life forms.

Core Principles and Challenges

Temporal calibration, or "molecular dating," infers the timescale of evolutionary history using genetic data and fossil or geological evidence. For bacteria like Marinisomatota, which lack a conventional fossil record, this presents unique challenges.

Key Challenges:

• Lack of Paleontological Proxies: Direct fossil evidence is extremely rare.
• Horizontal Gene Transfer (HGT): Pervasive HGT can obscure vertical phylogenetic signals used for dating.
• Rate Heterogeneity: Molecular evolutionary rates vary across lineages and over time, complicating clock models.
• Ancient Divergences: Deep nodes are sensitive to prior assumptions and model misspecification.

Opportunities:

• Genome-Scale Phylogenomics: Dense sampling of genes reduces stochastic error and helps identify core genes resistant to HGT.
• Geological Event Calibration: Using the ages of vicariance events (e.g., ocean basin formation, host lineage divergence for symbionts).
• Relaxed Clock Models: Bayesian methods (e.g., implemented in BEAST2, MCMCTree) account for rate variation among branches.
• Archaeal/ Eukaryotic Proxies: Calibrating based on co-evolution with datable hosts or environments.

Table 1: Common Geological and Biological Calibration Points for Bacterial Dating

Calibration Type	Example Event	Applicable to Marinisomatota?	Justification & Uncertainty
Great Oxidation Event (GOE)	Rise of atmospheric O~2~ ~2.4 Gya	Indirectly, for aerobic lineages	Provides a maximum age for oxygen-dependent metabolisms. Broad window (~2.3-2.5 Gya).
Host Divergence	Divergence of a eukaryotic host lineage	If symbiotic lifestyle is proven	Assumes co-divergence; risk of host-switching. Age derived from host fossil record.
Biomarker Fossils	Steranes from eukaryotes ~1.6 Gya	Indirectly, for associated communities	Provides minimum age for eukaryotic interaction.
Geographic Isolation	Closure of Isthmus of Panama ~3 Mya	For marine taxa with divided populations	Requires robust population genetic study to link vicariance to speciation.
Ancient Gene Duplication	Paralogue roots within gene families	Yes, for core metabolic genes	Requires clear orthology/paralogy delineation. Provides a minimum age.

Table 2: Comparison of Molecular Clock Software and Models

Software Package	Core Method	Key Feature	Best Suited For
BEAST2	Bayesian MCMC	Integrated relaxed clocks, flexible calibration priors (e.g., lognormal), user-friendly GUI (BEAUti).	Complex datasets with multiple calibrations, rate heterogeneity.
MCMCTree (PAML)	Bayesian MCMC	Efficient approximate likelihood, handles very large phylogenies.	Deep-time phylogenies with genome-scale data.
r8s	Penalized Likelihood	Fast, less computationally intensive than Bayesian methods.	Preliminary analyses, large trees under smooth rate variation.
TreePL	Penalized Likelihood	Highly optimized, can handle very large trees.	Phylogenies with 10,000+ tips where Bayesian is infeasible.

Detailed Experimental Protocol for aMarinisomatota-Focused Dating Analysis

Protocol: Time-Calibrated Phylogenomic Analysis Using BEAST2

Objective: To infer a time-scaled phylogeny for Marinisomatota and related Candidate Phyla Radiation (CPR) groups.

Step 1: Dataset Curation

• Genome Collection: Assemble a genomic dataset including publicly available Marinisomatota genomes, representative genomes from other CPR phyla, and outgroup taxa from well-established bacterial phyla (e.g., Terrabacteria).
• Core Genome Identification: Use OrthoFinder or similar to identify single-copy orthologous genes (SCGs) present in >90% of taxa.
• Alignment and Filtering: Align each SCG with MAFFT. Trim poorly aligned regions using trimAl (-automated1). Concatenate alignments into a supermatrix. Generate a partition file defining each gene.

Step 2: Substitution Model and Clock Model Selection

• Best-Fit Model: Determine the best-fit substitution model for each partition using ModelTest-NG or PartitionFinder2.
• Clock Testing: Perform a preliminary Bayesian analysis (without dates) in BEAST2 with a RandomLocalClock or RelaxedClockLogNormal model. Use Tracer to assess clock-likelihood and coefficient of variation—high variation supports a relaxed clock.

Step 3: Calibration Strategy Implementation

• Primary Calibration (Example): If any Marinisomatota lineage is inferred as an obligate symbiont of a marine protist with a fossil first appearance (e.g., 400 Mya), apply a lognormal prior to that node (mean in real space=400, offset=0, log StDev=0.5-1.0) to reflect uncertainty.
• Secondary Calibration: Use a previously published, well-justified age estimate for the divergence of CPR from other Bacteria (e.g., a conservative mean ~2.5 Gya) as a secondary, soft-bound calibration with a broad credible interval.

Step 4: BEAST2 Analysis Execution

• XML Configuration: Use BEAUti to set up the analysis: load alignment/partitions, select site and clock models (Relaxed Clock Log Normal), define tree prior (e.g., Birth-Death), and apply calibration priors on relevant nodes in the tree.
• MCMC Run: Run multiple independent MCMC chains for at least 100 million generations, sampling every 10,000. Ensure adequate chain mixing and ESS values >200 for all parameters (checked in Tracer).
• Post-Processing: Use LogCombiner to merge runs (discarding burn-in). Generate the maximum clade credibility time-tree with TreeAnnotator.

Step 5: Validation and Interpretation

• Perform a cross-validation analysis by sequentially removing one calibration point to assess its influence on node age estimates.
• Compare results with an alternative method (e.g., MCMCTree) to check for robustness.

Mandatory Visualizations

Bacterial Dating Workflow

Calibration Source Integration

The Scientist's Toolkit

Table 3: Essential Research Reagent Solutions for Phylogenomic Dating

Item / Software	Function / Purpose	Key Considerations
OrthoFinder	Identifies orthologous gene groups across genomes.	Critical for building a robust, HGT-minimized core genome dataset.
trimAl	Automatically trims spurious sequences/poorly aligned regions.	Improves alignment quality, reducing systematic error in divergence estimates.
PartitionFinder2 / ModelTest-NG	Selects best-fit nucleotide substitution model per partition.	Model accuracy improves branch length estimation, fundamental for dating.
BEAST2 Package	Bayesian evolutionary analysis for molecular dating.	Industry standard; requires careful configuration of priors and models.
Tracer	Diagnoses MCMC chain convergence and ESS.	Essential for validating the statistical reliability of dating results.
FigTree / IcyTree	Visualizes and annotates time-scaled phylogenetic trees.	Enables interpretation and presentation of node ages and credibility intervals.
Lognormal/Uniform Prior Densities (Conceptual)	Define probabilistic distributions for calibration nodes.	Lognormal priors are soft and realistic for most biological calibrations.
High-Performance Computing (HPC) Cluster	Provides computational resources for large phylogenomic analyses.	Dating analyses with genome-scale data are computationally intensive.

This whitepaper details a core methodological component of a broader thesis investigating the evolutionary history of the candidate phylum Marinisomatota. This recently described lineage, prevalent in marine subsurface and sediment niches, presents a unique opportunity to study microbial adaptation to extreme and oligotrophic environments. A central pillar of our phylogenomic research is the identification of genes under positive (diversifying) selection, which provides direct molecular evidence for adaptive evolution. This guide outlines the technical workflow for evolutionary rate analysis, specifically targeting genes that have been instrumental in the colonization and specialization of Marinisomatota within marine ecosystems.

Core Conceptual Framework: Evolutionary Rate Metrics

The detection of positive selection relies on quantifying the ratio (ω) of non-synonymous nucleotide substitutions per non-synonymous site (dN) to synonymous substitutions per synonymous site (dS). A ω > 1 indicates positive selection.

Table 1: Key Evolutionary Rate Metrics and Interpretation

Metric	Calculation	Interpretation	Value Indicating Positive Selection
dN	Non-synonymous substitutions / Non-synonymous sites	Rate of amino acid-changing mutations	N/A
dS	Synonymous substitutions / Synonymous sites	Rate of silent mutations (neutral baseline)	N/A
ω (dN/dS)	dN / dS	Selection pressure on protein	ω > 1

Detailed Experimental Protocol

Prerequisite: Genome and Ortholog Data Collection

• Source Material: High-quality metagenome-assembled genomes (MAGs) and/or isolate genomes of Marinisomatota and related outgroup taxa (e.g., other FCB group members).
• Objective: Construct a robust multiple sequence alignment for each protein-coding gene.

Protocol:

• Gene Prediction & Annotation: Use Prodigal v2.6.3 to predict open reading frames. Annotate functionally using eggNOG-mapper v2.1.12 against the COG and KEGG databases.
• Ortholog Identification: Perform an all-vs-all BLASTP (v2.13.0+) search with an E-value cutoff of 1e-10. Cluster genes into orthologous groups using OrthoFinder v2.5.5 with the MSA option (-M msa).
• Alignment and Cleaning: Align amino acid sequences for each orthogroup using MAFFT v7.505 (--auto). Back-translate to codon-aware nucleotide alignments using PAL2NAL v14. Poorly aligned regions are removed with trimAl v1.4.1 using the -automated1 heuristic.

Phylogenetic Tree Reconstruction

• Objective: Generate a species tree for branch-site model tests. Protocol:

• Concatenate alignments of single-copy universal orthologs (e.g., 120 marker genes).
• Construct a maximum-likelihood phylogeny using IQ-TREE v2.2.2.7 with ModelFinder Plus (-m MFP) and 1000 ultrafast bootstrap replicates (-B 1000).
• Root the tree using the outgroup taxa.

Detection of Positive Selection: Branch-Site Model

• Objective: Test if specific foreground branches (e.g., the stem lineage leading to Marinisomatota) show evidence of positive selection in a subset of sites within a gene.

Protocol (Using CODEML from PAML v4.10.7):

• Prepare Control File: Configure a codeml.ctl file specifying:
  • seqfile = cleaned codon alignment
  • treefile = Newick tree with foreground branch labeled
  • model = 2 (branch-site)
  • NSsites = 2
  • fix_omega = 0 (for alternative model, Alt) and 1 (for null model, Null)
  • omega = 1.5
• Run Models: Execute CODEML twice: once under the Null model (fix_omega = 1) and once under the Alternative model (fix_omega = 0).
• Likelihood Ratio Test (LRT): Extract log-likelihood values (lnL) from both runs. Calculate the test statistic: 2*(lnLAlt - lnLNull). This statistic follows a χ² distribution with 1 degree of freedom. A significant p-value (after correction for multiple testing, e.g., FDR < 0.05) rejects the null hypothesis, indicating positive selection on the foreground branch.
• Identify Sites: Under the significant Alternative model, use the Bayes Empirical Bayes (BEB) analysis to identify codon sites with posterior probability > 0.95 of being under positive selection.

Table 2: Example CODEML Results for a Hypothetical Marinisomatota Gene

Gene ID (Orthogroup)	Null Model lnL	Alt Model lnL	LRT Statistic	p-value (FDR adj.)	BEB Sites (PP>0.95)	Proposed Function
MSOG_00154	-3256.78	-3251.24	11.08	0.0009	12, 45, 178	TonB-dependent transporter
MSOG_00732	-4102.15	-4100.87	2.56	0.109 (ns)	N/A	DNA polymerase III

The Scientist's Toolkit: Research Reagent Solutions

Table 3: Essential Tools and Resources for Evolutionary Rate Analysis

Item	Function/Description	Example Tool/Resource (Version)
Genome Assembly/Prediction	Reconstruct and identify coding sequences from raw sequencing data.	Prodigal (v2.6.3), SPAdes (v3.15.5)
Orthology Inference	Define groups of genes descended from a single gene in the last common ancestor.	OrthoFinder (v2.5.5), Proteinortho (v6.1.2)
Sequence Alignment	Create accurate multiple sequence alignments for phylogenetic analysis.	MAFFT (v7.505), Clustal Omega (v1.2.4)
Phylogenetic Reconstruction	Infer evolutionary relationships among taxa.	IQ-TREE (v2.2.2.7), RAxML-NG (v1.2.0)
Selection Analysis Software	Perform codon-substitution model tests (dN/dS).	PAML/CODEML (v4.10.7), HyPhy (v2.5.52)
Multiple Testing Correction	Adjust p-values to control false discovery rate across many genes.	Benjamini-Hochberg procedure (statsmodels v0.14.0 in Python)
Visualization & Reporting	Visualize phylogenetic trees and generate publication-quality figures.	FigTree (v1.4.4), ggtree (R package, v3.6.2)

Visualization of Workflows

Conclusion

Phylogenomic analysis has fundamentally reshaped our understanding of the Marinisomatota phylum, precisely delineating its evolutionary history and relationships within the bacterial domain. By integrating robust methodological frameworks, overcoming analytical challenges, and employing rigorous validation, researchers can confidently map the genetic innovations that underpin this group's adaptation to marine ecosystems. The future of this field lies in leveraging these high-resolution evolutionary maps to guide functional studies and bioprospecting. The identified biosynthetic gene clusters and unique metabolic pathways, illuminated by their evolutionary context, present a promising frontier for the discovery of novel antimicrobials, enzymes, and bioactive compounds, directly impacting biomedical and clinical research pipelines.