Glypiated FGF directs contact-mediated bidirectional signaling to self-regulate tissue-specific dispersion How signaling proteins generate a multitude of information to organize tissue patterns is critical to understanding morphogenesis. In Drosophila, FGF produced in wing-disc cells regulates the development of the disc-associated air-sac-primordium (ASP). Here, we show that FGF is Glycosylphosphatidylinositol-anchored to the producing cell surface and that this modification both inhibits free FGF secretion and activates target-specific bidirectional FGF-FGFR signaling through cytonemes. FGF-source and recipient ASP cells extend cytonemes that present FGF and FGFR on their surfaces and reciprocally recognize each other over distance by contacting through CAM-like FGF-FGFR binding. Contact-mediated FGF-FGFR binding induces bidirectional signaling, which, in turn, promotes ASP and source cells to polarize cytonemes toward each other and reinforce signaling contacts. Subsequent un-anchoring of FGFR-bound-FGF from the source cell membrane dissociates cytoneme contacts and delivers FGF target-specifically to ASP cytonemes for paracrine functions. Thus, GPI-anchored FGF organizes both source and recipient cells and self-regulates its cytoneme-mediated tissue-specific dispersion and signaling.
Reuse of a wing venation gene-regulatory network in patterning the eyespot rings of butterflies Novel organismal traits might reuse ancestral gene-regulatory networks (GRNs) in their development, but data supporting this mechanism are still sparse. Here we propose that a conserved venation gene network, used in patterning insect wings, has been recruited to pattern the sharp and distinct rings of color in butterfly eyespot development. This network involves Decapentaplegic (Dpp), acting as a central morphogen, activating optix and spalt at different concentration thresholds. It also involves Spalt cross-regulating optix resulting in the formation of a sharp boundary between the two genes. First, using in situ hybridizations, we show that dpp is expressed in the eyespot centers during the early pupal stage and that CRISPR-Cas9 embryonic knockouts of dpp, optix, and spalt result in the loss of eyespots, loss of the orange ring, and loss of the black scales, respectively. Furthermore, using RNAi during the pupal stage, we show that dpp down-regulation results in the complete or partial loss of eyespots. Using CRISPR disruptions of optix or spalt, followed by immunostainings, we show that Spalt represses optix in the central region of the eyespot, limiting optix expression to a more peripheral ring. We further illustrate that optix disruptions alter scale pigmentation and scale ultrastructure as they transform orange to brown scales. Altogether, the regulatory interactions between dpp, spalt and optix mimic those found in earlier larval developmental stages, in the anterior compartment of an insects wing, where these genes play a role in venation patterning. The network similarities suggest that the venation GRN was co-opted to aid in the differentiation of the eyespot rings.
High performance single-cell gene regulatory network inference at scale: The Inferelator 3.0 MotivationGene regulatory networks define regulatory relationships between transcription factors and target genes within a biological system, and reconstructing them is essential for understanding cellular growth and function. Methods for inferring and reconstructing networks from genomics data have evolved rapidly over the last decade in response to advances in sequencing technology and machine learning. The scale of data collection has increased dramatically; the largest genome-wide gene expression datasets have grown from thousands of measurements to millions of single cells, and new technologies are on the horizon to increase to tens of millions of cells and above.

ResultsIn this work, we present the Inferelator 3.0, which has been significantly updated to integrate data from distinct cell types to learn context-specific regulatory networks and aggregate them into a shared regulatory network, while retaining the functionality of the previous versions. The Inferelator is able to integrate the largest single-cell datasets and learn cell-type specific gene regulatory networks. Compared to other network inference methods, the Inferelator learns new and informative Saccharomyces cerevisiae networks from single-cell gene expression data, measured by recovery of a known gold standard. We demonstrate its scaling capabilities by learning networks for multiple distinct neuronal and glial cell types in the developing Mus musculus brain at E18 from a large (1.3 million) single-cell gene expression dataset with paired single-cell chromatin accessibility data.

AvailabilityThe inferelator software is available on GitHub (https://github.com/flatironinstitute/inferelator) under the MIT license and has been released as python packages with associated documentation (https://inferelator.readthedocs.io/). Multi-omic analysis of the Arabidopsis clock activator mutant rve 4 6 8 reveals connections to carbohydrate metabolism and proteasome regulation Plants are able to sense changes in their light environments, such as the onset of day and night, as well as anticipate these changes in order to adapt and survive. Central to this ability is the plant circadian clock, a molecular circuit that precisely orchestrates plant cell processes over the course of a day. REVEILLE proteins (RVEs) are recently discovered members of the plant circadian circuitry that activate the evening complex and PRR genes to maintain regular circadian oscillation. The RVE 8 protein and its two homologs, RVE 4 and 6, have been shown to limit the length of the circadian period, with rve 4 6 8 triple-knockout plants possessing an elongated period along with increased leaf surface area, biomass, cell size and delayed flowering relative to wild-type Col-0 plants. Here, using a multi-omics approach consisting of phenomics, transcriptomics, proteomics, and metabolomics we draw novel connections between RVE8-like proteins and a number of core plant cell processes. In particular, we reveal that loss of RVE8-like proteins results in altered carbohydrate, organic acid and lipid metabolism, including a starch excess phenotype at dawn. We further demonstrate that rve 4 6 8 plants have lower levels of 20S proteasome subunits and possess significantly reduced proteasome activity, potentially explaining the increase in cell-size observed in RVE8-like mutants. Overall, this robust, multi-omic dataset, provides substantial new insights into the far reaching impact RVE8-like proteins have on the diel plant cell environment. CEP162 deficiency causes human retinal degeneration and reveals a dual role inciliogenesis and neurogenesis Defects in primary or motile cilia result in a variety of human pathologies, and retinal degeneration is frequently associated with these so-called ciliopathies. We show that homozygosity for a truncating variant in CEP162, a centrosome and microtubule-associated protein required for transition zone (TZ) assembly during ciliogenesis and neuronal differentiation in the retina, causes late-onset retinitis pigmentosa in 2 unrelated families. The mutant CEP162-E646R*5 protein is expressed and properly localized to the mitotic spindle but missing from the basal body in primary and photoreceptor cilia. This impairs recruitment of TZ components to the basal body and corresponds to complete loss of CEP162 function at the ciliary compartment, reflected by delayed formation of dysmorphic cilia. In contrast, rescue of increased cell death in the developing mouse retina after shRNA knockdown of Cep162 by expression of CEP162-E646R*5 indicates that the mutant retains its role for retinal neurogenesis. Human retinal degeneration thus results from specific loss of ciliary CEP162 function.

Graphical abstract

O_FIG O_LINKSMALLFIG WIDTH=134 HEIGHT=200 SRC="FIGDIR/small/469779v2_ufig1.gif" ALT="Figure 1">
View larger version (50K):
org.highwire.dtl.DTLVardef@83457eorg.highwire.dtl.DTLVardef@47775corg.highwire.dtl.DTLVardef@1ed643eorg.highwire.dtl.DTLVardef@5032d2_HPS_FORMAT_FIGEXP  M_FIG C_FIG Goal-centered representations in the human hippocampus Recent work in cognitive and systems neuroscience has suggested that the hippocampus might support planning, imagination, and navigation by forming "cognitive maps" that capture the structure of physical spaces, tasks, and situations. Critically, navigation involves planning within a context and disambiguating similar contexts to reach a goal. We examined hippocampal activity patterns in humans during a goal-directed navigation task to examine how contextual and goal information are incorporated in the construction and execution of navigational plans. Results demonstrate that, during planning, the hippocampus carries a context-specific representation of a future goal. Importantly, this effect could not be explained by stimulus or spatial information alone. During navigation, we observed reinstatement of activity patterns in the hippocampus ahead of participants required actions, which was strongest for behaviorally relevant points in the sequence. These results suggest that, rather than simply representing overlapping associations, hippocampal activity patterns are powerfully shaped by context and goals. Blebbistatin protects iPSC-CMs from hypercontraction and facilitates automated patch-clamp based electrophysiological study Recently, there have been great advances in cardiovascular channelopathy modeling and drug safety pharmacology using human induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs). The automated patch-clamp (APC) technique overcomes the disadvantages of the manual patch-clamp (MPC) technique, which is labor intensive and gives low output. However, the application of the APC platform is still limited in iPSC-CM based research, due to the difficulty in maintaining the high quality of single iPSC-CMs during dissociation and recording. In this study, we improved the method for single iPSC-CM preparation by applying 2.5 M blebbistatin (BB, an excitation-contraction coupling uncoupler) throughout APC procedures (dissociation, filtration, storage, and recording). Under non-BB buffered condition, iPSC-CMs in suspension showed a severe bleb-like morphology. However, BB-supplement led to significant improvements in morphology and INa recording, and we even obtained several CMs that showed spontaneous action potentials with typical morphology. Furthermore, APC faithfully recapitulated the single-cell electrophysiological phenotypes of iPSC-CMs derived from Brugada syndrome patients, as detected with MPC. Our study indicates that APC is capable of replacing MPC in the modeling of cardiac channelopathies using human iPSC-CMs by providing high-quality data with higher throughput. TeamTree analysis: a new approach to evaluate scientific production Advances in science and technology depend on the work of research teams and the publication of results through peer-reviewed articles representing a growing socio-economic resource. Current methods to mine the scientific literature regarding a field of interest focus on content, but the workforce credited by authorship remains largely unexplored, and appropriate measures of scientific production are debated. Here, a new bibliometric approach named TeamTree analysis is introduced that visualizes the development and composition of the workforce driving a field. A new citation-independent measure that scales with the H index estimates impact based on publication record, genealogical ties and collaborative connections. This author-centered approach complements existing tools to mine the scientific literature and to evaluate research across disciplines. Effects of substrate on shrimp growth, water quality and bacterial community in the biofloc system nursing Penaeus vannamei postlarvae This study aimed to investigate the effects of substrate on water quality, shrimp growth and bacterial community in the biofloc system with a salinity of 5{per thousand}. Two treatments, biofloc system with (sB) or without (nB) addition of elastic solid packing filler (nylon) as substrate, were set up. Penaeus vannamei postlarvae (PL, [~] stage 15) were stocked at a density of 4000 PL m-3 in each treatment with triplicates for a 28-days culture experiment, taking glucose as carbon source (C:N 15:1). Results showed that the survival rate (96.3{+/-}3.6%), FCR (0.76{+/-}0.06) and productivity (1.54{+/-}0.12 kg m-3) in sB treatment were significantly better than those in nB treatment (81.0{+/-}7.1%, 0.98{+/-}0.08 and 1.14{+/-}0.09 kg m-3, P<0.05). All water parameters were in the recommended ranges. Substrate showed significant effect on TAN, TSS, turbidity, biofloc volume, pH and carbonate alkalinity (P < 0.05). Actinobacteria (4.0-22.7%), Bacteroidetes (10.4-33.5%), Firmicutes (0.2-11.2%), Planctomycetes (4.0-14.9%) and Proteobacteria (29.4-59.0%) were the most dominant phyla for both treatments. However, the bacterial community in sB treatment showed to be significantly different from that in nB treatment (Jaccard distance 0.94{+/-}0.01, P=0.001). Substrate showed significant effects on Shannon, Heip, Pielou and Simpson index, as well as relative abundances of Actinobacteria, Bacteroidetes and Planctomycetes (P < 0.05). The results suggested that addition of substrate affected the shrimp growth, water quality and bacterial community in the biofloc system nursing P. vannamei PL with a 5{per thousand} salinity. Environment-dependence of the expression of mutational load and species' range limits Theoretical and empirical research on the causes of species range limits suggests the contribution of several intrinsic and extrinsic factors, with potentially complex interactions among them. An intrinsic factor proposed by recent theory is mutational load increasing towards range edges because of genetic drift. Furthermore, environmental quality may decline toward range edges and enhance the expression of load. Here we tested whether the expression of mutational load associated with range limits in the North American plant Arabidopsis lyrata is enhanced under stressful environmental conditions by comparing the performance of within- versus between-population crosses at common garden sites across the species distribution and beyond. Heterosis, reflecting the expression of load, increased with heightened estimates of genomic load and with environmental stress caused by warming, but the interaction was not significant. We conclude that range-edge populations suffer from a twofold genetic Allee effect caused by increased mutational load and stress-dependent load linked to general heterozygote deficiency, but no synergistic effect between them. Embryonic exposure to cannabidiol disrupts active neural circuits, an effect increased by Δ9-tetrahydrocannabinol and involving CB1R and CB2R in zebrafish In the light of legislative changes and the widespread use of cannabis as a recreational and medicinal drug, delayed effects of cannabis upon brief exposure during embryonic development are of high interest as early pregnancies often go undetected. Here, zebrafish embryos were exposed to cannabidiol (CBD) and {Delta}-9-tetrahydrocannabinol l (THC) until the end of gastrulation (1-10 hours post-fertilization) and analyzed later in development (4-5 days post-fertilization). In order to measure neural activity, we implemented CaMPARI (Calcium-Modulated Photoactivatable Ratiometric Integrator) and optimized the protocol for a 96-well format complemented by locomotor analysis. Our results revealed that neural activity was decreased by CBD more than THC. At higher doses, both cannabinoids could dramatically reduce neural activity and locomotor activity. Interestingly, the decrease was more pronounced when CBD and THC were combined. At the receptor level, CBD-mediated reduction of locomotor activity was partially prevented using cannabinoid type 1 and 2 receptor inhibitors. Overall, we report that CBD toxicity occurs via two cannabinoid receptors and is synergistically enhanced by THC exposure to negatively impact neural activity late in larval development. Future studies are warranted to reveal other cannabinoids and receptors involved in this pathway to understand the subsequent health implications of cannabis consumption on fetal development. ModEst - Precise estimation of genome size from NGS data Precise estimates of genome sizes are important parameters for both theoretical and practical biodiversity genomics. We present here a fast, easy-to-implement and precise method to estimate genome size from the number of bases sequenced and the mean sequence coverage. To estimate the latter, we take advantage of the fact that a precise estimation of the Poisson distribution parameter lambda is possible from truncated data, restricted to the part of the coverage distribution representing the true underlying distribution. With simulations we could show that reasonable genome size estimates can be gained even from low-coverage (10X), highly discontinuous genome drafts. Comparison of estimates from a wide range of taxa and sequencing strategies with flow-cytometry estimates of the same individuals showed a very good fit and suggested that both methods yield comparable, interchangeable results. Collective polymerase dynamics emerge from DNA supercoiling during transcription All biological processes ultimately come from physical interactions. The mechanical properties of DNA play a critical role in transcription. RNA polymerase can over or under twist DNA (referred to as DNA supercoiling) when it moves along a gene resulting in mechanical stresses in DNA that impact its own motion and that of other polymerases. For example, when enough supercoiling accumulates, an isolated polymerase halts and transcription stops. DNA supercoiling can also mediate non-local interactions between polymerases that shape gene expression fluctuations. Here, we construct a comprehensive model of transcription that captures how RNA polymerase motion changes the degree of DNA supercoiling which in turn feeds back into the rate at which polymerases are recruited and move along the DNA. Surprisingly, our model predicts that a group of three or more polymerases move together at a constant velocity and sustain their motion (forming what we call a polymeton) whereas one or two polymerases would have halted. We further show that accounting for the impact of DNA supercoiling on both RNA polymerase recruitment and velocity recapitulates empirical observations of gene expression fluctuations. Finally, we propose a mechanical toggle switch whereby interactions between genes are mediated by DNA twisting as opposed to proteins. Understanding the mechanical regulation of gene expression provides new insights into how endogenous genes can interact and informs the design of new forms of engineered interactions.

PACS numbers: The Functional Nanopore Screen: A Versatile High-throughput Assay to Study and Engineer Protein Nanopores in Escherichia coli Nanopores comprise a versatile class of membrane proteins that carry out a range of key physiological functions and are increasingly developed for different biotechnological applications. Yet, a capacity to study and engineer protein nanopores by combinatorial means has so far been hampered by a lack of suitable assays that combine sufficient experimental resolution with throughput. Addressing this technological gap, the Functional Nanopore (FuN) screen now provides a quantitative and dynamic read-out of nanopore assembly and function in E. coli. The assay is based on genetically-encoded fluorescent protein sensors that resolve the nanopore-dependent influx of Ca2+ across the inner membrane of E. coli. Illustrating its versatile capacity, the FuN screen is first applied to dissect the molecular features that underlie the assembly and stability of nanopores formed by the S2168 holin. In a subsequent step, nanopores are engineered by recombining the transmembrane module of S2168 with different ring-shaped oligomeric protein structures that feature defined hexa-, hepta- and octameric geometries. Library screening highlights substantial plasticity in the ability of the S2168 transmembrane module to oligomerize in alternative geometries while the functional properties of the resultant nanopores can be fine-tuned through the identity of the connecting linkers. Overall, the FuN screen is anticipated to facilitate both fundamental studies and complex nanopore engineering endeavors with many potential applications in biomedicine, biotechnology and synthetic biology. Glucose inhibits haemostasis and accelerates diet-induced hyperlipidaemia in zebrafish larvae Hyperglycaemia damages the microvasculature in part through the reduced recruitment of immune cells and interference with platelet signalling, leading to poor wound healing and accelerated lipid deposition in mammals. We investigated the utility of zebrafish larvae to model the effect of glucose on neutrophil and macrophage recruitment to a tail wound, wound-induced haemostasis, and chicken egg yolk feed challenge-induced hyperlipidaemia by supplementing larvae with exogenous glucose by immersion or injection. Neither method of glucose supplementation affected the recruitment of neutrophils and macrophages following tail transection. Glucose injection reduced thrombocyte retention and fibrin plug formation while only thrombocyte retention was reduced by glucose immersion following tail transection. We observed accelerated lipid accumulation in glucose-injected larvae challenged with high fat chicken egg yolk feeding. Our study identifies conserved and divergent effects of high glucose on inflammation, haemostasis, and hyperlipidaemia in zebrafish larvae compared to mammals. Single-dose AAV-based vaccine induces a high level of neutralizing antibodies against SARS-CoV-2 in rhesus macaques Coronavirus disease 2019 (COVID-19), which is triggered by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, continues to threaten global public health. Developing a vaccine that only requires single immunization but provides long-term protection for the prevention and control of COVID-19 is important. Here, we developed an adeno-associated virus (AAV)-based vaccine expressing a stable receptor-binding domain (SRBD) protein. The vaccine requires only a single shot but provides effective neutralizing antibodies (NAbs) over 598 days in rhesus macaques (Macaca mulatta). Importantly, our results showed that the NAbs were kept in high level and long lasting against authentic wild-type SARS-CoV-2, Beta, Delta and Omicron variants using plaque reduction neutralization test. Of note, although we detected pre-existing AAV2/9 antibodies before immunization, the vaccine still induced high and effective NAbs against COVID-19 in rhesus macaques. AAV-SRBD immune serum also efficiently inhibited the binding of ACE2 with RBD in the SARS-CoV-2 B.1.1.7 (Alpha), B.1.351 (Beta), P.1/P.2 (Gamma), B.1.617.2 (Delta), B.1.617.1/3(Kappa), and C.37 (Lambda) variants. Thus, these data suggest that the vaccine has great potential to prevent the spread of SARS-CoV-2. 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. The addition of FAIMS Increases Targeted Proteomics Sensitivity from FFPE Tumor Biopsies Mass spectrometry-based targeted proteomics allows objective protein quantitation of clinical biomarkers from a single section of formalin-fixed, paraffin-embedded (FFPE) tumor tissue biopsies. We combined high-field asymmetric waveform ion mobility spectrometry (FAIMS) and parallel reaction monitoring (PRM) to increase assay sensitivity. The modular nature of the FAIMS source allowed direct comparison of the performance of FAIMS-PRM to PRM. Limits of quantitation were determined by spiking synthetic peptides into a human spleen matrix. In addition, 20 clinical samples were analyzed using FAIMS-PRM and the quantitation of HER2 was compared with that obtained with the Ventana immunohistochemistry assay. FAIMS-PRM improved the overall signal-to-noise ratio over that from PRM and increased assay sensitivity in FFPE tissue analysis for four (HER2, EGFR, cMET, and KRAS) of five proteins of clinical interest. FAIMS-PRM enabled sensitive quantitation of basal HER2 expression in breast cancer samples classified as HER2 negative by immunohistochemistry. Furthermore, we determined the degree of FAIMS-dependent background reduction and showed that this correlated with an improved lower limit of quantitation with FAIMS. FAIMS-PRM is anticipated to benefit clinical trials in which multiple biomarker questions must be addressed and the availability of tumor biopsy samples is limited. Multiple lines of evidence combine to understand tree species coexistence Understanding drivers of species coexistence is a central challenge in ecology (1, 2). Coexistence cannot be observed directly, and while species co-occurrence in time and space is necessary for coexistence, it is not sufficient to prove coexistence (3, 4). Species exclusion from a region is potentially observable, but can take decades to occur, and still might just occur stochastically (5, 6). Thus, ecologists generally use theory to identify indirect observations that are indicative of mechanisms driving coexistence or exclusion (7-9). Various methods have been developed to indirectly infer coexistence, each of which requires different data, none of which are usually conclusive on their own, and are not typically combined (e.g. (9-12)). Here, we combine three different approaches which are normally used in isolation to build a case about a lack of coexistence of multiple hardwood species. Importantly, though each method has its flaws, we find that all three point to the same conclusion which strengthens the case for a lack of coexistence. First, in an experimental planting of three mature tree species we found no relationship between productivity and species diversity, which has been hypothesized to occur due to a lack of niche differences among species. Second, we used modern coexistence theory to calculate niche and fitness differences for each pair of species, which confirmed the high niche overlap among species and competitive dominance of one species. Third, we used the United States Department of Agriculture Forest Inventory and Analysis data to examine co-occurrence patterns of our species across thousands of natural forest stands and found that, these three species were distributed randomly throughout the USA, and do not significantly co-occur as would be expected if coexistence occurred naturally. Given that these independent methods agree; we take this as strong evidence about a lack of coexistence. Ustilago maydis telomere protein Pot1 harbors an extra N-terminal OB fold and regulates homology-directed DNA repair factors in a dichotomous and context-dependent manner The telomere G-strand binding protein Pot1 plays multifaceted roles in telomere maintenance and protection. We examined the structure and activities of Pot1 in Ustilago maydis, a fungal model that recapitulates key features of mammalian telomere regulation. Compared to the well-characterized primate and fission yeast Pot1 orthologs, UmPot1 harbors an extra N-terminal OB-fold domain (OB-N), which was recently shown to be present in most metazoans. UmPot1 binds directly to Rad51 and regulates the latters strand exchange activity. Deleting the OB-N domain, which is implicated in Rad51-binding, caused telomere shortening, suggesting that Pot1-Rad51 interaction facilitates telomere replication. Depleting Pot1 through transcriptional repression triggered growth arrest as well as rampant recombination, leading to multiple telomere aberrations. In addition, telomere repeat RNAs transcribed from both the G- and C-strand were dramatically up-regulated, and this was accompanied by elevated levels of telomere RNA-DNA hybrids. Telomere abnormalities of pot1-deficient cells were suppressed, and cell viability was restored by the deletion of genes encoding Rad51 or Brh2 (the BRCA2 ortholog), indicating that homology-directed repair (HDR) proteins are key mediators of telomere aberrations and cellular toxicity. Together, these observations underscore the complex physical and functional interactions between Pot1 and DNA repair factors, leading to context-dependent and dichotomous effects of HDR proteins on telomere maintenance and protection. Non-Newtonian Blood Rheology Impacts Left Atrial Stasis in Patient-Specific Simulations The lack of mechanically effective contraction of the left atrium (LA) during atrial fibrillation (AF) disturbs blood flow, increasing the risk of thrombosis and ischemic stroke. Thrombosis is most likely in the left atrial appendage (LAA), a small narrow sac where blood is prone to stagnate. Slow flow promotes the formation of erythrocyte aggregates in the LAA, also known as rouleaux, causing viscosity gradients that are usually disregarded in patient-specific simulations. To evaluate these non-Newtonian effects, we built atrial models derived from 4D computed tomography scans of patients and carried out computational fluid dynamics simulations using the Carreau-Yasuda constitutive relation. We examined six patients, three of whom had AF and LAA thrombosis or a history of transient ischemic attacks (TIAs). We modeled the effects of hematocrit and rouleaux formation kinetics by varying the parameterization of the Carreau-Yasuda relation and modulating non-Newtonian viscosity changes based on residence time. Comparing non-Newtonian and Newtonian simulations indicates that slow flow in the LAA increases blood viscosity, altering secondary swirling flows and intensifying blood stasis. While some of these effects can be subtle when examined using instantaneous metrics like shear rate or kinetic energy, they are manifested in the blood residence time, which accumulates over multiple heartbeats. Our data also reveal that LAA blood stasis worsens when hematocrit increases, offering a potential new mechanism for the clinically reported correlation between hematocrit and stroke incidence. In summary, we submit that hematocrit-dependent non-Newtonian blood rheology should be considered in calculating patient-specific blood stasis indices by computational fluid dynamics. Comparative Analysis of T Cell Spatial Proteomics and the Influence of HIV Expression As systems biology approaches to virology have become more tractable, highly studied viruses such as HIV can now be analyzed in new, unbiased ways, including spatial proteomics. We employed here a differential centrifugation protocol to fractionate Jurkat T cells for proteomic analysis by mass spectrometry; these cells contain inducible HIV-1 genomes, enabling us to look for changes in the spatial proteome induced by viral gene expression. Using these proteomics data, we evaluated the merits of several reported machine learning pipelines for classification of the spatial proteome and identification of protein translocations. From these analyses we found that classifier performance in this system was organelle-dependent, with Bayesian t-augmented Gaussian mixture modeling outperforming support vector machine (SVM) learning for mitochondrial and ER proteins, but underperforming on cytosolic, nuclear, and plasma membrane proteins by QSep analysis. We also observed a generally higher performance for protein translocation identification using a Bayesian model, BANDLE, on SVM-classified data. Comparative BANDLE analysis of cells induced to express the wild-type viral genome vs. cells induced to express a genome unable to express the accessory protein Nef identified known Nef-dependent interactors such as TCR signaling components and coatomer complex. Lastly, we found that SVM classification showed higher consistency and was less sensitive to HIV-dependent noise. These findings illustrate important considerations for studies of the spatial proteome following viral infection or viral gene expression and provide a reference for future studies of HIV-gene-dropout viruses. METAbolomics data Balancing with Over-sampling Algorithms (META-BOA): an online resource for addressing class imbalance MotivationClass imbalance, or unequal sample sizes between classes, is an increasing concern in machine learning for metabolomic and lipidomic data mining, which can result in overfitting for the over-represented class. Numerous methods have been developed for handling class imbalance, but they are not readily accessible to users with limited computational experience. Moreover, there is no resource that enables users to easily evaluate the effect of different over-sampling algorithms.

ResultsMETAbolomics data Balancing with Over-sampling Algorithms (META-BOA) is a web-based application that enables users to select between four different methods for class balancing, followed by data visualization and classification of the sample to observe the augmentation effects. META-BOA outputs a newly balanced dataset, generating additional samples in the minority class, according to the users choice of Synthetic Minority Over-sampling Technique (SMOTE), Borderline-SMOTE (BSMOTE), Adaptive Synthetic (ADASYN), or Random Over-Sampling Examples (ROSE). META-BOA further displays both principal component analysis (PCA) and t-distributed stochastic neighbor embedding (t-SNE) visualization of data pre- and post-over-sampling. Random forest classification is utilized to compare sample classification in both the original and balanced datasets, enabling users to select the most appropriate method for their analyses.

Availability and implementationMETA-BOA is available at https://complimet.ca/meta-boa.

Supplementary InformationSupplementary material is available at Bioinformatics online. Long term analysis of social structure: evidence of age-based consistent associations in male Alpine ibex Despite its recognized importance for understanding the evolution of animal sociality as well as for conservation, long term analysis of social networks of animal populations is still relatively uncommon. We investigated social network dynamics in males of a gregarious mountain ungulate (Alpine ibex, Capra ibex) over ten years focusing on groups, sub-groups and individuals, exploring the dynamics of sociality over different scales. Despite the social structure changing between seasons, the Alpine ibex population was highly cohesive: fission-fusion dynamics lead almost every male in the population to associate with each other male at least once. Nevertheless, we found that male Alpine ibex showed preferential associations that were maintained across seasons and years. Age seemed to be the most important factor driving preferential associations while other characteristics, such as social status, appeared less crucial. We also found that centrality measures were influenced by age and were also related to individual physical condition. The multi-scale and long-term frame of our study helped us show that ecological constrains, such as resource availability, may play a role in shaping associations in a gregarious species, but they cannot solely explain sociality and preferential association that are likely also to be driven by life-history linked physiological and social needs. Our results highlight the importance of long-term studies based on individually recognizable subjects to help us build on our understanding of the evolution of animal sociality. Structural basis for Nirmatrelvir in vitro efficacy against SARS-CoV-2 variants The COVID-19 pandemic continues to be a public health threat with emerging variants of SARS-CoV-2. Nirmatrelvir (PF-07321332) is a reversible, covalent inhibitor targeting the main protease (Mpro) of SARS-CoV-2 and the active protease inhibitor in PAXLOVID (nirmatrelvir tablets and ritonavir tablets). We evaluated the in vitro catalytic activity and in vitro potency of nirmatrelvir against the main protease (Mpro) of prevalent variants of concern (VOC) or variants of interest (VOI): Alpha (, B.1.1.7), Beta ({beta}, B.1.351), Delta ({delta}, B1.617.2), Gamma ({gamma}, P.1), Lambda ({lambda}, B.1.1.1.37/C37), Omicron (o, B.1.1.529) as well as the original Washington or wildtype strain. These VOC/VOI carry prevalent mutations at varying frequencies in the Mpro specifically for: , {beta},{gamma} (K90R), {lambda} (G15S) and o (P132H). In vitro biochemical enzymatic assay characterization of the enzyme kinetics of the mutant Mpros demonstrate that they are catalytically comparable to wildtype. Nirmatrelvir has similar potency against each mutant Mpro including P132H that is observed in the Omicron variant with a Ki of 0.635 nM as compared to a Ki of 0.933nM for wildtype. The molecular basis for these observations were provided by solution-phase structural dynamics and structural determination of nirmatrelvir bound to the o, {lambda} and {beta} Mpro at 1.63 - 2.09 [A] resolution. These in vitro data suggest that PAXLOVID has the potential to maintain plasma concentrations of nirmatrelvir many-fold times higher than the amount required to stop the SARS-CoV-2 VOC/VOI, including Omicron, from replicating in cells (1). The virucidal effects of 405 nm visible light on SARS-CoV-2 and influenza A virus Germicidal potential of specific wavelengths within the electromagnetic spectrum is an area of growing interest. While ultra-violet (UV) based technologies have shown satisfactory virucidal potential, the photo-toxicity in humans coupled with UV associated polymer degradation limit its use in occupied spaces. Alternatively, longer wavelengths with less irradiation energy such as visible light (405 nm) have largely been explored in the context of bactericidal and fungicidal applications. Such studies indicated that 405 nm mediated inactivation is caused by the absorbance of porphyrins within the organism creating reactive oxygen species which result in free radical damage to its DNA and disruption of cellular functions. The virucidal potential of visible-light based technologies has been largely unexplored and speculated to be ineffective given the lack of porphyrins in viruses. The current study demonstrated increased susceptibility of lipid-enveloped respiratory pathogens of importance such as SARS-CoV-2 (causative agent of COVID-19) as well as the influenza A virus to 405nm, visible light in the absence of exogenous photosensitizers indicating a potential porphyrin-independent alternative mechanism of visible light mediated viral inactivation. These results were obtained using less than expected irradiance levels which are generally safe for humans and commercially achievable. Our results support further exploration of the use of visible light technology for the application of continuous decontamination in occupied areas within hospitals and/or infectious disease laboratories, specifically for the inactivation of respiratory pathogens such as SARS-CoV-2 and Influenza A.