Studying the genetic loci related to the bovine leukemia virus using random forest method and genomic data

نویسندگان

1 Department of Animal Science, Faculty of Agricultural Science and Engineering, University College of Agriculture and Natural Resources (UTCAN), University of Tehran, Karaj, Iran

2 Department of Animal Science, Faculty of Agriculture, Tarbiat Modares University, Tehran, Iran

3 Department of Animal Science, Faculty of Agricultural Science and Engineering, University College of Agriculture and Natural Resources (UTCAN), University of Tehran, Karaj, Iran

4 Department of Animal Science, Faculty of Agricultural Science and Engineering, University College of Agriculture and Natural Resources (UTCAN), University of Tehran, Karaj, Iran

doi
10.22103/jlst.2025.25353.1634
چکیده

Bovine leukemia virus (BLV) is a causative agent of bovine leukosis, which, due to its long incubation period, can spread widely within a herd before clinical symptoms appear, causing significant economic losses. This study used a supervised machine learning method called random forest to identify genomic regions associated with BLV. The non-parametric nature of this method allows for the creation of predictive models without the need for initial statistical assumptions; whereas the standard Genome-wide association studies (GWAS) methods are usually based on single-variable hypothesis tests and cannot account for correlations resulting from connectivity imbalance or the combination of multiple markers. In this study, the genotyping data of 145 Holstein cows (77 BLV-positive, 68 healthy) after quality control by using the PLINK (v 1.02), which resulted in 23,910 Single nucleotide polymorphisms (SNPs) were analyzed. Random forest analyses on the mentioned data included three hyperparameters: mtry (0.5(p/3), (p/3), 2(p/3)), ntree (2000, 3000, 4000), and nodesize (5, 10, 15), where p is equal to the total number of SNPs (23,910). To find the best SNPs, the Mean Decrease Accuracy (MDA) index (> 1.89) was used which resulted in the selection of 50 SNPs. Genomic enrichment analyses showed that genes associated with the top 50 SNPs are predominantly involved in Positive Regulation, Intracellular Signaling, Apoptosis and Cell Death, Signal Transduction, Metabolic Processes, and Cell Differentiation and Development. In total, 82 genes were identified, including hub genes such as MYC, RABIF, IRS1, TRAPPC9, MAPK8, HTT, SNX9, BCLAF1, XRN1, and LSM6.