Jong-Joo Kim

Resolving the evolution of extant and extinct ruminants with high-throughput phylogenomics

The Pecorans (higher ruminants) are believed to have rapidly speciated in the Mid-Eocene, resulting in five distinct extant families: Antilocapridae, Giraffidae, Moschidae, Cervidae, and Bovidae. Due to the rapid radiation, the Pecoran phylogeny has proven difficult to resolve, and 11 of the 15 possible rooted phylogenies describing ancestral relationships among the Antilocapridae, Giraffidae, Cervidae, and Bovidae have each been argued as representations of the true phylogeny. Here we demonstrate that a genome-wide single nucleotide polymorphism (SNP) genotyping platform designed for one species can be used to genotype ancient DNA from an extinct species and DNA from species diverged up to 29 million years ago and that the produced genotypes can be used to resolve the phylogeny for this rapidly radiated infraorder. We used a high-throughput assay with 54,693 SNP loci developed for Bos taurus taurus to rapidly genotype 678 individuals representing 61 Pecoran species. We produced a highly resolved phylogeny for this diverse group based upon 40,843 genome-wide SNP, which is five times as many informative characters as have previously been analyzed. We also establish a method to amplify and screen genomic information from extinct species, and place Bison priscus within the Bovidae. The quality of genotype calls and the placement of samples within a well-supported phylogeny may provide an important test for validating the fidelity and integrity of ancient samples. Finally, we constructed a phylogenomic network to accurately describe the relationships between 48 cattle breeds and facilitate inferences concerning the history of domestication and breed formation.

Polar overdominant inheritance of a DLK1 polymorphism is associated with growth and fatness in pigs

The polar overdominance model of inheritance was proposed to explain the non-Mendelian expression of callipyge muscular hypertrophy in sheep. The callipyge locus (CLPG) maps to the distal portion of ovine Chromosome 18 within the DLK1– GTL2 region and corresponds to human Chromosome 14q32, where uniparental disomy (UPD) of the region is associated with multiple congenital anomalies, including growth retardation and obesity. We investigated the porcine DLK1– GTL2 region in a cross of two pig breeds to determine if the callipyge polar overdominance is present in another species. Analyses of the parental origin of DLK1 polymorphism in the F2 offspring found that paternal inheritance of DLK1 allele 2 and maternal inheritance of the allele 1 was significantly associated with decreased fat deposition and increased lean muscle mass, while the opposite parental inheritance of these alleles was associated with slower prenatal and postnatal growth. These results suggest that the polar overdominance mode of inheritance is present in the pig chromosomal region that is homologous to the CLPG locus in sheep. Further study in pigs can provide important insights into understanding the molecular regulation of imprinted genes that are associated with human UPD14 and sheep callipyge phenotypes.

Combined line-cross and half-sib QTL analysis of crosses between outbred lines.

Data from an F 2 cross between breeds of livestock are typically analysed by least squares line-cross or half-sib models to detect quantitative trait loci (QTL) that differ between or segregate within breeds. These models can also be combined to increase power to detect QTL, while maintaining the computational efficiency of least squares. Tests between models allow QTL to be characterized into those that are fixed (LC QTL), or segregating at similar (HS QTL) or different (CB QTL) frequencies in parental breeds. To evaluate power of the combined model, data wih various differences in QTL allele frequencies (FD) between parental breeds were simulated. Use of all models increased power to detect QTL. The line-cross model was the most powerful model to detect QTL for FD>0.6. The combined and half-sib models had similar power for FD<0.4. The proportion of detected QTL declared as LC QTL decreased with FD. The opposite was observed for HS QTL. The proportion of CB QTL decreased as FD deviated from 0.5. Accuracy of map position tended to be greatest for CB QTL. Models were applied to a cross of Berkshire and Yorkshire pig breeds and revealed 160 (40) QTL at the 5% chromosome (genome)-wise level for the 39 growth, carcass composition and quality traits, of which 72, 54, and 34 were declared as LC, HS and CB QTL. Fourteen CB QTL were detected only by the combined model. Thus, the combined model can increase power to detect QTL and mapping accuracy and enable characterization of QTL that segregate within breeds.

Genetic diversity, population structure and relationships in indigenous cattle populations of Ethiopia and Korean Hanwoo breeds using SNP markers

In total, 166 individuals from five indigenous Ethiopian cattle populations – Ambo (n = 27), Borana (n = 35), Arsi (n = 30), Horro (n = 36), and Danakil (n = 38) – were genotyped for 8773 single nucleotide polymorphism (SNP) markers to assess genetic diversity, population structure, and relationships. As a representative of taurine breeds, Hanwoo cattle (n = 40) were also included in the study for reference. Among Ethiopian cattle populations, the proportion of SNPs with minor allele frequencies (MAFs) ≥0.05 ranged from 81.63% in Borana to 85.30% in Ambo, with a mean of 83.96% across all populations. The Hanwoo breed showed the highest proportion of polymorphism, with MAFs ≥0.05, accounting for 95.21% of total SNPs. The mean expected heterozygosity varied from 0.370 in Danakil to 0.410 in Hanwoo. The mean genetic differentiation (FST; 1%) in Ethiopian cattle revealed that within individual variation accounted for approximately 99% of the total genetic variation. As expected, FST and Reynold genetic distance were greatest between Hanwoo and Ethiopian cattle populations, with average values of 17.62 and 18.50, respectively. The first and second principal components explained approximately 78.33% of the total variation and supported the clustering of the populations according to their historical origins. At K = 2 and 3, a considerable source of variation among cattle is the clustering of the populations into Hanwoo (taurine) and Ethiopian cattle populations. The low estimate of genetic differentiation (FST) among Ethiopian cattle populations indicated that differentiation among these populations is low, possibly owing to a common historical origin and high gene flow. Genetic distance, phylogenic tree, principal component analysis, and population structure analyses clearly differentiated the cattle population according to their historical origins, and confirmed that Ethiopian cattle populations are genetically distinct from the Hanwoo breed.

Genome Scan for Parent-of-Origin QTL Effects on Bovine Growth and Carcass Traits

Parent-of-origin effects (POE) such as genomic imprinting influence growth and body composition in livestock, rodents, and humans. Here, we report the results of a genome scan to detect quantitative trait loci (QTL) with POE on growth and carcass traits in Angus × Brahman cattle crossbreds. We identified 24 POE–QTL on 15 Bos taurus autosomes (BTAs) of which six were significant at 5% genome-wide (GW) level and 18 at the 5% chromosome-wide (CW) significance level. Six QTL were paternally expressed while 15 were maternally expressed. Three QTL influencing post-weaning growth map to the proximal end of BTA2 (linkage region of 0–9 cM; genomic region of 5.0–10.8 Mb), for which only one imprinted ortholog is known so far in the human and mouse genomes, and therefore may potentially represent a novel imprinted region. The detected QTL individually explained 1.4 ∼ 5.1% of each trait’s phenotypic variance. Comparative in silico analysis of bovine genomic locations show that 32 out of 1,442 known mammalian imprinted genes from human and mouse homologs map to the identified QTL regions. Although several of the 32 genes have been associated with quantitative traits in cattle, only two (GNAS and PEG3) have experimental proof of being imprinted in cattle. These results lend additional support to recent reports that POE on quantitative traits in mammals may be more common than previously thought, and strengthen the need to identify and experimentally validate cattle orthologs of imprinted genes so as to investigate their effects on quantitative traits.

Characterizations of Shell and Mantle Edge Pigmentation of a Pacific Oyster, Crassostrea gigas, in Korean Peninsula

The objectives of this study were to investigate color patterns of shell and mantle edge pigmentation of a Pacific oyster, C. gigas, and to estimate variance components of the two colors. A sample of 240 F0 oysters was collected from six aquaculture farms in Tongyeong, Korea to measure shell color and mantle edge pigmentation. Among the F0s, male and female individuals with black (white) shell and black (white) mantle edge were selected and mated to generate three F1 full-sib black (white) cross families (N = 265). Two and four F2 cross families (N = 286) were also produced from black and white F1 selected individuals, respectively. Variance component estimates due to residuals and families within color were obtained using SAS PROC VARCOMP procedures to estimate heritability of shell and mantle edge pigmentation. In the F0 generation, about 29% (11%) had black (white) color for both shell and mantle edge. However, in the F1 and F2 black (white) cross families, 75% (67%) and 100% (100%) of oysters had black (white) shell colors, and 59% (23%) and 79% (55%) had black (white) mantle edge, respectively. Spearman correlation coefficients between shell and mantle edge color were 0.25, 0.74, and 0.92 in F0, F1, and F2 generations, respectively, indicating that, with generations of selection process, an individual with black (white) shell color is more likely to have black (white) mantle edge pigmentation. This suggests that shell color could be a good indicator trait for mantle edge pigmentation if selection of both the colors is implemented for a couple of generations. Estimates of heritability were 0.41 and 0.77 for shell color and 0.27 and 0.08 for mantle edge pigmentation in the F1 and F2 generations, respectively, indicating that, in general, significant proportions of phenotypic variations for the shell and mantle edge colors are explained by genetic variations between individuals. These results suggest that the two color traits are inheritable and correlated, enabling effective selection on shell and mantle edge color.

Genome analysis of the domestic dog (Korean Jindo) by massively parallel sequencing.

Although pioneering sequencing projects have shed light on the boxer and poodle genomes, a number of challenges need to be met before the sequencing and annotation of the dog genome can be considered complete. Here, we present the DNA sequence of the Jindo dog genome, sequenced to 45-fold average coverage using Illumina massively parallel sequencing technology. A comparison of the sequence to the reference boxer genome led to the identification of 4 675 437 single nucleotide polymorphisms (SNPs, including 3 346 058 novel SNPs), 71 642 indels and 8131 structural variations. Of these, 339 non-synonymous SNPs and 3 indels are located within coding sequences (CDS). In particular, 3 non-synonymous SNPs and a 26-bp deletion occur in the TCOF1 locus, implying that the difference observed in cranial facial morphology between Jindo and boxer dogs might be influenced by those variations. Through the annotation of the Jindo olfactory receptor gene family, we found 2 unique olfactory receptor genes and 236 olfactory receptor genes harbouring non-synonymous homozygous SNPs that are likely to affect smelling capability. In addition, we determined the DNA sequence of the Jindo dog mitochondrial genome and identified Jindo dog-specific mtDNA genotypes. This Jindo genome data upgrade our understanding of dog genomic architecture and will be a very valuable resource for investigating not only dog genetics and genomics but also human and dog disease genetics and comparative genomics.

Effect of Population Reduction on mtDNA Diversity and Demographic History of Korean Cattle Populations.

The population sizes of three Korean indigenous cattle populations have been drastically reduced over the past decades. In this study, we examined the extent to which reduction in populations influenced genetic diversity, population structure and demographic history using complete mitochondrial DNA (mtDNA) control region sequences. The complete mtDNA control region was sequenced in 56 individuals from Korean Black (KB), Jeju Black (JEB) and Korean Brindle (BRI) cattle populations. We included 27 mtDNA sequences of Korean Brown (BRO) from the GenBank database. Haplotype diversity estimate for the total population was high (0.870) while nucleotide diversity was low (0.004). The KB showed considerably low nucleotide (π = 0.001) and haplotype (h = 0.368) diversities. Analysis of molecular variance revealed a low level of genetic differentiation but this was highly significant (p<0.001) among the cattle populations. Of the total genetic diversity, 7.6% was attributable to among cattle populations diversity and the rest (92.4%) to differences within populations. The mismatch distribution analysis and neutrality tests revealed that KB population was in genetic equilibrium or decline. Indeed, unless an appropriate breeding management practice is developed, inbreeding and genetic drift will further impoverish genetic diversity of these cattle populations. Rational breed development and conservation strategy is needed to safeguard these cattle population.

Evaluation of Single Nucleotide Polymorphisms (SNPs) Genotyped by the Illumina Bovine SNP50K in Cattle Focusing on Hanwoo Breed

In the present study, we evaluated the informativeness of SNPs genotyped by the Illumina Bovine SNP50K assay in different cattle breeds. To investigate these on a genome-wide scale, we considered 52,678 SNPs spanning the whole autosomal and X chromosomes in cattle. Our study samples consists of six different cattle breeds. Across the breeds approximately 72 and 6% SNPs were found polymorphic and fixed or close to fix in all the breeds, respectively. The variations in the average minor allele frequency (MAF) were significantly different between the breeds studied. The level of average MAF observed in Hanwoo was significantly lower than the other breeds. Hanwoo breed also displayed the lowest number of polymorphic SNPs across all the chromosomes. More importantly, this study indicated that the Bovine SNP50K assay will have reduced power for genome-wide association studies in Hanwoo as compared to other cattle breeds. Overall, the Bovine SNP50K assay described in this study offer a useful genotyping platform for mapping quantitative trait loci (QTLs) in the cattle breeds. The assay data represent a vast and generally untapped resource to assist the investigation of the complex production traits and the development of marker-assisted selection programs.

A Genome Wide Association Study on Age at First Calving Using High Density Single Nucleotide Polymorphism Chips in Hanwoo (Bos taurus coreanae)

Age at first calving is an important trait for achieving earlier reproductive performance. To detect quantitative trait loci (QTL) for reproductive traits, a genome wide association study was conducted on the 96 Hanwoo cows that were born between 2008 and 2010 from 13 sires in a local farm (Juk-Am Hanwoo farm, Suncheon, Korea) and genotyped with the Illumina 50K bovine single nucleotide polymorphism (SNP) chips. Phenotypes were regressed on additive and dominance effects for each SNP using a simple linear regression model after the effects of birth-year-month and polygenes were considered. A forward regression procedure was applied to determine the best set of SNPs for age at first calving. A total of 15 QTL were detected at the comparison-wise 0.001 level. Two QTL with strong statistical evidence were found at 128.9 Mb and 111.1 Mb on bovine chromosomes (BTA) 2 and 7, respectively, each of which accounted for 22% of the phenotypic variance. Also, five significant SNPs were detected on BTAs 10, 16, 20, 26, and 29. Multiple QTL were found on BTAs 1, 2, 7, and 14. The significant QTLs may be applied via marker assisted selection to increase rate of genetic gain for the trait, after validation tests in other Hanwoo cow populations.