IJSEM Track the topics, authors and articles important to you
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplementary Material
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via CrossRef
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Cadez, N.
Right arrow Articles by Smith, M. Th.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Cadez, N.
Right arrow Articles by Smith, M. Th.
Agricola
Right arrow Articles by Cadez, N.
Right arrow Articles by Smith, M. Th.
Int J Syst Evol Microbiol 56 (2006), 1157-1165; DOI  10.1099/ijs.0.64052-0
© 2006 International Union of Microbiological Societies

Phylogenetic placement of Hanseniaspora–Kloeckera species using multigene sequence analysis with taxonomic implications: descriptions of Hanseniaspora pseudoguilliermondii sp. nov. and Hanseniaspora occidentalis var. citrica var. nov.

Neza Cadez1, Peter Raspor1 and Maudy Th. Smith2

1 University of Ljubljana, Biotechnical Faculty, Department of Food Science and Technology, Jamnikarjeva 101, 1000 Ljubljana, Slovenia
2 Centraalbureau voor Schimmelcultures, PO Box 85167, 3508 AD Utrecht, The Netherlands

Correspondence
Neza Cadez
neza.cadez{at}bf.uni-lj.si


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
Two protein-coding genes, actin and translation elongation factor-1{alpha} (EF-1{alpha}), as well as two ribosomal gene regions, D1/D2 domains of the large subunit and both internal transcribed spacers including the 5.8S gene region, were evaluated regarding their usefulness for reconstruction of phylogenetic relationships in the Hanseniaspora–Kloeckera species group. This included analyses of sequence divergence values, heterogeneity of evolutionary rates and the reliability of the inferred trees. Both protein-coding genes showed greater capacities to resolve at the strain level and between the closely related species of Hanseniaspora–Kloeckera, compared with the ribosomal gene regions. However, to obtain a fully resolved and reliable phylogenetic tree that reflected the biological relationships it was necessary to combine three congruent sequence datasets. The novel species Hanseniaspora pseudoguilliermondii sp. nov. (type strain CBS 8772T) is described as a result of the application of various molecular approaches to delimit species. Furthermore, incongruent gene genealogies of genetically divergent strains of Hanseniaspora occidentalis, as determined by amplified fragment length polymorphism analysis and DNA–DNA reassociation measurements, indicated the presence of two novel varieties, H. occidentalis var. occidentalis (type strain CBS 2592T) and H. occidentalis var. citrica var. nov. (type strain CBS 6783T), which could be distinguished by habitat preference.


Abbreviations: AFLP, amplified fragment length polymorphism; EF-1{alpha}, translation elongation factor-1{alpha}; ITS, internal transcribed spacer; LSU, large subunit; RAPD, randomly amplified polymorphic DNA

The GenBank/EMBL/DDBJ accession numbers for the gene sequences determined in this study are listed in Table 1Go.

Figures showing saturation plots and phylogenetic analyses of sequences of actin and EF-1{alpha} genes and the D1/D2 regions of the LSU rDNA and ITS regions of Hanseniaspora–Kloeckera strains, and tables comparing dataset characteristics and characteristics inferred from the datasets using maximum-parsimony analyses and DNA–DNA relatedness values among strains of Hanseniaspora occidentalis are available as supplementary material in IJSEM Online.


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
The name Hanseniaspora was proposed for ascosporogenous apiculate yeasts by Zikes (1911)Go, although their common presence in fruit juices and fermenting musts had already led to the description of the apiculate species Saccharomyces apiculatus (Reess, 1870Go). The unstable ability of these yeasts to form ascospores contributed to several changes in nomenclature and definition of the genus (Lodder & Kreger-van Rij, 1952Go). With the introduction of molecular taxonomy, relationships between the species of Hanseniaspora and of the anamorphic genus Kloeckera were established on the basis of DNA base composition (Nakase & Komagata, 1970Go) and DNA–DNA relatedness (Meyer et al., 1978Go). The latter study represents the basis for the current classification, which was confirmed by several phylogenetic studies based on different parts of the ribosomal gene and two protein-coding genes (Yamada et al., 1992aGo; Boekhout et al., 1994Go; Esteve-Zarzoso et al., 2001Go; Kurtzman & Robnett, 2003Go). The genus Hanseniaspora currently comprises ten species, four of which (Hanseniaspora meyeri, Hanseniaspora clermontiae, Hanseniaspora lachancei and Hanseniaspora opuntiae) were recently circumscribed on the basis of DNA–DNA reassociation data (Smith, 1998Go; Cadez et al., 2003Go).

Phylogenetic placement of members of Hanseniaspora within the Saccharomyces clade revealed that the genus is divided into two clusters (Kurtzman & Robnett, 1998Go, 2003Go). The first cluster contains the species Hanseniaspora vineae, Hanseniaspora osmophila and Hanseniaspora occidentalis. Strains belonging to the latter species were found to be heterogeneous with respect to their randomly amplified polymorphic DNA (RAPD) profiles, electrophoretic karyotypes and restriction patterns of the internal transcribed spacer (ITS) regions (Cadez et al., 2002Go). The second cluster comprises Hanseniaspora valbyensis and four recently described species that are closely related to Hanseniaspora uvarum and Hanseniaspora guilliermondii (Cadez et al., 2003Go). The poorly resolved rDNA-based phylogenetic tree suggested that there were very close relationships among species of the H. uvarum–H. guilliermondii complex, which was not in agreement with the genetic relatedness determined by DNA–DNA reassociation analyses (Cadez et al., 2003Go). This supports the observation that the relationships between recently diverged species complexes are difficult to resolve on the basis of a single gene only (Rokas et al., 2003Go).

Several studies have demonstrated the phylogenetic usefulness of protein-coding genes for inferring relationships among different taxonomic units (Baldauf et al., 2000Go; Daniel et al., 2001Go; Belloch et al., 2000Go; O'Donnell et al., 1998Go). Housekeeping genes such as translation elongation factor-1{alpha} (EF-1{alpha}) and actin are frequently used as independently evolving markers because they do not appear to have undergone horizontal transfer or gene duplication. Moreover, the rates of sequence substitutions in exons are constrained by codon positions, which may determine the taxonomic level of analysis (Baldauf & Palmer, 1993Go; Kretzer & Bruns, 1999Go).

Because the species of Hanseniaspora–Kloeckera fall into several groups of closely related taxa separated by large phylogenetic distances (Cadez et al., 2003Go), the phylogenetic marker for resolving Hanseniaspora–Kloeckera species should have both the capacity to resolve the sister species and the potential to realistically reconstruct the phylogeny between the divergent lineages. Consequently, the aim of the study was to assess the phylogenetic usefulness of four phylogenetic markers by describing sequence divergence and by comparing the inter- and intraspecific phylogenetic relationships with those derived from genomic DNA–DNA relatedness and physiological properties.


    METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
Yeast strains.
The strains examined and their origins and accession numbers are listed in Table 1Go.


View this table:
[in this window]
[in a new window]
 
Table 1. Strains studied, their origin and GenBank/EMBL/DDBJ accession numbers

CBS, Centraalbureau voor Schimmelcultures, Utrecht, The Netherlands; DBVPG, Industrial Yeast Collection, Perugia, Italy; NCAIM, National Collection of Agricultural and Industrial Microorganisms, Budapest, Hungary.

 
Strains of H. occidentalis were characterized physiologically by using standard methods (Yarrow, 1998Go). Fermentation ability was determined after 3 and 7 days of incubation. Utilization of carbon compounds was tested in liquid media in tubes, with reciprocal shaking at 30 r.p.m. at 25 °C for 3 weeks. Assimilation of nitrogen compounds was examined by using the auxoanographic method after 1 week (Barnett et al., 2000Go).

DNA extraction, amplification and sequencing.
DNA was isolated from cultures grown on yeast-malt agar (0.3 % yeast extract, 0.3 % malt extract, 0.5 % peptone, 1 % glucose, 2 % agar), according to the method of Möller et al. (1992)Go.

PCR primers for the amplification of actin and EF-1{alpha} protein-coding regions were designed based on sequences of Saccharomyces cerevisiae (GenBank accession numbers NC_001138 and U51033, respectively). Additional primers were designed from generated sequences of the actin gene of H. vineae (Act-HV) and the EF-1{alpha} gene of H. meyeri (EF1-HM). Primer sequences, annealing temperatures and annealing sites and references are given in Table 2Go.


View this table:
[in this window]
[in a new window]
 
Table 2. Primers used

Tanneal, Annealing temperature. Position refers to S. cerevisiae sequence.

 
Amplifications were performed under the following conditions: an initial denaturing step of 5 min at 94 °C, followed by 35 cycles of 30 s at 94 °C, 30 s at primer-pair-specific annealing temperatures (Table 2Go), 2 min at 72 °C and termination with a final extension step of 7 min at 72 °C. The PCR products were purified using GenElute PCR Clean-up (Sigma–Aldrich), according to the manufacturer's instructions. The pGEM-T vector system was used to clone the PCR products of the protein-coding genes, following the procedure described in the manual (Promega). Plasmid DNA was isolated and purified for subsequent sequencing by using GenElute Plasmid Minipreps (Sigma–Aldrich). Nucleotide sequences were determined in both directions by using a Thermo Sequenase Primer Cycle sequencing kit (Amersham Biosciences) and labelled universal primers M13(–10) or M14R with an ALFexpress II Automatic DNA Sequencer (Amersham Pharmacia Biotech). The D1/D2 regions of the large subunit (LSU) rDNA and ITS regions were sequenced as described previously (Cadez et al., 2003Go). Sequences of both strands were aligned and edited with the ALFwin Sequence Analyser 2.10.06 software (Amersham Pharmacia Biotech).

Phylogenetic analyses.
The number of variable and parsimony-informative sites, the genomic DNA G+C contents and the mean uncorrected sequence divergences between the groups of taxa were calculated using the computer program MEGA 2.1 (Kumar et al., 2001Go).

Sequences were aligned automatically using the multiple sequence alignment program CLUSTAL X version 1.83 (Thompson et al., 1997Go). Positions where gaps existed in one of the sequences were manually excluded using BioEdit (Hall, 1999Go). All subsequent phylogenetic analyses employed the PAUP* 4.0b10 software package (Swofford, 2002Go). Most-parsimonious trees were generated by a heuristic search procedure with 1000 random addition replicates and tree bisection–reconnection branch swapping. Nucleotide sites were equally weighted. The stability of the branches was assessed by bootstrap analysis (Felsenstein, 1985Go) in which 1000 replicates were used. The homogeneity of the signal from the four datasets was tested with partition-homogeneity tests (Farris et al., 1995Go) implemented in PAUP* 4.0b10 with 1000 heuristic replicate searches and exclusion of the invariant characters.

For each dataset, a saturation plot was drawn to evaluate the degree of homoplasy and the heterogeneity of evolutionary rates and to choose an optimal outgroup species (Bush & Everett, 2001Go; Misof et al., 2001Go). The corrected evolutionary rates of nucleotide substitutions per site were plotted against the proportion of nucleotide changes per site (p distance). The corrected distances were calculated under the optimality criteria of maximum-likelihood using PAUP* 4.0b10 with corrections for multiple substitutions at the same site, substitutional rate biases and differences in evolutionary rates among sites for each dataset. These corrections were included in substitution models selected using the likelihood-ratio test implemented in MODELTEST 3.5 (Posada & Crandall, 1998Go). The general time reversible model with a gamma distribution parameter for among-site rate variations (GTR+{Gamma}) was found to be of adequate complexity for the actin gene dataset. The Tamura–Nei model was chosen for the EF-1{alpha} and D1/D2 datasets, with only the gamma distribution parameter (TN+{Gamma}) for the former dataset and the addition of the estimate of invariant sites (TN+I+{Gamma}) for the latter. The Hasegawa–Kishino–Yano 85 model with the gamma distribution parameter and the estimate of invariant sites (HKY+I+{Gamma}) was proposed for the ITS dataset. When sequences of closely related taxa are compared, a strong correlation between phyletic and p distances exists. With increasing distances between taxa, the degree of sequence divergence is underestimated due to parallel and backward substitutions and the regression curve starts to decline. When the slope of the curve approaches zero, mutational saturation is reached and the phylogenetic reconstruction is more likely to suffer from artefacts such as long-branch attraction (Felsenstein, 1978Go).

S. cerevisiae, Lachancea kluyveri and Pichia anomala were tested as outgroup species using the saturation plots. Because the levels of sequence divergence between S. cerevisiae, L. kluyveri and ingroup species overlapped, P. anomala was chosen as the outgroup.

AFLP assay.
Amplified fragment length polymorphism (AFLP) analysis was performed according to the procedure described by Vos et al. (1995)Go and adapted for fluorescence-based detection by Bandelj et al. (2004)Go. Genomic DNA was digested with the restriction enzymes MseI (1 U) and EcoRI (5 U) and linked to MseI and EcoRI adapters (50 and 5 pmol, respectively). Restricted and ligated DNA was pre-amplified using MseI-core and EcoRI-core primers. The PCR conditions for the pre-amplification step were 2 min at 72 °C, followed by 20 cycles of 30 s at 94 °C, 30 s at 56 °C and 2 min at 72 °C. Selective amplifications were performed with diluted (1 : 10, by vol.) pre-amplification mixture using three selective primer combinations (MseC/EcoAC, MseC/EcoCG and MseAC/EcoCG). The EcoAC and EcoCG primers were labelled for fluorescence detection with Cy5 dye. The PCR conditions were 2 min at 94 °C, 10 cycles of 30 s at 94 °C, 30 s at 66 °C (decreasing by 1 °C each step of the cycle) and 2 min at 72 °C, followed by 20 cycles of 30 s at 94 °C, 30 s at 56 °C and 2 min at 72 °C. The PCR products were prepared for electrophoresis by addition of an equal volume of formamide loading dye (98 % formamide, 10 mg blue dextran ml–1, 50 mM EDTA) and denatured by heating at 94 °C for 3 min. AFLP fragments were separated on 7.5 % polyacrylamide gels containing 7 M urea on ALF Express II (Amersham Pharmacia Biotech). ALF-Express sizer 50-500 (Amersham Pharmacia Biotech) was used as an external standard.

The AFLP profiles were combined in composite fingerprints using the Bionumerics software package (version 4.0; Applied Maths). Similarities between profiles were calculated using Pearson's correlation coefficient and cluster analysis was performed using the UPGMA algorithm.

Genomic DNA analysis.
For DNA extraction, the strains were grown for 2 days at 25 °C on a rotary shaker at 125 r.p.m. in 0.5 l yeast extract/malt extract broth (YM; Difco) using 1 l flat-bottom flasks. Isolation and purification of the DNA, determination of DNA base composition and DNA–DNA reassociation were performed according to procedures described previously (Golubev et al., 1989Go). Strains used for DNA–DNA reassociation are indicated in Table 1Go.


    RESULTS AND DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
Comparison of phylogenetic information in different datasets for resolving Hanseniaspora–Kloeckera species
Partial sequences of two nuclear protein-coding genes, actin and EF-1{alpha}, were generated and analysed phylogenetically for 26 strains of Hanseniaspora–Kloeckera species. The newly generated sequences were compared with the D1/D2 sequences of the LSU rDNA and ITS1–5.8S rDNA–ITS2 sequences. Sequences for the two regions of the ribosomal gene complex were obtained from previous studies in which they were used as phylogenetic markers for inferring relationships between Hanseniaspora–Kloeckera species (Kurtzman & Robnett, 1998Go; Cadez et al., 2003Go). We compared sequence divergence values, the heterogeneity of evolutionary rates and the reliability of the inferred trees from the four datasets of Hanseniaspora–Kloeckera strains. The results are presented in Supplementary Table S1 and Fig. S1 in IJSEM Online. The four datasets are influenced by functional and structural constraints and therefore differed in the amount and form of variation. For example, the ITS regions are under more relaxed structural constraints, which was evident from the highest mean sequence divergence between Hanseniaspora–Kloeckera species (0.253; Supplementary Table S1). On the other hand, the protein-coding datasets resolved better at the strain level (0.005; Supplementary Table S1) and between the closely related species of Hanseniaspora–Kloeckera (0.03; Supplementary Table S1). These results are consistent with the observations of Daniel & Meyer (2003)Go that synonymous substitutions in codons account for the high degree of sequence divergence between closely related taxa. Furthermore, the saturation plots (Supplementary Fig. S1 in IJSEM Online) indicated heterogeneity in the evolutionary rates for the ITS dataset (i.e. mean p distance, 0.12; SD=0.11; Supplementary Fig. S1, ITS, point 3). This had been already predicted from the fact that the number of substitutions in the ITS regions was incongruent with the DNA–DNA reassociation values for closely related species of Hanseniaspora (Cadez et al., 2003Go).

For all four datasets, we detected an overlap in the levels of sequence divergence between the distantly related Hanseniaspora–Kloeckera species and the levels of sequence divergence between the Hanseniaspora–Kloeckera species and the outgroup species of the Saccharomycetaceae (Supplementary Fig. S1 in IJSEM Online). Kurtzman (2003)Go observed a similar overlap of intra- and intergeneric distances for Hanseniaspora and four related genera calculated from a dataset of six genes. In view of these results, the reinstatement of the genus Kloeckeraspora to accommodate H. occidentalis, H. osmophila and H. vineae, as proposed by Yamada et al. (1992b)Go, could be justified. Nevertheless, we advocate the argument of Boekhout et al. (1994)Go against splitting the genus on the basis of genetic divergence only, because the two species groups share many similarities in morphology, physiology and ecology.

The phylogenetic usefulness of each dataset could only be determined after cladistic analyses (Damgaard & Cognato, 2003Go). Most-parsimonious trees based on single datasets and their characteristics are available as Supplementary Fig. S2 and in Supplementary Table S1, respectively, in IJSEM Online. In general, the trees showed some degree of congruence. However, the lack of divergence between the recently diverged H. uvarum–H. meyeri–H. clermontiae and H. guilliermondii–H. lachancei–H. opuntiae–Hanseniaspora sp. CBS 8772T species complexes resulted in a lack of statistical support for these clades in single genealogies. Therefore, we concatenated the datasets in order to increase the phylogenetic information. Prior to concatenation, we tested the congruence of the datasets with the partition homogeneity test (Farris et al., 1995Go). When all four datasets were combined, significant incongruence (P=0.016) was found and this resulted in a weakly supported phylogenetic tree. Therefore we tested different combinations of the concatenated datasets that resulted in the highest P values for the combined dataset without either ITS or EF-1{alpha} sequence data (P=0.35 and 0.24, respectively). Finally, parsimony analysis of the combined dataset of actin, D1/D2 and ITS sequences without EF-1{alpha} data produced a completely resolved tree with high statistical support (Fig. 1Go).


Figure 1
View larger version (13K):
[in this window]
[in a new window]
 
Fig. 1. Phylogeny of Hanseniaspora–Kloeckera species based on concatenated datasets of actin, D1/D2 and ITS gene sequences. Each species is represented by its type strain. The single most-parsimonious tree (tree length, 1169; consistency index=0.7707; retention index=0.8319) was constructed by heuristic search procedure under the parsimony criterion in PAUP*. Bootstrap percentages from 1000 replicates are shown. Pichia anomala was used as the outgroup. Bar, number of nucleotide substitutions. Interspecies DNA–DNA reassociation values are indicated by arrows. +/– show responses to key physiological tests for Hanseniaspora–Kloeckera species identification: 1, fermentation of sucrose; 2, growth with 0.01 % cycloheximide; 3, assimilation of 2-keto-D-gluconate; 4, growth at 37 °C; 5, growth at 30 °C.

 
With the aim of evaluating whether the inferred phylogenetic tree reflects the biological relationships, we compared the relationships with key physiological characteristics and DNA–DNA reassociation data of Hanseniaspora–Kloeckera species reported in previous studies (Meyer et al., 1978Go; Cadez et al., 2003Go). Fig. 1Go shows that clustering of the species is consistent with the current taxonomy. Clusters of closely related species (H. vineae–H. osmophila, H. guilliermondii–H. lachancei–H. opuntiae–Hanseniaspora sp. CBS 8772T, H. uvarum–H. meyeri–H. clermontiae and H. valbyensis–Kloeckera lindneri) generally share the same nutritional characteristics with two exceptions, H. lachancei and H. clermontiae. Furthermore, the congruence between the relationships determined from the concatenated datasets and DNA–DNA reassociation improved, in particular for H. lachancei, H. opuntiae and Hanseniaspora sp. CBS 8772T, which were placed inconsistently in the single genealogies (bootstrap support <50 %; Supplementary Fig. S2 in IJSEM Online). The latter species was not named in our previous study (Cadez et al., 2003Go) because its description would rest on a single isolate. However, based on additional molecular evidence of species delimitation we propose that strain CBS 8772T represents a novel species with the name Hanseniaspora pseudoguilliermondii sp. nov.

Conflicting gene genealogies within H. occidentalis
A comparison of multiple gene genealogies can be used to detect reproductive or geographical isolation between distinct populations. Incongruences between different gene genealogies may be indicative of genetic exchange among individuals within a species (Taylor et al., 2000Go). In order to clarify the species concept for H. occidentalis, we examined the consistency of phylogenetic relationships among H. occidentalis strains using three phylogenetic markers. In a previous study (Cadez et al., 2002Go), the six strains of H. occidentalis were found to be genetically heterogeneous, as shown by their chromosomal make-up, RAPD-PCR profiles and restriction patterns of the ITS regions.

First, we determined the genetic variability of 14 H. occidentalis strains by cluster analysis of the combined AFLP fingerprints. Fig. 2Go shows that the strains segregated into two major groups at a similarity level of 32 %. Interestingly, this grouping coincided with the source of isolation: the strains from group I were isolated from soil whereas the strains from group II were isolated from oranges or orange juices. An exception was strain CBS 2569, which was isolated from a fruit fly. Nevertheless, this strain showed a moderate AFLP similarity (42 %) to the rest of the strains in group I.


Figure 2
View larger version (20K):
[in this window]
[in a new window]
 
Fig. 2. UPGMA dendrogram of H. occidentalis strains based on combined AFLP fingerprints obtained with three primer combinations (MseC/EcoAC, MseC/EcoCG and MseAC/EcoCG). The distance between the strains was calculated as the Pearson's correlation coefficient.

 
The existence of two major groups was confirmed by DNA–DNA reassociation (Supplementary Table S2 in IJSEM Online). The reassociation values obtained when group I strains (CBS 2592T, CBS 280, CBS 282, CBS 283, CBS 284, CBS 2569 and CBS 2578) were paired with group II strains (CBS 6783T and CBS 6782) ranged from 60 to 78 %. The DNA–DNA relatedness values within group I, which included the type strain, ranged from 88 to 100 %. The mutual DNA–DNA relatedness of strains of group II was 99 %. Within group I, the lowest relatedness value (88 %) was observed between the type strain and strains CBS 2569 and CBS 282, which is in agreement with the AFLP fingerprinting results (Fig. 2Go). However, the interpretation of DNA–DNA reassociation values of 60–78 % determined between the major groups is ambiguous without further evidence of genetic isolation between the strains.

Phylogenetic trees inferred from the sequences of the actin gene, the EF-1{alpha} gene and the ITS regions of H. occidentalis strains are presented in Supplementary Fig. S3 in IJSEM Online. Only the topology of the actin tree (Supplementary Fig. S3a) recovered the relationships obtained from the AFLP fingerprints (Fig. 2Go), whereas the positions of strains CBS 2569 and CBS 9921 in the EF-1{alpha} and ITS trees were incongruent. For example, strain CBS 9921, a member of group II in the actin and ITS trees (Supplementary Fig. S3a, c), was placed in group I in the EF-1{alpha} tree (Supplementary Fig. S3b). On the other hand, strain CBS 2569 was basal to group I in the actin tree but was basal to group II in the EF-1{alpha} and ITS trees. The sequence divergence in the D1/D2 regions was too low (0–1 substitutions) for phylogenetic reconstruction.

From the results presented, we might predict that the two groups of H. occidentalis are isolated by habitat preference and may present an example of an ongoing speciation event. Therefore, we propose that the two divergent groups within H. occidentalis should be recognized as varieties, with the names H. occidentalis var. occidentalis for group I and H. occidentalis var. citrica for group II.

Taxonomy

Latin diagnosis of Hanseniaspora pseudoguilliermondii N. Cadez, P. Raspor et M. Th. Smith sp. nov.
In medio liquido post 48 horas 25 °C cellulae apiculatae, ovoideae vel elongatae, 2.2–8.7x1.6–4.2 µm, singulae vel binae; gemmatione bipolari repruducentes. Post unum mensem annulus tenuis et sedimentum formantur. In agaro farina Zeae maydis confecto pseudomycelium rudimentarium. In quoque asco 4, ascosporae petasiformes. Glucosum et cellobiosum fermentantur. Glucosum, cellobiosum, salicinum, arbutinum, glucono-{delta}-lactonum, 2-ketogluconatum, acidum gluconicum, ethylaminum, lysinum et cadaverinum assimilantur. Non assimilantur galactosum, L-sorbosum, D-glucosaminum, D-ribosum, D-xylosum, L-arabinosum, D-arabinosum, L-rhamnosum, sucrosum, maltosum, trehalosum, methyl {alpha}-D-glucosidum, melibiosum, lactosum, raffinosum, melezitosum, amylum solubile, glycerolum, erythritolum, ribitolum, xylitolum, L-arabinitolum, D-glucitolum, D-mannitolum, galactitolum, inositolum, acidum D-glucuronicum, acidum D-galacturonicum, acidum DL-lacticum, acidum succinicum, acidum citricum, methanolum et ethanolum. Maxima temperatura crescentiae 37 °C. Crescit in medio addito 10 % NaCl, 50 % glucoso et 0.1 % cycloheximido. G+C acidi deoxyribonucleati 31.5 mol%. Typus CBS 8772T (=NCAIM Y.741T) in collectione zymotica Centraalbureau voor Schimmelcultures, Trajectum ad Rhenum lyophilus depositus.

Description of Hanseniaspora pseudoguilliermondii N. Cadez, P. Raspor & M. Th. Smith sp. nov.
Etymology: the epithet pseudoguilliermondii is chosen because the species is similar to H. guilliermondii.

In YM liquid medium after 48 h at 25 °C cells are apiculate, ovoid to elongate, 2.2–8.7x1.6–4.2 µm, single or in pairs. Budding is bipolar. Sediment is present. After 1 month, a very thin ring and a sediment are formed. After 1 month at 25 °C a streak culture on malt agar is cream coloured, butyrous, smooth, glossy, flat to slightly raised at the centre, with an entire to slightly undulate margin. On cornmeal agar a rudimentary pseudomycelium is formed. Asci containing four hat-shaped ascospores are observed on 5 % Difco malt extract agar at 25 °C. Glucose and cellobiose are fermented. The following carbon compounds are assimilated: glucose, cellobiose, salicin, arbutin, glucono-{delta}-lactone, 2-ketogluconate and D-gluconate. The following are not assimilated: galactose, L-sorbose, D-glucosamine, D-ribose, D-xylose, L-arabinose, D-arabinose, L-rhamnose, sucrose, maltose, trehalose, methyl {alpha}-D-glucoside, melibiose, lactose, raffinose, melezitose, starch, glycerol, erythritol, ribitol, xylitol, L-arabinitol, D-glucitol, D-mannitol, galactitol, myo-inositol, D-glucuronate, D-galacturonate, DL-lactate, succinate, citrate, methanol and ethanol. Assimilation of nitrogen compounds is positive for ethylamine, lysine and cadaverine but negative for sodium nitrate and sodium nitrite. Growth at 37 °C is positive; growth at 40 °C is negative. Growth on YM agar with 10 % NaCl is positive. Growth on 50 % glucose (w/w) yeast extract agar is weak. Growth in the presence of 0.1 % cycloheximide is positive. Diazonium blue B reaction is negative. G+C content of DNA (Tm) is 31.5 mol%.

The type strain, CBS 8772T (=NCAIM Y.741T), was isolated from orange juice concentrate, GA, USA. It has been deposited at the Centraalbureau voor Schimmelcultures, Utrecht, The Netherlands.

Hanseniaspora occidentalis Smith var. occidentalis
Antonie van Leeuwenhoek (1974) 40, 441–444.

Type strain: ex-type isolate CBS 2592T (=ATCC 32053T) is a living strain in the CBS yeast collection.

Hanseniaspora occidentalis var. citrica N. Cadez, P. Raspor & M. Th. Smith var. nov.
Etymology: the epithet citrica of Citrus, the genus name of the host plant.

Type strain: ex-type isolate CBS 6783T is a living strain in the CBS yeast collection. Varietas a Hanseniaspora occidentalis var. citrica differt: trehalosum assimilantur (inferme). Typus: CBS 6783T in collectione zymotica Centraalbureau voor Schimmelcultures, Trajectum ad Rhenum lyophilus depositus.

This variety differs from Hanseniaspora occidentalis var. occidentalis by assimilation of trehalose (weak). The type strain, CBS 6783T, was isolated from orange juice from Italy. It has been deposited at the Centraalbureau voor Schimmelcultures, Utrecht, The Netherlands.


    ACKNOWLEDGEMENTS
 
We thank Ann Vaughan-Martini for providing cultures. Vincent Robert is acknowledged for valuable comments on the manuscript.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
Baldauf, S. L. & Palmer, J. D. (1993). Animals and fungi are each other's closest relatives: congruent evidence from multiple proteins. Proc Natl Acad Sci U S A 90, 11558–11562.[Abstract/Free Full Text]

Baldauf, S. L., Roger, A. J., Wenk-Siefert, I. & Doolittle, W. F. (2000). A kindom-level phylogeny of eukaryotes based on combined protein data. Science 290, 972–977.[Abstract/Free Full Text]

Bandelj, D., Jakse, J. & Javornik, B. (2004). Assessment of genetic variability of olive varieties by microsatellite and AFLP markers. Euphytica 136, 93–102.[CrossRef]

Barnett, J. A., Payne, R. W. & Yarrow, D. (2000). Yeasts: Characteristics and Identification. Cambridge: Cambridge University Press.

Belloch, C., Querol, A., Garcia, M. D. & Barrio, E. (2000). Phylogeny of the genus Kluyveromyces inferred from the mitochondrial cytochrome-c oxidase II gene. Int J Syst Evol Microbiol 50, 405–416.[Abstract]

Boekhout, T., Kurtzman, C. P. O'Donnell K. & Smith, M. T. (1994). Phylogeny of the yeast genera Hanseniaspora (anamorph Kloeckera), Dekkera (anamorph Brettanomyces), and Eeniella as inferred from partial 26S ribosomal DNA nucleotide sequences. Int J Syst Bacteriol 44, 781–786.[Abstract/Free Full Text]

Bush, R. M. & Everett, K. D. E. (2001). Molecular evolution of the Chlamydiaceae. Int J Syst Evol Microbiol 51, 203–220.[Abstract]

Cadez, N., Raspor, P., de Cock, A. W. A. M., Boekhout, T. & Smith, M. Th. (2002). Molecular identification and genetic diversity within species of the genera Hanseniaspora and Kloeckera. FEMS Yeast Res 1, 279–289.[Medline]

Cadez, N., Poot, G. A., Raspor, P. & Smith, M. Th. (2003). Hanseniaspora meyeri sp. nov., Hanseniaspora clermontiae sp. nov., Hanseniaspora lachancei sp. nov. and Hanseniaspora opuntiae sp. nov., novel apiculate yeast species. Int J Syst Evol Microbiol 53, 1671–1680.[Abstract/Free Full Text]

Damgaard, J. & Cognato, A. I. (2003). Sources of character conflict in a clade of water strides (Heteroptera: Gerridae). Cladistics 19, 512–526.[CrossRef]

Daniel, H.-M. & Meyer, W. (2003). Evaluation of ribosomal RNA and actin gene sequences for the identification of ascomycetous yeasts. Int J Food Microbiol 86, 61–78.[Medline]

Daniel, H.-M., Sorrell, T. C. & Meyer, W. (2001). Partial sequence analysis of the actin gene and its potential for studying the phylogeny of Candida species and their teleomorphs. Int J Syst Evol Microbiol 51, 1593–1606.[Abstract]

Esteve-Zarzoso, B., Peris-Torán, M. J., Ramón, D. & Querol, A. (2001). Molecular characterisation of Hanseniaspora species. Antonie van Leeuwenhoek 80, 85–92.[CrossRef][Medline]

Farris, J. S., Kallersjo, M., Kluge, A. G. & Bult, C. (1995). Testing significance of incongruence. Cladistics 10, 315–319.[CrossRef]

Felsenstein, J. (1978). Cases in which parsimony and compatibility methods will be positively misleading. Syst Zool 27, 401–410.[CrossRef]

Felsenstein, J. (1985). Confidence limits on phylogenies: an approach using the bootstrap. Evolution 39, 783–791.[CrossRef]

Gardes, M. & Bruns, T. D. (1993). ITS primers with enhanced specificity for basidiomycetes – application to the identification of mycorrhizae and rusts. Mol Ecol 2, 113–118.[Medline]

Golubev, W. I., Smith, M. Th., Poot, G. A. & Kock, J. L. F. (1989). Species delineation in the genus Nadsonia Sydow. Antonie van Leeuwenhoek 55, 369–382.[CrossRef][Medline]

Hall, T. A. (1999). BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser 41, 95–98.

Kretzer, A. M. & Bruns, T. D. (1999). Use of atp6 in fungal phylogenetics: an example from the Boletales. Mol Phylogenet Evol 13, 483–492.[CrossRef][Medline]

Kumar, S., Tamura, K., Jakobsen, I. B. & Nei, M. (2001). MEGA2: molecular evolutionary genetics analysis software. Bioinformatics 17, 1244–1245.[Abstract/Free Full Text]

Kurtzman, C. P. (2003). Phylogenetic circumscription of Saccharomyces, Kluyveromyces and other members of the Saccharomycetaceae, and the proposal of the new genera Lachancea, Nakaseomyces, Naumovia, Vanderwaltozyma and Zygotorulaspora. FEMS Yeast Res 4, 233–245.[CrossRef][Medline]

Kurtzman, C. P. & Robnett, C. J. (1998). Identification and phylogeny of ascomycetous yeasts from analysis of nuclear large subunit (26S) ribosomal DNA partial sequences. Antonie van Leeuwenhoek 73, 331–371.[CrossRef][Medline]

Kurtzman, C. P. & Robnett, C. J. (2003). Phylogenetic relationships among yeasts of the ‘Saccharomyces complex’ determined from multigene sequence analyses. FEMS Yeast Res 3, 417–432.[CrossRef][Medline]

Lodder, J. & Kreger-van Rij, N. J. W. (1952). The Yeasts, a Taxonomic Study. Amsterdam: North Holland.

Meyer, S. A., Smith, M. Th. & Simione, F. P., Jr (1978). Systematics of Hanseniaspora Zikes and Kloeckera Janke. Antonie van Leeuwenhoek 44, 79–96.[CrossRef][Medline]

Misof, B., Rickert, A. M., Buckley, T. R., Fleck, G. & Sauer, K. P. (2001). Phylogenetic signal and its decay in mitochondrial SSU and LSU rRNA gene fragments of Anisoptera. Mol Biol Evol 18, 27–37.[Abstract/Free Full Text]

Möller, E. M., Bahnweg, G., Sandermann, H. & Geiger, H. H. (1992). A simple and efficient protocol for isolation of high molecular weight DNA from filamentous fungi, fruit bodies, and infected plant tissues. Nucleic Acids Res 20, 6115–6116.[Free Full Text]

Nakase, T. & Komagata, K. (1970). Significance of DNA base composition in the classification of yeast genera Hanseniaspora and Kloeckera. J Gen Appl Microbiol 16, 241–250.

O'Donnell, K., Cigelnik, E. & Nirenberg, H. I. (1998). Molecular systematics and phylogeography of the Gibberella fujikuroi species complex. Mycologia 90, 465–493.[CrossRef]

Posada, D. & Crandall, K. A. (1998). MODELTEST: testing the model of DNA substitution. Bioinformatics 14, 817–818.[Abstract/Free Full Text]

Reess, M. (1870). Botanische Untersuchungen über die Alkoholgährungspilze. Leipzig: A. Felix (in German).

Rokas, A., Williams, B. L., King, N. & Carroll, S. B. (2003). Genome-scale approaches to resolving incongruence in molecular phylogenies. Nature 425, 798–804.[CrossRef][Medline]

Smith, M. Th. (1998). Hanseniaspora Zikes. In The Yeasts, a Taxonomic Study, 4th edn, pp. 214–220. Edited by C. P. Kurtzman & J. W. Fell. Amsterdam: Elsevier.

Swofford, D. L. (2002). PAUP*: Phylogenetic Analysis Using Parsimony (*and other methods), version 4.0b10. Sunderland, MA: Sinauer Associates.

Taylor, J. W., Jacobson, D. J., Kroken, S., Kasuga, T., Geiser, D. M., Hibbett, D. S. & Fisher, M. C. (2000). Phylogenetic species recognition and species concepts in fungi. Fungal Genet Biol 31, 21–32.[CrossRef][Medline]

Thompson, J. D., Gibson, T. J., Plewniak, F., Jeanmougin, F. & Higgins, D. G. (1997). The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res 25, 4876–4882.[Abstract/Free Full Text]

Vos, P., Hogers, R., Bleeker, M. & 7 other authors (1995). AFLP: a new technique for DNA fingerprinting. Nucleic Acids Res 23, 4407–4414.[Abstract/Free Full Text]

White, T. J., Bruns, T., Lee, S. & Taylor, J. (1990). Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In PCR Protocols: a Guide to Methods and Applications, pp. 315–322. Edited by N. Innis, D. Gelfand, J. Sninsky & T. White. London: Academic Press.

Yamada, Y., Maeda, K. & Banno, I. (1992a). The phylogenetic relationships of the Q6-equipped species in the teleomorphic apiculate yeast genera Hanseniaspora, Nadsonia, and Saccharomycodes based on the partial sequences of 18S and 26S ribosomal ribonucleic acids. J Gen Appl Microbiol 38, 585–596.

Yamada, Y., Maeda, K. & Banno, I. (1992b). An emendation of Kloeckeraspora Niehaus with the type species, Kloeckeraspora osmophila Niehaus and the proposals of two new combinations, Kloeckeraspora occidentalis and Kloeckeraspora vineae (Saccharomycetaceae). Bull Jpn Fed Cult Collect 8, 79–85.

Yarrow, D. (1998). Methods for the isolation, maintenance and identification of yeasts. In The Yeasts, a Taxonomic Study, 4th edn, pp. 77–100. Edited by C. P. Kurtzman & J. W. Fell. Amsterdam: Elsevier.

Zikes, H. (1911). Zur Nomenclaturfrage der Apiculatushefen. Zentralbl Bakteriol Parasitenkd Abt II 30, 145–149 (in German).





This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplementary Material
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via CrossRef
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Cadez, N.
Right arrow Articles by Smith, M. Th.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Cadez, N.
Right arrow Articles by Smith, M. Th.
Agricola
Right arrow Articles by Cadez, N.
Right arrow Articles by Smith, M. Th.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
INT J SYST EVOL MICROBIOL MICROBIOLOGY J GEN VIROL
J MED MICROBIOL ALL SGM JOURNALS