Phylogeny of Muntiacus (Cervidae) based on mitochondrial DNA restriction maps

The systematics of the Cervidae: a total evidence approach

Systematic relationships of cervids have been controversial for decades. Despite new input from molecular systematics, consensus could only be partially reached. The initial, gross (sub) classification based on morphology and comparative anatomy was mostly supported by molecular data. The rich fossil record of cervids has never been extensively tested in phylogenetic frameworks concerning potential systematic relationships of fossil cervids to extant cervids. The aim of this work was to investigate the systematic relationships of extant and fossil cervids using molecular and morphological characters and make implications about their evolutionary history based on the phylogenetic reconstructions. To achieve these objectives, molecular data were compiled consisting of five nuclear markers and the complete mitochondrial genome of 50 extant and one fossil cervids. Several analyses using different data partitions, taxon sampling, partitioning schemes, and optimality criteria were undertaken. In addition, the most extensive morphological character matrix for such a broad cervid taxon sampling was compiled including 168 cranial and dental characters of 41 extant and 29 fossil cervids. The morphological and molecular data were analysed in a combined approach and other comprehensive phylogenetic reconstructions. The results showed that most Miocene cervids were more closely related to each other than to any other cervids. They were often positioned between the outgroup and all other cervids or as the sister taxon to Muntiacini. Two Miocene cervids were frequently placed within Muntiacini. Plio- and Pleistocene cervids could often be affiliated to Cervini, Odocoileini or Capreolini. The phylogenetic analyses provide new insights into the evolutionary history of cervids. Several fossil cervids could be successfully related to living representatives, confirming previously assumed affiliations based on comparative morphology and introducing new hypotheses. New systematic relationships were observed, some uncertainties persisted and resolving systematics within certain taxa remained challenging.

Comparative sequence analyses reveal sites of ancestral chromosomal fusions in the Indian muntjac genome

Comparative mapping and sequencing was used to characterize the sites of ancestral chromosomal fusions in the Indian muntjac genome.

Background

Indian muntjac (Muntiacus muntjak vaginalis) has an extreme mammalian karyotype, with only six and seven chromosomes in the female and male, respectively. Chinese muntjac (Muntiacus reevesi) has a more typical mammalian karyotype, with 46 chromosomes in both sexes. Despite this disparity, the two muntjac species are morphologically similar and can even interbreed to produce viable (albeit sterile) offspring. Previous studies have suggested that a series of telocentric chromosome fusion events involving telomeric and/or satellite repeats led to the extant Indian muntjac karyotype.

Results

We used a comparative mapping and sequencing approach to characterize the sites of ancestral chromosomal fusions in the Indian muntjac genome. Specifically, we screened an Indian muntjac bacterial artificial-chromosome library with a telomere repeat-specific probe. Isolated clones found by fluorescence in situ hybridization to map to interstitial regions on Indian muntjac chromosomes were further characterized, with a subset then subjected to shotgun sequencing. Subsequently, we isolated and sequenced overlapping clones extending from the ends of some of these initial clones; we also generated orthologous sequence from isolated Chinese muntjac clones. The generated Indian muntjac sequence has been analyzed for the juxtaposition of telomeric and satellite repeats and for synteny relationships relative to other mammalian genomes, including the Chinese muntjac.

Conclusions

The generated sequence data and comparative analyses provide a detailed genomic context for seven ancestral chromosome fusion sites in the Indian muntjac genome, which further supports the telocentric fusion model for the events leading to the unusual karyotypic differences among muntjac species.

Population structure and gene flow among wild populations of the black muntjac (Muntiacus crinifrons) based on mitochondrial DNA control region sequences

The black muntjac (Muntiacus crinifrons) is a rare species endemic to China and its current distribution is confined to partial mountain ranges in eastern China. To assess the population structure and gene flow among the extant populations, we sequenced a portion of the mitochondrial control region (424 bp) for 47 samples collected from its current three large populations (Huangshan, Tiammushan and Suichang). A total of 18 unique haplotypes were defined based on 22 polymorphic sites. Of these haplotypes, 15 ones were population-specific and only one haplotype was shared among the three populations. Significant genetic differentiation was detected between Suichang and Huangshan populations (ø(ST)=0.1677, P<0.001) or between Suichang and Tianmushan populations (ø(ST)=0.2002, P<0.001), indicating that the Suichang population may be spatially structured from other two populations along matriline. The Mantel test revealed that this significant differentiation was not driven by geographic distance (P=0.14), implying that genetic divergence of Suichang population might result from recent human disturbances. Phylogenetic analyses suggested the mitochondrial control region haplotypes were split into two well divergent clades (Clade I and Clade II). Interestingly, the two distinct haplotype clades were found to coexist in Suichang area. The nested clade analysis revealed a significant phylogeographic structure among the black muntjac populations (total cladogram: chi2=18.68; P<0.001), which was inferred to result from past fragmentation followed by range expansion. The population expansion was supported by the analysis of mismatch distribution and the tests of neutrality. Therefore, we suggest that the coexistence of distinct haplotypes in Suichang population was induced by historical population expansion after fragmentation and that the current genetic differentiation should be attributed to the reduction of female-mediated gene flow due to recent habitat fragmentation and subsequent loss.

Microsatellite analysis of genetic variation and population subdivision for the black muntjac, Muntiacus crinifrons

The black muntjac (Muntiacus crinifrons) is a rare deer found only in a restricted region in east China. Recent studies of mitochondrial DNA diversity have shown a markedly low level of nucleotide diversity for the species, and the Suichang population was genetically differentiated from the two other populations, in Huangshan and Tianmushan mountains. In this study, we extended the analysis of genetic diversity and population subdivision for the black muntjac using data from 11 highly polymorphic nuclear DNA microsatellite loci. Contrary to the results based on mtDNA data, the microsatellite loci revealed that the black muntjac retained a rather high nuclear genetic diversity (overall average H (E) = 0.78). Nevertheless, both types of markers supported the idea that the extant black muntjac population is genetically disrupted (overall phi (ST) = 0.16 for mtDNA and overall F (ST) = 0.053 for microsatellite, both P < 0.001). The correlation between genetic differentiation and geographic distance was not significant (Mantel test; P > 0.05), implying that the patterns of genetic differentiation observed in this study might result from recent habitat fragmentation or loss. Based on the results from the mtDNA and nuclear DNA data sets, two management units were defined for the species, Huangshan/Tianmushan and Suichang. We also recommend that a new captive population be established with individuals from the Suichang region as a founder source.

High-density comparative BAC mapping in the black muntjac (Muntiacus crinifrons): molecular cytogenetic dissection of the origin of MCR 1p+4 in the X1X2Y1Y2Y3 sex chromosome system

The black muntjac (Muntiacus crinifrons, 2n = 8[female symbol]/9[male symbol]) is a critically endangered mammalian species that is confined to a narrow region of southeastern China. Male black muntjacs have an astonishing X1X2Y1Y2Y3 sex chromosome system, unparalleled in eutherian mammals, involving approximately half of the entire genome. A high-resolution comparative map between the black muntjac (M. crinifrons) and the Chinese muntjac (M. reevesi, 2n = 46) has been constructed based on the chromosomal localization of 304 clones from a genomic BAC (bacterial artificial chromosome) library of the Indian muntjac (M. muntjak vaginalis, 2n = 6[female symbol]/7[male symbol]). In addition to validating the chromosomal homologies between M. reevesi and M. crinifrons defined previously by chromosome painting, the comparative BAC map demonstrates that all tandem fusions that have occurred in the karyotypic evolution of M. crinifrons are centromere-telomere fusions. The map also allows for a more detailed reconstruction of the chromosomal rearrangements leading to this unique and complex sex chromosome system. Furthermore, we have identified 46 BAC clones that could be used to study the molecular evolution of the unique sex chromosomes of the male black muntjacs.

Mitochondrial DNA genetic diversity of black muntjac (Muntiacus crinifrons), an endangered species endemic to China

Genetic diversities based on the mtDNA control region were measured for both a wild population (n = 26) and a captive population (n = 18) of the black muntjac. In total, nine haplotypes were obtained from 44 samples. The wild population exhibited a low nucleotide diversity (pi = 0.00562), which suggests that the black muntjac had a small effective population size historically. In contrast to its low nucleotide diversity, haplotype diversity (h = 0.862) of the wild population was relatively high. Haplotype distribution among local samples shows a distinct difference. As anticipated because of the paucity of available founders, nucleotide diversity (pi = 0.00214) of the captive population was very low. Additionally, a high degree of haplotype identity and an obvious haplotype frequency bias was revealed in the captive population, which implies that the current breeding program should be readjusted to balance distributions of haplotypes, and some new founders should be introduced to the captive population to alleviate potential inbreeding depression.

A complete phylogeny of the whales, dolphins and even-toed hoofed mammals (Cetartiodactyla)

Despite the biological and economic importance of the Cetartiodactyla, the phylogeny of this clade remains controversial. Using the supertree approach of matrix representation with parsimony, we present the first phylogeny to include all 290 extant species of the Cetacea (whales and dolphins) and Artiodactyla (even-toed hoofed mammals). At the family-level, the supertree is fully resolved. For example, the relationships among the Ruminantia appear as (((Cervidae, Moschidae) Bovidae) (Giraffidae, Antilocapridae) Tragulidae). However, due to either lack of phylogenetic study or contradictory information, polytomies occur within the clades Sus, Muntiacus, Cervus, Delphinidae, Ziphiidae and Bovidae. Complete species-level phylogenies are necessary for both illustrating and analysing biological, geographical and ecological patterns in an evolutionary framework. The present species-level tree of the Cetartiodactyla provides the first opportunity to examine comparative hypotheses across entirely aquatic and terrestrial species within a single mammalian order.

Karyotypic evolution of a novel cervid satellite DNA family isolated by microdissection from the Indian muntjac Y-chromosome

A minilibrary was constructed from DOP-PCR products using microdissected Y-chromosomes of Indian muntjac as DNA templates. Two microclones designated as IM-Y4-52 and IM-Y5-7 were obtained from negative screening of all three cervid satellite DNAs (satellites I, II, and IV). These two microclones were 295 and 382 bp in size, respectively, and shared approximately 70% sequence homology. Southern blot analysis showed that the IM-Y4-52 clone was repetitive in nature with an approximately 0.32-kb register in HaeIII digest. Sequence comparison revealed no similarities to DNA sequences deposited in the GenBank database, suggesting that the microclone sequences were from a novel satellite DNA family designated as cervid satellite V. A subclone of an Indian muntjac BAC clone which screened positive for IM-Y4-52 had a 3,325-bp insert containing six intact monomers, four deleted monomers, and two partial monomers. The consensus sequence of the monomer was 328 bp in length and shared more than 80% sequence homology with every intact monomer. A zoo blot study using IM-Y4-52 as a probe showed that the strong hybridization with EcoRI digested male genomic DNA of Indian muntjac, Formosan muntjac, Chinese muntjac, sambar deer, and Chinese water deer. Female genomic DNA of Indian muntjac, Chinese water deer, and Formosan muntjac also showed positive hybridization patterns. Satellite V was found to specifically localize to the Y heterochromatin region of the muntjacs, sambar deer, and Chinese water deer and to chromosome 3 of Indian muntjac and the X-chromosome of Chinese water deer.

Estimation of the highest chromosome number of eukaryotes based on the minimum interaction theory

According to the minimum interaction theory, the chromosome evolution of eukaryotes proceeds as a whole toward increasing the chromosome number. This raises the following two questions: what was the starting chromosome number of eukaryotes and does the chromosome number increase infinitely? We attempted to provide a theoretical framework to resolve these questions. We propose that the species with n=2 observed in Protozoa, Platyhelminthes, Annelid, Algae, Fungi and higher plants would be chromosomal relicts conserving the karyotypes of ancestral eukaryotes. We also propose that the ideal highest number of eukaryotes (n(max)) can be given by an inverse of the minimum terminal interference distance (It(min)) in crossing-over (n(max)=100/It(min)). AsIt(min) =0.6 in mammals, n(max) approximately 166. On the other hand, the value estimated by computer simulations is somewhat lower with n(max)=133-138. Our arguments can be applied to other eukaryotes, if they have a localized centromere and the ratio of total synaptonemal complex/nuclear volume is comparable to that of mammals. We revealed that the index of gene shuffling per karyotypes (G) by means of the total number of gamete types with different gene combinations can be formulated asG =2(n+Fxi), where Fxi means interstitial chiasma frequency per cell corresponding to crossing-over mediated by the recombination nodule. The Fxi value increases in proportion to the n value in areas where n<40, but decreases gradually when n>40 and becomes zero when n>83. Therefore, in the ultimate karyotype with n(max)=166, FXi=0 andG =2(n)=2(166), where gene shuffling is guaranteed by the random orientation of chromosomes at the equatorial plate of meiotic metaphase I.

Integrative study on chromosome evolution of mammals, ants and wasps based on the minimum interaction theory

There is well-known evidence that in many eukaryotes, different species have different karyotypes (e.g. n=1-47 in ants and n=3-51 in mammals). Alternative (fusion and fission) hypotheses have been proposed to interpret this chromosomal diversity. Although the former has long been accepted, accumulating molecular genetics evidence seems to support the latter. We investigated this problem from a stochastic viewpoint using the Monte Carlo simulation method under the minimum interaction theory. We found that the results of simulations consistently interpreted the chromosomal diversity observed in mammals, ants and wasps, and concluded that chromosome evolution tends to evolve as a whole toward increasing chromosome numbers by centric fission. Accordingly, our results support the fission hypothesis. We discussed the process of chromosome evolution based on the latest theory of the molecular structure of chromosomes, and reconfirmed that the fission burst is the prime motive force in long-term chromosome evolution, and is effective in minimizing the genetic risks due to deleterious reciprocal translocations and in increasing the potential of genetic divergence. Centric fusion plays a biological role in eliminating heterochromatin (C-bands), but is only a local reverse flow in contrast to the previously held views.

Interstitial colocalization of two cervid satellite DNAs involved in the genesis of the Indian muntjac karyotype

A number of repetitive DNA clones were generated from PCR amplifications of Indian muntjac genomic DNA using primer sequences derived from a white tailed deer satellite II DNA sequence. One clone (Mmv-0.7) was characterized and shown to be a cervid satellite II DNA clone. Multiple colored FISH studies with cervid satellite I (C5) and this satellite II clone (Mmv-0.7) to Chinese muntjac metaphase chromosomes localized both satellite DNAs at the pericentromeric regions of all chromosomes except for chromosome 3 and the Y chromosome, whereas chromosome 3 exhibited pericentromeric satellite II DNA only. Where distinguishable, the pericentromeric satellite II signals appeared terminally oriented with respect to satellite I. Six pairs of Chinese muntjac autosomes had interstitial satellite I sites with four of these autosomal pairs (chromosomes 1, 2 and two other smaller autosomal pairs) also exhibiting interstitial satellite II signals. An interstitial site on the X chromosome was found to have satellite II signals. For the Indian muntjac chromosomes, FISH studies revealed a pericentromeric hybridization for satellites I and II as well as 27 distinct interstitial hybridization sites, each having at least one of the satellite DNAs. These data were used to more precisely define the chromosome fusion-associated breakpoints that presumably led to the formation of the present-day Indian muntjac karyotype. It further hints at the possibility that the Indian muntjac karyotype may have evolved directly from a 2n = 70 ancestral karyotype rather than from an intermediate 2n = 46 Chinese muntjac-like karyotype.

Rapid and parallel chromosomal number reductions in muntjac deer inferred from mitochondrial DNA phylogeny

Muntjac deer (Muntiacinae, Cervidae) are of great interest in evolutionary studies because of their dramatic chromosome variations and recent discoveries of several new species. In this paper, we analyze the evolution of karyotypes of muntjac deer in the context of a phylogeny which is based on 1,844-bp mitochondrial DNA sequences of seven generally recognized species in the muntjac subfamily. The phylogenetic results support the hypothesis that karyotypic evolution in muntjac deer has proceeded via reduction in diploid number. However, the reduction in number is not always linear, i.e., not strictly following the order: 46-->14/13-->8/9-->6/7. For example, Muntiacus muntjak (2n = 6/7) shares a common ancestor with Muntiacus feae (2n = 13/14), which indicates that its karyotype was derived in parallel with M. feae's from an ancestral karyotype of 2n >/= 13/14. The newly discovered giant muntjac (Muntiacus vuquangensis) may represent another parallel reduction lineage from the ancestral 2n = 46 karyotype. Our phylogenetic results indicate that the giant muntjac is relatively closer to Muntiacus reevesi than to other muntjacs and may be placed in the genus Muntiacus Analyses of sequence divergence reveal that the rate of change in chromosome number in muntjac deer is one of the fastest in vertebrates. Within the muntjac subfamily, the fastest evolutionary rate is found in the Fea's lineage, in which two species with different karyotypes diverged in around 0.5 Myr.

Zoo-fluorescence in situ hybridization analysis of human and Indian muntjac karyotypes (Muntiacus muntjak vaginalis) reveals satellite DNA clusters at the margins of conserved syntenic segments

Zoo-fluorescence in situ hybridization (FISH) with human whole chromosome-specific paint probes revealed extensive homoeologies between Indian muntjac (2n=6, 7 female, male) and human karyotypes (2n=46). Forty-two conserved syntenic segments, corresponding to all human chromosomes except the Y chromosome, produced a near-complete coverage of the muntjac complement and revealed margins of interspecific segmental homoeology. To test the hypothesis that interstitial satellite DNA loci, illuminated by a Chinese muntjac C5-satellite probe in Indian muntjac chromosome arms, mark ancestral fusion points (Lin CC, Sasi R, Fan YS, Chen Z-Q (1991) New evidence for tandem chromosome fusions in the karyotypic evolution of the Asian muntjacs. Chromosoma 101: 19-24), we combined Zoo-FISH with C5 satellite mapping. Twenty-six interstitial satellite DNA loci were detected in the haploid Indian muntjac genome and were found to co-localize with the margins of conserved human/Indian muntjac syntenic segments. These results were confirmed by two-colour FISH and are in accordance with the tandem fusion hypothesis for Indian muntjac chromosomes. Furthermore, conserved syntenic segment combinations detected in pig, cattle and Indian muntjac Zoo-FISH maps reveal ancestral artiodactyl chromosomes.

A reappraisal of the tandem fusion theory of karyotype evolution in Indian muntjac using chromosome painting

We have tested the tandem fusion hypothesis of the origin of the Indian muntjac karyotype (2n = 6/7) by using reciprocal chromosome painting between the Indian muntjac, Chinese muntjac (n = 46) and brown brocket deer (2n = 70 + 3B) with chromosome-specific paint probes derived from flow-sorted chromosomes of these three deer species. Our results have shown that the euchromatic blocks of all chromosome arms of the brown brocket deer have been conserved apparently unchanged in number and content in the Indian muntjac. While confirming the conservation in toto of most of Chinese muntjac euchromatin in the karyotype of the Indian muntjac, we demonstrate that the synteny of chromosomes 1, 2, 3, 4 and 5 of the Chinese muntjac has been disrupted by chromosome rearrangements other than fusions. This indicates that the present karyotype of the Indian muntjac cannot be reconstructed from the hypothetical Chinese muntjac-like 2n = 46 ancestral karyotype exclusively by chromosome fusions. Furthermore, we have shown that the breakpoints of these rearrangements appear to have occurred near to the fusion points formed during the origin of the 2n = 46 karyotype of the Chinese muntjac from a 2n = 70 karyotype, which is believed to be ancestral for the family Cervidae. Moreover, we substantiate that on the Indian muntjac chromosomes, the C5 probe, which is derived from the centromeric satellite sequences of the Chinese muntjac, maps to the putative fusion points determined by comparative chromosome painting and presumably represents the remnants of ancestral centromeric sequences.