Phylogenetic assessment of introns and SINEs within the Y chromosome using the cat family felidae as a species tree

Extensive Phylogenomic Discordance and the Complex Evolutionary History of the Neotropical Cat Genus Leopardus

Even in the genomics era, the phylogeny of Neotropical small felids comprised in the genus Leopardus remains contentious. We used whole-genome resequencing data to construct a time-calibrated consensus phylogeny of this group, quantify phylogenomic discordance, test for interspecies introgression, and assess patterns of genetic diversity and demographic history. We infer that the Leopardus radiation started in the Early Pliocene as an initial speciation burst, followed by another in its subgenus Oncifelis during the Early Pleistocene. Our findings challenge the long-held notion that ocelot (Leopardus pardalis) and margay (L. wiedii) are sister species and instead indicate that margay is most closely related to the enigmatic Andean cat (L. jacobita), whose whole-genome data are reported here for the first time. In addition, we found that the newly sampled Andean tiger cat (L. tigrinus pardinoides) population from Colombia associates closely with Central American tiger cats (L. tigrinus oncilla). Genealogical discordance was largely attributable to incomplete lineage sorting, yet was augmented by strong gene flow between ocelot and the ancestral branch of Oncifelis, as well as between Geoffroy's cat (L. geoffroyi) and southern tiger cat (L. guttulus). Contrasting demographic trajectories have led to disparate levels of current genomic diversity, with a nearly tenfold difference in heterozygosity between Andean cat and ocelot, spanning the entire range of variability found in extant felids. Our analyses improved our understanding of the speciation history and diversity patterns in this felid radiation, and highlight the benefits to phylogenomic inference of embracing the many heterogeneous signals scattered across the genome.

The roar of Rancho La Brea? Comparative anatomy of modern and fossil felid hyoid bones

Animal vocalization is broadly recognized as ecologically and evolutionarily important. In mammals, hyoid elements may influence vocalization repertoires because the hyoid apparatus anchors vocal tissues, and its morphology can be associated with variation in surrounding soft-tissue vocal anatomy. Thus, fossil hyoid morphology has the potential to shed light on vocalizations in extinct taxa. Yet, we know little about the hyoid morphology of extinct species because hyoid elements are rare in the fossil record. An exception is found in the Rancho La Brea tar pits in Los Angeles, California, where enough hyoids have been preserved to allow for quantitative analyses. The La Brea Tar Pits and Museum houses one of the largest and most diverse collections of carnivore fossils, including hyoid elements from the extinct felids Smilodon fatalis and Panthera atrox. Here, we found that extant members of Felinae (purring cats) and Panthera (roaring cats) showed characteristic differences in hyoid size and shape that suggest possible functional relationships with vocalization. The two extinct taxa had larger and more robust hyoids than extant felids, potentially reflecting the ability to produce lower frequency vocalizations as well as more substantial muscles associated with swallowing and respiration. Based on the shape of the hyoid elements, Panthera atrox resembled roaring cats, while Smilodon fatalis was quite variable and, contrary to suggestions from previous research, more similar overall to purring felids. Thus P. atrox may have roared and S. fatalis may have produced vocalizations similar to extant purring cats but at a lower frequency. Due to the confounding of vocalization repertoire and phylogenetic history in extant Felidae, we cannot distinguish between morphological signals related to vocalization behavior and those related to shared evolutionary history unrelated to vocalization.

Homoplasy of Retrotransposon Insertions in Toothed Whales

Retrotransposon insertion patterns facilitate a virtually homoplasy-free picture of phylogenetic history. Still, a few most likely random parallel insertions or deletions result in rare cases of homoplasy in primates. The following question arises: how frequent is retrotransposon homoplasy in other phylogenetic clades? Here, we derived genome insertion data of toothed whales to evaluate the extension of homoplasy in a representative laurasiatherian group. Among more than a thousand extracted and aligned retrotransposon loci, we detected 37 cases of precise parallel insertions in species that are separated by over more than 10 million years, a time frame which minimizes the effects of incomplete lineage sorting. We compared the phylogenetic signal of insertions with the flanking sequences of these loci to further exclude potential polymorphic loci derived by incomplete lineage sorting. We found that the phylogenetic signals of retrotransposon insertion patterns exhibiting true homoplasy differ from the signals of their flanking sequences. In toothed whales, precise parallel insertions account for around 0.18-0.29% of insertion cases, which is about 12.5 times the frequency of such insertions among Alus in primates. We also detected five specific deletions of retrotransposons on various lineages of toothed whale evolution, a frequency of 0.003%, which is slightly higher than such occurrences in primates. Overall, the level of retrotransposon homoplasy in toothed whales is still marginal compared to the phylogenetic diagnostic retrotransposon presence/absence signal.

Taxonomic revision of the pampas cat Leopardus colocola complex (Carnivora: Felidae): an integrative approach

The pampas cat Leopardus colocola has been subject to conflicting classifications over the years. Currently, one polytypic species with seven subspecies is recognized, but integrative taxonomic study for this debated group has never been done. Here, we combine the broadest morphological coverage of the pampas cat to date with molecular data and ecological niche models to clarify its species composition and test the validity of recently proposed subspecies. The multiple lines of evidence derived from morphology, molecular, biogeography and climatic niche datasets converged on the recognition of five monotypic species: L. braccatus, L. colocola, L. garleppi (including thomasi, budini, steinbachi, crespoi and wolffsohni as synonyms), L. munoai and L. pajeros (including crucina as synonym). These five species are morphologically diagnosable based on skin and skull traits, have evolved in distinct climatic niche spaces and were recovered in molecular species delimitation. Contrary to previous taxonomic arrangements, we do not recognize subspecies in pampas cats. To objectively define the two most controversial species, we designate neotypes for L. colocola and L. pajeros. The diversification of pampas cats is associated with Middle Pleistocene glaciations, but additional genetic samples from the central Andean region are still needed to conclusively reconstruct its evolutionary history.

Lynx canadensis (Carnivora: Felidae)

 Lynx canadensis Kerr, 1792, commonly called the Canada lynx, is a medium size felid and is the second largest of the four species in the genus Lynx. It is distributed throughout the boreal forest of most of Canada and Alaska and across portions of the northern United States. It prefers dense, regenerating coniferous forests with moderate canopy and understory cover. L. canadensis is a snowshoe hare specialist, and its ecology, morphology, and behavior closely reflect that of its main prey. It is listed as “Least Concern” by the International Union for Conservation of Nature and Natural Resources, is on Appendix II of the Convention on International Trade in Endangered Species of Wild Fauna and Flora, and its population size trend is considered stable. However, the status of United States subpopulations, being largely peripheral to the Canadian population, is more tenuous and the species is protected.

True Homoplasy of Retrotransposon Insertions in Primates

How reliable are the presence/absence insertion patterns of the supposedly homoplasy-free retrotransposons, which were randomly inserted in the quasi infinite genomic space? To systematically examine this question in an up-to-date, multigenome comparison, we screened millions of primate transposed Alu SINE elements for incidences of homoplasious precise insertions and deletions. In genome-wide analyses, we identified and manually verified nine cases of precise parallel Alu insertions of apparently identical elements at orthologous positions in two ape lineages and twelve incidences of precise deletions of previously established SINEs. Correspondingly, eight precise parallel insertions and no exact deletions were detected in a comparison of lemuriform primate and human insertions spanning the range of primate diversity. With an overall frequency of homoplasious Alu insertions of only 0.01% (for human-chimpanzee-rhesus macaque) and 0.02-0.04% (for human-bushbaby-lemurs) and precise Alu deletions of 0.001-0.002% (for human-chimpanzee-rhesus macaque), real homoplasy is not considered to be a quantitatively relevant source of evolutionary noise. Thus, presence/absence patterns of Alu retrotransposons and, presumably, all LINE1-mobilized elements represent indeed the virtually homoplasy-free markers they are considered to be. Therefore, ancestral incomplete lineage sorting and hybridization remain the only serious sources of conflicting presence/absence patterns of retrotransposon insertions, and as such are detectable and quantifiable. [Homoplasy; precise deletions; precise parallel insertions; primates; retrotransposons.].

Exploring Massive Incomplete Lineage Sorting in Arctoids (Laurasiatheria, Carnivora)

Freed from the competition of large raptors, Paleocene carnivores could expand their newly acquired habitats in search of prey. Such changing conditions might have led to their successful distribution and rapid radiation. Today, molecular evolutionary biologists are faced, however, with the consequences of such accelerated adaptive radiations, because they led to sequential speciation more rapidly than phylogenetic markers could be fixed. The repercussions being that current genealogies based on such markers are incongruent with species trees.Our aim was to explore such conflicting phylogenetic zones of evolution during the early arctoid radiation, especially to distinguish diagnostic from misleading phylogenetic signals, and to examine other carnivore-related speciation events. We applied a combination of high-throughput computational strategies to screen carnivore and related genomes in silico for randomly inserted retroposed elements that we then used to identify inconsistent phylogenetic patterns in the Arctoidea group, which is well known for phylogenetic discordances.Our combined retrophylogenomic and in vitro wet lab approach detected hundreds of carnivore-specific insertions, many of them confirming well-established splits or identifying and solving conflicting species distributions. Our systematic genome-wide screens for Long INterspersed Elements detected homoplasy-free markers with insertion-specific truncation points that we used to distinguish phylogenetically informative markers from conflicting signals. The results were independently confirmed by phylogenetic diagnostic Short INterspersed Elements. As statistical analysis ruled out ancestral hybridization, these doubly verified but still conflicting patterns were statistically determined to be genomic remnants from a time of ancestral incomplete lineage sorting that especially accompanied large parts of Arctoidea evolution.

The dynamic proliferation of CanSINEs mirrors the complex evolution of Feliforms

BACKGROUND

Repetitive short interspersed elements (SINEs) are retrotransposons ubiquitous in mammalian genomes and are highly informative markers to identify species and phylogenetic associations. Of these, SINEs unique to the order Carnivora (CanSINEs) yield novel insights on genome evolution in domestic dogs and cats, but less is known about their role in related carnivores. In particular, genome-wide assessment of CanSINE evolution has yet to be completed across the Feliformia (cat-like) suborder of Carnivora. Within Feliformia, the cat family Felidae is composed of 37 species and numerous subspecies organized into eight monophyletic lineages that likely arose 10 million years ago. Using the Felidae family as a reference phylogeny, along with representative taxa from other families of Feliformia, the origin, proliferation and evolution of CanSINEs within the suborder were assessed.

RESULTS

We identified 93 novel intergenic CanSINE loci in Feliformia. Sequence analyses separated Feliform CanSINEs into two subfamilies, each characterized by distinct RNA polymerase binding motifs and phylogenetic associations. Subfamily I CanSINEs arose early within Feliformia but are no longer under active proliferation. Subfamily II loci are more recent, exclusive to Felidae and show evidence for adaptation to extant RNA polymerase activity. Further, presence/absence distributions of CanSINE loci are largely congruent with taxonomic expectations within Feliformia and the less resolved nodes in the Felidae reference phylogeny present equally ambiguous CanSINE data. SINEs are thought to be nearly impervious to excision from the genome. However, we observed a nearly complete excision of a CanSINEs locus in puma (Puma concolor). In addition, we found that CanSINE proliferation in Felidae frequently targeted existing CanSINE loci for insertion sites, resulting in tandem arrays.

CONCLUSIONS

We demonstrate the existence of at least two SINE families within the Feliformia suborder, one of which is actively involved in insertional mutagenesis. We find SINEs are powerful markers of speciation and conclude that the few inconsistencies with expected patterns of speciation likely represent incomplete lineage sorting, species hybridization and SINE-mediated genome rearrangement.

AST: an automated sequence-sampling method for improving the taxonomic diversity of gene phylogenetic trees

A challenge in phylogenetic inference of gene trees is how to properly sample a large pool of homologous sequences to derive a good representative subset of sequences. Such a need arises in various applications, e.g. when (1) accuracy-oriented phylogenetic reconstruction methods may not be able to deal with a large pool of sequences due to their high demand in computing resources; (2) applications analyzing a collection of gene trees may prefer to use trees with fewer operational taxonomic units (OTUs), for instance for the detection of horizontal gene transfer events by identifying phylogenetic conflicts; and (3) the pool of available sequences is biased towards extensively studied species. In the past, the creation of subsamples often relied on manual selection. Here we present an Automated sequence-Sampling method for improving the Taxonomic diversity of gene phylogenetic trees, AST, to obtain representative sequences that maximize the taxonomic diversity of the sampled sequences. To demonstrate the effectiveness of AST, we have tested it to solve four problems, namely, inference of the evolutionary histories of the small ribosomal subunit protein S5 of E. coli, 16 S ribosomal RNAs and glycosyl-transferase gene family 8, and a study of ancient horizontal gene transfers from bacteria to plants. Our results show that the resolution of our computational results is almost as good as that of manual inference by domain experts, hence making the tool generally useful to phylogenetic studies by non-phylogeny specialists. The program is available at http://csbl.bmb.uga.edu/~zhouchan/AST.php.

Development of SNP markers identifying European wildcats, domestic cats, and their admixed progeny

Introgression can be an important evolutionary force but it can also lead to species extinction and as such is a crucial issue for species conservation. However, introgression is difficult to detect, morphologically as well as genetically. Hybridization with domestic cats (Felis silvestris catus) is a major concern for the conservation of European wildcats (Felis s. silvestris). The available morphologic and genetic markers for the two Felis subspecies are not sufficient to reliably detect hybrids beyond first generation. Here we present a single nucleotide polymorphism (SNP) based approach that allows the identification of introgressed individuals. Using high-throughput sequencing of reduced representation libraries we developed a diagnostic marker set containing 48 SNPs (Fst > 0.8) which allows the identification of wildcats, domestic cats, their hybrids and backcrosses. This allows assessing introgression rate in natural wildcat populations and is key for a better understanding of hybridization processes.

[Summary of phylogeny in family Felidae of Carnivora]

Felidae (cats) is one of the strict carnivorous groups in the order Carnivora, many of which are most familiar and spectacular to us. They are the top predators in the world. Thirty-six of 37 living cat species are considered as either "endangered" or "threatened". The relationships among species of the family Felidae, which evolved recently and rapidly, are difficult to resolve, and have been the subject of debate. Construction of a reliable Felidae phylogeny will be of evolutionarily significance and conservation value. In this paper, we summarized phylogeny of Felidae, including cytological, morphological and molecular evidence, and pointed out the existing phylogenetic problems. This review is expected to guide future researches of Felidae phylogeny, and to lay a theoretic foundation for the protection of this animal group.

Leopardus braccatus (Carnivora: Felidae)

Leopardus braccatus (Cope, 1889) is a small felid—not much larger than a domestic house cat—commonly called the Pantanal cat. No comprehensive surveys have been conducted to determine how many L. braccatus exist in the wild. It is found in humid, warm grasslands and wooded areas of extreme northwestern Argentina, southwestern and south- and north-central (newly reported ranges) Brazil, Paraguay, and Uruguay. Habitat loss and degradation are considered major threats throughout most of its range. It is uncommon in captivity and museum collections, listed with all Felidae under Appendix II of the Convention on International Trade in Endangered Species of Wild Fauna and Flora, and considered “Near Threatened” by the International Union for Conservation of Nature and Natural Resources in the pampas cat group ( = L. colocolo).

Carnivore-specific SINEs (Can-SINEs): distribution, evolution, and genomic impact

Short interspersed nuclear elements (SINEs) are a type of class 1 transposable element (retrotransposon) with features that allow investigators to resolve evolutionary relationships between populations and species while providing insight into genome composition and function. Characterization of a Carnivora-specific SINE family, Can-SINEs, has, has aided comparative genomic studies by providing rare genomic changes, and neutral sequence variants often needed to resolve difficult evolutionary questions. In addition, Can-SINEs constitute a significant source of functional diversity with Carnivora. Publication of the whole-genome sequence of domestic dog, domestic cat, and giant panda serves as a valuable resource in comparative genomic inferences gleaned from Can-SINEs. In anticipation of forthcoming studies bolstered by new genomic data, this review describes the discovery and characterization of Can-SINE motifs as well as describes composition, distribution, and effect on genome function. As the contribution of noncoding sequences to genomic diversity becomes more apparent, SINEs and other transposable elements will play an increasingly large role in mammalian comparative genomics.

Mitogenomic analysis of the genus Panthera

The complete sequences of the mitochondrial DNA genomes of Panthera tigris, Panthera pardus, and Panthera uncia were determined using the polymerase chain reaction method. The lengths of the complete mitochondrial DNA sequences of the three species were 16990, 16964, and 16773 bp, respectively. Each of the three mitochondrial DNA genomes included 13 protein-coding genes, 22 tRNA, two rRNA, one O(L)R, and one control region. The structures of the genomes were highly similar to those of Felis catus, Acinonyx jubatus, and Neofelis nebulosa. The phylogenies of the genus Panthera were inferred from two combined mitochondrial sequence data sets and the complete mitochondrial genome sequences, by MP (maximum parsimony), ML (maximum likelihood), and Bayesian analysis. The results showed that Panthera was composed of Panthera leo, P. uncia, P. pardus, Panthera onca, P. tigris, and N. nebulosa, which was included as the most basal member. The phylogeny within Panthera genus was N. nebulosa (P. tigris (P. onca (P. pardus, (P. leo, P. uncia)))). The divergence times for Panthera genus were estimated based on the ML branch lengths and four well-established calibration points. The results showed that at about 11.3 MYA, the Panthera genus separated from other felid species and then evolved into the several species of the genus. In detail, N. nebulosa was estimated to be founded about 8.66 MYA, P. tigris about 6.55 MYA, P. uncia about 4.63 MYA, and P. pardus about 4.35 MYA. All these estimated times were older than those estimated from the fossil records. The divergence event, evolutionary process, speciation, and distribution pattern of P. uncia, a species endemic to the central Asia with core habitats on the Qinghai-Tibetan Plateau and surrounding highlands, mostly correlated with the geological tectonic events and intensive climate shifts that happened at 8, 3.6, 2.5, and 1.7 MYA on the plateau during the late Cenozoic period.

On the phylogeny of Mustelidae subfamilies: analysis of seventeen nuclear non-coding loci and mitochondrial complete genomes

BACKGROUND

Mustelidae, as the largest and most-diverse family of order Carnivora, comprises eight subfamilies. Phylogenetic relationships among these Mustelidae subfamilies remain argumentative subjects in recent years. One of the main reasons is that the mustelids represent a typical example of rapid evolutionary radiation and recent speciation event. Prior investigation has been concentrated on the application of different mitochondrial (mt) sequence and nuclear protein-coding data, herein we employ 17 nuclear non-coding loci (>15 kb), in conjunction with mt complete genome data (>16 kb), to clarify these enigmatic problems.

RESULTS

The combined nuclear intron and mt genome analyses both robustly support that Taxidiinae diverged first, followed by Melinae. Lutrinae and Mustelinae are grouped together in all analyses with strong supports. The position of Helictidinae, however, is enigmatic because the mt genome analysis places it to the clade uniting Lutrinae and Mustelinae, whereas the nuclear intron analysis favors a novel view supporting a closer relationship of Helictidinae to Martinae. This finding emphasizes a need to add more data and include more taxa to resolve this problem. In addition, the molecular dating provides insights into the time scale of the origin and diversification of the Mustelidae subfamilies. Finally, the phylogenetic performances and limits of nuclear introns and mt genes are discussed in the context of Mustelidae phylogeny.

CONCLUSION

Our study not only brings new perspectives on the previously obscured phylogenetic relationships among Mustelidae subfamilies, but also provides another example demonstrating the effectiveness of nuclear non-coding loci for reconstructing evolutionary histories in a group that has undergone rapid bursts of speciation.

A molecular phylogeny of living primates

Comparative genomic analyses of primates offer considerable potential to define and understand the processes that mold, shape, and transform the human genome. However, primate taxonomy is both complex and controversial, with marginal unifying consensus of the evolutionary hierarchy of extant primate species. Here we provide new genomic sequence (∼8 Mb) from 186 primates representing 61 (∼90%) of the described genera, and we include outgroup species from Dermoptera, Scandentia, and Lagomorpha. The resultant phylogeny is exceptionally robust and illuminates events in primate evolution from ancient to recent, clarifying numerous taxonomic controversies and providing new data on human evolution. Ongoing speciation, reticulate evolution, ancient relic lineages, unequal rates of evolution, and disparate distributions of insertions/deletions among the reconstructed primate lineages are uncovered. Our resolution of the primate phylogeny provides an essential evolutionary framework with far-reaching applications including: human selection and adaptation, global emergence of zoonotic diseases, mammalian comparative genomics, primate taxonomy, and conservation of endangered species.

Advances in human biomedicine, including those focused on changes in genes triggered or disrupted in development, resistance/susceptibility to infectious disease, cancers, mechanisms of recombination, and genome plasticity, cannot be adequately interpreted in the absence of a precise evolutionary context or hierarchy. However, little is known about the genomes of other primate species, a situation exacerbated by a paucity of nuclear molecular sequence data necessary to resolve the complexities of primate divergence over time. We overcome this deficiency by sequencing 54 nuclear gene regions from DNA samples representing ∼90% of the diversity present in living primates. We conduct a phylogenetic analysis to determine the origin, evolution, patterns of speciation, and unique features in genome divergence among primate lineages. The resultant phylogenetic tree is remarkably robust and unambiguously resolves many long-standing issues in primate taxonomy. Our data provide a strong foundation for illuminating those genomic differences that are uniquely human and provide new insights on the breadth and richness of gene evolution across all primate lineages.

Against expectation: a short sequence with high signal elucidates cone snail phylogeny

A short (259 nucleotide) conserved intronic sequence (CIS) is surprisingly informative for delineating deep phylogenetic relationships in cone snails. Conus species previously have been assigned to clades based on the evidence from mitochondrial 12S and 16S rRNA gene sequences (1129 bp). Despite their length, these genes lack the phylogenetic information necessary to resolve the relationships among the clades. Here we show that the relationships can be inferred from just 46 sites in the very short CIS sequence (a portion of "intron 9" of the γ-glutamyl carboxylase gene). This is counterintuitive because in short sequences sampling error (noise) often drowns out phylogenetic signal. The intron 9 CIS is rich in synapomorphies that define the divergence patterns among eight clades of worm- and fish-hunting Conus, and it contains almost no homoplasy. Parsimony, maximum likelihood and Bayesian analyses of the combined sequences (mt rRNA+CIS) confirm most of the relationships among 23 Conus sequences. This phylogeny implies that fish-hunting behavior evolved at least twice during the history of Conus-once among New World species and independently in the Indo-Pacific clades.

Occurrence of Can-SINEs and intron sequence evolution supports robust phylogeny of pinniped carnivores and their terrestrial relatives

Investigating the dog genome we found 178965 introns with a moderate length of 200-1000 bp. A screening of these sequences against 23 different repeat libraries to find insertions of short interspersed elements (SINEs) detected 45276 SINEs. Virtually all of these SINEs (98%) belong to the tRNA-derived Can-SINE family. Can-SINEs arose about 55 million years ago before Carnivora split into two basal groups, the Caniformia (dog-like carnivores) and the Feliformia (cat-like carnivores). Genome comparisons of dog and cat recovered 506 putatively informative SINE loci for caniformian phylogeny. In this study we show how to use such genome information of model organisms to research the phylogeny of related non-model species of interest. Investigating a dataset including representatives of all major caniformian lineages, we analysed 24 randomly chosen loci for 22 taxa. All loci were amplifiable and revealed 17 parsimony-informative SINE insertions. The screening for informative SINE insertions yields a large amount of sequence information, in particular of introns, which contain reliable phylogenetic information as well. A phylogenetic analysis of intron- and SINE sequence data provided a statistically robust phylogeny which is congruent with the absence/presence pattern of our SINE markers. This phylogeny strongly supports a sistergroup relationship of Musteloidea and Pinnipedia. Within Pinnipedia, we see strong support from bootstrapping and the presence of a SINE insertion for a sistergroup relationship of the walrus with the Otariidae.

Molecular diversity and phylogeography of the Asian leopard cat, Felis bengalensis, inferred from mitochondrial and Y-chromosomal DNA sequences

To investigate genetic diversity and phylogeography of the Asian leopard cat (Felis bengalensis), mitochondrial DNA (mtDNA) sequences were determined for 39 individuals from various areas. Sequences combining the complete cytochrome b gene (1,140 bp) with the partial control region (646-810 bp) were classified into 24 haplotypes: 21 types from 21 animals, one from eight animals from Tsushima Islands, one from eight animals from Iriomote Island, and one from two animals from Southeast Asia. Phylogenetic trees of the 24 haplotypes clearly showed three clades: a Northern Lineage and Southern Lineages 1 and 2. The Northern Lineage consisted of animals from Tsushima Islands, the Korean Peninsula, the continental Far East, Taiwan, and Iriomote Island. Within the Northern Lineage, genetic contacts could have occurred between geographically neighboring populations before isolation by straits. Southern Lineage 1, comprising Southeast Asian animals, showed higher genetic diversity. Southern Lineage 2 had large genetic distances from other lineages. Within the control region, the Asian leopard cats shared two to four repetitive motifs, and the number of motifs and their constitution were highly variable among individuals. The motifs were polymorphic even within individuals and could be classified into 31 types. Finally, males of mtDNA Southern Lineage 1 had either of two types of the Y-chromosomal gene ZFY, whereas all males of Northern Lineage shared only one type. Our results indicate that the diversity of southern populations is higher and that genetic differentiation among northern local populations reflects past geographical isolation.

Confirming the phylogeny of mammals by use of large comparative sequence data sets

The ongoing generation of prodigious amounts of genomic sequence data from myriad vertebrates is providing unparalleled opportunities for establishing definitive phylogenetic relationships among species. The size and complexities of such comparative sequence data sets not only allow smaller and more difficult branches to be resolved but also present unique challenges, including large computational requirements and the negative consequences of systematic biases. To explore these issues and to clarify the phylogenetic relationships among mammals, we have analyzed a large data set of over 60 megabase pairs (Mb) of high-quality genomic sequence, which we generated from 41 mammals and 3 other vertebrates. All sequences are orthologous to a 1.9-Mb region of the human genome that encompasses the cystic fibrosis transmembrane conductance regulator gene (CFTR). To understand the characteristics and challenges associated with phylogenetic analyses of such a large data set, we partitioned the sequence data in several ways and utilized maximum likelihood, maximum parsimony, and Neighbor-Joining algorithms, implemented in parallel on Linux clusters. These studies yielded well-supported phylogenetic trees, largely confirming other recent molecular phylogenetic analyses. Our results provide support for rooting the placental mammal tree between Atlantogenata (Xenarthra and Afrotheria) and Boreoeutheria (Euarchontoglires and Laurasiatheria), illustrate the difficulty in resolving some branches even with large amounts of data (e.g., in the case of Laurasiatheria), and demonstrate the valuable role that very large comparative sequence data sets can play in refining our understanding of the evolutionary relationships of vertebrates.

Evolution of feline immunodeficiency virus in Felidae: implications for human health and wildlife ecology

Genetic analyses of feline immunodeficiency viruses provide significant insights on the worldwide distribution and evolutionary history of this emerging pathogen. Large-scale screening of over 3000 samples from all species of Felidae indicates that at least some individuals from most species possess antibodies that cross react to FIV. Phylogenetic analyses of genetic variation in the pol-RT gene demonstrate that FIV lineages are species-specific and suggest that there has been a prolonged period of viral-host co-evolution. The clinical effects of FIV specific to species other than domestic cat are controversial. Comparative genomic analyses of all full-length FIV genomes confirmed that FIV is host specific. Recently sequenced lion subtype E is marginally more similar to Pallas cat FIV though env is more similar to that of domestic cat FIV, indicating a possible recombination between two divergent strains in the wild. Here we review global patterns of FIV seroprevalence and endemnicity, assess genetic differences within and between species-specific FIV strains, and interpret these with patterns of felid speciation to propose an ancestral origin of FIV in Africa followed by interspecies transmission and global dissemination to Eurasia and the Americas. Continued comparative genomic analyses of full-length FIV from all seropositive animals, along with whole genome sequence of host species, will greatly advance our understanding of the role of recombination, selection and adaptation in retroviral emergence.

Unequal rates of Y chromosome gene divergence during speciation of the family Ursidae

Evolution of the bear family Ursidae is well investigated in terms of morphological, paleontological, and genetic features. However, several phylogenetic ambiguities occur within the subfamily Ursinae (the family Ursidae excluding the giant panda and spectacled bear), which may correlate with behavioral traits of female philopatry and male-biased dispersal which form the basis of the observed matriarchal population structure in these species. In the process of bear evolution, we investigate the premise that such behavioral traits may be reflected in patterns of variation among genes with different modes of inheritance: matrilineal mitochondrial DNA (mtDNA), patrilineal Y chromosome, biparentally inherited autosomes, and the X chromosome. In the present study, we sequenced 3 Y-linked genes (3,453 bp) and 4 X-linked genes (4,960 bp) and reanalyzed previously published sequences from autosome genes (2,347 bp) in ursid species to investigate differences in evolutionary rates associated with patterns of inheritance. The results describe topological incongruence between sex-linked genes and autosome genes and between nuclear DNA and mtDNA. In more ancestral branches within the bear phylogeny, Y-linked genes evolved faster than autosome and X-linked genes, consistent with expectations based on male-driven evolution. However, this pattern changes among branches leading to each species within the lineage of Ursinae whereby the evolutionary rates of Y-linked genes have fewer than expected substitutions. This inconsistency between more recent nodes of the bear phylogeny with more ancestral nodes may reflect the influences of sex-biased dispersal as well as molecular evolutionary characteristics of the Y chromosome, and stochastic events in species natural history, and phylogeography unique to ursine bears.

Phylogenetic relationships among cetaceans revealed by Y-chromosome sequences

The Y chromosome has recently come into the spotlight as a new and efficient genetic marker for tracing paternal lineages. We reconstructed cetacean phylogeny using a 1.7-kbp fragment of the non-recombining Y chromosome (NRY), including the SRY gene and a flanking non-coding region. The topology of the Y-chromosome tree is robust to various methods of analysis and exhibits high branch-support values, possibly due to the absence of recombination, small effective population size, and low homoplasy. The Y-chromosome tree indicates monophyly of each suborder, Mysticeti and Odontoceti, with high branch support values (BS> or =86%; PP> or =98%). In the Odontoceti clade, three superfamilies, Physeteroidea, Ziphioidea, and Delphinoidea, diverged soon after the split between Mysticeti and Odontoceti. Our analysis allows resolution of this rapid radiation and indicates that Physeteroidea is basal in the Odontoceti clade (BS, 99%; PP, 100%; MBS, 61%). The major split within the superfamily Delphinoidea is between the Delphinidae clade and the Monodontidae+ Phocoenidae clade. The phylogenetic relationships among delphinid species are ambiguous, probably because of the rapid radiation of this family. In the Mysticeti clade, the first major split is between Balaenidae and Balaenopteridae; within Balaenopteridae, a Balaenoptera acutorostrata+B. bonaerensis (minke whales) clade forms a sister clade with the other balaenopterid species. Megaptera novaeangliae is nested within Balaenoptera, making the latter paraphyletic. The low homoplasy exhibited by the Y-chromosome data presented here suggests that an extended data set incorporating longer sequences would provide better resolution of cetacean lower-level pylogeny.

Spliceosomal introns as tools for genomic and evolutionary analysis

Over the past 5 years, the availability of dozens of whole genomic sequences from a wide variety of eukaryotic lineages has revealed a very large amount of information about the dynamics of intron loss and gain through eukaryotic history, as well as the evolution of intron sequences. Implicit in these advances is a great deal of information about the structure and evolution of surrounding sequences. Here, we review the wealth of ways in which structures of spliceosomal introns as well as their conservation and change through evolution may be harnessed for evolutionary and genomic analysis. First, we discuss uses of intron length distributions and positions in sequence assembly and annotation, and for improving alignment of homologous regions. Second, we review uses of introns in evolutionary studies, including the utility of introns as indicators of rates of sequence evolution, for inferences about molecular evolution, as signatures of orthology and paralogy, and for estimating rates of nucleotide substitution. We conclude with a discussion of phylogenetic methods utilizing intron sequences and positions.

Initial sequence and comparative analysis of the cat genome

The genome sequence (1.9-fold coverage) of an inbred Abyssinian domestic cat was assembled, mapped, and annotated with a comparative approach that involved cross-reference to annotated genome assemblies of six mammals (human, chimpanzee, mouse, rat, dog, and cow). The results resolved chromosomal positions for 663,480 contigs, 20,285 putative feline gene orthologs, and 133,499 conserved sequence blocks (CSBs). Additional annotated features include repetitive elements, endogenous retroviral sequences, nuclear mitochondrial (numt) sequences, micro-RNAs, and evolutionary breakpoints that suggest historic balancing of translocation and inversion incidences in distinct mammalian lineages. Large numbers of single nucleotide polymorphisms (SNPs), deletion insertion polymorphisms (DIPs), and short tandem repeats (STRs), suitable for linkage or association studies were characterized in the context of long stretches of chromosome homozygosity. In spite of the light coverage capturing approximately 65% of euchromatin sequence from the cat genome, these comparative insights shed new light on the tempo and mode of gene/genome evolution in mammals, promise several research applications for the cat, and also illustrate that a comparative approach using more deeply covered mammals provides an informative, preliminary annotation of a light (1.9-fold) coverage mammal genome sequence.

Development of Y chromosome intraspecific polymorphic markers in the Felidae

Y chromosome haplotyping based on microsatellites and single nucleotide polymorphisms (SNPs) has proved to be a powerful tool for population genetic studies of humans. However, the promise of the approach is hampered in the majority of nonhuman mammals by the lack of Y-specific polymorphic markers. We were able to identify new male-specific polymorphisms in the domestic cat Felis catus and 6 additional Felidae species with a combination of molecular genetic and cytogenetic approaches including 1) identifying domestic cat male-specific microsatellites from markers generated from a male cat microsatellite-enriched genomic library, a flow-sorted Y cosmid library, or a Y-specific cat bacteria artificial chromosome (BAC) clone, (2) constructing microsatellite-enriched libraries from flow-sorted Y chromosomes isolated directly from focal wildcat species, and (3) screening Y chromosome conserved anchored tagged sequences primers in Felidae species. Forty-one male-specific microsatellites were identified, but only 6 were single-copy loci, consistent with the repetitive nature of the Y chromosome. Nucleotide diversity (pi) of Y-linked intron sequences (2.1 kbp) was in the range of 0 (tiger) to 9.95 x 10(-4) (marbled cat), and the number of SNPs ranged from none in the tiger to 7 in the Asian leopard cat. The Y haplotyping system described here, consisting of 4 introns (SMCY3, SMCY7, UTY11, and DBY7) and 1 polymorphic microsatellite (SMCY-STR), represents the first available markers for tracking intraspecific male lineage polymorphisms in Felidae species and promises to provide significant insights to evolutionary and population genetic studies of the species.

Evolution of the male-determining gene SRY within the cat family Felidae

In most placental mammals, SRY is a single-copy gene located on the Y chromosome and is the trigger for male sex determination during embryonic development. Here, we present comparative genomic analyses of SRY (705 bp) along with the adjacent noncoding 5' flank (997 bp) and 3' flank (948 bp) in 36 species of the cat family Felidae. Phylogenetic analyses indicate that the noncoding genomic flanks and SRY closely track species divergence. However, several inconsistencies are observed in SRY. Overall, the gene exhibits purifying selection to maintain function (omega = 0.815) yet SRY is under positive selection in two of the eight felid lineages. SRY has low numbers of nucleotide substitutions, yet most encode amino acid changes between species, and four different species have significantly altered SRY due to insertion/deletions. Moreover, fixation of nonsynonymous substitutions between sister taxa is not consistent and may occur rapidly, as in the case of domestic cat, or not at all over long periods of time, as observed within the Panthera lineage. The former resembles positive selection during speciation, and the latter purifying selection to maintain function. Thus, SRY evolution in cats likely reflects the different phylogeographic histories, selection pressures, and patterns of speciation in modern felids.

Using genomic data to unravel the root of the placental mammal phylogeny

The phylogeny of placental mammals is a critical framework for choosing future genome sequencing targets and for resolving the ancestral mammalian genome at the nucleotide level. Despite considerable recent progress defining superordinal relationships, several branches remain poorly resolved, including the root of the placental tree. Here we analyzed the genome sequence assemblies of human, armadillo, elephant, and opossum to identify informative coding indels that would serve as rare genomic changes to infer early events in placental mammal phylogeny. We also expanded our species sampling by including sequence data from >30 ongoing genome projects, followed by PCR and sequencing validation of each indel in additional taxa. Our data provide support for a sister-group relationship between Afrotheria and Xenarthra (the Atlantogenata hypothesis), which is in turn the sister-taxon to Boreoeutheria. We failed to recover any indels in support of a basal position for Xenarthra (Epitheria), which is suggested by morphology and a recent retroposon analysis, or a hypothesis with Afrotheria basal (Exafricoplacentalia), which is favored by phylogenetic analysis of large nuclear gene data sets. In addition, we identified two retroposon insertions that also support Atlantogenata and none for the alternative hypotheses. A revised molecular timescale based on these phylogenetic inferences suggests Afrotheria and Xenarthra diverged from other placental mammals approximately 103 (95-114) million years ago. We discuss the impacts of this topology on earlier phylogenetic reconstructions and repeat-based inferences of phylogeny.

Indel evolution of mammalian introns and the utility of non-coding nuclear markers in eutherian phylogenetics

Nuclear DNA intron sequences are increasingly used to investigate evolutionary relationships among closely related organisms. The phylogenetic usefulness of intron sequences at higher taxonomic levels has, however, not been firmly established and very few studies have used these markers to address evolutionary questions above the family level. In addition, the mechanisms driving intron evolution are not well understood. We compared DNA sequence data derived from three presumably independently segregating introns (THY, PRKC I and MGF) across 158 mammalian species. All currently recognized extant eutherian mammalian orders were included with the exception of Cingulata, Dermoptera and Scandentia. The total aligned length of the data was 6366 base pairs (bp); after the exclusion of autapomorphic insertions, 1511 bp were analyzed. In many instances the Bayesian and parsimony analyses were complementary and gave significant posterior probability and bootstrap support (>80) for the monophyly of Afrotheria, Euarchontoglires, Laurasiatheria and Boreoeutheria. Apart from finding congruent support when using these methods, the intron data also provided several indels longer than 3 bp that support, among others, the monophyly of Afrotheria, Paenungulata, Ferae and Boreoeutheria. A quantitative analysis of insertions and deletions suggested that there was a 75% bias towards deletions. The average insertion size in the mammalian data set was 16.49 bp +/- 57.70 while the average deletion was much smaller (4.47 bp +/- 14.17). The tendency towards large insertions and small deletions is highlighted by the observation that out of a total of 17 indels larger than 100 bp, 15 were insertions. The majority of indels (>60% of all events) were 1 or 2 bp changes. Although the average overall indel substitution rate of 0.00559 per site is comparable to that previously reported for rodents and primates, individual analyses among different evolutionary lineages provide evidence for differences in the formation rate of indels among the different mammalian groups.

Extensive Morphological Convergence and Rapid Radiation in the Evolutionary History of the Family Geoemydidae (Old World Pond Turtles) Revealed by SINE Insertion Analysis

The family Geoemydidae is one of three in the superfamily Testudinoidea and is the most diversified family of extant turtle species. The phylogenetic relationships in this family and among related families have been vigorously investigated from both morphological and molecular viewpoints. The evolutionary history of Geoemydidae, however, remains controversial. Therefore, to elucidate the phylogenetic relationships of Geoemydidae and related species, we applied the SINE insertion method to investigate 49 informative SINE loci in 28 species. We detected four major evolutionary lineages (Testudinidae, Batagur group, Siebenrockiella group, and Geoemyda group) in the clade Testuguria (a clade of Geoemydidae + Testudinidae). All five specimens of Testudinidae form a monophyletic clade. The Batagur group comprises five batagurines. The Siebenrockiella group has one species, Siebenrockiella crassicollis. The Geoemyda group comprises 15 geoemydines (including three former batagurines, Mauremys reevesii, Mauremys sinensis, and Heosemys annandalii). Among these four groups, the SINE insertion patterns were inconsistent at four loci, suggesting that an ancestral species of Testuguria radiated and rapidly diverged into the four lineages during the initial stage of its evolution. Furthermore, within the Geoemyda group we identified three evolutionary lineages, namely Mauremys, Cuora, and Heosemys. The Heosemys lineage comprises Heosemys, Sacalia, Notochelys, and Melanochelys species, and its monophyly is a novel assemblage in Geoemydidae. Our SINE phylogenetic tree demonstrates extensive convergent morphological evolution between the Batagur group and the three species of the Geoemyda group, M. reevesii, M. sinensis, and H. annandalii.

Analyses of carnivore microsatellites and their intimate association with tRNA-derived SINEs

Background

The popularity of microsatellites has greatly increased in the last decade on account of their many applications. However, little is currently understood about the factors that influence their genesis and distribution among and within species genomes. In this work, we analyzed carnivore microsatellite clones from GenBank to study their association with interspersed repeats and elucidate the role of the latter in microsatellite genesis and distribution.

Results

We constructed a comprehensive carnivore microsatellite database comprising 1236 clones from GenBank. Thirty-three species of 11 out of 12 carnivore families were represented, although two distantly related species, the domestic dog and cat, were clearly overrepresented. Of these clones, 330 contained tRNA^Lys-derived SINEs and 357 contained other interspersed repeats. Our rough estimates of tRNA SINE copies per haploid genome were much higher than published ones. Our results also revealed a distinct juxtaposition of AG and A-rich repeats and tRNA^Lys-derived SINEs suggesting their coevolution. Both microsatellites arose repeatedly in two regions of the insterspersed repeat. Moreover, microsatellites associated with tRNA^Lys-derived SINEs showed the highest complexity and less potential instability.

Conclusion

Our results suggest that tRNA^Lys-derived SINEs are a significant source for microsatellite generation in carnivores, especially for AG and A-rich repeat motifs. These observations indicate two modes of microsatellite generation: the expansion and variation of pre-existing tandem repeats and the conversion of sequences with high cryptic simplicity into a repeat array; mechanisms which are not specific to tRNA^Lys-derived SINEs. Microsatellite and interspersed repeat coevolution could also explain different distribution of repeat types among and within species genomes.

Finally, due to their higher complexity and lower potential informative content of microsatellites associated with tRNA^Lys-derived SINEs, we recommend avoiding their use as genetic markers.

The late Miocene radiation of modern Felidae: a genetic assessment

Modern felid species descend from relatively recent (<11 million years ago) divergence and speciation events that produced successful predatory carnivores worldwide but that have confounded taxonomic classifications. A highly resolved molecular phylogeny with divergence dates for all living cat species, derived from autosomal, X-linked, Y-linked, and mitochondrial gene segments (22,789 base pairs) and 16 fossil calibrations define eight principal lineages produced through at least 10 intercontinental migrations facilitated by sea-level fluctuations. A ghost lineage analysis indicates that available felid fossils underestimate (i.e., unrepresented basal branch length) first occurrence by an average of 76%, revealing a low representation of felid lineages in paleontological remains. The phylogenetic performance of distinct gene classes showed that Y-chromosome segments are appreciably more informative than mitochondrial DNA, X-linked, or autosomal genes in resolving the rapid Felidae species radiation.

A Nuclear DNA Phylogenetic Perspective on the Evolution of Echolocation and Historical Biogeography of Extant Bats (Chiroptera)

Bats (Order Chiroptera), the only mammals capable of powered flight and sophisticated laryngeal echolocation, represent one of the most species-rich and ubiquitous orders of mammals. However, phylogenetic relationships within this group are poorly resolved. A robust evolutionary tree of Chiroptera is essential for evaluating the phylogeny of echolocation within Chiroptera, as well as for understanding their biogeographical history. We generated 4 kb of sequence data from portions of four novel nuclear intron markers for multiple representatives of 17 of the 18 recognized extant bat families, as well as the putative bat family Miniopteridae. Three echolocation-call characters were examined by mapping them onto the combined topology: (1) high-duty cycle versus low-duty cycle, (2) high-intensity versus low-intensity call emission, and (3) oral versus nasal emission. Echolocation seems to be highly convergent, and the mapping of echolocation-call design onto our phylogeny does not appear to resolve the question of whether echolocation had a single or two origins. Fossil taxa may also provide insight into the evolution of bats; we therefore evaluate 195 morphological characters in light of our nuclear DNA phylogeny. All but 24 of the morphological characters were found to be homoplasious when mapped onto the supermatrix topology, while the remaining characters provided insufficient information to reconstruct the placement of the fossil bat taxa with respect to extant families. However, a morphological synapomorphy characterizing the Rhinolophoidea was identified and is suggestive of a separate origin of echolocation in this clade. Dispersal-Vicariance analysis together with a relaxed Bayesian clock were used to evaluate possible biogeographic scenarios that could account for the current distribution pattern of extant bat families. Africa was reconstructed as the center of origin of modern-day bat families.

Phylogenomics and the reconstruction of the tree of life

As more complete genomes are sequenced, phylogenetic analysis is entering a new era - that of phylogenomics. One branch of this expanding field aims to reconstruct the evolutionary history of organisms on the basis of the analysis of their genomes. Recent studies have demonstrated the power of this approach, which has the potential to provide answers to several fundamental evolutionary questions. However, challenges for the future have also been revealed. The very nature of the evolutionary history of organisms and the limitations of current phylogenetic reconstruction methods mean that part of the tree of life might prove difficult, if not impossible, to resolve with confidence.