Bioinformatics analysis of structure and function in the MRP gene family and its expression in response to various drugs in Bursaphelenchus xylophilus

Abstract

Genes homologous to members of the MRP gene family in Caenorhabditis elegans are important in drug resistance. To further explore the molecular mechanism of drug resistance in pine wood nematode (Bursaphelenchus xylophilus), we used bioinformatics approaches to analyze genomic data for B. xylophilus and identified Bx-MRP genes. We predicted the structure and function of the genes and encoded proteins. Using bioinformatics programs to predict and analyze various properties of the predicted proteins, including hydrophobicity, transmembrane regions, phosphorylation sites, and topologically isomeric structures, of these Bx-MRP genes, we determined that they function in transmembrane transport. From the results of RT-qPCR, the Bx-MRP family members confer significant differential resistance to different drug treatments. After treatment with different concentrations of emamectin benzoate, avermectin and matrine, the expression of each gene increased with increasing drug concentrations, indicating that the family members play a positive role in the regulation of multidrug resistance.

Keywords

Bursaphelenchus xylophilus / Bx-MRP gene family / Gene and protein structure / Bioinformatics / Multidrug stress

References

[1] Bailey TL, Johnson J, Grant CE, Noble WS. The MEME suite. Nucleic Acids Res, 2015, 43: W39-W49.

[2] Crooks GE, Hon G, Chandonia JM, Brenner SE. WebLogo: a sequence logo generator. Genome Res, 2004, 14: 1188-1190.

[3] Dermauw W, Van Leeuwen T. The ABC gene family in arthropods: comparative genomics and role in insecticide transport and resistance. Insect Biochem Mol, 2014, 45: 89-110.

[4] Domenichini A, Adamska A, Falasca M. ABC transporters as cancer drivers: potential functions in cancer development. Biochim Biophys Acta General Subjects, 2019, 1863: 52-60.

[5] Finn RD, Coggill P, Eberhardt RY, Eddy SR, Mistry J, Mitchell AL, Potter SC, Punta M, Qureshi M, Sangrador-Vegas A, Salazar GA, Tate J, Bateman A. The Pfam protein families database: towards a more sustainable future. Nucleic Acids Res, 2016, 44: D279-D285.

[6] Futai K. Pine Wood Nematode, Bursaphelenchus xylophilus. Annu Rev Phytopathol, 2013, 51: 61-83.

[7] Gameiro M, Silva R, Rocha-Pereira C, Carmo H, Carvalho F, Bastos MDL, Remiao F. Cellular models and in vitro assays for the screening of modulators of p-gp, MRP1 and BCRP. Molecules, 2017 22 4 600

[8] Hao X, Wang F, Ma L, Wang BW, Ma Y. Cloning and function of resistance gene of Bursaphelenchus xylophilus. J. North East For. Univ., 2018 46 89–92 97

[9] He ZL, Zhang HK, Gao SH, Lercher MJ, Chen WH, Hu SN. Evolview v2: an online visualization and management tool for customized and annotated phylogenetic trees. Nucleic Acids Res, 2016, 44: W236-W241.

[10] Hyde SC, Emsley P, Hartshorn MJ, Mimmack MM, Gileadi U, Pearce SR, Gallagher MP, Gill DR, Hubbard RE, Higgins CF. Structural model of ATP-binding proteins associated with cystic fibrosis, multidrug resistance and bacterial transport. Nature, 1990, 346: 362-365.

[11] Kikuchi T, Cotton JA, Dalzell JJ, Hasegawa K, Kanzaki N, McVeigh P, Takanashi T, Tsai IJ, Assefa SA, Cock PJA, Otto TD, Hunt M, Reid AJ, Sanchez-Flores A, Tsuchihara K, Yokoi T, Larsson MC, Miwa J, Maule AG, Sahashi N, Jones JT, Berriman M. Genomic insights into the origin of parasitism in the emerging plant pathogen Bursaphelenchus xylophilus. PLoS Pathog, 2011 7 9 e1002219

[12] Krogh A, Larsson B, von Heijne G, Sonnhammer EL. Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. J Mol Biol, 2001, 305: 567-580.

[13] Kwon T. Change of abundance of arthropods in pine forests caused by aerial insecticide spray. Arch Environ Contam Toxicol, 2008, 54: 92-106.

[14] Li YX, Meng FL, Deng X, Wang X, Feng YQ, Zhang W, Pan L, Zhang XY. Comparative transcriptome analysis of the pinewood nematode Bursaphelenchus xylophilus reveals the molecular mechanism underlying its defense response to host-derived-pinene. Int J Mol Sci, 2019 20 4 911

[15] Li YL, Zheng CY, Liu KC, Wu Y, Fan B, Han ZM. Transformation of multi-antibiotic resistant Stenotrophomonas maltophilia with GFP Gene to enable tracking its survival on pine trees. Forests, 2019 10 3 231

[16] Liu HB, Wu F, Wu XQ, Ye JR. Differential effects of rapamycin on Bursaphelenchus xylophilus with different virulence and differential expression of autophagy genes under stresses in nematodes. Acta Biochim Biophys Sin, 2019, 51: 254-262.

[17] Liu MJ, Hwang B, Jin CZ, Li WJ, Park DJ, Seo S, Kim CJ. Screening, isolation and evaluation of a nematicidal compound from actinomycetes against the pine wood nematode, Bursaphelenchus xylophilus. Pest Manag Sci, 2019, 75: 1585-1593.

[18] Ma YF, Wu ZH, Gao M, Loor JJ. Nuclear factor erythroid 2-related factor 2-antioxidant activation through the action of ataxia telangiectasia-mutated serine/threonine kinase is essential to counteract oxidative stress in bovine mammary epithelial cells. J Dairy Sci, 2018, 101: 5317-5328.

[19] Mota MM, Vieira P. Pine wilt disease: a worldwide threat to forest ecosystems, 2008, Berlin: Springer

[20] Omasits U, Ahrens CH, Mueller S, Wollscheid B. Protter: interactive protein feature visualization and integration with experimental proteomic data. Bioinformatics, 2014, 30: 884-886.

[21] Seddigh S. Comprehensive comparison of two protein family of P-ATPases (13A1 and 13A3) in insects. Comput Biol Chem, 2017, 68: 266-281.

[22] Seddigh S, Darabi M. Functional, structural, and phylogenetic analysis of mitochondrial cytochrome b (cytb) in insects. Mitochondrial DNA A, 2018, 29: 236-249.

[23] Soliman T, Mourits MCM, van der Werf W, Hengeveld GM, Robinet C, Lansink AGJM. Framework for modelling economic impacts of invasive species, applied to pine wood nematode in Europe. PLOS ONE, 2012 7 9 e45505

[24] Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol, 2011, 28: 2731-2739.

[25] Wang XR, Cheng X, Li YD, Zhang JA, Zhang ZF, Wu HR. Cloning arginine kinase gene and its RNAi in Bursaphelenchus xylophilus causing pine wilt disease. Eur J Plant Pathol, 2012, 134(3): 521-532.

[26] Wang F, Li DL, Ma L, Wang BW, Chen QL, Zhang RZ, Su D, Kang XY, Zhai W. Identification and RNAi of a Bursaphelenchus xylophilus dauer formation gene: Bx-DAF6. For Stud China, 2016, 38: 21-27.

[27] Wang BW, Ma L, Wang F, Wang BY, Hao X, Xu JY, Ma Y. Low temperature extends the lifespan of Bursaphelenchus xylophilus throuth the cGMP pathway. Int J Mol Sci, 2017 18 11 2320

[28] Wang BW, Hao X, Xu JY, Ma Y, Ma L. Transcriptome-Based analysis reveals a crucial role of BxGPCR17454 in low temperature response of pine wood nematode (Bursaphelenchus xylophilus). Int J Mol Sci, 2019 20 12 2898

[29] Wang BY, Wen RR, Ma L, Ma Y, Zhou TH. Molecular characterization and functional analysis of daf-16-2b gene in the pine wood nematode, Bursaphelenchus xylophilus. For Pathol, 2019 49 2 e12467

[30] Wilkins MR, Gasteiger E, Bairoch A, Sanchez JC, Williams KL, Appel RD, Hochstrasser DF. Protein identification and analysis tools in the ExPASy server. Methods Mol Biol, 1999, 112: 531-552.

[31] Yang XG, Jia LY, Wei L, Zuo WS, Song SP. Correlation of expression levels of multidrug resistance gene 1 (mdr1) mRNA, multidrug resistance-associated protein (MRP), amd P-glycoprotein (P-gp) with chemotherapy efficacy in malignant lymphomas. Zhonghua Yi Xue Za Zhi, 2002, 82: 1177-1179.

[32] Zhang HK, Gao SH, Lercher MJ, Hu SN, Chen WH. EvolView, an online tool for visualizing, annotating and managing phylogenetic trees. Nucleic Acids Res, 2012, 40: W569-W572.

[33] Zhao BG. Zhao BG, Futai K, Sutherland JR, Takeuchi Y. Pine wilt disease in China. Pine wilt disease, 2008, Tokyo: Springer