Highly Recombinant VGII Cryptococcus gattii Population Develops Clonal Outbreak Clusters through both Sexual Macroevolution and Asexual Microevolution

ABSTRACT An outbreak of the fungal pathogen Cryptococcus gattii began in the Pacific Northwest (PNW) in the late 1990s. This outbreak consists of three clonal subpopulations: VGIIa/major, VGIIb/minor, and VGIIc/novel. Both VGIIa and VGIIc are unique to the PNW and exhibit increased virulence. In this study, we sequenced the genomes of isolates from these three groups, as well as global isolates, and analyzed a total of 53 isolates. We found that VGIIa/b/c populations show evidence of clonal expansion in the PNW. Whole-genome sequencing provided evidence that VGIIb originated in Australia, while VGIIa may have originated in South America, and these were likely independently introduced. Additionally, the VGIIa outbreak lineage may have arisen from a less virulent clade that contained a mutation in the MSH2 ortholog, but this appears to have reverted in the VGIIa outbreak strains, suggesting that a transient mutator phenotype may have contributed to adaptation and evolution of virulence in the PNW outbreak. PNW outbreak isolates share genomic islands, both between the clonal lineages and with global isolates, indicative of sexual recombination. This suggests that VGII C. gattii has undergone sexual reproduction, either bisexual or unisexual, in multiple locales contributing to the production of novel, virulent subtypes. We also found that the genomes of two basal VGII isolates from HIV+ patients contain an introgression tract spanning three genes. Introgression substantially contributed to intra-VGII polymorphism and likely occurred through sexual reproduction with VGI. More broadly, these findings illustrate how both microevolution and sexual reproduction play central roles in the development of infectious outbreaks from avirulent or less virulent progenitors. IMPORTANCE Cryptococcus gattii is the causative agent responsible for ongoing infections in the Pacific Northwest of the United States and western Canada. The incidence of these infections increased dramatically in the 1990s and remains elevated. These infections are attributable to three clonal lineages of C. gattii, VGIIa, VGIIb, and VGIIc, with only VGIIa identified once previously in the Pacific Northwest prior to the start of the outbreak, albeit in a less virulent form. This study addresses the origin and emergence of this outbreak, using whole-genome sequencing and comparison of both outbreak and global isolates. We show that VGIIa arose mitotically from a less virulent clonal group, possibly via the action of a mutator phenotype, while VGIIb was likely introduced from Australia, and VGIIc appears to have emerged in the United States or in an undersampled locale via sexual reproduction. This work shows that multiple processes can contribute to the emergence of an outbreak. Cryptococcus gattii is the causative agent responsible for ongoing infections in the Pacific Northwest of the United States and western Canada. The incidence of these infections increased dramatically in the 1990s and remains elevated. These infections are attributable to three clonal lineages of C. gattii, VGIIa, VGIIb, and VGIIc, with only VGIIa identified once previously in the Pacific Northwest prior to the start of the outbreak, albeit in a less virulent form. This study addresses the origin and emergence of this outbreak, using whole-genome sequencing and comparison of both outbreak and global isolates. We show that VGIIa arose mitotically from a less virulent clonal group, possibly via the action of a mutator phenotype, while VGIIb was likely introduced from Australia, and VGIIc appears to have emerged in the United States or in an undersampled locale via sexual reproduction. This work shows that multiple processes can contribute to the emergence of an outbreak.

[1]  W. Meyer,et al.  Antifungal Drug Susceptibility and Phylogenetic Diversity among Cryptococcus Isolates from Dogs and Cats in North America , 2014, Journal of Clinical Microbiology.

[2]  Jacqueline E. Schein,et al.  Analysis of the Genome and Transcriptome of Cryptococcus neoformans var. grubii Reveals Complex RNA Expression and Microevolution Leading to Virulence Attenuation , 2014, PLoS genetics.

[3]  T. Myers,et al.  Microevolution During Serial Mouse Passage Demonstrates FRE3 as a Virulence Adaptation Gene in Cryptococcus neoformans , 2014, mBio.

[4]  Alexandros Stamatakis,et al.  RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies , 2014, Bioinform..

[5]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[6]  D. Vugia,et al.  Cryptococcus gattii in the United States: Genotypic Diversity of Human and Veterinary Isolates , 2013, PloS one.

[7]  P. Hedrick Adaptive introgression in animals: examples and comparison to new mutation and standing variation as sources of adaptive variation , 2013, Molecular ecology.

[8]  J. Heitman,et al.  Unisexual and Heterosexual Meiotic Reproduction Generate Aneuploidy and Phenotypic Diversity De Novo in the Yeast Cryptococcus neoformans , 2013, PLoS biology.

[9]  Rhys A. Farrer,et al.  Transmission of Hypervirulence Traits via Sexual Reproduction within and between Lineages of the Human Fungal Pathogen Cryptococcus gattii , 2013, PLoS genetics.

[10]  Leen Stougie,et al.  Ancient Dispersal of the Human Fungal Pathogen Cryptococcus gattii from the Amazon Rainforest , 2013, PloS one.

[11]  A. Fauci,et al.  Emerging Infectious Diseases: Threats to Human Health and Global Stability , 2013, PLoS pathogens.

[12]  A. Idnurm,et al.  Congenic Strains for Genetic Analysis of Virulence Traits in Cryptococcus gattii , 2013, Infection and Immunity.

[13]  J. Stajich,et al.  Complex history of the amphibian-killing chytrid fungus revealed with genome resequencing data , 2013, Proceedings of the National Academy of Sciences.

[14]  Rajesh Kunadharaju,et al.  Cryptococcus gattii, Florida, USA, 2011 , 2013, Emerging infectious diseases.

[15]  Helga Thorvaldsdóttir,et al.  Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration , 2012, Briefings Bioinform..

[16]  A. Idnurm,et al.  DNA Mutations Mediate Microevolution between Host-Adapted Forms of the Pathogenic Fungus Cryptococcus neoformans , 2012, PLoS pathogens.

[17]  J. Heitman,et al.  Gene Conversion Occurs within the Mating-Type Locus of Cryptococcus neoformans during Sexual Reproduction , 2012, PLoS genetics.

[18]  Pablo Cingolani,et al.  © 2012 Landes Bioscience. Do not distribute. , 2022 .

[19]  Paul Keim,et al.  Whole Genome Sequence Analysis of Cryptococcus gattii from the Pacific Northwest Reveals Unexpected Diversity , 2011, PloS one.

[20]  M. Nei,et al.  MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. , 2011, Molecular biology and evolution.

[21]  Gonçalo R. Abecasis,et al.  The variant call format and VCFtools , 2011, Bioinform..

[22]  Kristin L. Tangen,et al.  Genome Variation in Cryptococcus gattii, an Emerging Pathogen of Immunocompetent Hosts , 2011, mBio.

[23]  D. Carter,et al.  Sexual Reproduction of Cryptococcus gattii: a Population Genetics Perspective , 2011 .

[24]  Melissa A. Miller,et al.  Self-Mating in the Definitive Host Potentiates Clonal Outbreaks of the Apicomplexan Parasites Sarcocystis neurona and Toxoplasma gondii , 2010, PLoS genetics.

[25]  M. DePristo,et al.  The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. , 2010, Genome research.

[26]  Joseph Heitman,et al.  Emergence and Pathogenicity of Highly Virulent Cryptococcus gattii Genotypes in the Northwest United States , 2010, PLoS pathogens.

[27]  E. Sionov,et al.  Cryptococcus neoformans Overcomes Stress of Azole Drugs by Formation of Disomy in Specific Multiple Chromosomes , 2010, PLoS pathogens.

[28]  J. Heitman,et al.  Spores as Infectious Propagules of Cryptococcus neoformans , 2009, Infection and Immunity.

[29]  Yi Guan,et al.  Dating the emergence of pandemic influenza viruses , 2009, Proceedings of the National Academy of Sciences.

[30]  Gavin J. D. Smith,et al.  Origins and evolutionary genomics of the 2009 swine-origin H1N1 influenza A epidemic , 2009, Nature.

[31]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[32]  Richard Durbin,et al.  Sequence analysis Fast and accurate short read alignment with Burrows – Wheeler transform , 2009 .

[33]  T. G. Mitchell,et al.  Molecular evidence that the range of the Vancouver Island outbreak of Cryptococcus gattii infection has expanded into the Pacific Northwest in the United States. , 2009, The Journal of infectious diseases.

[34]  K. Wannemuehler,et al.  Estimation of the current global burden of cryptococcal meningitis among persons living with HIV/AIDS , 2009, AIDS.

[35]  Hiten D. Madhani,et al.  Systematic Genetic Analysis of Virulence in the Human Fungal Pathogen Cryptococcus neoformans , 2008, Cell.

[36]  Erik L. L. Sonnhammer,et al.  Kalign – an accurate and fast multiple sequence alignment algorithm , 2005, BMC Bioinformatics.

[37]  J. Heitman,et al.  Same-sex mating and the origin of the Vancouver Island Cryptococcus gattii outbreak , 2005, Nature.

[38]  Fran Lewitter,et al.  Intragenic tandem repeats generate functional variability , 2005, Nature Genetics.

[39]  T. Boekhout,et al.  A rare genotype of Cryptococcus gattii caused the cryptococcosis outbreak on Vancouver Island (British Columbia, Canada). , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[40]  L. Hoang,et al.  Cryptococcus neoformans infections at Vancouver Hospital and Health Sciences Centre (1997-2002): epidemiology, microbiology and histopathology. , 2004, Journal of medical microbiology.

[41]  Lawrence Corey,et al.  Epidemiology and outcome of mould infections in hematopoietic stem cell transplant recipients. , 2002, Clinical infectious diseases : an official publication of the Infectious Diseases Society of America.

[42]  A. Hehl,et al.  Success and Virulence in Toxoplasma as the Result of Sexual Recombination Between Two Distinct Ancestries , 2001, Science.

[43]  K. Crandall,et al.  TCS: a computer program to estimate gene genealogies , 2000, Molecular ecology.

[44]  W. Fitch,et al.  Predicting the evolution of human influenza A. , 1999, Science.

[45]  Tomas A. Prolla,et al.  Destabilization of tracts of simple repetitive DNA in yeast by mutations affecting DNA mismatch repair , 1993, Nature.

[46]  R. Reenan,et al.  Isolation and characterization of two Saccharomyces cerevisiae genes encoding homologs of the bacterial HexA and MutS mismatch repair proteins. , 1992, Genetics.