Linking individual differences in human primary visual cortex to contrast sensitivity around the visual field A central question in neuroscience is how the organization of cortical maps relates to perception, for which human primary visual cortex (V1) is an ideal model system. V1 nonuniformly samples the retinal image, with greater cortical magnification (surface area per degree of visual field) at the fovea than periphery and at the horizontal than vertical meridian. Moreover, the size and cortical magnification of V1 varies greatly across individuals. Here, we used fMRI and psychophysics in the same observers to quantify individual differences in V1 cortical magnification and contrast sensitivity at the four polar angle meridians. Across observers, the overall size of V1 and localized cortical magnification positively correlated with contrast sensitivity. Moreover, greater cortical magnification and higher contrast sensitivity at the horizontal than the vertical meridian were strongly correlated. These data reveal a link between cortical anatomy and visual perception at the level of individual observer and stimulus location.
Hierachical Resampling for Bagging in Multi-Study Prediction with Applications to Human Neurochemical Sensing We propose the "study strap ensemble," which combines advantages of two common approaches to fitting prediction models when multiple training datasets ("studies") are available: pooling studies and fitting one model versus averaging predictions from multiple models each fit to individual studies. The study strap ensemble fits models to bootstrapped datasets, or "pseudo-studies." These are generated by resampling from multiple studies with a hierarchical resampling scheme that generalizes the randomized cluster bootstrap. The study strap is controlled by a tuning parameter that determines the proportion of observations to draw from each study. When the parameter is set to its lowest value, each pseudo-study is resampled from only a single study. When it is high, the study strap ignores the multi-study structure and generates pseudo-studies by merging the datasets and drawing observations like a standard bootstrap. We empirically show the optimal tuning value often lies in between, and prove that special cases of the study strap draw the merged dataset and the set of original studies as pseudo-studies. We extend the study strap approach with an ensemble weighting scheme that utilizes information in the distribution of the covariates of the test dataset.

Our work is motivated by neuroscience experiments using real-time neurochemical sensing during awake behavior in humans. Current techniques to perform this kind of research require measurements from an electrode placed in the brain during awake neurosurgery and rely on prediction models to estimate neurotransmitter concentrations from the electrical measurements recorded by the electrode. These models are trained by combining multiple datasets that are collected in vitro under heterogeneous conditions in order to promote accuracy of the models when applied to data collected in the brain. A prevailing challenge is deciding how to combine studies or ensemble models trained on different studies to enhance model generalizability.

Our methods produce marked improvements in simulations and in this application. All methods are available in the studyStrap CRAN package.
Bone mineral density modeling via random field: normality, stationarity, sex and age dependence Background and ObjectiveCapturing the population variability of bone properties is of paramount importance to biomedical engineering. The aim of the present paper is to describe variability and correlations in bone mineral density with a spatial random field inferred from routine computed tomography data.

MethodsRandom fields were simulated by transforming pairwise uncorrelated Gaussian random variables into correlated variables through the spectral decomposition of an age-detrended correlation matrix. The validity of the random field model was demonstrated in the spatiotemporal analysis of bone mineral density. The similarity between the computed tomography samples and those generated via random fields was analyzed with the energy distance metric.

ResultsThe random field of bone mineral density was found to be approximately Gaussian/slightly left-skewed/strongly right-skewed at various locations. However, average bone density could be simulated well with the proposed Gaussian random field for which the energy distance, i.e., a measure that quantifies discrepancies between two distribution functions, is convergent with respect to the number of correlation eigenpairs.

ConclusionsThe proposed random field model allows the enhancement of computational biomechanical models with variability in bone mineral density, which could increase the usability of the model and provides a step forward in in-silico medicine. Spatial phylogenetics of Japanese ferns: Patterns, processes, and implications for conservation PremiseBiodiversity is often only measured with species richness. However, this ignores evolutionary history and is not sufficient for making conservation decisions. Here, we characterize multiple facets and drivers of biodiversity to understand how these relate to bioregions and conservation status in the ferns of Japan.

MethodsWe compiled a community dataset of 1,239 20 km x 20 km grid-cells including 672 taxa based on > 300,000 specimen records. We combined this with a phylogeny and functional traits to analyze taxonomic, phylogenetic, and functional diversity, and modeled biodiversity metrics in response to environmental factors and reproductive mode. Hierarchical clustering was used to delimit bioregions. Conservation status and threats were assessed by comparing the overlap of significantly diverse grid-cells with conservation zones and range maps of native Japanese deer.

ResultsTaxonomic richness was highest at mid-latitudes. Phylogenetic and functional diversity and phylogenetic endemism were highest in small southern islands. Relative phylogenetic and functional diversity were high at high and low latitudes, and low at mid-latitudes. Grid-cells were grouped into three (phylogenetic) or four (taxonomic) major bioregions. Temperature and apomixis were identified as drivers of biodiversity patterns. Conservation status was generally high for grid-cells with significantly high biodiversity, but the threat due to herbivory by deer was greater for taxonomic richness than other metrics.

ConclusionsOur integrative approach reveals previously undetected patterns and drivers of biodiversity in the ferns of Japan. Future conservation efforts should recognize that threats can vary by biodiversity metric and consider multiple multiple metrics when establishing conservation priorities. Analysis of science journalism reveals gender and regional disparities in coverage Science journalism is a critical way in which the public can remain informed and benefit from new scientific findings. Such journalism also shapes the publics view of the current state of scientific findings and legitimizes experts. Those covering science can only cite and quote a limited number of sources. Sources may be identified by the journalists research or by recommendations by other scientists. In both cases, biases may influence who is identified and ultimately included as an expert. To analyze possible biases in science journalism, we analyzed 22,001 non-research articles published by Nature. We chose to analyze Nature non-research articles since its research articles provide a natural comparator. Our analysis considered two possible sources of disparity: gender and name origin. To explore these sources of disparity, we extracted cited authors names as well as extracted names of quoted speakers. While citations and quotations within a piece do not reflect the entire information-gathering process, they can provide insight into the demographics of visible sources. We then used the extracted names to predict gender and name origin of the cited authors and speakers.

In order to appropriately quantify the level of difference, we must identify a suitable reference set for comparison. We chose first and last authors within primary research articles in Nature and a subset of Springer Nature articles in the same time period as our comparator. In our analysis, we found a skew towards male quotation in Nature science journalism-related articles. However, quotation is trending toward equal representation at a faster rate than first and last authorship in academic publishing. Interestingly, we found that the gender disparity in Nature quotes was column-dependent, with the "Career Features" column reaching gender parity. Our name origin analysis found a significant over-representation of names with predicted Celtic/English origin and under-representation of names with a predicted East Asian origin. This finding was observed both in extracted quotes and journal citations, but dampened in citations. Identifying and correcting for misspecifications in GWAS summary statistics and polygenic scores Publicly available genome-wide association studies (GWAS) summary statistics exhibit uneven quality, which can impact the validity of follow-up analyses. First, we present an overview of possible misspecifications that come with GWAS summary statistics. Then, in both simulations and real data analyses, we show that additional information such as imputation INFO scores, allele frequencies, and per-variant sample sizes in GWAS summary statistics can be used to detect possible issues and correct for misspecifications in the GWAS summary statistics. One important motivation for us is to improve the predictive performance of polygenic scores built from these summary statistics. Unfortunately, due to the lack of reporting standards for GWAS summary statistics, this additional information is not systematically reported. We also show that using well-matched LD references can improve model fit and translate into more accurate prediction. Finally, we discuss how to make polygenic score methods such as lassosum and LDpred2 more robust to these misspecifications to improve their predictive power. A novel Fiji/ImageJ plugin for the rapid analysis of blebbing cells When confined, cells have recently been shown to undergo a phenotypic switch to what has been termed, fast amoeboid (leader bleb-based) migration. However, as this is a nascent area of research, few tools are available for the rapid analysis of cell behavior. Here, we demonstrate that a novel Fiji/ImageJ-based plugin, Analyze_Blebs, can be used to quickly obtain cell migration parameters and morphometrics from time lapse images. As validation, we show that Analyze_Blebs can detect significant differences in cell migration and morphometrics, such as the largest bleb size, upon introducing different live markers of F-actin, including F-tractin and LifeAct tagged with green and red fluorescent proteins. We also demonstrate, using flow cytometry, that live markers increase total levels of F-actin. Furthermore, that F-tractin increases cell stiffness, which was found to correlate with a decrease in migration, thus reaffirming the importance of cell mechanics as a determinant of Leader Bleb-Based Migration (LBBM).

Highlight summaryA new plugin, Analyze_Blebs, enables the rapid analysis of cell migration and morphometrics of fast amoeboid cells. Morphometrics combined with cell stiffness are found to be predictive of confined, fast amoeboid (leader bleb-based), migration. Ultrahigh-Throughput Ambient MS: Direct Analysis at 22 Samples per Second by Infrared Matrix-Assisted Laser Desorption Electrospray Ionization Mass Spectrometry Infrared Matrix-Assisted Laser Desorption Electrospray Ionization (IR-MALDESI) mass spectrometry is an ambient-direct sampling method being developed for high-throughput, label-free, biochemical screening of large-scale compound libraries. Here, we report the development of an ultrahigh-throughput continuous motion IR-MALDESI sampling approach capable of acquiring data at rates up to 22.7 samples per second in a 384-well microtiter plate. At top speed, less than 1% analyte carryover is observed from well-to-well and signal intensity relative standard deviations (RSD) of 11.5% and 20.9% for 3 M 1-hydroxymidazolam and 12 M dextrorphan, respectively, are achieved. The ability to perform parallel kinetics studies on 384 samples with ~30s time resolution using an isocitrate dehydrogenase 1 (IDH1) enzyme assay is shown. Finally, we demonstrate the repeatability and throughput of our approach by measuring 115,200 samples from 300 microtiter plate reads consecutively over 5.54 hours with RSDs under 8.14% for each freshly introduced plate. Taken together, these results demonstrate the use of IR-MALDESI at sample acquisition rates that surpass other currently reported direct sampling mass spectrometry approaches used for high throughput compound screening.

For Table of Contents Only

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=125 SRC="FIGDIR/small/465730v3_ufig1.gif" ALT="Figure 1">
View larger version (49K):
org.highwire.dtl.DTLVardef@5f916corg.highwire.dtl.DTLVardef@4473b4org.highwire.dtl.DTLVardef@82e8forg.highwire.dtl.DTLVardef@128b57b_HPS_FORMAT_FIGEXP  M_FIG C_FIG 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. Risky choice: probability weighting explains Independence Axiom violations in monkeys Expected Utility Theory (EUT) provides axioms for maximizing utility in risky choice. The Independence Axiom (IA) is its most demanding axiom: preferences between two options should not change when altering both options equally by mixing them with a common gamble. We tested common consequence (CC) and common ratio (CR) violations of the IA over several months in thousands of stochastic choices using a large variety of binary option sets. Three monkeys showed consistently few outright Preference Reversals (8%) but substantial graded Preference Changes (46%) between the initial preferred gamble and the corresponding altered gamble. Linear Discriminant Analysis (LDA) indicated that gamble probabilities predicted most Preference Changes in CC (72%) and CR (88%) tests. The Akaike Information Criterion indicated that probability weighting within Cumulative Prospect Theory (CPT) explained choices better than models using Expected Value (EV) or EUT. Fitting by utility and probability weighting functions of CPT resulted in nonlinear and non-parallel indifference curves (IC) in the Marschak-Machina triangle and suggested IA non-compliance of models using EV or EUT. Indeed, CPT models predicted Preference Changes better than EV and EUT models. Indifference points in out-of-sample tests were closer to CPT-estimated ICs than EV and EUT ICs. Finally, while the few outright Preference Reversals may reflect the long experience of our monkeys, their more graded Preference Changes corresponded to those reported for humans. In benefitting from the wide testing possibilities in monkeys, our stringent axiomatic tests contribute critical information about risky decision-making and serve as basis for investigating neuronal decision mechanisms. Molecular diversity maintained by long-term balancing selection in mating loci defines multiple mating types in fungi Balancing selection, an evolutionary force that retains genetic diversity, has been detected in multiple genes and organisms, such as the sexual mating loci in fungi. However, to quantify the strength of balancing selection and define the mating-related genes require a large number of specimens. In tetrapolar basidiomycete fungi, sexual type is determined by two unlinked loci, MATA and MATB. Genes in both loci defines mating type identity, control successful mating and completion of the life cycle. These loci are usually highly diverse. Previous studies have speculated, based on culture crosses, that species of the non-model genus Trichaptum (Hymenochaetales, Basidiomycota) possess a tetrapolar mating system, with multiple alleles. Here, we sequenced a hundred and eighty specimens of three Trichaptum species. We characterized the chromosomal location of MATA and MATB, the molecular structure of MAT regions and their allelic richness. Our sequencing effort was sufficient to molecularly characterize multiple MAT alleles segregating before the speciation event of Trichaptum species. Our analyses suggested that long-term balancing selection has generated trans-species polymorphisms. Mating sequences were classified in different allelic classes based on an amino acid identity (AAI) threshold supported by phylogenetics. The inferred allelic information mirrored the outcome of in vitro crosses, thus allowing us to support the degree of allelic divergence needed for successful mating. Even with the high amount of divergence, key amino acids in functional domains are conserved. The observed allelic classes could potentially generate 14,560 different mating types. We conclude that the genetic diversity of mating in Trichaptum loci is due to long-term balancing selection, with limited recombination and duplication activity. Our large number of sequenced specimens highlighted the importance of sequencing multiple individuals from different species to detect the mating-related genes, the mechanisms generating diversity and the evolutionary forces maintaining them.

Author summaryFungi have complex mating systems, and basidiomycete fungi can encode thousands of mating types. Individuals with the same mating type cannot mate. This sexual system has evolved to facilitate sexual mating, increasing the chances to recombine into advantageous allelic combination and prune deleterious alleles. We explored the genomes of hundred and eighty specimens, combined with experimental mating studies of selected specimens, from a non-model organism (Trichaptum). We characterized the genomic regions controlling sex. The mating ability of the specimens confirmed the role of the mating alleles observed in the genomic data. The detailed analyses of many specimens allowed us to observe gene duplication and rearrangements within the mating loci, increasing the diversity within these loci. We supported previous suggestions of balancing selection in this region, an evolutionary force that maintains genomic diversity. These results supports that our fungal specimens are prone to outcross, which might facilitate the adaptation to new conditions. Inhibition of BCMA signaling alters protein translation and Multiple Myeloma growth Recent T Cell therapies have been effective in the treatment of hematological cancers. However, immunotoxicity and treatment relapse pose significant clinical challenges. Here, we revealed distinctive requirement for neutralizing TNF receptor ligands APRIL and BAFF in MM and DLBCL. Furthermore, we investigated the use of BCMA decoy receptor (sBCMA-Fc) as a therapeutic inhibitor of ARPIL and BAFF. While wild-type sBCMA-Fc successfully blocked APRIL signaling with picomolar binding affinity, inhibiting tumor growth in MM models, it lacked efficacy in inhibiting DLBCL progression due to its weak binding for BAFF. To expand the therapeutic utility of sBCMA-Fc, using a protein engineering approach, we generated an affinity enhanced mutant sBCMA-Fc fusion molecule (sBCMA-Fc V3) with 4-folds and 500-folds enhancement in binding to APRIL and BAFF respectively. The sBCMA-Fc V3 clone significantly enhanced antitumor activity against both MM and DLBCL. Importantly, sBCMA-Fc V3 was proven to be a viable clinical candidate by showing adequate toxicity profile and on-target mechanism of action in nonhuman primates.

SUMMARYThis study demonstrates the dichotomous function of APRIL and BAFF in MM and DLBCL, that can be safely targeted by an engineered fusion protein designed to trap APRIL and BAFF with ultra-high binding affinity. SERCA function is impaired in skeletal and cardiac muscles from young DBA/2J mdx mice The DBA/2J (D2) mdx mouse has emerged as a more severe model of Duchenne muscular dystrophy when compared to the traditional C57BL/10 (C57) mdx mouse. Here, we questioned whether sarco(endo)plasmic reticulum Ca2+-ATPase (SERCA) function would differ in muscles from young D2 and C57 mdx mice. In gastrocnemius muscles, both D2- and C57 mdx mice exhibited signs of impaired Ca2+ uptake, however, this was more severe in D2 mdx mice. Maximal SERCA activity was lowered only in D2 mdx gastrocnemius muscles and not C57 mdx muscles. Furthermore, in the left ventricle and diaphragm, Ca2+ uptake was impaired in C57 mdx muscles with lowered rates of Ca2+ uptake compared with C57 WT mice, whereas in muscles from D2 mdx mice, rates of Ca2+ uptake were unattainable due to the severe impairments in their ability to transport Ca2+. Overall, our study demonstrates that SERCA function is drastically impaired in young D2 mdx mice. A Hierarchy of Biomolecular Proportional-Integral-Derivative Feedback Controllers for Robust Perfect Adaptation and Dynamic Performance Proportional-Integral-Derivative (PID) feedback controllers have been the most widely used controllers in industry for almost a century due to their good performance, simplicity, and ease of tuning. Motivated by their success in various engineering disciplines, PID controllers recently found their way into synthetic biology, where the design of feedback molecular control systems has been identified as an important goal. In this paper, we consider the mathematical realization of PID controllers via biomolecular interactions. We propose an array of topologies that offer a compromise between simplicity and high performance. We first demonstrate that different Proportional-Integral (PI) controllers exhibit different capabilities for enhancing the dynamics and reducing variance (cell-to-cell variability). Next, we introduce several derivative controllers that are realized based on incoherent feedforward loops acting in a feedback configuration. Alternatively, we show that differentiators can be realized by placing molecular integrators in a negative feedback loop--an arrangement that can then be augmented by PI components to yield PID feedback controllers. We demonstrate that the derivative component can be exploited for enhancing system stability, dramatically increasing the molecular control systems dynamic performance, and for reducing the noise effect on the output of interest. The PID controller features are established through various deterministic and stochastic analyses as well as numerical simulations. Finally, we provide an experimental demonstration using a recently developed hybrid setup, the cyberloop, where the controller is implemented in silico to control a biological genetic circuit in vivo. The large array of novel biomolecular PID controllers introduced here forms a basis for the design and construction of advanced high-performance biomolecular control systems that robustly regulate the dynamics of living systems. Nanocrown electrodes for reliable and robust intracellular recording of cardiomyocytes and cardiotoxicity screening Drug-induced cardiotoxicity arises primarily when a compound alters the electrophysiological properties of cardiomyocytes. Features of intracellular action potentials (iAPs) are powerful biomarkers that predict proarrhythmic risks. However, the conventional patch clamp techniques for measuring iAPs are either laborious and low throughput or not suitable for measuring electrically connected cardiomyocytes. In the last decade, a number of vertical nanoelectrodes have been demonstrated to achieve parallel and minimally-invasive iAP recordings. Nanoelectrodes show great promise, but the large variability in success rate, signal strength, and the low throughput of device fabrication have hindered them from being broadly adopted for proarrhythmia drug assessment. In this work, we developed vertically-aligned and semi-hollow nanocrown electrodes that are mechanically robust and made through a scalable fabrication process. Nanocrown electrodes achieve >99% success rates in obtaining intracellular access through electroporation, allowing reliable and simultaneous iAP recordings from up to 57 human pluripotent stem-cell-derived cardiomyocytes (hPSC-CMs). The accuracy of nanocrown electrode recordings is validated by simultaneous patch clamp recording from the same cell. Nanocrown electrodes enable prolonged iAP recording for continual monitoring of the same cells upon the sequential addition of four to five incremental drug doses. In this way, the dose-response data is self-referencing, which avoids the cell-to-cell variations inherent to hPSC-CMs. We are hopeful that this technology development is a step towards establishing an iAP screening assay for preclinical evaluation of drug-induced arrhythmogenicity. Understanding the Spin Crossover Dynamical Effects of the Dioxygen Binding and Activation on HOD Enzyme For the cofactor-free 1-H-3-hydroxy-4-oxoquinaldine-2,4-dioxygenase (HOD), the dioxygen (O2) dependent steps are rate-limiting along with a spin state crossover to the singlet spin state. Here, the primary triplet O2 molecule activation on the 2-methyl-3-hydroxy-4(1H)-quinolone (MHQ) is investigated, and the catalytic role of the intersystem crossing effects is highlighted by directly comparing results from the Born-Oppenheimer dynamics and non-adiabatic surface hopping dynamics. This work confirms non-adiabatic dynamical effects are essential to modulate the O2 activation on the substrate MHQ. The time scale of the equilibration and conversion from triplet to singlet state should be in the range of a few hundreds of femtoseconds. We hope this work provides us a fresh look at the underlying physics of dioxygen activation reactions involving more than one spin state. Muscle activation strategies of the vastus lateralis according to sex AimDespite men exhibiting greater muscle strength and fatigibility than women, it remains unclear if there are sex-based differences in muscle recruitment strategies e.g. motor unit (MU) recruitment and modulation of firing rate (FR) at normalised forces and during progressive increases in force.

MethodsTwenty-nine healthy male and thirty-one healthy female participants (18-35 years) were studied. Intramuscular electromyography was used to record individual motor unit potentials (MUPs) and near fibre MUPs from the vastus lateralis (VL) during 10% and 25% maximum isometric voluntary contractions (MVC), and spike-triggered averaging was used to obtain motor unit number estimates (MUNE) of the VL. Multilevel mixed-effects linear regression models were used to investigate the effects of sex at each contraction level.

ResultsMen exhibited greater muscle strength (p<0.001) and size (p<0.001) than women, with no difference in force steadiness at 10% or 25% MVC. Women had smaller MUs and higher FR at 10% MVC (both p<0.02), similar to that at 25% MVC in MU size (p=0.062) and FR (p=0.031). However, both sexes showed similar increases in MU size and FR when moving from low-to mid-level contractions. There were no sex differences in any near fibre MUP parameters or in MUNE.

ConclusionIn the vastus lateralis, women produce muscle force via different neuromuscular recruitment strategies to men which is characterised by smaller MUs discharging at higher rates. However, similar strategies are employed to increase force production from low to moderate contractions. These findings of similar proportional increases between sexes support the use of mixed sex cohorts in studies of this nature.

Key pointsO_LIIncreases in muscle force production are mediated by motor unit (MU) recruitment, and MU firing rate (FR).
C_LIO_LIWomen are underrepresented in studies of human neuromuscular research and markedly differ to men in a number of aspects of neuromuscular function, yet little is known of the recruitment strategies of each.
C_LIO_LIHere we demonstrate men and women have similar vastus lateralis MU number estimates, yet women recruit smaller MUs with higher FR than men at normalised contraction levels. However, increases in force are achieved via similar trajectories of MU recruitment and MU FR in men and women.
C_LIO_LIAlthough men and women exhibit divergent neuromuscular recruitment strategies to achieve normalised forces, increases in force are achived similarly and support the inclusion of mixed sex cohorts in studies of this nature.
C_LI The activity and functions of soil microbial communities in the Finnish sub-Arctic vary across vegetation types Increased microbial activity in high-latitude soils due to climate change might lead to higher greenhouse gas (GHG) emissions. However, mechanisms of microbial GHG production and consumption in tundra soils are not thoroughly understood. We analyzed 116 soil metatranscriptomes from 73 sites in the Finnish sub-Arctic to investigate how the diversity and functional potential of bacterial and archaeal communities vary across vegetation types and soil layers. Soils differed in physicochemical conditions, with meadow soils being characterized by higher pH and low soil organic matter (SOM) and carbon/nitrogen ratio whereas dwarf shrub-dominated ecosystems with high SOM and low pH. Actinobacteria, Acidobacteria, Alphaproteobacteria, and Planctomycetes predominated all communities but there were significant differences on genus level between vegetation types, as plant polymer degrading groups were more active in shrub-dominated soils compared to meadows. Given that climate change scenarios predict expansion in dwarf shrubs at high latitudes, our results indicate that the rate of carbon turnover in tundra soils may increase in the future. Additionally, transcripts of methanotrophs were detected in the mineral layer of all soils, potentially moderating methane fluxes from deeper layers. In all, this study provides new insights into possible shifts in tundra microbial diversity and activity with climate change. Simpati: patient classifier identifies signature pathways based on similarity networks for the disease prediction BACKGROUNDPathway-based patient classification is a supervised learning task which implies a model learning pathways as features to predict the classes of patients. The counterpart of enrichment tools for the pathway analysis are fundamental methods for clinicians and biomedical scientists. They allow to find signature cellular functions which help to define and annotate a disease phenotype. They provide results which lead human experts to manually classify patients. It is a paradox that pathwaybased classifiers which natively resolve this objective are not strongly developed. They could simulate the human way of thinking, decipher hidden multivariate relationships between the deregulated pathways and the disease phenotype, and provide more information than a probability value. Instead, there are currently only two classifiers of such kind, they require a nontrivial hyperparameter tuning, are difficult to interpret and lack in providing new insights. There is the need of new classifiers which can provide novel perspectives about pathways, be easy to apply with different biological omics and produce new data enabling a further analysis of the patients.

RESULTSWe propose Simpati, an innovative and interpretable patient classifier based on pathway-specific patient similarity networks. The first classifier to adopt ad-hoc novel algorithms for such graph type. It standardizes the biological high-throughput dataset of patients profiles with a propagation algorithm that considers the interconnected nature of the cells molecules for inferring a new activity score. This allows Simpati to classify with dense, sparse, and non-homogenous omic data. Simpati organizes patients molecules in pathways represented by patient similarity networks for being interpretable, handling missing data and preserving the patient privacy. A network represents patients as nodes and a novel similarity measure determines how much every pair act co-ordinately in a pathway. Simpati detects signature biological processes based on how much the topological properties of the related networks separate the patient classes. In this step, it includes a new cohesive subgroup detection algorithm to handle patients not showing the same pathway activity as the other class members. An unknown patient is then classified by a unique recommender system which considers how much is similar to known patients and distant from being an outlier. Simpati outperforms previously published classifiers on five cancer datasets described with two biological omics, classifies well with sparse data, identifies more relevant pathways associated to the patients disease than the competitors and has the lowest computational requirements.

CONCLUSIONSimpati can serve as generic-purpose pathway-based classifier of patient classes. It provides signature pathways to unveil the altered biological mechanisms of a disease phenotype and to classify patients according to the learnt pathway-specific similarities. The signature condition and patient prediction can be deciphered considering the patient similarity networks which must reveal the members of a patient class more cohesive and similar than the non-members. Simpati divides the pathways in up and downinvolved. Upinvolved when the signaling cascades generated by the altered molecules of the disease patients impact stronger the pathway than the ones of the control class. We provide an R implementation, a graphical user interface and a visualization function for the patient similarity networks. The software is available at: https://github.com/LucaGiudice/Simpati Forecasting the effects of water regulation on the population viability of a threatened amphibian The regulation of river systems alters hydrodynamics and often reduces lateral connectivity between river channels and floodplains. For taxa such as frogs that rely on floodplain wetlands to complete their lifecycle, decreasing inundation frequency can reduce recruitment and increase the probability of local extinction. We virtually reconstructed the inundation patterns of wetlands under natural and regulated flow conditions and built stochastic population models to quantify the probability of local extinction under different inundation scenarios. Specifically, we explored the interplay of habitat size, inundation frequency, and successive dry years on the local extinction probability of the threatened southern bell frog Litoria raniformis in the Murray River floodplains of South Australia. We hypothesised that the changes to wetland inundation resulting from river regulation are a principal driver of L. raniformis declines in this semi-arid system.

Regulation has reduced the inundation frequency of essential habitats below critical thresholds for the persistence of many fresh water-dependent species. Successive dry years raise the probability of local extinction, and these effects are strongest in smaller wetlands. Larger wetlands and those with more frequent average inundation are less susceptible to these effects.

Elucidating these trends informs the prioritisation of treatment sites and the frequency of conservation interventions. Environmental water provision (through pumping or the operation of flow-regulating structures) is a promising tool to reduce the probability of breeding failure and local extinction. Our modelling approach can be used to prioritise the delivery of environmental water to L. raniformis and potentially many other frog species. Loss of the DYRK1A protein kinase results in reduction of ribosomal protein genes expression, ribosome mass and translation defects Ribosomal proteins (RPs) are evolutionary conserved proteins that are essential for protein translation. RP expression must be tightly regulated to ensure the appropriate assembly of ribosomes and to respond to the growth demands of cells. The elements regulating the transcription of RP genes (RPGs) have been characterized in yeast and Drosophila, yet how cells regulate the production of RPs in mammals is less well understood. Here, we show that a subset of RPG promoters is characterized by the presence of the palindromic TCTCGCGAGA motif and marked by the recruitment of the protein kinase DYRK1A. The presence of DYRK1A at these promoters is associated with enhanced binding of the TATA-binding protein, TBP, and it is negatively correlated with the binding of the GABP transcription factor, establishing at least two clusters of RPGs that could be coordinately regulated. However, DYRK1A silencing leads to a global reduction of RPGs, mRNA pointing at DYRK1A activities beyond those dependent on its chromatin association. Significantly, cells in which DYRK1A is depleted have reduced RP levels, fewer ribosomes, reduced global protein synthesis and smaller size. We therefore propose a novel role for DYRK1A in coordinating the expression of genes encoding RPs, thereby controlling cell growth in mammals. 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. Intestine-enriched apolipoprotein b orthologs are required for stem cell differentiation and regeneration in planarians Lipid metabolism plays an instructive role in regulating stem cell state and differentiation. However, the roles of lipid mobilization and utilization in stem cell-driven regeneration are unclear. Planarian flatworms readily restore missing tissue due to injury-induced activation of pluripotent somatic stem cells called neoblasts. Here, we identify two intestine-enriched orthologs of apolipoprotein b, apob-1 and apob-2, which mediate transport of neutral lipid stores from the intestine to target tissues including neoblasts, and are required for tissue homeostasis and regeneration. Inhibition of apob function by RNAi causes head regression and lysis in uninjured animals, and delays body axis re-establishment and regeneration of multiple organs in amputated fragments. Furthermore, apob RNAi causes expansion of the population of differentiating neoblast progeny and dysregulates expression of genes enriched in differentiating and mature cells in eight major cell type lineages. We conclude that intestine-derived lipids serve as a source of metabolites required for neoblast differentiation. Frequency modulated timer regulates mammalian hibernation Mammalian hibernators decrease basal metabolism and body temperature (Tb) to minimize energy expenditure in harsh seasons. During hibernation, Tb drops to low temperature (<10 {degrees}C) and remains constant for days, known as deep torpor in small mammalian hibernators. Spontaneous interbout arousals interrupt torpor bouts, when Tb recovers to euthermic state ~37 {degrees}C. Torpor-interbout arousal event repeats during several months of hibernation. However, little is known about mechanisms governing Tb fluctuation across torpor-interbout arousal cycles during hibernation. Recent improvement in data-logging techniques enables us to monitor Tb for more than hundred days with high precision, opening up new avenues for quantitative analysis to address the principle governing Tb fluctuation. Here, we analyzed Tb fluctuation across torpor-interbout arousal cycle of Syrian hamster, which can hibernate in response to chronic cold and short photoperiod under a laboratory condition, using generalized harmonic analysis and discovered a model with frequency modulation quantitatively reproducing Tb fluctuation. This analysis also identified an unexpectedly longer period of 120-430 days as the period that modulates another period of several days, generating Tb fluctuation for Syrian hamster. We propose that concerted action of two endogenous periods governs torpor-interbout arousal cycles during hibernation. Localization of signaling receptors maximizes cellular informationacquisition in spatially-structured natural environments Cells in natural environments like tissue or soil sense and respond to extracellular ligands with intricately structured and non-monotonic spatial distributions, sculpted by processes such as fluid flow and substrate adhesion. In this work, we show that spatial sensing and navigation can be optimized by adapting the spatial organization of signaling pathways to the spatial structure of the environment. We develop an information-theoretic framework for computing the optimal spatial organization of a sensing system for a given signaling environment. We find that receptor localization maximizes information acquisition in simulated natural contexts, including tissue and soil. Receptor localization extends naturally to produce a dynamic protocol for continuously redistributing signaling receptors, which when implemented using simple feedback, boosts cell navigation efficiency by 30-fold. Broadly, our work shows how cells can maximize the fidelity of information transfer by adapting the spatial organization of signaling molecules to the spatial structure of the environment. An Essential Role of UBXN3B in B Lymphopoiesis Hematopoiesis is finely regulated to enable timely production of the right numbers and types of mature immune cells to maintain tissue homeostasis. Dysregulated hematopoiesis may compromise antiviral immunity and/or exacerbate immunopathogenesis. Herein, we report an essential role of UBXN3B in maintenance of hematopoietic homeostasis and restriction of immunopathogenesis during respiratory viral infection. Ubxn3b deficient (Ubxn3b-/-) mice are highly vulnerable to SARS-CoV-2 and influenza A infection, characterized by more severe lung immunopathology, lower virus-specific IgG, significantly fewer B cells, but more myeloid cells than Ubxn3b+/+ littermates. This aberrant immune compartmentalization is recapitulated in uninfected Ubxn3b-/- mice. Mechanistically, UBXN3B controls precursor B-I (pre-BI) transition to pre-BII and subsequent proliferation in a cell-intrinsic manner, by maintaining BLNK protein stability and pre-BCR signaling. These results reveal an essential role of UBXN3B for the early stage of B cell development.