Phylogenetic relationships of horned lizards (Phrynosoma) based on nuclear and mitochondrial data: evidence for a misleading mitochondrial gene tree

It has proven remarkably difficult to obtain a well-resolved and strongly supported phylogeny for horned lizards (Phrynosoma) because of incongruence between morphological and mitochondrial DNA sequence data. We infer the phylogenetic relationships among all 17 extant Phrynosoma species using >5.1 kb of mtDNA (12S rRNA, 16S rRNA, ND1, ND2, ND4, Cyt b, and associated tRNA genes), and >2.2kb from three nuclear genes (RAG-1, BDNF, and GAPD) for most taxa. We conduct separate and combined phylogenetic analyses of these data using maximum parsimony, maximum likelihood, and Bayesian methods. The phylogenetic relationships inferred from the mtDNA data are congruent with previous mtDNA analyses based on fewer characters and provide strong support for most branches. However, we detected strong incongruence between the mtDNA and nuclear data using comparisons of branch support and Shimodaira-Hasegawa tests, with the (P. platyrhinos+P. goodei) clade identified as the primary source of this conflict. Our analysis of a P. mcalliixP. goodei hybrid suggests that this incongruence is caused by reticulation via introgressive hybridization. Our preferred phylogeny based on an analysis of the combined data (excluding the introgressed mtDNA data) provides a new framework for interpreting character evolution and biogeography within Phrynosoma. In the context of this improved phylogeny we propose a phylogenetic taxonomy highlighting four clades: (1) Tapaja, containing the viviparous short-horned lizards P. ditmarsi, P. hernandesi, P. douglasii, and P. orbiculare; (2) Anota, containing species with prominent cranial horns (P. solare, P. mcallii, and the P. coronatum group); (3) Doliosaurus, containing three species lacking antipredator blood-squirting (P. modestum, P. platyrhinos, and P. goodei); and (4) Brevicauda, containing two viviparous species with extremely short tails that lack blood-squirting (P. braconnieri and P. taurus).     

Horned lizard (Phrynosoma) phylogeny inferred from mitochondrial genes and morphological characters: understanding conflicts using multiple approaches

The genus Phrynosoma includes 13 species of North American lizards characterized by unique and highly derived morphologies and ecologies. Understanding interspecific relationships within this genus is essential for testing hypotheses about character evolution in this group. We analyzed mitochondrial ND4 and cytochrome b gene sequence data from all species of Phrynosoma in conjunction with a previously published dataset including 12S and 16S rRNA gene sequences and morphological characters. We used multiple phylogenetic methods and diagnostic tests for data combinability and taxonomic congruence to investigate the data in separate and combined analyses. Separate data partitions resulted in several well-supported lineages, but taxonomic congruence was lacking between topologies from separate and combined analyses. Partitioned Bremer support analyses also reveals conflict between data partitions in certain tree regions. When taxa associated with well-supported clades were removed from analyses, phylogenetic signal was lost. Combined, our results initially suggest conflict between data partitions, but further tests show the data are only appropriate for phylogenetic reconstruction of those parts of the topology that were well resolved. Nonetheless, our data analyses reveal five well-supported clades: (1) Phrynosoma ditmarsi and Phrynosoma hernandesi, (2) P. ditmarsi, P. hernandesi, and Phrynosoma douglasii, (3) P. ditmarsi, P. hernandesi, P. douglasii, and Phrynosoma orbiculare, (4) Phrynosoma mcallii and Phrynosoma platyrhinos, and (5) Phrynosoma braconnieri and Phrynosoma taurus.     

Comparative phylogeography of bagrid catfishes in Taiwan

A phylogeographic analysis of two bagrid catfishes in Taiwan was conducted using sequence data from a portion of the mitochondrial DNA (mtDNA) control region. For Pseudobagrus brevianalis, which is most probably endemic to Taiwan, a total of eight haplotypes were detected in 189 specimens from nine river systems covering its entire distribution range, from northern to central western areas of the island. Obvious genetic differentiation was observed among its populations (average F _ST = 0.753); in particular, the northernmost Tamsui River population was fixed for a single endemic haplotype. Nested clade phylogeographic analysis (NCPA) indicated that the dispersal center of mtDNA was the area around the Touchien River and Holong River, north to the Miaoli Plateau, in northwestern Taiwan, suggesting both northward and southward dispersal in this species. There was no evidence for the validity of P. taiwanensis, the nominal species described from Taiwan, morphologically similar to P. brevianalis. We confirmed that P. adiposalis was distributed discontinuously in three river systems; analysis of 42 specimens from the rivers indicated a total of four haplotypes and population differentiation (average F _ST = 0.876). Fixation into a largely differentiated haplotype in the northernmost Tamsui River population was also found in this species, but different processes for this phylogeographic pattern were implied for the two species. Comparison with P. ussuriensis, a widespread continental species morphologically similar to P. adiposalis, suggested the possibility that P. adiposalis is a group of local populations of P. ussuriensis in Taiwan. Two migration routes vs. random fixation scenarios for the population structure of P. adiposalis are discussed with information on other fishes and the geological history of the island.     

Phylogeography and population genetics of the maize stalk borer Busseola fusca (Lepidoptera, Noctuidae) in sub-Saharan Africa

The population genetics and phylogeography of African phytophagous insects have received little attention. Some, such as the maize stalk borer Busseola fusca, display significant geographic differences in ecological preferences that may be congruent with patterns of molecular variation. To test this, we collected 307 individuals of this species from maize and cultivated sorghum at 52 localities in West, Central and East Africa during the growing season. For all collected individuals, we sequenced a fragment of the mitochondrial cytochrome b. We tested hypotheses concerning the history and demographic structure of this species. Phylogenetic analyses and nested clade phylogeographic analyses (NCPA) separated the populations into three mitochondrial clades, one from West Africa, and two--Kenya I and Kenya II--from East and Central Africa. The similar nucleotide divergence between clades and nucleotide diversity within clades suggest that they became isolated at about the same time in three different refuges in sub-Saharan Africa and have similar demographic histories. The results of mismatch distribution analyses were consistent with the demographic expansion of these clades. Analysis of molecular variance (amova) indicated a high level of geographic differentiation at different hierarchical levels. NCPA suggested that the observed distribution of haplotypes at several hierarchical levels within the three major clades is best accounted for by restricted gene flow with isolation by distance. The domestication of sorghum and the introduction of maize in Africa had no visible effect on the geographic patterns observed in the B. fusca mitochondrial genome.     

南门峡裂腹鱼亚科鱼类形态相似种的分类学地位——形态趋同进化实例

受高原抬升所致的水系变迁及人类活动的影响,分布于南门峡河流的裂腹鱼亚科鱼类与黄河干支流种群间的基因交流受到长期限制。作为孤立小群体,探讨其分类学地位及其在小生境中的进化机制对了解青藏高原鱼类多样性和物种的形成、进化具有重要意义。本文采用聚合酶链式反应(PCR)和直接测序方法获得了南门峡裂腹鱼亚科鱼类(n=29)及其近缘种(n=19)共48个个体的线粒体DNA(mtDNA)细胞色素b(cytb)基因的全序列(1140bp),并以厚唇裸重鱼和尖裸鲤为外群构建了MP和Bayesian系统进化树。南门峡裂腹鱼亚科鱼类29个个体的序列经排序后,发现有100个(8·77%)多态性位点,共定义了16个单倍型,在系统进化树上分布于截然不同的两个族群中。其中5个单倍型(NMX3、6、7、13、15)与其近缘种花斑裸鲤和青海湖裸鲤形成单系群(MP99%,Bayesian98%),而其余11个单倍型(NMX1、2、4、5、8、9、10、11、12、14、16)与黄河干支流的黄河裸裂尻鱼形成另一个单系群(MP99%,Bayesian99%)。序列差异分析显示,分布于不同族群的南门峡裂腹鱼亚科鱼类之间存在较大的碱基差异(平均为7·42%),显示出种间差异水平,表明分布于南门峡河流的裂腹鱼亚科鱼类可能是花斑裸鲤和黄河裸裂尻鱼形态相似种的复合体。结合青藏高原隆升所致的气候环境变化和高原北部水系变迁的事件,推断形态趋同进化可能导致了南门峡河流裂腹鱼亚科鱼类形态相似种的共存,而小生境自然选择压力是引发适应性形态趋同进化的主要原因。     

Morphological correlates of ant eating in horned lizards (Phrynosoma)

The North American horned lizards ( Phrynosoma ) represent a morphologically specialized group of ant-eating lizards. Although variation in dietary fidelity is observed among the species, all appear to possess morphological specializations thought to be related to their ant-eating diets. Previous studies have examined morphological specialization in Phrynosoma , but they have not taken into account the phylogenetic relationships of its member species. In the present study, the morphological characteristics of the head, jaws and teeth that are thought to be important in prey capture and prey processing were examined to test whether variation in cranial morphology is associated with diet in lizards of the genus Phrynosoma . It is suggested that lizards of the genus Phrynosoma are indeed morphologically specialized and that ant-eating is associated with reduced dentition and an overall reduction in the robustness of morphological structures important in prey processing. Although this trend holds for the highly myrmecophagous species of Phrynosoma , a robust cranial morphology is apparent in the short-horned lizard clade ( Phrynosoma ditmarsi , Phrynosoma douglasii , Phrynosoma hernandesi , Phrynosoma orbiculare ), implying the ability to process a variety of dietary items. The present study suggests that additional feeding specializations exist within an already specialized clade (i.e. the short-horned lizard clade) and highlights the need for more detailed dietary and behavioural studies of feeding behaviour in this uniquely specialized group of lizards.  © 2006 The Linnean Society of London, Biological Journal of the Linnean Society , 2006, 89 , 13–24.     

Sea‐level changes and palaeo‐ranges: reconstruction of ancient shorelines and river drainages and the phylogeography of the Australian land crayfish Engaeus sericatus Clark (Decapoda: Parastacidae)

Historical sea levels have been influential in shaping the phylogeography of freshwater‐limited taxa via palaeodrainage and palaeoshoreline connections. In this study, we demonstrate an approach to phylogeographic analysis incorporating historical sea‐level information in a nested clade phylogeographic analysis (NCPA) framework, using burrowing freshwater crayfish as the model organism. Our study area focuses on the Bass Strait region of southeastern Australia, which is marine region encompassing a shallow seabed that has emerged as a land bridge during glacial cycles connecting mainland Australia and Tasmania. Bathymetric data were analysed using Geographical Information Systems (GIS) to delineate a palaeodrainage model when the palaeocoastline was 150 m below present‐day sea level. Such sea levels occurred at least twice in the past 500 000 years, perhaps more often or of larger magnitude within the last 10 million years, linking Victoria and Tasmania. Inter‐locality distance measures confined to the palaeodrainage network were incorporated into an NCPA of crayfish (Engaeus sericatus Clark 1936) mitochondrial 16S rDNA haplotypes. The results were then compared to NCPAs using present‐day river drainages and traditional great‐circle distance measures. NCPA inferences were cross‐examined using frequentist and Bayesian procedures in the context of geomorphological and historical sea‐level data. We found distribution of present‐day genetic variation in E. sericatus to be partly explained not only by connectivity through palaeodrainages but also via present‐day drainages or overland (great circle) routes. We recommend that future studies consider all three of these distance measures, especially for studies of coastally distributed species.     

THE ROLE OF INTROGRESSIVE HYBRIDIZATION IN THE EVOLUTION OF THE GILA ROBUSTA COMPLEX (TELEOSTEI: CYPRINIDAE)

Abstract The extent and impact of introgressive hybridization was examined in the Gila robusta complex of cyprinid fishes using mitochondrial DNA (mtDNA) sequence variation. Lower Colorado River basin populations of G. robusta, G. elegans, and G. cypha exhibited distinct mtDNAs, with only limited introgression of G. elegans into G. cypha. The impact of hybridization was significant in upper Colorado River basin populations; most upper basin fishes sampled exhibited only G. cypha mtDNA haplotypes, with some individuals exhibiting mtDNA from G. elegans. The complete absence of G. robusta mtDNA, even in populations of morphologically pure G. robusta, indicates extensive introgression that predates human influence. Analysis of the geographic distribution of variation identified two distinctive G. elegans lineages; however, the small number of individuals and localities sampled precluded a comprehensive analysis. Analysis of haplotype and population networks for G. cypha mtDNAs from 15 localities revealed low divergence among haplotypes; however, significant frequency differences among populations within and among drainages were found, largely attributable to samples in the Little Colorado River region. This structure was not associated with G. cypha and G. robusta, as morphotypes from the same location are more similar than conspecific forms in other locations. This indicates that morphological and mtDNA variation are affected by different evolutionary forces in Colorado River Gila and illustrates how both hybridization and local adaptation can play important roles in evolution.     

Genetic detection of multiple exotic water frog species in Belgium illustrates the need for monitoring and immediate action

The recent genetic screening of water frogs (genus Pelophylax) in Belgium has shown that the invasion of two water frog species from Eastern Europe and the Mediterranean region, P. ridibundus and P. cf. bedriagae is widespread. Possibly other exotic water frogs are invading and establishing themselves through commercial trade. We used a genetic identification approach to rapidly detect and identify morphologically cryptic exotic water frog species in a large number of populations throughout the northern part of Belgium. Among a total of 944 individuals, we found 506 non-indigenous specimens, seven of which belonged to species not recorded before with certainty in Belgium or neighbouring regions in the wild. One of them was identified genetically as the Iberian green frog (P. perezi), but was most likely a P. perezi × P. esculentus hybrid. Six individuals of the Levantine frog (P. bedriagae) were found in a pond in the vicinity of a pet shop where the species is sold. All other exotic frogs belonged to four different haplotypes of P. ridibundus, established in Belgium since c. 1970, and two haplotypes of P. cf. bedriagae, a poorly-known eastern Mediterranean sister species of the latter. Overall, our study underscores the extent of exotic water frog invasions associated with the pet trade. Although two of the exotic species were recorded in small numbers, their early detection is essential with regards to adequate control and eradication of invasive species.     

Using COI gene sequence to barcode two morphologically alike species: the cotton bollworm and the oriental tobacco budworm (Lepidoptera: Noctuidae)

Due to limited morphological difference, the two closely related sister species, the cotton bollworm, Helicoverpa armigera (Hübner) and the oriental tobacco budworm, H. assulta (Guenée) (Lepidoptera: Noctuidae), are very difficult to distinguish, especially at the larvae stage. Recently, DNA sequence has been widely used as a bio-barcode for species identification. In this study, we attempted to distinguish H. armigera and H. assulta using the mitochondrial cytochrome C oxidase subunit I gene (COI) gene sequence as the barcode. We determined a 658 bp segment of the COI gene for 28 individuals of H. armigera, 8 individuals of H. assulta, and 10 individuals of Mamestra brassicae (as the outgroup) in Yunnan Province, southwest of P. R. China, together with one H. assulta and two H. armigera reported sequences from GenBank. Twenty-three haplotypes were identified in all 49 samples. As expected, network analysis of the haplotypes of the three species presented a clustering pattern consistent with the respective species status. Haplotypes of the same species differed from each other by no more than three nucleotide substitutions. However, each haplotype of H. armigera differed from that of H. assulta by at least 22 nucleotide substitutions. Both species differed from M. brassicae by more than 50 nucleotide substitutions. 17 unique diagnostic nucleotides were also used to discriminate the two species. The finding of large COI sequence differences between H. armigera and H. assulta suggested that it could be used to distinguish the two morphologically alike species and be employed for quick species identification during pest control.     

Novel patterns of historical isolation, dispersal, and secondary contact across Baja California in the Rosy Boa (Lichanura trivirgata)

Mitochondrial DNA (mtDNA) sequence variation was examined in 131 individuals of the Rosy Boa (Lichanura trivirgata) from across the species range in southwestern North America. Bayesian inference and nested clade phylogeographic analyses (NCPA) were used to estimate relationships and infer evolutionary processes. These patterns were evaluated as they relate to previously hypothesized vicariant events and new insights are provided into the biogeographic and evolutionary processes important in Baja California and surrounding North American deserts. Three major lineages (Lineages A, B, and C) are revealed with very little overlap. Lineage A and B are predominately separated along the Colorado River and are found primarily within California and Arizona (respectively), while Lineage C consists of disjunct groups distributed along the Baja California peninsula as well as south-central Arizona, southward along the coastal regions of Sonora, Mexico. Estimated divergence time points (using a Bayesian relaxed molecular clock) and geographic congruence with postulated vicariant events suggest early extensions of the Gulf of California and subsequent development of the Colorado River during the Late Miocene-Pliocene led to the formation of these mtDNA lineages. Our results also suggest that vicariance hypotheses alone do not fully explain patterns of genetic variation. Therefore, we highlight the importance of dispersal to explain these patterns and current distribution of populations. We also compare the mtDNA lineages with those based on morphological variation and evaluate their implications for taxonomy.     

Polyphyly and species delimitation of Picea brachytyla (Pinaceae) based on population genetic data

Accurate species delimitation is fundamental to biodiversity conservation. The endangered spruce Picea brachytyla (Franch.) E. Pritz. was suggested to be polyphyletic based on a limited number of samples in previous studies. To evaluate polyphyly of P. brachytyla, we sampled 139 individuals from 16 populations across most of its distributional range, plus representatives of two related species, Picea likiangensis (Franch.) E. Pritz. and Picea wilsonii Mast. We sequenced 13 nuclear loci and three chloroplast and two mitochondrial loci for the following species delimitation. Phylogenetic analyses of nuclear loci grouped all individuals of P. brachytyla from Sichuan and Chongqing into one distinct lineage and those from Yunnan and Tibet (southern distribution) nested within the P. likiangensis species complex. Structure analyses confirmed this result. Networks of chloroplast DNA haplotypes similarly showed that P. brachytyla from the southern distribution nested within the P. likiangensis species complex, whereas haplotypes for the northern distribution comprised a separate and well-supported lineage. These results suggest that P. brachytyla from the southern distribution is a part of the P. likiangensis species complex and should be removed from P. brachytyla. Our study highlights the utility of population genetic evidence in delimitating endangered species and understanding the conservation status of such species.     

GENETIC VARIATION OF KOGIA SPP. WITH PRELIMINARY EVIDENCE FOR TWO SPECIES OF KOGIA SIMA

Concordance between mitochondrial DNA (mtDNA) markers and morphologically based species identifications was examined for the two currently recognized Kogia species. We sequenced 406 base pairs of the control region and 398 base pairs of the cytochrome b gene from 108 Kogia breviceps and 47 K. sima samples. As expecred, the two sister species were reciprocally monophyletic to each other in phylogenetic reconstructions, but within K. sima , we unexpectedly observed another reciprocally monophyletic relationship. The two K. sima clades resolved were phylogeographically concordant with all of the haplotypes in one clade observed solely among specimens sampled from the Atlantic Ocean and with those in the other clade observed solely among specimens sampled from the Indo-Pacific Ocean. These apparently allopatric clades were observed in all phylogenetic reconstructions using the maximum parsimony, maximum likelihood, and neighborjoining algorithms, with the mtDNA gene sequences analyzed separately and combined. The nucleotide diversity for the combined gene sequence haplotypes of the two K. sima clades resolved in our analyses was 0.58% and 1.03% for the Atlantic and Indo-Pacific, respectively, whereas for the two recognized sister species, nucleotide diversity was 1.65% and 4.02% for K. breviceps and K. sima , respectively. The combined gene sequence haplotypes have accumulated 44 fixed base pair differences between the two K. sima clades compared to 20 fixed base pair differences between the two recognized sister species. Although our results are consistent with species-level differences between the two K. sima clades, recognition of a third Kogia species awaits supporting evidence that these two apparently allopatric clades represent reproductively isolated groups of animals.     

Batbugs (Cimex pipistrelli group,Heteroptera: Cimicidae) are morphologically,but not genetically differentiated among bat hosts

Members of the family Cimicidae (Heteroptera) are obligate haematophagous ectoparasites. The Cimex pipistrelli species group parasitizes on bats, the likely ancestral hosts of the whole family. Based on morphology, it was suggested that three species of the group were present in the West‐Palaearctic region, although their validity remained a matter of discussion. Surprisingly, the status of these species has not been studied from the point of view of host specificity. We examined the diversification of the species group using morphological data, including the putative diagnostic characters, and sequences of one mitochondrial (cytochrome oxidase subunit I, COI) and four nuclear loci (internal transcribed spacer 2, 18S and 28s ribosomal genes and elongation factor 1 subunit α). This was carried out on a sample of 225 individuals from 69 bat roosts and 44 mist‐netted bats, altogether representing 12 bat species from 13 European countries and Lebanon. We revealed 27 mitochondrial haplotypes representing two distinct haplogroups and one outlying haplotype. The extent of morphological variability of specimens representing both haplogroups covers the range of characters reported for all three recognized species; therefore, the haplogroups clearly do not correspond to any described species. Also, the very limited variability found in the nuclear sequences of the cimicid bugs examined suggests that separate species do not exist in the region. We found considerable morphological differentiation among samples from different bat species, although individuals representing particular mitochondrial haplogroups often live sympatrically and on the same host species. It seems that batbugs are morphologically adapted to a particular bat host despite the low genetic structuring among individuals parasitizing different species of bats.     

Phylogeography of northern Dolly Varden Salvelinus malma (Salmoniformes: Salmonidae) from Asia and North America: An analysis based on the mitochondrial DNA genealogy

Spawning in habitats affected by Pleistocene glacial advances over most of its natural range, northern Dolly Varden Salvelinus malma malma typifies Arctic fauna distributed in northeastern Asia and northwestern North America. We reconstructed a genealogy of mtDNA haplotypes from 27 Alaskan and Asian populations to study the influence of historical events on the phylogeography and contemporary population genetic structure. Analysis of molecular variance partitioned most of the mtDNA variability to the intrapopulation component (72.5%) with much reduced differences between populations (21.1%) and regions (6.4%). Similar patterns of variation apparent from hierarchical diversity and nested clade phylogeographical analysis (NCPA) of mtDNA haplotypes identify weak spatial differentiation and low levels of divergence. These findings suggest (1) that demographic history has been influenced by historical range expansions and recent isolation by distance, (2) that present populations from Asia and North America were colonized from one main Beringian Refugium, and (3) that this taxon’s ancestral population probably experienced a bottleneck in the Beringian Refugium during the late Pleistocene (Wisconsin) glacial period. The genealogical and NCPA analyses, and mismatch distribution of S. m. malma mtDNA haplotypes do not confirm the assumptions about presence of the two refugia on the territories of the Beringian Land, in which allopatric S. m. malma ancestral populations evolved, and independent origin of the Sea of Okhotsk populations.     

DNA barcoding for identification of sand flies (Diptera: Psychodidae) in India

About 50 species of sand flies have been reported to be prevalent in India. We explored the utility of the DNA barcode approach towards species identification of these medically important insects. A total of 62 specimens belonging to seven morphologically identified species of two genera, Phlebotomus and Sergentomyia, collected from Puducherry Union Territory, Maharashtra and Rajasthan states of India were subjected to the analysis. Neighbor‐joining (NJ) analysis of DNA barcode sequences identified the individuals of seven morphological species into eight distinct species, as presented in the designed NJ tree. This methodology delineated morphologically identified species, S. bailyi, into two genetically isolated groups. Also, this study characterizes DNA barcodes of P. argentipes and P. papatasi, the vector species of leishmaniasis in India, for the first time.     

Monitoring of the Flat-Tailed Horned Lizard With Methods Incorporating Detection Probability

Abstract: Difficulty in monitoring the flat-tailed horned lizard (Phrynosoma mcallii) has led to controversy over its conservation status. The difficulty in detecting this species has discouraged large-scale estimates of abundance and led to uncertainty over whether the species exists in population sizes of sufficient size for long-term persistence. We incorporated detection probability into monitoring of this species using closed mark—recapture and distance-sampling methods. Density estimation from mark—recapture abundance estimates was improved using an estimate of the proportion of time lizards were on the plot. We estimated the probability of detection on the line for distance sampling and adjusted density estimates accordingly. We estimated the populations of the Yuha Basin Management Area in 2002 and the East Mesa Management Area, Imperial County, California, USA, in 2003 to be 25,514 (95% CI 14,444-38,970) and 42,619 (95% CI 23,161-67,639), respectively. Two estimates of detection probability on the line in distance sampling by different methods were 0.45 and 0.65. Density estimates derived from distance analyses for 3 East Mesa Management Area plots and the Yuha Basin Management Area were 1.55 per ha (95% CI 0.64-3.76) and 0.41 per ha (95% CI 0.22-0.7), respectively. These are the first large-scale estimates of abundance and density for P. mcallii.     

Molecular systematics of short-horned lizards: biogeography and taxonomy of a widespread species complex

We surveyed mitochondrial DNA (mtDNA) sequence variation in short-horned lizards (Phrynosoma douglasi) from throughout western North America and used these data to estimate an intraspecific phylogeny and to assess biogeographic scenarios underlying the geographic structure of lineages in this species. We sequenced 783 base pairs from 38 populations of P. douglasi and three putative outgroups (P. ditmarsi, P. orbiculare, P. platyrhinos). We detected high levels of nucleotide variation among populations and a spatial distribution of mtDNA lineages compatible with major geographic regions. The phylogenetic hypotheses best supported by the data suggest that P. douglasi, as currently described, is paraphyletic with respect to P. ditmarsi. Populations of P. douglasi from the Pacific Northwest (ID, CA, OR, WA) form a monophyletic group that is sister to the subsequent radiation of P. ditmarsi and other P. douglasi clades. These results suggest that divergences within this widespread species are fairly old. We focused on the genetic structure of populations of P. douglasi from a geographic perspective and interpreted the intraspecific phylogeny in light of geologic and climatic changes in western North America during the last 20 million years. The generally high levels of genetic variation found in these population comparisons are in accord with high levels of morphological variation in this species group; however, only in the Pacific Northwest region is there spatial congruence between these phylogenetic results and subspecific ranges based on previous morphological studies. We compared the evolutionary units delineated in this study with previously described subspecies of P. douglasi and evaluated the support (from morphology and mtDNA) for each population lineage in the phylogeny and the implications for the taxonomy of this group.     

On the specificity of avian blood parasites: revealing specific and generalist relationships between haemosporidians and biting midges

The study of host-parasite relationships involving vector-borne parasites requires understanding interactions between parasites and vectors. The capacity of haemosporidians to infect insects has clear evolutionary consequences for the transmission of diseases. Here, we investigated (i) the associations between blood parasites, biting midges and birds and (ii) the potential specificity between biting midge and haemosporidian haplotypes. A total of 629 parous biting midges Culicoides and 224 wild birds (belonging to seven species) from a locality of central Spain were individually examined for the presence of Haemoproteus and Plasmodium parasites by sequencing a fragment of cytochrome B. Biting midges were identified morphologically and characterized on the basis of a fragment of the cytochrome c oxidase (COI) gene. Overall, 12 Haemoproteus and three Plasmodium haplotypes were isolated and sequenced. Among them, 10 haplotypes were exclusively isolated from biting midges, three haplotypes only from birds and two haplotypes from both biting midges and birds. Biting midge haplotypes showed both specific and generalist relationships with Haemoproteus haplotypes but only generalist relationships with Plasmodium haplotypes. Several C. festivipennis and C. kibunesis haplotypes established significant coevolutionary links with Haemoproteus haplotypes. These results shed light on the specificity of interactions between vectors and blood parasites.