Genome-wide association analysis for susceptibility to infection by Mycobacterium avium ssp. paratuberculosis in US Holsteins.

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 and concurrently modifying management to reduce the likelihood of infection. The current study was motivated by the hypothesis that genetic variation in host susceptibility to MAP infection can be dissected and quantifiable associations with genetic markers identified. Two separate GWAS analyses were conducted, the first using 897 genotyped Holstein artificial insemination sires with phenotypes derived from incidence of MAP infection among daughters based on milk ELISA testing records. The second GWAS analysis was a case-control design using US Holstein cows phenotyped for MAP infection by serum ELISA or fecal culture tests. Cases included cows positive for either serum ELISA, fecal culture, or both. Controls consisted of animals negative for all tests conducted. A total of 376 samples (70 cases and 306 controls) from a University of Minnesota Johne's management demonstration project and 184 samples (76 cases and 108 controls) from a Michigan State University study were used. Medium-density (sires) and high-density (cows) genotype data were imputed to full genome sequence for the analyses. Marker-trait associations were analyzed using the single-step (ss)GWAS procedure implemented in the BLUPF90 suite of programs. Evidence of significant genomic contributions for susceptibility to MAP infection were observed on multiple chromosomes. Results were combined across studies in a meta-analysis, and increased support for genomic regions on BTA7 and BTA21 were observed. Gene set enrichment analysis suggested pathways for antigen processing and presentation, antimicrobial peptides and natural killer cell-mediated cytotoxicity are relevant to variation in host susceptibility to MAP infection, among others. Genomic prediction was evaluated using a 5-fold cross-validation, and moderate correlations were observed between genomic breeding value predictions and daughter averages (∼0.43 to 0.53) for MAP infection in testing data sets. These results suggest that genomic selection against susceptibility to MAP infection is feasible in Holstein cattle.

Herd management practices associated with bovine leukemia virus incidence rate in Michigan dairy farms.

The objective of this study was to identify associations between herd management practices and the incidence rate of bovine leukemia virus (BLV) infections in Michigan dairy herds. Previous management risk factor studies were of antibody prevalence rather than the rate of recent infections. Milk samples were collected from cohorts of cows on 112 Michigan dairy herds and tested for BLV using an antibody capture ELISA (n = 3849 cows). Cows were subsequently followed for an average of 21 months. Cows negative for anti-BLV antibodies and still present in their respective herds were retested by the same antibody capture ELISA to estimate within-herd incidence rates. The overall crude incidence rate was 1.46 infections per 100 cow-months at risk for the 1314 retested cows in 107 herds. The average within-herd incidence rate was 2.28 infections per 100 cow-months (range: 0 to 9.76 infections per 100 cow-months). A negative binomial regression model was used to identify herd management practices associated with the within-herd incidence rate. Results of the final multivariable model identified higher herd prevalence, milking frequency, needle reuse, as well as housing post-parturient cows separately, to be associated with increased incidence rate. Utilization of sand bedding for the lactating herd was found to be associated with decreased incidence rates. Results of this study suggest potential routes of BLV transmission which should be further investigated as disease control targets in ongoing control programs.

Controlling bovine leukemia virus in dairy herds by identifying and removing cows with the highest proviral load and lymphocyte counts.

The objective of this field trial was to reduce bovine leukemia virus (BLV) transmission and prevalence in commercial dairy herds using proviral load (PVL) and lymphocyte count (LC) measurements as indicators of the most infectious animals for culling or segregation. Bovine leukemia virus causes lymphoma in <5% of infected cattle, and increased lymphocyte counts (lymphocytosis) in about one-third. Recent research has shown that dairy cows infected with BLV have altered immune function associated with decreases in milk production and lifespan. Recent findings show that a minority of infected cattle have PVL concentrations in blood and other body fluids of over 1,000 times that of other infected cattle. In combination with a high LC, these animals are thought to be responsible for most transmission of BLV in a herd. Milk or blood samples from adult cows in our 3 Midwestern dairy farm field trials were tested semiannually with ELISA for BLV antibodies, and ELISA-positive cattle were then retested using a blood LC and a quantitative PCR test for PVL to identify the animals presumed to be most infectious. Herd managers were encouraged to consider PVL and LC status when making cull decisions, and to segregate cows with the highest PVL and LC from their BLV ELISA-negative herd mates where possible. After 2 to 2.5 yr of this intervention, the incidence risk of new infections decreased in all 3 herds combined, from 13.8 to 2.2, and the overall herd prevalence decreased in all 3 herds combined from 62.0 to 20.7%, suggesting that this approach can efficiently reduce BLV transmission as well as prevalence. This is encouraging, because a very low prevalence of BLV infection would make it economically feasible to cull the remaining ELISA-positive cattle, as was achieved in national eradication programs in other countries decades ago.

Differential gene expression between normal and pale, soft, and exudative turkey meat.

In response to high consumer demand, turkeys have been intensively selected for rapid growth rate and breast muscle mass and conformation. The success in breeding selection has coincided with an increasing incidence of pale, soft, and exudative (PSE) meat defect, especially in response to heat stress. We hypothesized that the underlying mechanism responsible for the development of PSE meat arises from differences in expression of several critical genes. The objective of this study was to determine differential gene expression between normal and PSE turkey meat using a 6K turkey skeletal muscle long oligonucleotide microarray. Breast meat samples were collected from Randombred Control Line 2 turkeys at 22 wk of age, and classified as normal or PSE primarily based on marinade uptake (high = normal, low = PSE). Total RNA was isolated from meat samples with the highest (normal, n = 6) and the lowest (PSE, n = 6) marinade uptake. Microarray data confirmation was conducted using quantitative real-time PCR. Selection of differentially expressed genes for pathway analysis was performed using a combination of fold change (FC) ranking (FC < -1.66, FC >1.66) and false discovery rate (<0.35) as criteria. The calcium signaling pathway was highlighted as the top canonical pathway associated with differential gene expression between normal and PSE turkey. Dramatic downregulation of fast-twitch myosin heavy chain coupled with upregulation of slow-twitch myosin and troponin C suggested a switch of skeletal muscle isoforms, which may alter muscle fiber arrangement and formation of actin-myosin complexes. Changes in expression of genes in the actin cytoskeleton signaling pathway also suggest altered structures of actin filaments that may affect cell motility as well as strength and flexibility of muscle cells. Substantial downregulation of pyruvate dehydrogenase kinase, isozyme 4 was observed in PSE samples, suggesting altered regulation of the aerobic metabolic pathway in the birds that developed PSE meat defect.

Versican, matrix Gla protein, and death-associated protein expression affect muscle satellite cell proliferation and differentiation.

Our previous transcriptional profiling study using a turkey skeletal muscle-specific oligonucleotide microarray revealed over 3,000 genes that were differentially expressed at 3 critical stages of muscle development: 18 d embryonic, 1 d posthatch, and 16 wk of age. The genes versican, matrix Gla protein (MGP), and death-associated protein (DAP) were selected to study for their potential effects on muscle satellite cell proliferation and differentiation, as their functions in other tissues are suggestive of possible key roles in the regulation of myogenesis and they are differentially expressed throughout muscle development in the turkey. Using small interfering RNA to knockdown the expression of these genes during proliferation and differentiation, each of the genes was found to differentially affect proliferation and differentiation. Versican and MGP predominantly affected proliferation with line effects, but later stages of differentiation were affected by the knockdown of versican and MGP. The underexpression of DAP inhibited myotube formation, which is a necessary stage in the development of muscle fibers. Without myotube development, muscle fiber formation will be inhibited or abolished. This is the first report that these genes with no previously documented functions with regard to muscle development play a critical role in muscle cell proliferation and differentiation.

Differential expression of calcium-regulating genes in heat-stressed turkey breast muscle is associated with meat quality.

Aberrant postmortem Ca(2+)-regulation in the early postmortem period is associated with the occurrence of inferior meat quality in turkeys, described as pale, soft, and exudative (PSE). The objective of the current study was to quantify expression of 4 candidate genes responsible for maintaining Ca(2+) homeostasis in turkey skeletal muscle as a function of heat stress: α and ß ryanodine receptors (RYR; Ca(2+)-release channels), the sarco/endoplasmic reticulum Ca(2+)-ATPase 1 (SERCA1), and the sarcoplasmic reticulum, Ca(2+)-storage protein calsequestrin (CASQ1). Two genetic lines of turkeys were used: a growth-selected commercial line and a randombred control line. Market-age birds were subjected to one of 5 heat stress treatments: no heat, 1 d, 3 d, 5 d, or 7 d of heat followed by 7 d of ambient temperature. Breast muscle samples were harvested and classified as normal or PSE using the meat quality parameters percentage of marinade uptake and percentage of cook loss. These parameters differed significantly by line, heat stress treatment, and meat quality status. Expression of candidate genes was measured using TaqMan quantitative PCR. Heat treatment was associated with significantly enhanced expression of αRYR, ßRYR, and CASQ1 in normal muscle from both lines. Conversely, mRNA abundance of these genes was reduced in PSE muscle from both lines and recovered or increased by 7 d + 7 d of rest. Genetic line differences were observed at several time points. Expression of SERCA1 in both normal and PSE samples from both lines was unchanged or trended downward with heat stress. Taken together, genetic line and heat-stress treatment affected the expression of important Ca(2+)-regulating genes in association with meat quality status. The data suggest that birds whose meat leads to PSE may fail to respond to heat stress appropriately due to a delay in the upregulation of the important calcium-regulating genes: αRYR, ßRYR, and CASQ1.

Differential expression of genes characterizing myofibre phenotype.

Skeletal muscle is composed of metabolically heterogeneous myofibres that exhibit high plasticity at both the morphological and transcriptional levels. The objective of this study was to employ microarray analysis to elucidate the differential gene expression between the tonic-'red' anterior latissimus dorsi (ALD) muscle, the phasic-'white' posterior latissimus dorsi (PLD) and 'mixed'-phenotype biceps femoris (BF) in 1-week-and 19-week-old male turkeys. A total of 170 differentially expressed genes were identified in the muscle samples analysed (P < 0.05). Gene GO analysis software was utilized to identify top gene networks and metabolic pathways involving differentially expressed genes. Quantitative real-time PCR for selected genes (BAT2D, CLU, EGFR and LEPROT) was utilized to validate the microarray data. The largest differences were observed between ALD and PLD muscles, in which 32 genes were over-expressed and 82 genes were under-expressed in ALD1-PLD1 comparison, and 70 genes were over-expressed and 70 under-expressed in ALD19-PLD19 comparison. The largest number of genes over-expressed in ALD muscles, as compared to other muscles, code for extracellular matrix proteins such as dystroglycan and collagen. The gene analysis revealed that phenotypically 'red' BF muscle has high expression of glycolytic genes usually associated with the 'white' muscle phenotype. Muscle-specific differences were observed in expression levels of genes coding for proteins involved in mRNA processing and translation regulation, proteosomal degradation, apoptosis and insulin resistance. The current findings may have large implications in muscle-type-related disorders and improvement of muscle quality in agricultural species.

Characterization of a 6K oligonucleotide turkey skeletal muscle microarray.

Consumer demand for lean, inexpensive meat products has driven the domestic turkey (Meleagris gallopavo) industry to unprecedented production; however, this has coincided with an increase in growth-induced myopathies and meat quality defects. With the aim of developing a new tool for the study of turkey growth and development at the muscle transcriptome level, a 6K oligonucleotide microarray was constructed, the Turkey Skeletal Muscle Long Oligo (TSKMLO) microarray. Skeletal muscle samples were collected at three critical stages in muscle development: 18-day embryo (hyperplasia), 1-day post-hatch (hypertrophy), and 16-week (market age) from two genetic lines of turkeys: RBC2, a line maintained without selection pressure, and F, a line selected from the RBC2 line for increased 16-week body weight. Oligonucleotides were designed from sequences obtained from skeletal muscle cDNA libraries from the three developmental stages. Several unique controls, including mismatch and distance controls and scrambled sequences, were designed for 30 genes. Quality control hybridizations were completed, confirming the validity and repeatability of the array. Control features were evaluated across two larger experiments comparing developmental stage within genetic line or genetic line within each developmental stage, totaling 70 arrays. Mismatch and scrambled sequences appeared to be useful controls of specific hybridization for most genes. In addition, quantitative real-time RT-PCR confirmed microarray results. This creation and assessment of the TSKMLO array provides a valuable community resource for the study of gene expression changes related to turkey muscle growth and development.