Using genome-wide measures of coancestry to maintain diversity and fitness in endangered and domestic pig populations

Conservation and breeding programs aim at maintaining the most diversity, thereby avoiding deleterious effects of inbreeding while maintaining enough variation from which traits of interest can be selected. Theoretically, the most diversity is maintained using optimal contributions based on many markers to calculate coancestries, but this can decrease fitness by maintaining linked deleterious variants. The heterogeneous patterns of coancestry displayed in pigs make them an excellent model to test these predictions. We propose methods to measure coancestry and fitness from resequencing data and use them in population management. We analyzed the resequencing data of Sus cebifrons, a highly endangered porcine species from the Philippines, and genotype data from the Pietrain domestic breed. By analyzing the demographic history of Sus cebifrons, we inferred two past bottlenecks that resulted in some inbreeding load. In Pietrain, we analyzed signatures of selection possibly associated with commercial traits. We also simulated the management of each population to assess the performance of different optimal contribution methods to maintain diversity, fitness, and selection signatures. Maximum genetic diversity was maintained using marker-by-marker coancestry, and least using genealogical coancestry. Using a measure of coancestry based on shared segments of the genome achieved the best results in terms of diversity and fitness. However, this segment-based management eliminated signatures of selection. We demonstrate that maintaining both diversity and fitness depends on the genomic distribution of deleterious variants, which is shaped by demographic and selection histories. Our findings show the importance of genomic and next-generation sequencing information in the optimal design of breeding or conservation programs.

Inference of gorilla demographic and selective history from whole-genome sequence data

Although population-level genomic sequence data have been gathered extensively for humans, similar data from our closest living relatives are just beginning to emerge. Examination of genomic variation within great apes offers many opportunities to increase our understanding of the forces that have differentially shaped the evolutionary history of hominid taxa. Here, we expand upon the work of the Great Ape Genome Project by analyzing medium to high coverage whole-genome sequences from 14 western lowland gorillas (Gorilla gorilla gorilla), 2 eastern lowland gorillas (G. beringei graueri), and a single Cross River individual (G. gorilla diehli). We infer that the ancestors of western and eastern lowland gorillas diverged from a common ancestor approximately 261 ka, and that the ancestors of the Cross River population diverged from the western lowland gorilla lineage approximately 68 ka. Using a diffusion approximation approach to model the genome-wide site frequency spectrum, we infer a history of western lowland gorillas that includes an ancestral population expansion of 1.4-fold around 970 ka and a recent 5.6-fold contraction in population size 23 ka. The latter may correspond to a major reduction in African equatorial forests around the Last Glacial Maximum. We also analyze patterns of variation among western lowland gorillas to identify several genomic regions with strong signatures of recent selective sweeps. We find that processes related to taste, pancreatic and saliva secretion, sodium ion transmembrane transport, and cardiac muscle function are overrepresented in genomic regions predicted to have experienced recent positive selection.

Comparative genomics reveals insights into avian genome evolution and adaptation

Birds are the most species-rich class of tetrapod vertebrates and have wide relevance across many research fields. We explored bird macroevolution using full genomes from 48 avian species representing all major extant clades. The avian genome is principally characterized by its constrained size, which predominantly arose because of lineage-specific erosion of repetitive elements, large segmental deletions, and gene loss. Avian genomes furthermore show a remarkably high degree of evolutionary stasis at the levels of nucleotide sequence, gene synteny, and chromosomal structure. Despite this pattern of conservation, we detected many non-neutral evolutionary changes in protein-coding genes and noncoding regions. These analyses reveal that pan-avian genomic diversity covaries with adaptations to different lifestyles and convergent evolution of traits.

Whole-genome analyses resolve early branches in the tree of life of modern birds

To better determine the history of modern birds, we performed a genome-scale phylogenetic analysis of 48 species representing all orders of Neoaves using phylogenomic methods created to handle genome-scale data. We recovered a highly resolved tree that confirms previously controversial sister or close relationships. We identified the first divergence in Neoaves, two groups we named Passerea and Columbea, representing independent lineages of diverse and convergently evolved land and water bird species. Among Passerea, we infer the common ancestor of core landbirds to have been an apex predator and confirm independent gains of vocal learning. Among Columbea, we identify pigeons and flamingoes as belonging to sister clades. Even with whole genomes, some of the earliest branches in Neoaves proved challenging to resolve, which was best explained by massive protein-coding sequence convergence and high levels of incomplete lineage sorting that occurred during a rapid radiation after the Cretaceous-Paleogene mass extinction event about 66 million years ago.

Molecular phylogeny and chromosomal evolution of Alcelaphini (Antilopinae)

Robertsonian (Rb) translocations, in particular centric fusions, are thought to play a primary role in evolution and speciation of the Bovidae family. However, Rb fusions are often polymorphic within species, being suggested as phylogenetically uninformative characters. This work studies chromosome variation in 72 captive individuals of 6 species of Alcelaphini (Antilopinae): The hartebeest (genus Alcelaphus), hirola (Beatragus), black and blue wildebeests (Connochaetes), and the topi and bontebok (Damaliscus). We infer the phylogenic relationships among Alcelaphini species and determine patterns of chromosomal evolution using G-banded karyotypes and complete mitochondrial genome sequences. The molecular phylogeny showed an early divergence of Connochaetes, followed by the split of Alcelaphus plus Beatragus + Damaliscus as sister taxa. Mitochondrial and chromosomal phylogenies only differed in the position of the critically endangered Beatragus, likely due to homoplasic chromosome characters. Patterns of chromosome evolution, reconstructed using a probabilistic approach, suggest that chromosome changes leading to speciation in Alcelaphini do not exclusively involve consecutive reduction of diploid number through centric fusion but also the losses and reversions of Rb translocations in Beatragus and Damaliscus lineages. Our results provide evidence that complex scenarios of chromosomal rearrangements can be detected in relatively recent-diverged bovids, as in this group of antelopes.

Characterization of Prdm9 in equids and sterility in mules

Prdm9 (Meisetz) is the first speciation gene discovered in vertebrates conferring reproductive isolation. This locus encodes a meiosis-specific histone H3 methyltransferase that specifies meiotic recombination hotspots during gametogenesis. Allelic differences in Prdm9, characterized for a variable number of zinc finger (ZF) domains, have been associated with hybrid sterility in male house mice via spermatogenic failure at the pachytene stage. The mule, a classic example of hybrid sterility in mammals also exhibits a similar spermatogenesis breakdown, making Prdm9 an interesting candidate to evaluate in equine hybrids. In this study, we characterized the Prdm9 gene in all species of equids by analyzing sequence variation of the ZF domains and estimating positive selection. We also evaluated the role of Prdm9 in hybrid sterility by assessing allelic differences of ZF domains in equine hybrids. We found remarkable variation in the sequence and number of ZF domains among equid species, ranging from five domains in the Tibetan kiang and Asiatic wild ass, to 14 in the Grevy’s zebra. Positive selection was detected in all species at amino acid sites known to be associated with DNA-binding specificity of ZF domains in mice and humans. Equine hybrids, in particular a quartet pedigree composed of a fertile mule showed a mosaic of sequences and number of ZF domains suggesting that Prdm9 variation does not seem by itself to contribute to equine hybrid sterility.

Diverse captive non-human primates with phytanic acid-deficient diets rich in plant products have substantial phytanic acid levels in their red blood cells

BACKGROUND

Humans and rodents with impaired phytanic acid (PA) metabolism can accumulate toxic stores of PA that have deleterious effects on multiple organ systems. Ruminants and certain fish obtain PA from the microbial degradation of dietary chlorophyll and/or through chlorophyll-derived precursors. In contrast, humans cannot derive PA from chlorophyll and instead normally obtain it only from meat, dairy, and fish products.

RESULTS

Captive apes and Old world monkeys had significantly higher red blood cell (RBC) PA levels relative to humans when all subjects were fed PA-deficient diets. Given the adverse health effects resulting from PA over accumulation, we investigated the molecular evolution of thirteen PA metabolism genes in apes, Old world monkeys, and New world monkeys. All non-human primate (NHP) orthologs are predicted to encode full-length proteins with the marmoset Phyh gene containing a rare, but functional, GA splice donor dinucleotide. Acox2, Scp2, and Pecr sequences had amino acid positions with accelerated substitution rates while Amacr had significant variation in evolutionary rates in apes relative to other primates.

CONCLUSIONS

Unlike humans, diverse captive NHPs with PA-deficient diets rich in plant products have substantial RBC PA levels. The favored hypothesis is that NHPs can derive significant amounts of PA from the degradation of ingested chlorophyll through gut fermentation. If correct, this raises the possibility that RBC PA levels could serve as a biomarker for evaluating the digestive health of captive NHPs. Furthermore, the evolutionary rates of the several genes relevant to PA metabolism provide candidate genetic adaptations to NHP diets.

Conservation genomics of threatened animal species

The genomics era has opened up exciting possibilities in the field of conservation biology by enabling genomic analyses of threatened species that previously were limited to model organisms. Next-generation sequencing (NGS) and the collection of genome-wide data allow for more robust studies of the demographic history of populations and adaptive variation associated with fitness and local adaptation. Genomic analyses can also advance management efforts for threatened wild and captive populations by identifying loci contributing to inbreeding depression and disease susceptibility, and predicting fitness consequences of introgression. However, the development of genomic tools in wild species still carries multiple challenges, particularly those associated with computational and sampling constraints. This review provides an overview of the most significant applications of NGS and the implications and limitations of genomic studies in conservation.

Macroevolutionary dynamics and historical biogeography of primate diversification inferred from a species supermatrix

Phylogenetic relationships, divergence times, and patterns of biogeographic descent among primate species are both complex and contentious. Here, we generate a robust molecular phylogeny for 70 primate genera and 367 primate species based on a concatenation of 69 nuclear gene segments and ten mitochondrial gene sequences, most of which were extracted from GenBank. Relaxed clock analyses of divergence times with 14 fossil-calibrated nodes suggest that living Primates last shared a common ancestor 71-63 Ma, and that divergences within both Strepsirrhini and Haplorhini are entirely post-Cretaceous. These results are consistent with the hypothesis that the Cretaceous-Paleogene mass extinction of non-avian dinosaurs played an important role in the diversification of placental mammals. Previous queries into primate historical biogeography have suggested Africa, Asia, Europe, or North America as the ancestral area of crown primates, but were based on methods that were coopted from phylogeny reconstruction. By contrast, we analyzed our molecular phylogeny with two methods that were developed explicitly for ancestral area reconstruction, and find support for the hypothesis that the most recent common ancestor of living Primates resided in Asia. Analyses of primate macroevolutionary dynamics provide support for a diversification rate increase in the late Miocene, possibly in response to elevated global mean temperatures, and are consistent with the fossil record. By contrast, diversification analyses failed to detect evidence for rate-shift changes near the Eocene-Oligocene boundary even though the fossil record provides clear evidence for a major turnover event ("Grande Coupure") at this time. Our results highlight the power and limitations of inferring diversification dynamics from molecular phylogenies, as well as the sensitivity of diversification analyses to different species concepts.

Incorporating indels as phylogenetic characters: impact for interfamilial relationships within Arctoidea (Mammalia: Carnivora)

Insertion and deletion events (indels) provide a suite of markers with enormous potential for molecular phylogenetics. Using many more indel characters than those in previous studies, we here for the first time address the impact of indel inclusion on the phylogenetic inferences of Arctoidea (Mammalia: Carnivora). Based on 6843 indel characters from 22 nuclear intron loci of 16 species of Arctoidea, our analyses demonstrate that when the indels were not taken into consideration, the monophyly of Ursidae and Pinnipedia tree and the monophyly of Pinnipedia and Musteloidea tree were both recovered, whereas inclusion of indels by using three different indel coding schemes give identical phylogenetic tree topologies supporting the monophyly of Ursidae and Pinnipedia. Our work brings new perspectives on the previously controversial placements among Arctoidea families, and provides another example demonstrating the importance of identifying and incorporating indels in the phylogenetic analyses of introns. In addition, comparison of indel incorporation methods revealed that the three indel coding methods are all advantageous over treating indels as missing data, given that incorporating indels produces consistent results across methods. This is the first report of the impact of different indel coding schemes on phylogenetic reconstruction at the family level in Carnivora, which indicates that indels should be taken into account in the future phylogenetic analyses.

Tissue sampling methods and standards for vertebrate genomics

The recent rise in speed and efficiency of new sequencing technologies have facilitated high-throughput sequencing, assembly and analyses of genomes, advancing ongoing efforts to analyze genetic sequences across major vertebrate groups. Standardized procedures in acquiring high quality DNA and RNA and establishing cell lines from target species will facilitate these initiatives. We provide a legal and methodological guide according to four standards of acquiring and storing tissue for the Genome 10K Project and similar initiatives as follows: four-star (banked tissue/cell cultures, RNA from multiple types of tissue for transcriptomes, and sufficient flash-frozen tissue for 1 mg of DNA, all from a single individual); three-star (RNA as above and frozen tissue for 1 mg of DNA); two-star (frozen tissue for at least 700 μg of DNA); and one-star (ethanol-preserved tissue for 700 μg of DNA or less of mixed quality). At a minimum, all tissues collected for the Genome 10K and other genomic projects should consider each species' natural history and follow institutional and legal requirements. Associated documentation should detail as much information as possible about provenance to ensure representative sampling and subsequent sequencing. Hopefully, the procedures outlined here will not only encourage success in the Genome 10K Project but also inspire the adaptation of standards by other genomic projects, including those involving other biota.

The fishes of Genome 10K

The Genome 10K project aims to sequence the genomes of 10,000 vertebrates, representing approximately one genome for each vertebrate genus. Since fishes (cartilaginous fishes, ray-finned fishes and lobe-finned fishes) represent more than 50% of extant vertebrates, it is planned to target 4,000 fish genomes. At present, nearly 60 fish genomes are being sequenced at various public funded labs, and under a Genome 10K and BGI pilot project. An additional 100 fishes have been identified for sequencing in the next phase of Genome 10K project.

Insights into hominid evolution from the gorilla genome sequence

Gorillas are humans’ closest living relatives after chimpanzees, and are of comparable importance for the study of human origins and evolution. Here we present the assembly and analysis of a genome sequence for the western lowland gorilla, and compare the whole genomes of all extant great ape genera. We propose a synthesis of genetic and fossil evidence consistent with placing the human-chimpanzee and human-chimpanzee-gorilla speciation events at approximately 6 and 10 million years ago (Mya). In 30% of the genome, gorilla is closer to human or chimpanzee than the latter are to each other; this is rarer around coding genes, indicating pervasive selection throughout great ape evolution, and has functional consequences in gene expression. A comparison of protein coding genes reveals approximately 500 genes showing accelerated evolution on each of the gorilla, human and chimpanzee lineages, and evidence for parallel acceleration, particularly of genes involved in hearing. We also compare the western and eastern gorilla species, estimating an average sequence divergence time 1.75 million years ago, but with evidence for more recent genetic exchange and a population bottleneck in the eastern species. The use of the genome sequence in these and future analyses will promote a deeper understanding of great ape biology and evolution.

A high density SNP array for the domestic horse and extant Perissodactyla: utility for association mapping, genetic diversity, and phylogeny studies

An equine SNP genotyping array was developed and evaluated on a panel of samples representing 14 domestic horse breeds and 18 evolutionarily related species. More than 54,000 polymorphic SNPs provided an average inter-SNP spacing of ∼43 kb. The mean minor allele frequency across domestic horse breeds was 0.23, and the number of polymorphic SNPs within breeds ranged from 43,287 to 52,085. Genome-wide linkage disequilibrium (LD) in most breeds declined rapidly over the first 50-100 kb and reached background levels within 1-2 Mb. The extent of LD and the level of inbreeding were highest in the Thoroughbred and lowest in the Mongolian and Quarter Horse. Multidimensional scaling (MDS) analyses demonstrated the tight grouping of individuals within most breeds, close proximity of related breeds, and less tight grouping in admixed breeds. The close relationship between the Przewalski's Horse and the domestic horse was demonstrated by pair-wise genetic distance and MDS. Genotyping of other Perissodactyla (zebras, asses, tapirs, and rhinoceros) was variably successful, with call rates and the number of polymorphic loci varying across taxa. Parsimony analysis placed the modern horse as sister taxa to Equus przewalski. The utility of the SNP array in genome-wide association was confirmed by mapping the known recessive chestnut coat color locus (MC1R) and defining a conserved haplotype of ∼750 kb across all breeds. These results demonstrate the high quality of this SNP genotyping resource, its usefulness in diverse genome analyses of the horse, and potential use in related species.

Comparative biology of mammalian telomeres: hypotheses on ancestral states and the roles of telomeres in longevity determination

Progressive telomere shortening from cell division (replicative aging) provides a barrier for human tumor progression. This program is not conserved in laboratory mice, which have longer telomeres and constitutive telomerase. Wild species that do/do not use replicative aging have been reported, but the evolution of different phenotypes and a conceptual framework for understanding their uses of telomeres is lacking. We examined telomeres/telomerase in cultured cells from > 60 mammalian species to place different uses of telomeres in a broad mammalian context. Phylogeny-based statistical analysis reconstructed ancestral states. Our analysis suggested that the ancestral mammalian phenotype included short telomeres (< 20 kb, as we now see in humans) and repressed telomerase. We argue that the repressed telomerase was a response to a higher mutation load brought on by the evolution of homeothermy. With telomerase repressed, we then see the evolution of replicative aging. Telomere length inversely correlated with lifespan, while telomerase expression co-evolved with body size. Multiple independent times smaller, shorter-lived species changed to having longer telomeres and expressing telomerase. Trade-offs involving reducing the energetic/cellular costs of specific oxidative protection mechanisms (needed to protect < 20 kb telomeres in the absence of telomerase) could explain this abandonment of replicative aging. These observations provide a conceptual framework for understanding different uses of telomeres in mammals, support a role for human-like telomeres in allowing longer lifespans to evolve, demonstrate the need to include telomere length in the analysis of comparative studies of oxidative protection in the biology of aging, and identify which mammals can be used as appropriate model organisms for the study of the role of telomeres in human cancer and aging.

Comparative biology of mammalian telomeres: hypotheses on ancestral states and the roles of telomeres in longevity determination

Progressive telomere shortening from cell division (replicative aging) provides a barrier for human tumor progression. This program is not conserved in laboratory mice, which have longer telomeres and constitutive telomerase. Wild species that do/do not use replicative aging have been reported, but the evolution of different phenotypes and a conceptual framework for understanding their uses of telomeres is lacking. We examined telomeres/telomerase in cultured cells from > 60 mammalian species to place different uses of telomeres in a broad mammalian context. Phylogeny-based statistical analysis reconstructed ancestral states. Our analysis suggested that the ancestral mammalian phenotype included short telomeres (< 20 kb, as we now see in humans) and repressed telomerase. We argue that the repressed telomerase was a response to a higher mutation load brought on by the evolution of homeothermy. With telomerase repressed, we then see the evolution of replicative aging. Telomere length inversely correlated with lifespan, while telomerase expression co-evolved with body size. Multiple independent times smaller, shorter-lived species changed to having longer telomeres and expressing telomerase. Trade-offs involving reducing the energetic/cellular costs of specific oxidative protection mechanisms (needed to protect < 20 kb telomeres in the absence of telomerase) could explain this abandonment of replicative aging. These observations provide a conceptual framework for understanding different uses of telomeres in mammals, support a role for human-like telomeres in allowing longer lifespans to evolve, demonstrate the need to include telomere length in the analysis of comparative studies of oxidative protection in the biology of aging, and identify which mammals can be used as appropriate model organisms for the study of the role of telomeres in human cancer and aging.

Impacts of the Cretaceous Terrestrial Revolution and KPg extinction on mammal diversification

Previous analyses of relations, divergence times, and diversification patterns among extant mammalian families have relied on supertree methods and local molecular clocks. We constructed a molecular supermatrix for mammalian families and analyzed these data with likelihood-based methods and relaxed molecular clocks. Phylogenetic analyses resulted in a robust phylogeny with better resolution than phylogenies from supertree methods. Relaxed clock analyses support the long-fuse model of diversification and highlight the importance of including multiple fossil calibrations that are spread across the tree. Molecular time trees and diversification analyses suggest important roles for the Cretaceous Terrestrial Revolution and Cretaceous-Paleogene (KPg) mass extinction in opening up ecospace that promoted interordinal and intraordinal diversification, respectively. By contrast, diversification analyses provide no support for the hypothesis concerning the delayed rise of present-day mammals during the Eocene Period.

A massively parallel sequencing approach uncovers ancient origins and high genetic variability of endangered Przewalski's horses

The endangered Przewalski's horse is the closest relative of the domestic horse and is the only true wild horse species surviving today. The question of whether Przewalski's horse is the direct progenitor of domestic horse has been hotly debated. Studies of DNA diversity within Przewalski's horses have been sparse but are urgently needed to ensure their successful reintroduction to the wild. In an attempt to resolve the controversy surrounding the phylogenetic position and genetic diversity of Przewalski's horses, we used massively parallel sequencing technology to decipher the complete mitochondrial and partial nuclear genomes for all four surviving maternal lineages of Przewalski's horses. Unlike single-nucleotide polymorphism (SNP) typing usually affected by ascertainment bias, the present method is expected to be largely unbiased. Three mitochondrial haplotypes were discovered—two similar ones, haplotypes I/II, and one substantially divergent from the other two, haplotype III. Haplotypes I/II versus III did not cluster together on a phylogenetic tree, rejecting the monophyly of Przewalski's horse maternal lineages, and were estimated to split 0.117–0.186 Ma, significantly preceding horse domestication. In the phylogeny based on autosomal sequences, Przewalski's horses formed a monophyletic clade, separate from the Thoroughbred domestic horse lineage. Our results suggest that Przewalski's horses have ancient origins and are not the direct progenitors of domestic horses. The analysis of the vast amount of sequence data presented here suggests that Przewalski's and domestic horse lineages diverged at least 0.117 Ma but since then have retained ancestral genetic polymorphism and/or experienced gene flow.

Sequence variation and genetic diversity in the giant panda

About 336-444 bp mitochondrial D-loop region and tRNA gene were sequenced for 40 individuals of the giant panda which were collected from Mabian, Meigu, Yuexi, Baoxing, Pingwu, Qingchuan, Nanping and Baishuijiang, respectively. 9 haplotypes were found in 21 founders. The results showed that the giant panda has low genetic variations, and that there is no notable genetic isolation among geographical populations. The ancestor of the living giant panda population perhaps appeared in the late Pleistocene, and unfortunately, might have suffered bottleneck attacks. Afterwards, its genetic diversity seemed to recover to some extent.

Human and great ape red blood cells differ in plasmalogen levels and composition

Background

Plasmalogens are ether phospholipids required for normal mammalian developmental, physiological, and cognitive functions. They have been proposed to act as membrane antioxidants and reservoirs of polyunsaturated fatty acids as well as influence intracellular signaling and membrane dynamics. Plasmalogens are particularly enriched in cells and tissues of the human nervous, immune, and cardiovascular systems. Humans with severely reduced plasmalogen levels have reduced life spans, abnormal neurological development, skeletal dysplasia, impaired respiration, and cataracts. Plasmalogen deficiency is also found in the brain tissue of individuals with Alzheimer disease.

Results

In a human and great ape cohort, we measured the red blood cell (RBC) levels of the most abundant types of plasmalogens. Total RBC plasmalogen levels were lower in humans than bonobos, chimpanzees, and gorillas, but higher than orangutans. There were especially pronounced cross-species differences in the levels of plasmalogens with a C16:0 moiety at the sn-1 position. Humans on Western or vegan diets had comparable total RBC plasmalogen levels, but the latter group showed moderately higher levels of plasmalogens with a C18:1 moiety at the sn-1 position. We did not find robust sex-specific differences in human or chimpanzee RBC plasmalogen levels or composition. Furthermore, human and great ape skin fibroblasts showed only modest differences in peroxisomal plasmalogen biosynthetic activity. Human and chimpanzee microarray data indicated that genes involved in plasmalogen biosynthesis show cross-species differential expression in multiple tissues.

Conclusion

We propose that the observed differences in human and great ape RBC plasmalogens are primarily caused by their rates of biosynthesis and/or turnover. Gene expression data raise the possibility that other human and great ape cells and tissues differ in plasmalogen levels. Based on the phenotypes of humans and rodents with plasmalogen disorders, we propose that cross-species differences in tissue plasmalogen levels could influence organ functions and processes ranging from cognition to reproduction to aging.

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.

Chromosome painting in Tragulidae facilitates the reconstruction of Ruminantia ancestral karyotype

Although Tragulidae, as the basal family in Ruminantia phylogenetic tree, is the key taxon for understanding the early chromosome evolution of extant ruminants, comparative molecular cytogenetic data on the tragulids are scarce. Here, we present the first genome-wide comparative map of the Java mouse deer (Tragulus javanicus, Tragulidae) revealed by chromosome painting with human and dromedary probes. Together with the published comparative maps of major representative cetartiodactyl species established with the same set of probes, our results allowed us to reconstruct a 2n = 48 Ruminantia ancestral karyotype, which is similar to the cetartiodactyl ancestral karyotype. The karyotype evolution of T. javanicus has involved multiple rearrangements, most of which appear to be apomorphic and have not found in karyotype evolution of pecoran species (i.e., Ruminantia excluding Tragulidae). The rate of chromosome evolution of the mouse deer was rather low-0.4 R/Ma, while the estimated tempo of chromosome changes on the lineages leading from Cetartiodactyla ancestor to Ruminantia and from Ruminantia to Pecora were roughly the same (about 1.2 R/Ma).