Cryptic pathogen-sugar interactions revealed by universal saturation transfer analysis Many host pathogen interactions such as human viruses (including non-SARS-coronaviruses) rely on attachment to host cell-surface glycans. There are conflicting reports about whether the Spike protein of SARS-CoV-2 binds to sialic acid commonly found on host cell-surface N-linked glycans. In the absence of a biochemical assay, the ability to analyze the binding of glycans to heavily- modified proteins and resolve this issue is limited. Classical Saturation Transfer Difference (STD) NMR can be confounded by overlapping sugar resonances that compound with known experimental constraints. Here we present  universal saturation transfer analysis (uSTA), an NMR method that builds on existing approaches to provide a general and automated workflow for studying protein-ligand interactions. uSTA reveals that B-origin-lineage-SARS-CoV-2 spike trimer binds sialoside sugars in an  end on manner and modelling guided by uSTA localises binding to the spike N-terminal domain (NTD). The sialylated-polylactosamine motif is found on tetraantennary human N-linked-glycoproteins in deeper lung and may have played a role in zoonosis. Provocatively, sialic acid binding is abolished by mutations in some subsequent SARS- CoV-2 variants-of-concern. A very high resolution cryo-EM structure confirms the NTD location and  end on mode; it rationalises the effect of NTD mutations and the structure-activity relationship of sialic acid analogues. uSTA is demonstrated to be a robust, rapid and quantitative tool for analysis of binding, even in the most demanding systems.

Extended AbstractThe surface proteins found on both pathogens and host cells mediate entry (and exit) and influence disease progression and transmission. Both types can bear host-generated post- translational modifications such as glycosylation that are essential for function but can confound biophysical methods used for dissecting key interactions. Several human viruses (including non- SARS-coronaviruses) attach to host cell-surface N-linked glycans that include forms of sialic acid (sialosides). There remains, however, conflicting evidence as to if or how SARS-associated coronaviruses might use such a mechanism. Here, we demonstrate quantitative extension of  saturation transfer protein NMR methods to a complete mathematical model of the magnetization transfer caused by interactions between protein and ligand. The method couples objective resonance-identification via a deconvolution algorithm with Bloch-McConnell analysis to enable a structural, kinetic and thermodynamic analysis of ligand binding beyond previously-perceived limits of exchange rates, concentration or system. Using an automated and openly available workflow this  universal saturation transfer analysis (uSTA) can be readily-applied in a range of even heavily-modified systems in a general manner to now obtain quantitative binding interaction parameters (KD, kEx). uSTA proved critical in mapping direct interactions between natural sialoside sugar ligands and relevant virus-surface attachment glycoproteins - SARS-CoV-2-spike and influenza-H1N1-haemagglutinin variants - by quantitating ligand signal in spectral regions otherwise occluded by resonances from mobile protein glycans (that also include sialosides). In B- origin-lineage-SARS-CoV-2 spike trimer  end on-binding to sialoside sugars was revealed contrasting with  extended surface-binding for heparin sugar ligands; uSTA-derived constraints used in structural modelling suggested sialoside-glycan binding sites in a beta-sheet-rich region of spike N-terminal domain (NTD). Consistent with this, uSTA-glycan binding was minimally- perturbed by antibodies that neutralize the ACE2-binding domain (RBD) but strongly disrupted in spike from the B1.1.7/alpha and B1.351/beta variants-of-concern, which possess hotspot mutations in the NTD. Sialoside binding in B-origin-lineage-NTD was unequivocally pinpointed by cryo-EM to a site that is created from residues that are notably deleted in variants (e.g. H69,V70,Y145 in alpha). An analysis of beneficial genetic variances in cohorts of patients from early 2020 suggests a model in which this site in the NTD of B-origin-lineage-SARS-CoV-2 (but not in alpha/beta-variants) may have exploited a specific sialylated-polylactosamine motif found on tetraantennary human N-linked-glycoproteins in deeper lung. Together these confirm a novel binding mode mediated by the unusual NTD of SARS-CoV-2 and suggest how it may drive virulence and/or zoonosis via modulation of glycan attachment. Since cell-surface glycans are widely relevant to biology and pathology, uSTA can now provide ready, quantitative, widespread analysis of complex, host-derived and post-translationally modified proteins with putative ligands relevant to disease even in previously confounding complex systems.
The physics of liquid-to-solid transitions in multi-domain protein condensates Many RNA-binding proteins (RBPs) that assemble into membraneless organelles, have a common architecture including disordered prion-like domain (PLD) and folded RNA-binding domain (RBD). An enrichment of PLD within the condensed phase gives rise to formation, on longer time scales, amyloid-like fibrils (aging). In this study, we employ coarse-grained Langevin dynamics simulations to explore the physical basis for the structural diversity in condensed phases of multi-domain RBPs. We discovered a highly cooperative first order transition between disordered structures and an ordered phase whereby chains of PLD organize in fibrils with high nematic orientational order. An interplay between homo-domain (PLD-PLD) and hetero-domain (PLD-RBD) interactions results in variety of structures with distinct spatial architectures. Interestingly, the different structural phases also exhibit vastly different intra-cluster dynamics of proteins, with diffusion coefficients 5 (disordered structures) to 50 times (ordered structures) lower than that of the dilute phase. Cooperativity of this liquid-solid transition makes fibril formation highly malleable to mutations or post-translational modifications. Our results provide a mechanistic understanding of how multi-domain RBPs could form assemblies with distinct structural and material properties.

Significance StatementAssembly of proteins and nucleic acids into dense, liquid-like pockets is associated with several key functions including stress response, gene-regulation, DNA-repair and RNA processing. Several RNA binding proteins such as FUS are known to form liquid-like condensates that progressively harden into more dynamically, solid-like structures, a phenomenon that gets accelerated by disease mutations. In this study, we discover the mechanistic origins of this transition and show that small mutational or posttranslational modifications could result in sharp disorder-order transitions that could characterize accelerated liquid-solid transition in disease mutants.
Trade-offs in biosensor optimization for dynamic pathway engineering Recent progress in synthetic biology allows the construction of dynamic control circuits for metabolic engineering. This technology promises to overcome many challenges encountered in traditional pathway engineering, thanks to their ability to self-regulate gene expression in response to bioreactor perturbations. The central components in these control circuits are metabolite biosensors that read out pathway signals and actuate enzyme expression. However, the construction of metabolite biosensors is a major bottleneck for strain design, and a key challenge is to understand the relation between biosensor dose-response curves and pathway performance. Here we employ multiobjective optimization to quantify performance trade-offs that arise in the design and calibration of metabolite biosensors. Our approach reveals strategies for tuning dose-response curves along an optimal trade-off between production flux and the cost of an increased expression burden on the host. We explore properties of control architectures built in the literature, and identify their advantages and caveats in terms of performance and robustness to growth conditions and leaky promoters. We demonstrate the optimality of a control circuit for glucaric acid production in Escherichia coli, which has been shown to increase titer by 2.5-fold as compared to static designs. Our results lay the groundwork for the automated design of control circuits for pathway engineering, with applications in the food, energy and pharmaceutical sectors. SARS-CoV-2 spike-glycoprotein processing at S1/S2 and S2' and shedding of the ACE2 viral receptor: roles of Furin and TMPRSS2 and implications for viral infectivity and cell-to-cell fusion The spike (S)-protein of SARS-CoV-2 binds ACE2 and requires proteolytic "priming" at PRRAR685{downarrow} into S1 and S2 (cleavage at S1/S2), and "fusion-activation" at a S2 site for viral entry. In vitro, Furin cleaved peptides mimicking the S1/S2 cleavage site more efficiently than at the putative S2, whereas TMPRSS2 inefficiently cleaved both sites. In HeLa cells Furin-like enzymes mainly cleaved at S1/S2 during intracellular protein trafficking, and S2 processing by Furin at KPSKR815{downarrow} was strongly enhanced by ACE2, but not for the optimized S2 KRRKR815{downarrow} mutant (S2), whereas individual/double KR815AA mutants were retained in the endoplasmic reticulum. Pharmacological Furin-inhibitors (Boston Pharmaceuticals, BOS-inhibitors) effectively blocked endogenous S-protein processing in HeLa cells. Furthermore, we show using pseudotyped viruses that while entry by a "pH-dependent" endocytosis pathway in HEK293 cells did not require Furin processing at S1/S2, a "pH-independent" viral entry in lung-derived Calu-3 cells was sensitive to inhibitors of Furin (BOS) and TMPRSS2 (Camostat). Consistently, these inhibitors potently reduce infectious viral titer and cytopathic effects, an outcome enhanced when both compounds were combined. Quantitative analyses of cell-to-cell fusion and spike processing revealed the key importance of the Furin sites for syncytia formation. Our assays showed that TMPRSS2 enhances fusion and proteolysis at S2 in the absence of cleavage at S1/S2, an effect that is linked to ACE2 shedding by TMPRSS2. Overall, our results indicate that Furin and TMPRSS2 play synergistic roles in generating fusion-competent S-protein, and in promoting viral entry, supporting the combination of Furin and TMPRSS2 inhibitors as potent antivirals against SARS-CoV-2.

IMPORTANCESARS-CoV-2 is the etiological agent of COVID-19 that resulted in >5 million deaths. The spike protein (S) of the virus directs infection of the lungs and other tissues by binding the angiotensin-converting enzyme 2 (ACE2) receptor. For effective infection, the S-protein is cleaved at two sites: S1/S2 and S2. Cleavage at S1/S2, induces a conformational change favoring the recognition of ACE2. The S2 cleavage is critical for cell-to-cell fusion and virus entry into host cells. Our study contributes to a better understanding of the dynamics of interaction between Furin and TMPRSS2 during SARS-CoV-2 entry and suggests that the combination of a non-toxic Furin inhibitor with a TMPRSS2 inhibitor could significantly reduce viral entry in lung cells, as evidenced by an average synergistic [~]95% reduction of viral infection. This represents a powerful novel antiviral approach to reduce viral spread in individuals infected by SARS-CoV-2 or future related coronaviruses. 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. Genetic interactions drive heterogeneity in causal variant effect sizes for gene expression and complex traits Despite the growing number of genome-wide association studies (GWAS), it remains unclear to what extent gene-by-gene and gene-by-environment interactions influence complex traits in humans. The magnitude of genetic interactions in complex traits has been difficult to quantify because GWAS are generally underpowered to detect individual interactions of small effect. Here, we develop a method to test for genetic interactions that aggregates information across all trait-associated loci. Specifically, we test whether SNPs in regions of European ancestry shared between European American and admixed African American individuals have the same causal effect sizes. We hypothesize that in African Americans, the presence of genetic interactions will drive the causal effect sizes of SNPs in regions of European ancestry to be more similar to those of SNPs in regions of African ancestry. We apply our method to two traits: gene expression in 296 African Americans and 482 European Americans in the Multi-Ethnic Study of Atherosclerosis (MESA) and low-density lipoprotein cholesterol (LDL-C) in 74K African Americans and 296K European Americans in the Million Veteran Program (MVP). We find significant evidence for genetic interactions in our analysis of gene expression; for LDL-C, we observe a similar point estimate although this is not significant, likely due to lower statistical power. These results suggest that gene-by-gene or gene-by-environment interactions modify the effect sizes of causal variants in human complex traits. Chromosome fusion affects genetic diversity and evolutionary turnover of functional loci, but consistently depends on chromosome size Major changes in chromosome number and structure are linked to a series of evolutionary phenomena, including intrinsic barriers to gene flow or suppression of recombination due to chromosomal rearrangements. However, chromosome rearrangements can also affect the fundamental dynamics of molecular evolution within populations by changing relationships between linked loci and altering rates of recombination. Here, we build chromosome-level assembly Eueides isabella and, together with the chromosome-level assembly of Dryas iulia, examine the evolutionary consequences of multiple chromosome fusions in Heliconius butterflies. These assemblies pinpoint fusion points on 10 of the 21 autosomal chromosomes and reveal striking differences in the characteristics of fused and unfused chromosomes. The ten smallest autosomes in D. iulia and E. isabella, which have each fused to a longer chromosome in Heliconius, have higher repeat and GC content, and longer introns than predicted by their chromosome length. Following fusion, these characteristics change to become more in line with chromosome length. The fusions also led to reduced diversity, which likely reflects increased background selection and selection against introgression between diverging populations, following a reduction in per-base recombination rate. We further show that chromosome size and fusion impact turnover rates of functional loci at a macroevolutionary scale. Together these results provide further evidence that chromosome fusion in Heliconius likely had dramatic effects on population level processes shaping rates of neutral and adaptive divergence. These effects may have impacted patterns of diversification in Heliconius, a classic example of an adaptive radiation. Genetic variation associated with plastic and homeostatic growth responses to drought in Arabidopsis O_LINatural genetic variation influences plant responses to environmental stressors. However, the extent to which such variation underlies plastic versus homeostatic response phenotypes has been little studied.
C_LIO_LIWe quantified the extent of drought-induced changes in leaf area in a set of Iberian Arabidopsis accessions and then performed association studies correlating variation in this phenotype with genomic and transcriptomic variation.
C_LIO_LIDrought-induced plastic reductions in relative leaf area typified accessions originating from productive environments. In contrast, homeostasis in relative leaf area typified accessions originating from unproductive environments. Genome-Wide Association Studies (GWAS), Transcriptome Wide Association Studies (TWAS), and expression GWAS (eGWAS) highlighted the importance of auxin-related processes and, in particular, the potential role of the SMALL AUXIN UP RNA 26 (SAUR26) gene in conferring leaf area plasticity. Homeostatic responses in relative leaf area were associated with a diverse gene set and positively associated with a higher intrinsic water use efficiency (WUEi) as confirmed in a TWAS metanalysis of previously published {delta}13C measurements.
C_LIO_LIWe have identified not only candidate "plasticity genes" but also candidate "homeostasis genes" controlling leaf area. Our results demonstrate the value of a combined GWAS, TWAS, and eGWAS approach to identify mechanisms underlying phenotypic responses to stress.
C_LI

One-sentence summaryInformation on phenotype, genotype, and transcript abundance is integrated to identify both plasticity and homeostasis genes and processes associated with local adaptation to drought stress in Arabidopsis accessions of the Iberian Peninsula. Relative contributions of sex hormones, sex chromosomes, and gonads to sex differences in tissue gene regulation Sex differences in physiology and disease in mammals result from the effects of three classes of factors that are inherently unequal in males and females: reversible (activational) effects of gonadal hormones, permanent (organizational) effects of gonadal hormones, and cell-autonomous effects of sex chromosomes, as well as genes driven by these classes of factors. Often, these factors act together to cause sex differences in specific phenotypes, but the relative contribution of each and the interactions among them remain unclear. Here, we used the Four Core Genotypes (FCG) mouse model with or without hormone replacement to distinguish the effects of each class of sex-biasing factors on transcriptome regulation in liver and adipose tissues. We found that the activational hormone levels have the strongest influence on gene expression, followed by the organizational gonadal sex effect and, lastly, sex chromosomal effect, along with interactions among the three factors. Tissue specificity was prominent, with a major impact of estradiol on adipose tissue gene regulation, and of testosterone on the liver transcriptome. The networks affected by the three sex-biasing factors include development, immunity and metabolism, and tissue-specific regulators were identified for these networks. Furthermore, the genes affected by individual sex-biasing factors and interactions among factors are associated with human disease traits such as coronary artery disease, diabetes, and inflammatory bowel disease. Our study offers a tissue-specific account of the individual and interactive contributions of major sex-biasing factors to gene regulation that have broad impact on systemic metabolic, endocrine, and immune functions. Experimental evolution in TP53 deficient human gastric organoids recapitulates tumorigenesis The earliest events during human tumor initiation are poorly characterized but may hold clues as to how to detect and prevent malignancy. Here we model occult pre-neoplasia by engineering TP53 deficiency, a common and early event in gastric cancer, into human gastric organoids. By performing experimental evolution in multiple clonally derived cultures over two years we define causal relationships between this common initiating genetic lesion and resulting phenotypes. TP53 loss elicited progressive aneuploidy, including copy number alterations and complex structural variants that are common in gastric cancers and which follow preferred temporal orders. Longitudinal single cell sequencing of TP53 deficient gastric organoids similarly indicates progression towards malignant transcriptional programs. Moreover, lineage tracing with expressed cellular barcodes demonstrates reproducible dynamics whereby initially rare subclones with shared transcriptional programs repeatedly attain clonal dominance. This powerful platform for experimental evolution exposes stringent selection, clonal interference and a striking degree of phenotypic convergence in pre-malignant epithelial organoids, implying that the earliest stages of tumorigenesis may be predictable while illuminating evolutionary constraints and barriers to malignant transformation.

Graphical abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=135 SRC="FIGDIR/small/487529v2_ufig1.gif" ALT="Figure 1">
View larger version (63K):
org.highwire.dtl.DTLVardef@122debborg.highwire.dtl.DTLVardef@4cd2daorg.highwire.dtl.DTLVardef@1a720d6org.highwire.dtl.DTLVardef@15ccaa7_HPS_FORMAT_FIGEXP  M_FIG C_FIG Circulating multimeric immune complexes drive immunopathology in COVID-19 A dysregulated immune response with high levels of SARS-CoV-2 specific IgG antibodies characterizes patients with severe or critical COVID-19. Although a robust IgG response is traditionally considered to be protective, excessive triggering of activating Fc-gamma-receptors (Fc{gamma}Rs) could be detrimental and cause immunopathology. Here, we document that patients who develop soluble circulating IgG immune complexes (sICs) during infection are subject to enhanced immunopathology driven by Fc{gamma}R activation. Utilizing cell-based reporter systems we provide evidence that sICs are predominantly formed prior to a specific humoral response against SARS-CoV-2. sIC formation, together with increased afucosylation of SARS-CoV-2 specific IgG eventually leads to an enhanced CD16 (Fc{gamma}RIII) activation of immune cells reaching activation levels comparable active systemic lupus erythematosus (SLE) disease. Our data suggest a vicious cycle of escalating immunopathology driven by an early formation of sICs in predisposed patients. These findings reconcile the seemingly paradoxical findings of high antiviral IgG responses and systemic immune dysregulation in severe COVID-19.

Clinical implicationsThe identification of sICs as drivers of an escalating immunopathology in predisposed patients opens new avenues regarding intervention strategies to alleviate critical COVID-19 progression.

Graphical abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=78 SRC="FIGDIR/small/449893v4_ufig1.gif" ALT="Figure 1">
View larger version (17K):
org.highwire.dtl.DTLVardef@b4616corg.highwire.dtl.DTLVardef@682d1aorg.highwire.dtl.DTLVardef@16946cborg.highwire.dtl.DTLVardef@a6ef7d_HPS_FORMAT_FIGEXP  M_FIG C_FIG A vicious cycle of immunopathology in COVID-19 patients is driven by soluble multimeric immune complexes (sICs). SARS-CoV-2 infection triggers sIC formation in prone individuals. Activation of Fc{gamma}RIII/CD16 expressing immune cells by sICs precedes a humoral response to SARS-CoV2 infection. sICs and infection add to IgG afucosylation, further enhancing Fc{gamma}RIII/CD16 activation by opsonized targets. High inflammation induces further sIC mediated immune cell activation ultimately leading to an escalating immunopathology. Comparison of hybrid learning and remote education in the implementation of the 'Adopt a Microorganism' methodology The Internet has changed the way teachers and students access information and build knowledge. The recent COVID-19 pandemic has created challenges for both teachers and students and a demand for new methodologies of remote learning. In the life sciences, mixing online content with practical activities represents an even greater challenge. In microbiology, the implementation of an active teaching methodology, the #Adopt project, based on the social network Facebook(R), represents an excellent option for connecting remote education with classroom activities. In 2020, the version applied in high school, "Adopt a Microorganism", was adapted to meet the demands of emergency remote education owing to the suppression of face-to-face activities caused by the pandemic. In the present study, we assessed how the change in methodology impacted the discourse richness of students from high school integrated with technical education in the Business Administration program of the Federal Institute of Sao Paulo, Sorocaba Campus. Three questionnaires related to the groups of microorganisms (Archaea, Bacteria, Virus, Fungi, and Protozoan) were applied. The students responses in the 2019 and 2020 classes were compared concerning content richness and multiplicity of concepts through the application of the Shannon diversity index, an approach that is generally used to assess biodiversity in different environments. The observed results suggest that remote learning provided students with a conceptual basis and richness of content equivalent to that achieved by students subjected to the hybrid teaching model. In conclusion, this study suggests that the #Adopt project methodology increases students discourse richness in microbiology even without face-to-face traditional classes. IRGM/Irgm1 Aggravates Progression of Atherosclerosis by Inducing Macrophage Apoptosis through the MAPK Signaling Pathway AimsAtherosclerosis underlies most cardiovascular diseases, among which acute coronary syndrome (ACS) caused by plaque rupture (PR) often leads to death. Immune-related GTPases (IRGM/Irgm1) have been extensively studied in inflammatory diseases, but their role in atherosclerosis is unclear. Determining how IRGM/Irgm1 promotes atherosclerotic plaque vulnerability will provide information for new biomarkers and/or therapeutic targets.

Methods and resultsWe identified ruptured and unruptured plaques by optical coherence tomography, and found that serum IRGM was highly expressed in patients with ST-segment elevation myocardial infarction in PR. We used ApoE-/-Irgm1+/+, ApoE-/-Irgm1+/- mice and chimeric mice to establish a model of advanced atherosclerosis. The results of pathological experiments showed that Irgm1 caused plaque necrosis. The ratio of neutral lipids and cholesterol crystals increases, while the content of collagen fibers decreases, aggravating the destabilization of atherosclerotic plaques. In vitro, we used multiple approaches to confirm that Irgm1 promotes macrophage apoptosis by promoting the production of reactive oxygen species and activating the MAPK signaling pathway.

ConclusionsIRGM may be a potential risk factor for PR. Mechanistic studies have shown that IRGM/Irgm1 contributes to the formation and rupture of fragile plaques. This is partly mediated by the induction of macrophage apoptosis via the MAPK signaling pathway. IRGM may offer new strategies for early treatment of ACS.

Translation viewOur findings indicate that IRGM/Irgm1 contributes to formation and rupture of vulnerable plaques. It suggests that IRGM may provide a new target for the early treatment of ACS. A unified mechanism for innate and learned visual landmark guidance in the insect central complex Insects can navigate efficiently in both novel and familiar environments, and this requires flexiblity in how they are guided by sensory cues. A prominent landmark, for example, can ellicit strong innate behaviours (attraction or menotaxis) but can also be used, after learning, as a specific directional cue as part of a navigation memory. However, the mechanisms that allow both pathways to co-exist, interact or override each other are largely unknown. Here we propose a model for the behavioural integration of innate and learned guidance based on the neuroanatomy of the central complex (CX), adapted to control landmark guided behaviours. We consider a reward signal provided either by an innate attraction to landmarks or a long-term visual memory in the mushroom bodies (MB) that modulates the formation of a local vector memory in the CX. Using an operant strategy for a simulated agent exploring a simple world containing a single visual cue, we show how the generated short-term memory can support both innate and learned steering behaviour. In addition, we show how this architecture is consistent with the observed effects of unilateral MB lesions in ants that cause a reversion to innate behaviour. We suggest the formation of a directional memory in the CX can be interpreted as transforming rewarding (positive or negative) sensory signals into a mapping of the environment that describes the geometrical attractiveness (or repulsion). We discuss how this scheme might represent an ideal way to combine multisensory information gathered during the exploration of an environment and support optimized cue integration.

1 Author summaryIn this paper, we modeled the neural pathway allowing insects to perform landmark guided behaviours using their internal compass. First, we observed the intrinsic property of the connectome, extracted from drosophila online database, between the internal compass neurons and the steering neurons to support an oriented behaviour towards a single landmark. Then, we proposed and evaluated an adaptation of the bees path integration neural circuit, to sustain flexible landmark guidance behaviours such as attraction or menotaxis. We showed the model ability to form a memory during the exploration of the local environment to support both innate or learned navigation behaviour using a single landmark in the environment. In addtion, we demonstrated the transformation of a simple goodness/badness signal, from innate or long-term memory pathways, into an oriented steering signal that could be applied to other sensory pathways. Furthermore, by reproducing lesion experiments in the mushroom bodies of wood ants we highlight the consistency of the model with biological observation. We then discuss the different emergent properties and the potential outcome that this local, and operant, memory supports. Electrophysiological correlates of mood and reward dynamics in human adolescents Despite its omnipresence in everyday interactions and its importance for mental health, mood and its neuronal underpinnings are poorly understood. Computational models can help identify parameters affecting self-reported mood during mood induction tasks. Here we test if computationally modelled dynamics of self-reported mood during monetary gambling can be used to identify trial-by-trial variations in neuronal activity. To this end, we shifted mood in healthy (N=24) and depressed (N=30) adolescents by delivering individually tailored reward prediction errors whilst recording magnetoencephalography (MEG) data. Following a pre-registered analysis, we hypothesize that expectation (defined by previous reward outcomes) would be predictive of beta-gamma oscillatory power (25-40Hz), a frequency shown to modulate to reward feedback. We also hypothesize that trial variations in the evoked response to the presentation of gambling options and in source localized responses to reward feedback. Through our multilevel statistical analysis, we found confirmatory evidence that beta-gamma power is positively related to reward expectation during mood shifts, with possible localized sources in the posterior cingulate cortex. We also confirmed reward prediction error to be predictive of trial-level variations in the response of the paracentral lobule and expectation to have an effect on the cerebellum after presentation of gambling options. To our knowledge, this is the first study to relate fluctuations in mood on a minute timescale to variations in neural oscillations with noninvasive electrophysiology.

Significance StatementBrain mechanisms underlying mood and its relationship with changes in reward contingencies in the environment are still elusive but could have a strong impact on our understanding and treatment of debilitating mood disorders. Building on a previously proposed computational mood model we use multilevel statistical models to find relationship between trial-by-trial variations in model components of mood and neural responses to rewards measured with non-invasive electrophysiology (MEG). Through confirmatory analysis we show that it is possible to observe relationships between trial variations in neural responses and computational parameters describing mood dynamics. Identifying the dynamics of mood and the neural processes it affects could pave the way for more effective neuromodulation treatments. Diversity of SARS-CoV-2 genome among various strain identified in Lucknow, Uttar Pradesh Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has emerged as a significant challenge worldwide. Rapid genome sequencing of SARS-CoV-2 is going on across the globe to detect mutations and genomic modifications in SARS-CoV-2. In this study, we have sequenced twenty-three SARS-CoV-2 positive samples collected during the first pandemic from the state of Uttar Pradesh, India. We observed a range of already reported mutations (2-22), including; D614G, L452R, Q613H, Q677H, T1027I in the S gene; S194L in the N gene; Q57H, L106F, T175I in the ORF3. Few unreported mutations such as P309S in the ORF1ab gene; T379I in the N gene; and L52F, V77I in the ORF3a gene were also detected. Phylogenetic genome analysis showed similarity with other SARS-CoV-2 viruses reported from Uttar Pradesh. The observed mutations may be associated with SARS-CoV-2 virus pathogenicity or disease severity. Mouse lemur, a new primate model to investigate cardiac pacemaker activity. BackgroundGray mouse lemurs (Microcebus murinus) are one of the smallest and most ancestral primates known. Their size falls in between that of mice and rats, and their genetic proximity to humans makes them an emerging model to study age-related neurodegeneration. Since Microcebus murinus replicates similar senescence processes of humans, it constitutes a useful model also for studying cardiovascular dysfunctions. However, its cardiac physiology is unknown. Thus, we investigated sinoatrial node (SAN) pacemaker activity of mouse lemurs, presenting the first characterization of heart automaticity in nonhuman primates.

Methods and ResultsWe recorded the heart automaticity in freely moving mouse lemurs and in their SAN tissues and pacemaker myocytes. Mouse lemurs have a heart rate (HR) in between those of mice and rats and a similar generation of electrical impulses in the SAN tissue. Conversely, their ventricular depolarization is more similar to large mammals. Mouse lemur SAN myocytes express the main pacemaker currents at densities similar to mice: the hyperpolarization-activated current (If) and the L-type (ICa,L) and T-type (ICa,T) calcium currents. Moreover, mouse lemur SAN myocytes display spontaneous ryanodine receptor (RyR)-dependent Ca2+ transients. Using muscle-derived stem cells (MDSCs) from the Microcebus murinus skeletal muscle, we also differentiated pacemaker-like (PML) cells showing spontaneous automaticity and expressing markers of native SAN myocytes (HCN4 and connexin-45).

ConclusionsIn conclusion, our characterization of heart automaticity in Microcebus murinus provides new opportunities for comparative cardio-physiology studies in primates and with humans, and for pharmacologic testing of molecules that could assist in modulating age-related dysfunctions of automaticity. Inter-microscope comparability of dental microwear texture data obtained from different optical profilometers. Dental microwear texture analysis has become a well-established approach for dietary inference and reconstruction of mammals and other tetrapods, both extant and extinct. As the amount of available data grows immensely, researchers could benefit from combining data gathered by others to perform meta-analyses. However, different devices used to capture three-dimensional surface scans for DMTA are known to produce variation even when measuring the same surface. Here we compare DMTA data of 36 guinea pigs that received different diets in a controlled feeding experiment, measured on a confocal disc-scanning and a confocal laser-scanning microscope. We are testing different pre-analysis filtering protocols to mitigate differences. We find inter-microscopes and filter-related differences for the majority of analysed 41 DMTA parameters. Certain microscope-specific filter routines resulted in less differences than other pre-analysis protocols. We further identify DMTA parameters which were stable regardless of microscope or pre-analysis treatment. Overall, the results obtained on both microscopes show the same dietary differentiation between guinea pig feeding groups, which supports that DMTA is a suitable method to obtain repeatable, objective dietary inferences. We finally propose a roadmap to enhance data exchange and inter-lab comparability and collaboration in the future. 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. Multiparametric profiling of cellular lipid metabolism in hypercholesterolemia Systematic insight into cellular dysfunctions can improve understanding of disease etiology, risk assessment and patient stratification. We present a multiparametric high-content imaging platform enabling quantification of low-density lipoprotein (LDL) uptake and lipid storage in cytoplasmic droplets of primary leukocyte subpopulations. We validated this platform with samples from 65 individuals with variable blood LDL-cholesterol (LDL-c) levels, including familial hypercholesterolemia (FH) and non-FH subjects. We integrated lipid storage data into a novel readout, lipid mobilization, measuring the efficiency with which cells deplete lipid reservoirs. Lipid mobilization correlated positively with LDL uptake and negatively with hypercholesterolemia and age, improving differentiation of individuals with normal and elevated LDL-c. Moreover, combination of cell-based readouts with a polygenic risk score for LDL-c explained hypercholesterolemia better than the genetic risk score alone. This platform provides functional insights into cellular lipid trafficking from a few mls of blood and is applicable to dissect metabolic disorders, such as hypercholesterolemia.

MotivationWe have limited information on how cellular lipid uptake and processing differ between individuals and influence the development of metabolic diseases, such as hypercholesterolemia. Available assays are labor intensive, require skilled personnel and are difficult to scale to higher throughput, making it challenging to obtain systematic functional cell-based data from individuals. To overcome this problem, we established a scalable automated analysis pipeline enabling reliable quantification of multiple cellular readouts, including lipid uptake, storage and mobilization, from different white blood cell populations. This approach provides new personalized insights into the cellular basis of hypercholesterolemia and obesity.

Graphical Abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=200 SRC="FIGDIR/small/440471v4_ufig1.gif" ALT="Figure 1">
View larger version (84K):
org.highwire.dtl.DTLVardef@9222c1org.highwire.dtl.DTLVardef@27e289org.highwire.dtl.DTLVardef@89ba82org.highwire.dtl.DTLVardef@33c9cd_HPS_FORMAT_FIGEXP  M_FIG C_FIG Relationship between shear wave velocity and muscle activation is inconsistent across different muscle types There is an increasing use of shear wave ultrasound elastography to quantify mechanical properties of muscles under various conditions such as changes muscle length and levels of activation in healthy and pathological muscle. However, little is known about the variability in shear wave velocity among muscles as most studies investigate one specific muscle. The purpose of this study was to determine if the relationship between SWV and muscle activation is consistent across muscles with different architectures: biceps brachii, tibialis anterior, and medial gastrocnemius, All measures were made at matching levels of activation and approximately at the optimal length for each muscle to control for length-dependent changes in the relationship between activation and force or stiffness. We also conducted a control experiment to determine how the passive force within a muscle alters the relationship between muscle activation and shear wave velocity. The relationship between shear wave velocity-squared and activation above 10% MVC differed across muscles, with biceps brachii and medial gastrocnemius showing a lower slope than tibialis anterior. Shear wave velocity-squared also differed between muscles at the shortest length (p<0.001) and the increase in shear wave velocity-squared with muscle lengthening also differed between muscle types (p = 0.005) Muscle-specific differences could not be explained by the architectural features such as pennation angle, fiber length, and physiological cross-sectional area. Our results demonstrate that there is not a unique relationship between muscle activation and shear wave velocity, highlighting the importance of understanding the many factors contributing to shear wave propagation in muscle before elastography can be used to make quantitative comparisons across muscle types. Newly developed ad hoc molecular assays show how eDNA can witness and anticipate the monk seal recolonization of central Mediterranean The monk seal Monachus monachus is the most endangered pinniped worldwide and the only one found in the Mediterranean, where its distribution and abundance have suffered a drastic decline in the last few decades. Data on its status are scattered due to both its rarity and evasiveness, and records are biased towards occasional, mostly coastal, encounters. Nowadays molecular techniques allow us to detect and quantify minute amounts of DNA traces released in the environment (eDNA) by any organism. We present three qPCR-assays targeting the monk seal mitogenome. The assays were soundly tested on an extensive and diversified sample set (n=73), including positive controls from Madeira breeding population collected during the peak of abundance, and two opportunistic Mediterranean eDNA-sample collections (offshore/coastal) from on-going projects. Monk seal DNA was detected in 47.2% and 66.7% of the samples collected in the Tyrrhenian from a ferry platform (2018-2019) and in the Pelagie archipelago -Strait of Sicily- (2020) respectively, anticipating (up to 2 year) visual observations occurred subsequently in proximity of the sampled areas. In the Tyrrhenian, detection occurrence increased between 2018 and 2019. Monk seal DNA recoveries were commoner in night-time ferry-samples, suggesting nocturnal predatory activity in pelagic waters. The proposed technique provides a non-invasive and yet highly-sensitive tool for defining the monk seal actual distribution and home range, its recovery rate and pinpoint coastal/offshore localities where prioritizing conservation, research, citizen science and education initiatives. Host factors drive the within-host competition between highly and low pathogenic H5N8 avian influenza viruses Highly pathogenic avian influenza viruses (HPAIV) emerge from low pathogenic avian influenza viruses (LPAIV) through the introduction of basic amino acids at the hemagglutinin (HA) cleavage site. Following viral evolution, the newly formed HPAIV likely represents a minority variant within the index host, predominantly infected with the LPAIV precursor. Using reverse-genetics engineered H5N8 viruses differing solely at the HA cleavage, we tested the hypothesis that the interaction between the minority HPAIV and the majority LPAIV could modulate the risk of HPAIV emergence and that the nature of the interaction could depend on the host species. In chickens, we observed that the H5N8LP increased H5N8HP replication and pathogenesis. By contrast, the H5N8LP antagonized H5N8HP replication and pathogenesis in ducks. Ducks mounted a more potent antiviral innate immune response than chickens against the H5N8LP, which correlated with H5N8HP inhibition. These data provide experimental evidence that HPAIV may be more likely to emerge in chickens than in ducks and underscore the importance of within-host viral variants interactions in viral evolution. SARS-CoV-2 variants show a gradual declining pathogenicity and pro-inflammatory cytokine spur, an increasing antigenic and anti-inflammatory cytokine induction, and rising structural protein instability Hyper-transmissibility with decreased disease severity are typical characteristics of Omicron variant. To understand this phenomenon, we used various bioinformatics approaches to analyze randomly selected genome sequences (one each) of the Gamma, Delta, and Omicron variants submitted to NCBI from 15 to 31 December 2021. We show that: (i) Pathogenicity of SARS-CoV-2 variants decreases in the order: Wuhan > Gamma > Delta > Omicron; however, the antigenic property follows the order: Omicron > Gamma > Wuhan > Delta. (ii) Omicron Spike RBD has lower pathogenicity but higher antigenicity than other variants. (iii) Decreased disease severity by Omicron variant may be due to its decreased pro-inflammatory and IL-6 stimulation and increased IFN-{gamma} and IL-4 induction efficacy. (iv) Mutations in N protein are associated with decreased IL-6 induction and human DDX21-mediated increased IL-4 production in Omicron. (v) Due to mutations, the stability of S, M, N, and E proteins decreases in the order: Omicron > Gamma > Delta > Wuhan. (vi) Stronger Spike RBD-hACE2 binding in Omicron is associated with increased transmissibility. However, the lowest stability of the Omicron Spike protein makes Spike RBD-hACE2 interaction weak for systemic infection and for causing severe disease. Finally (vii), the highest instability of Omicron E protein may also be associated with decreased viral maturation and low viral load leading to less severe disease and faster recovery. Our method may be used for other similar viruses, and these findings will contribute to the understanding of the dynamics of SARS-CoV-2 variants and the management of emerging variants. Sablefish (Anoplopoma fimbria Pallas, 1814) plasma biochemistry and hematology reference intervals including blood cell morphology Plasma biochemistry and hematology reference intervals are integral health assessment tools in all medical fields, including aquatic animal health. As sablefish (Anoplopoma fimbria) are becoming aquaculturally and economically more important, this manuscript provides essential reference intervals (RI) for their plasma biochemistry and hematology along with reference photomicrographs of blood cells in healthy, fasted sablefish. Blood cell morphology can differ between fish species. In addition, blood cell counts and blood chemistry can vary between fish species, demographics, water conditions, seasons, diets, and culture systems, which precludes the use of RIs from other fish species. For this study, blood was collected for plasma biochemistry and hematology analysis between June 20 and July 18, 2019, from healthy, yearling sablefish, hatched and reared in captivity on a commercial diet. Overnight fast of 16-18 hours did not sufficiently reduce lipids in the blood, which led to visible lipemia and frequent rupture of blood cells during analysis. Therefore, sablefish should be fasted for 24 to 36 hours before blood is collected to reduce hematology artifacts or possible reagent interference in plasma biochemistry analysis. Lymphocytes were the most dominant leukocytes (98%), while eosinophils were rare, and basophils were not detected in sablefish. Neutrophils were very large cells with Dohle bodies. In mammals and avian species, Dohle bodies are usually signs of toxic change from inflammation, but no such association was found in these fish. In conclusion, lipemia can interfere with sablefish blood analysis, and available removal methods should be evaluated as fasting for up to 36 h might not always be feasible. Also, more studies are required to establish RI for different developmental stages and rearing conditions. Multi-layered regulation of neuroectoderm differentiation by retinoic acid in a primitive streak-like context The formation of the primitive streak (PS) and the subsequent induction of neuroectoderm are hallmarks of gastrulation. Combining an in vitro reconstitution of this process based on mouse embryonic stem cells (mESCs) with a collection of knockouts in reporter mESC lines, we reassessed the contribution of retinoic acid (RA) signaling at early stages of neural commitment and its cross-talk with TGF{beta} and Wnt signaling inhibition. Single-cell RNA sequencing analysis captured the temporal unfolding of cell type diversification from epiblast and primitive streak-like cells up to the emergence of anterior and posterior neural fates. In conditions thought to lack RA synthesis, we discovered a hitherto unidentified residual RA production via a sensitive RA reporter. Genetic perturbations proved that the RA-degrading enzyme Cyp26a1 safeguard the developmental capabilities of the PS-like cells, limiting neural differentiation in wild type and in Chrd-/-Nog-/- or Dkk1-/- cells. Finally, the knockout of the three RAR receptors highlighted their function as negative regulators of loci critical for neural induction. Overall, we identified two mechanisms whereby components of the RA pathway can control the formation of neural progenitors in our PS-like context: a RA-dependent neural induction gated by RARs, and a receptor-mediated repression in the absence of ligand.