The metabolic domestication syndrome of budding yeast.

Cellular metabolism evolves through changes in the structure and quantitative states of metabolic networks. Here, we explore the evolutionary dynamics of metabolic states by focusing on the collection of metabolite levels, the metabolome, which captures key aspects of cellular physiology. Using a phylogenetic framework, we profiled metabolites in 27 populations of nine budding yeast species, providing a graduated view of metabolic variation across multiple evolutionary time scales. Metabolite levels evolve more rapidly and independently of changes in the metabolic network's structure, providing complementary information to enzyme repertoire. Although metabolome variation accumulates mainly gradually over time, it is profoundly affected by domestication. We found pervasive signatures of convergent evolution in the metabolomes of independently domesticated clades of Saccharomyces cerevisiae. Such recurring metabolite differences between wild and domesticated populations affect a substantial part of the metabolome, including rewiring of the TCA cycle and several amino acids that influence aroma production, likely reflecting adaptation to human niches. Overall, our work reveals previously unrecognized diversity in central metabolism and the pervasive influence of human-driven selection on metabolite levels in yeasts.

Plasticity and Stereotypic Rewiring of the Transcriptome Upon Bacterial Evolution of Antibiotic Resistance.

Bacterial evolution of antibiotic resistance frequently has deleterious side effects on microbial growth, virulence, and susceptibility to other antimicrobial agents. However, it is unclear how these trade-offs could be utilized for manipulating antibiotic resistance in the clinic, not least because the underlying molecular mechanisms are poorly understood. Using laboratory evolution, we demonstrate that clinically relevant resistance mutations in Escherichia coli constitutively rewire a large fraction of the transcriptome in a repeatable and stereotypic manner. Strikingly, lineages adapted to functionally distinct antibiotics and having no resistance mutations in common show a wide range of parallel gene expression changes that alter oxidative stress response, iron homeostasis, and the composition of the bacterial outer membrane and cell surface. These common physiological alterations are associated with changes in cell morphology and enhanced sensitivity to antimicrobial peptides. Finally, the constitutive transcriptomic changes induced by resistance mutations are largely distinct from those induced by antibiotic stresses in the wild type. This indicates a limited role for genetic assimilation of the induced antibiotic stress response during resistance evolution. Our work suggests that diverse resistance mutations converge on similar global transcriptomic states that shape genetic susceptibility to antimicrobial compounds.

The MCL1 inhibitor S63845 is tolerable and effective in diverse cancer models.

Avoidance of apoptosis is critical for the development and sustained growth of tumours. The pro-survival protein myeloid cell leukemia 1 (MCL1) is overexpressed in many cancers, but the development of small molecules targeting this protein that are amenable for clinical testing has been challenging. Here we describe S63845, a small molecule that specifically binds with high affinity to the BH3-binding groove of MCL1. Our mechanistic studies demonstrate that S63845 potently kills MCL1-dependent cancer cells, including multiple myeloma, leukaemia and lymphoma cells, by activating the BAX/BAK-dependent mitochondrial apoptotic pathway. In vivo, S63845 shows potent anti-tumour activity with an acceptable safety margin as a single agent in several cancers. Moreover, MCL1 inhibition, either alone or in combination with other anti-cancer drugs, proved effective against several solid cancer-derived cell lines. These results point towards MCL1 as a target for the treatment of a wide range of tumours.

Novel phylogenetic methods are needed for understanding gene function in the era of mega-scale genome sequencing.

Ongoing large-scale genome sequencing projects are forecasting a data deluge that will almost certainly overwhelm current analytical capabilities of evolutionary genomics. In contrast to population genomics, there are no standardized methods in evolutionary genomics for extracting evolutionary and functional (e.g. gene-trait association) signal from genomic data. Here, we examine how current practices of multi-species comparative genomics perform in this aspect and point out that many genomic datasets are under-utilized due to the lack of powerful methodologies. As a result, many current analyses emphasize gene families for which some functional data is already available, resulting in a growing gap between functionally well-characterized genes/organisms and the universe of unknowns. This leaves unknown genes on the 'dark side' of genomes, a problem that will not be mitigated by sequencing more and more genomes, unless we develop tools to infer functional hypotheses for unknown genes in a systematic manner. We provide an inventory of recently developed methods capable of predicting gene-gene and gene-trait associations based on comparative data, then argue that realizing the full potential of whole genome datasets requires the integration of phylogenetic comparative methods into genomics, a rich but underutilized toolbox for looking into the past.

Transcriptomic atlas of mushroom development reveals conserved genes behind complex multicellularity in fungi.

The evolution of complex multicellularity has been one of the major transitions in the history of life. In contrast to simple multicellular aggregates of cells, it has evolved only in a handful of lineages, including animals, embryophytes, red and brown algae, and fungi. Despite being a key step toward the evolution of complex organisms, the evolutionary origins and the genetic underpinnings of complex multicellularity are incompletely known. The development of fungal fruiting bodies from a hyphal thallus represents a transition from simple to complex multicellularity that is inducible under laboratory conditions. We constructed a reference atlas of mushroom formation based on developmental transcriptome data of six species and comparisons of >200 whole genomes, to elucidate the core genetic program of complex multicellularity and fruiting body development in mushroom-forming fungi (Agaricomycetes). Nearly 300 conserved gene families and >70 functional groups contained developmentally regulated genes from five to six species, covering functions related to fungal cell wall remodeling, targeted protein degradation, signal transduction, adhesion, and small secreted proteins (including effector-like orphan genes). Several of these families, including F-box proteins, expansin-like proteins, protein kinases, and transcription factors, showed expansions in Agaricomycetes, many of which convergently expanded in multicellular plants and/or animals too, reflecting convergent solutions to genetic hurdles imposed by complex multicellularity among independently evolved lineages. This study provides an entry point to studying mushroom development and complex multicellularity in one of the largest clades of complex eukaryotic organisms.

Unmatched Level of Molecular Convergence among Deeply Divergent Complex Multicellular Fungi.

Convergent evolution is pervasive in nature, but it is poorly understood how various constraints and natural selection limit the diversity of evolvable phenotypes. Here, we analyze the transcriptome across fruiting body development to understand the independent evolution of complex multicellularity in the two largest clades of fungi-the Agarico- and Pezizomycotina. Despite >650 My of divergence between these clades, we find that very similar sets of genes have convergently been co-opted for complex multicellularity, followed by expansions of their gene families by duplications. Over 82% of shared multicellularity-related gene families were expanding in both clades, indicating a high prevalence of convergence also at the gene family level. This convergence is coupled with a rich inferred repertoire of multicellularity-related genes in the most recent common ancestor of the Agarico- and Pezizomycotina, consistent with the hypothesis that the coding capacity of ancestral fungal genomes might have promoted the repeated evolution of complex multicellularity. We interpret this repertoire as an indication of evolutionary predisposition of fungal ancestors for evolving complex multicellular fruiting bodies. Our work suggests that evolutionary convergence may happen not only when organisms are closely related or are under similar selection pressures, but also when ancestral genomic repertoires render certain evolutionary trajectories more likely than others, even across large phylogenetic distances.

Gene family expansions and transcriptome signatures uncover fungal adaptations to wood decay.

Because they comprise some of the most efficient wood-decayers, Polyporales fungi impact carbon cycling in forest environment. Despite continuous discoveries on the enzymatic machinery involved in wood decomposition, the vision on their evolutionary adaptation to wood decay and genome diversity remains incomplete. We combined the genome sequence information from 50 Polyporales species, including 26 newly sequenced genomes and sought for genomic and functional adaptations to wood decay through the analysis of genome composition and transcriptome responses to different carbon sources. The genomes of Polyporales from different phylogenetic clades showed poor conservation in macrosynteny, indicative of genome rearrangements. We observed different gene family expansion/contraction histories for plant cell wall degrading enzymes in core polyporoids and phlebioids and captured expansions for genes involved in signalling and regulation in the lineages of white rotters. Furthermore, we identified conserved cupredoxins, thaumatin-like proteins and lytic polysaccharide monooxygenases with a yet uncharacterized appended module as new candidate players in wood decomposition. Given the current need for enzymatic toolkits dedicated to the transformation of renewable carbon sources, the observed genomic diversity among Polyporales strengthens the relevance of mining Polyporales biodiversity to understand the molecular mechanisms of wood decay.

Directed evolution of multiple genomic loci allows the prediction of antibiotic resistance.

Antibiotic development is frequently plagued by the rapid emergence of drug resistance. However, assessing the risk of resistance development in the preclinical stage is difficult. Standard laboratory evolution approaches explore only a small fraction of the sequence space and fail to identify exceedingly rare resistance mutations and combinations thereof. Therefore, new rapid and exhaustive methods are needed to accurately assess the potential of resistance evolution and uncover the underlying mutational mechanisms. Here, we introduce directed evolution with random genomic mutations (DIvERGE), a method that allows an up to million-fold increase in mutation rate along the full lengths of multiple predefined loci in a range of bacterial species. In a single day, DIvERGE generated specific mutation combinations, yielding clinically significant resistance against trimethoprim and ciprofloxacin. Many of these mutations have remained previously undetected or provide resistance in a species-specific manner. These results indicate pathogen-specific resistance mechanisms and the necessity of future narrow-spectrum antibacterial treatments. In contrast to prior claims, we detected the rapid emergence of resistance against gepotidacin, a novel antibiotic currently in clinical trials. Based on these properties, DIvERGE could be applicable to identify less resistance-prone antibiotics at an early stage of drug development. Finally, we discuss potential future applications of DIvERGE in synthetic and evolutionary biology.

Morphotype of bacteroids in different legumes correlates with the number and type of symbiotic NCR peptides.

In legume nodules, rhizobia differentiate into nitrogen-fixing forms called bacteroids, which are enclosed by a plant membrane in an organelle-like structure called the symbiosome. In the Inverted Repeat-Lacking Clade (IRLC) of legumes, this differentiation is terminal due to irreversible loss of cell division ability and is associated with genome amplification and different morphologies of the bacteroids that can be swollen, elongated, spherical, and elongated-branched, depending on the host plant. In Medicago truncatula, this process is orchestrated by nodule-specific cysteine-rich peptides (NCRs) delivered into developing bacteroids. Here, we identified the predicted NCR proteins in 10 legumes representing different subclades of the IRLC with distinct bacteroid morphotypes. Analysis of their expression and predicted sequences establishes correlations between the composition of the NCR family and the morphotypes of bacteroids. Although NCRs have a single origin, their evolution has followed different routes in individual lineages, and enrichment and diversification of cationic peptides has resulted in the ability to impose major morphological changes on the endosymbionts. The wide range of effects provoked by NCRs such as cell enlargement, membrane alterations and permeabilization, and biofilm and vesicle formation is dependent on the amino acid composition and charge of the peptides. These effects are strongly influenced by the rhizobial surface polysaccharides that affect NCR-induced differentiation and survival of rhizobia in nodule cells.

ELONGATED HYPOCOTYL 5 mediates blue light signalling to the Arabidopsis circadian clock.

Circadian clocks are gene networks producing 24-h oscillations at the level of clock gene expression that are synchronized to environmental cycles via light signals. The ELONGATED HYPOCOTYL 5 (HY5) transcription factor is a signalling hub acting downstream of several photoreceptors and is a key mediator of photomorphogenesis. Here we describe a mechanism by which light quality could modulate the pace of the circadian clock through governing abundance of HY5. We show that hy5 mutants display remarkably shorter period rhythms in blue but not in red light or darkness, and blue light is more efficient than red to induce accumulation of HY5 at transcriptional and post-transcriptional levels. We demonstrate that the pattern and level of HY5 accumulation modulates its binding to specific promoter elements of the majority of clock genes, but only a few of these show altered transcription in the hy5 mutant. Mathematical modelling suggests that the direct effect of HY5 on the apparently non-responsive clock genes could be masked by feedback from the clock gene network. We conclude that the information on the ratio of blue and red components of the white light spectrum is decoded and relayed to the circadian oscillator, at least partially, by HY5.

Comparative genomics reveals unique wood-decay strategies and fruiting body development in the Schizophyllaceae.

Agaricomycetes are fruiting body-forming fungi that produce some of the most efficient enzyme systems to degrade wood. Despite decades-long interest in their biology, the evolution and functional diversity of both wood-decay and fruiting body formation are incompletely known. We performed comparative genomic and transcriptomic analyses of wood-decay and fruiting body development in Auriculariopsis ampla and Schizophyllum commune (Schizophyllaceae), species with secondarily simplified morphologies, an enigmatic wood-decay strategy and weak pathogenicity to woody plants. The plant cell wall-degrading enzyme repertoires of Schizophyllaceae are transitional between those of white rot species and less efficient wood-degraders such as brown rot or mycorrhizal fungi. Rich repertoires of suberinase and tannase genes were found in both species, with tannases restricted to Agaricomycetes that preferentially colonize bark-covered wood, suggesting potential complementation of their weaker wood-decaying abilities and adaptations to wood colonization through the bark. Fruiting body transcriptomes revealed a high rate of divergence in developmental gene expression, but also several genes with conserved expression patterns, including novel transcription factors and small-secreted proteins, some of the latter which might represent fruiting body effectors. Taken together, our analyses highlighted novel aspects of wood-decay and fruiting body development in an important family of mushroom-forming fungi.

A highly precise and portable genome engineering method allows comparison of mutational effects across bacterial species.

Currently available tools for multiplex bacterial genome engineering are optimized for a few laboratory model strains, demand extensive prior modification of the host strain, and lead to the accumulation of numerous off-target modifications. Building on prior development of multiplex automated genome engineering (MAGE), our work addresses these problems in a single framework. Using a dominant-negative mutant protein of the methyl-directed mismatch repair (MMR) system, we achieved a transient suppression of DNA repair in Escherichia coli, which is necessary for efficient oligonucleotide integration. By integrating all necessary components into a broad-host vector, we developed a new workflow we term pORTMAGE. It allows efficient modification of multiple loci, without any observable off-target mutagenesis and prior modification of the host genome. Because of the conserved nature of the bacterial MMR system, pORTMAGE simultaneously allows genome editing and mutant library generation in other biotechnologically and clinically relevant bacterial species. Finally, we applied pORTMAGE to study a set of antibiotic resistance-conferring mutations in Salmonella enterica and E. coli. Despite over 100 million y of divergence between the two species, mutational effects remained generally conserved. In sum, a single transformation of a pORTMAGE plasmid allows bacterial species of interest to become an efficient host for genome engineering. These advances pave the way toward biotechnological and therapeutic applications. Finally, pORTMAGE allows systematic comparison of mutational effects and epistasis across a wide range of bacterial species.

Increased insulin-like growth factor 1 production by polyploid adipose stem cells promotes growth of breast cancer cells.

BACKGROUND: Adipose-tissue stem cells (ASCs) are subject of intensive research since their successful use in regenerative therapy. The drawback of ASCs is that they may serve as stroma for cancer cells and assist tumor progression. It is disquieting that ASCs frequently undergo genetic and epigenetic changes during their in vitro propagation. In this study, we describe the polyploidization of murine ASCs and the accompanying phenotypical, gene expressional and functional changes under long term culturing. METHODS: ASCs were isolated from visceral fat of C57BL/6 J mice, and cultured in vitro for prolonged time. The phenotypical changes were followed by microscopy and flow cytometry. Gene expressional changes were determined by differential transcriptome analysis and changes in protein expression were shown by Western blotting. The tumor growth promoting effect of ASCs was examined by co-culturing them with 4 T1 murine breast cancer cells. RESULTS: After five passages, the proliferation of ASCs decreases and cells enter a senescence-like state, from which a proportion of cells escape by polyploidization. The resulting ASC line is susceptible to adipogenic, osteogenic and chondrogenic differentiation, and expresses the stem cell markers CD29 and Sca-1 on an upregulated level. Differential transcriptome analysis of ASCs with normal and polyploid karyotype shows altered expression of genes that are involved in regulation of cancer, cellular growth and proliferation. We verified the increased expression of Klf4 and loss of Nestin on protein level. We found that elevated production of insulin-like growth factor 1 by polyploid ASCs rendered them more potent in tumor growth promotion in vitro. CONCLUSIONS: Our model indicates how ASCs with altered genetic background may support tumor progression.

Comprehensive DNA methylation study identifies novel progression-related and prognostic markers for cutaneous melanoma.

BACKGROUND: Cutaneous melanoma is the deadliest skin cancer, with an increasing incidence and mortality rate. Currently, staging of patients with primary melanoma is performed using histological biomarkers such as tumor thickness and ulceration. As disruption of the epigenomic landscape is recognized as a widespread feature inherent in tumor development and progression, we aimed to identify novel biomarkers providing additional clinical information over current factors using unbiased genome-wide DNA methylation analyses. METHODS: We performed a comprehensive DNA methylation analysis during all progression stages of melanoma using Infinium HumanMethylation450 BeadChips on a discovery cohort of benign nevi (n = 14) and malignant melanoma from both primary (n = 33) and metastatic (n = 28) sites, integrating the DNA methylome with gene expression data. We validated the discovered biomarkers in three independent validation cohorts by pyrosequencing and immunohistochemistry. RESULTS: We identified and validated biomarkers for, and pathways involved in, melanoma development (e.g., HOXA9 DNA methylation) and tumor progression (e.g., TBC1D16 DNA methylation). In addition, we determined a prognostic signature with potential clinical applicability and validated PON3 DNA methylation and OVOL1 protein expression as biomarkers with prognostic information independent of tumor thickness and ulceration. CONCLUSIONS: Our data underscores the importance of epigenomic regulation in triggering metastatic dissemination through the inactivation of central cancer-related pathways. Inactivation of cell-adhesion and differentiation unleashes dissemination, and subsequent activation of inflammatory and immune system programs impairs anti-tumoral defense pathways. Moreover, we identify several markers of tumor development and progression previously unrelated to melanoma, and determined a prognostic signature with potential clinical utility.

Specific Gene- and MicroRNA-Expression Pattern Contributes to the Epithelial to Mesenchymal Transition in a Rat Model of Experimental Colitis.

The aim of this study was to determine the gene- and microRNA-expression profile contributing to epithelial to mesenchymal transition in a rat model of experimental colitis. For this, inflammation was induced by injecting 2,4,6-trinitrobenzene sulphonic acid to the colon of male Wistar rats. Samples were taken from both inflamed and uninflamed regions of the same colon, total RNA was isolated, and the mRNA and microRNA expressions were monitored. We have determined that the expression of genes responsible for inducing mesenchymal phenotype, such as Egr1, Fgf2, Fgf7, Jak2, Notch2, Hif1α, Zeb2, Mmp9, Lox, and Vim, was all significantly induced in the inflamed regions of the affected colons while the epithelial marker E-cadherin (Cdh1) was downregulated. In contrast, the expression of microRNAs miR-192, miR-143, miR-375, miR-30a, miR-107, and miR-200b responsible for the regulation of the above mentioned genes was significantly downregulated in inflamed colon. Importantly, we detected moderate induction in the expression of five out of six tested microRNAs in the uninflamed regions. In summary, we identified numerous interacting genes and microRNAs with mutually exclusive expression pattern in inflamed regions of colitis-induced rats. These findings suggest that-among others-an important step in the epithelial to mesenchymal transition in experimental colitis is the dysregulated microRNA expression.

Precision medicine in cancer: challenges and recommendations from an EU-funded cervical cancer biobanking study.

BACKGROUND: Cervical cancer (CC) remains a leading cause of gynaecological cancer-related mortality worldwide. CC pathogenesis is triggered when human papillomavirus (HPV) inserts into the genome, resulting in tumour suppressor gene inactivation and oncogene activation. Collecting tumour and blood samples is critical for identifying these genetic alterations. METHODS: BIO-RAIDs is the first prospective molecular profiling clinical study to include a substantial biobanking effort that used uniform high-quality standards and control of samples. In this European Union (EU)-funded study, we identified the challenges that were impeding the effective implementation of such a systematic and comprehensive biobanking effort. RESULTS: The challenges included a lack of uniform international legal and ethical standards, complexities in clinical and molecular data management, and difficulties in determining the best technical platforms and data analysis techniques. Some difficulties were encountered by all investigators, while others affected only certain institutions, regions, or countries. CONCLUSIONS: The results of the BIO-RAIDs programme highlight the need to facilitate and standardise regulatory procedures, and we feel that there is also a need for international working groups that make recommendations to regulatory bodies, governmental funding agencies, and academic institutions to achieve a proficient biobanking programme throughout EU countries. This represents the first step in precision medicine.

Emergence and evolution of an international cluster of MDR Bacteroides fragilis isolates.

OBJECTIVES: The aim of this study was to examine the antibiotic resistance profiles, antibiotic resistance mechanisms and possible 'clonal' nature of some MDR Bacteroides fragilis strains that simultaneously harboured cfiA, nimB, IS1186 and IS4351. METHODS: Antibiotic susceptibilities were determined by Etests and antibiotic resistance genes and different genetic elements were detected by applying PCR methods. The environments of the cfiA and nimB genes were also determined by sequencing. The transferability of the cfiA, nimB and tet(Q) genes was tested by conjugation. The genetic relatedness of the test strains was tested by ERIC-PCR or PFGE. The complete genome sequences of two strains (B. fragilis BF8 and O:21) were determined by next-generation sequencing. RESULTS: Most of the seven B. fragilis strains tested displayed multidrug resistance phenotypes; five strains were resistant to at least five types of antibiotics. Besides the common genetic constitution, ERIC-PCR implied high genetic relatedness. Similarities in some of the antibiotic resistance mechanisms [carbapenems (cfiA) and metronidazole (nimB)] also confirmed their common origin, but some other resistance mechanisms {MLSB [erm(F)] and tetracycline [tet(Q)]} and PFGE typing revealed differences. In B. fragilis BF8 and O:21, erm(F) and tet(X) genes were found with IS4351 borders, thus constituting Tn4351. All the strains were tet(Q) positive and transferred this gene in conjugation experiments, but not the cfiA and nimB genes. CONCLUSIONS: An international cluster of MDR B. fragilis strains has been identified and characterized. This 'clone' may have emerged early in the evolution of division II B. fragilis strains, which was suggested by the low-complexity ERIC profiles and differences in the PFGE patterns.

A DNA methylation fingerprint of 1628 human samples.

Most of the studies characterizing DNA methylation patterns have been restricted to particular genomic loci in a limited number of human samples and pathological conditions. Herein, we present a compromise between an extremely comprehensive study of a human sample population with an intermediate level of resolution of CpGs at the genomic level. We obtained a DNA methylation fingerprint of 1628 human samples in which we interrogated 1505 CpG sites. The DNA methylation patterns revealed show this epigenetic mark to be critical in tissue-type definition and stemness, particularly around transcription start sites that are not within a CpG island. For disease, the generated DNA methylation fingerprints show that, during tumorigenesis, human cancer cells underwent a progressive gain of promoter CpG-island hypermethylation and a loss of CpG methylation in non-CpG-island promoters. Although transformed cells are those in which DNA methylation disruption is more obvious, we observed that other common human diseases, such as neurological and autoimmune disorders, had their own distinct DNA methylation profiles. Most importantly, we provide proof of principle that the DNA methylation fingerprints obtained might be useful for translational purposes by showing that we are able to identify the tumor type origin of cancers of unknown primary origin (CUPs). Thus, the DNA methylation patterns identified across the largest spectrum of samples, tissues, and diseases reported to date constitute a baseline for developing higher-resolution DNA methylation maps and provide important clues concerning the contribution of CpG methylation to tissue identity and its changes in the most prevalent human diseases.

DNA methylation map of mouse and human brain identifies target genes in Alzheimer's disease.

The central nervous system has a pattern of gene expression that is closely regulated with respect to functional and anatomical regions. DNA methylation is a major regulator of transcriptional activity, and aberrations in the distribution of this epigenetic mark may be involved in many neurological disorders, such as Alzheimer's disease. Herein, we have analysed 12 distinct mouse brain regions according to their CpG 5'-end gene methylation patterns and observed their unique epigenetic landscapes. The DNA methylomes obtained from the cerebral cortex were used to identify aberrant DNA methylation changes that occurred in two mouse models of Alzheimer's disease. We were able to translate these findings to patients with Alzheimer's disease, identifying DNA methylation-associated silencing of three targets genes: thromboxane A2 receptor (TBXA2R), sorbin and SH3 domain containing 3 (SORBS3) and spectrin beta 4 (SPTBN4). These hypermethylation targets indicate that the cyclic AMP response element-binding protein (CREB) activation pathway and the axon initial segment could contribute to the disease.

ContScout: sensitive detection and removal of contamination from annotated genomes.

Contamination of genomes is an increasingly recognized problem affecting several downstream applications, from comparative evolutionary genomics to metagenomics. Here we introduce ContScout, a precise tool for eliminating foreign sequences from annotated genomes. It achieves high specificity and sensitivity on synthetic benchmark data even when the contaminant is a closely related species, outperforms competing tools, and can distinguish horizontal gene transfer from contamination. A screen of 844 eukaryotic genomes for contamination identified bacteria as the most common source, followed by fungi and plants. Furthermore, we show that contaminants in ancestral genome reconstructions lead to erroneous early origins of genes and inflate gene loss rates, leading to a false notion of complex ancestral genomes. Taken together, we offer here a tool for sensitive removal of foreign proteins, identify and remove contaminants from diverse eukaryotic genomes and evaluate their impact on phylogenomic analyses.