Introducing 'identification probability' for automated and transferable assessment of metabolite identification confidence in metabolomics and related studies.

Methods for assessing compound identification confidence in metabolomics and related studies have been debated and actively researched for the past two decades. The earliest effort in 2007 focused primarily on mass spectrometry and nuclear magnetic resonance spectroscopy and resulted in four recommended levels of metabolite identification confidence - the Metabolite Standards Initiative (MSI) Levels. In 2014, the original MSI Levels were expanded to five levels (including two sublevels) to facilitate communication of compound identification confidence in high resolution mass spectrometry studies. Further refinement in identification levels have occurred, for example to accommodate use of ion mobility spectrometry in metabolomics workflows, and alternate approaches to communicate compound identification confidence also have been developed based on identification points schema. However, neither qualitative levels of identification confidence nor quantitative scoring systems address the degree of ambiguity in compound identifications in context of the chemical space being considered, are easily automated, or are transferable between analytical platforms. In this perspective, we propose that the metabolomics and related communities consider identification probability as an approach for automated and transferable assessment of compound identification and ambiguity in metabolomics and related studies. Identification probability is defined simply as 1/N, where N is the number of compounds in a reference library or chemical space that match to an experimentally measured molecule within user-defined measurement precision(s), for example mass measurement or retention time accuracy, etc. We demonstrate the utility of identification probability in an in silico analysis of multi-property reference libraries constructed from the Human Metabolome Database and computational property predictions, provide guidance to the community in transparent implementation of the concept, and invite the community to further evaluate this concept in parallel with their current preferred methods for assessing metabolite identification confidence.

DEIMoS: An Open-Source Tool for Processing High-Dimensional Mass Spectrometry Data.

We present DEIMoS: Data Extraction for Integrated Multidimensional Spectrometry, a Python application programming interface (API) and command-line tool for high-dimensional mass spectrometry data analysis workflows that offers ease of development and access to efficient algorithmic implementations. Functionality includes feature detection, feature alignment, collision cross section (CCS) calibration, isotope detection, and MS/MS spectral deconvolution, with the output comprising detected features aligned across study samples and characterized by mass, CCS, tandem mass spectra, and isotopic signature. Notably, DEIMoS operates on N-dimensional data, largely agnostic to acquisition instrumentation; algorithm implementations simultaneously utilize all dimensions to (i) offer greater separation between features, thus improving detection sensitivity, (ii) increase alignment/feature matching confidence among data sets, and (iii) mitigate convolution artifacts in tandem mass spectra. We demonstrate DEIMoS with LC-IMS-MS/MS metabolomics data to illustrate the advantages of a multidimensional approach in each data processing step.

Mass Spectrometry Adduct Calculator.

We describe the Mass Spectrometry Adduct Calculator (MSAC), an automated Python tool to calculate the adduct ion masses of a parent molecule. Here, adduct refers to a version of a parent molecule [M] that is charged due to addition or loss of atoms and electrons resulting in a charged ion, for example, [M + H]+. MSAC includes a database of 147 potential adducts and adduct/neutral loss combinations and their mass-to-charge ratios (m/z) as extracted from the NIST/EPA/NIH Mass Spectral Library (NIST17), Global Natural Products Social Molecular Networking Public Spectral Libraries (GNPS), and MassBank of North America (MoNA). The calculator relies on user-selected subsets of the combined database to calculate expected m/z for adducts of molecules supplied as formulas. This tool is intended to help researchers create identification libraries to collect evidence for the presence of molecules in mass spectrometry data. While the included adduct database focuses on adducts typically detected during liquid chromatography-mass spectrometry analyses, users may supply their own lists of adducts and charge states for calculating expected m/z. We also analyzed statistics on adducts from spectra contained in the three selected mass spectral libraries. MSAC is freely available at https://github.com/pnnl/MSAC.

Sex Chromosome Turnover in Bent-Toed Geckos (Cyrtodactylus).

Lizards and snakes (squamates) are known for their varied sex determining systems, and gecko lizards are especially diverse, having evolved sex chromosomes independently multiple times. While sex chromosomes frequently turnover among gecko genera, intrageneric turnovers are known only from Gekko and Hemidactylus. Here, we used RADseq to identify sex-specific markers in two species of Burmese bent-toed geckos. We uncovered XX/XY sex chromosomes in Cyrtodactylus chaunghanakwaensis and ZZ/ZW sex chromosomes in Cyrtodactylus pharbaungensis. This is the third instance of intrageneric turnover of sex chromosomes in geckos. Additionally, Cyrtodactylus are closely related to another genus with intrageneric turnover, Hemidactylus. Together, these data suggest that sex chromosome turnover may be common in this clade, setting them apart as exceptionally diverse in a group already known for diverse sex determination systems.

Simulations of Asymmetric Membranes Illustrate Cooperative Leaflet Coupling and Lipid Adaptability.

Biological membranes are composed of lipid bilayers that are often asymmetric with regards to the lipid composition and/or aqueous solvent they separate. Studying lipid asymmetry both experimentally and computationally is challenging. Molecular dynamics simulations of lipid bilayers with asymmetry are difficult due to finite system sizes and time scales accessible to simulations. Due to the very slow flip-flop rate for phospholipids, one must first choose how many lipids are on each side of the bilayer, but the resulting bilayer may be unstable (or metastable) due to differing tensile and compressive forces between leaflets. Here we use molecular dynamics simulations to investigate a number of different asymmetric membrane systems, both with atomistic and coarse-grained models. Asymmetries studied include differences in number of lipids, lipid composition (unsaturated and saturated tails and different headgroups), and chemical gradients between the aqueous phases. Extensive analysis of the bilayers' properties such as area per lipid, density, and lateral pressure profiles are used to characterize bilayer asymmetry. We also address how cholesterol (which flip-flops relatively quickly) influences membrane asymmetries. Our results show how each leaflet is influenced by the other and can mitigate the structural changes to the bilayer overall structure. Cholesterol can respond to changes in bilayer asymmetry to alleviate some of the effect on the bilayer structure, but that will alter its leaflet distribution, which in turn affects its chemical potential. Ionic imbalances are shown to have a modest change in bilayer structure, despite large changes in the electrostatic potential. Bilayer asymmetry can also induce a modest electrostatic potential across the membrane. Our results highlight the importance of membrane asymmetry on bilayer properties, the influence of lipid headgroups, tails and cholesterol on asymmetry, and the ability of lipids to adapt to different environments.

Escaping the evolutionary trap? Sex chromosome turnover in basilisks and related lizards (Corytophanidae: Squamata).

Most pleurodont lizard families (anoles, iguanas and their relatives), with the exception of the basilisks and casquehead lizards (family Corytophanidae), share homologous XX/XY sex chromosomes, syntenic with chicken chromosome 15. Here, we used a suite of methods (i.e. RADseq, RNAseq and qPCR) to identify corytophanid sex chromosomes for the first time. We reveal that all examined corytophanid species have partially degenerated XX/XY sex chromosomes, syntenic with chicken chromosome 17. Transcriptomic analyses showed that the expression of X-linked genes in the corytophanid, Basiliscus vittatus, is not balanced between the sexes, which is rather exceptional under male heterogamety, and unlike the dosage-balanced sex chromosomes in other well-studied XX/XY systems, including the green anole, Anolis carolinensis. Corytophanid sex chromosomes may represent a rare example of a turnover away from stable, differentiated sex chromosomes. However, because of poor phylogenetic resolution among pleurodont families, we cannot reject the alternative hypothesis that corytophanid sex chromosomes evolved independently from an unknown ancestral system.