Random forest algorithm reveals novel sites in HA protein that shift receptor binding preference of the H9N2 avian influenza virus

After filtering out redundant, incomplete, environmental source, and unclear information sequences, a total of 5,656 H9N2 HA gene sequences were obtained from the National Center for Biotechnology Information (NCBI), Global Initiative on Sharing All Influenza Data (GISAID), and Influenza Research Database (IRD) databases. HA sequences were divided into two datasets based on host information: avian-derived sequences (5,588) and mammal-derived sequences (68) and labeled accordingly.

Amino acid types were replaced with numerical values to transform the amino acid sequences into machine-readable vectors (Supplementary Table S1). We used the Random Forest Classifier method in Sklearn (1.0.22) to train the random forest classifier (Abraham et al., 2014), and in consideration of the large difference in the numbers of avian and non-avian sequences, we selected a balanced number of samples for training.