Environmental temperature is a strong driver of subspecies competition in the Drosophila microbiome

1. Problem Statement

Microbiome research has historically focused on taxonomic composition at the species level, largely overlooking intraspecific diversity (strain-level variation). While 16S rRNA profiling efficiently characterizes species diversity, it lacks the resolution to distinguish between conspecific strains, which can exhibit significant phenotypic and functional differences (e.g., pathogenic vs. commensal).

• The Gap: It is unclear how environmental factors, specifically temperature, drive the adaptation and competition of specific bacterial strains within a host microbiome.
• The Context: Lactiplantibacillus plantarum is a key nutritional symbiont in Drosophila, yet the functional roles of its different subspecies/clades within the fly gut remain uncharacterized.

2. Methodology

The authors employed a multi-faceted approach combining experimental evolution, comparative genomics, metagenomics, and phenotypic assays.

• Experimental Evolution System:

  • Utilized 20 replicate populations of Drosophila simulans derived from a natural Florida population.
  • Populations were maintained for >10 years (approx. 600+ generations) under two novel temperature regimes:
    • Hot: Fluctuating 28°C (day) / 18°C (night).
    • Cold: Fluctuating 20°C (day) / 10°C (night).
  • Additional isolates were sourced from independent experiments in Portugal (hot) and South Africa (constant cold).

• Genomic & Metagenomic Analysis:

  • Isolation & Sequencing: 33 L. plantarum isolates were cultured and sequenced (DNBSEQ).
  • Phylogeny: Constructed Maximum Likelihood trees using core genome SNPs and calculated Average Nucleotide Identity (ANI) to define clades.
  • Pangenomics: Analyzed accessory gene presence/absence and performed KEGG pathway enrichment to identify metabolic differences.
  • Metagenomic Tracking: Used competitive mapping of Pool-Seq reads against clade-specific pangenome graphs to track the relative abundance of clades over time in the host populations.

• Phenotypic Assays:

  • In Vitro Growth: Measured growth curves (growth rate, carrying capacity, inflection time) of 29 isolates in liquid MRS medium at hot temperatures.
  • In Vivo Mono-association: Inoculated axenic D. melanogaster and conventionally reared D. simulans with specific clades (U, H, C) or a mixture.
  • Fitness Metrics: Measured offspring number (eclosion) and developmental time under hot and cold conditions.

3. Key Contributions

• Discovery of Cryptic Diversity: Identified three distinct, sympatric clades of L. plantarum (designated U for unevolved, H for hot-evolved, and C for cold-evolved) co-existing in the Drosophila microbiome.
• Resolution of "Hidden" Adaptation: Demonstrated that environmental temperature drives selection at the subspecies level, leading to the dominance of specific clades that are indistinguishable by standard 16S profiling.
• Functional Decoupling: Showed that the clade dominating a specific temperature regime is not necessarily the one with the highest growth rate in liquid culture, highlighting the complexity of host-microbe interactions.
• Genotype-Dependent Symbiosis: Revealed that the fitness effect of L. plantarum on the host is clade-specific and can be detrimental (reduced fitness) depending on the strain and environmental context.

4. Key Results

A. Genomic Characterization of Clades

• Clade Definition: The three clades (U, H, C) showed ANI values between groups of 98.4–99.17% (below the ~99.5% species threshold), confirming they are distinct subspecies.
• Phylogenetic Origin: Phylogenetic analysis of 92 genomes (including public data) revealed that the three clades diverged before colonizing Drosophila, suggesting independent evolutionary histories rather than recent divergence within the host.
• Metabolic Divergence:
  • Clades H & C: Enriched in genes for chitobiose degradation (via hexosaminidase), allowing them to utilize chitin breakdown products from the fly exoskeleton/peritrophic matrix. This is hypothesized to be a key adaptation to the lab environment.
  • Clade U: Lacks chitobiose utilization but possesses unique genes for nitrate respiration (anaerobic adaptation) and guanylyl molybdenum cofactor synthesis.

B. Dynamics in Experimental Evolution

• Temperature-Driven Selection:
  • Initial State: Clade U was dominant (54.5%) in the ancestral population.
  • Hot Regime: Clade H rapidly increased to dominance; Clade U declined.
  • Cold Regime: Clade C rapidly increased to dominance; Clade U declined.
• Replicate Consistency: The shift was highly reproducible in cold regimes (all 10 replicates). Hot regimes showed more stochasticity and slower convergence, with some replicates temporarily dominated by Clade C before switching to H.

C. Phenotypic Discrepancies (In Vitro vs. In Vivo)

• Liquid Culture: Clade U exhibited the fastest growth rate and highest carrying capacity in MRS broth, outperforming Clade C.
• In Vivo Reality: Despite slower growth in liquid, Clade C dominated the cold fly populations. This suggests that factors other than simple growth rate (e.g., host interaction, death rate, or specific nutrient utilization like chitin) drive fitness in the gut.

D. Host Fitness Effects

• Detrimental Effects: In axenic D. melanogaster, inoculation with Clade C significantly reduced host fitness (fewer offspring, longer developmental time) compared to germ-free controls, particularly in the cold regime.
• Context Dependence: The negative effect was dose-dependent and less pronounced in non-axenic flies (likely due to competition from other microbiome members).
• Conclusion: The "mutualistic" nature of L. plantarum is not universal; it is clade-specific and can be parasitic under certain conditions.

5. Significance

• Redefining Microbiome Complexity: The study proves that intraspecific diversity in the Drosophila microbiome is orders of magnitude higher than previously thought, rivaling human gut strain richness.
• Mechanism of Adaptation: It provides a clear example of how environmental stressors (temperature) select for specific functional traits (e.g., chitin metabolism) at the subspecies level, driving rapid community restructuring.
• Ecological Theory: The findings challenge the assumption that "beneficial" symbionts are universally positive. Instead, the outcome of host-microbe interactions is highly context-dependent (genotype × environment).
• Methodological Impact: Highlights the necessity of moving beyond 16S rRNA sequencing to whole-genome metagenomics and strain-resolved analysis to fully understand microbial ecology and evolution.

In summary, this paper demonstrates that environmental temperature acts as a strong selective filter for bacterial subspecies, reshaping the microbiome through functional differentiation that is invisible to standard taxonomic surveys, with profound implications for host fitness and evolutionary adaptation.