DNA methylation associates with survival in non-metastatic clear cell renal cell carcinoma.

BACKGROUND: Clear cell renal cell carcinoma (ccRCC) is the most common subtype among renal cancer and is associated with poor prognosis if metastasized. Up to one third of patients with local disease at diagnosis will develop metastasis after nephrectomy, and there is a need for new molecular markers to identify patients with high risk of tumor progression. In the present study, we performed genome-wide promoter DNA methylation analysis at diagnosis to identify DNA methylation profiles associated with risk for progress. METHOD: Diagnostic tissue samples from 115 ccRCC patients were analysed by Illumina HumanMethylation450K arrays and methylation status of 155,931 promoter associated CpGs were related to genetic aberrations, gene expression and clinicopathological parameters. RESULTS: The ccRCC samples separated into two clusters (cluster A/B) based on genome-wide promoter methylation status. The samples in these clusters differed in tumor diameter (p < 0.001), TNM stage (p < 0.001), morphological grade (p < 0.001), and patients outcome (5 year cancer specific survival (pCSS5yr) p < 0.001 and cumulative incidence of progress (pCIP5yr) p < 0.001. An integrated genomic and epigenomic analysis in the ccRCCs, revealed significant correlations between the total number of genetic aberrations and total number of hypermethylated CpGs (R = 0.435, p < 0.001), and predicted mitotic age (R = 0.407, p < 0.001). We identified a promoter methylation classifier (PMC) panel consisting of 172 differently methylated CpGs accompanying progress of disease. Classifying non-metastatic patients using the PMC panel showed that PMC high tumors had a worse prognosis compared with the PMC low tumors (pCIP5yr 38% vs. 8%, p = 0.001), which was confirmed in non-metastatic ccRCCs in the publically available TCGA-KIRC dataset (pCIP5yr 39% vs. 16%, p < 0.001). CONCLUSION: DNA methylation analysis at diagnosis in ccRCC has the potential to improve outcome-prediction in non-metastatic patients at diagnosis.

Population genomics reveals structure at the individual, host-tree scale and persistence of genotypic variants of the undomesticated yeast Saccharomyces paradoxus in a natural woodland.

Genetic diversity in experimental, domesticated and wild populations of the related yeasts, Saccharomyces cerevisiae and Saccharomyces paradoxus, has been well described at the global scale. We investigated the population genomics of a local population on a small spatial scale to address two main questions. First, is there genomic variation in a S. paradoxus population at a spatial scale spanning centimetres (microsites) to tens of metres? Second, does the distribution of genomic variants persist over time? Our sample consisted of 42 S. paradoxus strains from 2014 and 43 strains from 2015 collected from the same 72 microsites around four host trees (Quercus rubra and Quercus alba) within 1 km2 in a mixed hardwood forest in southern Ontario. Six additional S. paradoxus strains recovered from adjacent maple and beech trees in 2015 are also included in the sample. Whole-genome sequencing and genomic SNP analysis revealed five differentiated groups (clades) within the sampled area. The signal of persistence of genotypes in their microsites from 2014 to 2015 was highly significant. Isolates from the same tree tended to be more related than strains from different trees, with limited evidence of dispersal between trees. In growth assays, one genotype had a significantly longer lag phase than the other strains. Our results indicate that different clades coexist at fine spatial scale and that population structure persists over at least a one-year interval in these wild yeasts, suggesting the efficacy of yearly sampling to follow longer term genetic dynamics in future studies.

Richard Paul Korf (1925-2016).

"With poetry, the tune is in the words themselves-and once you begin to hear it, it will stay with you." Richard P. Korf, notes to his narration of John Brown's Body.

DNA methylation status defines clinicopathological parameters including survival for patients with clear cell renal cell carcinoma (ccRCC).

Epigenetic alterations in the methylome have been associated with tumor development and progression in renal cell carcinoma (RCC). In this study, 45 tumor samples, 12 tumor-free kidney cortex tissues, and 24 peripheral blood samples from patients with clear cell RCC (ccRCC) were analyzed by genome-wide promoter-directed methylation arrays and related to clinicopathological parameters. Unsupervised hierarchical clustering separated the tumors into two distinct methylation groups (clusters A and B), where cluster B had higher average methylation and increased number of hypermethylated CpG sites (CpGs). Furthermore, tumors in cluster B had, compared with cluster A, a larger tumor diameter (p = 0.033), a higher morphologic grade (p < 0.001), a higher tumor-node-metastasis (TNM) stage (p < 0.001), and a worse prognosis (p = 0.005). Higher TNM stage was correlated to an increase in average methylation level (p = 0.003) and number of hypermethylated CpGs (p = 0.003), whereas a number of hypomethylated CpGs were mainly unchanged. However, the predicted age of the tumors based on methylation profile did not correlate with TNM stage, morphological grade, or methylation cluster. Differently methylated (DM) genes (n = 840) in ccRCC samples compared with tumor-free kidney cortex samples were predominantly hypermethylated and a high proportion were identified as polycomb target genes. The DM genes were overrepresented by transcription factors, ligands, and receptors, indicating functional alterations of significance for ccRCC progression. To conclude, increased number of hypermethylated genes was associated with increased TNM stage of the tumors. DNA methylation classification of ccRCC tumor samples at diagnosis can serve as a clinically applicable prognostic marker in ccRCC.

The underlying structure of adaptation under strong selection in 12 experimental yeast populations.

The aims of this study were to determine (i) whether adaptation under strong selection occurred through mutations in a narrow target of one or a few nucleotide sites or a broad target of numerous sites and (ii) whether the programs of adaptation previously observed from three experimental populations were unique or shared among populations that underwent parallel evolution. We used archived population samples from a previous study, representing 500 generations of experimental evolution in 12 populations under strong selection, 6 populations in a high-salt environment and 6 populations in a low-glucose environment. Each set of six populations included four with sexual reproduction and two with exclusively asexual reproduction. Populations were sampled as resequenced genomes of 115 individuals and as bulk samples from which frequencies of mutant alleles were estimated. In a high-salt environment, a broad target of 11 mutations within the proton exporter, PMA1, was observed among the six populations, in addition to expansions of the ENA gene cluster. This pattern was shared among populations that underwent parallel evolution. In a low-glucose environment, two programs of adaptation were observed. The originally observed pattern of mutation in MDS3/MKT1 in population M8 was a narrow target of a single nucleotide, unique to this population. Among the other five populations, the three mutations were shared in a broad target, sensing/signaling genes RAS1 and RAS2. RAS1/RAS2 mutations were not observed in the high-salt populations; PMA1 mutations were observed only in a high-salt environment.

Genomic analysis of the necrotrophic fungal pathogens Sclerotinia sclerotiorum and Botrytis cinerea.

Sclerotinia sclerotiorum and Botrytis cinerea are closely related necrotrophic plant pathogenic fungi notable for their wide host ranges and environmental persistence. These attributes have made these species models for understanding the complexity of necrotrophic, broad host-range pathogenicity. Despite their similarities, the two species differ in mating behaviour and the ability to produce asexual spores. We have sequenced the genomes of one strain of S. sclerotiorum and two strains of B. cinerea. The comparative analysis of these genomes relative to one another and to other sequenced fungal genomes is provided here. Their 38-39 Mb genomes include 11,860-14,270 predicted genes, which share 83% amino acid identity on average between the two species. We have mapped the S. sclerotiorum assembly to 16 chromosomes and found large-scale co-linearity with the B. cinerea genomes. Seven percent of the S. sclerotiorum genome comprises transposable elements compared to <1% of B. cinerea. The arsenal of genes associated with necrotrophic processes is similar between the species, including genes involved in plant cell wall degradation and oxalic acid production. Analysis of secondary metabolism gene clusters revealed an expansion in number and diversity of B. cinerea-specific secondary metabolites relative to S. sclerotiorum. The potential diversity in secondary metabolism might be involved in adaptation to specific ecological niches. Comparative genome analysis revealed the basis of differing sexual mating compatibility systems between S. sclerotiorum and B. cinerea. The organization of the mating-type loci differs, and their structures provide evidence for the evolution of heterothallism from homothallism. These data shed light on the evolutionary and mechanistic bases of the genetically complex traits of necrotrophic pathogenicity and sexual mating. This resource should facilitate the functional studies designed to better understand what makes these fungi such successful and persistent pathogens of agronomic crops.

Incipient speciation by divergent adaptation and antagonistic epistasis in yeast.

Establishing the conditions that promote the evolution of reproductive isolation and speciation has long been a goal in evolutionary biology. In ecological speciation, reproductive isolation between populations evolves as a by-product of divergent selection and the resulting environment-specific adaptations. The leading genetic model of reproductive isolation predicts that hybrid inferiority is caused by antagonistic epistasis between incompatible alleles at interacting loci. The fundamental link between divergent adaptation and reproductive isolation through genetic incompatibilities has been predicted, but has not been directly demonstrated experimentally. Here we empirically tested key predictions of speciation theory by evolving the initial stages of speciation in experimental populations of the yeast Saccharomyces cerevisiae. After replicate populations adapted to two divergent environments, we consistently observed the evolution of two forms of postzygotic isolation in hybrids: reduced rate of mitotic reproduction and reduced efficiency of meiotic reproduction. This divergent selection resulted in greater reproductive isolation than parallel selection, as predicted by the ecological speciation theory. Our experimental system allowed controlled comparison of the relative importance of ecological and genetic isolation, and we demonstrated that hybrid inferiority can be ecological and/or genetic in basis. Overall, our results show that adaptation to divergent environments promotes the evolution of reproductive isolation through antagonistic epistasis, providing evidence of a plausible common avenue to speciation and adaptive radiation in nature.

PRP8 inteins in species of the genus Botrytis and other ascomycetes.

The mobile elements termed inteins have a sporadic distribution in microorganisms. It is unclear how these elements are maintained. Inteins are intervening protein sequences that autocatalytically excise themselves from a precursor. Excision is a post-translational process referred to as 'protein splicing' in which the sequences flanking the intein are ligated, reforming the mature host protein. Some inteins contain a homing endonuclease domain (HEG) that is proposed to facilitate propagation of the intein element within a gene pool. We have previously demonstrated that the HEG of the PRP8 intein is highly active during meiosis in Botrytis cinerea. Here we analysed the Prp8 gene status in 21 additional Botrytis species to obtain insight into the mode of intein inheritance within the Botrytis lineage. Of the 21 species, 15 contained a PRP8 intein whereas six did not. The analysis was extended to closely related (Sclerotiniaceae) and distantly related (Ascomycota) taxa, focussing on evolutionary diversification of the PRP8 intein, including their possible acquisition by horizontal transfer and loss by deletion. Evidence was obtained for the occurrence of genetic footprints of previous intein occupation. There is no compelling evidence of horizontal transfer among species. Three distinct states of the Prp8 allele were identified, distributed over different orders within the Ascomycota: an occupied allele; an empty allele that was never occupied; an empty allele that was presumably previously occupied, from which the intein was precisely deleted. The presence of the genetic footprint identifies 20 species (including Neurospora crassa, Magnaporthe oryzae and Fusarium oxysporum) that previously contained the intein but have lost it entirely, while only 18 species (including Podospora anserina and Fusarium graminearum) appear never to have contained a PRP8 intein. The analysis indicates that inteins may be maintained in an equilibrium state.

Cellular effects and epistasis among three determinants of adaptation in experimental populations of Saccharomyces cerevisiae.

Epistatic interactions in which the phenotypic effect of an allele is conditional on its genetic background have been shown to play a central part in various evolutionary processes. In a previous study (J. B. Anderson et al., Curr. Biol. 20:1383-1388, 2010; J. R. Dettman, C. Sirjusingh, L. M. Kohn, and J. B. Anderson, Nature 447:585-588, 2007), beginning with a common ancestor, we identified three determinants of fitness as mutant alleles (each designated with the letter "e") that arose in replicate Saccharomyces cerevisiae populations propagated in two different environments, a low-glucose and a high-salt environment. In a low-glucose environment, MDS3e and MKT1e interacted positively to confer a fitness advantage. Also, PMA1e from a high-salt environment interacted negatively with MKT1e in a low-glucose environment, an example of a Dobzhansky-Muller incompatibility that confers reproductive isolation. Here we showed that the negative interaction between PMA1e and MKT1e is mediated by alterations in intracellular pH, while the positive interaction between MDS3e and MKT1e is mediated by changes in gene expression affecting glucose transporter genes. We specifically addressed the evolutionary significance of the positive interaction by showing that the presence of the MDS3 mutation is a necessary condition for the spread and fixation of the new mutations at the identical site in MKT1. The expected mutations in MKT1 rose to high frequencies in two of three experimental populations carrying MDS3e but not in any of three populations carrying the ancestral allele. These data show how positive and negative epistasis can contribute to adaptation and reproductive isolation.

Clinical and biological relevance of the transcriptomic-based prostate cancer metastasis subtypes MetA-C.

To improve treatment of metastatic prostate cancer, the biology of metastases needs to be understood. We recently described three subtypes of prostate cancer bone metastases (MetA-C), based on differential gene expression. The aim of this study was to verify the clinical relevance of these subtypes and to explore their biology and relations to genetic drivers. Freshly-frozen metastasis samples were obtained as hormone-naive (n = 17), short-term castrated (n = 21), or castration-resistant (n = 65) from a total of 67 patients. Previously published sequencing data from 573 metastasis samples were also analyzed. Through transcriptome profiling and sample classification based on a set of predefined MetA-C-differentiating genes, we found that most metastases were heterogeneous for the MetA-C subtypes. Overall, MetA was the most common subtype, while MetB was significantly enriched in castration-resistant samples and in liver metastases, and consistently associated with poor prognosis. By gene set enrichment analysis, the phenotype of MetA was described by high androgen response, protein secretion and adipogenesis, MetB by high cell cycle activity and DNA repair, and MetC by epithelial-to-mesenchymal transition and inflammation. The MetB subtype demonstrated single nucleotide variants of RB transcriptional corepressor 1 (RB1) and loss of 21 genes at chromosome 13, including RB1, but provided independent prognostic value to those genetic aberrations. In conclusion, a distinct set of gene transcripts can be used to classify prostate cancer metastases into the subtypes MetA-C. The MetA-C subtypes show diverse biology, organ tropism, and prognosis. The MetA-C classification may be used independently, or in combination with genetic markers, primarily to identify MetB patients in need of complementary therapy to conventional androgen receptor-targeting treatments.

Mutation spectra in autosomal dominant and recessive retinitis pigmentosa in northern Sweden.

Retinal degenerations represent a heterogeneous group of disorders affecting the function of the retina. The frequency of retinitis pigmentosa (RP) is 1/3500 worldwide, however, in northern Sweden it is 1/2000 due to limited migration and a 'founder' effect. In this study we identified genetic mechanisms underlying autosomal dominant and recessive RP present in northern Sweden. Several novel mutations unique for this region were found. In an autosomal recessive form of RP, Bothnia dystrophy caused by mutations in the RLBP1 gene, bi-allelic mutations R234W, M226K and compound heterozygosity, M226K+R234W was detected.In dominant form of RP mapped to 19q13.42 a 59 kb genomic deletion including the PRPF31 and three other genes was found.These data provide additional information on the molecular mechanisms of RP evolvement and in the future might be useful in development of therapeutic strategies. Identification of the disease-causing mutations allowed introducing molecular genetic testing of the patients and their families into the clinical practice.

Single nucleotide polymorphism-based diagnostic system for crop-associated Sclerotinia species.

A molecular diagnostic system using single nucleotide polymorphisms (SNPs) was developed to identify four Sclerotinia species: S. sclerotiorum (Lib.) de Bary, S. minor Jagger, S. trifoliorum Erikss., and the undescribed species Sclerotinia species 1. DNAs of samples are hybridized with each of five 15-bp oligonucleotide probes containing an SNP site midsequence unique to each species. For additional verification, hybridizations were performed using diagnostic single nucleotide substitutions at a 17-bp sequence of the calmodulin locus. The accuracy of these procedures was compared to that of a restriction fragment length polymorphism (RFLP) method based on Southern hybridizations of EcoRI-digested genomic DNA probed with the ribosomal DNA-containing plasmid probe pMF2, previously shown to differentiate S. sclerotiorum, S. minor, and S. trifoliorum. The efficiency of the SNP-based assay as a diagnostic test was evaluated in a blind screening of 48 Sclerotinia isolates from agricultural and wild hosts. One isolate of Botrytis cinerea was used as a negative control. The SNP-based assay accurately identified 96% of Sclerotinia isolates and could be performed faster than RFLP profiling using pMF2. This method shows promise for accurate, high-throughput species identification.

Evidence for alteration of fungal endophyte community assembly by host defense compounds.

* Plant defense compounds are common stressors encountered by endophytes. Fungi readily evolve tolerance to these compounds, yet few studies have addressed the influence of intraspecific variation in defense compound production on endophyte colonization. We compared the influence of defense compound production on the composition of fungal endophyte communities in replicated field experiments. * Maize (Zea mays) produces benzoxazinoids (BXs), compounds with antifungal byproducts persistent in the environment. Fungi were isolated from leaf and root tissue of two maize genotypes that produce BXs, and a natural mutant that does not. Isolates representing the species recovered were tested for tolerance to 2-benzoxazolinone (BOA), a toxic BX byproduct. * In seedling roots and mature leaves, the community proportion with low BOA tolerance was significantly greater in BX nonproducers than producers. Mean isolation frequency of Fusarium species was up to 35 times higher in mature leaves of BX producers than nonproducers. * Fungal species with relatively high tolerance to BOA are more abundant in BX producing than BX nonproducing maize. Production of BXs may increase colonization by Fusarium species in maize, including agents of animal toxicosis and yield-reducing disease in maize. Overall, results indicate that production of defense compounds can significantly alter endophyte community assembly.

Author Correction: Hybridization is a recurrent evolutionary stimulus in wild yeast speciation.

The original version of the Supplementary Information associated with this Article contained errors in Supplementary Figures 2, 12, 20 and 22. The HTML has been updated to include a corrected version of the Supplementary Information; the original incorrect versions of these Figures can be found as Supplementary Information associated with this Correction.

Hybridization is a recurrent evolutionary stimulus in wild yeast speciation.

Hybridization can result in reproductively isolated and phenotypically distinct lineages that evolve as independent hybrid species. How frequently hybridization leads to speciation remains largely unknown. Here we examine the potential recurrence of hybrid speciation in the wild yeast Saccharomyces paradoxus in North America, which comprises two endemic lineages SpB and SpC, and an incipient hybrid species, SpC*. Using whole-genome sequences from more than 300 strains, we uncover the hybrid origin of another group, SpD, that emerged from hybridization between SpC* and one of its parental species, the widespread SpB. We show that SpD has the potential to evolve as a novel hybrid species, because it displays phenotypic novelties that include an intermediate transcriptome profile, and partial reproductive isolation with its most abundant sympatric parental species, SpB. Our findings show that repetitive cycles of divergence and hybridization quickly generate diversity and reproductive isolation, providing the raw material for speciation by hybridization.

An integrated transcriptome analysis in T-cell acute lymphoblastic leukemia links DNA methylation subgroups to dysregulated TAL1 and ANTP homeobox gene expression.

Classification of pediatric T-cell acute lymphoblastic leukemia (T-ALL) patients into CIMP (CpG Island Methylator Phenotype) subgroups has the potential to improve current risk stratification. To investigate the biology behind these CIMP subgroups, diagnostic samples from Nordic pediatric T-ALL patients were characterized by genome-wide methylation arrays, followed by targeted exome sequencing, telomere length measurement, and RNA sequencing. The CIMP subgroups did not correlate significantly with variations in epigenetic regulators. However, the CIMP+ subgroup, associated with better prognosis, showed indicators of longer replicative history, including shorter telomere length (P = 0.015) and older epigenetic (P < 0.001) and mitotic age (P < 0.001). Moreover, the CIMP+ subgroup had significantly higher expression of ANTP homeobox oncogenes, namely TLX3, HOXA9, HOXA10, and NKX2-1, and novel genes in T-ALL biology including PLCB4, PLXND1, and MYO18B. The CIMP- subgroup, with worse prognosis, was associated with higher expression of TAL1 along with frequent STIL-TAL1 fusions (2/40 in CIMP+ vs 11/24 in CIMP-), as well as stronger expression of BEX1. Altogether, our findings suggest different routes for leukemogenic transformation in the T-ALL CIMP subgroups, indicated by different replicative histories and distinct methylomic and transcriptomic profiles. These novel findings can lead to new therapeutic strategies.

Divergent adaptation promotes reproductive isolation among experimental populations of the filamentous fungus Neurospora.

BACKGROUND: An open, focal issue in evolutionary biology is how reproductive isolation and speciation are initiated; elucidation of mechanisms with empirical evidence has lagged behind theory. Under ecological speciation, reproductive isolation between populations is predicted to evolve incidentally as a by-product of adaptation to divergent environments. The increased genetic diversity associated with interspecific hybridization has also been theorized to promote the development of reproductive isolation among independent populations. Using the fungal model Neurospora, we founded experimental lineages from both intra- and interspecific crosses, and evolved them in one of two sub-optimal, selective environments. We then measured the influence that initial genetic diversity and the direction of selection (parallel versus divergent) had on the evolution of reproductive isolation. RESULTS: When assayed in the selective environment in which they were evolved, lineages typically had greater asexual fitness than the progenitors and the lineages that were evolved in the alternate, selective environment. Assays for reproductive isolation showed that matings between lineages that were adapted to the same environment had greater sexual reproductive success than matings between lineages that were adapted to different environments. Evidence of this differential reproductive success was observed at two stages of the sexual cycle. For one of the two observed incompatibility phenotypes, results from genetic analyses were consistent with a two-locus, two-allele model with asymmetric (gender-specific), antagonistic epistasis. The effects of divergent adaptation on reproductive isolation were more pronounced for populations with greater initial genetic variation. CONCLUSION: Divergent selection resulted in divergent adaptation and environmental specialization, consistent with fixation of different alleles in different environments. When brought together by mating, these alleles interacted negatively and had detrimental effects on sexual reproductive success, in agreement with the Dobzhansky-Muller model of genetic incompatibilities. As predicted by ecological speciation, greater reproductive isolation was observed among divergent-adapted lineages than among parallel-adapted lineages. These results support that, given adequate standing genetic variation, divergent adaptation can indirectly cause the evolution of reproductive isolation, and eventually lead to speciation.

The importance of patient education in the treatment of chronic kidney disease.

The article by Finkelstein and colleagues is a seminal work in helping us to understand the educational needs of pre-end-stage renal disease (pre-ESRD) patients. There are currently major deficits in patients' perceived knowledge of their options for ESRD care. The Canadian system, with its more integrated and multidisciplinary approach to care, does a better job of communicating these options to patients. African Americans are less likely to be aware of their options.

Marker stability throughout 400 days of in vitro hyphal growth in the filamentous ascomycete, Sclerotinia sclerotiorum.

The stability of routinely used, population genetic markers through approximately 1 year of continuous laboratory growth was investigated in the common, plant pathogentic ascomycete Sclerotinia sclerotiorum. Given reports of accelerated mutation rates at higher temperatures, both a permissive temperature, 22 degrees C, and a temperature at the high end of tolerance, 30 degrees C, were employed. Because mycelial growth rate was tracked among mitotic lineages established for each strain, a subsidiary objective was addressed, testing the stability of a 30 degrees C-competent phenotype. Twelve laboratory strains of S. sclerotiorum, including the genome sequence isolate, 1980, were propagated serially for up to 400 days at 22 degrees C. Five of these strains were also propagated at 30 degrees C. No mutations were observed in mycelial compatibility groupings (MCGs), DNA fingerprints, alleles at 7 microsatellite loci, or alleles at 56 AFLP loci. All of these markers show variation in field populations, which are likely much larger and influenced by different and more stochastic environmental processes. In S. sclerotiorum, population genetic markers were stable over time through serial transfer and growth of laboratory strains at both 22 degrees C and 30 degrees C. The strain isolated after extended drought and capable of infecting plants at 28 degrees C demonstrated the stability of its high temperature-competent phenotype, in addition to its stable growth rate at 22 degrees C. This observation has implications for modeling pathogen tolerance or adaptation under conditions of environmental stochasticity, including climate warming.

Host-synthesized secondary compounds influence the in vitro interactions between fungal endophytes of maize.

Maize produces a suite of allelopathic secondary metabolites, the benzoxazinoids. 2,4-Dihydroxy-7-methoxy-2H-1,4-benzoxazin-3-one and 2,4-dihydroxy-2H-1,4-benzoxazin-3-one reside as glucosides in plant tissue and spontaneously degrade to 6-methoxy-2-benzoxazolinone (MBOA) and 2-benzoxazolinone (BOA) upon plant cell disruption. Several maize-associated fungi in the genus Fusarium can metabolize MBOA and BOA. BOA tolerance levels in 10 species of Fusarium and in the maize endophytes Nigrospora oryzae, Acremonium zeae, and Periconia macrospinosa were characterized. BOA tolerance ranged from 0.25 to 1.10 mg/ml among species. The influence of substrate alteration by one species on the subsequent growth of another species was assessed in the presence and absence of BOA. The colony area of the secondary colonizer in heterospecific interactions was compared to that in autospecific interactions (one isolate follows itself). In the presence of BOA, four of six secondary colonizers had greater growth (facilitation) when primary colonizers had higher BOA tolerance than the secondary colonizer. When the primary colonizer had lower tolerance than the secondary, three of six secondary colonizers were inhibited (competition) and three not significantly affected. In BOA-free medium, the number of isolates that were facilitated or inhibited was the same regardless of the tolerance level of the primary colonizer. Two of six secondary colonizers were facilitated, two inhibited, and two not significantly affected. This study provides some support for facilitation in stressful conditions under the Menge-Sutherland model. The results are not consistent with the corresponding prediction of competition in the absence of stress. The hypothesis drawn from these data is that in the presence of a toxin, fungal species that detoxify their substrate can enhance the colonization rate of less tolerant fungi.