Table 1 Primer Sequences, Peptide Templates and Positions

Table 1 Primer Sequences, Peptide Templates and Positions

Supplementary materials

Table S1 Information of the selected protein sequences of nitrate and nirtrous oxide reductases (coded by narG and nosZ) from GenBank, including the accession numbers and the names of strains

narG / nosZ
Accession no. a / SOURCE / Accession no. / SOURCE
1 / YP_001171435 / Pseudomonas stutzeri A1501 / NP_435889 / Sinorhizobium meliloti 1021
2 / BAD82906 / Pseudomonas sp. MT-1 / CAB53351 / Paracoccus denitrificans PD1222
3 / AAB50620 / Pseudomonas fluorescens / ABQ37945 / Bradyrhizobium sp. BTAi1
4 / NP_252564 / Pseudomonas aeruginosa PAO1 / CAA75425 / Achromobacter cycloclastes
5 / YP_001257334 / Brucella ovis ATCC 25840 / CAE27502 / Rhodopseudomonas palustris GA009
6 / YP_223631 / Brucella abortus biovar 1 str. 9-941 / YP_108223 / Methylobacterium sp. 4-46
7 / YP_983978 / Polaromonas naphthalenivorans CJ2 / YP_001059095 / Paracoccus denitrificans PD1222
8 / YP_741846 / Alkalilimnicola ehrlichei MLHE-1 / ZP_01848502 / Rhodopseudomonas palustris BisA53
9 / AAT47523 / Halomonas maura / CAA64426 / Pseudomonas sp. MT-1
10 / YP_310854 / Shigella sonnei Ss046 / YP_789757 / Pseudomonas aeruginosa PAO1
11 / NP_707133 / Shigella flexneri 2a str. 301 / NP_252082 / Pseudomonas aeruginosa PA7
12 / YP_402918 / Shigella dysenteriae Sd197 / ABR85628 / Burkholderia pseudomallei K96243
13 / NP_455739 / Salmonella enterica subsp. Enterica
serovar Typhi str. CT18 / YP_333650 / Burkholderia pseudomallei S13
14 / YP_960347 / Marinobacter aquaeolei VT8 / YP_102711 / Burkholderia mallei ATCC 23344
15 / ZP_01331210 / Burkholderia thailandensisE264

a The accession numbers were locus of the protein sequences of nitrate and nirtrous oxide reductases encoded by narG and nosZ from these bacteria.

Table S2 Bacterial strains used in this study to test the specificity of the narG andnosZprimers

Bacterial strains / narG
145F/773R / narG
571F/773R / nosZ
1126F/1884R / nosZ
1126qF/1381R / Source
Denitrifying strains with narG and nosZ
Pseudomonas stutzeri A1501 / + / + / + / + / CAAS
Pseudomonas aeruginosa ATCC9027 / + / + / + / + / GIM
Paracoccus denitrificans ATCC17741 / +- / +- / + / + / ACCC
Achromobacter xylosoxidans ATCC15173 / - / - / - / + / CGMCC
Ralstonia eutropha H16 / + / + / + / + / DSM
Ralstonia eutropha CH34 / + / + / + / + / DSM
Bordetella petrii DSM12804 / + / + / + / + / DSM
Thiobacillus denitrificans DSM12475 / + / + / + / - / DSM
Burkholderia thailandensis E264 / + / + / - / + / DSM
Chromohalobacter salexigensDSM3043 / + / + / + / + / DSM
Herbaspirillum frisingense GSF30 / + / + / + / + / DSM
Denitrifying strains with nosZ
Herbaspirillum seropedicae Z67 / - / - / DSM
Azospirillum lipoferum Sp59b / + / + / DSM
Azospirillum brasilense Sp7 / - / - / - / - / DSM
Azospirillum halopraeferens Au4 / - / - / - / + / DSM
Azospirillum doebereinerae GSF71 / - / - / + / + / DSM
Sinorhizobium meliloti1021 / - / - / + / +
Rhizobium meliloti sp. / - / - / + / + / IAE
Denitrifying strains with narG
Shigella flexneri GIM1.231 / + / + / GIM
Shigella dysenteriae GIM1.236 / + / + / GIM
Escherichia coli B7A / + / + / GIM
Ochrobactrum anthropi FZX-1 / + / + / GIM
Burkholderia soli DSM18235 / + / + / GIM
Firmicutes, Gram+
Staphylococcus epidermidissp. / + / + / - / - / IAE
Bacillus licheniformis ATCC14580 / + / + / - / - / DSM
Bacillus cereus sp. / + / + / - / - / IAE
Bacillus subtilis sp. / + / + / - / - / ISA
Corynebacterium glutamicum ATCC13032 / - / - / - / - / DSM
Pseudonocardia sp. / - / - / - / - / IAE
Actinomadura sp. / - / - / - / - / IAE

+, PCR product of the right size obtained; -, no PCR product obtained; + PCR product was cloned and sequenced and the sequence was deposited in the data banks. Pseudomonas stutzeri A1501 were provided by prof. Min Lin, The Chinese Academy of Agricultural Sciences (CAAS).

Table S3 Information of standard curves of the denitrifying genes

Gene / Standard curve / Linear range / R2 / Efficiency (%)
16S rRNA / Y=-3.60X+37.46 / 103~108 / 0.999 / 91
narG / Y=-3.83X+39.70 / 103~109 / 0.996 / 82
nirK / Y=-3.34X+38.33 / 104~108 / 0.995 / 99
nirS / Y=-2.82X+35.10 / 103~108 / 0.995 / 120
qnorB / Y=-3.13X+33.44 / 104~109 / 0.993 / 108
nosZ / Y=-3.39X+36.05 / 102~107 / 0.997 / 97

Table S4 The clonesand sequences information of the denitrifying genes

Fertilization
treatments / Sequenced clones / Clones similarity to the targeted gene / OTUsa
narG / qnorB / nosZ / narG / qnorB / nosZ / narG / qnorB / nosZ
NF / 80 / 80 / 85 / 73 / 76 / 78
UR / 89 / 80 / 80 / 86 / 76 / 64
BM / 80 / 78 / 84 / 71 / 78 / 75
BMR / 95 / 80 / 82 / 92 / 77 / 72
Total / 344 / 318 / 331 / 322 / 307 / 289 / 143 / 205 / 103

a Only the total OTUs from the four treatments of each gene were shown here as the OTUs from the four treatments overlapped.

Fig S1 Relative abundances of narG(a), qnorB (b) and nosZ(c)T-RFs from triplicate plots of long-term fertilization treatments. The relative abundance of T-RFs is given as a percentage of the total peak height. Fragment sizes within the graph indicate the sizes (bp) of the experimental T-RFs by T-RFLP.

Fig S2 Standard curves for 16S rRNA (a), narG(b), nirK(c), nirS(d), qnorB(e) andnosZ(f) assays obtained by plotting the concentration of plasmid DNA versus the threshold cycle. Red multiple sign present the amplification and quantity of soil samples.

Fig S3 Melt curves for 16S rRNA (a), narG(b), nirK(c), nirS(d), qnorB(e) andnosZ(f) genes. The Y-axis is the derivative of the slope of the signal intensity as temperature increases.

Fig S1a narG

Fig S1bqnorB

Fig S1cnosZ

Fig S2 a: 16S rRNA, b: narG, c: nirK, d: nirS, e: qnorB, f: nosZ

Fig S3 a: 16S rRNA, b: narG, c: nirK, d: nirS, e: qnorB, f: nosZ