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.

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.

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.

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.

Persistence of Resident and Transplanted Genotypes of the Undomesticated Yeast Saccharomyces paradoxus in Forest Soil.

One might expect yeasts in soil to be highly dispersed via water or insects, forming ephemeral, genetically heterogeneous populations subject to competition and environmental stochasticity. Here, we report persistence of genotypes of the yeast Saccharomyces paradoxus in space and time. Within 1 km2 in a mixed hardwood forest on scales from centimeters to tens of meters, we detected persistence over 3 years of native genotypes, identified by single nucleotide polymorphisms (SNPs) genome-wide, of the wild yeast Saccharomyces paradoxus growing around Quercus rubra and Quercus alba Yeasts were recovered by enrichment in ethanol-containing medium, which measures only presence or absence, not abundance. Additional transplantation experiments employed strains marked with spontaneous defects in the URA3 gene, which also confer resistance to 5-fluoroorotic acid (5FOA). Plating soil suspensions from transplant sites on 5FOA-containing medium permitted one-step quantification of yeast CFU, with no interference from other unmarked yeasts or microorganisms. After an initial steep decrease in abundance, the yeast densities fluctuated over time, increasing in association with rainfall and decreasing in association with drought. After 18 months, the transplanted yeasts remained in place on the nine sites. In vitro transplantation experiments into nonsterile soil in petri dishes showed similar patterns of persistence and response to moisture and drought. To determine whether Saccharomyces cerevisiae, not previously recovered from soils regionally, can persist in our cold climate sites, we transplanted marked S. cerevisiae alone and in mixture with S. paradoxus in the fall of 2017. Five months later, S. cerevisiae persisted to the same extent as S. paradoxusIMPORTANCESaccharomyces yeasts are intensively studied in biological research and in their domesticated roles in brewing and baking, and yet, remarkably little is known about their mode of life in forest soils. We report here that resident genotypes of the yeast S. paradoxus are persistent on a time scale of years in their microhabitats in forest soils. We also show that resident genotypes can be replaced by transplanted yeast genotypes. The high inoculum levels in experimental transplantations rapidly decreased over time, but the transplanted genotypes persisted at low abundance. We conclude that, in forest soils, Saccharomyces yeasts exist at very low abundance and that dispersal events are rare.

Liquid biopsies in lung cancer-time to implement research technologies in routine care?

Lung cancer is the leading cause of cancer mortality. A substantial progress in the understanding of lung cancer biology has resulted in several promising targeted therapies for advanced disease. Druggable targets today include point mutations such as EGFR, BRAF and re-arrangements in genes such as ALK and ROS1. Liquid biopsies collecting e.g., circulating tumor DNA (ctDNA) reflects overall tumor information and is not biased by analyzing of only a small fraction of the tumor and is always accessible in contrast to the lung cancer tissue. Technological advances in detection of low frequency mutation variants in ctDNA have made it the dominating liquid biopsy platform in terms of utility and sensitivity. Circulating DNA or RNA may possible be used to define populations with higher risk of developing lung cancer, thus reducing screening cohorts and increasing the positive predictive value of screening. Blood based-tests may also aid to identify genetic alterations several weeks prior to radiologically verified recurrence and may be of great value in the follow-up of lung cancer patients. Besides being an alternative to invasive biopsies in selected cases, liquid biopsies offer a unique possibility to monitor treatment response following medical treatment as well as treatment response and resistance development after targeted therapy, giving a possibility to modify the treatment after the genetic profile of the tumor. Ideally, genetic alterations found in ctDNA could be tracked in real-time discriminating between fast-growing life-threatening tumors from more indolent slow growing tumors or premalignant growth that are of no concern for the wellbeing of the patient. This review focuses on future perspectives of liquid biopsies in lung cancer care for different clinical settings and present current technological platforms for further discussion of possible strategies for implementation of liquid biopsies in lung cancer.

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.

Fungus Causing White-Nose Syndrome in Bats Accumulates Genetic Variability in North America with No Sign of Recombination.

Emerging fungal diseases of wildlife are on the rise worldwide, and the white-nose syndrome (WNS) epidemic in North American bats is a catastrophic example. The causal agent of WNS is a single clone of the fungus Pseudogymnoascus destructans. Early evolutionary change in this clonal population has major implications for disease ecology and conservation. Accumulation of variation in the fungus through mutation, and shuffling of variation through recombination, could affect the virulence and transmissibility of the fungus and the durability of what appears to be resistance arising in some bat populations. Our genome-wide analysis shows that the clonal population of P. destructans has expanded in size from a single genotype, has begun to accumulate variation through mutation, and presents no evidence as yet of genetic exchange among individuals. IMPORTANCE Since its discovery in 2006, the emerging infectious disease known as white-nose syndrome has killed millions of bats in North America, making it one of the most devastating wildlife epidemics in recorded history. We demonstrate that there has been as yet only spontaneous mutation across the North American population of P. destructans, and we find no indication of recombination. Thus, selective forces, which might otherwise impact pathogenic virulence, have so far had essentially no genetic variation on which to act. Our study confirmed the time of origin for the first and, thus far, only introduction of P. destructans to North America. This system provides an unprecedented opportunity to follow the evolution of a host-pathogen interaction unfolding in real time.

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.

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.

Genome Sequences of Three Phytopathogenic Species of the Magnaporthaceae Family of Fungi.

Magnaporthaceae is a family of ascomycetes that includes three fungi of great economic importance: Magnaporthe oryzae, Gaeumannomyces graminis var. tritici, and Magnaporthe poae. These three fungi cause widespread disease and loss in cereal and grass crops, including rice blast disease (M. oryzae), take-all disease in wheat and other grasses (G. graminis), and summer patch disease in turf grasses (M. poae). Here, we present the finished genome sequence for M. oryzae and draft sequences for M. poae and G. graminis var. tritici. We used multiple technologies to sequence and annotate the genomes of M. oryzae, M. poae, and G. graminis var. tritici. The M. oryzae genome is now finished to seven chromosomes whereas M. poae and G. graminis var. tritici are sequenced to 40.0× and 25.0× coverage respectively. Gene models were developed by the use of multiple computational techniques and further supported by RNAseq data. In addition, we performed preliminary analysis of genome architecture and repetitive element DNA.

Specific genomic aberrations predict survival, but low mutation rate in cancer hot spots, in clear cell renal cell carcinoma.

Detailed genetic profiling of clear cell renal cell carcinoma (ccRCC) has revealed genomic regions commonly affected by structural changes and a general genetic heterogeneity. VHL and PBRM1, both located at chromosome 3p, are 2 major genes mutated at high frequency but apart from these aberrations, the mutational landscape in ccRCC is largely undefined. Potential prognostic information given by the genomic changes appears to depend on the particular cohort studied. We analyzed a Swedish ccRCC cohort of 74 patients and found common changes (loss or gain occurring in >20% of the tumors) in 12 chromosomal regions (1p, 3p, 3q, 5q, 6q, 7p, 7q 8p, 9p, 9q, 10q, and 14q). A poor outcome was associated with gain of 7q and losses on 9p, 9q, and 14q. These aberrations were more frequent in metastasized tumors, suggesting alterations of genes important for tumor progression. Sequencing of 48 genes implicated in cancer revealed that only VHL, TP53, and PTEN were mutated at a noticeable frequency (51%, 9%, and 9%, respectively). Shorter relative telomere length (RTL) has been associated with loss of specific chromosomal regions in ccRCC tumors, but we could not verify this finding. However, a significantly lower tumor/nontumor (T/N) RTL ratio was detected for tumors with losses in 4q or 9p. In conclusion, poor outcome in ccRCC was associated with gain of 7q and loss on 9p, 9q, and 14q, whereas the mutation rate overall was low in a screen of cancer-associated genes.

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.

Genotype-phenotype correlations in Bothnia dystrophy caused by RLBP1 gene sequence variations.

PURPOSE: To evaluate phenotypes caused by different RLBP1 mutations in autosomal recessive retinitis pigmentosa of Bothnia type. METHODS: Compound heterozygotes for mutations in the RLBP1 gene [c.677T>A]+[c.700C>T] (p.M226K+p.R234W), n = 10, aged 7-84 years, and homozygotes c.677T>A (p.M226K), n = 2, aged 63 and 73 years, were studied using visual acuity (VA), low-contrast VA, visual fields (VFs) and optical coherence tomography (OCT). Retrospective VA and VFs, standardized dark adaptation and full-field electroretinograms (ERGs) were analysed and prolonged dark adaptometry and ERG (at 24 hr) were performed. RESULTS: Progressive decline of VA and VF areas was age-dependent. Retinal degenerative maculopathy, peripheral degenerative changes and retinitis punctata albescens (RPA) were present. Early retinal thinning in the central foveal, foveal (Ø 1 mm), and inner ring (Ø 3 mm) in the macular region, with homogenous, high-reflectance RPA changes, was visualized in and adjacent to the retinal pigment epithelium/choriocapillaris using OCT. Reduced dark adaptation and affected ERGs were present in all ages. Prolonged dark adaptation and ERG (at 24 hr), an increase in final threshold, and ERG rod and mixed rod/cone responses were found. CONCLUSIONS: The two RLBP1 genotypes presented a phenotypical and electrophysiological expression of progressive retinal disease similar to that previously described in homozygotes for the c.700C>T (p.R234W) RLBP1 mutation. The uniform phenotypical expression of RLBP1 mutations is relevant information for the disease and of importance in planning future treatment strategies.

Ocular phenotype of CORD5, an autosomal dominant retinal dystrophy associated with PITPNM3 p.Q626H mutation.

PURPOSE: To describe in detail the phenotype of CORD5 in two families segregating a mutation c.1878G>C (p.Q626H) in the PITPNM3 gene. METHODS: The study included 35 individuals from two different families of Swedish origin, all heterozygous for a PITPNM3 p.Q626H mutation. All participants underwent ophthalmological examination including kinetic perimetry, and in selected cases adaptometry, colour vision tests and optical coherence tomography (OCT). Electrophysiological studies were also performed. In some cases, the data were obtained from medical records. RESULTS: The majority of patients showed subnormal visual acuity and light sensitivity from childhood. Early signs of macular degeneration were also observed. There was a progressive decrease in visual acuity leading to legal blindness in early adulthood. Electrophysiological testing showed a progressive loss of photoreceptor function restricted mainly to the cones. OCT revealed decreased macular thickness with flattened and enlarged fovea. CONCLUSION: Our observations of the PITPNM3 p.Q626H mutation carriers confirm that CORD5 is a disease not to mix with other retinal degenerations mapped to 17p13. The results of our clinical evaluation so far indicate that CORD5 is characterized by predominant cone dysfunction without signs of general involvement of the retinal pigment epithelium. The rod system also seems to be unaffected. In this sense, CORD5 is different from other autosomal dominant CORDs where rod involvement is present to some degree in a late phase of the disease. Some intra- and inter-familial differences regarding the severity of the clinical picture were observed.

Evidence for a common toolbox based on necrotrophy in a fungal lineage spanning necrotrophs, biotrophs, endophytes, host generalists and specialists.

The Sclerotiniaceae (Ascomycotina, Leotiomycetes) is a relatively recently evolved lineage of necrotrophic host generalists, and necrotrophic or biotrophic host specialists, some latent or symptomless. We hypothesized that they inherited a basic toolbox of genes for plant symbiosis from their common ancestor. Maintenance and evolutionary diversification of symbiosis could require selection on toolbox genes or on timing and magnitude of gene expression. The genes studied were chosen because their products have been previously investigated as pathogenicity factors in the Sclerotiniaceae. They encode proteins associated with cell wall degradation: acid protease 1 (acp1), aspartyl protease (asps), and polygalacturonases (pg1, pg3, pg5, pg6), and the oxalic acid (OA) pathway: a zinc finger transcription factor (pac1), and oxaloacetate acetylhydrolase (oah), catalyst in OA production, essential for full symptom production in Sclerotinia sclerotiorum. Site-specific likelihood analyses provided evidence for purifying selection in all 8 pathogenicity-related genes. Consistent with an evolutionary arms race model, positive selection was detected in 5 of 8 genes. Only generalists produced large, proliferating disease lesions on excised Arabidopsis thaliana leaves and oxalic acid by 72 hours in vitro. In planta expression of oah was 10-300 times greater among the necrotrophic host generalists than necrotrophic and biotrophic host specialists; pac1 was not differentially expressed. Ability to amplify 6/8 pathogenicity related genes and produce oxalic acid in all genera are consistent with the common toolbox hypothesis for this gene sample. That our data did not distinguish biotrophs from necrotrophs is consistent with 1) a common toolbox based on necrotrophy and 2) the most conservative interpretation of the 3-locus housekeeping gene phylogeny--a baseline of necrotrophy from which forms of biotrophy emerged at least twice. Early oah overexpression likely expands the host range of necrotrophic generalists in the Sclerotiniaceae, while specialists and biotrophs deploy oah, or other as-yet-unknown toolbox genes, differently.

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.

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.