Next Article in Journal
The Low-Copy-Number Satellite DNAs of the Model Beetle Tribolium castaneum
Previous Article in Journal
Mimicked Mixing-Induced Heterogeneities of Industrial Bioreactors Stimulate Long-Lasting Adaption Programs in Ethanol-Producing Yeasts
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 14
Issue 5
10.3390/genes14050998
genes-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:
Article

Studying Genetic Diversity and Relationships between Mountainous Meihua Chickens Using Mitochondrial DNA Control Region

by
Bo Ran
1,2,
Wei Zhu
2,
Xiaoling Zhao
2,
Linxiang Li
3,
Zhixin Yi
3,
Miao Li
3,
Tao Wang
1 and
Diyan Li
1,*
1
School of Pharmacy, Chengdu University, Chengdu 610106, China
2
College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 611130, China
3
Bazhong Academy of Agricultural and Forestry Sciences, Bazhong 610066, China
*
Author to whom correspondence should be addressed.
Genes 2023, 14(5), 998; https://0-doi-org.brum.beds.ac.uk/10.3390/genes14050998
Submission received: 30 March 2023 / Revised: 18 April 2023 / Accepted: 26 April 2023 / Published: 28 April 2023
(This article belongs to the Section Animal Genetics and Genomics)
Browse Figures
Review Reports Versions Notes

Abstract

:
The Mountainous Meihua chicken is a unique regional germplasm resource from Tongjiang County, Bazhong City, China, but its genetic structure and evolutionary relationships with other native chicken breeds in the Sichuan region remain unclear. Here, we analyzed a total of 469 sequences, including 199 Mountainous Meihua chicken sequences generated in this study, together with 30 sequences representing 13 clades and 240 sequences from seven different Sichuan local chicken breeds downloaded from NCBI. These sequences were further used to analyze genetic diversity, patterns of population differentiation, and phylogenetic relationships between groups. We show that Mountainous Meihua chicken mtDNA sequences have high haplotypic and nucleotide diversity (0.876 and 0.012, respectively) and with a T bias that is suggestive of good breeding potential. Phylogenetic analysis showed that Mountainous Meihua chickens belong to clades A, B, E, and G and have a low affinity to other chicken breeds, with a moderate degree of differentiation. A non-significant Tajima’s D indicates that no demographic expansions occurred in the past. Finally, the four maternal lineages identified in Mountainous Meihua chicken showed unique genetic characteristics.

1. Introduction

The Bazhong Mountainous Meihua chicken (Galus galus domesticus) is a local genetic resource with ~2000 years of history in the Qinba mountain region [1]. In 2013, a general survey of germplasm resources showed that Mountainous Meihua chickens are primarily distributed in Nanjiang County, Bazhong City, in particular across 10 townships (e.g., Shaba, Huitan, and Hongyan of Nanjiang County) [2]. The unique shape, moderate growth speed, and excellent egg and meat quality make Mountainous Meihua chickens a highly valuable resource [3]. The breed is also resistant to rough feeding and possesses a good foraging ability, adaptability, and disease resistance [4], all of which are important traits for improving local breeds. In recent years, and thanks to the support of local municipal governments, seed conservation farms meant for planned seed selection and breeding have been established [5]. Now, there is one protected population, 68 protected families, and more than 2000 pure chickens of this breed.
As it is an important extracellular genetic information carrier, mitochondrial DNA (mtDNA) has low relative molecular weights, high copy numbers, maternal inheritance, and no/rare recombination events, which means that it is often used as a marker in many molecular studies [6,7,8,9]. It is a popular choice for phylogenetic, population genetic, and phylogeographic analyses [10]. The mtDNA D-loop region is a non-coding region rich in A/T bases, which can be recognized by specific molecules such as RNA polymerase and used as the starting sequence of transcription [11]. In the mitochondrial genome, coding genes are mainly used for the biodiversity research of different species, while non-coding D-loop regions are commonly used for genetic diversity analysis and maternal origin research due to their evolution speed being 5–10 times higher than that of other regions on mtDNA [12,13,14].
Using mtDNA information, previous studies indicated that domestic chickens independently originated from a group of red jungle poultry found in forests in Southeast Asia and India and gradually distributed throughout the world with continuous domestication and transmission, resulting in many new breeds [15]. Research on the origin of Chinese domestic chicken breeds revealed that the red jungle fowl from Thailand and its neighboring areas are considered to be the maternal ancestors of Chinese chicken breeds [16]. However, a series of independent domestication events occurred across Asia and the Indian subcontinent, supporting the theory of there being multiple origins of domestic chickens [17,18]. In 2006, Liu et al. [19] first proposed that domestic chickens in Eurasia are distributed among nine highly different mtDNA evolutionary branches (A–I) and pointed out different regions where different branches may originate. Later, other researchers proposed more branches and further determined the maternal origin of domestic chickens [20]. Analysis of the origin of Japanese indigenous chickens revealed that the wild red jungle fowl was domesticated in Southeast Asia and passed into China to form non-wild clades (A and B), which together form the foundation of Japanese indigenous chickens [21]. The research on the origins and transmission routes of domestic chickens in the Middle East and Africa supports the hypothesis that the Indian subcontinent is the maternal origin of these chickens. Additionally, ancient and modern trade networks contributed to the modern diversity of native chickens [22]. A growing number of studies on mtDNA D-loop regions have revealed genetic diversity within populations and inter-population relationships of some new local chicken breeds in various regions in recent years [23,24,25,26], which have provided important support for modern breeding or chicken genetic diversity protection plans.
Mountainous Meihua chickens have been genetically characterized in the past, in particular for phenotypes associated with growth and development and muscle and egg quality [3]. However, the evolutionary mechanisms of differentiation for this chicken breed have yet not been fully explored. Molecular technology can efficiently and accurately study the genetic differences among breeds of animals. Mitochondrial DNA (mtDNA) polymorphisms in various animal populations can help to clarify ancestral lineages or compare different populations, which has also been shown in chickens [27]. Here, we analyzed the partial mtDNA control region sequences of Mountainous Meihua chickens and seven other chicken breeds sampled in the Sichuan Province of China, with the goal to assess the levels of genetic diversity and population structure to uncover ancestral maternal lineages. Our results provide insights into the evolutionary history of Meihua chickens and other local breeds, which can help to improve this valuable genetic resource.

2. Materials and Methods

2.1. Ethical Statement

The collection of blood samples from Mountainous Meihua chickens was performed according to the Guidelines for Experimental Animals established by the Ministry of Science and Technology (Beijing, China, revised in March 2017). Ethical approval was provided by the Ethics Committee of the Sichuan Agricultural University (protocol number 2020202010).

2.2. Sample Collection and DNA Extraction

Blood samples were collected from 199 Mountainous Meihua chickens via wing veins. These chickens were all adult males in the Mountainous Meihua chickens breeding base in the Shuangquan Township, Tongjiang County. The blood samples were stored at −20 degrees to ensure the integrity of genomic DNA. These samples were then sent to the College of Animal Science and Technology at Chengdu Campus of Sichuan Agricultural University for genomic DNA and further genetic analysis. A proteinase-K/Phenol-chloroform standard procedure was used to extract genomic DNA, as previously described [28]. Isolated DNA samples were measured in terms of concentration and assessed in terms of quality using a NanoDrop 2000 ultramicroscopic spectrophotometer.

2.3. PCR Amplification and DNA Sequencing

The following forward (5′-AGGACTACGGCTTGAAAAGC-3′) [29] and reverse (5′-ATGTGCCTGACCGAGGAACCAG-3′) [19] primers were used for amplification of the mtDNA control region (approximately 500 bp) of Mountainous Meihua chicken. PCR amplification was performed using a 25 μL solution containing 12.5 μL of 2 × Taq PCR Master Mix, 1.25 μL of template genomic DNA (200 ng), 1.25 μL of each primer (10 pmol/μL), and 8.75 μL of ddH2O. After incubation at 95 °C for 3 min, 33 cycles were performed, as follows: 94 °C for 35 s, 58.4 °C for 35 s, and 72 °C for 1 min, with a final extension at 72 °C for 10 min. PCR products were purified using the gel recovery kit and sequenced using an ABI 3730 automated sequencer (Applied Biosystems, Carlsbad, CA, USA).

2.4. Data Analysis

The Nucleotide Basic Local Alignment Search Tool (BLASTN) was used to ensure the accuracy of the sequenced DNA fragments. The mitochondrial D-loop sequences of seven other Sichuan chicken breeds [30] were downloaded from NCBI (GenBank accession numbers: N381367 and MN381709), including 27 sequences for Pengxian yellow chickens (PX), 29 sequences for Jinyang silky chickens (JS), 30 sequences for Emei Black chickens (EM), 40 sequences for Jiuyuan black chickens (JY), 29 sequences for Muchuan silky chickens (MC), 39 sequences for Miyi chickens (MY), and 46 sequences for Shimian Caoke chickens (SM). The geographical distribution of the 8 chicken breeds is shown in Figure 1. The obtained 439 sequences were aligned using ClustalX 2.1 in the DNASTAR software package and truncated to the longest common sequence (482 bp) among all samples for further analysis. Taking the red jungle fowl (NC_007235) as the reference sequence, we analyzed genetic variation using DnaSP vs6.0 and defined haplotypes, extracted variation site information (V), and calculated the haplotype diversity (Hd), nucleotide diversity (π), the average number of nucleotide differences (K), Tajima’s D values, Nucleotide divergence (Dxy), net genetic distance (Da), Genetic differentiation index (Gst) and gene flow (Nst), haplotypes, and average no. of nucleotide differences. To evaluate sequence variation among and within populations, we conducted analysis of molecular variance (AMOVA) and pairwise fixation index (Fst) using Arlequin version 3.5 [31]. MEGA 6.0 [32] was used to analyze sequence base composition and establish phylogenetic relationships among the 8 chicken breeds. A total of 30 mitochondrial D-loop sequences downloaded from NCBI, representing a total of 13 clades, were used to analyze the maternal origin of the populations, following previous studies [20]. An unrooted neighbor-joining (NJ) tree was constructed using the Kimura 2-parameter model. Bootstrap values of the phylogenetic tree were estimated from 1000 replicates. A maximum parsimony median-joining (MJ) network was made of the 439 sequences using the software PopART [33].

3. Results

3.1. Sequence Variation and Population Diversity

The mtDNA of 439 individuals from eight chicken breeds was truncated to 482 bp. The observed base compositions of the mtDNA D-Loop in MH were 30.51% (T), 28.09% (C), 28.38% (A), and 13.02% (G) (Table 1), showing an obvious A/T bias. Other chicken breeds showed similar base composition signatures.
We identified 46 polymorphic sites, of which 39 were parsimony informative ones, and seven were singleton variable sites. These polymorphisms defined a total of 44 haplotypes (Figure 2), ranging from 5 (MC) to 23 (MH) in each breed. We also found 17 haplotypes shared among multiple populations (Table 2), in particular MY, SM, and MC. The haplotypes found in MH mainly include H2, H3, H5, H7, H10, H12, and H20, among which H3, H5, H7, H12, and H20 were shared haplotypes, accounting for 53.77% of the total number of MH samples analyzed. In contrast, H2 and H10 were unique MH haplotypes and accounted for 31.16% of the total samples.
Subsequently, we analyzed the overall genetic diversity indices in MH (Table 3) and contrasted them with the corresponding information available for other chicken breeds [30]. A total of 28 polymorphic sites were present in MH, of which 25 were parsimony informative sites and 3 were singletons. MH haplotype (0.876) and nucleotide (0.012) diversity, as well as the average number of nucleotide differences (5.936), were in the upper range of all chicken breeds studied. The estimated Tajima’s D across all breeds was not significant, indicating that no significant demographic expansions occurred in the past.
The genetic variation within and between population was quantified via AMOVA based on Kimura-2 parameter distances (Table 4) and showed that most genetic diversity was found within breeds (81.50%). Notwithstanding, we found that the differences between populations were highly significant (p < 0.01).

3.2. Nucleotide Divergence and Genetic Differentiation between Populations

We analyzed the average number of nucleotide substitutions per site (Dxy) and the net divergence (Da) between populations based on the mtDNA control region (Table 5). We estimated that the Dxy was between the eight chicken breeds ranged between 0.930 (PX) and 1.710 (JS), while the maximum Da was between MC and JY (0.735), and the minimum degree of differentiation was between JS and PX (−0.027). The Dxy between MH and other chicken breeds ranged from 1.212 to 1.652, with the Da ranging from 0.163 to 0.578.
We also estimated the coefficient of genetic differentiation (Gst) and gene flow (Nst) between eight populations (Table 5) and found the maximum degree of differentiation was between SM and PX and minimum degree of differentiation was between JS and PX. Furthermore, MH was the greatest genetic distance away from MC (Nst = 37.888) and had the most affinity with JS (Nst = 11.728).

3.3. Phylogenetic Relationships and Haplotype Network Using the mtDNA D-Loop

Phylogenetic analysis was conducted to estimate evolutionary genetic distances using the maximum composite likelihood method available in the MEGA6.0 software. Phylogenetic analysis (Figure 3) revealed 44 haplotypes in indigenous chickens that were divided across 5 (A, B, C, E, and G) of the 13 clades (clades A–I and W–Z), following the classification of Liu [19] and Miao [20]. The five clades, A, B, C, E, and G, contained 21, 6, 3, 9, and 5 haplotypes, respectively. MH chicken haplotypes were distributed across the A, B, E, and G clades. We further analyzed the haplotype composition in MH chicken breeds (Table 6), which showed that 29.65% of the samples belong to the A clade, 41.71% belong to the B clade, 27.14% belong to the E clade, and 1.50% belong to the G clade.
A maximum parsimony median-joining network was constructed based on the 44 haplotypes to further determine the phylogenetic relationship and genetic distances among the eight native chicken breeds. Network analysis depicts the minimum mutational distances among the haplotypes, and these are consistent with features observed in the neighbor-joining tree (Figure 4). By contrast, the haplotypes present within clades are closer to each other than those identified between clades are. Clades A and E are shared by all chicken breeds.

4. Discussion

Here, we presented an examination of the genetic diversity and maternal origin of Mountainous Meihua chickens (MH) and explored the genetic differentiation and evolutionary relationships between MH and seven other Sichuan chicken breeds. The analysis of the base composition of the mitochondrial D-loop indicated an A/T bias across all eight chicken breeds, which is consistent with studies in some other chicken breeds [25,34]. Most previous studies found that variable sites in the D-loop region were in 133–446 bp, as the mutation rate in this region is higher than it is in the other regions [35]. In our study, most of the genetic variation among the eight breeds was also located in this region. However, some variation sites found by Yin [36] were located at positions 711, 1214, 1215, and 1222 in the Lueyang Black-bone Chicken, and at positions 686, 792, and 1215 in the Jiangxi Anyi Tile-like Chicken [37]. As a result, longer sequences are likely to provide more comprehensive information on variation sites than shorter sequences will.
Studies have shown that commercial breeds have relatively low genetic diversity, as only a few sires produce many offspring [38]. Native chicken breeds, however, become good genetic resources for future breeding needs in order to create better varieties for commercial use due to their high level of genetic variation [39,40]. The ability of species to adapt to harsh environments and resist catching diseases can be improved by the genetic diversity within them [35]. Haplotype and nucleotide diversities are two key indicators to evaluate genetic variation, with higher values being indicative of higher genetic diversity [41]. In MH, we found haplotype and nucleotide diversities of 0.876 and 0.012, respectively, which are higher than those of most of local breeds here analyzed; it is important to note that the sample sizes of different breeds in our study are inconsistent. In order to properly compare the samples with different sizes, those conducting future studies should use the rarefaction method. Similar results were previously found when comparing MH with other indigenous chicken breeds from China using the same metrics [42,43]. These results could provide further support for improving breeding practices.
Tajima’s D values were not significant across all chicken breeds, indicating that neither balancing nor purifying selection occurred in these chicken breeds [30], and thus, they evolved neutrally. This is consistent with observations of native flocks in Thailand and the Philippines [17,44]. MH were raised as free-range backyard flocks for a long time, with episodes of random mating, and little to no control of breeding practices. These factors may have reduced artificial selective pressures, and thus, support neutral evolution. Further hierarchical analysis of molecular variance, AMOVA, indicated that the differentiation among maternal lineages in the eight chicken populations mainly resulted from variation within populations, which is similar to previous findings by Wani [7] and Teinlek [44].
Gst and Nst are commonly used to characterize genetic differentiation and gene exchanges between populations. Here, we found the highest degree of differentiation between MH and JY, and the lowest one between MH and MY. Similarly, while gene exchanges between MH and MC were frequent, those between MH and JS were rare. Our results showed large differences between MH and MY, and little diversification between MH and JY. In addition, we observed the highest Da value between MH and MC, and the lowest one between MH and EM. These differences may result from different living conditions, selection pressures, geographical isolation, and breeding histories [45,46,47,48,49,50].
According to Liu [19] and Miao [20], the mtDNA phylogenetic tree of domestic chicken and the red jungle fowl can be divided into 13 clades. Of these, clades A and B are widely distributed across the world and native to Yunnan and the surrounding areas. Clade C is endemic to the Huang-Huai River Basin in China, clade G is distributed in Yunnan and the adjacent areas, and clade E is common among Indian and European domestic chickens. A study on the matrilineal origin of Indonesian chickens showed that Indonesian chickens were distributed among haplogroups B, D, and E [51]. The phylogenetic results of a study of Bangladesh native chickens showed five distinct mitochondrial haplogroups (A, D, E, F, and I), suggesting that they have a more maternal origin [10]. A previous study indicated that the mtDNA D-loop sequences of seven Sichuan native chickens were clustered in clade B [30], but we found these groups include five clades, specifically A, B, C, E, and G. Clades A, B, and E were the dominant lineages in the MH group, with a frequency of 98.5%, with clade E representing 1.5%, which could have resulted from admixture with local Yunnan varieties. Interestingly, this clade is represented by two haplotypes. One is the most frequent clade, G haplotype, and the other one (detected only in MH) diverges from it, and it can be identified on the basis of one substitution. This could be indicative that this admixture happened a long time ago since the new haplotype detected only in MH had time to evolve. Unlike other chicken populations, a considerable part of the MH population (41.71%) belong to haplogroup B, showing characteristics that are different from those of other chicken breeds. In addition, maternal lineage sharing has been reported among different indigenous chicken breeds in various geographical locations [52,53,54]. Among the eight flocks studied, clades A and E were found to be in common, indicating a common origin.
In summary, our study suggests that Mountainous Meihua chickens have a relatively high level of genetic diversity and present a certain degree of genetic differentiation from seven other native breeds from China. The results also showed that there were four maternal lineages in this breed, specifically belonging to clades A, B, and E. Further genomic analyses should be conducted in the future to provide further understanding of this unique breed.

Author Contributions

Conceptualization, D.L. and T.W.; methodology, B.R.; formal analysis, W.Z.; investigation, W.Z. and M.L.; writing—original draft preparation, B.R., L.L., Z.Y. and W.Z.; writing—review and editing, X.Z.; visualization, L.L.; supervision, X.Z.; project administration, D.L. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the Sichuan Provincial Science and Technology Programs (Grant No. 2021YFN0112).

Institutional Review Board Statement

The collection of blood samples from Mountainous Meihua chickens was performed according to the Guidelines for Experimental Animals established by the Ministry of Science and Technology (Beijing, China, revised in March 2017). Ethical approval was pro-vided by the Ethics Committee of the Sichuan Agricultural University (protocol number 2020202010).

Informed Consent Statement

Not applicable.

Data Availability Statement

The mitochondrial DNA control region data were submitted to NCBI with the GenBank accession numbers OP455159–OP455357.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, Y.; Xu, J.-M.; Li, Q. High-quality local chicken breed-Qinba chicken. Sichuan Anim. Vet. Sci. 2016, 43, 40–41. [Google Scholar]
  2. Zhang, B.-W.; Li, J.; Ai, Z.-Y. The Protection Mechanism of the Geographical Indication Brand of Bashan native Chickens. Sichuan Anim. Vet. Sci. 2019, 46, 21–22. [Google Scholar]
  3. Deng, X.; Gui, T.-T.; Shu, G.; Liu, D.-S.; Li, M.; Pu, Z.-Q.; Zhang, G.; Sun, D.-F.; Tian, Y.-F.; Zhu, Q.; et al. Study on Genetic Characteristics of “Bazhong Mountainous Meihua Chicken”. China Poult. 2021, 43, 112–116. [Google Scholar]
  4. Du, H.; Sun, D.-F.; Yi, X.; Xu, D.; Liao, J.-H.; Wang, F.; Shu, G. Effects of Probiotics Combined with Fermentation Bed Culture Model on Growth Performance, Slaughter Performance, Immunologic Function and Intestinal Flora of Bazhong Mountain Plum Chicken. Feed. Ind. 2021, 42, 6–10. [Google Scholar]
  5. Zhang, Y.; Xu, J.-M.; Li, Q. Protection and Development of Qinba Chicken Breed Resources. China Anim. Ind. 2016, 16, 54–55. [Google Scholar]
  6. Lan, D.; Hu, Y.; Zhu, Q.; Liu, Y. Mitochondrial DNA study in domestic chicken. Mitochondrial DNA Part A DNA Mapp. Seq. Anal. 2017, 28, 25–29. [Google Scholar] [CrossRef]
  7. Wani, C.E.; Yousif, I.A.; Ibrahim, M.E.; Musa, H.H. Molecular Characterization of Sudanese and Southern Sudanese Chicken Breeds Using mtDNA D-Loop. Genet. Res. Int. 2014, 2014, 928420. [Google Scholar] [CrossRef] [PubMed]
  8. Quinn, T.W.; Wilson, A.C. Sequence evolution in and around the mitochondrial control region in birds. J. Mol. Evol. 1993, 37, 417–425. [Google Scholar] [CrossRef]
  9. Ishida, N.; Hasegawa, T.; Takeda, K.; Sakagami, M.; Onishi, A.; Inumaru, S.; Komatsu, M.; Mukoyama, H. Polymorphic sequence in the D-loop region of equine mitochondrial DNA. Anim. Genet. 1994, 25, 215–221. [Google Scholar] [CrossRef]
  10. Bhuiyan, M.S.; Chen, S.; Faruque, S.; Bhuiyan, A.K.; Beja-Pereira, A. Genetic diversity and maternal origin of Bangladeshi chicken. Mol. Biol. Rep. 2013, 40, 4123–4128. [Google Scholar] [CrossRef]
  11. Randi, E.; Lucchini, V. Organization and evolution of the mitochondrial DNA control region in the avian genus Alectoris. J. Mol. Evol. 1998, 47, 449–462. [Google Scholar] [CrossRef]
  12. Moore, W.S.; Defilippis, V.R. The Window of Taxonomic Resolution for Phylogenies Based on Mitochondrial Cytochrome b. In Avian Molecular Evolution and Systematics; Elsevier: Amsterdam, The Netherlands, 1997. [Google Scholar]
  13. Baker, A.J.; Marshall, H.D. Systematics. Mitochondrial Control Region Sequences as Tools for Understanding Evolution-ScienceDirect. In Avian Molecular Evolution and Systematics; Elsevier: Amsterdam, The Netherlands, 1997; pp. 51–82. [Google Scholar]
  14. Bao, W.; Chen, G.; Li, B.; Wu, X.; Shu, J.; Wu, S.; Xu, Q.; Weigend, S. Analysis of genetic diversity and phylogenetic relationships among red jungle fowls and Chinese domestic fowls. Sci. China Ser. C Life Sci. 2008, 51, 560–568. [Google Scholar] [CrossRef] [PubMed]
  15. Lyimo, C.M.; Weigend, A.; Msoffe, P.L.; Hocking, P.M.; Simianer, H.; Weigend, S. Maternal genealogical patterns of chicken breeds sampled in Europe. Anim. Genet. 2015, 46, 447–451. [Google Scholar] [CrossRef]
  16. Fu, Y.; Niu, D.; Luo, J.; Ruan, H.; He, G.Q.; Zhang, Y.P. Studies of the origin of Chinese domestic fowls. Acta Genet. Sin. 2001, 28, 411–417. [Google Scholar]
  17. Godinez, C.J.P.; Dadios, P.J.D.; Espina, D.M.; Matsunaga, M.; Nishibori, M. Population Genetic Structure and Contribution of Philippine Chickens to the Pacific Chicken Diversity Inferred From Mitochondrial DNA. Front. Genet. 2021, 12, 698401. [Google Scholar] [CrossRef] [PubMed]
  18. Kanginakudru, S.; Metta, M.; Jakati, R.D.; Nagaraju, J. Genetic evidence from Indian red jungle fowl corroborates multiple domestication of modern day chicken. BMC Evol. Biol. 2008, 8, 174. [Google Scholar] [CrossRef]
  19. Liu, Y.P.; Wu, G.S.; Yao, Y.G.; Miao, Y.W.; Luikart, G.; Baig, M.; Beja-Pereira, A.; Ding, Z.L.; Palanichamy, M.G.; Zhang, Y.P. Multiple maternal origins of chickens: Out of the Asian jungles. Mol. Phylogenetics Evol. 2006, 38, 12–19. [Google Scholar] [CrossRef] [PubMed]
  20. Miao, Y.W.; Peng, M.S.; Wu, G.S.; Ouyang, Y.N.; Yang, Z.Y.; Yu, N.; Liang, J.P.; Pianchou, G.; Beja-Pereira, A.; Mitra, B.; et al. Chicken domestication: An updated perspective based on mitochondrial genomes. Heredity 2013, 110, 277–282. [Google Scholar] [CrossRef]
  21. Oka, T.; Ino, Y.; Nomura, K.; Kawashima, S.; Kuwayama, T.; Hanada, H.; Amano, T.; Takada, M.; Takahata, N.; Hayashi, Y.; et al. Analysis of mtDNA sequences shows Japanese native chickens have multiple origins. Anim. Genet. 2007, 38, 287–293. [Google Scholar] [CrossRef]
  22. Al-Jumaili, A.S.; Boudali, S.F.; Kebede, A.; Al-Bayatti, S.A.; Essa, A.A.; Ahbara, A.; Aljumaah, R.S.; Alatiyat, R.M.; Mwacharo, J.M.; Bjornstad, G.; et al. The maternal origin of indigenous domestic chicken from the Middle East, the north and the horn of Africa. BMC Genet. 2020, 21, 30. [Google Scholar] [CrossRef]
  23. Sahu, P.K.; Das, B.; Sahoo, L.; Senapati, S.; Nayak, G.D. Genetic relationship and population structure of three Indian local chicken populations as revealed by mtDNA D-loop. Mitochondrial DNA Part A DNA Mapp. Seq. Anal. 2016, 27, 2986–2988. [Google Scholar] [CrossRef]
  24. Liu, L.L.; Yang, Y.Q.; Liu, Y.T.; Yang, N.; Xie, D.; Huang, Z. Complete mitochondrial genome sequence of Zhuxiang chicken (Gallus gallus. domesticus) and its phylogenetic analysis from D-loop region. Mitochondrial DNA B Resour. 2018, 3, 872–873. [Google Scholar] [CrossRef]
  25. Jia, X.X.; Tang, X.J.; Lu, J.X.; Fan, Y.F.; Chen, D.W.; Tang, M.J.; Gu, R.; Gao, Y.S. The investigation of genetic diversity and evolution of Daweishan Mini chicken based on the complete mitochondrial (mt)DNA D-loop region sequence. Mitochondrial DNA Part A DNA Mapp. Seq. Anal. 2016, 27, 3001–3004. [Google Scholar] [CrossRef]
  26. Peng, S.M.; Lin, Q.; Jiang, G.T.; Li, Y.H.; Dai, Q.Z.; He, X.; Yan, H.F. The complete mitochondrial genome of the Dongan black chicken and its phylogenetic analyses. Mitochondrial DNA B Resour. 2018, 3, 1127–1128. [Google Scholar] [CrossRef]
  27. Kawabe, K.; Worawut, R.; Taura, S.; Shimogiri, T.; Nishida, T.; Okamoto, S. Genetic Diversity of mtDNA D-loop Polymorphisms in Laotian Native Fowl Populations. Asian-Australas J. Anim. Sci. 2014, 27, 19–23. [Google Scholar] [CrossRef]
  28. Sambrook, J.E.; Maniatis, T.E.; Fritsch, E.F. Molecular Cloning: A Laboratory Manual; Cold Spring Harbor Laboratories: Cold Spring Harbor, NY, USA, 1989. [Google Scholar]
  29. Fumihito, A.; Miyake, T.; Sumi, S.; Takada, M.; Ohno, S.; Kondo, N. One subspecies of the red junglefowl (Gallus gallus gallus) suffices as the matriarchic ancestor of all domestic breeds. Proc. Natl. Acad. Sci. USA 1994, 91, 12505–12509. [Google Scholar] [CrossRef]
  30. Li, Q.; Zhang, P.; Li, M.; Li, L.; Hu, Y.; Mishra, S.K.; Guo, A.; Li, G.; Li, D.; Duan, Y. Genetic diversity and relationship of Dulong chickens using mitochondrial DNA control region. Mitochondrial DNA B Resour. 2019, 5, 275–280. [Google Scholar] [CrossRef] [PubMed]
  31. Excoffier, L.; Lischer, H.E. Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Mol. Ecol. Resour. 2010, 10, 564–567. [Google Scholar] [CrossRef]
  32. Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol. Biol. Evol. 2013, 30, 2725–2729. [Google Scholar] [CrossRef]
  33. Leigh, J.W.; Bryant, D. POPART: Full-Feature Software for Haplotype Network Construction. Methods Ecol. Evol. 2015, 6, 1110–1116. [Google Scholar] [CrossRef]
  34. Tang, X.-J.; Jia, X.-X.; Fan, Y.-F.; Gao, Y.-S.; Tang, M.-J.; Chen, D.-W.; Zhang, X.-Y.; Lu, J.-X.; Han, W. DNA Barcoding and Variety Identification of Four Chicken (Gallus gallus domesticus) Breeds Based on Mitochondrial DNA D-loop Region. J. Agric. Biotechnol. 2017, 25, 1653–1660. [Google Scholar]
  35. Gao, Y.S.; Jia, X.X.; Tang, X.J.; Fan, Y.F.; Lu, J.X.; Huang, S.H.; Tang, M.J. The genetic diversity of chicken breeds from Jiangxi, assessed with BCDO2 and the complete mitochondrial DNA D-loop region. PLoS ONE 2017, 12, e0173192. [Google Scholar] [CrossRef] [PubMed]
  36. Yin, Y.-J.; Du, W.-L.; Wang, L.; Zhang, T. Genetic Diversity of Lueyang Black-bone Chicken Based on Mitochondrial DNA D-loop Region Sequence. China Anim. Husb. Vet. Med. 2017, 44, 173–179. [Google Scholar]
  37. Gao, Y.-S.; Jia, X.-X.; Tang, X.-J.; Tang, M.-J.; Fan, Y.-F.; Lu, J.-X.; Gu, R.; Ge, Q.-L.; Su, Y.-J. The Genetic Diversity and Origin Analysis of Anyi Tile- like Chickens (Gallus gallus domestiaus) Based onMitochondrial DNAD-loop Sequence. J. Agric. Biotechnol. 2015, 25, 940–944. [Google Scholar]
  38. Mtileni, B.J.; Muchadeyi, F.C.; Maiwashe, A.; Groeneveld, E.; Groeneveld, L.F.; Dzama, K.; Weigend, S. Genetic diversity and conservation of South African indigenous chicken populations. J. Anim. Breed Genet. 2011, 128, 209–218. [Google Scholar] [CrossRef]
  39. Tadano, R.; Sekino, M.; Nishibori, M.; Tsudzuki, M. Microsatellite marker analysis for the genetic relationships among Japanese long-tailed chicken breeds. Poult. Sci. 2007, 86, 460–469. [Google Scholar] [CrossRef]
  40. Cabarles, J.C., Jr. Production potentials of native chickens (Gallus gallus domesticus L.) of Western Visayas, Philippines. Trop. Anim. Health Prod. 2013, 45, 405–410. [Google Scholar] [CrossRef]
  41. Chen, S.-Y.; Zhang, Y.-P. Genetic methods and applications in the study of the origin of domestic animals. Chin. Sci. Bull. 2006, 51, 2469–2475. [Google Scholar]
  42. Sha, Y.-Z.; Gao, C.-X.; Liu, M.-M.; Zhao, S.-G. Evaluation of the genetic diversity of six Chinese indigenous chickens. Asian-Australas. J. Anim. Sci. 2019, 33, 1566–1572. [Google Scholar] [CrossRef]
  43. Bao, W.-B.; Shu, J.-T.; Wang, C.-B.; Zhang, H.-X.; Chen, G.-H. Investigation on Genetic Diversity and Systematic Evolution in Chinese Domestic Fowls and Red Jungle Fowls by Analyzing the mtDNA Control Region. Acta V Eterinaria Zootech. Sin. 2008, 39, 1449–1459. [Google Scholar]
  44. Teinlek, P.; Siripattarapravat, K.; Tirawattanawanich, C. Genetic diversity analysis of Thai indigenous chickens based on complete sequences of mitochondrial DNA D-loop region. Asian-Australas J. Anim. Sci. 2018, 31, 804–811. [Google Scholar] [CrossRef] [PubMed]
  45. Liu, J.-B.; Ding, X.-Z.; Zeng, Y.-F.; Yue, Y.-J.; Guo, X.; Guo, T.-T.; Chu, M.; Wang, F.; Han, J.; Feng, R.; et al. Genetic Diversity and Phylogenetic Evolution of Tibetan Sheep Based on mtDNA D-Loop Sequences. PLoS ONE 2016, 11, e0159308. [Google Scholar] [CrossRef]
  46. Toalombo Vargas, P.A.; Leon, J.M.; Fiallos Ortega, L.R.; Martinez, A.; Villafuerte Gavilanes, A.A.; Delgado, J.V.; Landi, V. Deciphering the Patterns of Genetic Admixture and Diversity in the Ecuadorian Creole Chicken. Animals 2019, 9, 670. [Google Scholar] [CrossRef] [PubMed]
  47. Giontella, A.; Cardinali, I.; Pieramati, C.; Cherchi, R.; Biggio, G.P.; Achilli, A.; Silvestrelli, M.; Lancioni, H. A Genetic Window on Sardinian Native Horse Breeds through Uniparental Molecular Systems. Animals 2020, 10, 1544. [Google Scholar] [CrossRef] [PubMed]
  48. Xia, X.; Huang, G.; Wang, Z.; Sun, J.; Wu, Z.; Chen, N.; Lei, C.; Hanif, Q. Mitogenome Diversity and Maternal Origins of Guangxi Cattle Breeds. Animals 2019, 10, 19. [Google Scholar] [CrossRef]
  49. Lancioni, H.; Di Lorenzo, P.; Cardinali, I.; Ceccobelli, S.; Capodiferro, M.R.; Fichera, A.; Grugni, V.; Semino, O.; Ferretti, L.; Gruppetta, A.; et al. Survey of uniparental genetic markers in the Maltese cattle breed reveals a significant founder effect but does not indicate local domestication. Anim. Genet. 2016, 47, 267–269. [Google Scholar] [CrossRef]
  50. Kusliy, M.A.; Vorobieva, N.V.; Tishkin, A.A.; Makunin, A.I.; Druzhkova, A.S.; Trifonov, V.A.; Iderkhangai, T.O.; Graphodatsky, A.S. Traces of Late Bronze and Early Iron Age Mongolian Horse Mitochondrial Lineages in Modern Populations. Genes 2021, 12, 412. [Google Scholar] [CrossRef]
  51. Ulfah, M.; Perwitasari, D.; Jakaria, J.; Muladno, M.; Farajallah, A. Multiple maternal origins of Indonesian crowing chickens revealed by mitochondrial DNA analysis. Mitochondrial DNA Part A DNA Mapp. Seq. Anal. 2017, 28, 254–262. [Google Scholar] [CrossRef]
  52. Zang, H.; Wang, Y.; Yang, H.H.; He, P.L.; Pan, S.Q.; Geng, Z.Y.; Jin, S.H. Characterisation of the complete mitochondrial genome, genetic diversity and maternal origin of Huainan Partridge chicken. Br. Poult. Sci. 2021, 62, 320–327. [Google Scholar] [CrossRef]
  53. Hoque, R.; Seo, D.W.; Jo, C.; Gotoh, T.; Lee, J.H. Reconstruction of Phylogenetic Relationships of Korean Chickens Using Mitochondrial Dloop Sequences. J. Fac. Agric. Kyushu Univ. 2013, 58, 287–293. [Google Scholar]
  54. Liao, Y.; Mo, G.; Sun, J.; Wei, F.; Liao, D.J. Genetic diversity of Guangxi chicken breeds assessed with microsatellites and the mitochondrial DNA D-loop region. Mol. Biol. Rep. 2016, 43, 415–425. [Google Scholar] [CrossRef]
Genes 14 00998 g001 550
Figure 1. Distribution map of eight local chicken breeds in Sichuan. Note: MY: Miyi chickens; SM: Shimian Caoke chicken; EM: Emei Black chicken; MC: Muchuan silky chicken; PX: Pengxian yellow chicken; JS: Jinyang silky chicken; MH: Mountainous Meihua chicken; JY: Jiuyuan black chicken. The red dot Chengdu is capital of Sichuan province.
Figure 1. Distribution map of eight local chicken breeds in Sichuan. Note: MY: Miyi chickens; SM: Shimian Caoke chicken; EM: Emei Black chicken; MC: Muchuan silky chicken; PX: Pengxian yellow chicken; JS: Jinyang silky chicken; MH: Mountainous Meihua chicken; JY: Jiuyuan black chicken. The red dot Chengdu is capital of Sichuan province.
Genes 14 00998 g001
Genes 14 00998 g002 550
Figure 2. Variable information of the 8 Sichuan chicken breeds. Note: These numbers indicated mutation position. “n” represents the number of individuals per haplotype. Dot (.) represents the identical nucleotide with the reference sequence (NC_007235).
Figure 2. Variable information of the 8 Sichuan chicken breeds. Note: These numbers indicated mutation position. “n” represents the number of individuals per haplotype. Dot (.) represents the identical nucleotide with the reference sequence (NC_007235).
Genes 14 00998 g002
Genes 14 00998 g003 550
Figure 3. Phylogenetic tree of mitochondrial D-loop nucleotide sequences based on 74 haplotypes (30 from reference DNA Genbank representing 13 clades [20] and 44 from indigenous chicken breeds). The evolutionary history was inferred using the neighbor-joining method. The numeral at each branch indicates the bootstrap value of replications (1000 replicates). Bootstrap values lower than 50% are not shown.
Figure 3. Phylogenetic tree of mitochondrial D-loop nucleotide sequences based on 74 haplotypes (30 from reference DNA Genbank representing 13 clades [20] and 44 from indigenous chicken breeds). The evolutionary history was inferred using the neighbor-joining method. The numeral at each branch indicates the bootstrap value of replications (1000 replicates). Bootstrap values lower than 50% are not shown.
Genes 14 00998 g003
Genes 14 00998 g004 550
Figure 4. Maximum parsimony median-joining network of MH with other Sichuan local chicken breeds based on 482 bp control-region sequence. Node sizes are proportional to haplotype frequencies. The lines linking the nodes are proportional to the mutation steps. Different colors within each haplotype ring correspond to different chicken breeds. Five clades (A, B, C, E and G) were formed.
Figure 4. Maximum parsimony median-joining network of MH with other Sichuan local chicken breeds based on 482 bp control-region sequence. Node sizes are proportional to haplotype frequencies. The lines linking the nodes are proportional to the mutation steps. Different colors within each haplotype ring correspond to different chicken breeds. Five clades (A, B, C, E and G) were formed.
Genes 14 00998 g004
Table
Table 1. The base composition of mtDNA D-Loop sequence.
Table 1. The base composition of mtDNA D-Loop sequence.
BreedNumber of
Individuals		Total Number of BasesContent %
T(U)CAGA + TG + C
MH19948230.5128.0928.3813.0258.8941.11
JY4048230.4328.1028.3213.1558.7541.25
MY3948230.5927.8828.4113.1159.0041.00
SM4648230.4728.0728.3513.1158.8241.18
MC2948230.4828.0628.2413.2258.7241.28
PX2748230.3828.1528.2913.1858.6741.33
JS2948230.3728.0928.4113.1358.7941.21
EM3048230.4628.0828.3113.1558.7741.23
Table
Table 2. Distribution of 44 haplotypes in 8 Sichuan chicken breeds.
Table 2. Distribution of 44 haplotypes in 8 Sichuan chicken breeds.
MH(199)JY(40)MY(39)SM(46)MC(29)PX(27)JS(29)EM(30)
H1(2)H5(20)H5(4)H3(10)H5(4)H5(15)H3(1)H5(10)
H2(34)H6(5)H7(5)H7(22)H7(11)H7(1)H5(15)H7(10)
H3(45)H7(7)H12(2)H11(3)H12(1)H14(1)H20(6)H12(1)
H4(1)H24(3)H16(1)H12(6)H16(10)H20(7)H35(2)H16(1)
H5(17)H25(3)H20(2)H20(1)H20(3)H30(1)H37(1)H20(1)
H6(4)H26(1)H23(10)H30(2) H35(1)H38(1)H26(1)
H7(22)H27(1)H24(2)H33(2) H36(1)H39(1)H41(3)
H8(1) H26(3) H40(2)H42(1)
H9(1) H28(1) H43(1)
H10(28) H29(2) H44(1)
H11(2) H30(2)
H12(13) H31(2)
H13(1) H32(1)
H14(3) H33(1)
H15(1) H34(1)
H16(6)
H17(2)
H18(1)
H19(1)
H20(10)
H21(2)
H22(1)
H23(1)
Note: The parenthesized numbers are sample quantities, and the underlined numbers indicate shared haplotypes between populations.
Table
Table 3. Genetic diversity indices of mtDNA control region within 8 Sichuan chicken breeds.
Table 3. Genetic diversity indices of mtDNA control region within 8 Sichuan chicken breeds.
BreednVSsPsHHdπANDTajima’s D *
MH19928325230.8760.0125.9360.567
JY401871170.7090.0073.492−0.571
MY3926224150.9060.0157.3870.687
SM461901970.7410.0104.8170.367
MC291211150.7320.0073.1670.120
PX272081270.6410.0115.0540.289
JS29135880.7000.0136.0590.123
EM301789100.7860.0104.7840.393
Total43946739440.8950.0136.349−0.331
Note: * All of them were not statistically significant in this study (p-values > 0.1), and some of the data in this table came from our previous study [30]. Variable sites: V; singleton sites: Ss, parsimony informative sites: Ps; haplotypes: H; haplotype diversity: Hd; nucleotide diversity: π; average no. of nucleotide differences: AND.
Table
Table 4. Analysis of molecular variance (AMOVA) of 439 mtDNA D-loop sequences from 8 local chicken breeds in Sichuan province.
Table 4. Analysis of molecular variance (AMOVA) of 439 mtDNA D-loop sequences from 8 local chicken breeds in Sichuan province.
Sources of Variationd.f.Sum of SquaresVariance ComponentsPercentage of VariationFstp-Value
Among populations7221.6760.61549 Va18.50--
Within populations4311168.7872.71180 Vb81.50--
Total4381390.4623.32729-0.184980.00
Note: d.f., degree of freedom; Fst, fixation index; Va, among-population variance; Vb, between-individual within-population variance.
Table
Table 5. Genetic differentiation index (Gst), gene flow (Nst), nucleotide divergence (Dxy), and net genetic distance (Da) between each pair of chicken breeds.
Table 5. Genetic differentiation index (Gst), gene flow (Nst), nucleotide divergence (Dxy), and net genetic distance (Da) between each pair of chicken breeds.
MHJYMYSMMCPXJSEM
Coefficient of differentiation and gene flow
MH 18.53715.62618.34537.88813.41511.72812.640
JY4.793 27.88037.46451.4263.8604.0107.071
MY2.7116.564 8.39918.81118.03317.82214.482
SM3.48312.5306.804 10.45123.18522.53016.239
MC3.6428.9525.7325.601 36.30435.44028.150
PX4.7512.8518.37016.86512.470 −2.3960.693
JS4.2612.5397.44015.89811.815−1.036 2.065
EM3.0311.7533.9175.2643.3495.2854.972
Nucleotide divergence and net genetic distance
MH 0.2250.2580.2520.5780.1790.1670.163
JY1.212 0.4380.5210.7350.0360.0420.066
MY1.6521.576 0.1190.2570.2850.3030.215
SM1.3771.3911.396 0.0970.3130.3320.195
MC1.5301.4311.3610.932 0.4890.5280.324
PX1.3280.9301.5871.3461.349 −0.0270.007
JS1.4221.0411.7101.4691.4941.135 0.024
EM1.2840.9321.4881.1991.1561.0361.158
Note: The upper half is made up of the Gst (lower diagonal) and the Nst (upper diagonal). The lower half is made up of the Dxy (below diagonal) and the Da (above diagonal). All the values were enlarged 100 times, with the maximum and minimum values underlined. Relevant data of JY, MY, SM, MC, PX, JS, and EM came from our previous study [30].
Table
Table 6. Composition of haplogroups in MH chickens.
Table 6. Composition of haplogroups in MH chickens.
CladeABCEG
H1(2)H2(34)-H7(22)H21(2)
H5(17)H3(45)-H12(13)H23(1)
H6(4)H4(1)-H13(1)-
H10(28)H8(1)-H16(6)-
H11(2)H9(1)-H17(2)-
H14(3)H15(1)-H20(10)-
H18(1)----
H19(1)----
H22(1)----
Total59830543
Ratio(%)29.65%41.71%027.14%1.50%
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

Ran, B.; Zhu, W.; Zhao, X.; Li, L.; Yi, Z.; Li, M.; Wang, T.; Li, D. Studying Genetic Diversity and Relationships between Mountainous Meihua Chickens Using Mitochondrial DNA Control Region. Genes 2023, 14, 998. https://0-doi-org.brum.beds.ac.uk/10.3390/genes14050998

AMA Style

Ran B, Zhu W, Zhao X, Li L, Yi Z, Li M, Wang T, Li D. Studying Genetic Diversity and Relationships between Mountainous Meihua Chickens Using Mitochondrial DNA Control Region. Genes. 2023; 14(5):998. https://0-doi-org.brum.beds.ac.uk/10.3390/genes14050998

Chicago/Turabian Style

Ran, Bo, Wei Zhu, Xiaoling Zhao, Linxiang Li, Zhixin Yi, Miao Li, Tao Wang, and Diyan Li. 2023. "Studying Genetic Diversity and Relationships between Mountainous Meihua Chickens Using Mitochondrial DNA Control Region" Genes 14, no. 5: 998. https://0-doi-org.brum.beds.ac.uk/10.3390/genes14050998

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