title: Cation and anion channelrhodopsins: Sequence motifs and taxonomic distribution

abstract:
Cation and anion channelrhodopsins (CCRs and ACRs, respectively) primarily from two algal species, Chlamydomonas reinhardtii and Guillardia theta, have become widely used as optogenetic tools to control cell membrane potential with light. We mined algal and other protist polynucleotide sequencing projects and metagenomic samples to identify 75 channelrhodopsin homologs from three channelrhodopsin families, including one revealed in dinoflagellates in this study. We carried out electrophysiological analysis of 33 natural channelrhodopsin variants from different phylogenetic lineages and 10 metagenomic homologs in search of sequence determinants of ion selectivity, photocurrent desensitization, and spectral tuning in channelrhodopsins. Our results show that association of a reduced number of glutamates near the conductance path with anion selectivity depends on a wider protein context, because prasinophyte homologs with the identical glutamate pattern as in cryptophyte ACRs are cation-selective. Desensitization is also broadly context-dependent, as in one branch of stramenopile ACRs and their metagenomic homologs its extent roughly correlates with phylogenetic relationship of their sequences. Regarding spectral tuning, two prasinophyte CCRs exhibit red-shifted spectra to 585 nm, although their retinal-binding pockets do not match those of previously known similarly red-shifted channelrhodopsins. In cryptophyte ACRs we identified three specific residue positions in the retinal-binding pocket that define the wavelength of their spectral maxima. Lastly, we found that dinoflagellate rhodopsins with a TCP motif in the third transmembrane helix and a metagenomic homolog exhibit channel activity.

IMPORTANCEChannelrhodopsins are widely used in neuroscience and cardiology as research tools and are considered as prospective therapeutics, but their natural diversity and mechanisms remain poorly characterized. Genomic and metagenomic sequencing projects are producing an ever-increasing wealth of data, whereas biophysical characterization of the encoded proteins lags behind. In this study we used manual and automated patch clamp recording of representative members of four channelrhodopsin families including a family that we report in this study in dinoflagellates. Our results contribute to a better understanding of molecular determinants of ionic selectivity, photocurrent desensitization, and spectral tuning in channelrhodopsins.


title: Genetic analysis of coronary artery disease using tree-based automated machine learning informed by biology-based feature selection

abstract: Machine Learning (ML) approaches are increasingly being used in biomedical applications. Important challenges of ML include choosing the right algorithm and tuning the parameters for optimal performance. Automated ML (AutoML) methods, such as Tree-based Pipeline Optimization Tool (TPOT), have been developed to take some of the guesswork out of ML thus making this technology available to users from more diverse backgrounds. The goals of this study were to assess applicability of TPOT to genomics and to identify combinations of single nucleotide polymorphisms (SNPs) associated with coronary artery disease (CAD), with a focus on genes with high likelihood of being good CAD drug targets. We leveraged public functional genomic resources to group SNPs into biologically meaningful sets to be selected by TPOT. We applied this strategy to data from the UK Biobank, detecting a strikingly recurrent signal stemming from a group of 28 SNPs. Importance analysis of these uncovered functional relevance of the top SNPs to genes whose association with CAD is supported in the literature and other resources. Furthermore, we employed game-theory based metrics to study SNP contributions to individual level TPOT predictions and discover distinct clusters of well-predicted CAD cases. The latter indicates a promising approach towards precision medicine.

category: bioinformatics


title: PIMBA: a PIpeline for MetaBarcoding Analysis

abstract:
SummaryDNA metabarcoding is an emerging monitoring method capable of assessing biodiversity from environmental samples (eDNA). Advances in computational tools have been required due to the increase of Next-Generation Sequencing data. Tools for DNA metabarcoding analysis, such as MOTHUR, QIIME, Obitools, and mBRAVE have been widely used in ecological studies. However, some difficulties are encountered when there is a need to use custom databases. Here we present PIMBA, a PIpeline for MetaBarcoding Analysis, which allows the use of customized databases, as well as other reference databases used by the softwares mentioned here. PIMBA is an open-source and user-friendly pipeline that consolidates all analyses in just three command lines.

Availability and Implementationhttps://github.com/reinator/pimba

Contactrenato.renison@gmail.com and guilherme.oliveira@pq.itv.org.


title: Species richness determines C. difficile invasion outcome in synthetic human gut communities

abstract: Understanding the principles of colonization resistance of the gut microbiome to the pathogen Clostridioides difficile will enable the design of next generation defined bacterial therapeutics. We investigate the ecological principles of community resistance to C. difficile invasion using a diverse synthetic human gut microbiome. Our results show that species richness is a key determinant of C. difficile growth across a wide range of ecological contexts. Using a dynamic computational model, we demonstrate that C. difficile receives the largest number and magnitude of incoming negative interactions. We identify molecular mechanisms of inhibition including acidification of the environment and competition over glucose. We demonstrate that C. difficiles close relative Clostridium hiranonis strongly inhibits C. difficile via a pH-independent mechanism. While increasing the initial density of C. difficile can increase its abundance in the assembled community, the community context determines the maximum achievable C. difficile abundance. Our work suggests that the C. difficile inhibitory potential of defined bacterial therapeutics can be optimized by designing communities that feature a combination of mechanisms including species richness, environment acidification, and resource competition.

category: systems biology


title: Anterior claustrum cells are responsive during behavior but not passive sensory stimulation

abstract: The claustrum is uniquely positioned to communicate with almost all higher-order cortical areas through widespread and reciprocal anatomical projections, yet the in vivo functional properties of claustrum neurons are not well understood. Here we use microendoscope imaging in mice to measure activity in populations of genetically-labelled Gnb4+ claustrum neurons. We find that only a small fraction of cells in the anterior claustrum are responsive to visual or auditory stimuli when delivered under passive yet wakeful conditions. In contrast, during a visual behavioral task, the majority of cells in the anterior claustrum are strongly modulated, with separate and spatially intermingled cell populations showing either increases or decreases in activity relative to spontaneous levels. Our results suggest that the Gnb4+ cells in the anterior claustrum do not represent passively presented sensory stimuli; rather, these cells are strongly engaged during behavior associated with sensory-motor transformations.

category: neuroscience


title: Insular cortex corticotropin-releasing factor integrates stress signaling with social decision making.

abstract: Impairments in social cognition manifest in a variety of psychiatric disorders, making the neurobiological mechanisms underlying social decision making of particular translational importance. The insular cortex is consistently implicated in stress-related social and anxiety disorders, which are associated with diminished ability to make and use inferences about the emotions of others to guide behavior. We investigated how corticotropin releasing factor (CRF), a neuromodulator evoked by both self and social stressors, influenced the insula. In acute slices from male and female rats, CRF depolarized insular pyramidal neurons. In males, but not females, CRF suppressed presynaptic GABAergic inhibition leading to greater excitatory synaptic efficacy in a CRF receptor 1 (CRF1) and cannabinoid receptor 1 (CB1) dependent fashion. In males only, insular CRF increased social investigation, and CRF1 and CB1 antagonists interfered with social decision making. To investigate the molecular and cellular basis for the effect of CRF we examined insular CRF1 and CB1 mRNAs and found greater total insula CRF1 mRNA in females but greater CRF1 and CB1 mRNA colocalization in male insular cortex glutamatergic neurons which suggest complex, sex-specific organization of CRF and endocannabinoid systems. Together these results reveal a new sex-specific mechanism by which stress and affect contribute to social decision making.


title: Mycobacterium tuberculosis modifies cell wall carbohydrates during biofilm growth with a concomitant reduction in complement activation

abstract: There is an urgent need for drugs, new vaccines, and diagnostics for TB. It is recognised that research needed for the development of new vaccines for TB needs to be underpinned by understanding both the molecular and cellular mechanisms of host-pathogen interactions and how the immune response can be modulated to achieve protection with the use of a new vaccine for TB. Complement interacts with and orchestrates many aspects of the innate and adaptive immune responses and activation by Mycobacterium tuberculosis can be triggered by all three pathways. However, little is known about the contribution of each of these pathways during TB disease, particularly with respect to mycobacterial phenotype. There is strong evidence for extracellular communities of M. tuberculosis during TB disease (biofilms) that are found in the acellular rim of granulomas. These biofilms have been observed in cavities in lung resections from TB patients and are likely to be present in post-primary TB episodes in necrotic lesions. Our study aimed to understand more about the interactions between M. tuberculosis biofilms and complement activation, to determine which mycobacterial cell wall components are altered during biofilm growth, and how their alteration contributes to modulation of the complement response. We show that the lectin pathway has a reduced role compared to the classical pathway in initiating complement activation in biofilm bacteria. Analyses of the M. tuberculosis biofilm cell wall carbohydrate fractions revealed that there was reduced -glucan compared to planktonically-grown bacteria. Reduced C3b/iC3b deposition directly onto biofilm carbohydrates was observed which was consistent with both the observed reduction of C3b/iC3b deposition on biofilm bacilli and a reduction in the contribution of the lectin pathway in initiating complement activation on whole bacteria from biofilms, compared to planktonically-grown bacteria.

category: immunology


title: Innate triggering and antiviral effector functions of activin A

abstract:
First-line defence against viral infection is contingent upon rapid detection of conserved viral structural and genomic motifs by germline-encoded pattern recognition receptors, followed by activation of the type I IFN system and establishment of an intracellular antiviral state. Novel antiviral functions of bone morphogenetic protein and related activin cytokines, acting in conjunction with, and independently of, type I IFN, have recently been described. Activin A mediates multiple innate and adaptive immune functions - including antiviral effects. However, how such effects are mediated and how activin might be triggered by viral infection have not been defined. Here we addressed this in vivo and in vitro, in humans and mice.

Transcriptomic analyses delineated strikingly congruent patterns of gene regulation in hepatocytes stimulated with recombinant activin A and IFN in vitro. Activin A mRNA, encoded by INHBA, is induced upon activation of RIG-I, MDA5 and TLR7/8 viral nucleic acid sensors in vitro, across multiple cell lines and in human peripheral blood mononuclear cells. In vivo, infection of mice with influenza A also upregulated Inhba mRNA in the lung; this local upregulation of Inhba is retained in MAVS knockout mice, indicating a role for non-RIG-I-like receptors in its induction. Activin induction and signalling were also detectable in patients with chronic viral hepatitis.

Together, these data suggest Activin A is triggered in parallel with type I IFN responses and can trigger related antiviral effector functions. This model has implications for the development of targeted antiviral therapies, in addition to revealing novel facets of activin biology.


title: Capturing and analyzing pattern diversity: an example using the melanistic spotted patterns of leopard geckos

abstract: Animal color patterns are widely studied in ecology, evolution, and through mathematical modeling. Patterns may vary among distinct body parts such as the head, trunk or tail. As large amounts of photographic data is becoming more easily available, there is a growing need for general quantitative methods for capturing and analyzing the full complexity and details of pattern variation. Detailed information on variation in color pattern elements is necessary to understand how patterns are produced and established during development, and which evolutionary forces may constrain such a variation. Here, we develop an approach to capture and analyze variation in melanistic color pattern elements in leopard geckos. We use this data to study the variation among different body parts of leopard geckos and to draw inferences about their development. We compare patterns using 14 different indices such as the ratio of melanistic versus total area, the ellipticity of spots, and the size of spots and use these to define a composite distance between two patterns. Pattern presence/absence among the different body parts indicates a clear pathway of pattern establishment from the head to the back legs. Together with weak within-individual correlation between leg patterns and main body patterns, this suggests that pattern establishment in the head and tail may be independent from the rest of the body. We found that patterns vary greatest in size and density of the spots among body parts and individuals, but little in their average shapes. We also found a correlation between the melanistic patterns of the two front legs, as well as the two back legs, and also between the head, tail and trunk, especially for the density and size of the spots, but not their shape or inter-spot distance. Our data collection and analysis approach can be applied to other organisms to study variation in color patterns between body parts and to address questions on pattern formation and establishment in animals.

category: developmental biology


title: Comparing divergence landscapes from reduced-representation and whole-genome re-sequencing in the yellow-rumped warbler (Setophaga coronata) species complex

abstract: Researchers seeking to generate genomic data for non-model organisms are faced with a number of trade-offs when deciding which method to use. The selection of reduced representation approaches versus whole genome re-sequencing will ultimately affect the marker density, sequencing depth, and the number of individuals that can multiplexed. These factors can affect researchers ability to accurately characterize certain genomic features, such as landscapes of divergence--how FST varies across the genomes. To provide insight into the effect of sequencing method on the estimation of divergence landscapes, we applied an identical bioinformatic pipeline to three generations of sequencing data (GBS, ddRAD, and WGS) produced for the same system, the yellow-rumped warbler species complex. We compare divergence landscapes generated using each method for the myrtle warbler (Setophaga coronata coronata) and the Audubons warbler (S. c. auduboni), and for Audubons warblers with deeply divergent mtDNA resulting from mitochondrial introgression. We found that most high-FST peaks were not detected in the ddRAD dataset, and that while both GBS and WGS were able to identify the presence of large peaks, WGS was superior at a finer scale. Comparing Audubons warblers with divergent mitochondrial haplotypes, only WGS allowed us to identify small (10-20kb) regions of elevated differentiation, one of which contained the nuclear-encoded mitochondrial gene NDUFAF3. We calculated the cost per base pair for each method and found it was comparable between GBS and WGS, but significantly higher for ddRAD. These comparisons highlight the advantages of WGS over reduced representation methods when characterizing landscapes of divergence.

category: evolutionary biology


title: Genomes of Wolbachia endosymbionts from the human filarial parasites Mansonella perstans and Mansonella ozzardi

abstract:
The intracellular endosymbiotic proteobacteria Wolbachia have evolved across the phyla nematoda and arthropoda. In Wolbachia phylogeny, supergroup F is the only clade with members from both arthropod and filarial nematode hosts and therefore can provide unique insights into their evolution and biology. In this study, 4 new supergroup F Wolbachia genomes have been assembled using a metagenomic assembly and binning approach, wMoz and wMpe from the human filarial parasites Mansonella ozzardi and Mansonella perstans, and wOcae and wMoviF from the blue mason bee Osmia caerulescens and the sheep ked Melophagus ovinus respectively. A comprehensive phylogenomic analysis revealed two independent origins of filarial Wolbachia in supergroup F, likely from ancestral arthropod hosts, as well as a convergent pseudogenization and loss of the bacterioferritin gene across all filarial Wolbachia. The new genomes provide a valuable resource for further studies on symbiosis, evolution, and the discovery of new antibiotics to treat mansonellosis.

Data depositionRaw Illumina reads for the M. perstans isolate Mpe1 and the M. ozzardi isolate Moz2 have been deposited in the NCBI SRA database under BioProject accession numbers PRJNA666672 and PRJNA666671 respectively. The assembled genome and PGAP annotations for the Wolbachia wMpe and wMoz are available from the NCBI GenBank database with accession numbers JACZHU000000000 and JADAKK000000000 respectively. The annotated genomes for wOcae and wMoviF are available as GenBank accessions JAHDJT000000000 and JAHRBA000000000 respectively, deposited under BioProject accessions PRJNA244005 and PRJNA739984 respectively.

category: genomics


title: Multiple origins of nematode-Wolbachia symbiosis in supergroup F and convergent loss of bacterioferritin in filarial Wolbachia


title: High resolution ultrasonic neural modulation observedvia in vivo two-photon calcium imaging

abstract: Neural modulation plays a major role in delineating the circuit mechanisms and serves as the cornerstone of neural interface technologies. Among the various modulation mechanisms, ultrasound enables noninvasive label-free deep access to mammalian brain tissue. To date, most if not all ultrasonic neural modulation implementations are based on [~]1 MHz carrier frequency. The long acoustic wavelength results in a spatially coarse modulation zone, often spanning over multiple function regions. The modulation of one brain region is inevitably linked with the modulation of its neighboring regions. To significantly increase the spatial resolution, we explored the application of high-frequency ultrasound. To investigate the neuronal response at cellular resolutions, we developed a dual-modality system combining in vivo two-photon calcium imaging and focused ultrasound modulation. The studies show that the [~]30 MHz ultrasound can suppress the neuronal activity in awake mice at 100-micron scale spatial resolutions, paving the way for high-resolution ultrasonic neural modulation.

category: bioengineering


title: BubbleGun: Enumerating Bubbles and Superbubbles in Genome Graphs

abstract:
MotivationWith the fast development of third generation sequencing machines, de novo genome assembly is becoming a routine even for larger genomes. Graph-based representations of genomes arise both as part of the assembly process, but also in the context of pangenomes representing a population. In both cases, polymorphic loci lead to bubble structures in such graphs. Detecting bubbles is hence an important task when working with genomic variants in the context of genome graphs.

ResultsHere, we present a fast general-purpose tool, called BubbleGun, for detecting bubbles and superbubbles in genome graphs. Furthermore, BubbleGun detects and outputs runs of linearly connected bubbles and superbubbles, which we call bubble chains. We showcase its utility on de Bruijn graphs and compare our results to vgs snarl detection. We show that BubbleGun is considerably faster than vg especially in bigger graphs, where it reports all bubbles in less than 30 minutes on a human sample de Bruijn graph of around 2 million nodes.

AvailabilityBubbleGun is available and documented at https://github.com/fawaz-dabbaghieh/bubble_gun under MIT license.

Contactfawaz@hhu.de or tobias.marschall@hhu.de

Supplementary informationSupplementary data are available at Bioinformatics online.


title: Impact of increased membrane realism on conformational sampling of proteins

abstract: The realism and accuracy of lipid bilayer simulations through molecular dynamics (MD) is heavily dependent on the lipid composition. While the field is pushing towards implementing more heterogeneous and realistic membrane compositions, a lack of high-resolution lipidomic data prevents some membrane protein systems from being modeled with the highest level of realism. Given the additional diversity of real-world cellular membranes and protein-lipid interactions, it is still not fully understood how altering membrane complexity affects modeled membrane protein function or if it matters over long timescale simulations. This is especially true for organisms whose membrane environments have little to no computational study, such as the plant plasma membrane. Tackling these issues in tandem, a generalized, realistic, and asymmetric plant plasma with more than 10 different lipid species membrane is constructed herein. Classical MD simulations of pure membrane constructs were performed to evaluate how altering the compositional complexity of the membrane impacted the plant membrane properties. The apo form of a plant sugar transporter, OsSWEET2b, was inserted into membrane models where lipid diversity was calculated in either a size-dependent or -independent manner. An adaptive sampling simulation regime validated by Markov-state models was performed to capture the gating dynamics of OsSWEET2b in each of these membrane constructs. In comparison to previous OsSWEET2b simulations performed in a pure POPC bilayer, we confirm that simulations performed within a native-like membrane composition alter the stabilization of apo OsSWEET2b conformational states by ~1 kcal/mol. The free energy barriers of intermediate conformational states decrease when realistic membrane complexity is simplified, albeit roughly within sampling error, suggesting that protein-specific responses to membranes differ due to altered packing caused by compositional fluctuations. This work serves as a case study where a more realistic bilayer composition makes unbiased conformational sampling easier to achieve than with simplified bilayers.

category: biophysics


title: The SARS-CoV-2 replication-transcription complex is a priority target for broad-spectrum pan-coronavirus drugs

abstract: In the absence of effective treatment, COVID-19 is likely to remain a global disease burden. Compounding this threat is the near certainty that novel coronaviruses with pandemic potential will emerge in years to come. Pan-coronavirus drugs - agents active against both SARS-CoV-2 and other coronaviruses - would address both threats. A strategy to develop such broad-spectrum inhibitors is to pharmacologically target binding sites on SARS-CoV-2 proteins that are highly conserved in other known coronaviruses, the assumption being that any selective pressure to keep a site conserved across past viruses will apply to future ones. Here, we systematically mapped druggable binding pockets on the experimental structure of fifteen SARS-CoV-2 proteins and analyzed their variation across twenty-seven - and {beta}-coronaviruses and across thousands of SARS-CoV-2 samples from COVID-19 patients. We find that the two most conserved druggable sites are a pocket overlapping the RNA binding site of the helicase nsp13, and the catalytic site of the RNA-dependent RNA polymerase nsp12, both components of the viral replication-transcription complex. We present the data on a public web portal (https://www.thesgc.org/SARSCoV2_pocketome/) where users can interactively navigate individual protein structures and view the genetic variability of drug binding pockets in 3D.

category: pharmacology and toxicology


title: Genome-wide Association Study of Pancreatic Fat: The Multiethnic Cohort Adiposity Phenotype Study

abstract: Several studies have found associations between higher pancreatic fat content and adverse health outcomes, such as diabetes and the metabolic syndrome, but investigations into the genetic contributions to pancreatic fat are limited. This genome-wide association study, comprised of 804 participants with MRI-assessed pancreatic fat measurements, was conducted in the ethnically diverse Multiethnic Cohort-Adiposity Phenotype Study (MEC-APS). Two genetic variants reaching genome-wide significance, rs73449607 on chromosome 13q21.2 (Beta = -0.67, P = 4.50x10-8) and rs7996760 on chromosome 6q14 (Beta = -0.90, P = 4.91x10-8) were associated with percent pancreatic fat on the log scale. Rs73449607 was most common in the African American population (13%) and rs79967607 was most common in the European American population (6%). Rs73449607 was also suggestively associated with lower risk of type 2 diabetes (OR = 0.95, 95% CI = 0.89-1.00, P = 0.047) in the Population Architecture Genomics and Epidemiology (PAGE) Study and the DIAbetes Genetics Replication and Meta-analysis (DIAGRAM), which included substantial numbers of non-European ancestry participants (53,102 cases and 193,679 controls). Rs73449607 is located in an intergenic region between GSX1 and PLUT, and rs79967607 is in intron 1 of EPM2A. PLUT, a linkRNA, regulates transcription of an adjacent gene, PDX1, that controls beta-cell function in the mature pancreas, and EPM2A encodes the protein laforin, which plays a critical role in regulating glycogen production. If validated, these variants may suggest a genetic component for pancreatic fat and a common etiologic link between pancreatic fat and type 2 diabetes.

category: cell biology


title: Mast cell function in prostate inflammation, fibrosis, and smooth muscle cell dysfunction

abstract:
Intraurethral inoculation of mice with uropathogenic E. coli (CP1) results in prostate inflammation, fibrosis, and urinary dysfunction, recapitulating some but not all of the pathognomonic clinical features associated with benign prostatic hyperplasia (BPH) and lower urinary tract symptoms (LUTS). In both patients with LUTS and in CP1-infected mice, we observed increased numbers and activation of mast cells and elevated levels of prostate fibrosis. Therapeutic inhibition of mast cells and the histamine 1 receptor in the mouse model resulted in reduced mast cell activation in the prostate and significant alleviation of urinary dysfunction. Treated mice showed reduced prostate fibrosis, less infiltration of immune cells, and decreased inflammation. In addition, as opposed to symptomatic CP1-infected mice, treated mice showed reduced myosin light chain (MLC)-2 phosphorylation, a marker of prostate smooth muscle contraction. These results show that mast cells play a critical role in the pathophysiology of urinary dysfunction and may be an important therapeutic target for men with BPH/LUTS.

NEW AND NOTEWORTHYLUTS-associated BPH is derived from a combination of immune activation, extracellular matrix remodeling, hyperplasia, and smooth muscle cell contraction in prostates of men. Using a mouse model, we describe the importance of mast cells in regulating these multiple facets involved in the pathophysiology of LUTS. Mast cell inhibition alleviates both pathology and urinary dysfunction in this model suggesting the potential for mast cell inhibition as a therapeutic that prevents and reverses pathology and associated symptomology.


title: Avidity observed between a bivalent inhibitor and an enzyme monomer with a single active site

abstract:
Although myriad protein-protein interactions in nature use polyvalent binding, in which multiple ligands on one entity bind to multiple receptors on another, to date an affinity advantage of polyvalent binding has been demonstrated experimentally only in cases where the target receptor molecules are clustered prior to complex formation. Here, we demonstrate cooperativity in binding affinity (i.e., avidity) for a protein complex in which an engineered dimer of the amyloid precursor protein inhibitor (APPI), possessing two fully functional inhibitory loops, interacts with mesotrypsin, a soluble monomeric protein that does not self-associate or cluster spontaneously. We found that each inhibitory loop of the purified APPI homodimer was over three-fold more potent than the corresponding loop in the monovalent APPI inhibitor. This observation is consistent with a suggested mechanism whereby the two APPI loops in the homodimer simultaneously and reversibly bind two corresponding mesotrypsin monomers to mediate mesotrypsin dimerization. We propose a simple model for such dimerization that quantitatively explains the observed cooperativity in binding affinity. Binding cooperativity in this system reveals that the valency of ligands may affect avidity in protein-protein interactions including those of targets that are not surface-anchored and do not self-associate spontaneously. In this scenario, avidity may be explained by the enhanced concentration of ligand binding sites in proximity to the monomeric target, which may favor rebinding of the multiple ligand binding sites with the receptor molecules upon dissociation of the protein complex.

Impact statementLacham-Hartman et al. demonstrate enhancement of binding affinity through avidity in a complex between a bivalent ligand and a soluble monomeric target with a single binding site. Avidity effects have previously been demonstrated only for clustered receptor molecules presenting multiple binding sites. Our model may explain how polyvalent ligands can agonize or antagonize biological interactions involving nonclustered target molecules that are crucial for intra- and extracellular structural, metabolic, signaling, and regulatory pathways.

category: biochemistry


title: Mechanical influences on in silico tumor evolution

abstract:
Experimental insight and conceptual understanding of tumor growth are steadily growing and leading to new therapeutic interventions. Experiments and clinical studies are able to link single-cell properties to macroscopic tumor attributes. The development of cellular subpopulations in heterogeneous tumors can be understood as an evolutionary system with different cell types competing over both space and nutrients. However, to predict the growth trajectory and development of a tumor, fitness and trade-offs of cell properties in the context of the surroundings are required and often inaccessible. The optimum of the evolutionary trajectory provides a target for intervention, but can mostly not be identified. We propose that the optimal value of cellular properties is influenced by the tumor surrounding. Computational multiscale-modeling of tissue enables the observation of the trajectory of each cell while modeling the tumor surrounding. We model a 3D spheroid tumor and the fitness of individual cells and the evolutionary behavior of the tumor are quantified and linked to global parameters. Cell-cell adhesion and cell motility are two important mechanical properties for cell development and used as free parameters. Mechanical properties alone are able to drive the tumor towards low adhesion.We implement a dynamically changing nutrient surrounding representing the fluctuating blood-supply through blood vessel collapse and angiogenesis. We find that the evolutionary speed depends on the frequency of the fluctuations. We identify a frequency domain in which the evolutionary speed is significantly increased over a tumor with constant nutrient supply. The findings suggest that mechanically-induced fluctuations can accelerate tumor evolution.

Author summaryLimited space and nutrients together with competing cell types drive an evolutionary process inside tumors. This process selects for the fittest cell types and optimizes the growing behavior for its local surroundings. An expanding tumor exerts mechanical forces on its cells and its surroundings, leading to a fluctuating nutrient supply through collapsing blood vessels. Here, we observe the influence of a dynamically changing surrounding on the evolutionary behavior of heterogeneous tumors in a high-resolution computational model. We find that the evolutionary speed depends on the frequency of the fluctuations and a fitness advantage of low-adhesion cells.

category: cancer biology


title: Plant immune system activation is necessary for efficient interaction with auxin secreting beneficial bacteria

abstract: Plants continuously monitor the presence of microorganisms through their immune system to establish an adaptive response. Unlike immune recognition of pathogenic bacteria, how beneficial bacteria interact with the plant immune system is not well understood. Analysis of colonization of Arabidopsis thaliana by auxin producing beneficial bacteria revealed that activation of the plant immune system is necessary for efficient bacterial colonization and auxin secretion. A feedback loop is established in which bacterial colonization triggers an immune reaction and production of reactive oxygen species, which, in turn, stimulate auxin production by the bacteria. Auxin promotes bacterial survival and efficient root colonization, allowing the bacteria to inhibit fungal infection and promote plant health.

category: plant biology


title: Prostaglandin E2 induction by cytosolic Listeria monocytogenes in phagocytes is necessary for optimal T-cell priming

abstract:
Listeria monocytogenes is an intracellular bacterium that elicits robust CD8+ T-cell responses. Despite the ongoing development of L. monocytogenes-based platforms as cancer vaccines, our understanding of how L. monocytogenes drives robust CD8+ T-cell responses remains incomplete. One overarching hypothesis is that activation of cytosolic innate pathways is critical for immunity, as strains of L. monocytogenes that are unable to access the cytosol fail to elicit robust CD8+ T-cell responses and in fact inhibit optimal T-cell priming. Counterintuitively, however, activation of known cytosolic pathways, such as the inflammasome and type I IFN, lead to impaired immunity. Here, we describe a cytosol-dependent response that is critical for immunity to L. monocytogenes, namely production of prostaglandin E2 (PGE2) downstream of cyclooxygenase-2 (COX-2). Vacuole-constrained L. monocytogenes elicit reduced PGE2 production compared to wild-type strains in macrophages and dendritic cells ex vivo. In vivo, infection with wild-type L. monocytogenes leads to 10-fold increases in PGE2 production early during infection whereas vacuole-constrained strains fail to induce PGE2 over mock-immunized controls. Mice deficient in COX-2 specifically in Lyz2+ or CD11c+ cells produce less PGE2, suggesting these cell subsets contribute to PGE2 levels in vivo, while depletion of phagocytes with clodronate abolishes PGE2 production completely. Taken together, this work identifies the first known cytosol-dependent innate immune response critical for generating CD8+ T-cell responses to L. monocytogenes, suggesting that one reason cytosolic access is required to prime CD8+ T-cell responses may be due to induction of PGE2.

Author summaryL. monocytogenes is an intracellular bacterial pathogen that generates robust cell-mediated immune responses. Due to this robust induction, L. monocytogenes is used as both a model to understand how CD8+ T-cells are primed, as well as a platform for cancer immunotherapy vaccines. L. monocytogenes must enter the cytosol of an infected host cell to stimulate robust T-cell responses, however, which cytosolic innate pathway(s) contribute to T-cell priming remains unclear. Here, we define COX-2 dependent PGE2 production as the first cytosol-dependent innate immune response critical for immunity to L. monocytogenes. We found that ex vivo PGE2 production by macrophages and dendritic cells is partially dependent on cytosolic access, as vacuole-constrained strains of L. monocytogenes elicit reduced PGE2. In vivo, cytosolic access is essential for PGE2 production. L. monocytogenes elicits a 10-fold increase in PGE2 production, whereas strains of L. monocytogenes that cannot access the cytosol fail to elicit PGE2 compared to mock immunized mice. Furthermore, CD11c+ and Lyz2+ cells contribute to PGE2 production in vivo, as mice deficient in COX-2 in these cell subsets have impaired PGE2 production. Taken together, our work identifies the first known cytosol-dependent pathway that is critical for generating immunity to L. monocytogenes.


title: A deep-sea sulfate reducing bacterium directs the formation of zero-valent sulfur via sulfide oxidation

abstract: Zero-valent sulfur (ZVS) is a critical intermediate in the biogeochemical sulfur cycle. Up to date, sulfur oxidizing bacteria have been demonstrated to dominate the formation of ZVS. In contrast, formation of ZVS mediated by sulfate reducing bacteria (SRB) has been rarely reported. Here, we report for the first time that a typical sulfate reducing bacterium Desulfovibrio marinus CS1 directs the formation of ZVS via sulfide oxidation. In combination with proteomic analysis and protein activity assays, thiosulfate reductase (PhsA) and sulfide: quinone oxidoreductase (SQR) were demonstrated to play key roles in driving ZVS formation. In this process, PhsA catalyzed thiosulfate to form sulfide, which was then oxidized by SQR to form ZVS. Consistently, the expressions of PhsA and SQR were significantly up-regulated in strain CS1 when cultured in the deep-sea cold seep, strongly indicating strain CS1 might form ZVS in its real inhabiting niches. Notably, homologs of phsA and sqr widely distributed in the metagenomes of deep-sea SRB. Given the high abundance of SRB in cold seeps, it is reasonable to propose that SRB might greatly contribute to the formation of ZVS in the deep-sea environments. Our findings add a new aspect to the current understanding of the source of ZVS.

category: microbiology


title: Antimicrobial Resistance and Genomic Characterization of Salmonella Dublin Isolates in Cattle from the United States.

abstract: Salmonella enterica subspecies enterica serovar Dublin is a host-adapted serotype in cattle, associated with enteritis and systemic disease. While rare in humans, it can cause severe illness, including bacteremia, with hospitalization and death. In the United States, S. Dublin has become one of the most multidrug-resistant serotypes. The objective of this study was to characterize S. Dublin isolates from sick cattle by analyzing phenotypic and genotypic antimicrobial resistance (AMR) profiles, the presence of plasmids, and phylogenetic relationships. S. Dublin isolates (n=140) were selected from submissions to the NVSL for Salmonella serotyping (2014 - 2017) from 21 states. Isolates were tested for susceptibility against 14 class-representative antimicrobial drugs. Resistance profiles were determined using the ABRicate with Resfinder and NCBI databases, AMRFinder and PointFinder. Plasmids were detected using ABRicate with PlasmidFinder. Phylogeny was determined using vSNP. We found 98% of the isolates were resistant to more than 4 antimicrobials. Only 1 isolate was pan-susceptible and had no predicted AMR genes. All S. Dublin isolates were susceptible to azithromycin and meropenem. They showed 96% resistance to sulfonamides, 97% to tetracyclines, 95% to aminoglycosides and 85% to betalactams. The most common AMR genes were: sulf2 and tetA (98.6%), aph(3)-Ib and aph(6)-Id (96.4%), floR (94.3%), and blaCMY-2 (85.7%). All quinolone resistant isolates presented mutations in gyrA. Ten plasmid types were identified among all isolates with IncA/C2, IncX1, and IncFII(S) being the most frequent. The S. Dublin isolates show low genomic genetic diversity. This study provided antimicrobial susceptibility and genomic insight into S. Dublin clinical isolates from cattle in the U.S. Further sequence analysis integrating food and human origin S. Dublin isolates may provide valuable insight on increased virulence observed in humans.


title: Generation of functional ciliates cholangiocytes from human pluripotent stem cells

abstract: The derivation of mature functional cholangiocytes from human pluripotent stem cells (hPSCs) would provide a model for studying the pathogenesis of cholangiopathies and for developing novel therapies to treat them. Current differentiation protocols are not efficient and give rise to cholangiocytes that are not fully mature, limiting their therapeutic applications. Here, we describe a new strategy to generate functional hPSC-derived cholangiocytes that display many characteristics of mature bile duct cells including high levels of CFTR and the presence of primary cilia capable of sensing flow. With this level of maturation, these cholangiocytes are amenable for testing the efficacy of new cystic fibrosis drugs and for studying the role of cilia in cholangiocyte development and function. Transplantation studies showed that the mature cholangiocytes generate ductal structures in the liver of immunocompromised mice providing the first indication that it may be possible to develop cell-based therapies to restore bile duct function in patients with biliary disease.


title: Spatial heterogeneity enhance robustness of large multi-species ecosystems

abstract:
Understanding ecosystem stability and functioning is a long-standing goal in theoretical ecology, with one of the main tools being dynamical modelling of species abundances. With the help of dynamical population models limits to stability and regions of various ecosystem dynamics have been extensively mapped in terms of diversity (number of species), types of interactions, interaction strengths, varying interaction networks (for example plant-pollinator, food-web) and varying structures of these networks. Although it is apparent that ecosystems reside in and are affected by a spatial environment, local differences (spatial heterogeneity) is often excluded from studies mapping stability boundaries under the assumption of an average and equal amount of interaction for all individuals of a species. Here we show that extending the classic dynamical Generalised-Lotka-Volterra model into a connected space the boundaries of stability change. When viewing the ecosystem as a spatially heterogeneous whole, limits previously marking the end of stability can now be crossed without any remarkable change in species abundances and without loss of stability. Thus limits previously thought to mark catastrophic transitions are not critical due to the possibility of spatial heterogeneity within the system. In addition, we show that too much spatial fragmentation of ecosystem habitats acts destabilising and leads back to the stability boundaries found in spatially homogeneous ecosystems with average interactions. Thus, we conclude that spatially heterogeneous but connected systems are the most robust. In terms of ecosystem management, the risk of collapse or irreversible changes is lower in spatially heterogeneous systems, which real ecosystems are, and we should expect local changes in populations well in advance of system collapse. Although, too much fragmentation of an ecosystems available space can lead to a less robust system with higher risk of extinctions and collapse.

Author summaryOne of the major challenges facing humanity is the fragmentation of wildlife habitats and decline in biodiversity due to human need for resources and land-use practices. We need to find ways to combine human prosperity with biodiversity conservation. To achieve this a solid understanding of ecosystem stability and functioning is paramount. One way to gain such insight is to find limits when we expect species to go extinct or ecosystems to collapse by simulations of interacting species populations. Many such stability limits have been found theoretically the last decades, but for simplification of modelling, studies often exclude that ecosystems are spread out in space. By doing so, studies assume an average and equal amount of interaction for all individuals of a species. Here, we explicitly include space and thus allow for migration and spatial heterogeneity (local differences) in interactions. When modelling the ecosystem as a spatially heterogeneous whole, limits previously marking extinction or collapse can now be crossed without any remarkable change in species abundances. Thus, natural hindrances to migration improve ecosystem robustness and limits previously thought to mark catastrophic transitions are not critical due to the possibility of spatial heterogeneity within the system. In addition, we reconfirm that a large amount of fragmentation acts destabilising for an ecosystem. Thus, we conclude that the most robust ecosystems are spatially heterogeneous but connected.

category: ecology


title: Molecular phylogeny of Puerto Rico Bank dwarf geckos (Squamata: Sphaerodactylidae: Sphaerodactylus)

abstract: AO_SCPLOWBSTRACTC_SCPLOWThe genus Sphaerodactylus is a very species-rich assemblage of sphaerodactylid lizards that has undergone a level of speciation in parallel to that of the well-known Anolis lizards. Nevertheless, molecular phylogenetic research on this group consists of a handful of smaller studies of regional focus (e.g., western Puerto Rico, the Lesser Antilles) or large-scale analyses based on relatively limited sequence data. Few medium-scale multi-locus studies exist-- for example, studies that encompass an entire radiation on an island group. Building upon previous work done in Puerto Rican Sphaerodactylus, we performed multi-locus sampling of Sphaerodactylus geckos from across the Puerto Rico Bank. We then used these data for phylogeny estimation with near-complete taxon sampling. We focused on sampling the widespread nominal species S. macrolepis and in so doing, we uncovered a highly divergent and morphologically distinct lineage of Sphaerodactylus macrolepis from Puerto Rico, Culebra, and Vieques islands, which we recognize as S. grandisquamis (Stejneger, 1904) on the basis of molecular and morphological characters. S. grandisquamis co-occurs with S. macrolepis only on Culebra Island but is highly genetically differentiated and morphologically distinct. Sphaerodactylus macrolepis is now restricted to the eastern Puerto Rico Bank, from Culebra east through the Virgin Islands and including the topotypic population on St. Croix. We include additional discussion of the evolutionary history and historical biogeography of the Sphaerodactylus of the Puerto Rican Bank in the context of these new discoveries.


title: PRSS2 promotes tumor growth and progression by repressing Tsp-1 in the tumor microenvironment via binding to LRP1

abstract: In the earliest stages of tumor development, epithelial tumors (carcinomas) are physically confined to the area of the tissue in which they form. These nascent lesions (carcinomas in situ) are sequestered from the tissue parenchyma by the basement membrane. Within the tissue parenchyma lie a myriad of cell types comprised of fibroblasts, immune and inflammatory cells and endothelial cells. Upon invasion across the basement membrane and into the tissue parenchyma, tumors must manipulate the expression of pro- and anti-tumorigenic proteins such that pro-tumorigenic factors are produced in excess to anti-tumorigenic proteins. One such anti-tumorigenic protein is Thrombospondin-1 (Tsp-1). We have previously demonstrated that stimulation of Tsp-1 in the tumor microenvironment (TME) potently inhibits tumor growth and progression and in some cases induces tumor regression. Here, we identify a novel tumor-mediated mechanism to repress the expression of Tsp-1 in the TME via secretion of the serine protease PRSS2. We demonstrate that PRSS2 represses Tsp-1, not via its enzymatic activity, but by binding to low-density lipoprotein receptor-related protein 1 (LRP1). These findings describe a novel activity for PRSS2 through binding to LRP1 and represent a potential therapeutic strategy to treat cancer by blocking the PRSS2-mediated repression of Tsp-1.


title: PRSS2 stimulates tumor growth by remodeling the TME via repression of Tsp1


title: The proteome of remyelination is different from that of developmental myelination

abstract: Loss of myelin underlies the pathology of several neurological disorders of diverse etiology. CNS remyelination by adult oligodendrocyte progenitor cells (OPCs) can occur but it differs from developmental myelination carried out by neonatal OPCs. We asked whether the myelin proteome of remyelinated regions is changed. We compared the myelin proteome formed during development to the remyelination proteome attained after lysolecithin-induced demyelination in the mouse spinal cord. Mass-spectrometry analysis of iTRAQ labelled myelin protein lysates showed that the proteome of remyelination is different from that of developmental myelination, leading to profound changes in myelin protein content. Aside from known mediators of oligodendrocyte differentiation, we found proteome alterations included modulators of metabolism, cell signaling and actin cytoskeleton dynamics. Downregulating one candidate (FSCN1/Fascin1) was sufficient to partially hamper oligodendrocytes in-vitro. In summary, we identify the difference in the proteome of remyelinating oligodendrocytes as a novel potential contributor to the pathophysiology of demyelinating disorders, thus providing new potential therapeutic targets for future studies.


title: When makes you unique: temporality of the human brain fingerprint

abstract: The extraction of "fingerprints" from human brain connectivity data has become a new frontier in neuroscience. However, the time scales of human brain identifiability have not been addressed yet. In other words, what temporal features make our brains more "identifiable"? We here explore the dynamics of brain fingerprints (or brainprints) along two complementary axes: 1) what is the optimal time scale at which brainprints integrate sufficient information, 2) when best identification happens. Using dynamic identifiability, we show that the best identification emerges at longer time scales (~300s); however, short transient "bursts of identifiability" persist even when looking at shorter functional interactions. We find that these bursts of identifiability might be strongly associated with neuronal activity. Furthermore, we report evidence that different parts of connectome fingerprints relate to different time scales: i.e., more visual-somatomotor at short temporal windows, more frontoparietal-DMN driven by increasing temporal windows. Finally, using a meta-analytic approach, we show that there is a broad spectrum of associations between brainprints and behavior. At faster time scales, human brain fingerprints are linked to multisensory stimulation, eye movements, affective processing, visuospatial attention. At slower time scales instead, we find higher-cognitive functions, such as language and verbal semantics, awareness, declarative and working memory, social cognition. We hope that this first investigation of the temporality of the human brain fingerprint will pave the way towards a better understanding of what and when makes our brains unique.


title: Neural markers of speech comprehension: measuring EEG tracking of linguistic speech representations, controlling the speech acoustics

abstract:
When listening to speech, our brain responses time-lock to acoustic events in the stimulus. Recent studies have also reported that cortical responses track linguistic representations of speech. However, tracking of these representations is often described without controlling for acoustic properties. Therefore, the response to these linguistic representations might reflect unaccounted acoustic processing rather than language processing. Here, we evaluated the potential of several recently proposed linguistic representations as neural markers of speech comprehension. To do so, we investigated EEG responses to audiobook speech of 29 participants (22 [female]). We examined whether these representations contribute unique information over and beyond acoustic neural tracking and each other. Indeed, not all of these linguistic representations were significantly tracked after controlling for acoustic properties. However, phoneme surprisal, cohort entropy, word surprisal, and word frequency were all significantly tracked over and beyond acoustic properties. We also tested the generality of the associated responses by training on one story and testing on another. In general, the linguistic representations are tracked similarly across different stories spoken by different readers. These results suggests that these representations characterize processing of the linguistic content of speech.

Significance StatementFor clinical applications it would be desirable to develop a neural marker of speech comprehension derived from neural responses to continuous speech. Such a measure would allow for behaviour-free evaluation of speech understanding; this would open doors towards better quantification of speech understanding in populations from whom obtaining behavioral measures may be difficult, such as young children or people with cognitive impairments, to allow better targeted interventions and better fitting of hearing devices.


title: Ezh2 mediates epigenetic regulation of osteoclastogenesis and bone remodeling in mice

abstract: Osteoclasts derived from hematopoietic stem cells control bone resorption. Identifying novel molecules that can epigenetically regulate osteoclastogenesis has been an important basic and clinical issue. The polycomb group (PcG) protein enhancer of zeste homologue 2 (Ezh2), a histone lysine methyltransferase is associated with epigenetic regulation of numerous cellular processes, it is not yet clear on its involvement in bone cell development and homeostasis. Here, we crossed LysM-Cre mice with Ezh2flox/flox mice to delete Ezh2 in myeloid cell lineage mature macrophages. Conditional deletion of Ezh2 in macrophages resulted in significant increases in postnatal bone growth in the first 6 months of life, but tibia length and body weight gains were not different in knockout mice compared with their wild type controls. Significantly decreased osteoclastogenesis but increased bone mass without osteopetrosis were found in Ezh2 conditional knockout (CKO) mice. In contrast to female mice, one floxed Ezh2 gene copy recombinant with LysM-Cre+ (Ezh2flox/+LysM-Cre+) produced increased bone mass in young adult male mice compared with control mice (Ezh2flox/flox, LysM-Cre+ and wild type). Inflammatory milieu in bone was significantly lower in both male and female CKO mice compared with their respective controls. Deletion of Ezh2 in macrophages triggered increased gene expression of osteoclast suppressors, IRF8, MafB and Arg1 due to decreased Ezh2-induced trimethylation of H3K27me3. Conversely, NFATc1 and Cathepsin k expression were decreased. These findings suggest that pre-osteoclastic cell differentiation is under epigenetic control of osteoclast suppressive gene expression via an Ezh2-dependent mechanisms.


title: Effects of domestication on neophobia: A comparison between the domesticated Bengalese finch (Lonchura striata var. domestica) and its wild ancestor, the white-rumped munia (Lonchura striata)

abstract: Bengalese finches (Lonchura striata var. domestica) have more complex song traits than their wild ancestors, white-rumped munias (Lonchura striata). Domesticated finches are likely able to allocate more resources to song development rather than allocating resources to mechanisms aimed at coping with predation, which are no longer needed under domesticated conditions. Here, we aimed to examine the effects of changes in selection pressure due to domestication on the behaviour of Bengalese finches and to contemplate the possible evolutionary mechanisms underlying these changes. To do so, we compared neophobic responses to novel-object conditions as an assessment of reactions to potential predators. We studied groups of Bengalese finches and white-rumped munias and found that Bengalese finches were more likely to eat the food provided to them under novel-object conditions. Bengalese finches had a shorter latency time to eat, and this latency time was less affected by the novel object in the case of Bengalese finches compared to white-rumped munias. Therefore, Bengalese finches have reduced neophobic responses due to domestication. The behavioural strategies of white-rumped munias appear to be more suitable for natural environments, which include unpredictable risks, whereas Bengalese finches have likely adapted their behaviour to the conditions of artificial selection.

category: animal behavior and cognition


title: Phospholipidosis is a shared mechanism underlying the in vitro antiviral activity of many repurposed drugs against SARS-CoV-2

abstract:
Repurposing drugs as treatments for COVID-19 has drawn much attention. A common strategy has been to screen for established drugs, typically developed for other indications, that are antiviral in cells or organisms. Intriguingly, most of the drugs that have emerged from these campaigns, though diverse in structure, share a common physical property: cationic amphiphilicity. Provoked by the similarity of these repurposed drugs to those inducing phospholipidosis, a well-known drug side effect, we investigated phospholipidosis as a mechanism for antiviral activity. We tested 23 cationic amphiphilic drugs--including those from phenotypic screens and others that we ourselves had found--for induction of phospholipidosis in cell culture. We found that most of the repurposed drugs, which included hydroxychloroquine, azithromycin, amiodarone, and four others that have already progressed to clinical trials, induced phospholipidosis in the same concentration range as their antiviral activity; indeed, there was a strong monotonic correlation between antiviral efficacy and the magnitude of the phospholipidosis. Conversely, drugs active against the same targets that did not induce phospholipidosis were not antiviral. Phospholipidosis depends on the gross physical properties of drugs, and does not reflect specific target-based activities, rather it may be considered a confound in early drug discovery. Understanding its role in infection, and detecting its effects rapidly, will allow the community to better distinguish between drugs and lead compounds that more directly impact COVID-19 from the large proportion of molecules that manifest this confounding effect, saving much time, effort and cost.

One Sentence SummaryDrug-induced phospholipidosis is a single mechanism that may explain the in vitro efficacy of a wide-variety of therapeutics repurposed for COVID-19.


title: B.1.526 SARS-CoV-2 variants identified in New York City are neutralized by vaccine-elicited and therapeutic monoclonal antibodies

abstract: DNA sequence analysis recently identified the novel SARS-CoV-2 variant B.1.526 that is spreading at an alarming rate in the New York City area. Two versions of the variant were identified, both with the prevalent D614G mutation in the spike protein together with four novel point mutations and with an E484K or S477N mutation in the receptor binding domain, raising concerns of possible resistance to vaccine-elicited and therapeutic antibodies. We report that convalescent sera and vaccine-elicited antibodies retain full neutralizing titer against the S477N B.1.526 variant and neutralize the E484K version with a modest 3.5-fold decrease in titer as compared to D614G. The E484K version was neutralized with a 12-fold decrease in titer by the REGN10933 monoclonal antibody but the combination cocktail with REGN10987 was fully active. The findings suggest that current vaccines and therapeutic monoclonal antibodies will remain protective against the B.1.526 variants. The findings further support the value of wide-spread vaccination.


title: Isolation and Characterization of Cross-Neutralizing Coronavirus Antibodies from COVID-19+ Subjects

abstract: SARS-CoV-2 is one of three coronaviruses that have crossed the animal-to-human barrier in the past two decades. The development of a universal human coronavirus vaccine could prevent future pandemics. We characterized 198 antibodies isolated from four COVID19+ subjects and identified 14 SARS-CoV-2 neutralizing antibodies. One targeted the NTD, one recognized an epitope in S2 and twelve bound the RBD. Three anti-RBD neutralizing antibodies cross-neutralized SARS-CoV-1 by effectively blocking binding of both the SARS-CoV-1 and SARS-CoV-2 RBDs to the ACE2 receptor. Using the K18-hACE transgenic mouse model, we demonstrate that the neutralization potency rather than the antibody epitope specificity regulates the in vivo protective potential of anti-SARS-CoV-2 antibodies. The anti-S2 antibody also neutralized SARS-CoV-1 and all four cross-neutralizing antibodies neutralized the B.1.351 mutant strain. Thus, our study reveals that epitopes in S2 can serve as blueprints for the design of immunogens capable of eliciting cross-neutralizing coronavirus antibodies.


title: 3D visualization of SARS-CoV-2 infection and receptor distribution in Syrian hamster lung lobes display distinct spatial arrangements

abstract:
SARS-CoV-2 attaches to angiotensin-converting enzyme 2 (ACE2) to gain entry into cells after which the spike protein is cleaved by the transmembrane serine protease 2 (TMPRRS2) to facilitate viral-host membrane fusion. ACE2 and TMPRRS2 expression profiles have been analyzed at the genomic, transcriptomic, and single-cell RNAseq level, however, biologically relevant protein receptor organization in whole tissues is still poorly understood. To describe the organ-level architecture of receptor expression, related to the ability of ACE2 and TMPRRS2 to mediate infectivity, we performed a volumetric analysis of whole Syrian hamster lung lobes. Lung tissue of infected and control animals were stained using antibodies against ACE2 and TMPRRS2, combined with fluorescent spike protein and SARS-CoV-2 nucleoprotein staining. This was followed by light-sheet microscopy imaging to visualize expression patterns. The data demonstrates that infection is restricted to sites with both ACE2 and TMPRRS2, the latter is expressed in the primary and secondary bronchi whereas ACE2 is predominantly observed in the terminal bronchioles and alveoli. Conversely, infection completely overlaps at these sites where ACE2 and TMPRSS2 co-localize.

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=154 SRC="FIGDIR/small/435771v1_ufig1.gif" ALT="Figure 1">
View larger version (60K):
org.highwire.dtl.DTLVardef@150fc65org.highwire.dtl.DTLVardef@1ea6a56org.highwire.dtl.DTLVardef@eb3c43org.highwire.dtl.DTLVardef@1c1877e_HPS_FORMAT_FIGEXP  M_FIG C_FIG


title: Untangling the cell immune response dynamic for severe and critical cases of SARS-CoV-2 infection

abstract: COVID-19 is a global pandemic leading high death tolls worldwide day by day. Clinical evidence suggests that COVID-19 patients can be classified as non-severe, severe and critical cases. In particular, studies have highlighted the relationship between the lymphopenia and the severity of the illness, where CD8+ T cells have the lowest levels in critical cases. In this work, we aim to elucidate the key parameters that define the course of the disease deviating from severe to critical case. To this end, several mathematical models are proposed to represent the dynamic of the immune response in patients with SARS-CoV-2 infection. The best model had a good fit to reported experimental data, and in accordance with values found in the literature. Our results suggest that a rapid proliferation of CD8+ T cells is decisive in the severity of the disease.


title: Translocation of chloroplast NPR1 to the nucleus in retrograde signaling for adaptive response to salt stress in tobacco

abstract: Chloroplasts play a pivotal role in biotic and abiotic stress responses, accompanying changes in the cell reduction/oxidation (redox) state. Chloroplasts are an endosymbiotic organelle that sends retrograde signals to the nucleus to integrate with environmental changes. This study showed that salt stress causes the rapid accumulation of the nonexpressor of pathogenesis-related genes 1 (NPR1) protein, a redox-sensitive transcription coactivator that elicits many tolerance responses in chloroplasts and the nucleus. The transiently accumulated chloroplast NPR1 protein was translocated to the nucleus in a redox-dependent manner under salinity stress. In addition, immunoblotting and fluorescence image analysis showed that chloroplast-targeted NPR1-GFP fused with cTP (chloroplast transit peptide from RbcS) was localized in the nucleus during the responses to salt stress. Chloroplast functionality was essential for retrograde translocation, in which the stomules and cytoplasmic vesicles participated. Treatments with H2O2 and an ethylene precursor enhanced this retrograde translocation. Compared to each wild-type plant, retrograde signaling-related gene expression was severely impaired in the npr1-1 mutant in Arabidopsis, but enhanced transiently in the NPR1-Ox transgenic tobacco line. Therefore, NPR1 might be a retrograde signaling hub that improves a plants adaptability to changing environments.


title: miRador: a fast and precise tool for the prediction of plant miRNAs

abstract: Plant microRNAs (miRNAs) are short, non-coding RNA molecules that restrict gene expression via post-transcriptional regulation and function in several essential pathways including development, growth, and stress responses. Accurately identifying miRNAs in populations of small RNA (sRNA) sequencing libraries is a computationally intensive process which has resulted in the misidentification of inaccurately annotated miRNA sequences. In recent years, criteria for miRNA annotation have been refined to reduce these misannotations. Here, we describe miRador, a novel miRNA identification tool that utilizes the most up-to-date, community-established criteria for accurate identification of miRNAs in plants. We combine target prediction and Parallel Analysis of RNA Ends (PARE) data to assess the precision of the miRNAs identified by miRador. We compare miRador to other commonly used miRNA prediction tools and we find that miRador is at least as precise as other prediction tools while being significantly faster than other tools.


title: Deep Learning Predicts Patterns of Cardiotoxicity in a High-Content Screen Using Induced Pluripotent Stem Cell-Derived Cardiomyocytes

abstract:
Drug-induced cardiotoxicity and hepatotoxicity are major causes of drug attrition. To decrease late-stage drug attrition, pharmaceutical and biotechnology industries need to establish biologically relevant models that use phenotypic screening to predict drug-induced toxicity. In this study, we sought to rapidly detect patterns of cardiotoxicity using high-content image analysis with deep learning and induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs). We screened a library of 1280 bioactive compounds and identified those predicted to have cardiotoxic liabilities using a single-parameter score based on deep learning. Compounds with major predicted cardiotoxicity included DNA intercalators, ion channel blockers, epidermal growth factor receptor, cyclin-dependent kinase, and multi-kinase inhibitors. We also screened a diverse library of molecules with unknown targets and identified chemical frameworks with predicted cardiotoxic liabilities. By using this screening approach during target discovery and lead optimization, we can de-risk early-stage drug discovery. We show that the broad applicability of combining deep learning with iPSC technology is an effective way to interrogate cellular phenotypes and identify drugs that protect against diseased phenotypes and deleterious mutations.

GRAPHICAL ABSTRACT

O_FIG O_LINKSMALLFIG WIDTH=161 HEIGHT=200 SRC="FIGDIR/small/436666v1_ufig1.gif" ALT="Figure 1">
View larger version (42K):
org.highwire.dtl.DTLVardef@695707org.highwire.dtl.DTLVardef@1d0980forg.highwire.dtl.DTLVardef@1af30e7org.highwire.dtl.DTLVardef@b3a312_HPS_FORMAT_FIGEXP  M_FIG C_FIG CONTRIBUTION TO THE FIELDIn this article, Grafton and colleagues use induced pluripotent stem cell technology and deep learning to train a neural network capable of detecting patterns of cardiotoxicity. To identify bioactive and chemical classes that lead to cardiotoxicity, they combine the neural network with high-content screening of 2560 compounds. The methods described in this study can be used to de-risk early-stage drug development, triage hits, and identify drugs that protect against disease. This screening paradigm will serve as a useful resource for drug discovery and phenotypic interrogation of stem cells and stem cell-derived cell types.


title: Tick hazard in the South Downs National Park (UK): species, distribution, key locations for future interventions, site density, habitats

abstract:
BackgroundThe South Downs National Park (SDNP) is the UKs most visited National Park, and a foci of tick-borne Lyme disease. A range of human pathogens have been detected in UK ticks and related hosts, and the first presumed autochthonous cases of tick-borne encephalitis and babesiosis were recorded in 2019-20. SDNPs key objectives include conserving wildlife and encouraging enjoyment of the countryside, so interventions are needed that reduce hazard without negatively affecting ecosystem health. To be successful these require knowledge of site hazards, and we aimed to provide this to enable action.

MethodsBritish Deer Society volunteers submitted ticks removed from deer. Key potential intervention sites were selected and ticks collected by drag-sampling six 50 m2 transects per site, in most cases twice yearly for two years. Ticks were identified in-lab (sex, life stage, species), hazard measured as tick presence, Density of Ticks (all life stages, DOT), and Density of Nymphs (DON). Sites and habitat types were analysed for association with hazard. Distribution across SDNP was mapped in a Geographic Information System (GIS), by combining and comparing our fieldwork results with records from five other data sources (recent and historic).

Results87 Ixodes ricinus (all but one adults, 82%F) were removed from 14 deer (Dama dama n=10; Capreolus capreolus n=3; 1 not recorded; tick burden, 1-35) at 12 locations (commonly woodland). Five potential key intervention sites were identified and drag-sampled 2015-16, collecting 623 ticks (238 on-transects): 53.8% nymphs, 42.5% larvae, 3.7% adults (13M, 10F). Ticks were present on-transects at all sites drag-sampled (I.ricinus at three, Haemaphysalis punctata at two). The Mens (TM, the quietist site for human visitors) had the highest DOT at 30/300 m2 (DON=30/300 m2), followed by Queen Elizabeth Country Park (QECP, the busiest) at 22/300 m2 (12/300 m2), Cowdray Estate (CE) at 8/300 m2 (6/300 m2), and Seven Sisters Country Park (SSCP) at 1/300 m2 (1/300 m2). Ditchling Beacon Nature Reserve (DBNR) was sampled 2016 only (one adult H.punctata collected). Woodland had significantly higher hazard than grazed downland, but ticks were present at all downland sites drag-sampled. GIS mapping showed I.ricinus identified in 33/37 of SDNPs 10 km2 grid squares, Ixodes hexagonus 10/37, H.punctata 7/37, Dermacentor reticulatus 1/37.

ConclusionsMapping shows tick hazard is broadly distributed across SDNP. Ixodes ricinuswas most common, though the seeming range expansion of H.punctata is concerning, particularly as it seems to thrive better on grazed downland than I.ricinus. Site specific recommendations include: management of small high hazard plots with heavy visitor numbers (QECP); signage on post-visit precautions (all sites); repellent impregnated clothing for deerstalkers (CE); flock trials to control H.punctata (SSCP, DBNR). Further research at TM, which has high tick density, may contribute to knowledge on ecological dynamics underlying infection density, and the potential use of predator re-introduction/protection as a public health intervention. Ecological research on H.punctata would aid control. The SDNP Authority is ideally placed to link and champion site-based and regional policies to reduce hazard, whilst avoiding or reducing conflict between public health and ecosystem health.


title: Quantitative single cell 5hmC sequencing reveals non-canonical gene regulation by non-CG hydroxymethylation

abstract: Oxidative modification of 5-methylcytosine (5mC) generates 5-hydroxymethylcytosine (5hmC), a DNA modification that exhibits unique epigenetic regulatory functions and impacts diverse biological processes. However, standard single-nucleus/cell bisulfite sequencing methods cannot resolve the base ambiguity between 5mC and 5hmC to accurately measure cell-type specific epigenomic patterns and gene regulatory functions of 5hmC or true 5mC. Here, we develop single-nucleus 5hmC sequencing (snhmC-seq) for quantitative and unbiased 5hmC profiling in single cells by harnessing differential deaminase activity of APOBEC3A towards 5mC and chemically protected 5hmC. We used snhmC-seq to profile single nuclei from cryopreserved mouse brain samples to reveal epigenetic heterogeneity of 5hmC at single-cell resolution and uncovered a non-canonical gene regulatory role of genic 5hmC in non-CG context.


title: Integrating transcription factor abundance with chromatin accessibility in human erythroid lineage commitment

abstract: We present InTAC-seq, a method for simultaneous quantification of genome-wide open chromatin and intracellular protein abundance in fixed cells. Using InTAC we directly observe variation in chromatin accessibility and transcription factor motif occupancy driven by differences in transcription factor protein abundance. By purifying bone marrow progenitor cells based on GATA1 protein expression, we establish its role in both functional and epigenetic restriction of erythroid cell identity in human hematopoiesis.


title: A comparison of rotation-based scores for gene set analysis

abstract:
Gene-wise differential expression is usually the first major step in the statistical analysis of high-throughput data obtained from tech-niques such as microarrays or RNA-sequencing. The analysis at gene level is often complemented by the screening of the data in a broader biological context that considers as unit of analysis meaningful groups of genes that may have functions in certain biological processes.

Among the vast number of publications about gene set analysis, the rotation test for gene set analysis, also referred by roast, is a general sample randomization approach that maintains the integrity of the intra-gene set correlation structure in defining the null distribution of the test.

In this work we compare the performance of several enrichment score functions using such rotational approach for hypothesis testing. We find that computationally intensive measures based on Kolmogorov-Smirnov statistics fail to improve the rates of simpler measures of GSA like mean and maxmean scores. We also show the importance of ac-counting for the gene linear dependence structure of the testing set, which it is linked to the loss of effective signature size. In this regard, weighted statistics are introduced with the aim of maximizing the ef-fective signature size. These are found to out-power other usual scores in some simulations scenarios. The average of absolute values is found to be the most powerful score using both simulated and benchmarking data.

All tools are available in the roastgsa R package.

category: bioinformatics


title: RNA-Ribo Explorer: interactive mining and visualisation of Ribosome profiling data

abstract:
RNA translation has long been thought as a stable and uniform process by which a ribosome produces a protein encoded by the main Open Reading Frame (ORF) of an mRNA. Recently, growing evidence support incomplete correlation between RNA and protein abundance levels, the existence of alternative ORFs in numerous mammalian RNAs, and the involvement of ribosomes in gene expression regulation, thereby challenging previous views of translation. Ribosome profiling (aka Ribo-seq) has renewed the study of translation by enabling the mapping of translating ribosomes on the whole transcriptome using deep-sequencing.

Despite increasing use of Ribo-seq, recent review articles conclude that flexible, interactive tools for mining such data are missing. As Ribo-seq protocols still evolve, flexibility is highly desirable for the end-user. Here we describe RNA-Ribo-Explorer (RRE) a stand-alone tool that fills this gap. With RRE, one can explore read-count profiles of RNAs obtained after mapping, compare them between conditions, and visualize the profiles of individual RNAs. Importantly, the user can mine the data by defining queries that combine several criteria to detect interesting subsets of RNAs. For instance, one can ask RRE to find all RNAs whose translation of UTR region compared to that of the main ORF has changed between two conditions. This feature seems useful for finding candidate RNAs whose translation status or processing has changed across conditions.

RRE is a platform independent software and is freely available at https://gite.lirmm.fr/rivals/RRE/-/releases.


title: The synergic role of actomyosin architecture and biased detachment in muscle energetics: insights in cross bridge mechanism beyond the lever-arm swing

abstract: Muscle energetics reflects the ability of myosin motors to convert chemical energy into mechanical energy. How this process takes place remains one of the most elusive questions in the field. Here we combined experimental measurements of in vitro sliding velocity based on DNA-origami built filaments carrying myosins with different lever arm length and simulations based on a Monte-Carlo model which accounts for three basic components: (i) the geometrical hindrance, (ii) the mechano-sensing mechanism, and (iii) the biased kinetics for stretched or compressed motors. The model simulations showed that the geometrical hindrance due to acto-myosin spatial mismatching and the preferential detachment of compressed motors are synergic in generating the rapid increase in the ATP-ase rate from isometric to moderate velocities of contraction, thus acting as an energy-conservation strategy in muscle contraction. The velocity measurements on a DNA-origami filament that preserves the motors distribution showed that geometrical hindrance and biased detachment generate a non-zero sliding velocity even without rotation of the myosin lever-arm, which is widely recognized as the basic event in muscle contraction. Because biased detachment is a mechanism for the rectification of thermal fluctuations, in the Brownian-ratchet framework, we predict that it requires a non-negligible amount of energy to preserve the second law of thermodynamics. Taken together, our theoretical and experimental results elucidate non-conventional components in the chemo-mechanical energy transduction in muscle.


title: Comparative Multiomic Analysis Reveals Low T Cell Infiltration as the Primary Feature of Tobacco Use in HPV(+) Oropharyngeal Cancer

abstract:
PurposeTobacco use is an independent adverse prognostic feature in human papillomavirus (HPV)-associated oropharyngeal squamous cell carcinoma (OPSCC). Despite this, the biologic features associated with tobacco use have not been systematically investigated in this population. We sought to characterize the genomic and immunologic features of HPV(+) OPSCC associated with tobacco use and adverse oncologic outcomes.

Experimental DesignWhole exome sequencing of 47 primary HPV(+) OPSCC tumors was performed to investigate mutational differences associated with tobacco exposure. To characterize the tumor immune microenvironment (TIME), targeted mRNA hybridization was performed and immunohistochemical (IHC) staining was used to validate these findings.

ResultsLow expression of transcripts in a T cell-inflamed gene expression profile (TGEP) was associated with tobacco use at the time of diagnosis and lower overall and disease-free survival. Tobacco use was associated with an increased proportion of T>C substitutions and a lower proportion of mutational signatures typically observed in HPV(+) OPSCC tumors, but was not associated with increases in mutational burden or the rate of recurrent oncogenic mutations.

ConclusionsIn HPV(+) OPSCC, low T cell infiltration of primary tumors is associated with current tobacco use and worse oncologic outcomes. Rather than an increased mutational burden, tobaccos primary and clinically relevant association is immunosuppression of the primary TIME. An objective clinical assay like the TGEP, which quantifies immune infiltration of the primary TIME, may have value for HPV(+) OPSCC risk stratification in future clinical trials.


title: Viability markers for determination of desiccation tolerance and critical stages during dehydration in Selaginella species

abstract:
Plants can tolerate some degree of dehydration but below a threshold of water content most plants die. However, some species display specific physiological, biochemical, and molecular responses that allow survival to desiccation. Some of these responses are activated at critical stages during water loss and could represent the difference between desiccation tolerance (DT) and death. Here, we report the development of a simple and reproducible system to determine DT in Selaginella species. This system is based on the use of excised tissue (explants), exposed to a dehydration agent inside small containers, rather than whole plants making it faster, better controlled, and potential use under field conditions. We also report that the triphenyltetrazolium chloride test is a simple and accurate assay to determine tissue viability. The explant system is particularly useful to identify critical points during the dehydration process and was applied to identify novel desiccation-tolerant Selaginella species. Our data suggest that desiccation-sensitive Selaginella species have a change in viability when dehydrated to 40% RWC, indicating the onset of a critical condition at this water content. Comparative studies at these critical stages could provide a better understanding of DT mechanisms and unravel insights into the key responses to survive desiccation.

HighlightIn this article, we developed a simple and efficient system to determine desiccation tolerance and critical stages during the dehydration process in Selaginella that can be applied to other plant species.


title: A polyyne toxin produced by an antagonistic bacterium blinds and lyses a green microalga

abstract: Microalgae are key contributors to global carbon fixation and the basis of many food webs. In nature, their growth is often supported or suppressed by other microorganisms. The bacterium Pseudomonas protegens Pf-5 arrests the growth of the green alga Chlamydomonas reinhardtii, deflagellates the alga by the cyclic lipopeptide orfamide A, and alters its morphology. Using a combination of Raman microspectroscopy, genome mining and mutational analysis, we discovered a novel polyyne toxin we name protegencin that is secreted by P. protegens and penetrates algal cells to destroy their primitive visual system, the eyespot. Together with secreted orfamide A, protegencin prevents the phototactic behavior of C. reinhardtii needed to perform optimal photosynthesis. A protegencin-deficient biosynthetic mutant of P. protegens does not affect growth or eyespot carotenoids of C. reinhardtii. Thus, protegencin acts in a direct and destructive way, and reveals at least a two-pronged molecular strategy used by algicidal bacteria.


title: Turning plants from passive to active material: FERONIA and microtubules independently contribute to mechanical feedback

abstract: To survive, cells must constantly resist mechanical stress. In plants, this involves the reinforcement of cell walls, notably through microtubule-dependent cellulose deposition, and thus wall sensing. Several receptor-like kinases have been proposed to act as mechanosensors. Here we tested whether the microtubule response to stress acts downstream of known wall sensors. Using a multi-step screen with eleven mutant lines, we identify FERONIA as the primary candidate for controlling the microtubule response to stress. However, when performing mechanical perturbations, we show that the microtubule response to stress can be independent from FER. We reveal that the feronia phenotype can be partially rescued by reducing tensile stress levels. Conversely, in the absence of both microtubules and FER, cells swell and burst like soap bubbles. Altogether, this shows that the microtubule response to stress acts as an independent pathway to resist stress, in parallel to FER. We propose that both pathways are key components to turn plant cells from passive to active material.


title: Proteasome homeostasis is essential for a robust cauliflower mosaic virus infection

abstract: The ubiquitin-proteasome system (UPS) is essential for the maintenance and shifts in protein homeostasis, and thereby forms a founding pillar in virtually all cellular processes including plant immunity and viral infections. According to its importance in fine-tuning the complex plant immune response, proteasomal defects result in divergent outcomes including both resistance and susceptibility phenotypes in response to viruses. The final outcome will largely depend on the specific virus and its specific co-adaptation with the UPS as well as the immune system. Here, we show that cauliflower mosaic virus (CaMV) relies on the proteasome for robust infection. The proteasome system is induced during infection via SA and supports systemic accumulation of the virus as well as plant growth performance during infection. This establishes the UPS as a win-win pathway for the plant and the virus, and together with our demonstration of a proteasome-suppressing viral effector, the intimacy between the proteasome and CaMV is fortified.


title: RNA Pol IV has antagonistic parent-of-origin effects on Arabidopsis endosperm

abstract: Gene expression in endosperm - a seed tissue that mediates transfer of maternal resources to offspring - is under complex epigenetic control. We show here that plant-specific RNA Polymerase IV mediates parental control of endosperm gene expression. Pol IV is required for the production of small interfering RNAs that typically direct DNA methylation. We compared small RNAs, DNA methylation, and mRNAs in A. thaliana endosperm from reciprocal heterozygotes produced by crossing wild-type plants to Pol IV mutants. We find that maternally and paternally acting Pol IV have divergent effects on endosperm. Losses of maternal and paternal Pol IV impact sRNAs and DNA methylation at distinct genomic sites. Strikingly, maternally and paternally-acting Pol IV have antagonistic impacts on gene expression at some loci, divergently promoting or repressing endosperm gene expression. Antagonistic parent-of-origin effects have only rarely been described and are consistent with a gene regulatory system evolving under parental conflict.


title: Transmission potential of Floridian Aedes aegypti mosquitoes for dengue virus serotype 4: Implications for estimating local dengue risk

abstract:
Dengue virus serotype 4 (DENV-4) circulated in Aedes aegypti in 2016 and 2017 in Florida in the absence of human index cases, compelling a full assessment of local mosquito vector competence and DENV-4 risk. To better understand DENV-4 transmission risk in Florida, we used an expanded suite of tests to measure and compare the vector competence of both an established colony of Ae. aegypti (Orlando strain [ORL]) and a field-derived colony from Collier County, Florida in 2018 (COL) for a Haitian DENV-4 human field isolate and a DENV-4 laboratory strain (Philippines H241). We immediately noted that ORL saliva-positivity was higher for the field versus laboratory DENV-4 strains. In a subsequent comparison with the recent COL mosquito colony we also observed significantly higher midgut susceptibility of COL and ORL for the Haitian DENV-4 field strain, and significantly higher saliva-positivity rate for COL, although overall saliva virus titers were similar between the two. These data point to a potential midgut infection barrier for the DENV-4 laboratory strain for both mosquito colonies and that the marked difference in transmission potential estimates hinge on virus-vector combinations. Our study highlights the importance of leveraging an expanded suite of testing methods with emphasis on utilizing local mosquito populations and field relevant dengue serotypes and strains to accurately estimate transmission risk in a given setting.

ImportanceDENV-4 was found circulating in Florida (FL) Ae. aegypti mosquitoes in the absence of human index cases in the state (2016-2017). How DENV-4 was maintained locally is unclear, presenting a major gap in our understanding of DENV-4 public health risk. We determined the baseline arbovirus transmission potential of laboratory and field colonies of Ae. aegypti for both laboratory and field isolates of DENV-4. We observed high transmission potential of field populations of Ae. aegypti and evidence of higher vertical transmission of the DENV-4 field isolate, providing clues to the possible mechanism of undetected DENV-4 maintenance in the state. Our findings also move the field forward in the development of best practices for evaluating arbovirus vector competence, with evidence that transmission potential estimates vary depending on the mosquito-virus combinations. These data emphasize the poor suitability of lab-established virus strains and the high relevance of field-derived mosquito populations in estimating transmission risk.


title: Machine learning models exploring characteristic single-nucleotide signatures in Yellow Fever Virus

abstract:
Yellow fever virus (YFV) is the agent of the most severe mosquito-borne disease in the tropics. Recently, Brazil suffered major YFV outbreaks with a high fatality rate affecting areas where the virus has not been reported for decades, consisting of urban areas where a large number of unvaccinated people live. We developed a machine learning framework combining three different algorithms (XGBoost, random forest and regularized logistic regression). This method was applied to 56 YFV sequences from human infections and 27 from non-human primate (NHPs) infections to investigate the presence of genetic signatures possibly related to disease severity (in human related sequences) and differences in the PCR cycle threshold (Ct) values (in NHP related sequences). Our analyses reveal four non-synonymous single nucleotide variations (SNVs) on sequences from human infections, in proteins NS3 (E614D), NS4a (I69V), NS5 (R727G, V643A) and six non-synonymous SNVs on NHP sequences, in proteins E (L385F), NS1 (A171V), NS3 (I184V) and NS5 (N11S, I374V, E641D). We performed comparative protein structural analysis on these SNVs, describing possible impacts on protein function. Despite the fact that the dataset is limited in size and that this study does not consider virus-host interactions, our work highlights the use of machine learning as a versatile and fast initial approach to genomic data exploration.

ImportanceYellow fever is responsible for 29-60 thousand deaths annually in South America and Africa and is the most severe mosquito-borne disease in the tropics. Given the range of clinical outcomes and the availability of YFV genomic data, the use of machine learning analysis promises to be a powerful tool in the investigation of genetic signatures that could impact disease severity and its potential of being reintroduced in an urban transmission cycle. This can assist in the search for biomarkers of severity as well as help elucidating variations in hosts Ct value. This work aims to propose a relatively fast and inexpensive computational analysis framework, which can be used as a real-time, innitial strategy associated with genomic surveillance to identify a set of single nucleotide variants putatively related to biological and clinical characteristics being observed.


title: Nasonia segmentation is regulated by an ancestral insect segmentation regulatory network also present in flies

abstract:
Insect segmentation is a well-studied and tractable system with which to investigate the genetic regulation of development. Though insects segment their germband using a variety of methods, modelling work implies that a single gene regulatory network can underpin the two main types of insect segmentation. This means limited genetic changes are required to explain significant differences in segmentation mode between different insects. Evidence for this idea is limited to Drosophila melanogaster, Tribolium castaneum, and the spider Parasteatoda tepidariorum, and the nature of the gene regulatory network (GRN) underlying this model has not been tested. Some insects, for example Nasonia vitripennis and Apis mellifera segment progressively, a pattern not examined in studies of this segmentation model, producing stripes at different times throughout the embryo, but not from a segment addition zone.

Here we aim to understand the GRNs patterning Nasonia using a simulation-based approach. We found that an existing model of Drosophila segmentation (Clark, 2017) can be used to recapitulate Nasonias progressive segmentation, if provided with altered inputs in the form of expression of the timer genes N{nu}-caudal and N{nu}-odd paired. We also predict limited topological changes to the pair rule network. Together this implies that very limited changes to the Drosophila network are required to simulate Nasonia segmentation, despite the differences in segmentation modes, implying that Nasonia use a very similar version of an ancestral GRN also used by Drosophila.


title: DrosoPhyla: genomic resources for drosophilid phylogeny and systematics

abstract: The vinegar fly Drosophila melanogaster is a pivotal model for invertebrate development, genetics, physiology, neuroscience, and disease. The whole family Drosophilidae, which contains over 4000 species, offers a plethora of cases for comparative and evolutionary studies. Despite a long history of phylogenetic inference, many relationships remain unresolved among the groups and genera in the Drosophilidae. To clarify these relationships, we first developed a set of new genomic markers and assembled a multilocus data set of 17 genes from 704 species of Drosophilidae. We then inferred well-supported group and species trees for this family. Additionally, we were able to determine the phylogenetic position of some previously unplaced species. These results establish a new framework for investigating the evolution of traits in fruit flies, as well as valuable resources for systematics.


title: Bayestrat: Population Stratification Correction Using Bayesian Shrinkage Prior for Genetic Association Studies

abstract: In genetic association studies with common diseases, population stratification is a major source of confounding. Principle component regression (PCR) and linear mixed model (LMM) are two commonly used approaches to account for population stratification. Previous studies have shown that LMM can be interpreted as including all principle components (PCs) as random-effect covariates. However, including all PCs in LMM may inflate type I error in some scenarios due to redundancy, while including only a few pre-selected PCs in PCR may fail to fully capture the genetic diversity. Here, we propose a statistical method under the Bayesian framework, Bayestrat, that utilizes appropriate shrinkage priors to shrink the effects of non- or minimally confounded PCs and improve the identification of highly confounded ones. Simulation results show that Bayestrat consistently achieves lower type I error rates yet higher power, especially when the number of PCs included in the model is large. We also apply our method to two real datasets, the Dallas Heart Studies (DHS) and the Multi-Ethnic Study of Atherosclerosis (MESA), and demonstrate the superiority of Bayestrat over commonly used methods.

category: genetics