The order Chaetothyriales (Pezizomycotina, Ascomycetes) harbours obligatorily melanised fungi and includes numerous etiologic agents of chromoblastomycosis, phaeohyphomycosis and other diseases of vertebrate hosts. Diseases range from mild cutaneous to fatal cerebral or disseminated infections and affect humans and cold-blooded animals globally. In addition, Chaetothyriales comprise species with aquatic, rock-inhabiting, ant-associated, and mycoparasitic life-styles, as well as species that tolerate toxic compounds, suggesting a high degree of versatile extremotolerance. To understand their biology and divergent niche occupation, we sequenced and annotated a set of 23 genomes of main the human opportunists within the Chaetothyriales as well as related environmental species. Our analyses included fungi with diverse life-styles, namely opportunistic pathogens and closely related saprobes, to identify genomic adaptations related to pathogenesis. Furthermore, ecological preferences of Chaetothyriales were analysed, in conjuncture with the order-level phylogeny based on conserved ribosomal genes. General characteristics, phylogenomic relationships, transposable elements, sex-related genes, protein family evolution, genes related to protein degradation (MEROPS), carbohydrate-active enzymes (CAZymes), melanin synthesis and secondary metabolism were investigated and compared between species. Genome assemblies varied from 25.81 Mb (Capronia coronata) to 43.03 Mb (Cladophialophora immunda). The bantiana-clade contained the highest number of predicted genes (12817 on average) as well as larger genomes. We found a low content of mobile elements, with DNA transposons from Tc1/Mariner superfamily being the most abundant across analysed species. Additionally, we identified a reduction of carbohydrate degrading enzymes, specifically many of the Glycosyl Hydrolase (GH) class, while most of the Pectin Lyase (PL) genes were lost in etiological agents of chromoblastomycosis and phaeohyphomycosis. An expansion was found in protein degrading peptidase enzyme families S12 (serine-type D-Ala-D-Ala carboxypeptidases) and M38 (isoaspartyl dipeptidases). Based on genomic information, a wide range of abilities of melanin biosynthesis was revealed; genes related to metabolically distinct DHN, DOPA and pyomelanin pathways were identified. The MAT (MAting Type) locus and other sex-related genes were recognized in all 23 black fungi. Members of the asexual genera Fonsecaea and Cladophialophora appear to be heterothallic with a single copy of either MAT-1-1 or MAT-1-2 in each individual. All Capronia species are homothallic as both MAT1-1 and MAT1-2 genes were found in each single genome. The genomic synteny of the MAT-locus flanking genes (SLA2-APN2-COX13) is not conserved in black fungi as is commonly observed in Eurotiomycetes, indicating a unique genomic context for MAT in those species. The heterokaryon (het) genes expansion associated with the low selective pressure at the MAT-locus suggests that a parasexual cycle may play an important role in generating diversity among those fungi.
Three closely related thermally dimorphic pathogens are causal agents of major fungal diseases affecting humans in the Americas: blastomycosis, histoplasmosis and paracoccidioidomycosis. Here we report the genome sequence and analysis of four strains of the etiological agent of blastomycosis, Blastomyces, and two species of the related genus Emmonsia, typically pathogens of small mammals. Compared to related species, Blastomyces genomes are highly expanded, with long, often sharply demarcated tracts of low GC-content sequence. These GC-poor isochore-like regions are enriched for gypsy elements, are variable in total size between isolates, and are least expanded in the avirulent B. dermatitidis strain ER-3 as compared with the virulent B. gilchristii strain SLH14081. The lack of similar regions in related species suggests these isochore-like regions originated recently in the ancestor of the Blastomyces lineage. While gene content is highly conserved between Blastomyces and related fungi, we identified changes in copy number of genes potentially involved in host interaction, including proteases and characterized antigens. In addition, we studied gene expression changes of B. dermatitidis during the interaction of the infectious yeast form with macrophages and in a mouse model. Both experiments highlight a strong antioxidant defense response in Blastomyces, and upregulation of dioxygenases in vivo suggests that dioxide produced by antioxidants may be further utilized for amino acid metabolism. We identify a number of functional categories upregulated exclusively in vivo, such as secreted proteins, zinc acquisition proteins, and cysteine and tryptophan metabolism, which may include critical virulence factors missed before in in vitro studies. Across the dimorphic fungi, loss of certain zinc acquisition genes and differences in amino acid metabolism suggest unique adaptations of Blastomyces to its host environment. These results reveal the dynamics of genome evolution and of factors contributing to virulence in Blastomyces.
The fission yeast clade--comprising Schizosaccharomyces pombe, S. octosporus, S. cryophilus, and S. japonicus--occupies the basal branch of Ascomycete fungi and is an important model of eukaryote biology. A comparative annotation of these genomes identified a near extinction of transposons and the associated innovation of transposon-free centromeres. Expression analysis established that meiotic genes are subject to antisense transcription during vegetative growth, which suggests a mechanism for their tight regulation. In addition, trans-acting regulators control new genes within the context of expanded functional modules for meiosis and stress response. Differences in gene content and regulation also explain why, unlike the budding yeast of Saccharomycotina, fission yeasts cannot use ethanol as a primary carbon source. These analyses elucidate the genome structure and gene regulation of fission yeast and provide tools for investigation across the Schizosaccharomyces clade.
We have sequenced the genomes of 18 isolates of the closely related human pathogenic fungi Coccidioides immitis and Coccidioides posadasii to more clearly elucidate population genomic structure, bringing the total number of sequenced genomes for each species to 10. Our data confirm earlier microsatellite-based findings that these species are genetically differentiated, but our population genomics approach reveals that hybridization and genetic introgression have recently occurred between the two species. The directionality of introgression is primarily from C. posadasii to C. immitis, and we find more than 800 genes exhibiting strong evidence of introgression in one or more sequenced isolates. We performed PCR-based sequencing of one region exhibiting introgression in 40 C. immitis isolates to confirm and better define the extent of gene flow between the species. We find more coding sequence than expected by chance in the introgressed regions, suggesting that natural selection may play a role in the observed genetic exchange. We find notable heterogeneity in repetitive sequence composition among the sequenced genomes and present the first detailed genome-wide profile of a repeat-induced point mutation (RIP) process distinctly different from what has been observed in Neurospora. We identify promiscuous HLA-I and HLA-II epitopes in both proteomes and discuss the possible implications of introgression and population genomic data for public health and vaccine candidate prioritization. This study highlights the importance of population genomic data for detecting subtle but potentially important phenomena such as introgression.
Comparative analysis of multiple genomes in a phylogenetic framework dramatically improves the precision and sensitivity of evolutionary inference, producing more robust results than single-genome analyses can provide. The genomes of 12 Drosophila species, ten of which are presented here for the first time (sechellia, simulans, yakuba, erecta, ananassae, persimilis, willistoni, mojavensis, virilis and grimshawi), illustrate how rates and patterns of sequence divergence across taxa can illuminate evolutionary processes on a genomic scale. These genome sequences augment the formidable genetic tools that have made Drosophila melanogaster a pre-eminent model for animal genetics, and will further catalyse fundamental research on mechanisms of development, cell biology, genetics, disease, neurobiology, behaviour, physiology and evolution. Despite remarkable similarities among these Drosophila species, we identified many putatively non-neutral changes in protein-coding genes, non-coding RNA genes, and cis-regulatory regions. These may prove to underlie differences in the ecology and behaviour of these diverse species.
Here we report a high-quality draft genome sequence of the domestic dog (Canis familiaris), together with a dense map of single nucleotide polymorphisms (SNPs) across breeds. The dog is of particular interest because it provides important evolutionary information and because existing breeds show great phenotypic diversity for morphological, physiological and behavioural traits. We use sequence comparison with the primate and rodent lineages to shed light on the structure and evolution of genomes and genes. Notably, the majority of the most highly conserved non-coding sequences in mammalian genomes are clustered near a small subset of genes with important roles in development. Analysis of SNPs reveals long-range haplotypes across the entire dog genome, and defines the nature of genetic diversity within and across breeds. The current SNP map now makes it possible for genome-wide association studies to identify genes responsible for diseases and traits, with important consequences for human and companion animal health.
