title: Measurement of Klebsiella Intestinal Colonization Density to Assess Infection Risk

abstract:
BackgroundKlebsiella pneumoniae and closely related species K. variicola and K. quasipneumoniae are common causes of healthcare-associated infections, and patients frequently become infected with their intestinal colonizing strain. To assess the association between Klebsiella colonization density and subsequent infections, a case-control study was performed.

MethodsA multiplex qPCR assay was developed and validated to quantify Klebsiella (K. pneumoniae, K. variicola, and K. quasipneumoniae combined) relative to total bacterial DNA copies in rectal swabs. Cases of Klebsiella infection were identified based on clinical definitions and having a clinical culture isolate and preceding or co-incident colonization isolate with the same wzi capsular sequence type. Controls were colonized patients without subsequent infection and were matched 2:1 to cases based on age, sex, and rectal swab collection date. Quantitative PCR (qPCR) from rectal swab samples was used to measure the association between relative abundance (RA) of Klebsiella and subsequent infections.

ResultsKlebsiella RA by qPCR highly correlated with 16S sequencing ({rho}=0.79; P <.001). The median Klebsiella RA in the study group was 2.6% (interquartile range (IQR) 0.1-22.5, n=238), and was higher in cases (15.7%, IQR 0.93-52.6%, n=83) than controls (1.01%, IQR 0.02-12.8%; n=155; P <0.0001). After adjusting for multiple clinical covariates using inverse probability of treatment weighting, subjects with a Klebsiella RA > 22% had a 2.87-fold (1.64-5.03, P =0.0003) increased odds of infection compared to those with lower colonization density levels.

ConclusionsMeasurement of colonization density by qPCR could represent a novel approach to identify hospitalized patients at risk for Klebsiella infection.

ImportanceColonization by bacterial pathogens often precedes infection, and offers a window of opportunity to prevent these infections. Klebsiella colonization is significantly and reproducibly associated with subsequent infection, however factors that enhance or mitigate this risk in individual patients are unclear. This study developed an assay to measure the density of Klebsiella colonization, relative to total fecal bacteria, in rectal swabs from hospitalized patients. Applying this assay to 238 colonized patients, high Klebsiella density defined as >22% of total bacteria, was significantly associated with subsequent infection. Based on widely available polymerase chain reaction (PCR) technology, this type of assay could be deployed in clinical laboratories to identify patients at increased risk of Klebsiella infections. As novel therapeutics are developed to eliminate pathogens from the gut microbiome, a rapid Klebsiella colonization density assay could identify patients who would benefit from this type of infection prevention interventions.


title: Sec17/Sec18 can support membrane fusion without help from completion of SNARE zippering

abstract: Membrane fusion requires R-, Qa-, Qb-, and Qc-family SNAREs that zipper into RQaQbQc coiled coils, driven by the sequestration of apolar amino acids. Zippering has been thought to provide all the force driving fusion. Sec17/SNAP can form an oligomeric assembly with SNAREs with the Sec17 C-terminus bound to Sec18/NSF, the central region bound to SNAREs, and a crucial apolar loop near the N-terminus poised to insert into membranes. Though Sec17 aids zippering, we now report that Sec17 and Sec18 will drive robust fusion without requiring zippering completion. Zippering-driven fusion is blocked by deleting the C-terminal quarter of any Q-SNARE domain or by replacing the apolar amino acids of the Qa-SNARE which face the center of the 4-SNARE coiled coils with polar residues. These blocks to fusion, singly or combined, are bypassed by Sec17 and Sec18, and SNARE-dependent fusion is restored without help from completing zippering.

category: biochemistry


title: Crowdsourced Identification of Multi-Target Kinase Inhibitors for RET- and TAU- Based Disease: the Multi-Targeting Drug DREAM Challenge

abstract:
A continuing challenge in modern medicine is the identification of safer and more efficacious drugs. Precision therapeutics, which have one molecular target, have been long promised to be safer and more effective than traditional therapies. This approach has proven to be challenging for multiple reasons including lack of efficacy, rapidly acquired drug resistance, and narrow patient eligibility criteria. An alternative approach is the development of drugs that address the overall disease network by targeting multiple biological targets ( polypharmacology). Rational development of these molecules will require improved methods for predicting single chemical structures that target multiple drug targets. To address this need, we developed the Multi-Targeting Drug DREAM Challenge, in which we challenged participants to predict single chemical entities that target pro-targets but avoid anti-targets for two unrelated diseases: RET-based tumors and a common form of inherited Tauopathy. Here, we report the results of this DREAM Challenge and the development of two neural network-based machine learning approaches that were applied to the challenge of rational polypharmacology. Together, these platforms provide a potentially useful first step towards developing lead therapeutic compounds that address disease complexity through rational polypharmacology.

Author SummaryMany modern drugs are developed with the goal of modulating a single cellular pathway or target. However, many drugs are, in fact,  dirty; they target multiple cellular pathways or targets. This phenomenon is known as multi-targeting or polypharmacology. While some strive to develop  cleaner therapeutics that eliminate secondary targets, recent work has shown that multi-targeting therapeutics have key advantages for a variety of diseases. However, while multi-targeting drugs that affect a precisely-defined profile of targets may be more effective, it is difficult to computationally predict which molecules have desirable target profiles. Here, we report the results of a competitive crowdsourcing project (the Multi-Targeting Drug DREAM Challenge), where we challenged participants to predict chemicals that have desired target profiles for cancer and neurodegenerative disease.

category: pharmacology and toxicology


title: Crowdsourced identification of multi-target kinase inhibitors for RET- and TAU- based disease: the Multi-Targeting Drug DREAM Challenge


title: Fossil data do not support a long Pre-Cretaceous history of flowering plants

abstract: The origin of angiosperms is a classic macroevolutionary problem, because of their rapid rise in the Early Cretaceous fossil record, beginning about 139 Ma ago, and the conflict this creates with older crown-group ages based on molecular clock dating1. Silvestro et al.2 use a novel methodology to model past angiosperm diversity based on a Bayesian Brownian Bridge model of fossil finds assigned to extant families, concluding that a Cretaceous origin is vanishingly unlikely. However, their results strongly conflict with the known temporal distribution of angiosperm fossils, and, while we agree that statistical analysis aids interpretation of the fossil record, here we show the conclusions of Silvestro et al.2 are unsound.

category: paleontology


title: RaptGen: A variational autoencoder with profile hidden Markov model for generative aptamer discovery

abstract: Nucleic acid aptamers are generated by an in vitro molecular evolution method known as systematic evolution of ligands by exponential enrichment (SELEX). A variety of candidates is limited by actual sequencing data from an experiment. Here, we developed RaptGen, which is a variational autoencoder for in silico aptamer generation. RaptGen exploits a profile hidden Markov model decoder to represent motif sequences effectively. We showed that RaptGen embedded simulation sequence data into low-dimension latent space dependent on motif information. We also performed sequence embedding using two independent SELEX datasets. RaptGen successfully generated aptamers from the latent space even though they were not included in high-throughput sequencing. RaptGen could also generate a truncated aptamer with a short learning model. We demonstrated that RaptGen could be applied to activity-guided aptamer generation according to Bayesian optimization. We concluded that a generative method by RaptGen and latent representation are useful for aptamer discovery. Codes are available at https://github.com/hmdlab/raptgen.

category: bioinformatics


title: Early Dynamics of Chronic Myeloid Leukemia on Nilotinib Predicts Deep Molecular Response

abstract:
Chronic myeloid leukemia (CML) is a myeloproliferative disorder caused by the BCR-ABL1 tyrosine kinase1,2. ABL1-selective tyrosine kinase inhibitors (TKIs) including nilotinib have dramatically improved the prognosis of patients with CML3-7. The ultimate goal of CML treatment is likely to be TKI-free maintenance of deep molecular response (DMR), which is defined by quantitative measurement of BCR-ABL1 transcripts on the international scale (IS)8, and durable DMR is a prerequisite to reach this goal9. Thus, an algorithm to enable the early prediction of DMR achievement on TKI therapy is quite valuable. Here, we show that our mathematical framework based on a clinical trial dataset10 can accurately predict the response to frontline nilotinib. We found that our simple dynamical model can predict nilotinib response by using two common laboratory findings (observation values): IS11,12 and white blood cell (WBC) count. Furthermore, our proposed method identified patients who failed to achieve DMR with high fidelity according to the data collected only at three initial time points during nilotinib therapy. Since our model relies on the general properties of TKI response, our framework would be applicable to CML patients who receive frontline nilotinib or other TKIs in clinical practice.

Significance StatementChronic myeloid leukemia (CML) is a myeloproliferative disorder caused by the BCR-ABL1 tyrosine kinase. The goal of this treatment is the sequential achievement of deep molecular response (DMR). Tyrosine kinase inhibitors (TKIs) are effective in the reduction because they inhibit CML cell proliferation. However, because of individual differences in the TKI efficacy, some patients are unable to achieve DMR, so that early prediction of DMR reachability is necessary for personalized medicine. By combining time series analysis and mathematical modeling, we developed a mathematical method that accurately predicts patients who do not achieve DMR in the early stages of TKI administration. Our prediction method gives a basis of effective personalized treatments for CML patients.

category: systems biology


title: Gene loss during the transition to multicellularity

abstract:
Multicellular evolution is a major transition associated with momentous diversification of multiple lineages and increased developmental complexity. The volvocine algae comprise a valuable system for the study of this transition, as they span from unicellular to undifferentiated and differentiated multicellular morphologies despite their genomes being highly similar, suggesting multicellular evolution requires few genetic changes to undergo dramatic shifts in developmental complexity. Here, the evolutionary dynamics of six volvocine genomes were examined, where a gradual loss of genes was observed in parallel to the co-option of a few key genes. Protein complexes in the six species exhibited a high degree of novel interactions, suggesting that gene loss plays a role in evolutionary novelty. This finding was supported by gene network modeling, where gene loss outpaces gene gain in generating novel stable network states. These results suggest developmental complexity may be driven by gene loss, and not just by gene gain.

Significance StatementIncreased developmental complexity is thought to evolve mainly through genetic innovation and co-option. Comparison of the genomes of the closely related volvocines revealed that despite large changes in developmental complexity their genomes are exhibiting significant gene loss. Further, a burst of gene loss is shown to occur at the transition to undifferentiated multicellularity by gene inactivation and decay. This likely results in changes to protein-protein interactions within the cell, suggesting evolutionary novelty does not always require gene gain. Using empirical and modelling approaches we demonstrate gene loss can more easily produce protein network novelty than gene gain. We propose that gene loss can be a driver of biological innovation, allowing for reconfiguration of gene networks and differential use of existing functional repertoires.

category: evolutionary biology


title: Gene loss during a transition to multicellularity


title: A test statistic to quantify treelikeness in phylogenetics

abstract: Most phylogenetic analyses assume that the evolutionary history of an alignment (either that of a single locus, or of multiple concatenated loci) can be described by a single bifurcating tree, the so-called the treelikeness assumption. Treelikeness can be violated by biological events such as recombination, introgression, or incomplete lineage sorting, and by systematic errors in phylogenetic analyses. The incorrect assumption of treelikeness may then mislead phylogenetic inferences. To quantify and test for treelikeness in alignments, we develop a test statistic which we call the tree proportion. This statistic quantifies the proportion of the edge weights in a phylogenetic network that are represented in a bifurcating phylogenetic tree of the same alignment. We extend this statistic to a statistical test of treelikeness using a parametric bootstrap. We use extensive simulations to compare tree proportion to a range of related approaches. We show that tree proportion successfully identifies non-treelikeness in a wide range of simulation scenarios, and discuss its strengths and weaknesses compared to other approaches. The power of the tree-proportion test to reject non-treelike alignments can be lower than some other approaches, but these approaches tend to be limited in their scope and/or the ease with which they can be interpreted. Our recommendation is to test treelikeness of sequence alignments with both tree proportion and mosaic methods such as 3Seq. The scripts necessary to replicate this study are available at https://github.com/caitlinch/treelikeness


title: Pak1 kinase regulates ribonucleoprotein granules through Sts5 for polarized growth and the glucose starvation response

abstract: Fission yeast cells maintain a rod shape due to conserved signaling pathways that organize the cytoskeleton for polarized growth. We discovered a mechanism linking the conserved protein kinase Pak1 with cell shape through the RNA-binding protein Sts5. Pak1 (also called Shk1 and Orb2) prevents Sts5 association with P bodies by directly phosphorylating its intrinsically disordered region (IDR). Pak1 and the cell polarity kinase Orb6 both phosphorylate the Sts5 IDR but at distinct residues. Mutations preventing phosphorylation in the Sts5 IDR cause increased P body formation and defects in cell shape and polarity. Unexpectedly, when cells encounter glucose starvation, PKA signaling triggers Pak1 recruitment to stress granules with Sts5. Through retargeting experiments, we reveal that Pak1 localizes to stress granules to promote rapid dissolution of Sts5 upon glucose addition. Our work reveals a new role for Pak1 in regulating cell shape through ribonucleoprotein granules during normal and stressed growth conditions.

category: cell biology


title: Pak1 kinase controls cell shape through ribonucleoprotein granules


title: Multidimensional Single-Nuclei RNA-Seq Reconstruction of Adipose Tissue Reveals Adipocyte Plasticity Underlying Thermogenic Response

abstract:
Adipose tissue has been classified based on its morphology and function as white, brown, or beige / brite. It plays an essential role as a regulator of systemic metabolism through paracrine and endocrine signals. Recently, multiple adipocyte subtypes have been revealed using RNA sequencing technology, going beyond simply defined morphology but by their cellular origin, adaptation to metabolic stress, and plasticity. Here, we performed an in-depth analysis of publicly available single-nuclei RNAseq from adipose tissue and utilized a workflow template to characterize adipocyte plasticity, heterogeneity, and secretome profiles. The reanalyzed dataset led to the identification of different subtypes of adipocytes including three subpopulations of thermogenic adipocytes and provided a characterization of distinct transcriptional profiles along the adipocyte trajectory under thermogenic challenges. This study provides a useful resource for further investigations regarding mechanisms related to adipocyte plasticity and trans-differentiation.

HighlightsMultidimensional transcriptome analysis at single-nucleus resolution recovers nuclei of cell types in adipose tissue

Adaptative thermogenic response results in 3 distinct mature adipose cell types

Single-nuclei transcriptomic-based secretome analysis reveals adipose cell-type-specific genes

The in vivo trajectory of adipocyte plasticity for thermogenic response reveals sets of trans-differentiation genes

Graphic Abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=200 SRC="FIGDIR/small/431320v2_ufig1.gif" ALT="Figure 1">
View larger version (37K):
org.highwire.dtl.DTLVardef@291767org.highwire.dtl.DTLVardef@1bfa3b1org.highwire.dtl.DTLVardef@93e051org.highwire.dtl.DTLVardef@6bef5e_HPS_FORMAT_FIGEXP  M_FIG C_FIG


title: Phase transitions may explain why SARS-CoV-2 spreads so fast and why new variants are spreading faster

abstract: The novel coronavirus SARS CoV-2 responsible for the COVID-19 pandemic and SARS CoV-1 responsible for the SARS epidemic of 2002-2003 share an ancestor yet evolved to have much different transmissibility and global impact1. A previously developed thermodynamic model of protein conformations predicted that SARS CoV-2 is very close to a thermodynamic critical point, which makes it highly infectious but also easily displaced by a spike-based vaccine because there is a tradeoff between transmissibility and robustness2. The model identified a small cluster of four key mutations of SARS CoV-2 that promotes much stronger viral attachment and viral spreading. Here we apply the model to two new strains (B.1.1.7 and B.1.351)3 and predict, using no free parameters, how the new mutations can further enhance infectiousness.

category: biophysics


title: Eukaryotic translation initiation factor 2A protects pancreatic beta cells during endoplasmic reticulum stress while rescuing translation inhibition.

abstract: The endoplasmic reticulum (ER) stress-induced unfolded protein response (UPR) helps decide {beta} cell survival in diabetes. The alternative eukaryotic initiation factor 2A (EIF2A) has been proposed to mediate EIF2S1-independent translation during cellular stress and viral infection, but its role in {beta} cells is unknown. EIF2A abundance is high in human and mouse islets relative to other tissues, and both thapsigargin and palmitate significantly increased EIF2A mRNA and EIF2A protein levels in MIN6 cells, mouse islets and human islets. Knockdowns of EIF2A, the related factor EIF2D, or both EIF2A and EIF2D, were not sufficient to cause apoptosis. On the other hand, transient or stable EIF2A over-expression protected MIN6 cells, primary mouse islets, and human islets from ER stress-induced, caspase-3-dependent apoptosis. Mechanistically, EIF2A overexpression decreased ERN1 (also known as IRE1) expression in thapsigargin-treated MIN6 cells or human islets. In vivo, {beta} cell specific EIF2A viral overexpression reduced ER stress, improved insulin secretion, and abrogated hyperglycemia in Ins2Akita/WT mice. EIF2A overexpression significantly increased expression of genes involved in protein translation and reduced expression of pro-apoptotic genes (e.g. ALDH1A3). Remarkably, the decrease in global protein synthesis during UPR was prevented by EIF2A, despite ER stress-induced EIF2S1 phosphorylation. The protective effects of EIF2A were additive to those of ISRIB, a drug that counteracts the effects of EIF2S1 phosphorylation. Cells overexpressing EIF2A showed higher expression of translation factor EIF2B5, which may contribute to the lack of translational inhibition in these cells. We conclude that EIF2A is a novel target for {beta} cell protection and the circumvention of EIF2S1-mediated translational repression.

category: physiology


title: The Serine Protease HtrA2 mediates radiation-induced senescence in cancer cells

abstract:
Radiation therapy can induce cellular senescence in cancer cells leading to short-term tumor growth arrest, yet increased long-term recurrence. To better understand the molecular mechanisms involved, we developed a model of radiation-induced senescence in cultured cancer cells, which exhibited a typical senescent phenotype, including upregulation of p53 and its target p21, followed by sustained reduction in cellular proliferation, changes in cell size and cytoskeleton organization, and senescence-associated beta-galactosidase activity. A functional siRNA screen using a cell death-related library identified the mitochondrial Ser protease HtrA2 as necessary for senescence development. Mass spectrometry-based proteomic profiling of the senescent cells indicated downregulation of proteins involved in cell cycle progression and DNA repair, and upregulation of proteins associated with malignancy, while irradiation with HtrA2 inhibition upregulated cell proliferation components. In search of direct HtrA2 substrates following radiation, we determined that HtrA2 cleaves the intermediate filament protein vimentin, affecting its cytoplasmic organization. Ectopic expression of active cytosolic HtrA2 resulted in similar changes to vimentin filament assembly. Thus HtrA2, contributes to several hallmarks of senescence and is involved in the cytoskeletal reorganization that accompanies radiation-induced senescence.

SummaryHere the authors identify the Ser protease HtrA2 as a novel mediator of radiation-induced senescence, necessary for sustained proliferation arrest and reorganization of the vimentin filament network.

category: cancer biology


title: Targeting Diacylglycerol Lipase to Reduce AlcoholConsumption

abstract: Alcohol use disorder (AUD) is associated with substantial morbidity, mortality, and societal cost, and pharmacological treatment options for AUD are limited. The endogenous cannabinoid (eCB) signaling system is critically involved in reward processing and alcohol intake is positively correlated with release of the eCB ligand 2-Arachidonoylglycerol (2-AG) within reward neurocircuitry. Here we show that genetic and pharmacological inhibition of diacylglycerol lipase (DAGL), the rate limiting enzyme in the synthesis of 2-AG, reduces alcohol consumption in a variety of preclinical models ranging from a voluntary free-access model to aversion resistant-drinking, and dependence-like drinking induced via chronic intermittent ethanol vapor exposure in mice. DAGL inhibition also prevented ethanol-induced suppression of GABAergic transmission onto midbrain dopamine neurons, providing mechanistic insight into how DAGL inhibition could affect alcohol reward. Lastly, DAGL inhibition during either chronic alcohol consumption or protracted withdrawal was devoid of anxiogenic and depressive-like behavioral effects. These data suggest reducing 2-AG signaling via inhibition of DAGL could represent a novel approach to reduce alcohol consumption across the spectrum of AUD severity.

category: neuroscience


title: The Need for Sustainable Leadership in Academia - a German Case Study

abstract: Academic leaders are selected based on their publication record, citation index and acquisition of third party funding. However, heading a successful research team, also requires leadership skills. Despite the clear need, leadership development has been systematically neglected in the present academic system. At the same time, growing evidence suggests that leadership styles of academic supervisors can dramatically affect the mental health of academic employees as well as drive highly skilled researchers out of academia. Here, we assessed the current state of academic leadership in the German academic system by surveying 368 participants currently employed in academia in Germany. We report that 64% of current academic leaders did not feel prepared for their current position while 86% of participants expressed their interest in leadership development programs offered by their research institutions. Our results highlight the demand for leadership development programs in German academic institutions to ensure a more efficient academic system.

category: scientific communication and education


title: Dynamic predictive coding across the left fronto-temporal language hierarchy: Evidence from MEG, EEG and fMRI

abstract: We used MEG and ERPs to track the time-course and localization of evoked activity produced by expected, unexpected plausible and implausible words during incremental language comprehension. We interpret the full pattern of evoked responses and functional connectivity within a hierarchical predictive coding framework in which evoked activity reflects residual information (error) that is passed up and down the cortical hierarchy. Between 300-500ms, the three conditions produced progressively larger responses within left temporal cortex (lexico-semantic prediction error), while implausible inputs produced a selectively enhanced response within left inferior frontal cortex (event-level prediction error). Between 600-1000ms, unexpected plausible words activated left inferior frontal and middle temporal cortices (feedback activity that induced top-down shifts of event-level and lexico-semantic representations), while highly implausible inputs activated posterior fusiform and medial temporal cortices, supporting orthographic reanalysis and new learning. Therefore, predictive coding may provide a unifying theory that links language comprehension to other domains of cognition.


title: Predictive coding across the left fronto-temporal hierarchy during language comprehension


title: Stink bug Agonoscelis spp. (Heteroptera: Pentatomidae)-An Emerging Threat for Seed Production in Alfalfa Crop (Medicago sativa L.) and Their Successful Management

abstract: Forages are vital constituent for sustainable agriculture because they provide feed for animals that ultimately converted into human food. Alfalfa is one of the most important forages that has highest feeding value for livestock, and seed production of alfalfa seriously affected by several factors, but seed yield loss due to stink bug attack is more as compared to other factors. Studies were conducted to control stink bug by different insecticides at Fodder Research Institute, Sargodha, Punjab-Pakistan during 2016-17. The efficacy of ten insecticides viz., acephate, dimethoate, malathion, chlorpyriphos, bifenthrin, lambdacyhalothrin, deltamethrin, acetamiprid, imidacloprid and carbosulfan were tested against stink bug, Agonoscelis spp. (Heteroptera Pentatomidae) on alfalfa, Medicago sativa L. variety SGD-2002. The mortality of stink bug was recorded one, three, five, seven, ten and fifteen days after spray. Similarly the population of pollinators was recorded before and one, three and five days after spray. From this study it had been observed that acetamiprid (81.14 %) and acephate (80.65%) had the higest mortality of stink bug and proved to be the most effective insecticides against stink bug. By spray of insecticides the population of pollinators declined one day after spray, but it had been rehabilitated three days after spray. By chemical treatment against stink bug, seed yield increased from 28.05 Kg/acre (during last four year without chemical control of stink bug) to 116 Kg/acre in 2016-17 (with chemical control). From this study it has been concluded that chemicals can be used in integrated management program of alfalfa seed production.

category: zoology


title: Midbrain dopamine neurons provide teaching signals for goal-directed navigation

abstract:
In naturalistic environments, animals navigate in order to harvest rewards. Successful goal-directed navigation requires learning to accurately estimate location and select optimal state-dependent actions. Midbrain dopamine neurons are known to be involved in reward value learning1-13. They have also been linked to reward location learning, as they play causal roles in place preference14,15 and enhance spatial memory16-21. Dopamine neurons are therefore ideally placed to provide teaching signals for goal-directed navigation. To test this, we imaged dopamine neural activity as mice learned to navigate in a closed-loop virtual reality corridor and lick to report the reward location. Across learning, phasic dopamine responses developed to visual cues and trial outcome that resembled reward prediction errors and indicated the animals estimate of the reward location. We also observed the development of pre-reward ramping activity, the slope of which was modulated by both learning stage and task engagement. The slope of the dopamine ramp was correlated with the accuracy of licks in the next trial, suggesting that the ramp sculpted accurate location-specific action during navigation. Our results indicate that midbrain dopamine neurons, through both their phasic and ramping activity, provide teaching signals for improving goal-directed navigation.

HighlightsO_LIWe investigated midbrain dopamine activity in mice learning a goal-directed navigation task in virtual reality
C_LIO_LIPhasic dopamine signals reflected prediction errors with respect to subjective estimate of reward location
C_LIO_LIA slow ramp in dopamine activity leading up to reward location developed over learning and was enhanced with task engagement
C_LIO_LIPositive ramp slopes were followed by improved performance on subsequent trials, suggesting a teaching role during goal-directed navigation
C_LI


title: Testing hypotheses of marsupial brain size variation using phylogenetic multiple imputations and a Bayesian comparative framework

abstract: Considerable controversy exists about which hypotheses and variables best explain mammalian brain size variation. We use a new, high-coverage dataset of marsupial brain and body sizes, and the first phylogenetically imputed full datasets of 16 predictor variables, to model the prevalent hypotheses explaining brain size evolution using phylogenetically corrected Bayesian generalised linear mixed-effects modelling. Despite this comprehensive analysis, litter size emerges as the only significant predictor. Marsupials differ from the more frequently studied placentals in displaying much lower diversity of reproductive traits, which are known to interact extensively with many behavioural and ecological predictors of brain size. Our results therefore suggest that studies of relative brain size evolution in placental mammals may require targeted co-analysis or adjustment of reproductive parameters like litter size, weaning age, or gestation length. This supports suggestions that significant associations between behavioural or ecological variables with relative brain size may be due to a confounding influence of the extensive reproductive diversity of placental mammals.


title: Genomic architecture of 5S rDNA cluster and its variations within and between species

abstract: Ribosomal genes (rDNAs) are arranged in purely tandem repeats, preventing them from being reliably assembled onto chromosome. The uncertainty of rDNA genomic structure presents a significant barrier for studying their function and evolution. Here, we generate ultra-long Nanopore and short NGS reads to delineate the architecture and variation of the 5S rDNA cluster in the different strains of C. elegans and C. briggsae. We classify the individual rDNA units into 25 types based on the unique sequence variations in each unit of C. elegans (N2). We next perform manual assembly of the cluster using the long reads that carry these units, which led to an assembly of rDNA cluster consisting of up to 167 5S rDNA units. The ordering and copy number of various rDNA units are indicative of separation time between strains. Surprisingly, we observed a drastically lower level of variation in the 5S rDNA cluster in the C. elegans CB4856 and C. briggsae AF16 strains than C. elegans N2 strain, suggesting a unique mechanism in maintaining the rDNA cluster stability in the N2. Single-copy transgenes landed into the rDNA cluster shows the expected expression in the soma, supporting that rDNA genomic environment is transcriptionally compatible with RNA polymerase II. Delineating the structure and variation of rDNA cluster paves the way for its functional and evolutionary studies.

category: genomics


title: Surface Densities Prewet a Near-Critical Membrane

abstract: Recent work has highlighted roles for thermodynamic phase behavior in diverse cellular processes. Proteins and nucleic acids can phase separate into three-dimensional liquid droplets in the cytoplasm and nucleus and the plasma membrane of animal cells appears tuned close to a two-dimensional liquid-liquid critical point. In some examples, cytoplasmic proteins aggregate at plasma membrane domains, forming structures such as the post-synaptic density and diverse signaling clusters. Here we examine the physics of these surface densities, employing minimal simulations of co-acervating polymers coupled to an Ising membrane surface in conjunction with a complementary Landau theory. We argue that these surface densities are a novel phase reminiscent of pre-wetting, in which a molecularly thin three-dimensional liquid forms on a usually solid surface. However, in surface densities the solid surface is replaced by a membrane with an independent propensity to phase separate. We show that proximity to criticality in the membrane dramatically increases the parameter regime in which a pre-wetting-like transition occurs, leading to a broad region where coexisting surface phases can form even when a bulk phase is unstable. Our simulations naturally exhibit three surface phase coexistence even though both the membrane and the polymer bulk can only display two phase coexistence on their own. We argue that the physics of these surface densities enables diverse functions seen in Eukaryotic cells.


title: Impaired pattern separation in Tg2576 mice is associated with hyperexcitable dentate gyrus caused by Kv4.1 downregulation

abstract: Alzheimers disease (AD) is a progressive neurodegenerative disorder that causes memory loss. Most AD researches have focused on neurodegeneration mechanisms. Considering that neurodegenerative changes are not reversible, understanding early functional changes before neurodegeneration is critical to develop new strategies for early detection and treatment of AD. We found that Tg2576 mice exhibited impaired pattern separation at the early preclinical stage. Based on previous studies suggesting a critical role of dentate gyrus (DG) in pattern separation, we investigated functional changes in DG of Tg2576 mice. We found that granule cells in DG (DG-GCs) in Tg2576 mice showed increased action potential firing in response to long depolarizations and reduced 4-AP sensitive K+-currents compared to DG-GCs in wild-type (WT) mice. Among Kv4 family channels, Kv4.1 mRNA expression in DG was significantly lower in Tg2576 mice. We confirmed that Kv4.1 protein expression was reduced in Tg2576, and this reduction was restored by antioxidant treatment. Hyperexcitable DG and impaired pattern separation in Tg2576 mice were also recovered by antioxidant treatment. These results highlight the hyperexcitability of DG-GCs as a pathophysiologic mechanism underlying early cognitive deficits in AD and Kv4.1 as a new target for AD pathogenesis in relation to increased oxidative stress.


title: Relict groups of spiny frogs indicate Late Paleogene-Early Neogene trans-Tibet dispersal of thermophile faunal elements

abstract:
BackgroundThe Himalaya-Tibet orogen (HTO) presents an outstanding geologically active formation that contributed to, and fostered, modern Asian biodiversity. However, our concepts of the historical biogeography of its biota are far from conclusive, as are uplift scenarios for the different parts of the HTO. Here, we revisited our previously published data set of the tribe Paini extending it with sequence data from the most western Himalayan spiny frogs Allopaa and Chrysopaa and using them as an indirect indicator for the paleoecological development of Tibet.

MethodsWe obtained sequence data of two mitochondrial loci (16S rRNA, COI) and one nuclear marker (Rag1) from Allopaa samples from Kashmir Himalaya as well as Chrysopaa sequence data from the Hindu Kush available from GenBank to complement our previous data set. A Maximum likelihood and dated Bayesian gene tree were generated based on the concatenated data set. To resolve the inconsistent placement of Allopaa, we performed different topology tests.

ResultsConsistent with previous results, the Southeast Asian genus Quasipaa is sister to all other spiny frogs. The results further reveal a basal placement of Chrysopaa relative to Allopaa and Nanorana with an estimated age of ca. 26 Mya. Based on the topology tests, the phylogenetic position of Allopaa as a sister clade to Chaparana seems to be most likely, resulting in a paraphyletic genus Nanorana and a separation from the latter clade around 20 Mya. Both, the placements of Chrysopaa and Allopaa support the presence of basal Paini lineages in the far north western part of the HTO, which is diametrically opposite end of the HTO with respect to the ancestral area of spiny frogs in Southeast Asia. These striking distributional patterns can be most parsimoniously explained by trans-Tibet dispersal during the late Oligocene (subtropical Chrysopaa) respectively early Miocene (warm temperate Allopaa). Within spiny frogs, only members of the monophyletic Nanorana+Paa clade are adapted to the colder temperate climates, indicating that high-altitude environments did not dominate in the HTO before ca. 15 Mya. Our results are consistent with fossil records suggesting that large parts of Tibet were characterized by subtropical to warm temperate climates at least until the early Miocene.


title: Talaromyces marneffei promotes M2 polarization of human macrophages by downregulating SOCS3 expression and activating TLR9 pathway

abstract: Little is known about how Talaromyces marneffei, a thermally dimorphic fungus that causes substantial morbidity and mortality in Southeast Asia, evades the human immune system. Polarization of macrophages into fungal-inhibiting M1 and fungal-promoting M2 types has been shown to play an important role in the innate immune response against fungal pathogens. This mechanism has not been defined for T. marneffei. Here, we demonstrated that T. marneffei promotes its survival in human macrophages by inducing them towards M2 polarization. Our investigations of the mechanism revealed that T. marneffei infection led to SOCS3 protein degradation by inducing tyrosine phosphorylation, thereby relieving the inhibitory effect of SOCS3 on p-STAT6, a key factor for M2 polarization. Our SOCS3-overexpression experiments showed that SOCS3 is a positive regulator of M1 polarization and plays an important role in limiting M2 polarization. Furthermore, we found that inhibition of TLR9 pathway partially blocked T. marneffei-induced M2 polarization and significantly enhanced the killing activity of macrophages against T. marneffei. Collectively, these results reveal a novel mechanism by which T. marneffei evades the immune response of human macrophages.

category: microbiology


title: Glycine receptor α3K governs mobility and conductance of L/K splice variant heteropentamers

abstract:
Glycine receptors (GlyRs) are ligand-gated pentameric chloride channels in the central nervous system. GlyR-3 is a possible target for chronic pain treatment and temporal lobe epilepsy. Alternative splicing into K or L variants determines the subcellular fate and function of GlyR-3, yet it remains to be shown whether its different splice variants can functionally co-assemble, and what the properties of such heteropentamers would be. Here, we subjected GlyR-3 to a combined fluorescence microscopy and electrophysiology analysis. We employ masked Pearsons and dual-color spatiotemporal correlation analysis to prove that GlyR-3 splice variants heteropentamerize, adopting the mobility of the K variant. Fluorescence-based single-subunit counting experiments revealed a variable and concentration ratio dependent hetero-stoichiometry. Via single-channel on-cell patch clamp we show heteropentameric conductances resemble those of the 3K splice variant. Our data are compatible with a model where 3 heteropentamerization fine-tunes mobility and activity of GlyR 3 channels, which is important to understand and tackle 3 related diseases.

SummaryThe glycine receptor 3 is key to the central nervous systems physiology and involved in chronic pain and epilepsy. In this paper, Lemmens et al. reveal and functionally characterize 3 splice variant heteropentamerization via advanced single-molecule fluorescence image analysis.

DeclarationsO_ST_ABSFundingC_ST_ABSWe acknowledge the UHasselt Advanced Optical Microscopy Centre (AOMC). Prof. Em. Marcel Ameloot, the Research Foundation Flanders (FWO, project G0H3716N) and the province of Limburg (Belgium) (tUL Impuls II) are acknowledged for funding the microscopy hardware. V. Lemmens is grateful for a doctoral scholarship from the UHasselt (17DOC11BOF) and KU Leuven (C14/16/053) Special Research Funds (BOF).

Conflicts of interest / competing interestsNo conflicts of interest apply.

Ethics approvalNot applicable

Availability of data and materialAll data and material are available upon request.

Code availabilityFluctuation imaging and co-localization analyses were performed in the software package PAM [71]. The software is available as source code, requiring MATLAB to run, or as pre-compiled standalone distributions for Windows or MacOS at http://www.cup.uni-muenchen.de/pc/lamb/software/pam.html or hosted in Git repositories under http://www.gitlab.com/PAM-PIE/PAM and http://www.gitlab.com/PAM-PIE/PAMcompiled. A detailed user manual is available at http://pam.readthedocs.io.

Author contributionsConceptualization Meier J.C., Brone B. and Hendrix J.; Investigation and formal analysis Lemmens V. and Thevelein B.; Software development Hendrix J.; Writing the original draft Lemmens V., Thevelein B and Hendrix, J.; Review and editing by all authors.


title: Daily caffeinated soda intake is associated with impaired working memory and higher impulsivity in children

abstract:
Caffeinated soda contains two addictive substances, sugar and caffeine, and is the most preferred means of caffeine consumption among children. However, it remains unexplored if habitual caffeinated soda intake in childhood is associated with a higher risk of alcohol misuse in the future. Here, we investigated the neurocognitive correlates of daily consumption of caffeinated soda and examined whether daily caffeinated soda intake is associated with a higher risk of alcohol initiation in children. By applying machine learning and hierarchical regression to a large dataset from the Adolescent Brain Cognitive Development (ABCD) Study, we investigated the relationship between daily caffeinated soda intake and the well-known risk factors of substance misuse, including impaired working memory, high impulsivity, and altered reward processing. We next examined whether daily caffeinated soda intake could predict greater alcohol sipping in 12 months. The results showed that daily consumption of caffeinated soda was associated with impaired working memory, higher impulsivity, and altered reward processing, with increased amygdala activation during the emotional working memory task and hypoactivity in the nucleus accumbens and posterior cingulate cortex during reward processing tasks. Daily drinkers reported higher rates of alcohol sipping compared to non-drinkers after 12 months, and daily intake of caffeinated soda mediated the relationship between the neurocognitive risk factors and future alcohol sipping. These results suggest that daily caffeinated soda intake in childhood (1) is associated with impaired neurocognitive functioning, and (2) may work as a predictor for alcohol initiation, which has significant implications for public health recommendations.

Significance StatementIs daily caffeinated soda intake associated with higher risks of substance misuse? If so, which specific risk factors of substance misuse are involved? Moreover, can daily caffeinated soda intake predict future alcohol initiation? These are the questions that many parents are asking but surprisingly there is no clear guideline, and its potential vulnerabilities of substance misuse remain largely unknown. Here, by applying machine learning and hierarchical regression approaches, we demonstrate that daily caffeinated soda intake is associated with neurocognitive risk factors of substance misuse including impaired working memory, higher impulsivity, and altered reward processing. We also present that daily caffeinated soda consumption could predict higher alcohol sipping rates in 12 months. These results have significant implications for public health recommendations.


title: Daily caffeinated soda intake in children is associated with neurocognitive vulnerabilities of substance misuse and predicts higher risk of alcohol sipping in 12 months


title: A covalently crosslinked bioink for multi-materials drop-on-demand 3D bioprinting of three-dimensional cell cultures

abstract: In vitro three-dimensional (3D) cell models have been accepted to better recapitulate aspects of in vivo organ environment than 2D cell culture. Currently, the production of these complex in vitro 3D cell models with multiple cell types and microenvironments remains challenging and prone to human error. Here we report a versatile bioink comprised of a 4-arm PEG based polymer with distal maleimide derivatives as the main ink component and a bis-thiol species as the activator that crosslinks the polymer to form the hydrogel in less than a second. The rapid gelation makes the polymer system compatible with 3D bioprinting. The ink is combined with a drop-on-demand 3D bioprinting platform consisting of eight independently addressable nozzles and high-throughput printing logic for creating complex 3D cell culture models. The combination of multiple nozzles and fast printing logic enables the rapid preparation of many complex 3D structures comprising multiple hydrogel environments in the one structure in a standard 96-well plate format. The platform compatibility for biological applications was validated using pancreatic ductal adenocarcinoma cancer (PDAC) cells with their phenotypic responses controlled by tuning the hydrogel microenvironment.

category: bioengineering


title: Contribution of interneuron subtype-specific GABAergic signalling to emergent sensory processing in somatosensory whisker barrel cortex in mouse.

abstract: Mammalian neocortex is important for conscious processing of sensory information. Fundamental to this function is balanced glutamatergic and GABAergic signalling. Yet little is known about how this interaction arises in the developing forebrain despite increasing insight into early GABAergic interneuron (IN) circuits. To further study this, we assessed the contribution of specific INs to the development of sensory processing in the mouse whisker barrel cortex. Specifically we explored the role of INs in speed coding and sensory adaptation. In wild-type animals, both speed processing and adaptation were present as early as the layer 4 critical period of plasticity, and showed refinement over the period leading to active whisking onset. We then conditionally silenced action-potential-dependent GABA release in either somatostatin (SST) or vasoactive intestinal peptide (VIP) INs. These genetic manipulations influenced both spontaneous and sensory-evoked activity in an age and layer-dependent manner. Silencing SST+ INs reduced early spontaneous activity and abolished facilitation in sensory adaptation observed in control pups. In contrast, VIP+ IN silencing had an effect towards the onset of active whisking. Silencing either IN subtype had no effect on speed coding. Our results reveal how these IN subtypes differentially contribute to early sensory processing over the first few postnatal weeks.


title: Natural variation in fertility is correlated with species-wide levels of divergence in Caenorhabditis elegans

abstract: Life history traits underlie the fitness of organisms and are under strong natural selection in the face of environmental challenges. A new mutation that positively impacts a life history trait will likely increase in frequency and become fixed in a population (e.g. selective sweep). The identification of the beneficial alleles that underlie selective sweeps provides insights into the mechanisms that occurred during the evolution of species. In the global population of Caenorhabditis elegans, we previously identified selective sweeps that have drastically reduced chromosomal-scale genetic diversity in the species. Here, we measured the fertility (viable offspring) of a collection of wild C. elegans strains, including many recently isolated divergent strains from the Hawaiian islands and found that strains with larger swept genomic regions on multiple chromosomes have significantly higher fertility than strains that do not have evidence of the recent selective sweeps. We used genome-wide association (GWA) mapping to identify three quantitative trait loci (QTL) underlying the fertility variation. Additionally, we mapped previous fertility data of wild C. elegans strains and C. elegans recombinant inbred advanced intercross lines (RIAILs) that were grown in various conditions and detected eight QTL across the genome using GWA and linkage mappings. These QTL show the genetic complexity of life history traits such as fertility across this species. Moreover, the haplotype structure in each GWA QTL region revealed correlations with recent selective sweeps in the C. elegans population. North American and European strains had significantly higher fertility than most strains from Hawaii, a hypothesized origin of the C. elegans species, suggesting that beneficial alleles that cause increased fertility could underlie the selective sweeps during the worldwide expansion of C. elegans.

category: genetics


title: Retrograde Control of Cytosolic Translation Targets Synthesis of Plastid Localized Proteins and Nuclear Responses for Efficient Light Acclimation.

abstract: Canonical retrograde signaling is the transmission of information from organelles to the nucleus. Discrepancies between protein accumulation and transcript abundance in response to oxidative stress were suggestive of protein translation responding to retrograde signaling. Here we uncover multiple components of a translation-dependent retrograde signaling pathway that impact translation efficiency and gene expression, including the kinases, MPK6 and the SnRK1 subunit, AKIN10. Global ribosome foot-printing demonstrated rapid differential loading of 939 of transcripts from polyribosomes within 10 min after transfer from Low to High-light. Translationally regulated transcripts shared motifs in their 5-UTR that act as binding sites for RBPs such as GAPC. The Stress Associated Proteins 2 and 3 carry such motifs in their UTRs and interact with the calcium sensor Calmodulin-like 49, relocating to the nucleus to co-regulate a translation-dependent transcriptional response. Translation dependent retrograde signaling bifurcates into a direct translational circuit and a translation-reliant nuclear circuit synchronizing translation, nuclear and anterograde response pathways, which may serve as a just in time-provision of needed proteins to the plastids.

category: plant biology


title: Stress-induced glucocorticoids alter the rhythmic activity of Leydig cells depending on the circadian time

abstract:
BackgroundThe increasing amount of data points to the circadian timing system as an essential part of processes regulating androgen homeostasis. However, the relationship between stress response, timekeeping-, and steroidogenesis-related systems is unexplored.

ObjectiveThe purpose of the study was to analyze the stress-response of the testosterone-producing Leydig cells depending on the stressful events time.

Materials and methodsThe study was designed to follow the effects of 3-hour immobilization (IMO) applied at different periods during the day. The IMO performed once (1xIMO) or repeated in 10 consecutive days (10xIMO). Principal-component-analysis (PCA) followed the expression study of the clock and steroidogenic-related genes in Leydig cells.

ResultsBoth types of IMO in all investigated periods increased corticosterone and decreased testosterone blood level. Transcriptional analysis revealed different sensitivity to IMO events depending on the circadian time. The majority of steroidogenesis-related genes (Lhcgr, Cyp11a1, Cyp17a1, Hsd3b1/2) were down-regulated in the inactive but unchanged or even up-regulated in the active phase of the day. Both types of IMO potentiated the expression of clock elements Bmal1/BMAL1, Per1/PER1 regardless of the days stage and reduced Reverba in the inactive phase. The PCA confirmed a major shift, for both IMO-types, in the transcription of steroidogenesis and clock genes across the inactive/active phase. Further, the diurnal pattern of the glucocorticoid receptor (Nr3c1/GR) expression in Leydig cells was increased in the inactive phase due to 10xIMO. The observed time-dependent IMO-response of the Leydig cells correlated with different GR engagements.

DiscussionStress- and the circadian-system coordinatively shape Leydig cells physiology, assuming diverse GR engagement as a possible factor in mediating the diurnal effect of stress.

ConclusionThe Leydig cells stress-response depends on the time of the stressful situation, emphasizing the importance of circadian activity in supporting androgen homeostasis and male fertility.


title: Stress alters the activity of Leydig cells dependently on the diurnal time


title: Meta-analysis of orthogonal OMICs data from COVID-19 patients unveils prognostic markers and antiviral factors.

abstract: While our battle with the COVID-19 pandemic continues, a multitude of Omics data has been generated from patient samples in various studies, which remains to be translated. We conducted a meta-analysis of published transcriptome and proteome profiles of nasal swab and bronchioalveolar lavage fluid (BALF) samples of COVID-19 patients, to shortlist high confidence upregulated host factors. Subsequently, mRNA overexpression of selected genes was validated in nasal swab/BALF samples from a cohort of COVID-19 positive/negative, symptomatic/asymptomatic individuals. Analysis of these data revealed S100 family genes (S100A6, S100A8, S100A9, and S100P) as prognostic markers of COVID-19 disease. Furthermore, Thioredoxin gene (TXN) was identified as a significant upregulated host factor in our overlap analysis. An FDA-approved drug Auranofin, which inhibits Thioredoxin reduction, was found to mitigate SARS-CoV-2 replication in vitro and in vivo in the hamster challenge model. Overall, this study translates COVID-19 host response Big Data into potential clinical interventions.


title: Identification of COVID-19 prognostic markers and therapeutic targets through meta-analysis and validation of Omics data from nasopharyngeal samples


title: BRAF inhibitors reprogram cancer-associated fibroblasts to drive matrix remodeling and therapeutic escape in melanoma

abstract: The tumor stroma and its cellular components are known to play an important role in tumor resilience towards treatment. Here, we report a novel resistance mechanism in melanoma that is elicited by BRAF inhibitors (BRAFi)-induced non-canonical nuclear {beta}-catenin signaling in cancer-associated fibroblasts (CAFs). Our study reveals that BRAFi leads to an expanded CAF population with increased {beta}-catenin nuclear accumulation and enhances their biological properties. This CAF subpopulation is essential for melanoma cells to grow and resist to BRAFi/MEK inhibitors (MEKi). Mechanistically, BRAFi induces BRAF and CRAF dimerization and subsequently the activation of ERK signaling in CAFs, leading to the inactivation of the {beta}-catenin destruction complex. By RNA-Seq, periostin (POSTN) is identified as a major downstream effector of {beta} catenin in CAFs and can compensate for the loss of {beta}-catenin in CAFs in conferring melanoma cell BRAFi/MEKi resistance. Moreover, POSTN reactivates the ERK pathway, which is inhibited by BRAFi/MEKi, in melanoma cells through PI3K/AKT signaling. Our data underscore the roles of BRAFi-induced CAF reprogramming in matrix remodeling and therapeutic escape of melanoma cells and reveal POSTN as an important matrix target to eliminate BRAFi/MEKi resistance in melanoma.


title: A genetic program boosts mitochondrial function to power macrophage tissue invasion

abstract: Metabolic adaptation to changing demands underlies homeostasis. During inflammation or metastasis, cells leading migration into challenging environments require an energy boost, however what controls this capacity is unknown. We identify a previously unstudied nuclear protein, Atossa, as changing metabolism in Drosophila melanogaster immune cells to promote tissue invasion. Atossas vertebrate orthologs, FAM214A-B, can fully substitute for Atossa, indicating functional conservation from flies to mammals. Atossa increases mRNA levels of Porthos, an unstudied RNA helicase and two metabolic enzymes, LKR/SDH and GR/HPR. Porthos increases translation of a gene subset, including those affecting mitochondrial functions, the electron transport chain, and metabolism. Respiration measurements and metabolomics indicate that Atossa and Porthos powers up mitochondrial oxidative phosphorylation to produce sufficient energy for leading macrophages to forge a path into tissues. As increasing oxidative phosphorylation enables many crucial physiological responses, this unique genetic program may modulate a wide range of cellular behaviors beyond migration.


title: Conservation and convergence of genetic architecture in the adaptive radiation of Anolis lizards

abstract: The G matrix, which quantifies the genetic architecture of traits, is often viewed as an evolutionary constraint. However, G can evolve in response to selection and may also be viewed as a product of adaptive evolution. Convergent evolution of G in similar environments would suggest that G evolves adaptively, but it is difficult to disentangle such effects from phylogeny. Here, we use the adaptive radiation of Anolis lizards to ask whether convergence of G accompanies the repeated evolution of habitat specialists, or ecomorphs, across the Greater Antilles. We measured G in seven species representing three ecomorphs (trunk-crown, trunk- ground, and grass-bush). We found that the overall structure of G does not converge. Instead, the structure of G is well conserved and displays a phylogenetic signal consistent with Brownian motion. However, several elements of G showed signatures of convergence, indicating that some aspects of genetic architecture have been shaped by selection. Most notably, genetic correlations between limb traits and body traits were weaker in long-legged trunk-ground species, suggesting effects of recurrent selection on limb length. Our results demonstrate that common selection pressures may have subtle but consistent effects on the evolution of G, even as its overall structure remains conserved.


title: Membrane-binding mechanism of the EEA1 FYVE domain revealed by multi-scale molecular dynamics simulations

abstract: Early Endosomal Antigen 1 (EEA1) is a key protein in endosomal trafficking and is implicated in both autoimmune and neurological diseases. The C-terminal FYVE domain of EEA1 binds endosomal membranes, which contain phosphatidylinositol-3-phosphate (PI(3)P). Although it is known that FYVE binds PI(3)P specifically, it has not previously been described of how FYVE attaches and binds to endosomal membranes. In this study, we employed both coarse-grained (CG) and atomistic (AT) molecular dynamics (MD) simulations to determine how FYVE binds to PI(3)P-containing membranes. CG-MD showed that the dominant membrane binding mode resembles the crystal structure of EEA1 FYVE domain in complex with inositol-1,3-diphospate (PDB ID 1JOC). FYVE, which is a homodimer, binds the membrane via a hinge mechanism, where the C-terminus of one monomer first attaches to the membrane, followed by the C-terminus of the other monomer. The total binding energy is 70 kJ/mol, of which 50-60 kJ/mol stems from specific PI(3)P-interactions. By AT-MD, we could partition the binding mode into two types: (i) adhesion by electrostatic FYVE-PI(3)P interaction, and (ii) insertion of amphipathic loops. The AT simulations also demonstrated flexibility within the FYVE homodimer between the C-terminal heads and coiled-coil stem, allowing binding via a mechanism resembling that of a suction cup connected to a locally rigid stem via a flexible hinge.


title: Spatially confined sub-tumor microenvironments orchestrate pancreatic cancer pathobiology

abstract:
Pancreatic ductal adenocarcinoma (PDAC) remains resistant to most treatments and demonstrates a complex pathobiology. Here, we deconvolute regional heterogeneity in the human PDAC tumor microenvironment (TME), a long-standing obstacle, to define precise stromal contributions to PDAC progression. Large scale integration of histology-guided multiOMICs with clinical data sets and functional in vitro models uncovers two microenvironmental programs in PDAC that were anchored in fibroblast differentiation states. These sub-tumor microenvironments (subTMEs) co-occurred intratumorally and were spatially confined, producing patient-specific cellular and molecular heterogeneity associated with shortened patient survival. Each subTME was uniquely structured to support discrete aspects of tumor biology: reactive regions rich in activated fibroblast communities were immune-hot and promoted aggressive tumor progression while deserted regions enriched in extracellular matrix supported tumor differentiation yet were markedly chemoprotective. In conclusion, PDAC regional heterogeneity derives from biologically distinct reactive and protective TME elements with a defined, active role in PDAC progression.

Graphical Abstract & Key findings

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=191 SRC="FIGDIR/small/431890v1_ufig1.gif" ALT="Figure 1">
View larger version (73K):
org.highwire.dtl.DTLVardef@17adabdorg.highwire.dtl.DTLVardef@aee778org.highwire.dtl.DTLVardef@deae00org.highwire.dtl.DTLVardef@1718272_HPS_FORMAT_FIGEXP  M_FIG C_FIG O_LIPDAC regional heterogeneity originates in sub-tumor microenvironments (subTMEs)
C_LIO_LISubTMEs exhibit distinct immune phenotypes and CAF differentiation states
C_LIO_LIDifferent subTMEs are either tumor-promoting or chemoprotective
C_LIO_LIIntratumoral subTME co-occurrence links stromal heterogeneity to patient outcome
C_LI


title: Rapid Characterization of AAV gene therapy vectors by Mass Photometry

abstract: Recombinant adeno-associated viruses (rAAV) are extensively used as gene delivery vectors in clinical studies, and several rAAV based treatments have already been approved. Significant progress has been made in rAAV manufacturing, and large-scale vector production and purification methods have been developed. However, a better and more precise capsid characterization techniques are still needed to guarantee the purity and safety of the rAAV preparations. A recently developed single-molecule technique, mass photometry (MP), measures mass distributions of biomolecules with high resolution and sensitivity. Here we explore applications of MP for the characterization of capsid fractions. We demonstrate that MP is able to resolve and quantify not only empty and full-genome containing capsid populations, but also identify the partially packaged capsid impurities. MP data accurately measures full and empty capsid ratios, and can be used to estimate the size of the encapsidated genome. MP distributions provide information on sample heterogeneity and on the presence of aggregates. Current analytical techniques used to characterize rAAV preparations are susceptible to background signals, have limited accuracy, or are time-consuming and require a large amount of material. MP can analyze sub-picomole quantities of sample, and data can be obtained and analyzed within minutes. This method provides a simple, robust, and effective tool to monitor physical attributes of rAAV vectors.


title: Cultivation-based quantification and identification of bacteria at two hygienic key sides of domestic washing machines

abstract:
Detergent drawer and door seal represent important sites for microbial life in domestic washing machines. Interestingly, quantitative data on microbial contamination of these sites is scarce. Here, 10 domestic washing machines were swap-sampled for subsequent bacterial cultivation at four different sampling sites, each: detergent drawer, detergent drawer chamber as well as top and bottom part of the rubber door seal. The average bacterial load over all washing machines and sites was 2.1 {+/-} 1.0 x 104 CFU cm-2 (average {+/-} standard error of the mean (SEM)). The top part of the door seal showed the lowest contamination (11.1 {+/-} 9.2 x 101 CFU cm-2), probably due to less humidity.

Out of 212 isolates, 178 (84%) were identified on genus level and 118 (56%) on species level using MALDI biotyping, resulting in 29 genera and 40 identified species across all machines. The predominant bacterial genera were Staphylococcus and Micrococcus, which were found at all sites. 21 out of 40 species were classified as opportunistic pathogens, emphasizing the need for regular cleaning of the investigated sites.


title: Comparing optimized exoskeleton assistance of the hip, knee, and ankle in single and multi-joint configurations

abstract: Exoskeletons that assist the hip, knee, and ankle joints have begun to improve human mobility, particularly by reducing the metabolic cost of walking. However, direct comparisons of optimal assistance of these joints, or their combinations, have not yet been possible. Assisting multiple joints may be more beneficial than the sum of individual effects, because muscles often span multiple joints, or less effective, because single-joint assistance can indirectly aid other joints. In this study, we used a hip-knee-ankle exoskeleton emulator paired with human-in-the-loop optimization to find single-joint, two-joint, and whole-leg assistance that maximally reduced the metabolic cost of walking for three participants. Hip-only and ankle-only assistance reduced the metabolic cost of walking by 26% and 30% relative to walking in the device unassisted, confirming that both joints are good targets for assistance. Knee-only assistance reduced the metabolic cost of walking by 13%, demonstrating that effective knee assistance is possible. Two-joint assistance reduced the metabolic cost of walking by between 34% and 42%, with the largest improvements coming from hip-ankle assistance. Assisting all three joints reduced the metabolic cost of walking by 50%, showing that at least half of the metabolic energy expended during walking can be saved through exoskeleton assistance. Changes in kinematics and muscle activity indicate that single-joint assistance indirectly assisted muscles at other joints, such that the improvement from whole-leg assistance was smaller than the sum of its single-joint parts. Exoskeletons can assist the entire limb for maximum effect, but a single well-chosen joint can be more efficient when considering additional factors such as weight and cost.


title: Extracellular riboflavin induces anaerobic biofilm formation in Shewanella oneidensis

abstract:
Some microorganisms can respire with extracellular electron acceptors using an extended electron transport chain to the cell surface. These organisms apply flavin molecules as cofactors to facilitate one-electron transfer catalysed by the terminal reductases and as endogenous electron shuttles. In the model organism Shewanella oneidensis, riboflavin production and excretion triggers a specific biofilm formation response that is initiated at a specific threshold concentration, similar to canonical quorum sensing molecules. Riboflavin-mediated messaging is based on the overexpression of the gene encoding the putrescin decarboxylase speC which leads to posttranscriptional overproduction of proteins involved in biofilm formation. Using a model of growth-dependent riboflavin production under batch and biofilm growth conditions, the number of cells necessary to produce the threshold concentration per time was deduced. Furthermore, our results indicate that specific retention of riboflavin in the biofilm matrix leads to localized concentrations which by far exceed the necessary threshold value.

ImportanceFerric iron is the fourth most abundant element of the earth crust. It occurs at neutral pH in the form insoluble iron minerals. The dissimilatory reduction of these minerals is an import part of global geological cycles and is catalyzed by microorganisms through extended respiratory chains to the cell surface. Shewanella oneidensis is one of the best understood model organisms for this kind of extracellular respiration. Flavins are important for the reduction of extracellular electron acceptors by S. oneidensis. since they have a function as (I) cofactors of the terminal reductases and (II) electron shuttles. In this study we reveal that flavin molecules are further employed as quorum sensing molecules. They are excreted by the organisms in a growth dependent manner and lead to anaerobic biofilm formation as a specific response at a certain threshold concentration. Although we know multiple examples of quorum sensing mechanisms, the use of riboflavin was so far not described and at least in S. oneidensis proceeds via a new regulatory routine that proceeds on the trancriptomic and posttranscriptomic level.


title: Drivers of inter-population variation in the gut microbiomes of sister species of Phanaeus dung beetles

abstract: The microbiome plays key roles in host physiology and ecology, but how the microbiome varies among populations of hosts is not well understood. However, different abiotic and biotic selection pressures across a species range likely lead to variation in the microbiome. In addition, symbiotic microbiota may differ more between closely-related species in sympatry than in allopatry if selection favors the reduction of interspecific competition. We investigated variation in the maternally-transmitted, beneficial gut microbiomes of Phanaeus vindex and P. difformis, sister species of dung beetle that compete for the same resources in sympatry and occur across diverse climatic conditions. We sampled and sequenced bacterial/archaeal 16S rDNA from guts of P. difformis and P. vindex collected across 17 sympatric and allopatric sites. Gut microbial communities were best predicted by temperature and precipitation, cattle present at sites, and spatial relationships among sites. Contrary to our hypotheses, we did not find that the gut microbial communities of P. vindex and P. difformis differed more in sympatry than in allopatry, nor that P. vindex, the more broadly distributed of the two species, exhibited greater microbiome turnover among populations. However, the gut microbiome of P. vindex shifted more between sympatric and allopatric populations than did that of P. difformis, suggesting character displacement. While more research is needed, it is possible that differences in gut microbial communities allow P. vindex and P. difformis to partition their niches in sympatry. Our work argues for further exploration of the gut microbiomes potential role in niche partitioning and local adaptation.

category: ecology


title: Promoting axon regeneration by enhancing the non-coding function of the injury-responsive coding gene Gpr151

abstract: Gene expression profiling in response to nerve injury has been mainly focused on protein functions of coding genes to understand mechanisms of axon regeneration and to identify targets of potential therapeutics for nerve repair. However, the protein functions of several highly injury-induced genes including Gpr151 for regulating the regenerative ability remain unclear. Here we present an alternative approach focused on non-coding functions of the coding genes, which led to the identification of the non-coding function of Gpr151 RNA interacting with RNA- binding proteins such as CSDE1. Gpr151 promotes axon regeneration by the function of its 5- untranslated region (5UTR) and expression of an engineered form of the 5UTR improves regenerative capacity in vitro and in vivo in both sciatic nerve and optic nerve injury models. Our data suggest that searching injury-induced coding genes potentially functioning by their non- coding regions is required for the RNA-based gene therapy for improving axon regeneration.


title: Efficiency Improvements and Discovery of New Substrates for a SARS-CoV-2 Main Protease FRET Assay

abstract: The COVID-19 pandemic, caused by the SARS-CoV-2 virus, has a huge impact on the world. Although several vaccines have recently reached the market, the development of specific antiviral drugs against SARS-CoV-2 is an important additional strategy in fighting the pandemic. One of the most promising pharmacological targets is the viral main protease (Mpro). Here, we present an optimized biochemical assay procedure for SARS-CoV-2 Mpro. We have comprehensively investigated the influence of different buffer components and conditions on the assay performance, and characterized six FRET substrates with a 2-Abz/Tyr(3-NO2) FRET pair. The substrates 2-AbzSAVLQSGTyr(3-NO2)R-OH, a truncated version of the established DABCYL/EDANS FRET substrate, and a new substrate 2-AbzVVTLQSGTyr(3-NO2)R-OH are promising candidates for screening and inhibitor characterization. In the latter substrate, the incorporation of Val at the position P5 improved the catalytic efficacy. Based on the obtained results, we present here a reproducible, reliable assay protocol using highly affordable buffer components.

category: biochemistry


title: Notch-Jagged signaling complex defined by an interaction mosaic

abstract: The Notch signaling system links cellular fate to that of its neighbors, driving proliferation, apoptosis, and cell differentiation in metazoans, whereas dysfunction leads to debilitating developmental disorders and cancers. Other than a five-by-five domain complex, it is unclear how the 40 extracellular domains of the Notch1 receptor collectively engage the 19 domains of its canonical ligand Jagged1 to activate Notch1 signaling. Here, using cross-linking mass spectrometry (XL-MS), biophysical and structural techniques on the full extracellular complex and targeted sites,we identify five distinct regions, two on Notch1 and three on Jagged1, that form an interaction network.The Notch1 membrane-proximal regulatory region individually binds to the established Notch1 epidermal growth factor (EGF) 8-13 and Jagged1 C2-EGF3 activation sites, as well as to two additional Jagged1 regions, EGF 8-11 and cysteine-rich domain (CRD). XL-MS and quantitative interaction experiments show that the three Notch1 binding sites on Jagged1 also engage intramolecularly.These interactions, together with Notch1 and Jagged1 ectodomain dimensions and flexibility determined by small-angle X-ray scattering (SAXS), support the formation of backfolded architectures. Combined, the data suggest that critical Notch1 and Jagged1 regions are not distal, but engage directly to control Notch1 signaling, thereby redefining the Notch1-Jagged1 activation mechanism and indicating new routes for therapeutic applications.

category: molecular biology


title: Delivery of recombinant SARS-CoV-2 envelope protein into human cells

abstract: SARS-CoV-2 envelope protein (S2-E) is a conserved membrane protein that is essential to coronavirus assembly and budding. Here, we describe the recombinant expression and purification of S2-E into amphipol-class amphipathic polymer solutions. The physical properties of amphipols underpin their ability to solubilize and stabilize membrane proteins without disrupting membranes. Amphipol delivery of S2-E to pre-formed planar bilayers results in spontaneous membrane integration and formation of viroporin ion channels. Amphipol delivery of the S2-E protein to human cells results in membrane integration followed by retrograde trafficking to a location adjacent to the endoplasmic reticulum-to-Golgi intermediate compartment (ERGIC) and the Golgi, which are the sites of coronavirus replication. Delivery of S2-E to cells enables both chemical biological approaches for future studies of SARS-CoV-2 pathogenesis and development of "Trojan Horse" anti-viral therapies. This work also establishes a paradigm for amphipol-mediated delivery of membrane proteins to cells.


title: mRNA codon optimization on quantum computers

abstract: Reverse translation of polypeptide sequences to expressible mRNA constructs is a NP-hard combinatorial optimization problem. Each amino acid in the protein sequence can be represented by as many as six codons, and the process of selecting the combination that maximizes probability of expression is termed codon optimization. This work investigates the potential impact of leveraging quantum computing technology for codon optimization. An adiabatic quantum computer (AQC) is compared to a standard genetic algorithm (GA) programmed with the same objective function. The AQC is found to be competitive in identifying optimal solutions and future generations of AQCs may be able to outperform classical GAs. The utility of gate-based systems is also evaluated using a simulator resulting in the finding that while current generations of devices lack the hardware requirements, in terms of both qubit count and connectivity, to solve realistic problems, future generation devices may be highly efficient.

category: bioinformatics


title: Rice NIN-LIKE PROTEIN 3 plays a significant role in nitrogen use efficiency and grain yield under nitrate-sufficient conditions

abstract: Nitrogen (N) is an essential macronutrient for crop growth and yield, but excessive application of N fertilizer has caused serious environmental pollution and greatly increased the cost of agricultural production. One of the effective and economical solutions to this dilemma is to improve the N use efficiency (NUE) of crops. Although some components involved in regulating NUE have been identified, the underlying molecular mechanism remains largely elusive in rice. Here we report that the OsNLP3 (NIN-like protein 3) is an important regulator in NUE and grain yield under nitrate-sufficient conditions. Both NUE and grain yield were significantly improved by enhanced OsNLP3 expression in the field, but reduced in osnlp3 mutants. The expression of OsNLP3 responds to both nitrate and ammonium, while OsNLP3 nuclear retention is only induced by nitrate, not by ammonium. OsNLP3 regulates the expression of a set of genes involved in N transport and assimilation by directly binding to the nitrate-responsive cis-element in the promoters of these genes. Our study demonstrates that OsNLP3 is significant for the regulation of NUE and grain yield, particularly in nitrate-rich conditions, thus providing a candidate for improving NUE and grain yield in rice.


title: Competitive coordination of the dual roles of the Hedgehog co-receptor in homophilic adhesion and signal reception

abstract:
Hedgehog (Hh) signaling patterns embryonic tissues and contributes to homeostasis in adults. In Drosophila, Hh transport and signaling are thought to occur along a specialized class of actin-rich filopodia, termed cytonemes. Here, we report that Interference hedgehog (Ihog) not only forms a Hh receptor complex with Patched to mediate intracellular signaling, but Ihog also engages in trans-homophilic binding leading to cytoneme stabilization in a manner independent of its role as the Hh receptor. Both functions of Ihog (trans-homophilic binding for cytoneme stabilization and Hh binding for ligand sensing) involve a region of the first fibronectin repeat of the extracellular domain. Thus, the Ihog-Ihog interaction and the Hh-Ihog interaction cannot occur simultaneously for a single Ihog molecule. By combining experimental data and mathematical modeling, we determined that Hh-Ihog heterophilic interaction dominates and Hh can disrupt and displace Ihog molecules involved in trans-homophilic binding. Consequently, we proposed that the weaker Ihog-Ihog trans interaction promotes and stabilizes direct membrane contacts along cytonemes and that, as the cytoneme encounters secreted Hh ligands, the ligands trigger release of Ihog from trans Ihog-Ihog complex enabling transport or internalization of the Hh ligand-Ihog-Patched -receptor complex. Thus, the seemingly incompatible functions of Ihog in homophilic adhesion and ligand binding cooperate to assist Hh transport and reception along the cytonemes.

One-sentence summaryApparently incompatible functions of Ihog in homophilic adhesion and ligand binding cooperate in Hh transport and reception.


title: Impaired local intrinsic immunity to SARS-CoV-2 infection in severe COVID-19

abstract: Infection with SARS-CoV-2, the virus that causes COVID-19, can lead to severe lower respiratory illness including pneumonia and acute respiratory distress syndrome, which can result in profound morbidity and mortality. However, many infected individuals are either asymptomatic or have isolated upper respiratory symptoms, which suggests that the upper airways represent the initial site of viral infection, and that some individuals are able to largely constrain viral pathology to the nasal and oropharyngeal tissues. Which cell types in the human nasopharynx are the primary targets of SARS-CoV-2 infection, and how infection influences the cellular organization of the respiratory epithelium remains incompletely understood. Here, we present nasopharyngeal samples from a cohort of 35 individuals with COVID-19, representing a wide spectrum of disease states from ambulatory to critically ill, as well as 23 healthy and intubated patients without COVID-19. Using standard nasopharyngeal swabs, we collected viable cells and performed single-cell RNA-sequencing (scRNA-seq), simultaneously profiling both host and viral RNA. We find that following infection with SARS-CoV-2, the upper respiratory epithelium undergoes massive reorganization: secretory cells diversify and expand, and mature epithelial cells are preferentially lost. Further, we observe evidence for deuterosomal cell and immature ciliated cell expansion, potentially representing active repopulation of lost ciliated cells through coupled secretory cell differentiation. Epithelial cells from participants with mild/moderate COVID-19 show extensive induction of genes associated with anti-viral and type I interferon responses. In contrast, cells from participants with severe lower respiratory symptoms appear globally muted in their anti-viral capacity, despite substantially higher local inflammatory myeloid populations and equivalent nasal viral loads. This suggests an essential role for intrinsic, local epithelial immunity in curbing and constraining viral-induced pathology. Using a custom computational pipeline, we characterized cell-associated SARS-CoV-2 RNA and identified rare cells with RNA intermediates strongly suggestive of active replication. Both within and across individuals, we find remarkable diversity and heterogeneity among SARS-CoV-2 RNA+ host cells, including developing/immature and interferon-responsive ciliated cells, KRT13+ "hillock"-like cells, and unique subsets of secretory, goblet, and squamous cells. Finally, SARS-CoV-2 RNA+ cells, as compared to uninfected bystanders, are enriched for genes involved in susceptibility (e.g., CTSL, TMPRSS2) or response (e.g., MX1, IFITM3, EIF2AK2) to infection. Together, this work defines both protective and detrimental host responses to SARS-CoV-2, determines the direct viral targets of infection, and suggests that failed anti-viral epithelial immunity in the nasal mucosa may underlie the progression to severe COVID-19.

category: immunology


title: Predicted epitope binding core motifs and common amino acid polymorphisms distinguish protective- from susceptibility-associated HLA alleles in leishmaniasis

abstract: Susceptibility for leishmaniasis is largely dependent on genetic- and immune factors of the host. Despite the previously described association of human leukocyte antigen (HLA) gene cluster variants as genetic susceptibility factors, little is known on the mechanisms that mediate these associations. To characterize the functionality underpinning these associations between HLA and disease, we predicted the epitope binding repertoires for all known leishmaniasis-associated HLA variants collected in a thorough literature review. We identified several amino acid polymorphisms in the HLA sequences that distinguished protective-from risk-associated HLA-DRB1 alleles. Proteome-wide and multi-species T cell epitope binding predictions were carried out across these alleles, enabling us to map the effects on the epitope binding repertoires. The protective-associated HLA-DRB1 alleles were characterized by common binding core motifs, which map to the identified amino acid polymorphisms. These results strongly suggest that polymorphism in the HLA region, resulting in differential antigen presentation, affects the association between HLA and leishmaniasis disease development. Finally, we established a valuable open-access resource of putative epitopes, of which a set of 14 HLA-unrestricted strong-binding epitopes, conserved across species, were prioritized for further epitope discovery in the search for novel subunit-based vaccines.


title: Determining the pathogenicity of MUT gene variant by mini-gene splicing assay

abstract:
ObjectiveTo perform gene mutation splicing analysis on 12 suspected pathogenic MMUT/MUT gene mutation sites to verify the pathogenicity of the mutation.

MethodsWild-type and mutant minigenes were inserted into pcMINI vector, and total 5 wild-type recombinant vectors and 12 mutant recombinant vectors were constructed. The total RNA in 293T cells was extracted after the recombinant vectors were transfected into 293T cell line. Then PCR products were detected by agarose gel electrophoresis and analyzed by sequencing.

ResultsRT-PCR and sequencing results showed that among 12 mutations, 9 mutations (c.419T>C, 469G>T, c.470T>A, c.626dupC, c.693C>G, c.976A>G, c.1009T>C, c.1777G>T and c.1874A>C) did not affect the gene splicing, and the other 3 mutations (c.454C>T, c.421G>A and c.2125-3C>G) all affected mRNA splicing.

ConclusionIn recent case reports of MMUT/MUT gene mutation sites, variant of uncertain significance (VUS) variation is very common. In this study, the pathogenicity of three mutation sites is confirmed by the mini-gene method.


title: Improving stroke detection via Correlation Tensor MRI

abstract: Noninvasively detecting and characterizing modulations in cellular scale micro-architecture is a desideratum for contemporary neuroimaging. Diffusion MRI (dMRI) has become the mainstay methodology for probing microstructure, and, in ischemia, its contrasts have revolutionized stroke management. However, the biological underpinnings of the contrasts observed in conventional dMRI in general and in ischemia in particular are still highly debated since the markers only indirectly reporter on microstructure. Here, we present Correlation Tensor MRI (CTI), a method that rather than measuring diffusion, harnesses diffusion correlations as its source of contrast. We show that CTI can resolve the sources of diffusional kurtosis, which in turn, provide dramatically enhanced specificity and sensitivity towards ischemia. In particular, the sensitivity towards ischemia nearly doubles, both in grey matter (GM) and white matter (WM), and unique signatures for neurite beading, cell swelling, and edema are inferred from CTI. The enhanced sensitivity and specificity endowed by CTI bodes well for future applications in biomedicine, basic neuroscience, and in the clinic.


title: Unraveling micro-architectural modulations in neural tissue upon ischemia by Correlation Tensor MRI


title: Finding Druggable Sites in Proteins using TACTICS

abstract: Structure-based drug discovery efforts require knowledge of where drug-binding sites are located on target proteins. To address the challenge of finding druggable sites, we developed a machine-learning algorithm called TACTICS (Trajectory-based Analysis of Conformations To Identify Cryptic Sites), which uses an ensemble of molecular structures (such as molecular dynamics simulation data) as input. First, TACTICS uses k-means clustering to select a small number of conformations that represent the overall conformational heterogeneity of the data. Then, TACTICS uses a random forest model to identify potentially bindable residues in each selected conformation, based on protein motion and geometry. Lastly, residues in possible binding pockets are scored using fragment docking. As proof-of-principle, TACTICS was applied to the analysis of simulations of the SARS-CoV-2 main protease and methyltransferase and the Yersinia pestis aryl carrier protein. Our approach recapitulates known small-molecule binding sites and predicts the locations of sites not previously observed in experimentally determined structures. The TACTICS code is available at https://github.com/Albert-Lau-Lab/tactics_protein_analysis.


title: Structural insight into host plasma membrane association and assembly of HIV-1 Matrix protein

abstract: HIV-1 continues to be a global health concern since AIDS was first recognized by the World Health Organization (WHO). It is estimated that there were 38 million people infected with HIV-1 and 1.5 million deaths in 2019 alone. A better understanding of the details of the HIV late-stage life cycle, involving Pr55Gag attachment to the membrane for the further oligomerization to release virion, will provide us new avenues for potential treatment. Inositol hexakisphosphate (IP6) is an abundant endogenous cyclitol molecule and its binding was linked to the oligomerization of Pr55Gag via the MA domain. However, the binding site of IP6 on MA was unknown and the structural details of this interaction were missing. Here, we present three high-resolution crystal structures of the MA domain in complex with IP6 molecules to reveal its binding mode. Additionally, extensive Differential Scanning Fluorimetry analysis combined with cryo- and ambient-temperature X-ray crystallography and computational biology identify the key residues that participate in IP6 binding. Our data provide novel insights about the multilayered HIV-1 virion assembly process that involves the interplay of IP6 with PIP2, a phosphoinositide essential for the membrane binding of Pr55Gag. IP6 and PIP2 have neighboring alternate binding sites within the same highly basic region (residues 18-33). This indicates that IP6 and PIP2 bindings are not mutually exclusive and may play a key role in coordinating virion particles membrane localization. Based on our three different IP6-MA complex crystal structures, we propose a new model that involves the IP6 coordination of the oligomerization of outer MA and inner CA domain 2D layers during assembly and budding.


title: BTBD9 is a novel component of IGF signaling and regulates manganese-induced dopaminergic dysfunction

abstract: Restless legs syndrome (RLS) is a common neurological disorder associated with iron deficiency and dopaminergic (DAergic) neuronal dysfunction. BTBD9 is a genetic risk factor for RLS. However, its molecular function remains largely unknown. Here, we report the interaction between BTBD9, manganese (Mn) and insulin/insulin-like growth factor (IGF) signaling in Caenorhabditis elegans, mouse Neuro2a cells and humans. We found that elevated Mn downregulated BTBD9 mRNA levels; in turn, BTBD9 expression attenuated Mn-induced cellular stress and dopaminergic neurodegeneration. As Mn is a known co-factor for insulin receptor and IGF-1 receptor, which activates IGF signaling, we posited that BTBD9 negatively regulates IGF signaling. Our results showed that the protective effects of BTBD9 against Mn toxicity were dependent on the forkhead box O (FOXO) protein. Furthermore, BTBD9 overexpression significantly elevated FOXO level and decreased PKB level, while phosphoinositide-dependent kinase-1 (PDK1) level remained unchanged. We conclude that BTBD9 acts as a key component in the IGF signaling pathway. Meanwhile, the roles of Mn in DAergic neurotoxicity and regulating BTBD9 shed new light on the etiology of RLS.


title: Causal approach to environmental risks of seabed mining

abstract: Seabed mining is approaching the commercial mining phase across the worlds oceans. This rapid industrialization of seabed resource use is introducing new pressures to marine environments. The environmental impacts of such pressures should be carefully evaluated prior to permitting new activities, yet observational data is mostly missing. Here, we examine the environmental risks of seabed mining using a causal, probabilistic network approach. Drawing on a series of interviews with a multidisciplinary group of experts, we outline the cause-effect pathways related to seabed mining activities to inform quantitative risk assessments. The approach consists of (1) iterative model building with experts to identify the causal connections between seabed mining activities and the affected ecosystem components, and (2) quantitative probabilistic modelling to provide estimates of mortality of benthic fauna in the Baltic Sea. The model is used to evaluate alternative mining scenarios, offering a quantitative means to highlight the uncertainties around the impacts of mining. We further outline requirements for operationalizing quantitative risk assessments, highlighting the importance of a cross-disciplinary approach to risk identification. The model can be used to support permitting processes by providing a more comprehensive description of the potential environmental impacts of seabed resource use, allowing iterative updating of the model as new information becomes available.


title: GAT-FD: An Integrated MATLAB Toolbox for Graph Theoretical Analysis of Task-Related Functional Dynamics

abstract: Functional connectivity (FC) has been demonstrated to be varying over time during sensory and cognitive processes. Quantitative examinations of such variations can significantly advance our understanding on large-scale functional organizations and their topological dynamics that support normal brain functional connectome and can be altered in individuals with brain disorders. However, toolboxes that integrate the complete functions for analyzing task-related brain FC, functional network topological properties, and their dynamics, are still lacking. The current study has developed a MATLAB toolbox, the Graph Theoretical Analysis of Task-Related Functional Dynamics (GAT-FD), which consists of four modules for sliding-window analyses, temporal mask generation, estimations of network properties and dynamics, and result display, respectively. All the involved functions have been tested and validated using fMRI data collected from human subjects when performing a block-designed task. The results demonstrated that the GAT-FD allows for effective and quantitative evaluations of the functional network properties and their dynamics during the task period. As an open-source and user-friendly package, the GAT-FD and its detailed user manual are freely available at https://www.nitrc.org/projects/gat_fd and https://centers.njit.edu/cnnl/gat_fd/.


title: Genome-wide transcriptional responses of marine nematode Litoditis marina to hyposaline and hypersaline stresses

abstract: Maintenance of osmotic homeostasis is essential for all organisms, especially for marine animals in the ocean with 30{per thousand} salinity or higher. However, the underlying molecular mechanisms that how marine animals adapt to high salinity environment compared to their terrestrial relatives, remain elusive. Here, we investigated marine animals genome-wide transcriptional responses to salinity stresses using an emerging marine nematode model Litoditis marina. We found that the transthyretin-like family genes were significantly increased in both hyposaline and hypersaline conditions, while multiple neurotransmitter receptor and ion transporter genes were down-regulated in both conditions, suggesting the existence of conserved strategies for response to stressful salinity environments in L. marina. Unsaturated fatty acids biosynthesis related genes, neuronal related tubulins and intraflagellar transport genes were specifically up-regulated in hyposaline treated worms, while exhibited the opposite regulation in hypersaline condition. By contrast, cuticle related collagen genes were enriched and up-regulated for hypersaline response, interestingly, the expression of these collagen genes was significantly decreased in hyposaline condition. Given a wide range of salinity tolerance of the marine nematodes, this study and further genetic analysis of key gene(s) of osmoregulation in L. marina will likely provide important insights into biological evolution and environmental adaptation mechanisms in nematodes and other invertebrate animals in general.