Distinct repertoire development history during affinity maturation of antigen-responding VHH antibody Time-resolved antibody repertoire analysis is essential to elucidate the antibody maturation process during immunization. However, most analyses are merely snapshot comparisons of unimmunized and fully immunized animals. In this study, an alpaca was immunized with multiple antigen injections and its repertoire development was traced by high-throughput next-generation sequencing every week for 3 months. The use of alpaca-derived, single-domain antibody circumvents the problems associated with conventional antibodies including heavy-chain and light-chain matching and labor-intensive production. The sequences were processed into clusters and the antibodies were generated to determine whether the clusters included antigen-binding antibodies. For most of the antigen-responding clusters, repertoire development exhibited distinct time-dependent patterns compared with clusters unresponsive to antigens. Moreover, the sequences at the tips and root of the phylogenetic tree had strong and weak antigen-binding affinity, respectively. The foregoing characteristics could be exploited in the prediction of antigen-responding clusters and the physical properties of antibodies.

TeaserClusters of antigen-responding, single-domain antibodies exhibited distinct time-dependent repertoire development history
An AI-guided signature reveals the nature of the shared proximal pathways of host immune response in MIS-C and Kawasaki disease A significant surge in cases of multisystem inflammatory syndrome in children (MIS-C, also called Pediatric Inflammatory Multisystem Syndrome - PIMS) has been observed amidst the COVID-19 pandemic. MIS-C shares many clinical features with Kawasaki disease (KD), although clinical course and outcomes are divergent. We analyzed whole blood RNA sequences, serum cytokines, and formalin fixed heart tissues from these patients using a computational toolbox of two gene signatures, i.e., the 166-gene viral pandemic (ViP) signature, and its 20-gene severe (s)ViP subset that were developed in the context of SARS-CoV-2 infection and a 13-transcript signature previously demonstrated to be diagnostic for KD. Our analyses revealed that KD and MIS-C are on the same continuum of the host immune response as COVID-19. While both the pediatric syndromes converge upon an IL15/IL15RA-centric cytokine storm, suggestive of shared proximal pathways of immunopathogenesis, they diverge in other laboratory parameters and cardiac phenotypes. The ViP signatures also revealed unique targetable cytokine pathways in MIS-C, place MIS-C farther along in the spectrum in severity compared to KD and pinpoint key clinical (reduced cardiac function) and laboratory (thrombocytopenia and eosinopenia) parameters that can be useful to monitor severity.
The Ensemble of Gene Regulatory Networks at Mutation-Selection Balance The evolution of diverse phenotypes both involves and is constrained by molecular interaction networks. When these networks influence patterns of expression, we refer to them as gene regulatory networks (GRNs). Here, we develop a quasi-species model of GRN evolution. With this model, we prove that-across a broad spectrum of viability and mutation functions-the dynamics converge to a stationary distribution over GRNs. Next, we show from first principles how the frequency of GRNs at equilibrium will be proportional to each GRNs eigenvector centrality in the genotype network. Finally, we determine the structural characteristics of GRNs that are favored in response to a range of selective environments and mutational constraints. Our work connects GRN evolution to quasi-species models, and thus can provide a mechanistic explanation for the topology of GRNs experiencing various evolutionary forces. 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. Rho Kinase regulates neutrophil NET formation that is involved in UVB-induced skin inflammation ObjectiveUltraviolet B (UVB) is an important trigger of skin inflammation and lupus with leukocyte recruitment to inflamed skin. We recently reported the involvement of neutrophil NETosis in UVB-induced skin inflammation, and that NETotic nuclear envelope rupture is driven by PKC-mediated nuclear lamin B disassembly. To address the role of Actin cytoskeleton in NETosis, we investigated the effects of Rho kinase (ROCK) and its downstream actomyosin cytoskeletal networks on PKC nuclear translocation and NET formation, as well as their involvement in UVB-induced skin inflammation.

MethodsWe studied the dynamic changes of ROCK and actomyosin cytoskeletal networks during NETosis induction and their involvement in PKC nuclear translocation. Using mice with hematopoietic-specific ROCK1 deficiency, we investigated the effects of ROCK1 deficiency on NETosis, and its involvement in UVB-induced skin inflammation.

ResultsOur time course studies demonstrated the dynamic changes of actin polymerization and ROCK activation, support the role of actin cytoskeleton in nuclear translocation of cytosolic PKC in early stage of NETosis induction. Inhibition of actin polymerization or key molecules of the ROCK/MLCK/myosin pathway decreased PKC nuclear translocation and NET formation. Genetic deficiency of ROCK1, inhibited NETosis ex vivo and in vivo, decreased extracellular display of NET-associated IL-17A, TNF, IFN{gamma}, and IFN in inflamed skin, which were correlated with the ameliorated skin inflammation in UVB-irradiated mice with hematopoietic-specific ROCK1 deficiency.

ConclusionsROCK regulated NETosis through modulation of PKC nuclear translocation via actomyosin cytoskeletal networks in neutrophils. ROCK1 deficiency ameliorated UVB- induced skin inflammation by attenuation of NETosis and NET-associated cytokines. Cluster-specific gene markers enhance Shigella and Enteroinvasive Escherichia coli in silico serotyping Shigella and enteroinvasive Escherichia coli (EIEC) cause human bacillary dysentery with similar invasion mechanisms and share similar physiological, biochemical and genetic characteristics. The ability to differentiate Shigella and EIEC from each other is important for clinical diagnostic and epidemiologic investigations. The existing genetic signatures may not discriminate between Shigella and EIEC. However, phylogenetically, Shigella and EIEC strains are composed of multiple clusters and are different forms of E. coli. In this study, we identified 10 Shigella clusters, 7 EIEC clusters and 53 sporadic types of EIEC by examining over 17,000 publicly available Shigella/EIEC genomes. We compared Shigella and EIEC accessory genomes to identify the cluster-specific gene markers or marker sets for the 17 clusters and 53 sporadic types. The gene markers showed 99.63% accuracy and more than 97.02% specificity. In addition, we developed a freely available in silico serotyping pipeline named Shigella EIEC Cluster Enhanced Serotype Finder (ShigEiFinder) by incorporating the cluster-specific gene markers and established Shigella/EIEC serotype specific O antigen genes and modification genes into typing. ShigEiFinder can process either paired end Illumina sequencing reads or assembled genomes and almost perfectly differentiated Shigella from EIEC with 99.70% and 99.81% cluster assignment accuracy for the assembled genomes and mapped reads respectively. ShigEiFinder was able to serotype over 59 Shigella serotypes and 22 EIEC serotypes and provided a high specificity with 99.40% for assembled genomes and 99.38% for mapped reads for serotyping. The cluster markers and our new serotyping tool, ShigEiFinder (https://github.com/LanLab/ShigEiFinder), will be useful for epidemiologic and diagnostic investigations.

Impact statementThe differentiation of Shigella strains from enteroinvasive E. coli (EIEC) is important for clinical diagnosis and public health epidemiologic investigations. The similarities between Shigella and EIEC strains make this differentiation very difficult as both share common ancestries within E. coli. However, Shigella and EIEC are phylogenetically separated into multiple clusters, making high resolution separation using cluster specific genomic markers possible. In this study, we identified 17 Shigella or EIEC clusters including five that were newly identified through examination of over 17,000 publicly available Shigella and EIEC genomes. We further identified an individual or a set of cluster-specific gene markers for each cluster using comparative genomic analysis. These markers can then be used to classify isolates into clusters and were used to develop an in silico pipeline, ShigEiFinder (https://github.com/LanLab/ShigEiFinder) for accurate differentiation, cluster typing and serotyping of Shigella and EIEC from Illumina sequencing reads or assembled genomes. This study will have broad application from understanding the evolution of Shigella/EIEC to diagnosis and epidemiology.

Data summarySequencing data have been deposited at the National Center for Biotechnology Information under BioProject number PRJNA692536.

RepositoriesRaw sequence data are available from NCBI under the BioProject number PRJNA692536. Protective host-dependent antagonism among 	Pseudomonas in the Arabidopsis phyllosphere The plant microbiome is a rich biotic environment, comprising numerous taxa. The community structure of these colonizers is constrained by multiple factors, including host-microbe and microbe-microbe interactions, as well as the interplay between the two. While much can be learned from pairwise relationships between individual hosts and microbes, or individual microbes with themselves, the ensemble of interrelations between the host and microbial consortia may lead to different outcomes that are not easily predicted from the individual interactions. Their study can thus provide new insights into the complex relationship between plants and microbes. Of particular importance is how strain-specific such plant-microbe-microbe interactions are, and how they eventually affect plant health. Here, we test strain-level interactions in the phyllosphere between groups of co-existing commensal and pathogenic Pseudomonas among each other and with A. thaliana, by employing synthetic communities of genome-barcoded isolates. We found that commensal Pseudomonas prompted a host response leading to a selective inhibition of a specific pathogenic lineage, resulting in plant protection. The extent of plant protection, however, was dependent on plant genotype, indicating that these effects are host-mediated. There were similar genotype-specific effects on the microbe side, as we could pinpoint an individual Pseudomonas isolate as the predominant cause for this differential interaction. Collectively, our work highlights how within-species genetic differences on both the host and microbe side can have profound effects on host-microbe-microbe dynamics. The paradigm that we have established provides a platform for the study of host-dependent microbe-microbe competition and cooperation in the A. thaliana-Pseudomonas system. High-plex Multiomic Analysis in FFPE at Subcellular Level by Spatial Molecular Imaging The Spatial Molecular Imaging platform (CosMxTM SMI, NanoString, Seattle, WA) utilizes a high-plex in-situ imaging chemistry for both RNA and protein detection. This automated instrument provides 1000s of plex, at high sensitivity (1 to 2 copies/cell), very low error rate (0.0093 false calls/cell), and background ([~]0.04 counts/cell). The imaging system generates three-dimensional super-resolution localization of analytes (< 50 nm in the X-Y plane), at [~] 1 million cells per sample, four samples per run. Cell segmentation is antibody-based, compatible with FFPE samples. Multiomic data (980 RNAs, 108 proteins) were measured at subcellular resolution using FFPE tissues (non-small cell lung (NSCLC) and breast cancer) and allowed identification of over 18 distinct cell types, 10 unique tumor microenvironments, and 100 pairwise ligand-receptor interactions. Over 800,000 single cells and [~]260 million transcripts data are released into the public domain allowing extended data analysis by the entire spatial biology research community. Once-daily feeding is associated with better health in companion dogs: Results from the Dog Aging Project A variety of diets have been studied for possible anti-aging effects. In particular, studies of intermittent fasting and time-restricted feeding in laboratory rodents have found evidence of beneficial health outcomes. Companion dogs represent a unique opportunity to study diet in a large mammal that shares human environments. The Dog Aging Project has been collecting data on thousands of companion dogs of all different ages, sizes, and breeds since 2019. We leveraged this diverse cross-sectional dataset to investigate associations between feeding frequency and cognitive function (n = 10,474) as well as nine broad categories of health conditions (n = 24,238). Controlling for sex, age, breed, and other potential confounders, we found that dogs fed once daily rather than more frequently had lower mean scores on a cognitive dysfunction scale, and lower odds of having gastrointestinal, dental, orthopedic, kidney/urinary, and liver/pancreas disorders. Therefore, we find that once-daily feeding is associated with better health in multiple domains. Future research with longitudinal data can provide stronger evidence for a possible causal effect of feeding frequency on health in companion dogs. Kinetochore-mediated microtubule assembly and Augmin-dependent amplification drive k-fiber maturation in mammals Chromosome segregation in mammals relies on the maturation of a thick bundle of kinetochore-attached microtubules known as k-fibers. How k-fibers mature from initial kinetochore-microtubule attachments remains a fundamental question. By combining functional perturbations in Indian muntjac, a placental mammal with the lowest known chromosome number (n=3) and distinctively large kinetochores, with fixed- and subsecond live-cell super-resolution CH-STED nanoscopy, we uncovered the mechanism by which Augmin mediates k-fiber maturation. Augmin promoted kinetochore microtubule turnover by sustaining microtubule formation from kinetochores and poleward flux, regardless of pioneer centrosomal microtubules. Tracking of microtubule growth events within individual k-fibers revealed a wide angular dispersion, consistent with Augmin-mediated branched microtubule nucleation. Augmin depletion reduced the frequency of microtubule growth events within k-fibers and hampered normal repair after acute k-fiber injury by laser microsurgery. Our work underscores the contribution of microtubule formation from kinetochores and Augmin-mediated microtubule amplification for k-fiber maturation and spindle assembly in mammals. Missing the forest for the trees: Do seaweed ecosystems mitigate atmospheric CO2 emissions? Global seaweed carbon sequestration estimates are currently taken as the fraction of the net primary production (NPP) exported to the deep ocean. However, this perspective does not account for CO2 from the consumption of external subsidies. Here we clarify: i) the role of export relative to seaweed net ecosystem production (NEP) for a closed system and one more likely open to subsidies; ii) the importance of subsidies by compiling published estimates of NEP from seaweed-dominated ecosystems; and iii) discuss their impact on the global seaweed net carbon balance and other sequestration constraints as a mitigation service. Examples of seaweed NEP (n = 18) were sparse and variable. Nevertheless, the average NEP (-9.2mmol C m-2 day-1 SE {+/-} 11.6) suggested that seaweed ecosystems are a C source, becoming increasingly heterotrophic as their export is consumed. Critically, mitigation of greenhouse gas emissions was mixed relative to their replacement or baseline states, and where CO2 is supplied independently of organic metabolism and atmospheric exchange we caution a sole reliance on NEP or NPP. This will ensure a more accurate seaweed mitigation assessment, one that does exceed their capacity and is effective within a compliance and carbon trading scheme. Structural and mechanistic basis of reiterative transcription initiation Reiterative transcription initiation, observed at promoters that contain homopolymeric sequences at the transcription start site, generates RNA products having 5 sequences non-complementary to the DNA template. Here, using crystallography and cryo-EM to define structures, protein-DNA-photocrosslinking to map positions of RNAP leading and trailing edges relative to DNA, and single-molecule DNA nanomanipulation to assess RNAP-dependent DNA unwinding, we show that RNA extension in reiterative transcription initiation: (1) occurs without DNA scrunching, (2) involves a short, 2-3 bp RNA-DNA hybrid, and (3) generates RNA that exits RNAP through the portal by which scrunched nontemplate-strand DNA exits RNAP in standard transcription initiation. The results establish that, whereas RNA extension in standard transcription initiation proceeds through a scrunching mechanism, RNA extension in reiterative transcription initiation proceeds through a slippage mechanism, with slipping of RNA relative to DNA within a short RNA-DNA hybrid, and with extrusion of RNA from RNAP through an alternative RNA exit. 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. Transcriptional and Functional Activity of Hemangiosarcoma Support Bone Marrow Nurse Cell Ontogeny Hemangiosarcoma and angiosarcoma are soft tissue sarcomas of malignant blood vessel-forming cells in dogs and humans, respectively. These vasoformative sarcomas are aggressive and highly metastatic, with disorganized, irregular blood-filled vascular spaces. Our objective was to define molecular programs that support the niche, enabling progression of canine hemangiosarcoma and human angiosarcoma. Here, we show that the transcriptional landscape of canine hemangiosarcoma and human angiosarcoma included comparable angiogenic and inflammatory programs. Dog-in-mouse hemangiosarcoma xenografts recapitulated the vasoformative and highly angiogenic morphology and molecular characteristics of primary tumors. Blood vessels in the tumors were complex and disorganized, and they were lined by both donor and host cells, a trait that was not observed in xenografts from canine osteosarcoma and lymphoma. In some cases, the xenografted hemangiosarcoma cells created exuberant myeloid hyperplasia and gave rise to lymphoproliferative tumors of mouse origin. We did not uncover a definitive transmissible etiology, but our data indicate that transcriptional programs of hemangiosarcoma cells resemble those of hematopoietic nurse cells, and these malignant cells support expansion and differentiation of human hematopoietic progenitors. We conclude that canine hemangiosarcoma, and possibly human angiosarcoma, originate from nurse cells that make up the stromal bone marrow niche, and that these cells may also support the growth of hematopoietic tumors. Citation needed?Wikipedia bibliometrics during the first waveof 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. 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. Ribosome actuator via the SpoT/ppGpp pathway to mitigate gene overexpression burden Gene overexpression causes a burden on ribosomes, which decreases the expression of other genes and growth rate. Here, we engineer a genetic ribosome actuator that increases gene expression rate and growth rate by over 150% and 80%, respectively. The actuator expresses the hydrolysis domain of the SpoT enzyme (SpoTH) in a bacterial strain with high basal levels of ppGpp, thus increasing ribosome availability. We use this actuator in a feedforward control arrangement, wherein SpoTH is co-expressed with a gene of interest (GOI). Without SpoTH co-expression, activation of the GOI decreases the expression of a constitutive gene and growth rate by over 80% and 40%, respectively. By contrast, with SpoTH co-expression, the ribosome binding site (RBS) of SpoTH can be tuned to achieve inappreciable change of either the expression of the constitutive gene or growth rate when the GOI is activated to the same level. We finally employ the actuator to keep growth rate constant while expressing dCas9 to a level that would otherwise decrease growth rate by more than 40% due to toxicity. This actuator will be relevant in applications where high cellular protein concentrations are desired without growth defects or global reduction of gene expression. Time- and cell-dependent atypia and cell death are caused by progressive deficiency in DNA replication Pleiotropy caused by single-gene mutations is common and poorly understood. A zebrafish null mutant of DNA polymerase  subunit B, huli hutu (hht), evolves a complex pleiotropy associated with DNA damage and S phase arrest across multiple organ systems over 5-7 days, including nuclear atypia, a common cellular feature in human cancers and pre-cancers, in gastrointestinal organs, and nuclear fragmentation in the eye and brain. The pleiotropic pattern of hht phenotypes is explained by progressive loss of wild-type maternal pola2 function in homozygous mutant embryos whose pola2 mRNA becomes undetectable by 24 hours post-fertilization (hpf). Inhibition of DNA synthesis by aphidicolin or hydroxyurea in wild-type embryos from 24 hpf phenocopied the pleiotropic pattern of hht. These results are consistent with a model in which time-sensitive, reduced capacity for DNA synthesis results in cell death in fast-replicating cells, and nuclear atypia in tissues with fewer and larger cells. Patient-specific network connectivity combined with a next generation neural mass model to test clinical hypothesis of seizurepropagation Dynamics underlying epileptic seizures span multiple scales in space and time, therefore, understanding seizure mechanisms requires identifying the relations between seizure components within and across these scales, together with the analysis of their dynamical repertoire. In this view, mathematical models have been developed, ranging from single neuron to neural population.

In this study we consider a neural mass model able to exactly reproduce the dynamics of heterogeneous spiking neural networks. We combine the mathematical modelling with structural information from non-invasive brain imaging, thus building large-scale brain network models to explore emergent dynamics and test clinical hypothesis. We provide a comprehensive study on the effect of external drives on neuronal networks exhibiting multistability, in order to investigate the role played by the neuroanatomical connectivity matrices in shaping the emergent dynamics. In particular we systematically investigate the conditions under which the network displays a transition from a low activity regime to a high activity state, which we identify with a seizure-like event. This approach allows us to study the biophysical parameters and variables leading to multiple recruitment events at the network level. We further exploit topological network measures in order to explain the differences and the analogies among the subjects and their brain regions, in showing recruitment events at different parameter values.

We demonstrate, along the example of diffusion-weighted magnetic resonance imaging (MRI) connectomes of 20 healthy subjects and 15 epileptic patients, that individual variations in structural connectivity, when linked with mathematical dynamic models, have the capacity to explain changes in spatiotemporal organization of brain dynamics, as observed in network-based brain disorders. In particular, for epileptic patients, by means of the integration of the clinical hypotheses on the epileptogenic zone (EZ), i.e. the local network where highly synchronous seizures originate, we have identified the sequence of recruitment events and discussed their links with the topological properties of the specific connectomes. The predictions made on the basis of the implemented set of exact mean-field equations turn out to be in line with the clinical pre-surgical evaluation on recruited secondary networks. Receptor localization maximizes information acquisition in natural cell 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. Precartilage condensation during limb skeletogenesis occurs by tissue phase separation controlled by a bistable cell-state switch with suppressed oscillatory dynamics The tetrapod limb skeleton is initiated in unpatterned limb bud mesenchyme by the formation of precartilage condensations. Here, based on time-lapse videographic analysis of a forming condensation in a high-density culture of chicken limb bud mesenchyme, we observe a phase transition to a more fluidized state for cells within spatial compacted foci (protocondensations that will progress to condensations), as reflected in their spatial confinement, cell-substratum interaction and speed of motion. Previous work showed that galectin-8 and galectin-1A, two proteins of the galactoside-binding galectin family, are the earliest determinants of this process in the chicken limb bud, and that their interactions in forming skeletogenic patterns of condensations can be interpreted mathematically through a reaction-diffusion-adhesion framework. Based on this framework, we use an ordinary differential equation-based approach to analyze the core switching modality of the galectin reaction network and characterize the states of the network independent of the diffusive and adhesive arms of the patterning mechanism. We identify two steady states where the concentrations of both galectins are respectively, negligible, and very high. An explicit Lyapunov function shows that there are no periodic solutions. For sigmoidal galectin production terms, the model exhibits a bistable switch that arises from a monostable state via saddle-node bifurcation. Our model therefore predicts that the galectin network exists in low and high expression states separated in space or time without any intermediate states. This provides a causal basis for the observed outside vs. inside transition observed in the in vitro video data. We performed a quantitative analysis of the distribution of galectin-1A in cultures of condensing chick limb mesenchymal cells and found that the interior of the protocondensations had concentrations of this protein (compared to the immediate exterior) over and above that expected from its higher cell density, consistent with the models predictions. The galectin-based patterning network is thus suggested, on theoretical grounds, to incorporate a core switch independent of any spatial or temporal dynamics, that drives the chondrogenic cell state transition in limb skeletogenesis. Protein contact map prediction using multiple sequence alignment dropout and consistency learning for sequences with fewer homologs The prediction of protein contact map needs enough normalized number of effective sequence (Nf) in multiple sequence alignment (MSA). When Nf is small, the predicted contact maps are often not satisfactory. To solve this problem, we randomly selected a small part of sequence homologs for proteins with large Nf to generate MSAs with small Nf. From these MSAs, input features were generated and were passed through a consistency learning network, aiming to get the same results when using the features generated from the MSA with large Nf. The results showed that this method effectively improves the prediction accuracy of protein contact maps with small Nf. Elucidation of SARS-CoV-2 budding mechanisms through molecular dynamics simulations of M and E protein complexes SARS-CoV-2 and other coronaviruses pose major threats to global health, yet computational efforts to understand them have largely overlooked the process of budding, a key part of the coronavirus life cycle. When expressed together, coronavirus M and E proteins are sufficient to facilitate budding into the ER-Golgi intermediate compartment (ERGIC). To help elucidate budding, we ran atomistic molecular dynamics (MD) simulations using the Feig laboratorys refined structural models of the SARS-CoV-2 M protein dimer and E protein pentamer. Our MD simulations consisted of M protein dimers and E protein pentamers in patches of membrane. By examining where these proteins induced membrane curvature in silico, we obtained insights around how the budding process may occur. Multiple M protein dimers acted together to induce global membrane curvature through protein-lipid interactions while E protein pentamers kept the membrane planar. These results could eventually help guide development of antiviral therapeutics which inhibit coronavirus budding. Beta-cell specific Insr deletion promotes insulin hypersecretion and improves glucose tolerance prior to global insulin resistance Abstract

Insulin receptor (Insr) protein can be found at higher levels in pancreatic {beta}-cells than in most other tissues, but the consequences of {beta}-cell insulin resistance remain enigmatic. Ins1cre allele was used to delete Insr specifically in {beta}-cells of both female and male mice. Experimental mice were compared to Ins1cre-containing littermate controls at multiple ages and on multiple diets. RNA-seq of purified recombined {beta}-cells revealed transcriptomic consequences of Insr loss, which differed between female and male mice. Action potential and calcium oscillation frequencies were increased in Insr knockout {beta}- cells from female, but not male mice, whereas only male {beta}InsrKO mice had reduced ATP-coupled oxygen consumption rate and reduced expression of genes involved in ATP synthesis. Female {beta}InsrKO and {beta}InsrHET mice exhibited elevated insulin release in perifusion experiments, during hyperglycemic clamps, and following i.p. glucose challenge. Deletion of Insr did not alter {beta}-cell area up to 9 months of age, nor did it impair hyperglycemia-induced proliferation. Based on our data, we adapted a mathematical model to include {beta}-cell insulin resistance, which predicted that {beta}-cell Insr knockout would improve glucose tolerance depending on the degree of whole-body insulin resistance. Indeed, glucose tolerance was significantly improved in female {beta}InsrKO and {beta}InsrHET mice when compared to controls at 9, 21 and 39 weeks, and also in insulin-sensitive 4-week old males. We did not observe improved glucose tolerance in older male mice or in high fat diet-fed mice, corroborating the prediction that global insulin resistance obscures the effects of {beta}-cell specific insulin resistance. The propensity for hyperinsulinemia was associated with mildly reduced fasting glucose and increased body weight. We further validated our main in vivo findings using the Ins1-CreERT transgenic line and found that male mice had improved glucose tolerance 4 weeks after tamoxifen-mediated Insr deletion. Collectively, our data show that loss of {beta}-cell Insr contributes to glucose-induced hyperinsulinemia, thereby improving glucose homeostasis in otherwise insulin sensitive dietary and age contexts. Hydrophobic mismatch is a key factor in protein transport on the Tat pathway The twin-arginine translocation (Tat) pathway transports folded proteins across membranes in bacteria, thylakoid, plant mitochondria, and archaea. In most species, the active Tat machinery consists of three independent subunits, TatA, TatB and TatC. TatA and TatB from all bacterial species possess short transmembrane alpha-helices (TMHs), both of which are only fifteen residues long in E. coli. Such short TMHs cause a hydrophobic mismatch between Tat subunits and the membrane bilayer. Here, by modifying the length of the TMHs of E. coli TatA and TatB, we assess the functional importance of the hydrophobic mismatch in the Tat transport mechanism. Surprisingly, both TatA and TatB with as few as 11 residues in their respective TMHs are still able to insert into the membrane bilayer, albeit with a decline in membrane integrity. Three different assays, both qualitative and quantitative, were conducted to evaluate the Tat activity of the TMH length mutants. Our experiments indicate that the TMHs of TatA and TatB appear to be evolutionarily tuned to 15 amino acids, with activity dropping off with any modification of this length. We believe our study supports a model of Tat transport utilizing localized toroidal pores that form when the membrane bilayer is thinned to a critical threshold. In this context, the 15-residue length of the TatA and TatB TMHs can be seen as a compromise between the need for some hydrophobic mismatch to allow the membrane to reversibly reach the threshold thinness required for toroidal pore formation, and the permanently destabilizing effect of placing even shorter helices into these energy-transducing membranes. A phenomics approach for in vitro antiviral drug discovery BackgroundThe current COVID-19 pandemic has highlighted the need for new and fast methods to identify novel or repurposed therapeutic drugs. Here we present a method for untargeted phenotypic drug screening of virus-infected cells, combining Cell Painting with antibody-based detection of viral infection in a single assay. We designed an image analysis pipeline for segmentation and classification of virus-infected and non-infected cells, followed by extraction of morphological properties.

ResultsWe show that the methodology can successfully capture virus-induced phenotypic signatures of MRC-5 human lung fibroblasts infected with Human coronavirus 229E (CoV-229E). Moreover, we demonstrate that our method can be used in phenotypic drug screening using a panel of nine host- and virus-targeting antivirals. Treatment with effective antiviral compounds reversed the morphological profile of the host cells towards a non-infected state.

ConclusionsThe method can be used in drug discovery for morphological profiling of novel antiviral compounds on both infected and non-infected cells.