Compensatory mutations cause excess of antagonistic epistasis in RNA secondary structure folding

Background

The rate at which fitness declines as an organism's genome accumulates random mutations is an important variable in several evolutionary theories. At an intuitive level, it might seem natural that random mutations should tend to interact synergistically, such that the rate of mean fitness decline accelerates as the number of random mutations is increased. However, in a number of recent studies, a prevalence of antagonistic epistasis (the tendency of multiple mutations to have a mitigating rather than reinforcing effect) has been observed.

Results

We studied in silico the net amount and form of epistatic interactions in RNA secondary structure folding by measuring the fraction of neutral mutants as a function of mutational distance d. We found a clear prevalence of antagonistic epistasis in RNA secondary structure folding. By relating the fraction of neutral mutants at distance d to the average neutrality at distance d, we showed that this prevalence derives from the existence of many compensatory mutations at larger mutational distances.

Conclusions

Our findings imply that the average direction of epistasis in simple fitness landscapes is directly related to the density with which fitness peaks are distributed in these landscapes.

     

COGcollator: a web server for analysis of distant relationships between homologous protein families

Background

The Clusters of Orthologous Groups (COGs) of proteins systematize evolutionary related proteins into specific groups with similar functions. However, the available databases do not provide means to assess the extent of similarity between the COGs.

Aim

We intended to provide a method for identification and visualization of evolutionary relationships between the COGs, as well as a respective web server.

Results

Here we introduce the COGcollator, a web tool for identification of evolutionarily related COGs and their further analysis. We demonstrate the utility of this tool by identifying the COGs that contain distant homologs of (i) the catalytic subunit of bacterial rotary membrane ATP synthases and (ii) the DNA/RNA helicases of the superfamily 1.

Reviewers

This article was reviewed by Drs. Igor N. Berezovsky, Igor Zhulin and Yuri Wolf.

     

RSEARCH: Finding homologs of single structured RNA sequences

Background

For many RNA molecules, secondary structure rather than primary sequence is the evolutionarily conserved feature. No programs have yet been published that allow searching a sequence database for homologs of a single RNA molecule on the basis of secondary structure.

Results

We have developed a program, RSEARCH, that takes a single RNA sequence with its secondary structure and utilizes a local alignment algorithm to search a database for homologous RNAs. For this purpose, we have developed a series of base pair and single nucleotide substitution matrices for RNA sequences called RIBOSUM matrices. RSEARCH reports the statistical confidence for each hit as well as the structural alignment of the hit. We show several examples in which RSEARCH outperforms the primary sequence search programs BLAST and SSEARCH. The primary drawback of the program is that it is slow. The C code for RSEARCH is freely available from our lab's website.

Conclusion

RSEARCH outperforms primary sequence programs in finding homologs of structured RNA sequences.

     

Type 2 diabetes increases oocyte mtDNA mutations which are eliminated in the offspring by bottleneck effect

Background

Diabetes induces many complications including reduced fertility and low oocyte quality, but whether it causes increased mtDNA mutations is unknown.

Methods

We generated a T2D mouse model by using high-fat-diet (HFD) and Streptozotocin (STZ) injection. We examined mtDNA mutations in oocytes of diabetic mice by high-throughput sequencing techniques.

Results

T2D mice showed glucose intolerance, insulin resistance, low fecundity compared to the control group. T2D oocytes showed increased mtDNA mutation sites and mutation numbers compared to the control counterparts. mtDNA mutation examination in F1 mice showed that the mitochondrial bottleneck could eliminate mtDNA mutations.

Conclusions

T2D mice have increased mtDNA mutation sites and mtDNA mutation numbers in oocytes compared to the counterparts, while these adverse effects can be eliminated by the bottleneck effect in their offspring. This is the first study using a small number of oocytes to examine mtDNA mutations in diabetic mothers and offspring.

     

RiboTagger: fast and unbiased 16S/18S profiling using whole community shotgun metagenomic or metatranscriptome surveys

Background

Taxonomic profiling of microbial communities is often performed using small subunit ribosomal RNA (SSU) amplicon sequencing (16S or 18S), while environmental shotgun sequencing is often focused on functional analysis. Large shotgun datasets contain a significant number of SSU sequences and these can be exploited to perform an unbiased SSU--based taxonomic analysis.

Results

Here we present a new program called RiboTagger that identifies and extracts taxonomically informative ribotags located in a specified variable region of the SSU gene in a high-throughput fashion.

Conclusions

RiboTagger permits fast recovery of SSU-RNA sequences from shotgun nucleic acid surveys of complex microbial communities. The program targets all three domains of life, exhibits high sensitivity and specificity and is substantially faster than comparable programs.

     

Genetic mutations in influenza H3N2 viruses from a 2012 epidemic in Southern China

Background

An influenza H3N2 epidemic occurred throughout Southern China in 2012.

Methods

We analyzed the hemagglutinin (HA) and neuraminidase (NA) genes of influenza H3N2 strains isolated between 2011–2012 from Guangdong. Mutation sites, evolutionary selection, antigenic sites, and N-glycosylation within these strains were analyzed.

Results

The 2011–2012 Guangdong strains contained the HA-A214S, HA-V239I, HA-N328S, NA-L81P, and NA-D93G mutations, similar to those seen in the A/ Perth/16/2009 influenza strain. The HA-NSS_061–063 and NNS_160–162 glycosylation sites were prevalent among the 2011–2012 Guangdong strains but the NA-NRS_402–404 site was deleted. Antigenically, there was a four-fold difference between A/Perth/16/2009 -like strains and the 2011–2012 Guangdong strains.

Conclusion

Antigenic drift of the H3N2 subtype contributed to the occurrence of the Southern China influenza epidemic of 2012.

     

Deep learning for DNase I hypersensitive sites identification

Background

The DNase I hypersensitive sites (DHSs) are associated with the cis-regulatory DNA elements. An efficient method of identifying DHSs can enhance the understanding on the accessibility of chromatin. Despite a multitude of resources available on line including experimental datasets and computational tools, the complex language of DHSs remains incompletely understood.

Methods

Here, we address this challenge using an approach based on a state-of-the-art machine learning method. We present a novel convolutional neural network (CNN) which combined Inception like networks with a gating mechanism for the response of multiple patterns and longterm association in DNA sequences to predict multi-scale DHSs in Arabidopsis, rice and Homo sapiens.

Results

Our method obtains 0.961 area under curve (AUC) on Arabidopsis, 0.969 AUC on rice and 0.918 AUC on Homo sapiens.

Conclusions

Our method provides an efficient and accurate way to identify multi-scale DHSs sequences by deep learning.

     

Controlled dehydration improves the diffraction quality of two RNA crystals

Background

Post-crystallization dehydration methods, applying either vapor diffusion or humidity control devices, have been widely used to improve the diffraction quality of protein crystals. Despite the fact that RNA crystals tend to diffract poorly, there is a dearth of reports on the application of dehydration methods to improve the diffraction quality of RNA crystals.

Results

We use dehydration techniques with a Free Mounting System (FMS, a humidity control device) to recover the poor diffraction quality of RNA crystals. These approaches were applied to RNA constructs that model various RNA-mediated repeat expansion disorders.

Conclusion

The method we describe herein could serve as a general tool to improve diffraction quality of RNA crystals to facilitate structure determinations.

     

RNA binding protein Lin28A promotes osteocarcinoma cells progression by associating with the long noncoding RNA MALAT1

Objectives

To explore the effects of Lin28A on progression of osteocarcinoma (OS) cells.

Results

Lin28A mRNA and protein expressions were significantly increased in OS tissues compared with that in normal adjacent tissues. Expressions of Lin28A and long noncoding RNA MALAT1 were positively correlated. Patients with higher Lin28A expression had shorter overall survival. Moreover, Lin28A knockdown inhibited OS cells proliferation, migration, invasion and promoted cell apoptosis; Lin28A was found to harbor binding sites on MALAT1 sequences and associated with MALAT1, and increased MALAT1 stability and expression. Notably, the inhibition of Lin28A knockdown was attenuated or even reversed by MALAT1 overexpression.

Conclusions

RNA binding protein Lin28A could facilitate OS cells progression by associating with the long noncoding RNA MALAT1.

     

Reconstruction of ancestral RNA sequences under multiple structural constraints

Background

Secondary structures form the scaffold of multiple sequence alignment of non-coding RNA (ncRNA) families. An accurate reconstruction of ancestral ncRNAs must use this structural signal. However, the inference of ancestors of a single ncRNA family with a single consensus structure may bias the results towards sequences with high affinity to this structure, which are far from the true ancestors.

Methods

In this paper, we introduce achARNement, a maximum parsimony approach that, given two alignments of homologous ncRNA families with consensus secondary structures and a phylogenetic tree, simultaneously calculates ancestral RNA sequences for these two families.

Results

We test our methodology on simulated data sets, and show that achARNement outperforms classical maximum parsimony approaches in terms of accuracy, but also reduces by several orders of magnitude the number of candidate sequences. To conclude this study, we apply our algorithms on the Glm clan and the FinP-traJ clan from the Rfam database.

Conclusions

Our results show that our methods reconstruct small sets of high-quality candidate ancestors with better agreement to the two target structures than with classical approaches. Our program is freely available at: http://csb.cs.mcgill.ca/acharnement.

     

DTL reconciliation repair

Background

Maximum parsimony phylogenetic tree reconciliation is an important technique for reconstructing the evolutionary histories of hosts and parasites, genes and species, and other interdependent pairs. Since the problem of finding temporally feasible maximum parsimony reconciliations is NP-complete, current methods use either exact algorithms with exponential worst-case running time or heuristics that do not guarantee optimal solutions.

Results

We offer an efficient new approach that begins with a potentially infeasible maximum parsimony reconciliation and iteratively “repairs” it until it becomes temporally feasible.

Conclusions

In a non-trivial number of cases, this approach finds solutions that are better than those found by the widely-used Jane heuristic.

     

A modified yeast three-hybrid system enabling both positive and negative selections

Objective

To increase the reporter repertoire of the yeast three-hybrid system and introduce the option of negative selection.

Results

Two new versions of the yeast three-hybrid system were made by modifying the MS2 coat RNA-binding protein and fusing it to the Gal4 DNA-binding protein. This allows the use of Gal4 inducible reporters to measure RNA–protein interactions. We introduced two mutations, V29I and N55K into the tandem MS2 dimer and an 11 amino acid deletion to increase RNA–protein affinity and inhibit capsid formation. Introduction of these constructs into the yeast strains MaV204K and PJ69-2A (which contain more reporters than the conventional yeast three-hybrid strains L40-coat and YBZ-1) allows RNA–protein binding interactions with a wide range of affinities to be detected using histidine auxotrophy, and negative selection with 5-fluoroorotic acid.

Conclusion

This yeast three-hybrid system has advantages over previous versions as demonstrated by the increased dynamic range of detectable binding interactions using yeast survival assays and colony forming assays with multiple reporters using known RNA–protein interactions.

     

Identification of trunk mutations in gastric carcinoma: a case study

Background

Intratumor heterogeneity (ITH) poses an urgent challenge for cancer precision medicine because it can cause drug resistance against cancer target therapy and immunotherapy. The search for trunk mutations that are present in all cancer cells is therefore critical for each patient.

Case presentation

In this study, we aimed to evaluate the efficiency of multiregional sequencing for the identification of trunk mutations present in all regions of a tumor as a case study. We applied multiregional whole-exome sequencing (WES) to investigate the genetic heterogeneity and homogeneity of a case of gastric carcinoma. Approximately 83% of common missense mutations present in two samples and approximately 89% of common missense mutations present in three samples were trunk mutations. Notably, trunk mutations appeared to have higher variant allele frequencies (VAFs) than non-trunk mutations.

Conclusions

Our results indicate that small-scale multiregional sampling and subsequent screening of low VAF somatic mutations might be a cost-effective strategy for identifying the majority of trunk mutations in gastric carcinoma.

     

How the chemical features of molecules may have addressed the settlement of metabolic steps

Introduction

While the evolutionary adaptation of enzymes to their own substrates is a well assessed and rationalized field, how molecules have been originally selected in order to initiate and assemble convenient metabolic pathways is a fascinating, but still debated argument.

Objectives

Aim of the present study is to give a rationale for the preferential selection of specific molecules to generate metabolic pathways.

Methods

The comparison of structural features of molecules, through an inductive methodological approach, offer a reading key to cautiously propose a determining factor for their metabolic recruitment.

Results

Starting with some commonplaces occurring in the structural representation of relevant carbohydrates, such as glucose, fructose and ribose, arguments are presented in associating stable structural determinants of these molecules and their peculiar occurrence in metabolic pathways.

Conclusions

Among other possible factors, the reliability of the structural asset of a molecule may be relevant or its selection among structurally and, a priori, functionally similar molecules.

     

MetaBinG2: a fast and accurate metagenomic sequence classification system for samples with many unknown organisms

Background

Many methods have been developed for metagenomic sequence classification, and most of them depend heavily on genome sequences of the known organisms. A large portion of sequencing sequences may be classified as unknown, which greatly impairs our understanding of the whole sample.

Result

Here we present MetaBinG2, a fast method for metagenomic sequence classification, especially for samples with a large number of unknown organisms. MetaBinG2 is based on sequence composition, and uses GPUs to accelerate its speed. A million 100 bp Illumina sequences can be classified in about 1 min on a computer with one GPU card. We evaluated MetaBinG2 by comparing it to multiple popular existing methods. We then applied MetaBinG2 to the dataset of MetaSUB Inter-City Challenge provided by CAMDA data analysis contest and compared community composition structures for environmental samples from different public places across cities.

Conclusion

Compared to existing methods, MetaBinG2 is fast and accurate, especially for those samples with significant proportions of unknown organisms.

Reviewers

This article was reviewed by Drs. Eran Elhaik, Nicolas Rascovan, and Serghei Mangul.

     

Analysis of chicken anemia virus genome: evidence of intersubtype recombination

Background

Chicken anemia virus (CAV) is the causative agent of chicken infectious anemia. CAV putative intergenotypic recombinants have been reported previously. This fact is based on the previous classification of CAV sequences into three genotypes. However, it is unknown whether intersubtype recombination occurs between the recently reported four CAV genotypes and five subtypes of genome sequences.

Results

Phylogenetic analysis, together with a variety of computational recombination detection algorithms, was used to investigate CAV approximately full genomes. Statistically significant evidence of intersubtype recombination was detected in the parent-like and two putative CAV recombinant sequences. This event was shown to occur between CAV subgroup A1 and A2 sequences in the phylogenetic trees.

Conclusions

We revealed that intersubtype recombination in CAV genome sequences played a role in generating genetic diversity within the natural population of CAV.

     

Evaluation of several lightweight stochastic context-free grammars for RNA secondary structure prediction

Background

RNA secondary structure prediction methods based on probabilistic modeling can be developed using stochastic context-free grammars (SCFGs). Such methods can readily combine different sources of information that can be expressed probabilistically, such as an evolutionary model of comparative RNA sequence analysis and a biophysical model of structure plausibility. However, the number of free parameters in an integrated model for consensus RNA structure prediction can become untenable if the underlying SCFG design is too complex. Thus a key question is, what small, simple SCFG designs perform best for RNA secondary structure prediction?

Results

Nine different small SCFGs were implemented to explore the tradeoffs between model complexity and prediction accuracy. Each model was tested for single sequence structure prediction accuracy on a benchmark set of RNA secondary structures.

Conclusions

Four SCFG designs had prediction accuracies near the performance of current energy minimization programs. One of these designs, introduced by Knudsen and Hein in their PFOLD algorithm, has only 21 free parameters and is significantly simpler than the others.