Koonin Research Articles

Synthetic biology as a technoscience: The case of minimal genomes and essential genes

This article examines how minimal genome research mobilizes philosophical concepts such as minimality and essentiality. Following a historical approach the article aims to uncover what function this terminology plays and which problems are raised by them. Specifically, four historical moments are examined, linked to the work of Harold J. Morowitz, Mitsuhiro Itaya, Eugene Koonin and Arcady Mushegian, and J. Craig Venter. What this survey shows is a historical shift away from historical questions about life or descriptive questions about specific organisms towards questions that explore biological possibilities: what are possible forms of minimal genomes, regardless of whether they exist in nature?Moreover, it highlights a fundamental ambiguity at work in minimal genome research between a universality claim and a standardization claim: does a minimal genome refer to the minimal gene set for any organism whatsoever? Or does it refer rather to a gene set that will provide stable, robust and predictable behaviour, suited for biotechnological applications? Two diagnoses are proposed for this ambiguity: a philosophical diagnosis of how minimal genome research either misunderstands the ontology of biological entities or philosophically misarticulates scientific practice. Secondly, a historical diagnosis that suggests that this ambiguity is part of a broader shift towards technoscience.

Major Insights into Microbiology: An Interview with Luciano Marraffini.

Major Insights into Microbiology: An Interview with Luciano Marraffini.

From tumors to species: a SCANDAL hypothesis

Some tumor cells can evolve into transmissible parasites. Notable examples include the Tasmanian devil facial tumor disease, the canine transmissible venereal tumor and transmissible cancers of mollusks. We present a hypothesis that such transmissible tumors existed in the past and that some modern animal taxa are descendants of these tumors. We expect potential candidates for SCANDALs (speciated by cancer development animals) to be simplified relatives of more complex metazoans and have genomic alterations typical for cancer progression (such as deletions of universal apoptosis genes). We considered several taxa of simplified animals for our hypothesis: dicyemida, orthonectida, myxosporea and trichoplax. Based on genomic analysis we conclude that Myxosporea appear to be the most suitable candidates for a tumor ancestry. They are simplified parasitic cnidarians that universally lack major genes implicated in cancer progression including all genes with Caspase and BCL2 domains as well as any p53 and apoptotic protease activating factor – 1 (Apaf-1) homologs, suggesting the disruption of main apoptotic pathways in their early evolutionary history. Further comparative genomics and single-cell transcriptomic studies may be helpful to test our hypothesis of speciation via a cancerous stage.ReviewersThis article was reviewed by Eugene Koonin, Mikhail Gelfand and Gregory M Woods.

Darwinian selection of host and bacteria supports emergence of Lamarckian-like adaptation of the system as a whole

BackgroundThe relatively fast selection of symbiotic bacteria within hosts and the potential transmission of these bacteria across generations of hosts raise the question of whether interactions between host and bacteria support emergent adaptive capabilities beyond those of germ-free hosts.ResultsTo investigate possibilities for emergent adaptations that may distinguish composite host-microbiome systems from germ-free hosts, we introduce a population genetics model of a host-microbiome system with vertical transmission of bacteria. The host and its bacteria are jointly exposed to a toxic agent, creating a toxic stress that can be alleviated by selection of resistant individuals and by secretion of a detoxification agent (“detox”). We show that toxic exposure in one generation of hosts leads to selection of resistant bacteria, which in turn, increases the toxic tolerance of the host’s offspring. Prolonged exposure to toxin over many host generations promotes anadditional form of emergent adaptation due to selection of hosts based on detox produced by their bacterial community as a whole (as opposed to properties of individual bacteria).ConclusionsThese findings show that interactions between pure Darwinian selections of host and its bacteria can give rise to emergent adaptive capabilities, including Lamarckian-like adaptation of the host-microbiome system.ReviewersThis article was reviewed by Eugene Koonin, Yuri Wolf and Philippe Huneman.

Biology as a Constructive Physics

Yuri Manin’s approach to Zipf’s law (Kolmogorov complexity as energy) is applied to investigation of biological evolution. Model of constructive statistical mechanics where complexity is a contribution to energy is proposed to model genomics. Scaling laws in genomics are discussed in relation to Zipf’s law. This gives a model of Eugene Koonin’s Third Evolutionary Synthesis – physical model which should describe scaling in genomics.

Molecular mechanisms of the juvenile form of Batten disease: important role of MAPK signaling pathways (ERK1/ERK2, JNK and p38) in pathogenesis of the malady

BackgroundMutations in the CLN3 gene lead to so far an incurable juvenile-onset neuronal ceroid lipofuscinosis (JNCL) or Batten disease that starts at the age of 4–6 years with a progressive retinopathy leading to blindness. Motor disturbances, epilepsy and dementia manifest during several following years. Most JNCL patients carry the same 1.02-kb deletion in the CLN3 gene, encoding an unusual transmembrane protein, CLN3 or battenin.ResultsBased on data of genome-wide expression profiling in CLN3 patients with different rate of the disease progression [Mol. Med., 2011, 17: 1253–1261] and our bioinformatic analysis of battenin protein-protein interactions in neurons we propose that CLN3 can function as a molecular chaperone for some plasma membrane proteins, being crucially important for their correct folding in endoplasmic reticulum. Changes in spatial structure of these membrane proteins lead to transactivation of the located nearby receptors. Particularly, CLN3 interacts with a subunit of Na/K ATPase ATP1A1 which changes its conformation and activates the adjacent epidermal growth factor receptor (EGFR). As a result, a large amount of erroneously activated EGFR generates MAPK signal cascades (ERK1/ERK2, JNKs and p38) from cell surface eventually causing neurons’ death.ConclusionsMolecular mechanism of the juvenile form of Batten disease (JNCL), which is based on the excessive activation of signaling cascades in a time of the radical increase of neuronal membranes’ area in the growing brain, have been proposed and substantiated. The primary cause of this phenomenon is the defective function of the CLN3 protein that could not act properly as molecular chaperone for some plasma membrane proteins in the endoplasmic reticulum. The incorrect three-dimensional structure of at least one such protein, ATP1A1, leads to unregulated spontaneous and repetitive activation of the SRC kinase that transactivates EGFR with the subsequent uncontrolled launch of various MAPK cascades. Possible ways of treatment of patients with JNCL have been suggested.ReviewersThis article was reviewed by Konstantinos Lefkimmiatis, Eugene Koonin and Vladimir Poroikov.

ASXL gain-of-function truncation mutants: defective and dysregulated forms of a natural ribosomal frameshifting product?

BackgroundProgrammed ribosomal frameshifting (PRF) is a gene expression mechanism which enables the translation of two N-terminally coincident, C-terminally distinct protein products from a single mRNA. Many viruses utilize PRF to control or regulate gene expression, but very few phylogenetically conserved examples are known in vertebrate genes. Additional sex combs-like (ASXL) genes 1 and 2 encode important epigenetic and transcriptional regulatory proteins that control the expression of homeotic genes during key developmental stages. Here we describe an ~150-codon overlapping ORF (termed TF) in ASXL1 and ASXL2 that, with few exceptions, is conserved throughout vertebrates.ResultsConservation of the TF ORF, strong suppression of synonymous site variation in the overlap region, and the completely conserved presence of an EH[N/S]Y motif (a known binding site for Host Cell Factor-1, HCF-1, an epigenetic regulatory factor), all indicate that TF is a protein-coding sequence. A highly conserved UCC_UUU_CGU sequence (identical to the known site of +1 ribosomal frameshifting for influenza virus PA-X expression) occurs at the 5′ end of the region of enhanced synonymous site conservation in ASXL1. Similarly, a highly conserved RG_GUC_UCU sequence (identical to a known site of −2 ribosomal frameshifting for arterivirus nsp2TF expression) occurs at the 5′ end of the region of enhanced synonymous site conservation in ASXL2.ConclusionsDue to a lack of appropriate splice forms, or initiation sites, the most plausible mechanism for translation of the ASXL1 and 2 TF regions is ribosomal frameshifting, resulting in a transframe fusion of the N-terminal half of ASXL1 or 2 to the TF product, termed ASXL-TF. Truncation or frameshift mutants of ASXL are linked to myeloid malignancies and genetic diseases, such as Bohring-Opitz syndrome, likely at least in part as a result of gain-of-function or dominant-negative effects. Our hypothesis now indicates that these disease-associated mutant forms represent overexpressed defective versions of ASXL-TF.ReviewersThis article was reviewed by Laurence Hurst and Eugene Koonin.

On the global CRISPR array behavior in class I systems

BackgroundMuch effort is underway to build and upgrade databases and tools related to occurrence, diversity, and characterization of CRISPR-Cas systems. As microbial communities and their genome complements are unearthed, much emphasis has been placed on details of individual strains and model systems within the CRISPR-Cas classification, and that collection of information as a whole affords the opportunity to analyze CRISPR-Cas systems from a quantitative perspective to gain insight into distribution of CRISPR array sizes across the different classes, types and subtypes. CRISPR diversity, nomenclature, occurrence, and biological functions have generated a plethora of data that created a need to understand the size and distribution of these various systems to appreciate their features and complexity.ResultsBy utilizing a statistical framework and visual analytic techniques, we have been able to test several hypotheses about CRISPR loci in bacterial class I systems. Quantitatively, though CRISPR loci can expand to hundreds of spacers, the mean and median sizes are 40 and 25, respectively, reflecting rather modest acquisition and/or retention overall. Histograms uncovered that CRISPR array size displayed a parametric distribution, which was confirmed by a goodness-of fit test. Mapping the frequency of CRISPR loci on a standardized chromosome plot revealed that CRISPRs have a higher probability of occurring at clustered locations along the positive or negative strand. Lastly, when multiple arrays occur in a particular system, the size of a particular CRISPR array varies with its distance from the cas operon, reflecting acquisition and expansion biases.ConclusionsThis study establishes that bacterial Class I CRISPR array size tends to follow a geometric distribution; these CRISPRs are not randomly distributed along the chromosome; and the CRISPR array closest to the cas genes is typically larger than loci in trans. Overall, we provide an analytical framework to understand the features and behavior of CRISPR-Cas systems through a quantitative lens.ReviewersThis article was reviewed by Eugene Koonin (NIH-NCBI) and Uri Gophna (Tel Aviv University).

Demystification of animal symmetry: symmetry is a response to mechanical forces

ᅟSymmetry is an eye-catching feature of animal body plans, yet its causes are not well enough understood. The evolution of animal form is mainly due to changes in gene regulatory networks (GRNs). Based on theoretical considerations regarding fundamental GRN properties, it has recently been proposed that the animal genome, on large time scales, should be regarded as a system which can construct both the main symmetries – radial and bilateral – simultaneously; and that the expression of any of these depends on functional constraints. Current theories explain biological symmetry as a pattern mostly determined by phylogenetic constraints, and more by chance than by necessity. In contrast to this conception, I suggest that physical effects, which in many cases act as proximate, direct, tissue-shaping factors during ontogenesis, are also the ultimate causes – i.e. the indirect factors which provide a selective advantage – of animal symmetry, from organs to body plan level patterns. In this respect, animal symmetry is a necessary product of evolution. This proposition offers a parsimonious view of symmetry as a basic feature of the animal body plan, suggesting that molecules and physical forces act in a beautiful harmony to create symmetrical structures, but that the concert itself is directed by the latter.ReviewersThis article was reviewed by Eugene Koonin, Zoltán Varga and Michaël Manuel.

Expansion of the molecular and morphological diversity of Acanthamoebidae (Centramoebida, Amoebozoa) and identification of a novel life cycle type within the group.

BackgroundAcanthamoebidae is a “family” level amoebozoan group composed of the genera Acanthamoeba, Protacanthamoeba, and very recently Luapeleamoeba. This clade of amoebozoans has received considerable attention from the broader scientific community as Acanthamoeba spp. represent both model organisms and human pathogens. While the classical composition of the group (Acanthamoeba + Protacanthamoeba) has been well accepted due to the morphological and ultrastructural similarities of its members, the Acanthamoebidae has never been highly statistically supported in single gene phylogenetic reconstructions of Amoebozoa either by maximum likelihood (ML) or Bayesian analyses.ResultsHere we show using a phylogenomic approach that the Acanthamoebidae is a fully supported monophyletic group within Amoebozoa with both ML and Bayesian analyses. We also expand the known range of morphological and life cycle diversity found in the Acanthamoebidae by demonstrating that the amoebozoans “Protostelium” arachisporum, Dracoamoeba jormungandri n. g. n. sp., and Vacuolamoeba acanthoformis n.g. n.sp., belong within the group. We also found that “Protostelium” pyriformis is clearly a species of Acanthamoeba making it the first reported sporocarpic member of the genus, that is, an amoeba that individually forms a walled, dormant propagule elevated by a non-cellular stalk. Our phylogenetic analyses recover a fully supported Acanthamoebidae composed of five genera. Two of these genera (Acanthamoeba and Luapeleameoba) have members that are sporocarpic.ConclusionsOur results provide high statistical support for an Acanthamoebidae that is composed of five distinct genera. This study increases the known morphological diversity of this group and shows that species of Acanthamoeba can include spore-bearing stages. This further illustrates the widespread nature of spore-bearing stages across the tree of Amoebozoa.ReviewersThis article was reviewed by Drs. Eugene Koonin, Purificacion Lopez-Garcia and Sandra Baldauf. Sandra Baldauf was nominated by Purificacion Lopez-Garcia, an Editorial Board member.Electronic supplementary materialThe online version of this article (doi:10.1186/s13062-016-0171-0) contains supplementary material, which is available to authorized users.

Cpf1 nucleases demonstrate robust activity to induce DNA modification by exploiting homology directed repair pathways in mammalian cells.

BackgroundCpf1 nucleases have recently been repurposed for site-specific genome modification. Two members of the Cpf1 family, the AsCpf1 from Acidaminococcus sp. and the LbCpf1 from Lachnospiraceae bacterium were shown to induce higher indel frequencies than SpCas9 when examining four randomly-selected target sequences for each type of nuclease. Whether they are a real match for Cas9 nucleases, however, remains to be verified.ResultsHere, we used AsCpf1 and LbCpf1 to induce homology directed repair, either single strand annealing (SSA) or homologous recombination (HR), in N2a mouse neuroblastoma cells. Exploiting a plasmid that contains two GFP halves with overlapping sequences and exploring 20 targets, on all but one both nucleases consistently performed with above 10 % efficiency. Several Cas9 nucleases have been previously characterised in order to find an orthogonal counterpart for the most widely used promiscuous SpCas9. Here, we found that AsCpf1 and LbCpf1 might be better candidates than three of the best such counterparts: Cas9 from Staphylococcus aureus, from Streptococcus thermophilus and from Neisseria meningitidis, when assessed for inducing efficient SSA mediated repair in N2a cells. When tested on genomic targets exploiting HR, both nucleases were able to induce the integration of a donor cassette with 1000 bp-long homologous arms. We also generated plasmids that express these Cpf1 nucleases together with their cognate crRNAs and that are equipped with type IIS restriction enzyme sites to facilitate spacer cloning.ConclusionsOur results suggest that employing As- or LbCpf1 nuclease to induce homology directed repair in N2a cells, although is less effective at present than employing SpCas9, it is an equally or more effective tool than the most frequently used orthogonal Cas9 counterparts of SpCas9. These findings support the position of Cpf1 nucleases on the side of SpCas9 on the palette of effective genome engineering tools.ReviewersThis article was reviewed by Eugene Koonin, Haruhiko Siomi and Jean-Yves Masson.Electronic supplementary materialThe online version of this article (doi:10.1186/s13062-016-0147-0) contains supplementary material, which is available to authorized users.

Evolutionary consequences of polyploidy in prokaryotes and the origin of mitosis and meiosis.

BackgroundThe origin of eukaryote-specific traits such as mitosis and sexual reproduction remains disputable. There is growing evidence that both mitosis and eukaryotic sex (i.e., the alternation of syngamy and meiosis) may have already existed in the basal eukaryotes. The mating system of the halophilic archaeon Haloferax volcanii probably represents an intermediate stage between typical prokaryotic and eukaryotic sex. H. volcanii is highly polyploid, as well as many other Archaea. Here, we use computer simulation to explore genetic and evolutionary outcomes of polyploidy in amitotic prokaryotes and its possible role in the origin of mitosis, meiosis and eukaryotic sex.ResultsModeling suggests that polyploidy can confer strong short-term evolutionary advantage to amitotic prokaryotes. However, it also promotes the accumulation of recessive deleterious mutations and the risk of extinction in the long term, especially in highly mutagenic environment. There are several possible strategies that amitotic polyploids can use in order to reduce the genetic costs of polyploidy while retaining its benefits. Interestingly, most of these strategies resemble different components or aspects of eukaryotic sex. They include asexual ploidy cycles, equalization of genome copies by gene conversion, high-frequency lateral gene transfer between relatives, chromosome exchange coupled with homologous recombination, and the evolution of more accurate chromosome distribution during cell division (mitosis). Acquisition of mitosis by an amitotic polyploid results in chromosome diversification and specialization. Ultimately, it transforms a polyploid cell into a functionally monoploid one with multiple unique, highly redundant chromosomes. Specialization of chromosomes makes the previously evolved modes of promiscuous chromosome shuffling deleterious. This can result in selective pressure to develop accurate mechanisms of homolog pairing, and, ultimately, meiosis.ConclusionEmergence of mitosis and the first evolutionary steps towards eukaryotic sex could have taken place in the ancestral polyploid, amitotic proto-eukaryotes, as they were struggling to survive in the highly mutagenic environment of the Early Proterozoic shallow water microbial communities, through the succession of the following stages: (1) acquisition of high-frequency between-individual genetic exchange coupled with homologous recombination; (2) acquisition of mitosis, followed by rapid chromosome diversification and specialization; (3) evolution of homolog synapsis and meiosis. Additional evidence compatible with this scenario includes mass acquisition of new families of paralogous genes by the basal eukaryotes, and recently discovered correlation between polyploidy and the presence of histones in Archaea.ReviewerThis article was reviewed by Eugene Koonin, Uri Gophna and Armen Mulkidjanian. For the full reviews, please go to the Reviewers' comments section.

Computational recognition and analysis of hitherto uncharacterized nucleotide cyclase-like proteins in bacteria.

Evolutionary relationship between class III nucleotide cyclases and an uncharacterized set of bacterial proteins from Actinobacteria, Bacteroidetes and Proteobacteria has been recognized and analyzed. Detailed analyses of sequence and structural features resulted in the recognition of potential cyclase function conferring residues and presence of signature topological motif (βααββαβ) in the uncharacterized set of bacterial proteins. Lack of transmembrane domains and signal peptide cleavage sites is suggestive of their cytosolic subcellular localization. Furthermore, analysis on evolutionarily conserved gene clusters of the predicted nucleotide cyclase-like proteins and their evolutionary relationship with nucleotide cyclases suggest their participation in cellular signalling events. Our analyses suggest expansion of class III nucleotide cyclases.ReviewersThis article was reviewed by Eugene Koonin and Michael Gromiha.Electronic supplementary materialThe online version of this article (doi:10.1186/s13062-016-0130-9) contains supplementary material, which is available to authorized users.

Estimation of ribosome profiling performance and reproducibility at various levels of resolution.

BackgroundRibosome profiling (or Ribo-seq) is currently the most popular methodology for studying translation; it has been employed in recent years to decipher various fundamental gene expression regulation aspects.The main promise of the approach is its ability to detect ribosome densities over an entire transcriptome in high resolution of single codons. Indeed, dozens of ribo-seq studies have included results related to local ribosome densities in different parts of the transcript; nevertheless, the performance of Ribo-seq has yet to be quantitatively evaluated and reported in a large-scale multi-organismal and multi-protocol study of currently available datasets.ResultsHere we provide the first objective evaluation of Ribo-seq at the resolution of a single nucleotide(s) using clear, interpretable measures, based on the analysis of 15 experiments, 6 organisms, and a total of 612, 961 transcripts. Our major conclusion is that the ability to infer signals of ribosomal densities at nucleotide scale is considerably lower than previously thought, as signals at this level are not reproduced well in experimental replicates. In addition, we provide various quantitative measures that connect the expected error rate with Ribo-seq analysis resolution.ConclusionsThe analysis of Ribo-seq data at the resolution of codons and nucleotides provides a challenging task, calls for task-specific statistical methods and further protocol improvements. We believe that our results are important for every researcher studying translation and specifically for researchers analyzing data generated by the Ribo-seq approach.ReviewersThis article was reviewed by Dmitrij Frishman, Eugene Koonin and Frank Eisenhaber.Electronic supplementary materialThe online version of this article (doi:10.1186/s13062-016-0127-4) contains supplementary material, which is available to authorized users.

Lack of population differentiation patterns of previously identified putatively adaptive transposable element insertions at microgeographic scales.

BackgroundTransposable elements (TEs) play an important role in genome function and evolution. It has been shown that TEs are a considerable source of adaptive changes in the genome of Drosophila melanogaster. Specifically, footprints of selection at the DNA level, the presence of population differentiation patterns across environmental gradients, and detailed mechanistic and fitness analyses of a few candidate adaptive TEs pointed to the role of TEs in environmental adaptation. However, whether the population differentiation patterns observed at large geographic scales can be replicated at a microgeographic scale has never been assessed before.ResultsIn this work, we explored the population patterns of putatively adaptive TEs at a micro-spatial scale level. We compared the frequencies of TEs, previously identified as putatively adaptive and putatively neutral, in populations collected in opposite slopes of the Evolution Canyon at Mt. Carmel in Israel separated by 200 m on average. However, the differentiation patterns previously observed across large geographic distances (2000–2200 km) were not replicated at the microscale level of the Evolution Canyon populations.ConclusionTE insertions previously associated with D. melanogaster environmental adaptation at a macro scale level do not play such a role at the microscale level of the Evolution Canyon populations. However, these results do not exclude a role of TEs in microgeographic adaptation because the dataset analyzed in this work is restricted to TEs identified in a single North American strain and as such is highly biased and incomplete.ReviewersThis article was reviewed by Eugene Koonin, Limsoon Wong and Fyodor Kondrashov.

Sequence analysis reveals a conserved extension in the capping enzyme of the alphavirus supergroup, and a homologous domain in nodaviruses.

BackgroundMembers of the alphavirus supergroup include human pathogens such as chikungunya virus, hepatitis E virus and rubella virus. They encode a capping enzyme with methyltransferase-guanylyltransferase (MTase-GTase) activity, which is an attractive drug target owing to its unique mechanism. However, its experimental study has proven very difficult.ResultsWe examined over 50 genera of viruses by sequence analyses. Earlier studies showed that the MTase-GTase contains a “Core” region conserved in sequence. We show that it is followed by a long extension, which we termed “Iceberg” region, whose secondary structure, but not sequence, is strikingly conserved throughout the alphavirus supergroup. Sequence analyses strongly suggest that the minimal capping domain corresponds to the Core and Iceberg regions combined, which is supported by earlier experimental data. The Iceberg region contains all known membrane association sites that contribute to the assembly of viral replication factories. We predict that it may also contain an overlooked, widely conserved membrane-binding amphipathic helix. Unexpectedly, we detected a sequence homolog of the alphavirus MTase-GTase in taxa related to nodaviruses and to chronic bee paralysis virus. The presence of a capping enzyme in nodaviruses is biologically consistent, since they have capped genomes but replicate in the cytoplasm, where no cellular capping enzyme is present. The putative MTase-GTase domain of nodaviruses also contains membrane-binding sites that may drive the assembly of viral replication factories, revealing an unsuspected parallel with the alphavirus supergroup.ConclusionsOur work will guide the functional analysis of the alphaviral MTase-GTase and the production of domains for structure determination. The identification of a homologous domain in a simple model system, nodaviruses, which replicate in numerous eukaryotic cell systems (yeast, flies, worms, mammals, and plants), can further help crack the function and structure of the enzyme.ReviewersThis article was reviewed by Valerian Dolja, Eugene Koonin and Sebastian Maurer-Stroh.Electronic supplementary materialThe online version of this article (doi:10.1186/s13062-015-0050-0) contains supplementary material, which is available to authorized users.

"Off-Spotter": very fast and exhaustive enumeration of genomic lookalikes for designing CRISPR/Cas guide RNAs.

BackgroundCRISPR/Cas (Clustered Regularly Interspaced Short Palindromic Repeats/CRISPR associated nucleases) is a powerful component of the prokaryotic immune system that has been adapted for targeted genetic engineering in higher organisms. A key element of CRISPR/Cas is the “guide” RNA (gRNA) that is ~20 nucleotides (nts) in length and designed to be complementary to the intended target site. An integral requirement of the CRISPR/Cas system is that the target site be followed by a protospacer adjacent motif (PAM). Care needs to be exercised during gRNA design to avoid unintended (“off-target”) interactions.ResultsWe designed and implemented the Off-Spotter algorithm to assist with the design of optimal gRNAs. When presented with a candidate gRNA sequence and a PAM, Off-Spotter quickly and exhaustively identifies all genomic sites that satisfy the PAM constraint and are identical or nearly-identical to the provided gRNA. Off-Spotter achieves its extreme performance through purely algorithmic means and not through hardware accelerators such as graphical processing units (GPUs). Off-Spotter also allows the user to identify on-the-fly how many and which nucleotides of the gRNA comprise the “seed”. Off-Spotter’s output includes a histogram showing the number of potential off-targets as a function of the number of mismatches. The output also includes for each potential off-target the site’s genomic location, a human genome browser hyperlink to the corresponding location, genomic annotation in the vicinity of the off-target, GC content, etc.ConclusionOff-Spotter is very fast and flexible and can help in the design of optimal gRNAs by providing several PAM choices, a run-time definition of the seed and of the allowed number of mismatches, and a flexible output interface that allows sorting of the results, optional viewing/hiding of columns, etc. A key element of Off-Spotter is that it does not have a rigid definition of the seed: instead, the user can declare both the seed’s location and extent on-the-fly. We expect that this flexibility in combination with Off-Spotter’s speed and richly annotated output will enable experimenters to interactively and quickly explore different scenarios and gRNA possibilities.ReviewedThis article was reviewed by Dr Eugene Koonin and Dr Frank Eisenhaber.Electronic supplementary materialThe online version of this article (doi:10.1186/s13062-015-0035-z) contains supplementary material, which is available to authorized users.

Mammalian Argonaute-DNA binding?

When a field shares the consensus that a particular phenomenon does NOT occur, this may reflect extensive experimental investigations with negative outcomes, or may represent the “common sense” position based on current knowledge and established ways of thinking. The current consensus of the RNA field is that eukaryotic Argonaute (Ago) proteins employ RNA guides and target other RNAs. The alternative -- that eukaryotic Ago has biologically important interactions with DNA in vivo – has not been seriously considered, in part because the only role contemplated for DNA was as a guide strand, and in part because it did not seem plausible that any natural source of suitable DNAs exists in eukaryotic cells. However, eukaryotic Argonaute domains bind DNA in the test tube, and several articles report that small inhibitory double-stranded DNAs do have the ability to silence target RNAs in a sequence-dependent (though poorly characterized) manner. A search of the literature identified potential DNA binding partners for Ago, including (among others) single-stranded DNAs residing in extracellular vesicles, and cytoplasmic satellite-repeat DNA fragments that are associated with the plasma membrane and transcribed by Pol II. It is interesting to note that both cytoplasmic and extracellular vesicle DNA are expressed at greatly elevated levels in cancer cells relative to normal cells. In such a pathological scenario, if not under normal conditions, there may be appreciable binding of Ago to DNA despite its lower affinity compared to RNA. If so, DNA might displace Ago from binding to its normal partners (miRNAs, siRNAs and other short ncRNAs), disrupting tightly controlled post-transcriptional gene silencing processes that are vital to correct functioning of a normal cell. The possible contribution to cancer pathogenesis is a strong motivator for further investigation of Ago-DNA binding. More generally, this case underscores the need for better informatics tools to allow investigators to analyze the state of a given scientific question at a high-level and to identify possible new research directions.Reviewers: This article was reviewed by Eugene Koonin, Kira S. Makarova, Alexander Maxwell Burroughs (nominated by L Aravind), and Isidore Rigoutsos.Open peer review: Reviewed by Eugene Koonin, Kira S. Makarova, Alexander Maxwell Burroughs (nominated by L Aravind), and Isidore Rigoutsos. For the full reviews, please go to the Reviewers’ comments section.

The fundamental tradeoff in genomes and proteomes of prokaryotes established by the genetic code, codon entropy, and physics of nucleic acids and proteins.

BackgroundMutations in nucleotide sequences provide a foundation for genetic variability, and selection is the driving force of the evolution and molecular adaptation. Despite considerable progress in the understanding of selective forces and their compositional determinants, the very nature of underlying mutational biases remains unclear.ResultsWe explore here a fundamental tradeoff, which analytically describes mutual adjustment of the nucleotide and amino acid compositions and its possible effect on the mutational biases. The tradeoff is determined by the interplay between the genetic code, optimization of the codon entropy, and demands on the structure and stability of nucleic acids and proteins.ConclusionThe tradeoff is the unifying property of all prokaryotes regardless of the differences in their phylogenies, life styles, and extreme environments. It underlies mutational biases characteristic for genomes with different nucleotide and amino acid compositions, providing foundation for evolution and adaptation.ReviewersThis article was reviewed by Eugene Koonin, Michael Gromiha, and Alexander Schleiffer.Electronic supplementary materialThe online version of this article (doi:10.1186/s13062-014-0029-2) contains supplementary material, which is available to authorized users.

Microbes on the edge of Occam's razor.

Our body harbors hundreds of microbial species and contains many more bacterial than human cells. These microbes are not passive riders but rather a vital component of the organism. The human microbiota affects our health in multiple ways, both positively and negatively. One of the new attractive directions in microbiome biology is the “microbiome-brain axis”. Several groups of researchers have described the ability of the gut microbiota to communicate with the brain and thus modulate human behavior. These limited experimental data became the foundation of the “biomeme hypothesis” of possible microbial origin of some religious rituals that has recently appeared in Biology Direct. Here I propose a critical analysis of this hypothesis. I conclude that there is no evidence of the microbial origin of religious practices but there are strong indications of their psychological and social roots.ReviewersThis article was reviewed by Eugene Koonin, Neil R Smalheiser, Etienne Joly.