Tree invasions threaten the conservation potential and sustainability of U.S.rangelands O_LIRangelands of the United States provide ecosystem services that sustain biodiversity and rural economies. Native tree encroachment is a long-standing conservation challenge to these landscapes. Still, its impact is often overlooked due to the slow pace of tree invasions and the positive public perception of trees.
C_LIO_LIRecent innovations in satellite remote sensing permit the tracking of tree encroachment and the corresponding impact on herbaceous production. We analyzed tree cover change and herbaceous production across 17 western states from 1990 to 2019.
C_LIO_LIWe show that tree encroachment is a dominant change agent in western U.S. rangelands; tree cover has increased by 50% (77,323 km2) over 30 years, with more than 25% of U.S. rangelands experiencing sustained tree cover expansion. Since 1990, roughly 300 Tg of herbaceous biomass has been lost, totaling $5 billion in foregone revenue to small agricultural producers.
C_LIO_LISynthesis and applications: The impact of tree encroachment to rangeland loss is similar in magnitude to row-crop conversion, another primary threat to U.S. rangelands. Prioritizing conservation efforts to prevent tree encroachment can bolster ecosystem and economic sustainability, particularly among privately-owned lands threatened by land-use conversion.
C_LI
Integrated spatial models foster complementarity between monitoring programs in producing large-scale bottlenose dolphin indicators Over the last decades, large-scale ecological projects have emerged that require collecting ecological data over broad spatial and temporal coverage. Yet, obtaining relevant information about large-scale population dynamics from a single monitoring program is challenging, and often several sources of data, possibly heterogeneous, need to be integrated. In this context, integrated models combine multiple data types into a single analysis to quantify population dynamics of a targeted population. When working at large geographical scales, integrated spatial models have the potential to produce spatialised ecological estimates that would be difficult to obtain if data were analysed separately.

In this paper, we illustrate how spatial integrated modelling offers a relevant framework for conducting ecological inference at large scales. Focusing on the Mediterranean bottlenose dolphins (Tursiops truncatus), we combined 21,464 km of photo-identification boat surveys collecting spatial capture-recapture data with 24,624 km of aerial line-transect following a distance-sampling protocol. We analysed spatial capture-recapture data together with distance-sampling data to estimate abundance and density of bottlenose dolphins. We compared the performances of the distance sampling model and the spatial capture-recapture model fitted independently, to our integrated spatial model.

The outputs of our spatial integrated models inform bottlenose dolphin ecological status in the French Mediterranean Sea and provide ecological indicators that are required for regional scale ecological assessments like the EU Marine Strategy Framework Directive. We argue that integrated spatial models are widely applicable and relevant to conservation research and biodiversity assessment at large spatial scales.
ZAF -- The First Open Source Fully Automated Feeder for Aquatic Facilities In the past few decades, aquatic animals have become popular model organisms in biology, spurring a growing need for establishing aquatic facilities. Zebrafish are widely studied and relatively easy to culture using commercial systems. However, a challenging aspect of maintaining aquatic facilities is animal feeding, which is both time- and resource-consuming. We have developed an open-source fully automatic daily feeding system, Zebrafish Automatic Feeder (ZAF). ZAF is reliable, provides a standardized amount of food to every tank, is cost-efficient and easy to build. The advanced version, ZAF+, allows for the precise control of food distribution as a function of fish density per tank, and has a user-friendly interface. Both ZAF and ZAF+ are adaptable to any laboratory environment and facilitate the implementation of aquatic colonies. Here we provide all blueprints and instructions for building the mechanics, electronics, fluidics, as well as to setup the control software and its user-friendly graphical interface. Importantly, the design is modular and can be scaled to meet different user needs. Furthermore, our results show that ZAF and ZAF+ do not adversely affect zebrafish culture, enabling fully automatic feeding for any aquatic facility. Diverse stem-chondrichthyan oral structures and evidence for an independently acquired acanthodid dentition The teeth of sharks famously form a series of transversely-organised files with conveyor-belt replacement that are borne directly on the jaw cartilages, in contrast to the dermal plate-borne dentition of bony fishes that undergoes site-specific replacement. A major obstacle in understanding how this system evolved is the poorly understood relationships of the earliest chondrichthyans and the profusion of morphologically and terminologically diverse bones, cartilages, splints and whorls that they possess. Here we use tomographic methods to investigate mandibular structures in several early branching  acanthodian-grade stem-chondrichthyans. We show that the dentigerous jaw bones of disparate genera of ischnacanthids are united by a common construction, being growing bones with non-shedding dentition. Mandibular splints, which support the ventro-lateral edge of the Meckels cartilage in some taxa, are formed from dermal bone and may be an acanthodid synapomorphy. We demonstrate that the teeth of Acanthodopsis are borne directly on the mandibular cartilage and that this taxon is deeply nested within an edentulous radiation, representing an unexpected independent origin of teeth. Many or even all of the range of unusual oral structures may be apomorphic, but they should nonetheless be considered when building hypotheses of tooth and jaw evolution, both in chondrichthyans and more broadly. High-throughput variant detection using a color-mixing strategy Many diseases are related to multiple genetic alterations along a single gene. Probing for highly multiple (>10) variants in a single qPCR tube is not possible due to a limited number of fluorescence channels and one variant per channel, so many more tubes are needed. Here, we experimentally validate our novel color-mixing strategy that uses fluorescence combinations as digital color codes to probe multiple variants simultaneously. The color-mixing strategy relies on a simple intra-tube assay that can probe for 15 variants as part of an inter-tube assay that can probe for an exponentially increased number of variants. The color-mixing strategy is achieved using multiplex double-stranded toehold probes modified with fluorophores and quenchers; the probes are designed to be quenched or luminous after binding to wildtype or variant templates. We used the color-mixing strategy to probe for 21 pathogenic mutations in thalassemia and to distinguish between heterozygous and homozygous variants in 6 tubes, with a specificity of 99% and a sensitivity of 94%. To support tuberculosis diagnosis, we used the same strategy to simultaneously probe in Mycobacterium tuberculosis for rifampicin-resistance mutations occurring within one 81-bp region and one 48-bp region in rpoB gene, plus five isoniazid-resistance mutations in inhA and katG genes. The OPAQUE1/DISCORDIA2 myosin XI is required for phragmoplast guidance during asymmetric cell division in maize Asymmetric divisions produce cells with different fates and are critical for development. Here we show that a maize myosin XI protein, OPAQUE1 (O1), is necessary for asymmetric divisions during maize stomatal development. We analyzed stomatal precursor cells prior to and during asymmetric division to determine why o1 mutants have abnormal division planes. Cell polarization and nuclear positioning occur normally in the o1 mutant, and the future site of division is correctly specified. The defect in o1 occurs during late cytokinesis, when the plant-specific phragmoplast -- made of microtubules, actin and other proteins - forms the nascent cell plate. The phragmoplast becomes misguided and does not meet the previously established division site. Initial phragmoplast guidance is correct in o1. However, as phragmoplast expansion continues, phragmoplasts in o1 stomatal precursor cells become misguided and do not meet the cortex at the established division site. To understand how this myosin protein contributes to phragmoplast guidance, we identified O1-interacting proteins. Other myosins, specific actin-binding proteins, and maize kinesins related to the Arabidopsis thaliana division site markers PHRAGMOPLAST ORIENTING KINESINs (POKs) interact with O1. We propose that different myosins are important at multiple steps of phragmoplast expansion, and the O1 actin motor and POK-like microtubule motors work together to ensure correct late-stage phragmoplast guidance. Cell-specific Bioorthogonal Tagging of Glycoproteins Altered glycosylation is an undisputed corollary of cancer development. Understanding these alterations is paramount but hampered by limitations underlying cellular model systems. For instance, the intricate interactions between tumour and host cannot be adequately recapitulated in monoculture of tumour-derived cell lines. More complex co-culture models usually rely on sorting procedures for proteome analyses and rarely capture the details of protein glycosylation. Here, we report a strategy termed Bio-Orthogonal Cell line-specific Tagging of Glycoproteins (BOCTAG). Cells are equipped by transfection with an artificial biosynthetic pathway that transforms bioorthogonally tagged sugars into the corresponding nucleotide-sugars. Only transfected cells incorporate bioorthogonal tags into glycoproteins in the presence of non-transfected cells. We employ BOCTAG as an imaging technique and to annotate cell-specific glycosylation sites in mass spectrometry-glycoproteomics. We demonstrate application in co-culture and mouse models, allowing for profiling of the glycoproteome as an important modulator of cellular function. Fatal COVID-19 and non-COVID-19 Acute Respiratory Distress Syndrome is Associated with Incomplete Alveolar Type 1 Epithelial Cell Differentiation from the Transitional State Without Fibrosis ARDS due to COVID-19 and other etiologies results from injury to the alveolar epithelial cell (AEC) barrier resulting in noncardiogenic pulmonary edema, which causes acute respiratory failure; clinical recovery requires epithelial regeneration. During physiologic regeneration in mice, AEC2s proliferate, exit the cell cycle, and transiently assume a transitional state before differentiating into AEC1s; persistence of the transitional state is associated with pulmonary fibrosis in humans. It is unknown whether transitional cells emerge and differentiate into AEC1s without fibrosis in human ARDS and why transitional cells differentiate into AEC1s during physiologic regeneration but persist in fibrosis. We hypothesized that incomplete but ongoing AEC1 differentiation from transitional cells without fibrosis may underlie persistent barrier permeability and fatal acute respiratory failure in ARDS. Immunostaining of postmortem ARDS lungs revealed abundant transitional cells in organized monolayers on alveolar septa without fibrosis. They were typically cuboidal or partially spread, sometimes flat, and occasionally expressed AEC1 markers. Immunostaining and/or interrogation of scRNAseq datasets revealed that transitional cells in mouse models of physiologic regeneration, ARDS, and fibrosis express markers of cell cycle exit but only in fibrosis express a specific senescence marker. Thus, in severe, fatal early ARDS, AEC1 differentiation from transitional cells is incomplete, underlying persistent barrier permeability and respiratory failure, but ongoing without fibrosis; senescence of transitional cells may be associated with pulmonary fibrosis. DNA supercoiling-mediated collective behavior of co-transcribing RNA polymerases Multiple RNA polymerases (RNAPs) transcribing a gene have been known to exhibit collective group behavior, causing the transcription elongation rate to increase with the rate of transcription initiation. Such behavior has long been believed to be driven by a physical interaction or "push" between closely spaced RNAPs. However, recent studies have posited that RNAPs separated by longer distances may cooperate via the DNA segment under transcription. Here, we present a theoretical model incorporating the mechanical coupling between RNAP translocation and the torsional response of supercoiled DNA. Using stochastic simulations, we demonstrate long-range cooperation between co-transcribing RNAPs mediated by DNA supercoiling. We find that inhibiting transcription initiation can slow down the already recruited RNAPs, in agreement with recent experimental observations, and predict that the average transcription elongation rate varies non-monotonically with the rate of transcription initiation. We further show that while RNAPs transcribing neighboring genes oriented in tandem can cooperate, those transcribing genes in divergent or convergent orientations can act antagonistically, and that such behavior holds over a large range of intergenic separations. Our model makes testable predictions, revealing how the mechanical interplay between RNAPs and the DNA they transcribe can govern a key cellular process. Receptor-binding domain recombinant protein on alum-CpG induces broad protection against SARS-CoV-2 variants of concern We conducted preclinical studies in mice using a yeast-produced SARS-CoV-2 RBD subunit vaccine candidate formulated with aluminum hydroxide (alum) and CpG deoxynucleotides. This formulation is equivalent to the Corbevax vaccine that recently received emergency use authorization by the Drugs Controller General of India. We compared the immune response of mice vaccinated with RBD/alum to mice vaccinated with RBD/alum+CpG. We also evaluated mice immunized with RBD/alum+CpG and boosted with RBD/alum. Mice were immunized twice intramuscularly at a 21-day interval. Compared to two doses of the /alum formulation, the RBD/alum+CpG vaccine induced a stronger and more balanced Th1/Th2 cellular immune response, with high levels of neutralizing antibodies against the original Wuhan isolate of SARS-CoV-2 as well as the B.1.1.7 (Alpha), B. 1.351 (Beta), B. 1.617.2 and (Delta) variants. Neutralizing antibody titers against the B.1.1.529 (BA.1, Omicron) variant exceeded those in human convalescent plasma after Wuhan infection but were lower than against the other variants. Interestingly, the second dose did not benefit from the addition of CpG, possibly allowing dose-sparing of the adjuvant in the future. The data reported here reinforces that the RBD/alum+CpG vaccine formulation is suitable for inducing broadly neutralizing antibodies against SARS-CoV-2 including variants of concern. Stimulated oxygen transport in tissue by magnetic needle acupuncture AimTo understand physiology of magnetic needle acupuncture by investigating O2 transport in tissue during needle intervention.

MethodsO2 transport in tissue is modeled by utilizing COMSOL with magnetic needle inserted into muscle tissue in a 2D porous media. The damaged tissue has been mimicked by an extracted tissue block with 1st order O2 consumption rate. The convection-diffusion O2 transport in the damaged tissue has been further evaluated by varying magnetic flux density B0 of the needle (0-1 T), myoglobin concentration (0-1 mM), O2 tension (5-100 Torr), O2 consumption rates, tissue permeability (10-12-10-6 m2) and porosity (0.1-0.9).

Results1) Active O2 transport carried by interstitial flow is enhanced with the intervention of a magnetic needle by generating a high gradient magnetic field around the tip, which exerts a strong force (104 N/m3) on the diamagnetic interstitial fluid to accelerate the flow. 2) This interstitial flow can reach 30 m/s at B0 = 1T and strongly correlates to the needle B0 and tissue permeability. 3) The needle stirs the interstitial flow can pump O2 flux by 1-2 orders of magnitude compared to that without magnetic field. 4) The enhancement of active O2 transport by magnetic needle is site-specific to the tissue in the vicinity of the tip. This enhancement is more effective in edema condition with a high tissue permeability (>10-9 m2).

ConclusionsThe dramatic enhancement of O2 transport to restore the O2 mitochondria metabolism for dysfunctional muscle tissues is the fundamental physiological mechanism of magnetic needle acupuncture. Modeling miRNA-driven post-transcriptional regulation by using exon-intron split analysis (EISA) in pigs The contribution of microRNAs (miRNAs) to mRNA regulation has often been explored by post hoc selection of downregulated genes and determining whether they harbor binding sites for miRNAs of interest. This approach, however, does not discriminate whether these mRNAs are also downregulated at the transcriptional level. Here, we have characterized the transcriptional and post-transcriptional changes of mRNA expression in two porcine tissues: gluteus medius muscle of fasted and fed Duroc gilts and adipose tissue of lean and obese Duroc-Gottingen minipigs. Exon-intron split analysis (EISA) of RNA-seq data allowed us to identify downregulated mRNAs with high post-transcriptional signals in fed or obese states, and we assessed whether they harbor binding sites for upregulated miRNAs in any of these two physiological states. We found 26 downregulated mRNAs with high post-transcriptional signals in the muscle of fed gilts and 21 of these were predicted targets of upregulated miRNAs also in the fed state. For adipose tissue, 44 downregulated mRNAs in obese minipigs displayed high post-transcriptional signals, and 25 of these were predicted targets of miRNAs upregulated in the obese state. These results suggest that the contribution of miRNAs to mRNA repression is more prominent in the skeletal muscle system. Finally, we identified several genes that may play relevant roles in the energy homeostasis of the pig skeletal muscle (DKK2 and PDK4) and adipose (SESN3 and ESRRG) tissues. By differentiating transcriptional from post-transcriptional changes in mRNA expression, EISA provides a valuable view about the regulation of gene expression, complementary to canonical differential expression analyses. Different contribution of the monkey prefrontal and premotor dorsal cortex in decision making during a transitive inference task Several studies have reported similar neural modulations between brain areas of the frontal cortex, such as the dorsolateral prefrontal (DLPFC) and the premotor dorsal (PMd) cortex, in tasks requiring encoding of the abstract rules for selecting the proper action. Here we compared the neuronal modulation of the DLPFC and PMd of monkeys trained to choose the higher rank from a pair of abstract images (target item), selected from an arbitrarily rank-ordered set (A>B>C>D>E>F) in the context of a transitive inference task. Once acquired by trial-and-error, the ordinal relationship between pairs of adjacent images (i.e., A>B; B>C; C>D; D>E; E>F), monkeys were tested in indicating the ordinal relation between items of the list not paired during learning. During these decisions, we observed that the choice accuracy increased and the reaction time decreased as the rank difference between the compared items enhanced. This result is in line with the hypothesis that after learning, the monkeys built an abstract mental representation of the ranked items, where rank comparisons correspond to the items position comparison on this representation. In both brain areas, we observed higher neuronal activity when the target item appeared in a specific location on the screen with respect to the opposite position and that this difference was particularly enhanced at lower degrees of difficulty. By comparing the time evolution of the activity of the two areas, we observed that the neural encoding of target item spatial position occurred earlier in the DLPFC than in the PMd. Pericyte-to-endothelial cell signaling via vitronectin-integrin regulates blood-CNS barrier Endothelial cells of blood vessels of the central nervous system (CNS) constitute blood-CNS barriers. Barrier properties are not intrinsic to these cells; rather they are induced and maintained by CNS microenvironment. Notably, the abluminal surface of CNS capillaries are ensheathed by pericytes and astrocytes. However, extrinsic factors from these perivascular cells that regulate barrier integrity are largely unknown. Here, we establish vitronectin, an extracellular-matrix protein secreted by CNS pericytes, as a regulator of blood-CNS barrier function via interactions with its integrin receptor, 5 in endothelial cells. Genetic ablation of vitronectin or mutating vitronectin to prevent integrin binding as well as endothelial-specific deletion of integrin 5 causes barrier leakage. Furthermore, vitronectin-integrin 5 signaling maintains barrier integrity by actively inhibiting transcytosis in endothelial cells. These results demonstrate that signaling from perivascular cells to endothelial cells via ligand-receptor interactions is a key mechanism to regulate barrier permeability. Towards a mathematical understanding of colonization resistance in multispecies microbial communities Microbial community composition and dynamics are key to health and disease. Explaining the forces generating and shaping diversity in the microbial consortia making up our bodys defenses is a major aim of current research in microbiology. For this, tractable models are needed, that bridge the gap between observations of patterns and underlying mechanisms. While most microbial dynamics models are based on the Lotka-Volterra framework, we still do not have an analytic quantity for colonization resistance, by which a microbial systems fitness as a whole can be understood. Here, inspired by an epidemiological perspective, we propose a rather general modeling framework whereby colonization resistance can be clearly mathematically defined and studied. In our model, N similar species interact with each other through a co-colonization interaction network encompassing pairwise competition and cooperation, abstractly mirroring how organisms effectively modify their micro-scale environment in relation to others. This formulation relies on a generic notion of shared resources between members of a consortium, yielding explicit frequency-dependent dynamics among N species, in the form of a replicator equation, and offering a precise definition of colonization resistance. We demonstrate that colonization resistance arises and evolves naturally in a multispecies system as a collective quadratic term in a replicator equation, describing dynamic mean invasion fitness. Each pairwise invasion growth rate between two ecological partners, [Formula], is derived explicitly from species asymmetries and mean traits. This makes the systemic colonization resistance [Formula] also an emergent function of global mean-field parameters and trait variation architecture, weighted by the evolving relative abundances among species. In particular, if the underlying invasion fitness matrix {Lambda} displays species-specific  invasiveness or  invasibility structure, colonization resistance will be insensitive to mean micro-scale cooperation or competition. However, in general, colonization resistance depends on and may undergo critical transitions with changes in mean  environment, e.g. cooperation and growth level in a community. We illustrate several key links between our proposed measure of colonization resistance and invader success, including sensitivity to timing, and to the intrinsic pairwise invasion architecture of the resident community. Our simulations reveal that symmetric and invader-driven mutual invasion among resident species tend to maximize systemic colonization resistance to outsiders, when compared to resident-driven, anti-symmetric, almost anti-symmetric and random {Lambda} structures. We contend this modeling approach is a powerful new avenue to study, test and validate interaction networks and invasion topologies in diverse microbial consortia, and quantify analytically their role in colonization resistance, system function, and invasibility. Predictive uncertainty in mechanistic models of cellular processes calibrated to experimental data Mathematical models are often used to explore network-driven cellular processes from a systems perspective. However, a dearth of quantitative data suitable for model calibration leads to models with parameter unidentifiability and questionable predictive power. Here we introduce a Bayesian and Machine-Learning based Measurement Model approach to explore how quantitative and non-quantitative data constrain models of apoptosis execution within a missing data context. We find two orders of magnitude more ordinal (e.g. immunoblot) data are necessary to achieve accuracy comparable to quantitative (e.g. fluorescence) data. Notably, ordinal and nominal (e.g. immunostain) non-quantitative data synergize to reduce model uncertainty and improve accuracy. Further, model prediction accuracy and certainty strongly depend on rigorous data-driven formulations of the measurement, and the size and make-up of the datasets. Finally, we demonstrate the potential of a data-driven Measurement Model approach to identify model features that could lead to informative experimental measurements and improve model predictive power. An Immune Cell Atlas Reveals Dynamic COVID-19 Specific Neutrophil Programming Amenable to Dexamethasone Therapy SARS-CoV-2 is a novel coronavirus that causes acute respiratory distress syndrome (ARDS), death and long-term sequelae. Innate immune cells are critical for host defense but are also the primary drivers of ARDS. The relationships between innate cellular responses in ARDS resulting from COVID-19 compared to other causes of ARDS, such as bacterial sepsis is unclear. Moreover, the beneficial effects of dexamethasone therapy during severe COVID-19 remain speculative, but understanding the mechanistic effects could improve evidence-based therapeutic interventions. To interrogate these relationships, we developed an scRNA-Seq and plasma proteomics atlas (biernaskielab.ca/COVID_neutrophil). We discovered that compared to bacterial ARDS, COVID-19 was associated with distinct neutrophil polarization characterized by either interferon (IFN) or prostaglandin (PG) active states. Neutrophils from bacterial ARDS had higher expression of antibacterial molecules such as PLAC8 and CD83. Dexamethasone therapy in COVID patients rapidly altered the IFNactive state, downregulated interferon responsive genes, and activated IL1R2+ve neutrophils. Dexamethasone also induced the emergence of immature neutrophils expressing immunosuppressive molecules ARG1 and ANXA1, which were not present in healthy controls. Moreover, dexamethasone remodeled global cellular interactions by changing neutrophils from information receivers into information providers. Importantly, male patients had higher proportions of IFNactive neutrophils, a greater degree of steroid-induced immature neutrophil expansion, and increased mortality benefit compared to females in the dexamethasone era. Indeed, the highest proportion of IFNactive neutrophils was associated with mortality. These results define neutrophil states unique to COVID-19 when contextualized to other life-threatening infections, thereby enhancing the relevance of our findings at the bedside. Furthermore, the molecular benefits of dexamethasone therapy are also defined, and the identified pathways and plasma proteins can now be targeted to develop improved therapeutics. Enhanced nociceptor sprouting in a rabbit model of cerebral palsy Spastic cerebral palsy (CP) is a movement disorder marked by hypertonia and hyperreflexia, and the most prevalent comorbidity is pain. Since spinal nociceptive afferents contribute to both the sensation of painful stimuli as well as reflex circuits involved in movement, we investigated the relationship between prenatal hypoxia-ischemia (HI) injury which can cause CP, and possible changes in spinal nociceptive circuitry. To do this, we examined nociceptive afferents and mechanical and thermal sensitivity of New Zealand White rabbit kits after prenatal HI or a sham surgical procedure. As described previously, a range of motor deficits similar to spastic CP was observed in kits born naturally after HI (40 minutes at ~70-80% gestation). We found that HI caused an expansion of peptidergic afferents (marked by expression of calcitonin gene-related peptide; CGRP) in both the superficial and deep dorsal horn at postnatal day (P)5. Non-peptidergic nociceptive afferent arborization (labeled by isolectin B4; IB4) was unaltered in HI kits but overlap of the two populations (peptidergic and non-peptidergic nociceptors) was increased by HI. Density of glial fibrillary acidic protein (GFAP) was unchanged within spinal white matter regions important in nociceptive transmission at P5. We found that mechanical and thermal allodynia was present in HI kits even in the absence of motor deficits. These findings suggest that prenatal HI injury impacts spinal sensory pathways in addition to the more well-established disruptions to descending motor circuits. In conclusion, changes to spinal nociceptive circuitry could disrupt spinal reflexes and contribute to pain experienced by individuals with CP. Changes in Behavior and Salivary Serum Amyloid A in cows with Subclinical Mastitis The aim of this study was to identify detailed changes in behavior, and in salivary serum amyloid A (SAA), associated with subclinical mastitis. This included standard sickness behaviors (which form part of an adaptive response to conserve energy, minimize heat loss and facilitate recovery following infection and inflammation) and lesser-studied behaviors, that are biologically predicted to change at low-level infection, and therefore particularly relevant for detecting subclinical mastitis (here labelled  luxury behaviors). SAA is a physiological marker of systemic inflammation, with levels in milk and serum already known to increase during subclinical mastitis. We investigated whether the same was true for SAA in cow saliva. Data were collected for n = 34 commercial barn-housed dairy cows, tested in 17 matched-pairs. Each pair comprised a cow with subclinical mastitis (SCM) and a healthy control (CTRL), identified using somatic cell count (SCC) (SCM: SCC >200 x1000 cells/ml; CTRL: SCC <100 x1000 cells/ml). SCM cows were selected for study ad-hoc, at which point they were paired with a CTRL cow, based upon parity and calving date; consequently, the full data set was accrued over several months. Data were collected for each pair over 3 days: SCC (l4:00-l5:00h) Day 1; behavior (24h from 00:00h) Day 2; salivary serum amyloid-A (SAA) Day 3. We report, for the first time, that an increase in salivary SAA occurs during subclinical mastitis; SAA was higher in SCM cows and demonstrated a positive (weak) correlation with SCC. The behavioral comparisons revealed that SCM cows had reductions in activity (behavioral transitions and distance moved), social exploration, social reactivity (here: likelihood to move away/be displaced following receipt of agonism), performance of social grooming and head butts, and the receipt of agonistic non-contact challenges. In addition, SCM cows received more head swipes, and spent a greater proportion of time lying with their head on their flank than CTRL cows. SCM cows also displayed an altered feeding pattern; they spent a greater proportion of feeding time in direct contact with two conspecifics, and a lower proportion of feeding time at self-locking feed barriers, than CTRL cows. Behavioral measures were found to correlate, albeit loosely, with serum SAA in a direction consistent with predictions for sickness behavior. These included positive correlations with both lying duration and the receipt of all agonistic behavior, and negative correlations with feeding, drinking, the performance of all social and all agonistic behavior, and social reactivity. We conclude that changes in salivary SAA, social behavior, and activity offer potential in the detection of subclinical mastitis and recommend further investigation to substantiate and refine our findings. Critical re-evaluation of the slope factor of the operational model of agonism Background and PurposeAlthough being a relative term, agonist efficacy is a cornerstone in the proper assessment of agonist selectivity and signalling bias. The operational model of agonism (OMA) has become successful in the determination of agonist efficacies and ranking them. In 1985, Black et al. introduced the slope factor to the OMA to make it more flexible and allow for fitting steep as well as flat concentration-response curves.

Experimental ApproachFunctional analysis of OMA demonstrates that the slope factor implemented by Black et al. affects relationships among parameters of the OMA, leading to wrong estimates of operational efficacy and affinity of an agonist. In contrast, functional analysis of OMA with Hill coefficient as a slope factor has not this flaw.

Key ResultsFitting of the OMA with Black et al. slope factor to concentration-response curves of several models (homotropic allosteric modulation, substrate inhibition and non-competitive auto-inhibition) gives wrong estimates of operational efficacy and affinity. In contrast, fitting of the OMA modified by the Hill coefficient to concentration-response curves of these several models gives correct estimates of operational efficacy and affinity.

Conclusion and ImplicationsThe OMA modified by the Hill coefficient can be used for correct ranking of agonism and subsequent analysis, like quantification of signalling bias. TeamTree analysis: a new approach to evaluate scientific production Advances in science and technology depend on the work of research teams and the publication of results through peer-reviewed articles representing a growing socio-economic resource. Current methods to mine the scientific literature regarding a field of interest focus on content, but the workforce credited by authorship remains largely unexplored, and appropriate measures of scientific production are debated. Here, a new bibliometric approach named TeamTree analysis is introduced that visualizes the development and composition of the workforce driving a field. A new citation-independent measure that scales with the H index estimates impact based on publication record, genealogical ties and collaborative connections. This author-centered approach complements existing tools to mine the scientific literature and to evaluate research across disciplines. Integrative analysis reveals aged clonal B cells, microenvironment and c-Myc activation in the origin of age-related lymphoma While cancer is an age-related disease, most studies focus on genetically engineered younger mouse models. Here, we uncover how cancer develops as a consequence of the naturally aged immune system in mice. B-cell lymphoma frequently occurs in aged mice and is associated with increased cell size, splenomegaly, and a novel clonal B-cell population. Age-emergent B cells clonally expand outside of germinal centers driven by somatic mutations, activated c-Myc and hypermethylated promoters, and both genetically and epigenetically recapitulate human follicular and diffuse-large B-cell lymphomas. Mechanistically, mouse cancerous B cells originate from age-associated B cells, which are atypical memory B cells. Age-associated B cells secrete a spectrum of proinflammatory cytokines and activate paracrinally the expression of c-Myc in surrounding B cells. Although clonal B cells are a product of an aging microenvironment, they evolve being self-sufficient and support malignancy when transferred into young mice. Inhibition of mTOR and c-Myc attenuates premalignant changes in B cells during aging and emerges as a therapeutic strategy to delay the onset of age-related lymphoma. Together, we uncover how aging generates cancerous B cells, characterize a model that captures the origin of spontaneous cancer during aging and identify interventions that may postpone age-associated lymphoma. Non-Newtonian Blood Rheology Impacts Left Atrial Stasis in Patient-Specific Simulations The lack of mechanically effective contraction of the left atrium (LA) during atrial fibrillation (AF) disturbs blood flow, increasing the risk of thrombosis and ischemic stroke. Thrombosis is most likely in the left atrial appendage (LAA), a small narrow sac where blood is prone to stagnate. Slow flow promotes the formation of erythrocyte aggregates in the LAA, also known as rouleaux, causing viscosity gradients that are usually disregarded in patient-specific simulations. To evaluate these non-Newtonian effects, we built atrial models derived from 4D computed tomography scans of patients and carried out computational fluid dynamics simulations using the Carreau-Yasuda constitutive relation. We examined six patients, three of whom had AF and LAA thrombosis or a history of transient ischemic attacks (TIAs). We modeled the effects of hematocrit and rouleaux formation kinetics by varying the parameterization of the Carreau-Yasuda relation and modulating non-Newtonian viscosity changes based on residence time. Comparing non-Newtonian and Newtonian simulations indicates that slow flow in the LAA increases blood viscosity, altering secondary swirling flows and intensifying blood stasis. While some of these effects can be subtle when examined using instantaneous metrics like shear rate or kinetic energy, they are manifested in the blood residence time, which accumulates over multiple heartbeats. Our data also reveal that LAA blood stasis worsens when hematocrit increases, offering a potential new mechanism for the clinically reported correlation between hematocrit and stroke incidence. In summary, we submit that hematocrit-dependent non-Newtonian blood rheology should be considered in calculating patient-specific blood stasis indices by computational fluid dynamics. Integrated quality control of allele-specific copy numbers, mutations and tumour purity from cancer whole genome sequencing assays Cancer genomes contain thousands of somatic point mutations, chromosome copy alterations and more complex structural variants, which contribute to tumour growth and therapy response. Whole genome sequencing is a well established approach for somatic variant identification, but its broad application comes with complications, particularly in how proposed calls are quality assessed. To address this issue, we present CNAqc, a quantitative framework to quality control somatic mutations and allele-specific copy numbers, both in clonal and subclonal settings while accounting for variations in tumour purity, as commonly seen in bulk sampling. We test the model via extensive simulations, validate it using low-pass single-cell data, and apply it to 2778 single-sample PCAWG whole-genomes, 10 in-house multi-region whole-genomes and 48 TCGA whole-exomes. CNAqc is compatible with common bioinformatic pipelines and designed to support automated parameterization processes that are crucial in the era of large-scale whole genome sequencing. Shorter effective lifespan in laboratory populations of D. melanogaster can reduce sexual selection. The role of sexual selection in mediating levels of sexual conflict has been demonstrated in many experimental evolution studies on Drosophila spp. where competition among males for mating was the target of selection. Sexual selection has also been shown to affect the evolution of life-histories. However, the influence of divergent life-histories on reproductive strategies and, therefore, sexual selection and possibly sexual conflict, has been less well studied. We examined D. melanogaster populations selected for a short development time and early age at reproduction for changes in reproductive behaviour and traits that are proxies of sexual selection. We report a large reduction in reproductive competition experienced by the males of these populations, compared to ancestral populations that are not consciously selected for rapid development or early reproduction, potentially leading to reduced sexual selection. We show that rapidly developing and early reproducing populations have very low levels of mating in their lifetime (females are more or less monandrous), low courtship levels, shorter copulation duration, and longer time from eclosion to first mating, compared to the controls. These results are discussed in the context of the previously demonstrated reduction of inter-locus sexual conflict in these populations. We show that life-history strategies might have a large and significant impact on sexual selection, with each influencing the other and contributing to the complexities of adaptation.

Significance statementSexual conflict, often manifested as an arms-race between males and females trying to enhance their own reproductive success at some cost to the other, is of great evolutionary interest because it can maintain genetic variation in populations, prevent the independent optimization of male and female traits, and also promote speciation. Sexual selection, or variation in mating success, is well known to affect levels of sexual conflict. However, it is not so clear whether, and how, the regular evolution of life-histories also affects sexual selection. Here, we show that life-history evolution in fruit fly populations selected for traits not directly related to sexual conflict might, nevertheless, mediate the possible evolution of altered sexual conflict levels through effects on sexual selection. Populations that evolved to develop to adulthood fast, and reproduce relatively early in life, are shown to potentially experience less sexual selection, which can explain the low sexual conflict levels earlier observed in them. A thermodynamic model for water activity and redox potential in evolution and development Reactions involving water and oxygen are basic features of geological and biological processes. To understand how life interacts with its environment requires monitoring interactions with H2O and O2 at timescales relevant to not only organismal growth but also billions of years of geobiological evolution. Here, chemical transformations intrinsic to evolution and development were characterized by analyzing data from recent phylostratigraphic and proteomic studies. This two-stage analysis involves obtaining chemical metrics (carbon oxidation state and stoichiometric hydration state) from the elemental compositions of proteins followed by modeling the relative stabilities of target proteins against a proteomic background to infer thermodynamic parameters (oxygen fugacity, water activity, and virtual redox potential (Eh)). The main results of this study are a rise in carbon oxidation state of proteins spanning the time of the Great Oxidation Event, a rise in virtual redox potential that coincides with the likely emergence of aerobic metabolism, a rise in carbon oxidation state of proteins inferred from the transcriptome in late stages of Bacillus subtilis biofilm growth, and a drop in stoichiometric hydration state of the fruit fly developmental proteome at the same time as a drop in organismal water content. Stoichiometric hydration state also decreases for proteins with more recent gene ages and through stages of biofilm development, leading to predictions of higher hydration potentials at earlier time points. By building chemical representations of protein evolution and developmental proteomes, exciting new types of geobiochemical proxies can be developed.