Refine
Year of publication
Document Type
- Preprint (2316) (remove)
Has Fulltext
- yes (2316)
Keywords
- Kollisionen schwerer Ionen (33)
- heavy ion collisions (27)
- Deutsch (23)
- Quark-Gluon-Plasma (14)
- equation of state (13)
- QGP (12)
- heavy-ion collisions (11)
- Kongress (10)
- Syntax (10)
- quark-gluon plasma (10)
Institute
- Physik (1424)
- Frankfurt Institute for Advanced Studies (FIAS) (1006)
- Informatik (803)
- Medizin (176)
- Extern (82)
- Biowissenschaften (76)
- Ernst Strüngmann Institut (70)
- Mathematik (48)
- MPI für Hirnforschung (47)
- Psychologie (46)
RNA-sequencing analyses are often limited to identifying lowest p-value transcripts, which does not address polygenic phenomena. To overcome this limitation, we developed an integrative approach that combines large scale transcriptomic meta-analysis of patient brain tissues with single-cell sequencing data of CNS neurons, short RNA-sequencing of human male- and female-originated cell lines, and connectomics of transcription factor- and microRNA-interactions with perturbed transcripts. We used this pipeline to analyze cortical transcripts of schizophrenia and bipolar disorder patients. While these pathologies show massive transcriptional parallels, their clinically well-known sexual dimorphisms remain unexplained. Our method explicates the differences between afflicted men and women, and identifies disease-affected pathways of cholinergic transmission and gp130-family neurokine controllers of immune function, interlinked by microRNAs. This approach may open new perspectives for seeking biomarkers and therapeutic targets, also in other transmitter systems and diseases.
The SLC26 family of transporters maintains anion equilibria in all kingdoms of life. The family shares a 7 + 7 transmembrane segments inverted repeat architecture with the SLC4 and SLC23 families, but holds a regulatory STAS domain in addition. While the only experimental SLC26 structure is monomeric, SLC26 proteins form structural and functional dimers in the lipid membrane. Here we resolve the structure of an SLC26 dimer embedded in a lipid membrane and characterize its functional relevance by combining PELDOR distance measurements and biochemical studies with MD simulations and spin-label ensemble refinement. Our structural model reveals a unique interface different from the SLC4 and SLC23 families. The functionally relevant STAS domain exerts a stabilizing effect on regions central in this dimer. Characterization of heterodimers indicates that protomers in the dimer functionally interact. The combined structural and functional data define the framework for a mechanistic understanding of functional cooperativity in SLC26 dimers.
Autophagy is a highly conserved catabolic process through which defective or otherwise harmful cellular components are targeted for degradation via the lysosomal route. Regulatory pathways, involving post-translational modifications such as phosphorylation, play a critical role in controlling this tightly orchestrated process. Here, we demonstrate that TBK1 regulates autophagy by phosphorylating autophagy modifiers LC3C and GABARAP-L2 on surface-exposed serine residues (LC3C S93 and S96; GABARAP-L2 S87 and S88). This phosphorylation event impedes their binding to the processing enzyme ATG4 by destabilizing the complex. Phosphorylated LC3C/GABARAP-L2 cannot be removed from liposomes by ATG4 and are thus protected from ATG4-mediated premature removal from nascent autoph-agosomes. This ensures a steady coat of lipidated LC3C/GABARAP-L2 throughout the early steps in autophagosome formation and aids in maintaining a unidirectional flow of the autophagosome to the lysosome. Taken together, we present a new regulatory mechanism of autophagy, which influences the conjugation and de-conjugation of LC3C and GABARAP-L2 to autophagosomes by TBK1-mediated phosphorylation.
Molecular mechanisms of inorganic-phosphate release from the core and barbed end of actin filaments
(2023)
The release of inorganic phosphate (Pi) from actin filaments constitutes a key step in their regulated turnover, which is fundamental to many cellular functions. However, the molecular mechanisms underlying Pi release from both the core and barbed end of actin filaments remain unclear. Here, we combine cryo-EM with molecular dynamics simulations and in vitro reconstitution to demonstrate how actin releases Pi through a ‘molecular backdoor’. While constantly open at the barbed end, the backdoor is predominantly closed in filament-core subunits and only opens transiently through concerted backbone movements and rotameric rearrangements of residues close to the nucleotide binding pocket. This mechanism explains why Pi escapes rapidly from the filament end and yet slowly from internal actin subunits. In an actin variant associated with nemaline myopathy, the backdoor is predominantly open in filament-core subunits, resulting in greatly accelerated Pi release after polymerization and filaments with drastically shortened ADP-Pi caps. This demonstrates that the Pi release rate from F-actin is controlled by steric hindrance through the backdoor rather than by the disruption of the ionic bond between Pi and Mg2+ at the nucleotide-binding site. Our results provide the molecular basis for Pi release from actin and exemplify how a single, disease-linked point mutation distorts the nucleotide state distribution and atomic structure of the actin filament.
The interaction between the Heat Shock Proteins 70 and 40 is at the core of the ATPase regulation of the chaperone machinery that maintains protein homeostasis. However, the structural details of this fundamental interaction are still elusive and contrasting models have been proposed for the transient Hsp70/Hsp40 complexes. Here we combine molecular simulations based on both coarsegrained and atomistic models with co-evolutionary sequence analysis to shed light on this problem by focusing on the bacterial DnaK/DnaJ system. The integration of these complementary approaches resulted into a novel structural model that rationalizes previous experimental observations. We identify an evolutionary-conserved interaction surface formed by helix II of the DnaJ J-domain and a groove on lobe IIA of the DnaK nucleotide binding domain, involving the inter-domain linker.
Mosquito species belonging to the genus Aedes have attracted the interest of scientists and public health officers for their invasive species traits and efficient capacity of transmitting viruses affecting humans. Some of these species were brought outside their native range by human activities such as trade and tourism, and colonised new regions thanks to a unique combination of eco-physiological traits.
Considering mosquito physiological and behavioural traits to understand and predict the spatial and temporal population dynamics is thus a crucial step to develop strategies to mitigate the local densities of invasive Aedes populations.
Here, we synthesised the life cycle of four invasive Aedes species (Ae. aegypti, Ae. albopictus, Ae. japonicus and Ae. koreicus) in a single multi-scale stochastic modelling framework which we coded in the R package dynamAedes. We designed a stage-based and time-discrete stochastic model driven by temperature, photo-period and inter-specific larval competition that can be applied to three different spatial scales: punctual, local and regional. These spatial scales consider different degrees of spatial complexity and data availability, by accounting for both active and passive dispersal of mosquito species as well as for the heterogeneity of the input temperature data.
Our overarching aim was to provide a flexible, open-source and user-friendly tool rooted in the most updated knowledge on species biology which could be applied to the management of invasive Aedes populations as well as for more theoretical ecological inquiries.
The basidiomycete smut fungi are predominantly plant parasitic, causing severe losses in some crops. Most species feature a saprotrophic haploid yeast stage, and several smut fungi are only known from this stage, with some isolated from habitats without suitable hosts, e.g. from Antarctica. Thus, these species are generally believed to be apathogenic, but recent findings that some of these might have a plant pathogenic sexual counterpart, casts doubts on the validity of this hypothesis. Here, four Pseudozyma genomes were re-annotated and compared to published smut pathogens and the well-characterised effector gene Pep1 from these species was checked for its ability to complement a Pep1 deletion strain of Ustilago maydis. It was found that 113 high-confidence putative effector proteins were conserved among smut and Pseudozyma genomes. Among these were several validated effector proteins, including Pep1. By genetic complementation we show that Pep1 homologs from the supposedly apathogenic yeasts restore virulence in Pep1-deficient mutants Ustilago maydis. Thus, it is concluded that Pseudozyma species have retained a suite of effectors. This hints at the possibility that Pseudozyma species have kept an unknown plant pathogenic stage for sexual recombination or that these effectors have positive effects when colonising plant surfaces.
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.
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.
Vocal communication is essential to coordinate social interactions in mammals and it requires a fine discrimination of communication sounds. Auditory neurons can exhibit selectivity for specific calls, but how it is affected by preceding sounds is still debated. We tackled this using ethologically relevant vocalizations in a highly vocal mammalian species: Seba’s short-tailed bat. We show that cortical neurons present several degrees of selectivity for echolocation and distress calls. Embedding vocalizations within natural acoustic streams leads to stimulus-specific suppression of neuronal responses that changes sound selectivity in disparate manners: increases in neurons with poor discriminability in silence and decreases in neurons selective in silent settings. A computational model indicates that the observed effects arise from two forms of adaptation: presynaptic frequency specific adaptation acting in cortical inputs and stimulus unspecific postsynaptic adaptation. These results shed light into how acoustic context modulates natural sound discriminability in the mammalian cortex.
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.
Communication sounds are ubiquitous in the animal kingdom, where they play a role in advertising physiological states and/or socio-contextual scenarios. Distress sounds, for example, are typically uttered in distressful scenarios such as agonistic interactions. Here, we report on the occurrence of superfast temporal periodicities in distress calls emitted by bats (species Carollia perspicillata). Distress vocalizations uttered by this bat species are temporally modulated at frequencies close to 1.7 kHz, that is, ∼17 times faster than modulation rates observed in human screams. Fast temporal periodicities are represented in the bats’ brain by means of frequency following responses, and temporally periodic sounds are more effective in boosting the heart rate of awake bats than their demodulated versions. Altogether, our data suggest that bats, an animal group classically regarded as ultrasonic, can exploit the low frequency portion of the soundscape during distress calling to create spectro-temporally complex, arousing sounds.
Neural oscillations are at the core of important computations in the mammalian brain. Interactions between oscillatory activities in different frequency bands, such as delta (1-4 Hz), theta (4-8 Hz), or gamma (>30 Hz), are a powerful mechanism for binding fundamentally distinct spatiotemporal scales of neural processing. Phase-amplitude coupling (PAC) is one such plausible and well-described interaction, but much is yet to be uncovered regarding how PAC dynamics contribute to sensory representations. In particular, although PAC appears to have a major role in audition, the characteristics of coupling profiles in sensory and integration (i.e. frontal) cortical areas remain obscure. Here, we address this question by studying PAC dynamics in the frontal-auditory field (FAF; an auditory area in the bat frontal cortex) and the auditory cortex (AC) of the bat Carollia perspicillata. By means of simultaneous electrophysiological recordings in frontal and auditory cortices examining local-field potentials (LFPs), we show that the amplitude of gamma-band activity couples with the phase of low-frequency LFPs in both structures. Our results demonstrate that the coupling in FAF occurs most prominently in delta/high-gamma frequencies (1-4/75-100 Hz), whereas in the AC the coupling is strongest in the theta/low-gamma (2-8/25-55 Hz) range. We argue that distinct PAC profiles may represent different mechanisms for neuronal processing in frontal and auditory cortices, and might complement oscillatory interactions for sensory processing in the frontal-auditory cortex network.
In humans, screams have strong amplitude modulations (AM) at 30 to 150 Hz. These AM correspond to the acoustic correlate of perceptual roughness. In bats, distress calls can carry AMs, which elicit heart rate increases in playback experiments. Whether amplitude modulation occurs in fearful vocalisations of other animal species beyond humans and bats remains unknown. Here we analysed the AM pattern of rats’ 22-kHz ultrasonic vocalisations emitted in a fear conditioning task. We found that the number of vocalisations decreases during the presentation of conditioned stimuli. We also observed that AMs do occur in rat 22-kHz vocalisations. AMs are stronger during the presentation of conditioned stimuli, and during escape behaviour compared to freezing. Our results suggest that the presence of AMs in vocalisations emitted could reflect the animal’s internal state of fear related to avoidance behaviour.
Frontal areas of the mammalian cortex are thought to be important for cognitive control and complex behaviour. These areas have been studied mostly in humans, non-human primates and rodents. In this article, we present a quantitative characterization of response properties of a frontal auditory area responsive to sound in the bat brain, the frontal auditory field (FAF). Bats are highly vocal animals and they constitute an important experimental model for studying the auditory system. At present, little is known about neuronal sound processing in the bat FAF. We combined electrophysiology experiments and computational simulations to compare the response properties of auditory neurons found in the bat FAF and auditory cortex (AC) to simple sounds (pure tones). Anatomical studies have shown that the latter provide feedforward inputs to the former. Our results show that bat FAF neurons are responsive to sounds, however, when compared to AC neurons, they presented sparser, less precise spiking and longer-lasting responses. Based on the results of an integrate-and-fire neuronal model, we speculate that slow, low-threshold, synaptic dynamics could contribute to the changes in activity pattern that occur as information travels through cortico-cortical projections from the AC to the FAF.
Summary The auditory midbrain (inferior colliculus, IC) plays an important role in sound processing, acting as hub for acoustic information extraction and for the implementation of fast audio-motor behaviors. IC neurons are topographically organized according to their sound frequency preference: dorsal IC regions encode low frequencies while ventral areas respond best to high frequencies, a type of sensory map defined as tonotopy. Tonotopic maps have been studied extensively using artificial stimuli (pure tones) but our knowledge of how these maps represent information about sequences of natural, spectro-temporally rich sounds is sparse. We studied this question by conducting simultaneous extracellular recordings across IC depths in awake bats (Carollia perspicillata) that listened to sequences of natural communication and echolocation sounds. The hypothesis was that information about these two types of sound streams is represented at different IC depths since they exhibit large differences in spectral composition, i.e. echolocation covers the high frequency portion of the bat soundscape (> 45 kHz), while communication sounds are broadband and carry most power at low frequencies (20-25 kHz). Our results showed that mutual information between neuronal responses and acoustic stimuli, as well as response redundancy in pairs of neurons recorded simultaneously, increase exponentially with IC depth. The latter occurs regardless of the sound type presented to the bats (echolocation or communication). Taken together, our results indicate the existence of mutual information and redundancy maps at the midbrain level whose response cannot be predicted based on the frequency composition of natural sounds and classic neuronal tuning curves.
Transient receptor potential (TRP) ion channels are among the most well-studied classes of temperature-sensing molecules. Yet, the molecular mechanism and thermodynamic basis for the temperature sensitivity of TRP channels remains to this day poorly understood. One hypothesis is that the temperature-sensing mechanism can simply be described by a difference in heat capacity between the closed and open channel states. While such a two-state model may be simplistic it nonetheless has descriptive value, in the sense that it can be used to to compare overall temperature sensitivity between different channels and mutants. Here, we introduce a mathematical framework based on the two-state model to reliably extract temperature-dependent thermodynamic potentials and heat capacities from measurements of equilibrium constants at different temperatures. Our framework is implemented in an open-source data analysis package that provides a straightforward way to fit both linear and nonlinear van ‘t Hoff plots, thus avoiding some of the previous, potentially erroneous, assumptions when extracting thermodynamic variables from TRP channel electrophysiology data.
Animals extract behaviorally relevant signals from “noisy” environments. To investigate signal extraction, echolocating provides a rich system testbed. For orientation, bats broadcast calls and assign each echo to the corresponding call. When orienting in acoustically enriched environments or when approaching targets, bats change their spectro-temporal call design. Thus, to assess call adjustments that are exclusively meant to facilitate signal extraction in “noisy” environments, it is necessary to control for distance-dependent call changes. By swinging bats in a pendulum, we tested the influence of acoustic playback on the echolocation behavior of Carollia perspicillata. This paradigm evokes reproducible orientation behavior and allows a precise definition of the influence of the acoustic context. Our results show that bats dynamically switch between different adaptations to cope with sound-based navigation in acoustically contaminated environments. These dynamics of echolocation behavior may explain the large variety of adaptations that have been reported in the bat literature.
The development of super-resolution microscopy (SRM) has widened our understanding of biomolecular structure and function in biological materials. Imaging multiple targets within a single area would elucidate their spatial localization relative to the cell matrix and neighboring biomolecules, revealing multi-protein macromolecular structures and their functional co-dependencies. SRM methods are, however, limited to the number of suitable fluorophores that can be imaged during a single acquisition as well as the loss of antigens during antibody washing and restaining for organic dye multiplexing. We report the visualization of multiple protein targets within the pre- and postsynapse in 350-400 nm thick neuronal tissue sections using DNA-assisted single-molecule localization microscopy. Using antibodies labeled with short DNA oligonucleotides, multiple targets are visualized successively by sequential exchange of fluorophore-labeled complementary oligonucleotides present in the imaging buffer. The structural integrity of the tissue is maintained owing to only a single labelling step during sample preparation. Multiple targets are imaged using a single laser wavelength, minimizing chromatic aberration. This method proved robust for multi-target imaging in semi-thin tissue sections, paving the way towards structural cell biology with single-molecule super-resolution microscopy.
Understanding the nano-architecture of protein machines in diverse sub-cellular compartments remains a challenge despite rapid progress in super-resolution microscopy. While singlemolecule localization microscopy techniques allow the visualization and identification of cellular structures with near-molecular resolution, multiplex-labeling of tens of target proteins within the same sample has not yet been achieved routinely. However, single sample multiplexing is essential to detect patterns that threaten to get lost in multi-sample averaging. Here, we report maS3TORM (multiplexed automated serial staining stochastic optical reconstruction microscopy), a microscopy approach capable of fully automated 3D dSTORM imaging and solution exchange employing a re-staining protocol to achieve highly multiplexed protein localization within individual biological samples. We demonstrate 3D super-resolution images of 15 target proteins in single cultured cells and 16 targets in individual neuronal tissue samples with <10 nm localization precision. This allowed us to define novel nano-architectural features of protein distribution within the presynaptic nerve terminal.
Membrane receptors are central to cell-cell communication. Receptor clustering at the plasma membrane modulates physiological responses, and mesoscale receptor organization is critical for downstream signaling. Spatially restricted cluster formation of the neuropeptide Y2 hormone receptor (Y2R) was observed in vivo; however, the relevance of this confinement is not fully understood. Here, we controlled Y2R clustering in situ by a chelator nanotool. Due to the multivalent interaction, we observed a dynamic exchange in the microscale confined regions. Fast Y2R enrichment in clustered areas triggered a ligand-independent downstream signaling determined by an increase in cytosolic calcium, cell spreading, and migration. We revealed that the cell response to ligand-induced activation was amplified when cells were pre-clustered by the nanotool. Ligand-independent signaling by clustering differed from ligand-induced activation in the binding of arrestin-3 as downstream effector, which was recruited to the confined regions only in the presence of the ligand. This approach enables in situ clustering of membrane receptors and raises the possibility to explore different modalities of receptor activation.
Molecular recognition of M1-linked ubiquitin chains by native and phosphorylated UBAN domains
(2019)
Although the Ub-binding domain in ABIN proteins and NEMO (UBAN) is highly conserved, UBAN-containing proteins exhibit different Ub-binding properties, resulting in their diverse biological roles. Post-translational modifications further control UBAN domain specificity for poly-Ub chains. However, precisely, how the UBAN domain structurally confers such functional diversity remains poorly understood. Here we report crystal structures of ABIN-1 alone and in complex with one or two M1-linked di-Ub chains. ABIN-1 UBAN forms a homo-dimer that provides two symmetrical Ub-binding sites on either side of the coiled-coil structure. Moreover, crystal structures of ABIN1 UBAN in complex with di-Ub chains reveal a concentration-dependency of UBAN/di-Ub binding stoichiometry. Analysis of UBAN/M1-linked di-Ub binding characteristics indicates that phosphorylated S473 in OPTN and its corresponding phospho-mimetic residue in ABIN-1 (E484) are essential for high affinity interactions with M1-linked Ub chains. Also, a phospho-mimetic mutation of A303 in NEMO, corresponding to S473 of OPTN, increases binding affinity for M1-linked Ub chains. These findings are in line with the diverse physiological roles of UBAN domains, as phosphorylation of OPTN UBAN is required to enhance its binding to Ub during mitophagy.
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.
From loss to recovery: how to effectively assess chemosensory impairments during COVID-19 pandemic
(2021)
Chemosensory impairments have been established as a specific indicator of COVID-19. They affect most patients and may persist long past the resolution of respiratory symptoms, representing an unprecedented medical challenge. Since the SARS-CoV-2 pandemic started, we now know much more about smell, taste, and chemesthesis loss associated with COVID-19. However, the temporal dynamics and characteristics of recovery are still unknown. Here, capitalizing on data from the Global Consortium for Chemosensory Research (GCCR) crowdsourced survey, we assessed chemosensory abilities after the resolution of respiratory symptoms in participants diagnosed with COVID-19 during the first wave of the pandemic in Italy. This analysis led to the identification of two patterns of chemosensory recovery, limited (partial) and substantial, which were found to be associated with differential age, degrees of chemosensory loss, and regional patterns. Uncovering the self-reported phenomenology of recovery from smell, taste, and chemesthetic disorders is the first, yet essential step, to provide healthcare professionals with the tools to take purposeful and targeted action to address chemosensory disorders and its severe discomfort.
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.
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.
Mathematical modeling of the molecular switch of TNFR1-mediated signaling pathways using Petri nets
(2021)
The paper describes a mathematical model of the molecular switch of cell survival, apoptosis, and necroptosis in cellular signaling pathways initiated by tumor necrosis factor 1. Based on experimental findings in the current literature, we constructed a Petri net model in terms of detailed molecular reactions for the molecular players, protein complexes, post-translational modifications, and cross talk. The model comprises 118 biochemical entities, 130 reactions, and 299 connecting edges. Applying Petri net analysis techniques, we found 279 pathways describing complete signal flows from receptor activation to cellular response, representing the combinatorial diversity of functional pathways.120 pathways steered the cell to survival, whereas 58 and 35 pathways led to apoptosis and necroptosis, respectively. For 65 pathways, the triggered response was not deterministic, leading to multiple possible outcomes. Based on the Petri net, we investigated the detailed in silico knockout behavior and identified important checkpoints of the TNFR1 signaling pathway in terms of ubiquitination within complex I and the gene expression dependent on NF-κB, which controls the caspase activity in complex II and apoptosis induction.
Autophagosome biogenesis requires a localized perturbation of lipid membrane dynamics and a unique protein-lipid conjugate. Autophagy-related (ATG) proteins catalyze this biogenesis on cellular membranes, but the underlying molecular mechanism remains unclear. Focusing on the final step of the protein-lipid conjugation reaction, ATG8/LC3 lipidation, we show how membrane association of the conjugation machinery is organized and fine-tuned at the atomistic level. Amphipathic α-helices in ATG3 proteins (AHATG3) are found to have low hydrophobicity and to be less bulky. Molecular dynamics simulations reveal that AHATG3 regulates the dynamics and accessibility of the thioester bond of the ATG3∼LC3 conjugate to lipids, allowing covalent lipidation of LC3. Live cell imaging shows that the transient membrane association of ATG3 with autophagic membranes is governed by the less bulky- hydrophobic feature of AHATG3. Collectively, the unique properties of AHATG3 facilitate protein- lipid bilayer association leading to the remodeling of the lipid bilayer required for the formation of autophagosomes.
Molluscs are the second most species-rich phylum in the animal kingdom, yet only eleven genomes of this group have been published so far. Here, we present the draft genome sequence of the pulmonate freshwater snail Radix auricularia. Six whole genome shotgun libraries with different layouts were sequenced. The resulting assembly comprises 4,823 scaffolds with a cumulative length of 910 Mb and an overall read coverage of 72x. The assembly contains 94.6 % of a metazoan core gene collection, indicating an almost complete coverage of the coding fraction. The discrepancy of ~690 Mb compared to the estimated genome size of R. auricularia (1.6 Gb) results from a high repeat content of 70 % mainly comprising DNA transposons. The annotation of 17,338 protein coding genes was supported by the use of publicly-available transcriptome data. This draft will serve as starting point for further genomic and population genetic research in this scientifically important phylum.
It becomes more and more obvious that deregulation of host metabolism play an important role in SARS-CoV-2 pathogenesis with implication for increased risk of severe course of COVID-19. Furthermore, it is expected that COVID-19 patients recovered from severe disease may experience long-term metabolic disorders. Thereby understanding the consequences of SARS-CoV-2 infection on host metabolism can facilitate efforts for effective treatment option. We have previously shown that SARS-CoV-2-infected cells undergo a shift towards glycolysis and that 2-deoxy-D-glucose (2DG) inhibits SARS-CoV-2 replication. Here, we show that also pentose phosphate pathway (PPP) is remarkably deregulated. Since PPP supplies ribonucleotides for SARS-CoV-2 replication, this could represent an attractive target for an intervention. On that account, we employed the transketolase inhibitor benfooxythiamine and showed dose-dependent inhibition of SARS-CoV-2 in non-toxic concentrations. Importantly, the antiviral efficacy of benfooxythiamine was further increased in combination with 2DG.
Endogenous clocks enable organisms to adapt their physiology and behavior to daily variation in environmental conditions. Metabolic processes in cyanobacteria to humans are effected by the circadian clock, and its dysregulation causes metabolic disorders. In mouse and Drosophila were shown that the circadian clock directs translation of factors involved in ribosome biogenesis and synchronizes protein synthesis. However, the role of clocks in Drosophila neurogenesis and the potential impact of clock impairment on neural circuit formation and function is less understood. Here we demonstrate that light stimuli or circadian clock causes a defect in neural stem cell growth and proliferation accompanied by reduced nucleolar size. Further, we define that light and clock independently affect the InR/TOR growth regulatory pathway due to the effect on regulators of protein biosynthesis. Altogether, these data suggest that alterations in growth regulatory pathways induced by light and clock are associated with impaired neural development.
Under natural conditions, the visual system often sees a given input repeatedly. This provides an opportunity to optimize processing of the repeated stimuli. Stimulus repetition has been shown to strongly modulate neuronal-gamma band synchronization, yet crucial questions remained open. Here we used magnetoencephalography in 30 human subjects and find that gamma decreases across ~10 repetitions and then increases across further repetitions, revealing plastic changes of the activated neuronal circuits. Crucially, changes induced by one stimulus did not affect responses to other stimuli, demonstrating stimulus specificity. Changes partially persisted when the inducing stimulus was repeated after 25 minutes of intervening stimuli. They were strongest in early visual cortex and increased interareal feedforward influences. Our results suggest that early visual cortex gamma synchronization enables adaptive neuronal processing of recurring stimuli. These and previously reported changes might be due to an interaction of oscillatory dynamics with established synaptic plasticity mechanisms.
The firing pattern of ventral midbrain dopamine neurons is controlled by afferent and intrinsic activity to generate prediction error signals that are essential for reward-based learning. Given the absence of intracellular in vivo recordings in the last three decades, the subthreshold membrane potential events that cause changes in dopamine neuron firing patterns remain unknown. By establishing stable in vivo whole-cell recordings of >100 spontaneously active midbrain dopamine neurons in anaesthetized mice, we identified the repertoire of subthreshold membrane potential signatures associated with distinct in vivo firing patterns. We demonstrate that dopamine neuron in vivo activity deviates from a single spike pacemaker pattern by eliciting transient increases in firing rate generated by at least two diametrically opposing biophysical mechanisms: a transient depolarization resulting in high frequency plateau bursts associated with a reactive, depolarizing shift in action potential threshold; and a prolonged hyperpolarization preceding slower rebound bursts characterized by a predictive, hyperpolarizing shift in action potential threshold. Our findings therefore illustrate a framework for the biophysical implementation of prediction error and sensory cue coding in dopamine neurons by tuning action potential threshold dynamics.
Genome-wide CRISPR screens are becoming more widespread and allow the simultaneous interrogation of thousands of genomic regions. Although recent progress has been made in the analysis of CRISPR screens, it is still an open problem how to interpret CRISPR mutations in non-coding regions of the genome. Most of the tools concentrate on the interpretation of mutations introduced in gene coding regions. We introduce a computational pipeline that uses epigenomic information about regulatory elements for the interpretation of CRISPR mutations in non-coding regions. We illustrate our approach on the analysis of a genome-wide CRISPR screen in hTERT-RPE-1 cells and reveal novel regulatory elements that mediate chemoresistance against doxorubicin in these cells. We infer links to established and to novel chemoresistance genes. Our approach is general and can be applied on any cell type and with different CRISPR enzymes.
Several studies have probed perceptual performance at different times after a self-paced motor action and found frequency-specific modulations of perceptual performance phase-locked to the action. Such action-related modulation has been reported for various frequencies and modulation strengths. In an attempt to establish a basic effect at the population level, we had a relatively large number of participants (n=50) perform a self-paced button press followed by a detection task at threshold, and we applied both fixed- and random-effects tests. The combined data of all trials and participants surprisingly did not show any significant action-related modulation. However, based on previous studies, we explored the possibility that such modulation depends on the participant’s internal state. Indeed, when we split trials based on performance in neighboring trials, then trials in periods of low performance showed an action-related modulation at ≈17 Hz. When we split trials based on the performance in the preceding trial, we found that trials following a “miss” showed an action-related modulation at ≈17 Hz. Finally, when we split participants based on their false-alarm rate, we found that participants with no false alarms showed an action-related modulation at ≈17 Hz. All these effects were significant in random-effects tests, supporting an inference on the population. Together, these findings indicate that action-related modulations are not always detectable. However, the results suggest that specific internal states such as lower attentional engagement and/or higher decision criterion are characterized by a modulation in the beta-frequency range.
Several recent studies investigated the rhythmic nature of cognitive processes that lead to perception and behavioral report. These studies used different methods, and there has not yet been an agreement on a general standard. Here, we present a way to test and quantitatively compare these methods. We simulated behavioral data from a typical experiment and analyzed these data with several methods. We applied the main methods found in the literature, namely sine-wave fitting, the Discrete Fourier Transform (DFT) and the Least Square Spectrum (LSS). DFT and LSS can be applied both on the averaged accuracy time course and on single trials. LSS is mathematically equivalent to DFT in the case of regular, but not irregular sampling - which is more common. LSS additionally offers the possibility to take into account a weighting factor which affects the strength of the rhythm, such as arousal. Statistical inferences were done either on the investigated sample (fixed-effect) or on the population (random-effect) of simulated participants. Multiple comparisons across frequencies were corrected using False-Discovery-Rate, Bonferroni, or the Max-Based approach. To perform a quantitative comparison, we calculated Sensitivity, Specificity and D-prime of the investigated analysis methods and statistical approaches. Within the investigated parameter range, single-trial methods had higher sensitivity and D-prime than the methods based on the averaged-accuracy-time-course. This effect was further increased for a simulated rhythm of higher frequency. If an additional (observable) factor influenced detection performance, adding this factor as weight in the LSS further improved Sensitivity and D-prime. For multiple comparison correction, the Max-Based approach provided the highest Specificity and D-prime, closely followed by the Bonferroni approach. Given a fixed total amount of trials, the random-effect approach had higher D-prime when trials were distributed over a larger number of participants, even though this gave less trials per participant. Finally, we present the idea of using a dampened sinusoidal oscillator instead of a simple sinusoidal function, to further improve the fit to behavioral rhythmicity observed after a reset event.
Analyzing non-invasive recordings of electroencephalography (EEG) and magnetoencephalography (MEG) directly in sensor space, using the signal from individual sensors, is a convenient and standard way of working with this type of data. However, volume conduction introduces considerable challenges for sensor space analysis. While the general idea of signal mixing due to volume conduction in EEG/MEG is recognized, the implications have not yet been clearly exemplified. Here, we illustrate how different types of activity overlap on the level of individual sensors. We show spatial mixing in the context of alpha rhythms, which are known to have generators in different areas of the brain. Using simulations with a realistic 3D head model and lead field and data analysis of a large resting-state EEG dataset, we show that electrode signals can be differentially affected by spatial mixing by computing a sensor complexity measure. While prominent occipital alpha rhythms result in less heterogeneous spatial mixing on posterior electrodes, central electrodes show a diversity of rhythms present. This makes the individual contributions, such as the sensorimotor mu-rhythm and temporal alpha rhythms, hard to disentangle from the dominant occipital alpha. Additionally, we show how strong occipital rhythms rhythms can contribute the majority of activity to frontal channels, potentially compromising analyses that are solely conducted in sensor space. We also outline specific consequences of signal mixing for frequently used assessment of power, power ratios and connectivity profiles in basic research and for neurofeedback application. With this work, we hope to illustrate the effects of volume conduction in a concrete way, such that the provided practical illustrations may be of use to EEG researchers to in order to evaluate whether sensor space is an appropriate choice for their topic of investigation.
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.
In the course of global climate change, central Europe is experiencing more frequent and prolonged periods of drought. The drought years 2018 and 2019 affected European beeches (Fagus sylvatica L.) differently: even in the same stand, drought damaged trees neighboured healthy trees, suggesting that the genotype rather than the environment was responsible for this conspicuous pattern. We used this natural experiment to study the genomic basis of drought resistance with Pool-GWAS. Contrasting the extreme phenotypes identified 106 significantly associated SNPs throughout the genome. Most annotated genes with associated SNPs (>70%) were previously implicated in the drought reaction of plants. Non-synonymous substitutions led either to a functional amino acid exchange or premature termination. A SNP-assay with 70 loci allowed predicting drought phenotype in 98.6% of a validation sample of 92 trees. Drought resistance in European beech is a moderately polygenic trait that should respond well to natural selection, selective management, and breeding.
The gradual heterogeneity of climatic factors pose varying selection pressures across geographic distances that leave signatures of clinal variation in the genome. Separating signatures of clinal adaptation from signatures of other evolutionary forces, such as demographic processes, genetic drift, and adaptation to non-clinal conditions of the immediate local environment is a major challenge. Here, we examine climate adaptation in five natural populations of the harlequin fly Chironomus riparius sampled along a climatic gradient across Europe. Our study integrates experimental data, individual genome resequencing, Pool-Seq data, and population genetic modelling. Common-garden experiments revealed a positive correlation of population growth rates corresponding to the population origin along the climate gradient, suggesting thermal adaptation on the phenotypic level. Based on a population genomic analysis, we derived empirical estimates of historical demography and migration. We used an FST outlier approach to infer positive selection across the climate gradient, in combination with an environmental association analysis. In total we identified 162 candidate genes as genomic basis of climate adaptation. Enriched functions among these candidate genes involved the apoptotic process and molecular response to heat, as well as functions identified in other studies of climate adaptation in other insects. Our results show that local climate conditions impose strong selection pressures and lead to genomic adaptation despite strong gene flow. Moreover, these results imply that selection to different climatic conditions seems to converge on a functional level, at least between different insect species.
One of the major problems in evolutionary biology is to elucidate the relationships between historical events and the tempo and mode of lineage divergence. The development of relaxed molecular clock models and the increasing availability of DNA sequences resulted in more accurate estimations of taxa divergence times. However, finding the link between competing historical events and divergence is still challenging. Here we investigate assigning constrained-age priors to nodes of interest in a time-calibrated phylogeny as a means of hypothesis comparison. These priors are equivalent to historic scenarios for lineage origin. The hypothesis that best explains the data can be selected by comparing the likelihood values of the competing hypotheses, modelled with different priors. A simulation approach was taken to evaluate the performance of the prior-based method and to compare it with an unconstrained approach. We explored the effect of DNA sequence length and the temporal placement and span of competing hypotheses (i.e. historic scenarios) on selection of the correct hypothesis and the strength of the inference. Competing hypotheses were compared applying a posterior simulation analogue of the Akaike Information Criterion and Bayes factors (obtained after calculation of the marginal likelihood with three estimators: Harmonic Mean, Stepping Stone and Path Sampling). We illustrate the potential application of the prior-based method on an empirical data set to compare competing geological hypotheses explaining the biogeographic patterns in Pleurodeles newts. The correct hypothesis was selected on average 89% times. The best performance was observed with DNA sequence length of 3500-10000 bp. The prior-based method is most reliable when the hypotheses compared are not temporally too close. The strongest inferences were obtained when using the Stepping Stone and Path Sampling estimators. The prior-based approach proved effective in discriminating between competing hypotheses when used on empirical data. The unconstrained analyses performed well but it probably requires additional computational effort. Researchers applying this approach should rely only on inferences with moderate to strong support. The prior-based approach could be applied on biogeographical and phylogeographical studies where robust methods for historical inferences are still lacking.
Active transposable elements (TEs) may result in divergent genomic insertion and abundance patterns among conspecific populations. Upon secondary contact, such divergent genetic backgrounds can theoretically give rise to classical Dobzhansky-Muller incompatibilities (DMI), a way how TEs can contribute to the evolution of endogenous genetic barriers and eventually population divergence. We investigated whether differential TE activity created endogenous selection pressures among conspecific populations of the non-biting midge Chironomus riparius, focussing on a Chironomus-specific TE, the minisatellite-like Cla-element, whose activity is associated with speciation in the genus. Using an improved and annotated draft genome for a genomic study with five natural C. riparius populations, we found highly population-specific TE insertion patterns with many private insertions. A highly significant correlation of pairwise population FST from genome-wide SNPs with the FST estimated from TEs suggests drift as the major force driving TE population differentiation. However, the significantly higher Cla-element FST level due to a high proportion of differentially fixed Cla-element insertions indicates that segregating, i.e. heterozygous insertions are selected against. With reciprocal crossing experiments and fluorescent in-situ hybridisation of Cla-elements to polytene chromosomes, we documented phenotypic effects on female fertility and chromosomal mispairings that might be linked to DMI in hybrids. We propose that the inferred negative selection on heterozygous Cla-element insertions causes endogenous genetic barriers and therefore acts as DMI among C. riparius populations. The intrinsic genomic turnover exerted by TEs, thus, may have a direct impact on population divergence that is operationally different from drift and local adaptation.
Bears are iconic mammals with a complex evolutionary history. Natural bear hybrids and studies of few nuclear genes indicate that gene flow among bears may be more common than expected and not limited to the closely related polar and brown bears. Here we present a genome analysis of the bear family with representatives of all living species. Phylogenomic analyses of 869 mega base pairs divided into 18,621 genome fragments yielded a well-resolved coalescent species tree despite signals for extensive gene flow across species. However, genome analyses using three different statistical methods show that gene flow is not limited to closely related species pairs. Strong ancestral gene flow between the Asiatic black bear and the ancestor to polar, brown and American black bear explains numerous uncertainties in reconstructing the bear phylogeny. Gene flow across the bear clade may be mediated by intermediate species such as the geographically wide-spread brown bears leading to massive amounts of phylogenetic conflict. Genome-scale analyses lead to a more complete understanding of complex evolutionary processes. The increasing evidence for extensive inter-specific gene flow, found also in other animal species, necessitates shifting the attention from speciation processes achieving genome-wide reproductive isolation to the selective processes that maintain species divergence in the face of gene flow.
All giraffe (Giraffa) were previously assigned to a single species (G. Camelopardalis) and nine subspecies. However, multi-locus analyses of all subspecies have shown that there are four genetically distinct clades and suggest four giraffe species. This conclusion might not be fully accepted due to limited data and lack of explicit gene flow analyses. Here we present an extended study based on 21 independent nuclear loci from 137 individuals. Explicit gene flow analyses identify less than one migrant per generation, including between the closely related northern and reticulated giraffe. Thus, gene flow analyses and population genetics of the extended dataset confirm four genetically distinct giraffe clades and support four independent giraffe species. The new findings call for a revision of the IUCN classification of giraffe taxonomy. Three of the four species are threatened with extinction, mostly occurring in politically unstable regions, and as such, require the highest conservation support possible.
Welfare is the largest expenditure category in all advanced democracies. Consequently, much literature has studied partisan effects on total and policy-specific welfare expenditure. Yet, these results cannot be trusted: The methodological standard is to apply time-series cross-section-regressions to annual observation data. But governments hardly change annually. Thus, the number of observations is artificially inflated, leading to incorrect estimates. While this problem has recently been acknowledged, it has not been convincingly resolved. We propose Mixed-Effects Models as a solution, which allow decomposing variance into different levels and permit complex cross-classification data structures. We argue that Mixed-Effects models combine the strengths of existing methodological approaches while alleviating their weaknesses. Empirically, we study partisan effects on total and on disaggregated expenditure in 23 OECD-countries, 1960-2012, using several measures of party preferences.
The debate on effects of globalization on welfare states is extensive. Often couched in terms of a battle between the compensation and the efficiency thesis, the scholarly literature has provided contradictory arguments and findings. This article contributes to the scholarly debate by exploring in greater detail the micro-level foundations of compensation theory. More specifically, we distinguish between individual policy preferences for compensatory social policies (unemployment insurance) and human capital-focused social investment policies (education) and expect globalization to mainly affect demand for educational investment. A multi-level analysis of ISSP survey data provides empirical support for this hypothesis. This finding provides an important revision and extension of the classical analytical perspective of compensation theory, because it shows that citizens value the social investment function of the welfare state above and beyond simple compensation via social transfers. This might be particularly relevant in today's skill-centered knowledge economies.
We cannot imagine a political system without opposition. Despite this crucial position in politics, political science has largely neglected to study oppositions. Attempting to fill this gap, this article analyses the institutional opportunities of parliamentary oppositions. It offers a parsimonious framework by distinguishing two dimensions of opposition influence: Some institutions enable oppositions to control governments, while others offer opportunities to present alternatives. A comparison of oppositions’ opportunities in 21 democracies shows that countries fall into four groups along these dimensions: In majoritarian democracies, weak control mechanisms are countered by excellent opportunities to publicize alternatives. Consociational democracies are characterized by strong control mechanisms, but provide only weak opportunities to present alternatives. Moreover, in Southern Europe, control mechanisms and opportunities to present alternatives are weak, while both are pronounced in Nordic Europe. The results are summarized in three indices that can easily be applied in future research examining oppositions and their power.
Phylogenetic analyses of nuclear and mitochondrial genomes have shown that polar bears captured the mitochondrial genome of brown bears some 160,00 years ago. This hybridization event likely led to an extinction of the original polar bear mitochondrial genome. However, parts of the mitochondrial DNA occasionally integrates into the nuclear genome, forming pseudogenes called numts (nuclear mitochondrial integrations). Screening the polar bear genome for numts, we identified only 13 such integrations. Analyses of whole-genome sequences from additional polar bears, brown and American black bears as well as the giant panda indicates that the discovered numts entered the bear lineage before the initial ursid radiation some 14 million years ago. Our findings suggests a low integration rate of numts in the bear lineage and a complete loss of the original polar bear mitochondrial genome.
Background: Ever decreasing costs along with advances in sequencing and library preparation technologies enable even small research groups to generate chromosome-level assemblies today. Here we report the generation of an improved chromosome-level assembly for the Siamese fighting fish (Betta splendens) that was carried out during a practical university Master’s course. The Siamese fighting fish is a popular aquarium fish and an emerging model species for research on aggressive behaviour. We updated the current genome assembly by generating a new long-read nanopore-based assembly with subsequent scaffolding to chromosome-level using previously published HiC data.
Findings: The use of nanopore-based long-read data sequenced on a MinION platform (Oxford Nanopore Technologies) allowed us to generate a baseline assembly of only 1,276 contigs with a contig N50 of 2.1 Mbp, and a total length of 441 Mbp. Scaffolding using previously published HiC data resulted in 109 scaffolds with a scaffold N50 of 20.7 Mbp. More than 99% of the assembly is comprised in 21 scaffolds. The assembly showed the presence of 95.8% complete BUSCO genes from the Actinopterygii dataset indicating a high quality of the assembly.
Conclusion: We present an improved full chromosome-level assembly of the Siamese fighting fish generated during a university Master’s course. The use of ~35× long-read nanopore data drastically improved the baseline assembly in terms of continuity. We show that relatively in-expensive high-throughput sequencing technologies such as the long-read MinION sequencing platform can be used in educational settings allowing the students to gain practical skills in modern genomics and generate high quality results that benefit downstream research projects.
Borders and edges are salient and behaviourally relevant features for navigating the environment. The brain forms dedicated neural representations of environmental boundaries, which are assumed to serve as a reference for spatial coding. Here we expand this border coding network to include the retrosplenial cortex (RSC) in which we identified neurons that increase their firing near all boundaries of an arena. RSC border cells specifically encode walls, but not objects, and maintain their tuning in the absence of direct sensory detection. Unlike border cells in the medial entorhinal cortex (MEC), RSC border cells are sensitive to the animal’s direction to nearby walls located contralateral to the recorded hemisphere. Pharmacogenetic inactivation of MEC led to a disruption of RSC border coding, but not vice versa, indicating network directionality. Together these data shed light on how information about distance and direction of boundaries is generated in the brain for guiding navigation behaviour.
Improved integration of single cell transcriptome data demonstrated on heart failure in mice and men
(2023)
Biomedical research frequently uses murine models to study disease mechanisms. However, the translation of these findings to human disease remains a significant challenge. In order to improve the comparability of mouse and human data, we present a cross-species integration pipeline for single-cell transcriptomic assays.
The pipeline merges expression matrices and assigns clear orthologous relationships. Starting from Ensembl ortholog assignments, we allocated 82% of mouse genes to unique orthologs by using additional publicly available resources such as Uniprot, and NCBI databases. For genes with multiple matches, we employed the Needleman-Wunsch global alignment based on either amino acid or nucleotide sequence to identify the ortholog with the highest degree of similarity.
The workflow was tested for its functionality and efficiency by integrating scRNA-seq datasets from heart failure patients with the corresponding mouse model. We were able to assign unique human orthologs to up to 80% of the mouse genes, utilizing the known 17,492 orthologous pairs. Curiously, the integration process enabled the identification of both common and unique regulatory pathways between species in heart failure.
In conclusion, our pipeline streamlines the integration process, enhances gene nomenclature alignment and simplifies the translation of mouse models to human disease. We have made the OrthoIntegrate R-package accessible on GitHub (https://github.com/MarianoRuzJurado/OrthoIntegrate), which includes the assignment of ortholog definitions for human and mouse, as well as the pipeline for integrating single cells.
Recently, significant advances have been made by identifying the levels of synchronicity of the underlying dynamics of a given brain state. This research has demonstrated that unconscious dynamics tend to be more synchronous than those found in conscious states, which are more asynchronous. Here we go beyond this dichotomy to demonstrate that the different brain states are always underpinned by spatiotemporal chaos but with dissociable turbulent dynamics. We investigated human neuroimaging data from different brain states (resting state, meditation, deep sleep, and disorders of consciousness after coma) and were able to distinguish between them using complementary model-free and model-based measures of turbulent information transmission. Our model-free approach used recent advances describing a measure of information cascade across spatial scales using tools from turbulence theory. Complementarily, our model-based approach used exhaustive in silico perturbations of whole-brain models fitted to the empirical neuroimaging data, which allowed us to study the information encoding capabilities of the brain states. Overall, the current framework demonstrates that different levels of turbulent dynamics are fundamental for describing and differentiating between brain states.
BOLD signatures of sleep
(2019)
Sleep can be distinguished from wake by changes in brain electrical activity, typically assessed using electroencephalography (EEG). The hallmark of non-rapid-eye-movement sleep are two major EEG events: slow waves and spindles. Here we sought to identify possible signatures of sleep in brain hemodynamic activity, using simultaneous fMRI-EEG. We found that, during the transition from wake to sleep, blood-oxygen-level-dependent (BOLD) activity evolved from a mixed-frequency pattern to one dominated by two distinct oscillations: a low-frequency (~0.05Hz) oscillation prominent in light sleep and a high-frequency (~0.17Hz) oscillation in deep sleep. The two BOLD oscillations correlated with the occurrences of spindles and slow waves, respectively. They were detectable across the whole brain, cortically and subcortically, but had different regional distributions and opposite onset patterns. These spontaneous BOLD oscillations provide fMRI signatures of basic sleep processes, which may be employed to study human sleep at spatial resolution and brain coverage not achievable using EEG.
Summary: Understanding the role of short-interfering RNA (siRNA) in diverse biological processes is of current interest and often approached through small RNA sequencing. However, analysis of these datasets is difficult due to the complexity of biological RNA processing pathways, which differ between species. Several properties like strand specificity, length distribution, and distribution of soft-clipped bases are few parameters known to guide researchers in understanding the role of siRNAs. We present RAPID, a generic eukaryotic siRNA analysis pipeline, which captures information inherent in the datasets and automatically produces numerous visualizations as user-friendly HTML reports, covering multiple categories required for siRNA analysis. RAPID also facilitates an automated comparison of multiple datasets, with one of the normalization techniques dedicated for siRNA knockdown analysis, and integrates differential expression analysis using DESeq2. RAPID is available under MIT license at https://github.com/SchulzLab/RAPID. We recommend using it as a conda environment available from https://anaconda.org/bioconda/rapid.
Several studies suggested that transcription factor (TF) binding to DNA may be impaired or enhanced by DNA methylation. We present MeDeMo, a toolbox for TF motif analysis that combines information about DNA methylation with models capturing intra-motif dependencies. In a large-scale study using ChIP-seq data for 335 TFs, we identify novel TFs that are affected by DNA methylation. Overall, we find that CpG methylation decreases the likelihood of binding for the majority of TFs. For a considerable subset of TFs, we show that intra-motif dependencies are pivotal for accurately modelling the impact of DNA methylation on TF binding.
Background Enhancers play a fundamental role in orchestrating cell state and development. Although several methods have been developed to identify enhancers, linking them to their target genes is still an open problem. Several theories have been proposed on the functional mechanisms of enhancers, which triggered the development of various methods to infer promoter enhancer interactions (PEIs). The advancement of high-throughput techniques describing the three-dimensional organisation of the chromatin, paved the way to pinpoint long-range PEIs. Here we investigated whether including PEIs in computational models for the prediction of gene expression improves performance and interpretability.
Results We have extended our Tepic framework to include DNA contacts deduced from chromatin conformation capture experiments and compared various methods to determine PEIs using predictive modelling of gene expression from chromatin accessibility data and predicted transcription factor (TF) motif data. We found that including long-range PEIs deduced from both HiC and HiChIP data indeed improves model performance. We designed a novel machine learning approach that allows to prioritize TFs in distal loop and promoter regions with respect to their importance for gene expression regulation. Our analysis revealed a set of core TFs that are part of enhancer-promoter loops involving YY1 in different cell lines.
Conclusion: We show that the integration of chromatin conformation data improves gene expression prediction, underlining the importance of enhancer looping for gene expression regulation. Our general approach can be used to prioritize TFs that are involved in distal and promoter-proximal regulation using accessibility, conformation and expression data.
Interest in time-resolved connectivity in fMRI has grown rapidly in recent years. The most widely used technique for studying connectivity changes over time utilizes a sliding windows approach. There has been some debate about the utility of shorter versus longer windows, the use of fixed versus adaptive windows, as well as whether observed resting state dynamics during wakefulness may be predominantly due to changes in sleep state and subject head motion. In this work we use an independent component analysis (ICA)-based pipeline applied to concurrent EEG/fMRI data collected during wakefulness and various sleep stages and show: 1) connectivity states obtained from clustering sliding windowed correlations of resting state functional network time courses well classify the sleep states obtained from EEG data, 2) using shorter sliding windows instead of longer non-overlapping windows improves the ability to capture transition dynamics even at windows as short as 30 seconds, 3) motion appears to be mostly associated with one of the states rather than spread across all of them 4) a fixed tapered sliding window approach outperforms an adaptive dynamic conditional correlation approach, and 5) consistent with prior EEG/fMRI work, we identify evidence of multiple states within the wakeful condition which are able to be classified with high accuracy. Classification of wakeful only states suggest the presence of time-varying changes in connectivity in fMRI data beyond sleep state or motion. Results also inform about advantageous technical choices, and the identification of different clusters within wakefulness that are separable suggest further studies in this direction.
Mitochondrial NADH:ubiquinone oxidoreductase (complex I) is a 1 MDa membrane protein complex with a central role in energy metabolism. Redox-driven proton translocation by complex I contributes substantially to the proton motive force that drives ATP synthase. Several structures of complex I from bacteria and mitochondria have been determined but its catalytic mechanism has remained controversial. We here present the cryo-EM structure of complex I from Yarrowia lipolytica at 2.1 Å resolution, which reveals the positions of more than 1600 protein-bound water molecules, of which ∼100 are located in putative proton translocation pathways. Another structure of the same complex under steady-state activity conditions at 3.4 Å resolution indicates conformational transitions that we associate with proton injection into the central hydrophilic axis. By combining high-resolution structural data with site-directed mutagenesis and large-scale molecular dynamics simulations, we define details of the proton translocation pathways, and offer new insights into the redox-coupled proton pumping mechanism of complex I.
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.
Human behaviour is inextricably linked to the interaction of emotion and cognition. For decades, emotion and cognition were perceived as separable processes, yet with mutual interactions. Recently, this differen-tiation has been challenged by more integrative approaches, but without addressing the exact neurophysiological basis of their interaction. Here, we aimed to uncover neurophysiological mechanisms of emotion-cognition interaction. We used an emotional Flanker task paired with EEG/FEM beamforming in a large cohort (N=121) of healthy human participants, obtaining high temporal and fMRI-equivalent spatial resolution. Spatially, emotion and cognition processing overlapped in the right inferior frontal gyrus (rIFG), specifically in pars triangularis. Temporally, emotion and cognition processing overlapped during the transition from emotional to cognitive processing, with a stronger interaction in β-band power leading to worse behavioral performance. Despite functionally segregated subdivisions in rIFG, frequency-specific information flowed extensively within IFG and top-down to visual areas (V2, Precuneus) – explaining the behavioral interference effect. Thus, for the first time we here show the neural mechanisms of emotion-cognition interaction in space, time, frequency and information transfer with high temporal and spatial resolution, revealing a central role for β-band activity in rIFG. Our results support the idea that rIFG plays a broad role in both inhibitory control and emotional interference inhibition as it is a site of convergence in both processes. Furthermore, our results have potential clinical implications for understanding dysfunctional emotion-cognition interaction and emotional interference inhibition in psychiatric disor-ders, e.g. major depression and substance use disorder, in which patients have difficulties in regulating emotions and executing inhibitory control.
Summary statement When echolocating under demanding conditions e.g. noisy, narrow space, or cluttered environments, frugivorous bats adapt their call pattern by increasing the call rate within biosonar groups.
Abstract For orientation, echolocating bats emit biosonar calls and use echoes arising from call reflections. They often pattern their calls into groups which increases the rate of sensory feedback over time. Insectivorous bats emit call groups at a higher rate when orienting in cluttered compared to uncluttered environments. Frugivorous bats increase the rate of call group emission when they echolocate in noisy environments. Here, calls emitted by conspecifics potentially interfere with the bat’s biosonar signals and complicate the echolocation behavior. To minimize the information loss followed by signal interference, bats may profit from a temporally increased sensory acquisition rate, as it is the case for the call groups. In frugivorous bats, it remains unclear if call group emission represents an exclusive adaptation to avoid interference by signals from other bats or if it represents an adaptation that allows to orient under demanding environmental conditions. Here, we compared the emission pattern of the frugivorous bat Carollia perspicillata when the bats were flying in noisy versus silent, narrow versus wide or cluttered versus non-cluttered corridors. According to our results, the bats emitted larger call groups and they increased the call rate within the call groups when navigating in narrow, cluttered, or noisy environments. Thus, call group emission represents an adaptive behavior when the bats orient in complex environments.
Most mammals rely on the extraction of acoustic information from the environment in order to survive. However, the mechanisms that support sound representation in auditory neural networks involving sensory and association brain areas remain underexplored. In this study, we address the functional connectivity between an auditory region in frontal cortex (the frontal auditory field, FAF) and the auditory cortex (AC) in the bat Carollia perspicillata. The AC is a classic sensory area central for the processing of acoustic information. On the other hand, the FAF belongs to the frontal lobe, a brain region involved in the integration of sensory inputs, modulation of cognitive states, and in the coordination of behavioural outputs. The FAF-AC network was examined in terms of oscillatory coherence (local-field potentials, LFPs), and within an information theoretical framework linking FAF and AC spiking activity. We show that in the absence of acoustic stimulation, simultaneously recorded LFPs from FAF and AC are coherent in low frequencies (1-12 Hz). This “default” coupling was strongest in deep AC layers and was unaltered by acoustic stimulation. However, presenting auditory stimuli did trigger the emergence of coherent auditory-evoked gamma-band activity (>25 Hz) between the FAF and AC. In terms of spiking, our results suggest that FAF and AC engage in distinct coding strategies for representing artificial and natural sounds. Taken together, our findings shed light onto the neuronal coding strategies and functional coupling mechanisms that enable sound representation at the network level in the mammalian brain.
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.
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.
The ability to vocalize is ubiquitous in vertebrates, but neural networks leading to vocalization production remain poorly understood. Here we performed simultaneous, large scale, neuronal recordings in the frontal cortex and dorsal striatum (caudate nucleus) during the production of echolocation and non-echolocation calls in bats. This approach allows to assess the general aspects underlying vocalization production in mammals and the unique evolutionary adaptations of bat echolocation. Our findings show that distinct intra-areal brain rhythms in the beta (12-30 Hz) and gamma (30-80 Hz) bands of the local field potential can be used to predict the bats’ vocal output and that phase locking between spikes and field potentials occurs prior vocalization production. Moreover, the fronto-striatal network is differentially coupled in the theta-band during the production of echolocation and non-echolocation calls. Overall, our results present evidence for fronto-striatal network oscillations in motor action prediction in mammals.
Consistent individual differences in behaviour (animal personality) are widespread throughout the Animal Kingdom. This includes variation in risk-taking versus risk-averse behavioural tendencies. Variation in several personality dimensions is associated with distinct fitness consequences and thus, may become a target of natural and/or sexual selection. However, the link between animal personality and mate choice—as a major component of sexual selection—remains understudied. We asked (1) whether females and males of the livebearing fish Poecilia mexicana prefer risk-taking mating partners (directional mating preference), (2) or if their preferences are dependent on the choosing individual’s own personality type (assortative mating). We characterized each test subject for its risk-taking behaviour, assessed as the time to emerge from shelter and enter an unknown area. In dichotomous association preference tests, we offered two potential mating partners that differed in risk-taking behaviour but were matched for other phenotypic traits (body size, shape, and colouration). Females, but not males, exhibited a strong directional preference for risk-taking over risk-averse mating partners. At the same time, the strength of females’ preferences correlated positively with their own risk-taking scores. Our study is the first to demonstrate that a strong overall preference for risk-taking mating partners does not preclude effects of choosing individuals’ own personality type on (subtle) individual variation in mating preferences. More generally, two different preferences functions appear to interact to determine the outcome of individual mate choice decisions.
Cognition requires the dynamic modulation of effective connectivity, i.e. the modulation of the postsynaptic neuronal response to a given input. If postsynaptic neurons are rhythmically active, this might entail rhythmic gain modulation, such that inputs synchronized to phases of high gain benefit from enhanced effective connectivity. We show that visually induced gamma-band activity in awake macaque area V4 rhythmically modulates responses to unpredictable stimulus events. This modulation exceeded a simple additive superposition of a constant response onto ongoing gamma-rhythmic firing, demonstrating the modulation of multiplicative gain. Gamma phases leading to strongest neuronal responses also led to shortest behavioral reaction times, suggesting functional relevance of the effect. Furthermore, we find that constant optogenetic stimulation of anesthetized cat area 21a produces gamma-band activity entailing a similar gain modulation. As the gamma rhythm in area 21a did not spread backwards to area 17, this suggests that postsynaptic gamma is sufficient for gain modulation.
Synchronization has been implicated in neuronal communication, but causal evidence remains indirect. We used optogenetics to generate depolarizing currents in pyramidal neurons of cat visual cortex, emulating excitatory synaptic inputs under precise temporal control, while measuring spike output. Cortex transformed constant excitation into strong gamma-band synchronization, revealing the well-known cortical resonance. Increasing excitation with ramps increased the strength and frequency of synchronization. Slow, symmetric excitation profiles revealed hysteresis of power and frequency. Crucially, white-noise input sequences enabled causal analysis of network transmission, establishing that cortical resonance selectively transmits coherent input components. Models composed of recurrently coupled excitatory and inhibitory units uncovered a crucial role of feedback inhibition and suggest that hysteresis can arise through spike-frequency adaptation. The presented approach provides a powerful means to investigate the resonance properties of local circuits and probe how these properties transform input and shape transmission.
The gamma rhythm has been implicated in neuronal communication, but causal evidence remains indirect. We measured spike output of local neuronal networks and emulated their synaptic input through optogenetics. Opsins provide currents through somato-dendritic membranes, similar to synapses, yet under experimental control with high temporal precision. We expressed Channelrhodopsin-2 in excitatory neurons of cat visual cortex and recorded neuronal responses to light with different temporal characteristics. Sine waves of different frequencies entrained neuronal responses with a reliability that peaked for input frequencies in the gamma band. Crucially, we also presented white-noise sequences, because their temporal unpredictability enables analysis of causality. Neuronal spike output was caused specifically by the input’s gamma component. This gamma-specific transfer function is likely an emergent property of in-vivo networks with feedback inhibition. The method described here could reveal the transfer function between the input to any one and the output of any other neuronal group.
The plant family Brassicaceae includes some of the most studied hosts of plant microbiomes, targeting microbial diversity, community assembly rules, and effects on host performance. Compared to bacteria, eukaryotes in the brassicaceous microbiome remain understudied, especially under natural settings. Here, we assessed the impact of host identity and age on the assembly of fungal and oomycete root communities, using DNA metabarcoding of roots and associated soil of three annual co-habiting Brassicaceae collected at two time points. Our results showed that fungal communities are more diverse and structured than those of oomycetes. In both cases, plant identity and sampling time had little influence on community variation, whereas root/soil compartment had a strong effect by exerting control on the entry of soil microorganisms into the roots. The enrichment in roots of specific fungi suggests a specialization towards the asymptomatic colonization of plant tissues, which could be relevant to host’s fitness and health.
Wildfires are relatively rare in subarctic tundra ecosystems, but they can strongly change ecosystem properties. Short-term fire effects on subarctic tundra vegetation are well documented, but long-term vegetation recovery has been studied less. The frequency of tundra fires will increase with climate warming. Understanding the long-term effects of fire is necessary to predict future ecosystem changes.
We used a space-for-time approach to assess vegetation recovery after fire over more than four decades. We studied soil and vegetation patterns on three large fire scars (>44, 28 and 12 years old) in dry, lichen-dominated forest tundra in Western Siberia. On 60 plots, we determined soil temperature and permafrost thaw depth, sampled vegetation and measured plant functional traits. We assessed trends in NDVI to support the field-based results on vegetation recovery.
Soil temperature, permafrost thaw depth and total vegetation cover had recovered to pre-fire levels after >44 years, as well as total vegetation cover. In contrast, after >44 years, functional groups had not recovered to the pre-fire state. Burnt areas had lower lichen and higher bryophyte and shrub cover. The dominating shrub species, Betula nana, exhibited a higher vitality (higher specific leaf area and plant height) on burnt compared with control plots, suggesting a fire legacy effect in shrub growth. Our results confirm patterns of shrub encroachment after fire that were detected before in other parts of the Arctic and Subarctic. In the so far poorly studied Western Siberian forest tundra we demonstrate for the first time, long-term fire-legacies on the functional composition of relatively dry shrub- and lichen-dominated vegetation.
Seed harvesting from wild plant populations is key for ecological restoration, but may threaten the persistence of source populations. Consequently, several countries have set guidelines limiting the proportions of harvestable seeds. Here, we use high-resolution data from 298 plant species to model the demographic consequences of seed harvesting. We find that the current guidelines only protect some species, but are insufficient or overly restrictive for others. We show that the maximum possible fraction of seed harvesting is strongly associated with harvesting frequency and generation time of the target species, ranging from 100% in long-lived species to <1% in the most annuals. Our results provide quantitative basis to guide seed harvesting legislation based on species’ generation time and harvesting regime.
Animals sense ambient temperature so that they can adjust their behavior to the environment; they avoid noxious heat and coldness and stay within a survivable temperature range. C. elegans can sense temperature, memorize past cultivation temperature and navigate towards preferable temperature, for which a thermosensory neuron, AFD, is essential. AFD responds to temperature increase from the past cultivation temperature by increasing intracellular Ca2+ level. We aimed to reveal how AFD encodes and memorizes the information of temperature. Although cGMP synthesis is crucial for thermosensation by AFD, whether and how cGMP level temporally fluctuates in AFD remained elusive. We therefore monitored cGMP level in AFD and found that cGMP dynamically responded to temperature change in a manner dependent on past cultivation temperature. Given that cGMP dynamics is supposed to be upstream of Ca2+ dynamics, our results suggest that AFD’s memory is formed by simpler molecular mechanisms than previously expected from the Ca2+ dynamics. Moreover, we analyzed how guanylyl cyclases and phosphodiesterases, which synthesize and degrade cGMP, respectively, contributed to cGMP and Ca2+ dynamics and thermotaxis behavior.
Signal transfer of visual stimuli to V4 occurs in gamma-rhythmic, pulsed information packages
(2020)
Summary Selective visual attention allows the brain to focus on behaviorally relevant information while ignoring irrelevant signals. As a possible mechanism, routing by synchronization was proposed: neural populations sending attended signals align their gamma-rhythmic activities with receiving populations, such that spikes from the senders arrive at excitability peaks of the receivers, enhancing signal transfer. Conversely, the non-attended signals arrive unaligned to the receiver’s oscillation, reducing signal transfer. Therefore, visual signals should be transferred through periodically pulsed information packages, resulting in a modulation of the stimulus content within the receiver’s activity by its gamma phase and amplitude. To test this prediction, we quantified gamma phase-specific stimulus content within neural activity from area V4 of macaques performing a visual attention task. For the attended stimulus we find enhanced stimulus content reaching its maximum near excitability peaks, with effect magnitude increasing with oscillation amplitude, establishing a functional link between selective processing and gamma activity.
No disease modifying therapy is currently available for Parkinson’s disease (PD), the second most common neurodegenerative disease. The long non-motor prodromal phase of PD is a window of opportunity for early detection and intervention. However, we lack the pathophysiological understanding to develop selective biomarkers and interventions. By developing a mutant α-synuclein selective-overexpression mouse model of prodromal PD, we identified a cell-autonomous selective Kv4 channelopathy in dorsal motor nucleus of the vagus (DMV) neurons. This functional remodeling of intact DMV neurons leads to impaired pacemaker function in vitro and in vivo, which in turn reduces gastrointestinal motility which is a common, very early symptom of prodromal PD. We show for the first time a causal chain of events from α-synuclein via a biophysical dysfunction of specific neuronal populations to a clinically relevant prodromal symptom. These findings can facilitate the rational design of clinical biomarkers to identify people at risk for PD.
Substantia nigra dopamine (SN DA) neurons are progressively lost in Parkinson disease (PD). While the molecular and cellular mechanisms of their differential vulnerability and degeneration have been extensively studied, we still know very little about potential functional adaptations of those SN DA neurons that – at least for some time – manage to survive during earlier stages of PD. We utilized a partial lesion 6-OHDA mouse model to characterize initial electrophysiological impairments and chronic adaptations of surviving identified SN DA neurons, both in vivo and in vitro. Early after lesion (3 weeks), we detected a selective loss of in vivo burst firing in surviving SN DA neurons, which was accompanied by in vitro pacemaker instability. In contrast, late after lesion (>2 months), in vivo firing properties of surviving SN DA neurons had recovered in the presence of 2-fold accelerated pacemaking in vitro. Finally, we show that this chronic cell-autonomous adaptation in surviving SN DA neurons was mediated by Kv4.3 channel downregulation. Our study demonstrates substantial homeostatic plasticity of surviving SN DA neurons after a single-hit non-progressive lesion, which might contribute to the phenotype of initially surviving SN DA neurons in PD.
Parkinson disease (PD), one of the most common neurodegenerative disorder, is believed to be driven by toxic α-synuclein aggregates eventually resulting in selective loss of vulnerable neuron populations, prominent among them, nigrostriatal dopamine (DA) neurons in the lateral substantia nigra (l-SN). How α-synuclein aggregates initiate a pathophysiological cascade selectively in vulnerable neurons is still unclear. Here, we show that the exposure to low nanomolar concentrations of α-synuclein aggregates (i.e. fibrils) but not its monomeric forms acutely and selectively disrupted the electrical pacemaker function of the DA subpopulation most vulnerable in PD. This implies that only dorsolateral striatum projecting l-SN DA neurons were electrically silenced by α-synuclein aggregates, while the activity of neither neighboring DA neurons in medial SN projecting to dorsomedial striatum nor mesolimbic DA neurons in the ventral tegmental area (VTA) were affected. Moreover, we demonstrate functional K-ATP channels comprised of Kir6.2 subunit in DA neurons to be necessary to mediate this acute pacemaker disruption by α-synuclein aggregates. Our study thus identifies a molecularly defined target that quickly translates the presence of α-synuclein aggregates into an immediate impairment of essential neuronal function. This constitutes a novel candidate process how a protein-aggregation-driven sequence in PD is initiated that might eventually lead to selective neurodegeneration.
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.
Background Transposable elements (TEs) are an important source of genome plasticity across the tree of life. Accumulating evidence suggests that TEs may not be randomly distributed in the genome. Drift and natural selection are important forces shaping TE distribution and accumulation, acting directly on the TE element or indirectly on the host species. Fungi, with their multifaceted phenotypic diversity and relatively small genome size, are ideal models to study the role of TEs in genome evolution and their impact on the host’s ecological and life history traits. Here we present an account of all TEs found in a high-quality reference genome of the lichen-forming fungus Umbilicaria pustulata, a macrolichen species comprising two climatic ecotypes: Mediterranean and cold-temperate. We trace the occurrence of the newly identified TEs in populations along three replicated elevation gradients using a Pool-Seq approach, to identify TE insertions of potential adaptive significance.
Results We found that TEs cover 21.26 % of the 32.9 Mbp genome, with LTR Gypsy and Copia clades being the most common TEs. Out of a total of 182 TE copies we identified 28 insertions displaying consistent insertion frequency differences between the two host ecotypes across the elevation gradients. Most of the highly differentiated insertions were located near genes, indicating a putative function.
Conclusions This pioneering study into the content and climate niche-specific distribution of TEs in a lichen-forming fungus contributes to understanding the roles of TEs in fungal evolution. Particularly, it may serve as a foundation for assessing the impact of TE dynamics on fungal adaptation to the abiotic environment, and the impact of TE activity on the evolution and maintenance of a symbiotic lifestyle.
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.
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.
Afterimages result from a prolonged exposure to still visual stimuli. They are best detectable when viewed against uniform backgrounds and can persist for multiple seconds. Consequently, the dynamics of afterimages appears to be slow by their very nature. To the contrary, we report here that about 50% of an afterimage intensity can be erased rapidly—within less than a second. The prerequisite is that subjects view a rich visual content to erase the afterimage; fast erasure of afterimages does not occur if subjects view a blank screen. Moreover, we find evidence that fast removal of afterimages is a skill learned with practice as our subjects were always more effective in cleaning up afterimages in later parts of the experiment. These results can be explained by a tri-level hierarchy of adaptive mechanisms, as has been proposed by the theory of practopoiesis.
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.
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.
Cross-frequency coupling (CFC) has been proposed to coordinate neural dynamics across spatial and temporal scales. Despite its potential relevance for understanding healthy and pathological brain function, the standard CFC analysis and physiological interpretation come with fundamental problems. For example, apparent CFC can appear because of spectral correlations due to common non-stationarities that may arise in the total absence of interactions between neural frequency components. To provide a road map towards an improved mechanistic understanding of CFC, we organize the available and potential novel statistical/modeling approaches according to their biophysical interpretability. While we do not provide solutions for all the problems described, we provide a list of practical recommendations to avoid common errors and to enhance the interpretability of CFC analysis.
Untangling the cell immune response dynamic for severe and critical cases of SARS-CoV-2 infection
(2021)
COVID-19 is a global pandemic leading high death tolls worldwide day by day. Clinical evidence suggests that COVID-19 patients can be classified as non-severe, severe and critical cases. In particular, studies have highlighted the relationship between the lymphopenia and the severity of the illness, where CD8+ T cells have the lowest levels in critical cases. In this work, we aim to elucidate the key parameters that define the course of the disease deviating from severe to critical case. To this end, several mathematical models are proposed to represent the dynamic of the immune response in patients with SARS-CoV-2 infection. The best model had a good fit to reported experimental data, and in accordance with values found in the literature. Our results suggest that a rapid proliferation of CD8+ T cells is decisive in the severity of the disease.
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.
Abstract
Co-infections by multiple pathogens have important implications in many aspects of health, epidemiology and evolution. However, how to disentangle the contributing factors of the immune response when two infections take place at the same time is largely unexplored. Using data sets of the immune response during influenza-pneumococcal co-infection in mice, we employ here topological data analysis to simplify and visualise high dimensional data sets.
We identified persistent shapes of the simplicial complexes of the data in the three infection scenarios: single viral infection, single bacterial infection, and co-infection. The immune response was found to be distinct for each of the infection scenarios and we uncovered that the immune response during the co-infection has three phases and two transition points. During the first phase, its dynamics is inherited from its response to the primary (viral) infection. The immune response has an early (few hours post co-infection) and then modulates its response to finally react against the secondary (bacterial) infection. Between 18 to 26 hours post co-infection the nature of the immune response changes again and does no longer resembles either of the single infection scenarios.
Author summary
The mapper algorithm is a topological data analysis technique used for the qualitative analysis, simplification and visualisation of high dimensional data sets. It generates a low-dimensional image that captures topological and geometric information of the data set in high dimensional space, which can highlight groups of data points of interest and can guide further analysis and quantification.
To understand how the immune system evolves during the co-infection between viruses and bacteria, and the role of specific cytokines as contributing factors for these severe infections, we use Topological Data Analysis (TDA) along with an extensive semi-unsupervised parameter value grid search, and k-nearest neighbour analysis.
We find persistent shapes of the data in the three infection scenarios, single viral and bacterial infections and co-infection. The immune response is shown to be distinct for each of the infections scenarios and we uncover that the immune response during the co-infection has three phases and two transition points, a previously unknown property regarding the dynamics of the immune response during co-infection.
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.
During infection the SARS-CoV-2 virus fuses its viral envelope with cellular membranes of its human host. Initial contact with the host cell and membrane fusion are both mediated by the viral spike (S) protein. Proteolytic cleavage of S at the S2′ site exposes its 40 amino acid long fusion peptide (FP). Binding of the FP to the host membrane anchors the S2 domain of S in both the viral and the host membrane. The reorganization of S2 then pulls the two membranes together. Here we use molecular dynamics (MD) simulations to study the two core functions of the SARS-CoV-2 FP: to attach quickly to cellular membranes and to form an anchor strong enough to withstand the mechanical force during membrane fusion. In eight 10 μs-long MD simulations of FP in proximity to endosomal and plasma membranes, we find that FP binds spontaneously to the membranes and that binding proceeds predominantly by insertion of two short amphipathic helices into the membrane interface. Connected via a flexible linker, the two helices can bind the membrane independently, yet binding of one promotes the binding of the other by tethering it close to the target membrane. By simulating mechanical pulling forces acting on the C-terminus of the FP we then show that the bound FP can bear forces up to 250 pN before detaching from the membrane. This detachment force is more than ten-fold higher than an estimate of the force required to pull host and viral membranes together for fusion. We identify a fully conserved disulfide bridge in the FP as a major factor for the high mechanical stability of the FP membrane anchor. We conclude, first, that the sequential binding of two short amphipathic helices allows the SARS-CoV-2 FP to insert quickly into the target membrane, before the virion is swept away after shedding the S1 domain connecting it to the host cell receptor. Second, we conclude that the double attachment and the conserved disulfide bridge establish the strong anchoring required for subsequent membrane fusion. Multiple distinct membrane-anchoring elements ensure high avidity and high mechanical strength of FP-membrane binding.
Learning in the eyes: specific changes in gaze patterns track explicit and implicit visual learning
(2020)
What is the link between eye movements and sensory learning? Although some theories have argued for a permanent and automatic interaction between what we know and where we look, which continuously modulates human information- gathering behavior during both implicit and explicit learning, there exist surprisingly little evidence supporting such an ongoing interaction. We used a pure form of implicit learning called visual statistical learning and manipulated the explicitness of the task to explore how learning and eye movements interact. During both implicit exploration and explicit visual learning of unknown composite visual scenes, eye movement patterns systematically changed in accordance with the underlying statistical structure of the scenes. Moreover, the degree of change was directly correlated with the amount of knowledge the observers acquired. Our results provide the first evidence for an ongoing and specific interaction between hitherto accumulated knowledge and eye movements during both implicit and explicit learning.
VASP is a member of the Enabled/VASP protein family that is involved in cortical actin dynamics and may also contribute to the formation of gap junctions. In vessels, gap junctional coupling allows the transfer of signals along the vessel wall and coordinates vascular behavior. Moreover, VASP is reportedly a mediator of NO-induced inhibition of platelet aggregation. Therefore, we hypothesized that VASP exerts also important physiologic functions in arterioles. We examined the spread of vasodilations enabled by gap junctional coupling in endothelial cells as well as NO-induced arteriolar dilations in VASP-deficient mice by intravital microscopy of the microcirculation in a skeletal muscle in anesthetized mice. Conducted dilations were initiated by brief, locally confined stimulation of the arterioles with acetylcholine. The maximal diameters of the arterioles under study ranged from 30 to 40 μm. Brief stimulation with acetylcholine induced a short dilation at the local site that was also observed at remote, upstream sites without an attenuation of the amplitude up to a distance of 1.2 mm in control animals (wild-type). In contrast, remote dilations were reduced in VASP-deficient mice despite a similar local dilation indicating an impairment of conducted dilations. Superfusion of NOdonors induced a concentration-dependent dilation in wild-type mice. However, these dilations were slightly reduced in VASP-deficient animals. In contrast, dilations induced by the endothelial stimulator acetylcholine were fully preserved in VASP-deficient mice. In summary, this study suggests that VASP exerts critical functions in arteriolar diameter control. It is crucial for the conduction of dilator signals along the endothelial cell layer. The impairment possibly reflects a perturbed formation of gap junctions in the endothelial cell membrane. VASP also participates in the full dilatory potential of NOdonors although the effect of its deficiency is only subtle. In contrast, VASP is not required for dilations initiated by endothelial stimulation which are mediated in the murine microcirculation by an EDH-mechanism.
The SARS-CoV-2 pandemic has challenged researchers at a global scale. The scientific community’s massive response has resulted in a flood of experiments, analyses, hypotheses, and publications, especially in the field of drug repurposing. However, many of the proposed therapeutic compounds obtained from SARS-CoV-2 specific assays are not in agreement and thus demonstrate the need for a singular source of COVID-19 related information from which a rational selection of drug repurposing candidates can be made. In this paper, we present the COVID-19 PHARMACOME, a comprehensive drug-target-mechanism graph generated from a compilation of 10 separate disease maps and sources of experimental data focused on SARS-CoV-2 / COVID-19 pathophysiology. By applying our systematic approach, we were able to predict the synergistic effect of specific drug pairs, such as Remdesivir and Thioguanosine or Nelfinavir and Raloxifene, on SARS-CoV-2 infection. Experimental validation of our results demonstrate that our graph can be used to not only explore the involved mechanistic pathways, but also to identify novel combinations of drug repurposing candidates.
Plant communities provide floral resource-landscapes for pollinators. Yet, it is insufficiently understood how these landscapes shape pollinator-mediated interactions among multiple plant species. Here, we study how pollinators and the seed set of plants respond to the distribution of a floral resource (nectar sugar) in space and across plant species, inflorescences and flowering phenologies. In a global biodiversity hotspot, we quantified floral resource-landscapes on 27 sites of 4 ha comprising 127,993 shrubs of 19 species. Visitation rates of key bird pollinators strongly depended on the phenology of site-scale resource amounts. Seed set of focal plants increased with resources of conspecific neighbours and with site-scale resources, notably with heterospecific resources of lower quality (less sugar per inflorescence). Floral resources are thus a common currency determining how multiple plant species interact via pollinators. These interactions may alter conditions for species coexistence in plant communities and cause community-level Allee effects that promote extinction cascades.
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.
How much data do we need? Lower bounds of brain activation states to predict human cognitive ability
(2022)
Human functional brain connectivity can be temporally decomposed into states of high and low cofluctuation, defined as coactivation of brain regions over time. Despite their low frequency of occurrence, states of particularly high cofluctuation have been shown to reflect fundamentals of intrinsic functional network architecture (derived from resting-state fMRI) and to be highly subject-specific. However, it is currently unclear whether such network-defining states of high cofluctuation also contribute to individual variations in cognitive abilities – which strongly rely on the interactions among distributed brain regions. By introducing CMEP, an eigenvector-based prediction framework, we show that functional connectivity estimates from as few as 20 temporally separated time frames (< 3% of a 10 min resting-state fMRI scan) are significantly predictive of individual differences in intelligence (N = 281, p < .001). In contrast and against previous expectations, individual’s network-defining time frames of particularly high cofluctuation do not achieve significant prediction of intelligence. Multiple functional brain networks contribute to the prediction, and all results replicate in an independent sample (N = 831). Our results suggest that although fundamentals of person-specific functional connectomes can be derived from few time frames of highest brain connectivity, temporally distributed information is necessary to extract information about cognitive abilities from functional connectivity time series. This information, however, is not restricted to specific connectivity states, like network-defining high-cofluctuation states, but rather reflected across the entire length of the brain connectivity time series.
Probing the association between resting state brain network dynamics and psychological resilience
(2021)
Abstract
This study aimed at replicating a previously reported negative correlation between node flexibility and psychological resilience, i.e., the ability to retain mental health in the face of stress and adversity. To this end, we used multiband resting-state BOLD fMRI (TR = .675 sec) from 52 participants who had filled out three psychological questionnaires assessing resilience. Time-resolved functional connectivity was calculated by performing a sliding window approach on averaged time series parcellated according to different established atlases. Multilayer modularity detection was performed to track network reconfigurations over time and node flexibility was calculated as the number of times a node changes community assignment. In addition, node promiscuity (the fraction of communities a node participates in) and node degree (as proxy for time-varying connectivity) were calculated to extend previous work. We found no substantial correlations between resilience and node flexibility. We observed a small number of correlations between the two other brain measures and resilience scores, that were however very inconsistently distributed across brain measures, differences in temporal sampling, and parcellation schemes. This heterogeneity calls into question the existence of previously postulated associations between resilience and brain network flexibility and highlights how results may be influenced by specific analysis choices.
Author Summary We tested the replicability and generalizability of a previously proposed negative association between dynamic brain network reconfigurations derived from multilayer modularity detection (node flexibility) and psychological resilience. Using multiband resting-state BOLD fMRI data and exploring several parcellation schemes, sliding window approaches, and temporal resolutions of the data, we could not replicate previously reported findings regarding the association between node flexibility and resilience. By extending this work to other measures of brain dynamics (node promiscuity, degree) we observe a rather inconsistent pattern of correlations with resilience, that strongly varies across analysis choices. We conclude that further research is needed to understand the network neuroscience basis of mental health and discuss several reasons that may account for the variability in results.
Word familiarity and predictive context facilitate visual word processing, leading to faster recognition times and reduced neuronal responses. Previously, models with and without top-down connections, including lexical-semantic, pre-lexical (e.g., orthographic/ phonological), and visual processing levels were successful in accounting for these facilitation effects. Here we systematically assessed context-based facilitation with a repetition priming task and explicitly dissociated pre-lexical and lexical processing levels using a pseudoword familiarization procedure. Experiment 1 investigated the temporal dynamics of neuronal facilitation effects with magnetoencephalography (MEG; N=38 human participants) while Experiment 2 assessed behavioral facilitation effects (N=24 human participants). Across all stimulus conditions, MEG demonstrated context-based facilitation across multiple time windows starting at 100 ms, in occipital brain areas. This finding indicates context based-facilitation at an early visual processing level. In both experiments, we furthermore found an interaction of context and lexical familiarity, such that stimuli with associated meaning showed the strongest context-dependent facilitation in brain activation and behavior. Using MEG, this facilitation effect could be localized to the left anterior temporal lobe at around 400 ms, indicating within-level (i.e., exclusively lexical-semantic) facilitation but no top-down effects on earlier processing stages. Increased pre-lexical familiarity (in pseudowords familiarized utilizing training) did not enhance or reduce context effects significantly. We conclude that context based-facilitation is achieved within visual and lexical processing levels. Finally, by testing alternative hypotheses derived from mechanistic accounts of repetition suppression, we suggest that the facilitatory context effects found here are implemented using a predictive coding mechanism.
To characterize the left-ventral occipito-temporal cortex (lvOT) role during reading in a quantitatively explicit and testable manner, we propose the lexical categorization model (LCM). The LCM assumes that lvOT optimizes linguistic processing by allowing fast meaning access when words are familiar and filter out orthographic strings without meaning. The LCM successfully simulates benchmark results from functional brain imaging. Empirically, using functional magnetic resonance imaging, we demonstrate that quantitative LCM simulations predict lvOT activation across three studies better than alternative models. Besides, we found that word-likeness, which is assumed as input to LCM, is represented posterior to lvOT. In contrast, a dichotomous word/non-word contrast, which is assumed as the LCM’s output, could be localized to upstream frontal brain regions. Finally, we found that training lexical categorization results in more efficient reading. Thus, we propose a ventral-visual-stream processing framework for reading involving word-likeness extraction followed by lexical categorization, before meaning extraction.