Skip to main content
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Research and News from ASM Journals
    • mSphere of Influence: Commentaries from Early Career Microbiologists
    • Archive
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mSphere
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
mSphere
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Research and News from ASM Journals
    • mSphere of Influence: Commentaries from Early Career Microbiologists
    • Archive
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mSphere
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
Research Article | Ecological and Evolutionary Science

Dating the Cryptococcus gattii Dispersal to the North American Pacific Northwest

Chandler C. Roe, Jolene Bowers, Hanna Oltean, Emilio DeBess, Philippe J. Dufresne, Scott McBurney, David P. Overy, Bodo Wanke, Colleen Lysen, Tom Chiller, Wieland Meyer, George R. Thompson III, Shawn R. Lockhart, Crystal M. Hepp, David M. Engelthaler
J. Andrew Alspaugh, Editor
Chandler C. Roe
aTranslational Genomics Research Institute, Flagstaff, Arizona, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jolene Bowers
aTranslational Genomics Research Institute, Flagstaff, Arizona, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Hanna Oltean
bWashington State Department of Health, Shoreline, Washington, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Emilio DeBess
cPublic Health Division, Oregon Health Authority, Salem, Oregon, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Philippe J. Dufresne
dLaboratoire de Santé Publique du Québec, Institut National de Santé Publique du Québec, Charlottetown, Prince Edward Island, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Scott McBurney
eCanadian Wildlife Health Cooperative, Atlantic Veterinary College, University of Prince Edward Island, Charlottetown, Prince Edward Island, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
David P. Overy
eCanadian Wildlife Health Cooperative, Atlantic Veterinary College, University of Prince Edward Island, Charlottetown, Prince Edward Island, Canada
fAgriculture and Agrifood Canada, Ottawa Research and Development Centre, Ottawa, Ontario, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Bodo Wanke
gNational Institute of Infectious Diseases Evandro Chagas, FIOCRUZ, Rio de Janeiro, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Colleen Lysen
hMycotic Disease Branch, Centers for Disease Control and Prevention, Atlanta, Georgia, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Tom Chiller
hMycotic Disease Branch, Centers for Disease Control and Prevention, Atlanta, Georgia, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Wieland Meyer
iMolecular Mycology Research Laboratory, Center for Infectious Diseases and Microbiology, Marie Bashir Institute for Emerging Infections and Biosecurity, Westmead Clinical School, Sydney Medical School, The University of Sydney, Westmead, New South Wales, Australia
jWestmead Hospital, Westmead Institute for Medical Research, Westmead, New South Wales, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
George R. Thompson III
kUniversity of California, Davis, Davis, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Shawn R. Lockhart
hMycotic Disease Branch, Centers for Disease Control and Prevention, Atlanta, Georgia, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Crystal M. Hepp
lSchool of Informatics, Computing, and Cyber Systems, Northern Arizona University, Flagstaff, Arizona, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
David M. Engelthaler
aTranslational Genomics Research Institute, Flagstaff, Arizona, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for David M. Engelthaler
J. Andrew Alspaugh
Duke University Medical Center
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/mSphere.00499-17
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

The emergence of Cryptococcus gattii, previously regarded as a predominantly tropical pathogen, in the temperate climate of the North American Pacific Northwest (PNW) in 1999 prompted several questions. The most prevalent among these was the timing of the introduction of this pathogen to this novel environment. Here, we infer tip-dated timing estimates for the three clonal C. gattii populations observed in the PNW, VGIIa, VGIIb, and VGIIc, based on whole-genome sequencing of 134 C. gattii isolates and using Bayesian evolutionary analysis by sampling trees (BEAST). We estimated the nucleotide substitution rate for each lineage (1.59 × 10−8, 1.59 × 10−8, and 2.70 × 10−8, respectively) to be an order of magnitude higher than common neutral fungal mutation rates (2.0 × 10−9), indicating a microevolutionary rate (e.g., successive clonal generations in a laboratory) in comparison to a species’ slower, macroevolutionary rate (e.g., when using fossil records). The clonal nature of the PNW C. gattii emergence over a narrow number of years would therefore possibly explain our higher mutation rates. Our results suggest that the mean time to most recent common ancestor for all three sublineages occurred within the last 60 to 100 years. While the cause of C. gattii dispersal to the PNW is still unclear, our research estimates that the arrival is neither ancient nor very recent (i.e., <25 years ago), making a strong case for an anthropogenic introduction.

IMPORTANCE The recent emergence of the pathogenic fungus Cryptococcus gattii in the Pacific Northwest (PNW) resulted in numerous investigations into the epidemiological and enzootic impacts, as well as multiple genomic explorations of the three primary molecular subtypes of the fungus that were discovered. These studies lead to the general conclusion that the subtypes identified likely emerged out of Brazil. Here, we conducted genomic dating analyses to determine the ages of the various lineages seen in the PNW and propose hypothetical causes for the dispersal events. Bayesian evolutionary analysis strongly suggests that these independent fungal populations in the PNW are all 60 to 100 years old, providing a timing that is subsequent to the opening of the Panama Canal, which allowed for more direct shipping between Brazil and the western North American coastline, a possible driving event for these fungal translocation events.

INTRODUCTION

Combining epidemiology with microbial evolution analyses in a historical context is critical to understanding the nature of newly occurring infectious diseases. Diseases may emerge in a new region due to recent pathogen translocation events (e.g., West Nile virus in the United States in 1999 [1], Ebola virus in western Africa in 2014 [2], and Zika virus in Brazil in 2016 [3]). Many times “emerging diseases” are only emerging in our understanding of previously undetected endemic disease (e.g., Legionella in Philadelphia in 1976 [4] and hantavirus in the American Southwest in 1993 [5]). The appearance of Cryptococcus gattii in the Pacific Northwest (PNW) may represent both of these models (note that while the species nomenclature of C. gattii is currently under debate [6, 7], for consistency with previous and ongoing studies, we use the traditional C. gattii nomenclature herein, with a focus on the major molecular type VGII).

The source and timing of the emergence of C. gattii in the PNW have been a challenge to public health and mycology researchers since cryptococcosis seemingly first appeared in British Columbia in 1999 (8, 9). Early studies elucidated the highly clonal nature of the newly identified subtypes of the VGII major molecular type—VGIIa, VGIIb, and eventually, VGIIc—in the Oregon-Washington region (9, 10). More recent studies have identified the origin of C. gattii, including these subtypes, to be South America, likely Brazil (11–13), where C. gattii is endemic. However, their apparently sudden appearance, with novel phenotypes and relatively widespread nature, across the North American PNW remained a genomic enigma (12).

Here, we apply Bayesian analysis-based genomic dating to understand the timing of the PNW emergence(s) in order to better elucidate the causes of the C. gattii translocation events and dispersal to and within North America. In an effort to provide estimates of time to most recent common ancestor (TMRCA), we implemented an established dating method and combined that with public health surveillance and epidemiology, along with a historical understanding of global anthropogenic events, to establish a hypothesis of “how and when” for the dispersal of C. gattii.

RESULTS

A total of 134 Cryptococcus gattii whole genomes, 112 of which were previously published, representing both the global VGII lineages and the Pacific Northwest subtypes, were included in the maximum-likelihood tree (Fig. 1). This tree is based on 289,240 total single-nucleotide polymorphisms (SNPs) with 213,653 informative SNPs, using the TVM+ASC+G4 nucleotide substitution model (14). This SNP matrix had a quality breadth of coverage of 87.57% of the reference genome R265.

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

Maximum-likelihood phylogeny and geographic locations of the PNW Cryptococcus gattii VGII lineages. This tree is based on 289,240 total SNPs using the TVM+ASC+G4 nucleotide substitution model and includes 134 genomes that collectively cover 87.57% of the R265 BC human 2001 reference genome. Within this lineage, three clonal subclades exist, VGIIa, VGIIb, and VGIIc. This tree includes 1,000 bootstrap pseudoreplicates; nodes with bootstrap values less than 95% are denoted by red circles. Geographic representations of the three PNW clones are a generalized depiction of the primary impact region of each in the PNW.

The C. gattii subtypes in the PNW are highly clonal, with 412, 153, and 268 SNPs within VGIIa, VGIIb, and VGIIc, respectively (see Fig. S1 to S3 in the supplemental material), which is consistent with preliminary findings (9, 10) but unlike other C. gattii populations in the United States (15). The newly sequenced 2015/2016 isolates of VGIIa, VGIIb, and VGIIc all grouped within the PNW clades. Two 2015 autochthonous VGIIa isolates from Québec, Canada (16), that grouped together within the PNW clade, seemingly derived from an Oregon strain lineage, and one 2015 wildlife isolate from Nova Scotia (17) that grouped within the primary VGIIb PNW clade seemingly derived from a Washington lineage. The three PNW genotypes all continue to be clonal, with no evidence of recombination based on the pairwise homoplasy index (PHI) statistic.

FIG S1

Maximum-likelihood phylogeny of C. gattii VGIIa data set. The VGIIa tree is based on 412 SNPs using the K2P+ASC nucleotide substitution model and includes 29 genomes that collectively cover 96.9% of the R265 BC human 2001 reference genome. The PHI test did not find statistically significant evidence for recombination in the VGIIa SNP matrix (P = 1.0). Download FIG S1, JPG file, 1.1 MB.
Copyright © 2018 Roe et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

FIG S2

Maximum-likelihood phylogeny of C. gattii VGIIb data set. The VGIIb tree is based on 153 SNPs using the K2P+ASC nucleotide substitution model and includes 8 genomes with a quality breadth of coverage of the B11567 Nova Scotia deer VGIIb 2015 reference genome of 97.1%. The PHI test statistic, P = 1.0, indicates no significant evidence for recombination in the VGIIb SNP matrix. Download FIG S2, JPG file, 0.6 MB.
Copyright © 2018 Roe et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

FIG S3

Maximum-likelihood phylogeny of C. gattii VGIIc data set. The VGIIc tree included 268 SNPs using the K2P+ASC nucleotide substitution model. Here, 29 genomes had a quality breadth of coverage of the B11468 OR human 2016 reference genome of 96.4%. The PHI test again yielded P = 1.0. Download FIG S3, JPG file, 1.7 MB.
Copyright © 2018 Roe et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

The root-to-tip regressions identified various degrees of clocklike behavior among the three genotypes, with R2 values of 0.5971, 0.661, and 0.0745, suggesting that VGIIa and VGIIb have strong clocklike behavior while VGIIc has weak clocklike behavior (Fig. S4 to S6). However, because all three genotypes have positive regression slopes, molecular clock analyses are appropriate and reliable for mutation rate estimation (18). The best-fitting model from MEGA7 that was also available in the BEAST (Bayesian evolutionary analysis by sampling trees) software was implemented: the HKY model was applied to both the VGIIa and VGIIb data sets, while the best-fitting model for the VGIIc data set was the TN93 model.

FIG S4

Root-to-tip analyses to assess the temporal signal (the clocklike behavior) in the VGIIa data set of point mutation accumulation over time. The VGIIa SNPs illustrate strong clocklike behavior with an R2 of 0.5971. The 10,000-date randomization permutation distribution is included, showing a P value of 0.076, suggesting our regression analysis is on the cusp of being better than we would observe from random chance. Download FIG S4, JPG file, 0.2 MB.
Copyright © 2018 Roe et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

FIG S5

Root-to-tip analysis of the VGIIb group resulting in an R2 value of 0.661 and a positive slope, suggesting strong clocklike behavior. The 10,000-date randomization permutation distribution for the VGIIb data set had a P value of 0.268. Our observed R2 value is only better than 73% of what we observe from random chance. Download FIG S5, JPG file, 0.2 MB.
Copyright © 2018 Roe et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

FIG S6

The VGIIc regression plot is diffuse, with an R2 value of 0.074, suggesting weak clocklike behavior. The 10,000-date randomization permutation distribution for VGIIc produced a P value of 0.294. Our observed R2 value is only better than 70% of what we observe from random chance. Download FIG S6, JPG file, 0.2 MB.
Copyright © 2018 Roe et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

The 10,000-date randomization permutation testing on the VGIIa sample set produced a P value of 0.076, suggesting that the R2 value produced in the regression analysis was better than 9,240 date randomized regression analyses and that our results are significantly different from what would be expected from random chance (Fig. S4). The date randomization testing for VGIIb and VGIIc produced P values of 0.268 and 0.294, respectively, showing that our observed R2 values were not statistically different than random chance (Fig. S5 and S6).

The estimated mutation rates for VGIIa and VGIIc are extremely similar by BEAST calculations; the VGIIa mutation rate was calculated to be 1.59 × 10−8 SNPs per base per year (95% highest posterior density [HPD], 5.54 × 10−9 to 2.93 × 10−8), while that of VGIIc is 1.59 × 10−8 SNPs per base per year (95% HPD, 5.54 × 10−9 to 2.04 × 10−8). Even though VGIIa’s root-to-tip regression was clocklike compared to that of VGIIc, both genotypes exhibit low rates of evolutionary change. The estimated mutation rate for VGIIb was calculated to be nearly twice as high as those of the VGIIa and VGIIc lineages, with 2.70 × 10−8 SNPs per base per year (95% HPD, 6.75 × 10−9 to 5.23 × 10−8).

The time to most recent ancestor (TMRCA) was calculated for each sample set from the BEAST analysis; the mean TMRCA for the PNW VGIIa sample set was estimated at 87.99 years ago (95% HPD, 53.87 to 173.76), the mean TMRCA for the PNW VGIIb sample set was estimated at 81.43 years ago (95% HPD, 27.93, 187.95), and the mean TMRCA for the VGIIc sample set was estimated at 66.29 years ago (95% HPD, 26.91, 115.43) (Fig. 2).

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

Bayesian phylogenetic analyses of C. gattii samples. BEAST 1.8.4 was used to produce calibrated phylogenies with the mean TMRCA estimates, which were 87.99 years ago for the VGIIa data set (A), 81.43 years ago for the VGIIb data set (B), and 66.29 years ago for the VGIIc data set (C). The tips of the branches correspond to the year of sampling. Red nodes represent internal nodes with posterior probability support of >0.95.

DISCUSSION

Between 2005 and 2013, a total of 273 Cryptococcus gattii cases from both human and veterinary sources were reported in the PNW (19), with new cases still reported today. Understanding the emergence and continual evolution of this pathogen in a novel environment is critical to the understanding of the ongoing epidemiology of cryptococcal disease in this region and may be important to studying the evolution of other emerging pathogens. This unexpected and ongoing emergence of what was previously thought to be predominantly a tropical pathogen spurred several questions regarding the genetic and geographic origin of this fungus (12, 20) and, pointedly, the timing of the introduction of C. gattii into the PNW. Here, we applied a Bayesian approach to narrow in on a dispersal estimate.

Previous studies estimating the divergence time of C. neoformans and C. gattii applied a common neutral mutation rate of 2.0 × 10−9. We estimated the mean substitution rates for all three genotypes to be between 1.59 × 10−8 and 2.70 × 10−8, an order of magnitude higher than previously described rates (21). Several studies have shown microevolutionary rates (e.g., successive clonal generations in a laboratory) to be consistently higher than a species’ macroevolutionary rates (e.g., when using fossil records), which are likely affected by purifying selection that becomes evident over longer periods of time (22–27). The latter likely includes other nonclonal reproduction effects that would not impact the likely asexual life history of C. gattii in the PNW, where only one mating type has been identified (9). The clonal nature of the PNW emergences over a narrow number of years would explain the higher mutation rates.

Our results suggest that C. gattii has been in the PNW longer than previously hypothesized. The clinical emergence of specific major molecular types in 1999 (VGIIa and VGIIb) and 2005 (VGIIc) was likely an emergence in our understanding of previously arrived and dispersed fungi from Brazil. The polytomy effect of individual lineages dating back to the MRCA (i.e., founder populations) in each of the VGIIa, VGIIb, and VGIIc clades suggests an initial evolutionary bottleneck followed by an early intra-PNW dispersal event, and the subclades within these lineages demonstrate ongoing divergence events during their evolutionary histories on the continent. Unlike VGIIa and VGIIc, VGIIb was reported to have been introduced to the PNW on at least two occasions (12); in this analysis, we only included the members of the primary (i.e., major) PNW clade, which now includes the 2015 Nova Scotia isolate, as the other, minor clade contained only two isolate genomes with no new members identified since the original analysis (12). The increased attention following the initial case cluster of VGIIa on Vancouver Island in 1999 likely led to an increase of case findings and the subsequent finding of additional “novel” subtypes elsewhere in the region (e.g., VGIIc in Oregon) (9). Our common centurial dating for the seemingly independent emergences of the PNW C. gattii subtypes suggests a possible common mechanism for their translocation. A search for causes of transmission can now focus on events and means of transmission that occurred at the time of or shortly preceding the MRCA dates.

The Teddy Roosevelt effect.Given the commonality of just under 100 years for the PNW arrival of the VGIIa, -b, and -c subtypes and given their lack of known presence elsewhere on the North American continent (with three notable exceptions: a single Florida VGIIb isolate [28], the recently discovered Québec cluster of VGIIa [16, 29], which seems to have derived from the PNW VGIIa population, and the single VGIIb isolate from a deer in Nova Scotia [17], likely derived from the PNW VGIIb population), it is logical to hypothesize that these C. gattii subtypes arrived in the PNW region shortly into the 20th century. It is possible that populations arrived elsewhere on the continent and failed to become established; a similar scenario was seen with the concurrent spread of Yersinia pestis by boats from Asia to multiple U.S. ports around the turn of the 20th century that only became enzootic through the ports in California (30). What appears to be clear is that there was an outward spread of the PNW C. gattii populations from the PNW coastal regions to locales further inland (including across the country of Canada). These findings would suggest at a minimum that a translocation event or events occurred between Brazil and the PNW coast approximately 100 years ago. As shipping was minimal between Brazil and Pacific Coastal countries in South America until the opening of the Panama Canal, it follows that this historical event may have allowed for ships to carry cargo, and microbial stowaways, between these distant locations for the first time (Fig. 3).

FIG. 3
  • Open in new tab
  • Download powerpoint
FIG. 3

Geographic representation of the timing of dispersal of the VGII lineages to the PNW. ya, years ago. (Adapted from reference 60 with permission from National Geographic.)

The Panamanian isthmus acts as a land bridge between the American continents and as a natural barrier between the Pacific and Atlantic Ocean. The isthmus has played a crucial role in the geology, hydrology, biology, and climate in the region since the uprising of this intercontinental land bridge nearly 3 million years ago (31). The isthmus has also played a key role in the dispersal of microbes: notably, the fungal pathogen Coccidioides posadasii likely was translocated from North America to South America across the isthmus via mass movements of mammals during the later Great American Biotic Interchange events that occurred nearly a million years ago (32). After the final rising of the isthmus, a permanent land barrier was formed between the continents, separating the oceans, until just over 100 years ago. In 1914, after years of human toil, the Panama Canal, commissioned by Theodore Roosevelt, was completed. This allowed, for the first time in nearly 3 million years, water movement between the Atlantic and Pacific Oceans, subsequently allowing for cargo ships to move between the east and west coasts of the Western Hemisphere continents without having to go through dangerous waters at the southern point of South America. After the opening of the Panama Canal, shipping between ports of Brazil and western North America began in earnest.

Common goods transported out of Brazil at this time were Brazilian hardwood lumber, minerals, coffee, and rubber. While C. gattii is known to have been transported on (33) and to thrive on (34, 35) live trees and seeds, no contamination of lumber products has been documented previously. While we may be unable to ascertain a specific shipping product that may have carried the fungus from Brazil to the PNW and elsewhere, one hypothesis we propose is contaminated ballast water. Ballast water contamination is a well-understood problem in large vessels, and ballast has been known to contain and transport nonnative animals, algae, and microbes around the world (36). At a vessel’s source port, seawater is taken on by ships and stored in their ballast to provide vessel stability. This ballast water is released en route or at the destination port to correct for fluctuations in cargo weight (37). Cryptococcus has been shown to survive in vitro in saline (8), has been isolated from seawater (38, 39), and is known to cause significant morbidity and mortality in marine mammals (40–42), particularly in the PNW (43–45). Again, while only hypothetical, it is possible that a limited number of ships with contaminated ballast water could have transported these pathogens unknowingly to PNW ports and beyond. It is notable that the dominant PNW clade of VGIIb also shares a very recent common ancestor with multiple lineages in Southeast Asia. Again, a shipping vessel from Brazil carrying goods (and contaminants) to the PNW could also have caused the translocation to Asia at this time, given that the new transcontinental seaway would have allowed shipping to these markets as well.

An alternate and intriguing hypothesis for dispersal to Malaysia and elsewhere in Southeast Asia is the possible contamination of rubber tree (Hevea brasiliensis) seeds, which were infamously “stolen” and exported to Britain by Henry Wickham in 1876 to establish British-controlled rubber tree plantations in Ceylon (now Sri Lanka) and Malay (modern day Malaysia and Singapore) and, subsequently, to British colonies in Africa and India (46). While not previously shown to be associated with Hevea seeds, C. gattii contamination of other tree seeds, such as Eucalyptus spp., has been well documented (33). Brazil, attempting to maintain global control of the rubber trade, tried to prevent the shipment of such seeds, limiting other possible such exports out of Brazil; however, the Wickham seed transport likely resulted in the crash of Brazilian rubber beginning in 1910 (46). The 1876 seed transport timeline and global movement of Hevea seeds at this time could represent an indirect dispersal mechanism. Along these lines, the Brazil nut (Bertholletia excelsa) has been an additional plant-based export out of Brazil around the world for multiple centuries; however, intense exportations began following the Brazilian rubber crash in 1910 (47), providing an appropriate temporal connection to the C. gattii dispersal event. Extensive analyses of well-known ongoing fungal contamination of Brazil nuts during collection and processing have not identified the presence of Cryptococcus spp. (47), and there is no obvious route to localized environmental contamination in the PNW from these products.

While the prevalence of recombination along with the extent of temporal signal within a sample set plays a large role in estimating genome-wide evolutionary rates (18), these sample sets show no statistical evidence of recombination to hinder these analyses. However, the date randomization permutation testing showed there was little true temporal signal within the VGIIb and VGIIc sample sets. Given the narrow sampling dates and the limited number of samples in each data set, it is possible we are only capturing a glimpse of the temporal signal and, with a wider sampling time and more samples, it would be possible to narrow our MRCA estimates.

Whatever the cause of C. gattii dispersal to the PNW, it is clear that those populations are neither ancient nor very recent (i.e., <25 years ago) arrivals to the region. Dispersal in the last 100 years would strongly suggest anthropogenic causes, directly or indirectly, rather than movement by natural animal migrations, as suggested with Cryptococcus populations (48) and other environmental pathogenic fungi (32). The identification of VGIIa in a pet store bird (sample number B11566) and a separate, unrelated case in pet store worker (sample number B11569), both in Québec, are unexpected occurrences (29), and it is of interest to note that the bird reportedly originated from British Columbia and the infecting C. gattii strain clearly evolved from the PNW VGIIa population, suggesting a linkage to human transport of animals as the cause. Conversely, the appearance of a PNW-derived strain of VGIIb (sample number B11567) in a Nova Scotia yearling white-tailed deer (Odocoileus virginianus) (17) is less explainable. There is no other current evidence of long-distance dispersal of VGIIb across the North American continent, nor is there evidence of transport of animals from the PNW to the region, as the infected deer was part of a local deer population. White-tailed deer are nonmigratory, providing further evidence of the presence of VGIIb in Nova Scotia, and this is suggestive of a unique introduction/dispersal event in the region (17). The province of Nova Scotia is a peninsula that is surrounded by four major bodies of water, and the city of Halifax, throughout its history, has been an active sea port, one of the major international points of entry for shipping into Atlantic, Canada. More investigation of the local environment is necessary to understand the history of this incident and provide estimates of regional endemicity. As the North American populations of C. gattii continue to evolve and disperse, it will be useful to continually apply genomic dating to understand the nature of these events and the expanding impact of these fungi on human and veterinary health.

MATERIALS AND METHODS

Samples.A total of 66 genomes were included in the analysis of the three clonal lineages found in the PNW. The 22 new genomes from this study include 13 C. gattii recent isolates collected either in 2015 or 2016 from Washington, Oregon, and Canada, which were comprised of 7 VGIIa, 1 VGIIb, and 5 VGIIc isolates (Table 1). Previously published isolate genomes (12) were included for each genotype, including 22 VGIIa, 16 VGIIb, and 24 VGIIc genomes (Table 1). The VGII lineage tree was based off of 134 genomes, 112 of which were previously published (12), including the above-mentioned samples (Table 1). The recently described msh2 variant lineage of VGIIa (49) was not included in the genotype-specific genomic dating analyses (described below) due to possible confounding from the hypermutation phenotype, except for the NIH444 genome, which had no identifiable impacts on the total number of SNP mutations (12).

View this table:
  • View inline
  • View popup
TABLE 1

C. gattii VGII sample genomes used in this study, including 22 newly sequenced and 112 previously published genomes

Genomic DNA was extracted from the 22 isolates using the ZR fungal/bacterial DNA MiniPrep kit (Zymo Research), following the manufacturer’s instructions. DNA samples were fragmented by sonication and prepared for multiplexed, paired-end sequencing with a 700-bp insert using the library hyper-preparation kit with standard PCR library amplification (KAPA Biosystems) as previously described (32). Libraries were quantified using the KAPA library quantification kit (KAPA Biosystems) and sequenced to a read length of 300 bp on the MiSeq instrument (Illumina, Inc., San Diego, CA).

SNP matrix generation.Sequenced genomes were assembled de novo using SPAdes version 3.6.0 (50). Reference genomes for phylogenetic analyses were selected based on the assembly quality metric N50 and total assembly length (see Table S1 in the supplemental material). NASP (51) was used to generate an SNP matrix for each of the three C. gattii subtypes, as well as the VGII complex, using the selected reference genomes for each genotype and R265 for the VGII lineage analysis. In brief, sequencing reads were aligned to the reference genome using Novoalign (Novocraft Technologies Sdn Bhd) and SNPs were identified using the GATK Unified Genotyper ToolKit version 2.7-2 (52). NASP filtered SNP loci not present in every sample and with less than 10× coverage and less than a 90% consensus in any sample. Using NUCmer (53), duplicated regions identified within the reference genome were also removed.

TABLE S1

Reference genome assembly statistics, including number of contigs, number of contigs larger than 500 bases, N50, largest contig, and total size. Genomes were assembled using SPAdes. Download TABLE S1, XLSX file, 0.03 MB.
Copyright © 2018 Roe et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Phylogenetic analyses and root-to-tip regression.IQ-Tree version 1.3.10 (54) was used to infer the best-fitting nucleotide substitution model for each of the three SNP matrices and to produce maximum-likelihood trees with 1,000 nonparametric bootstrap pseudoreplicates for branch confidence. Trees were visualized in FigTree version 1.4.2, and SNP numbers placed on branches using phangorn (55). To test each SNP matrix for evidence of recombination, which can confound divergence-dating analyses (18), the PHI test was conducted using PHIPack (56). In order to assess the temporal signal of the data set, regression analysis implementing root-to-tip genetic distance as a function of the sample collection year was conducted using the software package TempEst version 1.5 (http://tree.bio.ed.ac.uk/software/tempest/). Using the determination coefficient, R2, a measure of clocklike behavior was assessed. In an effort to maximize R2, the best-fitting root was selected based on TempEst recommendations (18). Additionally, we performed 10,000 random permutations of the sampling dates over the sequences in an effort to evaluate the significance of our regression results (57).

Divergence time analyses.A Bayesian molecular clock using tip dating was implemented in the BEAST version 1.8.0 software package (58) to infer evolutionary rates and time to most recent common ancestor (TMRCA) for the three VGII genotypes. MEGA7 (14) was used for nucleotide substitution model selection for each genotype, utilizing the Bayesian information criterion results to determine the best-fitting models. As previously described (32), available nucleotide substitution models are limited in BEAST, and therefore, the best-fitting model from MEGA7 that was also available in the BEAST software was implemented. While the SNP matrices only included variable sites, we corrected for the invariant sites by specifying a Constant Patterns model in the BEAST XML file. For each separate analysis, we determined the numbers of constant A’s, C’s, T’s, and G’s and added them to the XML file. Additionally, “path and stepping stone” sampling marginal-likelihood estimators were used in order to determine the best-fitting clock and demographic model combinations (59). The statistical fits of 10 different clock and demographic model combinations were assessed using the log marginal likelihood (Table S2). In BEAST, four independent chains of 4 billion iterations each were run for all molecular clock and demographic model combinations for the VGIIa and VGIIc data sets. Convergence among the four chains for VGIIb analysis completed after 1 billion iterations.

TABLE S2

Metrics from BEAST analysis of VGIIa, VGIIb, and VGIIc SNP matrices. Metrics are provided for both relaxed lognormal clock and relaxed clock. Download TABLE S2, XLSX file, 0.04 MB.
Copyright © 2018 Roe et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Accession number(s).All new sequence data files were deposited in the NCBI Sequence Read Archive (BioProject accession number PRJNA388113).

ACKNOWLEDGMENTS

We thank all our colleagues from around the world who have contributed isolates to either the current or previously published Cryptococcus whole-genome studies.

This analysis was supported by the Centers for Disease Control and Prevention (contract no. 200201461029 to D.M.E.), NIH (grant number R21AI098059 to D.M.E., G.R.T., S.R.L. and W.M.), and an Australian NH&MRC grant (grant number APP1031943 to W.M.).

The findings and conclusions in the manuscript are those of the authors and do not necessarily represent the official views of the Centers for Disease Control and Prevention.

FOOTNOTES

    • Received October 27, 2017.
    • Accepted December 13, 2017.
  • Copyright © 2018 Roe et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

REFERENCES

  1. 1.↵
    1. Lanciotti RS,
    2. Roehrig JT,
    3. Deubel V,
    4. Smith J,
    5. Parker M,
    6. Steele K,
    7. Crise B,
    8. Volpe KE,
    9. Crabtree MB,
    10. Scherret JH,
    11. Hall RA,
    12. MacKenzie JS,
    13. Cropp CB,
    14. Panigrahy B,
    15. Ostlund E,
    16. Schmitt B,
    17. Malkinson M,
    18. Banet C,
    19. Weissman J,
    20. Komar N,
    21. Savage HM,
    22. Stone W,
    23. McNamara T,
    24. Gubler DJ
    . 1999. Origin of the West Nile virus responsible for an outbreak of encephalitis in the Northeastern United States. Science 286:2333–2337. doi:10.1126/science.286.5448.2333.
    OpenUrlAbstract/FREE Full Text
  2. 2.↵
    1. Spengler JR,
    2. Ervin ED,
    3. Towner JS,
    4. Rollin PE,
    5. Nichol ST
    . 2016. Perspectives on West Africa Ebola virus disease outbreak, 2013–2016. Emerg Infect Dis 22:956–963. doi:10.3201/eid2206.160021.
    OpenUrlCrossRef
  3. 3.↵
    1. Armstrong N,
    2. Hou W,
    3. Tang Q
    . 2017. Biological and historical overview of Zika virus. J Virol 6:1–8. doi:10.5501/wjv.v6.i1.1.
    OpenUrlCrossRef
  4. 4.↵
    1. Mercante JW,
    2. Morrison SS,
    3. Desai HP,
    4. Raphael BH,
    5. Winchell JM
    . 2016. Genomic analysis reveals novel diversity among the 1976 Philadelphia Legionnaires’ disease outbreak isolates and additional ST36 strains. PLoS One 11:e0164074. doi:10.1371/journal.pone.0164074.
    OpenUrlCrossRef
  5. 5.↵
    1. Hjelle B,
    2. Torres-Pérez F
    . 2010. Hantaviruses in the Americas and their role as emerging pathogens. Viruses 2:2559–2586. doi:10.3390/v2122559.
    OpenUrlCrossRefPubMedWeb of Science
  6. 6.↵
    1. Hagen F,
    2. Khayhan K,
    3. Theelen B,
    4. Kolecka A,
    5. Polacheck I,
    6. Sionov E,
    7. Falk R,
    8. Parnmen S,
    9. Lumbsch HT,
    10. Boekhout T
    . 2015. Recognition of seven species in the Cryptococcus gattii/Cryptococcus neoformans species complex. Fungal Genet Biol 78:16–48. doi:10.1016/j.fgb.2015.02.009.
    OpenUrlCrossRef
  7. 7.↵
    1. Kwon-Chung KJ,
    2. Bennett JE,
    3. Wickes BL,
    4. Meyer W,
    5. Cuomo CA,
    6. Wollenburg KR,
    7. Bicanic TA,
    8. Castañeda E,
    9. Chang YC,
    10. Chen J,
    11. Cogliati M,
    12. Dromer F,
    13. Ellis D,
    14. Filler SG,
    15. Fisher MC,
    16. Harrison TS,
    17. Holland SM,
    18. Kohno S,
    19. Kronstad JW,
    20. Lazera M,
    21. Levitz SM,
    22. Lionakis MS,
    23. May RC,
    24. Ngamskulrongroj P,
    25. Pappas PG,
    26. Perfect JR,
    27. Rickerts V,
    28. Sorrell TC,
    29. Walsh TJ,
    30. Williamson PR,
    31. Xu J,
    32. Zelazny AM,
    33. Casadevall A
    . 2017. The case for adopting the “species complex” nomenclature for the etiologic agents of cryptococcosis. mSphere 2:e00357-16. doi:10.1128/mSphere.00357-16.
    OpenUrlAbstract/FREE Full Text
  8. 8.↵
    1. Kidd SE,
    2. Chow Y,
    3. Mak S,
    4. Bach PJ,
    5. Chen H,
    6. Hingston AO,
    7. Kronstad JW,
    8. Bartlett KH
    . 2007. Characterization of environmental sources of the human and animal pathogen Cryptococcus gattii in British Columbia, Canada, and the Pacific Northwest of the United States. Appl Environ Microbiol 73:1433–1443. doi:10.1128/AEM.01330-06.
    OpenUrlAbstract/FREE Full Text
  9. 9.↵
    1. Byrnes EJ, III,
    2. Bildfell RJ,
    3. Frank SA,
    4. Mitchell TG,
    5. Marr KA,
    6. Heitman J
    . 2009. Molecular evidence that the range of the Vancouver Island outbreak of Cryptococcus gattii infection has expanded into the Pacific Northwest in the United States. J Infect Dis 199:1081–1086. doi:10.1086/597306.
    OpenUrlCrossRefPubMedWeb of Science
  10. 10.↵
    1. Gillece JD,
    2. Schupp JM,
    3. Balajee SA,
    4. Harris J,
    5. Pearson T,
    6. Yan Y,
    7. Keim P,
    8. DeBess E,
    9. Marsden-Haug N,
    10. Wohrle R,
    11. Engelthaler DM,
    12. Lockhart SR
    . 2011. Whole genome sequence analysis of Cryptococcus gattii from the Pacific Northwest reveals unexpected diversity. PLoS One 6:e28550. doi:10.1371/journal.pone.0028550.
    OpenUrlCrossRefPubMed
  11. 11.↵
    1. Hagen F,
    2. Ceresini PC,
    3. Polacheck I,
    4. Ma H,
    5. van Nieuwerburgh F,
    6. Gabaldón T,
    7. Kagan S,
    8. Pursall ER,
    9. Hoogveld HL,
    10. van Iersel LJJ,
    11. Klau GW,
    12. Kelk SM,
    13. Stougie L,
    14. Bartlett KH,
    15. Voelz K,
    16. Pryszcz LP,
    17. Castañeda E,
    18. Lazera M,
    19. Meyer W,
    20. Deforce D,
    21. Meis JF,
    22. May RC,
    23. Klaassen CHW,
    24. Boekhout T
    . 2013. Ancient dispersal of the human fungal pathogen Cryptococcus gattii from the Amazon rainforest. PLoS One 8:e71148. doi:10.1371/journal.pone.0071148.
    OpenUrlCrossRefPubMed
  12. 12.↵
    1. Engelthaler DM,
    2. Hicks ND,
    3. Gillece JD,
    4. Roe CC,
    5. Schupp JM,
    6. Driebe EM,
    7. Gilgado F,
    8. Carriconde F,
    9. Trilles L,
    10. Firacative C,
    11. Ngamskulrungroj P,
    12. Castañeda E,
    13. Lazera MDS,
    14. Melhem MSC,
    15. Pérez-Bercoff A,
    16. Huttley G,
    17. Sorrell TC,
    18. Voelz K,
    19. May RC,
    20. Fisher MC,
    21. Thompson GR,
    22. Lockhart SR,
    23. Keim P,
    24. Meyer W
    . 2014. Cryptococcus gattii in North American Pacific Northwest: whole-population genome analysis provides insights into species evolution and dispersal. mBio 5:e01464-14. doi:10.1128/mBio.01464-14.
    OpenUrlAbstract/FREE Full Text
  13. 13.↵
    1. Souto ACP,
    2. Bonfietti LX,
    3. Ferreira-Paim K,
    4. Trilles L,
    5. Martins M,
    6. Ribeiro-Alves M,
    7. Pham CD,
    8. Martins L,
    9. Dos Santos W,
    10. Chang M,
    11. Brito-Santos F,
    12. Santos DCS,
    13. Fortes S,
    14. Lockhart SR,
    15. Wanke B,
    16. Melhem MSC,
    17. Lazéra MS,
    18. Meyer W
    . 2016. Population genetic analysis reveals a high genetic diversity in the Brazilian Cryptococcus gattii VGII population and shifts the global origin from the Amazon rainforest to the semi-arid desert in the northeast of Brazil. PLoS Negl Trop Dis 10:e0004885. doi:10.1371/journal.pntd.0004885.
    OpenUrlCrossRef
  14. 14.↵
    1. Kumar S,
    2. Stecher G,
    3. Tamura K
    . 2016. MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets. Mol Biol Evol 33:1870–1874. doi:10.1093/molbev/msw054.
    OpenUrlCrossRefPubMed
  15. 15.↵
    1. Lockhart SR,
    2. Roe CC,
    3. Engelthaler DM
    . 2016. Whole-genome analysis of Cryptococcus gattii, Southeastern United States. Emerg Infect Dis 22:1098–1101. doi:10.3201/eid2206.151455.
    OpenUrlCrossRefPubMed
  16. 16.↵
    1. St-Pierre J,
    2. Dufresne PJ,
    3. Carignan A,
    4. Levesque E,
    5. Bernard F,
    6. Longtin J,
    7. LeBlanc L
    . 30 October 2017. Report of the first two human indigenous cases of Cryptococcus gattii infection in eastern Canada. Mycopathologia doi:10.1007/s11046-017-0215-8.
    OpenUrlCrossRef
  17. 17.↵
    1. Overy DP,
    2. McBurney S,
    3. Muckle A,
    4. Lund L,
    5. Lewis PJ,
    6. Strang R
    . 2016. Cryptococcus gattii VGIIb-like variant in white-tailed deer, Nova Scotia, Canada. Emerg Infect Dis 22:1131–1133. doi:10.3201/eid2206.160081.
    OpenUrlCrossRef
  18. 18.↵
    1. Duchêne S,
    2. Holt KE,
    3. Weill FX,
    4. Le Hello S,
    5. Hawkey J,
    6. Edwards DJ,
    7. Fourment M,
    8. Holmes EC
    . 2016. Genome-scale rates of evolutionary change in bacteria. Microb Genom 2:e000094. doi:10.1099/mgen.0.000094.
    OpenUrlCrossRef
  19. 19.↵
    1. Lockhart SR,
    2. Iqbal N,
    3. Harris JR,
    4. Grossman NT,
    5. DeBess E,
    6. Wohrle R,
    7. Marsden-Haug N,
    8. Vugia DJ
    . 2013. Cryptococcus gattii in the United States: genotypic diversity of human and veterinary isolates. PLoS One 8:e74737. doi:10.1371/journal.pone.0074737.
    OpenUrlCrossRefPubMed
  20. 20.↵
    1. Billmyre RB,
    2. Croll D,
    3. Li W,
    4. Mieczkowski P,
    5. Carter DA,
    6. Cuomo CA,
    7. Kronstad JW,
    8. Heitman J
    . 2014. Highly recombinant VGII Cryptococcus gattii population develops clonal outbreak clusters through both sexual macroevolution and asexual microevolution. mBio 5:e01494-14. doi:10.1128/mBio.01494-14.
    OpenUrlAbstract/FREE Full Text
  21. 21.↵
    1. Sharpton TJ,
    2. Neafsey DE,
    3. Galagan JE,
    4. Taylor JW
    . 2008. Mechanisms of intron gain and loss in Cryptococcus. Genome Biol 9:R24. doi:10.1186/gb-2008-9-1-r24.
    OpenUrlCrossRefPubMed
  22. 22.↵
    1. Parsons TJ,
    2. Muniec DS,
    3. Sullivan K,
    4. Woodyatt N,
    5. Alliston-Greiner R,
    6. Wilson MR,
    7. Berry DL,
    8. Holland KA,
    9. Weedn VW,
    10. Gill P,
    11. Holland MM
    . 1997. A high observed substitution rate in the human mitochondrial DNA control region. Nat Genet 15:363–368. doi:10.1038/ng0497-363.
    OpenUrlCrossRefPubMedWeb of Science
  23. 23.↵
    1. Howell N,
    2. Smejkal CB,
    3. Mackey DA,
    4. Chinnery PF,
    5. Turnbull DM,
    6. Herrnstadt C
    . 2003. The pedigree rate of sequence divergence in the human mitochondrial genome: there is a difference between phylogenetic and pedigree rates. Am J Hum Genet 72:659–670. doi:10.1086/368264.
    OpenUrlCrossRefPubMedWeb of Science
  24. 24.↵
    1. Santos C,
    2. Montiel R,
    3. Sierra B,
    4. Bettencourt C,
    5. Fernandez E,
    6. Alvarez L,
    7. Lima M,
    8. Abade A,
    9. Aluja MP
    . 2005. Understanding differences between phylogenetic and pedigree-derived mtDNA mutation rate: a model using families from the Azores Islands (Portugal). Mol Biol Evol 22:1490–1505. doi:10.1093/molbev/msi141.
    OpenUrlCrossRefPubMedWeb of Science
  25. 25.↵
    1. Gibbs AJ,
    2. Fargette D,
    3. García-Arenal F,
    4. Gibbs MJ
    . 2010. Time—the emerging dimension of plant virus studies. J Gen Virol 91:13–22. doi:10.1099/vir.0.015925-0.
    OpenUrlCrossRefPubMedWeb of Science
  26. 26.↵
    1. Ho SY,
    2. Lanfear R,
    3. Bromham L,
    4. Phillips MJ,
    5. Soubrier J,
    6. Rodrigo AG,
    7. Cooper A
    . 2011. Time-dependent rates of molecular evolution. Mol Ecol 20:3087–3101. doi:10.1111/j.1365-294X.2011.05178.x.
    OpenUrlCrossRefPubMedWeb of Science
  27. 27.↵
    1. Duchêne DA,
    2. Hua X,
    3. Bromham L
    . 2017. Phylogenetic estimates of diversification rate are affected by molecular rate variation. J Evol Biol 30:1884–1897. doi:10.1111/jeb.13148.
    OpenUrlCrossRef
  28. 28.↵
    1. Kunadharaju R,
    2. Choe U,
    3. Harris JR,
    4. Lockhart SR,
    5. Greene JN
    . 2013. Cryptococcus gattii, Florida, USA, 2011. Emerg Infect Dis 19:519–521. doi:10.3201/eid1903.121399.
    OpenUrlCrossRefPubMed
  29. 29.↵
    1. Maccolini ÉO,
    2. Dufresne PJ,
    3. Aschenbroich SA,
    4. McHale B,
    5. Fairbrother JH,
    6. Bédard C,
    7. Hébert JA
    . 2017. A disseminated Cryptococcus gattii VGIIa infection in a citron-crested cockatoo (Cacatua sulphurea citrinocristata) in Québec, Canada. J Avian Med Surg 31:142–151. doi:10.1647/2016-193.
    OpenUrlCrossRef
  30. 30.↵
    1. Kugeler KJ,
    2. Staples JE,
    3. Hinckley AF,
    4. Gage KL,
    5. Mead PS
    . 2015. Epidemiology of human plague in the United States, 1900–2012. Emerg Infect Dis 21:16–22. doi:10.3201/eid2101.140564.
    OpenUrlCrossRefPubMed
  31. 31.↵
    1. O’Dea A,
    2. Lessios HA,
    3. Coates AG,
    4. Eytan RI,
    5. Restrepo-Moreno SA,
    6. Cione AL,
    7. Collins LS,
    8. de Queiroz A,
    9. Farris DW,
    10. Norris RD,
    11. Stallard RF,
    12. Woodburne MO,
    13. Aguilera O,
    14. Aubry MP,
    15. Berggren WA,
    16. Budd AF,
    17. Cozzuol MA,
    18. Coppard SE,
    19. Duque-Caro H,
    20. Finnegan S,
    21. Gasparini GM,
    22. Grossman EL,
    23. Johnson KG,
    24. Keigwin LD,
    25. Knowlton N,
    26. Leigh EG,
    27. Leonard-Pingel JS,
    28. Marko PB,
    29. Pyenson ND,
    30. Rachello-Dolmen PG,
    31. Soibelzon E,
    32. Soibelzon L,
    33. Todd JA,
    34. Vermeij GJ,
    35. Jackson JBC
    . 2016. Formation of the Isthmus of Panama. Sci Adv 2:e1600883. doi:10.1126/sciadv.1600883.
    OpenUrlFREE Full Text
  32. 32.↵
    1. Engelthaler DM,
    2. Roe CC,
    3. Hepp CM,
    4. Teixeira M,
    5. Driebe EM,
    6. Schupp JM,
    7. Gade L,
    8. Waddell V,
    9. Komatsu K,
    10. Arathoon E,
    11. Logemann H,
    12. Thompson GR,
    13. Chiller T,
    14. Barker B,
    15. Keim P,
    16. Litvintseva AP
    . 2016. Local population structure and patterns of Western Hemisphere dispersal for Coccidioides spp., the fungal cause of valley fever. mBio 7:e00550-16. doi:10.1128/mBio.00550-16.
    OpenUrlAbstract/FREE Full Text
  33. 33.↵
    1. Chakrabarti A,
    2. Jatana M,
    3. Kumar P,
    4. Chatha L,
    5. Kaushal A,
    6. Padhye AA
    . 1997. Isolation of Cryptococcus neoformans var. gattii from Eucalyptus camaldulensis in India. J Clin Microbiol 35:3340–3342.
    OpenUrlAbstract/FREE Full Text
  34. 34.↵
    1. Hiremath SS,
    2. Chowdhary A,
    3. Kowshik T,
    4. Randhawa HS,
    5. Sun S,
    6. Xu J
    . 2008. Long-distance dispersal and recombination in environmental populations of Cryptococcus neoformans var. grubii from India. Microbiology 154:1513–1524. doi:10.1099/mic.0.2007/015594-0.
    OpenUrlCrossRefPubMedWeb of Science
  35. 35.↵
    1. Springer DJ,
    2. Mohan R,
    3. Heitman J
    . 2017. Plants promote mating and dispersal of the human pathogenic fungus Cryptococcus. PLoS One 12:e0171695. doi:10.1371/journal.pone.0171695.
    OpenUrlCrossRef
  36. 36.↵
    1. Xu J,
    2. Wickramarathne TL,
    3. Chawla NV,
    4. Grey EK,
    5. Steinhaeuser K,
    6. Keller RP,
    7. Drake JM,
    8. Lodge DM
    . 2014. Improving management of aquatic invasions by integrating shipping network, ecological, and environmental data, p 1699–1708. In. Proceedings of the 20th ACM SIGKDD International Conference on Knowledge Discovery and Data Mining—KDD 14. ACM Press, New York, NY.
  37. 37.↵
    1. DiBacco C,
    2. Humphrey DB,
    3. Nasmith LE,
    4. Levings CD
    . 2012. Ballast water transport of non-indigenous zooplankton to Canadian ports. ICES J Mar Sci 69:483–491. doi:10.1093/icesjms/fsr133.
    OpenUrlCrossRef
  38. 38.↵
    1. Kidd SE,
    2. Bach PJ,
    3. Hingston AO,
    4. Mak S,
    5. Chow Y,
    6. MacDougall L,
    7. Kronstad JW,
    8. Bartlett KH
    . 2007. Cryptococcus gattii dispersal mechanisms, British Columbia, Canada. Emerg Infect Dis 13:51–57. doi:10.3201/eid1301.060823.
    OpenUrlCrossRefPubMed
  39. 39.↵
    1. Chang CF,
    2. Lee CF,
    3. Lin KY,
    4. Liu SM
    . 2016. Diversity of yeasts associated with the sea surface microlayer and underlying water along the northern coast of Taiwan. Res Microbiol 167:35–45. doi:10.1016/j.resmic.2015.08.005.
    OpenUrlCrossRef
  40. 40.↵
    1. Miller WG,
    2. Padhye AA,
    3. van Bonn W,
    4. Jensen E,
    5. Brandt ME,
    6. Ridgway SH
    . 2002. Cryptococcosis in a bottlenose dolphin (Tursiops truncatus) caused by Cryptococcus neoformans var. gattii. J Clin Microbiol 40:721–724. doi:10.1128/JCM.40.2.721-724.2002.
    OpenUrlAbstract/FREE Full Text
  41. 41.↵
    1. Rotstein DS,
    2. West K,
    3. Levine G,
    4. Lockhart SR,
    5. Raverty S,
    6. Morshed MG,
    7. Rowles T
    . 2010. Cryptococcus gattiivgi in a spinner dolphin (Stenella longirostris) from Hawaii. J Zoo Wildl Med 41:181–183. doi:10.1638/2009-0145.1.
    OpenUrlCrossRefPubMed
  42. 42.↵
    1. Huckabone SE,
    2. Gulland FMD,
    3. Johnson SM,
    4. Colegrove KM,
    5. Dodd EM,
    6. Pappagianis D,
    7. Dunkin RC,
    8. Casper D,
    9. Carlson EL,
    10. Sykes JE,
    11. Meyer W,
    12. Miller MA
    . 2015. Coccidioidomycosis and other systemic mycoses of marine mammals stranding along the central California, USA coast: 1998–2012. J Wildl Dis 51:295–308. doi:10.7589/2014-06-143.
    OpenUrlCrossRef
  43. 43.↵
    1. Norman SA,
    2. Raverty S,
    3. Zabek E,
    4. Etheridge S,
    5. Ford JKB,
    6. Hoang LMN,
    7. Morshed M
    . 2011. Maternal–fetal transmission of Cryptococcus gattii in harbor porpoise. Emerg Infect Dis 17:304–305. doi:10.3201/eid1702.101232.
    OpenUrlCrossRefPubMed
  44. 44.↵
    1. Simeone CA,
    2. Gulland FMD,
    3. Norris T,
    4. Rowles TK
    . 2015. A systematic review of changes in marine mammal health in North America, 1972–2012: the need for a novel integrated approach. PLoS One 10:e0142105. doi:10.1371/journal.pone.0142105.
    OpenUrlCrossRef
  45. 45.↵
    1. Rosenberg JF,
    2. Haulena M,
    3. Hoang LMN,
    4. Morshed M,
    5. Zabek E,
    6. Raverty SA
    . 2016. Cryptococcus gattii type VGIIa infection in harbor seals (Phoca vitulina) in British Columbia, Canada. J Wildl Dis 52:677–681. doi:10.7589/2015-11-299.
    OpenUrlCrossRef
  46. 46.↵
    1. Musgrave T,
    2. Musgrave W
    . 2002. An empire of plants: people and plants that changed the world, p 173. Cassell, London, United Kingdom.
  47. 47.↵
    1. Freitas-Silva O,
    2. Venâncioa A
    . 2011. Brazil nuts: benefits and risks associated with contamination by fungi and mycotoxins. Food Res Int 44:1434–1440. doi:10.1016/j.foodres.2011.02.047.
    OpenUrlCrossRef
  48. 48.↵
    1. Dixit A,
    2. Carroll SF,
    3. Qureshi ST
    . 2009. Cryptococcus gattii: an emerging cause of fungal disease in North America. Interdiscip Perspect Infect Dis 2009:840452. doi:10.1155/2009/840452.
    OpenUrlCrossRefPubMed
  49. 49.↵
    1. Billmyre RB,
    2. Clancey SA,
    3. Heitman J
    . 2017. Natural mismatch repair mutations mediate phenotypic diversity and drug resistance in Cryptococcus deuterogattii. eLife 6:e28802. doi:10.7554/eLife.28802.
    OpenUrlCrossRef
  50. 50.↵
    1. Bankevich A,
    2. Nurk S,
    3. Antipov D,
    4. Gurevich AA,
    5. Dvorkin M,
    6. Kulikov AS,
    7. Lesin VM,
    8. Nikolenko SI,
    9. Pham S,
    10. Prjibelski AD,
    11. Pyshkin AV,
    12. Sirotkin AV,
    13. Vyahhi N,
    14. Tesler G,
    15. Alekseyev MA,
    16. Pevzner PA
    . 2012. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol 19:455–477. doi:10.1089/cmb.2012.0021.
    OpenUrlCrossRefPubMed
  51. 51.↵
    1. Sahl JW,
    2. Lemmer D,
    3. Travis J,
    4. Schupp JM,
    5. Gillece JD,
    6. Aziz M,
    7. Driebe EM,
    8. Drees KP,
    9. Hicks ND,
    10. Williamson CHD,
    11. Hepp CM,
    12. Smith DE,
    13. Roe C,
    14. Engelthaler DM,
    15. Wagner DM,
    16. Keim P
    . 2016. NASP: an accurate, rapid method for the identification of SNPs in WGS datasets that supports flexible input and output formats. Microb Genom 2:e000074. doi:10.1099/mgen.0.000074.
    OpenUrlCrossRef
  52. 52.↵
    1. McKenna A,
    2. Hanna M,
    3. Banks E,
    4. Sivachenko A,
    5. Cibulskis K,
    6. Kernytsky A,
    7. Garimella K,
    8. Altshuler D,
    9. Gabriel S,
    10. Daly M,
    11. DePristo MA
    . 2010. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res 20:1297–1303. doi:10.1101/gr.107524.110.
    OpenUrlAbstract/FREE Full Text
  53. 53.↵
    1. Delcher AL,
    2. Phillippy A,
    3. Carlton J,
    4. Salzberg SL
    . 2002. Fast algorithms for large-scale genome alignment and comparison. Nucleic Acids Res 30:2478–2483. doi:10.1093/nar/30.11.2478.
    OpenUrlCrossRefPubMedWeb of Science
  54. 54.↵
    1. Nguyen LT,
    2. Schmidt HA,
    3. von Haeseler A,
    4. Minh BQ
    . 2015. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol 32:268–274. doi:10.1093/molbev/msu300.
    OpenUrlCrossRefPubMed
  55. 55.↵
    1. Schliep KP
    . 2011. phangorn: phylogenetic analysis in R. Bioinformatics 27:592–593. doi:10.1093/bioinformatics/btq706.
    OpenUrlCrossRefPubMedWeb of Science
  56. 56.↵
    1. Bruen TC,
    2. Philippe H,
    3. Bryant D
    . 2006. A simple and robust statistical test for detecting the presence of recombination. Genetics 172:2665–2681. doi:10.1534/genetics.105.048975.
    OpenUrlAbstract/FREE Full Text
  57. 57.↵
    1. Murray GGR,
    2. Wang F,
    3. Harrison EM,
    4. Paterson GK,
    5. Mather AE,
    6. Harris SR,
    7. Holmes MA,
    8. Rambaut A,
    9. Welch JJ
    . 2016. The effect of genetic structure on molecular dating and tests for temporal signal. Methods Ecol Evol 7:80–89. doi:10.1111/2041-210X.12466.
    OpenUrlCrossRef
  58. 58.↵
    1. Drummond AJ,
    2. Suchard MA,
    3. Xie D,
    4. Rambaut A
    . 2012. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol 29:1969–1973. doi:10.1093/molbev/mss075.
    OpenUrlCrossRefPubMedWeb of Science
  59. 59.↵
    1. Baele G,
    2. Lemey P,
    3. Bedford T,
    4. Rambaut A,
    5. Suchard MA,
    6. Alekseyenko AV
    . 2012. Improving the accuracy of demographic and molecular clock model comparison while accommodating phylogenetic uncertainty. Mol Biol Evol 29:2157–2167. doi:10.1093/molbev/mss084.
    OpenUrlCrossRefPubMedWeb of Science
  60. 60.↵
    1. Wyeth NC
    . 1927. Map of discovery of the western hemisphere. Wall map. National Geographic, Evergreen, CO.
View Abstract
PreviousNext
Back to top
Download PDF
Citation Tools
Dating the Cryptococcus gattii Dispersal to the North American Pacific Northwest
Chandler C. Roe, Jolene Bowers, Hanna Oltean, Emilio DeBess, Philippe J. Dufresne, Scott McBurney, David P. Overy, Bodo Wanke, Colleen Lysen, Tom Chiller, Wieland Meyer, George R. Thompson III, Shawn R. Lockhart, Crystal M. Hepp, David M. Engelthaler
mSphere Jan 2018, 3 (1) e00499-17; DOI: 10.1128/mSphere.00499-17

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print
Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this mSphere article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Dating the Cryptococcus gattii Dispersal to the North American Pacific Northwest
(Your Name) has forwarded a page to you from mSphere
(Your Name) thought you would be interested in this article in mSphere.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Dating the Cryptococcus gattii Dispersal to the North American Pacific Northwest
Chandler C. Roe, Jolene Bowers, Hanna Oltean, Emilio DeBess, Philippe J. Dufresne, Scott McBurney, David P. Overy, Bodo Wanke, Colleen Lysen, Tom Chiller, Wieland Meyer, George R. Thompson III, Shawn R. Lockhart, Crystal M. Hepp, David M. Engelthaler
mSphere Jan 2018, 3 (1) e00499-17; DOI: 10.1128/mSphere.00499-17
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

Cryptococcus
genomics
molecular epidemiology
mycology

Related Articles

Cited By...

About

  • About mSphere
  • Board of Editors
  • Policies
  • For Reviewers
  • For the Media
  • Embargo Policy
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Author Warranty
  • Types of Articles
  • Getting Started
  • Ethics
  • Contact Us

Follow #mSphereJ

@ASMicrobiology

       

 

Website feedback

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

 

American Society for Microbiology
1752 N St. NW
Washington, DC 20036
Phone: (202) 737-3600

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

Online ISSN: 2379-5042