De Novo Sequencing Research Articles

Interaction of Salivary and Midgut Proteins of Helicoverpa armigera with Soybean Trypsin Inhibitor

Feeding of Helicoverpa armigera larvae on semi-synthetic diet containing Soybean trypsin inhibitor (STI) resulted in disappearance of STI sensitive protease in salivary and midgut protease extract. This might be due to in situ inhibition by dietary STI. STI was largely degraded within 1 h of incubation with total salivary protease (1:1). Degradation was relatively low in midgut proteases. STI interacting proteins were isolated from saliva and midgut extracts of larvae fed on STI supplemented diet using affinity column. Most of the isolated proteins showed caseinolytic activity in zymogram. Denovo sequencing data of seven different peptides selected from trypsin digested total protein showed similarity to chymotrypsinogen, serine protease, aminopeptidase N, peroxidase, hypothetical protein and muscle specific protein.

Feature-by-Feature – Evaluating De Novo Sequence Assembly

The whole-genome sequence assembly (WGSA) problem is among one of the most studied problems in computational biology. Despite the availability of a plethora of tools (i.e., assemblers), all claiming to have solved the WGSA problem, little has been done to systematically compare their accuracy and power. Traditional methods rely on standard metrics and read simulation: while on the one hand, metrics like N50 and number of contigs focus only on size without proportionately emphasizing the information about the correctness of the assembly, comparisons performed on simulated dataset, on the other hand, can be highly biased by the non-realistic assumptions in the underlying read generator. Recently the Feature Response Curve (FRC) method was proposed to assess the overall assembly quality and correctness: FRC transparently captures the trade-offs between contigs' quality against their sizes. Nevertheless, the relationship among the different features and their relative importance remains unknown. In particular, FRC cannot account for the correlation among the different features. We analyzed the correlation among different features in order to better describe their relationships and their importance in gauging assembly quality and correctness. In particular, using multivariate techniques like principal and independent component analysis we were able to estimate the “excess-dimensionality” of the feature space. Moreover, principal component analysis allowed us to show how poorly the acclaimed N50 metric describes the assembly quality. Applying independent component analysis we identified a subset of features that better describe the assemblers performances. We demonstrated that by focusing on a reduced set of highly informative features we can use the FRC curve to better describe and compare the performances of different assemblers. Moreover, as a by-product of our analysis, we discovered how often evaluation based on simulated data, obtained with state of the art simulators, lead to not-so-realistic results.

PEAKS DB: De Novo Sequencing Assisted Database Search for Sensitive and Accurate Peptide Identification

Many software tools have been developed for the automated identification of peptides from tandem mass spectra. The accuracy and sensitivity of the identification software via database search are critical for successful proteomics experiments. A new database search tool, PEAKS DB, has been developed by incorporating the de novo sequencing results into the database search. PEAKS DB achieves significantly improved accuracy and sensitivity over two other commonly used software packages. Additionally, a new result validation method, decoy fusion, has been introduced to solve the issue of overconfidence that exists in the conventional target decoy method for certain types of peptide identification software.

Secretome of the Coprophilous Fungus Doratomyces stemonitis C8, Isolated from Koala Feces

Coprophilous fungi inhabit herbivore feces, secreting enzymes to degrade the most recalcitrant parts of plant biomass that have resisted the digestive process. Consequently, the secretomes of coprophilous fungi have high potential to contain novel and efficient plant cell wall degrading enzymes of biotechnological interest. We have used one-dimensional and two-dimensional gel electrophoresis, matrix-assisted laser desorption ionization-time-of-flight tandem mass spectrometry (MALDI-TOF/TOF MS/MS), and quadrupole time-of-flight liquid chromatography-tandem mass spectrometry (Q-TOF LC-MS/MS) to identify proteins from the secretome of the coprophilous fungus Doratomyces stemonitis C8 (EU551185) isolated from koala feces. As the genome of D. stemonitis has not been sequenced, cross-species identification, de novo sequencing, and zymography formed an integral part of the analysis. A broad range of enzymes involved in the degradation of cellulose, hemicellulose, pectin, lignin, and protein were revealed, dominated by cellobiohydrolase of the glycosyl hydrolase family 7 and endo-1,4-β-xylanase of the glycosyl hydrolase family 10. A high degree of specialization for pectin degradation in the D. stemonitis C8 secretome distinguishes it from the secretomes of some other saprophytic fungi, such as the industrially exploited T. reesei. In the first proteomic analysis of the secretome of a coprophilous fungus reported to date, the identified enzymes provide valuable insight into how coprophilous fungi subsist on herbivore feces, and these findings hold potential for increasing the efficiency of plant biomass degradation in industrial processes such as biofuel production in the future.

Gapped Spectral Dictionaries and Their Applications for Database Searches of Tandem Mass Spectra

Generating all plausible de novo interpretations of a peptide tandem mass (MS/MS) spectrum (Spectral Dictionary) and quickly matching them against the database represent a recently emerged alternative approach to peptide identification. However, the sizes of the Spectral Dictionaries quickly grow with the peptide length making their generation impractical for long peptides. We introduce Gapped Spectral Dictionaries (all plausible de novo interpretations with gaps) that can be easily generated for any peptide length thus addressing the limitation of the Spectral Dictionary approach. We show that Gapped Spectral Dictionaries are small thus opening a possibility of using them to speed-up MS/MS searches. Our MS-Gapped-Dictionary algorithm (based on Gapped Spectral Dictionaries) enables proteogenomics applications (such as searches in the six-frame translation of the human genome) that are prohibitively time consuming with existing approaches. MS-Gapped-Dictionary generates gapped peptides that occupy a niche between accurate but short peptide sequence tags and long but inaccurate full length peptide reconstructions. We show that, contrary to conventional wisdom, some high-quality spectra do not have good peptide sequence tags and introduce gapped tags that have advantages over the conventional peptide sequence tags in MS/MS database searches.

Combinatorial analysis and algorithms for quasispecies reconstruction using next-generation sequencing

BackgroundNext-generation sequencing (NGS) offers a unique opportunity for high-throughput genomics and has potential to replace Sanger sequencing in many fields, including de-novo sequencing, re-sequencing, meta-genomics, and characterisation of infectious pathogens, such as viral quasispecies. Although methodologies and software for whole genome assembly and genome variation analysis have been developed and refined for NGS data, reconstructing a viral quasispecies using NGS data remains a challenge. This application would be useful for analysing intra-host evolutionary pathways in relation to immune responses and antiretroviral therapy exposures. Here we introduce a set of formulae for the combinatorial analysis of a quasispecies, given a NGS re-sequencing experiment and an algorithm for quasispecies reconstruction. We require that sequenced fragments are aligned against a reference genome, and that the reference genome is partitioned into a set of sliding windows (amplicons). The reconstruction algorithm is based on combinations of multinomial distributions and is designed to minimise the reconstruction of false variants, called in-silico recombinants.ResultsThe reconstruction algorithm was applied to error-free simulated data and reconstructed a high percentage of true variants, even at a low genetic diversity, where the chance to obtain in-silico recombinants is high. Results on empirical NGS data from patients infected with hepatitis B virus, confirmed its ability to characterise different viral variants from distinct patients.ConclusionsThe combinatorial analysis provided a description of the difficulty to reconstruct a quasispecies, given a determined amplicon partition and a measure of population diversity. The reconstruction algorithm showed good performance both considering simulated data and real data, even in presence of sequencing errors.

ADI pathway and histidine decarboxylation are reciprocally regulated in Lactobacillus hilgardii ISE 5211: proteomic evidence

Amine production by amino acid decarboxylation is a common feature that is used by lactic acid bacteria (LAB) to complement lactic fermentation, since it is coupled with a proton-extruding antiport system which leads to both metabolic energy production and the attenuation of intracellular acidity. Analogous roles are played in LAB by both malolactic fermentation (MLF) and the arginine deiminase (ADI) pathway. The present investigation was aimed at establishing reciprocal interactions between amino acid decarboxylation and the two above mentioned routes. The analyses were carried out on a Lactobacillus hilgardii strain (ISE 5211) that is able to decarboxylate histidine to histamine, which had previously been isolated from wine and whose complete genome is still unknown. The 2DE proteomic approach, followed by MALDI TOF-TOF and De Novo Sequencing, was used to study the protein expression levels. The experimental evidence has indicated that malate does not influence histidine decarboxylase (HDC) biosynthesis and that histidine does not affect the malolactic enzyme level. However, the expression of the ADI route enzymes, arginine deiminase and ornithine transcarbamylase, is down-regulated by histidine: this biosynthetic repression is more important (4-fold) in cultures that are not supplemented with arginine, but is also significant (2-fold) in an arginine supplemented medium that normally induces the ADI pathway. On the other hand, arginine partially represses HDC expression, but only when histidine and arginine are both present in the culture medium. This proteomic study has also pointed out a down-regulation exerted by histidine over sugar metabolism enzymes and a GroEL stress protein. These data, together with the reciprocal antagonism between arginine deimination and histidine decarboxylation, offer clue keys to the understanding of the accumulation of lactate, amine, ammonia and ethylcarbamate in wine, with consequent implications on different health risk controls.

Characterization of Quasispecies of Pandemic 2009 Influenza A Virus (A/H1N1/2009) by De Novo Sequencing Using a Next-Generation DNA Sequencer

Pandemic 2009 influenza A virus (A/H1N1/2009) has emerged globally. In this study, we performed a comprehensive detection of potential pathogens by de novo sequencing using a next-generation DNA sequencer on total RNAs extracted from an autopsy lung of a patient who died of viral pneumonia with A/H1N1/2009. Among a total of 9.4×106 40-mer short reads, more than 98% appeared to be human, while 0.85% were identified as A/H1N1/2009 (A/Nagano/RC1-L/2009(H1N1)). Suspected bacterial reads such as Streptococcus pneumoniae and other oral bacteria flora were very low at 0.005%, and a significant bacterial infection was not histologically observed. De novo assembly and read mapping analysis of A/Nagano/RC1-L/2009(H1N1) showed more than ×200 coverage on average, and revealed nucleotide heterogeneity on hemagglutinin as quasispecies, specifically at two amino acids (Gly172Glu and Gly239Asn of HA) located on the Sa and Ca2 antigenic sites, respectively. Gly239 and Asn239 on antigenic site Ca2 appeared to be minor amino acids compared with the highly distributed Asp239 in H1N1 HAs. This study demonstrated that de novo sequencing can comprehensively detect pathogens, and such in-depth investigation facilitates the identification of influenza A viral heterogeneity. To better characterize the A/H1N1/2009 virus, unbiased comprehensive techniques will be indispensable for the primary investigations of emerging infectious diseases.

De Novo Sequencing of Peptides on Single Resin Beads by MALDI-FTICR Tandem Mass Spectrometry

An efficient approach in combinatorial chemistry is the synthesis of one-bead-one-compound peptide libraries. In contrast to synthesis and functional screening, which is performed in a largely automated manner, structure determination has been frequently laborious and time-consuming. Here we report an approach for de novo sequencing of peptides on single beads by matrix-assisted laser desorption/ionization Fourier transform ion cyclotron resonance (MALDI-FTICR) tandem mass spectrometry, using a resin with a photolinker for solid-phase peptide synthesis. Upon sorting out single beads, an efficient sample preparation on the MALDI target was developed that enables fragmentation upon irradiation of the bead-matrix mixture with the ultraviolet (UV)-MALDI laser, with enhanced yield of sequence-specific fragment ions at increased laser energy. This approach is illustrated by sequence determinations of two peptides from a library with sequences varying in a single amino acid; the feasibility with tandem-MS procedures and fragment ion assignment was ascertained by sustained off-resonance irradiation/collision induced dissociation (SORI/CID) and infrared multiphoton dissociation (IRMPD) fragmentation.

N-Terminal Tagging Strategy for De Novo Sequencing of Short Peptides by ESI-MS/MS and MALDI-MS/MS

The major portion of skin secretory peptidome of the European Tree frog Hyla arborea consists of short peptides from tryptophyllin family. It is known that b-ions of these peptides undergo head-to-tail cyclization, forming a ring that can open, resulting in several linear forms. As a result, the spectrum contains multiple ion series, thus complicating de novo sequencing. This was observed in the Q-TOF spectrum of one of the tryptophyllins isolated from Hyla arborea; the sequence FLPFFP-NH(2) was established by Edman degradation and counter-synthesis. Though no rearrangements were observed in FTICR-MS and MALDI-TOF/TOF spectra, both of them were not suitable for mass-spectrometry sequencing due to the low sequence coverage. To obtain full amino acid sequence by mass spectrometry, three chemical modifications to N-terminal amino moiety were applied. They include acetylation and sulfobenzoylation of N-amino group and its transformation to 2,4,6-trimethylpyridinium by interaction with 2,4,6-trimethylpyrillium tetrafluoroborate. All three reagents block scrambling and provide spectra better than the intact peptide. Unfortunately, all of them also readily react with lysine side chain. Hence, all investigated procedures can be used to improve sequencing of short peptides, while acetylation is the recommended one. It shows excellent results, and it is plain and simple to perform. This is the procedure of choice for MS-sequencing of short peptides by manual or automatic algorithms.

Incorporating Biological Mass Spectrometry into Undergraduate Teaching Labs, Part 3: De Novo Peptide Sequencing Using Electrospray Tandem Mass Spectrometry

Tandem mass spectrometry (MS–MS) is an important technique routinely used to provide structural information about analytes of interest. One major application of tandem mass spectrometry is in the field of proteomics, which routinely uses MS–MS experiments to help identify peptides and proteins. In this experiment, students learn how to perform a tandem mass spectrometry experiment on a selected peptide ion using an ion-trap mass spectrometer and how to use the resulting data to sequence the peptide de novo. While the majority of protein and peptide identification done using mass spectrometry relies on database search algorithms, it is still important for a trained scientist to compare the computer results with the empirical data to make a judgment as to how well the two correlate. It is therefore important that students be exposed to peptide sequencing and learn how this can be done without a database to prepare them for research in the proteomics field. We have included several representative peptide fra...

Flipping NextGen: using biological systems to characterize NextGen sequencing technologies

Finally, de-novo sequencing (i.e. without a reference) would require a minimum of 1/2 of the sequence length to be unique in order to allow sufficient contig extension in the assembly process. For example, 40–50 bp reads are necessary for de-novo characterization of these systems uniquely defined by 20–25 bp reads. As of 2009, short read NextGen sequencing technologies have moved to 50 bp and beyond, ushering in what is expected to be the start of a revolution in genomics. Conclusion These results establish a lower bound on sequence length (x) required to sufficiently conduct re-sequencing, transcriptome, and de-novo sequencing projects. The asymptotic nature of the results also provides a guide for what percentage of the total space (y) we might expect to define in genomes/transcriptomes of similar size and complexity. from UT-ORNL-KBRIN Bioinformatics Summit 2009 Pikeville, TN, USA. 20–22 March 2009

De Novo Sequencing of RCB-1 to -3: Peptide Biomarkers from the Castor Bean PlantRicinus communis

Ricinus communis (also know as the castor bean plant) whose forbears escaped from suburban gardens or commercial cultivation grow wild in many countries. In temperate and tropical climates seeds will develop to maturity, and plants may be perennial. In Australia these plants have become widespread and are regarded as noxious weeds in many localities. The seeds of R. communis contain ricin, a protein toxin which can easily be extracted into an aqueous solution. Ricin is toxic by ingestion, inhalation, and injection. The history of terrorist and anarchist interest in the use of seeds from R. communis has driven the development of strategies for determination of cultivar and geographic location of the source of an extract of wild-grown castor bean seed. This forensic information is of considerable interest to law enforcement and intelligence organizations. During forensic studies of both the metabolome and proteome of extracts from eight specimens of six different cultivars of R. communis ("zanzibariensis" collected from Kenya and Tanzania, "gibsonii", "impala", "dehradun", "carmencita", and "sanguineus" collected from Spain and Tanzania), three peptide biomarkers (designated Ricinus communis biomarkers, or RCB) were identified in both the MALDI and electrospray LC-MS spectra. Two of these peptides (RCB-1 and RCB-2) were present in varying amounts in all cultivars, while RCB-3 was present only in the "carmencita" cultivar. The amino acid sequences of RCB-1 to -3 were determined using LC-MS(n) fragmentation and de novo sequencing on both the intact and the carbamidomethyl modified peptides. The connectivity of the two disulfide bonds that were present in all three RCB were determined using a strategy of partial reduction and differential alkylation using tris-(2-carboxyethyl)phosphine with N-ethylmaleimide to reduce and alkylate the most accessible disulfide bond, followed by reduction and alkylation of the remaining disulfide bond with dithiolthreitol and iodoacetamide. The possible functional role of RCB-1 to -3 in R. communis seeds is also discussed.

Sequencing of Single and Double Stranded RNA Oligonucleotides by Acid Hydrolysis and MALDI Mass Spectrometry

Treatment of RNA oligonucleotides with strong acids at pH 1-2 rapidly leads to hydrolysis of the phosphodiester bonds at the 5'-position of ribose. Analysis of the resulting degradation products by MALDI coupled to an Orbitrap high resolution mass spectrometer shows almost complete mass ladders from both sides of the nucleotides without interfering fragments from base losses or internal fragments. From the mass differences between adjacent peaks of a mass ladder, the sequence can be determined. Low cleavage efficiency at the termini leads to 2mers and 3mers which can be identified by MS/MS. In this way the complete sequences of different siRNA 21mer single and double strands could be verified. This simple and fast method can be applied for controlling sequences of synthetic oligomers, as well as for de-novo sequencing. Moreover, the method is applicable for localization and identification of RNA modifications as demonstrated using the examples of an oligonucleotide with phosphorothioate backbone and of one containing 2'-methoxy-ribose modifications.

Spectral Profiles, a Novel Representation of Tandem Mass Spectra and Their Applications for de Novo Peptide Sequencing and Identification

Despite many efforts in the last decade, the progress in de novo peptide sequencing has been slow with only 30-45% of all peptides correctly reconstructed. We argue that accurate full-length peptide sequencing may be an unattainable goal for some spectra and demonstrate how to accurately sequence gapped peptides instead. We further argue that gapped peptides are nearly as useful as full-length peptides for error-tolerant database searches. Gapped peptides occupy a niche between long but inaccurate full-length reconstructions and short but accurate peptide sequence tags. Our MS-Profile tool uses spectral profiles, a new representation of tandem mass spectra, to generate gapped peptides that are longer and more accurate than peptide sequence tags of length 3 traditionally used to speed up database searches in proteomics. In addition, spectral profiles also enable intuitive visualization of all high scoring de novo reconstructions of tandem mass spectra.

De Novo Sequencing of Unique Sequence Tags for Discovery of Post-Translational Modifications of Proteins

De novo sequencing is a spectrum analysis approach for mass spectrometry data to discover post-translational modifications in proteins; however, such an approach is still in its infancy and is still not widely applied to proteomic practices due to its limited reliability. In this work, we describe a de novo sequencing approach for the discovery of protein modifications based on identification of the proteome UStags (Shen, Y.; Tolić, N.; Hixson, K. K.; Purvine, S. O.; Pasa-Tolić, L.; Qian, W. J.; Adkins, J. N.; Moore, R. J.; Smith, R. D. Anal. Chem. 2008, 80, 1871-1882). The de novo information was obtained from Fourier-transform tandem mass spectrometry data for peptides and polypeptides from a yeast lysate, and the de novo sequences obtained were selected based on filter levels designed to provide a limited yet high quality subset of UStags. The DNA-predicted database protein sequences were then compared to the UStags, and the differences observed across or in the UStags (i.e., the UStags' prefix and suffix sequences and the UStags themselves) were used to infer possible sequence modifications. With this de novo-UStag approach, we uncovered some unexpected variances within several yeast protein sequences due to amino acid mutations and/or multiple modifications to the predicted protein sequences. To determine false discovery rates, two random (false) databases were independently used for sequence matching, and ~3% false discovery rates were estimated for the de novo-UStag approach. The factors affecting the reliability (e.g., existence of de novo sequencing noise residues and redundant sequences) and the sensitivity of the approach were investigated and described. The combined de novo-UStag approach complements the UStag method previously reported by enabling the discovery of new protein modifications.

Accurate Mass as a Bioinformatic Parameter in Data-to-Knowledge Conversion: Fourier Transform Ion Cyclotron Resonance Mass Spectrometry for Peptide De Novo Sequencing

With the availability of ultra-precise mass spectrometric biomolecular data, the accurate mass is becoming a physical quantity of high interest for bioinformatics tools and strategies. Fourier transform ion cyclotron resonance mass spectrometry with electrospray ionization or matrix-assisted laser desorption/ionization sources now allows the easy determination of amino acid composition of medium size, unknown peptides when employing combinatorial calculation of hypothetical parent and fragment ion masses. This new method, which in a second step, allows the reliable de-novo sequencing of completely unknown peptides ["composition-based sequencing (CBS)"(1)] appears to open a wide new field of bioanalytical investigation. It has been shown that even unspecifically cleaved proteins can be identified easily and reliably when accurate mass evaluation is combined with protein database search tools.(2) It is quite clear that, while the nominal mass of a peptide has obviously no useful correlation to biomolecular structure, the accurate mass, instead, has a strong and direct correlation to structure that so far has not been exploited considerably by bioinformatic tools. It has already become obvious that accurate mass evaluation is going to become a central goal for bioinformatics strategies in the near future.(3-11) Strategies for extracting structural, and even functional, information out of accurate mass values of biomolecules still have to be developed. Examples and prospects of accurate mass evaluation in bioinformatics for the field of proteomics are outlined in the following.

Lookup Peaks: A Hybrid of de Novo Sequencing and Database Search for Protein Identification by Tandem Mass Spectrometry

A powerful technique for peptide and protein identification is tandem mass spectrometry followed by database search using a program such as SEQUEST or Mascot. These programs, however, become slow and lose sensitivity when allowing nonspecific cleavages or peptide modifications. De novo sequencing and hybrid methods such as sequence tagging offer speed and robustness for wider searches, yet these approaches require better spectra with more complete and consecutive fragmentation and, hence, are less sensitive to low-abundance peptides. Here we describe a new hybrid method that retains the sensitivity of pure database search. The method uses a small amount of de novo analysis to identify likely b- and y-ion peaks--"lookup peaks"--that can then be used to extract candidate peptides from the database, with the number of candidates tunable to fit a computing budget. We describe a program called ByOnic that implements this method, and we benchmark ByOnic on several data sets, including one of mouse blood plasma spiked with low concentrations of recombinant human proteins. We demonstrate that ByOnic is more sensitive than sequence tagging and, indeed, more sensitive than the three most popular pure database search tools--SEQUEST, Mascot, and X!Tandem--on both the peptide and protein levels. On the mouse plasma samples, ByOnic consistently found spiked proteins missed by the other tools.

Improving de Novo Sequencing of Peptides Using a Charged Tag and C-Terminal Digestion

An improved method for peptide de novo sequencing by MALDI mass spectrometry is presented. The method couples a charge derivatization reaction with C-terminal digestion to modify tryptic peptides. The charge derivatization attaches a fixed charge group onto the N-termini of peptides, and the enzymatic digestion after the derivatization step removes C-terminal basic amino acid residues such as arginine and lysine. The fragmentation of the modified peptide(s) under low-energy CID conditions (MALDI Q-TOF mass spectrometer) yields a simplified yet complete ion series of the peptide sequence. The validity of the method is demonstrated by the results from several model protein digests, where peptide sequences were correctly deduced either manually or through an automated sequencing program.

Production and physicochemical characterization of acidocin D20079, a bacteriocin produced by Lactobacillus acidophilus DSM 20079

Lactobacillus acidophilus DSM 20079 is the producer of a novel bacteriocin termed acidocin D20079. In this paper, a partial sequence of this peptide is determined, together with data on its secondary structure. A modification of the MRS-growth medium (replacing the detergent Tween 80 with oleic acid), was shown to improve the production level of the peptide by one order of magnitude, as well as to stabilize the activity level. Addition of a detergent (Tween 20, less interfering in mass spectrometric analysis), was however necessary for solubilization of the purified acidocin D20079. Digestion of the peptide followed by de-novo sequencing of generated fragments, allowed determination of a partial sequence consisting of 39 of the totally estimated 65 residues. Acidocin D20079 has a high content of glycine residues, hydrophobic residues, and acidic residues. No modified amino acids were found. Edman degradation, and C-terminal sequencing failed, suggesting that the peptide may be cyclic, and a novel member of class IIc bacteriocins. Circular dichroism spectroscopy and secondary structure prediction showed random coil conformation in aqueous solution, but secondary structure was induced in the presence of sodium-dodecyl sulfate. The data could be fitted assuming 2–13% of the residues to be in α-helix and 23–27% of the residues to be in β-strand conformation. This indicates that a membrane/membrane-mimicking hydrocarbon–water interface induces an active conformation.