Search

Bovine leukemia virus (BLV) is the etiological agent of enzootic bovine leukosis and causes a persistent infection that can leave cattle with no symptoms. Many countries have been able to successfully eradicate BLV through improved detection and management methods. However, with the increasing novel molecular detection methods there have been few efforts to standardize these results at global scale. This study aimed to determine the interlaboratory accuracy and agreement of 11 molecular tests in detecting BLV. Each qPCR/ddPCR method varied by target gene, primer design, DNA input and chemistries. DNA samples were extracted from blood of BLV-seropositive cattle and lyophilized to grant a better preservation during shipping to all participants around the globe. Twenty nine out of 44 samples were correctly identified by the 11 labs and all methods exhibited a diagnostic sensitivity between 74 and 100%. Agreement amongst different assays was linked to BLV copy numbers present in samples and the characteristics of each assay (i.e., BLV target sequence). Finally, the mean correlation value for all assays was within the range of strong correlation. This study highlights the importance of continuous need for standardization and harmonization amongst assays and the different participants. The results underscore the need of an international calibrator to estimate the efficiency (standard curve) of the different assays and improve quantitation accuracy. Additionally, this will inform future participants about the variability associated with emerging chemistries, methods, and technologies used to study BLV. Altogether, by improving tests performance worldwide it will positively aid in the eradication efforts., (© 2024. The Author(s).)

Bovine leukemia virus (BLV) is the causative agent of enzootic bovine leucosis (EBL), which has been reported worldwide. The expression of viral structural proteins: surface glycoprotein (gp51) and three core proteins - p15 (matrix), p24 (capsid), and p12 (nucleocapsid) induce a strong humoral and cellular immune response at first step of infection. CD4+ T-cell activation is generally induced by bovine leukocyte antigen (BoLA) region- positive antigen-presenting cells (APC) after processing of an exogenous viral antigen. Limited data are available on the BLV epitopes from the core proteins recognized by CD4+ T-cells. Thus, immunoinformatic analysis of Gag sequences obtained from 125 BLV isolates from Poland, Canada, Pakistan, Kazakhstan, Moldova and United States was performed to identify the presence of BoLA-DRB3 restricted CD4+ T-cell epitopes. The 379 15-mer overlapping peptides spanning the entire Gag sequence were run in BoLA-DRB3 allele-binding regions using a BoLA-DRB- peptide binding affinity prediction algorithm. The analysis identified 22 CD4+ T-cell peptide epitopes of variable length ranging from 17 to 22 amino acids. The predicted epitopes interacted with 73 different BoLA-DRB3 alleles found in BLV-infected cattle. Importantly, two epitopes were found to be linked with high proviral load in PBMC. A majority of dominant and subdominant epitopes showed high conservation across different viral strains, and therefore could be attractive targets for vaccine development., (© 2023. The Author(s).)

Small non-coding RNAs, such as microRNAs (miRNA) and tRNA-derived fragments (tRF), are known to be involved in post-transcriptional gene regulation. Research has provided evidence that small RNAs may influence immune development in calves. Bovine leukosis is a disease in cattle caused by Bovine Leukemia Virus (BLV) that leads to increased susceptibility to opportunistic pathogens. No research has addressed the potential influence that a maternal BLV infection may have on gene regulation through the differential expression of miRNAs or tRFs in progeny. Blood samples from 14-day old Holstein calves born to BLV-infected dams were collected. Antibodies for BLV were assessed using ELISA and levels of BLV provirus were assessed using qPCR. Total RNA was extracted from whole blood samples for small RNA sequencing. Five miRNAs (bta-miR-1, bta-miR-206, bta-miR-133a, bta-miR-133b, and bta-miR-2450d) and five tRFs (tRF-36-8JZ8RN58X2NF79E, tRF-20-0PF05B2I, tRF-27-W4R951KHZKK, tRF-22-S3M8309NF, and tRF-26-M87SFR2W9J0) were dysregulated in calves born to BLV-infected dams. The miRNAs appear to be involved in the gene regulation of immunological responses and muscle development. The tRF subtypes and parental tRNA profiles in calves born to infected dams appear to be consistent with previous publications in adult cattle with BLV infection. These findings offer insight into how maternal BLV infection status may impact the development of offspring.

Objective: Determine the association between bovine leukemia virus (BLV) status and fertility in beef cows. BLV-status was defined using 3 different testing strategies (ELISA-, quantitative polymerase chain reaction- [qPCR], and high proviral load [PVL]-status). Fertility was defined as the overall probability of pregnancy as well as the probability of becoming pregnant in the first 21 days of the breeding season., Animals: Convenience sample of 2,820 cows from 43 beef herds., Procedures: The association of BLV-status with the probability of becoming pregnant was evaluated with a multivariable logistic regression analysis that used pregnancy status as a binary outcome, herd nested within ranch as a random effect, and BLV-status (ELISA-, qPCR-, and PVL-status as separate models) and potential covariates (eg, age, Body Condition Score [BCS] category, and interactions) as fixed effects., Results: Raw data revealed that 55% (1,552/2,820) of cows were classified as BLV-positive by ELISA, and 95.3% (41/43) of herds had a least 1 ELISA-positive cow. Classification as BLV ELISA-positive was positively associated with the probability of being pregnant; however, when qPCR or PVL were used to classify BLV-status, there was no association with the probability of being pregnant. None of the methods of classifying BLV-status were associated with the probability of becoming pregnant in the first 21 days of the breeding season., Clinical Relevance: This study did not find evidence that testing beef cows for BLV-status using ELISA, qPCR, or a cut-off of 0.9 PVL and removing test-positive cows will improve cowherd fertility as described by the probability of becoming pregnant during the breeding season or becoming pregnant during the first 21 days of the breeding season.

Objective: Determine bovine leukemia virus (BLV) seroprevalence of adult female cattle in Eastern Kansas beef herds and the proviral load (PVL) of those cattle found to be ELISA positive., Animals: Convenience sample of 2,845 cows from 44 beef herds., Procedures: BLV serostatus was determined using an ELISA antibody test (gp-51; IDEXX). BLV quantitative PCR (qPCR) status and PVL were determined utilizing a qPCR test (SS1 qPCR test; CentralStar Laboratories). The association of age, herd size, and body condition score (BCS) category on the probability of being BLV positive was evaluated with a multiple variable logistic regression analysis that used BLV status as a binary outcome, herd nested within ranch as a random effect, and BCS, herd size, and age category as fixed effects., Results: Forty-two of 44 herds had at least 1 BLV ELISA-positive cow (95.5% herd seroprevalence). Overall, 1,564 of the 2,845 cows were BLV ELISA positive (55.0% individual animal prevalence). No association between BLV ELISA status and herd size or BCS was identified. When evaluated by age, the model-adjusted probability of being BLV ELISA positive was lowest for heifers (1 year of age, first parity) and increased until 5 to 6 years of age. Of the 1,564 ELISA-positive animals, 838 were qPCR positive (53.6%). The model-adjusted probability of being qPCR positive was not associated with age, herd size, or BCS category., Clinical Relevance: This study indicated that BLV-seropositive status both as a herd classification and individual animal classification was very common in this population. Because the percentage of BLV-seropositive cows varied between herds and by age, this study provides evidence that it is essential for investigators to control for herd and age in any analysis of the association of BLV serostatus and health and production outcomes of interest. Some BLV ELSIA-seropositive cows were classified as BLV negative by qPCR, and risk factors may differ between classification status by ELISA and qPCR.

Bovine leukemia virus (BLV) causes enzootic bovine leukosis, the most common neoplastic disease in cattle. Previous work estimates that 78% of US beef operations and 38% of US beef cattle are seropositive for BLV. Infection by BLV in a herd is an economic concern for producers as evidence suggests that it causes an increase in cost and a subsequent decrease in profit to producers. Studies investigating BLV in dairy cattle have noted disease resistance or susceptibility, measured by a proviral load (PVL) associated with specific alleles of the bovine leukocyte antigen (BoLA) DRB3 gene. This study aims to investigate the associations between BoLA DRB3 alleles and BLV PVL in beef cattle. Samples were collected from 157 Midwest beef cows. BoLA DRB3 alleles were identified and compared with BLV PVL. One BoLA DRB3 allele, *026:01 , was found to be associated with high PVL in relation to the average of the sampled population. In contrast, two alleles, *033:01 and *002:01 , were found to be associated with low PVL. This study provides evidence of a relationship between BoLA DRB3 alleles and BLV PVL in US beef cows.

Characterization of the bovine leukocyte antigen (BoLA) DRB3 gene has shown that specific alleles associate with susceptibility or resilience to the progression of bovine leukemia virus (BLV), measured by proviral load (PVL). Through surveillance of multi-farm BLV eradication field trials, we observed differential phenotypes within seropositive cows that persist from months to years. We sought to develop a multiplex next-generation sequencing workflow (NGS-SBT) capable of genotyping 384 samples per run to assess the relationship between BLV phenotype and two BoLA genes. We utilized longitudinal results from milk ELISA screening and subsequent blood collections on seropositive cows for PVL determination using a novel BLV proviral load multiplex qPCR assay to phenotype the cows. Repeated diagnostic observations defined two distinct phenotypes in our study population, ELISA-positive cows that do not harbor detectable levels of provirus and those who do have persistent proviral loads. In total, 565 cows from nine Midwest dairy farms were selected for NGS-SBT, with 558 cows: 168 BLV susceptible (ELISA-positive/PVL-positive) and 390 BLV resilient (ELISA-positive/PVL-negative) successfully genotyped. Three BoLA-DRB3 alleles, including one novel allele, were shown to associate with disease resilience, *009:02 , *044:01 , and *048:02 were found at rates of 97.5%, 86.5%, and 90.3%, respectively, within the phenotypically resilient population. Alternatively, DRB3*015:01 and *027:03 , both known to associate with disease progression, were found at rates of 81.1% and 92.3%, respectively, within the susceptible population. This study helps solidify the immunogenetic relationship between BoLA-DRB3 alleles and BLV infection status of these two phenotypic groupings of US dairy cattle.

Bovine leukemia virus (BLV) causes Enzootic Bovine Leukosis (EBL), a persistent life-long disease resulting in immune dysfunction and shortened lifespan in infected cattle, severely impacting the profitability of the US dairy industry. Our group has found that 94% of dairy farms in the United States are infected with BLV with an average in-herd prevalence of 46%. This is partly due to the lack of clinical presentation during the early stages of primary infection and the elusive nature of BLV transmission. This study sought to validate a near-complete genomic sequencing approach for reliability and accuracy before determining its efficacy in characterizing the sequence identity of BLV proviral genomes collected from a pilot study made up of 14 animals from one commercial dairy herd. These BLV-infected animals were comprised of seven adult dam/daughter pairs that tested positive by ELISA and qPCR. The results demonstrate sequence identity or divergence of the BLV genome from the same samples tested in two independent laboratories, suggesting both vertical and horizontal transmission in this dairy herd. This study supports the use of Oxford Nanopore sequencing for the identification of viral SNPs that can be used for retrospective genetic contact tracing of BLV transmission.

Enzootic Bovine Leukosis (EBL) caused by the bovine leukemia virus (BLV) has been eradicated in over 20 countries. In contrast, the U.S. and many other nations are experiencing increasing prevalence in the absence of efforts to control transmission. Recent studies have shown that BLV infection in dairy cattle has a greater impact beyond the long-recognized lymphoma development that occurs in <5% of infected cattle. Like other retroviruses, BLV appears to cause multiple immune system disruptions, affecting both cellular and humoral immunity, which are likely responsible for increasingly documented associations with decreased dairy production and decreased productive lifespan. Realization of these economic losses has increased interest in controlling BLV using technology that was unavailable decades ago, when many nations eradicated BLV via traditional antibody testing and slaughter methods. This traditional control is not economically feasible for many nations where the average herd antibody prevalence is rapidly approaching 50%. The ELISA screening of cattle with follow-up testing via qPCR for proviral load helps prioritize the most infectious cattle for segregation or culling. The efficacy of this approach has been demonstrated in at least four herds. Breeding cattle for resistance to BLV disease progression also appears to hold promise, and several laboratories are working on BLV vaccines. There are many research priorities for a wide variety of disciplines, especially including the need to investigate the reports linking BLV and human breast cancer.

Acute infections of bovine viral diarrhea virus (BVDV) lead to a range of clinical presentations. Laboratory tests for detection depend on collection of samples during a short viremia. Acutely infected animals remain largely undiagnosed. Transfer RNA halves (tsRNAs) are hypothesized to function like microRNAs to regulate gene expression during an immune response. The objective of this study was to identify tsRNAs in cattle that had been challenged with a non-cytopathic field strain of BVDV. Colostrum-deprived neonatal Holstein calves were either challenged with BVDV (n=5) or mock challenged (n=4). Sera was collected prior to challenge and days 4, 9, and 16 post challenge. RNA was extracted and read counts of small non-coding RNAs were assessed using next-generation sequencing. A total of 87,838,207 reads identified 41 different tsRNAs. Two 5' tsRNAs, tsRNAProAGG and tsRNAValAAC, differed across time. Two 5' tsRNAs, tsRNAGlyCCC and tsRNAGlyGCC, differed between treatment groups across time. Four days post challenge, 5' tsRNAGlyCCC and tsRNAGlyGCC were significantly lower in the challenged group than the control group. Further studies are needed to identify the importance and function of 5' tsRNAGlyCCC and tsRNAGlyGCC in serum samples of cattle challenged with BVDV.

Bovine leukemia virus (BLV) affects cattle health and productivity worldwide, causing abnormal immune function and immunosuppression. Transfer RNA fragments (tRFs) are known to be involved in inhibition of gene expression and have been associated with stress and immune response, tumor growth, and viral infection. The objective of this study was to identify tRFs associated with antibody response to BLV in Holstein cattle. Sera from 14 animals were collected to establish IgG reactivity to BLV by ELISA. Seven animals were seropositive (positive group) and seven were seronegative (negative group) for BLV exposure. Leukocytes from each animal were collected and tRFs were extracted for sequencing. tRF5 GlnCTG , tRF5 GlnTTG , and tRF5 HisGTG , were significantly different between seropositive and seronegative groups ( P < 0.0067). In all cases the positive group had a lower number of normalized sequences for tRFs when compared to the negative group. Result suggests that tRF5s could potentially be used as biomarkers to establish exposure of cattle to BLV.

Bovine viral diarrhea virus (BVDV) is an RNA virus that is often associated with respiratory disease in cattle. MicroRNAs have been proposed as indicators of exposure to respiratory pathogens. The objective of this study was to identify microRNAs in cattle that had been challenged with a non-cytopathic field strain of BVDV. Five colostrum deprived neonate Holstein calves were inoculated with BVDV (challenged) and 4 were mock challenged (control). Serum from all calves was collected at four different times: prior to challenge (day 0) and at 4, 9, and 16 days post-challenge. RNA was extracted from sera, and expression, via read counts, of small non-coding RNAs were obtained using next-generation sequencing. A total of 905,861 sequences identified 427 microRNAs. Sixty-two microRNAs had >1,000 total reads across all samples. Bta-miR-339a, bta-miR-185, bta-miR-486, Bta-miR-92a, bta-miR-30e-5p, bta-let-7c, and bta-miR-2284x were significantly different ( P < 0.05) across time regardless of challenge status. Bta-miR-423-5p ( P = 0.008) and bta-miR-151-3p ( P = 0.005) were significantly different between challenged and control animals across time. In challenged animals, bta-miR-423-5p peaked in number of reads by day 4 and steadily declined from day 4 to day 16. In control animals, bta-miR-423-5p declined from day 0 to day 9 and increased in number by day 16. By day 16, both challenged and control animals had similar levels of bta-miR-423-5p, and these levels were similar to day 0 levels. Bta-miR-151-3p peaked at day 9 in challenged animals, while control animals decreased across time. By day 16, the number of reads of bta-miR-151-3p were similar between challenged and control animals. The level in challenged animals had returned to day 0 levels by day 16, whereas the levels for control animals was significantly lower ( P = 0.006) than day 0. Further studies are needed to establish if bta-miR-423-5p or bta-miR-151-3p could be used as a biomarker for exposure to BVDV.

Background: Grazing mammals rely on their ruminal microbial symbionts to convert plant structural biomass into metabolites they can assimilate. To explore how this complex metabolic system adapts to the host animal's diet, we inferred a microbiome-level metabolic network from shotgun metagenomic data., Results: Using comparative genomics, we then linked this microbial network to that of the host animal using a set of interface metabolites likely to be transferred to the host. When the host sheep were fed a grain-based diet, the induced microbial metabolic network showed several critical differences from those seen on the evolved forage-based diet. Grain-based (e.g., concentrate) diets tend to be dominated by a smaller set of reactions that employ metabolites that are nearer in network space to the host's metabolism. In addition, these reactions are more central in the network and employ substrates with shorter carbon backbones. Despite this apparent lower complexity, the concentrate-associated metabolic networks are actually more dissimilar from each other than are those of forage-fed animals. Because both groups of animals were initially fed on a forage diet, we propose that the diet switch drove the appearance of a number of different microbial networks, including a degenerate network characterized by an inefficient use of dietary nutrients. We used network simulations to show that such disparate networks are not an unexpected result of a diet shift., Conclusion: We argue that network approaches, particularly those that link the microbial network with that of the host, illuminate aspects of the structure of the microbiome not seen from a strictly taxonomic perspective. In particular, different diets induce predictable and significant differences in the enzymes used by the microbiome. Nonetheless, there are clearly a number of microbiomes of differing structure that show similar functional properties. Changes such as a diet shift uncover more of this type of diversity.

By mapping translated metagenomic reads to a microbial metabolic network, we show that ruminal ecosystems that are rather dissimilar in their taxonomy can be considerably more similar at the metabolic network level. Using a new network bi-partition approach for linking the microbial network to a bovine metabolic network, we observe that these ruminal metabolic networks exhibit properties consistent with distinct metabolic communities producing similar outputs from common inputs. For instance, the closer in network space that a microbial reaction is to a reaction found in the host, the lower will be the variability of its enzyme copy number across hosts. Similarly, these microbial enzymes that are nearby to host nodes are also higher in copy number than are more distant enzymes. Collectively, these results demonstrate a widely expected pattern that, to our knowledge, has not been explicitly demonstrated in microbial communities: namely that there can exist different community metabolic networks that have the same metabolic inputs and outputs but differ in their internal structure., (© The Author(s) 2015. Published by Oxford University Press on behalf of Nucleic Acids Research.)

The callipyge phenotype in sheep involves substantial postnatal muscle hypertrophy and other changes to carcass composition. A single nucleotide polymorphism in the DLK1-DIO3 imprinted gene cluster alters gene expression of the paternal allele-specific protein-coding genes and several maternal allele-specific long noncoding RNA and microRNA when the mutation is inherited in cis. The inheritance pattern of the callipyge phenotype is polar overdominant because muscle hypertrophy only occurs in heterozygous animals that inherit a normal maternal allele and the callipyge SNP on the paternal allele (+/C). We examined the changes of gene expression of four major transcripts from the DLK1-DIO3 cluster and four myosin isoforms during the development of muscle hypertrophy in the semimembranosus as well as in the supraspinatus that does not undergo hypertrophy. The homozygous (C/C) animals had an intermediate gene expression pattern for the paternal allele-specific genes and two myosin isoforms, indicating a biological activity that was insufficient to change muscle mass. Transcriptome analysis was conducted by RNA sequencing in the four callipyge genotypes. The data show that homozygous animals (C/C) have lower levels of gene expression at many loci relative to the other three genotypes. A number of the downregulated genes are putative targets of the maternal allele-specific microRNA with gene ontology, indicating regulatory and cell signaling functions. These results suggest that the trans-effect of the maternal noncoding RNA and associated miRNA is to stabilize the expression of a number of regulatory genes at a functional, but low level to make the myofibers of homozygous (C/C) lambs less responsive to hypertrophic stimuli of the paternal allele-specific genes., (© 2014 Stichting International Foundation for Animal Genetics.)

We performed a genome-wide association study for Warner-Bratzler shear force (WBSF), a measure of meat tenderness, by genotyping 3360 animals from five breeds with 54 790 BovineSNP50 and 96 putative single-nucleotide polymorphisms (SNPs) within μ-calpain [HUGO nomenclature calpain 1, (mu/I) large subunit; CAPN1] and calpastatin (CAST). Within- and across-breed analyses estimated SNP allele substitution effects (ASEs) by genomic best linear unbiased prediction (GBLUP) and variance components by restricted maximum likelihood under an animal model incorporating a genomic relationship matrix. GBLUP estimates of ASEs from the across-breed analysis were moderately correlated (0.31-0.66) with those from the individual within-breed analyses, indicating that prediction equations for molecular estimates of breeding value developed from across-breed analyses should be effective for genomic selection within breeds. We identified 79 genomic regions associated with WBSF in at least three breeds, but only eight were detected in all five breeds, suggesting that the within-breed analyses were underpowered, that different quantitative trait loci (QTL) underlie variation between breeds or that the BovineSNP50 SNP density is insufficient to detect common QTL among breeds. In the across-breed analysis, CAPN1 was followed by CAST as the most strongly associated WBSF QTL genome-wide, and associations with both were detected in all five breeds. We show that none of the four commercialized CAST and CAPN1 SNP diagnostics are causal for associations with WBSF, and we putatively fine-map the CAPN1 causal mutation to a 4581-bp region. We estimate that variation in CAST and CAPN1 explains 1.02 and 1.85% of the phenotypic variation in WBSF respectively., (© 2012 The Authors, Animal Genetics © 2012 Stichting International Foundation for Animal Genetics.)

Background: Genomic selection is a recently developed technology that is beginning to revolutionize animal breeding. The objective of this study was to estimate marker effects to derive prediction equations for direct genomic values for 16 routinely recorded traits of American Angus beef cattle and quantify corresponding accuracies of prediction., Methods: Deregressed estimated breeding values were used as observations in a weighted analysis to derive direct genomic values for 3570 sires genotyped using the Illumina BovineSNP50 BeadChip. These bulls were clustered into five groups using K-means clustering on pedigree estimates of additive genetic relationships between animals, with the aim of increasing within-group and decreasing between-group relationships. All five combinations of four groups were used for model training, with cross-validation performed in the group not used in training. Bivariate animal models were used for each trait to estimate the genetic correlation between deregressed estimated breeding values and direct genomic values., Results: Accuracies of direct genomic values ranged from 0.22 to 0.69 for the studied traits, with an average of 0.44. Predictions were more accurate when animals within the validation group were more closely related to animals in the training set. When training and validation sets were formed by random allocation, the accuracies of direct genomic values ranged from 0.38 to 0.85, with an average of 0.65, reflecting the greater relationship between animals in training and validation. The accuracies of direct genomic values obtained from training on older animals and validating in younger animals were intermediate to the accuracies obtained from K-means clustering and random clustering for most traits. The genetic correlation between deregressed estimated breeding values and direct genomic values ranged from 0.15 to 0.80 for the traits studied., Conclusions: These results suggest that genomic estimates of genetic merit can be produced in beef cattle at a young age but the recurrent inclusion of genotyped sires in retraining analyses will be necessary to routinely produce for the industry the direct genomic values with the highest accuracy.

Next generation sequencing platforms have democratized genome sequencing. Large genome centers are no longer required to produce genome sequences costing millions. A few lanes of paired-end sequence on an Illumina Genome Analyzer, costing < $10,000, will produce more sequence than generated only a few years ago to produce the human and cow assemblies. The de novo assembly of large numbers of short reads into a high-quality whole-genome sequence is now technically feasible and will allow the whole genome sequencing and assembly of a broad spectrum of ruminant species. Next-generation sequencing instruments are also proving very useful for transcriptome or resequencing projects in which the entire RNA population produced by a tissue, or the entire genomes of individual animals are sequenced, and the produced reads are aligned to a reference assembly. We have used this strategy to examine gene expression differences in tissues from cattle differing in feed efficiency, to perform genome-wide single nucleotide polymorphism discovery for the construction of ultrahigh-density genotyping assays, and in combination with genome-wide association analysis, for the identification of mutations responsible for Mendelian diseases. The new 800K SNP bovine genotyping assays possess the resolution to map trait associations to the locations of individual genes and the 45 million polymorphisms identified in > 180X genome sequence coverage on over 200 animals can be queried to identify the polymorphisms present within positional candidate genes. These new tools should rapidly allow the identification of genes and mutations underlying variation in cattle production and reproductive traits.

Callipyge sheep exhibit extreme postnatal muscle hypertrophy in the loin and hindquarters as a result of a single nucleotide polymorphism (SNP) in the imprinted DLK1-DIO3 domain on ovine chromosome 18. The callipyge SNP up-regulates the expression of surrounding transcripts when inherited in cis without altering their allele-specific imprinting status. The callipyge phenotype exhibits polar overdominant inheritance since only paternal heterozygous animals have muscle hypertrophy. Two studies were conducted profiling gene expression in lamb muscles to determine the down-stream effects of over-expression of paternal allele-specific DLK1 and RTL1 as well as maternal allele-specific MEG3, RTL1AS and MEG8, using Affymetrix bovine expression arrays. A total of 375 transcripts were differentially expressed in callipyge muscle and 25 transcripts were subsequently validated by quantitative PCR. The muscle-specific expression patterns of most genes were similar to DLK1 and included genes that are transcriptional repressors or affect feedback mechanisms in beta-adrenergic and growth factor signaling pathways. One gene, phosphodiesterase 7A had an expression pattern similar to RTL1 expression indicating a biological activity for RTL1 in muscle. Only transcripts that localize to the DLK1-DIO3 domain were affected by inheritance of a maternal callipyge allele. Callipyge sheep are a unique model to study over expression of both paternal allele-specific genes and maternal allele-specific non-coding RNA with an accessible and nonlethal phenotype. This study has identified a number of genes that are regulated by DLK1 and RTL1 expression and exert control on postnatal skeletal muscle growth. The genes identified in this model are primary candidates for naturally regulating postnatal muscle growth in all meat animal species, and may serve as targets to ameliorate muscle atrophy conditions including myopathic diseases and age-related sarcopenia.