Community profiling and gene expression of fungal assimilatory nitrate reductases in agricultural soil.

Although fungi contribute significantly to the microbial biomass in terrestrial ecosystems, little is known about their contribution to biogeochemical nitrogen cycles. Agricultural soils usually contain comparably high amounts of inorganic nitrogen, mainly in the form of nitrate. Many studies focused on bacterial and archaeal turnover of nitrate by nitrification, denitrification and assimilation, whereas the fungal role remained largely neglected. To enable research on the fungal contribution to the biogeochemical nitrogen cycle tools for monitoring the presence and expression of fungal assimilatory nitrate reductase genes were developed. To the ~100 currently available fungal full-length gene sequences, another 109 partial sequences were added by amplification from individual culture isolates, representing all major orders occurring in agricultural soils. The extended database led to the discovery of new horizontal gene transfer events within the fungal kingdom. The newly developed PCR primers were used to study gene pools and gene expression of fungal nitrate reductases in agricultural soils. The availability of the extended database allowed affiliation of many sequences to known species, genera or families. Energy supply by a carbon source seems to be the major regulator of nitrate reductase gene expression for fungi in agricultural soils, which is in good agreement with the high energy demand of complete reduction of nitrate to ammonium.

Introduction

Agricultural soils normally encounter regular nitrate (NO3−) flushes, either directly from nitrate fertilization or from ammonium (NH4+) or urea, which are rapidly converted to NO3− by nitrification. NO3− can subsequently be lost from agricultural systems via leaching to the groundwater or via denitrification to nitrous oxide (N2O) or dinitrogen (N2) to the atmosphere (Zumft, 1997). Both processes cause environmental problems and, moreover, consumption of NO3−-contaminated drinking water poses health risks to humans and animals and N2O is a potent greenhouse gas (for a recent summary on biogeochemical nitrogen cycles and their impact on the environment, see Galloway ).

Both plants and phylogenetically diverse microorganisms have high NO3− assimilation capacities (Glass ; Stitt ; Inselsbacher ). Incorporation of nitrogen from NO3− into microbial biomass leads to immobilization and therefore reduces losses through leaching or denitrification. The ability to assimilate NO3− is widespread among microbes. In the fungal kingdom, NO3− utilization is restricted to the Dikarya, that is, the Ascomycota and the Basidiomycota (Slot ; Slot and Hibbett, 2007) and certain members from the Mucorales (Sarbhoy, 1965). NO3− assimilation is however not ubiquitous in the Dikarya with many nonutilizing species found especially among the yeasts (Lodder, 1970).

In many fungi, NO3− assimilation genes are clustered and horizontal gene transfer (HGT) from an oomycete to an ancestor of the Dikarya has been proposed (Slot ; Slot and Hibbett, 2007). The presence of a nitrate reductase gene in the genome of Mucor circinelloides could, however, be indicative of an earlier acquisition of the NO3− assimilation cluster by fungi. Within the Dikarya, a secondary HGT event was proposed from a basidiomycetous donor to the ascomycetous genus Trichoderma (Slot and Hibbett, 2007).

Regulation of fungal NO3− assimilation has been extensively studied (Marzluf, 1997; Siverio, 2002) and is today best understood in the model ascomycete Aspergillus nidulans (Hynes, 1973; Cove, 1979; Punt ; Muro-Pastor , 2004; Narendja ; Berger , 2008; Bernreiter ; Schinko ). In general, NO3− uptake and assimilation in fungi is repressed in the presence of preferred nitrogen sources (Wong ), but differences occur at the level of repression. Although NH4+ or glutamine completely repress uptake and utilization of any alternative nitrogen source in the model ascomycetes Neurospora crassa (Premakumar ) and A. nidulans (Arst and Cove, 1973), a more relaxed mode of action was found in several Fusarium species (Celar, 2003) and in the yeast Candida nitratophila (Ali and Hipkin, 1986). Similar variations on the same theme are seen at the pathway-specific level: although the presence of NO3− (or nitrite) is obligatory for full induction of uptake and assimilation in N. crassa and A. nidulans (Hawker ; Marzluf, 1997), this is not the case in the ectomycorrhizal basidiomycetes Hebeloma cylindrosporum (Jargeat , 2003) and Laccaria bicolor (Kemppainen ) and the yeasts C. nitratophila and Sporobolomyces roseus (Ali and Hipkin, 1985).

In good agreement with studies on model organisms under laboratory conditions, assimilatory nitrate reductase activity in aerated soil slurries amended with glucose was found to be strongly inhibited by NH4+ and glutamine (McCarty and Bremner, 1992a, 1992b). On the other hand, high NO3− immobilization rates were found in grassland (Davidson ) and arable (Inselsbacher ) soils even in the presence of NH4+. Whereas the studies by McCarty and Bremner, (1992a, 1992b) employed measurement of NO3− removal in aerated soil slurries amended with glucose and NO3−, the study by Inselsbacher measured incorporation of nitrogen-15 into microbial biomass in a soil microcosm system (Inselsbacher ). Under these experimental conditions, microsites might occur where the repressive effect of NH4+ on NO3− assimilation is alleviated because of contrasting soil mobility of the two inorganic nitrogen forms (Davidson ; Inselsbacher ).

To characterize assimilatory nitrate reductases in environmental samples, tools and data are available for bacterial (Allen ; He ), archaeal (Alcántara-Hernández ) and basidiomycetous (Nygren ) genes. In contrast, the analogous tools and databases have not yet been available for ascomycetes, which constitute the dominant fungal group in many agricultural and grassland soils (Klaubauf ). In this study we developed degenerate primers that allowed the amplification of nitrate reductase encoding genes (euknr) from fungi. The deduced protein sequences are consequently termed EukNR. This naming should avoid confusion with the bacterial nomenclature, where nar is used for the respiratory, nap for the periplasmic and nas for the assimilatory nitrate reductases (Richardson and Watmough, 1999). EukNR is structurally not related to bacterial nitrate reductases but rather to pro- and eukaryotic sulphite oxidases (Campbell, 2001). The main focus in primer design has been put on Ascomycota, and within the Ascomycota on the subphylum Pezizomycotina, which constitute the main group in agricultural soils (Klaubauf ). The newly developed primers have been used on individual culture isolates to improve the reference database of euknr genes. In addition, total DNA isolated from agricultural soils has also been used as template in an attempt to gain a better insight into the fungal community able to utilize NO3−. Finally, gene expression of fungal euknr has been studied in a selected agricultural soil, revealing the most prominent members of the soil fungal community expressing nitrate reductase mRNA.

Materials and methods

Fungal strains and growth conditions

A complete list of fungal strains used in this study is given in Supplementary Table 1. All cultures were maintained at room temperature on malt extract agar (Merck, Darmstadt, Germany) except for Bionectria coronata CBS 696.93 and Bionectria grammicospora CBS 209.93, which were maintained on potato dextrose agar (Roth, Karlsruhe, Germany).

DNA extraction, PCR and sequencing

DNA extraction from axenically growing culture isolates was performed with the DNeasy Plant Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer's instructions. DNA from soils Maissau (soil M) and Purkersdorf (soil P) was taken from Klaubauf . DNA from experiments with soil from Niederschleinz (soil N) was eluted after RNA extraction (see below) with the RNA PowerSoil DNA Elution Accessory Kit (MO BIO, Carlsbad, CA, USA). For field site and soil sampling description, see Inselsbacher .

Primers for amplification of fungal euknr were designed based on multiple alignments of available full-size EukNR sequences. Primer details are given in Table 1 and Supplementary Table 2. Positions of primers in the euknr gene are shown in Figure 1. The routine amplification protocol employed a first PCR step with primer pair niaD01F/niaD04R and a nested PCR with primer pairs niaD15F/niaD12R (Ascomycota) or niaD15F/niaD13R (Basidiomycota), although from pure culture isolates single-round amplification with primer pairs niaD15F/niaD12R or niaD15F/niaD13R usually worked well. From Bionectria ochroleuca and Myxotrichaceae, longer PCR products were generated with more specific primers to improve phylogenetic placement. Primer variants with T7- or M13fwd-sequences attached to the 5′-end were used for direct sequencing of PCR products, as sequencing with highly degenerate primers often resulted in poor-quality reads. The PCR mixtures were as follows: 1 × GoTaq Green Master Mix (Promega, Mannheim, Germany), 2 μ of each primer, 0.5 mg ml–1 bovine serum albumin, 1.25 m additional MgCl2 and template. The cycling conditions were as follows: initial denaturation at 95 °C for 2 min and 30 s followed by 35 cycles of 94 °C for 20 s, 52 °C for 20 s and 72 °C for 1 min followed by a final extension at 72 °C for 5 s.

Table 1

Primers for routine amplification of fungal euknr from pure culture isolates and environmental samples

Primer	Sequence	Expected PCR product size
niaD01F	5′-GTNTGYGCNGGNAA-3′
niaD04R	5′-GTNGGRTGYTCRAA-3′	niaD01F+niaD04R: 1 kb
niaD15F	5′-GGNAAYMGNMGNAARGARCARAA-3′
niaD12R	5′-AACCANGGRTTRTTCATCATNCC-3′	niaD15R+niaD12R: 0.7–1 kb
niaD13R	5′-GGTRGCGTTCCAGTACATRTC-3′	niaD15R+niaD13R: 0.8–1 kb

Figure 1

Position of primers for amplification of euknr genes. Positions of primers for amplification of fungal euknr genes are shown on a schematic map of the EukNR domain structure. Forward primers are indicated above and reverse primers below the map; niaD primers are from this study, and nar primers from Nygren are shown for comparison. For primer details, see Table 1 and Supplementary Table 2. Regions used for phylogenetic studies are indicated at the bottom. Region A was used for construction of the backbone tree and for phylogenetic placement of soil genomic and cDNA clones. Region B was used to obtain improved phylogenetic resolution of horizontal gene transfer events in the Trichoderma/Bionectria-clade and in the Myxotrichaceae (see text and Figure 2). Region C was used by Nygren in their study of EukNR from ectomycorrhizal basidiomycetes.

Sequencing was done by AGOWA GmbH (Berlin, Germany). Accession numbers are HQ115646–HQ115737 for rRNA gene cluster sequences for identification of pure culture isolates, HQ234766–HQ234874 for partial euknr from pure culture isolates (see Supplementary Table 1), HQ243445–HQ243656 for partial euknr amplified from soil genomic DNA and HQ288852–HQ288891 for partial euknr amplified from soil complementary DNA (cDNA).

Microcosm experiment

The experimental system for analysing soil fungal RNA expression consisted of microcosms filled with agricultural soil N as described in Inselsbacher , 2010). Soil was amended with 1% (w/w on a fresh weight basis) of a cellulose/starch/xylan mixture before setup of the experiment. Treatment A received no additional organic fertilizer, treatment B received 0.05% (w/w on a fresh weight basis) Biovin (Intertrest, Guntramsdorf, Austria) and treatment C received 0.05% autoclaved Biovin. Biovin is an organic fertilizer produced from humified grapevine pomace. Each treatment was replicated ten times (five replicates for each harvest time point). Microcosms were equilibrated to 62% water-filled pore space, fertilized twice with potassium nitrate (KNO3) and planted with winter barley as described in Inselsbacher , 2010). Sampling time points for soil RNA isolation and determination of soil nitrogen pools corresponded to days 1 and 6 in Inselsbacher , that is, 1 and 6 days after planting the pregerminated seedlings and the second fertilization. Only microcosms with fully developed winter barley seedlings were used for further analyses. At the first sampling day, four A, five B and four C treatments and at the sixth day, four A, four B and five C treatments had been fully developed and were used for further analyses. Soil chemical analyses were done as described in Inselsbacher .

RNA isolation and cDNA synthesis

At harvest the whole content of one microcosm was transferred to a petri-dish, plant shoots were removed and plant roots were carefully separated from the soil. Soil was mixed and an aliquot (2 g fresh weight) was used for RNA isolation with the PowerSoil Total RNA Isolation Kit (MO BIO) according to the manufacturer's instructions with a minor modification: the organic phase and the interphase were re-extracted by addition of 2.5 ml Bead Solution, 0.25 ml SR1 and 0.8 ml SR2, followed by thorough vortexing and centrifugation. The aqueous layers from the first and second extraction were combined. Amounts of buffers were doubled for the subsequent steps until loading onto the RNA capture columns. Integrity of RNA was checked on a conventional 1.5% agarose gel, and the concentration determined spectrophotometrically with the Quant-iT RNA BR Kit (Invitrogen, Paisley, UK). From each replicate, 8 μl of total RNA (∼0.1 μg) was treated with 1 U DNase I (Invitrogen) in a volume of 10 μl for 15 min at room temperature. Reactions were terminated by addition of 1 μl 25 m EDTA and heating to 65 °C for 10 min. DNase-I-treated RNA was directly subjected to reverse transcription with random hexamer primers with RevertAid H− First Strand cDNA Synthesis Kit (Fermentas, St Leon-Rot, Germany) in a volume of 20 μl according to the manufacturer's instructions. Single-stranded cDNA was diluted to 100 μl and directly used for PCR with euknr primers as described above. Control reactions included PCR with untreated RNA and RNA after DNase I treatment as templates and PCR with primers EF1α983F/EF1α1567R for amplification of fungal translation elongation factor 1α (for primer details, see http://aftol1.biology.duke.edu /pub/primers/viewPrimers). No PCR products could be obtained from RNA after DNase I treatment but before reverse transcription.

At sampling time point day 1, nested PCR on cDNA with primers designed for ascomycetous euknr gave products of the expected size in two, four and three replicates from treatments A, B and C, respectively. PCR products from day 1 were pooled per treatment (A, B and C) before cloning. As from day 6, euknr could only be detected with a weak signal in 1 out of 13 RNA samples, this time point was omitted from further analyses.

Cloning

PCR products were cloned in plasmid pTZ57R/T (Fermentas) according to the manufacturer's instructions. Cloning of euknr-PCR products from pure culture isolates was only necessary in one case (strain Phialophora mustea KSS12_F02), which showed overlapping signals upon direct sequencing of PCR products. In all, 12 randomly selected clones were picked and sequenced from both ends with vector primers. For Cadophora finlandica, a genomic clone containing a nearly full-length euknr including the promoter region was isolated from an existing library (Gorfer ) by PCR amplification with euknr and plasmid primers (Acc. Nr. HQ234780).

For analysis of fungal nitrate reductase genes in agricultural soils, euknr was amplified from genomic DNA from soils M and P with primer pairs designed for ascomycetous and basidiomycetous euknr, respectively. PCRs were run in triplicate and pooled before cloning. From each of the resulting four libraries, 96 clones were randomly picked. For soil N from the microcosm experiment, only the primer pair designed for the ascomycetous euknr was used. From each of the resulting three libraries, 72 clones were randomly picked. For analysis of euknr cDNA pools in soil N, only the primer pair designed for the ascomycetous euknr was used. Amplification was done separately on the replicates from harvest day 1 and pooled per treatment before cloning. From the pools A, B and C, 85, 46 and 39 clones were picked, respectively.

Clones were first screened by restriction fragment length polymorphism with the restriction enzyme BsuRI (Fermentas isoschizomere of HaeIII) and up to five randomly selected clones from each restriction fragment length polymorphism pattern were sequenced by vector primers from both ends.

Sequence analysis

Forward and reverse sequence reads were assembled using Vector NTI Advance 10 and 11 (Invitrogen). Introns were detected based on conserved intron positions, conserved boundary sequences and deduced amino-acid sequence after in silico splicing of the introns. Edited nucleotide sequences were translated to amino-acid sequences by Vector NTI. Multiple alignments of protein sequences (MUSCLE; Edgar, 2004) and construction of phylogenetic trees were done with molecular evolutionary genetics analysis version 4 (MEGA4; Tamura ) and MEGA5 (Tamura ). Maximum likelihood trees were calculated with the WAG+G model (Whelan and Goldman, 2001) and 500 bootstrap replicates.

Chimaeric reads from genomic and cDNA libraries were detected by Bellerophon (Huber ) and removed before further analyses. For affiliation of soil genomic and cDNA clones, sequences were clustered with MOTHUR (Schloss ) at 96% identity (furthest neighbour) at the protein level. At this cutoff, separation of most species was possible, although some species (for example, in the genus Tetracladium) still clustered together whereas, for example, Fusarium oxysporum is divided into two subgroups. With these settings, clones with identical BsuRI restriction patterns consistently clustered together.

Phylogenetic hypothesis testing

Four replicates each of constrained and unconstrained trees were calculated with RaxML 7.0.3 (Stamatakis, 2006) using the PROTMIXWAGF setting, an implementation of the WAG+G model of amino-acid sequence evolution (Whelan and Goldman, 2001), as selected with MultiPhyl v1.0.6 (Keane ). Single-site likelihoods were calculated in RaxML using the PROTGAMMAWAGF setting and imported into CONSEL 0.1i (Shimodaira and Hasegawa, 2001) to compare the tree topologies by a series of statistical tests, including the approximately unbiased test of phylogenetic tree selection (Shimodaira, 2002).

Statistical analyses

Differences in soil nitrogen pools between soils were tested by one-way analysis of variance. Cluster analysis (nearest neighbour method with squared Euclidian distance metric) of euknr communities in soils M and P amplified with different primer pairs was performed to investigate the similarity between and within soil samples. All statistical analyses were run with StatGraphics version 5 (StatPoint Technologies, Inc., Warrenton, VA, USA).

Results

The phylogenetic tree of fungal EukNR sequences

The fungal EukNR tree is based on published full-length or nearly full-length sequences, mainly from fungal genome sequencing projects, but also from early studies on structure and function of EukNR. An overview of these sequences is given in Supplementary Table 3. The full-length data set consists of three Phytophthora sequences, which represent the putative oomycete donor of horizontally transferred fungal EukNR (Slot and Hibbett, 2007), 1 sequence from Mucoromycotina, 9 sequences from Basidiomycota and 72 sequences from Ascomycota. Additionally, one sequence from an Oryza sativa EST library (AK110249; Kikuchi ) was included. This sequence obviously encodes a EukNR from a fungus related to Ustilago maydis, which had apparently infected the analysed rice cultivar. No corresponding sequence was found in the rice genome, excluding the possibility that the sequence was obtained by HGT. Phylogenetic analysis of the full-length data set provides evidence that euknr was transferred from an oomycete donor to a common ancestor of Mucorales and Dikarya (95% bootstrap support for placing M. circinelloides EukNR sister to dikaryan EukNR; see Supplementary Figure 1). The rare but evenly spread distribution of nitrate assimilation in the Mucorales (cf Domsch ) suggests frequent losses of euknr from mucoralean genomes, as for example in Rhizopus oryzae (Ma ). More sequence data from the basal fungal lineages are needed for a reliable reconstruction of the gene gains and losses that finally led to today's distribution of euknr in the fungal kingdom.

As many fungal species commonly found in agricultural soils have no close relatives with completely sequenced genomes, additional euknr sequences were generated from pure culture isolates. To this end, 115 ascomycetes and 3 basidiomycetes were analysed for the presence of euknr by PCR. From 104 ascomycetes and 3 basidiomycetes, euknr could be amplified. All PCR products showed highest similarity to fungal nitrate reductases, and not to any other proteins also containing the Molybdenum cofactor like sulphite oxidases. In only two cases, gene duplication of euknr was found: for B. ochroleuca Rd0801 (see below) and Phialophora mustea KSS12_F02. EukNRs from P. mustea KSS12_F02 have 78.8% sequence identity and cluster together in a sub-branch of the Hypocreales (see Supplementary Figure 2). In both copies, essential residues for binding of Molybdenum cofactor are conserved. This pattern is indicative for a gene duplication in a recent ancestor of P. mustea KSS12_F02.

Overall, there is good congruence of the fungal EukNR tree with widely accepted phylogenies (Hibbett ). In the EukNR tree based on the full-length sequences, all orders appear monophyletic, although the Hypocreales (excluding Trichoderma, see Slot and Hibbett, 2007) and Onygenales lack bootstrap support (see Supplementary Figure 1). The EukNR tree based on the partial Molybdenum cofactor-binding domain (region A in Figure 1) provides increased phylogenetic coverage because of inclusion of additional taxa, but shows decreased resolution at certain nodes (Supplementary Figure 2).

In two cases, indications of HGT within the Dikarya could be found. Although acquisition of a nitrate assimilation gene cluster in the genus Trichoderma from a basidiomycete related to Ustilago via HGT was already previously described (Slot and Hibbett, 2007), the origin and destination of this cluster could now be narrowed down. According to our phylogenetic data, the donor belonged to the Ustilaginaceae and the acceptor was a common ancestor of Trichoderma/Hypocrea and Bionectria (Figure 2). As one out of seven analysed B. ochroleuca strains contained two euknr genes—one ascomycetous and one basidiomycetous—it seems likely that HGT occurred before loss of the autochthonous euknr. Presence of an ascomycetous euknr in the remaining B. ochroleuca strains cannot be excluded, but in none of the completely sequenced Trichoderma genomes an ascomycetous euknr could be found.

Figure 2

Horizontal gene transfer of euknr. The EukNR phylogenetic tree provides evidence for horizontal gene transfer of euknr from a basidiomycete related to the Ustilaginaceae to a common ancestor of the ascomycetes Trichoderma and Bionectria (Hypocreales; Slot and Hibbett, 2007 and this study) and from an ascomycete related to the Eurotiales and Onygenales to an ancestor of the Myxotrichaceae, which are closely related to the Helotiales (this study). Tree branches affected by gene transfer events are highlighted in bold. The arrows indicate different positioning of the involved branches between accepted phylogenies and the EukNR tree (shown here). The phylogenetic tree is based on the partial sequence of the Molybdenum cofactor (MoCo)-binding domain and the partial cytochrome b5 domain (region B in Figure 1) of available full-length or nearly full-length EukNR sequences and partial sequences from Trichoderma, Bionectria and the Myxotrichaceae. Where appropriate, branches were collapsed at the ordinal level (grey triangles) to reduce complexity of the image. Bootstrap support for branches above 50% is shown.

Another HGT event seems to have occurred within the Ascomycota. EukNR from the Myxotrichaceae groups with the Eurotiales and Onygenales. Formerly the genera Geomyces, Pseudogymnoascus and Oidiodendron were placed in the Onygenales (Currah, 1985), although phylogenetic studies recently moved this family close to the Helotiales within the Leotiomycetes (Tsuneda and Currah, 2004; Wang , 2006b). There is, however, excellent bootstrap support for placement of EukNR from G. pannorum L01, Geomyces sp. NG_p41, P. roseus NG_p50 and O. cerealis NG_p39 sister to the Eurotiales and Onygenales, and additional support comes from intron numbers and positions in the amplified euknr fragment: although all analysed euknr genes in the subclass Eurotiomycetidae and in the family Myxotrichaceae show three introns at conserved positions within the amplified fragment, no intron is found in the corresponding region in the Helotiales (see Supplementary Figure 3). Constrained trees, where EukNR from Myxotrichaceae was forced to group with the Helotiales, were rejected by all topological tests applied, including the conservative weighted Shimodaira–Hasegawa test (Shimodaira and Hasegawa, 2001).

Fungal euknr genes in agricultural soils

Agricultural soil fungal communities with the potential to assimilate NO3− were studied by direct amplification of fungal euknr genes from total soil DNA. The following soils from Lower Austria were analysed: Maissau (soil M), Purkersdorf (soil P) and Niederschleinz (soil N). For details on treatment of soil N, see Materials and methods. Soils M and P were analysed in parallel with primers designed for ascomycetous and basidiomycetous euknr, respectively. Surprisingly, both primer pairs amplified mainly ascomycetous euknr genes from soils M and P. Only one basidiomycetous euknr gene affiliated to Coprinellus was detected in soil M. Libraries amplified from the same soil had a significantly lower Euclidian distance (that is, were significantly more similar) than libraries amplified from the two different soils. Therefore, data were merged to two soil-specific data sets. Based on results from soils M and P, only the primer pair designed for the ascomycetous euknr was applied to soil N. Although some sequences from the genomic libraries could be clearly assigned to a distinct species or genus, the majority of sequences could only be assigned to higher taxonomic ranks (see Figure 3 and Table 2). From sequences that showed clear affiliation at the ordinal level, the majority belonged to the Sordariales, Hypocreales and Helotiales. In soil N, however, a large number of clones (39%) belonged to the Peronosporales, and one sequence clustered with M. circinelloides. These results indicate that the primers for amplification of ascomycetous euknr also show a certain degree of unspecificity.

Figure 3

Phylogenetic placement of soil euknr clones. The phylogenetic tree is based on the partial Molybdenum cofactor (MoCo)-binding domain of EukNR (region A in Figure 1). Phytophthora spp. were used as outgroup, branches lacking environmental sequences were collapsed at common phylogenetic levels. The upper part of the tree is collapsed in (a) and shown fully expanded in (b). Environmental sequences are labelled as follows: genomic clone soil M, genomic clone soil P, genomic clone soil N, cDNA clone soil N. Bootstrap support for branches above 50% is shown. Filled green circle, genomic clone soil M; filled yellow square, genomic clone soil P; filled blue triangle, genomic clone soil N; open blue triangle, cDNA clone soil N.

Table 2

Tentative phylogenetic affiliation and relative abundance (%) of euknr clones amplified from agricultural soil DNA

Ordera	Speciesb	Agricultural soilc
		M	P	Ntot	Nfun
Pezizomycotina i.s.	Pezizomycotina sp. 1a	2.1
	Pezizomycotina sp. 1b		3.0
	Pezizomycotina sp. 2		20.1
	Pezizomycotina sp. 3	1.4
	Pezizomycotina sp. 4		0.6
	Pezizomycotina sp. 5		5.3
	Pezizomycotina sp. 6		3.6
	Pezizomycotina sp. 7		0.6
	Pezizomycotina sp. 8		5.3
Capnodiales	Cladosporium herbarum complex	0.7
	Devriesia sp. NG_p52	2.1
Pleosporales	Pleosporales NG_p49			1.7	2.9
	Pleosporales sp. 1	0.7
Chaetothyriales	Exophiala salmonis		1.2
Leotiomycetes i.s.	Geomyces sp.			0.6	1.0
Helotiales	Tetracladium marchalianum complex			1.2	1.9
	Tetracladium sp. 1		13.0
	Tetracladium sp. 2	3.5
Sordariomycetes i.s.	Sordariomycetes sp. 1			1.7	2.9
	Sordariomycetes sp. 2	19.7
	Sordariomycetes sp. 3	5.6
	Sordariomycetes sp. 4	29.6
Hypocreales	Fusarium sp. 1		2.4
	Fusarium sp. 2			1.2	1.9
	Fusarium sp. 3			0.6	1.0
	Hypocreales NG_p26			1.2	1.9
	Hypocreales sp. 1a			0.6	1.0
	Hypocreales sp. 1b			5.2	8.6
Microascales	Doratomyces sp.		21.3
Phyllachorales	Gibellulopsis nigrescens	7.0		5.8	9.5
Sordariales	Lasiosphaeriaceae sp. 1	7.7	4.1
	Lasiosphaeriaceae sp. 2	1.4	4.7
	Lasiosphaeriaceae sp. 3		3.0
	Lasiosphaeriaceae sp. 4			1.7	2.9
	Lasiosphaeriaceae sp. 5	16.9		17.3	28.6
	Lasiosphaeriaceae sp. 6		2.4
	Lasiosphaeriaceae WF151		4.7
	Sordariales sp. 1			4.6	7.6
	Trichocladium asperum related 1			5.2	8.6
	Trichocladium asperum related 2			11.6	19.0
	Trichocladium sp. 1	0.7
Xylariales	Arthrinium sp. 1a		0.6
	Arthrinium sp. 1b		4.1
Agaricales	Coprinellus sp.	0.7
Mucoromycotina i.s.	Mucoromycotina sp. 1			0.6	1.0
Peronosporales	Peronosporales sp. 1			38.7
Stramenopiles i.s.	Stramenopiles sp. 1			0.6
Number (n)d		142	169	173	105
		100.0	100.0	100.0	100.0

Species are sorted according to class and order, with Ascomycota on the top followed by Basidiomycota (Coprinellus sp. 1), Mucoromycotina and the non-fungal phylum Stramenopiles (Peronosporales sp. 1 and Stramenopiles sp.).

Species epithets like sp. 1, 2 and so on indicate different species of unknown affiliation, where 1a and 1b denote closely related species, whereas epithets like NG_p52 indicates identity at the 96 % threshold level (protein sequences) to insufficiently identified species from the reference database (see Supplementary Table 1).

Agricultural soils from Maissau (M), Purkersdorf (P) and Niederschleinz (N). Data for soil N treatments A, B and C (see Materials and methods) were merged and once calculated for all clones including non-fungal sequences (Ntot) and once for fungal sequences only (Nfung).

Absolute number of clones analysed per library.

As not unusual for fungal communities, the euknr gene pool is very uneven, with the three most dominant species accounting in all three agricultural soils for >50% of the clones.

Fungal euknr gene expression in agricultural soil

Expression of ascomycetous euknr genes was examined in agricultural soil N. As no fungal euknr gene expression could be detected in soil N even after fertilization with KNO3 (from experiments described in Inselsbacher ), a second set of experiments was conducted where additional carbohydrates were supplied to the soil microbial community. Soil N was amended with a mixture of cellulose, starch and xylan and fertilized with KNO3. Treatments A, B and C differed by addition of organic fertilizer. Organic fertilizer at the applied dose had no effects on soil pools of NO3–, NH4+ and microbial nitrogen at both sampling time points (data not shown). There was, however, a major effect of carbohydrate addition on soil nitrogen pools: NH4+ content and microbial nitrogen increased slightly, whereas NO3– levels dropped tremendously (Table 3). Accordingly, no fungal euknr gene expression could be detected in soil N without carbohydrate addition, whereas gene expression could consistently be detected in carbohydrate-amended soil N at sampling day 1. At sampling day 6, fungal euknr gene expression could be detected in only 1 out of 13 replicates. Analysis of the fungal community actively expressing euknr at day 1 revealed five species from the ascomycetes and one from the Stramenopiles. The latter was only found in treatment A, but at a very high number. Although four out of the five fungal euknr cDNAs were only detected in one of the three treatments, cDNAs affiliated to Trichocladium asperum related 2 were found in all three treatments. Clone abundances in cDNA libraries were always higher than in genomic libraries. A comparative overview of the results is shown in Table 4 and Figure 3.

Table 3

Nitrogen pools and fungal euknr gene expression in microcosm soils

	Soil Na	Soil N + CHb	Factorc	P-valued
Day 1e
NH4+ nitrogenf	1.08±0.31	2.22±0.55	× 2	<0.0001
NO3− nitrogeng	67.7±12.5	0.35±0.062	÷200	<0.0001
Microbial nitrogenh	45.7±20.5	198±28.7	× 4	<0.0001
ni	15	13
tefj	15	12
euknrk	0	9
Day 6
NH4+ nitrogen	0.93±0.31	1.85±0.40	× 2	<0.0001
NO3− nitrogen	5.58±3.83	0.19±0.063	÷30	<0.0001
Microbial nitrogen	45.6±17.4	172±19.5	× 4	<0.0001
n	10	14
tef	10	14
euknr	0	1

Agricultural soil from Niederschleinz (N) fertilized with KNO3 equivalent to 36 μg nitrogen per gram dry weight; for experimental details, see Materials and methods. Data for nitrogen pools in soil N without carbohydrate addition were taken from Inselsbacher et al. (2010).

As above, but with addition of carbohydrates (mixture of cellulose, starch and xylan).

Approximate x-fold change in nitrogen pools between experiments without and with carbohydrate addition.

Statistical significance of the differences between treatments (one-way analysis of variance (ANOVA)).

Sampling day.

Plant-available ammonium nitrogen in the soil in μg g−1 dry weight (means±s.d.).

Plant-available nitrate nitrogen in the soil in μg g−1 dry weight (means±s.d.).

Microbial nitrogen in the soil in μg g−1 dry weight (means±s.d.).

Number of microcosms analysed.

Number of microcosms where fungal tef (translation elongation factor 1α) gene expression could be detected.

Number of microcosms where fungal euknr gene expression could be detected.

Table 4

Tentative phylogenetic affiliation of fungal euknr complementary (c)DNA and relative expression levels in microcosms with agricultural soil N (soil from Niederschleinz) after addition of potassium nitrate (KNO3) and carbohydrates

		Aa			B			C
	gb	cc	c/gd	g	c	c/g	g	c	c/g
Fusarium sp. 2				6 (2)				69 (25)
Lasiosphaeriaceae sp. 4	3 (1)			3 (1)			3 (1)	3 (1)	1.0
Lasiosphaeriaceae sp. 5	21 (7)	51 (24)	2.4	40 (14)			25 (9)
Trichocladium asperum rel. 1				14 (5)	13 (6)	0.9	11 (4)
Trichocladium asperum rel. 2	21 (7)	49 (23)	2.2	14 (5)	87 (39)	6.1	22 (8)	28 (10)	1.3
Stramenopiles sp.2		- (30)
Ascomycotae	100 (33)	100 (47)		100 (35)	100 (45)		100 (36)	100 (36)

Treatments differed in application of organic fertilizer (A: no organic fertilizer; B: 0.05% Biovin; C: 0.05% autoclaved Biovin). For further details on the experimental setup, see Materials and methods.

Relative clone abundances (rounded to integer) and absolute clone numbers (in parentheses) in the genomic libraries. Only ascomycetous euknr sequences were considered for calculation of the relative abundances and only species for which euknr cDNA could be found at least in one of the three treatments are shown.

Relative clone abundances (rounded to integer) and absolute clone numbers (in parentheses) in the cDNA libraries. Only ascomycetous euknr sequences were considered for calculation of the relative abundances.

Ratios of relative clone abundances in genomic and cDNA libraries (c/g) were calculated as approximations for expression levels. No relative expression levels were calculated for euknr from Fusarium sp. 2, which was found at the genomic level in treatment B and at the cDNA level in treatment C. Stramenopiles sp. 2 was only found in the cDNA library from treatment A, but in no genomic library.

Number of ascomycetous euknr sequences in the respective libraries, on which the calculation of the relative abundances was based.

Discussion

The current data set for the fungal EukNR phylogenetic tree contains 191 sequences, 82 previously published and 109 newly obtained, and covers all major ascomycete orders found in arable soils (Klaubauf ). There are, however, still substantial gaps. No sequences can currently be assigned to soil clone group I, a group of uncultured fungi at the base of the Ascomycota (Porter ) that can occur at high frequency, for example, in grassland soil (Klaubauf ). As these fungi have not been yet cultivated, it is currently unknown whether or not they can utilize NO3−. From the Taphrinomycotina there are a few completely sequenced genomes, but none of them contains a putative euknr. Within the subphylum Pezizomycotina the class Lecanoromycetes of lichen-forming fungi is not yet represented in the EukNR tree. Based on the available data it seems likely that several events of gains and losses led to today's distribution of euknr in the fungal kingdom: primarily obtained by HGT from an oomycetous donor to an ancestor of the Dikarya (Slot and Hibbett, 2007) or even earlier to a common ancestor of Mucoromycotina and Dikarya (this study, see Supplementary Figure 1), euknr was subsequently frequently lost, mainly in fungi with a yeast-like form in their life cycle. Many species from not only the Saccharomycotina and the Taphrinomycotina, but also basidiomycetous yeasts like Cryptococcus or the dimorphic ascomycete Lecythophora sp. NG_p46 (see Supplementary Table 1) are unable to assimilate nitrate. Good evidence for secondary HGT events is available in two cases: (1) the genera Trichoderma/Hypocrea- and Bionectria-obtained euknr most likely from an ancestor closely related to the Ustilaginaceae (Slot and Hibbett, 2007this study), presumably before they lost their corresponding autochthonous gene, as evidenced by the presence of two copies of euknr in one strain from B. ochroleuca; and (2) the Myxotrichaceae obtained euknr from an ancestor within the Eurotiomycetidae. For other genera with unusual placement of EukNR in the phylogenetic tree (like for example, Doratomyces sp.), bootstrap support is currently lacking to make statements about HGT events and more data are needed. With these exceptions, clustering of most taxa according to established phylogenies (Hibbett ) is found in the available EukNR data set. Good correspondence was found between number and position of introns and phylogenetic affiliation, with six introns at conserved positions for all thus far analysed members in the Eurotiales and Onygenales (class Eurotiomycetes) and one intron for all analysed members in the Hypocreales and the Sordariales (class Sordariomycetes). The newly detected HGT event from a Eurotiomycete to an ancestor of the Myxotrichaceae is also supported by conserved intron positions (see Supplementary Figure 3).

As evidenced by the high number of sequences from soil clone libraries that could only be assigned to higher taxonomic levels, expansion of the database of fungal EukNR sequences is fundamental for an accurate affiliation of sequences from environmental samples. The need for a database with good coverage of EukNR from all major fungal lineages is further documented by the common occurrence of HGT events in the evolution of fungal nitrate assimilation. The availability of EukNR sequences beyond the borders of the Molybdenum cofactor-binding domain would moreover be important in data analysis of metagenomic and metatranscriptomic studies.

The relatively high number of agricultural EukNR sequences belonging to the orders Sordariales, Hypocreales and Helotiales is in good agreement with the structure of the fungal communities in these soils as analysed by ITS/LSU (internal transcribed spacer and partial large subunit rRNA gene) libraries (Klaubauf ). The dominance of these orders was even more pronounced at the euknr gene expression level in soil N, where all ascomycetous euknr sequences belonged to the Sordariales or Hypocreales. Although many sequences at both the community (=genomic DNA) and gene expression (=cDNA) level could only be loosely affiliated to known families, orders and so on, some sequences could clearly be affiliated to genera or even species. Noteworthy, a set of euknr sequences showing consistently high expression levels in soil N amended with NO3− and carbohydrates was closely related to T. asperum. Fungi related to T. asperum were detected in two out of five investigated agricultural soils, in one of which it was the most abundant species (Klaubauf ). although fungi related to T. asperum were not detected in native soil N, they were probably present below detection limit and thrived under the applied conditions, that is, carbohydrate addition and nitrate fertilization. The availability of cultured isolates opens up the possibility to study in more detail the regulation of euknr gene expression of autochthonous fungi in agroecosystems. Regulation of euknr gene expression was thus far mainly studied in response to nitrogen sources (Cove, 1979; Marzluf, 1997; Berger , 2008; Kemppainen ; Schinko ) and rarely in response to carbon sources (Hynes, 1973; Choudary and Ramananda Rao, 1982). The current data suggest that probably carbon source availability acts as the primary signal. In soil N without addition of carbohydrates, NO3− levels remain high, and expression levels of fungal euknr are undetectable. Fungi are, however, not metabolically inactive under these conditions as evidenced by expression of tef (Table 3). After addition of carbohydrates, fungal euknr gene expression is turned on and NO3− levels drop well below levels of available NH4+. Concomitantly, microbial nitrogen increases, suggesting that at least part of the NO3− is assimilated by fungi for incorporation into biomass. After prolonged incubation (day 6), euknr gene expression levels drop, probably because of the absence of induction by NO3− and increased competition for NO3− by the developing plants (Inselsbacher ). In a comparative study, organic farming employing high organic matter inputs had a significantly positive effect on nitrate assimilation by soil microbes compared with conventional farming (Burger and Jackson, 2003). The dominant role of carbon source on fungal euknr gene expression is in good agreement with previous reports that show carbon source dependence of nitrate reductase activity in Candida utilis (Choudary and Ramananda Rao, 1982) and a recent meta-analysis identifying the ratio of organic carbon to NO3− as the main regulator for the fate of NO3− in terrestrial and aquatic ecosystems (Taylor and Townsend, 2010). NO3− reduction is an energy-demanding process requiring strong supply of reduced nicotinamide adenine dinucleotide phosphate from the pentose phosphate pathway in fungi. Consequently, several enzymes of the pentose phosphate pathway are co-regulated with the nitrate assimilation genes in A. nidulans (Hankinson, 1974; Hankinson and Cove, 1974; Schinko ).

The presented data will certainly contribute to our understanding of the role of fungi in biogeochemical nitrogen cycles. The availability of partial euknr sequences from a broad taxonomic range of Ascomycota allows incorporation of euknr into functional analytical tools like the GeoChip (He ). Furthermore, it will improve phylogenetic affiliation of euknr sequence tags from metatranscriptomic approaches. It must be highlighted that the current data set contains euknr sequences from several orders without completely sequenced genomes like the Botryosphaeriales, Chaetosphaeriales, Chaetothyriales, Diaporthales, Dothideales, Microascales and Xylariales. Future work will further expand the ascomycetous data set and explore groups herein largely neglected like the Basidiomycota or Mucoromycotina. Moreover, the availability of a comprehensive euknr data set will allow the development of efficient tools for quantification of gene expression by real-time PCR techniques in a broad range of environmental samples.