Metabolic Adaptation, a Specialized Leaf Organ Structure and Vascular Responses to Diurnal N2 Fixation by Nostoc azollae Sustain the Astonishing Productivity of Azolla Ferns without Nitrogen Fertilizer

Biosynthesis and differential spatial distribution of the 3-deoxyanthocyanidins apigenidin and luteolinidin at the interface of a plant-cyanobacteria symbiosis exposed to cold

Aquatic ferns of the genus Azolla (Azolla) form highly productive symbioses with filamentous cyanobacteria fixing N2 in their leaf cavities, Nostoc azollae. Stressed symbioses characteristically turn red due to 3-deoxyanthocyanidin (DA) accumulation, rare in angiosperms and of unknown function. To understand DA accumulation upon cold acclimation and recovery, we integrated laser-desorption-ionization mass-spectrometry-imaging (LDI-MSI), a new Azolla filiculoides genome-assembly and annotation, and dual RNA-sequencing into phenotypic analyses of the symbioses. Azolla sp. Anzali recovered even when cold-induced DA-accumulation was inhibited by abscisic acid. Cyanobacterial filaments generally disappeared upon cold acclimation and Nostoc azollae transcript profiles were unlike those of resting stages formed in cold-resistant sporocarps, yet filaments re-appeared in leaf cavities of newly formed green fronds upon cold-recovery. The high transcript accumulation upon cold acclimation of AfDFR1 encoding a flavanone 4-reductase active in vitro suggested that the enzyme of the first step in the DA-pathway may regulate accumulation of DAs in different tissues. However, LDI-MSI highlighted the necessity to describe metabolite accumulation beyond class assignments as individual DA and caffeoylquinic acid metabolites accumulated differentially. For example, luteolinidin accumulated in epithelial cells, including those lining the leaf cavity, supporting a role for the former in the symbiotic interaction during cold acclimation.

The crane fly glycosylated triketide δ-lactone cornicinine elicits akinete differentiation of the cyanobiont in aquatic Azolla fern symbioses

The restriction of plant-symbiont dinitrogen fixation by an insect semiochemical had not been previously described. Here we report on a glycosylated triketide δ-lactone from Nephrotoma cornicina crane flies, cornicinine, that causes chlorosis in the floating-fern symbioses from the genus Azolla. Only the glycosylated trans-A form of chemically synthesized cornicinine was active: 500 nM cornicinine in the growth medium turned all cyanobacterial filaments from Nostoc azollae inside the host leaf-cavities into akinetes typically secreting CTB-bacteriocins. Cornicinine further inhibited akinete germination in Azolla sporelings, precluding re-establishment of the symbiosis during sexual reproduction. It did not impact development of the plant Arabidopsis thaliana or several free-living cyanobacteria from the genera Anabaena or Nostoc but affected the fern host without cyanobiont. Fern-host mRNA sequencing from isolated leaf cavities confirmed high NH4-assimilation and proanthocyanidin biosynthesis in this trichome-rich tissue. After cornicinine treatment, it revealed activation of Cullin-RING ubiquitin-ligase-pathways, known to mediate metabolite signaling and plant elicitation consistent with the chlorosis phenotype, and increased JA-oxidase, sulfate transport and exosome formation. The work begins to uncover molecular mechanisms of cyanobiont differentiation in a seed-free plant symbiosis important for wetland ecology or circular crop-production today, that once caused massive CO2 draw-down during the Eocene geological past.

Comparative impact of replacing fish meal with azolla on growth, water quality, and physiology of red tilapia fingerlings at varying salinities

A 210-day experiment to assess the efficacy of substituting azolla plant powder at levels of 0, 20, 40, and 60% for fish meal on red tilapia fingerlings (RTF, initial weight of 18.23 ± 0.12 g) performance under salinity levels of 5, 18, and 28ppt. Among the various conditions, RTF-fed 20% azolla at 28 and 5ppt salinity showcased the highest specific growth rate (SGR), whereas the lowest SGR was observed in fish-fed 60% azolla at 5ppt salinity. Upon azolla incorporation, noteworthy elevations in phytoplankton, zooplankton, dissolved oxygen (DO), pH, NH3, and NO3 were noted and conversely, azolla introduction led to decreased NH4 and NO2 concentrations in all salinity levels. Further, a significant (p < 0.05) interaction between azolla levels and water salinity (S×A) significantly impacted the hematological parameters of RTF. The highest levels of superoxide dismutase (SOD), catalase (CAT), and total protein (TP) were found in RTF-fed 20% azolla at 28ppt salinity, while the lowest CAT and TP levels occurred in RTF-fed 60% azolla at 5ppt salinity. The highest aspartate aminotransferase (AST) and alanine aminotransferase (ALT) levels were recorded in the RTF group fed 60% azolla at 5ppt salinity, with the lowest values seen in the group given 20% azolla at 28ppt salinity. RTF fed a 20% azolla diet at 18ppt salinity exhibited the highest lysozyme value, in contrast to the lowest value observed in the RTF group fed the control diet at 18ppt salinity. In conclusion, this study recommends the utilization of azolla at inclusion levels ranging from 20 to 40%, as it has the potential to notably enhance the immune system and elevate the survival rate of RTF.

Harnessing biological nitrogen fixation in plant leaves

The importance of biological nitrogen fixation (BNF) in securing food production for the growing world population with minimal environmental cost has been increasingly acknowledged. Leaf surfaces are one of the biggest microbial habitats on Earth, harboring diverse free-living N2-fixers. These microbes inhabit the epiphytic and endophytic phyllosphere and contribute significantly to plant N supply and growth. Here, we summarize the contribution of phyllosphere-BNF to global N cycling, evaluate the diversity of leaf-associated N2-fixers across plant hosts and ecosystems, illustrate the ecological adaptation of N2-fixers to the phyllosphere, and identify the environmental factors driving BNF. Finally, we discuss potential BNF engineering strategies to improve the nitrogen uptake in plant leaves and thus sustainable food production.

Impact of High Light Intensity and Low Temperature on the Growth and Phenylpropanoid Profile of Azolla filiculoides

Exposure to high light intensity (HL) and cold treatment (CT) induces reddish pigmentation in Azolla filiculoides, an aquatic fern. Nevertheless, how these conditions, alone or in combination, influence Azolla growth and pigment synthesis remains to be fully elucidated. Likewise, the regulatory network underpinning the accumulation of flavonoids in ferns is still unclear. Here, we grew A. filiculoides under HL and/or CT conditions for 20 days and evaluated the biomass doubling time, relative growth rate, photosynthetic and non-photosynthetic pigment contents, and photosynthetic efficiency by chlorophyll fluorescence measurements. Furthermore, from the A. filiculoides genome, we mined the homologs of MYB, bHLH, and WDR genes, which form the MBW flavonoid regulatory complex in higher plants, to investigate their expression by qRT-PCR. We report that A. filiculoides optimizes photosynthesis at lower light intensities, regardless of the temperature. In addition, we show that CT does not severely hamper Azolla growth, although it causes the onset of photoinhibition. Coupling CT with HL stimulates the accumulation of flavonoids, which likely prevents irreversible photoinhibition-induced damage. Although our data do not support the formation of MBW complexes, we identified candidate MYB and bHLH regulators of flavonoids. Overall, the present findings are of fundamental and pragmatic relevance to Azolla's biology.

Evolutionary genomic insights into cyanobacterial symbioses in plants

Photosynthesis, the ability to fix atmospheric carbon dioxide, was acquired by eukaryotes through symbiosis: the plastids of plants and algae resulted from a cyanobacterial symbiosis that commenced more than 1.5 billion years ago and has chartered a unique evolutionary path. This resulted in the evolutionary origin of plants and algae. Some extant land plants have recruited additional biochemical aid from symbiotic cyanobacteria; these plants associate with filamentous cyanobacteria that fix atmospheric nitrogen. Examples of such interactions can be found in select species from across all major lineages of land plants. The recent rise in genomic and transcriptomic data has provided new insights into the molecular foundation of these interactions. Furthermore, the hornwort Anthoceros has emerged as a model system for the molecular biology of cyanobacteria-plant interactions. Here, we review these developments driven by high-throughput data and pinpoint their power to yield general patterns across these diverse symbioses.

Moroccan Lagoon Microbiomes

Lagoons are fragile marine ecosystems that are considerably affected by anthropogenic pollutants. We performed a spatiotemporal characterization of the microbiome of two Moroccan lagoons, Marchica and Oualidia, both classified as Ramsar sites, the former on the Mediterranean coast and the latter on the Atlantic coast. We investigated their microbial diversity and abundance using 16S rRNA amplicon- and shotgun-based metagenomics approaches during the summers of 2014 and 2015. The bacterial microbiome was composed primarily of Proteobacteria (25–53%, 29–29%), Cyanobacteria (34–12%, 11–0.53%), Bacteroidetes (24–16%, 23–43%), Actinobacteria (7–11%, 13–7%), and Verrucomicrobia (4–1%, 15–14%) in Marchica and Oualidia in 2014 and 2015, respectively. Interestingly, 48 strains were newly reported in lagoon ecosystems, while eight unknown viruses were detected in Mediterranean Marchica only. Statistical analysis showed higher microbial diversity in the Atlantic lagoon than in the Mediterranean lagoon and a robust relationship between alpha diversity and geographic sampling locations. This first-ever metagenomics study on Moroccan aquatic ecosystems enriched the national catalog of marine microorganisms. They will be investigated as candidates for bioindication properties, biomonitoring potential, biotechnology valorization, biodiversity protection, and lagoon health assessment.

Nitrogen fixation: a poorly understood process along the freshwater-marine continuum

Although N₂ fixation is a major component of the global N cycle and has been extensively studied in open-ocean and terrestrial ecosystems, rates and ecological dynamics remain virtually unknown for the inland and coastal aquatic ecosystems (lakes, wetlands, rivers, streams, estuaries) that connect terrestrial and marine biomes. This is due to the diversity of these habitats, as well as the traditional paradigm that N₂ fixation rates were low to nonexistent, and therefore not important, in these ecosystems. We identify three major research themes to advance understanding of aquatic N₂ fixation: 1) the biological diversity of diazotrophs and variability of N₂ fixation rates, 2) the ecological stoichiometry of N₂ fixation, and 3) the upscaling of N₂ fixation rates from genes to ecosystems. Coordinating research across these areas will advance limnology and oceanography by fully integrating N₂ fixation into ecological dynamics of aquatic ecosystems from local to global scales.

Perspectives on the utilization of Azolla-Anabaena system as feed supplement

The nitrogen fixing aquatic pteridophyte Azolla is one of the fastest growing nitrogen-fixing plants and it is used as a potential source for high rate biomass production. Azolla has the ability to fix atmospheric nitrogen at cheaper and faster rates due to the presence of a symbiotic cyanobacterium Anabaena azollae. Therefore, the ability to fix atmospheric nitrogen is important from an agricultural perspective. However, Azolla is gaining popularity as feed supplement for cattle, poultry and fish. Further, the ease of cultivation and favourable nutrient composition make Azolla an important feed supplement. This review focuses on the perspectives of Azolla as feed supplement.

Light and Temperature Shape the Phenylpropanoid Profile of Azolla filiculoides Fronds

Azolla is a genus of floating freshwater ferns. By their high growth and N₂ fixation rates, Azolla species have been exploited for centuries by populations of South-east Asia as biofertilizers in rice paddies. The use of Azolla species as a sustainable plant material for diverse applications, such as feeding, biofuel production, and bioremediation, has encountered a growing interest over the last few years. However, high levels of feed deterrent flavonoids in their fronds have discouraged the use of these ferns as a sustainable protein source for animal consumption. Additionally, information on how and to what extent environmental determinants affect the accumulation of secondary metabolites in these organisms remains poorly understood. Moving from these considerations, here, we investigated by an untargeted metabolomics approach the profiles of phenylpropanoid compounds in the fronds of Azolla filiculoides sampled under control and pigment-inducing stress conditions. In parallel, we assayed the expression of essential structural genes of the phenylpropanoid pathway by quantitative RT-PCR. This study provides novel information concerning A. filiculoides phenylpropanoid compounds and their temporal profiling in response to environmental stimuli. In particular, we show that besides the already known 3-deoxyanthocyanidins, anthocyanidins, and proanthocyanidins, this fern can accumulate additional secondary metabolites of outstanding importance, such as chemoattractants, defense compounds, and reactive oxygen species (ROS) scavengers, and crucial as dietary components for humans, such as dihydrochalcones, stilbenes, isoflavones, and phlobaphenes. The findings of this study open an opportunity for future research studies to unveil the interplay between genetic and environmental determinants underlying the elicitation of the secondary metabolites in ferns and exploit these organisms as sustainable sources of beneficial metabolites for human health.

Azolla ferns testify: seed plants and ferns share a common ancestor for leucoanthocyanidin reductase enzymes

Questions about in vivo substrates for proanthocyanidin (PA) biosynthesis and condensation have not been resolved and wide gaps in the understanding of transport and biogenesis in 'tannosomes' persist. Here we examined the evolution of PA biosynthesis in ferns not previously reported, asking what PAs are synthesised and how. Chemical and gene-expression analyses were combined to characterise PA biosynthesis, leveraging genome annotation from the floating fern Azolla filiculoides. In vitro assay and phylogenomics of PIP-dehydrogenases served to infer the evolution of leucoanthocyanidin reductase (LAR). Sporophyte-synthesised (epi)catechin polymers, averaging only seven subunits, accumulated to 5.3% in A. filiculoides, and 8% in A. pinnata biomass dry weight. Consistently, a LAR active in vitro was highly expressed in A. filiculoides. LAR, and paralogous fern WLAR-enzymes with differing substrate binding sites, represent an evolutionary innovation of the common ancestor of fern and seed plants. The specific ecological niche of Azolla ferns, a floating plant-microbe mat massively fixing CO₂ and N₂ , shaped their metabolism in which PA biosynthesis predominates and employs novel fern LAR enzymes. Characterisation of in vivo substrates of these LAR, will help to shed light on the recently assigned and surprising dual catalysis of LAR from seed plants.

Far-Red Light-Induced Azolla filiculoides Symbiosis Sexual Reproduction: Responsive Transcripts of Symbiont Nostoc azollae Encode Transporters Whilst Those of the Fern Relate to the Angiosperm Floral Transition

Water ferns of the genus Azolla and the filamentous cyanobacteria Nostoc azollae constitute a model symbiosis that enabled the colonization of the water surface with traits highly desirable for the development of more sustainable crops: their floating mats capture CO₂ and fix N₂ at high rates using light energy. Their mode of sexual reproduction is heterosporous. The regulation of the transition from the vegetative phase to the spore forming phase in ferns is largely unknown, yet a prerequisite for Azolla domestication, and of particular interest as ferns represent the sister lineage of seed plants. Sporocarps induced with far red light could be crossed so as to verify species attribution of strains from the Netherlands but not of the strain from the Anzali lagoon in Iran; the latter strain was assigned to a novel species cluster from South America. Red-dominated light suppresses the formation of dissemination stages in both gametophyte- and sporophyte-dominated lineages of plants, the response likely is a convergent ecological strategy to open fields. FR-responsive transcripts included those from MIKC^C homologues of CMADS1 and miR319-controlled GAMYB transcription factors in the fern, transporters in N. azollae, and ycf2 in chloroplasts. Loci of conserved microRNA (miRNA) in the fern lineage included miR172, yet FR only induced miR529 and miR535, and reduced miR319 and miR159. Phylogenomic analyses of MIKC^C TFs suggested that the control of flowering and flower organ specification may have originated from the diploid to haploid phase transition in the homosporous common ancestor of ferns and seed plants.

TMEA: A Thermodynamically Motivated Framework for Functional Characterization of Biological Responses to System Acclimation

The objective of gene set enrichment analysis (GSEA) in modern biological studies is to identify functional profiles in huge sets of biomolecules generated by high-throughput measurements of genes, transcripts, metabolites, and proteins. GSEA is based on a two-stage process using classical statistical analysis to score the input data and subsequent testing for overrepresentation of the enrichment score within a given functional coherent set. However, enrichment scores computed by different methods are merely statistically motivated and often elusive to direct biological interpretation. Here, we propose a novel approach, called Thermodynamically Motivated Enrichment Analysis (TMEA), to account for the energy investment in biological relevant processes. Therefore, TMEA is based on surprisal analysis, which offers a thermodynamic-free energy-based representation of the biological steady state and of the biological change. The contribution of each biomolecule underlying the changes in free energy is used in a Monte Carlo resampling procedure resulting in a functional characterization directly coupled to the thermodynamic characterization of biological responses to system perturbations. To illustrate the utility of our method on real experimental data, we benchmark our approach on plant acclimation to high light and compare the performance of TMEA with the most frequently used method for GSEA.

Understanding transport processes in lichen, Azolla-cyanobacteria, ectomycorrhiza, endomycorrhiza, and rhizobia-legume symbiotic interactions

Intimate interactions between photosynthetic and non-photosynthetic organisms require the orchestrated transfer of ions and metabolites between species. We review recent progress in identifying and characterizing the transport proteins involved in five mutualistic symbiotic interactions: lichens, Azolla–cyanobacteria, ectomycorrhiza, endomycorrhiza, and rhizobia–legumes. This review focuses on transporters for nitrogen and carbon and other solutes exchanged in the interactions. Their predicted functions are evaluated on the basis of their transport mechanism and prevailing transmembrane gradients of H ⁺ and transported substrates. The symbiotic interactions are presented in the assumed order from oldest to most recently evolved.

Jasmonic and salicylic acid response in the fern Azolla filiculoides and its cyanobiont

Plants sense and respond to microbes utilizing a multilayered signalling cascade. In seed plants, the phytohormones jasmonic and salicylic acid (JA and SA) are key denominators of how plants respond to certain microbes. Their interplay is especially well-known for tipping the scales in plants' strategies of dealing with phytopathogens. In non-angiosperm lineages, the interplay is less well understood, but current data indicate that it is intertwined to a lesser extent and the canonical JA/SA antagonism appears to be absent. Here, we used the water fern Azolla filiculoides to gain insights into the fern's JA/SA signalling and the molecular communication with its unique nitrogen fixing cyanobiont Nostoc azollae, which the fern inherits both during sexual and vegetative reproduction. By mining large-scale sequencing data, we demonstrate that Azolla has most of the genetic repertoire to produce and sense JA and SA. Using qRT-PCR on the identified biosynthesis and signalling marker genes, we show that Azolla is responsive to exogenously applied SA. Furthermore, exogenous SA application influenced the abundance and gene expression of Azolla's cyanobiont. Our data provide a framework for JA/SA signalling in ferns and suggest that SA might be involved in Azolla's communication with its vertically inherited cyanobiont.

Growing Azolla to produce sustainable protein feed: the effect of differing species and CO2 concentrations on biomass productivity and chemical composition

BACKGROUND

Since available arable land is limited and nitrogen fertilizers pollute the environment, cropping systems ought to be developed that do not rely on them. Here we investigate the rapidly growing, N₂ -fixing Azolla/Nostoc symbiosis for its potential productivity and chemical composition to determine its potential as protein feed.

RESULTS

In a small production system, cultures of Azolla pinnata and Azolla filiculoides were continuously harvested for over 100 days, yielding an average productivity of 90.0-97.2 kg dry weight (DW) ha^-1  d^-1 . Under ambient CO₂ levels, N₂ fixation by the fern's cyanobacterial symbionts accounted for all nitrogen in the biomass. Proteins made up 176-208 g kg^-1 DW (4.9 × total nitrogen), depending on species and CO₂ treatment, and contained more essential amino acids than protein from soybean. Elevated atmospheric CO₂ concentrations (800 ppm) significantly boosted biomass production by 36-47%, without decreasing protein content. Choice of species and CO₂ concentrations further affected the biomass content of lipids (79-100 g kg^-1 DW) and (poly)phenols (21-69 g kg^-1 DW).

CONCLUSIONS

By continuous harvesting, high protein yields can be obtained from Azolla cultures, without the need for nitrogen fertilization. High levels of (poly)phenols likely contribute to limitations in the inclusion rate of Azolla in animal diets and need further investigation. © 2018 The Authors. Journal of the Science of Food and Agriculture published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Fern genomes elucidate land plant evolution and cyanobacterial symbioses

Ferns are the closest sister group to all seed plants, yet little is known about their genomes other than that they are generally colossal. Here, we report on the genomes of Azolla filiculoides and Salvinia cucullata (Salviniales) and present evidence for episodic whole-genome duplication in ferns-one at the base of 'core leptosporangiates' and one specific to Azolla. One fern-specific gene that we identified, recently shown to confer high insect resistance, seems to have been derived from bacteria through horizontal gene transfer. Azolla coexists in a unique symbiosis with N2-fixing cyanobacteria, and we demonstrate a clear pattern of cospeciation between the two partners. Furthermore, the Azolla genome lacks genes that are common to arbuscular mycorrhizal and root nodule symbioses, and we identify several putative transporter genes specific to Azolla-cyanobacterial symbiosis. These genomic resources will help in exploring the biotechnological potential of Azolla and address fundamental questions in the evolution of plant life.

Azolla along a phosphorus gradient: biphasic growth response linked to diazotroph traits and phosphorus-induced iron chlorosis

Azolla spp., a water fern often used for phytoremediation, is a strong phosphorus (P) accumulator due to its high growth rate and N₂ fixing symbionts (diazotrophs). It is known that plant growth is stimulated by P, but the nature of the interactive response of both symbionts along a P gradient, and related changes in growth-limiting factors, are unclear. We determined growth, and N and P sequestration rates of Azolla filiculoides in N-free water at different P concentrations. The growth response appeared to be biphasic and highest at levels ≥10 P µmol l⁻¹. Diazotrophic N sequestration increased upon P addition, and rates were three times higher at high P than at low P. At 10 µmol P l⁻¹, N sequestration rates reached its maximum and A. filiculoides growth became saturated. Due to luxury consumption, P sequestration rates increased until 50 µmol P l⁻¹. At higher P concentrations (≥50 µmol l⁻¹), however, chlorosis occurred that seems to be caused by iron- (Fe-), and not by N-deficiency. We demonstrate that traits of the complete symbiosis in relation to P and Fe availability determine plant performance, stressing the role of nutrient stoichiometry. The results are discussed regarding Azolla’s potential use in a bio-based economy.

Is there foul play in the leaf pocket? The metagenome of floating fern Azolla reveals endophytes that do not fix N2 but may denitrify

Dinitrogen fixation by Nostoc azollae residing in specialized leaf pockets supports prolific growth of the floating fern Azolla filiculoides. To evaluate contributions by further microorganisms, the A. filiculoides microbiome and nitrogen metabolism in bacteria persistently associated with Azolla ferns were characterized. A metagenomic approach was taken complemented by detection of N₂ O released and nitrogen isotope determinations of fern biomass. Ribosomal RNA genes in sequenced DNA of natural ferns, their enriched leaf pockets and water filtrate from the surrounding ditch established that bacteria of A. filiculoides differed entirely from surrounding water and revealed species of the order Rhizobiales. Analyses of seven cultivated Azolla species confirmed persistent association with Rhizobiales. Two distinct nearly full-length Rhizobiales genomes were identified in leaf-pocket-enriched samples from ditch grown A. filiculoides. Their annotation revealed genes for denitrification but not N₂ -fixation. ¹⁵ N₂ incorporation was active in ferns with N. azollae but not in ferns without. N₂ O was not detectably released from surface-sterilized ferns with the Rhizobiales. N₂ -fixing N. azollae, we conclude, dominated the microbiome of Azolla ferns. The persistent but less abundant heterotrophic Rhizobiales bacteria possibly contributed to lowering O₂ levels in leaf pockets but did not release detectable amounts of the strong greenhouse gas N₂ O.

Is there foul play in the leaf pocket? The metagenome of floating fern Azolla reveals endophytes that do not fix N2 but may denitrify

Dinitrogen fixation by Nostoc azollae residing in specialized leaf pockets supports prolific growth of the floating fern Azolla filiculoides. To evaluate contributions by further microorganisms, the A. filiculoides microbiome and nitrogen metabolism in bacteria persistently associated with Azolla ferns were characterized. A metagenomic approach was taken complemented by detection of N₂ O released and nitrogen isotope determinations of fern biomass. Ribosomal RNA genes in sequenced DNA of natural ferns, their enriched leaf pockets and water filtrate from the surrounding ditch established that bacteria of A. filiculoides differed entirely from surrounding water and revealed species of the order Rhizobiales. Analyses of seven cultivated Azolla species confirmed persistent association with Rhizobiales. Two distinct nearly full-length Rhizobiales genomes were identified in leaf-pocket-enriched samples from ditch grown A. filiculoides. Their annotation revealed genes for denitrification but not N₂ -fixation. ¹⁵ N₂ incorporation was active in ferns with N. azollae but not in ferns without. N₂ O was not detectably released from surface-sterilized ferns with the Rhizobiales. N₂ -fixing N. azollae, we conclude, dominated the microbiome of Azolla ferns. The persistent but less abundant heterotrophic Rhizobiales bacteria possibly contributed to lowering O₂ levels in leaf pockets but did not release detectable amounts of the strong greenhouse gas N₂ O.