Assembly Software Research Articles

Helsinki VideoMEG Project: Augmenting magnetoencephalography with synchronized video recordings

The primary goal of the Helsinki VideoMEG Project is to enable magnetoencephalography (MEG) practitioners to record and analyze the video of the subject during an MEG experiment jointly with the MEG data. The project provides:•Hardware assembly instructions and software for setting up video and audio recordings of the participant synchronized to MEG data acquisition.•Basic software tools for analyzing video and audio together with the MEG data.The resulting setup allows reliable recording of video and audio from the subject in various real-world usage scenarios. The Helsinki VideoMEG Project allowed successful establishment of video-MEG facilities in four different MEG laboratories in Finland, Sweden and the United States.

Transcriptome sequencing reveals high isoform diversity in the ant Formica exsecta.

Transcriptome resources for social insects have the potential to provide new insight into polyphenism, i.e., how divergent phenotypes arise from the same genome. Here we present a transcriptome based on paired-end RNA sequencing data for the ant Formica exsecta (Formicidae, Hymenoptera). The RNA sequencing libraries were constructed from samples of several life stages of both sexes and female castes of queens and workers, in order to maximize representation of expressed genes. We first compare the performance of common assembly and scaffolding software (Trinity, Velvet-Oases, and SOAPdenovo-trans), in producing de novo assemblies. Second, we annotate the resulting expressed contigs to the currently published genomes of ants, and other insects, including the honeybee, to filter genes that have annotation evidence of being true genes. Our pipeline resulted in a final assembly of altogether 39,262 mRNA transcripts, with an average coverage of >300X, belonging to 17,496 unique genes with annotation in the related ant species. From these genes, 536 genes were unique to one caste or sex only, highlighting the importance of comprehensive sampling. Our final assembly also showed expression of several splice variants in 6,975 genes, and we show that accounting for splice variants affects the outcome of downstream analyses such as gene ontologies. Our transcriptome provides an outstanding resource for future genetic studies on F. exsecta and other ant species, and the presented transcriptome assembly can be adapted to any non-model species that has genomic resources available from a related taxon.

RNA‐Seq of Agropyron mongolicum Keng in response to drought stress

AbstractAgropyron mongolicum Keng is the predominant perennial forage crop in the steppes of the Inner Mongolian Plateau, China. It is known to be tolerant to drought, low temperature and salinity, and is also a valuable resource for wheat improvement. However, the corresponding molecular mechanism underlying the adaptation to water shortage in A. mongolicum remains uncharacterized. In this study, the molecular responses of A. mongolicum were compared during drought stress. Two libraries prepared from control and drought‐stressed plants were sequenced via the Illumina HiSeq 2500 platform. A total of 44733 unigenes were obtained using de novo assembly software, of which 31944 had BLAST hits in the NCBI non‐redundant database. 22659 unigenes were assigned to gene ontology (GO), 8451 unigenes to clusters of orthologous groups (COG), 5178 unigenes to kyoto encyclopedia of genes and genomes (KEGG), and 20 850 unigenes to Swissprot database. There were significant differences in expression amongst a total of 2171 unigenes under drought‐stress conditions. These DEGs included 1300 GO terms and 83 KEGG pathways. Using next generation sequencing (NGS) technology, the transcriptome of A. mongolicum was studied; this provided important information for further gene expression and functional genomics studies, confirmed our results on the molecular mechanisms of drought response in A. mongolicum, and provides a useful resource for breeding of drought‐resistant crops.

Critical Assessment of Metagenome Interpretation\u2014a benchmark of metagenomics software

The Critical Assessment of Metagenome Interpretation (CAMI) community initiative presents results from its first challenge, a rigorous benchmarking of software for metagenome assembly, binning and taxonomic profiling.

Comparative genomic analysis of Mycobacterium neoaurum MN2 and MN4 substrate and product tolerance.

The microbial bioconversion of sterols can afford valuable steroid precursors, such as 4-androstene-3,17-dione (AD) and androsta-1,4-diene-3,17-dione (ADD). The Mycobacterium neoaurum MN4 mutant strain can produce AD in high yield and can tolerate a higher concentration of the substrate phytosterol than the parent strain M. neoaurum MN2. In order to further investigate the mechanisms underlying the enhanced substrate and product tolerance, we performed a genomic analysis of the MN2 and MN4 strains. The genomes were sequenced using a high-throughput approach and analyzed using software for genome assembly, gene prediction and functional annotation, KEGG (Kyoto Encyclopedia of Genes and Genomes) annotation, COG (cluster of orthologous) group cluster analysis, GO cluster analysis, and SNP detection and annotation. Based on comparative genomics, 184 mutations were identified in MN4, the average variant rate of 1 variant every 27,249 bases, with a TS/TV value of 0.5877 and missense mutations in one key sterol degradation genes (ChoM1) and four side chain degradation genes that encode enzymes catalysing β-oxidation. The results suggest the high AD yield might be due to mutation of ChoM and genes encoding FadE, FadB and FadA β-oxidation enzymes. This study provides a theoretical basis for further functional genomics analysis and heterologous production of M. neoaurum MN2 secondary metabolites.

OMSim: a simulator for optical map data.

MotivationThe Bionano Genomics platform allows for the optical detection of short sequence patterns in very long DNA molecules (up to 2.5 Mbp). Molecules with overlapping patterns can be assembled to generate a consensus optical map of the entire genome. In turn, these optical maps can be used to validate or improve de novo genome assembly projects or to detect large-scale structural variation in genomes. Simulated optical map data can assist in the development and benchmarking of tools that operate on those data, such as alignment and assembly software. Additionally, it can help to optimize the experimental setup for a genome of interest. Such a simulator is currently not available.ResultsWe have developed a simulator, OMSim, that produces synthetic optical map data that mimics real Bionano Genomics data. These simulated data have been tested for compatibility with the Bionano Genomics Irys software system and the Irys-scaffolding scripts. OMSim is capable of handling very large genomes (over 30 Gbp) with high throughput and low memory requirements.Availability and implementationThe Python simulation tool and a cross-platform graphical user interface are available as open source software under the GNU GPL v2 license (http://www.bioinformatics.intec.ugent.be/omsim).Supplementary information Supplementary data are available at Bioinformatics online.

Asset condition assessment of municipal drinking water, wastewater and stormwater systems – Challenges and directions forward

The technical bases for the development of asset management plans (AMPs) for municipal water infrastructure are rapidly evolving, with regard to both available software and strategies for data assembly. The challenges and assistance available from an array of AMPs in the technical literature are reviewed, including the capabilities and suitability of some data acquisition technologies. Specific strategies which are important for municipalities while developing effective AMPs are also described. The review concludes there is extensive literature on ‘how’ successful AMPs can be implemented for drinking water and wastewater systems in Canada, but guidance for stormwater is lacking. Guidance on assessing the condition of existing assets, how to define and maintain acceptable levels of service, and the availability of software tools to assist in asset management exist for the former, but not the latter. Moving forward, opportunities to better link both grey (traditional) and green (non-traditional) stormwater, wastewater and drinking water infrastructure are described.

Development of a Feature and Template-Assisted Assembler and Application to the Analysis of a Foot-and-Mouth Disease Virus Genotyping Microarray.

Several RT-PCR and genome sequencing strategies exist for the resolution of Foot-and-Mouth Disease virus (FMDV). While these approaches are relatively straightforward, they can be vulnerable to failure due to the unpredictable nature of FMDV genome sequence variations. Sequence independent single primer amplification (SISPA) followed by genotyping microarray offers an attractive unbiased approach to FMDV characterization. Here we describe a custom FMDV microarray and a companion feature and template-assisted assembler software (FAT-assembler) capable of resolving virus genome sequence using a moderate number of conserved microarray features. The results demonstrate that this approach may be used to rapidly characterize naturally occurring FMDV as well as an engineered chimeric strain of FMDV. The FAT-assembler, while applied to resolving FMDV genomes, represents a new bioinformatics approach that should be broadly applicable to interpreting microarray genotyping data for other viruses or target organisms.

Genome Sequencing of Steroid-Producing Bacteria with Illumina Technology.

Illumina technology is widely used for bacterial whole-genome sequencing due to its simplicity, cheapness, reliability, and abundant software for manipulation with raw data. Illumina technology belongs to a second generation of whole genome sequencing that yields great amount of short reads for genome regions. Genomic DNA is fragmented to short pieces. DNA fragments are amplified for signal increasing, and are read using sequencing-by-synthesis. Millions of short reads up to 100-300 bp in length are assembled in continuous sequences. Mate-pair technology allows resolving a long repeat.Here, we describe the principles of standard and mate-pair library preparation from DNA samples, library quality control, sequencing with MiSeq instrument and following data bioinformatics treatment. Software for genome assembly and completion are listed that assemble, map, annotate, visualize, edit and allow doing other manipulations with genomic sequences. The whole genomes sequencing of the steroid-producing Actinobacteria using these protocols is exemplified.

Comparison of bacterial genome assembly software for MinION data and their applicability to medical microbiology

Translating the Oxford Nanopore MinION sequencing technology into medical microbiology requires on-going analysis that keeps pace with technological improvements to the instrument and release of associated analysis software. Here, we use a multidrug-resistant Enterobacter kobei isolate as a model organism to compare open source software for the assembly of genome data, and relate this to the time taken to generate actionable information. Three software tools (PBcR, Canu and miniasm) were used to assemble MinION data and a fourth (SPAdes) was used to combine MinION and Illumina data to produce a hybrid assembly. All four had a similar number of contigs and were more contiguous than the assembly using Illumina data alone, with SPAdes producing a single chromosomal contig. Evaluation of the four assemblies to represent the genome structure revealed a single large inversion in the SPAdes assembly, which also incorrectly integrated a plasmid into the chromosomal contig. Almost 50 %, 80 % and 90 % of MinION pass reads were generated in the first 6, 9 and 12 h, respectively. Using data from the first 6 h alone led to a less accurate, fragmented assembly, but data from the first 9 or 12 h generated similar assemblies to that from 48 h sequencing. Assemblies were generated in 2 h using Canu, indicating that going from isolate to assembled data is possible in less than 48 h. MinION data identified that genes responsible for resistance were carried by two plasmids encoding resistance to carbapenem and to sulphonamides, rifampicin and aminoglycosides, respectively.

De novo transcriptome sequencing facilitates genomic resource generation in Tinospora cordifolia.

Tinospora cordifolia is known for its medicinal properties owing to the presence of useful constituents such as terpenes, glycosides, steroids, alkaloids, and flavonoids belonging to secondary metabolism origin. However, there is little information available pertaining to critical genomic elements (ESTs, molecular markers) necessary for judicious exploitation of its germplasm. We employed 454 GS-FLX pyrosequencing of entire transcripts and altogether ∼25K assembled transcripts or Expressed sequence tags (ESTs) were identified. As the interest in T. cordifolia is primarily due to its secondary metabolite constituents, the ESTs pertaining to terpenoids biosynthetic pathway were identified in the present study. Additionally, several ESTs were assigned to different transcription factor families. To validate our transcripts dataset, the novel EST-SSR markers were generated to assess the genetic diversity among germplasm of T. cordifolia. These EST-SSR markers were found to be polymorphic and the dendrogram based on dice similarity index revealed three distinct clustering of accessions. The present study demonstrates effectiveness in using both NEWBLER and MIRA sequence read assembler software for enriching transcript-dataset and thus enables better exploitation of EST resources for mining candidate genes and designing molecular markers.

MO-FG-202-04: Gantry-Resolved Linac QA for VMAT: A Comprehensive and Efficient System Using An Electronic Portal Imaging Device

Purpose: To automate gantry-resolved linear accelerator (linac) quality assurance (QA) for volumetric modulated arc therapy (VMAT) using an electronic portal imaging device (EPID). Methods: A QA system for VMAT was developed that uses an EPID, frame-grabber assembly and in-house developed image processing software. The system relies solely on the analysis of EPID image frames acquired without the presence of a phantom. Images were acquired at 8.41 frames per second using a frame grabber and ancillary acquisition computer. Each image frame was tagged with a gantry angle from the linac's on-board gantry angle encoder. Arc-dynamic QA plans were designed to assess the performance of each individual linac component during VMAT. By analysing each image frame acquired during the QA deliveries the following eight machine performance characteristics were measured as a function of gantry angle: MLC positional accuracy, MLC speed constancy, MLC acceleration constancy, MLC-gantry synchronisation, beam profile constancy, dose rate constancy, gantry speed constancy, dose-gantry angle synchronisation and mechanical sag. All tests were performed on a Varian iX linear accelerator equipped with a 120 leaf Millennium MLC and an aS1000 EPID (Varian Medical Systems, Palo Alto, CA, USA). Results: Machine performance parameters were measured as a function of gantry angle using EPID imaging and compared to machine log files and the treatment plan. Data acquisition is currently underway at 3 centres, incorporating 7 treatment units, at 2 weekly measurement intervals. Conclusion: The proposed system can be applied for streamlined linac QA and commissioning for VMAT. The set of test plans developed can be used to assess the performance of each individual components of the treatment machine during VMAT deliveries as a function of gantry angle. The methodology does not require the setup of any additional phantom or measurement equipment and the analysis is fully automated to allow for regular routine testing.

Analysis of the magnetic field for MR damper working in shear mode

Seismic control, restorative and mechanical vibrations speak about a wide range of research that is emerging rapidly. Magneto-rheological dampers have as of late turned into an object of concentrated studies because of their interesting physical features and their ability to control damping in mechanical systems. In this paper an MR damper is analysed as a 2D model to be axisymmetric and finite element analysis using ANSYS is done to simulate the distribution of the magnetic field in the damper for a given current. Then the design for assembly (DFA) software tool is used for estimating cost of assembly of the final product.

De novo assembly and characterization of farmed blue fox (Alopex lagopus) global transcriptome using Illumina paired-end sequencing.

The blue fox (Alopex lagopus), a coat-color variant of the Arctic fox, is a domesticated fur-bearing mammal. In the present study, transcriptome data generated from a pool of nine different tissues were obtained with Illumina HiSeq2500 paired-end sequencing technology. After filtering from raw reads, 32,358,290 clean reads were assembled into 161,269 transcripts and 97,252 unigenes by the Trinity fragment assembly software. Of the assembled unigenes, 37,967 were annotated in the National Center for Biotechnology Information (NCBI) Non-Redundant (NR) protein database and 26,264 in the Swiss-Prot database. Among the annotated unigenes, 24,839 and 24,267 were assigned using the Gene Ontology (GO) and euKaryotic Orthologous Groups (KOG) databases, respectively. Altogether, 17,057 unigenes were mapped onto 227 pathways using the Kyoto Encyclopedia of Genes and Genomes database. In addition, 6394 simple sequence repeats were identified by examining 12,965 unigenes (>1 kb), which could contribute to the development of molecular markers. This study generated transcriptome data for the blue fox that will promote further progress in expression profiling studies, and provide a good annotation basis for genomic studies.

Genome sequencing and analysis of a granulovirus isolated from the Asiatic rice leafroller, Cnaphalocrocis medinalis.

The complete genome of Cnaphalocrocis medinalis granulovirus (CnmeGV) from a serious migratory rice pest, Cnaphalocrocis medinalis (Lepidoptera: Pyralidae), was sequenced using the Roche 454 Genome Sequencer FLX system (GS FLX) with shotgun strategy and assembled by Roche GS De Novo assembler software. Its circular double-stranded genome is 111,246 bp in size with a high A+T content of 64.8% and codes for 118 putative open reading frames (ORFs). It contains 37 conserved baculovirus core ORFs, 13 unique ORFs, 26 ORFs that were found in all Lepidoptera baculoviruses and 42 common ORFs. The analysis of nucleotide sequence repeats revealed that the CnmeGV genome differs from the rest of sequenced GVs by a 23 kb and a 17kb gene block inversions, and does not contain any typical homologous region (hr) except for a region of non-hr-like sequence. Chitinase and cathepsin genes, which are reported to have major roles in the liquefaction of the hosts, were not found in the CnmeGV genome, which explains why CnmeGV infected insects do not show the phenotype of typical liquefaction. Phylogenetic analysis, based on the 37 core baculovirus genes, indicates that CnmeGV is closely related to Adoxophyes orana granulovirus. The genome analysis would contribute to the functional research of CnmeGV, and would benefit to the utilization of CnmeGV as pest control reagent for rice production.

Tiling Assembly: a new tool for reference annotation-independent transcript assembly and novel gene identification by RNA-sequencing.

Annotation of the rice (Oryza sativa) genome has evolved significantly since release of its draft sequence, but it is far from complete. Several published transcript assembly programmes were tested on RNA-sequencing (RNA-seq) data to determine their effectiveness in identifying novel genes to improve the rice genome annotation. Cufflinks, a popular assembly software, did not identify all transcripts suggested by the RNA-seq data. Other assembly software was CPU intensive, lacked documentation, or lacked software updates. To overcome these shortcomings, a heuristic ab initio transcript assembly algorithm, Tiling Assembly, was developed to identify genes based on short read and junction alignment. Tiling Assembly was compared with Cufflinks to evaluate its gene-finding capabilities. Additionally, a pipeline was developed to eliminate false-positive gene identification due to noise or repetitive regions in the genome. By combining Tiling Assembly and Cufflinks, 767 unannotated genes were identified in the rice genome, demonstrating that combining both programmes proved highly efficient for novel gene identification. We also demonstrated that Tiling Assembly can accurately determine transcription start sites by comparing the Tiling Assembly genes with their corresponding full-length cDNA. We applied our pipeline to additional organisms and identified numerous unannotated genes, demonstrating that Tiling Assembly is an organism-independent tool for genome annotation.

The Streptomyces leeuwenhoekii genome: de novo sequencing and assembly in single contigs of the chromosome, circular plasmid pSLE1 and linear plasmid pSLE2.

BackgroundNext Generation DNA Sequencing (NGS) and genome mining of actinomycetes and other microorganisms is currently one of the most promising strategies for the discovery of novel bioactive natural products, potentially revealing novel chemistry and enzymology involved in their biosynthesis. This approach also allows rapid insights into the biosynthetic potential of microorganisms isolated from unexploited habitats and ecosystems, which in many cases may prove difficult to culture and manipulate in the laboratory. Streptomyces leeuwenhoekii (formerly Streptomyces sp. strain C34) was isolated from the hyper-arid high-altitude Atacama Desert in Chile and shown to produce novel polyketide antibiotics.ResultsHere we present the de novo sequencing of the S. leeuwenhoekii linear chromosome (8 Mb) and two extrachromosomal replicons, the circular pSLE1 (86 kb) and the linear pSLE2 (132 kb), all in single contigs, obtained by combining Pacific Biosciences SMRT (PacBio) and Illumina MiSeq technologies. We identified the biosynthetic gene clusters for chaxamycin, chaxalactin, hygromycin A and desferrioxamine E, metabolites all previously shown to be produced by this strain (J Nat Prod, 2011, 74:1965) and an additional 31 putative gene clusters for specialised metabolites. As well as gene clusters for polyketides and non-ribosomal peptides, we also identified three gene clusters encoding novel lasso-peptides.ConclusionsThe S. leeuwenhoekii genome contains 35 gene clusters apparently encoding the biosynthesis of specialised metabolites, most of them completely novel and uncharacterised. This project has served to evaluate the current state of NGS for efficient and effective genome mining of high GC actinomycetes. The PacBio technology now permits the assembly of actinomycete replicons into single contigs with >99 % accuracy. The assembled Illumina sequence permitted not only the correction of omissions found in GC homopolymers in the PacBio assembly (exacerbated by the high GC content of actinomycete DNA) but it also allowed us to obtain the sequences of the termini of the chromosome and of a linear plasmid that were not assembled by PacBio. We propose an experimental pipeline that uses the Illumina assembled contigs, in addition to just the reads, to complement the current limitations of the PacBio sequencing technology and assembly software.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-1652-8) contains supplementary material, which is available to authorized users.

Bioinformatic Suggestions on MiSeq-Based Microbial Community Analysis.

Recent sequencing technology development has revolutionized fields of microbial ecology. MiSeq-based microbial community analysis allows us to sequence more than a few hundred samples at a time, which is far more cost-effective than pyrosequencing. The approach, however, has not been preferably used owing to computational difficulties of processing huge amounts of data as well as known Illumina-derived artefact problems with amplicon sequencing. The choice of assembly software to take advantage of paired-end sequencing and methods to remove Illumina artefacts sequences are discussed. The protocol we suggest not only removed erroneous reads, but also dramatically reduced computational workload, which allows even a typical desktop computer to process a huge amount of sequence data generated with Illumina sequencers. We also developed a Web interface (http://biotech.jejunu.ac.kr/ ~abl/16s/) that allows users to conduct fastq-merging and mothur batch creation. The study presented here should provide technical advantages and supports in applying MiSeq-based microbial community analysis.

QuorUM: An Error Corrector for Illumina Reads

MotivationIllumina Sequencing data can provide high coverage of a genome by relatively short (most often 100 bp to 150 bp) reads at a low cost. Even with low (advertised 1%) error rate, 100 × coverage Illumina data on average has an error in some read at every base in the genome. These errors make handling the data more complicated because they result in a large number of low-count erroneous k-mers in the reads. However, there is enough information in the reads to correct most of the sequencing errors, thus making subsequent use of the data (e.g. for mapping or assembly) easier. Here we use the term “error correction” to denote the reduction in errors due to both changes in individual bases and trimming of unusable sequence. We developed an error correction software called QuorUM. QuorUM is mainly aimed at error correcting Illumina reads for subsequent assembly. It is designed around the novel idea of minimizing the number of distinct erroneous k-mers in the output reads and preserving the most true k-mers, and we introduce a composite statistic π that measures how successful we are at achieving this dual goal. We evaluate the performance of QuorUM by correcting actual Illumina reads from genomes for which a reference assembly is available.ResultsWe produce trimmed and error-corrected reads that result in assemblies with longer contigs and fewer errors. We compared QuorUM against several published error correctors and found that it is the best performer in most metrics we use. QuorUM is efficiently implemented making use of current multi-core computing architectures and it is suitable for large data sets (1 billion bases checked and corrected per day per core). We also demonstrate that a third-party assembler (SOAPdenovo) benefits significantly from using QuorUM error-corrected reads. QuorUM error corrected reads result in a factor of 1.1 to 4 improvement in N50 contig size compared to using the original reads with SOAPdenovo for the data sets investigated.AvailabilityQuorUM is distributed as an independent software package and as a module of the MaSuRCA assembly software. Both are available under the GPL open source license at http://www.genome.umd.edu.Contact gmarcais@umd.edu.

An improved fruit transcriptome and the identification of the candidate genes involved in fruit abscission induced by carbohydrate stress in litchi.

Massive young fruit abscission usually causes low and unstable yield in litchi (Litchi chinensis Sonn.), an important fruit crop cultivated in tropical and subtropical areas. However, the molecular mechanism of fruit drop has not been fully characterized. This study aimed at identification of molecular components involved in fruitlet abscission in litchi, for which reference genome is not available at present. An improved de novo transcriptome assembly was firstly achieved by using an optimized assembly software, Trinity. Using improved transcriptome assembly as reference, digital transcript abundance (DTA) profiling was performed to screen and identify candidate genes involved in fruit abscission induced by girdling plus defoliation (GPD), a treatment significantly decreased the soluble sugar contents causing carbohydrate stress to fruit. Our results showed that the increasing fruit abscission rate after GPD treatment was associated with higher ethylene production and lower glucose levels in fruit. A total of 2,771 differentially expressed genes were identified as GPD-responsive genes, 857 of which were defined by GO and KEGG enrichment analyses as the candidate genes involved in fruit abscission process. These genes were involved in diverse metabolic processes and pathways, including carbohydrate metabolism, plant hormone synthesis, and signaling, transcription factor activity and cell wall modification that were rapidly induced in the early stages (within 2 days after treatment). qRT-PCR was used to explore the expression pattern of 15 selected candidate genes in the abscission zone, pericarp, and seed, which confirmed the accuracy of our DTA data. More detailed information for different functional categories was also analyzed. This study profiled the gene expression related to fruit abscission induced by carbohydrate stress at whole transcriptome level and thus provided a better understanding of the regulatory mechanism of young fruit abscission in litchi.