M&E Papers in Press. Published online on July 13, 2010 doi:10.1264/jsme2.ME10108

1 Revised ME10108

2

3 Molecular Characterization of Fungal Communities in Non-Tilled, Cover-Cropped

4 Upland Rice Field Soils

5

1†‡* 1† 1, 2 6 TOMOYASU NISHIZAWA , ZHAORIGETU , MASAKAZU KOMATSUZAKI , YOSHINORI

2 3 1, 2 7 SATO , NOBUHIRO KANEKO , and HIROYUKI OHTA

8

9 1Ibaraki University College of Agriculture, 3-21-1 Chuou, Ami-machi, Ibaraki 300-0393, Japan;

10 2Institute for Global Change Adaptation Science, Ibaraki University, 2-1-1 Bunkyo, Mito,

11 Ibaraki 310-8512, Japan; and 3Graduate School of Environment and Information Science,

12 Yokohama National University, 79-7 Tokiwadai, Hdogaya-ku, Yokohama,Proofs Kanagawa 240-8501,

13 Japan 14 View 15 (Received February 4, 2010 – Accepted June 16, 2010)

16

17 * Corresponding author. E-mail: [email protected]; Tel: +81–29–888–8684; Fax: 18 +81–29–888–8525.Advance 19 † T.N. and Z. contributed equally to this study.

20 ‡ Present address: Department of Applied Biological Chemistry, Graduate School of Agricultural

21 and Life Sciences, The University of Tokyo, 1–1–1, Yayoi, Bunkyo-ku, Tokyo 113–8657, Japan

22

23 Running headline: Soil Fungal Community in Upland Rice Field

24

25

1

Copyright 2010 by the Japanese Society of Microbial Ecology / the Japanese Society of Soil Microbiology 26 Abstract

27 This study aimed to characterize soil fungal communities in upland rice fields

28 managed with tillage/non-tillage and winter cover-cropping (hairy vetch and cereal rye)

29 practices, using PCR-based molecular methods. The study plots were maintained as

30 upland fields for 5 years and the soils sampled in the second and fifth years were

31 analyzed using T-RFLP (terminal restriction fragment length polymorphism) profiling

32 and clone libraries with the internal transcribed spacer (ITS) region and domain 1 (D1)

33 of the fungal large-subunit (fLSU) rRNA (D1fLSU) as the target DNA sequence. From

34 the 2nd-year-sample, 372 cloned sequences of fungal ITS-D1fLSU were obtained and

35 clustered into 80 nonredundant fungal OTUs (operational taxonomic units) in 4 fungal 36 phyla. The T-RFLP profiling was performed with the 2nd- andProofs 5th-year-samples and the 37 major T-RFs (terminal restriction fragments) were identified using a theoretical

38 fragment analysis of the ITS-D1fLSU clones. These molecular analyses showed that the

39 fungal community was influenced more stronglyView by the cover-cropping than tillage

40 practices. Moreover, the non-tilled, cover-cropped soil was characterized by a

41 predominance of sp. in the phylum . We provided a genetic

42 database of the fungal ITS-D1fLSUs in the differently managed soils of upland rice fields.

43 Advance

44 Key words: clone library, Cryptococcus, soil fungal community, T-RFLP, ribosomal

45 RNA gene

46

47

2

Copyright 2010 by the Japanese Society of Microbial Ecology / the Japanese Society of Soil Microbiology 48 Introduction

49 Agricultural management practices of non-tillage and cover-cropping have a great

50 influence on soil carbon and nitrogen dynamics through soil microbial activity (15, 40).

51 In most well-aerated, upland fields, fungi account for the majority of the total microbial

52 population (3). In addition to the saprophytic property, evidence of fungal nitrous

53 oxide-producing denitrification has been revealed in grassland soil using a combination

54 of the substrate-induced respiration inhibition method and labeling with a stable isotope

55 of nitrogen (18). Chemical inhibitors of fungi and bacteria decreased the flux of nitrous

56 oxide by 89% and 23%, respectively, suggesting that the fungal community was

57 involved in the nitrous oxide production. Prior to this finding, experiments in vivo had 58 shown that nitrous oxide was the main end product of fungal denitrificationProofs (30, 36). 59 Recently, Zhaorigetu et al. (40) reported a positive relationship between fungal biomass

60 and nitrous oxide emissions in an upland rice field. It was also shown that the fungal

61 biomass in the surface soil (0-10 cm) incrVieweased with the adoption of conservation

62 management systems such as non-tilled fallow and cover-cropping (hairy vetch and rye)

63 practices. However, the soil fungal community associated with the production of nitrous

64 oxide still remains indistinct in upland rice field soil.

65 To assessAdvance microbial diversify, small-s ubunit (SSU) and/or large-subunit (LSU) rRNA

66 genes are often used as markers (1). Sequencing of either LSU-rRNA or internal

67 transcribed spacer (ITS) regions has been used for the identification of fungal species

68 (12, 28, 29, 32). Molecular methods such as fungal rRNA gene and ITS region-based

69 clone library and phylogenetic analyses were utilized for identifying the dominant fungi

70 and elucidating the diversity of fungal communities in environmental samples (2, 24, 25,

71 27, 37). Terminal-restriction fragment length polymorphism (T-RFLP) profiling (19, 21)

3

Copyright 2010 by the Japanese Society of Microbial Ecology / the Japanese Society of Soil Microbiology 72 using an automated DNA sequencer with fluorescence-labeling nucleotides is a

73 powerful tool for monitoring the diversity in microbial communities (5, 7, 10, 11, 14, 23,

74 27, 31).

75 The objective of this study was to characterize the fungal community of nitrogen

76 oxide-production-enhanced upland soils maintained with cover-cropped and non-tillage

77 practices in comparison with those of conventional tillage soils. PCR-based T-RFLP

78 profiling and a clone library based on the fungal ITS–LSU-rRNA region in upland rice

79 fields were used.

80

81 Materials and Methods 82 Field management and soil sampling Proofs 83 The study plots were located at the Field Science Center of Ibaraki University

84 College of Agriculture, Ibaraki, Japan. The soil at the sites was a humic allophane soil

85 (Haludands, Haplic Andosols). The experimeViewnt was conducted over a period of 5 years,

86 from October 2003 until October 2007. Hairy vetch (Vicia villosa Roth.) and cereal rye

87 (Secale cereale L.) were planted every October from 2003 to 2007. Seeds of field rice

88 (Oryza sativa L.) were planted with a drill in April of each year. A detailed explanation

89 of the Advanceexperimental plots has been given previously (23, 40). In the present study, four

90 soils were used: plow-tilled fallow (PF), non-tilled fallow (NF), non-tilled hairy

91 vetch-cover cropped (NH), and non-tilled rye-cover cropped (NR) soil. Soil samples

92 were taken on May 14, 2004 and May 15, 2007. Five sub-samples (10 g fresh weight

93 per sub-sample) per depth of between 0 and 10 cm were collected from three plots using

94 sterile implements and then thoroughly mixed. The 50 g of soil was kept at -20˚C in

95 darkness before the genetic analysis.

4

Copyright 2010 by the Japanese Society of Microbial Ecology / the Japanese Society of Soil Microbiology 96

97 DNA extraction and manipulation

98 Total DNA from soil samples (0.5–1.0 g wet weight) was isolated using an Isoil for

99 Beads Beating kit (Nippon Gene, Tokyo, Japan) with skim milk powder (Wako, Osaka,

100 Japan) as described previously (23). Escherichia coli DH5α (Takara Bio, Otsu, Japan)

101 was used as a host for recombinant plasmids and grown at 37˚C in Luria-Bertani’s (LB)

102 broth or LB agar. As an antibiotic, ampicillin sodium salt was added at a final

103 concentration of 75 µg mL-1 when necessary. The pGEM-T-Easy vector system

104 (Promega, Madison, WI, USA) was used for cloning according to instructions.

105 106 PCR amplification of the internal transcribed spacer (ITS) sequenceProofs and domain 1 (D1) 107 of the large-subunit (LSU) rRNA gene

108 The fungal ITS-D1fLSU rRNA region in soil DNA (0.1 µg) was amplified using ITS1

109 (5’- TCCGTAGGTGAACCTGCGG -3’)View and LR21 (5’- ACTTCAAGCGTTTCCCTTT

110 -3’, position 424-393 in Saccharomyces cerevisiae large subunit rRNA). The reactions

111 were performed in a Takara PCR thermal cycler Personal (Takara Bio) with Ex Taq 112 (TakaraAdvance Bio) under the following conditions: 2 min at 95˚C, followed by 30 cycles of 113 95˚C (30 s), 54˚C (45 s), and 72˚C (1.5 min), with a final cycle of 72˚C for 3 min.

114 For the construction of the clone library, the ITS-D1fLSU rRNA gene region amplified

115 using DNA of PF, NF, NH, and NR soils sampled in 2004 was cloned and the PCR

116 DNA fragments were ligated into the pGEM-T-easy vector.

117

118 Sequencing of the ITS region and computer analysis of DNA sequences

119 Nucleotide sequences were determined with a BigDye Terminator Cycle Sequencing

5

Copyright 2010 by the Japanese Society of Microbial Ecology / the Japanese Society of Soil Microbiology 120 Ready Reaction Kit (Applied Biosystems, Foster City, CA, USA) according to the

121 instructions supplied and with a PE Applied Biosystems Automated DNA Sequencer

122 (model 3130xl, Applied Biosystems). The double-stranded DNA sequences were

123 assembled and analyzed using Genetyx (version 11.0) and Genetyx-ATSQ (version 4.0)

124 software (Genetyx, Tokyo, Japan), respectively, and compared with similar DNA

125 sequences retrieved from the DDBJ/EMBL/GenBank databases using the

126 NCBI-BLAST and DDBJ-BLAST programs. The phylogenetic analysis was carried out

127 using the multiple-sequence alignment tool from Clustal W (35).

128

129 Terminal-restriction fragment length polymorphism (T-RFLP)-based profiling 130 The ITS regions were amplified by PCR from the total genomicProofs soil DNA with the 131 primers QLR21 (5’- CACTTCAAGCGTTTCCCTTT -3’), labeled with quenching

132 fluorescence, and ITS1. The 5’-end fluorescence-labeled primer (16) was purchased

133 from J–Bio21 (Tsukuba, Japan). The PCRView mixture (30 µL) was prepared by combining

134 0.1 µg of template soil DNA, 1.0 µL of 10 pmol µL-1 primers, Takara Ex Taq, dNTPs,

135 and 3 µL of optimized 10×Ex buffer (Takara Bio Inc.) in a PCR cycler. The PCR of ITS 136 for T-RFLPAdvance profiling was carried out under the following conditions: 2 min at 95˚C, 137 followed by 30 cycles of 95˚C (30 s), 54˚C (45 s), and 72˚C (1.5 min).

138 T-RFLP was carried out as described previously (23). Aliquots (5 µL) of the

139 amplified ITS-D1fLSU fragment were separately digested with AluI, HaeIII, MboI, MspI,

140 and Sau3AI (Takara Bio) according to the manufacturer’s instructions. The precise

141 lengths of T-RFs (terminal restriction fragments) from the amplified ITS-D1fLSU

142 fragments were determined on a 3130xl PE DNA Sequencer (Applied Biosystems). The

143 purified T-RF DNA (2 µL) was mixed with 15 µL of Hi-Di formamide and 0.1 µL of

6

Copyright 2010 by the Japanese Society of Microbial Ecology / the Japanese Society of Soil Microbiology 144 DNA standard LIZ®500 (Applied Biosystems) for standardization. This mixture was

145 then denatured at 96˚C for 2 min and immediately chilled on ice prior to electrophoresis

146 with the automated DNA sequencer in the GeneScan mode. The lengths of fluorescently

147 labeled T-RFs were determined by comparison with internal standards using

148 GeneMapper software (version 3.7; Applied Biosystems). T-RFLP profiling of the

149 fungal community structure in soil samples was obtained with peaks ranging from 50

150 bases to 650 bases. The definition of T-RFs was peaks with a fluorescence threshold of

151 more than 30.

152

153 Analysis of the diversity of data sets 154 Diversity was analyzed using different definitions of operationalProofs taxonomic units 155 (OTUs). The conventional affiliation of a fungal species (OTU) was based o the

156 ITS-D1fLSU rRNA sequence within the GenBank database. Due to the computational

157 challenges associated with estimating diversViewity indices in a previous study (23), the

158 observed richness and diversity were calculated with the Chao1 (6) and

159 Shannon-Wiener indices utilizing the web-based bioinformatics platform, FastGroupII

160 (39) as a group with 97% identity in ITS-D1fLSU rRNA region sequence, respectively.

161 StatisticalAdvance calculations were made with Microsoft Excel using the two-sided, unpaired

162 Student's t-test. The similarity of data sets (relative abundance of the peak height and

163 fragment length) of T-RFs profiles was examined using nonmetric multidimensional

164 scaling (NMDS) based on the free software package ‘R’ (version 2.8.0) (26).

165

166

7

Copyright 2010 by the Japanese Society of Microbial Ecology / the Japanese Society of Soil Microbiology 167 Sequence accession numbers

168 All cloned sequences were deposited in the DNA Database of Japan (DDBJ) under

169 accession numbers AB520238 to AB520609.

170

171 Results

172 Cloning of the ITS and D1 domain of the large subunit rRNA gene

173 To investigate the fungal community structure, the primer pair ITS1/LR21 was used

174 to amplify a long DNA fragment (850-950 bp) comprising the whole ITS region and

175 approximately 300 bp of the LSU rRNA gene containing the domain 1 (D1)

176 region (D1fLSU) (Fig. 1). Simultaneous amplification of both the conserved and variable 177 regions makes it possible to infer the phylogenetic position ofProofs the fungus at a higher 178 taxonomic level. Tedersoo et al. (33) described that additional sequencing of the LSU

179 rRNA genes could improve the determination of some unknown species, at least at the

180 family or order level. View

181 We observed an unexpected irregular peak in the T-RFLP profile when total DNA of

182 the fungal wild-type strain was used for amplification with the fluorescence

183 labeled-ITS1F/ITS4 (9, 38) (data not shown). We therefore decided to utilize the

184 ITS1/LR21Advance as a stable primer set to avoid irregular peaks as a result of amplification

185 with the ITS4 primer (Fig. 1).

186 The fungal diversity in samples from four differently managed soils was evaluated by

187 including the ITS and D1fLSU (ITS- D1fLSU) libraries for each soil type. PCR

188 amplifications were successfully performed for samples of four differently managed

189 soils. A total of 372 clones (PF, 86 clones; NF, 98; NH, 93; NR, 95) were sequenced. To

190 identify the closest match of cloned sequences, an NCBI-BLAST search was carried out

8

Copyright 2010 by the Japanese Society of Microbial Ecology / the Japanese Society of Soil Microbiology 191 and the fungal ITS-D1fLSU database shown in Table S1 was identified. The taxonomic

192 analysis revealed that most sequences could be assigned to 4 phyla or

193 unclassified/uncultured fungal clones of soil samples that were unevenly distributed

194 throughout the PF, NF, NH, and NR plots. According to this result, the ITS-D1fLSU

195 rRNA region was suitable for characterization of the assemblage of soil fungi. From the

196 appearance rates of clones, it was suggested that the diversity of fungi in the soil was

197 dependent on the soil management practices (Table 1). In PF soil, 54.7% of the cloned

198 sequences were identified in unculturable fungal clones. Ascomycota accounted for

199 26.7% of the total PF clone library. In NF soil, unculturable fungal clones and

200 Basidiomycota cloned sequences accounted for 45.9% and 33.7% of the total, 201 respectively. On the other hand, the NH and NR soil samples Proofstaken from plots with 202 cover-cropping management consisted primarily of Basidiomycota (NH, 49.4%; NR,

203 68.4%). Moreover, in NH soil, Ascomycota accounted for 34.4% of the clones. The

204 uncultured clones (IU-FSC Fun71 and Fun66)View were found in the four libraries. The ratio

205 of Basidiomycota was increased by the non-tilled fallow management. In particular,

206 sequences of the non-mycorrhizal fungal species, the Cryptococcus sp. (OTUs: IU-FSC

207 Fun22 to Fun28) of the phylum Basidiomycota, were predominant in the clone libraries

208 of the AdvanceNF, NH, and NR soils (Table S1). It was also confirmed that Basidiomycota was

209 predominant in non-tilled, cover-cropped soils and uncultured fungal clones dominated

210 in the fallow-treated soil.

211

212 Diversity of identified sequences

213 In order to analyze the fungal diversity associated with the different soil management

214 practices, the proportion of eukaryotic microbial diversity represented by the clone

9

Copyright 2010 by the Japanese Society of Microbial Ecology / the Japanese Society of Soil Microbiology 215 libraries was estimated by rarefaction analysis combined with nonlinear regression. The

216 diversity of the ITS-D1fLSU sequences was analyzed using the FastGroupII

217 bioinformatics tool (Fig. 2). Although the rarefaction curves from all soils failed to

218 plateau, the rarefaction analyses showed that the NF soil had greater fungal richness

219 than the other soils (Fig. 2). In contrast to the results above for the soils subjected to

220 different management practices, the comparative rarefaction results were corroborated

221 by nonparametric-based estimates of total richness and diversity. The Chao1 total

222 richness estimate for the PF soil was 63 ribotypes (OTUs). The Chao1 estimate for the

223 NF soil showed the highest diversity with 87 OTUs, whereas the NH and NR soil

224 samples had total richness values of 65 and 81 OTUs, respectively. The coverage of the 225 number of OTUs analyzed relative to the expected maximum Proofsnumber of ribotypes for 226 ITS-D1fLSU ranged from 38% to 57%. Additionally, the Shannon-Wiener index, a

227 measurement of biodiversity used to estimate the number of species and their

228 homogeneity across a given plot values, wereView very similar for the four differently

229 managed soils (PF, 2.9; NH, 2.6; NR, 2.7), with the exception of 3.6 with the NF library.

230 These molecular analyses clarified that non-tilled fallow soil had remarkable fungal

231 diversity.

232 Advance

233 Identification of the predominant T-RFs combined with in silico T-RFs

234 To investigate the fungal community in the soils managed with non-tillage and

235 cover-cropping practices in the upland field, fungal ITS-D1fLSU PCR-based T-RFLP

236 analyses were carried out (Fig. S1). It seemed that the fungal T-RFLP profiles had their

237 own characteristic pattern of restriction enzyme-digested ITS-D1fLSUs. The size of the

238 experimental T-RFs (T-RFexps) on T-RFLP profiling was determined by the in silico

10

Copyright 2010 by the Japanese Society of Microbial Ecology / the Japanese Society of Soil Microbiology 239 T-RF of the fungal ITS-D1fLSU cloned sequence (T-RFseqs). Furthermore, to determine

240 the dominant T-RFexps, we used the database of T-RFseqs shown in Table S1, which was

241 determined by an analysis of the clone libraries. This database included OTUs chosen

242 from more than 5 cloned sequences from the PF, NF, NH, and NR soils (Table S1).

243 In our previous study, a terminal restriction fragment length method utilizing a

244 quenching fluorescence primer was developed and its accuracy was confirmed (23). The

245 112-base and 143-base AluI T-RFexps were contained in T-RF profiles of four differently

246 managed soils (data not shown). This result was consistent with the results for in silico

247 AluI T-RFseq (Table S1), indicating that AluI is not suitable for comparing soil fungal

248 communities. Four HaeIII T-RFexps [102-, 258-, 298-, and 445-base] were mainly 249 observed in the four soils sampled in 2004 (Fig. S1A). The dominantProofs HaeIII 102-base 250 T-RF exp was shown to correspond precisely to the IU-FSC of the phylum Ascomycota 251 and IU-FSC Fun63 and Fun70 of unculturedView fungus (Table S1). The HaeIII 258-base 252 and 298-base T-RFexps were consistent with the HaeIII 259-base and 299-base T-RFseqs

253 of IU-FSC Fun29 and Fun39, respectively. The HaeIII 445-base T-RFexp corresponded

254 to the HaeIII 447-base T-RFseq of IU-FSC Fun66b. In the MboI/Sau3AI-digested

255 T-RFLP profiles, the MboI/Sau3AI 68/69-base and 425/426-base T-RFexps were

256 identifiedAdvance in the 04PF (PF soil sampled in 2004), 04NF, and 04NH soils, respectively

257 (Fig. S1B). According to the clone library analysis, the dominant MboI/Sau3AI

258 68/69-base T-RFexp corresponded to the 72-base T-RFseq of IU-FSC Fun04, Fun07,

259 Fun09, and Fun13 among the phylum Ascomycota and IU-FSC Fun63 and Fun70 of

260 uncultured fungal clones. The dominant MboI/Sau3AI 425/426-base T-RFexp was shown

261 to correspond to the T-RFseq of IU-FSC Fun59, Fun66ab, and Fun71. The dominant

262 MspI 103-base T-RFexps were identified as the 104-base T-RFseq of IU-FSC Fun07 and

11

Copyright 2010 by the Japanese Society of Microbial Ecology / the Japanese Society of Soil Microbiology 263 Fun66c (Fig S1C and Table S1). The MspI 197-base T-RFexps was consistent with the

264 199-base T-RFseq of IU-FSC Fun70. In the case of 04NH soil, two T-RFexps [190- and

265 505/507-base] were observed and similar to the T-RFseq of IU-FSC Fun24, Fun25a,

266 Fun28b, and Fun38a and IU-FSC Fun04 and Fun09ab, respectively (Fig. S1C and Table

267 S1). According to these results, the in silico analysis of clone libraries enabled us to

268 determine the major T-RFs in the fungal T-RFLP profiling.

269

270 Comparisons of T-RFLP profiling

271 The soil fungal community structure was assessed in the soils treated with different

272 management practices between 2004 and 2007 based on the identification of major 273 T-RF exps using the 2004 data from T-RFLP profiling. The HaeProofsIII, MboI/Sau3AI and 274 MspI-digested T-RFLP profiles showed the change in diversity of the fungal community

275 among the different soil conditions from 2004 to 2007 as judged with the two-sided,

276 unpaired Student's t-test (P≤0.05) (Table View2). Continuous non-tillage fallow rye

277 cover-cropping decreased the number of T-RFexps in the soil fungal community

278 compared to non-tillage fallow hairy vetch cover-cropping.

279 In the HaeIII T-RFLP profiling, the fungal community apparently differed in T-RFs of

280 aroundAdvance 400- and 600-bases between 04NH and 07NH (NH soil sampled in 2007) (Fig.

281 S1A). In the detailed comparison of MboI/Sau3AI T-RFLP profiles, two MboI/Sau3AI

282 T-RF exps [356/357- and 425-base] mainly appeared in the T-RFLP profiles of the 07PF,

283 07NF, 07NH, and 07NR soils. Likewise, the dominant MboI/Sau3AI 428/430-base

284 T-RF exp appeared in the 04NR soil. The major T-RFexp cluster consisted of around 420

285 bases in the 2004 soil. However, another T-RFexp cluster of around 350 bases occurred

286 in the 2007 soil. A change in the fungal community structure was very likely from the

12

Copyright 2010 by the Japanese Society of Microbial Ecology / the Japanese Society of Soil Microbiology 287 comparison of the MboI/Sau3AI-digested T-RFexps profiles with the 04NF (Fig. S1B). In

288 the MspI T-RFLP profiling, around 190- and 510-base MspI T-RFexps were clustered in

289 the 04NH and 04 NR soils. In particular, a convergence of 190-base T-RFexps occurred

290 in the 04NH and 04NR soils (Fig. S1C and Table 2). The occurrence of 190-base

291 MspI-T-RFexp remained until 2007.

292

293 Discussion

294 In the present study, the fungal T-RFLP profiles showed a characteristic pattern of

295 restriction enzyme-digested ITS-D1fLSUs and the ITS-D1fLSU rRNA-based clone library

296 analysis elucidated the predominant T-RFs in T-RFLP profiling. Additionally, a 297 comparative analysis using the data of T-RF profiles in samplesProofs of 2004 and 2007 298 revealed that continuous soil management practices resulted in a change in diversity of

299 the fungal community (Fig. S1 and Table 2). To estimate the similarity in fungal

300 communities among the differently managedView soils, NMDS was performed based on the

301 results of T-RFLP profiling and clone libraries. First of all, a HaeIII-, MboI/Sau3AI-, or

302 MspI-T-RFLP profiling based-NMDS was carried out (data not shown). NMDS of the

303 MboI/Sau3AI-T-RFLP profiling showed that the four fungal communities were highly

304 similar,Advance indicating that it was quite difficult to understand the similarity or dissimilarity

305 of fungal community. Therefore, NMDS analysis of the combined HaeIII- and

-3 306 MspI-T-RFexp profiles was performed (stress value, 8.0×10 ) (Fig. 3A). NMDS

307 combined with T-RFexps profiles showed that the effect of tillage and non-tillage

308 management practices was minor because of the close distance among the profiles of PF

309 and NF from 2004 to 2007. Regarding the NR soil, a significant effect of rye

310 cover-cropping practice on the soil fungal community structure was clearly confirmed

13

Copyright 2010 by the Japanese Society of Microbial Ecology / the Japanese Society of Soil Microbiology 311 by NMDS (Fig. 3A). Recently, it was reported that the degradation (for 5.5 months) of

312 wheat straw residue in soil had almost no influence on fungal community structure

313 according to a statistical analysis of fungal ITS-based DNA fingerprints (4). Likewise,

314 as for the clone library, T-RFseq data (fragment lengths and ratio of clone number) after

315 digestion by HaeIII and MspI in silico were used for NMDS. A similar fungal

316 community structure between hairy vetch- and rye-cover cropped soils appeared in the

317 NMDS analysis (stress value, 1.1×10-14; Fig. 3B). Our results indicated the cover crop

318 residue, especially rye straw residue, on non-tilled fallow soil might strongly influence

319 the fungal community at the surface of soil.

320 A total of 152 cloned sequences of the phylum Basidiomycota were found in the four 321 clone libraries of soil sampled in 2004. Some 95 clones were Proofsrelated to the genus 322 Cryptococcus, indicating that OTUs belonging to the phylum Basidiomycota,

323 Cryptococcus predominated in the non-tilled soils. A detailed search of IU-FSC

324 Fun22-Fun29 was carried out using the GenBaViewnk database. Eight OTUs were divided

325 into three clades: 1) IU-FSC Fun24, Fun26, Fun27, and Fun28 were Basidiomycota,

326 Agaricomycotina, , , mitosporic Tremellales, Cryptococcus; 327 2) IU-FSCAdvance Fun23 and Fun29 were Basidiomycota, Agaricomycotina, Tremellomycetes, 328 Cystofilobasidiales, mitosporic Cystofilobasidiales, Cryptococcus; 3) IU-FSC Fun22

329 and Fun25 were Basidiomycota, Agaricomycotina, Tremellomycetes, Filobasidiales,

330 mitosporic Filobasidiales, Cryptococcus; respectively. In particular, IU-FSC Fun28 (41

331 clones), with a similar sequence to Cryptococcus victoriae (34), was dominant in the

332 non-tilled fallow and non-tilled cover-cropping soil. Based on the T-RFLP profiling of

333 NR soils sampled in 2004 and 2007, the characteristic pattern of HaeIII-MboI-MspI

334 [(428/430)-(429/430)-(190/191)] T-RFexps corresponding to IU-FSC Fun28b (40 clones)

14

Copyright 2010 by the Japanese Society of Microbial Ecology / the Japanese Society of Soil Microbiology 335 was observed (Fig. S1 and Table S1). These results indicate that the contribution of

336 cover-cropping of rye enhanced this trend. The clade Tremellales is well-known to be

337 made up of human pathogenic strains such as C. neoformans, C. albidus, and C.

338 laurentii (8, 32). Therefore, a phylogenetic tree was constructed among Tremellales, and

339 Cryptococcus containing the OTUs IU-FSC Fun24, 26, 27, and 28 (Fig. 4).

340 Phylogenetic analysis based on the neighbor-joining method using CLUSTAL W

341 showed that three of the four OTUs (i.e., all except for IU-FSC Fun24) were integrated

342 into one group. Obviously, these clones were distinguished from the human pathogenic

343 Cryptococcus strains (Fig. 4). The clones of the putative Cryptococcus sp. strains

344 (IU-FSC Fun24, 26, 27, and 28) were almost identical with those of Cryptococcus sp. 345 from soils (13, 17, 20, 22, 27). This suggests that clones of soilProofs Cryptococcus sp. are 346 widespread throughout the world. However, the role of Cryptococcus in soil has not

347 been clarified yet. With respect to nitrous oxide-producing activity in Cryptococcus

348 strains of Basidiomycota, Cryptococcus gastricusView JCM 3691 in the clade Filobasidiales

349 did not produce nitrous oxide or nitric oxide in pure cultures (36). The isolation of

350 Cryptococcus strains from the NR soil is now in progress in our laboratory.

351 In conclusion, our results showed the diversity of soil-surface fungal communities in

352 uplandAdvance rice field by a combination of T-RFLP profiling and clone library analysis based

353 on fungal ITS-D1fLSU sequences and revealed that soil management techniques such as

354 non-tilling and cover-cropping were involved in determining the diversity and dynamics

355 of the soil-surface fungal community. Furthermore, we provided a genetic database of

356 the fungal ITS-D1fLSUs in the different agricultural soils of the upland rice field.

357

358

15

Copyright 2010 by the Japanese Society of Microbial Ecology / the Japanese Society of Soil Microbiology 359 Acknowledgements

360 This work was supported by a grant from the Global Environment Research Fund of the

361 Ministry of the Environment, Japan (F-073) and in part by MEXT through Special

362 Coordination Funds for Promoting Science and Technology as a part of the “Sustainable

363 Agriculture Practices to Mitigate and Adapt to Global Warming” project undertaken by

364 the Institute for Global Change Adaptation Science, Ibaraki University.

365

366 References

367 1. Amann, R., and W., Ludwig. 2000. Ribosomal RNA-targeted nucleic acid probes

368 for studies in microbial ecology. FEMS Microbiol. Rev. 24:555–565. 369 2. Anderson. I.C., B.A. Bastias, D.R. Genney, P.I. Parkin, andProofs J.W.G. Cairney. 2007. 370 Basidiomycete fungal communities in Australian sclerophyll forest soil are altered

371 by repeated prescribed burning. Mycol. Res. 111:482–486.

372 3. Anderson, J.P.E., and K.H. Domsch. View1975. Measurement of bacterial and fungal

373 contributions to respiration of selected agricultural and forest soils. Can. J.

374 Microbiol. 21:314–322.

375 4. Bastian, F., L. Bouziri, B. Nicolardot, and L. Ranjard. 2009. Impact of wheat straw

376 decompositionAdvance on successional patterns of soil microbial community structure. Soil

377 Biol. Biochem. 41:262–275.

378 5. Brodie, E., S. Edwards, and N. Clipson. 2003. Soil fungal community structure in a

379 temperate upland grassland soil. FEMS Microbiol. Ecol. 45:105–114.

380 6. Chao, A. 1984. Non-parametric estimation of the number of classes in a population.

381 Scand. J. Statist. 11:783–791.

382 7. Euringer, K., and T. Lueders. 2008. An optimised PCR/T-RFLP fingerprinting

16

Copyright 2010 by the Japanese Society of Microbial Ecology / the Japanese Society of Soil Microbiology 383 approach for the investigation of protistan communities in groundwater

384 environments. J. Microbiol. Methods 75:262–268.

385 8. Fonseca, A., G. Scorzetti, and J.W. Fell. 2000. Diversity in the yeast Cryptococcus

386 albidus and related species as revealed by ribosomal DNA sequence analysis. Can.

387 J. Microbiol. 46:7–27.

388 9. Gardes, M., and T.D. Bruns. 1993. ITS primers with enhanced specificity for

389 basidiomycetes-application to the identification of mycorrhizae and rusts. Mol.

390 Ecol. 2:113–118.

391 10. Hartmann, M., and F. Widmer. 2006. Community structure analyses are more

392 sensitive to differences in soil bacterial communities than anonymous diversity 393 indices. Appl. Environ. Microbiol. 72:7804–7812. Proofs 394 11. Hartmann, M., and F. Widmer. 2008. Reliability for detecting composition and

395 changes of microbial communities by T-RFLP genetic profiling. FEMS Microbiol.

396 Ecol. 63:249–260. View

397 12. Horton, T.R., and T.D. Bruns. 2001. The molecular revolution in ectomycorrhizal

398 ecology: peeking into the black-box. Mol. Ecol. 10:1855–1871.

399 13. Hoshino, Y.T., and N. Matsumoto. 2007. Changes in fungal community structure in

400 bulkAdvance soil and spinach rhizosphere soil after chemical fumigation as revealed by 18S

401 rDNA PCR-DGGE. Soil Sci. Plant Nutr. 53:40–55.

402 14. Kasel, S., L.T. Bennett, and J. Tibbits. 2008. Land use influences soil fungal

403 community composition across central Victoria, south-eastern Australia. Soil Biol.

404 Biochem. 40:1724–1732.

405 15. Komatsuzaki, M., and H. Ohta. 2007. Soil management practices for sustainable

406 agro-ecosystems. Sustain. Sci. 2:103–120.

17

Copyright 2010 by the Japanese Society of Microbial Ecology / the Japanese Society of Soil Microbiology 407 16. Kurata, S., T. Kanagawa, K. Yamada, M. Torimura, T. Yokomaku, Y. Kamagata, and

408 R. Kurane. 2001. Fluorescent quenching-based quantitative detection of specific

409 DNA/RNA using a BODIPY® FL-labeled probe or primer. Nucleic Acids Res.

410 29:e34.

411 17. Landeweert, R., P. Leeflang, T.W. Kuyper, E. Hoffland, A. Rosling, K. Wernars, and

412 E. Smit. 2003. Molecular identification of ectomycorrhizal mycelium in soil

413 horizons. Appl. Environ. Microbiol. 69:327–333.

414 18. Laughlin, R.J., and R.J. Stevens. 2002. Evidence for fungal dominance of

415 denitrification and codenitrification in a grassland soil. Soil Sci. Soc. Am. J.

416 66:1540–1548. 417 19. Liu, W.-T., T.L. Marsh, H. Cheng, and L.J. Forney. 1997. ProofsCharacterization of 418 microbial diversity by determining terminal restriction fragment length

419 polymorphism of genes encoding 16S rRNA. Appl. Environ. Microbiol.

420 63:4516–4522. View

421 20. Lynch, M.D.J., and R.G. Thorn. 2006. Diversity of basidiomycetes in Michigan

422 agricultural soils. Appl. Environ. Microbiol. 72:7050–7056.

423 21. Marsh, T.L. 1999. Terminal restriction fragment length polymorphism (T-RFLP):

424 AnAdvance emerging method for characterizing diversity among homologous populations

425 of amplification products. Cur. Opin. Microbiol. 2:323–327.

426 22. Midgley, D.J., J.A. Saleeba, M.I. Stewart, A.E. Simpson, and P.A. McGee. 2007.

427 Molecular diversity of soil basidiomycete communities in northern-central New

428 South Wales, Australia. Mycol. Res. 111:370–378.

429 23. Nishizawa, T., M. Komatsuzaki, N. Kaneko, and H. Ohta. 2008. Archaeal diversity

430 of upland rice field soils assessed by the terminal restriction fragment length

18

Copyright 2010 by the Japanese Society of Microbial Ecology / the Japanese Society of Soil Microbiology 431 polymorphism method combined with real time quantitative-PCR and a alone

432 library analysis. Microbes Environ. 23:237–243.

433 24. O’Brien, H.E., J.L. Parrent, J.A. Jackson, J.-M. Moncalvo, and R.Vilgalys. 2005.

434 Fungal community analysis by large-scale sequencing of environmental samples.

435 Appl. Environ. Microbiol. 71:5544–5550.

436 25. Porter, T.M., C.W. Schadt, L. Rizvi, A.P. Martin, S.K. Schmidt, L. Scott-Denton, R.

437 Vilgalys, and J.M. Moncalvo. 2008. Widespread occurrence and phylogenetic

438 placement of a soil clone group adds a prominent new branch to the fungal tree of

439 life. Mol. Phylogenet. Evol. 46:635–644.

440 26. R Development Core Team. 2008. R: A language and environment for statistical 441 computing. R Foundation for Statistical Computing, Vienna,Proofs Austria. ISBN 442 3-900051-07-0, URL http://www.R-project.org.

443 27. Robinson, C.H., T.M. Szaro, A.D. Izzo, I.C. Anderson, P.I. Parkin, and T.D. Bruns.

444 2009. Spatial distribution of fungal communitieViews in a coastal grassland soil. Soil

445 Biol. Biochem. 41:414–416.

446 28. Schadt, C.W., A.P. Martin, D.A. Lipson, and S.K. Schmidt. 2003. Seasonal

447 dynamics of previously unknown fungal lineages in tundra soils. Science

448 301:1359–1361.Advance

449 29. Scorzetti, G., J.W. Fell, A. Fonseca, and A. Statzell-Tallman. 2003. Systematics of

450 basidiomycetous yeasts: A comparison of large subunit D1/D2 and internal

451 transcribed spacer rDNA regions. FEMS Yeast Res. 2:495–517.

452 30. Shoun, H., D.-H. Kim, H. Uchiyama, and J. Sugiyama. 1992. Denitrification by

453 fungi. FEMS Microbiol. Lett. 94:277–282.

454 31. Singh, B.K., L. Nazaries, S. Munro, I.C. Anderson, and C.D. Campbell. 2006. Use

19

Copyright 2010 by the Japanese Society of Microbial Ecology / the Japanese Society of Soil Microbiology 455 of multiplex terminal restriction fragment length polymorphism for rapid and

456 simultaneous analysis of different components of the soil microbial community.

457 Appl. Environ. Microbiol. 72:7278–7285.

458 32. Sugita, T., and A. Nishikawa. 2004. Molecular taxonomy and identification of

459 pathogenic fungi based on DNA sequence analysis. Jpn. J. Med. Mycol. 45:55–58

460 (in Japanese).

461 33. Tedersoo, L., U. Kõljalg, N. Hallenberg, and K.-H. Larsson. 2003. Fine scale

462 distribution of ectomycorrhizal fungi and roots across substrate layers including

463 coarse woody debris in a mixed forest. New Phytol. 159:153–165.

464 34. Thomas-Hall, S., K. Watson, and G. Scorzetti. 2002. Cryptococcus statzelliae sp. 465 nov. and three novel strains of Cryptococcus victoriae, yeastsProofs isolated form 466 Antarctic soils. Int. J. Syst. Evol. Microbiol. 52:2303–2308.

467 35. Thompson, J.D., D.G. Higgins, and T.J. Gibson. 1994. CLUSTAL W: Improving the

468 sensitivity of progressive multiple sequeViewnce alignment through sequence weighting,

469 position specific gap penalties and weight matrix choice. Nucleic Acids Res.

470 22:4673–4680.

471 36. Tsuruta, S., N. Takaya, L. Zhang, H. Shoun, K. Kimura, M. Hamamoto, and T.

472 Nakase.Advance 1998. Denitrification by yeasts and occurrence of cytochrome P450nor in

473 Trichosporon cutaneum. FEMS Microbiol. Lett. 168:105–110.

474 37. Weber, S.D., A. Hofmann, M. Pilhofer, G. Wanner, R. Agerer, W. Ludwig, K.-H.

475 Schleifer, and J. Fried. 2009. The diversity of fungi in aerobic sewage granules

476 assessed by 18S rRNA gene and ITS sequence analyses. FEMS Microbiol. Ecol.

477 68:246–254.

478 38. White, T.J., T.D. Bruns, S. Lee, and J. Taylor. 1990. Amplification and direct

20

Copyright 2010 by the Japanese Society of Microbial Ecology / the Japanese Society of Soil Microbiology 479 sequencing of fungal ribosomal RNA genes for phylogenetics, p.315–322. In T.J.

480 White, J.J. Sninsky, D.H. Gelfand, and M.A. Innin (ed.), PCR Protocols a Guide to

481 Methods and Applications Academic Press, San Diego, USA.

482 39. Yu, Y., M. Breitbart, P. McNairnie, and F. Rohwer. 2006. FastGroupII: A web-based

483 bioinformatics platform for analyses of large 16S rDNA libraries. BMC

484 Bioinformatics 7:57.

485 40. Zhaorigetu, M. Komatsuzaki, Y. Sato, and H. Ohta. 2008. Relationships between

486 fungal biomass and nitrous oxide emission in upland rice soils under no tillage and

487 cover-cropping systems. Microbes Environ. 23:201–208.

488 489 Proofs View

Advance

21

Copyright 2010 by the Japanese Society of Microbial Ecology / the Japanese Society of Soil Microbiology 490 Figure legends

491 Fig. 1. Representation of the primer position in a eukaryotic (yeast and filamentous

492 fungi) ribosomal RNA gene. The loci and direction of the primers are indicated by

493 arrows. ITS, internal transcribed spacer; LSU, large-subunit rRNA; SSU, small-subunit

494 rRNA. The plus (+) and minus (–) signs indicate a successful and unsuccessful analysis,

495 respectively. Primer positions of ITS1, LR21, QLR21, QITS1F (nucleotides sequence:

496 5’- CTTGGTCATTTAGAGGAAGTAA -3’) and ITS4 (5’-

497 TCCTCCGCTTATTGATATGC-3’) used for T-RFLP show in the ribosomal RNA gene

498 region.

499 500 Fig. 2. Rarefaction curves of the diversity of the ITS-D1fLSU regionProofs obtained for soil 501 management. Operational taxonomic units (OTUs) were assigned using a 3%

502 dissimilarity cut-off with sequence data of the ITS-D1fLSU region clone library. PF,

503 plow-tilled fallow; NF, non-tilled fallow; ViewNH, non-tilled fallow cover crop (hairy

504 vetch); NR, non-tilled fallow cover crop (rye).

505

506 Fig. 3. Nonmetric multidimensional scaling (NMDS) of terminal-restriction fragments

507 of T-RFLPAdvance profiling (A) and fungal clone libraries (B) from study sites. NMDS was

508 preformed using T-RFexps and T-RFseqs profiles combined HaeIII and MspI enzymes,

509 respectively. The two-dimensional stress values were 1.1×10-14 (A) and 8.0×10-3 (B),

510 respectively. PF, plow-tilled fallow (circle); NF, non-tilled fallow (triangle); NH,

511 non-tilled fallow cover crop of hairy vetch (rhombus); NR, non-tilled fallow cover crop

512 of rye (square). Years 2004 and 2007 are represented as 04 and 07, respectively. Closed

513 and open marks indicate the years 2004 and 2007, respectively. Dotted arrows indicate

22

Copyright 2010 by the Japanese Society of Microbial Ecology / the Japanese Society of Soil Microbiology 514 the movement of the soil fungal community structure.

515

516 Fig. 4. Phylogenetic relationship of the fungal internal transcribed spacer sequence of

517 the characterized clone library of non-tilled fallow managed soils. The method of

518 neighbor-joining was carried out using the package Clustal W (35). The sequences

519 obtained in the present study are shown in bold. Bootstrap values greater than 50 are

520 shown at branch points. Cystofilobasidium bisporidii CBS 6347 was used as an

521 out-group. The scale bar indicates 10% sequence divergence. The accession number

522 deposited in GenBank database is indicated in parentheses.

523 524 Proofs View

Advance

23

Copyright 2010 by the Japanese Society of Microbial Ecology / the Japanese Society of Soil Microbiology 525 Table 1. Appearance rate of sequenced clones in the different soil management 526 practices 527 Soil management methodsa) 528 Phylogenetic relations 04PF 04NF 04NH 04NR 529 Phylum 530 Ascomycota 26.7 % 4.1 % 34.4 % 8.4 % 531 Basidiomycota 10.5 % 33.7 % 49.5 % 68.4 % 532 Chytridiomycota 2.3 % 12.2 % 2.2 % 0 % 533 Glomeromycota 0 % 3.1 % 3.2 % 2.2 % 534 Unclassified fungi 5.8 % 1.0 % 1.0 % 4.2 % 535 Uncultured fungal clone 54.7 % 45.9 % 9.7 % 16.8 % 536 Total sequenced number 86 98 93 95 537 a)PF, plow-tilled fallow; NF, non-tilled fallow; NH, non-tilled with hairy vetch; NR,

538 non-tilled with cereal rye. Years 2004 is represented as 04. 539 Proofs View

Advance

24

Copyright 2010 by the Japanese Society of Microbial Ecology / the Japanese Society of Soil Microbiology 540 Table 2. Number of T-RFs detected by fungal T-RFLP profiling of upland soils under 541 different treatmentsa) 542 Restriction Number of T-RFs (2004, May / 2007, May) 543 enzymes in soil managementb) 544 PF NF* NH* NR* 545 HaeIII 22 / 33 26 / 6 53 / 14 27 / 9 546 MboI/Sau3AI 51 / 39 37 / 9 50 / 14 22 / 8 547 MspI 51 / 35 41 / 9 53 / 9 29 / 4 a) 548 ITS-D1fLSU regions were amplified under 30 cycles using 0.1-µg environmental DNA. 549 The number of T-RFs of between 50 and 650 bases was counted. 550 b)PF, plow-tilled fallow; NF, non-tilled fallow; NH, non-tilled with hairy vetch; NR, 551 non-tilled with cereal rye. 552 Asterisk indicates a significant difference from T-RF numbers in 2004 (P≤0.05). 553 Proofs View

Advance

25

Copyright 2010 by the Japanese Society of Microbial Ecology / the Japanese Society of Soil Microbiology ITS1 LR21/QLR21Proofs QITS1F ITS4 SSU 5.8S LSU Ribosomal RNA ITS1View ITS2 D1 D2 Primer set T-RFLP profiling QITS1F, ITS4 ITS1, QLR21

Advance Fig. 1. Nishizawa et al.

Copyright 2010 by the Japanese Society of Microbial Ecology / the Japanese Society of Soil Microbiology 60

50 NF

40 PF OTUs NR 30 Proofs NH 20 Number of Number 10

0 View 0 20 40 60 80 100 120 Number of clones analyzed

Advance Fig. 2. Nishizawa et al.

Copyright 2010 by the Japanese Society of Microbial Ecology / the Japanese Society of Soil Microbiology AB 4 10.2 . 0 Proofs 0

02 NF04 NH04 NMDS axis 2 NMDS axis -2 PF07 NH07 View -0.2 -0.1 -4 3 -0. -6 -4 -2 0 2 4 6 -0.2 0 0.2 0.4 NMDS axis 1 NMDS axis 1

Advance Fig. 3. Nishizawa et al.

Copyright 2010 by the Japanese Society of Microbial Ecology / the Japanese Society of Soil Microbiology Cystofilobasidium bisporidii CBS 6347 (AF444299) 100 Cryptococcus nyarrowii CBS 8804 (AY006481) C116 (NR,2 clones; NH, 6 clones) IU-FSC Fun24 65 Cryptococcus sp. CBS 7743 (AF408419) Cryptococcus allantoinivorans CBS 9604 (AY315664) Cryptococcus aureus G7a (DQ640764) 67 Cryptococcus cellulolyticus CBS 8294 (A(AF444442)F444442) 50 Cryptococcus anemochorus CBS 10258 (DQ830986) 63 Cryptococcus rajasthanensis 15L (AM262325) 68 Cryptococcus laurentii NRRL Y-2536 (FJ153214) 91 Cryptococcus nemorosus VKM Y-2906 (AF472628) 91 Cryptococcus perniciosus VKM Y-2905 (AF472627) 96 98 Cryptococcus heveanensis CBS 569 (AF444301) Cryptococcus bestiolae (AY917101) 87 Cryptococcus dejecticola (AY917103) Cryptococcus flavus CBS 331 (AF444338) 89 Cryptococcus podzolicus CBS 6819 (AF444321) Cryptococcus statzellii 87a (AY029342) 94 100 C087 (NH, 1 clone) IU-FSC Fun26Proofs Cryptococcus surugaensis JCM 11903 (AB100440) 67 63 Cryptococcus luteolus ATCC 32044 (EU252551) 100 Cryptococcus zeae HB1207T (AJ965481) Cryptococcus fagi CBS 9964 (DQ054534) 99 Cryptococcus skinneri CBS 5029 (AF444305) Cryptococcus dimennae ATCC 22024 (EU266559) Cryptococcus taibaiensis AS 2.2444 (AY557601) 76 View Cryptococcus foliicola AS2.2471 (AY557600) 98 74 C053 (NH, 1 clone) IU-FSC Fun27 89 Cryptococc u s tephrensis CBS 8935 (D(DQ000318)Q000318) 95 Cryptococcus carnescens CBS 10755 (EU149786) Cryptococcus victoriae HB1221 (AM160647) 62 D071 (NR, 2 clones) C065 (NF, 1 clone; NH, 9 clones; NR, 25 clones) 50 C104 (NH, 1 clone) IU-FSC Fun28 C122 (NH, 2 clones) C058 (NH, 1 clone) Advance0.1

Fig. 4. Nishizawa et al.

Copyright 2010 by the Japanese Society of Microbial Ecology / the Japanese Society of Soil Microbiology