Genome‐wide association analysis and gene set enrichment analysis with SNP data identify genes associated with 305‐day milk yield in Holstein dairy cows

Milk production traits, such as 305‐day milk yield (305MY), have been under direct selection to improve production in dairy cows. Over the past 50 years, the average milk yield has nearly doubled, and over 56% of the increase is attributable to genetic improvement. As such, additional improvements in milk yield are still possible as new loci are identified. The objectives of this study were to detect SNPs and gene sets associated with 305MY in order to identify new candidate genes contributing to variation in milk production. A population of 781 primiparous Holstein cows from six central Washington dairies with records of 305MY and energy corrected milk were used to perform a genome‐wide association analysis (GWAA) using the Illumina BovineHD BeadChip (777 962 SNPs) to identify QTL associated with 305MY (P < 1.0 × 10^?5). A gene set enrichment analysis with SNP data (GSEA‐SNP) was performed to identify gene sets (normalized enrichment score > 3.0) and leading edge genes (LEGs) influencing 305MY. The GWAA identified three QTL comprising 34 SNPs and 30 positional candidate genes. In the GSEA‐SNP, five gene sets with 58 unique and 24 shared LEGs contributed to 305MY. Identification of QTL and LEGs associated with 305MY can provide additional targets for genomic selection to continue to improve 305MY in dairy cattle.     

Validating genome‐wide associated signals for clinical mastitis in German Holstein cattle

A validation study for six genomic regions previously identified by a genome‐wide association study for somatic cell score was conducted with data of clinical mastitis in German Holstein cattle. Out of 10 tested SNPs, five on chromosomes 6, 13 and 19 were significantly associated with clinical mastitis (P < 0.05). Three SNPs on chromosomes 6 and 19 had the same direction of effect as those previously reported in the initial genome‐wide association study for somatic cell score. The other two SNPs on chromosome 13 had opposite effects. As well as validating associations within known QTL from previous studies, e.g. chromosomes 6 and 19, novel loci on chromosome 13 were confirmed. Promising candidate genes are, for example: deoxycytidine kinase, immunoglobulin J chain, vitamin D binding protein, forkhead box K2, sodium/hydrogen exchanger 8 and cytoplasmic nuclear factor of activated T‐cells 2. Our confirmation study provides additional evidence for the functional role of the linked genomic regions to immune response. This information can be used as a basis for further functional studies for those potential genes.     

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.     

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.     

Effect of FASN gene on milk yield and milk composition in the Chinese Holstein dairy population

Fatty acid synthase (FASN) is a multifunctional protein that catalyzes de novo synthesis of fatty acids in cells. It plays a key role in the lipid biosynthesis as well as in the general metabolism of all living animals. We herein investigated polymorphisms of FASN. As a result, six single nucleotide polymorphisms (SNPs) were found and then genotyped in 752 Chinese Holstein cows. It was found that g.17924A>G was non‐synonymous, g.13965 C>T, g.16907 T>C and g.18663T>C were synonymous mutations and two other two SNPs, g.8948 C>T (ss491228481) and g.14439T>C (rs133498277), were in intronic sequences of the gene. All such identified SNPs were found to be associated with milk yield and composition traits (P = 0.0441 to <0.0001). Significant additive and allele substitution effects were observed for three yield traits at all six loci as well (P < 0.05 to <0.01). Complete linkage disequilibrium among the five SNPs, with the exception of g.8948 C>T, was observed.     

Antimicrobial resistance and toxin gene profiles of Staphylococcus aureus strains from Holstein milk

Isolation of Staphylococcus aureus (Staph. aureus) from Holstein milk samples with mastitis and nonmastitis was conducted to estimate its prevalence, antimicrobial resistance and toxin genes. A total of 353 milk samples were collected from three Chinese Holstein herds. Fifty‐three Staph. aureus isolates collected from 29 Staph. aureus‐positive samples were characterized via antimicrobial susceptibility, toxin genes and Pulsed‐field Gel Electrophoresis (PFGE) profiles. The prevalence of Staph. aureus was 4·0–9·5% in mastitic and 7·3–11·5% in nonmastitic samples in the analysed herds. Approximately 61·0% of Staph. aureus strains isolated from mastitis cows were resistant to ≥10 antimicrobials compared with 0% of isolates with nonmastitis. The most frequently observed super antigenic toxin gene was pvl (41·5%) followed by seh + pvl (13·2%). We did not find mecA‐positive methicillin‐resistant Staph. aureus (MRSA) strains, while mecA‐negative MRSA strains were identified in the three herds. PFGE results suggested potential transmission of Staph. aureus strains in different farms. These results open new insights into Staph. aureus transmission and antimicrobial resistance of Holstein dairy cows and into developing strategies for udder health improvement of dairy cattle.     

Genetic association between bovine NOD2 polymorphisms and infection by Mycobacterium avium subsp. paratuberculosis in Holstein‐Friesian cattle

Nucleotide‐Binding Oligomerization Domain 2 (NOD2) has been reported to be a candidate gene for Mycobacterium avium subsp. paratuberculosis (MAP) infection in a Bos taurus × Bos indicus mixed breed based on a genetic association with the c.2197T>C single nucleotide polymorphism (SNP). Nevertheless, this SNP has also been reported to be monomorphic in the B. taurus species. In the present work, 18 SNPs spanning the bovine NOD2 gene have been analysed in a genetic association study of two independent populations of Holstein‐Friesian cattle. We found that the C allele of SNP c.*1908C>T, located in the 3′‐UTR region of the gene, is significantly more frequent in infected animals than in healthy ones, which supports the idea that the bovine NOD2 gene plays a role in susceptibility to MAP infection. However, in silico analyses of the NOD2 nucleotide sequence did not yield definitive data about a possible direct effect of SNP c.*1908C>T on susceptibility to infection and led us to consider its linkage disequilibrium with the causative variant. A more exhaustive genetic association study including all putative, functional SNPs from this gene and subsequent functional analyses needs to be conducted to achieve a more complete understanding of how different variants of NOD2 may affect susceptibility to MAP infection in cattle.     

Genome‐wide association study for semen traits of the bulls in Chinese Holstein

A genome‐wide association study (GWAS) was performed to identify markers and candidate genes for five semen traits in the Holstein bull population in China. The analyzed dataset consisted of records from 692 bulls from eight bull stations; each bull was genotyped using the Illumina BovineSNP50 BeadChip. Association tests between each trait and the 41 188 informative high‐quality SNPs were achieved with gapit software. In total, 19 suggestive significant SNPs, partly located within the reported QTL regions or within or close to the reported candidate genes, associated with five semen traits were detected. By combining our GWAS results with the biological functions of these genes, eight novel promising candidate genes, including ETNK1, PDE3A, PDGFRB, CSF1R, WT1, DSCAML1, SOD1 and RUNX2, were identified that potentially relate to semen traits. Our findings may provide a basis for further research on the genetic mechanism of semen traits and marker‐assisted selection of such traits in Holstein bulls.     

SNPs in SNCA,MCCC1, DLG2, GBF1 and MBNL2 are associated with Parkinson's disease in southern Chinese population

Numerous single nucleotide polymorphisms (SNPs), which have been identified as susceptibility factors for Parkinson's disease (PD) as per genome‐wide association studies, have not been fully characterized for PD patients in China. This study aimed to replicate the relationship between 12 novel SNPs of 12 genes and PD risk in southern Chinese population. Twelve SNPs of 12 genes were detected in 231 PD patients and 249 controls, using the SNaPshot technique. Meta‐analysis was used to assess heterogeneity of effect sizes between this study and published data. The impact of SNPs on gene expression was investigated by analysing the SNP‐gene association in the expression quantitative trait loci (eQTL) data sets. rs8180209 of SNCA (allele model: P = .047, OR = 0.77; additive model: P = .047, OR = 0.77), rs2270968 of MCCC1 (dominant model: P = .024, OR = 1.52), rs7479949 of DLG2 (recessive model; P = .019, OR = 1.52), rs10748818 of GBF1 (additive model: P < .001, OR = 0.37), and rs4771268 of MBNL2 (recessive model: P = .003, OR = 0.48) were replicated to be significantly associated with the increased risk of PD. Noteworthy, a meta‐analysis of previous studies suggested rs8180209, rs2270968, rs7479949 and rs4771268 were in line with those of our cohort. Our study replicated five novel functional SNPs in SNCA, MCCC1, DLG2, GBF1 and MBNL2 could be associated with increased risk of PD in southern Chinese population.     

Genome-Wide Association Study of Susceptibility to Infection by Mycobacterium avium Subspecies paratuberculosis in Holstein Cattle

Paratuberculosis, or Johne''s disease, is a chronic, granulomatous, gastrointestinal tract disease of cattle and other ruminants caused by the bacterium Mycobacterium avium, subspecies paratuberculosis (MAP). Control of Johne''s disease is based on programs of testing and culling animals positive for infection with MAP while concurrently modifying management to reduce the likelihood of infection. The current study is motivated by the hypothesis that genetic variation in host susceptibility to MAP infection can be dissected and quantifiable associations with genetic markers identified. For this purpose, a case-control, genome-wide association study was conducted using US Holstein cattle phenotyped for MAP infection using a serum ELISA and/or fecal culture test. Cases included cows positive for either serum ELISA, fecal culture or both. Controls consisted of animals negative for the serum ELISA test or both serum ELISA and fecal culture when both were available. Controls were matched by herd and proximal birth date with cases. A total of 856 cows (451 cases and 405 controls) were used in initial discovery analyses, and an additional 263 cows (159 cases and 104 controls) from the same herds were used as a validation data set. Data were analyzed in a single marker analysis controlling for relatedness of individuals (GRAMMAR-GC) and also in a Bayesian analysis in which multiple marker effects were estimated simultaneously (GenSel). For the latter, effects of non-overlapping 1 Mb marker windows across the genome were estimated. Results from the two discovery analyses were generally concordant; however, discovery results were generally not well supported in analysis of the validation data set. A combined analysis of discovery and validation data sets provided strongest support for SNPs and 1 Mb windows on chromosomes 1, 2, 6, 7, 17 and 29.     

Identification of genetic associations of ECHS1 gene with milk fatty acid traits in dairy cattle

Our previous genome‐wide association study identified 83 genome‐wide significant SNPs and 20 novel promising candidate genes for milk fatty acids in Chinese Holstein. Among them, the enoyl‐CoA hydratase, short chain 1 (ECHS1) and enoyl‐CoA hydratase and 3‐hydroxyacyl CoA dehydrogenase (EHHADH) genes were located near two SNPs and one SNP respectively, and they play important roles in fatty acid metabolism pathways. We herein validated whether the two genes have genetic effects on milk fatty acid traits in dairy cattle. By re‐sequencing the full‐length coding region, partially adjacent introns and 3000 bp up/downstream flanking sequences, we identified 12 SNPs in ECHS1: two in exons, four in the 3′ flanking region and six in introns. The g.25858322C>T SNP results in an amino acid replacement from leucine to phenylalanine and changes the secondary structure of the ECHS1 protein, and single‐locus association analysis showed that it was significantly associated with three milk fatty acids (P = 0.0002–0.0013). The remaining 11 SNPs were found to be significantly associated with at least one milk fatty acid (P = <0.0001–0.0040). Also, we found that two haplotype blocks, consisting of nine and two SNPs respectively, were significantly associated with eight milk fatty acids (P = <0.0001–0.0125). However, none of polymorphisms was observed in the EHHADH gene. In conclusion, our findings are the first to indicate that the ECHS1 gene has a significant genetic impact on long‐chain unsaturated and medium‐chain saturated fatty acid traits in dairy cattle, although the biological mechanism is still undetermined and requires further in‐depth validation.     

Candidate gene analysis of osteochondrosis in Spanish Purebred horses

Equine osteochondrosis (OC) is a frequent developmental orthopaedic disease with high economic impact on the equine industry and may lead to premature retirement of the animal as a result of chronic pain and lameness. The genetic background of OC includes different genes affecting several locations; however, these genetic associations have been tested in only one or few populations, lacking the validation in others. The aim of this study was to identify the genetic determinants of OC in the Spanish Purebred horse breed. For that purpose, we used a candidate gene approach to study the association between loci previously implicated in the onset and development of OC in other breeds and different OC locations using radiographic data from 144 individuals belonging to the Spanish Purebred horse breed. Of the 48 polymorphisms analysed, three single nucleotide polymorphisms (SNPs) located in the FAF1, FCN3 and COL1A2 genes were found to be associated with different locations of OC lesions. These data contribute insights into the complex gene networks underlying the multifactorial disease OC, and the associated SNPs could be used in a marker‐assisted selection strategy to improve horse health, welfare and competitive lifespan.     

Polymorphisms in the AKT1 and AKT2 genes and oesophageal squamous cell carcinoma risk in an Eastern Chinese population

Ethnic Han Chinese are at high risk of developing oesophageal squamous cell carcinoma (ESCC). Aberrant activation of the AKT signalling pathway is involved in many cancers, including ESCC. Some single nucleotide polymorphisms (SNPs) in genes involved in this pathway may contribute to ESCC susceptibility. We selected five potentially functional SNPs in AKT1 (rs2494750, rs2494752 and rs10138277) and AKT2 (rs7254617 and rs2304186) genes and investigated their associations with ESCC risk in 1117 ESCC cases and 1096 controls in an Eastern Chinese population. None of individual SNPs exhibited an association with ESCC risk. However, the combined analysis of three AKT1 SNPs suggested that individuals carrying one of AKT1 variant genotypes had a decreased ESCC risk [adjusted odds ratio (OR) = 0.60, 95% CI = 0.42–0.87]. Further stratified analysis found that AKT1 rs2294750 SNP was associated with significantly decreased ESCC risk among women (adjusted OR = 0.63, 95% CI = 0.43–0.94) and non‐drinkers (OR = 0.79, 95% CI = 0.64–0.99). Similar protective effects on women (adjusted OR = 0.56, 95% CI = 0.37–0.83) and non‐drinker (adjusted OR = 0.75, 95% CI = 0.60–0.94) were also observed for the combined genotypes of AKT1 SNPs. Consistently, logistic regression analysis indicated significant gene–gene interactions among three AKT1 SNPs (P < 0.015). A three‐AKT1 SNP haplotype (C‐A‐C) showed a significant association with a decreased ESCC risk (adjusted OR = 0.70, 95% CI = 0.52–0.94). Multifactor dimensionality reduction analysis confirmed a high‐order gene–environment interaction in ESCC risk. Overall, we found that three AKT1 SNPs might confer protection against ESCC risk; nevertheless, these effects may be dependent on other risk factors. Our results provided evidence of important gene–environment interplay in ESCC carcinogenesis.     

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.     

Combining next‐generation sequencing and progeny testing for rapid identification of induced recessive and dominant mutations in maize M2 individuals

Molecular identification of mutant alleles responsible for certain phenotypic alterations is a central goal of genetic analyses. In this study we describe a rapid procedure suitable for the identification of induced recessive and dominant mutations applied to two Zea mays mutants expressing a dwarf and a pale green phenotype, respectively, which were obtained through pollen ethyl methanesulfonate (EMS) mutagenesis. First, without prior backcrossing, induced mutations (single nucleotide polymorphisms, SNPs) segregating in a (M₂) family derived from a heterozygous (M₁) parent were identified using whole‐genome shotgun (WGS) sequencing of a small number of (M₂) individuals with mutant and wild‐type phenotypes. Second, the state of zygosity of the mutation causing the phenotype was determined for each sequenced individual by phenotypic segregation analysis of the self‐pollinated (M₃) offspring. Finally, we filtered for segregating EMS‐induced SNPs whose state of zygosity matched the determined state of zygosity of the mutant locus in each sequenced (M₂) individuals. Through this procedure, combining sequencing of individuals and Mendelian inheritance, three and four SNPs in linkage passed our zygosity filter for the homozygous dwarf and heterozygous pale green mutation, respectively. The dwarf mutation was found to be allelic to the an1 locus and caused by an insertion in the largest exon of the AN1 gene. The pale green mutation affected the nuclear W2 gene and was caused by a non‐synonymous amino acid exchange in encoded chloroplast DNA polymerase with a predicted deleterious effect. This coincided with lower cpDNA levels in pale green plants.     

A genome‐wide association study in Caucasian women suggests the involvement of HLA genes in the severity of facial solar lentigines

Solar lentigines are a common feature of sun‐induced skin ageing. Little is known, however, about the genetic factors contributing to their development. In this genome‐wide association study, we aimed to identify genetic loci associated with solar lentigines on the face in 502 middle‐aged French women. Nine SNPs, gathered in two independent blocks on chromosome 6, exhibited a false discovery rate below 25% when looking for associations with the facial lentigine score. The first block, in the 6p22 region, corresponded to intergenic SNPs and also exhibited a significant association with forehead lentigines (P = 1.37 × 10^?8). The second block, within the 6p21 HLA region, was associated with decreased HLA‐C expression according to several eQTL databases. Interestingly, these SNPs were also in high linkage disequilibrium with the HLA‐C*0701 allele (r² = 0.95). We replicated an association recently found by GWAS in the IRF4 gene. Finally, a complementary study on 44 selected candidate SNPs revealed novel associations in the MITF gene. Overall, our results point to several mechanisms involved in the severity of facial lentigines, including HLA/immunity and the melanogenesis pathway.     

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.     

Polymorphisms of the bovine luteinizing hormone/choriogonadotropin receptor (LHCGR) gene and its association with superovulation traits

The major limitation to the development of embryo transfer technique in cattle is the highly variable between individuals in ovulatory response to FSH-induced superovulation. The objective of this study was to identify a predictor to forecast superovulation response on the basis of associations between superovulation performance and gene polymorphism, variation in the bovine luteinizing hormone/choriogonadotropin receptor (LHCGR) gene was investigated using PCR-single-strand conformational (PCR-SSCP) and DNA sequencing. Four single nucleotide polymorphisms (SNPs) of G51656T, A51703G, A51726G and G51737A were identified at the intron 9 of the LHCGR gene in 171 Chinese Holstein cows treated for superovulation, and evaluated its associations with superovulatory response. Association analysis showed that these four SNPs had significant effects on the total number of ova (TNO) (P < 0.05). Moreover, the A51703G and A51726G polymorphisms significantly associated with the number of transferable embryos (NTE) (P < 0.05). In addition, significant additive effect on TNO was detected in polymorphisms of G51656T (P < 0.05) and A51703G (P < 0.01), and the A51703G polymorphism also had significant additive effects on NTE (P < 0.01). These results indicate that LHCGR gene is a potential marker for superovulation response and can be used to predict the most appropriate dose of FSH for superovulation in Chinese Holstein cows.