• Genes and Genomics
• /
• v.40 no.12
• /
• pp.1373-1382
• /
• 2018

Hanwoo and Chikso are classified as Korean native cattle breeds that are currently registered with the Food and Agriculture Organization. However, there is still a lack of genomic studies to compare Hanwoo to Chikso populations. The objective of this study was to perform genome-wide analysis of Hanwoo and Chikso populations, investigating the genetic relationships between these two populations. We genotyped a total of 319 cattle including 214 Hanwoo and 105 Chikso sampled from Gangwon Province Livestock Technology Research Institute, using the Illumina Bovine SNP50K Beadchip. After performing quality control on the initially generated datasets, we assessed linkage disequilibrium patterns for all the possible SNP pairs within 1 Mb apart. Overall, average $r^2$ values in Hanwoo (0.048) were lower than Chikso (0.074) population. The genetic relationship between the populations was further assured by the principal component analysis, exhibiting clear clusters in each of the Hanwoo and Chikso populations, respectively. Overall heterozygosity for Hanwoo (0.359) was slightly higher than Chikso (0.345) and inbreeding coefficient was also a bit higher in Hanwoo (-0.015) than Chikso (-0.035). The average $F_{ST}$ value was 0.036 between Hanwoo and Chikso, indicating little genetic differentiation between those two breeds. Furthermore, we found potential selection signatures including LRP1B and NTRK2 genes that might be implicated with meat and reproductive traits in cattle. In this study, the results showed that both Hanwoo and Chikso populations were not under severe level of inbreeding. Although the principal component analysis exhibited clear clusters in each of the populations, we did not see any clear evidence that those two populations are highly differentiated each other.

• The Korean Journal of Physiology and Pharmacology
• /
• v.23 no.5
• /
• pp.367-379
• /
• 2019

Although atopic dermatitis (AD) is known to be a representative skin disorder, it also affects the systemic immune response. In a recent study, myoblasts were shown to be involved in the immune regulation, but the roles of muscle cells in AD are poorly understood. We aimed to identify the relationship between mitochondria and atopy by genome-wide analysis of skeletal muscles in mice. We induced AD-like symptoms using house dust mite (HDM) extract in NC/Nga mice. The transcriptional profiles of the untreated group and HDM-induced AD-like group were analyzed and compared using microarray, differentially expressed gene and functional pathway analyses, and protein interaction network construction. Our microarray analysis demonstrated that immune response-, calcium handling-, and mitochondrial metabolism-related genes were differentially expressed. In the Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology pathway analyses, immune response pathways involved in cytokine interaction, nuclear factor-kappa B, and T-cell receptor signaling, calcium handling pathways, and mitochondria metabolism pathways involved in the citrate cycle were significantly upregulated. In protein interaction network analysis, chemokine family-, muscle contraction process-, and immune response-related genes were identified as hub genes with many interactions. In addition, mitochondrial pathways involved in calcium signaling, cardiac muscle contraction, tricarboxylic acid cycle, oxidation-reduction process, and calcium-mediated signaling were significantly stimulated in KEGG and Gene Ontology analyses. Our results provide a comprehensive understanding of the genome-wide transcriptional changes of HDM-induced AD-like symptoms and the indicated genes that could be used as AD clinical biomarkers.

• Asian-Australasian Journal of Animal Sciences
• /
• v.32 no.4
• /
• pp.494-500
• /
• 2019

Objective: Feed consumption contributes a large percentage for total production costs in the poultry industry. Detecting genes associated with feeding traits will be of benefit to improve our understanding of the molecular determinants for feed efficiency. The objective of this study was to identify candidate genes associated with feed conversion ratio (FCR) via genomewide association study (GWAS) using sequence data imputed from single nucleotide polymorphism (SNP) panel in a Chinese indigenous chicken population. Methods: A total of 435 Chinese indigenous chickens were phenotyped for FCR and were genotyped using a 600K SNP genotyping array. Twenty-four birds were selected for sequencing, and the 600K SNP panel data were imputed to whole sequence data with the 24 birds as the reference. The GWAS were performed with GEMMA software. Results: After quality control, 8,626,020 SNPs were used for sequence based GWAS, in which ten significant genomic regions were detected to be associated with FCR. Ten candidate genes, ubiquitin specific peptidase 44, leukotriene A4 hydrolase, ETS transcription factor, R-spondin 2, inhibitor of apoptosis protein 3, sosondowah ankyrin repeat domain family member D, calmodulin regulated spectrin associated protein family member 2, zinc finger and BTB domain containing 41, potassium sodium-activated channel subfamily T member 2, and member of RAS oncogene family were annotated. Several of them were within or near the reported FCR quantitative trait loci, and others were newly reported. Conclusion: Results from this study provide valuable prior information on chicken genomic breeding programs, and potentially improve our understanding of the molecular mechanism for feeding traits.

• Asian-Australasian Journal of Animal Sciences
• /
• v.32 no.7
• /
• pp.913-921
• /
• 2019

Objective: The objectives of this study were to compare identified informative regions through two genome-wide association study (GWAS) approaches and determine the accuracy and bias of the direct genomic value (DGV) for milk production traits in Korean Holstein cattle, using two genomic prediction approaches: single-step genomic best linear unbiased prediction (ss-GBLUP) and Bayesian Bayes-B. Methods: Records on production traits such as adjusted 305-day milk (MY305), fat (FY305), and protein (PY305) yields were collected from 265,271 first parity cows. After quality control, 50,765 single-nucleotide polymorphic genotypes were available for analysis. In GWAS for ss-GBLUP (ssGWAS) and Bayes-B (BayesGWAS), the proportion of genetic variance for each 1-Mb genomic window was calculated and used to identify informative genomic regions. Accuracy of the DGV was estimated by a five-fold cross-validation with random clustering. As a measure of accuracy for DGV, we also assessed the correlation between DGV and deregressed-estimated breeding value (DEBV). The bias of DGV for each method was obtained by determining regression coefficients. Results: A total of nine and five significant windows (1 Mb) were identified for MY305 using ssGWAS and BayesGWAS, respectively. Using ssGWAS and BayesGWAS, we also detected multiple significant regions for FY305 (12 and 7) and PY305 (14 and 2), respectively. Both single-step DGV and Bayes DGV also showed somewhat moderate accuracy ranges for MY305 (0.32 to 0.34), FY305 (0.37 to 0.39), and PY305 (0.35 to 0.36) traits, respectively. The mean biases of DGVs determined using the single-step and Bayesian methods were $1.50{\pm}0.21$ and $1.18{\pm}0.26$ for MY305, $1.75{\pm}0.33$ and $1.14{\pm}0.20$ for FY305, and $1.59{\pm}0.20$ and $1.14{\pm}0.15$ for PY305, respectively. Conclusion: From the bias perspective, we believe that genomic selection based on the application of Bayesian approaches would be more suitable than application of ss-GBLUP in Korean Holstein populations.

• Proceedings of the Korean Society of Crop Science Conference
• /
• 2022.10a
• /
• pp.202-202
• /
• 2022

Recently days, soybean production in paddy field is increasing, from 4,422 ha in 2016 to 10,658 ha in 2021 in Korea. It is easy for Phytophthora stem and root rot (PSR) occurring in paddy field condition, when it is poorly drained soils with a high clay content, and temporary flooding and ponding. Therefore PSR resistant soybean cultivar is required. The objective of this study is to identify QTL region and candidate genes relating to PSR resistance of the race in main soybean cultivation area in Korea. 210 soybean materials including cultivars and germplasm were used for inoculation and genome-wide association study (GWAS). Inoculation was conducted using stem-scar method with 2 replications in 2-year for the race 3053 from Kimje and 3617 from Andong. 210 materials were genotyped with Soya SNP 180K chip, and structure analysis and association mapping were conducted with QTLMAX V2. The results of inoculation showed that survival ratio ranged from 0% to 96.7% and mean 9.7% for 3053 and ranged from 0% to 100% and mean 7.6% for 3617. Structure analysis showed linkage disequillibrium (LD) was decayed below r²=0.5 at 335kb of SNP distance. Significant SNPs (LOD>7.0) were identified in Chr 1, 2, 3, 4, 5, 11, 14, 15 for 3053 and Chr 1, 2, 3, 7, 10, 14 for 3617. Especially, LD blocks (AX-90455181;15,056,628bp~AX-90475572;15,298,872bp) in Chr 2 for 3053 and 3067 were duplicated. 29 genes were identified on these genetic regions including Glyma.02gl47000 relating to ribosome recycling factor and defense response to fungus in Soybase.

• Molecules and Cells
• /
• v.41 no.11
• /
• pp.953-963
• /
• 2018

The stepwise development of T cells from a multipotent precursor is guided by diverse mechanisms, including interactions among lineage-specific transcription factors (TFs) and epigenetic changes, such as DNA methylation and hydroxymethylation, which play crucial roles in mammalian development and lineage commitment. To elucidate the transcriptional networks and epigenetic mechanisms underlying T-cell lineage commitment, we investigated genome-wide changes in gene expression, DNA methylation and hydroxymethylation among populations representing five successive stages of T-cell development (DN3, DN4, DP, $CD4^+$, and $CD8^+$) by performing RNA-seq, MBD-seq and hMeDIP-seq, respectively. The most significant changes in the transcriptomes and epigenomes occurred during the DN4 to DP transition. During the DP stage, many genes involved in chromatin modification were up-regulated and exhibited dramatic changes in DNA hydroxymethylation. We also observed 436 alternative splicing events, and approximately 57% (252) of these events occurred during the DP stage. Many stage-specific, differentially methylated regions were observed near the stage-specific, differentially expressed genes. The dynamic changes in DNA methylation and hydroxymethylation were associated with the recruitment of stage-specific TFs. We elucidated interactive networks comprising TFs, chromatin modifiers, and DNA methylation and hope that this study provides a framework for the understanding of the molecular networks underlying T-cell lineage commitment.

• Asian-Australasian Journal of Animal Sciences
• /
• v.28 no.7
• /
• pp.926-935
• /
• 2015

The efficiency of genome-wide association analysis (GWAS) depends on power of detection for quantitative trait loci (QTL) and precision for QTL mapping. In this study, three different strategies for GWAS were applied to detect QTL for carcass quality traits in the Korean cattle, Hanwoo; a linkage disequilibrium single locus regression method (LDRM), a combined linkage and linkage disequilibrium analysis (LDLA) and a $BayesC{\pi}$ approach. The phenotypes of 486 steers were collected for weaning weight (WWT), yearling weight (YWT), carcass weight (CWT), backfat thickness (BFT), longissimus dorsi muscle area, and marbling score (Marb). Also the genotype data for the steers and their sires were scored with the Illumina bovine 50K single nucleotide polymorphism (SNP) chips. For the two former GWAS methods, threshold values were set at false discovery rate <0.01 on a chromosome-wide level, while a cut-off threshold value was set in the latter model, such that the top five windows, each of which comprised 10 adjacent SNPs, were chosen with significant variation for the phenotype. Four major additive QTL from these three methods had high concordance found in 64.1 to 64.9Mb for Bos taurus autosome (BTA) 7 for WWT, 24.3 to 25.4Mb for BTA14 for CWT, 0.5 to 1.5Mb for BTA6 for BFT and 26.3 to 33.4Mb for BTA29 for BFT. Several candidate genes (i.e. glutamate receptor, ionotropic, ampa 1 [GRIA1], family with sequence similarity 110, member B [FAM110B], and thymocyte selection-associated high mobility group box [TOX]) may be identified close to these QTL. Our result suggests that the use of different linkage disequilibrium mapping approaches can provide more reliable chromosome regions to further pinpoint DNA makers or causative genes in these regions.

• Proceedings of the Korean Society of Crop Science Conference
• /
• 2017.06a
• /
• pp.145-145
• /
• 2017

Foxglove aphid, Aulacorthum solani (Kaltenbach), is a Hemipteran insect that infected a wide variety of plants worldwide and caused serious yield losses in crops. The objective of this study was to identify the putative genes to foxglove aphid resistance in wild soybean, PI 366121 (Glycine soja Sieb. and Zucc.). One hundred and forty-one F4:8 recombinant inbred lines developed from a cross between susceptible variety, Williams 82 and foxglove aphid resistance wild soybean, PI 366121 were used. The two type of resistance response, antibiosis and antixenosis resistance were evaluated through choice and no-choice test, graded by the degree of total plant damage and primary infestation leaf damage; a genome-wide molecular linkage map was constructed with 29,898 single-nucleotide polymorphism markers utilizing a Axiom(R) 180K soyaSNP array. Using inclusive composite interval mapping analysis for foxglove aphid resistance, one major candidate QTL on chromosome 7 was identified. The major QTL on chromosome 7 showed both antixenosis and antibiosis resistance responses. The newly identified major QTL was consistent with previously reported QTL, Raso2, which showed around 5 times narrow down interval range with 8 candidate genes. Furthermore, total 1,115 soybean varieties including Glycine soja and Glycine max were exposed to germplasm screening, and 31 varieties, which showed significant antibiosis type foxglove aphid resistance were identified. This result could be useful in breeding for new foxglove aphid resistant soybean cultivars and developing novel insecticides.

• Journal of Microbiology and Biotechnology
• /
• v.25 no.4
• /
• pp.511-520
• /
• 2015

Triclosan, the widely used biocide, specifically targets enoyl-acyl carrier protein reductase (ENR) in the bacterial fatty acid synthesis system. Although the fish pathogen Aeromonas salmonicida subsp. salmonicida exhibits triclosan resistance, the nature of this resistance has not been elucidated. Here, we aimed to characterize the triclosan resistance of A. salmonicida subsp. salmonicida causing furunculosis. The fosmid library of triclosan-resistant A. salmonicida subsp. salmonicida was constructed to select a fosmid clone showing triclosan resistance. With the fosmid clone showing triclosan resistance, a subsequent secondary library search resulted in the selection of subclone pTSR-1. DNA sequence analysis of pTSR-1 revealed the presence of a chromosomal-borne fabV-encoding ENR homolog. The ENR of A. salmonicida (FabVas) exhibited significant homology with previously known FabV, including the catalytic domain YX₍₈₎K. fabVas introduction into E. coli dramatically increased its resistance to triclosan. Heterologous expression of FabVas might functionally replace the triclosan-sensitive FabI in vivo to confer E. coli with triclosan resistance. A genome-wide search for fabVas homologs revealed the presence of an additional fabV gene (fabVas2) paralog in A. salmonicida strains and the fabVas orthologs from other gram-negative fish pathogens. Both of the potential FabV ENRs expressed similarly with or without triclosan supplement. This is the first report about the presence of two potential FabV ENRs in a single pathogenic bacterium. Our result suggests that triclosan-resistant ENRs are widely distributed in various bacteria in nature, and the wide use of this biocide can spread these triclosan-tolerant ENRs among fish pathogens and other pathogenic bacteria.

• Tuberculosis and Respiratory Diseases
• /
• v.82 no.2
• /
• pp.126-132
• /
• 2019

Background: The development of lung cancer results from the interaction between genetic mutations and dynamic epigenetic alterations, although the exact mechanisms are not completely understood. Changes in DNA methylation may be a promising biomarker for early detection and prognosis of lung cancer. We evaluated the serial changes in genome-wide DNA methylation patterns in blood samples of lung cancer patients. Methods: Blood samples were obtained for three consecutive years from three patients (2 years before, 1 year before, and after lung cancer detection) and from three control subjects (without lung cancer). We used the MethylationEPIC BeadChip method, which covers the 850,000 bp cytosine-phosphate-guanine (CpG) site, to conduct an epigenome-wide analysis. Significant differentially methylated regions (DMRs) were identified using p-values <0.05 in a correlation test identifying serial methylation changes and serial increase or decrease in ${\beta}$ value above 0.1 for three consecutive years. Results: We found three significant CpG sites with differentially methylated ${\beta}$ values and 7,105 CpG sites with significant correlation from control patients without lung cancer. However, there were no significant DMRs. In contrast, we found 11 significant CpG sites with differentially methylated ${\beta}$ values and 10,562 CpG sites with significant correlation from patients with lung cancer. There were two significant DMRs: cg21126229 (RNF212) and cg27098574 (BCAR1). Conclusion: This study revealed DNA methylation changes that might be implicated in lung cancer development. The DNA methylation changes may be the possible candidate target regions for the early detection and prevention of lung cancer.