the sacculated forestomachs of kangaroos and wallabies ( family macropodidae ) commonly contain large numbers of nematodes belonging to the strongylid sub - family cloacininae .
for example , vendl and beveridge ( 2014 ) reported means and ranges in numbers of cloacinine nematodes in the stomachs of the red - necked wallaby , macropus rufogriseus , the eastern grey kangaroo , macropus giganteus and the swamp wallaby , wallabia bicolor , as 60,800 ( 2000210,000 ) , 20,500 ( 700079,000 ) and 20,000 ( 300058 , 000 ) , respectively .
these data support earlier published figures for high intensities of infection in the western grey kangaroo , macropus fuliginosus , the red kangaroo , macropus rufus , and m.  giganteus ( arundel et  al .
although there is considerable species diversity in the stomach - inhabiting cloacinine nematodes ( spratt et  al . , 1991 ) ,
for instance , rugopharynx australis dominates the gastric helminth community of m.  rufus , with mean and maximum burdens of 47,000 and 266,000 nematodes , respectively ( arundel et  al . , 1979 ) .
in spite of the numerical significance of this genus in the gastric helminth communities of macropodids , continuing taxonomic studies are needed , as it is unlikely that all species have yet been described and the presence of cryptic species could potentially complicate the interpretation of ecological data published to date .
the genus rugopharynx was revised by beveridge ( 1982 ) , who recognised nine species .
however , beveridge ( 1982 ) noted that r.  australis was potentially a complex of a number of species , differentiable only by very minor and often overlapping morphological characteristics .
multilocus enzyme electrophoretic ( mee ) studies , combined with morphological evidence , indicated the existence of two new species , r.  sigma and rugopharynx mawsonae , both formerly confused with rugopharynx zeta ( chilton et  al .
1994 ) , while an additional mee investigation of r.  australis by chilton et  al .
subsequently , beveridge and chilton ( 1999 ) split r.  australis into 10 species and resurrected r.  alpha as a valid species .
the latter revision was based on extrapolating the minor morphological differences identified in samples included in the earlier electrophoretic study across the entire species complex . in spite of this progress , there has been no attempt to independently verify the validity of the species erected to date using molecular methods .
a closely related genus , rugonema was erected by beveridge ( 1999 ) for specimens formerly referred to as r.  australis occurring in the stomach of the black - gloved wallaby , macropus irma , from western australia , based on morphological differences in the labial collar .
because of this close association with rugopharynx , rugonema labiatum was also included in the present study . given the prevalence and abundance of this nematode genus in the stomachs of kangaroos and wallabies , and the diversity of species currently recognised based on minor morphological criteria , this study was undertaken to attempt to establish species boundaries within the genus based on sequences of the first and second internal transcribed spacers ( its-1 and its-2 ) of nuclear ribosomal dna .
these molecular target regions have proved to be highly informative for the specific identification of a range of strongylid nematodes , including taxa within the cloacininae ( chilton , 2004 ) , the subfamily to which rugopharynx belongs .
nematodes were obtained from the stomachs of a range of kangaroos and wallabies ( fig .  1 ; table  1 ) , which had been collected as fresh road - kills or from road - kills frozen prior to examination .
host nomenclature follows van dyck and strahan ( 2008 ) . in instances where a nematode species occurred across a large geographical area
, an attempt was made to include samples from different geographical regions of australia , particularly any occurring on the island of tasmania ( fig .  1 ) .
australian state names are abbreviated as : nsw , new south wales ; qld , queensland ; sa , south australia ; tas , tasmania ; vic , victoria ; wa , western australia .
nematodes were washed in saline , frozen in liquid nitrogen and stored at  80  c until morphological and molecular studies were undertaken .
nematodes were then thawed , the head and tail of each worm removed , fixed in lactophenol and mounted permanently in polyvinyl lactophenol as voucher specimens .
nematodes were identified according to previous descriptions ( beveridge , 1982 , beveridge et  al . , 1994 , chilton et  al . , 1993 , beveridge and chilton , 1999 , appan et  al . , 2004 ) .
voucher specimens have been deposited in the south australian museum ( sam ) , adelaide ( table  1 ) .
species within the genus rugopharynx were divided into three groups based on the morphology of the buccal capsule ( table  1 ) .
nematodes with a simple , cylindrical buccal capsule were designated as type i ( fig . 
2 ) , while those with a buccal capsule divided into two sections were designated as type ii .
these buccal capsules were further subdivided into species with a buccal capsule divided in the mid region ( r.  epsilon ) ( type iia ) and those in which the division occurred in the anterior quarter ( r.  rufogrisea ) ( type iib ) .
species with a buccal capsule divided into three segments were designated as type iii . in some species , such as r.  spratti and r.  tau ,
the potential division of the buccal capsule is subtle ; hence , in these cases , species were classified as type i. genomic dna was isolated from the mid - body part of each nematode using a small - scale sodium - dodecyl - sulphate / proteinase k extraction method ( gasser et  al . , 1993 ) , followed by purification
the region of rdna comprising the its-1 , 5.8s rrna gene , its-2 and flanking sequences (= its+ ) was amplified by the pcr using primers nc16 ( forward ; 5-agttcaatcgcaatggctt-3 ) and nc2 ( reverse ; 5-ttagtttcttttcctccgct-3 ) .
pcr was performed in a 50  l volume for 30 cycles at 94  c for 30  s ( denaturation ) , 55  c for 30  s ( annealing ) and 72  c for 30  s ( extension ) , followed by one cycle at 72  c for 5  min ( final extension ) .
amplicons were purified using mini - columns ( using wizard pcr - preps , promega ) , and the its+ sequenced in both directions using the same primers ( separately ) as used for pcr .
the sequences generated in the present study have been deposited in the genbank database ( accession numbers ln906946-ln906995 ; table  1 ) .
sequences were initially aligned using the program muscle ( edgar , 2004 ) , and alignments adjusted manually using the program mesquite v.2.75 ( maddison and maddison , 2011 ) .
phylogenetic analyses of the aligned sequence data were conducted by bayesian inference ( bi ) using monte carlo markov chain analysis in the program mrbayes v.3.2.2 ( ronquist and huelsenbeck , 2003 ) .
the likelihood parameters set for the bi analysis of sequence data were based on the akaike information criteria test in jmodeltest v.2.1.5 ( posada , 2008 ) .
the number of substitutions was set at 6 , with a gamma - distribution . for the tree ,
posterior probability ( pp ) values were calculated by running 2,000,000 generations with four simultaneous tree - building chains .
the standard deviation of split frequencies was < 0.01 , and the potential scale reduction factor approached one .
for each analysis , a 50%-majority rule consensus tree was constructed based on the final 75% of trees produced by bi .
phylogenetic analyses of the its+ sequence data were also conducted using the neighbor - joining ( nj ) and maximum parsimony ( mp ) methods in paup * 4.0b10 ( swofford , 1999 ) . for the mp analyses ,
heuristic searches were carried out with random addition of sequences ( n  =  100 ) , tree - bisection - reconstruction ( tbr ) branch swapping , the multrees option in effect , maxtrees set at 2000 and saving all equally parsimonious trees .
the tree length ( l ) , consistency index , excluding uninformative characters ( ci ) and the retention index ( ri ) were recorded .
bootstrap analyses ( 1000 replicates ) were conducted to determine the relative support for clades in the consensus trees ; nodal support was expressed as a percentage . the its+ sequences of labiostrongylus australis ( genbank accession numbers aj308403-aj308411 ) were used as the outgroup in the phylogenetic analyses , because this nematode species belongs to a related tribe within the same subfamily ( i.e. cloacininae ) ( see lichtenfels , 1980 ) to which rugopharynx also belongs ( beveridge , 1982 , chilton et  al . , 1997 ) and because it has been relatively well studied genetically ( chilton et  al . , 2009b ) .
a consensus tree depicting the phylogenetic relationships of the nematodes was compared with a composite tree of the hosts .
the host phylogenetic tree was based primarily on the molecular analyses of meredith et  al .
( 2008 ) ; however , as some species ( e.g. , thylogale billardierii , macropus dorsalis , petrogale herberti , petrogale assimilis and petrogale inornata ) were not included in that study , the tree was modified to incorporate these hosts based on the studies of campeau - ploquin et  al .
, 2001 , cardillo et  al . , 2004 and eldridge et  al .
m. ( osphranter ) bernardus and m. ( o. ) antilopinus , from which no species of rugopharynx have been reported .
not all species of petrogale and thylogale are included in the phylogenetic tree , only those included in the present study .
the length of the its+ region , excluding flanking regions was 744777  bp for all taxa within the genus rugopharynx and rugonema labiatum . the length of the its-1 sequences ranged from 373  bp ( rugopharynx mawsonae , r.  sigma from thylogale thetis and rugonema labiatum ) to 387  bp ( r.  mu , r.  rufogrisea and some individuals of r.  australis ) , whereas the its-2 sequences were much shorter , ranging from 216  bp ( r. longibursaris , r. omega and r.  zeta ) to 313  bp ( r.  alpha and r.  mawsonae ) ( supplementary table  1 ) .
the length ( 153  bp ) and nucleotide sequence of the 5.8s rrna gene was the same for all morphospecies within the genus rugopharynx and rugonema labiatum .
there were five fixed differences ( i.e. where there are no shared nucleotides at an alignment position ) between these two taxa ; two in the its-1 and three in the its-2 .
the magnitude of fixed differences in its+ sequence among morphospecies within the genus rugopharynx ranged from eight ( i.e. between r.  omega and r.  longibursaris ) to 84 ( i.e. between r.  alpha and r.  rufogrisea ) .
only one fixed difference in its-1 sequence was detected between r.  setonicis and rugopharynx rho compared with 10 fixed differences between these two morphospecies for the its-2 , while the lowest number of fixed differences in its-2 sequence among morphospecies ( i.e. five ) was between r.  omega and r.  longibursaris .
for some morphospecies , there was variation among individuals in the dna sequences of the its-1 and/or its-2 ( see supplementary tables  28 ) .
three variable positions ( two in the its-1 and one in the its-2 ) were detected in the aligned sequences of r.  mu individuals from hosts ( w.  bicolor ) collected in new south wales and victoria ; however , none of these mutations represented a fixed difference ( supplementary table  2 ) .
similarly , there were no fixed differences in its+ sequence among three specimens of r.  mawsonae from the same host ( m.  dorsalis ) , or among two specimens of r.  pi from two host species ( m.  rufogriseus and m. parryi ) , even though variations were detected at nine ( four in the its-1 and five in the its-2 ) and three ( two in its-1 and one in the its-2 ) alignment positions , respectively . a similar pattern of intraspecific variation was found for three specimens of r.  macropodis from two host species ( m.  giganteus and m.  fuliginosus ) , except for a single nucleotide insertion in the its-1 sequence for one of two individuals ( i.e. f736 ) collected from m. giganteus .
specimens of r.  longibursaris from m.  rufogriseus collected from launceston ( tasmania ) and emu flat ( new south wales ) differed in its+ sequence at eight alignment positions ( five in the its-1 and three in the its-2 ) , whereas no genetic variation occurred in the dna sequence of r.  spratti , specimens of which were collected from the same host species and localities as r.  longiburaris . similarly , there were no fixed differences in its+ sequence of five r.  epsilon individuals collected on the mainland ( i.e. from m.  rufogriseus and w.  bicolor ) and in tasmania ( from m.  rufogriseus ) , except for a single insertion in the its-1 sequence of specimen f77 collected in new south wales ( supplementary table  3 ) .
in contrast , fixed differences in both the its-1 and its-2 sequences were detected among individuals of four morphospecies , r.  rufogrisea , r.  australis , r.  zeta and r.  rho , collected from different host species and/or geographical regions ( see supplementary tables  48 , respectively ) . in the case of r.  rufogrisea , although there were no fixed nucleotide differences in its+ sequences between specimens from m.  rufogriseus collected in new south wales ( f84 ) and tasmania ( f720 ) , there were seven ( of 12 ) variable nucleotide positions ( two in the its-1 and five in the its-2 ) when compared to a specimen from m.  parryi collected in queensland .
similarly , there were 52 variable nucleotide positions in its+ between r.  sigma specimens from thylogale stigmatica and those from t.  thetis , most ( 25 in its-1 and 24 in its-2 ) representing putative fixed differences between nematodes from the two host species ; 20 of these differences ( five in its-1 and 15 in its-2 ) represented 1 - 10  bp nucleotides .
in addition , nucleotide variation was detected at 20 positions in its+ sequences among individuals of r.  zeta collected from different species of petrogale in queensland .
there was only one fixed nucleotide difference in its-1 sequence between specimens of r.  zeta and those from p.  herberti and p.  inornata , whereas these nematodes differed unequivocally at 13 ( 11 in its-1 and two in its-2 ) of 20 variable positions when compared with the its+ sequence from a specimen of r.  zeta from p.  assimilis .
genetic variation was also detected at 52 nucleotide positions in the its+ sequences among seven individuals of r.  australis .
this morphospecies could be separated into two groups ( clades ) based on their its+ sequences .
nematodes in clade 1 had fixed differences at 38 ( 14 in its-1 and 24 in its-2 ) alignment positions when compared to those in clade 2 .
one eastern grey kangaroo ( m.  giganteus ) collected from new south wales contained taxa from both clades 1 and 2 ( i.e. specimens f751 and f754 ) .
likewise , in the case of r.  rho , individuals from macropus eugenii from south australia ( f724 ) and western australia ( f905 ) had almost identical its+ sequences , which differed at 16 positions ( eight in its-1 and eight in its-2 ) from specimens from m.  fuliginosus from south australia ( f780 ) and western australia ( f910 ) , and from m.  irma ( g114 ) from western australia .
eighteen ( five in its-1 and 13 in its-2 ) of these differences represented 14 nucleotides .
phylogenetic analyses were conducted to determine whether individual morphospecies from different host species and/or geographical regions represented monophyletic assemblages .
the its+ sequences were aligned over 824 positions ( 404 for its-1 , 153 for 5.8s rdna and 267 for its-2 ) , 175 of which were informative for the mp analysis .
the topology of the strict consensus tree of the mp analysis ( not shown ) , based on 1217 equally most parsimonious trees ( l  =  579 , ci  =  0.58 , and ri  =  0.83 ) , was very similar to that produced from the bi ( fig .  3 ) and the nj analyses ( fig .  4 ) . in all phylogenetic analyses ,
there was no support for all specimens of r.  australis representing a monophyletic assemblage .
however , there was absolute support ( pp  =  1.0 in the bi analysis and bs  =  100% in both the nj and mp trees ) for the separation of r.  australis into two clades .
similarly , there was no support for r.  sigma from the two host species ( t.  stigmatica and t.  thetis ) forming a monophyletic assemblage .
there was strong support ( pp  =  0.999 and 1 , and bs  =  95100% ) for the separation of r.  rho into two clades , one containing individuals from m.  eugenii and the other containing individuals from m.  fuliginosus and m.  irma . however , there was support in the bi analysis ( pp  =  0.994 ) for a sister taxa relationship for the two clades of r.  rho , and absolute support ( pp  =  1 ) for a sister relationship of the two clades with r.  setonicis .
in contrast , in the mp and nj analyses , there was no support for a sister taxon relationship between the two clades of r.  rho .
however , there was strong support ( bs  =  90% ) in the nj analysis for an assemblage comprising r.  setonicis and the two clades of r.  rho . there was also strong support ( pp  =  1.0 in the bi analysis and bs  =  9697% in the mp and nj analyses ) for r.  zeta from the three host species , with support ( bs  =  7095% , pp  =  0.831 ) for r.  zeta from p.  inornata and p.  herberti forming a clade to the exclusion of r.  zeta from p.  assimilis .
similarly , there was absolute support ( pp  =  1.0 ; bs  =  100% ) for r.  rufogrisea representing a monophyletic clade , and support ( bs  =  85100% ; pp  =  0.882 ) for r. 
rufogrisea collected from m.  rufogriseus in tasmania and new south wales forming a clade to the exclusion of r.  rufogrisea collected from m.  parryi in queensland . in the phylogenetic analysis of its+ sequence data for the morphospecies ,
there was support in the bi and mp analyses ( i.e. pp  =  0.973 ; bs  =  86% ) for a clade consisting of r.  macropodis , r.  rosemariae , r.  rufogrisea and r.  theta forming a clade with respect to the other species .
there was also some support ( i.e. pp  =  0.995 ; bs  =  7175% ) for a clade containing r.  longibursaris , r.  omega , r.  tau and r.  spratii ( fig . 
comparison of the phylogeny of the nematodes , derived from the analysis of its+ sequence data , with that currently available for the host species ( fig .  5 )
this comparison revealed that a simple buccal capsule ( i ) appears to be the plesiomorphic state ( as for r.  alpha ) .
type i buccal capsules occurred in all clades , being the exclusive buccal capsule type in one clade , but was mixed with type ii capsules in a second clade and mixed with type iii capsules in the third clade .
to date , species of the genus rugopharynx found in the stomachs of macropodid marsupials have been identified solely on the basis of morphological criteria ( beveridge , 1982 ) or the combination of morphological and mee data ( chilton et  al . , 1993 ,
1994 ) . in describing new species within the r.  australis complex , beveridge and chilton ( 1999 ) relied on limited mee data ( chilton et  al . ,
1996 ) and extrapolated from this base in describing nine new species based on small but potentially significant morphological characters .
the present study represents the first detailed examination of the genus using dna sequence data and therefore represents the first test of the validity of the range of species currently recognised either exclusively on morphological grounds or on the basis of morphological and mee data .
the current molecular analyses included most of the currently known species , apart from r.  longispicularis and r.  petrogale . both of these species occur in hosts , such as the parma wallaby , macropus parma , and the brush - tailed rock wallaby , petrogale penicillata , which are currently considered to be rare or vulnerable species ( maynes , 2008 , eldridge and close , 2008 )
; therefore , it is challenging to obtain parasite material from these host species for molecular studies .
no species of rugopharynx are known to occur in m. ( n. ) agilis , m.(o . ) antilopinus and m. ( o. ) bernardus ( spratt et  al . , 1991 ) . in the current study ,
all of the species of rugopharynx presently recognised based on morphological differences ( i.e. morphospecies ) , in some cases with supporting mee data , had a unique set of its+ sequences , thereby providing additional evidence in support of their validity .
fixed differences in the its-1 and its-2 sequences between or among morphospecies were limited ( e.g. , 0.3% between r.  setonicis and r.  rho , and 2.3% between r.  omega and r.  longibursaris , for the its-1 and its-2 , respectively ) , and data from a single pair of ribosomal dna spacers may not always provide unequivocal support for specific status ( nadler and prez - ponce de lon , 2011 ) . however , given the reliability of its-1 and its-2 sequences for identifying and distinguishing nematode species from macropodids to date ( e.g. , chilton et  al .
, 2009a , chilton et  al . , 2012 ) , the evidence presented here is relatively strong . in some instances ,
the morphological and genetic differentiation is supported by the occurrence of formerly cryptic species in the same host individual . in the case of r.  australis and r.  macropodis (
i.e. previously included within r.  australis ) , the former occurring in kangaroos in arid environments and the latter in areas of higher rainfall ( beveridge and chilton , 1999 ) , both species were found at one intermediate location in victoria ( pine plains station ) ( beveridge and chilton , 1999 ) , indicating genetic isolation and providing further support for the validity of the two species .
the data presented here suggest that additional cryptic ( i.e. genetically distinct but morphologically similar ) species of rugopharynx remain to be described
. specimens of r.  australis occurred in two quite distinct clades , although some specimens from the two clades were collected from the same individual host ( m.  giganteus ) . the magnitude of fixed differences in its-1 and its-2 sequences between members of the two r.  australis clades ( 3.6% and 10.4% , respectively ) was greater than that between related morphospecies ( e.g. , 0.8% and 2.3% , respectively , between r.  omega and r.  longibursaris ) .
furthermore , there was no support for the two r.  australis clades forming a monophyletic assemblage .
examination of the voucher specimens involved in this study suggested that the two specimens were differentiable based on spicule lengths , with those from m.  giganteus from trangie , nsw ( f751 ) being 2.0  mm , and those of the additional specimen from the same individual host ( f754 ) being only 1.28  mm long ; the spicules of a similar specimen from m.  fuliginosus from hattah lakes , victoria ( f784 ) , were 1.36  mm long .
beveridge and chilton ( 1999 ) gave the spicule lengths of r.  australis as 1.441.95  mm , which virtually encompasses the range of the specimens used in this study .
the spicule lengths of the most ( morphologically ) similar species to r.  australis , ( i.e. r.  macropodis ) are 1.141.23  mm ( beveridge and chilton , 1999 ) .
in addition , beveridge and chilton ( 1999 ) noted significant differences in bursal morphology within this species .
therefore , it appears that r.  australis , as currently defined , is a composite of at least two species . in the case of r. 
sigma , there was limited variation ( i.e. two fixed differences in its-2 but none in its-1 ) in the its+ sequences of two specimens from t.  stigmatica , collected 1200  km from one another in queensland , whereas they had 49 fixed differences ( 25 in its-1 and 24 in its-2 ) when compared with r.  sigma from t.  thetis .
this magnitude of sequence difference ( 6.5% and 10.3% for its-1 and its-2 , respectively ) exceeded that among many morphospecies within the genus .
furthermore , specimens f794 from t.  thetis and g384 from t.  stigmatica were collected at the same locality ( lamington national park , qld ) , thus being in sympatry .
the phylogenetic analyses also showed that r.  sigma from the two host species did not form a monophyletic clade , providing additional support that they represent cryptic species .
examination of the female voucher specimen of r.  sigma from t.  thetis ( f794 ) indicates a tail length of 0.21  mm compared with 0.390.45  mm for specimens from t.  stigmatica and the distance of the vulva from the posterior end as 0.30  mm compared with 0.600.70  mm in specimens from t.  stigmatica ( see chilton et  al . , 1993 ) .
hence , there appear to be morphological features supporting the molecular differences for these specimens .
( 2000 ) concluded that the helminth communities of these two host species in southern queensland , where they are sympatric , were essentially similar .
specimens of r.  zeta from p.  inornata and p.  herberti formed a strongly supported clade to the exclusion of specimens from p.  assimilis .
these three closely related species of rock wallabies have parapatric distributions along the eastern coast of queensland ( potter et  al . , 2012 ) .
( 2009a ) examined three species of cloacinine nematodes , cloacina caenis , c.  pearsoni and c.  robertsi , which occur in these related rock wallaby species and demonstrated genetic differences between them , suggesting the existence of cryptic species .
there was one fixed difference in its+ between specimens of r.  zeta from p.  inornata and p.  herberti ( but none in the its-2 ) , whereas they had 13 fixed differences ( 11 in its-1 and 2 in its-1 ) when compared to the r.  zeta from p.  assimilis .
the magnitude of sequence difference in its-1 ( 2.8% ) exceeds that detected between r.  omega and r.  longibursaris , whereas the 0.9% sequence differences in the its-2 is less than between these two morphospecies .
therefore , additional molecular investigations are required to test the hypothesis that r.  zeta represents a species complex .
the current data suggest the existence of one species in p.  assimilis and a second species in both p.  inornata and p.  herberti .
rugopharynx zeta also occurs in a number of related species of rock wallabies ( petrogale mareeba , petrogale sharmani and p.  penicillata ) ( spratt et  al . ,
1991 ) and specimens from these additional hosts would need to be included in future studies .
rugopharynx rho from m.  fuliginosus from both south australia and western australia , together with m.  irma from western australia , formed a clade distinct from the same morphospecies obtained from m.  eugenii in south australia and western australia .
there was also no support in the mp and nj analyses for these two clades forming a monophyletic assemblage .
as m.  fuliginosus and m.  eugenii are sympatric at both localities ( van dyck and strahan , 2008 ) , the data suggest that cryptic species may also exist within this taxon . in western australia ,
m.  irma and m.  fuliginosus are sympatric ( van dyck and strahan , 2008 ) and the occurrence of this species in m.  irma may have resulted from host switching . consequently the data presented here suggest that additional cryptic species may exist within r.  australis , r.  rho , r. sigma and r.  zeta .
there were no fixed differences in its+ sequence between specimens of r.  rufogrisea from m.  rufogriseus collected on the mainland of australia and the island state of tasmania , which is consistent with the findings for other morphospecies ( e.g. , r.  epsilon and r.  spratti ) that parasitise m.  rufogriseus .
however , the magnitude of the fixed sequence differences between r.  rufogrisea from m.  rufogriseus and m.  parryi ( 0.5% and 2.1% in its-1 and its-2 , respectively ) is very similar to that between r.  omega and r.  longibursaris , two other species that parasitise m.  rufogriseus ( beveridge , 1982 ) .
the current analysis also provides some insight into host specificity within the genus , although many of the species examined in this study appear to be moderately host specific , occurring in one or two host species ( i.e. r.  alpha , r.  chi , r.  delta , r.  longibursaris , r.  mawsonae , r.  mu , r.  omega , r.  pi , r.  spratti , r.  tau and r.  theta ) , other species included in this study appear to have a wide host range .
thus , r.  australis was identified in m.  rufus ( the type host for the species ) , as well as in m.  fuliginosus , m.  giganteus , macropus robustus and m.  dorsalis .
specimens from the first four host species were collected in arid or semiarid regions of australia ( table  1 ) where these host species are sympatric and r.  australis is a common parasite in each of them ( arundel et  al . , 1979 , beveridge et  al . , 1998 ) .
however , r.  australis is an uncommon parasite of m.  dorsalis ( see beveridge et  al . , 1998 )
and the specimen collected here was in an area in which m.  dorsalis is sympatric with kangaroo species commonly parasitized by this nematode .
analyses of its+ sequence data of r.  epsilon specimens from different host species ( i.e. m.  rufogriseus and w.  bicolor ) suggest that it represents a single species with a broad host range , given that r.  epsilon parasitizes a variety of macropodid hosts ( spratt et  al . , 1991 )
. therefore , examination of additional specimens of r.  epsilon from other host species needs to be studied to test this proposal .
the genus rugonema was erected ( beveridge , 1999 ) for a single species of nematode , ru .
labiatum , from the stomach of m.  irma which resembled rugopharynx but in which the labial collar , instead of forming an annulus , was prominently four - lobed , similar to the genus wallabinema . however , wallabinema lacks a striated buccal capsule and differs in the morphology of the oesophagus .
the sole distinguishing morphological feature of rugonema , that is the four lip - like lobes of the labial collar , is an autapomorphy within the tribe pharyngostrongylinea .
based on the molecular data and a reconsideration of its morphological differentiation , rugonema is here made a synonym of rugopharynx with its sole species becoming rugopharynx labiatum ( beveridge , 1999 ) n. comb .
this species was most genetically similar to r.  pi and belonged to a clade that also included r.  mu .
beveridge ( 1982 ) divided rugopharynx into three groups based on the morphology of the buccal capsule , with either a simple cylindrical buccal capsule , a bilobed or a trilobed buccal capsule .
the group with bilobed buccal capsules consisted of r.  epsilon and r.  rufogrisea , with the indentation occurring in the mid region in r.  epsilon and in the anterior quarter in r.  rufogrisea .
however , among the new species described by beveridge and chilton ( 1999 ) , some were difficult to allocate to a particular group ( i.e. , r.  tau and r.  petrogale ) because the division of the buccal capsule was subtle ( r.  tau ) or variable ( r.  petrogale ) .
the mapping of the morphological data on to the consensus nematode phylogenetic tree ( fig . 
5 ) suggest that the simple buccal capsule is the plesiomorphic state within the genus and that bilobed buccal capsules have evolved independently .
the findings of the present study also suggest that the trilobed buccal capsules are a derived character .
the lack of resolution in the cladogram prevents more detailed conclusions from being drawn on the evolution of buccal capsule shapes within the genus .
there appears to be no co - evolutionary relationship between nematodes and hosts . based on molecular evidence ,
setonix diverged within the macropodine lineage about 10 million years ago , while thylogale and petrogale are sister taxa to the clade that contains macropus and wallabia , their estimated time of divergence being about eight million years ( meredith et  al . , 2008 ) .
the largest clade of the molecular phylogenetic tree of the nematodes contains species from each of these macropodine genera apart from setonix .
beveridge and chilton ( 2001 ) undertook a morphological phylogenetic analysis of the r.  australis complex , which produced a completely unresolved tree , and consequently these authors concluded that there was no obvious co - evolutionary association with hosts .
comparable studies of other cloacinine genera also supported the hypothesis that evolution within this group of nematodes was primarily by host switching ( beveridge and chilton , 2001 ) , a hypothesis concordant with the data presented above .
( 2016 ) have provided molecular evidence for host switching in the related cloacinine genus cyclostrongylus . in summary ,
the molecular data presented here support the earlier morphological studies of the genus rugopharynx and provide additional evidence that the species of rugopharynx currently established , many of them based on minor morphological differences and mee data , are indeed valid .
the study has also revealed the existence of additional cryptic species within the genus that need to be characterised morphologically .
comparisons of buccal capsule morphology with the phylogenetic tree provided some insights into the evolution of more complex buccal capsules , but there were no obvious co - evolutionary associations with hosts , the data instead suggesting a pattern of host switching in the evolution of the genus .
 sequences of the internal transcribed spacers of nuclear ribosomal dna ( its-1 and its-2 ) were determined for species of the genus rugopharynx and rugonema labiatum , nematodes from the stomachs of macropodid marsupials . 
 phylogenetic analyses of the aligned sequence data were conducted . 
 the relationships provided molecular support for all species currently recognised , some of which are based on minor morphological differences and on multilocus enzyme electrophoretic data , but also indicated that additional , cryptic species exist within the genus . 
 in addition , the genus rugonema is placed as a synonym of rugopharynx , its sole species becoming rugopharynx labiatum n. comb . 
 the molecular data provided some insights into the evolution of complex buccal capsule morphologies within the genus , but there was no evidence of co - evolution between the macropodid hosts and their parasites . 
Introduction
Materials and methods
Results
Discussion
Conflicts of interest