For this study, a comprehensive set of 30 genomes was obtained from the NCBI database, encompassing diverse strains from the four Frankia genera, including both uncultured microsymbionts from Protofrankia and Frankia genera and their closely related cultivated counterparts (Table 1). As per the Minimum Information about a MAG (MIMAG) standards (Bowers et al., 2017), MAGs are classified as ‘high-quality’ if they exhibit >90% completeness and less than 5% contamination. In this study, the completeness of the used MAGs ranges from 82.45 to 97.39%, and contamination levels range from 0.25 to 5.2%. These values indicate that the MAGs meet the criteria necessary for conducting comparative analyses with isolate-assembled genomes (IAGs), thereby substantiating rigorous downstream functional analysis and ensuring the validity of potential conclusions.

Pairwise Mash distances were calculated using Mash v2.1 (Ondov et al., 2016). Digital DNA–DNA hybridization (dDDH) values and confidence intervals were calculated using GGDC version 3.0 (Meier-Kolthoff et al., 2013). Pairwise average nucleotide identity (ANI) were calculated at https://www.ezbiocloud.net (Yoon et al., 2017).

For the presence of insertion sequences, ISfinder (Siguier et al., 2006, 2014) was employed. Pseudogene prediction was conducted utilizing Pseudofinder (Syberg-Olsen et al., 2022), and CRISPR elements were identified using the CRISPR Recognition Tool (Bland et al., 2007). Data concerning plasmid-like and virus-like components within a genome were sourced from IMG-M (Chen et al., 2023), incorporating geNomad for the automated identification of both virus-like and plasmid-like entities. Putative horizontal gene transfer (HGT) events were identified using data from IMG-M (Chen et al., 2023).

Pangenome analysis which permitted to categorize genes into core genes include strict-core (>99 to 100% prevalence) and soft-core genes (>95 to >96% prevalence), and cloud genes, was performed using Panaroo v1.3.4 (Tonkin-Hill et al., 2020) with default settings (sequence identity threshold 95%, protein family sequence identity threshold 70%, length difference threshold 95%). Gene gain and loss events were assessed using Panstripe (Tonkin-Hill et al., 2023) alongside the phylogenetic tree from the core gene alignment constructed using IQ-TREE (Nguyen et al., 2015), along with the Panaroo gene presence-absence matrix.

To investigate the evolution of gene content, ancestral gene numbers and events were estimated using the Dollo parsimony method implemented in COUNT v9.1106 program, which does not require an out-group species for the phylogenetic tree (Csurös, 2010). The final tree was generated using the web application chiplot (Xie et al., 2023).

To explore genomic differences in uncultured and cultured strains, a detailed comparative analysis was conducted using Orthovenn3 using the OrthoFinder algorithm, with an e-value of 1e⁻² and an inflation value of 1.50, was employed for orthologous analyses and GO term enrichment analysis was also performed using Orthovenn3 (Sun et al., 2023). Pairwise non-synonymous and synonymous substitution ratios (dN/dS) for each single-copy orthologous gene (OG) were estimated with the Codeml module of the PAML v4.10.0 package using the following settings (model = 0 and codon frequency = 2) (Yang, 2007).

Functional annotations of uncultured Frankia and Protofrankia, along with their closely related cultured counterparts, were carried out utilizing the Clusters of Orthologous Genes (COG) database (Galperin et al., 2019). The completeness of metabolic pathways was assessed using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (Kanehisa et al., 2002; Touchon and Roche, 2007). The KEGG metabolic pathways were examined through BlastKOALA (Kanehisa et al., 2016), aligning against the prokaryotes database to elucidate pathway functionalities and completeness.

Cell Wall and Capsule, Dormancy and Sporulation, and Regulation and Cell signaling, were obtained using SEED server for Rapid Annotation Subsystem Technology server (RAST) (Overbeek et al., 2014). The heatmaps were created with ComplexHeatMap package (Gu, 2022) on RStudio.

Genomes of uncultured and cultured strains were submitted to the PhenDB web server, ¹ using recommended default settings with a balanced accuracy cutoff of 0.75 and a prediction confidence cutoff of 0.6 (Feldbauer et al., 2015) were employed to predict obligate or facultative intracellular lifestyles. GapMind was used to detect gene pathways for essential carbon sources (Price et al., 2021) and amino acid biosynthesis (Price et al., 2020). Functional traits extraction was mainly performed using MicroTrait, an R package streamlining the extraction of fitness traits from microbial genome sequences (Karaoz and Brodie, 2022). This tool employs profile hidden Markov models (profile-HMM) and logical operations to predict and map the protein family content within genome sequences to fitness traits. The microTrait framework depends on meticulously curated HMMs, known as microtrait-hmms, which capture the sequence diversity of protein families sourced from IMG/M. These HMMs undergo benchmarking using KEGG orthologs for trusted cutoff (TC) scores. Software dependencies include HMMER (Potter et al., 2018) and Prodigal (Hyatt et al., 2010), with additional features utilizing Infernal (Nawrocki and Eddy, 2013), tRNAscan-SE (Schattner et al., 2005) and bedtools (Quinlan, 2014). Essential data dependencies encompass microtrait-hmm (gene level profile-HMM database) and dbCAN (Huang et al., 2018)-HMMdb (domain level profile-HMM database for Carbohydrate-active enzymes). Additionally, secondary metabolite clusters were determined using AntiSMASH 7.0 (Blin et al., 2023).

The 560 single-copy gene phylogenomic trees (Figure 1A) revealed a robust phylogenetic radiation within the Frankiaceae family supported by robust bootstraps and posterior probability values. This analysis clearly delineated the four genera—Pseudofrankia, Protofrankia, Parafrankia, and Frankia (Gtari, 2022). Mash-based analyses, utilizing MinHash sketches and Jaccard indices for shared k-mers, enabled cohesive clustering of the 30 genomes (Supplementary Figure S1), ensuring a robust alignment of type strain species classifications within their respective genera. Significantly, the most elevated similarities observed for Frankia umida Ag45_Mut15^T and Candidatus Frankia nodulisporulans members were 97.15% with AgTrS, 96.80% with AgUmASH1, and 96.82% with AgUmASt1 (Supplementary Figure S1). Digital DNA:DNA hybridization (dDDH) between Candidatus Frankia nodulisporulans members and Ag45_Mut15^T yielded 73.7% ([70.7–76.5]), designating AgTrS, AgUmASH1, and AgUmASt1 as members of Frankia umida species (Chun et al., 2018). Candidatus Protofrankia datiscae Dg1 and Protofrankia symbiont of Coriaria nepalensis Dg1_Cn_nod exhibited genomic similarities of 96.72 and 96.65%, respectively, with Protofrankia coriariae BMG5.1^T (Supplementary Figure S2). dDDH values ranged from 66.5% ([63.5–69.3]) to 67.5% ([64.5–70.3]), classifying them as a borderline taxon. Candidatus Frankia alpina strains AiOr and AvVan displayed values of 95.58 and 95.86%, respectively, with Frankia gtarii strain Agncl-4^T. However, dDDH values of 56.90% ([54.1–59.6%]) for AiOr and 4.50% ([51.8–57.2%]) for AvVan suggest genomic distinctions.

This observation aligns with findings in microbial ecology, where the uncultivability of certain taxa is not solely attributed to evolutionary divergence but may involve a complex interplay of ecological adaptations and genomic features (Stewart, 2012). It challenges the conventional notion that uncultured microbes represent entirely novel lineages, emphasizing the need for a nuanced understanding of microbial diversity and ecology (Connon and Giovannoni, 2002). The potential genomic continuity hints at shared ancestry, possibly driven by environmental factors influencing microbial evolution (Shade et al., 2014).

As stated previously, an obligate intracellular lifestyle may lead to genome reduction, making microorganisms uncultivable in the lab (Casadevall, 2008; Rodríguez-Gijón et al., 2022). Genomic comparisons between cultured and uncultured Frankia strains revealed consistent genome size differences (0.924–3.427 Mb) (Figure 1A), a trend also observed among cultivated strains associated with host plant speciation (Normand et al., 2007; Tisa et al., 2016). Insertion sequences, horizontal gene transfer, viral interactions, plasmid dynamics, and CRISPR systems collectively drive genome evolution (Casadevall, 2008; Sloan and Moran, 2013; Brito, 2021; Arnold et al., 2022). These elements play a role in rearranging genomes, introducing novel genetic material, and the CRISPR systems serve as a defense mechanism against integration (Nussenzweig and Marraffini, 2020). Diverse adaptive strategies were evident in gene density (ratio of the number of genes per number of base pairs) variations, with some uncultured Frankia strains exhibiting high gene densities (Figure 1A), possibly indicating specialized adaptations (Nakabachi et al., 2006; Kuo et al., 2009; Kuo and Ochman, 2010). Pseudogenes identified in the entire set of 30 genomes predominantly displayed features such as incompleteness, frameshifts, multiple issues, internal stop codons, with fewer instances of ambiguous residues (Supplementary Table S1). The heightened presence of pseudogenes in Dg2 implies continuous evolutionary dynamics, likely influenced by specific environmental pressures. Furthermore, the KEGG distribution of pseudogenes primarily associated them with “carbohydrate metabolism,” followed by “Genetic information processing,” and “Unclassified genetic information processing” (Supplementary Figure S3). The examination of gene events in the 30 analyzed Frankiaceae genomes unveiled a predominant sequence of occurrences, primarily involving loss, followed by duplication, and subsequently, gain (Figure 1A). Notably, in the case of uncultured Frankia, the total number of gene gains and losses was comparatively lower than in their cultured counterparts (Figure 2A). Specifically, the counts range from 1,000 to 2,000 for uncultured Frankia, whereas they range from 3,000 to 4,000 for cultured strains. Conversely, in the Protofrankia group, cultured strains exhibit similar or lower counts (1,000–2,000) compared to their uncultured counterparts (1,000–5,000) (Figure 3A).

Gene-gain events are contingent on the rates of horizontal gene transfer and gene duplication, while gene-loss events hinge on the rates of gene inactivation and deletion (Bordenstein and Reznikoff, 2005). The occurrence of various horizontal gene transfers events, regardless of the cultivation status, underscores the inherent ability of Frankia strains, whether cultured or uncultured, to facilitate the exchange of genetic material. This phenomenon is intricately linked to the symbionts’ ability to establish contact not only with conspecific strains but also with other bacterial species, a possibility that extends even within the confines of actinorhizal nodules (Ghodhbane-Gtari et al., 2021).

The pangenome categorizes genes into three groups: core genes, shared universally across all strains (subdivided into core and soft core genes); accessory genes, found in subsets of strains (subdivided into shell and cloud genes); and unique/singleton genes, exclusive to individual strains (Tettelin et al., 2005). Core genes are essential for basic functions (Tettelin et al., 2005; Kim et al., 2020), while accessory and unique genes reflect adaptations to environments, hosts, or lifestyles (Zhang and Sievert, 2014; Nelson and Stegen, 2015), which may influence their propensity for growth under axenic conditions.

Our investigation of 30 Frankiaceae genomes revealed a diverse genomic landscape with 37 universally present core genes, 17 soft core genes in a substantial subset, and 4,559 shell genes (Figure 1B). Additionally, 100,164 cloud genes indicate a dynamic genetic reservoir within this taxonomic group, forming a pangenome of 104,777 genes that offers insights into genetic diversity, adaptability, and evolutionary dynamics.

Analysis of Frankia and Protofrankia revealed distinct genetic landscapes. The pangenome analysis of Frankia strains revealed a total of 19,414 genes (Figure 2B). Among these, 588 core genes were present in 99–100% of the strains, indicating their essential role. The pangenome includes 7,076 shell genes (15–95% presence), contributing to functional diversity, and 11,750 cloud genes (less than 15% presence), highlighting extensive genetic variability.

The pangenome analysis of Protofrankia strains revealed a total of 18,411 genes (Figure 3B). Among these, 194 core genes were present in 99–100% of the strains, indicating their essential role. The pangenome includes 12,949 shell genes (15–95% presence), contributing to functional diversity, and 5,268 cloud genes (less than 15% presence), highlighting extensive genetic variability.

No soft-core genes (95–99% presence) were found in the analyzed Frankia and Protofrankia strains. Soft-core genes are particularly relevant in analyses involving draft genomes as some genes may be missing in draft genomes due to their lower assembly quality than completely assembled genomes (Nelson and Stegen, 2015). The difference in core genes (588 vs. 194) between Frankia and Protofrankia can be attributed to several factors. Sampling bias is a significant factor, as core gene numbers vary with the diversity and quantity of strains analyzed. Protofrankia likely shows a smaller core genome due to fewer sequenced strains and closer genetic relationships compared to Frankia (Nouioui et al., 2014). Methodological differences also play a crucial role, with core genome estimates significantly impacted by genome incompleteness, fragmentation, and contamination, which are more prevalent in metagenome-assembled genomes (MAGs) compared to single isolate assembled genomes (Li and Yin, 2022). Additionally, genomic diversity affects core genome size and composition; Frankia species likely possess a larger core genome due to greater genomic diversification across varied environments, unlike Protofrankia species (Gtari, 2022).

Analyses of all three sets—Frankia, Protofrankia, and Frankiaceae—have revealed open pangenomes (Figures 1C, 2C,D, 3C,D), as indicated by the continuous discovery of new gene families with the addition of new genomes (Costa et al., 2020). The configuration of a pangenome, whether open or closed, is intricately tied to the lifestyle of bacterial species, with sympatric species thriving in interactive, communal environments demonstrating open pangenomes characterized by a continuous influx of new genes through horizontal gene transfer (Medini et al., 2005). Conversely, allopatric species, residing in isolated settings, typically exhibit smaller, closed pangenomes (Medini et al., 2005; Georgiades and Raoult, 2011).

Three distinct types of evolutionary selection processes, based on codon substitution, include positive selection (dN/dS > 1), purifying selection (dN/dS < 1), and neutral selection (dN/dS = 1) (Spielman and Wilke, 2015). The average dN/dS ratios for genes in both Frankia and Protofrankia, whether cultured or uncultured, were consistently found to be less than 1 (Figures 2E, 3E). This pattern suggests that all analyzed genes are subject to purifying selection, signifying the preservation of functional integrity and evolutionary stability across the analyzed strains (Sarkar and Sen, 2022).

The absence of specific genes in uncultivated microsymbionts, as revealed by Orthovenn3 analysis, highlights potential factors contributing to the challenges associated with their cultivability. In cultivated Frankia strains, 521 unique genes, with functions spanning transcription regulation, fatty acid biosynthesis, and defense responses, suggest enhanced capabilities in cellular regulation, metabolic activities, and defense mechanisms (Tables 2, 3; Supplementary Figure S4). The molecular functions of these genes, including oxidoreductase and ATPase activities, underscore their vital roles in cellular processes. Similarly, Protofrankia exhibits 825 exclusive genes, emphasizing its distinct genetic repertoire (Table 3; Supplementary Figure S5). These genes were predominantly associated with transcription regulation and transmembrane transport, indicating potential adaptations to environmental cues and nutrient uptake strategies. The inclusion of functions related to sporulation, phospholipid transport, oxidoreductase, ATPase activities, DNA binding, and metal ion binding in Protofrankia’s exclusive gene set suggests a diverse array of cellular processes. The absence of these genes in uncultivated strains may signify a less ability compared to cultivated strains to regulate gene expression, perform essential metabolic functions, and adapt to diverse environmental conditions, all of which are crucial for successful laboratory cultivation.

Bacteria with unconventional cell wall structures can pose challenges for cultivation (Van Teeseling et al., 2015; Markova, 2020; Megrian et al., 2022). Apart from the “Exoenzymes regulatory protein AepA in lipid-linked oligosaccharide synthesis cluster,” absent in certain uncultured Frankia strains (AiOr, AvVan, AgTrS, and AgUmASH1), and the “Apolipoprotein N-acyltransferase (EC 2.3.1.-) in lipid-linked oligosaccharide synthesis cluster,” lacking in strain AiOr, all components of the “Cell wall and Capsule pathway” were detected in uncultured strains of Frankia and Protofrankia compared to their cultivated counterparts (Supplementary Figures S6, S7). This finding suggests that the uncultured Frankia and Protofrankia strains are capable of generating intact Cell Wall and Capsule structures.

Dynamics of dormancy and sporulation, contributing to our understanding of the cultivability of specific strains (Kell and Young, 2000; Wörmer et al., 2019). Except for the “Spore pigment biosynthetic cluster in Actinomycetes. Polyketide hydroxylase WhiE VIII,” absent in strain Dg2, all components of the “Dormancy and Sporulation, no subcategory” are detectable in uncultured strains of Frankia and Protofrankia compared to their cultured counterparts. This indicates functional integrity within this pathway across the uncultured strains (Supplementary Figures S6, S7).

Dearth of regulatory elements in obligate lifestyle stems from the loss of structural genes during reductive evolution is a consequence of adapting to narrow ecological niches (Wilcox et al., 2003). Regulation and Cell Signaling, both cultured and uncultured Frankia strains share cAMP signaling components and exhibit consistent expression of essential elements (Supplementary Figure S6). The cAMP signaling, may indicate common regulatory features, and disruptions in these components can influence cultivation success (Vartoukian et al., 2010). The Stringent Response, Programmed Cell Death, and Toxin-Antitoxin Systems are prevalent in all analyzed genomes which contribute to bacterial behavior (Barer and Harwood, 1999; Prozorov and Danilenko, 2010). Uncultured Protofrankia, and Frankia strains along with cultured related strains share all “Regulatory and cell Signaling” features (Supplementary Figure S7). Common enzymatic activities in the cAMP pathway and consistent presence of regulatory proteins underscore similarities in regulatory mechanisms across these strains. Shared functionalities in the stringent response, (p)ppGpp metabolism, play a crucial role in bacterial survival under varying conditions (Fagen et al., 2014), and sporulation-related proteins further highlight commonalities in these pathways (Galperin et al., 2012; Browne et al., 2016).

The absence or deregulation of specific central metabolic pathways is known to hinder the cultivability of certain bacterial taxa (Pande and Kost, 2017; Bodor et al., 2020; Zhao et al., 2023). Analyzing metabolic pathway data from uncultured Frankia and Protofrankia strains alongside their closely related cultured counterparts reveals that crucial central metabolic pathways are largely intact across these organisms (Supplementary Figure S8). Key pathways such as glycolysis (Embden-Meyerhof pathway), gluconeogenesis, the citrate cycle (TCA cycle), and the pentose phosphate pathway exhibit high completeness levels. Furthermore, pathways involved in pyruvate oxidation, fatty acid biosynthesis, and various amino acid biosynthesis pathways are also well-represented with substantial completeness (Supplementary Figure S9).

These findings indicate a diverse metabolic repertoire uncultured Frankia and Protofrankia compered to their cultured counterparts, enabling efficient utilization of diverse carbon sources and synthesis of essential biomolecules critical for growth and survival. The presence of these pathways underscores the metabolic flexibility of these bacterial strains, likely facilitating their adaptation to different environmental conditions and ecological niches. However, notable exceptions are observed in the uncultured Protofrankia Dg2 strain, particularly in pathways such as the Entner-Doudoroff pathway and certain aspects of nucleotide sugar biosynthesis, which display partial completeness. These observations suggest specific metabolic adaptations or limitations in these pathways within the Protofrankia Dg2 strain.

In contrast to obligate symbionts and pathogens that rely on hosts for nutrients, saprophytic and facultative symbiotic microorganisms can utilize decaying organic matter (Ray et al., 2009; Salje, 2021). Their efficient decomposition involves rich enzymatic sets for breakdown, membrane transporters, and uptake through central metabolic pathways, providing nutrients and energy for cellular functions (Cao et al., 2023). Comparing cultured and uncultured strains in complex carbohydrate depolymerization and resource acquisition pathways reveals shared features correlating with similar nutritional requirements (Figures 4A, 5A). Both groups exhibit proficiency in breaking down complex carbohydrates (cellulose, chitin, heteromannan, mixed linkage glucan, xylan, and xyloglucan) and degrading simple compounds (proteins, EMP pathway carbohydrates, maltose, galactose, mannose, trehalose, and glycerol). For substrate uptake, both groups show similarities in transporting aromatic acids, monosaccharides, polysaccharides, dicarboxylates, monocarboxylates, free amino acids, inorganic ions, metal ions, amides, ammonium, nucleotides, organophosphorus compounds, sulfur compounds, siderophores, and vitamin B. Both groups display balanced assimilation of C1 compounds via CO2 fixation through the dicarboxylate hydroxybutyrate pathway, methanotrophy utilizing formaldehyde oxidation or assimilation with serine, assimilatory nitrate reduction, nitrogen fixation, ammonia assimilation to aspartate, glutamate, or glutamine, assimilatory sulfate reduction, and phosphatase-mediated phosphorus compound assimilation. They also exhibit the ability to ferment pyruvate to short-chain fatty acids (SCFAs) and use the oxidative pentose phosphate pathway for aerobic respiration, sharing the presence of the electron transport chain (ETC) complex III. In terms of phototrophy, both cultured and uncultured strains produce the accessory carotenoid pigments neurosporene and lycopene. However, specific pathways like the aerobic respiration nonoxidative pentose phosphate pathway and the glyoxylate cycle are found in various cultured strains but are absent in some uncultured strains. Additionally, the electron transport chain ETC complex IV is present in cultured but missing in all uncultured strains (Figures 4A, 5A). These distinctions highlight the enhanced metabolic diversity and adaptability of cultured strains under controlled conditions.

While host tissues can impose chemical and/or physical stresses, notably through the host immune system’s response to microbe-associated molecular patterns (Schwartzman and Ruby, 2016; Bénard et al., 2020), an obligate lifestyle can be regarded as an adaptive response to narrow ecological niches, providing a buffered environment conducive to maintaining homeostasis. This adaptation involves the regulation of pH, temperature, and osmotic stress, enabling the host to efficiently mitigate these environmental challenges.

The comparison of stress tolerance between cultured and uncultured Frankia strains showed extensive shared behaviors (Figure 4A). Both groups excel in scavenging reactive oxygen species (ROS) and maintaining redox balance, indicating strong defenses against oxidative stress (score 1). Cultured strains, however, exhibit enhanced biofilm formation and EPS biosynthesis/export compared to uncultured strains, potentially increasing their ability to form structured communities. Under heat stress, both groups demonstrate robust protein protection and repair mechanisms involving ATP-dependent proteases. Cultured strains also display superior cold-sensing machinery, highlighting their adaptability to temperature fluctuations. In desiccation, osmotic, and salt stress, both groups show similar proficiency in compatible solute synthesis and transport, suggesting shared strategies for osmoregulation. All strains possess redox-sensing systems for pH adaptation, crucial for survival in diverse environmental conditions. Cultured Protofrankia strains exhibit higher enzyme counts in ROS scavenging, suggesting distinct oxidative stress response mechanisms compared to uncultured strains (Figure 5A). While biofilm formation varies, both groups demonstrate comparable abilities in protecting, repairing, and degrading denatured proteins (score 4–5), emphasizing their resilience under stress conditions. Temperature-specific, osmotic, salt stress responses, and pH stress responses show variations among strains but consistent patterns within each group. Importantly, all strains share capabilities in ROS scavenging, redox-sensing, and response to oxygen limitation and envelope stress, including phage resistance, enhancing their resilience in challenging environments.

It is important to note that cultured Protofrankia strains have showcased a notable preference for alkalophilic conditions, setting them apart from other cultured members of the Frankiaceae family (Gtari et al., 2015). This investigation further delineates these distinctions, revealing that Frankia strains generally display a heightened capacity for the consumption of excessive cytoplasmic protons when compared to Protofrankia. Significantly, specific enzyme systems implicated in pH stress responses exhibit variations between the two genera. Frankia strains, in particular, tend to favor misfolded proteins ATP-dependent proteases, decarboxylase, and antiporter systems, in contrast to Protofrankia. These discerning findings underscore the unique adaptive strategies employed by Frankia and Protofrankia in response to pH stress, potentially elucidating the alkalophilic preference observed in Protofrankia.

Secondary metabolites, distinct from primary metabolites, are characterized by their low molecular mass and do not play direct roles in the growth, development, or reproduction of the producing organism. However, they are often produced during the late growth phase of the microorganisms and play crucial roles in adaptation, competition, and defense mechanisms, conferring fitness advantages (O’Brien and Wright, 2011). These compounds are versatile in their biological functions, participating in inter-microbial warfare through processes such as antibiotic biosynthesis and resistance (Perry et al., 2022) and bacteriocin production (Riley and Wertz, 2002). Additionally, they can act as mediators of cross-species mutualism (O’Brien and Wright, 2011), or even participate in both processes, as seen with siderophores (Kramer et al., 2020). It’s crucial to note that obligate symbiotic bacteria, such as those residing within root nodules, may produce a distinct set of secondary metabolites compared to facultative symbtiotic and cultivable counterparts (Pinto-Carbó et al., 2016). This difference arises from their specialized lifestyle, which prioritizes mutualistic interactions within host environments rather than competitive survival in broader ecological settings. As a result, these bacteria may focus on producing metabolites that support their symbiotic relationship with the host plant, such as compounds involved in nitrogen fixation or modulation of plant physiology, rather than defensive or competitive compounds typically found in free-living bacteria (Gtari et al., 2015).

The analysis of biosynthetic gene clusters in both cultivated and uncultivated Frankia strains revealed that both groups exhibit the presence of polyketide synthases (t2PKS, t3PKS, otherKS), non-ribosomal peptide synthetases (NRPS), and genes related to terpene, siderophore, and lassopeptide biosynthesis, suggesting a common potential for secondary metabolite production (Figure 4A). However, cultured strains show higher counts in specific clusters, including t1PKS1, NRPS, and terpene genes, indicating potential metabolic advantages in these pathways. Conversely, the absence of specific clusters, such as lassopeptide and betalactone, in uncultivated strains suggests potential ecological adaptations and differences in secondary metabolite production. For cultivated and uncultivated Protofrankia strains, shared polyketide synthases (t2PKS, otherKS), NRPS, and genes related to terpene and betalactone biosynthesis are observed (Figure 5A). The presence of siderophore biosynthetic genes suggests a common strategy for iron acquisition. Distinctions, such as lassopeptide genes in cultivated strains and unique distributions of t1PKS1, t3PKS, and lanthipeptide genes in specific uncultivated strains, imply adaptations to specific ecological niches.