writeFasta2

Write sequences in fasta format to file 
This function writes sequences from character vector as fasta formatted file (from <a href="https://www.uniprot.org">UniProt</a>) 
Line-headers are based on names of elements of input vector <code>prot</code>.
This function also allows comparing the main vector of sequences with a reference vector <code>ref</code> to check if any of the sequences therein are truncated.

Data analysis of proteomics experiments by mass spectrometry is supported by this collection of functions mostly dedicated to the analysis of (bottom-up) quantitative (XIC) data.
Fasta-formatted proteomes (eg from UniProt Consortium <doi:10.1093/nar/gky1049>) can be read with automatic parsing and multiple annotation types (like species origin, abbreviated gene names, etc) extracted.
Initial results from multiple software for protein (and peptide) quantitation can be imported (to a common format):
MaxQuant (Tyanova et al 2016 <doi:10.1038/nprot.2016.136>), Dia-NN (Demichev et al 2020 <doi:10.1038/s41592-019-0638-x>),
Fragpipe (da Veiga et al 2020 <doi:10.1038/s41592-020-0912-y>), ionbot (Degroeve et al 2021 <doi:10.1101/2021.07.02.450686>),
MassChroq (Valot et al 2011 <doi:10.1002/pmic.201100120>),
OpenMS (Strauss et al 2021 <doi:10.1038/nmeth.3959>), ProteomeDiscoverer (Orsburn 2021 <doi:10.3390/proteomes9010015>),
Proline (Bouyssie et al 2020 <doi:10.1093/bioinformatics/btaa118>), AlphaPept (preprint Strauss et al <doi:10.1101/2021.07.23.453379>)
and Wombat-P (Bouyssie et al 2023 <doi:10.1021/acs.jproteome.3c00636>.
Meta-data provided by initial analysis software and/or in sdrf format can be integrated to the analysis.
Quantitative proteomics measurements frequently contain multiple NA values, due to physical absence of given peptides in some samples, limitations in sensitivity or other reasons.
Help is provided to inspect the data graphically to investigate the nature of NA-values via their respective replicate measurements
and to help/confirm the choice of NA-replacement algorithms.
Meta-data in sdrf-format (Perez-Riverol et al 2020 <doi:10.1021/acs.jproteome.0c00376>) or similar tabular formats can be imported and included.
Missing values can be inspected and imputed based on the concept of NA-neighbours or other methods.
Dedicated filtering and statistical testing using the framework of package 'limma' <doi:10.18129/B9.bioc.limma> can be run, enhanced by multiple rounds of NA-replacements to provide robustness towards rare stochastic events.
Multi-species samples, as frequently used in benchmark-tests (eg Navarro et al 2016 <doi:10.1038/nbt.3685>, Ramus et al 2016 <doi:10.1016/j.jprot.2015.11.011>), can be run with special options considering
such sub-groups during normalization and testing. Subsequently, ROC curves (Hand and Till 2001 <doi:10.1023/A:1010920819831>) can be constructed to compare multiple analysis approaches.
As detailed example the data-set from Ramus et al 2016 <doi:10.1016/j.jprot.2015.11.011>) quantified by MaxQuant, ProteomeDiscoverer,
and Proline is provided with a detailed analysis of heterologous spike-in proteins.

Wolfgang Raffelsberger

wrProteo

Proteomics Data Analysis Functions

writeFasta2 function

<dl><dt>prot</dt>
<dd>(character) vector of sequenes, names will be used for fasta-header</dd>
<dt>fileNa</dt>
<dd>(character) name (and path) for file to be written</dd>
<dt>ref</dt>
<dd>(character) optional/additional set of (reference-) sequences (only for comparison to <code>prot</code>), length of proteins from <code>prot</code> will be checked to mark truncated proteins by '_tru'</dd>
<dt>lineLength</dt>
<dd>(integer, length=1) number of sequence characters per line (default 60, should be &gt;1 and &lt;10000)</dd>
<dt>eol</dt>
<dd>(character) the character(s) to print at the end of each line (row); for example, eol = "\r\n" will produce Windows' line endings on a Unix-alike OS</dd>
<dt>truSuf</dt>
<dd>(character) suffix to be added for sequences found truncated when comparing with <code>ref</code></dd>
<dt>silent</dt>
<dd>(logical) suppress messages</dd>
<dt>debug</dt>
<dd>(logical) supplemental messages for debugging</dd>
<dt>callFrom</dt>
<dd>(character) allows easier tracking of messages produced</dd></dl>

Arguments

Write sequences in fasta format to file 
This function writes sequences from character vector as fasta formatted file (from <a href='https://www.uniprot.org'>UniProt</a>) 
Line-headers are based on names of elements of input vector <code>prot</code>.
This function also allows comparing the main vector of sequences with a reference vector <code>ref</code> to check if any of the sequences therein are truncated.

Write sequences in fasta format to file
  
This function writes sequences from character vector as fasta formatted file (from UniProt) 
Line-headers are based on names of elements of input vector prot.
This function also allows comparing the main vector of sequences with a reference vector ref to check if any of the sequences therein are truncated. — writeFasta2

<dl>

<dt>prot</dt>
<dd>(character) vector of sequenes, names will be used for fasta-header</dd>


<dt>fileNa</dt>
<dd>(character) name (and path) for file to be written</dd>


<dt>ref</dt>
<dd>(character) optional/additional set of (reference-) sequences (only for comparison to <code>prot</code>), length of proteins from <code>prot</code> will be checked to mark truncated proteins by '_tru'</dd>


<dt>lineLength</dt>
<dd>(integer, length=1) number of sequence characters per line (default 60, should be &gt;1 and &lt;10000)</dd>


<dt>eol</dt>
<dd>(character) the character(s) to print at the end of each line (row); for example, eol = "\r\n" will produce Windows' line endings on a Unix-alike OS</dd>


<dt>truSuf</dt>
<dd>(character) suffix to be added for sequences found truncated when comparing with <code>ref</code></dd>


<dt>silent</dt>
<dd>(logical) suppress messages</dd>


<dt>debug</dt>
<dd>(logical) supplemental messages for debugging</dd>


<dt>callFrom</dt>
<dd>(character) allows easier tracking of messages produced</dd>

</dl>

Write sequences in fasta format to file
 
This function writes sequences from character vector as fasta formatted file (from <a href='https://www.uniprot.org'>UniProt</a>) 
Line-headers are based on names of elements of input vector <code>prot</code>.
This function also allows comparing the main vector of sequences with a reference vector <code>ref</code> to check if any of the sequences therein are truncated.

Write sequences in fasta format to file
 
This function writes sequences from character vector as fasta formatted file (from <a href="https://www.uniprot.org">UniProt</a>) 
Line-headers are based on names of elements of input vector <code>prot</code>.
This function also allows comparing the main vector of sequences with a reference vector <code>ref</code> to check if any of the sequences therein are truncated.

Description

Usage

Value

Arguments

Details

See Also

Examples