Comparative Evolutionary Analysis Research Articles

Building a three-dimensional model of rat albumin molecule by homology modeling

Albumin is known as being able to cleave ester bonds in organophosphates (OP) and serve as a major means of OP detoxification due to its large amount in the blood serum. Most of in vivo toxicological studies are conducted on rodents, mice and rats. To adequately extrapolate results of ligand/albumin interaction studies to a human organism, it is necessary to carry out a comparative analysis of the rat (RSA) and human (HSA) albumin structure by in silico methods. X-ray diffraction analysis of RSA has not been done as yet, and the RSA structure is still undetermined. The aim of this study was to build a three-dimensional (3D) model of RSA by homology modeling. As templates, 14 RSA homologs were selected from the Protein Data Bank. A 3D model of the RSA molecule was built using the RSA primary sequence and 3D models of the selected templates. The 3D model geometry was improved by the energy minimization method. The quality of the model was evaluated by 9 parameters. According to this evaluation, in 97% of protein structures from the comparative dataset these parameters were worse than in the obtained model. It was concluded that this model is suitable for further research, specifically, for in silico studies of the RSA/OP interaction and comparative evolutionary analysis of RSA and HSA.

Phylogenomic Synteny Network Analysis of MADS-Box Transcription Factor Genes Reveals Lineage-Specific Transpositions, Ancient Tandem Duplications, and Deep Positional Conservation.

Conserved genomic context provides critical information for comparative evolutionary analysis. With the increase in numbers of sequenced plant genomes, synteny analysis can provide new insights into gene family evolution. Here, we exploit a network analysis approach to organize and interpret massive pairwise syntenic relationships. Specifically, we analyzed synteny networks of the MADS-box transcription factor gene family using 51 completed plant genomes. In combination with phylogenetic profiling, several novel evolutionary patterns were inferred and visualized from synteny network clusters. We found lineage-specific clusters that derive from transposition events for the regulators of floral development (APETALA3 and PI) and flowering time (FLC) in the Brassicales and for the regulators of root development (AGL17) in Poales. We also identified two large gene clusters that jointly encompass many key phenotypic regulatory Type II MADS-box gene clades (SEP1, SQUA, TM8, SEP3, FLC, AGL6, and TM3). Gene clustering and gene trees support the idea that these genes are derived from an ancient tandem gene duplication that likely predates the radiation of the seed plants and then expanded by subsequent polyploidy events. We also identified angiosperm-wide conservation of synteny of several other less studied clades. Combined, these findings provide new hypotheses for the genomic origins, biological conservation, and divergence of MADS-box gene family members.

Comparative genomics reveals contraction in olfactory receptor genes in bats

Gene loss and gain during genome evolution are thought to play important roles in adaptive phenotypic diversification. Among mammals, bats possess the smallest genomes and have evolved the unique abilities of powered flight and laryngeal echolocation. To investigate whether gene family evolution has contributed to the genome downsizing and phenotypic diversification in this group, we performed comparative evolutionary analyses of complete proteome data for eight bat species, including echolocating and non-echolocating forms, together with the proteomes of 12 other laurasiatherian mammals. Our analyses revealed extensive gene loss in the most recent ancestor of bats, and also of carnivores (both >1,000 genes), although this gene contraction did not appear to correlate with the reduction in genome size in bats. Comparisons of highly dynamic families suggested that expansion and contraction affected genes with similar functions (immunity, response to stimulus) in all laurasiatherian lineages. However, the magnitude and direction of these changes varied greatly among groups. In particular, our results showed contraction of the Olfactory Receptor (OR) gene repertoire in the last common ancestor of all bats, as well as that of the echolocating species studied. In contrast, non-echolocating fruit bats showed evidence of expansion in ORs, supporting a “trade-off” between sensory modalities.

Gene expression profiling during the embryo-to-larva transition in the giant red sea urchin Mesocentrotus franciscanus.

In echinoderms, major morphological transitions during early development are attributed to different genetic interactions and changes in global expression patterns that shape the regulatory program for the specification of embryonic territories. In order more thoroughly to understand these biological and molecular processes, we examined the transcriptome structure and expression profiles during the embryo‐to‐larva transition of a keystone species, the giant red sea urchin Mesocentrotus franciscanus. Using a de novo assembly approach, we obtained 176,885 transcripts from which 60,439 (34%) had significant alignments to known proteins. From these transcripts, ~80% were functionally annotated allowing the identification of ~2,600 functional, structural, and regulatory genes involved in developmental process. Analysis of expression profiles between gastrula and pluteus stages of M. franciscanus revealed 791 differentially expressed genes with 251 GO overrepresented terms. For gastrula, up‐regulated GO terms were mainly linked to cell differentiation and signal transduction involved in cell cycle checkpoints. In the pluteus stage, major GO terms were associated with phosphoprotein phosphatase activity, muscle contraction, and olfactory behavior, among others. Our evolutionary comparative analysis revealed that several of these genes and functional pathways are highly conserved among echinoids, holothuroids, and ophiuroids.

Molecular Basis of ABHD5 Lipolysis Activation

Alpha-beta hydrolase domain-containing 5 (ABHD5), the defective gene in human Chanarin-Dorfman syndrome, is a highly conserved regulator of adipose triglyceride lipase (ATGL)-mediated lipolysis that plays important roles in metabolism, tumor progression, viral replication, and skin barrier formation. The structural determinants of ABHD5 lipolysis activation, however, are unknown. We performed comparative evolutionary analysis and structural modeling of ABHD5 and ABHD4, a functionally distinct paralog that diverged from ABHD5 ~500 million years ago, to identify determinants of ABHD5 lipolysis activation. Two highly conserved ABHD5 amino acids (R299 and G328) enabled ABHD4 (ABHD4 N303R/S332G) to activate ATGL in Cos7 cells, brown adipocytes, and artificial lipid droplets. The corresponding ABHD5 mutations (ABHD5 R299N and ABHD5 G328S) selectively disrupted lipolysis without affecting ATGL lipid droplet translocation or ABHD5 interactions with perilipin proteins and ABHD5 ligands, demonstrating that ABHD5 lipase activation could be dissociated from its other functions. Structural modeling placed ABHD5 R299/G328 and R303/G332 from gain-of-function ABHD4 in close proximity on the ABHD protein surface, indicating they form part of a novel functional surface required for lipase activation. These data demonstrate distinct ABHD5 functional properties and provide new insights into the functional evolution of ABHD family members and the structural basis of lipase regulation.

Systematic Analyses and Prediction of Human Drug Side Effect Associated Proteins from the Perspective of Protein Evolution.

Identification of various factors involved in adverse drug reactions in target proteins to develop therapeutic drugs with minimal/no side effect is very important. In this context, we have performed a comparative evolutionary rate analyses between the genes exhibiting drug side-effect(s) (SET) and genes showing no side effect (NSET) with an aim to increase the prediction accuracy of SET/NSET proteins using evolutionary rate determinants. We found that SET proteins are more conserved than the NSET proteins. The rates of evolution between SET and NSET protein primarily depend upon their noncomplex (protein complex association number = 0) forming nature, phylogenetic age, multifunctionality, membrane localization, and transmembrane helix content irrespective of their essentiality, total druggability (total number of drugs/target), m-RNA expression level, and tissue expression breadth. We also introduced two novel terms—killer druggability (number of drugs with killing side effect(s)/target), essential druggability (number of drugs targeting essential proteins/target) to explain the evolutionary rate variation between SET and NSET proteins. Interestingly, we noticed that SET proteins are younger than NSET proteins and multifunctional younger SET proteins are candidates of acquiring killing side effects. We provide evidence that higher killer druggability, multifunctionality, and transmembrane helices support the conservation of SET proteins over NSET proteins in spite of their recent origin. By employing all these entities, our Support Vector Machine model predicts human SET/NSET proteins to a high degree of accuracy (∼86%).

East Polynesian Islands as Models of Cultural Divergence: The Case of Rapa Nui and Rapa Iti

ABSTRACTThe diverse island societies of East Polynesia are well-suited as models for comparative evolutionary analysis. Settled ca. 750 BP by a common ancestral population, colonists of the remote corners of the Pacific shared a pool of cultural traits that included commensal species, language, technology, and other cultural practices. Following colonization however, island populations diverged in language, subsistence practices, degree of territoriality, settlement patterns, investment and forms of monumental architecture, and social organization. Driven by historical circumstances and varied environmental conditions, this divergence presents evolutionary case studies of alternative paths of cultural change. One explanatory approach to this evolutionary divergence involves isolating the critical ecological parameters that likely constrained and shaped the diverse history of island populations. Here, we offer a comparative evolutionary analysis that explores the divergent histories of two marginal East Polynesian islands: Rapa Nui and Rapa Iti.

Comparative evolutionary study on transitions of national innovation systems in the Baltic Region

The paper provides an evolutionary study on the transition of national innovation systems (NIS) within the catching up economy context in the Baltic Region. The regional and institutional context provides an exceptional chance to look at the evolution of economies from planned economy towards free market, and reaching advanced high economy profile over 25 years (International Monetary Fund, 2015). However, the authors argue, that economy transition does not simultaneously lead to the technological upgrading and sophistication of NIS. The comparative evolutionary data analysis in three Baltic States, namely Estonia, Latvia and Lithuania, that reached advanced economy profile in 2015 (IMF) is used to illustrate the statement. The period of study encompasses ten years of membership in European Union (2006-2015), and aims at providing the national innovation system transition trends from the evolutionary perspective. The authors aim to provide the reader with the broad overview of transition patterns taken. Empirical findings demonstrate the differences in national innovation system transition, whereas Estonia demonstrates balanced development of national innovation system, and Lithuania together with Latvia demonstrate asymmetries and stagnating development despite successful economic transition.

Comparative evolutionary study on transitions of national innovation systems in the Baltic Region

The paper provides an evolutionary study on the transition of national innovation systems (NIS) within the catching up economy context in the Baltic Region. The regional and institutional context provides an exceptional chance to look at the evolution of economies from planned economy towards free market, and reaching advanced high economy profile over 25 years (International Monetary Fund, 2015). However, the authors argue, that economy transition does not simultaneously lead to the technological upgrading and sophistication of NIS. The comparative evolutionary data analysis in three Baltic States, namely Estonia, Latvia and Lithuania, that reached advanced economy profile in 2015 (IMF) is used to illustrate the statement. The period of study encompasses ten years of membership in European Union (2006-2015), and aims at providing the national innovation system transition trends from the evolutionary perspective. The authors aim to provide the reader with the broad overview of transition patterns taken. Empirical findings demonstrate the differences in national innovation system transition, whereas Estonia demonstrates balanced development of national innovation system, and Lithuania together with Latvia demonstrate asymmetries and stagnating development despite successful economic transition.

Genetic diversity among eight Dendrolimus species in Eurasia (Lepidoptera: Lasiocampidae) inferred from mitochondrial COI and COII, and nuclear ITS2 markers.

BackgroundMoths of genus Dendrolimus (Lepidoptera: Lasiocampidae) are among the major pests of coniferous forests worldwide. Taxonomy and nomenclature of this genus are not entirely established, and there are many species with a controversial taxonomic position. We present a comparative evolutionary analysis of the most economically important Dendrolimus species in Eurasia.ResultsOur analysis was based on the nucleotide sequences of COI and COII mitochondrial genes and ITS2 spacer of nuclear ribosomal genes. All known sequences were extracted from GenBank. Additional 112 new sequences were identified for 28 specimens of D. sibiricus, D. pini, and D. superans from five regions of Siberia and the Russian Far East to be able to compare the disparate data from all previous studies. In total, 528 sequences were used in phylogenetic analysis. Two clusters of closely related species in Dendrolimus were found. The first cluster includes D. pini, D. sibiricus, and D. superans; and the second, D. spectabilis, D. punctatus, and D. tabulaeformis. Species D. houi and D. kikuchii appear to be the most basal in the genus.ConclusionGenetic difference among the second cluster species is very low in contrast to the first cluster species. Phylogenetic position D. tabulaeformis as a subspecies was supported. It was found that D. sibiricus recently separated from D. superans. Integration of D. sibiricus mitochondrial DNA sequences and the spread of this species to the west of Eurasia have been established as the cause of the unjustified allocation of a new species: D. kilmez. Our study further clarifies taxonomic problems in the genus and gives more complete information on the genetic structure of D. pini, D. sibiricus, and D. superans.Electronic supplementary materialThe online version of this article (doi:10.1186/s12863-016-0463-5) contains supplementary material, which is available to authorized users.

Mining the Cicer arietinum genome for the mildew locus O (Mlo) gene family and comparative evolutionary analysis of the Mlo genes from Medicago truncatula and some other plant species.

The mildew locus O (Mlo) gene family is ubiquitous in land plants. Some members of this gene family are involved in negative regulation of powdery mildew resistance, while others are involved in several other biological functions. Mlo proteins have characteristic seven transmembrane domains and a calmodulin-binding domain at their C-termini, and are associated with plasma membrane. The Mlo gene family has been studied in several economically important cereals, but little information is available on this gene family in the important legumes, Medicago truncatula Gaertn. and Cicer arietinum L. We carried out a comprehensive and comparative investigation of the Mlo gene family in these two species using the genome sequences available at the M. truncatula genome database (Mt v4.0) and NCBI (C. arietinum). A genome-wide homology-based search using Arabidopsis Mlo proteins as query identified 16 MtMlo (M. truncatula Mlo) and 14 CarMlo (C. arietinum Mlo) genes. The MtMlo and CarMlo genes had comparable gene structure, protein sequence and topology. Their chromosomal locations indicated the occurrence of extensive reorganization in the genomes of the two species after their divergence from the common ancestor. A multiple sequence alignment of 53 Mlo proteins from these two and several other species showed a highly conserved sequence block of seven amino acids, viz., L-ETPTW, towards their N-termini. The evolutionary phylogenetic analysis grouped the MtMlo and CarMlo members into four clusters, and most of the MtMlo and CarMlo members formed one-to-one ortholog pairs. The ka/ks analyses indicated that the MtMlo and CarMlo genes are subjected to intense purifying selection.

Aquerium: A web application for comparative exploration of domain-based protein occurrences on the taxonomically clustered genome tree.

Gene duplication and loss are major driving forces in evolution. While many important genomic resources provide information on gene presence, there is a lack of tools giving equal importance to presence and absence information as well as web platforms enabling easy visual comparison of multiple domain-based protein occurrences at once. Here, we present Aquerium, a platform for visualizing genomic presence and absence of biomolecules with a focus on protein domain architectures. The web server offers advanced domain organization querying against the database of pre-computed domains for ∼26,000 organisms and it can be utilized for identification of evolutionary events, such as fusion, disassociation, duplication, and shuffling of protein domains. The tool also allows alternative inputs of custom entries or BLASTP results for visualization. Aquerium will be a useful tool for biologists who perform comparative genomic and evolutionary analyses. The web server is freely accessible at http://aquerium.utk.edu. Proteins 2016; 85:72-77. © 2016 Wiley Periodicals, Inc.

Functional alterations due to amino acid changes and evolutionary comparative analysis of ARPKD and ADPKD genes

A targeted customized sequencing of genes implicated in autosomal recessive polycystic kidney disease (ARPKD) phenotype was performed to identify candidate variants using the Ion torrent PGM next-generation sequencing. The results identified four potential pathogenic variants in PKHD1 gene [c.4870C>T, p.(Arg1624Trp), c.5725C>T, p.(Arg1909Trp), c.1736C>T, p.(Thr579Met) and c.10628T>G, p.(Leu3543Trp)] among 12 out of 18 samples. However, one variant c.4870C>T, p.(Arg1624Trp) was common among eight patients. Some patient samples also showed few variants in autosomal dominant polycystic kidney disease (ADPKD) disease causing genes PKD1 and PKD2 such as c.12433G>A, p.(Val4145Ile) and c.1445T>G, p.(Phe482Cys), respectively. All causative variants were validated by capillary sequencing and confirmed the presence of a novel homozygous variant c.10628T>G, p.(Leu3543Trp) in a male proband. We have recently published the results of these studies (Edrees et al., 2016). Here we report for the first time the effect of the common mutation p.(Arg1624Trp) found in eight samples on the protein structure and function due to the specific amino acid changes of PKHD1 protein using molecular dynamics simulations. The computational approaches provide tool predict the phenotypic effect of variant on the structure and function of the altered protein. The structural analysis with the common mutation p.(Arg1624Trp) in the native and mutant modeled protein were also studied for solvent accessibility, secondary structure and stabilizing residues to find out the stability of the protein between wild type and mutant forms. Furthermore, comparative genomics and evolutionary analyses of variants observed in PKHD1, PKD1, and PKD2 genes were also performed in some mammalian species including human to understand the complexity of genomes among closely related mammalian species. Taken together, the results revealed that the evolutionary comparative analyses and characterization of PKHD1, PKD1, and PKD2 genes among various related and unrelated mammalian species will provide important insights into their evolutionary process and understanding for further disease characterization and management.

Comparative genome and evolutionary analysis of naturally occurring Beilong virus in brown and black rats

Recently, we reported the presence of Beilong virus in spleen and kidney samples of brown rats and black rats, suggesting that these rodents could be natural reservoirs of Beilong virus. In this study, four genomes of Beilong virus from brown rats and black rats were sequenced. Similar to the Beilong virus genome sequenced from kidney mesangial cell line culture, those of J-virus from house mouse and Tailam virus from Sikkim rats, these four genomes from naturally occurring Beilong virus also contain the eight genes (3′-N-P/V/C-M-F-SH-TM-G-L-5′). In these four genomes, the attachment glycoprotein encoded by the G gene consists of 1046 amino acids; but for the original Beilong virus genome sequenced from kidney mesangial cell line, the G CDS was predicted to be prematurely terminated at position 2205 (TGG→TAG), resulting in a 734-amino-acid truncated G protein. This phenomenon of a lack of nonsense mutation in naturally occurring Beilong viruses was confirmed by sequencing this region of 15 additional rodent samples. Phylogenetic analyses showed that the cell line and naturally occurring Beilong viruses were closely clustered, without separation into subgroups. In addition, these viruses were further clustered with J-virus and Tailam virus, with high bootstrap supports of >90%, forming a distinct group in Paramyxoviridae. Brown rats and black rats are natural reservoirs of Beilong virus. Our results also supports that the recently proposed genus, Jeilongvirus, should encompass Beilong virus, J-virus and Tailam virus as members.

Retinal transcriptome sequencing sheds light on the adaptation to nocturnal and diurnal lifestyles in raptors.

Owls (Strigiformes) represent a fascinating group of birds that are the ecological night-time counterparts to diurnal raptors (Accipitriformes). The nocturnality of owls, unusual within birds, has favored an exceptional visual system that is highly tuned for hunting at night, yet the molecular basis for this adaptation is lacking. Here, using a comparative evolutionary analysis of 120 vision genes obtained by retinal transcriptome sequencing, we found strong positive selection for low-light vision genes in owls, which contributes to their remarkable nocturnal vision. Not surprisingly, we detected gene loss of the violet/ultraviolet-sensitive opsin (SWS1) in all owls we studied, but two other color vision genes, the red-sensitive LWS and the blue-sensitive SWS2, were found to be under strong positive selection, which may be linked to the spectral tunings of these genes toward maximizing photon absorption in crepuscular conditions. We also detected the only other positively selected genes associated with motion detection in falcons and positively selected genes associated with bright-light vision and eye protection in other diurnal raptors (Accipitriformes). Our results suggest the adaptive evolution of vision genes reflect differentiated activity time and distinct hunting behaviors.

Visions of the past and dreams of the future in the Orient: the Irano-Turanian region from classical botany to evolutionary studies.

Ever since the 19th century, the immense arid lands of the Orient, now called the Irano-Turanian (IT) floristic region, attracted the interest of European naturalists with their tremendous plant biodiversity. Covering approximately 30% of the surface of Eurasia (16000000 km2 ), the IT region is one of the largest floristic regions of the world. The IT region represents one of the hotspots of evolutionary and biological diversity in the Old World, and serves as a source of xerophytic taxa for neighbouring regions. Moreover, it is the cradle of the numerous species domesticated in the Fertile Crescent. Over the last 200 years, naturalists outlined different borders for the IT region. Yet, the delimitation and evolutionary history of this area remain one of the least well-understood fields of global biogeography, even though it is crucial to explaining the distribution of life in Eurasia. No comprehensive review of the biogeographical delimitations nor of the role of geological and climatic changes in the evolution of the IT region is currently available. After considering the key role of floristic regions in biogeography, we review the history of evolving concepts about the borders and composition of the IT region over the past 200 years and outline a tentative circumscription for it. We also summarise current knowledge on the geological and climatic history of the IT region. We then use this knowledge to generate specific evolutionary hypotheses to explain how different geological, palaeoclimatic, and ecological factors contributed to range expansion and contraction, thus shaping patterns of speciation in the IT region over time and space. Both historical and ecological biogeography should be applied to understand better the floristic diversification of the region. This will ultimately require evolutionary comparative analyses based on integrative phylogenetic, geological, climatic, ecological, and species distribution studies on the region. Furthermore, an understanding of evolutionary and ecological processes will play a major role in regional planning for protecting biodiversity of the IT region in facing climatic change. With this review, we aim to introduce the IT floristic region to a broader audience of evolutionary, ecological and systematic biologists, thus promoting cutting-edge research on this area and raising awareness of this vast and diverse, yet understudied, part of the world.

Bioinformatics Approach in Plant Genomic Research.

The advance in genomics technology leads to the dramatic change in plant biology research. Plant biologists now easily access to enormous genomic data to deeply study plant high-density genetic variation at molecular level. Therefore, fully understanding and well manipulating bioinformatics tools to manage and analyze these data are essential in current plant genome research. Many plant genome databases have been established and continued expanding recently. Meanwhile, analytical methods based on bioinformatics are also well developed in many aspects of plant genomic research including comparative genomic analysis, phylogenomics and evolutionary analysis, and genome-wide association study. However, constantly upgrading in computational infrastructures, such as high capacity data storage and high performing analysis software, is the real challenge for plant genome research. This review paper focuses on challenges and opportunities which knowledge and skills in bioinformatics can bring to plant scientists in present plant genomics era as well as future aspects in critical need for effective tools to facilitate the translation of knowledge from new sequencing data to enhancement of plant productivity.

Behavioral responses to injury and death in wild Barbary macaques (Macaca sylvanus).

The wounding or death of a conspecific has been shown to elicit varied behavioral responses throughout thanatology. Recently, a number of reports have presented contentious evidence of epimeletic behavior towards the dying and dead among non-human animals, a behavioral trait previously considered uniquely human. Here, we report on the behavioral responses of Barbary macaques, a social, non-human primate, to the deaths of four group members (one high-ranking adult female, one high-ranking adult male, one juvenile male, and one female infant), all caused by road traffic accidents. Responses appeared to vary based on the nature of the death (protracted or instant) and the age class of the deceased. Responses included several behaviors with potential adaptive explanations or consequences. These included exploration, caretaking (guarding, carrying, and grooming), and proximity to wounded individuals or corpses, and immediate as well as longer-lasting distress behaviors from other group members following death, all of which have been reported in other non-human primate species. These observations add to a growing body of comparative evolutionary analysis of primate thanatology and help to highlight the multifaceted impacts of human-induced fatalities on an endangered and socially complex primate.

A high-quality carrot genome assembly provides new insights into carotenoid accumulation and asterid genome evolution.

We report a high-quality chromosome-scale assembly and analysis of the carrot (Daucus carota) genome, the first sequenced genome to include a comparative evolutionary analysis among members of the euasterid II clade. We characterized two new polyploidization events, both occurring after the divergence of carrot from members of the Asterales order, clarifying the evolutionary scenario before and after radiation of the two main asterid clades. Large- and small-scale lineage-specific duplications have contributed to the expansion of gene families, including those with roles in flowering time, defense response, flavor, and pigment accumulation. We identified a candidate gene, DCAR_032551, that conditions carotenoid accumulation (Y) in carrot taproot and is coexpressed with several isoprenoid biosynthetic genes. The primary mechanism regulating carotenoid accumulation in carrot taproot is not at the biosynthetic level. We hypothesize that DCAR_032551 regulates upstream photosystem development and functional processes, including photomorphogenesis and root de-etiolation.

The Genome of Haemoproteus tartakovskyi and Its Relationship to Human Malaria Parasites.

The phylogenetic relationships among hemosporidian parasites, including the origin of Plasmodium falciparum, the most virulent malaria parasite of humans, have been heavily debated for decades. Studies based on multiple-gene sequences have helped settle many of these controversial phylogenetic issues. However, denser taxon sampling and genome-wide analyses are needed to confidently resolve the evolutionay relationships among hemosporidian parasites. Genome sequences of several Plasmodium parasites are available but only for species infecting primates and rodents. To root the phylogenetic tree of Plasmodium, genomic data from related parasites of birds or reptiles are required. Here, we use a novel approach to isolate parasite DNA from microgametes and describe the first genome of a bird parasite in the sister genus to Plasmodium, Haemoproteus tartakovskyi. Similar to Plasmodium parasites, H. tartakovskyi has a small genome (23.2 Mb, 5,990 genes) and a GC content (25.4%) closer to P. falciparum (19.3%) than to Plasmodium vivax (42.3%). Combined with novel transcriptome sequences of the bird parasite Plasmodium ashfordi, our phylogenomic analyses of 1,302 orthologous genes demonstrate that mammalian-infecting malaria parasites are monophyletic, thus rejecting the repeatedly proposed hypothesis that the ancestor of Laverania parasites originated from a secondary host shift from birds to humans. Genes and genomic features previously found to be shared between P. falciparum and bird malaria parasites, but absent in other mammal malaria parasites, are therefore signatures of maintained ancestral states. We foresee that the genome of H. tartakovskyi will open new directions for comparative evolutionary analyses of malarial adaptive traits.