Determining evolutionary origin and phylogenetic relationships of mallard-like ducks of Oceania, greater Indonesia, and the Philippines with ddRAD-seq data.

AIM: We aim to determine the evolutionary origins and population genetics of mallard-like ducks of Oceania, greater Indonesia, and the Philippines. LOCATION: Oceania, greater Indonesia, and the Philippines. TAXON: Mallard (Anas platyrhynchos), Pacific black duck (A. superciliosa spp.), and Philippine duck (A. luzonica) METHODS: Thousands of nuclear ddRAD-seq loci and the mitochondrial DNA control region were assayed across individuals representative of each species' range. We assessed population structure and phylogenetic relationships, as well as estimated demographic histories to reconstruct the biogeographical history of each species. RESULTS: Philippine and Pacific black ducks represent unique genetic lineages that diverged from the mallard 1-2 million years ago. We find no support for the Philippine duck representing a hybrid species as once posited; however, their low levels of genetic diversity requires further attention. We find a lack of substructure among Philippine ducks. However, we found pronounced differentiation between subspecies of Pacific black ducks, especially between A. s. superciliosa from New Zealand and A. s. rogersi from Australia, Papua New Guinea, and Timor-Leste, Indonesia. Anas superciliosa pelewensis gave mixed results; individuals from the Solomon Islands were differentiated from the other subspecies, but those from the island of Aunu'u, American Samoa, were genetically more similar to A. s. rogersi than A. s. pelewensis samples from the Solomon Islands. Finally, we find limited evidence of interspecific gene flow at evolutionary scales, and mallard introgression among contemporary samples. MAIN CONCLUSIONS: Mallard-like ducks radiated across Oceania, greater Indonesia, and the Philippines within the last 2 million years. Only the Pacific black duck showed unique sub-structuring that largely followed known sub-species ranges, except for A. s. pelewensis. We posit that the high interrelatedness among Solomon Island samples suggests that their genetic distinctiveness may simply be the result of high levels of genetic drift. In contrast, we conclude that mainland Australian Pacific black ducks were the most likely source for the recent colonization of American Samoa. As a result, our findings suggest that either the A. s. pelewensis subspecies designations and/or its geographical range may require re-evaluation. Continued re-evaluation of evolutionary and taxonomic relationships is necessary when attempting to reconstruct and understand biogeographical histories, with important implications towards any attempts to implement conservation strategies.

Sustained plumage divergence despite weak genomic differentiation and broad sympatry in sister species of Australian woodswallows (Artamus spp.).

Plumage divergence can function as a strong premating barrier when species come into secondary contact. When it fails to do so, the results are often genome homogenization and phenotypic hybrids at the zone of contact. This is not the case in the largely sympatric masked woodswallow and white-browed woodswallow species (Passeriformes: Artamidae: Artamus spp) complex in Australia where phenotypic integrity is sustained despite no discernible mitochondrial structure in earlier work. This lack of structure may suggest recent divergence, ongoing gene flow or both, and phenotypic hybrids are reported albeit rarely. Here, we further assessed the population structure and differentiation across the species' nuclear genomes using ddRAD-seq. As found in the mitochondrial genome, no structure or divergence within or between the two species was detected in the nuclear genome. This coarse sampling of the genome nonetheless revealed peaks of differentiation around the genes SOX5 and AXIN1. Both are involved in the Wnt/ß-catenin signalling pathway, which regulates feather development. Reconstruction of demographic history and estimation of parameters supports a scenario of secondary contact. Our study informs how divergent plumage morphs may arise and be sustained despite whole-genome homogenization and reveals new candidate genes potentially involved in plumage divergence.

Adaptive introgression of the beta-globin cluster in two Andean waterfowl.

Introgression of beneficial alleles has emerged as an important avenue for genetic adaptation in both plant and animal populations. In vertebrates, adaptation to hypoxic high-altitude environments involves the coordination of multiple molecular and cellular mechanisms, including selection on the hypoxia-inducible factor (HIF) pathway and the blood-O2 transport protein hemoglobin (Hb). In two Andean duck species, a striking DNA sequence similarity reflecting identity by descent is present across the ~20 kb ß-globin cluster including both embryonic (HBE) and adult (HBB) paralogs, though it was yet untested whether this is due to independent parallel evolution or adaptive introgression. In this study, we find that identical amino acid substitutions in the ß-globin cluster that increase Hb-O2 affinity have likely resulted from historical interbreeding between high-altitude populations of two different distantly-related species. We examined the direction of introgression and discovered that the species with a deeper mtDNA divergence that colonized high altitude earlier in history (Anas flavirostris) transferred adaptive genetic variation to the species with a shallower divergence (A. georgica) that likely colonized high altitude more recently possibly following a range shift into a novel environment. As a consequence, the species that received these ß-globin variants through hybridization might have adapted to hypoxic conditions in the high-altitude environment more quickly through acquiring beneficial alleles from the standing, hybrid-origin variation, leading to faster evolution.

ddRAD-seq data reveal significant genome-wide population structure and divergent genomic regions that distinguish the mallard and close relatives in North America.

Recently evolved species typically share genetic variation across their genomes due to incomplete lineage sorting and/or ongoing gene flow. Given only subtle allele frequency differences at most loci and the expectation that divergent selection may affect only a tiny fraction of the genome, distinguishing closely related species based on multi-locus data requires substantial genomic coverage. In this study, we used ddRAD-seq to sample the genomes of five recently diverged, New World "mallards" (Anas spp.), a group of dabbling duck species characterized by diagnosable phenotypic differences but minimal genetic differentiation. With increased genomic sampling, we aimed to characterize population structure within this group and identify genomic regions that may have experienced divergent selection during speciation. We analyzed 3,017 autosomal ddRAD-seq loci and 177 loci from the Z-chromosome. In contrast to previous studies, the ddRAD-seq data were sufficient to assign individuals to their respective species or subspecies and to generate estimates of gene flow in a phylogenetic framework. We find limited evidence of contemporary gene flow between the dichromatic mallard and several monochromatic taxa, but find evidence for historical gene flow between some monochromatic species pairs. We conclude that the overall genetic similarity of these taxa likely reflects retained ancestral polymorphism rather than recent and extensive gene flow. Thus, despite recurring cases of hybridization in this group, our results challenge the current dogma predicting the genetic extinction of the New World monochromatic dabbling ducks via introgressive hybridization with mallards. Moreover, ddRAD-seq data were sufficient to identify previously unknown outlier regions across the Z-chromosome and several autosomal chromosomes that may have been involved in the diversification of species in this recent radiation.

Persistence of an endangered native duck, feral mallards, and multiple hybrid swarms across the main Hawaiian Islands.

Interspecific hybridization is recognized as an important process in the evolutionary dynamics of both speciation and the reversal of speciation. However, our understanding of the spatial and temporal patterns of hybridization that erode versus promote species boundaries is incomplete. The endangered, endemic koloa maoli (or Hawaiian duck, Anas wyvilliana) is thought to be threatened with genetic extinction through ongoing hybridization with an introduced congener, the feral mallard (A. platyrhynchos). We investigated spatial and temporal variation in hybrid prevalence in populations throughout the main Hawaiian Islands, using genomic data to characterize population structure of koloa, quantify the extent of hybridization, and compare hybrid proportions over time. To accomplish this, we genotyped 3,308 double-digest restriction-site-associated DNA (ddRAD) loci in 425 putative koloa, mallards, and hybrids from populations across the main Hawaiian Islands. We found that despite a population decline in the last century, koloa genetic diversity is high. There were few hybrids on the island of Kaua'i, home to the largest population of koloa. By contrast, we report that sampled populations outside of Kaua'i can now be characterized as hybrid swarms, in that all individuals sampled were of mixed koloa × mallard ancestry. Further, there is some evidence that these swarms are stable over time. These findings demonstrate spatial variation in the extent and consequences of interspecific hybridization, and highlight how islands or island-like systems with small population sizes may be especially prone to genetic extinction when met with a congener that is not reproductively isolated.

Host community heterogeneity and the expression of host specificity in avian haemosporidia in the Western Cape, South Africa.

Similar patterns of parasite prevalence in animal communities may be driven by a range of different mechanisms. The influences of host heterogeneity and host-parasite interactions in host community assemblages are poorly understood. We sampled birds at 27 wetlands in South Africa to compare four hypotheses explaining how host community heterogeneity influences host specificity in avian haemosporidia communities: the host-neutral hypothesis, the super-spreader hypothesis, the host specialist hypothesis and the heterogeneity hypothesis. A total of 289 birds (29%) were infected with Plasmodium, Haemoproteus and/or Leucocytozoon lineages. Leucocytozoon was the most diverse and generalist parasite genus, and Plasmodium the most conservative. The host-neutral and host specialist hypotheses received the most support in explaining prevalence by lineage (Leucocytozoon) and genus (Plasmodium and Haemoproteus), respectively. We observed that haemosporidian prevalence was potentially amplified or reduced with variation in host and/or parasitic taxonomic levels of analysis. Our results show that Leucocytozoon host abundance and diversity was influential to parasite prevalence at varying taxonomic levels, particularly within heterogeneous host communities. Furthermore, we note that prevalent mechanisms of infection can potentially act as distinct roots for shaping communities of avian haemosporidia.

Becoming pure: identifying generational classes of admixed individuals within lesser and greater scaup populations.

Estimating the frequency of hybridization is important to understand its evolutionary consequences and its effects on conservation efforts. In this study, we examined the extent of hybridization in two sister species of ducks that hybridize. We used mitochondrial control region sequences and 3589 double-digest restriction-associated DNA sequences (ddRADseq) to identify admixture between wild lesser scaup (Aythya affinis) and greater scaup (A. marila). Among 111 individuals, we found one introgressed mitochondrial DNA haplotype in lesser scaup and four in greater scaup. Likewise, based on the site-frequency spectrum from autosomal DNA, gene flow was asymmetrical, with higher rates from lesser into greater scaup. However, using ddRADseq nuclear DNA, all individuals were assigned to their respective species with >0.95 posterior assignment probability. To examine the power for detecting admixture, we simulated a breeding experiment in which empirical data were used to create F1 hybrids and nine generations (F2-F10) of backcrossing. F1 hybrids and F2, F3 and most F4 backcrosses were clearly distinguishable from pure individuals, but evidence of admixed histories was effectively lost after the fourth generation. Thus, we conclude that low interspecific assignment probabilities (0.011-0.043) for two lesser and nineteen greater scaup were consistent with admixed histories beyond the F3 generation. These results indicate that the propensity of these species to hybridize in the wild is low and largely asymmetric. When applied to species-specific cases, our approach offers powerful utility for examining concerns of hybridization in conservation efforts, especially for determining the generational time until admixed histories are effectively lost through backcrossing.

Speciation genomics and a role for the Z chromosome in the early stages of divergence between Mexican ducks and mallards.

Speciation is a continuous and dynamic process, and studying organisms during the early stages of this process can aid in identifying speciation mechanisms. The mallard (Anas platyrhynchos) and Mexican duck (A. [p.] diazi) are two recently diverged taxa with a history of hybridization and controversial taxonomy. To understand their evolutionary history, we conducted genomic scans to characterize patterns of genetic diversity and divergence across the mitochondrial DNA (mtDNA) control region, 3523 autosomal loci and 172 Z-linked sex chromosome loci. Between the two taxa, Z-linked loci (ΦST  = 0.088) were 5.2 times more differentiated than autosomal DNA (ΦST  = 0.017) but comparable to mtDNA (ΦST  = 0.092). This elevated Z differentiation deviated from neutral expectations inferred from simulated data that incorporated demographic history and differences in effective population sizes between marker types. Furthermore, 3% of Z-linked loci, compared to <0.1% of autosomal loci, were detected as outlier loci under divergent selection with elevated relative (ΦST ) and absolute (dXY ) estimates of divergence. In contrast, the ratio of Z-linked and autosomal differentiation among the seven Mexican duck sampling locations was close to 1:1 (ΦST  = 0.018 for both markers). We conclude that between mallards and Mexican ducks, divergence at autosomal markers is largely neutral, whereas greater divergence on the Z chromosome (or some portions thereof) is likely the product of selection that has been important in speciation. Our results contribute to a growing body of literature indicating elevated divergence on the Z chromosome and its likely importance in avian speciation.

Mito-nuclear discord in six congeneric lineages of Holarctic ducks (genus Anas).

Many species have Holarctic distributions that extend across Europe, Asia and North America. Most genetics research on these species has examined only mitochondrial (mt) DNA, which has revealed wide variance in divergence between Old World (OW) and New World (NW) populations, ranging from shallow, unstructured genealogies to deeply divergent lineages. In this study, we sequenced 20 nuclear introns to test for concordant patterns of OW-NW differentiation between mtDNA and nuclear (nu) DNA for six lineages of Holarctic ducks (genus Anas). Genetic differentiation for both marker types varied widely among these lineages (idiosyncratic population histories), but mtDNA and nuDNA divergence within lineages was not significantly correlated. Moreover, compared with the association between mtDNA and nuDNA divergence observed among different species, OW-NW nuDNA differentiation was generally lower than mtDNA divergence, at least for lineages with deeply divergent mtDNA. Furthermore, coalescent estimates indicated significantly higher rates of gene flow for nuDNA than mtDNA for four of the six lineages. Thus, Holarctic ducks show prominent mito-nuclear discord between OW and NW populations, and we reject differences in sorting rates as the sole cause of the within-species discord. Male-mediated intercontinental gene flow is likely a leading contributor to this discord, although selection could also cause increased mtDNA divergence relative to weak nuDNA differentiation. The population genetics of these ducks contribute to growing evidence that mtDNA can be an unreliable indicator of stage of speciation and that more holistic approaches are needed for species delimitation.

Phylogenetics of a recent radiation in the mallards and allies (Aves: Anas): inferences from a genomic transect and the multispecies coalescent.

Reconstructing species trees by incorporating information from many independent gene trees reduces the confounding influence of stochastic lineage sorting. Such analyses are particularly important for taxa that share polymorphisms due to incomplete lineage sorting or introgressive hybridization. We investigated phylogenetic relationships among 14 closely related taxa from the mallard (Anas spp.) complex using the multispecies coalescent and 20 nuclear loci sampled from a genomic transect. We also examined how treating recombining loci and hybridizing species influences results by partitioning the data using various protocols. In general, topologies were similar among the various species trees, with major clades consistently composed of the same taxa. However, relationships among these clades and among taxa within clades changed among partitioned data sets. Posterior support generally decreased when filtering for recombination, whereas excluding mallards (Anas platyrhynchos) increased posterior support for taxa known to hybridize with them. Furthermore, branch lengths decreased substantially for recombination-filtered data. Finally, concordance between nuclear and morphometric topologies conflicted with those in the mitochondrial tree, particularly with regard to the placement of the Hawaiian duck (A. wyvilliana), Philippine duck (A. luzonica), and two spot-billed ducks (A. zonorhyncha and A. poecilorhyncha). These results demonstrate the importance of maximizing sequence length and taxon sampling when inferring taxonomic relationships that are confounded by extensive allele sharing.

Heteropatric speciation in a duck, Anas crecca.

Heteropatric differentiation is a mode of speciation with gene flow in which divergence occurs between lineages that are in sympatry and allopatry at different times during cyclic spatial movements. Empirical evidence suggests that heteropatric differentiation may prove to be common among seasonally migratory organisms. We examined genetic differentiation between the sedentary Aleutian Islands population of green-winged teal (Anas crecca-nimia) and its close migratory relative, the Eurasian, or Old World (OW), Anas c. crecca population, a portion of which passes through the range of nimia during its seasonal migrations. We also examined its relationship with the parapatric North American, New World (NW), A. c. carolinensis population. Sequence data from eight nuclear introns and the mtDNA control region showed that the nimia-crecca divergence occurred much more recently than the deeper crecca-carolinensis split (~83 000 years vs. ~1.1 Myr). Despite considerable spatial overlap between crecca and nimia during seasonal migration, three key predictions of heteropatric differentiation are supported: significant genetic divergence (overall mean Φst  = 0.07), low gene flow (2Ne m ~ 1.8), and an effective population size in nimia that is not especially low (Ne  ~ 80 000 individuals). Similar levels of gene flow have come into nimia from carolinensis, but no detectable nuclear gene flow has gone out of nimia into either OW (crecca) or NW (carolinensis) populations. We infer that adaptations of these populations to local optima in different places (e.g. each matching their reproductive effort to different resource blooms) promote genetic isolation and divergence despite periods of sympatry between them, as the heteropatric model predicts.

Host associations, biogeography, and phylogenetics of avian malaria in southern African waterfowl.

The relevance of spatial variation in the environment and host communities for parasite community composition is poorly documented, creating a need for additional case studies from which general principles can be developed. Avian malaria in southern African waterfowl has not previously been studied. As a first step towards documenting and understanding its biogeography, we used PCR and molecular sequencing techniques to analyse 454 blood samples from Afrotropical ducks from 5 different locations (spread around the subregion) for avian malaria. Fifty-five blood samples were positive for one or more genera of haematozoa. The regional infection rate across all sites and sampling periods was 12·1%. Nine individuals carried dual infections containing multiple haematozoa. Fifteen different cytochrome b haplotypes among 52 positives (3 samples failed to sequence) and 61 total sequences were found. Eleven haplotypes closely matched Plasmodium, whereas 4 were more similar to Haemoproteus. Five distinct haematozoan clades were identified. Haemoproteus parasites appeared to be more host-specific than Plasmodium, which occurred at every sampling location and in every host species examined. There were no significant differences in overall parasite prevalence attributable to either site or species, although Plasmodium and Haemoproteus occurrences differed by site-species combination and the borderline significance of our test for between-site variation (P < 0·06) implied that with a larger sample size, differences in parasite prevalence among locations might be detectable.

Genetic structure of the endangered northeastern bulrush (Scirpus ancistrochaetus) in Pennsylvania, USA, using information from RAPD and SNPs.

The northeastern bulrush, Scirpus ancistrochaetus, is a federally endangered wetland plant species found primarily in Pennsylvania, USA. Data on the population genetic structure of this species are needed by conservation managers to prioritize conservation efforts. In this study, we used two genetic marker systems to examine diversity and structure of this species in populations throughout Pennsylvania. The first simple and inexpensive approach utilized RAPD primers; our second, more detailed approach relied on DNA sequencing of single nucleotide polymorphisms. We found genetic variation using both RAPDs and sequencing and found some overlap in information between the two methods, including clusters of related populations and the identification of a genetically unique population. Future studies will seek to examine variation across the full geographic range of the species.

Behavioural vs. molecular sources of conflict between nuclear and mitochondrial DNA: the role of male-biased dispersal in a Holarctic sea duck.

Genetic studies of waterfowl (Anatidae) have observed the full spectrum of mitochondrial (mt) DNA population divergence, from apparent panmixia to deep, reciprocally monophyletic lineages. Yet, these studies often found weak or no nuclear (nu) DNA structure, which was often attributed to male-biased gene flow, a common behaviour within this family. An alternative explanation for this 'conflict' is that the smaller effective population size and faster sorting rate of mtDNA relative to nuDNA lead to different signals of population structure. We tested these alternatives by sequencing 12 nuDNA introns for a Holarctic pair of waterfowl subspecies, the European goosander (Mergus merganser merganser) and the North American common merganser (M. m. americanus), which exhibit strong population structure in mtDNA. We inferred effective population sizes, gene flow and divergence times from published mtDNA sequences and simulated expected differentiation for nuDNA based on those histories. Between Europe and North America, nuDNA Ф(ST) was 3.4-fold lower than mtDNA Ф(ST) , a result consistent with differences in sorting rates. However, despite geographically structured and monophyletic mtDNA lineages within continents, nuDNA Ф(ST) values were generally zero and significantly lower than predicted. This between- and within-continent contrast held when comparing mtDNA and nuDNA among published studies of ducks. Thus, male-mediated gene flow is a better explanation than slower sorting rates for limited nuDNA differentiation within continents, which is also supported by nonmolecular data. This study illustrates the value of quantitatively testing discrepancies between mtDNA and nuDNA to reject the null hypothesis that conflict simply reflects different sorting rates.

Multilocus coalescent analysis of haemoglobin differentiation between low- and high-altitude populations of crested ducks (Lophonetta specularioides).

Hypoxia is a key factor determining survival, and haemoglobins are targets of selection in species native to high-altitude regions. We studied population genetic structure and evaluated evidence for local adaptation in the crested duck (Lophonetta specularioides). Differentiation, gene flow and time since divergence between highland and lowland populations were assessed for three haemoglobin genes (α(A) , α(D) , ß(A) ) and compared to seven reference loci (six autosomal introns and mtDNA). Four derived amino acid replacements were found in the globin genes that had elevated Φ(ST) values between the Andean highlands and Patagonian lowlands. A single ß(A) -globin polymorphism at a site known to influence O(2) affinity was fixed for different alleles in the two populations, whereas three α(A) - and α(D) -globin polymorphisms exhibited high heterozygosity in the highlands but not in the lowlands. Coalescent analyses supported restricted gene flow for haemoglobin alleles and mitochondrial DNA but nonzero gene flow for the introns. Simulating genetic data under a drift-migration model of selective neutrality, the ß(A) -globin fell outside the 95% confidence limit of simulated data, suggesting that directional selection is maintaining different variants in the contrasting elevational environments, thereby restricting migration of ß(A) -globin alleles. The α(A) - and α(D) -globins, by contrast, did not differ from the simulated values, suggesting that variants in these genes are either selectively neutral, or that the effects of selection could not be differentiated from background levels of population structure and linkage disequilibrium. This study illustrates the combined effects of selection and population history on inferring levels of population divergence for a species distributed across an altitudinal gradient in which selection for hypoxia resistance has likely played an important role.

A parapatric propensity for breeding precludes the completion of speciation in common teal (Anas crecca, sensu lato).

Speciation is a process in which genetic drift and selection cause divergence over time. However, there is no rule dictating the time required for speciation, and even low levels of gene flow hinder divergence, so that taxa may be poised at the threshold of speciation for long periods of evolutionary time. We sequenced mitochondrial DNA (mtDNA) and eight nuclear introns (nuDNA) to estimate genomic levels of differentiation and gene flow between the Eurasian common teal (Anas crecca crecca) and the North American green-winged teal (Anas crecca carolinensis). These ducks come into contact in Beringia (north-eastern Asia and north-western North America) and have probably done so, perhaps cyclically, since the Pliocene-Pleistocene transition, ~2.6 Ma, when they apparently began diverging. They have diagnosable differences in male plumage and are 6.9% divergent in the mtDNA control region, with only 1 of 58 crecca and 2 of 86 carolinensis having haplotypes grouping with the other. Two nuclear loci were likewise strongly structured between these teal (Φ(st) ≥ 0.35), but six loci were undifferentiated or only weakly structured (Φ(st) = 0.0-0.06). Gene flow between crecca and carolinensis was ~1 individual per generation in both directions in mtDNA, but was asymmetrical in nuDNA, with ~1 and ~20 individuals per generation immigrating into crecca and carolinensis, respectively. This study illustrates that species delimitation using a single marker oversimplifies the complexity of the speciation process, and it suggests that even with divergent selection, moderate levels of gene flow may stall the speciation process short of completion.

Gene flow in the face of countervailing selection: adaptation to high-altitude hypoxia in the betaA hemoglobin subunit of yellow-billed pintails in the Andes.

When populations become locally adapted to contrasting environments, alleles that have high fitness in only one environment may be quickly eliminated in populations adapted to other environments, such that gene flow is partly restricted. The stronger the selection, the more rapidly immigrant alleles of lower fitness will be eliminated from the population. However, gene flow may continue to occur at unlinked loci, and adaptive divergence can proceed in the face of countervailing gene flow if selection is strong relative to migration (s > m). We studied the population genetics of the major hemoglobin genes in yellow-billed pintails (Anas georgica) experiencing contrasting partial pressures of oxygen in the Andes of South America. High gene flow and weak population subdivision were evident at seven putatively neutral loci in different chromosomal linkage groups. In contrast, amino acid replacements (Ser-beta13, Ser-beta116, and Met-beta133) in the betaA hemoglobin subunit segregated by elevation between lowland and highland populations with significantly elevated F(ST). Migration rates for the betaA subunit alleles were approximately 17-24 times smaller than for five unlinked reference loci, the alphaA hemoglobin subunit (which lacks amino acid replacements) and the mitochondrial DNA control region. The betaA subunit alleles of yellow-billed pintails were half as likely to be transferred downslope, from the highlands to the lowlands, than in the opposite direction upslope. We hypothesize that migration between the lowlands and highlands is restricted by local adaptation, and the betaA hemoglobin subunit is a likely target of selection related to high-altitude hypoxia; however, gene flow may be sufficiently high to retard divergence at most unlinked loci. Individuals homozygous for lowland alleles may have relatively little difficulty dispersing to the highlands initially but may experience long-term fitness reduction. Individuals homozygous for highland genotypes, in contrast, would be predicted to have difficulty dispersing to the lowlands if their hemoglobin alleles confer high oxygen affinity, predicted to result in chronic erythrocytosis at low elevation. Heterozygous individuals may have a dispersal advantage if their hemoglobin has a wider range of function due to the presence of multiple protein isoforms with a mixture of different oxygen affinities.

Prevalence and Coinfection of Three Tick-Borne Pathogens in Questing Adult Blacklegged Ticks Ixodes scapularis (Vilas County, Wisconsin).

Introduction: In North America, the blacklegged tick (Ixodes scapularis) is a vector of several human pathogens, and tick-borne disease incidence is increasing. Objectives: We estimated the prevalence of questing blacklegged ticks vectoring three zoonotic pathogens in Vilas County, Wisconsin. Materials and Methods: We collected 461 adult blacklegged ticks and used PCR to screen for the presence of pathogens that cause Lyme disease (Borrelia burgdorferi), human granulocytic anaplasmosis (HGA, Anaplasma phagocytophilum), and babesiosis (Babesia microti). Results: We found that 52.5% of ticks carried at least one pathogen. The estimated infection prevalence in the tick population was 17.4% (Lyme disease), 14.3% (HGA), and 6.5% (babesiosis). Multiple pathogens were present in 14.3% of ticks surveyed. Conclusion: About half of questing ticks tested in this study carried at least one zoonotic pathogen. Coinfection was common in our study area and, if not properly recognized, leads to greater risk of underdiagnosis.

Signatures of high-altitude adaptation in the major hemoglobin of five species of andean dabbling ducks.

Hypoxia is one of the most important factors affecting survival at high altitude, and the major hemoglobin protein is a likely target of selection. We compared population genetic structure in the alphaA and betaA hemoglobin subunits (HBA2 and HBB) of five paired lowland and highland populations of Andean dabbling ducks to unlinked reference loci. In the hemoglobin genes, parallel amino acid replacements were overrepresented in highland lineages, and one to five derived substitutions occurred at external solvent-accessible positions on the alpha and beta subunits, at alpha(1)beta(1) intersubunit contacts, or in close proximity to inositolpentaphosphate (IPP) binding sites. Coalescent analyses incorporating the stochasticity of drift and mutation indicated that hemoglobin alleles were less likely to be transferred between highland and lowland populations than unlinked alleles at five other loci. Amino acid replacements that were overrepresented in the highlands were rarely found within lowland populations, suggesting that alleles segregating at high frequency in the highlands may be maladaptive in the lowlands and vice versa. Most highland populations are probably nonmigratory and locally adapted to the Altiplano, but gene flow for several species may be sufficiently high to retard divergence at unlinked loci. Heterozygosity was elevated in the alphaA or betaA subunits of highland populations exhibiting high gene flow between the southern lowlands and the highlands and in highland species that disperse seasonally downslope to midelevation environments from the central Andean plateau. However, elevated heterozygosity occurred more frequently in the alphaA subunit but not simultaneously in both subunits, suggesting that selection may be more constrained by epistasis in the betaA subunit. Concordant patterns among multiple species with different evolutionary histories and depths of historical divergence and gene flow suggest that the major hemoglobin genes of these five dabbling duck species have evolved adaptively in response to high-altitude hypoxia in the Andes.

Old divergence and restricted gene flow between torrent duck (Merganetta armata) subspecies in the Central and Southern Andes.

AIM: To investigate the structure and rate of gene flow among populations of habitat-specialized species to understand the ecological and evolutionary processes underpinning their population dynamics and historical demography, including speciation and extinction. LOCATION: Peruvian and Argentine Andes. TAXON: Two subspecies of torrent duck (Merganetta armata). METHODS: We sampled 156 individuals in Peru (M. a. leucogenis; Chillón River, n = 57 and Pachachaca River, n = 49) and Argentina (M. a. armata; Arroyo Grande River, n = 33 and Malargüe River, n = 17), and sequenced the mitochondrial DNA (mtDNA) control region to conduct coarse and fine-scale demographic analyses of population structure. Additionally, to test for differences between subspecies, and across genetic markers with distinct inheritance patterns, a subset of individuals (Peru, n = 10 and Argentina, n = 9) was subjected to partial genome resequencing, obtaining 4,027 autosomal and 189 Z-linked double-digest restriction-associated DNA sequences. RESULTS: Haplotype and nucleotide diversities were higher in Peru than Argentina across all markers. Peruvian and Argentine subspecies showed concordant species-level differences (ΦST mtDNA = 0.82; ΦST autosomal = 0.30; ΦST Z chromosome = 0.45), including no shared mtDNA haplotypes. Demographic parameters estimated for mtDNA using IM and IMa2 analyses, and for autosomal markers using ∂a∂i (isolation-with-migration model), supported an old divergence (mtDNA = 600,000 years before present (ybp), 95% HPD range = 1.2 Mya to 200,000 ybp; and autosomal ∂a∂i = 782,490 ybp), between the two subspecies, characteristic of deeply diverged lineages. The populations were well-differentiated in Argentina but moderately differentiated in Peru, with low unidirectional gene flow in each country. MAIN CONCLUSIONS: We suggest that the South American Arid Diagonal was preexisting and remains a current phylogeographic barrier between the ranges of the two torrent duck subspecies, and the adult territoriality and breeding site fidelity to the rivers define their population structure.