Association of TLR4 polymorphisms with Mycobacterium avium subspecies paratuberculosis infection status in Canadian Holsteins

Mycobacterium avium ssp. paratuberculosis (MAP) causes chronic enteritis in cattle that results in substantial financial losses to the cattle industry worldwide. Given that susceptibility to MAP infection is determined in part by genetics, marker‐assisted selection may help in the breeding of animals that are more resistant to MAP infection. The toll‐like receptor 4 gene (TLR4) was selected as a potential candidate gene because of its role in innate immunity and its involvement in MAP recognition and infection. The objective of this study, therefore, was to identify associations between TLR4 polymorphisms and susceptibility to MAP infection in Canadian Holstein cows. Two biologically relevant SNPs, including c.‐226G>C in the 5′‐untranslated region and the non‐synonymous SNP c.2021C>T in the potential TIR domain, were selected for an association analysis with MAP infection status in 409 Canadian Holsteins. The haplotype C‐T from these combined SNPs yielded significant association with susceptibility to MAP infection, supporting the involvement of TLR4 in susceptibility to MAP infection.     

A cohort study of racing performance in Japanese Thoroughbred racehorses using genome information on ECA18

Using 1710 Thoroughbred racehorses in Japan, a cohort study was performed to evaluate the influence of genotypes at four single nucleotide polymorphisms (SNPs) on equine chromosome 18 (ECA18), which were associated in a previous genome‐wide association study for racing performance with lifetime earnings and performance rank. In males, both g.65809482T>C and g.65868604G>T were related to performance rank (P = 0.005). In females, g.65809482T>C (P = 1.76E‐6), g.65868604G>T (P = 6.81E‐6) and g.66493737C>T (P = 4.42E‐5) were strongly related to performance rank and also to lifetime earnings (P < 0.05). When win‐race distance (WRD) among all winning racehorses and best race distance (BRD) among elite racehorses were considered as the phenotypes, significant associations (P < 0.001) were observed for all four SNPs. The favourable race distance of both elite (BRD) and novice racehorses (WRD) was also associated with genotypes in the ECA18 region, indicating the presence of a gene in this region influencing optimum race distance in Thoroughbred racehorses. Therefore, the association with performance rank is likely due to the bias in the race distances. The location of the SNPs within and proximal to the gene encoding myostatin (MSTN) strongly suggests that regulation of the MSTN gene affects racing performance. In particular, the g.65809482T>C, g.65868604G>T and g.66493737C>T SNPs, or their combinations, may be genetic diagnostic markers for racing performance indicators such as WRD and BRD.     

Bovine PGLYRP1 polymorphisms and their association with resistance to Mycobacterium avium ssp. paratuberculosis

Mycobacterium avium ssp. paratuberculosis (MAP) causes a chronic, granulomatous inflammatory condition of the intestines in ruminants and wild-type species. It causes significant economic losses to the dairy and beef industries owing to reduced productivity, premature culling and mortality. Bovine peptidoglycan recognition protein 1 is an important pattern recognition molecule that is capable of directly killing microorganisms. The goal of this study was to identify single nucleotide polymorphisms (SNPs) in the gene encoding bovine peptidoglycan recognition protein 1 and to assess their association with susceptibility to MAP infection in dairy cattle. Blood and milk samples were collected from Holsteins in Southwestern and Eastern Ontario and tested for MAP infection using blood and milk ELISAs. A resource population consisting of 197 infected (S/P > 0.25) and 242 healthy (S/P < 0.10) cattle was constructed. Sequencing of pooled DNA was used to identify three SNPs (c.102G>C, c.480G>A and c.625C>A) that were genotyped in the resource population. Statistical analysis was performed using a logistic regression model fitting the additive and dominance effects of each SNP in the model. SNP c.480G>A (P = 0.054) was found to be associated with susceptibility to MAP infection. Cows with a copy of the major allele 'G' at this locus had an odds ratio of 1.51 (95% CI: 0.99-2.31) for being infected with MAP.     

Whole‐Genome association analysis of susceptibility to paratuberculosis in Holstein cattle

The objective of this study was to identify genetic markers and genomic regions associated with susceptibility to Mycobacterium avium ssp. paratuberculosis (MAP) infection in Holstein cattle. Associated single nucleotide polymorphisms (SNPs) were identified by genotyping 521 MAP‐infected Holstein cows and comparing SNP allele frequencies of these infected cows with allele frequencies estimated from specific reference populations. Reference population allele frequency estimates used Holstein sire genotype data and were weighted estimates based on sire usage within the population in question. The 521 infected cows were 233 and 288 cows from two resource populations of approximately 5000 cows each, collected independently. Population 1 was comprised primarily of daughters of twelve Holstein artificial insemination sires used heavily within the US dairy cattle population. Samples were obtained from 300 co‐operating commercial dairy herds throughout the US and were tested by both MAP faecal culture and blood‐enzyme‐linked immunosorbent assay (ELISA). Population 2 consisted of dairy cattle from six co‐operating dairy herds in Wisconsin, with all animals in the herds tested by blood enzyme‐linked immunosorbent assay (ELISA) for MAP infection. Genotyping was performed with the Illumina Bovine SNP50 Bead Chip, providing genotypes for 35 772 informative SNPs. Data from the two resource populations were analysed both in separate and combined analyses. The most significant autosomal markers from the individual and combined analyses (n = 197, nominal P < 0.001) were used in a stepwise logistic regression analysis to identify a set of 51 SNPs that could be used as a predictor of genetics for Holstein cattle susceptibility to MAP infection.     

The challenge of separating signatures of local adaptation from those of isolation by distance and colonization history: The case of two white pines

Accurately detecting signatures of local adaptation using genetic‐environment associations (GEAs) requires controlling for neutral patterns of population structure to reduce the risk of false positives. However, a high degree of collinearity between climatic gradients and neutral population structure can greatly reduce power, and the performance of GEA methods in such case is rarely evaluated in empirical studies. In this study, we attempted to disentangle the effects of local adaptation and isolation by environment (IBE) from those of isolation by distance (IBD) and isolation by colonization from glacial refugia (IBC) using range‐wide samples in two white pine species. For this, SNPs from 168 genes, including 52 candidate genes for growth and phenology, were genotyped in 133 and 61 populations of Pinus strobus and P. monticola, respectively. For P. strobus and using all 153 SNPs, climate (IBE) did not significantly explained among‐population variation when controlling for IBD and IBC in redundancy analyses (RDAs). However, 26 SNPs were significantly associated with climate in single‐locus GEA analyses (Bayenv2 and LFMM), suggesting that local adaptation took place in the presence of high gene flow. For P. monticola, we found no evidence of IBE using RDAs and weaker signatures of local adaptation using GEA and F_ST outlier tests, consistent with adaptation via phenotypic plasticity. In both species, the majority of the explained among‐population variation (69 to 96%) could not be partitioned between the effects of IBE, IBD, and IBC. GEA methods can account differently for this confounded variation, and this could explain the small overlap of SNPs detected between Bayenv2 and LFMM. Our study illustrates the inherent difficulty of taking into account neutral structure in natural populations and the importance of sampling designs that maximize climatic variation, while minimizing collinearity between climatic gradients and neutral structure.     

Association of FTO Gene Variants With Adiposity in African‐American Adolescents

The prevalence of obesity continues to increase significantly, with the largest rise in the African‐American adolescents. Genetic contributions to obesity are being identified with the advent of genome‐wide association studies (GWAS). Specifically, variants of the fat mass and obesity associated (FTO) gene have been associated with obesity in populations of European descent. The studies in African Americans have been inconclusive. To further evaluate the association of the FTO gene and adiposity in African Americans, we genotyped 47 single‐nucleotide polymorphisms (SNPs), including seven SNPs previously reported to be significant in the literature in a cohort consisting of 561 non‐Hispanic white and 497 African‐American individuals. Analysis of our data showed 17 SNPs to be associated with BMI Z‐score (BMI‐Z) in our study population. The strongest association was found in the African Americans. The most significant SNP was rs8057044, which was associated with BMI‐Z in the African Americans (P = 0.00054). SNP rs9939609 was found to be significant in the non‐Hispanic white population (P = 0.028). Our data confirm the association between FTO and adiposity suggesting that FTO is a childhood obesity susceptibility gene. Our data also identify a novel SNP of the FTO gene (rs8057044) that is associated with measures of adiposity in the African‐American population.     

Genetic polymorphisms of the bovine Fatty acid binding protein 4 gene are significantly associated with marbling and carcass weight in Hanwoo (Korean Cattle)

The objective of this study was to investigate an association between polymorphisms in the FABP4 gene and phenotypic variation for marbling and carcass weight (CWT) in a population of Hanwoo steers. We re‐sequenced 4.3 kb of the FABP4 gene region in 24 Hanwoo bulls and identified 16 SNPs and 1 microsatellite polymorphism. Of these 16 SNPs, three SNPs [g.2774G>C (intron I), g.3473A>T (intron II) and g.3631G>A (exon III, creating a p.Met >Val amino acid substitution)] were genotyped in 583 steers to assess their association with carcass traits. The g.3473A allele showed a significant increasing effect on CWT (P = 0.01) and the g.3631G allele was associated with higher marbling score (P = 0.006). One haplotype of these three SNPs (CAG) was significantly associated with CWT (P = 0.02) and marbling score (P = 0.05) and could potentially be of value for marker assisted selection in Hanwoo cattle. The CAG haplotype effect for CWT was larger (11.14 ± 5.03 kg) than the largest single locus effect of g.3473A>T (5.01 ± 2.2 kg).     

Single‐nucleotide Polymorphism of CD36 Locus and Obesity in European Adolescents

CD36 is a membrane receptor with a wide variety of functions, including the regulation of energy metabolism, fat storage, and adipocyte differentiation. To assess the relationship between CD36 gene single‐nucleotide polymorphisms (SNPs) and obesity in adolescents, we evaluated the relationship between CD36 SNPs and the risk of obesity in a case–control study composed of 307 obese (age = 15.0 ± 1.1 years) and 339 normal‐weight adolescents (age = 14.6 ± 1.1 years). To validate the results, we assessed the relation between the same SNPs and percentage of body fat (BF%) and BMI in 1,151 European adolescents (age = 14.8 ± 1.4 years). SNPs with a minor allele frequency >0.10 were selected to tag CD36. Genotyping was performed on an Illumina system. Four SNPs (rs3211867, rs3211883, rs3211908, and rs1527483) were associated with increased risk of obesity in the case–control study (odds ratio (OR) (95% confidence interval)): 1.96 (1.26–3.04], P = 0.003; 1.73 (1.16–2.59), P = 0.007; 2.42 (1.47–4.01), P = 0.0005 and 1.95 (1.25–3.05), P = 0.003, respectively). The same four SNPs were associated with higher BMI (P < 0.05) and BF% (P < 0.04) in the validation study. Further analyses identified a haplotype (frequency: 0.05) carrying the minor allele of these SNPs as being associated with obesity (OR: 2.28; P = 0.0008) in the case–control study and with excess adiposity (i.e., higher BF% (P = 0.03) and BMI (P = 0.04)) in the validation study. Our data suggest that genetic variability at the CD36 gene locus could be associated with body weight variability in European adolescents but these findings require replication.     

Associations between newly discovered polymorphisms in the Bos taurusgrowth hormone receptor gene and performance traits in Holstein–Friesian dairy cattle

Variations in the growth hormone receptor (GHR) gene sequence are associated with performance traits in cattle. For example, the single nucleotide polymorphism (SNP) F279Y in transmembrane exon 8 has a strong association with milk yield. In this study, 32 previously unreported, putative novel SNPs (31 in the 5′ non‐coding region) were identified by resequencing ～19 kb of the GHR gene in genomic DNA from 22 cattle of multiple breeds. A population of 848 Holstein–Friesian AI sires was subsequently genotyped for the 32 putative novel SNPs and seven published SNPs (including F279Y, one in exon 1A promoter and five in exon 10). Associations between each segregating SNP and genetic merit for performance were quantified in the 848 Holstein–Friesians using weighted animal linear mixed models. Six of the published SNPs and seven of the novel SNPs were associated with at least one of the traits – milk yield, fat yield, protein yield, fat percentage, protein percentage, somatic cell score, calving interval, survival and growth and size traits. Even when the allelic substitution effect (P < 0.001) of F279Y was accounted for, the allelic substitution effect of one of the novel SNPs (GHR4.2) in the 5′ non‐coding region of GHR was associated with a lactation milk yield of 37.46 kg (P < 0.001). GHR4.2 and F279Y were not in linkage disequilibrium (r² = 0.00, D’ = 0.04) in the 848 Holstein–Friesians, indicating that their association with milk yield was independent.     

Identification of single nucleotide polymorphisms in the bovine follicle‐stimulating hormone receptor and effects of genotypes on superovulatory response traits

In dairy cows, there is evidence that failure to respond to superovulation protocols is a heritable trait. In women, genotyping for the p.N680S single nucleotide polymorphism (SNP) in the follicle‐stimulating hormone receptor (FSHR) gene may help identify poor responders before ovarian stimulation is initiated. Our objectives were to identify SNPs in the coding region of the bovine FSHR gene and to investigate the effect of FSHR genotypes on superovulatory response in Holstein cattle. Sequencing of FSHR exons 1–10 revealed seven SNPs. Three were non‐synonymous mutations (c.337C>G, c.871A>G and c.1973C>G). SNP c.337C>G encodes for a proline‐to‐alanine (p.Pro113Ala) amino acid replacement in the extracellular ligand‐binding domain of the receptor. PCR‐RFLP analyses showed that homozygous GG Holstein cows present a higher percentage of viable embryos, whereas GG and CG animals have less unfertilised oocytes. SNP c.871A>G results in an isoleucine‐to‐valine (p.Ile291Val) modification, and homozygous AA animals present lower embryo yield after superovulatory treatments. SNP c.1973C>G corresponds to a threonine‐to‐serine (p.The658Ser) modification in the intracellular carboxyl‐terminal domain of the FSHR protein, and homozygous GG Holstein cows were associated with a lower embryo yield and a higher percentage of unfertilised oocytes. Our results suggest that specific alleles of the bovine FSHR gene are associated with variations in embryo yield and in the number of unfertilised oocytes.     

Osteogenic and neurogenic differentiation of EGFP gene transfected neural stem cells derived from the brain of porcine fetuses at intermediate and late gestational age

The aims of this study were (i) to determine whether NSCs (neural stem cells) could be isolated from the brain of porcine fetuses at intermediate and late gestational age and (ii) to determine if these stem cells could be differentiated in vitro into osteogenic and neurogenic lineages following transfection with a reporter gene, EGFP (enhanced green fluorescence protein). The NSCs were isolated from the brains of porcine fetuses at intermediate and late gestational age and transfected with EGFP gene using lipofection. The transfected NSCs cells were induced to differentiate into cells of osteogenic and neurogenic lineages. Markers associated with NSCs and their osteogenic and neurogenic derivatives were tested by PCR. The results demonstrated that NSCs could be isolated from the brain of porcine fetus at intermediate and late gestational age and that transfected NSCs expressed EGFP and could be induced to differentiate in vitro. NSCs expressed CD‐90, Hes1, Oct4, Sox2 and Nestin, while following differentiation cells expressed markers for osteogenic (osteocalcin and osteonectin) and neurogenic cells such as astrocyte [GFAP (glial fibrillary acidic protein)], oligodendrocyte [GALC (galactosylceramide)] and neuron [NF (neurofilament), ENO2 (enolase 2) and MAP (microtubule‐associated protein)].     

Association of the variants in the BUD13‐ZNF259 genes and the risk of hyperlipidaemia

The single nucleotide polymorphisms (SNPs) in the BUD13 homolog (BUD13) and zinc finger protein 259 (ZNF259) genes have been associated with one or more serum lipid traits in the European populations. However, little is known about such association in the Chinese populations. Our objectives were to determine the association of the BUD13/ZNF259 SNPs and their haplotypes with hypercholesterolaemia (HCH)/hypertriglyceridaemia (HTG) and to identify the possible gene–gene interactions among these SNPs. Genotyping of 6 SNPs was performed in 634 hyperlipidaemic and 547 normolipidaemic participants. The ZNF259 rs2075290, ZNF259 rs964184 and BUD13 rs10790162 SNPs were significantly associated with serum lipid levels in both HCH and non‐HCH populations (P < 0.008–0.001). On single locus analysis, only BUD13 rs10790162 was associated with HCH (OR: 2.23, 95% CI: 1.05, 4.75, P = 0.015). The G‐G‐A‐A‐C‐C haplotype, carrying rs964184‐G‐allele, was associated with increased risk of HCH (OR: 1.35, 95% CI: 1.10, 1.66, P = 0.005) and HTG (OR: 1.75, 95% CI: 1.39, 2.21, P = 0.000). The A‐C‐G‐G‐C‐C and A‐C‐A‐G‐T‐C haplotypes, carrying rs964184‐C‐allele, were associated with reduced risk of HCH (OR: 0.77, 95% CI: 0.61, 0.99, P = 0.039 and OR: 0.66, 95% CI: 0.47, 0.94, P = 0.021 respectively). On multifactor dimensionality reduction analyses, the two‐ to three‐locus models showed a significant association with HCH and HTG (P < 0.01–0.001). The BUD13/ZNF259 SNPs, which were significant in the European populations, are also replicable in the Southern Chinese population. Moreover, inter‐locus interactions may exist among these SNPs. However, further functional studies are required to clarify how these SNPs and genes actually affect the serum lipid levels.     

Evaluation of PAX3 genetic variants and nevus number

The presence of a high nevus number is the strongest phenotypic predictor of melanoma risk. Here, we describe the results of a three‐stage study directed at identifying risk variants for the high nevus phenotype. At the first stage, 263 melanoma cases from Barcelona were genotyped for 223 single‐nucleotide polymorphisms (SNPs) in 39 candidate genes. Seven SNPs in the PAX3 gene were found to be significantly associated with nevus number under the additive model. Next, the associations for seven PAX3 variants were evaluated in 1217 melanoma cases and 475 controls from Leeds; and in 3054 healthy twins from TwinsUK. Associations with high nevus number were detected for rs6754024 (P values < 0.01) in the Barcelona and Leeds datasets and for rs2855268 (P values < 0.01) in the Barcelona and the TwinsUK sets. Associations (P values < 0.001) in the opposite direction were detected for rs7600206 and rs12995399 in the Barcelona and TwinsUK sets. This study suggests that SNPs in PAX3 are associated with nevus number, providing support for PAX3 as a candidate nevus gene. Further studies are needed to examine the role of PAX3 in melanoma susceptibility.     

A case-control study between the gene polymorphisms of polyunsaturated fatty acids metabolic rate-limiting enzymes and coronary artery disease in a Chinese Han population

To investigate the association between the polymorphisms of fatty acid desaturase 1 (FADS1), fatty acid desaturase 2 (FADS2) and elongation of very long chain fatty acids like 2 (ELOVL2) gene and coronary artery disease (CAD) in a Chinese Han population. Three single nucleotide polymorphisms (SNPs) from these genes were genotyped using PCR-based restriction fragment length polymorphism analysis in 199 CAD cases and 192 controls of Han Chinese origin. rs174556 in the FADS1 gene showed allelic (P=0.002) and genotypic (P=0.030) association with the disease, while there was no disease association for the other two SNPs. The frequency of rs174556 minor allele (T) was significantly higher in the case group than the control group. The trans phase gene–gene interaction analysis showed that the combined genotype of rs174556 (T/T) and rs3756963 (T/T) was weakly associated with the disease (P=0.043). rs174556 in the FADS1 gene is very likely to be associated with CAD in the Chinese Han population.     

Association analyses of porcine SERPINE1 reveal sex-specific effects on muscling, growth, fat accretion and meat quality

The serpin peptidase inhibitor, clade E (nexin, plasminogen activator inhibitor type 1), member 1 (SERPINE1) gene encodes plasminogen activator inhibitor type 1 (PAI), which is the major physiological inhibitor of tissue‐type and urokinase‐type plasminogen activators and plays a role in obesity and insulin resistance in women but not in men. We detected SNP FN396538:g.566G>A in intron 3 and a non‐synonymous substitution NM_213910:c.612A>G in exon 3 (p.Ile159Val) and mapped the gene to position 8.4 cM on the linkage map of chromosome 3. Association analyses were conducted on the 12th–15th generation of the Meishan × Large White (MLW) cross (n = 565), with records for weight at the end of test, lifetime daily gain, test time daily gain, loin depth and backfat depth, as well as on a European wild boar × Meishan (W × M) F₂ population (n = 333) with 47 traits recorded for carcass composition and meat quality. Analyses performed across the entire MLW population or in the male animals did not show any trait significantly associated with the loci studied. In female animals, both SNPs were associated with loin depth at nominal P < 0.05 with adjusted P values equal to 0.051 (g.566) and 0.057 (c.612). Differences between homozygotes were up to 0.65 SD. In the entire W × M population and female animals, SERPINE1 was significantly associated at adjusted P < 0.05 in descending order with muscling, growth and fat accretion and in male animals with meat quality (R‐value). In the studied populations, allele effects were in opposite directions, which implies that the SNPs are markers that are in linkage disequilibrium with a causative mutation.     

Galling and reversion disease incidence in a range of blackcurrant genotypes, differing in resistance to the blackcurrant gall mite (Cecidophyopsis ribis) and blackcurrant reversion disease

An assessment was made of the response of 10 blackcurrant genotypes, differing in resistance to the agent of blackcurrant reversion disease (BRD) and to its gall mite vector, in field trials in Scotland and Finland over 5 yr. At each location, the 10 genotypes were planted in plots containing infector plants with high inoculum levels of the two organisms. In Scotland, the infector plants contained large numbers of gall mites (Cecidophyopsis ribis) and were infected with the European (E) form of BRD; in Finland, infector plants contained a different species of gall mite (C. spicata) and the severe Russian form (R) of BRD. As expected, at both sites, almost all plants of cvs Ben Alder, Ben Lomond, Ben Tirran, Ojebyn and an SCRI selection F4/1/66, which are susceptible to gall mite and BRD, became infested with each of these organisms. However, in Scotland but not in Finland, 58% of cv. Ojebyn plants were affected by BRD. The cv. Foxendown, which contains gene Ce that confers apparent immunity to C. ribis, was free from galls and failed to develop distinctive BRD symptoms at both sites. The cvs Rus and Neosypajuscajaija, which contain gene P which is reported to confer resistance to C. ribis, were infested more slowly by mites than the mite-susceptible genotypes and showed a smaller number of galls per plant. Also, they were infected with BRD more slowly than some mite-susceptible genotypes although by the end of the experiment, most plants were affected by BRD. All plants of cvs Golubka and Ben Gairn, which are resistant to the agent of BRD, remained free from distinct BRD symptoms at both sites despite the fact that plants contained galls. These data indicate the superiority of gene Ce over gene P for resistance to gall mites with the added benefit that the virtual immunity to gall mites provided by gene Ce provides a high level of protection against infection with BRD. The relative merits of these different forms of resistance to gall mites and BRD in blackcurrant are discussed in relation to present control methods for these two organisms and in the light of recent findings of different species of eriophyid mites on Ribes species.     

Genotype‐by‐sequencing of the plant‐pathogenic fungi Pyrenophora teres and Sphaerulina musiva utilizing Ion Torrent sequence technology

Genetic and genomics tools to characterize host–pathogen interactions are disproportionately directed to the host because of the focus on resistance. However, understanding the genetics of pathogen virulence is equally important and has been limited by the high cost of de novo genotyping of species with limited marker data. Non‐resource‐prohibitive methods that overcome the limitation of genotyping are now available through genotype‐by‐sequencing (GBS). The use of a two‐enzyme restriction‐associated DNA (RAD)‐GBS method adapted for Ion Torrent sequencing technology provided robust and reproducible high‐density genotyping of several fungal species. A total of 5783 and 2373 unique loci, ‘sequence tags’, containing 16 441 and 9992 single nucleotide polymorphisms (SNPs) were identified and characterized from natural populations of Pyrenophora teres f. maculata and Sphaerulina musiva, respectively. The data generated from the P. teres f. maculata natural population were used in association mapping analysis to map the mating‐type gene to high resolution. To further validate the methodology, a biparental population of P. teres f. teres, previously used to develop a genetic map utilizing simple sequence repeat (SSR) and amplified fragment length polymorphism (AFLP) markers, was re‐analysed using the SNP markers generated from this protocol. A robust genetic map containing 1393 SNPs on 997 sequence tags spread across 15 linkage groups with anchored reference markers was generated from the P. teres f. teres biparental population. The robust high‐density markers generated using this protocol will allow positional cloning in biparental fungal populations, association mapping of natural fungal populations and population genetics studies.