Dajeong Lim

Whole-genome resequencing of Hanwoo (Korean cattle) and insight into regions of homozygosity

BackgroundHanwoo (Korean cattle), which originated from natural crossbreeding between taurine and zebu cattle, migrated to the Korean peninsula through North China. Hanwoo were raised as draft animals until the 1970s without the introduction of foreign germplasm. Since 1979, Hanwoo has been bred as beef cattle. Genetic variation was analyzed by whole-genome deep resequencing of a Hanwoo bull. The Hanwoo genome was compared to that of two other breeds, Black Angus and Holstein, and genes within regions of homozygosity were investigated to elucidate the genetic and genomic characteristics of Hanwoo.ResultsThe Hanwoo bull genome was sequenced to 45.6-fold coverage using the ABI SOLiD system. In total, 4.7 million single-nucleotide polymorphisms and 0.4 million small indels were identified by comparison with the Btau4.0 reference assembly. Of the total number of SNPs and indels, 58% and 87%, respectively, were novel. The overall genotype concordance between the SNPs and BovineSNP50 BeadChip data was 96.4%. Of 1.6 million genetic differences in Hanwoo, approximately 25,000 non-synonymous SNPs, splice-site variants, and coding indels (NS/SS/Is) were detected in 8,360 genes. Among 1,045 genes containing reliable specific NS/SS/Is in Hanwoo, 109 genes contained more than one novel damaging NS/SS/I. Of the genes containing NS/SS/Is, 610 genes were assigned as trait-associated genes. Moreover, 16, 78, and 51 regions of homozygosity (ROHs) were detected in Hanwoo, Black Angus, and Holstein, respectively. ‘Regulation of actin filament length’ was revealed as a significant gene ontology term and 25 trait-associated genes for meat quality and disease resistance were found in 753 genes that resided in the ROHs of Hanwoo. In Hanwoo, 43 genes were located in common ROHs between whole-genome resequencing and SNP chips in BTA2, 10, and 13 coincided with quantitative trait loci for meat fat traits. In addition, the common ROHs in BTA2 and 16 were in agreement between Hanwoo and Black Angus.ConclusionsWe identified 4.7 million SNPs and 0.4 million small indels by whole-genome resequencing of a Hanwoo bull. Approximately 25,000 non-synonymous SNPs, splice-site variants, and coding indels (NS/SS/Is) were detected in 8,360 genes. Additionally, we found 25 trait-associated genes for meat quality and disease resistance among 753 genes that resided in the ROHs of Hanwoo. These findings will provide useful genomic information for identifying genes or casual mutations associated with economically important traits in cattle.

Neuronal Genes for Subcutaneous Fat Thickness in Human and Pig Are Identified by Local Genomic Sequencing and Combined SNP Association Study

Obesity represents a major global public health problem that increases the risk for cardiovascular or metabolic disease. The pigs represent an exceptional biomedical model related to energy metabolism and obesity in humans. To pinpoint causal genetic factors for a common form of obesity, we conducted local genomic de novo sequencing, 18.2 Mb, of a porcine QTL region affecting fatness traits, and carried out SNP association studies for backfat thickness and intramuscular fat content in pigs. In order to relate the association studies in pigs to human obesity, we performed a targeted genome wide association study for subcutaneous fat thickness in a cohort population of 8,842 Korean individuals. These combined association studies in human and pig revealed a significant SNP located in a gene family with sequence similarity 73, member A (FAM73A) associated with subscapular skin-fold thickness in humans (rs4121165, GC-corrected p-value  = 0.0000175) and with backfat thickness in pigs (ASGA0029495, p-value  = 0.000031). Our combined association studies also suggest that eight neuronal genes are responsible for subcutaneous fat thickness: NEGR1, SLC44A5, PDE4B, LPHN2, ELTD1, ST6GALNAC3, ST6GALNAC5, and TTLL7. These results provide strong support for a major involvement of the CNS in the genetic predisposition to a common form of obesity.

Generation and analysis of large-scale expressed sequence tags (ESTs) from a full-length enriched cDNA library of porcine backfat tissue

BackgroundGenome research in farm animals will expand our basic knowledge of the genetic control of complex traits, and the results will be applied in the livestock industry to improve meat quality and productivity, as well as to reduce the incidence of disease. A combination of quantitative trait locus mapping and microarray analysis is a useful approach to reduce the overall effort needed to identify genes associated with quantitative traits of interest.ResultsWe constructed a full-length enriched cDNA library from porcine backfat tissue. The estimated average size of the cDNA inserts was 1.7 kb, and the cDNA fullness ratio was 70%. In total, we deposited 16,110 high-quality sequences in the dbEST division of GenBank (accession numbers: DT319652-DT335761). For all the expressed sequence tags (ESTs), approximately 10.9 Mb of porcine sequence were generated with an average length of 674 bp per EST (range: 200–952 bp). Clustering and assembly of these ESTs resulted in a total of 5,008 unique sequences with 1,776 contigs (35.46%) and 3,232 singleton (65.54%) ESTs. From a total of 5,008 unique sequences, 3,154 (62.98%) were similar to other sequences, and 1,854 (37.02%) were identified as having no hit or low identity (<95%) and 60% coverage in The Institute for Genomic Research (TIGR) gene index of Sus scrofa. Gene ontology (GO) annotation of unique sequences showed that approximately 31.7, 32.3, and 30.8% were assigned molecular function, biological process, and cellular component GO terms, respectively. A total of 1,854 putative novel transcripts resulted after comparison and filtering with the TIGR SsGI; these included a large percentage of singletons (80.64%) and a small proportion of contigs (13.36%).ConclusionThe sequence data generated in this study will provide valuable information for studying expression profiles using EST-based microarrays and assist in the condensation of current pig TCs into clusters representing longer stretches of cDNA sequences. The isolation of genes expressed in backfat tissue is the first step toward a better understanding of backfat tissue on a genomic basis.

Genome-wide association study identifies major loci for carcass weight on BTA14 in Hanwoo (Korean cattle).

This genome-wide association study (GWAS) was conducted to identify major loci that are significantly associated with carcass weight, and their effects, in order to provide increased understanding of the genetic architecture of carcass weight in Hanwoo. This genome-wide association study identified one major chromosome region ranging from 23 Mb to 25 Mb on chromosome 14 as being associated with carcass weight in Hanwoo. Significant Bonferroni-corrected genome-wide associations (P<1.52×10−6) were detected for 6 Single Nucleotide Polymorphic (SNP) loci for carcass weight on chromosome 14. The most significant SNP was BTB-01280026 (P = 4.02×10−11), located in the 25 Mb region on Bos taurus autosome 14 (BTA14). The other 5 significant SNPs were Hapmap27934-BTC-065223 (P = 4.04×10−11) in 25.2 Mb, BTB-01143580 (P = 6.35×10−11) in 24.3 Mb, Hapmap30932-BTC-011225 (P = 5.92×10−10) in 24.8 Mb, Hapmap27112-BTC-063342 (P = 5.18×10−9) in 25.4 Mb, and Hapmap24414-BTC-073009 (P = 7.38×10−8) in 25.4 Mb, all on BTA 14. One SNP (BTB-01143580; P = 6.35×10−11) lies independently from the other 5 SNPs. The 5 SNPs that lie together showed a large Linkage disequilibrium (LD) block (block size of 553 kb) with LD coefficients ranging from 0.53 to 0.89 within the block. The most significant SNPs accounted for 6.73% to 10.55% of additive genetic variance, which is quite a large proportion of the total additive genetic variance. The most significant SNP (BTB-01280026; P = 4.02×10−11) had 16.96 kg of allele substitution effect, and the second most significant SNP (Hapmap27934-BTC-065223; P = 4.04×10−11) had 18.06 kg of effect on carcass weight, which correspond to 44% and 47%, respectively, of the phenotypic standard deviation for carcass weight in Hanwoo cattle. Our results demonstrated that carcass weight was affected by a major Quantitative Trait Locus (QTL) with a large effect and by many SNPs with small effects that are normally distributed.

The genome landscape of indigenous African cattle

BackgroundThe history of African indigenous cattle and their adaptation to environmental and human selection pressure is at the root of their remarkable diversity. Characterization of this diversity is an essential step towards understanding the genomic basis of productivity and adaptation to survival under African farming systems.ResultsWe analyze patterns of African cattle genetic variation by sequencing 48 genomes from five indigenous populations and comparing them to the genomes of 53 commercial taurine breeds. We find the highest genetic diversity among African zebu and sanga cattle. Our search for genomic regions under selection reveals signatures of selection for environmental adaptive traits. In particular, we identify signatures of selection including genes and/or pathways controlling anemia and feeding behavior in the trypanotolerant N’Dama, coat color and horn development in Ankole, and heat tolerance and tick resistance across African cattle especially in zebu breeds.ConclusionsOur findings unravel at the genome-wide level, the unique adaptive diversity of African cattle while emphasizing the opportunities for sustainable improvement of livestock productivity on the continent.

Use of a bovine genome array to identify new biological pathways for beef marbling in Hanwoo (Korean Cattle)

BackgroundMarbling (intramuscular fat) is a valuable trait that impacts on meat quality and an important factor determining price of beef in the Korean beef market. Animals that are destined for this high marbling market are fed a high concentrate ration for approximately 30 months in the Korean finishing farms. However, this feeding strategy leads to inefficiencies and excessive fat production. This study aimed to identify candidate genes and pathways associated with intramuscular fat deposition on highly divergent marbling phenotypes in adult Hanwoo cattle.ResultsBovine genome array analysis was conducted to detect differentially expressed genes (DEGs) in m. longissimus with divergent marbling phenotype (marbling score 2 to 7). Three data-processing methods (MAS5.0, GCRMA and RMA) were used to test for differential expression (DE). Statistical analysis identified 21 significant transcripts from at least two data-processing methods (P < 0.01). All 21 differentially expressed genes were validated by real-time PCR. Results showed a high concordance in the gene expression fold change between the microarrays and the real time PCR data. Gene Ontology (GO) and pathway analysis demonstrated that some genes (ADAMTS4, CYP51A and SQLE) over expressed in high marbled animals are involved in a protein catabolic process and a cholesterol biosynthesis process. In addition, pathway analysis also revealed that ADAMTS4 is activated by three regulators (IL-17A, TNFα and TGFβ1). QRT-PCR was used to investigate gene expression of these regulators in muscle with divergent intramuscular fat contents. The results demonstrate that ADAMTS4 and TGFβ1 are associated with increasing marbling fat. An ADAMTS4/TGFβ1 pathway seems to be associated with the phenotypic differences between high and low marbled groups.ConclusionsMarbling differences are possibly a function of complex signaling pathway interactions between muscle and fat. These results suggest that ADAMTS4, which is involved in connective tissue degradation, could play a role in an important biological pathway for building up marbling in cattle. Moreover, ADAMTS4 and TGFβ 1could potentially be used as an early biological marker for marbling fat content in the early stages of growth.

Gene expression profiling of chicken primordial germ cell ESTs

BackgroundGerm cells are the only cell type that can penetrate from one generation to next generation. At the early embryonic developmental stages, germ cells originally stem from primordial germ cells, and finally differentiate into functional gametes, sperm in male or oocyte in female, after sexual maturity. This study was conducted to investigate a large-scale expressed sequence tag (EST) analysis in chicken PGCs and compare the expression of the PGC ESTs with that of embryonic gonad.ResultsWe constructed 10,851 ESTs from a chicken cDNA library of a collection of highly separated embryonic PGCs. After chimeric and problematic sequences were filtered out using the chicken genomic sequences, there were 5,093 resulting unique sequences consisting of 156 contigs and 4,937 singlets. Pearson chi-square tests of gene ontology terms in the 2nd level between PGC and embryonic gonad set showed no significance. However, digital gene expression profiling using the Audics test showed that there were 2 genes expressed significantly with higher number of transcripts in PGCs compared with the embryonic gonads set. On the other hand, 17 genes in embryonic gonads were up-regulated higher than those in the PGC set.ConclusionOur results in this study contribute to knowledge of mining novel transcripts and genes involved in germline cell proliferation and differentiation at the early embryonic stages.

Whole-Genome Resequencing Analysis of Hanwoo and Yanbian Cattle to Identify Genome-Wide SNPs and Signatures of Selection.

Over the last 30 years, Hanwoo has been selectively bred to improve economically important traits. Hanwoo is currently the representative Korean native beef cattle breed, and it is believed that it shared an ancestor with a Chinese breed, Yanbian cattle, until the last century. However, these two breeds have experienced different selection pressures during recent decades. Here, we whole-genome sequenced 10 animals each of Hanwoo and Yanbian cattle (20 total) using the Illumina HiSeq 2000 sequencer. A total of approximately 3.12 and 3.07 billion sequence reads were mapped to the bovine reference sequence assembly (UMD 3.1) at an average of approximately 10.71- and 10.53-fold coverage for Hanwoo and Yanbian cattle, respectively. A total of 17,936,399 single nucleotide polymorphisms (SNPs) were yielded, of which 22.3% were found to be novel. By annotating the SNPs, we further retrieved numerous nonsynonymous SNPs that may be associated with traits of interest in cattle. Furthermore, we performed whole-genome screening to detect signatures of selection throughout the genome. We located several promising selective sweeps that are potentially responsible for economically important traits in cattle; the PPP1R12A gene is an example of a gene that potentially affects intramuscular fat content. These discoveries provide valuable genomic information regarding potential genomic markers that could predict traits of interest for breeding programs of these cattle breeds.

Heat shock protein B1 and its regulator genes are negatively correlated with intramuscular fat content in the longissimus thoracis muscle of Hanwoo (Korean cattle) steers.

In previous proteomic studies, heat shock protein β 1 (HSPB1) was detected as a candidate protein related to meat quality in cattle. This study sought to determine if its gene expression was associated with intramuscular fat content in the longissimus thoracis muscle of Korean cattle (Hanwoo). Tissue from two groups of 10 steers each, low-marbling (mean intramuscular fat content, 7.4 ± 1.5%) and high-marbling (23.5 ± 2.8%), were used for immunoblotting, real-time PCR, and statistical analyses. HSPB1 expression in both mRNA and protein was shown to be negatively related to intramuscular fat content (P < 0.05). Pathway analysis found two genes, TNF receptor superfamily member 6 (FAS) and angiotensinogen (AGT), that were regulators of the HSPB1 gene. The expression of the two genes showed a negative correlation with intramuscular fat content (P < 0.05). These results suggest that HSPB1, FAS, and AGT may be good candidate genes associated with intramuscular fat content in the longissimus muscle of Korean cattle.

Whole-genome resequencing analyses of five pig breeds, including Korean wild and native, and three European origin breeds

Pigs have been one of the most important sources of meat for humans, and their productivity has been substantially improved by recent strong selection. Here, we present whole-genome resequencing analyses of 55 pigs of five breeds representing Korean native pigs, wild boar and three European origin breeds. 1,673.1 Gb of sequence reads were mapped to the Swine reference assembly, covering ∼99.2% of the reference genome, at an average of ∼11.7-fold coverage. We detected 20,123,573 single-nucleotide polymorphisms (SNPs), of which 25.5% were novel. We extracted 35,458 of non-synonymous SNPs in 9,904 genes, which may contribute to traits of interest. The whole SNP sets were further used to access the population structures of the breeds, using multiple methodologies, including phylogenetic, similarity matrix, and population structure analysis. They showed clear population clusters with respect to each breed. Furthermore, we scanned the whole genomes to identify signatures of selection throughout the genome. The result revealed several promising loci that might underlie economically important traits in pigs, such as the CLDN1 and TWIST1 genes. These discoveries provide useful genomic information for further study of the discrete genetic mechanisms associated with economically important traits in pigs.