Statistical Uncertainty in Paleoclimate Proxy Reconstructions.

A quantitative analysis of any environment older than the instrumental record relies on proxies. Uncertainties associated with proxy reconstructions are often underestimated, which can lead to artificial conflict between different proxies, and between data and models. In this paper, using ordinary least squares linear regression as a common example, we describe a simple, robust and generalizable method for quantifying uncertainty in proxy reconstructions. We highlight the primary controls on the magnitude of uncertainty, and compare this simple estimate to equivalent estimates from Bayesian, nonparametric and fiducial statistical frameworks. We discuss when it may be possible to reduce uncertainties, and conclude that the unexplained variance in the calibration must always feature in the uncertainty in the reconstruction. This directs future research toward explaining as much of the variance in the calibration data as possible. We also advocate for a "data-forward" approach, that clearly decouples the presentation of proxy data from plausible environmental inferences.

Mesoscopic Dynamical Differences from Quantum State Preparation in a Bose-Hubbard Trimer.

Conventional wisdom is that quantum effects will tend to disappear as the number of quanta in a system increases, and the evolution of a system will become closer to that described by mean-field classical equations. In this Letter we combine newly developed theoretical and experimental techniques to propose and analyze an experiment using a Bose-Hubbard trimer where the opposite is the case. We find that differences in the preparation of a centrally evacuated trimer can lead to readily observable differences in the subsequent dynamics which increase with system size. Importantly, these differences can be detected by the simple measurements of atomic number.

Identifying a Superfluid Reynolds Number via Dynamical Similarity.

The Reynolds number provides a characterization of the transition to turbulent flow, with wide application in classical fluid dynamics. Identifying such a parameter in superfluid systems is challenging due to their fundamentally inviscid nature. Performing a systematic study of superfluid cylinder wakes in two dimensions, we observe dynamical similarity of the frequency of vortex shedding by a cylindrical obstacle. The universality of the turbulent wake dynamics is revealed by expressing shedding frequencies in terms of an appropriately defined superfluid Reynolds number, Re(s), that accounts for the breakdown of superfluid flow through quantum vortex shedding. For large obstacles, the dimensionless shedding frequency exhibits a universal form that is well-fitted by a classical empirical relation. In this regime the transition to turbulence occurs at Re(s)≈0.7, irrespective of obstacle width.

Onsager-Kraichnan condensation in decaying two-dimensional quantum turbulence.

Despite the prominence of Onsager's point-vortex model as a statistical description of 2D classical turbulence, a first-principles development of the model for a realistic superfluid has remained an open problem. Here we develop a mapping of a system of quantum vortices described by the homogeneous 2D Gross-Pitaevskii equation (GPE) to the point-vortex model, enabling Monte Carlo sampling of the vortex microcanonical ensemble. We use this approach to survey the full range of vortex states in a 2D superfluid, from the vortex-dipole gas at positive temperature to negative-temperature states exhibiting both macroscopic vortex clustering and kinetic energy condensation, which we term an Onsager-Kraichnan condensate (OKC). Damped GPE simulations reveal that such OKC states can emerge dynamically, via aggregation of small-scale clusters into giant OKC clusters, as the end states of decaying 2D quantum turbulence in a compressible, finite-temperature superfluid. These statistical equilibrium states should be accessible in atomic Bose-Einstein condensate experiments.

Altered functional B cell subset populations in patients with chronic fatigue syndrome compared to healthy controls.

Chronic fatigue syndrome (CFS) is a heterogeneous disorder of unknown aetiology characterized by disabling fatigue, headaches, sleep disturbance and several other symptoms. The onset of CFS may follow a viral infection or period of stress. Patients with CFS do not have hypogammaglobulinaemia, predisposition to recurrent bacterial infections or symptoms of autoimmunity. To date, defects in B cell numbers or function have not been shown in the literature. However, treatment with anti-B cell therapy using Rituximab has recently shown benefit to CFS patients. We therefore postulated that patients with CFS had a subtle humoral immune dysfunction, and performed extended B cell immunophenotyping. We undertook a detailed characterization of the proportions of the different B cell subsets in 33 patients with CFS fulfilling the Canadian and Fukada criteria for CFS and compared these with 24 age- and gender-matched healthy controls (HC). CFS patients had greater numbers of naive B cells as a percentage of lymphocytes: 6·3 versus 3·9% in HC (P = 0·034), greater numbers of naive B cells as a percentage of B cells: 65 versus 47% in controls (P = 0·003), greater numbers of transitional B cells: 1·8 versus 0·8% in controls (P = 0·025) and reduced numbers of plasmablasts: 0·5 versus 0·9% in controls (P = 0·013). While the cause of these changes is unclear, we speculate whether they may suggest a subtle tendency to autoimmunity.

Characteristics of two-dimensional quantum turbulence in a compressible superfluid.

Fluids subjected to suitable forcing will exhibit turbulence, with characteristics strongly affected by the fluid's physical properties and dimensionality. In this work, we explore two-dimensional (2D) quantum turbulence in an oblate Bose-Einstein condensate confined to an annular trapping potential. Experimentally, we find conditions for which small-scale stirring of the condensate generates disordered 2D vortex distributions that dissipatively evolve toward persistent currents, indicating energy transport from small to large length scales. Simulations of the experiment reveal spontaneous clustering of same-circulation vortices and an incompressible energy spectrum with k(-5/3) dependence for low wave numbers k. This work links experimentally observed vortex dynamics with signatures of 2D turbulence in a compressible superfluid.

Quantum tunneling of a vortex between two pinning potentials.

A vortex can tunnel between two pinning potentials in an atomic Bose-Einstein condensate on a time scale of the order of 1s under typical experimental conditions. This makes it possible to detect the tunneling experimentally. We calculate the tunneling rate by phenomenologically treating vortices as charged particles moving in an inhomogeneous magnetic field. The obtained results are in close agreement with numerical simulations based on the stochastic c-field theory.

Chronic fatigue syndrome, the immune system and viral infection.

The chronic fatigue syndrome (CFS), as defined by recent criteria, is a heterogeneous disorder with a common set of symptoms that often either follows a viral infection or a period of stress. Despite many years of intense investigation there is little consensus on the presence, nature and degree of immune dysfunction in this condition. However, slightly increased parameters of inflammation and pro-inflammatory cytokines such as interleukin (IL) 1, IL6 and tumour necrosis factor (TNF) α are likely present. Additionally, impaired natural killer cell function appears evident. Alterations in T cell numbers have been described by some and not others. While the prevalence of positive serology for the common herpes viruses appears no different from healthy controls, there is some evidence of viral persistence and inadequate containment of viral replication. The ability of certain herpes viruses to impair the development of T cell memory may explain this viral persistence and the continuation of symptoms. New therapies based on this understanding are more likely to produce benefit than current methods.

Observation of vortex dipoles in an oblate Bose-Einstein condensate.

We report experimental observations and numerical simulations of the formation, dynamics, and lifetimes of single and multiply charged quantized vortex dipoles in highly oblate dilute-gas Bose-Einstein condensates (BECs). We nucleate pairs of vortices of opposite charge (vortex dipoles) by forcing superfluid flow around a repulsive Gaussian obstacle within the BEC. By controlling the flow velocity we determine the critical velocity for the nucleation of a single vortex dipole, with excellent agreement between experimental and numerical results. We present measurements of vortex dipole dynamics, finding that the vortex cores of opposite charge can exist for many seconds and that annihilation is inhibited in our trap geometry. For sufficiently rapid flow velocities, clusters of like-charge vortices aggregate into long-lived multiply charged dipolar flow structures.

Direct Observation of the Dynamics of Single-Cell Metabolic Activity during Microbial Diauxic Growth.

Population-level analyses are rapidly becoming inadequate to answer many of biomedical science and microbial ecology's most pressing questions. The role of microbial populations within ecosystems and the evolutionary selective pressure on individuals depend fundamentally on the metabolic activity of single cells. Yet, many existing single-cell technologies provide only indirect evidence of metabolic specialization because they rely on correlations between transcription and phenotype established at the level of the population to infer activity. In this study, we take a top-down approach using isotope labels and secondary ion mass spectrometry to track the uptake of carbon and nitrogen atoms from different sources into biomass and directly observe dynamic changes in anabolic specialization at the level of single cells. We investigate the classic microbiological phenomenon of diauxic growth at the single-cell level in the model methylotroph Methylobacterium extorquens In nature, this organism inhabits the phyllosphere, where it experiences diurnal changes in the available carbon substrates, necessitating an overhaul of central carbon metabolism. We show that the population exhibits a unimodal response to the changing availability of viable substrates, a conclusion that supports the canonical model but has thus far been supported by only indirect evidence. We anticipate that the ability to monitor the dynamics of anabolism in individual cells directly will have important applications across the fields of ecology, medicine, and biogeochemistry, especially where regulation downstream of transcription has the potential to manifest as heterogeneity that would be undetectable with other existing single-cell approaches.IMPORTANCE Understanding how genetic information is realized as the behavior of individual cells is a long-term goal of biology but represents a significant technological challenge. In clonal microbial populations, variation in gene regulation is often interpreted as metabolic heterogeneity. This follows the central dogma of biology, in which information flows from DNA to RNA to protein and ultimately manifests as activity. At present, DNA and RNA can be characterized in single cells, but the abundance and activity of proteins cannot. Inferences about metabolic activity usually therefore rely on the assumption that transcription reflects activity. By tracking the atoms from which they build their biomass, we make direct observations of growth rate and substrate specialization in individual cells throughout a period of growth in a changing environment. This approach allows the flow of information from DNA to be constrained from the distal end of the regulatory cascade and will become an essential tool in the rapidly advancing field of single-cell metabolism.

Oxygen isotope effects during microbial sulfate reduction: applications to sediment cell abundances.

The majority of anaerobic biogeochemical cycling occurs within marine sediments. To understand these processes, quantifying the distribution of active cells and gross metabolic activity is essential. We present an isotope model rooted in thermodynamics to draw quantitative links between cell-specific sulfate reduction rates and active sedimentary cell abundances. This model is calibrated using data from a series of continuous culture experiments with two strains of sulfate reducing bacteria (freshwater bacterium Desulfovibrio vulgaris strain Hildenborough, and marine bacterium Desulfovibrio alaskensis strain G-20) grown on lactate across a range of metabolic rates and ambient sulfate concentrations. We use a combination of experimental sulfate oxygen isotope data and nonlinear regression fitting tools to solve for unknown kinetic, step-specific oxygen isotope effects. This approach enables identification of key isotopic reactions within the metabolic pathway, and defines a new, calibrated framework for understanding oxygen isotope variability in sulfate. This approach is then combined with porewater sulfate/sulfide concentration data and diagenetic modeling to reproduce measured 18O/16O in porewater sulfate. From here, we infer cell-specific sulfate reduction rates and predict abundance of active cells of sulfate reducing bacteria, the result of which is consistent with direct biological measurements.

Numerical method for evolving the dipolar projected Gross-Pitaevskii equation.

We describe a method for evolving the projected Gross-Pitaevskii equation (PGPE) for an interacting Bose gas in a harmonic-oscillator potential, with the inclusion of a long-range dipolar interaction. The central difficulty in solving this equation is the requirement that the field is restricted to a small set of prescribed modes that constitute the low-energy c -field region of the system. We present a scheme, using a Hermite-polynomial-based spectral representation, which precisely implements this mode restriction and allows an efficient and accurate solution of the dipolar PGPE. We introduce a set of auxiliary oscillator states to perform a Fourier transform necessary to evaluate the dipolar interaction in reciprocal space. We extensively characterize the accuracy of our approach and derive Ehrenfest equations for the evolution of the angular momentum.

Patterns of sulfur isotope fractionation during microbial sulfate reduction.

Studies of microbial sulfate reduction have suggested that the magnitude of sulfur isotope fractionation varies with sulfate concentration. Small apparent sulfur isotope fractionations preserved in Archean rocks have been interpreted as suggesting Archean sulfate concentrations of <200 µm, while larger fractionations thereafter have been interpreted to require higher concentrations. In this work, we demonstrate that fractionation imposed by sulfate reduction can be a function of concentration over a millimolar range, but that nature of this relationship depends on the organism studied. Two sulfate-reducing bacteria grown in continuous culture with sulfate concentrations ranging from 0.1 to 6 mm showed markedly different relationships between sulfate concentration and isotope fractionation. Desulfovibrio vulgaris str. Hildenborough showed a large and relatively constant isotope fractionation ((34) εSO 4-H2S ≅ 25‰), while fractionation by Desulfovibrio alaskensis G20 strongly correlated with sulfate concentration over the same range. Both data sets can be modeled as Michaelis-Menten (MM)-type relationships but with very different MM constants, suggesting that the fractionations imposed by these organisms are highly dependent on strain-specific factors. These data reveal complexity in the sulfate concentration-fractionation relationship. Fractionation during MSR relates to sulfate concentration but also to strain-specific physiological parameters such as the affinity for sulfate and electron donors. Previous studies have suggested that the sulfate concentration-fractionation relationship is best described with a MM fit. We present a simple model in which the MM fit with sulfate concentration and hyperbolic fit with growth rate emerge from simple physiological assumptions. As both environmental and biological factors influence the fractionation recorded in geological samples, understanding their relationship is critical to interpreting the sulfur isotope record. As the uptake machinery for both sulfate and electrons has been subject to selective pressure over Earth history, its evolution may complicate efforts to uniquely reconstruct ambient sulfate concentrations from a single sulfur isotopic composition.

Numerical method for the stochastic projected Gross-Pitaevskii equation.

We present a method for solving the stochastic projected Gross-Pitaevskii equation (SPGPE) for a three-dimensional weakly interacting Bose gas in a harmonic-oscillator trapping potential. The SPGPE contains the challenge of both accurately evolving all modes in the low-energy classical region of the system, and evaluating terms from the number-conserving scattering reservoir process. We give an accurate and efficient procedure for evaluating the scattering terms using a Hermite-polynomial based spectral-Galerkin representation, which allows us to precisely implement the low-energy mode restriction. Stochastic integration is performed using the weak semi-implicit Euler method. We extensively characterize the accuracy of our method, finding a faster-than-expected rate of stochastic convergence. Physical consistency of the algorithm is demonstrated by considering thermalization of initially random states.

Spatial and temporal variability of biomarkers and microbial diversity reveal metabolic and community flexibility in Streamer Biofilm Communities in the Lower Geyser Basin, Yellowstone National Park.

Detailed analysis of 16S rRNA and intact polar lipids (IPLs) from streamer biofilm communities (SBCs), collected from geochemically similar hot springs in the Lower Geyser Basin, Yellowstone National Park, shows good agreement and affirm that IPLs can be used as reliable markers for the microbial constituents of SBCs. Uncultured Crenarchaea are prominent in SBS, and their IPLs contain both glycosidic and mixed glyco-phospho head groups with tetraether cores, having 0-4 rings. Archaeal IPL contributions increase with increasing temperature and comprise up to one-fourth of the total IPL inventory at >84 °C. At elevated temperatures, bacterial IPLs contain abundant glycosidic glycerol diether lipids. Diether and diacylglycerol (DAG) lipids with aminopentanetetrol and phosphatidylinositol head groups were identified as lipids diagnostic of Aquificales, while DAG glycolipids and glyco-phospholipids containing N-acetylgycosamine as head group were assigned to members of the Thermales. With decreasing temperature and concomitant changes in water chemistry, IPLs typical of phototrophic bacteria, such as mono-, diglycosyl, and sulfoquinovosyl DAG, which are specific for cyanobacteria, increase in abundance, consistent with genomic data from the same samples. Compound-specific stable carbon isotope analysis of IPL breakdown products reveals a large isotopic diversity among SBCs in different hot springs. At two of the hot springs, 'Bison Pool' and Flat Cone, lipids derived from Aquificales are enriched in (13) C relative to biomass and approach values close to dissolved inorganic carbon (DIC) (approximately 0‰), consistent with fractionation during autotrophic carbon fixation via the reversed tricarboxylic acid pathway. At a third site, Octopus Spring, the same Aquificales-diagnostic lipids are 10‰ depleted relative to biomass and resemble stable carbon isotope values of dissolved organic carbon (DOC), indicative of heterotrophy. Other bacterial and archaeal lipids show a similar variance, with values resembling the DIC or DOC pool or a mixture thereof. This variance cannot be explained by hot spring chemistry or temperature alone, but instead, we argue that intermittent input of exogenous organic carbon can result in metabolic shifts of the chemotrophic communities from autotrophy to heterotrophy and vice versa.

Revisiting the dissimilatory sulfate reduction pathway.

Sulfur isotopes in the geological record integrate a combination of biological and diagenetic influences, but a key control on the ratio of sulfur isotopes in sedimentary materials is the magnitude of isotope fractionation imparted during dissimilatory sulfate reduction. This fractionation is controlled by the flux of sulfur through the network of chemical reactions involved in sulfate reduction and by the isotope effect associated with each of these chemical reactions. Despite its importance, the network of reactions constituting sulfate reduction is not fully understood, with two principle networks underpinning most isotope models. In this study, we build on biochemical data and recently solved crystal structures of enzymes to propose a revised network topology for the flow of sulfur through the sulfate reduction metabolism. This network is highly branched and under certain conditions produces results consistent with the observations that motivated previous sulfate reduction models. Our revised network suggests that there are two main paths to sulfide production: one that involves the production of thionate intermediates, and one that does not. We suggest that a key factor in determining sulfur isotope fractionation associated with sulfate reduction is the ratio of the rate at which electrons are supplied to subunits of Dsr vs. the rate of sulfite delivery to the active site of Dsr. This reaction network may help geochemists to better understand the relationship between the physiology of sulfate reduction and the isotopic record it produces.

Proteorhodopsin photosystem gene expression enables photophosphorylation in a heterologous host.

Proteorhodopsins (PRs) are retinal-containing proteins that catalyze light-activated proton efflux across the cell membrane. These photoproteins are known to be globally distributed in the ocean's photic zone, and they are found in a diverse array of Bacteria and Archaea. Recently, light-enhanced growth rates and yields have been reported in at least one PR-containing marine bacterium, but the physiological basis of light-activated growth stimulation has not yet been determined. To describe more fully PR photosystem genetics and biochemistry, we functionally surveyed a marine picoplankton large-insert genomic library for recombinant clones expressing PR photosystems in vivo. Our screening approach exploited transient increases in vector copy number that significantly enhanced the sensitivity of phenotypic detection. Two genetically distinct recombinants, initially identified by their orange pigmentation, expressed a small cluster of genes encoding a complete PR-based photosystem. Genetic and biochemical analyses of transposon mutants verified the function of gene products in the photopigment and opsin biosynthetic pathways. Heterologous expression of six genes, five encoding photopigment biosynthetic proteins and one encoding a PR, generated a fully functional PR photosystem that enabled photophosphorylation in recombinant Escherichia coli cells exposed to light. Our results demonstrate that a single genetic event can result in the acquisition of phototrophic capabilities in an otherwise chemoorganotrophic microorganism, and they explain in part the ubiquity of PR photosystems among diverse microbial taxa.