title: Charged sequence motifs increase affinity towards liquid-liquid phase separation

abstract: Protein phase separation is a major governing factor in multiple cellular processes, such as those concerning RNA metabolism and RNA-binding proteins. Despite many key observations, the exact structural characteristics of proteins involved in the process are still not fully deciphered. In this work we show that proteins harbouring sequence regions with specific charged residue patterns are significantly associated with liquid-liquid phase separation. In particular, regions with repetitive arrays of alternating charges show the strongest association, whereas segments with generally high charge density and single -helices also show detectable but weaker connections.


title: Atria: An Ultra-fast and Accurate Trimmer for Adapter and Quality Trimming

abstract:
BackgroundAs Next Generation Sequencing takes a dominant role in terms of output capacity and sequence length, adapters attached to the reads and low-quality bases hinder the performance of downstream analysis directly and implicitly, such as producing false-positive single nucleotide polymorphisms (SNP), and generating fragmented assemblies. A fast trimming algorithm is in demand to remove adapters precisely, especially in read tails with relatively low quality.

FindingsWe present a trimming program named Atria. Atria matches the adapters in paired reads and finds possible overlapped regions with a super-fast and carefully designed byte-based matching algorithm (O(n) time with O(1) space). Atria also implements multi-threading in both sequence processing and file compression and supports single-end reads.

ConclusionsAtria performs favorably in various trimming and runtime benchmarks of both simulated and real data with other cutting-edge trimmers. We also provide an ultra-fast and lightweight byte-based matching algorithm. The algorithm can be used in a broad range of short-sequence matching applications, such as primer search and seed scanning before alignment.

Availability & ImplementationThe Atria executables, source code, and benchmark scripts are available at https://github.com/cihga39871/Atria under the MIT license.

category: bioinformatics


title: CardIAP: Calcium images analyzer web application

abstract:
MotivationIonic calcium (Ca2+) plays the role of the second messenger in eukaryotic cells associated with cellular functions of regulation of the cell cycle, such as transport, motility, gene expression, and metabolism (Permyakov and Kretsinger, 2009). The use of fluorometric techniques in isolated cells, loaded with Ca2+ sensitive fluorescent probes allows the quantitative measurement of dynamic events that occur in living, functioning cells. The Cardiomyocytes Images Analyzer Application (CardIAP) covers the need for tools to analyze and retrieve information from confocal microscopy images, in a systematic, accurate, and fast way.

ResultsHere we present the CardIAP web app, an automated method for the identification of spatio-temporal patterns in a calcium fluorescence imaging sequence. Through this tool, users can analyze single or multiple Ca2+ transients from confocal line-scan images and obtain quantitative information on the dynamic response of the stimulated myocyte.

Our web application also allows the user the extraction of data on calcium dynamics in downloadable tables and plots, simplifying the calculation of the alternation and discordance indices and their classification. CardIAP could assist in studying the underlying mechanisms of anomalous calcium release phenomena.

Availability and implementationCardIAP is an open-source app, entirely developed in Python, which can be freely accessed and used at http://cardiap.herokuapp.com/.


title: MOLECULAR VIDEOGAMING: SUPER-RESOLVED TRAJECTORY-BASED NANOCLUSTERING ANALYSIS USING SPATIO-TEMPORAL INDEXING

abstract: Single-molecule localization microscopy (SMLM) techniques are emerging as vital tools to unravel the nanoscale world of living cells. However, current analysis methods primarily focus on defining spatial nanoclusters based on detection density, but neglect important temporal information such as cluster lifetime and recurrence in "hotspots" on the plasma membrane. Spatial indexing is widely used in videogames to effectively detect interactions between moving geometric objects. Here, we use the R-tree spatial indexing algorithm to perform SMLM data analysis and determine whether the bounding boxes of individual molecular trajectories overlap, as a measure of their potential membership in nanoclusters. Extending the spatial indexing into the time dimension allows unique resolution of spatial nanoclusters into multiple spatiotemporal clusters. We have validated this approach using synthetic and SMLM-derived data. Quantitative characterization of recurring nanoclusters allowed us to demonstrate that both syntaxin1a and Munc18-1 molecules transiently cluster in hotspots on the neurosecretory plasma membrane, offering unprecedented insights into the dynamics of these protein which are essential to neuronal communication. This new analytical tool, named Nanoscale Spatiotemporal Indexing Clustering (NASTIC), has been implemented as a free and open-source Python graphic user interface.

category: biophysics


title: Combined free energy calculation and machine learning methods for understanding ligand unbinding kinetics

abstract: The determination of drug residence times, which define the time an inhibitor is in complex with its target, is a fundamental part of the drug discovery process. Synthesis and experimental measurements of kinetic rate constants are, however, expensive, and time-consuming. In this work, we aimed to obtain drug residence times computationally. Furthermore, we propose a novel algorithm to identify molecular design objectives based on ligand unbinding kinetics. We designed an enhanced sampling technique to accurately predict the free energy profiles of the ligand unbinding process, focusing on the free energy barrier for unbinding. Our method first identifies unbinding paths determining a corresponding set of internal coordinates (IC) that form contacts between the protein and the ligand, it then iteratively updates these interactions during a series of biased molecular-dynamics (MD) simulations to reveal the ICs that are important for the whole of the unbinding process. Subsequently, we performed finite temperature string simulations to obtain the free energy barrier for unbinding using the set of ICs as a complex reaction coordinate. Importantly, we also aimed to enable further design of drugs focusing on improved residence times. To this end, we developed a supervised machine learning (ML) approach with inputs from unbiased "downhill" trajectories initiated near the transition state (TS) ensemble of the string unbinding path. We demonstrate that our ML method can identify key ligand-protein interactions driving the system through the TS. Some of the most important drugs for cancer treatment are kinase inhibitors. One of these kinase targets is Cyclin Dependent Kinase 2 (CDK2), an appealing target for anticancer drug development. Here, we tested our method using two different CDK2 inhibitors for potential further development of these compounds. We compared the free energy barriers obtained from our calculations with those observed in available experimental data. We highlighted important interactions at the distal ends of the ligands that can be targeted for improved residence times. Our method provides a new tool to determine unbinding rates, and to identify key structural features of the inhibitors that can be used as starting points for novel design strategies in drug discovery.


title: Structural origins of Escherichia coli RNA polymerase open promoter complex stability

abstract: The first step of gene expression in all organisms requires opening the DNA duplex to expose one strand for templated RNA synthesis. In Escherichia coli, promoter DNA sequence fundamentally determines how fast the RNA polymerase (RNAP) forms "open" complexes (RPo), whether RPo persists for seconds or hours, and how quickly RNAP transitions from initiation to elongation. These rates control promoter strength in vivo but their structural origins remain largely unknown. Here we use cryo-electron microscopy to determine structures of RPo formed de novo at three promoters with widely differing lifetimes at 37{degrees}C: {lambda}PR (t1/2 [~] 10 hours), T7A1 (t1/2 [~] 4 minutes), and a point mutant in {lambda}PR ({lambda}PR-5C) (t1/2 [~] 2 hours). Two distinct RPo conformers are populated at {lambda}PR, likely representing productive and unproductive forms of RPo observed in solution studies. We find that changes in the sequence and length of DNA in the transcription bubble just upstream of the start site (+1) globally alter the network of DNA-RNAP interactions, base stacking, and strand order in the single-stranded DNA of the transcription bubble; these differences propagate beyond the bubble to upstream and downstream DNA. After expanding the transcription bubble by one base (T7A1), the nontemplate-strand "scrunches" inside the active site cleft; the template-strand bulges outside the cleft at the upstream edge of the bubble. The structures illustrate how limited sequence changes trigger global alterations in the transcription bubble that modulate RPo lifetime and affect the subsequent steps of the transcription cycle.


title: Interaction of hypoxia and nicotine acetylcholine receptor signaling network reveals a novel mechanism for lung adenocarcinoma progression in never-smokers

abstract: High incident of lung cancer among never smokers and their disease pathogenesis is an unexplained phenomenon. We have analyzed 1727 lung cancer patient data to understand the impact of smoking on overall survival of lung cancer patients and have observed a difference of only 47 days between smokers and never smokers in adenocarcinoma patients suggesting that the disease is equally fatal in never-smokers irrespective of gender. In this study, we have investigated the possible collaboration between the nAChR and hypoxia signaling pathway to elucidate a mechanism of disease progression in never-smokers. We report a previously unidentified increase in both acetylcholine and nAChR-7 levels in non small cell lung cancer cells in hypoxia. Similar increase in ubiquitously expressed nAChR-7 transcripts was also observed in other cancer lines. A direct binding of HIF-1 with the hypoxia response element (HRE) present at -48 position preceding the transcriptional start site in nAChR-7 promoter region was established. Significantly, the increased acetylcholine levels in hypoxia drove a feedback loop via modulation of PI3K/AKT pathway to stabilize HIF-1 in hypoxia. Further, Bungarotoxin, an antagonist of nAChR-7 significantly reversed hypoxia mediated metastasis and induction of HIF-1 in these cells. Our study gives a plausible explanation for the equally worse prognosis of lung adenocarcinoma in never-smokers wherein the nAChR signaling is enhanced in hypoxia by acetylcholine, in the absence of nicotine.

category: cancer biology


title: Flexible Path Planning through Vicarious Trial and Error

abstract: AO_SCPLOWBSTRACTC_SCPLOWFlexible planning is necessary for reaching goals and adapting when conditions change. We introduce a biologically plausible path planning model that learns its environment, rapidly adapts to change, and plans efficient routes to goals. Unlike prior models of hippocampal replay, our model addresses the decision-making process when faced with uncertainty. We tested the model in simulations of human and rodent navigation in mazes. Like the human and rat, the model was able to generate novel shortcuts, and take detours when familiar routes were blocked. Similar to rodent hippocampus recordings, the neural activity of the model resembles neural correlates of Vicarious Trial and Error (VTE) during early learning or during uncertain conditions. Similar to rodent studies, after learning, the neural activity resembles forward replay or preplay predicting a future route, and VTE activity decreases. We suggest that VTE, in addition to weighing possible outcomes, is a way in which an organism may gather information for future use.

category: neuroscience


title: Horizontal orientation facilitates pollinator attraction and rain avoidance in radially symmetrical flowers

abstract:
PremiseFloral angle, such as upward, horizontal, and downward orientation are known to evolve under both biotic and abiotic agents to enhance pollination success in zoophilious plants. Adaptive significance of horizontal orientation in radially symmetrical (actinomorphic) flowers under biotic and abiotic selection pressures were largely unknown, although those in bilaterally symmetrical flowers have been well studied.

MethodsUsing experimentally angle changed flowers, we examined the effects of flower angle on pollinator behaviors, pollination success and rain avoidance in a population of insect-pollinated Platycodon grandiflorus. We further investigated the frequency and amount of precipitation in the flowering season and pollen damage by water in this species. Main results: Horizontally oriented flowers received more visitations and pollen grains on the stigma in male and/or female phases than downward and/or upward oriented flowers and avoided pollen damage by rainfall compared to upward oriented flowers. The pollen germination experiment showed that approximately 30% of pollen grains burst in distilled water, thus pollen damage by rainfall was potentially serious in P. garndiflorus.

ConclusionIn this study, our field experiments revealed that upward flowers cannot avoid damage from rainfall during the flowering period whereas both upward and downward flowers suffered from pollinator limitation in female success. Thus, horizontal flower orientation is suggested to be adaptive in this insect-pollinated actinomorphic species which blooms in the rainy season.

category: ecology


title: The X-linked splicing regulator MBNL3 has been co-opted to restrict placental growth in eutherians

abstract: The eutherian placenta is a major site for parental genetic conflict. Here, we identify the X-linked Mbnl3 gene as a novel player in this dispute. Mbnl3 belongs to an RNA binding protein family whose members regulate alternative splicing and other aspects of RNA metabolism in association with cellular differentiation. We find that, in eutherians, Mbnl3 has become specifically expressed in placenta and has undergone accelerated sequence evolution leading to changes in its RNA binding specificities. Although its molecular roles are partly redundant with those of Mbnl2, Mbnl3 has also acquired novel biological functions. In particular, whereas Mbnl2;Mbnl3 double knockout mice display severe placental maturation defects leading to strong histological and functional abnormalities, Mbnl3 knockout alone results in increased placental growth and favors placental and fetal resource allocation during limiting conditions.

category: evolutionary biology


title: Effects of age and surface instability on the control of the center of mass

abstract: During standing, posture can be controlled by accelerating the Center of Mass (CoM) through shifting the center of pressure (CoP) within the base of support by applying ankle moments ("CoP mechanism"), or through the "counter-rotation mechanism", i.e., changing the angular momentum of segments around the CoM to change the direction of the ground reaction force. Postural control develops over the lifespan; at both the beginning and the end of the lifespan adequate postural control appears more challenging. In this study, we aimed to assess mediolateral balance performance and the related use of the postural control mechanisms in children, older adults and young adults when standing on different (unstable) surfaces. Sixteen pre-pubertal children (6-9y), 17 young adults (18-24y) and eight older adults (65-80y) performed bipedal upright standing trials of 16 seconds on a rigid surface and on three balance boards that could freely move in the frontal plane, varying in height (15-19 cm) of the surface of the board above the point of contact with the floor. Full body kinematics (16 segments, 48 markers, using SIMI 3D-motion analysis system (GmbH) and DeepLabCut and Anipose) were retrieved. Performance related outcome measures, i.e., the number of trials with balance loss and the Root Mean Square (RMS) of the time series of the CoM acceleration, the contributions of the CoP mechanism and the counter-rotation mechanism to the CoM acceleration in the frontal plane and selected kinematic measures, i.e. the orientation of the board and the head and the Mean Power Frequency (MPF) of the balance board orientation and the CoM acceleration were determined. Balance loss only occurred when standing on the highest balance board, twice in one older adult once in one young adult. In children and older adults, the RMS of the CoM accelerations were larger, corresponding to poorer balance performance. Across age groups and conditions, the contribution of the CoP mechanism to the total CoM acceleration was much larger than that of the counter-rotation mechanisms, ranging from 94%-113% vs 23%-38% (with totals higher than 100% indicating opposite effects of both mechanisms). Deviations in head orientation were small compared to deviations in balance board orientation. We hypothesize that the CoP mechanism is dominant, since the counter-rotation mechanism would conflict with stabilizing the orientation of the head in space.

category: bioengineering


title: ShareLoc - an open platform for sharing localization microscopy data

abstract: Novel insights and more powerful analytical tools can emerge from the reanalysis of existing data sets, especially via machine learning methods. Despite the widespread use of single molecule localization microscopy (SMLM) for super-resolution bioimaging, the underlying data are often not publicly accessible. We developed ShareLoc (https://shareloc.xyz), an open platform designed to enable sharing, easy visualization and reanalysis of SMLM data. We discuss its features and show how data sharing can improve the performance and robustness of SMLM image reconstruction by deep learning.


title: BITC induces cardiomyocyte proliferation and heart regeneration

abstract: Cardiomyocyte proliferation is an evolutionarily conserved mechanism that supports cardiac regeneration in vertebrates. Mammalian cardiomyocytes are arrested from the cell cycle shortly after birth, and therefore mammals lose the ability to regenerate injured myocardium for the rest of their lives. Pharmacological induction of cardiomyocyte proliferation has gained a lot of interest in recent years, as researchers strive to achieve heart tissue regeneration. Here we show that a small chemical, benzyl isothiocyanate (BITC), induced cardiomyocyte proliferation through activation of the cyclin-dependent kinase (CDK) pathway. BITC treatment also allowed heart regeneration in the infarcted neonatal heart, even after the regeneration period in mice. Furthermore, administration of BITC to adult mice in parallel with mild hypoxia (10% O2) induced cell cycle reentry and tissue regeneration in the adult heart. Our findings thus suggest that pharmacological activation of the CDK pathway using BITC, concurrently with the activation of hypoxia-related signaling pathways, may be a promising approach to inducing cardiac regeneration in patients with heart disease.

category: pharmacology and toxicology


title: Large dataset of octocoral mitochondrial genomes provides new insights into mt-mutS evolution and function

abstract: All studied octocoral mitochondrial genomes contain a gene from the MutS family, whose members code for proteins involved in DNA mismatch repair, other types of DNA repair, meiotic recombination, and other functions. Although mutS homologues are found in all domains of life as well as viruses, octocoral mt-mutS is the only such gene encoded in an organellar genome. While the function of mtMutS is not known, its domain architecture, conserved sequence, and presence of some characteristic residues suggest its involvement in mitochondrial DNA repair. This inference is supported by exceptionally low rates of mt-sequence evolution observed in octocorals. Previous studies of mt-mutS have been limited by the small number of octocoral mt-genomes available. We utilized sequence-capture data from the recent Quattrini et al. study to assemble complete mitochondrial genomes for 97 species of octocorals. Combined with sequences publicly available in GenBank, this resulted in a dataset of 184 complete mitochondrial genomes, which we used to re-analyze the conservation and evolution of mt-mutS. We discovered the first case of mt-mutS loss among octocorals in one of the two Pseudoanthomastus sp. assembled from Quattrini et al. data. This species displayed accelerated rate and and changed patterns of nucleotide substitutions in mt-genome, which we argue provide additional evidence for the role of mtMutS in DNA repair. In addition, we found accelerated mt-sequence evolution in the presence of mt-mutS in several octocoral lineages. This accelerated evolution did not appear to be the result of relaxed selection pressure and did not entail changes in patterns of nucleotide substitutions. Overall, our results support previously reported patterns of conservation in mt-mutS and suggest that mtMutS is involved in DNA repair in octocoral mitochondria. They also indicate that the presence of mt-mutS contributes to, but does not fully explain, the low rates of sequence evolution in octocorals


title: A thermodynamic model for interpreting tryptophan excitation-energy-dependent fluorescence spectra provides insight into protein conformational sampling and stability

abstract: It is now over thirty years since Demchenko and Ladokhin first posited the potential of the tryptophan red edge excitation shift (REES) effect to capture information on protein molecular dynamics. Whilst there have been many key efforts in the intervening years, a biophysical thermodynamic model to quantify the relationship between the REES effect and protein flexibility has been lacking. Without such a model the full potential of the REES effect cannot be realized. Here, we present a thermodynamic model of the protein REES effect that captures information on protein conformational flexibility, even with proteins containing multiple tryptophan residues. Our study incorporates exemplars at every scale, from tryptophan in solution, single tryptophan peptides to multi-tryptophan proteins, with examples including a structurally disordered peptide, de novo designed enzyme, human regulatory protein, therapeutic monoclonal antibody in active commercial development, and a mesophilic and hyperthermophilic enzyme. Combined, our model and data suggest a route forward for the experimental measurement of the protein REES effect and point to the potential for integrating bimolecular simulation with experimental data to yield novel insights.


title: Angiotensin-converting enzyme governs endogenous opioid signaling and synaptic plasticity in nucleus accumbens

abstract: Angiotensin-converting enzyme (ACE) regulates blood pressure by cleaving angiotensin I to produce angiotensin II. In the brain, ACE is expressed at uniquely high levels in the striatonigral pathway, but its central function remains poorly understood. We find that ACE degrades an unconventional enkephalin heptapeptide, Met-enkephalin-Arg-Phe, in the nucleus accumbens of mice. ACE inhibition enhanced mu opioid receptor activation by Met-enkephalin-Arg-Phe, causing a cell type-specific long-term depression of glutamate release onto medium spiny projection neurons expressing the Drd1 dopamine receptor. Systemic ACE inhibition was not intrinsically rewarding, but decreased the conditioned place preference caused by fentanyl administration, and enhanced reciprocal social interaction. Our results raise the enticing prospect that central ACE inhibition can boost endogenous opioid signaling for clinical benefit, while mitigating risk of addiction.


title: Multiple motors cooperate to establish and maintain acentrosomal spindle bipolarity in C. elegans oocyte meiosis

abstract: While centrosomes organize spindle poles during mitosis, oocyte meiosis can occur in their absence. Spindles in human oocytes frequently fail to maintain bipolarity and consequently undergo chromosome segregation errors, making it important to understand mechanisms that promote acentrosomal spindle stability. To this end, we have optimized the auxin-inducible degron system in C. elegans to remove factors from pre-formed oocyte spindles within minutes and assess effects on spindle structure. This approach revealed that dynein is required to maintain the integrity of acentrosomal poles; removal of dynein from bipolar spindles caused pole splaying, and when coupled with a monopolar spindle induced by depletion of kinesin-12 motor KLP-18, dynein depletion led to a complete dissolution of the monopole. Surprisingly, we went on to discover that following monopole disruption, individual chromosomes were able to reorganize local microtubules and re-establish a miniature bipolar spindle that mediated chromosome segregation. This revealed the existence of redundant microtubule sorting forces that are undetectable when KLP-18 and dynein are active. We found that the kinesin-5 family motor BMK-1 provides this force, uncovering the first evidence that kinesin-5 contributes to C. elegans meiotic spindle organization. Altogether, our studies have revealed how multiple motors are working synchronously to establish and maintain bipolarity in the absence of centrosomes.

category: cell biology


title: The function of α-Catenin mechanosensing in tissue morphogenesis

abstract: -catenin couples the cadherin-catenin complex to the actin cytoskeleton. The mechanosensitive -catenin M region undergoes conformational changes upon application of force to recruit binding partners. Here, we took advantage to the tension landscape in the Drosophila embryo to define three different states of -catenin mechanosensing in support of cell adhesion. Low, medium, and high tension contacts showed -catenin M region-dependent low, medium, and high levels of Vinculin and Ajuba recruitment. In contrast, Afadin/Canoe acts in parallel to -catenin at bicellular low and medium tension junctions, but requires an interaction with -catenin for its tension-sensitive enrichment at high-tension tricellular junctions. Individual M region domains make complex contributions to cell adhesion through their impact on binding partner recruitment, and redundancies with the function of Afadin/Canoe. Our data argue that -catenin and its interaction partners are part of a cooperative and partially redundant, mechanoresponsive network that supports AJs remodelling during morphogenesis.


title: Distal regulation, silencers and a shared combinatorial syntax are hallmarks of animal embryogenesis

abstract: The chromatin environment plays a central role in regulating developmental gene expression in metazoans. Yet, the basal regulatory landscape of metazoan embryogenesis is unknown. Here, we generate chromatin accessibility profiles for six embryonic, plus larval and adult stages in the sponge Amphimedon queenslandica. These profiles are reproducible within stages, reflect histone modifications, and identify transcription factor (TF) binding sequence motifs predictive of cis-regulatory elements during embryogenesis in other metazoans, but not the unicellular relative Capsaspora. Motif analysis of chromatin accessibility profiles across Amphimedon embryogenesis identifies three major developmental periods. As in bilaterian embryogenesis, early development in Amphimedon involves activating and repressive chromatin in regions both proximal and distal to transcription start sites. Transcriptionally repressive elements ( silencers) are prominent during late embryogenesis. They coincide with an increase in cis-regulatory regions harbouring metazoan TF binding motifs, and an increase in the expression of metazoan-specific genes. Changes in chromatin state and gene expression in Amphimedon suggest the conservation of distal enhancers, dynamically silenced chromatin, and TF-DNA binding specificity in animal embryogenesis.

category: developmental biology


title: Visceral mesoderm signaling regulates assembly position and function of the Drosophila testis niche

abstract: Tissue homeostasis often requires a properly placed niche to support stem cells. Morphogenetic processes that position a niche are just being described. For the Drosophila testis we recently showed that pro niche cells, specified at disparate positions during early gonadogenesis, must assemble into one collective at the anterior of the gonad. We now find that Slit and FGF signals emanating from adjacent visceral mesoderm regulate assembly. In response to signaling, niche cells express islet, which we find is also required for niche assembly. Without signaling, niche cells specified furthest from the anterior are unable to migrate, remaining dispersed. Function of such niches is severely disrupted, with niche cells evading cell cycle quiescence, compromised in their ability to signal the incipient stem cell pool, and failing to orient stem cell divisions properly. Our work identifies both extrinsic signaling and intrinsic responses required for proper assembly and placement of the testis niche.


title: The pseudobranch of jawed vertebrates is a mandibular arch-derived gill

abstract: The pseudobranch is a gill-like epithelial elaboration that sits behind the jaw of most fishes. This structure was classically regarded as a vestige of the ancestral gill-arch like condition of the gnathostome jaw. However, more recently, hypotheses of jaw evolution by transformation of a gill arch have been challenged, and the pseudobranch has alternatively been considered a specialised derivative of the second (hyoid) pharyngeal arch. Here, we demonstrate by cell lineage tracing in a cartilaginous fish, the skate (Leucoraja erinacea), that the pseudobranch does, in fact, derive from the mandibular arch, and that it shares gene expression features and cell types with gills. We also show that the mandibular arch pseudobranch is supported by a spiracular cartilage that is patterned by a shh-expressing epithelial signalling centre. This closely parallels the condition seen in the gill arches, where cartilaginous appendages called branchial rays supporting the respiratory lamellae of the gills are patterned by a shh-expressing gill arch epithelial ridge (GAER). Taken together, these findings support serial homology of the pseudobranch and gills, and an ancestral origin of gill arch-like anatomical features from the gnathostome mandibular arch.


title: SMAD4 target genes are part of a transcriptional network that integrates the response to BMP and SHH signaling during early limb bud patterning

abstract: SMAD4 regulates gene expression in response to BMP and TGF{beta} signal transduction and is required for diverse morphogenetic processes, but its target genes have remained largely elusive. Here, we use an epitope-tagged Smad4 allele for ChIP-seq analysis together with transcriptome analysis of wild-type and mouse forelimb buds lacking Smad4 in the mesenchyme. This analysis identifies the SMAD4 target genes during establishment of the feedback signaling system and establishes that SMAD4 predominantly mediates BMP signal-transduction during early limb bud development. Unexpectedly, the initial analysis reveals that the expression of cholesterol biosynthesis enzymes is precociously down-regulated and intracellular cholesterol levels reduced in Smad4-deficient limb bud mesenchymal progenitors. The SMAD4 target GRNs includes genes, whose expression in the anterior limb bud is up-regulated by interactions of SMAD4 complexes with enhancers active in the anterior mesenchyme. This reveals a predominant function of SMAD4 in up-regulating target gene expression in the anterior limb bud mesenchyme. Analysis of differentially expressed genes that are shared between Smad4- and Shh-deficient limb buds corroborates the positive role of SMAD4 in transcriptional regulation of anterior genes and reveals a repressive effect on posterior genes that are positively regulated by SHH signaling. This analysis uncovers the overall opposing effects of SMAD4-mediated BMP and SHH signalling on transcriptional regulation during early limb bud development. In summary, this analysis indicates that during early digit patterning and limb bud outgrowth, the anterior/proximal and proximo/distal expression dynamics of co-regulated genes are controlled by distinct and contrasting trans-regulatory inputs from SHH and SMAD4-mediated BMP signal transduction.


title: Distinct proteostasis states drive pharmacologic chaperone susceptibility for Cystic Fibrosis Transmembrane Conductance Regulator misfolding mutants

abstract: Pharmacological chaperones represent a class of therapeutic compounds for treating protein misfolding diseases. One of the most prominent examples is the FDA-approved pharmacological chaperone lumacaftor (VX-809), which has transformed cystic fibrosis (CF) therapy. CF is a fatal disease caused by mutations in the cystic fibrosis transmembrane conductance regulator (CFTR). VX-809 corrects folding of F508del CFTR, the most common patient mutation, yet F508del exhibits only mild VX-809 response. In contrast, rarer mutations P67L and L206W are hyper-responsive to VX-809, while G85E is non-responsive. Despite the clinical success of VX-809, the mechanistic origin for the distinct susceptibility of mutants remains unclear. Here, we use interactomics to characterize the impact of VX-809 on proteostasis interactions of P67L and L206W and compare these to F508del and G85E. We determine hyper-responsive mutations P67L and L206W exhibit decreased interactions with proteasomal, and autophagy degradation machinery compared to F508del and G85E. We then show inhibiting the proteasome attenuates P67L and L206W VX-809 response, and inhibiting the lysosome attenuates F508del VX-809 response. Our data suggests a previously unidentified but required role for protein degradation in VX-809 correction. Furthermore, we present an approach for identifying proteostasis characteristics of mutant-specific therapeutic response to pharmacological chaperones.


title: Digital Spatial Profiling of Collapsing Glomerulopathy

abstract: Collapsing glomerulopathy is a histologically distinct variant of focal and segmental glomerulosclerosis that presents with heavy proteinuria and portends a poor prognosis. Collapsing glomerulopathy can be triggered by viral infections such as HIV and SARS-CoV-2. Transcriptional profiling of collapsing glomerulopathy lesions is difficult since only a few glomeruli may exhibit this histology within a kidney biopsy and the mechanisms driving this heterogeneity are unknown. Therefore, we used recently developed digital spatial profiling (DSP) technology which permits quantification of mRNA at the level of individual glomeruli. Using DSP, we profiled 1,852 transcripts in glomeruli from HIV and SARS-CoV-2 infected patients with biopsy confirmed collapsing glomerulopathy. The increased resolution of DSP uncovered heterogeneity in glomerular transcriptional profiles that were missed in early laser capture microdissection studies of pooled glomeruli. Focused validation using immunohistochemistry and RNA in situ hybridization showed good concordance with DSP results. Therefore, DSP represents a powerful method to dissect transcriptional programs of pathologically discernible kidney lesions.

category: pathology


title: A thermostable oral SARS-CoV-2 vaccine induces mucosal and protective immunity

abstract: An ideal protective vaccine against SARS-CoV-2 should not only be effective in preventing disease, but also in preventing virus transmission. It should also be well accepted by the population and have a simple logistic chain. To fulfill these criteria, we developed a thermostable, orally administered vaccine that can induce a robust mucosal neutralizing immune response. We used our platform based on retrovirus-derived enveloped virus-like particles (e-VLPs) harnessed with variable surface proteins (VSPs) from the intestinal parasite Giardia lamblia, affording them resistance to degradation and the triggering of robust mucosal cellular and antibody immune responses after oral administration. We made e-VLPs expressing various forms of the SARS-CoV-2 Spike protein (S), with or without membrane protein (M) expression. We found that prime-boost administration of VSP-decorated e-VLPs expressing a pre-fusion stabilized form of S and M triggers robust mucosal responses against SARS-CoV-2 in mice and hamsters, which translate into complete protection from a viral challenge. Moreover, they dramatically boosted the IgA mucosal response of intramuscularly injected vaccines. We conclude that our thermostable orally administered e-VLP vaccine could be a valuable addition to the current arsenal against SARS-CoV-2, in a stand-alone prime-boost vaccination strategy or as a boost for existing vaccines.

category: immunology


title: Heat efficiently inactivates coronaviruses inside vehicles

abstract: Heat is an established method to inactivate coronaviruses, and there is utility in using heat to reduce viral load on common touch points in vehicles exposed to a person shedding SARS-CoV-2. As SARS-CoV-2 is a Biosafety level (BSL)-3 pathogen, real world testing of heat as a sanitation method for public and private vehicles becomes a challenge, requiring a surrogate coronavirus that can be handled safely outside of a BSL-3 facility. In this study, we used Bovine Coronavirus (BCoV) as a surrogate for SARS-CoV-2 to test the efficacy of heat-based betacoronavirus inactivation. In vitro, a 30-minute exposure to 56{degrees}C completely inactivated BCoV in solution, and a 15-minute exposure reduced recovery of BCoV >1000-fold. When heated to 56{degrees}C for 15 minutes, the infectivity of BCoV spotted and dried on typical porous and non-porous automobile interior materials was reduced by 99 - 99.99%. When BCoV was spotted and dried on hard plastic (seat) material placed inside an out of service transit bus, 56{degrees}C heat for 30 minutes reduced BCoV infectivity 85 - 99.5%. Thus, 56{degrees}C is an accessible, rapid, and effective method to inactivate coronaviruses inside motor vehicles.

category: microbiology


title: SARS-CoV-2 triggers DNA damage response in Vero E6 cells

abstract: The novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus responsible for the current COVID-19 pandemic and has now infected more than 200 million people with more than 4 million deaths globally. Recent data suggest that symptoms and general malaise may continue long after the infection has ended in recovered patients, suggesting that SARS-CoV-2 infection has profound consequences in the host cells. Here we report that SARS-CoV-2 infection can trigger a DNA damage response (DDR) in African green monkey kidney cells (Vero E6). We observed a transcriptional upregulation of the Ataxia telangiectasia and Rad3 related protein (ATR) in infected cells. In addition, we observed enhanced phosphorylation of CHK1, a downstream effector of the ATR DNA damage response, as well as H2AX. Strikingly, SARS-CoV-2 infection lowered the expression of TRF2 shelterin-protein complex, and reduced telomere lengths in infected Vero E6 cells. Thus, our observations suggest SARS-CoV-2 may have pathological consequences to host cells beyond evoking an immunopathogenic immune response.


title: Biparental age effects in the burying beetle Nicrophorus vespilloides

abstract: Parental age at reproduction influences offspring size and survival by affecting prenatal and postnatal conditions in a wide variety of species, including humans. However, most investigations into this manifestation of ageing focus upon maternal age effects; the effects of paternal age and interactions between maternal and paternal age are often neglected. Furthermore, even when maternal age effects are studied, pre- and postnatal effects are confounded. Using a cross-fostered experimental design, we investigated the joint effects of paternal and pre- and postnatal maternal ages on numerous offspring outcomes in a laboratory population of a species of burying beetle, Nicrophorus vespilloides. When we correct our tests for significance for multiple comparisons, we found no clear evidence for any parental effect senescence acting on egg size, larval weight, or larval survival. Nor did we find a statistical effect of paternal or egg producer age on the outcomes of foster mothers as measured by weight change experienced during caregiving. These findings are consistent with recent negative results reported in a similar study of N. vespilloides maternal age effects while also expanding these to other offspring traits and to paternal age effects. We discuss how the peculiar life history of this species may promote selection to resist the evolution of parental age effects, and how this might have influenced our ability to detect senescence.


title: STRIDE: accurately decomposing and integrating spatial transcriptomics using single-cell RNA sequencing

abstract: The recent advances in spatial transcriptomics have brought unprecedented opportunities to understand the cellular heterogeneity in the spatial context. However, the current limitations of spatial technologies hamper the exploration of cellular localizations and interactions at single-cell level. Here, we present spatial transcriptomics deconvolution by topic modeling (STRIDE), a computational method to decompose cell types from spatial mixtures by leveraging topic profiles trained from single-cell transcriptomics. STRIDE accurately estimated the cell-type proportions and showed balanced specificity and sensitivity compared to existing methods. We demonstrate STRIDEs utility by applying it to different spatial platforms and biological systems. Deconvolution by STRIDE not only mapped rare cell types to spatial locations but also improved the identification of spatial localized genes and domains. Moreover, topics discovered by STRIDE were associated with cell-type-specific functions, and could be further used to integrate successive sections and reconstruct the three-dimensional architecture of tissues. Taken together, STRIDE is a versatile and extensible tool for integrated analysis of spatial and single-cell transcriptomics and is publicly available at https://github.com/wanglabtongji/STRIDE.


title: Feasibility of Administering Human Pancreatic Cancer Chemotherapy in Spontaneous Pancreatic Cancer Mouse Model

abstract:
BackgroundBoth modified FOLFIRINOX (mFFX) and gemcitabine/nab-paclitaxel chemotherapy regimens have been shown to improve clinical outcomes in patients with pancreatic cancer, and are often used interchangeably as the standard of care. Preclinical studies often do not use these regimens, since administering these multiagent approaches can be difficult. In this study, we assessed the feasibility of administering these two chemotherapy regimens in spontaneous pancreatic tumors using KPC mice with the ultimate goal of advancing preclinical studies.

MethodsKPC mice were created by breeding KrasLSL-G12D/+ to Trp53fl/fl;Ptf1Cre/+, resulting in KrasLSL-G12D/+;p53fl/+;Ptf1Cre/+ mice. At 14 weeks of age, mice were palpated for spontaneous tumor growth that was verified using ultrasounds. Mice with tumors under 15 mm in diameter were used. The mice were assigned to one of seven treatment regimens: 1 cycle of mFOLFIRINOX (FFX X1), 2 cycles of mFOLFIRINOX (FFX X2), 1 cycle of mFOLFIRINOX with 40Gy SBRT (FFX SBRT), 1 cycle of gemcitabine/nab-paclitaxel (GEM/AB X1), 2 cycles of gemcitabine/nab-paclitaxel (GEM/AB X2), 2 cycles of gemcitabine/nab-paclitaxel with 40Gy SBRT (GEM/AB SBRT), or no chemotherapy treatment (control).

ResultsIn total, 92 mice were included. The median OS in the FFX X2 group was slightly longer that the median OS in the FFX X1 group (15 days vs 11 days, P=0.003). Mice in the GEM/AB X2 group had longer OS when compared to mice in the GEM/AB X1 group (33.5 vs 13 days, P=0.001). Mice treated with chemotherapy survived longer than untreated control animals (median OS: 6.5 days, P <0.001). Moreover, in mice treated with chemotherapy, mice that received 2 cycles of GEM/AB X2 had the longest survival, while the FFX X1 group had the poorest OS (P <0.001). Lastly, chemotherapy followed by consolidative SBRT trended towards increased local control and survival.

ConclusionsWe demonstrate the utility and feasibility of clinicall relevant mFOLFIRINOX and gemcitabine/nab-paclitaxel in preclinical models of pancreatic cancer.


title: Feasibility of Administering Human Pancreatic Cancer Chemotherapy in a Spontaneous Pancreatic Cancer Mouse Model


title: (Un)common space in infant neuroimaging studies: a systematic review of infant templates

abstract: In neuroimaging, spatial normalization is an important step that maps an individuals brain onto a template brain permitting downstream statistical analyses. Yet, in infant neuroimaging, there remain several technical challenges that have prevented the establishment of a standardized template for spatial normalization. Thus, many different approaches are used in the literature. To quantify the popularity and variability of these approaches in infant neuroimaging studies, we performed a systematic review of infant MRI studies from 2000 to 2020. Here, we present results from 833 studies meeting inclusion criteria. Studies were classified into 1) processing data in single subject space, 2) using a predefined, or "off the shelf", template, 3) creating a study specific template or 4) using a hybrid of these methods. We found that across the studies in the systematic review, single subject space was the most used (no common space). This was the most used common space for DWI and structural MRI studies while fMRI studies preferred off the shelf atlases. We found a pattern such that more recently published studies are more commonly using off the shelf atlases. When considering special populations, preterm studies most used single subject space while, when no special populations were being analyzed, an off the shelf template was most common. The most used off the shelf templates were the UNC Infant Atlases (26.1%). Using a systematic review of infant neuroimaging studies, we highlight a lack of an established "standard" template brain in these studies.


title: Parkinson's Disease-Associated, Sex-specific Changes in DNA Methylation at PARK7 (DJ-1), ATXN1, SLC17A6, NR4A2, and PTPRN2 in Cortical Neurons

abstract: Evidence for epigenetic regulation playing a role in Parkinsons disease (PD) is growing, particularly for DNA modifications. Approximately 90% of PD cases are due to a complex interaction between age, genes, and environmental factors, and epigenetic marks are thought to mediate the relationship between aging, genetics, the environment, and disease risk. To date, there are a small number of published genome-wide studies of DNA modifications in PD, but none accounted for cell-type or sex in their analyses. Given the hetereogeneity of bulk brain tissue samples and known sex differences in PD risk, progression, and severity, these are critical variables to account for. In this first genome-wide analysis of DNA methylation in an enriched neuronal population from PD post-mortem parietal cortex, we report sex-specific PD-associated methylation changes in PARK7 (DJ-1), SLC17A6 (VGLUT2), PTPRN2 (IA-2{beta}), NR4A2 (NURR1), and other genes involved in developmental pathways, neurotransmitter packaging and release, and axon and neuron projection guidance.


title: Functional characterization of the biogenic amine transporter system on human macrophages

abstract: Monocyte-derived macrophages are key players in tissue homeostasis and disease regulated by a variety of signaling molecules. Recent literature has highlighted the ability for biogenic amines to regulate macrophage functions, but the mechanisms governing biogenic amine signaling on and around immune cells remains nebulous. In the central nervous system, biogenic amine transporters are regarded as the master regulators of neurotransmitter signaling. While we and others have shown macrophages express these transporters, relatively little is known of their function on these cells. To address these knowledge gaps, we interrogated the function of norepinephrine (NET) and dopamine (DAT) transporters on human monocyte-derived macrophages. We found that both NET and DAT are present and can uptake substrate from the extracellular space at baseline. Not only was DAT expressed in cultured macrophages, but it was also detected in a subset of intestinal macrophages in situ. Surprisingly, we discovered a NET-independent, DAT-mediated immuno-modulatory mechanism in response to lipopolysaccharide (LPS). LPS induced reverse transport of dopamine through DAT, engaging autocrine/paracrine signaling loop that regulated the macrophage response. Removing this signaling loop enhanced the pro-inflammatory response to LPS. Finally, we found that this DAT-immune axis was disrupted in disease. Collectively, our data introduce a novel role for DAT in the regulation of innate immunity during health and disease.


title: The gut microbial signature of gestational diabetes mellitus and the association with diet intervention

abstract:
Gestational diabetes mellitus (GDM) is a high-risk pregnancy complication that is associated with metabolic disorder phenotypes, such as abnormal blood glucose and obesity. The link between microbiota and diet management contributes to metabolic homeostasis in GDM. Therefore, it is crucial to understand the structure of the gut microbiota in GDM and to explore the effect of dietary management on the microbiota structure. In this study, we analyzed the composition of the gut microbiota between 27 GDM and 30 healthy subjects at two time points using Illumina HiSeq 2500 platform. The taxonomy analyses suggested that the overall bacteria clustered by diabetes status, rather than diet intervention. Of particular interest, the phylum Acidobacteria in GDM was significantly increased, and positively correlated with blood glucose levels. Moreover, Partial least-squares discriminant analysis (PLS-DA) revealed that certain genera in the phyla Firmicutes, Bacteroidetes, Proteobacteria, and Lentisphaerae characterized the GDM gut microbiota. Correlation analysis indicated that blood glucose levels and BMI index were correlated with the relative abundance of SCFAS-producing genera. Through the comparison between the GDM and healthy samples with or without diet intervention, we discovered that the role of short-term diet management in GDM processes is associated with the change in the Firmicutes/Bacteroidetes ratio and some specific taxa, rather than an alternative gut microbial pattern. Our study have important implications for understanding the beneficial effects of diet intervention on the specific gut microbiota and thus possibly their metabolism in pregnant women with GDM.

ImportanceUnderstanding the composition and dynamics of the gut microbiota in GDM women under diet intervention is important because there may be opportunities for preventive strategies. We examined the relationships between GDM gut microbiota at two times before and after the diet intervention during second trimester of pregnancy and clinical characteristics in cohort of GDM women. We found that short-term diet management in GDM processes is associated with changes in the Firmicutes/Bacteroidetes ratio and some specific taxa rather than an alternative gut microbial pattern. Our study highlights the importance of considering diet intervention as the rescue of microbial dysfunction of GDM disease and can serve as a strategy for early prevention in future study.


title: Targeting Stem-loop 1 of the SARS-CoV-2 5'UTR to suppress viral translation and Nsp1 evasion

abstract: SARS-CoV-2 is a highly pathogenic virus that evades anti-viral immunity by interfering with host protein synthesis, mRNA stability, and protein trafficking. The SARS-CoV-2 nonstructural protein 1 (Nsp1) uses its C-terminal domain to block the mRNA entry channel of the 40S ribosome to inhibit host protein synthesis. However, how SARS-CoV-2 circumvents Nsp1-mediated suppression for viral protein synthesis and if the mechanism can be targeted therapeutically remain unclear. Here we show that N- and C-terminal domains of Nsp1 coordinate to drive a tuned ratio of viral to host translation, likely to maintain a certain level of host fitness while maximizing replication. We reveal that the SL1 region of the SARS-CoV-2 5 UTR is necessary and sufficient to evade Nsp1-mediated translational suppression. Targeting SL1 with locked nucleic acid antisense oligonucleotides (ASOs) inhibits viral translation and makes SARS-CoV-2 5 UTR vulnerable to Nsp1 suppression, hindering viral replication in vitro at a nanomolar concentration. Thus, SL1 allows Nsp1 to switch infected cells from host to SARS-CoV-2 translation, presenting a therapeutic target against COVID-19 that is conserved among immune-evasive variants. This unique strategy of unleashing a virus own virulence mechanism against itself could force a critical trade off between drug resistance and pathogenicity.

category: molecular biology


title: A Process-Based Model with Temperature, Water, and Lab-derived Data Improves Predictions of Daily Mosquito Density

abstract: While the number of human cases of mosquito-borne diseases has increased in North America in the last decade, accurate modeling of mosquito population density has remained a challenge. Longitudinal mosquito trap data over the many years needed for model calibration is relatively rare. In particular, capturing the relative changes in mosquito abundance across seasons is necessary for predicting the risk of disease spread as it varies from year to year. We developed a process-based mosquito population model that captures life-cycle egg, larva, pupa, adult stages, and diapause for Culex pipiens and Culex restuans mosquito populations. Mosquito development through these stages is a function of time, temperature, daylight hours, and aquatic habitat availability. The time-dependent parameters are informed by both laboratory studies and mosquito trap data from the Greater Toronto Area. The model incorporates city-wide water-body gauge and precipitation data as a proxy for aquatic habitat. This approach accounts for the nonlinear interaction of temperature and aquatic habitat variability on the mosquito life stages. We demonstrate that the full model predicts the yearly variations in mosquito populations better than a statistical model using the same data sources. This improvement in modeling mosquito abundance can help guide interventions for reducing mosquito abundance in mitigating mosquito-borne diseases like the West Nile virus.


title: Epidemiological and ecological consequences of virus manipulation of host and vector in plant virus transmission

abstract:
Many plant viruses are transmitted by insect vectors. Transmission can be described as persistent or non-persistent depending on rates of acquisition, retention, and inoculation of virus. Much experimental evidence has accumulated indicating vectors can prefer to settle and/or feed on infected versus noninfected host plants. For persistent transmission, vector preference can also be conditional, depending on the vectors own infection status. Since viruses can alter host plant quality as a resource for feeding, infection potentially also affects vector population dynamics. Here we use mathematical modelling to develop a theoretical framework addressing the effects of vector preferences for landing, settling and feeding - as well as potential effects of infection on vector population density - on plant virus epidemics. We explore the consequences of preferences that depend on the host (infected or healthy) and vector (viruliferous or nonviruliferous) phenotypes, and how this is affected by the form of transmission, persistent or non-persistent. We show how different components of vector preference have characteristic effects on both the basic reproduction number and the final incidence of disease. We also show how vector preference can induce bistability, in which the virus is able to persist even when it cannot invade from very low densities. Feedbacks between plant infection status, vector population dynamics and virus transmission potentially lead to very complex dynamics, including sustained oscillations. Our work is supported by an interactive interface https://plantdiseasevectorpreference.herokuapp.com/. Our model reiterates the importance of coupling virus infection to vector behaviour, life history and population dynamics to fully understand plant virus epidemics.

AUTHOR SUMMARYPlant virus diseases - which cause devastating epidemics in plant populations worldwide - are most often transmitted by insect vectors. Recent experimental evidence indicates how vectors do not choose between plants at random, but instead can be affected by whether plants are infected (or not). Virus infection can cause plants to "smell" different, because they produce different combinations of volatile chemicals, or "taste" different, due to chemical changes in infected tissues. Vector reproduction rates can also be affected when colonising infected versus uninfected plants. Potential effects on epidemic spread through a population of plants are not yet entirely understood. There are also interactions with the mode of virus transmission. Some viruses can be transmitted after only a brief probe by a vector, whereas others are only picked up after an extended feed on an infected plant. Furthermore there are differences in how long vectors remain able to transmit the virus. This ranges from a matter of minutes, right up to the entire lifetime of the insect, depending on the plant-virus-vector combination under consideration. Here we use mathematical modelling to synthesise all this complexity into a coherent theoretical framework. We illustrate our model via an online interface https://plantdiseasevectorpreference.herokuapp.com/.


title: Comprehensive analysis of DNA replication timing in genetic diseases and gene knockouts identifies MCM10 as a novel regulator of the replication program

abstract: Cellular proliferation depends on the accurate and timely replication of the genome. Several genetic diseases are caused by mutations in key DNA replication genes; however, it remains unclear whether these genes influence the normal program of DNA replication timing. Similarly, the factors that regulate DNA replication dynamics are poorly understood. To systematically identify trans-acting modulators of replication timing, we profiled replication in 184 cell lines from three cell types, encompassing 60 different gene knockouts or genetic diseases. Through a rigorous approach that considers the background variability of replication timing, we concluded that most samples displayed normal replication timing. However, mutations in two genes showed consistently abnormal replication timing. The first gene was RIF1, a known modulator of replication timing. The second was MCM10, a highly conserved member of the pre-replication complex. MCM10 mutant cells demonstrated replication timing variability comprising 46% of the genome and at different locations than RIF1 knockouts. Replication timing alterations in MCM10-mutant cells was predominantly comprised of replication initiation defects. Taken together, this study demonstrates the remarkable robustness of the human replication timing program and reveals MCM10 as a novel modulator of DNA replication timing.

category: genetics


title: Genome assembly and annotation of the tambaqui (Colossoma macropomum): an emblematic fish of the Amazon River basin

abstract: Colossoma macropomum known as "tambaqui" is the largest Characiformes fish in the Amazon River Basin and a leading species in Brazilian aquaculture and fisheries. Good quality meat and great adaptability to culture systems are some of its remarkable farming features. To support studies into the genetics and genomics of the tambaqui, we have produced the first high-quality genome for the species. We combined Illumina and PacBio sequencing technologies to generate a reference genome, assembled with 39X coverage of long reads and polished to a QV=36 with 130X coverage of short reads. The genome was assembled into 1,269 scaffolds to a total of 1,221,847,006 bases, with a scaffold N50 size of 40 Mb where 93% of all assembled bases were placed in the largest 54 scaffolds that corresponds to the diploid karyotype of the tambaqui. Furthermore, the NCBI Annotation Pipeline annotated genes, pseudogenes, and non-coding transcripts using the RefSeq database as evidence, guaranteeing a high-quality annotation. A Genome Data Viewer for the tambaqui was produced which benefits any groups interested in exploring unique genomic features of the species. The availability of a highly accurate genome assembly for tambaqui provides the foundation for novel insights about ecological and evolutionary facets and is a helpful resource for aquaculture purposes.

category: genomics


title: Radiotherapy induces specific miRNA expression profiles in glioblastoma exosomes

abstract:
BackgroundGlioblastoma (GBM) is the most aggressive and frequent primary brain tumor during adulthood. One of the major treatments is the association of surgery and a combination of chemo and radiotherapies. Despite its immediate efficiency, it fails to prevent the cancer recurrence in the irradiated area due to radioresistance mechanisms.

MicroRNAs (miRNAs or miR) are small non-coding, single strand RNA molecules encoding to various specific genes and able to regulate their expression and induce the tumor cell survival leading to radioresistance. Small extracellular vesicles (EVs), or exosomes released by tumor cells in tumor microenvironment and blood circulation are able to transport and diffuse miRNAs and affect the microenvironment by spreading the miRNAs, which drive radioresistance.

AimsO_LITo identify the variations of miRNAs expression induced by irradiation in human glioblastoma U87-MG cells and their secreted exosomes collected in supernatants.
C_LIO_LITo analyze the miRNAs variations in EVs-derived from the plasma of patients during radiotherapy, in order to identify a miRNA signature induced by radiotherapy in a liquid biopsy.
C_LI

Materiel and methodsU87-MG cells were cultured on plates and exposed to irradiation. miRNAs analyzes were performed in cells and in EVs isolated from cell supernatants to determine miRNAs expressions both in cells and in secreted exosomes before and after irradiation.

Plasma-derived EVs were collected from 4 glioblastoma patients before and after surgery and radiotherapy treatments.

ConclusionThe analysis of miRNAs expression profiles in both GBM cells and their derived EVs revealed that miR profile changes after irradiation. However, the number of similar miR between cells or EVs, following cell irradiation, was restricted to 3 miRs alone suggesting that the irradiation-induced changes in the miR profile in the cells and their EVs are not closely linked. In this context, the miR profile in EVs from patients plasma was investigated to establish a potential link with the miRNAs profile observed in EVs from irradiated cells and to assess its relationship with the response to radiotherapy. Three miRs (different from those identified in cells) were common between EVs derived from cells and patients derived-exosomes. These miRs detected in circulating EVs could provide a specific and reliable signature in response to ionizing radiation, which could be useful for monitoring the effectiveness of radiotherapy. Further experiments on a larger patients population with clinical data could also help to define whether this signature might have a prognostic value on the response to radiotherapy.


title: Specialization of the photoreceptor transcriptome by Srrm3-dependent microexons is required for outer segment maintenance and vision

abstract: Retinal photoreceptors differ in their transcriptomic profiles from other neuronal subtypes, likely as a reflection of their unique cellular morphology and function in the detection of light thorough the ciliary outer segment. We discovered a new layer of this molecular specialization by revealing that the vertebrate retina expresses the largest number of tissue-enriched microexons of all tissue types. A subset of these microexons is included exclusively in photoreceptor transcripts, particularly in genes involved in cilia biogenesis and in vesicle-mediated transport. This microexon program is regulated by Srrm3, a paralog of the neural microexon regulator Srrm4. Despite both proteins positively regulate retina microexons in vitro, only Srrm3 is highly expressed in mature photoreceptors and its deletion in zebrafish results in widespread downregulation of microexon inclusion, severe photoreceptor alterations and blindness. These results shed light into photoreceptors transcriptomic specialization and functionality, uncovering new cell type-specific roles for Srrm3 and microexons with implication for retinal diseases.

category: systems biology


title: TAAC - TMS Adaptable Auditory Control: a universal tool to mask TMS click

abstract:
BackgroundCoupling transcranial magnetic stimulation with electroencephalography (TMS-EEG) allows recording the EEG response to a direct, non-invasive cortical perturbation. However, obtaining a genuine TMS-evoked EEG potential requires controlling for several confounds, among which a main source is represented by the auditory evoked potentials (AEPs) associated to the TMS discharge noise (TMS click). This contaminating factor can be in principle prevented by playing a masking noise through earphones.

New methodHere we release TMS Adaptable Auditory Control (TAAC), a highly flexible, open-source, Matlab(R)-based interface that generates in real-time customized masking noises. TAAC creates noises starting from the stimulator-specific TMS click and tailors them to fit the individual, subject-specific click perception by mixing and manipulating the standard noises in both time and frequency domains.

ResultsWe showed that TAAC allows us to provide standard as well as customized noises able to effectively and safely mask the TMS click.

Comparison with existing methodsHere, we showcased two customized noises by comparing them to two standard noises previously used in the TMS literature (i.e., a white noise and a noise generated from the stimulator-specific TMS click only). For each, we quantified the Sound Pressure Level (SPL; measured by a Head and Torso Simulator - HATS) required to mask the TMS click in a population of 20 healthy subjects. Both customized noises were effective at safe (according to OSHA and NIOSH safety guidelines), lower SPLs with respect to standard noises.

ConclusionsAt odds with previous methods, TAAC allows creating effective and safe masking noises specifically tailored on each TMS device and subject. The combination of TAAC with tools for the real-time visualization of TEPs can help control the influence of auditory confounds also in non-compliant patients. Finally, TAAC is a highly flexible and open-source tool, so it can be further extended to meet different experimental requirements.


title: Urgency disrupts cognitive control of human action

abstract: Intelligent behavior requires cognitive control, the ability to act directed by goals despite competing action tendencies triggered by stimuli in the environment. For eye movements, it has recently been discovered that cognitive control is briefly disrupted in urgent situations (Salinas et al., 2019). In a time-window before an urgent response, participants could not help but look at a suddenly appearing visual stimulus, even though their goal was to look away from it. Urgency seemed to provoke a new visual-oculomotor phenomenon: A period in which saccadic eye movements are dominated by external stimuli, and uncontrollable by current goals. This period was assumed to arise from brain mechanisms controlling eye movements and spatial attention, such as those of the frontal eye field. Here, we show that the phenomenon is more general than previously thought. We found that urgency disrupted cognitive control also in well-investigated manual tasks, so that responses were dominated by goal-conflicting stimulus features. This dominance of behavior followed established trial-to-trial signatures of cognitive control that replicate across a variety of tasks. Thus together, these findings reveal that urgency temporarily impairs cognitive control in general, not only at brain mechanisms controlling eye movements.


title: A virus-encoded microRNA contributes to evade innate immune response during SARS-CoV-2 infection

abstract: SARS-CoV-2 infection results in impaired interferon response in severe COVID-19 patients. However, how SARS-CoV-2 interferes with host immune response is incompletely understood. Here, we sequenced small RNAs from SARS-CoV-2-infected human cells and identified a micro-RNA (miRNA) encoded in a recently evolved region of the viral genome. We show that the virus-encoded miRNA produces two miRNA isoforms in infected cells by the enzyme Dicer and they are loaded into Argonaute proteins. Moreover, the predominant miRNA isoform targets the 3UTR of interferon-stimulated genes and represses their expression in a miRNA-like fashion. Finally, the two viral miRNA isoforms were detected in nasopharyngeal swabs from COVID-19 patients. We propose that SARS-CoV-2 employs a virus-encoded miRNA to hijack the host miRNA machinery and evade the interferon-mediated immune response.


title: rbioacc: an R-package to analyse toxicokinetic data

abstract:
O_LI rbioacc is an  R-package dedicated to the analysis of experimental data collected from bioaccumulation tests during which organisms are exposed to a chemical (exposure phase) and then put into a clean media (depuration phase). Internal concentrations are regularly measured over time all along the experiment.
C_LIO_LI rbioacc provides ready-to-use functions to visualize and fully analyze such data. Under a Bayesian framework, this package fits a generic one-compartment toxicokinetic (TK) model automatically built from the data. It provides TK parameter estimates (appropriate uptake and elimination rates) and bioaccumulation metrics (e.g., BCF, BSAF, BMF). All parameter estimates, bioaccumulation metrics as well as predictions of internal concentrations into organisms are delivered with their uncertainty.
C_LIO_LIBioaccumulation metrics are primarily provided in support of environmental risk assessment, in full compliance with regulatory requirements required to approve marketing applications of chemical substances.
C_LIO_LIThis paper gives brief worked examples of the use of  rbioacc from data collected through standard bioaccumulation tests, and publicly available within the scientific literature. These examples constitute step-by-step user-guides to analyze any new data set, uploaded in the right format.
C_LI


title: `rbioacc`: an R-package to analyse toxicokinetic data


title: 'rbioacc': an R-package to analyse toxicokinetic data


title: Malayan kraits (Bungarus candidus) show affinity to anthropogenic structures in a human dominated landscape

abstract: Animal movement can impact human-wildlife conflict through a variety of features: increased movement can lead to greater chance of encounter, remaining still can lead to greater or lesser detection, and activity can modulate their impact on humans. Here we assess the movement and space use of the highly venomous and medically important Malayan krait (Bungarus candidus) on a suburban university campus. We radio- tracked 14 kraits for an average of 114 days (min: 19, max: 218), with an average of 106 fixes (min: 21, max: 229). We assessed movement pathways and activity with dynamic Brownian Bridge Movement models, evaluated patterns of revisitation to identify site fidelity, and determined whether individuals selected for particular land-use types with Integrated Step Selection Functions. Most individuals displayed some level of attraction to buildings (n = 10) and natural areas (n = 12); we identified a similar unambiguous pattern of attraction to buildings and natural areas at the population level (of our sample). Snakes remained in shelter sites for long periods of time (max = 94 days) and revisited sites on average every 15.45 days. Over 50% of fixes were within human settlements and 37.1% were associated with buildings. We found generally seasonal patterns of activity, with higher activity in wet seasons (when classes are typically in session on campus), and lower activity in the hot season (when there are a number of short breaks causing students to leave campus). These results show frequent proximity between Malayan kraits and humans at Suranaree University of Technology, Thailand; thereby, suggesting a near constant potential for human- wildlife conflict. Despite the fact that no snakebites from this species occurred at the university during our study period, substantial education and awareness training should be considered to raise awareness to ensure continued coexistence on campus.


title: Remotely detected plant function in two midwestern prairie grassland experiments reveals belowground processes

abstract: Imaging spectroscopy provides the opportunity to incorporate leaf and canopy optical data into ecological studies, but the extent to which remote sensing of vegetation can enhance the study of belowground processes is not well understood. In grassland systems, aboveground and belowground vegetation quantity and quality are coupled, and both influence belowground microbial processes and nutrient cycling, providing a potential link between remote sensing and belowground processes. We hypothesized that ecosystem productivity, and the chemical, structural and phylogenetic-functional composition of plant communities would be detectable with remote sensing and could be used to characterize belowground plant and soil processes in two grassland biodiversity experiments--the BioDIV experiment at Cedar Creek Ecosystem Science Reserve in Minnesota and the Wood River Nature Conservancy experiment in Nebraska. Specifically, we tested whether aboveground vegetation chemistry and productivity, as detected from airborne sensors, predict soil properties microbial processes and community composition. Imaging spectroscopy data were used to map aboveground biomass and green vegetation cover, functional traits and phylogenetic-functional community composition of vegetation. We examined the relationships between the image-derived variables and soil carbon and nitrogen concentration, microbial community composition, biomass and extracellular enzyme activity, and soil processes, including net nitrogen mineralization. In the BioDIV experiment--which has low overall diversity and productivity despite high variation in each--belowground processes were driven mainly by variation in the amount of organic matter inputs to soils. As a consequence, soil respiration, microbial biomass and enzyme activity, and fungal and bacterial composition and diversity were significantly predicted by remotely sensed vegetation cover. In contrast, at Wood River, where plant diversity and productivity were consistently higher, remotely sensed functional, chemical and phylogenetic composition of vegetation predicted belowground extracellular enzyme activity, microbial biomass, and net nitrogen mineralization rates, while aboveground biomass did not. The strong, contrasting associations between the quantity and chemistry of aboveground inputs with belowground soil processes and properties provide a basis for using imaging spectroscopy to understand belowground processes across productivity gradients in grassland systems. However, a mechanistic understanding of how above and belowground components interact among different ecosystems remains critical to extending these results broadly.


title: An adjuvanted SARS-CoV-2 RBD nanoparticle elicits neutralizing antibodies and fully protective immunity in aged mice

abstract: Development of affordable and effective vaccines that can also protect vulnerable populations such as the elderly from COVID-19-related morbidity and mortality is a public health priority. Here we took a systematic and iterative approach by testing several SARS-CoV-2 protein antigens and adjuvants to identify a combination that elicits neutralizing antibodies and protection in young and aged mice. In particular, SARS-CoV-2 receptorbinding domain (RBD) displayed as a protein nanoparticle (RBD-NP) was a highly effective antigen, and when formulated with an oil-in-water emulsion containing Carbohydrate fatty acid MonoSulphate derivative (CMS) induced the highest levels of cross-neutralizing antibodies compared to other oil-in-water emulsions or AS01B. Mechanistically, CMS induced antigen retention in the draining lymph node (dLN) and expression of cytokines, chemokines and type I interferon-stimulated genes at both injection site and dLN. Overall, CMS:RBD-NP is effective across multiple age groups and is an exemplar of a SARS-CoV-2 subunit vaccine tailored to the elderly.


title: Distinct, opposing functions for CFIm59 and CFIm68 in mRNA alternative polyadenylation of Pten and in the PI3K/Akt signalling cascade

abstract: The precise maintenance of PTEN dosage is crucial for tumor suppression across a wide variety of cancers. Post-transcriptional regulation of Pten heavily relies on regulatory elements encoded by its 3UTR. We previously reported the important diversity of 3UTR isoforms of Pten mRNAs produced through alternative polyadenylation (APA). Here, we reveal the direct regulation of Pten APA by the mammalian cleavage factor I (CFIm) complex, which in turn contributes to PTEN protein dosage. CFIm consists of the UGUA-binding CFIm25 and APA regulatory subunits CFIm59 or CFIm68. Deep sequencing analyses of perturbed (KO and KD) cell lines uncovered the differential regulation of Pten APA by CFIm59 and CFIm68 and further revealed that their divergent functions have widespread impact for APA in transcriptomes. Differentially regulated genes include numerous factors within the phosphoinositide 3-kinase (PI3K)/protein kinase B (Akt) signalling pathway that PTEN counter-regulates. We further reveal a stratification of APA dysregulation among a subset of PTEN-driven cancers, with recurrent alterations among PI3K/Akt pathway genes regulated by CFIm. Our results refine the transcriptome selectivity of the CFIm complex in APA regulation, and the breadth of its impact in PTEN-driven cancers.

category: biochemistry


title: B cell-intrinsic requirement for WNK1 kinase in T cell-dependent antibody responses

abstract:
Migration and adhesion play critical roles in B cells, regulating recirculation between lymphoid organs, migration within lymphoid tissue and interaction with CD4+ T cells. However, there is limited knowledge of how B cells integrate chemokine receptor and integrin signaling with B cell activation to generate efficient humoral responses. Here we show that the WNK1 kinase, a regulator of migration and adhesion, is essential in B cells for T-dependent antibody responses. We demonstrate that WNK1 transduces signals from the BCR, CXCR5 and CD40, and using intravital imaging we show that WNK1 regulates migration of naive and activated B cells, and their interactions with T cells. Unexpectedly, we show that WNK1 is required for BCR- and CD40-induced proliferation, acting through the OXSR1 and STK39 kinases, and for efficient B cell-T cell collaboration in vivo. Thus, WNK1 is critical for humoral immune responses, by regulating B cell migration, adhesion and T cell-dependent activation.

SummaryThe WNK1 kinase is essential in B cells for T-dependent antibody responses because it is activated by signaling from BCR, CXCR5 and CD40 and regulates B cell migration, adhesion, T-dependent activation, and differentiation into germinal center B cells and plasma cells.


title: SkewDB: A comprehensive database of GC and a 10 other skews for over 28,000 chromosomes and plasmids

abstract:
GC skew denotes the relative excess of G nucleotides over C nucleotides on the leading versus the lagging replication strand of eubacteria. While the effect is small, typically around 2.5%, it is robust and pervasive. GC skew and the analogous TA skew are a localized deviation from Chargaffs second parity rule, which states that G and C, and T and A occur with (mostly) equal frequency even within a strand.

Most bacteria also show the analogous TA skew. Different phyla show different kinds of skew and differing relations between TA and GC skew.

This article introduces an open access database (https://skewdb.org) of GC and 10 other skews for over 28,000 chromosomes and plasmids. Further details like codon bias, strand bias, strand lengths and taxonomic data are also included. The SkewDB database can be used to generate or verify hypotheses. Since the origins of both the second parity rule, as well as GC skew itself, are not yet satisfactorily explained, such a database may enhance our understanding of microbial DNA.


title: SkewDB: A comprehensive database of GC and 10 other skews for over 28,000 chromosomes and plasmids


title: Trans-cellular tunnels induced by the fungal pathogen Candida albicans facilitate invasion through successive epithelial cells without host damage

abstract: The opportunistic fungal pathogen Candida albicans is normally commensal, residing in the mucosa of most healthy individuals. In susceptible hosts, its filamentous hyphal form can invade epithelial layers leading to superficial or severe systemic infection. Invasion is mainly intracellular, though it causes no apparent damage to host cells. We investigated the invasive lifestyle of C. albicans in vitro using live-cell imaging and the damage-sensitive reporter galectin-3. Quantitative single cell analysis showed that invasion can result in host membrane breaching at different stages of invasion and cell death, or in traversal of host cells without membrane breaching. Membrane labelling and three-dimensional "volume" electron microscopy revealed that hyphae can traverse several host cells within trans-cellular tunnels that are progressively remodelled and may undergo  inflations linked to host glycogen stores, possibly during nutrient uptake. Thus, C. albicans invade epithelial tissues by either inflicting or avoiding host damage, the latter facilitated by trans-cellular tunnelling.


title: Variant surface glycoprotein expression in human African trypanosomiasis patients

abstract:
Trypanosoma brucei gambiense, an extracellular protozoan parasite, is the primary causative agent of human African trypanosomiasis. T. b. gambiense is endemic to West and Central Africa, where it is transmitted by the bite of infected tsetse flies. In the bloodstream of an infected host, the parasite evades antibody recognition by altering the Variant Surface Glycoprotein (VSG) that forms a dense coat on its cell surface through a process known as antigenic variation. Each VSG has a variable N-terminal domain exposed to the host and a less variable C-terminal domain that is partially hidden from host antibodies. Our lab developed VSG-seq, a targeted RNA-seq method, to study VSG gene expression in T. brucei. Studies using VSG-seq to characterize antigenic variation in a mouse model have revealed marked diversity in VSG expression within parasite populations, but this finding has not yet been validated in a natural human infection. Here, we used VSG-seq to analyze VSGs expressed in the blood of 12 patients infected with T. b. gambiense. The number of VSGs identified per patient ranged from one to 14; notably, two VSGs were shared by more than one patient. Analysis of expressed VSG N-terminal domain types revealed that 82% of expressed VSGs encoded a type B N-terminus, a bias not seen in datasets from other T. brucei subspecies. C-terminal types in T. b. gambiense infection were also restricted. These results demonstrate a bias either in the underlying VSG repertoire of T. b. gambiense or in this subspecies expression of VSGs during infection. This work demonstrates the feasibility of using VSG-seq to study antigenic variation in human infections and highlights the importance of understanding VSG repertoires in the field.

Author SummaryHuman African Trypanosomiasis is a neglected tropical disease primarily caused by the extracellular parasite Trypanosoma brucei gambiense. To avoid elimination by the host, these parasites repeatedly replace their dense surface coat of Variant Surface Glycoprotein (VSG). Despite the important role of VSGs in prolonging infection, VSG expression during natural human infections is poorly understood. A better understanding of natural VSG gene expression dynamics can clarify the mechanisms that T. brucei uses to alter its VSG coat and improve how trypanosomiasis is diagnosed in humans. We analyzed the expressed VSGs detected in the blood of patients with trypanosomiasis. Our findings indicate that parasites in both natural and experimental infections express diverse sets of VSGs.


title: Trypanosomal variant surface glycoprotein expression in human African trypanosomiasis patients


title: The Xer activation factor of TLC{Phi} expands the possibilities for Xer recombination

abstract: Many mobile elements take advantage of the highly-conserved chromosome dimer resolution system of bacteria, Xer. They participate in the transmission of antibiotic resistance and pathogenicity determinants. In particular, the toxin-linked cryptic satellite phage (TLC{Phi}) plays an essential role in the continuous emergence of new toxigenic clones of the Vibrio cholerae strain at the origin of the ongoing 7th cholera pandemic. The Xer machinery is composed of two chromosomally-encoded tyrosine recombinases, XerC and XerD. They resolve chromosome dimers by adding a crossover between sister copies of a specific 28 base pair site of bacterial chromosomes, dif. The activity of XerD depends on a direct contact with a cell division protein, FtsK, which spatially and temporally constrains the process. TLC{Phi} encodes for a XerD-activation factor (XafT), which drives the integration of the phage into the dif site of the primary chromosome of V. cholerae independently of FtsK. However, XerD does not bind to the attachment site (attP) of TLC{Phi}, which raised questions on the integration process. Here, we compared the integration efficiency of thousands of synthetic mini-TLC{Phi} plasmids harbouring different attP sites and assessed their stability in vivo. In addition, we compared the efficiency with which XafT and the XerD activation domain of FtsK drive recombination reactions in vitro. Taken together, our results suggest that XafT promotes the formation of synaptic complexes between canonical Xer recombination sites and imperfect sites.


title: Widespread roles of Trypanosoma brucei ATR in nuclear genome function and transmission are linked to R-loops

abstract: Inheritance of aberrant chromosomes can compromise genome integrity and affect cellular fitness. In eukaryotes, surveillance pathways and cell cycle checkpoints monitor for aberrant DNA transmission and the ATR kinase, a regulator of the DNA damage response, plays a pivotal role. Prior work revealed that ATR acts during antigenic variation in Trypanosoma brucei mammal-infective life cycle forms and that its loss is lethal, but how widely ATR operates in genome maintenance is largely unknown. Here, we show that after prolonged ATR depletion by RNAi T. brucei continues to synthesise DNA and enters new rounds of cell division, despite increased genome damage. Furthermore, we detect defective chromosome segregation,  micronuclei formation and disruption of the nuclear architecture. RNA-seq revealed that loss of ATR affects the expression of nearly half the genes in the genome, including both RNA Polymerase I and II transcription. Using ChIP-seq of yH2A and DRIP-seq, we reveal overlapping signals for genome damage and R-loops after ATR depletion in all intergenic regions. In addition, we report reduced R-loop levels and accumulation of yH2A signal within centromeres. Together, our data indicates widespread roles of ATR in T. brucei, including differing roles in R-loop homeostasis during multigene transcription and in chromosome segregation.


title: A broad-specificity O-glycoprotease that enables improved analysis of glycoproteins and glycopeptides containing intact complex O-glycans

abstract: Analysis of mucin type O-glycans linked to serine/threonine of glycoproteins is technically challenging, in part, due to a lack of effective enzymatic tools that enable their analysis. Recently, several O-glycan-specific endoproteases that can cleave the protein adjacent to the appended glycan have been described. Despite significant progress in understanding the biochemistry of these enzymes, known O-glycoproteases have specificity constrains, such as inefficient cleavage of glycoproteins bearing sialylated O-glycans, high selectivity for certain type of glycoproteins or protein sequence bias, that limit their analytical application. In this study, we examined the capabilities of an immunomodulating metalloprotease (IMPa) from Pseudomonas aeruginosa. The peptide substrate sequence selectivity and its impact on IMPa activity was interrogated using an array of synthetic peptides and their glycoforms. We show that IMPa has no specific P1 residue preference and can tolerate most amino acids at the P1 position, except aspartic acid. The enzyme does not cleave between two adjacent O-glycosites, indicating that O-glycosylated serine/threonine is not allowed at position P1. Glycopeptides with as few as two amino acids on either side of an O-glycosite were specifically cleaved by IMPa. Finally, IMPa efficiently cleaved peptides and proteins carrying sialylated and asialylated O-glycans of varying complexity. We present the use of IMPa in a one-step O-glycoproteomics workflow for glycoprofiling of individual purified glycoproteins granulocyte colony-stimulating factor (G-CSF) and receptor-type tyrosine-protein phosphatase C (CD45) without the need for glycopeptide enrichment. In these examples, IMPa enabled identification of O-glycosites and the range of complex O-glycan structures at each site.


title: Simultaneous Triple-Parametric Optical Mapping of Transmembrane Potential, Intracellular Calcium and NADH for Cardiac Physiology Assessment

abstract: Investigation of the complex relationships and dependencies of multiple cellular processes that govern cardiac physiology and pathophysiology requires simultaneous dynamic assessment of multiple parameters. In this study, we introduce triple-parametric optical mapping to simultaneously image metabolism, electrical excitation, and calcium signaling from the same field of view and demonstrate its application in the field of drug testing and cardiovascular research. We applied this metabolism-excitation-contraction coupling (MECC) methodology to test the effects of blebbistatin, 4-aminopyridine and verapamil on cardiac physiology. While blebbistatin and 4-aminopyridine alter multiple aspects of cardiac function suggesting off-target effects, the effects of verapamil were on-target and it altered only one of ten tested parameters. Triple-parametric optical mapping was also applied during ischemia and reperfusion; and we identified that metabolic changes precede the effects of ischemia on cardiac electrophysiology.


title: Diversity and dynamics of the CRISPR-Cas systems associated with Bacteroides fragilis in human population

abstract: CRISPR-Cas systems are adaptive immune systems commonly found in prokaryotes that provide sequence-specific defense against invading mobile genetic elements (MGEs). The memory of these immunological encounters are stored in CRISPR arrays, where spacer sequences record the identity and history of past invaders. Analyzing such CRISPR arrays provide insights into the dynamics of CRISPR-Cas systems and the adaptation of their host bacteria to rapidly changing environments such as the human gut. In this study, we utilized 601 Bacteroides fragilis genome isolates from 12 healthy individuals, 6 of which include longitudinal observations, and 222 available B. fragilis reference genomes to update the understanding of B. fragilis CRISPR-Cas dynamics and their differential activities. Analysis of longitudinal genomic data showed that some CRISPR array structures remained relatively stable over time whereas others involved radical spacer acquisition during some periods, and diverse CRISPR arrays (associated with multiple isolates) co-existed in the same individuals with some persisted over time. Furthermore, features of CRISPR adaptation, evolution, and microdynamics were highlighted through an analysis of host-MGE network, such as modules of multiple MGEs and hosts, reflecting complex interactions between B. fragilis and its invaders mediated through the CRISPR-Cas systems. This work demonstrates the power of using culture-based population genomics to reveal the activities and evolution of the CRISPR-Cas systems associated with gut bacteria in human population. We made available of all annotated CRISPR-Cas systems and their target MGEs, and their interaction network as a web resource at https://omics.informatics.indiana.edu/CRISPRone/Bfragilis.

category: bioinformatics


title: Investigation of gating in outer membrane porins provides new perspectives on antibiotic resistance mechanisms

abstract: Gram-negative bacteria pose a serious public health concern, primarily due to a higher frequency of antibiotic resistance conferred to them as a result of low permeability of their outer membrane (OM). Antibiotics capable of traversing the OM typically permeate through OM porins; thus, understanding the permeation properties of these porins is instrumental to the development of new antibiotics. A common macroscopic feature of many OM porins is their ability to transition between functionally distinct open and closed states that regulate transport properties and rate. To obtain a molecular basis for these processes, we performed tens of microseconds of molecular dynamics simulations of E. coli OM porin, OmpF. We observed that large-scale motion of the internal loop, L3, leads to widening and narrowing of the pore, suggesting its potential role in gating. Furthermore, Markov state analysis revealed multiple energetically stable conformations of L3 corresponding to open and closed states of the porin. Dynamics between these functional states occurs on the time scale of tens of microseconds and are mediated by the movement of highly conserved acidic residues of L3 to form H-bonds with opposing sides of the barrel wall of the pore. To validate our mechanism, we mutated key residues involved in the gating process that alter the H-bond pattern in the open/closed states and performed additional simulations. These mutations shifted the dynamic equilibrium of the pore towards open or closed states. Complementarily, the mutations favoring the open/closed states lead to increased/decreased accumulation of multiple antibiotics in our whole-cell accumulation assays. Notably, porins containing one of the mutations favoring the closed state has previously been found in antibiotic resistant bacterial strains. Overall, our 180 {micro}s of simulation data (wild type and mutants) with concerted experiments suggests that regulation of the dynamic equilibrium between open and closed states of OM porins could be a mechanism by which Gram-negative bacteria acquire antibiotic resistance.

category: biophysics


title: HIF-dependent expression of creatine kinase brain isoform (CKB) promotes breast cancer metastasis, and blocking creatine kinase activity with cyclocreatine impairs invasion and improves the efficacy of conventional chemotherapies

abstract: The oxygen-responsive Hypoxia Inducible Factor (HIF)-1 promotes several steps of the metastatic cascade. A hypoxic gene signature is enriched in triple negative breast cancers (TNBCs) and correlates with poor patient survival. Since inhibiting the HIF transcription factors with small molecules is challenging, we sought to identify genes downstream of HIF-1 that could be targeted to block invasion and metastasis. Creatine kinase brain isoform (CKB) was identified as a highly differentially expressed gene in a screen of HIF-1 wild type and knockout mammary tumor cells derived from a transgenic model of metastatic breast cancer. CKB is a cytosolic enzyme that reversibly catalyzes the phosphorylation of creatine, generating phosphocreatine (PCr) in the forward reaction, and regenerating ATP in the reverse reaction. Creatine kinase activity is inhibited by the creatine analog cyclocreatine (cCr). Loss and gain of function genetic approaches were used in combination with cCr therapy to define the contribution of CKB expression or creatine kinase activity to cell proliferation, migration, invasion, and metastasis in ER-negative breast cancers. CKB was necessary for cell invasion in vitro and strongly promoted tumor growth and metastasis in vivo. Similarly, cyclocreatine therapy repressed cell migration, cell invasion, formation of invadopodia, and lung metastasis. Moreover, in common TNBC cell line models, the addition of cCr to conventional cytotoxic chemotherapy agents was either additive or synergistic to repress tumor cell growth.

category: cancer biology


title: HIF-dependent expression of creatine kinase brain isoform (CKB) promotes breast cancer metastasis, whereas cyclocreatine therapy impairs invasion and improves the efficacy of conventional chemotherapies