Natural proteome diversity links aneuploidy tolerance to protein turnover.

Accessing the natural genetic diversity of species unveils hidden genetic traits, clarifies gene functions and allows the generalizability of laboratory findings to be assessed. One notable discovery made in natural isolates of Saccharomyces cerevisiae is that aneuploidy-an imbalance in chromosome copy numbers-is frequent1,2 (around 20%), which seems to contradict the substantial fitness costs and transient nature of aneuploidy when it is engineered in the laboratory3-5. Here we generate a proteomic resource and merge it with genomic1 and transcriptomic6 data for 796 euploid and aneuploid natural isolates. We find that natural and lab-generated aneuploids differ specifically at the proteome. In lab-generated aneuploids, some proteins-especially subunits of protein complexes-show reduced expression, but the overall protein levels correspond to the aneuploid gene dosage. By contrast, in natural isolates, more than 70% of proteins encoded on aneuploid chromosomes are dosage compensated, and average protein levels are shifted towards the euploid state chromosome-wide. At the molecular level, we detect an induction of structural components of the proteasome, increased levels of ubiquitination, and reveal an interdependency of protein turnover rates and attenuation. Our study thus highlights the role of protein turnover in mediating aneuploidy tolerance, and shows the utility of exploiting the natural diversity of species to attain generalizable molecular insights into complex biological processes.

Ancient and recent origins of shared polymorphisms in yeast.

Shared genetic polymorphisms between populations and species can be ascribed to ancestral variation or to more recent gene flow. Here, we mapped shared polymorphisms in Saccharomyces cerevisiae and its sister species Saccharomyces paradoxus, which diverged 4-6 million years ago. We used a dense map of single-nucleotide diagnostic markers (mean distance 15.6 base pairs) in 1,673 sequenced S. cerevisiae isolates to catalogue 3,852 sequence blocks (≥5 consecutive markers) introgressed from S. paradoxus, with most being recent and clade-specific. The highly diverged wild Chinese S. cerevisiae lineages were depleted of introgressed blocks but retained an excess of individual ancestral polymorphisms derived from incomplete lineage sorting, perhaps due to less dramatic population bottlenecks. In the non-Chinese S. cerevisiae lineages, we inferred major hybridization events and detected cases of overlapping introgressed blocks across distinct clades due to either shared histories or convergent evolution. We experimentally engineered, in otherwise isogenic backgrounds, the introgressed PAD1-FDC1 gene pair that independently arose in two S. cerevisiae clades and revealed that it increases resistance against diverse antifungal drugs. Overall, our study retraces the histories of divergence and secondary contacts across S. cerevisiae and S. paradoxus populations and unveils a functional outcome.

Single nucleotide polymorphisms associated with wine fermentation and adaptation to nitrogen limitation in wild and domesticated yeast strains.

For more than 20 years, Saccharomyces cerevisiae has served as a model organism for genetic studies and molecular biology, as well as a platform for biotechnology (e.g., wine production). One of the important ecological niches of this yeast that has been extensively studied is wine fermentation, a complex microbiological process in which S. cerevisiae faces various stresses such as limited availability of nitrogen. Nitrogen deficiencies in grape juice impair fermentation rate and yeast biomass production, leading to sluggish or stuck fermentations, resulting in considerable economic losses for the wine industry. In the present work, we took advantage of the "1002 Yeast Genomes Project" population, the most complete catalogue of the genetic variation in the species and a powerful resource for genotype-phenotype correlations, to study the adaptation to nitrogen limitation in wild and domesticated yeast strains in the context of wine fermentation. We found that wild and domesticated yeast strains have different adaptations to nitrogen limitation, corroborating their different evolutionary trajectories. Using a combination of state-of-the-art bioinformatic (GWAS) and molecular biology (CRISPR-Cas9) methodologies, we validated that PNP1, RRT5 and PDR12 are implicated in wine fermentation, where RRT5 and PDR12 are also involved in yeast adaptation to nitrogen limitation. In addition, we validated SNPs in these genes leading to differences in fermentative capacities and adaptation to nitrogen limitation. Altogether, the mapped genetic variants have potential applications for the genetic improvement of industrial yeast strains.

Telomere-to-telomere assemblies of 142 strains characterize the genome structural landscape in Saccharomyces cerevisiae.

Pangenomes provide access to an accurate representation of the genetic diversity of species, both in terms of sequence polymorphisms and structural variants (SVs). Here we generated the Saccharomyces cerevisiae Reference Assembly Panel (ScRAP) comprising reference-quality genomes for 142 strains representing the species' phylogenetic and ecological diversity. The ScRAP includes phased haplotype assemblies for several heterozygous diploid and polyploid isolates. We identified circa (ca.) 4,800 nonredundant SVs that provide a broad view of the genomic diversity, including the dynamics of telomere length and transposable elements. We uncovered frequent cases of complex aneuploidies where large chromosomes underwent large deletions and translocations. We found that SVs can impact gene expression near the breakpoints and substantially contribute to gene repertoire evolution. We also discovered that horizontally acquired regions insert at chromosome ends and can generate new telomeres. Overall, the ScRAP demonstrates the benefit of a pangenome in understanding genome evolution at population scale.

Disadvantaged Economic Conditions and Stricter Border Rules Shape Afghan Refugees' Ethnobotany: Insights from Kohat District, NW Pakistan.

The study of migrants' ethnobotany can help to address the diverse socio-ecological factors affecting temporal and spatial changes in local ecological knowledge (LEK). Through semi-structured and in-depth conversations with ninety interviewees among local Pathans and Afghan refugees in Kohat District, NW Pakistan, one hundred and forty-five wild plant and mushroom folk taxa were recorded. The plants quoted by Afghan refugees living inside and outside the camps tend to converge, while the Afghan data showed significant differences with those collected by local Pakistani Pathans. Interviewees mentioned two main driving factors potentially eroding folk plant knowledge: (a) recent stricter border policies have made it more difficult for refugees to visit their home regions in Afghanistan and therefore to also procure plants in their native country; (b) their disadvantaged economic conditions have forced them to engage more and more in urban activities in the host country, leaving little time for farming and foraging practices. Stakeholders should foster the exposure that refugee communities have to their plant resources, try to increase their socio-economic status, and facilitate both their settling outside the camps and their transnational movement for enhancing their use of wild plants, ultimately leading to improvements in their food security and health status.

Telomeres are shorter in wild Saccharomyces cerevisiae isolates than in domesticated ones.

Telomeres are ribonucleoproteins that cap chromosome-ends and their DNA length is controlled by counteracting elongation and shortening processes. The budding yeast Saccharomyces cerevisiae has been a leading model to study telomere DNA length control and dynamics. Its telomeric DNA is maintained at a length that slightly varies between laboratory strains, but little is known about its variation at the species level. The recent publication of the genomes of over 1,000 S. cerevisiae strains enabled us to explore telomere DNA length variation at an unprecedented scale. Here, we developed a bioinformatic pipeline (YeaISTY) to estimate telomere DNA length from whole-genome sequences and applied it to the sequenced S. cerevisiae collection. Our results revealed broad natural telomere DNA length variation among the isolates. Notably, telomere DNA length is shorter in those derived from wild rather than domesticated environments. Moreover, telomere DNA length variation is associated with mitochondrial metabolism, and this association is driven by wild strains. Overall, these findings reveal broad variation in budding yeast's telomere DNA length regulation, which might be shaped by its different ecological life-styles.

miRNAs Copy Number Variations Repertoire as Hallmark Indicator of Cancer Species Predisposition.

Aging is one of the hallmarks of multiple human diseases, including cancer. We hypothesized that variations in the number of copies (CNVs) of specific genes may protect some long-living organisms theoretically more susceptible to tumorigenesis from the onset of cancer. Based on the statistical comparison of gene copy numbers within the genomes of both cancer-prone and -resistant species, we identified novel gene targets linked to tumor predisposition, such as CD52, SAT1 and SUMO. Moreover, considering their genome-wide copy number landscape, we discovered that microRNAs (miRNAs) are among the most significant gene families enriched for cancer progression and predisposition. Through bioinformatics analyses, we identified several alterations in miRNAs copy number patterns, involving miR-221, miR-222, miR-21, miR-372, miR-30b, miR-30d and miR-31, among others. Therefore, our analyses provide the first evidence that an altered miRNAs copy number signature can statistically discriminate species more susceptible to cancer from those that are tumor resistant, paving the way for further investigations.

Genome-Wide Association Study Reveals Host Factors Affecting Conjugation in Escherichia coli.

The emergence and dissemination of antibiotic resistance threaten the treatment of common bacterial infections. Resistance genes are often encoded on conjugative elements, which can be horizontally transferred to diverse bacteria. In order to delay conjugative transfer of resistance genes, more information is needed on the genetic determinants promoting conjugation. Here, we focus on which bacterial host factors in the donor assist transfer of conjugative plasmids. We introduced the broad-host-range plasmid pKJK10 into a diverse collection of 113 Escherichia coli strains and measured by flow cytometry how effectively each strain transfers its plasmid to a fixed E. coli recipient. Differences in conjugation efficiency of up to 2.7 and 3.8 orders of magnitude were observed after mating for 24 h and 48 h, respectively. These differences were linked to the underlying donor strain genetic variants in genome-wide association studies, thereby identifying candidate genes involved in conjugation. We confirmed the role of fliF, fliK, kefB and ucpA in the donor ability of conjugative elements by validating defects in the conjugation efficiency of the corresponding lab strain single-gene deletion mutants. Based on the known cellular functions of these genes, we suggest that the motility and the energy supply, the intracellular pH or salinity of the donor affect the efficiency of plasmid transfer. Overall, this work advances the search for targets for the development of conjugation inhibitors, which can be administered alongside antibiotics to more effectively treat bacterial infections.

Domestication reprogrammed the budding yeast life cycle.

Domestication of plants and animals is the foundation for feeding the world human population but can profoundly alter the biology of the domesticated species. Here we investigated the effect of domestication on one of our prime model organisms, the yeast Saccharomyces cerevisiae, at a species-wide level. We tracked the capacity for sexual and asexual reproduction and the chronological life span across a global collection of 1,011 genome-sequenced yeast isolates and found a remarkable dichotomy between domesticated and wild strains. Domestication had systematically enhanced fermentative and reduced respiratory asexual growth, altered the tolerance to many stresses and abolished or impaired the sexual life cycle. The chronological life span remained largely unaffected by domestication and was instead dictated by clade-specific evolution. We traced the genetic origins of the yeast domestication syndrome using genome-wide association analysis and genetic engineering and disclosed causative effects of aneuploidy, gene presence/absence variations, copy number variations and single-nucleotide polymorphisms. Overall, we propose domestication to be the most dramatic event in budding yeast evolution, raising questions about how much domestication has distorted our understanding of the natural biology of this key model species.

Aborting meiosis allows recombination in sterile diploid yeast hybrids.

Hybrids between diverged lineages contain novel genetic combinations but an impaired meiosis often makes them evolutionary dead ends. Here, we explore to what extent an aborted meiosis followed by a return-to-growth (RTG) promotes recombination across a panel of 20 Saccharomyces cerevisiae and S. paradoxus diploid hybrids with different genomic structures and levels of sterility. Genome analyses of 275 clones reveal that RTG promotes recombination and generates extensive regions of loss-of-heterozygosity in sterile hybrids with either a defective meiosis or a heavily rearranged karyotype, whereas RTG recombination is reduced by high sequence divergence between parental subgenomes. The RTG recombination preferentially arises in regions with low local heterozygosity and near meiotic recombination hotspots. The loss-of-heterozygosity has a profound impact on sexual and asexual fitness, and enables genetic mapping of phenotypic differences in sterile lineages where linkage analysis would fail. We propose that RTG gives sterile yeast hybrids access to a natural route for genome recombination and adaptation.

On the Trail of an Ancient Middle Eastern Ethnobotany: Traditional Wild Food Plants Gathered by Ormuri Speakers in Kaniguram, NW Pakistan.

An ethnobotanical field study focusing on traditional wild food botanical taxa was carried out in Kaniguram, South Waziristan, Pakistan, among Ormur (or Burki or Baraki) peoples, which represent a diasporic minority group, as well as among the surrounding Pashtuns. Through sixty semi-structured interviews, fifty-two wild food plants (taxa) were recorded, and they were primarily used raw as snacks and cooked as vegetables. Comparative analysis found a remarkable overlap of the quoted plant uses between the two studied groups, which may reflect complex socio-cultural adaptations Ormur speakers faced. Ormur people retain a rich knowledge of anthropogenic weeds and the phytonyms reveal important commonalities with Persian and Kurdish phytonyms, which may indicate their possible horticultural-driven human ecological origin from the Middle East. Some novel or rare food uses of Cirsiumarvense, Nannorrhops ritchiana, Periploca aphylla, Perovskia atriplicifolia, Viscum album,Oxalis corniculata and Withania coagulans were documented. Since the Ormuri language represents a moribund language, still spoken by only a few thousand speakers in NW Pakistan and Afghanistan, it is recommended that the traditional bio-cultural and gastronomical heritage of this minority group be appropriately protected and bolstered in future rural development programs.

Deconvolution of intergenic polymorphisms determining high expression of Factor H binding protein in meningococcus and their association with invasive disease.

Neisseria meningitidis is a strictly human pathogen and is the major cause of septicemia and meningitis worldwide. Factor H binding protein (fHbp) is a meningococcal surface-exposed lipoprotein that binds the human Complement factor H allowing the bacterium to evade the host innate immune response. FHbp is also a key antigen in two vaccines against N. meningitidis serogroup B. Although the fHbp gene is present in most circulating meningococcal strains, level of fHbp expression varies among isolates and has been correlated to differences in promoter sequences upstream of the gene. Here we elucidated the sequence determinants that control fHbp expression in globally circulating strains. We analyzed the upstream fHbp intergenic region (fIR) of more than 5800 strains representative of the UK circulating isolates and we identified eleven fIR sequence alleles which represent 88% of meningococcal strains. By engineering isogenic recombinant strains where fHbp expression was under the control of each of the eleven fIR alleles, we confirmed that the fIR sequence determines a specific and distinct level of expression. Moreover, we identified the molecular basis for variation in expression through polymorphisms within key regulatory regions that are known to affect fHbp expression. We experimentally established three expression groups, high-medium-low, that correlated directly with the susceptibility to killing mediated by anti-fHbp antibodies and the ability of the meningococcal strain to survive within human serum. By using this sequence classification and information about the variant, we predicted fHbp expression in the panel of UK strains and we observed that strains with higher expressing fIR alleles are more likely associated with invasive disease. Overall, our findings can contribute to understand and predict vaccine coverage mediated by fHbp as well as to shed light on the role of this virulence factor in determining an invasive phenotype.

Lager Yeast Design Through Meiotic Segregation of a Saccharomyces cerevisiae × Saccharomyces eubayanus Hybrid.

Yeasts in the lager brewing group are closely related and consequently do not exhibit significant genetic variability. Here, an artificial Saccharomyces cerevisiae × Saccharomyces eubayanus tetraploid interspecies hybrid was created by rare mating, and its ability to sporulate and produce viable gametes was exploited to generate phenotypic diversity. Four spore clones obtained from a single ascus were isolated, and their brewing-relevant phenotypes were assessed. These F1 spore clones were found to differ with respect to fermentation performance under lager brewing conditions (15°C, 15 °Plato), production of volatile aroma compounds, flocculation potential and temperature tolerance. One spore clone, selected for its rapid fermentation and acetate ester production was sporulated to produce an F2 generation, again comprised of four spore clones from a single ascus. Again, phenotypic diversity was introduced. In two of these F2 clones, the fermentation performance was maintained and acetate ester production was improved relative to the F1 parent and the original hybrid strain. Strains also performed well in comparison to a commercial lager yeast strain. Spore clones varied in ploidy and chromosome copy numbers, and faster wort fermentation was observed in strains with a higher ploidy. An F2 spore clone was also subjected to 10 consecutive wort fermentations, and single cells were isolated from the resulting yeast slurry. These isolates also exhibited variable fermentation performance and chromosome copy numbers, highlighting the instability of polyploid interspecific hybrids. These results demonstrate the value of this natural approach to increase the phenotypic diversity of lager brewing yeast strains.

Meiotic Cells Counteract Programmed Retrotransposon Activation via RNA-Binding Translational Repressor Assemblies.

Retrotransposon proliferation poses a threat to germline integrity. While retrotransposons must be activated in developing germ cells in order to survive and propagate, how they are selectively activated in the context of meiosis is unclear. We demonstrate that the transcriptional activation of Ty3/Gypsy retrotransposons and host defense are controlled by master meiotic regulators. We show that budding yeast Ty3/Gypsy co-opts binding sites of the essential meiotic transcription factor Ndt80 upstream of the integration site, thereby tightly linking its transcriptional activation to meiotic progression. We also elucidate how yeast cells thwart Ty3/Gypsy proliferation by blocking translation of the retrotransposon mRNA using amyloid-like assemblies of the RNA-binding protein Rim4. In mammals, several inactive Ty3/Gypsy elements are undergoing domestication. We show that mammals utilize equivalent master meiotic regulators (Stra8, Mybl1, Dazl) to regulate Ty3/Gypsy-derived genes in developing gametes. Our findings inform how genes that are evolving from retrotransposons can build upon existing regulatory networks during domestication.

A yeast living ancestor reveals the origin of genomic introgressions.

Genome introgressions drive evolution across the animal1, plant2 and fungal3 kingdoms. Introgressions initiate from archaic admixtures followed by repeated backcrossing to one parental species. However, how introgressions arise in reproductively isolated species, such as yeast4, has remained unclear. Here we identify a clonal descendant of the ancestral yeast hybrid that founded the extant Saccharomyces cerevisiae Alpechin lineage5, which carries abundant Saccharomyces paradoxus introgressions. We show that this clonal descendant, hereafter defined as a 'living ancestor', retained the ancestral genome structure of the first-generation hybrid with contiguous S. cerevisiae and S. paradoxus subgenomes. The ancestral first-generation hybrid underwent catastrophic genomic instability through more than a hundred mitotic recombination events, mainly manifesting as homozygous genome blocks generated by loss of heterozygosity. These homozygous sequence blocks rescue hybrid fertility by restoring meiotic recombination and are the direct origins of the introgressions present in the Alpechin lineage. We suggest a plausible route for introgression evolution through the reconstruction of extinct stages and propose that genome instability allows hybrids to overcome reproductive isolation and enables introgressions to emerge.

CRISpy-Pop: A Web Tool for Designing CRISPR/Cas9-Driven Genetic Modifications in Diverse Populations.

CRISPR/Cas9 is a powerful tool for editing genomes, but design decisions are generally made with respect to a single reference genome. With population genomic data becoming available for an increasing number of model organisms, researchers are interested in manipulating multiple strains and lines. CRISpy-pop is a web application that generates and filters guide RNA sequences for CRISPR/Cas9 genome editing for diverse yeast and bacterial strains. The current implementation designs and predicts the activity of guide RNAs against more than 1000 Saccharomyces cerevisiae genomes, including 167 strains frequently used in bioenergy research. Zymomonas mobilis, an increasingly popular bacterial bioenergy research model, is also supported. CRISpy-pop is available as a web application (https://CRISpy-pop.glbrc.org/) with an intuitive graphical user interface. CRISpy-pop also cross-references the human genome to allow users to avoid the selection of guide RNAs with potential biosafety concerns. Additionally, CRISpy-pop predicts the strain coverage of each guide RNA within the supported strain sets, which aids in functional population genetic studies. Finally, we validate how CRISpy-pop can accurately predict the activity of guide RNAs across strains using population genomic data.

An Out-of-Patagonia migration explains the worldwide diversity and distribution of Saccharomyces eubayanus lineages.

Population-level sampling and whole-genome sequences of different individuals allow one to identify signatures of hybridization, gene flow and potential molecular mechanisms of environmental responses. Here, we report the isolation of 160 Saccharomyces eubayanus strains, the cryotolerant ancestor of lager yeast, from ten sampling sites in Patagonia along 2,000 km of Nothofagus forests. Frequency of S. eubayanus isolates was higher towards southern and colder regions, demonstrating the cryotolerant nature of the species. We sequenced the genome of 82 strains and, together with 23 available genomes, performed a comprehensive phylogenetic analysis. Our results revealed the presence of five different lineages together with dozens of admixed strains. Various analytical methods reveal evidence of gene flow and historical admixture between lineages from Patagonia and Holarctic regions, suggesting the co-occurrence of these ancestral populations. Analysis of the genetic contribution to the admixed genomes revealed a Patagonian genetic origin of the admixed strains, even for those located in the North Hemisphere. Overall, the Patagonian lineages, particularly the southern populations, showed a greater global genetic diversity compared to Holarctic and Chinese lineages, in agreement with a higher abundance in Patagonia. Thus, our results are consistent with a likely colonization of the species from peripheral glacial refugia from South Patagonia. Furthermore, fermentative capacity and maltose consumption resulted negatively correlated with latitude, indicating better fermentative performance in northern populations. Our genome analysis, together with previous reports in the sister species S. uvarum suggests that a S. eubayanus ancestor was adapted to the harsh environmental conditions of Patagonia, a region that provides the ecological conditions for the diversification of these ancestral lineages.

Discordant evolution of mitochondrial and nuclear yeast genomes at population level.

BACKGROUND: Mitochondria are essential organelles partially regulated by their own genomes. The mitochondrial genome maintenance and inheritance differ from the nuclear genome, potentially uncoupling their evolutionary trajectories. Here, we analysed mitochondrial sequences obtained from the 1011 Saccharomyces cerevisiae strain collection and identified pronounced differences with their nuclear genome counterparts. RESULTS: In contrast with pre-whole genome duplication fungal species, S. cerevisiae mitochondrial genomes show higher genetic diversity compared to the nuclear genomes. Strikingly, mitochondrial genomes appear to be highly admixed, resulting in a complex interconnected phylogeny with a weak grouping of isolates, whereas interspecies introgressions are very rare. Complete genome assemblies revealed that structural rearrangements are nearly absent with rare inversions detected. We tracked intron variation in COX1 and COB to infer gain and loss events throughout the species evolutionary history. Mitochondrial genome copy number is connected with the nuclear genome and linearly scale up with ploidy. We observed rare cases of naturally occurring mitochondrial DNA loss, petite, with a subset of them that do not suffer the expected growth defect in fermentable rich media. CONCLUSIONS: Overall, our results illustrate how differences in the biology of two genomes coexisting in the same cells can lead to discordant evolutionary histories.

Transcriptional Activity and Protein Levels of Horizontally Acquired Genes in Yeast Reveal Hallmarks of Adaptation to Fermentative Environments.

In the past decade, the sequencing of large cohorts of Saccharomyces cerevisiae strains has revealed a landscape of genomic regions acquired by Horizontal Gene Transfer (HGT). The genes acquired by HGT play important roles in yeast adaptation to the fermentation process, improving nitrogen and carbon source utilization. However, the functional characterization of these genes at the molecular level has been poorly attended. In this work, we carried out a systematic analysis of the promoter activity and protein level of 30 genes contained in three horizontally acquired regions commonly known as regions A, B, and C. In three strains (one for each region), we used the luciferase reporter gene and the mCherry fluorescent protein to quantify the transcriptional and translational activity of these genes, respectively. We assayed the strains generated in four different culture conditions; all showed low levels of transcriptional and translational activity across these environments. However, we observed an increase in protein levels under low nitrogen culture conditions, suggesting a possible role of the horizontally acquired genes in the adaptation to nitrogen-limited environments. Furthermore, since the strains carrying the luciferase reporter gene are null mutants for the horizontally acquired genes, we assayed growth parameters (latency time, growth rate, and efficiency) and the fermentation kinetics in this set of deletion strains. The results showed that single deletion of 20 horizontally acquired genes modified the growth parameters, whereas the deletion of five of them altered the maximal CO2 production rate (Vmax). Interestingly, we observed a correlation between growth parameters and Vmax for an ORF within region A, encoding an ortholog to a thiamine (vitamin B1) transporter whose deletion decreased the growth rate, growth efficiency, and CO2 production. Altogether, our results provided molecular and phenotypic evidence highlighting the importance of horizontally acquired genes in yeast adaptation to fermentative environments.

Comparison of Phylogenetic Tree Topologies for Nitrogen Associated Genes Partially Reconstruct the Evolutionary History of Saccharomyces cerevisiae.

Massive sequencing projects executed in Saccharomyces cerevisiae have revealed in detail its population structure. The recent "1002 yeast genomes project" has become the most complete catalogue of yeast genetic diversity and a powerful resource to analyse the evolutionary history of genes affecting specific phenotypes. In this work, we selected 22 nitrogen associated genes and analysed the sequence information from the 1011 strains of the "1002 yeast genomes project". We constructed a total evidence (TE) phylogenetic tree using concatenated information, which showed a 27% topology similarity with the reference (REF) tree of the "1002 yeast genomes project". We also generated individual phylogenetic trees for each gene and compared their topologies, identifying genes with similar topologies (suggesting a shared evolutionary history). Furthermore, we pruned the constructed phylogenetic trees to compare the REF tree topology versus the TE tree and the individual genes trees, considering each phylogenetic cluster/subcluster within the population, observing genes with cluster/subcluster topologies of high similarity to the REF tree. Finally, we used the pruned versions of the phylogenetic trees to compare four strains considered as representatives of S. cerevisiae clean lineages, observing for 15 genes that its cluster topologies match 100% the REF tree, supporting that these strains represent main lineages of yeast population. Altogether, our results showed the potential of tree topologies comparison for exploring the evolutionary history of a specific group of genes.