Sequence Analysis Tools Research Articles

Evaluation of the Abbott RealTime HCV genotype II plus RUO (PLUS) assay with reference to core and NS5B sequencing

BackgroundHCV genotyping remains a critical tool for guiding initiation of therapy and selecting the most appropriate treatment regimen. Current commercial genotyping assays may have difficulty identifying 1a, 1b and genotype 6. ObjectiveTo evaluate the concordance for identifying 1a, 1b, and genotype 6 between two methods: the PLUS assay and core/NS5B sequencing. Study designThis study included 236 plasma and serum samples previously genotyped by core/NS5B sequencing. Of these, 25 samples were also previously tested by the Abbott RealTime HCV GT II Research Use Only (RUO) assay and yielded ambiguous results. The remaining 211 samples were routine genotype 1 (n=169) and genotype 6 (n=42). Genotypes obtained from sequence data were determined using a laboratory-developed HCV sequence analysis tool and the NCBI non-redundant database. ResultsAgreement between the PLUS assay and core/NS5B sequencing for genotype 1 samples was 95.8% (162/169), with 96% (127/132) and 95% (35/37) agreement for 1a and 1b samples respectively. PLUS results agreed with core/NS5B sequencing for 83% (35/42) of unselected genotype 6 samples, with the remaining seven “not detected” by the PLUS assay. Among the 25 samples with ambiguous GT II results, 15 were concordant by PLUS and core/NS5B sequencing, nine were not detected by PLUS, and one sample had an internal control failure. ConclusionsThe PLUS assay is an automated method that identifies 1a, 1b and genotype 6 with good agreement with gold-standard core/NS5B sequencing and can aid in the resolution of certain genotype samples with ambiguous GT II results.

StructAnalyzer - a tool for sequence versus structure similarity analysis.

In the world of RNAs and proteins, similarities at the level of primary structures of two comparable molecules usually correspond to structural similarities at the tertiary level. In other words, measures of sequence and structure similarities are in general correlated - a high value of sequence similarity imposes a high value of structural similarity. However, important exceptions that stay in contrast to this general rule can be identified. It is possible to find similar structures with very different sequences, as well as similar sequences with very different structures. In this paper, we focus our attention on the latter case and propose a tool, called StructAnalyzer, supporting analysis of relations between the sequence and structure similarities. Recognition of tertiary structure diversity of molecules with very similar primary structures may be the key for better understanding of mechanisms influencing folding of RNAs or proteins, and as a result for better understanding of their function. StructAnalyzer allows exploration and visualization of structural diversity in relation to sequence similarity. We show how this tool can be used to screen RNA structures in Protein Data Bank (PDB) for sequences with structural variants.

Introduction on Using the FastPCR Software and the Related Java Web Tools for PCR and Oligonucleotide Assembly and Analysis.

This chapter introduces the FastPCR software as an integrated tool environment for PCR primer and probe design, which predicts properties of oligonucleotides based on experimental studies of the PCR efficiency. The software provides comprehensive facilities for designing primers for most PCR applications and their combinations. These include the standard PCR as well as the multiplex, long-distance, inverse, real-time, group-specific, unique, overlap extension PCR for multi-fragments assembling cloning and loop-mediated isothermal amplification (LAMP). It also contains a built-in program to design oligonucleotide sets both for long sequence assembly by ligase chain reaction and for design of amplicons that tile across a region(s) of interest. The software calculates the melting temperature for the standard and degenerate oligonucleotides including locked nucleic acid (LNA) and other modifications. It also provides analyses for a set of primers with the prediction of oligonucleotide properties, dimer and G/C-quadruplex detection, linguistic complexity as well as a primer dilution and resuspension calculator. The program consists of various bioinformatical tools for analysis of sequences with the GC or AT skew, CG% and GA% content, and the purine-pyrimidine skew. It also analyzes the linguistic sequence complexity and performs generation of random DNA sequence as well as restriction endonucleases analysis. The program allows to find or create restriction enzyme recognition sites for coding sequences and supports the clustering of sequences. It performs efficient and complete detection of various repeat types with visual display. The FastPCR software allows the sequence file batch processing that is essential for automation. The program is available for download at http://primerdigital.com/fastpcr.html , and its online version is located at http://primerdigital.com/tools/pcr.html .

An Alarmingly High Proportion of HIV-1 Isolates Carrying Mutations Corresponding to Resistance to Antiretroviral Drugs among HIV-Positive High-Risk Groups in Central Vietnam: a Substudy of the National Sentinel Survey.

We studied the prevalence of HIV drug resistance among high-risk groups such as injecting drug users (IDUs), female sex workers (FSWs), and men having sex with men (MSM) in central Vietnam. We used HIV-positive blood samples from 2012-2013 sentinel surveillance surveys. Study subjects were screened for HIV infection by standardized screening assays, and the HIV-positive samples were further tested for HIV viral load and drug-resistance mutations (DRMs) by in-house assays. DRMs were identified using the Stanford University online sequence analysis tool. Their risk behaviors were also investigated. During the study period, 6,016 (high-risk) subjects were screened, and 97 tested positive for HIV infection (IDUs: n = 63, 3.0%; FSWs: n = 24, 0.9%; and MSM: n = 10, 1.0%). Ninety-two of the 97 samples (45 from 2012 and 47 from 2013) were available for further testing. HIV RNA was detected in 56 (60.9%) of the 92 samples, and drug resistance genotyping was successfully performed on 40 (71.4%) samples. All these isolates were subtype CRF01_AE, except for 1 (2.5%) IDU whose HIV belongs to subtype B. Thirteen individuals (32.5%) were carrying HIV with at least 1 DRM: 9 IDUs, 1 FSW, and 3 MSM. Thus, HIV seroprevalence among high-risk individuals in central Vietnam is low, but a high proportion of drug resistant HIV-1 isolates is observed in the high-risk group.

MIToS.jl: mutual information tools for protein sequence analysis in the Julia language.

MIToS is an environment for mutual information analysis and a framework for protein multiple sequence alignments (MSAs) and protein structures (PDB) management in Julia language. It integrates sequence and structural information through SIFTS, making Pfam MSAs analysis straightforward. MIToS streamlines the implementation of any measure calculated from residue contingency tables and its optimization and testing in terms of protein contact prediction. As an example, we implemented and tested a BLOSUM62-based pseudo-count strategy in mutual information analysis. The software is totally implemented in Julia and supported for Linux, OS X and Windows. It's freely available on GitHub under MIT license: http://mitos.leloir.org.ar . diegozea@gmail.com or cmb@leloir.org.ar. Supplementary data are available at Bioinformatics online.

The European Classical Swine Fever Virus Database: Blueprint for a Pathogen-Specific Sequence Database with Integrated Sequence Analysis Tools.

Molecular epidemiology has become an indispensable tool in the diagnosis of diseases and in tracing the infection routes of pathogens. Due to advances in conventional sequencing and the development of high throughput technologies, the field of sequence determination is in the process of being revolutionized. Platforms for sharing sequence information and providing standardized tools for phylogenetic analyses are becoming increasingly important. The database (DB) of the European Union (EU) and World Organisation for Animal Health (OIE) Reference Laboratory for classical swine fever offers one of the world’s largest semi-public virus-specific sequence collections combined with a module for phylogenetic analysis. The classical swine fever (CSF) DB (CSF-DB) became a valuable tool for supporting diagnosis and epidemiological investigations of this highly contagious disease in pigs with high socio-economic impacts worldwide. The DB has been re-designed and now allows for the storage and analysis of traditionally used, well established genomic regions and of larger genomic regions including complete viral genomes. We present an application example for the analysis of highly similar viral sequences obtained in an endemic disease situation and introduce the new geographic “CSF Maps” tool. The concept of this standardized and easy-to-use DB with an integrated genetic typing module is suited to serve as a blueprint for similar platforms for other human or animal viruses.

Investigating Triticeae anther gene promoter activity in transgenic Brachypodium distachyon.

In this report, we demonstrate that Brachypodium distachyon could serve as a relatively high throughput in planta functional assay system for Triticeae anther-specific gene promoters. There remains a vast gap in our knowledge of the promoter cis-acting elements responsible for the transcriptional regulation of Triticeae anther-specific genes. In an attempt to identify conserved cis-elements, 14 pollen-specific and 8 tapetum-specific Triticeae putative promoter sequences were analyzed using different promoter sequence analysis tools. Several cis-elements were found to be enriched in these sequences and their possible role in gene expression regulation in the anther is discussed. Despite the fact that potential cis-acting elements can be identified within putative promoter sequence datasets, determining whether particular promoter sequences can in fact direct proper tissue-specific and developmental gene expression still needs to be confirmed via functional assays preferably performed in closely related plants. Transgenic functional assays with Triticeae species remain challenging and Brachypodium distachyon may represent a suitable alternative. The promoters of the triticale pollen-specific genes group 3 pollen allergen (PAL3) and group 4 pollen allergen (PAL4), as well as the tapetum-specific genes chalcone synthase-like 1 (CHSL1), from wheat and cysteine-rich protein 1 (CRP1) from triticale were fused to the green fluorescent protein gene (GFP) and analyzed in transgenic Brachypodium. This report demonstrates that this model species could serve to accelerate the functional analysis of Triticeae anther-specific gene promoters.

VSEARCH: a versatile open source tool for metagenomics.

BackgroundVSEARCH is an open source and free of charge multithreaded 64-bit tool for processing and preparing metagenomics, genomics and population genomics nucleotide sequence data. It is designed as an alternative to the widely used USEARCH tool (Edgar, 2010) for which the source code is not publicly available, algorithm details are only rudimentarily described, and only a memory-confined 32-bit version is freely available for academic use.MethodsWhen searching nucleotide sequences, VSEARCH uses a fast heuristic based on words shared by the query and target sequences in order to quickly identify similar sequences, a similar strategy is probably used in USEARCH. VSEARCH then performs optimal global sequence alignment of the query against potential target sequences, using full dynamic programming instead of the seed-and-extend heuristic used by USEARCH. Pairwise alignments are computed in parallel using vectorisation and multiple threads.ResultsVSEARCH includes most commands for analysing nucleotide sequences available in USEARCH version 7 and several of those available in USEARCH version 8, including searching (exact or based on global alignment), clustering by similarity (using length pre-sorting, abundance pre-sorting or a user-defined order), chimera detection (reference-based or de novo), dereplication (full length or prefix), pairwise alignment, reverse complementation, sorting, and subsampling. VSEARCH also includes commands for FASTQ file processing, i.e., format detection, filtering, read quality statistics, and merging of paired reads. Furthermore, VSEARCH extends functionality with several new commands and improvements, including shuffling, rereplication, masking of low-complexity sequences with the well-known DUST algorithm, a choice among different similarity definitions, and FASTQ file format conversion. VSEARCH is here shown to be more accurate than USEARCH when performing searching, clustering, chimera detection and subsampling, while on a par with USEARCH for paired-ends read merging. VSEARCH is slower than USEARCH when performing clustering and chimera detection, but significantly faster when performing paired-end reads merging and dereplication. VSEARCH is available at https://github.com/torognes/vsearch under either the BSD 2-clause license or the GNU General Public License version 3.0.DiscussionVSEARCH has been shown to be a fast, accurate and full-fledged alternative to USEARCH. A free and open-source versatile tool for sequence analysis is now available to the metagenomics community.

The Papillomavirus Episteme: a major update to the papillomavirus sequence database.

The Papillomavirus Episteme (PaVE) is a database of curated papillomavirus genomic sequences, accompanied by web-based sequence analysis tools. This update describes the addition of major new features. The papillomavirus genomes within PaVE have been further annotated, and now includes the major spliced mRNA transcripts. Viral genes and transcripts can be visualized on both linear and circular genome browsers. Evolutionary relationships among PaVE reference protein sequences can be analysed using multiple sequence alignments and phylogenetic trees. To assist in viral discovery, PaVE offers a typing tool; a simplified algorithm to determine whether a newly sequenced virus is novel. PaVE also now contains an image library containing gross clinical and histopathological images of papillomavirus infected lesions. Database URL: https://pave.niaid.nih.gov/.

PseKRAAC: a flexible web server for generating pseudo K-tuple reduced amino acids composition.

Freely available on the web at http://bigdata.imu.edu.cn/psekraac CONTACTS: yczuo@imu.edu.cn or imu.hema@foxmail.com or yanglei_hmu@163.comSupplementary information: Supplementary data are available at Bioinformatics online.

The Effects of Trivalent Lanthanide Cationization on the Electron Transfer Dissociation of Acidic Fibrinopeptide B and its Analogs.

Electrospray ionization (ESI) on mixtures of acidic fibrinopeptide B and two peptide analogs with trivalent lanthanide salts generates [M + Met + H](4+), [M + Met](3+), and [M + Met -H](2+), where M = peptide and Met = metal (except radioactive promethium). These ions undergo extensive and highly efficient electron transfer dissociation (ETD) to form metallated and non-metallated c- and z-ions. All metal adducted product ions contain at least two acidic sites, which suggest attachment of the lanthanide cation at the side chains of one or more acidic residues. The three peptides undergo similar fragmentation. ETD on [M + Met + H](4+) leads to cleavage at every residue; the presence of both a metal ion and an extra proton is very effective in promoting sequence-informative fragmentation. Backbone dissociation of [M + Met](3+) is also extensive, although cleavage does not always occur between adjacent glutamic acid residues. For [M + Met - H ](2+), a more limited range of product ions form. All lanthanide metal peptide complexes display similar fragmentation except for europium (Eu). ETD on [M + Eu - H](2+) and [M + Eu](3+) yields a limited amount of peptide backbone cleavage; however, [M + Eu + H](4+) dissociates extensively with cleavage at every residue. With the exception of the results for Eu(III), metallated peptide ion formation by ESI, ETD fragmentation efficiencies, and product ion formation are unaffected by the identity of the lanthanide cation. Adduction with trivalent lanthanide metal ions is a promising tool for sequence analysis of acidic peptides by ETD. Graphical Abstract ᅟ.

Agronomic management options for sustaining chickpea yield under climate change scenario

The impact of future climate change on the chickpea productivity was studied using the sequence analysis tool of DSSAT V 4.5 to simulate fallow-chickpea rotation at four locations viz Anantapur, Kurnool, Kadapa and Prakasam of Andhra Pradesh State. The results indicated that as compared to baseline climate, the climate change to be anticipated by 2069 (Mid –century period) would decrease the yield of chickpea by 4.3 to 18.6 per cent across various locations tested. Yield benefits obtained based on the simulation study from various adaptation options revealed that advancing the sowing window by one fortnight and application of one critical irrigation at 60 DAS found to be beneficial in increasing chickpea yields under climate change scenario.

ANTICALIgN: visualizing, editing and analyzing combined nucleotide and amino acid sequence alignments for combinatorial protein engineering.

ANTIC ALIGN: is an interactive software developed to simultaneously visualize, analyze and modify alignments of DNA and/or protein sequences that arise during combinatorial protein engineering, design and selection. ANTIC ALIGN: combines powerful functions known from currently available sequence analysis tools with unique features for protein engineering, in particular the possibility to display and manipulate nucleotide sequences and their translated amino acid sequences at the same time. ANTIC ALIGN: offers both template-based multiple sequence alignment (MSA), using the unmutated protein as reference, and conventional global alignment, to compare sequences that share an evolutionary relationship. The application of similarity-based clustering algorithms facilitates the identification of duplicates or of conserved sequence features among a set of selected clones. Imported nucleotide sequences from DNA sequence analysis are automatically translated into the corresponding amino acid sequences and displayed, offering numerous options for selecting reading frames, highlighting of sequence features and graphical layout of the MSA. The MSA complexity can be reduced by hiding the conserved nucleotide and/or amino acid residues, thus putting emphasis on the relevant mutated positions. ANTIC ALIGN: is also able to handle suppressed stop codons or even to incorporate non-natural amino acids into a coding sequence. We demonstrate crucial functions of ANTIC ALIGN: in an example of Anticalins selected from a lipocalin random library against the fibronectin extradomain B (ED-B), an established marker of tumor vasculature. Apart from engineered protein scaffolds, ANTIC ALIGN: provides a powerful tool in the area of antibody engineering and for directed enzyme evolution.

DeepSNVMiner: a sequence analysis tool to detect emergent, rare mutations in subsets of cell populations.

Background. Massively parallel sequencing technology is being used to sequence highly diverse populations of DNA such as that derived from heterogeneous cell mixtures containing both wild-type and disease-related states. At the core of such molecule tagging techniques is the tagging and identification of sequence reads derived from individual input DNA molecules, which must be first computationally disambiguated to generate read groups sharing common sequence tags, with each read group representing a single input DNA molecule. This disambiguation typically generates huge numbers of reads groups, each of which requires additional variant detection analysis steps to be run specific to each read group, thus representing a significant computational challenge. While sequencing technologies for producing these data are approaching maturity, the lack of available computational tools for analysing such heterogeneous sequence data represents an obstacle to the widespread adoption of this technology.Results. Using synthetic data we successfully detect unique variants at dilution levels of 1 in a 1,000,000 molecules, and find DeeepSNVMiner obtains significantly lower false positive and false negative rates compared to popular variant callers GATK, SAMTools, FreeBayes and LoFreq, particularly as the variant concentration levels decrease. In a dilution series with genomic DNA from two cells lines, we find DeepSNVMiner identifies a known somatic variant when present at concentrations of only 1 in 1,000 molecules in the input material, the lowest concentration amongst all variant callers tested.Conclusions. Here we present DeepSNVMiner; a tool to disambiguate tagged sequence groups and robustly identify sequence variants specific to subsets of starting DNA molecules that may indicate the presence of a disease. DeepSNVMiner is an automated workflow of custom sequence analysis utilities and open source tools able to differentiate somatic DNA variants from artefactual sequence variants that likely arose during DNA amplification. The workflow remains flexible such that it may be customised to variants of the data production protocol used, and supports reproducible analysis through detailed logging and reporting of results. DeepSNVMiner is available for academic non-commercial research purposes at https://github.com/mattmattmattmatt/DeepSNVMiner.

Advances in Molecular Serotyping and Subtyping of Escherichia coli.

Escherichia coli plays an important role as a member of the gut microbiota; however, pathogenic strains also exist, including various diarrheagenic E. coli pathotypes and extraintestinal pathogenic E. coli that cause illness outside of the GI-tract. E. coli have traditionally been serotyped using antisera against the ca. 186 O-antigens and 53 H-flagellar antigens. Phenotypic methods, including bacteriophage typing and O- and H- serotyping for differentiating and characterizing E. coli have been used for many years; however, these methods are generally time consuming and not always accurate. Advances in next generation sequencing technologies have made it possible to develop genetic-based subtyping and molecular serotyping methods for E. coli, which are more discriminatory compared to phenotypic typing methods. Furthermore, whole genome sequencing (WGS) of E. coli is replacing established subtyping methods such as pulsed-field gel electrophoresis, providing a major advancement in the ability to investigate food-borne disease outbreaks and for trace-back to sources. A variety of sequence analysis tools and bioinformatic pipelines are being developed to analyze the vast amount of data generated by WGS and to obtain specific information such as O- and H-group determination and the presence of virulence genes and other genetic markers.

An update on IgBLAST, a tool for Ig and TCR sequence analysis

Abstract The variable domain of B and T cell receptor is encoded by multiple genes, including the variable (V) gene, the diversity (D) gene and the joining (J) gene. Previously we developed an analysis tool, IgBLAST, for immunoglobulin sequences (http://www.ncbi.nlm.nih.gov/igblast/). Here we describe our recent update on this tool. The new functionalities include capability to analyze T cell receptor sequences (TR), annotating the CDR3 region, extending alignment to 5′ end of the V region, more flexibility in search criteria and providing summarized information to aid large-scale sequence analysis. We use examples to demonstrate these new functionalities and their implications in Ig and TCR sequence analysis.

Efficient Nucleic Acid Extraction and 16S rRNA Gene Sequencing for Bacterial Community Characterization.

There is a growing appreciation for the role of microbial communities as critical modulators of human health and disease. High throughput sequencing technologies have allowed for the rapid and efficient characterization of bacterial communities using 16S rRNA gene sequencing from a variety of sources. Although readily available tools for 16S rRNA sequence analysis have standardized computational workflows, sample processing for DNA extraction remains a continued source of variability across studies. Here we describe an efficient, robust, and cost effective method for extracting nucleic acid from swabs. We also delineate downstream methods for 16S rRNA gene sequencing, including generation of sequencing libraries, data quality control, and sequence analysis. The workflow can accommodate multiple samples types, including stool and swabs collected from a variety of anatomical locations and host species. Additionally, recovered DNA and RNA can be separated and used for other applications, including whole genome sequencing or RNA-seq. The method described allows for a common processing approach for multiple sample types and accommodates downstream analysis of genomic, metagenomic and transcriptional information.

Efficient Nucleic Acid Extraction and 16S rRNA Gene Sequencing for Bacterial Community Characterization

There is a growing appreciation for the role of microbial communities as critical modulators of human health and disease. High throughput sequencing technologies have allowed for the rapid and efficient characterization of bacterial communities using 16S rRNA gene sequencing from a variety of sources. Although readily available tools for 16S rRNA sequence analysis have standardized computational workflows, sample processing for DNA extraction remains a continued source of variability across studies. Here we describe an efficient, robust, and cost effective method for extracting nucleic acid from swabs. We also delineate downstream methods for 16S rRNA gene sequencing, including generation of sequencing libraries, data quality control, and sequence analysis. The workflow can accommodate multiple samples types, including stool and swabs collected from a variety of anatomical locations and host species. Additionally, recovered DNA and RNA can be separated and used for other applications, including whole genome sequencing or RNA-seq. The method described allows for a common processing approach for multiple sample types and accommodates downstream analysis of genomic, metagenomic and transcriptional information.

Management of genetic epilepsies: From empirical treatment to precision medicine

Despite the over 20 antiepileptic drugs (AEDs) now licensed for epilepsy treatment, seizures can be effectively controlled in about ⿼70% of patients. Thus, epilepsy treatment is still challenging in about one third of patients and this may lead to a severe medically, physically, and socially disabling condition. However, there is clear evidence of heterogeneity of response to existing AEDs and a significant unmet need for effective intervention. A number of studies have shown that polymorphisms may influence the poor or inadequate therapeutic response as well as the occurrence of adverse effects. In addition, the new frontier of genomic technologies, including chromosome microarrays and next-generation sequencing, improved our understanding of the genetic architecture of epilepsies. Recent findings in some genetic epilepsy syndromes provide insights into mechanisms of epileptogenesis, unrevealing the role of a number of genes with different functions, such as ion channels, proteins associated to the vesical synaptic cycle or involved in energy metabolism. The rapid progress of high-throughput genomic sequencing and corresponding analysis tools in molecular diagnosis are revolutionizing the practice and it is a fact that for some monogenic epilepsies the molecular confirmation may influence the choice of the treatment. Moreover, the novel genetic methods, that are able to analyze all known genes at a reasonable price, are of paramount importance to discover novel therapeutic avenues and individualized (or precision) medicine.

DRAMP: a comprehensive data repository of antimicrobial peptides.

The growing problem of antibiotic-resistant microorganisms results in an urgent need for substitutes to conventional antibiotics with novel modes of action and effective activities. Antimicrobial peptides (AMPs), produced by a wide variety of living organisms acting as a defense mechanism against invading pathogenic microbes, are considered to be such promising alternatives. AMPs display a broad spectrum of antimicrobial activity and a low propensity for developing resistance. Therefore, a thorough understanding of AMPs is essential to exploit them as antimicrobial drugs. Considering this, we developed a comprehensive user-friendly data repository of antimicrobial peptides (DRAMP), which holds 17349 antimicrobial sequences, including 4571 general AMPs, 12704 patented sequences and 74 peptides in drug development. Entries in the database have detailed annotations, especially detailed antimicrobial activity data (shown as target organism with MIC value) and structure information. Annotations also include accession numbers crosslinking to Pubmed, Swiss-prot and Protein Data Bank (PDB). The website of the database comes with easy-to-operate browsing as well as searching with sorting and filtering functionalities. Several useful sequence analysis tools are provided, including similarity search, sequence alignment and conserved domain search (CD-Search). DRAMP should be a useful resource for the development of novel antimicrobial peptide drugs.