MiteFinderII: a novel tool to identify miniature inverted-repeat transposable elements hidden in eukaryotic genomes

Background

Miniature inverted-repeat transposable element (MITE) is a type of class II non-autonomous transposable element playing a crucial role in the process of evolution in biology. There is an urgent need to develop bioinformatics tools to effectively identify MITEs on a whole genome-wide scale. However, most of currently existing tools suffer from low ability to deal with large eukaryotic genomes.

Methods

In this paper, we proposed a novel tool MiteFinderII, which was adapted from our previous algorithm MiteFinder, to efficiently detect MITEs from genomics sequences. It has six major steps: (1) build K-mer Index and search for inverted repeats; (2) filtration of inverted repeats with low complexity; (3) merger of inverted repeats; (4) filtration of candidates with low score; (5) selection of final MITE sequences; (6) selection of representative sequences.

Results

To test the performance, MiteFinderII and three other existing algorithms were applied to identify MITEs on the whole genome of oryza sativa. Results suggest that MiteFinderII outperforms existing popular tools in terms of both specificity and recall. Additionally, it is much faster and more memory-efficient than other tools in the detection.

Conclusion

MiteFinderII is an accurate and effective tool to detect MITEs hidden in eukaryotic genomes. The source code is freely accessible at the website: https://github.com/screamer/miteFinder.

     

PSAT: A web tool to compare genomic neighborhoods of multiple prokaryotic genomes

Background  

The conservation of gene order among prokaryotic genomes can provide valuable insight into gene function, protein interactions, or events by which genomes have evolved. Although some tools are available for visualizing and comparing the order of genes between genomes of study, few support an efficient and organized analysis between large numbers of genomes. The Prokaryotic Sequence homology Analysis Tool (PSAT) is a web tool for comparing gene neighborhoods among multiple prokaryotic genomes.     

Genomic suppression subtractive hybridization as a tool to identify differences in mycorrhizal fungal genomes

Understanding the ecology of methanogens in natural and engineered environments is a prerequisite to predicting or managing methane emissions. In this study, a novel high-throughput fingerprint method was developed for determining methanogen diversity and relative abundance within environmental samples. The method described here, designated amplicon length heterogeneity PCR of the mcrA gene (LH-mcrA), is based on the natural length variation in the mcrA gene. The mcrA gene encodes the alpha-subunit of the methyl-coenzyme M reductase, which is involved in the terminal step of methane production by methanogens. The methanogenic communities from stored swine and dairy manures were distinct from each other. To validate the method, methanogenic communities in a plug flow-type bioreactor (PFBR) treating swine manure were characterized using LH-mcrA method and correlated to mcrA gene clone libraries. The diversity and relative abundance of the methanogenic groups were assessed. Methanobrevibacter, Methanosarcinaceae, Methanoculleus, Methanogenium, Methanocorpusculum and one unidentified group were assigned to particular LH-mcrA amplicons. Particular phylotypes related to Methanoculleus were predominant in the last compartment of the PFBR where the bulk of methane was produced. LH-mcrA method was found to be a reliable, fast and cost-effective alternative for diversity assessment of methanogenic communities in microbial systems.     

Intraspecific phylogenetic congruence among multiple symbiont genomes

Eukaryotes often form intimate endosymbioses with prokaryotic organisms. Cases in which these symbionts are transmitted cytoplasmically to host progeny create the potential for co-speciation or congruent evolution among the distinct genomes of these partners. If symbionts do not move horizontally between different eukaryotic hosts, strict phylogenetic congruence of their genomes is predicted and should extend to relationships within a single host species. Conversely, even rare 'host shifts' among closely related lineages should yield conflicting tree topologies at the intraspecific level. Here, we investigate the historical associations among four symbiotic genomes residing within an aphid host: the mitochondrial DNA of Uroleucon ambrosiae aphids, the bacterial chromosome of their Buchnera bacterial endosymbionts, and two plasmids associated with Buchnera. DNA sequence polymorphisms provided a significant phylogenetic signal and no homoplasy for each data set, yielding completely and significantly congruent phylogenies for these four genomes and no evidence of horizontal transmission. This study thus provides the first evidence for strictly vertical transmission and 'co-speciation' of symbiotic organisms at the intraspecific level, and represents the lowest phylogenetic level at which such coevolution has been demonstrated. These results may reflect the obligate nature of this intimate mutualism and indicate opportunities for adaptive coevolution among linked symbiont genomes.     

Java GUI for InterProScan (JIPS): A tool to help process multiple InterProScans and perform ortholog analysis

Background  

Recent, rapid growth in the quantity of available genomic data has generated many protein sequences that are not yet biochemically classified. Thus, the prediction of biochemical function based on structural motifs is an important task in post-genomic analysis. The InterPro databases are a major resource for protein function information. For optimal results, these databases should be searched at regular intervals, since they are frequently updated.     

MGView: an alignment and visualization tool to enhance gap closure of microbial genomes

Gap closure is a challenging phase in microbial random shotgun genome sequencing projects, particularly since genome assemblies are often complicated by the presence of repeat elements, insertion sequences and other similar factors that contribute to sequence misassemblies. While it is well recognized that the conservation of genetic information between microbial genomes, combined with the exponential increase in available microbial sequences, can be exploited to increase the efficiency of gap closure, we lack the computational tools to aid in this process. We describe here a new tool, MGView, which was developed to create a graphical depiction of the alignment of a set of microbial contigs against a completed microbial genome. The results of our assembly of the Staphylococcus aureus RF122 genome show that MGView enables a considerable reduction in time and economic cost associated with closure. Together, the results also show that the application of MGView not only enables a reduction in fold-coverage requirements of the random shotgun sequence phase, but also provides interesting insights into differences in gene content and organization between finished and unfinished microbial genomes.     

Interactions among multiple genomes: tsetse, its symbionts and trypanosomes

Insect-borne diseases exact a high public health burden and have a devastating impact on livestock and agriculture. To date, control has proved to be exceedingly difficult. One such disease that has plagued sub-Saharan Africa is caused by the protozoan African trypanosomes (Trypanosoma species) and transmitted by tsetse flies (Diptera: Glossinidae). This presentation describes the biology of the tsetse fly and its interactions with trypanosomes as well as its symbionts. Tsetse can harbor up to three distinct microbial symbionts, including two enterics (Wigglesworthia glossinidia and Sodalis glossinidius) as well as facultative Wolbachia infections, which influence host physiology. Recent investigations into the genome of the obligate symbiont Wigglesworthia have revealed characteristics indicative of its long co-evolutionary history with the tsetse host species. Comparative analysis of the commensal-like Sodalis with free-living enterics provides examples of adaptations to the host environment (physiology and ecology), reflecting genomic tailoring events during the process of transitioning into a symbiotic lifestyle. From an applied perspective, the extensive knowledge accumulated on the genomic and developmental biology of the symbionts coupled with our ability to both express foreign genes in these microbes in vitro and repopulate tsetse midguts with these engineered microbes now provides a means to interfere with the host physiological traits which contribute to vector competence promising a novel tool for disease management.     

MUTAGEN: multi-user tool for annotating genomes

SUMMARY: MUTAGEN is a free prokaryotic annotation system. It offers the advantages of genome comparison, graphical sequence browsers, search facilities and open-source for user-specific adjustments. The web-interface allows several users to access the system from standard desktop computers. The Sulfolobus acidocaldarius genome, and several plasmids and viruses have so far been analysed and annotated using MUTAGEN. AVAILABILITY: MUTAGEN is released as open-source software under GPL. The code is available for download and/or contribution at http://dac.molbio.ku.dk/bioinformatics/MUTAGEN/     

MSOAR 2.0: Incorporating tandem duplications into ortholog assignment based on genome rearrangement

Background  

Ortholog assignment is a critical and fundamental problem in comparative genomics, since orthologs are considered to be functional counterparts in different species and can be used to infer molecular functions of one species from those of other species. MSOAR is a recently developed high-throughput system for assigning one-to-one orthologs between closely related species on a genome scale. It attempts to reconstruct the evolutionary history of input genomes in terms of genome rearrangement and gene duplication events. It assumes that a gene duplication event inserts a duplicated gene into the genome of interest at a random location (i.e., the random duplication model). However, in practice, biologists believe that genes are often duplicated by tandem duplications, where a duplicated gene is located next to the original copy (i.e., the tandem duplication model).     

CHSMiner: a GUI tool to identify chromosomal homologous segments

Background  

The identification of chromosomal homologous segments (CHS) within and between genomes is essential for comparative genomics. Various processes including insertion/deletion and inversion could cause the degeneration of CHSs.     

PFAAT version 2.0: A tool for editing,annotating, and analyzing multiple sequence alignments

Background  

By virtue of their shared ancestry, homologous sequences are similar in their structure and function. Consequently, multiple sequence alignments are routinely used to identify trends that relate to function. This type of analysis is particularly productive when it is combined with structural and phylogenetic analysis.     

MMG: a probabilistic tool to identify submodules of metabolic pathways

Motivation: A fundamental task in systems biology is the identificationof groups of genes that are involved in the cellular responseto particular signals. At its simplest level, this often reducesto identifying biological quantities (mRNA abundance, enzymeconcentrations, etc.) which are differentially expressed intwo different conditions. Popular approaches involve using t-teststatistics, based on modelling the data as arising from a mixturedistribution. A common assumption of these approaches is thatthe data are independent and identically distributed; however,biological quantities are usually related through a complex(weighted) network of interactions, and often the more pertinentquestion is which subnetworks are differentially expressed,rather than which genes. Furthermore, in many interesting cases(such as high-throughput proteomics and metabolomics), onlyvery partial observations are available, resulting in the needfor efficient imputation techniques. Results: We introduce Mixture Model on Graphs (MMG), a novelprobabilistic model to identify differentially expressed submodulesof biological networks and pathways. The method can easily incorporateinformation about weights in the network, is robust againstmissing data and can be easily generalized to directed networks.We propose an efficient sampling strategy to infer posteriorprobabilities of differential expression, as well as posteriorprobabilities over the model parameters. We assess our methodon artificial data demonstrating significant improvements overstandard mixture model clustering. Analysis of our model resultson quantitative high-throughput proteomic data leads to theidentification of biologically significant subnetworks, as wellas the prediction of the expression level of a number of enzymes,some of which are then verified experimentally. Availability: MATLAB code is available from http://www.dcs.shef.ac.uk/~guido/software.html Contact: guido{at}dcs.shef.ac.uk Supplementary information: Supplementary data are availableat Bioinformatics online. Associate Editor: Jonathan Wren     

Xenopatients 2.0: Reprogramming the epigenetic landscapes of patient-derived cancer genomes

In the science-fiction thriller film Minority Report, a specialized police department called “PreCrime” apprehends criminals identified in advance based on foreknowledge provided by 3 genetically altered humans called “PreCogs”. We propose that Yamanaka stem cell technology can be similarly used to (epi)genetically reprogram tumor cells obtained directly from cancer patients and create self-evolving personalized translational platforms to foresee the evolutionary trajectory of individual tumors. This strategy yields a large stem cell population and captures the cancer genome of an affected individual, i.e., the PreCog-induced pluripotent stem (iPS) cancer cells, which are immediately available for experimental manipulation, including pharmacological screening for personalized “stemotoxic” cancer drugs. The PreCog-iPS cancer cells will re-differentiate upon orthotopic injection into the corresponding target tissues of immunodeficient mice (i.e., the PreCrime-iPS mouse avatars), and this in vivo model will run through specific cancer stages to directly explore their biological properties for drug screening, diagnosis, and personalized treatment in individual patients. The PreCog/PreCrime-iPS approach can perform sets of comparisons to directly observe changes in the cancer-iPS cell line vs. a normal iPS cell line derived from the same human genetic background. Genome editing of PreCog-iPS cells could create translational platforms to directly investigate the link between genomic expression changes and cellular malignization that is largely free from genetic and epigenetic noise and provide proof-of-principle evidence for cutting-edge “chromosome therapies” aimed against cancer aneuploidy. We might infer the epigenetic marks that correct the tumorigenic nature of the reprogrammed cancer cell population and normalize the malignant phenotype in vivo. Genetically engineered models of conditionally reprogrammable mice to transiently express the Yamanaka stemness factors following the activation of phenotypic copies of specific cancer diseases might crucially evaluate a “reprogramming cure” for cancer. A new era of xenopatients 2.0 generated via nuclear reprogramming of the epigenetic landscapes of patient-derived cancer genomes might revolutionize the current personalized translational platforms in cancer research.