ASU Scholarship Showcase

Cyanobacteria are considered good models for biohydrogen production because they are relatively simple organisms with a demonstrable ability to generate H₂ under certain physiological conditions. However, most produce only little H₂, revert readily to H₂ consumption, and suffer from hydrogenase sensitivity to O₂. Strains of the cyanobacteria Lyngbya aestuarii and Microcoleus chthonoplastes obtained from marine intertidal cyanobacterial mats were recently found to display much better H₂ production potential. Because of their ecological origin in environments that become quickly anoxic in the dark, we hypothesized that this differential ability may have evolved to serve a role in the fermentation of the photosynthate. Here we show that, when forced to ferment internal substrate, these cyanobacteria display desirable characteristics of physiological H₂ production. Among them, the strain L. aestuarii BL J had the fastest specific rates and attained the highest H₂ concentrations during fermentation of photosynthate, which proceeded via a mixed acid fermentation pathway to yield acetate, ethanol, lactate, H₂, CO₂, and pyruvate. Contrary to expectations, the H₂ yield per mole of glucose was only average compared to that of other cyanobacteria. Thermodynamic analyses point to the use of electron donors more electronegative than NAD(P)H in Lyngbya hydrogenases as the basis for its strong H₂ production ability. In any event, the high specific rates and H₂ concentrations coupled with the lack of reversibility of the enzyme, at the expense of internal, photosynthetically generated reductants, makes L. aestuarii BL J and/or its enzymes, a potentially feasible platform for large-scale H₂ production.

Biological Soil Crusts (BSCs) are organosedimentary assemblages comprised of microbes and minerals in topsoil of terrestrial environments. BSCs strongly impact soil quality in dryland ecosystems (e.g., soil structure and nutrient yields) due to pioneer species such as Microcoleus vaginatus; phototrophs that produce filaments that bind the soil together, and support an array of heterotrophic microorganisms. These microorganisms in turn contribute to soil stability and biogeochemistry of BSCs. Non-cyanobacterial populations of BSCs are less well known than cyanobacterial populations. Therefore, we attempted to isolate a broad range of numerically significant and phylogenetically representative BSC aerobic heterotrophs. Combining simple pre-treatments (hydration of BSCs under dark and light) and isolation strategies (media with varying nutrient availability and protection from oxidative stress) we recovered 402 bacterial and one fungal isolate in axenic culture, which comprised 116 phylotypes (at 97% 16S rRNA gene sequence homology), 115 bacterial and one fungal. Each medium enriched a mostly distinct subset of phylotypes, and cultivated phylotypes varied due to the BSC pre-treatment. The fraction of the total phylotype diversity isolated, weighted by relative abundance in the community, was determined by the overlap between isolate sequences and OTUs reconstructed from metagenome or metatranscriptome reads. Together, more than 8% of relative abundance of OTUs in the metagenome was represented by our isolates, a cultivation efficiency much larger than typically expected from most soils. We conclude that simple cultivation procedures combined with specific pre-treatment of samples afford a significant reduction in the culturability gap, enabling physiological and metabolic assays that rely on ecologically relevant axenic cultures.

Background: The extracellular sunscreen scytonemin is the most common and widespread indole-alkaloid among cyanobacteria. Previous research using the cyanobacterium Nostoc punctiforme ATCC 29133 revealed a unique 18-gene cluster (NpR1276 to NpR1259 in the N. punctiforme genome) involved in the biosynthesis of scytonemin. We provide further genomic characterization of these genes in N. punctiforme and extend it to homologous regions in other cyanobacteria.

Results: Six putative genes in the scytonemin gene cluster (NpR1276 to NpR1271 in the N. punctiforme genome), with no previously known protein function and annotated in this study as scyA to scyF, are likely involved in the assembly of scytonemin from central metabolites, based on genetic, biochemical, and sequence similarity evidence. Also in this cluster are redundant copies of genes encoding for aromatic amino acid biosynthetic enzymes. These can theoretically lead to tryptophan and the tyrosine precursor, p-hydroxyphenylpyruvate, (expected biosynthetic precursors of scytonemin) from end products of the shikimic acid pathway. Redundant copies of the genes coding for the key regulatory and rate-limiting enzymes of the shikimic acid pathway are found there as well. We identified four other cyanobacterial strains containing orthologues of all of these genes, three of them by database searches (Lyngbya PCC 8106, Anabaena PCC 7120, and Nodularia CCY 9414) and one by targeted sequencing (Chlorogloeopsis sp. strain Cgs-089; CCMEE 5094). Genomic comparisons revealed that most scytonemin-related genes were highly conserved among strains and that two additional conserved clusters, NpF5232 to NpF5236 and a putative two-component regulatory system (NpF1278 and NpF1277), are likely involved in scytonemin biosynthesis and regulation, respectively, on the basis of conservation and location. Since many of the protein product sequences for the newly described genes, including ScyD, ScyE, and ScyF, have export signal domains, while others have putative transmembrane domains, it can be inferred that scytonemin biosynthesis is compartmentalized within the cell. Basic structural monomer synthesis and initial condensation are most likely cytoplasmic, while later reactions are predicted to be periplasmic.

Conclusion: We show that scytonemin biosynthetic genes are highly conserved among evolutionarily diverse strains, likely include more genes than previously determined, and are predicted to involve compartmentalization of the biosynthetic pathway in the cell, an unusual trait for prokaryotes.

N₂ fixation and ammonia oxidation (AO) are the two most important processes in the nitrogen (N) cycle of biological soil crusts (BSCs). We studied the short-term response of acetylene reduction assay (ARA) rates, an indicator of potential N₂ fixation, and AO rates to temperature (T, -5°C to 35°C) in BSC of different successional stages along the BSC ecological succession and geographic origin (hot Chihuahuan and cooler Great Basin deserts). ARA in all BSCs increased with T until saturation occurred between 15 and 20°C, and declined at 30–35°C. Culture studies using cyanobacteria isolated from these crusts indicated that the saturating effect was traceable to their inability to grow well diazotrophically within the high temperature range. Below saturation, temperature response was exponential, with Q₁₀ significantly different in the two areas (~ 5 for Great Basin BSCs; 2–3 for Chihuahuan BSCs), but similar between the two successional stages. However, in contrast to ARA, AO showed a steady increase to 30–35°C in Great Basin, and Chihuhuan BSCs showed no inhibition at any tested temperature. The T response of AO also differed significantly between Great Basin (Q₁₀ of 4.5–4.8) and Chihuahuan (Q₁₀ of 2.4–2.6) BSCs, but not between successional stages. Response of ARA rates to T did not differ from that of AO in either desert. Thus, while both processes scaled to T in unison until 20°C, they separated to an increasing degree at higher temperature. As future warming is likely to occur in the regions where BSCs are often the dominant living cover, this predicted decoupling is expected to result in higher proportion of nitrates in soil relative to ammonium. As nitrate is more easily lost as leachate or to be reduced to gaseous forms, this could mean a depletion of soil N over large landscapes globally.

Endolithic microbial communities are prominent features of intertidal marine habitats, where they colonize a variety of substrates, contributing to their erosion. Almost 2 centuries worth of naturalistic studies focused on a few true-boring (euendolithic) phototrophs, but substrate preference has received little attention. The Isla de Mona (Puerto Rico) intertidal zone offers a unique setting to investigate substrate specificity of endolithic communities since various phosphate rock, limestone and dolostone outcrops occur there. High-throughput 16S rDNA genetic sampling, enhanced by targeted cultivation, revealed that, while euendolithic cyanobacteria were dominant operational taxonomic units (OTUs), the communities were invariably of high diversity, well beyond that reported in traditional studies and implying an unexpected metabolic complexity potentially contributed by secondary colonizers. While the overall community composition did not show differences traceable to the nature of the mineral substrate, we detected specialization among particular euendolithic cyanobacterial clades towards the type of substrate they excavate but only at the OTU phylogenetic level, implying that close relatives have specialized recurrently into particular substrates. The cationic mineral component was determinant in this preference, suggesting the existence in nature of alternatives to the boring mechanism described in culture that is based exclusively on transcellular calcium transport.

The sensitivity of Earth’s wetlands to observed shifts in global precipitation and temperature patterns and their ability to produce large quantities of methane gas are key global change questions. We present a microwave satellite-based approach for mapping fractional surface water (FW) globally at 25-km resolution. The approach employs a land cover-supported, atmospherically-corrected dynamic mixture model applied to 20+ years (1992–2013) of combined, daily, passive/active microwave remote sensing data. The resulting product, known as Surface Water Microwave Product Series (SWAMPS), shows strong microwave sensitivity to sub-grid scale open water and inundated wetlands comprising open plant canopies. SWAMPS’ FW compares favorably (R² = 91%–94%) with higher-resolution, global-scale maps of open water from MODIS and SRTM-MOD44W. Correspondence of SWAMPS with open water and wetland products from satellite SAR in Alaska and the Amazon deteriorates when exposed wetlands or inundated forests captured by the SAR products were added to the open water fraction reflecting SWAMPS’ inability to detect water underneath the soil surface or beneath closed forest canopies. Except for a brief period of drying during the first 4 years of observation, the inundation extent for the global domain excluding the coast was largely stable. Regionally, inundation in North America is advancing while inundation is on the retreat in Tropical Africa and North Eurasia. SWAMPS provides a consistent and long-term global record of daily FW dynamics, with documented accuracies suitable for hydrologic assessment and global change-related investigations.

The filamentous, non-heterocystous cyanobacterium Lyngbya aestuarii is an important contributor to marine intertidal microbial mats system worldwide. The recent isolate L. aestuarii BL J, is an unusually powerful hydrogen producer. Here we report a morphological, ultrastructural, and genomic characterization of this strain to set the basis for future systems studies and applications of this organism. The filaments contain circa 17 μm wide trichomes, composed of stacked disk-like short cells (2 μm long), encased in a prominent, laminated exopolysaccharide sheath. Cellular division occurs by transversal centripetal growth of cross-walls, where several rounds of division proceed simultaneously. Filament division occurs by cell self-immolation of one or groups of cells (necridial cells) at the breakage point. Short, sheath-less, motile filaments (hormogonia) are also formed. Morphologically and phylogenetically L. aestuarii belongs to a clade of important cyanobacteria that include members of the marine Trichodesmiun and Hydrocoleum genera, as well as terrestrial Microcoleus vaginatus strains, and alkalyphilic strains of Arthrospira. A draft genome of strain BL J was compared to those of other cyanobacteria in order to ascertain some of its ecological constraints and biotechnological potential.

The genome had an average GC content of 41.1%. Of the 6.87 Mb sequenced, 6.44 Mb was present as large contigs (>10,000 bp). It contained 6515 putative protein-encoding genes, of which, 43% encode proteins of known functional role, 26% corresponded to proteins with domain or family assignments, 19.6% encode conserved hypothetical proteins, and 11.3% encode apparently unique hypothetical proteins. The strain's genome reveals its adaptations to a life of exposure to intense solar radiation and desiccation. It likely employs the storage compounds, glycogen, and cyanophycin but no polyhydroxyalkanoates, and can produce the osmolytes, trehalose, and glycine betaine. According to its genome, BL J strain also has the potential to produce a plethora of products of biotechnological interest such as Curacin A, Barbamide, Hemolysin-type calcium-binding toxin, the suncreens scytonemin, and mycosporines, as well as heptadecane and pentadecane alkanes. With respect to hydrogen production, initial comparisons of the genetic architecture and sequence of relevant genes and loci, and a comparative model of protein structure of the NiFe bidirectional hydrogenase, did not reveal conspicuous differences that could explain its unusual hydrogen producing capacity.

A warming climate is altering land-atmosphere exchanges of carbon, with a potential for increased vegetation productivity as well as the mobilization of permafrost soil carbon stores. Here we investigate land-atmosphere carbon dioxide (CO₂) cycling through analysis of net ecosystem productivity (NEP) and its component fluxes of gross primary productivity (GPP) and ecosystem respiration (ER) and soil carbon residence time, simulated by a set of land surface models (LSMs) over a region spanning the drainage basin of Northern Eurasia. The retrospective simulations cover the period 1960–2009 at 0.5° resolution, which is a scale common among many global carbon and climate model simulations. Model performance benchmarks were drawn from comparisons against both observed CO₂ fluxes derived from site-based eddy covariance measurements as well as regional-scale GPP estimates based on satellite remote-sensing data.

The site-based comparisons depict a tendency for overestimates in GPP and ER for several of the models, particularly at the two sites to the south. For several models the spatial pattern in GPP explains less than half the variance in the MODIS MOD17 GPP product. Across the models NEP increases by as little as 0.01 to as much as 0.79 g C m^-2 yr^-2, equivalent to 3 to 340 % of the respective model means, over the analysis period. For the multimodel average the increase is 135 % of the mean from the first to last 10 years of record (1960–1969 vs. 2000–2009), with a weakening CO₂ sink over the latter decades. Vegetation net primary productivity increased by 8 to 30 % from the first to last 10 years, contributing to soil carbon storage gains. The range in regional mean NEP among the group is twice the multimodel mean, indicative of the uncertainty in CO₂ sink strength.

The models simulate that inputs to the soil carbon pool exceeded losses, resulting in a net soil carbon gain amid a decrease in residence time. Our analysis points to improvements in model elements controlling vegetation productivity and soil respiration as being needed for reducing uncertainty in land-atmosphere CO₂ exchange. These advances will require collection of new field data on vegetation and soil dynamics, the development of benchmarking data sets from measurements and remote-sensing observations, and investments in future model development and intercomparison studies.

Soil temperature (T_s) change is a key indicator of the dynamics of permafrost. On seasonal and interannual timescales, the variability of T_s determines the active-layer depth, which regulates hydrological soil properties and biogeochemical processes. On the multi-decadal scale, increasing T_s not only drives permafrost thaw/retreat but can also trigger and accelerate the decomposition of soil organic carbon. The magnitude of permafrost carbon feedbacks is thus closely linked to the rate of change of soil thermal regimes. In this study, we used nine process-based ecosystem models with permafrost processes, all forced by different observation-based climate forcing during the period 1960–2000, to characterize the warming rate of T_s in permafrost regions. There is a large spread of T_s trends at 20 cm depth across the models, with trend values ranging from 0.010 ± 0.003 to 0.031 ± 0.005 °C yr^-1. Most models show smaller increase in T_s with increasing depth. Air temperature (Tsub>a) and longwave downward radiation (LWDR) are the main drivers of T_s trends, but their relative contributions differ amongst the models. Different trends of LWDR used in the forcing of models can explain 61 % of their differences in T_s trends, while trends of T_a only explain 5 % of the differences in T_s trends. Uncertain climate forcing contributes a larger uncertainty in T_s trends (0.021 ± 0.008 °C yr^-1, mean ± standard deviation) than the uncertainty of model structure (0.012 ± 0.001 °C yr^-1), diagnosed from the range of response between different models, normalized to the same forcing. In addition, the loss rate of near-surface permafrost area, defined as total area where the maximum seasonal active-layer thickness (ALT) is less than 3 m loss rate, is found to be significantly correlated with the magnitude of the trends of T_s at 1 m depth across the models (R = −0.85, P = 0.003), but not with the initial total near-surface permafrost area (R = −0.30, P = 0.438). The sensitivity of the total boreal near-surface permafrost area to T_s at 1 m is estimated to be of −2.80 ± 0.67 million km²°C^-1. Finally, by using two long-term LWDR data sets and relationships between trends of LWDR and T_s across models, we infer an observation-constrained total boreal near-surface permafrost area decrease comprising between 39 ± 14  ×  10³ and 75 ± 14  ×  10³km²yr^-1 from 1960 to 2000. This corresponds to 9–18 % degradation of the current permafrost area.

A realistic simulation of snow cover and its thermal properties are important for accurate modelling of permafrost. We analyse simulated relationships between air and near-surface (20  cm) soil temperatures in the Northern Hemisphere permafrost region during winter, with a particular focus on snow insulation effects in nine land surface models, and compare them with observations from 268 Russian stations. There are large cross-model differences in the simulated differences between near-surface soil and air temperatures (ΔT; 3 to 14 °C), in the sensitivity of soil-to-air temperature (0.13 to 0.96 °C °^C-1), and in the relationship between ΔT and snow depth. The observed relationship between ΔT and snow depth can be used as a metric to evaluate the effects of each model's representation of snow insulation, hence guide improvements to the model's conceptual structure and process parameterisations. Models with better performance apply multilayer snow schemes and consider complex snow processes. Some models show poor performance in representing snow insulation due to underestimation of snow depth and/or overestimation of snow conductivity. Generally, models identified as most acceptable with respect to snow insulation simulate reasonable areas of near-surface permafrost (13.19 to 15.77 million  km²). However, there is not a simple relationship between the sophistication of the snow insulation in the acceptable models and the simulated area of Northern Hemisphere near-surface permafrost, because several other factors, such as soil depth used in the models, the treatment of soil organic matter content, hydrology and vegetation cover, also affect the simulated permafrost distribution.