Difference between revisions of "MSigDB v3.1 Release Notes"
m |
m |
||
| Line 2: | Line 2: | ||
<h2>Gene Symbols Update</h2> | <h2>Gene Symbols Update</h2> | ||
<p> Gene sets consist of a large variety of gene identifiers, called <strong>original members</strong> here. To use gene sets by GSEA and other querying tools, original members have to be converted to a common universal kind of gene identifiers. Previous releases of MSigDB used human gene symbols for this purpose. Researchers prefer working with gene symbols because they can easily recognize, remember and put them in the context of their work. Unfortunately, a gene usually has multiple different symbols. Conversely, the same symbol can often refer to a number of different genes. Finally, gene symbols change frequently. To overcome these issues, we chose Entrez Gene IDs as robust universal identifiers (called <strong>ezid members</strong> here. Entrez Gene IDs uniquely identify human genes and never change. For convenience, we continue displaying gene sets as made from human gene symbols. However, the symbols are now unambiguously derived from the corresponding human Entrez Gene IDs. For non-human original members, we first convert them to the organism-specific Entrez Gene IDs and then seek their orthologous counterparts as human Entrez Gene IDs. For this, we rely on a collection of Bioconductor Annotation packages and internal lookup tables.</p> | <p> Gene sets consist of a large variety of gene identifiers, called <strong>original members</strong> here. To use gene sets by GSEA and other querying tools, original members have to be converted to a common universal kind of gene identifiers. Previous releases of MSigDB used human gene symbols for this purpose. Researchers prefer working with gene symbols because they can easily recognize, remember and put them in the context of their work. Unfortunately, a gene usually has multiple different symbols. Conversely, the same symbol can often refer to a number of different genes. Finally, gene symbols change frequently. To overcome these issues, we chose Entrez Gene IDs as robust universal identifiers (called <strong>ezid members</strong> here. Entrez Gene IDs uniquely identify human genes and never change. For convenience, we continue displaying gene sets as made from human gene symbols. However, the symbols are now unambiguously derived from the corresponding human Entrez Gene IDs. For non-human original members, we first convert them to the organism-specific Entrez Gene IDs and then seek their orthologous counterparts as human Entrez Gene IDs. For this, we rely on a collection of Bioconductor Annotation packages and internal lookup tables.</p> | ||
| − | <p>We have updated gene symbols for all sets and collections according to gene_info.gz and gene_history.gz files downloaded from <a href="http://www.ncbi.nlm.nih.gov/gene">Entrez Gene</a> FTP site on November, 15, 2011.</p> | + | <p>We have updated gene symbols for all sets and collections according to <tt>gene_info.gz</tt> and <tt>gene_history.gz</tt> files downloaded from <a href="http://www.ncbi.nlm.nih.gov/gene">Entrez Gene</a> FTP site on November, 15, 2011.</p> |
<h2>Gene Sets Update</h2> | <h2>Gene Sets Update</h2> | ||
Revision as of 14:54, 2 October 2012
This page describes changes in Release 3.1 of the Molecular Signatures Database (MSigDB)
Contents
Gene Symbols Update
Gene sets consist of a large variety of gene identifiers, called original members here. To use gene sets by GSEA and other querying tools, original members have to be converted to a common universal kind of gene identifiers. Previous releases of MSigDB used human gene symbols for this purpose. Researchers prefer working with gene symbols because they can easily recognize, remember and put them in the context of their work. Unfortunately, a gene usually has multiple different symbols. Conversely, the same symbol can often refer to a number of different genes. Finally, gene symbols change frequently. To overcome these issues, we chose Entrez Gene IDs as robust universal identifiers (called ezid members here. Entrez Gene IDs uniquely identify human genes and never change. For convenience, we continue displaying gene sets as made from human gene symbols. However, the symbols are now unambiguously derived from the corresponding human Entrez Gene IDs. For non-human original members, we first convert them to the organism-specific Entrez Gene IDs and then seek their orthologous counterparts as human Entrez Gene IDs. For this, we rely on a collection of Bioconductor Annotation packages and internal lookup tables.
We have updated gene symbols for all sets and collections according to gene_info.gz and gene_history.gz files downloaded from <a href="http://www.ncbi.nlm.nih.gov/gene">Entrez Gene</a> FTP site on November, 15, 2011.
Gene Sets Update
The following describes the changes made to the gene set collections for MSigDB v3.1.
Size Filtering
After mapping to human Entrez Gene IDs, sets with fewer than five genes were not included in the v3.0 release.
C1: Positional gene sets (-60)
- 60 gene sets have been deprecated due to small size (fewer than five human Entrez Gene IDs).
No other changes were made in the C1 gene sets. For a description of this collection, see the <a href="http://www.broad.mit.edu/gsea/msigdb/collections.jsp">Browse Collections</a> page.
C2: Curated gene sets (+1,380)
The C2 collection consists of gene sets collected from various sources such as online pathway databases, publications in PubMed, and knowledge of domain experts. Gene sets in this collection have been extensively revised and expanded.
Note that all the gene set names for C2 have changed. Many of the names used in v2.5 were confusing or wrong, so these have been clarified or corrected. For CGP, the new naming convention is that all gene set names begin with the surname of the first author of the source paper. For CP, the names now begin with the contributor organization.
- CGP: chemical and genetic perturbations (2,392 gene sets). See <a href="http://www.broadinstitute.org/cancer/software/gsea/wiki/index.php/Msigdb_mapping_v2.5_to_v3">this page</a> for information about MSigDB 2.5 gene sets that have been renamed, retired, recombined, or replaced in the MSigDB 3.0 release. All these gene sets have been verified against the original sources. During the reviewing process, we have:
- renamed gene sets to follow consistent conventions throughout the whole collection
- wrote new, enhanced, brief descriptions according to consistent conventions throughout the whole collection
- validated and corrected, if necessary, every attribute for each existing gene set
- added the exact source of the gene set (e.g., Table 1)
- added GEO or ArrayExpress ID when available
- added links to human Entrez Gene entries and PubChem Compound entries as appropriate
- used the original gene identifiers as reported in the source paper (not all gene sets did this originally)
- resolved cases of redundant gene sets
- CP: canonical pathways (880 gene sets).
- We have deprecated all gene sets:
- from GenMAPP gene sets because the majority of them in the previous release are based on KEGG or GO information that we already have
- from GO in this collection because they are already represented by C5
- based on NetAffx annotations because these are largely based on GO and thus are already represented by C5
- with untraceable origins
- We have replaced all existing BioCarta and KEGG gene sets with updated versions from these resources.
- In collaboration with Reactome, we added 430 new canonical pathway sets
- To reduce redundancy in canonical pathways from BioCarta, KEGG, and Reactome, we developed and applied the following filters:
- Source priority: KEGG > Reactome > BioCarta
- Size priority: keep the set with the smaller size
- Name length priority: keep the set with the shorter name
- External ID priority: keep the set with the smaller ID (applied to Reactome sets only)
- For convenience, we have organized gene sets from BioCarta, KEGG, and Reactome as separate, third-level divisions within C2:CP
- We have deprecated all gene sets:
C3: Motif gene sets (-1)
- Thanks to a sharp user, we fixed an error in the description of the gene set "V$NRF2_01".
- All uncategorized gene sets in this collection have been assigned to the TFT subcollection.
- One gene set in the MIR subcollection has been deprecated due to small size (less than five human Entrez Gene IDs).
- No other changes were made in the C3 gene sets.
Gene sets in the C3 collection consist of genes sharing a cis-regulatory motif. This collection contains the following two subcollections:
Transcription factor targets (TFT)
These sets share upstream cis-regulatory motifs which can function as potential transcripton factor binding sites. We used two approaches to generate these gene sets.
We extracted 460 mammalian transcriptional regulatory motifs from v7.4 TRANSFAC database. We then generated the motif gene sets consisting of the inferred target genes for each motif. Every such set consists of human genes whose promoters (defined as regions -2kb to +2kb around transcription start site) contain at least one instance of the motif. We named these sets by the corresponding TRANSFAC matrix identifiers, e.g., V$MIF1_01. The set’s full description is the TRANSFAC entry for the matching matrix, in a format described here.
Motif gene sets of ‘conserved instances’ consist of the inferred target genes for each motif m of 174 upstream motifs highly conserved among five mammalian species. The motifs are catalogued in Xie, et al. (2005, Nature 434, 338–345) and represent potential transcription factor binding sites. Each motif gene set consists of all human genes whose promoters (defined as regions -kb to +2kb around transcription start site) contained at least one conserved instance of motif m. If the motif’s sequence matched a transcription factor binding site documented in the TRANSFAC database (see above), then we appended the name of the TRANSFAC binding matrix to the motif sequence in the gene name, e.g.: MOTIFSEQ_FOO, where MOTIFSEQ is the sequence of motif m and FOO is the TRANSFAC matrix name (e.g., V$MIF1_01). The set’s full description in this case is the TRANSFAC entry for the matching matrix. If the motif’s sequence matched no transcription factor binding site from TRANSFAC v.7.4, then we named the set as MOTIFSEQ_UNKNOWN where MOTIFSEQ is the sequence of motif m.
microRNA Targets (MIR)
These gene sets consist of the inferred target gene for each motif m of 221 3'-UTR motifs highly conserved among five mammalian species. The motifs are catalogued catalogued in http://www.ncbi.nlm.nih.gov/pubmed/15735639 Xie, et al. (2005, Nature 434, 338–345)] and represent potential microRNA binding sites. Each motif gene set consists of all genes whose 3’-UTR contained at least one conserved instance of motif m.
C4: Computational gene sets (-2)
- Two gene sets in the CM subcollection have been deprecated due to small size (fewer than five human Entrez Gene IDs).
No other changes were made in the C4 gene sets. For a description of this collection, see the <a href="http://www.broad.mit.edu/gsea/msigdb/collections.jsp">Browse Collections</a> page.
Cancer gene neighborhoods (CGN)
Starting with a curated list of 380 cancer-associated genes (Brentani et al. 2003, Proc. Natl. Acad. Sci. USA 100, 13418-13423), Subramanian, Tamayo et al. (2005 Proc. Natl. Acad. Sci. USA 102, 15545-15550) mined four expression compendia for correlated gene sets. Gene neighborhoods with fewer than 25 genes at a Pearson correlation threshold of 0.8 were omitted yielding 427 sets.
Gene set names indicate the corresponding expression compendia and the seed cancer-associated genes.
C5: Gene Ontology gene sets
- Names of 71 gene sets have been changed by removing pairs of consecutive underscore characters ('_').
No other changes were made in the C5 gene sets. For a description of this collection, see the <a href="http://www.broad.mit.edu/gsea/msigdb/collections.jsp">Browse Collections</a> page.
For more information
For complete descriptions of all collections or to download the updated gene sets, go to the <a href="http://www.broad.mit.edu/gsea/msigdb/collections.jsp">Browse Collections</a> page.
Other Updates
Changes in organism annotation
We applied consistent organism names throughout MSigDB gene sets. For example, 'Human' was renamed to 'Homo sapiens'. Pathways reported as 'Human' which, according to the original records, are not present in 'Homo sapiens' were deprecated. We inspected annotations of every gene set and replaced 'Generic' organism name with the appropriate correct organism name. For genes in the C3 collection, we replaced 'Human,Mouse,Rat,Dog' organism name with the more correct 'Homo sapiens' name. While the majority of motifs themselves are conserved across the five species, human genome sequences were used to identify genes containing these and other motifs.
Changes in the database XML file format
We have added new attributes to the XML database file. See <a href="http://www.broadinstitute.org/cancer/software/gsea/wiki/index.php/MSigDB_XML_description">this page</a> for more information.
More options with GMT files
Traditionally, all gene sets in MSigDB use human gene symbols as standard gene identifiers. GMT files with gene sets based on human gene symbols contain the word symbols in their names. For standard GSEA analysis, no change is expected: just continue using *.symbols.gmt files as previously. For the first time, this release contains additional kinds of GMT files, with gene sets using:
- original gene identifiers: exactly as they appear in the sources of gene sets; these files contain the word orig in their names
Since original identifiers are from a variety of platforms, we do not recommend using them for routine GSEA analysis. Instead, these files should serve for uses other than standard GSEA analysis. - human Entrez Gene IDs; these files contain the word entrez in their names
Because these files are not identical to the standard *.symbols.gmt files, GSEA done with these files might produce slightly different results compared to the standard *.symbols.gmt files.
For standard GSEA analysis, no change is expected: just continue using .*.symbols.gmt files as previously. Alternatively, you can now choose GMT files with genes rendered as human Entrez Gene IDs.
Compute Overlaps Error Corrected
Several users noted an error in calculations of p-values in the Compute Overlaps tool at the MSigDB web site. Thanks to these reports, we have corrected the error.
MSigDB v2.5 Files
The MSigDB v2.5 files are archived and are still available for download on the <a href="http://www.broadinstitute.org/gsea/downloads.jsp">Downloads</a> page. To load them into GSEA, see the GSEA 2.0.7 Release Notes.
