The files in the subdirectories of
extdata support the examples in the package documentation and vignettes.
abundance contain protein abundance and microbial occurrence data:
stress is a data frame listing proteins identified in selected proteomic stress response experiments. The names of proteins begin at row 3, and columns are all the same length (padded as necessary at the bottom by
NAs). Names correspond to ordered locus names (for Sce) or gene names (for Eco). The column names identify the experiments, the first row contains the name of the organism (Sce or Eco) and the third row has the reference key for the source of the data (listed in
AA03.csv has reference abundances for 71 proteins taken from Fig. 3 of Anderson and Anderson, 2002 (as corrected in Anderson and Anderson, 2003). The columns with data taken from these sources are type (hemoglobin, plasma, tissue, or interleukin), description (name used in the original figure), log10(pg/ml) (upper limit of abundance interval shown in Anderson and Anderson, 2003, log10 of concentration in pg/ml). The additional columns are data derived from a search of the SWISS-PROT/UniProtKB database based on the descriptions of the proteins: name (nominal UniProtKB name for this protein), name2 (other UniProtKB names(s) that could apply to the protein), and note (notes based on searching for a protein of this description). The amino acid compositions of all proteins whose names are not NA are included in
abbrv column for the proteins contains the description given by Anderson and Anderson, 2003, followed by (in parentheses) the UniProtKB accession number. Annotated initiator methionines (e.g. for ferritin, myoglobin, ENOG), signal peptides or propeptides were removed from the proteins (except where they are not annotated in UniProtKB: IGHG1, IGHA1, IGHD, MBP). In cases were multiple isoforms are present in UniProtKB (e.g. Albumin) only the first isoform was taken. In the case of C4 Complement (CO4A) and C5 Complement (CO5), the amino acid composition of only the alpha chains are listed. In the case of the protein described as iC3b, the amino acid sequence is taken to be that of Complement C3c alpha' chain fragment 1 from CO3, and is given the name CO3.C3c. The non-membrane (soluble) chains of TNF-binding protein (TNR1A) and TNF-alpha (TNFA) were used. Rantes, MIP-1 beta and MIP-1 alpha were taken from C-C motif chemokines (CCL5, CCL4, CCL3 respectively). C-peptide was taken from the corresponding annotation for insulin and here is named INS.C. See
read.expr for examples that use this file.
ISR+08.csv has columns excerpted from Additional File 2 of Ishihama et al. (2008) for protein abundances in E. coli cytosol. The columns in this file are ID (Swiss-Prot ID), accession (Swiss-Prot accession), emPAI (exponentially modified protein abundance index), copynumber (emPAI-derived copy number/cell), GRAVY (Kyte-Doolittel), FunCat (FunCat class description), PSORT (PSORT localisation), ribosomal (yes/no). See
read.expr for examples that use this file.
yeastgfp.csv.xz Has 28 columns; the names of the first five are
GFP visualized?, and
abundance. The remaining columns correspond to the 23 subcellular localizations considered in the YeastGFP project (Huh et al., 2003 and Ghaemmaghami et al., 2003) and hold values of either
F for each protein. yeastgfp.csv was downloaded on 2007-02-01 from http://yeastgfp.ucsf.edu using the Advanced Search, setting options to download the entire dataset and to include localization table and abundance, sorted by orf number. See
yeastgfp for examples that use this file.
microbes.csv has data for microbial occurrence (i.e. relative enrichement) in colorectal cancer and normal tissue. The file is from the Supporting Information of Dick (2016). This file is used by
bison contain BLAST results and taxonomic information for a metagenome:
bisonP... are partial tabular BLAST results for proteins in the Bison Pool Environmental Genome. Protein sequences predicted in the metagenome were downloaded from the Joint Genome Institute's IMG/M system on 2009-05-13. The target database for the searches was constructed from microbial protein sequences in National Center for Biotechnology Information (NCBI) RefSeq database version 57, representing 7415 microbial genomes. The ‘blastall’ command was used with the default setting for E value cuttoff (10.0) and options to make a tabular output file consisting of the top 20 hits for each query sequence. The function
read.blast was used to extract only those hits with E values less than or equal to 1e-5 and with sequence similarity (percent identity) at least 30 percent, and to keep only the first hit for each query sequence. The function
write.blast was used to save partial BLAST files (only selected columns). The files provided with CHNOSZ contain the first 5,000 hits for each sampling site at Bison Pool, representing between about 7 to 15 percent of the first BLAST hits after similarity and E value filtering.
gi.taxid.txt.xz is a table that lists the sequence identifiers (gi numbers) that appear in the example BLAST files (see above), together with the corresponding taxon ids used in the NCBI databases. This file is not a subset of the complete ‘gi_taxid_prot.dmp.gz’ available at ftp://ftp.ncbi.nih.gov/pub/taxonomy/ but instead is a subset of ‘gi.taxid.txt’ generated from the RefSeq release catalog using ‘gencat.sh’ in the
refseq directory. See
id.blast for an example that uses this file and the BLAST files described above.
cpetc contain heat capacity and other thermodynamic data and parameters:
PM90.csv Heat capacities of four unfolded aqueous proteins taken from Privalov and Makhatadze, 1990. Temperature in °C is in the first column, and heat capacities of the proteins in J mol^-1 K^-1 in the remaining columns. See
ionize.aa for an example that uses this file.
RH95.csv Heat capacity data for iron taken from Robie and Hemingway, 1995. Temperature in Kelvin is in the first column, heat capacity in J K^-1 mol^-1 in the second. See
subcrt for an example that uses this file.
RT71.csv pH titration measurements for unfolded lysozyme (LYSC_CHICK) taken from Roxby and Tanford, 1971. pH is in the first column, net charge in the second. See
ionize.aa for an example that uses this file.
SOJSH.csv Experimental equilibrium constants for the reaction NaCl(aq) = Na+ + Cl- as a function of temperature and pressure taken from Fig. 1 of Shock et al., 1992. Data were extracted from the figure using g3data (http://www.frantz.fi/software/g3data.php). See
water for an example that uses this file.
V.CH4.HWM96.csv Apparent molar heat capacities and volumes of CH4 in dilute aqueous solutions reported by Hnedkovsky and Wood, 1997 and Hnedkovsky et al., 1996. See
EOSregress for examples that use these files.
SC10_Rainbow.csv Values of temperature (°C), pH and logarithms of activity of CO2, H2, NH4+, H2S and CH4 for mixing of seawater and hydrothermal fluid at Rainbow field (Mid-Atlantic Ridge), taken from Shock and Canovas, 2010.
SS98_Fig5b.csv Values of logarithm of fugacity of O2 and pH as a function of temperature for mixing of seawater and hydrothermal fluid, digitized from Figs. 5a and b of Shock and Schulte, 1998.
rubisco.csv UniProt IDs for Rubisco, ranges of optimal growth temperature of organisms, domain and name of organisms, and URL of reference for growth temperature, from Dick, 2014.
fasta contain protein sequences:
EF-Tu.aln consists of aligned sequences (394 amino acids) of elongation factor Tu (EF-Tu). The sequences correspond to those taken from UniProtKB for ECOLI (Escherichia coli), THETH (Thermus thermophilus) and THEMA (Thermotoga maritima), and reconstructed ancestral sequences taken from Gaucher et al., 2003 (maximum likelihood bacterial stem and mesophilic bacterial stem, and alternative bacterial stem)
rubisco.csv Sequences of Rubisco obtained from UniProt (see Dick, 2014).
protein contain protein composition data for model organisms.
read.expr for examples that use these files.
Data frame of amino acid composition of 6716 proteins from the Saccharomyces Genome Database (SGD).
Values in the first three columns are the
ORF names of proteins,
GENE names. The remaining twenty columns (
VAL) contain the numbers of the respective amino acids in each protein.
The sources of data for Sce.csv are the files protein_properties.tab and SGD_features.tab (for the gene names), downloaded from http://www.yeastgenome.org on 2013-08-24.
Amino acid compositions of 4407 proteins in Escherichia coli strain K12.
Format is the one used
thermo$protein, with columns protein holding the gene name, organism set to ECOLI, and
abbrv holding the UniProt ID.
The source of data is the file ECOLI.fas downloaded from the HAMAP (High-quality Automated and Manual Annotation of microbial Proteomes system) FTP site (Gattiker et al., 2003) on 2010-09-25 (old URL:
These two files contain amino acid compositions of metagenomically encoded proteins, averaged together according to functional annotation (DS11) or taxonomic affiliation (DS13).
The data are from Dick and Shock, 2011 and 2013.
They are used in the vignette Hot-spring proteins in CHNOSZ.
Overall protein compositions of microbial species reported to be positively or negatively enriched in colorectal cancer.
This file is taken from Dick, 2016.
It is used by
refseq contain code and results of processing NCBI Reference Sequences (RefSeq) for microbial proteins, using RefSeq release 61 of 2013-09-09:
README.txt Instructions for producing the data files.
gencat.sh Bash script to extract microbial protein records from the RefSeq catalog.
gi.taxid.txt Output from above. The complete file is too large to distribute with CHNOSZ, but a portion is included in
extdata/bison to support processing example BLAST files for the Bison Pool metagenome (based on RefSeq 57, 2013-01-08).
mkfaa.sh Combine the contents of .faa.gz files into a single FASTA file (to use e.g. for making a BLAST database).
protein.refseq.R Calculate average amino acid composition of all proteins for each organism identified by a taxonomic ID.
trim_refseq.R Keep only selected organism names (reduces number of taxa from 6758 to 779, helps to control package size).
protein_refseq.csv.xz Output from above. See example in
taxid.names.R Generate a table of scientific names for the provided taxids. Requires the complete
nodes.dmp from NCBI taxonomy files.
taxid_names.csv.xz Output from above.
NOTE: For backward compatibility with the example BLAST files for the Bison Pool metagenome, the packaged file merges records for taxids found in either RefSeq 57 or 61.
Certain taxids in release 57 were not located in the current RefSeq catalog, probably related to the transition to the “WP” multispecies accessions (ftp://ftp.ncbi.nlm.nih.gov/refseq/release/announcements/WP-proteins-06.10.2013.pdf).
See example for
taxonomy contain example taxonomic data files:
nodes.dmp are excerpts of the taxonomy files available on the NCBI ftp site (ftp://ftp.ncbi.nih.gov/pub/taxonomy/taxdump.tar.gz, accessed 2010-02-15). These example files contain only the entries for Escherichia coli K-12, Saccharomyces cerevisiae, Homo sapiens, Pyrococcus furisosus and Methanocaldococcus jannaschii (taxids 83333, 4932, 9606, 186497, 243232) and the higher-ranking nodes (genus, family, etc.) in the respective lineages. See
taxonomy for examples that use this file.
thermo contain additional thermodynamic data and group additivity definitions:
OBIGT-2.csv contains supplementary thermodynamic data in the same format as the primary database in
data/OBIGT.csv. Data for some entries in the primary database are taken from different literature sources in this file. The default action of
add.obigt is to add the contents of this file to CHNOSZ's working database in
diagram and the code of
anim.TCA for examples that use this file.
obigt_check.csv contains the results of running
check.obigt to check the internal consistency of entries in the primary and supplementary databases.
RH98_Table15.csv Group stoichiometries for high molecular weight crystalline and liquid organic compounds taken from Table 15 of Richard and Helgeson, 1998. The first three columns have the
formula and physical
state (cr or liq). The remaining columns have the numbers of each group in the compound; the names of the groups (columns) correspond to species in
thermo$obigt. The compound named 5a(H),14a(H)-cholestane in the paper has been changed to 5a(H),14b(H)-cholestane here to match the group stoichiometry given in the table. See
RH2obigt for a function that uses this file.
DLEN67.csv Standard Gibbs energies of formation, in kcal/mol, from Dayhoff et al., 1967, for nitrogen (N2) plus 17 compounds shown in Fig. 2 of Dayhoff et al., 1964, at 300, 500, 700 and 1000 K.
Anderson, N. L. and Anderson, N. G. (2002) The human plasma proteome: History, character and diagnostic prospects. Molecular and Cellular Proteomics 1, 845–867. http://dx.doi.org/10.1074/mcp.R200007-MCP200
Anderson, N. L. and Anderson, N. G. (2003) The human plasma proteome: History, character and diagnostic prospects (Vol. 1 (2002) 845-867). Molecular and Cellular Proteomics 2, 50. http://dx.doi.org/10.1074/mcp.A300001-MCP200
Baas Becking, L. G. M., Kaplan, I. R. and Moore, D. (1960) Limits of the natural environment in terms of pH and oxidation-reduction potentials. Journal of Geology 68(3), 243–284. http://www.jstor.org/stable/30059218
Benson, S. W. and Buss, J. H. (1958) Additivity rules for the estimation of molecular properties. Thermodynamic properties. J. Chem. Phys. 29, 546–572. http://dx.doi.org/10.1063/1.1744539
Dayhoff, M. O. and Lippincott, E. R. and Eck, R. V. (1964) Thermodynamic Equilibria In Prebiological Atmospheres. Science 146, 1461–1464. http://dx.doi.org/10.1126/science.146.3650.1461
Dayhoff, M. O. and Lippincott, E. R., Eck, R. V. and Nagarajan (1967) Thermodynamic Equilibrium In Prebiological Atmospheres of C, H, O, N, P, S, and Cl. Report SP-3040, National Aeronautics and Space Administration. http://ntrs.nasa.gov/search.jsp?R=19670017966
Dick, J. M. (2014) Average oxidation state of carbon in proteins. J. R. Soc. Interface 11, 20131095. http://dx.doi.org/10.1098/rsif.2013.1095
Dick, J. M. (2016) Proteomic indicators of oxidation and hydration state in colorectal cancer. PeerJ 4:e2238. http://dx.doi.org/10.7717/peerj.2238
Dick, J. M. and Shock, E. L. (2011) Calculation of the relative chemical stabilities of proteins as a function of temperature and redox chemistry in a hot spring. PLoS ONE 6, e22782. http://dx.doi.org/10.1371/journal.pone.0022782
Dick, J. M. and Shock, E. L. (2013) A metastable equilibrium model for the relative abundance of microbial phyla in a hot spring. PLoS ONE 8, e72395. http://dx.doi.org/10.1371/journal.pone.0072395
Domalski, E. S. and Hearing, E. D. (1993) Estimation of the thermodynamic properties of C-H-N-O-S-Halogen compounds at 298.15 K J. Phys. Chem. Ref. Data 22, 805–1159. http://dx.doi.org/10.1063/1.555927
Gattiker, A., Michoud, K., Rivoire, C., Auchincloss, A. H., Coudert, E., Lima, T., Kersey, P., Pagni, M., Sigrist, C. J. A., Lachaize, C., Veuthey, A.-L., Gasteiger, E. and Bairoch, A. (2003) Automatic annotation of microbial proteomes in Swiss-Prot. Comput. Biol. Chem. 27, 49–58. http://dx.doi.org/10.1016/S1476-9271(02)00094-4
Gaucher, E. A., Thomson, J. M., Burgan, M. F. and Benner, S. A (2003) Inferring the palaeoenvironment of ancient bacteria on the basis of resurrected proteins. Nature 425(6955), 285–288.
Ghaemmaghami, S., Huh, W., Bower, K., Howson, R. W., Belle, A., Dephoure, N., O'Shea, E. K. and Weissman, J. S. (2003) Global analysis of protein expression in yeast. Nature 425(6959), 737–741.
Huh, W. K., Falvo, J. V., Gerke, L. C., Carroll, A. S., Howson, R. W., Weissman, J. S. and O'Shea, E. K. (2003) Global analysis of protein localization in budding yeast. Nature 425(6959), 686–691.
HAMAP system. HAMAP FTP directory, ftp://ftp.expasy.org/databases/hamap/
Hnedkovsky, L., Wood, R. H. and Majer, V. (1996) Volumes of aqueous solutions of CH4, CO2, H2S, and NH3 at temperatures from 298.15 K to 705 K and pressures to 35 MPa. J. Chem. Thermodyn. 28, 125–142. http://dx.doi.org/10.1006/jcht.1996.0011
Hnedkovsky, L. and Wood, R. H. (1997) Apparent molar heat capacities of aqueous solutions of CH4, CO2, H2S, and NH3 at temperatures from 304 K to 704 K at a pressure of 28 MPa. J. Chem. Thermodyn. 29, 731–747. http://dx.doi.org/10.1006/jcht.1997.0192
Ishihama, Y., Schmidt, T., Rappsilber, J., Mann, M., Hartl, F. U., Kerner, M. J. and Frishman, D. (2008) Protein abundance profiling of the Escherichia coli cytosol. BMC Genomics 9:102. http://dx.doi.org/10.1186/1471-2164-9-102
Joint Genome Institute (2007) Bison Pool Environmental Genome. Protein sequence files downloaded from IMG/M (http://img.jgi.doe.gov/cgi-bin/m/main.cgi?section=FindGenomes&page=findGenomes)
Privalov, P. L. and Makhatadze, G. I. (1990) Heat capacity of proteins. II. Partial molar heat capacity of the unfolded polypeptide chain of proteins: Protein unfolding effects. J. Mol. Biol. 213, 385–391. http://dx.doi.org/10.1016/S0022-2836(05)80198-6
Richard, L. and Helgeson, H. C. (1998) Calculation of the thermodynamic properties at elevated temperatures and pressures of saturated and aromatic high molecular weight solid and liquid hydrocarbons in kerogen, bitumen, petroleum, and other organic matter of biogeochemical interest. Geochim. Cosmochim. Acta 62, 3591–3636. http://dx.doi.org/10.1016/S0016-7037(97)00345-1
Robie, R. A. and Hemingway, B. S. (1995) Thermodynamic Properties of Minerals and Related Substances at 298.15 K and 1 Bar (10^5 Pascals) Pressure and at Higher Temperatures. U. S. Geol. Surv., Bull. 2131, 461 p. http://www.worldcat.org/oclc/32590140
Roxby, R. and Tanford, C. (1971) Hydrogen ion titration curve of lysozyme in 6 M guanidine hydrochloride. Biochemistry 10, 3348–3352. http://dx.doi.org/10.1021/bi00794a005
SGD project. Saccharomyces Genome Database, http://www.yeastgenome.org
Shock, E. and Canovas, P. (2010) The potential for abiotic organic synthesis and biosynthesis at seafloor hydrothermal systems. Geofluids 10, 161–192. http://dx.doi.org/10.1111/j.1468-8123.2010.00277.x
Shock, E. L., Oelkers, E. H., Johnson, J. W., Sverjensky, D. A. and Helgeson, H. C. (1992) Calculation of the thermodynamic properties of aqueous species at high pressures and temperatures: Effective electrostatic radii, dissociation constants and standard partial molal properties to 1000 °C and 5 kbar. J. Chem. Soc. Faraday Trans. 88, 803–826. http://dx.doi.org/10.1039/FT9928800803
Shock, E. L. and Schulte, M. D. (1998) Organic synthesis during fluid mixing in hydrothermal systems. J. Geophys. Res. 103, 28513–28527. http://dx.doi.org/10.1029/98JE02142
YeastGFP project. Yeast GFP Fusion Localization Database, http://yeastgfp.ucsf.edu; Current location: http://yeastgfp.yeastgenome.org