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 study for energy balance and fat/protein ratio in German Holstein bull dams

The metabolic status of cows is important to health and fertility, especially in early lactation, and energy balance (EB) and fat/protein ratio (FPR) are considered as appropriate indicators for metabolic disorders. The aim of this study was to detect SNPs (single nucleotide polymorphisms) associated with EB and FPR in German Holstein bull dams belonging to the research herd Karkendamm. Bull dams were genotyped using the Illumina Bovine SNP 50K Bead chip^® comprising 54 001 SNPs. A total of 43 593 SNPs and 586 (EB) and 668 (FPR) bull dams passed the quality control criteria. Phenotypes were deregressed breeding values estimated via random regression animal models for lactation days 11, 20, 30, and 42 for EB and FPR. Whole‐genome association analyses were carried out fitting principal components as covariates to adjust for genetic substructure. Permutation tests were applied to estimate genome‐wise significance. Across all observed lactation days, 19 SNPs located in four different intervals on chromosomes 1, 14, 16, and 27 were detected. For EB, seven markers across four chromosomes were identified. There was no overlap between markers associated with FPR and EB. SNPs associated with FPR were mostly located in QTL regions for milk production traits, especially in the region of DGAT1, whereas SNPs associated with EB mainly showed positional relationships to previously described QTL regions affecting functional traits, especially fertility.     

LINE-1 and EPAS1 DNA methylation associations with high-altitude exposure

Recent discoveries indicate a genetic basis for high-altitude adaptation among human groups who have resided at high altitude for millennia, including Andeans, Tibetans, and Ethiopians. Yet, genetics alone does not explain the extent of variation in altitude-adaptive phenotypes. Current and past environments may also play a role, and one way to determine the effect of the environment is through the epigenome. To characterize if Andean adaptive responses to high altitude have an epigenetic component, we analyzed DNA methylation of the promoter region of EPAS1 and LINE-1 repetitive element among 572 Quechua individuals from high- (4,388 m) and low-altitude (0 m) in Peru. Participants recruited at high altitude had lower EPAS1 DNA methylation and higher LINE-1 methylation. Altitude of birth was associated with higher LINE-1 methylation, not with EPAS1 methylation. The number of years lived at high altitude was negatively associated with EPAS1 methylation and positively associated with LINE-1 methylation. We found four one-carbon metabolism SNPs (MTHFD1 rs2236225, TYMS rs502396, FOLH1 rs202676, GLDC rs10975681) that cumulatively explained 11.29% of the variation in average LINE-1 methylation. And identified an association between LINE-1 methylation and genome-wide SNP principal component 1 that distinguishes European from Indigenous American ancestry suggesting that European admixture decreases LINE-1 methylation. Our results indicate that both current and lifetime exposure to high-altitude hypoxia have an effect on EPAS1 and LINE-1 methylation among Andean Quechua, suggesting that epigenetic modifications may play a role in high-altitude adaptation.     

Genes with expression levels correlating to drip loss prove association of their polymorphism with water holding capacity of pork

Six genes that were known to exhibit expression levels that are correlated to drip loss BVES, SLC3A2, ZDHHC5, CS, COQ9, and EGFR have been for candidate gene analysis. Based on in silico analysis SNPs were detected, confirmed by sequencing, and used for genotyping. The SNPs were genotyped in about 1,800 animals from six pig populations including commercial herds of Pietrain (PI) and German Landrace (DL), different commercial herds of Pietrain × (German Large White × German Landrace) (PIF1(a/b/c)), and one experimental F₂-population Duroc × Pietrain (DUPI). Comparative and genetic mapping established the location of BVES on SSC1, of SLC3A2 and ZDHHC5 on SSC2, of CS on SSC5, of COQ9 on SSC6 and of EGFR on SSC9, respectively, coinciding with QTL regions for carcass and meat quality traits. BVES, SLC3A2, and CS revealed association at least with drip loss and with several other measures of water holding capacity (WHC). Moreover, COQ9 and EGFR were associated with several meat quality traits such as meat color and/or thawing loss. This study reveals statistic evidence in addition to the functional relationship of these genes to WHC previously evidenced by expression analysis. This study reveals positional and genetic statistical evidence for a link of genetic variation at these loci or close to them and promotes those six candidate genes as functional and/or positional candidate genes for meat quality traits.     

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.     

A genome‐wide association study for the incidence of persistent bovine viral diarrhea virus infection in cattle

Bovine viral diarrhea viruses (BVDV) comprise a diverse group of viruses that cause disease in cattle. BVDV may establish both transient and persistent infections depending on the developmental stage of the animal at exposure. The objective was to determine whether genomic regions harboring single nucleotide polymorphisms (SNPs) could be associated with the presence or absence of persistent BVDV infection. A genome‐wide association approach based on 777 000 SNP markers was used. Samples of animals identified as positive (n = 1200) or negative (n = 1200) for the presence of BVDV in skin samples (n = 1200) were used. DNA samples were combined in 24 pools (100 animals per pool). One SNP, significant at the 5 percent genome‐wide level (P = 9.41 × 10^?8), was detected on chromosome 14, located at position 80 675 176 bp. Fifteen SNPs, residing on chromosomes 1, 2, 6, 8, 10, 15 and 18, were moderately associated (P < 1 × 10^?5) with persistent BVDV infection. Results show that genes harboring or neighboring significant SNPs are involved in leucopenia, signal transduction, RNA splicing and DNA methylation processes.     

The EPAS1 gene influences the aerobic–anaerobic contribution in elite endurance athletes

Absract EPAS1 is a gene involved in complex oxygen sensing. It is expressed in microvascular endothelial cells, lung epithelial cells, cardiac myocytes and the brain. An association study was undertaken comparing elite endurance athletes classified into two groups according to a power–time model of performance intensity: power–time-maximum (PT-MAX; N=242, event duration 50 s to 10 min) and power–time–steady state (PT-SS; N=151, event duration ~2–10 h), with normal controls (N=444) using 12 SNPs across EPAS1. Ordinal regression analysis of allele frequencies revealed significant differences at SNPs 2 and 3 (P=0.01). Haplotype analysis revealed the presence of haplotypes involving SNPs 2–5 that significantly differentiated (P<0.05) the groups based on an ordinal ranking using the power–time classification. These same haplotypes differentiated the PT-MAX group in which a significant decrease in a haplotype (F: G-C-C-G; OR=0.57, P=0.02, 95% CI 0.36–0.92) and increase in a second haplotype (G: A-T-G-G; OR=1.75, P=0.03, 95% CI 1.05–2.91) was observed compared to controls. The PT-SS group was differentiated from the PT-MAX group by a third haplotype (H: A-T-G-A; OR=0.46, P=0.04, 95% CI 0.22–0.96). Since EPAS1 has a role as a sensor capable of integrating cardiovascular function, energetic demand, muscle activity and oxygen availability into physiological adaptation, we propose that DNA variants in EPAS1 influence the relative contribution of aerobic and anaerobic metabolism and hence the maximum sustainable metabolic power for a given event duration.     

Genetic dissection of reproductive performance of dairy cows under heat stress

Heat stress negatively impacts the reproductive performance of dairy cows. The main objective of this study was to dissect the genetic basis underlying dairy cow fertility under heat stress conditions. Our first goal was to estimate genetic components of cow conception across lactations considering heat stress. Our second goal was to reveal individual genes and functional gene-sets that explain a cow’s ability to conceive under thermal stress. Data consisted of 74 221 insemination records on 13 704 Holstein cows. Multitrait linear repeatability test-day models with random regressions on a function of temperature–humidity index values were used for the analyses. Heritability estimates for cow conception under heat stress were around 2–3%, whereas genetic correlations between general and thermotolerance additive genetic effects were negative and ranged between −0.35 and −0.82, indicating an unfavorable relationship between cows’ ability to conceive under thermo-neutral vs. thermo-stress conditions. Whole-genome scans identified at least six genomic regions on BTA1, BTA10, BTA11, BTA17, BTA21 and BTA23 associated with conception under thermal stress. These regions harbor candidate genes such as BRWD1, EXD2, ADAM20, EPAS1, TAOK3, and NOS1, which are directly implicated in reproductive functions and cellular response to heat stress. The gene-set enrichment analysis revealed functional terms related to fertilization, developmental biology, heat shock proteins and oxidative stress, among others. Overall, our findings contribute to a better understanding of the genetics underlying the reproductive performance of dairy cattle under heat stress conditions and point out novel genomic strategies for improving thermotolerance and fertility via marker-assisted breeding.     

Lack of association between genetic polymorphisms in ROBO1, MRPL19/C2ORF3 and THEM2 with Developmental Dyslexia

Developmental Dyslexia (DD) is a heritable, complex genetic disorder characterized by specific impairment in reading and writing ability that is substantially below the expected reading ability given the person's chronological age, measured intelligence and age-appropriate education. More than ten susceptible genes have been identified for DD. A Single Nucleotide Polymorphism (SNP) of these genes was found to be associated with various phenotypes of DD. To identify the role of SNPs of four candidate genes namely, MRPL19/C2ORF3, ROBO1 and THEM2 in an Indian population, we genotyped eight SNPs of these genes in 157 children with DD and 212 normal readers using a MassARRAY technique with a MALDI-TOF MS analyzer. Power analysis of some of these SNPs showed > 80% of power. Chi-square test, Odds Ratios (ORs), 95% Confidence Intervals (CIs) and Bonferroni's correction were applied to identify the significance of the genotyped SNPs and haplotypes. Our study failed to show any association of SNPs and haplotypes of these genes with DD in an Indian population.     

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.     

Analysis of dyslexia candidate genes in the Raine cohort representing the general Australian population

Several genes have been suggested as dyslexia candidates. Some of these candidate genes have been recently shown to be associated with literacy measures in sample cohorts derived from the general population. Here, we have conducted an association study in a novel sample derived from the Australian population (the Raine cohort) to further investigate the role of dyslexia candidate genes. We analysed markers, previously reported to be associated with dyslexia, located within the MRPL19/C2ORF3, KIAA0319, DCDC2 and DYX1C1 genes in a sample of 520 individuals and tested them for association with reading and spelling measures. Association signals were detected for several single nucleotide polymorphisms (SNPs) within DYX1C1 with both the reading and spelling tests. The high linkage disequilibrium (LD) we observed across the DYX1C1 gene suggests that the association signal might not be refined by further genetic mapping.     

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.     

Genome‐wide association study for pigmentation traits in Chinese Holstein population

With the Illumina BovineSNP50K BeadChip, we performed a genome‐wide association study (GWAS) for two pigmentation traits in a Chinese Holstein population: proportion of black (PB) and teat colour (TC). A case–control design was used. Cases were the cows with PB <0.30 (n = 129) and TC <2 points (n = 140); controls were those with PB >0.90 (n = 58) and TC >4 points (n = 281). The RM test of roadtrips (version 1.2) was applied to detect SNPs for the two traits with 42 883 and 42 741 SNPs respectively. A total of nine and 12 genome‐wide significant (P < 0.05) SNPs associated with PB and TC respectively were identified. Of these, two SNPs for PB were located within the KIT and IGFBP7 genes, and the other four SNPs were 23~212 kb away from the PDGFRA gene on BTA6; nine SNPs associated with TC were located within or 21~78.8 kb away from known genes on chromosomes 4, 11, 22, 23 and 24. By combing through our GWAS results and the biological functions of the genes, we suggest that the KIT, IGFBP7, PDGFRA, MITF, ING3 and WNT16 genes are promising candidates for PB and TC in Holstein cattle, providing a basis for further investigation on the genetic mechanism of pigmentation formation.     

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.     

Helicobacter pylori infection in connective tissue disorders is associated with high levels of antibodies to mycobacterial hsp65 but not to human hsp60

Background. To investigate whether the Helicobacter pylori status influences levels of antibodies against mycobacterial heat shock protein (hsp) 65 and human hsp60 in systemic autoimmune diseases and to study the concentration of anti‐H. pylori antibodies in autoimmune patients and healthy controls. Materials and Methods. Antibodies against human heat‐shock protein hsp60, mycobacterial heat‐shock protein hsp65 were analyzed by ELISA. Anti‐Helicobacter antibodies were determined by enzyme immunoassay. Results. There was a markedly higher prevalence of H. pylori infection in undifferentiated connective tissue disease (82%) (n = 33) and systemic sclerosis (78%) (n = 55) but not in systemic lupus erythematosus (n = 49), polymyositis/dermatomyositis (n = 14), rheumatoid arthritis (n = 21) or primary Raynaud's syndrome (n = 26) compared with controls (59%) (n = 349). In autoimmune diseases H. pylori infection was associated with elevated levels of antihsp65 (p = .008) but not of antihsp60. Anti‐hsp65 levels were significantly higher in H. pylori‐infected (n = 129) than in uninfected patients (n = 69) (p = .0007). Conclusions. These findings indicate that in autoimmune diseases the infection with the H. pylori bacterium is associated with increased concentration of antimycobacterial hsp65.     

Genetic risk variants for dyslexia on chromosome 18 in a German cohort

Dyslexia is characterized by impaired reading and spelling. The disorder has a prevalence of about 5% in Germany, and a strong hereditary component. Several loci are thought to be involved in the development of dyslexia. Scerri et al. identified eight potential dyslexia‐associated single nucleotide polymorphisms (SNPs) in seven genes on chromosome 18 in an English‐speaking population. Here, we present an association analysis that explores the relevance of these SNPs in a German population comprising 388 dyslexia cases and 364 control cases. In case–control analysis, three nominal SNP associations were replicated. The major alleles of NEDD4L‐rs12606138 and NEDD4L‐rs8094327 were risk associated [odds ratio (OR) = 1.35, 95% confidence interval (CI) = 1.0–1.7, P‐value = 0.017 and OR = 1.39, 95% CI = 1.1–1.7, P‐value = 0.007, respectively], and both SNPs were in strong linkage disequilibrium (r² = 0.95). For MYO5B‐rs555879, the minor allele was risk associated (OR = 1.31, 95% CI = 1.1–1.6, P‐value = 0.011). The combined analysis of SNP sets using set enrichment analysis revealed a study‐wide significant association for three SNPs with susceptibility for dyslexia. In summary, our results substantiate genetic markers in NEDD4L and MYO5B as risk factors for dyslexia and provide first evidence that the relevance of these markers is not restricted to the English language .