Nuclear DNA contents, rDNAs, and karyotype evolution in subgenus Vicia: III. The heterogeneous section Hypechusa

Summary. Nuclear DNA contents, automated karyotype analyses, and sequences of internal transcribed spacers from ribosomal genes have been determined in the species belonging to section Hypechusa of the subgenus Vicia. Karyomorphological results and phylogenetic data generated from the comparison of rDNA (genes coding for rRNA) sequences showed that sect. Hypechusa is not monophyletic; however, some monophyletic units are apparent (one including Vicia galeata, V. hyrcanica, V. noeana, and V. tigridis, another including V. assyriaca, V. hybrida, V. melanops, V. mollis, and V. sericocarpa), which partly correspond to morphology-based infrasectional groups. The relationships among these species and the species in sections Faba, Narbonensis, Bithynicae, and Peregrinae have been also investigated. Correspondence and reprints: Dipartimento di Agrobiologia ed Agrochimica, Università della Tuscia, via San Camillo de Lellis, 01100 Viterbo, Italy.     

Phylogenetic relationships between and within sections Hypechusa, Narbonensis and Peregrinae of genus Vicia (Fabaceae) based on evidence from isozymes and morphology

Previous studies of vetch species belonging to sections Hypechusa, Narbonensis and Peregrinae have revealed controversies and open questions that deserve revisiting with the use of contemporary approaches. The paper extends the previous work of Leht and Jaaska (2002) by describing variation of 16 isozymes and 86 morphological characters among 17 Vicia species belonging to sections Hypechusa, Narbonensis and Peregrinae of the type subgenus in comparison with V. bithynica of section Pseudolathyrus. Three more isozymes and 12 new morphological characters are included to improve the resolution of phylogenetic relationships. The maximum parsimony trees based on isozymes and morphology revealed essentially the same four clusters of species that correspond to series Hypechusa and Hyrcanicae of section Hypechusa, to sections Narbonensis and Peregrinae of current taxonomy, except that the isozyme tree placed V. lutea of series Hypechusa into a separate clade. Isozyme variation in V. mollis, V. noeana, V. sericocarpa and V. ciliatula of section Hypechusa is described for the first time. Twenty-three species-specific isozymes that may be used for species identification by analysis of seeds are described. Two new isozymes that distinguish V. serratifolia from V. narbonensis are described that further support their specific status. Series Hypechusa and Hyrcanicae of section Hypechusa, sections Narbonensis and Peregrinae of traditional taxonomy are monophyletic groups that deserve sectional status within the genus Vicia. The phylogenetic position and sectional placement of V. lutea warrants further study with independent sets of DNA characters.     

Cladistic and phenetic analysis of relationships in Vicia subgenus Vicia (Fabaceae) by morphology and isozymes

 Variation of 80 multistate morphological characters and isozymes encoded by 13 loci among 23 vetch species of the type subgenus of the genus Vicia in comparison with V. dumetorum, V. pisiformis and V. sylvatica of the subgenus Cracca is described and analyzed with cladistic parsimony and phenetic neighbour-joining methods by using two different ways of coding. Morphological analyses showed the subgenus Vicia monophyletic and revealed subgroups in a general agreement with traditionally recognized sections, except showing V. faba nested within section Narbonensis and ambiguity in the position of V. lathyroides and V. bithynica. Parsimony analysis of orthozymes as presence/absence characters revealed in the subgenus two basic monophyletic clades: 1) V. faba and three species of the section Peregrinae, V. michauxii, V. aintabensis and V. peregrina, in one subclade linked with species of the Narbonensis and Hyperchusa sections together with V. pisiformis of subgenus Cracca in a second subclade; 2) species belonging to sections Vicia, Sepium, Pseudolathyrus and Lathyroides together with V. sylvatica of the subgenus Cracca. Neighbour-joining analysis of orthozymes revealed the same two basic groups, differing only in the relative position of some species in subclusters. Both isozyme analyses showed paraphyly of the subgenera Vicia and Cracca. Parsimony analysis of orthozymes as character states of isozymes yielded a largely unresolved strict consensus cladogram of 209 most parsimonious trees, and reweighting of characters failed to produce a stable tree. Phylogenetic congruence and discordance among morphological and isozyme analyses, coding ways, homoplasy and weighting of characters are discussed. Received November 20, 2001 Accepted January 31, 2002     

Molecular phylogeny of Russulaceae (Basidiomycetes; Russulales) inferred from the nucleotide sequences of nuclear large subunit rDNA

Phylogenetic relationships of the genera Russula and Lactarius were investigated using sequence data from the nuclear-encoded large subunit ribosomal DNA (LSU rDNA). Ninety-five sequences belonging to the genera Russula and Lactarius, including 31 sequences from the databases, were used in this study. Analysis of the LSU rDNA region indicated that Russulaceae was divided into six groups (group A–F) in the neighbor-joining (NJ) tree. Lactarius consisted of one large clade (group A). Therefore, this genus was found to be monophyletic. However, the monophyly of genus Russula remained unclear. The genus Russula consisted of five groups in the NJ tree. Group B includes sects. Plorantes and Archaeinae (Heim), and group C includes sects. Delicoarchaeae and Russula in the NJ tree. Neither of the two groups formed a single clade in the most parsimonius (MP) tree. Group D includes many taxa having colored spore prints and amyloid in suprahilar plage of spores in sect. Russula and sect. Rigidae. Group E consists of only sect. Compactae and is further divided into three subclades, represented by R. densifolia, R. nigricans, and R. subnigricans, respectively. Group F contains sects. Rigidae, Ingratae, and Pelliculariae. Sect. Compactae and sect. Plorantes should not be as closely related as previously supposed. Russula earlei may be placed in sect. Archaeinae Heim. Russula flavida (subsect. Amoeninae) is placed in sect. Russula with R. aurea with a high bootstrap value (99%). The nuclear LSU rDNA region is a useful tool in recognization of species of Russulaceae and may provide information concerning phylogenetic relationships between the genera Russula and Lactarius.     

Microarray-based survey of repetitive genomic sequences in Vicia spp.

A modified DNA microarray-based technique was devised for preliminary screening of short fragment genomic DNA libraries from three Vicia species (V. melanops, V. narbonensis, and V. sativa) to isolate representative highly abundant DNA sequences that show different distribution patterns among related legume species. The microarrays were sequentially hybridized with labeled genomic DNAs of thirteen Vicia and seven other Fabaceae species and scored for hybridization signals of individual clones. The clones were then assigned to one of the following groups characterized by hybridization to: (1) all tested species, (2) most of the Vicia and Pisum species, (3) only a few Vicia species, and (4) preferentially a single Vicia species. Several clones from each group, 65 in total, were sequenced. All Group I clones were identified as rDNA genes or fragments of chloroplast genome, whereas the majority of Group II clones showed significant homologies to retroelement sequences. Clones in Groups III and IV contained novel dispersed repeats with copy numbers 10²–10⁶/1C and two genus-specific tandem repeats. One of these belongs to the VicTR-B repeat family, and the other clone (S12) contains an amplified portion of the rDNA intergenic spacer. In situ hybridization using V. sativa metaphase chromosomes revealed the presence of the S12 sequences not only within rDNA genes, but also at several additional loci. The newly identified repeats, as well as the retroelement-like sequences, were characterized with respect to their abundance within individual genomes. Correlations between the repeat distributions and the current taxonomic classification of these species are discussed.     

Isozyme evidence on the specific distinctness and phylogenetic position of <Emphasis Type="Italic">Vicia incisa</Emphasis> (Fabaceae)

Vicia incisa is a taxonomically controversial species that has been also treated as a subspecies of V. sativa because of a great morphological similarity. The phylogenetic position of V. incisa is uncertain because various DNA markers have provided contradictory results. Isozymes of V. incisa encoded by 15 loci and resolved with the use of polyacrylamide gel electrophoresis (PAGE) are described and compared with those of seven related species belonging to sections Vicia, Sepium, Lathyroides and Pseudolathyrus in order to get new evidence about its taxonomic rank and phylogenetic position. Phylogenetic relationships are analyzed with maximum parsimony and neighbour joining methods. Vicia incisa is shown to differ from all three subspecies of V. sativa including, sativa, cordata and nigra, by alternate variants of ten isozymes out of 15 analysed. Instead, V. incisa has much more similarity to V. grandiflora and V. sepium by sharing eight isozyme variants which differ from the subspecies of V. sativa. The most parsimony and neighbour joining analyses of the isozyme variation placed V. incisa as basally linked to the V. grandiflora and V. sepium species couple in the clade of section Sepium (= sect. Atossa), while the subspecies of V. sativa together with V. lathyroides formed a separate clade of section Vicia. The isozyme data provide further support to the species status of V. incisa.     

Physical mapping of 18S-5.8S-26S and 5S ribosomal RNA gene families in three important vetches (Vicia species) and their allied taxa constituting three species complexes

The most-important vetch species, Vicia narbonensis (narbon vetch, section Faba), Vicia villosa (hairy vetch, section Cracca) and Vicia sativa (common vetch, section Vicia) and their close relatives (often difficult to circumscribe into distinct taxa) constitute respectively, Narbonensis, Villosa and Sativa species complexes in the genus Vicia. The distribution of the 18S-5.8S-26S (18S-26S) and 5S ribosomal RNA (rRNA) gene families on the chromosomes of 19 (2n=2x=10,12,14) of the 24 species and subspecies belonging to the three species complexes, and Vicia bithynica (2n=12, section Faba) and Vicia hybrida (2n=12, section Hypechusa) was studied by fluorescence in situ hybridization (FISH) with pTa 71 (18S-26S rDNA) and pTa 794 (5S rDNA) DNA clones. Computer – aided chromosome analysis was performed on the basis of chromosome length, the arm-length ratio and the position of the hybridization signals. The positions of the four (2+2) signals of the two rRNA gene families were similar between each of the three, as well as two subspecies of V. narbonensis and Vicia johannis, respectively. Two major 18S-26S rDNA loci were found in the nucleolus organiser regions (NORs) of each of the species except V. hybrida, where it was present in two out of four SAT chromosomes. In addition to major NORs, two minor loci have been physically mapped at the centromeric regions of chromosomes of group 1 in Vicia amphicarpa, Vicia macrocarpa and V. sativa, and two NORs of group 5 in V. hybrida, and on the long arms of group 4 in V. bithynica. Two or four 5S rDNA loci, observed in the short arms of groups 2–4 and 5, and 18S-26S rDNA loci were located in different chromosomes of all the species within the Narbonensis and Villosa species complexes, and Vicia angustifolia of the Sativa species complex. In the remaining six species of the Sativa species complex, and V. bithynica and V. hybrida, the two or four 5S rDNA sites were present in chromosomes which harbor 18S-26S rRNA genes. The tandemly repeated 5S rDNA sites, located at the proximal part of the long arm of groups 3–5, were diagnostic for V. angustifolia, Vicia cordata, Vicia incisa, V. macrocarpa, Vicia nigra and V. sativa of the Sativa species complex. In V. amphicarpa of the same complex, the tandem repeats were located at the distal part of the long arms of group 3. Variability in the number, size and location of two ribosomal DNA probes could generally distinguish species within the Narbonensis and Sativa species complex, V. bithynica and V. hybrida. With respect to the four species of the Villosa species complex the karyotypes could not be identified individually on the basis of the distribution of two ribosomal gene families in three out of seven pairs of chromosomes. Received: 18 October 2000 / Accepted: 20 March 2001     

Cytology of Vicia species. X. Karyotype evolution and phylogenetic implication in Vicia species of the sections Atossa, Microcarinae, Wiggersia and Vicia

Automated karyotype analyses and sequence of rDNA spacers have been analysed for the species belonging to sections Atossa, Microcarinae, Wiggersia and Vicia. Karyomorphological parameters, based on Rec, Syi and TF% indices, have been determined and evidenced that, in term of symmetry, the karyotype of Vicia lathyroides was the most asymmetric one. A multivariate analysis using 34 karyological parameters, in addition to the symmetry indices, has been carried out and the corresponding dendrogram of linkage distances showed six different groups. Molecular investigations on the inclusive group in study by employing ITS DNA sequences indicated a different pattern of relationships. The cladistic analysis combining the molecular data set with karyological parameters evidenced that the species of sections Vicia and Atossa join closely to each other in a paraphyletic group, which includes the monophyletic section Wiggersia. Therefore, our karyological and molecular data provide information about the phylogenetic position of the analysed species inside the subgenus Vicia and are discussed in relation to previous results obtained by morphology, isozymes and ribosomal genes analyses.     

Phylogeny of Carpinus and subfamily Coryloideae (Betulaceae) based on chloroplast and nuclear ribosomal sequence data

Phylogenetic studies were conducted for Carpinus and the subfamily Coryloideae (Betulaceae) using sequences of the chloroplast matK gene, the trnL-trnF region (trnL intron, and trnL [UAA] 3' exon-trnF [GAA] intergenic spacer) and the psbA-trnH intergenic spacer, and the nuclear ribosomal ITS regions. The combined analyses of the three chloroplast regions suggest that Coryloideae is monophyletic; Ostryopsis is sister to the Carpinus - Ostrya clade; Corylus is monophyletic and sister to the Ostrya - Carpinus - Ostryopsis clade; Ostrya is paraphyletic; and within Carpinus, species of sect. Carpinus from eastern Asia form a monophyletic group, whereas the positions of C. betulus from Europe and C. caroliniana from eastern North America are unresolved within the Carpinus clade. The cpDNA tree generated in this study is largely congruent with the previously published ITS results, but the ITS tree places Carpinus sect. Distegocarpus as sister to the Ostrya - Carpinus sect. Carpinus clade. Future work is needed to examine the relationships within the Ostrya - Carpinus clade, evaluate the generic status of Ostrya, and test the phylogenetic position of Ostryopsis.     

Phylogenetic study of benthic, spine-bearing prorocentroids, including Prorocentrum fukuyoi sp. nov.

Species of prorocentroid dinoflagellates are common in marine benthic sediment and epibenthic habitats, as well as in planktonic habitats. Marine planktonic prorocentroids typically possess a small spine in the apical region. In this study, we describe a new, potentially widely distributed benthic species of Prorocentrum, P. fukuyoi sp. nov., from tidal sand habitats in several sites in Australia and from central Japan. This species was found to possess an apical spine or flange and was sister species to P. emarginatum. We analyzed the phylogeny of the group including this new species, based on large subunit (LSU) rDNA sequences. The genus contained a high level of divergence in LSU rDNA, in some cases among sister taxa. P. fukuyoi and P. emarginatum were found to be most closely related to a clade of generally planktonic taxa. Several morphological features may constitute more informative synapomorphies than habitat in distinguishing clades of prorocentroid species.     

Pollen morphology of some taxa of Vicia L. subgenus Vicia (Fabaceae) from Turkey

In this study, the pollen morphology of 11 taxa belonging to Atossa (Alef.) Asch. & Graebner, Hypechusa (Alef.) Asch. & Graebner, Peregrinae Kupicha, Wiggersia (Alef.) Maxted, Vicia L. and Narbonensis (Radzhi) Maxted sections of the genus Vicia L. subgenus Vicia (Fabeae, Fabaceae) naturally growing in Turkey has been studied using Light Microscopy (LM) and Scanning Electron Microscopy (SEM) to evaluate the taxonomic relevance of pollen characters. Twelve morphometric characters are analysed with one-way analysis of variance and Tukey’s honestly significant difference test for multiple comparisons. Of the taxa studied V. narbonensis var. narbonensis (sect. Narbonensis) has the largest pollen grains (P = 51.98 μm × E = 30.52 μm) and V. lathyroides (sect. Wiggersia) has the smallest pollen grains (P = 27.71 μm × E = 20.14 μm). The pollen grains are subprolate to perprolate (P/E = 1.16–2.11), but the prolate shape occurs in the majority of the taxa. The regular pollen grains of all taxa are trizonocolporate, isopolar, and released in monads. Ornamentation of the mesocolpium is psilate-perforate in V. lathyroides (sect. Wiggersia), reticulate-rugulate in V. narbonensis var. narbonensis (sect. Narbonensis), (micro)reticulate in V. sericocarpa var. sericocarpa (sect. Hypechusa), V. sativa subsp. sativa (sect. Vicia) and V. grandiflora var. grandiflora (sect. Vicia), and reticulate-perforate in the remaining taxa. The apocolpium and colpus area are psilate or perforate in all taxa except V. sericocarpa var. sericocarpa (sect. Hypechusa) and V. sativa subsp. sativa (sect. Vicia), which exhibit the obscurely reticulate-perforate pattern. Several palynological features have taxonomic importance in sectional level: polar axis, equatorial diameter, pollen shape (P/E ratio), colpus length, colpus width, size of pori, porus length/width ratio, lumina diameter, muri thickness and ornamentation. The results also indicate that pollen characters can be useful in distinguishing the examined taxa.     

Nuclear DNA contents, rDNAs, chromatin organization, and karyotype evolution inVicia sect,faba

Summary Automated karyotype analyses, nuclear DNA contents, and sequences of rDNA internal transcribed spacers of the nine species inVicia sect.faba are reported. As karyomorphological parameters are used to evaluate the karyotype evolution, so the determination of the heterochromatin by Feulgen absorption at different thresholds of optical density provided further evidence on the chromatin organization withinVicia sect.faba. The comparison of sequences of rDNA spacers has enabled the definition of the phylogenetic relationships between the analyzed species.     

Phylogeny of Oedogoniales (Chlorophyceae, Chlorophyta) inferred from 18S rDNA sequences with emphasis on the relationships in the genus Oedogonium based on ITS-2 sequences

The phylogeny of Oedogoniales was investigated by using nuclear 18S rDNA sequences. Results showed that the genus Oedocladium, as a separated clade, was clustered within the clade of Oedogonium; whereas the genus Bulbochaete was in a comparatively divergent position to the other two genera. The relationship among the species of Oedogonium was discussed, focusing on ITS-2 phylogeny analyzed combining with some morphological characteristics. Our results showed that all the dioecious nannandrous taxa involved in this study were resolved into one clade, while all the monocious taxa were clustered into another clade as a sister group to the former. The report also suggests that the dioecious macrandrous taxa form a paraphyly and could be more basally situated than the dioecious nannandrous and the monoecious taxa by means of molecular phylogeny and morphotype investigations.     

Cytological and molecular characterization of Vicia esdraelonensis Warb. & Eig: a rare taxon

Summary. Vicia esdraelonensis, a rare taxon belonging to section Hypechusa of subgenus Vicia, was recovered and analyzed by cytological, karyological, and molecular methods, with the aim of both characterizing this species and furthering our knowledge of the phylogeny of subgenus Vicia. Automated karyotype analysis, nuclear DNA content, and chromatin organization were determined by the Feulgen reaction, as well as chromosome banding after double staining with 4′,6-diamidino-2-phenylindole (DAPI) and chromomycin A3. The chromosome number and the nuclear DNA content were in agreement with the values of the species of section Hypechusa. The GC- and AT-rich preferential sites were determined by chromomycin A3 and DAPI staining. Karyomorphological parameters indicated that V. esdraelonensis is in an intermediate position in the spatial representation of the species of section Hypechusa on the basis of symmetry indices, as well as in the dendrogram of linkage distance constructed on 37 chromosome parameters. Molecular data based on internal transcribed spacer sequences show that V. esdraelonensis can doubtlessly be included in section Hypechusa and document its closeness to V. noeana. A cladistic analysis combining the molecular data set with karyological characters is also reported. Karyological, cytological, and molecular data allow characterization of the V. esdraelonensis genome and provide information about the phylogenetic position of this species within the Hyrcanicae series of section Hypechusa. Correspondence and reprints: Dipartimento di Biologia, Università di Pisa, Via L. Ghini 5, 56126 Pisa, Italy.     

Taxonomic relationships between V. faba and its relatives based on nuclear and mitochondrial RFLPs and PCR analysis

Summary The taxonomic relationships between 52 accessions of 12 Vicia species and three accessions of Lathyrus were examined using nuclear RFLP- and PCR-generated data. Two hundred and sixty informative restriction fragments or amplification products were analysed by single linkage analysis, average cluster analysis, and the Wagner parsimony method. Dendrograms constructed from each type of analysis showed similar overall topologies and could be divided into three parts corresponding respectively to the Lathyrus outgroup, the species grouped in the section Faba/narbonemis complex, and the species belonging to the sections Hypechusa and Peregrinae. With few exceptions, the majority of accessions belonging to one species grouped together before branching to other species. An analysis of mitochondrial DNA phenotypes was both consistent with and complemented the results from the nuclear data. Overall, the species relationships show a good correlation with the classification of Maxted et al. but suggest that V. faba is more closely aligned to species from the sections Hypechusa and Peregrinae than to those in the narbonensis complex. In addition, the position of two new species, V. kalakhensis and V. eristaloides, as members of the narbonensis complex was supported by the molecular data, which also allowed a preliminary classification for recently collected Vicia accessions.     

Phylogenetic relationships among annual and perennial species of the genus <Emphasis Type="Italic">Cicer</Emphasis> as inferred from ITS sequences of nuclear ribosomal DNA

The cladistic analysis of the DNA sequences of the internal transcribed spacers of ribosomal cistrons (ITS1 and ITS2) for 20 species of Cicer L. (among which all the annuals), shows that various sections of the genus are not monophyletic. Annual species do not form a clade: C. arietinum, in fact, is closely related to both C. echinospermum and C. reticulatum, whereas C. bijugum, C. judaicum, and C. pinnatifidum form a separate clade. The annual C. cuneatum is sister group to the perennial C. canariense and both are archaic species within the genus. C. yamashitae is, on the contrary, the only annual species belonging to a group of perennials, within which close relationships are evident between C. graecum and C. montbretii as well as among a group of mainly Asian species.     

New clade of silicified bolidophytes that belong to Triparma (Bolidophyceae,Stramenopiles)

The small phytoplankton genus Triparma belongs to the class Bolidophyceae and contains two distinct forms: silicified species and naked flagellated species (formerly Bolidomonas). Recent studies showed that four silicified species/strains (Triparma laevis f. inornata, T. laevis f. longispina, T. strigata, and T. aff. verrucosa) belong to a single clade that is paraphyletic, because it also contains an unclassified flagellated strain, and is sister to a flagellated species, T. eleuthera. In this study, we isolated and characterized two new strains of silicified species to test the phylogenetic unity of silicified bolidophytes. The isolates were identified as T. retinervis strains because they possessed fine areolation on the cell wall. 18S rDNA and rbcL phylogenetic analyses demonstrated that T. retinervis formed a new silicified clade that is sister to the flagellated species T. pacifica. This reveals that there are at least two distinct clades including both silicified and flagellated Triparma species.     

The phylogenetic relationships of Cynoglottis (Boraginaceae- Boragineae) inferred from ITS, 5.8S and trnL sequences

A molecular phylogenetic analysis of Cynoglottis was performed to evaluate previous hypotheses based on non-molecular evidence concerning the position of this genus within Boraginaceae tribe Boragineae. ITS-5.8S and trnL_UAA sequences from the nuclear and chloroplast non-coding genomes were obtained for four Cynoglottis taxa and selected members of the related genera Anchusa, Anchusella, Gastrocotyle, Brunnera and Pentaglottis. Cynoglottis is monophyletic, but neither trnL nor ITS support a close relationship with Brunnera, unlike previously supposed on morphological grounds. Brunnera is, instead, related to the southwestern European monotypic genus Pentaglottis, with which it forms a basal clade. ITS-5.8S sequences show that Anchusa thessala, a southeastern European annual species of Anchusa subg. Buglossellum, is sister to Cynoglottis and that the two taxa form a clade which also includes the Balkan endemic Gastrocotyle macedonica. Species of Anchusa subg. Anchusa form a separate lineage with high bootstrap support, suggesting that this heterogeneous genus is paraphyletic with respect to Cynoglottis. ITS sequences also discriminate between the Balkan-Apenninic diploid C. barrelieri and the Anatolian tetraploid C. chetikiana, albeit with low support. The molecular results are discussed in the light of karyological, morphological and chorological aspects.This work has been supported by M.I.U.R. 40% 2003 and the University of Firenze.     

Molecular phylogeny of black flies (Diptera: Simuliidae) from Thailand, using ITS2 rDNA

The sequences of the second internal transcribed spacer (ITS2) of ribosomal DNA (rDNA) were determined for 40 black fly species from Thailand, belonging to 4 subgenera of the genus Simulium, namely Gomphostilbia (12 species), Nevermannia (5 species), Montisimulium (1 species), Simulium sensu stricto (21 species), and an unknown subgenus with one species (Simulium baimaii). The length of the ITS2 ranged from 247 to 308 bp. All black fly species had high AT content, ranging from 71 to 83.8%. Intraindividual variation (clonal variation) occurred in 13 species, ranging from 0.3 to 1.1%. Large intrapopulation and interpopulation heterogeneities exist in S. feuerboni from the same and different locations in Doi Inthanon National Park, northern Thailand. Phylogenetic relationships among 40 black fly species were examined using PAUP (version 4.0b10) and MrBAYS (version 3.0B4). The topology of the trees revealed two major monophyletic clades. The subgenus Simulium and Simulium baimaii were placed in the first monophyletic clade, whereas the subgenera Nevermannia + Montisimulium were placed as the sister group to the subgenus Gomphostilbia in the second monophyletic clade. Our results suggest that S. baimaii belongs to the malyschevi-group or variegatum-group in the subgenus Simulium. The molecular phylogeny generally agrees with existing morphology-based phylogenies.