Phylogeography of the South China field mouse (Apodemus draco) on the southeastern Tibetan Plateau reveals high genetic diversity and glacial refugia

The southeastern margin of the Tibetan Plateau (SEMTP) is a particularly interesting region due to its topographic complexity and unique geologic history, but phylogeographic studies that focus on this region are rare. In this study, we investigated the phylogeography of the South China field mouse, Apodemus draco, in order to assess the role of geologic and climatic events on the Tibetan Plateau in shaping its genetic structure. We sequenced mitochondrial cytochrome b (cyt b) sequences in 103 individuals from 47 sampling sites. In addition, 23 cyt b sequences were collected from GenBank for analyses. Phylogenetic, demographic and landscape genetic methods were conducted. Seventy-six cyt b haplotypes were found and the genetic diversity was extremely high (π = 0.0368; h = 0.989). Five major evolutionary clades, based on geographic locations, were identified. Demographic analyses implied subclade 1A and subclade 1B experienced population expansions at about 0.052-0.013 Mya and 0.014-0.004 Mya, respectively. The divergence time analysis showed that the split between clade 1 and clade 2 occurred 0.26 Mya, which fell into the extensive glacial period (EGP, 0.5-0.17 Mya). The divergence times of other main clades (2.20-0.55 Mya) were congruent with the periods of the Qingzang Movement (3.6-1.7 Mya) and the Kun-Huang Movement (1.2-0.6 Mya), which were known as the most intense uplift events in the Tibetan Plateau. Our study supported the hypothesis that the SEMTP was a large late Pleistocene refugium, and further inferred that the Gongga Mountain Region and Hongya County were glacial refugia for A. draco in clade 1. We hypothesize that the evolutionary history of A. draco in the SEMTP primarily occurred in two stages. First, an initial divergence would have been shaped by uplift events of the Tibetan Plateau. Then, major glaciations in the Pleistocene added complexity to its demographic history and genetic structure.     

Molecular phylogeny of Eastern Pacific porcelain crabs, genera Petrolisthes and Pachycheles, based on the mtDNA 16S rDNA sequence: phylogeographic and systematic implications

Porcelain crabs, genera Petrolisthes and Pachycheles, are diverse and abundant members of the eastern Pacific near-shore decapod crustacean community. Morphology-based taxonomic analyses of these crabs have determined groupings of affiliated species, but phylogenetic relationships remain unknown. We used sequence data from the mitochondrial 16S rRNA gene of 46 species of eastern Pacific porcelain crabs to perform phylogenetic analyses by distance and parsimony methods. Our results are used to compare the taxonomic significance of morphological and molecular characters, to examine sequence divergence rates of crab 16S rRNA genes, and to analyze the phylogeographic history of these crabs. Our phylogenetic trees indicate that the genus Petrolisthes is divided into two main clades, reflecting morphological features. One clade contains primarily tropical species, and the other contains species from throughout the eastern Pacific, as well as species in the genera Allopetrolisthes and Liopetrolisthes. Phylogenetic trees of Pachycheles suggest an antitropical distribution; north and south temperate species form one clade and tropical species form a second clade. Sequence divergence rates of the 16S rRNA gene from three pairs of geminate species can be used to date divergence times, and we discuss porcelain crab phylogeographic patterns in relation to paleogeographic events.     

The complete mitochondrial genome of Accipiter virgatus and evolutionary history of the pseudo-control regions in Falconiformes

The complete mitochondrial genome sequence of Accipiter virgatus was determined. This mt-genome was 17,952 bp in length and consisted of 22 tRNA genes, 13 protein-coding genes, 2 rRNA genes, one control region (CR) and one pseudo-control region (CCR). Phylogenetic analyses of 14,644 bp of mitochondrial DNA (12 protein-coding genes, 2 rRNAs and 22 tRNAs) revealed the phylogenetic position of Cathartidae (Cathartes aura) was more closer to Ciconiidae (storks) than Accipitridae. To investigate the divergence times of the CCRs in Falconiformes, detailed analyses of the noncoding regions (CR and CCR) were performed. We found the recently reported novel gene order in Falconiformes had multiple independent origins and hence cannot be used to infer phylogenetic lineages. Indeed, the molecular clock suggested the CCR in Falconidae emerged about 65.4 million years (Mya), while that in Pandionidae–Accipitridae clade emerged about 19.16 Mya. The intra-genomic homology between the noncoding regions was detected in Spilornis cheela, which supporting the duplication hypothesis. Furthermore, the structure of CCR should be featured by a region containing tandem repeats as two definitely separated clusters of tandem repeats were found. The findings presented here should be considered in future phylogenetic and evolutionary studies targeting the pseudo-control regions of all Falconiformes species.     

Dating of divergences within the Rattus genus phylogeny using whole mitochondrial genomes

The timing and order of divergences within the genus Rattus have, to date, been quite speculative. In order to address these important issues we sequenced six new whole mitochondrial genomes from wild-caught specimens from four species, Rattus exulans, Rattus praetor, Rattus rattus and Rattus tanezumi. The only rat whole mitochondrial genomes available previously were all from Rattus norvegicus specimens. Our phylogenetic and dating analyses place the deepest divergence within Rattus at approximately 3.5 million years ago (Mya). This divergence separates the New Guinean endemic R. praetor lineage from the Asian lineages. Within the Asian/Island Southeast Asian clade R. norvegicus diverged earliest at approximately 2.9Mya. R. exulans and the ancestor of the sister species R. rattus and R. tanezumi subsequently diverged at approximately 2.2Mya, with R. rattus and R. tanezumi separating as recently as approximately 0.4Mya. Our results give both a better resolved species divergence order and diversification dates within Rattus than previous studies.     

A phylogenetic study of the genus Apodemus by sequencing the mitochondrial DNA control region

The sequences of the mitochondrial DNA control region (D-loop) and flanking tRNA genes (about 1000 bp) of 20 samples of wood mice (genus Apodemus ) were analyzed in order to clarify the relationships between different species belonging to the genus. The phylogenetic trees obtained using different methods showed similar topologies with distinct Karstomys ( Apodemus epimelas and Apodemus mystacinus ) and Sylvaemus ( Apodemus alpicola , Apodemus flavicollis , Apodemus hermonensis , Apodemus sylvaticus and Apodemus uralensis ) subtrees. Within Sylvaemus all species appeared to be closely related to each other, probably as result of a bush-like radiation event. Nevertheless, A. hermonensis seemed to be the first diverging branch followed by A. sylvaticus ; A. alpicola and A. flavicollis appeared to be very closely related. Three individuals of uncertain taxonomical status were included in the analysis: hypotheses as to their status are discussed. Further phylogenetic analysis was carried out combining the D-loop sequences of part of the samples of certain taxonomical status with 12S rRNA and cytochrome b sequences obtained by other researchers. Furthermore, I present a structure analysis of the D-loop in Apodemus as compared other rodent species.     

Phylogeny of the genus Apodemus with a special emphasis on the subgenus Sylvaemus using the nuclear IRBP gene and two mitochondrial markers: cytochrome b and 12S rRNA

Phylogenetic relationships among 17 extant species of Murinae, with special reference to the genus Apodemus, were investigated using sequence data from the nuclear protein-coding gene IRBP (15 species) and the two mitochondrial genes cytochrome b and 12S rRNA (17 species). The analysis of the three genes does not resolve the relationships between Mus, Apodemus, and Rattus but separates Micromys from these three genera. The analysis of the two mitochondrial regions supported an association between Apodemus and Tokudaia and indicated that these two genera are more closely related to Mus than to Rattus or Micromys. Within Apodemus, the mitochondrial data sets indicated that 8 of the 9 species analyzed can be sorted into two main groups: an Apodemus group, with A. agrarius, semotus, and peninsulae, and a Sylvaemus group, with uralensis, flavicollis, alpicola, sylvaticus, and hermonensis. The position of Apodemus mystacinus is ambiguous and might be either included in Sylvaemus or considered a distinct subgenus, Karstomys, more closely related to Sylvaemus than to Apodemus. Estimation of the divergence time for these taxa suggests a separation between 7 and 8 My ago for the three groups (mystacinus and the two subgenera Apodemus and Sylvaemus). Within each subgenus, divergence times are between 5.4 and 6 My for Apodemus and between 2.2 and 3.5 My for Sylvaemus and mystacinus.     

Temporal diversification of Mesoamerican cichlid fishes across a major biogeographic boundary

The Mexican Neovolcanic Plateau sharply divides the vertebrate fauna of Mesoamerica where the climate of both the neotropics and temperate North America gradually blend. Only a few vertebrate groups such as the Heroine cichlids, distributed from South America to the Rio Grande in North America, are found both north and south of the Neovolcanic Plateau. To better understand the geography and temporal diversification of cichlids at this geologic boundary, we used mitochondrial DNA sequences of the cytochrome b (cyt b) gene to reconstruct the relationships of 52 of the approximately 80 species of Heroine cichlids in Mesoamerica. Our analysis suggests several cichlids in South America should be considered as part of the Mesoamerican Heroine clade because they and the cichlids north of the Isthmus of Panama are clearly supported as monophyletic with respect to all other Neotropical cichlids. We also recovered a group containing species in Paratheraps+Paraneetroplus+Vieja as the sister clade to Herichthys. Herichthys is the only cichlid clade north of the Mexican Plateau and it is monophyletic. Non-parametric rate smoothing of cichlid cyt b sequence resulted in an estimated divergence time of approximately 6 million years for Herichthys. This temporal diversification is concordant with divergence times estimated for anurans in the genus Bufo, a group that exhibits a similar geographic distribution. Our results indicate the 5-million-year-old extension of the Mexican Neovolcanic Plateau to the Gulf Coast of Mexico has strongly influenced the current transition between the vertebrate faunas of the Neotropics and Nearctic.     

Phylogenetic relationships in genus Niviventer (Rodentia: Muridae) in China inferred from complete mitochondrial cytochrome b gene

Chinese species of the genus Niviventer, predominantly distributed in the southeastern Tibetan Plateau and in Taiwan, are a diverse group and have not yet received a thorough molecular phylogenetic analysis. Here, we reconstructed the phylogenetic relationships of 32 specimens representing nine Chinese species of Niviventer, based on sequences of the complete mitochondrial cytochrome b gene. Maximum parsimony, maximum likelihood and Bayesian analysis resulted in three consistent trees, each supported by high bootstrap values. The results showed that the Niviventer species included here are monophyletic. The nine species were classified into three distinct clades: clade A with Niviventer brahma, N. confucianus, N. coxingi, N. culturatus, N. eha and N. fulvescens; clade B with N. andersoni and N. excelsior; clade C with N. cremoriventer. Our results also suggested that N. culturatus should be a valid species rather than a subspecies of N. confucianus. Divergence times among species were calibrated according to the middle-late Pleistocene (1.2-0.13 Mya) fossil records of N. confucianus. The results demonstrated that the first radiation event of the genus Niviventer occurred in early Pleistocene (about 1.66 Mya), followed by the divergence of clades A and B at about 1.46 Mya. Most of the extant Niviventer species appeared during early to middle Pleistocene (about 1.29-0.67 Mya). These divergence times are coincidental with the last uplift events of the Tibetan Plateau, Kun-Huang movement, Pleistocene glaciations and the vicariant formation of Taiwan Strait. Consequently geographical events and Pleistocene glaciations have played a great role in the diversification of Niviventer.     

Phylogeny of Paramysis (Crustacea: Mysida) and the origin of Ponto-Caspian endemic diversity: resolving power from nuclear protein-coding genes

The Ponto-Caspian (Black and Caspian seas) brackish-water fauna represents a special case of the endemic diversification in world's ancient lakes; it also involves a hotspot of continental diversity in the predominantly marine mysid crustaceans. We explored the origins and history of the mysid diversification in a phylogenetic analysis of some 20 endemic Ponto-Caspian species mainly of the genus Paramysis and their marine congeners, using sequences of two nuclear protein-coding genes, two nuclear rRNA genes, the mitochondrial COI gene and morphological data. A nearly completely resolved phylogeny was recovered, with no indication of rapid diversification bursts. Deep divergences were found among the main endemic clades, attesting to a long independent faunal history in the continental Paratethys waters. The current marine Paramysis species make a monophyletic cluster secondarily derived from the continental Paratethyan (Ponto-Caspian) Paramysis ancestors. The good phylogenetic resolution was mainly due to the two nuclear protein-coding genes, opsin and EPRS, here for the first time applied to peracarid systematics. In contrast, 'conventional' mtDNA and nuclear rRNA genes provided poor topological resolution and weak congruence of divergence rates. The two nuclear protein-coding genes had more congruent rates of evolution, and were about 10-15 times slower than the mitochondrial COI gene.     

Resolution of phylogenetic relationships among recently evolved species as a function of amount of DNA sequence: an empirical study based on woodpeckers (Aves: Picidae)

Synonymous substitutions in the 13 mitochondrial encoded protein genes form a large pool of characters that should approach the ideal for phylogenetic analysis of being independently and identically distributed. Pooling sequences from multiple mitochondrial protein-coding genes should result in statistically more powerful estimates of relationships among species that diverged sufficiently recently that most nucleotide substitutions are synonymous. Cytochrome oxidase I (COI) was sequenced for woodpecker species for which cytochrome b (cyt b) sequences were available. A pairing-design test based on the normal distribution indicated that cyt b evolves more rapidly than COI when all nucleotides are compared but their rates are equal for synonymous substitutions. Nearly all of the phylogenetically informative substitutions among woodpeckers are synonymous. Statistical support for relationships, as measured by bootstrap proportions, increased as the number of nucleotides increased from 1047 (cyt b) to 1512 (COI) to 2559 nucleotides (aggregate data set). Pseudo-bootstrap replicates showed the same trend and increasing the amount of sequence beyond the actual length of 2559 nucleotides to 5120 (2x) resulted in stronger bootstrap support, even though the amount of phylogenetic information was the same. However, the amount of sequence required to resolve an internode depends on the length of the internode and its depth in the phylogeny.     

Complete nucleotide sequences and gene organization of mitochondrial genome of Bufo gargarizans

The complete mitochondrial genome of Bufo gargarizans was sequenced using overlapping polymerase chain reaction (PCR) amplicons (GenBank Accession No. DQ275350). The genome is 17,277 base pairs in length, containing 13 protein-coding genes (ATP6, ATP8, COI-III, ND1-6, ND4L, Cyt b), 2 ribosomal RNAs (12S rRNA and 16S rRNA), 22 transfer RNAs and a putative control region. We analyzed the sequence using bioinformatics methods comparing the obtained mtDNA sequence with other toads and frogs. Based on the concatenated nucleotide sequences of protein-coding genes, we constructed a phylogenetic tree with maximum likelihood (ML) and maximum parsimony (MP) methods and discussed the phylogenetic relationships among 11 species of Anura.     

Phylogenetic systematics and biogeography of the poison frogs: evidence from mitochondrial DNA sequences

Regions of the mitochondrial genome were sequenced and analysed in representative species of poison frogs, in order to investigate phylogenetic relationships within the family Dendrobatidae. Mitochondrial DNA (mfDNA) fragments from three gene regions; cytochrome b, 16S ribosomal RNA (rRNA), and 12S rRNA, provided 1198 base pairs of DNA sequence and 589 informative sites. Phylogenetic analysis using parsimony was used to infer the evolutionary relationships among the species in the survey. Our analysis supported previous partitions of species into the genera Epipedobates, Phyllobates and Dendrobates , with two exceptions; Epipedobates (Allobates) femoralis was placed outside the clade containing the other toxic dendrobatids, and Minyobates minutus was placed within the genus Dendrobates. Genetic distances estimated between all pairs of taxa using the Kimura 2-parameter model indicated substantial genetic divergence between species, particularly those found in Amazonia. Time of divergence estimates were highly variable depending on gene region, but even the lowest estimates were inconsistent with the Pleistocene Refugia hypothesis.     

Genetic divergence and speciation in lowland and montane peruvian poison frogs

Amazonia is famous for high biodiversity, and the highlands of the transition zone between the Andes and the lowlands of the Amazon basin show particularly high species diversity. Hypotheses proposed to explain the high levels of diversity in the highlands include repeated parapatric speciation across ecological gradients spanning the transition zone, repeated allopatric speciation across geographic barriers between the highlands and lowlands, divergence across geographic barriers within the transition zone, and simple lineage accumulation over long periods of time. In this study, we investigated patterns of divergence in frogs of the genus Epipedobates (family Dendrobatidae) using phylogenetic and biogeographic analyses of divergence in mitochondrial DNA (1778 aligned positions from genes encoding cyt b, 12S and 16S rRNA for 60 Epipedobates and 11 outgroup specimens) and coloration (measured for 18 specimens representing nine species in Epipedobates). The majority of phenotypic and species diversity in the poison frog genus Epipedobates occurs in the transition zone, although two morphologically conserved members of the genus are distributed across the lowlands of the Amazon basin. Phylogenetic analysis reveals that there is a single highland clade derived from an ancestral colonization event in northern Peru by a population of lowland ancestry. Epipedobates trivittatus, a widespread Amazonian species, is a member of the highland clade that reinvaded the lowlands. Comparative analyses of divergence in coloration and mtDNA reveals that divergence in coloration among populations and species in the highlands has been accelerated relative to the lowlands. This suggests a role for selection in the divergence of coloration among populations and species.     

Phylogenetic relationships of leaf monkeys (Presbytis; Colobinae) based on cytochrome b and 12S rRNA genes

Little is known about the classification and phylogenetic relationships of the leaf monkeys (Presbytis). We analyzed mitochondrial DNA sequences of cytochrome b (Cyt b) and 12S rRNA to determine the phylogenetic relationships of the genus Presbytis. Gene fragments of 388 and 371 bp of Cyt b and 12S rRNA, respectively, were sequenced from samples of Presbytis melalophos (subspecies femoralis, siamensis, robinsoni, and chrysomelas), P. rubicunda and P. hosei. The genus Trachypithecus (Cercopithecidae) was used as an outgroup. The Cyt b NJ and MP phylogeny trees showed P. m. chrysomelas to be the most primitive, followed by P. hosei, whereas 12S rRNA tree topology only indicated that these two species have close relationships with the other members of the genus. In our analysis, chrysomelas, previously classified as a subspecies of P. melalophos, was not included in either the P. m. femoralis clade or the P. m. siamensis clade. Whether or not there should be a separation at the species level remains to be clarified. The tree topologies also showed that P. m. siamensis is paraphyletic with P. m. robinsoni, and P. m. femoralis with P. rubicunda, in two different clades. Cyt b and 12S rRNA are good gene candidates for the study of phylogenetic relationships at the species level. However, the systematic relationships of some subspecies in this genus remain unclear.     

Differential geographic patterns of mitochondrial DNA variation in two sympatric species of Japanese wood mice, Apodemus speciosus and A. argenteus

We examined the gene sequences of mitochondrial cytochrome b (cyt b) in two Japanese wood mouse species, Apodemus speciosus (n = 89) and A. argenteus (n = 46), which are distributed on the four main islands of Japan (Hokkaido, Honshu, Shikoku, and Kyushu) and on the small islands surrounding them. Apodemus speciosus, the larger of the two species, showed substantial genetic variation, with a maximum of 3% sequence divergence, and remarkable phylogenetic subdivision with two major clades. One clade represents haplotypes from a central region, including Honshu, Shikoku, Kyushu, and their adjacent islands; the other clade includes haplotypes from Hokkaido and the peripheral islands, forming four subclades: a) Hokkaido, b) Sado Island, c) Satsunan Islands, and d) the Izu Islands. Sequence divergence among the four subclades was 1.0 to 1.5%, implying that A. speciosus colonized these geographic regions 0.2 to 0.3 million years ago, assuming a substitution rate of 2.4% per million years. The population on the Izu Islands has preserved haplotypes that are distinct from those in any other region, providing good evidence for the natural colonization of the volcanic islands of the Izu Islands. The cyt b sequence variation had no relation to the karyotypic dimorphism for the eastern (2n = 48) and western (2n = 46) geographic groups, between which a strict border exists at central Honshu. On the other hand, Apodemus argenteus, the smaller of the two species, showed a similar level of sequence divergence (maximum of 3%) but no substantial geographic differentiation: populations in Hokkaido, Sado, and Yakushima shared similar haplotypes with each of the central populations, suggesting that genetic exchanges occurred among the localities in the last 0.15 million years. The apparent genetic structure of the mitochondrial DNA found in the A. speciosus population might be caused solely by long-term existence in insular regions, presumably due to ecological superiority relative to A. argenteus.     

Molecular phylogeny and divergence times of drosophilid species

The phylogenetic relationships and divergence times of 39 drosophilid species were studied by using the coding region of the Adh gene. Four genera--Scaptodrosophila, Zaprionus, Drosophila, and Scaptomyza (from Hawaii)--and three Drosophila subgenera--Drosophila, Engiscaptomyza, and Sophophora--were included. After conducting statistical analyses of the nucleotide sequences of the Adh, Adhr (Adh-related gene), and nuclear rRNA genes and a 905-bp segment of mitochondrial DNA, we used Scaptodrosophila as the outgroup. The phylogenetic tree obtained showed that the first major division of drosophilid species occurs between subgenus Sophophora (genus Drosophila) and the group including subgenera Drosophila and Engiscaptomyza plus the genera Zaprionus and Scaptomyza. Subgenus Sophophora is then divided into D. willistoni and the clade of D. obscura and D. melanogaster species groups. In the other major drosophilid group, Zaprionus first separates from the other species, and then D. immigrans leaves the remaining group of species. This remaining group then splits into the D. repleta group and the Hawaiian drosophilid cluster (Hawaiian Drosophila, Engiscaptomyza, and Scaptomyza). Engiscaptomyza and Scaptomyza are tightly clustered. Each of the D. repleta, D. obscura, and D. melanogaster groups is monophyletic. The splitting of subgenera Drosophila and Sophophora apparently occurred about 40 Mya, whereas the D. repleta group and the Hawaiian drosophilid cluster separated about 32 Mya. By contrast, the splitting of Engiscaptomyza and Scaptomyza occurred only about 11 Mya, suggesting that Scaptomyza experienced a rapid morphological evolution. The D. obscura and D. melanogaster groups apparently diverged about 25 Mya. Many of the D. repleta group species studied here have two functional Adh genes (Adh-1 and Adh-2), and these duplicated genes can be explained by two duplication events.      

Molecular phylogeny and biogeography of caecilians from Southeast Asia (Amphibia, Gymnophiona, Ichthyophiidae), with special reference to high cryptic species diversity in Sundaland

We investigated the phylogenetic relationships and estimated the history of species diversification and character evolution in two ichthyophiid genera: Caudacaecilia and Ichthyophis. We estimated the phylogenetic relationships of 67 samples from 33 localities in Southeast Asia from 3840-bp sequences of the mitochondrial 12S rRNA, 16S rRNA, and cyt b genes using Bayesian inference, maximum likelihood, and maximum parsimony methods. The Southeast Asian samples formed a well-supported clade differentiated from a South Asian sample. The Southeast Asian clade was divided into two subclades, one containing samples from South China, Indochina, Malay Peninsula, and Java. The other consisted of samples from Borneo and the Philippines. Neither Caudacaecilia nor Ichthyophis was monophyletic, nor did samples with or without light stripes lateral to the body form clades. We found several distinct sympatric lineages and undescribed species, especially from Sundaland.     

基于线粒体DNA序列探讨斑头鱼分类地位

测定了斑头鱼Agrammus agrammus和大泷六线鱼Hexagrammos otakii的线粒体COⅠ、Cyt b和16S rRNA基因的部分序列,结合从GenBank中获得的六线鱼属3种的同源序列,以单鳍多线鱼Pleurogrammus monopterygius为外群,运用邻接法(NJ)、最大简约法(MP)和贝叶斯法(BI)构建了分子系统树。同时联合了3个基因片段序列,运用贝叶斯联合模型综合探讨了六线鱼类的系统发育关系。结果显示:除16S rRNA基因外,其余2个基因片段以及联合模型所构建的系统树拓扑结构完全一致,即斑头鱼与大泷六线鱼亲缘关系最近,应归为六线鱼属,拉丁学名应为Hexagrammos agrammus。Cyt b基因片段序列分析结果显示,斑头鱼和大泷六线鱼分歧时间约为175万年。结合形态学研究资料,支持将斑头鱼归为六线鱼属的观点,斑头鱼和大泷六线鱼亲缘关系最近,属于六线鱼科中分化较晚的种类。     

Molecular phylogenetic analysis of circumnuclear hemoproteids (Haemosporida: Haemoproteidae) of sylviid birds, with a description of Haemoproteus parabelopolskyi sp. nov

Haemoproteus spp., with circumnuclear gametocytes and tentatively belonging to Haemoproteus belopolskyi, are widespread and prevalent in warblers belonging to the Sylviidae, with numerous mitochondrial cytochrome b (cyt b) lineages detected among them.We sampled the hemoproteids from 6 species of warblers adjacent to the Baltic Sea. Parasites were identified to species based on morphology of their gametocytes, and a segment of the parasite's cyt b gene was sequenced. Sixteen mitochondrial cyt b lineages of hemoproteids with circumnuclear gametocytes were recorded. Two clades of lineages (clade A in species of Acrocephalus and Hippolais and clade B in species of Sylvia) with sequence divergence between their lineages >5% are distinguished in the phylogenetic tree. Within the clades A and B, the genetic distance between the lineages is < or = 3.9 and < and = 2.8%, respectively. We compared the morphology of gametocytes of 3 lineages (hHIICT1, hMW1, and hSYAT2) in detail. The lineages hHIICTI and hMW1 (clade A) belong to the morphospecies H. belopolskyi. Parasites of the lineage hSYAT2 (clade B) are described as a new species Haemoproteus parabelopolskyi, which can be readily distinguished from H. belopolskyi by the significantly smaller nuclei of its macrogametocytes. Lineages closely related to H. belopolskyi and H. parabelopolskyi are identified. The sequence divergence between lineages of these 2 morphospecies ranges between 5.3 and 8.1%. It seems probable that avian Haemoproteus spp. with a genetic differentiation of > or =5% in mitochondrial cyt b gene might be morphologically differentiated at the stage of gametocytes. This study establishes the value of both PCR and morphology in identification of avian hemoproteids.