Next Article in Journal
In Vitro Effects of Photon Beam and Carbon Ion Radiotherapy on the Perineural Invasion of Two Cell Lines of Neurotropic Tumours
Previous Article in Journal
Late-Onset Psoriatic Arthritis: Are There Any Distinct Characteristics? A Retrospective Cohort Data Analysis
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Life
Volume 13
Issue 3
10.3390/life13030793
life-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Brief Report

Population Structure and Genetic Diversity Analysis of “Yufen 1” H Line Chickens Using Whole-Genome Resequencing

by
Cong Liu
1,†,
Dongxue Wang
1,†,
Yuehua He
1,
Wenjie Liang
1,
Wenting Li
1,2,3,
Kejun Wang
1,
Donghua Li
1,2,3,
Zhuanjian Li
1,2,3,
Yadong Tian
1,2,3,
Xiangtao Kang
1,2,3,* and
Guirong Sun
1,2,3,*
1
College of Animal Science and Technology, Henan Agricultural University, Zhengzhou 450001, China
2
Henan Key Laboratory for Innovation and Utilization of Chicken Germplasm Resources, Henan Agricultural University, Zhengzhou 450001, China
3
The Shennong Laboratory, Zhengzhou 450002, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Life 2023, 13(3), 793; https://doi.org/10.3390/life13030793
Submission received: 7 February 2023 / Revised: 13 March 2023 / Accepted: 13 March 2023 / Published: 15 March 2023
(This article belongs to the Section Animal Science)
Browse Figure
Versions Notes

Abstract

:
The effective protection and utilization of poultry resources depend on an accurate understanding of the genetic diversity and population structure. The breeding of the specialized poultry lineage “Yufen 1”, with its defined characteristics, was approved by the China Poultry Genetic Resource Committee in 2015. Thus, to investigate the relationship between the progenitor H line and other poultry breeds, the genetic diversity and population structure of “Yufen 1” H line (YF) were investigated and compared with those of 2 commercial chicken breeds, the ancestor breed Red Jungle Fowls, and 11 Chinese indigenous chicken breeds based on a whole-genome resequencing approach using 8,112,424 SNPs. The genetic diversity of YF was low, and the rate of linkage disequilibrium decay was significantly slower than that of the other Chinese indigenous breeds. In addition, it was shown that the YF population was strongly selected during intensive breeding and that genetic resources have been seriously threatened, which highlights the need to establish a systematic conservation strategy as well as utilization techniques to maintain genetic diversity within YF. Moreover, a principal component analysis, a neighbor-joining tree analysis, a structure analysis, and genetic differentiation indices indicated that YF harbors a distinctive genetic resource with a unique genetic structure separate from that of Chinese indigenous breeds at the genome level. The findings provide a valuable resource and the theoretical basis for the further conservation and utilization of YF.

1. Introduction

China is among the countries with the richest poultry genetic resources worldwide [1]. With the domestication of chickens, various types of poultry breeds have been formed through both natural and artificial selection, accumulating considerable genetic and phenotypic variations [2]. In this context, Chinese indigenous poultry breeds have excellent characteristics, such as early sexual maturity, a strong resistance to adversity, a rough feeding resistance, good meat quality, and a low frequency of harmful genes [3]. However, exotic breeds are widely utilized in poultry production as a way to meet market and economic needs. Although the extensive utilization and hybridization of exotic poultry breeds have improved production performance and economic efficiency, they have also negatively impacted the resources of indigenous breeds, especially those in the southeast coastal areas of China, which are on the verge of extinction [4]. In fact, it is estimated that a considerable proportion of indigenous poultry breeds (21.3% worldwide) are endangered due to breeding strategies that are deemed inadequate for the genetic sustainability of indigenous breeds [5]. Therefore, it is particularly urgent to study the genetic diversity and the relationships of indigenous poultry breeds, since this type of assessment can provide the basis for developing more effective resource conservation, development, and utilization strategies.
In recent years, with the reduction in sequencing costs, whole-genome sequencing technology has been applied to several aspects of agricultural practices, including molecular detection and breed conservation [6,7]. Numerous studies have been conducted on genetic diversity and systematic relationships based on various molecular markers at the regional scale and worldwide [8,9]. Sun J et al. [10] assessed the genetic diversity of indigenous poultry from Guangxi provinces using resequencing and found that most of the local populations could be characterized by a higher genetic diversity and lower differentiation compared to commercial chicken breeds. These authors revealed that Guangxi chickens were not significantly affected by recent inbreeding, indicating the effectiveness of breed conservation efforts for indigenous chicken populations in Guangxi. By conducting a population genetic analysis of seven Chinese indigenous chicken breeds in the context of global breeds based on 600 k SNPs, Chen L et al. [11] found that the Dongxiang blue-eggshell and Chongren Partridge (CR) breeds displayed a markedly reduced genetic diversity and a signature of admixture with European commercial breeds, which may help in establishing an efficient conservation program.
In total, 115 different chicken breeds are considered indigenous to China, which was shown to possess great advantages for genetic breeding. However, breeding methods in many parts of China still require considerable development, and several excellent poultry breeds have not yet been efficiently utilized. The matching line of the “Yufen 1” layer was authorized by the National Commission on Livestock and Poultry Genetic Resource in 2015. A closed breeding method was used to develop line H by crossing barred plumaged-original Gushi chicken (GS) with an egg-laying brown-shelled Babcock B-380 grandparent line C. As the progenitor, line H is characterized by early sexual maturity, a fast plumage line, a high egg production, and a Columbian feather pattern [12]. However, knowledge on the genomic structure of the “Yufen 1” H line chicken (YF) is still scant.
Herein, whole-genome sequencing (WGS) data were obtained from ten YF chickens from the poultry germplasm resource farm of Henan Agricultural University. By combining the published whole-genome resequencing data of 2 commercial chicken breeds, 11 Chinese indigenous chicken breeds, and 1 wild chicken breed, the genetic structure of YF representatives was thoroughly evaluated using genetic diversity and population structure analyses. The results discussed constitute scientific and effective data for the development and conservation of YF genetic resources.

2. Materials and Methods

2.1. Sample Collection and WGS Analysis

A total of 154 individuals from 15 chicken breeds were included in the present study (Supplementary Table S1). Ten blood samples of YF chickens (half males and females, unrelated per breed) were collected from the poultry germplasm resource farm of Henan Agricultural University, China. Genomic DNA was extracted from chicken blood using a TianGen DNA Kit. Paired-end libraries with a ~500 bp insert size were constructed and then subjected to sequencing using the BGISEQ-500 platform to generate 150 bp paired-end reads (BGI Genomics Co., Ltd., Cambridge, MA, USA). The genome resequencing data of the remaining 14 breeds were obtained from previously published data [13,14,15]; the breeds include Huanglang (HL), Huaixiang (HX), Huaibei-Ma (HBM), Hetian (HT), Huxu (HUXU), Jianghan (JH), Ningdu (ND), Wenchang (WC), Wuhua (WH), Yao (YAO), Gushi (GS), Red Jungle Fowl (RJF), White Leghorn (WLH), and Rhode Island Red (RIR). The accession numbers are PRJNA597842, PRJCA004227, PRJNA482210, and PRJEB30270.

2.2. Read Mapping and SNP Calling

Clean reads were aligned to the Gallus GRCg6a reference genome using the BWA-MEM alignment algorithm built into BWA (v0.7.17) [16]. Subsequent quality control processes were performed using the SortSam and MarkDuplicates in Picard Tools (v.1.56) [16], followed by variant calling, merging, and filtering using the Genome Analysis Toolkit (GATK v4.1.7.0) [17,18]. Filtering was conducted considering the criteria of QD < 2.0 || MQ < 40.0 || FS > 60.0 || SOR > 3.0 || MQRankSum < −12.5, thereby excluding reads with segregation distortions or sequencing errors. In addition, SNPs were retained using Vcftools (v0.1.16) [19] for subsequent analyses according to the following criteria: a minor allele frequency > 0.1 and a maximum miss rate > 0.8. The Ensemble genome database and SNPEff (v4.1) [20] programs were used to obtain information about SNP annotation. Additionally, the R package CMplot was used to visualize the distribution of the SNPs across the chromosome.

2.3. Genomic Diversity Analysis

To determine the genetic diversity among all chicken populations, the average minor allele frequency (MAF), observed heterozygosity (Ho), expected heterozygosity (He), and runs of homozygosity (ROHs) were estimated using PLINK (v1.90). The parameters for the ROH analysis were as follows: a sliding window of 50 SNP slides along the chromosome to estimate homozygosity; each sliding window allowed no more than 1 heterozygote, with no more than 5 missing SNPs, a minimum length of ROH of 100 kb, a minimum density of 1 SNP/50 kb, and a maximum gap between consecutive SNPs of 1000 kb. To estimate individual genomic inbreeding coefficients (FROH) using the ROH data, the length of the genome covered by ROH was divided by the total chicken autosomal genome length covered by SNPs (960,796.788 kb in the present study) [21]. The average genome-wide nucleotide polymorphisms (π) and genetic differentiation (FST) were calculated with a 10 kb sliding window and 5 kb stepwise increments using Vcftools (v0.1.16) [22]. GONE [23] was used to estimate the current generation’s effective population size for all populations. Studies have reported that the total linkage map of the chicken genome ranges from 2600 to 3800 cM [24], and the total genome size is approximately 1100 Mb [25]. A genetic distance of approximately 3 cm is equal to a physical distance of 1 Mb; hence, we used the parameter of 1 Mb = 3 cM for estimation. A linkage disequilibrium (LD) decay analysis was performed using PopLDdecay (v3.41) [26].

2.4. Population Genetic Structure and Gene Flow Analysis

Prior to performing a population genetic analysis, the SNPs were identified using LD-based pruned data with PLINK (v1.90), and the parameter was indep—pairwise 25 5 0.2. Based on the pruned SNP data, a principal component analysis (PCA) was conducted using PLINK (v1.90) to explore the genomic repartition in 15 chicken breeds. A neighbor-joining (NJ) tree of the studied populations was constructed based on pairwise distances using PHYLIP (v3.697) [27] and MEGA (v7.0) [28]. The population structure was inferred by applying the model-based clustering algorithm implemented in ADMIXTURE (v1.3.0) software [29]. TreeMix software (v1.13) [30] was used to establish a maximum likelihood tree for a gene flow analysis of the 15 chicken breeds with RJF as the outgroup. Migration events from 1 to 6 were set, and the corresponding residual matrix was generated using the options”-k 500” and “-noss”.

3. Results

3.1. SNP Calling and Annotation

The samples were sequenced at a genome coverage depth from 5.17× to 32.63×. After quality control and the filtering of unqualified SNPs, a total of 8,112,424 autosomal SNPs were obtained for statistical analyses. After annotating the obtained SNPs, most of the variations were found in introns (29.84%), intergenic regions (10.18%), and transcripts (40.41%) (Supplementary Table S2). A density distribution analysis of the SNPs on each chromosome revealed that SNPs were more frequently distributed on large chromosomes and less frequently distributed on small chromosomes, and the SNP distribution was positively correlated with chromosome length (Supplementary Figure S1). In general, the SNPs appeared to be evenly distributed on each chromosome, except for uneven distribution at the telomeres of some chromosomes.

3.2. Genetic Diversity and LD Analysis

The genetic diversity parameters are shown in Table 1. Overall, the genetic diversity indices of YF (Ho = 0.2766, He = 0.2630, MAF = 0.1932, π = 0.0024, and Ne = 70) were between those of commercial and indigenous breeds, and the indigenous breeds showed a higher genetic diversity than the commercial chicken breeds. The genetic diversity indices of RJF were lower (Ho = 0.2724, He = 0.2903, MAF = 0.2155, π = 0.0026, and Ne = 43). The highest average FROH was found in WLH (0.4427) and YF (0.3273), whereas the lowest FROH was identified in HL (0.0284). The overall level of species diversity can be determined by comparing the degree of LD decay among populations. Thus, an LD decay analysis showed that WLH had the slowest LD decay rate, followed by RIR, YF, and the other Chinese indigenous breeds. Conspicuously, the decay rate of YF was significantly slower than that of the Chinese indigenous chicken breeds and RJF (Supplementary Figure S2). To obtain greater insights into the genetic relationship between YF and the other chicken breeds, relatedness was investigated by calculating pairwise FST values (Table 2). It was revealed that the FST values of YF and the Chinese indigenous chicken breeds ranged from 0.1339 (JH) to 0.1472 (ND) with moderate genetic differentiation. In contrast, the FST values between YF and the commercial chicken breeds ranged from 0.2066 (RIR) to 0.2636 (WLH), showing a high genetic differentiation. Within the Chinese indigenous chicken populations, the lowest FST values (0.0056) were found between YAO and HX, suggesting that these chicken breeds have a considerable shared genomic background.

3.3. Population Genetic Structure and Gene Flow Analysis

PCA showed that the contribution rates of the first two principal components (PCs) to explain the variance in the data were 16.53% and 13.30%, respectively (Figure 1A). PC1 isolated WLH from all the other chicken breeds, whereas PC2 mainly contributed to the separation of YF and RIR from the other indigenous chicken breeds, WLH and RJF. The PCA results were mainly consistent with those obtained for the NJ tree constructed based on whole-genome polymorphic SNPs (Figure 1B). Of note, the NJ tree revealed the formation of 3 genetic clusters: cluster 1 comprised WLH, cluster 2 comprised YF and RIR, and cluster 3 comprised the remaining 11 Chinese indigenous chicken breeds and RJF. Furthermore, in order to determine the historical admixture patterns of the chicken populations, an ADMIXTURE analysis was conducted with K values ranging from 2 to 15. At K = 2, a genetic divergence first occurred between the commercial and noncommercial chicken breeds. When K was increased, WLH (K = 3), RIR (K = 3), RJF (K = 4), and YF (K = 5) were progressively assigned to a distinct cluster (Figure 1C). When K = 5, the cross-validation error reached the lowest value (Supplementary Figure S3 and Supplementary Table S3). To further investigate the historical split and admixture in the studied chicken breeds, a TreeMix analysis was conducted to infer the maximum likelihood population (ML) tree and potential migration events for the chicken breeds (Figure 1D,E). It was found that the inferred migration edges under the assumption of three migration events returned the lowest residuals (Supplementary Figure S4), resulting in a ML tree that best fit the data obtained in the present study. The ML tree explained 80.0% of the variance when considering three migrations. Based on this ML tree, YF and RIR were nonetheless found in a large clade. Taken together, the results from the gene flow analysis suggest that the occurrence of three migration events were as follows: RJF –> WLH, RIR –> JH, and GS –> RIR and YF.

4. Discussion

An analysis of genetic diversity and population structure using genetic data can improve or drive animal selection toward maintaining genetic diversity [31]. In this study, genetic diversity, population structure, and genetic differentiation were investigated in commercial chicken breeds, Chinese indigenous breeds, RJF, and YF populations based on whole-genome resequencing data. Furthermore, the genetic distinctiveness of YF chickens was described, which provides useful information for maintaining biodiversity [32,33].
The genetic diversity indices (Ho, He, MAF, π, and Ne) of the two commercial chicken breeds (WLH and RIR) were the lowest compared with those of the indigenous chicken breeds, indicating a lower genetic diversity among the commercial breeds. In addition, the slowest rate of LD decay was found among the two commercial chicken breeds, which could be explained by the fact that both breeds experienced intensive selection for egg production. This hypothesis can be also confirmed by the high FROH values, indicating that these two breeds underwent inbreeding to a considerably large extent [11]. Compared with the other indigenous breeds, YF had lower genetic diversity indices (Ho, He, MAF, π, and Ne), a slower rate of LD decay, and a higher degree of inbreeding. These findings are consistent with the current breeding process of YF. Line H was bred for at least twelve generations and underwent fierce selection and purification, aiming to improve individual growth, reproductive performance, and appearance. Moreover, it was observed that inbreeding occurred among YF populations; thus, the genetic diversity in YF was relatively low. These results highlight the urgent need to preserve the biodiversity in YF, and special attention should be given to the conservation of YF. RJF is the domestic chicken ancestor, and the genetic diversity within RJF should be relatively high [34]. Surprisingly, it was found that the genetic diversity in RJF was even lower than that in most indigenous breeds in this study. The observed low genetic diversity in RJF was mainly due to the RJF sampling scheme, which included closely related individuals from the same subgroup [15], with a high degree of inbreeding (FROH = 0.1696) between the sampled populations. In recent years, increased human activity has led to a decrease in the number and distribution of RJF, which resulted in a decrease in genetic diversity within this breed [35]. Egg harvesting and hunting have been considered the greatest threats to RJF populations [35,36]. Concerning the genetic differentiation indices, the highest values were found between each commercial breed and the Chinese indigenous chickens, and the results are in accordance with those reported in a previous study [1]. In particular, the level of genetic differentiation between YF and the commercial chicken breeds was high, whereas that between YF and the indigenous breeds was moderate. Collectively, these results indicate that YF has unique genetic characteristics, and it is worth further exploring its breed-specific molecular markers.
An effective population size and an appropriate selection and mating strategy are necessary components of species conservation plans. Studies have shown that the appropriate effective population size not only maintains population fitness but also facilitates species monitoring and management [37]. Indeed, through suitable mating plans, specific genetic conservation strategies might be applied to animals with the goal of increasing genetic variability and controlling inbreeding [38]. In a previous study, the genomic diversity among populations of three Chinese indigenous chicken breeds was assessed based on genome-wide SNPs, which revealed that the genetic diversity of conserved populations in which random mating prevails was high, and the inbreeding coefficient was below 0.1 [39]. This result indicated that random mating largely avoided inbreeding and can be considered a solid strategy for improving conservation schemes. Considering the current status of genetic diversity, particular attention should be given to YF, as it showed a lower genomic variability and a higher inbreeding coefficient. Therefore, an elaborate conservation program should be established for YF to avoid further inbreeding and to maintain its genetic diversity. In addition, strategies for building and utilizing genomic resource banks for conservation breeding can also greatly protect the genetic diversity of species [40].
Furthermore, the population structure and historical admixture patterns were assessed for the 15 different chicken breeds. Based on PCA and an NJ tree analysis, each individual of YF clustered together, showing a consistent genetic relationship. Interestingly, the ML tree of the gene flow analysis showed that YF and RIR clustered in a single clade, which was also confirmed by PCA and the NJ tree analysis, thereby suggesting a possible genetic relationship between these two breeds. We believe that the close genetic relationship between YF and RIR is mainly due to their low genetic diversity and drift, which further emphasizes the urgency and importance of conserving the genetic resources of YF. The migration events not only indicated the genetic contribution of the wild chicken breed RJF to the commercial breed WLH but also suggested gene exchange between the commercial chicken breed RIR and the indigenous breeds GS, YF, and JH. Finally, an ADMIXTURE analysis revealed that, at K = 5, WLH, RIR, RJF, and the other 12 indigenous chicken breeds could be distinguished from one another, and YF was clearly assigned to a separate gene pool. These results confirm the reliability of the classification of YF as a separate breed.

5. Conclusions

In conclusion, the genetic diversity and population structure of YF, two commercial chicken breeds, indigenous breeds, and the ancestor breed RJF were comprehensively studied based on genome-wide SNPs. The results show that the genetic diversity of YF was low, which may have undergone intensive selective pressure. Additionally, this study demonstrated at the whole-genome level that YF is a precious genetic resource and helped us to understand the genetic structure of YF. Based on the data described herein, conservation strategies should be implemented to maintain genetic diversity in YF, and the potential of this genetic resource should be further explored. In fact, we can control inbreeding through appropriate mating programs and develop in situ and off-site conservation programs to safeguard genetic variability and achieve specific genetic conservation plans for this breed. Moreover, the present data provide useful information about the genomics of Chinese and Western poultry breeds, which will likely assist in the development of a national project for the conservation and utilization of these breeds.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/life13030793/s1, Supplementary Figure S1: The location of SNPs along chromosomes. The x-axis represents the chromosome position (Mb), and the y-axis denotes the chromosomes. The number of SNPs present in each 1 Mb genome block is expressed via colors as indicated in the key; Supplementary Figure S2: LD decay in different populations within 300 kb; Supplementary Figure S3: The cross-validation error for K = 2 to 15; Supplementary Figure S4: Residual matrix of each migration event scenario (from 1 to 6 migration events, corresponding to A–F, respectively); the ML tree with 3 migration edges showing the smallest residuals; Supplementary Table S1: List of chicken breeds used for analyses; Supplementary Table S2: Distribution of variants in functional regions of the chicken genome; Supplementary Table S3: The cross-validation error values for K = 2 to 15.

Author Contributions

Conceptualization, Z.L. and Y.T.; data curation, K.W., D.W., Y.H., and W.L. (Wenjie Liang); methodology, W.L. (Wenting Li) and D.L.; formal analysis: C.L.; writing—original draft, C.L.; writing—review and editing, G.S. and X.K. All authors have read and agreed to the published version of the manuscript.

Funding

This project was supported by the Key Research Project of the Shennong Laboratory (SN01-2022-05), the Scientific Studio of Zhongyuan Scholars (30601985), Major Projects of Henan Province (221100110200), and the Special Fund for Genetic Breeding Post of National Laying Hens Industrial Technology System (NO. CARS-40-K04).

Institutional Review Board Statement

The animal experiments were performed according to the Guide for the Care and Use of Laboratory Animals (Ministry of Science and Technology, China, 2004). The study proposal was approved by the Institutional Animal Care and Use Committee (IACUC) of Henan Agricultural University (approval number: 11–0085).

Informed Consent Statement

Not applicable.

Data Availability Statement

All the sequence data and SNP data generated in this study were deposited in the National Genomics Data Center (https://bigd.big.ac.cn (accessed on 18 October 2022)) with the accession codes PRJCA012580 and PRJCA013041.

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. Luo, W.; Luo, C.; Wang, M.; Guo, L.; Nie, Q. Genome diversity of Chinese indigenous chicken and the selective signatures in Chinese gamecock chicken. Sci. Rep. 2020, 10, 14532. [Google Scholar] [CrossRef] [PubMed]
  2. Changsheng, N.; Pedro, A.; Yaxiong, J.; Haigang, B.; Zhonghua, N.; Lujiang, Q. Genome-Wide Single-Nucleotide Polymorphism Data Unveil Admixture of Chinese Indigenous Chicken Breeds with Commercial Breeds. Genome Biol. Evol. 2019, 11, 1847–1856. [Google Scholar]
  3. Zhang, C.; Lin, D.; Wang, Y.; Peng, D.; Li, H.; Fei, J.; Chen, K.; Yang, N.; Hu, X.; Zhao, Y.; et al. Widespread introgression in Chinese indigenous chicken breeds from commercial broiler. Evol. Appl. 2019, 12, 610–621. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Sha, Y.; Gao, C.; Liu, M.; Zhao, S. Evaluation of the genetic diversity of six Chinese indigenous chickens. Asian-Australas. J. Anim. Sci. 2020, 33, 1566–1572. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Scherf, B.D.; Fao, R.; Pilling, D. The Second Report on the State of the World’s Animal Genetic Resources for Food and Agriculture; Food and Agriculture Organization: Rome, Italy, 2015. [Google Scholar]
  6. Cheng, J.; Zhao, H.; Chen, N.; Cao, X.; Chen, H. Population structure, genetic diversity, and selective signature of Chaka sheep revealed by whole genome sequencing. BMC Genom. 2020, 21, 520. [Google Scholar] [CrossRef]
  7. Li, D.; Sun, G.; Zhang, M.; Cao, Y.; Kang, X. Breeding history and candidate genes responsible for black skin of Xichuan black-bone chicken. BMC Genom. 2020, 21, 511. [Google Scholar] [CrossRef] [PubMed]
  8. David, D.; Hubert, A.S.; Kossigan, D.G.; Chrysostome, C.C.A.A.M.; Olorounto, K.D. Genetic diversity and phylogenetic relationships within local pigs in southern Benin. Trop. Anim. Health Prod. 2021, 53, 434. [Google Scholar]
  9. Yang, Z.; Zou, L.; Sun, T.; Xu, W.; Zeng, L.; Jia, Y.; Jiang, J.; Deng, J.; Yang, X. Genome-Wide Association Study Using Whole-Genome Sequencing Identifies a Genomic Region on Chromosome 6 Associated With Comb Traits in Nandan-Yao Chicken. Front. Genet. 2021, 12, 682501. [Google Scholar] [CrossRef] [PubMed]
  10. Sun, J.; Chen, T.; Zhu, M. Whole-genome sequencing revealed genetic diversity and selection of Guangxi indigenous chickens. PLoS ONE 2022, 17, e0250392. [Google Scholar] [CrossRef]
  11. Chen, L.; Wang, X.; Cheng, D.; Chen, K.; Fan, Y.; Wu, G.; You, J.; Liu, S.; Mao, H.; Ren, J. Population genetic analyses of seven Chinese indigenous chicken breeds in a context of global breeds. Anim. Genet. 2019, 50, 82–86. [Google Scholar] [CrossRef] [Green Version]
  12. Wang, X.; Li, D.; Song, S.; Zhang, Y.; Kang, X. Combined transcriptomics and proteomics forecast analysis for potential genes regulating the Columbian plumage color in chickens. PLoS ONE 2019, 14, e0210850. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Wang, K.; Hu, H.; Tian, Y.; Li, J.; Kang, X. The chicken pan-genome reveals gene content variation and a promoter region deletion in IGF2BP1 affecting body size. Molecular Biology and Evolution. Mol. Biol. Evol. 2021, 38, 5066–5081. [Google Scholar] [CrossRef] [PubMed]
  14. Huang, X.; Otecko, N.O.; Peng, M.; Weng, Z.; Du, B. Genome-wide genetic structure and selection signatures for color in 10 traditional Chinese yellow-feathered chicken breeds. BMC Genom. 2020, 21, 316. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Qanbari, S.; Rubin, C.J.; Maqbool, K.; Weigend, S.; Weigend, A.; Geibel, J.; Kerje, S.; Wurmser, C.; Peterson, A.T.; Brisbin, I.L., Jr.; et al. Genetics of adaptation in modern chicken. PLoS Genet. 2019, 15, e1007989. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Heng, L. Toward better understanding of artifacts in variant calling from high-coverage samples. Bioinformatics 2014, 30, 2843–2851. [Google Scholar]
  17. Heng, L.; Richard, D. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 2009, 25, 1754–1760. [Google Scholar]
  18. Lin, Y.L.; Chang, P.C.; Hsu, C.; Hung, M.Z.; Chien, Y.H.; Hwu, W.L.; Lai, F.; Lee, N.C. Comparison of GATK and DeepVariant by trio sequencing. Sci. Rep. 2022, 12, 1809. [Google Scholar] [CrossRef] [PubMed]
  19. Bayer, P.E.; Gill, M.; Danilevicz, M.F.; Edwards, D. Producing High-Quality Single Nucleotide Polymorphism Data for Genome-Wide Association Studies. Methods Mol. Biol. 2022, 2481, 153–159. [Google Scholar] [PubMed]
  20. Cingolani, P.; Platts, A.; Wang, L.L.; Coon, M.; Nguyen, T.; Wang, L.; Land, S.J.; Lu, X. Ruden DM: A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff. Fly 2012, 6, 80–92. [Google Scholar] [CrossRef] [Green Version]
  21. Mcquillan, R.; Leutenegger, A.L.; Abdelrahman, R.; Franklin, C.S.; Pericic, M.; Baraclauc, L.; Smolejnarancic, N.; Janicijevic, B.; Polasek, O.; Tenesa, A. Runs of homozygosity in European populations. Am. J. Hum. Genet. 2008, 83, 359–372. [Google Scholar] [CrossRef] [Green Version]
  22. Danecek, P.; Auton, A.; Abecasis, G.; Albers, C.A.; Banks, E.; Depristo, M.A.; Handsaker, R.E.; Lunter, G.; Marth, G.T.; Sherry, S.T. The variant call format and VCFtools. Bioinformatics 2011, 27, 2156–2158. [Google Scholar] [CrossRef] [PubMed]
  23. Santiago, E.; Novo, I.; Pardiñas, A.F.; Saura, M.; Wang, J.; Caballero, A. Recent demographic history inferred by high-resolution analysis of linkage disequilibrium. Mol. Biol. Evol. 2020, 37, 3642–3653. [Google Scholar] [CrossRef] [PubMed]
  24. Groenen, M.A.; Cheng, H.H.; Bumstead, N.; Benkel, B.F.; Briles, W.E.; Burke, T.; Burt, D.W.; Crittenden, L.B.; Dodgson, J.; Hillel, J. A consensus linkage map of the chicken genome. Genome Res. 2000, 10, 137–147. [Google Scholar] [PubMed]
  25. Warren, W.C.; Hillier, L.W.; Tomlinson, C.; Minx, P.; Kremitzki, M.; Graves, T.; Markovic, C.; Bouk, N.; Pruitt, K.D.; Thibaud-Nissen, F. A new chicken genome assembly provides insight into avian genome structure. G3 Genes Genomes Genet. 2017, 7, 109–117. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Zhang, C.; Dong, S.S.; Xu, J.Y.; He, W.M.; Yang, T.L. PopLDdecay: A fast and effective tool for linkage disequilibrium decay analysis based on variant call format files. Bioinformatics 2019, 35, 1786–1788. [Google Scholar] [CrossRef]
  27. Felsenstein, J. PHYLIP-Phylogeny Inference Package (Version 3.2). Cladistics 1989, 5, 164–166. [Google Scholar]
  28. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [Green Version]
  29. Alexander, D.H.; Lange, K. Enhancements to the ADMIXTURE algorithm for individual ancestry estimation. BMC Bioinform. 2011, 12, 246. [Google Scholar] [CrossRef] [Green Version]
  30. Pickrell, J.K.; Pritchard, J.K.; Tang, H. Inference of population splits and mixtures from genome-wide allele frequency data. PLoS Genet. 2012, 8, e1002967. [Google Scholar] [CrossRef] [Green Version]
  31. Wang, Y.; Dong, R.; Li, X.; Cui, C.; Yu, G. Analysis of the Genetic Diversity and Family Structure of the Licha Black Pig Population on Jiaodong Peninsula, Shandong Province, China. Animals 2022, 12, 1045. [Google Scholar] [CrossRef]
  32. Strillacci, M.G.; Cozzi, M.C.; Gorla, E.; Mosca, F.; Schiavini, F.; Román-Ponce, S.; López, F.R.; Schiavone, A.; Marzoni, M.; Cerolini, S. Genomic and genetic variability of six chicken populations using single nucleotide polymorphism and copy number variants as markers. Animal 2017, 11, 737–745. [Google Scholar] [CrossRef] [PubMed]
  33. Cendron, F.; Perini, F.; Mastrangelo, S.; Tolone, M.; Cassandro, M. Genome-Wide SNP Analysis Reveals the Population Structure and the Conservation Status of 23 Italian Chicken Breeds. Animals 2020, 10, 1441. [Google Scholar] [CrossRef] [PubMed]
  34. Zhang, J.; Nie, C.; Li, X.; Ning, Z.; Qu, L. Genome-Wide Population Genetic Analysis of Commercial, Indigenous, Game, and Wild Chickens Using 600K SNP Microarray Data. Front. Genet. 2020, 11, 543294. [Google Scholar] [CrossRef] [PubMed]
  35. Setianto, J.; Zain, B.; Sutriyono; Prakoso, H. Domestication of red jungle fowl: A case study of the red jungle fowl chicks procurement by the communities in Central Bengkulu, Indonesia. Biodiversitas 2017, 18, 183–189. [Google Scholar] [CrossRef]
  36. Akrim, F.; Awan, M.S.; Mahmood, T.; Anjum, M.Z.; Qasim, S.; Khalid, J.; Andleeb, S. Threats to Red Junglefowl (Gallus gallus murghi) in Deva Vatala National Park, District Bhimber, Azad Jammu and Kashmir, Pakistan. Sciencedomain 2015, 6, 59–65. [Google Scholar] [CrossRef]
  37. Meuwissen, T.; Woolliams, J.A. Effective sizes of livestock populations to prevent a decline in fitness. Theor. Appl. Genet. 1994, 89, 1019–1026. [Google Scholar] [CrossRef] [PubMed]
  38. Zanetti, E.; Marchi, M.D.; Dalvit, C.; Cassandro, M. Genetic characterization of local Italian breeds of chickens undergoing in situ conservation. Poult. Sci. 2010, 89, 420–427. [Google Scholar] [CrossRef]
  39. Mengmeng, Z.; Wei, H.; Hui, T.; Guohui, L.; Minjie, Z.; Ran, X.; Yijun, L.; Tao, Y.; Wenting, L.; Jianmin, Z. Genomic diversity dynamics in conserved chicken populations are revealed by genome-wide SNPs. BMC Genom. 2018, 19, 598. [Google Scholar]
  40. Ballou, J.D.; Lacy, R.C.; Traylor-Holzer, K.; Bauman, K.; Ivy, J.A.; Asa, C. Strategies for establishing and using genome resource banks to protect genetic diversity in conservation breeding programs. Zoo Biol. 2022. Online ahead of print. [Google Scholar] [CrossRef]
Life 13 00793 g001 550
Figure 1. Population structure and gene flow analysis of all chicken individuals. (A) First and second principal components from a principal component analysis of all populations; (B) neighbor-joining tree for all individual chickens. (C) ADMIXTURE analysis for K = 2 up to 6. The same color indicates groups with the same ancestry. (D) Maximum likelihood tree with 3 migration events. Migration events are shown as colored arrows, shaded according to their weight. Horizontal branch lengths are proportional to the amount of genetic drift that has occurred at each branch. The scale bar shows 10 times the average standard error of the entries in the sample covariance matrix; (E) residual fit from the maximum likelihood tree in (D).
Figure 1. Population structure and gene flow analysis of all chicken individuals. (A) First and second principal components from a principal component analysis of all populations; (B) neighbor-joining tree for all individual chickens. (C) ADMIXTURE analysis for K = 2 up to 6. The same color indicates groups with the same ancestry. (D) Maximum likelihood tree with 3 migration events. Migration events are shown as colored arrows, shaded according to their weight. Horizontal branch lengths are proportional to the amount of genetic drift that has occurred at each branch. The scale bar shows 10 times the average standard error of the entries in the sample covariance matrix; (E) residual fit from the maximum likelihood tree in (D).
Life 13 00793 g001
Table
Table 1. Genetic diversity indices.
Table 1. Genetic diversity indices.
BreedsN 1Breed AssignmentHo 2He 3FROH 4Π 5Ne 6MAF 7
WLH 810Commercial breeds (Italy)0.1332 0.1757 0.4427 0.0016 110.1279
RIR 910Commercial breeds (America)0.1915 0.2128 0.3096 0.0019 590.1582
GS 1010Native breeds (Henan province, China)0.3129 0.3206 0.0793 0.0028 1200.2372
HBM 1110Native breeds (Anhui province, China)0.2992 0.3056 0.0469 0.0027 1670.2273
HL 1210Native breeds (Hunan province, China)0.3042 0.3174 0.0284 0.0028 4970.2353
HT 1310Native breeds (Fujian province, China)0.3019 0.3095 0.0472 0.0028 3060.2301
HUXU 149Native breeds (Guangdong province, China)0.2794 0.3072 0.0424 0.0027 1830.2283
HX 1510Native breeds (Guangdong province, China)0.2991 0.3133 0.0405 0.0028 2620.2325
JH 1610Native breeds (Hubei province, China)0.2817 0.3064 0.0784 0.0027 520.2278
ND 1710Native breeds (Jiangxi province, China)0.2891 0.3055 0.0540 0.0027 2620.2269
WC 1810Native breeds (Hainan province, China)0.2936 0.3134 0.0374 0.0028 2160.2324
WH 1910Native breeds (Anhui province, China)0.3033 0.3078 0.0352 0.0027 520.2288
YAO 2010Native breeds (Guangxi province, China)0.3026 0.3136 0.0422 0.0028 2640.2328
RJF 2115Wild breeds (Thailand)0.2724 0.2903 0.1696 0.0026 430.2155
YF 2210Native breeds (Henan province, China)0.2766 0.2630 0.3273 0.0024 700.1932
1 N = the number of animals per breed, 2 Ho = observed heterozygosity, 3 He = expected heterozygosity, 4 FROH = inbreeding coefficients, 5 π = nucleotide polymorphisms, 6 Ne = effective population size, 7 MAF = minor allele frequency, 8 WLH = White Leghorn, 9 RIR = Rhode Island Red, 10 GS = Gushi, 11 HBM = Huaibei-Ma, 12 HL = Huanglang, 13 HT = Hetian, 14 HUXU = Huxu, 15 HX = Huaixiang, 16 JH = Jianghan, 17 ND = Ningdu, 18 WC = Wenchang, 19 WH = Wuhua, 20 YAO = Yao, 21 RJF = Red Jungle Fowl, 22 YF = “Yufen 1” H line.
Table
Table 2. Genetic differentiation index between populations.
Table 2. Genetic differentiation index between populations.
BreedsWLH 1RIR 2GS 3HBM 4HL 5HT 6HUXU 7HX 8JH 9ND 10WC 11WH 12YAO 13RJF 14YF 15
WLH 10 0.3006 0.2232 0.2239 0.2039 0.2123 0.2043 0.2010 0.2079 0.2164 0.1989 0.2092 0.2020 0.2127 0.2636
RIR 2 0 0.1648 0.1648 0.1580 0.1650 0.1619 0.1584 0.1480 0.1712 0.1593 0.1654 0.1589 0.1947 0.2066
GS 3 0 0.0278 0.0272 0.0429 0.0460 0.0437 0.0405 0.0424 0.0428 0.0492 0.0437 0.0947 0.1404
HBM 4 0 0.0193 0.0365 0.0405 0.0386 0.0333 0.0350 0.0371 0.0438 0.0390 0.0921 0.1411
HL 5 0 0.0166 0.0193 0.0149 0.0262 0.0151 0.0135 0.0240 0.0163 0.0683 0.1347
HT 6 0 0.0307 0.0278 0.0417 0.0254 0.0244 0.0355 0.0288 0.0761 0.1427
HUXU 7 0 0.0205 0.0398 0.0320 0.0170 0.0293 0.0221 0.0746 0.1372
HX 8 0 0.0390 0.0274 0.0166 0.0272 0.0056 0.0696 0.1368
JH 9 0 0.0410 0.0375 0.0452 0.0393 0.0903 0.1339
ND 10 0 0.0253 0.0349 0.0291 0.0757 0.1472
WC 11 0 0.0249 0.0179 0.0650 0.1378
WH 12 0 0.0280 0.0781 0.1432
YAO 13 0 0.0706 0.1372
RJF 14 0 0.1780
YF 15 0
1 WLH = White Leghorn, 2 RIR = Rhode Island Red, 3 GS = Gushi, 4 HBM = Huaibei-Ma, 5 HL = Huanglang, 6 HT = Hetian, 7 HUXU = Huxu, 8 HX = Huaixiang, 9 JH = Jianghan, 10 ND = Ningdu, 11 WC = Wenchang, 12 WH = Wuhua, 13 YAO = Yao, 14 RJF = Red Jungle Fowl, 15 YF = “Yufen 1” H line.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Liu, C.; Wang, D.; He, Y.; Liang, W.; Li, W.; Wang, K.; Li, D.; Li, Z.; Tian, Y.; Kang, X.; et al. Population Structure and Genetic Diversity Analysis of “Yufen 1” H Line Chickens Using Whole-Genome Resequencing. Life 2023, 13, 793. https://doi.org/10.3390/life13030793

AMA Style

Liu C, Wang D, He Y, Liang W, Li W, Wang K, Li D, Li Z, Tian Y, Kang X, et al. Population Structure and Genetic Diversity Analysis of “Yufen 1” H Line Chickens Using Whole-Genome Resequencing. Life. 2023; 13(3):793. https://doi.org/10.3390/life13030793

Chicago/Turabian Style

Liu, Cong, Dongxue Wang, Yuehua He, Wenjie Liang, Wenting Li, Kejun Wang, Donghua Li, Zhuanjian Li, Yadong Tian, Xiangtao Kang, and et al. 2023. "Population Structure and Genetic Diversity Analysis of “Yufen 1” H Line Chickens Using Whole-Genome Resequencing" Life 13, no. 3: 793. https://doi.org/10.3390/life13030793

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop