Refine
Year of publication
- 2022 (330) (remove)
Document Type
- Preprint (330) (remove)
Language
- English (330)
Has Fulltext
- yes (330)
Is part of the Bibliography
- no (330)
Keywords
Institute
- Physik (190)
- Frankfurt Institute for Advanced Studies (FIAS) (144)
- Informatik (102)
- Medizin (28)
- Ernst Strüngmann Institut (25)
- Biowissenschaften (22)
- Senckenbergische Naturforschende Gesellschaft (8)
- MPI für Hirnforschung (7)
- Psychologie (6)
- Biochemie, Chemie und Pharmazie (5)
Long non-coding RNAs (lncRNAs) can act as regulatory RNAs which, by altering the expression of target genes, impact on the cellular phenotype and cardiovascular disease development. Endothelial lncRNAs and their vascular functions are largely undefined. Deep RNA-Seq and FANTOM5 CAGE analysis revealed the lncRNA LINC00607 to be highly enriched in human endothelial cells. LINC00607 was induced in response to hypoxia, arteriosclerosis regression in non-human primates and also in response to propranolol used to induce regression of human arteriovenous malformations. siRNA knockdown or CRISPR/Cas9 knockout of LINC00607 attenuated VEGF-A-induced angiogenic sprouting. LINC00607 knockout in endothelial cells also integrated less into newly formed vascular networks in an in vivo assay in SCID mice. Overexpression of LINC00607 in CRISPR knockout cells restored normal endothelial function. RNA- and ATAC-Seq after LINC00607 knockout revealed changes in the transcription of endothelial gene sets linked to the endothelial phenotype and in chromatin accessibility around ERG-binding sites. Mechanistically, LINC00607 interacted with the SWI/SNF chromatin remodeling protein BRG1. CRISPR/Cas9-mediated knockout of BRG1 in HUVEC followed by CUT&RUN revealed that BRG1 is required to secure a stable chromatin state, mainly on ERG-binding sites. In conclusion, LINC00607 is an endothelial-enriched lncRNA that maintains ERG target gene transcription by interacting with the chromatin remodeler BRG1.
Several clinically used drugs are derived from microorganisms that often produce them via non-ribosomal peptide synthetases (NRPS), giant megasynthases that activate and connect individual amino acids in an assembly line fashion. Since NRPS are not restricted to the incorporation of the 20 proteinogenic amino acids, their efficient manipulation would allow the biotechnological generation of several different peptides including linear, cyclic and further modified derivatives. Here we describe a detailed phylogenetic analysis of several bacterial NRPS that led to the identification of a new recombination breakpoint within the thiolation (T) domain important in natural NRPS evolution. From this an evolutionary-inspired eXchange Unit between T domains (XUT) approach was developed, which allows the assembly of NRPS fragments over a broad range of GC contents, protein similarities, and extender unit specificities, as was shown for the specific production of a proteasome inhibitor, designed and assembled from five different NRPS fragments.
Many clinically used drugs are derived from or inspired by bacterial natural products that often are biosynthesised via non-ribosomal peptide synthetases (NRPS), giant megasynthases that activate and join individual amino acids in an assembly line fashion. Since NRPS are not limited to the incorporation of the 20 proteinogenic amino acids, their efficient manipulation would allow the biotechnological generation of complex peptides including linear, cyclic and further modified natural product analogues, e.g. to optimise natural product leads. Here we describe a detailed phylogenetic analysis of several bacterial NRPS that led to the identification of a new recombination breakpoint within the thiolation (T) domain that is important for natural NRPS evolution. From this, an evolution-inspired eXchange Unit between T domains (XUT) approach was developed which allows the assembly of NRPS fragments over a broad range of GC contents, protein similarities, and extender unit specificities, as demonstrated for the specific production of a proteasome inhibitor designed and assembled from five different NRPS fragments.
HER2 belongs to the ErbB sub-family of receptor tyrosine kinases and regulates cellular proliferation and growth. Different from other ErbB receptors, HER2 has no known ligand. Activation occurs through heterodimerization with other ErbB receptors and their cognate ligands. This suggests several possible activation paths of HER2 with ligand-specific, differential response, which so far remained unexplored. Using single-molecule tracking and the diffusion profile of HER2 as a proxy for activity, we measured the activation strength and temporal profile in live cells. We found that HER2 is strongly activated by EGFR-targeting ligands EGF and TGFα, yet with a distinguishable temporal fingerprint. The HER4-targeting ligands EREG and NRGβ1 showed weaker activation of HER2, a preference for EREG and a delayed response to NRGβ1. Our results indicate a selective ligand response of HER2 that may serve as a regulatory element. Our experimental approach is easily transferable to other membrane receptors targeted by multiple ligands.
Highlights
HER2 exhibits heterogeneous motion in the plasma membrane
The fraction of immobile HER2 correlates with phosphorylation levels
Diffusion properties serve as proxies for HER2 activation
HER2 exhibits ligand-specific activation strength and temporal profiles
Mechanisms by which specific histone modifications regulate distinct gene regulatory networks remain little understood. We investigated how H3K79me2, a modification catalyzed by DOT1L and previously considered a general transcriptional activation mark, regulates gene expression in mammalian cardiogenesis. Early embryonic cardiomyocyte ablation of Dot1l revealed that H3K79me2 does not act as a general transcriptional activator, but rather regulates highly specific gene regulatory networks at two critical cardiogenic junctures: left ventricle patterning and postnatal cardiomyocyte cell cycle withdrawal. Mechanistic analyses revealed that H3K79me2 in two distinct domains, gene bodies and regulatory elements, synergized to promote expression of genes activated by DOT1L. Surprisingly, these analyses also revealed that H3K79me2 in specific regulatory elements contributed to silencing genes usually not expressed in cardiomyocytes. As DOT1L mutants had increased numbers of postnatal mononuclear cardiomyocytes and prolonged cardiomyocyte cell cycle activity, controlled inhibition of DOT1L might be a strategy to promote cardiac regeneration post-injury.
Although vaccines are currently used to control the coronavirus disease 2019 (COVID-19) pandemic, treatment options are urgently needed for those who cannot be vaccinated and for future outbreaks involving new severe acute respiratory syndrome coronavirus virus 2 (SARS-CoV-2) strains or coronaviruses not covered by current vaccines. Thus far, few existing antivirals are known to be effective against SARS-CoV-2 and clinically successful against COVID-19.
As part of an immediate response to the COVID-19 pandemic, a high-throughput, high content imaging–based SARS-CoV-2 infection assay was developed in VeroE6-eGFP cells and was used to screen a library of 5676 compounds that passed phase 1 clinical trials. Eight candidates (nelfinavir, RG-12915, itraconazole, chloroquine, hydroxychloroquine, sematilide, remdesivir, and doxorubicin) with in vitro anti–SARS-CoV-2 activity in VeroE6-eGFP and/or Caco-2 cell lines were identified. However, apart from remdesivir, toxicity and pharmacokinetic data did not support further clinical development of these compounds for COVID-19 treatment.
The change in allele frequencies within a population over time represents a fundamental process of evolution. By monitoring allele frequencies, we can analyze the effects of natural selection and genetic drift on populations. To efficiently track time-resolved genetic change, large experimental or wild populations can be sequenced as pools of individuals sampled over time using high-throughput genome sequencing (called the Evolve & Resequence approach, E&R). Here, we present a set of experiments using hundreds of natural genotypes of the model plant Arabidopsis thaliana to showcase the power of this approach to study rapid evolution at large scale. First, we validate that sequencing DNA directly extracted from pools of flowers from multiple plants -- organs that are relatively consistent in size and easy to sample -- produces comparable results to other, more expensive state-of-the-art approaches such as sampling and sequencing of individual leaves. Sequencing pools of flowers from 25-50 individuals at ∼40X coverage recovers genome-wide frequencies in diverse populations with accuracy r > 0.95. Secondly, to enable analyses of evolutionary adaptation using E&R approaches of plants in highly replicated environments, we provide open source tools that streamline sequencing data curation and calculate various population genetic statistics two orders of magnitude faster than current software. To directly demonstrate the usefulness of our method, we conducted a two-year outdoor evolution experiment with A. thaliana to show signals of rapid evolution in multiple genomic regions. We demonstrate how these laboratory and computational Pool-seq-based methods can be scaled to study hundreds of populations across many climates.
The brains of black 6 mice (Mus musculus) and Seba’s short-tailed bats (Carollia perspicillata) weigh roughly the same and share the mammalian neocortical laminar architecture. Bats have highly developed sonar calls and social communication and are an excellent neuroethological animal model for auditory research. Mice are olfactory and somatosensory specialists and are used frequently in auditory neuroscience, particularly for their advantage of standardization and genetic tools. Investigating their potentially different general auditory processing principles would advance our understanding of how the ecological needs of a species shape the development and function of the mammalian nervous system. We compared two existing datasets, recorded with linear multichannel electrodes down the depth of the primary auditory cortex (A1) while awake, across both species while presenting repetitive stimulus trains with different frequencies (∼5 and ∼40 Hz). We found that while there are similarities between cortical response profiles in bats and mice, there was a better signal to noise ratio in bats under these conditions, which allowed for a clearer following response to stimuli trains. This was most evident at higher frequency trains, where bats had stronger response amplitude suppression to consecutive stimuli. Phase coherence was far stronger in bats during stimulus response, indicating less phase variability in bats across individual trials. These results show that although both species share cortical laminar organization, there are structural differences in relative depth of layers. Better signal to noise ratio in bats could represent specialization for faster temporal processing shaped by their individual ecological niches.
Determining the phase structure of Quantum Chromodynamics (QCD) and its Equation of State (EOS) at densities and temperatures realized inside neutron stars and their mergers is a long-standing open problem. The holographic V-QCD framework provides a model for the EOS of dense and hot QCD, which describes the deconfinement phase transition between a dense baryonic and a quark matter phase. We use this model in fully general relativistic hydrodynamic (GRHD) simulations to study the formation of quark matter and the emitted gravitational wave signal of binary systems that are similar to the first ever observed neutron star merger event GW170817.
The purpose of the paper is to initiate the development of the theory of Newton Okounkov bodies of curve classes. Our denition is based on making a fundamental property of NewtonOkounkov bodies hold also in the curve case: the volume of the NewtonOkounkov body of a curve is a volume-type function of the original curve. This construction allows us to conjecture a new relation between NewtonOkounkov bodies, we prove it in certain cases.
Motivation Expert curation to differentiate between functionally diverged homologs and those that may still share a similar function routinely relies on the visual interpretation of domain architecture changes. However, the size of contemporary data sets integrating homologs from hundreds to thousands of species calls for alternate solutions. Scoring schemes to evaluate domain architecture similarities can help to automatize this procedure, in principle. But existing schemes are often too simplistic in the similarity assessment, many require an a-priori resolution of overlapping domain annotations, and those that allow overlaps to extend the set of annotations sources cannot account for redundant annotations. As a consequence, the gap between the automated similarity scoring and the similarity assessment based on visual architecture comparison is still too wide to make the integration of both approaches meaningful.
Results Here, we present FAS, a scoring system for the comparison of multi-layered feature architectures integrating information from a broad spectrum of annotation sources. Feature architectures are represented as directed acyclic graphs, and redundancies are resolved in the course of comparison using a score maximization algorithm. A benchmark using more than 10,000 human-yeast ortholog pairs reveals that FAS consistently outperforms existing scoring schemes. Using three examples, we show how automated architecture similarity assessments can be routinely applied in the benchmarking of orthology assignment software, in the identification of functionally diverged orthologs, and in the identification of entries in protein collections that most likely stem from a faulty gene prediction.
Tree bark constitutes ideal habitat for microbial communities, because it is a stable substrate, rich in micro-niches. Bacteria, fungi, and terrestrial microalgae together form microbial communities, which in turn support more bark-associated organisms, such as mosses, lichens, and invertebrates, thus contributing to forest biodiversity. We have a limited understanding of the diversity and biotic interactions of the bark-associated microbiome, as investigations have mainly focussed on agriculturally relevant systems and on single taxonomic groups. Here we implemented a multi-kingdom metabarcoding approach to analyse diversity and community structure of the green algal, bacterial, and fungal components of the bark-associated microbial communities of beech, the most common broadleaved tree of Central European forests. We identified the most abundant taxa, hub taxa, and co-occurring taxa. We found that tree size (as a proxy for age) is an important driver of community assembly, suggesting that environmental filtering leads to less diverse fungal and algal communities over time. Conversely, forest management intensity had negligible effects on microbial communities on bark. Our study suggests the presence of undescribed, yet ecologically meaningful taxa, especially in the fungi, and highlights the importance of bark surfaces as a reservoir of microbial diversity. Our results constitute a first, essential step towards an integrated framework for understanding microbial community assembly processes on bark surfaces, an understudied habitat and neglected component of terrestrial biodiversity. Finally, we propose a cost-effective sampling strategy to study bark-associated microbial communities across large spatial or environmental scales.
The pitfalls of measuring representational similarity using representational similarity analysis
(2022)
A core challenge in cognitive and brain sciences is to assess whether different biological systems represent the world in a similar manner. Representational Similarity Analysis (RSA) is an innovative approach to address this problem and has become increasingly popular across disciplines ranging from artificial intelligence to computational neuroscience. Despite these successes, RSA regularly uncovers difficult-to-reconcile and contradictory findings. Here, we demonstrate the pitfalls of using RSA and explain how contradictory findings arise due to false inferences about representational similarity based on RSA-scores. In a series of studies that capture increasingly plausible training and testing scenarios, we compare neural representations in computational models, primate cortex and human cortex. These studies reveal two problematic phenomena that are ubiquitous in current research: a “mimic” effect, where confounds in stimuli can lead to high RSA-scores between provably dissimilar systems, and a “modulation effect”, where RSA-scores become dependent on stimuli used for testing. Since our results bear on a number of influential findings and the inferences drawn by current practitioners in a wide range of disciplines, we provide recommendations to avoid these pitfalls and sketch a way forward to a more solid science of representation in cognitive systems.
The pitfalls of measuring representational similarity using representational similarity analysis
(2022)
A core challenge in neuroscience is to assess whether diverse systems represent the world similarly. Representational Similarity Analysis (RSA) is an innovative approach to address this problem and has become increasingly popular across disciplines from machine learning to computational neuroscience. Despite these successes, RSA regularly uncovers difficult-to-reconcile and contradictory findings. Here we demonstrate the pitfalls of using RSA to infer representational similarity and explain how contradictory findings arise and support false inferences when left unchecked. By comparing neural representations in primate, human and computational models, we reveal two problematic phenomena that are ubiquitous in current research: a “mimic” effect, where confounds in stimuli can lead to high RSA scores between provably dissimilar systems, and a “modulation effect”, where RSA-scores become dependent on stimuli used for testing. Since our results bear on existing findings and inferences, we provide recommendations to avoid these pitfalls and sketch a way forward.
Some pitfalls of measuring representational similarity using Representational Similarity Analysis
(2022)
A core challenge in cognitive and brain sciences is to assess whether different biological systems represent the world in a similar manner. Representational Similarity Analysis (RSA) is an innovative approach that addresses this problem by looking for a second-order isomorphisim in neural activation patterns. This innovation makes it easy to compare latent representations across individuals, species and computational models, and accounts for its popularity across disciplines ranging from artificial intelligence to computational neuroscience. Despite these successes, using RSA has led to difficult-to-reconcile and contradictory findings, particularly when comparing primate visual representations with deep neural networks (DNNs): even though DNNs have been shown to learn and behave in vastly different ways to humans, comparisons based on RSA have shown striking similarities in some studies. Here, we demonstrate some pitfalls of using RSA and explain how contradictory findings can arise due to false inferences about representational similarity based on RSA-scores. In a series of studies that capture increasingly plausible training and testing scenarios, we compare neural representations in computational models, primate cortex and human cortex. These studies reveal two problematic phenomena that are ubiquitous in current research: a “mimic effect”, where confounds in stimuli can lead to high RSA-scores between provably dissimilar systems, and a “modulation effect”, where RSA-scores become dependent on stimuli used for testing. Since our results bear on a number of influential findings, such as comparisons made between human visual representations and those of primates and DNNs, we provide recommendations to avoid these pitfalls and sketch a way forward to a more solid science of representation in cognitive systems.
We deal with the reconstruction of inclusions in elastic bodies based on monotonicity methods and construct conditions under which a resolution for a given partition can be achieved. These conditions take into account the background error as well as the measurement noise. As a main result, this shows us that the resolution guarantees depend heavily on the Lamé parameter μ and only marginally on λ.
Holography has provided valuable insights into the time evolution of strongly coupled gauge theories in a fixed spacetime. However, this framework is insufficient if this spacetime is dynamical. We present a scheme to evolve a four-dimensional, strongly interacting gauge theory coupled to four-dimensional dynamical gravity in the semiclassical regime. As in previous work, we use holography to evolve the quantum gauge theory stress tensor, whereas the four-dimensional metric evolves according to Einstein's equations coupled to the expectation value of the stress tensor. The novelty of our approach is that both the boundary and the bulk spacetimes are constructed dynamically, one time step at a time. We focus on Friedmann-Lemaître-Robertson-Walker geometries and evolve far-from-equilibrium initial states that lead to asymptotically expanding, flat or collapsing Universes.
Quantitative MRI maps of human neocortex explored using cell type-specific gene expression analysis
(2022)
Quantitative MRI (qMRI) allows extraction of reproducible and robust parameter maps. However, the connection to underlying biological substrates remains murky, especially in the complex, densely packed cortex. We investigated associations in human neocortex between qMRI parameters and neocortical cell types by comparing the spatial distribution of the qMRI parameters longitudinal relaxation rate (R1), effective transverse relaxation rate (R2∗), and magnetization transfer saturation (MTsat) to gene expression from the Allen Human Brain Atlas, then combining this with lists of genes enriched in specific cell types found in the human brain. As qMRI parameters are magnetic field strength-dependent, the analysis was performed on MRI data at 3T and 7T. All qMRI parameters significantly covaried with genes enriched in GABA- and glutamatergic neurons, i.e. they were associated with cytoarchitecture. The qMRI parameters also significantly covaried with the distribution of genes enriched in astrocytes (R2∗ at 3T, R1 at 7T), endothelial cells (R1 and MTsat at 3T), microglia (R1 and MTsat at 3T, R1 at 7T), and oligodendrocytes (R1 at 7T). These results advance the potential use of qMRI parameters as biomarkers for specific cell types.
Neural computations emerge from recurrent neural circuits that comprise hundreds to a few thousand neurons. Continuous progress in connectomics, electrophysiology, and calcium imaging require tractable spiking network models that can consistently incorporate new information about the network structure and reproduce the recorded neural activity features. However, it is challenging to predict which spiking network connectivity configurations and neural properties can generate fundamental operational states and specific experimentally reported nonlinear cortical computations. Theoretical descriptions for the computational state of cortical spiking circuits are diverse, including the balanced state where excitatory and inhibitory inputs balance almost perfectly or the inhibition stabilized state (ISN) where the excitatory part of the circuit is unstable. It remains an open question whether these states can co-exist with experimentally reported nonlinear computations and whether they can be recovered in biologically realistic implementations of spiking networks. Here, we show how to identify spiking network connectivity patterns underlying diverse nonlinear computations such as XOR, bistability, inhibitory stabilization, supersaturation, and persistent activity. We established a mapping between the stabilized supralinear network (SSN) and spiking activity which allowed us to pinpoint the location in parameter space where these activity regimes occur. Notably, we found that biologically-sized spiking networks can have irregular asynchronous activity that does not require strong excitation-inhibition balance or large feedforward input and we showed that the dynamic firing rate trajectories in spiking networks can be precisely targeted without error-driven training algorithms.
The exploration of hot and dense nuclear matter: Introduction to relativistic heavy-ion physics
(2022)
This article summarizes our present knowledge about nuclear matter at the highest energy densities and its formation in relativistic heavy ion collisions. We review what is known about the structure and properties of the quark-gluon plasma and survey the observables that are used to glean information about it from experimental data.
We investigate the impact of non-Hermiticity on the thermodynamic properties of interacting fermions by examining bilinear extensions to the 3+1 dimensional SU(2)-symmetric Nambu--Jona-Lasinio (NJL) model of quantum chromodynamics at finite temperature and chemical potential. The system is modified through the anti-PT-symmetric pseudoscalar bilinear ψ¯γ5ψ and the PT-symmetric pseudovector bilinear iBνψ¯γ5γνψ, introduced with a coupling g. Beyond the possibility of dynamical fermion mass generation at finite temperature and chemical potential, our findings establish model-dependent changes in the position of the chiral phase transition and the critical end-point. These are tunable with respect to g in the former case, and both g and |B|/B0 in the latter case, for both lightlike and spacelike fields. Moreover, the behavior of the quark number, entropy, pressure and energy densities signal a potential fermion or antifermion excess compared to the standard NJL model, due to the pseudoscalar and pseudovector extension respectively. In both cases regions with negative interaction measure I=ϵ−3p are found. Future indications of such behaviors in strongly interacting fermion systems, for example in the context of neutron star physics, may point toward the presence of non-Hermitian contributions. These trends provide a first indication of curious potential mechanisms for producing non-Hermitian baryon asymmetry. In addition, the formalism described in this study is expected to apply more generally to other Hamiltonians with four-fermion interactions and thus the effects of the non-Hermitian bilinears are likely to be generic.
We investigate the impact of non-Hermiticity on the thermodynamic properties of interacting fermions by examining bilinear extensions to the 3+1 dimensional SU(2)-symmetric Nambu--Jona-Lasinio (NJL) model of quantum chromodynamics at finite temperature and chemical potential. The system is modified through the anti-PT-symmetric pseudoscalar bilinear ψ¯γ5ψ and the PT-symmetric pseudovector bilinear iBνψ¯γ5γνψ, introduced with a coupling g. Beyond the possibility of dynamical fermion mass generation at finite temperature and chemical potential, our findings establish model-dependent changes in the position of the chiral phase transition and the critical end-point. These are tunable with respect to g in the former case, and both g and |B|/B0 in the latter case, for both lightlike and spacelike fields. Moreover, the behavior of the quark number, entropy, pressure, and energy densities signal a potential fermion or antifermion excess compared to the standard NJL model, due to the pseudoscalar and pseudovector extension respectively. In both cases regions with negative interaction measure I=ϵ−3p are found. Future indications of such behaviors in strongly interacting fermion systems, for example in the context of neutron star physics, may point toward the presence of non-Hermitian contributions. These trends provide a first indication of curious potential mechanisms for producing non-Hermitian baryon asymmetry. In addition, the formalism described in this study is expected to apply more generally to other Hamiltonians with four-fermion interactions and thus the effects of the non-Hermitian bilinears are likely to be generic.
Chern numbers can be calculated within a frame of vortex fields related to phase conventions of a wave function. In a band protected by gaps the Chern number is equivalent to the total number of flux carrying vortices. In the presence of topological defects like Dirac cones this method becomes problematic, in particular if they lack a well-defined winding number. We develop a scheme to include topological defects into the vortex field frame. A winding number is determined by the behavior of the phase in reciprocal space when encircling the defect's contact point. To address the possible lack of a winding number we utilize a more general concept of winding vectors. We demonstrate the usefulness of this ansatz on Dirac cones generated from bands of the Hofstadter model.
Chern numbers can be calculated within a frame of vortex fields related to phase conventions of a wave function. In a band protected by gaps the Chern number is equivalent to the total number of flux carrying vortices. In the presence of topological defects like Dirac cones this method becomes problematic, in particular if they lack a well-defined winding number. We develop a scheme to include topological defects into the vortex field frame. A winding number is determined by the behavior of the phase in reciprocal space when encircling the defect's contact point. To address the possible lack of a winding number we utilize a more general concept of winding vectors. We demonstrate the usefulness of this ansatz on Dirac cones generated from bands of the Hofstadter model.
Abstract
Natural plant populations often harbour substantial heritable variation in DNA methylation. However, a thorough understanding of the genetic and environmental drivers of this epigenetic variation requires large-scale and high-resolution data, which currently exist only for a few model species. Here, we studied 207 lines of the annual weed Thlaspi arvense (field pennycress), collected across a large latitudinal gradient in Europe and propagated in a common environment. By screening for variation in DNA sequence and DNA methylation using whole-genome (bisulfite) sequencing, we found significant epigenetic population structure across Europe. Average levels of DNA methylation were strongly context-dependent, with highest DNA methylation in CG context, particularly in transposable elements and in intergenic regions. Residual DNA methylation variation within all contexts was associated with genetic variants, which often co-localized with annotated methylation machinery genes but also with new candidates. Variation in DNA methylation was also significantly associated with climate of origin, with methylation levels being higher in warmer regions and lower in more variable climates. Finally, we used variance decomposition to assess genetic versus environmental associations with differentially methylation regions (DMRs). We found that while genetic variation was generally the strongest predictor of DMRs, the strength of environmental associations increased from CG to CHG and CHH, with climate-of-origin as the strongest predictor in about one third of the CHH DMRs. In summary, our data show that natural epigenetic variation in Thlaspi arvense is significantly associated with both DNA sequence and environment of origin, and that the relative importance of the two factors strongly depends on the sequence context of DNA methylation. T. arvense is an emerging biofuel and winter cover crop; our results may hence be relevant for breeding efforts and agricultural practices in the context of rapidly changing environmental conditions.
Author Summary: Variation within species is an important level of biodiversity, and it is key for future adaptation. Besides variation in DNA sequence, plants also harbour heritable variation in DNA methylation, and we want to understand the evolutionary significance of this epigenetic variation, in particular how much of it is under genetic control, and how much is associated with the environment. We addressed these questions in a high-resolution molecular analysis of 207 lines of the common plant field pennycress (Thlaspi arvense), which we collected across Europe, propagated under standardized conditions, and sequenced for their genetic and epigenetic variation. We found large geographic variation in DNA methylation, associated with both DNA sequence and climate of origin. Genetic variation was generally the stronger predictor of DNA methylation variation, but the strength of environmental association varied between different sequence contexts. Climate-of-origin was the strongest predictor in about one third of the differentially methylated regions in the CHH context, which suggests that epigenetic variation may play a role in the short-term climate adaptation of pennycress. As pennycress is currently being domesticated as a new biofuel and winter cover crop, our results may be relevant also for agriculture, particularly in changing environments.
The mammalian frontal and auditory cortices are important for vocal behaviour. Here, using local field potential recordings, we demonstrate for the first time that the timing and spatial pattern of oscillations in the fronto-auditory cortical network of vocalizing bats (Carollia perspicillata) predict the purpose of vocalization: echolocation or communication. Transfer entropy analyses revealed predominantly top-down (frontal-to-auditory cortex) information flow during spontaneous activity and pre-vocal periods. The dynamics of information flow depended on the behavioural role of the vocalization and on the timing relative to vocal onset. Remarkably, we observed the emergence of predominantly bottom-up (auditory-to-frontal cortex) information transfer patterns specific echolocation production, leading to self-directed acoustic feedback. Electrical stimulation of frontal areas selectively enhanced responses to echolocation sounds in auditory cortex. These results reveal unique changes in information flow across sensory and frontal cortices, potentially driven by the purpose of the vocalization in a highly vocal mammalian model.
More than 75% of surface and secreted proteins are modified by covalent addition of complex sugars through N- and O-glycosylation. Unlike proteins, glycans do not typically adopt specific secondary structures and remain very mobile, influencing protein dynamics and interactions with other molecules. Glycan conformational freedom impairs complete structural elucidation of glycoproteins. Computer simulations may be used to model glycan structure and dynamics. However, such simulations typically require thousands of computing hours on specialized supercomputers, thus limiting routine use. Here, we describe a reductionist method that can be implemented on personal computers to graft ensembles of realistic glycan conformers onto static protein structures in a matter of minutes. Using this open-source pipeline, we reconstructed the full glycan cover of SARS-CoV-2 Spike protein (S-protein) and a human GABAA receptor. Focusing on S-protein, we show that GlycoSHIELD recapitulates key features of extended simulations of the glycosylated protein, including epitope masking, and provides new mechanistic insights on N-glycan impact on protein structural dynamics.
The human brain achieves visual object recognition through multiple stages of nonlinear transformations operating at a millisecond scale. To predict and explain these rapid transformations, computational neuroscientists employ machine learning modeling techniques. However, state-of-the-art models require massive amounts of data to properly train, and to the present day there is a lack of vast brain datasets which extensively sample the temporal dynamics of visual object recognition. Here we collected a large and rich dataset of high temporal resolution EEG responses to images of objects on a natural background. This dataset includes 10 participants, each with 82,160 trials spanning 16,740 image conditions. Through computational modeling we established the quality of this dataset in five ways. First, we trained linearizing encoding models that successfully synthesized the EEG responses to arbitrary images. Second, we correctly identified the recorded EEG data image conditions in a zero-shot fashion, using EEG synthesized responses to hundreds of thousands of candidate image conditions. Third, we show that both the high number of conditions as well as the trial repetitions of the EEG dataset contribute to the trained models’ prediction accuracy. Fourth, we built encoding models whose predictions well generalize to novel participants. Fifth, we demonstrate full end-to-end training of randomly initialized DNNs that output M/EEG responses for arbitrary input images. We release this dataset as a tool to foster research in visual neuroscience and computer vision.
NAD is a coenzyme central to metabolism that was also found to serve as a 5’-terminal cap of bacterial and eukaryotic RNA species. The presence and functionality of NAD-capped RNAs (NAD-RNAs) in the archaeal domain remain to be characterized in detail. Here, by combining LC-MS and NAD captureSeq methodology, we quantified the total levels of NAD-RNAs and determined the identity of NAD-RNAs in the two model archaea, Sulfolobus acidocaldarius and Haloferax volcanii. A complementary differential RNA-Seq (dRNA-Seq) analysis revealed that NAD transcription start sites (NAD-TSS) correlate with well-defined promoter regions and often overlap with primary transcription start sites (pTSS). The population of NAD-RNAs in the two archaeal organisms shows clear differences, with S. acidocaldarius possessing more capped small non-coding RNAs (sncRNAs) and leader sequences. The NAD-cap did not prevent 5’→3’ exonucleolytic activity by the RNase Saci-aCPSF2. To investigate enzymes that facilitate the removal of the NAD-cap, four Nudix proteins of S. acidocaldarius were screened. None of the recombinant proteins showed NAD decapping activity. Instead, the Nudix protein Saci_NudT5 showed activity after incubating NAD-RNAs at elevated temperatures. Hyperthermophilic environments promote the thermal degradation of NAD into the toxic product ADPR. Incorporating NAD into RNAs and the regulation of ADPR-RNA decapping by Saci_NudT5 is proposed to provide additional layers of maintaining stable NAD levels in archaeal cells.
Importance: This study reports the first characterization of 5’-terminally modified RNA molecules in Archaea and establishes that NAD-RNA modifications, previously only identified in the other two domains of life, are also prevalent in the archaeal model organisms Sulfolobus acidocaldarius and Haloferax volcanii. We screened for NUDIX hydrolases that could remove the NAD-RNA cap and showed that none of these enzymes removed NAD modifications, but we discovered an enzyme that hydrolyzes ADPR-RNA. We propose that these activities influence the stabilization of NAD and its thermal degradation to potentially toxic ADPR products at elevated growth temperatures.
Intrinsically disordered regions (IDRs) are essential for membrane receptor regulation but often remain unresolved in structural studies. TRPV4, a member of the TRP vanilloid channel family involved in thermo- and osmosensation, has a large N-terminal IDR of approximately 150 amino acids. With an integrated structural biology approach, we analyze the structural ensemble of the TRPV4 IDR and identify a network of regulatory elements that modulate channel activity in a hierarchical lipid-dependent manner through transient long-range interactions. A highly conserved autoinhibitory patch acts as a master regulator by competing with PIP2 binding to attenuate channel activity. Molecular dynamics simulations show that loss of the interaction between the PIP2-binding site and the membrane reduces the force exerted by the IDR on the structured core of TRPV4. This work demonstrates that IDR structural dynamics are coupled to TRPV4 activity and highlights the importance of IDRs for TRP channel function and regulation.
The most basic behavioural states of animals can be described as active or passive. However, while high-resolution observations of activity patterns can provide insights into the ecology of animal species, few methods are able to measure the activity of individuals of small taxa in their natural environment. We present a novel approach in which the automated VHF radio-tracking of small vertebrates fitted with lightweight transmitters (< 0.2 g) is used to distinguish between active and passive behavioural states.
A dataset containing > 3 million VHF signals was used to train and test a random forest model in the assignment of either active or passive behaviour to individuals from two forest-dwelling bat species (Myotis bechsteinii (n = 50) and Nyctalus leisleri (n = 20)). The applicability of the model to other taxonomic groups was demonstrated by recording and classifying the behaviour of a tagged bird and by simulating the effect of different types of vertebrate activity with the help of humans carrying transmitters. The random forest model successfully classified the activity states of bats as well as those of birds and humans, although the latter were not included in model training (F-score 0.96–0.98).
The utility of the model in tackling ecologically relevant questions was demonstrated in a study of the differences in the daily activity patterns of the two bat species. The analysis showed a pronounced bimodal activity distribution of N. leisleri over the course of the night while the night-time activity of M. bechsteinii was relatively constant. These results show that significant differences in the timing of species activity according to ecological preferences or seasonality can be distinguished using our method.
Our approach enables the assignment of VHF signal patterns to fundamental behavioural states with high precision and is applicable to different terrestrial and flying vertebrates. To encourage the broader use of our radio-tracking method, we provide the trained random forest models together with an R-package that includes all necessary data-processing functionalities. In combination with state-of-the-art open-source automated radio-tracking, this toolset can be used by the scientific community to investigate the activity patterns of small vertebrates with high temporal resolution, even in dense vegetation.
Temporal anticipation is a fundamental process underlying complex neural functions such as associative learning, decision-making, and motor-preparation. Here we study event anticipation in its simplest form in human participants using magnetoencephalography. We distributed events in time according to different probability density functions and presented the stimuli separately in two different sensory modalities. We found that the temporal dynamics in right parietal cortex correlate with reaction times to anticipated events. Specifically, after an event occurred, event probability was represented in right parietal activity, hinting at a functional role of event-related potential component P300 in temporal expectancy. The results are consistent across both visual and auditory modalities. The right parietal cortex seems to play a central role in the processing of event probability density. Overall, this work contributes to the understanding of the neural processes involved in the anticipation of events in time.
Using the notion of a root datum of a reductive group G we propose a tropical analogue of a principal G-bundle on a metric graph. We focus on the case G=GLn, i.e. the case of vector bundles. Here we give a characterization of vector bundles in terms of multidivisors and use this description to prove analogues of the Weil--Riemann--Roch theorem and the Narasimhan--Seshadri correspondence. We proceed by studying the process of tropicalization. In particular, we show that the non-Archimedean skeleton of the moduli space of semistable vector bundles on a Tate curve is isomorphic to a certain component of the moduli space of semistable tropical vector bundles on its dual metric graph.
Candida boidinii NAD+-dependent formate dehydrogenase (CbFDH) has gained significant attention for its potential applications in the production of biofuels and various industrial chemicals from inorganic carbon dioxide. The present study reports the atomic X-ray crystal structures of the wild-type CbFDH at cryogenic and ambient temperatures as well as Val120Thr mutant at cryogenic temperature determined at the Turkish Light Source "Turkish DeLight". The structures reveal new hydrogen bonds between Thr120 and water molecules in the mutant CbFDH's active site, suggesting increased stability of the active site and more efficient electron transfer during the reaction. Further experimental data is needed to test these hypotheses. Collectively, our findings provide invaluable insights into future protein engineering efforts that could potentially enhance the efficiency and effectiveness of CbFDH.
Background Overweight and decreased physical fitness are highly prevalent in schizophrenia, represent a major risk factor for cardio-vascular diseases and decrease the patients’ life expectancies. It is thus important to understand the underlying mechanisms that link psychopathology and weight gain. We hypothesize that the dopaminergic reward system plays an important role in this.
Methods: We analyzed the seed-based functional connectivity (FC) of the ventral tegmental area (VTA) in a group of schizophrenic patients (n = 32) and age- as well as gender matched healthy controls (n = 27). We then correlated the resting-state results with physical fitness parameters, obtained in a fitness test, and psychopathology.
Results: The seed-based connectivity analysis revealed decreased functional connections between the VTA and the anterior cingulate cortex (ACC), as well as the dorsolateral prefrontal cortex and increased functional connectivity between the VTA and the middle temporal gyrus in patients compared to healthy controls. The decreased FC between the VTA and the ACC of the patient group could further be associated with increased body fat and negatively correlated with the overall physical fitness. We found no significant correlations with psychopathology.
Conclusion: Although we did not find significant correlations with psychopathology, we could link decreased physical fitness and high body fat with dysconnectivity between the VTA and the ACC in schizophrenia. These findings demonstrate that a dysregulated reward system is not just responsible for symptomatology in schizophrenia but is also involved in comorbidities and could pave the way for future lifestyle therapy interventions.
The COVID-19 pandemic and the associated prevention measures did not only impact on the transmission of COVID-19 but also on the spread of other infectious diseases in an unprecedented natural experiment. Here, we analysed the transmission patterns of 22 different infectious diseases during the COVID-19 pandemic in England. Our results show that the COVID-19 prevention measures generally reduced the spread of pathogens that are transmitted via the air and the faecal-oral route. Moreover, the COVID-19 prevention measures resulted in the sustained suppression of vaccine-preventable infectious diseases also after the removal of restrictions, while non-vaccine preventable diseases displayed a rapid rebound. Despite concerns that a lack of exposure to common pathogens may affect population immunity and result in large outbreaks by various pathogens post-COVID-19, only four of the 22 investigated diseases and disease groups displayed higher post-than pre-pandemic levels without an obvious causative relationship. Notably, this included chickenpox for which an effective vaccine is available but not used in the UK, which provides strong evidence supporting the inclusion of the chickenpox vaccination into the routine vaccination schedule in the UK. In conclusion, our findings provide unique, novel insights into the impact of non-pharmaceutical interventions on the spread of a broad range of infectious diseases.
The NVX-CoV2373-vaccine has recently been licensed, although data on vaccine-induced humoral and cellular immunity towards the parental strain and variants of concern (VOCs) in comparison to dual-dose mRNA-regimens are limited. In this observational study including 66 participants, we show that NVX-CoV2373-induced IgG-levels were lower than after vaccination with BNT162b2 or mRNA-1273 (n=22 each, p=0.006). Regardless of the vaccine and despite different IgG-levels, neutralizing activity towards VOCs was highest for Delta, followed by BA.2 and BA.1. Interestingly, spike-specific CD8 T-cell levels after NVX-CoV2373-vaccination were significantly lower and were detectable in 3/22 (14%) individuals only. In contrast, spike-specific CD4 T-cells were induced in 18/22 (82%) individuals. However, CD4 T-cell levels were lower (p<0.001), had lower CTLA-4 expression (p<0.0001) and comprised less multifunctional cells co-expressing IFNγ, TNFαα and IL-2 (p=0.0007) as compared to mRNA-vaccinated individuals. Unlike neutralizing antibodies, NVX-CoV2373-induced CD4 T cells cross-reacted to all tested VOCs from Alpha to Omicron, which may hold promise to protect from severe disease.
Electroencephalography (EEG) has been used for decades to identify neurocognitive processes related to intelligence. Evidence is accumulating for associations with neural markers of higher-order cognitive processes (e.g., working memory); however, whether associations are specific to complex processes or also relate to earlier processing stages remains unclear. Addressing these issues has implications for improving our understanding of intelligence and its neural correlates. The mismatch negativity (MMN) is an event-related brain potential (ERP) that is elicited when, within a series of frequent standard stimuli, rare deviant stimuli are presented. As stimuli are typically presented outside the focus of attention, the MMN is suggested to capture automatic pre-attentive discrimination processes. However, the MMN and its relation to intelligence has largely only been studied in the auditory domain, thus preventing conclusions about the involvement of automatic discrimination processes in humans’ dominant sensory modality vision. Electroencephalography was recorded from 50 healthy participants during a passive visual oddball task that presented simple sequence violations as well as deviations within a more complex hidden pattern. Signed area amplitudes and fractional area latencies of the visual mismatch negativity (vMMN) were calculated with and without Laplacian transformation. Correlations between vMMN and intelligence (Raven’s Advanced Progressive Matrices) were of negligible to small effect sizes, differed critically between measurement approaches, and Bayes Factors provided anecdotal to substantial evidence for the absence of an association. We discuss differences between the auditory and visual MMN, the implications of different measurement approaches, and offer recommendations for further research in this evolving field.
Effective spectral functions of the ρ meson are reconstructed by considering the lifetimes inside different media using the hadronic transport SMASH (Simulating Many Accelerated Strongly-interacting Hadrons). Due to inelastic scatterings, resonance lifetimes are dynamically shortened (collisional broadening), even though the employed approach assumes vacuum resonance properties. Analyzing the ρ meson lifetimes allows to quantify an effective broadening of the decay width and spectral function, which is important in order to distinguish dynamical effects from additional genuine medium modifications to the spectral functions, indicating e.g. an onset of chiral symmetry restoration. The broadening of the spectral function in a thermalized system is shown to be consistent with other theoretical calculations. The effective ρ meson spectral function is also presented for the dynamical evolution of heavy-ion collisions, finding a clear correlation of the broadening to system size, which is explained by an observed dependence of the width on the local hadron density. Furthermore, the difference in the results between the thermal system and full collision dynamics is explored, which may point to non-equilibrium effects.
Reactive oxygen species (ROS) are constant by-products of aerobic life. In excess, ROS lead to cytotoxic protein aggregates, which are a hallmark of ageing in animals and linked to age-related pathologies in humans. Acylamino acid-releasing enzymes (AARE) are bifunctional serine proteases, acting on oxidized proteins. AARE are found in all domains of life, albeit under different names, such as acylpeptide hydrolase (APEH/ACPH), acylaminoacyl peptidase (AAP), or oxidized protein hydrolase (OPH). In humans, AARE malfunction is associated with age-related pathologies, while their function in plants is less clear. Here, we provide a detailed analysis of AARE genes in the plant lineage and an in-depth analysis of AARE localization and function in the moss Physcomitrella and the angiosperm Arabidopsis. AARE loss-of-function mutants have not been described for any organism so far. We generated and analysed such mutants and describe a connection between AARE function, aggregation of oxidized proteins and plant ageing, including accelerated developmental progression and reduced life span. Our findings complement similar findings in animals and humans, and support a unified concept of ageing.
Reactive oxygen species (ROS) are constant by-products of aerobic life. In excess, ROS lead to cytotoxic protein aggregates, which are a hallmark of ageing in animals and linked to age-related pathologies in humans. Acylamino acid-releasing enzymes (AARE) are bifunctional serine proteases, acting on oxidized proteins. AARE are found in all domains of life, albeit under different names, such as acylpeptide hydrolase (APEH/ACPH), acylaminoacyl peptidase (AAP), or oxidized protein hydrolase (OPH). In humans, AARE malfunction is associated with age-related pathologies, while their function in plants is less clear. Here, we provide a detailed analysis of AARE genes in the plant lineage and an in-depth analysis of AARE localization and function in the moss Physcomitrella and the angiosperm Arabidopsis. AARE loss-of-function mutants have not been described for any organism so far. We generated and analysed such mutants and describe a connection between AARE function, aggregation of oxidized proteins and plant ageing, including accelerated developmental progression and reduced life span. Our findings complement similar findings in animals and humans, and suggest a unified concept of ageing may exist in different life forms.
Protein oxidation results from the reaction of amino-acid side chains with reactive oxygen species (ROS) and is partly irreversible. In non-photosynthetic tissues, mitochondria are a main source of ROS, whereas plastids are the major source in photosynthetic tissues. Oxidized proteins suffer from decreased structural integrity and even loss of function, and their accumulation leads to cytotoxic aggregates. In mammals, aggregate formation correlates with aging and is linked to several age-related pathologies. Mammalian proteolytic pathways for clearance of oxidized proteins are under intensive research, while mechanistic insights into this process in plants is scarce. Acylamino acid-releasing (AARE) enzymes are ATP-independent serine proteases, presumably acting on oxidized proteins and operating in a dual exo-/endopeptidase mode. They are found in all domains of life. Here, we investigated AARE enzymes in the moss Physcomitrella and the angiosperm Arabidopsis and identified three homologous nuclear genes in Physcomitrella (PpAARE1-3) and a single nuclear gene in Arabidopsis (AtAARE). Surprisingly, we observed triple localization of the proteins AtAARE and PpAARE1 to plastids, mitochondria and the cytosol in vivo, likely conserved across the plant lineage. This represents an ATP-independent possibility for degradation of oxidized proteins in the major source organelles of ROS in plants, which is distinct to mammals. Combinatorial knockout plants and protein interaction analysis revealed specific interactions of the moss AARE isoforms and functions in progressive aging. Analysis of an AtAARE T-DNA mutant further suggests the evolutionary conservation of AARE function in age-related development.
Background Eukaryotic gene expression is controlled by cis-regulatory elements (CREs) including promoters and enhancers which are bound by transcription factors (TFs). Differential expression of TFs and their putative binding sites on CREs cause tissue and developmental-specific transcriptional activity. Consolidating genomic data sets can offer further insights into the accessibility of CREs, TF activity, and thus gene regulation. However, the integration and analysis of multi-modal data sets are hampered by considerable technical challenges. While methods for highlighting differential TF activity from combined ChIP-seq and RNA-seq data exist, they do not offer good usability, have limited support for large-scale data processing, and provide only minimal functionality for visual result interpretation.
Results We developed TF-Prioritizer, an automated java pipeline to prioritize condition-specific TFs derived from multi-modal data. TF-Prioritizer creates an interactive, feature-rich, and user-friendly web report of its results. To showcase the potential of TF-Prioritizer, we identified known active TFs (e.g., Stat5, Elf5, Nfib, Esr1), their target genes (e.g., milk proteins and cell-cycle genes), and newly classified lactating mammary gland TFs (e.g., Creb1, Arnt).
Conclusion TF-Prioritizer accepts ChIP-seq and RNA-seq data, as input and suggests TFs with differential activity, thus offering an understanding of genome-wide gene regulation, potential pathogenesis, and therapeutic targets in biomedical research.
Targeted protein degradation is a drug modality represented by compounds that recruit a target to an E3 ubiquitin ligase to promote target ubiquitination and proteasomal degradation. Historically, the field distinguishes monovalent degraders from bifunctional degraders (PROTACs) that connect target and ligase via separate binding ligands joined via a linker1–4. Here, we elucidate the mechanism of action of a PROTAC-like degrader of the transcriptional coactivator BRD4, composed of a BRD4 ligand linked to a ligand for the E3 ligase CRL4DCAF15. Using orthogonal CRISPR/Cas9 screens we identify the degrader activity is independent of DCAF15, and relies on a different CRL4 substrate receptor, DCAF16. We demonstrate an intrinsic affinity between BRD4 and DCAF16, which is dependent on the tandem bromodomains of BRD4 and further increased by the degrader without physically engaging DCAF16 in isolation. Structural characterization of the resulting ternary complex reveals both BRD4 bromodomains are bivalently engaged in cis by the degrader and are bound to DCAF16 through several interfacial BRD4-DCAF16 and degrader-DCAF16 contacts. Our findings demonstrate that intramolecularly bridging domains can confer glue-type stabilization of intrinsic target-E3 interactions, and we propose this as a general strategy to modulate the surface topology of target proteins to nucleate co-opting of E3 ligases or other cellular effector proteins for effective proximity-based pharmacology.
Glutathione (GSH) is the main determinant of intracellular redox potential and participates in multiple cellular signaling pathways. Achieving a detailed understanding of intracellular GSH trafficking and regulation depends on the development of tools to map GSH compartmentalization and intra-organelle fluctuations. Herein, we present a new GSH sensing platform, TRaQ-G, for live-cell imaging. This small-molecule/protein hybrid sensor possesses a unique reactivity turn-on mechanism that ensures that the small molecule is only sensitive to GSH in the desired location. Furthermore, TRaQ-G can be fused to a fluorescent protein of choice to give a ratiometric response. Using TRaQ-G-mGold, we demonstrated that the nuclear and cytosolic GSH pools are independently regulated during cell proliferation. We also used this sensor, in combination with roGFP, to quantify redox potential and GSH concentration simultaneously in the endoplasmic reticulum. Finally, by exchanging the fluorescent protein, we created a near-infrared, targetable and quantitative GSH sensor.
Tracking influenza a virus infection in the lung from hematological data with machine learning
(2022)
The tracking of pathogen burden and host responses with minimal-invasive methods during respiratory infections is central for monitoring disease development and guiding treatment decisions. Utilizing a standardized murine model of respiratory Influenza A virus (IAV) infection, we developed and tested different supervised machine learning models to predict viral burden and immune response markers, i.e. cytokines and leukocytes in the lung, from hematological data. We performed independently in vivo infection experiments to acquire extensive data for training and testing purposes of the models. We show here that lung viral load, neutrophil counts, cytokines like IFN-γ and IL-6, and other lung infection markers can be predicted from hematological data. Furthermore, feature analysis of the models shows that blood granulocytes and platelets play a crucial role in prediction and are highly involved in the immune response against IAV. The proposed in silico tools pave the path towards improved tracking and monitoring of influenza infections and possibly other respiratory infections based on minimal-invasively obtained hematological parameters.
We have investigated the systematic differences introduced when performing a Bayesian-inference analysis of the equation of state of neutron stars employing either variable- or constant-likelihood functions. The former have the advantage that it retains the full information on the distributions of the measurements, making an exhaustive usage of the data. The latter, on the other hand, have the advantage of a much simpler implementation and reduced computational costs. In both approaches, the EOSs have identical priors and have been built using the sound-speed parameterization method so as to satisfy the constraints from X-ray and gravitationalwaves observations, as well as those from Chiral Effective Theory and perturbative QCD. In all cases, the two approaches lead to very similar results and the 90%-confidence levels are essentially overlapping. Some differences do appear, but in regions where the probability density is extremely small and are mostly due to the sharp cutoff set on the binary tidal deformability Λ˜ ≤ 720 employed in the constant-likelihood analysis. Our analysis has also produced two additional results. First, a clear inverse correlation between the normalized central number density of a maximally massive star, nc,TOV/ns, and the radius of a maximally massive star, RTOV. Second, and most importantly, it has confirmed the relation between the chirp mass Mchirp and the binary tidal deformability Λ˜. The importance of this result is that it relates a quantity that is measured very accurately, Mchirp, with a quantity that contains important information on the micro-physics, Λ˜. Hence, once Mchirp is measured in future detections, our relation has the potential of setting tight constraints on Λ˜.
Viewpoint effects on object recognition interact with object-scene consistency effects. While recognition of objects seen from “accidental” viewpoints (e.g., a cup from below) is typically impeded compared to processing of objects seen from canonical viewpoints (e.g., the string-side of a guitar), this effect is reduced by meaningful scene context information. In the present study we investigated if these findings established by using photographic images, generalise to 3D models of objects. Using 3D models further allowed us to probe a broad range of viewpoints and empirically establish accidental and canonical viewpoints. In Experiment 1, we presented 3D models of objects from six different viewpoints (0°, 60°, 120°, 180° 240°, 300°) in colour (1a) and grayscaled (1b) in a sequential matching task. Viewpoint had a significant effect on accuracy and response times. Based on the performance in Experiments 1a and 1b, we determined canonical (0°-rotation) and non-canonical (120°-rotation) viewpoints for the stimuli. In Experiment 2, participants again performed a sequential matching task, however now the objects were paired with scene backgrounds which could be either consistent (e.g., a cup in the kitchen) or inconsistent (e.g., a guitar in the bathroom) to the object. Viewpoint interacted significantly with scene consistency in that object recognition was less affected by viewpoint when consistent scene information was provided, compared to inconsistent information. Our results show that viewpoint-dependence and scene context effects generalize to depth rotated 3D objects. This supports the important role object-scene processing plays for object constancy.
In a dynamic environment, the already limited information that human working memory can maintain needs to be constantly updated to optimally guide behaviour. Indeed, previous studies showed that working memory representations are continuously being transformed during delay periods leading up to a response. This goes hand-in-hand with the removal of task-irrelevant items. However, does such removal also include veridical, original stimuli, as they were prior to transformation? Here we aimed to assess the neural representation of task-relevant transformed representations, compared to the no-longer-relevant veridical representations they originated from. We applied multivariate pattern analysis to electroencephalographic data during maintenance of orientation gratings with and without mental rotation. During maintenance, we perturbed the representational network by means of a visual impulse stimulus, and were thus able to successfully decode veridical as well as imaginary, transformed orientation gratings from impulse-driven activity. On the one hand, the impulse response reflected only task-relevant (cued), but not task-irrelevant (uncued) items, suggesting that the latter were quickly discarded from working memory. By contrast, even though the original cued orientation gratings were also no longer task-relevant after mental rotation, these items continued to be represented next to the rotated ones, in different representational formats. This seemingly inefficient use of scarce working memory capacity was associated with reduced probe response times and may thus serve to increase precision and flexibility in guiding behaviour in dynamic environments.
Background: Blood donation saves lives. Provided they are in good health, male volunteers can donate as often as six times per year from the age of 18 into their late sixties. The burden of blood donation is very unevenly distributed, with a small minority of altruistic individuals providing this critical resource. While the consequences of persistent iron depletion in blood donors have been studied in the context of cancer and coronary heart disease, potential effects of the erythropoietic stress from repetitive large-volume phlebotomy remain unexplored. We sought to investigate if and how repeated blood donations affect the clonal composition of the hematopoietic stem and progenitor cell (HSPC) compartment.
Methods: 105 healthy, male individuals with an extensive blood donation history (median of 120 donations per donor; median age of 66 yrs.) were screened for the presence of clonal hematopoiesis (CH) using a sequencing panel covering 141 genes commonly mutated in human myeloid neoplasms. The control cohort consisted of 103 healthy, male donors with a median of 5 donations per donor and a median age of 63. Donors positive for CH were subsequently studied longitudinally. The pathogenicity of detected variants was compared using established scoring systems. Finally, to assess the functional consequences of blood-donation induced CH, selected CH mutations were introduced by CRISPR-mediated editing into HSPCs from human cord blood (CB) or bone marrow (BM). The effect of these mutations was tested under different stress stimuli using functional ex vivo long-term culture initiating cells (LTC-IC) assays.
Results: Compared to the control cohort, frequent donors were significantly more likely to have mutations in genes encoding for epigenetic modifiers (44.7 vs. 22.3 %), most specifically in the two genes most commonly mutated in CH, DNMT3A and TET2 (35.2 vs. 20.3 %). However, no difference in the variant allele frequency (VAF) of detected mutations was found between the groups. Longitudinal analysis revealed that the majority of the mutations remained at a stable VAF over an observation period of approximately one year. Three DNMT3A variants from the frequent donor cohort were introduced into healthy HSPCs and functionally analyzed: All expanded in response to EPO, but none responded to LPS or IFNγ stimulation. This contrasted with the leukemogenic DNMT3A R882H mutation, which did not expand in the presence of EPO but instead responded strongly to inflammatory stimuli.
Conclusions: Frequent whole blood donation is associated with a higher prevalence of CH driven by mutations in genes encoding for epigenetic modifiers, with DNMT3A and TET2 being the most common. This increased CH prevalence is not associated with a higher pathogenicity of the associated variants and is likely a result of the selection of clones with improved responsiveness to EPO under the condition of bleeding stress. Our data show that even highly frequent blood donations over many years is not increasing the risk for malignant clones further underscoring the safety of repetitive blood donations. To our knowledge, this is the first CH study analyzing a cohort of individuals known for their superior health and survival, able to donate blood until advanced age. Thus, our analysis possibly identified mutations associated with beneficial outcomes, rather than a disease condition, such as mutations in DNMT3A that mediated the improved expansion of HSPCs in EPO enriched environments. Our data support the notion of ongoing Darwinian evolution in humans at the somatic stem cell level and present EPO as one of the environmental factors to which HSPCs with specific mutations may respond with superior fitness.
Selective attention implements preferential routing of attended stimuli, likely through increasing the influence of the respective synaptic inputs on higher-area neurons. As the inputs of competing stimuli converge onto postsynaptic neurons, presynaptic circuits might offer the best target for attentional top-down influences. If those influences enabled presynaptic circuits to selectively entrain postsynaptic neurons, this might lead to selective routing. Indeed, when two visual stimuli induce two gamma rhythms in V1, only the gamma induced by the attended stimulus entrains gamma in V4. Here, we modeled this selective entrainment with a Dynamic Causal Model for Cross-Spectral Densities and found that it can be explained by attentional modulation of intrinsic V1 connections. Specifically, local inhibition was decreased in the granular input layer and increased in the supragranular output layer of the V1 circuit that processed the attended stimulus. Thus, presynaptic attentional influences and ensuing entrainment were sufficient to mediate selective routing.
RNA-binding proteins (RBPs) control every RNA metabolic process by multiple protein-RNA and protein-protein interactions. Their roles have largely been analyzed by crude mutations, which abrogate multiple functions at once and likely impact the structural integrity of the large messenger ribonucleoprotein particle (mRNP) assemblies, these proteins often function in. Using UV-induced RNA-protein crosslinking and subsequent mass spectrometric analysis, we first identified more than 100 in vivo RNA crosslinks in 16 nuclear mRNP components in S. cerevisiae. For functional analysis, we chose Npl3, for which we determined crosslinks in its two RNA recognition motifs (RRM) and in the flexible linker region connecting the two. Using NMR and structural analyses, we show that both RRM domains and the linker uniquely contribute to RNA recognition. Interestingly, mutations in these regions cause different phenotypes, indicating distinct functions of the different RNA-binding domains of Npl3. Notably, the npl3-Linker mutation strongly impairs recruitment of several mRNP components to chromatin and incorporation of further mRNP components into nuclear mRNPs, establishing a function of Npl3 in nuclear mRNP assembly. Taken together, we determined the specific function of the RNA-binding activity of the nuclear mRNP component Npl3, an approach that can be applied to many RBPs in any RNA metabolic process.
Vaccines are central to controlling the coronavirus disease 2019 (COVID-19) pandemic but the durability of protection is limited for currently approved COVID-19 vaccines. Further, the emergence of variants of concern (VoCs) that evade immune recognition has reduced vaccine effectiveness, compounding the problem. Here, we show that a single dose of a murine cytomegalovirus (MCMV)-based vaccine, which expresses the spike (S) protein of the virus circulating early in the pandemic (MCMVS), protects highly susceptible K18-hACE2 mice from clinical symptoms and death upon challenge with a lethal dose of D614G SARS-CoV-2. Moreover, MCMVS vaccination controlled two immune-evading VoCs, the Beta (B.1.135) and the Omicron (BA.1) variants in BALB/c mice, and S-specific immunity was maintained for at least 5 months after immunization, where neutralizing titers against all tested VoCs were higher at 5-months than at 1-month post-vaccination. Thus, cytomegalovirus (CMV)-based vector vaccines might allow for long-term protection against COVID-19.
The free energy of TAP-solutions for the SK-model of mean field spin glasses can be expressed as a nonlinear functional of local terms: we exploit this feature in order to contrive abstract REM-like models which we then solve by a classical large deviations treatment. This allows to identify the origin of the physically unsettling quadratic (in the inverse of temperature) correction to the Parisi free energy for the SK-model, and formalizes the true cavity dynamics which acts on TAP-space, i.e. on the space of TAP-solutions. From a non-spin glass point of view, this work is the first in a series of refinements which addresses the stability of hierarchical structures in models of evolving populations.
The Kitaev material α-RuCl3 is among the most prominent candidates to host a quantum spin-liquid state endowed with fractionalized excitations. Recent experimental and theoretical investigations have separately revealed the importance of both the magnetoelastic coupling and the magnetic anisotropy, in dependence of the applied magnetic field direction. In this combined theoretical and experimental research, we investigate the anisotropic magnetic and magnetoelastic properties for magnetic fields applied along the main crystallographic axes as well as for fields canted out of the honeycomb plane. We found that the magnetostriction anisotropy is unusually large compared to the anisotropy of the magnetization, which is related to the strong magnetoelastic Γ′˜-type coupling in our \textit{ab-initio} derived model. We observed large, non-symmetric magnetic anisotropy for magnetic fields canted out of the honeycomb ab-plane in opposite directions, namely towards the +c∗ or −c∗ axes, respectively. The observed directional anisotropy is explained by considering the relative orientation of the magnetic field with respect to the co-aligned RuCl6 octahedra. Magnetostriction measurements in canted fields support this non-symmetric magnetic anisotropy, however these experiments are affected by magnetic torque effects. Comparison of theoretical predictions with experimental findings allow us to recognize the significant contribution of torque effects in experimental setups where α-RuCl3 is placed in canted magnetic fields.
The phase diagram of the (1+1)-dimensional Gross-Neveu model is reanalyzed for (non-)zero chemical potential and (non-)zero temperature within the mean-field approximation. By investigating the momentum dependence of the bosonic two-point function, the well-known second-order phase transition from the Z2 symmetric phase to the so-called inhomogeneous phase is detected. In the latter phase the chiral condensate is periodically varying in space and translational invariance is broken. This work is a proof of concept study that confirms that it is possible to correctly localize second-order phase transition lines between phases without condensation and phases of spatially inhomogeneous condensation via a stability analysis of the homogeneous phase. To complement other works relying on this technique, the stability analysis is explained in detail and its limitations and successes are discussed in context of the Gross-Neveu model. Additionally, we present explicit results for the bosonic wave-function renormalization in the mean-field approximation, which is extracted analytically from the bosonic two-point function. We find regions -- a so-called moat regime -- where the wave function renormalization is negative accompanying the inhomogeneous phase as expected.
The phase diagram of the (1+1)-dimensional Gross-Neveu model is reanalyzed for (non-)zero chemical potential and (non-)zero temperature within the mean-field approximation. By investigating the momentum dependence of the bosonic two-point function, the well-known second-order phase transition from the Z2 symmetric phase to the so-called inhomogeneous phase is detected. In the latter phase the chiral condensate is periodically varying in space and translational invariance is broken. This work is a proof of concept study that confirms that it is possible to correctly localize second-order phase transition lines between phases without condensation and phases of spatially inhomogeneous condensation via a stability analysis of the homogeneous phase. To complement other works relying on this technique, the stability analysis is explained in detail and its limitations and successes are discussed in context of the Gross-Neveu model. Additionally, we present explicit results for the bosonic wave-function renormalization in the mean-field approximation, which is extracted analytically from the bosonic two-point function. We find regions -- a so-called moat regime -- where the wave function renormalization is negative accompanying the inhomogeneous phase as expected.
The family of scaffold attachment factor B (SAFB) proteins comprises three members and was first identified as binders of the nuclear matrix/scaffold. Over the past two decades, SAFBs were shown to act in DNA repair, mRNA/(l)ncRNA processing, and as part of protein complexes with chromatin-modifying enzymes. SAFB proteins are approximately-100-kDa-sized dual nucleic acid-binding proteins with dedicated domains in an otherwise largely unstructured context, but whether and how they discriminate DNA- and RNA-binding has remained enigmatic. We here provide the SAFB2 DNA- and RNA-binding SAP and RRM domains in their functional boundaries and use solution NMR spectroscopy to ascribe DNA- and RNA-binding functions. We give insight into their target nucleic acid preferences and map the interfaces with respective nucleic acids on sparse data-derived SAP and RRM domain structures. Further, we provide evidence that the SAP domain exhibits intra-domain dynamics and a potential tendency to dimerise, which may expand its specifically targeted DNA sequence range. Our data provide a first molecular basis of and a starting point towards deciphering DNA- and RNA-binding functions of SAFB2 on the molecular level and serve a basis for understanding its localization to specific regions of chromatin and its involvement in the processing of specific RNA species.
Background: Driven by globalization, urbanization and climate change, the distribution range of invasive vector species has expanded to previously colder ecoregions. To reduce health-threatening impacts on humans, insect vectors are extensively studied. Population genomics can reveal the genomic basis of adaptation and help to identify emerging trends of vector expansion.
Results: By applying whole genome analyses and genotype-environment associations to populations of the main dengue vector Ae. aegypti, sampled along an altitudinal temperature gradient in Nepal (200- 1300m), we identify adaptive traits and describe the species’ genomic footprint of climate adaptation to colder ecoregions. We found two clusters of differentiation with significantly different allele frequencies in genes associated to climate adaptation between the highland population (1300m) and all other lowland populations (≤ 800 m). We revealed non-synonymous mutations in 13 of the candidate genes associated to either altitude, precipitation or cold tolerance and identified an isolation-by-environment differentiation pattern.
Conclusion: Other than the expected gradual differentiation along the altitudinal gradient, our results reveal a distinct genomic differentiation of the highland population. This finding either indicates a differential invasion history to Nepal or local high-altitude adaptation explaining the population’s phenotypic cold tolerance. In any case, this highland population can be assumed to carry pre-adapted alleles relevant for the species’ invasion into colder ecoregions worldwide that way expanding their climate niche.
Objects that are congruent with a scene are recognised more efficiently than objects that are incongruent. Further, semantic integration of incongruent objects elicits a stronger N300/N400 EEG component. Yet, the time course and mechanisms of how contextual information supports access to semantic object information is unclear. We used computational modelling and EEG to test how context influences semantic object processing. Using representational similarity analysis, we established that EEG patterns dissociated between objects in congruent or incongruent scenes from around 300 ms. By modelling semantic processing of objects using independently normed properties, we confirm that the onset of semantic processing of both congruent and incongruent objects is similar (∼150 ms). Critically, after ∼275 ms, we discover a difference in the duration of semantic integration, lasting longer for incongruent compared to congruent objects. These results constrain our understanding of how contextual information supports access to semantic object information.
Respiratory complex I in mitochondria and bacteria catalyzes the transfer of electrons from NADH to quinone (Q). The free energy available from the reaction is used to pump protons and to establish a membrane proton electrochemical gradient, which drives ATP synthesis. Even though several high-resolution structures of complex I have been resolved, how Q reduction is linked with proton pumping, remains unknown. Here, microsecond long molecular dynamics (MD) simulations were performed on Yarrowia lipolytica complex I structures where Q molecules have been resolved in the ~30 Å long Q tunnel. MD simulations of several different redox/protonation states of Q reveal the coupling between the Q dynamics and the restructuring of conserved loops and ion pairs. Oxidized quinone stabilizes towards the N2 FeS cluster, a binding mode not previously described in Yarrowia lipolytica complex I structures. On the other hand, reduced (and protonated) species tend to diffuse towards the Q binding sites closer to the tunnel entrance. Mechanistic and physiological relevance of these results are discussed.
Neuroscience studies in non-human primates (NHP) often follow the rule of thumb that results observed in one animal must be replicated in at least one other. However, we lack a statistical justification for this rule of thumb, or an analysis of whether including three or more animals is better than including two. Yet, a formal statistical framework for experiments with few subjects would be crucial for experimental design, ethical justification, and data analysis. Also, including three or four animals in a study creates the possibility that the results observed in one animal will differ from those observed in the others: we need a statistically justified rule to resolve such situations. Here, I present a statistical framework to address these issues. This framework assumes that conducting an experiment will produce a similar result for a large proportion of the population (termed ‘representative’), but will produce spurious results for a substantial proportion of animals (termed ‘outliers’); the fractions of ‘representative’ and ‘outliers’ animals being defined by a prior distribution. I propose a procedure in which experimenters collect results from M animals and accept results that are observed in at least N of them (‘N-out-of-M’ procedure). I show how to compute the risks α (of reaching an incorrect conclusion) and β (of failing to reach a conclusion) for any prior distribution, and as a function of N and M. Strikingly, I find that the N-out-of-M model leads to a similar conclusion across a wide range of prior distributions: recordings from two animals lowers the risk α and therefore ensures reliable result, but leaves a large risk β; and recordings from three animals and accepting results observed in two of them strikes an efficient balance between acceptable risks α and β. This framework gives a formal justification for the rule of thumb of using at least two animals in NHP studies, suggests that recording from three animals when possible markedly improves statistical power, provides a statistical solution for situations where results are not consistent between all animals, and may apply to other types of studies involving few animals.
Multiple resistance and pH adaptation (Mrp) cation/proton antiporters are essential for growth of a variety of halophilic and alkaliphilic bacteria under stress conditions. Mrp-type antiporters are closely related to the membrane domain of respiratory complex I. We determined the structure of the Mrp antiporter from Bacillus pseudofirmus by electron cryo-microscopy at 2.2 Å resolution. The structure resolves more than 99% of the sidechains of the seven membrane subunits MrpA to MrpG plus 360 water molecules, including ∼70 in putative ion translocation pathways. Molecular dynamics simulations based on the high-resolution structure revealed details of the antiport mechanism. We find that switching the position of a histidine residue between three hydrated pathways in the MrpA subunit is critical for proton transfer that drives gated transmembrane sodium translocation. Several lines of evidence indicate that the same histidine-switch mechanism operates in respiratory complex I.
Our lives (and deaths) have by now been dominated for two years by COVID-19, a pandemic that has caused hundreds of millions of disease cases, millions of deaths, trillions in economic costs, and major restrictions on our freedom. Here we suggest a novel tool for controlling the COVID-19 pandemic. The key element is a method for a population-scale PCR-based testing, applied on a systematic and repeated basis. For this we have developed a low cost, highly sensitive virus-genome-based test. Using Germany as an example, we demonstrate by using a mathematical model, how useful this strategy could have been in controlling the pandemic. We show using real-world examples how this might be implemented on a mass scale and discuss the feasibility of this approach.
Mollusca is the second-largest animal phylum with over 100,000 species among eight distinct taxonomic classes. Across 1000 living species in the class Polyplacophora, chitons have a relatively constrained morphology but with some notable deviations. Several genera possess “shell eyes”, true eyes with a lens and retina that are embedded within the dorsal shells, which represent the most recent evolution of animal eyes. The phylogeny of major chiton clades is mostly well established, in a set of superfamily and higher-level taxa supported by various approaches including multiple gene markers, mitogenome-phylogeny and phylotranscritomic approaches as well as morphological studies. However, one critical lineage has remained unclear: Schizochiton was controversially suggested as a potential independent origin of chiton shell eyes. Here, with the draft genome sequencing of Schizochiton incisus (superfamily Schizochitonoidea) plus assembly of transcriptome data from other polyplacophorans, we present phylogenetic reconstructions using both mitochondrial genomes and phylogenomic approaches with multiple methods. Phylogenetic trees from mitogenomic data are inconsistent, reflecting larger scale confounding factors in molluscan mitogenomes. A consistent robust topology was generated with protein coding genes using different models and methods. Our results support Schizochitonoidea is a sister group to other Chitonoidea in Chitonina, in agreement with established classification. This suggests that the earliest origin of shell eyes is in Schizochitonoidea, which were also gained secondarily in other genera in Chitonoidea. Our results have generated a holistic review of the internal relationship within Polyplacophora, and a better understanding on the evolution of Polyplacophora.
In meditation practices that involve focused attention to a specific object, novice practitioners often experience moments of distraction (i.e., mind wandering). Previous studies have investigated the neural correlates of mind wandering during meditation practice through Electroencephalography (EEG) using linear metrics (e.g., oscillatory power). However, their results are not fully consistent. Since the brain is known to be a chaotic/nonlinear system, it is possible that linear metrics cannot fully capture complex dynamics present in the EEG signal. In this study, we assess whether nonlinear EEG signatures can be used to characterize mind wandering during breath focus meditation in novice practitioners. For that purpose, we adopted an experience sampling paradigm in which 25 participants were iteratively interrupted during meditation practice to report whether they were focusing on the breath or thinking about something else. We compared the complexity of EEG signals during mind wandering and breath focus states using three different algorithms: Higuchi’s fractal dimension (HFD), Lempel-Ziv complexity (LZC), and Sample entropy (SampEn). Our results showed that EEG complexity was generally reduced during mind wandering relative to breath focus states. We conclude that EEG complexity metrics are appropriate to disentangle mind wandering from breath focus states in novice meditation practitioners, and therefore, they could be used in future EEG neurofeedback protocols to facilitate meditation practice.
When speech is too fast, the tracking of the acoustic signal along the auditory pathway deteriorates, leading to suboptimal speech segmentation and decoding of speech information. Thus, speech comprehension is limited by the temporal constraints of the auditory system. Here we ask whether individual differences in auditory-motor coupling strength in part shape these temporal constraints. In two behavioral experiments, we characterize individual differences in the comprehension of naturalistic speech as function of the individual synchronization between the auditory and motor systems and the preferred frequencies of the systems. Obviously, speech comprehension declined at higher speech rates. Importantly, however, both higher auditory-motor synchronization and higher spontaneous speech motor production rates were predictive of better speech-comprehension performance. Furthermore, performance increased with higher working memory capacity (Digit Span) and higher linguistic, model-based sentence predictability – particularly so at higher speech rates and for individuals with high auditory-motor synchronization. These findings support the notion of an individual preferred auditory– motor regime that allows for optimal speech processing. The data provide evidence for a model that assigns a central role to motor-system-dependent individual flexibility in continuous speech comprehension.
In this work, we study for the first time the thermal evolution of twin star pairs, i.e., stars that present the same mass but different radius and compactness. We collect available equations of state that give origin to a second branch of stable compact stars with quarks in their core. For each equation of state, we investigate the particle composition inside stars and how differently each twin evolves over time, which depends on the central density/pressure and consequent crossing of the threshold for the Urca cooling process. We find that, although the general stellar thermal evolution depends on mass and particle composition, withing one equation of state, only twin pairs that differ considerably on compactness can be clearly distinguished by how they cool down.
Although new advances in neuroscience allow the study of vocal communication in awake animals, substantial progress in the processing of vocalizations has been made from brains of anaesthetized preparations. Thus, understanding how anaesthetics affect neuronal responses is of paramount importance. Here, we used electrophysiological recordings and computational modelling to study how the auditory cortex of bats responds to vocalizations under anaesthesia and in wakefulness. We found that multifunctional neurons that process echolocation and communication sounds were affected by ketamine anaesthesia in a manner that could not be predicted by known anaesthetic effects. In wakefulness, acoustic contexts (preceding echolocation or communication sequences) led to stimulus-specific suppression of lagging sounds, accentuating neuronal responses to sound transitions. However, under anaesthesia, communication contexts (but not echolocation) led to a global suppression of responses to lagging sounds. Such asymmetric effect was dependent on the frequency composition of the contexts and not on their temporal patterns. We constructed a neuron model that could replicate the data obtained in vivo. In the model, anaesthesia modulates spiking activity in a channel-specific manner, decreasing responses of cortical inputs tuned to high-frequency sounds and increasing adaptation in the respective cortical synapses. Combined, our findings obtained in vivo and in silico reveal that ketamine anaesthesia does not reduce uniformly the neurons’ responsiveness to low and high frequency sounds. This effect depends on combined mechanisms that unbalance cortical inputs and ultimately affect how auditory cortex neurons respond to natural sounds in anaesthetized preparations.
A biophysical regulator of inhibitory integration and learning in mesolimbic dopamine neurons
(2022)
Midbrain dopamine neurons are essential for flexible control of adaptive behaviors. DA neurons that project to different target regions have unique biophysical properties, and it is thought that this diversity reflects functional specialization. This assumption implies the presence of specific genetic determinants with precise impacts on behavior. We tested this general hypothesis by homing in on one particular biophysical mechanism, Kv4 channel inactivation, using a combination of molecular, proteomic, electrophysiological, computational, and behavioral approaches. We demonstrate that KChIP4a, a singular Kv4 β-subunit splice variant, prolongs hyperpolarization-rebound delays selectively in dopamine neurons projecting to the nucleus accumbens core, shifts the integration of inhibitory inputs and, in turn, selectively regulates learning from negative prediction-errors. Our results reveal a highly specialized, gene-to-behavior mechanistic chain that is only operative in a particular dopaminergic subsystem, illuminating how molecularly defined biophysical switches are employed for neuron subtype-specific information processing in the brain.
The knowledge that brain functional connectomes are both unique and reliable has enabled behaviourally relevant inferences at a subject level. However, it is unknown whether such “fingerprints” persist under altered states of consciousness. Ayahuasca is a potent serotonergic psychedelic which elicits a widespread dysregulation of functional connectivity. Used communally in religious ceremonies, its shared use may highlight relevant novel interactions between mental state and FC inherency. Using 7T fMRI, we assessed resting-state static and dynamic FCs for 21 Santo Daime members after collective ayahuasca intake in an acute, within-subject study. Here, connectome fingerprinting revealed a shared functional space, accompanied by a spatiotemporal reallocation of keypoint edges. Importantly, we show that interindividual differences in higher-order FCs motifs are relevant to experiential phenotypes, given that they can predict perceptual drug effects. Collectively, our findings offer an example as to how individualised connectivity markers can be used to trace a subject’s functional connectome across altered states of consciousness.
In this article we provide a stack-theoretic framework to study the universal tropical Jacobian over the moduli space of tropical curves. We develop two approaches to the process of tropicalization of the universal compactified Jacobian over the moduli space of curves -- one from a logarithmic and the other from a non-Archimedean analytic point of view. The central result from both points of view is that the tropicalization of the universal compactified Jacobian is the universal tropical Jacobian and that the tropicalization maps in each of the two contexts are compatible with the tautological morphisms. In a sequel we will use the techniques developed here to provide explicit polyhedral models for the logarithmic Picard variety.
Intraspecific genomic variability affects a species’ adaptive potential towards climatic conditions. Variation in gene content across populations and environments may point at genomic adaptations to specific environments. The lichen symbiosis, a stable association of fungal and photobiont partners, offers an excellent system to study environmentally driven gene content variation. Many species have remarkable environmental tolerances, and often form populations in different climate zones. Here we combine comparative and population genomics to assess the presence and absence of genes in high elevation and low elevation genomes of two lichenized fungi of the genus Umbilicaria. The two species have non-overlapping ranges, but occupy similar climatic niches in North America (U. phaea) and Europe (U. pustulata): high elevation populations are located in the cold temperate zone and low elevation populations in the Mediterranean zone. We assessed gene content variation along replicated elevation gradients in each of the two species, based on a total of 2050 individuals across 26 populations. Specifically, we assessed shared orthologs across species within the same climate zone, and tracked which genes increase or decrease in abundance within populations along elevation. In total, we found 16 orthogroups with shared orthologous genes in genomes at low elevation and 13 at high elevation. Coverage analysis revealed one ortholog that is exclusive to genomes at low elevation. Conserved domain search revealed domains common to the protein kinases (PKs) superfamily. We traced the discovered ortholog in populations along five replicated elevation gradients on both continents. The protein kinase gene linearly declined in abundance with increasing elevation, and was absent in the highest populations. We consider the parallel loss of an ortholog in two species and in two geographic settings a rare find, and a step forward in understanding the genomic underpinnings of climatic tolerances in lichenized fungi. In addition, the tracking of gene content variation provides a widely applicable framework for retrieving biogeographical determinants of gene presence/absence patterns. Our work provides insights into gene content variation of lichenized fungi in relation to climatic gradients, suggesting a new research direction with implications for understanding evolutionary trajectories of complex symbioses in relation to climatic change.
Moving in synchrony to external rhythmic stimuli is an elementary function that humans regularly engage in. It is termed “sensorimotor synchronization” and it is governed by two main parameters, the period and the phase of the movement with respect to the external rhythm. There has been an extensive body of research on the characteristics of these parameters, primarily once the movement synchronization has reached a steady-state level. Particular interest has been shown about how these parameters are corrected when there are deviations for the steady-state level. However, little is known about the initial “tuning-in” interval, when one aligns the movement to the external rhythm from rest. The current work investigates this “tuning-in” period for each of the four limbs and makes various novel contributions in the understanding of sensorimotor synchronization. The results suggest that phase and period alignment appear to be separate processes. Phase alignment involves limb-specific somatosensory memory in the order of minutes while period alignment has very limited memory usage. Phase alignment is the primary task but then the brain switches to period alignment where it spends most its resources. In overall this work suggests a central, cognitive role of period alignment and a peripheral, sensorimotor role of phase alignment.
Research on psychopathy has so far been largely limited to the investigation of high-level processes, such as emotion perception and regulation. In the present work, we investigate whether psychopathy has an effect on the estimation of fundamental physical parameters, which are computed in the brain during early stages of sensory processing. We employed a simple task in which participants had to estimate their interpersonal distance from a moving avatar and stop it at a given distance. The face expression of the avatars were positive, negative, or neutral. Participants carried out the task online on their home computers. We measured the psychopathy level via a self-report questionnaire. Regardless of the degree of psychopathy, the facial expression of the avatars showed no effect on distance estimation. Our results show that individuals with a high degree of psychopathy underestimate distance of approaching avatars significantly less (let the avatar approach them significantly closer) than did participants with a lesser degree of psychopathy. Moreover, participants who scored high in Self-Centered Impulsivity underestimate the distance to approaching avatars significantly less (let the avatar approach closer) than participants with a low score. Distance estimation is considered an automatic process performed at early stages of visual processing. Therefore, our results imply that psychopathy affects basic early sensory processes, such as feature extraction, in the visual cortex.
The neural mechanisms that unfold when humans form a large group defined by an overarching context, such as audiences in theater or sports, are largely unknown and unexplored. This is mainly due to the lack of availability of a scalable system that can record the brain activity from a significantly large portion of such an audience simultaneously. Although the technology for such a system has been readily available for a long time, the high cost as well as the large overhead in human resources and logistic planning have prohibited the development of such a system. However, during the recent years reduction in technology costs and size have led to the emergence of low-cost, consumer-oriented EEG systems, developed primarily for recreational use. Here by combining such a low-cost EEG system with other off-the-shelve hardware and tailor-made software, we develop in the lab and test in a cinema such a scalable EEG hyper-scanning system. The system has a robust and stable performance and achieves accurate unambiguous alignment of the recorded data of the different EEG headsets. These characteristics combined with small preparation time and low-cost make it an ideal candidate for recording large portions of audiences.
Neuronal hyperexcitability is a feature of Alzheimer’s disease (AD). Three main mechanisms have been proposed to explain it: i), dendritic degeneration leading to increased input resistance, ii), ion channel changes leading to enhanced intrinsic excitability, and iii), synaptic changes leading to excitation-inhibition (E/I) imbalance. However, the relative contribution of these mechanisms is not fully understood. Therefore, we performed biophysically realistic multi-compartmental modelling of excitability in reconstructed CA1 pyramidal neurons of wild-type and APP/PS1 mice, a well-established animal model of AD. We show that, for synaptic activation, the excitability promoting effects of dendritic degeneration are cancelled out by excitability decreasing effects of synaptic loss. We find an interesting balance of excitability regulation with enhanced degeneration in the basal dendrites of APP/PS1 cells potentially leading to increased excitation by the apical but decreased excitation by the basal Schaffer collateral pathway. Furthermore, our simulations reveal that three additional pathomechanistic scenarios can account for the experimentally observed increase in firing and bursting of CA1 pyramidal neurons in APP/PS1 mice. Scenario 1: increased excitatory burst input; scenario 2: enhanced E/I ratio and scenario 3: alteration of intrinsic ion channels (IAHP down-regulated; INap, INa and ICaT up-regulated) in addition to enhanced E/I ratio. Our work supports the hypothesis that pathological network and ion channel changes are major contributors to neuronal hyperexcitability in AD. Overall, our results are in line with the concept of multi-causality and degeneracy according to which multiple different disruptions are separately sufficient but no single disruption is necessary for neuronal hyperexcitability.
The ability to extract regularities from the environment is arguably an adaptive characteristic of intelligent systems. In the context of speech, statistical learning is thought to be an important mechanism for language acquisition. By considering individual differences in speech auditory-motor synchronization, an independent component analysis of fMRI data revealed that the neural substrates of statistical word form learning are not fully shared across individuals. While a network of auditory and superior pre/motor regions is universally activated in the process of learning, a fronto-parietal network is instead additionally and selectively engaged by some individuals, boosting their performance. Furthermore, interfering with the use of this network via articulatory suppression (producing irrelevant speech during learning) normalizes performance across the entire sample. Our work provides novel insights on language-related statistical learning and reconciles previous contrasting findings, while highlighting the need to factor in fundamental individual differences for a precise characterization of cognitive phenomena.
Research points to neurofunctional differences underlying fluent speech production in stutterers and non-stutterers. There has been considerably less work focusing on the processes that underlie stuttered speech, primarily due to the difficulty of reliably eliciting stuttering in the unnatural contexts associated with neuroimaging experiments. We used magnetoencephalography (MEG) to test the hypothesis that stuttering events result from global motor inhibition–a “freeze” response typically characterized by increased beta power in nodes of the action-stopping network. We leveraged a novel clinical interview to develop participant-specific stimuli in order to elicit a comparable amount of stuttered and fluent trials. Twenty-nine adult stutterers participated. The paradigm included a cue prior to a go signal, which allowed us to isolate processes associated with stuttered and fluent trials prior to speech initiation. During this pre-speech time window, stuttered trials were associated with greater beta power in the right pre-supplementary motor area, a key node in the action-stopping network, compared to fluent trials. Beta power in the right pre-supplementary area was related to a clinical measure of stuttering severity. We also found that anticipated words identified independently by participants were stuttered more often than those generated by the researchers, which were based on the participants’ reported anticipated sounds. This suggests that global motor inhibition results from stuttering anticipation. This study represents the largest comparison of stuttered and fluent speech to date. The findings provide a foundation for clinical trials that test the efficacy of neuromodulation on stuttering. Moreover, our study demonstrates the feasibility of using our approach for eliciting stuttering during MEG and functional magnetic resonance imaging experiments so that the neurobiological bases of stuttered speech can be further elucidated.
Speech imagery (the ability to generate internally quasi-perceptual experiences of speech) is a fundamental ability linked to cognitive functions such as inner speech, phonological working memory, and predictive processing. Speech imagery is also considered an ideal tool to test theories of overt speech. The study of speech imagery is challenging, primarily because of the absence of overt behavioral output as well as the difficulty in temporally aligning imagery events across trials and individuals. We used magnetoencephalography (MEG) paired with temporal-generalization-based neural decoding and a simple behavioral protocol to determine the processing stages underlying speech imagery. We monitored participants’ lip and jaw micromovements during mental imagery of syllable production using electromyography. Decoding participants’ imagined syllables revealed a sequence of task-elicited representations. Importantly, participants’ micromovements did not discriminate between syllables. The decoded sequence of neuronal patterns maps well onto the predictions of current computational models of overt speech motor control and provides evidence for hypothesized internal and external feedback loops for speech planning and production, respectively. Additionally, the results expose the compressed nature of representations during planning which contrasts with the natural rate at which internal productions unfold. We conjecture that the same sequence underlies the motor-based generation of sensory predictions that modulate speech perception as well as the hypothesized articulatory loop of phonological working memory. The results underscore the potential of speech imagery, based on new experimental approaches and analytical methods, and further pave the way for successful non-invasive brain-computer interfaces.
We study the μ-μ45-T phase diagram of the 2+1-dimensional Gross-Neveu model, where μ denotes the ordinary chemical potential, μ45 the chiral chemical potential and T the temperature. We use the mean-field approximation and two different lattice regularizations with naive chiral fermions. An inhomogeneous phase at finite lattice spacing is found for one of the two regularizations. Our results suggest that there is no inhomogeneous phase in the continuum limit. We show that a chiral chemical potential is equivalent to an isospin chemical potential. Thus, all results presented in this work can also be interpreted in the context of isospin imbalance.
A prototype of a new type of calorimeter has been designed and constructed, based on a silicon-tungsten sampling design using pixel sensors with digital readout. It makes use of the Alpide MAPS sensor developed for the ALICE ITS upgrade. A binary readout is possible due to the pixel size of ≈30×30μm2. This prototype has been successfully tested with cosmic muons and with test beams at DESY and the CERN SPS. We report on performance results obtained at DESY, showing good energy resolution and linearity, and compare to detailed MC simulations. Also shown are preliminary results of the high-energy performance as measured at the SPS. The two-shower separation capabilities are discussed.
Brookshire (2022) claims that previous analyses of periodicity in detection performance after a reset event suffer from extreme false-positive rates. Here we show that this conclusion is based on an incorrect implemention of a null-hypothesis of aperiodicity, and that a correct implementation confirms low false-positive rates. Furthermore, we clarify that the previously used method of shuffling-in-time, and thereby shuffling-in-phase, cleanly implements the null hypothesis of no temporal structure after the reset, and thereby of no phase locking to the reset. Moving from a corresponding phase-locking spectrum to an inference on the periodicity of the underlying process can be accomplished by parameterizing the spectrum. This can separate periodic from non-periodic components, and quantify the strength of periodicity.
Coronavirus disease 2019 (COVID-19) is caused by the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) and can affect multiple organs, among which is the circulatory system. Inflammation and mortality risk markers were previously detected in COVID-19 plasma and red blood cells (RBCs) metabolic and proteomic profiles. Additionally, biophysical properties, such as deformability, were found to be changed during the infection. Based on such data, we aim to better characterize RBC functions in COVID-19. We evaluate the flow properties of RBCs in severe COVID-19 patients admitted to the intensive care unit by using in vitro microfluidic techniques and automated methods, including artificial neural networks, for an unbiased RBC analysis. We find strong flow and RBC shape impairment in COVID-19 samples and demonstrate that such changes are reversible upon suspension of COVID-19 RBCs in healthy plasma. Vice versa, healthy RBCs immediately resemble COVID-19 RBCs when suspended in COVID-19 plasma. Proteomics and metabolomics analyses allow us to detect the effect of plasma exchanges on both plasma and RBCs and demonstrate a new role of RBCs in maintaining plasma equilibria at the expense of their flow properties. Our findings provide a framework for further investigations of clinical relevance for therapies against COVID-19 and possibly other infectious diseases.
Coronavirus disease 2019 (COVID-19) is caused by the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) and can affect multiple organs, among which is the circulatory system. Inflammation and mortality risk markers were previously detected in COVID-19 plasma and red blood cells (RBCs) metabolic and proteomic profiles. Additionally, biophysical properties, such as deformability, were found to be changed during the infection. Based on such data, we aim to better characterize RBC functions in COVID-19. We evaluate the flow properties of RBCs in severe COVID-19 patients admitted to the intensive care unit by using in vitro microfluidic techniques and automated methods, including artificial neural networks, for an unbiased RBC analysis. We find strong flow and RBC shape impairment in COVID-19 samples and demonstrate that such changes are reversible upon suspension of COVID-19 RBCs in healthy plasma. Vice versa, healthy RBCs immediately resemble COVID-19 RBCs when suspended in COVID-19 plasma. Proteomics and metabolomics analyses allow us to detect the effect of plasma exchanges on both plasma and RBCs and demonstrate a new role of RBCs in maintaining plasma equilibria at the expense of their flow properties. Our findings provide a framework for further investigations of clinical relevance for therapies against COVID-19 and possibly other infectious diseases.
After myocardial infarction in the adult heart the remaining, non-infarcted tissue adapts to compensate the loss of functional tissue. This adaptation requires changes in gene expression networks, which are mostly controlled by transcription regulating proteins. Long non-coding transcripts (lncRNAs) are now recognized for taking part in fine-tuning such gene programs. We identified and characterized the cardiomyocyte specific lncRNA Sweetheart RNA (Swhtr), an approximately 10 kb long transcript divergently expressed from the cardiac core transcription factor coding gene Nkx2-5. We show that Swhtr is dispensable for normal heart development and function, but becomes essential for the tissue adaptation process after myocardial infarction. Re-expressing Swhtr from an exogenous locus rescues the Swhtr null phenotype. Genes depending on Swhtr after cardiac stress are significantly occupied, and therefore most likely regulated by NKX2-5. Our results indicate a synergistic role for Swhtr and the developmentally essential transcription factor NKX2-5 in tissue adaptation after myocardial injury.
The discovery of superconductivity in layered vanadium-based kagome metals AV3Sb5 (A: K, Rb, Cs) has added a new family of materials to the growing class of possible unconventional superconductors. However, the nature of the superconducting pairing in these materials remains elusive. We present a microscopic theoretical study of the leading superconducting instabilities on the kagome lattice based on spin- and charge-fluctuation mediated Cooper pairing. The applied methodology includes effects of both on-site and nearest-neighbor repulsive Coulomb interactions. Near the upper van Hove filling -- relevant for the AV3Sb5 materials -- we find a rich phase diagram with several pairing symmetries being nearly degenerate. In particular, while a substantial fraction of the phase diagram is occupied by a spin-singlet order parameter transforming as a two-dimensional irreducible representation of the point group, several nodal spin-triplet pairing states remain competitive. We compute the band and interaction parameter-dependence of the hierarchy of the leading superconducting instabilities, and determine the detailed momentum dependence of the resulting preferred gap structures. Crucially, for moderate values of the interaction parameters, the individual pairing states depend strongly on momentum and exhibit multiple nodes on the Fermi surface. We discuss the properties of these superconducting gap structures in light of recent experimental developments of the AV3Sb5 materials.
The electrical and computational properties of neurons in our brains are determined by a rich repertoire of membrane-spanning ion channels and elaborate dendritic trees. However, the precise reason for this inherent complexity remains unknown. Here, we generated large stochastic populations of biophysically realistic hippocampal granule cell models comparing those with all 15 ion channels to their reduced but functional counterparts containing only 5 ion channels. Strikingly, valid parameter combinations in the full models were more frequent and more stable in the face of perturbations to channel expression levels. Scaling up the numbers of ion channels artificially in the reduced models recovered these advantages confirming the key contribution of the actual number of ion channel types. We conclude that the diversity of ion channels gives a neuron greater flexibility and robustness to achieve target excitability.
From early to middle childhood, brain regions that underlie memory consolidation undergo profound maturational changes. However, there is little empirical investigation that directly relates age-related differences in brain structural measures to the memory consolidation processes. The present study examined system-level memory consolidations of intentionally studied object-location associations after one night of sleep (short delay) and after two weeks (long delay) in normally developing 5-to-7-year-old children (n = 50) and young adults (n = 39). Behavioural differences in memory consolidation were related to structural brain measures. Our results showed that children, in comparison to young adults, consolidate correctly learnt object-location associations less robustly over short and long delay. Moreover, using partial least squares correlation method, a unique multivariate profile comprised of specific neocortical (prefrontal, parietal, and occipital), cerebellar, and hippocampal subfield structures was found to be associated with variation in short-delay memory consolidation. A different multivariate profile comprised of a reduced set of brain structures, mainly consisting of neocortical (prefrontal, parietal, and occipital), and selective hippocampal subfield structures (CA1-2 and subiculum) was associated with variation in long-delay memory consolidation. Taken together, the results suggest that multivariate structural pattern of unique sets of brain regions are related to variations in short- and long-delay memory consolidation across children and young adults.
RESEARCH HIGHLIGHTS
* Short- and long-delay memory consolidation is less robust in children than in young adults
* Short-delay brain profile comprised of hippocampal, cerebellar, and neocortical brain regions
* Long-delay brain profile comprised of neocortical and selected hippocampal brain regions.
* Brain profiles differ between children and young adults.
Hydrogen is a promising fuel in a carbon-neutral economy, and many efforts are currently undertaken to produce hydrogen. One of the challenges is to store and transport the highly explosive gas in a safe and easy way. One option that is intensively analyzed by chemists and biologists is the conversion of hydrogen and CO2 to formic acid, the liquid organic hydrogen carrier. Here, we demonstrate for the first time that a bio-based system, using Acetobacterium woodii as the biocatalyst, allows multiple cycles of bi-directional hydrogenation of CO2 to formic acid in one bioreactor. The process was kept running over 2 weeks producing and oxidizing 330 mM formic acid in total. Unwanted side-product formation of acetic acid was prevented through metabolic engineering of the organism. The demonstrated process design can be considered as a future “bio-battery” for the reversible storage of electrons in the form of H2 in formic acid, a versatile compound.
Gasdermin-D (GSDMD) is the ultimate effector of pyroptosis, a form of programmed cell death associated with pathogen invasion and inflammation. After proteolytic cleavage by caspases activated by the inflammasome, the GSDMD N-terminal domain (GSDMDNT) assembles on the inner leaflet of the plasma membrane and induces the formation of large membrane pores. We use atomistic molecular dynamics simulations to study GSDMDNT monomers, oligomers, and rings in an asymmetric plasma membrane mimetic. We identify distinct interaction motifs of GSDMDNT with phosphatidylinositol-4,5-bisphosphate (PI(4,5)P2) and phosphatidylserine (PS) head-groups and describe differential lipid binding between the pore and prepore conformations. Oligomers are stabilized by shared lipid binding sites between neighboring monomers acting akin to double-sided tape. We show that already small GSDMDNT oligomers form stable, water-filled and ion-conducting membrane pores bounded by curled beta-sheets. In large-scale simulations, we resolve the process of pore formation by lipid detachment from GSDMDNT arcs and lipid efflux from partial rings. We find that that high-order GSDMDNT oligomers can crack under the line tension of 86 pN created by an open membrane edge to form the slit pores or closed GSDMDNT rings seen in experiment. Our simulations provide a detailed view of key steps in GSDMDNT-induced plasma membrane pore formation, including sublytic pores that explain nonselective ion flux during early pyroptosis.
The heterotetrameric human transfer RNA (tRNA) splicing endonuclease (TSEN) catalyzes the excision of intronic sequences from precursor tRNAs (pre-tRNAs)1. Mutations in TSEN and its associated RNA kinase CLP1 are linked to the neurodegenerative disease pontocerebellar hypoplasia (PCH)2–8. The three-dimensional (3D) assembly of TSEN/CLP1, the mechanism of substrate recognition, and the molecular details of PCH-associated mutations are not fully understood. Here, we present cryo-electron microscopy structures of human TSEN with intron-containing pre-tRNATyrgta and pre-tRNAArgtct. TSEN exhibits broad structural homology to archaeal endonucleases9 but has evolved additional regulatory elements that are involved in handling and positioning substrate RNA. Essential catalytic residues of subunit TSEN34 are organized for the 3’ splice site which emerges from a bulge-helix configuration. The triple-nucleotide bulge at the intron/3’-exon boundary is stabilized by an arginine tweezer motif of TSEN2 and an interaction with the proximal minor groove of the helix. TSEN34 and TSEN54 define the 3’ splice site by holding the tRNA body in place. TSEN54 adapts a bipartite fold with a flexible central region required for CLP1 binding. PCH-associated mutations are located far from pre-tRNA binding interfaces explaining their negative impact on structural integrity of TSEN without abrogating its catalytic activity in vitro10. Our work defines the molecular framework of pre-tRNA recognition and cleavage by TSEN and provides a structural basis to better understand PCH in the future.
The new variant of concern (VOC) of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), Omicron (B.1.1.529), is genetically very different from other VOCs. We compared Omicron with the preceding VOC Delta (B.1.617.2) and the wildtype strain (B.1) with respect to their interactions with the antiviral type I interferon (IFN-alpha/beta) response in infected cells. Our data indicate that Omicron has gained an elevated capability to suppress IFN-beta induction upon infection and to better withstand the antiviral state imposed by exogenously added IFN-alpha.
Precisely estimating event timing is essential for survival, yet temporal distortions are ubiquitous in our daily sensory experience. Here, we tested whether the relative position, relative duration and relative distance in time of two sequentially-organized events —standard S, with constant duration, and comparison C, varying trial-by-trial— are causal factors in generating temporal distortions. We found that temporal distortions emerge when the first event is shorter than the second event. Importantly, a significant interaction suggests that a longer ISI helps counteracting such serial distortion effect only the constant S is in first position, but not if the unpredictable C is in first position. These results suggest the existence of a perceptual bias in perceiving ordered event durations, mechanistically contributing to distortion in time perception. We simulated our behavioral results with a Bayesian model and replicated the finding that participants disproportionately expand first-position dynamic (unpredictable) short events. Our results clarify the mechanics generating time distortions by identifying a hitherto unknown duration-dependent encoding inefficiency in human serial temporal perception, akin to a strong prior that can be overridden for highly predictable sensory events but unfolds for unpredictable ones.
TriMem: a parallelized hybrid Monte Carlo software for efficient simulations of lipid membranes
(2022)
Lipid membranes are integral building blocks of living cells and perform a multitude of biological functions. Currently, molecular simulations of cellular-scale membrane structures at atomic resolution are nearly impossible, due to their size, complexity, and the large times-scales required. Instead, elastic membrane models are used to simulate membrane topologies and transitions between them, and to infer their properties and functions. Unfortunately, efficiently parallelized open-source simulation code to do so has been lacking. Here, we present TriMem, a parallel hybrid Monte Carlo simulation engine for triangulated lipid membranes. The kernels are efficiently coded in C++ and wrapped with Python for ease-of-use. The parallel implementation of the energy and gradient calculations and of Monte Carlo flip moves of edges in the triangulated membrane enable us to simulate also large and highly curved sub-cellular structures. For validation, we reproduce phase diagrams of vesicles with varying surface-to-volume ratios and area difference. The software can tackle a range of membrane remodelling processes on sub-cellular and cellular scales. Additionally, extensive documentation make the software accessible to the broad biophysics and computational cell biology communities.
A candidate gene cluster for the bioactive natural product gyrophoric acid in lichen-forming fungi
(2022)
Natural products of lichen-forming fungi are structurally diverse and have a variety of medicinal properties. Despite this, they a have limited implementation in industry, because the corresponding genes remain unknown for most of the natural products. Here we implement a long-read sequencing and bioinformatic approach to identify the biosynthetic gene cluster of the bioactive natural product gyrophoric acid (GA). Using 15 high-quality genomes representing nine GA-producing species of the lichen-forming fungal genus Umbilicaria, we identify the most likely GA cluster and investigate cluster gene organization and composition across the nine species. Our results show that GA clusters are promiscuous within Umbilicaria, with only three genes that are conserved across species, including the PKS gene. In addition, our results suggest that the same cluster codes for different but structurally similar NPs, i.e., GA, umbilicaric acid and hiascic acid, bringing new evidence that lichen metabolite diversity is also generated through regulatory mechanisms at the molecular level. Ours is the first study to identify the most likely GA cluster, and thus provides essential information to open new avenues for biotechnological approaches to producing and modifying GA and similar lichen-derived compounds. We show that bioinformatics approaches are useful in linking genes and potentially associated natural products. Genome analyses help unlocking the pharmaceutical potential of organisms such as lichens, which are biosynthetically diverse but slow growing, and difficult to cultivate due to their symbiotic nature.
Genome mining as a biotechnological tool for the discovery of novel biosynthetic genes in lichens
(2022)
The ever-increasing demand for novel drugs highlights the need for bioprospecting unexplored taxa for their biosynthetic potential. Lichen-forming fungi (LFF) are a rich source of natural products but their implementation in pharmaceutical industry is limited, mostly because the genes corresponding to a majority of their natural products is unknown. Furthermore, it is not known to what extent these genes encode structurally novel molecules. Advance in next-generation sequencing technologies has expanded the range of organisms that could be exploited for their biosynthetic potential. In this study, we mine the genomes of nine lichen-forming fungal species of the genus Umbilicaria for biosynthetic genes, and categorize the BGCs as “associated product structurally known”, and “associated product putatively novel”. We found that about 25-30% of the biosynthetic genes are divergent when compared to the global database of BGCs comprising of 1,200,000 characterized biosynthetic genes from planta, bacteria and fungi. Out of 217 total BGCs, 43 were only distantly related to known BGCs, suggesting they encode structurally and functionally unknown natural products. Clusters encoding the putatively novel metabolic diversity comprise PKSs (30), NRPSs (12) and terpenes (1). Our study emphasizes the utility of genomic data in bioprospecting microorganisms for their biosynthetic potential and in advancing the industrial application of unexplored taxa. We highlight the untapped structural metabolic diversity encoded in the lichenized fungal genomes. To the best of our knowledge, this is the first investigation identifying genes coding for NPs with potentially novel therapeutic properties in LFF.
Under temperature or pressure tuning, tetragonal EuPd2Si2 is known to undergo a valence transition from nearly divalent to nearly trivalent Eu accompanied by a volume reduction. Albeit intensive work, its origin is not yet completely understood. Here, we investigate the mechanism of the valence transition under volume compression by density functional theory calculations (DFT). Our analysis suggests that the transition is a consequence of an enhanced c−f hybridization between localized Eu 4f states and itinerant conduction states (Eu 5d, Pd 4d, and Si 3p) where the interplay of the electronic bandwidth, crystal field environment, Coulomb repulsion, Hund's coupling and spin-orbit coupling plays a crucial role for the transition to happen. The change in the electronic structure is intimately related to the volume reduction where Eu-Pd(Si) bond lengths shorten. In a next step we compare our DFT results to surface-sensitive photoemission data in which the mixed-valent properties of EuPd2Si2 are reflected in a simultaneous observation of divalent and trivalent signals from the Eu 4f shell.
Spatial attention increases both inter-areal synchronization and spike rates across the visual hierarchy. To investigate whether these attentional changes reflect distinct or common mechanisms, we performed simultaneous laminar recordings of identified cell classes in macaque V1 and V4. Enhanced V4 spike rates were expressed by both excitatory neurons and fast-spiking interneurons, and were most prominent and arose earliest in time in superficial layers, consistent with a feedback modulation. By contrast, V1-V4 gamma-synchronization reflected feedforward communication and surprisingly engaged only fast-spiking interneurons in the V4 input layer. In mouse visual cortex, we found a similar motif for optogenetically identified inhibitory-interneuron classes. Population decoding analyses further indicate that feedback-related increases in spikes rates encoded attention more reliably than feedforward-related increases in synchronization. These findings reveal distinct, cell-type-specific feedforward and feedback pathways for the attentional modulation of inter-areal synchronization and spike rates, respectively.