Eui-Soo Kim

Genomic adaptation of admixed dairy cattle in East Africa

Dairy cattle in East Africa imported from the U.S. and Europe have been adapted to new environments. In small local farms, cattle have generally been maintained by crossbreeding that could increase survivability under a severe environment. Eventually, genomic ancestry of a specific breed will be nearly fixed in genomic regions of local breeds or crossbreds when it is advantageous for survival or production in harsh environments. To examine this situation, 25 Friesians and 162 local cattle produced by crossbreeding of dairy breeds in Kenya were sampled and genotyped using 50K SNPs. Using principal component analysis (PCA), the admixed local cattle were found to consist of several imported breeds, including Guernsey, Norwegian Red, and Holstein. To infer the influence of parental breeds on genomic regions, local ancestry mapping was performed based on the similarity of haplotypes. As a consequence, it appears that no genomic region has been under the complete influence of a specific parental breed. Nonetheless, the ancestry of Holstein-Friesians was substantial in most genomic regions (>80%). Furthermore, we examined the frequency of the most common haplotypes from parental breeds that have changed substantially in Kenyan crossbreds during admixture. The frequency of these haplotypes from parental breeds, which were likely to be selected in temperate regions, has deviated considerably from expected frequency in 11 genomic regions. Additionally, extended haplotype homozygosity (EHH) based methods were applied to identify the regions responding to recent selection in crossbreds, called candidate regions, resulting in seven regions that appeared to be affected by Holstein-Friesians. However, some signatures of selection were less dependent on Holsteins-Friesians, suggesting evidence of adaptation in East Africa. The analysis of local ancestry is a useful approach to understand the detailed genomic structure and may reveal regions of the genome required for specialized adaptation when combined with methods for searching for the recent changes of haplotype frequency in an admixed population.

Exploring the Genetic Signature of Body Size in Yucatan Miniature Pig

Since being domesticated about 10,000–12,000 years ago, domestic pigs (Sus scrofa domesticus) have been selected for traits of economic importance, in particular large body size. However, Yucatan miniature pigs have been selected for small body size to withstand high temperature environment and for laboratory use. This renders the Yucatan miniature pig a valuable model for understanding the evolution of body size. We investigate the genetic signature for selection of body size in the Yucatan miniature pig. Phylogenetic distance of Yucatan miniature pig was compared to other large swine breeds (Yorkshire, Landrace, Duroc and wild boar). By estimating the XP-EHH statistic using re-sequencing data derived from 70 pigs, we were able to unravel the signatures of selection of body size. We found that both selections at the level of organism, and at the cellular level have occurred. Selection at the higher levels include feed intake, regulation of body weight and increase in mass while selection at the molecular level includes cell cycle and cell proliferation. Positively selected genes probed by XP-EHH may provide insight into the docile character and innate immunity as well as body size of Yucatan miniature pig.

Exploring evidence of positive selection reveals genetic basis of meat quality traits in Berkshire pigs through whole genome sequencing.

BackgroundNatural and artificial selection following domestication has led to the existence of more than a hundred pig breeds, as well as incredible variation in phenotypic traits. Berkshire pigs are regarded as having superior meat quality compared to other breeds. As the meat production industry seeks selective breeding approaches to improve profitable traits such as meat quality, information about genetic determinants of these traits is in high demand. However, most of the studies have been performed using trained sensory panel analysis without investigating the underlying genetic factors. Here we investigate the relationship between genomic composition and this phenotypic trait by scanning for signatures of positive selection in whole-genome sequencing data.ResultsWe generated genomes of 10 Berkshire pigs at a total of 100.6 coverage depth, using the Illumina Hiseq2000 platform. Along with the genomes of 11 Landrace and 13 Yorkshire pigs, we identified genomic variants of 18.9 million SNVs and 3.4 million Indels in the mapped regions. We identified several associated genes related to lipid metabolism, intramuscular fatty acid deposition, and muscle fiber type which attribute to pork quality (TG, FABP1, AKIRIN2, GLP2R, TGFBR3, JPH3, ICAM2, and ERN1) by applying between population statistical tests (XP-EHH and XP-CLR). A statistical enrichment test was also conducted to detect breed specific genetic variation. In addition, de novo short sequence read assembly strategy identified several candidate genes (SLC25A14, IGF1, PI4KA, CACNA1A) as also contributing to lipid metabolism.ConclusionsResults revealed several candidate genes involved in Berkshire meat quality; most of these genes are involved in lipid metabolism and intramuscular fat deposition. These results can provide a basis for future research on the genomic characteristics of Berkshire pigs.

Identification of signatures of selection for intramuscular fat and backfat thickness in two Duroc populations

Intramuscular fat (IMF) content is an important trait affecting the quality of pork. Two Duroc populations, one under positive selection for IMF and the other selected for decreased backfat but under stabilizing selection for IMF, were used to identify signatures of selection associated with IMF using 60,000 single-nucleotide polymorphism data. The effects of selection were analyzed between 2 lines or groups representing selected and control animals within each population using a discriminant analysis of principal components and Wright’s fixation index (FST). Moreover, extended haplotype homozygosity-based approaches were used to examine the changes in haplotype frequency due to recent selection. Each statistical method identified 10–20 selection signatures. A few haplotype-based signatures of selection agreed with results from a genome-wide association study (GWAS), while FST measures showed a better agreement with GWAS results. Agreement of marker-trait associations and signatures of selection was limited, and further examination will be necessary to understand the effect of selection on IMF and why some regions identified by GWAS did not appear to respond to the selection practiced. The genes in 21 consensus selection signatures were examined. Several genes with an effect on overall fatness were identified, but further research is needed to assess whether or not some of them could have a specific effect on IMF.

Genomic footprints of dryland stress adaptation in Egyptian fat-tail sheep and their divergence from East African and western Asia cohorts

African indigenous sheep are classified as fat-tail, thin-tail and fat-rump hair sheep. The fat-tail are well adapted to dryland environments, but little is known on their genome profiles. We analyzed patterns of genomic variation by genotyping, with the Ovine SNP50K microarray, 394 individuals from five populations of fat-tail sheep from a desert environment in Egypt. Comparative inferences with other East African and western Asia fat-tail and European sheep, reveal at least two phylogeographically distinct genepools of fat-tail sheep in Africa that differ from the European genepool, suggesting separate evolutionary and breeding history. We identified 24 candidate selection sweep regions, spanning 172 potentially novel and known genes, which are enriched with genes underpinning dryland adaptation physiology. In particular, we found selection sweeps spanning genes and/or pathways associated with metabolism; response to stress, ultraviolet radiation, oxidative stress and DNA damage repair; activation of immune response; regulation of reproduction, organ function and development, body size and morphology, skin and hair pigmentation, and keratinization. Our findings provide insights on the complexity of genome architecture regarding dryland stress adaptation in the fat-tail sheep and showcase the indigenous stocks as appropriate genotypes for adaptation planning to sustain livestock production and human livelihoods, under future climates.

Determination of Genetic Structure and Signatures of Selection in Three Strains of Tanzania Shorthorn Zebu, Boran and Friesian Cattle by Genome-Wide SNP Analyses

Background More than 90 percent of cattle in Tanzania belong to the indigenous Tanzania Short Horn Zebu (TSZ) population which has been classified into 12 strains based on historical evidence, morphological characteristics, and geographic distribution. However, specific genetic information of each TSZ population has been lacking and has caused difficulties in designing programs such as selection, crossbreeding, breed improvement or conservation. This study was designed to evaluate the genetic structure, assess genetic relationships, and to identify signatures of selection among cattle of Tanzania with the main goal of understanding genetic relationship, variation and uniqueness among them. Methodology/Principal findings The Illumina Bos indicus SNP 80K BeadChip was used to genotype genome wide SNPs in 168 DNA samples obtained from three strains of TSZ cattle namely Maasai, Tarime and Sukuma as well as two comparative breeds; Boran and Friesian. Population structure and signatures of selection were examined using principal component analysis (PCA), admixture analysis, pairwise distances (FST), integrated haplotype score (iHS), identical by state (IBS) and runs of homozygosity (ROH). There was a low level of inbreeding (F~0.01) in the TSZ population compared to the Boran and Friesian breeds. The analyses of FST, IBS and admixture identified no considerable differentiation between TSZ trains. Importantly, common ancestry in Boran and TSZ were revealed based on admixture and IBD, implying gene flow between two populations. In addition, Friesian ancestry was found in Boran. A few common significant iHS were detected, which may reflect influence of recent selection in each breed or strain. Conclusions Population admixture and selection signatures could be applied to develop conservation plan of TSZ cattle as well as future breeding programs in East African cattle.

Effect of a c-MYC Gene Polymorphism (g.3350G>C) on Meat Quality Traits in Berkshire.

c-MYC (v-myelocytomatosis viral oncogene homologue) is a transcription factor that plays important role in many biological process including cell growth and differentiation, such as myogenesis and adipogenesis. In this study, we aimed to detect MYC gene polymorphisms, their genotype frequencies and to determine associations between these polymorphisms and meat quality traits in Berkshire pigs. We identified a single nucleotide polymorphism (SNP) in intron 2 of MYC gene by Sanger sequencing, i.e., g.3350G>C (rs321898326), that is only found in Berkshire pigs, but not in other breeds including Duroc, Landrace, and Yorkshire pigs that were used in this study. Genotypes of total 378 Berkshire pigs (138 sows and 240 boars) were determined using Hha I restriction enzyme digestion after polymerase chain reaction. Observed allele frequencies of GG, GC, and CC genotypes were 0.399, 0.508, and 0.093 respectively. Statistical analysis indicated that the g.3350G>C polymorphism was significantly associated with pH45min and cooking loss (p<0.05), suggesting that g.3350G>C SNP can be used for pre-selection of pH45min and cooking loss traits in Berkshire pigs.

Characterization of the acute heat stress response in gilts: III. Genome-wide association studies of thermotolerance traits in pigs

Heat stress is one of the limiting factors negatively affecting pig production, health, and fertility. Characterizing genomic regions responsible for variation in HS tolerance would be useful in identifying important genetic factor(s) regulating physiological responses to HS. In the present study, we performed genome-wide association analyses for respiration rate (RR), rectal temperature (TR), and skin temperature (TS) during HS in 214 crossbred gilts genotyped for 68,549 single nucleotide polymorphisms (SNP) using the Porcine SNP 70K BeadChip. Considering the top 0.1% smoothed phenotypic variances explained by SNP windows, we detected 26, 26, 21, and 14 genes that reside within SNPs explaining the largest proportion of variance (top 25 SNP windows) and associated with change in RR (ΔRR) from thermoneutral (TN) conditions to HS environment, as well as the change in prepubertal TR (ΔTR), change in postpubertal ΔTR, and change in TS (ΔTS), respectively. The region between 28.85 Mb and 29.10 Mb on chromosome 16 explained about 0.05% of the observed variation for ΔRR. The growth hormone receptor (GHR) gene resides in this region and is associated with the HS response. The other important candidate genes associated with ΔRR (PAIP1, NNT, and TEAD4), ΔTR (LIMS2, TTR, and TEAD4), and ΔTS (ERBB4, FKBP1B, NFATC2, and ATP9A) have reported roles in the cellular stress response. The SNP explaining the largest proportion of variance and located within and in the vicinity of genes were related to apoptosis or cellular stress and are potential candidates that underlie the physiological response to HS in pigs.

Association of a Single Nucleotide Polymorphism (SNP) in Porcine MYC Gene with Economic Traits in Pigs

MYC (v-myelocytomatosis viral oncogene homologue) is a regulator gene that encodes for a nuclear phosphoprotein. Porcine MYC gene was mapped on chromosome SSC 4p13 and is associated with a variety of functions such as cell proliferation and cell growth. MYC expression is coupled to a multitude of physiological processes and is regulated by hormones, growth factors, cytokines, lymphokines, nutritional status, development and differentiation. MYC is also involved in myogenesis, muscle hyperplasia and adipogenesis. In this study, we investigated SNPs in MYC gene and their association with economic traits in Duroc, Landrace and Yorkshire populations. We detected a single point mutation in exon 3 of porcine MYC gene as a change of T to C at 906 base (amino acid position 302, nonsynomous mutation of alanine) in MYC-N domain. MYC mutation (T906C) was significantly associated with age at 90 kg in these breeds, signifying that this mutation can serve as a selection marker for growth traits in pigs.