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

An osteocalcin-deficient mouse strain without endocrine abnormalities

Osteocalcin (OCN), the most abundant noncollagenous protein in the bone matrix, is reported to be a bone-derived endocrine hormone with wide-ranging effects on many aspects of physiology, including glucose metabolism and male fertility. Many of these observations were made using an OCN-deficient mouse allele (Osc– ) in which the 2 OCN-encoding genes in mice, Bglap and Bglap2, were deleted in ES cells by homologous recombination. Here we describe mice with a new Bglap and Bglap2 double-knockout (dko) allele (Bglap/2p.Pro25fs17Ter) that was generated by CRISPR/Cas9-mediated gene editing. Mice homozygous for this new allele do not express full-length Bglap or Bglap2 mRNA and have no immunodetectable OCN in their serum. FTIR imaging of cortical bone in these homozygous knockout animals finds alterations in the collagen maturity and carbonate to phosphate ratio in the cortical bone, compared with wild-type littermates. However, μCT and 3-point bending tests do not find differences from wild-type littermates with respect to bone mass and strength. In contrast to the previously reported OCN-deficient mice with the Osc− allele, serum glucose levels and male fertility in the OCN-deficient mice with the Bglap/ 2pPro25fs17Ter allele did not have significant differences from wild-type littermates. We cannot explain the absence of endocrine effects in mice with this new knockout allele. Possible explanations include the effects of each mutated allele on the transcription of neighboring genes, or differences in genetic background and environment. So that our findings can be confirmed and extended by other interested investigators, we are donating this new Bglap and Bglap2 double-knockout strain to the Jackson Laboratories for academic distribution

Investigating Transcriptional Regulation Within Bone Development: Characterization of HRPT2/CDC73 In Vivo and the Effects of Ascorbic Acid in Osteoblast Differentiation In Vitro

The integrity of the mechanisms that control gene transcription during development and in post-natal life is essential to maintain tissue homeostasis and impede the development of genetic diseases such as cancer. Inheritance of a defective hyperpatahyroidism 2 (HRPT2) allele, an essential regulator of gene transcription, predisposes individuals to a constellation of symptoms ranging from endocrine abnormalities to parathyroid adenomas and jaw bone tumors (HPT-JT). In order to elucidate the function of the HRPT2 gene and the pathogenesis that results upon spontaneous inactivation in familial cases of parathyroid cancer and HPT-JT, mouse models were generated that allow for deletion of Hrpt2 within different stages and tissues during development. We have used the Hrpt2 flox mouse model to delete Hrpt2 in mesenchymal progenitor cells as well as committed, terminally differentiated osteoblasts and osteocytes. Whereas loss of Hrpt2 in mesenchymal progenitors was embryonic lethal, genetic deletion of Hrpt2 in mature bone forming cells led to increased bone mass and bone strength. However, Hrpt2 conditional knockout bones had increased cortical porosity and osteocyte apoptosis associated with increased osteoblast specific gene expression. This work suggests that Hrpt2 is required for cell proliferation and differentiation and acts as a transcriptional repressor in terminally differentiated cell types. Control of gene transcription defines cell identity and fate. Ascorbic acid (AA, also known as vitamin C) is an essential vitamin for humans and is well known for its role in collagen synthesis. AA acts as a cofactor for TET enzymes, which hydroxylate methylated cytosines. Here, we characterize how 7 days of AA treatment causes changes in gene transcription, 5-hydroxymethylcytosine deposition, and the active chromatin marks H3K4me3 and H3k27ac in MC3T3-E1 murine pre-osteoblasts cells, initiating cell differentiation and expression of the osteoblast phenotype. Though 5hmC deposition was not specific for only highly expressed genes, it was highly enriched at transcriptional start sites and CpG islands. While H3K4me3 was mostly unchanged, H3K27ac was predictive of driving gene expression. This work suggests that AA causes dramatic changes to the epigenome through epigenetic modifiers to impact cell differentiation

MicroRNAs Encoded by Bovine Leukemia Virus (BLV) Are Associated with Reduced Expression of B Cell Transcriptional Regulators in Dairy Cattle Naturally Infected with BLV

Bovine leukemia virus (BLV) is estimated to infect over 83% of dairy herds and over 40% of all dairy cows in the United States. While, BLV only causes leukemia in a small proportion of animals, research indicates that BLV+ cattle exhibit reduced milk production and longevity that is distinct from lymphoma development. It is hypothesized that BLV negatively affects production by interfering with cattle immunity and increasing the risk of secondary infections. In particular, BLV+ cows demonstrate reduced circulating levels of both antigen-specific and total IgM. This study investigated possible mechanisms by which BLV could interfere with the production of IgM in naturally infected cattle. Specifically, total plasma IgM and the expression of genes IGJ, BLIMP1, BCL6, and PAX5 in circulating IgM+ B cells were measured in 15 naturally infected BLV+ and 15 BLV− cows. In addition, BLV proviral load (PVL) (a relative measurement of BLV provirus integrated into host DNA) and the relative expression of BLV TAX and 5 BLV microRNAs (miRNAs) were characterized and correlated to the expression of selected endogenous genes. BLV+ cows exhibited lower total plasma IgM and lower expression of IGJ, BLIMP1, and BCL6. While, BLV TAX and BLV miRNAs failed to correlate with IGJ expression, both BLV TAX and BLV miRNAs exhibited negative associations with BLIMP1 and BCL6 gene expression. The results suggest a possible transcriptional pathway by which BLV interferes with IgM production in naturally infected cattle

Dairy Cows Naturally Infected with Bovine Leukemia Virus Exhibit Abnormal B- and T-Cell Phenotypes after Primary and Secondary Exposures to Keyhole Limpet Hemocyanin

Bovine leukemia virus (BLV) is a retrovirus that is highly prevalent in US dairy herds: over 83% are BLV infected and the within-herd infection rate can be almost 50% on average. While BLV is known to cause lymphosarcomas, only 5% or fewer infected cattle will develop lymphoma; this low prevalence of cancer has historically not been a concern to dairy producers. However, more recent research has found that BLV+ cows without lymphoma produce less milk and have shorter lifespans than uninfected herdmates. It has been hypothesized that BLV infection interferes with normal immune function in infected cattle, and this could lead to reduced dairy production. To assess how naturally infected BLV+ cows responded to a primary and secondary immune challenge, 10 BLV+ and 10 BLV− cows were injected subcutaneously with keyhole limpet hemocyanin (KLH) and dimethyldioctadecylammonium bromide. B- and T-cell responses were characterized over the following 28 days. A total of 56 days after primary KLH exposure, cows were re-injected with KLH and B- and T-cell responses were characterized again over the following 28 days. BLV+ cows produced less KLH-specific IgM after primary immune stimulation; demonstrated fewer CD45R0+ B cells, altered proportions of CD5+ B cells, altered expression of CD5 on CD5+ B cells, and reduced MHCII surface expression on B cells ex vivo; exhibited reduced B-cell activation in vitro; and displayed an increase in BLV proviral load after KLH exposure. In addition, BLV+ cows had a reduced CD45R0+γδ+ T-cell population in the periphery and demonstrated a greater prevalence of IL4-producing T cells in vitro. All together, our results demonstrate that both B- and T-cell immunities are disrupted in BLV+ cows and that antigen-specific deficiencies can be detected in BLV+ cows even after a primary immune exposure

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