Genome-wide scan for commons SNPs affecting bovine leukemia virus infection level in dairy cattle

Direct genetic effects, maternal genetic effects, and maternal genetic sensitivity on prenatal heat stress for calf diseases and corresponding genomic loci in German Holsteins

The aim of this study was to infer the effects of heat stress (HS) of dams during late gestation on direct and maternal genetic parameters for pneumonia (PNEU, 112,563 observations), diarrhea (DIAR, 176,904 observations), and omphalitis (OMPH, 176,872 observations) in Holstein calves kept in large-scale co-operator herds. The genotype dataset included 41,135 SNPs from 19,247 male and female cattle. Temperature-humidity indices (THI) during the last 8 wk of pregnancy were calculated, using the climate data from the nearest public weather station for each herd. Heat load effects were considered for average weekly THI larger than 60. Phenotypically, regression coefficients of calf diseases on prenatal THI during the last 8 wk of gestation were estimated in 8 consecutive runs. The strongest detrimental effects of prenatal HS on PNEU and DIAR were identified for the last week of pregnancy (wk 1). Thus, only wk 1 was considered in ongoing genetic and genomic analyses. In an advanced model considering prenatal HS, random regression coefficients on THI in wk 1 nested within maternal genetic effects (maternal slope effects for heat load) were considered as parameters to infer maternal sensitivity in response to prenatal THI alterations. Direct heritabilities from the advanced model ranged from 0.10 (THI 60) to 0.08 (THI 74) for PNEU and were close to 0.16 for DIAR. Maternal heritabilities for PNEU increased from 0.03 to 0.10 along the THI gradient. For DIAR, the maternal heritability was largest (0.07) at the minimum THI (THI = 60) and decreased to 0.05 at THI 74. Genetic correlations smaller than 0.80 for PNEU and DIAR recorded at THI 60 with corresponding diseases at THI 74 indicated genotype by climate interactions for maternal genetic effects. Genome-wide associations studies were performed using de-regressed proofs of genotyped sires for direct genetic, maternal genetic, and maternal slope effects. Thirty suggestive and 2 significant SNPs were identified from the GWAS. Forty-three genes located close to the suggestive SNPs (±100 kb) were annotated as potential candidate genes. Three biological processes were inferred on the basis of the these genes, addressing the negative regulation of the viral life cycle, innate immune response, and protein ubiquitination. Hence, the genetics of prenatal heat stress mechanisms are associated with immune physiology and disease resistance mechanisms.

Single-Nucleotide Polymorphism on Spermatogenesis Associated 16 Gene-Coding Region Affecting Bovine Leukemia Virus Proviral Load

Bovine leukemia virus (BLV) is an etiological agent of malignant lymphoma in cattle and is endemic in many cattle-breeding countries. Thus, the development of cattle genetically resistant to BLV is desirable. The purpose of this study was to identify novel single-nucleotide polymorphisms (SNPs) related to resistance to BLV. A total of 146 DNA samples from cattle with high BLV proviral loads (PVLs) and 142 samples from cattle with low PVLs were used for a genome-wide association study (GWAS). For the verification of the GWAS results, an additional 1342 and 456 DNA samples from BLV-infected Japanese Black and Holstein cattle, respectively, were used for an SNP genotyping PCR to compare the genotypes for the identified SNPs and PVLs. An SNP located on the spermatogenesis associated 16 (SPATA16)-coding region on bovine chromosome 1 was found to exceed the moderate threshold (p < 1.0 × 10−5) in the Additive and Dominant models of the GWAS. The SNP genotyping PCR revealed that the median values of the PVL were 1278 copies/50 ng of genomic DNA for the major homozygous, 843 for the heterozygous, and 621 for the minor homozygous genotypes in the Japanese Black cattle (p < 0.0001). A similar tendency was also observed in the Holstein cattle. We found that cattle with the minor allele for this SNP showed 20−25% lower PVLs. Although the mechanisms through which this SNP impacts the PVL remain unknown, we found a novel SNP related to BLV resistance located on the SPATA16 gene-coding region on bovine chromosome 1.

Genetic Influence on Frequencies of Myeloid-Derived Cell Subpopulations in Mouse

Differences in frequencies of blood cell subpopulations were reported to influence the course of infections, atopic and autoimmune diseases, and cancer. We have discovered a unique mouse strain B10.O20 containing extremely high frequency of myeloid-derived cells (MDC) in spleen. B10.O20 carries 3.6% of genes of the strain O20 on the C57BL/10 genetic background. It contains much higher frequency of CD11b⁺Gr1⁺ cells in spleen than both its parents. B10.O20 carries O20-derived segments on chromosomes 1, 15, 17, and 18. Their linkage with frequencies of blood cell subpopulations in spleen was tested in F₂ hybrids between B10.O20 and C57BL/10. We found 3 novel loci controlling MDC frequencies: Mydc1, 2, and 3 on chromosomes 1, 15, and 17, respectively, and a locus controlling relative spleen weight (Rsw1) that co-localizes with Mydc3 and also influences proportion of white and red pulp in spleen. Mydc1 controls numbers of CD11b⁺Gr1⁺ cells. Interaction of Mydc2 and Mydc3 regulates frequency of CD11b⁺Gr1⁺ cells and neutrophils (Gr1⁺Siglec-F^- cells from CD11b⁺ cells). Interestingly, Mydc3/Rsw1 is orthologous with human segment 6q21 that was shown previously to determine counts of white blood cells. Bioinformatics analysis of genomic sequence of the chromosomal segments bearing these loci revealed polymorphisms between O20 and C57BL/10 that change RNA stability and genes’ functions, and we examined expression of relevant genes. This identified potential candidate genes Smap1, Vps52, Tnxb, and Rab44. Definition of genetic control of MDC can help to personalize therapy of diseases influenced by these cells.

Taurine and Indicine Haplotype Representation in Advanced Generation Individuals From Three American Breeds

Development of the American Breeds of beef cattle began in the 1920s as breeders and U. S. Experiment Station researchers began to create Bos taurus taurus × Bos taurus indicus hybrids using Brahman as the B. t. indicus source. By 1954, U.S. Breed Associations had been formed for Brangus (5/8 Angus × 3/8 Brahman), Beefmaster (½ Brahman × ¼ Shorthorn × ¼ Hereford), and Santa Gertrudis (5/8 Shorthorn × 3/8 Brahman). While these breeds were developed using mating designs expected to create base generation animals with the required genome contributions from progenitor breeds, each association has now registered advanced generation animals in which selection or drift may have caused the realized genome compositions to differ from initial expected proportions. The availability of high-density SNP genotypes for 9,161 Brangus, 3,762 Beefmaster, and 1,942 Santa Gertrudis animals allowed us to compare the realized genomic architectures of breed members to the base generation expectations. We used RFMix to estimate local ancestry and identify genomic regions in which the proportion of Brahman ancestry differed significantly from a priori expectations. For all three breeds, lower than expected levels of Brahman composition were found genome-wide, particularly in early-generation animals where we demonstrate that selection on beef production traits was likely responsible for the taurine enrichment. Using a proxy for generation number, we also contrasted the genomes of early- and advanced-generation animals and found that the indicine composition of the genome has increased with generation number likely due to selection on adaptive traits. Many of the most-highly differentiated genomic regions were breed specific, suggesting that differences in breeding objectives and selection intensities exist between the breeds. Global ancestry estimation is commonly performed in admixed animals to control for stratification in association studies. However, local ancestry estimation provides the opportunity to investigate the evolution of specific chromosomal segments and estimate haplotype effects on trait variation in admixed individuals. Investigating the genomic architecture of the American Breeds not only allows the estimation of indicine and taurine genome proportions genome-wide, but also the locations within the genome where either taurine or indicine alleles confer a selective advantage.

Genome-wide selection of discriminant SNP markers for breed assignment in indigenous sheep breeds

The assignment of an individual to the true population of origin is one of the most important applications of genomic data for practical use in animal breeding. The aim of this study was to develop a statistical method and then, to identify the minimum number of informative SNP markers from high-throughput genotyping data that would be able to trace the true breed of unknown samples in indigenous sheep breeds. The total numbers of 217 animals were genotyped using Illumina OvineSNP50K BeadChip in Zel, Lori-Bakhtiari, Afshari, Moqani, Qezel and a wild-type Iranian sheep breed. After SNP quality check, the principal component analysis (PCA) was used to determine how the animals allocated to the groups using all genotyped markers. The results revealed that the first principal component (PC1) separated out the two domestic and wild sheep breeds, and all domestic breeds were separated from each other for PC2. The genetic distance between different breeds was calculated using FST and Reynold methods and the results showed that the breeds were well differentiated. A statistical method was developed using the stepwise discriminant analysis (SDA) and the linear discriminant analysis (LDA) to reduce the number of SNPs for discriminating 6 different Iranian sheep populations and K-fold cross-validation technique was employed to evaluate the potential of a selected subset of SNPs in assignment success rate. The procedure selected reduced pools of markers into 201 SNPs that were able to exactly discriminate all sheep populations with 100% accuracy. Moreover, a discriminate analysis of principal components (DAPC) developed using 201 linearly independent SNPs revealed that these markers were able to assign all individuals into true breed. Finally, these 201 identified SNPs were successfully used in an independent out-group breed consisting of 96 samples of Baluchi sheep breed and the results indicated that these markers are able to correctly allocate all unknown samples to true population of origin. In general, the results of this study indicated that the combined use of the SDA and LDA techniques represents an efficient strategy for selecting a reduced pool of highly discriminant markers.

Exonic SNP in MHC-DMB2 is associated with gene expression and humoral immunity in Japanese quails

The DMB2 gene is widely expressed at high levels in avian. This gene plays an important role in humoral immunity. The aim of this study was to investigate the effects of 361 G > C Single nucleotide polymorphism (SNP) on DMB2 protein structure and gene expression to determine how the 361 G > C SNP affects humoral immune response in Japanese quails. 0.2 mL of 5% sheep red blood cell (SRBC) was injected into breast muscle of 130 Japanese quails on 28 days. After DNA extraction, PCR was carried out to amplify a 333-base pair DNA fragment from the exon 2 of DMB2 gene. The pattern of all samples was determined through RFLP technique. PCR-RFLP results identified two alleles segregating (C, G) as three genotypes (CC, CG and GG) in Japanese Quails. The antibody response to SRBC with CC genotype was significantly higher than the CG and GG genotypes (P < 0.01). In silico analysis showed that the 361 G > C SNP has no effect on the physicochemical properties and 3D structure. The results of RT-qPCR indicated that the effect of genotype on gene expression is significant, so that the expression of CC genotype is more than CG and GG genotype. It can be inferred that the 361 G > C SNP in the exon 2 of MHC-DMB2 gene is not desirable. This mutation decreases humoral immune response by reducing DMB2 gene expression.

Specific antiviral effect of violaceoid E on bovine leukemia virus

Bovine leukemia virus (BLV) infection has spread worldwide causing significant economic losses in the livestock industry. In countries with a high prevalence of BLV, minimizing economic losses is challenging; thus, research into various countermeasures is important for improving BLV control. Because anti-BLV drugs have not been developed, the present study explored a promising chemical compound with anti-BLV activity. Initially, screening of a chemical compound library revealed that violaceoid E (vioE), which is isolated from fungus, showed antiviral activity. Further analysis demonstrated that the antiviral effect of vioE inhibited transcriptional activation of BLV. Cellular thermal shift assay and pulldown assays provided evidence for a direct interaction between vioE and the viral transactivator protein, Tax. These data indicate that interference with Tax-dependent transcription could be a novel target for development of anti-BLV drugs. Therefore, it is suggested that vioE is a novel antiviral compound against BLV.

Genomic inbreeding and selection signatures in the local dairy breed Icelandic Cattle

Icelandic Cattle is the only dairy cattle breed native to Iceland. It currently numbers ca. 26 000 breeding females. We used 50k genotypes of over 8000 Icelandic Cattle to estimate genomic and pedigree-based inbreeding and to detect selection signatures using the integrated haplotype score. We used 47 Icelandic bulls genotyped with a 770k SNP chip to estimate LD decay for comparison with other Nordic dairy cattle breeds. We detected ROHs on the autosomes and computed ROH-based autosomal inbreeding coefficients. Average inbreeding coefficients according to pedigree and ROHs were 0.0621 and 0.101. Effective population sizes for the years 2009-2017 according to pedigree, ROHs, genomic relationship matrix, excess of homozygosity and individual increase in inbreeding were 81, 65, 60, 58 and 92 respectively. We identified three regions and six candidate genes that showed a signature of selection according to the integrated haplotype score (P < 0.05) on chromosomes 1, 16 and 23. The LD structure of Icelandic Cattle is shaped by a long period of isolation and a small founder population. The estimate of LD at distances closer than 0.3 Mb is lower in Icelandic Cattle than in Danish Jersey, but is higher than in Danish Holstein and Red Nordic Dairy Cattle breeds. Our findings show that inbreeding rates in Icelandic Cattle currently are sustainable according to FAO guidelines, and our results do not indicate severe historical inbreeding.

Expression-based analysis of genes related to single nucleotide polymorphism hits associated with bovine leukemia virus proviral load in Argentinean dairy cattle

In dairy cattle infected with bovine leukemia virus (BLV), the proviral load (PVL) level is directly related to the viral transmission from infected animals to their healthy herdmates. Two contrasting phenotypic groups can be identified when assessing PVL in peripheral blood of infected cows. A large number of reports point to bovine genetic variants (single nucleotide polymorphisms) as one of the key determinants underlying PVL level. However, biological mechanisms driving BLV PVL profiles and infection progression in cattle have not yet been elucidated. In this study, we evaluated whether a set of candidate genes affecting BLV PVL level according to whole genome association studies are differentially expressed in peripheral blood mononuclear cells derived from phenotypically contrasting groups of BLV-infected cows. During a 10-mo-long sampling scheme, 129 Holstein cows were phenotyped measuring anti-BLV antibody levels, PVL quantification, and white blood cell subpopulation counts. Finally, the expression of 8 genes (BOLA-DRB3, PRRC2A, ABT1, TNF, BAG6, BOLA-A, LY6G5B, and IER3) located within the bovine major histocompatibility complex region harboring whole genome association SNP hits was evaluated in 2 phenotypic groups: high PVL (n = 7) and low PVL (n = 8). The log₂ initial fluorescence value (N₀) transformed mean expression values for the ABT1 transcription factor were statistically different in high- and low-PVL groups, showing a higher expression of the ABT1 gene in low-PVL cows. The PRRC2A and IER3 genes had a significant positive (correlation coefficient = 0.61) and negative (correlation coefficient = -0.45) correlation with the lymphocyte counts, respectively. Additionally, the relationships between gene expression values and lymphocyte counts were modeled using linear regressions. Lymphocyte levels in infected cows were better explained (coefficient of determination = 0.56) when fitted a multiple linear regression model using both PRRC2A and IER3 expression values as independent variables. The present study showed evidence of differential gene expression between contrasting BLV infection phenotypes. These genes have not been previously related to BLV pathobiology. This valuable information represents a step forward in understanding the BLV biology and the immune response of naturally infected cows under a commercial milk production system. Efforts to elucidate biological mechanisms leading to BLV infection progression in cows are valuable for BLV control programs. Further studies integrating genotypic data, global transcriptome analysis, and BLV progression phenotypes are needed to better understand the BLV-host interaction.

Risk factor analysis of bovine leukemia virus infection in dairy cattle in Egypt

Identification of the risk factors associated with Enzootic bovine leukosis (EBL) is essential for the adoption of potentially prevention strategies. Accordingly, our objectives were to determine the geographic distribution of Bovine Leukemia Virus (BLV) infection and identify the risk factors associated with cow-level BLV infection in the Egyptian dairy cattle. A cross-sectional study was conducted on 1299 mixed breed cows distributed over four provinces in the Nile Delta of Egypt in 2018. The randomly selected cows on each farm were serologically tested for BLV, and the cow's information was obtained from the farm records. Four variables (geographic location, herd size, number of parities, and age) were used for risk analysis. A total of 230 serum samples (17.7 %) were serologically positive for BLV. The highest prevalence of BLV infection was associated with parity (OR = 3.4, 95 %CI 2.4-4.9) with 80 % probability of being BLV-positive at parity ≥5, followed by herd size (OR = 1.8, 95 %CI 1.4-2.2). However, geographic location seems to have no impact on the prevalence of BLV infection in Egypt. Our findings strongly indicate that the intensive surveillance and effective prevention strategies against BLV infection in Egypt should be provided to multiparous cows with ≥5 parities and live in large farm with more than 200 cows.

PRMT5 Is Required for Bovine Leukemia Virus Infection In Vivo and Regulates BLV Gene Expression, Syncytium Formation, and Glycosylation In Vitro

Bovine leukemia virus (BLV) is the causative agent of enzootic bovine leukosis, which is the most common neoplastic disease of cattle and is closely related to human T-cell leukemia viruses. We investigated the role of a new host protein, PRMT5, in BLV infection. We found that PRMT5 is overexpressed only in BLV-infected cattle with a high proviral load, but not in those with a low proviral load. Furthermore, this upregulation continued to the lymphoma stage. PRMT5 expression was upregulated in response to experimental BLV infection; moreover, PRMT5 upregulation began in an early stage of BLV infection rather than after a long period of proviral latency. Second, siRNA-mediated PRMT5 knockdown enhanced BLV gene expression at the transcript and protein levels. Additionally, a selective small-molecule inhibitor of PRMT5 (CMP5) enhanced BLV gene expression. Interestingly, CMP5 treatment, but not siRNA knockdown, altered the gp51 glycosylation pattern and increased the molecular weight of gp51, thereby decreasing BLV-induced syncytium formation. This was supported by the observation that CMP5 treatment enhanced the formation of the complex type of N-glycan more than the high mannose type. In conclusion, PRMT5 overexpression is related to the development of BLV infection with a high proviral load and lymphoma stage and PRMT5 inhibition enhances BLV gene expression. This is the first study to investigate the role of PRMT5 in BLV infection in vivo and in vitro and to reveal a novel function for a small-molecule compound in BLV-gp51 glycosylation processing.

Genetic variant of SPARC gene and its association with growth traits in Chinese cattle

SPARC is a cysteine-rich acidic secreted protein. It is a non-collagen component of bone, which is widely distributed in humans and animals and plays an important role. SPARC has been found in a variety of human cancers (breast cancer, stomach cancer, ovarian cancer, etc.) and diabetes-related research. Especially the muscle and fat metabolism are closely related. In this study, we used a DNA pool to detect a new SNP site (g.12454T  >  C). A total of 616 samples of four breeds of Qinchuan cattle (QC, n = 176 ), Xianan cattle (XN, n = 160 ), Pinan cattle (PN, n = 136 ) and Jiaxian cattle (JX, n = 144 ) were analyzed and identified with ARMS-PCR. In addition, we correlated SNP with growth traits and showed significant correlation with growth traits such as rump length, hip width, and body length ( p < 0.05 ). Moreover, we tested the SPARC gene expression level in different tissues belonging to XN adult cattle ( n = 3 ) and found its high expression in muscle tissues (relative to the kidney). Further, we found the SNP is able to increase the SPARC expression level in skeletal muscle ( n = 12 ). According to statistical data, this SNP site may be applied to a molecular marker of an early marker-assisted selection for early growth of beef 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.

Single nucleotide polymorphism rs110861313 in the intergenic region of chromosome 23 is associated with the development of leukosis in the Russian Black Pied cattle

In recent years, using a genome-wide association study (GWAS), a number of single nucleotide polymorphisms (SNPs) have been suggested to be associated with susceptibility to leukemia in cattle. However, all studies have been done with purebred Holstein cows and their hybrids. In this regard, it is important to confirm the functional role of polymorphisms previously identified in a GWAS study in Russian cattle breeds. The aim of this study was to verify the association between rs110861313 in the intergenic region of bovine chromosome 23 and leukemia in the Russian Black Pied cattle. Based on the levels of bovine leukemia virus (BLV)-specific antibodies detected in serum using serodiagnostic techniques, animals were divided into three groups: healthy animals (n = 115), asymptomatic virus carriers (n = 145) and animals with leukemia (n = 107). Genotyping of rs110861313 was carried out using polymerase chain reaction followed by analysis of restriction fragment length polymorphisms. A significant decrease in the frequency of the A/A genotype (11.2 %) was revealed in animals with persistent lymphocytosis compared to virus carriers (27.6 %) (p < 0.002). At the same time, the frequency of animals with the C/C genotype in animals with persistent lymphocytosis (41.1 %) was significantly higher than that of virus carriers (21.4 %) (p < 0.001). In this case, asymptomatic virus carriers can be considered a more suitable control than healthy animals that have not been in contact with the virus. According to bioinformatics analysis, resistance to BLV can be due to the presence of the transcription factor FOXM1 binding site in the region of rs110861313. FOXM1 is expressed in immune cells and can potentially affect the expression of the neighboring genes (LY6G5B, GPANK1, ABHD16A, LY6G6F, LY6G6E, CSNK2B, ApoM). Thus, we found that SNP rs110861313 in the intergenic region of bovine chromosome 23 is associated with the development of leukemia following BLV infection in the Russian Black Pied cattle.

A point mutation to the long terminal repeat of bovine leukemia virus related to viral productivity and transmissibility

It is important to establish the molecular basis of the high transmissibility of bovine leukemia virus (BLV) to develop new methods of preventing viral transmission. Hence, the aim of this study was to determine whether some strains had transmission advantages. First, we determined the whole BLV genome sequences of all 34 BLV-infected cows from one farm. Phylogenetic analysis divided strains into 26 major and 8 minor strains. The major strains dominantly spread independent of host factor, bovine leucocyte antigen. Further analysis, with molecular clones, associated transmissibility with viral productivity in vitro. In addition, the two groups could be classified by group-specific mutations. The reverse genetic approach demonstrated that a spontaneous mutation at nucleotide 175 of the BLV genome, which is located in the viral promoter region, could alter viral productivity by changing viral transactivation, suggesting that BLV transmissibility is affected by a spontaneous mutation associated with viral productivity.

Bovine leukemia virus proviral load is more strongly associated with bovine major histocompatibility complex class II DRB3 polymorphism than with DQA1 polymorphism in Holstein cow in Japan

Bovine leukemia virus (BLV) causes enzootic bovine leukosis and is closely related to the human T-lymphotropic virus. Bovine major histocompatibility complex (BoLAs) are used extensively as markers of disease and immunological traits in cattle. For BLV diagnosis, proviral load is a major diagnosis index for the determination of disease progression and transmission risk. Therefore, we investigated the frequency of BoLA-DRB3 alleles, BoLA-DQA1 alleles, and haplotypes of BoLA class II isolated from the heads of 910 BLV-infected cows out of 1290 cows assessed from BLV-positive farms, in a nationwide survey from 2011 to 2014 in Japan. Our aim was to identify BoLA class II polymorphisms associated with the BLV proviral load in the Holstein cow. The study examined 569 cows with a high proviral load and 341 cows with a low proviral load. Using the highest odds ratio (OR) as a comparison index, we confirmed that BoLA-DRB3 was the best marker for determining which cow spread the BLV (OR 13.9 for BoLA-DRB3, OR 11.5 for BoLA-DQA1, and OR 6.2 for BoLA class II haplotype). In addition, DRB3*002:01, *009:02, *012:01, *014:01, and *015:01 were determined as BLV provirus associated alleles. BoLA-DRB3*002:01, *009:02, and *014:01 were determined as resistant alleles (OR > 1), and BoLA-DRB3*012:01 and *015:01 were determined as susceptible alleles (OR < 1). In this study, we showed that BoLA-DRB3 was a good marker for determining which cow spread BLV, and we found not only one resistant allele (BoLA-DRB3*009:02), but also two other disease-resistant alleles and two disease-susceptible alleles. This designation of major alleles as markers of susceptibility or resistance can allow the determination of the susceptibility or resistance of most cows to disease. Overall, the results of this study may be useful in eliminating BLV from farms without having to separate cows into several cowsheds.

Prioritization of genes associated with the pathogenesis of leukosis in cattle

Selection by means of genetic markers is a promising approach to the eradication of infectious diseases in farm animals, especially in the absence of effective methods of treatment and prevention. Bovine leukemia virus (BLV) is spread throughout the world and represents one of the biggest problems for the livestock production and food security in Russia. However, recent genome-wide association studies have shown that sensitivity/resistance to BLV is polygenic. The aim of this study was to create a catalog of cattle genes and genes of other mammalian species involved in the pathogenesis of BLV-induced infection and to perform gene prioritization using bioinformatics methods. Based on manually collected information from a range of open sources, a total of 446 genes were included in the catalog of cattle genes and genes of other mammals involved in the pathogenesis of BLV-induced infection. The following criteria were used to prioritize 446 genes from the catalog: (1) the gene is associated with leukemia according to a genome-wide association study; (2) the gene is associated with leukemia according to a case-control study; (3) the role of the gene in leukemia development has been studied using knockout mice; (4) protein-protein interactions exist between the gene-encoded protein and either viral particles or individual viral proteins; (5) the gene is annotated with Gene Ontology terms that are overrepresented for a given list of genes; (6) the gene participates in biological pathways from the KEGG or REACTOME databases, which are over-represented for a given list of genes; (7) the protein encoded by the gene has a high number of protein-protein interactions with proteins encoded by other genes from the catalog. Based on each criterion, a rank was assigned to each gene. Then the ranks were summarized and an overall rank was determined. Prioritization of 446 candidate genes allowed us to identify 5 genes of interest (TNF,LTB,BOLA-DQA1,BOLA-DRB3,ATF2), which can affect the sensitivity/resistance of cattle to leukemia.