avian haemosporidians belong to a diverse group of vector - borne , protozoan blood parasites that infect a broad range of avian species throughout north america and around the world ( valkinas , 2005 ) .
representative species from the genera haemoproteus , plasmodium , and leucocytozoon have been detected on every continent except antarctica ( beadell et  al . , 2009 ) , infecting both wild and domestic birds .
infections by these parasites can vary among taxa and across habitat types ( greiner et  al .
, 1975 , white et  al . , 1978 , durrant et  al . ,
2006 ) , and prevalence in avian host populations may be determined , at least in part , by exposure to viable dipteran vectors ( bennett et  al . , 1992 ) .
other factors such as age , host suitability , and a multitude of environmental variables may also be important determinants of parasite infection ( martnez - abran et  al . , 2004 ) .
rapid warming in northern regions has raised concerns about environmentally - driven changes in the distribution and prevalence of disease in birds and other wildlife ( van hemert et  al . , 2014 ) .
many vector - borne parasites , including plasmodium and other haemosporidians , are temperature - sensitive and are projected to expand in response to climate warming in the arctic ( loiseau et  al .
heavy infection by haemosporidian parasites can cause mild to severe pathogenic effects in various avian species ( ots and hrak , 1998 , palinauskas et  al .
, 2008 ) and exposure of nave hosts to blood parasites may result in dramatic population - level consequences , such as have occurred among native birds in hawaii ( lapointe et  al . , 2012 ) .
recent studies have proposed an increase in haemosporidian infections , particularly plasmodium , among alaskan bird populations as a result of shifts in temperature , vegetation cover , and vector populations ( loiseau et  al . , 2012 ,
however , limited baseline data currently exist for many avian species , particularly year - round residents , making it difficult to assess potential risks and evaluate future patterns of change .
species of birds from the order galliformes ( family phasianidae ) have adapted to the cold climates of arctic and sub - arctic regions and inhabit a broad range of habitats across the northern hemisphere ( aldrich , 1963 , braun and willers , 1967 ) .
given their status as year - round residents , grouse and ptarmigan ( subfamily tetraonidae ) provide excellent model species for studying the prevalence and diversity of locally - transmitted blood parasites in alaska as well as the relationships between climate variables and parasite infection .
four species of grouse and three species of ptarmigan occur throughout alaska 's many biomes and each plays an important role as a prey species for predators and game species for subsistence and recreational hunters .
all species of tetraonids in alaska are year - round residents and many of their populations overlap geographically , but occupy vastly different habitats .
spruce grouse ( falcipennis canadensis ) and ruffed grouse ( bonasa umbellus ) are typically associated with boreal and mixed deciduous forests , respectively , whereas rock ptarmigan ( lagopus muta ) , willow ptarmigan ( lagopus lagopus ) , and white - tailed ptarmigan ( lagopus leucura ) inhabit more sparsely vegetated arctic and subarctic tundra and subalpine environments ( aldrich , 1963 ) .
sooty grouse ( dendrogapus fuligninosus ) populations in alaska are found in the temperate coastal rainforests of southeastern alaska and sharp - tailed grouse ( tympanuchus phasianellus ) tend to prefer open grassy habitat with shrub - like brush in interior alaska ( aldrich , 1963 , dublin and taras , 2005 ) .
grouse and ptarmigan populations are secondary prey species to many predators and undergo cyclical population fluctuations ( holmstad et  al . ,
some researchers have hypothesized that parasitic infection may contribute to these population fluctuations , either as a direct cause of mortality or via a secondary effect such as decreasing mobility or delayed flushing response ( fallis , 1945 , holmstad et  al . , 2006 , skirnisson et  al . , 2012 ) .
early microscopic studies from north america determined that most tetraonid species were infected by at least one genus of haemosporidian parasite , with many individuals showing diverse infections of multiple parasite genera at moderate to high prevalence ( fallis , 1945 , stabler et  al .
, 1967a , stabler et  al . , 1967b , bennett and inder , 1972 , mahrt , 1981 , forbes et  al . , 1994 ) .
although research on haemosporidian infections in alaskan tetraonids has been extremely limited , two historical studies focused on grouse and ptarmigan species within specific regions of the state .
( 1967a ) sampled rock ptarmigan in interior alaska using microscopy and found 88% prevalence of leucocytozoon parasites and no evidence of haemoproteus or plasmodium infection .
( 1967b ) collected blood smears from spruce grouse in the southwestern and southcentral regions of alaska , where they detected 80% and 60% prevalence of leucocytozoon and haemoproteus , respectively , and no plasmodium .
molecular - based detection methods have become more widely used in the past decade , allowing for detailed description of parasite genetic diversity , including identification of unique parasite lineages in each genus ( bensch et  al . , 2000 , hellgren et  al .
since the advent of these methods , however , there has been little research on the genetic diversity of haemosporidians in grouse and ptarmigan worldwide ( sato et  al . , 2007 ) , and
none in alaska . to address the lack of current knowledge about blood parasite infection in grouse and ptarmigan species in alaska and to provide information on the genetic diversity of parasite
lineages that infect this group of resident birds , our objectives for this study were to : ( 1 ) to estimate the prevalence of haemosporidian infections ( genera haemoproteus , plasmodium , and leucocytozoon ) in grouse and ptarmigan species throughout the state of alaska ; ( 2 ) obtain genetic information to estimate the diversity of identified haemoproteus , plasmodium , and leucoctyzoon parasite lineages ; and ( 3 ) evaluate the relationship between prevalence of blood parasite infection and geographic region , age , and species of the host .
our results will provide some of the first genetic data on haemosporidian lineages infecting tetraonid species and add new information about blood parasite prevalence and distribution in resident avian hosts from northern regions .
wings from hunter - harvested grouse and ptarmigan were voluntarily submitted to the alaska department of fish and game through their grouse and ptarmigan wing collection program from october 2012 to november 2014 .
we received 459 samples from four species of grouse ( ruffed grouse , spruce grouse , sharp - tailed grouse , and sooty grouse ) and three species of ptarmigan ( willow ptarmigan , rock ptarmigan , and white - tailed ptarmigan ) from multiple geographic regions of alaska ( fig .  1 ) .
all collections occurred between august and april , primarily during the non - breeding season .
wings were stored individually at  20 c until tissue extraction and subsequent genetic analysis .
we extracted dna from approximately 500  mg of muscle tissue of frozen wings using a dneasy blood and tissue extraction kit ( qiagen , valencia , ca ) following the manufacturer 's protocol . in order to confirm the viability of each dna extraction , a 695 base pair ( bp ) fragment of the avian mitochondrial dna ( mtdna ) cytochrome oxidase
i ( coi ) gene was amplified following protocols described by kerr et  al .
each sample was considered viable if bands were visible when pcr product was visualized on 0.8% agarose gels stained with gel red nucleic acid gel stain ( biotium , hayward , ca ) . each extracted dna sample that proved viable via our coi positive control (
n  =  459 ) was subsequently screened for the presence of haemoproteus , plasmodium , and leucocytozoon parasites using a nested - pcr protocol described by hellgren et  al .
( 2004 ) , which allows for simultaneous identification of infections from all three parasite genera .
a 479 bp fragment of haemosporidian mtdna cytochrome b gene was bi - directionally sequenced for all positive samples using identical primers from pcr .
positive pcr products were purified with exo - sapit ( usb inc . , cleveland , oh ) according to the manufacturer 's protocol , and sequencing was conducted using big dye terminator v3.1 mix ( applied biosystems , foster city , ca ) and analyzed on an abi 3730xl automated dna sequencer ( applied biosystems , foster city , ca ) .
sequence data were cleaned up and edited using sequencher 5.0.1 software ( gene codes corp . , ann arbor , mi ) .
raw sequences were assigned to one of three parasite genera ( haemoproteus , plasmodium , or leucocytozoon ) using the nucleotide blast function available on the national center for biotechnology information ( ncbi ) website .
samples that did not produce double - stranded dna sequence during at least one of the pcr reactions or failed to be assigned by blast were considered negative . to determine the relative frequency and genetic distances for all unique haemosporidian haplotypes in grouse and ptarmigan samples
, we created a median - joining minimum spanning network using network 4.6.1 ( fig .  2 ; bandelt et  al . , 1999 ) .
additionally , we compared haemosporidian haplotypes from our samples to previously identified lineages available ( as of 27 november , 2015 ) on the genbank and malavi databases ( bensch et  al . , 2009 ) using their basic local alignment search tool function ( blast ) to determine if they matched lineages reported in other avian hosts . in order to provide insight into patterns of host or geographic conservatism of these parasites ,
phylogenetic analyses were performed on sequence data obtained from our samples and reference sequences collected from genbank and malavi public databases ( bensch et  al . , 2009 ) .
we included a single representative of each unique haemosporidian mtdna haplotype identified in our study ( n  =  23 ; genbank accession numbers ku257614-ku257636 ) , as well as representative haplotypes from previous studies conducted on haemosporidians in alaskan avian species ( ramey et  al .
, 2012 , oakgrove et  al . , 2014 , smith et  al . , 2015 ) , lineages isolated from black fly vectors in colorado ( murdock et  al . , 2015 ) , and representative lineages that have been previously identified in birds from the phasianidae family obtained from the malavi database , to examine geographic or host - specific patterns .
all sequences were aligned and cropped to a final length of 413 bp to match the length of the shortest lineage in our phylogeny ; any lineages shorter than this or those that contained ambiguous bases from co - infection by multiple haemosporidian lineages were excluded from further analysis .
, 2012 ) using a general time reversible model and a gamma distribution for among - site variation ( gtr  +  g ) .
four heated chains were used and the first 25,000 sampled trees were discarded as burn - in . our analysis was replicated multiple times to ensure consistent results with each analysis running for a minimum of 3.0    10 generations or until the split frequencies of the posterior probability standard deviation stabilized below 0.01 .
trees were sampled every 1000 generations and rooted with mammalian plasmodium ( plasmodium yoelii : genbank accession number ay099051 ) , based on a phylogeny described in perkins and schall ( 2002 ) . to investigate geographic patterns of haemosporidian prevalence , we identified six geographic regions for analysis by grouping together multiple game management units and sub - units based on proximity and separation by major geographical features ( fig . 
1 ) : southcoastal ( kenai peninsula and southeastern alaska ; gmus 1c , 1d , 2 , 7 , 15a , 15b , and 15c ) , southcentral ( anchorage area and the matanuska - susitna valley ; gmus 13a , 13d , 14a , 14c , 16a , and 16b ) , southwestern ( bristol bay , alaska peninsula , and eastern aleutian islands ; 9d , 9e , and 17c ) , southern interior ( south side of alaska range ; gmus 12 , 13b , and 13e ) , northern interior ( north side of alaska range ; gmus 20a-20e and 25c ) , and seward peninsula ( gmu 22c ) .
we used generalized linear models ( proc genmod with binomial distribution of errors and logit link function ; sas institute , version 9.4 ) to investigate the relationship between parasite infection status ( dependent , binary variable ) and the following explanatory variables associated with capture of each bird : region , age ( adult [ > 1 year of age ] versus juvenile [ < 1 year of age ] ) , and species group ( grouse versus ptarmigan ) .
we also included interaction terms between region and species group and between age and species group as we hypothesized that prevalence in each species group could vary geographically and by age .
however , we did not expect age - related patterns of infection to vary by region so this term was not included .
we also did not include sex in our models because sex was not known for many of the samples .
for plasmodium , because no infections were detected in ptarmigan , species group was not included and only region and age were considered in our models .
we evaluated the relative support for models in each candidate set using akaike 's information criterion corrected for sample size ( aicc ) ( burnham and anderson , 2002 ) . within each candidate model
set , we then calculated the akaike weight ( w ) for each model , which is the relative likelihood of the model , given the data and the models being considered ( burnham and anderson , 2002 ) . in each analysis
, we used an all - subsets modeling approach , in which the full model included all additive combinations of the independent variables and interactions ; all possible subsets of that model were included in the candidate model set .
we present least - squares means for leucocytozoon and haemoproteus prevalence by species group and region , calculated from best - supported generalized linear models . for plasmodium ,
we report least - squares means for age and region from single - factor models including grouse only .
of the 459 grouse and ptarmigan wings screened , 342 individual samples were identified as positive for at least one genus of haemosporidian parasite , with an overall apparent prevalence of 74.5% ( 2.1% ) . from these 342 infected individuals , we detected a total of 374 individual haemosporidian infections consisting of leucocytozoon ( n  =  241 ) , haemoproteus ( n  =  94 ) , and plasmodium ( n  =  39 ; table  1 ) .
the majority of individuals had single infections ( 90.6%    1.6% ; n  =  310 ) , with a smaller number having co - infections of leucocytozoon and haemoproteus ( 5.3%    1.2% ; n  =  18 ) or leucocytozoon and plasmodium ( 4.1%    1.1% ; n  =  14 ) .
apparent prevalence of leucocytozoon ( 52.5%    2.3% ) was the highest among the three genera , followed by haemoproteus ( 20.5%    1.9% ) , and plasmodium ( 8.5%    1.3% ) .
leucocytozoon infections were represented across all seven species of bird hosts sampled , ranging from 21% to 83% prevalence depending upon species ( table  1 ) .
plasmodium parasites were detected in only three of the grouse species sampled ( spruce , ruffed , and sharp - tailed grouse ) , and no ptarmigan species ( table  1 ) .
our duplicate screening of each sample confirmed 423 consistent screening results for leucocytozoon and 431 for haemoproteus / plasmodium , the latter of which were not differentiated since they are identified with a single molecular primer set .
this concordance of 92.2% and 93.9% for leucocytozoon and haemoproteus / plasmodium , respectively , is comparable with detection probabilities reported by other studies using the same molecular protocol ( ramey et  al . , 2015 ) .
among ptarmigan , apparent prevalence of leucocytozoon infection was relatively high in all 6 regions sampled , ranging from 47.1% ( 12.1% ; n  =  17 ) in southwestern alaska to 96.3% ( 3.6% ; n  =  27 ) on the seward peninsula ( table  2 ) .
prevalence of leucocytozoon was generally lower in grouse , ranging from 2.0% ( 1.9% ; n  =  51 ) in southwestern alaska to 47.1% ( 12.1% ; n  =  17 ) in the southern interior ( table  2 ) .
similar geographic patterns were observed for leucocytozoon across the two species groups , with the exception of the seward peninsula , for which we had no grouse samples .
apparent prevalence of haemoproteus was highest among grouse in southwestern alaska ( 98.0%    1.9% ) and lowest in the southern interior ( 5.9%    5.7% ; table  2 ) .
the highest prevalence of haemoproteus infection among ptarmigan also occurred in southwestern alaska but was only 11.8% ( 7.8% ; table  2 ) . among grouse , plasmodium prevalence was highest in southcentral alaska ( 39.2%    6.8% ) and was higher among adults ( 26.6%    2.8% ) than juveniles ( 13.4%    2.4% ; table  2 ) .
analysis of haemosporidian mtdna sequences revealed a total of 23 unique haplotypes ( fig .  2 ; table  3 ) .
haplotypes were identified as haemoproteus ( n  =  6 ) , plasmodium ( n  =  3 ) , and leucocytozoon ( n  =  14 ) .
haplotype frequency varied from detection in only a single individual to detection in 75 separate individuals of multiple species ( table  3 ) .
a total of 15 haplotypes were isolated to a single host species , some of which included multiple individuals sampled across large geographic distances .
the remaining eight haplotypes were detected in multiple host species , with four of these detected in both grouse and ptarmigan hosts and the remaining four restricted to one species group ( table  3 ) .
comparison to haemosporidian lineages on malavi and genbank public databases showed that three of our haplotypes were identical to lineages previously identified in other avian hosts ( table  3 ) .
akgph03 ,  which was isolated from four different tetraonid species ( table  3 ) , was identical to lineage cygnus01 , and leucocytozoon haplotype 
akgpp01  detected in our study , which was isolated from three of the four species of grouse sampled , was identical to the bt7 lineage found on the malavi database .
the 92 separate haemosporidian lineages in our bayesian tree clustered by genus as expected , with high posterior probability support for all of the larger clades ( fig .  3 ) .
leucocytozoon haplotypes detected in our study were more diverse than those of haemoproteus or plasmodium .
analysis with reference sequences grouped most of our leucocytozoon haplotypes by host family , with one clade of several haplotypes deeply diverged from the others and grouped with sequences obtained from black flies and one lineage previously identified in phasianid species .
haemoproteus and plasmodium haplotypes from our study showed less host conservatism than those of leucocytozoon .
our phylogeny provided support for haemoproteus clades comprised of our sequences and those detected in waterfowl from alaska .
additionally , all of our plasmodium haplotypes were grouped in a single clade with lineages isolated from tundra swans and mallards in alaska ( fig . 
the model receiving the most support ( w  =  0.32 ) for explaining variation in leucocytozoon prevalence included region , species group , and a region*species group interaction term ( table  4 ) , although models that also included age ( aicc  =  0.22 ; w  =  0.28 ) and an age*species group interaction ( aicc  =  0.13 ; w  =  0.30 ) had almost equal support ( table  4 ) . for haemoproteus ,
the best - supported model ( w  =  0.62 ) included region , species group , and a region*species group interaction .
there was moderate support ( w  =  0.23 ) for a model that also included age ( table  4 ) .
plasmodium was not detected in any ptarmigan samples ( making prevalence inestimable ) so we modeled the likelihood of plasmodium infection in grouse only .
the best - supported model for explaining plasmodium prevalence in grouse included region and age ( w  =  0.99 ) .
our study constitutes the most comprehensive survey of haemosporidian parasites in alaskan tetraoninds to date , and has revealed a relatively high overall prevalence of haemoproteus , plasmodium , and leucocytozoon ( 74.5%    2.1% ) .
prevalence of all three genera varied by region , with differing patterns of infection between grouse and ptarmigan .
surprisingly , we detected no plasmodium infections in ptarmigan , even across the relatively broad geographic area sampled .
apparent prevalence of leucocytozoon was high across most regions and host species , particularly among ptarmigan ( table  2 ) .
prevalence ranged from 47% to 96% in ptarmigan and , with the exception of southwestern alaska , exceeded 30% in grouse .
these results suggest that tetraonid populations across alaska are exposed to relatively high concentrations of suitable vectors ( table  2 ) .
consistently higher leucocytozoon prevalence in ptarmigan compared to grouse in each geographic region is interesting given the respective habitat selection of these birds .
ptarmigan , particularly rock and white - tailed ptarmigan , typically inhabit alpine environments in the regions from which we received samples , and grouse tend to occur in more lowland regions closer to streams and boreal forests ( aldrich , 1963 , dublin and taras , 2005 ) .
leucocytozoon parasites are transmitted to vertebrate hosts exclusively by hematophagus black flies ( simuliidae ) , except for leucocytozoon caulleryi , which is transmitted by biting midges ( culicoides ) and is specific to domestic chickens ( gallus gallus ) ( valkinas , 2005 ) .
simuliids are typically found in the highest density near moving water , where they undergo their larval stages attached to rocks ( borror and delong , 1954 ) .
leucocytozoon lovati , a parasite specific to tetraoninds and the only species previously described in ptarmigan , is spread by simuliid species that are known to undergo larval stages in smaller , slow moving streams or the outflows of impounded water bodies , such as beaver ponds ( currie , 1997 ) .
further investigation into parasite - vector relationships for leucocytozoon in alaska is needed to determine the underlying factors responsible for differences in prevalence between grouse and ptarmigan .
in contrast to the pattern observed with leucocytozoon , prevalence of haemoproteus was generally higher in grouse than in ptarmigan .
overall , infection with haemoproteus parasites was less common than with leucocytozoon in most species and regions .
one notable exception was the remarkably high prevalence of haemoproteus ( 98.0% ) in grouse from the southwestern region , where leucocytozoon infection in this species group was relatively rare ( 2.0% ; table  2 ) .
1967a ) also detected substantially higher prevalence of haemoproteus parasites in spruce grouse sampled from lake aleknagik , which is near where the majority of our southwestern samples originated , than from spruce grouse collected on the kenai peninsula of alaska .
additionally , tundra swans ( cygnus columbianus ) sampled in the bristol bay lowlands had a higher prevalence of haemoproteus than swans sampled in other regions of alaska ( ramey et  al . , 2012 ) .
haemoproteus parasites are primarily transmitted by biting midges of the genus culicoides , although a small number of these parasites use hippoboscid flies ( hippoboscidae ) as vectors ( valkinas , 2005 ) .
culicoides require semi - aquatic or very moist habitats in order to reproduce , and often form large swarms near the area where they are hatched , but can also travel great distances through both wing- and wind - propelled dispersal ( borkent , 2004 ) .
the similarity of our results to previous studies suggests increased opportunities for transmission of haemoproteus in southwestern alaska , particularly around bristol bay .
it is possible that local habitat and environmental factors in this region support larger quantities of biting midges , leading to increased transmission of haemoproteus infections .
additionally , habitat associations of grouse with lowland areas and boreal forest habitats may help to explain the notably higher prevalence of haemoproteus in grouse than ptarmigan across all regions .
plasmodium infections were restricted to grouse species , with not a single infection detected in any of the three species of ptarmigan .
the majority of these infections occurred in spruce grouse from the southcentral region , which is largely dominated by boreal forest habitat .
this habitat is suitable for species of culicidae , the family of mosquitoes that transmit plasmodium , and may support relatively high densities of potential vectors in this area .
plasmodium is also the most temperature - sensitive of the three genera of haemosporidian parasites we screened for in this study ( valkinas , 2005 ) , and the lowland areas of the southcentral region have summertime temperatures that are conducive to completion of the parasite 's life cycle ( wilkinson et  al . ,
we detected plasmodium in a smaller proportion of grouse samples from other regions , including the northern interior , which is similar in latitude to the most northerly detections of locally - acquired plasmodium infections reported in passerines ( oakgrove et  al .
age was also an important factor in determining prevalence of plasmodium in grouse , with adults infected at a higher frequency than juveniles ( table  2 ) .
, 2007 , cosgrove et  al . , 2008 , wilkinson et  al . , 2016 ) and may be explained in part by the longer temporal period to which adults are exposed to potential vectors ( valkinas , 2005 ) .
the apparent absence of plasmodium in ptarmigan is notable , especially given the broad geographic area sampled , which included regions where plasmodium has been detected in a diversity of other avian hosts ( loiseau et  al . , 2012 , reeves et  al . , 2015 , wilkinson et  al . , 2016 ) .
willow ptarmigan are commonly found in areas of dense vegetation in lowland areas and valley bottoms ( hannon et  al . , 1998 ) , where temperature and wind conditions are conducive to high vector density .
additionally , although the alpine environments typical of white - tailed and rock ptarmigan may have a lower abundance of mosquitoes than forest or shrub habitats , mosquitoes do occur even at higher elevations ( m. smith , personal observation ) .
the complete absence of plasmodium in our alaskan ptarmigan samples suggests that factors such as host behavior , local habitat selection , or other host life history traits are important determinants of infection status .
interestingly , to our knowledge , no other studies of blood parasites in ptarmigan have detected plasmodium infections , including in japan ( hagihara et  al . , 2004 ) ,
norway ( holmstad and skorping , 1998 ) , iceland ( skirnisson et  al . , 2012 ) , and canada ( mahrt , 1981 ) .
it is unclear whether the higher contemporary prevalence of plasmodium we observed among spruce grouse in southcentral alaska represents a true increase or if differences in detection methods ( pcr in our study versus microscopy in stabler et  al .
other studies , dating back more than 70 years , have detected plasmodium in grouse from more southerly locations in canada , including british columbia and ontario , ( e.g. , fallis , 1945 , forbes et  al . , 1994 ) .
the potential for an increase in plasmodium prevalence associated with climate warming has been suggested in other studies of blood parasites in alaska and warrants additional research ( loiseau et  al . , 2012 , wilkinson et  al . , 2016 ) .
molecular analysis of historical blood samples from tetraonids across a larger geographic area , if available , would provide greater inference about potential long - term changes in haemosporidian prevalence .
it should also be noted that our collection of samples during the non - breeding season may reflect a lower apparent prevalence than samples collected during summer , when vectors are active , especially for plasmodium , for which parasitemia drops significantly during the chronic infection phase , making detection difficult ( valkinas , 2005 ) .
similar to other studies of haemosporidians in alaska , co - infections with leucocytozoon and haemoproteus were the most common among grouse and ptarmigan , with fewer co - infections of leucocytozoon and plasmodium . however , we detected a lower prevalence of co - infection compared to other alaskan species ( oakgrove et  al .
, 2015 , reeves et  al . , 2015 ) . while we did not detect any haemoproteus / plasmodium co - infections
, this finding is n't surprising as it is rare to detect these type of co - infections using molecular methods .
haemoproteus infections typically produce much higher parasitemia in avian hosts than plasmodium ( valkinas , 2005 ) , making it difficult to detect plasmodium in a co - infected individual . future work involving microscopic examination of blood smears should be conducted to determine dynamics of co - infections between haemoproteus and plasmodium parasites .
comparison of haemosporidian haplotypes identified in our study revealed differing levels of host specificity for haemoproteus , plasmodium , and leucocytozoon parasites .
leucocytozoon lineages in our phylogeny were grouped largely by host family . the majority of our haplotypes clustered into a single , monophyletic group with other lineages isolated from phasianid hosts and one representative lineage from black flies ( fig .  3 ) .
the remaining leucocytozoon haplotypes clustered into another , highly diverged clade containing haplotype akgpl06 , which was identical to previously identified lineage lamut01 .
lamut01 , one of the only leucocytozoon lineages isolated from ptarmigan species , was detected in japanese rock ptarmigan ( lagopus mutua japonicas ) and identified as leucocytozoon lovati ( sato et  al . , 2007 ) .
although there are four species of leucocytozoon known to infect phasianid hosts , the only species that has been previously identified in ptarmigan is l.  lovati .
previous studies have observed this species in ptarmigan populations around the world , typically at relatively high prevalence ( stabler et  al .
, 1967a , mahrt , 1981 , murata et  al . , 2007 ) .
2 ) , combined with the level of divergence between the two clades in our phylogeny ( fig . 
cryptic speciation of leucocytozoon species has been suggested by several other studies ( sehgal et  al . , 2006 , reeves et  al . , 2015 ) , although the lack of genetic analyses of leucocytozoon species in tetraonids may explain why only one species has previously been reported in ptarmigan .
microscopic examination of blood smears from ptarmigan would be beneficial to provide identification of parasite morphospecies .
haemoproteus haplotypes clustered into two , tightly grouped clades consisting of either tetraonid or alaskan waterfowl lineages ( fig .  3 ) .
lineage akgph03 , detected in four tetraonid species ( table  3 ) , was identical to cygnus01 .
cygnus01 has been detected in anseriform hosts , including tundra swans ( cygnus columbianus ) in minnesota and alaska ( ricklefs and fallon , 2002 ) , northern pintails ( anas acuta ) in california ( ramey et  al . ,
2012 ) , and cinnamon teal ( anas cyanoptera ) in south america ( smith and ramey , 2015 ) .
this lineage had previously only been documented in waterfowl , and is associated with haemoproteus nettionis , which is thought to be specific to anseriform hosts .
smith and ramey ( 2015 ) showed that haemoproteus parasites in waterfowl were relatively host specific , even intercontinentally .
however , the fact that the haemoproteus lineages we detected in tetraonids clustered so tightly with those from waterfowl and the identification of identical lineages in the two families suggests more complex host - parasite relationships than previously thought .
although the identification of identical parasite lineages in both waterfowl and tetraonids is not necessarily indicative of infection by the same parasite species , this novel finding provides incentive for further study of haemoproteus parasites in alaska .
the addition of microscopic methods would allow for morphological identification of these parasites to definitively determine if diverse hosts are being parasitized by the same haemoproteus species .
our phylogeny grouped all plasmodium haplotypes from our study , as well as plasmodium lineages from alaskan waterfowl , into a single , monophyletic clade .
akgpp01  had a 100% identity match to lineage bt7 ( also described as p43 in genbank ) .
this plasmodium lineage displays generalist behavior and has been documented in 37 bird species to date , spanning 11 families and four orders ( bensch et  al . , 2009 ) . in alaska ,
bt7 has been reported in waterfowl ( reeves et  al . , 2015 ) , black - capped chickadees ( poecile atricapillus ; 
 wilkinson et  al . , 2016
plasmodium species have historically been known to display generalist behavior in avian hosts ( bensch et  al . , 2000 , waldenstrm et  al .
the fact that plasmodium haplotypes from our study were more genetically diverged from other lineages isolated from phasianidae hosts than those from diverse alaskan bird species supports this trend of generalist behavior in plasmodium parasites .
results from our study of haemosporidian parasites in alaskan tetraonids revealed variation in prevalence of haemoproteus , plasmodium , and leucocytozoon by species and geographic region across alaska .
the absence of plasmodium infection in ptarmigan species and the markedly higher prevalence of leucocytozoon in ptarmigan and haemoproteus in grouse were notable findings .
future work involving finer - scale sampling would allow for better resolution of regional and local environmental factors that may influence parasite prevalence .
additionally , study of parasite - vector relationships for haemosporidians in alaska is needed to understand the mechanisms responsible for apparent differences in haemosporidian infections in grouse and ptarmigan throughout the state .
phylogenetic assessment of haemosporidian haplotypes showed a diverse array of parasite lineages infecting tetraonid species , with apparent generalist behavior by haemoproteus and plasmodium parasites , and host conservatism in leucocytozoon lineages .
several haplotypes in our study were identical to haemosporidian lineages previously detected in other host families , including waterfowl and passerines .
the most essential step for future work involving diversity of haemosporidians in alaskan tetraonids would be the addition of microscopic examination of blood smears to molecular and phylogenetic analysis .
doing so would allow researchers to match parasite morphospecies to genetic lineages and determine the level of host sharing or specificity among these blood parasites .
results from this study provide important baseline information for future research examining potential effects of climate warming on haemosporidian infections in avian species from alaska and other northern regions . given the status of grouse and ptarmigan as year - round residents , long - term monitoring of haemoproteus , plasmodium , and
leucocytozoon infections in alaskan tetraonids will allow for assessment of local changes in parasite prevalence and distribution .
additionally , the genetic data provided here offer some of the first information on haemosporidian lineages in tetraonid species and are useful for phylogenetic comparison in grouse and ptarmigan worldwide .
 projections related to future climate warming indicate the potential for an increase in the distribution and prevalence of blood parasites in northern regions . 
 however , baseline data are lacking for resident avian host species in alaska . 
 grouse and ptarmigan occupy a diverse range of habitat types throughout the northern hemisphere and are among the most well - known and important native game birds in north america . 
 information regarding the prevalence and diversity of haemosporidian parasites in tetraonid species is limited , with few recent studies and an almost complete lack of genetic data . to better understand the genetic diversity of haemosporidian parasites in alaskan tetraonids and to determine current patterns of geographic range and host specificity , we used molecular methods to screen 459 tissue samples collected from grouse and ptarmigan species across multiple regions of alaska for infection by leucocytozoon , haemoproteus , and plasmodium blood parasites . 
 infections were detected in 342 individuals , with overall apparent prevalence of 53% for leucocytozoon , 21% for haemoproteus , and 9% for plasmodium . 
 parasite prevalence varied by region , with different patterns observed between species groups ( grouse versus ptarmigan ) . 
 leucocytozoon was more common in ptarmigan , whereas haemoproteus was more common in grouse . 
 we detected plasmodium infections in grouse only . 
 analysis of haemosporidian mitochondrial dna cytochrome b sequences revealed 23 unique parasite haplotypes , several of which were identical to lineages previously detected in other avian hosts . 
 phylogenetic analysis showed close relationships between haplotypes from our study and those identified in alaskan waterfowl for haemoproteus and plasmodium parasites . 
 in contrast , leucocytozoon lineages were structured strongly by host family . 
 our results provide some of the first genetic data for haemosporidians in grouse and ptarmigan species , and provide an initial baseline on the prevalence and diversity of blood parasites in a group of northern host species . 
Introduction
Materials and methods
Results
Discussion