Influence of Maternal BLV Infection on miRNA and tRF Expression in Calves.

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.

No evidence for a negative association between bovine leukemia virus status and fertility in Kansas beef herds: a cross-sectional study.

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.

Cross-sectional study to describe bovine leukemia virus herd and within-herd ELISA prevalence and bovine leukemia virus proviral load of convenience-sampled Kansas beef cow-calf herds.

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.

Identification of BoLA Alleles Associated with BLV Proviral Load in US Beef Cows.

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.

Phenotypic Selection of Dairy Cattle Infected with Bovine Leukemia Virus Demonstrates Immunogenetic Resilience through NGS-Based Genotyping of BoLA MHC Class II Genes.

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.

Tracing Viral Transmission and Evolution of Bovine Leukemia Virus through Long Read Oxford Nanopore Sequencing of the Proviral Genome.

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.

Current Developments in the Epidemiology and Control of Enzootic Bovine Leukosis as Caused by Bovine Leukemia Virus.

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.

Analysis of tRNA halves (tsRNAs) in serum from cattle challenged with bovine viral diarrhea virus.

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.

Association of Transfer RNA Fragments in White Blood Cells With Antibody Response to Bovine Leukemia Virus in Holstein Cattle.

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. tRF5GlnCTG, tRF5GlnTTG, and tRF5HisGTG, 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.

Circulating MicroRNAs in Serum from Cattle Challenged with Bovine Viral Diarrhea Virus.

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.

Diet shifts provoke complex and variable changes in the metabolic networks of the ruminal microbiome.

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.

The players may change but the game remains: network analyses of ruminal microbiomes suggest taxonomic differences mask functional similarity.

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.

Accuracies of genomic breeding values in American Angus beef cattle using K-means clustering for cross-validation.

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.

Effect of DLK1 and RTL1 but not MEG3 or MEG8 on muscle gene expression in Callipyge lambs.

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.