SringInteractome

From Wikili
Jump to: navigation, search

StringInteractome est intégré dans Gscope.

A partir d'une liste de gènes, StringInteractoime interroge la base de données Postgresql String et crée 3 réseaux d'interactions (WithAny, Sandwich et QueryOnly.

Il crée tous les fichiers nécessaires à Cytoscape et surtout en web à CytoscapeWeb qui sait afficher les réseaux en interactif dans le navigateur


Utilisation en ligne de commande

Il suffit de donner un fichier contenant une liste de gènes ou de probesets (NM, affy, nom de gène, etc.), par exemple le fichier String_essai.txt avec une ligne d'entête (s'il n'y en a pas on mettra data1 data2, etc...). NB: la ligne d'entête doit commencer par #

 #accessnum	genename
 NM_177470	Acaa2
 NM_010678	Aff3


setgscoperr
gscope puts StringInteractome {FichierDepart RepertoireParent Org {SeuilString ""} {From ""} {WithoutTextMiningOnly ""}}
# (for example, see details described below :)
gscope puts StringInteractome inputFile . Homo_sapiens 700

ou directement en http (ici je donne un exemple avec RepertoireParent qui est '/ici/par/la/SameAsFile' (SameAsFile tel quel dans le texte), et l'orga Homo_sapiens (avec '_' merci)

http://lbgi.igbmc.fr/~ripp/cgi-bin/GscopeServer?EVImm&StringInteractome&/home/ripp/lg.lst&/ici/par/la/SameAsFile&Homo_sapiens
  • FichierDepart (est ici String_essai.txt)
  • RepertoireParent va contenir un nouveau répertoire
    • appelé String_taxid_date qui est créé automatiquement
    • ou si RepertoireParent est de la forme /bidule/machin/SameAsFile on créera /bidule/machin/FichierDepart (sans .txt ... ici String_essai)
  • Org est soit Mus_musculus ou "Mus musculus" ou 10090, Homo_sapiens ou "Homo sapiens" ou 9606, Rattus_norvegicus etc.
  • SeuilString par défaut 700
  • From est FromNM, FromAffy ou FromTest (défaut). Si FromTest on teste chaque valeur pour détecter si c'est NM ou Affy etc.
  • WithoutTextMiningOnly on ne prend pas ceux qui ne sont trouvés QUE PAR textmining

StringInteractome va créer le répertoire RepertoireParent/String_10090_20091213_1655 (ou RepertoireParent/String_essai) qui contiendra

 AttributeNode		contient PreferredName, StringId, et toutes les colonnes du FichierDepart. Permet de colorer les noeuds dans Cytoscape selon les colonnes de ce fichier (exemple: colorer selon l'expression). Par défaut, contient une colonne appelée QueryColor qui permettra de colorer les gènes de départ
 Identifiers		contient les identifiants String correspondant aux gènes du FichierDepart
 QueryOnly.string	network avec uniquement les gènes du FichierDepart
 Sandwich.string	network avec les gènes du FichierDepart et ceux possédant au moins deux interactions avec nos gènes de départ
 String_essai.tgz	le tar gzip de tout ça.
 String_essai.txt	le FichierDepart
 String_essai_context.txt le contexte dans lequel on a créé cette liste
 Summary.txt		contient pour chaque query: Query;StringId;PreferredName;toutes les colonnes du FichierDepart
 WithAny.parser         contient pour chaque query: le nombre d'interactions, le nombre d'interactions uniquement avec les gènes du FichierDepart, le % d'interactions uniquement avec les gènes du FichierDepart et la liste de ces gènes
 WithAny.string         network avec les gènes du FichierDepart et toutes leurs interactions

PHP StringInteractome

tout ça est accessible depuis PHP.

  • par RrFun::Wscope("EVImm", "StringInteractome", $params);
 params est un array de la forme (on voit ici les valeurs prises par défaut)
       $params["org"]       = "Mus_musculus";  (puisque Laëtitia refuse de mettre un défaut dans StringInteractome !!!)
       $params["seuil"]     = "";
       $params["from"]      = "";
       $params["wtmo"]      = "";
       $params["repparent"] = "home/ripp/GenoretWww/ImAnno/ImAnnoEar/String";

Dans touts les cas on crée un tar gzip .tgz qui peut être téléchargé.

Cytoscape

On a maintenant tout ce qu'il faut pour entrer dans Cytoscape

  1. File/import network from table WithAny.string
  2. File/import network from table Sandwich.string
  3. File/import network from table QueryOnly.string
  4. import attribute AttributeNode

Quand on fait import network il faut :

  1. select file ...
  2. cocher Show text file import options
    • Delimiter Tab (enlever Space)
    • Show all entries
    • Transfert first line as attribute names
  1. Source Interaction prendre column 2 (PNInitial)
  2. Target Interaction prendre column 4 (PNCopain)
  3. Import
  4. Close la fenêtre de baratin

Quand on fait import Attributes

  1. Select table ...
  2. Show text file import options
    • Delimiter Tab (c'est déjà bon)
    • Show all entries
    • Transfert first line as attribute names
  1. Import
  2. Close la fenêtre de baratin

Avec ça on peut colorier nos gènes

  1. Cliquer sur l'onglet WizMapper
  2. Double-cliquer Node Color (il saute en haut)
  3. Select Value QueryColor
  4. Mapping type Discrete Mapper
  5. Cliquer sur le 1 puis dans sa case vide à droite
  6. Il apparaît trois petits points cliquer dessus pour choisir une couleur

Bonne chance !

              Raymond