Gabriela Aguileta

Assessing the Performance of Single-Copy Genes for Recovering Robust Phylogenies

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.

Comparative genomics of emerging pathogens in the Candida glabrata clade

BackgroundCandida 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.ResultsOur 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.ConclusionsRemarkably, 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.Sequence Accession NumbersNakaseomyces delphensis: CAPT01000001 to CAPT01000179Candida bracarensis: CAPU01000001 to CAPU01000251Candida nivariensis: CAPV01000001 to CAPV01000123Candida castellii: CAPW01000001 to CAPW01000101Nakaseomyces bacillisporus: CAPX01000001 to CAPX01000186

Fungal evolutionary genomics provides insight into the mechanisms of adaptive divergence in eukaryotes

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 genomic 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.

FUNYBASE: a FUNgal phYlogenomic dataBASE

BackgroundThe 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.DescriptionFUNYBASE 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.ConclusionFUNYBASE 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.

Ancient Trans-Specific Polymorphism at Pheromone Receptor Genes in Basidiomycetes

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.

High Variability of Mitochondrial Gene Order among Fungi

From their origin as an early alpha proteobacterial endosymbiont to their current state as cellular organelles, large-scale genomic reorganization has taken place in the mitochondria of all main eukaryotic lineages. So far, most studies have focused on plant and animal mitochondrial (mt) genomes (mtDNA), but fungi provide new opportunities to study highly differentiated mtDNAs. Here, we analyzed 38 complete fungal mt genomes to investigate the evolution of mtDNA gene order among fungi. In particular, we looked for evidence of nonhomologous intrachromosomal recombination and investigated the dynamics of gene rearrangements. We investigated the effect that introns, intronic open reading frames (ORFs), and repeats may have on gene order. Additionally, we asked whether the distribution of transfer RNAs (tRNAs) evolves independently to that of mt protein-coding genes. We found that fungal mt genomes display remarkable variation between and within the major fungal phyla in terms of gene order, genome size, composition of intergenic regions, and presence of repeats, introns, and associated ORFs. Our results support previous evidence for the presence of mt recombination in all fungal phyla, a process conspicuously lacking in most Metazoa. Overall, the patterns of rearrangements may be explained by the combined influences of recombination (i.e., most likely nonhomologous and intrachromosomal), accumulated repeats, especially at intergenic regions, and to a lesser extent, mobile element dynamics.

Microsatellite loci to recognize species for the cheese starter and contaminating strains associated with cheese manufacturing.

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 P. chrysogenum 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 synonymized 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.

Phylo-MCOA: A Fast and Efficient Method to Detect Outlier Genes and Species in Phylogenomics Using Multiple Co-inertia Analysis

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.

Sex in Penicillium: Combined phylogenetic and experimental approaches

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.

Chaos of Rearrangements in the Mating-Type Chromosomes of the Anther-Smut Fungus Microbotryum lychnidis-dioicae

Sex chromosomes in plants and animals and fungal mating-type chromosomes often show exceptional genome features, with extensive suppression of homologous recombination and cytological differentiation between members of the diploid chromosome pair. Despite strong interest in the genetics of these chromosomes, their large regions of suppressed recombination often are enriched in transposable elements and therefore can be challenging to assemble. Here we show that the latest improvements of the PacBio sequencing yield assembly of the whole genome of the anther-smut fungus, Microbotryum lychnidis-dioicae (the pathogenic fungus causing anther-smut disease of Silene latifolia), into finished chromosomes or chromosome arms, even for the repeat-rich mating-type chromosomes and centromeres. Suppressed recombination of the mating-type chromosomes is revealed to span nearly 90% of their lengths, with extreme levels of rearrangements, transposable element accumulation, and differentiation between the two mating types. We observed no correlation between allelic divergence and physical position in the nonrecombining regions of the mating-type chromosomes. This may result from gene conversion or from rearrangements of ancient evolutionary strata, i.e., successive steps of suppressed recombination. Centromeres were found to be composed mainly of copia-like transposable elements and to possess specific minisatellite repeats identical between the different chromosomes. We also identified subtelomeric motifs. In addition, extensive signs of degeneration were detected in the nonrecombining regions in the form of transposable element accumulation and of hundreds of gene losses on each mating-type chromosome. Furthermore, our study highlights the potential of the latest breakthrough PacBio chemistry to resolve complex genome architectures.