Nitrogen fixation by native Bradyrhizobia in symbiosis with Lupinus mariae-josephae requires a T3SS encoding a NopE-like effector

Several bradyrhizobial isolates from L. mariae-josephae root nodules [1] contain a type III secretion system (T3SS) within a cluster of about 30 genes. Among those genes, ttsI codes for the transcriptional activator of the system. Mutation of ttsI resulted in the formation of white, non-fixing nodules with the natural legume host, L. mariae-josephae. The T3SS cluster also contains a gene coding for a NopE-like protein. NopE proteins have been demonstrated to be effectors in the Bradyrhizobium-soybean symbiosis [2] and belong to a small group of poorly characterized proteins from plant-associated bacteria that contain one or two autocleavage motifs known as DUF1521 (Schirrmeister et al. 2011). The amino acid sequence of a NopE-like protein in the L. mariae-josephae strain LmjC contains just one autocatalytic motif. This is unlike NopE1 and NopE2 proteins secreted by the T3SS of B. japonicum, that contain two motifs [3]. The autocleavage of LmjC NopE protein was analyzed after expression in E. coli and purification. Two protein fragments of the predicted sizes appeared in the presence of Ca2+, Cu2+, Cd2+, Zn2+ and Mn2+ cations. In contrast, autocleavage did not take place in the presence of Ni2+, Co2+ or Mg2+. Site-directed mutagenesis of the DUF1521 motif in LmjC NopE abolished self-cleavage in vitro. Symbiotic competence of a NopE- mutant with the L. mariae-josephae host was not affected. Possible roles of NopE are discussed.

Identification and functional characterization of Rhizobium leguminosarum bv. viciae genetic systems involved in nickel homeostasis.

A collection of Rhizobium leguminosarum bv. viciae strains isolated from ultramafic and contaminated soils in Italy and Germany, respectively, was analyzed for resistance to nickel and cobalt ions. These assays led to the identification of strain UPM1137, which is able to grow at high concentrations of nickel and cobalt. In order to identify genetic systems involved in the homeostasis to these metals, a random mutagenesis was carried out in UPM1137 by inserting a Tn5-derivative minitransposon. As a result 4313 transconjugants were obtained, being 39 of them (0.90%) unable to grow at 1.5 mM NiCl2. The identification of the transposon insertion site in these mutants showed that the disrupted genes encode proteins belonging to different functional categories, where the secreted and membrane proteins were the most numerous. The analysis of heavy metal resistance and phenotypes in symbiotic and free –living cells will define the contribution of these genes to metal homeostasis.

Azotobacter vinelandii nitrogenase: “Kinetics of nif gene expression and insights into the roles of FdxN and NifQ in FeMo-co biosynthesis”

The Molybdenum-nitrogenase is responsible for most biological nitrogen fixation activity (BNF) in the biosphere. Due to its great agronomical importance, it has been the subject of profound genetic and biochemical studies. The Mo nitrogenase carries at its active site a unique iron-molybdenum cofactor (FeMoco) that consists of an inorganic 7 Fe, 1 Mo, 1 C, 9 S core coordinated to the organic acid homocitrate. Biosynthesis of FeMo-co occurs outside nitrogenase through a complex and highly regulated pathway involving proteins acting as molecular scaffolds, metallocluster carriers or enzymes that provide substrates in appropriate chemical forms. Specific expression regulatory factors tightly control the accumulation levels of all these other components. Insertion of FeMo-co into a P-cluster containing apo-NifDK polypeptide results in nitrogenase reconstitution. Investigation of FeMo-co biosynthesis has uncovered new radical chemistry reactions and new roles for Fe-S clusters in biology.

Biogeochemical indicators of elevated nitrogen deposition in semiarid Mediterranean ecosystems

Nitrogen (N) deposition has doubled the natural N inputs received by ecosystems through biological N fixation and is currently a global problem that is affecting the Mediterranean regions. We evaluated the existing relationships between increased atmospheric N deposition and biogeochemical indicators related to soil chemical factors and cryptogam species across semiarid central, southern, and eastern Spain. The cryptogam species studied were the biocrust-forming species Pleurochaete squarrosa (moss) and Cladonia foliacea (lichen). Sampling sites were chosen in Quercus coccifera (kermes oak) shrublands and Pinus halepensis (Aleppo pine) forests to cover a range of inorganic N deposition representative of the levels found in the Iberian Peninsula (between 4.4 and 8.1 kg N ha(-1) year(-1)). We extended the ambient N deposition gradient by including experimental plots to which N had been added for 3 years at rates of 10, 20, and 50 kg N ha(-1) year(-1). Overall, N deposition (extant plus simulated) increased soil inorganic N availability and caused soil acidification. Nitrogen deposition increased phosphomonoesterase (PME) enzyme activity and PME/nitrate reductase (NR) ratio in both species, whereas the NR activity was reduced only in the moss. Responses of PME and NR activities were attributed to an induced N to phosphorus imbalance and to N saturation, respectively. When only considering the ambient N deposition, soil organic C and N contents were positively related to N deposition, a response driven by pine forests. The PME/NR ratios of the moss were better predictors of N deposition rates than PME or NR activities alone in shrublands, whereas no correlation between N deposition and the lichen physiology was observed. We conclude that integrative physiological measurements, such as PME/NR ratios, measured on sensitive species such as P. squarrosa, can provide useful data for national-scale biomonitoring programs, whereas soil acidification and soil C and N storage could be useful as additional corroborating ecosystem indicators of chronic N pollution.

Análisis genómico y funcional de los sistemas de Quorum Sensing en Rhizobium leguminosarum

Rhizobium leguminosarum (Rl) es una alfa-proteobacteria capaz de establecer una simbiosis diazotrófica con distintas leguminosas. A pesar de la importancia de esta simbiosis en el balance global del ciclo del nitrógeno, muy pocos genomas de rhizobios han sido secuenciados, que aporten nuevos conocimientos relacionados con las características genéticas que contribuyen a importantes procesos simbióticos. Únicamente tres secuencias completas de Rl han sido publicadas: Rl bv. viciae 3841 y dos genomas de Rl bv. trifolii (WSM1325 y WSM2304), ambos simbiontes de trébol. La secuencia genómica de Rlv UPM791 se ha determinado por medio de secuenciación 454. Este genoma tiene un tamaño aproximado de 7.8 Mb, organizado en un cromosoma y 5 replicones extracromosómicos, que incluyen un plásmido simbiótico de 405 kb. Este nuevo genoma se ha analizado en relación a las funciones simbióticas y adaptativas en comparación con los genomas completos de Rlv 3841 y Rl bv. trifolii WSM1325 y WSM2304. Mientras que los plásmidos pUPM791a y b se encuentran conservados, el plásmido simbiótico pUPM791c exhibe un grado de conservación muy bajo comparado con aquellos descritos en las otras cepas de Rl. Uno de los factores implicados en el establecimiento de la simbiosis es el sistema de comunicación intercelular conocido como Quorum Sensing (QS). El análisis del genoma de Rlv UPM791 ha permitido la identificación de dos sistemas tipo LuxRI mediados por señales de tipo N-acyl-homoserina lactonas (AHLs). El análisis mediante HPLC-MS ha permitido asociar las señales C6-HSL, C7-HSL y C8-HSL al sistema rhiRI, codificado en el plásmido simbiótico; mientras que el sistema cinRI, localizado en el cromosoma, produce 3OH-C14:1-HSL. Se ha identificado una tercera sintasa (TraI) codificada en el plásmido simbiótico, pero su regulador correspondiente se encuentra truncado debido a un salto de fase. Adicionalmente, se han encontrado tres reguladores de tipo LuxR-orphan que no presentan una sintasa LuxI asociada. El efecto potencial de las señales tipo AHL se ha estudiado mediante una estrategia de quorum quenching, la cual interfiere con los sistemas de QS de la bacteria. Esta estrategia está basada en la introducción del gen aiiA de Bacillus subtilis, que expresa constitutivamente una enzima lactonasa degradadora de AHLs. Para llevar a cabo el análisis en condiciones simbióticas, se ha desarrollado un sistema de doble marcaje que permite la identificación basado en los marcadores gusA y celB, que codifican para una enzima β–glucuronidasa y una β–galactosidasa termoestable, respectivamente. Los resultados obtenidos indican que Rlv UPM791 predomina sobre la cepa Rlv 3841 para la formación de nódulos en plantas de guisante. La baja estabilidad del plásmido que codifica para aiiA, no ha permitido obtener una conclusión definitiva sobre el efecto de la lactonasa AiiA en competitividad. Con el fin de analizar el significado y la regulación de la producción de moléculas señal tipo AHL, se han generado mutantes defectivos en cada uno de los dos sistemas de QS. Se ha llevado a cabo un análisis detallado sobre la producción de AHLs, formación de biofilm y simbiosis con plantas de guisante, veza y lenteja. El efecto de las deleciones de los genes rhiI y rhiR en Rlv UPM791 es más drástico en ausencia del plásmido pUPM791d. Mutaciones en cinI o cinRIS muestran tanto ausencia de señales, como producción exclusivamente de las de bajo peso molecular, respectivamente, producidas por el sistema rhiRI. Estas mutaciones mostraron un efecto importante en simbiosis. El sistema rhiRI se necesita para un comportamiento simbiótico normal. Además, mutantes cinRIS generaron nódulos blancos e ineficientes, mientras que el mutante cinI fue incapaz de producir nódulos en ninguna de las leguminosas utilizadas. Dicha mutación resultó en la inestabilización del plasmido simbiótico por un mecanismo dependiente de cinI que no ha sido aclarado. En general, los resultados obtenidos indican la existencia de un modelo de regulación dependiente de QS significativamente distinto a los que se han descrito previamente en otras cepas de R. leguminosarum, en las cuales no se había observado ningún fenotipo relevante en simbiosis. La regulación de la producción de AHLs Rlv UPM791 es un proceso complejo que implica genes situados en los plásmidos UPM791c y UPM791d, además de la señal 3-OH-C14:1-HSL. Finalmente, se ha identificado un transportador de tipo RND, homologo a mexAB-oprM de P. aeruginosa e implicado en la extrusión de AHLs de cadena larga. La mutación he dicho transportador no tuvo efectos apreciables sobre la simbiosis. ABSTRACT Rhizobium leguminosarum (Rl) is a soil alpha-proteobacterium that establishes a diazotrophic symbiosis with different legumes. Despite the importance of this symbiosis to the global nitrogen cycling balance, very few rhizobial genomes have been sequenced so far which provide new insights into the genetic features contributing to symbiotically relevant processes. Only three complete sequences of Rl strains have been published: Rl bv. viciae 3841, harboring six plasmids (7.75 Mb) and two Rl bv. trifolii (WSM1325 and WSM2304), both clover symbionts, harboring 5 and 4 plasmids, respectively (7.41 and 6.87 Mb). The genomic sequence of Rlv UPM791 was undertaken by means of 454 sequencing. Illumina and Sanger reads were used to improve the assembly, leading to 17 final contigs. This genome has an estimated size of 7.8 Mb organized in one chromosome and five extrachromosomal replicons, including a 405 kb symbiotic plasmid. Four of these plasmids are already closed, whereas there are still gaps in the smallest one (pUPM791d) due to the presence of insertion elements and repeated sequences, which difficult the assembly. The annotation has been carried out thanks to the Manatee pipeline. This new genome has been analyzed as regarding symbiotic and adaptive functions in comparison to the Rlv 3841 complete genome, and to those from Rl bv. trifolii strains WSM1325 and WSM2304. While plasmids pUPM791a and b are conserved, the symbiotic plasmid pUPM791c exhibited the lowest degree of conservation as compared to those from the other Rl strains. One of the factors involved in the symbiotic process is the intercellular communication system known as Quorum Sensing (QS). This mechanism allows bacteria to carry out diverse biological processes in a coordinate way through the production and detection of extracellular signals that regulate the transcription of different target genes. Analysis of the Rlv UPM791 genome allowed the identification of two LuxRI-like systems mediated by N-acyl-homoserine lactones (AHLs). HPLC-MS analysis allowed the adscription of C6-HSL, C7-HSL and C8-HSL signals to the rhiRI system, encoded in the symbiotic plasmid, whereas the cinRI system, located in the chromosome, produces 3OH-C14:1-HSL, previously described as “bacteriocin small”. A third synthase (TraI) is encoded also in the symbiotic plasmid, but its cognate regulator TraR is not functional due to a fameshift mutation. Three additional LuxR orphans were also found which no associated LuxI-type synthase. The potential effect of AHLs has been studied by means of a quorum quenching approach to interfere with the QS systems of the bacteria. This approach is based upon the introduction into the strains Rl UPM791 and Rl 3841 of the Bacillus subtilis gene aiiA expressing constitutively an AHL-degrading lactonase enzyme which led to virtual absence of AHL even when AiiA-expressing cells were a fraction of the total population. No significant effect of AiiA-mediated AHL removal on competitiveness for growth in solid surface was observed. For analysis under symbiotic conditions we have set up a two-label system to identify nodules produced by two different strains in pea roots, based on the markers gusA and celB, encoding a β–glucuronidase and a thermostable β–galactosidase enzymes, respectively. The results obtained show that Rlv UPM791 outcompetes Rlv 3841 for nodule formation in pea plants, and that the presence of the AiiA plasmid does not significantly affect the relative competitiveness of the two Rlv strains. However, the low stability of the pME6863 plasmid, encoding aiiA, did not lead to a clear conclusion about the AiiA lactonase effect on competitiveness. In order to further analyze the significance and regulation of the production of AHL signal molecules, mutants deficient in each of the two QS systems were constructed. A detailed analysis of the effect of these mutations on AHL production, biofilm formation and symbiosis with pea, vetch and lentil plants has been carried out. The effect of deletions on Rlv UPM791 rhiI and rhiR genes is more pronounced in the absence of plasmid pUPM791d, as no signal is detected in UPM791.1, lacking this plasmid. Mutations in cinI or cinRIS show either no signals, or only the small ones produced by the rhiRI system, suggesting that cinR might be regulating the rhiRI system. These mutations had a strong effect on symbiosis. Analysis of rhi mutants revealed that rhiRI system is required for normal symbiotic performance, as a drastic reduction of symbiotic fitness is observed when rhiI is deleted, and rhiR is essential for nitrogen fixation in the absence of plasmid pUPM791d. Furthermore, cinRIS mutants resulted in white and inefficient nodules, whereas cinI mutant was unable to form nodules on any legume tested. The latter mutation is associated to the instabilization of the symbiotic plasmid through a mechanism still uncovered. Overall, the results obtained indicate the existence of a model of QS-dependent regulation significantly different to that previously described in other R. leguminosarum strains, where no relevant symbiotic phenotype had been observed. The regulation of AHL production in Rlv UPM791 is a complex process involving the symbiotic plasmid (pUPM791c) and the smallest plasmid (pUPM791d), with a key role for the 3-OH-C14:1-HSL signal. Finally, we made a search for potential AHL transporters in Rlv UPM791 genome. These signals diffuse freely across membranes, but in the case of the long-chain AHLs an active efflux system might be required, as it has been described for C12-HSL in the case of Pseudomonas aeruginosa. We have identified a putative AHL transporter of the RND family homologous to P. aeruginosa mexAB-oprM. A mutant strain deficient in this transporter has been generated, and TLC analysis shows absence of 3OH-C14:1-HSL in its supernatant. This deficiency was complemented by the reintroduction of an intact copy of the genes via plasmid transfer. The mutation in mexAB genes had no significant effects on the symbiotic performance of R. leguminosarum bv. viciae.

Nitrogen deposition alters nitrogen cycling and reduces soil carbon content in low-productivity semiarid Mediterranean ecosystems.

Anthropogenic N deposition poses a threat to European Mediterranean ecosystems. We combined data from an extant N deposition gradient (4.3–7.3 kg N ha−1 yr−1) from semiarid areas of Spain and a field experiment in central Spain to evaluate N deposition effects on soil fertility, function and cyanobacteria community. Soil organic N did not increase along the extant gradient. Nitrogen fixation decreased along existing and experimental N deposition gradients, a result possibly related to compositional shifts in soil cyanobacteria community. Net ammonification and nitrification (which dominated N-mineralization) were reduced and increased, respectively, by N fertilization, suggesting alterations in the N cycle. Soil organic C content, C:N ratios and the activity of β-glucosidase decreased along the extant gradient in most locations. Our results suggest that semiarid soils in low-productivity sites are unable to store additional N inputs, and that are also unable to mitigate increasing C emissions when experiencing increased N deposition.

Medicago truncatula Natural Resistance Associated Macrophage Protein1 is required for iron uptake by rhizobia infected nodule cells

Iron is critical for symbiotic nitrogen fixation (SNF) as a key component ofmultiple ferroproteins involved in this biological process. In the model legume Medicago truncatula, iron is delivered by the vasculature to the infection/maturation zone (zone II) of the nodule, where it is released to the apoplast. From there, plasma membrane iron transporters move it into rhizobia-containing cells, where iron is used as the cofactor of multiple plant and rhizobial proteins (e.g. plant leghemoglobin and bacterial nitrogenase). MtNramp1 (Medtr3g088460) is the M. truncatula Natural Resistance-Associated Macrophage Protein family member, with the highest expression levels in roots and nodules. Immunolocalization studies indicate that MtNramp1 is mainly targeted to the plasma membrane. A loss-of-function nramp1 mutant exhibited reduced growth compared with the wild type under symbiotic conditions, but not when fertilized with mineral nitrogen. Nitrogenase activity was low in the mutant, whereas exogenous iron and expression of wild-type MtNramp1 in mutant nodules increased nitrogen fixation to normal levels. These data are consistent with a model in which MtNramp1 is the main transporter responsible for apoplastic iron uptake by rhizobia-infected cells in zone II.

Iron distribution through the developmental stages of Medicago truncatula nodules.

Paramount to symbiotic nitrogen fixation (SNF) is the synthesis of a number of metalloenzymes that use iron as a critical component of their catalytical core. Since this process is carried out by endosymbiotic rhizobia living in legume root nodules, the mechanisms involved in iron delivery to the rhizobia-containing cells are critical for SNF. In order to gain insight into iron transport to the nodule, we have used synchrotron-based X-ray fluorescence to determine the spatio-temporal distribution of this metal in nodules of the legume Medicago truncatula with hitherto unattained sensitivity and resolution. The data support a model in which iron is released from the vasculature into the apoplast of the infection/differentiation zone of the nodule (zone II). The infected cell subsequently takes up this apoplastic iron and delivers it to the symbiosome and the secretory system to synthesize ferroproteins. Upon senescence, iron is relocated to the vasculature to be reused by the shoot. These observations highlight the important role of yet to be discovered metal transporters in iron compartmentalization in the nodule and in the recovery of an essential and scarce nutrient for flowering and seed production.

MtNramp1 mediates iron import in rhizobia-infected Medicago truncatula cells.

Symbiotic nitrogen fixation is a process that requires relatively high quantities of iron provided by the host legume. Using synchrotron-based X-ray fluorescence, we have determined that this iron is released from the vasculature into the apoplast of zone II of M. truncatula nodules. This overlaps with the distribution of MtNramp1, a plasma membrane iron importer. The importance of MtNramp1 in iron transport for nitrogen fixation is indicated by the 60% reduction of nitrogenase activity observed in knock-down lines, most likely due to deficient incorporation of this essential metal cofactor at the necessary levels.

Functional characterization of Rhizobium leguminosarum bv. viciae DmeRF, a cation diffusion facilitator system involved in nickel and cobalt resistance

In prokaryotes, nickel is an essential element participating in the structure of enzymes involved in multiple cellular processes. Nickel transport is a challenge for microorganisms since, although essential, high levels of this metal inside the cell are toxic. For this reason, bacteria have developed high-affinity nickel transporters as well as nickel-specific detoxification systems. Ultramafic soils, and soils contaminated with heavy metals are excellent sources of nickel resistant bacteria. Molecular analysis of strains isolated in the habitats has revealed novel genetic systems involved in adaptation to such hostile conditions.

Identification and analysis of a nickel-inducible cation diffusion facilitator-NiCoT combined system in Rhizobium leguminosarum bv. viciae

Nickel, like other transition metals, can be toxic to cells even at moderate concentration (low microM range) by displacing essential metals from their native binding sites or by generating reactive oxygen species that cause oxidative DNA damage. For this reason, cells have evolved mechanisms to deal with excess nickel. Efflux systems include members of the Resistance-Nodulation-cell Division (RND) protein family, P-type ATPases, cation diffusion facilitators (CDF) and other resistance factors. Nickel-specific exporters have been characterized in Cupravidus metallidurans, Helicobacter pylori, Achromobacter xylosoxidans, Serratia marcenses and Escherichia coli.

Metagenomic Anlaysis of microsymbiont selection by the legume plant host

Rhizobium leguminosarum bv.viciae is able to establish nitrogen-fixing symbioses with legumes of the genera Pisum, Lens, Lathyrus and Vicia. Classic studies using trap plants (Laguerre et al., Young et al.) provided evidence that different plant hosts are able to select different rhizobial genotypes among those available in a given soil. However, these studies were necessarily limited by the paucity of relevant biodiversity markers. We have now reappraised this problem with the help of genomic tools. A well-characterized agricultural soil (INRA Bretennieres) was used as source of rhizobia. Plants of Pisum sativum, Lens culinaris, Vicia sativa and V. faba were used as traps. Isolates from 100 nodules were pooled, and DNA from each pool was sequenced (BGI-Hong Kong; Illumina Hiseq 2000, 500 bp PE libraries, 100 bp reads, 12 Mreads). Reads were quality filtered (FastQC, Trimmomatic), mapped against reference R. leguminosarum genomes (Bowtie2, Samtools), and visualized (IGV). An important fraction of the filtered reads were not recruited by reference genomes, suggesting that plant isolates contain genes that are not present in the reference genomes. For this study, we focused on three conserved genomic regions: 16S-23S rDNA, atpD and nodDABC, and a Single Nucleotide Polymorphism (SNP) analysis was carried out with meta / multigenomes from each plant. Although the level of polymorphism varied (lowest in the rRNA region), polymorphic sites could be identified that define the specific soil population vs. reference genomes. More importantly, a plant-specific SNP distribution was observed. This could be confirmed with many other regions extracted from the reference genomes (data not shown). Our results confirm at the genomic level previous observations regarding plant selection of specific genotypes. We expect that further, ongoing comparative studies on differential meta / multigenomic sequences will identify specific gene components of the plant-selected genotypes

Population genomics of host specificity in Rhizobium leguminosarum bv. viciae

Legumes establish a root-nodule symbiosis with soil bacteria collectively known as rhizobia. This symbiosis allows legumes to benefit from the nitrogen fixation capabilities of rhizobia and thus to grow in the absence of any fixed nitrogen source. This is especially relevant for Agriculture, where intensive plant growth depletes soils of useable, fixed nitrogen sources. One of the main features of the root nodule symbiosis is its specificity. Different rhizobia are able to nodulate different legumes. Rhizobium leguminosarum bv. viciae is able to establish an effective symbiosis with four different plant genera (Pisum, Lens, Vicia, Lathyrus), and any given isolate will nodulate any of the four plant genera. A population genomics study with rhizobia isolated from P. sativum, L. culinaris, V. sativa or V. faba, all originating in the same soil, showed that plants select specific genotypes from those available in that soil. This was demonstrated at the genome-wide level, but also for specific genes. Accelerated mesocosm studies with successive plant cultures provided additional evidence on this plant selection and on the nature of the genotypes selected. Finally, representatives from the major rhizobial genotypes isolated from these plants allowed characterization of the size and nature of the respective pangenome and specific genome compartments. These were compared to the different genotypes ?symbiotic and non-symbiotic?present in rhizobial populations isolated directly from the soil without plant intervention.

Estimation of nitrogen budgets for contrasting catchments at the landscape scale.

A comprehensive assessment of nitrogen (N) flows at the landscape scale is fundamental to understand spatial interactions in the N cascade and to inform the development of locally optimised N management strategies. To explore these interactions, complete N budgets were estimated for two contrasting hydrological catchments (dominated by agricultural grassland vs. semi-natural peat-dominated moorland), forming part of an intensively studied landscape in southern Scotland. Local scale atmospheric dispersion modelling and detailed farm and field inventories provided high resolution estimations of input fluxes. Direct agricultural inputs (i.e. grazing excreta, N2 fixation, organic and synthetic fertiliser) accounted for most of the catchment N inputs, representing 82% in the grassland and 62% in the moorland catchment, while atmospheric deposition made a significant contribution, particularly in the moorland catchment, contributing 38% of the N inputs. The estimated catchment N budgets highlighted areas of key uncertainty, particularly N2 exchange and stream N export. The resulting N balances suggest that the study catchments have a limited capacity to store N within soils, vegetation and groundwater. The "catchment N retention", i.e. the amount of N which is either stored within the catchment or lost through atmospheric emissions, was estimated to be 13% of the net anthropogenic input in the moorland and 61% in the grassland catchment. These values contrast with regional scale estimates: Catchment retentions of net anthropogenic input estimated within Europe at the regional scale range from 50% to 90%, with an average of 82% (Billen et al., 2011). This study emphasises the need for detailed budget analyses to identify the N status of European landscapes.

Evaluation of cover crops based on characteristics of nitrogen use, and aerial and root growth

La caracterización de los cultivos cubierta (cover crops) puede permitir comparar la idoneidad de diferentes especies para proporcionar servicios ecológicos como el control de la erosión, el reciclado de nutrientes o la producción de forrajes. En este trabajo se estudiaron bajo condiciones de campo diferentes técnicas para caracterizar el dosel vegetal con objeto de establecer una metodología para medir y comparar las arquitecturas de los cultivos cubierta más comunes. Se estableció un ensayo de campo en Madrid (España central) para determinar la relación entre el índice de área foliar (LAI) y la cobertura del suelo (GC) para un cultivo de gramínea, uno de leguminosa y uno de crucífera. Para ello se sembraron doce parcelas con cebada (Hordeum vulgare L.), veza (Vicia sativa L.), y colza (Brassica napus L.). En 10 fechas de muestreo se midieron el LAI (con estimaciones directas y del LAI-2000), la fracción interceptada de la radiación fotosintéticamente activa (FIPAR) y la GC. Un experimento de campo de dos años (Octubre-Abril) se estableció en la misma localización para evaluar diferentes especies (Hordeum vulgare L., Secale cereale L., x Triticosecale Whim, Sinapis alba L., Vicia sativa L.) y cultivares (20) en relación con su idoneidad para ser usadas como cultivos cubierta. La GC se monitorizó mediante análisis de imágenes digitales con 21 y 22 muestreos, y la biomasa se midió 8 y 10 veces, respectivamente para cada año. Un modelo de Gompertz caracterizó la cobertura del suelo hasta el decaimiento observado tras las heladas, mientras que la biomasa se ajustó a ecuaciones de Gompertz, logísticas y lineales-exponenciales. Al final del experimento se determinaron el C, el N y el contenido en fibra (neutrodetergente, ácidodetergente y lignina), así como el N fijado por las leguminosas. Se aplicó el análisis de decisión multicriterio (MCDA) con objeto de obtener un ranking de especies y cultivares de acuerdo con su idoneidad para actuar como cultivos cubierta en cuatro modalidades diferentes: cultivo de cobertura, cultivo captura, abono verde y forraje. Las asociaciones de cultivos leguminosas con no leguminosas pueden afectar al crecimiento radicular y a la absorción de N de ambos componentes de la mezcla. El conocimiento de cómo los sistemas radiculares específicos afectan al crecimiento individual de las especies es útil para entender las interacciones en las asociaciones, así como para planificar estrategias de cultivos cubierta. En un tercer ensayo se combinaron estudios en rhizotrones con extracción de raíces e identificación de especies por microscopía, así como con estudios de crecimiento, absorción de N y 15N en capas profundas del suelo. Las interacciones entre raíces en su crecimiento y en el aprovisionamiento de N se estudiaron para dos de los cultivares mejor valorados en el estudio previo: uno de cebada (Hordeum vulgare L. cv. Hispanic) y otro de veza (Vicia sativa L. cv. Aitana). Se añadió N en dosis de 0 (N0), 50 (N1) y 150 (N2) kg N ha-1. Como resultados del primer estudio, se ajustaron correctamente modelos lineales y cuadráticos a la relación entre la GC y el LAI para todos los cultivos, pero en la gramínea alcanzaron una meseta para un LAI>4. Antes de alcanzar la cobertura total, la pendiente de la relación lineal entre ambas variables se situó en un rango entre 0.025 y 0.030. Las lecturas del LAI-2000 estuvieron correlacionadas linealmente con el LAI, aunque con tendencia a la sobreestimación. Las correcciones basadas en el efecto de aglutinación redujeron el error cuadrático medio del LAI estimado por el LAI-2000 desde 1.2 hasta 0.5 para la crucífera y la leguminosa, no siendo efectivas para la cebada. Esto determinó que para los siguientes estudios se midieran únicamente la GC y la biomasa. En el segundo experimento, las gramíneas alcanzaron la mayor cobertura del suelo (83-99%) y la mayor biomasa (1226-1928 g m-2) al final del mismo. Con la mayor relación C/N (27-39) y contenido en fibra digestible (53-60%) y la menor calidad de residuo (~68%). La mostaza presentó elevadas GC, biomasa y absorción de N en el año más templado en similitud con las gramíneas, aunque escasa calidad como forraje en ambos años. La veza presentó la menor absorción de N (2.4-0.7 g N m-2) debido a la fijación de N (9.8-1.6 g N m-2) y escasa acumulación de N. El tiempo térmico hasta alcanzar el 30% de GC constituyó un buen indicador de especies de rápida cubrición. La cuantificación de las variables permitió hallar variabilidad entre las especies y proporcionó información para posteriores decisiones sobre la selección y manejo de los cultivos cubierta. La agregación de dichas variables a través de funciones de utilidad permitió confeccionar rankings de especies y cultivares para cada uso. Las gramíneas fueron las más indicadas para los usos de cultivo de cobertura, cultivo captura y forraje, mientras que las vezas fueron las mejor como abono verde. La mostaza alcanzó altos valores como cultivo de cobertura y captura en el primer año, pero el segundo decayó debido a su pobre actuación en los inviernos fríos. Hispanic fue el mejor cultivar de cebada como cultivo de cobertura y captura, mientras que Albacete como forraje. El triticale Titania alcanzó la posición más alta como cultiva de cobertura, captura y forraje. Las vezas Aitana y BGE014897 mostraron buenas aptitudes como abono verde y cultivo captura. El MCDA permitió la comparación entre especies y cultivares proporcionando información relevante para la selección y manejo de cultivos cubierta. En el estudio en rhizotrones tanto la mezcla de especies como la cebada alcanzaron mayor intensidad de raíces (RI) y profundidad (RD) que la veza, con valores alrededor de 150 cruces m-1 y 1.4 m respectivamente, comparados con 50 cruces m-1 y 0.9 m para la veza. En las capas más profundas del suelo, la asociación de cultivos mostró valores de RI ligeramente mayores que la cebada en monocultivo. La cebada y la asociación obtuvieron mayores valores de densidad de raíces (RLD) (200-600 m m-3) que la veza (25-130) entre 0.8 y 1.2 m de profundidad. Los niveles de N no mostraron efectos claros en RI, RD ó RLD, sin embargo, el incremento de N favoreció la proliferación de raíces de veza en la asociación en capas profundas del suelo, con un ratio cebada/veza situado entre 25 a N0 y 5 a N2. La absorción de N de la cebada se incrementó en la asociación a expensas de la veza (de ~100 a 200 mg planta-1). Las raíces de cebada en la asociación absorbieron también más nitrógeno marcado de las capas profundas del suelo (0.6 mg 15N planta-1) que en el monocultivo (0.3 mg 15N planta-1). ABSTRACT Cover crop characterization may allow comparing the suitability of different species to provide ecological services such as erosion control, nutrient recycling or fodder production. Different techniques to characterize plant canopy were studied under field conditions in order to establish a methodology for measuring and comparing cover crops canopies. A field trial was established in Madrid (central Spain) to determine the relationship between leaf area index (LAI) and ground cover (GC) in a grass, a legume and a crucifer crop. Twelve plots were sown with either barley (Hordeum vulgare L.), vetch (Vicia sativa L.), or rape (Brassica napus L.). On 10 sampling dates the LAI (both direct and LAI-2000 estimations), fraction intercepted of photosynthetically active radiation (FIPAR) and GC were measured. A two-year field experiment (October-April) was established in the same location to evaluate different species (Hordeum vulgare L., Secale cereale L., x Triticosecale Whim, Sinapis alba L., Vicia sativa L.) and cultivars (20) according to their suitability to be used as cover crops. GC was monitored through digital image analysis with 21 and 22 samples, and biomass measured 8 and 10 times, respectively for each season. A Gompertz model characterized ground cover until the decay observed after frosts, while biomass was fitted to Gompertz, logistic and linear-exponential equations. At the end of the experiment C, N, and fiber (neutral detergent, acid and lignin) contents, and the N fixed by the legumes were determined. Multicriteria decision analysis (MCDA) was applied in order to rank the species and cultivars according to their suitability to perform as cover crops in four different modalities: cover crop, catch crop, green manure and fodder. Intercropping legumes and non-legumes may affect the root growth and N uptake of both components in the mixture. The knowledge of how specific root systems affect the growth of the individual species is useful for understanding the interactions in intercrops as well as for planning cover cropping strategies. In a third trial rhizotron studies were combined with root extraction and species identification by microscopy and with studies of growth, N uptake and 15N uptake from deeper soil layers. The root interactions of root growth and N foraging were studied for two of the best ranked cultivars in the previous study: a barley (Hordeum vulgare L. cv. Hispanic) and a vetch (Vicia sativa L. cv. Aitana). N was added at 0 (N0), 50 (N1) and 150 (N2) kg N ha-1. As a result, linear and quadratic models fitted to the relationship between the GC and LAI for all of the crops, but they reached a plateau in the grass when the LAI > 4. Before reaching full cover, the slope of the linear relationship between both variables was within the range of 0.025 to 0.030. The LAI-2000 readings were linearly correlated with the LAI but they tended to overestimation. Corrections based on the clumping effect reduced the root mean square error of the estimated LAI from the LAI-2000 readings from 1.2 to less than 0.50 for the crucifer and the legume, but were not effective for barley. This determined that in the following studies only the GC and biomass were measured. In the second experiment, the grasses reached the highest ground cover (83- 99%) and biomass (1226-1928 g/m2) at the end of the experiment. The grasses had the highest C/N ratio (27-39) and dietary fiber (53-60%) and the lowest residue quality (~68%). The mustard presented high GC, biomass and N uptake in the warmer year with similarity to grasses, but low fodder capability in both years. The vetch presented the lowest N uptake (2.4-0.7 g N/m2) due to N fixation (9.8-1.6 g N/m2) and low biomass accumulation. The thermal time until reaching 30% ground cover was a good indicator of early coverage species. Variable quantification allowed finding variability among the species and provided information for further decisions involving cover crops selection and management. Aggregation of these variables through utility functions allowed ranking species and cultivars for each usage. Grasses were the most suitable for the cover crop, catch crop and fodder uses, while the vetches were the best as green manures. The mustard attained high ranks as cover and catch crop the first season, but the second decayed due to low performance in cold winters. Hispanic was the most suitable barley cultivar as cover and catch crop, and Albacete as fodder. The triticale Titania attained the highest rank as cover and catch crop and fodder. Vetches Aitana and BGE014897 showed good aptitudes as green manures and catch crops. MCDA allowed comparison among species and cultivars and might provide relevant information for cover crops selection and management. In the rhizotron study the intercrop and the barley attained slightly higher root intensity (RI) and root depth (RD) than the vetch, with values around 150 crosses m-1 and 1.4 m respectively, compared to 50 crosses m-1 and 0.9 m for the vetch. At deep soil layers, intercropping showed slightly larger RI values compared to the sole cropped barley. The barley and the intercropping had larger root length density (RLD) values (200-600 m m-3) than the vetch (25-130) at 0.8-1.2 m depth. The topsoil N supply did not show a clear effect on the RI, RD or RLD; however increasing topsoil N favored the proliferation of vetch roots in the intercropping at deep soil layers, with the barley/vetch root ratio ranging from 25 at N0 to 5 at N2. The N uptake of the barley was enhanced in the intercropping at the expense of the vetch (from ~100 mg plant-1 to 200). The intercropped barley roots took up more labeled nitrogen (0.6 mg 15N plant-1) than the sole-cropped barley roots (0.3 mg 15N plant-1) from deep layers.