ABSTRACTParasitic chytrid fungi have emerged as a significant threat to amphibian species worldwide, necessitating the development of techniques to isolate these pathogens into sterile culture for research purposes. However, early methods of isolating chytrids from their hosts relied on killing amphibians. We modified a pre-existing protocol for isolating chytrids from infected animals to use toe clips and biopsies from toe webbing rather than euthanizing hosts, and distributed the protocol to interested researchers worldwide as part of the BiodivERsA project RACE – here called the RML protocol. In tandem, we developed a lethal procedure for isolating chytrids from tadpole mouthparts. Reviewing a database of use a decade after their inception, we find that these methods have been widely applied across at least 5 continents, 23 countries and in 62 amphibian species, and have been successfully used to isolate chytrids in remote field locations. Isolation of chytrids by the non-lethal RML protocol occured in 18% of attempts with 207 fungal isolates and three species of chytrid being recovered. Isolation of chytrids from tadpoles occured in 43% of attempts with 334 fungal isolates of one species (Batrachochytrium dendrobatidis) being recovered. Together, these methods have resulted in a significant reduction and refinement of our use of threatened amphibian species and have improved our ability to work with this important group of emerging fungal pathogens.

AbstractCandida auris is an emergent fungal pathogen of rising public health concern due to increasing reports of outbreaks in healthcare settings and resistance to multiple classes of antifungal drugs. While distantly related to the more common pathogens C. albicans and C. glabrata, C. auris is closely related to three rarely observed and often multidrug-resistant species, C. haemulonii, C. duobushaemulonii and C. pseudohaemulonii. Here, we generated and analyzed near complete genome assemblies and RNA-Seq-guided gene predictions for isolates from each of the four major C. auris clades and for C. haemulonii, C. duobushaemulonii and C. pseudohaemulonii. Our analyses mapped seven chromosomes and revealed chromosomal rearrangements between C. auris clades and related species. We found conservation of genes involved in mating and meiosis and identified both MTLa and MTLα C. auris isolates, suggesting the potential for mating between clades. Gene conservation analysis highlighted that many genes linked to drug resistance and virulence in other pathogenic Candida species are conserved in C. auris and related species including expanded families of transporters and lipases, as well as mutations and copy number variants in ERG11 that confer drug resistance. In addition, we found genetic features of the emerging species that likely underlie differences in virulence and drug response between these and other Candida species, including genes involved in cell wall structure. To begin to characterize the species-specific genes important for antifungal response, we profiled the gene expression of C. auris in response to voriconazole and amphotericin B and found induction of several transporters and metabolic regulators that may play a role in drug resistance. This study provides a comprehensive view of the genomic basis of drug resistance, potential for mating, and virulence in this emerging fungal clade.

AbstractCryptococcus gattii is a pathogenic yeast of humans and other animals, which causes disease predominantly in immunocompetent hosts. Infection begins when aerosolized yeast or spores enter the body, triggering an immune response, including engulfment by macrophages. To understand the early transcriptional signals in both the yeast and its mammalian host, we performed a time-course dual RNA-seq experiment for four lineages of C. gattii (VGI-IV) interacting with mouse macrophages at 1hr, 3hr and 6hr post infection. Comparison of in vitro to ex vivo gene expression indicates lineage VGII is transcriptionally divergent to non-VGII lineages, including differential expression of genes involved in capsule synthesis, capsule attachment and ergosterol production. Various paralogs demonstrate sub-functionalisation between lineages including an upregulation of capsule biosynthesis-related gene CAP2, and downregulation of CAP1 in VGIII. Isolates also compensate for lineage-specific gene-losses by over-expression of genetically similar paralogs, including an over-expression of capsule gene CAS3 in VGIV having lost CAS31. Differential expression of one in five C. gattii genes was detected following co-incubation with mouse macrophages; all isolates showed high induction of oxidative-reduction functions and a downregulation of capsule attachment genes. We also show that VGII switches expression of two laccase paralogs (from LAC1 to LAC2) during co-incubation of macrophages. Finally, we found that mouse macrophages respond to all four lineages of C. gattii by upregulating FosB/Jun/Egr1 regulatory proteins at early time points. This study highlights the evolutionary breadth of expression profiles amongst the lineages of C. gattii and the diversity of transcriptional responses at this host-pathogen interface.ImportanceThe transcriptional profiles of related pathogens and their response to host induced stresses underpin their pathogenicity. Expression differences between related pathogens during host interaction can indicate when and how these genes contribute to virulence, ultimately informing new and improved treatment strategies for those diseases. In this paper, we compare the transcriptional profiles of five isolates representing four lineages of C. gattii in rich media. Our analyses identified key processes including cell capsule, ergosterol production and melanin that are differentially expressed between lineages, and we find that VGII has the most distinct profile in terms of numbers of differentially expressed genes. All lineages have also undergone sub-functionalisation for various paralogs including capsule biosynthesis and attachment genes. Most genes appeared down-regulated during co-incubation with macrophages, with the largest decrease observed for capsule attachment genes, which appears coordinated with a stress response, as all lineages also upregulated oxidative stress response genes. Furthermore, VGII upregulated many genes that are linked to ergosterol biosynthesis and switched expression of the laccase LAC1 to LAC2 ex vivo. Finally, we saw a pronounced increase in the FosB/Jun/Egr1 regulatory proteins at early time points in bone marrow derived macrophages, marking a role in the host response to C. gattii. This work highlights the dynamic roles of key C. gattii virulence genes in response to macrophages.

ABSTRACT
We discovered a new lineage of the globally important fungal pathogen Cryptococcus gattii on the basis of analysis of six isolates collected from three locations spanning the Central Miombo Woodlands of Zambia, Africa. All isolates were from environments (middens and tree holes) that are associated with a small mammal, the African hyrax. Phylogenetic and population genetic analyses confirmed that these isolates form a distinct, deeply divergent lineage, which we name VGV. VGV comprises two subclades (A and B) that are capable of causing mild lung infection with negligible neurotropism in mice. Comparing the VGV genome to previously identified lineages of C. gattii revealed a unique suite of genes together with gene loss and inversion events. However, standard URA5 restriction fragment length polymorphism (RFLP) analysis could not distinguish between VGV and VGIV isolates. We therefore developed a new URA5 RFLP method that can reliably identify the newly described lineage. Our work highlights how sampling understudied ecological regions alongside genomic and functional characterization can broaden our understanding of the evolution and ecology of major global pathogens.
IMPORTANCE Cryptococcus gattii is an environmental pathogen that causes severe systemic infection in immunocompetent individuals more often than in immunocompromised humans. Over the past 2 decades, researchers have shown that C. gattii falls within four genetically distinct major lineages. By combining field work from an understudied ecological region (the Central Miombo Woodlands of Zambia, Africa), genome sequencing and assemblies, phylogenetic and population genetic analyses, and phenotypic characterization (morphology, histopathological, drug-sensitivity, survival experiments), we discovered a hitherto unknown lineage, which we name VGV (variety gattii five). The discovery of a new lineage from an understudied ecological region has far-reaching implications for the study and understanding of fungal pathogens and diseases they cause.

Candida auris is an emergent multidrug-resistant fungal pathogen causing increasing reports of outbreaks. While distantly related to C. albicans and C. glabrata, C. auris is closely related to rarely observed and often multidrug-resistant species from the C. haemulonii clade. Here, we analyze near complete genome assemblies for the four C. auris clades and three related species, and map intra- and inter-species rearrangements across the seven chromosomes. Using RNA-Seq-guided gene predictions, we find that most mating and meiosis genes are conserved and that clades contain either the MTLa or MTLα mating loci. Comparing the genomes of these emerging species to those of other Candida species identifies genes linked to drug resistance and virulence, including expanded families of transporters and lipases, as well as mutations and copy number variants in ERG11. Gene expression analysis identifies transporters and metabolic regulators specific to C. auris and those conserved with related species which may contribute to differences in drug response in this emerging fungal clade.

Globalized infectious diseases are causing species declines worldwide, but their source often remains elusive. We used whole-genome sequencing to solve the spatiotemporal origins of the most devastating panzootic to date, caused by the fungus Batrachochytrium dendrobatidis, a proximate driver of global amphibian declines. We traced the source of B. dendrobatidis to the Korean peninsula, where one lineage, BdASIA-1, exhibits the genetic hallmarks of an ancestral population that seeded the panzootic. We date the emergence of this pathogen to the early 20th century, coinciding with the global expansion of commercial trade in amphibians, and we show that intercontinental transmission is ongoing. Our findings point to East Asia as a geographic hotspot for B. dendrobatidis biodiversity and the original source of these lineages that now parasitize amphibians worldwide.

Cryptococcus gattii is a pathogenic yeast of humans and other animals which causes disease predominantly in immunocompetent hosts. Infection begins when aerosolized yeast or spores enter the body, triggering an immune response, including engulfment by macrophages. To understand the early transcriptional signals in both the yeast and its mammalian host, we performed a time-course dual-transcriptome sequencing (RNA-seq) experiment for four lineages of C. gattii (lineages VGI to IV) interacting with mouse macrophages at 1, 3, and 6 h postinfection. Comparisons of in vitro to ex vivo gene expression levels indicated that lineage VGII is transcriptionally divergent from non-VGII lineages, including differential expression of genes involved in capsule synthesis, capsule attachment, and ergosterol production. Several paralogous genes demonstrated subfunctionalization between lineages, including upregulation of capsule biosynthesis-related gene CAP2 and downregulation of CAP1 in VGIII. Isolates also compensate for lineage-specific gene losses by overexpression of genetically similar paralogs, including overexpression of capsule gene CAS3 in VGIV, which have lost the CAS31 gene. Differential expression of one in five C. gattii genes was detected following coincubation with mouse macrophages; all isolates showed high induction of oxidative-reduction functions and downregulation of capsule attachment genes. We also found that VGII switches expression of two laccase paralogs (from LAC1 to LAC2) during coincubation of macrophages. Finally, we found that mouse macrophages respond to all four lineages of C. gattii by upregulating FosB/Jun/Egr1 regulatory proteins at early time points. This report highlights the evolutionary breadth of expression profiles among the lineages of C. gattii and the diversity of transcriptional responses at this host-pathogen interface.IMPORTANCE the transcriptional profiles of related pathogens and their responses to host-induced stresses underpin their pathogenicity. Expression differences between related pathogens during host interaction can indicate when and how these genes contribute to virulence, ultimately informing new and improved treatment strategies for those diseases. In this paper, we compare the transcriptional profiles of five isolates representing four lineages of C. gattii in rich media. Our analyses identified key processes, including those involving cell capsule, ergosterol production, and melanin, that are differentially expressed between lineages, and we found that VGII has the most distinct profile in terms of numbers of differentially expressed genes. All lineages have also undergone subfunctionalization for several paralogs, including capsule biosynthesis and attachment genes. Most genes appeared downregulated during coincubation with macrophages, with the largest decrease observed for capsule attachment genes, which appeared to be coordinated with a stress response, as all lineages also upregulated oxidative stress response genes. Furthermore, VGII upregulated many genes that are linked to ergosterol biosynthesis and switched from expression of the laccase LAC1 to expression of LAC2 ex vivo Finally, we saw a pronounced increase in the FosB/Jun/Egr1 regulatory proteins at early time points in bone marrow-derived macrophages, marking a role in the host response to C. gattii This work highlights the dynamic roles of key C. gattii virulence genes in response to macrophages.

To understand the evolutionary pathways that lead to emerging infections of vertebrates, here we explore the genomic innovations that allow free-living chytrid fungi to adapt to and colonize amphibian hosts. Sequencing and comparing the genomes of two pathogenic species of Batrachochytrium to those of close saprophytic relatives reveals that pathogenicity is associated with remarkable expansions of protease and cell wall gene families, while divergent infection strategies are linked to radiations of lineage-specific gene families. By comparing the host-pathogen response to infection for both pathogens, we illuminate the traits that underpin a strikingly different immune response within a shared host species. Our results show that, despite commonalities that promote infection, specific gene-family radiations contribute to distinct infection strategies. The breadth and evolutionary novelty of candidate virulence factors that we discover underscores the urgent need to halt the advance of pathogenic chytrids and prevent incipient loss of biodiversity.

BACKGROUND: Ortholog prediction and synteny visualization across whole genomes are valuable methods for detecting and representing a range of evolutionary processes such as genome expansion, chromosomal rearrangement, and chromosomal translocation. Few standalone methods are currently available to visualize synteny across any number of annotated genomes. RESULTS: Here, I present a Synteny Imaging tool (Synima) written in Perl, which uses the graphical features of R. Synima takes orthologues computed from reciprocal best BLAST hits or OrthoMCL, and DAGchainer, and outputs an overview of genome-wide synteny in PDF. Each of these programs are included with the Synima package, and a pipeline for their use. Synima has a range of graphical parameters including size, colours, order, and labels, which are specified in a config file generated by the first run of Synima - and can be subsequently edited. Synima runs quickly on a command line to generate informative and publication quality figures. Synima is open source and freely available from https://github.com/rhysf/Synima under the MIT License. CONCLUSIONS: Synima should be a valuable tool for visualizing synteny between two or more annotated genome assemblies.

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.

Fever is a universal response to infection, and opportunistic pathogens such as Candida albicans have evolved complex circuitry to sense and respond to heat. Here we harness RNA-seq and ChIP-seq to discover that the heat shock transcription factor, Hsf1, binds distinct motifs in nucleosome-depleted promoter regions to regulate heat shock genes and genes involved in virulence in C. albicans. Consequently, heat shock increases C. albicans host cell adhesion, damage and virulence. Hsf1 activation depends upon the molecular chaperone Hsp90 under basal and heat shock conditions, but the effects are opposite and in part controlled at the level of Hsf1 expression and DNA binding. Finally, we demonstrate that Hsp90 regulates global transcription programs by modulating nucleosome levels at promoters of stress-responsive genes. Thus, we describe a mechanism by which C. albicans responds to temperature via Hsf1 and Hsp90 to orchestrate gene expression and chromatin architecture, thereby enabling thermal adaptation and virulence.

Emerging fungal pathogens cause an expanding burden of disease across the animal kingdom, including a rise in morbidity and mortality in humans. Yet, we currently have only a limited repertoire of available therapeutic interventions. A greater understanding of the mechanisms of fungal virulence and of the emergence of hypervirulence within species is therefore needed for new treatments and mitigation efforts. For example, over the past decade, an unusual lineage of Cryptococcus gattii, which was first detected on Vancouver Island, has spread to the Canadian mainland and the Pacific Northwest infecting otherwise healthy individuals. The molecular changes that led to the development of this hypervirulent cryptococcal lineage remain unclear. To explore this, we traced the history of similar microevolutionary events that can lead to changes in host range and pathogenicity. Here, we detail fine-resolution mapping of genetic differences between two highly related Cryptococcus gattii VGIIc isolates that differ in their virulence traits (phagocytosis, vomocytosis, macrophage death, mitochondrial tubularization and intracellular proliferation). We identified a small number of single site variants within coding regions that potentially contribute to variations in virulence. We then extended our methods across multiple lineages of C. gattii to study how selection is acting on key virulence genes within different lineages.This article is part of the themed issue 'Tackling emerging fungal threats to animal health, food security and ecosystem resilience'.

UNLABELLED: Cryptococcus gattii is a fungal pathogen of humans, causing pulmonary infections in otherwise healthy hosts. To characterize genomic variation among the four major lineages of C. gattii (VGI, -II, -III, and -IV), we generated, annotated, and compared 16 de novo genome assemblies, including the first for the rarely isolated lineages VGIII and VGIV. By identifying syntenic regions across assemblies, we found 15 structural rearrangements, which were almost exclusive to the VGI-III-IV lineages. Using synteny to inform orthology prediction, we identified a core set of 87% of C. gattii genes present as single copies in all four lineages. Remarkably, 737 genes are variably inherited across lineages and are overrepresented for response to oxidative stress, mitochondrial import, and metal binding and transport. Specifically, VGI has an expanded set of iron-binding genes thought to be important to the virulence of Cryptococcus, while VGII has expansions in the stress-related heat shock proteins relative to the other lineages. We also characterized genes uniquely absent in each lineage, including a copper transporter absent from VGIV, which influences Cryptococcus survival during pulmonary infection and the onset of meningoencephalitis. Through inclusion of population-level data for an additional 37 isolates, we identified a new transcontinental clonal group that we name VGIIx, mitochondrial recombination between VGII and VGIII, and positive selection of multidrug transporters and the iron-sulfur protein aconitase along multiple branches of the phylogenetic tree. Our results suggest that gene expansion or contraction and positive selection have introduced substantial variation with links to mechanisms of pathogenicity across this species complex. IMPORTANCE: the genetic differences between phenotypically different pathogens provide clues to the underlying mechanisms of those traits and can lead to new drug targets and improved treatments for those diseases. In this paper, we compare 16 genomes belonging to four highly differentiated lineages of Cryptococcus gattii, which cause pulmonary infections in otherwise healthy humans and other animals. Half of these lineages have not had their genomes previously assembled and annotated. We identified 15 ancestral rearrangements in the genome and over 700 genes that are unique to one or more lineages, many of which are associated with virulence. In addition, we found evidence for recent transcontinental spread, mitochondrial genetic exchange, and positive selection in multidrug transporters. Our results suggest that gene expansion/contraction and positive selection are diversifying the mechanisms of pathogenicity across this species complex.

Sequence alignments form the basis for many comparative and population genomic studies. Alignment tools provide a range of accuracies dependent on the divergence between the sequences and the alignment methods. Despite widespread use, there is no standard method for assessing the accuracy of a dataset and alignment strategy after resequencing. We present a framework and tool for determining the overall accuracies of an input read dataset, alignment and SNP-calling method providing an isolate in that dataset has a corresponding, or closely related reference sequence available. In addition to this tool for comparing False Discovery Rates (FDR), we include a method for determining homozygous and heterozygous positions from an alignment using binomial probabilities for an expected error rate. We benchmark this method against other SNP callers using our FDR method with three fungal genomes, finding that it was able achieve a high level of accuracy. These tools are available at http://cfdr.sourceforge.net/.

Chytridiomycosis is a cutaneous infection of amphibians caused by the chytridiomycete fungal pathogen Batrachochytrium dendrobatidis (Bd). Despite being in a phylum not known for pathogenicity in vertebrates, Bd is now recognized as a primary driver of amphibian declines. Data show that this novel pathogen emerged in the 20th century to colonize amphibians worldwide. Such rapid emergence of a previously unrecognized pathogen illustrates many aspects of emerging fungal infections that threaten human health, namely long-distance human-mediated dispersal, multihost reservoirs, and altered virulence. In order to combat Bd, new tools have been developed to track its global spread and to analyze in parallel whole-genome diversity. This article details how such tools have applications to tracking and managing human fungal infections.

Many plant pathogens, including those in the lineage of the Irish potato famine organism Phytophthora infestans, evolve by host jumps followed by specialization. However, how host jumps affect genome evolution remains largely unknown. To determine the patterns of sequence variation in the P. infestans lineage, we resequenced six genomes of four sister species. This revealed uneven evolutionary rates across genomes with genes in repeat-rich regions showing higher rates of structural polymorphisms and positive selection. These loci are enriched in genes induced in planta, implicating host adaptation in genome evolution. Unexpectedly, genes involved in epigenetic processes formed another class of rapidly evolving residents of the gene-sparse regions. These results demonstrate that dynamic repeat-rich genome compartments underpin accelerated gene evolution following host jumps in this pathogen lineage.

Illumina's Genome Analyzer generates ultra-short sequence reads, typically 36 nucleotides in length, and is primarily intended for resequencing. We tested the potential of this technology for de novo sequence assembly on the 6 Mbp genome of Pseudomonas syringae pv. syringae B728a with several freely available assembly software packages. Using an unpaired data set, velvet assembled >96% of the genome into contigs with an N50 length of 8289 nucleotides and an error rate of 0.33%. EDENA generated smaller contigs (N50 was 4192 nucleotides) and comparable error rates. SSAKE and VCAKE yielded shorter contigs with very high error rates. Assembly of paired-end sequence data carrying 400 bp inserts produced longer contigs (N50 up to 15 628 nucleotides), but with increased error rates (0.5%). Contig length and error rate were very sensitive to the choice of parameter values. Noncoding RNA genes were poorly resolved in de novo assemblies, while >90% of the protein-coding genes were assembled with 100% accuracy over their full length. This study demonstrates that, in practice, de novo assembly of 36-nucleotide reads can generate reasonably accurate assemblies from about 40 x deep sequence data sets. These draft assemblies are useful for exploring an organism's proteomic potential, at a very economic low cost.