Plants are associated with a complex microbiota that contributes to
nutrient acquisition, plant growth, and plant defense. Nitrogen-fixing
microbial associations are efficient and well characterized in legumes but
are limited in cereals, including maize. We studied an indigenous landrace
of maize grown in nitrogen-depleted soils in the Sierra Mixe region of
Oaxaca, Mexico. This landrace is characterized by the extensive development
of aerial roots that secrete a carbohydrate-rich mucilage. Analysis of the
mucilage microbiota indicated that it was enriched in taxa for which many
known species are diazotrophic, was enriched for homologs of genes encoding
nitrogenase subunits, and harbored active nitrogenase activity as assessed
by acetylene reduction and 15N2 incorporation assays. Field experiments in
Sierra Mixe using 15N natural abundance or 15N-enrichment assessments over
5 years indicated that atmospheric nitrogen fixation contributed 29%–82% of
the nitrogen nutrition of Sierra Mixe maize.
Author summary
Nitrogen is an essential nutrient for plants, and for many nonlegume crops,
the requirement for nitrogen is primarily met by the use of inorganic
fertilizers. These fertilizers are produced from fossil fuel by
energy-intensive processes that are estimated to use 1% to 2% of the total
global energy supply and produce an equivalent share of greenhouse gases.
Because maize (*Zea mays L*.) is a significant recipient of nitrogen
fertilization, a research goal for decades has been to identify or engineer
mechanisms for biological fixation of atmospheric nitrogen in association
with this crop. We hypothesized that isolated indigenous landraces of maize
grown using traditional practices with little or no fertilizer might have
evolved strategies to improve plant performance under low-nitrogen nutrient
conditions. Here, we show that for one such maize landrace grown in
nitrogen-depleted fields near Oaxaca, Mexico, 29%–82% of the plant nitrogen
is derived from atmospheric nitrogen. High levels of nitrogen fixation are
supported, at least in part, by the abundant production of a sugar-rich
mucilage associated with aerial roots that provides a home to a complex
nitrogen-fixing microbiome.
Figures
[image: Table 2]
[image: Fig 1]
[image: Fig 2]
[image: Fig 3]
[image: Fig 4]
[image: Fig 5]
[image: Table 1]
[image: Table 2]
[image: Fig 1]
[image: Fig 2]
[image: Fig 3]
*Citation: *Van Deynze A, Zamora P, Delaux P-M, Heitmann C, Jayaraman D,
Rajasekar S, et al. (2018) Nitrogen fixation in a landrace of maize is
supported by a mucilage-associated diazotrophic microbiota. PLoS Biol
16(8): e2006352. https://doi.org/10.1371/journal.pbio.2006352
*Academic Editor: *Eric Kemen, University of Tübingen, Germany
*Funding: *Mars, Incorporated http://www.mars.com/global. The research was
funded by an unrestricted gift and a grant to ABB. The funder had no role
in study design, data collection and analysis, decision to publish, or
preparation of the manuscript. The research was funded by grants to ABB
from BioN2, Incorporated. The funder had no role in study design, data
collection and analysis, decision to publish, or preparation of the
manuscript. Preliminary research was supported by a Dissertation Research
Grant from UC Mexus to KS. The funder had no role in study design, data
collection and analysis, decision to publish, or preparation of the
manuscript.
*Competing interests: * Howard-Yana Shapiro is affiliated with Mars,
Incorporated, one of the research sponsors. Author Cristobal Heitmann was
unable to confirm authorship or contributions himself, and this was carried
out collectively by the other co-authors.
*Abbreviations: *%Ndfa, percent of nitrogen derived from the atmosphere;
%Ndiff, percent total nitrogen difference; ARA, acetylene reduction assay;
CCRC, Complex Carbohydrate Center Research Center; GC/MS, gas
chromatography/mass spectrometry; IRMS, isotope-ratio mass spectrometry;
NPGS, National Plant Germplasm System; PCoA, principal component analysis;
TMS, trimethylsilyl
Introduction
Triplett suggested that it may be interesting to survey primitive maize
landraces from the areas of maize origin to identify maize diazotrophic
endophytes [5
<http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.2006352#pbio.2006352.ref005>].
Estrada and colleagues [23
<http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.2006352#pbio.2006352.ref023>]
followed this suggestion and examined a landrace of maize in the Sierra
Mixe region of Oaxaca, Mexico, and isolated a nitrogen-fixing endophyte
from the resident maize landrace. The isolate was tentatively identified as
a new species of *Burkholderia*, but the contribution of atmospheric
dinitrogen (N2) to the nitrogen economy of the plant was not tested. This
group also reported the isolation of a similar endophyte from field-grown
teosinte plants and speculated that the *Burkholderia* strain might have
formed a primitive symbiosis with teosinte that persisted during
domestication of maize.
Fig 1. Physiological features of Sierra Mixe maize.
(A) The transition between juvenile and adult phases (black arrow) occurs 5
weeks after planting in Sierra Mixe maize (black bars) and in the tall
maize heirloom Hickory King (gray bars). (B) Number of aerial roots
observed on Sierra Mixe maize and Hickory King after 14 weeks of growth in
the field in Madison, United States of America. Error bars represent
standard errors; an asterisk indicates a significant difference between
Sierra Mixe maize and Hickory King (Student *t* test, *P* < 0.01). (Data at
DOI: 10.6084/m9.figshare.6534545
<https://doi.org/10.6084/m9.figshare.6534545>).
The aerial roots of Sierra Mixe maize (left) secrete large quantities of
mucilage between 3 and 6 months after planting. The mucilage is
carbohydrate rich, with the composition dominated by arabinose, fucose, and
galactose (side panel).
Fig 3. DNA sequencing–based characterization of the microbiome of Sierra
Mixe maize.
(A) Samples clustered by PCoA on Bray-Curtis dissimilarity distance matrix.
Bacterial communities were assessed using PCR amplification and sequencing
of rRNA genes. Each point corresponds to an individual sample. Permanova
tests run using *adonis* in *vegan* revealed mucilage samples were
statistically distinct (e.g., mucilage versus rhizosphere *P* = 0.002,
mucilage versus roots *P* = 0.03, mucilage versus aerial roots *P* = 0.03).
(B) Metagenomic sequencing–based analysis of homologs of core *nif* genes (
*nifH*, *nifD*, *nifE*, *nifK*, *nifN*, and *nifB*) and alternate
nitrogenase (*anfG*/*vnfG*). Metagenomic samples were searched for homologs
by mapping reads on reference *nif* trees. The number of hits was
normalized to an estimate of the number of bacterial genes in the
metagenomic sample (measured using the number of hits to the RecA hidden
Markov model). *The lowest abundance of a core *nif* gene in mucilage and
rhizosphere libraries. **Only core *nif* gene hit in stem library. anfG and
vnfG alternate nitrogenase were observed only in mucilage library. PCoA,
principal component analysis.
Fig 4. Nitrogenase and N2 fixation activity in mucilage produced by Sierra
Mixe maize.
(A) Mucilage of various Sierra Mixe maize lines collected in Sierra Mixe or
field-grown plants in Madison, USA, and of teosinte display strong
acetylene reduction activity. (-, no acetylene; +, 10% acetylene).
Asterisks indicate significant differences (**P* < 0.05; ***P* < 0.01,
Mann-Whitney test). (B) Nitrogen fixation in Sierra Mixe maize mucilage by
15N2 assimilation. Mucilage collected from Sierra Mixe maize grown in
Sierra Mixe was incubated in gas-tight vials filled with 15N2 or 14N2 gas
for 70 hours at 37 °C. 15N (atom % excess) was determined by IRMS. (C and
D) *H*. *seropedicae* and *A*. *brasilense* display acetylene reduction
activity when added to nonfixing mucilage, whereas the same mucilage
supplemented with sterile medium (-) or the same bacteria without mucilage
(-) do not. (E) Oxygen concentration at 8 mm inside of the mucilage. Means
and standard errors are shown. Different letters indicate statistically
supported groups (Kruskal-Wallis test). (Data at DOI:
10.6084/m9.figshare.6534545 <https://doi.org/10.6084/m9.figshare.6534545>).
IRMS, isotope-ratio mass spectrometry.
As with Sierra Mixe maize, a wild relative, *Z*. *mays* ssp. *mexicana*
(teosinte), also produced extensive aerial roots (S5 Fig
<http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.2006352#pbio.2006352.s005>)
but much smaller amounts of secreted mucilage. To test the ability of the
teosinte mucilage to support nitrogen fixation, we collected mucilage from
several plants of teosinte and measured endogenous nitrogenase activity
using ARA. Acetylene reduction was readily observed in teosinte mucilage (Fig
4A
<http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.2006352#pbio-2006352-g004>),
suggesting that production of mucilage that supports nitrogen fixation by
an associated nitrogen-fixing microbiota may be an ancient trait of maize
and potentially introgressed from *Z*. *mays* ssp. *mexicana* into the
Sierra Mixe landrace postdomestication.
To assess the mucilage characteristics that support nitrogen fixation, we
tested the ability of 2 phylogenetically distinct nitrogen-fixing bacteria
to reduce acetylene when inoculated in the mucilage collected from aerial
roots of Sierra Mixe maize. Before the experiment, the mucilage was frozen
for 2 weeks at –80 °C and thawed, which abolished endogenous nitrogenase
activity (Fig 4C
<http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.2006352#pbio-2006352-g004>).
Two nitrogen-fixing bacteria, *H*. *seropedicae*, and *A*. *brasilense*,
showed readily detectable ARA activity when added to the mucilage (Fig 4C
and 4D
<http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.2006352#pbio-2006352-g004>).
Bacterial nitrogenase is O2-sensitive and needs to be protected by a
low-oxygen (<5%) environment or physiological protective mechanisms, as
well as an abundant carbon source to derive energy for this process [35
<http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.2006352#pbio.2006352.ref035>].
To determine if mucilage could fulfill these requirements, we measured the
free-oxygen concentration in the mucilage of Sierra Mixe maize and teosinte
at the depth of 8 mm and found it to be <5%, indicating that the mucilage
can provide a microaerobic environment compatible with nitrogen fixation
for these bacteria [36
<http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.2006352#pbio.2006352.ref036>]
(Fig 4D
<http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.2006352#pbio-2006352-g004>).
Mucilage is also comprised of complex sugars that may be catabolized to
provide free sugars—mainly arabinose, fucose, and galactose—capable of
supporting bacterial growth and nitrogen fixation. To determine whether
these properties of the mucilage are sufficient to support nitrogen
fixation, we created an artificial medium mimicking these mucilage
properties by using a low-N medium, solidified with 0.2% agar, that reduced
oxygen concentration to levels almost as low as those found in the mucilage
(S6A Fig
<http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.2006352#pbio.2006352.s006>)
and supplemented with a mix of free sugars corresponding to the composition
of the fully hydrolyzed mucilage carbohydrates. *H*. *seropedicae*, *A*.
*brasilense*, and *Burkholderia unamae* (S6B, S6C and S6D Fig
<http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.2006352#pbio.2006352.s006>)
showed significant ARA activity in this reconstructed mucilage, indicating
that the low O2 and free sugars provided by the aerial root mucilage are
sufficient to support nitrogen fixation by these diazotrophs.
Based on ARA, we can conclude that mucilage from Sierra Mixe maize harbors
native diazotrophs and can also support the N2-fixing activity of the
exogenously inoculated diazotrophs, *H*. *seropedicae*, *A*. *brasilense*,
and *B*. *unamae*. However, these data did not demonstrate that the aerial
roots had the capacity to take up and assimilate the fixed N. To test
whether atmospheric N2 that was fixed by mucilage-associated diazotrophs
could be transferred to and utilized by the Sierra Mixe maize, a more
direct 15N2 gas–enrichment experiment was used. Aerial roots, along with
their generated mucilage, inoculated with *A*. *brasilense* Sp7 exhibited
significant 15N2 gas incorporation in comparison to the 14N2 gas-treated
roots (Fig 5
<http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.2006352#pbio-2006352-g005>),
and isotope-ratio mass spectrometry (IRMS) analysis confirmed significant
enrichment of 15N in chlorophyll (converted to pheophytin for analysis) of
these roots compared to the negative controls (Fig 5
<http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.2006352#pbio-2006352-g005>
).
[image: thumbnail]
<http://journals.plos.org/plosbiology/article/figure/image?size=medium&id=info:doi/10.1371/journal.pbio.2006352.g005>
Download:
Fig 5. Analysis of Sierra Mixe maize samples for 15N2 enrichment in
mucilage, aerial roots, and pheophytin from aerial roots.
Mucilage was generated from aerial roots and inoculated with *A*.
*brasilense* Sp7. 15N2 gas was pumped in and, after incubation mucilage,
was separated from the aerial roots. Mucilage alone and aerial roots alone
were subjected to analysis by IRMS. Results revealed a significant
enrichment of 15N2 in mucilage alone and aerial roots alone (left y-axis).
Since the inoculated aerial roots may contain *A*. *brasilense* Sp7
attached to the surface, we extracted pheophytin from these aerial roots to
test 15N2 incorporation in pheophytin. Results revealed a significant
enrichment of 15N2 in these aerial roots, indicating that aerial roots are
indeed the sites for transfer of fixed nitrogen to the plants (right
y-axis). 14N2 samples were used as negative controls. *n* = 4 (aerial roots
and mucilage), *n* = 3 (pheophytin). The asterisk (*) indicates a
statistically significant difference (*p* < 0.05). (Data at DOI:
10.6084/m9.figshare.6534545 <https://doi.org/10.6084/m9.figshare.6534545>)
IRMS, isotope-ratio mass spectrometry.
(A) δ15N values from Sierra Mixe maize, a conventional maize variety (Maiz
Blanco Conasupo), and reference plants grown in Sierra Mixe, 3 months after
planting. Values are given as mean and standard deviation. Statistical
comparisons were made between the means of the reference plants, Maiz
Blanco Conasupo, and *Z*. *mays* S. Mixe, using Student *t* tests (*p* <
0.05). Different letters indicate statistically supported groups. (B) δ15N
values from Sierra Mixe maize plants grown in Fields 1 and 2 in Sierra Mixe
during 2011 and 2012 at 2, 3, 4, 5, and 6 months after planting. Values are
given as mean and standard deviation. Statistical comparisons were made
between the reference plants mean and *Z*. *mays* S. Mixe means at each
time point using Student *t* tests (*p* < 0.05). Different letters indicate
statistically supported groups. Reference plants are listed in S3 Table
<http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.2006352#pbio.2006352.s010>.
(Data at DOI: 10.6084/m9.figshare.6534545
<https://doi.org/10.6084/m9.figshare.6534545>).
The method of using natural abundance 15N and other species as reference
plants to calculate %Ndfa is potentially limiting because of differences in
root and shoot growth and phenology of reference and test plants and the
limited range of δ15N found in soils. An alternative method, 15N
enrichment, is similar to the natural abundance methods but enriches soil 15N
by the addition of 15N fertilizer, thereby increasing the difference
between soil and atmospheric δ15N and assay sensitivity. Several direct
comparisons have indicated that both methods can give comparable but
slightly different results [38
<http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.2006352#pbio.2006352.ref038>
,40
<http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.2006352#pbio.2006352.ref040>
].
Another method of estimating N2 fixation is the “N Difference” method,
which determines the difference between the total N content of an N2-fixing
plant and the total N content of a reference nonfixing plant. Total N is
calculated by multiplying total N content (%) in a specific plant sample
and the total biomass (kg/ha or kg) produced by the plant. The %Ndiff is
calculated as described in Materials and methods.
In 2016 and 2017, we assessed atmospheric nitrogen fixation using the
15N-isotope-enrichment
method (1%–10% enrichment) at 3 vegetative growth stages—V9, V12, and
Tassel—[41
<http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.2006352#pbio.2006352.ref041>]
in a random complete block design (5 replicates) in 3 low-N Sierra Mixe
fields and at Tassel stage in 2017 in the same fields using the same design
(Table 2
<http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.2006352#pbio-2006352-t002>).
Field 3 (with a history of 0–1 year of maize production) yielded
significant differences in Atom% 15N only for SM2 at the Tassel stage in
2016, but shoot N was significantly different in both 2016 and 2017. In
2016, Sierra Mixe maize landrace varieties exhibited significantly lower
Atom%15N levels than the reference plants in Field 4, (with a history of
1–2 years of maize production) at Tassel and at V9 in 2016, and at Tassel
in both years for shoot N. In both 2016 and 2017, Sierra Mixe maize
landraces in Field 5 (over 3 years of continuous maize production)
exhibited significantly lower Atom%15N and shoot N at all stages sampled,
indicating a significant level of atmospheric nitrogen fixation in those
experiments. The calculated %Ndfa ranged between 31% and 55%, and Ndiff
ranged from 29%–82%. The correlation between Ndfa and Ndiff was 0.55 and
0.44 (*P* < 0.01) across locations in 2016 and 2017, respectively.
Significant measures of N2 fixation were detected in 4 of 6 experiments by
Atom% δ15N N determinations (Ndfa) and in 6 of 6 experiments by Ndiff
determinations. Root and shoots exhibited significant differences in
biomass, height, and stem diameter between control hybrids and local
landraces (S4 Table
<http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.2006352#pbio.2006352.s011>).
Soil analyses (S5 Table
<http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.2006352#pbio.2006352.s012>)
showed that all sites were depleted in soil nitrogen yet produced a crop
greater than 2,000 kg/ha. It is possible that differences between
microbiota associated with the 3 fields in different stages of crop
rotation (0 to >4 years of continuous maize production) account for
differences observed among fields and between years, but further research
is needed to answer this question.
[image: thumbnail]
<http://journals.plos.org/plosbiology/article/figure/image?size=medium&id=info:doi/10.1371/journal.pbio.2006352.t002>
Download:
Table 2. 15N-isotope-enrichment determinations in field trials.
Percent Ndfa and Ndiff were calculated for Sierra Mixe varieties when Atom%
15N excess or Shoot N for the whole plant were significantly different from
reference varieties as assessed by ANOVA and single-degree-of-freedom
contrasts (*p* = 0.05), respectively. Values followed by asterisks are
significantly different from the reference varieties based on
single-degree-of-freedom contrasts (*p* < 0.05). Percent Ndfa and Ndiff
were not calculated for reference (dashes). (Data at DOI:
10.6084/m9.figshare.6534545 <https://doi.org/10.6084/m9.figshare.6534545>).
We have demonstrated that the mucilage associated with the aerial roots of
Sierra Mixe maize can support a complex diazotrophic microbiota enriched
for homologs of genes encoding nitrogenase subunits that harbor active
nitrogenase activity, and that nitrogen is transferred efficiently from the
nitrogen-fixing bacteria to the host plant tissues. Collectively, over
several years and locations, the 15N natural abundance and 15N enrichment
results of 2 selections of a Sierra Mixe indigenous maize landrace suggest
that its nitrogen nutrition when grown in its native environment is
partially fulfilled by fixation of atmospheric nitrogen. Nitrogen fixation
is a particularly difficult phenotype to evaluate, as all the techniques
available are prone to artifacts and can give different estimates for
nitrogen fixation [42
<http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.2006352#pbio.2006352.ref042>
,43
<http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.2006352#pbio.2006352.ref043>].
In this study, we used several classical techniques with Sierra Mixe maize
and have shown over multiple locations and years that Sierra Mixe maize can
fix nitrogen at rates (29%–82%) not previously reported, to our knowledge,
in maize.
The genetic basis of the trait or the source of the microbial inoculum,
which may be either environmental or seed-borne, are unresolved. The
observation that a teosinte species (*Z*. *mays* ssp. *mexicana*) also
exhibits a similar diazotrophic activity in aerial root–associated mucilage
suggests that this is an ancient trait that may have been introgressed and
amplified in the Sierra Mixe landrace. It will be important, in the future,
to determine the genetic basis of the trait, the identity of associated
microbial diazotrophs, and the mechanisms of microbial recruitment. This
research, together with other published research [5
<http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.2006352#pbio.2006352.ref005>
,18
<http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.2006352#pbio.2006352.ref018>
,23
<http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.2006352#pbio.2006352.ref023>],
suggests new avenues for research into potentially novel mechanisms of
biological N2 fixation in maize. This could have a significant impact on
maize crop productivity and nitrogen use efficiency, particularly in
regions of the world where agriculture is characterized by poor soil
nutrition.
Materials and methods
Plant material
Sierra Mixe maize seeds were obtained in Sierra Mixe region of Oaxaca,
Mexico, from an open pollinated population. *Z*. *mays* ssp. *mexicana*
(teosinte), LH123HT, Mo17, PHG39, LH82, PH207, and B73 seeds were obtained
from USDA National Plant Germplasm System (NPGS) (accessions Ames 8083,
PI601079, PI558532, PI600981, PI601170, PI601005, and PI550473,
respectively). Maize line Hickory King was obtained from Victory Seeds
(accession 3140041). Tornado-F21 and H-377 were obtained from Semillas
Ceres in Oaxaca. SM1 and SM2 are different selections from the Sierra Mixe
landrace. SM1 is a uniform population based on kernel size, shape, and
color, and SM2 is a heterogeneous population representing the landrace.
Biological materials were accessed and utilized under an Access and Benefit
Sharing Agreement between the Sierra Mixe community and BioN2, Inc., and
with permission from the Mexican government. An internationally recognized
certificate of compliance under the Nagoya Protocol
(ABSCH-IRCC-MX-207343-3) has been issued for such activities.
Bacterial strains and media
*A*. *brasilense* Sp7 and *B*. *unamae* MTI-641 were kindly provided by Dr.
G. Alexandre (University of Tennessee, Knoxville, USA) and Dr. A. Hirsch
(University of California, Los Angeles, USA), respectively. *H*.
*seropedicae* Z152 (ATCC 35894) was provided by the ATCC ( http://www.atcc.org/). Bacteria were grown in liquid BSE medium.
Sample collection
The rhizosphere and plant tissues that include stem, leaf, aerial roots,
underground roots, and mucilage of Sierra Mixe maize were sampled during
seasons 2010, 2011, and 2012 from Fields 1 and 2 in Sierra Mixe. For plant
endophyte analysis, tissues were surface sterilized by rinsing with mqH2O,
shaken gently in 70% ethanol for 5 minutes, placed into 1% hypochlorite
bleach, gently stirred for 10 minutes, rinsed 3 times in mqH20, and dried
in a laminar flow cabinet. Roots and stems were also dissected to remove
epidermal tissues before extraction. For seed endophyte analysis, embryo
and endosperm of Sierra Mixe, Hickory King, and B73 were withdrawn from the
seeds by hand using a razor blade in a laminar flow cabinet and were
collected in 1.5 ml sterile microcentrifuge tubes. For this study, the
rhizosphere was defined as a layer of soil covering the outer surface of
the root system that could be washed from roots in a buffer/detergent
solution. Roots were first separated and shaken to remove loosely adhering
soil. All soils and plant material samples were used immediately for DNA
extraction. Soil fertility analysis, which included physical parameters,
soil reaction, and salinity, was performed in AgroLab (Pachuma, Mexico).
Plant phenotyping
The number of nodes with aerial roots was monitored weekly (greenhouse) or
after 14 weeks (field). The total number of aerial roots was quantified
after 14 weeks. The disappearance of leaf wax and appearance of trichomes
were monitored weekly to determine the transition between juvenile stage
and adult stage in Sierra Mixe and Hickory King maize.
Greenhouse and field experiments, Madison, USA
For experiments in the greenhouse, seeds of Sierra Mixe of Fields 1 and 2
and Hickory King were surface sterilized and germinated as described
previously. After 1 week, the seedlings are transplanted in 40-liter pots
filled with a mix of sand and perlite (v:v) and grown in a high-ceiling
greenhouse at the Biotron facility (University of Wisconsin, Madison, USA).
Plants were watered twice a day for 2 minutes with half-strength of
Hoagland solution. For experiments in the field, 3 independent plots of 20
plants per genotype were planted, with 3 border rows (B73) between each
genotype. Sierra Mixe and Teosinte plants that were grown in Madison for
the ARA were planted in the same field at the same time. This experiment
was replicated in 3 different field plots.
Mucilage glycosyl composition
Glycosyl composition analysis was performed by combined gas
chromatography/mass spectrometry (GC/MS) of the per-*O*-trimethylsilyl
(TMS) derivatives of the monosaccharide methyl glycosides produced from the
sample by acidic methanolysis. Methyl glycosides were first prepared from
dry mucilage samples by methanolysis in 1 M HCl in methanol at 80 °C (18–22
hours), followed by re-*N*-acetylation with pyridine and acetic anhydride
in methanol (for detection of amino sugars). The samples were then
per-*O*-trimethylsilylated
by treatment with Tri-Sil (Pierce) at 80 °C (0.5 hours). GC/MS analysis of
the TMS methyl glycosides was performed on an HP 6890 GC interfaced to a
5975b MSD, using an All Tech EC-1 fused silica capillary column (30 m ×
0.25 mm ID). The analysis was performed at the Complex Carbohydrate Center
Research Center (CCRC) of University of Georgia, Athens, USA.
DNA extraction
DNA (80–150 ng μl−1) was extracted from 100 mg of the rhizosphere plant
tissues using a DNA isolation kit (Mo Bio Laboratories, Carlsbad, USA). The
PCR control for microbial DNA isolation was performed on 16S rRNA genes.
PCR was performed using eubacterial primers 27F
(5′-AGAGTTTGATCCTGGCTCAG-3′) and 1492R (5′-GGTTACCTTGTTACGACTT-3′), and the
product was approximately a 1,450 bp fragment. Amplification was carried
out with 1 μM of each primer in 3 mM MgCl2, 20 μM of each dNTP, 1.25 units
of *Taq* polymerase (Promega) in a total volume of 20 μl of 1X reaction
buffer (Promega, Madison, USA). PCR conditions included an initial
denaturation at 95 °C for 3 minutes followed by 35 cycles of denaturation
at 94 °C for 1 minute, annealing at 56 °C for 1 minute, and elongation at
72 °C for 1.5 minutes, with a final elongation at 72 °C for 7 minutes. DNA
was resolved using an agarose gel run at 100 V for 30 minutes for analysis
of total DNA and amplification of PCR products, respectively. Gels were
visualized by ethidium bromide staining under UV light in a gel
documentation system. The products obtained were purified with a NucleoSpin
Gel extraction kit (Clontech, Palo Alto, USA).
Illumina-based 16S rRNA gene sequencing
16S rRNA gene PCR and sequencing of rhizosphere and plant tissues were
carried out using the Caporaso protocol [45
<http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.2006352#pbio.2006352.ref045>].
We extended the Caporaso approach [46
<http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.2006352#pbio.2006352.ref046>]
to a dual barcode scheme for each sample and replaced the Golay barcodes
with a different set of Illumina-compatible barcodes that were designed to
balance base composition and tolerate up to 4 sequencing errors in barcode
sequences. The forward primer used was
AATGATACGGCGACCACCGAGATCTACAC[Barcode]TATGGTAATTGTGTGCCAGCMGCCGCGGTAA, and
the reverse primer was
CAAGCAGAAGACGGCATACGAGAT[Barcode]AGTCAGTCAGCCGGACTACHVGGGTWTCTAAT. The
barcodes used were designed to allow pooling of multiple samples within a
single MiSeq run. Ten cycles of PCR with barcoded primers were performed at
low annealing temperature (55 °C); samples were then pooled and cleaned
using a Qiagen column to remove the unincorporated primers. At this stage,
an additional 10 or 20 cycles of PCR were performed on the pool using the
Illumina paired-end flowcell primers with a higher annealing temperature
(65 °C). The resulting PCR product was subjected to QC with an Agilent
Bioanalyzer and estimated concentration using KAPA Biosystems qPCR kit. The
samples were diluted to the appropriate loading concentration for a MiSeq
run, spiked with 25% phiX control library, and sequenced using an Illumina
MiSeq instrument with the manufacturer’s standard 150 nucleotides
paired-end dual-index sequencing protocol and the custom sequencing primers.
The Illumina sequences were obtained from 2 MiSeq runs (2X150 bp
paired-end) and were demultiplexed using a custom script ( https://figshare.com/s/04997ae7f7d18b53174a). The 822,804 reads thus
obtained were trimmed to a Phred-equivalent of 20 and filtered for adaptor
contamination using BBDuk (of the BBTool packages https://sourceforge.net/projects/bbmap/). Because of the low quality of the
reverse mate pair, reads of the forward mate were used in the analysis. The
reads were then preprocessed and analyzed using DADA2 [47
<http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.2006352#pbio.2006352.ref047>].
Prescribed standard filtering parameters were used, such as PhiX
contamination check and removal of reads with more than 2 errors or
ambiguous bases or with an expected error greater than 2. Chimeras were
identified and removed using the removeBimeraDenovo function of DADA2. The
clean reads were then collapsed into sequence variants and classified using
RDP training set (version 14). The sequence variants that were classified
as chloroplast or mitochondria were removed from further analyses. From
samples with library size ranging from 84 to 20,597 reads (mean of 4,703),
995 unique sequence variants were identified.
For ARA with the mucilage, 2 ml of freshly collected mucilage from 1 or 2
aerial roots (Sierra Mixe) or several plants (teosinte) grown in the field
were introduced in 14.5 ml vials (Wheaton, Millville, USA) that were
tightly closed. For ARA with added bacteria, *A*. *brasilense*, *H*.
*seropedicae*, and *B*. *unamae* were grown in BSE medium for 48 hours.
Then, bacteria were collected by centrifugation (5 minutes, 4,000 × *g*)
and suspended in the Fahraeus medium. The 14.5 ml vials, each containing 5
ml of mucilage previously stored for several months at –20 °C to reduce
endogenous nitrogen-fixing bacteria, were inoculated with the bacterial
suspension at a final OD600nm = 0.01. Control tubes were prepared either
without bacteria or with 5 ml of Fahraeus medium instead of mucilage. Then,
850 μl of acetylene (Airgas) was injected into each vial. OD600nm was
measured for each tube after 72 hours. For both conditions, controls
without acetylene were performed in parallel. For ARA with aerial roots, 1
aerial root without mucilage was introduced in each 14.5 ml vial (10
replicates). One ml of acetylene (Airgas) was injected into each vial. For
ARA with seedlings, Sierra Mixe seedlings were inoculated with *A*.
*brasilense* and grown for 3 weeks. Plants were then transferred to 500 ml
jars, and 50 ml of acetylene was injected in each jar. For ARA with
underground roots, pieces of roots (about 10 cm long) were collected from
plants grown in pots and introduced into 14.5 ml vials (3 replicates).
Ethylene quantification was made by injecting 1 ml of the air phase,
sampled after 72 hours, on a gas chromatography (GC-2010 Shimadzu) equipped
with a Rt-Alumina BOND/KCL column (Restek).
Mucilage 15N2 assimilation
The enrichment of mucilage in 15N atom was achieved by removing 4 ml of
headspace gas and replacing it with 4 ml of either 15N2 (Sigma-Aldrich) or
14N2 nitrogen gas directly into a vial containing 1.0 mL of mucilage.
Mucilage was collected from Sierra Mixe maize plants grown in Sierra Mixe
and stored at 4 °C for up to 2 weeks between sampling and the determination
of 15N2 assimilation. The mucilage samples were incubated at 37 °C for 0
and 70 hours in the presence of 15N2. 15N2 assimilation was stopped by
freezing the mucilage samples at −20 °C. The samples were then freeze-dried
and weighed. The 15N2 analysis in the mucilage samples was performed at the
UC Davis Stable Isotope Facility (Davis, USA) and the UW-Madison Soil
Science Facility (Madison, USA). Statistical analysis was performed using
SYSTAT version 10 (Chicago, USA).
Measurement of free-oxygen concentration
For measurement in collected mucilage, 2 ml of mucilage was introduced in a
15 ml tube. The probe (robust oxygen mini probe, Pyroscience) was
introduced 8 mm deep in the mucilage and oxygen measurements performed
until stabilization of the signal was observed. Control corresponds to
free-oxygen concentration in the liquid Fahraeus medium. One-point
calibration was made in aerated water, as advised by the manufacturer.
15N2 gas–enrichment experiments
Aerial roots were collected from Sierra Mixe maize grown at the Biotron
greenhouse facility (University of Wisconsin, Madison, USA). Mucilage was
generated from each of these aerial roots by incubating them in 5 ml of
water at room temperature for 48 hours. Mucilage, along with the aerial
roots, was inoculated with *A*. *brasilense* Sp7. Then, 10%–15% (v/v) of 15N
2 gas was pumped into the vials, and the samples were incubated at 30 °C
for 48 hours. After incubation of mucilage alone, or aerial roots alone,
pheophytin extracted from these aerial roots was subjected to IRMS
analysis. To obtain pheophytin, chlorophyll was extracted from aerial roots
and converted to pheophytin by acid treatment, following as described [59
<http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.2006352#pbio.2006352.ref059>].
14N2-treated mucilage and aerial roots were used as negative controls.
15N natural abundance
The proportion (%) of nitrogen derived from biological nitrogen fixation
(%Ndfa) was estimated from the 15N natural abundance (expressed in delta
units, ‰) of the Sierra Mixe maize (δ15Nfixing plant) and that of the
reference plant species (δ15Nref). In each of the 2011 and 2012 field
seasons in Sierra Mixe, 90–114 individual maize samples (depending on the
year) and 270 reference plant samples, representing 8–10 species (depending
on the year) of non-nitrogen-fixing plants, were analyzed. For the single
time point in 2010, 12 individual maize samples and 33 reference plant
samples, representing 8 species of non-nitrogen-fixing plants, were
analyzed. The reference plant species in the field were identified using
universal 18S PCR analysis from DNA sampled using FTA Plant Saver card (GE
Life Sciences, Pittsburg, USA) simultaneously with tissue samples collected
for 15N analysis. PCRs were performed using Sigma’s Extract-N-Amp Plant PCR
Kit according to the manufacturer for sequencing and BLASTN comparison. For
15N natural abundance, the third-youngest leaf of Sierra Mixe maize or
reference plants was collected from Field 3 and 4 from the second to the
sixth month postplanting and analyzed for N-isotope composition. Total
organic nitrogen was determined by Kjeldahl digestion followed by steam
distillation. Analysis for natural 15N abundance was carried out as
described by Bremer and van Kessel [38
<http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.2006352#pbio.2006352.ref038>].
The 15N analysis was performed at UCD Stable Isotope Facility ( http://stableisotopefacility.ucdavis.edu/13cand15n.html). The percentage of
nitrogen derived from nitrogen fixation (%Ndfa) was calculated as
follows: where
“δ15N” is stable nitrogen isotopes, “ref” is the value from non-N-fixing
reference plants, “fixing plant” is Sierra Mixe maize, and “B” is the 15N
abundance in the air, assumed to be 0.0‰.
15N-enrichment field experiments
In 2016, 3 locations were chosen: Field 3, land that had not been planted
to crops for over 10 years; Field 4, land that had maize for 1 year; and
Field 5, land with continuous maize. A randomized complete block design
trial was established at each site, with 5 replicates with 4 varieties.
Each plot consisted of 6 matas surrounded by a common border of SM2 on all
sides and a double border on outside rows. A mata is the traditional
planting design in the Sierra Mixe region, similar to a hill plot in which
multiple plants are seeded together. Each mata was planted with 5 seeds and
thinned to 3 seeds for 18 plants per plot. Matas were planted in a grid 80
× 100 cm from each other. 15N was applied at a dose of 0.36 grams per plot
in a liquid solution, reaching the desired enrichment of 1% for all 3
fields, with 50 ml added per mata.
In 2017, the same 3 locations were planted with same design and entries,
except Field 3 consisted only of H377 and SM2 entries. In all 3 fields, 15N
was applied at a dose of 0.95 grams per plot, reaching the desired
enrichment of greater than 1%. A solution was spread evenly over each plot
using a garden watering can, such that the whole experimental area received
an equal amount of enriched 15N. Plants were covered with plastic bags at
the time of application (V5) to ensure that 15N was not directly applied to
the leaves and that the 15N was uniformly available to all plants.
Soil samples were taken from a 0–60 cm depth in each plot, blended, and
sent for analysis at UC Davis Soil lab. Means were calculated across
locations. In 2016, at V9 and V12, 1 mata (3 plants) was sampled; and at
Tassel, 4 matas (12 plants) were sampled. In 2017, a single sampling of 6
matas (18 plants) was sampled at Tassel. For each sampling, plants were dug
out to include all roots. Because of the high rainfall (2,100 mm
concentrated in the growing season from June to October), roots were
shallow for both reference and test varieties. Each plant was photographed,
and data were recorded for the number of plants, plant height, the total
fresh weight of shoots, and roots and stem diameter. Whole plants were
chopped, ground, subsampled, and dried in an oven to record total dry
weight for shoots and roots. Well-blended subsamples were taken and shipped
to Davis to measure Total N and 15N. Total nitrogen and Atom% 15N were
determined for each plot at each time point for the shoot. Total organic
nitrogen was determined by Kjeldahl digestion followed by steam
distillation. The 15N analysis was performed at UCD Stable Isotope Facility
(http://stableisotopefacility.ucdavis.edu/13cand15n.html).
Nitrogen difference (%*NDiff*) method was calculated as
The average of Tornado F21 and H377 was used as the reference to calculate
%Ndfa and %Ndiff. For Ndiff calculations, the area harvested was adjusted
based on matas harvested at each time point.
Field data analyses
Data were analyzed using the R lme4 package. Data were checked for outliers
and subjected to ANOVA and mean separation using Least Significant
Difference (*p* = 0.05) for each location. LSDs were calculated only when
ANOVA F-tests were significant at *p* = 0.05. For %Ndfa and %Ndiff,
single-degree-of-freedom contrasts were calculated to compare test
varieties to the mean of the reference varieties (*P* > 0.05). Pearson
correlation coefficients were calculated between %Ndfa and %Ndiff.
Supporting information
pbio.2006352.s001.eps
sorry, we can't preview this file
fig*share*
<https://figshare.com/articles/Nitrogen_fixation_in_a_landrace_of_maize_is_supported_by_a_mucilage-associated_diazotrophic_microbiota/6941555>
*1 / 12*
Box plot indicating the alpha diversity as calculated by Phyloseq using (A)
Simpson index and (B) Shannon index.
The heat map depicts the abundances of (A) 1,000 most abundant SVs and (B)
20 most abundant SVs in the dataset that were transformed using
variance-stabilizing transformation in DESeq2. The libraries are Mucilage
(Blue; OLMC00, OLMD00, OLMV00, OLMX00), Aerial Root (Gray; OLAR00, OLAR02,
OLAR04, OLAR05), Aerial Root with Mucilage (Pink; OLAR01, OLAR03), Stem
(Green; OLST00, OLST01, OLST02, OLST03), Underground Root (Brown; OLUR01,
OLUR02, OLUR03), and Rhizosphere (Magenta; OXRZ11, OXRZ12, OXRZ13, OXRZ21,
OXRZ22, OXRZ23, OXRZ32, OXRZ31, OXRZ33).
(A) Number of nodes with aerial roots and (B) number of aerial roots
observed on teosinte, Sierra Mixe maize, and Hickory King after 14 weeks.
Bar = standard error of the mean. Different letters indicate statistically
supported groups according to the Kruskal-Wallis test.
(A) Oxygen measured at 3 depths in Fahraeus medium with (black bars) or
without (gray bars) 0.2% agar. (B) Effect of the different sugars present
in the mucilage on the ability of *H*. *seropedicae*, (C) *A*. *brasilense*,
and (D) *B*. *unamae* to reduce acetylene.
>From each location, 6 leaf samples were randomly sampled from Sierra Mixe
maize plants and 6 leaf samples from each of 2 reference plants. The third
emergent leaf of each maize plant was sampled. Reference plants were
selected from the most abundant weed species within each sample location,
and from a plant family (Asteraceae and Ranunculaceae) that is neither
actinorhizal nor leguminous nor has members known to associate with
diazotrophic bacteria. δ15N was determined for each plant sampled, and
%Ndfa was calculated for Sierra Mixe maize according to the equation 2 in [
19
<http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.2006352#pbio.2006352.ref019>].
Values are given as mean and s.e. Different letters indicate statistically
supported groups (one-way ANOVA, *P* < 0.05). %Ndfa, percent of nitrogen
derived from the atmosphere.
This paper is dedicated to the life and memory of Cristobal Heitmann.
Cristobal’s energy and enthusiasm was a major catalyst in completing this
research. He died tragically while the manuscript was under review. We
thank Vicente Vasquez and Maria del Refugio Vasquez for assistance in
developing the program in Mexico; Carmen Ortega and Saulon Zamora for
sample collection and field trial assistance; Shawn Kaeppler, Natalia de
Leon, and Jillian Foerster for field assistance; Nguyet Dao for
microbial 15N-fixation
assays; Harry Read for processing 15N2-enrichment samples; John Zhang for
assistance with library construction; and Armando Garcia-Llanos for
assistance in sample processing. We thank the Comisiriado of the Sierra
Mixe, Mexico, for their support and access to community genetic resources.
References
1. 1. Bulgarelli D, Rott M, Schlaeppi K, van Themaat EVL, Ahmadinejad N,
Assenza F, et al. Revealing structure and assembly cues for Arabidopsis
root-inhabiting bacterial microbiota. Nature. 2012;488(7409):91–5.
pmid:22859207
- View Article <https://doi.org/10.1038/nature11336>
- PubMed/NCBI <http://www.ncbi.nlm.nih.gov/pubmed/22859207>
- Google Scholar