Search

This paper provides an overview of the United States (US) Department of Energy's (DOE's) Energy Exascale Earth System Model version 2 (E3SMv2) fully coupled regionally refined model (RRM) and documents the overall atmosphere, land, and river results from the Coupled Model Intercomparison Project 6 (CMIP6) DECK (Diagnosis, Evaluation, and Characterization of Klima) and historical simulations - a first-of-its-kind set of climate production simulations using RRM. The North American (NA) RRM (NARRM) is developed as the high-resolution configuration of E3SMv2 with the primary goal of more explicitly addressing DOE's mission needs regarding impacts to the US energy sector facing Earth system changes. The NARRM features finer horizontal resolution grids centered over NA, consisting of 25→100¯km atmosphere and land, a 0.125 river-routing model, and 14→60¯km ocean and sea ice. By design, the computational cost of NARRM is 1/43× of the uniform low-resolution (LR) model at 100¯km but only 1/4¯10¯%-20¯% of a globally uniform high-resolution model at 25¯km. A novel hybrid time step strategy for the atmosphere is key for NARRM to achieve improved climate simulation fidelity within the high-resolution patch without sacrificing the overall global performance. The global climate, including climatology, time series, sensitivity, and feedback, is confirmed to be largely identical between NARRM and LR as quantified with typical climate metrics. Over the refined NA area, NARRM is generally superior to LR, including for precipitation and clouds over the contiguous US (CONUS), summertime marine stratocumulus clouds off the coast of California, liquid and ice phase clouds near the North Pole region, extratropical cyclones, and spatial variability in land hydrological processes. The improvements over land are related to the better-resolved topography in NARRM, whereas those over ocean are attributable to the improved air-sea interactions with finer grids for both atmosphere and ocean and sea ice. Some features appear insensitive to the resolution change analyzed here, for instance the diurnal propagation of organized mesoscale convective systems over CONUS and the warm-season land-atmosphere coupling at the southern Great Plains. In summary, our study presents a realistically efficient approach to leverage the fully coupled RRM framework for a standard Earth system model release and high-resolution climate production simulations.

Low-grade, sustained inflammation in white adipose tissue (WAT) characterizes obesity and coincides with type 2 diabetes mellitus (T2DM). However, pharmacological targeting of inflammation lacks durable therapeutic effects in insulin-resistant conditions. Through a computational screen, we discovered that the FDA-approved rheumatoid arthritis drug auranofin improved insulin sensitivity and normalized obesity-associated abnormalities, including hepatic steatosis and hyperinsulinemia in mouse models of T2DM. We also discovered that auranofin accumulation in WAT depleted inflammatory responses to a high-fat diet without altering body composition in obese wild-type mice. Surprisingly, elevated leptin levels and blunted beta-adrenergic receptor activity achieved by leptin receptor deletion abolished the antidiabetic effects of auranofin. These experiments also revealed that the metabolic benefits of leptin reduction were superior to immune impacts of auranofin in WAT. Our studies uncover important metabolic properties of anti-inflammatory treatments and contribute to the notion that leptin reduction in the periphery can be accomplished to treat obesity and T2DM.

This work documents version two of the Department of Energy's Energy Exascale Earth System Model (E3SM). E3SMv2 is a significant evolution from its predecessor E3SMv1, resulting in a model that is nearly twice as fast and with a simulated climate that is improved in many metrics. We describe the physical climate model in its lower horizontal resolution configuration consisting of 110 km atmosphere, 165 km land, 0.5° river routing model, and an ocean and sea ice with mesh spacing varying between 60 km in the mid-latitudes and 30 km at the equator and poles. The model performance is evaluated with Coupled Model Intercomparison Project Phase 6 Diagnosis, Evaluation, and Characterization of Klima simulations augmented with historical simulations as well as simulations to evaluate impacts of different forcing agents. The simulated climate has many realistic features of the climate system, with notable improvements in clouds and precipitation compared to E3SMv1. E3SMv1 suffered from an excessively high equilibrium climate sensitivity (ECS) of 5.3 K. In E3SMv2, ECS is reduced to 4.0 K which is now within the plausible range based on a recent World Climate Research Program assessment. However, a number of important biases remain including a weak Atlantic Meridional Overturning Circulation, deficiencies in the characteristics and spectral distribution of tropical atmospheric variability, and a significant underestimation of the observed warming in the second half of the historical period. An analysis of single-forcing simulations indicates that correcting the historical temperature bias would require a substantial reduction in the magnitude of the aerosol-related forcing.

Although macroautophagy/autophagy deficiency causes degenerative diseases, the deletion of essential autophagy genes in adipocytes paradoxically reduces body weight. Brown adipose tissue (BAT) plays an important role in body weight regulation and metabolic control. However, the key cellular mechanisms that maintain BAT function remain poorly understood. in this study, we showed that global or brown adipocyte-specific deletion of pink1, a Parkinson disease-related gene involved in selective mitochondrial autophagy (mitophagy), induced BAT dysfunction, and obesity-prone type in mice. Defective mitochondrial function is among the upstream signals that activate the NLRP3 inflammasome. NLRP3 was induced in brown adipocyte precursors (BAPs) from pink1 knockout (KO) mice. Unexpectedly, NLRP3 induction did not induce canonical inflammasome activity. Instead, NLRP3 induction led to the differentiation of pink1 KO BAPs into white-like adipocytes by increasing the expression of white adipocyte-specific genes and repressing the expression of brown adipocyte-specific genes. nlrp3 deletion in pink1 knockout mice reversed BAT dysfunction. Conversely, adipose tissue-specific atg7 KO mice showed significantly lower expression of Nlrp3 in their BAT. Overall, our data suggest that the role of mitophagy is different from general autophagy in regulating adipose tissue and whole-body energy metabolism. Our results uncovered a new mitochondria-NLRP3 pathway that induces BAT dysfunction. The ability of the nlrp3 knockouts to rescue BAT dysfunction suggests the transcriptional function of NLRP3 as an unexpected, but a quite specific therapeutic target for obesity-related metabolic diseases.Abbreviations: ACTB: actin, beta; BAPs: brown adipocyte precursors; BAT: brown adipose tissue; BMDMs: bone marrow-derived macrophages; CASP1: caspase 1; CEBPA: CCAAT/enhancer binding protein (C/EBP), alpha; ChIP: chromatin immunoprecipitation; EE: energy expenditure; HFD: high-fat diet; IL1B: interleukin 1 beta; ITT: insulin tolerance test; KO: knockout; LPS: lipopolysaccharide; NLRP3: NLR family, pyrin domain containing 3; PINK1: PTEN induced putative kinase 1; PRKN: parkin RBR E3 ubiquitin protein ligase; RD: regular diet; ROS: reactive oxygen species; RT: room temperature; UCP1: uncoupling protein 1 (mitochondrial, proton carrier); WT: wild-type.

This study provides an overview of the coupled high-resolution Version 1 of the Energy Exascale Earth System Model (E3SMv1) and documents the characteristics of a 50-year-long high-resolution control simulation with time-invariant 1950 forcings following the HighResMIP protocol. In terms of global root-mean-squared error metrics, this high-resolution simulation is generally superior to results from the low-resolution configuration of E3SMv1 (due to resolution, tuning changes, and possibly initialization procedure) and compares favorably to models in the CMIP5 ensemble. Ocean and sea ice simulation is particularly improved, due to better resolution of bathymetry, the ability to capture more variability and extremes in winds and currents, and the ability to resolve mesoscale ocean eddies. The largest improvement in this regard is an ice-free Labrador Sea, which is a major problem at low resolution. Interestingly, several features found to improve with resolution in previous studies are insensitive to resolution or even degrade in E3SMv1. Most notable in this regard are warm bias and associated stratocumulus deficiency in eastern subtropical oceans and lack of improvement in El Niño. Another major finding of this study is that resolution increase had negligible impact on climate sensitivity (measured by net feedback determined through uniform +4K prescribed sea surface temperature increase) and aerosol sensitivity. Cloud response to resolution increase consisted of very minor decrease at all levels. Large-scale patterns of precipitation bias were also relatively unaffected by grid spacing.

This work documents the first version of the U.S. Department of Energy (DOE) new Energy Exascale Earth System Model (E3SMv1). We focus on the standard resolution of the fully coupled physical model designed to address DOE mission-relevant water cycle questions. Its components include atmosphere and land (110-km grid spacing), ocean and sea ice (60 km in the midlatitudes and 30 km at the equator and poles), and river transport (55 km) models. This base configuration will also serve as a foundation for additional configurations exploring higher horizontal resolution as well as augmented capabilities in the form of biogeochemistry and cryosphere configurations. The performance of E3SMv1 is evaluated by means of a standard set of Coupled Model Intercomparison Project Phase 6 (CMIP6) Diagnosis, Evaluation, and Characterization of Klima simulations consisting of a long preindustrial control, historical simulations (ensembles of fully coupled and prescribed SSTs) as well as idealized CO2 forcing simulations. The model performs well overall with biases typical of other CMIP-class models, although the simulated Atlantic Meridional Overturning Circulation is weaker than many CMIP-class models. While the E3SMv1 historical ensemble captures the bulk of the observed warming between preindustrial (1850) and present day, the trajectory of the warming diverges from observations in the second half of the twentieth century with a period of delayed warming followed by an excessive warming trend. Using a two-layer energy balance model, we attribute this divergence to the model's strong aerosol-related effective radiative forcing (ERFari+aci = −1.65 W/m2) and high equilibrium climate sensitivity (ECS = 5.3 K).

Specific metabolic underpinnings of androgen receptor (AR)-driven growth in prostate adenocarcinoma (PCa) are largely undefined, hindering the development of strategies to leverage the metabolic dependencies of this disease when hormonal manipulations fail. Here we show that the mitochondrial pyruvate carrier (MPC), a critical metabolic conduit linking cytosolic and mitochondrial metabolism, is transcriptionally regulated by AR. Experimental MPC inhibition restricts proliferation and metabolic outputs of the citric acid cycle (TCA) including lipogenesis and oxidative phosphorylation in AR-driven PCa models. Mechanistically, metabolic disruption resulting from MPC inhibition activates the eIF2α/ATF4 integrated stress response (ISR). ISR signaling prevents cell cycle progression while coordinating salvage efforts, chiefly enhanced glutamine assimilation into the TCA, to regain metabolic homeostasis. We confirm that MPC function is operant in PCa tumors in-vivo using isotopomeric metabolic flux analysis. In turn, we apply a clinically viable small molecule targeting the MPC, MSDC0160, to pre-clinical PCa models and find that MPC inhibition suppresses tumor growth in hormone-responsive and castrate-resistant conditions. Collectively, our findings characterize the MPC as a tractable therapeutic target in AR-driven prostate tumors.

Chronic inflammation accompanies obesity and limits subcutaneous white adipose tissue (WAT) expandability, accelerating the development of insulin resistance and type 2 diabetes mellitus. MicroRNAs (miRNAs) influence expression of many metabolic genes in fat cells, but physiological roles in WAT remain poorly characterized. Here, we report that expression of the miRNA miR-30a in subcutaneous WAT corresponds with insulin sensitivity in obese mice and humans. To examine the hypothesis that restoration of miR-30a expression in WAT improves insulin sensitivity, we injected adenovirus (Adv) expressing miR-30a into the subcutaneous fat pad of diabetic mice. Exogenous miR-30a expression in the subcutaneous WAT depot of obese mice coupled improved insulin sensitivity and increased energy expenditure with decreased ectopic fat deposition in the liver and reduced WAT inflammation. High-throughput proteomic profiling and RNA-Seq suggested that miR-30a targets the transcription factor STAT1 to limit the actions of the proinflammatory cytokine interferon-γ (IFN-γ) that would otherwise restrict WAT expansion and decrease insulin sensitivity. We further demonstrated that miR-30a opposes the actions of IFN-γ, suggesting an important role for miR-30a in defending adipocytes against proinflammatory cytokines that reduce peripheral insulin sensitivity. Together, our data identify a critical molecular signaling axis, elements of which are involved in uncoupling obesity from metabolic dysfunction.

The precise molecular alterations driving castration-resistant prostate cancer (CRPC) are not clearly understood. Using a novel network-based integrative approach, here, we show distinct alterations in the hexosamine biosynthetic pathway (HBP) to be critical for CRPC. Expression of HBP enzyme glucosamine-phosphate N-acetyltransferase 1 (GNPNAT1) is found to be significantly decreased in CRPC compared with localized prostate cancer (PCa). Genetic loss-of-function of GNPNAT1 in CRPC-like cells increases proliferation and aggressiveness, in vitro and in vivo. This is mediated by either activation of the PI3K-AKT pathway in cells expressing full-length androgen receptor (AR) or by specific protein 1 (SP1)-regulated expression of carbohydrate response element-binding protein (ChREBP) in cells containing AR-V7 variant. Strikingly, addition of the HBP metabolite UDP-N-acetylglucosamine (UDP-GlcNAc) to CRPC-like cells significantly decreases cell proliferation, both in-vitro and in animal studies, while also demonstrates additive efficacy when combined with enzalutamide in-vitro. These observations demonstrate the therapeutic value of targeting HBP in CRPC.

The GraphBLAS standard (GraphBlas.org) is being developed to bring the potential of matrix-based graph algorithms to the broadest possible audience. Mathematically, the GraphBLAS defines a core set of matrix-based graph operations that can be used to implement a wide class of graph algorithms in a wide range of programming environments. This paper provides an introduction to the mathematics of the GraphBLAS. Graphs represent connections between vertices with edges. Matrices can represent a wide range of graphs using adjacency matrices or incidence matrices. Adjacency matrices are often easier to analyze while incidence matrices are often better for representing data. Fortunately, the two are easily connected by matrix multiplication. A key feature of matrix mathematics is that a very small number of matrix operations can be used to manipulate a very wide range of graphs. This composability of a small number of operations is the foundation of the GraphBLAS. A standard such as the GraphBLAS can only be effective if it has low performance overhead. Performance measurements of prototype GraphBLAS implementations indicate that the overhead is low.

BackgroundMutations in microtubule-regulating genes are associated with disorders of neuronal migration and microcephaly. Regulation of centriole length has been shown to underlie the pathogenesis of certain ciliopathy phenotypes. Using a next-generation sequencing approach, we identified mutations in a novel centriolar disease gene in a kindred with an embryonic lethal ciliopathy phenotype and in a patient with primary microcephaly.Methods and resultsWhole exome sequencing data from a non-consanguineous Caucasian kindred exhibiting mid-gestation lethality and ciliopathic malformations revealed two novel non-synonymous variants in CENPF, a microtubule-regulating gene. All four affected fetuses showed segregation for two mutated alleles [IVS5-2A>C, predicted to abolish the consensus splice-acceptor site from exon 6; c.1744G>T, p.E582X]. In a second unrelated patient exhibiting microcephaly, we identified two CENPF mutations [c.1744G>T, p.E582X; c.8692 C>T, p.R2898X] by whole exome sequencing. We found that CENP-F colocalised with Ninein at the subdistal appendages of the mother centriole in mouse inner medullary collecting duct cells. Intraflagellar transport protein-88 (IFT-88) colocalised with CENP-F along the ciliary axonemes of renal epithelial cells in age-matched control human fetuses but did not in truncated cilia of mutant CENPF kidneys. Pairwise co-immunoprecipitation assays of mitotic and serum-starved HEKT293 cells confirmed that IFT88 precipitates with endogenous CENP-F.ConclusionsOur data identify CENPF as a new centriolar disease gene implicated in severe human ciliopathy and microcephaly related phenotypes. CENP-F has a novel putative function in ciliogenesis and cortical neurogenesis.

The analysis of graphs has become increasingly important to a wide range of applications. Graph analysis presents a number of unique challenges in the areas of (1) software complexity, (2) data complexity, (3) security, (4) mathematical complexity, (5) theoretical analysis, (6) serial performance, and (7) parallel performance. Implementing graph algorithms using matrix-based approaches provides a number of promising solutions to these challenges. The GraphBLAS standard (istcbigdata.org/GraphBlas) is being developed to bring the potential of matrix based graph algorithms to the broadest possible audience. The GraphBLAS mathematically defines a core set of matrix-based graph operations that can be used to implement a wide class of graph algorithms in a wide range of programming environments. This paper provides an introduction to the GraphBLAS and describes how the GraphBLAS can be used to address many of the challenges associated with analysis of graphs.

This paper presents results from a 40 year Weather Research and Forecasting (WRF) based dynamical downscaling experiment performed at 12 km horizontal grid spacing, centered on the state of California, and forced by a 1° × 1.25° finite-volume current-climate Community Climate System Model ver. 3 (CCSM3) simulation. In-depth comparisons between modeled and observed regional-average precipitation, 2 m temperature, and snowpack are performed. The regional model reproduces the spatial distribution of precipitation quite well, but substantially overestimates rainfall along windward slopes. This is due to strong overprediction of precipitation intensity; precipitation frequency is actually underpredicted by the model. Moisture fluxes impinging on the coast seem to be well-represented over California, implying that precipitation bias is caused by processes internal to WRF. Positive-definite moisture advection and use of the Grell cumulus parameterization result in some decrease in precipitation bias, but other sources are needed to explain the full bias magnitude. Surface temperature is well simulated in all seasons except summer, when overly-dry soil moisture results in a several degree warm bias in both CCSM3 and WRF. Additionally, coastal temperatures appear to be too warm due to a coastal sea surface temperature bias inherited from CCSM3. Modeled snowfall/snowmelt agrees quite well with observations, but snow water equivalent is found to be much too low due to monthly reinitialization of all regional model fields from CCSM3 values.