Genetic association between NOD2 polymorphism and infection status by Mycobacterium avium ssp. paratuberculosis in German Holstein cattle

The aim of the study was the analysis of the nucleotide‐binding oligomerization domain containing 2 (NOD2, formerly CARD15) as a candidate gene for Mycobacterium avium ssp. paratuberculosis infection in cattle. Eleven SNPs in the NOD2 gene were identified, and finally, four SNPs were included in a case–control study using 324 German Holstein cows tested for paratuberculosis using fecal culture and ELISA. The SNP (GenBank) AY518738S04:g.521G>A in exon 4 showed a significant association between the fecal culture status of the animals and NOD2 allele variants. The other three SNPs showed no associations in German Holstein cows.     

Meta-analysis of two genome-wide association studies of bovine paratuberculosis

Background

Bovine paratuberculosis (ParaTB) also known as Johne''s disease, is a contagious fatal disease resulting from infection by Mycobacterium avium subspecies paratuberculosis (MAP). Previous studies have identified loci associated with ParaTB using different measurements to define cases and controls. The objective of this study was to combine the data from two recent studies to identify genetic loci associated with MAP tissue infection and humoral immune response, defined by MAP ELISA-positive cattle, by comparing cases and control animals for one or both measures of infection.

Methodology/Principal Findings

The two populations used for the association analyses were a cohort of MAP tissue infected animals and control Holstein cows from the USA and the second cohort composed of ELISA-positive and ELISA-negative Holstein cows from Italy. Altogether 1190 cattle were genotyped with the Illumina BovineSNP50 BeadChip. SNP markers were removed if the minor allele frequency <0.01 or genotyping failure was >5%. Animals were removed with >5% genotyping failure. Whole genome association analyses were conducted with the GRAMMAR-CG method using two different definitions of control populations.

Conclusion/Significance

The analyses identified several loci (P<5 e-05) associated with ParaTB, defined by positive ELISA and presence of bacteria in tissue compared to ELISA and tissue negative animals, on chromosomes 1, 12 and 15 and one unassigned SNP. These results confirmed associations on chromosome 12 and the unassigned SNP with ParaTB which had been found in the Italian population alone. Furthermore, several additional genomic regions were found associated with ParaTB when ELISA and tissue positive animals were compared with tissue negative samples. These loci were on chromosomes 1, 6, 7, 13, 16, 21,23 and 25 (P<5 e-05). The results clearly indicate the importance of the phenotype definition when seeking to identify markers associated with different disease responses.     

Genome-wide association study of <Emphasis Type="Italic">Mycobacterium avium</Emphasis> subspecies <Emphasis Type="Italic">Paratuberculosis</Emphasis> infection in Chinese Holstein

Background

Paratuberculosis is a contagious, chronic and enteric disease in ruminants, which is caused by Mycobacterium avium subspecies paratuberculosis (MAP) infection, resulting in enormous economic losses worldwide. There is currently no effective cure for MAP infection or a vaccine, it is thus important to explore the genetic variants that contribute to host susceptibility to infection by MAP, which may provide a better understanding of the mechanisms of paratuberculosis and benefit animal genetic improvement. Herein we performed a genome-wide association study (GWAS) to identify genomic regions and candidate genes associated with susceptibility to MAP infection in dairy cattle.

Results

Using Illumina Bovine 50?K (54,609 SNPs) and GeneSeek HD (138,893 SNPs) chips, two analytical approaches were performed, GRAMMAR-GC and ROADTRIPS in 937 Chinese Holstein cows, among which individuals genotyped by the 50?K chip were imputed to HD SNPs with Beagle software. Consequently, 15 and 11 significant SNPs (P?<?5?×?10^??5) were identified with GRAMMAR-GC and ROADTDRIPS, respectively. A total of 10 functional genes were in proximity to (i.e., within 1?Mb) these SNPs, including IL4, IL5, IL13, IRF1, MyD88, PACSIN1, DEF6, TDP2, ZAP70 and CSF2. Functional enrichment analysis showed that these genes were involved in immune related pathways, such as interleukin, T cell receptor signaling pathways and inflammatory bowel disease (IBD), implying their potential associations with susceptibility to MAP infection. In addition, by examining the publicly available cattle QTLdb, a previous QTL for MAP was found to be overlapped with one of regions detected currently at 32.5?Mb on BTA23, where the TDP2 gene was anchored.

Conclusions

In conclusion, we identified 26 SNPs located on 15 chromosomes in the Chinese Holstein population using two GWAS strategies with high density SNPs. Integrated analysis of GWAS, biological functions and the reported QTL information helps to detect positional candidate genes and the identification of regions associated with susceptibility to MAP traits in dairy cattle.

     

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.     

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 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.     

Single nucleotide polymorphisms,haplotypes and combined genotypes of lactoferrin gene and their associations with mastitis in Chinese Holstein cattle

Lactoferrin (Lf) is naturally produced by the mammary gland, having biological functions of antibacterial and anti-inflammatory activities. To investigate whether the Lf gene is associated with mastitis in dairy cattle, a DNA sequencing approach was used to identify single nucleotide polymorphisms (SNPs) in the gene. Three previously reported SNPs in the 5′ flanking region and one novel SNP in exon1 of Lf gene were identified. A total of 353 individuals from Holstein cattle populations were genotyped for their SNPs using Created Restriction Site PCR (CRS-PCR) and PCR-RFLP methods. Twenty-two and nineteen combinations of three SNPs (g.3440T>G, g.3879_3880insG, and g.4432T>C) and another three SNPs (g.3429G>A, g.3440T>G, g.3879_3880insG) were observed, respectively. The result of haplotype analysis of four SNPs showed that fourteen different haplotypes were identified. Two major haplotypes (GECB and GECA) occurred with a frequency of 22.5 and 18.5% in the study population, respectively. Statistical analyses revealed no significant association between one single SNP of Lf gene and SCS, whereas significant associations between their combined genotypes of three SNPs, haplotype and SCS. Combined genotype EFCDBB and GGEFDD with the lowest SCS were favorable for the mastitis resistance. They may be used as a possible candidate for marker-assisted selection in dairy cattle breeding program.     

Association of single nucleotide polymorphisms in the ANKRA2 and CD180 genes with bovine respiratory disease and presence of Mycobacterium avium subsp. paratuberculosis(1)

The objective was to determine whether single nucleotide polymorphisms (SNPs) in the ANKRA2 and CD180 genes are associated with incidence of bovine respiratory disease (BRD) and presence of Mycobacterium avium subsp. paratuberculosis (MAP) in cattle. Two independent populations were used. The first population (BRD‐affected; N = 90) was composed of 31 half‐sib progeny, from a Brahman × Angus sire, that were treated for BRD. Untreated offspring from the sire were selected to serve as controls. The second population (MAP‐infected) of 330 animals of unknown parentage was evaluated for the presence of MAP in ileocecal lymph node and classified as positive or negative. Markers in both genes were assessed for association in these two populations. In the BRD‐affected population, five SNPs in the ANKRA2 gene were significantly associated (P < 0.05), and two SNPs were highly associated (P < 0.01) with incidence of BRD. In addition, two SNPs in the CD180 gene were found to be associated with this trait. In the MAP‐infected population, one SNP in the ANKRA2 gene was significantly associated (P < 0.05) with the presence or absence of MAP, and a SNP in the CD180 gene was highly associated (P < 0.01) with the trait. Haplotypes, using significant markers, showed a positive association with both incidence of BRD (P = 0.0001) and with the presence of MAP (P = 0.0032). Markers in the ANKRA2 and CD180 genes are associated with the ability of the animal to cope with pathogens.     

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.     

Mycoflora of cervicovaginal fluids in dairy cows with or without reproductive disorders

Fungal infection of reproduction system of dairy cattle have not been received much more attention. The aim of this study was to determine the fungal infection of Holstein dairy cows with reproductive disorders or healthy. Fungal isolates of cervicovaginal fluids of 176 Holstein dairy cows were collected by using the double rod swabs for cervix and the sterilized cotton swabs for discharges of vagina. They were evaluated for fungal infections. The treatment group included 70 dairy cows with reproductive diseases, such as abortion, repeat breeder, endometritis, metritis, retention of fetal membrane, dystocia, cervicitis, and vaginitis. The control groups were included 42 healthy non-pregnant cows and 64 pregnant cows. Isolates of fungi were obtained from cervix and vagina of 27.1% and 28.6% of treatment group, 26.7% and 31.2% of pregnant cows, and 33.3% and 21.4% of healthy non-pregnant cows, respectively, indicating no significant differences. It is showed that the cervix and vagina of the treatment group have been infected by six different mycotic isolates. However, the cervix and vagina of pregnant and non-pregnant healthy cows in control group were infected with 5, 6 and 5, 4 different fungal agents, respectively. Penicillium and yeast were the most common isolated agents. Regarding to the result of this study, it is concluded that fungal infections can occur in cervicovaginal cavity of Holstein dairy cows with or without reproductive diseases.     

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.     

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.     

Genomic breeding values,SNP effects and gene identification for disease traits in cow training sets

Holstein Friesian cow training sets were created according to disease incidences. The different datasets were used to investigate the impact of random forest (RF) and genomic BLUP (GBLUP) methodology on genomic prediction accuracies. In addition, for further verifications of some specific scenarios, single‐step genomic BLUP was applied. Disease traits included the overall trait categories of (i) claw disorders, (ii) clinical mastitis and (iii) infertility from 80 741 first lactation Holstein cows kept in 58 large‐scale herds. A subset of 6744 cows was genotyped (50K SNP panel). Response variables for all scenarios were de‐regressed proofs (DRPs) and pre‐corrected phenotypes (PCPs). Initially, all sick cows were allocated to the testing set, and healthy cows represented the training set. For the ongoing cow allocation schemes, the number of sick cows in the training set increased stepwise by moving 10% of the sick cows from the testing to the training set in each step. The size of training and testing sets was kept constant by replacing the same number of cows in the testing set with (randomly selected) healthy cows from the training set. For both the RF and GBLUP methods, prediction accuracies were larger for DRPs compared to PCPs. For PCPs as a response variable, the largest prediction accuracies were observed when the disease incidences in training sets reflected the disease incidence in the whole population. A further increase in prediction accuracies for some selected cow allocation schemes (i.e. larger prediction accuracies compared to corresponding scenarios with RF or GBLUB) was achieved via single‐step GBLUP applications. Correlations between genome‐wide association study SNP effects and RF importance criteria for single SNPs were in a moderate range, from 0.42 to 0.57, when considering SNPs from all chromosomes or from specific chromosome segments. RF identified significant SNPs close to potential positional candidate genes: GAS1, GPAT3 and CYP2R1 for clinical mastitis; SPINK5 and SLC26A2 for laminitis; and FGF12 for endometritis.     

Identification of complex vertebral malformation carriers in Holstein cattle in south China

Complex vertebral malformation (CVM) is a recently described monogenic autosomal recessive hereditary defect of Holstein dairy cattle that causes premature birth, aborted fetuses and stillborn calves. Guanine is substituted by thymine (G>T) in the solute carrier family 35 member A3 gene (SLC35A3). A valine is changed to a phenylalanine at position 180 of uridine 5'-diphosphate-N-acetyl-glucosamine transporter protein. CVM is expected to occur in many countries due to the widespread use of sire semen. We developed a created restriction site PCR (CRS-PCR) method to diagnose CVM in dairy cows. This was tested on 217 cows and 125 bulls selected randomly from a Holstein cattle population in south China. Five Holstein cows and five Holstein bulls were identified to be CVM carriers; the percentages of CVM carriers were estimated to be 2.3, 4.0 and 2.9% in the cows, bulls and entire Holstein cattle sample, respectively.     

Haplotype analysis of TLR4 gene and its effects on milk somatic cell score in Chinese commercial cattle

Bovine mastitis is a very complex and common disease of dairy cattle and a major source of economic losses to the dairy industry worldwide. In this study, the bovine TLR4 was taken as a candidate gene for mastitis resistance. This study aimed to analyze the associations of single nucleotide polymorphisms (SNP) or haplotype and somatic cell score (SCS) in 404 Chinese commercial dairy cattle including Chinese Holstein, Sanhe cattle and Chinese Simmental breeds. The polymerase chain reaction and sequencing methods were used for detecting genotype and allele frequency distribution of the two SNPs (rs8193062, rs8193064), statistical results showed that T allele at rs8193062 and C allele at rs8193064 were the predominate alleles. Moreover, six SNPs, including two SNPs (rs8193062, rs8193064) and four SNPs (rs8193060, rs8193069, rs29017188, rs8193046) which were chosen according the polymorphism level for the same cattle populations in previous studies, were used for haplotype analysis, the results revealed that twenty-one haplotypes were found in the mentioned animals, of which, Hap1 (30.5 %) and Hap2 (30.4 %) were the most common haplotypes. Hap2, Hap4 and Hap12 might negatively effect on milk SCS, whereas Hap13 might positively effect on milk SCS. The results in this study might assist in marker assisted selection and provided some reference to be implemented in breeding programs to improve the mastitis resistance of dairy cattle.     

Bovine CLEC7A genetic variants and their association with seropositivity in Johne's disease ELISA

Mycobacterium avium ssp. paratuberculosis (MAP) infection in cattle causes significant economic losses to the dairy and beef industries resulting from reduced productivity, premature culling and mortality. Bovine Dectin-1, an important pattern recognition molecule that is able to generate a proinflammatory response by acting alongside Toll like receptor (TLR) signaling, is known to co-operate with TLR2 to specifically activate a macrophage proinflammatory response against mycobacterial infections. Therefore, the goal of this study was to identify single nucleotide polymorphisms (SNPs) in the gene encoding bovine Dectin-1 (CLEC7A) and to assess their association with susceptibility to MAP infection in dairy cattle. Blood and milk samples, collected from commercial dairy operations, were tested for MAP infection using blood and milk ELISAs and a resource population consisting of 197 infected and 242 healthy cattle was constructed. Pooled DNA was used for sequencing and eight single nucleotide polymorphisms (SNPs) were identified. Identified SNPs were genotyped on the resource population using the iPLEX MassARRAY system and statistical analysis was performed using logistic regression fitting the additive and dominance effects of each SNP in the model. Out of a total of eight identified SNPs, five were successfully genotyped, and three out of these five SNPs were found to be in complete linkage. Statistical analysis revealed a strong association between a non-synonymous SNP c.589A>G (p = 0.008), and MAP infection status of the resource population inferred by seropositivity in MAP antibody specific ELISAs. This SNP c.589A>G was located in the geneic region that encodes the carbohydrate recognition domain of bovine Dectin-1. Therefore, further investigation of its functional relevance is warranted.     

A post‐GWAS confirming the SCD gene associated with milk medium‐ and long‐chain unsaturated fatty acids in Chinese Holstein population

The stearoyl‐CoA desaturase (delta‐9‐desaturase) gene encodes a key enzyme in the cellular biosynthesis of monounsaturated fatty acids. In our initial genome‐wide association study (GWAS) of Chinese Holstein cows, 19 SNPs fell in a 1.8‐Mb region (20.3–22.1 Mb) on chromosome 26 underlying the SCD gene and were highly significantly associated with C14:1 or C14 index. The aims of this study were to verify whether the SCD gene has significant genetic effects on milk fatty acid composition in dairy cattle. By resequencing the entire coding region of the bovine SCD gene, a total of six variations were identified, including three coding variations (g.10153G>A, g.10213T>C and g.10329C>T) and three intronic variations (g.6926A>G, g.8646G>A and g.16158G>C). The SNP in exon 3, g.10329C>T, was predicted to result in an amino acid replacement from alanine (GCG) to valine (GTG) in the SCD protein. An association study for 16 milk fatty acids using 346 Chinese Holstein cows with accurate phenotypes and genotypes was performed using the mixed animal model with the proc mixed procedure in sas 9.2. All six detected SNPs were revealed to be associated with six medium‐ and long‐chain unsaturated fatty acids (P = 0.0457 to P < 0.0001), specifically for C14:1 and C14 index (P = 0.0005 to P < 0.0001). Subsequently, strong linkage disequilibrium (D′ = 0.88–1.00) was observed among all six SNPs in SCD and the five SNPs (rs41623887, rs109923480, rs42090224, rs42092174 and rs42091426) within the 1.8‐Mb region identified in our previous GWAS, indicating that the significant association of the SCD gene with milk fatty acid content traits reduced the observed significant 1.8‐Mb chromosome region in GWAS. Haplotype‐based analysis revealed significant associations of the haplotypes encompassing the six SCD SNPs and one SNP (rs109923480) in a GWAS with C14:1, C14 index, C16:1 and C16 index (P = 0.0011 to P < 0.0001). In summary, our findings provide replicate evidence for our previous GWAS and demonstrate that variants in the SCD gene are significantly associated with milk fatty acid composition in dairy cattle, which provides clear evidence for an increased understanding of milk fatty acid synthesis and enhances opportunities to improve milk‐fat composition in dairy cattle.     

A genome‐wide association study suggests new candidate genes for milk production traits in Chinese Holstein cattle

A genome‐wide association study (GWAS) was conducted on 15 milk production traits in Chinese Holstein. The experimental population consisted of 445 cattle, each genotyped by the GGP (GeneSeek genomic profiling)‐BovineLD V3 SNP chip, which had 26 151 public SNPs in its manifest file. After data cleaning, 20 326 SNPs were retained for the GWAS. The phenotypes were estimated breeding values of traits, provided by a public dairy herd improvement program center that had been collected once a month for 3 years. Two statistical models, a fixed‐effect linear regression model and a mixed‐effect linear model, were used to estimate the association effects of SNPs on each of the phenotypes. Genome‐wide significant and suggestive thresholds were set at 2.46E‐06 and 4.95E‐05 respectively. The two statistical models concurrently identified two genome‐wide significant (P < 0.05) SNPs on milk production traits in this Chinese Holstein population. The positional candidate genes, which were the ones closest to these two identified SNPs, were EEF2K (eukaryotic elongation factor 2 kinase) and KLHL1 (kelch like family member 1). These two genes could serve as new candidate genes for milk yield and lactation persistence, yet their roles need to be verified in further function studies.     

Selection and use of SNP markers for animal identification and paternity analysis in U.S. beef cattle

DNA marker technology represents a promising means for determining the genetic identity and kinship of an animal. Compared with other types of DNA markers, single nucleotide polymorphisms (SNPs) are attractive because they are abundant, genetically stable, and amenable to high-throughput automated analysis. In cattle, the challenge has been to identify a minimal set of SNPs with sufficient power for use in a variety of popular breeds and crossbred populations. This report describes a set of 32 highly informative SNP markers distributed among 18 autosomes and both sex chromosomes. Informativity of these SNPs in U.S. beef cattle populations was estimated from the distribution of allele and genotype frequencies in two panels: one consisting of 96 purebred sires representing 17 popular breeds, and another with 154 purebred American Angus from six herds in four Midwestern states. Based on frequency data from these panels, the estimated probability that two randomly selected, unrelated individuals will possess identical genotypes for all 32 loci was 2.0 × 10⁻¹³ for multi-breed composite populations and 1.9 × 10⁻¹⁰ for purebred Angus populations. The probability that a randomly chosen candidate sire will be excluded from paternity was estimated to be 99.9% and 99.4% for the same respective populations. The DNA immediately surrounding the 32 target SNPs was sequenced in the 96 sires of the multi-breed panel and found to contain an additional 183 polymorphic sites. Knowledge of these additional sites, together with the 32 target SNPs, allows the design of robust, accurate genotype assays on a variety of high-throughput SNP genotyping platforms.     

Higher non-return rate associated with Mycobacterium avium subspecies paratuberculosis infection at early stage in Holstein dairy cows

The effects of infection by Mycobacterium avium paratuberculosis (Map) on dairy cows are poorly documented and quite controversial. This retrospective study aimed at quantifying the variation in non-return to service of Holstein dairy cows according to their Map-infection status. Three different statuses were defined based on both individual and herd tests results: ELISA positive cow, all tests negative cow in a negative herd and all tests negative cow in a positive herd. Whatever the age at Map testing, the status was attributed to a cow from its first lactation onwards. Non-return to service was determined at 200 days after first and second services. The study was performed from 1999 to 2007 on 185,950 AI from 48,914 cows in early stage of the infection in 1069 herds by logistic regression controlling for known factors influencing non-return rate. Non-return rate was higher for infected cows compared to negative cows from negative herds (RR of 1.10 or +3.9 points of % of non-return rate). The effect was significant for parities 1 and 2 (RR of 1.11 and 1.12, respectively) but not for higher ones. This effect was lower when comparing positive cows to negative cows in the same herds but relative risks were still above 1. The hypothesis that the effect of Map on non-return depends upon the stage of infection is formulated.