The diversity ofMalassezia yeasts confirmed by rRNA sequence and nuclear DNA comparisons

One hundred and fourMalassezia strains (52 isolated from humans and 52 from animals) were compared using large subunit (LSU) ribosomal RNA sequence similarity and nuclear DNA complementarity. Eight groups of strains were recognized as genetically distinct species. Each taxon was confirmed by a homogeneous mole % GC and percentages of DNA/DNA reassociations higher than 85%. The non-lipid-dependentMalassezia yeasts were maintained as the unique taxonM. pachydermatis. In contrast, lipid-dependent strains were shown to be distributed among seven species:M. furfur, M. sympodialis andM. species 1–5. These taxa matched remarkably well with morphological and serological differences documented by previous investigators. The LSU rRNA sequences allowed a further intraspecific resolution with most of genomic taxa represented by several closely related sequences:M. pachydermatis counted up to seven sequences,M. furfur four sequences,M. species 1 comprised three sequences andM. species 2 andM. species 5 two sequences. Three species,M. sympodialis, M. species 3 andM. species 4, displayed a unique type of sequence. Thus, the present report demonstrates the usefulness of sequencing for both taxonomic and epidemiological purposes.     

DNA homologies among homothallic, pseudo-homothallic and heterothallic species of Neurospora.

Summary Highly purified DNAs from three homothallic speciesNeurospora africana, N. dodgei andN. lineolata; three reference strains representing authentic heterothallic species,N. crassa, N. intermedia andN. sitophila; and two strains of pseudo-homothallic speciesN. tetrasperma were characterized by spectrophotometry and DNA reassociation using hydroxyapatite chromatography. All of these known species are closely related on the basis of DNA characteristics such as base composition and thermal denaturation profiles of major DNA components. Minor components of ascospore DNA was, however, only 5–7% of total DNA instead of 15–20% minor component DNA shown by mycelial DNA. Species belonging to same group were not distinguishable morphologically, but all of these species were distinguishable by DNA:DNA homology studies. Greater DNA homology was noticed between DNAs of heterothallic species and DNAs of pseudohomothallic species than DNA of true homothallic species. Difference on DNA-nucleiotide sequences among homothallic species was very little. Pseudo-homothallic speciesN. tetrasperma was found to be distinctly different from homothallic species but closer to heterothallic species based on such studies.Supported in part by a contract No. E(40-1)4182 with the U.S. Energy Research and Development Administration. We are grateful to Departments of Oncology and Radiotherapy, College of Medicine, Howard University, Washington, D.C. for providing us with material assistance     

Molecular phylogeny of Lactarius volemus and its allies inferred from the nucleotide sequences of nuclear large subunit rDNA

The phylogenetic relationships of Lactarius volemus and its relatives were investigated using the nucleotide sequences of the nuclear-encoded large subunit ribosomal DNA (LSU rDNA). Thirty-six sequences from L. volemus, L. corrugis, and L. hygrophoroides, including three sequences obtained from the GenBank database, were used in this study. Samples studied were divided into four major subclades (A–D) in both neighbor-joining (NJ) and maximum-parsimony (MP) trees. Lactarius volemus and L. corrugis formed one large clade in both NJ and MP trees (bootstrap value, 100%), which was divided into three subclades (A–C). Subclade A included three clusters of L. volemus strains, i.e., velvet, red, and Chinese types. Subclade B included the common and red types of L. corrugis. Subclade C included the common and yellow types of L. volemus. Subclade D included only one type of L. hygrophoroides. An analysis of the fatty acid composition supported the divisions found in the molecular analysis. Analyses of nucleotide sequence, fatty acid composition, morphological characteristics, and the taste of the fruiting bodies all led us to conclude that the common, velvet, red, and Chinese types of L. volemus, and the common and red types of L. corrugis, may each belong to different species, subspecies, or varieties. Further studies with more material from a wide range of regions are required to conduct taxonomic revision of these types. The LSU rDNA region may be a useful tool to investigate phylogenetic relationships within the section Dulces of the genus Lactarius.     

Characterization of <Emphasis Type="Italic">Penicillium</Emphasis> Species by Ribosomal DNA Sequencing and BOX,ERIC and REP-PCR Analysis

The genus Penicillium is one of the largest and widely distributed fungal genera described to date. As a result, its taxonomic classification and species discrimination within this genus has become complicated. In this study, 52 isolates that belonged to the Penicillum genus and other related genera were characterized using two DNA-based methods: (i) analysis of the nucleotide sequences of internal transcribed spacers in ribosomal DNA and (ii) analysis of DNA fingerprints that were generated by polymerase chain reactions with specific primers for enterobacterial repetitive intergenic consensus (ERIC) and repetitive extragenic palindromic (REP) sequences, and BOX elements. Using both methods, Penicillium species were discriminated from other fungal genera. Furthermore, Penicillium species that include strains which are used as biocontrol agents, such as P. glabrum, P. purpurogenum, and P. oxalicum, could be distinguished from other Penicillium species using these techniques. Based on our findings, we propose that a polyphasic approach that includes analysis of the nucleotide sequences of ribosomal DNA and detecting the presence of highly conserved, repeated nucleotide sequences can be used to determine the genetic relationships between different Penicillium species. Furthermore, we propose that our results can be used as a start point to develop a strategy to monitor the environmental presence of particular strains of Penicillium species when they are used as biocontrol agents.     

Toxic Pseudo-nitzschia multistriata (Bacillariophyceae) from the Gulf of Naples: morphology,toxin analysis and phylogenetic relationships with other Pseudo-nitzschia species

The genus Pseudo-nitzschia includes several species capable of producing domoic acid, the causative agent of Amnesic Shellfish Poisoning. Some of these species have been recorded frequently in the Gulf of Naples. For one of the species, P. multistriata, which has been recurrently found in our sampling area since 1995, this is the first report for European waters. Here we provide further details on the fine structure of this species. Pseudo-nitzschia multistriata was the only one found to produce domoic acid among all the Pseudo-nitzschia species from the Gulf of Naples, and this finding raises the number of potentially toxic species in this genus to nine. Phylogenetic relationships among several Pseudo-nitzschia species were assessed using the hypervariable domains (D1–D3) of the large subunit (LSU) rDNA. The match between the phylogeny obtained and important taxonomic characters used in this genus are discussed. Results show that P. multistriata clusters with wider species lacking a central larger interspace in the raphe. Close genetic relationships were determined between P. fraudulenta and P. subfraudulenta, and between P. pungens and P. multiseries. Genetic differences among these pairs of species are comparable to those among isolates of P. pseudodelicatissima from the Gulf of Naples, indicating high intraspecific genetic diversity of Pseudo-nitzschia species in the relatively conserved LSU region. This could explain the problematic results obtained when testing a match between species-specific Pseudo-nitzschia LSU probes and our sequences.     

Characteristics and utility of nuclear-encoded large-subunit ribosomal gene sequences in phylogenetic studies of red algae

Primer sequences are described for amplifying and sequencing a large fragment (approximately 2500 b.p.) of the nuclear-encoded large-subunit ribosomal RNA gene (LSU) from red algae. In comparison to RuBisCo large-subunit gene (rbcL) and nuclear-encoded small-subunit ribosomal RNA gene (SSU) sequence data, LSU sequence data was intermediate in the number of phylogenetically informative positions and sequence divergence. Parsimony analysis of LSU sequences for 16 Gelidiales species resolved some nodes unresolved in rbcL and SSU parsimony trees. An analysis of LSU sequences from 13 species of red algae classified in 11 orders suggests that this gene may be useful in studies of higher-level relationships of red algae.     

New taxa refresh the phylogeny and classification of pleurostomatid ciliates (Ciliophora,Litostomatea)

A high diversity of pleurostomatid ciliates has been discovered in the last decade, and their systematics needs to be improved in the light of new findings concerning their morphology and molecular phylogeny. In this work, a new genus, Protolitonotus gen. n., and two new species, Protolitonotus magnus sp. n. and Protolitonotus longus sp. n., were studied. Furthermore, 19 novel nucleotide sequences of SSU rDNA, LSU rDNA and ITS1‐5.8S‐ITS2 were collected to determine the phylogenetic relationships and systematic positions of the pleurostomatid ciliates in this study. Based on both molecular and morphological data, the results demonstrated that: (i) as disclosed by the sequence analysis of SSU rDNA, LSU rDNA and ITS1‐5.8S‐ITS2, Protolitonotus gen. n. is sister to all other pleurostomatids and thus represents an independent lineage and a separate family, Protolitonotidae fam. n., which is defined by the presence of a semi‐suture formed by the right somatic kineties near the dorsal margin of the body; (ii) the families Litonotidae and Kentrophyllidae are both monophyletic based on both SSU rDNA and LSU rDNA sequences, whereas Amphileptidae are non‐monophyletic in trees inferred from SSU rDNA sequences; and (iii) the genera Loxophyllum and Kentrophyllum are both monophyletic, whereas Litonotus is non‐monophyletic based on SSU rDNA analyses. ITS1‐5.8S‐ITS2 sequence data were used for the phylogenetic analyses of pleurostomatids for the first time; however, species relationships were less well resolved than in the SSU rDNA and LSU rDNA trees. In addition, a major revision to the classification of the order Pleurostomatida is suggested and a key to its families and genera is provided.     

Phylogenetic relationships of Steinernema Travassos, 1927 (Nematoda: Cephalobina: Steinernematidae) based on nuclear, mitochondrial and morphological data

Entomopathogenic nematodes of the genus Steinernema are lethal parasites of insects that are used as biological control agents of several lepidopteran, dipteran and coleopteran pests. Phylogenetic relationships among 25 Steinernema species were estimated using nucleotide sequences from three genes and 22 morphological characters. Parsimony analysis of 28S (LSU) sequences yielded a well-resolved phylogenetic hypothesis with reliable bootstrap support for 13 clades. Parsimony analysis of mitochondrial DNA sequences (12S rDNA and cox 1 genes) yielded phylogenetic trees with a lower consistency index than for LSU sequences, and with fewer reliably supported clades. Combined phylogenetic analysis of the 3-gene dataset by parsimony and Bayesian methods yielded well-resolved and highly similar trees. Bayesian posterior probabilities were high for most clades; bootstrap (parsimony) support was reliable for approximately half of the internal nodes. Parsimony analysis of the morphological dataset yielded a poorly resolved tree, whereas total evidence analysis (molecular plus morphological data) yielded a phylogenetic hypothesis consistent with, but less resolved than trees inferred from combined molecular data. Parsimony mapping of morphological characters on the 3-gene trees showed that most structural features of steinernematids are highly homoplastic. The distribution of nematode foraging strategies on these trees predicts that S. hermaphroditum, S. diaprepesi and S. longicaudum (US isolate) have cruise forager behaviours.     

LSU rDNA-BASED RFLP ASSAYS FOR DISCRIMINATING SPECIES AND STRAINS OF ALEXANDRIUM (DINOPHYCEAE)1

In a previous study large-subunit ribosomal RNA gene (LSU rDNA) sequences from the marine dinoflagellates Alexandrium tamarense (Lebour) Balech, A. catenella (Whedon et Kofoid) Balech, A. fundyense Balech, A. affine (Fukuyo et Inoue) Balech, A. minutum Halim, A. lusitanicum Balech, and A. andersoni Balech were compared to assess inter- and intraspecific relationships. Many cultures compared in that study contained more than one class of LSU rDNA. Sequencing pooled clones of rDNA from single cultures revealed length heterogeneities and sequence ambiguities. This complicated sequence comparisons because multiple rDNA clones from a single culture had to be sequenced individually to document the different classes of molecules present in that culture. A further complication remained as to whether or not the observed intraculture sequence variations were reliable genetic markers or were instead artifacts of the polymerase chain reaction (PCR) amplification, cloning, and/or sequencing methods employed. The goals of the present study were to test the accuracy of Alexandrium LSU rDNA sequences using restriction fragment-length polymorphism (RFLP) analysis and to devise RFLP-based assays for discriminating among representatives of that group. Computer-assisted examination of the sequences allowed us to identify a set of restriction enzymes that were predicted to reveal species, strain, and intraculture LSU rDNA heterogeneities. All groups identified by sequencing were revealed independently and repeatedly by RFLP analysis of PCR-amplified material. Five ambiguities and one length heterogeneity, each of which ascribes a unique group of Alexandrium species or strains, were confirmed by restriction digests. Observed intraculture LSU rDNA heterogeneities were not artifacts of cloning and sequencing but were instead a good representation of the spectrum of molecules amplified during PCR reactions. Intraculture LSU rDNA heterogeneities thus serve as unique genetic markers for particular strains of Alexandrium, particularly those of A. tamarense, A. catenella, and A. fundyense. However, some of these “signature heterogeneities” represented a smaller portion of PCR product than was expected given acquired sequences. Other deviations from predicted RFLP patterns included incomplete digestions and appearance of spurious products. These observations indicate that the diversity of sequences in PCR product pools were greater than that observed by cloning and sequencing. The RFLP tests described here are useful tools for characterizing Alexandrium LSU rDNA to define the evolutionary lineage of cultures and are applicable at a fraction of the time, cost, and labor required for sequencing.     

Phylogenetic relationship withinFusarium sambucinum Fuckel sensu lato,determined from ribosomal RNA sequences

Partial ribosomal RNA nucleotide sequences were determined for 11 strains ofFusarium sambucinum Fuckelsensu lato to assess by molecular genetic means, Nirenberg's recent morphotaxonomic interpretation which split the species into three distinct taxa:F. sambucinum sensu stricto, F. torulosum, and one other species, as yet unnamed (Fusarium species nova). Four sequence patterns were identified among the 11 strains. Two sequences that varied at one site were found among strains ofF. sambucinum, strains ofF. torulosum andFusarium sp. nov. showed no intraspecific variation. Interspecific comparisons revealed nucleotide sequence differences of 3–9 substitutions in the ca. 240 nucleotide rRNA segment examined. Although interspecific differences are not large in terms of percent nucleotide substitution, they are much larger than the observed intraspecific variation and support the morphological interpretation distinguishing three taxa. When the data were analysed using parsimony and bootstrapping, the three taxon tree was well supported. The phylogenetic arrangement of these strains is congruent with secondary metabolite profile similarities.     

Molecular Analysis of Ribosomal RNA Gene of Red Tide Algae Obtained from the Seto Inland Sea

Eleven clones from five species of the planktonic microalgae, (Chattonella antiqua, Chattonella marina, Heterosigma akashiwo, Alexandrium catenella, and Scrippsiella trochoidea), which were collected from the Seto Inland Sea in Japan and from Thailand, were subjected to nucleotide sequence analysis of the D1/D2 domain of the large subunit (LSU) of their ribosomal RNA genes. After amplification by polymerase chain reaction using degenerated primers, whole-nucleotide sequences for the D1/D2 domains of the LSU rRNA gene of 11 microalgae were analyzed. Phylogenic tree analysis using these nucleotide sequences showed each species located in a cluster corresponding to its morphological classification. The nucleotide sequence data for Chattonella spp. suggest that multiple clones of both Chattonella antiqua and Chattonella marina are present in the Seto Inland Sea and that red tide blooms of Chattonella spp. in different years may have contained different clones. Received September 6, 1999; accepted December 16, 1999.     

Studies on the diversity of the distinct phylogenetic lineage encompassing Glomus claroideum and Glomus etunicatum

Morphological and molecular characters were analysed to investigate diversity within isolates of the Glomus claroideum/Glomus etunicatum species group in the genus Glomus. The inter- and intra-isolate sequence diversity of the large subunit (LSU) rRNA gene D2 region of eight isolates of G. claroideum and G. etunicatum was studied using PCR-single strand conformational polymorphism (SSCP)-sequencing. In addition, two isolates recently obtained from Southern China were included in the analysis to allow for a wider geographic screening. Single spore DNA isolation confirmed the magnitude of gene diversity found in multispore DNA extractions. An apparent overlap of spore morphological characters was found between G. claroideum and G. etunicatum in some isolates. Analysis of the sequence frequencies in all G. etunicatum and G. claroideum isolates (ten) showed that four LSU D2 sequences, representing 32.1% of the clones analysed for multispore extraction (564) were found to be common to both species, and those sequences were the most abundant in four of the ten isolates analysed. The frequency of these sequences ranged between 23.2% and 87.5% of the clones analysed in each isolate. The implications for the use of phenotypic characters to define species in arbuscular mycorrhizal fungi are discussed. The current position of G. claroideum/G.etunicatum in the taxonomy of the Glomeromycota is also discussed.     

Comparison of primary and secondary 26S rRNA structures in twoTetrahymena species: Evidence for a strong evolutionary and structural constraint in expansion segments

Summary We have determined the nucleotide sequence of the 26S large subunit (LSU) rRNA genes for twoTetrahymena species,T. thermophila andT. pyriformis. The inferred rRNA sequences are presented in their most probable secondary structures based on compensatory mutations, energy, and conservation criteria. The majority of the nucleotide changes between the twoTetrahymena LSU rRNAs and the positions of a relatively large deletion and of the processing cleavage sites resulting in the generation of the hidden break are all located within the so-called divergent domains or expansion segments. These are regions within the common core of secondary structure where expansions have taken place during the evolution of the rRNA of higher eukaryotes.The dispensable nature of some of the expansion segments has been taken as evidence of their non-functionality. However, our data show that a considerable selective constraint has operated to presesrve the secondary structure of these segments. Especially in the case of the D2 and D8 segments, the presence of a considerable number of compensatory base changes suggests that the secondary structure of these regions is of functional importance. Alternatively, these expansion segments may have maintained characteristic folding patterns because only such structures are being tolerated within otherwise functionally important regions.     

Phylogenetic analysis of Glomeromycota by partial LSU rDNA sequences

We analyzed the large subunit ribosomal RNA (rRNA) gene [LSU ribosomal DNA (rDNA)] as a phylogenetic marker for arbuscular mycorrhizal (AM) fungal taxonomy. Partial LSU rDNA sequences were obtained from ten AM fungal isolates, comprising seven species, with two new primers designed for Glomeromycota LSU rDNA. The sequences, together with 58 sequences available from the databases, represented 31 AM fungal species. Neighbor joining and parsimony analyses were performed with the aim of evaluating the potential of the LSU rDNA for phylogenetic resolution. The resulting trees indicated that Archaeosporaceae are a basal group in Glomeromycota, Acaulosporaceae and Gigasporaceae belong to the same clade, while Glomeraceae are polyphyletic. The results support data obtained with the small subunit (SSU) rRNA gene, demonstrating that the LSU rRNA gene is a useful molecular marker for clarifying taxonomic and phylogenetic relationships in Glomeromycota.     

Polynucleotide sequence relatedness among motileAeromonas species

Relatedness among 55 motileAeromonas strains was assessed by determining the extent of reassociation in heterologous DNA preparations. S1 nuclease and diethylaminoethyl-cellulose filters were used to separate reassociated nucleotide sequences from nonreassociated sequences and to determine the thermal stability of related nucleotide sequences. The motile aeromonads would presently consist of three species:A. hydrophila (type strain ATCC 7966),A. caviae (type strain ATCC 15468), andA. sobria (type strain CIP 7433). These three species are clearly separated on the basis of both biochemical characteristics and similarities in DNA. Each of these three species contains more than one hybridization group: three groups inA. hydrophila; two groups inA. caviae; and at least two groups inA. sobria. DNA hybridization groups are biochemically similar within each species. When additional data is available, these separate DNA hybridization groups may merit designation. Any decision to delineate species in terms of DNA hybridization groups must await a phenotypic basis for their separation.     

Phylogenetic study of Cryptococcus laurentii mycocinogenic strains

Previous studies reported Cryptococcus laurentii strains to secrete mycocins of different inhibition range. These mycocinogenic strains were studied by molecular methods to verify their phylogenetic relationships. Based on the combined analysis of LSU and ITS sequence data, these strains were confirmed to belong to the species C. laurentii. Therefore, different strains of the same basidiomycetous yeast C. laurentii secrete a distinct mycocin. This trait has not been demonstrated for basidiomycetous yeasts before. Additionally, this study provides several DNA-barcodes (LSU, ITS and RPB1) and reports their variability for this species.     

Torulaspora indica a novel yeast species isolated from coal mine soils

Four yeast strains (APSS 805, APSS 806, APSS 815 and AP-18) belonging to a novel Torulaspora species were isolated from coal mine soils of Singareni in Andhra Pradesh state, India. Another strain (PBA-22) was isolated from agricultural field soil from Gujarat state, India. The vegetative cells of all these strains were round, haploid and produced asci by conjugation between independent cells or mother cell and bud, with rough ascospores, suggesting their possible relation to ascomycetous yeast genus Torulaspora. Phylogenetic analysis of the D1/D2 domain of the large subunit (LSU) rRNA gene and Internal Transcribed Spacer (ITS) regions revealed that, among the five strains, three viz. APSS 805, APSS 806 and APSS 815 have identical sequences. The other two strains (AP-18 and PBA-22) differed from the other three strains in less than 1% nucleotide substitutions in the combined D1/D2 domain and ITS sequences, indicating that all of them (five strains) may belong to the same species. These five strains were closely related to Torulaspora globosa, but showed more than 3–7% sequence divergence from T. globosa and all other species in the genus Torulaspora in the combined sequence analysis of D1/D2 domain and ITS region of rRNA gene. In addition, these strains also showed distinct microsatellite finger-printing pattern from related species and differed in several physiological responses suggesting that these strains belong to a novel species of Torulaspora. We propose to name these strains as Torulaspora indica sp. nov., and designate APSS 805^T = MTCC 9772 ^T = CBS 12408 ^T as the type strain of this novel species. The Mycobank number of the novel species is MB 563738.     

Genetic diversity of Paramecium species on the basis of multiple loci analysis and ITS secondary structure models

Among ciliates, Paramecium has become a privileged model for the study of “species problem” particularly in the case of the “Paramecium aurelia complex” that has been intensely investigated. Despite extensive studies, the taxonomy of Paramecium is still challenging. The major problem is an uneven sampling of Paramecium with relatively few representatives of each species. To investigate species from the less discovered region (Pakistan), 10 isolates of Paramecium species including a standing-alone FT8 strain previously isolated by some of us were subjected to molecular characterization. Fragments of 18S recombinant DNA (rDNA), ITS1-5.8S-ITS2-5′LSU rDNA, cytochrome c oxidase subunit II, and hsp70 genes were used as molecular markers for phylogenetic analysis of particular isolates. The nucleotide sequences of polymerase chain reaction products of all markers were compared with the available sequences of relevant markers of other Paramecium species from GenBank. Phylogenetic trees based on all molecular markers showed that all the nine strains had a very close relationship with Paramecium primaurelia except for the FT8 strain. FT8 consistently showed its unique position in comparison to all other species in the phylogenetic trees. Available sequences of internal transcribed spacer 1 (ITS1) and ITS2 and some other ciliate sequences from GenBank were used for the construction of secondary models. Two highly conserved helices supported by compensatory base changes among all ciliates of ITS2 secondary structures were found similar to other eukaryotes. Therefore, the most conserved 120 to 180 base pairs regions were identified for their comparative studies. We found that out of the three helices in ITS1 structure, helix B was more conserved in Paramecium species. Despite various substitutions in the primary sequence, it was observed that secondary structures of ITS1 and ITS2 could be helpful in interpreting the phylogenetic relationships both at species as well as at generic level.     

Relationships among didelphid marsupials based on sequence variation in the mitochondrial cytochrome B gene

Variation in the mitochondrial cytochrome b gene (nucleotide and amino acid sequences) is evaluated for 9 genera and 15 species of American opossums in the family Didelphidae, using the American caenolestid rat opossumLestoros and the New Guinean peroryctid bandicootEchimypera as outgroups. Phylogenetic analyses (parsimony and distance) strongly support the monophyly of the Didelphidae and delineate two major clades; (1)Didelphis andPhilander are strongly aligned sister taxa, withMetachirus weakly but consistently associated with them, and (2)Marmosa plusMicoureus, withMonodelphis falling outside that pair. The generaMarmosops, Caluromys, andGlironia exhibit varied relationships, depending upon the method of analysis and data (DNA or amino acid sequences) used, but generally are placed individually or in combinations near or at the base of the didelphid radiation. Some aspects of these relationships are consistent with current taxonomic views, but others are in marked contrast. Specifically, a clade comprised of the mouse opossumsMarmosa, Micoureus, andMarmosops is strongly rejected by log-likelihood analysis, contrary to expectations from some current classifications. Also, the woolly opossumsCaluromys andGlironia also do not form a sister-taxon relationship, as suggested by their placement in a subfamily separate from the remaining didelphids examined. However, such a relationship cannot be rejected from log-likelihood analyses. The relationships suggested fromcyt-b sequences are strongly concordant with those based on DNA-DNA hybridization analyses. In addition to systematic and phylogenetic properties, molecular evolution of the didelphid cytochrome b gene sequence is characterized according to nucleotide bias and rate differentials at each codon position and across the entire sequence.To whom correspondence should be addressed.     

Utility of mitochondrial‐encoded cytochrome c oxidase I gene for phylogenetic analysis and species identification of the planktonic diatom genus Skeletonema

Small subunit (SSU) and large subunit (LSU) rDNA sequences have been commonly used to delineate the taxonomy and biogeography of the planktonic diatom genus Skeletonema, but the genes occur as multiple copies and are therefore not suitable for barcoding purposes. Here, we analyzed phylogenetic relationships of Skeletonema using the mitochondrial‐encoded cytochrome c oxidase I gene (cox1), as well as partial LSU rDNA (D1–D3) and SSU rDNA, to identify the factors that define species and to evaluate the utility of these three markers for this taxon. Twelve Skeletonema species were divided into six clades, I–VI, each of which comprised the same species by the three markers: clades I (S. japonicum, S. grethae, S. pseudocostatum, and S. tropicum), II (S. menzelii), III (S. dohrnii and S. marinoi), IV (S. costatum, S. potamos, and S. subsalsum), V (S. grevillei), and VI (S. ardens). However, the branching order among these clades was incongruent among the markers. In clade III, six S. marinoi strains had identical cox1 sequences. These S. marinoi strains branched along with S. dohrnii, except for strains from the Gulf of Naples, with high support in cox1. Species delimitation between S. dohrnii and S. marinoi was therefore not supported. In clade IV, S. costatum and S. subsalsum were robustly clustered, with S. potamos as a sister clade in the cox1 tree, not in the LSU and SSU trees. In clade II, cox1 also confirmed that S. menzelii includes three subclades potentially distinguishable from each other by morphological features. Cox1 proved to be the most useful marker for the identification of Skeletonema species because it gave a tree with highly supported clades, has sufficient variation within and among species, encodes a protein in a single copy, and requires relatively few primers.