miércoles, 18 de mayo de 2016

Use of Population Genetics to Assess the Ecology, Evolution, and Population Structure of Coccidioides - Volume 22, Number 6—June 2016 - Emerging Infectious Disease journal - CDC

Use of Population Genetics to Assess the Ecology, Evolution, and Population Structure of Coccidioides - Volume 22, Number 6—June 2016 - Emerging Infectious Disease journal - CDC

Volume 22, Number 6—June 2016


Use of Population Genetics to Assess the Ecology, Evolution, and Population Structure ofCoccidioides

Marcus M. Teixeira and Bridget M. BarkerComments to Author 
Author affiliations: Translational Genomics Research Institute, Flagstaff, Arizona, USA (M.M. Teixeira, B.M. Barker)Northern Arizona University, Flagstaff (B.M. Barker); University of Arizona, Tucson, Arizona, USA (B.M. Barker)


During the past 20 years, a general picture of the genetic diversity and population structure of Coccidioides, the causal agent of coccidioidomycosis (Valley fever), has emerged. The genus consists of 2 genetically diverse species, C. immitis and C. posadasii, each of which contains 1 or more distinct populations with limited gene flow. Genotypic data indicate that C. immitis is divided into 2 subpopulations (central and southern California populations) and C. posadasii is divided into 3 subpopulations (Arizona, Mexico, and Texas/South America populations). However, admixture within and among these populations and the current paucity of environmental isolates limit our understanding of the population genetics ofCoccidioides. We assessed population structure of Coccidioides in Arizona by analyzing 495 clinical and environmental isolates. Our findings confirm the population structure as previously described and indicate a finer scale population structure in Arizona. Environmental isolates appear to have higher genetic diversity than isolates from human patients.
Coccidioides immitis and C. posadasii are the only 2 species recognized within the genus Coccidioides (1). These fungi are endemic to arid or semi-arid regions of the Americas. Both species cause the disease coccidioidomycosis (Valley fever), which is contracted by dogs, humans, and other mammals living in or visitingCoccidioides-endemic areas (2,3). Infection is acquired through inhalation of air-dispersed arthroconidia (asexual single-cell fungal propagules). When a mammalian host inhales these conidia, a switch from polar to isotropic growth is initiated, resulting in the development of a specialized infectious structure called a spherule (4). Within 4 to 5 days, the mature spherules disrupt, releasing potentially hundreds of endospores, each of which are capable of developing into a new spherule (5). This cycle continues until the host’s immune system represses fungal propagation or the fungus goes quiescent (4). If infection is not controlled, it can disseminate to other organs and tissues and is capable of crossing the blood–brain barrier and causing meningitis, which is fatal if untreated (6). Approximately 40% of infections are symptomatic (4).
The geographic distribution of C. immitis was thought to be restricted to central and southern California (7). However, the range extends south into Baja California and east into Arizona, and recent work shows this species was also found in eastern Washington (8,9), at Dinosaur National Monument in Utah (10), and in a patient in Colombia with no travel history (11). The species C. posadasii is present in Arizona, with its range extending into Utah, Texas, and Mexico and dispersed populations in Central and South America (1215). C. immitis and C. posadasii probably co-occur in nature, given that both species have been isolated from patients in San Diego and Mexico and hybrid strains have been identified (1,16). Environmental sampling and recovery of isolates would be more helpful in confirming this hypothesis than using isolates derived from patients.
One approach to assessing genetic diversity in fungal populations is to develop microsatellite markers (17,18). Microsatellites are short (1–6 bp) tandem repeats, which are found throughout eukaryotic genomes and are thought to be evolving under neutrality in fungi (19). These markers have been useful in population genetics studies that compare genotypes among closely related fungal species or populations (17,2023). Here we focus on the genotyping ofCoccidioides strains from various origins by combining multiple studies in a meta-analysis and by using population genetics to clarify the causative agents of coccidioidomycosis.
Because coccidioidomycosis is increasing and disease severity is highly variable, defining genotypic distribution is important for monitoring outbreaks and determining whether increased pathogenicity is an emerging trait (24). Previous analysis showed that a single clone did not cause the rise in infection rates in Arizona; rather, each isolate recovered from a patient was unique (25). Thus, the question remains: why is coccidioidomycosis on the rise? It has been hypothesized that climate change, changes in human susceptibility, changes in reporting, or a result of the interaction of these factors, overlaid with high genetic variation and the possibility that Coccidioides can colonize new hosts and new environments, are some of the factors responsible (3,26,27). We aimed to answer 4 main questions: 1) if the subpopulation structure previously proposed has support when a larger dataset is analyzed by using multiple methods; 2) if there is evidence for population structure within Arizona; 3) if environmental isolates from Arizona are distinct from Arizona human host isolates; and 4) if patient data confound population structure because of incorrect identification of the point source of infection.


In total, we compiled data from 66 soil-derived isolates retrieved by mouse passage in Tucson, Arizona (28); 141 isolates from Arizona patients with Valley fever (25); 106 C. posadasii and 62 C. immitis isolates from a broad geographic range (1); and 266 clinical C. posadasii isolates (human and veterinary) newly analyzed for this study (Technical Appendix 1 Table). Of these 641 isolates, 22 were removed from final analysis for failure to amplify >2 of the 9 loci.
DNA Extraction
To extract DNA, we placed ≈0.2 g of mycelia in a 2-mL screw-cap tube containing 0.5-mm–diameter sterile glass beads (BioSpec, Bartlesville, OK, USA) and 1 ml of lysis buffer (50 mmol/L Tris-HCl [pH 7.5], 100 mmol/L EDTA [pH 8.0], 100 mmol/L NaCl, 0.5% sodium dodecyl sulfate, and 100 mmol/L β-mercaptoethanol) and subjected it to mechanical disruption by vortexing on a flat 12-tube holder (MoBio, Carlsbad, CA, USA) at 3,700 rpm for 10 minutes. Samples were incubated at 65°C for 60 minutes and centrifuged at 8,000 rpm for 5 minutes. We extracted nucleic acids from the supernatant with buffered phenol:chloroform:isoamyl alcohol pH 8.0 (25:24:1) and again with cholorform:isoamyl (24:1) and precipitated from the aqueous layer with 0.6 volumes of isopropyl alcohol. We washed the pellets twice with ethanol and resuspended them in 150 µL of double-distilled H2O. DNA concentration was determined on NanoDrop 1000 spectrophotometer (Thermo Scientific, Wilmington, DE, USA) and was diluted to 20 ng/µL.
Multilocus Microsatellite Typing Markers and PCR
To genotype isolates, we used 9 microsatellite primers developed for phylogenetic analysis and tested for concordance in Coccidioides (17,18,25). All microsatellite fragments were first denatured for 2 min at 96°C, followed by 30 amplification cycles (30 s at 94°C, 30 s at 55°C, and 1 min at 72°C) and 1 extension cycle of 5 min at 72°C with 2.5× Hotmaster mix (Eppendorf, New York, NY, USA). One primer from each set was end-labeled with a fluorescent tag (either NED dye [ABI, Shirley, NY, USA] or FAM or HEX [Eurogentec, Seraing, Belgium). Primer concentrations were 200 nmol/L each per reaction, and 100 ng of DNA was used for each reaction.
Fragment Analysis
We grouped microsatellite fragments from each isolate into 3 sets of 3 fragments and labeled 1 primer set in each grouping with HEX, FAM, or NED. Pooled PCR products were separated on an ABI 3730 DNA Analyzer (Applied Biosystems, Foster City, CA, USA) at the University of Arizona Genomic Analysis and Technology Core sequencing facility, using a ROX-labeled ladder (Invitrogen, Carlsbad, CA, USA) for sizing. Chromatographs were read in Genotyper (ABI, Shirley, NY, USA), and a single peak was scored (Coccidioides is haploid). No evidence of multiple peaks was detected. Microsatellites were amplified and analyzed at least twice to verify their size. To compare our isolates to those described in published data, we analyzed the microsatellite sizes from a subset of previously analyzed isolates on our ABI 3730 system (Technical Appendix 1 Table, duplicates tab). Calibration was necessary to compare the published microsatellite sizes to our data (Technical Appendix 1 Table, correction tab).
Population Analyses
We tabulate data from the Genotyper program maintained them in a spreadsheet (Technical Appendix 1 Table). Files were checked for duplicates and clone–correction checked using GenAlEx 6.501 (29). We found identical isolates from multiple isolations from the same patient and from isolates collected from the same soil site. Any samples that were missing >3 loci were eliminated from the final dataset. Locations were incorporated into a nexus file containing 619 isolates. We assigned locations based on the isolation/hospital origin as follows: Phoenix, Yuma, and Tucson (Arizona); San Diego and San Joaquin Valley (California); Texas; Mexico; and South America (Brazil, Argentina, and Venezuela).
We analyzed microsatellite matrices by using STRUCTURE 2.3.4 (Pritchard Laboratory, Stanford University, Stanford, CA, USA) to determine population structure within Coccidioides (30). The running length of burn-in period was 100,000 repetitions with 1,000,000 Markov chain Monte Carlo repetitions. Default settings in STRUCTURE 2.3.4 were as follows: the admixture model was used to infer α along with the previous sampling location information model (LOCPRIOR) (30). We used CLUMPP, a cluster matching and permutation program (https://web.stanford.edu/group/rosenberglab/clumpp.html), to define populations within the STRUCTURE algorithm. K is the number of significant populations in each main group. Allele frequencies were assumed to be correlated among populations, assuming that there are different Fst values for different subpopulations, the previous mean of Fst for populations is 0.01, and λ is constant at 1.0. Ten runs for each k from 1 to 10 were performed, and results were analyzed using Evanno’s method implemented in StructureHARVESTER (31). We generated a consensual STRUCTURE plot from the admixture values using the Clustering Markov Packager Across K (CLUMPAK) (http://www.clumpak.tau.ac.il) and built final plots with STRUCTURE PLOT (32,33).
We also inferred Coccidioides population splits and mixtures trees using a statistical model related to common ancestors through a graph of ancestral populations via TreeMix software (Pritchard Laboratory) (34). In brief, we inferred a population tree on the basis of microsatellite data for each of the identified populations in STRUCTURE (Technical Appendix 2 Table). Migration events were placed on admixed edges, which are correlated with the degree of ancestry for each population and represents unidirectional gene flow between populations. Horizontal branch lengths are proportional to the accumulated genetic drift (drift parameter) from each population that was placed in a given branch. The drift parameter measures the variance in allele frequency that changes along each population of the tree. We also analyzed the same data were by using Nei’s unbiased genetic distance estimate (Table 1), to complete a principal coordinate analysis (PCoA) (Table 2) in GENALEX 6.501 (http://www.biology-assets.anu.edu.au/GenAlEx/Welcome.html) (29). We documented allele frequencies, private alleles, and haploid diversity calculations (Table 3) for Arizona samples (Technical Appendix 2 Table).


Combining Data from Multiple Sources
Thumbnail of Results of principal coordinate analysis and STRUCTURE analyses of Coccidioides spp. populations. A) Principal coordinate analysis using Nei’s unbiased genetic distance estimates supports 3 main groupings: C. immitis, C. posadasii TX/SA/MX, and C. posadasii Arizona (see also Technical Appendix 3 Figure 2). The greatest separation occurs between species and is reflected in principal coordinate 1 (93.92% of variance). Color-coding for populations: lime green, San Joaquin Valley (SJV);
Figure 1. Results of principal coordinate analysis and STRUCTURE analyses ofCoccidioides spp. populations. A) Principal coordinate analysis using Nei’s unbiased genetic distance estimates supports 3 main groupings: C. immitis, C. posadasii TX/SA/MX,...
We documented microsatellite frequencies (Technical Appendix 3 Figure 1Technical Appendix 1 Table). Three loci (GAC2, 621.1, and ACJ) had low diversity in C. posadasii, and these same loci were variable in C. immitis. Three loci showed the opposite pattern (K01, K03, and K07) and had low diversity for C. immitis and are variable in C. posadasii. Three loci (K09, GA1, and GA37) were diverse for both species. These results were similar to those of earlier reports (1). We merged datasets were merged for analysis (Technical Appendix 1Table) and analyzed isolates from both published datasets (Technical Appendix 1, duplicates tab). Manual corrections of 1 or 2 bp were needed because of slight variations among machines and ladders (Technical Appendix 1, corrections tab).
Population Structure of Coccidioides Subspecies
STRUCTURE analysis based on 619 isolates revealed 3 Coccidioides populations for C. immitis (n = 61), C. posadasii Mexico/Texas/South America (n = 63), and C. posadasii Arizona (n = 495) (Technical Appendix 3Figure 2). We detected low gene flow between the 3 major populations as observed by unique bar plots for each of these populations and observed gene flow between C. posadasii Mexico/Texas/South America and C. posadasiiArizona and between C. immitis and C. posadasii Mexico/Texas/South America (Technical Appendix 3 Figure 2). The population tree displays the 3 main populations and population assignments for each isolate along the bar plots (Technical Appendix 3 Figure 2). PCoA analysis using Nei’s unbiased genetic distance estimates revealed 3 main groupings when considering all data (Figure 1, panel A). Principal component (PC) 1 explains 93.92% of the variation, mainly attributable to variation between the species (Eigen value 1.202). PC2 explains 3.95% of the variation, reflecting the subpopulation structure in both species (Eigen value 0.051). PC3 explains 1.44% of variation and further separates Mexico from Arizona (Table 2).
Population Structure within C. immitis Population
Thumbnail of Population tree of Coccidioides subspecies population splits and mixtures. Maximum-likelihood population tree and presence of gene flow between diverged Coccidioides populations were inferred by using TreeMix software and microsatellites data (34). Direction of arrow indicates migration or gene flow based on admixture models; migration weights are shaded according their importance, supporting gene flow from a soil-derived population (AZSOIL) recovered from animal passage to a clinic
Figure 2. Population tree ofCoccidioides subspecies population splits and mixtures. Maximum-likelihood population tree and presence of gene flow between divergedCoccidioides populations were inferred by using TreeMix software and microsatellites data (<>
Thumbnail of Coccidioides subspecies distribution for North, Central, and South America. The frequency of assignment for each Coccidioides population was plotted in a pie chart for each location, and numbers of isolates from each location are displayed. For example, patients from Mexico were infected with isolates from Texas, San Diego, and Mexico populations, as determined by analysis with STRUCTURE. Because each of the patients’ location is the hospital, no fine-scale population is assessed.
Figure 3Coccidioidessubspecies distribution for North, Central, and South America. The frequency of assignment for eachCoccidioides population was plotted in a pie chart for each location, and numbers of isolates from...
Results of PCoA analysis strongly indicated population structure within C. immitis, separating San Joaquin Valley (SJV) from San Diego and Mexico (SDMX) isolates (Figure 1, panel A). STRUCTURE analysis also indicates a strong population subdivision within C. immitis (Figure 1, panel B). According to the optimal number of clusters determined by using StructureHARVESTER, the SJV and SDMX isolates are clustered into 2 different populations (k = 2) (Figure 1, panel B; Technical Appendix 3 Figure 3). Bar plots show that limited gene flow was observed between these subpopulations; however, the bar plots also indicated that the C. immitis isolates 17TX, 4SD, 8SD, 4M3, and 8M3 share alleles from both populations (Technical Appendix 3 Figure 3). The population tree indicates that the C. immitis SJV population has a migration event from the C. immitis SDMX population (Figure 2). The isolate population distribution frequency for C. immitis reveals differences between SDMX and SJV populations (Figure 3).
Population Structure within C. posadasii Mexico/Texas/South America Population
For C. posadasii Mexico/Texas/South America population, we detected 2 optimal clusters (Technical Appendix 3 Figure 3), 1 including Texas/South America isolates and 1 constituting isolates from Mexico (Figure 1). The Mexico isolates display a high level of hybridization between 2 different populations as well as within C. immitis(Technical Appendix 3 Figures 1, 3). The migration event from C. immitis SDMX to C. posadasii Mexico and the more basal C. immitis to C. posadasii Tucson migration event implicate the Sonoran desert as a convergent source of multiple Coccidioides genotypes and possible center of origin of the genus (Figures 23).
Population Structure within C. posadasii arizona Population
Population structure analysis of 495 separate fungal isolates suggests at least 3 different C. posadasiisubpopulations in Arizona, in agreement with PCoA data (Figure 1). Clinical samples from Yuma and Phoenix (designated PHOENIX) and Tucson patients (designated TUCSON) fall in 2 different populations according to STRUCTURE (Figure 1, panel B). All environmental samples and some veterinary/clinical samples from Tucson, Phoenix, and Yuma regions constitute a third population (designated AZSOIL) apart from the TUCSON and PHOENIX populations (Figure 1Technical Appendix 3 Figure 3). We detected high level of admixture in the Arizona population, suggesting gene flow between 3 populations. However, the presence of private alleles for different loci within each of the 3 Arizona populations supports genetic isolation (Technical Appendix 2 Table, AZ_PAL tab). Structure plots of AZSOIL, PHOENIX, and TUCSON populations contain isolates with genotypes from all 3 populations (Technical Appendix 3 Figure 3). AZSOIL and TUCSON populations arose from the same geographic origin (Figure 3). The population tree (Figure 2) supports a migration event from AZSOIL to TUCSON. The AZSOIL is placed nearer to the ancestral branch for Arizona subpopulations. In addition, a low number of clinical isolates clustered with AZSOIL, leading us to consider variable pathogenicity or host specificity (Figure 2). We propose that a mammalian host or its close microenvironment (e.g., mammal burrows) could contribute to increased fitness of a virulent phenotype. Thus, the environmental reservoir could play a role in the emergence of pathogenic strains.
Clinical Isolates
Data obtained from genotyping human patient isolates might lead to incorrect estimates of population structure. Two C. immitis were found in patients in Phoenix hospitals, and both patients had confirmed travel to California; however, we analyzed only 1 because the other did not meet our cutoff criteria (25). A Texas patient isolate was determined to be C. immitis (1). Patients from China, Switzerland, and Colorado (1 patient from each) and 7 California patients were infected with C. posadasii (1). One of the widely used laboratory C. posadasii strains (Silveira) was isolated from a patient with coccidioidomycosis diagnosed in California. In northern Mexico (including Baja California) and southern Mexico (Michoacán state), many strains are genotyped as C. immitis but have evidence of hybridization with C. posadasii and signatures of introgression (16). Less is known about the prevalence of introgression found in the C. posadasii Mexico population. For Arizona isolates newly analyzed for this study, no C. immitis were identified (Technical Appendix 1 Table).


Multiple methods and previous reports show that there are 2 species within Coccidioides defined as C. posadasii and C. immitis (1,7). Within species, C. posadasii contains the 2 main populations of Texas/South America/Mexico and Arizona, and within C. immitis, 2 populations are suggested, SJV and SDMX, supported by our data and previous reports (Figure 1) (7). Gene flow between C. immitis populations is not abolished, as exemplified by the admixture isolates 17TX and 22SD (Technical Appendix 3 Figure 3). STRUCTURE analysis suggests that C. posadasii Arizona and Texas/Mexico/South America populations are highly differentiated, with few isolates sharing genotypes among them (Technical Appendix 3 Figure 3). Additionally, divergence between Mexico and South America/Texas is evident, such that they are evolving independently (Figure 1).
Within the Arizona population, we observed 3 clusters: PHOENIX, TUCSON, and AZSOIL (Figure 1). PHOENIX consistently groups separately from TUCSON and AZSOIL, which might reflect differences in ecology between Arizona upland (Tucson) and the Lower Colorado Valley (Phoenix and Yuma) or variation in pathogenicity among hosts. Variation in mean soil temperature, precipitation, natural hosts, and vegetation could exert differential selection pressure on the fungus in the environment (35,36). In addition, according to the population tree, the AZSOIL subpopulation appears to be basal within Arizona. The migration event from AZSOIL to TUCSON might reflect selection of more pathogenic genotypes because only ≈40% of infections are symptomatic (4), and even fewer of these would result in severe disease where the isolate would be collected from the patient (Figure 2). This leads us to propose that the AZSOIL subpopulation reflects greater diversity than the TUCSON and PHOENIX subpopulations and that this greater diversity might be driven by selection of certain pathogenic strains in humans.
Moreover, our soil sampling reflects diversity at only 7 locations in and around Tucson, and all samples were collected with a single year, whereas the patient isolates from Tucson were collected over a period of 30 years. These soil isolates were obtained using a highly sensitive murine model of coccidioidomycosis. Not all mice had evidence of illness, and infection was only realized upon necropsy. Thus, we might have selected for infectious strains, but we believe we captured diversity in pathogenesis. This assumption would suggest that we have underestimated diversity in the environment. Diversity at some soil locations was high (i.e., multiple genotypes were recovered), whereas other sites were clonal, or we only recovered a single colony. Thus, it was surprising to find higher unbiased genetic diversity in AZSOIL (0.607, ±0.078 SE) than in TUCSON (0.584, ±0.085 SE) (Table 3). Patient isolates can provide information on a coarse level, but finer-scale mapping of geographic and population boundaries will require environmental sampling and analysis of genotypes. Our data suggest that environmental isolates reflect a broad diversity of genotypes and only a subset may be capable of causing severe disease in humans. A primary concern with our analysis is the precise location of the isolate origin. Few environmental isolates of Coccidioides exist, and methods to obtain them for genotypic analysis are currently inadequate (9,10,28,36).
Admixtures were found in the Arizona population, and gene flow was observed between the three defined subpopulations (Technical Appendix 3Figure 3); however, the presence of private alleles within each of those subpopulations and high genetic distance supports genetic isolation (Tables 13). The same was observed for the Mexico isolates nested in the Texas/Mexico/South America population. Additionally, the presence of private alleles and high diversity within C. immitis suggests that our results are not affected by oversampling in the Arizona subpopulation (Table 3). Because a sexual life cycle has not been observed, questions related to frequency, timing, and directionality of genetic exchange remain to be explored experimentally. Additional multilocus microsatellite types might be needed to support populations or could be resolved by using whole genome sequence comparison. Our data support previous work identifying the same main populations (1,7) and can be further tested with additional single nucleotide polymorphisms identified using whole genome sequence comparison.
Questions remain about the population biology of Coccidioides. The spatial and temporal distribution of individual genotypes, the amount of spatial overlap between the 2 species, and population boundaries within each species are still unclear (9,10,16). Overlap between species is likely, because of the identification of both C. immitis and C. posadasii recovered among patients in San Diego and northern Mexico and the observation of hybridization and introgression (1,7,16). This work shows that analyzing a large number of patient isolates and assigning regional population information reveals the potential for population structure within Arizona, at a much finer scale than previously thought (Figure 1). Thus, genetic differences and population subdivision among isolates and populations are likely greater than has been shown to date. The question of which population is basal to the C. posadasiilineage remains unanswered, and greater efforts to explore genotypic variation in Texas, Mexico, and Central and South America are needed.
Understanding the ecology of Coccidioides has been a longstanding goal (3,4). We used multiple methods to understand population genetics and determine population structure (29,30). However, environmental isolates must be more deeply explored by using direct fungal isolation (not passaging in mice) or high-coverage metagenomic sequencing, so that a specific location can be assigned to each isolate and potential for greater genetic diversity in the environment could be specifically tested. Surveying human patient isolates will continue to be valuable to track new outbreaks, such as the current coccidioidomycosis cases in Washington State (9).
Investigating the ecology and distribution of genotypes within and among populations of a pathogen is important for monitoring outbreaks, determining variance in virulence, and predicting disease progression (37). Correlating disease severity with pathogen genotype by using genome-wide association studies might assist in identifying genetic-based differences in virulence (38). Monitoring disease progression and response to antifungal therapy in animal models of coccidioidomycosis with more than a few well-characterized laboratory strains might provide information that could assist with better treatment options (39). Finally, a better understanding of ecologic and environmental factors that influence the growth and reproduction of the organism will assist in predicting and preventing exposure to the pathogen (28).
Dr. Teixeira is a postdoctoral fellow for the Northern Arizona Center for Valley Fever Research in the Pathogen Genomics Division at the Translational Genomics Research Institute in Flagstaff, Arizona. His research interests are population genetics and genomics of human fungal pathogens, primarily fungi in the order Onygenales
Dr. Barker is director of the Northern Arizona Center for Valley Fever Research and assistant professor in the Pathogen Genomics Division at the Translational Genomics Research Institute in Flagstaff, Arizona. Her research interests are population genetics and genomics of Coccidioides, development of clinical assays for early detection of Coccidioides Valley fever, and understanding the interplay of host and pathogen factors associated with Coccidioides infection.


Thanks to the Valley Fever Center for Excellence, John Galgiani, M. Lourdes Lewis, and Marc Orbach for technical support and isolates. Thanks to Matthew Fisher and Kelsea Jewell for providing raw data and DNA.
This publication is based on data reported in B.M.B.’s dissertation work. B.M.B. was supported by the Graduate Interdisciplinary Program in Genetics at the University of Arizona. Deidentified information on physician and diagnostic laboratory–reported cases was provided by the Arizona Department of Health Services. The reanalysis of data and preparation of the manuscript was supported by internal startup funding from the Translational Genomics Institute and a National Institutes of Health/National Institute of Allergy and Infectious Diseases K22 grant awarded to B.M.B.


  1. Fisher MCKoenig GLWhite TJTaylor JWMolecular and phenotypic description of Coccidioides posadasii sp. nov., previously recognized as the non-California population of Coccidioides immitis. Mycologia2002;94:7384DOIPubMed
  2. Cairns LBlythe DKao APappagianis DKaufman LKobayashi JOutbreak of coccidioidomycosis in Washington state residents returning from Mexico. Clin Infect Dis2000;30:614DOIPubMed
  3. Nguyen CBarker BMHoover SNix DEAmpel NMFrelinger JARecent advances in our understanding of the environmental, epidemiological, immunological, and clinical dimensions of coccidioidomycosis. Clin Microbiol Rev2013;26:50525DOIPubMed
  4. Lewis ERBowers JRBarker BMDust devil: the life and times of the fungus that causes Valley fever. PLoS Pathog2015;11:e1004762DOIPubMed
  5. Huppert MSun SHHarrison JLMorphogenesis throughout saprobic and parasitic cycles of Coccidioides immitis. Mycopathologia.1982;78:10722DOIPubMed
  6. Chiller TMGalgiani JNStevens DACoccidioidomycosis. [viii. ]Infect Dis Clin North Am2003;17:4157DOIPubMed
  7. Fisher MCKoenig GLWhite TJSan-Blas GNegroni RAlvarez IGBiogeographic range expansion into South America by Coccidioides immitismirrors New World patterns of human migration. Proc Natl Acad Sci U S A2001;98:455862DOIPubMed
  8. Marsden-Haug NHill HLitvintseva APEngelthaler DMDriebe EMRoe CCCoccidioides immitis identified in soil outside of its known range—Washington, 2013. MMWR Morb Mortal Wkly Rep2014;63:450 .PubMed
  9. Litvintseva APMarsden-Haug NHurst SHill HGade LDriebe EMValley fever: finding new places for an old disease: Coccidioides immitis found in Washington State soil associated with recent human infection. Clin Infect Dis2015;60:e13DOIPubMed
  10. Johnson SMCarlson ELFisher FSPappagianis DDemonstration of Coccidioides immitis and Coccidioides posadasii DNA in soil samples collected from Dinosaur National Monument, Utah. Med Mycol2014;52:6107 and. DOIPubMed
  11. Canteros CEVélez HAToranzo AISuárez-Alvarez RTobón Jimenez AMPMolecular identification of Coccidioides immitis in formalin-fixed, paraffin-embedded (FFPE) tissues from a Colombian patient. Med Mycol2015;53:5207DOIPubMed
  12. Whiston ETaylor JWGenomics in Coccidioides: insights into evolution, ecology, and pathogenesis. Med Mycol2014;52:14955DOIPubMed
  13. Duarte-Escalante EZúñiga GFrías-De-León MGCanteros CCastañón-Olivares LRReyes-Montes MRAFLP analysis reveals high genetic diversity but low population structure in Coccidioides posadasii isolates from Mexico and Argentina. BMC Infect Dis2013;13:411DOIPubMed
  14. Brilhante RSde Lima RARibeiro JFde Camargo ZPCastelo-Branco DSGrangeiro TBGenetic diversity of Coccidioides posadasii from Brazil. Med Mycol2013;51:4327DOIPubMed
  15. Campins HCoccidioidomycosis in South America. A review of its epidemiology and geographic distribution. Mycopathol Mycol Appl.1970;41:2534DOIPubMed
  16. Neafsey DEBarker BMSharpton TJStajich JEPark DJWhiston EPopulation genomic sequencing of Coccidioides fungi reveals recent hybridization and transposon control. Genome Res2010;20:93846DOIPubMed
  17. Fisher MCKoenig GWhite TJTaylor JWA test for concordance between the multilocus genealogies of genes and microsatellites in the pathogenic fungus Coccidioides immitis. Mol Biol Evol2000;17:116474DOIPubMed
  18. Fisher MCWhite TJTaylor JWPrimers for genotyping single nucleotide polymorphisms and microsatellites in the pathogenic fungus Coccidioides immitis. Mol Ecol1999;8:10824DOIPubMed
  19. Lim SNotley-McRobb LLim MCarter DAA comparison of the nature and abundance of microsatellites in 14 fungal genomes. Fungal Genet Biol.2004;41:102536DOIPubMed
  20. Fisher MCDE Hoog S, Akom NV. A highly discriminatory multilocus microsatellite typing (MLMT) system for Penicillium marneffei. Mol Ecol Notes.2004;4:5158DOIPubMed
  21. Fisher MCAanensen Dde Hoog SVanittanakom NMultilocus microsatellite typing system for Penicillium marneffei reveals spatially structured populations. J Clin Microbiol2004;42:50659DOIPubMed
  22. Taylor MLHernández-García LEstrada-Bárcenas DSalas-Lizana RZancopé-Oliveira RMGarcía de la Cruz SGenetic diversity of Histoplasma capsulatum isolated from infected bats randomly captured in Mexico, Brazil, and Argentina, using the polymorphism of (GA)(n) microsatellite and its flanking regions. Fungal Biol. 2012;116:308–17.
  23. Matute DRSepulveda VEQuesada LMGoldman GHTaylor JWRestrepo AMicrosatellite analysis of three phylogenetic species ofParacoccidioides brasiliensis. J Clin Microbiol2006;44:21537DOIPubMed
  24. Barker BMJewell KAKroken SOrbach MJThe population biology of coccidioides: epidemiologic implications for disease outbreaks. Ann N Y Acad Sci2007;1111:14763DOIPubMed
  25. Jewell KCheshier RCage GDGenetic diversity among clinical Coccidioides spp. isolates in Arizona. Med Mycol2008;46:44955DOIPubMed
  26. Hector RFRutherford GWTsang CAErhart LMMcCotter OAnderson SMThe public health impact of coccidioidomycosis in Arizona and California. Int J Environ Res Public Health2011;8:115073DOIPubMed
  27. Fisher MCKoenig GLWhite TJTaylor JWPathogenic clones versus environmentally driven population increase: analysis of an epidemic of the human fungal pathogen Coccidioides immitis. J Clin Microbiol2000;38:80713 .PubMed
  28. Barker BMTabor JAShubitz LFPerrill ROrbach MJDetection and phylogenetic analysis of Coccidioides posadasii in Arizona soil samples. Fungal Ecol2012;5:16376DOI
  29. Peakall RSmouse PEGenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research—an update. Bioinformatics.2012;28:25379DOIPubMed
  30. Hubisz MJFalush DStephens MPritchard JK. Inferring weak population structure with the assistance of sample group information. Mol Ecol Resour. 2009;9:1322–32.
  31. Earl DvonHoldt B. STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour. 2012;4:359–61.
  32. Kopelman NMMayzel JJakobsson MRosenberg NAMayrose I. Clumpak: a program for identifying clustering modes and packaging population structure inferences across K. Mol Ecol Resour. 2015;15:1179–91.
  33. Ramasamy RKRamasamy SBindroo BBNaik VG. STRUCTURE PLOT: a program for drawing elegant STRUCTURE bar plots in user friendly interface. Springerplus. 2014;3:431.
  34. Pickrell JKPritchard JKInference of population splits and mixtures from genome-wide allele frequency data. PLoS Genet2012;8:e1002967.DOIPubMed
  35. Baptista-Rosas RCHinojosa ARiquelme MEcological niche modeling of Coccidioides spp. in western North American deserts. Ann N Y Acad Sci.2007;1111:3546DOIPubMed
  36. Fisher FSBultman MWJohnson SMPappagianis DZaborsky ECoccidioides niches and habitat parameters in the southwestern United States: a matter of scale. Ann N Y Acad Sci2007;1111:4772DOIPubMed
  37. Litvintseva APBrandt MEMody RKLockhart SRInvestigating fungal outbreaks in the 21st century. PLoS Pathog2015;11:e1004804DOIPubMed
  38. Muller LALucas JEGeorgianna DRMcCusker JHGenome-wide association analysis of clinical vs. nonclinical origin provides insights intoSaccharomyces cerevisiae pathogenesis. Mol Ecol2011;20:408597DOIPubMed
  39. Thompson GR IIIStevens DAClemons KVFierer JJohnson RHSykes JCall for a California coccidioidomycosis consortium to face the top ten challenges posed by a recalcitrant regional disease. Mycopathologia2015;179:19DOIPubMed



Technical Appendices

Suggested citation for this article: Teixeira MM, Parker BM. Use of population genetics to assess the ecology, evolution, and population structure ofCoccidioides, Arizona, USA. Emerg Infect Dis. 2016 Jun [date cited]. http://dx.doi.org/10.3201/eid2206.151565

DOI: 10.3201/eid2206.151565

No hay comentarios:

Publicar un comentario