Identifying genetic loci controlling neonatal passive transfer of immunity using a hybrid genotyping strategy

Colostrum intake is critical to a piglet's survival and can be measured by precipitating out the γ‐immunoglobulins from serum with ammonium sulfate (immunocrit). Genetic analysis of immunocrits on 5312 piglets indicated that the heritabilities (se) for direct and maternal effects were 0.13 (0.06) and 0.53 (0.08) respectively. To identify QTL for direct genetic effects, piglets with the highest and lowest immunocrits from 470 litters were selected. Six sets of DNA pools were created based on sire of the litter. These 12 DNA pools were applied to Illumina Porcine SNP60 BeadChips. Normalized X and Y values were analyzed. Three different SNP selection methods were used: deviation of the mean from high vs. low pools, the deviation adjusted for variance based on binomial theory and ANOVA. The 25 highest ranking SNPs were selected from each evaluation for further study along with 12 regions selected based on a five‐SNP window approach. Selected SNPs were individually genotyped in the 988 piglets included in pools as well as in 524 piglets that had intermediate immunocrits. Association analyses were conducted fitting an animal model using the estimated genetic parameters. Nineteen SNPs were nominally associated (P < 0.01) with immunocrit values, of which nine remained significant (P < 0.05) after Bonferroni correction, located in 16 genomic regions on 13 chromosomes. In conclusion, the pooling strategy reduced the cost to scan the genome by more than 80% and identified genomic regions associated with a piglet's ability to acquire γ‐immunoglobulin from colostrum. Each method to rank SNPs from the pooled analyses contributed unique validated markers, suggesting that multiple analyses will reveal more QTL than a single analysis.     

Effect of swine leukocyte antigen‐DQA gene variation on diarrhea in Large White,Landrace, and Duroc piglets

Piglet diarrhea is one of the most common factors that affects the benefits of the swine industry. Although recent studies have shown that exon 2 of SLA‐DQA is associated with piglet resistance to diarrhea, contributions of genetic variation in the additional exon coding regions of this gene remain unclear. Here, we investigated variation in exons 1, 3 and 4 of the SLA‐DQA gene and evaluated their effects on diarrheal infection in 425 suckling piglets. No variation was identified in exon 1. In exon 3, there were eight alleles detected, generated by 14 single nucleotide polymorphisms (SNPs) and three nucleotide deletions, eight SNPs being newly identified. Four allele sequences and three SNPs were identified in exon 4, only one SNP being newly identified. Statistical analysis showed that the genotypes of exon 3 are significantly associated with piglet diarrhea; indeed, genotypes DQA*wb01/wb02 and wb04/wb05 are clearly associated with resistance to piglet diarrhea, as they have the lowest probabilities of infection (P < 0.05). However, no significant association was found between the genotypes of exon 4 and diarrhea (P > 0.05). These results provide important new information concerning the level of genetic diversity at the SLA‐DQA locus and suggest that further genetic association studies of piglet diarrhea resistance should include analyses of both exons 2 and 3 of this locus.     

Genetic markers that influence feed efficiency phenotypes also affect cattle temperament as measured by flight speed

Flight speed is a predictive indicator of cattle temperament and is associated with feed efficiency phenotypes. Genetic markers associated with both traits may assist with selection of calmer animals with improved economic value. A preliminary genome‐wide association study determined chromosomal regions on BTA9, and 17 were associated with flight speed. The genes quaking (QKI), glutamate receptor, ionotropic, AMPA 2 (GRIA2) and glycine receptor β (GLRB) were identified in these regions as potential functional candidates. Beef steers (n = 1057) were genotyped with SNPs located within and flanking these genes. One SNP located near QKI and one near GRIA2 were nominally associated with flight speed (P ≤ 0.05) although neither was significant after Bonferroni correction. Several studies have shown a correlation between flight speed and feed intake or gain; therefore, we also analyzed SNPs on BTA6:38–39 Mb known to be associated with average daily gain (ADG) and average daily feed intake (ADFI) for association with flight speed. Several SNPs on BTA6 were associated with flight speed (P ≤ 0.005), and three were significant after Bonferroni correction. These results suggest that the genes tested are unlikely to contribute to flight speed variation for our cattle population, but SNPs on BTA6 associated with ADG and ADFI may influence temperament. Use of these markers to select for economically important feed efficiency phenotypes may produce cattle with more desirable temperaments.     

Identification of candidate genes and mutations in QTL regions for immune responses in chicken

There are two categories of immune responses – innate and adaptive immunity – both having polygenic backgrounds and a significant environmental component. In our study, adaptive immunity was represented by the specific antibody response toward keyhole limpet hemocyanin (KLH); innate immunity was represented by natural antibodies toward lipopolysaccharide (LPS) and lipoteichoic acid (LTA). Defining genetic bases of immune responses leads from defining quantitative trait loci (QTL) toward a single mutation responsible for variation in the phenotypic trait. The goal of the reported study was to define candidate genes and mutations for the immune traits of interest in chicken by performing an association study of SNPs located in candidate genes defined in QTL regions. Candidate genes and SNPs in QTL regions were selected in silico. SNP association was based on a custom SNP panel, GoldenGate genotyping assay (Illumina) and two statistical models: random mixed model and CAR score. The most significant SNP for immune response toward KLH was located in the JMJD6 gene located on GGA18. Four SNPs in candidate genes FOXJ1 (GGA18), EPHB1 (GGA9), PTGER4 (GGAZ) and PRKCB (GGA14) showed association with natural antibodies for LPS. A single SNP in ITGB4 (GGA18) was associated with natural antibodies for LTA. All associated SNPs mentioned above showed additive effects.     

New resources for genetic studies in Populus nigra: genome‐wide SNP discovery and development of a 12k Infinium array

Whole genome resequencing of 51 Populus nigra (L.) individuals from across Western Europe was performed using Illumina platforms. A total number of 1 878 727 SNPs distributed along the P. nigra reference sequence were identified. The SNP calling accuracy was validated with Sanger sequencing. SNPs were selected within 14 previously identified QTL regions, 2916 expressional candidate genes related to rust resistance, wood properties, water‐use efficiency and bud phenology and 1732 genes randomly spread across the genome. Over 10 000 SNPs were selected for the construction of a 12k Infinium Bead‐Chip array dedicated to association mapping. The SNP genotyping assay was performed with 888 P. nigra individuals. The genotyping success rate was 91%. Our high success rate was due to the discovery panel design and the stringent parameters applied for SNP calling and selection. In the same set of P. nigra genotypes, linkage disequilibrium throughout the genome decayed on average within 5–7 kb to half of its maximum value. As an application test, ADMIXTURE analysis was performed with a selection of 600 SNPs spread throughout the genome and 706 individuals collected along 12 river basins. The admixture pattern was consistent with genetic diversity revealed by neutral markers and the geographical distribution of the populations. These newly developed SNP resources and genotyping array provide a valuable tool for population genetic studies and identification of QTLs through natural‐population based genetic association studies in P. nigra.     

A genome‐wide association study of heat stress‐associated SNPs in catfish

Heat tolerance is a complex and economically important trait for catfish genetic breeding programs. With global climate change, it is becoming an increasingly important trait. To better understand the molecular basis of heat stress, a genome‐wide association study (GWAS) was carried out using the 250 K catfish SNP array with interspecific backcross progenies, which derived from crossing female channel catfish with male F1 hybrid catfish (female channel catfish × male blue catfish). Three significant associated SNPs were detected by performing an EMMAX approach for GWAS. The SNP located on linkage group 14 explained 12.1% of phenotypical variation. The other two SNPs, located on linkage group 16, explained 11.3 and 11.5% of phenotypical variation respectively. A total of 14 genes with heat stress related functions were detected within the significant associated regions. Among them, five genes—TRAF2, FBXW5, ANAPC2, UBR1 and KLHL29— have known functions in the protein degradation process through the ubiquitination pathway. Other genes related to heat stress include genes involved in protein biosynthesis (PRPF4 and SYNCRIP), protein folding (DNAJC25), molecule and iron transport (SLC25A46 and CLIC5), cytoskeletal reorganization (COL12A1) and energy metabolism (COX7A2, PLCB1 and PLCB4) processes. The results provide fundamental information about genes and pathways that is useful for further investigation into the molecular mechanisms of heat stress. The associated SNPs could be promising candidates for selecting heat‐tolerant catfish lines after validating their effects on larger and various catfish populations.     

Identifying chromosomal selection‐sweep regions in facial eczema selection‐line animals using an ovine 50K‐SNP array

Facial eczema (FE) is a hepato‐mycotoxicosis found mainly in New Zealand sheep and cattle. When genetics was found to be a factor in FE susceptibility, resistant and susceptible selection lines of Romney sheep were established to enable further investigations of this disease trait. Using the Illumina OvineSNP50 BeadChip, we conducted a selection‐sweep experiment on these FE genetic lines. Two analytical methods were used to detect selection signals, namely the Peddrift test (Dodds & McEwan, 1997) and fixation index F_ST (Weir & Hill, 2002). Of 50 975 single nucleotide polymorphism (SNP) markers tested, there were three that showed highly significant allele frequency differences between the resistant and susceptible animals (Peddrift nominal P < 0.000001). These SNP loci are located on chromosomes OAR1, OAR11 and OAR12 that coincide precisely with the three highest genomic F_ST peaks. In addition, there are nine less significant Peddrift SNPs (nominal P ≤ 0.000009) on OAR6 (n = 2), OAR9 (n = 2), OAR12, OAR19 (n = 2), OAR24 and OAR26. In smoothed F_ST (five‐SNP moving average) plots, the five most prominent peaks are on OAR1, OAR6, OAR7, OAR13 and OAR19. Although these smoothed F_ST peaks do not coincide with the three most significant Peddrift SNP loci, two (on OAR6 and OAR19) overlap with the set of less significant Peddrift SNPs above. Of these 12 Peddrift SNPs and five smoothed F_ST regions, none is close to the FE candidate genes catalase and ABCG2; however, two on OAR1 and one on OAR13 fall within suggestive quantitative trait locus regions identified in a previous genome screen experiment. The present studies indicated that there are at least eight genomic regions that underwent a selection sweep in the FE lines.     

Detection of quantitative trait loci for resistance to gastrointestinal nematode infections in Creole goats

This study aimed to identify regions of the genome affecting resistance to gastrointestinal nematodes in a Creole goat population naturally exposed to a mixed nematode infection (Haemonchus contortus, Trichostrongylus colubriformis and Oesophagostomum columbianum) by grazing on irrigated pasture. A genome‐wide quantitative trait loci (QTL) scan was performed on 383 offspring from 12 half‐sib families. A total of 101 microsatellite markers were genotyped. Traits analysed were faecal egg count (FEC), packed cell volume (PCV), eosinophil count and bodyweight (BW) at 7 and 11 months of age. Levels of activity of immunoglobulin A (IgA) and activity of immunoglobulin E (IgE) anti‐Haemonchus contortus L3 crude extracts and adult excretion/secretion products (ESPs) were also analysed. Using interval mapping, this study identified 13 QTL for parasite resistance. Two QTL linked with FEC were found on chromosomes 22 and 26. Three QTL were detected on chromosomes 7, 8 and 14 for eosinophil counts. Three QTL linked with PCV were identified on chromosomes 5, 9 and 21. A QTL for BW at 7 months of age was found on chromosome 6. Lastly, two QTL detected on chromosomes 3 and 10 were associated with IgE anti‐L3, and IgE anti‐ESP was linked with two QTL on chromosomes 1 and 26. This study is the first to have identified regions of the genome linked with nematode resistance in a goat population using a genome scan. These results provide useful tools for the understanding of parasite resistance in small ruminants.     

A candidate gene association study for nine economically important traits in Italian Holstein cattle

We genotyped 58 single nucleotide polymorphisms (SNPs) in 25 candidate genes in about 800 Italian Holstein sires. Fifty‐six (minor allele frequency >0.02) were used to evaluate their association with single traits: milk yield (MY), milk fat yield (FY), milk protein yield (PY), milk fat percentage (FP), milk protein percentage (PP), milk somatic cell count (MSCC); and complex indexes: longevity, fertility and productivity–functionality type (PFT), using deregressed proofs, after adjustment for familial relatedness. Thirty‐two SNPs were significantly associated (proportion of false positives <0.05) with different traits: 16 with MSCC, 15 with PY, 14 with MY, 12 with PFT, eight with longevity, eight with FY, eight with PP, five with FP and two with fertility. In particular, a SNP in the promoter region of the PRLR gene was associated with eight of nine traits. DGAT1 polymorphisms were highly associated with FP and FY. Casein gene markers were associated with several traits, confirming the role of the casein gene cluster in affecting milk yield, milk quality and health traits. Other SNPs in genes located on chromosome 6 were associated with PY, PP, PFT, MY (PPARGC1A) and MSCC (KIT). This latter association may suggest a biological link between the degree of piebaldism in Holstein and immunological functions affecting somatic cell count and mastitis resistance. Other significant SNPs were in the ACACA, CRH, CXCR1, FASN, GH1, LEP, LGB (also known as PAEP), MFGE8, SRC, TG, THRSP and TPH1 genes. These results provide information that can complement QTL mapping and genome‐wide association studies in Holstein.     

Signatures of selection in sheep bred for resistance or susceptibility to gastrointestinal nematodes

Background

Gastrointestinal nematodes are one of the most serious causes of disease in domestic ruminants worldwide. There is considerable variation in resistance to gastrointestinal nematodes within and between sheep breeds, which appears to be due to underlying genetic diversity. Through selection of resistant animals, rapid genetic progress has been demonstrated in both research and commercial flocks. Recent advances in genome sequencing and genomic technologies provide new opportunities to understand the ovine host response to gastrointestinal nematodes at the molecular level, and to identify polymorphisms conferring nematode resistance.

Results

Divergent lines of Romney and Perendale sheep, selectively bred for high and low faecal nematode egg count, were genotyped using the Illumina® Ovine SNP50 BeadChip. The resulting genome-wide SNP data were analysed for selective sweeps on loci associated with resistance or susceptibility to gastrointestinal nematode infection. Population differentiation using F_ST and Peddrift revealed sixteen regions, which included candidate genes involved in chitinase activity and the cytokine response. Two of the sixteen regions identified were contained within previously identified QTLs associated with nematode resistance.

Conclusions

In this study we identified fourteen novel regions associated with resistance or susceptibility to gastrointestinal nematodes. Results from this study support the hypothesis that host resistance to internal nematode parasites is likely to be controlled by a number of loci of moderate to small effects.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-637) contains supplementary material, which is available to authorized users.     

Genome-wide scan reveals important additive and non-additive genetic effects associated with resistance to Haemonchus contortus in Florida Native sheep

Florida Native sheep is among the sheep breeds best adapted to humid and hot climatic conditions such as those of Florida, USA, and have shown a superior ability to regulate nematode burdens. This is one of the oldest sheep breeds in North America and is an endangered species. To ensure genetic diversity and long-term survival of the breed, protection of the current genetic stock is critical and conservation efforts are required to promote its breeding and production. The objective of the present study was to investigate the importance of additive and non-additive genetic effects on resistance to natural Haemonchus contortus infections in Florida Native sheep using a whole genome scan. A total of 200 sheep were evaluated in the present study. Phenotypic records included faecal egg count (FEC, eggs/gram), FAMACHA® score, packed cell volume (PCV, %), body condition score and average daily gain (ADG, kg). Sheep were genotyped using the GGP Ovine 50K SNP chip and 45.2 k single nucleotide polymorphism (SNP) markers spanning the entire genome were available for quality control procedures. Mixed models were used to analyse the response variables and included the identity by state matrix to control for population structure. Bonferroni correction was used to control for multiple testing and a second arbitrary threshold (0.1 × 10^?3) was used. Fifteen SNPs with additive and non-additive genetic effects and located in Ovis aries chromosome OAR1, 2, 3, 6, 8, 10, 11, 12, 13 and 21 were associated with FEC, FAMACHA® score, PCV and ADG. These SNPs could be potential genetic markers for resistance to natural H. contortus exposure in Florida Native sheep.     

Genome‐wide detection of genetic loci and candidate genes for teat number and body conformation traits at birth in Chinese Sushan pigs

Body conformation at birth and teat number are economically important traits in the pig industry, as these traits are usually explored to evaluate the growth and reproductive potential of piglets. To detect genetic loci and candidate genes for these traits, we performed a GWAS on 269 pigs from a recently developed Chinese breed (Sushan) using 38  128 informative SNPs on the Affymetrix Porcine SNP 55K Array. In total, we detected one genome‐wide significant (P = 1.31e‐6) SNP for teat number on chromosome X and 15 chromosome‐wide significant SNPs for teat number, body weight, body length, chest circumference and cannon circumference at birth on chromosomes 1, 3, 4, 6, 7, 9, 10, 13, 14, 15, 17 and 18. The most significant SNP had an additive effect of 0.74 × total teat number, explaining 20% of phenotypic variance. Five significant SNPs resided in the previously reported quantitative trait loci for these traits and seven significant SNPs had a pleiotropic effect on multiple traits. Intriguingly, 12 of the genes nearest to the significant SNPs are functionally related to body conformation and teat number traits, including SPRED2, MKX, TMSB4X and ESR1. GO analysis revealed that candidate genes proximal to the significant SNPs were enriched in the G‐protein coupled receptor and steroid hormone‐mediated signaling pathway. Our findings shed light on the genetic basis of the measured traits and provide molecular markers especially for the genetic improvement of teat number in Sushan and related pigs.     

Genetic variation at KIT locus may predispose to melanoma

As loss of KIT frequently occurs in melanoma progression, we hypothesized that KIT is implicated in predisposition to melanoma (MM). Thus, we sequenced the KIT coding region in 112 familial MM cases and 143 matched controls and genotyped tag single‐nucleotide polymorphisms (SNPs) in two cohorts of melanoma patients and matched controls. Five rare KIT substitutions, all predicted possibly or probably deleterious, were identified in five patients, but none in controls [RR = 2.26 (1.26–2.26)]. Expressed in melanocyte lines, three substitutions inhibited KIT signaling. Comparison with exomes database (7020 alleles) confirmed a significant excess of rare deleterious KIT substitutions in patients. Additionally, a common SNP, rs2237028, was associated with MM risk, and 6 KIT variants were associated with nevus count. Our data strongly suggest that rare KIT substitutions predispose to melanoma and that common variants at KIT locus may also impact nevus count and melanoma risk.     

Identification of new single nucleotide polymorphisms affecting total number born and candidate genes related to ovulation rate in Chinese Erhualian pigs

The Chinese Erhualian pig has the highest record for litter size in the world. However, the genetic mechanism of its high prolificacy remains poorly understood. In our study, large phenotypic variations in litter size were found among Erhualian sows. Significant differences in total number born (TNB) and corpora lutea numbers were observed between sows with high and low estimated breeding values (EBVs) for TNB. To identify single nucleotide polymorphisms (SNPs) associated with TNB, a selective genomic scan was conducted on 18 sows representing the top 10% and 18 sows representing the bottom 10% of EBVs of 177 sows using Illumina Porcine SNP60 genotype data. Genome‐wide fixation coefficient (F_ST) values were calculated for each SNP between the high‐ and low‐EBV groups. A total of 154 SNPs were significantly differentiated loci between the two groups. Of the top 10 highest F_ST SNPs, rs81399474, rs81400131 and rs81405013 on SSC8 and rs81434499 and rs81434489 on SSC 12 corresponded to previously reported QTL for litter size. The other five SNPs, rs81367039 on SSC2, rs80891106 on SSC7, rs81477883 on SSC12 and rs80938898 and rs80971725 on SSC14, appeared to be novel QTL for TNB. Significant associations between rs81399474 on SSC8 and TNB were confirmed in 313 Erhualian sows. Forty genes were identified around the top 10 highest F_ST SNPs, of which UCHL1, adjacent to rs81399474, and RPS6KB1 and CLTC, adjacent to rs81434499, have been reported to affect the ovulation rate in pig. The findings can advance understanding of the genetic variations in litter size of pigs.     

Genome‐wide association study identifies an NLR gene that confers partial resistance to Magnaporthe oryzae in rice

Because of the frequent breakdown of major resistance (R) genes, identification of new partial R genes against rice blast disease is an important goal of rice breeding. In this study, we used a core collection of the Rice Diversity Panel II (C‐RDP‐II), which contains 584 rice accessions and are genotyped with 700 000 single‐nucleotide polymorphism (SNP) markers. The C‐RDP‐II accessions were inoculated with three blast strains collected from different rice‐growing regions in China. Genome‐wide association study identified 27 loci associated with rice blast resistance (LABRs). Among them, 22 LABRs were not associated with any known blast R genes or QTLs. Interestingly, a nucleotide‐binding site leucine‐rich repeat (NLR) gene cluster exists in the LABR12 region on chromosome 4. One of the NLR genes is highly conserved in multiple partially resistant rice cultivars, and its expression is significantly up‐regulated at the early stages of rice blast infection. Knockout of this gene via CRISPR‐Cas9 in transgenic plants partially reduced blast resistance to four blast strains. The identification of this new non‐strain specific partial R gene, tentatively named rice blast Partial Resistance gene 1 (PiPR1), provides genetic material that will be useful for understanding the partial resistance mechanism and for breeding durably resistant cultivars against blast disease of rice.     

Genome‐wide association reveals the locus responsible for four‐horned ruminant

Phenotypic variability in horn characteristics, such as their size, number and shape, offers the opportunity to elucidate the molecular basis of horn development. The objective of this study was to map the genetic determinant controlling the production of four horns in two breeds, Jacob sheep and Navajo‐Churro, and examine whether an eyelid abnormality occurring in the same populations is related. Genome‐wide association mapping was performed using 125 animals from the two breeds that contain two‐ and four‐horned individuals. A case–control design analysis of 570 712 SNPs genotyped with the ovine HD SNP Beadchip revealed a strong association signal on sheep chromosome 2. The 10 most strongly associated SNPs were all located in a region spanning Mb positions 131.9–132.6, indicating the genetic architecture underpinning the production of four horns is likely to involve a single gene. The closest genes to the most strongly associated marker (OAR2_132568092) were MTX2 and the HOXD cluster, located approximately 93 Kb and 251 Kb upstream respectively. The occurrence of an eyelid malformation across both breeds was restricted to polled animals and those carrying more than two horns. This suggests the eyelid abnormality may be associated with departures from the normal developmental production of two‐horned animals and that the two conditions are developmentally linked. This study demonstrated the presence of separate loci responsible for the polled and four‐horned phenotypes in sheep and advanced our understanding of the complexity that underpins horn morphology in ruminants.     

Identification of an Interaction between VWF rs7965413 and Platelet Count as a Novel Risk Marker for Metabolic Syndrome: An Extensive Search of Candidate Polymorphisms in a Case-Control Study

Although many single nucleotide polymorphisms (SNPs) have been identified to be associated with metabolic syndrome (MetS), there was only a slight improvement in the ability to predict future MetS by the simply addition of SNPs to clinical risk markers. To improve the ability to predict future MetS, combinational effects, such as SNP—SNP interaction, SNP—environment interaction, and SNP—clinical parameter (SNP × CP) interaction should be also considered. We performed a case-control study to explore novel SNP × CP interactions as risk markers for MetS based on health check-up data of Japanese male employees. We selected 99 SNPs that were previously reported to be associated with MetS and components of MetS; subsequently, we genotyped these SNPs from 360 cases and 1983 control subjects. First, we performed logistic regression analyses to assess the association of each SNP with MetS. Of these SNPs, five SNPs were significantly associated with MetS (P < 0.05): LRP2 rs2544390, rs1800592 between UCP1 and TBC1D9, APOA5 rs662799, VWF rs7965413, and rs1411766 between MYO16 and IRS2. Furthermore, we performed multiple logistic regression analyses, including an SNP term, a CP term, and an SNP × CP interaction term for each CP and SNP that was significantly associated with MetS. We identified a novel SNP × CP interaction between rs7965413 and platelet count that was significantly associated with MetS [SNP term: odds ratio (OR) = 0.78, P = 0.004; SNP × CP interaction term: OR = 1.33, P = 0.001]. This association of the SNP × CP interaction with MetS remained nominally significant in multiple logistic regression analysis after adjustment for either the number of MetS components or MetS components excluding obesity. Our results reveal new insight into platelet count as a risk marker for MetS.     

Sequence analysis of three canine adipokine genes revealed an association between TNF polymorphisms and obesity in Labrador dogs

Obesity is an emerging health problem in purebred dogs. Due to their crucial role in energy homeostasis control, genes encoding adipokines are considered candidate genes, and their variants may be associated with predisposition to obesity. Searching for polymorphism was carried out in three adipokine genes (TNF, RETN and IL6). The study was performed on 260 dogs, including lean (n = 109), overweight (n = 88) and obese (n = 63) dogs. The largest cohort was represented by Labrador Retrievers (n = 136). Altogether, 24 novel polymorphisms were identified: 12 in TNF (including one missense SNP), eight in RETN (including one missense SNP) and four in IL6. Distributions of five common SNPs (two in TNF, two in RETN and one in IL6) were further analyzed with regard to body condition score. Two SNPs in the non‐coding parts of TNF (c.‐40A>C and c.233+14G>A) were associated with obesity in Labrador dogs. The obtained results showed that the studied adipokine genes are highly polymorphic and two polymorphisms in the TNF gene may be considered as markers predisposing Labrador dogs to obesity.     

Validation of genetic polymorphisms on BTA14 associated with carcass trait in a commercial Hanwoo population

The objective of this study was to validate the association of significant SNPs identified from a previous genome‐wide association study with carcass weight (CWT) in a commercial Hanwoo population. We genotyped 13 SNPs located on BTA14 in 867 steers from Korea Hanwoo feedlot bulls. Of these 13 SNPs, five SNPs, namely rs29021868, rs110061498, rs109546980, rs42404006 and rs42303720, were found to be significantly associated (P < 0.001) with CWT. These five significant markers spanned the 24.3 to 29.4 Mb region of BTA14. The most significant marker (rs29021868) for CWT in this study had a 13.07 kg allele substitution effect and accounted for 2.4% of the additive genetic variance in the commercial Hanwoo population. The SNP marker rs109546980 was found to be significantly associated with both CWT (P < 0.001) and eye muscle area (P < 0.001) and could potentially be exploited for marker‐assisted selection in Hanwoo cattle. We also genotyped the ss319607402 variation, which maps to intron2 of PLAG1 gene and which is already reported to be associated with height, to identify any significant association with carcass weight; however, no such association was observed in this Hanwoo commercial population.