Application of genomic big data to analyze the genetic diversity and population structure of Korean domestic chickens.

Genetic diversity analysis is crucial for maintaining and managing genetic resources. Several studies have examined the genetic diversity of Korean domestic chicken (KDC) populations using microsatellite markers, but it is difficult to capture the characteristics of the whole genome in this manner. Hence, this study analyzed the genetic diversity of several KDC populations using high-density single nucleotide polymorphism (SNP) genotype data. We examined 935 birds from 21 KDC populations, including indigenous and adapted Korean native chicken (KNC), Hyunin and Jeju KDC, and Hanhyup commercial KDC populations. A total of 212,420 SNPs of 21 KDC populations were used for calculating genetic distances and fixation index, and for ADMIXTURE analysis. As a result of the analysis, the indigenous KNC groups were genetically closer and more fixed than the other groups. Furthermore, Hyunin and Jeju KDC were similar to the indigenous KNC. In comparison, adapted KNC and Hanhyup KDC populations derived from the same original species were genetically close to each other, but had different genetic structures from the others. In conclusion, this study suggests that continuous evaluation and management are required to prevent a loss of genetic diversity in each group. Basic genetic information is provided that can be used to improve breeds quickly by utilizing the various characteristics of native chickens.

Genome-wide association study revealed the genomic regions associated with skin pigmentation in an Ogye x White Leghorn F2 chicken population.

Skin color in chickens is an economically important trait that determines the first impression of a consumer toward a broiler and can ultimately affect consumer choice in the market. Therefore, identification of genomic regions associated with skin color is crucial for increasing the sales value of chickens. Although previous studies have attempted to reveal the genetic markers associated with the skin coloration in chickens, most were limited to investigations of candidate genes, such as melanin-related genes, and focused on case/control studies based on a single or small population. In this study, we performed a genome-wide association study (GWAS) on 770 F2 intercrosses produced by an experimental population of 2 chicken breeds, namely Ogye and White Leghorns, with different skin colors. The GWAS demonstrated that the L* value among the 3 skin color traits is highly heritable, and the genomic regions located on 2 chromosomes (20 and Z) were detected to harbor SNPs significantly associated with the skin color trait, accounting for most of the total genetic variance. Particular genomic regions spanning a ∼2.94 Mb region on GGA Z and a ∼3.58 Mb region on GGA 20 were significantly associated with skin color traits, and in these regions, certain candidate genes, including MTAP, FEM1C, GNAS, and EDN3, were found. Our findings could help elucidate the genetic mechanisms underlying chicken skin pigmentation. Furthermore, the candidate genes can be used to provide a valuable breeding strategy for the selection of specific chicken breeds with ideal skin coloration.

Estimating genetic diversity and population structure of 22 chicken breeds in Asia using microsatellite markers.

OBJECTIVE: Estimating the genetic diversity and structures, both within and among chicken breeds, is critical for the identification and conservation of valuable genetic resources. In chickens, microsatellite (MS) marker polymorphisms have previously been widely used to evaluate these distinctions. Our objective was to analyze the genetic diversity and relationships among 22 chicken breeds in Asia based on allelic frequencies. METHODS: We used 469 genomic DNA samples from 22 chicken breeds from eight Asian countries (South Korea, KNG, KNB, KNR, KNW, KNY, KNO; Laos, LYO, LCH, LBB, LOU; Indonesia, INK, INS, ING; Vietnam, VTN, VNH; Mongolia, MGN; Kyrgyzstan, KGPS; Nepal, NPS; Sri Lanka, SBC) and three imported breeds (RIR, Rhode Island Red; WLG, White Leghorn; CON, Cornish). Their genetic diversity and phylogenetic relationships were analyzed using 20 MS markers. RESULTS: In total, 193 alleles were observed across all 20 MS markers, and the number of alleles ranged from 3 (MCW0103) to 20 (LEI0192) with a mean of 9.7 overall. The NPS breed had the highest expected heterozygosity (Hexp, 0.718±0.027) and polymorphism information content (PIC, 0.663±0.030). Additionally, the observed heterozygosity (Hobs) was highest in LCH (0.690±0.039), whereas WLG showed the lowest Hexp (0.372±0.055), Hobs (0.384±0.019), and PIC (0.325±0.049). Nei's DA genetic distance was the closest between VTN and VNH (0.086), and farthest between KNG and MGN (0.503). Principal coordinate analysis showed similar results to the phylogenetic analysis, and three axes explained 56.2% of the variance (axis 1, 19.17%; 2, 18.92%; 3, 18.11%). STRUCTURE analysis revealed that the 22 chicken breeds should be divided into 20 clusters, based on the highest ΔK value (46.92). CONCLUSION: This study provides a basis for future genetic variation studies and the development of conservation strategies for 22 chicken breeds in Asia.

Haplogroup Classification of Korean Cattle Breeds Based on Sequence Variations of mtDNA Control Region.

Many studies have reported the frequency and distribution of haplogroups among various cattle breeds for verification of their origins and genetic diversity. In this study, 318 complete sequences of the mtDNA control region from four Korean cattle breeds were used for haplogroup classification. 71 polymorphic sites and 66 haplotypes were found in these sequences. Consistent with the genetic patterns in previous reports, four haplogroups (T1, T2, T3, and T4) were identified in Korean cattle breeds. In addition, T1a, T3a, and T3b sub-haplogroups were classified. In the phylogenetic tree, each haplogroup formed an independent cluster. The frequencies of T3, T4, T1 (containing T1a), and T2 were 66%, 16%, 10%, and 8%, respectively. Especially, the T1 haplogroup contained only one haplotype and a sample. All four haplogroups were found in Chikso, Jeju black and Hanwoo. However, only the T3 and T4 haplogroups appeared in Heugu, and most Chikso populations showed a partial of four haplogroups. These results will be useful for stable conservation and efficient management of Korean cattle breeds.

Identification of Single Nucleotide Polymorphism Marker and Association Analysis of Marbling Score in Fas Gene of Hanwoo.

The Fas (APO-1, TNFRSF6) gene known as a member of the tumor necrosis factor receptor superfamily was selected for DNA marker development in Korean cattle. It is a cell membrane protein and mediates programmed cell death (apoptosis). We discovered single nucleotide polymorphisms (SNPs) within Fas gene in order to develop novel DNA markers related to economical traits at the genomic level. The sequences of whole exon and 1 kb range of both front and back of the gene were determined by direct-sequencing methods using 24 cattle. A total of 55 SNPs were discovered and we selected 31 common polymorphic sites considering their allele frequencies, haplotype-tagging status and linkage disequilibrium (LD) for genotyping in larger-scale subjects. The SNPs were confirmed genotype through the SNaPshot method (n = 274) and were examined for a possible genetic association between Fas polymorphisms and marbling score. So, the SNPs that were identified significant are g.30256G>C, g.31474C>A, g.31940A>G, and g.32982G>A. These results suggest that SNPs of Fas gene were associated with intramuscular fat content of meat quality traits in Korean cattle.

Mutations in calpastatin and µ-calpain are associated with meat tenderness, flavor and juiciness in Hanwoo (Korean cattle): molecular modeling of the effects of substitutions in the calpastatin/µ-calpain complex.

The objective of this study was to evaluate the effects of seven single nucleotide polymorphisms (SNPs) in Calpain 1 and Calpastatin genes previously associated with meat tenderness attributes in other cattle breeds in Korean Hanwoo cattle. The Hanwoo resource population was used to study association of 7 SNPs with beef tenderness, flavor, juiciness, intramuscular fat and shear force. In this association study, CAST:c.182A>G (+0.14, P=0.04) and CAST:c.1985G>C (-0.12, P=0.02) had significant effects on juiciness, but no effects on other traits. In contrast, CAPN1:c.1589G>A was associated with meat tenderness (P=0.01) and juiciness (P=0.04). The CAPN1:c.1589G>A (Val530Ile) SNP marker displayed significant effect on the meat tenderness score which is strongly supported by molecular modeling of the CAPN1:c.1589G>A (Val530Ile) variant that inhibits CAST protein from binding more strongly than the wild-type protein, which may explain its effect on meat tenderness.

Changes in expression of insulin signaling pathway genes by dietary fat source in growing-finishing pigs.

This study investigated changes in gene expression by dietary fat source, i.e., beef tallow, soybean oil, olive oil, and coconut oil (each 3% in feed), in both male and female growing-finishing pigs. Real-time PCR was conducted on seven genes (insulin receptor; INSR, insulin receptor substrate; IRS, phosphatidylinositol (3,4,5)-triphosphate; PIP3, 3-phosphoinositide-dependent protein kinase-1; PDK1, protein kinase B; Akt, forkhead box protein O1; FOXO1 and cGMP-inhibited 3', 5'-cyclic phosphodiesterase; PDE3) located upstream of the insulin signaling pathway in the longissimus dorsi muscle (LM) of pigs. The INSR, IRS, PIP3, and PDE3 genes showed significantly differential expression in barrow pigs. Expression of the PIP3 and FOXO1 genes was significantly different among the four dietary groups in gilt pigs. In particular, the PIP3 gene showed the opposite expression pattern between barrow and gilt pigs. These results show that dietary fat source affected patterns of gene expression according to animal gender. Further, the results indicate that the type of dietary fat affects insulin signaling-related gene expression in the LM of pigs. These results can be applied to livestock production by promoting the use of discriminatory feed supplies.

ERα/E2 signaling suppresses the expression of steroidogenic enzyme genes via cross-talk with orphan nuclear receptor Nur77 in the testes.

Estrogen receptor alpha (ERα) has been reported to affect steroidogenesis in testicular Leydig cells, but its molecular mechanism remains unclear. Here, we investigate the effect of estrogen and ERα on Nur77, a major transcription factor that regulates the expression of steroidogenic enzyme genes. In MA-10 Leydig cells, estradiol (E2) treatment, and interestingly ERα overexpression, suppressed the cAMP-induced and Nur77-activated promoter activity of steroidogenic enzyme genes via the suppression of Nur77 transactivation. ERα physically interacted with Nur77 and inhibited its DNA binding activity. In addition, ERα/E2 signaling decreased Nur77 protein levels. Consistent with the above results, the testicular testosterone level was higher in Leydig cell-specific ERα knock-out mice (ERα(flox/flox)Cyp17iCre) than in wild-type mice (ERα(flox/flox)). Taken together, these results suggest that ERα/E2 signaling controls the Nur77-mediated expression of steroidogenic enzyme genes in Leydig cells. These findings may provide a mechanistic explanation for the local regulation of testicular steroidogenesis by estrogenic compounds and ERα.

Porcine LMNA Is a Positional Candidate Gene Associated with Growth and Fat Deposition.

Crosses between Korean and Landrace pigs have revealed a large quantitative trait loci (QTL) region for fat deposition in a region (89 cM) of porcine chromosome 4 (SSC4). To more finely map this QTL region and identify candidate genes for this trait, comparative mapping of pig and human chromosomes was performed in the present study. A region in the human genome that corresponds to the porcine QTL region was identified in HSA1q21. Furthermore, the LMNA gene, which is tightly associated with fat augmentation in humans, was localized to this region. Radiation hybrid (RH) mapping using a Sus scrofa RH panel localized LMNA to a region of 90.3 cM in the porcine genome, distinct from microsatellite marker S0214 (87.3 cM). Two-point analysis showed that LMNA was linked to S0214, SW1996, and S0073 on SSC4 with logarithm (base 10) of odds scores of 20.98, 17.78, and 16.73, respectively. To clone the porcine LMNA gene and to delineate the genomic structure and sequences, including the 3'untranslated region (UTR), rapid amplification of cDNA ends was performed. The coding sequence of porcine LMNA consisted of 1,719 bp, flanked by a 5'UTR and a 3'UTR. Two synonymous single nucleotide polymorphisms (SNPs) were identified in exons 3 and 7. Association tests showed that the SNP located in exon 3 (A193A) was significantly associated with weight at 30 wks (p<0.01) and crude fat content (p<0.05). This association suggests that SNPs located in LMNA could be used for marker-assisted selection in pigs.

Molecular cloning and expression of steroidogenic acute regulatory protein from bullfrog (Rana catesbeiana).

Steroidogenic acute regulatory protein (StAR) transfers cholesterol from the outer mitochondrial membrane to the inner membrane where the cytochrome P450 side chain cleavage enzyme (P450scc) resides. This process is the rate-limiting step in steroidogenesis. StAR cDNAs have been cloned and characterized from a range of different species. To investigate the role of StAR in the amphibian system, we cloned a full-length StAR cDNA from bullfrog (Rana catesbeiana) using reverse transcription polymerase chain reaction (RT-PCR) in conjunction with rapid amplification of cDNA ends (RACE). The putative full-length bullfrog StAR (bfStAR) cDNA was 1862 base pairs (bp) in length, and the longest open reading frame (ORF) encoded a protein of 284 amino acids. Amino acid sequence comparison showed that amphibian StAR has a high degree of sequence identity, ranging from 62% to 98%, with StAR proteins of other species. Similar to other species, bfStAR contained two conserved domains, the mitochondrial targeting domain and cholesterol-binding domain, in the N-terminus and C-terminus of the protein, respectively. Northern blot analysis and RT-PCR indicated that StAR mRNA is expressed in the gonads and adrenal gland. Transfection of green monkey kidney (COS-1) cells with an expression construct for bfStAR revealed that it encoded 34 and 27kDa proteins that were recognized by antiserum raised against the human StAR-related lipid transfer (START) domain.

Effects of butyltin compounds on follicular steroidogenesis in the bullfrog(Rana catesbeiana).

The effects of butyltin compounds on follicular steroidogenesis in amphibians were examined using ovarian follicles of Rana catesbeiana. Isolated follicles were cultured for 18h in the presence and absence of frog pituitary homogenate (FPH) or various steroid precursors, and the steroid levels in the follicles or culture media were measured by radioimmunoassay (RIA). Among the butyltin compounds, tributyltin (TBT) strongly inhibited the FPH-induced synthesis of pregnenolone (P(5)), progesterone (P(4)) and testosterone (T). It also inhibited the conversion of P(5)-P(4) and T to estradiol-17ß(E(2)) and it partially suppressed the conversion of androstenedione (AD) to T, but not P(4) to 17α-hydroxyprogesterone (17α-OHP(4)). A high concentration of dibutyltin (DBT) also inhibited steroidogenesis by the follicles while monobutyltin (MBT) and tetrabutylin (TeBT) had no effect. These results suggest that the initial step of steroidogenesis (P(5) synthesis) and enzymes such as 3ß-HSD, 17ß-HSD and aromatase are inhibited by TBT or DBT. However, 17α-hydroxylase was not suppressed by TBT or the other butyltin compounds.