dirofilaria immitis , a filaroid nematode that causes dirofilariasis or heartworm disease , is arguably the most important parasite that infects dogs in north america ( wolstenholme et  al . , 2015 ) . this mosquito - borne nematode affects 30 other mammal species including cats , wild canids , felids and occasionally humans ( lee et  al . , 2010 , the disease has worldwide distribution but is more prevalent in temperate , tropical and sub - tropical regions of the world ( simn et  al . , 2009 ) . for over 2 decades , prevention of the disease has relied upon a single class of drug , the macrocyclic lactones ( mls ) . concern about possible loss of efficacy ( loe ) to this class of drug was reported in 2005 ( hampshire , 2005 ) . recent studies have unequivocally confirmed ml resistance in d.  immitis ( bourguinat et  al . , 2011a , bourguinat et  al . , 2015 , pulaski et  al . , 2014 ) . the atp binding cassette ( abc ) transporters have been implicated in ml resistance ( see review lespine et  al . , specific studies on the role of abc transporters and ml resistance include studies in caenorhabditis elegans ( james and davey , 2009 , ardelli and prichard , 2013 , bygarski et  al . , 2014 ) , haemonchus contortus ( blackhall et  al . , 1998 ; xu et  al . , 1998 prichard and roulet , 2007 ) , onchocerca volvulus ( eng and prichard , 2005 , ardelli and prichard , 2007 , bourguinat et  al . 2014 ) , teladorsagia circumcincta ( dicker et  al . , 2011 ) , and for benzimidazole resistance in h.  contortus ( blackhall et  al . , 2008 ) , triclabendazole resistance in fasciola hepatica ( mottier et  al . , 2006 ) and praziquantel resistance in schistosoma mansoni ( messerli et  al . , 2009 ) . in d.  immitis , a diplotypic genotype in dim - pgp-11 ( abc - b class of the transporter superfamily ) gene was found to correlate strongly with a loss of response to ivermectin ( bourguinat et  al . members of the abc transporter superfamily , particularly p - glycoproteins ( pgps ) and multidrug resistance proteins ( mrps ) have broad substrate specificity ( e.g. , antibiotics , antimalarials , herbicides , antifungals and other chemotherapeutic agents ) ( higgins , 2007 ) . abc transporters , as efflux pumps , influence the bioavailability and disposition of drugs through active efflux of compound out of cells , thereby contributing to a phenomenon called multidrug resistance ( mdr ) ( lespine et  al . , 2012 , ardelli , 2013 ) . overexpression of these transporters has been suggested as a mdr conferring mechanism ( ardelli , 2013 ) . in nematodes , multidrug resistance is suggested to be conferred by members from subfamily b and c , though many reports specifically implicate the role of pgps in the resistance mechanism ( prichard et  al . , 2012 
 lespine et  al . , 2012 ) a study in b.  malayi has shown that transporters from subfamily a and g may also have a role in resistance , as they were also overexpressed along with pgps and mrps following treatment with ivermectin and moxidectin ( stitt et  al . although abc transporters may play a role in developing and maintaining anthelmintic resistance , it is believed that their main physiological function is to protect neurons and other tissues in which they are expressed , from a broad spectrum of toxins ( prichard and roulet , 2007 , ardelli , 2013 ) . unlike vertebrates , which have a vascular system contained within epithelial cells which line the blood capillaries and other blood vessels where pgps and mrps are commonly expressed , the number of efflux transporters is greater in nematodes , compared to mammals , as they may have evolved this diversity as part of their protective mechanism for a variety of different tissues and specific cells . abc transporters can regulate anthelmintic efficacy by modulating the intracellular or intra cell membrane concentration of lipophilic drugs . therefore , a deletion or inhibition of such transporters may partially overcome resistance ( lespine et  al . , 2012 , ardelli , 2013 , greenberg , 2013a , greenberg , 2014a ) . for example , deletion or disruption of mammalian and nematode pgp genes leads to increased sensitivity to ivermectin ( schinkel et  al . , 1994 , 2013 ) . also , substrates of pgp that are competitive inhibitors ( e.g. , verapamil ) , or that block the transport function directly , have restored anthelmintic sensitivity ( ardelli and prichard , 2013 ) . in b.  malayi , abc transporter inhibitors have potentiated sensitivity to ivermectin ( tompkins et  al . , 2011 ) . similar reversal mechanisms have been shown for h.  contortus resistant to ivermectin ( molento and prichard , 1999 , lifschitz et  al . inhibitors of mrp , such as mk571 and buthionine sulfoxamine , have also reversed resistance ( prichard et  al . , 2012 ) . in nematodes , abc transporters are believed to be essential for the survival of nematodes as they are known to be involved in a wide range of processes . for example , in c.  elegans , they play critical roles in apoptotic cell corpse removal ( wu and horvitz , 1998 ) , dauer formation ( yabe et  al . , 2005 ) , rnai ( sundaram et  al . , 2006 ) , directed sperm motility ( kubagawa et  al . , 2006 ) , apart from resistance to toxins ( broeks et  al . 2005 ) , and resistance to pathogens ( mahajan - miklos et  al . , 1999 ) . taking all of this into consideration , abc transporter genes may prove to be attractive therapeutic targets for new or repurposed mdr reversing agents , and/or for new anthelmintics ( lespine et  al . , 2012 , ardelli , 2013 , greenberg , 2014b ) . knowledge of allelic variants in abc transporter genes may be useful information to have during drug design . knowledge about polymorphism in these genes may help ensure that a new drug will be active against all the allelic forms of these targets . some polymorphism in mdr transporter genes may influence their level of expression and/or the functional characteristics ( gerloff , 2004 ) , and such knowledge may be useful to understand the mechanism of resistance . also , heterogeneity of a drug target may determine the drug selection process ( prichard , 2001 ) and knowledge of polymorphism in drug transporters may help to anticipate resistance development . in order to address these possibilities , the objective of this study was to identify putative abc transporter genes / subunits from subfamilies a , b , c and g in d.  immitis , by comparative genomic approaches , and to analyze the genes for single nucleotide polymorphisms ( snps ) . the genome of d.  immitis is not fully annotated ( bourguinat et  al . , 2015 ) and therefore a complete inventory of abc transporter genes is not yet available . in a recent study ( bourguinat et  al . , 2016 ) , eight abc transporter genes including a pseudogene from d.  immitis were described . to identify the remaining , un - annotated homologs of abc transporter genes in d.  immitis , we followed a similar approach as previously described ( mani et  al . , 2016 ) . briefly , all the nucleotide sequences that encode genes or subunits of abc transporter subfamilies ( a , b , c and g ) of all nematodes were extracted from available databases such as ncbi ( http://www.ncbi.nlm.nih.gov/ ) , wormbase ( http://www.wormbase.org/#01-23-6 ) , broad institute ( https://www.broadinstitute.org/ ) and nembase4 ( http://www.nematodes.org/nembase4/ ) . these sequences were then used in the nucleotide blastt server v2.2 ( http://nematodes.org/genomes/dirofilaria_immitis/ ) to indicate putative protein type . these blast searches were helpful to locate each of the putative abc transporter genes in the scaffolds of the d.  immitis nuclear genome ( version ndi.2.2.2 ) ( http://salmo.bio.ed.ac.uk/cgi-bin/gbrowse/gbrowse/ndi.2.2.2/ ) . whole genome data from diverse pooled samples of 122 adult d.  immitis , from 17 dogs from the usa , grand canary , grenada and italy , and approximately 32,000 d.  immitis microfilariae from 4 dogs from usa , described previously ( bourguinat et  al . , 2015 , mani et  al . , 2016 ) were used . more information on the d.  immitis samples and their classification into ml susceptible ( sus ) and loss of efficacy ( loe ) populations is given in the supplementary table st1 . the adult worms were mixed sex populations characterized as susceptible to macrocyclic lactone heartworm preventives , while the microfilariae were termed loss of efficacy samples as it was suspected that the dogs may have become infected despite being on macrocyclic lactone heartworm preventives , as previously described ( bourguinat et  al . , 2015 ) . the method followed for synchronized file generation has been described elsewhere ( mani et  al . , 2016 ) . d.  immitis nuclear genome v2.2 ( http://nematodes.org/genomes/dirofilaria_immitis/ ) was used as the reference genome . bam files that consisted of aligned sequences from each population against the reference genome were generated . the program popoolation2 was used to generate a synchronized file that contained nucleotide frequencies for each population at every base in the reference genome . this information was obtained for each gene in a concise format after filtering for base quality . previously 8 abc - b transporter genes in d.  immitis had been identified , named and classified based on phylogenetic analysis ( bourguinat et  al . , 2016 ) . for the rest , the naming was either based on the closest homolog in c.  elegans , identified during a blastp search , or as given during annotation of the d.  immitis ndi.2.2.2 nuclear genome . fifteen abc transporter genes / subunits of d.  immitis were classified , based on the classifications of abc transporter homologs in c.  elegans ( sheps et  al . , 2004 ) and b.  malayi ( ardelli et  al . , 2010 ) . a similar methodology , as described by mani et  al . ( 2016 ) was followed to identify the location of a snp within a gff ( generic feature format ) file that contained annotations of the nuclear genome of d.  immitis ndi.2.2.2 was used to call whether a locus was in an exonic or intronic region of the gene . the type of snp , whether synonymous ( ssnps ) or non - synonymous ( nssnps ) ( missense or nonsense ) , and the respective amino acid changes were assessed using sequencher software 4.10.1 ( gene codes corporation , ann arbor , mi , usa ) . transmembrane helices of these membrane proteins were predicted using online web tools such as phobius ( http://phobius.sbc.su.se/ ) and hmmtop ( http://www.enzim.hu/hmmtop/index.php ) . for those of the abc transporter genes whose transmembrane domain and atp - binding cassette domain were previously predicted ( bourguinat et  al . , 2016 ) , the topological structures were displayed using the topo2 transmembrane protein display software ( http://www.sacs.ucsf.edu/cgi-bin/open-topo2.py ) . abc transporters with potential drug target capabilities were highlighted , based on similar filtering methodology used for b.  malayi ( kumar et  al . , 2007 ) and d.  immitis ( godel et  al . , 2012 , mani et  al . , 2016 ) . one such filter was evidence of detrimental effects of its homolog during rnai in c.  elegans . such detrimental effects could be , for example , embryonic lethal / larval arrest , shortened life span , locomotion variant , organism development variant , sluggish / fainter , slow growth , egg size defective , pharyngeal pumping variant , as known from the wormbase website ( http://legacy.wormbase.org/ ) . the other type of filter used was the absence of a blastp hit with an e - value below 10 , in the predicted proteomes of homo sapiens and canis lupus familiaris . to study the impact of identified snps , especially non - synonymous snp ( nssnps ) , on the secondary structure of an abc transporter , a secondary structure prediction algorithm called psipred ( http://bioinf.cs.ucl.ac.uk/psipred/ ) ( mcguffin et  al . a total of 169 nucleotide sequences that could encode abc transporter genes or subunits from subfamilies a , b , c and g were extracted from available databases such as ncbi , nembase4 , wormbase and broad institute . the number of nucleotide sequences , of related nematodes , extracted for analysis is given in table  1 ( complete information of sequences is available in supplementary data file s1 ) . of the total nucleotide sequences extracted , 13 of them belonged to abc - a , 84 to abc - b , 29 to abc - c and 43 to abc - g encoding genes / subunits . by the complementary genomic approaches , we have identified 17 genes or subunits that encode 15 unique abc transporter genes ( see table  2 ) . out of the 15 , eight abc transporter genes have already been described ( bourguinat et  al . , 2016 ) . these genes , identified on the basis of pcr amplification , sequencing and bioinformatics tools , included three pgp encoding genes ( dim - pgp-3 , dim - pgp-10 , dim - pgp-11 ) , a pseudogene and two half - transporter genes ( dim - haf-1 , dim - haf-4 ) from the abc - b subfamily , and two half transporter genes ( dim - haf-5.1 , dim - haf-5.2 ) from the abc - c subfamily ( bourguinat et  al . , 2016 ) . we identified 7 more abc transporter genes , two ( dim - abt-2 , dim - abt-4 ) from sub - family a , three ( dim - mrp-1 , dim - mrp-5 , dim - mrp-7 ) from sub - family c and two ( dim - wht-4 , dim - wht-8 ) abc transporter genes from sub - family g. some of the abc transporter genes identified by bourguinat et  al . ( 2016 ) were different from those of their respective annotations in the d.  immitis genome ( http://nematodes.org/genomes/dirofilaria_immitis/ ) . for example , dim - pgp-3 , dim - pgp-10 , dim - haf-1 , dim - haf-4 and dim - haf-5 were differently annotated in the genome . in such cases , every polymorphic position within those genes has been called , based on the annotations by bourguinat et  al . ( 2016 ) , as the authors had used bioinformatic tools and pcr amplification to confirm each gene . for the rest of the abc transporter genes , the position of every snp was called as per gene annotations in the d.  immitis genome . the blastp search - filter , although applied to help identify nematode - specific abc transporters ( and their associated polymorphisms ) , did not affect any of the results on polymorphism . out of 17 abc transporter genes / subunits that encoded 127,275 base pairs ( bp ) , sequence coverage was only available for 119,033 bp . this number of base pairs , that covered 0.14% of 84.2  mb sized nuclear assembly of d.  immitis , was analyzed for polymorphism . as previously described ( mani et  al . , 2016 ) , a position was called polymorphic if there was a 20% difference between the frequency of 2 or more nucleotides . using this 20% threshold level , we identified 231 snps in 119,033 bp with a nucleotide diversity rate of 0.19% ( see table  3 ) . of the total snps , 67 were present across all of the populations , while 89 and 75 were specific to sus and loe populations , respectively . we identified 193 snps in intron regions , with an intronic snp rate of 1/412 bp , and the remaining 38 snps were present in exons , with snp rate of 1/1041 bp . twenty - one of the snps in exons caused changes in amino acids ( see table  3 ) . polymorphic rates calculated for the pooled sus populations , from different locations , are given in table  4 . three full - length transporter pgps , a pseudogene and 2 half transporters from abc - b subfamily of d.  immitis were investigated for polymorphism . together , they contained 59 snps within 55,588 bps covered with an exonic snp rate of 1/1457 and an intronic snp rate of 1/797 . dim - pgp-3 , dim - pgp-10 , dim - pgp-11 and a pseudogene were identified in scaffolds ndi.2.2.scaf00046 , -48 ,  4 and  41 respectively , as previously described ( bourguinat et  al . , 2016 ) . these authors identified half - transporter genes dim - haf-1 and dim - haf-4 in scaffolds ndi.2.2.scaf00023 and ndi.2.2.scaf00101 , respectively . in dim - pgp-3 a polymorphism at the 645th codon of the gene changed aliphatic amino acid isoleucine to nucleophilic amino acid threonine ( table  5 ) . this snp was predicted to be located at the linker region of the transporter ( see supplementary data file s3 , fig . two synonymous types of polymorphism occurred at amino acid positions 322 and 834 of the gene . dim - pgp-10 was the least polymorphic ( 0.05% ) among the pgp genes with 8 snps . of the two snps localized in exon regions , one at amino acid position 176 caused an amino acid variant with either serine if the corresponding codon was tct or phenylalanine if the corresponding codon was ttt ( table  5 ) . this locus was predicted to be located at the cytoplasmic region between tm2 and tm3 ( see supplementary data file s3 , fig . a second polymorphism ( gctggt ) identified at 317 codon of the protein was identified with a317 g , predicted to be present at the fifth transmembrane region of the protein ( see supplementary data file s3 , fig . dim - pgp-11 ( bourguinat et  al . , 2016 ) previously implicated in loss of response to ivermectin ( bourguinat et  al . , 2011a , bourguinat et  al . , of those , 8 were only detected in sus and 5 were only observed in loe populations . two ( r230p , k1203r ) of five exonic snps in dim - pgp-11 were non - synonymous types of polymorphism ( table  5 ) . the snp k1203r , found previously as an exonic part of the diplotypic  gg - gg  genotype of d.  immitis microfilariae that correlated with loe response to ivermectin ( bourguinat et  al . , 2011b ) , was present predominantly in loe pooled populations as expected ( see supplementary data file s2 ) . two adjacently placed polymorphic sites at the 230 position of the amino acid were known . a second position change at the 230 codon caused a change from arginine to proline whereas the third position nucleotide change did not affect the amino acid coded . an abc transporter - encoding gene identified in scaffold ndi.2.2.scaf00041 was called as a pseudogene as it contained multiple genetic lesions that could prevent translation to a functional protein ( bourguinat et  al . , 2016 ) . though the gene appears to be non - functional , polymorphism in pseudogenes can be important ( macphee et  al . , 2002 , bourguinat et  al . , 2016 ) . polymorphic analysis of the pseudogene allowed us to identify 20 snps , out of which , three were in exonic regions . a deletion type of polymorphism at codon 135 caused a stop codon at the downstream flanking region across the populations . at position 293 of the pseudogene , a polymorphic variant with either an alanine or valine was present ( see supplementary data file s2 ) . among the half abc - b transporter genes , dim - haf-4 had two snps in the intron regions of the gene , one exclusively in the loe populations and the other ( deletion polymorphism ) across all populations . a single intronic polymorphic locus in the dim - haf-3 gene was found to be present only in the sus populations ( see supplementary data file s2 ) . three mrp genes ( dim - pgp-2 , dim - pgp-5 , dim - pgp-7 ) and two half - transporter genes ( dim - haf-5.1 , dim - haf-5.2 ) ( bourguinat et  al . , 2016 ) from the abc - c subfamily , that together covered 44,689 bp of the ndi2.2 d.  immitis genome , were used for snp investigation . with a total of 75 and 20 snps in the intron and exon regions of these genes , respectively , an overall polymorphic rate of 0.2% was determined . a homolog of the mrp-1 gene , identified in scaffold ndi.2.2.scaf00010 , had 38 snps , out of which 22 were specific to sus populations and 8 were specific to loe populations . two of the five exonic snps caused amino acid change , namely t24p and l1185f , whereas the remaining three snps were synonymous types of polymorphism . the snp l1185f caused a change in predicted secondary structure of the protein ( table  6 , fig .  1 ) . dim - mrp-5 had 30 snps with estimated polymorphic rate of 0.25% ; the most polymorphic amongst the mrp genes studied . eleven polymorphic sites were specific to sus populations and another 14 were specific to loe populations . seven out of the eight exonic snps at amino acids positions 411 , 412 , 634 , 638 , 640 , 1230 and 1295 of mrp-5 , were identified to cause amino acid changes ( table  6 ) , out of which , the snp a411 g was predicted to cause a change in secondary structure ( fig .  1 ) . dim - mrp-7 in scaffold ndi.2.2.scaf00004 was identified with 16 snps , four of them were in exonic regions . all four exonic snps were non - synonymous , namely , w258 g , k263n , l280s and g1100c ( table  6 ) . as shown in fig .  1 , two of the snps , k263n and g1100c would cause a change in secondary structure , as predicted by psipred . one out of three in the exon regions was non - synonymous , causing a change in amino acid from serine ( agc ) to arginine ( aac ) at position 61 of the predicted amino acid sequence . this snp is anticipated to be in the cytoplasmic region , before the tm1 ( see supplementary data file s3 , fig . two other snps , at positions 84 and 126 of the amino acid sequence , were synonymous types of polymorphism . two abc - g transporter genes ( dim - wht-4 , dim - wht-8 ) together covered 17,912 bp of the genome . with 70 snps in introns and 5 in exons , they had an overall polymorphic rate of 0.4% . the wht-4 gene , which spanned 12,560 bp in the scaffold ndi.2.2.scaf00589 , had at least 69 snps , with a calculated polymorphic rate of 0.54% , which was the highest among all the abc transporter genes studied . among the snps identified , 33 of them were only in loe populations and 14 of them were specific to sus populations . all of the 5 snps in exons were synonymous , of which 4 occurred only in sus populations . in another wht gene ( dim - wht-8 ) , identified in scaffold ndi.2.2.scaf00252 , all 6 snps discovered were in intron regions ( supplementary data file s2 ) and 5 were specific to sus populations . multidrug transporters such as p - glycoproteins are members of the abc transporter superfamily of efflux transporters and they have been associated with drug resistance in helminths , including nematodes ( lespine et  al . , 2012 , ardelli , 2013 ) and trematodes ( mottier et  al . , 2006 , greenberg , 2013b ) . these transporters may influence anthelmintic drug concentrations at the target sites of the parasite and the co - administration of mdr inhibitors ( e.g. , verapamil , valspodar ) with anthelmintics has enhanced efficacy in  vivo ( lespine et  al . , 2012 ) . in addition , these transporters are likely to play key roles in worm physiological functions such as excretion , reproduction and modulation of host responses ( greenberg , 2014a , greenberg , 2014b ) . therefore , abc transporters may prove to be attractive therapeutic targets for agents that can act as inhibitors , or as anti - nematode agents in their own right . recently , there has been an interest to exploit these transporters for resistance reversal . for resistance mediated by mdr transporters , different strategies for reversal new compounds may either act as mdr inhibitors , directly inhibiting efflux , or as inhibitors of transcriptional regulation that may modulate gene expression at the transcription level or at the level of mrna stability , at the post - transcriptional level ( lespine et  al . , knowledge of the polymorphism of abc transporter genes of d.  immitis , especially from sub - families a , b , c and g , implicated in drug resistance ( ardelli , 2013 ) , may be relevant and provide useful information during transport inhibitor or reversing agent drug design . such knowledge might also be helpful to assess potential resistance development . to study polymorphism , we wished to identify all putative abc transporter genes / subunits of d.  immitis by comparative genomic approaches . we followed this approach as the genome of d.  immitis was not fully annotated ( bourguinat et  al . , 2015 ) . we found that the genome of d.  immitis contains at least 15 abc transporter genes ; 2 abt genes from sub - family a , 3 pgp , 1 pseudogene and 2 half transporters from sub - family b , 3 mrp , 2 half transporters from subfamily c , and 2 white - like genes from subfamily g. ardelli ( 2013 ) predicted 33 abc transporter genes ( including 8 pgp , 8 haf and 5 mrp genes ) in b.  malayi , a closely related filarial nematode ; whereas another filarial nematode , o.  volvulus was reported to have 7 abc transporter genes ( including 1 pgp , 3 haf and 1 mrp genes ) ( ardelli , 2013 ) . the polymorphic rate was calculated for each sus pooled population and it was found that the rate was similar between samples from the usa , grand canary and grenada . however , a low polymorphic rate was determined in samples from italy and this can be explained by the fact that those worms were from a single dog . this hypo - variability of abc transporter genes is in agreement with the finding of low genetic variability among ion channel genes ( mani et  al . , 2016 ) . moreover , polymorphism studied at the microsatellite level ( belanger et  al . , 2011 ) and at the whole genome level ( godel et  al . , 2012 ) also suggested low genetic variability among d.  immitis populations . the latter study , however , used two single d.  immitis worms , one each from italy and the usa to draw this conclusion . this trend of hypovariability would be consistent with recent common ancestry of european and american heartworm populations due to migration . d.  immitis pgps could be of interest as targets for multidrug resistance reversal , based on evidence that overexpression of pgps plays a role in resistance in other nematode species and that resistance mediated by pgp genes in a related nematode can be reversed after the introduction of inhibitors . a significant increase in the expression profiles of cel - pgp-11 , cel - haf-1 , cel - mrp-1 , cel - mrp-5 , cel - mrp-7 after ivm treatment ( yan et  al . , 2012 , ardelli and prichard , 2013 , bygarski et  al . , 2014 ) was reported in an ivm - resistant strain , compared to the wild - type bristol n2 strain . sensitivity to ivermectin was increased by mdr inhibitors in b.  malayi ( tompkins et  al . p - glycoproteins of d.  immitis could also be interesting targets for new anthelmintic drugs , considering the physiological roles their homologs play in other nematodes , including c.  elegans . deletion of some of the pgp genes had a measurable effect on phenotypes of c.  elegans ( janssen et  al . , 2013 ) . for example , c.  elegans deletion strains pgp-3 and pgp-10 were found to have significantly higher levels of pharyngeal pump rate than the wild - type strain , suggesting that these genes affected feeding in these strains . also , the pgp-11 gene is known to be essential for larval development of c.  elegans ( janssen et  al . , 2013 ) . these findings suggest that pgp-34 , pgp-10 , and pgp-11 of c.  elegans have phenotypic effects on resistance , feeding and development of the parasite . if we suppose similar gene product activities for dim - pgp-3 , dim - pgp-10 , dim - pgp-11 , in d.  immitis , then these genes could be of potential interest as targets for candidate drugs . six abc - b transporter genes ( dim - pgp-3 , dim - pgp-10 , dim - pgp-11 , 1 pseudogene , dim - haf-1 , dim - haf-4 ) , previously identified and named ( bourguinat et  al . , 2016 ) were found to have at least 59 snps , an overall polymorphic rate of 0.1% . in dim - pgp-3 , a missense variant ( i645 t ) with an aliphatic amino acid , isoleucine , or a nucleophilic amino acid , threonine , predicted to be in the linker region of the protein ( supplementary data file s3 , fig . it has been shown that the linker region of abc transporters mediates ubiquitination and control of protein turnover ( klling and losko , 1997 , sato et  al . , 2009 ) , two processes implicated in mechanisms of drug resistance in plasmodium ( hunt et  al . , 2007 , the pgp-3 homolog in c.  elegans , cel - pgp-3 , found to be expressed in the intestine and excretory cell ( zhao et  al . , 2004 ) , has a role in feeding of the parasite ( bygarski et  al . , 2014 ) . considering all this , the pgp-3 of d.  immitis could be a potential target for a reversal mechanism . in the dim - pgp-10 gene , 8 snps , including 2 exonic snps , were identified . in c.  elegans , pgp-10 is expressed in the intestine ( zhao et  al . , 2004 ) but its role is not yet known . further studies need to be done to understand the implications of the snps identified in dim - pgp-10 . another pgp gene of d.  immitis , dim - pgp-11 , was found to be very polymorphic as it contained at least 18 snps . of them , five were exonic including one of interest , k1203r , found mainly in loe populations . the same snp ( k1203r ( aagagg ) ) in the dim - pgp-11 gene has been captured in different forms by two separate studies , validating our approach to identify snps in abc transporter genes . it was identified as position 11 that corresponded to the 1st gg of the gg - gg genotype of a partial dim - pgp-11 sequence ( genbank accession number hm596853 ) ( bourguinat et  al . , 2011b ) . in another study , it was known as locus_48,992_b , identified by comparing read frequencies of four sus and four loe d.  immitis populations ( bourguinat et  al . , 2015 ) . the nbds and walker a and walker b sites are very conserved across abc - b proteins , and recently it has been shown ( verhalen et  al . , 2017 ) that changing e1197 to 1197q ; which makes the nbd region just adjacent to the walker b site more polar and more basic , impairs atp turnover in mammalian p - glycoprotein , and thus likely transport function . a similar change occurring with the k1203r snp , found in dim - pgp-11 , similarly introduces a more polar , very basic group adjacent to the walker b site which may modulate atp turnover and transport activity . alignment of this region in a number of organisms indicates that k1203 is highly conserved . thus the mutation to 1203r may represent an important change affecting the functioning of the dim - pgp-11 transporter . in the pgp-11 of d.  immitis , two additional snps at the second and third position of codon 230 were identified , with the former causing an amino acid change from arginine ( basic amino acid ) to proline ( cyclic amino acid ) . this non - synonymous snp was predicted to be present at the cytoplasmic region between tm4 and tm5 , immediately after the tm4 region ( see supplementary data file s3 , fig . the crystal structure of a eukaryotic homolog , cmabcb1 ( from cyanidioschyzon merolae ) revealed that the tm4 helix acts as a  gatekeeper  of the substrate entrance pore and is therefore critical for efficient drug binding and transport ( kodan et  al . , 2014 ) . in the mouse pgp , it was suggested that tm4 forms one of the portals open to the inner leaflet for drug entry and so is believed to play an important role in selecting drugs from the lipid bilayer ( li et  al . , 2014 ) . if we assume a similar scenario , any snp present in and around the tm4 helix of pgp-11 of d.  immitis , may play a role in substrate selection , binding and transport ; therefore the snp r230p in dim - pgp-11 maybe of interest during drug design . another nssnp , k1203r , predicted to be present immediately before the walker b motif ( illde ) in the second nucleotide binding domain of the transporter ( see supplementary data file s3 , fig . s3c ) may be of possible interest for further study as this polymorphic locus was implicated in the loe phenotype to ivm treatment ( bourguinat et  al . , 2011b ) . d.  immitis pgp-11 is only 36% identical to pgp-1 of c.  elegans for which a higher resolution ( 3.4 a ) crystal structure is available ( jin et  al . , 2012 ) and has been shown to bind ml anthelmintics ( david et  al . , 2016 ) and therefore caution is required in using the structural knowledge of cel - pgp-1 to compute functional implications of identified snps in dim - pgp-11 . in c.  elegans , the pgp-11 gene is known to be expressed in the intestine and excretory cell ( zhao et  al . , 2004 ) . moreover , the homolog of pgp-11 is essential for survival in c.  elegans , as the gene is required for larval development ( janssen et  al . , 2013 ) . in d.  immitis , a diplotypic gg - gg genotype of the gene strongly correlated with loss of ivermectin response phenotype of mfs ( bourguinat et  al . , 2011a , bourguinat et  al . , 2011b ) taken together , the pgp-11 of d.  immitis may be a gene of potential interest for drug discovery . their protective role in worms is well known as they were found to reduce the toxicity of heavy metals ( broeks et  al . , 1996 ) and mls ( lespine et  al . , 2012 ) . there is an increased expression of mrps in wild - type ( lespine et  al . , 2012 ) and ivermectin - resistant strains ( james and davey , 2009 ) of c.  elegans after exposure to mls . immitis , namely dim - mrp-1 , dim - mrp-5 and dim - mrp-7 may be interesting drug targets . in the mrp-1 gene , a polymorphism at codon 24 changed the amino acid from threonine ( acc ) to proline ( ccc ) . another snp with either a leucine or phenylalanine at position 1185 of the amino acid sequence was identified . the homolog of mrp-1 , in c.  elegans , regulates dauer larva formation ( yabe et  al . , 2005 ) in addition to its protective role against heavy metals ( broeks et  al . , 1996 ) . therefore dim - mrp-1 ( and its associated snp ) may be an interesting target for drug design . in c.  elegans , the likely role of the mrp-5 protein , as a heme exporter , was shown and targeted deletion of mrp-5 in the intestine caused embryonic lethality in the mutant worms ( korolnek et  al . , 2014 ) . mrp-5 of c.  elegans is thus essential for viability and therefore could be an interesting drug target . however , dim - mrp-5 was found to be highly polymorphic , with at least 30 snps ( eight of them in the exon regions ) and so there is a higher chance of resistance development against a drug targeting this gene . in dim - mrp-7 of d.  immitis , 4 snps were identified in exonic regions , all of which resulted in non - synonymous types of polymorphism . some of the snps identified within mrp genes , namely l1185f ( mrp-1 ) , a411 g ( mrp-5 ) , k263n and g1100c ( mrp-7 ) are predicted to cause a change in secondary structure ( see fig .  1 ) . those snps might have interesting effects on their respective protein function and therefore further studies are needed to understand their implications . in the dim - mrp-7 ( abcc10 ) homolog of humans , two snps ( one in intron 4 and another in exon 12 ) were associated with renal tubular dysfunction ( kathawala et  al . , we believe that some of the snps identified , including snps in non - coding sequence , may have effects on gene expression , mrna stability , function , substrate specificity , or drug resistance . we are aware that some polymorphic sites may have been missed in this study , due to lack of coverage during sequencing , small sample size and limited sample diversity , or due to the use of the 20% threshold value for snp calling . further studies need to be done to confirm each snp in the field and to understand its implications . information on the degree of heterogeneity of d.  immitis has been limited and this study adds to that information . the study lists all genes detected from a , b , c and g sub - families of abc transporters . the merit of the study lies in the diverse large number of samples ( 122 adult worms and 32,000 mfs ) used for snp genotyping , with 231 snps identified that may be of relevance for mdr inhibitor or antinematodal drug design . this study also provides an updated inventory of abc transporter genes in d.  immitis .
dirofilaria immitis , a filarial nematode , causes dirofilariasis in dogs , cats and occasionally in humans . 
 prevention of the disease has been mainly by monthly use of the macrocyclic lactone ( ml ) endectocides during the mosquito transmission season . 
 recently , ml resistance has been confirmed in d.  immitis and therefore , there is a need to find new classes of anthelmintics . 
 one of the mechanisms associated with ml resistance in nematodes has been the possible role of atp binding cassette ( abc ) transporters in reducing drug concentrations at receptor sites . 
 abc transporters , mainly from sub - families b , c and g , may contribute to multidrug resistance ( mdr ) by active efflux of drugs out of the cell . 
 gene products of abc transporters may thus serve as the targets for agents that may modulate susceptibility to drugs , by inhibiting drug transport . 
 abc transporters are believed to be involved in a variety of physiological functions critical to the parasite , such as sterol transport , and therefore may also serve as the target for drugs that can act as anthelmintics on their own . 
 knowledge of polymorphism in these abc transporter genes in nematode parasites could provide useful information for the process of drug design . 
 we have identified 15 abc transporter genes from sub - families a , b , c and g , in d.  immitis , by comparative genomic approaches and analyzed them for polymorphism . 
 whole genome sequencing data from four ml susceptible ( sus ) and four loss of efficacy ( loe ) pooled populations were used for single nucleotide polymorphism ( snp ) genotyping . out of 231 snps 
 identified in those 15 abc transporter genes , 89 and 75 of them were specific to the sus or loe populations , respectively . 
 a few of the snps identified may affect gene expression , protein function , substrate specificity or resistance development and may be useful for transporter inhibitor / anthelmintic drug design , or in order to anticipate resistance development .