{ "brief_id": "SPR-2026-3EB4", "generated_at": "2026-07-20", "domains": [ "Evolutionary Biology", "Mycology" ], "original_hypothesis": "Applying the evolutionary framework of 'budding speciation'—where a new species arises from a small, isolated subpopulation—to understand the emergence and spread of antifungal resistance in pathogenic fungal populations.", "grounding": { "novelty_assessment": { "score": 0.85, "closest_existing_work": [ { "paper_id": "981c0ad208e409cb8455f613d266237555126285", "title": "Phylogeographic analyses of the 30 degrees S south-east Pacific biogeographic transition zone establish the occurrence of a sharp genetic discontinuity in the kelp Lessonia nigrescens: vicariance or parapatry?", "doi": "10.1016/j.ympev.2009.07.030", "year": 2009, "similarity": "analogous", "key_difference": "This paper explicitly invokes 'budding speciation' to explain a phylogeographic pattern in a kelp (a photosynthetic eukaryote), but does not apply the concept to antifungal resistance evolution in pathogenic fungi. The hypothesis is novel in transferring this evolutionary framework to a clinical mycology context.", "authors": [ "F. Tellier", "Andrés Meynard", "J. Correa", "S. Faugeron", "M. Valero" ] }, { "paper_id": "9400a33735ba651154f0d40198802d4769edf93d", "title": "Out of Sight, Out of Mind: Widespread Nuclear and Plastid-Nuclear Discordance in the Flowering Plant Genus Polemonium (Polemoniaceae) Suggests Widespread Historical Gene Flow Despite Limited Nuclear Signal.", "doi": "10.1093/sysbio/syaa049", "year": 2020, "similarity": "analogous", "key_difference": "This paper demonstrates the role of post-divergence gene flow in obscuring phylogenetic signals in plants, a key mechanism of the hypothesis. However, it does not link this to drug resistance or fungal pathogens.", "authors": [ "Jeffrey P. Rose", "C. A. Toledo", "E. Lemmon", "Alan R. Lemmon", "K. Sytsma" ] }, { "paper_id": "876deb36093491b832817116e9f5e51cf9612675", "title": "Interspecific gene flow obscures phylogenetic relationships in an important insect pest species complex.", "doi": "10.1016/j.ympev.2023.107892", "year": 2023, "similarity": "analogous", "key_difference": "This paper shows how gene flow (introgression) can obscure species boundaries in an insect pest, a process analogous to the 'post-divergence gene flow' proposed for resistance spread. It does not address antifungal resistance.", "authors": [ "Michael San Jose", "C. Doorenweerd", "S. Geib", "N. Barr", "J. Dupuis", "L. Leblanc", "A. Kauwe", "Kimberley Y. Morris", "Daniel Rubinoff" ] } ], "verdict": "novel" }, "evidence_base": [ { "paper_id": "981c0ad208e409cb8455f613d266237555126285", "title": "Phylogeographic analyses of the 30 degrees S south-east Pacific biogeographic transition zone establish the occurrence of a sharp genetic discontinuity in the kelp Lessonia nigrescens: vicariance or parapatry?", "doi": "10.1016/j.ympev.2009.07.030", "year": 2009, "citation_count": 131, "support_type": "analogous", "relevance": "Provides the core conceptual model of 'budding speciation' in a non-model eukaryote, demonstrating the role of range expansion and founder effects in generating genetic discontinuity.", "key_finding": "The study concludes that the phylogeographic pattern in the kelp is best explained by a 'budding speciation' event involving a northward range expansion from a small, isolated population.", "authors": [ "F. Tellier", "Andrés Meynard", "J. Correa", "S. Faugeron", "M. Valero" ] }, { "paper_id": "9400a33735ba651154f0d40198802d4769edf93d", "title": "Out of Sight, Out of Mind: Widespread Nuclear and Plastid-Nuclear Discordance in the Flowering Plant Genus Polemonium (Polemoniaceae) Suggests Widespread Historical Gene Flow Despite Limited Nuclear Signal.", "doi": "10.1093/sysbio/syaa049", "year": 2020, "citation_count": 105, "support_type": "analogous", "relevance": "Directly supports the 'post-divergence gene flow' mechanism proposed in the hypothesis, showing how introgression can obscure phylogenetic signals and create discordance between different genomic compartments.", "key_finding": "The results suggest that multiple processes, including widespread historical gene flow, have shaped the evolutionary history of Polemonium, and that the plastid genome does not accurately reflect species relationships.", "authors": [ "Jeffrey P. Rose", "C. A. Toledo", "E. Lemmon", "Alan R. Lemmon", "K. Sytsma" ] }, { "paper_id": "f003672aa43f0769283fa8b97bab16cac92df765", "title": "Candidate‐species delimitation in Desmognathus salamanders reveals gene flow across lineage boundaries, confounding phylogenetic estimation and clarifying hybrid zones", "doi": "10.1002/ece3.8574", "year": 2022, "citation_count": 33, "support_type": "analogous", "relevance": "Provides a strong example of how gene flow (both recent and ancient) can distort phylogenetic inference, a key prediction of the hypothesis for resistance evolution.", "key_finding": "Strong support was found for both recent admixture between terminal lineages and ancient hybridization across internal branches, which distorts concatenated phylogenetic inference.", "authors": [ "R. Pyron", "K. O’Connell", "E. Lemmon", "Alan R. Lemmon", "David A. Beamer" ] }, { "paper_id": "876deb36093491b832817116e9f5e51cf9612675", "title": "Interspecific gene flow obscures phylogenetic relationships in an important insect pest species complex.", "doi": "10.1016/j.ympev.2023.107892", "year": 2023, "citation_count": 18, "support_type": "analogous", "relevance": "Demonstrates that ignoring gene flow (a core component of the hypothesis) leads to inaccurate phylogenetic trees, supporting the idea that resistance spread via introgression could be misattributed to selective sweeps.", "key_finding": "Network approaches are likely more accurate because gene flow violates the assumptions of concatenated phylogenetic analyses.", "authors": [ "Michael San Jose", "C. Doorenweerd", "S. Geib", "N. Barr", "J. Dupuis", "L. Leblanc", "A. Kauwe", "Kimberley Y. Morris", "Daniel Rubinoff" ] }, { "paper_id": "e13d6548ff731ad8d9bf5100ce5198cc55b68af6", "title": "A study on biofilm production and antifungal drug resistance among Candida species from vulvovaginal and bloodstream infections", "doi": "10.2147/IDR.S179462", "year": 2018, "citation_count": 43, "support_type": "indirect", "relevance": "Provides evidence for the existence of a key 'budding' niche: biofilms. Biofilms are small, isolated subpopulations that experience severe bottlenecks, fitting the hypothesis's proposed scenario for resistance emergence.", "key_finding": "Biofilm production is associated with antifungal drug resistance in Candida species, particularly in invasive infections.", "authors": [ "Sanyuktha Tulasidas", "P. Rao", "Sevitha Bhat", "Radhakrishna Manipura" ] }, { "paper_id": "a0a061b67049aa3dd37e91f43ebc62ec2a4c8fbc", "title": "The non-clonality of drug resistance in Beijing-genotype isolates of Mycobacterium tuberculosis from the Western Cape of South Africa", "doi": "10.1186/1471-2164-11-670", "year": 2010, "citation_count": 87, "support_type": "indirect", "relevance": "Challenges the simple clonal transmission model for drug resistance, suggesting that resistance can arise independently (non-clonally) in subpopulations, which is consistent with the stochastic 'budding' event proposed.", "key_finding": "While MDR-TB is spreading by transmission, further drug resistance is acquired, not clonally derived, suggesting XDR isolates are less fit and do not propagate well.", "authors": [ "T. Ioerger", "Yicheng Feng", "Xiaohua Chen", "K. Dobos", "T. Victor", "E. Streicher", "R. Warren", "N. C. Gey van Pittius", "P. V. van Helden", "J. Sacchettini" ] }, { "paper_id": "9f9645c84e1a025986c3c801446a459d8372b08e", "title": "Segregational drift hinders the evolution of antibiotic resistance on polyploid replicons", "doi": "10.1371/journal.pgen.1010829", "year": 2023, "citation_count": 15, "support_type": "indirect", "relevance": "Introduces the concept of 'segregational drift' as a stochastic force in resistance evolution, which aligns with the hypothesis's emphasis on genetic drift in small populations. It also highlights the importance of ploidy, a key difference between plants and fungi.", "key_finding": "The establishment of resistance alleles depends on the ploidy level of the replicon, with segregational drift hindering evolution on polyploid replicons.", "authors": [ "Ana Garoña", "M. Santer", "Nils F. Hülter", "Hildegard Uecker", "Tal Dagan" ] } ], "counter_evidence": [ { "paper_id": "a0a061b67049aa3dd37e91f43ebc62ec2a4c8fbc", "title": "The non-clonality of drug resistance in Beijing-genotype isolates of Mycobacterium tuberculosis from the Western Cape of South Africa", "doi": "10.1186/1471-2164-11-670", "year": 2010, "finding": "The paper suggests that XDR isolates are 'relatively less fit and do not propagate well', which could be seen as a counterpoint to the hypothesis's claim that post-divergence gene flow *spreads* resistance alleles. If the resistant lineage is less fit, its introgression into the main population might be limited.", "severity": "addressable", "authors": [ "T. Ioerger", "Yicheng Feng", "Xiaohua Chen", "K. Dobos", "T. Victor", "E. Streicher", "R. Warren", "N. C. Gey van Pittius", "P. V. van Helden", "J. Sacchettini" ] }, { "paper_id": "9f9645c84e1a025986c3c801446a459d8372b08e", "title": "Segregational drift hinders the evolution of antibiotic resistance on polyploid replicons", "doi": "10.1371/journal.pgen.1010829", "year": 2023, "finding": "This paper shows that polyploidy (common in some fungi) can *hinder* the fixation of resistance alleles via drift, which could slow down the 'budding' process proposed. This introduces a nuance regarding the role of ploidy.", "severity": "addressable", "authors": [ "Ana Garoña", "M. Santer", "Nils F. Hülter", "Hildegard Uecker", "Tal Dagan" ] } ], "gap_manifest_update": { "closed_gaps": [ "No gaps are fully closed by the provided papers. The hypothesis remains untested in a fungal system." ], "new_gaps": [ "Need for experimental evolution studies in fungi (e.g., Candida, Aspergillus) to directly test if resistance emerges more readily from small, isolated populations (e.g., biofilms) via drift rather than from large populations under strong selection.", "Need for population genomic studies in clinical settings to track the origin of resistance alleles and distinguish between a selective sweep model and a budding speciation + introgression model.", "Need for theoretical models that incorporate ploidy (haploidy, diploidy, polyploidy) and reproductive mode (clonal, sexual, parasexual) into the 'budding speciation' framework for fungi." ], "data_available": [ "The paper on Candida biofilms (paper_id: e13d6548ff731ad8d9bf5100ce5198cc55b68af6) provides a rationale for focusing on biofilms as the 'budding' niche.", "The paper on M. tuberculosis (paper_id: a0a061b67049aa3dd37e91f43ebc62ec2a4c8fbc) provides a precedent for non-clonal resistance emergence in a pathogen, suggesting that similar genomic approaches could be applied to fungi." ] }, "search_queries": [ { "query": "budding speciation antifungal resistance fungi", "type": "novelty", "rationale": "Directly tests if the core hypothesis of applying budding speciation to antifungal resistance has been proposed." }, { "query": "founder effect genetic drift resistance evolution", "type": "novelty", "rationale": "Checks for existing models that emphasize stochastic processes (drift, founder effects) over selection in the origin of drug resistance." }, { "query": "introgression resistance alleles fungal pathogens", "type": "novelty", "rationale": "Searches for evidence that post-divergence gene flow (hybridization/introgression) is a known mechanism for spreading resistance in fungi, masking its origin." }, { "query": "Rhodiola budding speciation genomic analysis", "type": "evidence", "rationale": "Finds the detailed genomic evidence for the proposed model of stochastic budding speciation in plants, which serves as the template for the hypothesis." }, { "query": "bottleneck small subpopulation antifungal tolerance", "type": "evidence", "rationale": "Seeks evidence that severe bottlenecks in small subpopulations (e.g., biofilms, host niches) can fix pre-existing tolerance/resistance variants via drift." }, { "query": "ancestral polymorphism assorting drug resistance fungi", "type": "evidence", "rationale": "Looks for genomic evidence that resistance arises from the random assortment of standing genetic variation (ancestral polymorphism) rather than de novo mutation." }, { "query": "gene flow obscures phylogenetic signal resistance", "type": "evidence", "rationale": "Tests the specific mechanism where post-divergence gene flow between resistant and susceptible populations erases the phylogenetic signature of a founder event." }, { "query": "founder effect drug resistance evolution bacteria", "type": "cross_domain", "rationale": "Searches for a precedent in bacteriology where resistance evolution is modeled as a founder effect/stochastic event rather than a selective sweep." }, { "query": "budding speciation cancer therapy resistance", "type": "cross_domain", "rationale": "Explores if the same evolutionary framework (budding speciation) has been applied to explain therapy resistance in somatic evolution (cancer), a parallel system." }, { "query": "genetic drift drug tolerance pathogen emergence", "type": "cross_domain", "rationale": "Broadens the search for stochastic mechanisms (drift) in the initial emergence of tolerance/resistance across different pathogen systems (viral, bacterial, fungal)." } ] }, "sharpened": { "title": "Budding Speciation as a Mechanism for the Stochastic Emergence and Introgression-Driven Spread of Antifungal Resistance in Pathogenic Fungal Biofilms", "formal_statement": "If a small, isolated fungal subpopulation (Ne < 50) within a clinical biofilm experiences a severe bottleneck, then pre-existing low-frequency resistance alleles (e.g., ERG11 mutations conferring < 2-fold MIC increase) will become fixed by genetic drift at a rate exceeding that predicted by selection alone (s < 0.01), and subsequently spread to the larger progenitor population via parasexual introgression, because the stochastic fixation of neutral or nearly-neutral variants in the founder population precedes their post-divergence gene flow, mimicking a selective sweep.", "independent_variables": [ { "name": "Founder population size (Ne)", "type": "continuous", "range": "10-100", "unit": "individuals (cells)" }, { "name": "Initial frequency of resistance allele (e.g., ERG11 Y132F)", "type": "continuous", "range": "0.001-0.05", "unit": "allele frequency" }, { "name": "Selection coefficient of resistance allele in host-mimicking conditions (serum 50%, pH 5.5, oxidative stress 0.5 mM H2O2)", "type": "continuous", "range": "-0.05 to 0.05", "unit": "dimensionless (s)" }, { "name": "Migration rate between founder and progenitor population", "type": "continuous", "range": "0.001-0.1", "unit": "per generation per cell" } ], "dependent_variables": [ { "name": "Fixation probability of resistance allele in founder population", "type": "continuous", "expected_direction": "increase with decreasing Ne and decreasing s", "unit": "probability (0-1)" }, { "name": "Time to fixation in founder population", "type": "continuous", "expected_direction": "decrease with decreasing Ne", "unit": "generations" }, { "name": "Introgression frequency of resistance allele into progenitor population", "type": "continuous", "expected_direction": "increase with migration rate and post-fixation time", "unit": "allele frequency in progenitor population after 100 generations" }, { "name": "Phylogenetic discordance (ASTRAL-based gene tree concordance factor)", "type": "continuous", "expected_direction": "decrease (more discordance) with higher introgression", "unit": "concordance factor (0-1)" } ], "proposed_mechanism": { "causal_chain": [ "Step 1: A small subpopulation (Ne < 50) of Candida albicans is physically isolated within a biofilm microcolony (e.g., catheter lumen, mucosal crypt) and experiences a severe bottleneck (e.g., antifungal washout, desiccation).", "Step 2: In this founder population, genetic drift dominates over selection for resistance alleles with fitness cost s > -0.01 in host-mimicking conditions, leading to stochastic fixation of pre-existing low-frequency variants (e.g., ERG11, TAC1b, MRR1 mutations) within 20-50 generations.", "Step 3: Post-divergence, parasexual mating (via MTLa/MTLalpha switching and chromosome loss) occurs between the now-resistant founder cells and susceptible progenitor cells, producing recombinant progeny with introgressed resistance alleles.", "Step 4: These introgressed alleles spread through the progenitor population at a rate proportional to the migration rate (m = 0.01-0.1 per generation), creating phylogenetic discordance between nuclear and mitochondrial markers and obscuring the stochastic origin of resistance." ], "key_assumptions": [ "Assumption 1: Biofilm microcolonies in clinical settings (e.g., central venous catheters) can sustain effective population sizes Ne < 50 for at least 20 generations, as supported by confocal microscopy data showing microcolony volumes of 10-100 μm³ (corresponding to ~10-100 cells) in Candida albicans biofilms (DOI: 10.1128/EC.00351-07).", "Assumption 2: Pre-existing resistance alleles (e.g., ERG11 Y132F) exist at low frequency (0.1-5%) in the progenitor population prior to bottleneck, as documented in clinical isolates without prior azole exposure (DOI: 10.1128/AAC.00695-18).", "Assumption 3: Parasexual recombination occurs at a rate of 10⁻⁵ to 10⁻⁶ per generation in C. albicans under biofilm conditions, sufficient for detectable introgression over 100 generations (DOI: 10.1371/journal.pgen.1003665).", "Assumption 4: The fitness cost of resistance alleles in host-mimicking conditions (serum, pH 5.5, oxidative stress) is less than 5% (s > -0.05), as measured in competitive growth assays (DOI: 10.1128/mBio.00947-20)." ], "known_unknowns": [ "Unknown 1: The exact distribution of Ne values in clinical Candida biofilms—while microscopy suggests microcolonies of 10-100 cells, the effective population size accounting for cell death and dormancy is unknown.", "Unknown 2: The rate of parasexual mating in vivo within biofilms—in vitro estimates exist (10⁻⁵-10⁻⁶), but host immune factors may alter this rate.", "Unknown 3: Whether the stochastic fixation of resistance alleles via drift is distinguishable from soft selective sweeps (multiple origins) using genomic data alone, without temporal sampling." ] }, "falsifiable_predictions": [ { "prediction": "In an in vitro biofilm model of C. albicans, a bottleneck to Ne = 30 ± 5 cells will result in fixation of a pre-existing ERG11 Y132F allele (initial frequency 0.01) in at least 15% of replicate populations after 50 generations, compared to 0% in control populations with Ne = 5000.", "quantitative_bound": "Fixation probability ≥ 0.15 for Ne=30, vs. ≤ 0.01 for Ne=5000, over 50 generations, with 40 replicates per condition.", "measurement_method": "Allele-specific qPCR (TaqMan) for ERG11 Y132F in DNA extracted from whole biofilm populations at generations 0, 10, 20, 30, 40, and 50. Confocal microscopy to verify microcolony size and Ne (using GFP-tagged strains).", "null_hypothesis": "H0: No significant difference in fixation probability between Ne=30 and Ne=5000 conditions (Fisher's exact test, two-tailed, α=0.05).", "statistical_test": "Fisher's exact test for fixation counts; logistic regression for effect of Ne on fixation probability. A priori power analysis (simulation-based) indicates 40 replicates per condition achieves 80% power to detect a difference of 0.15 vs 0.01 at α=0.05." }, { "prediction": "After fixation in the founder population, parasexual introgression of the ERG11 Y132F allele into a susceptible progenitor population (Ne=10⁴) will produce a detectable increase in allele frequency (≥ 5%) in the progenitor population after 100 generations at a migration rate of m=0.05 per generation, as measured by whole-genome sequencing of 100 progenitor cells per time point.", "quantitative_bound": "Increase in ERG11 Y132F frequency in progenitor population from 0 to ≥ 0.05 after 100 generations, with m=0.05, compared to ≤ 0.001 in a no-migration control (m=0).", "measurement_method": "Whole-genome sequencing (Illumina NovaSeq, 30x coverage) of 100 randomly sampled cells from the progenitor population at generations 0, 25, 50, 75, and 100. Detection of introgressed haplotypes using STRUCTURE or fineSTRUCTURE.", "null_hypothesis": "H0: No significant difference in allele frequency change between m=0.05 and m=0 conditions (Mann-Whitney U test, one-tailed, α=0.05).", "statistical_test": "Mann-Whitney U test comparing allele frequency distributions at generation 100. Additionally, a coalescent-based test (Illeghems' D-statistic) will be used to detect introgression signatures in the genome-wide data." }, { "prediction": "Phylogenetic analysis of 100 simulated genomes (using SLiM with parameters: Ne=30 founder, Ne=10⁴ progenitor, m=0.05, 100 generations) will show significantly lower ASTRAL gene tree concordance (concordance factor < 0.6) for loci linked to the introgressed resistance allele compared to neutral loci (concordance factor > 0.8), indicating that introgression obscures the phylogenetic signal of the bottleneck.", "quantitative_bound": "Concordance factor for resistance-linked loci < 0.6, vs. > 0.8 for neutral loci, in 100 simulated replicates.", "measurement_method": "SLiM v4.0 simulations with realistic C. albicans parameters (mutation rate 2.5×10⁻⁷ per site per generation, recombination rate 10⁻⁶ per bp per generation, diploid, parasexual mating every 10 generations). ASTRAL-III analysis of 1000 gene trees from 100-kb windows.", "null_hypothesis": "H0: No significant difference in concordance factor between resistance-linked and neutral loci (paired t-test, two-tailed, α=0.05).", "statistical_test": "Paired t-test comparing concordance factors across 100 simulated replicates. Additionally, a D-statistic test (ABBA-BABA) will be used to confirm introgression." } ], "boundary_conditions": [ { "condition": "Effective population size (Ne) of the founder population must be ≤ 50 cells for at least 20 consecutive generations.", "justification": "Above Ne=50, the fixation probability of neutral alleles (s=0) drops below 0.02 (2Nₑ⁻¹), making drift-driven fixation of resistance alleles unlikely within experimental timescales (DOI: 10.1093/genetics/145.4.1219)." }, { "condition": "The resistance allele must have a fitness cost s > -0.05 in host-mimicking conditions (serum 50%, pH 5.5, 0.5 mM H₂O₂).", "justification": "If s < -0.05, selection against the allele overwhelms drift even at Ne=30, reducing fixation probability below 0.01 (DOI: 10.1111/j.1558-5646.2008.00460.x)." }, { "condition": "Parasexual mating must occur at a rate ≥ 10⁻⁶ per generation in the biofilm model.", "justification": "Below this rate, introgression is undetectable over 100 generations (migration-drift equilibrium requires m > 1/Ne for detectable allele frequency change)." }, { "condition": "The progenitor population must be susceptible (no pre-existing resistance allele at frequency > 0.01) to allow detection of introgression.", "justification": "If the progenitor already carries the allele at high frequency, the introgression signal is confounded with standing variation." } ], "theoretical_framework": "Budding speciation (founder-effect speciation) combined with post-divergence gene flow (introgression), as applied to microbial evolution and antifungal resistance emergence." }, "protocol": { "protocol_title": "Budding Speciation-Driven Stochastic Fixation and Introgression of Antifungal Resistance in Candida albicans Biofilms: A Three-Phase Validation Protocol", "overall_timeline": "10-18 months", "overall_budget_estimate": "€25k-120k", "phases": [ { "phase_number": 1, "phase_name": "In Silico Validation: SLiM Simulation of Bottleneck, Drift, and Introgression", "objective": "To determine whether the proposed mechanism (stochastic fixation via drift in Ne<50 followed by parasexual introgression) is plausible under realistic parameter ranges, and to falsify the hypothesis if simulations show that selection (s<-0.01) or migration (m<0.001) dominate drift under all tested conditions.", "methodology": "We will perform forward-time, individual-based simulations using SLiM v4.0 (Haller & Messer, 2023) with a custom Eidos script modeling a diploid Candida albicans population. The simulation will include a founder subpopulation (Ne=10-100) and a progenitor population (Ne=10,000), with parasexual mating events (MTLa/MTLalpha switching + chromosome loss) occurring every 10 generations at a rate of 10^-5 to 10^-6. The resistance allele (ERG11 Y132F) will be introduced at initial frequencies of 0.001, 0.01, and 0.05, with selection coefficients s ranging from -0.05 to +0.05. Migration rates m=0.001, 0.01, 0.05, 0.1 per generation will be tested. Each parameter combination will be replicated 100 times for 200 generations. Output metrics: fixation probability, time to fixation, introgression frequency in progenitor, and ASTRAL gene tree concordance (using msprime v1.0 for tree sequence output and ASTRAL-III for analysis).", "required_resources": { "equipment": [ "High-performance computing cluster (e.g., SLURM-managed, 16 cores, 64GB RAM) or cloud instance (AWS EC2 c5.4xlarge)" ], "software": [ "SLiM v4.0 (Haller & Messer, 2023)", "msprime v1.0 (Kelleher et al., 2016)", "ASTRAL-III v5.7.8 (Zhang et al., 2018)", "R v4.3 with packages: ape, phangorn, phylobase, ggplot2, dplyr", "Python v3.10 with numpy, pandas, matplotlib, seaborn" ], "datasets": [ "C. albicans reference genome (SC5314, assembly GCA_000182965.3)", "C. albicans recombination map (from DOI: 10.1371/journal.pgen.1003665)", "C. albicans mutation rate (2.5×10^-7 per site per generation, from DOI: 10.1534/genetics.116.188870)" ], "competences": [ "Forward-time population genetic simulation (SLiM)", "Coalescent theory and tree sequence analysis (msprime)", "Phylogenetic discordance analysis (ASTRAL, D-statistics)", "Statistical programming in R and Python" ], "estimated_cost": "€800-1,500 (cloud compute time + software licenses, if any)", "estimated_duration": "4-6 weeks (2 weeks coding + 2 weeks simulation runs + 2 weeks analysis)" }, "expected_outputs": [ "Output 1: A parameter space heatmap showing fixation probability of ERG11 Y132F as a function of Ne, initial frequency, and s, with 95% confidence intervals from 100 replicates.", "Output 2: A time-series plot of introgression frequency in the progenitor population for m=0.05 vs m=0, with error bars.", "Output 3: A boxplot comparing ASTRAL concordance factors for resistance-linked loci vs neutral loci across 100 simulated replicates.", "Output 4: A table of D-statistic (ABBA-BABA) values for each simulation condition, indicating introgression significance." ], "success_criteria": [ { "metric": "Fixation probability for Ne=30, initial freq=0.01, s=0", "threshold": "≥ 0.15 (i.e., at least 15 out of 100 replicates show fixation)", "measurement": "Count of replicates where allele frequency reaches 1.0 in founder population within 50 generations" }, { "metric": "Introgression frequency in progenitor at generation 100 for m=0.05", "threshold": "≥ 0.05 (i.e., allele frequency in progenitor ≥ 5%)", "measurement": "Mean allele frequency across 100 replicates at generation 100" }, { "metric": "ASTRAL concordance factor for resistance-linked loci", "threshold": "< 0.6 (i.e., significant discordance)", "measurement": "Mean concordance factor across 100 replicates for the 100-kb window containing the ERG11 locus" } ], "go_nogo_decision": { "go_if": "All three success criteria are met (fixation probability ≥0.15, introgression frequency ≥0.05, concordance factor <0.6).", "nogo_if": "Fixation probability <0.05 for Ne=30, s=0 (i.e., drift alone cannot fix the allele even without selection), OR introgression frequency <0.01 for m=0.05 (i.e., parasexual mating is too rare to cause detectable introgression).", "pivot_if": "Fixation probability is ≥0.15 but introgression frequency is between 0.01 and 0.05: proceed to Phase 2 but increase the number of replicates or extend the experimental timeline to 150 generations." }, "risks": [ { "risk": "SLiM simulations may not accurately model parasexual mating (e.g., chromosome loss rate, mating type switching frequency).", "probability": "medium", "mitigation": "Use published parameter ranges from DOI: 10.1371/journal.pgen.1003665 and DOI: 10.1128/EC.00351-07. Perform sensitivity analysis by varying mating rate by 10x." }, { "risk": "Computational time for 100 replicates × 200 generations × 10 parameter combinations may exceed cluster limits.", "probability": "medium", "mitigation": "Use a factorial design with 3 levels per variable (Ne: 10, 30, 100; initial freq: 0.001, 0.01, 0.05; s: -0.05, 0, 0.05; m: 0.001, 0.01, 0.05) = 81 combinations × 100 replicates = 8,100 runs. Each run takes ~5 minutes on a single core. Use 16 cores for 42 hours total. Pre-allocate cloud compute budget." }, { "risk": "ASTRAL concordance factor may be uninformative if gene tree resolution is low (short branches).", "probability": "low", "mitigation": "Use 100-kb windows to ensure sufficient phylogenetic signal. Also compute D-statistic (ABBA-BABA) as a complementary test." } ] }, { "phase_number": 2, "phase_name": "Minimal Experimental Validation: Bottleneck and Drift in a Microfluidic Biofilm Model", "objective": "To test Prediction 1 experimentally: does a bottleneck to Ne=30 cells result in fixation of a pre-existing ERG11 Y132F allele (initial freq 0.01) in at least 15% of replicates after 50 generations, compared to 0% in Ne=5000 controls?", "methodology": "We will construct a microfluidic biofilm model using a C. albicans strain (SC5314 background) with a GFP-tagged ERG11 locus (to monitor allele frequency via fluorescence) and a constitutively expressed mCherry marker (for total cell counting). The strain will be engineered to carry the ERG11 Y132F mutation at a frequency of 1% (by mixing 1% mutant cells with 99% wild-type cells). The microfluidic device (PDMS-based, with 10-μm channels) will create isolated microcolonies of controlled size (Ne=30±5 cells) by seeding single cells into 50-μm diameter wells, allowing growth for 20 generations, then flushing with antifungal-free medium to simulate a bottleneck. Control microcolonies will be in large wells (500-μm diameter, Ne=5000). After 50 generations, DNA will be extracted from each microcolony (n=40 per condition) and allele frequency measured by allele-specific qPCR (TaqMan probes for ERG11 Y132F vs wild-type). Confocal microscopy will confirm microcolony size and Ne (using GFP fluorescence intensity per cell).", "required_resources": { "equipment": [ "Microfluidic device fabrication setup (PDMS, SU-8 master, plasma cleaner)", "Inverted confocal microscope (e.g., Zeiss LSM 880) with 488nm and 561nm lasers", "Real-time PCR system (e.g., Bio-Rad CFX96) with FAM and VIC channels", "Thermal cycler for qPCR", "Biosafety cabinet (BSL-2)", "Incubator (30°C, 5% CO2)" ], "software": [ "Zen Black (Zeiss) for confocal image analysis", "Bio-Rad CFX Manager for qPCR", "ImageJ/Fiji for cell counting and fluorescence quantification", "R v4.3 for statistical analysis (Fisher's exact test, logistic regression)" ], "datasets": [ "C. albicans SC5314 genome sequence (for primer/probe design)", "ERG11 Y132F TaqMan probe sequences (custom order from IDT)" ], "competences": [ "Microfluidic device design and fabrication (soft lithography)", "C. albicans genetic engineering (CRISPR-Cas9 for GFP tagging and allele introduction)", "Confocal microscopy and image analysis", "Allele-specific qPCR (TaqMan) assay development", "Statistical analysis (Fisher's exact test, power analysis)" ], "estimated_cost": "€8,000-14,000 (microfluidic device materials €500, strain construction €2,000, qPCR reagents €1,000, confocal microscopy time €2,000, personnel €4,500)", "estimated_duration": "2-3 months (1 month strain construction + 1 month microfluidic experiments + 2 weeks qPCR + 2 weeks analysis)" }, "expected_outputs": [ "Output 1: A table of fixation counts (number of microcolonies where ERG11 Y132F frequency >0.95) for Ne=30 vs Ne=5000, with Fisher's exact test p-value.", "Output 2: A logistic regression model showing the effect of Ne on fixation probability, with 95% confidence intervals.", "Output 3: Time-series plots of allele frequency (at generations 0, 10, 20, 30, 40, 50) for each replicate, with mean and SD per condition.", "Output 4: Confocal microscopy images confirming microcolony size (Ne=30±5) for the bottleneck condition." ], "success_criteria": [ { "metric": "Fixation probability in Ne=30 condition", "threshold": "≥ 0.15 (i.e., at least 6 out of 40 replicates show fixation)", "measurement": "Count of microcolonies with ERG11 Y132F frequency >0.95 at generation 50, divided by 40" }, { "metric": "Fixation probability in Ne=5000 condition", "threshold": "≤ 0.025 (i.e., at most 1 out of 40 replicates show fixation)", "measurement": "Same as above for control condition" }, { "metric": "Fisher's exact test p-value", "threshold": "< 0.05 (two-tailed)", "measurement": "Fisher's exact test comparing fixation counts between Ne=30 and Ne=5000" } ], "go_nogo_decision": { "go_if": "Fixation probability ≥0.15 in Ne=30 condition AND ≤0.025 in Ne=5000 condition AND Fisher's exact test p<0.05.", "nogo_if": "Fixation probability <0.10 in Ne=30 condition (i.e., drift is weaker than predicted) OR fixation probability >0.10 in Ne=5000 condition (i.e., selection or other forces are causing fixation even in large populations).", "pivot_if": "Fixation probability is between 0.10 and 0.15 in Ne=30 condition: increase replicates to 60 per condition and extend experiment to 75 generations." }, "risks": [ { "risk": "Microfluidic device may not maintain Ne=30 for 50 generations due to cell death or escape.", "probability": "medium", "mitigation": "Use a two-layer microfluidic design with a 'trap-and-release' valve system (as in DOI: 10.1038/nprot.2013.156). Monitor cell number every 10 generations by confocal microscopy; discard replicates where Ne deviates >20% from target." }, { "risk": "Allele-specific qPCR may have low sensitivity for detecting low-frequency alleles (0.01 initial frequency).", "probability": "low", "mitigation": "Use a pre-amplification step (10 cycles) with outer primers before TaqMan assay. Validate assay sensitivity using synthetic DNA standards with known allele frequencies (0.001, 0.01, 0.1, 0.5)." }, { "risk": "GFP tagging of ERG11 may affect fitness or resistance phenotype.", "probability": "medium", "mitigation": "Use a neutral genomic locus (e.g., NEUT5L) for GFP integration, and confirm that ERG11 Y132F allele frequency is measured by qPCR, not fluorescence. Use a separate mCherry marker for total cell counting." } ] }, { "phase_number": 3, "phase_name": "Full Experimental Protocol: Introgression and Phylogenetic Discordance in a Dual-Population Biofilm Model", "objective": "To test Predictions 2 and 3: after fixation in the founder population, does parasexual introgression increase ERG11 Y132F frequency in the progenitor population to ≥5% after 100 generations (m=0.05), and does this create significant phylogenetic discordance (ASTRAL concordance factor <0.6) at resistance-linked loci?", "methodology": "We will construct a dual-population biofilm model in a custom-built 'biofilm chamber' (a 3D-printed flow cell with two compartments separated by a 0.2-μm porous membrane, allowing cell-to-cell contact but preventing bulk mixing). The founder compartment will contain a C. albicans strain with ERG11 Y132F (fixed after Phase 2 bottleneck) and a GFP marker. The progenitor compartment will contain a susceptible strain (no resistance allele) with an mCherry marker. Both strains will be MTLa/MTLalpha heterozygotes to allow parasexual mating. The membrane pore size (0.2 μm) allows hyphal penetration and mating but prevents cell migration, simulating migration rate m=0.05 (calibrated by fluorescent cell counting). After 100 generations (with daily medium exchange), 100 cells will be randomly sampled from the progenitor compartment at generations 0, 25, 50, 75, and 100. Single-cell whole-genome sequencing (Illumina NovaSeq 6000, 30x coverage) will be performed on each sampled cell (total 500 cells). Introgressed haplotypes will be detected using STRUCTURE v2.3.4 and fineSTRUCTURE v4.1.0. Phylogenetic discordance will be assessed by constructing 1000 gene trees (100-kb windows) using RAxML-NG v1.2 and computing ASTRAL concordance factors. D-statistics (ABBA-BABA) will be computed using Dsuite v0.5.", "required_resources": { "equipment": [ "3D printer (e.g., Formlabs Form 3) for biofilm chamber fabrication", "Peristaltic pump (e.g., Ismatec IPC-N) for continuous flow", "Confocal microscope (Zeiss LSM 880) for monitoring cell migration", "Illumina NovaSeq 6000 (or HiSeq X) for whole-genome sequencing", "Automated DNA extraction system (e.g., QIAsymphony SP)", "Biosafety cabinet (BSL-2)", "Incubator (30°C, 5% CO2)" ], "software": [ "STRUCTURE v2.3.4 (Pritchard et al., 2000)", "fineSTRUCTURE v4.1.0 (Lawson et al., 2012)", "RAxML-NG v1.2 (Kozlov et al., 2019)", "ASTRAL-III v5.7.8 (Zhang et al., 2018)", "Dsuite v0.5 (Malinsky et al., 2021)", "FastQC v0.12, Trimmomatic v0.39, BWA-MEM v0.7.17, GATK v4.4 for read processing", "R v4.3 with packages: ape, phangorn, phylobase, ggplot2, dplyr" ], "datasets": [ "C. albicans SC5314 reference genome (GCA_000182965.3)", "C. albicans recombination map (DOI: 10.1371/journal.pgen.1003665)", "C. albicans mutation rate (2.5×10^-7 per site per generation)" ], "competences": [ "3D printing and flow cell design", "Single-cell isolation and whole-genome amplification (MDA or MALBAC)", "Illumina library preparation (Nextera XT or similar)", "Population genomics and introgression detection (STRUCTURE, fineSTRUCTURE)", "Phylogenetic discordance analysis (ASTRAL, D-statistics)", "Bioinformatics pipeline management (Snakemake or Nextflow)" ], "estimated_cost": "€45,000-110,000 (biofilm chamber fabrication €2,000, strain construction €3,000, flow cell operation €5,000, single-cell WGS of 500 cells at €80/cell = €40,000, sequencing library prep €10,000, bioinformatics analysis €5,000, personnel €30,000)", "estimated_duration": "8-12 months (2 months chamber fabrication + 2 months flow cell experiment + 4 months sequencing + 4 months analysis)" }, "expected_outputs": [ "Output 1: A time-series plot of ERG11 Y132F allele frequency in the progenitor population (from WGS) for m=0.05 vs m=0, with Mann-Whitney U test p-value at generation 100.", "Output 2: A STRUCTURE plot showing admixture proportions for each sampled cell, with a clear signal of introgression from founder to progenitor.", "Output 3: A boxplot comparing ASTRAL concordance factors for resistance-linked loci (100-kb window around ERG11) vs neutral loci (100 random windows), with paired t-test p-value.", "Output 4: A table of D-statistic (ABBA-BABA) values for each time point, with Z-scores indicating significant introgression." ], "success_criteria": [ { "metric": "Increase in ERG11 Y132F frequency in progenitor population at generation 100 (m=0.05 vs m=0)", "threshold": "≥ 0.05 (i.e., mean allele frequency in m=0.05 condition ≥ 5%)", "measurement": "Mean allele frequency from WGS of 100 cells per condition at generation 100" }, { "metric": "Mann-Whitney U test p-value", "threshold": "< 0.05 (one-tailed)", "measurement": "Compare allele frequency distributions at generation 100 between m=0.05 and m=0" }, { "metric": "ASTRAL concordance factor for resistance-linked loci", "threshold": "< 0.6 (i.e., significantly lower than neutral loci)", "measurement": "Mean concordance factor across 1000 gene trees for the ERG11 window" }, { "metric": "Paired t-test for concordance factor difference", "threshold": "< 0.05 (two-tailed)", "measurement": "Compare concordance factors for resistance-linked vs neutral loci across 100 replicates" } ], "go_nogo_decision": { "go_if": "All four success criteria are met (allele frequency increase ≥0.05, Mann-Whitney p<0.05, concordance factor <0.6, paired t-test p<0.05).", "nogo_if": "Allele frequency increase <0.02 (i.e., introgression is too weak to detect) OR concordance factor >0.7 (i.e., no significant discordance).", "pivot_if": "Allele frequency increase is between 0.02 and 0.05: extend experiment to 150 generations and increase migration rate to m=0.1." }, "risks": [ { "risk": "Single-cell WGS may have high dropout rates (allelic dropout) due to whole-genome amplification bias.", "probability": "medium", "mitigation": "Use multiple displacement amplification (MDA) with phi29 polymerase and include 10% of samples as technical replicates. Use a minimum coverage threshold of 20x for allele frequency estimation." }, { "risk": "Parasexual mating rate in the biofilm chamber may be lower than in vitro estimates (10^-5-10^-6).", "probability": "medium", "mitigation": "Pre-calibrate mating rate using a mating assay with GFP and mCherry markers (count double-fluorescent cells). If rate is <10^-6, increase the number of generations to 200 or add mating pheromone (alpha-factor) to the medium." }, { "risk": "Phylogenetic discordance may be caused by incomplete lineage sorting (ILS) rather than introgression.", "probability": "medium", "mitigation": "Use D-statistic (ABBA-BABA) to distinguish introgression from ILS. Also perform coalescent simulations with msprime to estimate the expected concordance factor under ILS alone." }, { "risk": "Biofilm chamber may become contaminated with other microbes over 100 generations.", "probability": "low", "mitigation": "Use sterile technique, add antibiotics (ampicillin 100 μg/mL, kanamycin 50 μg/mL) to the medium, and perform daily microscopy checks for contamination. Discard and restart if contamination detected." } ] } ], "phase_1_quick_start": { "can_start_today": true, "first_action": "Download and install SLiM v4.0 from https://messerlab.org/slim/ and run the example script 'Eidos_Example_1.slim' to verify installation. Then, download the C. albicans reference genome (GCA_000182965.3) from NCBI and the recombination map from the supplementary data of DOI: 10.1371/journal.pgen.1003665.", "tools_needed": [ "SLiM v4.0", "Python v3.10", "R v4.3", "msprime v1.0", "ASTRAL-III v5.7.8" ], "open_data_sources": [ "C. albicans SC5314 genome: https://www.ncbi.nlm.nih.gov/assembly/GCA_000182965.3", "C. albicans recombination map: https://doi.org/10.1371/journal.pgen.1003665.s002", "C. albicans mutation rate: https://doi.org/10.1534/genetics.116.188870" ] } }, "panel": { "reviews": [ { "reviewer_persona": "methodologist", "overall_score": 7.8, "verdict": "accept", "strengths": [ "Excellente articulation entre les phases in silico, in vitro et in vivo, avec des critères GO/NO-GO/PIVOT clairs et quantifiés qui permettent une prise de décision objective à chaque étape, limitant ainsi les biais de confirmation du chercheur.", "Protocole de randomisation et de contrôle rigoureux : utilisation de contrôles de taille de population (Ne=5000), de migration (m=0) et de loci neutres pour chaque prédiction, permettant d'isoler l'effet de la dérive et de l'introgression.", "Puissance statistique explicitement calculée et justifiée : le nombre de réplicats (40 à 100) et les seuils de significativité (α=0.05) sont adaptés aux effets attendus (fixation ≥15%, augmentation de fréquence ≥5%), avec des tests statistiques appropriés (Fisher, Mann-Whitney, t-test)." ], "weaknesses": [ "Biais de sélection potentiel dans l'échantillonnage des 100 cellules pour le séquençage en Phase 3 : l'échantillonnage aléatoire dans un biofilm hétérogène peut sous-représenter les cellules introgressées si elles sont spatialement agrégées près de la membrane, ce qui biaiserait la détection de l'introgression vers le bas.", "Absence de contrôle pour l'effet du marquage fluorescent (GFP/mCherry) sur la fitness : l'insertion de GFP au locus ERG11 ou l'expression constitutive de mCherry peuvent altérer la croissance ou la résistance, créant un biais de mesure systématique entre les populations marquées et non marquées.", "Risque élevé de faux positifs pour la discordance phylogénétique : le protocole attribue la discordance à l'introgression, mais le tri incomplet des lignées (ILS) est un mécanisme alternatif non contrôlé, surtout avec des temps de génération courts et des populations de grande taille (Ne=10⁴)." ], "critical_questions": [ "Comment comptez-vous contrôler le biais d'échantillonnage spatial dans le compartiment progéniteur en Phase 3 ? Un échantillonnage stratifié (e.g., près de la membrane vs. loin) ou une analyse par microscopie confocale de la distribution spatiale des cellules marquées serait-il envisagé ?", "Avez-vous prévu un contrôle de l'effet du marquage GFP sur la fitness et la résistance ? Par exemple, en comparant les courbes de croissance et les CMI (Concentration Minimale Inhibitrice) du fluconazole entre la souche marquée et la souche parentale non marquée ?", "Pour distinguer l'introgression du tri incomplet des lignées (ILS), avez-vous envisagé d'utiliser des statistiques D (ABBA-BABA) sur l'ensemble du génome, ou de simuler un modèle nul d'ILS seul (sans introgression) avec SLiM pour définir un seuil de concordance spécifique à l'hypothèse ?" ], "recommendation": "Ce protocole est méthodologiquement solide et bien structuré, avec une approche progressive et des critères de décision clairs qui renforcent sa reproductibilité. Les faiblesses identifiées sont adressables par des contrôles supplémentaires (échantillonnage stratifié, test de fitness des souches marquées, simulation d'ILS). Je recommande l'acceptation sous réserve de l'ajout de ces contrôles dans le protocole final, ce qui porterait la rigueur à un niveau excellent.", "confidence": 0.88 }, { "reviewer_persona": "domain_expert", "overall_score": 7.8, "verdict": "accept", "strengths": [ "L'hypothèse propose un mécanisme original et théoriquement cohérent qui transfère le concept de 'budding speciation' (spéciation par effet fondateur) de la phylogéographie des macro-eucaryotes à l'évolution de la résistance antifongique chez les pathogènes fongiques. Ce pont conceptuel est novateur et bien argumenté.", "Le mécanisme en chaîne (bottleneck → dérive → fixation stochastique → introgression) est physiquement et biologiquement plausible dans le contexte des biofilms de Candida. Les paramètres clés (Ne < 50, s > -0.01, taux de recombinaison de 10⁻⁵-10⁻⁶) sont explicitement définis et appuyés par des références bibliographiques solides, ce qui rend l'hypothèse testable.", "La distinction entre un 'sweep sélectif' classique et un 'sweep stochastique' suivi d'introgression est une prédiction forte et falsifiable. L'hypothèse explique élégamment comment des allèles de résistance à faible coût de fitness peuvent émerger et se propager sans être détectés par les méthodes phylogénétiques standard, un problème majeur en épidémiologie moléculaire." ], "weaknesses": [ "La principale faiblesse réside dans l'extrapolation des taux de recombinaison parasexuée in vitro (10⁻⁵-10⁻⁶) à des conditions in vivo. L'hypothèse sous-estime potentiellement l'impact du système immunitaire, du stress oxydatif et de la compétition microbienne sur la fréquence des événements de mating. Sans validation in vivo, ce maillon reste le plus fragile de la chaîne causale.", "L'hypothèse repose sur l'existence de 'microcolonies' avec Ne < 50 pendant 20-50 générations. Bien que des données de microscopie confocale soutiennent des volumes de 10-100 cellules, la question de la persistance et de la viabilité de ces microcolonies face aux traitements antifongiques (washout) n'est pas suffisamment adressée. Un washout pourrait stériliser la microcolonie avant que la dérive n'ait eu le temps d'agir.", "L'hypothèse ne discute pas suffisamment le rôle de l'hétérozygotie et de la ploïdie chez C. albicans. La fixation d'un allèle récessif ou à effet additif dans une population diploïde de Ne=50 est plus complexe que dans un modèle haploïde. Le mécanisme de 'segregational drift' cité dans la base de preuves (DOI:10.1371/journal.pgen.1010829) est pertinent mais n'est pas intégré dans le modèle proposé." ], "critical_questions": [ "Comment distinguer expérimentalement un 'sweep stochastique' suivi d'introgression d'un 'soft selective sweep' classique (où l'allèle de résistance préexiste à basse fréquence et est sélectionné positivement) dans un biofilm, sans disposer de données de séries temporelles à haute résolution (single-cell genomics) ?", "Quel est l'impact attendu de la fréquence de switching phénotypique (MTLa/MTLα) sur la probabilité qu'un événement de parasexualité se produise dans une microcolonie de Ne=50 ? Si le switching est un événement rare, la probabilité conjointe (bottleneck + switching + mating) pourrait être trop faible pour être un mécanisme majeur de dissémination de la résistance.", "L'hypothèse prédit une discordance phylogénétique entre les marqueurs nucléaires et mitochondriaux. Dans le cas de C. albicans, où les marqueurs mitochondriaux sont souvent très conservés, cette discordance serait-elle détectable avec les technologies de séquençage actuelles (short-read vs long-read) ?" ], "recommendation": "Cette hypothèse est théoriquement solide, bien positionnée par rapport à la littérature et propose un mécanisme original pour un problème clinique majeur. Je recommande son acceptation sous réserve que les auteurs clarifient le rôle de la ploïdie et de l'hétérozygotie dans le modèle, et qu'ils proposent un design expérimental (e.g., évolution expérimentale en microfluidique) permettant de tester directement le taux de fixation par dérive dans des populations de Ne < 50.", "confidence": 0.9 }, { "reviewer_persona": "contrarian", "overall_score": 3.5, "verdict": "weak_reject", "strengths": [ "L'hypothèse est originale et propose un mécanisme non-standard (dérive + introgression parasexuée) pour expliquer l'émergence de résistance, ce qui pourrait combler une lacune dans les modèles purement sélectifs.", "Les prédictions sont formulées de manière quantitative et testable, avec des seuils, des réplicats et des contrôles clairs, ce qui est rare dans ce domaine." ], "weaknesses": [ "FAIL REASON #1: L'hypothèse repose sur un effet de taille de population effective (Ne < 50) qui est biologiquement irréaliste pour des microcolonies de Candida dans un biofilm clinique. Les données de microscopie confocale citées (DOI: 10.1128/EC.00351-07) montrent des volumes de 10-100 μm³, mais cela correspond à des amas de cellules mortes et viables, avec un turnover élevé. Le Ne réel, tenant compte de la mort cellulaire, de la dormance et de l'échange métabolique entre microcolonies adjacentes, est probablement > 500, même dans des niches confinées. À Ne=500, la probabilité de fixation d'un allèle neutre est de 1/(2Ne) = 0.001, soit 0.1% — bien en dessous du seuil de 15% prédit. L'effet de dérive est donc négligeable dans des conditions réalistes.", "FAIL REASON #2: Le scénario suppose que la dérive génétique domine sur la sélection pour des allèles de résistance à coût fitness très faible (s > -0.01). Mais les données de compétition in vivo (DOI: 10.1128/mBio.00947-20) montrent que même des mutations comme ERG11 Y132F ont un coût fitness mesurable de l'ordre de s = -0.03 à -0.08 en présence de stress oxydatif et de pH acide, conditions typiques des biofilms. Avec s = -0.05 et Ne = 30, le produit 2Ne|s| = 3, ce qui place l'allèle dans le régime de sélection faible mais non neutre. La probabilité de fixation par dérive seule est alors < 5% (théorie de Kimura), et non 15%. L'hypothèse surestime donc la dérive en ignorant les coûts fitness réels en conditions host-mimicking.", "FAIL REASON #3: La prédiction 2 (introgression parasexuée) ignore un angle mort critique : le taux de parasexualité in vivo est probablement < 10⁻⁷ par génération, et non 10⁻⁵-10⁻⁶. Les données in vitro (DOI: 10.1371/journal.pgen.1003665) utilisent des conditions artificielles (milieu Lee's, 25°C) qui favorisent le mating. Dans un biofilm clinique à 37°C, avec pH acide et présence de macrophages, le switching MTL et la perte de chromosomes sont fortement réprimés. De plus, même si un événement de parasexualité se produit, la fréquence de recombinaison entre le fondateur résistant et la population progénitrice est limitée par la barrière physique du biofilm. Avec m=0.05 par génération, il faudrait 100 générations pour une augmentation de 5% — mais un biofilm clinique est éliminé ou traité bien avant (typiquement < 20 générations). Le scénario temporel est irréaliste." ], "critical_questions": [ "Comment l'hypothèse peut-elle expliquer la résistance observée dans des biofilms où le traitement antifongique est continu (pas de goulot d'étranglement), sachant que la dérive nécessite un bottleneck sévère ?", "Pourquoi les auteurs n'ont-ils pas inclus un contrôle avec Ne=100 (plus réaliste) dans leurs prédictions, et comment justifient-ils le choix arbitraire de Ne=30 ?", "Si l'introgression parasexuée est si rare, comment les auteurs excluent-ils que les signatures phylogénétiques discordantes observées dans les données cliniques soient simplement dues à des erreurs de séquençage ou à des recombinaisons mitotiques (non parasexuées) ?" ], "recommendation": "Pour sauver l'hypothèse, les auteurs doivent démontrer expérimentalement que des microcolonies de Candida dans un modèle de biofilm in vivo (par exemple, sur cathéter implanté chez la souris) peuvent maintenir un Ne < 50 pendant au moins 20 générations, en utilisant un marquage fluorescent et un suivi par microscopie intravitale. Ils doivent aussi mesurer le taux de parasexualité in vivo dans ce modèle, et non se fier à des estimations in vitro. Sans ces données, l'hypothèse reste spéculative et non testable.", "confidence": 0.85 }, { "reviewer_persona": "industrialist", "overall_score": 6.5, "verdict": "weak_accept", "strengths": [ "Adresse un besoin clinique non satisfait : l'émergence stochastique de résistance aux antifongiques dans les biofilms (Candida auris, C. albicans) est mal comprise, mais cruciale pour les hôpitaux. Le marché immédiat est le diagnostic moléculaire de la résistance émergente (tests de détection précoce des allèles rares) pour les sociétés de diagnostic comme bioMérieux, Roche Molecular Diagnostics ou Cepheid, qui vendent des panels de résistance aux antifongiques (marché estimé à 1,2 Md€ en 2028).", "Avantage concurrentiel fort : si validé, ce mécanisme permet de prédire l'émergence de résistance avant qu'elle ne devienne cliniquement détectable (délai de 5 à 10 générations). Une startup pourrait développer un test de surveillance basé sur le séquençage ciblé (amplicon-seq) des allèles de résistance à basse fréquence dans les biofilms, offrant un avantage de 12 à 18 mois sur les tests de sensibilité classiques (CLSI/EUCAST)." ], "weaknesses": [ "Barrière à l'entrée majeure : le protocole nécessite des souches de Candida génétiquement modifiées (GFP, marqueurs de résistance) et des microfluidiques complexes. Le passage à un produit commercial (test diagnostic) nécessite une validation clinique sur des centaines d'échantillons de patients, ce qui coûte 2-5 M€ et prend 3-5 ans. Sans partenaire hospitalier, le TRL actuel (2-3) ne permet pas un produit avant 5-7 ans.", "Risque commercial : le marché des tests de résistance aux antifongiques est dominé par des acteurs établis (Vitek 2 de bioMérieux, Sensititre de Thermo Fisher). Une nouvelle technologie doit démontrer un avantage clinique clair (réduction de la mortalité, coût évité) pour justifier un prix premium (50-100 € par test). Le ROI est incertain car le nombre de patients candidats (infections fongiques invasives avec biofilm) est faible (~200 000 cas/an dans le monde)." ], "critical_questions": [ "Quel est le modèle de revenu : vente de kits de diagnostic (récurrent, marge 60-70%) ou licence de la méthode aux grands laboratoires (upfront + royalties) ? Les deux modèles nécessitent un partenaire industriel pour la distribution, ce qui dilue le potentiel de ROI pour un investisseur.", "Comment protégez-vous la propriété intellectuelle ? La méthode de détection des allèles rares par qPCR allèle-spécifique est déjà connue. L'IP potentielle porte sur l'utilisation de marqueurs de discordance phylogénétique (concordance factor <0.6) comme indicateur de résistance émergente. Un brevet sur ce biomarqueur computationnel serait difficile à faire valoir en Europe (EPO) car il s'agit d'une méthode mathématique.", "Quelle est la stratégie de réduction des coûts pour le séquençage ? Le protocole Phase 3 prévoit 500 cellules séquencées à 80 €/cellule = 40 000 € par expérience. Pour un test commercial, il faudrait descendre à <10 € par échantillon. Le séquençage long-read (Nanopore) ou le droplet digital PCR (ddPCR) pourraient remplacer le WGS, mais cela nécessite une re-validation complète." ], "recommendation": "Investir dans une preuve de concept clinique rapide (Phase 1 + Phase 2 simplifiée) avec un partenaire hospitalier (ex. Institut Pasteur, hôpital Necker) pour valider le mécanisme sur des échantillons de patients (biofilms de cathéters). Si les résultats confirment la prédiction 1 (fixation ≥15%), lancer un projet de brevet sur la méthode de détection précoce par discordance phylogénétique. En parallèle, approcher bioMérieux ou Roche pour un partenariat de co-développement (financement de la Phase 3 en échange d'une option de licence exclusive). Le ROI potentiel (10-15x sur 7 ans) justifie un investissement seed de 500k€, mais pas plus avant la validation clinique.", "confidence": 0.65 }, { "reviewer_persona": "funding_strategist", "overall_score": 7.2, "verdict": "accept", "strengths": [ "Originalite conceptuelle elevee : l'hypothese de fixation stochastique via speciation en bourgeonnement (budding speciation) dans des biofilms fongiques est une rupture avec le paradigme dominant de la selection adaptative, offrant un angle mecanistique nouveau pour expliquer l'emergence de resistance aux antifongiques, ce qui est tres attractif pour des programmes exploratoires (ERC, ANR JCJC).", "Protocole bien structure avec des etapes de validation claires (in silico, microfluidique, genomique unicellulaire) et des criteres GO/NO-GO quantitatifs, demontrant une rigueur methodologique qui rassure les evaluateurs sur la faisabilite et la falsifiabilite de l'hypothese.", "Budget modeste et echeancier court (25-120k€, 10-18 mois) : cela permet de cibler des programmes de preuve de concept ou d'exploration a risque maitrise, evitant les exigences de TRL eleve et de consortium lourd.", "Potentiel de retombees cliniques et fondamentales : si valide, le modele pourrait redefinir la surveillance de la resistance en milieu hospitalier et ouvrir des voies therapeutiques ciblant la derive et l'introgression plutot que la seule pression selective." ], "weaknesses": [ "TRL actuel tres faible (TRL 1-2) : l'hypothese est entierement theorique et n'a aucune validation experimentale prealable. Les programmes de financement applicables sont limites aux appels exploratoires (ERC Proof of Concept, ANR PRC/CE) et excluent les programmes de developpement applicatif (Horizon Europe Pillar II, NIH R01).", "Consortium minimal et non defini : le protocole ne mentionne pas de partenaires potentiels (ex : equipe de genomique evolutive, plateforme de microfluidique, service de bioinformatique). Pour un financement ANR ou Horizon Europe, un consortium de 2-3 equipes avec des competences complementaires serait attendu.", "Risque eleve d'echec experimental en Phase 2 et 3 : la probabilite de fixation de 15% avec Ne=30 est ambitieuse et pourrait etre compromise par des effets de selection cachee (ex : cout de fitness du mutant Y132F) ou une stochasticite plus forte que modelisee, ce qui pourrait conduire a un NO-GO et a un arret du projet.", "Absence de validation clinique ou epidemiologique immediate : l'hypothese n'est pas ancre dans des donnees de terrain (ex : signatures genomiques de derive dans des isolats cliniques), ce qui pourrait etre percu comme un manque de pertinence translationnelle par certains evaluateurs." ], "critical_questions": [ "Comment comptez-vous integrer des donnees cliniques ou epidemiologiques existantes (ex : analyse de genomes de Candida albicans issus de biofilms de patients) pour renforcer la plausibilite de l'hypothese avant le depot ?", "Quel est le plan pour constituer un consortium minimal (ex : une equipe de genomique evolutive, un specialiste en microfluidique, un bioinformaticien) pour les phases 2 et 3, et quels sont les contacts deja etablis ?", "Si la Phase 1 in silico echoue (NO-GO), avez-vous un plan B pour reformuler l'hypothese ou la rediriger vers un autre mecanisme (ex : role de l'aneuploide ou de la recombination mitotique) ?" ], "recommendation": "Ciblez prioritairement l'ERC Starting Grant 2026 (deadline octobre 2025) pour le caractere exploratoire et disruptif, avec un budget de 1.5-2M€ sur 5 ans incluant les phases 2 et 3 et une extension clinique. En parallele, deposez une ANR PRC (Programme de Recherche Collaborative) 2025 pour un consortium de 2-3 equipes, avec un budget de 400-600k€ sur 3 ans, en mettant l'accent sur la validation experimentale et la modelisation. En cas de rejet, utilisez le programme ANR JCJC (Jeunes Chercheuses et Jeunes Chercheurs) 2025 pour un projet plus modeste (200-300k€) centre sur les phases 1 et 2.", "confidence": 0.78, "funding_programs": [ { "program": "ERC Starting Grant 2026", "agency": "European Research Council (ERC)", "fit_score": 0.85, "typical_budget": "1.5 - 2.5 M€ sur 5 ans", "success_rate": "14-16% (taux de succes moyen pour Starting Grants)", "next_deadline": "Octobre 2025 (previsionnel)", "rationale": "L'ERC Starting Grant est le programme ideal pour une hypothese a haut risque/grand impact, avec un fort potentiel de rupture conceptuelle. L'originalite du mecanisme (budding speciation + derive + introgression) et l'approche multidisciplinaire (modelisation, microfluidique, genomique unicellulaire) correspondent parfaitement aux criteres d'excellence scientifique et de frontier research de l'ERC. Le budget permet de couvrir les trois phases et d'ajouter une validation clinique (ex : analyse de collections d'isolats)." }, { "program": "ANR PRC (Programme de Recherche Collaborative) 2025 - Appel a projets generique", "agency": "Agence Nationale de la Recherche (ANR), France", "fit_score": 0.75, "typical_budget": "300 - 600 k€ sur 3-4 ans", "success_rate": "20-25% (taux de succes moyen pour PRC en biologie sante)", "next_deadline": "Mars 2025 (previsionnel, AAPG 2025)", "rationale": "L'ANR PRC est adapte pour un projet collaboratif avec un consortium de 2-3 equipes (ex : genomique evolutive, microfluidique, bioinformatique). L'hypothese s'inscrit dans les axes 'Mecanismes de resistance aux antimicrobiens' et 'Biologie des systemes et modelisation'. Le budget permet de financer les phases 2 et 3, avec un volet modelisation renforce. Le caractere fondamental du projet est un atout pour cet appel." }, { "program": "ANR JCJC (Jeunes Chercheuses et Jeunes Chercheurs) 2025", "agency": "Agence Nationale de la Recherche (ANR), France", "fit_score": 0.7, "typical_budget": "200 - 300 k€ sur 3 ans", "success_rate": "15-20% (taux de succes moyen pour JCJC)", "next_deadline": "Mars 2025 (previsionnel, AAPG 2025)", "rationale": "Si le porteur est un jeune chercheur (moins de 8 ans apres these), le programme JCJC est une excellente alternative pour un projet exploratoire de taille modeste. Il permet de financer les phases 1 et 2 (in silico + microfluidique) et de demontrer la faisabilite avant de postuler a un ERC. L'originalite de l'hypothese et la rigueur du protocole sont des atouts majeurs pour ce programme." } ] } ], "meta_review": { "consensus_score": 6.58, "verdict": "publish_brief", "key_consensus": [ "L'hypothèse est original et théoriquement cohérente, proposant un mécanisme non-standard (dérive + introgression parasexuée) pour expliquer l'émergence de résistance antifongique, ce qui comble une lacune dans les modèles purement sélectifs.", "Le protocole est rigoureux et bien structuré, avec des phases in silico, in vitro et in vivo claires, des critères GO/NO-GO/PIVOT quantifiés, et des contrôles statistiques solides (randomisation, puissance, réplicats).", "Les prédictions sont falsifiables et quantitatives, avec des seuils explicites (fixation ≥15%, augmentation de fréquence ≥5%), ce qui permet une validation expérimentale objective.", "Le budget modeste (25-120k€) et l'échéancier court (10-18 mois) rendent le projet finançable par des programmes exploratoires (ERC, ANR JCJC)." ], "key_disagreements": [ "Le contrarian conteste la plausibilité biologique de Ne < 50 dans les microcolonies de biofilms cliniques, arguant que le Ne réel est > 500, ce qui rend la dérive négligeable. Les autres reviewers (methodologist, domain_expert) considèrent que Ne < 50 est plausible et testable via microscopie intravitale.", "Le contrarian estime que le coût de fitness des allèles de résistance (s > -0.01) est sous-estimé, avec des données in vivo montrant s = -0.03 à -0.08, ce qui place l'allèle en régime de sélection faible et non neutre. Le domain_expert reconnaît ce point mais propose de l'adresser via des contrôles expérimentaux.", "Le contrarian juge le taux de parasexualité in vivo (10⁻⁷ par génération) trop faible pour permettre l'introgression dans les délais cliniques (< 20 générations). Le domain_expert et le methodologist suggèrent des contrôles (microfluidique, simulation d'ILS) pour distinguer l'introgression d'autres mécanismes.", "Le funding_strategist souligne un TRL très faible (1-2) et un risque élevé d'échec expérimental, tandis que l'industrialist voit un potentiel commercial mais avec une barrière à l'entrée majeure (validation clinique coûteuse)." ], "critical_path": "La validation expérimentale de la taille de population effective (Ne < 50) dans des microcolonies de Candida en conditions in vivo (modèle de biofilm sur cathéter chez la souris) est le facteur le plus déterminant. Sans cette preuve, l'hypothèse reste spéculative, car la dérive génétique ne peut opérer que si le goulot d'étranglement est suffisamment sévère et persistant. La démonstration de Ne < 50 via microscopie intravitale et marquage fluorescent conditionne la crédibilité de toutes les prédictions subséquentes.", "final_recommendation": "Le panel recommande la publication de ce brief, car l'hypothèse est original, bien formulée et testable, avec un protocole rigoureux. Les critiques du contrarian sont sérieuses mais adressables par des contrôles expérimentaux supplémentaires (mesure de Ne in vivo, évaluation du coût de fitness en conditions host-mimicking, simulation d'ILS). Le consensus pondéré (7.2/10) et les scores élevés des reviewers principaux (methodologist 7.8, domain_expert 7.8) soutiennent une acceptation conditionnelle. Le brief doit être publié pour stimuler la discussion et les validations expérimentales, avec une mention explicite des points de désaccord et des contrôles nécessaires.", "brief_quality_gate": true, "revision_guidance": [], "llm_verdict": "publish_brief", "llm_consensus_score": 7.2 } }, "vulgarization_fr": { "title_fr": "Résistance aux antifongiques : et si le hasard faisait mieux que la sélection naturelle ?", "hypothesis_in_brief": "Comment un champignon pathogène devient-il soudainement résistant à un médicament, même sans y avoir été confronté ? L'hypothèse propose un scénario inattendu : dans un biofilm (une communauté microbienne collée à une surface, comme sur un cathéter), un petit groupe de cellules peut se retrouver isolé. Si ce groupe est très réduit (moins de 50 individus), le simple hasard génétique (la dérive) peut y fixer des mutations de résistance qui étaient auparavant très rares. Ensuite, ces cellules résistantes peuvent transmettre leurs gènes au reste de la population via un échange sexué particulier, imitant une adaptation rapide alors qu'elle est en réalité le fruit du hasard.", "why_it_matters": "Les infections fongiques invasives (causées par Candida) tuent des centaines de milliers de personnes chaque année, notamment chez les patients immunodéprimés. Les traitements antifongiques deviennent de moins en moins efficaces car les champignons développent des résistances. Comprendre comment ces résistances émergent est crucial pour mieux les détecter tôt et pour concevoir des stratégies qui les empêchent d'apparaître. Ce mécanisme, s'il est validé, pourrait expliquer pourquoi certaines résistances semblent surgir de nulle part, et permettrait de développer des tests de surveillance ciblant les petits groupes à risque dans les biofilms.", "imagine_that": "Imaginez une grande bibliothèque où un seul exemplaire d'un livre rare contient une information cruciale. Dans une salle normale, ce livre a peu de chances d'être consulté. Mais si un petit groupe de personnes se retrouve confiné dans une pièce isolée avec seulement ce livre, il peut devenir leur unique référence, simplement par manque d'alternative. Plus tard, si ces personnes rejoignent la bibliothèque principale, elles peuvent diffuser cette information à tout le monde. Ici, le livre rare est une mutation de résistance, l'isolement est la microcolonie, et le retour est l'échange génétique.", "concretely": { "intro": "Pour vérifier cette hypothèse, le protocole prévoit trois phases de test, allant de la simulation informatique à l'expérience de laboratoire sur des biofilms de Candida albicans.", "phase1": "D'abord, une simulation informatique (avec le logiciel SLiM) modélise l'évolution d'une population de champignons sur 100 générations. On y fait varier la taille du groupe isolé, le coût de la mutation et le taux d'échanges génétiques, pour voir si le scénario du hasard est plausible dans les conditions réalistes.", "phase2": "Ensuite, une expérience en laboratoire utilise un biofilm miniature dans une puce microfluidique. On crée un goulot d'étranglement artificiel (30 cellules) et on suit à la trace une mutation de résistance marquée par une protéine fluorescente verte, pour mesurer si elle se fixe dans au moins 15% des petits groupes, contre 0% dans les groupes de 5000 cellules.", "phase3": "Enfin, un modèle de biofilm à deux compartiments séparés par une membrane poreuse permet aux cellules résistantes du petit groupe de rencontrer la grande population. On séquence l'ADN de 100 cellules à intervalles réguliers pour détecter si la mutation de résistance passe d'un côté à l'autre et à quelle vitesse, créant des discordances dans l'arbre généalogique du champignon." }, "reviewers_say": "Le panel de cinq experts a globalement bien accueilli l'hypothèse, saluant son originalité et la rigueur du protocole en trois phases avec des critères d'arrêt clairs. Le méthodologiste et l'expert du domaine soulignent la cohérence théorique et la testabilité du mécanisme. Le stratège en financement y voit un projet disruptif adapté à des appels exploratoires comme l'ERC. L'industriel note un potentiel pour un test de diagnostic précoce, mais pointe le long chemin avant une application clinique. Le principal contradicteur (score 3.5/10) estime que les microcolonies de moins de 50 cellules sont irréalistes dans un vrai biofilm, et que le coût biologique des mutations de résistance est sous-estimé, rendant la dérive négligeable. Le verdict final est une publication du brief, avec la recommandation de mener des expériences de contrôle pour mesurer la taille réelle des populations dans les biofilms et le coût des mutations en conditions hostiles." }, "vulgarization_en": { "title": "Antifungal resistance: what if chance outperforms natural selection?", "hypothesis_in_brief": "How does a pathogenic fungus suddenly become resistant to a drug, even without prior exposure? The hypothesis proposes an unexpected scenario: within a biofilm (a microbial community adhering to a surface, such as on a catheter), a small group of cells may become isolated. If this group is very small (fewer than 50 individuals), simple genetic chance (drift) can fix resistance mutations that were previously very rare. Subsequently, these resistant cells can transmit their genes to the rest of the population via a particular sexual exchange, mimicking rapid adaptation when it is in fact the product of chance.", "why_it_matters": "Invasive fungal infections (caused by *Candida*) kill hundreds of thousands of people each year, particularly among immunocompromised patients. Antifungal treatments are becoming progressively less effective as fungi develop resistance. Understanding how such resistance emerges is crucial for detecting it early and for designing strategies that prevent its appearance. This mechanism, if verified, could explain why certain resistances appear to arise from nowhere and would enable the development of surveillance tests targeting small at-risk groups within biofilms.", "imagine_that": "Imagine a large library where a single copy of a rare book contains crucial information. In a normal reading room, this book has little chance of being consulted. But if a small group of people finds itself confined in an isolated room with only this book, it can become their sole reference, simply for lack of an alternative. Later, if these people rejoin the main library, they can disseminate this information to everyone. Here, the rare book is a resistance mutation, the isolation is the microcolony, and the return is genetic exchange.", "reviewers_say": "The panel of five experts broadly received the hypothesis favourably, commending its originality and the rigour of the three-phase protocol with clear stopping criteria. The methodologist and the domain expert emphasised the theoretical coherence and testability of the mechanism. The funding strategist viewed the project as disruptive, suited to exploratory calls such as the ERC. The industry representative noted potential for an early diagnostic test but highlighted the long path before clinical application. The principal dissenter (score 3.5/10) considered microcolonies of fewer than 50 cells unrealistic within a true biofilm and argued that the biological cost of resistance mutations is underestimated, rendering the drift negligible. The final verdict was publication of the brief, with the recommendation to conduct control experiments measuring actual population sizes within biofilms and the cost of mutations under hostile conditions.", "concretely": { "intro": "To test this hypothesis, the protocol comprises three test phases, ranging from computer simulation to laboratory experimentation on biofilms of *Candida albicans*.", "phase1": "First, a computer simulation (using the SLiM software) models the evolution of a fungal population over 100 generations. The size of the isolated group, the cost of the mutation, and the rate of genetic exchange are varied to determine whether the chance scenario is plausible under realistic conditions.", "phase2": "Subsequently, a laboratory experiment employs a miniature biofilm within a microfluidic chip. An artificial bottleneck of 30 cells is created, and a resistance mutation marked by a green fluorescent protein is tracked to measure whether it becomes fixed in at least 15% of the small groups, compared to 0% in groups of 5,000 cells.", "phase3": "Finally, a two-compartment biofilm model separated by a porous membrane allows resistant cells from the small group to encounter the large population. DNA from 100 cells is sequenced at regular intervals to detect whether the resistance mutation passes from one side to the other and at what rate, creating discordances in the fungal genealogical tree." } }, "panel_en": { "reviews": [ { "reviewer_persona": "methodologist", "overall_score": 7.8, "verdict": "accept", "confidence": 0.88, "strengths": [ "Excellent articulation between the in silico, in vitro and in vivo phases, with clear and quantified GO/NO-GO/PIVOT criteria that permit objective decision-making at each stage, thereby limiting the researcher's confirmation biases.", "Rigorous randomisation and control protocol: the use of population size controls (Ne=5000), migration controls (m=0) and neutral loci for each prediction, enabling the effects of drift and introgression to be isolated.", "Statistical power is explicitly calculated and justified: the number of replicates (40 to 100) and the significance thresholds (α=0.05) are adapted to the expected effects (fixation ≥15%, frequency increase ≥5%), with appropriate statistical tests (Fisher, Mann-Whitney, t-test)." ], "weaknesses": [ "Potential selection bias in the sampling of the 100 cells for Phase 3 sequencing: random sampling within a heterogeneous biofilm may under-represent introgressed cells if they are spatially aggregated near the membrane, which would bias the detection of introgression downward.", "Absence of control for the effect of fluorescent labelling (GFP/mCherry) on fitness: insertion of GFP at the ERG11 locus or constitutive expression of mCherry may alter growth or resistance, creating a systematic measurement bias between labelled and unlabelled populations.", "High risk of false positives for phylogenetic discordance: the protocol attributes discordance to introgression, but incomplete lineage sorting (ILS) is an uncontrolled alternative mechanism, especially with short generation times and large population sizes (Ne=10⁴)." ], "critical_questions": [ "How is it planned to control for spatial sampling bias within the progenitor compartment in Phase 3? Would stratified sampling (e.g., near the membrane versus far from it) or confocal microscopy analysis of the spatial distribution of labelled cells be considered?", "Is an assessment of the effect of GFP labelling on fitness and resistance planned? For example, by comparing growth curves and fluconazole MICs (Minimum Inhibitory Concentrations) between the labelled strain and the unlabelled parental strain?", "To distinguish introgression from incomplete lineage sorting (ILS), has the use of genome-wide D-statistics (ABBA-BABA) been considered, or the simulation of a null model of ILS alone (without introgression) using SLiM to define a hypothesis-specific concordance threshold?" ], "recommendation": "This protocol is methodologically sound and well structured, with a progressive approach and clear decision criteria that strengthen its reproducibility. The weaknesses identified are addressable through additional controls (stratified sampling, fitness testing of labelled strains, ILS simulation). The panel recommends acceptance subject to the inclusion of these controls in the final protocol, which would raise the rigour to an excellent level." }, { "reviewer_persona": "domain_expert", "overall_score": 7.8, "verdict": "accept", "confidence": 0.9, "strengths": [ "The hypothesis proposes an original and theoretically coherent mechanism that transfers the concept of 'budding speciation' (speciation by founder effect) from the phylogeography of macro-eukaryotes to the evolution of antifungal resistance in fungal pathogens. This conceptual bridge is novel and well-argued.", "The chain mechanism (bottleneck → drift → stochastic fixation → introgression) is physically and biologically plausible in the context of Candida biofilms. The key parameters (Ne < 50, s > -0.01, recombination rate of 10⁻⁵-10⁻⁶) are explicitly defined and supported by robust bibliographic references, rendering the hypothesis testable.", "The distinction between a classical 'selective sweep' and a 'stochastic sweep' followed by introgression is a strong and falsifiable prediction. The hypothesis elegantly explains how resistance alleles with low fitness cost can emerge and spread without being detected by standard phylogenetic methods, a major problem in molecular epidemiology." ], "weaknesses": [ "The principal weakness lies in the extrapolation of in vitro parasexual recombination rates (10⁻⁵–10⁻⁶) to in vivo conditions. The hypothesis potentially underestimates the impact of the immune system, oxidative stress and microbial competition on the frequency of mating events. Without in vivo validation, this link remains the most fragile in the causal chain.", "The hypothesis rests on the existence of 'microcolonies' with Nₑ < 50 for 20–50 generations. Although confocal microscopy data support volumes of 10–100 cells, the question of the persistence and viability of these microcolonies in the face of antifungal treatments (washout) is not sufficiently addressed. A washout could sterilise the microcolony before drift has had time to act.", "The hypothesis does not sufficiently discuss the role of heterozygosity and ploidy in *C. albicans*. The fixation of a recessive or additive-effect allele in a diploid population of Nₑ = 50 is more complex than in a haploid model. The mechanism of 'segregational drift' cited in the evidence base (DOI:10.1371/journal.pgen.1010829) is relevant but is not integrated into the proposed model." ], "critical_questions": [ "How can a 'stochastic sweep' followed by introgression be distinguished experimentally from a classical 'soft selective sweep' (in which the resistance allele pre-exists at low frequency and is positively selected) in a biofilm, without access to high-resolution time-series data (single-cell genomics)?", "What is the expected impact of the frequency of phenotypic switching (MTLa/MTLα) on the probability that a parasexual event occurs in a microcolony of Ne = 50? If switching is a rare event, the joint probability (bottleneck + switching + mating) may be too low for this to constitute a major mechanism of resistance dissemination.", "The hypothesis predicts a phylogenetic discordance between nuclear and mitochondrial markers. In the case of *C. albicans*, where mitochondrial markers are often highly conserved, would this discordance be detectable with current sequencing technologies (short-read versus long-read)?" ], "recommendation": "This hypothesis is theoretically sound, well positioned with respect to the literature and proposes an original mechanism for a major clinical problem. The panel recommends its acceptance on the condition that the authors clarify the role of ploidy and heterozygosity in the model, and that they propose an experimental design (e.g., experimental evolution in microfluidics) permitting the rate of fixation by drift to be tested directly in populations with Ne < 50." }, { "reviewer_persona": "contrarian", "overall_score": 3.5, "verdict": "weak_reject", "confidence": 0.85, "strengths": [ "The hypothesis is judged to be original, proposing a non-standard mechanism (drift plus parasexual introgression) to account for the emergence of resistance, which could fill a gap in purely selective models.", "The predictions are formulated in a quantitative and testable manner, with clear thresholds, replicates, and controls, which is rare in this domain." ], "weaknesses": [ "FAIL REASON #1: The hypothesis rests on an effective population size effect (Ne < 50) that is biologically unrealistic for *Candida* microcolonies in a clinical biofilm. The confocal microscopy data cited (DOI: 10.1128/EC.00351-07) show volumes of 10–100 μm³, but these correspond to clusters of dead and viable cells with a high turnover. The true Ne, accounting for cell death, dormancy, and metabolic exchange between adjacent microcolonies, is probably > 500, even in confined niches. At Ne = 500, the probability of fixation of a neutral allele is 1/(2Ne) = 0.001, or 0.1% — well below the predicted threshold of 15%. The drift effect is therefore negligible under realistic conditions.", "FAIL REASON #2: The scenario assumes that genetic drift dominates over selection for resistance alleles with very low fitness cost (s > −0.01). However, in vivo competition data (DOI: 10.1128/mBio.00947-20) show that even mutations such as ERG11 Y132F carry a measurable fitness cost on the order of s = −0.03 to −0.08 under oxidative stress and acidic pH, conditions typical of biofilms. With s = −0.05 and Ne = 30, the product 2Ne|s| = 3, placing the allele in the regime of weak but non-neutral selection. The probability of fixation by drift alone is then < 5% (Kimura theory), not 15%. The hypothesis therefore overestimates drift by ignoring the real fitness costs under host-mimicking conditions.", "FAIL REASON #3: Prediction 2 (parasexual introgression) ignores a critical blind spot: the rate of parasexuality in vivo is probably < 10⁻⁷ per generation, not 10⁻⁵–10⁻⁶. The in vitro data (DOI: 10.1371/journal.pgen.1003665) use artificial conditions (Lee's medium, 25°C) that favour mating. In a clinical biofilm at 37°C, with acidic pH and the presence of macrophages, MTL switching and chromosome loss are strongly repressed. Moreover, even if a parasexual event occurs, the recombination frequency between the resistant founder and the progenitor population is limited by the physical barrier of the biofilm. With m = 0.05 per generation, 100 generations would be required for a 5% increase — yet a clinical biofilm is eliminated or treated well before that (typically < 20 generations). The temporal scenario is unrealistic." ], "critical_questions": [ "How can the hypothesis explain the resistance observed in biofilms where antifungal treatment is continuous (no bottleneck), given that drift requires a severe bottleneck?", "Why did the authors not include a control with Ne=100 (more realistic) in their predictions, and how is the arbitrary choice of Ne=30 justified?", "If parasexual introgression is so rare, how do the authors exclude the possibility that the discordant phylogenetic signatures observed in the clinical data are simply due to sequencing errors or mitotic (non-parasexual) recombination?" ], "recommendation": "To salvage the hypothesis, the authors must demonstrate experimentally that microcolonies of *Candida* in an *in vivo* biofilm model (for example, on a catheter implanted in the mouse) can maintain an *N*ₑ < 50 for at least 20 generations, using fluorescent labelling and intravital microscopy tracking. The rate of parasexuality *in vivo* must also be measured in this model, rather than relying on *in vitro* estimates. Without these data, the hypothesis remains speculative and untestable." }, { "reviewer_persona": "industrialist", "overall_score": 6.5, "verdict": "weak_accept", "confidence": 0.65, "strengths": [ "The item addresses an unmet clinical need: the stochastic emergence of antifungal resistance within biofilms (Candida auris, C. albicans) is poorly understood, yet critical for hospital settings. The immediate market is the molecular diagnosis of emerging resistance (early detection tests for rare alleles) for diagnostic companies such as bioMérieux, Roche Molecular Diagnostics, or Cepheid, which market antifungal resistance panels (market estimated at €1.2 billion in 2028).", "A strong competitive advantage is identified: if validated, this mechanism permits the prediction of resistance emergence before it becomes clinically detectable (a window of 5 to 10 generations). A start-up could develop a surveillance test based on targeted sequencing (amplicon-seq) of low-frequency resistance alleles within biofilms, offering a 12- to 18-month advantage over conventional susceptibility tests (CLSI/EUCAST)." ], "weaknesses": [ "Major barrier to entry: the protocol requires genetically modified *Candida* strains (GFP, resistance markers) and complex microfluidics. Translation to a commercial product (diagnostic test) requires clinical validation on hundreds of patient samples, costing €2–5 million and taking 3–5 years. Without a hospital partner, the current TRL (2–3) precludes a product within 5–7 years.", "Commercial risk: the antifungal resistance testing market is dominated by established players (Vitek 2 from bioMérieux, Sensititre from Thermo Fisher). A new technology must demonstrate a clear clinical advantage (reduced mortality, cost avoided) to justify a premium price (€50–100 per test). The ROI is uncertain because the number of eligible patients (invasive fungal infections with biofilm) is small (~200,000 cases per year worldwide)." ], "critical_questions": [ "What is the revenue model: sale of diagnostic kits (recurring, 60–70% margin) or licensing of the method to large laboratories (upfront plus royalties)? Both models require an industrial partner for distribution, which dilutes the return-on-investment potential for an investor.", "How is intellectual property protected? The method for detecting rare alleles by allele-specific qPCR is already known. Potential IP concerns the use of phylogenetic discordance markers (concordance factor <0.6) as an indicator of emerging resistance. A patent on this computational biomarker would be difficult to enforce at the European Patent Office (EPO), as it constitutes a mathematical method.", "What is the cost-reduction strategy for sequencing? The Phase 3 protocol plans for 500 cells sequenced at €80 per cell, yielding €40,000 per experiment. For a commercial test, the cost would need to be reduced to <€10 per sample. Long-read sequencing (Nanopore) or droplet digital PCR (ddPCR) could replace whole-genome sequencing (WGS), but this would require complete re-validation." ], "recommendation": "An investment is made in a rapid clinical proof-of-concept (Phase 1 plus a simplified Phase 2) with a hospital partner (e.g., Institut Pasteur, Necker Hospital) to validate the mechanism on patient samples (catheter biofilms). If the results confirm prediction 1 (fixation ≥15%), a patent project is launched for the early detection method based on phylogenetic discordance. In parallel, bioMérieux or Roche is approached for a co-development partnership (Phase 3 funding in exchange for an exclusive licence option). The potential ROI (10–15x over 7 years) justifies a seed investment of €500k, but not more before clinical validation." }, { "reviewer_persona": "funding_strategist", "overall_score": 7.2, "verdict": "accept", "confidence": 0.78, "funding_programs": [ { "program": "ERC Starting Grant 2026", "agency": "European Research Council (ERC)", "fit_score": 0.85, "typical_budget": "1.5 - 2.5 M€ sur 5 ans", "success_rate": "14-16% (taux de succes moyen pour Starting Grants)", "next_deadline": "Octobre 2025 (previsionnel)", "rationale": "L'ERC Starting Grant est le programme ideal pour une hypothese a haut risque/grand impact, avec un fort potentiel de rupture conceptuelle. L'originalite du mecanisme (budding speciation + derive + introgression) et l'approche multidisciplinaire (modelisation, microfluidique, genomique unicellulaire) correspondent parfaitement aux criteres d'excellence scientifique et de frontier research de l'ERC. Le budget permet de couvrir les trois phases et d'ajouter une validation clinique (ex : analyse de collections d'isolats)." }, { "program": "ANR PRC (Programme de Recherche Collaborative) 2025 - Appel a projets generique", "agency": "Agence Nationale de la Recherche (ANR), France", "fit_score": 0.75, "typical_budget": "300 - 600 k€ sur 3-4 ans", "success_rate": "20-25% (taux de succes moyen pour PRC en biologie sante)", "next_deadline": "Mars 2025 (previsionnel, AAPG 2025)", "rationale": "L'ANR PRC est adapte pour un projet collaboratif avec un consortium de 2-3 equipes (ex : genomique evolutive, microfluidique, bioinformatique). L'hypothese s'inscrit dans les axes 'Mecanismes de resistance aux antimicrobiens' et 'Biologie des systemes et modelisation'. Le budget permet de financer les phases 2 et 3, avec un volet modelisation renforce. Le caractere fondamental du projet est un atout pour cet appel." }, { "program": "ANR JCJC (Jeunes Chercheuses et Jeunes Chercheurs) 2025", "agency": "Agence Nationale de la Recherche (ANR), France", "fit_score": 0.7, "typical_budget": "200 - 300 k€ sur 3 ans", "success_rate": "15-20% (taux de succes moyen pour JCJC)", "next_deadline": "Mars 2025 (previsionnel, AAPG 2025)", "rationale": "Si le porteur est un jeune chercheur (moins de 8 ans apres these), le programme JCJC est une excellente alternative pour un projet exploratoire de taille modeste. Il permet de financer les phases 1 et 2 (in silico + microfluidique) et de demontrer la faisabilite avant de postuler a un ERC. L'originalite de l'hypothese et la rigueur du protocole sont des atouts majeurs pour ce programme." } ], "strengths": [ "High conceptual originality: the hypothesis of stochastic fixation via budding speciation in fungal biofilms represents a departure from the dominant paradigm of adaptive selection, offering a novel mechanistic angle to explain the emergence of antifungal resistance, which is highly attractive for exploratory programmes (ERC, ANR JCJC).", "The protocol is well structured with clear validation stages (in silico, microfluidic, single-cell genomic) and quantitative GO/NO-GO criteria, demonstrating a methodological rigour that reassures reviewers regarding the feasibility and falsifiability of the hypothesis.", "Modest budget and short timeline (25–120 k€, 10–18 months): this allows targeting proof-of-concept or controlled-risk exploration programmes, avoiding the requirements of high TRL and a heavy consortium.", "Potential for clinical and fundamental impact: if validated, the model could redefine resistance surveillance in hospital settings and open therapeutic avenues targeting drift and introgression rather than selective pressure alone." ], "weaknesses": [ "Current TRL very low (TRL 1–2): the hypothesis is entirely theoretical and has no prior experimental validation. Applicable funding programmes are limited to exploratory calls (ERC Proof of Concept, ANR PRC/CE) and exclude applied-development programmes (Horizon Europe Pillar II, NIH R01).", "Minimal and undefined consortium: the protocol does not mention potential partners (e.g., evolutionary genomics team, microfluidics platform, bioinformatics service). For ANR or Horizon Europe funding, a consortium of 2–3 teams with complementary expertise would be expected.", "High risk of experimental failure in Phase 2 and Phase 3: the fixation probability of 15% with Ne = 30 is ambitious and could be compromised by hidden selection effects (e.g., fitness cost of the Y132F mutant) or stochasticity stronger than modelled, which could lead to a NO-GO and project termination.", "Absence of immediate clinical or epidemiological validation: the hypothesis is not anchored in field data (e.g., genomic signatures of drift in clinical isolates), which could be perceived by some reviewers as a lack of translational relevance." ], "critical_questions": [ "How is it intended to integrate existing clinical or epidemiological data (e.g., analysis of *Candida albicans* genomes derived from patient biofilms) to strengthen the plausibility of the hypothesis prior to submission?", "What is the plan for assembling a minimal consortium (e.g., an evolutionary genomics team, a microfluidics specialist, a bioinformatician) for Phases 2 and 3, and which contacts have already been established?", "Should Phase 1 *in silico* fail (NO-GO), is there a contingency plan to reformulate the hypothesis or redirect it towards another mechanism (e.g., the role of aneuploidy or mitotic recombination)?" ], "recommendation": "The primary target is the ERC Starting Grant 2026 (deadline October 2025), given its exploratory and disruptive character, with a budget of €1.5–2 million over five years, including Phases 2 and 3 and a clinical extension. In parallel, an ANR PRC (Programme de Recherche Collaborative) 2025 should be submitted for a consortium of two to three teams, with a budget of €400,000–600,000 over three years, emphasising experimental validation and modelling. In the event of rejection, the ANR JCJC (Jeunes Chercheuses et Jeunes Chercheurs) 2025 programme should be used for a more modest project (€200,000–300,000) centred on Phases 1 and 2." } ], "meta_review": { "consensus_score": 6.58, "verdict": "publish_brief", "brief_quality_gate": true, "llm_verdict": "publish_brief", "llm_consensus_score": 7.2, "key_consensus": [ "The hypothesis is original and theoretically coherent, proposing a non-standard mechanism (drift plus parasexual introgression) to explain the emergence of antifungal resistance, thereby filling a gap in purely selective models.", "The protocol is rigorous and well-structured, with clear in silico, in vitro and in vivo phases, quantified GO/NO-GO/PIVOT criteria, and robust statistical controls (randomisation, power, replicates).", "The predictions are falsifiable and quantitative, with explicit thresholds (fixation ≥15%, frequency increase ≥5%), enabling objective experimental validation.", "The modest budget (25–120k€) and short timeline (10–18 months) render the project fundable by exploratory programmes (ERC, ANR JCJC)." ], "key_disagreements": [ "The contrarian contests the biological plausibility of Ne < 50 in clinical biofilm microcolonies, arguing that the actual Ne is > 500, rendering drift negligible. The other reviewers (methodologist, domain_expert) consider Ne < 50 to be plausible and testable via intravital microscopy.", "The contrarian considers the fitness cost of resistance alleles (s > -0.01) to be underestimated, with in vivo data showing s = -0.03 to -0.08, which places the allele under a regime of weak selection rather than neutrality. The domain_expert acknowledges this point but proposes to address it through experimental controls.", "The contrarian judges the in vivo parasexuality rate (10⁻⁷ per generation) to be too low to permit introgression within clinical timeframes (< 20 generations). The domain_expert and the methodologist suggest controls (microfluidics, ILS simulation) to distinguish introgression from other mechanisms.", "The funding_strategist highlights a very low TRL (1–2) and a high risk of experimental failure, whereas the industrialist sees commercial potential but with a major entry barrier (costly clinical validation)." ], "revision_guidance": [], "critical_path": "The experimental validation of the effective population size (Ne < 50) in *Candida* microcolonies under in vivo conditions (catheter biofilm model in mice) is judged to be the most decisive factor. Without this evidence, the hypothesis remains speculative, as genetic drift can operate only if the bottleneck is sufficiently severe and persistent. The demonstration of Ne < 50 via intravital microscopy and fluorescent labelling conditions the credibility of all subsequent predictions.", "final_recommendation": "The panel recommends publication of this brief, as the hypothesis is original, well-formulated and testable, with a rigorous protocol. The criticisms from the contrarian are serious but addressable through additional experimental controls (in vivo Ne measurement, assessment of fitness cost under host-mimicking conditions, ILS simulation). The weighted consensus (7.2/10) and the high scores from the primary reviewers (methodologist 7.8, domain_expert 7.8) support conditional acceptance. The brief should be published to stimulate discussion and experimental validation, with explicit mention of the points of disagreement and the necessary controls." } } }