Silicon cyling along the land-ocean continuum

The alteration of the global environment by human activities is so widespread that scientists argue we've entered a new geologic epoch known as the Anthropocene. This dissertation examines the impact of human activities on biogeochemical cycling at the land-sea interface. I focus primarily on the role of land use/land cover (LULC) and coastal nutrient enrichment on silicon (Si) cycling in New England rivers and salt marshes. On land, Si is taken up by vegetation, improving plant fitness and protecting plants from a variety of environmental stressors. In aquatic systems, diatoms, the dominant type of phytoplankton in coastal temperate waters, require Si to survive. My research demonstrates that LULC is an important driver of Si export to coastal systems, accounting for 40-70% of the variability of riverine fluxes. Developed watersheds export significantly (p=0.03) more Si than their forested counterparts, which I hypothesize is due to less vegetated cover, a known Si sink, in developed watersheds. Building on this, I calculated the amount of Si fixed by land plants globally (84 Tmol yr-1) and the percent (55%) of global terrestrial net primary production that can be attributed to active Si-accumulating organisms. Next, I created the first complete salt marsh Si budget by quantifying tidal creek fluxes and net Si accumulation in a relatively undisturbed low-nutrient salt marsh. Further, comparing this Si accumulation to that of a high-nutrient marsh revealed that the high-nutrient marsh contained significantly (p<0.05) more Si within the sediments, roots, and porewater. Combining my original data from six New England salt marshes with published values, I quantify the mode of Si accumulation (rejective, passive, or active) by Spartina grasses and the environmental controls on such accumulation. Finally, using radionuclide tracers 137Cs and 210Pb, I calculated vertical accretion rates of five salt marshes and compared these values to historical measurements. I found that accretion rates have slowed and this deceleration is driven, in part, by a decrease in organic matter accumulation. Together, this dissertation improves our knowledge of Si cycling in terrestrial and aquatic ecosystems, and identifies previously unrecognized ways in which humans are perturbing biogeochemical cycles at the land-sea interface

Does elevated CO2 alter silica uptake in trees?

© The Author(s), 2015. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Frontiers in Plant Science 5 (2015): 793, doi:10.3389/fpls.2014.00793.Human activities have greatly altered global carbon (C) and Nitrogen (N) cycling. In fact, atmospheric concentrations of carbon dioxide (CO2) have increased 40% over the last century and the amount of N cycling in the biosphere has more than doubled. In an effort to understand how plants will respond to continued global CO2 fertilization, long-term free-air CO2 enrichment experiments have been conducted at sites around the globe. Here we examine how atmospheric CO2 enrichment and N fertilization affects the uptake of silicon (Si) in the Duke Forest, North Carolina, a stand dominated by Pinus taeda (loblolly pine), and five hardwood species. Specifically, we measured foliar biogenic silica concentrations in five deciduous and one coniferous species across three treatments: CO2 enrichment, N enrichment, and N and CO2 enrichment. We found no consistent trends in foliar Si concentration under elevated CO2, N fertilization, or combined elevated CO2 and N fertilization. However, two-thirds of the tree species studied here have Si foliar concentrations greater than well-known Si accumulators, such as grasses. Based on net primary production values and aboveground Si concentrations in these trees, we calculated forest Si uptake rates under control and elevated CO2 concentrations. Due largely to increased primary production, elevated CO2 enhanced the magnitude of Si uptake between 20 and 26%, likely intensifying the terrestrial silica pump. This uptake of Si by forests has important implications for Si export from terrestrial systems, with the potential to impact C sequestration and higher trophic levels in downstream ecosystems.This research was supported in part by the Sloan Foundation in a fellowship to Robinson W. Fulweiler. The Duke Forest FACE was supported by his study was supported by the US Department of Energy (Grant No. DE-FG02-95ER62083) through the Office of Biological and Environmental Research (BER) and its National Institute for Global Environmental Change (NIGEC), Southeast Regional Center (SERC) at the University of Alabama, and by the US Forest Service through both the Southern Global Climate Change Program and the Southern Research Station. Adrien C. Finzi acknowledges ancillary support from the US NSF (DEB0236356)

Soil warming accelerates biogeochemical silica cycling in a temperate forest.

© The Author(s), 2019. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Gewirtzman, J., Tang, J., Melillo, J. M., Werner, W. J., Kurtz, A. C., Fulweiler, R. W., & Carey, J. C. Soil warming accelerates biogeochemical silica cycling in a temperate forest. Frontiers in Plant Science, 10, (2019): 1097, doi:10.3389/fpls.2019.01097.Biological cycling of silica plays an important role in terrestrial primary production. Soil warming stemming from climate change can alter the cycling of elements, such as carbon and nitrogen, in forested ecosystems. However, the effects of soil warming on the biogeochemical cycle of silica in forested ecosystems remain unexplored. Here we examine long-term forest silica cycling under ambient and warmed conditions over a 15-year period of experimental soil warming at Harvard Forest (Petersham, MA). Specifically, we measured silica concentrations in organic and mineral soils, and in the foliage and litter of two dominant species (Acer rubrum and Quercus rubra), in a large (30 × 30 m) heated plot and an adjacent control plot (30 × 30 m). In 2016, we also examined effects of heating on dissolved silica in the soil solution, and conducted a litter decomposition experiment using four tree species (Acer rubrum, Quercus rubra, Betula lenta, Tsuga canadensis) to examine effects of warming on the release of biogenic silica (BSi) from plants to soils. We find that tree foliage maintained constant silica concentrations in the control and warmed plots, which, coupled with productivity enhancements under warming, led to an increase in total plant silica uptake. We also find that warming drove an acceleration in the release of silica from decaying litter in three of the four species we examined, and a substantial increase in the silica dissolved in soil solution. However, we observe no changes in soil BSi stocks with warming. Together, our data indicate that warming increases the magnitude of silica uptake by vegetation and accelerates the internal cycling of silica in in temperate forests, with possible, and yet unresolved, effects on the delivery of silica from terrestrial to marine systems.This research was supported by the National Science Foundation (NSF PLR-1417763 to JT), the Geological Society of America (Stephen G. Pollock Undergraduate Research Grant to JG), the Institute at Brown for Environment and Society, and the Marine Biological Laboratory. Sample analysis and Fulweiler’s involvement were supported by Boston University and a Bullard Fellowship from Harvard University. The soil warming experiment was supported by the National Science Foundation (DEB-0620443) and Department of Energy (DE-FC02-06-ER641577 and DE-SC0005421)

Evaluation of simulated soil carbon dynamics in Arctic-Boreal ecosystems

© The Author(s), 2020. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Huntzinger, D. N., Schaefer, K., Schwalm, C., Fisher, J. B., Hayes, D., Stofferahn, E., Carey, J., Michalak, A. M., Wei, Y., Jain, A. K., Kolus, H., Mao, J., Poulter, B., Shi, X., Tang, J., & Tian, H. Evaluation of simulated soil carbon dynamics in Arctic-Boreal ecosystems. Environmental Research Letters, 15(2), (2020): 025005, doi:10.1088/1748-9326/ab6784.Given the magnitude of soil carbon stocks in northern ecosystems, and the vulnerability of these stocks to climate warming, land surface models must accurately represent soil carbon dynamics in these regions. We evaluate soil carbon stocks and turnover rates, and the relationship between soil carbon loss with soil temperature and moisture, from an ensemble of eleven global land surface models. We focus on the region of NASA's Arctic-Boreal vulnerability experiment (ABoVE) in North America to inform data collection and model development efforts. Models exhibit an order of magnitude difference in estimates of current total soil carbon stocks, generally under- or overestimating the size of current soil carbon stocks by greater than 50 PgC. We find that a model's soil carbon stock at steady-state in 1901 is the prime driver of its soil carbon stock a hundred years later—overwhelming the effect of environmental forcing factors like climate. The greatest divergence between modeled and observed soil carbon stocks is in regions dominated by peat and permafrost soils, suggesting that models are failing to capture the frozen soil carbon dynamics of permafrost regions. Using a set of functional benchmarks to test the simulated relationship of soil respiration to both soil temperature and moisture, we find that although models capture the observed shape of the soil moisture response of respiration, almost half of the models examined show temperature sensitivities, or Q10 values, that are half of observed. Significantly, models that perform better against observational constraints of respiration or carbon stock size do not necessarily perform well in terms of their functional response to key climatic factors like changing temperature. This suggests that models may be arriving at the right result, but for the wrong reason. The results of this work can help to bridge the gap between data and models by both pointing to the need to constrain initial carbon pool sizes, as well as highlighting the importance of incorporating functional benchmarks into ongoing, mechanistic modeling activities such as those included in ABoVE.This work was supported by NASA'S Arctic Boreal Vulnerability Experiment (ABoVE; https://above.nasa.gov); NNN13D504T. Funding for the Multi-scale synthesis and Terrestrial Model Intercomparison Project (MsTMIP; https://nacp.ornl.gov/MsTMIP.shtml) activity was provided through NASA ROSES Grant #NNX10AG01A. Data management support for preparing, documenting, and distributing model driver and output data was performed by the Modeling and Synthesis Thematic Data Center at Oak Ridge National Laboratory (MAST-DC; https://nacp.ornl.gov), with funding through NASA ROSES Grant #NNH10AN681. Finalized MsTMIP data products are archived at the ORNL DAAC (https://daac.ornl.gov). We also acknowledge the modeling groups that provided results to MsTMIP. The synthesis of site-level soil respiration, temperature, and moisture data reported in Carey et al 2016a, 2016b) was funded by the US Geological Survey (USGS) John Wesley Powell Center for Analysis and Synthesis Award G13AC00193. Additional support for that work was also provided by the USGS Land Carbon Program. JBF carried out the research at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration. California Institute of Technology. Government sponsorship acknowledged

Directly Measured Denitrification Reveals Oyster Aquaculture and Restored Oyster Reefs Remove Nitrogen at Comparable High Rates

Coastal systems are increasingly impacted by over-enrichment of nutrients, which has cascading effects for ecosystem functioning. Oyster restoration and aquaculture are both hypothesized to mitigate excessive nitrogen (N) loads via benthic denitrification. The degree to which these management activities perform similar functions for removing N, however, has not been extensively examined in New England, a place where nutrient runoff is high and increasing oyster (Crassostrea virginica) restoration and aquaculture activity is taking place. Here, we use a novel in situ methodology to directly measure net N2 and O2 fluxes across the sediment-water interface in a shallow (~1 m) coastal pond in southern Rhode Island. We collected data seasonally during 2013 and 2014 at restored oyster reefs, oyster aquaculture, oyster cultch (shell), and bare sediment. Restored oyster reefs and aquaculture had the highest mean (±SE) denitrification rates, 581.9 (±164.2) and 346 (±168.6) μmol N2−N m−2 h−1, respectively, and are among the highest recorded for oyster-dominated environments. Denitrification rates at sites with oyster cultch were 60.9 (±44.3) μmol N2−N m−2 h−1, which is substantially less than the sites with active oysters but still more than 50% higher than denitrification rates measured in bare sediment (24.4 ± 10.1 μmol N2–N m−2 h−1). The increase in denitrification rates at treatments, however, varied by season and the greatest rates for restored reefs were in the fall. Overall, the greatest aggregate denitrification rates occurred in the fall. Sediment oxygen demand (SOD) followed similar patterns but with greater overall rates in the summer, and displayed a strong linear relationship with denitrification (R2 = 0.93). Our results demonstrate that habitats associated with live oysters have higher net denitrification rates and that oyster reef restoration and oyster aquaculture may provide similar benefits to the ecosystem in terms of N removal. However, gas fluxes may also be affected where three-dimensional structure is introduced via oyster shell cultch and this appears to be seasonally-dependent. These data will be important for managers as they incorporate oysters into nutrient reduction strategies and consider system-level trade-offs in services provided by oyster reef restoration and aquaculture activities

Directly Measured Denitrification Reveals Oyster Aquaculture and Restored Oyster Reefs Remove Nitrogen at Comparable High Rates

Coastal systems are increasingly impacted by over-enrichment of nutrients, which has cascading effects for ecosystem functioning. Oyster restoration and aquaculture are both hypothesized to mitigate excessive nitrogen (N) loads via benthic denitrification. The degree to which these management activities perform similar functions for removing N, however, has not been extensively examined in New England, a place where nutrient runoff is high and increasing oyster (Crassostrea virginica) restoration and aquaculture activity is taking place. Here, we use a novel in situ methodology to directly measure net N2 and O2 fluxes across the sediment-water interface in a shallow (~1 m) coastal pond in southern Rhode Island. We collected data seasonally during 2013 and 2014 at restored oyster reefs, oyster aquaculture, oyster cultch (shell), and bare sediment. Restored oyster reefs and aquaculture had the highest mean (±SE) denitrification rates, 581.9 (±164.2) and 346 (±168.6) μmol N2−N m−2 h−1, respectively, and are among the highest recorded for oyster-dominated environments. Denitrification rates at sites with oyster cultch were 60.9 (±44.3) μmol N2−N m−2 h−1, which is substantially less than the sites with active oysters but still more than 50% higher than denitrification rates measured in bare sediment (24.4 ± 10.1 μmol N2–N m−2 h−1). The increase in denitrification rates at treatments, however, varied by season and the greatest rates for restored reefs were in the fall. Overall, the greatest aggregate denitrification rates occurred in the fall. Sediment oxygen demand (SOD) followed similar patterns but with greater overall rates in the summer, and displayed a strong linear relationship with denitrification (R2 = 0.93). Our results demonstrate that habitats associated with live oysters have higher net denitrification rates and that oyster reef restoration and oyster aquaculture may provide similar benefits to the ecosystem in terms of N removal. However, gas fluxes may also be affected where three-dimensional structure is introduced via oyster shell cultch and this appears to be seasonally-dependent. These data will be important for managers as they incorporate oysters into nutrient reduction strategies and consider system-level trade-offs in services provided by oyster reef restoration and aquaculture activities

Pan-cancer Alterations of the MYC Oncogene and Its Proximal Network across the Cancer Genome Atlas

Although theMYConcogene has been implicated incancer, a systematic assessment of alterations ofMYC, related transcription factors, and co-regulatoryproteins, forming the proximal MYC network (PMN),across human cancers is lacking. Using computa-tional approaches, we define genomic and proteo-mic features associated with MYC and the PMNacross the 33 cancers of The Cancer Genome Atlas.Pan-cancer, 28% of all samples had at least one ofthe MYC paralogs amplified. In contrast, the MYCantagonists MGA and MNT were the most frequentlymutated or deleted members, proposing a roleas tumor suppressors.MYCalterations were mutu-ally exclusive withPIK3CA,PTEN,APC,orBRAFalterations, suggesting that MYC is a distinct onco-genic driver. Expression analysis revealed MYC-associated pathways in tumor subtypes, such asimmune response and growth factor signaling; chro-matin, translation, and DNA replication/repair wereconserved pan-cancer. This analysis reveals insightsinto MYC biology and is a reference for biomarkersand therapeutics for cancers with alterations ofMYC or the PMN

Pan-Cancer Analysis of lncRNA Regulation Supports Their Targeting of Cancer Genes in Each Tumor Context

Long noncoding RNAs (lncRNAs) are commonly dys-regulated in tumors, but only a handful are known toplay pathophysiological roles in cancer. We inferredlncRNAs that dysregulate cancer pathways, onco-genes, and tumor suppressors (cancer genes) bymodeling their effects on the activity of transcriptionfactors, RNA-binding proteins, and microRNAs in5,185 TCGA tumors and 1,019 ENCODE assays.Our predictions included hundreds of candidateonco- and tumor-suppressor lncRNAs (cancerlncRNAs) whose somatic alterations account for thedysregulation of dozens of cancer genes and path-ways in each of 14 tumor contexts. To demonstrateproof of concept, we showed that perturbations tar-geting OIP5-AS1 (an inferred tumor suppressor) andTUG1 and WT1-AS (inferred onco-lncRNAs) dysre-gulated cancer genes and altered proliferation ofbreast and gynecologic cancer cells. Our analysis in-dicates that, although most lncRNAs are dysregu-lated in a tumor-specific manner, some, includingOIP5-AS1, TUG1, NEAT1, MEG3, and TSIX, synergis-tically dysregulate cancer pathways in multiple tumorcontexts

Genomic, Pathway Network, and Immunologic Features Distinguishing Squamous Carcinomas

This integrated, multiplatform PanCancer Atlas study co-mapped and identified distinguishing molecular features of squamous cell carcinomas (SCCs) from five sites associated with smokin

Spatial Organization and Molecular Correlation of Tumor-Infiltrating Lymphocytes Using Deep Learning on Pathology Images

Beyond sample curation and basic pathologic characterization, the digitized H&E-stained images of TCGA samples remain underutilized. To highlight this resource, we present mappings of tumorinfiltrating lymphocytes (TILs) based on H&E images from 13 TCGA tumor types. These TIL maps are derived through computational staining using a convolutional neural network trained to classify patches of images. Affinity propagation revealed local spatial structure in TIL patterns and correlation with overall survival. TIL map structural patterns were grouped using standard histopathological parameters. These patterns are enriched in particular T cell subpopulations derived from molecular measures. TIL densities and spatial structure were differentially enriched among tumor types, immune subtypes, and tumor molecular subtypes, implying that spatial infiltrate state could reflect particular tumor cell aberration states. Obtaining spatial lymphocytic patterns linked to the rich genomic characterization of TCGA samples demonstrates one use for the TCGA image archives with insights into the tumor-immune microenvironment