Clement K. M. Tsui

Genome and transcriptome analyses of the mountain pine beetle-fungal symbiont Grosmannia clavigera, a lodgepole pine pathogen

In western North America, the current outbreak of the mountain pine beetle (MPB) and its microbial associates has destroyed wide areas of lodgepole pine forest, including more than 16 million hectares in British Columbia. Grosmannia clavigera (Gc), a critical component of the outbreak, is a symbiont of the MPB and a pathogen of pine trees. To better understand the interactions between Gc, MPB, and lodgepole pine hosts, we sequenced the ∼30-Mb Gc genome and assembled it into 18 supercontigs. We predict 8,314 protein-coding genes, and support the gene models with proteome, expressed sequence tag, and RNA-seq data. We establish that Gc is heterothallic, and report evidence for repeat-induced point mutation. We report insights, from genome and transcriptome analyses, into how Gc tolerates conifer-defense chemicals, including oleoresin terpenoids, as they colonize a host tree. RNA-seq data indicate that terpenoids induce a substantial antimicrobial stress in Gc, and suggest that the fungus may detoxify these chemicals by using them as a carbon source. Terpenoid treatment strongly activated a ∼100-kb region of the Gc genome that contains a set of genes that may be important for detoxification of these host-defense chemicals. This work is a major step toward understanding the biological interactions between the tripartite MPB/fungus/forest system.

The Faces of Fungi database: fungal names linked with morphology, phylogeny and human impacts

Taxonomic names are key links between various databases that store information on different organisms. Several global fungal nomenclural and taxonomic databases (notably Index Fungorum, Species Fungorum and MycoBank) can be sourced to find taxonomic details about fungi, while DNA sequence data can be sourced from NCBI, EBI and UNITE databases. Although the sequence data may be linked to a name, the quality of the metadata is variable and generally there is no corresponding link to images, descriptions or herbarium material. There is generally no way to establish the accuracy of the names in these genomic databases, other than whether the submission is from a reputable source. To tackle this problem, a new database (FacesofFungi), accessible at www.facesoffungi.org (FoF) has been established. This fungal database allows deposition of taxonomic data, phenotypic details and other useful data, which will enhance our current taxonomic understanding and ultimately enable mycologists to gain better and updated insights into the current fungal classification system. In addition, the database will also allow access to comprehensive metadata including descriptions of voucher and type specimens. This database is user-friendly, providing links and easy access between taxonomic ranks, with the classification system based primarily on molecular data (from the literature and via updated web-based phylogenetic trees), and to a lesser extent on morphological data when molecular data are unavailable. In FoF species are not only linked to the closest phylogenetic representatives, but also relevant data is provided, wherever available, on various applied aspects, such as ecological, industrial, quarantine and chemical uses. The data include the three main fungal groups (Ascomycota, Basidiomycota, Basal fungi) and fungus-like organisms. The FoF webpage is an output funded by the Mushroom Research Foundation which is an NGO with seven directors with mycological expertise. The webpage has 76 curators, and with the help of these specialists, FoF will provide an updated natural classification of the fungi, with illustrated accounts of species linked to molecular data. The present paper introduces the FoF database to the scientific community and briefly reviews some of the problems associated with classification and identification of the main fungal groups. The structure and use of the database is then explained. We would like to invite all mycologists to contribute to these web pages.

Revision of lignicolous Tubeufiaceae based on morphological reexamination and phylogenetic analysis

In this paper we revisit the family Tubeufiaceae with notes on genera that we have re-examined where possible. Generic type specimens of Acanthophiobolus, Kamalomyces, Podonectria, Thaxteriella and Thaxteriellopsis were re-examined, described and illustrated and shown to belong to Tubeufiaceae. Notes are provided on Acanthostigma, Chaetosphaerulina, Thaxterina and Tubeufia, which are retained in Tubeufiaceae; however, we were unable to locate the types of these genera during the time frame of this study. Allonecte is excluded from the Tubeufiaceae, as the ascospores are fusiform-ellipsoidal, grey-brown and 1-septate and the asci are cylindrical, all of which are features more typical of Pleosporaceae, where it is transferred. Byssocallis has yellow to orange ascomata and clavate ascospores which is atypical of Tubeufiaceae. Thus its taxonomic status needs to be reevaluated. Lentendraeopsis has an endophytic habit, cylindro-clavate asci and two-celled ascospores more typical of Pleosporales, where it is transferred. Taphrophila has small ascomata, a thin peridium, branching setae around the apex of the ascomata, clavate to saccate asci and lacks pseudoparaphyses. These are features atypical of the Tubeufiaceae, and Taphrophila should be placed in the Dothideomycetes incertae cedis. Twelve new collections of Tubeufiaceae from Thailand were isolated, and their DNA was extracted. The sequence data of LSU, SSU and ITS rDNA were amplified and analyzed using parsimony and likelihood methods. The results of phylogenetic analysis was used to establish the inter-generic relationships in Tubeufiaceae. Thaxteriellopsis lignicola, epitypified in this investigation, is a sister taxon in the family Tubeufiaceae based on phylogenetic analysis of rRNA sequence data. Chlamydotubeufia is introduced as a new genus based on the production of dictyochlamydosporous anamorphs, including two new species. Three new species, one each in Acanthostigma, Tubeufia and Thaxteriella are also described and illustrated. The phylogenetic placement of these genera is also discussed.

Population structure and migration pattern of a conifer pathogen, Grosmannia clavigera, as influenced by its symbiont, the mountain pine beetle

We investigated the population structure of Grosmannia clavigera (Gc), a fungal symbiont of the mountain pine beetle (MPB) that plays a crucial role in the establishment and reproductive success of this pathogen. This insect–fungal complex has destroyed over 16 million ha of lodgepole pine forests in Canada, the largest MPB epidemic in recorded history. During this current epidemic, MPB has expanded its range beyond historically recorded boundaries, both northward and eastward, and has now reached the jack pine of Alberta, potentially threatening the Canadian boreal forest. To better understand the dynamics between the beetle and its fungal symbiont, we sampled 19 populations in western North America and genotyped individuals from these populations with eight microsatellite markers. The fungus displayed high haplotype diversity, with over 250 unique haplotypes observed in 335 single spore isolates. Linkage equilibria in 13 of the 19 populations suggested that the fungus reproduces sexually. Bayesian clustering and distance analyses identified four genetic clusters that corresponded to four major geographical regions, which suggested that the epidemic arose from multiple geographical sources. A genetic cluster north of the Rocky Mountains, where the MPB has recently become established, experienced a population bottleneck, probably as a result of the recent range expansion. The two genetic clusters located north and west of the Rocky Mountains contained many fungal isolates admixed from all populations, possibly due to the massive movement of MPB during the epidemic. The general agreement in north–south differentiation of MPB and G. clavigera populations points to the fungal pathogen’s dependence on the movement of its insect vector. In addition, the patterns of diversity and the individual assignment tests of the fungal associate suggest that migration across the Rocky Mountains occurred via a northeastern corridor, in accordance with meteorological patterns and observation of MPB movement data. Our results highlight the potential of this pathogen for both expansion and sexual reproduction, and also identify some possible barriers to gene flow. Understanding the ecological and evolutionary dynamics of this fungus–beetle association is important for the modelling and prediction of MPB epidemics.

Unequal Recombination and Evolution of the Mating-Type (MAT) Loci in the Pathogenic Fungus Grosmannia clavigera and Relatives

Sexual reproduction in fungi is regulated by the mating-type (MAT) locus where recombination is suppressed. We investigated the evolution of MAT loci in eight fungal species belonging to Grosmannia and Ophiostoma (Sordariomycetes, Ascomycota) that include conifer pathogens and beetle symbionts. The MAT1-2 idiomorph/allele was identified from the assembled and annotated Grosmannia clavigera genome, and the MAT locus is flanked by genes coding for cytoskeleton protein (SLA) and DNA lyase. The synteny of these genes is conserved and consistent with other members in Ascomycota. Using sequences from SLA and flanking regions, we characterized the MAT1-1 idiomorph from other isolates of G. clavigera and performed dotplot analysis between the two idiomorphs. Unexpectedly, the MAT1-2 idiomorph contains a truncated MAT1-1-1 gene upstream of the MAT1-2-1 gene that bears the high-mobility-group domain. The nucleotide and amino acid sequence of the truncated MAT1-1-1 gene is similar to its homologous copy in the MAT1-1 idiomorph in the opposite mating-type isolate, except that positive selection is acting on the truncated gene and the alpha(α)-box that encodes the transcription factor has been deleted. The MAT idiomorphs sharing identical gene organization were present in seven additional species in the Ophiostomatales, suggesting that the presence of truncated MAT1-1-1 gene is a general pattern in this order. We propose that an ancient unequal recombination event resulted in the ancestral MAT1-1-1 gene integrated into the MAT1-2 idiomorph and surviving as the truncated MAT1-1-1 genes. The α-box domain of MAT1-1-1 gene, located at the same MAT locus adjacent to the MAT1-2-1 gene, could have been removed by deletion after recombination due to mating signal interference. Our data confirmed a 1:1 MAT/sex ratio in two pathogen populations, and showed that all members of the Ophiostomatales studied here including those that were previously deemed asexual have the potential to reproduce sexually. This ability can potentially increase genetic variability and can enhance fitness in new, ecological niches.

Rapid identification and detection of pine pathogenic fungi associated with mountain pine beetles by padlock probes.

Fifteen million hectares of pine forests in western Canada have been attacked by the mountain pine beetle (Dendroctonus ponderosae; MPB), leading to devastating economic losses. Grosmannia clavigera and Leptographium longiclavatum, are two fungi intimately associated with the beetles, and are crucial components of the epidemic. To detect and discriminate these two closely related pathogens, we utilized a method based on ligase-mediated nucleotide discrimination with padlock probe technology, and signal amplification by hyperbranched rolling circle amplification (HRCA). Two padlock probes were designed to target species-specific single nucleotide polymorphisms (SNPs) located at the inter-generic spacer 2 region and large subunit of the rRNA respectively, which allows discrimination between the two species. Thirty-four strains of G. clavigera and twenty-five strains of L. longiclavatum representing a broad geographic origin were tested with this assay. The HRCA results were largely in agreement with the conventional identification based on morphology or DNA-based methods. Both probes can also efficiently distinguish the two MPB-associated fungi from other fungi in the MPB, as well as other related fungi in the order Ophiostomatales. We also tested this diagnostic method for the direct detection of these fungi from the DNA of MPB. A nested PCR approach was used to enrich amplicons for signal detection. The results confirmed the presence of these two fungi in MPB. Thus, the padlock probe assay coupled with HRCA is a rapid, sensitive and reproducible method for the identification and detection of these ophiostomatoid fungi.

Asexual propagation of a virulent clone complex in a human and feline outbreak of sporotrichosis.

ABSTRACT Sporotrichosis is one of the most frequent subcutaneous fungal infections in humans and animals caused by members of the plant-associated, dimorphic genus Sporothrix. Three of the four medically important Sporothrix species found in Brazil have been considered asexual as no sexual stage has ever been reported in Sporothrix schenckii, Sporothrix brasiliensis, or Sporothrix globosa. We have identified the mating type (MAT) loci in the S. schenckii (strain 1099-18/ATCC MYA-4821) and S. brasiliensis (strain 5110/ATCC MYA-4823) genomes by using comparative genomic approaches to determine the mating type ratio in these pathogen populations. Our analysis revealed the presence of a MAT1-1 locus in S. schenckii while a MAT1-2 locus was found in S. brasiliensis representing genomic synteny to other Sordariomycetes. Furthermore, the components of the mitogen-activated protein kinase (MAPK)-pheromone pathway, pheromone processing enzymes, and meiotic regulators have also been identified in the two pathogens, suggesting the potential for sexual reproduction. The ratio of MAT1-1 to MAT1-2 was not significantly different from 1:1 for all three Sporothrix species, but the population of S. brasiliensis in the outbreaks originated from a single mating type. We also explored the population genetic structure of these pathogens using sequence data of two loci to improve our knowledge of the pattern of geographic distribution, genetic variation, and virulence phenotypes. Population genetics data showed significant population differentiation and clonality with a low level of haplotype diversity in S. brasiliensis isolates from different regions of sporotrichosis outbreaks in Brazil. In contrast, S. schenckii isolates demonstrated a high degree of genetic variability without significant geographic differentiation, indicating the presence of recombination. This study demonstrated that two species causing the same disease have contrasting reproductive strategies and genetic variability patterns.

Giardia spp. Are Commonly Found in Mixed Assemblages in Surface Water, as Revealed by Molecular and Whole-Genome Characterization.

ABSTRACT Giardia is the most common parasitic cause of gastrointestinal infections worldwide, with transmission through surface water playing an important role in various parts of the world. Giardia duodenalis (synonyms: G. intestinalis and G. lamblia), a multispecies complex, has two zoonotic subtypes, assemblages A and B. When British Columbia (BC), a western Canadian province, experienced several waterborne giardiasis outbreaks due to unfiltered surface drinking water in the late 1980s, collection of isolates from surface water, as well as from humans and beavers (Castor canadensis), throughout the province was carried out. To better understand Giardia in surface water, 71 isolates, including 29 from raw surface water samples, 29 from human giardiasis cases, and 13 from beavers in watersheds from this historical library were characterized by PCR. Study isolates also included isolates from waterborne giardiasis outbreaks. Both assemblages A and B were identified in surface water, human, and beavers samples, including a mixture of both assemblages A and B in waterborne outbreaks. PCR results were confirmed by whole-genome sequencing (WGS) for one waterborne outbreak and supported the clustering of human, water, and beaver isolates within both assemblages. We concluded that contamination of surface water by Giardia is complex, that the majority of our surface water isolates were assemblage B, and that both assemblages A and B may cause waterborne outbreaks. The higher-resolution data provided by WGS warrants further study to better understand the spread of Giardia.

Single-nucleotide polymorphism discovery in Leptographium longiclavatum, a mountain pine beetle-associated symbiotic fungus, using whole-genome resequencing

Single‐nucleotide polymorphisms (SNPs) are rapidly becoming the standard markers in population genomics studies; however, their use in nonmodel organisms is limited due to the lack of cost‐effective approaches to uncover genome‐wide variation, and the large number of individuals needed in the screening process to reduce ascertainment bias. To discover SNPs for population genomics studies in the fungal symbionts of the mountain pine beetle (MPB), we developed a road map to discover SNPs and to produce a genotyping platform. We undertook a whole‐genome sequencing approach of Leptographium longiclavatum in combination with available genomics resources of another MPB symbiont, Grosmannia clavigera. We sequenced 71 individuals pooled into four groups using the Illumina sequencing technology. We generated between 27 and 30 million reads of 75 bp that resulted in a total of 1, 181 contigs longer than 2 kb and an assembled genome size of 28.9 Mb (N50 = 48 kb, average depth = 125x). A total of 9052 proteins were annotated, and between 9531 and 17 266 SNPs were identified in the four pools. A subset of 206 genes (containing 574 SNPs, 11% false positives) was used to develop a genotyping platform for this species. Using this roadmap, we developed a genotyping assay with a total of 147 SNPs located in 121 genes using the Illumina® Sequenom iPLEX Gold. Our preliminary genotyping (success rate = 85%) of 304 individuals from 36 populations supports the utility of this approach for population genomics studies in other MPB fungal symbionts and other fungal nonmodel species.

Colonization history, host distribution, anthropogenic influence and landscape features shape populations of white pine blister rust, an invasive alien tree pathogen.

White pine blister rust is caused by the fungal pathogen Cronartium ribicola J.C. Fisch (Basidiomycota, Pucciniales). This invasive alien pathogen was introduced into North America at the beginning of the 20th century on pine seedlings imported from Europe and has caused serious economic and ecological impacts. In this study, we applied a population and landscape genetics approach to understand the patterns of introduction and colonization as well as population structure and migration of C. ribicola. We characterized 1,292 samples of C. ribicola from 66 geographic locations in North America using single nucleotide polymorphisms (SNPs) and evaluated the effect of landscape features, host distribution, and colonization history on the structure of these pathogen populations. We identified eastern and western genetic populations in North America that are strongly differentiated. Genetic diversity is two to five times higher in eastern populations than in western ones, which can be explained by the repeated accidental introductions of the pathogen into northeastern North America compared with a single documented introduction into western North America. These distinct genetic populations are maintained by a barrier to gene flow that corresponds to a region where host connectivity is interrupted. Furthermore, additional cryptic spatial differentiation was identified in western populations. This differentiation corresponds to landscape features, such as mountain ranges, and also to host connectivity. We also detected genetic differentiation between the pathogen populations in natural stands and plantations, an indication that anthropogenic movement of this pathogen still takes place. These results highlight the importance of monitoring this invasive alien tree pathogen to prevent admixture of eastern and western populations where different pathogen races occur.