Search

Ecological and genetic factors have influenced the composition of the human microbiome during our evolutionary history. We analyzed the oral microbiota of the Agta, a hunter-gatherer population where part of its members is adopting an agricultural diet. We show that age is the strongest factor modulating the microbiome, likely through immunosenescence as there is an increase of pathogenicity with age. Biological and cultural processes generate sexual dimorphism in the oral microbiome. A small subset of oral bacteria is influenced by the host genome, linking host collagen genes to bacterial biofilm formation. Our data also suggests that shifting from a fish/meat to a rice-rich diet transforms their microbiome, mirroring the Neolithic transition. All these factors have implications in the epidemiology of oral diseases. Thus, the human oral microbiome is multifactorial, and shaped by various ecological and social factors that modify the oral environment.

Chronic inflammation that affects primarily metabolic organs, such as white adipose tissue (WAT), is considered as a major cause of human obesity-associated co-morbidities. However, the molecular mechanisms initiating this inflammation in WAT are poorly understood. By combining transcriptomics, ChIP-seq and modeling approaches, we studied the global early and late responses to a high-fat diet (HFD) in visceral (vWAT) and subcutaneous (scWAT) AT, the first being more prone to obesity-induced inflammation. HFD rapidly triggers proliferation of adipocyte precursors within vWAT. However, concomitant antiadipogenic signals limit vWAT hyperplastic expansion by interfering with the differentiation of proliferating adipocyte precursors. Conversely, in scWAT, residing beige adipocytes lose their oxidizing properties and allow storage of excessive fatty acids. This phase is followed by tissue hyperplastic growth and increased angiogenic signals, which further enable scWAT expansion without generating inflammation. Our data indicate that scWAT and vWAT differential ability to modulate adipocyte number and differentiation in response to obesogenic stimuli has a crucial impact on the different susceptibility to obesity-related inflammation of these adipose tissue depots. Electronic supplementary material The online version of this article (10.1007/s00018-020-03485-z) contains supplementary material, which is available to authorized users.

HLA-G is a promiscuous immune checkpoint molecule. The HLA-G gene presents substantial nucleotide variability in its regulatory regions. However, it encodes a limited number of proteins compared to classical HLA class I genes. We characterized the HLA-G genetic variability in 4640 individuals from 88 different population samples across the globe by using a state-of-the-art method to characterize polymorphisms and haplotypes from high-coverage next-generation sequencing data. We also provide insights regarding the HLA-G genetic diversity and a resource for future studies evaluating HLA-G polymorphisms in different populations and association studies. Despite the great haplotype variability, we demonstrated that: (1) most of the HLA-G polymorphisms are in introns and regulatory sequences, and these are the sites with evidence of balancing selection, (2) linkage disequilibrium is high throughout the gene, extending up to HLA-A, (3) there are few proteins frequently observed in worldwide populations, with lack of variation in residues associated with major HLA-G biological properties (dimer formation, interaction with leukocyte receptors). These observations corroborate the role of HLA-G as an immune checkpoint molecule rather than as an antigen-presenting molecule. Understanding HLA-G variability across populations is relevant for disease association and functional studies., We acknowledge Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq/Brazil, #302060/2019-7), Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES/Brazil, Finance code: 001), INCTC/FAPESP (Grant #465539/2014-9) for financial support, and FAPESP/Brazil (#Grant 2017/19223-0).

The haploid ascomycete Diplocarpon rosae is the causal agent of black spot disease on roses, a widespread and devastating disease in the outdoor landscape. In this study, we established a Eurasian collection of 77 monoconidial strains of D. rosae: 50 strains collected on cultivated roses in Europe and Asia, and 27 strains on wild roses in Kazakhstan. To provide tools to describe its biology and to study its genetic diversity, we sequenced two strains of D. rosae using Illumina paired-end technology. The genome sizes of these two strains were estimated at 31.1 and 35.2 Mb, which are two times smaller than the genome size of the unique strain previously published. A BUSCO analysis confirmed a genome duplication of the strain previously sequenced and partial gene duplication of strains analysed in this study. Using the two genome sequences, 27 polymorphic microsatellite markers were identified. Polymorphism analysis of the 77 strains revealed a strong genetic differentiation between strains from cultivated and wild roses, and a lower diversity within the fungal population from cultivated roses compared to the population from wild roses. Pathogenicity of 10 strains was evaluated on 9 rose cultivars inoculated in the greenhouse. Disease scoring allowed the classification of strains into three groups and the characterization of resistance of rose cultivars. Good correlation observed between resistance scoring in greenhouse conditions and in the field indicates that pathogenicity assays in controlled conditions could be very useful in the near future to rapidly characterize the resistance of new rose varieties to black spot disease., This work was supported by the CASDAR no. C-2014-06 from the French Ministry of Agriculture, Agrifood and Forestry. The authors thank Markus Linde for kindly supplying German infected leaves, Célia Calmels for microsatellite analysis, Sandrine Pierre for the contribution to assays on whole plants, Phenotic and Experimental unit Horticole for plant management in the greenhouse and in the field, respectively, Thibault Leroy for reads assembly of the Kazakh strain, and Biogenouest for read assembly of the German strain. All microsatellite data were obtained on the ANAN platform (M. Bahut). This work was financially supported by SFR 4207 Quasav.

Identifying the genes underlying adaptation, their distribution in genomes, and the evolutionary forces shaping genomic diversity, are key challenges in evolutionary biology. Very few studies have investigated the abundance and distribution of selective sweeps in species with high-quality reference genomes, outside of a handful of model species. Pathogenic fungi are tractable eukaryote models for investigating the genomics of adaptation. By sequencing 53 genomes of two species of anther-smut fungi and mapping them against a high-quality reference genome, we showed that selective sweeps were abundant and scattered throughout the genome in one species, affecting near 17% of the genome, but much less numerous and in different genomic regions in its sister species, where they left footprints in only 1% of the genome. Polymorphism was negatively correlated with linkage disequilibrium levels in the genomes, consistent with recurrent positive and/or background selection. Differential expression in planta and in vitro, and functional annotation, suggested that many of the selective sweeps were probably involved in adaptation to the host plant. Examples include glycoside hydrolases, pectin lyases and an extracellular membrane protein with CFEM domain. This study thus provides candidates genes for being involved in plant - pathogen interaction (effectors), that have remained elusive for long in this otherwise well-studied system. Their identification will foster future functional and evolutionary studies, in the plant and in the anther-smut pathogens, being model species of natural plant-pathogen associations. In addition, our results suggest that positive selection can have a pervasive impact in shaping genomic variability in pathogens and selfing species, broadening our knowledge of the occurrence and frequency of selective events in natural populations. This article is protected by copyright. All rights reserved.

While the extent and impact of horizontal transfers in prokaryotes are widely acknowledged, their importance to the eukaryotic kingdom is unclear and thought by many to be anecdotal. Here we report multiple recent transfers of a huge genomic island between Penicillium spp. found in the food environment. Sequencing of the two leading filamentous fungi used in cheese making, P. roqueforti and P. camemberti, and comparison with the penicillin producer P. rubens reveals a 575 kb long genomic island in P. roqueforti—called Wallaby—present as identical fragments at non-homologous loci in P. camemberti and P. rubens. Wallaby is detected in Penicillium collections exclusively in strains from food environments. Wallaby encompasses about 250 predicted genes, some of which are probably involved in competition with microorganisms. The occurrence of multiple recent eukaryotic transfers in the food environment provides strong evidence for the importance of this understudied and probably underestimated phenomenon in eukaryotes., Horizontal gene transfers are known to play an important role in prokaryote evolution but their impact and prevalence in eukaryotes is less clear. Here, the authors sequence the genomes of cheese making fungi P. roqueforti and P. camemberti, and provide evidence for recent horizontal transfers of a large genomic region.

About 100,000 species of fungi have been described so far, of which a high percentage obtain nutrients by living in close association with other organisms, mainly plants. Many fungi are pathogenic and can have important impact on human health or lead to severe economic losses due to infected crops or to animal diseases. In this chapter, we describe the main pathogenic fungi, parasitizing humans, animals, and plants, and having important consequences on human health or human activities. We focus on some examples of emerging fungal diseases, of humans, animals, and plants. We then review (1) the modern molecular tools used for epidemiology and population genetics of fungal pathogens, the types of markers most useful, and the different types of analyses that can be performed to unravel their mating systems and dispersal; and (2) the advances since 2000s in fungal genomics, in particular the insights that have been gained so far regarding the pathogenic lifestyles.

Parallels have been drawn between the evolution of nonrecombining regions in fungal mating-type chromosomes and animal and plant sex chromosomes, particularly regarding the stages of recombination cessation forming evolutionary strata of allelic divergence. Currently, evidence and explanations for recombination cessation in fungi are sparse, and the presence of evolutionary strata has been examined in a minimal number of fungal taxa. Here, the basidiomycete genus Microbotryum was used to determine the history of recombination cessation for loci on the mating-type chromosomes. Ancestry of linkage with mating type for 13 loci was assessed across 20 species by a phylogenetic method. No locus was found to exhibit trans-specific polymorphism for alternate alleles as old as the mating pheromone receptor, indicating that ages of linkage to mating type varied among the loci. The ordering of loci in the ancestry of linkage to mating type does not agree with their previously proposed assignments to evolutionary strata. This study suggests that processes capable of influencing divergence between alternate alleles may act at loci in the nonrecombining regions (e.g., gene conversion) and encourages further work to dissect the evolutionary processes acting upon genomic regions that determine mating compatibility.

Phylogenies are fundamental to comparative biology as they help to identify independent events on which statistical tests rely. Two groups of phylogenetic comparative methods (PCMs) can be distinguished: those that take phy- logenies into account by introducing explicit models of evolution and those that only consider phylogenies as a statistical constraint and aim at partitioning trait values into a phylogenetic component (phylogenetic inertia) and one or multiple specific components related to adaptive evolution. The way phylogenetic information is incorporated into the PCMs de- pends on the method used. For the first group of methods, phylogenies are converted into variance-covariance matrices of traits following a given model of evolution such as Brownian motion (BM). For the second group of methods, phyloge- nies are converted into distance matrices that are subsequently transformed into Euclidean distances to perform principal coordinate analyses. Here, we show that simply taking the elementwise square root of a distance matrix extracted from a phylogenetic tree ensures having a Euclidean distance matrix. This is true for any type of distances between species (patristic or nodal) and also for trees harboring multifurcating nodes. Moreover, we illustrate that this simple transforma- tion using the square root imposes less geometric distortion than more complex transformations classically used in the literature such as the Cailliez method. Given the Euclidean nature of the elementwise square root of phylogenetic distance matrices, the positive semidefinitiveness of the phylogenetic variance-covariance matrix of a trait following a BM model, or related models of trait evolution, can be established. In that way, we build a bridge between the two groups of statistical methods widely used in comparative analysis. These results should be of great interest for ecologists and evolutionary biologists performing statistical analyses incorporating phylogenies. (Biodiversity estimation; coevolution; cospeciation; functional diversity; multivariate analysis; phylogeny.)

Several cases of speciation by hybridization have been reported in fungi, mostly involving recent hybridization between closely related species. In the basidiomycete genus Microbotryum by contrast some species were suspected to have arisen by hybridization between moderately distant species. In particular two species, M. lagerheimii and M. silenes-acaulis, had different placements in phylogenetic trees depending on the genes considered. Microbotryum species exhibit bipolar heterothallism, and here we analyzed se- quences of the two alternate pheromone receptors to obtain further insights on the occurrence of hybridization. Indeed because mating-type loci are always heterozygous homoploid hybrid speciation should leave a permanent footprint at the mating- type locus by retaining the alternate alleles from their respective parental species. The trees obtained with each of the two pheromone receptors were well resolved, and the species relationships were in agreement with published phylogenies. Fungal pher- omone receptor genes of basidiomycetes thus appear useful for phylogenetic studies, although it may not be true for the homobasidiomycetes where duplica- tions of these genes have occurred. Furthermore an incongruence between the phylogenies of the two pheromone receptors was found for one species, M. lagerheimii, as previously observed between other nuclear genes. However additional species analyzed here revealed that the incongruence involved the whole clade including both M. lagerheimii and the Microbotryum species parasitizing Lychnis flos-cucucli. The ancestor of these species thus possibly arose via hybridization between distant ancestral lineages, although further studies should address alternative hypotheses, such as chance events during lineage sorting.

Numerous genes in diverse organisms have been shown to be under positive selection, especially genes involved in reproduction, adaptation to contrasting environments, hybrid inviability, and host-pathogen interactions. Looking for genes under positive selection in pathogens has been a priority in efforts to investigate coevolution dynamics and to develop vaccines or drugs. To elucidate the functions involved in host specialization, here we aimed at identifying candidate sequences that could have evolved under positive selection among closely related pathogens specialized on different hosts. For this goal, we sequenced c. 17,000-32,000 ESTs from each of four Microbotryum species, which are fungal pathogens responsible for anther smut disease on host plants in the Caryophyllaceae. Forty-two of the 372 predicted orthologous genes showed significant signal of positive selection, which represents a good number of candidate genes for further investigation. Sequencing 16 of these genes in 9 additional Microbotryum species confirmed that they have indeed been rapidly evolving in the pathogen species specialized on different hosts. The genes showing significant signals of positive selection were putatively involved in nutrient uptake from the host, secondary metabolite synthesis and secretion, respiration under stressful conditions and stress response, hyphal growth and differentiation, and regulation of expression by other genes. Many of these genes had transmembrane domains and may therefore also be involved in pathogen recognition by the host. Our approach thus revealed fruitful and should be feasible for many non-model organisms for which candidate genes for diversifying selection are needed.

International audience; Nuclear disasters at Chernobyl and Fukushima provide examples of effects of acute ionizing radiation on mutations that can affect the fitness and distribution of species. Here, we investigated the prevalence of Microbotryum lychnidis-dioicae, a pollinator-transmitted fungal pathogen of plants causing anther-smut disease in Chernobyl, its viability, fertility and karyotype variation, and the accumulation of nonsynonymous mutations in its genome. We collected diseased flowers of Silene latifolia from locations ranging by more than two orders of magnitude in background radiation, from 0.05 to 21.03 lGy/h. Disease prevalence decreased significantly with increasing radiation level, possibly due to lower pollinator abundance and altered pollinator behaviour. Viability and fertility, measured as the budding rate of haploid sporidia following meiosis from the diploid teliospores, did not vary with increasing radiation levels and neither did karyotype overall structure and level of chromosomal size heterozygosity. We sequenced the genomes of twelve samples from Chernobyl and of four samples collected from uncontaminated areas and analysed alignments of 6068 predicted genes, corresponding to 1.04 3 10 7 base pairs. We found no dose-dependent differences in substitution rates (neither dN, dS, nor dN/dS). Thus, we found no significant evidence of increased deleterious mutation rates at higher levels of background radiation in this plant pathogen. We even found lower levels of nonsynonymous substitution rates in contaminated areas compared to control regions, suggesting that purifying selection was stronger in contaminated than uncontaminated areas. We briefly discuss the possibilities for a mechanistic basis of radio resistance in this nonmelanized fungus.

A comprehensive dataset of 62 beta globin gene sequences from various vertebrates was compiled to test the molecular clock and to estimate dates of gene duplications. We found that evolution of the beta globin family of genes is not clock-like, a result that is at odds with the common use of this family as an example of a constant rate of evolution over time. Divergence dates were estimated either with or without assuming the molecular clock, and both analyses produced similar date estimates, which are also in general agreement with estimates reported previously. In addition we report date estimates for seven previously unexamined duplication events within the beta globin family. Despite multiple sources of rate variation, the average rate across the beta globin phylogeny yielded reasonable estimates of divergence dates in most cases. Exceptions were cases of gene conversion, where it appears to have led to underestimates of divergence dates. Our results suggest (i) the major duplications giving rise to the paralogous beta globin genes are associated with significant evolutionary rate variation among gene lineages; and (ii) genes arising from more recent gene duplications (e.g., tandem duplications within lineages) do not appear to differ greatly in rate. We believe this pattern reflects a complex interplay of evolutionary forces where natural selection for diversifying paralogous functions and lineage-specific effects contribute to rate variation on a long-term basis, while gene conversion tends to increase sequence similarity. Gene conversion effects appear to be stronger on recent gene duplicates, as their sequences are highly similar. Lastly, phylogenetic analyses do not support a previous report that avian globins are members of a relic lineage of omega globins.

Different models of gene family evolution have been proposed to explain the mechanism whereby gene copies created by gene duplications are maintained and diverge in function. Ohta proposed a model which predicts a burst of nonsynonymous substitutions following gene duplication and the preservation of duplicates through positive selection. An alternative model, the duplication-degeneration-complementation (DDC) model, does not explicitly require the action of positive Darwinian selection for the maintenance of duplicated gene copies, although purifying selection is assumed to continue to act on both copies. A potential outcome of the DDC model is heterogeneity in purifying selection among the gene copies, due to partitioning of subfunctions which complement each other. By using the d(N)/ d(S) (omega) rate ratio to measure selection pressure, we can distinguish between these two very different evolutionary scenarios. In this study we investigated these scenarios in the beta-globin family of genes, a textbook example of evolution by gene duplication. We assembled a comprehensive dataset of 72 vertebrate beta-globin sequences. The estimated phylogeny suggested multiple gene duplication and gene conversion events. By using different programs to detect recombination, we confirmed several cases of gene conversion and detected two new cases. We tested evolutionary scenarios derived from Ohta's model and the DDC model by examining selective pressures along lineages in a phylogeny of beta-globin genes in eutherian mammals. We did not find significant evidence for an increase in the omega ratio following major duplication events in this family. However, one exception to this pattern was the duplication of gamma-globin in simian primates, after which a few sites were identified to be under positive selection. Overall, our results suggest that following gene duplications, paralogous copies of beta-globin genes evolved under a nonepisodic process of functional divergence.

© 2015 Perlin et al. This is an Open Access article distributed under the terms of the Creative Commons Attribution License/n(http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium,/nprovided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http:///ncreativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated. Background: The genus Microbotryum includes plant pathogenic fungi afflicting a wide variety of hosts with anther smut disease. Microbotryum lychnidis-dioicae infects Silene latifolia and replaces host pollen with fungal spores, exhibiting biotrophy and necrosis associated with altering plant development. Results: We determined the haploid genome sequence for M. lychnidis-dioicae and analyzed whole transcriptome data from plant infections and other stages of the fungal lifecycle, revealing the inventory and expression level of genes that facilitate pathogenic growth. Compared to related fungi, an expanded number of major facilitator superfamily transporters and secretory lipases were detected; lipase gene expression was found to be altered by exposure to lipid compounds, which signaled a switch to dikaryotic, pathogenic growth. In addition, while enzymes to digest cellulose, xylan, xyloglucan, and highly substituted forms of pectin were absent, along with depletion of peroxidases and superoxide dismutases that protect the fungus from oxidative stress, the repertoire of glycosyltransferases and of enzymes that could manipulate host development has expanded. A total of 14 % of the genome was categorized as repetitive sequences. Transposable elements have accumulated in mating-type chromosomal regions and were also associated across the genome with gene clusters of small secreted proteins, which may mediate host interactions. Conclusions: The unique absence of enzyme classes for plant cell wall degradation and maintenance of enzymes that break down components of pollen tubes and flowers provides a striking example of biotrophic host adaptation.

International audience; Horizontal gene transfers (HGT), i.e., the transmission of genetic material between species not directly attributable to meiotic gene exchange, have long been acknowledged as a major driver of prokaryotic evolution and is increasingly recognized as an important source of adaptation in eukaryotes. In fungi in particular , many convincing examples of HGT have been reported to confer selective advantages on the recipient fungal host, either promoting fungal pathogenicity on plants or increasing their toxicity by the acquisition of secondary metabolic clusters, resulting in adaptation to new niches and in some cases eventually even in speciation. These horizontal gene transfers involve single genes, complete metabolic pathways or even entire chromosomes. A recent study has uncovered multiple recent horizontal transfers of a 575 kb ge-nomic island in cheese Penicillium fungi, representing ca. 2% of the Penicillium roqueforti's genome, that may confer selective advantage in the competing cheese environment where bacteria and fungi occur. Novel phylogenomic methods are being developed, revealing massive HGT among fungi. Altogether, these recent studies indicate that HGT is a crucial mechanism of rapid adaptation, even among eukaryotes. A well-documented example of pathogenicity acquisition in fungi involves the recent transfer of the toxin-coding gene ToxA, from Stagonospora nodorum, a fungus pathogen of wheat, to Pyrenophora tritici-repentis, a distant fungal species having thereby acquired the ability to infect wheat. Pyrenophora has also acquired numerous virulence factors

BACKGROUND: Candida glabrata follows C. albicans as the second or third most prevalent cause of candidemia worldwide. These two pathogenic yeasts are distantly related, C. glabrata being part of the Nakaseomyces, a group more closely related to Saccharomyces cerevisiae. Although C. glabrata was thought to be the only pathogenic Nakaseomyces, two new pathogens have recently been described within this group: C. nivariensis and C. bracarensis. To gain insight into the genomic changes underlying the emergence of virulence, we sequenced the genomes of these two, and three other non-pathogenic Nakaseomyces, and compared them to other sequenced yeasts. RESULTS: Our results indicate that the two new pathogens are more closely related to the non-pathogenic N. delphensis than to C. glabrata. We uncover duplications and accelerated evolution that specifically affected genes in the lineage preceding the group containing N. delphensis and the three pathogens, which may provide clues to the higher propensity of this group to infect humans. Finally, the number of Epa-like adhesins is specifically enriched in the pathogens, particularly in C. glabrata. CONCLUSIONS: Remarkably, some features thought to be the result of adaptation of C. glabrata to a pathogenic lifestyle, are present throughout the Nakaseomyces, indicating these are rather ancient adaptations to other environments. Phylogeny suggests that human pathogenesis evolved several times, independently within the clade. The expansion of the EPA gene family in pathogens establishes an evolutionary link between adhesion and virulence phenotypes. Our analyses thus shed light onto the relationships between virulence and the recent genomic changes that occurred within the Nakaseomyces. This work was supported by funding from the European Research Council/nunder the European Union’s Seventh Framework Programme (FP/2007- 2013)/ERC Grant Agreement n.310325, a Grant from the Qatar National Research Fund grant (NPRP 5-298-3-086), and by a grant from the Spanish Ministry of Economy and Competitiveness (BIO2012-37161). CF’s research is/nfunded in part by an “Attractivité” grant from the University Paris Sud. GA is/na recipient of a Marie Curie grant (FP7-PEOPLE-2010-IEF-No.274223). SB, HD and RA are recipients of, respectively, a shared post-doctoral grant and a Ph. D. grant, from the Région Ile-de-France’s DIM Malinf program. JAC was supported by the Ph.D. Program in Computational Biology of the Instituto Gulbenkian de Ciência, Portugal (sponsored by Fundação Calouste Gulbenkian, Siemens SA, and Fundação para a Ciência e Tecnologia; SFRH/BD/33528/2008)

International audience; Fungi are ideal model organisms for dissecting the genomic bases of adaptive divergence in eukaryotes. They have simple morphologies and small genomes, occupy contrasting, well-identified ecological niches and tend to have short generation times, and many are amenable to experimental approaches. Fungi also display diverse lifestyles , from saprotrophs to pathogens or mutualists, and they play extremely important roles in both ecosystems and human activities, as wood decayers, mycorrhizal fungi, lichens, endophytes, plant and animal pathogens, and in fermentation or drug production. We review here recent insights into the patterns and mechanisms of adaptive divergence in fungi, including sources of divergence, genomic variation and, ultimately , speciation. We outline the various ecological sources of divergent selection and genomic changes, showing that gene loss and changes in gene expression and in geno-mic architecture are important adaptation processes, in addition to the more widely recognized processes of amino acid substitution and gene duplication. We also review recent findings regarding the interspecific acquisition of genomic variation and suggesting an important role for introgression, hybridization and horizontal gene transfers (HGTs). We show that transposable elements can mediate several of these genomic changes, thus constituting important factors for adaptation. Finally, we review the consequences of divergent selection in terms of speciation, arguing that genetic incompatibilities may not be as widespread as generally thought and that pleiotropy between adaptation and reproductive isolation is an important route of speciation in fungal pathogens.

International audience; Tests of selective neutrality based on comparative and population genetic data target different time scales of selection. Investigating both time scales on the same genes has rarely been done and yet can provide insights into histories of selection. In pathogens, such a comparison can elucidate whether the same genes experience recurrent positive selection across specialization events onto novel hosts and selection in the short term, e.g. by coevolution with a current host. Here, we investigated whether the genes showing signs of positive selection (based on d(N)/d(S) ratios) between pathogens specialized on different hosts also exhibited footprints of recent positive selection. We sequenced 58 DNA fragments in the model fungal plant pathogen Microbotryum lychnidis-dioicae infecting Silene latifolia. Eleven focal genes were chosen because they showed a significant signal of positive selection in a previous study and putative functions likely involved in host-parasite interactions. In addition, 47 control genes were randomly chosen among genes showing no evidence of positive selection. We used a population genetics approach to search for signatures of recent selection, controlling for population structure and demographic history. We found footprints of purifying selection in all focal genes and of recent positive selection in two of them. Signs of purifying selection were also found at the same eleven genes in a closely related species, Microbotryum silenes-dioicae, infecting Silene dioica. These results suggest that genes experiencing episodes of adaptive diversification during host shifts may subsequently be under strong functional constraint, although some can remain under positive selection in the short term.

The rapid evolution of particular genes is essential for the adaptation of pathogens to new hosts and new environments. Powerful methods have been developed for detecting targets of selection in the genome. Here we used divergence data to compare genes among four closely related fungal pathogens adapted to different hosts to elucidate the functions putatively involved in adaptive processes. For this goal, ESTs were sequenced in the specialist fungal pathogens Botrytis tulipae and Botrytis ficariarum, and compared with genome sequences of Botrytis cinerea and Sclerotinia sclerotiorum, responsible for diseases on over 200 plant species. A maximum likelihood-based analysis of 642 predicted orthologs detected 21 genes showing footprints of positive selection. These results were validated by resequencing nine of these genes in additional Botrytis species, showing they have also been rapidly evolving in other related species. Twenty of the 21 genes had not previously been identified as pathogenicity factors in B. cinerea, but some had functions related to plant-fungus interactions. The putative functions were involved in respiratory and energy metabolism, protein and RNA metabolism, signal transduction or virulence, similarly to what was detected in previous studies using the same approach in other pathogens. Mutants of B. cinerea were generated for four of these genes as a first attempt to elucidate their functions.

Full genome data sets are currently being explored on a regular basis to infer phylogenetic trees, but there are often discordances among the trees produced by different genes. An important goal in phylogenomics is to identify which individual gene and species produce the same phylogenetic tree and are thus likely to share the same evolutionary history. On the other hand, it is also essential to identify which genes and species produce discordant topologies and therefore evolve in a different way or represent noise in the data. The latter are outlier genes or species and they can provide a wealth of information on potentially interesting biological processes, such as incomplete lineage sorting, hybridization, and horizontal gene transfers. Here, we propose a new method to explore the genomic tree space and detect outlier genes and species based on multiple co-inertia analysis (MCOA), which efficiently captures and compares the similarities in the phylogenetic topologies produced by individual genes. Our method allows the rapid identification of outlier genes and species by extracting the similarities and discrepancies, in terms of the pairwise distances, between all the species in all the trees, simultaneously. This is achieved by using MCOA, which finds successive decomposition axes from individual ordinations (i.e., derived from distance matrices) that maximize a covariance function. The method is freely available as a set of R functions. The source code and tutorial can be found online at http://phylomcoa.cgenomics.org.

Publisher Summary This chapter discusses the main pathogenic fungi parasitizing humans, animals, and plants, and having important consequences on human health or human activities. It discusses the modern molecular tools used for epidemiology and population genetics of fungal pathogens. The two major groups, Ascomycota and Basidiomycota, are the most common causes of invasive fungal infections in humans. Fungal pathogens can cause severe disease in wildlife species and plants. They reduce crop yield and lower product quality by attacking cultivated plants and their products. Our primary food production is at risk due to emerging crop diseases. Epidemics caused by invasive pathogens have been repeatedly reported to alter natural ecosystems. Molecular epidemiologic studies to study fungal pathogens use techniques such as MLST, multilocus strain typing and MLMT, multilocus microsatellite typing. Gene-sequence data from studies using genealogical concordance phylogenetic recognition (GCPSR) can easily be converted to strain sequence types (STs) and stored in web databases. Population genetics is also needed to understand fungal diseases. It has provided important insights for some fungal pathogens on their mating systems, dispersal, and population structure. It offers the opportunity to forecast the emergence of genotypes, populations, or species with detrimental characteristics for human affairs. This chapter reveals that comparative genomic studies in plant pathogenic and symbiotic fungi have brought many insights into the evolution of the pathogenic lifestyle, in particular into the mechanisms of virulence and host adaptations. Regarding epidemiology, it is concluded that molecular methods have much to offer to the study of fungal pathogens, allowing elucidation of ecological and microevolutionary processes. The chapter states that microsatellite markers in particular are very powerful tools and should be more widely used for population studies in fungi, despite the technical challenges of their isolation in this kingdom.

UMR BGPI Equipe 5; International audience; We report the development of 17 microsatellite markers in the cheese fungi Penicillium camemberti and P. roqueforti, using an enrichment protocol. Polymorphism and cross-amplification were explored using 23 isolates of P. camemberti, 26 isolates of P. roqueforti and 2 isolates of each of the A chtysogenum and P. nalgiovense species. used to produce meat fermented products. The markers appeared useful for differentiating species, both using their amplification sizes and the sequences of their flanking regions. The microsatellite locus PC4 was particularly suitable for distinguishing contaminant species closely related to P. camemberti and for clarifying the phylogenetic relationship of this species with its supposed ancestral form, P. commune. We analyzed 22 isolates from different culture collections assigned to the morphospecies P commune, most of them occurring as food spoilers, mainly from the cheese environment. None of them exhibited identical sequences with the ex-type isolate of the species P. commune. They were instead distributed into two other distinct lineages, corresponding to the old species P. fuscoglaucum and P. biforme, previously synonymised respectively with P. commune and P. camemberti. The ex-type isolate of P. commune was strictly identical to P. camemberti at all the loci examined. P. caseifulvum, a non toxinogenic species described as a new candidate for cheese fermentation, also exhibited sequences identical to P. camemberti. The microsatellite locus PC4 may therefore be considered as a useful candidate for the barcode of these economically important species

We studied the mode of reproduction and its evolution in the fungal subgenus Penicillium Biverticillium using phylogenetic and experimental approaches. We sequenced mating type (MAT) genes and nuclear DNA fragments in sexual and putatively asexual species. Examination of the concordance between individual trees supported the recognition of the morphological species. MAT genes were detected in two putatively asexual species and were found to evolve mostly under purifying selection, although high substitution rates were detected at some sites in some clades. The first steps of sexual reproduction could be induced under controlled conditions in one of the two species, although no mature cleistothecia were produced. Altogether, these findings suggest that the asexual Penicillium species may have lost sex only very recently and/or that the MAT genes are involved in other functions. An ancestral state reconstruction analysis indicated several events of putative sex loss in the genus. Alternatively, it is possible that the supposedly asexual Penicillium species may have retained a cryptic sexual stage.

BGPI : équipe 5; International audience; The ongoing coevolutionary struggle between hosts and pathogens, with hosts evolving to escape pathogen infection and pathogens evolving to escape host defences, can generate an ‘arms race’, i.e., the occurrence of recurrent selective sweeps that each favours a novel resistance or virulence allele that goes to fixation. Host–pathogen coevolution can alternatively lead to a ‘trench warfare’, i.e., balancing selection, maintaining certain alleles at loci involved in host–pathogen recognition over long time scales. Recently, technological and methodological progress has enabled detection of footprints of selection directly on genes, which can provide useful insights into the processes of coevolution. This knowledge can also have practical applications, for instance development of vaccines or drugs. Here we review the methods for detecting genes under positive selection using divergence data (i.e., the ratio of nonsynonymous to synonymous substitution rates, dN/dS). We also review methods for detecting selection using polymorphisms, such as methods based on FST measures, frequency spectrum, linkage disequilibrium and haplotype structure. In the second part, we review examples where targets of selection have been identified in pathogens using these tests. Genes under positive selection in pathogens have mostly been sought among viruses, bacteria and protists, because of their paramount importance for human health. Another focus is on fungal pathogens owing to their agronomic importance. We finally discuss promising directions in pathogen studies, such as detecting selection in non-coding regions.

In the majority of sexual organisms, reproduction occurs almost exclusively through the combination of distinct and alternate forms, called sexes or mating types. In some fungi, there can be dozens to hundreds of alternate alleles that determine compatible mating types. Such extensive polymorphism is expected to be maintained by balancing selection, and in extreme cases may give rise to trans-specific polymorphism. Here, we analyzed sequences of two pheromone receptors in the Microbotryum fungal species complex (Basidiomycota), which has only two alternate mating types. Several lines of evidence strongly suggest that the pheromone receptors are two allelic sequences acting to determine the alternate A1 and A2 mating types required for mating in Microbotryum. Phylogenetic trees of pheromone receptors in the Microbotryum species complex indicated a trans-specific polymorphism: the Microbotryum sequences from a given mating type were all more similar to the pheromone receptors of distantly related classes of fungi than to the alternate pheromone receptor in the Microbotryum species. A phylogenetic tree built using other known pheromone receptors from basidiomycetes showed that trans-specific polymorphism is widespread. The pheromone receptor alleles from Microbotryum appeared as the oldest, being at least 370 million years old. This represents the oldest known trans-specific polymorphism known in any organism so far, which may be due to the existence of sex chromosomes, obligate sexuality, mitochondrial inheritance linked to the mating type, and a highly selfing mating system in Microbotryum.

Approximately 100,000 species of fungi have been described so far, of which, a high percentage obtain nutrients by living in close association with other organisms (mainly plants), being pathogens or symbionts (i.e. commensalists or mutualists). At the genomic level, an association between broad-scale genetic changes and the emergences of the parasitic and symbiotic life style in fungi has been proposed. Although comparative genomic studies in fungi are still in the early stages, they have enormous potential to improve our understanding of how such life styles evolve. In this chapter, we review the main characteristics of genome evolution in fungi, particularly in species that are pathogenic or symbiotic with plants, by focusing on the mechanisms involved in host interactions. We address the following topics in relation to the pathogenic and symbiotic lifestyles in fungi: the evolution of genome organization, chromosomal rearrangements, evolution of gene families and clusters, suppression of recombination around mating type loci, rapidly evolving genes, horizontal transfer, hybridization, transposable elements, telomeres, introns, mitochondrial genomes, and we finally compare genome evolution between pathogenic bacteria and fungi.

International audience; Phylogenies involving nonmodel species are based on a few genes, mostly chosen following historical or practical criteria. Because gene trees are sometimes incongruent with species trees, the resulting phylogenies may not accurately reflect the evolutionary relationships among species. The increase in availability of genome sequences now provides large numbers of genes that could be used for building phylogenies. However, for practical reasons only a few genes can be sequenced for a wide range of species. Here we asked whether we can identify a few genes, among the single-copy genes common to most fungal genomes, that are sufficient for recovering accurate and well-supported phylogenies. Fungi represent a model group for phylogenomics because many complete fungal genomes are available. An automated procedure was developed to extract single-copy orthologous genes from complete fungal genomes using a Markov Clustering Algorithm (Tribe-MCL). Using 21 complete, publicly available fungal genomes with reliable protein predictions, 246 single-copy orthologous gene clusters were identified. We inferred the maximum likelihood trees using the individual orthologous sequences and constructed a reference tree from concatenated protein alignments. The topologies of the individual gene trees were compared to that of the reference tree using three different methods. The performance of individual genes in recovering the reference tree was highly variable. Gene size and the number of variable sites were highly correlated and significantly affected the performance of the genes, but the average substitution rate did not. Two genes recovered exactly the same topology as the reference tree, and when concatenated provided high bootstrap values. The genes typically used for fungal phylogenies did not perform well, which suggests that current fungal phylogenies based on these genes may not accurately reflect the evolutionary relationships among species. Analyses on subsets of species showed that the phylogenetic performance did not seem to depend strongly on the sample. We expect that the best-performing genes identified here will be very useful for phylogenetic studies of fungi, at least at a large taxonomic scale. Furthermore, we compare the method developed here for finding genes for building robust phylogenies with previous ones and we advocate that our method could be applied to other groups of organisms when more complete genomes are available.

Background The increasing availability of fungal genome sequences provides large numbers of proteins for evolutionary and phylogenetic analyses. However the heterogeneity of data, including the quality of genome annotation and the difficulty of retrieving true orthologs, makes such investigations challenging. The aim of this study was to provide a reliable and integrated resource of orthologous gene families to perform comparative and phylogenetic analyses in fungi. Description FUNYBASE is a database dedicated to the analysis of fungal single-copy genes extracted from available fungal genomes sequences, their classification into reliable clusters of orthologs, and the assessment of their informative value for phylogenetic reconstruction based on amino acid sequences. The current release of FUNYBASE contains two types of protein data: (i) a complete set of protein sequences extracted from 30 public fungal genomes and classified into clusters of orthologs using a robust automated procedure, and (ii) a subset of 246 reliable ortholog clusters present as single copy genes in 21 fungal genomes. For each of these 246 ortholog clusters, phylogenetic trees were reconstructed based on their amino acid sequences. To assess the informative value of each ortholog cluster, each was compared to a reference species tree constructed using a concatenation of roughly half of the 246 sequences that are best approximated by the WAG evolutionary model. The orthologs were classified according to a topological score, which measures their ability to recover the same topology as the reference species tree. The full results of these analyses are available on-line with a user-friendly interface that allows for searches to be performed by species name, the ortholog cluster, various keywords, or using the BLAST algorithm. Examples of fruitful utilization of FUNYBASE for investigation of fungal phylogenetics are also presented. Conclusion FUNYBASE constitutes a novel and useful resource for two types of analyses: (i) comparative studies can be greatly facilitated by reliable clusters of orthologs across sets of user-defined fungal genomes, and (ii) phylogenetic reconstruction can be improved by identifying genes with the highest informative value at the desired taxonomic level.

The globin family of genes and proteins has been a recurrent object of study for many decades. This interest has generated a vast amount of knowledge. However it has also created an inconsistent and confusing nomenclature, due to the lack of a systematic approach to naming genes and failure to reflect the phylogenetic relationships among genes of the gene family. To alleviate the problems with the existing system, here we propose a standardized nomenclature for the alpha and beta globin family of genes, based on a phylogenetic analysis of vertebrate alpha and beta globins, and following the Guidelines for Human Gene Nomenclature.

BGPI : équipe 4 et 5; Deciphering the genetic bases of pathogen adaptation to its host is a key question in ecology and evolution. To understand how the fungus Magnaporthe oryzae adapts to different plants, we sequenced eight M. oryzae isolates differing in host specificity (rice, foxtail millet, wheat and goosegrass), and one Magnaporthe grisea isolate specific of crabgrass. Analysis of Magnaporthe genomes revealed small variation in genome sizes [39-43 Mb] and gene content [12,283-14,781 genes] between isolates. The whole set of Magnaporthe genes comprised 14,966 shared families, 63% of which included genes present in all the 9 M. oryzae genomes. The evolutionary relationships among Magnaporthe isolates were inferred using 6,878 single-copy orthologs. The resulting genealogy was mostly bifurcating among the different host-specific lineages, but reticulate inside the rice lineage. We detected traces of introgression from a non-rice genome in the rice reference 70-15 genome. Among M. oryzae isolates and host specific lineages, the genome composition in terms of frequencies of genes putatively involved in pathogenicity (effectors, secondary metabolism, cazome) was conserved. However, 529 shared families were found only in non-rice lineages, whereas the rice lineage possessed 86 specific families absent from the non-rice genomes. Our results confirmed that the host specificity of M. oryzae isolates was associated with a divergence between lineages without major gene flow and that, despite the strong conservation of gene families between lineages, adaptation to different hosts, especially to rice, was associated with the presence of a small number of specific gene families. All information was gathered in a public database (http://genome.jouy.inra.fr/gemo).