RESEARCH ARTICLE

# Population structure analysis of Yeonsan Ogye using microsatellite markers

Sung Hyun Cho1,#, Seung-Sook Lee2,#, Prabuddha Manjula1, Minjun Kim1, Seung Hwan Lee1, Jun Heon Lee1,*, Dongwon Seo1,*
1Division of Animal and Dairy Science, Chungnam National University, Daejeon 34134, Korea
2Jisan Farm, Nonsan 29009, Korea
*Corresponding author: Jun Heon Lee, Division of Animal and Dairy Science, Chungnam National University, Daejeon 34134, Korea. Tel: +82-42-821-5779, E-mail: junheon@cnu.ac.kr
*Corresponding author: Dongwon Seo, Division of Animal and Dairy Science, Chungnam National University, Daejeon 34134, Korea. Tel: +82-42-821-5878, E-mail: seotuna@gmail.com

#These authors contributed equally to this work.

© Copyright 2020 Korean Society of Animal Science and Technology. This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received: Sep 10, 2020; Revised: Sep 25, 2020; Accepted: Sep 27, 2020

Published Online: Nov 30, 2020

## Abstract

The Yeonsan Ogye (YO) chicken is a natural heritage of Korea, characterized by black feathers, skin, bones, eyes, and comb. The purebred of YO population has been reared under the natural mating system with no systematic selection and breeding plan. The purpose of this study was to identify the genetic diversity and find the optimal number of population sub-division using 12 polymorphic microsatellite (MS) markers to construct a pedigree-based breeding plan for the YO population. A total of 509 YO birds were used for this study. Genetic diversity and population structure analysis were conducted based on the MS marker genotype information. The overall average polymorphic information content value and expected heterozygosity of the population were 0.586, and 0.642, respectively. The K-mean cluster analysis based on the genetic distance result confirmed that the current YO population can be divided into three ancestry groups. Individuals in each group were evaluated based on their genetic distance to identify the potential candidates for a future breeding plan. This study concludes that a future breeding plan with known pedigree information of selected founder animals, which holds high genetic diversity, could be the best strategy to ensure the conservation of the Korean YO chicken population.

Keywords: Genetic diversity study; K-mean clustering; Microsatellite marker; Yeonsan Ogye chicken

## INTRODUCTION

Knowledge of population genetic diversity is important for effective management and the development of sustainable conservation and improvement strategies. DNA polymorphism analysis is a powerful tool for analyzing genetic variation within and among individuals, families, populations, and other taxa. Various types of molecular genetic markers can find in the eukaryotic genome, such as microsatellite (MS), single-nucleotide polymorphism (SNP), restriction fragment length polymorphism (RFLP), and amplified fragment length polymorphism (AFLP) markers have been developed and utilized in population genetic studies. In addition, DNA profiling methods, which are based on the polymerase chain reaction (PCR) [1] are independent with environmental factors, can provide more accurate and useful information for the genomic evidence.

MS markers are typically genotyped by PCR amplification, followed by the determination of allele size by either gel electrophoresis or capillary electrophoresis, which involves laser-based detection of fluorescence dye modified oligonucleotides [2]. MS markers have also been demonstrated to be useful for estimating genetic relatedness and diversity in several livestock species including, cattle, sheep, and chickens, particularly in indigenous breeds, inbred strains, and commercial lines [36]. Based on their degree of polymorphism and genome coverage, several MS markers have been recommended by the Measurement of Domestic Animal Diversity (MoDAD) for application in diversity studies [7]. Various population genetic parameters have been described to measure genetic variation within and among populations using MS markers, including observed (Ho) and expected heterozygosity (He), polymorphic information content (PIC), number of alleles (Na) and allele richness per locus, genetic distance, and phylogeny. Simple and affordable genotyping methods and analytical tools offering easy interpretation have led to the extensive application of MS markers in genetic diversity studies to implements the conservation program for various native chicken breeds.

In 1980, the Yeonsan Ogye (YO) chicken was designated a natural monument of Korea (no. 265). YO is characterized by black feathers, skin, bones, eyes, and crest. Among native chicken breeds in South Korea, YO has very low economic value, owing to its slow breeding potential and hatching rate. However, the importance of YO as a natural cultural asset must be considered. Korea’s Cultural Heritage Administration has established guidelines to maintain a population of at least 1,000 YO chickens to prevent its extinction. To maintain genetic heterogeneity, it is essential to prevent inbreeding depression. Previous genomic analyses have shown that the YO population has maintained its genetic diversity and a suitable inbreeding coefficient despite the absence of a systematic management system. However, the current YO breeding group possesses no pedigree information to maintain a conservation population of over 1,000 animals via natural mating and artificial hatching. In each generation, male and female birds were selected based on the phenotypic traits (qualitative measures) such as completely black color for feather, comb, and shank, without considering their genetic diversity or the inbreeding ratio. Therefore, it is necessary to establish a breeding system for YO population that prevents inbreeding depression likely happen in the future and to maintain genetic diversity. Such efforts may prevent the loss of this nationally important YO genetic resource. This study aimed to assess the genetic diversity of current YO population and find the optimal number of clusters (the groups of birds that differ by genetic distance, and have low relatedness) using 19 polymorphic microsatellite (MS) markers and to construct a pedigree-based breeding plan for the YO population.

## MATERIALS AND METHODS

Sample collection and DNA extraction

In the present study, a complete sample strategy was implemented to collect the sample of YO population maintained at a private farm in Nonsan city, Korea. Approximately, 2 mL of whole blood samples were collected from 509 birds (> 5 months old males and females) and were stored in tubes containing ethylenediaminetetraacetic acid (EDTA) to prevent coagulation before DNA extraction step. At the same time, all the birds were tagged with unique electronic radio-frequency identification (RFID) tags and recorded to facilitate the identification of each individual (Fig. 1). Blood samples were subjected to genomic DNA (gDNA) extraction using the PrimePrep Genomic DNA Extraction Kit from Blood (Genet Bio, Daejeon, Korea). To assess the quality of the extracted gDNA, its concentration and purity were measured using a NanoDrop 2000c spectrophotometer. The extracted gDNA samples were stored at –20°C until the experiment began.

Fig. 1. Yeonsan Ogye (YO) chicken with radio-frequency identification (RFID) tags on its legs.
Microsatellite marker amplification test and marker selection

Initially, a total of 19 highly polymorphic MS markers were tested to investigate their polymorphic status in the YO breed. This set of markers were selected from the list of MS markers recommended by Food and Agriculture Organization of the United Nations (FAO) and a previous study of Seo et al. [8]. Polymorphism level of 19 markers was initially evaluated using single-plex PCR and agarose gel electrophoresis. Of these, 12 highly polymorphic markers were selected for multiplex PCR and genotyping. The 3’ end of each marker was modified to add four types of fluorescent substances (FAM, NED, VIC, and PET) for multiplex PCR and fragment analysis. To confirm the genotype (Table 1), a single PCR was also performed using HS Prime Taq Premix (Genet Bio) on YO test samples. The following amplification conditions were used in each PCR amplification: pre-denaturation at 95°C for 10 min, followed by denaturation, annealing, and extension processes repeated 30 times at 95°C for 30 s, 60°C for 30 s, and 72°C for 1 min, respectively, followed by a final extension at 72°C for 60 min. Fragment sizes of the PCR amplification products was then obtained by fragment analysis using the GeneScan-500 LIZ (Applied Biosystems, Foster City, CA, USA) size standard marker and Genetic Analyzer (GA)3730xl (Applied Biosystems). Allele size was assigned using the GeneMapper ver. 3.7 software (Applied Biosystems).

Table 1. Primer information for the 19 candidate microsatellite (MS) markers
No Marker Chr Dye Primer (5’-3’) Allele size range
1 ADL0259 9 VIC F CTCATTGCAGAGGAAGTTCT 107–129
R GTAATGGAGGATGCTCAGGT
2 ADL0268 1 FAM F CTCCACCCCTCTCAGAACTA 105–117
R CAACTTCCCATCTACCTACT
3 ADL0293 17 NED F GTAATCTAGAAACCCCATCT 105–119
R ACATACCGCAGTCTTTGTTC
4 ADL0304 18 NED F GGGGAGGAACTCTGGAAATG 137–159
R CCTCATGCTTCGTGCTTTTT
5 ADL0317 4 FAM F AGTTGGTTTCAGCCATCCAT 178–204
R CCCAGAGCACACTGTCACTG
6 GCT0016 9 PET F TCCAAGGTTCTCCAGTTC 109–125
R GGCATAAGGATAGCAACAG
7 LEI0074 26 VIC F GACCTGGTCCTGACATGGGTG 224–240
R GTTTGCTGATTAGCCATCGCG
8 LEI0094 4 VIC F GATCTCACCAGTATGAGCTGC 254–280
R TCTCACACTGTAACACAGTGC
9 LEI0141 2 FAM F CGCATTTGATGCATAACACATG 220–242
R AAGGCAAACTCAGCTGGAACG
10 MCW0029 5 PET F GTGGACACCCATTTGTACCCTATG 139–189
R CATGCAATTCAGGACCGTGCA
11 MCW0087 2 FAM F ATTTCTGCAGCCAACTTGGAG 267–283
R CTCAGGCAGTTCTCAAGAACA
12 MCW0104 13 NED F TAGCACAACTCAAGCTGTGAG 189–225
R AGACTTGCACAGCTGTGTACC
13 MCW0123 14 FAM F CCACTAGAAAAGAACATCCTC 79–89
R GGCTGATGTAAGAAGGGATGA
14 MCW0145 1 VIC F ACTTTATTCTCCAAATTTGGCT 181–211
R AAACACAATGGCAACGGAAAC
15 MCW0228 10 PET F GATCTCTGCATTACAAGCATG 221–239
R TTGCTGACCTGCTCATGCAAG
16 MCW0264 2 PET F CTTACTTTTCACGACAGAAGC 224–240
R AGACTGAGTCACACTCGTAAG
17 MCW0330 17 NED F TGGACCTCATCAGTCTGACAG 254–286
R AATGTTCTCATAGAGTTCCTGC
18 ROS0013 5 NED F TGCTGCTCCTGGRAAATTG 220–242
R GAAAAGCCATGGAGGAATCA
19 ROS0083 13 VIC F CATTACAGCTCAGTGTTGGCA 109–129
R TTGCAAGTGCTCTCCCATC

The bold markers are the selected 12 microsatellite markers for the population study.

Chr, chromosome number; Dye, fluorescent dye.

Population genetic analyses and breeding plans

To obtain the genetic parameter for 12 MS markers, we performed basic population-level statistical and genetic diversity analyses using the ‘adegenet’ package in the R software [9]. The total number of alleles, polymorphic information content (PIC), Ho and He were calculated.

Multi-dimensional scaling (MDS) plot and phylogenetic tree analyses were performed using individual genetic distance information derived using the Euclidean distance method and genetic distances between clustered groups were identified based on Nei’s standard genetic distance values obtained using the ‘poppr’ R package [10], which can be represented as follows:

${D}_{Ne{i}^{\prime }s}=-\text{ln}\left(\frac{\Sigma \text{\hspace{0.17em}}pop{1}_{i}\cdot pop{2}_{i}}{\sqrt{\left(\Sigma \text{\hspace{0.17em}}pop{1}_{i}^{2}\right)\left(\Sigma \text{\hspace{0.17em}}pop{2}_{i}^{2}\right)}}\right)$

Where pop1i and pop2i are the allele frequencies at the ith locus in populations 1 and 2, respectively.

To obtain the genetic structure in current YO population (i.e., for optimal YO group division), k-means clustering was performed using the ‘NbClust’ R software package, and the obtained grouping information, and was heatmap analyzed using the ‘pheatmap’ R software package to visualize the genetic distance between identified groups.

To prepare a mating plan that will minimizes the inbreeding within the YO group, the possible breeding groups were identify by cluster analysis and genetic distance data.

## RESULTS AND DISCUSSION

For effective management and conservation of YO chicken as a natural heritage, it is necessary to separate breeding groups and maintain the pedigree by establishing an all-in-all-out system for generation replacement. However, due to the lack of pedigree and breeding record information of current YO population, we have distinguished several groups of YO based on their genetic distance using genotype information from 12 MS markers to establish an initial breeding plan.

Selection of highly polymorphic microsatellite markers for Yeonsan Ogye population

A total of 19 MS markers were tested to identify highly polymorphic marker combinations in YO chicken population. The 12 selected MS markers (ADL0259, ADL0268, ADL017, LEI0074, LEI0094, LEI0141, MCW0029, MCW0087, MCW0127, MCW0264, MCW0330, and ROS0083) contained an average of five alleles per locus. The Ho, He, and PIC of each marker in the YO population are listed in Table 2.

Table 2. Allele distribution, polymorphism information content (PIC), and observed (Ho) and expected heterozygosity (He) of the 12 microsatellite markers
Marker Set Identified allele size PIC Ho He Allele number
ADL0259 1 107 109 111 113 127 0.549 0.595 0.617 5
ADL0268 2 102 108 110 112 114 0.736 0.727 0.773 5
ADL0317 1 182 186 188 196 198 0.308 0.389 0.356 5
LEI0074 2 227 229 235 237 239 243 0.684 0.715 0.727 6
LEI0094 1 252 262 264 266 278 282 0.686 0.766 0.734 6
LEI0141 2 222 230 234 244 0.623 0.251 0.686 4
MCW0029 2 149 151 161 163 181 185 187 0.650 0.697 0.693 7
MCW0087 1 267 271 273 277 281 283 0.726 0.633 0.767 6
MCW0123 1 80 86 90 0.479 0.544 0.541 3
MCW0264 2 230 232 238 240 244 0.403 0.448 0.487 5
MCW0330 1 268 276 286 0.496 0.536 0.587 3
ROS0083 2 109 117 119 123 129 0.691 0.760 0.730 5

The average PIC of the selected markers was 0.586, and the average He was 0.642. Hinderbrand et al. [11] reported that MS markers with PIC > 0.5 and He > 0.6 are suitable for population genetic studies. Thus, all the selected markers except ADL0317 (PIC ≥ 0.4) had sufficient PIC for individual identification and lineage classification in the YO population. Two markers, ADL0268 and MCW0087 had exceptionally high PIC levels, of 7.0 or higher.

The markers selected in this study had higher Ho and He values than those identified in a previous study that have used 25 MS markers [12]. However, these results may have caused by the relatively large number of samples used (509 samples), compared with the previous study (27 samples). Notably, the combination of markers selected in this study provides sufficient polymorphism to distinguish individuals and breeds with a relatively small number of markers. Since all 12 markers were included in previous genetic diversity analysis of eight chicken varieties [13], we can indirectly infer the genetic distance and relationship between Ogye and various other chicken breeds in Korea. Choi et al. [13] analyzed the genotypes of 20–30 chickens in different chicken groups including hybrid population to obtain a combination of 12 markers with PIC > 0.7 and a very low probability of identity (PI). These results are consistent with the findings of the present study, suggesting that our 12 selected markers may be sufficient for explaining genetic diversity and distances within the YO population.

Diversity and clustering analyses using 12 microsatellite markers

The result of MDS analysis confirmed that all birds were evenly distributed, without bias. This finding was expected because the samples used in this study were obtained exclusively from one large population. Therefore, these chickens have shared the common alleles (Fig. 2).

Fig. 2. Multi-dimensional scaling analysis results obtained using 12 selected microsatellite markers, indicating no clear differentiation among birds.

To design a breeding plan that minimizes inbreeding in a population without existing lineage information, it is important to calculate genetic distances using genotype information, and mating birds from different groups that identified to have relatively far genetic distance, while maintaining a new lineage record. Allele frequency for each bird and group is calculated to determine the genetic distance among each bird and group [14]. We performed a k-means clustering analysis to determine the optimal number of YO groups using genotype information obtained from the selected MS markers. By gradually increasing the number of clusters from 2 to 15, we determined the optimal number of groups with the smallest within-group variance and the largest between-group variance. This analysis results showed that the optimal k-value of YO groups was three (Fig. 3). In addition, based on the MDS, and phylogenetic tree analysis results, we further determined that the YO population could be separated into three groups (Figs. 4 and 5). These results imply that the YO chicken population is maintained as three ancestral groups. However, each group maintained a similar and close genetic distance from 0.16 to 0.17, and within the group, the genetic distance was very close with an average of 0.694, confirmed that their genetic similarities were high.

Fig. 3. K-means cluster analysis of the YO population. The best fitted cluster was k = 3 value. YO, Yeonsan Ogye.
Fig. 4. The results for the phylogenetic tree analysis using YO population. This indicates three populations can be identified in this population. pop, populations; YO, Yeonsan Ogye.
Fig. 5. MDS plots for K = 3 (A) and K = 5 (B) groups in the YO population. Each group is identified by the shape and color of the spot. pop, populations; MDS, multi-dimensional scaling; YO, Yeonsan Ogye.

If the YO population is divided based on the three suggested groups (with lowest within-group variance and highest between-group variance), the genetic distance between the suggested groups did not much differ as indicated in the Fig. 6A, suggest that mating among three-ancestry groups would not be effective and sustainable. Furthermore, it is very difficult to maintain a sufficient number of generations for future replacement by keeping only three groups which may also lead to high inbreeding and low genetic diversity in the future.

Fig. 6. Heatmap of genetic distances using 3 and 5 groups in the YO population. Blue and red indicate closer and farther genetic distances, respectively. pop, populations; YO, Yeonsan Ogye.

These limitations can be overcome by either increasing number of ancestry group or immigration. In the context of increasing ancestry group, we try to increase the number of ancestry groups up to 5 (k = 5) and analyses the genetic distance between groups. As summarized in heatmap in Figs. 5B, 6B, and 7 the genetic distance between groups was varied and range between 0.119 (pop1 vs pop2) and 0.261 (pop1 vs pop3). It was clear that the increasing number of groups possibly result in more genetically distant ancestry populations. Maintaining a sufficient number of the breeding groups at the initial stage is more important as it provides several options to decide the breeding groups to maintain the effective breeding plan. Having five base populations would be sufficient to increase the number of breeding groups that possess genetically distal birds. Rotational breeding among these five groups is more effective instead of keeping only three groups with low genetic distance.

Fig. 7. Represent MDS plot of K = 5 groups with additional dimensions. (A) Clustering of 5 groups at dimension 1 and 2, (B) clustering of 5 groups at dimension 1 and 3, and (C) clustering of 5 groups at dimension 1 and 4. pop, populations; MDS, multi-dimensional scaling; YO, Yeonsan Ogye.

In general, inbreeding has been reported to have a greater impact on population extinct than natural selection. In addition, inbreeding reduces the heterogeneity of the population and increases homozygosity, making it susceptible to recessive genetic diseases. The heterogeneity of a specific population becomes a measure of genetic diversity, and it is essential to maintain the population by adapting to various environmental changes and hereditary diseases. However, YO has been continued strong selection to maintain its unique appearance as a natural monument while maintaining a closed chicken population. In an environment where pedigrees are not recorded, the influence of phenotypic selections can have a great influence on the gene pool, and there is a possibility that the appearance selection of a certain standard may cause inbreeding by selecting only a single ancestry group. Our recent study of Major Histocompatibility Complex B region variation showed that YO population maintain in National Institute of Animal Science (NIAS) has unique genetic variation compare to other chicken populations [15]. MHC-B linked microsatellite marker LEI0258 diversity in the present study population showed that they have 9 alleles with varied frequencies (unpublished data). Keeping a diverse population will be improved the adaptive genetic variation and benefits their immunogenetic capacity to respond to the diseases in a free-range production system.

Therefore, this study provided sufficient information to understand the current genetic diversity of the YO population and the number of ancestry groups that can be used to increase or maintain the diverse population for future improvement.

Fortunately, 20–30 years ago, when the restoration project of Korean native chicken was being carried out, some groups of YO were delivered to Seoul National University and the National Institute of Animal Science in Chungbuk, Chungnam, and Gyeongbuk to maintain the purebred population. According to Roh et al. [12], each of these Ogye populations has the same phenotypic characteristics but different genetic diversity due to different population size and breeding strategies. Whereas the Gyeongbuk population is different from the other three populations because this population might be maintained as a crossbred despite their similar phenotypic characteristics. In particular, there is a difference between the four groups of YO that was not strongly selected by their appearance. Investigating the genetic diversity of these populations may find new ancestry groups to find clues that could increase genetic diversity. In addition, It has been reported that crossbreeding between genetically close lineages can alleviate the increase in effective population size and inbreeding [16]. Thus, as a second strategy, mutual exchange of population for breeding and effectively population management would be a successful conservation strategy for the YO chicken population.

## CONCLUSION

The results of this study indicate that 12 MS markers represent a viable approach for determining genetic diversity and provide basic information for preparing a pedigree to establish an appropriate YO breeding plan. Based on genetic distances, we have identified three clusters of YO chickens and further divided into five clusters to implement the systematic breeding plan, in which rotational mating between these clusters is considered more favorable due to their different genetic distance. To sustain the unique phenotypic characteristics of YO chickens as a national monument of Korea, it also is necessary to consider the genetic information of individual birds in each breeding group and possible exchange of breeding materials among the available resource populations.

## Funding sources

This work was supported by research fund of Chungnam National University (2020-0537-01).

Not applicable.

## Availability of data and material

Upon reasonable request, the datasets of this study can be available from the corresponding author.

## Authors’ contributions

Conceptualization: Seo D, Lee JH, Lee SH.

Data curation: Cho SH, Lee SS, Kim M.

Formal analysis: Cho SH, Lee SS, Seo D.

Methodology: Cho SH, Lee SS, Kim M.

Software: Cho SH, Manjula P, Kim M.

Validation: Lee JH, Seo D.

Investigation: Cho SH, Lee SS.

Writing - original draft: Cho SH, Lee SS.

Writing - review & editing: Cho SH, Lee SS, Manjula P, Lee SH, Seo D, Lee JH.

## Ethics approval and consent to participate

This article does not require IRB/IACUC approval because there are no human and animal participants.

## REFERENCES

1.

Karp A, Kresovich S, Bhat KV, Ayad WG, Hodgkin T. Molecular tools in plant genetic resources conservation: a guide to the technologies. Rome: International Plant Genetic Resources. 1997IPGRI Technical Bulletin No. 2

2.

Spathis R, Lum JK. An updated validation of Promega’s PowerPlex® 16 system: high throughput databasing under reduced PCR volume conditions on Applied Biosystem’s 96 capillary 3730xl DNA Analyzer. J Forensic Sci. 2008; 53:1353-7

3.

Weigend S, Romanov MN. Current strategies for the assessment and evaluation of genetic diversity in chicken resources. World’s Poult Sci J. 2001; 57:275-88

4.

Zhang X, Leung F, Chan D, Yang G, Wu C. Genetic diversity of Chinese native chicken breeds based on protein polymorphism, randomly amplified polymorphic DNA, and microsatellite polymorphism. Poult Sci. 2002; 81:1463-72

5.

Dessie T, Tadesse M, Yami A, Peters KJ. Village chicken production systems in Ethiopia: 1. flock characteristics and performance. Livest Res Rural Dev. 2003; :15-9

6.

Sheriff O, Alemayehu K. Genetic diversity studies using microsatellite markers and their contribution in supporting sustainable sheep breeding programs: a review. Cogent Food Agric. 2018; :4-1459062

7.

ISAG/FAO Standing Committee. Secondary guidelines for development of national farm animal genetic resources management plans. Measurement of domestic animal diversity (MoDAD): recommended microsatellite markers. Rome: Food and Agriculture Organization of the United Nations. 2004

8.

Seo D, Hoque M, Choi N, Sultana H, Park H, Heo K, et al. Discrimination of Korean native chicken lines using fifteen selected microsatellite markers. Asian-Australas J Anim Sci. 2013; 26:316-22

9.

Jombart T, Ahmed I. adegenet 1.3-1: new tools for the analysis of genome-wide SNP data. Bioinformatics. 2011; 27:3070-1

10.

Kamvar Z, Tabima J, Grünwald N. Poppr: an R package for genetic analysis of populations with clonal, partially clonal, and/or sexual reproduction. PeerJ. 2014; 2:e281

11.

Hildebrand CE, Torney DC, Wagner RP. Informativeness of polymorphic DNA markers.In In: Cooper NG, editor.editor The human genome project: deciphering the blueprint of heredity. Mill Valley, CA: University Science Books. 1994; p p. 100-2

12.

Roh HJ, Kim KW, Lee JW, Jeon DY, Kim SC, Jeon IS, et al. Genetic diversity and relationship of Ogye population in Korea using 25 microsatellite markers. Korean J Poult Sci. 2018; 45:229-36

13.

Choi NR, Seo DW, Jemaa SB, Sultana H, Heo KN, Jo C, et al. Discrimination of the commercial Korean native chicken population using microsatellite markers. J Anim Sci Technol. 2015; 57:1-8

14.

Takezaki N, Nei M. Genetic distances and reconstruction of phylogenetic trees from microsatellite DNA. Genetics. 1996; 144:389-99

15.

Manjula P, Fulton JE, Seo D, Lee JH. Major histocompatibility complex B variability in Korean native chicken breeds. Poult Sci. 2020; 99:4704-13

16.

Piyasatian N, Kinghorn B. Balancing genetic diversity, genetic merit and population viability in conservation programmes. J Anim Breed Genet. 2003; 120:137-49

##### 2019 JCR Impact Factor: 1.685

The 2019 Journal Citation Reports (JCR) was announced, and the impact factor of JAST was determined to be 1.685.

We would like to ask for your continued interest and support in our journal.

Thank you.

JAST Editorial Office

I don't want to open this window for a day.