BoLA-DRB3 Polymorphism is Associated with Differential Susceptibility to Bovine Leukemia Virus-Induced Lymphoma and Proviral Load

Bovine leukemia virus (BLV) is the causative agent of enzootic bovine leucosis. However, less than 5% of BLV-infected cattle will develop lymphoma, suggesting that, in addition to viral infection, host genetic polymorphisms might play a role in disease susceptibility. Bovine leukocyte antigen (BoLA)-DRB3 is a highly polymorphic gene associated with BLV proviral load (PVL) susceptibility. Due to the fact that PVL is positively associated with disease progression, it is believed that controlling PVL can prevent lymphoma development. Thus, many studies have focused on the relationship between PVL and BoLA-DRB3. Despite this, there is little information regarding the relationship between lymphoma and BoLA-DRB3. Furthermore, whether or not PVL-associated BoLA-DRB3 is linked to lymphoma-associated BoLA-DRB3 has not been clarified. Here, we investigated whether or not lymphoma-associated BoLA-DRB3 is correlated with PVL-associated BoLA-DRB3. We demonstrate that two BoLA-DRB3 alleles were specifically associated with lymphoma resistance (*010:01 and *011:01), but no lymphoma-specific susceptibility alleles were found; furthermore, two other alleles, *002:01 and *012:01, were associated with PVL resistance and susceptibility, respectively. In contrast, lymphoma and PVL shared two resistance-associated (DRB3*014:01:01 and *009:02) BoLA-DRB3 alleles. Interestingly, we found that PVL associated alleles, but not lymphoma associated alleles, are related with the anti-BLV gp51 antibody production level in cows. Overall, our study is the first to demonstrate that the BoLA-DRB3 polymorphism confers differential susceptibility to BLV-induced lymphoma and PVL.


Introduction
Viral load in chronic infections with viruses, such as hepatitis B virus (HBV), hepatitis C virus (HCV), human T cell leukemia virus type 1 (HTLV-1), and human immunodeficiency virus type 1 (HIV-1), has been reported to determine the likelihood of pathogenesis and disease progression [1][2][3][4]. For retroviruses, whose genome integrates with the host genome, proviral load (PVL) is an important risk factor of viruss-associated disease prediction [5,6]. Bovine leukemia virus (BLV) is closely related to HTLV-1 and is the causative agent of enzootic bovine leukosis (EBL), a disease that is characterized by long-term symptoms, including persistent lymphocytosis (PL), which may culminate in B-cell lymphosarcoma [7,8]. Several studies indicate that BLV PVL is associated with BLV-related disease progression [9][10][11][12][13]. However, only 5% of infected cattle progress to develop lymphoma, suggesting that in addition to viral infection, host genetic polymorphisms might play a role in disease susceptibility.
The major histocompatibility complex (MHC), a highly polymorphic gene set, plays a crucial role in antigen presentation and immune responsiveness [14][15][16], and thus, it is associated with numerous infectious diseases. In cattle, the MHC system is known as bovine leukocyte antigen (BoLA). Several studies have identified genetic variations in BoLA-DRB3, a functionally important locus and the most highly polymorphic BoLA class II locus in cattle. To date, 330 DRB3 alleles have been registered in the Immuno Polymorphism Database (IPD)-MHC database (https://www.ebi.ac. uk/ipd/mhc/group/BoLA/). The BoLA-DRB3 polymorphism influences susceptibility to BLV-induced lymphoma [17][18][19], and to PVL [20][21][22]. As PVL is positively related to lymphoma development, it is possible that lymphoma-associated BoLA-DRB3 is consistent with PVL-associated BoLA-DRB3. However, the consistency of the above association has not been studied yet.
Indeed, lymphoma development and viral replication depend on different cellular mechanisms, potentially leading to the differential susceptibility of lymphoma and PVL to BoLA-DRB3. It has been reported that proviral integration and BLV proteins are required for initial cell transformation [7,23]. However, the host immune system can remove transformed cells by lymphocyte activation via MHC molecules [24]. Due to the fact that MHC class II alleles affect antigen presentation and MHC expression levels in cancer cells [25,26], it is reasonable to hypothesize that MHC class II alleles would bind to peptides derived from viral or tumor antigens, and that the resulting complex would be recognized by CD4+ T cells. Consequently, some BoLA-DRB3 alleles might specifically bind with the processed viral antigen, while others might specifically recognize the tumor antigens. Thus, it is likely that different BoLA-DRB3 alleles are specifically associated with BLV-induced lymphoma and PVL. Consistent with this, PVL does not always correlate with lymphoma development, as many infected cows with a high PVL do not develop lymphoma. On the contrary, attenuated BLV-infected sheep were found to exhibit significantly lower PVL, but still developed lymphoma [27], suggesting that lymphoma and PVL may induce different susceptibilities, depending on different BoLA-DRB3 polymorphisms. In this study, using asymptomatic and lymphoma Holstein cows randomly collected in a nationwide survey in Japan, we demonstrated that BoLA-DRB3 polymorphism is associated with differential susceptibility to BLV-induced lymphoma and PVL.

Sample Collection and Diagnosis
Blood samples from 611 BLV-infected but clinically normal Holstein cows (asymptomatic cows; information summarized in Tables S1 and S6) and 221 BLV-infected Holstein cows with lymphoma (lymphoma cows; information summarized in Table S2) were randomly collected in a nationwide Viruses 2020, 12, 352 3 of 15 survey across Japan (32 prefectures out of 47), and the genomic DNA and plasma from peripheral blood were isolated. The subclinical stage of BLV infection was diagnosed according to the lymphocyte count (cells/µL) and the age of each cow (≤ 8,500 = normal and ≥ 13,000 = lymphocytosis for cows aged 2-3 years; ≤ 5,500 = normal and ≥ 7,500 = lymphocytosis for cows aged ≥ 6 years). Asymptomatic cows were defined as BLV-infected but clinically and hematologically normal cattle; PL cows were defined as BLV-infected but clinically normal cattle showing with an increase in the number of apparently normal B lymphocytes. Subsequently, lymphoma was diagnosed by both gross and histological observation and by detecting atypical mononuclear cells in the slaughterhouse. In this study, PL cases were excluded and used only samples from asymptomatic cows and lymphoma cows.

Detection of Anti-BLV gp51 Antibody by Enzyme-Linked Immunosorbent Assay (ELISA)
The anti-BLV gp51 antibody was measured with an anti-BLV antibody ELISA Kit (JNC, Tokyo, Japan), according to the manufacturer's instructions. Two-fold serial dilutions of plasma samples starting at 1:16 were tested by the ELISA Kit. The OD value in each DRB3 group was compared at each dilution.

Association Study and Statistical Analysis
An association study based on Fisher's exact test was performed by comparing the allele and genotype frequencies between asymptomatic and lymphoma cows or low PVL and high PVL cows. The results were penalized with the Benjamini-Hochberg (BH) procedure to correct for the false positive rate. Each allele or genotype was ranked based on their p-value starting from the smallest one. The BH value was calculated based on the equation p-value rank / total allele (genotype) number × 0.05). The alleles or genotypes with p-value < BH value and odds ratio (OR) < 1 were categorized as resistance alleles. In contrast, those with p-value < BH value and OR > 1 were defined as susceptibility alleles or genotypes. The association of cow mean age and birth location with lymphoma or PVL was evaluated by the Mann-Whitney U test and Tukey's multiple comparison test, respectively. When we confirmed the association between age and lymphoma or PVL, we performed logistic regression analyses to adjust for age. To evaluate a multiplicative interaction between BoLA-DRB3 alleles, we introduced the interaction term in a logistic regression model as conditional analysis [32]. We assessed the significance level of the association study by applying a Bonferroni correction according to the number of assessed alleles (adjusted p < 0.05). All calculations were performed using R software (version 3.5.0, R Foundation for Statistical Computing, Vienna, Austria).

PVL Is Not Fully Correlated with Lymphoma Development
The PVL of 250 asymptomatic cows (Table S1) ranged from 5-120,482 copies/10 5 cells (mean: 9401 copies/10 5 cells), while in 221 lymphoma cows (Table S2), the PVL ranged from 28-1,960,674 copies/10 5 cells (mean 99,522 copies/10 5 cells; Figure 1). This difference suggested that animals with a high BLV PVL were at a higher risk of developing lymphoma. Our previous report indicated that cows with a PVL of greater than 14,000 copies/10 5 cells secreted BLV into nasal mucus [30], and BLV provirus was detected in milk samples from cows when the PVL in blood samples was higher than 10,000 copies/10 5 cells [33]. These results suggest that a PVL around 10,000 copies/10 5 cells in blood might be an indicator of efficient BLV spreading within the whole body and thus this is a relatively high number. Therefore, a BLV PVL of 10,000 copies/10 5 cells was set as a threshold to distinguish between high-PVL (HPVL) and low-PVL (LPVL) cows (Figure 1), which is also in line with our previous study [19]. Consistently, in lymphoma and asymptomatic cows, the mean PVL was found to be above and below this threshold, respectively ( Figure 1). However, 62 HPVL cows remained asymptomatic, whereas 37 LPVL cows developed lymphoma (

PVL Is Not Fully Correlated with Lymphoma Development
The PVL of 250 asymptomatic cows (Table S1) ranged from 5-120,482 copies/10 5 cells (mean: 9401 copies/10 5 cells), while in 221 lymphoma cows (Table S2), the PVL ranged from 28-1,960,674 copies/10 5 cells (mean 99,522 copies/10 5 cells; Figure 1). This difference suggested that animals with a high BLV PVL were at a higher risk of developing lymphoma. Our previous report indicated that cows with a PVL of greater than 14,000 copies/10 5 cells secreted BLV into nasal mucus [30], and BLV provirus was detected in milk samples from cows when the PVL in blood samples was higher than 10,000 copies/10 5 cells [33]. These results suggest that a PVL around 10,000 copies/10 5 cells in blood might be an indicator of efficient BLV spreading within the whole body and thus this is a relatively high number. Therefore, a BLV PVL of 10,000 copies/10 5 cells was set as a threshold to distinguish between high-PVL (HPVL) and low-PVL (LPVL) cows (Figure 1), which is also in line with our previous study [19]. Consistently, in lymphoma and asymptomatic cows, the mean PVL was found to be above and below this threshold, respectively ( Figure 1). However, 62 HPVL cows remained asymptomatic, whereas 37 LPVL cows developed lymphoma ( Table 1), indicating that lymphoma development is not fully correlated with the PVL. This could be because BLV-induced lymphoma and BLV PVLs are associated with different BoLA-DRB3 alleles.  (Table S1) and 221 lymphoma (Table S2) cows in a nationwide survey in Japan. BLV infection was analyzed using BLV-CoCoMo-qPCR-2. The red line represents a BLV PVL of 10,000 copies/10 5 cell, which was set as the threshold between high-and low-PVL cows.
To exclude the effect from other potential factors that might associate with lymphoma development such as cow origin and age, we then applied multivariable logistic regression to adjust for the effect of these potential factors. Our calculation indicates that age showed a significant association with disease susceptibility (p = 6.56 × 10 −6 ). However, no significant difference was observed between location and lymphoma susceptibility (p = 0.182). Therefore, we performed the logistic regression analysis adjusted by age in the only association study of lymphoma. After studying the association between each BoLA-DRB3 allele and lymphoma susceptibility, we conducted a stepwise conditional analysis with respect to the top-associated BoLA-DRB3 alleles (Table S4). In Table 3, a conditional analysis of DRB3*011:01 revealed an independent association with DRB3*009:02 and DRB3*010:01. A subsequent conditional analysis regarding DRB3*009:02 and DRB3*011:01 revealed an independent association with  (Table S3); 13 out of 32 alleles with frequency > 1% are shown. The X-axis shows the allele name and the Y-axis shows allele frequency (%) for each BoLA-DRB3 allele. The Benjamini-Hochberg (BH) procedure was performed to adjust the false positive rate. Alleles with a p-value < BH value were defined as susceptibility (S) with an odds ratio (OR) > 1 and as resistance (R) with an OR < 1. BH value = (p-value rank/total allele number) × 0.05.
To exclude the effect from other potential factors that might associate with lymphoma development such as cow origin and age, we then applied multivariable logistic regression to adjust for the effect of these potential factors. Our calculation indicates that age showed a significant association with disease susceptibility (p = 6.56 × 10 −6 ). However, no significant difference was observed between location and lymphoma susceptibility (p = 0.182). Therefore, we performed the logistic regression analysis adjusted by age in the only association study of lymphoma. After studying the association between each BoLA-DRB3 allele and lymphoma susceptibility, we conducted a stepwise conditional analysis with respect to the top-associated BoLA-DRB3 alleles (Table S4). In Table 3, a conditional analysis of DRB3*011:01 revealed an independent association with DRB3*009:02 and DRB3*010:01. A subsequent conditional analysis regarding DRB3*009:02 and DRB3*011:01 revealed an independent association with DRB3*010:01 and DRB3*014:01:01. Next, the subsequent conditional analysis of DRB3*009:02, DRB3*010:01, and DRB3*011:01 revealed an independent association with DRB3*014:01:01. After conditioning DRB3*009:02, DRB3*010:01, DRB3*011:01, and DRB3*014:01:01, no significant association locus was observed. We then conducted a multivariate regression analysis incorporating the four associated BoLA-DRB3 alleles (DRB3*009:02, DRB3*010:01, DRB3*011:01, and DRB3*014:01:01). We identified that all were independently associated with lymphoma resistance (Table 3). Table 3. Logistic regression analysis-based association study of BoLA-DRB3 alleles in asymptomatic and lymphoma cows after adjustments for age. For the genotype association study, genotypes with frequencies > 1% are shown in Figure 3 (complete genotype frequencies are summarized in Table S5). However, no genotypes reached statistical significance in terms of their association with lymphoma development after BH correction to adjust the false discovery rate (Table 4). For the genotype association study, genotypes with frequencies > 1% are shown in Figure 3 (complete genotype frequencies are summarized in Table S5). However, no genotypes reached statistical significance in terms of their association with lymphoma development after BH correction to adjust the false discovery rate (Table 4).  (Table S4); a total of 31 out of 94 genotypes with frequency > 1% are shown. The X-axis and Y-axis show the genotype name and frequency (%) for each BoLA-DRB3 genotype, respectively.  (Table S4); a total of 31 out of 94 genotypes with frequency > 1% are shown. The X-axis and Y-axis show the genotype name and frequency (%) for each BoLA-DRB3 genotype, respectively.

Association Study of BoLA-DRB3 with PVL
Subsequently, to determine the association between BoLA-DRB3 and BLV PVL, we selected an additional 361 asymptomatic cows, in addition to the original 250 asymptomatic cows (used in Figures 1-3). A total of 611 asymptomatic cows were then divided into the HPVL group (n = 294; Table S6) and LPVL group (n = 317; Table S6). The frequencies of BoLA-DRB3 alleles from LPVL cows and HPVL cows were calculated by Fisher's exact test, and p-values and ORs were estimated for each allele (Table S7). The analysis of allele frequencies ( Figure 4) and association (Table 5)  consistent with previous findings [20]. In addition, DRB3*012:01 (OR = 3.84) was identified as a susceptibility allele.  (Table S6). Total of 12 out of 26 alleles with frequency > 1% are shown.   (Table S6). Total of 12 out of 26 alleles with frequency > 1% are shown.

Differential Susceptibility of BoLA-DRB3 Polymorphisms to Lymphoma and PVL
We compared the effect of BoLA-DRB3 on cow susceptibility to lymphoma and PVL, based on the multivariable logistic regression analysis in Figures 2 and 4. Several different types of BoLA-DRB3 alleles were found to be associated with BLV-induced lymphoma and BLV PVL ( Figure 6). There were two lymphoma resistance alleles, DRB3*010:01, and DRB3*011:01, but no susceptibility alleles were identified. In addition, one allele associated with PVL resistance, DRB3*002:01, and one PVL susceptibility allele, DRB3*012:01, were found. Two resistance alleles, DRB3*009:02 and DRB3*014:01:01, were commonly identified in lymphoma and BLV PVL.
Viruses 2020, 12, x FOR PEER REVIEW 11 of 16 Figure 6. Summary of the differences in BoLA-DRB3 allele-associated proviral load (PVL) and lymphoma susceptibility, based on the logistic regression association study results. R, resistance; S, susceptibility.

BoLA-DRB3 Polymorphisms Are Associated with anti-BLV Antibody Production Levels
Finally, we tried to link the potential biological functions with BoLA-DRB3 polymorphisms. Previously, it has been demonstrated that PVL resistance and susceptibility alleles are associated with anti-BLV antibody production levels [34]. Here, we hypothesized that only PVL-associated alleles, but not lymphoma-associated alleles, would be related to viral antigen-induced immune responses.
To test this, we compared the anti-BLV antibody (anti-gp51) production level between the cows with the PVL susceptibility allele (DRB3*012:01), PVL resistance allele (DRB3*002:01), PVL/lymphoma resistance allele (DRB3*009:02), and lymphoma-specific resistance allele (DRB3*011:01), as shown in the summary in Figure 6. Two-fold serially-diluted plasma samples were tested by ELISA and the OD value in each BoLA-DRB3 group was compared at the dilution of 1:2048. Cows with the PVL resistance allele (PVL resistance group and PVL/lymphoma resistance group) had significantly lower anti-gp51 production levels compared to those in cows carrying the PVL susceptibility allele (p = 0.006 and p = 0.012 respectively; Figure 7). However, cows with the lymphoma specific resistance allele did not have explicitly and significantly different levels of anti-gp51, compared to those in animals with the PVL susceptibility allele (p = 0.474), suggesting that the lymphoma specific associated allele has a lesser effect on anti-gp51 production than the PVL associated allele. Figure 6. Summary of the differences in BoLA-DRB3 allele-associated proviral load (PVL) and lymphoma susceptibility, based on the logistic regression association study results. R, resistance; S, susceptibility.

BoLA-DRB3 Polymorphisms Are Associated with anti-BLV Antibody Production Levels
Finally, we tried to link the potential biological functions with BoLA-DRB3 polymorphisms.
Previously, it has been demonstrated that PVL resistance and susceptibility alleles are associated with anti-BLV antibody production levels [34]. Here, we hypothesized that only PVL-associated alleles, but not lymphoma-associated alleles, would be related to viral antigen-induced immune responses.
To test this, we compared the anti-BLV antibody (anti-gp51) production level between the cows with the PVL susceptibility allele (DRB3*012:01), PVL resistance allele (DRB3*002:01), PVL/lymphoma resistance allele (DRB3*009:02), and lymphoma-specific resistance allele (DRB3*011:01), as shown in the summary in Figure 6. Two-fold serially-diluted plasma samples were tested by ELISA and the OD value in each BoLA-DRB3 group was compared at the dilution of 1:2048. Cows with the PVL resistance allele (PVL resistance group and PVL/lymphoma resistance group) had significantly lower anti-gp51 production levels compared to those in cows carrying the PVL susceptibility allele (p = 0.006 and p = 0.012 respectively; Figure 7). However, cows with the lymphoma specific resistance allele did not have explicitly and significantly different levels of anti-gp51, compared to those in animals with the PVL susceptibility allele (p = 0.474), suggesting that the lymphoma specific associated allele has a lesser effect on anti-gp51 production than the PVL associated allele.

Discussion
In the present study, by using both Fisher's exact test and multivariable logistic regression analysis, we showed for the first time that the susceptibility to BLV-induced lymphoma and PVL is affected by different BoLA-DRB3 polymorphisms. For example, two BoLA-DRB3 alleles, DRB3*010:01, and DRB3*011:01, were found to be associated with resistance to lymphoma but not to PVL. In addition, DRB3*002:01 was specifically associated with PVL resistance. In contrast, we found that DRB3*009:02 was common between lymphoma and PVL resistance, in line with the reciprocal association between PVL and lymphoma development. Thus, we might conclude that host polymorphisms at the BoLA-DRB3 locus are an important factor in both PVL and lymphoma development; interestingly, the PVL-associated BoLA-DRB3 allele did not show a major correlation with the lymphoma-associated BoLA-DRB3 allele. The potential reason for the differential susceptibility to PVL and Lymphoma might be the differential immune response, which depends on BoLA-DRB3 polymorphisms, as we found that the level of anti-BLV antibody is related to the PVL-associated allele but not the lymphoma-associated allele, suggesting that lymphoma specific-associated alleles have a lesser effect on anti-BLV antibody production than PVL-associated alleles.
Some discrepancies were observed between the association based on Fisher's exact test and multivariable logistic regression analysis. For example, DRB3*015:01 was indicated as a lymphoma susceptibility allele by Fisher's exact test, but not by multivariable logistic regression analysis. The inconsistencies can also be found in DRB3*012:01 and DRB3*014:01:01. This is a common problem between these two statistical methods, as Fisher's exact test includes all factors, such as environmental and genetically factors, that together influence BoLA-DRB3 polymorphisms. In contrast, multivariable logistic regression can adjust for the effect of other associated factors, such as age, for the lymphoma association study. Therefore, further experiments for allele functional confirmation are needed.
It is not clear how most BLV-infected cattle do not develop bovine leukosis. The following findings in our study may help solve this issue: (i) In asymptomatic cows, two of the major BoLA-DRB3 alleles, DRB3*011:01 (22%) and DRB3*010:01 (12%), were significantly associated with lymphoma resistance, but were unrelated to PVL. This may explain why some HPVL cows remained asymptomatic. (ii) Susceptibility alleles specifically associated with lymphoma were absent in Japanese Holstein cows, suggesting that malignant transformation requires other factors, besides BoLA-DRB3 polymorphism. For instance, the deregulation of lymphocyte homeostasis is known to lead to leukemia [35]. Provirus integration close to cancer-driver sites and transcriptionally active regions may affect host gene expression [23,36,37]. The viral accessory proteins Tax and G4 also play a crucial role in cell transformation [38]. Besides, p53 mutation [39,40] and tumor necrosis factor-α polymorphisms [41] are also related to lymphoma development.
The major function of MHC class II molecules is to present antigens for T cells to activate the adaptive immune response. It is possible that PVL-specific BoLA-DRB3 and lymphoma-specific BoLA-DRB3 recognize different antigens and thus the subsequent immune response targeting the virus or tumor cells, respectively. To link BoLA-DRB3 polymorphisms with their biological functions, we tested the anti-BLV gp51 antibody level in cows with the PVL-associated BoLA-DRB3 allele and lymphoma-associated BoLA-DRB3 allele. Interestingly, significantly different anti-gp51 levels were found between cows with PVL resistance and PVL susceptibility alleles. However, no significant difference was found in cows carrying lymphoma-specific resistance alleles compared to those in cows carrying PVL susceptibility alleles. This result is in line with the hypothesis that proteins encoded by PVL-associated and lymphoma-associated BoLA-DRB3 alleles recognize different antigens and thus trigger different subsequent immune responses. Similar to that in a previous study, we found that cows with PVL resistance alleles exhibit significantly lower anti-gp51 levels than cows with PVL susceptibility alleles [34]. This is probably because we measured anti-gp51 levels in steady state virus infections, and thus, the antibody concentration might change and correlate with the viral titer in cows. As a result, in PVL resistance cows, which are associated with low viral expression levels, a low anti-BLV antibody level would be detected. In addition to humoral immunity, whether BoLA-DRB3 polymorphisms are associated with effects on cytotoxic T lymphocytes needs further study.
In the PVL association study, DRB3*015:01 has been reported as a PVL susceptibility allele in our and other studies [19,34]. However, we found only one PVL susceptibility allele, DRB3*012:01, in the current investigation. This difference might be due to the sample collection bias and also the statistical analysis method. For the lymphoma association study, it has been reported that the BoLA-DRB3*018:02, DRB3*032:02, and DRB3*009:01 alleles are associated with the susceptibility to BLV-induced lymphoma, whereas DRB3*001:01 and DRB3*011:01 are involved in lymphoma resistance in Iranian Holstein cows [17]. In the present study, DRB3*011:01 was confirmed to be a resistance allele. The other identified alleles were different from those previously reported, suggesting that regional genetic variations may exist. Indeed, ethnicity-related differences in the frequency of human MHC alleles have been observed [42,43]. Furthermore, allelic diversity in the BoLA locus between cattle breeds has been previously demonstrated [44][45][46][47]. This variability is strongly influenced by selective pressures such as exposure to infectious diseases and breed origin. Therefore, the association between the BoLA-DRB3 locus and the resistance or susceptibility to BLV-induced lymphoma, as well as the regulation of PVL, should be further explored in different countries and in distinct cow breeds.
In conclusion, we have demonstrated for the first time that BLV-induced lymphoma and PVL are associated with different BoLA-DRB3 alleles in Holstein cows in Japan. Although BLV infects cattle worldwide, effective treatments and vaccines are not available. Consequently, breed selection based on BoLA-DRB3 polymorphism is a promising strategy to reduce the burden of BLV-induced lymphoma. Contrarily, the sporadic inconsistency between PVL and terminal diseases might be a common phenomenon, due to host genetic polymorphisms during different infectious viral diseases. Indeed, partial inconsistency between PVL and the related symptoms was also observed in HTLV-1-infected patients [48]. As BLV is closely related to HTLV-1, the consistency between the susceptibility of host genetic polymorphisms with PVL and HTLV-1-related symptoms is worth confirming.
Supplementary Materials: The following are available online at http://www.mdpi.com/1999-4915/12/3/352/s1, Table S1. Sample information of asymptomatic cows in BoLA-DRB3 and BLV-induced lymphoma association study. Table S2. Sample information of lymphoma cows in BoLA-DRB3 and BLV-induced lymphoma association study. Table S3. Association study of BoLA-DRB3 alleles between asymptomatic cows and lymphoma cows. Table S4. Results of the stepwise conditional analysis of BoLA-DRB3 alleles in asymptomatic and lymphoma cows. Table S5. Association study of BoLA-DRB3 genotypes between Asymptomatic cows and lymphoma cows. Table S6. Sample information of LPVL and HPVL cows in PVL association study. Table S7. Association study of BoLA-DRB3 alleles between LPVL and HPVL cows. Table S8. Results of the stepwise conditional analysis of BoLA-DRB3 alleles in LPVL and HPVL cows. Table S9. Association study of BoLA-DRB3 genotypes between LPVL and HPVL cows.