Rooting depth and xylem vulnerability are independent woody plant traits jointly selected by aridity, seasonality, and water table depth.

Evolutionary radiations of woody taxa within arid environments were made possible by multiple trait innovations including deep roots and embolism-resistant xylem, but little is known about how these traits have coevolved across the phylogeny of woody plants or how they jointly influence the distribution of species. We synthesized global trait and vegetation plot datasets to examine how rooting depth and xylem vulnerability across 188 woody plant species interact with aridity, precipitation seasonality, and water table depth to influence species occurrence probabilities across all biomes. Xylem resistance to embolism and rooting depth are independent woody plant traits that do not exhibit an interspecific trade-off. Resistant xylem and deep roots increase occurrence probabilities in arid, seasonal climates over deep water tables. Resistant xylem and shallow roots increase occurrence probabilities in arid, nonseasonal climates over deep water tables. Vulnerable xylem and deep roots increase occurrence probabilities in arid, nonseasonal climates over shallow water tables. Lastly, vulnerable xylem and shallow roots increase occurrence probabilities in humid climates. Each combination of trait values optimizes occurrence probabilities in unique environmental conditions. Responses of deeply rooted vegetation may be buffered if evaporative demand changes faster than water table depth under climate change.

Areas of global importance for conserving terrestrial biodiversity, carbon and water.

To meet the ambitious objectives of biodiversity and climate conventions, the international community requires clarity on how these objectives can be operationalized spatially and how multiple targets can be pursued concurrently. To support goal setting and the implementation of international strategies and action plans, spatial guidance is needed to identify which land areas have the potential to generate the greatest synergies between conserving biodiversity and nature's contributions to people. Here we present results from a joint optimization that minimizes the number of threatened species, maximizes carbon retention and water quality regulation, and ranks terrestrial conservation priorities globally. We found that selecting the top-ranked 30% and 50% of terrestrial land area would conserve respectively 60.7% and 85.3% of the estimated total carbon stock and 66% and 89.8% of all clean water, in addition to meeting conservation targets for 57.9% and 79% of all species considered. Our data and prioritization further suggest that adequately conserving all species considered (vertebrates and plants) would require giving conservation attention to ~70% of the terrestrial land surface. If priority was given to biodiversity only, managing 30% of optimally located land area for conservation may be sufficient to meet conservation targets for 81.3% of the terrestrial plant and vertebrate species considered. Our results provide a global assessment of where land could be optimally managed for conservation. We discuss how such a spatial prioritization framework can support the implementation of the biodiversity and climate conventions.

The Importance of Protein Phosphorylation for Signaling and Metabolism in Response to Diel Light Cycling and Nutrient Availability in a Marine Diatom.

Diatoms are major contributors to global primary production and their populations in the modern oceans are affected by availability of iron, nitrogen, phosphate, silica, and other trace metals, vitamins, and infochemicals. However, little is known about the role of phosphorylation in diatoms and its role in regulation and signaling. We report a total of 2759 phosphorylation sites on 1502 proteins detected in Phaeodactylum tricornutum. Conditionally phosphorylated peptides were detected at low iron (n = 108), during the diel cycle (n = 149), and due to nitrogen availability (n = 137). Through a multi-omic comparison of transcript, protein, phosphorylation, and protein homology, we identify numerous proteins and key cellular processes that are likely under control of phospho-regulation. We show that phosphorylation regulates: (1) carbon retrenchment and reallocation during growth under low iron, (2) carbon flux towards lipid biosynthesis after the lights turn on, (3) coordination of transcription and translation over the diel cycle and (4) in response to nitrogen depletion. We also uncover phosphorylation sites for proteins that play major roles in diatom Fe sensing and utilization, including flavodoxin and phytotransferrin (ISIP2A), as well as identify phospho-regulated stress proteins and kinases. These findings provide much needed insight into the roles of protein phosphorylation in diel cycling and nutrient sensing in diatoms.

Multiplexed Knockouts in the Model Diatom Phaeodactylum by Episomal Delivery of a Selectable Cas9.

Marine diatoms are eukaryotic microalgae that play significant ecological and biogeochemical roles in oceans. They also have significant potential as organismal platforms for exploitation to address biotechnological and industrial goals. In order to address both modes of research, sophisticated molecular and genetic tools are required. We presented here new and improved methodologies for introducing CRISPR-Cas9 to the model diatom Phaeodactylum tricornutum cells and a streamlined protocol for genotyping mutant cell lines with previously unknown phenotypes. First, bacterial-conjugation was optimized for the delivery of Cas9 by transcriptionally fusing Cas9 to a selectable marker by the 2A peptide. An episome cloning strategy using both negative and positive selection was developed to streamline CRISPR-episome assembly. Next, cell line picking and genotyping strategies, that utilize manual sequencing curation, TIDE sequencing analysis, and a T7 endonuclease assay, were developed to shorten the time required to generate mutants. Following this new experimental pipeline, both single-gene and two-gene knockout cell lines were generated at mutagenesis efficiencies of 48% and 25%, respectively. Lastly, a protocol for precise gene insertions via CRISPR-Cas9 targeting was developed using particle-bombardment transformation methods. Overall, the novel Cas9 episome design and improved genotyping methods presented here allow for quick and easy genotyping and isolation of Phaeodactylum mutant cell lines (less than 3 weeks) without relying on a known phenotype to screen for mutants.

A regional-scale assessment of using metabolic scaling theory to predict ecosystem properties.

Metabolic scaling theory (MST) is one of ecology's most high-profile general models and can be used to link size distributions and productivity in forest systems. Much of MST's foundation is based on size distributions following a power law function with a scaling exponent of -2, a property assumed to be consistent in steady-state ecosystems. We tested the theory's generality by comparing actual size distributions with those predicted using MST parameters assumed to be general. We then used environmental variables and functional traits to explain deviation from theoretical expectations. Finally, we compared values of relative productivity predicted using MST with a remote-sensed measure of productivity. We found that fire-prone heath communities deviated from MST-predicted size distributions, whereas fire-sensitive rainforests largely agreed with the theory. Scaling exponents ranged from -1.4 to -5.3. Deviation from the power law assumption was best explained by specific leaf area, which varies along fire frequency and moisture gradients. While MST may hold in low-disturbance systems, we show that it cannot be applied under many environmental contexts. The theory should remain general, but understanding the factors driving deviation from MST and subsequent refinements is required if it is to be applied robustly across larger scales.

Evolution and regulation of nitrogen flux through compartmentalized metabolic networks in a marine diatom.

Diatoms outcompete other phytoplankton for nitrate, yet little is known about the mechanisms underpinning this ability. Genomes and genome-enabled studies have shown that diatoms possess unique features of nitrogen metabolism however, the implications for nutrient utilization and growth are poorly understood. Using a combination of transcriptomics, proteomics, metabolomics, fluxomics, and flux balance analysis to examine short-term shifts in nitrogen utilization in the model pennate diatom in Phaeodactylum tricornutum, we obtained a systems-level understanding of assimilation and intracellular distribution of nitrogen. Chloroplasts and mitochondria are energetically integrated at the critical intersection of carbon and nitrogen metabolism in diatoms. Pathways involved in this integration are organelle-localized GS-GOGAT cycles, aspartate and alanine systems for amino moiety exchange, and a split-organelle arginine biosynthesis pathway that clarifies the role of the diatom urea cycle. This unique configuration allows diatoms to efficiently adjust to changing nitrogen status, conferring an ecological advantage over other phytoplankton taxa.

Nitrate Reductase Knockout Uncouples Nitrate Transport from Nitrate Assimilation and Drives Repartitioning of Carbon Flux in a Model Pennate Diatom.

The ecological prominence of diatoms in the ocean environment largely results from their superior competitive ability for dissolved nitrate (NO3-). To investigate the cellular and genetic basis of diatom NO3- assimilation, we generated a knockout in the nitrate reductase gene (NR-KO) of the model pennate diatom Phaeodactylum tricornutum In NR-KO cells, N-assimilation was abolished although NO3- transport remained intact. Unassimilated NO3- accumulated in NR-KO cells, resulting in swelling and associated changes in biochemical composition and physiology. Elevated expression of genes encoding putative vacuolar NO3- chloride channel transporters plus electron micrographs indicating enlarged vacuoles suggested vacuolar storage of NO3- Triacylglycerol concentrations in the NR-KO cells increased immediately following the addition of NO3-, and these increases coincided with elevated gene expression of key triacylglycerol biosynthesis components. Simultaneously, induction of transcripts encoding proteins involved in thylakoid membrane lipid recycling suggested more abrupt repartitioning of carbon resources in NR-KO cells compared with the wild type. Conversely, ribosomal structure and photosystem genes were immediately deactivated in NR-KO cells following NO3- addition, followed within hours by deactivation of genes encoding enzymes for chlorophyll biosynthesis and carbon fixation and metabolism. N-assimilation pathway genes respond uniquely, apparently induced simultaneously by both NO3- replete and deplete conditions.

Using single nuclei for RNA-seq to capture the transcriptome of postmortem neurons.

A protocol is described for sequencing the transcriptome of a cell nucleus. Nuclei are isolated from specimens and sorted by FACS, cDNA libraries are constructed and RNA-seq is performed, followed by data analysis. Some steps follow published methods (Smart-seq2 for cDNA synthesis and Nextera XT barcoded library preparation) and are not described in detail here. Previous single-cell approaches for RNA-seq from tissues include cell dissociation using protease treatment at 30 °C, which is known to alter the transcriptome. We isolate nuclei at 4 °C from tissue homogenates, which cause minimal damage. Nuclear transcriptomes can be obtained from postmortem human brain tissue stored at -80 °C, making brain archives accessible for RNA-seq from individual neurons. The method also allows investigation of biological features unique to nuclei, such as enrichment of certain transcripts and precursors of some noncoding RNAs. By following this procedure, it takes about 4 d to construct cDNA libraries that are ready for sequencing.

Inactivation of Phaeodactylum tricornutum urease gene using transcription activator-like effector nuclease-based targeted mutagenesis.

Diatoms are unicellular photosynthetic algae with promise for green production of fuels and other chemicals. Recent genome-editing techniques have greatly improved the potential of many eukaryotic genetic systems, including diatoms, to enable knowledge-based studies and bioengineering. Using a new technique, transcription activator-like effector nucleases (TALENs), the gene encoding the urease enzyme in the model diatom, Phaeodactylum tricornutum, was targeted for interruption. The knockout cassette was identified within the urease gene by PCR and Southern blot analyses of genomic DNA. The lack of urease protein was confirmed by Western blot analyses in mutant cell lines that were unable to grow on urea as the sole nitrogen source. Untargeted metabolomic analysis revealed a build-up of urea, arginine and ornithine in the urease knockout lines. All three intermediate metabolites are upstream of the urease reaction within the urea cycle, suggesting a disruption of the cycle despite urea production. Numerous high carbon metabolites were enriched in the mutant, implying a breakdown of cellular C and N repartitioning. The presented method improves the molecular toolkit for diatoms and clarifies the role of urease in the urea cycle.

High-resolution DNA melt-curve analysis for cost-effective mass screening of pairwise species interactions.

Ecological studies of pairwise interactions are constrained by the methods available for rapid species identification of the interacting organisms. The resolution of data required to characterize species interaction networks at multiple spatio-temporal scales can be intensive, and therefore laborious and costly to collect. We explore the utility of high-resolution DNA melt-curve analysis (HRM) as a rapid species identification method. An approach was developed to identify organisms at the pairwise interaction level, with particular application to cryptic species interactions that are traditionally difficult to study. Here, we selected a challenging application; to identify the presence/absence of pathogenic fungi (Sporothrix inflata, Ophiostoma nigrocarpum and Ophiostoma galeiforme) transported by bark beetle vectors (Hylastes ater and Hylurgus ligniperda). The technique was able to distinguish between different species of DNA within a single, pooled sample. In test applications, HRM was effective in the mass screening and identification of pathogenic fungal species carried by many individual bark beetle vectors (n = 455 beetles screened) across large geographic scales. For two of the fungal species, there was no difference in the frequency of association with either of their vectors, but for the third fungal species there was a shift in vector-pathogen associations across locations. This technique allows rapid, mass screening and characterization of species interactions at a fraction of the time and cost of traditional methods. It is anticipated that this method can be readily applied to explore other cryptic species interactions, or other studies requiring rapid generation of large data sets and/or high-throughput efficiency.

An experimental test of insect-mediated colonisation of damaged Pinus radiata trees by sapstain fungi.

Vector-pathogen dynamics play a central role in understanding tree health and forest dynamics. There is substantial evidence that bark beetles act as spore vectors for many species of fungi that cause 'sapstain' discolouration of damaged trees and timber. However, the direct quantitative link between vector-mediated spore dispersal and subsequent sapstain colonisation of wood is not fully understood. Here, we used caged versus uncaged experimental logs to test whether the exclusion of bark beetles quantitatively alters the distribution and intensity of sapstain fungal spread within damaged trees. Using generalised linear mixed models, we tested the effect of bark beetle exclusion on sapstain intensity within and among cut logs at two plantation forest sites. Overall, sapstain was found on all logs regardless of caging treatment, indicating that sapstain colonisation can occur (to some degree) without arthropod vectors, probably via wind, rain-splash and, potentially, latent endophytic development. This was supported by the dominance of Diplodia pinea in fungal isolations taken from trees felled at the site, as this fungal species is known to disperse independently of bark beetles. However, the intensity of sapstain within and among experimental logs was significantly greater in uncaged than in caged logs, where beetle colonisation was significantly greater. This appeared to be driven by a significant within-log association between the intensity of staining and the intensity of beetle, and other arthropod, tunnelling and feeding activities. Taken together, these results strongly suggest that the dominant mechanism underlying the role of bark beetles in sapstain development in this study system is not vector-mediated spore dispersal, per se, but rather the facilitation of spore entry and hyphal development through tunnelling and feeding activities. We discuss the implications of these findings for forest management and the effective salvage-harvest of trees damaged by stochastic climate events such as storm and fire damage.

Elimination of manganese(II,III) oxidation in Pseudomonas putida GB-1 by a double knockout of two putative multicopper oxidase genes.

Bacterial manganese(II) oxidation impacts the redox cycling of Mn, other elements, and compounds in the environment; therefore, it is important to understand the mechanisms of and enzymes responsible for Mn(II) oxidation. In several Mn(II)-oxidizing organisms, the identified Mn(II) oxidase belongs to either the multicopper oxidase (MCO) or the heme peroxidase family of proteins. However, the identity of the oxidase in Pseudomonas putida GB-1 has long remained unknown. To identify the P. putida GB-1 oxidase, we searched its genome and found several homologues of known or suspected Mn(II) oxidase-encoding genes (mnxG, mofA, moxA, and mopA). To narrow this list, we assumed that the Mn(II) oxidase gene would be conserved among Mn(II)-oxidizing pseudomonads but not in nonoxidizers and performed a genome comparison to 11 Pseudomonas species. We further assumed that the oxidase gene would be regulated by MnxR, a transcription factor required for Mn(II) oxidation. Two loci met all these criteria: PputGB1_2447, which encodes an MCO homologous to MnxG, and PputGB1_2665, which encodes an MCO with very low homology to MofA. In-frame deletions of each locus resulted in strains that retained some ability to oxidize Mn(II) or Mn(III); loss of oxidation was attained only upon deletion of both genes. These results suggest that PputGB1_2447 and PputGB1_2665 encode two MCOs that are independently capable of oxidizing both Mn(II) and Mn(III). The purpose of this redundancy is unclear; however, differences in oxidation phenotype for the single mutants suggest specialization in function for the two enzymes.

Evaluation of GFP tag as a screening reporter in directed evolution of a hyperthermophilic beta-glucosidase.

By applying a directed evolution methodology specific enzymatic characteristics can be enhanced, but to select mutants of interest from a large mutant bank, this approach requires high throughput screening and facile selection. To facilitate such primary screening of enhanced clones, an expression system was tested that uses a green fluorescent protein (GFP) tag from Aequorea victoria linked to the enzyme of interest. As GFP's fluorescence is readily measured, and as there is a 1:1 molar correlation between the target protein and GFP, the concept proposed was to determine whether GFP could facilitate primary screening of error-prone PCR (EPP) clones. For this purpose a thermostable beta-glucosidase (BglA) from Fervidobacterium sp. was used as a model enzyme. A vector expressing the chimeric protein BglA-GFP-6XHis was constructed and the fusion protein purified and characterized. When compared to the native proteins, the components of the fusion displayed modified characteristics, such as enhanced GFP thermostability and a higher BglA optimum temperature. Clones carrying mutant BglA proteins obtained by EPP, were screened based on the BglA/GFP activity ratio. Purified tagged enzymes from selected clones resulted in modified substrate specificity.

Genomic insights into Mn(II) oxidation by the marine alphaproteobacterium Aurantimonas sp. strain SI85-9A1.

Microbial Mn(II) oxidation has important biogeochemical consequences in marine, freshwater, and terrestrial environments, but many aspects of the physiology and biochemistry of this process remain obscure. Here, we report genomic insights into Mn(II) oxidation by the marine alphaproteobacterium Aurantimonas sp. strain SI85-9A1, isolated from the oxic/anoxic interface of a stratified fjord. The SI85-9A1 genome harbors the genetic potential for metabolic versatility, with genes for organoheterotrophy, methylotrophy, oxidation of sulfur and carbon monoxide, the ability to grow over a wide range of O(2) concentrations (including microaerobic conditions), and the complete Calvin cycle for carbon fixation. Although no growth could be detected under autotrophic conditions with Mn(II) as the sole electron donor, cultures of SI85-9A1 grown on glycerol are dramatically stimulated by addition of Mn(II), suggesting an energetic benefit from Mn(II) oxidation. A putative Mn(II) oxidase is encoded by duplicated multicopper oxidase genes that have a complex evolutionary history including multiple gene duplication, loss, and ancient horizontal transfer events. The Mn(II) oxidase was most abundant in the extracellular fraction, where it cooccurs with a putative hemolysin-type Ca(2+)-binding peroxidase. Regulatory elements governing the cellular response to Fe and Mn concentration were identified, and 39 targets of these regulators were detected. The putative Mn(II) oxidase genes were not among the predicted targets, indicating that regulation of Mn(II) oxidation is controlled by other factors yet to be identified. Overall, our results provide novel insights into the physiology and biochemistry of Mn(II) oxidation and reveal a genome specialized for life at the oxic/anoxic interface.

Geomicrobiology of manganese(II) oxidation.

Mn(II)-oxidizing microbes have an integral role in the biogeochemical cycling of manganese, iron, nitrogen, carbon, sulfur, and several nutrients and trace metals. There is great interest in mechanistically understanding these cycles and defining the importance of Mn(II)-oxidizing bacteria in modern and ancient geochemical environments. Linking Mn(II) oxidation to cellular function, although still enigmatic, continues to drive efforts to characterize manganese biomineralization. Recently, complexed-Mn(III) has been shown to be a transient intermediate in Mn(II) oxidation to Mn(IV), suggesting that the reaction might involve a unique multicopper oxidase system capable of a two-electron oxidation of the substrate. In biogenic and abiotic synthesis experiments, the application of synchrotron-based X-ray scattering and spectroscopic techniques has significantly increased our understanding of the oxidation state and relatively amorphous structure (i.e. delta-MnO(2)-like) of biogenic oxides, providing a new blueprint for the structural signature of biogenic Mn oxides.

Improved catalytic efficiency and active site modification of 1,4-beta-D-glucan glucohydrolase A from Thermotoga neapolitana by directed evolution.

Thermotoga neapolitana 1,4-beta-d-glucan glucohydrolase A preferentially hydrolyzes cello-oligomers, such as cellotetraose, releasing single glucose moieties from the reducing end of the cello-oligosaccharide chain. Using directed evolution techniques of error-prone PCR and mutant library screening, a variant glucan glucohydrolase has been isolated that hydrolyzes the disaccharide, cellobiose, at a 31% greater rate than its wild type (WT) predecessor. The mutant library, expressed in Escherichia coli, was screened at 85 degrees C for increased hydrolysis of cellobiose, a native substrate rather than a chromogenic analog, using a continuous, thermostable coupled enzyme assay. The V(max) for the mutant was 108 +/- 3 units mg(-1), whereas that of the WT was 75 +/- 2 units mg(-1). The K(m) for both proteins was nearly the same. The k(cat) for the new enzyme increased by 31% and its catalytic efficiency (k(cat)/K(m)) for cellobiose also rose by 31% as compared with the parent. The nucleotide sequence of two positive clones and two null clones identified 11 single base shifts. The nucleotide transition in the most active clone caused an isoleucine to threonine amino acid substitution at position 170. Structural models for I170T and WT proteins were derived by sequence homology with Protein Data Bank code 1BGA from Paenibacillus polymyxa. Analysis of the WT and I170T model structures indicated that the substitution in the mutant enzyme repositioned the conserved catalytic residue Asn-163 and reconfigured entry to the active site.

Thermostable continuous coupled assay for measuring glucose using glucokinase and glucose-6-phosphate dehydrogenase from the marine hyperthermophile Thermotoga maritima.

A novel, thermostable adaptation of the coupled-enzyme assay for monitoring glucose concentrations was developed for an optimal temperature of 85 degrees C. This is the first report of a thermostable glucostat from a marine hyperthermophile. The continuous assay, using glucokinase (Glk) and glucose-6-phosphate dehydrogenase (Gpd) from Thermotoga maritima, demonstrated robust activity over a range of temperatures (75-90 degrees C) and pH values (6.8- 8.5). Purified glucokinase had a monomeric molecular mass of 33.8kDa while that of glucose-6-phosphate dehydrogenase (D-glucose 6-phosphate:NADP oxidoreductase) was 57.5kDa. The high-temperature assay provided a method for directly assaying the activity of another hyperthermophilic enzyme, 1,4-beta-D-glucan glucohydrolase (GghA) from Thermotoga neapolitana. To provide a benchmark for protein-engineering experiments involving GghA, a three-enzyme continuous assay (performed at 85 degrees C), linking wild-type GghA, Glk, and Gpd, measured glucose produced from GghA's hydrolysis of cellobiose, one of GghA's secondary substrates. The assay established the kinetic behavior of wild-type GghA toward cellobiose and was used to screen for changes in the catalytic efficiency of variant GghA(s) induced by random mutagenesis. The assay's development will allow high-throughput screening of other thermostable glucose-producing enzymes, including those applicable to commercial biomass conversion.