Open-ST: High-resolution spatial transcriptomics in 3D.

Spatial transcriptomics (ST) methods unlock molecular mechanisms underlying tissue development, homeostasis, or disease. However, there is a need for easy-to-use, high-resolution, cost-efficient, and 3D-scalable methods. Here, we report Open-ST, a sequencing-based, open-source experimental and computational resource to address these challenges and to study the molecular organization of tissues in 2D and 3D. In mouse brain, Open-ST captured transcripts at subcellular resolution and reconstructed cell types. In primary head-and-neck tumors and patient-matched healthy/metastatic lymph nodes, Open-ST captured the diversity of immune, stromal, and tumor populations in space, validated by imaging-based ST. Distinct cell states were organized around cell-cell communication hotspots in the tumor but not the metastasis. Strikingly, the 3D reconstruction and multimodal analysis of the metastatic lymph node revealed spatially contiguous structures not visible in 2D and potential biomarkers precisely at the 3D tumor/lymph node boundary. All protocols and software are available at https://rajewsky-lab.github.io/openst.

The underappreciated diversity of bile acid modifications.

The repertoire of modifications to bile acids and related steroidal lipids by host and microbial metabolism remains incompletely characterized. To address this knowledge gap, we created a reusable resource of tandem mass spectrometry (MS/MS) spectra by filtering 1.2 billion publicly available MS/MS spectra for bile-acid-selective ion patterns. Thousands of modifications are distributed throughout animal and human bodies as well as microbial cultures. We employed this MS/MS library to identify polyamine bile amidates, prevalent in carnivores. They are present in humans, and their levels alter with a diet change from a Mediterranean to a typical American diet. This work highlights the existence of many more bile acid modifications than previously recognized and the value of leveraging public large-scale untargeted metabolomics data to discover metabolites. The availability of a modification-centric bile acid MS/MS library will inform future studies investigating bile acid roles in health and disease.

An Integrated TCGA Pan-Cancer Clinical Data Resource to Drive High-Quality Survival Outcome Analytics.

For a decade, The Cancer Genome Atlas (TCGA) program collected clinicopathologic annotation data along with multi-platform molecular profiles of more than 11,000 human tumors across 33 different cancer types. TCGA clinical data contain key features representing the democratized nature of the data collection process. To ensure proper use of this large clinical dataset associated with genomic features, we developed a standardized dataset named the TCGA Pan-Cancer Clinical Data Resource (TCGA-CDR), which includes four major clinical outcome endpoints. In addition to detailing major challenges and statistical limitations encountered during the effort of integrating the acquired clinical data, we present a summary that includes endpoint usage recommendations for each cancer type. These TCGA-CDR findings appear to be consistent with cancer genomics studies independent of the TCGA effort and provide opportunities for investigating cancer biology using clinical correlates at an unprecedented scale.

Grow now, pay later: When should a bacterium go into debt?

Microbes grow in a wide variety of environments and must balance growth and stress resistance. Despite the prevalence of such trade-offs, understanding of their role in nonsteady environments is limited. In this study, we introduce a mathematical model of "growth debt," where microbes grow rapidly initially, paying later with slower growth or heightened mortality. We first compare our model to a classical chemostat experiment, validating our proposed dynamics and quantifying Escherichia coli's stress resistance dynamics. Extending the chemostat theory to include serial-dilution cultures, we derive phase diagrams for the persistence of "debtor" microbes. We find that debtors cannot coexist with nondebtors if "payment" is increased mortality but can coexist if it lowers enzyme affinity. Surprisingly, weak noise considerably extends the persistence of resistance elements, pertinent for antibiotic resistance management. Our microbial debt theory, broadly applicable across many environments, bridges the gap between chemostat and serial dilution systems.

Costs of being a diet generalist for the protist predator Dictyostelium discoideum.

Consumers range from specialists that feed on few resources to generalists that feed on many. Generalism has the clear advantage of having more resources to exploit, but the costs that limit generalism are less clear. We explore two understudied costs of generalism in a generalist amoeba predator, Dictyostelium discoideum, feeding on naturally co-occurring bacterial prey. Both involve costs of combining prey that are suitable on their own. First, amoebas exhibit a reduction in growth rate when they switched to one species of prey bacteria from another compared to controls that experience only the second prey. The effect was consistent across all six tested species of bacteria. These switching costs typically disappear within a day, indicating adjustment to new prey bacteria. This suggests that these costs are physiological. Second, amoebas usually grow more slowly on mixtures of prey bacteria compared to the expectation based on their growth on single prey. There were clear mixing costs in three of the six tested prey mixtures, and none showed significant mixing benefits. These results support the idea that, although amoebas can consume a variety of prey, they must use partially different methods and thus must pay costs to handle multiple prey, either sequentially or simultaneously.

Monitoring resilience in bursts.

The possibility to anticipate critical transitions through detecting loss of resilience has attracted attention in many fields. Resilience indicators rely on the mathematical concept of critical slowing down, which means that a system recovers more slowly from external perturbations when it gets closer to tipping point. This decrease in recovery rate can be reflected in rising autocorrelation and variance in data. To test whether resilience is changing, resilience indicators are often calculated using a moving window in long, continuous time series of the system. However, for some systems, it may be more feasible to collect several high-resolution time series in short periods of time, i.e., in bursts. Resilience indicators can then be calculated to detect a change of resilience between such bursts. Here, we compare the performance of both methods using simulated data and showcase the possible use of bursts in a case study using mood data to anticipate depression in a patient. With the same number of data points, the burst approach outperformed the moving window method, suggesting that it is possible to downsample the continuous time series and still signal an upcoming transition. We suggest guidelines to design an optimal sampling strategy. Our results imply that using bursts of data instead of continuous time series may improve the capacity to detect changes in resilience. This method is promising for a variety of fields, such as human health, epidemiology, or ecology, where continuous monitoring can be costly or unfeasible.

Spatial synchrony cascades across ecosystem boundaries and up food webs via resource subsidies.

Cross-ecosystem subsidies are critical to ecosystem structure and function, especially in recipient ecosystems where they are the primary source of organic matter to the food web. Subsidies are indicative of processes connecting ecosystems and can couple ecological dynamics across system boundaries. However, the degree to which such flows can induce cross-ecosystem cascades of spatial synchrony, the tendency for system fluctuations to be correlated across locations, is not well understood. Synchrony has destabilizing effects on ecosystems, adding to the importance of understanding spatiotemporal patterns of synchrony transmission. In order to understand whether and how spatial synchrony cascades across the marine-terrestrial boundary via resource subsidies, we studied the relationship between giant kelp forests on rocky nearshore reefs and sandy beach ecosystems that receive resource subsidies in the form of kelp wrack (detritus). We found that synchrony cascades from rocky reefs to sandy beaches, with spatiotemporal patterns mediated by fluctuations in live kelp biomass, wave action, and beach width. Moreover, wrack deposition synchronized local abundances of shorebirds that move among beaches seeking to forage on wrack-associated invertebrates, demonstrating that synchrony due to subsidies propagates across trophic levels in the recipient ecosystem. Synchronizing resource subsidies likely play an underappreciated role in the spatiotemporal structure, functioning, and stability of ecosystems.

Toward a metabolic theory of catchments: Scaling of water and carbon fluxes with size.

Allometric scaling relations are widely used to link biological processes to body size in nature. Several studies have shown that such scaling laws hold also for natural ecosystems, including individual trees and forests, riverine metabolism, and river network organization. However, the derivation of scaling laws for catchment-scale water and carbon fluxes has not been achieved so far. Here, we focus on scaling relations of catchment green metabolism, defined as the set of ecohydrological and biogeochemical processes through which vegetation assemblages in catchments maintain their structure and react to the surrounding environment. By revising existing plant size-density relationships and integrating them across large-scale domains, we show that the ecohydrological fluxes occurring at the catchment scale are invariant with respect to the above-ground vegetation biomass per unit area of the basin, while they scale linearly with catchment size. We thus demonstrate that the sublinear scaling of plant metabolism results in an isometric scaling at catchment and regional scales. Deviations from such predictions are further shown to collapse onto a common distribution, thus incorporating natural fluctuations due to resource limitations into a generalized scaling theory. Results from scaling arguments are supported by hyperresolution ecohydrological simulations and remote sensing observations.

Using global remote camera data of a solitary species complex to evaluate the drivers of group formation.

The social system of animals involves a complex interplay between physiology, natural history, and the environment. Long relied upon discrete categorizations of "social" and "solitary" inhibit our capacity to understand species and their interactions with the world around them. Here, we use a globally distributed camera trapping dataset to test the drivers of aggregating into groups in a species complex (martens and relatives, family Mustelidae, Order Carnivora) assumed to be obligately solitary. We use a simple quantification, the probability of being detected in a group, that was applied across our globally derived camera trap dataset. Using a series of binomial generalized mixed-effects models applied to a dataset of 16,483 independent detections across 17 countries on four continents we test explicit hypotheses about potential drivers of group formation. We observe a wide range of probabilities of being detected in groups within the solitary model system, with the probability of aggregating in groups varying by more than an order of magnitude. We demonstrate that a species' context-dependent proclivity toward aggregating in groups is underpinned by a range of resource-related factors, primarily the distribution of resources, with increasing patchiness of resources facilitating group formation, as well as interactions between environmental conditions (resource constancy/winter severity) and physiology (energy storage capabilities). The wide variation in propensities to aggregate with conspecifics observed here highlights how continued failure to recognize complexities in the social behaviors of apparently solitary species limits our understanding not only of the individual species but also the causes and consequences of group formation.

Effective green treatment of sewage sludge from Fenton reactions: Utilizing MoS2 for sustainable resource recovery.

Effectively managing sewage sludge from Fenton reactions in an eco-friendly way is vital for Fenton technology's viability in pollution treatment. This study focuses on sewage sludge across various treatment stages, including generation, concentration, dehydration, and landfill, and employs chemical composite MoS2 to facilitate green resource utilization of all types of sludge. MoS2, with exposed Mo4+ and low-coordination sulfur, enhances iron cycling and creates an acidic microenvironment on the sludge surface. The MoS2-modified iron sludge exhibits outstanding (>95%) phenol and pollutant degradation in hydrogen peroxide and peroxymonosulfate-based Fenton systems, unlike unmodified sludge. This modified sludge maintains excellent Fenton activity in various water conditions and with multiple anions, allowing extended phenol degradation for over 14 d. Notably, the generated chemical oxygen demand (COD) in sludge modification process can be efficiently eliminated through the Fenton reaction, ensuring effluent COD compliance and enabling eco-friendly sewage sludge resource utilization.

Nutritional control of developmental processes.

Nutritional and metabolic cues are integral to animal development. Organisms use them both as sustenance and environmental indicators, fueling, informing and influencing developmental decisions. Classical examples, such as the Warburg effect, clearly illustrate how genetic programs control metabolic changes. However, the way that nutrition and metabolism can also modulate or drive genetic programs to instruct developmental trajectories is much more elusive, owing to several difficulties including uncoupling permissive and instructive functions. Here, we discuss recent advancements in the field that highlight the developmental role of nutritional and metabolic cues across multiple levels of organismal complexity.

Resource theory of quantum scrambling.

Quantum chaos has become a cornerstone of physics through its many applications. One trademark of quantum chaotic systems is the spread of local quantum information, which physicists call scrambling. In this work, we introduce a mathematical definition of scrambling and a resource theory to measure it. We also describe two applications of this theory. First, we use our resource theory to provide a bound on magic, a potential source of quantum computational advantage, which can be efficiently measured in experiment. Second, we also show that scrambling resources bound the success of Yoshida's black hole decoding protocol.

Trade-offs between cost and information in cellular prediction.

Living cells can leverage correlations in environmental fluctuations to predict the future environment and mount a response ahead of time. To this end, cells need to encode the past signal into the output of the intracellular network from which the future input is predicted. Yet, storing information is costly while not all features of the past signal are equally informative on the future input signal. Here, we show for two classes of input signals that cellular networks can reach the fundamental bound on the predictive information as set by the information extracted from the past signal: Push-pull networks can reach this information bound for Markovian signals, while networks that take a temporal derivative can reach the bound for predicting the future derivative of non-Markovian signals. However, the bits of past information that are most informative about the future signal are also prohibitively costly. As a result, the optimal system that maximizes the predictive information for a given resource cost is, in general, not at the information bound. Applying our theory to the chemotaxis network of Escherichia coli reveals that its adaptive kernel is optimal for predicting future concentration changes over a broad range of background concentrations, and that the system has been tailored to predicting these changes in shallow gradients.

Tradeoffs in bacterial physiology determine the efficiency of antibiotic killing.

Antibiotic effectiveness depends on a variety of factors. While many mechanistic details of antibiotic action are known, the connection between death rate and bacterial physiology is poorly understood. A common observation is that death rate in antibiotics rises linearly with growth rate; however, it remains unclear how other factors, such as environmental conditions and whole-cell physiological properties, affect bactericidal activity. To address this, we developed a high-throughput assay to precisely measure antibiotic-mediated death. We found that death rate is linear in growth rate, but the slope depends on environmental conditions. Growth under stress lowers death rate compared to nonstressed environments with similar growth rate. To understand stress's role, we developed a mathematical model of bacterial death based on resource allocation that includes a stress-response sector; we identify this sector using RNA-seq. Our model accurately predicts the minimal inhibitory concentration (MIC) with zero free parameters across a wide range of growth conditions. The model also quantitatively predicts death and MIC when sectors are experimentally modulated using cyclic adenosine monophosphate (cAMP), including protection from death at very low cAMP levels. The present study shows that different conditions with equal growth rate can have different death rates and establishes a quantitative relation between growth, death, and MIC that suggests approaches to improve antibiotic efficacy.

Microbial population dynamics decouple growth response from environmental nutrient concentration.

How the growth rate of a microbial population responds to the environmental availability of chemical nutrients and other resources is a fundamental question in microbiology. Models of this response, such as the widely used Monod model, are generally characterized by a maximum growth rate and a half-saturation concentration of the resource. What values should we expect for these half-saturation concentrations, and how should they depend on the environmental concentration of the resource? We survey growth response data across a wide range of organisms and resources. We find that the half-saturation concentrations vary across orders of magnitude, even for the same organism and resource. To explain this variation, we develop an evolutionary model to show that demographic fluctuations (genetic drift) can constrain the adaptation of half-saturation concentrations. We find that this effect fundamentally differs depending on the type of population dynamics: Populations undergoing periodic bottlenecks of fixed size will adapt their half-saturation concentrations in proportion to the environmental resource concentrations, but populations undergoing periodic dilutions of fixed size will evolve half-saturation concentrations that are largely decoupled from the environmental concentrations. Our model not only provides testable predictions for laboratory evolution experiments, but it also reveals how an evolved half-saturation concentration may not reflect the organism's environment. In particular, this explains how organisms in resource-rich environments can still evolve fast growth at low resource concentrations. Altogether, our results demonstrate the critical role of population dynamics in shaping fundamental ecological traits.

Resource competition can explain simplicity in microbial community assembly.

Predicting the composition and diversity of communities is a central goal in ecology. While community assembly is considered hard to predict, laboratory microcosms often follow a simple assembly rule based on the outcome of pairwise competitions. This assembly rule predicts that a species that is excluded by another species in pairwise competition cannot survive in a multispecies community with that species. Despite the empirical success of this bottom-up prediction, its mechanistic origin has remained elusive. In this study, we elucidate how this simple pattern in community assembly can emerge from resource competition. Our geometric analysis of a consumer-resource model shows that trio community assembly is always predictable from pairwise outcomes when one species grows faster than another species on every resource. We also identify all possible trio assembly outcomes under three resources and find that only two outcomes violate the assembly rule. Simulations demonstrate that pairwise competitions accurately predict trio assembly with up to 100 resources and the assembly of larger communities containing up to twelve species. We then further demonstrate accurate quantitative prediction of community composition using the harmonic mean of pairwise fractions. Finally, we show that cross-feeding between species does not decrease assembly rule prediction accuracy. Our findings highlight that simple community assembly can emerge even in ecosystems with complex underlying dynamics.

Expanding the BonnMu sequence-indexed repository of transposon induced maize (Zea mays L.) mutations in dent and flint germplasm.

The BonnMu resource is a transposon tagged mutant collection designed for functional genomics studies in maize. To expand this resource, we crossed an active Mutator (Mu) stock with dent (B73, Co125) and flint (DK105, EP1, and F7) germplasm, resulting in the generation of 8064 mutagenized BonnMu F2-families. Sequencing of these Mu-tagged families revealed 425 924 presumptive heritable Mu insertions affecting 36 612 (83%) of the 44 303 high-confidence gene models of maize (B73v5). On average, we observed 12 Mu insertions per gene (425 924 total insertions/36 612 affected genes) and 53 insertions per BonnMu F2-family (425 924 total insertions/8064 families). Mu insertions and photos of seedling phenotypes from segregating BonnMu F2-families can be accessed through the Maize Genetics and Genomics Database (MaizeGDB). Downstream examination via the automated Mutant-seq Workflow Utility (MuWU) identified 94% of the presumptive germinal insertion sites in genic regions and only a small fraction of 6% inserting in non-coding intergenic sequences of the genome. Consistently, Mu insertions aligned with gene-dense chromosomal arms. In total, 42% of all BonnMu insertions were located in the 5' untranslated region of genes, corresponding to accessible chromatin. Furthermore, for 38% of the insertions (163 843 of 425 924 total insertions) Mu1, Mu8 and MuDR were confirmed to be the causal Mu elements. Our publicly accessible European BonnMu resource has archived insertions covering two major germplasm groups, thus facilitating both forward and reverse genetics studies.

Functional analysis of regA paralog rlsD in Volvox carteri.

Volvox carteri is an excellent system for investigating the origins of cell differentiation because it possesses just two cell types, reproductive gonidia and motile somatic cells, which evolved relatively recently. The somatic phenotype depends on the regA gene, which represses cell growth and reproduction, preventing cells expressing it from growing large enough to become gonidia. regA encodes a putative transcription factor and was generated in an undifferentiated ancestor of V. carteri through duplication of a progenitor gene whose ortholog in V. carteri is named rlsD. Here we analyze the function of rlsD through knockdown, overexpression, and RNA-seq experiments, to gain clues into the function of a member of an understudied putative transcription factor family and to obtain insight into the origins of cell differentiation in the volvocine algae. rlsD knockdown was lethal, while rlsD overexpression dramatically reduced gonidial growth. rlsD overexpression led to differential expression of approximately one-fourth of the genome, with repressed genes biased for those typically overexpressed in gonidia relative to somatic cells, and upregulated genes biased toward expression in soma, where regA expression is high. Notably, rlsD overexpression affects accumulation of transcripts for genes/Pfam domains involved in ribosome biogenesis, photosynthetic light harvesting, and sulfate generation, functions related to organismal growth, and responses to resource availability. We also found that in the wild type, rlsD expression is induced by light deprivation. These findings are consistent with the idea that cell differentiation in V. carteri evolved when a resource-responsive, growth-regulating gene was amplified, and a resulting gene duplicate was co-opted to repress growth in a constitutive, spatial context.

From single- to multi-omics: future research trends in medicinal plants.

Medicinal plants are the main source of natural metabolites with specialised pharmacological activities and have been widely examined by plant researchers. Numerous omics studies of medicinal plants have been performed to identify molecular markers of species and functional genes controlling key biological traits, as well as to understand biosynthetic pathways of bioactive metabolites and the regulatory mechanisms of environmental responses. Omics technologies have been widely applied to medicinal plants, including as taxonomics, transcriptomics, metabolomics, proteomics, genomics, pangenomics, epigenomics and mutagenomics. However, because of the complex biological regulation network, single omics usually fail to explain the specific biological phenomena. In recent years, reports of integrated multi-omics studies of medicinal plants have increased. Until now, there have few assessments of recent developments and upcoming trends in omics studies of medicinal plants. We highlight recent developments in omics research of medicinal plants, summarise the typical bioinformatics resources available for analysing omics datasets, and discuss related future directions and challenges. This information facilitates further studies of medicinal plants, refinement of current approaches and leads to new ideas.

Trade-offs between the instantaneous growth rate and long-term fitness: Consequences for microbial physiology and predictive computational models.

Microbial systems biology has made enormous advances in relating microbial physiology to the underlying biochemistry and molecular biology. By meticulously studying model microorganisms, in particular Escherichia coli and Saccharomyces cerevisiae, increasingly comprehensive computational models predict metabolic fluxes, protein expression, and growth. The modeling rationale is that cells are constrained by a limited pool of resources that they allocate optimally to maximize fitness. As a consequence, the expression of particular proteins is at the expense of others, causing trade-offs between cellular objectives such as instantaneous growth, stress tolerance, and capacity to adapt to new environments. While current computational models are remarkably predictive for E. coli and S. cerevisiae when grown in laboratory environments, this may not hold for other growth conditions and other microorganisms. In this contribution, we therefore discuss the relationship between the instantaneous growth rate, limited resources, and long-term fitness. We discuss uses and limitations of current computational models, in particular for rapidly changing and adverse environments, and propose to classify microbial growth strategies based on Grimes's CSR framework.