Search

In recent years it became clear that in eukaryotic genome evolution gene loss is prevalent over gene gain. However, the absence of genes in an annotated genome is not always equivalent to the loss of genes. Due to sequencing issues, or incorrect gene prediction, genes can be falsely inferred as absent. This implies that loss estimates are overestimated and, more generally, that falsely inferred absences impact genomic comparative studies. However, reliable estimates of how prevalent this issue is are lacking. Here we quantified the impact of gene prediction on gene loss estimates in eukaryotes by analysing 209 phylogenetically diverse eukaryotic organisms and comparing their predicted proteomes to that of their respective six-frame translated genomes. We observe that 4.61% of domains per species were falsely inferred to be absent for Pfam domains predicted to have been present in the last eukaryotic common ancestor. Between phylogenetically different categories this estimate varies substantially: for clade-specific loss (ancestral loss) we found 1.30% and for species-specific loss 16.88% to be falsely inferred as absent. For BUSCO 1-to-1 orthologous families, 18.30% were falsely inferred to be absent. Finally, we showed that falsely inferred absences indeed impact loss estimates, with the number of losses decreasing by 11.78%. Our work strengthens the increasing number of studies showing that gene loss is an important factor in eukaryotic genome evolution. However, while we demonstrate that on average inferring gene absences from predicted proteomes is reliable, caution is warranted when inferring species-specific absences., Competing Interests: The authors have declared that no competing interests exist.

Comparative genomics has proven a fruitful approach to acquire many functional and evolutionary insights into core cellular processes. Here it is argued that in order to perform accurate and interesting comparative genomics, one first and foremost has to be able to recognize, postulate, and revise different evolutionary scenarios. After all, these studies lack a simple protocol, due to different proteins having different evolutionary dynamics and demanding different approaches. The authors here discuss this challenge from a practical (what are the observations?) and conceptual (how do these indicate a specific evolutionary scenario?) viewpoint, with the aim to guide investigators who want to analyze the evolution of their protein(s) of interest. By sharing how the authors draft, test, and update such a scenario and how it directs their investigations, the authors hope to illuminate how to execute molecular evolution studies and how to interpret them. Also see the video abstract here https://youtu.be/VCt3l2pbdbQ., (© The Authors. Bio Essays Published by WILEY Periodicals, Inc.)

Hotspots of rapid genome evolution hold clues about human adaptation. We present a comparative analysis of nine whole-genome sequenced primates to identify high-confidence targets of positive selection. We find strong statistical evidence for positive selection in 331 protein-coding genes (3%), pinpointing 934 adaptively evolving codons (0.014%). Our new procedure is stringent and reveals substantial artefacts (20% of initial predictions) that have inflated previous estimates. The final 331 positively selected genes (PSG) are strongly enriched for innate and adaptive immunity, secreted and cell membrane proteins (e.g. pattern recognition, complement, cytokines, immune receptors, MHC, Siglecs). We also find evidence for positive selection in reproduction and chromosome segregation (e.g. centromere-associated CENPO, CENPT), apolipoproteins, smell/taste receptors and mitochondrial proteins. Focusing on the virus-host interaction, we retrieve most evolutionary conflicts known to influence antiviral activity (e.g. TRIM5, MAVS, SAMHD1, tetherin) and predict 70 novel cases through integration with virus-human interaction data. Protein structure analysis further identifies positive selection in the interaction interfaces between viruses and their cellular receptors (CD4-HIV; CD46-measles, adenoviruses; CD55-picornaviruses). Finally, primate PSG consistently show high sequence variation in human exomes, suggesting ongoing evolution. Our curated dataset of positive selection is a rich source for studying the genetics underlying human (antiviral) phenotypes. Procedures and data are available at https://github.com/robinvanderlee/positive-selection., (© The Author(s) 2017. Published by Oxford University Press on behalf of Nucleic Acids Research.)

The intraflagellar transport (IFT) complex is an integral component of the cilium, a quintessential organelle of the eukaryotic cell. The IFT system consists of three subcomplexes [i.e., intraflagellar transport (IFT)-A, IFT-B, and the BBSome], which together transport proteins and other molecules along the cilium. IFT dysfunction results in diseases collectively called ciliopathies. It has been proposed that the IFT complexes originated from vesicle coats similar to coat protein complex (COP) I, COPII, and clathrin. Here we provide phylogenetic evidence for common ancestry of IFT subunits and α, β', and ε subunits of COPI, and trace the origins of the IFT-A, IFT-B, and the BBSome subcomplexes. We find that IFT-A and the BBSome likely arose from an IFT-B-like complex by intracomplex subunit duplication. The distribution of IFT proteins across eukaryotes identifies the BBSome as a frequently lost, modular component of the IFT. Significantly, loss of the BBSome from a taxon is a frequent precursor to complete cilium loss in related taxa. Given the inferred late origin of the BBSome in cilium evolution and its frequent loss, the IFT complex behaves as a "last-in, first-out" system. The protocoatomer origin of the IFT complex corroborates involvement of IFT components in vesicle transport. Expansion of IFT subunits by duplication and their subsequent independent loss supports the idea of modularity and structural independence of the IFT subcomplexes.

The TOR kinase is a major regulator of growth in eukaryotes. Many components of the TOR pathway are implicated in cancer and metabolic diseases in humans. Analysis of the evolution of TOR and its pathway may provide fundamental insight into the evolution of growth regulation in eukaryotes and provide a practical framework on which experimental evidence can be compared between species. Here we performed phylogenetic analyses on the components of the TOR pathway and determined their point of invention. We find that the two TOR complexes and a large part of the TOR pathway originated before the Last Eukaryotic Common Ancestor and form a core to which new inputs have been added during animal evolution. In addition, we provide insight into how duplications and sub-functionalization of the S6K, RSK, SGK and PKB kinases shaped the complexity of the TOR pathway. In yeast we identify novel AGC kinases that are orthologous to the S6 kinase. These results demonstrate how a vital signaling pathway can be both highly conserved and flexible in eukaryotes.

The members of the Ras-like superfamily of small GTP-binding proteins are molecular switches that are in general regulated in time and space by guanine nucleotide exchange factors and GTPase activating proteins. The Ras-like G-proteins Ras, Rap and Ral are regulated by a variety of guanine nucleotide exchange factors that are characterized by a CDC25 homology domain. Here we study the evolution of the Ras pathway by determining the evolutionary history of CDC25 homology domain coding sequences. We identified CDC25 homology domain coding sequences in animals, fungi and a wide range of protists, but not in plants. This suggests that the CDC25 homology domain originated in or before the last eukaryotic ancestor but was subsequently lost in plant. We provide evidence that at least seven different ancestral Ras guanine nucleotide exchange factors were present in the ancestor of fungi and animals. Differences between present day fungi and animals are the result of loss of ancestral Ras guanine nucleotide exchange factors early in fungal and animal evolution combined with lineage specific duplications and domain acquisitions. In addition, we identify Ral guanine exchange factors and Ral in early diverged fungi, dating the origin of Ral signaling back to before the divergence of animals and fungi. We conclude that the Ras signaling pathway evolved by gradual change as well as through differential sampling of the ancestral CDC25 homology domain repertoire by both fungi and animals. Finally, a comparison of the domain composition of the Ras guanine nucleotide exchange factors shows that domain addition and diversification occurred both prior to and after the fungal-animal split.

The formation of proteins into stable protein complexes plays a fundamental role in the operation of the cell. The study of the degree of evolutionary conservation of protein complexes between species and the evolution of protein-protein interactions has been hampered by lack of comprehensive coverage of the high-throughput (HTP) technologies that measure the interactome. We show that new high-throughput datasets on protein co-purification in yeast have a substantially lower false negative rate than previous datasets when compared to known complexes. These datasets are therefore more suitable to estimate the conservation of protein complex membership than hitherto possible. We perform comparative genomics between curated protein complexes from human and the HTP data in Saccharomyces cerevisiae to study the evolution of co-complex memberships. This analysis revealed that out of the 5,960 protein pairs that are part of the same complex in human, 2,216 are absent because both proteins lack an ortholog in S. cerevisiae, while for 1,828 the co-complex membership is disrupted because one of the two proteins lacks an ortholog. For the remaining 1,916 protein pairs, only 10% were never co-purified in the large-scale experiments. This implies a conservation level of co-complex membership of 90% when the genes coding for the protein pairs that participate in the same protein complex are also conserved. We conclude that the evolutionary dynamics of protein complexes are, by and large, not the result of network rewiring (i.e. acquisition or loss of co-complex memberships), but mainly due to genomic acquisition or loss of genes coding for subunits. We thus reveal evidence for the tight interrelation of genomic and network evolution.