4
$\begingroup$

I currently have two sets of gene expression data. The first is a dataframe of genes identified by Annotation id CG numbers (for example "CG10005") in one column and a numerical variable of interest associated with each of those genes in the other (the important part is that the genes are identified by the CG number). The second dataset is this study from multiple species where genes are identified by "GLEANR" ID's like this one: "dsim_GLEANR_10060".

I would like to extract all of the values from this study for the D.melanogaster and D.simulans data, then link that to the Annotation ID's in the first dataset (to connect my variable of interest to the sex bias in expression for each species). The trouble I am having is connecting the orthologous genes within the second study to each other, and then connecting them to the CG identifiers in the first study.

Does anyone have any suggestions on resources which can be used to get this all to link up?

Improve this question
$\endgroup$

2 Answers 2

Reset to default
2
$\begingroup$

This is the FlyBase page for the example gene: Dsim\GD10095. There, you have a section "orthologs", linking to OrthoDB. So my suggestion is: Find the list of synonyms for D. simulans on FlyBase (perhaps here?), download the Drosophila section of OrthoDB, and finally find the 1:1 orthologs.

Improve this answer
$\endgroup$
3
  • $\begingroup$ Thanks, I am looking at this now - I assume each gene has a unique FBgn and their orthologs are identified by the ODB6_OG_ID (?) and the species identified by DROME, DROSI, DROSE... taking just those three variables I can create a data frame of each ODB6_OG_ID in the first column and columns for each species giving its orthologous FBgn ID for each species of interest? But how do I "Find the 1:1 orthologs"? $\endgroup$
    – rg255
    Commented Mar 8, 2013 at 10:30
  • $\begingroup$ Sometimes you have a strain-specific gene duplication, which you probably want to exclude to keep things simple. So you can check for each species if there's more than one gene associated with an ODB6_OG_ID, and remove those. $\endgroup$
    – Michael Kuhn
    Commented Mar 8, 2013 at 11:18
  • $\begingroup$ I got this all working after a good ol' fashioned wrestle with R today - thanks a lot for the help! $\endgroup$
    – rg255
    Commented Mar 11, 2013 at 20:02
0
$\begingroup$

You could also use biomart for this purposes.

Click on the "ID Converter" tab in the "Tools" section (on the left).

Improve this answer
$\endgroup$

You must log in to answer this question.

Not the answer you're looking for? Browse other questions tagged .