- Research article
- Open Access
Characterization and expression profiling of ATP-binding cassette transporter genes in the diamondback moth, Plutella xylostella (L.)
BMC Genomics volume 17, Article number: 760 (2016)
ATP-binding cassette (ABC) transporters are one of the major transmembrane protein families found in all organisms and play important roles in transporting a variety of compounds across intra and extra cellular membranes. In some species, ABC transporters may be involved in the detoxification of substances such as insecticides. The diamondback moth, Plutella xylostella (L.), a destructive pest of cruciferous crops worldwide, is an important species to study as it is resistant to many types of insecticides as well as biological control Bacillus thuringiensis toxins.
A total of 82 ABC genes were identified from our published P. xylostella genome, and grouped into eight subfamilies (ABCA-H) based on phylogenetic analysis. Genes of subfamilies ABCA, ABCC and ABCH were found to be expanded in P. xylostella compared with those in Bombyx mori, Manduca sexta, Heliconius melpomene, Danaus plexippus, Drosophila melanogaster, Tetranychus urticae and Homo sapiens. Phylogenetic analysis indicated that many of the ABC transporters in P. xylostella are orthologous to the well-studied ABC transporter genes in the seven other species. Transcriptome- and qRT-PCR-based analysis elucidated physiological effects of ABC gene expressions of P. xylostella which were developmental stage- and tissue-specific as well as being affected by whether or not the insects were from an insecticide-resistant strain. Two ABCC and one ABCA genes were preferentially expressed in midgut of the 4th-instar larvae of a susceptible strain (Fuzhou-S) suggesting their potential roles in metabolizing plant defensive chemicals. Most of the highly expressed genes in insecticide-resistant strains were also predominantly expressed in the tissues of Malpighian tubules and midgut.
This is the most comprehensive study on identification, characterization and expression profiling of ABC transporter genes in P. xylostella to date. The diversified features and expression patterns of this gene family may be associated with the evolutionary capacity of this species to develop resistance to a wide range of insecticides and biological toxins. Our findings provide a solid foundation for future functional studies on specific ABC transporter genes in P. xylostella, and for further understanding of their physiological roles and regulatory pathways in insecticide resistance.
ATP-binding cassette (ABC) transporters constitute one of the largest transmembrane protein families, which is widespread across all species. The first ABC transporter was found in prokaryotes, and the first cloned and characterized human ABC transporter member was ABCB1, which confers multidrug resistance (MDR) to cancer cells, preventing the accumulation of chemotherapeutic drugs . According to their functions, the ABC proteins can be divided into three categories: importers, exporters and non-transport proteins . Importers are only found in prokaryotes. In eukaryotes, ABC transporters are exporters and involved in the excretion of drugs, and endogenous and exogenous toxins. The third class of ABC proteins are apparently not related to molecule transport, but rather acting as ion channels, regulators of ion channels and receptors, and in some cases involved in DNA repair, ribosome assembly and translation .
ABC transporters are generally composed of four core domains with two being nucleotide-binding domains (NBDs) that bind and hydrolyze ATP and two transmembrane domains (TMDs) mediating translocation of the respective substrate . Reflecting the diversity of substrates handled by ABC transporters, TMDs are much more diverse than NBDs, whose sequences are highly conserved in order to perform their roles as ATP hydrolyzing enzymes . Half-transporters consist of one NBD and one TMD and need to form homo- or heterodimers to be functional . The mode of action of ABC transporters is called “ATP-switch” transport cycle [7, 8]. The ABC family is classified into eight subfamilies, annotated A to H according to their sequence similarity and domain conservation . The H subfamily appears to be present in all insects, mites, the slime mould Dictyostelium, and zebrafish, but is absent in genomes of plants, worms, yeasts, and mammalian species [9, 10]. In humans, many ABC proteins have been characterized with special functions and mutations of ABC genes can cause or contribute to a series of genetic disorders .
ABC transporters have recently been documented in insects as a family of detoxification-involved proteins , complementing the activity of another three classes of major metabolic enzymes: cytochrome P450 monooxygenases (P450s), glutathione S-transferases (GSTs) and carboxylesterases (COEs). Detoxification of toxins occurs in three phases with the first phase being associated predominantly with P450s . In phase II, GSTs are dominant and known to be linked to resistance development to most classes of insecticides . A series of transporters, including members of ABC transporters, are involved in phase III, which aims at the elimination of products generated during phases I and II . Some ABC members of subfamilies B, C and G are involved in resistance to xenobiotics including insecticides . Epis et al.  report that combining the insecticide permethrin with the ABC transporter inhibitor leads to greater Anopheles stephensi mortality than when using permethrin alone, demonstrating the importance of ABC transporters in insecticide resistance. Besides their detoxification roles, RNAi-mediated knockdown of some ABC genes in Tribolium castaneum results in a series of abnormal developmental phenotypes, such as growth arrest, eye pigmentation defects, abnormal cuticle formation, egg-laying and egg-hatching defects, and mortality .
The insect pest Plutella xylostella is a cosmopolitan Lepidoptera that almost exclusively feeds on cruciferous plants . Due to its short life cycle and capacity to rapidly develop insecticide resistance, P. xylostella is difficult to control [19, 20]. The species is the first to be reported resistant to dichlorodiphenyltrichloroethane (DDT) in the 1950s  and Bacillus thuringiensis (Bt) toxins in the 1990s . Bt resistance of P. xylostella is associated with ABCC2 alone  or in combination with ABCC3  or ABCG1 . In addition, the silencing of an ABCH1 gene results in the death of larvae and pupae . Expression of ABC genes is found to be more frequently up-regulated than that of GSTs, COEs or P450s in insecticide-resistant larvae of P. xylostella, suggesting a potential detoxification role of ABC transporters .
In this study, based on the ABC transporter genes (PxABCs) previously identified from the P. xylostella genome , we further characterized the gene structure and motifs, and performed phylogenetic analysis using P. xylostella, Bombyx mori, Manduca sexta, Heliconius melpomene, Danaus plexippus, Drosophila melanogaster, Tetranychus urticae, and Homo sapiens to further understand the evolutionary relationships among the eight subfamilies identified in this study. In addition, we carried out transcriptome- and qRT-PCR-based expression profiling of the ABC transporter genes in different developmental stages, tissues, and insecticide-susceptible and resistant strains of P. xylostella.
Results and discussion
Identification and grouping of the PxABCs
Based on our previous work on annotation of PxABCs in the P. xylostella genome , we identified 82 ABC transporter genes (Table 1 and Additional file 1) and 19 ABC fragments (Additional file 2). The 19 ABC fragments had homology to ABC transporters of other insects, but lacked the highly conserved NBDs of canonical ABC proteins . P. xylostella ABC transporter genes were grouped into the eight subfamilies (A-H) (Additional file 3). The number of genes in each subfamily greatly varied, ranging from one gene in ABCE to 21 in ABCC (Table 2). The ABCC subfamily was further divided into two groups with one group highly similar to the ABCB subfamily, which was also found in the other Lepidoptera, B. mori .
Characterization of the PxABCs and their motifs
The 82 PxABCs were dispersed on 59 scaffolds, 40 of which were found being individually located on different scaffolds. The remaining PxABCs were clustered on 19 scaffolds with each containing two or three genes, suggesting tandem duplication of these genes. The length of most predicted ABC transporters ranged from 124 to 2,714 amino acids (aa), with two exceptionally long genes containing 3,796 and 4,008 aa. The corresponding exon numbers ranged from 2 to 73 (Table 1), showing high structural complexity (Additional file 3).
The NBDs of ABC transporters generally contain seven highly conserved, but not invariant, motifs including Walker A, Walker B, ABC signature, A-loop, Q-loop, D-loop and H-loop . The Walker motifs A and B show that the ATP binding sites  and the “signature sequence” (also called C-loop or LSGGQ motif) are exclusive to ABC transporters, contributing to the formation of a composite catalytic site . Additional motifs such as the A-loop, Q-loop, D-loop and H-loop (also called the switch loop) only possess a single conserved aa residue , which are key amino acids in the catalytic cycle. We identified in P. xylostella six of these seven conserved motifs, with A-loop being absent (Fig. 1). However, the nature of such structural invariance in NBDs among arthropods remains poorly understood .
Subfamily-based comparison of ABC transporters
In P. xylostella, 15 ABCA genes were identified, including one 3-gene cluster on scaffold 295, two 2-gene clusters on scaffolds 188 and 856, and eight genes that were individually located on different scaffolds. Apparently, gene tandem duplication resulted in high protein diversity of ABCA. ABCA subfamily harbored the two longest ABC transporter genes, Px013614 (3,796 aa) and Px017838 (4,008 aa), which were similar in length to B. mori  and T. urticae . In Saccharomyces cerevisiae (yeast), no ABCA member has been identified  (Table 2). ABCA transporters in human are characterized as full-transporters , however, we found their structures are variable in arthropods. The ABCA transporters of T. urticae and D. melanogaster were all full-transporters, while the moths and butterflies contained both full- and half-transporters (Additional file 4). The P. xylostella ABCA subfamily comprised of nine full- and six half-transporters. Phylogenetic analysis of the ABCA subfamily revealed a Lepidoptera-specific clade with a distinct expansion in P. xylostella, which was orthologous to a human-specific clade (Fig. 2). A same orthologous relationship was also observed between another Lepidoptera-specific clade and a mite-specific clade (Fig. 2).
ABCA transporters in human are involved in lipid transport and metabolism . For example, the ABCA1 protein is mutated in a recessive disorder (Tangier disease) characterized by a defect in biogenesis of high density lipoprotein (HDL) . It is reported that HsABCA4 is linked to the Stargardt disease because of its role in retinal integrity, HsABCA3 is involved in lung surfactant production, and HsABCA12 participates in keratinization processes in the skin . However, the roles of ABCA transporter genes in arthropods are currently unclear. They may share a function related to lipid trafficking processes based on the high conservation of their structure. RNAi-mediated knockdown of TcABCA-9A or TcABCA-9B results in ~30 % mortality of T. castaneum , implying their significant roles in insects.
ABCB subfamily contains both full- and half-transporters . In the P. xylostella genome, we identified seven full- and seven half-transporters. The phylogenetic analysis of full-transporters showed that H. sapiens and T. urticae were located in a separate clade from insects (Fig. 3a). The phylogenetic tree of most ABCB half-transporters showed obvious orthologs and high bootstrap values among the eight species (Fig. 3b), indicating that they were orthologous and evolutionary divergent. Our results suggest that full-transporters may have evolved through lineage-specific duplication, while half-transporters may be evolutionarily conserved in metazoan species [10, 28, 33].
In human, the full-transporter ABCB1 plays important roles in cancer development by contributing to MDR . Expression of the ABCB1 gene of P. xylostella, PxPgp1, was found being up-regulated in the strain subjected to insecticide abamectin . In Heliothis virescens, increased expression of the human Pgp orthologs is associated with resistance to thiodicarb . Mutations in bile salt export pump gene (BSEP, HsABCB11) and multidrug-resistant gene 3 (MDR3, HsABCB4) in human can result in progressive familial intrahepatic cholestasis . Based on the phylogenetic tree (Fig. 3a), HsABCB4 and HsABCB11, both being full-transporters, seemed to be functioning in a human-specific manner. Human ABCB6, 7, 8 and 10 are mitochondrial half-transporters involved in iron metabolism and transport of Fe/S protein precursors . Their orthologs in arthropods are clear (Fig. 3b), implying the similar roles as in humans . However, the human-specific ABCB half-transporters (Fig. 3b), ABCB2, 3 and 9, are involved in antigen processing .
In P. xylostella, ABCC transporters formed the largest subfamily with 21 members (Table 1). There were one 3-gene cluster on scaffold 137, two 2-gene clusters on scaffolds 21 and 321, and 14 genes that were individually located on different scaffolds, suggesting tandem duplication of those clustered genes. Human ABCC transporters are all full-transporters, but in insects both full- and half-transporters can be found [3, 15, 28]. It showed similar structural characteristics to what we have found in ABCA subfamily of the eight species. The P. xylostella ABCC subfamily consisted of eight full- and 13 half-transporters. The phylogenetic tree of ABCC transporters showed that lineage-specific members were significantly duplicated in T. urticae, while this kind of duplication was inconspicuous in other species, suggesting varying evolutionary divergence of ABCC subfamily among the studied species (Fig. 4).
ABCC transporters are regarded as multidrug-resistance associated proteins (MRPs) due to their ability to extrude drugs with broad specificity. In D. melanogaster adults, MRPs may play a role in secretion of methotrexate by the Malpighian tubules , which is a drug used for treating cancers and auto-immune disorders . Labbe et al. propose that ABCC transporters can be involved in xenobiotic efflux of moths . Similar to ABCB4, ABCC2 is also involved in bile acids, phospholipids, and bilirubin export and related to progressive familial intrahepatic cholestasis . ABCC7 (cystic fibrosis transmembrane conductance regulator, CFTR) is an important ABC transporter involved in ion transport in human, and those who have mutations on functional CFTR may suffer diseases [40, 41]. However, no orthologs of human ABCC7 were found in P. xylostella and the other arthropod species.
Although mutations in ABCC2 have been implicated in Bt resistance in H. virescens, P. xylostella, Trichoplusia ni, and B. mori [23, 42–44], species-specific patterns of cross-resistance to Bt toxins are varied. This suggests that different mutations in ABCC2 may have occurred influencing the types of Bt resistance in a toxin binding site-dependent manner . Recently, ABCC3 ortholog in P. xylostella (Px008999) has been reported to also be associated with Bt resistance in insects , while its ortholog in human (ABCC3) functions as a marker for MDR in non-small cell lung cancer .
ABCC transporters also have cell-surface receptor activity, such as sulfonylurea receptor (SUR) . Among the eight species studied, P. xylostella and M. sexta had no ortholog of SUR (Fig. 4). In arthropods, SUR has been proposed as the direct target site for benzoylphenylureas (BPUs), a group of chitin synthesis inhibitors . The insecticidal activity of BPUs in P. xylostella  provides a clue for further exploring its SUR gene. Recent research however suggests that SUR may not be the only receptor for chitin synthesis as shown in D. melanogaster embryos . Knocking down SUR ortholog in T. castaneum has no effect on chitin synthesis .
The ABCD subfamily solely consists of half-transporters with the topology TMD–NBD, except for some plant representatives . All three ABCD members in P. xylostella were half-transporters. ABCD subfamily was present in all studied taxa and comprised of two to five members (Table 2). The ABCDs showed four distinct groups with HsABCD4 standing in its own branch and the three PxABCDs being clustered into two Lepidoptera-specific clades (Fig. 5). Their high orthologous relationships indicate that they are evolutionarily conserved in metazoan species.
ABCD transporters are peroxisomal membrane-located and are involved in importing fatty acids and/or fatty acyl-CoAs into peroxisome with specific metabolic and developmental functions in various organisms . However, the role of ABCD transporters in arthropods is still unclear .
ABCE and ABCF
ABCE proteins showed a high similarity (>50 % of aa identity) among the studied species (Fig. 6). The numbers of ABCE and ABCF were highly conserved, with most eukaryotes (including P. xylostella) having one ABCE and three ABCF genes (Table 2). Phylogenetic analysis showed that they had clear orthologous relationships (Fig. 6).
The structure of ABCE and ABCF are quite distinct from other ABC transporters due to their lack of TMDs and only containing a pair of linked NBDs . Therefore, they are not related to transport. Initially, ABCE was identified as an RNase L inhibitor in human . In human and yeast, ABCE proteins contribute in translation initiation . ABCF is similar to ABCE and relates to ribosome biogenesis and translational control. RNAi-mediated knockdown of members in the ABCE and ABCF subfamilies in the flour beetle T. castaneum results in 100 % mortality in penultimate larvae, showing fundamental roles of ABCE and ABCF in biological processes of the insects . In P. xylostella, the ABCE gene (Px007660) exhibited the highest expression level compared to other PxABCs based on the P. xylostella genome  and transcriptome . The same result is also reported in B. mori .
In P. xylostella, we found a total of 19 ABCG transporters, with 18 of them being half-transporters and only one full-transporter (Px007949). Gene duplication appeared to have occurred several times among ABCG genes of P. xylostella, which was evidenced by the location of gene paralogs on scaffolds 3, 8, 82, 131 and 282 (Table 1). Most ABCG transporters are half-transporters and they need to form homo- or hetero-dimers to perform the transport function .
In the phylogenetic tree of ABCG, all the orthologs of Dmwhite, Dmbrown and Dmscarlet were identified in Lepidoptera (Fig. 7). They are the most well-known ABCG genes in D. melanogaster, and encode for proteins to transport guanine or tryptophan precursors of the red and brown eye color pigments influencing the development of compound and simple eyes [56, 57]. These orthologs are also present in other insects, such as Anopheles gambiae , Bactrocera dorsalis  and Ceratitis capitata , suggesting that white, brown and scarlet are highly conserved protein transporters for eye color pigments in insects. Although several P. xylostella ABCG proteins have similar sequences with the white, brown and scarlet proteins in D. melanogaster, it is unclear whether they are involved in eye color pigments. It has been reported that Pxwhite is associated with Bt resistance in P. xylostella . In A. stephensi, analysis shows that the ortholog of Dmscarlet, AnstABCG4, exhibits a ten-fold increased expression when treated with permethrin after 48 h compared to untreated control .
In human, ABCG proteins are involved in lipid transport across membranes. The most intensively studied ABCG in human is ABCG2 (breast cancer resistance protein, BCRP), which acts as a MRP transporting anticancer drugs and a series of substrates [61–63]. In our study, no human ABCG2 orthologs were identified in Lepidoptera (Fig. 7), which were identified in D. melanogaster (E23) and T. urticae. D. melanogaster E23 (DmCG3327) is a 20-hydroxyecdysone (20E) primary response ABC transporter, which can suppress 20E-mediated gene activation . In B. mori, five ABCG genes including BmABC005226, BmABC005203, BmABC005202, BmABC010555 and BmABC010557, are 20E responsive genes although they are not orthologous to E23 . The orthologous relationships of HsABCG5 and HsABCG8 with the other species showed their evolutionarily conserved function (Fig. 7), the translated proteins of which form a functional heterodimer that is coordinately regulated by cholesterol .
We identified 6 ABCH genes in P. xylostella and they were all half-transporters. One of them (Px004510) tended to share high sequence similarity with ABCF subfamily (Additional file 3). ABCH subfamily was first identified with the sequencing of D. melanogaster genome . The structure of ABCH proteins is most closely related to subfamily ABCG. A previous study using five insect species with available genomes suggests that insect ABCHs may have originated from a common ancestral copy . Our phylogenetic analysis further supports this hypothesis (Fig. 8). Outside the Insecta, however, the T. urticae ABCHs seem not to conform to this rule (Fig. 8).
Down-regulation of an ABCH1 gene can cause high mortality of larvae and pupae of P. xylostella . Similarly, RNAi-mediated knockdown of TcABCH-9C results in desiccation and 100 % mortality in T. castaneum , and knockdown of ABCH gene CG9990 in D. melanogaster is lethal [65, 66]. It may be inferred that ABCH is a newly evolved subfamily with specific functions in the development of certain species, but this needs to be investigated.
Stage- and tissue-specific expression of the PxABCs
In order to develop an understanding of the physiological functions of ABC transporters in P. xylostella, the expression of 82 PxABCs was profiled at different developmental stages and with different tissues of the susceptible strain (Fuzhou-S, SS) (Fig. 9 and Additional file 13), based on the P. xylostella genome and transcriptome datasets [27, 55]. The RPKM values of the 82 ABC genes were clustered into three separate clades. The first clade contained 22 ABC genes with higher expression levels in various tissues and developmental stages, including two ABCDs, one ABCE and two ABCFs, which showed particularly high expression levels. This indicated that while ABCD, ABCE and ABCF were numerically small subfamilies (Table 2), they might be of some biological importance, though currently unclear. Within the second clade, there were 25 ABC genes exhibiting low expressions with most RPKM values being > 1.
The third clade contained 35 ABC genes with variable expression patterns. We found that Px002415 (ABCC3), Px002416 (ABCC2) and Px008256 (ABCA7) were relatively highly expressed in midgut of 4th instar and all larval stages compared with the other tissues and stages. Previous studies have shown that ABCC2 proteins play important physiological roles in resistance to Bt toxins [23, 42, 43]. Since midgut is an important organ involved with food digestion and detoxification, and larva is the main feeding stage of P. xylostella, we suspect that these genes may also play important roles in metabolizing plant secondary compounds. Many ABC genes were preferentially expressed in Malpighian tubule and most of them belonged to ABCB and ABCC subfamilies (Fig. 9). This suggests that ABC transporter genes of these subfamilies may also be involved in xenobiotic detoxification because Malpighian tubules make up an important excretory organ for transporting wastes .
Strain-specific expression profiles of the PxABCs
In order to understand the potential function of the P. xylostella ABC proteins in insecticide resistance, expressions of the P. xylostella ABC genes were profiled for larvae of two insecticide-resistant strains (fipronil (FRS) and chlorpyrifos (CRS)) (Fig. 10 and Additional file 14), based on the genome and transcriptome data and compared with the susceptible strain (SS) [27, 55]. The majority of the P. xylostella ABC genes were up-regulated in FRS and/or CRS, with a few genes being down-regulated. The ABCE gene (Px007660) showed the highest expression level in all three strains, suggesting its fundamental role in various physiological processes of the cell. At least, five ABCC genes including Px002418 (ABCC4), Px008999 (ABCC13), Px009835 (ABCC17), Px002416 (ABCC2), and Px002419 (ABCC5) were up-regulated in FRS and CRS. Among them, ABCC genes are proposed to be involved in insecticide resistance , and ABCC2 is linked to Bt resistance .
In our case, one of the ABCH genes (Px014955) was mostly up-regulated in FRS and CRS, and weakly expressed in SS. Four genes including Px002784 (ABCC7), Px013614 (ABCA11), Px013659 (ABCA12) and Px015888 (ABCC21) showed up-regulated expression in CRS, but not in FRS, suggesting that the expression of these genes could be induced by specific insecticides. There were eight PxABCs exhibiting low expression (RPKM value < 1) in all three strains. Three of the PxABC genes (ABCA14, ABCB4 and ABCG4) had no differential expression in the three strains, suggesting that they might not be related to the resistance of chlorpyrifos and fipronil. Px006766 (ABCF3) and Px013659 (ABCA12) were down-regulated in FRS and CRS, while Px013614 (ABCA11), Px008371 (ABCG11) and Px004725 (ABCG5) exhibited up-regulated expression in both insecticide-resistant strains. This may indicate that not all ABC transporter genes are responsible for xenobiotics excretion, and a strictly transcriptional regulation of ABCs in response to external challenges is required for balancing transmembrane substrate transport.
qRT-PCR validation of expressions
Eighteen ABC transporter genes with up- or down-regulated expressions in the two insecticide-resistant strains (FRS and CRS) of P. xylostella were selected for experimental validation using qRT-PCR. Quantitative expression of the PxABCs at different developmental stages of the susceptible strain (Additional file 15) and two insecticide-resistant strains (Additional file 16) overall confirmed the stage- and strain-specific expression patterns as profiled using the P. xylostella genome and transcriptome datasets (Figs. 9 and 10). The results provide significant clues for further studies on functions of some specific PxABCs in host plant adaptation and insecticide resistance development.
Our work presents the most comprehensive study to date on identification, characterization and expression profiling of ABC transporters in the genome of P. xylostella, integrating evolutionary and physiological aspects. A comparison of ABC transporters from seven arthropod species and human provides an overview of this vital gene family. The variances in genomic features and expression patterns of the genes reflect evolutionary and functional diversification of the ABC transporters. Some genes are preferentially expressed in larvae, midgut, and Malpighian tubules, which may be involved in detoxification of plant defense chemicals. Moreover, some of the ABC transporter genes are up-regulated in two insecticide-resistant strains compared to the susceptible strain, suggesting their involvement in detoxification of the chemical insecticides through transportation. Our work offers a solid foundation for further research on gene-specific functions of the ABC transporters, which is essential to better understand the molecular and genetic mechanisms involved in detoxification in P. xylostella.
Experimental DBM strains
The experimental population of P. xylostella was derived from a susceptible strain (Fuzhou-S) that was collected from a vegetable field of Fuzhou (26.08°N, 119.28°E) in 2004 and used for genome sequencing . Since then this initial population has been reared on potted radish seedlings (Raphanus sativus L.) at 25 ± 1 °C, 65 ± 5%RH and L:D = 16:8 h in a separate greenhouse without exposure to insecticides. The two insecticide-resistant strains (CRS and FRS) were selected from this susceptible strain, and detailed in the published P. xylostella transcriptome .
Identification of lepidopteran ABC transporters
To identify the ABC transporter genes in the P. xylostella genome (http://iae.fafu.edu.cn/DBM/), a systematic BLASTP search was performed using arthropod ABC transporter protein sequences available from NCBI as queries, with the cutoff set at e-value < e−20. Candidate ABC transporter sequences were submitted to the NCBI protein database to search for ABC transporter domains. We used the online FGENESH and FGENESH+ programs (http://linux1.softberry.com/berry.phtml) to predict the gene structures of ABC transporter genes on their genomic DNA sequences. Finally, we used program Pfam (http://pfam.xfam.org/) to identify the NBD and TMD structure. The same method was applied to identify ABC transporters from up to date genome versions of Manduca sexta (ftp://ftp.bioinformatics.ksu.edu/pub/Manduca/OGS2/), Danaus plexippus (http://monarchbase.umassmed.edu/resource.html), and Heliconius melpomene (http://www.butterflygenome.org/node/4).
Sequence alignment and phylogenetic tree construction
To assign the ABC transporter genes to specific subfamilies, the complete aa of NBDs of all P. xylostella ABC transporters were aligned by MUSCLE algorithm. Phylogenetic trees were constructed with the Maximum Likelihood method using MEGA-CC (7.0.18) for Linux users. Bootstrap analysis with 1,000 replicates was used to evaluate the significance of the nodes. Poisson correction aa model and all sites were used for the tree reconstruction. Comparison analyses were conducted among P. xylostella, B. mori, M. sexta, H. melpomene, D. plexippus, D. melanogaster, T. urticae and H. sapiens for each of the ABC transporter subfamilies separately using the full-length aa sequences and the same method as what phylogenetic analysis of P. xylostella ABC transporters used.
Examination of gene structure and motifs
The P. xylostella ABC genes were mapped onto corresponding scaffolds of the genome sequence assembly version 2 . Gene structure of the P. xylostella ABC genes was visualized using the online tool Gene Structure Display Server (http://gsds.cbi.pku.edu.cn/). We used the MEME software (http://meme.sdsc.edu) to identify the conserved motifs in the 82 PxABC sequences using the following parameters: number of repetitions = any, maximum number of motifs = 10, and optimum motif width = 3 to 10 residues.
Gene expression profiling of the PxABCs
The RNA-seq data of the P. xylostella ABC genes were downloaded from the published database (http://iae.fafu.edu.cn/DBM/). Expressions of the PxABCs were profiled in different developmental stages and tissues of the susceptible strain, including eggs (within 24 h after oviposition), 1st ~ 4th-instar larvae, pupae, male and female adults, fat body, hemolymph, Malpighian tubules, silk glands, salivary glands and midgut of the 4th-instar larvae, head of the 4th-instar larvae and male/female adults, and individuals between 3rd-instar larvae of SS, and each of the insecticide-resistant strains (FRS and CRS). Each of the RPKM values was transformed into base-2 logarithm, and the expression profiling of ABC transporter genes was generated and visualized by pheatmap package of the R program (https://www.r-project.org/, version 3.2.5) using the similarity metric of Euclidean distance and clustering method of complete linkage.
Validation of the gene expression by qRT-PCR
To confirm expression patterns of the PxABCs based on P. xylostella transcriptome, qRT-PCR analysis was performed using SYBR-green fluorescence with gene-specific primers. The PCR products were examined by dissociation curve analysis after the PCR reaction to confirm the specific detection of target transcripts by the qRT-PCR analysis. Three independent biological replicates were included for qRT-PCR, each of which had three technical replicates. The first-strand cDNA was synthesized from total RNA using the reverse transcriptase kit from Promega (Madison, WI). The qRT-PCR reactions were prepared using the SYBR SELECT MASTER MIX FOR CFX from Invitrogen (Carlsbad, CA) following manufacturer’s instructions and run on a CFX96 Touch™ Real-Time PCR Detection System (Bio-Rad, USA), following the program: 95 °C for 3 min; 45 cycles of 95 °C for 15 s, and 57 °C for 35 s, and a final melt curve at 60 °C for 5 s to 95 °C with 0.5 °C increments. The P. xylostella ribosomal protein L32 (RPL32) gene (GenBank acc. no. AB180441) was used as an internal reference. Standard curves were generated by 5-fold dilutions of the cDNA templates. The 2−ΔCt method was used to analyze the relative values of mRNA expression.
ABC transporter of A. stephensi
Breast cancer resistance protein
ABC transporter of B. mori
Cystic fibrosis transmembrane conductance regulator
ABC transporter of H. sapiens
Multidrug-resistance associated protein
ABC transporter of P. xylostella
P-glycoprotein 1 of P. xylostella
Quantitative real-time polymerase chain reaction
Reads per kilobase per millionread
ABC transporter of T. castaneum
Gottesman MM, Fojo T, Bates SE. Multidrug resistance in cancer: role of ATP–dependent transporters. Nat Rev Cancer. 2002;2(1):48–58.
Saurin W, Hofnung M, Dassa E. Getting in or out: early segregation between importers and exporters in the evolution of ATP-binding cassette (ABC) transporters. J Mol Evol. 1999;48(1):22–41.
Dean M, Hamon Y, Chimini G. The human ATP-binding cassette (ABC) transporter superfamily. J Lipid Res. 2001;42(7):1007–17.
Rees DC, Johnson E, Lewinson O. ABC transporters: the power to change. Nat Rev Mol Cell Biol. 2009;10(3):218–27.
Seeger MA, van Veen HW. Molecular basis of multidrug transport by ABC transporters. Biochim Biophys Acta-Proteins Proteom. 2009;1794(5):725–37.
Hollenstein K, Dawson RJ, Locher KP. Structure and mechanism of ABC transporter proteins. Curr Opin Struct Biol. 2007;17(4):412–8.
Higgins CF, Linton KJ. The ATP switch model for ABC transporters. Nat Struct Mol Biol. 2004;11(10):918–26.
Linton KJ. Structure and function of ABC transporters. Physiology. 2007;22(2):122–30.
Dean M, Annilo T. Evolution of the ATP-binding cassette (ABC) transporter superfamily in vertebrates. Annu Rev Genomics Hum Genet. 2005;6:123–42.
Dermauw W, Osborne EJ, Clark RM, Grbić M, Tirry L, Van Leeuwen T. A burst of ABC genes in the genome of the polyphagous spider mite Tetranychus urticae. BMC Genomics. 2013;14(1):317.
Wenzel JJ, Piehler A, Kaminski WE. ABC A-subclass proteins: gatekeepers of cellular phospho-and sphingolipid transport. Front Biosci. 2006;12:3177–93.
Labbe R, Caveney S, Donly C. Genetic analysis of the xenobiotic resistance-associated ABC gene subfamilies of the Lepidoptera. Insect Mol Biol. 2011;20(2):243–56.
Feyereisen R. Insect P450 enzymes. Annu Rev Entomol. 1999;44:507–33.
Ranson H, Hemingway J. Mosquito glutathione transferases. Methods Enzymol. 2005;401:226–41.
Chahine S, O’Donnell MJ. Interactions between detoxification mechanisms and excretion in Malpighian tubules of Drosophila melanogaster. J Exp Biol. 2011;214(3):462–8.
Epis S, Porretta D, Mastrantonio V, Comandatore F, Sassera D, Rossi P, et al. ABC transporters are involved in defense against permethrin insecticide in the malaria vector Anopheles stephensi. Parasite Vector. 2014;7(1):349.
Broehan G, Kroeger T, Lorenzen M, Merzendorfer H. Functional analysis of the ATP-binding cassette (ABC) transporter gene family of Tribolium castaneum. BMC Genomics. 2013;14(1):6.
Talekar NS, Shelton AM. Biology, ecology, and management of the diamondback moth. Annu Rev Entomol. 1993;38(1):275–301.
Sarfraz M, Keddie BA. Conserving the efficacy of insecticides against Plutella xylostella (L.) (Lep., Plutellidae). J Appl Entomol. 2005;129(3):149–57.
Furlong MJ, Wright DJ, Dosdall LM. Diamondback moth ecology and management: problems, progress, and prospects. Annu Rev Entomol. 2013;58:517–41.
Ankersmit GW. DDT-resistance in Plutella maculipennis (Curt.) (Lep.) in Java. B Entomol Res. 1953;44(3):421–5.
Heckel DG, Gahan LJ, Liu YB, Tabashnik BE. Genetic mapping of resistance to Bacillus thuringiensis toxins in diamondback moth using biphasic linkage analysis. Proc Natl Acad Sci U S A. 1999;96(15):8373–7.
Baxter SW, Badenes-Pérez FR, Morrison A, Vogel H, Crickmore N, Kain W, et al. Parallel evolution of Bacillus thuringiensis toxin resistance in Lepidoptera. Genetics. 2011;189(2):675–9.
Guo ZJ, Kang S, Chen DF, Wu QJ, Wang SL, Xie W, et al. MAPK signaling pathway alters expression of midgut ALP and ABCC genes and causes resistance to Bacillus thuringiensis Cry1Ac toxin in diamondback moth. PLoS Genet. 2015;11(4):e1005124.
Guo ZJ, Kang S, Zhu X, Xia JX, Wu QJ, Wang SL, et al. Down-regulation of a novel ABC transporter gene (Pxwhite) is associated with Cry1Ac resistance in the diamondback moth, Plutella xylostella (L.). Insect Biochem Mol Biol. 2015;59:30–40.
Guo ZJ, Kang S, Zhu X, Xia JX, Wu QJ, Wang SL, et al. The novel ABC transporter ABCH1 is a potential target for RNAi-based insect pest control and resistance management. Sci Rep. 2015;5:13728.
You MS, Yue Z, He WY, Yang XH, Yang G, Xie M, et al. A heterozygous moth genome provides insights into herbivory and detoxification. Nat Genet. 2013;45(2):220–5.
Liu SM, Zhou S, Tian L, Guo EN, Luan YX, Zhang JZ, et al. Genome-wide identification and characterization of ATP-binding cassette transporters in the silkworm. Bombyx mori BMC Genomics. 2011;12(1):491.
Walker JE, Saraste M, Runswick MJ, Gay NJ. Distantly related sequences in the alpha-and beta-subunits of ATP synthase, myosin, kinases and other ATP-requiring enzymes and a common nucleotide binding fold. EMBO J. 1982;1(8):945–51.
Dermauw W, Van Leeuwen T. The ABC gene family in arthropods: comparative genomics and role in insecticide transport and resistance. Insect Biochem Mol Biol. 2014;45:89–110.
Taglicht D, Michaelis S. Saccharomyces cerevisiae ABC proteins and their relevance to human health and disease. Method Enzymol. 1998;292:130–62.
Santamarina-Fojo S, Remaley AT, Neufeld EB, Brewer Jr HB. Regulation and intracellular trafficking of the ABCA1 transporter. J Lipid Res. 2001;42(9):1339–45.
Sturm A, Cunningham P, Dean M. The ABC transporter gene family of Daphnia pulex. BMC Genomics. 2009;10(1):170.
Tian L, Yang J, Hou W, Xu B, Xie W, Wang S, et al. Molecular cloning and characterization of a P-glycoprotein from the diamondback moth, Plutella xylostella (Lepidoptera: Plutellidae). Int J Mol Sci. 2013;14(11):22891–905.
Lanning CL, Fine RL, Corcoran JJ, Ayad HM, Rose RL, Abou-Donia MB. Tobacco budworm P-glycoprotein: biochemical characterization and its involvement in pesticide resistance. Biochim Biophys Acta-Gen Subjects. 1996;1291(2):155–62.
Chen HL, Liu YJ, Chen HL, Wu SH, Ni YH, Ho MC, et al. Expression of hepatocyte transporters and nuclear receptors in children with early and late-stage biliary atresia. Pediatr Res. 2008;63(6):667–73.
Abele R, Tampé R. Modulation of the antigen transport machinery TAP by friends and enemies. FEBS Lett. 2006;580(4):1156–63.
Chahine S, O’Donnell MJ. Physiological and molecular characterization of methotrexate transport by Malpighian tubules of adult Drosophila melanogaster. J Insect Physiol. 2009;55(10):927–35.
Affleck JG, Neumann K, Wong L, Walker VK. The effects of methotrexate on Drosophila development, female fecundity, and gene expression. Toxicol Sci. 2006;89(2):495–503.
Pignatti PF, Bombieri C, Marigo C, Benetazzo M, Luisetti M. Increased incidence of cystic fibrosis gene mutations in adults with disseminated bronchiectasis. Hum Mol Genet. 1995;4(4):635–9.
Cohn JA, Friedman KJ, Noone PG, Knowles MR, Silverman LM, Jowell PS. Relation between mutations of the cystic fibrosis gene and idiopathic pancreatitis. N Engl J Med. 1998;339(10):653–8.
Gahan LJ, Pauchet Y, Vogel H, Heckel DG. An ABC transporter mutation is correlated with insect resistance to Bacillus thuringiensis Cry1Ac toxin. PLoS Genet. 2010;6(12):e1001248.
Atsumi S, Miyamoto K, Yamamoto K, Narukawa J, Kawai S, Sezutsu H, et al. Single amino acid mutation in an ATP-binding cassette transporter gene causes resistance to Bt toxin Cry1Ab in the silkworm, Bombyx mori. Proc Natl Acad Sci U S A. 2012;109(25):E1591–8.
Heckel DG. Learning the ABCs of Bt: ABC transporters and insect resistance to Bacillus thuringiensis provide clues to a crucial step in toxin mode of action. Pestic Biochem Physiol. 2012;104(2):103–10.
Hernández-Martínez P, Hernández-Rodríguez CS, Krishnan V, Crickmore N, Escriche B, Ferré J. Lack of Cry1Fa binding to the midgut brush border membrane in a resistant colony of Plutella xylostella moths with a mutation in the ABCC2 locus. Appl Environ Microbiol. 2012;78(18):6759–61.
Zhao YB, Lu HL, Yan A, Yang YM, Meng QW, Sun LC, et al. ABCC3 as a marker for multidrug resistance in non-small cell lung cancer. Sci Rep. 2013;3:3120.
Dean M, Rzhetsky A, Allikmets R. The human ATP-binding cassette (ABC) transporter superfamily. Genome Res. 2001;11:1156–66.
Abo-Elghar GE, Fujiyoshi P, Matsumura F. Significance of the sulfonylurea receptor (SUR) as the target of diflubenzuron in chitin synthesis inhibition in Drosophila melanogaster and Blattella germanica. Insect Biochem Mol Biol. 2004;34(8):743–52.
Li Y, Qin Y, Yang N, Sun Y, Yang X, Sun R, et al. Studies on insecticidal activities and action mechanism of novel benzoylphenylurea candidate NK-17. PLoS One. 2013;8(6):e66251.
Meyer F, Flötenmeyer M, Moussian B. The sulfonylurea receptor SUR is dispensable for chitin synthesis in Drosophila melanogaster embryos. Pest Manag Sci. 2013;69(10):1136–40.
Theodoulou FL, Holdsworth M, Baker A. Peroxisomal ABC transporters. FEBS Lett. 2006;580(4):1139–55.
Kerr ID. Sequence analysis of twin ATP binding cassette proteins involved in translational control, antibiotic resistance, and ribonuclease L inhibition. Biochem Biophys Res Commun. 2004;315(1):166–73.
Zhou A, Hassel BA, Silverman RH. Expression cloning of 2-5A-dependent RNAase A uniquely regulated mediator of interferon action. Cell. 1993;72:753–65.
Chen ZQ, Dong J, Ishimura A, Daar I, Hinnebusch AG, Dean M. The essential vertebrate ABCE1 protein interacts with eukaryotic initiation factors. J Biol Chem. 2006;281(11):7452–7.
He W, You M, Vasseur L, Yang G, Xie M, Cui K, et al. Developmental and insecticide-resistant insights from the de novo assembled transcriptome of the diamondback moth, Plutella xylostella. Genomics. 2012;99(3):169–77.
Morgan TH. Sex limited inheritance in Drosophila. Science. 1910;32(812):120–2.
Mackenzie SM, Brooker MR, Gill TR, Cox GB, Howells AJ, Ewart GD. Mutations in the white gene of Drosophila melanogaster affecting ABC transporters that determine eye colouration. Biochim Biophys Acta-Biomembranes. 1999;1419(2):173–85.
Besansky NJ, Bedell JA, Benedict MQ, Mukabayire O, Hilfiker D, Collins FH. Cloning and characterization of the white gene from Anopheles gambiae. Insect Mol Biol. 1995;4(4):217–31.
Handler AM, McCombs SD. The piggyBac transposon mediates germ-line transformation in the Oriental fruit fly and closely related elements exist in its genome. Insect Mol Biol. 2000;9(6):605–12.
Gomulski LM, Pitts RJ, Costa S, Saccone G, Torti C, Polito LC, et al. Genomic organization and characterization of the white locus of the Mediterranean fruitfly, Ceratitis capitata. Genetics. 2001;157(3):1245–55.
Tarr PT, Tarling EJ, Bojanic DD, Edwards PA, Baldan A. Emerging new paradigms for ABCG transporters. Biochim Biophys Acta-Mol Cell Biol L. 2009;1791(7):584–93.
Kerr ID, Haider AJ, Gelissen IC. The ABCG family of membrane-associated transporters: you don’t have to be big to be mighty. Br J Pharmacol. 2011;164(7):1767–79.
Fukuda Y, Schuetz JD. ABC transporters and their role in nucleoside and nucleotide drug resistance. Biochem Pharmacol. 2012;83(8):1073–83.
Hock T, Cottrill T, Keegan J, Garza D. The E23 early gene of Drosophila encodes an ecdysone-inducible ATP-binding cassette transporter capable of repressing ecdysone-mediated gene activation. Proc Natl Acad Sci U S A. 2000;97(17):9519–24.
Mummery-Widmer JL, Yamazaki M, Stoeger T, Novatchkova M, Bhalerao S, Chen D, et al. Genome-wide analysis of Notch signalling in Drosophila by transgenic RNAi. Nature. 2009;458(7241):987–92.
Zhang S, Feany MB, Saraswati S, Littleton JT, Perrimon N. Inactivation of Drosophila Huntingtin affects long-term adult functioning and the pathogenesis of a Huntington’s disease model. Dis Mod Mech. 2009;2(5–6):247–66.
The work was supported by National Natural Science Foundation of China (Nos. 31320103922, 31230061 and 31301677), National Key Project of Fundamental Scientific Research (“973” Programs, No. 2011CB100404) in China, Natural Science Foundation of Fujian Proince, China (2014 J01086), and Outstanding Youth Fellowships for WH at FAFU (xjq201403). GMG is supported by the National Thousand Talents Program in China and the Advanced Talents of SAEFA, and LV by the Minjiang Scholar Program in Fujian Province (PRC) and the Advanced Talents of State Administration of Foreign Experts Affairs (PRC).
Availability of data and materials
All sequences used for this study are available in respective public databases and provided in additional supporting files. The phylogenetic data for Figs. 2, 3, 4, 5, 6, 7 and 8 were uploaded to TreeBase and the original NEXUS files can be downloaded from the DBM-DB (http://iae.fafu.edu.cn/DBM/family/PxABCs.php).
WQ, WC and MZ designed and/or performed experiments. WQ drafted the manuscript. WQ, XM and WH carried out the bioinformatics analysis. MY supervised the project. MY, WQ, XM, WH, GMG and LV interpreted the data and critically revised the manuscript. All authors have read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
P. xylostella is not protected under any legislation in China, as a protected or endangered species regulating or restricting its collection. Neither specific permits are required for collecting the specimens from the field nor is animal ethics approval required for work with this invertebrate.
Protein sequences of the 82 ABC transporters of P. xylostella. (FA 17 kb)
Protein sequences of the 19 ABC transporter fragments of P. xylostella. (FA 78 kb)
Phylogenetic tree and gene structure of ABCs in P. xylostella. The left side of the figure presents the phylogenetic relationship among ABC transporters in P. xylostella and the right, the corresponding gene structures. The green bars represent exons, and black lines represent introns. The digits (0, 1 and 2) on each of the black lines represent three different phases of the introns existed in eukaryotic genes. The length of genes is measured with a scale on the bottom right. Protein sequences used for phylogenetic analysis are provided in the Additional file 1. (FA 30 kb)
Protein sequences of the ABCA transporters of the eight species. (XLSX 25 kb)
Protein sequences of the full ABCB transporters of the eight species. (XLSX 13 kb)
Protein sequences of the half ABCB transporters of the eight species. (DOCX 163 kb)
Protein sequences of the ABCC transporters of the eight species. (DOCX 273 kb)
Protein sequences of the ABCD transporters of the eight species. (FA 80 kb)
Protein sequences of the ABCE transporters of the eight species. (FA 7 kb)
Protein sequences of the ABCF transporters of the eight species. (DOCX 1284 kb)
Protein sequences of the ABCG transporters of the eight species. (FA 126 kb)
Protein sequences of the ABCH transporters of the seven species except H. sapiens. (FA 54 kb)
Log2 RPKM values of 82 PxABCs at different developmental stages and in different tissues of the susceptible P. xylostella. (FA 27 kb)
Log2 RPKM values of 82 PxABCs in larvae of the susceptible strain and two insecticide (chlorpyrifos and fipronil)-resistant strains of P. xylostella. (FA 172 kb)
Expression profiles of 18 selected ABC genes at different developmental stages of the susceptible P. xylostella as determined by qRT-PCR. X axis: samples; Y axis: relative expression values; E: eggs; 11: the 1st-instar larvae; 21–22: first and second day of the 2nd-instar larvae; 31–32: first and third day of the 3rd-instar larvae; 41–44: each day of the 4th-instar larvae; P1-p4: each day of pupa stage; AM: virgin male adult; AF: virgin female adult; A*M: post-mating male adult; A*F: post-mating female adult. (FA 12 kb)
Expression profiles of 18 selected ABC genes in the susceptible strain and two insecticide (chlorpyrifos and fipronil)-resistant strains of P. xylostella as determined by qRT-PCR. X axis: samples; Y axis: relative expression values; 3-rd: third day of the 3rd-instar larvae; 4-th: third day of the 4th-instar larvae; Pupae: third day of the pupa stage; Males: virgin male adults; Females: virgin female adults; SS: susceptible strain; CRS: chlorpyrifos-resistant strains; FRS: fipronil-resistant strain. (FA 5 kb)
About this article
Cite this article
Qi, W., Ma, X., He, W. et al. Characterization and expression profiling of ATP-binding cassette transporter genes in the diamondback moth, Plutella xylostella (L.). BMC Genomics 17, 760 (2016). https://doi.org/10.1186/s12864-016-3096-1
- ABC transporter
- Phylogenetic analysis
- Transcriptome analysis
- Insecticide resistance