Clonality Despite Sex: The Evolution of Host-Associated Sexual Neighborhoods in the Pathogenic Fungus


Molecular genetic approaches typically detect recombination in microbes regardless of assumed asexuality. However, genetic data have shown the AIDS-associated pathogen Penicillium marneffei to have extensive spatial genetic structure at local and regional scales, and although there has been some genetic evidence that a sexual cycle is possible, this haploid fungus is thought to be genetically, as well as morphologically, asexual in nature because of its highly clonal population structure. Here we use comparative genomics, experimental mixed-genotype infections, and population genetic data to elucidate the role of recombination in natural populations of P. marneffei. Genome wide comparisons reveal that all the genes required for meiosis are present in P. marneffei, mating type genes are arranged in a similar manner to that found in other heterothallic fungi, and there is evidence of a putatively meiosis-specific mutational process. Experiments suggest that recombination between isolates of compatible mating types may occur during mammal infection. Population genetic data from 34 isolates from bamboo rats in India, Thailand and Vietnam, and 273 isolates from humans in China, India, Thailand, and Vietnam show that recombination is most likely to occur across spatially and genetically limited distances in natural populations resulting in highly clonal population structure yet sexually reproducing populations. Predicted distributions of three different spatial genetic clusters within P. marneffei overlap with three different bamboo rat host distributions suggesting that recombination within hosts may act to maintain population barriers within P. marneffei.


Published in the journal: . PLoS Pathog 8(10): e32767. doi:10.1371/journal.ppat.1002851
Category: Research Article
doi: https://doi.org/10.1371/journal.ppat.1002851

Summary

Molecular genetic approaches typically detect recombination in microbes regardless of assumed asexuality. However, genetic data have shown the AIDS-associated pathogen Penicillium marneffei to have extensive spatial genetic structure at local and regional scales, and although there has been some genetic evidence that a sexual cycle is possible, this haploid fungus is thought to be genetically, as well as morphologically, asexual in nature because of its highly clonal population structure. Here we use comparative genomics, experimental mixed-genotype infections, and population genetic data to elucidate the role of recombination in natural populations of P. marneffei. Genome wide comparisons reveal that all the genes required for meiosis are present in P. marneffei, mating type genes are arranged in a similar manner to that found in other heterothallic fungi, and there is evidence of a putatively meiosis-specific mutational process. Experiments suggest that recombination between isolates of compatible mating types may occur during mammal infection. Population genetic data from 34 isolates from bamboo rats in India, Thailand and Vietnam, and 273 isolates from humans in China, India, Thailand, and Vietnam show that recombination is most likely to occur across spatially and genetically limited distances in natural populations resulting in highly clonal population structure yet sexually reproducing populations. Predicted distributions of three different spatial genetic clusters within P. marneffei overlap with three different bamboo rat host distributions suggesting that recombination within hosts may act to maintain population barriers within P. marneffei.

Introduction

Hypotheses of globally continuous populations and strict clonality in putatively asexual microbial pathogens are rarely supported [1], [2], [3], [4], [5]. Instead, genetic approaches detect recombination in microbes regardless of assumed asexuality, and pathogens are surprisingly promiscuous despite strong population genetic structure [6], [7], [8]. In some eukaryotic pathogens spatial structuring is readily attributable to dispersal limitations [9], [10], but many fungal pathogens of humans display extensive spatial population genetic structure despite their ability to disperse via aerosolized spores [11], [12], [13], [14], [15]. Examples that cause extensive morbidity include Cryptococcus neoformans, Coccidioides sp., Histoplasma capsulatum and Penicillium marneffei. These fungi are maintained in natural environmental reservoirs that might contribute to structured populations via local adaptation, and they are thought to be largely clonal. However, fungi have many different mating systems that encompass asexual propagation t7hrough multiple forms of sexual and parasexual recombination, and clonal structure may be arrived at via very different mechanisms [16], [17]. Evidence suggests that population structure in fungal pathogens is strongly influenced by host distributions and extrinsic geographic boundaries [12], [18], [19], [20], [21]. Therefore, the interplay between mating systems, population structure, and host adaptation is a central question underpinning the evolutionary epidemiology of fungal pathogens.

Heterothallic mating systems in fungi require physical contact between two isolates containing opposite mating types at the mating-type locus (MAT) in order to undergo sexual reproduction. If no mating partners are present, then sexual reproduction does not occur and the fungus reproduces asexually (but see Lin et al. [22] for evidence of same-sex mating in the otherwise heterothallic fungus Cryptococcus neoformans). In this case, the relative capacities of fungal lineages to disperse and co-occupy environmental niches can drive population-level recombination rates. If strains of opposite mating type do not equally penetrate environments then species recombination rates may be reduced to levels nearing complete asexuality [23]. Previously, it has been shown that the HIV-associated emerging pathogen Penicillium marneffei shows extensive spatial genetic structure at local and regional scales across Thailand [XPATH ERROR: unknown variable "start2".], [24]. Although there has been some genetic evidence that a sexual cycle is possible in P. marneffei, this haploid fungus is thought to be genetically, as well as morphologically, asexual within these populations [25].

In this study, we use comparative genomics, experimental approaches, and population genetic data to identify the role of sexual recombination in maintaining spatial and genetic structure in this infection. We attempt to answer 4 specific questions: 1) Does the P. marneffei genome show evidence of sex? 2) How are populations of P. marneffei genetically structured? 3) Can population structure be reconciled with sex? 4) Do spatial or host factors correlate with population structure and sex? We use comparative genomics to identify genes linked to mating and genomic signatures of mutation bias associated with meiosis, and we experimentally detect recombination in vivo. We expand our collection of population genetic data across southeast Asia to include mating type data and 34 isolates from bamboo rats in India, Thailand and Vietnam, and 273 isolates from humans in China, India, Thailand, and Vietnam. Together these data form a mosaic that reveals some physical and genetic underpinnings of mating in P. marneffei that are linked to patterns of genetic diversity across its known endemic range.

Materials and Methods

Genomic analyses

We used 84 sexual cycle genes (Table S1) to blast against the NCBI genome sequences NZ_AAHF00000000 (A. fumigatus), NZ_ABAR00000000 (P. marneffei), and NZ_ABAS00000000 (T. stipitatus). We screened transposon families for substitution bias by making BLAST based alignments to determine the dominant form of a functional integrase gene in each family. We counted the type of substitution based on differences of alleles as low as 70% identical to the dominant intact type. We compared gene sequences of the genomic region between slaB and apnB (the genes that flank the MAT idiomorph in related fungi) of strains FRR2161 and FRR3842.

Isolates, MLMT barcoding and MAT discrimination

We acquired 307 isolates of P. marneffei from humans and vertebrate hosts (bamboo rats), covering the known global range of the fungus. Our study obtained 273 epidemiologically unlinked human isolates of P. marneffei from HIV-AIDS patients covering the time-period 1959 to 2005. Of these isolates, 258 were georeferenced to either the broad geographical region of collection or the patients home address. The remaining 15 isolates were recovered from patients whose infections were diagnosed in non-endemic regions, and no accurate geographical origin could be assigned. In addition to human isolates of P. marneffei, we obtained 34 isolates from the bamboo rats species Rhizomys pruinosis (n = 3), R. sumatrensis (n = 13), R. sinensis (n = 1) and Cannomys badius (n = 17). We also include the type isolate for P. marneffei Segretain et al. ATCC 18224, CBS 388.87, isolated from R. sinensis in 1959 [26]. All isolates were cultured on Sabouraud's agar and DNA extracted as previously described [27]. Subsequently, isolates were genotyped at 21 microsatellite loci using the methods described in Fisher et al. [11], [27].

The presence within each of isolate of the MAT1-1 α box and MAT1-2 high mobility group idiomorphs was determined using the PCR protocol detailed by Woo et al. [25].

Population genetic structure and spatial analyses

Genotypes were analysed using GenAlex 6.0 [28] to determine allelic diversity, genotypic diversity and spatial correlation across regions and the global distribution of P. marneffei. We used the package adegenet and its dependencies in R to conduct spatial PCA and DAPC analyses [29]. To compare our inferred results against a model of a single continuous population structured by a dingle dispersal kernel and mutation rate we used the coalescent based program IBDsim [30]. Additional distribution data for bamboo rats were collected from specimen databases AMNH, FMNH, NMNH, and the GBIF. We used the bioclim layers 1–21 at 30 sec from the world clim database in MAXENT to generate predicted distributions for bamboo rats and P. marneffei genetic clusters. We measured distributional overlap using Schoener's D and a resampling approach [31]. We compared the relative overlap of genetic clusters to host distributions by generating null distributions of D based on resampling of R. sumatrensis and Cannomys (Text S1). Possible parental distances were compared to null distributions generated by choosing isolates randomly that met the genetic distance criteria from the population that met the parental criteria.

Co-infection analysis

Five co-housed outbred CD-1 male mice (16–18 g) were inoculated intranasally with 107 spores suspended in 40 µl of PBS. Conidia from two isolates, PM9, a MAT 1–2 isolate from Thailand, and the type strain ATCC18824 (FR2161) were mixed in a 1∶1 ratio to form the inoculum (S9). Serial dilutions of homogenized saline samples were plated (no later than 6 hours after they were removed from the mice) on Sabouraud agar. Colonies were counted after 4 days in 27°C. Individual colonies used for DNA extraction and subsequent genotyping as before [11]. Isolate genotypes were compared to the initial genotypes of the inoculum and genotypes differing from inoculum were confirmed via DNA sequencing.

Ethics statement

All the clinical studies from which isolates are available were approved by the Wellcome Trust ethics committees at the study sites, in the UK and by the regulatory authorities of the countries involved. All patients or their next of kin gave written informed consent and all patient data are anonymised. This work strictly complied with the animal regulations and guidelines under UK law and was approved by Imperial College's Ethical Review Process (ERP) Committee and the British Home Office. All murine work was carried out in a Biosafety level 3 secure animal facility under licensed approval from the British Home Office.

Results/Discussion

Genomic evidence for sex in P. marneffei

Sexual reproduction leaves an imprint on fungal genomes by maintaining genes required for mating and by generating patterns of mutation and recombination restricted to meiotic processes [32], [33], [34], [35], [36]. Successful mating in fungi requires that a genome contains a functioning series of interconnected genetic pathways [37]. Using a comparative genomic approach we assessed the presence and functionality in P. marneffei of genes known to be involved in sexual development in fungi. First, comparing between strains FRR2161 and FRR3842 revealed that the region between genes slaB and apnB resembled other fungal mating type idiomorphs. A region of complete dissimilarity was flanked by regions that were nearly identical between the strains (Fig. 1A and B). We found homologs for nearly all of the genes needed for a complete sexual cycle in yeast to be present and putatively functional (Table S1). Those genes not detected in P. marneffei were also not detected in Talaromyces stipitatus, a fungus with a complete sexual cycle, and most genes that were absent in those two fungi were also missing in the recently demonstrated heterothallic fungus Aspergillus fumigatus. Although these sex-related genes may be conserved to function in processes other than mating, their presence suggests that P. marneffei has preserved the ability to complete a sexual cycle. We detected another genomic signature of a functional sexual cycle, a type of mutation bias associated with meiosis. Repeat induced point mutation (RIP), a process by which some fungi silence genes involved in mobile genetic element function by preferentially mutating repeated sequences within their genomes, is associated with meiosis [38]. This process results in skewed base pair distributions due to the induced mutations. Using an approach similar to that of Clutterbuck [39], we found evidence of an excess of sliding windows with zero AG and CT dinucleotides and mutation bias in P. marneffei transposon family Ty-1 with a skew towards G to A and C to T transitions (Fig. 2). We also detected a putatively functional RID gene (Locus ID PMAA079888), the only conserved gene so far implicated in RIP [40], [41], [42]. Although the RID gene and the observed mutation bias can be explained by several factors including those acting during mitosis, they point towards a RIP or RIP-like process that is generally considered a feature of sexually reproducing fungi and an overall genomic pattern consistent with sex. Although these genomic signatures could represent relics from a sexual past rather than ongoing sexual recombination within P. marneffei, in the related human pathogenic fungus Aspergillus fumigatus, the discovery of mating type genes and evidence of RIP heralded the eventual description of a full sexual life cycle [35], [43], [44].

Fig. 1. Mating type idiomorph of P. marneffei.
Mating type idiomorph of <i>P. marneffei</i>.
A) A dot plot showing similarity between strains FRR2161 and FRR3482 similarity across the 7 kb covering MAT loci. B) Gene cartoon showing arrangement of genes in the MAT idiomorphs. Grey regions are areas of high homology, and black lines are small segments of rearrangement. Coding regions of MAT alpha and HMG genes are shown in black, flanking genes with high homology are shown in grey, and hypothetical genes with low homology are shown in white. Arrows indicate the direction of transcription.

Fig. 2. Biased mutation in P. marneffei.
Biased mutation in <i>P. marneffei</i>.
Nucleotide bias in P. marneffei. A) AG dinucleotide frequency across the entire genome of P. marneffei shows a non-normal distribution consistent with patterns of RIP-like process B) Comparison between the types of mutations from the inferred ancestral state shows a bias in transitions and towards G to A and C to T mutations in particular.

Spatially correlated genetic diversity

Microsatellite allelic diversity was high overall and within localities (Table 1). With the exception of Thai Central and Thai South, populations assigned a priori by locality were significantly differentiated from one another by Wright's FST [45]; this metric ascertains the proportion of genetic variance among geographical regions relative to the total variance. FST values near zero mean that populations are not distinct and variation is shared equally within and between them, while higher FST values mean that more genetic differences occur between populations compared to those within populations. The China and Taiwan populations were most different from the other a priori populations (Table 2). Phylogenetic analysis revealed associations between sampling area and the occurrence of phylogenetic clustering (Fig. 3). Using discriminant analysis of principal components (DAPC) to identify genetic clusters [46], [47], we assigned individuals to 3 clusters based on the Bayesian information criterion. The clusters show some spatial association. Cluster 1 is composed mostly of isolates from central and southern Thailand, Cluster 2 of isolates from China, and Cluster 3 of isolates from northern Thailand (Fig. 3, Figure S1). As expected, we observed a strong pattern of spatial genetic correlation (r2 = 0.41, p<0.01). We also detected significant ‘global’ genetic structure (positive correlation between spatial and genetic distance) but no ‘local’ genetic structure (negative correlation) using spatial principal coordinate analysis [46], [47]. To test for a homogenous neutral process of genetic differentiation we used a spatially explicit coalescent-based simulation of isolation by distance generated with IBDsim [30] to simulate a uniform dispersal/mutation process across our exact sampling scheme. This uniform genetic structure was then compare against our recovered spatial genetic pattern. By controlling for the spatial distribution of our sample sites we are able to determine if the apparent genetic clustering is simply an artifactual product of clustered sampling and a single uniform process of genetic clustering. Because the hypothesised parameter space is nearly infinite, we concentrated on dispersal scenarios that most closely resembled the spatial genetic correlation present in our data, namely, the strength of spatial genetic correlation at the smallest spatial scale and the decay rate of the correlation. Although the simulated datasets largely overlapped with our recovered data we observed important departures between the two. The simulated datasets had a single peak in spatial genetic correlation at the smallest spatial scale and a decay in correlation dependent on the dispersal kernel, a feature of all single population isolation-by-distance models, but the observed pattern had additional peaks in certain distance classes that disrupted the uniform decay (Fig S2 and S3). One peak was composed of distances between individuals belonging to the outer edges of clusters 1 and 3. Another major peak comes at the spatial scale where the outer edges of Clusters 1 and 3 contact with Cluster 2. These results differ from previous results that observed different rates of decay for spatial genetic correlation [24], a feature that probably owes to the limited geographic scope and power of the earlier study. Our data now suggest that the observed genetic clusters are not the result of a process of uniform decay with geographic distance, and that other factors are also driving the heterogeneity observed in our dataset.

Fig. 3. Correspondence between genetic diversity and spatial location.
Correspondence between genetic diversity and spatial location.
A) Unrooted neighborjoining tree of genetic distances between isolates labelled with the colour each isolation locality. B) Genetic clusters inferred using DAPC and isolates coloured according to cluster. C) Spatial distribution of DAPC clusters.

Tab. 1. Local population diversity.
Local population diversity.
n, number of samples; MT, number of microsatellite types ( = haplotypes); Na, mean number of alleles; Ne, mean effective number of alleles ; He, mean haploid genetic diversity; Pa, mean number of private alleles. Standard errors are shown in parentheses for Na, Ne, He and Pa.

Tab. 2. Local population differentiation.
Local population differentiation.
FST, the proportion of variance among geographical regions relative to the total variance. Pairwise population values of FST are below the diagonal with significant values (p<0.001 from 999 permutations) in bold. Pairwise Nei's genetic distances are shown above the diagonal.

Recombination and linkage disequilibrium

Linkage disequilibrium was high throughout the sample with an overall of 0.113 (Table S2). We determined the relative frequency of mutation to recombination using the single locus variant approach applied by Fisher et al. [24] and found a mutation to recombination frequency of 0.083 suggesting that mutation is up to 12 times less frequent than recombination across the whole population. When restricted to only bamboo rat isolates, all single locus variants would be due to recombination, while for human-only isolates the ratios are unchanged in comparison to the entire dataset. Average fungal microsatellite mutation rates have been inferred from between 2.80×10−6 and 2.50×10−5 mutations per generation [48], making the inferred recombination rate in P. marneffei between 3×10−5 and 2.50×10−4, a rate about half that observed in wild yeast [49]. This approach only detects single locus recombination events, which may be a minority in eukaryotic populations, while it should detect virtually all mutations that have not otherwise been masked by recombination. However, the method could be strongly biased towards inferring recombination due to convergent mutations in microsatellite length. The measure of minimal recombination (RM), which represents the minimum number of recombination events necessary to explain alleles failing the four gamete test [50] given the arrangement of the alleles in a contig, showed that recombination did occur within contigs (Figure S4). Complete clonality and complete panmixia are rejected for P. marneffei, but similar to previous results the inferred levels of clonality remain among the highest observed for fungi [24]. To explain the high level of clonal structure either recombination must be rare or it must occur largely between closely related individuals.

Distribution of MAT genes

The entire sample population of P. marneffei showed a distribution of mating types that was significantly skewed (p = 0.02 or p = 0.04 when clone corrected) from a ratio of 1∶1 in favour of an overabundance of MAT1-1 alleles, but some local populations were skewed towards MAT1-2 alleles (Table 1). Two of the genetic clusters inferred by DAPC were skewed towards more MAT1-1 alleles, but MAT genes within the central cluster did not differ from a 1∶1 ratio (Table 3). As predicted in work prior to the discovery of MAT loci in P. marneffei, highly skewed MAT ratios would be expected in a predominately asexual population [24]. On one hand, in the absence of sex and selection, MAT genes at an initial frequency of 0.5 are expected to be fixed in a population on average by ln 2(Ne) generations. Alternatively, in a completely sexual population without selection associated with a mating type, MAT alleles would be maintained at frequencies near 0.5 with very limited variance because all individuals in each generation will possess MAT alleles according to a binomial distribution, and there is no opportunity for drift beyond a single generation. MAT allele counts can be used to represent the reduction in effective population size caused by drift in MAT ratios [51], [52], but this assumes a fully sexual population. When sex is limited, the average allele frequencies for populations that do not lose sex and become fixed remain 0.5, but the variance in MAT allele frequency depends on population size and the frequency of sex. Based only on the differences in MAT allele frequencies between clusters and an intrinsic restriction on sex, P. marneffei would have an intrinsic upper bound of sexual recombination frequency at less than 4.5% given a modest population size of 1000. This small level of sex could explain the highly skewed ratio of MAT alleles in Cluster 3 and still accommodate the 1∶1 ratio in Cluster 1 while avoiding any fixation of MAT alleles. However, if the intrinsic sexual recombination rate explained the distribution of MAT alleles it would predict equal frequencies of clone detection across populations. Instead, percent clonality tracks the MAT allele skew, suggesting that sexual recombination in Cluster 3 is reduced relative to Cluster 1 (Table 3). We do note, however, that MAT allele frequencies would not be informative about where sex occurs if unisexual mating occurs in P. marneffei as is known in C. neoformans [22].

Tab. 3. Mating type counts within clusters and the percent of isolates with MLMT identical to at least one other isolate in the sample.
Mating type counts within clusters and the percent of isolates with MLMT identical to at least one other isolate in the sample.
Ne(MAT) is the estimated effective population size based on mating type frequency within the cluster.

Sexual neighbourhoods

Given that recombination occurs in P. marneffei, we wanted to determine the geographic and genetic scope of sex. Out of 43 clonal groups inferred with EBURST [53], six contained both MAT alleles, and four otherwise genetically identical multi-locus microsatellite types contained both MAT alleles (Fig. 4). Otherwise genetically identical isolates that differ only at mating type have also been detected in Cryptococcus gattii populations [54], [55], [56]. Clones with both MAT alleles represent the smallest possible genetic scale of sex, and are unequivocal evidence for recombination. To detect the more divergent recombination events we defined putative recombinants as any genotype that had no unique alleles, yet differed from the most similar genotype for at least three loci. Putative parents or ancestral parents were defined as all isolates that together could complete the multi-locus microsatellite type (MLMT) of the recombinant genotype (Figure S5). This allows us to compare between observed distances of maximal observed recombination against a null hypothesis that any two isolates could recombine. We identified 11 potential recombinants with 54 possible parent genotype combinations. Geographic distance between putative parents was shorter, 382 km, and genetic similarity higher, 60.06% identical, than random potential parents drawn from the entire population, 675 km (p = 0.005) and 49.29% identical (p = 0.025) respectively.

Fig. 4. EBURST diagram showing clones with mating types.
EBURST diagram showing clones with mating types.
Mating types are shown for each clone composed of more than a single isolate. MAT1-1-1 isolates are labelled with a red 1, and MAT1-2-1 isolates are labelled with a blue 2. The four genotypes that have both mating types are labelled by Microsatellite Type (MT) number.

The scale of recombination determines the efficacy of adaptation and the adaptive potential of populations. Although recombination across large distances allows generation of greater genetic diversity and more rapid spread of advantageous alleles, it disrupts locally advantageous combinations reducing local ecological genetic correlation. When sex is limited to small geographic distances it can reinforce local adaptation, and when limited to smaller genetic distances can reinforce genomic coadaptation. Together these effects can promote ecological speciation [57], [58], [59]. When ecological adaptation acts to reinforce genetic differentiation, strong correlations between key ecological factors and population distributions will exist [60], [61], [62], [63], [64]. To assess the possibility that ecological adaptation drives population differentiation in P. marneffei we used MAXENT [64] to predict overlap between the ecological niches of the genetic clusters (Fig. 5). Cluster 1, the cluster with the ratio of MAT alleles nearest 1∶1, had the widest predicted range and overlapped with the entire predicted range of Cluster 3, including the predicted range that was not sampled in Myanmar. Cluster 1, Cluster 2, and Cluster 3 isolates are all found in bamboo rats, but 16 of 17 samples from Cannomys badius and 13 of 14 unique genotypes were from cluster 1 and distributed among India, Thai Central, Thai North, and Thai South sampling localities. None of the Cluster 1 isolates were among the 13 recovered from Rhizomys sumatrensis, which were all in Cluster 3. Cannomys badius is relatively more abundant in the western portion of the range of P. marneffei. The predicted distributions of bamboo rats were similar to the IUCN species ranges and had overlap with P. marneffei distribution. Although our spatial sample of Cluster 2 was geographically restricted it was entirely within the distribution of Rhizomys sinensis, a species that has been shown to consistently harbour clinically relevant P. marneffei [65]. The distributions for R. sumatrensis overlapped with Cannomys and Clusters 1 and 3 (Fig. 5). However, using ENMTools [66] to account for sampling error we found that Cluster 1 predicted distributions overlapped more with the Cannomys distribution than R. sumatrensis distribution, and Cluster 3 similarly overlapped more with R. sumatrensis than Cannomys distribution (S8). This observed range overlap supports a host specific effect on P. marneffei population structure.

Fig. 5. Predicted distributions of P. marneffei genetic clusters and bamboo rat species.
Predicted distributions of <i>P. marneffei</i> genetic clusters and bamboo rat species.
A) Niche model distribution of Cluster 1. B) Niche model distribution of Cluster 2. C) Niche model distribution of Cluster 3. D) IUCN distribution of Cannomys badius. E) IUCN distribution of Rhizomys sinensis. F) IUCN distribution of Rhizomys sumatrensis. G) Niche model distribution of C. badius. H) Niche modeldistribution of R. sinensis. I) Predicted distribution of R. sumatrensis.

Hosts may structure populations of pathogenic fungi in many ways, including by providing an environment in which recombination can occur and by acting as a selective filter on population genetic diversity [67], [68]. We used a murine inhalation model of co-infection with genetically distinct strains to investigate the effect of host infection on P. marneffei (Text S2). Isolates of different mating types were used for experimental co-infection of 5 mice. Subsequent culture after 15 days from the livers showed a strong bias towards recovery of the MAT1-1 genotype for each of the mice. However, in two mice, genotypes of 4 isolates recovered from co-infections also revealed infrequent transfer of alleles between isolates of different mating type and genetic cluster (Table 4), suggesting that recombination may be possible across genetic barriers if multiple strains are within a host. In a smaller but similar in vitro experiment we did not observe significant bias towards MAT1-1, and from our scan of partial genotypes we did not recover any recombinants (S9, Table 4). We do not rule out regular recombination outside of hosts, but in the context of our spatial genetic evidence, the result of experimental infections indicate that hosts may play an important role in the development of sexual neighborhoods in populations of P. marneffei. However, the evolution of that role may involve restricted mating with or without host adaptation and remains to be explored.

Tab. 4. Counts of isolate genotypes recovered from murine co-infection and <i>in vitro</i> co-inoculation with ATCC18224 and PM9.
Counts of isolate genotypes recovered from murine co-infection and &lt;i&gt;in vitro&lt;/i&gt; co-inoculation with ATCC18224 and PM9.

Causes and consequences of sexually driven clonal structure

Asexual spores are common in vitro and likely a feature of natural P. marneffei populations, but sexual recombination may be an unexpectedly common occurrence in natural populations. The evidence supports the occurrence of recombination and perhaps even frequent sex, yet the natural populations remain strongly clonal and spatially structured. Although many mycologists might perceive this as a paradox because clonality is usually used as a proxy for asexuality, many fungi, including key pathogens, also employ same clone mating or sibling mating [7], [16], [69]. Three key hypotheses could explain the perceived clonality in P. marneffei; 1) Spatially restricted dispersal keeps individuals in contact with only closely related individuals; 2) Genetic incompatibility between dissimilar individuals restricts sex to genetically similar individuals; 3) Local adaptation restricts the ability of dissimilar genotypes to penetrate habitats ensuring mating between genetically similar individuals. All three are likely to be partially correct. Although the genetic evidence shows spatial limitations to effective dispersal, the physical dispersal of airborne conidia is not likely to be a limiting factor, and four genetically identical clones are dispersed across distances over 800 km. We have little information about the effect of genetic similarity on mating success in P. marneffei, but genetic restrictions on successful recombination are present in some plant pathogens [63], [70], [71] and should not be completely discounted. Local adaptation is not fully supported by ecological niche models that show overlap between distinct genetic clusters, but there is limited evidence of host specialisation. A key question unanswered in all of these hypotheses is why have sex at all?

Previous work has focused on the consequences of selectively neutral loss of sex in P. marneffei [72], but the persistence of a sexual cycle in P. marneffei despite abundant asexual reproduction in the lab suggests that there is a selective advantage for sex not associated with the advantages of greater adaptive potential provided by outcrossing. One major consequence of sexual clonality is release from Muller's Ratchet compared to asexuality [51]. Large numbers of haploid offspring and wide dispersal maximize environmental exposure of genets and increase the efficiency of purging deleterious alleles [73], [74]. However, the sexual process itself may also reduce the accumulation of deleterious mutations independent of recombinational effects [75]. Among the close relatives of P. marneffei in the subgenus Biverticillium, outcrossing has not been shown to occur, but self-fertility is common [76], [77], and the distribution of mating systems in the subgenus suggests that inbreeding may not reduce the evolutionary longevity of this group [78].

Another compelling scenario favouring sex recognises the opportunity presented by mating itself for dramatic shifts in morphology and physiology. A predominant view in the fungal literature is that sex occurs in otherwise mostly clonal fungi in response to stressful conditions [16], [79]. Sexually produced spores are often viable for long periods of time and are resistant to extreme environmental conditions [78], [80], [81], [82]. Regardless of the costs and benefits of recombination, P. marneffei might withstand stress by mating when it would otherwise not survive. If that were the case, recombination in P. marneffei might be clustered in space or time where or when stress occurs. Unfortunately, little is known about the natural ecology of P. marneffei, and any conditions that might allow mating to occur are unknown. Isolates are commonly recovered from bamboo rats, yet the epizoology of the fungus is poorly known including unknown routes of infection and unknown course and outcomes of the zoonosis. In other dimorphic fungi including Histoplasma and Blastomyces, mating in natural populations is also poorly known, but in these species it mating has long been studied outside of the host and at lower temperatures in vitro [83], [84]. Nevertheless, the association with small mammals may be the best starting point for a search for the natural sexual niche of P. marneffei.

Cryptic mating and inbreeding in P. marneffei has some parallels with other fungal and non-fungal eukaryotic pathogens [7]. There is growing support for high inbreeding in addition to asexual reproduction in Leishmania brasiliensis [85]. Experimental data support a role for within-vector recombination and it is thought that this restricted recombination in Leishmania results in sexual neighbourhoods of pathogen genotypes with high differentiation at multiple spatial scales [86]. In the malaria parasite Plasmodium falciparum, high inbreeding has been linked to faster emergence of drug resistance in some low infection intensity regions, but inbreeding has also emerged as a general property of P. falciparum populations regardless of infection intensity [87], [88]. In Toxoplasma gondii and Sarcocystis neurona, ‘clonal’ emergence is enabled by high rates of selfing and results in spatially structured populations [89], [90]. The most common fungal infection of humans, Candida albicans, also undergoes same-sex mating that can facilitate inbreeding and the ability to accelerate evolution via sex within clonal populations [91]. In C. neoformans, multiple ecological niches where recombination between strains with opposite mating type occurs have been found [92], [93]. However, most mating in C. neoformans globally is likely to occur between strains of the same mating type, and although there is as yet no indication of what factors alter the probability of this kind of inbreeding across natural populations, it may have facilitated the global emergence of a single mating type and highly clonal populations [15], [22], [94], [95]. Although some of the mechanisms underlying these inbreeding eukaryotic pathogens remain mysterious and likely differ between organisms, there is an emerging consistent pattern of clonality resulting from inbreeding rather than strictly asexual propagation even in the absence of a recognized sexual stage.

Supporting Information

Attachment 1

Attachment 2

Attachment 3

Attachment 4

Attachment 5

Attachment 6

Attachment 7

Attachment 8

Attachment 9


Zdroje

1. de WitR, BouvierT (2006) ‘Everything is everywhere, but, the environment selects’; what did Baas Becking and Beijerinck really say? Environ Microbiol 8: 755–758.

2. FinlayBJ (2002) Global dispersal of free-living microbial eukaryote species. Science 296: 1061–1063.

3. GreenJL, BohannanBJM, WhitakerRJ (2008) Microbial biogeography: From taxonomy to traits. Science 320: 1039–1043.

4. RydholmC, SzakacsG, LutzoniF (2006) Low genetic variation and no detectable population structure in Aspergillus fumigatus compared to closely related Neosartorya species. Eukaryot Cell 5: 650–657.

5. HenkDA, EagleCE, BrownK, Van Den BergMA, DyerPS, et al. (2011) Speciation despite globally overlapping distributions in Penicillium chrysogenum: the population genetics of Alexander Fleming's lucky fungus. Mol Ecol 4288–4301.

6. HeitmanJ (2006) Sexual reproduction and the evolution of microbial pathogens. Curr Biol 16: R711–R725.

7. HeitmanJ (2010) Evolution of eukaryotic microbial pathogens via covert sexual Reproduction. Cell Host Microbe 8: 86–99.

8. SunS, HeitmanJ (2011) Is sex necessary? BMC Biol 9: 56.

9. BarrettLG, ThrallPH, BurdonJJ, LindeCC (2008) Life history determines genetic structure and evolutionary potential of host-parasite interactions. Trends Ecol Evol 23: 678–685.

10. BruyndonckxN, HenryI, ChristeP, KerthG (2009) Spatio-temporal population genetic structure of the parasitic mite Spinturnix bechsteini is shaped by its own demography and the social system of its bat host. Mol Ecol 18: 3581–3592.

11. FisherMC, AanensenD, de HoogS, VanittanakomN (2004) Multilocus microsatellite typing system for Penicillium marneffei reveals spatially structured populations. J Clin Microbol 42: 5065–5069.

12. KasugaT, WhiteTJ, KoenigG, McewenJ, RestrepoA, et al. (2003) Phylogeography of the fungal pathogen Histoplasma capsulatum. Mol Ecol 12: 3383–3401.

13. MeeceJK, AndersonJL, FisherMC, HenkDA, SlossBL, et al. (2011) Population genetic structure of clinical and environmental isolates of Blastomyces dermatitidis, based on 27 polymorphic microsatellite markers. App Environ Microb 77: 5123–5131.

14. NeafseyDE, BarkerBM, SharptonTJ, StajichJE, ParkDJ, et al. (2010) Population genomic sequencing of Coccidioides fungi reveals recent hybridization and transposon control. Genome Res 20: 938–946.

15. SimwamiSP, KhayhanK, HenkDA, AanensenDM, BoekhoutT, et al. (2011) Low diversity Cryptococcus neoformans Variety grubii multilocus sequence types from Thailand are consistent with an ancestral african origin. PLoS Pathog 7: e1001343.

16. LeeSC, NiM, LiWJ, ShertzC, HeitmanJ (2010) The Evolution of Sex: a perspective from the fungal kingdom. Microbiol Mol Biol R 74: 298–340.

17. BilliardS, Lopez-VillavicencioM, DevierB, HoodME, FairheadC, et al. (2011) Having sex, yes, but with whom? Inferences from fungi on the evolution of anisogamy and mating types. Biol Rev 86: 421–442.

18. FisherMC, KoenigGL, WhiteTJ, San-BlasG, NegroniR, et al. (2001) Biogeographic range expansion into South America by Coccidioides immitis mirrors New World patterns of human migration. Proc Natl Acad Sci USA 98: 4558–4562.

19. LitvintsevaAP, LinX, TempletonI, HeitmanJ, MitchellTG (2007) Many globally isolated AD hybrid strains of Cryptococcus neoformans originated in Africa. PLoS Pathog 3: 1109–1117.

20. KasugaT, TaylorJW, WhiteTJ (1999) Phylogenetic relationships of varieties and geographical groups of the human pathogenic fungus Histoplasma capsulatum Darling. J Clin Microbiol 37: 653–663.

21. TaylorML, Chavez-TapiaCB, Rojas-MartinezA, Reyes-MontesMD, del ValleMB, et al. (2005) Geographical distribution of genetic polymorphism of the pathogen Histoplasma capsulatum isolated from infected bats, captured in a central zone of Mexico. FEMS Immunol Med Mic 45: 451–458.

22. LinXR, HullCM, HeitmanJ (2005) Sexual reproduction between partners of the same mating type in Cryptococcus neoformans. Nature 434: 1017–1021.

23. LitvintsevaAP, MarraRE, NielsenK, HeitmanJ, VilgalysR, et al. (2003) Evidence of sexual recombination among Cryptococcus neoformans serotype A isolates in sub-Saharan Africa. Eukaryot Cell 2: 1162–1168.

24. FisherMC, HanageWP, de HoogS, JohnsonE, SmithMD, et al. (2005) Low effective dispersal of asexual genotypes in heterogeneous landscapes by the endemic pathogen Penicillium marneffei. PLoS Pathog 1: 159–165.

25. WooPCY, ChongKTK, TseH, CaiJJ, LauCCY, et al. (2006) Genomic and experimental evidence for a potential sexual cycle in the pathogenic thermal dimorphic fungus Penicillium marneffei. FEBS Lett 580: 4976–4977.

26. SegretainG (1959) Penicillium marneffei n.sp., agent of a mycosis of the reticuloendothelial system. Mycopathologia 11: 327–353.

27. FisherMC, de HoogGS, VannittanakomN (2004) A highly discriminatory Multilocus Microsatellite Typing System (MLMT) for Penicillium marneffei. Mol Ecol Notes 5: 231–234.

28. PeakallR, SmousePE (2006) GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol Ecol Notes 6: 288–295.

29. JombartT (2008) adegenet: a R package for the multivariate analysis of genetic markers. Bioinformatics 24: 1403–1405.

30. WarrenDL, GlorRE, TurelliM (2008) Environmental niche equivalency versus conservatism: quantitative approaches to niche evolution. Evolution 62: 2868–2883.

31. LebloisR, EstoupA, RoussetF (2009) IBDSim: a computer program to simulate genotypic data under isolation by distance. Mol Ecol Resour 9: 107–109.

32. ForcheA, AlbyK, SchaeferD, JohnsonAD, BermanJ, et al. (2008) The parasexual cycle in Candida albicans provides an alternative pathway to meiosis for the formation of recombinant strains. PLoS Biol 6: 1084–1097.

33. GalaganJE, SelkerEU (2004) RIP: the evolutionary cost of genome defense. Trends Genet 20: 417–423.

34. KeeneyS, GirouxCN, KlecknerN (1997) Meiosis-specific DNA double-strand breaks are catalyzed by Spo11, a member of a widely conserved protein family. Cell 88: 375–384.

35. PoggelerS (2002) Genomic evidence for mating abilities in the asexual pathogen Aspergillus fumigatus. Curr Genet 42: 153–160.

36. TzungKW, WilliamsRM, SchererS, FederspielN, JonesT, et al. (2001) Genomic evidence for a complete sexual cycle in Candida albicans. Proc Natl Acad Sci USA 98: 3249–3253.

37. Dyer PS (2007) Sexual reproduction and significance of MAT in the Aspergilli. In: Sex in fungi: molecular determination and evolutionary implications. Heitman J, Kronstad, J.W, Taylor, J.W, and Casselton, L.A., editors. Washington, DC: ASM Press. pp. 123–142.

38. SelkerEU (1990) Premeiotic instability of repeated sequences in Neurospora crassa. Annu Rev Genet 24: 579–613.

39. ClutterbuckAJ (2011) Genomic evidence of repeat-induced point mutation (RIP) in filamentous ascomycetes. Fungal Genet Biol 48: 306–326.

40. BraumannI, van den BergM, KempkenF (2008) Repeat induced point mutation in two asexual fungi, Aspergillus niger and Penicillium chrysogenum. Curr Genet 53: 287–297.

41. FreitagM, WilliamsRL, KotheGO, SelkerEU (2002) A cytosine methyltransferase homologue is essential for repeat-induced point mutation in Neurospora crassa. Proc Natl Acad Sci USA 99: 8802–8807.

42. MalagnacF, WendelB, GoyonC, FaugeronG, ZicklerD, et al. (1997) A gene essential for de novo methylation and development in ascobolus reveals a novel type of eukaryotic DNA methyltransferase structure. Cell 91: 281–290.

43. O'GormanCM, FullerHT, DyerPS (2009) Discovery of a sexual cycle in the opportunistic fungal pathogen Aspergillus fumigatus. Nature 457: 471–U475.

44. NeuvegliseC, SarfatiJ, LatgeJP, ParisS (1996) Afut1, a retrotransposon-like element from Aspergillus fumigatus. Nucleic Acids Res 24: 1428–1434.

45. WrightS (1965) The interpretation of population structure by F-Statistics with special regard to systems of mating. Evolution 19: 395–420.

46. JombartT, DevillardS, DufourAB, PontierD (2008) Revealing cryptic spatial patterns in genetic variability by a new multivariate method. Heredity 101: 92–103.

47. JombartT, DrayS, DufourAB (2009) Finding essential scales of spatial variation in ecological data: a multivariate approach. Ecography 32: 161–168.

48. DettmanJR, TaylorJW (2004) Mutation and evolution of microsatellite loci in neurospora. Genetics 168: 1231–1248.

49. RuderferDM, PrattSC, SeidelHS, KruglyakL (2006) Population genomic analysis of outcrossing and recombination in yeast. Nat Genet 38: 1077–1081.

50. HudsonRR, KaplanNL (1985) Statistical Properties of the Number of Recombination Events in the History of a Sample of DNA-Sequences. Genetics 111: 147–164.

51. WrightS (1931) Evolution in Mendelian Populations. Genetics 16: 97–159.

52. LeslieJF, KleinKK (1996) Female fertility and mating type effects on effective population size and evolution in filamentous fungi. Genetics 144: 557–567.

53. FeilEJ, LiBC, AanensenDM, HanageWP, SprattBG (2004) eBURST: inferring patterns of evolutionary descent among clusters of related bacterial genotypes from multilocus sequence typing data. J Bacteriol 186: 1518–1530.

54. BoversM, HagenF, KuramaeEE, BoekhoutT (2008) Six monophyletic lineages identified within Cryptococcus neoformans and Cryptococcus gattii by multi-locus sequence typing. Fungal Genet Biol 45: 400–421.

55. FraserJA, GilesSS, WeninkEC, Geunes-BoyerSG, WrightJR, et al. (2005) Same-sex mating and the origin of the Vancouver Island Cryptococcus gattii outbreak. Nature 437: 1360–1364.

56. SaulN, KrockenbergerM, CarterD (2008) Evidence of recombination in mixed-mating-type and alpha-only populations of Cryptococcus gattii sourced from single Eucalyptus tree hollows. Eukaryot Cell 7: 727–734.

57. KaweckiTJ, EbertD (2004) Conceptual issues in local adaptation. Ecol Lett 7: 1225–1241.

58. RundleHD, NosilP (2005) Ecological speciation. Ecol Lett 8: 336–352.

59. SchluterD, ConteGL (2009) Genetics and ecological speciation. Proc Natl Acad Sci USA 106: 9955–9962.

60. EganSP, NosilP, FunkDJ (2008) Selection and genomic differentiation during ecological speciation: Isolating the contributions of host association via a comparative genome scan of Neochlamisus bebbianae leaf beetles. Evolution 62: 1162–1181.

61. HendryAP, NosilP, RiesebergLH (2007) The speed of ecological speciation. Funct Ecol 21: 455–464.

62. ViaS (2002) The ecological genetics of speciation. Am Nat 159: S1–S7.

63. GiraudT, RefregierG, Le GacM, de VienneDM, HoodME (2008) Speciation in fungi. Fungal Genet Biol 45: 791–802.

64. PhillipsSJ, AndersonRP, SchapireRE (2006) Maximum entropy modeling of species geographic distributions. Ecol Model 190: 231–259.

65. LiangL, CaoCW, WangWJ, LuoH, HuangSB, et al. (2011) Common Reservoirs for Penicillium marneffei Infection in Humans and Rodents, China. Emerg Infect Dis 17: 209–214.

66. WarrenDL, GlorRE, TurelliM (2010) ENMTools: a toolbox for comparative studies of environmental niche models. Ecography 33: 607–611.

67. GladieuxP, VerckenE, FontaineMC, HoodME, JonotO, et al. (2011) Maintenance of fungal pathogen species that are specialized to different hosts: allopatric divergence and introgression through secondary contact. Mol Biol Evol 28: 459–471.

68. HullCM, RaisnerRM, JohnsonAD (2000) Evidence for mating of the “asexual” yeast Candida albicans in a mammalian host. Science 289: 307–310.

69. HsuehYP, HeitmanJ (2008) Orchestration of sexual reproduction and virulence by the fungal mating-type locus. Curr Opin Microbiol 11: 517–524.

70. GiraudT, YocktengR, Lopez-VillavicencioM, RefregierG, HoodME (2008) Mating system of the anther smut fungus Microbotryum violaceum: selfing under heterothallism. Eukaryot Cell 7: 765–775.

71. KohnLM (2005) Mechanisms of fungal speciation. Annu Rev Phytopathol 43: 279–308.

72. Fisher MC (2007) The evolutionary implications of an asexual lifestyle manifested by Penicillium marneffei. In: Sex in fungi: molecular determination and evolutionary implications. Heitman J, Kronstad JW, Taylor JW, Casselton LA, editors. Washington, DC: ASM Press. pp. 201–212.

73. CharlesworthD, WillisJH (2009) The genetics of inbreeding depression. Nat Rev Genet 10: 783–796.

74. Goddard MR (2007) Why bother with sex? Answers from experiments with yeast and other organisms. In: Sex in fungi: molecular determination and evolutionary implications. Heitman J, Kronstad, J.W, Taylor, J.W, and Casselton, L.A., editors. Washington, DC: ASM Press. pp. 489–506.

75. BruggemanJ, DebetsAJM, WijngaardenPJ, deVisserJAGM, HoekstraRF (2003) Sex slows down the accumulation of deleterious mutations in the homothallic fungus Aspergillus nidulans. Genetics 164: 479–485.

76. Lopez-VillavicencioM, AguiletaG, GiraudT, de VienneDM, LacosteS, et al. (2010) Sex in Penicillium: Combined phylogenetic and experimental approaches. Fungal Genet Biol 47: 693–706.

77. GeiserDM (2009) Sexual structures in Aspergillus: morphology, importance and genomics. Med Mycolo 47 Suppl 1: S21–26.

78. NelsonMA, MetzenbergRL (1992) Sexual development genes of Neurospora Crassa. Genetics 132: 149–162.

79. BahnYS, XueCY, IdnurmA, RutherfordJC, HeitmanJ, et al. (2007) Sensing the environment: lessons from fungi. Nature Reviews Microbiology 5: 57–69.

80. WangP, PerfectJR, HeitmanJ (2000) The G-protein beta subunit GPB1 is required for mating and haploid fruiting in Cryptococcus neoformans. J Mol Cell Biol 20: 352–362.

81. GrishkanI, KorolAB, NevoE, WasserSP (2003) Ecological stress and sex evolution in soil microfungi. Proc Biol Sci 270: 13–18.

82. Alonso-MongeR, RomanE, AranaDM, PlaJ, NombelaC (2009) Fungi sensing environmental stress. Clin Microbiol Infec 15: 17–19.

83. Kwon-ChungKJ (1972) Emmonsiella capsulata: perfect state of Histoplasma capsulatum. Science 177: 368–369.

84. McDonoughES, LewisAL (1967) Blastomyces dermatitidis: production of the sexual stage. Science 156: 528–529.

85. RougeronV, De MeeusT, HideM, WaleckxE, BermudezH, et al. (2009) Extreme inbreeding in Leishmania braziliensis. Proc Natl Acad Sci USA 106: 10224–10229.

86. RougeronV, BanulsAL, CarmeB, SimonS, CouppieP, et al. (2011) Reproductive strategies and population structure in Leishmania: substantial amount of sex in Leishmania viannia guyanensis. Mol Ecol 20: 3116–3127.

87. DyeC, WilliamsBG (1997) Multigenic drug resistance among inbred malaria parasites. Proc R Soc Lond B Biol Sci 264: 61–67.

88. RazakandrainibeFG, DurandP, KoellaJC, De MeeusT, RoussetF, et al. (2005) “Clonal” population structure of the malaria agent Plasmodium falciparum in high-infection regions. Proc Natl Acad Sci USA 102: 17388–17393.

89. DubeyJP, VelmuruganGV, RajendranC, YabsleyMJ, ThomasNJ, et al. (2011) Genetic characterisation of Toxoplasma gondii in wildlife from North America revealed widespread and high prevalence of the fourth clonal type. International J Parasitol 41: 1139–1147.

90. WendteJM, MillerMA, LambournDM, MagargalSL, JessupDA, et al. (2010) Self-mating in the definitive host potentiates clonal outbreaks of the apicomplexan parasites Sarcocystis neurona and Toxoplasma gondii. PLoS Genet 6: e1001261.

91. AlbyK, SchaeferD, BennettRJ (2009) Homothallic and heterothallic mating in the opportunistic pathogen Candida albicans. Nature 460: 890–893.

92. NielsenK, De ObaldiaAL, HeitmanJ (2007) Cryptococcus neoformans mates on pigeon guano: Implications for the realized ecological niche and globalization. Eukaryot Cell 6: 949–959.

93. XueCY, TadaY, DoingXN, HeitmanJ (2007) The human fungal pathogen Cryptococcus can complete its sexual cycle during a pathogenic association with plants. Cell Host Microbe 1: 263–273.

94. BuiT, LinXR, MalikR, HeitmanJ, CarterD (2008) Isolates of Cryptococcus neoformans from Infected Animals Reveal Genetic Exchange in Unisexual, alpha Mating Type Populations. Eukaryot Cell 7: 1771–1780.

95. LinXR (2009) Cryptococcus neoformans: Morphogenesis, infection, and evolution. Infect Genet Evol 9: 401–416.

Štítky
Hygiena a epidemiologie Infekční lékařství Laboratoř

Článek vyšel v časopise

PLOS Pathogens


2012 Číslo 10
Nejčtenější tento týden
Nejčtenější v tomto čísle
Kurzy

Zvyšte si kvalifikaci online z pohodlí domova

Aktuální možnosti diagnostiky a léčby litiáz
nový kurz
Autoři: MUDr. Tomáš Ürge, PhD.

Střevní příprava před kolonoskopií
Autoři: MUDr. Klára Kmochová, Ph.D.

Závislosti moderní doby – digitální závislosti a hypnotika
Autoři: MUDr. Vladimír Kmoch

Aktuální možnosti diagnostiky a léčby AML a MDS nízkého rizika
Autoři: MUDr. Natália Podstavková

Jak diagnostikovat a efektivně léčit CHOPN v roce 2024
Autoři: doc. MUDr. Vladimír Koblížek, Ph.D.

Všechny kurzy
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se