title: massNet: integrated processing and classification of spatially resolved mass spectrometry data using deep learning for rapid tumor delineation

abstract:
MotivationMass spectrometry imaging (MSI) provides rich biochemical information in a label-free manner and therefore holds promise to substantially impact current practice in disease diagnosis. However, the complex nature of MSI data poses computational challenges in its analysis. The complexity of the data arises from its large size, high dimensionality, and spectral non-linearity. Preprocessing, including peak picking, has been used to reduce raw data complexity, however peak picking is sensitive to parameter selection that, perhaps prematurely, shapes the downstream analysis for tissue classification and ensuing biological interpretation.

ResultsWe propose a deep learning model, massNet, that provides the desired qualities of scalability, non-linearity, and speed in MSI data analysis. This deep learning model was used, without prior preprocessing and peak picking, to classify MSI data from a mouse brain harboring a patient-derived tumor. The massNet architecture established automatically learning of predictive features, and automated methods were incorporated to identify peaks with potential for tumor delineation. The models performance was assessed using cross-validation, and the results demonstrate higher accuracy and a 174-fold gain in speed compared to the established classical machine learning method, support vector machine.

Availability and ImplementationThe code is publicly available on GitHub.


title: Restoring vascular endothelial autophagic flux reduces atherosclerotic lesions

abstract: Atherosclerotic lesions preferentially develop in arterial areas exposed to low shear stress, where endothelial cells express a pro-inflammatory, apoptotic, and senescent phenotype. Autophagy is a lysosomal mechanism that recycles damaged organelles and protein aggregates to maintain cellular homeostasis. Stimulation of autophagy in high shear stress conditions is an atheroprotective process. Conversely, endothelial cells exposed to atheroprone low shear stress present a defective autophagic flux, which favors a pro-inflammatory phenotype and the formation of atherosclerotic lesions. Since an efficient autophagic flux is dependent on -tubulin acetylation, which is reduced under low shear stress, we hypothesized that increasing -tubulin acetylation could restore adequate levels of autophagy in endothelial cells exposed to low shear stress. We found that blocking Histone Deacetylase 6 (HDAC6) activity, either by pharmacological inhibition (Tubastatin-A) or genetic approaches (shHDAC6), raised levels of acetylated -tubulin, as well as LC3-II/I ratio, LC3 punctae area and autophagic flux in cultured endothelial cells exposed to low shear stress. This effect was associated with a reduced expression of inflammatory markers (Intercellular adhesion molecule-1 (ICAM-1), Vascular cell Adhesion Protein-1 (VCAM-1) and Monocyte Chemoattractant Protein-1 (MCP-1)) in Tumor Necrosis Factor-alpha (TNF-)-stimulated cells. We observed increased endothelial autophagic flux in the aortic arch of the HDAC6-/-/ApoE-/- mice. Subsequently, atherosclerotic plaque size was significantly reduced in the atheroprone areas of chimeric HDAC6-/-/ApoE-/- mice, transplanted with HDAC6+/+/ApoE-/- bone marrow, when compared to HDAC6+/+/ApoE-/- littermate controls. Taken together, these results indicate that targeting -tubulin acetylation, via HDAC6-inhibition, may be an interesting strategy to restore endothelial autophagic flux and to promote an atheroprotective endothelial phenotype despite unfavorable shear stress conditions.

category: pathology


title: Cardioluminescence in transgenic zebrafish embryos: a calcium imaging tool to study drug effects and pathological modeling

abstract: The zebrafish embryo has emerged as an excellent model in cardiovascular research. The existing techniques to monitor Ca2+ in the heart based on fluorescent Ca2+ biosensors are limited due to phototoxicity and photobleaching. To overcome these issues, we have used bioluminescence. We generated a transgenic line expressing GFP-Aequorin in the heart, Tg(cmlc2:GA), and optimized an in vivo aequorin reconstitution protocol to improve the luminescence capacity. This allowed imaging Ca2+ in long duration recordings in embryos of 3 to 5 days post-fertilization. The analogs diacetyl h-coelenterazine and f-coelenterazine enhanced the light output and signal-to-noise ratio from the embryos. With this cardioluminescence model, we monitored the time-averaged Ca2+ levels and beat-to-beat Ca2+ oscillations. Changes in Ca2+ levels were observed by incubation with BayK8644, an L-type Ca2+ channel agonist, the channel blocker nifedipine, and {beta}-adrenergic blocker propranolol. Treatment of zebrafish embryos with terfenadine for 24 hours has been proposed as a model of heart failure. Tg(cmlc2:GA) embryos treated with terfenadine showed a 2:1 atrioventricular block and a decrease in the ventricular Ca2+ levels.

category: physiology


title: Relationships of diversity and evenness in adaptation strategies of the effect of protective coloration of animals

abstract: This study aims to investigate the relationship between the diversity and uniformity of the protective coloration of animals. The protective coloration of animals is the result of the formation of adaptation strategies. The aim of this study is to identify the colorimetric parameters that characterize the adaptive properties of the protective coloration of animals in relation to their natural habitat. The relevance of this study is related to the importance of the tasks of remote collection of information about various biological species; in particular, animals that carry dangerous infections. This is all the more important in light of the current epidemic problems facing humanity. The revealed colorimetric characteristics can be used to describe digital images and their further recognition. They can be classifying features in machine learning algorithms. We give an example of digital image processing of a duck (Anas platyrhynchos). The analysis of the distribution of the found colorimetric parameter makes it possible to identify the image segments corresponding to the Anas platyrhynchos. The results thus obtained are also of interest from the point of view of the problem of the adaptive role of the diversity of biological systems in relation to the problem of the mechanisms of functioning of the protective coloration of animals.

category: systems biology


title: Host adaptation through hybridization: Genome analysis of triticale powdery mildew reveals unique combination of lineage-specific effectors

abstract: The emergence of new fungal pathogens through hybridization represents a serious challenge for agriculture. Hybridization between the wheat mildew (Blumeria graminis f.sp. tritici) and rye mildew (B.g. f.sp. secalis) pathogens have led to the emergence of a new mildew form (B.g. f.sp. triticale) growing on triticale, a man-made amphiploid crop derived from crossing rye and wheat which was originally resistant to the powdery mildew disease. The identification of the genetic basis of host-adaptation in triticale mildew has been hampered by the lack of a reference genome. Here we report the 141.4 Mb reference assembly of B.g. triticale isolate THUN-12 derived from long-read sequencing and genetic map-based scaffolding. All eleven B.g. triticale chromosomes were assembled from telomere-to-telomere and revealed that 19.7% of the hybrid genome was inherited from the rye mildew parental lineage. We identified lineage-specific regions in the hybrid, inherited from the rye or wheat mildew parental lineages, that harbour numerous bona fide candidate effectors. We propose that the combination of lineage-specific effectors in the hybrid genome is crucial for host-adaptation, allowing the fungus to simultaneously circumvent the immune systems contributed by wheat and rye in the triticale crop. In line with this we demonstrate the functional transfer of the SvrPm3 effector from wheat to triticale mildew, a virulence effector that specifically suppresses resistance of the wheat Pm3 allelic series. This transfer is the likely underlying cause for the observed poor effectiveness of several Pm3 alleles against triticale mildew and exemplifies the negative implications of pathogen hybridizations on resistance breeding.

category: genomics


title: Phylodynamic insights on the early spread of the COVID-19 pandemic and the efficacy of intervention measures

abstract:
We performed phylodynamic analyses of all available SARS-CoV-2 genomes from the early phase of the COVID-19 pandemic--combined with a novel dataset on contemporary global air-travel volume--to assess the efficacy of public-health measures on viral geographic spread. Globally, viral dispersal rates are significantly correlated with air-travel volume, and widespread international air-travel bans imposed against China by early February coincide with a significant reduction in geographic viral spread. In North America, the efficacy of this travel ban was temporary, possibly due to the lack of both containment measures against other infected regions and domestic mitigation measures. By contrast, in China, domestic mitigation measures were correlated with a long-term reduction in viral spread, despite repeated international introductions. Our study supports a role for both targeted international containment and domestic mitigation measures as critical components of a more comprehensive public-health strategy to mitigate future outbreaks caused by the emergence of novel SARS-CoV-2 variants.

One sentence summaryPhylodynamic analyses reveal that variation in rates of early geographic spread of COVID-19 are correlated with intervention measures.

category: evolutionary biology


title: Single Cell Atlas of Human Putamen Reveals Disease Specific Changes in Synucleinopathies: Parkinson's Disease and Multiple System Atrophy

abstract:
Understanding disease biology at a cellular level from disease specific tissues is imperative for effective drug development for complex neurodegenerative diseases. We profiled 87,086 nuclei from putamen tissue of healthy controls, Parkinsons Disease (PD), and Multiple System Atrophy (MSA) subjects to construct a comprehensive single cell atlas. Although both PD and MSA are manifestations of -synuclein protein aggregation, we observed that both the diseases have distinct cell-type specific changes. We see a possible expansion and activation of microglia and astrocytes in PD compared to MSA and controls. Contrary to PD microglia, we found absence of upregulated unfolded protein response in MSA microglia compared to controls. Differentially expressed genes in major cell types are enriched for genes associated with PD-GWAS loci. We found altered expression of major neurodegeneration associated genes -- SNCA, MAPT, LRRK2, and APP -- at cell-type resolution. We also identified disease associated gene modules using a network biology approach. Overall, this study creates an interactive atlas from synucleinopathies and provides major cell-type specific disease insights.

Link to interactive atlas will be made available at the time of publication.


title: Pharmacologically induced weight loss is associated with distinct gut microbiome changes in obese rats

abstract:
BackgroundObesity, metabolic disease and some psychiatric conditions are associated with changes to relative abundance of bacterial species and specific genes in the fecal microbiome. Little is known about the impact of pharmacologically induced weight loss on distinct gut microbiome species and their respective gene programs in obese individuals.

ResultsUsing shotgun metagenomics, the composition of the microbiome was obtained for two cohorts of obese female Wistar rats (n=10-12, total of 82) maintained on a high fat diet before and after a 42-day treatment with a panel of four investigatory or approved anti-obesity drugs (tacrolimus/FK506, bupropion, naltrexone and sibutramine), alone or in combination. We found that sibutramine treatment induced consistent weight loss through reducing food intake. Decreased weight loss in sibutramine-treated rats was associated with changes to the gut microbiome that included increased beta-diversity, increased Bacteroides/Firmicutes ratio and increased relative abundance of multiple Bacteroides species. In addition, the relative abundance of multiple genes was found to be differentially abundant, including significant reductions in components of flagellum and genes involved in flagellum assembly.

ConclusionsThis study provides a large resource comprising complete shotgun metagenomics datasets of the fecal microbiome coupled with weight change and food intake at day 3, day 15 and day 42 from 82 obese rats treated with a range of compounds used for weight loss, which are available to the community for detailed analysis. Furthermore, by conducting a detailed analysis of the microbiome associated with sibutramine-induced weight loss, we have identified multiple weight-loss associated microbial taxa and pathways. These include a reduction in components of flagellum and the flagellum assembly pathway that points to a potential role of sibutramine-induced weight-loss on regulating bacterially driven anti-inflammatory responses.


title: Squeegee: de-novo identification of reagent and laboratory induced microbial contaminants in low biomass microbiomes

abstract: Computational analysis of host-associated microbiomes has opened the door to numerous discoveries relevant to human health and disease. However, contaminant sequences in metagenomic samples can potentially impact the interpretation of findings reported in microbiome studies, especially in low biomass environments. Our hypothesis is that contamination from DNA extraction kits or sampling lab environments will leave taxonomic "bread crumbs" across multiple distinct sample types, allowing for the detection of microbial contaminants when negative controls are unavailable. To test this hypothesis we implemented Squeegee, a de novo contamination detection tool. We tested Squeegee on simulated and real low biomass metagenomic datasets. On the low biomass samples, we compared Squeegee predictions to experimental negative control data and show that Squeegee accurately recovers known contaminants. We also analyzed 749 metagenomic datasets from the Human Microbiome Project and identified likely previously unreported kit contamination. Collectively, our results highlight that Squeegee can identify microbial contaminants with high precision. Squeegee is open-source and available at: https://gitlab.com/treangenlab/squeegee

category: bioinformatics


title: iProMix: A decomposition model for studying the function of ACE2 based on bulk proteogenomic data for coronavirus pathogenesis

abstract: Both SARS-CoV and SARS-CoV-2 use ACE2 receptors to enter epithelial cells in lung and many other tissues to cause human diseases. Genes and pathways that regulate ACE2 may facilitate/inhibit viral entry and replication, and genes and pathways that are controlled by ACE2 may be perturbed during infection, both affecting disease severity and outcomes. It is critical to understand how genes and pathways are associated with ACE2 in epithelial cells by leveraging proteomic data, but an accurate large-scale proteomic profiling at cellular resolution is not feasible at current stage. Therefore, we propose iProMix, a novel framework that decomposes bulk tissue proteomic data to identify epithelial cell component specific associations between ACE2 and other proteins. Unlike existing decomposition based association analyses, iProMix allows both predictors and outcomes to be impacted by cell type composition of the tissue and accounts for the impacts of decomposition variations and errors on hypothesis tests. It also builds in the functions to improve cell type estimation if estimates from existing literature are unsatisfactory. Simulations demonstrated that iProMix has well-controlled false discovery rate and large power in non-asymptotic settings with both correctly and mis-specified cell-type composition. We applied iProMix to the 110 adjacent normal tissue samples of patients with lung adenocarcinoma from Clinical Proteomic Tumor Analysis Consortium, and identified that interferon  and{gamma} pathways were most significantly associated with ACE2 protein abundances in epithelial cells. Interestingly, the associations were sex-specific that the positive associations were only observed in men, while in women the associations were negative.


title: Generic solving of one-compartment toxicokinetic models

abstract:
This paper gives the full analytical solution of the generic set of ordinary differential equations that define one-compartment toxicokinetic models. These models describe uptake and elimination processes taking place within living organisms when exposed to chemical substances. The models solved in this paper consider living organisms as a unique compartment, into which a parent compound enters via several possible exposure routes and from which it is eliminated as well as its potential metabolites. Benefiting from generic solutions of one-compartment toxicokinetic models is particularly useful when fitting them to experimental data, facilitating the writing of the inference algorithms leading to parameter estimates. Additionally, these models are of crucial interest in environmental risk assessment for the calculation of bioaccumulation metrics as required by regulators in support of decision making when they evaluate dossiers for marketing authorisation of active substances.

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=129 SRC="FIGDIR/small/442956v2_ufig1.gif" ALT="Figure 1">
View larger version (18K):
org.highwire.dtl.DTLVardef@1e8381org.highwire.dtl.DTLVardef@14b957forg.highwire.dtl.DTLVardef@143dc14org.highwire.dtl.DTLVardef@25c478_HPS_FORMAT_FIGEXP  M_FIG O_FLOATNOGraphical abstractC_FLOATNO C_FIG

category: pharmacology and toxicology


title: A hybrid PDE-ABM model for viral dynamics with application to SARS-CoV-2 and influenza

abstract:
We propose a hybrid partial differential equation - agent-based (PDE-ABM) model to describe the spatio-temporal viral dynamics in a cell population. The virus concentration is considered as a continuous variable and virus movement is modelled by diffusion, while changes in the states of cells (i.e. healthy, infected, dead) are represented by a stochastic agent-based model. The two subsystems are intertwined: the probability of an agent getting infected in the ABM depends on the local viral concentration, and the source term of viral production in the PDE is determined by the cells that are infected.

We develop a computational tool that allows us to study the hybrid system and the generated spatial patterns in detail. We systematically compare the outputs with a classical ODE system of viral dynamics, and find that the ODE model is a good approximation only if the diffusion coefficient is large.

We demonstrate that the model is able to predict SARS-CoV-2 infection dynamics, and replicate the output of in vitro experiments. Applying the model to influenza as well, we can gain insight into why the outcomes of these two infections are different.


title: Multi-site enzymes as a mechanism for bistability in reaction networks

abstract: Here, we focus on a common class of enzymes that have multiple substrate-binding sites (multi-site enzymes), and analyse their capacity to generate bistable dynamics in the reaction systems that they are embedded in. Using mathematical techniques, we show that the inherent binding and catalysis reactions arising from multiple substrate-enzyme complexes creates a potential for bistable dynamics in a reaction system. We construct a generic model of an enzyme with n substrate binding sites and derive an analytical solution for the steady state concentration of all enzyme-substrate complexes. By studying these expressions, we obtain a mechanistic understanding for bistability and derive parameter combinations that guarantee bistability and show how changing specific enzyme kinetic parameters and enzyme levels can lead to bistability in reaction systems involving mjulti-site enzymes. Thus, the presented findings provide a biochemical and mathematical basis for predicting and engineering bistability in multi-site enzymes.


title: Independent dynamics of slow, intermediate, and fast intracranial EEG spectral activities during human memory formation

abstract: A wide spectrum of brain rhythms are engaged throughout the human cortex in cognitive functions. How the rhythms of various low and high frequencies are spatiotemporally coordinated across the human brain during memory processing is inconclusive. They can either be coordinated together across a wide range of the frequency spectrum or induced in specific bands. We used a large dataset of human intracranial electroencephalography (iEEG) to parse the spatiotemporal dynamics of spectral activities induced during formation of verbal memories. Encoding of words for subsequent free recall activated slow theta, intermediate alpha and beta, and fast gamma frequency power in discrete cortical sites. A majority of the electrode sites recorded activity in only one of these frequencies, except for the visual cortex where spectral power was induced across multiple bands. Each frequency band showed characteristic dynamics of the induced power specific to cortical area and hemisphere. The power of the low, intermediate, and fast activities propagated in distinct spatiotemporal patterns across the visual, temporal and prefrontal cortical areas as the words were presented for encoding. Our results suggest anatomically and temporally distributed spectral activities in the formation of human memory.

category: neuroscience


title: Subpopulations of neurons in lOFC encode previous and current rewards at time of choice

abstract: 1Studies of neural dynamics in lateral orbitofrontal cortex (lOFC) have shown that subsets of neurons that encode distinct aspects of behavior, such as value, may project to common downstream targets. However, it is unclear whether reward history, which may subserve lOFCs well-documented role in learning, is represented by functional subpopulations in lOFC. Previously, we analyzed neural recordings from rats performing a value-based decision-making task, and we documented trial-by-trial learning that required lOFC (Constantinople et al., 2019). Here we characterize functional subpopulations of lOFC neurons during behavior, including their encoding of task variables. We found five distinct clusters of lOFC neurons, either based on clustering of their trial-averaged peristimulus time histograms (PSTHs), or a feature space defined by their average conditional firing rates aligned to different task variables. We observed weak encoding of reward attributes, but stronger encoding of reward history, the animals left or right choice, and reward receipt across all clusters. Only one cluster, however, encoded the animals reward history at the time shortly preceding the choice, suggesting a possible role in integrating previous and current trial outcomes at the time of choice. This cluster also exhibits qualitatively similar responses to identified corticostriatal projection neurons in a recent study (Hirokawa et al., 2019), and suggests a possible role for subpopulations of lOFC neurons in mediating trial-by-trial learning.


title: The influence of the respiratory cycle on reaction times in sensory-cognitive paradigms

abstract: Behavioural and electrophysiological studies point to widespread influence of the state of respiration on brain activity and cognitive performance. Still, the prevalence and relevance of such respiratory-behavioural relations in typical sensory-cognitive tasks remain unclear. We here used a battery of six tasks probing sensory detection, discrimination and short-term memory to address the questions of whether and by how much behaviour covaries with the respiratory cycle. Our results show that participants tended to align their respiratory cycle to the experimental paradigm. Furthermore, their reaction times, but not so much their response accuracy, consistently and significantly covaried with the respiratory cycle, and this effect was strongest when analyzed contingent on the respiratory state at participants responses. The respective effect sizes where comparable to those seen in many typical neurocognitive experimental manipulations. These results support a prominent relation between respiration and sensory-cognitive function and suggest that sensation is intricately linked to rhythmic bodily or interoceptive functions.


title: Capturing the nature of events and event context using Hierarchical Event Descriptors (HED)

abstract: Event-related data analysis plays a central role in EEG and MEG (MEEG) and other neuroimaging modalities such as fMRI. Choices about which events to report and how to annotate their full natures significantly influence the value, reliability, and reproducibility of neuroimaging datasets for further analysis and meta- or mega-analysis. A powerful annotation strategy using the new third-generation formulation of the Hierarchical Event Descriptors (HED) framework and tools (hedtags.org) combines robust event description with details of experiment design and metadata in a human-readable as well as machine-actionable form, making event annotation relevant to the full range of neuroimaging and other time series data. This paper considers the event design and annotation process using as a case study the well-known multi-subject, multimodal dataset of Wakeman and Henson made available by its authors as a Brain Imaging Data Structure (BIDS) dataset (bids.neuroimaging.io). We propose a set of best practices and guidelines for event annotation integrated in a natural way into the BIDS metadata file architecture, examine the impact of event design decisions, and provide a working example of organizing events in MEEG and other neuroimaging data. We demonstrate how annotations using HED can document events occurring during neuroimaging experiments as well as their interrelationships, providing machine-actionable annotation enabling automated within- and across-experiment analysis and comparisons. We discuss the evolution of HED software tools and have made an accompanying HED-annotated BIDS-formatted edition of the MEEG data of the Wakeman and Henson dataset (openneuro.org, ds003645).


title: Theta phase mediates deliberate action switching in human Supplementary Motor Areas

abstract: The ability to deliberately overwrite ongoing automatic actions is a necessary feature of adaptive behavior. It has been proposed that the supplementary motor areas (SMAs) operate as a controller that orchestrates the switching between automatic and deliberate processes by inhibiting ongoing behaviors and so facilitating the execution of alternative ones. In addition, previous studies support the involvement of SMAs theta waves (4-9 Hz) in cognitive control. However, the exact role of such oscillatory dynamics and their contribution to the control of action are not fully understood. To investigate the mechanisms by which the SMAs support direct control of deliberate behavior, we recorded intracranial electroencephalography (iEEG) activity in humans performing a motor sequence task. Subjects had to perform a "change of plans" motor task requiring habitual movements to be overwritten at unpredictable moments. We found that SMAs were exclusively active during trials that demand action reprogramming in response to the unexpected cue but were silent during automatic action execution. Importantly, SMAs activity was characterized by a distinct temporal pattern, expressed in a stereotypical phase alignment of theta oscillations. More specifically, single trial motor performance was correlated with the trial contribution to the global inter-trial phase coherence, with higher coherence associated with faster trials. In addition, theta phase modulated the amplitude of gamma oscillations, with higher cross-frequency coupling in faster trials. Our results suggest that within frontal cortical networks, theta oscillations could encode a control signal that promotes the execution of deliberate actions.


title: Infection, recovery and re-infection of farmed mink with SARS-CoV-2

abstract:
Mink, on a farm with about 15,000 animals, became infected with SARS-CoV-2. Over 75% of tested animals were positive for SARS-CoV-2 RNA in throat swabs and 100% of tested animals were seropositive. The virus responsible had a deletion of nucleotides encoding residues H69 and V70 within the spike protein gene. The infected mink recovered and after free-testing of the mink, the animals remained seropositive. During follow-up studies, after a period of more than 2 months without virus detection, over 75% of tested animals scored positive again for SARS-CoV-2 RNA. Whole genome sequencing showed that the virus circulating during this re-infection was most closely related to the virus identified in the first outbreak on this farm but additional sequence changes had occurred. Animals had much higher levels of anti-SARS-CoV-2 antibodies after re-infection than at free-testing. Thus, following recovery from an initial infection, seropositive mink rapidly became susceptible to re-infection by SARS-CoV-2.

Article Summary LineFollowing widespread infection with SARS-CoV-2 of mink on a farm, all tested animals had seroconverted and the farm was then tested free of infection; however, less than 3 months later, a further round of infection affected more than 75% of tested animals.

category: microbiology


title: Subtle differences in the pathogenicity of SARS-CoV-2 variants of concern B.1.1.7 and B.1.351 in rhesus macaques

abstract: The emergence of several SARS-CoV-2 variants has caused global concerns about increased transmissibility, increased pathogenicity, and decreased efficacy of medical countermeasures. Animal models can be used to assess phenotypical changes in the absence of confounding factors that affect observed pathogenicity and transmissibility data in the human population. Here, we studied the pathogenicity of variants of concern (VOC) B.1.1.7 and B.1.351 in rhesus macaques and compared it to a recent clade B.1 SARS-CoV-2 isolate containing the D614G substitution in the spike protein. The B.1.1.7 VOC behaved similarly to the D614G with respect to clinical disease, virus shedding and virus replication in the respiratory tract. Inoculation with the B.1.351 isolate resulted in lower clinical scores in rhesus macaques that correlated with lower virus titers in the lungs, less severe histologic lung lesions and less viral antigen detected in the lungs. We observed differences in the local innate immune response to infection. In bronchoalveolar lavages, cytokines and chemokines were upregulated on day 4 in animals inoculated with D614G and B.1.1.7 but not in those inoculated with B.1.351. In nasal samples, we did not detect upregulation of cytokines and chemokines in D614G or B.1.351-inoculated animals. However, cytokines and chemokines were upregulated in the noses of B.1.1.7-inoculated animals. Taken together, our comparative pathogenicity study suggests that ongoing circulation under diverse evolutionary pressure favors transmissibility and immune evasion rather than an increase in intrinsic pathogenicity.


title: iCodon: ideal codon design for customized gene expression

abstract: Messenger RNA (mRNA) stability substantially impacts steady-state gene expression levels in a cell. mRNA stability, in turn, is strongly affected by codon composition in a translation-dependent manner across species, through a mechanism termed codon optimality. We have developed iCodon (www.iCodon.org), an algorithm for customizing mRNA expression through the introduction of synonymous codon substitutions into the coding sequence. iCodon is optimized for four vertebrate transcriptomes: mouse, human, frog, and fish. Users can predict the mRNA stability of any coding sequence based on its codon composition and subsequently generate more stable (optimized) or unstable (deoptimized) variants encoding for the same protein. Further, we show that codon optimality predictions correlate with expression levels using fluorescent reporters and endogenous genes in human cells and zebrafish embryos. Therefore, iCodon will benefit basic biological research, as well as a wide range of applications for biotechnology and biomedicine.

category: molecular biology


title: GLOBAL LAND-USE AND LAND-COVER DATA FOR ECOLOGISTS: HISTORICAL, CURRENT AND FUTURE SCENARIOS

abstract: Land-use land-cover (LULC) data are important predictors of species occurrence and biodiversity threat. Although there are LULC datasets available under current conditions, there is a lack of such data under historical and future climatic conditions. This hinders, for example, projecting niche and distribution models under global change scenarios at different time scenarios. The Land Use Harmonization Project (LUH2) is a global terrestrial dataset at 0.25{degrees} spatial resolution that provides LULC data from 850 to 2300 for 12 LULC state classes. The dataset, however, is compressed in a file format (NetCDF) that is incompatible for many analyses and requires layer extractions and transformations that are intractable for most researchers. Here we selected and transformed the LUH2 in a standard GIS format data to make it more user-friendly. We provide LULC for every year from 850 to 2100, and from 2015 on, the LULC dataset is provided under two Shared Socioeconomic Pathways (SSP2-4.5 and SSP5-8.5). We provide two types of files for each year: separate files with continuous values for each of the 12 LULC state classes, and a single categorical file with all state classes combined. To create the categorical layer, we assigned the state with the highest value in a given pixel among the 12 continuous data. LUH2 predicts a pronounced decrease in primary forest, particularly noticeable in the Amazon, the Brazilian Atlantic Forest, the Congo Basin and the boreal forests, an equally pronounced increase in secondary forest and non-forest lands, and in croplands in the Brazilian Atlantic Forest and sub-Saharan Africa. The final dataset provides LULC data for 1251 years that will be of interest for macroecology, ecological niche modeling, global change analysis, and other applications in ecology and conservation.

category: ecology


title: GLOBAL LAND-USE AND LAND-COVER DATA: HISTORICAL, CURRENT AND FUTURE SCENARIOS


title: Single Nuclei Sequencing of Human Putamen Oligodendrocytes Reveals Altered Heterogeneity and Disease-Associated Changes in Parkinson's Disease and Multiple System Atrophy

abstract: The role of oligodendrocytes in neurodegenerative diseases remains incompletely understood and largely unexplored at the single cell level. We profiled 87,086 single nuclei from human brain putamen region for healthy control, Parkinsons Disease (PD), and Multiple System Atrophy (MSA). Oligodendrocyte lineage cells were the dominant cell-type in the putamen with oligodendrocyte subpopulations clustered by transcriptomic variation found to exhibit diverse functional enrichment patterns, and this oligodendrocyte heterogeneity was altered in a disease-specific way. Among profiled oligodendrocyte subpopulations, differences in expression of SNCA, HAPLN2, MAPT, APP, and OPALIN were observed for PD and MSA compared with healthy controls. Intriguingly, greater activation of unfolded protein response pathway gene expression was observed in PD nuclei versus MSA. Using network analysis, we then identified specific PD- and MSA-correlated gene co-expression modules enriched with disease relevant pathways; the PD-correlated module was significantly enriched for Parkinsons Disease GWAS loci (p = 0.01046). Our analysis provides a broader understanding of oligodendrocyte heterogeneity and reveals distinctive oligodendrocyte pathological alterations associated with PD and MSA which may suggest potential novel therapeutic targets and new strategies for disease modification.


title: A phylogeny-based metric for estimating changes in transmissibility from recurrent mutations in SARS-CoV-2

abstract:
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged in late 2019 and spread globally to cause the COVID-19 pandemic. Despite the constant accumulation of genetic variation in the SARS-CoV-2 population, there was little evidence for the emergence of significantly more transmissible lineages in the first half of 2020. Starting around November 2020, several more contagious and possibly more virulent  Variants of Concern (VoCs) were reported in various regions of the world. These VoCs share some mutations and deletions that haven arisen recurrently in distinct genetic backgrounds. Here, we build on our previous work modelling the association of mutations to SARS-CoV-2 transmissibility and characterise the contribution of individual recurrent mutations and deletions to estimated viral transmissibility. We then assess how patterns of estimated transmissibility in all SARS-CoV-2 clades have varied over the course of the COVID-19 pandemic by summing transmissibility estimates for all individual mutations carried by any sequenced genome analysed. Such an approach recovers the Delta variant (21A) as the most transmissible clade currently in circulation, followed by the Alpha variant (20I). By assessing transmissibility over the time of sampling, we observe a tendency for estimated transmissibility within clades to slightly decrease over time in most clades. Although subtle, this pattern is consistent with the expectation of a decay in transmissibility in mainly non-recombining lineages caused by the accumulation of weakly deleterious mutations. SARS-CoV-2 remains a highly transmissible pathogen, though such a trend could conceivably play a role in the turnover of different global viral clades observed over the pandemic so far.

CaveatsO_LIThis work is not about the severity of disease. We do not analyse the severity of disease. We do not present any evidence that SARS-CoV-2 has decreased in severity.
C_LIO_LILineage replacement dynamics are affected by many factors. The trend we recover for a decrease in inferred transmissibility of a clade over time is a small effect. We caution against over-interpretation. This result would not affect the management of the SARS-CoV-2 pandemic: for example, we make no claims about any impact on the efficacy of particular non-pharmaceutical interventions (NPIs).
C_LIO_LIOur phylogeny-based method to infer changes in estimated transmissibility due to recurrent mutations and deletions makes a number of simplifying assumptions. These may not all be valid. The consistent trend for the slight decrease we report might be due to an as-yet-unidentified systematic bias.
C_LI


title: Self-assembly of tessellated tissue sheets by growth and collision

abstract: Tissues do not exist in isolation-they interact with other tissues within and across organs. While cell-cell interactions have been intensely investigated, less is known about tissue-tissue interactions. Here, we studied collisions between monolayer tissues with different geometries, cell densities, and cell types. First, we determine rules for tissue shape changes during binary collisions and describe complex cell migration at tri-tissue boundaries. Next, we demonstrate that genetically identical tissues displace each other based solely on cell density gradients, and present a physical model of tissue interactions that allows us to estimate the bulk modulus of the tissues from collision dynamics. Finally, we introduce TissEllate, a design tool for self-assembling complex tessellations from arrays of many tissues, and we use cell sheet engineering techniques to transfer these composite tissues like cellular films. Overall, our work provides insight into the mechanics of tissue collisions, harnessing them to engineer tissue composites as designable living materials.

category: bioengineering


title: Accurate prediction of protein torsion angles using evolutionary signatures and recurrent neural network

abstract: The amino acid sequence of a protein contains all the necessary information to specify its shape, which dictates its biological activities. However, it is challenging and expensive to experimentally determine the three-dimensional structure of proteins. The backbone torsion angles, as an important structural constraint, play a critical role in protein structure prediction, and accurately predicting the angles can considerably advance the tertiary structure prediction by accelerating efficient sampling of the large conformational space for low energy structures. On account of the rapid growth of protein databases and striking breakthroughs in deep learning algorithms, computational advances allow us to extract knowledge from large-scale data to address key biological questions. Here we propose evolutionary signatures that are computed from protein sequence profiles, and a deep neural network, termed ESIDEN, that adopts a straightforward architecture of recurrent neural networks with a small number of learnable parameters. The proposed ESIDEN is validated on three benchmark datasets, including D2020, TEST2016/2018, and CASPs datasets. On the D2020, using the combination of the four novel features and basic features, the ESIDEN achieves the mean absolute error (MAE) of 15.7 and 19.8 for{phi} and{psi} , respectively. Comparing to the best-so-far methods, we show that the ESIDEN significantly improves the angle{psi} by the MAE decrements of more than 3.5 degrees on both TEST2016 and TEST2018 and achieves better MAE of the angle{phi} by decrements of at least 0.3 degrees although it adopts simple architecture and fewer learnable parameters. On fifty-nine template-free modeling targets, the ESIDEN achieves high accuracy by reducing the MAEs by 0.6 and more than 2.3 degrees on average for the torsion angles{phi} and{psi} in the CASPs, respectively. Using the predicted torsion angles, we infer the tertiary structures of four representative template-free modeling targets that achieve high precision with regard to the root-mean-square deviation and TM-score by comparing them to the native structures. The results demonstrate that the ESIDEN can make accurate predictions of the torsion angles by leveraging the evolutionary signatures. The proposed evolutionary signatures would be also used as alternative features in predicting residue-residue distance, protein structure, and protein-ligand binding sites. Moreover, the high-precision torsion angles predicted by the ESIDEN can be used to accurately infer protein tertiary structures, and the ESIDEN would potentially pave the way to improve protein structure prediction.


title: Unilateral Resection of Both Cortical Visual Pathways Alters Action but not Perception

abstract: The human cortical visual system consists of two major pathways, a ventral pathway that subserves perception and a dorsal pathway that subserves visuomotor control. These pathways follow dissociable developmental trajectories, and, accordingly, might be differentially susceptible to neurodevelopmental disorders or injuries. Previous studies have found that children with cortical resections of the ventral visual pathway retain largely normal visuoperceptual abilities. Whether visually guided actions, supported by computations carried out by the dorsal pathway, follow a similar pattern remains unknown. To address this question, we examined visuoperceptual and visuomotor behaviors in a pediatric patient, TC, who underwent a cortical resection that included portions of the left ventral and dorsal pathways. We collected data when TC used her right and left hands to perceptually estimate the width blocks that varied in width and length, and, separately, to grasp the same blocks. TCs perceptual estimation performance was comparable to that of controls, independent of the hand used. In contrast, relative to controls, she showed reduced visuomotor sensitivity to object shape and this was more evident when she grasped the objects with her contralesional right hand. These results provide evidence for a striking difference in the reorganization profiles of the two visual pathways. This difference supports the notion that the two pathways exhibit differential susceptibility to neurodevelopmental disorders.


title: Multiple mechanisms for overcoming lethal over-initiation of DNA replication

abstract: DNA replication is a highly regulated process that is primarily controlled at the step of initiation. In the gram-positive bacterium Bacillus subtilis the replication initiator DnaA, is regulated by YabA, which inhibits cooperative binding at the origin. Mutants lacking YabA have increased and asynchronous initiation. We found that under conditions of rapid growth, the dnaA1 mutation that causes replication over-initiation, was synthetic lethal with a deletion of yabA. We isolated several classes of suppressors of the lethal phenotype of the {Delta}yabA dnaA1 double mutant. Some suppressors (dnaC, cshA) caused a decrease in replication initiation. Others (relA, nrdR) stimulate replication elongation. One class of suppressors decreased levels of the replicative helicase, DnaC, thereby limiting replication initiation. We found that decreased levels of helicase were sufficient to decrease replication initiation under fast growth conditions. Our results highlight the multiple mechanisms cells use to regulate DNA replication.


title: Palbociclib-mediated cell cycle arrest can occur in the absence of the CDK inhibitors p21 and p27

abstract: The use of CDK4/6 inhibitors in the treatment of a wide range of cancers is an area of ongoing investigation. Despite their increasing clinical use, there is limited understanding of the determinants of sensitivity and resistance to these drugs. Recent data has cast doubt on how CDK4/6 inhibitors arrest proliferation, provoking renewed interest in the role(s) of CDK4/6 in driving cell proliferation. As the use of CDK4/6 inhibitors in cancer therapies becomes more prominent, an understanding of their effect on the cell cycle becomes more urgent. Here, we investigate the mechanism of action of CDK4/6 inhibitors in promoting cell cycle arrest. Two main models explain how CDK4/6 inhibitors cause G1 cell cycle arrest, which differ in their dependence on the CDK inhibitor proteins p21 and p27. We have used live and fixed single-cell quantitative imaging, with inducible degradation systems, to address the roles of p21 and p27 in the mechanism of action of CDK4/6 inhibitors. We find that CDK4/6 inhibitors can initiate and maintain a cell cycle arrest without p21 or p27. This work clarifies our current understanding of the mechanism of action of CDK4/6 inhibitors and has implications for cancer treatment and patient stratification.

category: cell biology


title: Metagenome-genome-wide association studies reveal human genetic impact on the oral microbiome

abstract: The oral microbiota contains billions of microbial cells, which could contribute to diseases in a number of body sites. Challenged by eating, drinking and dental hygiene on a daily basis, the oral microbiota is regarded as highly dynamic. Here, we report significant human genomic associations with the oral metagenome from more than 1,915 individuals, for both the tongue dorsum and saliva. We identified five genetic loci associated with oral microbiota at study-wide significance (p < 3.16 x 10-11). Four of the five associations were well replicated in an independent cohort of 1,439 individuals: rs1196764 at APPL2 with Prevotella jejuni, Oribacterium uSGB 3339 and Solobacterium uSGB 315; rs3775944 at the serum uric acid transporter SLC2A9 with Oribacterium uSGB 1215, Oribacterium uSGB 489 and Lachnoanaerobaculum umeaense; rs4911713 near OR11H1 with species F0422 uSGB 392; and rs36186689 at LOC105371703 with Eggerthia. Further analyses confirmed 84% (386/455 for tongue dorsum) and 85% (391/466 for saliva) of genetic-microbiota associations including 6 genome-wide significant associations mutually validated between the two niches. Human genetics accounted for at least 10% of oral microbiome compositions between individuals. Machine learning models showed that polygenetic risk score dominated over oral microbiome in predicting predisposing risk of dental diseases such as dental calculus and gingival bleeding. These findings indicate that human genetic differences are one explanation for a stable or recurrent oral microbiome in each individual.

category: genetics


title: Ozone sensitivity of diverse maize genotypes is associated with differences in gene regulation, not gene content

abstract: The maize pangenome has demonstrate large amounts of presence/absence variation and it has been hypothesized that presence/absence variation contributes to stress response. To uncover whether the observed genetic variation in physiological response to elevated ozone (a secondary air pollutant that causes significant crop yield losses) concentration is due to variation in genic content, and/or variation in gene expression, we examine the impact of sustained elevated ozone concentration on the leaf tissue from 5 diverse maize inbred genotypes (B73, Mo17, Hp301, C123, NC338). Analysis of long reads from the transcriptomes of the 10 conditions found expressed genes in the leaf are part of the shared genome, with 94.5% of expressed genes from syntenic loci. Quantitative analysis of short reads from 120 plants (twelve from each condition) found limited transcriptional response to sustained ozone stress in the ozone resistant B73 genotype (151 genes), while more than 3,300 genes were significantly differentially expressed in the more sensitive NC338 genotype. The genes underpinning the divergence of B73 from the other 4 genotypes implicates ethylene signaling consistent with some findings in Arabidopsis. For the 82 of the 83 genes differentially expressed among all 5 genotypes and the 788 of 789 genes differentially expressed in 4 genotypes (excluding B73) in sensitivity to ozone is associated with oxidative stress tolerance being associated with a weaker response to a reactive oxygen species (ROS) signal and suggests that genetic variation in downstream processes is key to ozone tolerance.


title: NERD-seq: A novel approach of Nanopore direct RNA sequencing that expands representation of non-coding RNAs

abstract: The new next-generation sequencing platforms by Oxford Nanopore Technologies for direct RNA sequencing (direct RNA-seq) allow for an in-depth and comprehensive study of the epitranscriptome by enabling direct base calling of RNA modifications. Non-coding RNAs constitute the most frequently documented targets for RNA modifications. However, the current standard direct RNA-seq approach is unable to detect many of these RNAs. Here we present NERD-seq, a sequencing approach which enables the detection of multiple classes of non-coding RNAs excluded by the current standard approach. Using total RNA from a tissue with high known transcriptional and non-coding RNA activity in mouse, the brain hippocampus, we show that, in addition to detecting polyadenylated coding and non-coding transcripts as the standard approach does, NERD-seq is able to significantly expand the representation for other classes of RNAs such as snoRNAs, snRNAs, scRNAs, srpRNAs, tRNAs, rRFs and non-coding RNAs originating from LINE L1 elements. Thus, NERD-seq presents a new comprehensive direct RNA-seq approach for the study of epitranscriptomes in brain tissues and beyond.


title: Ca2+-mediated higher-order assembly of b0,+AT-rBAT is a key step for system b0,+ biogenesis and cystinuria

abstract: Cystinuria is a genetic disorder characterized by overexcretion of dibasic amino acids and cystine, which causes recurrent kidney stones and occasionally severe kidney failure. Mutations of the two responsible proteins, rBAT and b0,+AT, which comprise system b0,+, are linked to type I and non-type I cystinuria respectively and they exhibit distinct phenotypes due to protein trafficking defects or catalytic inactivation. Although recent structural insights into human b0,+AT-rBAT suggested a model for transport-inactivating mutations, the mechanisms by which type I mutations trigger trafficking deficiencies are not well understood. Here, using electron cryo-microscopy and biochemistry, we discover that Ca2+-mediated higher-order assembly of system b0,+ is the key to its trafficking on the cell surface. We show that Ca2+ stabilizes the interface between two rBAT molecules to mediate super-dimerization, and this in turn facilitates the N-glycan maturation of system b0,+. A common cystinuria mutant T216M and mutations that disrupt the Ca2+ site in rBAT cause the loss of higher-order assemblies, resulting in protein trafficking deficiency. Mutations at the super-dimer interface reproduce the mis-trafficking phenotype, demonstrating that super-dimerization is essential for cellular function. Cell-based transport assays confirmed the importance of the Ca2+ site and super-dimerization, and additionally suggested which residues are involved in cationic amino acid recognition. Taken together, our results provide the molecular basis of type I cystinuria and serve as a guide to develop new therapeutic strategies against it. More broadly, our findings reveal an unprecedented link between transporter oligomeric assembly and trafficking diseases in general.

category: biochemistry


title: Endothelial restoration of CAD GWAS gene PLPP3 by nanomedicine suppresses YAP/TAZ activity and reduces atherosclerosis in vivo

abstract: Genome-wide association studies (GWAS) have suggested new molecular mechanisms in vascular cells driving atherosclerotic diseases such as coronary artery disease (CAD) and ischemic stroke (IS). Nevertheless, a major challenge to develop new therapeutic approaches is to spatiotemporally manipulate these GWAS-identified genes in specific vascular tissues in vivo. YAP (Yes-associated protein) and TAZ (transcriptional coactivator with PDZ-binding motif) have merged as critical transcriptional regulators in cells responding to biomechanical stimuli, such as in athero-susceptible endothelial cells activated by disturbed flow (DF). The molecular mechanisms by which DF activates while unidirectional flow (UF) inactivates YAP/TAZ remain incompletely understood. Recent studies demonstrated that DF and genetic predisposition (risk allele) of CAD/IS locus 1p32.2 converge to reduce phospholipid phosphatase 3 (PLPP3) expression in vascular endothelium. Restoration of endothelial PLPP3 in vivo, although remains challenging and unexplored, is hypothesized to reduce atherosclerosis. We devised a nanomedicine system integrating nanoparticles and Cdh5 promoter-driven plasmids to successfully restore PLPP3 expression in activated endothelium, resulting in suppressed YAP/TAZ activity and reduced DF-induced atherosclerosis in mice. Mechanistically, our studies discovered a molecular paradigm by which CAD/IS GWAS gene PLPP3 inactivates YAP/TAZ by reducing lysophosphatidic acid (LPA)-induced myosin II and ROCK in endothelium under UF. These results highlight a new mechanistic link between GWAS and YAP/TAZ mechano-regulation and moreover, establish a proof of concept of vascular wall-based therapies employing targeted nanomedicine to manipulate CAD/IS GWAS genes in vivo.


title: Phase separation in fluids with many interacting components

abstract: Fluids in natural systems, like the cytoplasm of a cell, often contain thousands of molecular species that are organized into multiple coexisting phases that enable diverse and specific functions. How interactions between numerous molecular species encode for various emergent phases is not well understood. Here we leverage approaches from random matrix theory and statistical physics to describe the emergent phase behavior of fluid mixtures with many species whose interactions are drawn randomly from an underlying distribution. Through numerical simulation and stability analyses, we show that these mixtures exhibit staged phase separation kinetics and are characterized by multiple coexisting phases at equilibrium with distinct compositions. Random-matrix theory predicts the number of existing phases at equilibrium, validated by simulations with diverse component numbers and interaction parameters. Surprisingly, this model predicts an upper bound on the number of phases, derived from dynamical considerations, that is much lower than the limit from the Gibbs phase rule, which is obtained from equilibrium thermodynamic constraints. Using a biophysically motivated model of pairwise interactions between components, we design ensembles that encode either linear or non-monotonic scaling relationships between number of components and co-existing phases, which we validate through simulation and theory. Finally, inspired by parallels in biological systems, we show that including nonequilibrium turnover of components through chemical reactions can tunably modulate the number of coexisting phases at steady-state without changing overall fluid composition. Together, our study provides a model framework that describes the emergent dynamical and steady-state phase behavior of liquid-like mixtures with many interacting constituents.

category: biophysics


title: Bergamottin a CYP3A inhibitor found in grapefruit juice inhibits prostate cancer cell growth by downregulating androgen receptor signaling causing cell cycle block and apoptosis

abstract: Prostate cancer is the second leading cause of cancer related death in American men. Several therapies have been developed to treat advanced prostate cancer, but these therapies often have severe side effects. To improve the outcome with fewer side effects we focused on the furanocoumarin bergamottin, a natural product found in grapefruit juice and a potent CYP3A inhibitor. Our recent studies have shown that CYP3A5 inhibition can block androgen receptor (AR) signaling, critical for prostate cancer growth. We observed that bergamottin reduces prostate cancer (PC) cell growth by decreasing both total and nuclear AR (AR activation) reducing downstream AR signaling. Bergamottins role in reducing AR activation was confirmed by confocal microscopy studies and reduction in PSA levels. Further studies revealed that bergamottin promotes cell cycle block and accumulates G0/G1 cells. The cell cycle block was accompanied with reduction in cyclin D, cyclin B, CDK4, P-cdc2 (Y15) and P-wee1 (S642). We also observed that bergamottin triggers apoptosis in prostate cancer cell lines as evident by TUNEL staining and PARP cleavage. Our data suggest that bergamottin may be used as an adjunctive nutritional supplement to suppress prostate cancer growth and is of relevance to AA patients carrying wild type CYP3A5 often presenting aggressive disease.

category: cancer biology


title: Tools for efficient analysis of neurons in a 3D reference atlas of whole mouse spinal cord

abstract: Spinal neurons are highly heterogeneous in location, transcriptional identity and function. To understand their contributions to sensorimotor circuits, it is essential to map the positions of identified subsets of neurons in relation to others throughout the spinal cord (SC), but we lack tools for whole SC sample preparation, imaging and in toto analysis. To overcome this problem, we have (1) designed scaffolds (SpineRacks) that facilitate efficient and ordered cryo-sectioning of the entire SC in a single block, (2) constructed a 3D reference atlas of adult mouse SC and (3) developed software (SpinalJ) to register images of sections and for standardized analysis of cells and projections in atlas space. We have verified mapping accuracies for known neurons and demonstrated the usefulness of this platform to reveal unknown neuronal distributions. Together, these tools provide high-throughput analyses of whole mouse SC and enable direct comparison of 3D spatial information between animals and studies.


title: Ongoing neural oscillations predict the post-stimulus outcome of closed loop auditory stimulation during slow-wave sleep.

abstract:
The large slow oscillation (SO, 0.5-2Hz) that characterises slow-wave sleep is crucial to memory consolidation and other physiological functions. Manipulating slow oscillations can enhance sleep and memory, as well as benefitting the immune system. Closed-loop auditory stimulation (CLAS) has been demonstrated to increase the SO amplitude and to boost fast sleep spindle activity (11-16Hz). Nevertheless, not all such stimuli are effective in evoking SOs, even if they are precisely phase-locked. Here, we studied whether it is possible to use ongoing activity patterns to determine which oscillations to stimulate in order to effectively enhance SOs or SO-locked spindle activity. To this end, we trained classifiers using the morphological characteristics of the ongoing SO, as measured by electroencephalography (EEG), to predict whether stimulation would lead to a benefit in terms of the resulting SO and spindle amplitude. Separate classifiers were trained using trials from spontaneous control and stimulated datasets, and we evaluated their performance by applying them to held-out data both within and across conditions. We were able to predict both when large SOs will occur spontaneously, and whether a phase-locked auditory click will effectively enlarge them with an accuracy of ~70%. We were also able to predict when stimulation would elicit spindle activity with an accuracy of ~60%. Finally, we evaluate the importance of the various SO features used to make these predictions. Our results offer new insight into SO and spindle dynamics and provide a new method for online optimisation of stimulation.

HIGHLIGHTS- Random forest classifiers can predict spontaneous and stimulated SOs and spindle amplitudes.
- Morphological wave features predicted the response of SOs and spindles to CLAS.
- SO amplitude during the click is the main predictor for post-stimulus SO amplitude.
- Prediction of spindle activity did not differ in accuracy for stimulated vs spontaneous data.


title: Role of hypocretin in the medial preoptic area in the regulation of sleep, maternal behavior and body temperature of lactating rats.

abstract: The hypocretins (HCRT), also known as orexin, includes two neuroexcitatory peptides, HCRT-1 and HCRT-2 (orexin A y B, respectively), synthesized by neurons located in the postero-lateral hypothalamus, whose projections and receptors are widely distributed throughout the brain, including the medial preoptic area (mPOA). HCRT have been associated with a wide range of physiological functions including sleep-wake cycle, maternal behavior and body temperature, all regulated by the mPOA. Previously we showed that HCRT in the mPOA facilitates certain active maternal behaviors, while the blockade of HCRT-R1 increased the time spent in nursing. As mother rats mainly sleep while they nurse, we hypothesize that HCRT in the mPOA of lactating rats reduce sleep and nursing, while the intra-mPOA administration of the dual orexin receptor antagonist (DORA) would generate the opposite effect. Therefore, the aim of this study was to determine the role of HCRT within the mPOA, in the regulation and integration of the sleep-wake cycle, maternal behavior and body temperature of lactating rats. To evaluate this idea, we assessed the sleep-wake states, maternal behavior and body temperature of lactating rats following microinjections of HCRT-1 (100 and 200 M) and DORA (5mM) into the mPOA. As expected, our data shows that HCRT-1 in mPOA promoted wakefulness and a slightly increase in body temperature, whereas DORA increased both NREM and REM sleep along with nursing and milk ejection. Taken together, our results strongly suggest that the reduction of the endogenous HCRT within the mPOA of lactating rats is important to promote sleep, nursing and milk ejection.


title: Discoidin Domain Receptor 2 regulates AT1 receptor expression in Angiotensin II-stimulated cardiac fibroblasts via fibronectin-dependent Integrin-β1 signalling

abstract: Recent reports on the cardioprotective effects of fibronectin inhibition following myocardial injury suggest a largely unexplored role for the extracellular matrix (ECM) glycoprotein in cardiac fibroblast function. We probed the molecular regulation and functional implications of fibronectin gene expression in cardiac fibroblasts exposed to Angiotensin II, a major pro-fibrotic factor in the myocardium. Using gene knockdown and over-expression approaches, western blotting and promoter pull-down assay, we show that collagen type I-activated Discoidin Domain Receptor 2 (DDR2) mediates Angiotensin II-stimulated transcriptional up-regulation of fibronectin by Yes-activated Protein in cardiac fibroblasts. Further, siRNA-mediated fibronectin knockdown attenuated Angiotensin II-stimulated expression of collagen type I and anti-apoptotic cIAP2, and enhanced susceptibility to apoptosis. Importantly, an obligate role for fibronectin was observed in Angiotensin II-stimulated expression of AT1R, the Angiotensin II receptor, which would link ECM signaling and Angiotensin II signaling in cardiac fibroblasts. Moreover, conditioned medium from DDR2- or fibronectin-silenced cardiac fibroblasts reduced AT1R expression in H9c2 cardiomyoblasts. The regulatory role of fibronectin in Angiotensin II-stimulated cIAP2, collagen type I and AT1R expression was mediated by Integrin-{beta}1-integrin-linked kinase signaling. In vivo, we observed modestly reduced basal levels of AT1R in DDR2-null mouse myocardium, associated with the previously reported reduction in myocardial Integrin-{beta}1 levels. The role of fibronectin, downstream of DDR2, could be a critical determinant of cardiac fibroblast-mediated wound healing following myocardial injury. In summary, our findings suggest a complex mechanism of regulation of cardiac fibroblast function involving two major extracellular matrix proteins, collagen type I and fibronectin, and their receptors, DDR2 and Integrin-{beta}1.


title: Discoidin Domain Receptor 2 regulates AT1R expression in Angiotensin II-stimulated cardiac fibroblasts via fibronectin-dependent Integrin-β1 signaling


title: Meta-analysis and multi-omics to elucidate pathogenic mechanisms of age-related knee osteoarthritis

abstract: Increased mechanistic insight into the pathogenesis of knee osteoarthritis (KOA) is needed to develop efficacious disease-modifying treatments. Though age-related pathogenic mechanisms are most relevant to the great majority of KOA seen clinically, the bulk of our mechanistic understanding of KOA has been derived using surgically induced post-traumatic OA (PTOA) models. Here, we took an integrated approach of meta-analysis and multi-omics to elucidate pathogenic mechanisms of age-related KOA in murine model. Protein-level data together with transcriptomic profiling revealed inflammation, autophagy, and cellular senescence as primary hallmarks of age-related KOA. Importantly, the molecular profiles of aged cartilage were unique from those in PTOA, with only 1% overlap between the two. At the nexus of aging hallmarks, Advanced Glycation End-Product (AGE)/Receptor for AGE emerged as intrinsically linked to age-related KOA. This pathway was further validated by mass spectrometry. Collectively, these findings implicate dysregulation of AGE-RAGE signaling as a key driver of age-related KOA.


title: poreCov - an easy to use, fast, and robust workflow for SARS-CoV-2 genome reconstruction via nanopore sequencing

abstract: In response to the SARS-CoV-2 pandemic, a highly increased sequencing effort has been established worldwide to track and trace ongoing viral evolution. Technologies such as nanopore sequencing via the ARTIC protocol are used to reliably generate genomes from raw sequencing data as a crucial base for molecular surveillance. However, for many labs that perform SARS-CoV-2 sequencing, bioinformatics is still a major bottleneck, especially if hundreds of samples need to be processed in a recurring fashion. Pipelines developed for short-read data cannot be applied to nanopore data. Therefore, specific long-read tools and parameter settings need to be orchestrated to enable accurate genotyping and robust reference-based genome reconstruction of SARS-CoV-2 genomes from nanopore data. Here we present poreCov, a highly parallel workflow written in Nextflow, using containers to wrap all the tools necessary for a routine SARS-CoV-2 sequencing lab into one program. The ease of installation, combined with concise summary reports that clearly highlight all relevant information, enables rapid and reliable analysis of hundreds of SARS-CoV-2 raw sequence data sets or genomes. poreCov is freely available on GitHub under the GNUv3 license: github.com/replikation/poreCov.


title: Quantifying Central Metabolic Fluxes in Human Platelets Using Metabolic Flux Analysis

abstract: Platelet metabolism is linked to platelet hyper- and hypoactivity in numerous human diseases. Developing a detailed understanding of the link between metabolic shifts and platelet activation state is integral to improving human health. Here, we show the first application of isotopically nonstationary 13C metabolic flux analysis to quantitatively measure carbon fluxes in both resting and thrombin activated platelets. Resting platelets primarily metabolize glucose to lactate via glycolysis, while acetate is oxidized to fuel the tricarboxylic acid cycle. Upon activation with thrombin, a potent platelet agonist, platelets increase their uptake of glucose 3-fold. This results in an absolute increase in flux throughout central metabolism, but when compared to resting platelets they redistribute carbon dramatically. Activated platelets decrease relative flux to the oxidative pentose phosphate pathway and TCA cycle from glucose and increase relative flux to lactate. These results provide the first report of reaction-level carbon fluxes in platelets and allow us to distinguish metabolic fluxes with much higher resolution than previous studies.


title: Isotopically Nonstationary 13C Metabolic Flux Analysis in Resting and Activated Human Platelets


title: Reversion of antibiotic resistance in drug-resistant bacteria using non-steroidal anti-inflammatory drug benzydamine

abstract: Antimicrobial resistance has been a growing concern that gradually undermines our tradition treatment regimen. The fact that few antibacterial drugs with new scaffolds or targets have been approved in the past two decades aggravates this crisis. Repurposing previously approved drugs as potent antibiotic adjuvants offers a cost-effective strategy to mitigate the development of resistance and tackle the increasing infections by multidrug-resistant (MDR) bacteria. Herein, we found that benzydamine, a widely used non-steroidal anti-inflammatory drug in clinic, remarkably potentiated broad-spectrum antibiotic-tetracyclines activity against a panel of clinical important resistant pathogens, including MRSA, VRE, MCRPEC and tet(X)-positive Gram-negative bacteria. Further mechanistically experiments showed that benzydamine dissipated membrane potential ({Delta}{Psi}) in both Gram-positive and negative bacteria, which in turn upregulated the transmembrane proton gradient ({Delta}pH) and promoted the uptake of tetracyclines. Additionally, benzydamine exacerbated the oxidative stress by triggering the production of ROS and suppressing GAD system-mediated oxidative defensive. This mode of action explains the great bactericidal activity of the doxycycline-benzydamine combination against different metabolic states of bacteria including persister cells. As a proof-of-concept, the in vivo efficacy of this combination therapy was evidenced in multiple animal infection models. These findings revealed that benzydamine is a promising tetracycline antibiotics adjuvant and has the potential to address life-threatening infections by MDR bacteria.


title: Complete and high-quality genomes of novel microbial species from a meromictic lake using a workflow combining long- and short-read sequencing platforms

abstract:
BackgroundMost of Earths bacteria have yet to be cultivated. The metabolic and functional potentials of these uncultivated microorganisms thus remain mysterious, and the metagenome-assembled genome (MAG) approach is the most robust method for uncovering these potentials. However, MAGs discovered by conventional metagenomic assembly and binning methods are usually highly fragmented genomes with heterogeneous sequence contamination, and this affects the accuracy and sensitivity of genomic analyses. Though the maturation of long-read sequencing technologies provides a good opportunity to fix the problem of highly fragmented MAGs as mentioned above, the methods error-prone nature causes severe problems of long-read-alone metagenomics. Hence, methods are urgently needed to retrieve MAGs by a combination of both long- and short-read technologies to advance genome-centric metagenomics.

ResultsIn this study, we combined Illumina and Nanopore data to develop a new workflow to reconstruct 233 MAGs--six novel bacterial orders, 20 families, 66 genera, and 154 species--from Lake Shunet, a secluded meromictic lake in Siberia. Those new MAGs were underrepresented or undetectable in other MAGs studies using metagenomes from human or other common organisms or habitats. Using this newly developed workflow and strategy, the average N50 of reconstructed MAGs greatly increased 10-40-fold compared to when the conventional Illumina assembly and binning method were used. More importantly, six complete MAGs were recovered from our datasets, five of which belong to novel species. We used these as examples to demonstrate many novel and intriguing genomic characteristics discovered in these newly complete genomes and proved the importance of high-quality complete MAGs in microbial genomics and metagenomics studies.

ConclusionsThe results show that it is feasible to apply our workflow with a few additional long reads to recover numerous complete and high-quality MAGs from short-read metagenomes of high microbial diversity environment samples. The unique features we identified from five complete genomes highlight the robustness of this method in genome-centric metagenomic research. The recovery of 154 novel species MAGs from a rarely explored lake greatly expands the current bacterial genome encyclopedia and broadens our knowledge by adding new genomic characteristics of bacteria. It demonstrates a strong need to recover MAGs from diverse unexplored habitats in the search for microbial dark matter.


title: Hepatitis B virus polymerase restricts LINE-1 retrotransposition

abstract: Long interspersed element-1 (LINE-1, L1) retrotransposon composes about 17% of the human genome. However, genetic and biochemical interactions between L1 and hepatitis B virus (HBV) remain poorly understood. In this study, we found that HBV restricts L1 mobility without inhibiting the L1 promoter activity. Notably, HBV polymerase (Pol) strongly inhibited L1 retrotransposition in a reverse transcriptase (RT)-independent manner. Indeed, the ribonuclease H (RNase H) domain was essential for inhibition of L1 retrotransposition. L1 ORF1p RNA-binding protein predominantly localized into cytoplasmic RNA granule termed P-body. However, HBV Pol sequestered L1 ORF1p from P-body and colocalized with L1 ORF1p in cytoplasm, when both proteins were co-expressed. Altogether, HBV Pol seems to restrict L1 mobility through a sequestration of L1 ORF1p from P-body. Thus, these results suggest a novel function or activity of HBV Pol in regulation of L1 retrotransposition.


title: Inhibitory concentrations of ciprofloxacin induce an adaptive response promoting the intracellular survival of Salmonella Typhimurium

abstract:
Antimicrobial resistance (AMR) is a pressing global health crisis, which has been fuelled by the sustained use of certain classes of antimicrobials, including fluoroquinolones. While the genetic mutations responsible for decreased fluoroquinolone (ciprofloxacin) susceptibility are known, the implications of ciprofloxacin exposure on bacterial growth, survival, and interactions with host cells are not well described. Aiming to understand the influence of inhibitory concentrations of ciprofloxacin in vitro, we subjected three clinical isolates of S. Typhimurium to differing concentrations of ciprofloxacin, dependent on their minimum inhibitory concentrations (MIC), and assessed the impact on bacterial growth, morphology, and transcription. We further investigated the differential morphology and transcription that occurred following ciprofloxacin exposure and measured the ability of ciprofloxacin-treated bacteria to invade and replicate in host cells. We found that ciprofloxacin-exposed S. Typhimurium are able to recover from inhibitory concentrations of ciprofloxacin, and that the drug induces specific morphological and transcriptional signatures associated with the bacterial SOS response, DNA repair, and intracellular survival. In addition, ciprofloxacin-treated S. Typhimurium have increased capacity for intracellular replication in comparison to untreated organisms. These data suggest that S. Typhimurium undergoes an adaptive response under ciprofloxacin perturbation that promotes cellular survival, a consequence that may justify more measured use of ciprofloxacin for Salmonella infections. The combination of multiple experimental approaches provides new insights into the collateral effects that ciprofloxacin and other antimicrobials have on invasive bacterial pathogens.

ImportanceAntimicrobial resistance is a critical concern in global health. In particular, there is rising resistance to fluoroquinolones, such as ciprofloxacin, a first-line antimicrobial for many Gram-negative pathogens. We investigated the adaptive response of clinical isolates of Salmonella Typhimurium to ciprofloxacin, finding that the bacteria adapt in short timespans to high concentrations of ciprofloxacin in a way that promotes intracellular survival during early infection. Importantly, by studying three clinically relevant isolates, we were able to show that individual isolates respond differently to ciprofloxacin, and for each isolate, there was a heterogeneous response under ciprofloxacin treatment. The heterogeneity that arises from ciprofloxacin exposure may drive survival and proliferation of Salmonella during treatment and lead to drug resistance.


title: Longitudinal changes in diet cause repeatable and largely reversible shifts in gut microbial communities of laboratory mice and are observed across segments of the entire intestinal tract

abstract: Dietary changes are known to alter the composition of the gut microbiome. However, it is less understood how repeatable and reversible these changes are and how diet-switches affect the microbiota in the various segments of the gastrointestinal tract. Here, a treatment group of conventionally raised laboratory mice is subjected to two periods of western diet (WD) interrupted by a period of standard diet (SD) of the same duration. Beta-diversity analyses show that diet-induced microbiota changes are largely reversible (q = 0.1501; PERMANOVA, weighted-UniFrac comparison of the treatment-SD group to the control-SD group) and repeatable (q = 0.032; PERMANOVA, weighted-UniFrac comparison of both WD treatments). Furthermore, we report that diet switches alter the gut microbiota composition along the length of the intestinal tract in a segment-specific manner, leading to gut segment-specific Firmicutes/Bacteroidota ratios. We identified prevalent and distinct Amplicon Sequencing Variants (ASVs), particularly in genera of the recently described Muribaculaceae, along the gut as well as ASVs that are differentially abundant between segments of treatment and control groups. Overall, this study provides insights into the reversibility of diet-induced microbiota changes and highlights the importance to expand sampling efforts beyond the collections of fecal samples to characterize diet-dependent and segment-specific microbiome differences.


title: Network analysis of large phospho-signalling datasets: application to Plasmodium-erythrocyte interactions

abstract: Phosphorylation based signalling is a complicated and intertwined series of pathways critical to all domains of life. This interconnectivity, though essential to life, makes understanding and decoding the interactions difficult. Large datasets of phosphorylation interactions through the activity of kinases on their numerous effectors are now being generated, however interpretation of the network environment remains challenging. In humans, many phosphorylation interactions have been identified across published works to form the known phosphorylation interaction network. We overlayed phosphorylation datasets onto this network which provided information to each of the connections. To analyse the datasets now mapped into a network, we designed a pathway analysis that uses random walks to identify chains of phosphorylation events occurring much more or much less frequently than expected. This analysis highlights pathways of phosphorylation that work synergistically, providing a rapid interpretation of the most critical pathways in a given dataset. Here we used datasets of human red blood cells infected with the notable stages of Plasmodium falciparum asexual development. The analysis identified several known signalling interactions, and additional interactions which could form the basis of numerous future studies. The network analysis designed here is widely applicable to any comparative phosphorylation dataset across infection and disease and can provide a rapid and reliable analysis to guide validation studies.


title: FliW and CsrA govern flagellin (FliC) synthesis and play pleiotropic roles in virulence and physiology of Clostridioides difficile R20291

abstract:
Clostridioides difficile is a Gram-positive, spore-forming, and toxin-producing anaerobe that can cause nosocomial antibiotic-associated intestinal disease. In C. difficile, the expression of flagellar genes is coupled to toxin gene regulation and bacterial colonization and virulence. The flagellin FliC is responsible for pleiotropic gene regulation during in vivo infection. However, how fliC expression is regulated is unclear. In Bacillus subtilis, flagellin homeostasis and motility are coregulated by flagellar assembly factor FliW, Flagellin Hag (FliC homolog), and CsrA (Carbon storage regulator A), which is referred to as partner-switching mechanism "FliW-CsrA-Hag". In this study, we characterized FliW and CsrA functions by deleting or overexpressing fliW, csrA, and fliW-csrA in C. difficile R20291. We showed that both fliW deletion or csrA overexpression in R20291, and csrA complementation in R20291{Delta}WA (fliW-csrA codeletion) dramatically decreased FliC production, however, fliC gene transcription was unaffected. While suppression of fliC translation by csrA overexpression was mostly relieved when fliW was coexpressed, and no significant difference in FliC production was detected when only fliW was complemented in R20291{Delta}WA. Further, loss of fliW led to increased biofilm formation, cell adhesion, toxin production, and pathogenicity in a mouse model of C. difficile infection (CDI), while fliW-csrA codeletion decreased toxin production and mortality in vivo. Taken together, these data suggest that CsrA negatively modulates fliC expression and FliW indirectly affects fliC expression through inhibition of CsrA post-transcriptional regulation, which seems similar to the "FliW-CsrA-Hag" switch in B. subtilis. Our data also suggest that "FliW-CsrA-fliC/FliC" can regulate many facets of C. difficile R20291 pathogenicity.

IMPORTANCEC. difficile flagellin FliC is associated with toxin gene expression, bacterial colonization and virulence, and is also involved in pleiotropic gene regulation during in vivo infection. However, how fliC expression is regulated remains unclear. In light of "FliW-CsrA-Hag" switch coregulation mechanism reported in B. subtilis, we showed that fliW and csrA play an important role in flagellin synthesis which affects C. difficile motility directly. Our data also suggest that FliW-CsrA-fliC/FliC" can regulate many facets of C. difficile R20291 pathogenicity. These findings further aid us in understanding the virulence regulation in C. difficile.


title: Host bioenergetic parameters reveal cytotoxicity of anti-tuberculosis drugs undetected using conventional viability assays.

abstract: High attrition rates in tuberculosis (TB) drug development have been largely attributed to safety, which is likely due to the use of endpoint assays measuring cell viability to detect drug cytotoxicity. In drug development of cancer, metabolic and neurological disorders, and antibiotics, cytotoxicity is increasingly being assessed using extracellular flux (XF) analysis, which measures cellular bioenergetic metabolism in real-time. Here, we adopt the XF platform to investigate the cytotoxicity of drugs currently used in TB treatment on the bioenergetic metabolism of HepG2 cells, THP-1 macrophages, and human monocyte derived macrophages (hMDM). We found that the XF analysis reveals earlier drug-induced effects on the cells bioenergetic metabolism prior to cell death, measured by conventional viability assays. Furthermore, each cell type has a distinct response to drug treatment, suggesting that more than one cell type should be considered to examine cytotoxicity in TB drug development. Interestingly, chemically unrelated drugs with different modes of action on Mycobacterium tuberculosis have similar effects on the bioenergetic parameters of the cells, thus, discouraging the prediction of potential cytotoxicity based on chemical structure and mode of action of new chemical entities. The clustering of the drug-induced effects on the hMDM bioenergetic parameters are reflected in the clustering of the effects of the drugs on cytokine production in hMDMs, demonstrating concurrence between the effects of the drugs on the metabolism and functioning of the macrophages. These findings can be used as a benchmark to establish XF analysis as a new tool to assay cytotoxicity in TB drug development.


title: Visualization of trigeminal ganglion sensory neuronal signaling regulated by Cdk5

abstract: The mechanisms underlying facial and oral pain are still incompletely understood, posing major therapeutic challenges. Cyclin-dependent kinase 5 (Cdk5) is a key neuronal kinase involved in pain signaling. However, the regulatory roles of Cdk5 in orofacial pain signaling and the possibility of therapeutic intervention at the level of mouse trigeminal ganglion primary neurons remain elusive. In this study, we used optimized intravital imaging to directly compare trigeminal neuronal activities after mechanical, thermal, and chemical stimulation. We then tested whether facial inflammatory pain in mice could be alleviated by the Cdk5 inhibitor peptide TFP5. We demonstrated regulation of total Ca2+ intensities by Cdk5 activity using transgenic and knockout mouse models. In mice with orofacial inflammation, application of TFP5 specifically decreased total Ca2+ intensities in response to noxious stimuli. It also alleviated inflammation-induced allodynia by inhibiting activation of trigeminal peripheral sensory neurons. Cdk5 inhibitors may provide promising non-opioid candidates for pain treatment.


title: Spatial-temporal targeted and non-targeted surveys to assess microbiological composition of drinking water in Puerto Rico following Hurricane Maria.

abstract:
Loss of basic utilities, such as drinking water and electricity distribution, were sustained for months in the aftermath of Hurricane Marias (HM) landfall in Puerto Rico (PR) in September 2017. The goal of this study was to assess if there was deterioration in biological quality of drinking water due to these disruptions. This study characterized the microbial composition of drinking water following HM across nine drinking water systems (DWSs) in PR and utilized an extended temporal sampling campaign to determine if changes in the drinking water microbiome were indicative of HM associated disturbance followed by recovery. In addition to monitoring water chemistry, the samples were subjected to culture independent targeted and non-targeted microbial analysis including quantitative PCR (qPCR) and genome-resolved metagenomics. The qPCR results showed that residual disinfectant was the major driver of bacterial concentrations in tap water with marked decrease in concentrations from early to late sampling timepoints. While Mycobacterium avium and Pseudomonas aeruginosa were not detected in any sampling locations and timepoints, genetic material from Leptospira and Legionella pneumophila were transiently detected in a few sampling locations. The majority of metagenome assembled genomes (MAGs) recovered from these samples were not associated with pathogens and were consistent with bacterial community members routinely detected in DWSs. Further, whole metagenome-level comparisons between drinking water samples collected in this study with samples from other full-scale DWS indicated no significant deviation from expected community membership of the drinking water microbiome. Overall, our results suggest that disruptions due to HM did not result in significant and sustained deterioration of biological quality of drinking water at our study sites.

Graphical Abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=86 SRC="FIGDIR/small/442998v1_ufig1.gif" ALT="Figure 1">
View larger version (29K):
org.highwire.dtl.DTLVardef@100c1e9org.highwire.dtl.DTLVardef@bae402org.highwire.dtl.DTLVardef@1cf41bborg.highwire.dtl.DTLVardef@1e46105_HPS_FORMAT_FIGEXP  M_FIG C_FIG


title: Changes in wheat rhizosphere microbiota in response to chemical inputs, plant genotype and phenotypic plasticity

abstract: Since modern wheat varieties are grown with chemical inputs, we ignore if changes observed in rhizosphere microorganisms between ancient and modern varieties are due to i) breeding-induced changes in plant genotype, ii) modifications of the environment via synthetic chemical inputs, or (iii) phenotypic plasticity, defined as the interaction between the genotype and the environment. In the field, we evaluated the effects of various wheat varieties (modern and ancient) grown with or without chemical inputs (N-fertilizer, fungicide and herbicide together) in a crossed factorial design. We analysed rhizosphere bacteria and fungi by amplicons sequencing and mycorrhizal association by microscopic observations. When considered independently of plant genotype, chemical inputs were responsible for an increase in dominance for bacteria and decrease in evenness for bacteria and fungi. Independently of inputs, modern varieties had richer and more even bacterial communities compared to ancient varieties. Phenotypic plasticity had a significant effect: bacterial and fungal diversity decreased when inputs were applied in ancient varieties but not in modern ones. Mycorrhiza were more abundant in modern than ancient varieties, and less abundant when using chemical inputs. Although neglected, phenotypic plasticity is important to understand the evolution of plant-microbiota associations and a relevant target in breeding programs.


title: Spatial point pattern analysis identifies mechanisms shaping the skin parasite landscape in Leishmania donovani infection

abstract: Increasing evidence suggests that infectiousness of hosts carrying parasites of the Leishmania donovani complex, the causative agents of visceral leishmaniasis, is linked to parasite repositories in the host skin. However, a detailed understanding of the dispersal and dispersion of these obligatory-intracellular parasites and their host phagocytes in the skin is lacking. Using endogenously fluorescent parasites as a proxy, we apply image analysis combined with spatial point pattern models borrowed from ecology to characterize dispersion of parasitized myeloid cells (including ManR+ and CD11c+ cells) and predict dispersal mechanisms in a previously described immunodeficient model of L. donovani infection. Our results suggest that after initial seeding of infection in the skin, heavily parasite-infected myeloid cells are found in patches that resemble innate granulomas. Spread of parasites from these initial patches subsequently occurs through infection of recruited myeloid cells, ultimately leading to self-propagating networks of patch clusters. This combination of imaging and ecological pattern analysis to identify mechanisms driving the skin parasite landscape offers new perspectives on myeloid cell behavior following parasitism by L. donovani and may also be applicable to elucidating the behavior of other intracellular tissue-resident pathogens and their host cells.


title: Alphavirus and flavivirus infectivity is disrupted by components of the bacterial cell wall

abstract:
The impact of the host microbiota on arbovirus infections is currently not well understood. Arboviruses are viruses transmitted through the bites of infected arthropods, predominantly mosquitoes or ticks. The first site of arbovirus inoculation is the biting site in the host skin, which is colonized by a complex microbial community that could possibly influence arbovirus infection. We demonstrated that pre-incubation of arboviruses with certain components of the bacterial cell wall, including lipopolysaccharides (LPS) of some Gram-negative bacteria and lipoteichoic acids or peptidoglycan of certain Gram-positive bacteria, significantly reduced arbovirus infectivity in vitro. This inhibitory effect was observed for arboviruses of different virus families, including chikungunya virus of the Alphavirus genus and Zika virus of the Flavivirus genus, showing that this is a broad phenomenon. A modest inhibitory effect was observed following incubation with a panel of heat-inactivated bacteria, including bacteria residing on the skin. No viral inhibition was observed after pre-incubation of cells with LPS. Furthermore, a virucidal effect of LPS on viral particles was noticed by electron microscopy. Therefore, the main inhibitory mechanism seems to be due to a direct effect on the virus particles. Together, these results suggest that bacteria are able to decrease the infectivity of alphaviruses and flaviviruses.

ImportanceDuring the past decades the world has experienced a vast increase in epidemics of alphavirus and flavivirus infections. These viruses can cause severe diseases such as hemorrhagic fever, encephalitis and arthritis. Several alpha- and flaviviruses, such as chikungunya virus, Zika virus and dengue virus, are significant global health threats because of their high disease burden, their widespread (re-)emergence and the lack of (good) anti-arboviral strategies. Despite the clear health burden, alphavirus and flavivirus infection and disease are not fully understood. A knowledge gap in the interplay between the host and the arbovirus is the potential interaction with host skin bacteria. Therefore, we studied the effect of (skin) bacteria and bacterial cell wall components on alphavirus and flavivirus infectivity in cell culture. Our results show that certain bacterial cell wall components markedly reduced viral infectivity by directly interacting with the virus particle.


title: Perturbation of alphavirus and flavivirus infectivity by components of the bacterial cell wall


title: Learning dynamics by computational integration of single cell genomic and lineage information

abstract: A goal of single cell genome-wide profiling is to reconstruct dynamic transitions during cell differentiation, disease onset, and drug response. Single cell assays have recently been integrated with lineage tracing, a set of methods that identify cells of common ancestry to establish bona fide dynamic relationships between cell states. These integrated methods have revealed unappreciated cell dynamics, but their analysis faces recurrent challenges arising from noisy, dispersed lineage data. Here, we develop coherent, sparse optimization (CoSpar) as a robust computational approach to infer cell dynamics from single-cell genomics integrated with lineage tracing. CoSpar is robust to severe down-sampling and dispersion of lineage data, which enables simpler, lower-cost experimental designs and requires less calibration. In datasets representing hematopoiesis, reprogramming, and directed differentiation, CoSpar identifies fate biases not previously detected, predicting transcription factors and receptors implicated in fate choice. Documentation and detailed examples for common experimental designs are available at https://cospar.readthedocs.io/.

category: developmental biology


title: Extracting socio-spatial networks from photo-ID data using multilevel multinomial models

abstract:
O_LIEstimating the impacts of anthropogenic disturbances requires an understanding of the habitat use patterns of individuals within a population. This is especially the case when disturbances are localized within a populations spatial range, as variation in habitat-use within a population can drastically alter the distribution of impacts.
C_LIO_LIHere, we illustrate the potential for multilevel multinomial models to generate spatial networks from capture-recapture data, a common data source use in wildlife studies to monitor population dynamics and habitat use. These spatial networks capture which regions of a populations spatial distribution share similar/dissimilar individual usage patterns, and can be especially useful for detecting structured habitat use within the populations spatial range.
C_LIO_LIUsing simulations and 18 years of capture-recapture data from St. Lawrence Estuary (SLE) beluga, we show that this approach can successfully estimate the magnitude of similarities/dissimilarities in individual usage patterns across sectors, and identify sectors that share similar individual usage patterns that differ from other sectors, i.e., structured habitat use. In the case of SLE beluga, this method identified multiple clusters of individuals, each preferentially using restricted areas within their summer range of the SLE.
C_LIO_LISynthesis and applications. Multilevel multinomial models can be effective at estimating spatial structure in habitat use within wildlife populations sampled by capture-recapture of individuals. Our finding of a structured habitat use within the SLE beluga summer range has direct implications for estimating individual exposures to localized stressors, such as underwater noise from shipping or other activities.
C_LI


title: Multilevel multinomial models as a tool to extract population spatial networks from capture-recapture data


title: Extracting spatial networks from capture-recapture data reveal structured habitat use within a marine mammal's spatial range