The updated phylogenies of the phasianidae based on combined data of nuclear and mitochondrial DNA

A high-quality assembly reveals genomic characteristics, phylogenetic status, and causal genes for leucism plumage of Indian peafowl

Background

The dazzling phenotypic characteristics of male Indian peafowl (Pavo cristatus) are attractive both to the female of the species and to humans. However, little is known about the evolution of the phenotype and phylogeny of these birds at the whole-genome level. So far, there are no reports regarding the genetic mechanism of the formation of leucism plumage in this variant of Indian peafowl.

Results

A draft genome of Indian peafowl was assembled, with a genome size of 1.05 Gb (the sequencing depth is 362×), and contig and scaffold N50 were up to 6.2 and 11.4 Mb, respectively. Compared with other birds, Indian peafowl showed changes in terms of metabolism, immunity, and skeletal and feather development, which provided a novel insight into the phenotypic evolution of peafowl, such as the large body size and feather morphologies. Moreover, we determined that the phylogeny of Indian peafowl was more closely linked to turkey than chicken. Specifically, we first identified that PMEL was a potential causal gene leading to the formation of the leucism plumage variant in Indian peafowl.

Conclusions

This study provides an Indian peafowl genome of high quality, as well as a novel understanding of phenotypic evolution and phylogeny of Indian peafowl. These results provide a valuable reference for the study of avian genome evolution. Furthermore, the discovery of the genetic mechanism for the development of leucism plumage is both a breakthrough in the exploration of peafowl plumage and also offers clues and directions for further investigations of the avian plumage coloration and artificial breeding in peafowl.

Exploring potentialities of avian genomic research in Nepalese Himalayas

Nepal, a small landlocked country in South Asia, holds about 800 km of Himalayan Mountain range including the Earth's highest mountain. Within such a mountain range in the north and plain lowlands in the south, Nepal provides a habitat for about 9% of global avian fauna. However, this diversity is underrated because of the lack of enough studies, especially using molecular tools to quantify and understand the distribution patterns of diversity. In this study, we reviewed the studies in the last two decades (2000‒2019) that used molecular methods to study the biodiversity in Nepal to examine the ongoing research trend and focus. Although Nepalese Himalaya has many opportunities for cutting-edge molecular research, our results indicated that the rate of genetic/genomic studies is much slower compared to the regional trends. We found that genetic research in Nepal heavily relies on resources from international institutes and that too is mostly limited to research on species monitoring, distribution, and taxonomic validations. Local infrastructures to carry out cutting-edge genomic research in Nepal are still in their infancy and there is a strong need for support from national/international scientists, universities, and governmental agencies to expand such genomic infrastructures in Nepal. We particularly highlight avian fauna as a potential future study system in this region that can be an excellent resource to explore key biological questions such as understanding eco-physiology and molecular basis of organismal persistence to changing environment, evolutionary processes underlying divergence and speciation, or mechanisms of endemism and restrictive distribution of species.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1186/s40657-021-00290-5.

When good mitochondria go bad: Cyto-nuclear discordance in landfowl (Aves: Galliformes)

Mitochondrial sequences were among the first molecular data collected for phylogenetic studies and they are plentiful in DNA sequence archives. However, the future value of mitogenomic data in phylogenetics is uncertain, because its phylogenetic signal sometimes conflicts with that of the nuclear genome. A thorough understanding of the causes and prevalence of cyto-nuclear discordance would aid in reconciling different results owing to sequence data type, and provide a framework for interpreting megaphylogenies when taxa which lack substantial nuclear data are placed using mitochondrial data. Here, we examine the prevalence and possible causes of cyto-nuclear discordance in the landfowl (Aves: Galliformes), leveraging 47 new mitogenomes assembled from off-target reads recovered as part of a target-capture study. We evaluated two hypotheses, that cyto-nuclear discordance is "genuine" and a result of biological processes such as incomplete lineage sorting or introgression, and that cyto-nuclear discordance is an artifact of inaccurate mitochondrial tree estimation (the "inaccurate estimation" hypothesis). We identified seven well-supported topological differences between the mitogenomic tree and trees based on nuclear data. These well-supported topological differences were robust to model selection. An examination of sites suggests these differences were driven by small number of sites, particularly from third-codon positions, suggesting that they were not confounded by convergent directional selection. Hence, the hypothesis of genuine discordance was supported.

The draft genome of the blood pheasant (Ithaginis cruentus): Phylogeny and high-altitude adaptation

The blood pheasant (Ithaginis cruentus), the only species in the genus Ithaginis, lives in an extremely inhospitable high-altitude environment, coping with hypoxia and ultraviolet (UV) radiation. To further investigate the phylogeny of Phasianidae species based on complete genomes and understand the molecular genetic mechanisms of the high-altitude adaptation of the blood pheasant, we de novo assembled and annotated the complete genome of the blood pheasant. The blood pheasant genome size is 1.04 Gb with scaffold N50 of 10.88 Mb. We identified 109.92 Mb (10.62%) repetitive elements, 279,037 perfect microsatellites, and 17,209 protein-coding genes. The phylogenetic tree of Phasianidae based on whole genomes revealed three highly supported major clades with the blood pheasant included in the "erectile clade." Comparative genomics analysis showed that many genes were positively selected in the blood pheasant, which was associated with response to hypoxia and/or UV radiation. More importantly, among these positively selected genes (PSGs) which were related to high-altitude adaptation, sixteen PSGs had blood pheasant-specific missense mutations. Our data and analysis lay solid foundation to the study of Phasianidae phylogeny and provided new insights into the potential adaptation mechanisms to the high altitude employed by the blood pheasant.

Comparison of the Chinese bamboo partridge and red Junglefowl genome sequences highlights the importance of demography in genome evolution

BACKGROUND

Recent large-scale whole genome sequencing efforts in birds have elucidated broad patterns of avian phylogeny and genome evolution. However, despite the great interest in economically important phasianids like Gallus gallus (Red Junglefowl, the progenitor of the chicken), we know little about the genomes of closely related species. Gallus gallus is highly sexually dichromatic and polygynous, but its sister genus, Bambusicola, is smaller, sexually monomorphic, and monogamous with biparental care. We sequenced the genome of Bambusicola thoracicus (Chinese Bamboo Partridge) using a single insert library to test hypotheses about genome evolution in galliforms. Selection acting at the phenotypic level could result in more evidence of positive selection in the Gallus genome than in Bambusicola. However, the historical range size of Bambusicola was likely smaller than Gallus, and demographic effects could lead to higher rates of nonsynonymous substitution in Bambusicola than in Gallus.

RESULTS

We generated a genome assembly suitable for evolutionary analyses. We examined the impact of selection on coding regions by examining shifts in the average nonsynonymous to synonymous rate ratio (dN/dS) and the proportion of sites subject to episodic positive selection. We observed elevated dN/dS in Bambusicola relative to Gallus, which is consistent with our hypothesis that demographic effects may be important drivers of genome evolution in Bambusicola. We also demonstrated that alignment error can greatly inflate estimates of the number of genes that experienced episodic positive selection and heterogeneity in dN/dS. However, overall patterns of molecular evolution were robust to alignment uncertainty. Bambusicola thoracicus has higher estimates of heterozygosity than Gallus gallus, possibly due to migration events over the past 100,000 years.

CONCLUSIONS

Our results emphasized the importance of demographic processes in generating the patterns of variation between Bambusicola and Gallus. We also demonstrated that genome assemblies generated using a single library can provide valuable insights into avian evolutionary history and found that it is important to account for alignment uncertainty in evolutionary inferences from draft genomes.

A novel mitochondrial genome of Arborophila and new insight into Arborophila evolutionary history

The lineage of the Bar-backed Partridge (Arborophila brunneopectus) was investigated to determine the phylogenetic relationships within Arborophila as the species is centrally distributed within an area covered by the distributions of 22 South-east Asian hill partridge species. The complete mitochondrial genome (mitogenome) of A. brunneopectus was determined and compared with four other hill partridge species mitogenomes. NADH subunit genes are radical in hill partridge mitogenomes and contain the most potential positive selective sites around where variable sites are abundant. Together with 44 other mitogenomes of closely related species, we reconstructed highly resolved phylogenetic trees using maximum likelihood (ML) and Bayesian inference (BI) analyses and calculated the divergence and dispersal history of Arborophila using combined datasets composed of their 13-protein coding sequences. Arborophila is reportedly be the oldest group in Phasianidae whose ancestors probably originated in Asia. A. rufipectus shares a closer relationship with A. ardens and A. brunneopectus compared to A. gingica and A. rufogularis, and such relationships were supported and profiled by NADH dehydrogenase subunit 5 (ND5). The intragenus divergence of all five Arborophila species occurred in the Miocene (16.84~5.69 Mya) when there were periods of climate cooling. We propose that these cooling events in the Miocene forced hill partridges from higher to lower altitudes, which led to geographic isolation and speciation. We demonstrated that the apparently deleterious +1 frameshift mutation in NADH dehydrogenase subunit 3 (ND3) found in all Arborophila is an ancient trait that has been eliminated in some younger lineages, such as Passeriformes. It is unclear of the biological advantages of this elimination for the relevant taxa and this requires further investigation.

Complete mitochondrial genome of Francolinus pintadeanus (Galliformes: Phasianidae)

The complete mitochondrial genome (mitogenome) of Francolinus pintadeanus, consisting of 16,693 bp, was determined and analyzed. It displayed as typical genome organization as other Galliformes mitogenomes: 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes and one control region. The phylogenetic relationships of 25 Phasianidae species and three Odontophoridae species as outgroup using maximum likelihood and Bayesian inference methods based on a concatenated dataset from mitogenomes were analyzed. The results reveal that F. pintadeanus had a close relationship with Gallus gallus gallus/Bambusicola thoracica, then this clade formed a sister group with Pavo muticus/Argusianus argus.

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.

Effects of Supplementary Feeding on the Breeding Ecology of the Buff-Throated Partridge in a Tibetan Sacred Site, China

Our goal was to document effects of year-round supplemental feeding on breeding ecology of the Buff-throated Partridge, Tetraophasis szechenyii, within a Tibetan sacred site. We evaluated effects of supplemental feeding used as religious/cultural practices which could potentially aid conservation of endangered phasianids. We compared fed breeding groups to neighboring nonfed groups. Fed groups initiated first clutches significantly earlier than nonfed groups. Earlier laying groups within fed and nonfed groups showed significantly lower hatching rates than later groups; however, fed groups showed significantly higher hatching rates than nonfed groups laying in the same period. Earlier laying increased opportunities to renest. All six fed groups with clutch failures renested compared to only one of five nonfed groups with clutch failures. Fed female breeders showed significantly greater investment in their young with larger clutches and larger eggs, which likely increased survivability of early hatchlings. We observed no predation on birds at feeding sites and recorded only four cases of predation on incubating females, which showed no detectable difference between fed and nonfed groups. Ground-nesting birds typically face high risks of predation. Ten of the 48 groups nested in trees, which occurs in few phasianid species. Tree nests showed significantly higher hatching rates compared to ground nests; however, we found no significant difference in tree nesting between fed and nonfed groups. This partridge is one of four gallinaceous species with cooperative breeding. Breeding groups with helpers had significantly greater reproductive success than single pairs, and fed female breeders with helpers laid bigger eggs than single pairs. Comparing annual reproductive output per group, fed groups not only produced significantly more independent young (≥150 days post-hatching), their young hatched significantly earlier, which likely have greater reproductive value over later hatched young of nonfed groups. Supplemental feeding year-round is likely what enabled the successes of the fed partridges.

Comparative mitochondrial genomics and phylogenetic relationships of the Crossoptilon species (Phasianidae, Galliformes)

BACKGROUND

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.

RESULTS

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.

CONCLUSIONS

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.