# Stochastic emergence of antifungal resistance via population bottleneck-driven fixation and subsequent introgression in pathogenic fungi

## Metadata

- **SPORE ID**: SPR-2026-2D9D
- **Domaines**: Evolutionary Biology x Mycology
- **Date de generation**: 2026-04-19
- **Panel consensus score**: 6.2/10
- **Novelty score**: 0.85/1.0
- **Panel verdict**: publish_brief

## Abstract

If a pathogenic fungal population experiences a severe bottleneck (effective population size Nₑ = 50-200) in an isolated niche, then genetic drift will fix pre-existing resistance alleles with low selective coefficients (|s| < 0.01) at a probability ≥ 0.05, significantly exceeding expectations under selection alone. Subsequent gene flow via homologous recombination (rate ≥ 10⁻⁸ events/bp/generation) will introgress these alleles into the larger progenitor population, creating phylogenetic reticulations detectable by discordance ≥ 15% between gene trees.

The proposed mechanism involves 4 causal steps: (1) A subpopulation of a pathogenic fungus (e.g., Candida auris) becomes isolated in -> (2) In this bottleneck, genetic drift dominates over selection for alleles with |s|  -> (3) Subsequent contact and gene flow (homologous recombination rate ≥ 10⁻⁸) between  -> (4) Introgression creates phylogenetic network reticulations, with discordance ≥ 15%.

Literature grounding on 6 verified references yields a novelty score of 0.85 (novel). A 3-phase experimental protocol (budget: €20k-80k, timeline: 8-14 months) is proposed, starting with in silico validation. A panel of 5 expert reviewers reached a consensus score of 6.2/10.

## 1. Hypothese et mecanisme propose

### 1.1 Formulation formelle

If a pathogenic fungal population experiences a severe bottleneck (effective population size Nₑ = 50-200) in an isolated niche, then genetic drift will fix pre-existing resistance alleles with low selective coefficients (|s| < 0.01) at a probability ≥ 0.05, significantly exceeding expectations under selection alone. Subsequent gene flow via homologous recombination (rate ≥ 10⁻⁸ events/bp/generation) will introgress these alleles into the larger progenitor population, creating phylogenetic reticulations detectable by discordance ≥ 15% between gene trees.

### 1.2 Variables

**Variables independantes :**

| Variable | Type | Plage | Unite |
|----------|------|-------|-------|
| Bottleneck effective population size (Nₑ) | continuous | 50-200 | number of individuals |
| Pre-bottleneck allele frequency of resistance variants | continuous | 0.001-0.05 | frequency |
| Selection coefficient (s) of resistance allele in absence of drug | continuous | -0.02 to +0.01 | dimensionless |
| Homologous recombination rate | continuous | 10⁻⁹ to 10⁻⁷ | events/base pair/generation |

**Variables dependantes :**

| Variable | Type | Direction attendue | Unite |
|----------|------|-------------------|-------|
| Probability of allele fixation in bottleneck population | continuous | increase | probability (0-1) |
| Time to allele fixation in progenitor population post-introgression | continuous | decrease | generations |
| Phylogenetic discordance between resistance locus and housekeeping genes | continuous | increase | percentage of discordant nodes |
| Minimum inhibitory concentration (MIC) distribution in progenitor population | continuous | non-monotonic | μg/mL |

### 1.3 Chaine causale

1. Step 1: A subpopulation of a pathogenic fungus (e.g., Candida auris) becomes isolated in a niche (e.g., biofilm, specific host tissue) with reduced effective population size (Nₑ = 50-200).
1. Step 2: In this bottleneck, genetic drift dominates over selection for alleles with |s| < 0.01, leading to stochastic fixation of pre-existing resistance variants (e.g., ERG11 mutations) present at low frequency (0.001-0.05).
1. Step 3: Subsequent contact and gene flow (homologous recombination rate ≥ 10⁻⁸) between the bottleneck-derived lineage and the larger progenitor population introgresses the fixed resistance allele.
1. Step 4: Introgression creates phylogenetic network reticulations, with discordance ≥ 15% between trees built from resistance loci versus neutral genomic regions.

**Hypotheses cles :**

- Resistance alleles exist at low frequencies in ancestral populations due to mutation-selection balance or fluctuating selection.
- Bottleneck events are common in clinical/ecological settings (e.g., transmission between hosts, biofilm formation).
- Homologous recombination occurs at sufficient rates to transfer large genomic regions containing resistance alleles.
- The fitness cost of resistance (s in drug-free conditions) is small enough (|s| < 0.02) that drift can overcome selection in small populations.

**Inconnues identifiees :**

- The exact distribution of fitness effects (DFE) for resistance mutations in relevant in vivo conditions.
- The quantitative relationship between Nₑ in clinical niches and measurable census population sizes.
- The rate of recombination between divergent fungal lineages in natural and host environments.

### 1.4 Conditions aux limites

- **Bottleneck size Nₑ > 500** — Above this threshold, selection dominates over drift for alleles with |s| > 0.001, making stochastic fixation negligible
- **Selection coefficient |s| > 0.02 in the absence of drug** — Strongly deleterious alleles are purged even in small populations; strongly beneficial alleles fix deterministically
- **Recombination rate < 10⁻¹⁰ events/bp/generation** — Below this rate, introgression is negligible compared to new mutation as a source of variation
- **Population subdivision with FST < 0.25 between subpopulations** — High gene flow prevents the establishment of significant allele frequency differences via drift

### 1.5 Cadre theorique

Nearly Neutral Theory of Molecular Evolution (Ohta) combined with Founder Effect Speciation (Mayr) and Network Phylogenetics

## 2. Etat de l'art et positionnement

### 2.1 Travaux les plus proches

- **[2010] The existence of species rests on a metastable equilibrium between inbreeding and outbreeding. An essay on the close relationship between speciation, inbreeding and recessive mutations** — [10.1186/1745-6150-6-62](https://doi.org/10.1186/1745-6150-6-62)
  - Similarite: conceptual
  - Difference cle: Joly (2010) proposes a general model of speciation via 'budding' from an ancestral stock, driven by inbreeding. The SPORE hypothesis specifically applies this evolutionary framework to the emergence of antifungal resistance in pathogenic fungi, proposing a detailed mechanism involving genetic drift in small, isolated subpopulations (e.g., biofilms, host niches) and subsequent spread via introgression, which is not addressed in Joly's paper.
- **[2024] Microbial species and intraspecies units exist and are maintained by ecological cohesiveness coupled to high homologous recombination** — [10.1038/s41467-024-53787-0](https://doi.org/10.1038/s41467-024-53787-0)
  - Similarite: mechanistic
  - Difference cle: Conrad et al. (2024) describe how microbial species and intraspecies units are maintained by ecological cohesiveness and recombination. While this supports the idea of gene flow (recombination) shaping populations, it does not specifically model the initial emergence of a trait (like resistance) via a stochastic budding event from a small founder population, which is the core of the SPORE hypothesis.

### 2.2 Base de preuves

- **[2010] The existence of species rests on a metastable equilibrium between inbreeding and outbreeding. An essay on the close relationship between speciation, inbreeding and recessive mutations** — [10.1186/1745-6150-6-62](https://doi.org/10.1186/1745-6150-6-62)
  - Type: analogous | Citations: 22
  - Provides the foundational evolutionary concept of 'budding speciation' where new species arise from small, inbreeding groups derived from an ancestral stock. This is the core model being transferred to mycology.
- **[2024] Microbial species and intraspecies units exist and are maintained by ecological cohesiveness coupled to high homologous recombination** — [10.1038/s41467-024-53787-0](https://doi.org/10.1038/s41467-024-53787-0)
  - Type: indirect | Citations: 14
  - Supports the 'post-divergence gene flow' component of the hypothesis by demonstrating that high homologous recombination is a key mechanism maintaining microbial species and intraspecies units, which aligns with the idea of resistance alleles spreading via introgression.
- **[2022] Genomic Diversity across Candida auris Clinical Isolates Shapes Rapid Development of Antifungal Resistance In Vitro and In Vivo** — [10.1128/mbio.00842-22](https://doi.org/10.1128/mbio.00842-22)
  - Type: indirect | Citations: 64
  - Demonstrates that high genomic diversity in C. auris populations provides a substrate for rapid, multiple-pathway evolution of resistance. This pre-existing variation is a prerequisite for the 'random fixation of pre-existing genetic variants' proposed in the budding model.
- **[2025] Genomic epidemiology and antifungal resistance of emerging Trichophyton indotineae in China** — [10.1080/22221751.2025.2498571](https://doi.org/10.1080/22221751.2025.2498571)
  - Type: tangential | Citations: 8
  - While focused on tracing the spread of an already-resistant lineage, its use of genomic epidemiology on a global collection of isolates is the type of study that could, in principle, be used to detect patterns consistent with a budding origin (e.g., identifying a small, ancestral subpopulation). It does not test the hypothesis directly.
- **[2023] Two chromosome-level genome assemblies of Rhodiola shed new light on genome evolution in rapid radiation and evolution of the biosynthetic pathway of salidroside.** — [10.1111/tpj.16501](https://doi.org/10.1111/tpj.16501)
  - Type: tangential | Citations: 11
  - The paper is cited in the hypothesis as an example of detailed genomic analysis of budding speciation in plants (Rhodiola). It provides the methodological analogy (comparative genomics) but does not discuss budding speciation or resistance evolution.

### 2.3 Contre-preuves et limitations connues

- **[addressable] Cytotoxicity and resistance evolution of a novel antifungal carbon nanoparticle**
  - Shows that yeast populations can rapidly adapt to novel antifungal carbon dots through laboratory evolution experiments. This demonstrates strong, direct selection can drive resistance evolution quickly in a large population, which is the traditional model the SPORE hypothesis seeks to complement/reframe.
- **[minor] Genomic Diversity across Candida auris Clinical Isolates Shapes Rapid Development of Antifungal Resistance In Vitro and In Vivo**
  - Demonstrates rapid evolution of resistance under strong antifungal selection in experimental settings, supporting a classic selective sweep model. This does not contradict the budding hypothesis but shows an alternative, likely common, pathway.

### 2.4 Evaluation de nouveaute

- **Score**: 0.85/1.0
- **Verdict**: novel

## 3. Predictions falsifiables

| # | Prediction | Borne quantitative | Methode | H0 | Test statistique |
|---|-----------|-------------------|---------|-----|-----------------|
| 1 | In experimental bottlenecks of Nₑ = 100, the fixation probability of a resistance allele with s = -0.005 will be ≥ 0.05, exceeding the probability under strong selection (s > 0.1) in large populations. | Fixation probability 0.05-0.10 vs. <0.001 under deterministic selection | Fluctuation assay with whole-population sequencing across 100 replicate bottleneck lineages | H0: Fixation probability in bottlenecks is ≤ fixation probability under selection in large populations (s > 0.1). | One-tailed Fisher's exact test, alpha=0.05, power=0.8 (requires n=100 replicates per condition) |
| 2 | In phylogenomic analyses of resistant clinical isolates, gene trees containing resistance alleles will show ≥ 15% topological discordance with species trees based on 100 conserved housekeeping genes. | 15-40% discordance in nodes with bootstrap support >70% | Maximum likelihood phylogenies from whole-genome sequences, measured by Robinson-Foulds distance | H0: Topological discordance between resistance gene trees and species trees is ≤ 5%. | Permutation test of tree distances with 1000 replicates, alpha=0.01 |
| 3 | Following experimental introgression, the time to reach 50% frequency of a resistance allele in a large recipient population (Nₑ = 10⁴) will be reduced by ≥ 30% compared to de novo mutation under identical selection pressure (s = 0.05). | Reduction from 1000 to ≤700 generations | Long-term evolution experiment with periodic whole-genome sequencing and recombination quantification via SNP linkage decay | H0: Time to 50% frequency via introgression is ≥ time via de novo mutation. | Log-rank test for frequency trajectories, alpha=0.05 |

## 4. Protocole experimental

**Timeline globale**: 8-14 months
**Budget global**: €20k-80k

### 4.1 Phase 1 — In Silico Validation

**Objectif**: Determine if genetic drift in simulated bottlenecks (Nₑ=50-200) can fix low-frequency resistance alleles (|s|<0.01) with probability ≥0.05, and if resulting phylogenetic patterns match clinical isolate data.

**Methodologie**: 1. Wright-Fisher forward simulations using SLiM v4.2 to model bottleneck dynamics. Parameters: Nₑ=50-200, s=[-0.02, +0.01], initial allele frequency=0.001-0.05. 1000 replicates per parameter combination. 2. Analysis of publicly available Candida auris/Trichophyton indotineae genomes from NCBI SRA (BioProjects PRJNA485481, PRJNA949202) using IQ-TREE2 for phylogenies and PhyloNet for network analysis. Calculate Robinson-Foulds distances between gene trees (ERG11, FKS1) and species trees (100 housekeeping genes).

- Cout: €500-2000 (cloud computing + software licenses)
- Duree: 4-6 weeks
- Equipement: High-performance computing cluster or cloud instance (AWS EC2 r6i.xlarge)
- Logiciels: SLiM v4.2, IQ-TREE2, PhyloNet, R v4.3 with ape, phytools, ggtree packages, bcftools/samtools

**Criteres de succes :**

- Fixation probability for s=-0.005, Nₑ=100: ≥0.05 (vs. <0.001 for s>0.1 in large populations)
- Phylogenetic discordance in clinical isolates: ≥15% nodes discordant (bootstrap>70%)

- **GO**: Both criteria met: simulations show P(fixation)≥0.05 AND clinical data show ≥15% discordance
- **NO-GO**: Simulations show P(fixation)<0.01 OR clinical discordance <5% with high confidence (p<0.01)
- **PIVOT**: Only one criterion met: reconsider parameter ranges or focus on specific fungal species showing signal

### 4.2 Phase 2 — Minimal Experimental Validation

**Objectif**: Empirically test if bottleneck of Nₑ~100 in Candida auris leads to stochastic fixation of a marked resistance allele with known small fitness cost.

**Methodologie**: 1. Engineer a C. auris strain (B8441 background) with a neutral marker (mCherry) linked to a resistance allele (fluconazole resistance via ERG11^Y132F) with measured s=-0.005 in drug-free medium. 2. Perform 100 replicate bottleneck experiments: inoculate 100 cells from a large population (initial allele frequency=0.01) into 100 separate wells, grow to saturation (~20 generations), sequence entire population (Illumina MiSeq 2x150bp) to determine fixation. 3. Control: large population (Nₑ>10⁴) replicates under identical conditions.

- Cout: €8k-15k (strain engineering, sequencing, consumables)
- Duree: 2-3 months
- Equipement: Biosafety Level 2 lab, Microplate reader (Tecan Spark), Illumina MiSeq, Microbial incubator, Flow cytometer (for mCherry sorting)
- Logiciels: bcl2fastq, breseq for variant calling, R for statistical analysis

**Criteres de succes :**

- Fixation probability in bottlenecks (Nₑ=100): ≥0.05
- Difference from large population controls: Odds ratio ≥5 (Fisher's exact test p<0.05)

- **GO**: Bottleneck fixation probability ≥0.05 AND significantly higher than large population controls (p<0.05)
- **NO-GO**: No fixation events in bottlenecks OR no difference from controls (p>0.1)
- **PIVOT**: Fixation occurs but at lower probability (0.01-0.04): adjust bottleneck size or test allele with different s

### 4.3 Phase 3 — Full Experimental Protocol

**Objectif**: Validate complete causal chain: bottleneck fixation → introgression via recombination → spread in recipient population with phylogenetic discordance.

**Methodologie**: 1. Evolve 10 independent bottleneck lineages from Phase 2 where resistance fixed. 2. Co-culture each bottleneck lineage with a genetically distinct, susceptible progenitor strain (marked with GFP) under conditions promoting recombination (biofilm, host-mimicking media). Measure recombination rates via loss of marker linkage using whole-genome sequencing (PacBio HiFi). 3. Track introgression in large recipient populations (Nₑ=10⁴) under fluconazole selection (0.5×MIC). Sample every 50 generations for 500 generations, sequence pools to measure allele frequency trajectories. 4. Construct phylogenies from evolved isolates (10× coverage Illumina) for resistance loci vs. housekeeping genes.

- Cout: €50k-100k+ (long-term culturing, extensive sequencing, personnel)
- Duree: 9-12 months
- Equipement: PacBio Sequel IIe for recombination detection, Illumina NextSeq 2000 for time-series sequencing, Continuous culture devices (chemostats or morbidostats), Confocal microscope for biofilm imaging
- Logiciels: pbmm2, pbsv for PacBio analysis, GATK for variant calling, ClonalFrameML for recombination detection, BEAST2 for phylogenetic dating

**Criteres de succes :**

- Time to 50% frequency with introgression: ≥30% reduction vs. de novo mutation control
- Phylogenetic discordance in evolved populations: ≥15% discordant nodes (bootstrap>70%)
- Recombination rate: ≥10⁻⁸ events/bp/generation

- **GO**: All three criteria met: faster spread via introgression, measurable recombination, and phylogenetic discordance
- **NO-GO**: Introgression provides no speed advantage (≤10% reduction) OR recombination rate <10⁻¹⁰
- **PIVOT**: Partial success: recombination occurs but no discordance—focus on population genetic signatures rather than phylogenetics

### 4.4 Quick start : comment demarrer aujourd'hui

- **Peut demarrer maintenant**: Oui
- **Premiere action**: Download Candida auris genome data from NCBI SRA (PRJNA485481) using sratoolkit and run IQ-TREE2 on ERG11 gene sequences to build initial gene tree.
- **Outils**: sratoolkit, fastq-dump, Trimmomatic, SPAdes, prokka, IQ-TREE2
- **Donnees ouvertes**: NCBI SRA (PRJNA485481, PRJNA949202), Candida Genome Database, MycoCosm

## 5. Analyse d'impact

### 5.1 Impact scientifique

Novelty score: 0.85/1.0 (novel)

### 5.2 Applications industrielles et marche

**Score industriel**: 4.5/10

**Forces :**
- Adresse un besoin critique non satisfait dans la lutte antifongique : la prediction de l'emergence de resistances 'stochastiques' echappant aux modeles deterministes, visant directement les departements R&D des grands groupes pharmaceutiques (Pfizer, Gilead, MSD) et les agences de sante publique (CDC, ECDC).
- Avantage competitif potentiel dans la surveillance epidemiologique : une plateforme predictive basee sur ce mecanisme pourrait identifier les souches a 'haut risque d'introgression' avant leur emergence clinique, offrant un delai d'alerte de 6-18 mois sur les methodes actuelles de sequencage reactif.

**Faiblesses :**
- Marche adressable trop niche et difficile a monetiser directement. Le client final (hopital) ne paie pas pour des modeles evolutifs, et le transfert vers un kit diagnostique ou un logiciel SaaS rencontre une competition farouche d'acteurs etablis (QIAGEN, bioMerieux) avec des barrieres reglementaires (FDA/CE-IVD) insurmontables a court terme.
- Le TRL est extremement bas (2-3). La voie de commercialisation la plus realiste est la vente de service de R&D ou de licence de propriete intellectuelle a un grand pharma, un modele a faible marge et a cycle de vente long (>3 ans), avec un ROI incertain face aux priorites internes des acheteurs potentiels.

**Recommandation**: Reorienter l'hypothese vers un service de recherche sous contrat (CRO) pour les departements de resistance antifongique des grands pharmas, en utilisant le protocole valide comme preuve de concept. Parallellement, developper un algorithme predictif minimal (MVP) et le proposer en acces libre pour etablir un standard academique et capturer l'expertise, en visant a long terme une licence exclusive pour un developpeur de logiciel de surveillance epidemiologique (ex : BioNumerics).

### 5.3 Opportunites de financement

| Programme | Agence | Fit | Budget type | Taux succes |
|-----------|--------|-----|-------------|-------------|
| ERC Starting Grant / ERC Proof of Concept | European Research Council | 0.7 | €1.5M pour StG ; €150k pour PoC | ~13% (StG) |
| ANR AAPG Générique - Jeunes Chercheuses et Jeunes Chercheurs (JCJC) | Agence Nationale de la Recherche (France) | 0.85 | €150k - €300k | ~15-20% |
| EMBO Installation Grants / EMBO Scientific Exchange Grants | European Molecular Biology Organization | 0.65 | €50k - €150k (Installation) ; montants plus faibles (Exchange) | Variable, compétitif |

- **ERC Starting Grant / ERC Proof of Concept** (European Research Council): L'hypothèse est suffisamment risquée, fondamentale et novatrice pour un ERC. Le cœur du projet (dérive génétique, fixation, introgression) est de la science de base d'excellence. Une Starting Grant pourrait financer l'ensemble du protocole étendu. Un Proof of Concept serait adapté si un chercheur ERC en poste voulait développer cet axe.
- **ANR AAPG Générique - Jeunes Chercheuses et Jeunes Chercheurs (JCJC)** (Agence Nationale de la Recherche (France)): Fit parfait pour le budget et la maturité. L'appel JCJC vise explicitement des projets exploratoires et risqués portés par des chercheurs en début de carrière. Le protocole proposé (20k-80k) correspond à une sous-partie d'un projet JCJC, qui permettrait de financer intégralement les phases 2 et 3, avec des moyens pour du séquençage et du personnel.
- **EMBO Installation Grants / EMBO Scientific Exchange Grants** (European Molecular Biology Organization): L'EMBO finance l'excellence en biologie moléculaire, y compris la biologie évolutive et fongique. Un Installation Grant (pour chercheurs revenant en Europe) pourrait couvrir le budget de démarrage. Un Exchange Grant serait idéal pour initier une collaboration avec un laboratoire expert en génomique des populations de champignons pathogènes, étape clé pour bâtir le futur consortium.

**Recommandation financement**: Cibler d'abord un financement de démarrage pour atteindre un TRL 3-4 avec des données expérimentales solides sur un modèle fongique. En parallèle, constituer un consortium exploratoire avec un modélisateur, un généticien des populations fongiques et un clinicien/microbiologiste hospitalier. Présenter ensuite le projet comme une recherche fondamentale avec un potentiel de transformation pour la surveillance et la gestion de la résistance.

## 6. Panel Review Summary

| Reviewer | Score | Verdict | Point cle |
|----------|-------|---------|-----------|
| methodologist | 8.2/10 | accept | Le protocole adopte une approche progressive et prudente, avec une phase in silico pour valider les prédictions théoriques avant de s'engager dans des expériences coûteuses et longues. Ce design séquentiel avec des critères GO/NO-GO explicites est une excellente pratique méthodologique. |
| domain_expert | 7.0/10 | weak_accept | L'hypothèse intègre avec succès trois cadres théoriques bien établis mais souvent distincts — la théorie quasi neutre (Ohta), la dynamique de l'effet fondateur et des goulots d'étranglement (Mayr), et la phylogénétique de réseau — en un modèle mécaniste cohérent pour un problème appliqué pressant (la résistance aux antifongiques). Cette synthèse est intellectuellement convaincante. |
| contrarian | 4.5/10 | weak_reject | L'hypothèse établit élégamment un pont entre la théorie de la génétique des populations (dérive dans les goulots d'étranglement) et un problème clinique pressant (résistance aux antifongiques), en proposant une voie non adaptative souvent négligée. |
| industrialist | 4.5/10 | weak_reject | Adresse un besoin critique non satisfait dans la lutte antifongique : la prediction de l'emergence de resistances 'stochastiques' echappant aux modeles deterministes, visant directement les departements R&D des grands groupes pharmaceutiques (Pfizer, Gilead, MSD) et les agences de sante publique (CDC, ECDC). |
| funding_strategist | 6.5/10 | weak_accept | Hypothèse mécanistique originale qui relie la génétique des populations, l'écologie évolutive et la santé publique, un angle très apprécié dans les appels interdisciplinaires. |

**Consensus score**: 6.2/10
**Verdict final**: publish_brief

### 6.1 Consensus

- L'hypothèse propose une synthèse intellectuelle originale et valable entre la théorie de la dérive génétique, les effets de goulot d'étranglement et l'évolution de la résistance antifongique, méritant une investigation plus poussée.
- Le protocole expérimental proposé est bien structuré, avec une approche séquentielle et des critères GO/NO-GO clairs, représentant une bonne pratique méthodologique.
- Les préoccupations majeures portent sur la plausibilité biologique de deux paramètres clés : le taux de recombinaison supposé (≥ 10⁻⁸) et la distribution des effets fitness (|s| < 0.01) pour les allèles de résistance. Ces points sont identifiés comme des risques potentiels de rupture de la chaîne causale.

### 6.2 Points de desaccord

- L'évaluation de la maturité et du potentiel de l'hypothèse : le méthodologiste et l'expert du domaine recommandent l'acceptation/faible acceptation pour tester les prédictions, tandis que le contrarien et l'industriel recommandent le rejet, jugeant les étapes clés (introgression) peu plausibles ou non commercialisables.
- L'interprétation des faiblesses : pour certains, ce sont des points à affiner dans le protocole ; pour le contrarien, ce sont des « raisons d'échec » (fail reasons) qui invalident le mécanisme central sans validation expérimentale préalable.
- La finalité : vision académique de test d'une hypothèse fondamentale vs. attente d'une application translationnelle ou commerciale immédiate.

### 6.3 Critical path

La validation expérimentale directe de l'étape d'introgression par recombinaison homologue. Le succès ou l'échec de l'hypothèse repose entièrement sur la démonstration que cette étape est biologiquement plausible et mesurable dans des systèmes fongiques pathogènes pertinents, avec des taux de recombinaison et des coûts fitness réalistes.

**Recommandation finale**: Le panel reconnaît la créativité et la rigueur méthodologique de la proposition, mais le scepticisme sur la plausibilité biologique du mécanisme central, particulièrement l'étape d'introgression, est substantiel et bien fondé. À cette deuxième itération, sans preuve préliminaire solide pour étayer les paramètres critiques (DFE, taux de recombinaison), l'hypothèse ne remplit pas le seuil de confiance nécessaire pour une publication en l'état. Il est recommandé d'abandonner cette formulation spécifique et de recentrer les efforts sur la collecte de données empiriques concernant les effets fitness des allèles de résistance et les taux de recombinaison dans des niches à goulot d'étranglement, avant de proposer un nouveau modèle mécanistique.

## 7. Gap Manifest residuel

### 7.1 Data gaps

- Lack of direct empirical studies in mycology that test for the signature of a 'budding' event (e.g., severe bottleneck, drift, then introgression) in the evolutionary history of a known antifungal resistance trait. All provided papers on fungal resistance describe patterns consistent with selection or spread, not the stochastic origin proposed.
- No population genetic models in the provided literature formally adapt the plant-based 'budding speciation' model for fungal life histories (e.g., haploidy, clonality, mixed mating systems). The competence gap remains open.

### 7.2 Competence gaps

- Phase 1: Population genetics simulation
- Phase 1: Phylogenomic analysis
- Phase 1: Bioinformatics scripting (Python/R)
- Phase 2: Fungal genetics (CRISPR-Cas9 editing)
- Phase 2: Microbial evolution experiment design
- Phase 2: NGS library prep and analysis
- Phase 3: Long-term evolution experiments
- Phase 3: Recombination rate quantification
- Phase 3: Population genomic analysis
- Phase 3: Phylogenetic network modeling

### 7.3 Epistemic gaps

- The exact distribution of fitness effects (DFE) for resistance mutations in relevant in vivo conditions.
- The quantitative relationship between Nₑ in clinical niches and measurable census population sizes.
- The rate of recombination between divergent fungal lineages in natural and host environments.

## References

[1] E. Joly (2010). *The existence of species rests on a metastable equilibrium between inbreeding and outbreeding. An essay on the close relationship between speciation, inbreeding and recessive mutations*. DOI: [10.1186/1745-6150-6-62](https://doi.org/10.1186/1745-6150-6-62)
[2] R. Conrad, Catherine E. Brink, et al. (2024). *Microbial species and intraspecies units exist and are maintained by ecological cohesiveness coupled to high homologous recombination*. DOI: [10.1038/s41467-024-53787-0](https://doi.org/10.1038/s41467-024-53787-0)
[3] Laura S. Burrack, R. Todd, et al. (2022). *Genomic Diversity across Candida auris Clinical Isolates Shapes Rapid Development of Antifungal Resistance In Vitro and In Vivo*. DOI: [10.1128/mbio.00842-22](https://doi.org/10.1128/mbio.00842-22)
[4] Wenting Xie, Guanzhao Liang, et al. (2025). *Genomic epidemiology and antifungal resistance of emerging Trichophyton indotineae in China*. DOI: [10.1080/22221751.2025.2498571](https://doi.org/10.1080/22221751.2025.2498571)
[5] Dan Zhang, Xiao-Ying Liu, et al. (2023). *Two chromosome-level genome assemblies of Rhodiola shed new light on genome evolution in rapid radiation and evolution of the biosynthetic pathway of salidroside.*. DOI: [10.1111/tpj.16501](https://doi.org/10.1111/tpj.16501)
[6] Sijan Poudel Sharma, Suraj Paudyal, et al. (?). *Cytotoxicity and resistance evolution of a novel antifungal carbon nanoparticle*. DOI: [10.1101/2024.02.11.579833](https://doi.org/10.1101/2024.02.11.579833)

## Annexes

### A. Detailed Reviewer Reports

#### Methodologist

- **Score**: 8.2/10 | **Verdict**: accept | **Confidence**: 0.88
- **Strengths**: Le protocole adopte une approche progressive et prudente, avec une phase in silico pour valider les prédictions théoriques avant de s'engager dans des expériences coûteuses et longues. Ce design séquentiel avec des critères GO/NO-GO explicites est une excellente pratique méthodologique.; La conception expérimentale intègre des contrôles appropriés à chaque phase (ex: populations larges en Phase 2, contrôle par mutation de novo en Phase 3) et utilise des réplicats biologiques en nombre élevé (100 lignées de goulot d'étranglement), ce qui renforce la validité interne et la robustesse des conclusions.; L'hypothèse est formulée de manière précise et falsifiable, avec des prédictions quantitatives claires et des méthodes de mesure adéquates (probabilité de fixation, distances de Robinson-Foulds, suivi de fréquence allélique). La distinction entre les effets de la dérive et de la sélection est bien articulée.
- **Weaknesses**: La puissance statistique pour la Phase 2, bien que raisonnable (100 réplicats), pourrait être marginale pour détecter une probabilité de fixation de 0.05. Avec P=0.05, on s'attend à ~5 événements de fixation sur 100. Une fluctuation aléatoire pourrait facilement donner 2-3 événements, menant à un faux NO-GO. Une analyse de puissance préalable plus formelle est nécessaire.; Le risque de biais de mesure dans la Phase 3 est sous-estimé. Le suivi de la fréquence allélique par séquençage de pools tous les 50 générations peut manquer des dynamiques transitoires ou être biaisé par des différences d'efficacité d'amplification. Un contrôle par séquençage d'isolats individuels à des points clés serait un complément utile.; Le protocole repose fortement sur la mesure précise du coefficient de sélection (s = -0.005) pour l'allèle de résistance marqué. Une erreur dans cette mesure (risque noté [medium]) invaliderait l'interprétation centrale de la Phase 2. La méthode pour mesurer 's' n'est pas décrite et devrait être détaillée, avec une estimation de l'intervalle de confiance.
- **Questions**: Pour la Phase 2, quel est le plan d'analyse statistique exact pour comparer la probabilité de fixation entre les goulots d'étranglement et les populations larges ? Un test exact de Fisher sur une table de contingence 2x2 (fixé vs non-fixé ; goulot vs large) est-il prévu, et sa puissance a-t-elle été calculée compte tenu du faible nombre d'événements attendus ?; Dans la Phase 3, comment allez-vous isoler et quantifier spécifiquement l'effet de l'introgression par recombinaison homologue, par rapport à d'autres mécanismes potentiels de transfert génétique (ex: conversion génique, transfert horizontal de chromosomes) ou à une sélection secondaire sur des mutations compensatrices ? Le design avec marqueurs liés est bon, mais des contrôles additionnels pourraient être nécessaires.; Le critère de discordance phylogénétique (≥15% des nœuds) est-il biologiquement pertinent et statistiquement robuste ? Avez-vous prévu une analyse de sensibilité pour ce seuil, et comment allez-vous distinguer une discordance due à l'introgression d'une discordance due à une évolution convergente ou à une saturation mutationnelle sur le gène de résistance ?
- **Recommendation**: Je recommande l'acceptation du protocole avec des révisions mineures. Le design est globalement solide, novateur et méthodologiquement rigoureux. Avant de commencer la Phase 2, il est impératif de : 1) réaliser une analyse de puissance statistique formelle et d'envisager d'augmenter le nombre de réplicats à 150-200 pour la condition goulot d'étranglement ; 2) détailler et valider le protocole de mesure du coefficient de sélection (s) avec réplicats techniques et biologiques. La Phase 3 devrait inclure un plan de contrôle supplémentaire pour vérifier le mécanisme de recombinaison homologue (par exemple, en invalidant un gène clé de la recombinaison dans une lignée témoin).

#### Domain Expert

- **Score**: 7.0/10 | **Verdict**: weak_accept | **Confidence**: 0.8
- **Strengths**: L'hypothèse intègre avec succès trois cadres théoriques bien établis mais souvent distincts — la théorie quasi neutre (Ohta), la dynamique de l'effet fondateur et des goulots d'étranglement (Mayr), et la phylogénétique de réseau — en un modèle mécaniste cohérent pour un problème appliqué pressant (la résistance aux antifongiques). Cette synthèse est intellectuellement convaincante.; Le mécanisme proposé est biologiquement plausible et s'accorde avec l'écologie fongique connue. La reconnaissance des niches cliniques et écologiques (biofilms, tissus spécifiques) comme sites potentiels de goulots d'étranglement sévères est judicieuse et témoigne d'une compréhension de la biologie des populations fongiques allant au-delà des simples modèles in vitro.; Le cadre quantitatif (Nₑ = 50-200, |s| < 0,01, taux de recombinaison ≥ 10⁻⁸) fournit des prédictions testables et fait sortir l'hypothèse d'un modèle purement verbal. La prédiction d'une discordance phylogénétique ≥ 15 % constitue un résultat concret et falsifiable pour les études génomiques.
- **Weaknesses**: La base de preuves est davantage analogique que directe. L'article fondateur (Joly, 2010) traite d'un modèle général de spéciation, et non de l'évolution de la résistance chez les champignons. Les autres citations fournissent un contexte étayant (diversité, recombinaison, épidémiologie) mais ne constituent pas des preuves de la chaîne causale proposée. L'hypothèse manque actuellement de soutien empirique direct issu de systèmes fongiques.; Les paramètres clés sont présentés comme des postulats mais constituent des inconnues critiques. La Distribution of Fitness Effects (DFE) des mutations de résistance aux antifongiques in vivo est mal caractérisée ; certaines mutations de résistance peuvent avoir |s| >> 0,01, rendant improbable une fixation médiée par la dérive. De même, le taux de recombinaison entre lignées divergentes en conditions naturelles est souvent bien inférieur à 10⁻⁸, ce qui pourrait rompre l'étape d'introgression.; L'hypothèse peut sous-estimer le rôle de la sélection directe. En contexte clinique, les goulots d'étranglement (par exemple, la transmission) coïncident souvent avec l'exposition aux médicaments (prophylaxie ou traitement), créant un scénario où la sélection pour la résistance est forte et simultanée. Distinguer les effets de la dérive d'une sélection de fond faible dans ces environnements complexes constitue un défi méthodologique majeur.
- **Questions**: Étant donné que de nombreuses mutations de résistance aux antifongiques (par exemple, dans ERG11) sont connues pour conférer des coûts de fitness significativement supérieurs à 1 % in vitro, comment le postulat central de |s| < 0,01 peut-il être justifié pour les allèles pertinents dans les environnements de niche proposés ? Pouvez-vous fournir des estimations issues de la littérature ou des données préliminaires ?; L'hypothèse postule que l'introgression par recombinaison homologue constitue une étape finale clé. Cependant, de nombreux champignons pathogènes comme Candida auris présentent des structures de population majoritairement clonales, avec une recombinaison limitée entre les principaux clades. Comment le mécanisme se déroulerait-il chez de telles espèces, ou l'hypothèse est-elle spécifique aux pathogènes qui se recombinent ?; Comment distinguer opérationnellement la signature de ce modèle « goulot d'étranglement–fixation–introgression » de modèles alternatifs, tels que (a) la sélection sur la variation préexistante dans une grande population, ou (b) la mutation de novo et la sélection au cours de l'événement de goulot d'étranglement lui-même, à l'aide de jeux de données génomiques standard (par exemple, une collection d'isolats cliniques) ?
- **Recommendation**: L'hypothèse présente une synthèse novatrice et théoriquement solide qui mérite d'être étudiée. Je recommande de poursuivre avec un programme de recherche ciblé pour tester ses postulats centraux. La priorité devrait être accordée à : 1) Quantifier les effets sur la fitness des allèles de résistance dans des conditions pertinentes imitant l'hôte et sans médicament, afin d'affiner la DFE. 2) Développer des simulations de génomique des populations pour établir les signaux phylogénétiques attendus (discordance de réseau) sous ce modèle par rapport aux alternatives. 3) Le séquençage ciblé de populations fongiques issues de niches putatives de goulot d'étranglement (par exemple, des biofilms sur un hôte unique au cours du temps) afin de rechercher des preuves directes de décalages de fréquence allélique induits par la dérive.

#### Contrarian

- **Score**: 4.5/10 | **Verdict**: weak_reject | **Confidence**: 0.85
- **Strengths**: L'hypothèse établit élégamment un pont entre la théorie de la génétique des populations (dérive dans les goulots d'étranglement) et un problème clinique pressant (résistance aux antifongiques), en proposant une voie non adaptative souvent négligée.; Les prédictions proposées sont formulées de manière quantitative et testables en principe, dépassant les observations purement corrélationnelles pour proposer des expériences spécifiques et des signatures bio-informatiques.
- **Weaknesses**: FAIL REASON #1: Le taux de recombinaison proposé (≥ 10⁻⁸ événements/bp/génération) est probablement insuffisant pour une introgression efficace dans le laps de temps proposé. De nombreux champignons pathogènes, dont Candida auris, présentent des structures de population majoritairement clonales, avec une recombinaison rare et restreinte. La probabilité qu'un allèle rare, ayant dérivé dans une minuscule sous-population, se recombine avec succès dans un fond génétique progéniteur génétiquement distinct avant d'être à nouveau perdu est extrêmement faible, ce qui fait de la chaîne causale allant de l'étape 3 à l'étape 4 le point de rupture critique.; FAIL REASON #2: La distribution supposée du coût de fitness (|s| < 0,01) est optimiste et probablement irréaliste pour de nombreux allèles de résistance dans un environnement *sans médicament*, qui est la condition pertinente pour que la dérive l'emporte sur la sélection. Les mutations connues d'ERG11 portent souvent un coût de fitness mesurable in vitro. Si le coût réel est s = -0,05, et non -0,005, la dérive dans un Nₑ de 200 devient négligeable, et l'allèle est éliminé même dans le goulot d'étranglement.; FAIL REASON #3: La discordance phylogénétique prédite (≥ 15 %) est un motif post-hoc attendu de nombreux processus (par exemple, évolution convergente, transfert horizontal de gènes, tri de lignées incomplet). Sa détection ne valide pas le mécanisme de goulot d'étranglement-introgression proposé. Il s'agit d'une prédiction soumise à des facteurs de confusion ; l'hypothèse manque d'une signature unique et falsifiable permettant de distinguer son mécanisme spécifique de ces alternatives.
- **Questions**: Comment prévoyez-vous de distinguer empiriquement la signature de votre goulot d'étranglement proposé comme moteur de l'introgression de l'explication bien plus parcimonieuse d'une mutation *de novo* suivie d'une expansion clonale sous sélection médicamenteuse, en particulier compte tenu des faibles taux de recombinaison ?; Votre mécanisme repose sur l'existence d'allèles de résistance « préexistants » à des fréquences de 0,001-0,05. Quelles preuves génomiques directes existent pour une telle variation présente dans les populations cliniques pour les cibles antifongiques clés, étant donné que même un séquençage profond d'échantillons de patients peut ne pas capturer les variants présents dans une future sous-population soumise à un goulot d'étranglement ?
- **Recommendation**: Pour être crédible, les auteurs doivent d'abord démontrer la faisabilité de l'étape d'introgression. Cela nécessite un système d'évolution expérimentale dans lequel deux lignées marquées (l'une portant un allèle de résistance ayant dérivé) sont co-cultivées sans sélection, avec des taux de recombinaison et des fréquences d'introgression mesurés directement par séquençage. Sans validation empirique de cette étape clé, limitante, l'hypothèse demeure un scénario théorique intéressant mais non étayé.

#### Industrialist

- **Score**: 4.5/10 | **Verdict**: weak_reject | **Confidence**: 0.75
- **Strengths**: Adresse un besoin critique non satisfait dans la lutte antifongique : la prediction de l'emergence de resistances 'stochastiques' echappant aux modeles deterministes, visant directement les departements R&D des grands groupes pharmaceutiques (Pfizer, Gilead, MSD) et les agences de sante publique (CDC, ECDC).; Avantage competitif potentiel dans la surveillance epidemiologique : une plateforme predictive basee sur ce mecanisme pourrait identifier les souches a 'haut risque d'introgression' avant leur emergence clinique, offrant un delai d'alerte de 6-18 mois sur les methodes actuelles de sequencage reactif.
- **Weaknesses**: Marche adressable trop niche et difficile a monetiser directement. Le client final (hopital) ne paie pas pour des modeles evolutifs, et le transfert vers un kit diagnostique ou un logiciel SaaS rencontre une competition farouche d'acteurs etablis (QIAGEN, bioMerieux) avec des barrieres reglementaires (FDA/CE-IVD) insurmontables a court terme.; Le TRL est extremement bas (2-3). La voie de commercialisation la plus realiste est la vente de service de R&D ou de licence de propriete intellectuelle a un grand pharma, un modele a faible marge et a cycle de vente long (>3 ans), avec un ROI incertain face aux priorites internes des acheteurs potentiels.
- **Questions**: Quel est le business model concret ? Vendez-vous un algorithme (SaaS aux labos academiques), un service d'analyse contractuelle (CRO), ou visez-vous un partenariat de co-developpement avec un acteur du diagnostic ? Chaque voie exige des competences et un capital radicalement differents.; Quelle est la strategie de protection IP ? Le mecanisme evolutif lui-meme n'est pas brevetable. Seul un protocole experimental specifique ou un algorithme predictif unique pourrait l'etre, offrant une barriere a l'entree faible face a des concurrents mieux capitalises.
- **Recommendation**: Reorienter l'hypothese vers un service de recherche sous contrat (CRO) pour les departements de resistance antifongique des grands pharmas, en utilisant le protocole valide comme preuve de concept. Parallellement, developper un algorithme predictif minimal (MVP) et le proposer en acces libre pour etablir un standard academique et capturer l'expertise, en visant a long terme une licence exclusive pour un developpeur de logiciel de surveillance epidemiologique (ex : BioNumerics).

#### Funding Strategist

- **Score**: 6.5/10 | **Verdict**: weak_accept | **Confidence**: 0.8
- **Strengths**: Hypothèse mécanistique originale qui relie la génétique des populations, l'écologie évolutive et la santé publique, un angle très apprécié dans les appels interdisciplinaires.; Protocole bien structuré avec des critères GO/NO-GO clairs, démontrant une approche par étapes et une gestion des risques adaptée à un projet exploratoire.
- **Weaknesses**: TRL actuel très bas (TRL 1-2). Le budget et la timeline proposés correspondent à une preuve de concept préliminaire, insuffisante pour un appel collaboratif majeur type Horizon Europe.; Le narratif actuel est trop focalisé sur la validation technique. Pour un appel compétitif, il faut renforcer l'impact sociétal (résistance aux antifongiques = menace sanitaire mondiale) et le lien avec les données cliniques/épidémiologiques existantes.
- **Questions**: Quelle est la preuve préliminaire, au-delà de la simulation, que ces événements de goulot d'étranglement se produisent dans des niches cliniques ou environnementales réelles (biofilms, infection chronique, traitement antifongique) ?; Le consortium idéal manque d'un partenaire clinique ou de santé publique pour valider la pertinence des souches et le modèle écologique. Qui pourrait remplir ce rôle pour renforcer la proposition ?
- **Recommendation**: Cibler d'abord un financement de démarrage pour atteindre un TRL 3-4 avec des données expérimentales solides sur un modèle fongique. En parallèle, constituer un consortium exploratoire avec un modélisateur, un généticien des populations fongiques et un clinicien/microbiologiste hospitalier. Présenter ensuite le projet comme une recherche fondamentale avec un potentiel de transformation pour la surveillance et la gestion de la résistance.

### B. Semantic Scholar Search Queries

- [novelty] `budding speciation fungal pathogens` — Direct search for the core hypothesis linking the evolutionary model to mycology.
- [novelty] `founder effect antifungal resistance evolution` — Tests if the specific mechanism of stochastic origin via bottlenecks is documented.
- [novelty] `genetic drift azole resistance emergence` — Searches for evidence of non-selective, stochastic processes in resistance evolution.
- [novelty] `introgression antifungal resistance spread` — Tests the 'post-divergence gene flow' component of the hypothesis for spreading alleles.
- [evidence] `Candida biofilm population bottleneck genomics` — Seeks evidence for small, isolated subpopulations (biofilms) with founder effects in fungi.
- [evidence] `standing genetic variation antifungal resistance` — Looks for pre-existing polymorphisms that could be randomly fixed via drift in a bud.
- [evidence] `fungal hybridization drug resistance transfer` — Seeks evidence for gene flow/introgression as a mechanism for spreading resistance alleles.
- [evidence] `clonal expansion vs selection resistance` — Investigates alternative models to selective sweeps in resistance emergence.
- [cross_domain] `Rhodiola budding speciation genomics` — Finds the foundational plant studies on the model to understand its genetic mechanisms.
- [cross_domain] `peripatric speciation microbial pathogens` — Searches for application of similar speciation models (small isolate origin) in other microbes like bacteria.
- [cross_domain] `founder effect antibiotic resistance evolution` — Looks for a direct precedent of the proposed mechanism in bacterial systems, a closely related field.

---

*Generated by SPORE (Systeme de Production d'Opportunites de Recherche par Exploration) on 2026-04-19*