Multilocus phylogenetic study of the Scheffersomyces yeast clade and characterization of the N-terminal region of xylose reductase gene.

Many of the known xylose-fermenting (X-F) yeasts are placed in the Scheffersomyces clade, a group of ascomycete yeasts that have been isolated from plant tissues and in association with lignicolous insects. We formally recognize fourteen species in this clade based on a maximum likelihood (ML) phylogenetic analysis using a multilocus dataset. This clade is divided into three subclades, each of which exhibits the biochemical ability to ferment cellobiose or xylose. New combinations are made for seven species of Candida in the clade, and three X-F taxa associated with rotted hardwood are described: Scheffersomyces illinoinensis (type strain NRRL Y-48827(T) = CBS 12624), Scheffersomyces quercinus (type strain NRRL Y-48825(T) = CBS 12625), and Scheffersomyces virginianus (type strain NRRL Y-48822(T) = CBS 12626). The new X-F species are distinctive based on their position in the multilocus phylogenetic analysis and biochemical and morphological characters. The molecular characterization of xylose reductase (XR) indicates that the regions surrounding the conserved domain contain mutations that may enhance the performance of the enzyme in X-F yeasts. The phylogenetic reconstruction using XYL1 or RPB1 was identical to the multilocus analysis, and these loci have potential for rapid identification of cryptic species in this clade.

Introduction

D-xylose is a five-carbon backbone molecule of the hemicellulose component of plant cell walls and is one of the most abundant renewable carbon resources on Earth. Some bacteria and certain fungi, including fewer than twenty-five species of the more than 1500 described ascomycete yeasts, share the ability to produce ethanol by the fermentation of D-xylose [1]. In order to ferment D-xylose, yeasts express xylose reductase (XR), xylitol dehydrogenase (XDH), and xylulose kinase (XK) to convert D-xylose to D-xylulose-5-phosphate; D-xylulose-5-phosphate is then incorporated into the pentose phosphate pathway to be catalyzed to ethanol [2]–[4].

Xylose fermentation has been the focus of several studies in order to identify differences in the catabolic rate between strains [2], [5]-[14]. Overexpression, homologous and heterologous expression, and direct mutagenesis of genes involved in D-xylose assimilation and fermentation have only modestly enhanced the quantity of ethanol production by yeasts due to several metabolic constraints; these include rate of regeneration of the cofactor NADP(H) required by XR and XDH, repression by glucose, and anaerobic respiration regulatory control [15]-[17]. These studies have resulted in the present understanding of the biochemical pathway, but the main goal of bioengineering yeasts capable of fermenting D-xylose at a high rate to be used at industrial scales has not yet been achieved. Consequently, recent research has been focused on the discovery of new X-F yeasts, e.g. Spathaspora passalidarum and Candida jeffriesii [18], and Spathaspora arborariae and other taxa [19], [20], from the guts of lignicolous beetles and rotted wood, niches from which a number of X-F yeasts have been isolated [21]-[23].

Although, X-F yeasts appear scattered throughout the Saccharomycotina, the yeasts that have been reported to exhibit the highest rate of xylose fermentation under certain conditions are members of the Scheffersomyces clade [9], [21], [24], and for this reason fermentative ability has been intensively studied in this clade [9], [25]-[51]. Few studies, however, have undertaken clarification of the phylogenetic relationships among the yeasts of the clade [1], [52]-[58]. Therefore, a robust, well-supported phylogeny including as many taxa as possible is necessary to clarify the phylogenetic relationships among the Scheffersomyces clade members along with a comparison of the nucleotide mutations and enzymatic activity of the XR to understand the importance of the biochemical ability in the speciation process of this yeast clade.

In order to distinguish the species in the Scheffersomyces clade we used BLAST searches, biochemical and morphological characterization, and a multilocus phylogenetic analysis that included the traditional SSU and LSU markers, the orthologous RPB1, and the recently proposed ITS barcoding region for fungi [59]. We present a taxonomic revision of the Scheffersomyces clade, trace the nucleotide differences in the XYL1, and propose three new species of X-F yeasts, Scheffersomyces illinoinensis, Scheffersomyces quercinus, and Scheffersomyces virginianus associated with rotted hardwoods, Carya illinoinensis (pecan), Quercus nigra (water oak), and Quercus virginiana (live oak). We performed a molecular study of the XYL1 that codifies XR in certain members of the Scheffersomyces clade.

Materials and Methods

Yeast Isolation and Culture

Partially decayed logs and fallen branches of Carya illinoinensis (pecan, 30 cm diam), Quercus nigra (water oak, approximately 10 cm diam), and Quercus virginiana (live oak, approximately 10 cm diam) were collected from Pecan Drive, Saint Gabriel, Ascension Parish, Louisiana, and LSU Burden Center and the corner of Highland Road and S. Stadium Drive on the LSU campus, Baton Rouge, East Baton Rouge, Louisiana, respectively, between Sep and Oct 2007. The wood samples were divided into approximately 1 cm2 samples, and each sample was placed in a 1.5 mL microcentrifuge tube with 1 mL of sterile water. Tubes were vortexed for 30 s and a 100 µL aliquot was plated on acidified yeast medium agar [60]. Plates were incubated for 3 d at 25 C, and single colonies were isolated and streaked 4 times to obtain pure cultures.

Molecular Identification

Genomic DNA was extracted using a Wizard® Genomic DNA purification kit (Promega). The concentration, integrity, and purity of total DNA extracted were confirmed by gel electrophoresis in 0.8% agarose in 0.5 × Tris-Borate-EDTA (TBE) buffer. Initial rapid identification was carried out by PCR amplification and sequencing of the LSU (D1/D2 region ∼600 bp) rRNA gene for use in BLAST searches [1], [54]. In order to increase the robustness of the phylogenetic analyses, PCR amplifications of the small subunit (SSU ∼1.6 Kbp) and internal transcribed spacers 1 and 2 (ITS ∼500 bp) of the rRNA marker were carried out in addition to the D1/D2 region [61], [62] and the SSU rRNA gene was amplified using the combination of primers NS1 (forward) (5′-GTAGTCATATGCTTGTCTC-3′) and NS8 (reverse) (5′-TCCGCAGGTTCACCTACGGA-3′); ITS1-LSU markers were amplified using the combination of primers ITS1 (forward) (5′-TCCGTAGGTGAACCTGCGG-3′) and LR3 (reverse) (5′-CCGTGTTTCAAGACGGG-3′) in a PCR reaction with 20 µg of total DNA, 0.5 mM DTPs, 2.5 mM MgSO4, 0.8 µM of each primer set, 1 × PCR buffer, and 1 U of Taq polymerase (Promega) in 25 µL of final volume. The PCR amplification protocol included 5 min of DNA pre-denaturation at 95 C followed by 35 cycles of 1 min of DNA denaturation at 95 C, 45 s of primer annealing at 55 C, and a 2 min extension at 72 C, and 10 min final PCR extension.

In addition, another nuclear locus, RNA polymerase II, was used in the phylogenetic analysis. A fragment of ∼700 bp of the subunit I (RPB1) gene was amplified by the primer pair RPB1-Af (forward) 5′-GARTGYCCDGGDCAYTTYGG-3′ and RPB1-Cr (reverse) 5′-CCNGCDATNTCRTTRTCCATRTA-3′ [63], [64]. The PCR reaction was performed using 100 µg of total DNA, 0.6 mM DTPs, 2.5 mM MgSO4, 1 µM of each primer, 1 × PCR buffer, and 1.5 U of Taq polymerase (Promega) in 35 µL total final volume of reaction. The PCR amplification program included 5 min of DNA pre-denaturation at 95 C followed by 35 cycles of 1 min of DNA denaturation at 95 C, 45 s of primer annealing at 55 C, and 2 min of extension at 72 C, and 10 min final PCR extension.

The purified PCR products were sequenced in both directions by Beckman Coulter Genomics (Danvers, MA). Each molecular marker was sequenced on three independent occasions in order to avoid nucleotide differences due to sequencing errors.

RAPD-PCR Fingerprinting

RAPD–PCR analysis was performed with the oligonucleotide primer CDU (5′-GCGATCCCCA-3′) [65]-[68]. Aliquots of 25 µL of amplified product were analyzed by electrophoresis on 1.8% agarose gel in 1× TBE buffer with 1 × SYBR® Safe DNA Gel Stain (Invitrogen, Grand Island, NY) at 70 V for 80 min. DNA fragments were visualized with a UV-light transilluminator and photographed using a Polaroid system.

Morphological, Biochemical and Physiological Characteristics

The yeast standard description based on phenotypic characters was executed following standardized protocols [1], [69], [70].

Xylose Reductase (XR) Molecular Studies

The ∼600 bp fragment of the XYL1 was amplified using the following degenerate primers: XYL1-forward (5′-GGTYTTYGGMTGYTGGAARSTC-3′) and XYL1-reverse (5′-AAWGATTGWGGWCCRAAWGAWGA-3′) designed in this study, in a PCR reaction with 100 µg of total DNA, 0.4 mM DTPs, 4 mM MgSO4, 1 µM of each primer, 1× PCR buffer, and 1U of Taq polymerase (Promega) in a final volume of 25 µL. The PCR amplification program included 5 min of DNA pre-denaturization at 95 C followed by 35 cycles of 1 min of denaturization at 95 C, 1 min primer annealing at 57 C, 1 min extension at 72 C, and 10 min final PCR extension. The purified PCR products were sequenced as described above.

Phylogenetic Analyses

Contig sequences and sequencing manipulations were performed with Se-AL v2.01a11 (http://tree.bio.ed.ac.uk/software/seal/) and MESQUITE v2.74 [71]. The sequence alignments were carried out using the online interface MAFFT v6.859 (http://mafft.cbrc.jp/alignment/software/) with different advanced alignment strategies per locus: LSU and RPB1, global homology (G-INS-i); ITS, one conserved domain (L-INS-i); and SSU, secondary structure of RNA (Q-INS-i). In particular ITS loci were realigned using the software SATé v2.1.2 [72], and ambiguous sequence alignment ends were eliminated in all the alignments. Maximum likelihood (ML) phylogenic inference was performed in RAxML-VI-HPC [73] using a partitioned multilocus matrix (each partition for SSU, ITS1, 5.8 S, ITS2, LSU, and RPB1) under a general time reversible model with a gamma distribution of site rate variation (GTRGAMMA), and ML support was estimated using 1000 bootstrap replicates. The sequences for RPB1 and XR were obtained using SEQUIN v11.0 (http://www.ncbi.nlm.nih.gov/Sequin/) with the alternative yeast nuclear codon bias [45], [74]. Phylogenetic analysis using the amino acid matrix was performed in RAxML-VI-HPC under the VT AA selection model [75] and ML support was estimated using 1000 bootstrap replicates. Tree editing was done with FigTree v1.3.1 software (http://tree.bio.ed.ac.uk/software/figtree/).

Results

Diversity of Yeasts Isolated from Rotted Wood

The initial rapid molecular identification of the 29 yeast strains isolated from the wood samples using the D1/D2 LSU region confirmed the presence of several closely related species classified in the Scheffersomyces (11 isolates), Sugiyamaella (10 isolates), Trichomonascus (5 isolates), Meyerozyma (2 isolates), and Candida tanzawaensis (1 isolate) clades (Fig. 1). Initial biochemical characterization of fermentation abilities performed on all the isolates showed that only strains classified in the Scheffersomyces clade had the ability to ferment xylose.

Figure 1

ML tree based of the D1/D2 LSU region using a 606-character matrix for yeast species isolated from the wood samples (in bold).

Schizosaccharomyces pombe was used as an outgroup taxon (in grey). X-F, xylose-fermenting yeasts. Numbers above each branch refer to bootstrap values out of 1000 repetitions. ML score -11353.90.

New Xylose-fermenting Yeasts

Species delimitation of S. quercinus, S. virginianus, and S. illinoinensis, was based on phylogenetic placement, nucleotide differences in the rRNA markers (SSU, ITS and LSU), RPB1 and XYL1, and biochemical assay differences compared with their closest relatives (Table 1, 2 and 3, Fig. 2). In addition, we obtained different banding patterns using CDU RAPID-PCR fingerprinting primers (Fig. 3). These fingerprinting primers have been used previously to discriminate among cryptic yeast species [66], [67], [76].

Table 1

GenBank accession numbers of the nucleotide sequences used in this study. Sequences generated in this work shown in bold.

Species	Codes	SSU	ITS	LSU	RPB1	XYL1
C. bolitotheri T	NRRL Y-27587, CBS 9832, BG 00-8-15-1-1	AY242142	FJ623599	AY242249	JN804828	-
C. terraborum T	NRRL Y-27573, CBS 9826, BG 02-7-15-019A-2-1	AY426956	FJ623596	AY309810	JN804831	-
C. panamericana T	NRRL Y-27567, CBS 9834, BG 01-7-26-006B-2-1	AY426960	FJ623601	AY309872	JN804835	-
S. coipomoensis T	NRRL Y-17651, ATCC 58904, CBS 8178	HQ651931	HQ652070	HQ651966	EU344070	-
S. lignicola T	ATCC MYA-4674, CBS 10610, BCC 7733	AY845351	HQ652074	AY845350	-	-
S. ergatensis T	NRRL Y-17652, ATCC 22589, CBS 6248	AB013524	EU343826	U45746	EU344098	JQ436926
S. insectosa T	NRRL Y-12854, ATCC 66611, CBS 4286	AB013583	HQ652064	FM200041	JN804842	JQ235697
S. lignosus T	NRRL Y-12856, ATCC 58779, CBS 4705	HQ651941	JN943262	U45772	JN804837	JQ235693
S. segobiensis T	NRRL Y-11571, ATCC 58375, CBS 6857	AB054288	DQ409166	U45742	EF599429	JQ436925
L. elongisporus T	NRRL YB-4239, ATCC 11503, CBS 2605	HQ876033	HQ876042	HQ876050	AY653537	-
C. tropicalis T	NRRL Y-12968, ATCC 4563, CBS 616	EU348785	AB437068	U45749	-	-
S. queiroziae T	NRRL Y-48722, UFMG-CLM 5.1, CBS 11853	-	HM566445	HM566445	-	-
S. gosingicus T	CBS-11433, BCRC 23194, SJ7S11	HQ876040	HQ999978	HQ999955	-	-
S. spartinae T	NRRL Y-7322, ATCC 18866, CBS 6059	FJ153139	HQ876044	U45764	-	-
S. stipitis T	NRRL Y-7124, ATCC 58376, CBS 5773	AB054280	JN943257	U45741	JN804841	JQ235696
Scheffersomyces sp.	NRRL Y-48762, CBS 12363, UFMG HMD-26.3	-	JF826438	JF826438	-	-
S. shehatae T	NRRL Y-12858, CBS 5813, ATCC 34887	AB013582	JN943264	AF178049	JQ436927	JQ235691
S. quercinus T	NRRL Y-48825T, CBS 12625, W07-09-15-1-3-2	JN940981	JN943260	JN703957	JN804838	JQ008829
S. virginianus T	NRRL Y-48822T, CBS 12626, W07-10-04-4-6-2	JN940969	JN943259	JN703958	JN804839	JQ235695
S. illinoinensis T	NRRL Y-48827T, CBS 12624, W07-11-15-9-2-1	JN940968	JN943261	JN703959	JN804840	JQ235694

Table 2

Nucleotide differences and percentages of homology between the new xylose-fermenting yeasts and the type cultures of closest relatives, S. shehatae or S. stipitis.

Species	SSU	ITS	D1/D2 LSU	RPB1	XYL1
S. insectosa T	100%	98% (5 n)	100%	94% (34 n)	92% (44 n)
S. lignosus T	100%	98% (5 n)	100%	94% (37 n)	91% (46 n)
S. quercinus T	98% (27n)	99% (1 n)	98% (7 n)	98% (10 n)	96% (11 n)
S. virginianus T	97% (54n)	99% (1 n)	99% (4 n)	98% (10 n)	97% (8 n)
S. illinoinensis T	99% (3n)	99% (3 n)	100%	95% (28 n)	97% (15 n)

Table 3

Differences in physiological reactions of S. quercinus, S. illinoinensis, and S. virginianus and their closest relatives*.

Species	Rhammose	Galactitol	Lactate	Cadaverine	Glucosamine	Tryptophan	Melibiose	Inulin
S. insectosa	-	-	-	+	-	-	-	-
S. lignosus	-	+	w	+	-	-	-	-
S. quercinus	-	+,f	+	+	+	+	-	+
S. shehatae	-	-	-	+	-	-	-	-
S. virginianus	w	w	-	-	+	+	+	+
S. stipitis	+	-	+/−	+	-	-	-	-
S. illinoinensis	+	+/d	+/d	w	+	+	w	w

Biochemical assay results of previously described species were complied from Kurtzman (1990) and Barnett et al. (2000). Abbreviations for reaction results: +, positive; -, negative; d, delayed positive; w, weak positive.

Figure 2

ML tree based on a multilocus dataset using a 3488-character matrix for Scheffersomyces clade.

Candida tropicalis was used as an outgroup taxon (in grey). C-F and X-F, cellobiose- and xylose-fermenting yeasts respectively. Numbers above each branch refer to bootstrap values out of 1000 repetitions. ML score is -13300.52.

Figure 3

Characterization of S. quercinus, S. illinoinensis, and S. virginianus using RAPID-PCR CDU fingerprinting primers.

Species Description of Scheffersomyces quercinus H. Urbina & M. Blackw. sp. nov

MycoBank accession number MB563719, (Fig. 4, a-b).

Figure 4

Morphological characterization of S. quercinus, budding cells (a-b); S. virginianus, budding cells (c-d); and S. illinoinensis, budding cells (e-f); grown at 25 C in YM broth and V8 agar at 7 d days respectively, bar 10 µm.

After 7 d growth in YM broth at 25 C, cells are subglobose (5–8×5–7.5 µm), and occur singly, in pairs, or in chains. Pseudohyphae are present; true hyphae are absent. After 7 days on YM agar at 25 C, colonies are cream-colored with pale-pinkish perimeter on some older colonies, smooth, flat, and/or with scattered filaments at the margin. After 10 d of Dalmau plate culture on corn meal agar at 25 C, true hyphae are present. Aerobic growth is white, shiny, and smooth with filamentous margin. Asci and ascospores are not observed on YM or V8 agar. Diazonium blue B reaction is negative. See Table 4 for physiological characterization.

Table 4

Biochemical characterization of S. quercinus, S. illinoinensis, and S. virginianus.

		S. quercinus	S. illinoinensis	S. virginianus
	Fermentation
F1	D-Glucose	+	+	+, f
F2	D-Galactose	+	+	+
F3	Maltose	w	+	+
F4	α-Methyl-D-glucoside	-	-	-
F5	Sucrose	-	-	-
F6	α,α- Trehalose	+	+	+
F7	Melibose	-	w	-
F8	Lactose	-	-	-
F9	Cellobiose	-	-	-
F10	Melezitose	-	-	-
F11	Raffinose	-	-	-
F12	Inulin	-	-	-
F13	Starch	-	-	-
F14	D-Xylose	+	+	+
	Assimilation
C1	D-Glucose	+, f	+, f	+
C2	D-Galactose	+, f	+, f	+, f
C3	L-Sorbose	+,d	-	w
C4	D-Glucosamine	w	-	+
C5	D-Ribose	+,d	+,d	+, d
C6	D-Xylose	+	+	+
C7	L-Arabinose	+,d	+	+
C8	D-Arabinose	w	+	w
C9	L-Rhamnose	-	+	w
C10	Sucrose	+	+	+
C11	Maltose	+	+, f	+
C12	Trehalose	+	+, f	+
C13	α-Methyl-D-glucoside	+	+	+
C14	Cellobiose	+	+	+
C15	Salicin	+	+	+
C16	Arbutin	+	+	+
C17	Melibiose	-	+	+
C18	Lactose	-	-	+
C19	Raffinose	-	-	w
C20	Melezitose	-	+	+, d
C21	Inulin	w	w	+
C22	Soluble Starch	+	+	+
C23	Glycerol	+	+	+
C24	Erythritol	+	+	+
C25	Ribitol	+	+	+
C26	Xylitol	+, d	+	w
C27	L-Arabinitol	+, d	+,d	-
C28	D-Glucitol	+	+	+
C29	D-Mannitol	+	+	+
C30	Galactitiol	+, f	+, f	w
C31	myo-Inositol	-	-	w
C32	D-Glucono-1,5-lactone	+	+	+
C33	2-Keto-D-gluconate	+	+,f	+, f
C34	5-Keto-D-gluconate	?	?	?
C35	D-Gluconate	+	+,d	+
C36	D-Glucuronate	-	-	-
C38	DL-Lactate	+	+,d	-
C39	Succinate	+	+	+
C40	Citrate	+, d	+	+
C41	Methanol	-	-	-
C42	Ethanol	+, f	+,f	+, f
C43	Propane 1,2 diol	+, d	w	w
C44	Butane 2,3 diol	w	-	w
C45	Quinic acid	-	-	w
C46	D-Glucarate	-	-	-
	Temperature
T1	30°C	+	+	+
T2	35°C	+, f	+,f	+
T3	37°C	-	+,d	-
	Osmotic pressure
O1	0.01% Cycloheximide	+, f	+, f	+, f
O2	0.1% Cycloheximide	+, f	+, f	+,f
O3	1% Acetic Acid	-	-	-
O4	50% D-Glucose	+	+	+
O5	60% D-Glucose	w	+	w
O6	10% NaCl	+	+	w
O7	16% NaCl	-	-	-
	Nitrogen assimilation
N1	Nitrate	-	+, d	+
N2	Nitrite	-	+, d	-
N3	Ethylamine	+	+	+
N4	L-Lysine	+, d	-	+
N5	Cadaverine	+	w	-
N6	Creatine	-	-	-
N7	Creatinine	-	-	-
N8	D-Glucosamine	+	+	+
N9	Imidazole	-	-	-
N10	D-Tryptophan	+	+	+
	Vitamins
V1	w/o Vitamins	-	-	-
V2	w/o myo-Inositiol	+, d	w	+
V3	w/o Pantothenate	+	+	+
V4	w/o Biotin	+	-	-
V5	w/o Thiamin	+	+	+
V6	w/o Biotin & Thiamin	-	-	w
V7	w/o Pyridoxine	+	+	+
V8	w/o Pyrid. &Thiam	+	+	+
V9	w/o Niacin	+	+	+
V10	w/o PABA	+	+	+

Abbreviations for reaction results: +, positive; -, negative; d, delayed positive; w, weak positive.

Type strain

NRRL Y-48825T ( = CBS 12625; W07-09-15-1-3-2) is preserved as a lyophilized preparation in the Agricultural Research Service Culture Collection (NRRL), Peoria, Illinois, USA. The strain was isolated from rotted wood (Quercus niger), collected on 15 Nov 2007 at Louisiana State University, Burden Research Plantation, Baton Rouge, East Baton Rouge Parish, Louisiana, USA, N30°40′59″-W91°10′29″.

Etymology

The species name quercinus (N.L. gen. n.) refers to the genus of the substrate, Quercus niger, from which this species was isolated.

Species Description of Scheffersomyces virginianus H. Urbina & M. Blackw. sp. nov

MycoBank accession number MB 563720 (Fig. 4, c-d).

After 7 d growth in YM broth at 25 C, cells are globose to ellipsoidal (5.5–10×4–6.5 µm), and occur singly, in pairs, in short chains, or in small clusters. Pseudohyphae are present; true hyphae are absent. After 7 d on YM agar at 25 C, cream-colored to light pink with abundant filaments at margin. After 10 d Dalmau plate culture on corn meal agar at 25 C, pseudohyphae are present; septate hyphae are absent. Aerobic growth is white, shiny, and smooth with filaments at margin. Asci and ascospores are not observed on YM or V8 agar. Diazonium blue B reaction is negative. See Table 4 for physiological characterization.

Type Strain

NRRL Y-48822T ( = CBS 12626; W07-10-04-4-6-2) is preserved as a lyophilized preparation in the Agricultural Research Service Culture Collection (NRRL), Peoria, Illinois, USA. The strain was isolated from rotted wood (Quercus virginiana) collected on 15 Sep 2007, Highland Rd. and S. Stadium Dr., Louisiana State University, Baton Rouge, East Baton Rouge Parish, Louisiana, USA, N30°40′90″-W91°17′60″.

Etymology

The species name virginianus (N.L. gen. n.) refers to the species of the substrate, Quercus virginiana, from which this species was isolated.

Species Description of Scheffersomyces illinoinensis H. Urbina & M. Blackw. sp. nov

MycoBank accession number MB563457 (Fig. 4, e-f).

After 7 d growth in YM broth at 25 C, cells are globose (5−7.5 µm), and occur singly, in pairs, or in short chains. Pseudohyphae and hyphae are not present. After 7 d on YM agar at 25 C, colonies are cream-colored, smooth, and flat with a smooth margin. After 10 d Dalmau plate culture on corn meal agar at 25 C, pseudohyphae are present; septate hyphae are absent. Asci and ascospores are not observed on YM and V8 agar. Diazonium blue B reaction is negative. See Table 4 for physiological characterization.

NRRL Y-48827T ( = CBS 12624; W07-11-15-9-2-1) is preserved as a lyophilized preparation in the Agricultural Research Service Culture Collection (NRRL), Peoria, Illinois, USA. The strain was isolated from rotted wood (Carya illinoinensis) collected on 15 Nov 2007 at 810 Pecan Dr., St. Gabriel, Ascension Parish, Louisiana, USA, N30°30′56″-W91°10′30″.

The species name illinoinensis (N.L. gen. n.) refers to the species of the substrate, Carya illinoinensis, from which this species was isolated.

Multilocus Phylogenetic Study

The phylogenetic placement of S. quercinus, S. virginianus, and S. illinoinensis was based on ML analysis results of a concatenated nucleotide dataset containing 3488 characters (Fig. 2). The Scheffersomyces clade was divided into three subclades: 1) the early diverging S. spartinae and S. gosingicus subclade, 2) the cellobiose-fermenting S. ergatensis subclade, and 3) the largest, xylose-fermenting S. stipitis subclade to which the three new species belong.

Addition of more taxa and more molecular data in phylogenetic analyses has helped to define monophyletic clades, among a number of genera now recognized as polyphyletic. One problematic taxon, Pichia, was based primarily on ascospore shape. Hat-shaped ascospores, however, are found among several distant clades of yeasts as well as distantly related members of the Pezizomycotina. Scheffersomyces was proposed recently for species in the Pichia stipitis clade [55]. The genus included the type species, S. stipitis, and S. ergatensis and S. spartinae [55]. We propose additional new combinations in the genus Scheffersomyces by including clade members that previously were described as asexual species of the polyphyletic genus Candida [77].

(C. Ramírez & A. González) H. Urbina & M. Blackw. comb. nov.

MycoBank accession number: MB563714.

Basionym: Candida coipomoensis C. Ramírez & A. González, Mycopathologia 88:84, 1984.

(C.F. Lee) H. Urbina & M. Blackw. comb. nov.

MycoBank accession number: MB563799.

Basionym: Candida gosingica C.F. Lee, Int J Syst Evol Microbiol 61:670, 2011.

(Kurtzman) H. Urbina & M. Blackw. comb. nov.

MycoBank accession number: MB 564799.

Basionym: Candida shehatae var. insectosa Kurtzman, Antonie van Leeuwenhoek 57:218, 1990. Synonym: Candida insectosa (Kurtzman) Kurtzman, in Kurtzman & M. Suzuki, Mycoscience 51:10, 2010.

(Jindam., Limtong, Yongman., Tuntir., Potach., H. Kawas & Nakase) H. Urbina & M. Blackw. comb. nov.

MycoBank accession number: MB 563803.

Basionym: Candida lignicola Jindam., Limtong, Yongman., Tuntir., Potach., H. Kawas & Nakase, FEMS Yeast Res 7:1412, 2007.

(Kurtzman) H. Urbina & M. Blackw. comb. nov.

MycoBank accession number: MB 563801.

Basionym: Candida shehatae var. lignosa Kurtzman, Antonie van Leeuwenhoek 57:218, 1990. Synonym: Candida lignosa (Kurtzman) Kurtzman, in Kurtzman & M. Suzuki, Mycoscience 51:10, 2010.

(R.O. Santos, R.M. Cadete, Badotti, A. Mouro, Wallheim, F.C.O. Gomes, Stambuk, Lachance & C.A. Rosa) H. Urbina & M. Blackw. comb. nov.

MycoBank accession number: MB563717.

Basionym: Candida queiroziae R.O. Santos, R.M. Cadete, Badotti, A. Mouro, Wallheim, F.C.O. Gomes, Stambuk, Lachance & C.A. Rosa, Antonie van Leeuwenhoek 99:639. 2011.

(H.R. Buckley & van Uden) H. Urbina & M. Blackw. comb. nov.

MycoBank accession number: MB563716.

Basionym: Candida shehatae H.R. Buckley & van Uden, Mycopathologia Mycol. appl. 34:297, 1966.

Xylose reductase among X-F yeast members of the S. stipitis subclade

We amplified a ∼600 bp PCR product of XYL1 from X-F and non X-F yeasts tested (Fig. 5). The translated protein sequences at the N-terminal region have the conserved amino acids 49-D, 51-A, and 54-Y, described as part of the catalytic GX3DXAX2Y domain; the LX8DX4H and the GX3GXG domains, and the amino acids 83-K, 132-P, and 167-K that form the xylose-binding pocket previously reported [78], [79] (Fig. 6). Conserved amino acid substitutions at positions 64, 170, 194 (I  =  V); 65 (N  =  K); 69 (K  =  D); 70 (D  =  E); 80, 178 (I  =  L); 107 (L  =  V); 127 (N  =  E); 149 (L  =  I), 150 (A  =  L), 174 (T  =  P), 179 (Y  =  L), 184 (S  =  G), 186 (E  =  T), 188 (P  =  K), and 204 (S  =  P) and the percentage of conserved substitutions usually were higher than 50% at all sites compared with the non X-F fermenting yeast S. ergatensis (Fig. 6).

Figure 5

ML consensus tree based on XYL1 and the putative XR from of X-F members in S. stipitis subclade.

Scheffersomyces ergatensis was used as an outgroup taxon (in grey). Numbers above each branch refer to bootstrap values out of 1000 repetitions. ML scores are -1829.90 (DNA data) and -1284.79 (amino acid data).

Figure 6

Multiple sequence alignment of the N-terminal region of XR.

Identical residues are indicated by dots. Arrows indicate amino acids that constitute the active sites, stars indicate amino acids that form the xylose-binding pocket, and brackets indicate conserved domains.

Discussion

Yeasts Sampled from Rotted Wood

Most of the isolates from live oak and water oak were members of the Sugiyamaella clade: Candida boreocaroliniensis, Candida lignohabitans, and S. smithiae. The cosmopolitan genus Sugiyamaella is comprised of yeasts reported primarily from wood and frass of lignicolous beetles. It is of interest that, unlike other yeasts from these habitats, they are unable to ferment D-xylose [80], [81]. We also isolated other non X-F species Candida athensensis, Trichomonascus petasosporus, and a close relative of Candida anneliseae, ascomycete yeasts that previously were found associated with fungus-feeding beetles collected in Panama and the USA [56].

The yeast strains isolated from pecan wood, on the contrary, were dominated by the new X-F yeast species, S. illinoinensis, a close relative of S. stipitis. Only this wood was inhabited by Odontotaenius disjunctus (Passalidae, Coleoptera) in our study. Zhang et al. [82] recognized the relationship between S. stipitis and the lignicolous beetle commonly found inhabiting decayed hardwoods in the southeastern US, so the phylogenetic placement of these closely related species is consistent with the previous findings.

Other members of the Scheffersomyces clade (S. ergatensis, S. shehatae, and S. stipitis) have been reported frequently from associations with the gut of lignicolous beetles, including, not only the passalid beetle O. disjunctus but also beetles in the families Cerambycidae, Lucanidae, Buprestidiae, and Tenebrionidae [22], [83]. It is likely that the common gut yeasts are efficient at digesting components of the host diet, resisting toxic secondary metabolites, and adapting to the gut physiological environment (low oxygen and high carbon dioxide concentrations and extreme pH variation), characteristics that give them a greater chance to be horizontally transmitted to progeny. Consequently, in each host generation the symbiotic yeasts may be exposed to bottlenecks and positive selection driven by the host beetles, and these selective pressures increase when changes in the host diet occur. These evolutionary processes could have favored rapid speciation with morphological and other traditional characters lagging behind molecular changes in the Scheffersomyces yeast members. Evidence that supports our hypothesis is: 1) S. shehatae, S. lignicola, and S. insectosa, often found in association with insects, are indistinguishable by morphology and some molecular markers (e.g. SSU and D1/D2 LSU); 2) branch lengths are constrained in the phylogenetic tree (Fig. 2); and 3) gut morphology is modified to enhance the horizontal transmission of gut yeasts across generations, e.g. the posterior hindgut region of O. disjunctus is colonized mainly by filamentous yeasts attached by a holdfast [84]; in addition mycetomes occur in lignicolous cerambycid beetles colonized exclusively by closely related species of S. shehatae [22].

Kurtzman [58] recognized Candida shehatae var. shehatae, var. insectosa, and var. lignicola based on biochemical assays, a single nucleotide difference in the D1/D2 LSU region, and identical SSU rRNA among the varieties. He suggested that in order to understand the phylogenetic relationship between varieties of C. shehatae, analyses including ITS rRNA should be included. More recently, these varieties were raised to species level based on their distinctive electrokariotype profiles and the reinterpretation of the D1/D2 LSU locus [55], [57]. Yeast species often have been underestimated on the basis of only LSU and SSU rRNA data [85]-[87]. Therefore, the addition of more molecular markers in phylogenetic analyses has been used to increase the power of species recognition (see next section).

Phylogenetic Study of the Scheffersomyces Clade

Lachance et al. [88] proposed a method for species delimitation in yeasts based on parsimony networks. In our experience, the results obtained by using this method are difficult to interpret for several reasons: 1) the species are plotted as isolated entities with little information on phylogenetic relationship among members; 2) the method does not allow inclusion of a model of nucleotide selection for the analysis; and 3) node support values are lacking.

Because of these disadvantages, we implemented, instead, a multilocus phylogenetic analysis based on ML, commonly applied to fungi. We incorporated additional loci (ITS and RPB1) in order to increase the robustness of the phylogenetic analysis in the study of the Scheffersomyces members. We followed the recommendations of the Assembling the Fungal Tree of Life (AFTOL http://aftol.org/) research group in searching for orthologous genes in fungal genomes [89]. In addition results obtained by Schoch et al. [90] showed that RPB1, RPB2, and TEF1 (elongation factor 1 alpha) are more phylogenetically informative compared to rRNA genes in Ascomycota. More recently, Schoch et al. [59] proposed the use of ITS as a barcode gene for fungi, although RPB1 had higher species discriminatory power than ITS in Saccharomycotina. The authors [59] also pointed out that datasets containing combinations of at least three molecular markers (SSU, LSU, ITS or RPB1) showed the highest probability of correct identification for all Fungi.

As mentioned in the results section, the Scheffersomyces clade, including all known taxa reported in the literature, is comprised of three subclades: 1) an early diverged subclade of S. spartinae, the only member that lacks the ability to ferment both xylose and cellobiose, and S. gosingicus, a cellobiose-fermenting (C-F) yeast isolated from soil in southeastern Asia [91]; 2) a C-F subclade that includes Scheffersomyces sp. UFMG-HMD-26.3 (not yet validly published [20]), S. coipomoensis, S. ergatensis, S. lignicola, and S. queiroziae; and 3) the X-F subclade containing S. lignosus, S. illinoinensis, S. insectosa, S. quercinus, S. segobiensis, S. stipitis, S. shehatae, and S. virginianus (Fig. 2).

The phylogenetic outcome also suggests that the common ancestor of the Scheffersomyces clade may have shown the ability to ferment D-xylose and cellobiose. This hypothesis is supported by the phylogenetic analyses based on a multilocus dataset that places S. gosingicus, a C-F species, in the earliest derived subclade, and the results of the phylogenetic study showing XR presence in all Scheffersomyces clade members studied (Fig. 5). Moreover, the same tree topology of the multilocus analysis was recovered by using either XYL1 or XR, suggesting that X-F ability might have played a fundamental role in the speciation process of the X-F subclade (Figs. 2, 5). We did not include the single copy gene XYL1 in the multilocus phylogenetic analysis because the orthology of this locus has not been confirmed across several yeast taxa.

The ability to ferment both wood components (D-xylose and cellobiose) is exhibited by only a few yeasts: Ogataea polymorpha [92], Brettanomyces naardenensis [1], and Spathaspora passalidarum [18], and in these species fermentative abilities are weak or delayed. On the contrary, Scheffersomyces clade members exhibit only one or the other fermentative ability. The loss of fermentation capability could be a consequence of becoming more efficient in carrying out the fermentation of fewer sugars. In particular, the fermentation of cellobiose has an antagonistic effect against the fermentation of D-xylose, because during the extracellular fermentation of cellobiose, units of glucose are released by β-glucosidase, a sugar that has higher affinity for the pentose membrane transporter rather than D-xylose. Consequently, glucose is first incorporated into the cells to be fermented [93].

XR in the Scheffersomyces Clade Members

The metabolic constraint in the regeneration of the cofactor NAD+ to NADH+ has been described as a major constraint in the fermentation of D-xylose, therefore most studies have focused on the molecular characterization of the conserved domains involved in the uptake of this cofactor that is present in the C-terminal region of the XR, but few studies have been done on the characterization of the N-terminal region.

In the study in the N-terminal region of the XR of the X-F members of the S. stipitis subclade we were able to identify all of the conserved domains and amino acids in both X-F and non X-F yeasts. These findings indicated that the ability to ferment xylose does not rely solely on the presence of these conserved regions. We also found several biased nucleotide mutations that maintain the same polarity as the codified amino acid, and only the mutations on the residues 174, 179, 184, 188, and 204 showed a change in amino acid polarity. These mutations were found mainly surrounding conserved domains in comparison with the amino acid sequence of the non X-F yeast S. ergatensis (Fig. 6). The mutations could generate structural modifications that allow the fermentation of xylose in the Scheffersomyces clade members, results that could be supported by performing direct mutagenesis studies on the amino acids to characterize their individual roles in the performance of XR.

Although the X-F yeasts were dispersed throughout the Saccharomycotina, as we mentioned in the introduction, the Debaryomycetaceae includes the largest number of X-F yeasts, such as species of Scheffersomyces and Spathaspora found in association with lignicolous insects [21], [24]. This phylogenetic placement of the X-F yeasts supports two alternatives: the X-F ability was the result of convergent evolution in ascomycete yeasts, or the X-F ability was present in the earliest common ancestor in the Saccharomycotina and has been retained mainly in the yeasts associated with lignicolous habitats. Several independent lines of evidence favor the second premise: 1) classical biochemical studies have determined that the ability to assimilate D-xylose and xylitol is common throughout the yeasts; 2) several relatively early diverging yeasts, B. naardenensis, O. polymorpha, and P. tannophilus, ferment xylose [1], [92]; 3) many studies characterizing yeast diversity indicate that most X-F yeasts are associated with wood and the gut of lignicolous insects [21]; and 4) more recently, a report of the genome sequences of a diverse group of X-F and non X-F yeasts confirmed the presence of xylose genes in all of them [45].

Yeasts as a group, are known for their biochemical versatility in utilizing a wide variety of carbon sources. Individual strains, however, may be characterized by specific physiological profiles depending on their life style and environment. Wood substrates and the gut of lignicolous insects previously were unexplored environments for the isolation of X-F yeasts. The findings of this study further support the hypothesis that X-F yeasts and yeasts in the Sugiyamaella clade are common inhabitants of the wood substrates. The Scheffersomyces clade is comprised mainly of cellobiose- and D-xylose-fermenting yeasts isolated from distant geographical regions and associated with wood and insects that feed on plant tissues. The amino acid modifications present in the putative XR of X-F yeasts in the S. stipitis subclade, could be responsible for the enhanced rate of fermentation shown by the members of this clade. The addition of ITS and RPB1 loci in the phylogenetic studies on the Scheffersomyces clade dramatically increased the support of the phylogenetic relationships of the members. We have used the primers for XYL1 designed for this study across several yeast species (data not show), and they could be used to help understand how these genes have evolved in the members of Saccharomycotina. The phylogenetic reconstruction using only XYL1 or RPB1 was similar to the multilocus analysis, and these loci have potential for rapid identification of cryptic species in this clade.