Partial molecular characterization of different proviral strains of bovine leukemia virus

Molecular investigation of possible relationships concerning bovine leukemia virus and breast cancer

Worldwide, breast cancer has an eminent morbidity and mortality rate, as it is a neoplastic disease among females. The query of the prospective danger of bovine leukemia virus (BLV) to humans is an old but exceedingly topical focus of scientific debate. The objective of the current study was to determine the possible relationship between BLV and breast cancer. A total of 2710 formalin-fixed paraffin-embedded (FFPE) breast cancer samples were selected regardless of the age, ethnicity, or municipality origin of the subjects. The presence of BLV in human breast cancer was determined through nested PCR by amplifying tax and gag genes followed by partial sequencing. Homology was confirmed by using the online BLAST Tool. BLV genes were found to be positive in 26.8% (728/2710) of the samples from breast cancer patients and 10% (10/80) of the samples without cancer (negative control). The results indicated a correlation between the presence of the BLV gene and breast cancer (odds ratio = 0.3889; confidence interval = 1,18; p = 0.0029). The current findings suggest a possible link between BLV and human breast carcinoma. Therefore, screening cattle herds and milk products is suggested to reduce the viral transmission risk to humans.

Lack of association between amino acid sequences of the bovine leukemia virus envelope and varying stages of infection in dairy cattle

We collected 724 blood samples from dairy cattle from six Mexican states, and tested them for the presence of antibodies against BLV using a commercial ELISA test. Our study groups consisted of 32 samples: 12 asymptomatic cows, 12 cows with lymphocytosis and 8 samples of tumor tissue of the abomasum and heart of cattle with lymphoma. We designed three pairs of primers to amplify the complete BLV env gene, and obtained a fragment of 1548 nucleotides in length with the sequenced products. According to the phylogenetic tree we constructed to identify the viral genotype, 96.87 % of the sequences grouped into genotype 1, while a single sample from a cow with lymphocytosis (3.13 %) was associated with genotype 3 sequences. The similarity between the Mexican BLV sequences ranged from 0.985-1.00. In addition, the proportion of non-synonymous and synonymous mutations indicated negative selection. We did not identify any conserved residues in the viral protein sequences that could be related to BLV infection stage in cattle. Proviral quantification was performed using quantitative polymerase chain reaction, and we used Mood´s median test as statistical analysis. We found no significant association between proviral load and phase of infection. The sequences showed high similarity without any association between BLV surface glycoprotein and the different infection stages, nor differences in the proviral load. BLV genotype 1 was identified as prevalent in the studied samples, and for the first time in Mexico, we identified BLV genotype 3 in cattle.

Early events following bovine leukaemia virus infection in calves with different alleles of the major histocompatibility complex DRB3 gene

Cattle maintaining a low proviral load (LPL) status after bovine leukaemia virus (BLV) infection have been recognized as BLV controllers and non-transmitters to uninfected cattle in experimental and natural conditions. LPL has been associated with host genetics, mainly with the BoLA class II DRB3 gene. The aim of this work was to study the kinetics of BLV and the host response in Holstein calves carrying different BoLA-DRB3 alleles. Twenty BLV-free calves were inoculated with infected lymphocytes. Two calves were maintained uninfected as controls. Proviral load, total leukocyte and lymphocyte counts, anti-BLVgp51 titres and BLVp24 expression levels were determined in blood samples at various times post-inoculation. The viral load peaked at 30 days post-inoculation (dpi) in all animals. The viral load decreased steadily from seroconversion (38 dpi) to the end of the study (178 dpi) in calves carrying a resistance-associated allele (*0902), while it was maintained at elevated levels in calves with *1501 or neutral alleles after seroconversion. Leukocyte and lymphocyte counts and BLVp24 expression did not significantly differ between genetic groups. Animals with < 20 proviral copies/30 ng of DNA at 178 dpi or < 200 proviral copies at 88 dpi were classified as LPL, while calves with levels above these limits were considered to have high proviral load (HPL) profiles. All six calves with the *1501 allele progressed to HPL, while LPL was attained by 6/7 (86%) and 2/6 (33%) of the calves with the *0902 and neutral alleles, respectively. One calf with both *0902 and *1501 developed LPL. This is the first report of experimental induction of the LPL profile in cattle.

Bovine leukemia virus long terminal repeat variability: identification of single nucleotide polymorphisms in regulatory sequences

Background

Limited data are available on the incidence of variations in nucleotide sequences of long terminal repeat (LTR) regions of Bovine Leukemia Virus (BLV). Consequently, the possible impact of SNPs on BLV LTR function are poorly elucidated. Thus, a detailed and representative study of full-length LTR sequences obtained from sixty-four BLV isolates from different geographical regions of Poland, Moldova, Croatia, Ukraine and Russia were analyzed for their genetic variability.

Methods

Overlap extension PCR, sequencing and Bayesian phylogenetic reconstruction of LTR sequences were performed. These analyses were followed by detailed sequence comparison, estimation of genetic heterogeneity and identification of transcription factor binding site (TFBS) modifications.

Results

Phylogenetic analysis of curated LTR sequences and those available in the GenBank database reflected the acknowledged env gene classification of BLV into 10 genotypes, and further clustered analysed sequences into three genotypes - G4, G7 and G8. Additional molecular studies revealed the presence of 97 point mutations distributed at 89 positions throughout all 64 LTR sequences. The highest rate of variability was noted in U3 and U5 subregions. However, the variability in regulatory sequences (V_R) was assessed as lower than the variability within non-regulatory sequences (V_NR) for both, U3 and U5 subregions. In contrast, V_R value for R subregion, as well as for the total LTR, was higher than the V_NR suggesting the existence of positive selection. Twelve unique SNPs for these LTR sequences localized in regulatory and non-regulatory elements were identified. The presence of different types of substitutions lead to the abrogation of present or to the creation of additional TFBS.

Conclusion

This study represents the largest study of LTR genetic variability of BLV field isolates from Eastern part of Europe. Phylogenetic analysis of LTRs supports the clustering BLV variants based on their geographic origin. The SNP screening showed variations modifying LTR regulatory sequences, as well as altering TFBS. These features warrant further exploration as they could be related to proviral load and distinctive regulation of BLV transcription and replication.

Electronic supplementary material

The online version of this article (10.1186/s12985-018-1062-z) contains supplementary material, which is available to authorized users.

Heterogeneity determination of bovine leukemia virus genome in Santa Catarina state, Brazil

ABSTRACT: Bovine leukemia virus (BLV) is a member of Retroviridae family, genus Deltaretrovirus, and the main viral agent responsible for economic loses in dairy herds. Some studies have been carried out about BLV genotypes, and at least seven genotypes were found out in samples of different regions of the world. The objective of this study was to identify BLV samples from seropositive dairy cattle in Santa Catarina state, Brazil, using molecular techniques. Blood samples were collected (454) from dairy cattle from 31 different farms, and serology using agar gel immunodiffusion test (AGID) was performed. After that, 191 seropositive samples were submitted to DNA extraction, and in 77 samples the polymerase chain reaction (PCR) for amplification of a 440 bp fragment of the env gene was performed. Nineteen DNA samples were subjected to restriction fragment length polymorphism (RFLP) analysis by digestion of the PCR fragment by five restriction endonucleases - BamHI, HaeIII, Tru9I, TaqI, and MwoI. It was found 42% seropositive animals (191/454) and 68% positives of the farms (21/31). The PCR showed 80.5% (62/77) of animals positive. The RFLP analysis identified five different genotypes dispersed by Santa Catarina state, with the highest prevalence for genotype X (47.4%). Overall, our results identified the viral genotypes present in dairy cattle and the prevalence of new variants in representative farms from Santa Catarina state.

Stable infection of a bovine mammary epithelial cell line (MAC-T) with bovine leukemia virus (BLV)

Bovine leukemia virus (BLV) is a retrovirus that affects cattle causing a lymphoproliferative disease. BLV infection has been associated with misbalance of the immune response causing a higher incidence of other infections. Mastitis is one of the most important conditions that affect milk production in cattle. The aim of this study was to stably infect a bovine mammary epithelial cell line (MAC-T). MAC-T cell line was successfully infected with BLV and the infection was confirmed by nested PCR, qPCR, immunocytochemistry, western blot and transmission electron microscopy. This is the first report of a bovine mammary epithelial cell line stably infected with BLV. This new cell line could be used as an in vitro model to study the effect of BLV on the immune response in the mammary gland and the relationship with other agents causing mastitis.

Genetic Diversity of Bovine Major Histocompatibility Complex Class II DRB3 locus in cattle breeds from Asia compared to those from Africa and America

Genetic polymorphisms and diversity of BoLA-DRB3.2 are essential because of DRB3 gene's function in innate immunity and its association with infectious diseases resistance or tolerance in cattle. The present study was aimed at assessing the level of genetic diversity of DRB3 in the exon 2 (BoLA-DRB3.2) region in African, American and Asian cattle breeds. Amplification of exon 2 in 174 cattle revealed 15 haplotypes. The breeds with the highest number of haplotypes were Brangus (10), Sokoto Gudali (10) and Dajal (9), while the lowest number of haplotypes were found in Holstein and Sahiwal with 4 haplotypes each. Medium Joining network obtained from haplotypic data showed that all haplotypes condensed around a centric area and each sequence (except in H-3, H-51 and H-106) representing almost a specific haplotype. The BoLA-DRB3.2 sequence analyses revealed a non-significant higher rate of non-synonymous (dN) compared to synonymous substitutions (dS). The ratio of dN/dS substitution across the breeds were observed to be greater than one suggesting that variation at the antigen-binding sites is under positive selection; thus increasing the chances of these breeds to respond to wide array of pathogenic attacks. An analysis of molecular variance revealed that 94.01 and 5.99% of the genetic variation was attributable to differences within and among populations, respectively. Generally, results obtained suggest that within breed genetic variation across breeds is higher than between breeds. This genetic information will be important for investigating the relationship between BoLADRB3.2 and diseases in various cattle breeds studied with attendant implication on designing breeding programs that will aim at selecting individual cattle that carry resistant alleles.

Computational analysis of envelope glycoproteins from diverse geographical isolates of bovine leukemia virus identifies highly conserved peptide motifs

BACKGROUND

Bovine leukemia virus (BLV) is a deltaretrovirus infecting bovine B cells and causing enzootic bovine leucosis. The SU or surface subunit, gp51, of its envelope glycoprotein is involved in receptor recognition and virion attachment. It contains the major neutralizing and CD4+ and CD8+ T cell epitopes found in naturally infected animals. In this study, we aimed to determine global variation and conservation within gp51 in the context of developing an effective global BLV vaccine.

RESULTS

A total of 256 sequences extracted from the NCBI database and collected in different parts of the world, were studied to identify conserved segments along the env gene sequences that encode the gp51 protein. Using the MEME server and the conserved DNA Region module for analysis within DnaSP, we identified six conserved segments, referred to as A-F, and five semi-conserved segments, referred to as G-K. The amino acid conservation ranged from 98.8 to 99.8% in conserved segments A to F, while segments G to K had 89.6-95.2% conserved amino acid sequence. Selection analysis of individual segments revealed that residues of conserved segments had undergone purifying selection, whereas, particular residues in the semi-conserved segments are currently undergoing positive selection, specifically at amino acid positions 48 in segment K, 74 in segment G, 82 in segment I, 133 and 142 in segment J, and residue 291 in segment H. Each of the codons for these six residues contain the most highly variable nucleotides within their respective semi-conserved segments.

CONCLUSIONS

The data described here show that the consensus amino acid sequence constitutes a strong candidate from which a global vaccine can be derived for use in countries where eradication by culling is not economically feasible. The most conserved segments overlap with amino acids in known immunodeterminants, specifically in epitopes D-D', E-E', CD8+ T-cell epitopes, neutralizing domain 1 and CD4+ T-cell epitopes. Two of the segments reported here represent unique segments that do not overlap with previously identified antigenic determinants. We propose that evidence of positive selection in some residues of the semi-conserved segments suggests that their variation is involved in viral strategy to escape immune surveillance of the host.

A new genotype of bovine leukemia virus in South America identified by NGS-based whole genome sequencing and molecular evolutionary genetic analysis

BACKGROUND

Bovine leukemia virus (BLV) is a member of retroviridae family, together with human T cell leukemia virus types 1 and 2 (HTLV-1 and -2) belonging to the genes deltaretrovirus, and infects cattle worldwide. Previous studies have classified the env sequences of BLV provirus from different geographic locations into eight genetic groups. To investigate the genetic variability of BLV in South America, we performed phylogenetic analyses of whole genome and partial env gp51 sequences of BLV strains isolated from Peru, Paraguay and Bolivia, for which no the molecular characteristics of BLV have previously been published, and discovered a novel BLV genotype, genotype-9, in Bolivia.

RESULTS

In Peru and Paraguay, 42.3 % (139/328) and over 50 % (76/139) of samples, respectively, were BLV positive. In Bolivia, the BLV infection rate was up to 30 % (156/507) at the individual level. In Argentina, 325/420 samples were BLV positive, with a BLV prevalence of 77.4 % at the individual level and up to 90.9 % at herd level. By contrast, relatively few BLV positive samples were detected in Chile, with a maximum of 29.1 % BLV infection at the individual level. We performed phylogenetic analyses using two different approaches, maximum likelihood (ML) tree and Bayesian inference, using 35 distinct partial env gp51 sequences from BLV strains isolated from Peru, Paraguay, and Bolivia, and 74 known BLV strains, representing eight different BLV genotypes from various geographical locations worldwide. The results indicated that Peruvian and Paraguayan BLV strains were grouped into genotypes-1, -2, and -6, while those from Bolivia were clustered into genotypes-1, -2, and -6, and a new genotype, genotype-9. Interestingly, these results were confirmed using ML phylogenetic analysis of whole genome sequences obtained by next generation sequencing of 25 BLV strains, assigned to four different genotypes (genotypes-1, -2, -6, and -9) from Peru, Paraguay, and Bolivia. Comparative analyses of complete genome sequences clearly showed some specific substitutions, in both structural and non-structural BLV genes, distinguishing the novel genotype-9 from known genotypes.

CONCLUSIONS

Our results demonstrate widespread BLV infection in South American cattle and the existence of a new BLV genotype-9 in Bolivia. We conclude that at least seven BLV genotypes (genotypes-1, -2, -4, -5, -6, -7, and -9) are circulating in South America.

Animal models of bovine leukemia virus and human T-lymphotrophic virus type-1: insights in transmission and pathogenesis

Bovine leukemia virus (BLV) and human T-lymphotrophic virus type-1 (HTLV-1) are related retroviruses associated with persistent and lifelong infections and a low incidence of lymphomas within their hosts. Both viruses can be spread through contact with bodily fluids containing infected cells, most often from mother to offspring through breast milk. Each of these complex retroviruses contains typical gag, pol, and env genes but also unique, nonstructural proteins encoded from the pX region. These nonstructural genes encode the Tax and Rex regulatory proteins, as well as novel proteins essential for viral spread in vivo. Improvements in the molecular tools to test these viral determinants in cellular and animal models have provided new insights into the pathogenesis of each virus. Comparisons of BLV and HTLV-1 provide insights into mechanisms of spread and tumor formation, as well as potential approaches to therapeutic intervention against the infections.