How we pay attention in naturalistic visual search settings Research on attentional control has largely focused on single senses and the importance of behavioural goals in controlling attention. However, everyday situations are multisensory and contain regularities, both likely influencing attention. We investigated how visual attentional capture is simultaneously impacted by top-down goals, the multisensory nature of stimuli, and the contextual factors of stimulis semantic relationship and temporal predictability. Participants performed a multisensory version of the Folk et al. (1992) spatial cueing paradigm, searching for a target of a predefined colour (e.g. a red bar) within an array preceded by a distractor. We manipulated: 1) stimulis goal-relevance via distractors colour (matching vs. mismatching the target), 2) stimulis multisensory nature (colour distractors appearing alone vs. with tones), 3) the relationship between the distractor sound and colour (arbitrary vs. semantically congruent) and 4) the temporal predictability of distractor onset. Reaction-time spatial cueing served as a behavioural measure of attentional selection. We also recorded 129-channel event-related potentials (ERPs), analysing the distractor-elicited N2pc component both canonically and using a multivariate electrical neuroimaging framework. Behaviourally, arbitrary target-matching distractors captured attention more strongly than semantically congruent ones, with no evidence for context modulating multisensory enhancements of capture. Notably, electrical neuroimaging of surface-level EEG analyses revealed context-based influences on attention to both visual and multisensory distractors, in how strongly they activated the brain and type of activated brain networks. For both processes, the context-driven brain response modulations occurred long before the N2pc time-window, with topographic (network-based) modulations at ~30ms, followed by strength-based modulations at ~100ms post-distractor onset. Our results reveal that both stimulus meaning and predictability modulate attentional selection, and they interact while doing so. Meaning, in addition to temporal predictability, is thus a second source of contextual information facilitating goal-directed behaviour. More broadly, in everyday situations, attention is controlled by an interplay between ones goals, stimulis perceptual salience, meaning and predictability. Our study calls for a revision of attentional control theories to account for the role of contextual and multisensory control.
Brain expressed FKBP5 delineates a therapeutic subtype of severe mental illness Deducing genes capable of classifying biologically-distinct psychiatric subtypes, and their targets for treatment, is a priority approach in psychiatry. FKBP5 is one such gene with strong evidence of utility to delineate a trans-diagnostic psychiatric subtype. Yet how brain-expressed FKBP5 is affected in psychiatric disorders in humans is not fully understood and critical for propelling FKBP5-targeting treatment development. We performed a large-scale postmortem study (n=895) of FKBP5 using dorsolateral prefrontal cortex samples derived from individuals with severe psychiatric disorders with a comprehensive battery of bulk/single-cell omics and histological analyses. We observed consistently heightened FKBP5 mRNA and protein in psychopathology, moderated by genotype and age, and accompanied by DNA methylation changes in key enhancers. These effects were most prominent in superficial-layer pyramidal cells. Heightened FKBP5 was also differentially associated with downstream pathways according to age, being specifically associated with synaptic transmission in early adulthood, and neuroinflammation and neurodegeneration in later life.
Nucleus size and its effect on the chromatin structure in living cells. DNA-architectural proteins play a major role in organization of chromosomal DNA in living cells by packaging it into chromatin, whose spatial conformation is determined by an intricate interplay between the DNA-binding properties of architectural proteins and physical constraints applied to the DNA by a tight nuclear space. Yet, the exact effects of the cell nucleus size on DNA-protein interactions and chromatin structure currently remain obscure. Furthermore, there is even no clear understanding of molecular mechanisms responsible for the nucleus size regulation in living cells. To find answers to these questions, we developed a general theoretical framework based on a combination of polymer field theory and transfer-matrix calculations, which showed that the nucleus size is mainly determined by the difference between the surface tensions of the nuclear envelope and the endoplasmic reticulum membrane as well as the osmotic pressure exerted by cytosolic macromolecules on the nucleus. In addition, the model demonstrated that the cell nucleus functions as a piezoelectric element, changing its electrostatic potential in a size-dependent manner. This effect has been found to have a profound impact on stability of nucleosomes, revealing a previously unknown link between the nucleus size and chromatin structure. Overall, our study provides new insights into the molecular mechanisms responsible for regulation of the cell nucleus size, as well as the potential role of nuclear organization in shaping the cell response to environmental cues.

SIGNIFICANCE STATEMENTThe cell nucleus plays a central role in the life of eukaryotic cells, providing the highest level of control of intracellular processes. Depending on the stage of the cell cycle and / or surrounding environment, the size of the cell nucleus may undergo changes that are believed to cause chromatin reorganization, affecting gene transcription. At present, however, there is no clear understanding of the molecular mechanisms that may be responsible for such regulation, whose exact effect on chromatin structure remains unclear. In this study, by developing an advanced computational approach, we explore these issues from a physical perspective, revealing previously unknown mechanisms contributing to organization of the cell nucleus and chromatin. Penalized Logistic Regression Analysis for Genetic Association Studies of Binary Phenotypes IntroductionIncreasingly, logistic regression methods for genetic association studies of binary phenotypes must be able to accommodate data sparsity, which arises from unbalanced case-control ratios and/or rare genetic variants. Sparseness leads to maximum likelihood estimators (MLEs) of log-OR parameters that are biased away from their null value of zero and tests with inflated type 1 errors. Different penalized-likelihood methods have been developed to mitigate sparse-data bias. We study penalized logistic regression using a class of log-F priors indexed by a shrinkage parameter m to shrink the biased MLE towards zero.

MethodsWe propose a two-step approach to the analysis of a genetic association study: first, a set of variants that show evidence of association with the trait is used to estimate m; and second, the estimated m is used for log-F -penalized logistic regression analyses of all variants using data augmentation with standard software. Our estimate of m is the maximizer of a marginal likelihood obtained by integrating the latent log-ORs out of the joint distribution of the parameters and observed data. We consider two approximate approaches to maximizing the marginal likelihood: (i) a Monte Carlo EM algorithm (MCEM) and (ii) a Laplace approximation (LA) to each integral, followed by derivative-free optimization of the approximation.

ResultsWe evaluate the statistical properties of our proposed two-step method and compared its performance to other shrinkage methods by a simulation study. Our simulation studies suggest that the proposed log-F -penalized approach has lower bias and mean squared error than other methods considered. We also illustrate the approach on data from a study of genetic associations with "super senior" cases and middle aged controls.

Discussion/ConclusionWe have proposed a method for single rare variant analysis with binary phenotypes by logistic regression penalized by log-F priors. Our method has the advantage of being easily extended to correct for confounding due to population structure and genetic relatedness through a data augmentation approach. Dose-dependent thresholds of dexamethasone destabilize CAR T-cell treatment efficacy Chimeric antigen receptor (CAR) T-cell therapy is potentially an effective targeted immunotherapy for glioblastoma, yet there is presently little known about the efficacy of CAR T-cell treatment when combined with the widely used anti-inflammatory and immunosuppressant glucocorticoid, dexamethasone. Here we present a mathematical model-based analysis of three patient-derived glioblastoma cell lines treated in vitro with CAR T-cells and dexamethasone. Advanced in vitro experimental cell killing assay technologies allow for highly resolved temporal dynamics of tumor cells treated with CAR T-cells and dexamethasone, making this a valuable model system for studying the rich dynamics of nonlinear biological processes with translational applications. We model the system as a nonautonomous, two-species predator-prey interaction of tumor cells and CAR T-cells, with explicit time-dependence in the clearance rate of dexamethasone. Using time as a bifurcation parameter, we show that (1) dexamethasone destabilizes coexistence equilibria between CAR T-cells and tumor cells in a dose-dependent manner and (2) as dexamethasone is cleared from the system, a stable coexistence equilibrium returns in the form of a Hopf bifurcation. With the model fit to experimental data, we demonstrate that high concentrations of dexamethasone antagonizes CAR T-cell efficacy by exhausting, or reducing the activity of CAR T-cells, and by promoting tumor cell growth. Finally, we identify a critical threshold in the ratio of CAR T-cell death to CAR T-cell proliferation rates that predicts eventual treatment success or failure that may be used to guide the dose and timing of CAR T-cell therapy in the presence of dexamethasone in patients.

Author summaryBioengineering and gene-editing technologies have paved the way for advance immunotherapies that can target patient-specific tumor cells. One of these therapies, chimeric antigen receptor (CAR) T-cell therapy has recently shown promise in treating glioblastoma, an aggressive brain cancer often with poor patient prognosis. Dexamethasone is a commonly prescribed anti-inflammatory medication due to the health complications of tumor associated swelling in the brain. However, the immunosuppressant effects of dexamethasone on the immunotherapeutic CAR T-cells are not well understood. To address this issue, we use mathematical modeling to study in vitro dynamics of dexamethasone and CAR T-cells in three patient-derived glioblastoma cell lines. We find that in each cell line studied there is a threshold of tolerable dexamethasone concentration. Below this threshold, CAR T-cells are successful at eliminating the cancer cells, while above this threshold, dexamethasone critically inhibits CAR T-cell efficacy. Our modeling suggests that in the presence of high dexamethasone reduced CAR T-cell efficacy, or increased exhaustion, can occur and result in CAR T-cell treatment failure. Citation needed? Wikipedia and the COVID-19 pandemic BackgroundWith the COVID-19 pandemics outbreak, millions flocked to Wikipedia for updated information. Amid growing concerns regarding an "infodemic", ensuring the quality of information is a crucial vector of public health. Investigating if and how Wikipedia remained up to date and in line with science is key to formulating strategies to counter misinformation. Using citation analyses, we asked: which sources informed Wikipedias COVID-19-related articles before and during the pandemics first wave (January-May 2020).

ResultsWe found that coronavirus-related articles referenced trusted media sources and high-quality academic research. Moreover, despite a surge in COVID-19 preprints, Wikipedia had a clear preference for open-access studies published in respected journals and made little use of preprints. Building a timeline of English COVID-19 articles from 2001-2020 revealed a nuanced trade-off between quality and timeliness. It further showed how preexisting articles on key topics related to the virus created a framework for integrating new knowledge. Supported by a rigid sourcing policy, this "scientific infrastructure" facilitated contextualization and regulated the influx of new information. Lastly, we constructed a network of DOI-Wikipedia articles, which showed the shifting landscape of pandemic-related knowledge on Wikipedia and how academic citations create a web of shared knowledge supporting topics like COVID-19 vaccine development.

ConclusionsUnderstanding how scientific research interacts with the digital knowledge-sphere during the pandemic provides insight into how Wikipedia can facilitate access to science. It also reveals how, aided by what we term its "citizen encyclopedists", it successfully fended off COVID-19 disinformation and how this unique model may be deployed in other contexts. Distinct mutations and lineages of SARS-CoV-2 virus in the early phase of COVID-19 pandemic and subsequent global expansion A novel coronavirus, SARS-CoV-2, has caused over 190 million cases and over 4 million deaths worldwide since it occurred in December 2019 in Wuhan, China. Here we conceptualized the temporospatial evolutionary and expansion dynamics of SARS-CoV-2 by taking a series of cross-sectional view of viral genomes from early outbreak in January 2020 in Wuhan to early phase of global ignition in early April, and finally to the subsequent global expansion by late December 2020. Based on the phylogenetic analysis of the early patients in Wuhan, Wuhan/WH04/2020 is supposed to be a more appropriate reference genome of SARS-CoV-2, instead of the first sequenced genome Wuhan-Hu-1. By scrutinizing the cases from the very early outbreak, we found a viral genotype from the Seafood Market in Wuhan featured with two concurrent mutations (i.e. M type) had become the overwhelmingly dominant genotype (95.3%) of the pandemic one year later. By analyzing 4,013 SARS-CoV-2 genomes from different continents by early April, we were able to interrogate the viral genomic composition dynamics of initial phase of global ignition over a timespan of 14-week. 11 major viral genotypes with unique geographic distributions were also identified. WE1 type, a descendant of M and predominantly witnessed in western Europe, consisted a half of all the cases (50.2%) at the time. The mutations of major genotypes at the same hierarchical level were mutually exclusive, which implying that various genotypes bearing the specific mutations were propagated during human-to-human transmission, not by accumulating hot-spot mutations during the replication of individual viral genomes. As the pandemic was unfolding, we also used the same approach to analyze 261,323 SARS-CoV-2 genomes from the world since the outbreak in Wuhan (i.e. including all the publicly available viral genomes) in order to recapitulate our findings over one-year timespan. By 25 December 2020, 95.3% of global cases were M type and 93.0% of M-type cases were WE1. In fact, at present all the four variants of concern (VOC) are the descendants of WE1 type. This study demonstrates the viral genotypes can be utilized as molecular barcodes in combination with epidemiologic data to monitor the spreading routes of the pandemic and evaluate the effectiveness of control measures. Moreover, the dynamics of viral mutational spectrum in the study may help the early identification of new strains in patients to reduce further spread of infection, guide the development of molecular diagnosis and vaccines against COVID-19, and help assess their accuracy and efficacy in real world at real time. Pseudomonas aeruginosa C-terminal processing protease CtpA assembles into a hexameric structure that requires activation by a spiral-shaped lipoprotein binding partner Pseudomonas aeruginosa CtpA is a carboxyl terminal-processing protease that partners with the outer membrane lipoprotein LbcA to degrade at least five cell wall-associated proteins, four of which are cell wall hydrolases. This activity plays an important role in supporting P. aeruginosa virulence in a mouse model of acute pneumonia. However, almost nothing is known about the molecular mechanisms underlying CtpA and LbcA function. Here, we used structural analysis to show that CtpA alone assembles into an inactive hexamer comprising a trimer of dimers, which limits its substrate access and prevents nonspecific degradation. The adaptor protein LbcA is a right-handed open spiral with 11 tetratricopeptide repeats, which might wrap around a substrate to deliver it to CtpA for degradation. By structure-guided mutagenesis and functional assays, we also showed that the interfaces of the CtpA trimer-of-dimers, and an N-terminal helix of LbcA, are important for LbcA-mediated substrate degradation by CtpA both in vitro and in vivo. This work improves our understanding of the molecular mechanism of a CTP within the C-terminal processing peptidase-3 group.

IMPORTANCECarboxyl-terminal processing proteases (CTPs) are found in all three domains of life. In bacteria, some CTPs have been associated with virulence, raising the possibility that they could be theraputic targets. However, relatively little is known about their molecular mechanisms of action. In Pseudomonas aeruginosa, CtpA supports virulence by working in complex with the outer membrane lipoprotein LbcA to degrade cell wall hydrolases. Here, we report structure-function analyses of CtpA and LbcA, which reveals that CtpA assembles into an inactive hexamer comprising a trimer of dimers. LbcA is monomeric, with an N-terminal region important for binding to and activating CtpA, followed by a spiral structure composed of 11 tetratricopetide repeats, which could wrap around a substrate for delivery to CtpA. This work provides the first structure of a CTP-3 group member, revealing a unique mutimeric arrangement and insight into how this important proteolytic system functions. Crowdsourced identification of multi-target kinase inhibitors for RET- and TAU- based disease: the Multi-Targeting Drug DREAM Challenge A continuing challenge in modern medicine is the identification of safer and more efficacious drugs. Precision therapeutics, which have one molecular target, have been long promised to be safer and more effective than traditional therapies. This approach has proven to be challenging for multiple reasons including lack of efficacy, rapidly acquired drug resistance, and narrow patient eligibility criteria. An alternative approach is the development of drugs that address the overall disease network by targeting multiple biological targets ( polypharmacology). Rational development of these molecules will require improved methods for predicting single chemical structures that target multiple drug targets. To address this need, we developed the Multi-Targeting Drug DREAM Challenge, in which we challenged participants to predict single chemical entities that target pro-targets but avoid anti-targets for two unrelated diseases: RET-based tumors and a common form of inherited Tauopathy. Here, we report the results of this DREAM Challenge and the development of two neural network-based machine learning approaches that were applied to the challenge of rational polypharmacology. Together, these platforms provide a potentially useful first step towards developing lead therapeutic compounds that address disease complexity through rational polypharmacology.

Author SummaryMany modern drugs are developed with the goal of modulating a single cellular pathway or target. However, many drugs are, in fact,  dirty; they target multiple cellular pathways or targets. This phenomenon is known as multi-targeting or polypharmacology. While some strive to develop  cleaner therapeutics that eliminate secondary targets, recent work has shown that multi-targeting therapeutics have key advantages for a variety of diseases. However, while multi-targeting drugs that affect a precisely-defined profile of targets may be more effective, it is difficult to computationally predict which molecules have desirable target profiles. Here, we report the results of a competitive crowdsourcing project (the Multi-Targeting Drug DREAM Challenge), where we challenged participants to predict chemicals that have desired target profiles for cancer and neurodegenerative disease. Unevolved proteins from modern and prebiotic amino acids manifest distinct structural profiles Natural proteins represent numerous but tiny structure/function islands in a vast ocean of possible protein sequences, most of which has not been explored by either biological evolution or research. Recent studies have suggested this uncharted sequence space possesses surprisingly high structural propensity, but development of an understanding of this phenomenon has been awaiting a systematic high-throughput approach.

Here, we designed, prepared, and characterized two combinatorial protein libraries consisting of randomized proteins, each 105 residues in length. The first library constructed proteins from the entire canonical alphabet of 20 amino acids. The second library used a subset of only 10 residues (A,S,D,G,L,I,P,T,E,V) that represent a consensus view of plausibly available amino acids through prebiotic chemistry. Our study shows that compact conformations resistant to proteolysis are (i) abundant (up to 40%) in random sequence space, (ii) independent of general Hsp70 chaperone system activity, and (iii) not granted solely by "late" and complex amino acid additions. The Hsp70 chaperone system effectively increases solubility and refoldability of the canonical alphabet but has only a minor impact on the "early" library. The early alphabet proteins are inherently more soluble and refoldable, possibly assisted by the cell-like environment in which these assays were performed.

Our work indicates that both early and modern amino acids are predisposed to supporting protein structure (either in forms of oligomers or globular/molten globule structures) and that protein structure may not be a unique outcome of evolution. Environment dependent benefits of inter-individual variation in honey bee recruitment Inter-individual differences in behaviour within the members of a social group can affect the groups productivity. In eusocial insects, individual differences amongst workers in a colony play a central role division of labour and task allocation. Extensive empirical and theoretical work has highlighted variation in response thresholds as a proximate mechanism underlying individual behavioural differences and hence division of labour. However, other response parameters, like response probability and intensity, can affect these differences. In this study, we first extended a previously published agent-based model on honey bee foraging to understand the relative importance of response (dance) probability and response (dance) intensity in the task of recruitment. Comparing variation obtained from the simulations with previously published empirical data, we found that response intensity plays a more important role than probability in producing consistent inter-individual differences in recruitment behaviour. We then explored the benefits provided by this individual variation in recruitment behaviour to the colonys collective foraging effort under different environmental conditions. Our results revealed that individual variation leads to a greater energetic yield per forager, but only when food is abundant. Our study highlights the need to consider all response parameters while studying division of labour and adds to the growing body of evidence linking individual variation in behavioural responses to the success of social groups. Topaz1, an essential gene for murine spermatogenesis, down-regulates the expression of numerous testis-specific long non-coding RNAs. Spermatogenesis involves coordinated processes, including meiosis, to produce functional gametes. We previously reported Topaz1 as a germ cell-specific gene highly conserved in vertebrates. Topaz1 knockout males are sterile with testes that lack haploid germ cells because of meiotic arrest after prophase I. To better characterize Topaz1-/- testes, we used RNA-sequencing analyses at two different developmental stages (P16 and P18). The absence of TOPAZ1 disturbed the expression of genes involved in microtubule and/or cilium mobility, which was consistent with testicular histology showing the disruption of microtubules and centrosomes. Moreover, a quarter of P18 dysregulated genes are long non-coding RNAs (lncRNAs), and three of them are testis-specific and located in spermatocytes, their expression starting between P11 and P15. The suppression of one of them, 4939463O16Rik, did not alter fertility although sperm parameters were disturbed and sperm concentration fell. The transcriptome of P18-4939463O16Rik-/- testes was altered and the molecular pathways affected included microtubule-based processes, the regulation of cilium movement and spermatogenesis. The absence of TOPAZ1 protein or 4930463O16Rik produced the same enrichment clusters in mutant testes despite a contrasted phenotype on male fertility. In conclusion, TOPAZ1 appeared to stabilize the expression of numerous lncRNAs. Its suppression is not essential for fertility but required during the terminal differentiation of male gametes.

Author SummaryThe Topaz1 gene was initially characterized during the initiation of meiosis in the sheep fetal ovary. In order to determine its function, a KO of the murine gene was performed. In this species, only males were sterile and spermatogenesis was blocked before the first meiotic division. Here, we show that cytoskeletal elements are markedly disturbed in mutant testes, indicating that these elements play an important function in spermatogenesis. While the mitotic spindle of spermatogonia was normal, the meiotic spindle of spermatocytes was hemi-spindle-shaped and the homologous chromosome pairs could position themselves on the equatorial plate. In addition, lncRNAs account for 25% of genes whose expression in testes varies significantly in the absence of Topaz1. This suggests a key role for these factors in spermatogenesis. Largely testis-specific, they may be involved in spermatogenesis and play a more or less critical role in mouse fertility, which probably also depends on their redundancies. The renal inflammatory network of nephronophthisis STRUCTURED ABSTRACTO_ST_ABSBACKGROUNDC_ST_ABSThe majority of genetic kidney disease leading to kidney failure is caused by mutations in ciliary genes. How cilia malfunction leads to progressive kidney damage is poorly understood, but recent evidence links ciliopathy genes to CCL2 dependent macrophage recruitment in autosomal dominant polycystic kidney disease (ADPKD), the most studied renal ciliopathy. Whether or not renal inflammation is involved in other renal ciliopathies is unclear.

METHODSWe combined mice models with kidney biopsies and renal epithelial cells sampled from human urine to characterize the renal inflammatory network of nephronophthisis (NPH), the most frequent renal ciliopathy in children.

RESULTSIn human, mutations in cilia genes involved in NPH enhance urine excretion of the chemokine CCL2, causing abnormal macrophage recruitment in kidney tissues from NPH patients. Differing from ADPKD, inactivating Ccl2 specifically in mouse tubular cells does not rescue the NPH phenotype, suggesting that other inflammatory mediators are involved. Using transcriptional data from 2 NPH models, we identify a set of pro-inflammatory cytokines upregulated in this disease, independently of CCL2. The majority of detectable transcripts from this set are specifically upregulated in kidney cells from NPH patients. In line with the function of these cytokines, NPH kidneys show disproportionate neutrophils and T cells infiltrates compared to healthy subject or hypertensive and diabetic chronic kidney disease patients.

CONCLUSIONSThis study reveals that inflammation is a central aspect in human NPH and delineates a specific set of inflammatory mediators that regulates immune cell recruitment in human NPH.

SIGNIFICANCE STATEMENTMutations in genes encoding primary cilia proteins are the leading cause of genetic kidney failure. In autosomal dominant polycystic kidney disease (ADPKD), deregulated cilia signaling leads to kidney infiltration by macrophages through the chemokine CCL2. Little is known about renal inflammation in nephronophthisis (NPH), the most frequent pediatric renal ciliopathy. Using NPH mice models, tissues and cells from NPH patients, we unveil renal inflammation as preeminent feature of NPH. Remarkably, the renal inflammatory evoked by ciliary gene mutations in NPH does not overlap with ADPKD: it is CCL2 independent, involves a prominent recruitment of neutrophils and T cells and a specific cytokine signature. This unforeseen findings strengthen the link between primary cilia and renal inflammation. Induced antigen-binding polyreactivity in human serum IgA Previous studies have shown that polyreactive antibodies play an important role in the frontline defense against the dissemination of pathogens in the pre-immune host. Interestingly, antigen-binding polyreactivity can not only be inherent, but also acquired post-translationally. The ability of individual monoclonal IgG and IgE antibodies to acquire polyreactivity following contact with various agents that destabilize protein structure (urea, low pH) or have pro-oxidative potential (heme, ferrous ions) has been studied in detail. However, to the best of our knowledge this property of human IgA has previously been described only cursorily. In the present study pooled human serum IgA and two human monoclonal IgA antibodies were exposed to buffers with acidic pH, to free heme or to ferrous ions, and the antigen-binding behavior of the native and modified IgA to viral and bacterial antigens were compared using immunoblot and ELISA. We observed a dose-dependent increase in reactivity to several bacterial extracts and to pure viral antigens. This newly described property of IgA may have therapeutic potential as has already been shown for pooled IgG with induced polyreactivity. Differences in gait stability and acceleration characteristics between healthy young and older females Our aim was to evaluate differences in gait acceleration intensity, variability and stability of feet and trunk between older females and young females using inertial sensors. Twenty older females (OF; mean age 68.4, SD 4.1 years) and eighteen young females (YF; mean age 22.3, SD 1.7 years) were asked to walk straight for 100 meters at their preferred speed, while wearing inertial sensors on heels and lower back. We calculated spatiotemporal measures, foot and trunk acceleration characteristics and their variability, as well as trunk stability using the local divergence exponent (LDE). Two-way analysis of variance (including the factors foot and age), Students t-test, and Mann-Whitney U test were used to compare statistical differences of measures between groups. Cohens d effects were calculated for each variable. Foot maximum vertical acceleration and amplitude, trunk-foot vertical acceleration attenuation, as well as their variability were significantly smaller in OF than in YF. In contrast, trunk mediolateral acceleration amplitude, maximum vertical acceleration, and amplitude, as well as their variability were significantly larger in OF than in YF. Moreover, OF showed lower stability (i.e., higher LDE values) in mediolateral acceleration, mediolateral and vertical angular velocity of the trunk. Even though we measured healthy older females, these participants showed lower vertical foot accelerations with higher vertical trunk acceleration, lower trunk-foot vertical acceleration attenuation, less gait stability, and more variability of the trunk, and hence, were more likely to fall. These findings suggest that instrumented gait measurements may help for early detection of changes or impairments in gait performance, even before this can be observed by clinical eye or gait speed. Multiscale X-ray study of Bacillus subtilis biofilms reveals interlinked structural hierarchy and elemental heterogeneity Biofilms are multicellular microbial communities that encase themselves in an extracellular matrix (ECM) of secreted biopolymers and attach to surfaces and interfaces. Bacterial biofilms are detrimental in hospital and industrial settings, but they can be beneficial in agricultural contexts. An essential property of biofilms that grants them with increased survival relative to planktonic cells is phenotypic heterogeneity; the division of the biofilm population into functionally distinct subgroups of cells. Phenotypic heterogeneity in biofilms can be traced to the cellular level, however, the molecular structures and elemental distribution across whole biofilms as well as possible linkages between them remain unexplored. Mapping X-ray diffraction (XRD) across intact biofilms in time and space, we revealed the dominant structural features in Bacillus subtilis biofilms, stemming from matrix components, spores and water. By simultaneously following the X-ray fluorescence (XRF) signal of biofilms and isolated matrix components, we discovered that the ECM preferentially binds calcium ions over other metal ions, specifically, zinc, manganese and iron. These ions, remaining free to flow below macroscopic wrinkles that act as water channels, eventually accumulate and lead to sporulation. The possible link between ECM properties, regulation of metal ion distribution and sporulation across whole intact biofilms unravels the importance of molecular-level heterogeneity in shaping biofilm physiology and development.

Significance StatementBiofilms are multicellular soft microbial communities that are able to colonize synthetic surfaces as well as living organisms. To survive sudden environmental changes and efficiently share their common resources, cells in a biofilm divide into subgroups with distinct functions, leading to phenotypic heterogeneity. Here, by studying intact biofilms by synchrotron X-ray diffraction and fluorescence, we revealed correlations between biofilm macroscopic architectural heterogeneity and the spatio-temporal distribution of extracellular matrix, spores, water and metal ions. Our findings demonstrate that biofilm heterogeneity is not only affected by local genetic expression and cellular differentiation, but also by passive effects resulting from the physicochemical properties of the molecules secreted by the cells, leading to differential distribution of nutrients that propagates through macroscopic length scales. Building a Knowledge Graph to Enable Precision Medicine Developing personalized diagnostic strategies and targeted treatments requires a deep understanding of disease biology and the ability to dissect the relationship between molecular and genetic factors and their phenotypic consequences. However, such knowledge is fragmented across publications, non-standardized research repositories, and evolving ontologies describing various scales of biological organization between genotypes and clinical phenotypes. Here, we present PrimeKG, a precision medicine-oriented knowledge graph that provides a holistic view of diseases. PrimeKG integrates 20 high-quality resources to describe 17,080 diseases with 4,050,249 relationships representing ten major biological scales, including disease-associated protein perturbations, biological processes and pathways, anatomical and phenotypic scales, and the entire range of approved and experimental drugs with their therapeutic action, considerably expanding previous efforts in disease-rooted knowledge graphs. In addition, PrimeKG supports artificial intelligence analyses of how drugs might target disease-associated molecular perturbations by containing an abundance of 'indications', 'contradictions', and 'off-label use' drug-disease edges lacking in other knowledge graphs. We accompany PrimeKG's graph structure with text descriptions of clinical guidelines to enable multimodal analyses. Metacommunities from bacteria to birds: stronger environmental effects in mediterranean than in tropical ponds. The metacommunity concept provides a theoretical framework that aims at explaining organism distributions by a combination of environmental filtering, dispersal and drift. With the development of statistical tools to quantify and partially isolate the role of each of these processes, empirical metacommunity studies have multiplied worldwide. However, few works attempt a multi-taxon approach and even fewer compare two distant biogeographical regions using the same methodology. Under this framework, we tested the expectation that temperate (mediterranean-climate) pond metacommunities would be more influenced by environmental and spatial processes than tropical ones, because of stronger environmental gradients and greater isolation of waterbodies.

We surveyed 30 tropical and 32 mediterranean temporary ponds from Costa Rica and Spain, respectively, and obtained data on 49 environmental variables (including limnological, hydrogeomorphological, biotic, climatic, and landscape variables). We characterized the biological communities of Bacteria and Archaea (from both the water column and the sediments), phytoplankton, zooplankton, benthic invertebrates, amphibians and birds, and estimated the relative role of space and environment on metacommunity organization for each group and region, by means of variation partitioning using Generalized Additive Models (GAMs).

Environmental selection was important in both tropical and mediterranean ponds, but markedly stronger in the latter, probably due to their larger limnological heterogeneity. Spatialized environment and pure spatial effects were greater in the tropics, related to higher climatic heterogeneity and dispersal processes (e.g. restriction, surplus) acting at different scales. The variability between taxonomic groups in spatial and environmental contributions was very wide. Effects on passive and active dispersers were similar within regions but different across regions, with higher environmental effects in mediterranean active dispersers. The residual (unexplained) variation was larger in tropical pond metacommunities, suggesting a higher role for stochastic processes and/or effects of biotic interactions in the tropics. Overall, these results provide support, for a wide variety of organisms related to aquatic habitats, for the classical view of stronger abiotic niche constraints in temperate areas compared to the tropics. Pan-cancer HLA Gene-mediated Tumor Immunogenicity and Immune Evasion Human Leukocyte Antigen (HLA) expression contributes to the activation of anti-tumor immunity through interactions with T cell receptors. However, pan-cancer HLA-mediated immunogenicity and immunoediting mechanisms has not been systematically studied. In a retrospective analysis of 33 tumor types from the Cancer Genome Atlas, we uncovered HLA class I and class II differential expression, which outperformed traditional clinical metrics in predicting patient survival. We also characterized the distribution of HLA supertypes across cancers and showed that patients with high HLA allelic diversity and gene expression had better prognosis. Immune microenvironments with varied survival outcomes could be predicted using a neural network model trained on HLA expression data. Furthermore, we identified a subset of tumors which upregulated HLA class I but not class II genes and exploited HLA-mediated escape strategies. Our results suggest the potential of using HLA expression data to predict patient prognosis. Taken together, we emphasize the crucial role of HLA upregulation in shaping prolonged anti-tumor immunity.

SynopsisIn a retrospective analysis of 11080 patients of 33 TCGA cancer types, we showed HLA class I and class II differential expression shape various immune microenvironments and tumor immunoediting mechanisms, predicting tumor immunogenicity and patient survival. The evolution of stomatal traits along the trajectory towards C4 photosynthesis C4 photosynthesis optimizes plants carbon and water relations, allowing high photosynthetic rate with low stomatal conductance. Stomata have long been believed as a part of C4 syndrome. However, it remains unclear how stomata traits evolved along the path from C3 to C4. Stomatal patterning was examined in Flaveria genus, a model for studying C4 evolution. Comparative, transgenic and semi-in-vitro experiments were used to study molecular basis that underlies stomatal traits along C4 evolution. Novel results: the evolution from C3 to C4 species through intermediate species is accompanied by a stepwise rather than an abrupt change in the stomatal traits. The initial change occurs near Type II and dramatic change occurs at the C4-like species. On the road to C4, stomata become less in number but bigger in size and changes in stomatal density dominates changes in maximum stomatal conductance (gsmax). The reduction of FSTOMAGEN expression underlies altered gsmax between Flaveria species with different photosynthetic pathways and likely occurs in other C4 lineages. Our study provides insight into the pattern, mechanism and role of stomatal evolution along the road towards C4. This work highlights the stomatal traits in the current C4 evolutionary model and the co-evolution of photosynthetic pathway and stomata. Morphology of immatures of the thelytokous ant, Monomorium triviale Wheeler (Formicidae: Myrmicinae: Solenopsidini) with descriptions of the extraordinary last-instar queen larvae The ant genus Monomorium is one of the most species-rich but taxonomically problematic groups in the hyperdiverse subfamily Myrmicinae. An East Asian species, M. triviale Wheeler, produces both reproductive queens and sterile workers via obligate thelytokous parthenogenesis. Here, we describe the immature forms of M. triviale based on light and scanning electron microscopy observations, with a note on the striking caste dimorphism in the last larval instar. The last-instar queen larvae were easily recognized by their large size, "aphaenogastroid" body shape, and rows of doorknob-like tubercles on the lateral and dorsal body surface. This type of queen-specific structure has not been found in ants in general, let alone congeneric species found in Japan. In stark contrast to the queen larvae, worker larvae showed a "pheidoloid" body shape and a body surface similar to other ants. The worker larvae were estimated to have three instars, consistent with previously described congeners. The pupae of both castes had no cocoon, a characteristic commonly described in other Myrmicinae species. In total, the developmental period from egg to adult worker averaged 59 days under 25{degrees}C. We discuss possible functions of the tubercles of queen larvae based on previous studies. An Evolutionary Process Without Variation and Selection Natural selection successfully explains how organisms accumulate adaptive change despite that traits acquired over a lifetime are eliminated at the end of each generation. However, in some domains that exhibit cumulative, adaptive change--e.g., cultural evolution, and earliest life--acquired traits are retained; these domains do not face the problem that Darwins theory was designed to solve. Lack of transmission of acquired traits occurs when germ cells are protected from environmental change, due to a self-assembly code used in two distinct ways: (i) actively interpreted during development to generate a soma, and (ii) passively copied without interpretation during reproduction to generate germ cells. Early life and cultural evolution appear not to involve a self-assembly code used in these two ways. We suggest that cumulative, adaptive change in these domains is due to a lower-fidelity evolutionary process, and model it using Reflexively Autocatalytic and Foodset-generated networks. We refer to this more primitive evolutionary process as Self-Other Reorganisation (SOR) because it involves internal self-organising and self-maintaining processes within entities, as well as interaction between entities. SOR encompasses learning but in general operates across groups. We discuss the relationship between SOR and Lamarckism, and illustrate a special case of SOR without variation. RNA polymerase II is recruited to DNA double-strand breaks for dilncRNA transcription in Drosophila DNA double-strand breaks are among the most toxic lesions that can occur in a genome and their faithful repair is thus of great importance. Recent findings have uncovered local transcription that initiates at the break and forms a non-coding transcript, called damage-induced long non-coding RNA (dilncRNA), which helps to coordinate the DNA transactions necessary for repair. We provide nascent RNA sequencing-based evidence that RNA polymerase II transcribes the dilncRNA in Drosophila and that this is more efficient for DNA breaks in an intron-containing gene, consistent with the higher damage-induced siRNA levels downstream of an intron. The spliceosome thus stimulates recruitment of RNA polymerase II to the break, rather than merely promoting the annealing of sense and antisense RNA to form the siRNA precursor. In contrast, RNA polymerase III nascent RNA libraries did not contain reads corresponding to the cleaved loci and selective inhibition of RNA polymerase III did not reduce the yield of damage-induced siRNAs. Finally, the damage-induced siRNA density was unchanged downstream of a T8 sequence, which terminates RNA polymerase III transcription. We thus found no evidence for a participation of RNA polymerase III in dilncRNA transcription in cultured Drosophila cells. 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. Co-infection of Candidatus Piscichlamydia trichopodus (order Chlamydiales) and Henneguya sp (Myxosporea, Myxobolidae) in snakeskin gourami Trichopodus pectoralis (Regan 1910) The present study describes a simultaneous infection of a novel Chlamydia-like organism (CLO) with a Myxozoa parasite, Henneguya sp. in snakeskin gourami Trichopodus pectoralis in Thailand. A new CLO is proposed "Candidatus Piscichlamydia trichopodus" (CPT) based on 16S rRNA phylogenetic analysis. Systemic intracellular CPT infection was confirmed by histological examination, in situ hybridization, PCR assay, and sequencing of 16S rRNA. This novel pathogen belongs to the order Chlamydiales but differs in certain aspects from other species. The histopathological changes associated with CPT infection were different from the typical pathological lesions of epitheliocystis caused by previously known CLO. Unlike other CLO, CPT localized in the connective tissue rather than in the epithelial cells and formed smaller clumps of intracellular bacteria that stained dark blue with hematoxylin. On the other hand, typical myxospores of the genus Henneguya with tails were observed in the gill sections. Infection with Henneguya sp. resulted in extensive destruction of the gill filaments, most likely leading to respiratory distress. Due to the frequency of co-infections and the unavailability of culture methods for CLO and Henneguya sp, it was difficult to determine which pathogens were directly responsible for the associated mortality. However, co-infections may increase the negative impact on the host and the severity of the disease. Given the commercial importance of the snakeskin gourami and its significant aquaculture potential, the findings of this study are important for further studies on disease prevention. Spatiotemporal transcriptomic atlas of mouse organogenesis using DNA nanoball patterned arrays Spatially resolved transcriptomic technologies are promising tools to study cell fate decisions in a physical microenvironment, which is fundamental for enhancing our knowledge of mammalian development. However, the imbalance between resolution, transcript capture and field of view of current methodologies precludes their systematic application to analyze relatively large and three-dimensional mid- and late-gestation mammalian embryos. Here, we combined DNA nanoball (DNB) patterned arrays and tissue RNA capture to create SpaTial Enhanced REsolution Omics-sequencing (Stereo-seq). This approach allows transcriptomic profiling of large histological sections with high resolution and sensitivity. We have applied Stereo-seq to study the kinetics and directionality of transcriptional variation in a time course of mouse organogenesis. We used this information to gain insight into the molecular basis of regional specification, neuronal migration and differentiation in the developing brain. Furthermore, we mapped the expression of a panel of developmental disease-related loci on our global transcriptomic maps to define the spatiotemporal windows of tissue vulnerability. Our panoramic atlas constitutes an essential resource to investigate longstanding questions concerning normal and abnormal mammalian development.