No abstract

Bacterial plasmids play a large role in allowing bacteria to adapt to changing environments and can pose a significant risk to human health if they confer virulence and antimicrobial resistance (AMR). Plasmids differ significantly in the taxonomic breadth of host bacteria in which they can successfully replicate, this is commonly referred to as ‘host range’ and is usually described in qualitative terms of ‘narrow’ or ‘broad’. Understanding the host range potential of plasmids is of great interest due to their ability to disseminate traits such as AMR through bacterial populations and into human pathogens. We developed the MOB-suite to facilitate characterization of plasmids and introduced a whole-sequence-based classification system based on clustering complete plasmid sequences using Mash distances (https://github.com/phac-nml/mob-suite). We updated the MOB-suite database from 12 091 to 23 671 complete sequences, representing 17 779 unique plasmids. With advances in new algorithms for rapidly calculating average nucleotide identity (ANI), we compared clustering characteristics using two different distance measures – Mash and ANI – and three clustering algorithms on the unique set of plasmids. The plasmid nomenclature is designed to group highly similar plasmids together that are unlikely to have multiple representatives within a single cell. Based on our results, we determined that clusters generated using Mash and complete-linkage clustering at a Mash distance of 0.06 resulted in highly homogeneous clusters while maintaining cluster size. The taxonomic distribution of plasmid biomarker sequences for replication and relaxase typing, in combination with MOB-suite whole-sequence-based clusters have been examined in detail for all high-quality publicly available plasmid sequences. We have incorporated prediction of plasmid replication host range into the MOB-suite based on observed distributions of these sequence features in combination with known plasmid hosts from the literature. Host range is reported as the highest taxonomic rank that covers all of the plasmids which share replicon or relaxase biomarkers or belong to the same MOB-suite cluster code. Reporting host range based on these criteria allows for comparisons of host range between studies and provides information for plasmid surveillance.

The 18.5 kDa myelin basic protein (MBP), the most abundant splice isoform in adult mammalian myelin, is a multifunctional, intrinsically disordered protein involved in the development and compaction of the myelin sheath in the central nervous system. A highly conserved central segment comprises a membrane-anchoring amphipathic α-helix followed by a proline-rich segment that represents a ligand for SH3 domain-containing proteins. Here, we have determined using solution nuclear magnetic resonance spectroscopy the structure of a 36-residue peptide fragment of MBP (murine 18.5 kDa residues S72-S107, denoted the α2-peptide) comprising these two structural motifs, in association with dodecylphosphocholine (DPC) micelles. The structure was calculated using CS-ROSETTA (version 1.01) because the nuclear Overhauser effect restraints were insufficient for this protein. The experimental studies were complemented by molecular dynamics simulations of a corresponding 24-residue peptide fragment (murine 18.5 kDa residues E80-G103, denoted the MD-peptide), also in association with a DPC micelle in silico. The experimental and theoretical results agreed well with one another, despite the independence of the starting structures and analyses, both showing membrane association via the amphipathic α-helix, and a sharp bend in the vicinity of the Pro93 residue (murine 18.5 kDa sequence numbering). Overall, the conformations elucidated here show how the SH3 ligand is presented to the cytoplasm for interaction with SH3 domain-containing proteins such as Fyn and contribute to our understanding of myelin architecture at the molecular level.

The classic isoforms of myelin basic protein (MBP) are essential for the formation and maintenance of myelin in the central nervous system of higher vertebrates. The protein is involved in all facets of the development, compaction, and stabilization of the multilamellar myelin sheath, and also interacts with cytoskeletal and signaling proteins. The predominant 18.5-kDa isoform of MBP is an intrinsically-disordered protein that is a candidate auto-antigen in the human demyelinating disease multiple sclerosis. A highly-conserved central segment within classic MBP consists of a proline-rich region (murine 18.5-kDa sequence –T92-P93-R94-T95-P96-P97-P98-S99–) containing a putative SH3-ligand, adjacent to a region that forms an amphipathic α-helix (P82-I90) upon interaction with membranes, or under membrane-mimetic conditions. The T92 and T95 residues within the proline-rich region can be post-translationally modified through phosphorylation by mitogen-activated protein (MAP) kinases. Here, we have investigated the structure of the α-helical and proline-rich regions in dilute aqueous buffer, and have evaluated the effects of phosphorylation at T92 and T95 on the stability and dynamics of the α-helical region, by utilizing four 36-residue peptides (S72–S107) with differing phosphorylation status. Nuclear magnetic resonance spectroscopy reveals that both the α-helical as well as the proline-rich regions are disordered in aqueous buffer, whereas they are both structured in a lipid environment (cf., Ahmed et al., Biochemistry 51, 7475-9487, 2012). Thermodynamic analysis of trifluoroethanol-titration curves monitored by circular dichroism spectroscopy reveals that phosphorylation, especially at residue T92, impedes formation of the amphipathic α-helix. This conclusion is supported by molecular dynamics simulations, which further illustrate that phosphorylation reduces the folding reversibility of the α-helix upon temperature perturbation and affect the global structure of the peptides through altered electrostatic interactions. The results support the hypothesis that the central conserved segment of MBP constitutes a molecular switch in which the conformation and/or intermolecular interactions are mediated by phosphorylation/dephosphorylation at T92 and T95.

We propose a minimal protocol for exhaustive genome-wide association interaction analysis that involves screening for epistasis over large-scale genomic data combining strengths of different methods and statistical tools. The different steps of this protocol are illustrated on a real-life data application for Alzheimer's disease (AD) (2259 patients and 6017 controls from France). Particularly, in the exhaustive genome-wide epistasis screening we identified AD-associated interacting SNPs-pair from chromosome 6q11.1 (rs6455128, the KHDRBS2 gene) and 13q12.11 (rs7989332, the CRYL1 gene) (p = 0.006, corrected for multiple testing). A replication analysis in the independent AD cohort from Germany (555 patients and 824 controls) confirmed the discovered epistasis signal (p = 0.036). This signal was also supported by a meta-analysis approach in 5 independent AD cohorts that was applied in the context of epistasis for the first time. Transcriptome analysis revealed negative correlation between expression levels of KHDRBS2 and CRYL1 in both the temporal cortex (β = −0.19, p = 0.0006) and cerebellum (β = −0.23, p < 0.0001) brain regions. This is the first time a replicable epistasis associated with AD was identified using a hypothesis free screening approach.