Genomic analysis across 53 canine cancer types reveals novel mutations and high clinical actionability potential

A genomic understanding of the oncogenic processes and individual variability of human cancer has steadily fueled improvement in patient outcomes over the past 20 years. Mutations within tumour tissues are routinely assessed through clinical genomic diagnostic assays by academic and commercial laboratories to facilitate diagnosis, prognosis and effective treatment stratification. The application of genomics has unveiled a wealth of mutation-based biomarkers in canine cancers, suggesting that the transformative principles that have revolutionized human cancer medicine can be brought to bear in veterinary oncology. To advance clinical genomics and genomics-guided medicine in canine oncology, we have developed and validated a canine cancer next-generation sequencing gene panel for the identification of multiple mutation types in clinical specimens. With this panel, we examined the genomic landscapes of 828 tumours from 813 dogs, spanning 53 cancer types. We identified 7856 alterations, encompassing copy number variants, single nucleotide variants, indels and internal tandem duplications. Additionally, we evaluated the clinical utility of these alterations by incorporating a biomarker framework from comprehensive curation of primary canine literature and inferences from human cancer genomic biomarker literature and clinical diagnostics. Remarkably, nearly 90% of the cases exhibited mutations with diagnostic, prognostic or therapeutic implications. Our work represents a thorough assessment of genomic landscapes in a large cohort of canine cancers, the first of its kind for its comprehensive inclusion of multiple mutation types and structured annotation of biomarkers, demonstrating the clinical potential of leveraging mutation-based biomarkers in veterinary oncology.

Novel genomic prognostic biomarkers for dogs with cancer

BACKGROUND

Growing evidence from dogs and humans supports the abundance of mutation-based biomarkers in tumors of dogs. Increasing the use of clinical genomic diagnostic testing now provides another powerful data source for biomarker discovery.

HYPOTHESIS

Analyzed clinical outcomes in dogs with cancer profiled using SearchLight DNA, a cancer gene panel for dogs, to identify mutations with prognostic value.

ANIMALS

A total of 127 cases of cancer in dogs were analyzed using SearchLight DNA and for which clinical outcome information was available.

METHODS

Clinical data points were collected by medical record review. Variables including mutated genes, mutations, signalment, and treatment were fitted using Cox proportional hazard models to identify factors associated with progression-free survival (PFS). The log-rank test was used to compare PFS between patients receiving and not receiving targeted treatment before first progression.

RESULTS

Combined genomic and outcomes analysis identified 336 unique mutations in 89 genes across 26 cancer types. Mutations in 6 genes (CCND1, CCND3, SMARCB1, FANCG, CDKN2A/B, and MSH6) were significantly associated with shorter PFS. Dogs that received targeted treatment before first progression (n = 45) experienced significantly longer PFS compared with those that did not (n = 82, P = .01). This significance held true for 29 dogs that received genomically informed targeted treatment compared with those that did not (P = .05).

CONCLUSION AND CLINICAL IMPORTANCE

We identified novel mutations with prognostic value and demonstrate the benefit of targeted treatment across multiple cancer types. These results provide clinical evidence of the potential for genomics and precision medicine in dogs with cancer.

Standing in the canine precision medicine knowledge gap: Improving annotation of canine cancer genomic biomarkers through systematic comparative analysis of human cancer mutations in COSMIC

The accrual of cancer mutation data and related functional and clinical associations have revolutionised human oncology, enabling the advancement of precision medicine and biomarker-guided clinical management. The catalogue of cancer mutations is also growing in canine cancers. However, without direct high-powered functional data in dogs, it remains challenging to interpret and utilise them in research and clinical settings. It is well-recognised that canine and human cancers share genetic, molecular and phenotypic similarities. Therefore, leveraging the massive wealth of human mutation data may help advance canine oncology. Here, we present a structured analysis of sequence conservation and conversion of human mutations to the canine genome through a 'caninisation' process. We applied this analysis to COSMIC, the Catalogue of Somatic Mutations in Cancer, the most prominent human cancer mutation database. For the project's initial phase, we focused on the subset of the COSMIC data corresponding to Cancer Gene Census (CGC) genes. A total of 670 canine orthologs were found for 721 CGC genes. In these genes, 365 K unique mutations across 160 tumour types were converted successfully to canine coordinates. We identified shared putative cancer-driving mutations, including pathogenic and hotspot mutations and mutations bearing similar biomarker associations with diagnostic, prognostic and therapeutic utility. Thus, this structured caninisation of human cancer mutations facilitates the interpretation and annotation of canine mutations and helps bridge the knowledge gap to enable canine precision medicine.

Using evolutionary constraint to define novel candidate driver genes in medulloblastoma

Current knowledge of cancer genomics remains biased against noncoding mutations. To systematically search for regulatory noncoding mutations, we assessed mutations in conserved positions in the genome under the assumption that these are more likely to be functional than mutations in positions with low conservation. To this end, we use whole-genome sequencing data from the International Cancer Genome Consortium and combined it with evolutionary constraint inferred from 240 mammals, to identify genes enriched in noncoding constraint mutations (NCCMs), mutations likely to be regulatory in nature. We compare medulloblastoma (MB), which is malignant, to pilocytic astrocytoma (PA), a primarily benign tumor, and find highly different NCCM frequencies between the two, in agreement with the fact that malignant cancers tend to have more mutations. In PA, a high NCCM frequency only affects the BRAF locus, which is the most commonly mutated gene in PA. In contrast, in MB, >500 genes have high levels of NCCMs. Intriguingly, several loci with NCCMs in MB are associated with different ages of onset, such as the HOXB cluster in young MB patients. In adult patients, NCCMs occurred in, e.g., the WASF-2/AHDC1/FGR locus. One of these NCCMs led to increased expression of the SRC kinase FGR and augmented responsiveness of MB cells to dasatinib, a SRC kinase inhibitor. Our analysis thus points to different molecular pathways in different patient groups. These newly identified putative candidate driver mutations may aid in patient stratification in MB and could be valuable for future selection of personalized treatment options.

Leveraging base-pair mammalian constraint to understand genetic variation and human disease

Thousands of genomic regions have been associated with heritable human diseases, but attempts to elucidate biological mechanisms are impeded by an inability to discern which genomic positions are functionally important. Evolutionary constraint is a powerful predictor of function, agnostic to cell type or disease mechanism. Single-base phyloP scores from 240 mammals identified 3.3% of the human genome as significantly constrained and likely functional. We compared phyloP scores to genome annotation, association studies, copy-number variation, clinical genetics findings, and cancer data. Constrained positions are enriched for variants that explain common disease heritability more than other functional annotations. Our results improve variant annotation but also highlight that the regulatory landscape of the human genome still needs to be further explored and linked to disease.

Leveraging Base Pair Mammalian Constraint to Understand Genetic Variation and Human Disease

Although thousands of genomic regions have been associated with heritable human diseases, attempts to elucidate biological mechanisms are impeded by a general inability to discern which genomic positions are functionally important. Evolutionary constraint is a powerful predictor of function that is agnostic to cell type or disease mechanism. Here, single base phyloP scores from the whole genome alignment of 240 placental mammals identified 3.5% of the human genome as significantly constrained, and likely functional. We compared these scores to large-scale genome annotation, genome-wide association studies (GWAS), copy number variation, clinical genetics findings, and cancer data sets. Evolutionarily constrained positions are enriched for variants explaining common disease heritability (more than any other functional annotation). Our results improve variant annotation but also highlight that the regulatory landscape of the human genome still needs to be further explored and linked to disease.

Genomic tumor analysis provides clinical guidance for the management of diagnostically challenging cancers in dogs

OBJECTIVE

To evaluate the diagnostic, prognostic, and therapeutic utility of a cancer genomic diagnostic assay (SearchLight DNA; Vidium Animal Health) for diagnostically ambiguous cancer cases.

ANIMALS

69 privately owned dogs with ambiguous cancer diagnoses and for which the genomic assay was performed.

PROCEDURES

Genomic assay reports generated between September 28, 2020, and July 31, 2022, for dogs with malignancy or suspected malignancy were reviewed to determine the assay's clinical utility defined as providing diagnostic clarity, prognostic information, and/or therapeutic options.

RESULTS

Genomic analysis provided diagnostic clarity in 37 of 69 cases (54%; group 1) and therapeutic and/or prognostic information in 22 of the remaining 32 cases (69%; group 2) for which the diagnosis remained elusive. Overall, the genomic assay was clinically useful in 86% (59/69) of cases.

CLINICAL RELEVANCE

To our knowledge, this was the first study to evaluate the multifaceted clinical utility of a single cancer genomic test in veterinary medicine. Study findings supported the use of tumor genomic testing for dogs with cancer, particularly those that are diagnostically ambiguous and therefore inherently challenging to manage. This evidence-driven genomic assay provided diagnostic guidance, prognostic support, and therapeutic options for most patients with an unclear cancer diagnosis that would otherwise have an unsubstantiated clinical plan. Furthermore, 38% (26/69) of samples were easily obtained aspirates. Sample factors (sample type, percentage of tumor cells, and number of mutations) did not influence diagnostic yield. Our study demonstrated the value of genomic testing for the management of canine cancer.

PIK3CA is recurrently mutated in canine mammary tumors, similarly to in human mammary neoplasia

Biological features of neoplastic disease affecting mammary gland tissue are shared between canines and humans. Research performed in either species has translational value and early phase clinical trials performed in canines with spontaneous disease could be informative for human trials. The purpose of this study was to investigate the somatic genetic aberrations occurring in canine mammary neoplasia by exome capture and next generation sequencing. Based on 55 tumor-normal pairs we identified the PIK3CA gene as the most commonly mutated gene in canine mammary tumors, with 25% of samples carrying mutations in this gene. A recurrent missense mutation was identified, p.H1047R, which is homologous to the human PIK3CA hotspot mutation found in different types of breast neoplasia. Mutations homologous to other known human mutation hotspots such as the PIK3CA p.E545K and the KRAS p.G12V/D were also identified. We identified copy number aberrations affecting important tumor suppressor and oncogenic pathways including deletions affecting the PTEN tumor suppressor gene. We suggest that activation of the KRAS or PIK3CA oncogenes or loss of the PTEN suppressor gene may be important for mammary tumor development in dogs. This data endorses the conservation of cancer across species and the validity of studying cancer in non-human species.

Whole-genome sequencing of multiple isolates of Puccinia triticina reveals asexual lineages evolving by recurrent mutations

The wheat leaf rust fungus, Puccinia triticina Erikss., is a worldwide pathogen of tetraploid durum and hexaploid wheat. Many races of P. triticina differ for virulence to specific leaf rust resistance genes and are found in most wheat-growing regions of the world. Wheat cultivars with effective leaf rust resistance exert selection pressure on P. triticina populations for virulent race types. The objectives of this study were to examine whole-genome sequence data of 121 P. triticina isolates and to gain insight into race evolution. The collection included isolates comprising of many different race phenotypes collected worldwide from common and durum wheat. One isolate from wild wheat relative Aegilops speltoides and two from Ae. cylindrica were also included for comparison. Based on 121,907 informative variants identified relative to the reference Race 1-1 genome, isolates were clustered into 11 major lineages with 100% bootstrap support. The isolates were also grouped based on variation in 1311 predicted secreted protein genes. In gene-coding regions, all groups had high ratios of nonsynonymous to synonymous mutations and nonsense to readthrough mutations. Grouping of isolates based on two main variation principle components for either genome-wide variation or variation just within the secreted protein genes, indicated similar groupings. Variants were distributed across the entire genome, not just within the secreted protein genes. Our results suggest that recurrent mutation and selection play a major role in differentiation within the clonal lineages.

Whole-genome sequencing of glioblastoma reveals enrichment of non-coding constraint mutations in known and novel genes

BACKGROUND

Glioblastoma (GBM) has one of the worst 5-year survival rates of all cancers. While genomic studies of the disease have been performed, alterations in the non-coding regulatory regions of GBM have largely remained unexplored. We apply whole-genome sequencing (WGS) to identify non-coding mutations, with regulatory potential in GBM, under the hypothesis that regions of evolutionary constraint are likely to be functional, and somatic mutations are likely more damaging than in unconstrained regions.

RESULTS

We validate our GBM cohort, finding similar copy number aberrations and mutated genes based on coding mutations as previous studies. Performing analysis on non-coding constraint mutations and their position relative to nearby genes, we find a significant enrichment of non-coding constraint mutations in the neighborhood of 78 genes that have previously been implicated in GBM. Among them, SEMA3C and DYNC1I1 show the highest frequencies of alterations, with multiple mutations overlapping transcription factor binding sites. We find that a non-coding constraint mutation in the SEMA3C promoter reduces the DNA binding capacity of the region. We also identify 1776 other genes enriched for non-coding constraint mutations with likely regulatory potential, providing additional candidate GBM genes. The mutations in the top four genes, DLX5, DLX6, FOXA1, and ISL1, are distributed over promoters, UTRs, and multiple transcription factor binding sites.

CONCLUSIONS

These results suggest that non-coding constraint mutations could play an essential role in GBM, underscoring the need to connect non-coding genomic variation to biological function and disease pathology.

Comparative Genomics Reveals Shared Mutational Landscape in Canine Hemangiosarcoma and Human Angiosarcoma

Angiosarcoma is a highly aggressive cancer of blood vessel-forming cells with few effective treatment options and high patient mortality. It is both rare and heterogenous, making large, well-powered genomic studies nearly impossible. Dogs commonly suffer from a similar cancer, called hemangiosarcoma, with breeds like the golden retriever carrying heritable genetic factors that put them at high risk. If the clinical similarity of canine hemangiosarcoma and human angiosarcoma reflects shared genomic etiology, dogs could be a critically needed model for advancing angiosarcoma research. We assessed the genomic landscape of canine hemangiosarcoma via whole-exome sequencing (47 golden retriever hemangiosarcomas) and RNA sequencing (74 hemangiosarcomas from multiple breeds). Somatic coding mutations occurred most frequently in the tumor suppressor TP53 (59.6% of cases) as well as two genes in the PI3K pathway: the oncogene PIK3CA (29.8%) and its regulatory subunit PIK3R1 (8.5%). The predominant mutational signature was the age-associated deamination of cytosine to thymine. As reported in human angiosarcoma, CDKN2A/B was recurrently deleted and VEGFA, KDR, and KIT recurrently gained. We compared the canine data to human data recently released by The Angiosarcoma Project, and found many of the same genes and pathways significantly enriched for somatic mutations, particularly in breast and visceral angiosarcomas. Canine hemangiosarcoma closely models the genomic landscape of human angiosarcoma of the breast and viscera, and is a powerful tool for investigating the pathogenesis of this devastating disease. IMPLICATIONS: We characterize the genomic landscape of canine hemangiosarcoma and demonstrate its similarity to human angiosarcoma.

SETD2 Is Recurrently Mutated in Whole-Exome Sequenced Canine Osteosarcoma

Osteosarcoma is a debilitating bone cancer that affects humans, especially children and adolescents. A homologous form of osteosarcoma spontaneously occurs in dogs, and its differential incidence observed across breeds allows for the investigation of tumor mutations in the context of multiple genetic backgrounds. Using whole-exome sequencing and dogs from three susceptible breeds (22 golden retrievers, 21 Rottweilers, and 23 greyhounds), we found that osteosarcoma tumors show a high frequency of somatic copy-number alterations (SCNA), affecting key oncogenes and tumor-suppressor genes. The across-breed results are similar to what has been observed for human osteosarcoma, but the disease frequency and somatic mutation counts vary in the three breeds. For all breeds, three mutational signatures (one of which has not been previously reported) and 11 significantly mutated genes were identified. TP53 was the most frequently altered gene (83% of dogs have either mutations or SCNA in TP53), recapitulating observations in human osteosarcoma. The second most frequently mutated gene, histone methyltransferase SETD2, has known roles in multiple cancers, but has not previously been strongly implicated in osteosarcoma. This study points to the likely importance of histone modifications in osteosarcoma and highlights the strong genetic similarities between human and dog osteosarcoma, suggesting that canine osteosarcoma may serve as an excellent model for developing treatment strategies in both species.Significance: Canine osteosarcoma genomics identify SETD2 as a possible oncogenic driver of osteosarcoma, and findings establish the canine model as a useful comparative model for the corresponding human disease. Cancer Res; 78(13); 3421-31. ©2018 AACR.

Comparative Genomics Integrated with Association Analysis Identifies Candidate Effector Genes Corresponding to Lr20 in Phenotype-Paired Puccinia triticina Isolates from Australia

Leaf rust is one of the most common and damaging diseases of wheat, and is caused by an obligate biotrophic basidiomycete, Puccinia triticina (Pt). In the present study, 20 Pt isolates from Australia, comprising 10 phenotype-matched pairs with contrasting pathogenicity for Lr20, were analyzed using whole genome sequencing. Compared to the reference genome of the American Pt isolate 1-1 BBBD Race 1, an average of 404,690 single nucleotide polymorphisms (SNPs) per isolate was found and the proportion of heterozygous SNPs was above 87% in the majority of the isolates, demonstrating a high level of polymorphism and a high rate of heterozygosity. From the genome-wide SNPs, a phylogenetic tree was inferred, which consisted of a large clade of 15 isolates representing diverse presumed clonal lineages including 14 closely related isolates and the more diverged isolate 670028, and a small clade of five isolates characterized by lower heterozygosity level. Principle component analysis detected three distinct clusters, corresponding exactly to the two major subsets of the small clade and the large clade comprising all 15 isolates without further separation of isolate 670028. While genome-wide association analysis identified 302 genes harboring at least one SNP associated with Lr20 virulence (p < 0.05), a Wilcoxon rank sum test revealed that 36 and 68 genes had significant (p < 0.05) and marginally significant (p < 0.1) differences in the counts of non-synonymous mutations between Lr20 avirulent and virulent groups, respectively. Twenty of these genes were predicted to have a signal peptide without a transmembrane segment, and hence identified as candidate effector genes corresponding to Lr20. SNP analysis also implicated the potential involvement of epigenetics and small RNA in Pt pathogenicity. Future studies are thus warranted to investigate the biological functions of the candidate effectors as well as the gene regulation mechanisms at epigenetic and post-transcription levels. Our study is the first to integrate phenotype-genotype association with effector prediction in Pt genomes, an approach that may circumvent some of the technical difficulties in working with obligate rust fungi and accelerate avirulence gene identification.

Genome Diversity, Recombination, and Virulence across the Major Lineages of Paracoccidioides

The Paracoccidioides genus includes two species of thermally dimorphic fungi that cause paracoccidioidomycosis, a neglected health-threatening human systemic mycosis endemic to Latin America. To examine the genome evolution and the diversity of Paracoccidioides spp., we conducted whole-genome sequencing of 31 isolates representing the phylogenetic, geographic, and ecological breadth of the genus. These samples included clinical, environmental and laboratory reference strains of the S1, PS2, PS3, and PS4 lineages of P. brasiliensis and also isolates of Paracoccidioides lutzii species. We completed the first annotated genome assemblies for the PS3 and PS4 lineages and found that gene order was highly conserved across the major lineages, with only a few chromosomal rearrangements. Comparing whole-genome assemblies of the major lineages with single-nucleotide polymorphisms (SNPs) predicted from the remaining 26 isolates, we identified a deep split of the S1 lineage into two clades we named S1a and S1b. We found evidence for greater genetic exchange between the S1b lineage and all other lineages; this may reflect the broad geographic range of S1b, which is often sympatric with the remaining, largely geographically isolated lineages. In addition, we found evidence of positive selection for the GP43 and PGA1 antigen genes and genes coding for other secreted proteins and proteases and lineage-specific loss-of-function mutations in cell wall and protease genes; these together may contribute to virulence and host immune response variation among natural isolates of Paracoccidioides spp. These insights into the recent evolutionary events highlight important differences between the lineages that could impact the distribution, pathogenicity, and ecology of Paracoccidioides. IMPORTANCE Characterization of genetic differences between lineages of the dimorphic human-pathogenic fungus Paracoccidioides can identify changes linked to important phenotypes and guide the development of new diagnostics and treatments. In this article, we compared genomes of 31 diverse isolates representing the major lineages of Paracoccidioides spp. and completed the first annotated genome sequences for the PS3 and PS4 lineages. We analyzed the population structure and characterized the genetic diversity among the lineages of Paracoccidioides, including a deep split of S1 into two lineages (S1a and S1b), and differentiated S1b, associated with most clinical cases, as the more highly recombining and diverse lineage. In addition, we found patterns of positive selection in surface proteins and secreted enzymes among the lineages, suggesting diversifying mechanisms of pathogenicity and adaptation across this species complex. These genetic differences suggest associations with the geographic range, pathogenicity, and ecological niches of Paracoccidioides lineages.

Genetic and phenotypic intra-species variation in Candida albicans

Candida albicans is a commensal fungus of the human gastrointestinal tract and a prevalent opportunistic pathogen. To examine diversity within this species, extensive genomic and phenotypic analyses were performed on 21 clinical C. albicans isolates. Genomic variation was evident in the form of polymorphisms, copy number variations, chromosomal inversions, subtelomeric hypervariation, loss of heterozygosity (LOH), and whole or partial chromosome aneuploidies. All 21 strains were diploid, although karyotypic changes were present in eight of the 21 isolates, with multiple strains being trisomic for Chromosome 4 or Chromosome 7. Aneuploid strains exhibited a general fitness defect relative to euploid strains when grown under replete conditions. All strains were also heterozygous, yet multiple, distinct LOH tracts were present in each isolate. Higher overall levels of genome heterozygosity correlated with faster growth rates, consistent with increased overall fitness. Genes with the highest rates of amino acid substitutions included many cell wall proteins, implicating fast evolving changes in cell adhesion and host interactions. One clinical isolate, P94015, presented several striking properties including a novel cellular phenotype, an inability to filament, drug resistance, and decreased virulence. Several of these properties were shown to be due to a homozygous nonsense mutation in the EFG1 gene. Furthermore, loss of EFG1 function resulted in increased fitness of P94015 in a commensal model of infection. Our analysis therefore reveals intra-species genetic and phenotypic differences in C. albicans and delineates a natural mutation that alters the balance between commensalism and pathogenicity.

Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement

Advances in modern sequencing technologies allow us to generate sufficient data to analyze hundreds of bacterial genomes from a single machine in a single day. This potential for sequencing massive numbers of genomes calls for fully automated methods to produce high-quality assemblies and variant calls. We introduce Pilon, a fully automated, all-in-one tool for correcting draft assemblies and calling sequence variants of multiple sizes, including very large insertions and deletions. Pilon works with many types of sequence data, but is particularly strong when supplied with paired end data from two Illumina libraries with small e.g., 180 bp and large e.g., 3–5 Kb inserts. Pilon significantly improves draft genome assemblies by correcting bases, fixing mis-assemblies and filling gaps. For both haploid and diploid genomes, Pilon produces more contiguous genomes with fewer errors, enabling identification of more biologically relevant genes. Furthermore, Pilon identifies small variants with high accuracy as compared to state-of-the-art tools and is unique in its ability to accurately identify large sequence variants including duplications and resolve large insertions. Pilon is being used to improve the assemblies of thousands of new genomes and to identify variants from thousands of clinically relevant bacterial strains. Pilon is freely available as open source software.

Comparative genomic analysis of human fungal pathogens causing paracoccidioidomycosis

Paracoccidioides is a fungal pathogen and the cause of paracoccidioidomycosis, a health-threatening human systemic mycosis endemic to Latin America. Infection by Paracoccidioides, a dimorphic fungus in the order Onygenales, is coupled with a thermally regulated transition from a soil-dwelling filamentous form to a yeast-like pathogenic form. To better understand the genetic basis of growth and pathogenicity in Paracoccidioides, we sequenced the genomes of two strains of Paracoccidioides brasiliensis (Pb03 and Pb18) and one strain of Paracoccidioides lutzii (Pb01). These genomes range in size from 29.1 Mb to 32.9 Mb and encode 7,610 to 8,130 genes. To enable genetic studies, we mapped 94% of the P. brasiliensis Pb18 assembly onto five chromosomes. We characterized gene family content across Onygenales and related fungi, and within Paracoccidioides we found expansions of the fungal-specific kinase family FunK1. Additionally, the Onygenales have lost many genes involved in carbohydrate metabolism and fewer genes involved in protein metabolism, resulting in a higher ratio of proteases to carbohydrate active enzymes in the Onygenales than their relatives. To determine if gene content correlated with growth on different substrates, we screened the non-pathogenic onygenale Uncinocarpus reesii, which has orthologs for 91% of Paracoccidioides metabolic genes, for growth on 190 carbon sources. U. reesii showed growth on a limited range of carbohydrates, primarily basic plant sugars and cell wall components; this suggests that Onygenales, including dimorphic fungi, can degrade cellulosic plant material in the soil. In addition, U. reesii grew on gelatin and a wide range of dipeptides and amino acids, indicating a preference for proteinaceous growth substrates over carbohydrates, which may enable these fungi to also degrade animal biomass. These capabilities for degrading plant and animal substrates suggest a duality in lifestyle that could enable pathogenic species of Onygenales to transfer from soil to animal hosts.

Comparative genomics reveals mobile pathogenicity chromosomes in Fusarium

Fusarium species are among the most important phytopathogenic and toxigenic fungi. To understand the molecular underpinnings of pathogenicity in the genus Fusarium, we compared the genomes of three phenotypically diverse species: Fusarium graminearum, Fusarium verticillioides and Fusarium oxysporum f. sp. lycopersici. Our analysis revealed lineage-specific (LS) genomic regions in F. oxysporum that include four entire chromosomes and account for more than one-quarter of the genome. LS regions are rich in transposons and genes with distinct evolutionary profiles but related to pathogenicity, indicative of horizontal acquisition. Experimentally, we demonstrate the transfer of two LS chromosomes between strains of F. oxysporum, converting a non-pathogenic strain into a pathogen. Transfer of LS chromosomes between otherwise genetically isolated strains explains the polyphyletic origin of host specificity and the emergence of new pathogenic lineages in F. oxysporum. These findings put the evolution of fungal pathogenicity into a new perspective.