Generators of phenotypic diversity in the evolution of pathogenic microorganisms

Comparative Genomic Analyses of Colletotrichum lindemuthianum Pathotypes with Different Virulence Levels and Lifestyles

Colletotrichum lindemuthianum is the most frequent pathogenic fungus of the common bean Phaseolus vulgaris. This filamentous fungus employs a hemibiotrophic nutrition/infection strategy, which is characteristic of many Colletotrichum species. Due to host-pathogen coevolution, C. lindemuthianum includes pathotypes with a diversity of virulence against differential common bean varieties. In this study, we performed comparative genomic analyses on three pathotypes with different virulence levels and a non-pathogenic pathotype, isolated from different geographical areas in Mexico. Our results revealed large genomes with high transposable element contents that have undergone expansions, generating intraspecific diversity. All the pathotypes exhibited a similar number of clusters of orthologous genes (COGs) and Gene Ontology (GO) terms. TFomes contain families that are typical in fungal genomes; however, they show different contents between pathotypes, mainly in transcription factors with the fungal-specific TF and Zn2Cys6 domains. Peptidase families mainly contain abundant serine peptidases, metallopeptidases, and cysteine peptidases. In the secretomes, the number of genes differed between the pathotypes, with a high percentage of candidate effectors. Both the virulence gene and CAZyme gene content for each pathotype was abundant and diverse, and the latter was enriched in hemicellulolytic enzymes. We provide new insights into the nature of intraspecific diversity among C. lindemuthianum pathotypes and the origin of their ability to rapidly adapt to genetic changes in its host and environmental conditions.

Fungal and fungal-like diversity present in ornithogenically influenced maritime Antarctic soils assessed using metabarcoding

We assessed soil fungal and fungal-like diversity using metabarcoding in ornithogenically influenced soils around nests of the bird species Phalacrocorax atriceps, Macronectes giganteus, Pygoscelis antarcticus, and Pygoscelis adelie on the South Shetland Islands, maritime Antarctic. A total of 1,392,784 fungal DNA reads was obtained and assigned to 186 amplicon sequence variants (ASVs). The dominant fungal phylum was Ascomycota, followed by Basidiomycota, Chytridiomycota, Blastocladiomycota, Rozellomycota, Mortierellomycota, Monoblepharomycota, Aphelidiomycota, Basidiobolomycota, Mucoromycota, and the fungal-like Oomycota (Stramenopila), in rank order. Antarctomyces sp., Blastocladiomycota sp., Pseudogymnoascus pannorum, Microascaceae sp., Mortierella sp., Lobulomycetales sp., Sordariomycetes sp., Fungal sp., Rhizophydiales sp., Pseudeurotiaceae sp., Chytridiomycota sp. 1, Filobasidiella sp., Tausonia pullulans, Betamyces sp., and Leucosporidium sp. were the most abundant assigned taxa. The fungal assemblages present in the different ornithogenically influenced soils displayed different diversity indices. However, in general, we detected high fungal diversity and few taxa shared between the samples. Despite the polyextreme environmental conditions experienced in these Antarctic soils, the metabarcoding approach detected a rich and complex fungal community dominated by saprophytes, but with some pathogenic taxa also present. The community was dominated by psychrophilic and psychrotolerant taxa, some apparently endemic to Antarctica, and those identified only at higher taxonomic levels, which may represent currently undescribed fungi. The mycobiome detected included taxa characterized by different ecological roles, including saprotrophic, human- and animal-associated, phytopathogenic, mutualistic, and cosmopolitan. These fungi may potentially be dispersed by birds or in the air column over great distances, including between different regions within Antarctica and from South America, Africa, and Oceania.

Differential Carbon Catabolite Repression and Hemicellulolytic Ability among Pathotypes of Colletotrichum lindemuthianum against Natural Plant Substrates

Colletotrichum lindemuthianum is a phytopathogenic fungus that causes anthracnose in common beans (Phaseolus vulgaris) and presents a great diversity of pathotypes with different levels of virulence against bean varieties worldwide. The purpose of this study was to establish whether pathotypic diversity is associated with differences in the mycelial growth and secretion of plant-cell-wall-degrading enzymes (PCWDEs). We evaluated growth, hemicellulase and cellulase activity, and PCWDE secretion in four pathotypes of C. lindemuthianum in cultures with glucose, bean hypocotyls and green beans of P. vulgaris, and water hyacinth (Eichhornia crassipes). The results showed differences in the mycelial growth, hemicellulolytic activity, and PCWDE secretion among the pathotypes. Glucose was not the preferred carbon source for the best mycelial growth in all pathotypes, each of which showed a unique PCWDE secretion profile, indicating different levels of carbon catabolite regulation (CCR). The pathotypes showed a high differential hemicellulolytic capacity to degrade host and water hyacinth tissues, suggesting CCR by pentoses and that there are differences in the absorption and metabolism of different monosaccharides and/or disaccharides. We propose that different levels of CCR could optimize growth in different host tissues and could allow for consortium behavior in interactions with bean crops.

Population Genetic Structure of the Rubber Tree Powdery Mildew Pathogen (Erysiphe quercicola) from China

In order to manage agricultural pathogens, it is crucial to understand the population structure underlying epidemics. Rubber tree powdery mildew, caused by Erysiphe quercicola, is a serious threat to rubber plantations worldwide, especially in subtropical environments including all rubber tree-growing regions in China. However, the population structure of the pathogen is uncertain. In this study, 16 polymorphic microsatellite markers were used to genotype powdery mildew samples from the main rubber tree-growing regions including Yunnan (YN), Hainan (HN), western Guangdong (WG), and eastern Guangdong (EG). YN had higher genotypic diversity (Simpson's indices), genotypic evenness, Nei's gene diversity, allelic richness, and private allelic richness than the other regions. Cluster analysis, discriminant analysis of principal components, pairwise divergence, and shared multilocus genotype analyses all showed that YN differed significantly from the other regions. The genetic differentiation was small among the other three regions (HN, WG, and EG). Analysis of molecular variance indicated that the variability among regions accounted for 22.37% of the total variability. Genetic differentiation was significantly positively correlated (Rxy = 0.772, P = 0.001) with geographic distance. Linkage equilibrium analysis suggested possible occurrence of sexual recombination although asexual reproduction predominates in E. quercicola. The results suggested that although significant genetic differentiation of E. quercicola occurred between YN and the other regions, pathogen populations from the other three regions lacked genetic differentiation.

Whole genome resequencing reveal patterns of genetic variation within Colletotrichum acutatum species complex from rubber trees in China

The Colletotrichum acutatum species complex possesses a diverse number of important traits, such as a wide host range and host preference, different modes of reproduction, and different strategies of host infection. Research using comparative genomics has attempted to find correlations between these traits. Here, we used multi-locus techniques and gene genealogical concordance analysis to investigate the phylogenetic relationships and taxonomic status of the Colletotrichum acutatum species complex using field isolates obtained from rubber trees. The results revealed that the dominant species was C. australisinense, followed by C. bannaense, while strain YNJH17109 was identified as C. laticiphilum. The taxonomic status of strains YNLC510 and YNLC511 was undetermined. Using whole-genome single nucleotide polymorphism data to analyze population structure, 18 strains of C. australisinense were subsequently divided into four populations, one of which was derived from an admixture of two populations. In addition, the strains LD1687, GD1628, and YNLC516, did not belong to any populations, and were considered to be admixtures of two or more populations. A split decomposition network analysis also provided evidence for genetic recombination within Colletotrichum acutatum species complex from rubber trees in China. Overall, a weak phylogeographic sub-structure was observed. Analysis also revealed significant differences in morphological characters and levels of virulence between populations.

Genetic Analysis of Colletotrichum siamense Populations from Different Hosts and Counties in Hainan, China, Using Microsatellite Markers

Colletotrichum siamense was demonstrated as the dominant species among Colletotrichum spp. that infected rubber tree, areca palm, and coffee in Hainan, China. However, the extent of genetic differentiation within the species C. siamense in relation to geographical regions and host species is not known. In this study, 112 C. siamense isolates were genotyped with 12 microsatellite markers. In total, there were 99 multilocus genotypes. Results from permutational multivariate analysis of variance and analysis of molecular variance indicated that there was no significant genetic differentiation between fungal populations with respect to host, location (county), and year. Discriminant analysis of principal components and STRUCTURE analysis showed that C. siamense isolates grouped into three clusters; further analysis confirmed that there were significant (P < 0.001) genetic differences among the three clusters. However, each cluster had isolates from different hosts, counties, or years, supporting the lack of genetic differentiation with respect to host, county, and year. Statistical analyses of allelic associations indicated some evidence for recombination within the populations defined on the basis of host or county. The present findings provide insights into the genetic structure of C. siamense on the three perennial host species in Hainan and suggest that the disease on these three crops can be effectively considered as one disease and, hence, needs to be controlled simultaneously in mixed plantations.

Sterylglucosides in Fungi

Sterylglucosides (SGs) are sterol conjugates widely distributed in nature. Although their universal presence in all living organisms suggests the importance of this kind of glycolipids, they are yet poorly understood. The glycosylation of sterols confers a more hydrophilic character, modifying biophysical properties of cell membranes and altering immunogenicity of the cells. In fungi, SGs regulate different cell pathways to help overcome oxygen and pH challenges, as well as help to accomplish cell recycling and other membrane functions. At the same time, the level of these lipids is highly controlled, especially in wild-type fungi. In addition, modulating SGs metabolism is becoming a novel tool for vaccine and antifungal development. In the present review, we bring together multiple observations to emphasize the underestimated importance of SGs for fungal cell functions.

Drug Resistance and Evolvability in an Emerging Human Fungal Pathogen

Over the past decade, Candida auris has emerged as a highly transmissible human fungal pathogen. Because of its ability to transmit between patients in hospitals and its ability to rapidly develop drug resistance, C. auris presents unique challenges. However, at a genetic and genomic level we still understand relatively little about how drug resistance develops in this pathogen. Burrack et al. use experimental evolution and whole-genome sequencing to identify mutations correlated with fluconazole resistance in C. auris. They identify interesting genomic features, including highly plastic subtelomeric regions and whole chromosomal and segmental aneuploidies. Excitingly, they also identify the first example of a hypermutator strain in C. auris. In comparison with the model human fungal pathogen Candida albicans, C. auris is more likely to undergo mutation and less likely to undergo copy number variation in response to drug selection, which may be linked to differences in base ploidy level.

Role of Sexual Reproduction in the Evolution of the Wheat Stripe Rust Fungus Races in China

Experimental and population genetic approaches have reshaped our view of how fungal pathogens reproduce, with consequences for our understanding of fungal invasions. Puccinia striiformis f. sp. tritici, the causal agent of stripe rust, poses a severe threat to wheat production worldwide. The sexual stage of P. striiformis f. sp. tritici was discovered >10 years ago, but how it affects the evolution of the pathogen, especially the emergence of the new virulent races, remains largely unknown. Here, using population genetic analyses, we demonstrate that sexual reproduction plays an important role in the evolution of P. striiformis f. sp. tritici races in China, specifically the newly emerged and devastating race virulent to resistance gene Yr26, which is widely used in China and exerts strong selective pressure on the pathogen population. Association analysis identified six genes encoding secreted proteins as candidates for virulence on wheat cultivars carrying the Yr26 resistance gene. Our results highlight the important role of sexual reproduction and selection exerted by hosts in the emergence of new virulent races in China.

Genome Analysis of the Broad Host Range Necrotroph Nalanthamala psidii Highlights Genes Associated With Virulence

Guava wilt disease is caused by the fungus Nalanthamala psidii. The wilt disease results in large-scale destruction of orchards in South Africa, Taiwan, and several Southeast Asian countries. De novo assembly, annotation, and in-depth analysis of the N. psidii genome were carried out to facilitate the identification of characteristics associated with pathogenicity and pathogen evolution. The predicted secretome revealed a range of CAZymes, proteases, lipases and peroxidases associated with plant cell wall degradation, nutrient acquisition, and disease development. Further analysis of the N. psidii carbohydrate-active enzyme profile exposed the broad-spectrum necrotrophic lifestyle of the pathogen, which was corroborated by the identification of putative effectors and secondary metabolites with the potential to induce tissue necrosis and cell surface-dependent immune responses. Putative regulatory proteins including transcription factors and kinases were identified in addition to transporters potentially involved in the secretion of secondary metabolites. Transporters identified included important ABC and MFS transporters involved in the efflux of fungicides. Analysis of the repetitive landscape and the detection of mechanisms linked to reproduction such as het and mating genes rendered insights into the biological complexity and evolutionary potential of N. psidii as guava pathogen. Hence, the assembly and annotation of the N. psidii genome provided a valuable platform to explore the pathogenic potential and necrotrophic lifestyle of the guava wilt pathogen.

Fungi as mediators linking organisms and ecosystems

Fungi form a major and diverse component of most ecosystems on Earth. They are both micro and macroorganisms with high and varying functional diversity as well as great variation in dispersal modes. With our growing knowledge of microbial biogeography, it has become increasingly clear that fungal assembly patterns and processes differ from other microorganisms such as bacteria, but also from macroorganisms such as plants. The success of fungi as organisms and their influence on the environment lies in their ability to span multiple dimensions of time, space, and biological interactions, that is not rivalled by other organism groups. There is also growing evidence that fungi mediate links between different organisms and ecosystems, with the potential to affect the macroecology and evolution of those organisms. This suggests that fungal interactions are an ecological driving force, interconnecting different levels of biological and ecological organisation of their hosts, competitors, and antagonists with the environment and ecosystem functioning. Here we review these emerging lines of evidence by focusing on the dynamics of fungal interactions with other organism groups across various ecosystems. We conclude that the mediating role of fungi through their complex and dynamic ecological interactions underlie their importance and ubiquity across Earth's ecosystems.

Fungicide Resistance Evolution and Detection in Plant Pathogens: Plasmopara viticola as a Case Study

The use of single-site fungicides to control plant pathogens in the agroecosystem can be associated with an increased selection of resistance. The evolution of resistance represents one of the biggest challenges in disease control. In vineyards, frequent applications of fungicides are carried out every season for multiple years. The agronomic risk of developing fungicide resistance is, therefore, high. Plasmopara viticola, the causal agent of grapevine downy mildew, is a high risk pathogen associated with the development of fungicide resistance. P. viticola has developed resistance to most of the fungicide classes used and constitutes one of the most important threats for grapevine production. The goals of this review are to describe fungicide resistance evolution in P. viticola populations and how to conduct proper monitoring activities. Different methods have been developed for phenotyping and genotyping P. viticola for fungicide resistance and the different phases of resistance evolution and life cycles of the pathogen are discussed, to provide a full monitoring toolkit to limit the spread of resistance. A detailed revision of the available tools will help in shaping and harmonizing the monitoring activities between countries and organizations.

Draft genomes of Perkinsus olseni and Perkinsus chesapeaki reveal polyploidy and regional differences in heterozygosity

Perkinsus spp. parasites have significant impact on aquaculture and wild mollusc populations. We sequenced the genomes of five monoclonal isolates of Perkinsus olseni and one Perkinsus chesapeaki from international sources. Sequence analysis revealed similar levels of repetitive sequence within species, a polyploid genome structure, and substantially higher heterozygosity in Oceanian-sourced isolates. We also identified tandem replication of the rRNA transcriptional unit, with high strain variation. Characterized gene content was broadly similar amongst all Perkinsus spp. but P. olseni Oceanian isolates contained an elevated number of genes compared to other P. olseni isolates and cox3 could not be identified in any Perkinsus spp. sequence. Phylogenetics and average nucleotide identity scans were consistent with all P. olseni isolates being within one species. These are the first genome sequences generated for both P. olseni and P. chesapeaki and will allow future advances in diagnostic design and population genomics of these important aquatic parasites.

Animal Models of Cryptococcus neoformans in Identifying Immune Parameters Associated With Primary Infection and Reactivation of Latent Infection

Cryptococcus species are environmental fungal pathogens and the causative agents of cryptococcosis. Infection occurs upon inhalation of infectious particles, which proliferate in the lung causing a primary infection. From this primary lung infection, fungal cells can eventually disseminate to other organs, particularly the brain, causing lethal meningoencephalitis. However, in most cases, the primary infection resolves with the formation of a lung granuloma. Upon severe immunodeficiency, dormant cryptococcal cells will start proliferating in the lung granuloma and eventually will disseminate to the brain. Many investigators have sought to study the protective host immune response to this pathogen in search of host parameters that keep the proliferation of cryptococcal cells under control. The majority of the work assimilates research carried out using the primary infection animal model, mainly because a reactivation model has been available only very recently. This review will focus on anti-cryptococcal immunity in both the primary and reactivation models. An understanding of the differences in host immunity between the primary and reactivation models will help to define the key host parameters that control the infections and are important for the research and development of new therapeutic and vaccine strategies against cryptococcosis.

Genetic Diversity and Differentiation in Phoma betae Populations on Table Beet in New York and Washington States

Phoma betae is an important seedborne pathogen of table beet worldwide that is capable of causing foliar, root, and damping-off diseases. Ten microsatellite and mating type markers were developed to investigate the genetics of P. betae populations in table beet root crops in New York and in table beet seed crops in Washington, from where table beet seed is predominantly sourced. The markers were used to characterize 175 isolates comprising five P. betae populations (two from New York and three from Washington), and they were highly polymorphic with an allelic range of 4 to 33 and an average of 11.7 alleles per locus. All populations had high genotypic diversity (Simpson's complement index = 0.857 to 0.924) and moderate allelic diversity (Nei's unbiased gene diversity = 0.582 to 0.653). Greater differentiation observed between populations from the two states compared with populations within the same state suggested that an external inoculum source, such as windblown ascospores, may be homogenizing the populations. However, most genetic diversity (87%) was among individual isolates within populations (pairwise index of population differentiation = 0.127; P = 0.001), suggesting that local within-field inoculum source(s), such as infested field debris or infected weeds, may also be important in initiating disease outbreaks. Standardized index of association, proportion of compatible pairs of loci, and mating type ratio calculations showed evidence for a mixed reproduction mode in all populations. These findings could be useful in designing more effective management strategies for diseases caused by P. betae in table beet production.

Mitotic Recombination and Rapid Genome Evolution in the Invasive Forest Pathogen Phytophthora ramorum

Invasive alien species often have reduced genetic diversity and must adapt to new environments. Given the success of many invasions, this is sometimes called the genetic paradox of invasion. Phytophthora ramorum is invasive, limited to asexual reproduction within four lineages, and presumed clonal. It is responsible for sudden oak death in the United States, sudden larch death in Europe, and ramorum blight in North America and Europe. We sequenced the genomes of 107 isolates to determine how this pathogen can overcome the invasion paradox. Mitotic recombination (MR) associated with transposons and low gene density has generated runs of homozygosity (ROH) affecting 2,698 genes, resulting in novel genotypic diversity within the lineages. One ROH enriched in effectors was fixed in the NA1 lineage. An independent ROH affected the same scaffold in the EU1 lineage, suggesting an MR hot spot and a selection target. Differences in host infection between EU1 isolates with and without the ROH suggest that they may differ in aggressiveness. Non-core regions (not shared by all lineages) had signatures of accelerated evolution and were enriched in putative pathogenicity genes and transposons. There was a striking pattern of gene loss, including all effectors, in the non-core EU2 genome. Positive selection was observed in 8.0% of RxLR and 18.8% of Crinkler effector genes compared with 0.9% of the core eukaryotic gene set. We conclude that the P. ramorum lineages are diverging via a rapidly evolving non-core genome and that the invasive asexual lineages are not clonal, but display genotypic diversity caused by MR.IMPORTANCE Alien species are often successful invaders in new environments, despite the introduction of a few isolates with a reduced genetic pool. This is called the genetic paradox of invasion. We found two mechanisms by which the invasive forest pathogen causing sudden oak and sudden larch death can evolve. Extensive mitotic recombination producing runs of homozygosity generates genotypic diversity even in the absence of sexual reproduction, and rapid turnover of genes in the non-core, or nonessential portion of genome not shared by all isolates, allows pathogenicity genes to evolve rapidly or be eliminated while retaining essential genes. Mitotic recombination events occur in genomic hot spots, resulting in similar ROH patterns in different isolates or groups; one ROH, independently generated in two different groups, was enriched in pathogenicity genes and may be a target for selection. This provides important insights into the evolution of invasive alien pathogens and their potential for adaptation and future persistence.

Hemizygosity Enables a Mutational Transition Governing Fungal Virulence and Commensalism

Candida albicans is a commensal fungus of human gastrointestinal and reproductive tracts, but also causes life-threatening systemic infections. The balance between colonization and pathogenesis is associated with phenotypic plasticity, with alternative cell states producing different outcomes in a mammalian host. Here, we reveal that gene dosage of a master transcription factor regulates cell differentiation in diploid C. albicans cells, as EFG1 hemizygous cells undergo a phenotypic transition inaccessible to "wild-type" cells with two functional EFG1 alleles. Notably, clinical isolates are often EFG1 hemizygous and thus licensed to undergo this transition. Phenotypic change corresponds to high-frequency loss of the functional EFG1 allele via de novo mutation or gene conversion events. This phenomenon also occurs during passaging in the gastrointestinal tract with the resulting cell type being hypercompetitive for commensal and systemic infections. A "two-hit" genetic model therefore underlies a key phenotypic transition in C. albicans that enables adaptation to host niches.

Pheromone independent unisexual development in Cryptococcus neoformans

The fungus Cryptococcus neoformans can undergo a-α bisexual and unisexual reproduction. Completion of both sexual reproduction modes requires similar cellular differentiation processes and meiosis. Although bisexual reproduction generates equal number of a and α progeny and is far more efficient than unisexual reproduction under mating-inducing laboratory conditions, the α mating type dominates in nature. Population genetic studies suggest that unisexual reproduction by α isolates might have contributed to this sharply skewed distribution of the mating types. However, the predominance of the α mating type and the seemingly inefficient unisexual reproduction observed under laboratory conditions present a conundrum. Here, we discovered a previously unrecognized condition that promotes unisexual reproduction while suppressing bisexual reproduction. Pheromone is the principal stimulus for bisexual development in Cryptococcus. Interestingly, pheromone and other components of the pheromone pathway, including the key transcription factor Mat2, are not necessary but rather inhibitory for Cryptococcus to complete its unisexual cycle under this condition. The inactivation of the pheromone pathway promotes unisexual reproduction despite the essential role of this pathway in non-self-recognition during bisexual reproduction. Nonetheless, the requirement for the known filamentation regulator Znf2 and the expression of hyphal or basidium specific proteins remain the same for pheromone-dependent or independent sexual reproduction. Transcriptome analyses and an insertional mutagenesis screen in mat2Δ identified calcineurin being essential for this process. We further found that Znf2 and calcineurin work cooperatively in controlling unisexual development in this fungus. These findings indicate that Mat2 acts as a repressor of pheromone-independent unisexual development while serving as an activator for a-α bisexual development. The bi-functionality of Mat2 might have allowed it to act as a toggle switch for the mode of sexual development in this ubiquitous eukaryotic microbe.

A Population Genomics Approach to Assessing the Genetic Basis of Within-Host Microevolution Underlying Recurrent Cryptococcal Meningitis Infection

Recurrence of meningitis due to Cryptococcus neoformans after treatment causes substantial mortality in HIV/AIDS patients across sub-Saharan Africa. In order to determine whether recurrence occurred due to relapse of the original infecting isolate or reinfection with a different isolate weeks or months after initial treatment, we used whole-genome sequencing (WGS) to assess the genetic basis of infection in 17 HIV-infected individuals with recurrent cryptococcal meningitis (CM). Comparisons revealed a clonal relationship for 15 pairs of isolates recovered before and after recurrence showing relapse of the original infection. The two remaining pairs showed high levels of genetic heterogeneity; in one pair we found this to be a result of infection by mixed genotypes, while the second was a result of nonsense mutations in the gene encoding the DNA mismatch repair proteins MSH2, MSH5, and RAD5. These nonsense mutations led to a hypermutator state, leading to dramatically elevated rates of synonymous and nonsynonymous substitutions. Hypermutator phenotypes owing to nonsense mutations in these genes have not previously been reported in C. neoformans, and represent a novel pathway for rapid within-host adaptation and evolution of resistance to first-line antifungal drugs.

Exploring the genomic diversity of black yeasts and relatives (Chaetothyriales, Ascomycota)

The order Chaetothyriales (Pezizomycotina, Ascomycetes) harbours obligatorily melanised fungi and includes numerous etiologic agents of chromoblastomycosis, phaeohyphomycosis and other diseases of vertebrate hosts. Diseases range from mild cutaneous to fatal cerebral or disseminated infections and affect humans and cold-blooded animals globally. In addition, Chaetothyriales comprise species with aquatic, rock-inhabiting, ant-associated, and mycoparasitic life-styles, as well as species that tolerate toxic compounds, suggesting a high degree of versatile extremotolerance. To understand their biology and divergent niche occupation, we sequenced and annotated a set of 23 genomes of main the human opportunists within the Chaetothyriales as well as related environmental species. Our analyses included fungi with diverse life-styles, namely opportunistic pathogens and closely related saprobes, to identify genomic adaptations related to pathogenesis. Furthermore, ecological preferences of Chaetothyriales were analysed, in conjuncture with the order-level phylogeny based on conserved ribosomal genes. General characteristics, phylogenomic relationships, transposable elements, sex-related genes, protein family evolution, genes related to protein degradation (MEROPS), carbohydrate-active enzymes (CAZymes), melanin synthesis and secondary metabolism were investigated and compared between species. Genome assemblies varied from 25.81 Mb (Capronia coronata) to 43.03 Mb (Cladophialophora immunda). The bantiana-clade contained the highest number of predicted genes (12 817 on average) as well as larger genomes. We found a low content of mobile elements, with DNA transposons from Tc1/Mariner superfamily being the most abundant across analysed species. Additionally, we identified a reduction of carbohydrate degrading enzymes, specifically many of the Glycosyl Hydrolase (GH) class, while most of the Pectin Lyase (PL) genes were lost in etiological agents of chromoblastomycosis and phaeohyphomycosis. An expansion was found in protein degrading peptidase enzyme families S12 (serine-type D-Ala-D-Ala carboxypeptidases) and M38 (isoaspartyl dipeptidases). Based on genomic information, a wide range of abilities of melanin biosynthesis was revealed; genes related to metabolically distinct DHN, DOPA and pyomelanin pathways were identified. The MAT (MAting Type) locus and other sex-related genes were recognized in all 23 black fungi. Members of the asexual genera Fonsecaea and Cladophialophora appear to be heterothallic with a single copy of either MAT-1-1 or MAT-1-2 in each individual. All Capronia species are homothallic as both MAT1-1 and MAT1-2 genes were found in each single genome. The genomic synteny of the MAT-locus flanking genes (SLA2-APN2-COX13) is not conserved in black fungi as is commonly observed in Eurotiomycetes, indicating a unique genomic context for MAT in those species. The heterokaryon (het) genes expansion associated with the low selective pressure at the MAT-locus suggests that a parasexual cycle may play an important role in generating diversity among those fungi.

Population Genetic Structure of Rhizoctonia oryzae-sativae from Rice in Latin America and Its Adaptive Potential to Emerge as a Pathogen on Urochloa Pastures

The fungus Rhizoctonia oryzae-sativae is an important pathogen that causes the aggregated sheath spot disease on rice. In this study, we investigated the genetic structure of rice-adapted populations of R. oryzae-sativae sampled from traditional rice-cropping areas from the Paraíba Valley, São Paulo, Brazil, and from Meta, in the Colombian Llanos, in South America. We used five microsatellite loci to measure population differentiation and infer the pathogen's reproductive system. Gene flow was detected among the three populations of R. oryzae-sativae from lowland rice in Brazil and Colombia. In contrast, a lack of gene flow was observed between the lowland and the upland rice populations of the pathogen. Evidence of sexual reproduction including low clonality, Hardy-Weinberg equilibrium within loci and gametic equilibrium between loci, indicated the predominance of a mixed reproductive system in all populations. In addition, we assessed the adaptive potential of the Brazilian populations of R. oryzae-sativae to emerge as a pathogen to Urochloa spp. (signalgrass) based on greenhouse aggressiveness assays. The Brazilian populations of R. oryzae-sativae were probably only incipiently adapted as a pathogen to Urochloa spp. Comparison between R_ST and Q_ST showed the predominance of diversifying selection in the divergence between the two populations of R. oryzae-sativae from Brazil.

Is Predominant Clonal Evolution a Common Evolutionary Adaptation to Parasitism in Pathogenic Parasitic Protozoa, Fungi, Bacteria, and Viruses?

We propose that predominant clonal evolution (PCE) in microbial pathogens be defined as restrained recombination on an evolutionary scale, with genetic exchange scarce enough to not break the prevalent pattern of clonal population structure. The main features of PCE are (1) strong linkage disequilibrium, (2) the widespread occurrence of stable genetic clusters blurred by occasional bouts of genetic exchange ('near-clades'), (3) the existence of a "clonality threshold", beyond which recombination is efficiently countered by PCE, and near-clades irreversibly diverge. We hypothesize that the PCE features are not mainly due to natural selection but also chiefly originate from in-built genetic properties of pathogens. We show that the PCE model obtains even in microbes that have been considered as 'highly recombining', such as Neisseria meningitidis, and that some clonality features are observed even in Plasmodium, which has been long described as panmictic. Lastly, we provide evidence that PCE features are also observed in viruses, taking into account their extremely fast genetic turnover. The PCE model provides a convenient population genetic framework for any kind of micropathogen. It makes it possible to describe convenient units of analysis (clones and near-clades) for all applied studies. Due to PCE features, these units of analysis are stable in space and time, and clearly delimited. The PCE model opens up the possibility of revisiting the problem of species definition in these organisms. We hypothesize that PCE constitutes a major evolutionary strategy for protozoa, fungi, bacteria, and viruses to adapt to parasitism.

Population Genetics of Hirsutella rhossiliensis, a Dominant Parasite of Cyst Nematode Juveniles on a Continental Scale

Hirsutella rhossiliensis is a parasite of juvenile nematodes, effective against a diversity of plant-parasitic nematodes. Its global distribution on various nematode hosts and its genetic variation for several geographic regions have been reported, while the global population genetic structure and factors underlying patterns of genetic variation of H. rhossiliensis are unclear. In this study, 87 H. rhossiliensis strains from five nematode species (Globodera sp., Criconemella xenoplax, Rotylenchus robustus, Heterodera schachtii, and Heterodera glycines) in Europe, the United States, and China were investigated by multilocus sequence analyses. A total of 280 variable sites (frequency, 0.6%) at eight loci and six clustering in high accordance with geographic populations or host nematode-associated populations were identified. Although H. rhossiliensis is currently recognized as an asexual fungus, recombination events were frequently detected. In addition, significant genetic isolation by geography and nematode hosts was revealed. Overall, our analyses showed that recombination, geographic isolation, and nematode host adaptation have played significant roles in the evolutionary history of H. rhossiliensis IMPORTANCE: H. rhossiliensis has great potential for use as a biocontrol agent to control nematodes in a sustainable manner as an endoparasitic fungus. Therefore, this study has important implications for the use of H. rhossiliensis as a biocontrol agent and provides interesting insights into the biology of this species.

Gene Network Polymorphism Illuminates Loss and Retention of Novel RNAi Silencing Components in the Cryptococcus Pathogenic Species Complex

RNAi is a ubiquitous pathway that serves central functions throughout eukaryotes, including maintenance of genome stability and repression of transposon expression and movement. However, a number of organisms have lost their RNAi pathways, including the model yeast Saccharomyces cerevisiae, the maize pathogen Ustilago maydis, the human pathogen Cryptococcus deuterogattii, and some human parasite pathogens, suggesting there may be adaptive benefits associated with both retention and loss of RNAi. By comparing the RNAi-deficient genome of the Pacific Northwest Outbreak C. deuterogattii strain R265 with the RNAi-proficient genomes of the Cryptococcus pathogenic species complex, we identified a set of conserved genes that were lost in R265 and all other C. deuterogattii isolates examined. Genetic and molecular analyses reveal several of these lost genes play roles in RNAi pathways. Four novel components were examined further. Znf3 (a zinc finger protein) and Qip1 (a homolog of N. crassa Qip) were found to be essential for RNAi, while Cpr2 (a constitutive pheromone receptor) and Fzc28 (a transcription factor) are involved in sex-induced but not mitosis-induced silencing. Our results demonstrate that the mitotic and sex-induced RNAi pathways rely on the same core components, but sex-induced silencing may be a more specific, highly induced variant that involves additional specialized or regulatory components. Our studies further illustrate how gene network polymorphisms involving known components of key cellular pathways can inform identification of novel elements and suggest that RNAi loss may have been a core event in the speciation of C. deuterogattii and possibly contributed to its pathogenic trajectory.

Genome instability and mutations provoked by transposon movement are counteracted by novel defense mechanisms in organisms as diverse as fungi, plants, and mammals. In the human fungal pathogen Cryptococcus neoformans, an RNAi silencing pathway operates to defend the genome against mobile elements and transgene repeats. RNAi silencing pathways are conserved in the Cryptococcus pathogenic species complex and are mediated by canonical RNAi components. Surprisingly, several of these components are missing from all analyzed C. deuterogattii VGII strains, the molecular type responsible for the North American Pacific Northwest outbreak. To identify novel components of the RNAi pathways, we surveyed the reference genomes of C. deuterogattii, C. gattii, C. neoformans, and C. deneoformans. We identified 14 otherwise conserved genes missing in R265, including the RDP1, AGO1, and DCR1 canonical RNAi components, and focused on four potentially novel RNAi components: ZNF3, QIP1, CPR2, and FZC28. We found that Znf3 and Qip1 are both required for mitotic- and sex-induced silencing, while Cpr2 and Fzc28 contribute to sex-induced but not mitosis-induced silencing. Our studies reveal elements of RNAi pathways that operate to defend the genome during sexual development and vegetative growth and illustrate the power of network polymorphisms to illuminate novel components of biological pathways.

Stochastic developmental variation, an epigenetic source of phenotypic diversity with far-reaching biological consequences

This article reviews the production of different phenotypes from the same genotype in the same environment by stochastic cellular events, nonlinear mechanisms during patterning and morphogenesis, and probabilistic self-reinforcing circuitries in the adult life. These aspects of phenotypic variation are summarized under the term 'stochastic developmental variation' (SDV) in the following. In the past, SDV has been viewed primarily as a nuisance, impairing laboratory experiments, pharmaceutical testing, and true-to-type breeding. This article also emphasizes the positive biological effects of SDV and discusses implications for genotype-to-phenotype mapping, biological individuation, ecology, evolution, and applied biology. There is strong evidence from experiments with genetically identical organisms performed in narrowly standardized laboratory set-ups that SDV is a source of phenotypic variation in its own right aside from genetic variation and environmental variation. It is obviously mediated by molecular and higher-order epigenetic mechanisms. Comparison of SDV in animals, plants, fungi, protists, bacteria, archaeans, and viruses suggests that it is a ubiquitous and phylogenetically old phenomenon. In animals, it is usually smallest for morphometric traits and highest for life history traits and behaviour. SDV is thought to contribute to phenotypic diversity in all populations but is particularly relevant for asexually reproducing and genetically impoverished populations, where it generates individuality despite genetic uniformity. In each generation, SDV produces a range of phenotypes around a well-adapted target phenotype, which is interpreted as a bet-hedging strategy to cope with the unpredictability of dynamic environments. At least some manifestations of SDV are heritable, adaptable, selectable, and evolvable, and therefore, SDV may be seen as a hitherto overlooked evolution factor. SDV is also relevant for husbandry, agriculture, and medicine because most pathogens are asexuals that exploit this third source of phenotypic variation to modify infectivity and resistance to antibiotics. Since SDV affects all types of organisms and almost all aspects of life, it urgently requires more intense research and a better integration into biological thinking.

Potential of agricultural fungicides for antifungal drug discovery

While it is true that only a small fraction of fungal species are responsible for human mycoses, the increasing prevalence of fungal diseases has highlighted an urgent need to develop new antifungal drugs, especially for systemic administration. This contribution focuses on the similarities between agricultural fungicides and drugs. Inorganic, organometallic and organic compounds can be found amongst agricultural fungicides. Furthermore, fungicides are designed and developed in a similar fashion to drugs based on similar rules and guidelines, with fungicides also having to meet similar criteria of lead-likeness and/or drug-likeness. Modern approved specific-target fungicides are well-characterized entities with a proposed structure-activity relationships hypothesis and a defined mode of action. Extensive toxicological evaluation, including mammalian toxicology assays, is performed during the whole discovery and development process. Thus modern agrochemical research (design of modern agrochemicals) comes close to drug design, discovery and development. Therefore, modern specific-target fungicides represent excellent lead-like structures/models for novel drug design and development.

Parasexual Ploidy Reduction Drives Population Heterogeneity Through Random and Transient Aneuploidy in Candida albicans

The opportunistic pathogen Candida albicans has a large repertoire of mechanisms to generate genetic and phenotypic diversity despite the lack of meiosis in its life cycle. Its parasexual cycle enables shifts in ploidy, which in turn facilitate recombination, aneuploidy, and homozygosis of whole chromosomes to fuel rapid adaptation. Here we show that the tetraploid state potentiates ploidy variation and drives population heterogeneity. In tetraploids, the rate of losing a single heterozygous marker [loss of heterozygosity (LOH)] is elevated ∼30-fold higher than the rate in diploid cells. Furthermore, isolates recovered after selection for LOH of one, two, or three markers were highly aneuploid, with a broad range of karyotypes including strains with a combination of di-, tri-, and tetrasomic chromosomes. We followed the ploidy trajectories for these tetraploid- and aneuploid-derived isolates, using a combination of flow cytometry and double-digestion restriction-site-associated DNA analyzed with next-generation sequencing. Isolates derived from either tetraploid or aneuploid isolates predominately resolved to a stable euploid state. The majority of isolates reduced to the conventional diploid state; however, stable triploid and tetraploid states were observed in ∼30% of the isolates. Notably, aneuploid isolates were more transient than tetraploid isolates, resolving to a euploid state within a few passages. Furthermore, the likelihood that a particular isolate will resolve to the same ploidy state in replicate evolution experiments is only ∼50%, supporting the idea that the chromosome loss process of the parasexual cycle is random and does not follow trajectories involving specific combinations of chromosomes. Together, our results indicate that tetraploid progenitors can produce populations of progeny cells with a high degree of genomic diversity, from altered ploidy to homozygosis, providing an excellent source of genetic variation upon which selection can act.

Limited Sexual Reproduction and Quick Turnover in the Population Genetic Structure of Phytophthora infestans in Fujian, China

The mating system plays an important role in the spatiotemporal dynamics of pathogen populations through both its direct and indirect impact on the generation and distribution of genetic variation. Here, we used a combination of microsatellite and phenotypic markers to investigate the spatiotemporal distribution of genetic variation in Phytophthora infestans isolates collected from Fujian, China and to determine the role of sexual reproduction in the dynamics. Although the pathogen populations in this region were dominated by self-fertile genotypes, sexual reproduction only occurred occasionally and its contributions to the population genetic structure of P. infestans and epidemics of late blight in the region were limited. Only 49 genotypes were detected among the 534 isolates assayed and the pathogen populations displayed significant heterozygosity excess. Hierarchical analysis revealed that 21.42% of genetic variation was attributed to the difference among sampling years while only 4.45% was attributed to the difference among locations, suggesting temporal factors play a more important role in the population genetic dynamics of P. infestans than spatial factors in this region. We propose that clonal reproduction, combined with founder effects and long distance dispersal of sporangia, is responsible for the observed pattern of spatiotemporal dynamics in P. infestans.

Genetic differentiation and recombination among geographic populations of the fungal pathogen Colletotrichum truncatum from chili peppers in China

Colletotrichum truncatum is an extremely important fungal pathogen. It can cause diseases both in humans and in over 460 plant species. However, little is known about its genetic diversity within and among populations. One of the major plant hosts of C. truncatum is pepper, and China is one of the main pepper-producing countries in the world. Here, we propose the hypotheses that geography has a major influence on the relationships among populations of C. truncatum in China and that infections in different populations need to be managed differently. To test these hypotheses, we obtained and analyzed 266 C. truncatum isolates from 13 regions representing the main pepper-growing areas throughout China. The analysis based on nine microsatellite markers identified high intrapopulation genetic diversity, evidence of sexual recombination, and geographic differentiation. The genetic differentiation was positively correlated with geographic distance, with the southern and northern China populations grouped in two distinct clusters. Interestingly, isolates collected from the pepper-breeding center harbored the most private alleles. The results suggest that the geographic populations of C. truncatum on peppers in China are genetically differentiated and should be managed accordingly. Our study also provides a solid foundation from which to further explore the global genetic epidemiology of C. truncatum in both plants and humans.

Fungal model systems and the elucidation of pathogenicity determinants

Fungi have the capacity to cause devastating diseases of both plants and animals, causing significant harvest losses that threaten food security and human mycoses with high mortality rates. As a consequence, there is a critical need to promote development of new antifungal drugs, which requires a comprehensive molecular knowledge of fungal pathogenesis. In this review, we critically evaluate current knowledge of seven fungal organisms used as major research models for fungal pathogenesis. These include pathogens of both animals and plants; Ashbya gossypii, Aspergillus fumigatus, Candida albicans, Fusarium oxysporum, Magnaporthe oryzae, Ustilago maydis and Zymoseptoria tritici. We present key insights into the virulence mechanisms deployed by each species and a comparative overview of key insights obtained from genomic analysis. We then consider current trends and future challenges associated with the study of fungal pathogenicity.

How often do they have sex? A comparative analysis of the population structure of seven eukaryotic microbial pathogens

The model of predominant clonal evolution (PCE) proposed for micropathogens does not state that genetic exchange is totally absent, but rather, that it is too rare to break the prevalent PCE pattern. However, the actual impact of this “residual” genetic exchange should be evaluated. Multilocus Sequence Typing (MLST) is an excellent tool to explore the problem. Here, we compared online available MLST datasets for seven eukaryotic microbial pathogens: Trypanosoma cruzi, the Fusarium solani complex, Aspergillus fumigatus, Blastocystis subtype 3, the Leishmania donovani complex, Candida albicans and Candida glabrata. We first analyzed phylogenetic relationships among genotypes within each dataset. Then, we examined different measures of branch support and incongruence among loci as signs of genetic structure and levels of past recombination. The analyses allow us to identify three types of genetic structure. The first was characterized by trees with well-supported branches and low levels of incongruence suggesting well-structured populations and PCE. This was the case for the T. cruzi and F. solani datasets. The second genetic structure, represented by Blastocystis spp., A. fumigatus and the L. donovani complex datasets, showed trees with weakly-supported branches but low levels of incongruence among loci, whereby genetic structuration was not clearly defined by MLST. Finally, trees showing weakly-supported branches and high levels of incongruence among loci were observed for Candida species, suggesting that genetic exchange has a higher evolutionary impact in these mainly clonal yeast species. Furthermore, simulations showed that MLST may fail to show right clustering in population datasets even in the absence of genetic exchange. In conclusion, these results make it possible to infer variable impacts of genetic exchange in populations of predominantly clonal micro-pathogens. Moreover, our results reveal different problems of MLST to determine the genetic structure in these organisms that should be considered.

Cryptosporidium,Giardia, Cryptococcus, Pneumocystis genetic variability: cryptic biological species or clonal near-clades?

An abundant literature dealing with the population genetics and taxonomy of Giardia duodenalis, Cryptosporidium spp., Pneumocystis spp., and Cryptococcus spp., pathogens of high medical and veterinary relevance, has been produced in recent years. We have analyzed these data in the light of new population genetic concepts dealing with predominant clonal evolution (PCE) recently proposed by us. In spite of the considerable phylogenetic diversity that exists among these pathogens, we have found striking similarities among them. The two main PCE features described by us, namely highly significant linkage disequilibrium and near-clading (stable phylogenetic clustering clouded by occasional recombination), are clearly observed in Cryptococcus and Giardia, and more limited indication of them is also present in Cryptosporidium and Pneumocystis. Moreover, in several cases, these features still obtain when the near-clades that subdivide the species are analyzed separately (“Russian doll pattern”). Lastly, several sets of data undermine the notion that certain microbes form clonal lineages simply owing to a lack of opportunity to outcross due to low transmission rates leading to lack of multiclonal infections (“starving sex hypothesis”). We propose that the divergent taxonomic and population genetic inferences advanced by various authors about these pathogens may not correspond to true evolutionary differences and could be, rather, the reflection of idiosyncratic practices among compartmentalized scientific communities. The PCE model provides an opportunity to revise the taxonomy and applied research dealing with these pathogens and others, such as viruses, bacteria, parasitic protozoa, and fungi.

Micropathogen species definition is extremely difficult, since concepts applied to higher organisms (the biological species concept) are inadequate. In particular, the pathogens here surveyed have given rise to long-lasting controversies about their species status and that of the genotypes that subdivide them. The population genetic approach based on the predominant clonal evolution (PCE) concept proposed by us could bring simple solutions to these controversies, since it permits the description of clearly defined evolutionary entities (clonal multilocus genotypes and near-clades [incompletely isolated clades]) that could be the basis for species description, if the concerned specialists find it justified for applied research. The PCE model also provides a convenient framework for applied studies (molecular epidemiology, vaccine and drug design, clinical research) dealing with these pathogens and others.

Origins of eukaryotic sexual reproduction

Sexual reproduction is a nearly universal feature of eukaryotic organisms. Given its ubiquity and shared core features, sex is thought to have arisen once in the last common ancestor to all eukaryotes. Using the perspectives of molecular genetics and cell biology, we consider documented and hypothetical scenarios for the instantiation and evolution of meiosis, fertilization, sex determination, uniparental inheritance of organelle genomes, and speciation.