Identification of transcriptional biomarkers by RNA-sequencing for improved detection of β(2)-agonists abuse in goat skeletal muscle

In this paper, high-throughput RNA-sequencing (RNA-seq) was used to search for transcriptional biomarkers for β(2)-agonists. In combination with drug mechanisms, a smaller group of genes with higher detection accuracy was screened out. Unknown samples were first predicted by this group of genes, and liquid chromatograph tandem mass spectrometer (LC-MS/MS) was applied to positive samples to validate the biomarkers. The results of principal component analysis (PCA), hierarchical cluster analysis (HCA) and discriminant analysis (DA) indicated that the eight genes screened by high-throughput RNA-seq were able to distinguish samples in the experimental group and control group. Compared with the nine genes selected from an earlier literature, 17 genes including these nine genes were proven to have a more satisfactory effect, which validated the accuracy of gene selection by RNA-seq. Then, six key genes were selected from the 17 genes according to the variable importance in projection (VIP) value of greater than 1. The test results using the six genes and 17 genes were similar, revealing that the six genes were critical genes. By using the six genes, three positive samples possibly treated with drugs were screened out from 25 unknown samples through DA and partial least squares discriminant analysis (PLS-DA). Then, the three samples were verified by a standard method, and mapenterol was detected in a sample. Therefore, the six genes can be used as biomarkers to detect β(2)-agonists. Compared with the previous study, accurate detection of β(2)-agonists abuse using six key genes is an improvement method, which show great significance in the monitoring of β(2)-agonists abuse in animal husbandry.