Phylogeny-based Analysis Research Articles

Revision of the full stereochemistry of telomycin

The full stereochemistry of telomycin, an exciting nonribosomal cyclic depsipeptide antibiotic, has been corrected by phylogeny-based bioinformatic analysis of its nonribosomal peptide synthetases sequence and rigorous chiral amino acid analysis.

Genome size and GC content of myxomycetes

More than 1272 myxomycetes species have been described, accounting for more than half of all Amoebozoa species. However, the genome size of only three myxomycetes species has been reported. Therefore, we used flow cytometry to present an extensive survey and a phylogeny-based analysis of genome size and GC content evolution in 144 myxomycetes species. The genome size of myxomycetes ranged from 18.7 Mb to 470.3 Mb, and the GC content ranged from 38.7% to 70.1%. Bright-spored clade showed larger genome sizes and more intra-order genome size variations than the dark-spored clade. GC content and genome size were positively correlated in both bright-spored and dark-spored clades, and spore size was positively correlated with genome size and GC content in the bright-spored clade. We provided the first genome size data set in Myxomycetes, and our results will provide helpful information for future Myxomycetes studies, such as genome sequencing.

PhyloMed: a phylogeny-based test of mediation effect in microbiome

Microbiome data from sequencing experiments contain the relative abundance of a large number of microbial taxa with their evolutionary relationships represented by a phylogenetic tree. The compositional and high-dimensional nature of the microbiome mediator challenges the validity of standard mediation analyses. We propose a phylogeny-based mediation analysis method called PhyloMed to address this challenge. Unlike existing methods that directly identify individual mediating taxa, PhyloMed discovers mediation signals by analyzing subcompositions defined on the phylogenic tree. PhyloMed produces well-calibrated mediation test p-values and yields substantially higher discovery power than existing methods.

Hooking the scientific community on thorny-headed worms: interesting and exciting facts, knowledge gaps and perspectives for research directions on Acanthocephala.

Although interest in Acanthocephala seems to have reached only a small community of researchers worldwide, we show in this opinion article that this group of parasites is composed of excellent model organisms for studying key questions in parasite molecular biology and cytogenetics, evolutionary ecology, and ecotoxicology. Their shared ancestry with free-living rotifers makes them an ideal group to explore the origins of the parasitic lifestyle and evolutionary drivers of host shifts and environmental transitions. They also provide useful features in the quest to decipher the proximate mechanisms of parasite-induced phenotypic alterations and better understand the evolution of behavioral manipulation. From an applied perspective, acanthocephalans' ability to accumulate contaminants offers useful opportunities to monitor the impacts - and evaluate the possible mitigation - of anthropogenic pollutants on aquatic fauna and develop the environmental parasitology framework. However, exploring these exciting research avenues will require connecting fragmentary knowledge by enlarging the taxonomic coverage of molecular and phenotypic data. In this opinion paper, we highlight the needs and opportunities of research on Acanthocephala in three main directions: (i) integrative taxonomy (including non-molecular tools) and phylogeny-based comparative analysis; (ii) ecology and evolution of life cycles, transmission strategies and host ranges; and (iii) environmental issues related to global changes, including ecotoxicology. In each section, the most promising ideas and developments are presented based on selected case studies, with the goal that the present and future generations of parasitologists further explore and increase knowledge of Acanthocephala.

Identification of diverse novel genomoviruses in gut of wild birds

Genomoviridae is a virus family belonging to circular replication associated protein encoding ssDNA (CRESS-DNA) viruses, which have diverse genomic architecture and are widely distributed among different ecosystems. In this study, we characterized 39 novel genomovirus genomes including 3 from wild egrets, 9 from wild cranes, as well as 27 from wild finches in three different types of cloacal swabs of wild bird sampled in the station of Xinqing bird ringing, Maoershan Nature Reserve in Heilongjiang Province, Changbai Mountain in Jilin Province, and Hangzhou Wetland Park in Zhejiang Province, China. Here, using a Rep sequence phylogeny-based analysis, we divided the 39 genomoviruses into 8 genera within the family Genomoviridae, including Gemycircularvirus (n=20), Gemykibivirus (n=3), Gemygorvirus (n=2), Gemyvongvirus (n=2), Gemykolovirus (n=3), Gemykrogvirus (n=6), Gemytondvirus (n=1), Gemyduguivirus (n=1) and one unclassified genomovirus. The 39 genomovirus genomes belong to 36 species (27 of which are new) based on the currently accepted genomovirus pairwise nucleotide sequence identity species demarcation threshold of 78%. Overall, the research enriches our knowledge of CRESS-DNA viral diversity in China and emphasizes the prevalence of genomoviruses in nature.

Functional morphology and post-larval development of the buccal complex in Eubranchus rupium (Nudibranchia: Aeolidida: Fionidae)

Nudibranch molluscs represent an interesting model group to study the evolution of feeding apparatus and feeding modes, being characterized by specialized buccal complex in combination with extremely diverse dietary preferences and multiply prey shifts in evolutionary history. However, the plasticity of the buccal complex morphology in response to diet and specific feeding modes remains understudied. Here we study the general morphology and ontogenesis of the buccal complex in Eubranchus rupium (Nudibranchia: Fionidae). Specific goals are to provide a detailed description of buccal structures morphology in post-larval stages, suggest the feeding mechanism and discuss the phylogenetic value of the morphological characteristics of buccal armature within the genus Eubranchus. Methods included in vivo observations of the feeding process for E. rupium, light microscopic methods, scanning electron microscopy (SEM), transmission electron microscopy (TEM) and confocal laser scanning microscopy (cLSM). According to our results, E. rupium is a mechanical driller, boring holes in hydrozoan perisarc and sucking internal content. The mechanical drilling is supplied by functionally uniserial radula with plate-like laterals teeth of exclusively supportive function and by massive buccal musculature. Comparative phylogeny-based analysis suggests that the drilling feeding mode is common for the genus Eubranchus and indicates radular characters may have a high phylogenetic signal. The buccal complex morphology and feeding mode were found to be similar in both adults and post-metamorphic specimens, its general structures occur even in settled veligers. Juveniles and adults compete for food source, but the competition is smoothed due to characteristics of prey species growth and life cycle.

Predicting the chemical space of fungal polyketides by phylogeny-based bioinformatics analysis of polyketide synthase-nonribosomal peptide synthetase and its modification enzymes

Fungal polyketide synthase (PKS)–nonribosomal peptide synthetase (NRPS) hybrids are key enzymes for synthesizing structurally diverse hybrid natural products (NPs) with characteristic biological activities. Predicting their chemical space is of particular importance in the field of natural product chemistry. However, the unexplored programming rule of the PKS module has prevented prediction of its chemical structure based on amino acid sequences. Here, we conducted a phylogenetic analysis of 884 PKS–NRPS hybrids and a modification enzyme analysis of the corresponding biosynthetic gene cluster, revealing a hidden relationship between its genealogy and core structures. This unexpected result allowed us to predict 18 biosynthetic gene cluster (BGC) groups producing known carbon skeletons (number of BGCs; 489) and 11 uncharacterized BGC groups (171). The limited number of carbon skeletons suggests that fungi tend to select PK skeletons for survival during their evolution. The possible involvement of a horizontal gene transfer event leading to the diverse distribution of PKS–NRPS genes among fungal species is also proposed. This study provides insight into the chemical space of fungal PKs and the distribution of their biosynthetic gene clusters.

Description of soybean aphid (Aphis glycines Matsumura) mitochondrial genome and comparative mitogenomics of Aphididae (Hemiptera: Sternorrhyncha)

The complete mitochondrial genome of the soybean aphid (Aphis glycines Matsumura), a major agricultural pest in the world, is described for the first time, which consists of 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, as well as a large repeat region between tRNA-Glu and tRNA-Phe, and an AT-rich control region. The 17,954 bp mtgenome is the largest reported from the family Aphididae, and its gene order follows the ancestral insect mtgenome except for the repeat region, which contains a 195 bp unit repeated 11.9 times, representing the highest reported repeats among the known aphid mtgenomes to date. A new molecular phylogeny of Aphidae is reconstructed based on all available aphid mtgenomes, and it is shown that the mtgenome data can robustly resolve relationships at the subfamily level, but do not have sufficient phylogenetic information to resolve deep relationships. A phylogeny-based comparative analysis of mtgenomes has been performed to investigate the evolution of the repeat region between tRNA-Glu and tRNA-Phe. So far, among aphids, 13 species are known to have this repeat region of variable lengths, and a phylogenetic analysis of the repeat region shows that a large proportion of the sequences are conserved across the phylogeny, suggesting that the repeat region evolved in the most recent common ancestor of Aphidinae and Eriosomatinae, and that it has gone through numerous episodes of lineage-specific losses and expansions. Combined together, this study provides novel insights into how the repeat regions have evolved within aphids.

Estimating Age-Dependent Extinction: Contrasting Evidence from Fossils and Phylogenies.

The estimation of diversification rates is one of the most vividly debated topics in modern systematics, with considerable controversy surrounding the power of phylogenetic and fossil-based approaches in estimating extinction. Van Valen’s seminal work from 1973 proposed the “Law of constant extinction,” which states that the probability of extinction of taxa is not dependent on their age. This assumption of age-independent extinction has prevailed for decades with its assessment based on survivorship curves, which, however, do not directly account for the incompleteness of the fossil record, and have rarely been applied at the species level. Here, we present a Bayesian framework to estimate extinction rates from the fossil record accounting for age-dependent extinction (ADE). Our approach, unlike previous implementations, explicitly models unobserved species and accounts for the effects of fossil preservation on the observed longevity of sampled lineages. We assess the performance and robustness of our method through extensive simulations and apply it to a fossil data set of terrestrial Carnivora spanning the past 40 myr. We find strong evidence of ADE, as we detect the extinction rate to be highest in young species and declining with increasing species age. For comparison, we apply a recently developed analogous ADE model to a dated phylogeny of extant Carnivora. Although the phylogeny-based analysis also infers ADE, it indicates that the extinction rate, instead, increases with increasing taxon age. The estimated mean species longevity also differs substantially, with the fossil-based analyses estimating 2.0 myr, in contrast to 9.8 myr derived from the phylogeny-based inference. Scrutinizing these discrepancies, we find that both fossil and phylogeny-based ADE models are prone to high error rates when speciation and extinction rates increase or decrease through time. However, analyses of simulated and empirical data show that fossil-based inferences are more robust. This study shows that an accurate estimation of ADE from incomplete fossil data is possible when the effects of preservation are jointly modeled, thus allowing for a reassessment of Van Valen’s model as a general rule in macroevolution.

Sequence-based taxonomic framework for the classification of uncultured single-stranded DNA viruses of the family Genomoviridae.

With the advent of metagenomics approaches, a large diversity of known and unknown viruses has been identified in various types of environmental, plant, and animal samples. One such widespread virus group is the recently established family Genomoviridae which includes viruses with small (∼2–2.4 kb), circular ssDNA genomes encoding rolling-circle replication initiation proteins (Rep) and unique capsid proteins. Here, we propose a sequence-based taxonomic framework for classification of 121 new virus genomes within this family. Genomoviruses display ∼47% sequence diversity, which is very similar to that within the well-established and extensively studied family Geminiviridae (46% diversity). Based on our analysis, we establish a 78% genome-wide pairwise identity as a species demarcation threshold. Furthermore, using a Rep sequence phylogeny-based analysis coupled with the current knowledge on the classification of geminiviruses, we establish nine genera within the Genomoviridae family. These are Gemycircularvirus (n = 73), Gemyduguivirus (n = 1), Gemygorvirus (n = 9), Gemykibivirus (n = 29), Gemykolovirus (n = 3), Gemykrogvirus (n = 3), Gemykroznavirus (n = 1), Gemytondvirus (n = 1), Gemyvongvirus (n = 1). The presented taxonomic framework offers rational classification of genomoviruses based on the sequence information alone and sets an example for future classification of other groups of uncultured viruses discovered using metagenomics approaches.

Comparative analysis reveals migratory swallows (Hirundinidae) have less pointed wings than residents

The correlation between migration and wing pointedness, a pattern generally attributed to the need to reduce drag during powered straight-line flight, is well established in avian ecomorphology. However, most studies investigating this pattern have focused on relatively rounded-wing taxa, which employ different flight modes during foraging and migration. The basic assumption that migrants have comparatively pointed wings has not been questioned by analysing taxa with similar migratory and foraging flight. This study examines the correlation between migration and wing shape in swallows (Hirundinidae), a family with relatively pointed wings in which foraging flight resembles migratory flight. Using a phylogeny-based analysis, we compare the wing shape of species pairs with varying migratory habits in eight swallow genera. Surprisingly, migratory swallows have less pointed wings than sedentary species, and wing pointedness declines linearly with increasing migratory distance. This study represents the first published result documenting a reversal of the correlation between migration and wing pointedness found in other avian taxa. Interpreting this reversal requires a more nuanced understanding of wing ecomorphology; we hypothesize that inclement weather conditions on the breeding grounds and/or the roosting habits of migratory swallows create conflicting selective pressures that increase the cost of wing pointedness in migratory swallows.

Evolution of genome size and chromosome number in the carnivorous plant genus Genlisea (Lentibulariaceae), with a new estimate of the minimum genome size in angiosperms.

Some species of Genlisea possess ultrasmall nuclear genomes, the smallest known among angiosperms, and some have been found to have chromosomes of diminutive size, which may explain why chromosome numbers and karyotypes are not known for the majority of species of the genus. However, other members of the genus do not possess ultrasmall genomes, nor do most taxa studied in related genera of the family or order. This study therefore examined the evolution of genome sizes and chromosome numbers in Genlisea in a phylogenetic context. The correlations of genome size with chromosome number and size, with the phylogeny of the group and with growth forms and habitats were also examined. Nuclear genome sizes were measured from cultivated plant material for a comprehensive sampling of taxa, including nearly half of all species of Genlisea and representing all major lineages. Flow cytometric measurements were conducted in parallel in two laboratories in order to compare the consistency of different methods and controls. Chromosome counts were performed for the majority of taxa, comparing different staining techniques for the ultrasmall chromosomes. Genome sizes of 15 taxa of Genlisea are presented and interpreted in a phylogenetic context. A high degree of congruence was found between genome size distribution and the major phylogenetic lineages. Ultrasmall genomes with 1C values of <100 Mbp were almost exclusively found in a derived lineage of South American species. The ancestral haploid chromosome number was inferred to be n = 8. Chromosome numbers in Genlisea ranged from 2n = 2x = 16 to 2n = 4x = 32. Ascendant dysploid series (2n = 36, 38) are documented for three derived taxa. The different ploidy levels corresponded to the two subgenera, but were not directly correlated to differences in genome size; the three different karyotype ranges mirrored the different sections of the genus. The smallest known plant genomes were not found in G. margaretae, as previously reported, but in G. tuberosa (1C ≈ 61 Mbp) and some strains of G. aurea (1C ≈ 64 Mbp). Genlisea is an ideal candidate model organism for the understanding of genome reduction as the genus includes species with both relatively large (∼1700 Mbp) and ultrasmall (∼61 Mbp) genomes. This comparative, phylogeny-based analysis of genome sizes and karyotypes in Genlisea provides essential data for selection of suitable species for comparative whole-genome analyses, as well as for further studies on both the molecular and cytogenetic basis of genome reduction in plants.

Inter- and intra-specific pan-genomes of Borrelia burgdorferi sensu lato: genome stability and adaptive radiation

BackgroundLyme disease is caused by spirochete bacteria from the Borrelia burgdorferi sensu lato (B. burgdorferi s.l.) species complex. To reconstruct the evolution of B. burgdorferi s.l. and identify the genomic basis of its human virulence, we compared the genomes of 23 B. burgdorferi s.l. isolates from Europe and the United States, including B. burgdorferi sensu stricto (B. burgdorferi s.s., 14 isolates), B. afzelii (2), B. garinii (2), B. “bavariensis” (1), B. spielmanii (1), B. valaisiana (1), B. bissettii (1), and B. “finlandensis” (1).ResultsRobust B. burgdorferi s.s. and B. burgdorferi s.l. phylogenies were obtained using genome-wide single-nucleotide polymorphisms, despite recombination. Phylogeny-based pan-genome analysis showed that the rate of gene acquisition was higher between species than within species, suggesting adaptive speciation. Strong positive natural selection drives the sequence evolution of lipoproteins, including chromosomally-encoded genes 0102 and 0404, cp26-encoded ospC and b08, and lp54-encoded dbpA, a07, a22, a33, a53, a65. Computer simulations predicted rapid adaptive radiation of genomic groups as population size increases.ConclusionsIntra- and inter-specific pan-genome sizes of B. burgdorferi s.l. expand linearly with phylogenetic diversity. Yet gene-acquisition rates in B. burgdorferi s.l. are among the lowest in bacterial pathogens, resulting in high genome stability and few lineage-specific genes. Genome adaptation of B. burgdorferi s.l. is driven predominantly by copy-number and sequence variations of lipoprotein genes. New genomic groups are likely to emerge if the current trend of B. burgdorferi s.l. population expansion continues.

Molecular evolution of the pigmentation gene melanocortin-1 receptor in rodents

Adaptive variation in the melanocortin 1-receptor gene (MC1R), a key locus in melanogenesis, has been identified in some species of rodents. However, in others, MC1R has no causative role in pigmentation phenotypes despite their coat color variation. In this study, we characterized the rates and patterns of MC1R nucleotide and amino acid sequence evolution and, particularly, selective pressures in the separated domains of the protein using a comparative analysis of 43 species representing three major lineages of rodents with variable coat colors. We found high amino acid variation (44% of sites) throughout the protein. Most substitutions were observed in extracellular and transmembrane domains; the intracellular segment was conserved across species. Pairwise non-synonymous substitutions did not vary significantly in different domains among the rodent lineages - i.e., variation was not associated with phylogenetic distance. Phylogeny-based likelihood analysis suggested that purifying selection has mostly shaped the evolutionary course of MC1R. However, a high proportion of sites (27%) were under relaxation of functional constraints (ω = 0.38), and four sites (3, 14, 26, and 251) clearly evolved under positive selection (ω ≅ 2.9). Thus, our data indicate a high proportion of sites evolving under relaxed evolutionary constraints, which might indicate the evolvability of the system in the generation of adaptive changes in specific taxa in rodent lineages.

Evolution of bract development and B‐class MADS box gene expression in petaloid bracts of Cornus s. l. (Cornaceae)

Despite increasing interest in the molecular mechanisms of floral diversity, few studies have investigated the developmental and genetic bases of petaloid bracts. This study examined morphological patterns of bract initiation and expression patterns of B-class MADS-box genes in bracts of several Cornus species. We suggest that petaloid bracts in this genus may not share a single evolutionary origin. Developmental pathways of bracts and spatiotemporal expression of B-class genes in bracts and flowers were examined for four closely related dogwood species. Divergent morphological progressions and gene expression patterns were found in the two sister lineages with petaloid bracts, represented by Cornus florida and Cornus canadensis. Phylogeny-based analysis identified developmental and gene expression changes that are correlated with the evolution of petaloid bracts in C. florida and C. canadensis. Our data support the existence of independent evolutionary origins of petaloid bracts in C. canadensis and C. florida. Additionally, we suggest that functional transference within B-class gene families may have contributed to the origin of bract petaloidy in C. florida. However, the underlying mechanisms of petaloid bract development likely differ between C. florida and C. canadensis. In the future this hypothesis can be tested by functional analyses of Cornus B-class genes.

Genome-wide screening reveals the emergence and divergence of RTK homologues in basal Metazoan Hydra magnipapillata

Receptor tyrosine kinases (RTKs) are key components of cell-cell signalling required for growth and development of multicellular organisms. It is therefore likely that the divergence of RTKs and associated components played a significant role in the evolution of multicellular organisms. We have carried out the present study in hydra, a diploblast, to investigate the divergence of RTKs after parazoa and before emergence of triploblast phyla. The domain-based screening using Hidden Markov Models (HMMs) for RTKs in Genomescan predicted gene models of the Hydra magnipapillata genome resulted in identification of 15 RTKs. These RTKs have been classified into eight families based on domain architecture and homology. Only 5 of these RTKs have been previously reported and a few of these have been partially characterized. A phylogeny-based analysis of these predicted RTKs revealed that seven subtype duplications occurred between 'parazoan-eumetazoan split' and 'diploblast-triploblast split' in animal phyla. These results suggest that most of the RTKs evolved before the radiata-bilateria divergence during animal evolution.

Phylogeny‐based developmental analyses illuminate evolution of inflorescence architectures in dogwoods (Cornus s. l., Cornaceae)

• Inflorescence architecture is important to angiosperm reproduction, but our knowledge of the developmental basis underlying the evolution of inflorescence architectures is limited. Using a phylogeny-based comparative analysis of developmental pathways, we tested the long-standing hypothesis that umbel evolved from elongated inflorescences by suppression of inflorescence branches, while head evolved from umbels by suppression of pedicels. • The developmental pathways of six species of Cornus producing different inflorescence types were characterized by scanning electron microscopy (SEM) and histological analysis. Critical developmental events were traced over the molecular phylogeny to identify evolutionary changes leading to the formation of umbels and heads using methods accounting for evolutionary time and phylogenetic uncertainty. • We defined 24 developmental events describing the developmental progression of the different inflorescence types. The evolutionary transition from paniculate cymes to umbels and heads required alterations of seven developmental events occurring at different evolutionary times. • Our results indicate that heads and umbels evolved independently in Cornus from elongated forms via an umbellate dichasium ancestor and this process involved several independent changes. Our findings shed novel insights into head and umbel evolution concealed by outer morphology. Our work illustrates the importance of combining developmental and phylogenetic data to better define morphological evolutionary processes.

A Comparative Analysis of the Morphology and Evolution of Permanent Sperm Depletion in Spiders

Once thought to be energetically cheap and easy to produce, empirical work has shown that sperm is a costly and limited resource for males. In some spider species, there is behavioral evidence that sperm are permanently depleted after a single mating. This extreme degree of mating investment appears to co-occur with other reproductive strategies common to spiders, e.g. genital mutilation and sexual cannibalism. Here we corroborate that sperm depletion in the golden orb-web spider Nephila clavipes is permanent by uncovering its mechanistic basis using light and electron microscopy. In addition, we use a phylogeny-based statistical analysis to test the evolutionary relationships between permanent sperm depletion (PSD) and other reproductive strategies in spiders. Male testes do not produce sperm during adulthood, which is unusual in spiders. Instead, spermatogenesis is nearly synchronous and ends before the maturation molt. Testis size decreases as males approach their maturation molt and reaches its lowest point after sperm is transferred into the male copulatory organs (pedipalps). As a consequence, the amount of sperm available to males for mating is limited to the sperm contained in the pedipalps, and once it is used, males lose their ability to fertilize eggs. Our data suggest that PSD has evolved independently at least three times within web-building spiders and is significantly correlated with the evolution of other mating strategies that limit males to monogamy, including genital mutilation and sexual cannibalism. We conclude that PSD may be an energy-saving adaptation in species where males are limited to monogamy. This could be particularly important in web-building spiders where extreme sexual size dimorphism results in large, sedentary females and small, searching males who rarely feed as adults and are vulnerable to starvation. Future work will explore possible energetic benefits and the evolutionary lability of PSD relative to other mate-limiting reproductive behaviors.

A New Family of Ty1-copia-Like Retrotransposons Originated in the Tomato Genome by a Recent Horizontal Transfer Event

Rider is a novel and recently active Ty1-copia-like retrotransposon isolated from the T3238fer mutant of tomato. Structurally, it is delimited by a duplication of target sites and contains two long terminal direct repeats and an internal open reading frame, which encodes a Ty1-copia-type polyprotein with characteristic protein domains required for retrotransposition. The family of Rider elements has an intermediate copy number and is scattered in the chromosomes of tomato. Rider family members in the tomato genome share high sequence similarity, but different structural groups were identified (full-size elements, deletion derivatives, and solo LTRs). Southern blot analysis in Solanaceae species showed that Rider was a Lycopersicon-specific element. Sequence analysis revealed that among other plants, two Arabidopsis elements (named as Rider-like 1 and Rider-like 2) are most similar to Rider in both the coding and noncoding regions. RT-PCR analysis indicates that Rider is constitutively expressed in tomato plants. The phylogeny-based parsimony analysis and the sequence substitution analyses of these data suggest that these Rider-like elements originated from a recent introgression of Rider into the tomato genome by horizontal transfer 1-6 million years ago. Considering its transcriptional activity and the recent insertion of the element into at least two genes, Rider is a recently active retrotransposon in the tomato genome.

Ancient divergence of animal protein tyrosine kinase genes demonstrated by a gene family tree including choanoflagellate genes

Animal-specific gene families involved in cell–cell communication and developmental control comprise many subfamilies with distinct domain structures and functions. They diverged by subfamily-generating duplications and domain shufflings before the parazoan–eumetazoan split. Here, we have cloned 40 PTK cDNAs from choanoflagellates, Monosiga ovata, Stephanoeca diplocostata and Codosiga gracilis, the closest relatives to animals. A phylogeny-based analysis of PTKs revealed that 40 out of 47 subfamilies analyzed have unique domain structures and are possibly generated independently in animal and choanoflagellate lineages by domain shufflings. Seven cytoplasmic subfamilies showed divergence before the animal–choanoflagellate split originated by both duplications and shufflings.