- Research article
- Open Access
Comparative mitochondrial genomics and phylogenetic relationships of the Crossoptilon species (Phasianidae, Galliformes)
BMC Genomics volume 16, Article number: 42 (2015)
Phasianidae is a family of Galliformes containing 38 genera and approximately 138 species, which is grouped into two tribes based on their morphological features, the Pheasants and Partridges. Several studies have attempted to reconstruct the phylogenetic relationships of the Phasianidae, but many questions still remain unaddressed, such as the taxonomic status and phylogenetic relationships among Crossoptilon species. The mitochondrial genome (mitogenome) has been extensively used to infer avian genetic diversification with reasonable resolution. Here, we sequenced the entire mitogenomes of three Crossoptilon species (C. harmani, C. mantchuricum and C. crossoptilon) to investigate their evolutionary relationship among Crossoptilon species.
The complete mitogenomes of C. harmani, C. mantchuricum and C. crossoptilon are 16682 bp, 16690 bp and 16680 bp in length, respectively, encoding a standard set of 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and a putative control region. C. auritum and C. mantchuricum are more closely related genetically, whereas C. harmani is more closely related to C. crossoptilon. Crossoptilon has a closer relationship with Lophura, and the following phylogenetic relationship was reconstructed: ((Crossoptilon + Lophura) + (Phasianus + Chrysolophus)). The divergence time between the clades C. harmani-C. crossoptilon and C. mantchuricum-C. auritum is consistent with the uplift of the Tibetan Plateau during the Tertiary Pliocene. The Ka/Ks analysis showed that atp8 gene in the Crossoptilon likely experienced a strong selective pressure in adaptation to the plateau environment.
C. auritum with C. mantchuricum and C. harmani with C. crossoptilon form two pairs of sister groups. The genetic distance between C. harmani and C. crossoptilon is far less than the interspecific distance and is close to the intraspecific distance of Crossoptilon, indicating that C. harmani is much more closely related to C. crossoptilon. Our mito-phylogenomic analysis supports the monophyly of Crossoptilon and its closer relationship with Lophura. The uplift of Tibetan Plateau is suggested to impact the divergence between C. harmani-C. crossoptilon clade and C. mantchuricum-C. auritum clade during the Tertiary Pliocene. Atp8 gene in the Crossoptilon species might have experienced a strong selective pressure for adaptation to the plateau environment.
Crossoptilon, belonging to Phasianidae in Galliformes, is a rare but important genus endemic to China. The four previous recognised Crossoptilon species are C. harmani, C. mantchuricum, C. crossoptilon and C. auritum [1,2]. C. harmani is only distributed in a typical alpine and taiga habitats at elevations of 3500–3900 meters of southeastern Tibet. C. mantchuricum is mainly distributed in Lvliang Mountains of Shanxi Prov., Xiaowutai Mountains of Hebei Prov., Dongling Mountains of Beijing and some other local areas . C. crossoptilon is only present in western Sichuan Prov., northwestern Yunnan Prov., southeastern Qinghai Prov. and eastern Tibet at high mountain areas [3,4]. It is common at high altitudes (3000–4300 m) and exists in coniferous forests, mixed broadleaf-conifer forests and alpine scrubs . C. auritum is only encountered in the mountainous regions of Qinghai, Gansu and Sichuan provinces and Ningxia Hui Autonomous Region .
The taxonomic status and phylogenetic position of C. harmani are controversial [4-10]. Whereas some studies have proposed that it is an independent species [1,2,6,11], others have suggested that it is a subspecies of C. crossoptilon [4,5,12-14]. The C. mantchuricum taxon is believed to have diverged first, and C. crossoptilon and C. auritum are considered more closely related than the other species [7,9,10]. Other studies have suggested that C. crossoptilon is a relatively original species and that C. auritum and C. mantchuricum are more closely related genetically [1,4-6,8,12,14].
Previous studies on Crossoptilon species have focused on their distribution, ecology and behaviour. Some studies have used single gene or partial sequences in reconstructing their phylogenetic relationships , but a comprehensive study of the entire mitogenomes has rarely been attempted. To reconstruct the phylogenetic relationships among Crossoptilon and other Phasianidae species, the entire mitogenomes of C. harmani, C. mantchuricum and C. crossoptilon were sequenced, and thus four species of the genus Crossoptilon were thoroughly compared. Furthermore, the estimated divergence times of the four Crossoptilon species were calculated, and a Ka/Ks analysis was used to estimate the adaptation of the Crossoptilon mitogenome to different environments.
Sample collection and DNA extraction
Samples of C. harmani were collected from the Nagqu area in Tibet, China, and voucher specimen was deposited in Institute of Zoology, Shaanxi Normal University; C. mantchuricum and C. crossoptilon were collected from the Beijing Zoo, and samples were obtained from bird specimen collections of the National Zoological Museum, Institute of Zoology, Chinese Academy of Sciences. The collection was under the permit from Forestry Department and conformed to the National Wildlife Conservation Law in China. No living animal experiments were conducted in the current research. All the samples were preserved in 100% ethanol and stored at −20°C. The total genomic DNA was extracted from the liver/muscle tissue using the standard phenol/chloroform method .
PCR amplification and sequencing
The PCRs were performed under the following conditions, 2 min initial denaturation at 93°C, 40 cycles of: 10 s denaturation at 92°C, 30 s annealing at 58–53°C, 10 min elongation at 68°C in the preliminary 20 cycles, and 10 s denaturation at 92°C, 30 s annealing at 53°C, and 10 min elongation at 68°C with 20 s per cycle added to the elongation step in the succeeding 20 cycles, and finally an extension for 7 min at 68°C. The amplifications were performed in 15 μL reactions containing 2.4 μL of 2.5 mM dNTPs, 2.1 μL of each primer at 10 μΜ, 1.5 μl of 10× LA PCR Buffer I (Mg2+-free), 1.5 μL of 25 mM MgCl2, 1 μL of DNA template, 0.18 μL of 5 U/μl LA Taq polymerase (Takara, Dalian, China) and 4.22 μL ddH2O.
The sequences of the primers used for PCR amplification and sequencing of the mitochondrial genes were obtained from Sorenson (2003)  with minor changes (Additional file 1). The mitogenome of C. harmani was amplified in seven parts and the gaps were bridged using other adjacent primers. The PCR products were purified using the DNA Agarose Gel Extraction Kit (Bioteke, Beijing, China) after separation by electrophoresis on a 0.8% agarose gel. After separation and purification, the PCR products were sequenced by Sangon Biotech (Shanghai) Co., LTD. using the Primer-Walking method. The C. mantchuricum and C. crossoptilon mitogenomes were amplified in 12 or 13 fragments, and sequencing was performed using the Illumina Hiseq2000 high-throughput sequencing system of Shenzhen Huada Gene Technology Co., LTD. The gaps in the assembly after high-throughput sequencing were filled in by direct sequencing using the ABI 3730 DNA sequencer by Sangon Biotech (Shanghai) Co., LTD., using adjacent PCR primers.
Gene identification and genome analyses
The Staden sequence analysis package  was used for sequence assembly and annotation of C. harmani mitogenome. The complete genome assemblies for the mitogenomes of C. mantchuricum and C. crossoptilon were performed using the SOAP de novo software. Most tRNA genes were identified using tRNAscan-SE 1.21  under the ‘tRNAscan only’ search mode, with the vertebrate mitochondrial genetic code and ‘mito/chloroplast’ source. The protein-coding genes (PCGs), rRNA genes and the remaining putative tRNA genes that were not identified by tRNAscan-SE were identified by sequence comparison with other Galliformes species. The rrnS secondary structure of Crossoptilon was predicted based on the structure of Gallus gallus and Anas platyrhynchos obtained from the Comparative RNA Web (CRW) , and Pseudopodoces humilis (now as Parus humilis) structure . The rrnL secondary structure was predicted based on the structure of Xenopus laevis obtained from the CRW database, Bos taurus  and P. humilis structures . The RNAstructure software was used to identify and draw potential secondary structures in the single-stranded control region. The nucleotide compositions of the mitogenomes and amino acid information were analysed using MEGA 4.1 .
Along with the entire mitogenomes obtained in this study, 42 Galliformes sequences were used in the phylogenetic analysis, including two outgroups (Numida meleagris and Alectura lathami). The DNA sequences of the other species used in the phylogenetic analyses were downloaded from GenBank (the accession numbers and key information are shown in Additional file 2). The tRNA and rRNA genes and the CR were individually aligned using ClustalX 1.83  with the default settings. All 13 protein-coding genes were translated into amino acids, and then aligned using MEGA 4.1  with default parameters for each gene, and finally retranslated into nucleotide sequences.
Phylogenetic analyses of Phasianidae
Datasets containing 13 protein-coding genes (PCG) and all 37 genes plus the control region (mitogenome) were used to study the phylogenetic relationships within Phasianidae. Phylogenetic analysis based on nucleotide sequences was performed using PAUP*4.0b10 for the Maximum parsimony (MP) method , RAxML-7.0.3 for maximum likelihood (ML)  and MrBayes 3.1.2 for Bayesian inference (BI) . For ML and BI analyses, models of the concatenated nucleotide sequences datasets were assessed independently using AICc in MrModeltest2.2 . The best fit model GTR + I + G was chosen for the likelihood and Bayesian analyses. A consensus tree was generated for MP analysis under the majority rule. The reliability of the clades in the phylogenetic trees was assessed by bootstrap probabilities (BSP) computed using 1000 replicates, with random addition for each bootstrap replicate. The 1000 replicates bootstrap support was also performed in the ML analysis. Bayesian analysis with Markov Chain Monte Carlo sampling was run for 1000000 generations saving a tree every 100 generations, with one cold and three heated chains, and the burn-in time was determined by the time to convergence of the likelihood scores. The Bayesian posterior probabilities (BPP) were estimated on a 50% majority rule consensus tree of the remaining trees.
We examined the performance of individual genes and datasets based on nucleotides; PBS analyses were performed in the program combining TreeRot.v3  and PAUP*4.0b10 . The following datasets were used for analysis: 13 protein-coding genes, rrnS, rrnL, CR partitions, the first, second and third codons of PCG, the three tRNA gene cluster (IQM, WANCY and HSL), ATP (atp6 + atp8), COX (cox1 + cox2 + cox3) and NADH (nad1 + nad2 + nad3 + nad4 + nad4L + nad5 + nad6).
The MEGA 4.1  was used to calculate the pairwise genetic distance for four Crossoptilon species with default parameters. The mitogenome aligned data and four single genes (nad2, CR, cytb and rrnS) obtained from GenBank (Additional file 3) were used to calculate the genetic distances. These genes were aligned singly, and be adjusted to consistent sequence lengths manually.
Divergence time estimates focused on Crossoptilon
Along with the PCG dataset obtained in this study, 42 Galliformes sequences were used to estimate the divergence time of the Crossoptilon species. The divergence time of Crossoptilon species was calculated using the Bayesian procedure implemented in BEAST v. 1.7.2 [29,30]. A relaxed clock was used with rates complying with a log-normal distribution . The GTR + I + G model and a Yule prior were used in the analysis. The calibration points were based on the fossil records showing that stem Numididae-Phasianidae split at 50–54 Mya (million years ago) ; Arborophila rufipectus diverged from the other lineages in the Galliformes around 39 Mya ; Coturnix-Gallus split at 35 Mya [34-36]. The results of runs of 10 million generations were used after a burn-in of 100.
Ka and Ks analysis
To better understand the evolution at the DNA level and the role of selection in the four Crossoptilon species, we calculated the nonsynonymous and synonymous substitution rates using Kaks_calculator 2.0  for six groups [C. harmani-C. mantchuricum (C.har-C.man), C. mantchuricum-C. crossoptilon (C.man-C.cro), C. harmani-C. crossoptilon (C.har-C.cro), C. harmani-C. auritum (C.har-C.aur), C. mantchuricum-C. auritum (C.man-C.aur), and C. crossoptilon-C. auritum (C.cro-C.aur)]. The ratio of nonsynonymous substitution rate (Ka) to synonymous substitution rate (Ks) is widely used as an indicator of selective pressure at the sequence level among different species. It is commonly accepted that Ka > Ks, Ka = Ks, and Ka < Ks generally indicate positive selection, neutral mutation, and negative selection, respectively [38,39]. To calculate Ka, Ks and Ka/Ks, a model averaging method was selected. This method includes 14 different models for calculation and derived the average values for Ka, Ks, and Ka/Ks . The genetic code selected was the ‘vertebrate mitochondrial code’. To further study the selective pressure acted on each protein-coding gene in the genus Crossoptilon, CodeML in PAMLX software  was used to find sites under strong selective pressure. The secondary structure analysis of amino acid was performed by using an online software TOPCONS .
Comparison of the two sequencing methods used in this study revealed that high-throughput sequencing has greater coverage and accuracy relative to standard sequencing, albeit with higher cost. Using standard sequencing, we obtained 62 effective sequences with 2.63-fold coverage. In contrast, the high-throughput sequencing yielded effective assembly data with sequence depths (X) of 7604.35 for C. mantchuricum and 7810.82 for C. crossoptilon after filtering out some reads, such as low-quality or adapter-sequence-polluted reads.
The complete mitogenomes of C. harmani, C. mantchuricum and C. crossoptilon are 16682 bp, 16690 bp and 16680 bp in length, respectively. The annotated genomes have been deposited in the GenBank database (accession numbers: KP259806-KP259808). The mitochondrial genes of three Crossoptilon species are coded on the H-strand, except for one protein-coding gene (nad6) and eight tRNA genes (trnQ, trnA, trnN, trnC, trnY, trnS (UCN), trnP, and trnE) (Table 1). Consistent with previous results on the Phasianidae mitogenomes, the mitochondrial sequences of four Crossoptilon species (C. harmani, C. mantchuricum, C. crossoptilon and C. auritum) are biased toward adenine and thymine (54.3%, 54.2%, 54.3%, and 54.2%, respectively). The GC skew in Crossoptilon mitogenomes is similar, with parts of the CR and rrnS containing significantly higher skews than other regions. The CR, trnF, rrnS, trnV, rrnL and cox1 regions have strong GC skews.
The 13 protein-coding genes of Crossoptilon genomes are similar to most other Phasianidae species, with nad5 and atp8 being the longest and shortest genes, respectively. The total length of the PCGs in each Crossoptilon species is 11358 bp after removing termination codons, containing approximately 3786 codons. The A + T content of the 13 PCGs is 53.4% in C. harmani and C. crossoptilon, and 53.3% in C. auritum and C. mantchuricum. Analysis of the base composition at each codon position of the concatenated PCGs shows that the second codon position has a higher A + T content (57.9% in C. auritum and C. mantchuricum, 58.1% in C. harmani, and 58.0% in C. crossoptilon, respectively) than the first and third codon positions. The amino acid frequencies in the Crossoptilon PCGs are similar, with Leu significantly more frequent than other amino acids. The different codon positions have the same base distributions, with C and G as the most and least frequent bases in the third codon, respectively.
All PCGs in Crossoptilon species start with the typical ATG codon, except the nad5 gene in C. mantchuricum and the cox1 gene in four Crossoptilon species, which start with GTG. Four types of stop codons are used by the coding genes, including TAA and TAG for most genes, AGG for cox1 in Crossoptilon and nad6 in C. harmani and C. crossoptilon, and an incomplete stop codon T- for cox3, nad4 and nad2 in four Crossoptilon species, respectively.
Similar to previously sequenced mitogenomes, the genomes sequenced in this study contain 2 rRNA genes encoding the small and large rRNA subunits, which are located between trnF and trnL(UUR) and separated by the trnV gene. The lengths and A + T contents of the rrnS and rrnL in the Crossoptilon genomes are within the range observed in other Phasianidae species. The rrnS contains three domains with 46 predicted stems, and the rrnL contains six domains with 59 stems (Additional files 4 and 5). The secondary structures of the rrnS in C. crossoptilon and C. harmani are identical; the rrnL secondary structures differ by only 2 bp in length. However, there are many differences in the rrnL secondary structures among the four Crossoptilon species, specifically in the loop near stem 44, which contains several substitutions and indels.
The complete mitogenome sequence contains 22 interspersed tRNA genes. All the tRNA sequences have the potential to fold into typical cloverleaf secondary structures except for trnS(AGY), which lacks the DHU arm (Additional file 6). The secondary structures of trnF, trnL(UUR), trnQ, trnM, trnC, trnY, trnS(UCN), trnH, trnL(CUN), trnT, and trnP are relatively conserved. However, the structure of trnS(AGY) is different in C. auritum relative to the other three Crossoptilon species (Additional file 6); the structure in C. harmani, C. mantchuricum and C. crossoptilon contain two additional bases (G and A) in the amino acid acceptor arm. The most frequent mismatch is G-U; other mismatches are also present, including U-U in trnM and trnG, C-C in trnD and trnL(UUR), C-U in trnF and trnI, and A-C in trnS(AGY). The C. crossoptilon and C. harmani tRNA structures are almost identical, with the exception of 1 bp. The A + T content of the tRNA genes is 57.6% in C. harmani and C. crossoptilon, 57.4% in C. mantchuricum and C. auritum.
The nucleotide composition of the control region in the Crossoptilon species has a bias against G, which is common in the mitogenome sense strand of vertebrates . The control region of Crossoptilon is located in the conserved position between trnE and trnF and the length (1146 nucleotides) is conserved in all four species. The control region contains three domains: the ETAS Domain I (nt 1–312), Central Conserved Domain II (nt 313–780) and CSB Domain III (nt 781–1146). The A + T content of the three domains is similar in all the Crossoptilon species; Domains I, II, and III have higher contents of C and A, T, and A, respectively, relative to other bases. Domain III contains a higher percentage of A than Domain II, and Domain II has the highest G content among the three domains. The distribution of variable sites and conserved sites suggests that Domain II has relatively more conserved sites and Domain I has more variable sites compared to other domains.
The ETAS Domain I can be divided into two parts: part A, from nt 1–163 and part B, from nt 164–312. There are two conserved blocks in part A, ETAS1 (nt 64–126) and ETAS2 (nt 124–163), which are similar to motifs previously identified in other avian and mammalian species. The first block is perfectly conserved among the Crossoptilon species and has sequence similarity to the “goose hairpin” described in some phasianids and Anas species [43-46] (Additional file 7). The secondary structure of this hairpin is determined by a stem of seven complementary Cs/Gs and a loop containing a TCCC motif also present in mammalian control region , which was associated experimentally with termination of H strands . There are two copies of TCCC located at nt 22–25 and 183–186. The second block, which is perfectly conserved among the Crossoptilon species, has sequence similarity to the mammalian TASs, including the highly conserved motif GTGCAT, which is present in all sequenced Phasianids and Anseriforms. The motif GYRCAT (Y = C/T; R = A/G) is widespread in Domain I of some mammalian control region , and it has been duplicated in the R1 repeats of many species, including the opossum  and several rodents [50,51]. Its functional importance is suggested by both comparative  and experimental  data. Single-stranded Domain I can form potential secondary cloverleaf structures. The first 190 nucleotides of Domain I of the Crossoptilon species can form stable secondary structures (Additional file 7). The configurations of the cloverleaves are similar in C. harmani and C. crossoptilon and vary slightly in C. mantchuricum and C. auritum; the stem-loop structures of the cloverleaves always include the “goose hairpin”. The control region of C. harmani and C. crossoptilon have similar structures with five stems, while C. mantchuricum and C. auritum have only four.
The Central Conserved Domain II of the control region of the Crossoptilon species includes 468 conserved nucleotides, from the hyper-variable site at position 313 at the end of the ETAS domain to the beginning of the putative OH sequence block at position 781. Several conserved blocks of the central conserved domain are similar to the F (nt 350–377), E (nt 395–414), D (nt 455–479), and C (nt 504–529) boxes of vertebrates and other avian species [43,52], and a bird similarity box is also present.
The CSB Domain III is highly variable and has sequences similar to mammalian CSB1. A poly(C) sequence (nt 781–792), similar to the OH of mammals, maps just a few nucleotides downstream from the putative CSB1 (nt 803–828). However, it is difficult to identify sequences corresponding to mammalian CSB2 and CSB3. The secondary structures of CSB1 and the putative goose hairpin in Crossoptilon are consistent with those of G. gallus (Additional file 7). The bidirectional LSP/HSP promoters (nt 982–1003)  are almost perfectly conserved among the Crossoptilon species. A stable hairpin (nt 1004–1017) that is rich in poly (T) and poly (A) strings is located immediately upstream of the promoters.
Phylogeny and divergence time of the Crossoptilon species
The final combined PCG dataset has 11376 characters after alignment. For this dataset, parsimony analysis shows a length of 27412 steps, with CI =0.313, RI = 0.499. The topologies between ML and BI trees of PCG dataset were consistent (Figure 1). For the mitogenome dataset, there was no difference in topology between the ML and BI trees of Galliformes (Additional file 8); however, both trees differed from the MP tree (Additional file 8). According to the phylogenetic results, the monophyly of Crossoptilon was strongly supported in MP, ML and BI analyses. Within Phasianidae, the topology ((Crossoptilon + Lophura) + (Phasianus + Chrysolophus)) was formed in most trees. The sister-group relationship between Crossoptilon and Lophura was supported with bootstrap values 69 and 100 in MP and ML trees, and posterior probabilities 1.00 in BI tree of PCG dataset.
The genetic distances of the mitogenomes (Additional file 9) are identical (0.026) for C. auritum and C. harmani, C. auritum and C. crossoptilon, C. harmani and C. mantchuricum, and C. mantchuricum and C. crossoptilon. However, the genetic distances between C. harmani and C. crossoptilon, and C. auritum and C. mantchuricum are only 0.001 and 0.002, respectively. The genetic distances based on four single genes show that the genetic distances of C. harmani-C. crossoptilon and C. auritum-C. mantchuricum are much smaller. This value is similar to the intraspecific level and far less than the interspecific level (Additional file 9).
According to the estimated timescale obtained from the phylogenetic tree containing consistent topology with PCG-ML/BI trees, the C. harmani-C. crossoptilon and C. mantchuricum-C. auritum splits occurred at 3.21 Mya (95% highest posterior-probability density (HPD) = 2.36-4.23 Mya) (Figure 2). The divergence time for C. harmani and C. crossoptilon is approximately 0.11 Mya (95% HPD = 0.05-0.18 Mya) and for C. mantchuricum and C. auritum, 0.18 Mya (95% HPD = 0.1-0.29 Mya) (Figure 2).
Nonsynonymous and synonymous substitution
The analysis of variable sites in each protein-coding gene showed that the group C.har-C.cro and C.man-C.aur contained less variable sites, while other four groups included more (Additional file 10). The percentages of variable sites were higher in nad3 and nad6 genes among the groups C.har-C.man, C.man-C.cro, C.har-C.aur and C.cro-C.aur, while the percentages were lower in atp8 and atp6 genes. The Ka/Ks value of atp8 in the five groups is far greater than 1 except for C.har-C.cro, and atp6 and cytb in C.har-C.cro (Table 2), which show a strong positive selection. The Ka/Ks values of the other genes analysed in Crossoptilon species were less than 1, which shows a purifying selection (Table 2). Furthermore, the P-value (Fisher exact test) is much less than 0.001, except for atp8 in C.har-C.man, C.man-C.cro, C.har-C.aur and C.cro-C.aur (0.558274 or 0.289204, respectively, which are obviously higher than 0.05) and nad3 in C.har-C.cro (0.154), indicating that the difference is significant.
Both atp6 and cytb genes have one different base in the C.har-C.cro group, i.e., nt 245 in atp6 gene (T base or Phe in C. harmani, while C base or Ser in C. crossoptilon) and nt 709 in cytb gene (C base or Leu in C. harmani, while T base or Phe in C. crossoptilon). The atp8 gene in C.man-C.aur group also contains one different base (nt 124), G base in C. mantchuricum, while A base in C. auritum, and conrespondingly Val in C. mantchuricum, while Met in C. auritum. Further study showed that four amino acid sites in atp8 (T at 11 position, I at 12 position, S at 39 position and V at 42 position) in the genus Crossoptilon were positively selected sites, based on Bayes Empirical Bayes (BEB) analysis. This four different amino acid sites (T, I, S and V) corresponded to nucleotide sequences of ACT, ATC, AGC and GTA in C. harmani and C. crossoptilon, ATT, ACC, AAC and GTA in C. mantchuricum, ATT, ACC, AAC and ATA in C. auritum. The further secondary structure analysis of the amino acid sequence of atp8 showed that T at 11 position and I at 12 position located at helix (IN- > OUT) domain in transmembrane (TM) protein, while S at 39 position and V at 42 position located at outside domain.
PBS analyses were performed to better understand the contributions of different parts of the mitogenomes to the genome phylogeny based on the mitogenome-ML tree. The relationships among the PBS value, length, singleton sites (S), parsimony informative sites (Pi), variable sites (V) and conserved sites (C) in different partitions are also studied. Ranking individual protein-coding genes by their respective contribution to the total PBS values shows that some genes, such as nad5 and nad4 provide a higher contribution compared to other markers, while atp8 and nad4L contribute less. But analysis also shows that PBS values are roughly correlated with gene lengths. The variable sites are closely related to parsimony informative sites, which follow the same trend. The third codon has a high PBS value but with less conserved sites, while the tRNA gene cluster IQM has a lower PBS value with more conserved sites compared to other partitions.
The third codon has a high PBS value and therefore contributes most to the mitogenome-ML tree. The NADH genes have secondary high PBS, which might be the effect of the third codon in the relatively long sequence. To further understand the contribution of the 3rd codon, we reconstructed the phylogenetic tree using this single partition.
In the phylogenetic trees based on the third codon (Additional file 8), the structure of the ML tree is similar to the BI with the exception of Phasianus. The positions of Pucrasia macrolopha, Bonasa bonasia, Meleagris gallopavo and Polyplectron bicalcaratum in the 3rd-MP tree are different compared to the BI tree. The 3rd-BI tree is identical to the mitogenome-ML with relative high support except for the position of Polyplectron bicalcaratum.
Phylogenetic relationship and divergence time of the Crossoptilon species
In this study, the phylogenetic trees based on different datasets support the monophyly of Crossoptilon and the close genetic relationships of C. auritum and C. mantchuricum [4,12,54,55], C. harmani and C. crossoptilon [1,4,8,14,55]. Our analysis of the gene sequences of Crossoptilon based on different datasets confirms that Crossoptilon is the sister group to Lophura [55,56], indicating their closer relationship.
The genetic distance is rather small among some Crossoptilon species, for instance, between C. harmani and C. crossoptilon, C. auritum and C. mantchuricum. Morphological, ecological and behavioural studies have revealed the high similarity between C. harmani and C. crossoptilon, but with different plumage coloration . Their distributions are overlapping, and C. harmani routinely hybridises with C. c. drouynii , which indicates their reproductive compatibility. Furthermore, the distribution of C. harmani is limited between the Himalayan and Nyenchen Tanglha ranges, which provided the geographical isolation required for the formation of a subspecies or species. Therefore, the taxonomic status and its relationship to C. crossoptilon are much questionable, which need to be re-evaluated by multiple markers and population genetics in the further studies.
Based on their significant differences in many aspects such as morphology and behaviour, C. auritum and C. mantchuricum are commonly considered to be independent species. However, the genetic distance of 0.002 indicates that they have a much closer genetic relationship. Consistent with this observation, Tsam et al. (2003)  observed that the interspecific differentiation between C. mantchuricum and C. auritum (0.18%) is less than the degree of differentiation among the C. crossoptilon subspecies.
Some Phasianidae species were found in the Pliocene epoch, and several modern taxa had already appeared in the Quaternary Pleistocene. The earliest fossils of Crossoptilon were found in the Cenozoic Zhoukoudian strata of Beijing, and in the Pleistocene strata at Yanjinggou in Wanxian of Sichuan Prov. [57,58]. The evolution of Crossoptilon mainly occurred corresponding to the Tibetan Plateau intensive uplift during the Tertiary Pliocene and the alternation of the glacial and interglacial periods in the Quaternary, which created profound and complex changes in the geographical and ecological environments [59,60] and was thought to have greatly affected its topography and avian species diversification [61,62]. Based on our results, the spliting time between C. harmani-C. crossoptilon and C. mantchuricum-C. auritum is consistent with the uplift during the Tertiary Pliocene and the fossil record in the Sichuan Wanxian Yanjinggou Pleistocene strata, which is approximate to the divergence times estimated by Jiang et al. (2014)  (3.78 Mya, 95% HPD = 1.17-6.56 Mya). Our study shows that C. mantchuricum and C. auritum are relatively ancient groups, while C. harmani and C. crossoptilon diverged later; these observations are not consistent with the Crossoptilon evolution hypothesis proposed by Lu et al. (1998) .
Lu et al. (1998)  suggested that the ancestor of C. crossoptilon and C. harmani originated in Sichuan Prov., Yunnan Prov., and the Tibetan border area. They noted that the clade containing the C. harmani plesiomorphy migrated to the hinterland plateau and subsequently differentiated into the ancestor of C. harmani and C. auritum. The ancestor of C. auritum dispersed in the northern plateau, and differentiated into C. auritum and C. mantchuricum. Based on the divergence time and distribution characteristics of Crossoptilon species, our results support the ancestor of Crossoptilon was initially distributed in Sichuan Prov., Yunnan Prov., and the Tibetan border area. The Crossoptilon first diverged into the ancestor of C. harmani-C. crossoptilon and the ancestor of C. auritum-C. mantchuricum. Considering the geographical distribution of Crossoptilon [1,64-66], we propose that the ancestor of C. harmani-C. crossoptilon might migrated to platform of the plateau, and then splitted into two lineages during the last uplift of the Tibetan Plateau , one adapted the high altitude environment to form C. harmani in east Tibet, and the other evoluted into C. crossoptilon, while the ancestor of C. auritum-C. mantchuricum dispersed to the northern plateau and further differentiated nearly the same period with C. harmani-C. crossoptilon.
Selective pressure on protein-coding genes in Crossoptilon species
The Ka/Ks values in most protein-coding genes are less than 1 (Ka is lower than Ks) with P-values less than 0.001, which indicates that these genes in Crossoptilon are under purifying selection. The Ka/Ks values of each gene among C.har-C.man, C.man-C.cro, C.har-C.aur and C.cro-C.aur are similar. The Ka/Ks value of the atp8 gene is 50 in five groups except for C.har-C.cro, which reveals the strong positive selection. However, only the P-value for C.man-C.aur is far less than 0.01 (P = 0), which indicates significant differences. In contrast, the Ka/Ks values of the atp6 and cytb genes in C.har-C.cro are 50 with P-values less than 0.001, which also indicates strong positive selection. Previous studies have shown that atp6 gene is highly variable in the mitogenomes of humans living in extremely cold areas [68-71], the atp6 gene faces strong selection pressure with increasing altitude , and the variation of atp6 gene in Artemia tibetiana is the result of adaptation to the cold and hypoxic environment of the plateau . The variation in atp6 may reflect the adaptation of C. harmani to the plateau environment; cytb might also play an important role in this process. The atp8 gene in C. mantchuricum and C. auritum may also have experienced a strong selection in adaption to the plateau environment.
However, the further study showed that only four amino acid sites of atp8 in the Crossoptilon are positively selected sites, and accordingly this gene may have experienced strong selective pressure for plateau adaptation. The corresponding TM-helix (IN- > OUT) or outside position might be important for the Crossoptilon adaptation. After discarding groups (C.har-C.cro and C.man-C.aur) and genes (cox1, atp8 and nad4l) with obviously different Ka/Ks values to others, the average value (0.0431) is similar with the strongly locomotive group (0.04) in the study of Shen et al. (2009) , which indicates the Crossoptilon accumulates fewer nonsynonymous mutations and is corresponding to their strong running ability.
In summary, the mitogenomes of the Crossoptilon species contain the same gene arrangements and similar compositions including base contents and secondary structures. According to the phylogenetic trees, C. auritum with C. mantchuricum, and C. harmani with C. crossoptilon form two pairs of sister groups. Crossoptilon have a closer relationship with Lophura. Based on the genetic distances, C. harmani is more closely related with C. crossoptilon, and is the most recent diverged descendent of Crossoptilon as a result of the plateau adaptation. According to the molecular dating results, the divergence time between C. harmani-C. crossoptilon and C. mantchuricum-C. auritum is consistent with the uplift of the Tibetan plateau and the subsequently climate change during the Tertiary Pliocene. The Ka/Ks analysis showed that atp8 gene in the Crossoptilon species may have experienced strong selection for the plateau environment adaptation.
- atp6 and atp8 :
ATP synthase subunits 6 and 8
- cytb :
- cox1-3 :
Cytochrome c oxidase subunit 1 to 3
NADH dehydrogenase, subunit 1–6
- nad4L :
NADH dehydrogenase, subunit 4L
- rrnS and rrnL :
Small and large ribosomal RNA (rRNA) subunits
Lu X, Zheng GM, Gu BY. A preliminary investigation on taxonomy, distribution and evolutionary relationship of the eared pheasants, Crossoptilon. Acta Zool Sin. 1998;44(2):131–7.
Cheng TH. A complete checklist of species and subspecies of the Chinese birds. Beijing: Science Press; 1994. p. 31.
Lu TC. The rare and endangered wild chicken in Chain. Fuzhou: Fujian Science and Technology Press; 1991.
Cheng TH, Tan YK, Lu TC, Tang CG, Bao GJ, Li FL. Fauna of China. Aves, Vol. IV. Galliformes. Beijing: Science Press; 1978. p. 106–7.
Johnsgard PA. The Pheasants of the World. New York: Oxford University Press; 1986.
Shi XD, Zhang ZW, Liu LY. Karyo types and G-banding patterns of three eared-pheasant (Crossoptilon) species. Acta Zool Sin. 2001;47(3):280–4.
Gan YL, Lu TC, Liu RS, He FQ, Lu CL, Gan YL, et al. Observation on scanning electron microscope of eggshell of C. mantchuricumis, C. crossoptilon and C. auritum endemic pheasants in China. Acta Zool Sin. 1992;38(2):124–7.
Zheng GM, Zhang W, Zhao XR. A comparative research on the hind limb muscles of eared pheasants, Crossoptilon. In: Gao W, editor. The Study of Birds in China. Beijing: Science Press; 1991. p. 14–6.
Lu TC, Liu RS, He FQ, Lu CL, Li GY. Ecology and systematic relationship of three species of the genus Crossoptilon. Sichuan J Zool. 1989;8:21–3.
Liu RS, Guo YJ, Li FL, Hou LH. Study on the relationship among three species of the genus Crossoptilon by electraofocusing technique. Acta Zool Sin. 1985;31(2):206–13.
Ludlow F, Kinnear NB. The birds of south-eastern Tibet. Ibis. 1944;86:348–89.
Tsam CDM, Rao G, Ji JG, Suo LCR, Wan QH, Fang SG. Taxonomic status of Crossoptilon harmani and a phylogenetic study of the genus crossoptilon. Acta Zool Sin. 2003;28(2):173–9.
Rothschild L. On the avifauna of Yunnan, with critical notes. Novit Zool. 1926;33:189–343.
Delacour J. The Pheasant of the World. 2nd ed. London: World Pheasant Association and Spur Publications; 1977.
Wu AP, Ding W, Zhang ZW, Zhan XJ. Phylogenetic relationship of the avian genus Crossoptilon. Acta Zool Sin. 2005;51(5):898–902.
Yang C, Lei FM, Huang Y. Sequencing and Analysis of the Complete Mitochondrial Genome of Pseudopodoces humilis (Aves, Paridae). Zool Res. 2010;31(4):333–44.
Sorenson MD. Avian mtDNA primers. 2003. Available from: http://people.bu.edu/msoren/primers.html.
Staden R, Beal KF, Bonfield JK. The Staden package, 1998. Methods Mol Biol. 2000;132:115–30.
Lowe TM, Eddy SR. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997;25:955–64.
Cannone JJ, Subramanian S, Schnare MN, Collett JR, D'Souza LM, Du Y, et al. The Comparative RNA Web (CRW) site: an online database of comparative sequence and structure information for ribosomal, intron, and other RNAs. BMC Bioinformatics. 2002;3(2):1–31.
Burk A, Douzery EJP, Springer MS. The secondary structure of mammalian mitochondrial 16S rRNA molecules: refinements based on a comparative phylogenetic approach. J Mamm Evol. 2002;9(3):225–52.
Tamura K, Dudley J, Nei M, Kumar S. MEGA4: molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol. 2007;24(8):1596–9.
Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG. The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997;24:4876–82.
Swofford DL. PAUP*. Phylogenetic Analysis Using Parsimony (*and Other Methods), Version 4. Sunderland, MA: Sinauer Associates; 2003.
Stamatakis A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinform. 2006;22:2688–90.
Ronquist F, Huelsenbeck JP. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003;19:1572–4.
Nylander JAA. MrModeltest v2. Program distributed by the author. Evolutionary Biology Centre, Uppsala University;2004.
Sorenson MD, Franzosa EA. TreeRot, version 3. Boston: Boston University, MA; 2007.
Drummond AJ, Rambaut A. BEAST: bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007;7:214.
Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29(8):1969–73.
Drummond AJ, Ho SYW, Phillips MJ, Rambaut A. Relaxed phylogenetics and dating with confidence. PLoS Biol. 2006;4:e88.
Dyke GJ, Gulas BE, Crowe TM. Suprageneric relationships of galliform birds (Aves, Galliformes): a cladistic analysis of morphological characters. Zool J Linnean Soc. 2003;137:227–44.
Crowe TM, Bowie RCK, Bloomer P, Mandiwana TG, Hedderson TAJ, Randi E, et al. Phylogenetics, biogeography and classification of, and character evolution in, gamebirds (Aves: Galliformes): effects of character exclusion, data partitioning and missing data. Cladistics. 2006;22:495–532.
Tordoff HB, Macdonald JR. A new bird (family Cracidae) from the early Oligocene of South Dakota. Auk. 1957;74:174–84.
Brodkorb P. Catalogue of fossil birds, part 2 (Anseriformes through Galliformes). Bull Florida State Mus Biol Sci. 1964;8:195–335.
Mourer-Chauviré C. The Galliformes (Aves) of the Phosphorites du Quercy (France): systematics and biogeography. Natur Hist Mus Los Angeles County Sci Ser. 1992;36:67–95.
Zhang Z, Li J, Zhao XQ, Wang J, Wong GK, Yu J. KaKs_Calculator: calculating Ka and Ks through model selection and model averaging. Genomics Proteomics Bioinforma. 2006;4:259–63.
Nei M, Kumar S. Molecular Evolution and Phylogenetics. New York: Oxford University Press; 2000.
Yang Z. Computational Molecular Evolution. New York: Oxford University Press; 2006.
Xu B, Yang Z. PAMLX: a graphical user interface for PAML. Mol Biol Evol. 2013;30:2723–4.
Bernsel A, Viklund H, Hennerdal A, Elofsson A. TOPCONS: consensus prediction of membrane protein topology. Nucleic Acids Res. 2009;37(W):W465–468.
Wolstenholme DR. Animal mitochondrial DNA: structure and evolution. Int Rev Cytol. 1992;141:173–216.
Quinn TW. The genetic legacy of Mother Goose–phylogeographic patterns of lesser snow goose Chen caerulescens caerulescens maternal lineages. Mol Ecol. 1992;1:105–17.
Fumihito A, Miyake T, Sumi S-I, Takada M, Ohno S, Kondo N. One subspecies of the red junglefowl (Gallus gallus gallus) suffices as the matriarchic ancestor of all domestic breeds. Proc Natl Acad Sci U S A. 1994;91:12505–9.
Fumihito A, Miyake T, Takada M, Ohno S, Kondo N. The genetic link between the Chinese bamboo partridge (Bambusicola thoracica) and the chicken and junglefowls of the genus Gallus. Proc Natl Acad Sci U S A. 1995;92:11053–6.
Ramirez V, Savoie P, Morais R. Molecular characterization and evolution of a duck mitochondrial genome. J Mol Evol. 1993;37:296–310.
Douzery E, Randi E. The mitochondrial control region of Cervidae: evolutionary patterns and phylogenetic contents. Mol Biol Evol. 1997;14:1154–66.
Dufresne C, Mignotte F, Guéride M. The presence of tandem repeats and the initiation of replication in rabbit mitochondrial DNA. Eur J Biochem. 1996;235:593–600.
Gemmell NJ, Western PS, Watson JM, Marshall-Graves JA. Evolution of the mammalian mitochondrial control region—comparisons of control region sequences between monotreme and therian mammals. Mol Biol Evol. 1996;13:798–808.
Stewart DT, Baker AJ. Patterns of sequence variation in the mitochondrial D-loop region of shrews. Mol Biol Evol. 1994;11:9–21.
Fumagalli L, Taberlet P, Favre L, Hausser J. Origin and evolution of homologous repeated sequences in the mitochondrial DNA control region of shrews. Mol Biol Evol. 1996;13:31–46.
Sbisà E, Tanzariello F, Reyes A, Pesole G, Saccone C. Mammalian mitochondrial D-loop region structural analysis: identification of new conserved sequences and their functional and evolutionary implications. Gene. 1997;205:125–40.
L’Abbé DL, Duhaime JF, Lang BF, Morais R. The transcription of DNA in chicken mitochondria initiates from one major bidirectional promoter. J Biol Chem. 1991;266:10844–50.
Tang CZ. The analysis of system classification and geographical distribution of Crossoptilon. Acta Zool Sin. 1998;23:86–92.
Wang N, Kimball RT, Braun EL, Liang B, Zhang ZW. Assessing phylogenetic relationships among Galliformes: a multigene phylogeny with expanded taxon sampling in Phasianidae. PLoS One. 2013;8(5):e64312.
Shen YY, Dai K, Cao X, Murphy RW, Shen XJ, Zhang YP. The Updated Phylogenies of the Phasianidae Based on Combined Data of Nuclear and Mitochondrial DNA. PLoS One. 2014;9(4):e95786.
Wetmore A. Fossil birds from Mongolia and China. Am Aus Novit. 1934;711:1–16.
Hou LH. Avian fossils of pleistocene from Zhoukoudian, China. Vertebrata Pal Asiatica. 1982;20:366–8.
Yang YC, Li BY, Yin ZS, Zhang QS, Wang FB, Jing K, et al. Geomorphology of Xizang (Tibet). Beijing: Science Press; 1983. p. 1–238.
Li BY, Wang FB. Basic characteristics of landforms in the northwest Yunnan and southwest Sichuan area. In: Chinese Academy of Sciences, editor. The Comprehensive Scientific Expedition to the Qinghai-Xizang Plateau. Studies in Qinghai-Xizang (Tibet) Plateau Special Issue of Hengduan Mountains Scientific Expedition (II). Beijing: Beijing Science & Technology Press; 1986. p. 175–83.
Li JJ, Shi YF, Li BY. Uplift of the Qinghai–Xizang (Tibet) Plateau and Global Change. Lanzhou: Lanzhou University Press; 1995.
Lei FM, Qu YH, Song G. Species diversification and phylogeographical patterns of birds in response to the uplift of the Qinghai-Tibet Plateau and Quaternary glaciations. Curr Zool. 2014;60(2):149–61.
Jiang LC, Wang GC, Peng R, Peng QK, Zou FD. Phylogenetic and molecular dating analysis of Taiwan Blue Pheasant (Lophura swinhoii). Gene. 2014;539(1):21–9.
Wang FL, Chen JM, Lai RX. Studies on the ancient and modern geographical distribution of Brown-eared Pheasants. J Shanxi Univ. 1985;3:86–92.
Zhang CA, Ding CQ. The distribution pattern of the Galliformes in China. Acta Zool Sin. 2008;33(2):317–23.
Wu MX, Wu JG, Kuang MS, Heng T. Relationship between Geographic Distribution of Endemic Birds and Climatic Factors in China. Res Env Sci. 2011;24(4):409–20.
Li JJ, Fang XM. Uplift of the Tibetan Plateau and environmental changes. Chinese Sci Bull. 1999;44:2117–24.
Elson JL, Turnbull DM, Howell N. Comparative genomics and the evolution of human mitochondrial DNA: assessing the effects of selection. Am J Hum Genet. 2004;74:229–38.
Mishmar D, Ruiz-Pesini E, Golik P, Macaulay V, Clark AG, Hosseini S, et al. Natural selection shaped regional mtDNA variation in humans. Proc Natl Acad Sci U S A. 2003;100:171–6.
Coskun PE, Ruiz-Pesini E, Wallace DC. Control region mtDNA variants: longevity, climatic adaptation, and a forensic conundrum. Proc Natl Acad Sci U S A. 2003;100:2174–6.
Bhopal RS, Rafnsson SB. Could mitochondrial efficiency explain the susceptibility to adiposity, metabolic syndrome, diabetes and cardiovascular diseases in South Asian populations? Int J Epidemiol. 2009;38:1072–81.
Gu ML, Wang YJ, Shi L, Zhang YB, Chu JY. Comparison on mitochondrial atp6, atp8 and cytb genes between Chinese Tibetans in three different zones: detecting the signature of natural selection on mitochondrial genome. Hereditas (Beijing). 2009;31:147–52.
Zhang HX, Luo QB, Sun J, Liu F, Wu G, Yu J, et al. Mitochondrial genome sequences of Artemia tibetiana and A. urmiana: assessing molecular changes for high plateau adaptation. Sci China Life Sci. 2013;56:440–52.
Shen YY, Shi P, Sun YB, Zhang YP. Relaxation of selective constraints on avian mitochondrial DNA following the degeneration of flight ability. Genome Res. 2009;19(10):1760–5.
This research was supported by State Key Program of NSFC (31330073), the Strategic Priority Research Program of the Chinese Academy of Sciences (Grant No. XDB13020300, XDA05080703), the Ministry of Science and Technology of the People’s Republic of China (MOST Grant No. 2011FY120200-3) and the Key Laboratory of the Zoological Systematics and Evolution of the Chinese Academy of Sciences (No. O529YX5105).
The authors declare that they have no competing interests.
YH, FML and XJL conceived and designed the experimental plan. XJL carried out the experimental design, and drafted the manuscript, tables, and figures. YH and FML analyzed and interpreted the sequence data, and extensively revised the manuscript. All authors read and approved the final manuscript.
PCR primers used in this study. The mitogenome of C. harmani was amplified by seven parts, that is L1263x-H2891x, L2260x-H6681x, L5758x-H7122x, L6615x-H10884x, L10635x-H13563x, L13040x-H16064x and L15725x-H1530x. The C. mantchuricum mitogenome was amplified by 12 fragments, that is L1263x-H2891x, L2260x-H4644x, L3803x-H6681x, L5758x-H7122x, L6615x-H8121x, L7525x-H8628x, L8386x-H10884x, L10635x-H12344x, L11458x-H13563x, L13040x-H15646b, L14080x-H16064x and L15413b-H1530x. The C. crossoptilon mitogenome was amplified by 13 fragments, that is L1263x-H2891x, L2260x-H4644x, L3803x-H6681x, L5758x-H7122x, L6615x-H8121x, L7525x-H8628x, L8386x-H10884x, L10635x-H12344x, L11458x-H13563x, L12156b-H14127b, L13525b-H15049x, L14080x-H16064x and L15725x-H1530x. The gaps were completed by other adjacent primers.
Species studied in the phylogenetic analysis.
The single genes ( nad2 , CR, cytb and rrnS ) used to calculate genetic distances among Crossoptilon species.
The predicted secondary structure of rrnS in C. harmani . Nucleotide differences across all other three species are labled (C. mantchuricum: C. man, C. crossoptilon: C. cro, C. auritum: C. aur).
The predicted secondary structure of rrnL in C. harmani . Nucleotide differences across all other three species are labled (C. mantchuricum: C. man, C. crossoptilon: C. cro, C. auritum: C. aur).
The secondary structure of all tRNAs in C. harmani . Nucleotide differences across all other three species are labled (C. mantchuricum: C. man, C. crossoptilon: C. cro, C. auritum: C. aur).
The putative secondary structure of the control region in Crossoptilon . Note: A: nt 1–190 in C. crossoptilon and C. harmani; B: nt 1–190 in C. auritum and C. mantchuricum.
The MP, ML and BI trees of different datasets. Less than 50% bootstrap values were omitted. Notes: (1) results based on PCG dataset; (2) results based on mitogenome dataset; (3) results based on the third codons dataset (3rd). (A) MP tree; (B) ML tree; (C) BI tree.
Interspecific and intraspecific genetic distances of the Crossoptilon species based on mitogenome dataset and single genes. Note: C. harmani: C.har, C. mantchuricum: C. man, C. crossoptilon: C. cro, C. auritum: C. aur.
The variable sites of protein-coding genes among six Crossoptilon groups. Note: C. harmani: C.har, C. mantchuricum: C. man, C. crossoptilon: C. cro, C. auritum: C. aur.
About this article
Cite this article
Li, X., Huang, Y. & Lei, F. Comparative mitochondrial genomics and phylogenetic relationships of the Crossoptilon species (Phasianidae, Galliformes). BMC Genomics 16, 42 (2015). https://doi.org/10.1186/s12864-015-1234-9
- Mitochondrial genome
- Divergence time