Research Article
Print
Research Article
Multi-locus analysis supports the taxonomic validity of Arborophila gingica guangxiensis Fang Zhou & Aiwu Jiang, 2008
expand article infoDe Chen, Qiong Liu, Jiang Chang§, Aiwu Jiang|, Fang Zhou, Yanyun Zhang, Zhengwang Zhang
‡ Beijing Normal University, Beijing, China
§ Chinese Research Academy of Environmental Sciences, Beijing, China
| College of Forestry, Guangxi University, Nanning, China
¶ Guangxi University, Nanning, China
Open Access

Abstract

The taxonomic status of subspecies has long been debated, especially in conservation biology. Some proposed subspecies must be evolutionarily distinct to be considered conservation units. White-necklaced Partridge (Arborophila gingica) comprises two subspecies, A. g. gingica and A. g. guangxiensis. A. g. guangxiensis, restricted to three isolated small areas in Guangxi, China, with limited population sizes, is a newly discovered subspecies based on recently identified geographic and phenotypic differences between A. g. gingica; however, evidence is lacking that can effectively identify whether the subspecies is evolutionarily distinct. Here, three mitochondrial DNA segments and four nuclear introns were used to test whether the two subspecies are reciprocally monophyletic, which has been proposed as an objective method to evaluate evolutionary distinctiveness. The results indicate that the two subspecies are genetically divergent and form reciprocal monophyletic groups. Therefore, this study further supports the taxonomic validity and distinctiveness of A. g. guangxiensis and suggests that this subspecies be considered as a conservation unit.

Keywords

Conservation unit, evolutionary distinctiveness, hill partridge, mitochondrial DNA, nuclear introns, monophyletic groups

Introduction

The taxonomic status of subspecies has long been debated (Wilson and Brown 1953; Mayr 1982), especially in resource-limited conservation biology (Zink 2004; Phillimore and Owens 2006). Zink (2004) proposed that subspecies must be evolutionarily distinct to be considered conservation units. However, a recent global analysis showed that only 36% of traditional avian subspecies can be distinguishable by mitochondrial DNA (mtDNA) (Phillimore and Owens 2006). Therefore, assessing the validity of subspecies before proposing conservation efforts may be a practical way to effectively protect biodiversity with limited resources.

The near threatened (NT) White-necklaced Partridge (Arborophila gingica) (BirdLife International 2012), also known as the Collared or Rickett’s hill partridge, is a small partridge endemic to the southern Chinese forests of Hunan, Jiangxi, Zhejiang, Fujian, Guangdong and Guangxi (Fig. 1, Cheng 1978; Zheng 2011). Although the distribution area appears extensive, populations are severely fragmented and continuously declining (Zhou 1996; BirdLife International 2012), except for the Fujian population (He et al. 2007). This species was believed to be monotypic (Johnsgard 1988; Madge and McGowan 2002) until Zhou and Jiang (2008) discovered that the populations in north and central Guangxi differed from other southeastern populations in the coloration of the forehead, which is chestnut instead of white. Zhou and Jiang (2008) described these populations as a new subspecies, A. g. guangxiensis, a proposal subsequently accepted by several avian checklists (Zheng 2011; Dickinson and Remsen 2013; del Hoyo and Collar 2014; Gill and Donsker 2015). However, except for its diagnostic forehead coloration, no other plumage differences are known, and body weight, body length, wing, culmen, tarsus and tail show no obvious differences between the two subspecies (Zhou and Jiang 2008). A. g. guangxiensis is only found in three isolated small areas of north (Jiuwanshan Mountain and Sijianshan Mountain) and central (Damingshan Mountain) Guangxi (Fig. 1, Zhou and Jiang 2008). The estimated population is about 600 to 1000 individuals, thus undoubtedly A. g. guangxiensis suffers more threats than A. g. gingica (Taylor 2011). However, the subspecies status of A. g. guangxiensis is solely dependent on geographic and phenotypic differences, lacking molecular evidence that can effectively identify whether this isolated subspecies is evolutionarily distinct and can be considered a conservation unit (Zink 2004). Indeed, some operational criteria for subspecies recognition require that subspecies are both phenotypically distinct and correlate with evolutionary independence according to population genetic structure (Braby et al. 2012).

Figure 1.

Map of southeast China showing the distribution area of A. gingica. The purple area represents the distribution of the nominate subspecies A. g. gingica according to del Hoyo and Collar (2014), with the dark blue dot indicating the sampling site in Wuyishan Mountain. The orange dots represent the three isolated populations of A. g. guangxiensis (Zhou and Jiang 2008), with the sampling site in Jiuwanshan Mountain.

Recently, molecular systematics has become one of the most vigorous disciplines to assist in avian taxonomy (Fjeldså 2013). Specifically, mtDNA has been extensively used at various taxonomic levels (Zink and Barrowclough 2008). Zink (2004) advocated that subspecies should be reciprocally monophyletic in mtDNA gene trees to document the evolutionary distinctiveness of subspecies. However, studies that are based solely on mtDNA have been debated because population differentiation relies on the accumulated signals from many genes and mtDNA only represents a single locus (Rubinoff and Holland 2005; Edwards and Bensch 2009). Therefore, a reasonable strategy for phylogenetic analysis is to combine mtDNA with nuclear DNA (nuDNA) sequences (Corl and Ellegren 2013).

Here, three mtDNA segments and four nuclear introns of White-necklaced Partridge were combined to conduct a series of phylogenetic analyses and test whether A. g. guangxiensis and A. g. gingica form reciprocally monophyletic groups. Furthermore, times of divergence within A. gingica, between A. gingica and its closest relative, were investigated, and attempts to identify possible drivers of the diversification process were made.

Methods

Sampling, DNA extraction, PCR and sequencing

Seven individuals of A. g. guangxiensis were sampled from Jiuwanshan National Nature Reserve, Guangxi, and three individuals of A. g. gingica from Wuyishan National Nature Reserve, Jiangxi (Fig. 1). Previous studies indicated that the sister species of A. gingica was A. rufogularis (Wang et al. 2013); therefore, we used one individual of A. rufogularis from Tongbiguan National Nature Reserve in Yunnan as an outgroup. All samples were taken from live birds (blood or feather). Permissions for blood or feather sampling were granted by the regional forestry departments. Total DNA was extracted using a TIANamp Blood Genomic DNA Extraction Kit (TIANGEN BIOTECH CO, BEIJING, CHINA).

We amplified three mtDNA segments, cytochrome oxidase subunit 1 (COI), cytochrome b (CYTB) and NADH dehydrogenase subunit 2 (ND2); and four nuclear introns, aldolase b intron 6 (ALDOB), fibrinogen intron 5 (FGB), glyceraldehyde 3-phosphate dehydrogenase intron 11 (G3PDH) and ovomucoid intron G (OVOG) using the primers listed in Suppl. material 1. Both strands of each PCR product were sequenced by BGI-BEIJING. The sequences were visually proofread to the original chromatograms and were also checked against published DNA sequences. Each sequence was then assembled using MEGA v6 (Tamura et al. 2013). Then, we aligned the sequences using MUSCLE (Edgar 2004) implemented in MEGA v6 (Tamura et al. 2013) to obtain seven partitions, all sequences obtained from this study were submitted to GenBank (KU057820KU057877). Each nuclear partition was then phased (Stephens and Donnelly 2003) in DNASP v5.10 (Librado and Rozas 2009) to resolve the haplotypes of diploid nuclear sequences. Finally, we assembled the seven partitions into a complete matrix, an mtDNA matrix and a nuDNA matrix.

Phylogenetic analysis

The best-fitting nucleotide substitution model for each partition was selected using the Akaike Information Criterion with JMODELTEST v2.1.7 (Darriba et al. 2012). The mean genetic distances between and within subspecies were calculated in MEGA v6 (Tamura et al. 2013) using the Kimura two-parameter (K2P) model (with A. rufogularis removed); and standard error estimates were obtained by a bootstrap procedure (1000 replicates). Partitioned maximum likelihood (ML) analyses were conducted in GARLI v2.0 (Bazinet et al. 2014) using the best-fitting nucleotide substitution model for each partition. The subtree pruning and regrafting tree-searching method was used, and bootstrap values (BS) were calculated with 1,000 replicates. Partitioned Bayesian Inference (BI) was performed in BEAST v1.8.0 (Drummond and Rambaut 2007) with the best-fitting nucleotide substitution model for each partition (similar to Divergence time estimates, see below for details).

Divergence time estimates

First, we performed molecular clock tests in MEGA v6 (Tamura et al. 2013). The results showed that each partition was clock-like. Therefore we used the strict clock model for each partition. It is believed that a species tree analysis using combined mtDNA, Z-linked (ALDOB) and autosomal (FGB, G3PDH and OVOG) loci can substantially improve the resolution of the tree (Corl and Ellegren 2013). Therefore, we performed a species tree analysis using the complete matrix in *BEAST (Heled and Drummond 2010) implemented in BEAST v1.8.0 (Drummond and Rambaut 2007), with a fixed molecular rate of 2.38% for CYTB (average molecular rate for Galliform birds, Weir and Schluter (2008)) to estimate the molecular rates of the other loci. The ESS value was verified to be greater than 200 in TRACER v1.5 (Rambaut and Drummond 2009) to confirm that the chains had reached apparent stationarity. The final analysis was run for 100 million generations with trees sampled every 1,000 generations. TreeAnnotator v1.8.0 was then used to discard the first 20% of trees and to generate the consensus tree with Bayesian posterior probability.

Results

The complete matrix was 4750 base pairs (bp) in length, including 2861 bp of mtDNA sequence data, and 1889 bp of nuclear intron sequence data. Exclude outgroup, there were 18 variable and 13 informative sites in mtDNA, and 24 variable and 19 informative sites in nuDNA (after phasing). The genetic distance between the two subspecies was higher in mtDNA (0.0038) than in nuDNA (0.0028), and in nuclear introns the genetic distance within subspecies partially overlapped with that between subspecies (Table 1) due to some shared haplotypes (data not shown).

Phylogenetic analyses of the complete matrix and mtDNA matrix showed that A. g. guangxiensis and A. g. gingica formed monophyletic groups, with relatively high support (Fig. 2). However, analyses of the nuDNA matrix and separate analyses of each of the nuclear introns failed to recover the monophyletic relationships between A. g. guangxiensis and A. g. gingica (Suppl. material 2), and the support values are extremely low (data not shown).

Divergence time estimates from the species tree showed that the two subspecies A. g. guangxiensis and A. g. gingica diverged approximately 0.11 (0.05–0.19) mya (million years ago), whereas the divergence between A. gingica and A. rufogularis occurred 2.02 (0.91–2.91) mya (Fig. 2).

Figure 2.

Phylogenetic consensus trees from the mtDNA matrix and complete data matrix. Node values above the branches represented the BI posterior probability and ML bootstrap support. Values below the branches represent the divergence times (median) and 95% highest posterior density (HPD) between lineage groups, note that the divergence times in the multi-locus tree were estimated by species tree analysis. The last number in tip labels in the multi-locus tree represent the two haplotypes phased from diploid nuclear sequences.

Mean genetic distances (K2P) between and within subspecies.

Mean distance mtDNA ALDOB FGB G3PDH OVOG nuDNA
Within guangxiensis 0.0013
±0.0005
0.0017
±0.0009
0.0032
±0.0012
0.0037
±0.0017
0.0032
±0.0015
0.0029
±0.0007
Within gingica 0.0012
±0.0005
0.0013
±0.0012
0.0021
±0.0012
0.0009
±0.0009
0.0019
±0.0013
0.0016
±0.0007
Between subspecies 0.0038
±0.0009
0.0017
±0.0011
0.0033
±0.0012
0.0024
±0.0011
0.0038
±0.0018
0.0028
±0.0007

Discussion

This study documents genetic differentiation between A. g. guangxiensis and A. g. gingica. The phylogenetic analyses based on mtDNA indicate that A. g. guangxiensis and A. g. gingica form reciprocal monophyletic groups (Fig. 2), which meets the criterion that subspecies should be monophyletic in mtDNA to demonstrate evolutionary distinctiveness (Zink 2004). Monophyly was also supported by the multi-locus tree (Fig. 2).

However, although the A. g. guangxiensis clade received high support in the mtDNA tree (Fig. 2), nuDNA trees failed to recover the two subspecies as monophyletic groups (Suppl. material 2). This difference might be explained by the longer coalescence time of nuDNA due to its larger effective population size than mtDNA (Moore 1995), so that in recently diverged taxa lineage sorting would be complete for mtDNA but not yet for nuDNA (Zink and Barrowclough 2008). Between A. g. guangxiensis and A. g. gingica, mtDNA haplotypes were fully sorted whereas both taxa had a few shared nuDNA haplotypes. Furthermore, the combined mtDNA and nuDNA tree showed that A. g. guangxiensis and A. g. gingica formed reciprocally monophyletic groups (Fig. 2). Although the monophyly was mainly resolved by mtDNA (Zink and Barrowclough 2008), our results indicate that the two subspecies already exhibit recognizable divergences in nuDNA haplotype frequency, although the divergence was not complete.

In general, molecular phylogenetic study often reveals non-monophyly of avian subspecies (Zink 2004, Phillimore and Owens 2006), which may be due to incorrect taxonomy or rapid divergence. Among Chinese birds, non-monophyly has been documented in some of the subspecies of Lophura nycthemera (Dong et al. 2013), Charadrius alexandrinus (Rheindt et al. 2011) and Garrulax chinensis (Wu et al. 2012), and all subspecies of Motacilla alba (Li et al. 2015) and Leucosticte brandti (Sangster et al. 2015). Thus, the congruent divergence of morphological and molecular markers in A. gingica contrasts with the divergence patterns observed in several other avian species. However, many tropical and subtropical subspecies have been shown to be monophyletic, and sometimes highly divergent, underscoring the necessity of phylogeographic study for taxonomy within species (e.g. Song et al. 2009; Irestedt et al. 2013).

The divergence between A. gingica and A. rufogularis in southwest China (Chen et al. 2015) occurred approximately 2.02 mya when there was a major uplift of the Yunnan-Guizhou Plateau during the Plio-Pleistocene boundary (1.8–2.6 mya) (Zhang and Fang 2012). The uplift may have promoted geographical isolation in many species during this period (Qu et al. 2015), including Stachyridopsis ruficeps (Liu et al. 2012) and Aegithalos concinnus (Dai et al. 2011). The dramatic climatic cooling during the Plio-Pleistocene boundary may have resulted in altitudinal shifts in montane species (Hewitt 2000). These two events may have resulted in the divergence of A. gingica from A. rufogularis, and that between several other species (Liu et al. 2012).

Our results suggest that the divergence between A. g. guangxiensis and A. g. gingica occurred 0.11 (0.05–0.19) mya, during or after the penultimate glaciation (0.13–0.42 mya). We speculate that A. g. guangxiensis and A. g. gingica might have had separate refugia during the glaciation, inducing population differentiation. This Pleistocene refugia scenario has been proposed for several bird species in southeast China, including Tragopan caboti (Dong et al. 2010) and Alcippe morrisonia (Song et al. 2009).

In any case, geographical isolation has likely played a role in population differentiation. A. g. guangxiensis and A. g. gingica are currently separated by the karst basin in central Guangxi. This area also represents unfavorable habitat for some montane species, including Gorsachius magnificus (Hu and Liu 2014) and S. ruficeps (Liu et al. 2012), perhaps due to the large portion of limestone in the karst basin. Early modern human activities dating back to 0.14 mya have been discovered in this area (Shen et al. 2002), and these activities might have interrupted gene flow between bird populations (Zhou and Jiang 2008). These isolation hypotheses may also have affected differentiation between A. g. guangxiensis and A. g. gingica.

The estimated temporal diversification and historical biogeography of A. gingica proposed here is based on a small dataset and thus should ideally be substantiated with additional data. To better explore the underlying diversification process (e.g. speciation-with-gene-flow, Nosil 2008), more sampling using additional nuclear loci is needed (Edwards and Bensch 2009). In addition, ecological niche modelling may help to identify the potential distribution of both subspecies and the main environmental variables which determine the range of each subspecies (Hu and Liu 2014).

Conclusion

Our study demonstrates that the newly found subspecies A. g. guangxiensis and nominate A. g. gingica formed reciprocal monophyletic groups in a multi-locus molecular phylogenetic analyses. The allopatric distribution of A. g. guangxiensis and A. g. gingica and a single diagnostic morphological difference underscore the distinctiveness of these two taxa (Zhou and Jiang 2008). The total body of evidence thus meets the traditional requirement that subspecies are geographically non-overlapping and phenotypically divergent (Wilson and Brown 1953; Mayr and Ashlock 1991) and meets the modern requirements that subspecies are either genetically distinct (Zink 2004), diagnosable (Remsen 2010) or both (Braby et al. 2012). Therefore, our results further support the taxonomic validity of A. g. guangxiensis and we suggest that this subspecies should be considered as a conservation unit.

Acknowledgments

This study was funded by the National Natural Science Foundation of China (No. 31272296), the National Basic Research Program of China (2011BAZ03186) and the China Postdoctoral Science Foundation (2015M580060). We thank Dr Lu Dong and Dr Yue Sun for their help with the collection of samples, Mr Yinong Liu for his help with figure preparation, and Prof. Shou-Hsien Li for his help with manuscript’s revision. We also thank Dr Geoffrey Davison and another anonymous reviewer for their comments on the manuscript, and especially thanks to the editor Dr George Sangster for his comments, edits, and patience.

References

  • Bazinet AL, Zwickl DJ, Cummings MP (2014) A gateway for phylogenetic analysis powered by grid computing featuring GARLI 2.0. Systematic Biology 63: 812–818. doi: 10.1093/sysbio/syu031
  • Braby MF, Eastwood R, Murray N (2012) The subspecies concept in butterflies: has its application in taxonomy and conservation biology outlived its usefulness? Biological Journal of the Linnean Society 106: 699–716. doi: 10.1111/j.1095-8312.2012.01909.x
  • Chen D, Chang J, Li S, Liu Y, Liang W, Zhou F, Yao C, Zhang Z (2015) Was the exposed continental shelf a long-distance colonization route in the ice age? The Southeast Asia origin of Hainan and Taiwan partridges. Molecular Phylogenetics and Evolution 83: 167–173. doi: 10.1016/j.ympev.2014.11.006
  • Cheng T (1978) Fauna Sinica: Aves. Galliformes. Vol, 4. Science Press, Academia Sinica, Beijing, China, 92–93. [In Chinese]
  • Corl A, Ellegren H (2013) Sampling strategies for species trees: The effects on phylogenetic inference of the number of genes, number of individuals, and whether loci are mitochondrial, sex-linked, or autosomal. Molecular Phylogenetics and Evolution 67: 358–366. doi: 10.1016/j.ympev.2013.02.002
  • Dai C, Zhao N, Wang W, Lin C, Gao B, Yang X, Zhang Z, Lei F (2011) Profound climatic effects on two East Asian Black-throated tits (Ave: Aegithalidae), revealed by ecological niche models and phylogeographic analysis. PLoS ONE 6: e29329. doi: 10.1371/journal.pone.0029329
  • Darriba D, Taboada GL, Doallo R, Posada D (2012) jModelTest 2: more models, new heuristics and parallel computing. Nature Methods 9: 772. doi: 10.1038/nmeth.2109
  • del Hoyo J, Collar NJ (2014) HBW and Birdlife International Illustrated Checklist of the Birds of the World. Volume 1: Non-passerines. Lynx Edicions, Barcelona, 88–89.
  • Dickinson EC, Remsen JV (2013) The Howard and Moore Complete Checklist of the Birds of the World. 4th. Edition. Aves Press, Eastbourne, UK, 31–32.
  • Dong L, Heckel G, Liang W, Zhang Y (2013) Phylogeography of Silver Pheasant (Lophura nycthemera L.) across China: aggregate effects of refugia, introgression and riverine barriers. Molecular Ecology 22: 3376–3390. doi: 10.1111/mec.12315
  • Dong L, Zhang J, Sun Y, Liu Y, Zhang Y, Zheng G (2010) Phylogeographic patterns and conservation units of a vulnerable species, Cabot’s tragopan (Tragopan caboti), endemic to southeast China. Conservation Genetics 11: 2231–2242. doi: 10.1007/s10592-010-0108-6
  • Drummond AJ, Rambaut A (2007) BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evolutionary Biology 7: 214. doi: 10.1186/1471-2148-7-214
  • Edgar RC (2004) MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Research 32: 1792–1797. doi: 10.1093/nar/gkh340
  • Edwards SV, Bensch S (2009) Looking forwards or looking backwards in avian phylogeography? A comment on Zink and Barrowclough 2008. Molecular Ecology 18: 2930–2933. doi: 10.1111/j.1365-294X.2009.04270.x
  • Fjeldså J (2013) Avian classification in flux. Lynx Edition, Barcelona, 77–146.
  • He F, Lin Z, Jiang H, Lin J (2007) The occurrence of the White-necklaced Partridge in Fujian of SE China. Chinese Journal of Zoology 42: 147–148. [In Chinese]
  • Heled J, Drummond AJ (2010) Bayesian inference of species trees from multilocus data. Molecular Biology and Evolution 27: 570–580. doi: 10.1093/molbev/msp274
  • Hewitt G (2000) The genetic legacy of the Quaternary ice ages. Nature 405: 907–913. doi: 10.1038/35016000
  • Hu J, Liu Y (2014) Unveiling the conservation biogeography of a data-deficient endangered bird species under climate change. PLoS ONE 9: e84529. doi: 10.1371/journal.pone.0084529
  • Irestedt M, Fabre P-H, Batalha-Filho H, Jønsson KA, Roselaar CS, Sangster G, Ericson PGP (2013) The spatio-temporal colonization and diversification across the Indo-Pacific by a ‘great speciator’ (Aves, Erythropitta erythrogaster). Proceedings of the Royal Society of London B 280: 20130309. doi: 10.1098/rspb.2013.0309
  • Johnsgard PA (1988) The quails, partridges, and francolins of the world. Oxford University Press, Oxford, UK, 221 pp.
  • Li X, Dong F, Lei F, Alström P, Zhang R, Ödeen A, Fjeldså J, Ericson PGP, Zou F, Yang X (2015) Shaped by uneven Pleistocene climate: mitochondrial phylogeographic pattern and population history of White Wagtail Motacilla alba (Aves: Passeriformes). Journal of Avian Biology 47. doi: 10.1111/jav.00826
  • Librado P, Rozas J (2009) DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25: 1451–1452. doi: 10.1093/bioinformatics/btp187
  • Liu H, Wang W, Song G, Qu Y, Li S, Fjeldså J, Lei F (2012) Interpreting the process behind endemism in China by integrating the phylogeography and ecological niche models of the Stachyridopsis ruficeps. PLoS ONE 7: e46761. doi: 10.1371/journal.pone.0046761
  • Madge S, McGowan P (2002) Pheasants, partridges and grouse. Christopher Helm, London, UK, 257 pp.
  • Mayr E (1982) Of what use are subspecies? Auk 99: 593–595.
  • Mayr E, Ashlock PD (1991) Principles of systematic zoology. McGraw-Hill, New York, USA, 43–45.
  • Moore WS (1995) Inferring phylogenies from mtDNA variation: mitochondrial-gene trees versus nuclear-gene trees. Evolution 49: 718–726. doi: 10.2307/2410325
  • Phillimore AB, Owens IPF (2006) Are subspecies useful in evolutionary and conservation biology? Proceedings of the Royal Society of London B 273: 1049–1053. doi: 10.1098/rspb.2005.3425
  • Qu Y, Song G, Gao B, Quan Q, Ericson PGP, Lei F (2015) The influence of geological events on the endemism of East Asian birds studied through comparative phylogeography. Journal of Biogeography 42: 179–192. doi: 10.1111/jbi.12407
  • Remsen JV (2010) Subspecies as a meaningful taxonomic rank in avian classification. Ornithological Monographs 67: 62–78. doi: 10.1525/om.2010.67.1.62
  • Rheindt FE, Székely T, Edwards SV, Lee PLM, Burke T, Kennerley PR, Bakewell DN, Alrashidi M, Kosztolányi A, Weston MA, Liu W, Lei W, Shigeta Y, Javed S, Zefania S, Küpper C (2011) Conflict between Genetic and Phenotypic Differentiation: The Evolutionary History of a ‘Lost and Rediscovered’ Shorebird. PLoS ONE 6: e26995. doi: 10.1371/journal.pone.0026995
  • Rubinoff D, Holland BS (2005) Between two extremes: mitochondrial DNA is neither the panacea nor the nemesis of phylogenetic and taxonomic inference. Systematic Biology 54: 952–961. doi: 10.1080/10635150500234674
  • Sangster G, Roselaar CS, Irestedt M, Ericson PGP (2015) Sillem’s Mountain Finch Leucosticte sillemi is a valid species of rosefinch (Carpodacus, Fringillidae). Ibis 158. doi: 10.1111/ibi.12323
  • Song G, Qu Y, Yin Z, Li S, Liu N, Lei F (2009) Phylogeography of the Alcippe morrisonia (Aves: Timaliidae): long population history beyond late Pleistocene glaciations. BMC Evolutionary Biology 9: 143. doi: 10.1186/1471-2148-9-143
  • Shen G, Wang W, Wang Q, Zhao J, Collerson K, Zhou C, Tobias PV (2002) U-Series dating of Liujiang hominid site in Guangxi, Southern China. Journal of Human Evolution 43: 817–829. doi: 10.1006/jhev.2002.0601
  • Stephens M, Donnelly P (2003) A comparison of Bayesian methods for haplotype reconstruction from population genotype data. The American Journal of Human Genetics 73: 1162–1169. doi: 10.1086/379378
  • Tamura K, Stecher G, Peterson D, Filipski A, Kumar S (2013) MEGA6: Molecular Evolutionary Genetics Analysis Version 6.0. Molecular Biology and Evolution 30: 2725–2729. doi: 10.1093/molbev/mst197
  • Wang N, Kimball RT, Braun EL, Liang B, Zhang Z (2013) Assessing phylogenetic relationships among Galliformes: a multigene phylogeny with expanded taxon sampling in Phasianidae. PLoS ONE 8: e64312. doi: 10.1371/journal.pone.0064312
  • Wilson EO, Brown WL (1953) The subspecies concept and its taxonomic application. Systematic Zoology 2: 97–111. doi: 10.2307/2411818
  • Wu Y, Huang J, Zhang M, Luo S, Zhang Y, Lei F, Sheldon FH, Zou F (2012) Genetic divergence and population demography of the Hainan endemic Black-throated Laughingthrush (Aves: Timaliidae, Garrulax chinensis monachus) and adjacent mainland subspecies. Molecular Phylogenetics and Evolution 65: 482–489. doi: 10.1016/j.ympev.2012.07.005
  • Zhang L, Fang X (2012) Chinese Ancient Geography: The formation of Chinese natural environment. Science Press, Beijing, China, 79 pp. [In Chinese]
  • Zheng G (2011) A Checklist on the Classification and Distribution of the Birds of China (Second Edition). Science Press, Beijing, China, 55–56. [In Chinese]
  • Zhou F (1996) The distribution of White-browed hill partridge (Arborophila gingica) in Guangxi. Acta Zoologica Sinica 42: 143. [In Chinese]
  • Zhou F, Jiang A (2008) A new subspecies of the Collard (sic) hill-partridge. Acta Zootaxonomica Sinica 33: 802–806. [In Chinese]
  • Zink RM (2004) The role of subspecies in obscuring avian biological diversity and misleading conservation policy. Proceedings of the Royal Society of London B 271: 561–564. doi: 10.1098/rspb.2003.2617
  • Zink RM, Barrowclough GF (2008) Mitochondrial DNA under siege in avian phylogeography. Molecular Ecology 17: 2107–2121. doi: 10.1111/j.1365-294X.2008.03737.x