Combining stereo‐video monitoring and physiological trials to estimate reef fish metabolic demands in the wild

Organismal metabolic rates (MRs) are the basis of energy and nutrient fluxes through ecosystems. In the marine realm, fishes are some of the most prominent consumers. However, their metabolic demand in the wild (field MR [FMR]) is poorly documented, because it is challenging to measure directly. Here, we introduce a novel approach to estimating the component of FMR associated with voluntary activity (i.e., the field active MR [AMRfield]). Our approach combines laboratory‐based respirometry, swimming speeds, and field‐based stereo‐video systems to estimate the activity of individuals. We exemplify our approach by focusing on six coral reef fish species, for which we quantified standard MR and maximum MR (SMR and MMR, respectively) in the laboratory, and body sizes and swimming speeds in the field. Based on the relationships between MR, body size, and swimming speeds, we estimate that the activity scope (i.e., the ratio between AMRfield and SMR) varies from 1.2 to 3.2 across species and body sizes. Furthermore, we illustrate that the scaling exponent for AMRfield varies across species and can substantially exceed the widely assumed value of 0.75 for SMR. Finally, by scaling organismal AMRfield estimates to the assemblage level, we show the potential effect of this variability on community metabolic demand. Our approach may improve our ability to estimate elemental fluxes mediated by a critically important group of aquatic animals through a non‐destructive, widely applicable technique.     

Importance of the description of light interception in crop growth models

Canopy light interception determines the amount of energy captured by a crop, and is thus critical to modeling crop growth and yield, and may substantially contribute to the prediction uncertainty of crop growth models (CGMs). We thus analyzed the canopy light interception models of the 26 wheat (Triticum aestivum) CGMs used by the Agricultural Model Intercomparison and Improvement Project (AgMIP). Twenty-one CGMs assume that the light extinction coefficient (K) is constant, varying from 0.37 to 0.80 depending on the model. The other models take into account the illumination conditions and assume either that all green surfaces in the canopy have the same inclination angle (θ) or that θ distribution follows a spherical distribution. These assumptions have not yet been evaluated due to a lack of experimental data. Therefore, we conducted a field experiment with five cultivars with contrasting leaf stature sown at normal and double row spacing, and analyzed θ distribution in the canopies from three-dimensional canopy reconstructions. In all the canopies, θ distribution was well represented by an ellipsoidal distribution. We thus carried out an intercomparison between the light interception models of the AgMIP–Wheat CGMs ensemble and a physically based K model with ellipsoidal leaf angle distribution and canopy clumping (KellC). Results showed that the KellC model outperformed current approaches under most illumination conditions and that the uncertainty in simulated wheat growth and final grain yield due to light models could be as high as 45%. Therefore, our results call for an overhaul of light interception models in CGMs.     

Demographic effects of a megafire on a declining prairie grouse in the mixed‐grass prairie

Recent studies have documented benefits of small, prescribed fire and wildfire for grassland‐dependent wildlife, such as lesser prairie‐chickens (Tympanuchus pallidicintus), but wildlife demographic response to the scale and intensity of megafire (wildfire >40,000 ha) in modern, fragmented grasslands remains unknown. Limited available grassland habitat makes it imperative to understand if increasing frequency of megafires could further reduce already declining lesser prairie‐chicken populations, or if historical evolutionary interactions with fire make lesser prairie‐chickens resilient. To evaluate lesser prairie‐chicken demographic response to megafires, we compared lek counts, nest density, and survival rates of adults, nests, and chicks before (2014–2016) and after (2018–2020) a 2017 megafire in the mixed‐grass prairie of Kansas, USA (Starbuck fire ~254,000 ha). There was a 67% decline in attending males on leks post‐fire and a 57% decline in occupied leks post‐fire. Despite population declines as indicated by lek counts, adult female breeding season survival (S^) was similar pre‐ (S^ = 0.65 ± 0.08 [SE]) and post‐fire (0.61 ± 0.08), as was chick survival (pre‐fire: 0.23 ± 0.07; post‐fire: 0.27 ± 0.11). Nest survival appeared lower post‐fire (pre‐fire: 0.38 ± 0.06; post‐fire: 0.20 ± 0.06), but did not differ at the 95% confidence interval. Nest density of marked females declined 73% in areas burned by megafire. Although lesser prairie‐chickens persisted in the study area and we documented minimal effects on most demographic rates, reduced lesser prairie‐chicken abundance and reproductive output suggests full recovery may take >3 years. Increased propensity for megafire resulting from suppression of smaller fires, compounded by climate change and woody encroachment, may impose a short‐term (3–5 year) threat to already declining lesser prairie‐chicken populations.     

Dominance Genetic Variation Contributes Little to the Missing Heritability for Human Complex Traits

For human complex traits, non-additive genetic variation has been invoked to explain “missing heritability,” but its discovery is often neglected in genome-wide association studies. Here we propose a method of using SNP data to partition and estimate the proportion of phenotypic variance attributed to additive and dominance genetic variation at all SNPs (hSNP2 and δSNP2) in unrelated individuals based on an orthogonal model where the estimate of hSNP2 is independent of that of δSNP2. With this method, we analyzed 79 quantitative traits in 6,715 unrelated European Americans. The estimate of δSNP2 averaged across all the 79 quantitative traits was 0.03, approximately a fifth of that for additive variation (average hSNP2 = 0.15). There were a few traits that showed substantial estimates of δSNP2, none of which were replicated in a larger sample of 11,965 individuals. We further performed genome-wide association analyses of the 79 quantitative traits and detected SNPs with genome-wide significant dominance effects only at the ABO locus for factor VIII and von Willebrand factor. All these results suggest that dominance variation at common SNPs explains only a small fraction of phenotypic variation for human complex traits and contributes little to the missing narrow-sense heritability problem.     

Patch utilization and flower visitations by wild bees in a honey bee‐dominated,grassland landscape

Understanding habitat needs and patch utilization of wild and managed bees has been identified as a national research priority in the United States. We used occupancy models to investigate patterns of bee use across 1030 transects spanning a gradient of floral resource abundance and richness and distance from apiaries in the Prairie Pothole Region (PPR) of the United States. Estimates of transect use by honey bees were nearly 1.0 during our 3.5‐month sampling period, suggesting honey bees were nearly ubiquitous across transects. Wild bees more frequently used transects with higher flower richness and more abundant flowers; however, the effect size of the native flower abundance covariate (β^native = 3.90 ± 0.65 [1SE]) was four times greater than the non‐native flower covariate (β^non‐native = 0.99 ± 0.17). We found some evidence that wild bee use was lower at transects near commercial apiaries, but the effect size was imprecise (β^distance = 1.4 ± 0.81). Honey bees were more frequently detected during sampling events with more non‐native flowers and higher species richness but showed an uncertain relationship with native flower abundance. Of the 4039 honey bee and flower interactions, 85% occurred on non‐native flowers, while only 43% of the 738 wild bee observations occurred on non‐native flowers. Our study suggests wild bees and honey bees routinely use the same resource patches in the PPR but often visit different flowering plants. The greatest potential for resource overlap between honey bees and wild bees appears to be for non‐native flowers in the PPR. Our results are valuable to natural resource managers tasked with supporting habitat for managed and wild pollinators in agroecosystems.     

A Variable Sampling Interval Synthetic Xbar Chart for the Process Mean

The usual practice of using a control chart to monitor a process is to take samples from the process with fixed sampling interval (FSI). In this paper, a synthetic X¯ control chart with the variable sampling interval (VSI) feature is proposed for monitoring changes in the process mean. The VSI synthetic X¯ chart integrates the VSI X¯ chart and the VSI conforming run length (CRL) chart. The proposed VSI synthetic X¯ chart is evaluated using the average time to signal (ATS) criterion. The optimal charting parameters of the proposed chart are obtained by minimizing the out-of-control ATS for a desired shift. Comparisons between the VSI synthetic X¯ chart and the existing X¯, synthetic X¯, VSI X¯ and EWMA X¯ charts, in terms of ATS, are made. The ATS results show that the VSI synthetic X¯ chart outperforms the other X¯ type charts for detecting moderate and large shifts. An illustrative example is also presented to explain the application of the VSI synthetic X¯ chart.     

Higher ventilatory responses during and after passive walking-like leg movement in older individuals

Background

Minute ventilation (V·E) during walking has been shown to be higher in older individuals than in young individuals, but the mechanisms underlying the higher ventilatory response is unclear. Central command and peripheral neural reflex are important neural control mechanisms underlying ventilatory response during exercise. Passive leg movement has been used to exclude the influence of central command due to the lack of voluntary activation of muscles. The aim of the present study was to compare the ventilatory response during and after passive walking-like leg movement (PWM) in young and older individuals.

Methods

Eight young subjects (20 ± 2 years) and seven older subjects (70 ± 1 years) participated in this study. Subjects spent 7 minutes in a quiet standing (QS) position. Thereafter, they performed 14-minute rhythmic PWM at 1 Hz and this was followed by 7 minutes of QS.

Results

V·E values during pre-PWM QS were calculated as 1-minute averages using data obtained between 5 and 6 minutes. V·E values at pre-PWM QS in the young and older groups were 8.4 ± 2.1 and 7.5 ± 1.2 l/minute, respectively. V·E values increased significantly at the first minute of PWM to 11.4 ± 2.2 and 10.4 ± 2.5 l/minute in the young and older groups, respectively (P <0.001). In the young group, V·E at the last minute of PWM (9.2 ± 2.0 l/minute) was not significantly different from that at pre-PWM QS due to a decline in V·E, whereas V·E at the last minute of PWM in the older group (9.4 ± 2.2 l/minute) was still significantly higher (P <0.01). On the other hand, V·E at the first minute of post-PWM QS (7.2 ± 1.8 l/minute) was significantly lower than that during pre-PWM QS in the young group (P <0.05) but not in the older group.

Conclusions

Ventilatory response during and after PWM is higher in older individuals than in young individuals. This may be associated with a mechanism(s) other than central command. Our findings may explain part of the higher V·E response while walking in older individuals.     

Beta-Adrenergic Stimulation Maintains Cardiac Function in Serca2 Knockout Mice

Previous studies on Serca2 knockout (KO) mice showed that cardiac function is sustained in vivo for several weeks after knockout, whereas SERCA protein levels decrease and calcium dynamics are significantly impaired. In this study, we reconcile observed cellular and organ level contractile function using a cardiac multiscale model. We identified and quantified the changes in cellular function that are both consistent with observations and able to compensate for the decrease in SERCA. Calcium transients were used as input for multiscale computational simulations to predict whole-organ response. Although this response matched experimental pressure-volume (PV) measurements in healthy mice, the reduced magnitude calcium transients observed in KO cells were insufficient to trigger ventricular ejection. To replicate the effects of elevated catecholamine levels observed in vivo, cells were treated with isoproterenol. Incorporation of the resulting measured β-adrenergically stimulated calcium transients into the model resulted in a close match with experimental PV loops. Changes in myofilament properties, when considered in isolation, were not able to increase tension development to levels consistent with measurements, further confirming the necessity of a high β-adrenergic state. Modeling additionally indicated that increased venous return observed in the KO mice helps maintain a high ejection fraction via the Frank-Starling effect. Our study shows that increased β-adrenergic stimulation is a potentially highly significant compensatory mechanism by which cardiac function is maintained in Serca2 KO mice, producing the increases in both systolic and diastolic calcium, consistent with the observed contractile function observed in experimental PV measurements.     

Error-based or target-based? A unified framework for learning in recurrent spiking networks

The field of recurrent neural networks is over-populated by a variety of proposed learning rules and protocols. The scope of this work is to define a generalized framework, to move a step forward towards the unification of this fragmented scenario. In the field of supervised learning, two opposite approaches stand out, error-based and target-based. This duality gave rise to a scientific debate on which learning framework is the most likely to be implemented in biological networks of neurons. Moreover, the existence of spikes raises the question of whether the coding of information is rate-based or spike-based. To face these questions, we proposed a learning model with two main parameters, the rank of the feedback learning matrix R and the tolerance to spike timing τ_⋆. We demonstrate that a low (high) rank R accounts for an error-based (target-based) learning rule, while high (low) tolerance to spike timing promotes rate-based (spike-based) coding. We show that in a store and recall task, high-ranks allow for lower MSE values, while low-ranks enable a faster convergence. Our framework naturally lends itself to Behavioral Cloning and allows for efficiently solving relevant closed-loop tasks, investigating what parameters (R,τ⋆) are optimal to solve a specific task. We found that a high R is essential for tasks that require retaining memory for a long time (Button and Food). On the other hand, this is not relevant for a motor task (the 2D Bipedal Walker). In this case, we find that precise spike-based coding enables optimal performances. Finally, we show that our theoretical formulation allows for defining protocols to estimate the rank of the feedback error in biological networks. We release a PyTorch implementation of our model supporting GPU parallelization.     

The Impact of Prior Information on Estimates of Disease Transmissibility Using Bayesian Tools

The basic reproductive number (R₀) and the distribution of the serial interval (SI) are often used to quantify transmission during an infectious disease outbreak. In this paper, we present estimates of R₀ and SI from the 2003 SARS outbreak in Hong Kong and Singapore, and the 2009 pandemic influenza A(H1N1) outbreak in South Africa using methods that expand upon an existing Bayesian framework. This expanded framework allows for the incorporation of additional information, such as contact tracing or household data, through prior distributions. The results for the R₀ and the SI from the influenza outbreak in South Africa were similar regardless of the prior information (R^0 = 1.36–1.46, μ^ = 2.0–2.7, μ^ = mean of the SI). The estimates of R₀ and μ for the SARS outbreak ranged from 2.0–4.4 and 7.4–11.3, respectively, and were shown to vary depending on the use of contact tracing data. The impact of the contact tracing data was likely due to the small number of SARS cases relative to the size of the contact tracing sample.     

The catalytic mechanism of the mitochondrial methylenetetrahydrofolate dehydrogenase/cyclohydrolase (MTHFD2)

Methylenetetrahydrofolate dehydrogenase/cyclohydrolase (MTHFD2) is a new drug target that is expressed in cancer cells but not in normal adult cells, which provides an Achilles heel to selectively kill cancer cells. Despite the availability of crystal structures of MTHFD2 in the inhibitor- and cofactor-bound forms, key information is missing due to technical limitations, including (a) the location of absolutely required Mg²⁺ ion, and (b) the substrate-bound form of MTHFD2. Using computational modeling and simulations, we propose that two magnesium ions are present at the active site whereby (i) Arg233, Asp225, and two water molecules coordinate MgA2+, while MgA2+ together with Arg233 stabilize the inorganic phosphate (P_i); (ii) Asp168 and three water molecules coordinate MgB2+, and MgB2+ further stabilizes P_i by forming a hydrogen bond with two oxygens of P_i; (iii) Arg201 directly coordinates the P_i; and (iv) through three water-mediated interactions, Asp168 contributes to the positioning and stabilization of MgA2+, MgB2+ and P_i. Our computational study at the empirical valence bond level allowed us also to elucidate the detailed reaction mechanisms. We found that the dehydrogenase activity features a proton-coupled electron transfer with charge redistribution connected to the reorganization of the surrounding water molecules which further facilitates the subsequent cyclohydrolase activity. The cyclohydrolase activity then drives the hydration of the imidazoline ring and the ring opening in a concerted way. Furthermore, we have uncovered that two key residues, Ser197/Arg233, are important factors in determining the cofactor (NADP⁺/NAD⁺) preference of the dehydrogenase activity. Our work sheds new light on the structural and kinetic framework of MTHFD2, which will be helpful to design small molecule inhibitors that can be used for cancer treatment.     

A Bayesian method for identifying associations between response variables and bacterial community composition

Determining associations between intestinal bacteria and continuously measured physiological outcomes is important for understanding the bacteria-host relationship but is not straightforward since abundance data (compositional data) are not normally distributed. To address this issue, we developed a fully Bayesian linear regression model (BRACoD; Bayesian Regression Analysis of Compositional Data) with physiological measurements (continuous data) as a function of a matrix of relative bacterial abundances. Bacteria can be classified as operational taxonomic units or by taxonomy (genus, family, etc.). Bacteria associated with the physiological measurement were identified using a Bayesian variable selection method: Stochastic Search Variable Selection. The output is a list of inclusion probabilities (p^) and coefficients that indicate the strength of the association (β^included) for each bacterial taxa. Tests with simulated communities showed that adopting a cut point value of p^ ≥ 0.3 for identifying included bacteria optimized the true positive rate (TPR) while maintaining a false positive rate (FPR) of ≤ 5%. At this point, the chances of identifying non-contributing bacteria were low and all well-established contributors were included. Comparison with other methods showed that BRACoD (at p^ ≥ 0.3) had higher precision and a higher TPR than a commonly used center log transformed LASSO procedure (clr-LASSO) as well as higher TPR than an off-the-shelf Spike and Slab method after center log transformation (clr-SS). BRACoD was also less likely to include non-contributing bacteria that merely correlate with contributing bacteria. Analysis of a rat microbiome experiment identified 47 operational taxonomic units that contributed to fecal butyrate levels. Of these, 31 were positively and 16 negatively associated with butyrate. Consistent with their known role in butyrate metabolism, most of these fell within the Lachnospiraceae and Ruminococcaceae. We conclude that BRACoD provides a more precise and accurate method for determining bacteria associated with a continuous physiological outcome compared to clr-LASSO. It is more sensitive than a generalized clr-SS algorithm, although it has a higher FPR. Its ability to distinguish genuine contributors from correlated bacteria makes it better suited to discriminating bacteria that directly contribute to an outcome. The algorithm corrects for the distortions arising from compositional data making it appropriate for analysis of microbiome data.     

Recombination Facilitates Adaptive Evolution in Rhizobial Soil Bacteria

Homologous recombination is expected to increase natural selection efficacy by decoupling the fate of beneficial and deleterious mutations and by readily creating new combinations of beneficial alleles. Here, we investigate how the proportion of amino acid substitutions fixed by adaptive evolution (α) depends on the recombination rate in bacteria. We analyze 3,086 core protein-coding sequences from 196 genomes belonging to five closely related species of the genus Rhizobium. These genes are found in all species and do not display any signs of introgression between species. We estimate α using the site frequency spectrum (SFS) and divergence data for all pairs of species. We evaluate the impact of recombination within each species by dividing genes into three equally sized recombination classes based on their average level of intragenic linkage disequilibrium. We find that α varies from 0.07 to 0.39 across species and is positively correlated with the level of recombination. This is both due to a higher estimated rate of adaptive evolution and a lower estimated rate of nonadaptive evolution, suggesting that recombination both increases the fixation probability of advantageous variants and decreases the probability of fixation of deleterious variants. Our results demonstrate that homologous recombination facilitates adaptive evolution measured by α in the core genome of prokaryote species in agreement with studies in eukaryotes.     

Angle of Inclination of Tank-Treading Red Cells: Dependence on Shear Rate and Suspending Medium

Red cells suspended in solutions much more viscous than blood plasma assume an almost steady-state orientation when sheared above a threshold value of shear rate. This orientation is a consequence of the motion of the membrane around the red cell called tank-treading. Observed along the undisturbed vorticity of the shear flow, tank-treading red cells appear as slender bodies. Their orientation can be quantified as an angle of inclination (θ) of the major axis with respect to the undisturbed flow direction. We measured θ using solution viscosities (η₀) and shear rates (γ˙) covering one and three orders of magnitude, respectively. At the lower values of η₀, θ was almost independent of γ˙. At the higher values of η₀, θ displayed a maximum at intermediate shear rates. The respective maximal values of θ increased by ∼10° from 10.7 to 104 mPas. After accounting for the absent membrane viscosity in models by using an increased cytoplasmic viscosity, their predictions of θ agree qualitatively with our data. Comparison of the observed variation of θ at constant γ˙ with model results suggests a change in the reference configuration of the shear stiffness of the membrane.     

Viscoelastic Transient of Confined Red Blood Cells

The unique ability of a red blood cell to flow through extremely small microcapillaries depends on the viscoelastic properties of its membrane. Here, we study in vitro the response time upon flow startup exhibited by red blood cells confined into microchannels. We show that the characteristic transient time depends on the imposed flow strength, and that such a dependence gives access to both the effective viscosity and the elastic modulus controlling the temporal response of red cells. A simple theoretical analysis of our experimental data, validated by numerical simulations, further allows us to compute an estimate for the two-dimensional membrane viscosity of red blood cells, ηmem2D ∼ 10⁻⁷ N⋅s⋅m⁻¹. By comparing our results with those from previous studies, we discuss and clarify the origin of the discrepancies found in the literature regarding the determination of ηmem2D, and reconcile seemingly conflicting conclusions from previous works.     

Oxygen uptake,heart rate,perceived exertion,and integrated electromyogram of the lower and upper extremities during level and Nordic walking on a treadmill

The purpose of this study was to characterize responses in oxygen uptake ( V·O2), heart rate (HR), perceived exertion (OMNI scale) and integrated electromyogram (iEMG) readings during incremental Nordic walking (NW) and level walking (LW) on a treadmill. Ten healthy adults (four men, six women), who regularly engaged in physical activity in their daily lives, were enrolled in the study. All subjects were familiar with NW. Each subject began walking at 60 m/min for 3 minutes, with incremental increases of 10 m/min every 2 minutes up to 120 m/min V·O2 , V·E and HR were measured every 30 seconds, and the OMNI scale was used during the final 15 seconds of each exercise. EMG readings were recorded from the triceps brachii, vastus lateralis, biceps femoris, gastrocnemius, and tibialis anterior muscles. V·O2 was significantly higher during NW than during LW, with the exception of the speed of 70 m/min (P < 0.01). V·E and HR were higher during NW than LW at all walking speeds (P < 0.05 to 0.001). OMNI scale of the upper extremities was significantly higher during NW than during LW at all speeds (P < 0.05). Furthermore, the iEMG reading for the VL was lower during NW than during LW at all walking speeds, while the iEMG reading for the BF and GA muscles were significantly lower during NW than LW at some speeds. These data suggest that the use of poles in NW attenuates muscle activity in the lower extremities during the stance and push-off phases, and decreases that of the lower extremities and increase energy expenditure of the upper body and respiratory system at certain walking speeds.