Identification of new major histocompatibility complex-B Haplotypes in Bangladesh native chickens.

OBJECTIVE: The major histocompatibility complex in chicken demonstrates a great range of variations within varities, breeds, populations and that can eventually influence their immuneresponses. The preset study was conducted to understand the major histocompatibility complex-B (MHC-B) variability in five major populations of Bangladesh native chicken: Aseel, Hilly, Junglefowl, Non-descript Deshi, and Naked Neck. METHODS: These five major populations of Bangladesh native chicken were analyzed with a subset of 89 single nucleotide polymorphisms (SNPs) in the high-density MHC-B SNP panel and Kompetitive Allele-Specific polymerase chain reaction genotyping was applied. To explore haplotype diversity within these populations, the results were analyzed both manually and computationally using PHASE 2.1 program. The phylogenetic investigations were also performed using MrBayes program. RESULTS: A total of 136 unique haplotypes were identified within these five Bangladesh chicken populations, and only one was shared (between Hilly and Naked Neck). Phylogenetic analysis showed no distinct haplotype clustering among the five populations, although they were shared in distinct clades; notably, the first clade lacked Naked Neck haplotypes. CONCLUSION: The present study discovered a set of unique MHC-B haplotypes in Bangladesh chickens that could possibly cause varied immune reponses. However, further investigations are required to evaluate their relationships with global chicken populations.

Major histocompatibility complex genes exhibit a potential immunological role in mixed Eimeria-infected broiler cecum analyzed using RNA sequencing.

OBJECTIVE: This study was conducted to investigate the differential expression of the major histocompatibility complex (MHC) gene region in Eimeria-infected broiler. METHODS: We profiled gene expression of Eimeria-infected and uninfected ceca of broilers sampled at 4, 7, and 21 days post-infection (dpi) using RNA sequencing. Differentially expressed genes (DEGs) between two sample groups were identified at each time point. DEGs located on chicken chromosome 16 were used for further analysis. Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis was conducted for the functional annotation of DEGs. RESULTS: Fourteen significant (false discovery rate <0.1) DEGs were identified at 4 and 7 dpi and categorized into three groups: MHC-Y class I genes, MHC-B region genes, and non-MHC genes. In Eimeria-infected broilers, MHC-Y class I genes were upregulated at 4 dpi but downregulated at 7 dpi. This result implies that MHC-Y class I genes initially activated an immune response, which was then suppressed by Eimeria. Of the MHC-B region genes, the DMB1 gene was upregulated, and TAP-related genes significantly implemented antigen processing for MHC class I at 4 dpi, which was supported by KEGG pathway analysis. CONCLUSION: This study is the first to investigate MHC gene responses to coccidia infection in chickens using RNA sequencing. MHC-B and MHC-Y genes showed their immune responses in reaction to Eimeria infection. These findings are valuable for understanding chicken MHC gene function.

Application of next-generation sequencing for the high-resolution typing of MHC-B in Korean native chicken.

The major histocompatibility complex-B (MHC-B) region of chicken is crucially important in their immunogenesis and highly diverse among different breeds, lines, and even populations. Because it determines the resistance/susceptibility to numerous infectious diseases, it is important to analyze this genomic region, particularly classical class I and II genes, to determine the variation and diversity that ultimately affect antigen presentation. This study investigated five lines of indigenous Korean native chicken (KNC) and the Ogye breed using next-generation sequencing (NGS) data with Geneious Prime-based assembly and variant calling with the Genome Analysis Toolkit (GATK) best practices pipeline. The consensus sequences of MHC-B (BG1-BF2) were obtained for each chicken line/breed and their variants were analyzed. All of the Korean native chicken lines possessed an excessive number of variants, including an ample amount of high-impact variants that provided useful information regarding modified major histocompatibility complex molecules. The study confirmed that next-generation sequencing techniques can effectively be used to detect MHC variabilities and the KNC lines are highly diverse for the MHC-B region, suggesting a substantial divergence from red junglefowl.

Time-series transcriptome analysis identified differentially expressed genes in broiler chicken infected with mixed Eimeria species.

Coccidiosis caused by the Eimeria species is a highly problematic disease in the chicken industry. Here, we used RNA sequencing to observe the time-dependent host responses of Eimeria-infected chickens to examine the genes and biological functions associated with immunity to the parasite. Transcriptome analysis was performed at three time points: 4, 7, and 21 days post-infection (dpi). Based on the changes in gene expression patterns, we defined three groups of genes that showed differential expression. This enabled us to capture evidence of endoplasmic reticulum stress at the initial stage of Eimeria infection. Furthermore, we found that innate immune responses against the parasite were activated at the first exposure; they then showed gradual normalization. Although the cytokine-cytokine receptor interaction pathway was significantly operative at 4 dpi, its downregulation led to an anti-inflammatory effect. Additionally, the construction of gene co-expression networks enabled identification of immunoregulation hub genes and critical pattern recognition receptors after Eimeria infection. Our results provide a detailed understanding of the host-pathogen interaction between chicken and Eimeria. The clusters of genes defined in this study can be utilized to improve chickens for coccidiosis control.

A retroviral insertion in the tyrosinase (TYR) gene is associated with the recessive white plumage color in the Yeonsan Ogye chicken.

The recessive white (locus c) phenotype observed in chickens is associated with three alleles (recessive white c, albino ca , and red-eyed white cre ) and causative mutations in the tyrosinase (TYR) gene. The recessive white mutation (c) inhibits the transcription of TYR exon 5 due to a retroviral sequence insertion in intron 4. In this study, we genotyped and sequenced the insertion in TYR intron 4 to identify the mutation causing the unusual white plumage of Yeonsan Ogye chickens, which normally have black plumage. The white chickens had a homozygous recessive white genotype that matched the sequence of the recessive white type, and the inserted sequence exhibited 98% identity with the avian leukosis virus ev-1 sequence. In comparison, brindle and normal chickens had the homozygous color genotype, and their sequences were the same as the wild-type sequence, indicating that this phenotype is derived from other mutation(s). In conclusion, white chickens have a recessive white mutation allele. Since the size of the sample used in this study was limited, further research through securing additional samples to perform validation studies is necessary. Therefore, after validation studies, a selection system for conserving the phenotypic characteristics and genetic diversity of the population could be established if additional studies to elucidate specific phenotype-related genes in Yeonsan Ogye are performed.

Comparison of Selection Signatures between Korean Native and Commercial Chickens Using 600K SNP Array Data.

Korean native chickens (KNCs) comprise an indigenous chicken breed of South Korea that was restored through a government project in the 1990s. The KNC population has not been developed well and has mostly been used to maintain purebred populations in the government research institution. We investigated the genetic features of the KNC population in a selection signal study for the efficient improvement of this breed. We used 600K single nucleotide polymorphism data sampled from 191 KNCs (NG, 38; NL, 29; NR, 52; NW, 39; and NY, 33) and 54 commercial chickens (Hy-line Brown, 10; Lohmann Brown, 10; Arbor Acres, 10; Cobb, 12; and Ross, 12). Haplotype phasing was performed using EAGLE software as the initial step for the primary data analysis. Pre-processed data were analyzed to detect selection signals using the 'rehh' package in R software. A few common signatures of selection were identified in KNCs. Most quantitative trait locus regions identified as candidate regions were associated with traits related to reproductive organs, eggshell characteristics, immunity, and organ development. Block patterns with high linkage disequilibrium values were observed for LPP, IGF11, LMNB2, ERBB4, GABRB2, NTM, APOO, PLOA1, CNTN1, NTSR1, DEF3, CELF1, and MEF2D genes, among regions with confirmed selection signals. NL and NW lines contained a considerable number of selective sweep regions related to broilers and layers, respectively. We recommend focusing on improving the egg and meat traits of KNC NL and NW lines, respectively, while improving multiple traits for the other lines.

High Levels of Genetic Variation in MHC-Linked Microsatellite Markers from Native Chicken Breeds.

The major histocompatibility complex (MHC) is a highly polymorphic gene region that regulates cellular communication in all specific immune responses. In this study, we investigated 11 microsatellite (MS) markers in the MHC-B region of chicken populations from four countries: Sri Lanka, Bangladesh, South Korea, and Nigeria. The MS markers were divided into two sets. Set 1 included five novel MS markers, which we assessed using 192 samples from 21 populations. Set 2 included six previously reported markers, which we assessed using 881 samples from 29 populations. The Set 1 MS markers had lower polymorphism (polymorphic information content (PIC) < 0.5) than the Set 2 markers (PIC = 0.4-0.9). In all populations, the LEI0258 marker was the most polymorphic, with a total of 38 alleles (PIC = 0.912, expected heterozygosity (He) = 0.918). Local populations from Sri Lanka, Bangladesh, and Nigeria had higher allele diversity and more haplotypes for Set 2 MS markers than Korean and commercial populations. The Sri Lankan Karuwalagaswewa village population had the highest MHC diversity (mean allele number = 8.17, He = 0.657), whereas the white leghorn population had the lowest (mean allele number = 2.33, He = 0.342). A total of 409 haplotypes (89 shared and 320 unique), with a range of 4 (Rhode Island red) to 46 (Karuwalagaswewa village (TA)), were identified. Among the shared haplotypes, the B21-like haplotype was identified in 15 populations. The genetic relationship observed in a neighbour-joining tree based on the DA distance agreed with the breeding histories and geographic separations. The results indicated high MHC diversity in the local chicken populations. The difference in the allelic pattern among populations presumably reflects the effects of different genotypes, environments, geographic variation, and breeding policies in each country. The selection of MHC allele in domestic poultry can vary due to intensification of poultry production. Preserved MHC diversity in local chicken provides a great opportunity for future studies that address the relationships between MHC polymorphisms and differential immune responses.

Identification of Target Chicken Populations by Machine Learning Models Using the Minimum Number of SNPs.

A marker combination capable of classifying a specific chicken population could improve commercial value by increasing consumer confidence with respect to the origin of the population. This would facilitate the protection of native genetic resources in the market of each country. In this study, a total of 283 samples from 20 lines, which consisted of Korean native chickens, commercial native chickens, and commercial broilers with a layer population, were analyzed to determine the optimal marker combination comprising the minimum number of markers, using a 600 k high-density single nucleotide polymorphism (SNP) array. Machine learning algorithms, a genome-wide association study (GWAS), linkage disequilibrium (LD) analysis, and principal component analysis (PCA) were used to distinguish a target (case) group for comparison with control chicken groups. In the processing of marker selection, a total of 47,303 SNPs were used for classifying chicken populations; 96 LD-pruned SNPs (50 SNPs per LD block) served as the best marker combination for target chicken classification. Moreover, 36, 44, and 8 SNPs were selected as the minimum numbers of markers by the AdaBoost (AB), Random Forest (RF), and Decision Tree (DT) machine learning classification models, which had accuracy rates of 99.6%, 98.0%, and 97.9%, respectively. The selected marker combinations increased the genetic distance and fixation index (Fst) values between the case and control groups, and they reduced the number of genetic components required, confirming that efficient classification of the groups was possible by using a small number of marker sets. In a verification study including additional chicken breeds and samples (12 lines and 182 samples), the accuracy did not significantly change, and the target chicken group could be clearly distinguished from the other populations. The GWAS, PCA, and machine learning algorithms used in this study can be applied efficiently, to determine the optimal marker combination with the minimum number of markers that can distinguish the target population among a large number of SNP markers.

Major histocompatibility complex B variability in Korean native chicken breeds.

Adaptive genetic variations have direct influences on the fitness traits of the animal. The major histocompatibility complex B (MHC-B) region is responsible for adaptive and innate immune responses in chickens. In native Korean chicken breeds, no information on serologically defined B haplotypes is available. We investigated the MHC-B diversity in 5 restored lines of Korean native chicken and Ogye chicken breeds using a recently described MHC-B single-nucleotide polymorphism (SNP) panel and the MHC-linked LEI0258 variable number of tandem repeat marker. High SNP haplotype diversity was observed in Korean native chicken breeds with an average of 9.7 MHC-B SNP haplotypes per line. The total number of haplotypes ranged from 6 to 12 per line, and population-specific haplotypes ranged from 3 to 4. A total of 41 BSNP haplotypes, including 26 novel population-specific haplotypes and 15 common haplotypes, were reported over all populations. The 15 common haplotypes included 7 novel and 8 previously reported standard haplotypes. Selection and breeding evidence supports the observation of common haplotypes between the Korean native chicken and exotic breeds. Similarly, the LEI0258 marker showed allele variation, between 193 bp and 474 bp having 5 to 8 alleles per population. Some of these alleles (193, 249, 309, and 443 bp) were shared and more frequently observed. Comparison between SNP haplotypes and LEI0258 allele sizes for the same samples showed that some LEI0258 allele sizes correspond to more than one BSNP haplotype. The use of the MHC-B SNP panel greatly enhances the identification of MHC diversity compared with the sole use of the LEI0258 marker in native chicken populations.

Genetic diversity of MHC-B in 12 chicken populations in Korea revealed by single-nucleotide polymorphisms.

This study used a single-nucleotide polymorphism (SNP) panel to characterise the diversity in the major histocompatibility complex B region (MHC-B) in 12 chicken populations in Korea. Samples were genotyped for 96 MHC-B SNPs using an Illumina GoldenGate genotyping assay. The MHC-B SNP haplotypes were predicted using 58 informative SNPs and a coalescence-based Bayesian algorithm implemented by the PHASE program and a manual curation process. In total, 117 haplotypes, including 24 shared and 93 unique haplotypes, were identified. The unique haplotype numbers ranged from 0 in Rhode Island Red to 32 in the Korean native commercial chicken population 2 ("Hanhyup-3ho"). Population and haplotype principal component analysis (PCA) indicated no clear population structure based on the MHC haplotypes. Three haplotype clusters (A, B, C) segregated in these populations highlighted the relationship between the haplotypes in each cluster. The sequences from two clusters (B and C) overlapped, whereas the sequences from the third cluster (A) were very different. Overall, native breeds had high genetic diversity in the MHC-B region compared with the commercial breeds. This highlights their immune capabilities and genetic potential for resistance to many different pathogens.

Genetic diversity and population structure of indigenous chicken of Bangladesh using microsatellite markers.

OBJECTIVE: The objectives of this study were to investigate the genetic diversity, population structure and relatedness among the five chicken populations of Bangladesh using microsatellite markers. METHODS: A total of 161 individuals representing 5 chicken populations (non-descript Deshi [ND], naked neck [NN], hilly [HI], Aseel [AS], and red jungle fowl [JF]) were included in this study to investigate genetic diversity measures, population structure, genetic distance and phylogenetic relationships. Genotyping was performed using 16 selected polymorphic microsatellite markers distributed across 10 chromosomes. RESULTS: The average observed and expected heterozygosity, mean number of alleles and polymorphic information content were found to be 0.67±0.01, 0.70±0.01, 10.7 and 0.748, respectively in the studied populations. The estimated overall fixation index across the loci (F), heterozygote deficiency within (FIS) and among (FIT) chicken populations were 0.04±0.02, 0.05 and 0.16, respectively. Analysis of molecular variance analysis revealed 88.07% of the total genetic diversity was accounted for within population variation and the rest 11.93% was incurred with population differentiation (FST). The highest pairwise genetic distance (0.154) was found between ND and AS while the lowest distance was between JF and AS (0.084). Structure analysis depicted that the studied samples can be categorized into four distinct types or varieties (ΔK = 3.74) such as ND, NN, and HI where AS and JF clustered together as an admixed population. The Neighbor-Joining phylogenetic tree and discriminant analysis of principal component also showed close relatedness among three chicken varieties namely AS, HI, and JF. CONCLUSION: The results reflected that indigenous chicken of Bangladesh still possess rich genetic diversity but weak differentiation among the studied populations. This finding provides some important insight on genetic diversity measures that could support the designing and implementing of future breeding plans for indigenous chickens of Bangladesh.

Mitochondrial DNA variation and phylogeography of native Mongolian goats.

OBJECTIVE: Mongolia is one of a few countries that supports over 25 million goats, but genetic diversity, demographic history, and the origin of goat populations in Mongolia have not been well studied. This study was conducted to assess the genetic diversity, phylogenetic status and population structure of Mongolian native goats, as well as to discuss their origin together with other foreign breeds from different countries using hypervariable region 1 (HV1) in mtDNA. METHODS: In this study, we examined the genetic diversity and phylogenetic status of Mongolian native goat populations using a 452 base-pair long fragment of HVI of mitochondrial DNA from 174 individuals representing 12 populations. In addition, 329 previously published reference sequences from different regions were included in our phylogenetic analyses. RESULTS: Investigated native Mongolian goats displayed relatively high genetic diversities. After sequencing, we found a total of 109 polymorphic sites that defined 137 haplotypes among investigated populations. Of these, haplotype and nucleotide diversities of Mongolian goats were calculated as 0.997±0.001 and 0.0283±0.002, respectively. These haplotypes clearly clustered into four haplogroups (A, B, C, and D), with the predominance of haplogroup A (90.8%). Estimates of pairwise differences (Fst) and the analysis of molecular variance values among goat populations in Mongolia showed low genetic differentiation and weak geographical structure. In addition, Kazakh, Chinese (from Huanghuai and Leizhou), and Arabian (Turkish and Baladi breeds) goats had smaller genetic differentiation compared to Mongolian goats. CONCLUSION: In summary, we report novel information regarding genetic diversity, population structure, and origin of Mongolian goats. The findings obtained from this study reveal that abundant haplogroups (A to D) occur in goat populations in Mongolia, with high levels of haplotype and nucleotide diversity.

Population structure analysis of Yeonsan Ogye using microsatellite markers.

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.

Sequence characterization and polymorphism of melanocortin 1 receptor gene in some goat breeds with different coat color of Mongolia.

OBJECTIVE: Extension and Agouti loci play a key role for proportions of eumelanin and pheomelanin in determining coat color in several species, including goat. Mongolian goats exhibit diverse types of coat color phenotypes. In this study, investigation of the melanocortin 1 receptor (MC1R) coding region in different coat colors in Mongolian goats was performed to ascertain the presence of the extension allele. METHODS: A total of 105 goat samples representing three goat breeds were collected for this study from middle Mongolia. A 938 base pair (bp) long coding region of the MC1R gene was sequenced for three different breeds with different coat colors (Gobi Gurwan Saikhan: complete black, Zalaa Jinstiin Tsagaan: complete white, Mongolian native goat: admixture of different of coat colors). The genotypes of these goats were obtained from analyzing and comparing the sequencing results. RESULTS: A total of seven haplotypes defined by five substitution were identified. The five single nucleotide polymorphisms included two synonymous mutations (c.183C>T and c.489G>A) and three missense (non-synonymous) mutations (c.676A>G, c.748T>G, and c.770T>A). Comparison of genotypes frequencies of two common missense mutions using chi-sqaure (x2) test revealed significant differences between coat color groups (p<0.001). A logistic regression analysis additionally suggested highly significant association between genotypes and variation of black versus white uniform combination. Alternatively, most investigated goats (60.4%) belonged to H2 (TGAGT) haplotype. CONCLUSION: According to the findings obtained in this study on the investigated coat colors, mutations in MC1R gene may have the crucial role for determining eumelanin and pheomelanin phenotypes. Due to the complication of coat color phenotype, more detailed investigation needed.

Genetic diversity and phylogenetic relationship analyzed by microsatellite markers in eight Indonesian local duck populations.

OBJECTIVE: At least eight local duck breeds have been recognized and documented as national germplasm of Indonesia so far. It is necessary to genetically characterize the local duck breeds for aiding conservation and future improvement strategies. Thus, this study was carried out to assess genetic diversity and phylogenetic relationship of eight local duck populations of Indonesia using microsatellite markers. METHODS: In total, 240 individuals (30 individuals each population) from Alabio (AL), Bayang (BY), Magelang (MG), Mojosari (MJ), Pegagan (PG), Pitalah (PT), Rambon (RM), and Turi (TR) duck populations were genotyped using 22 microsatellite markers. RESULTS: The results showed a moderate level of genetic diversity among populations, with a total of 153 alleles detected over all loci and populations, ranging from 3 to 22 alleles per locus. Observed (Ho) and expected heterozygosity (He), as well as polymorphism information content over all loci and populations were 0.440, 0.566, and 0.513, respectively. Heterozygote deficiency in the overall populations (FIT = 0.237), was partly due to the heterozygote deficiency within populations (FIS = 0.114) and moderate level of genetic differentiation among populations (FST = 0.137). The most diverse population was MG (He = 0.545) and the least diverse population was AL (He = 0.368). The majority of populations were relatively in heterozygote deficiency (except AL), due to inbreeding. The genetic distances, phylogenetic trees, and principal coordinates analysis concluded that the populations can be grouped into two major clusters, resulting AL, MG, and MJ in one cluster separated from the remaining populations. CONCLUSION: The present study revealed a considerable genetic diversity of studied populations and thus, proper management strategies should be applied to preserve genetic diversity and prevent loss of alleles.

Identification of quantitative trait loci for the fatty acid composition in Korean native chicken.

OBJECTIVE: Fatty acid composition is one of the most important meat quality traits because it can contribute to functional, sensorial, and nutritional factors. In this study, quantitative trait locus (QTL) analyses for fatty acid composition traits were investigated in thigh and breast meat of Korean native chicken (KNC). METHODS: In total, 18 fatty acid composition traits were investigated from each meat sample using 83 parents, and 595 F1 chicks of 20 week old. Genotype assessment was performed using 171 informative DNA markers on 26 autosomes. The KNC linkage map was constructed by CRI-MAP software, which calculated genetic distances, with map orders between markers. The half-sib and full-sib QTL analyses were performed using GridQTL and SOLAR programs, respectively. RESULTS: In total, 30 QTLs (12 in the thigh and 18 in the breast meat) were detected by the half-sib analysis and 7 QTLs (3 in the thigh and 4 in the breast meat) were identified by the full-sib analysis. CONCLUSION: With further verification of the QTL regions using additional markers and positional candidate gene studies, these results can provide valuable information for determining causative mutations affecting the fatty acid composition of KNC meat. Moreover, these findings may aid in the selection of birds with favorable fatty acid composition traits.

Estimation of heritability and genetic correlation of body weight gain and growth curve parameters in Korean native chicken.

OBJECTIVE: This study estimated the genetic parameters for body weight gain and growth curve parameter traits in Korean native chicken (KNC). METHODS: A total of 585 F1 chickens were used along with 88 of their F0 birds. Body weights were measured every 2 weeks from hatching to 20 weeks of age to measure weight gain at 2-week intervals. For each individual, a logistic growth curve model was fitted to the longitudinal growth dataset to obtain three growth curve parameters (α, asymptotic final body weight; ß, inflection point; and γ, constant scale that was proportional to the overall growth rate). Genetic parameters were estimated based on the linear-mixed model using a restricted maximum likelihood method. RESULTS: Heritability estimates of body weight gain traits were low to high (0.057 to 0.458). Heritability estimates for α, ß, and γ were 0.211±0.08, 0.249±0.09, and 0.095±0.06, respectively. Both genetic and phenotypic correlations between weight gain traits ranged from -0.527 to 0.993. Genetic and phenotypic correlation between the growth curve parameters and weight gain traits ranged from -0.968 to 0.987. CONCLUSION: Based on the results of this study population, we suggest that the KNC could be used for selective breeding between 6 and 8 weeks of age to enhance the overall genetic improvement of growth traits. After validation of these results in independent studies, these findings will be useful for further optimization of breeding programs for KNC.

POU class 1 homeobox 1 gene polymorphisms associated with growth traits in Korean native chicken.

OBJECTIVE: POU class 1 homeobox 1 (POU1F1) mediates growth hormone expression and activity by altering transcription, eventually resulting in growth rate variations. Therefore, we aimed to identify chicken POU1F1 polymorphisms and evaluate the association between single nucleotide polymorphisms (SNPs) and growth-related traits, and logistic growth curve parameter traits (α, ß, and γ). METHODS: Three SNPs (M_1 to M_3) were used to genotype 585 F1 and 88 F0 birds from five Korean native chicken lines using a polymerase chain reaction-restriction fragment length polymorphism method. RESULTS: Single marker analyses and traits association analyses showed that M_2 was significantly associated with body weight at two weeks, weight gain from hatch to 2 weeks, and weight gain from 16 to 18 weeks (p<0.05). M_3 was significantly associated with weight gain from 14 to 16 weeks and from 16 to 18 weeks, and asymptotic body weight (α) (p<0.05). No traits were associated with M_1. The POU1F1 haplogroups were significantly associated with weight gain from 14 to 16 weeks (p = 0.020). Linkage disequilibrium test and Haploview analysis shown one main haploblock between M_2 and M_3 SNP. CONCLUSION: Thus, POU1F1 significantly affects the growth of Korean native chickens and their growth curve traits.