A phylogenomic supermatrix of Galliformes (Landfowl) reveals biased branch lengths.

Building taxon-rich phylogenies is foundational for macroevolutionary studies. One approach to improve taxon sampling beyond individual studies is to build supermatricies of publicly available data, incorporating taxa sampled across different studies and utilizing different loci. Most existing supermatrix studies have focused on loci commonly sequenced with Sanger technology ("legacy" markers, such as mitochondrial data and small numbers of nuclear loci). However, incorporating phylogenomic studies into supermatrices allows problem nodes to be targeted and resolved with considerable amounts of data, while improving taxon sampling with legacy data. Here we estimate phylogeny from a galliform supermatrix which includes well-known model and agricultural species such as the chicken and turkey. We assembled a supermatrix comprising 4500 ultra-conserved elements (UCEs) collected as part of recent phylogenomic studies in this group and legacy mitochondrial and nuclear (intron and exon) sequences. Our resulting phylogeny included 88% of extant species and recovered well-accepted relationships with strong support. However, branch lengths, which are particularly important in down-stream macroevolutionary studies, appeared vastly skewed. Taxa represented only by rapidly evolving mitochondrial data had high proportions of missing data and exhibited long terminal branches. Conversely, taxa sampled for slowly evolving UCEs with low proportions of missing data exhibited substantially shorter terminal branches. We explored several branch length re-estimation methods with particular attention to terminal branches and conclude that re-estimation using well-sampled mitochondrial sequences may be a pragmatic approach to obtain trees suitable for macroevolutionary analysis.

Ongoing green peafowl protection in China.

The green peafowl (Pavo muticus) is one of the most threatened pheasants in the world. In China, it is widely appreciated for its beauty as well as historical and cultural value, but current populations number less than 500 individuals. Recently, Tang and colleagues reported in Science that the green peafowl is likely to become extinct due to the construction of the Jiasajiang Level 1 Hydropower Station within the Red River Upstream District (RRUD) and thus called for a stop to this project (Tang et al., 2019). According to our recent surveys, however, this species is still extant in 22 counties of Yunnan Province, China, among which, only two within the RRUD have been predicted to be affected by floods from the hydropower station. Therefore, the conclusion that the species will likely go extinct in China upon completion of the dam is unwarranted. In fact, construction of the Jiasajiang Level 1 Hydropower Station was stopped in August 2017. The main challenge for green peafowl conservation is that over 65% of the population occurs outside of protected reserves in China. Fortunately, the Chinese government has adopted an Ecological Redline (ERL) strategy to achieve ecological civilization plans, thus bringing new hope to the conservation of green peafowls both inside and outside of protected reserves. As a top conservation priority for China, the government is fully committed to conserving this peafowl.

Symmetry breaking in the embryonic skin triggers directional and sequential plumage patterning.

The development of an organism involves the formation of patterns from initially homogeneous surfaces in a reproducible manner. Simulations of various theoretical models recapitulate final states of natural patterns, yet drawing testable hypotheses from those often remains difficult. Consequently, little is known about pattern-forming events. Here, we surveyed plumage patterns and their emergence in Galliformes, ratites, passerines, and penguins, together representing the three major taxa of the avian phylogeny, and built a unified model that not only reproduces final patterns but also intrinsically generates shared and varying directionality, sequence, and duration of patterning. We used in vivo and ex vivo experiments to test its parameter-based predictions. We showed that directional and sequential pattern progression depends on a species-specific prepattern: an initial break in surface symmetry launches a travelling front of sharply defined, oriented domains with self-organising capacity. This front propagates through the timely transfer of increased cell density mediated by cell proliferation, which controls overall patterning duration. These results show that universal mechanisms combining prepatterning and self-organisation govern the timely emergence of the plumage pattern in birds.

Mitochondrial DNA variation of the ruffed grouse (Bonasa umbellus).

OBJECTIVE: The ruffed grouse, Bonasa umbellus, is broadly distributed across North America and displays considerable taxonomic diversity. Except for a genetic study of some western populations of ruffed grouse, nothing is known about genetic variation in other regions of Canada and the United States. Our objective is to examine patterns of mitochondrial DNA (mtDNA) variation in the ruffed grouse across western, central, and eastern parts of its distribution. We compare patterns of mtDNA variation to those characterized by morphology and ecology. Additionally, we evaluate the demographic history of the species based on mitochondrial haplotype diversity. RESULTS: Patterns of mtDNA variation revealed geographic subdivision, with populations of ruffed grouse subdivided into 3 to 4 genetically distinct groups. This subdivision partially coincided with the ranges of described subspecies. Behavioral traits prohibiting long-distance movement and barriers to dispersal in response to physiography and unsuitable habitat help explain these patterns of subdivision. Historically, the ruffed grouse probably experienced a population expansion, possibly in response to changes during the Pleistocene.

The Draft Genome of the Endangered Sichuan Partridge (Arborophila rufipectus) with Evolutionary Implications.

The Sichuan partridge (Arborophila rufipectus, Phasianidae, Galliformes) is distributed in south-west China, and classified as endangered grade. To examine the evolution and genomic features of Sichuan partridge, we de novo assembled the Sichuan partridge reference genome. The final draft assembly consisted of approximately 1.09 Gb, and had a scaffold N50 of 4.57 Mb. About 1.94 million heterozygous single-nucleotide polymorphisms (SNPs) were detected, 17,519 protein-coding genes were predicted, and 9.29% of the genome was identified as repetitive elements. A total of 56 olfactory receptor (OR) genes were found in Sichuan partridge, and conserved motifs were detected. Comparisons between the Sichuan partridge genome and chicken genome revealed a conserved genome structure, and phylogenetic analysis demonstrated that Arborophila possessed a basal phylogenetic position within Phasianidae. Gene Ontology (GO) enrichment analysis of positively selected genes (PSGs) in Sichuan partridge showed over-represented GO functions related to environmental adaptation, such as energy metabolism and behavior. Pairwise sequentially Markovian coalescent analysis revealed the recent demographic trajectory for the Sichuan partridge. Our data and findings provide valuable genomic resources not only for studying the evolutionary adaptation, but also for facilitating the long-term conservation and genetic diversity for this endangered species.

De novo assembly of the Indian blue peacock (Pavo cristatus) genome using Oxford Nanopore technology and Illumina sequencing.

BACKGROUND: The Indian peafowl (Pavo cristanus) is native to South Asia and is the national bird of India. Here we present a draft genome sequence of the male blue peacock using Illumina and Oxford Nanopore technology (ONT). RESULTS: ONT sequencing gave ∼2.3-fold sequencing coverage, whereas Illumina generated 150-base pair paired-end sequence data at 284.6-fold coverage from 5 libraries. Subsequently, we generated a 0.915-gigabase pair de novo assembly of the peacock genome with a scaffold N50 of 0.23 megabase pairs (Mb). We predict that the peacock genome contains 23,153 protein-coding genes and 75.3 Mb (7.33%) of repetitive sequences. CONCLUSIONS: We report a high-quality assembly of the peacock genome using a hybrid approach of sequences generated by both Illumina and ONT. The long-read chemistry generated by ONT was useful for addressing challenges related to de novo assembly, particularly at regions containing repetitive sequences spanning longer than the read length, and which could not be resolved with only short-read-based assembly. Contig assembly of Illumina short reads gave an N50 of 1,639 bases, whereas with ONT, the N50 increased by >9-fold to 14,749 bases. The initial contig assembly based on Illumina sequencing reads alone gave 685,241 contigs. Further scaffolding on assembled contigs using both Illumina and ONT sequencing reads resulted in a final assembly of 15,025 super-scaffolds, with an N50 of ∼0.23 Mb. Ninety-five percent of proteins predicted by homology matched with those in a public repository, verifying the completeness of our assembly. Like other phylogenetic studies of avian conserved genes, we found P. cristatus to be most closely related to Gallus gallus, followed by Meleagris gallopavo and Anas platyrhynchos. Compared with the recently published peacock genome assembly, the current, superior, hybrid assembly has greater sequencing depth, fewer non-ATGC sequences, and fewer scaffolds.

A guinea fowl genome assembly provides new evidence on evolution following domestication and selection in galliformes.

The helmeted guinea fowl Numida meleagris belongs to the order Galliformes. Its natural range includes a large part of sub-Saharan Africa, from Senegal to Eritrea and from Chad to South Africa. Archaeozoological and artistic evidence suggest domestication of this species may have occurred about 2,000 years BP in Mali and Sudan primarily as a food resource, although villagers also benefit from its capacity to give loud alarm calls in case of danger, of its ability to consume parasites such as ticks and to hunt snakes, thus suggesting its domestication may have resulted from a commensal association process. Today, it is still farmed in Africa, mainly as a traditional village poultry, and is also bred more intensively in other countries, mainly France and Italy. The lack of available molecular genetic markers has limited the genetic studies conducted to date on guinea fowl. We present here a first-generation whole-genome sequence draft assembly used as a reference for a study by a Pool-seq approach of wild and domestic populations from Europe and Africa. We show that the domestic populations share a higher genetic similarity between each other than they do to wild populations living in the same geographical area. Several genomic regions showing selection signatures putatively related to domestication or importation to Europe were detected, containing candidate genes, most notably EDNRB2, possibly explaining losses in plumage coloration phenotypes in domesticated populations.

Genomic regions of speciation and adaptation among three species of grouse.

Understanding the molecular basis of adaption is one of the central goals in evolutionary biology and when investigated across sister species it can provide detailed insight into the mechanisms of speciation. Here, we sequence the genomes of 34 individuals from three closely related grouse species in order to uncover the genomic architecture of speciation and the genes involved in adaptation. We identify 6 regions, containing 7 genes that show lineage specific signs of differential selection across the species. These genes are involved in a variety of cell processes ranging from stress response to neural, gut, olfactory and limb development. Genome wide neutrality test statistics reveal a strong signal of population expansion acting across the genomes. Additionally, we uncover a 3.5 Mb region on chromosome 20 that shows considerably lower levels of differentiation across the three grouse lineages, indicating possible action of uniform selection in this region.

The first draft genome of Lophophorus: A step forward for Phasianidae genomic diversity and conservation.

The monal genus (Lophophorus) is a branch of Phasianidae and its species inhabit the high-altitude mountains of the Qinghai-Tibet Plateau. The Chinese monal, L. lhuysii, is a threatened endemic bird of China that possesses high-altitude adaptability, diversity of plumage color and potentially low reproductive life history. This is the first study to describe the monal genome using next generation sequencing technology. The Chinese monal genome size is 1.01 Gb, with 16,940 protein-coding genes. Gene annotation yielded 100.93 Mb (9.97%) repeat elements, 785 ncRNA, 5,465,549 bp (0.54%) SSR and 15,550 (92%) genes in public databases. Compared to other birds and mammals, the genome evolution analysis showed numerous expanded gene families and positive selected genes involved in high-altitude adaptation, especially related to the adaptation of low temperature and hypoxia. Consequently, this gene data can be used to investigate the molecular evolution of high-altitude adaptation in future bird research. Our first published genome of the genus Lophophorus will be integral for the study of monal population genetic diversity and conservation, genomic evolution and Galliformes species differentiation in the Qinghai-Tibetan Plateau.

Human-modified biogeographic patterns and conservation in game birds: The dilemma of the black francolin (Francolinus francolinus, Phasianidae) in Pakistan.

The ever-increasing human-mediated wildlife reshuffling is raising concern for the conservation of biodiversity. The loss of biological distinctiveness among regions lessens the genetic diversity and consequently the evolutionary potential of local biotas to tackle present-day global change and human disturbance. This process may be sometimes cryptic unless investigated by means of a molecular approach. In this respect, game birds are a paradigmatic case. The black francolin (Francolinus francolinus, Phasianidae) is a medium-sized galliform whose distribution range stretches from Cyprus to the Gulf of Bengal. Six morphologic subspecies are known, with three of which occurring in Pakistan, where the species is heavily hunted and used as pet for chirping competitions. We genotyped 98 samples (feathers) at both the entire mitochondrial DNA Control Region gene and nine microsatellite loci to get a deeper insight into the genetic diversity of the black francolin in Pakistan in order to offer cogent recommendations for its conservation management. We identified several mtDNA lineages that were consistent with the currently described subspecies/taxonomy whose pattern of co-occurrence is compatible with the geological history and the faunal movement routes of the region under study. However, the biparentally inherited microsatellites returned a quite discordant picture of an extensive, sex-biased genetic mixing due to the intensive relocations of already overharvested male individuals for chirping competitions. Our results indicated that the genetic integrity of the black francolin in Pakistan could be seriously at risk and call for monitoring and limiting its trade other than enhancing the public awareness of the importance of local biodiversity resources.

The periodic coloration in birds forms through a prepattern of somite origin.

The periodic stripes and spots that often adorn animals' coats have been largely viewed as self-organizing patterns, forming through dynamics such as Turing's reaction-diffusion within the developing skin. Whether preexisting positional information also contributes to the periodicity and orientation of these patterns has, however, remained unclear. We used natural variation in colored stripes of juvenile galliform birds to show that stripes form in a two-step process. Autonomous signaling from the somite sets stripe position by forming a composite prepattern marked by the expression profile of agouti Subsequently, agouti regulates stripe width through dose-dependent control of local pigment production. These results reveal that early developmental landmarks can shape periodic patterns upstream of late local dynamics, and thus constrain their evolution.

Crumble analysis of the historic sympatric distribution between Dendrortyx macroura and D. barbatus (Aves: Galliformes).

In Mexico, the Long-tailed Wood-Partridge (Dendrortyx macroura) is distributed in the mountains of the Trans-Mexican Volcanic Belt, Sierra Madre del Sur and Sierra Norte de Oaxaca; while the Bearded Wood-Partridge (D. barbatus) is distributed in the Sierra Madre Oriental (SMO). There is a controversial overlap in distribution (sympatry) between these two species (on the Cofre de Perote and Pico de Orizaba volcanoes, SMO and Sierra Norte de Oaxaca), based on the ambiguity and current lack of information regarding the distribution of these two species. In order to disentangle the possible presence of both species in the area of sympatry, we conducted a crumble analysis of the historic knowledge regarding the geographic distribution of both species, based on a review of scientific literature, database records, the specimen examination (in ornithological collections), field work and a reconstruction of the distribution range based on Ecological Niche Modeling. Our results support the presence of only one of these two species in the overlapping area, rejecting the existence of such an area of sympatry between the two species. We discuss alternative hypotheses that could explain the historically reported distribution pattern: 1) an error in the single existing historical record; 2) a possible local extinction of the species and 3) the past existence of interspecific competition that has since been resolved under the principle of competitive exclusion. We propose that the Santo Domingo River in northern Oaxaca and western slope of the Sierra Madre Oriental, mark the distribution limits between these species.

How do seemingly non-vagile clades accomplish trans-marine dispersal? Trait and dispersal evolution in the landfowl (Aves: Galliformes).

Dispersal ability is a key factor in determining insular distributions and island community composition, yet non-vagile terrestrial organisms widely occur on oceanic islands. The landfowl (pheasants, partridges, grouse, turkeys, quails and relatives) are generally poor dispersers, but the Old World quail (Coturnix) are a notable exception. These birds evolved small body sizes and high-aspect-ratio wing shapes, and hence are capable of trans-continental migrations and trans-oceanic colonization. Two monotypic partridge genera, Margaroperdix of Madagascar and Anurophasis of alpine New Guinea, may represent additional examples of trans-marine dispersal in landfowl, but their body size and wing shape are typical of poorly dispersive continental species. Here, we estimate historical relationships of quail and their relatives using phylogenomics, and infer body size and wing shape evolution in relation to trans-marine dispersal events. Our results show that Margaroperdix and Anurophasis are nested within the Coturnix quail, and are each 'island giants' that independently evolved from dispersive, Coturnix-like ancestral populations that colonized and were subsequently isolated on Madagascar and New Guinea. This evolutionary cycle of gain and loss of dispersal ability, coupled with extinction of dispersive taxa, can result in the false appearance that non-vagile taxa somehow underwent rare oceanic dispersal.

Conservation genetics of the capercaillie in Poland - Delineation of conservation units.

The capercaillie (Tetrao urogallus) is one of Poland's most endangered bird species, with an estimated population of 380-500 individuals in four isolated areas. To study these natural populations in Poland further, more than 900 non-invasive genetic samples were collected, along with samples from 59 birds representing large, continuous populations in Sweden and Russia; and from two centres in Poland breeding capercaillie. Microsatellite polymorphism at nine loci was then analysed to estimate within-population genetic diversity and genetic differentiation among populations. The results confirmed that isolation of populations and recent decreases in their sizes have reduced genetic diversity among capercaillie in Poland, with all the country's natural populations found to be experiencing the genetic after-effects of demographic bottlenecks. The results of analyses of genetic differentiation and structure further suggest the presence of a 'lowland' cluster (encompassing birds of the Augustowska and Solska Primaeval Forests in Poland, and of Sweden and Russia), and a Carpathian cluster. Capercaillie from Sweden and Russia are also found to differ markedly. The Polish lowland populations seem more closely related to birds from Scandinavia. Our genetic analysis also indicates that the stocks at breeding centres are of a high genetic diversity effectively reflecting the origins of founder individuals, though identification of ancestry requires further study in the case of some birds. Overall, the results sustain the conclusion that the Polish populations of capercaillie from the Carpathians and the lowlands should be treated as independent Management Units (MUs). This is to say that the breeding lines associated with these two sources should be maintained separately at breeding centres. The high level of genetic differentiation of birds from the Solska Primaeval Forest suggests that this population should also be assigned the status of independent MU.

Historical relationships of three enigmatic phasianid genera (Aves: Galliformes) inferred using phylogenomic and mitogenomic data.

The phylogeny of the Phasianidae (pheasants, partridges, and allies) has been studied extensively. However, these studies have largely ignored three enigmatic genera because of scarce DNA source material and limited overlapping phylogenetic data: blood pheasants (Ithaginis), snow partridges (Lerwa), and long-billed partridges (Rhizothera). Thus, phylogenetic positions of these three genera remain uncertain in what is otherwise a well-resolved phylogeny. Previous studies using different data types place Lerwa and Ithaginis in similar positions, but the absence of overlapping data means the relationship between them could not be inferred. Rhizothera was originally described in the genus Perdix (true partridges), although a partial cytochrome b (CYB) sequence suggests it is sister to Pucrasia (koklass pheasant). To identify robust relationships among Ithaginis, Lerwa, Rhizothera, and their phasianid relatives, we used 3692 ultra-conserved element (UCE) loci and complete mitogenomes from 19 species including previously hypothesized relatives of the three focal genera and representatives from all major phasianid clades. We used DNA extracted from historical specimen toepads for species that lacked fresh tissue in museum collections. Maximum likelihood and multispecies coalescent UCE analyses strongly supported Lerwa sister to a large clade which included Ithaginis at its base, and also including turkey, grouse, typical pheasants, tragopans, Pucrasia, and Perdix. Rhizothera was also in this clade, sister to a diverse group comprising Perdix, typical pheasants, Pucrasia, turkey and grouse. Mitogenomic genealogies differed from UCEs topologies, supporting a sister relationship between Ithaginis and Lerwa rather than a grade. The position of Rhizothera using mitogenomes depended on analytical choices. Unpartitioned and codon-based analyses placed Rhizothera sister to a tragopan clade, whereas a partitioned DNA model of the mitogenome was congruent with UCE results. In all mitogenome analyses, Pucrasia was sister to a clade including Perdix and the typical pheasants with high support, in contrast to UCEs and published nuclear intron data. Due to the strong support and consistent topology provided by all UCE analyses, we have identified phylogenetic relationships of these three enigmatic, poorly-studied, phasianid taxa.

Recovering the Genetic Identity of an Extinct-in-the-Wild Species: The Puzzling Case of the Alagoas Curassow.

The conservation of many endangered taxa relies on hybrid identification, and when hybrids become morphologically indistinguishable from the parental species, the use of molecular markers can assign individual admixture levels. Here, we present the puzzling case of the extinct in the wild Alagoas Curassow (Pauxi mitu), whose captive population descends from only three individuals. Hybridization with the Razor-billed Curassow (P. tuberosa) began more than eight generations ago, and admixture uncertainty affects the whole population. We applied an analysis framework that combined morphological diagnostic traits, Bayesian clustering analyses using 14 microsatellite loci, and mtDNA haplotypes to assess the ancestry of all individuals that were alive from 2008 to 2012. Simulated data revealed that our microsatellites could accurately assign an individual a hybrid origin until the second backcross generation, which permitted us to identify a pure group among the older, but still reproductive animals. No wild species has ever survived such a severe bottleneck, followed by hybridization, and studying the recovery capability of the selected pure Alagoas Curassow group might provide valuable insights into biological conservation theory.

Phylogeography of chukar partridge (Alectoris chukar) in China based on mtDNA control region.

The chukar partridge (Alectoris chukar) is distributed in north and northwest of China, in mountainous areas that were heavily affected by cyclic climate and landscape changes during the last Pleistocene glaciations. Some partridge populations have colonized and expanded their present ranges only after deglaciation and recent deforestation by human. Consequently, partridges from different areas could be genetically differentiated. In this study, a 1152-1154 bp portion of the mitochondrial DNA control region were analyzed for all 279 specimens collected from 28 populations through their distribution in China. A total number of 91 haplotypes were defined by 113 variable sites. The mean haplotype diversity and nucleotide diversity were 0.939 ± 0.008 and 0.0030 ± 0.0017, respectively, for chukar. Haplotype diversity among the 28 populations varied from 0.600 ± 0.215 (HH) to 1.000 ± 0.272 (SBC), and nucleotide diversity ranged from 0.0006 ± 0.0007 (HJ) to (0.0071 ± 0.0041) (HEG). The nucleotide diversity of chukar was descending from west to east and a possible historical scenario might be that A. chukar inhabited central Asia, and then dispersed eastward. Phylogenetic analyses indicated that there was no distinct phylogeographic structure in chukar populations in China. Haplotype network of chukar was star-like with some common haplotypes shared by different samples came from different populations as center. Both unstructured phylogeographic tree and star-like haplotype network are signatures of population expansion.

Rapid and recent diversification of curassows, guans, and chachalacas (Galliformes: Cracidae) out of Mesoamerica: Phylogeny inferred from mitochondrial, intron, and ultraconserved element sequences.

The Cracidae (curassows, guans, and chachalacas) include some of the most spectacular and endangered Neotropical bird species. They lack a comprehensive phylogenetic hypothesis, hence their geographic origin and the history of their diversification remain unclear. We present a species-level phylogeny of Cracidae inferred from a matrix of 430 ultraconserved elements (UCEs; at least one species sampled per genus) and eight more variable loci (introns and mtDNA; all available species). We use this phylogeny along with probabilistic biogeographic modeling to test whether Gondwanan vicariance, ancient dispersal to South America, ancient dispersal from South America, or massive global cooling isolated cracids in the Neotropics. Contrary to previous estimates that extant cracids diversified in the Cretaceous, our fossil-calibrated divergence time estimates instead support that crown Cracidae originated in the late Miocene. Species-rich genera Crax, Penelope, and Ortalis began diversifying as recently as 3Mya. Biogeographic reconstructions indicate that modern cracids originated in Mesoamerica and were isolated from a widespread Laurasian ancestor, consistent with the massive global cooling hypothesis. Current South American diversity is the result of multiple colonization events following uplift of the Panamanian Isthmus, coupled with rapid diversification and evolution of secondary sympatry. Of the four major cracid lineages (curassows, chachalacas, typical guans, horned guan), the only lineage that has failed to colonize and diversify South America is the unique horned guan (Oreophasis derbianus), which is sister to curassows and chachalacas rather than typical guans.

Complete sequence determination and analysis of the Guangxi Phasianus colchicus (Galliformes: Phasianidae) mitochondrial genome.

The complete mitochondrial genome sequence of the Guangxi Phasianus colchicus was measured by PCR-based methods and analyzed in detail. Our research findings reveal that the entire mitochondrial genome of P. colchicus is a circular molecule of length 16,687 bp (GenBank accession number: KT364526). The contents of A, T, C, and G in the mitochondrial genome were found to be 30.64%, 25.29%, 30.81% and 13.26%, respectively. The complete mitochondrial genome of the P. colchicus is a typical structure, including 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes, and 1 control region (D-loop region). This complete mitochondrial genome sequence provides essential information in understanding phylogenetic relationships among Galliformes mitochondrial genomes.

Identification of CR1 retroposons in Arborophila rufipectus and their application to Phasianidae phylogeny.

Chicken repeat 1 (CR1), a member of non-LTR retroposon, is an important phylogenetic marker in avian systematics. In this study, we reported several characteristics of CR1 elements in a draft genome of Arborophila rufipectus (Sichuan partridge). According to the analyses of RepeatMasker, approximately 254 966 CR1 elements were identified in A. rufipectus, covering 6.7% of the genome. Subsequently, we selected eighteen novel CR1 elements by comparing the chicken genome, turkey genome and assembled A. rufipectus scaffolds. Here, a combined data set comprising of 22 CR1 loci, mitochondrial genomes and eight unlinked introns was analysed to infer the evolutionary relationships of twelve Phasianidae species. The applicability of CR1 sequences for inferring avian phylogeny relative to mtDNA and intron sequences was investigated as well. Our results elucidated the position of A. rufipectus in Phasianidae with robust supports that it presented a sister clade to Arborophila ardens/Arborophila brunneopectus, and implied that genus Arborophila was in a basal phylogenetic position within Phasianidae and a phylogenetic affinity between Meleagris gallopavo and Pucrasia macrolopha. Therefore, this work not only resolved some of the confounding relationships among Phasianidae, but also suggested CR1 sequences could provide powerful complementary data for phylogeny reconstruction.