title: Idea paper: An envelope model of ecological disturbance

abstract: Disturbance is common in natural ecosystems, but increasingly defines them. While there are many descriptions for the dynamics of an ecosystems response to disturbance, there are few descriptions for the dynamics of the disturbance itself. I describe a novel application of a model based on the production of amplitude envelopes in acoustics and electronic music synthesis, with varying parameters Attack, Decay, Sustain, and Release (ADSR). I show that varying the parameters of the ADSR model is sufficient to produce and vary the qualitative disturbance regimes described by previous authors, and is capable of producing dynamics not previously considered. I tested the utility of the ADSR model by applying it to a logistic growth model. I found that manipulating the attack and release parameters of the ADSR model changes the population dynamics estimated by these models. This implies that responses to disturbance are determined not only by the resilience and resistance of the ecological system, but also the dynamics of the disturbance itself. My hope is that the ADSR model will prove useful to researchers in either describing disturbances in long-term ecological data, or in producing disturbances for simulations or experiments.


title: Brain regulatory program predates central nervous system evolution

abstract: Understanding if bilaterian centralized nervous systems (CNS) evolved once or multiple times has been debated for over a century. Recent efforts determined that the nerve chords found in bilaterian CNSs likely evolved independently, but the origin(s) of brains remains debatable. Developing brains are regionalized by stripes of gene expression along the anteroposterior axis. Gene homologs are expressed in the same relative order in disparate species, which has been interpreted as evidence for homology. However, regionalization programs resemble anteroposterior axial patterning programs, which supports an alternative model by which conserved expression in brains arose convergently through the independent co-option of deeply conserved axial patterning programs. To begin resolving these hypotheses, we sought to determine when the neurogenic role for axial programs evolved. Here we show that the nerve net in the cnidarian Nematostella vectensis and bilaterian brain are regionalized by the same molecular programs, which indicates nervous system regionalization predates the emergence of bilaterians and CNSs altogether. This argues that shared regionalization mechanisms are insufficient to support the homology of brains and supports the notion that axial programs were able to be co-opted multiple times during evolution of brains.

category: evolutionary biology


title: TIMSCONVERT: A workflow to convert trapped ion mobility data to open data formats

abstract:
SummaryAdvances in mass spectrometry instrumentation have led to the development of mass spectrometers with ion mobility separation (IMS) capabilities and dual source instrumentation, but the current software ecosystem lacks interoperability with downstream data analysis using open-source software/pipelines. Here, we present TIMSCONVERT, a data conversion workflow from timsTOF fleX MS raw data files to size conscious mzML and imzML formats with minimal preprocessing to allow for compatibility with downstream data analysis tools, which we showcase with several examples using data acquired across different experiments and acquisition modalities on the timsTOF fleX.

Availability and ImplementationTIMSCONVERT and its documentation can be found at https://github.com/gtluu/timsconvert and is available as a standalone command line interface, Nextflow workflow, and online in the Global Natural Products Social (GNPS) platform (https://proteomics2.ucsd.edu/ProteoSAFe/index.jsp?params={%22workflow%22%3A%20%22TIMSCONVERT%22}).

ContactMingxun Wang, miw023@ucsd.edu

Supplementary InformationSupplementary data are available at Bioinformatics online.

category: bioinformatics


title: Video-Based Deep Learning to Detect Dyssynergic Defecation with 3D High-Definition Anorectal Manometry

abstract:
BACKGROUND AND AIMSEvaluation for dyssynergia is the most common reason that gastroenterologists refer patients for anorectal manometry, because dyssynergia is amenable to biofeedback by physical therapists. High-definition anorectal manometry (3D-HDAM) is a promising technology to evaluate anorectal physiology, but adoption remains limited by its sheer complexity. We developed a 3D-HDAM deep learning algorithm to evaluate for dyssynergia.

METHODSSpatial-temporal data were extracted from consecutive 3D-HDAM studies performed between 2018-2020 at a tertiary institution. The technical procedure and gold standard definition of dyssynergia were based on the London consensus, adapted to the needs of 3D-HDAM technology. Three machine learning models were generated: (1) traditional machine learning informed by conventional anorectal function metrics, (2) deep learning, and (3) a hybrid approach. Diagnostic accuracy was evaluated using bootstrap sampling to calculate area-under-the-curve (AUC). To evaluate overfitting, models were validated by adding 502 simulated defecation maneuvers with diagnostic ambiguity.

RESULTS302 3D-HDAM studies representing 1,208 simulated defecation maneuvers were included (average age 55.2 years; 80.5% women). The deep learning model had comparable diagnostic accuracy (AUC=0.91 [95% confidence interval 0.89-0.93]) to traditional (AUC=0.93[0.92-0.95]) and hybrid (AUC=0.96[0.94-0.97]) predictive models in training cohorts. However, the deep learning model handled ambiguous tests more cautiously than other models; the deep learning model was more likely to designate an ambiguous test as inconclusive (odds ratio=4.21[2.78-6.38]) versus traditional/hybrid approaches.

CONCLUSIONSBy considering complex spatial-temporal information beyond conventional anorectal function metrics, deep learning on 3D-HDAM technology may enable gastroenterologists to reliably identify and manage dyssynergia in broader practice.


title: Comparative analysis of gene prediction tools for viral genome annotation

abstract: The number of newly available viral genomes and metagenomes has increased exponentially since the development of high throughput sequencing platforms and genome analysis tools. Bioinformatic annotation pipelines are largely based on open reading frame (ORF) calling software, which identifies genes independently of the sequence taxonomical background. Although ORF-calling programs provide a rapid genome annotation, they can misidentify ORFs and start codons; errors that might be perpetuated and propagated over time. This study evaluated the performance of multiple ORF-calling programs for viral genome annotation against the complete RefSeq viral database. Programs outputs varied when considering the viral nucleic acid type versus the viral host. According to the number of ORFs, Prodigal and Metaprodigal were the most accurate programs for DNA viruses, while FragGeneScan and Prodigal generated the most accurate outputs for RNA viruses. Similarly, Prodigal outperformed the benchmark for viruses infecting prokaryotes, and GLIMMER and GeneMarkS produced the most accurate annotations for viruses infecting eukaryotes. When the coordinates of the ORFs were considered, Prodigal scored high for all scenarios except for RNA viruses, where GeneMarkS generated the most reliable results. Overall, the quality of the coordinates predicted for RNA viruses was poorer than for DNA viruses, suggesting the need for improved ORF-calling programs to deal with RNA viruses. Moreover, none of the ORF-calling programs reached 90% accuracy for annotation of DNA viruses. Any automatic annotation can still be improved by manual curation, especially when the presence of ORFs is validated with wet-lab experiments. However, our evaluation of the current ORF-calling programs is expected to be useful for the improvement of viral genome annotation pipelines and highlights the need for more expression data to improve the rigor of reference genomes.


title: Buffering of transcription rate by mRNA half-life is a conserved feature of Rett syndrome models

abstract: Models of MECP2 dysfunction in Rett syndrome (RTT) assume that transcription rate changes directly correlate with altered steady-state mRNA levels. However, limited evidence suggests that transcription rate changes are buffered by poorly understood compensatory post-transcriptional mechanisms. Here we measure transcription rate and mRNA half-life changes in RTT patient neurons using RATE-seq, and reinterpret nuclear and whole-cell RNAseq from Mecp2 mice. Genes are dysregulated by changing transcription rate only or half-life only and are buffered when both are changed. We utilized classifier models to understand the direction of transcription rate changes based on gene-body DNA sequence, and combined frequencies of three dinucleotides were better predictors than contributions by CA and CG. MicroRNA and RNA-Binding Protein (RBP) motifs were enriched in 3UTRs of genes with half-life changes. Motifs for nuclear localized RBPs were enriched on buffered genes with increased transcription rate. Our findings identify post-transcriptional mechanisms in humans and mice that alter half-life only or buffer transcription rate changes when a transcriptional modulator gene is mutated in a neurodevelopmental disorder.

category: neuroscience


title: Two Stages of Speech Envelope Tracking in Human Auditory Cortex Modulated by Speech Intelligibility

abstract: The envelope is essential for speech perception. Recent studies have shown that cortical activity can track the acoustic envelope. However, whether the tracking strength reflects the extent of speech intelligibility processing remains controversial. Here, using stereo-electroencephalogram (sEEG) technology, we directly recorded the activity in human auditory cortex while subjects listened to either natural or noise-vocoded speech. These two stimuli have approximately identical envelopes, but the noise-vocoded speech does not have speech intelligibility. We found two stages of envelope tracking in auditory cortex: an early high-{gamma} (60-140 Hz) power stage (delay {approx} 49 ms) that preferred the noise-vocoded speech, and a late {theta} (4-8 Hz) phase stage (delay {approx} 178 ms) that preferred the natural speech. Furthermore, the decoding performance of high-{gamma} power was better in primary auditory cortex than in non-primary auditory cortex, consistent with its short tracking delay. We also found distinct lateralization effects: high-{gamma} power envelope tracking dominated left auditory cortex, while {theta} phase showed better decoding performance in right auditory cortex. In sum, we suggested a functional dissociation between high-{gamma} power and {theta} phase: the former reflects fast and automatic processing of brief acoustic features, while the latter correlates to slow build-up processing facilitated by speech intelligibility.


title: Peptide Aggregation Induced Immunogenic Rupture (PAIIR)

abstract: Under the influence of stress and membrane damage, cells undergo immunogenic cell death (ICD), which involves the release of damage associated molecular patterns (DAMPs), natural adjuvants for enhancing an immune response. In the presence of an antigen, released DAMPs can determine the type and magnitude of the immune response, and therefore the longevity and efficacy of an antigen-specific immunity. In the last decade, the immune response effect of ICD has been shown, yet there is no tool that can induce controlled ICD with predictable results, regardless of the cell type. We designed a peptide-based tool, called [II], for controlled damage to cell membrane to induce ICD and DAMPs release. Herein we describe a series of experiments that determine that the mechanism of action of [II] includes a caspase-dependent ICD and subsequent release of immune stimulating DAMPs, on various cell types. Moreover, we tested the hypothesis that controlled DAMP release via [II] in vivo was associated with enhancement of antigen-specific adaptive immunity with influenza hemagglutinin (HA) subunit vaccine. HA and [II] showed significantly higher HA specific IgG1 and IgG2a antibodies, compared to HA-only immunized mice, while the peptide itself did not elicit antibodies. In this paper, we demonstrate the first peptide-aggregation induced immunogenic rupture (PAIIR) approach as vaccine adjuvants for increasing both humoral and cellular immunity. In consideration of its ability to enhance IgG2a responses that are associated with heterosubtypic influenza virus protection, PAIIR is a promising adjuvant to promote universal protection upon influenza HA vaccination.

category: bioengineering


title: A qualitative mathematical model of immunocompetence with applications to SARS-CoV-2 immunity

abstract: A qualitative mathematical model of the notion of immunocompetence is developed, based on the formalism of Memory Evolutive Systems (MES), from which, immunocompetence is defined as an emergent structure of a higher order arising from the signal networks that are established between effector cells and molecules of the immune response in the presence of a given antigen. In addition, a possible mechanism of functorial nature is proposed, which may explain how immunocompetence is achieved in an organism endowed with innate and adaptive components of its immune system. Finally, a practical method to measure the immunocompetence status is established, using elements of the theory of small random graphs and taking into account the characteristics of the immune networks, established through transcriptional studies, of patients with severe COVID-19 and healthy patients, assuming that both types of patients were vaccinated with an effective biological against SARS-CoV-2.

category: systems biology


title: Cell Surface SARS-CoV-2 Nucleocapsid Protein Modulates Innate and Adaptive Immunity

abstract: SARS-CoV-2 nucleocapsid protein (N) induces strong antibody and T cell responses. Although considered to be localized in the cytosol, we readily detect N on the surface of live cells. N released by SARS-CoV-2 infected cells or N-expressing transfected cells binds to neighboring cells by electrostatic high-affinity binding to heparan sulfate and heparin, but not other sulfated glycosaminoglycans. N binds with high affinity to 11 human chemokines, including CXCL12{beta}, whose chemotaxis of leukocytes is inhibited by N from SARS-CoV-2, SARS-CoV-1, and MERS CoV. Anti-N Abs bound to the surface of N expressing cells activate Fc receptor-expressing cells. Our findings indicate that cell surface N manipulates innate immunity by sequestering chemokines and can be targeted by Fc expressing innate immune cells. This, in combination with its conserved antigenicity among human CoVs, advances its candidacy for vaccines that induce cross-reactive B and T cell immunity to SARS-CoV-2 variants and other human CoVs, including novel zoonotic strains.

category: microbiology


title: Present and future distribution of bat hosts of sarbecoviruses: implications for conservation and public health

abstract: Global changes in response to human encroachment into natural habitats and carbon emissions are driving the biodiversity extinction crisis and increasing disease emergence risk. Host distributions are one critical component to identify areas at risk of spillover, and bats act as reservoirs of diverse viruses. We developed a reproducible ecological niche modelling pipeline for bat hosts of SARS-like viruses (subgenus Sarbecovirus), given that since SARS-CoV-2 emergence several closely-related viruses have been discovered and sarbecovirus-host interactions have gained attention. We assess sampling biases and model bats current distributions based on climate and landscape relationships and project future scenarios. The most important predictors of species distribution were temperature seasonality and cave availability. We identified concentrated host hotspots in Myanmar and projected range contractions for most species by 2100. Our projections indicate hotspots will shift east in Southeast Asia in >2 {degrees}C hotter locations in a fossil-fueled development future. Hotspot shifts have implications for conservation and public health, as loss of population connectivity can lead to local extinctions, and remaining hotspots may concentrate near human populations.

category: ecology


title: Structural models of SARS-CoV-2 Omicron variant in complex with ACE2 receptor or antibodies suggest altered binding interfaces

abstract: There is enormous ongoing interest in characterizing the binding properties of the SARS-CoV-2 Omicron Variant of Concern (VOC) (B.1.1.529), which continues to spread towards potential dominance worldwide. To aid these studies, based on the wealth of available structural information about several SARS-CoV-2 variants in the Protein Data Bank (PDB) and a modeling pipeline we have previously developed for tracking the ongoing global evolution of SARS-CoV-2 proteins, we provide a set of computed structural models (henceforth models) of the Omicron VOC receptor-binding domain (omRBD) bound to its corresponding receptor Angiotensin-Converting Enzyme (ACE2) and a variety of therapeutic entities, including neutralizing and therapeutic antibodies targeting previously-detected viral strains. We generated bound omRBD models using both experimentally-determined structures in the PDB as well as machine learningbased structure predictions as starting points. Examination of ACE2-bound omRBD models reveals an interdigitated multi-residue interaction network formed by omRBD-specific substituted residues (R493, S496, Y501, R498) and ACE2 residues at the interface, which was not present in the original Wuhan-Hu-1 RBD-ACE2 complex. Emergence of this interaction network suggests optimization of a key region of the binding interface, and positive cooperativity among various sites of residue substitutions in omRBD mediating ACE2 binding. Examination of neutralizing antibody complexes for Barnes Class 1 and Class 2 antibodies modeled with omRBD highlights an overall loss of interfacial interactions (with gain of new interactions in rare cases) mediated by substituted residues. Many of these substitutions have previously been found to independently dampen or even ablate antibody binding, and perhaps mediate antibody-mediated neutralization escape (e.g., K417N). We observe little compensation of corresponding interaction loss at interfaces when potential escape substitutions occur in combination. A few selected antibodies (e.g., Barnes Class 3 S309), however, feature largely unaltered or modestly affected protein-protein interfaces. While we stress that only qualitative insights can be obtained directly from our models at this time, we anticipate that they can provide starting points for more detailed and quantitative computational characterization, and, if needed, redesign of monoclonal antibodies for targeting the Omicron VOC Spike protein. In the broader context, the computational pipeline we developed provides a framework for rapidly and efficiently generating retrospective and prospective models for other novel variants of SARS-CoV-2 bound to entities of virological and therapeutic interest, in the setting of a global pandemic.

category: biophysics


title: Structure of the 5' untranslated region in SARS-CoV-2 genome and its specific recognition by innate immune system via the human oligoadenylate synthase 1

abstract: 2-5-Oligoadenylate synthetase 1 (OAS1) is one of the key enzymes driving the innate immune system response to SARS-CoV-2 infection whose activity has been related to COVID-19 severity. In particular, OAS1 is a sensor of endogenous RNA that triggers the 2-5 oligoadenylate/RNase L pathway in response to viral infections, ultimately activating the RNA-Lyase which cleaves endogenous and exogenous RNA impeding the viral maturation. Upon SARS-CoV-2 infection, OAS1 is responsible for the recognition of viral RNA and has been shown to possess a particularly high sensitivity for the 5-untranslated (5-UTR) RNA region, which is organized in a double-strand stem loop motif (SLA). Since the structure of the nucleic acid/protein complex has not been resolved, here we report its structure obtained by molecular modeling, including enhanced sampling approaches. We also pinpoint that the SL1 region enhances the interaction network with the enzyme, promoting specific hydrogen bonds, absent in normal double strand RNA fragments, hence rationalizing the high affinity for OAS1.


title: Transcriptional and Cytotoxic Responses of Human Intestinal Organoids to Interferon Types I, II, and III

abstract: The three types of interferon (IFN) have roles in antimicrobial immunity and inflammation that must be properly balanced to maintain tissue homeostasis. For example, IFNs are elevated in the context of inflammatory bowel disease and may synergize with inflammatory cytokines such as tumor necrosis factor alpha (TNF) to promote tissue damage. Prior studies suggest that in mouse intestinal epithelial cells (IECs), type III IFNs are preferentially produced during viral infections and are less cytotoxic than type I IFN. Here, we generated human IEC organoid lines from biopsies of ileum, ascending colon, and sigmoid colon of three healthy subjects to establish the baseline responses of normal human IECs to types I, II, and III IFN. We found that all IFN types elicited responses that were qualitatively consistent across intestinal biopsy sites. However, IFN types differed in magnitude of STAT1 phosphorylation and identity of genes in their downstream transcriptional programs. Specifically, there was a core transcriptional module shared by IFN types, but types I and II IFN stimulated unique transcriptional modules beyond this core gene signature. The transcriptional modules of type I and II IFN included pro-apoptotic genes, and expression of these genes correlated with potentiation of TNF cytotoxicity. These data define the response profiles of healthy human IEC organoids across IFN types, and suggest that cytotoxic effects mediated by TNF in inflamed tissues may be amplified by a simultaneous high-magnitude IFN response.

category: immunology


title: Acetyltransferases GCN5 and PCAF are required for B cell maturation

abstract: B lymphocyte development has two DNA recombination processes: V(D)J recombination of the immunoglobulin (Igh) gene variable region and class switching of the Igh constant regions from IgM to IgG, IgA, or IgE. V(D)J recombination is required for successful maturation of B cells from pro-B to pre-B to immature-B and then to mature B cells in the bone marrow. CSR occurs outside of the bone marrow when mature B cells migrate to peripheral lymphoid organs, such as spleen and lymph nodes. Both V(D)J recombination and CSR depend on an open chromatin state that makes DNA accessible to specific enzymes, recombination activating gene (RAG), and activation-induced cytidine deaminase (AID). Acetyltransferases GCN5 and PCAF possess redundant functions acetylating histone H3 lysine 9 (H3K9). Here, we generated a mouse model that lacks both GCN5 and PCAF in B cells. We found that double-deficient mice possess low levels of mature B cells in the bone marrow and peripheral organs, an accumulation of pro-B cells in bone marrow, and reduced CSR levels. We conclude that both GCN5 and PCAF are required for B cell development in vivo.


title: Acetyltransferases GCN5 and PCAF are required for B lymphocyte maturation in mice


title: Cell transcriptomic atlas of the non-human primate Macaca fascicularis

abstract: Studying tissue composition and function in non-human primates (NHP) is crucial to understand the nature of our own species. Here, we present a large-scale single-cell and single-nucleus transcriptomic atlas encompassing over one million cells from 43 tissues from the adult NHP Macaca fascicularis. This dataset provides a vast, carefully annotated, resource to study a species phylogenetically close to humans. As proof of principle, we have reconstructed the cell-cell interaction networks driving Wnt signalling across the body, mapped the distribution of receptors and co-receptors for viruses causing human infectious diseases and intersected our data with human genetic disease orthologous coordinates to identify both expected and unexpected associations. Our Macaca fascicularis cell atlas constitutes an essential reference for future single-cell studies in human and NHP.

category: cell biology


title: Conformational variation in enzyme catalysis: A structural study on catalytic residues

abstract: Conformational variation in catalytic residues can be captured as alternative snapshots in enzyme crystal structures. Addressing the question of whether active site flexibility is an intrinsic and essential property of enzymes for catalysis, we present a comprehensive study on the 3D variation of active sites of 925 enzyme families, using explicit catalytic residue annotations from the Mechanism and Catalytic Site Atlas and structural data from the Protein Data Bank. Through weighted pairwise superposition of the functional atoms of active sites, we captured structural variability at single-residue level and examined the geometrical changes as ligands bind or as mutations occur. We demonstrate that catalytic centres of enzymes can be inherently rigid or flexible to various degrees according to the function they perform, and structural variability most often involves a subset of the catalytic residues, usually those not directly involved in the formation or cleavage of bonds. Moreover, data suggest that 2/3 of active sites are flexible, and in half of those, flexibility is only observed in the side chain. The goal of this work is to characterise our current knowledge of the extent of flexibility at the heart of catalysis and ultimately place our findings in the context of the evolution of catalysis as enzymes evolve new functions and bind different substrates.


title: High performance of a GPU-accelerated variant calling tool in genome data analysis

abstract: Rapid advances in next-generation sequencing (NGS) have facilitated ultralarge population and cohort studies that utilized whole-genome sequencing (WGS) to identify DNA variants that may impact gene function. Massive sequencing data require highly efficient bioinformatics tools to complete read alignment and variant calling as the fundamental analysis. Multiple software and hardware acceleration strategies have been developed to boost the analysis speed. This study comprehensively evaluated the germline variant calling of a GPU-based acceleration tool, BaseNumber, using WGS datasets from several sources, including gold-standard samples from the Genome in a Bottle (GIAB) project and the Golden Standard of China Genome (GSCG) project, resequenced GSCG samples, and 100 in-house samples from the China Deafness Genetics Consortium (CDGC) project. Sequencing data were analyzed on the GPU server using BaseNumber, the variant calling outputs of which were compared to the reference VCF or the results generated by the Burrows-Wheeler Aligner (BWA) + Genome Analysis Toolkit (GATK) pipeline on a generic CPU server. BaseNumber demonstrated high precision (99.32%) and recall (99.86%) rates in variant calls compared to the standard reference. The variant calling outputs of the BaseNumber and GATK pipelines were very similar, with a mean F1 of 99.69%. Additionally, BaseNumber took only 23 minutes on average to analyze a 48X WGS sample, which was 215.33 times shorter than the GATK workflow. The GPU-based BaseNumber provides a highly accurate and ultrafast variant calling capability, significantly improving the WGS analysis efficiency and facilitating time-sensitive tests, such as clinical WGS genetic diagnosis, and sheds light on the GPU-based acceleration of other omics data analyses.


title: Contrastive Cycle Adversarial Autoencoders for Single-cell Multi-omics Alignment and Integration

abstract: We have entered the multi-omics era, and we can measure cells from different aspects. When dealing with such multi-omics data, the first step is to determine the correspondence among different omics. In other words, we should match data from different spaces corresponding to the same object. This problem is particularly challenging in the single-cell multi-omics scenario because such data are very sparse with extremely high dimensions. Secondly, matched single-cell multi-omics data are rare and hard to collect. Furthermore, due to the limitations of the experimental environment, the data are usually highly noisy. To promote the single-cell multi-omics research, we overcome the above challenges, proposing a novel framework to align and integrate single-cell RNA-seq data and single-cell ATAC-seq data. Our approach can efficiently map the above data with high sparsity and noise from different spaces to a low-dimensional manifold in a unified space, making the downstream alignment and integration straightforward. Compared with the other state-of-the-art methods, our method performs better on both simulated and real single-cell data. On the real data, the performance improvement on accuracy over the previous methods is up to 55.7% regarding scRNA-seq and scATAC-seq data integration. Downstream trajectory inference analysis shows that our tool can transfer the labels from scRNA-seq to scATAC-seq with very high accuracy, which indicates our methods effectiveness.


title: Stabilization of Microbial Communities by Responsive Phenotypic Switching

abstract: Clonal microbes can switch between different phenotypes and recent theoretical work has shown that stochastic switching between these subpopulations can stabilize microbial communities. This phenotypic switching need not be stochastic, however, but can also be in response to environmental factors, both biotic and abiotic. Here, motivated by the bacterial persistence phenotype, we explore the ecological effects of such responsive switching by analyzing phenotypic switching in response to competing species. We show how the stability of microbial communities with responsive switching differs generically from that of communities with stochastic switching only. To understand this effect, we go on to analyse simple two-species models. Combining exact results and numerical simulations, we extend the classical stability results for models of two competing species without phenotypic variation to the case where one of the two species switches, stochastically and responsively, between two phenotypes. In particular, we show that responsive switching can stabilize coexistence even when stochastic switching on its own does not affect the stability of the community.


title: Landscape of epithelial-mesenchymal plasticity as an emergent property of coordinated teams in regulatory networks

abstract: Elucidating the principles of cellular decision-making is of fundamental importance. These decisions are often orchestrated by underlying regulatory networks. While we understand the dynamics of simple network motifs, how do large networks lead to a limited number of phenotypes, despite their complexity, remains largely elusive. Here, we investigate five different networks governing epithelial-mesenchymal plasticity and identified a latent design principles in their topology that limits their phenotypic repertoire - the presence of two "teams" of nodes engaging in a mutually inhibitory feedback loop, forming a toggle switch. These teams are specific to these networks and directly shape the phenotypic landscape and consequently the frequency and stability of terminal phenotypes vs. the intermediary ones. Our analysis reveals that network topology alone can contain information about phenotypic distributions it can lead to, thus obviating the need to simulate them. We unravel topological signatures that can drive canalization of cell-fates during diverse decision-making processes.


title: Continual reelin signaling by the prime neurogenic niche of the adult brain

abstract: During development reelin sets the pace of neocortical neurogenesis enabling in turn newborn neurons to migrate, but whether and, if so, how reelin signaling affects the adult neurogenic niches remains uncertain. We show that reelin signaling, resulting in Dab1 phosphorylation, occurs in the ependymal-subependymal zone (EZ/SEZ) of the lateral ventricles where, along with its associated rostral migratory stream (RMS), the highest density of functional ApoER2 accumulates. Mice deficient for reelin, ApoER2 or Dab1 exhibit enlarged ventricles and dysplastic RMS. Moreover, while the conditional ablation of Dab1 in neural progenitor cells (NPCs) enlarges the ventricles and impairs neuroblasts clearance from the SEZ, the transgenic misexpression of reelin in NPCs of reelin-deficient mice normalizes the ventricular lumen and the density of ependymal cilia, ameliorating in turn neuroblasts migration; consistently, intraventricular infusion of reelin reroutes neuroblasts. These results demonstrate that reelin signaling persists sustaining the germinal niche of the lateral ventricles and influencing neuroblasts migration in the adult brain.


title: FLAMs: A self-replicating ex vivo model of alveolar macrophages for functional genetic studies

abstract: Alveolar macrophages (AMs) are tissue resident cells in the lungs derived from the fetal liver that maintain lung homeostasis and respond to inhaled stimuli. While the importance of AMs is undisputed, they remain refractory to standard experimental approaches and high-throughput functional genetics as they are challenging to isolate and rapidly lose AM properties in standard culture. This limitation hinders our understanding of key regulatory mechanisms that control AM maintenance and function. Here, we describe the development of a new model, fetal liver-derived alveolar-like macrophages (FLAMs), which maintains cellular morphologies, expression profiles, and functional mechanisms similar to murine AMs. FLAMs combine treatment with two key cytokines for AM maintenance, GM-CSF and TGF{beta}. We leveraged the long-term stability of FLAMs to develop functional genetic tools using CRISPR-Cas9-mediated gene editing. Targeted editing confirmed the role of AM-specific gene Marco and the IL-1 receptor Il1r1 in modulating the AM response to crystalline silica. Furthermore, a genome-wide knockout library using FLAMs identified novel genes required for surface expression of the AM marker Siglec-F, most notably those related to the peroxisome. Taken together, our results suggest that FLAMs are a stable, self-replicating model of AM function that enables previously impossible global genetic approaches to define the underlying mechanisms of AM maintenance and function.


title: The Omicron variant is highly resistant against antibody-mediated neutralization - implications for control of the COVID-19 pandemic

abstract: The rapid spread of the SARS-CoV-2 Omicron variant suggests that the virus might become globally dominant. Further, the high number of mutations in the viral spike-protein raised concerns that the virus might evade antibodies induced by infection or vaccination. Here, we report that the Omicron spike was resistant against most therapeutic antibodies but remained susceptible to inhibition by Sotrovimab. Similarly, the Omicron spike evaded neutralization by antibodies from convalescent or BNT162b2-vaccinated individuals with 10- to 44-fold higher efficiency than the spike of the Delta variant. Neutralization of the Omicron spike by antibodies induced upon heterologous ChAdOx1/BNT162b2-vaccination or vaccination with three doses of BNT162b2 was more efficient, but the Omicron spike still evaded neutralization more efficiently than the Delta spike. These findings indicate that most therapeutic antibodies will be ineffective against the Omicron variant and that double immunization with BNT162b2 might not adequately protect against severe disease induced by this variant.


title: Quantification of nuclear transport inhibition by SARS-CoV-2 ORF6 using a broadly applicable live-cell dose-response pipeline

abstract: SARS coronavirus ORF6 inhibits the classical nuclear import pathway to antagonize host antiviral responses. Several models were proposed to explain its inhibitory function, but quantitative measurement is needed for model evaluation and refinement. We report a broadly applicable live-cell method for calibrated dose-response characterization of the nuclear transport alteration by a protein of interest. Using this method, we found that SARS-CoV-2 ORF6 is [~]5 times more potent than SARS-CoV-1 ORF6 in inhibiting bidirectional nuclear transport, due to differences in the NUP98-binding C-terminal region that is required for the inhibition. The N-terminal region was also required, but its membrane binding function was dispensable, since loss of the inhibitory function due to N-terminal truncation was rescued by forced oligomerization using a soluble construct. Based on these data, we propose that the hydrophobic N-terminal region drives oligomerization of ORF6 to multivalently cross-link the FG domains of NUP98 at the nuclear pore complex.


title: PTRN-1/CAMSAP and NOCA-2/NINEIN are required for microtubule polarity in Caenorhabditis elegans dendrites

abstract: The neuronal microtubule cytoskeleton is key to establish axon-dendrite polarity. Dendrites are characterized by the presence of minus-end out microtubules, however the mechanisms that organize these microtubules minus-end out is still poorly understood. Here, we characterized the role of two microtubule minus-end related proteins in this process in Caenorhabditis elegans, the microtubule minus-end stabilizing protein CAMSAP (PTRN-1) and a NINEIN homologue (NOCA-2). We found that CAMSAP and NINEIN function in parallel to mediate microtubule organization in dendrites. During dendrite outgrowth, RAB-11 positive vesicles localized to the dendrite tip function as a microtubule organizing center (MTOC) to nucleate microtubules. In the absence of either CAMSAP or NINEIN, we observed a low penetrance MTOC vesicles mis-localization to the cell body, and a nearly fully penetrant phenotype in double mutant animals. This suggests that both proteins are important for localizing the MTOC vesicles to the growing dendrite tip to organize microtubules minus-end out. Whereas NINEIN localizes to the MTOC vesicles where it is important for the recruitment of the microtubule nucleator {gamma}-tubulin, CAMSAP localizes around the MTOC vesicles and is co-translocated forward with the MTOC vesicles upon dendritic growth. Together, these results indicate that microtubule nucleation from the MTOC vesicles and microtubule stabilization are both important to localize the MTOC vesicles distally to organize dendritic microtubules minus-end out.


title: Antimicrobial resistance determinants in silage

abstract: Animal products may play a role in developing and spreading antimicrobial resistance in several ways. On the one hand, residues of antibiotics not adequately used in animal farming can enter the human body via food. But resistant bacteria may also be present in animal products, which can transfer the antimicrobial resistance genes (ARG) to the bacteria in the consumers body by horizontal gene transfer. As previous studies have shown that fermented foods have a meaningful ARG content, it is indicated that such genes may also be present in silage used as mass feed in the cattle sector. In our study, we aspired to answer what ARGs occur in silage and what mobility characteristics they have? For this purpose, we have analyzed bioinformatically 52 freely available deep sequenced silage samples from shotgun metagenome next-generation sequencing. A total of 17 perfect matched ARGs occurred 55 times in the samples. More than half of these ARGs are mobile because they can be linked to integrative mobile genetic elements, prophages or plasmids. Our results point to a neglected but substantial ARG source in the food chain.

category: genomics


title: Mutations in RBD of SARS-CoV-2 Omicron variant result stronger binding to human ACE2 protein

abstract: The COVID-19 pandemic caused by the SARS-CoV-2 virus has led to more than 270 million infections and 5.3 million of deaths worldwide. Several major variants of SARS-CoV-2 have emerged and posed challenges in controlling the pandemic. The recently occurred Omicron variant raised serious concerns about reducing the efficacy of vaccines and neutralization antibodies due to its vast mutations. We have modelled the complex structure of the human ACE2 protein and the receptor binding domain (RBD) of Omicron Spike protein (S-protein), and conducted atomistic molecular dynamics simulations to study the binding interactions. The analysis shows that the Omicron RBD binds more strongly to the human ACE2 protein than the original strain. The mutations at the ACE2-RBD interface enhance the tight binding by increasing hydrogen bonding interaction and enlarging buried solvent accessible surface area.


title: Mutations on RBD of SARS-CoV-2 Omicron variant result in stronger binding to human ACE2 receptor


title: SARS-CoV2 variant-specific replicating RNA vaccines protect from disease and pathology and reduce viral shedding following challenge with heterologous SARS-CoV2 variants of concern

abstract: Despite mass public health efforts, the SARS-CoV2 pandemic continues as of late-2021 with resurgent case numbers in many parts of the world. The emergence of SARS-CoV2 variants of concern (VoC) and evidence that existing vaccines that were designed to protect from the original strains of SARS-CoV-2 may have reduced potency for protection from infection against these VoC is driving continued development of second generation vaccines that can protect against multiple VoC. In this report, we evaluated an alphavirus-based replicating RNA vaccine expressing Spike proteins from the original SARS-CoV-2 Alpha strain and recent VoCs delivered in vivo via a lipid inorganic nanoparticle. Vaccination of both mice and Syrian Golden hamsters showed that vaccination induced potent neutralizing titers against each homologous VoC but reduced neutralization against heterologous challenges. Vaccinated hamsters challenged with homologous SARS-CoV2 variants exhibited complete protection from infection. In addition, vaccinated hamsters challenged with heterologous SARS-CoV-2 variants exhibited significantly reduced shedding of infectious virus. Our data demonstrate that this vaccine platform elicits significant protective immunity against SARS-CoV2 variants and supports continued development of this platform.


title: SARS-CoV-2 Spike triggers barrier dysfunction and vascular leak via integrins and TGF-β signaling

abstract: Severe COVID-19 is associated with epithelial and endothelial barrier dysfunction within the lung as well as in distal organs. While it is appreciated that an exaggerated inflammatory response is associated with barrier dysfunction, the triggers of this pathology are unclear. Here, we report that cell-intrinsic interactions between the Spike (S) glycoprotein of SARS-CoV-2 and epithelial/endothelial cells are sufficient to trigger barrier dysfunction in vitro and vascular leak in vivo, independently of viral replication and the ACE2 receptor. We identify an S-triggered transcriptional response associated with extracellular matrix reorganization and TGF-{beta} signaling. Using genetic knockouts and specific inhibitors, we demonstrate that glycosaminoglycans, integrins, and the TGF-{beta} signaling axis are required for S-mediated barrier dysfunction. Our findings suggest that S interactions with barrier cells are a contributing factor to COVID-19 disease severity and offer mechanistic insight into SARS-CoV-2 triggered vascular leak, providing a starting point for development of therapies targeting COVID-19 pathogenesis.


title: Improved neutralization of the SARS-CoV-2 Omicron variant after Pfizer-BioNTech BNT162b2 COVID-19 vaccine boosting

abstract: In late November 2021, the World Health Organization declared the SARS-CoV-2 lineage B.1.1.529 the fifth variant of concern, Omicron. This variant has acquired 15 mutations in the receptor binding domain of the spike protein, raising concerns that Omicron could evade naturally acquired and vaccine-derived immunity. We utilized an authentic virus, multicycle neutralisation assay to demonstrate that sera collected 1, 3, and 6 months post-two doses of Pfizer-BioNTech BNT162b2 has a limited ability to neutralise SARS-CoV-2. However, four weeks after a third dose, neutralizing antibody titres are boosted. Despite this increase, neutralising antibody titres are reduced 4-fold for Omicron compared to lineage A.2.2 SARS-CoV-2.


title: Following the track: accuracy and reproducibility of predation assessment on artificial caterpillars

abstract: Experimental studies of biotic interactions in real field conditions are essential to understand the structure and functioning of ecological networks. The use of artificial caterpillars to mimic actual prey availability is generally seen as a standard approach to compare the activity and diversity of predators along environmental gradients. Yet, even with standardized material and procedures, biases may still affect data from multiple observers with different expertise. We used pictures of artificial caterpillars with or without various predation marks, in an online survey that was targeted for the participants of the project, to evaluate the reliability of predation marks identification made by non-scientists and by scientists with and without previous experience in predation mark identification. Pictures of artificial caterpillars displayed typical marks left by birds, mammals and arthropods, as well as non-predation marks ( false positive). 357 respondents scanned 7140 pictures of these pictures. Self-declared scientists were more confident and accurate in their observations than non-scientists, but the differences in correct identifications among scientists and non-scientists were low. Self-declared scientists with previous experience were also more accurate than scientists without previous experience, while there were no differences in self-confidence among scientists with and without previous experience. Accuracy in predation mark identification did not differ among types of predators, but respondents were more keen to identify marks left by birds or mammals than arthropods. Our results have practical implications for the design of multi-observer projects relying on artificial caterpillars as a proxy to assess predation intensity, in particular in the context of citizen science.


title: Epistatic models predict mutable sites in SARS-CoV-2 proteins and epitopes

abstract:
The emergence of new variants of SARS-CoV-2 is a major concern given their potential impact on the transmissibility and pathogenicity of the virus as well as the efficacy of therapeutic interventions. Here, we predict the mutability of all positions in SARS-CoV-2 protein domains to forecast the appearance of unseen variants. Using sequence data from other coronaviruses, pre-existing to SARS-CoV-2, we build statistical models that do not only capture amino-acid conservation but more complex patterns resulting from epistasis. We show that these models are notably superior to conservation profiles in estimating the already observable SARS-CoV-2 variability. In the receptor binding domain of the spike protein, we observe that the predicted mutability correlates well with experimental measures of protein stability and that both are reliable mutability predictors (ROC AUC ~0.8). Most interestingly, we observe an increasing agreement between our model and the observed variability as more data become available over time, proving the anticipatory capacity of our model. When combined with data concerning the immune response, our approach identifies positions where current variants of concern are highly overrepresented. These results could assist studies on viral evolution, future viral outbreaks and, in particular, guide the exploration and anticipation of potentially harmful future SARS-CoV-2 variants.

Significance statementDuring the COVID pandemic, new SARS-CoV-2 variants emerge and spread, some being of major concern due to their increased infectivity or their capacity to reduce vaccine efficiency. Anticipating new mutations, which might give rise to new variants, would be of great interest. Here we construct sequence models predicting how mutable SARS-CoV-2 positions are, using a single SARS-CoV-2 sequence and databases of other coronaviruses. Predictions are tested against available mutagenesis data and the observed variability of SARS-CoV-2 proteins. Interestingly, our predictions agree increasingly with observations, as more SARS-CoV-2 sequences become available. Combining predictions with immunological data, we find a clear overrepresentation of mutations in current variants of concern. The approach may become relevant for potential outbreaks of future viral diseases.


title: No evidence of paralogous loci or new bona fide microRNAs in telomere to telomere (T2T) genomic data

abstract: The telomere to telomere (T2T) genome project discovered and mapped [~]240 million additional base pairs of primarily telomeric and centromeric reads. Much of this sequence was comprised of satellite sequences and large segmental duplications. We evaluated the extent to which human bona fide microRNAs (miRNAs) may be found in additional paralogous genomic loci or if previously undescribed microRNAs are present in these newly sequenced regions of the human genome. New genomic regions of the T2T project spanning [~]240 million bp of sequence were obtained and evaluated by blastn for the human miRNAs contained in MirGeneDB2.0 (N=556) and miRBase (N = 1917) along with all species of MirGeneDB2.0 miRNAs (N=10,899). Additionally, bowtie was used to compare unmapped reads from >4,000 primary cell samples to the new T2T sequence. Based on sequence and structure, no bona fide miRNAs were identified. Ninety-seven miRNAs of questionable authenticity (frequently known repeat elements) were identified from the miRBase dataset across the newly described regions of the human genome. These 97 represent only 51 miRNA families due to paralogy of highly similar miRNAs such as 24 members of the hsa-mir-548 family. Altogether, this data strongly supports our having identified widely expressed bona fide miRNAs in the human genome and move us further toward the completion of human miRNA discovery.


title: Personalized models of breast cancer desmoplasia reveal biomechanical determinants of drug penetration

abstract: Breast cancer desmoplasia heterogeneity contributes to high disease mortality due to discrepancies in treatment efficacy between patients. Personalized in vitro breast cancer models can be used for high throughput testing and ranking of therapeutic strategies to normalize the aberrant microenvironment in a patient-specific manner. Here, tumoroids assembled from patient-derived cells cultured in microphysiological systems including perfusable microvasculature reproduce key aspects of stromal and vascular dysfunction. Increased hyaluronic acid and collagen deposition, loss of vascular glycocalyx and reduced perfusion, and elevated interstitial fluid pressure in the models result in impaired drug distribution to tumor cells. We demonstrate the application of these personalized models as tools to rank molecular therapies for the normalization of the tumoroid microenvironment and to discover new therapeutic targets such as IL8 and CD44, which may ultimately improve drug efficacy in breast cancer patients.


title: EraSOR: Erase Sample Overlap in polygenic score analyses

abstract:
BackgroundPolygenic risk score (PRS) analyses are now routinely applied in biomedical research, with great hope that they will aid in our understanding of disease aetiology and contribute to personalized medicine. The continued growth of multi-cohort genome-wide association studies (GWASs) and large-scale biobank projects has provided researchers with a wealth of GWAS summary statistics and individual-level data suitable for performing PRS analyses. However, as the size of these studies increase, the risk of inter-cohort sample overlap and close relatedness increases. Ideally sample overlap would be identified and removed directly, but this is typically not possible due to privacy laws or consent agreements. This sample overlap, whether known or not, is a major problem in PRS analyses because it can lead to inflation of type 1 error and, thus, erroneous conclusions in published work.

ResultsHere, for the first time, we report the scale of the sample overlap problem for PRS analyses by generating known sample overlap across sub-samples of the UK Biobank data, which we then use to produce GWAS and target data to mimic the effects of inter-cohort sample overlap. We demonstrate that inter-cohort overlap results in a significant and often substantial inflation in the observed PRS-trait association, coefficient of determination (R2) and false-positive rate. This inflation can be high even when the absolute number of overlapping individuals is small if this makes up a notable fraction of the target sample. We develop and introduce EraSOR (Erase Sample Overlap and Relatedness), a software for adjusting inflation in PRS prediction and association statistics in the presence of sample overlap or close relatedness between the GWAS and target samples. A key component of the EraSOR approach is inference of the degree of sample overlap from the intercept of a bivariate LD score regression applied to the GWAS and target data, making it powered in settings where both have sample sizes over 1,000 individuals. Through extensive benchmarking using UK Biobank and HapGen2 simulated genotype-phenotype data, we demonstrate that PRSs calculated using EraSOR-adjusted GWAS summary statistics are robust to inter-cohort overlap in a wide range of realistic scenarios and are even robust to high levels of residual genetic and environmental stratification.

ConclusionThe results of all PRS analyses for which sample overlap cannot be definitively ruled out should be considered with caution given high type 1 error observed in the presence of even low overlap between base and target cohorts. Given the strong performance of EraSOR in eliminating inflation caused by sample overlap in PRS studies with large (>5k) target samples, we recommend that EraSOR be used in all future such PRS studies to mitigate the potential effects of inter-cohort overlap and close relatedness.


title: Black chromatin is indispensable for accurate simulations of Drosophila melanogaster chromatin structure.

abstract:
The interphase chromatin structure is extremely complex, precise and dynamic. Experimental methods can only show the frequency of interaction of the various parts of the chromatin. Therefore, it is extremely important to develop theoretical methods to predict the chromatin structure. In this publication, we describe the necessary factors for the effective modeling of the chromatin structure in Drosophila melanogaster. We also compared Monte Carlo with Molecular Dynamic methods. We showed that incorporating black, non-reactive chromatin is necessary for successfully prediction of chromatin structure, while the loop extrusion model or using Hi-C data as input are not essential for the basic structure reconstruction.

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=68 SRC="FIGDIR/small/472204v1_ufig1.gif" ALT="Figure 1">
View larger version (10K):
org.highwire.dtl.DTLVardef@65d59aorg.highwire.dtl.DTLVardef@1aa9cbforg.highwire.dtl.DTLVardef@18f5cborg.highwire.dtl.DTLVardef@b1034d_HPS_FORMAT_FIGEXP  M_FIG C_FIG


title: Acoustic camera system for measuring ultrasound communication in mice

abstract: To investigate molecular, genetic, and neural mechanisms underlying social behaviors and their deficits in neuropsychiatric disorders, social communication via ultrasonic vocalizations (USVs) in mice has received considerable attention as a powerful experimental model. The advances in sound localization technology have facilitated the analysis of vocal interactions between multiple mice. However, existing sound localization systems are built around distributed-microphone arrays, which require a special recording arena and long processing time. Here we report a novel acoustic camera system, USVCAM, which enables simpler and faster USV localization and assignment. The system comprises recently developed USV segmentation algorithms with a modification for overlapping vocalizations that results in high accuracy. Using USVCAM, we analyzed USV communications in a conventional home cage, and demonstrated novel vocal interactions in female ICR mice under a resident-intruder paradigm. The extended applicability and usability of USVCAM may facilitate future studies investigating normal and abnormal vocal communication and social behaviors, as well as the underlying physiological mechanisms.


title: Localized, time-dependent responses of rat cranial bone to repeated mild traumatic brain injuries

abstract: While it is well-established that bone responds dynamically to mechanical loading, the effects of mild traumatic brain injury (mTBI) on cranial bone composition are unclear. We hypothesized that repeated mTBI (rmTBI) would change the microstructure of cranial bones, without gross skull fractures. To address this, young adult female Piebald Viral Glaxo rats received sham, 1x, 2x or 3x closed-head mTBIs delivered at 24h intervals, using a weight-drop device custom built for reproducible impact. Skull bones were collected at 2 or 10 weeks after the final injury/sham procedure, imaged by micro computed tomography and analyzed at predetermined regions of interest. In the interparietal bone, proximal to the injury site, modest increases in bone thickness was observed at 2 weeks, particularly following 3x mTBI. By 10 weeks, 2x mTBI induced a robust increase in the volume and thickness of the interparietal bone, alongside a corresponding decrease in the volume of marrow cavities in the diploe region. In contrast, neither parietal nor frontal skull samples were affected by rmTBI. Our findings demonstrate time- and location-dependent effects of rmTBI on cranial bone structure, highlighting a need to consider microstructural alterations to cranial bone when assessing the consequences of rmTBI.


title: Adaptive discrimination of antigen risk by predictive coding in immune system

abstract: The immune system discriminates between harmful and harmless antigens based on past experiences; however, the underlying mechanism is largely unknown. From the viewpoint of machine learning, the learning system predicts the observation and updates the prediction based on prediction error, a process known as  predictive coding. Here, we modeled the population dynamics of T cells by adopting the concept of predictive coding; helper and regulatory T cells predict the antigen amount and excessive immune response, respectively. Their prediction error signals, possibly via cytokines, induce their differentiation to memory T cells. Through numerical simulations, we found that the immune system identifies antigen risks depending on the concentration and input rapidness of the antigen. Further, our model reproduced history-dependent discrimination, as in allergy onset and subsequent therapy. Together, this study provided a novel framework to improve our understanding of how the immune system adaptively learns the risks of diverse antigens.


title: Dehalobium species implicated in 2,3,7,8-tetrachloro-p-dioxin dechlorination in the contaminated sediments of Sydney Harbour Estuary

abstract: Polychlorinated dibenzo-p-dioxins and furans (PCDD/F) are some of the most environmentally recalcitrant and toxic compounds. They are naturally occurring and by-products of anthropogenic activity. Sydney Harbour Estuary (Sydney, Australia), is heavily contaminated with PCDD/F. Analysis of sediment cores revealed that the contamination source in Homebush Bay continues to have one of the highest levels of PCDD/F contamination in the world (5207 pg WHO-TEQ g-1) with >50% of the toxicity attributed to 2,3,7,8-tetrachlorodibenzo-p-dioxin (2,3,7,8-TCDD) the most toxic and concerning of the PCDD/F congeners. Comparison of congener profiles at the contamination source with surrounding bays and historical data provided evidence for the attenuation of 2,3,7,8-TCDD and other congeners at the source. This finding was supported by the detection of di-, mono- and unchlorinated dibenzo-p-dioxin. Microbial community analysis of sediments by 16S amplicon sequencing revealed an abundance of lineages from the class Dehalococcoidia (up to 15% of the community), including the genus Dehalobium (up to 0.5%). Anaerobic seawater enrichment cultures using perchloroethene as a more amenable growth substrate enriched only the Dehalobium population by more than six-fold. The enrichment culture then proved capable of reductively dechlorinating 2,3,7,8-TCDD to 2,3,7-TCDD and octachlorodibenzo-p-dibenzodioxin to hepta and hexa congeners. This work is the first to show microbial reductive dehalogenation of 2,3,7,8-TCDD with a bacterium from outside the Dehalococcoides genus, and one of only a few that demonstrates PCDD/F degradation in a marine environment.


title: Acyl-CoA Thioesterases; a rheostat that controls activated fatty acids modulates dengue virus serotype 2 replication

abstract: During infection with dengue viruses (DENVs), the lipid landscape within host cells is significantly altered to assemble membrane platforms that support viral replication and particle assembly. Fatty acyl-CoAs are key intermediates in the biosynthesis of complex lipids that form these membranes. They also function as key signaling lipids in the cell. Here, we carried out loss of function studies on acyl-CoA thioesterases (ACOTs), a family of enzymes that hydrolyze fatty acyl-CoAs to free fatty acids and coenzyme A, to understand their influence on the lifecycle of DENVs. Loss of function of the type I ACOTs 1 (cytoplasmic) and 2 (mitochondrial) together significantly increased DENV serotype 2 (DENV2) viral replication and infectious particle release. However, isolated knockdown of mitochondrial ACOT2 significantly decreased DENV2 protein translation, genome replication, and infectious virus release. Furthermore, loss of ACOT7 function, a mitochondrial type II ACOT, similarly suppressed DENV2. As ACOT1 and ACOT2 are splice variants, these data suggest that location (cytosol and mitochondria, respectively) rather than function of these proteins may account for the differences in DENV2 infection phenotype. Additionally, loss of mitochondrial ACOT2 and ACOT7 expression also altered the expression of several ACOTs located in multiple organelle compartments within the cell highlighting a complex relationship between ACOTs in the DENV2 virus lifecycle.


title: Microbial associations with Microscopic Colitis

abstract:
BACKGROUND AND OBJECTIVEMicroscopic colitis is a relatively common cause of chronic diarrhea and may be linked to luminal factors. Given the essential role of the microbiome in human gut health, analysis of microbiome changes associated with microscopic colitis could provide insights into the development of the disease.

METHODSWe enrolled patients who underwent colonoscopy for diarrhea. An experienced pathologist classified patients as having microscopic colitis (n=52) or controls (n=153). Research biopsies were taken from the ascending and descending colon, and the microbiome was characterized with Illumina sequencing. We analyzed the associations between microscopic colitis and microbiome with a series of increasingly complex models adjusted for a range of demographic and health factors.

RESULTSWe found that alpha-diversity was significantly lower in microscopic colitis cases compared to controls in the descending colon microbiome. In the descending colon, a series of models that adjusted for an increasing number of co-variates found taxa significantly associated with microscopic colitis, including Proteobacteria that was enriched in cases and Collinsella enriched in controls. While the alpha-diversity and taxa were not significantly associated with microscopic colitis in the ascending colon microbiome, the inference p-values based on ascending and descending microbiomes were highly correlated.

CONCLUSIONOur study demonstrates an altered microbiome in microscopic colitis cases compared to controls. Because both the cases and controls had diarrhea, we have identified candidate taxa that could be mechanistically responsible for the development of microscopic colitis independent of changes to the microbial community caused by diarrhea.

Significance of this studyO_ST_ABSWhat is already known about this subject?C_ST_ABSO_LIMicroscopic colitis is a common cause of chronic diarrhea. The exact etiology of microscopic colitis is unknown, but the gut microbiome is considered to play an important role.
C_LIO_LIPrior studies have reported microbial changes associated with microscopic colitis, but the results are not consistent. Prior studies have generally involved small numbers of patients, fecal samples, and comparison to healthy controls.
C_LI

What are the new findings?O_LIThe microbiome was altered in microscopic colitis patients.
C_LIO_LIMicroscopic colitis was associated with lower alpha-diversity, increase of inflammation related taxa and decrease of Collinsella in the descending colon microbiome independent of diarrhea status.
C_LIO_LIThe changes of microbial features were consistent between ascending and descending colon microbiomes, but with smaller changes in the ascending colon microbiome.
C_LI

How might it impact on clinical practice in the foreseeable future?O_LIThe altered microbial communities in patients with microscopic colitis suggest that the patients might benefit from prevention or treatment with live biotherapeutic products.
C_LI


title: Stressed out underground?: Transcriptomics illuminates the genome-wide response to heat stress in surface and subterranean diving beetles

abstract: Subterranean habitats are environmentally stable with respect to temperature, humidity, and the absence of light. The transition to a subterranean lifestyle might therefore be expected to cause considerable shifts in an organisms physiology; here, we investigate how subterranean colonisation affects thermal tolerance. Subterranean organisms might be at an increased risk of decline in the face of global temperature rises, but robust data on the fauna is lacking, particularly at the molecular level. In this study we compare the heat shock response of two species of diving beetle in the genus Paroster: one surface-dwelling (P. nigroadumbratus), the other restricted to a single aquifer (P. macrosturtensis). P. macrosturtensis has been previously established as having a lower thermal tolerance compared to surface-dwelling relatives, but the genomic basis of this difference is unknown. By sequencing transcriptomes of experimentally heat-shocked individuals and performing differential expression analysis, we demonstrate both species can mount a heat shock response at high temperatures (35C), in agreement with past survival experiments. However, the genes involved in these responses differ between species, and far greater genes are differentially expressed in the surface species, which may explain its more robust response to heat stress. In contrast, the subterranean species significantly upregulated the heat shock protein gene Hsp68 in the experimental setup under conditions it would likely encounter in nature (25C), suggesting it may be more sensitive to ambient stressors, e.g. handling. The results presented here contribute to an emerging narrative concerning weakened thermal tolerances in obligate subterranean organisms at the molecular level.


title: genomicSimulation: fast R functions for stochastic simulation of breeding programs

abstract: Simulation tools are key to designing and optimising breeding programs that are multi-year, high-effort endeavours. Tools that operate on real genotypes and integrate easily with other analysis software are needed to guide users to crossing decisions that best balance genetic gains and diversity to maintain gains in the future. This paper presents genomicSimulation, a fast and flexible tool for the stochastic simulation of crossing and selection on real genotypes. It is fully written in C for high execution speeds, has minimal dependencies, and is available as an R package for integration with Rs broad range of analysis and visualisation tools. Comparisons of a simulated recreation of a breeding program to the real data shows that the tools simulated offspring correctly show key population features, such as linkage disequilibrium patterns and genomic relationships. Both versions of genomicSimulation are freely available on GitHub: The R package version at https://github.com/vllrs/genomicSimulation/ and the C library version at https://github.com/vllrs/genomicSimulationC/

category: genetics


title: Protective effect of peroxisome proliferator-activated receptor- in regulating the early inflammation response to extended hepatectomy from a comparative transcriptomic analysis

abstract: The regulation mechanism of small-for-size syndrome remains unclear. Thus, we aimed to analyze the molecular profiles following extended hepatectomy and identify the therapeutic target. Major hepatectomy and extended hepatectomy were performed in the rat model, and the remnant livers were obtained dynamically for the high-throughput transcriptome analysis to identify the differentially expressed genes (DEGs). The general framework for weighted gene co-expression network analysis (WGCNA) was employed to explore the expression patterns of DEGs. As result, WGCNA identified 10 distinct gene co-expression modules according to the correlation between module eigengene and different postoperative time-points. The magenta module (gene count: 289) and the lightcyan module (gene count: 484) were found positively correlated with major hepatectomy instead of extended hepatectomy. In the lightcyan module, peroxisome proliferator-activated receptor- (PPAR) was selected and found the down-regulation in the remnant liver following extended (marginal) hepatectomy in rats and humans. Besides, administration of PPAR agonist attenuated hepatic inflammation injury while PPAR antagonist increased liver inflammation injury after extended hepatectomy in rats, marked by the significantly changed aminotransferases, tumor necrosis factor- and interleukin-6 levels in the plasm, and histological Suzuki criteria. Consequently, DEGs and their molecular profiles after extended hepatectomy were identified, and PPAR might be a potential therapy target for small-for-size syndrome.


title: Multiscale mechanical models of DNA nanotubes: A theoretical basis for the design and tuning of DNA nanocarriers

abstract: DNA nanostructures are one of potential candidates for drug carriers due to their good biocompatibility and non-specificity in vivo. A reliable prediction about mechanical properties of artificial DNA structures is desirable to improve the efficiency of DNA drug carriers, however there is only a handful of information on mechanical functionalities of DNA nanotubes (DNTs). This paper focuses on quantifying the multiscale correlations among DNT deformation, packaging conditions and surrounding factors to tune mechanical properties of DNTs. By combining WLC statistical mechanics model, Parsegians mesoscopic liquid crystal model and Eulers continuum beam theory, we developed a multiscale DNA-frame model; then theoretically characterize the initial packed states of DNTs for the first time, and reveal the diversity mechanism in mechanical properties of DNTs induced by interchain interactions and initial packed states. Moreover, the study of parameters, such as packaging conditions and environmental factors, provides a potential control strategy for tuning mechanical properties of DNTs. These conclusions provide a theoretical basis for accurately controlling the property and deformation of DNT in various DNT dynamic devices, such as DNA nanocarriers.


title: The Prion Protein Octarepeat Domain Forms Transient β-sheet Structures Upon Residue-Specific Cu(II) and Zn(II) Binding

abstract: Misfolding of the cellular prion protein (PrPC) is associated with the development of fatal neurodegenerative diseases called transmissible spongiform encephalopathies (TSEs). Metal ions appear to play a crucial role in the protein misfolding, and metal imbalance may be part of TSE pathologies. PrPC is a combined Cu(II) and Zn(II) metal binding protein, where the main metal binding site is located in the octarepeat (OR) region. Here, we used biophysical methods to characterize Cu(II) and Zn(II) binding to the isolated OR region. Circular dichroism (CD) spectroscopy data suggest that the OR domain binds up to four Cu(II) ions or two Zn(II) ions. Upon metal binding, the OR region seems to adopt a transient antiparallel {beta}-sheet hairpin structure. Fluorescence spectroscopy data indicates that under neutral conditions, the OR region can bind both Cu(II) and Zn(II) ions, whereas under acidic conditions it binds only Cu(II) ions. Molecular dynamics simulations suggest that binding of both metal ions to the OR region results in formation of {beta}-hairpin structures. As formation of {beta}-sheet structures is a first step towards amyloid formation, we propose that high concentrations of either Cu(II) or Zn(II) ions may have a pro-amyloid effect in TSEs.


title: Dynamics of CTCF and cohesin mediated chromatin looping revealed by live-cell imaging

abstract: Animal genomes are folded into loops and topologically associating domains (TADs) by CTCF and cohesin, but whether these loops are stable or dynamic is unknown. Here, we directly visualize chromatin looping at the Fbn2 TAD in mouse embryonic stem cells using super-resolution live-cell imaging and quantify looping dynamics by Bayesian inference. Our results are consistent with cohesin-mediated loop extrusion in cells, and with CTCF both stopping and stabilizing cohesin. Surprisingly, the Fbn2 loop is both rare and dynamic, with a looped fraction of ~3-6.5% and a median loop lifetime of ~10-30 minutes. Instead of a stable loop, our results establish a highly dynamic view of TADs and loops where the Fbn2 TAD exists predominantly in a partially extruded conformation. This dynamic and quantitative view of TADs may facilitate a mechanistic understanding of their functions.


title: Biocatalytic Formation of Novel Polyesters with para-Hydroxyphenyl groups in the Backbone - Engineering Cupriavidus necator for production of high-performance materials from CO2 and electricity

abstract:
Synthetic materials are integral components of consumables and durable goods and indispensable in our modern world. Polyesters are the most versatile bulk- and specialty-polymers, but their production is not sustainable, and their fate at end-of-life of great concern. Bioplastics are highly regarded alternatives but have shortcomings in material properties and commercial competitiveness with conventional synthetic plastics. These constraints have limited the success in global markets. Enabling bio-production of advanced bioplastics with superior properties from waste-derived feedstocks could change this.

We have created microbial cell factories that can produce a range of aliphatic and aromatic polyesters. A {Delta}phaC1 mutant of Cupriavidus necator H16 was complemented with hydroxyacyl-CoA transferases from either Clostridium propionicum (pct540) or Clostridium difficile (hadA), respectively. These were combined with a mutant PHA synthase (phaC1437) from Pseudomonas sp. MBEL 6-19, which rescued the PHA- phenotype of the knock-out mutant and allowed polymerization of various hydroxy carboxylates, including phloretic acid. This is the first-time, incorporation of an aromatic ring in the backbone of a biological polyester was achieved. Polymers contain para-hydroxyphenyl subunits are structurally analogous to synthetic aromatic polyesters like PET and high-strength polyarylates.

In a further advance, the transgenic strain was cultivated in a bio-electrochemical system under autotrophic conditions, enabling synthesis of aromatic bio-polyesters from H2 and O2 generated in situ, while assimilating CO2. Follow-up elementary flux-mode analysis established the feasibility of de novo production of twenty different polyesters from five different carbon- and energy-sources. This comprehensive study opens the door to sustainable bio-production of high-performance thermoplastics and thermosets.

Significance statementNew biomaterials can facilitate transition to a carbon-neutral chemical industry and a circular economy while at the same time preventing accumulation of plastic wastes in the environment. This can be accomplished by developing "drop-in" replacements for existing fossil carbon-based plastics. To that end, this work demonstrates for the first time that biocatalytic polymerization can incorporate an aromatic para-hydroxyphenyl carbonic acid into the backbone of a bio-polyester. This was accomplished using a genetically engineered microbial cell factory that assimilates carbon dioxide using hydrogen that is produced electrochemically in situ. A bio-electrochemical system eliminates the need for external supply of oxyhydrogen and avoids explosive mixtures, opening the door for sustainable production of biomaterials analogous to aromatic bulk-polyesters such as PET and high-performance "liquid-crystal polymers".

Highlights- Biocatalytic formation of aromatic polyesters with structural analogy to polyarylates
- Production of novel PHAs from CO2 and H2+O2 produced in situ in a bio-electrochemical system
- Expression-level of PHA synthase and molecular weight of polyesters are inversely correlated
- In silico design and analysis of pathways towards novel PHAs from different carbon-sources

category: synthetic biology


title: In vitro validation of neoantigen prediction algorithm for developing personalized cancer vaccine therapy

abstract:
BackgroundThe development of personalized neoantigen-based therapeutic cancer vaccines relies on computational algorithm-based pipelines. One of the critical issues in the pipeline is obtaining higher positive predictive value (PPV) performance, i.e., how many are immunogenic when selecting the top 5 to 20 candidate neoepitopes for the vaccination. We attempted to test the PPV of a neoepitope prediction algorithm Neopepsee.

MethodsSix breast cancer patients and patient-derived xenografts from three lung cancer patients and their paired peripheral blood samples were subjected to whole-exome and RNA sequencing. Neoantigen was predicted using two different algorithms (Neopepsee and pVACseq). Response of induced memory T cells to neopeptide candidates was evaluated by IFN-{gamma} Enzyme-linked immune absorbent spot (ELISpot) assays of peripheral blood mononuclear cell (PBMC) from three HLA-matched donors. Positive ELISpot response to a candidate peptide in at least 2 of 3 donor PBMC was regarded as an immunogenic response.

ResultsNeopepsee predicted 159 HLA-A matched neoepitope candidates out of 898 somatic mutations in nine patients (six breast and three lung cancer patients), whereas pVACseq predicted 84 HLA-A matched candidates. A total of 26 neopeptide candidates overlapped between the two predicted candidate pools. Among the candidates, 28 (20%, 28/ 137) and 15 (20%, 15/ 75) were positive by ELISpot assay, respectively. Among 26 overlapped candidates, 20 could be tested, and 7 of them (35%) were validated by ELISpot. Neopepsee identified at least one neoepitope in 7 of 9 patients (range 0-6), compared to 6 by pVACseq (range 0-5).

ConclusionAs suggested by Tumor Neoantigen Selection Alliance (TESLA), our results demonstrate low PPV of individual prediction models as well as the complementary nature of the Neopepsee and pVACseq and may help design neoepitope targeted cancer vaccines. Our data contribute a significant addition to the database of tested neoepitope candidates that can be utilized to further train and improve the prediction algorithms.


title: Age dependent changes in circulating Tfh cells influence the development of functional antibodies to malaria in children.

abstract: T-follicular helper (Tfh) cells are key drivers of antibodies that protect from malaria. However, little is known regarding the host and parasite factors that influence Tfh and functional antibody development. Here, we use samples from a large cross-sectional study of children residing in an area of high malaria transmission in Uganda to characterize Tfh cells and functional antibodies to multiple parasites stages. We identify a dramatic re-distribution of the Tfh cell compartment with age that is independent of malaria exposure, with Th2-Tfh cells predominating in early childhood, while Th1-Tfh cell gradually increase to adult levels over the first decade of life. Functional antibody acquisition is age-dependent and hierarchical acquired based on parasite stage, with merozoite responses followed by sporozoite and gametocyte antibodies. Antibodies were boosted in children with current infection, and were higher in females. The children with the very highest antibody levels had increased Tfh cell activation and proliferation, consistent with a key role of Tfh cells in antibody development. Together, these data reveal a complex relationship between the circulating Tfh compartment, antibody development and protection from malaria.


title: Disruption of c-di-GMP signaling networks unlocks cryptic expression of secondary metabolites during biofilm growth in Burkholderia pseudomallei

abstract:
The regulation and production of secondary metabolites during biofilm growth of Burkholderia spp. is not well understood. To learn more about the crucial role and regulatory control of cryptic molecules produced during biofilm growth, we disrupted c-di-GMP signaling in Burkholderia pseudomallei, a soil-borne bacterial saprophyte and the etiologic agent of melioidosis. Our approach to these studies combined transcriptional profiling with genetic deletions that targeted key c-di-GMP regulatory components to characterize responses to changes in temperature. Mutational analyses and conditional expression studies of c-di-GMP genes demonstrates their contribution to phenotypes such as biofilm formation, colony morphology, motility, and expression of secondary metabolite biosynthesis when grown as a biofilm at different temperatures. RNA-seq analysis was performed at varying temperatures in a {Delta}II2523 mutant background that is responsive to temperature alterations resulting in hypo- and hyper-biofilm forming phenotypes. Differential regulation of genes was observed for polysaccharide biosynthesis, secretion systems, and nonribosomal peptide and polyketide synthase (NRPS/PKS) clusters in response to temperature changes. Deletion mutations of biosynthetic gene clusters (BGCs) clusters 2, 11, 14 (syrbactin), and 15 (malleipeptin) in wild-type and {Delta}II2523 backgrounds also reveals the contribution of these BGCs to biofilm formation and colony morphology in addition to inhibition of Bacillus subtilis and Rhizoctonia solani. Our findings suggest that II2523 impacts the regulation of genes that contribute to biofilm formation and competition.

Characterization of cryptic BGCs under differing environmental conditions will allow for a better understanding of the role of secondary metabolites in the context of biofilm formation and microbe-microbe interactions.

ImportanceBurkholderia pseudomallei is a saprophytic bacterium residing in the environment that switches to a pathogenic lifestyle during infection of a wide range of hosts. The environmental cues that serve as the stimulus to trigger this change are largely unknown. However, it is well established that the cellular level of c-di-GMP, a secondary signal messenger, controls the switch from growth as planktonic cells to growth as a biofilm. Disrupting the signaling mediated by c-di-GMP allows for a better understanding of the regulation and the contribution of the surface associated and secreted molecules that contribute to the various lifestyles of this organism. The genome of B. pseudomallei also encodes cryptic biosynthetic gene clusters predicted to encode small molecules that potentially contribute to growth as a biofilm, adaptation, and interactions with other organisms. A better understanding of the regulation of these molecules is crucial to understanding how this versatile pathogen alters its lifestyle.


title: Insights into the role of lipoteichoic acids in Bacillus subtilis- a new function for MprF

abstract: In Bacillus subtilis, the cell is protected from the environment by a cell envelope, which comprises of layers of peptidoglycan that maintain the cell shape and anionic teichoic acids polymers whose biological function remains unclear. In B. subtilis, loss of all Class A Penicillin-Binding Proteins (aPBPs) which function in peptidoglycan synthesis is conditionally lethal. Here we show that this lethality is associated with an alteration of the lipoteichoic acids (LTA) and the accumulation of the major autolysin LytE in the cell wall. We provide the first evidence that the length and abundance of LTA acts to regulate the cellular level of LytE. Importantly, we identify a novel function for the aminoacyl-phosphatidylglycerol synthase MprF which acts to modulate LTA biosynthesis in B. subtilis and in the pathogen Staphylococcus aureus. This finding has implications for our understanding of antimicrobial peptide resistance (particularly daptomycin) in clinically relevant bacteria and MprF-associated virulence in pathogens, such as methicillin resistant S. aureus.


title: Scalable batch-correction method for integrating large-scale single-cell transcriptomes

abstract: Integration of the evolving large-scale single-cell transcriptomes requires scalable batch-correction approaches. Here we propose a simple batch-correction method that is scalable for integrating super large-scale single-cell transcriptomes from diverse sources. The core idea of the method is encoding batch information of each cell as a trainable parameter and added to its expression profile; subsequently, a contrastive learning approach is used to learn feature representation of the additive expression profile. We demonstrate the scalability of the proposed method by integrating 18 million cells obtained from the Human Cell Atlas. Our benchmark comparisons with current state-of-the-art single-cell integration methods demonstrated that our method could achieve comparable data alignment and cluster preservation. Our study would facilitate the integration of super large-scale single-cell transcriptomes. The source code is available at https://github.com/xilinshen/Fugue.


title: Single-molecule sizing through nano-cavity confinement

abstract: An approach relying on nano-cavity confinement is developed in this paper for the sizing of nanoscale particles and single biomolecules in solution. The approach, termed nano-cavity diffusional sizing (NDS), measures particle residence times within fluidic nano-cavities to determine their hydrodynamic radii. Using theoretical modelling and simulation, we show that the residence time of particles within nano-cavities above a critical timescale depends on the diffusion coefficient of the particle, which allows estimation of the particles size. We demonstrate this approach experimentally through measurement of particle residence times within nano-fluidic cavities using single-molecule confocal microscopy. Our data show that the residence times scale linearly with the sizes of nanoscale colloids, protein aggregates and single DNA oligonucleotides. NDS thus constitutes a new single molecule optofluidic approach that allows rapid and quantitative sizing of nanoscale objects for potential application in nanobiotechnology, biophysics, and clinical diagnostics.


title: Yield is negatively correlated with nucleotide-binding leucine-rich repeat gene content in soybean

abstract: The availability of increasing quantities of crop pangenome data permits the detailed association of gene content with agronomic traits. Here, we investigate disease resistance gene content of diverse soybean cultivars and report a significant negative correlation between the number of NLR resistance (R) genes and yield. We find no association between R-genes with seed weight, oil or protein content, and we find no correlation between yield and the number of RLK, RLP genes, or the total number of genes. These results suggest that recent yield improvement in soybean may be partially associated with the selective loss of NLR genes. Three quarters of soybean NLR genes do not show presence/absence variation, limiting the ability to select for their absence, and so the deletion or disabling of select NLR genes may support future yield improvement.

category: plant biology


title: Bacteria-derived peptidoglycan triggers an NF-kB dependent response in Drosophila gustatory neurons

abstract: Probing the external world is essential for eukaryotes to distinguish beneficial from pathogenic microorganisms. If it is clear that this task falls to the immune cells, recent work shows that neurons can also detect microbes, although the molecules and mechanisms involved are less characterized. In Drosophila, detection of bacteria-derived peptidoglycan by pattern recognition receptor (PRR) of the PGRP family expressed in immune cells, triggers NF-{kappa}B/IMD dependent signaling. We show here that one PGRP protein, called PGRP-LB, is expressed in some probosciss bitter taste neurons. In vivo calcium imaging reveals that the PGRP/IMD pathway is cell-autonomously required in these neurons to transduce the PGN signal. We finally show that NF-{kappa}B/IMD pathway activation in bitter neurons influences fly behavior. This demonstrates that flies use the same bacterial elicitor and signaling module to sense bacterial presence via the peripheral nervous system and trigger an anti-bacterial response in immune-competent cells.


title: Rapid generation of ventral spinal cord-like astrocytes from human iPSCs for modeling non-cell autonomous mechanisms of lower motor neuron disease

abstract: Astrocytes play important roles in the function and survival of neuronal cells. Dysfunctions of astrocytes are associated with numerous disorders and diseases of the nervous system, including motor neuron diseases such as amyotrophic lateral sclerosis (ALS). Human induced pluripotent stem cell (iPSC)-based approaches are becoming increasingly important for the study of the mechanisms underlying the involvement of astrocytes in non-cell autonomous processes of motor neuron degeneration in ALS. These studies must account for the molecular and functional diversity among astrocytes in different regions of the brain and spinal cord. It is essential that the most pathologically-relevant astrocyte preparations are used when investigating non-cell autonomous mechanisms of either upper or lower motor neuron degeneration in ALS. In this context, the main aim of this study was to establish conditions enabling rapid and robust generation of physiologically-relevant ventral spinal cord-like astrocytes that would provide an enhanced experimental model for the study of lower motor neuron degeneration in ALS. Neural progenitor cells with validated caudal and ventral features were derived from human iPSCs and differentiated into astrocytes, which were then characterized by examining morphology, markers of ventral spinal cord astrocytes, spontaneous and induced calcium transients, and astrogliosis markers. Efficient and streamlined generation of human iPSC-derived astrocytes with molecular and biological properties similar to physiological astrocytes in the ventral spinal cord was achieved. These induced astrocytes express markers of mature ventral spinal cord astrocytes, exhibit spontaneous and ATP-induced calcium transients, and lack signs of overt activation. Human iPSC- derived astrocytes with ventral spinal features offer advantages over more generic astrocyte preparations for the study of both ventral spinal cord astrocyte biology and the involvement of astrocytes in mechanisms of lower motor neuron degeneration in ALS.


title: Molecular Epidemiology and Microbiological Characteristics of Cryptococcus gattii VGII isolates from China

abstract:
Cryptococcus gattii (C. gattii) is a fungal pathogen that once caused an outbreak of cryptococcosis on Vancouver Island, and had spread worldwide, while few data were available in China. In this study, seven clinical isolates of C. gattii VGII were collected from 19 hospitals, Multi-locus Sequence Typing (MLST) analysis and whole-genome sequencing (WGS) was performed, and combined with published data for phylogenetic analysis. In addition, in vitro antifungal susceptibility testing, phenotypic analysis, and in vivo virulence studies were performed, subsequently, histopathological analysis of lung tissue was performed. C.gattii VGII infected patients were mainly immunocompetent male, and most of them had symptoms of central nervous system (CNS) involvement. MLST results showed that isolates from china exhibited high genetic diversity, and sequence type (ST) 7 was the major ST among the isolates. Some clinical isolates showed a close phylogenetic relationship with strains from Australia and South America. All clinical isolates did not show resistance to antifungal drugs. In addition, there was no correlation between virulence factors (temperature, melanin production, and capsule size) and virulence while in vivo experiments showed significant differences in virulence among strains. Lung fungal burden and damage to lung tissue correlated with virulence, and degree of damage to lung tissue in mice may highlight differences in virulence. Our work seeks to provide useful data for molecular epidemiology, antifungal susceptibility, and virulence differences of C. gattii VGII in China.

Author SummaryCryptococcus gattii is one of the life-threatening fungal pathogens that can infect immunocompetent individuals. C.gattii can be invaded through the respiratory tract and can be spread to the brain, causing meningoencephalitis and even death. In 1999, an outbreak of C. gattii VGII occurred on Vancouver Island in Canada and spread in the Pacific Northwest region, which attracted more attention. Compared with other nations, the data of C gattii VGII in China was limited. To better understand the molecular epidemiology and microbiological characteristics of C.gattii VGII isolates from China, Bioinformatics analysis and in vivo and in vitro experiments were performed for clinical strains. Multi-locus Sequence Typing and whole-genome sequencing reveal the genetic diversity of the isolates and genetic relationships with isolates from other countries. We also investigated the antifungal drug susceptibility and virulence differences of the isolates. The authors expect that our work can provide useful data for the prevention and control of C. gattii VGII in China and to draw attention to C. gattii.


title: Integrative interactome and ubiquitinome analyses reveal multiple regulatory pathways targeted by ICP0 in HSV-1 infected neuronal cells

abstract:
Herpes simplex virus 1 (HSV-1) is a neurotropic virus that can undergo both productive and latent infection in neurons. ICP0 is an HSV-1 E3 ubiquitin ligase crucial for productive infection and reactivation from latency. However, its targets have not been systematically investigated in neuronal cells. After confirming the importance of ICP0 in HSV-1 neuronal replication using an ICP0-null virus, we identified many ICP0-interacting proteins in infected neuronal and non-neuronal cells by mass-spectrometry-based interactome analysis. Co-immunoprecipitation assays validated ICP0 interactions with ACOT8, C1QBP, OTUD4, SNX9 and VIM in both Neuro-2a and 293T cells. Overexpression and knockdown experiments showed that SNX9 restricted replication of the ICP0-null but not wild-type virus in Neuro-2a cells. Ubiquitinome analysis by immunoprecipitating the trypsin digested ubiquitin reminant followed by mass spectrometry identified numerous candidate ubiquitination substrates of ICP0 in infected Neuro-2a cells, among which OTUD4 and VIM were novel substrates confirmed to be ubiquitinated by transfected ICP0 in Neuro-2a cells despite no evidence of their degradation by ICP0. Expression of OTUD4 was induced independently of ICP0 during HSV-1 infection. Overexpressed OTUD4 enhanced type I interferon expression during infection with the ICP0-null but not wild-type virus. In summary, by combining two proteomic approaches followed by confirmatory and functional experiments, we identified and validated multiple novel targets of ICP0 in neuronal cells, and revealed potential restrictive activities of SNX9 and OTUD4 as well as ICP0-dependent antagonism of these activities.

Author SummaryHerpes simplex virus 1 (HSV-1) establishes latent infection in neurons. ICP0 is known for its critical role in antagonizing cellular restrictive functions thereby initiating productive infection. It has been demonstrated to be important for both acute infection and reactivation from latency in neurons. However, little is known about its targets in neuronal cells. Here we combined two proteomic approaches, interactome and ubiquitinome analyses, to integratively identify interaction partners and substrates of ICP0 in HSV-1 infected neuronal cells. The results identified many novel targets as well as confirming previously reported ones. We also further validated some of the binding interactions and ubiquitin modifications. Functional studies revealed that the ICP0-interacting protein SNX9 restricted HSV-1 replication and the ICP0 substrate OTUD4 was induced to enhance type I interferon expression during HSV-1 neuronal infection. Moreover, the activities of these proteins appeared to be antagonized by ICP0-dependent mechanisms. This study provided comprehensive insight into ICP0 targets in neuronal cells and might prompt further investigation into the newly identified targets of ICP0.