A major bottleneck in our understanding of the molecular underpinnings of life is the assignment of function to proteins. While molecular experiments provide the most reliable annotation of proteins, their relatively low throughput and restricted purview have led to an increasing role for computational function prediction. However, assessing methods for protein function prediction and tracking progress in the field remain challenging. We conducted the second critical assessment of functional annotation (CAFA), a timed challenge to assess computational methods that automatically assign protein function. We evaluated 126 methods from 56 research groups for their ability to predict biological functions using Gene Ontology and gene-disease associations using Human Phenotype Ontology on a set of 3681 proteins from 18 species. CAFA2 featured expanded analysis compared with CAFA1, with regards to data set size, variety, and assessment metrics. To review progress in the field, the analysis compared the best methods from CAFA1 to those of CAFA2. The top-performing methods in CAFA2 outperformed those from CAFA1. This increased accuracy can be attributed to a combination of the growing number of experimental annotations and improved methods for function prediction. The assessment also revealed that the definition of top-performing algorithms is ontology specific, that different performance metrics can be used to probe the nature of accurate predictions, and the relative diversity of predictions in the biological process and human phenotype ontologies. While there was methodological improvement between CAFA1 and CAFA2, the interpretation of results and usefulness of individual methods remain context-dependent.
Whole-genome duplications have shaped the genomes of several vertebrate, plant, and fungal lineages. Earlier studies have focused on establishing when these events occurred and on elucidating their functional and evolutionary consequences, but we still lack sufficient understanding of how genome duplications first originated. We used phylogenomics to study the ancient genome duplication occurred in the yeast Saccharomyces cerevisiae lineage and found compelling evidence for the existence of a contemporaneous interspecies hybridization. We propose that the genome doubling was a direct consequence of this hybridization and that it served to provide stability to the recently formed allopolyploid. This scenario provides a mechanism for the origin of this ancient duplication and the lineage that originated from it and brings a new perspective to the interpretation of the origin and consequences of whole-genome duplications. Genome duplication is a major evolutionary process that has shaped the genomes of several eukaryotic lineages including vertebrates, plants, and fungi. The sequencing of the baker's yeast Saccharomyces cerevisiae in the 1990s revealed the presence of conserved blocks of duplicated genes, indicating an ancestral duplication of the entire genome. Subsequent work has clarified when this event occurred and what genomic rearrangements followed, but the underlying mechanistic origin of such a large-scale event remains poorly understood. Here we used a large-scale phylogenetic approach to examine the individual evolutionary histories of all yeast genes and assessed the time at which each duplication occurred. This survey revealed evidence for an ancient hybridization event between two ancestral species in the lineage in which the whole-genome duplication had occurred. We further characterize this hybridization event and the properties of the putative parental species. We propose that the whole-genome duplication was a direct consequence of this hybridization, providing a means by which the initially sterile hybrid could regain fertility. This scenario provides a mechanistic understanding of the origin of the ancient yeast whole-genome duplication and brings a radically different perspective on the interpretation of the origin and evolutionary consequences of whole-genome duplications in eukaryotic lineages.
The budding yeast Saccharomyces cerevisiae can be found in the wild and is also frequently associated with human activities. Despite recent insights into the phylogeny of this species, much is still unknown about how evolutionary processes related to anthropogenic niches have shaped the genomes and phenotypes of S. cerevisiae. To address this question, we performed population-level sequencing of 82 S. cerevisiae strains from wine, flor, rum, dairy products, bakeries, and the natural environment (oak trees). These genomic data enabled us to delineate specific genetic groups corresponding to the different ecological niches and revealed high genome content variation across the groups. Most of these strains, compared with the reference genome, possessed additional genetic elements acquired by introgression or horizontal transfer, several of which were population-specific. In addition, several genomic regions in each population showed evidence of nonneutral evolution, as shown by high differentiation, or of selective sweeps including genes with key functions in these environments (e.g., amino acid transport for wine yeast). Linking genetics to lifestyle differences and metabolite traits has enabled us to elucidate the genetic basis of several niche-specific population traits, such as growth on galactose for cheese strains. These data indicate that yeast has been subjected to various divergent selective pressures depending on its niche, requiring the development of customized genomes for better survival in these environments. These striking genome dynamics associated with local adaptation and domestication reveal the remarkable plasticity of the S. cerevisiae genome, revealing this species to be an amazing complex of specialized populations.
Phylogenetic trees representing the evolutionary relationships of homologous genes are the entry point for many evolutionary analyses. For instance, the use of a phylogenetic tree can aid in the inference of orthology and paralogy relationships, and in the detection of relevant evolutionary events such as gene family expansions and contractions, horizontal gene transfer, recombination or incomplete lineage sorting. Similarly, given the plurality of evolutionary histories among genes encoded in a given genome, there is a need for the combined analysis of genome-wide collections of phylogenetic trees (phylomes). Here, we introduce a new release of PhylomeDB (http://phylomedb.org), a public repository of phylomes. Currently, PhylomeDB hosts 120 public phylomes, comprising >1.5 million maximum likelihood trees and multiple sequence alignments. In the current release, phylogenetic trees are annotated with taxonomic, protein-domain arrangement, functional and evolutionary information. PhylomeDB is also a major source for phylogeny-based predictions of orthology and paralogy, covering >10 million proteins across 1059 sequenced species. Here we describe newly implemented PhylomeDB features, and discuss a benchmark of the orthology predictions provided by the database, the impact of proteome updates and the use of the phylome approach in the analysis of newly sequenced genomes and transcriptomes.
Background: The Mediterranean olive tree (Olea europaea subsp. europaea) was one of the first trees to be domesticated and is currently of major agricultural importance in the Mediterranean region as the source of olive oil. The molecular bases underlying the phenotypic differences among domesticated cultivars, or between domesticated olive trees and their wild relatives, remain poorly understood. Both wild and cultivated olive trees have 46 chromosomes (2n). Findings: A total of 543 Gb of raw DNA sequence from whole genome shotgun sequencing, and a fosmid library containing 155,000 clones from a 1,000+ year-old olive tree (cv. Farga) were generated by Illumina sequencing using different combinations of mate-pair and pair-end libraries. Assembly gave a final genome with a scaffold N50 of 443 kb, and a total length of 1.31 Gb, which represents 95 % of the estimated genome length (1.38 Gb). In addition, the associated fungus Aureobasidium pullulans was partially sequenced. Genome annotation, assisted by RNA sequencing from leaf, root, and fruit tissues at various stages, resulted in 56,349 unique protein coding genes, suggesting recent genomic expansion. Genome completeness, as estimated using the CEGMA pipeline, reached 98.79 %. Conclusions: The assembled draft genome of O. europaea will provide a valuable resource for the study of the evolution and domestication processes of this important tree, and allow determination of the genetic bases of key phenotypic traits. Moreover, it will enhance breeding programs and the formation of new varieties.
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.