void gdfFreeGeneList(struct gdfGene **pList)
/* Free a whole list of genes. */
{
struct gdfGene *gene, *next;
gene = *pList;
while (gene != NULL)
{
next = gene->next;
gdfFreeGene(gene);
gene = next;
}
*pList = NULL;
}
boolean getWormGeneExonDna(char *name, DNA **retDna)
/* Get the DNA associated with a gene, without introns. */
{
struct gdfGene *g;
struct slName *syn = NULL;
long lstart, lend;
int start, end;
int dnaSize;
DNA *dna;
int i;
struct gdfDataPoint *pt = NULL;
struct wormGdfCache *gdfCache;
struct dyString *dy = newDyString(1000);
/* Translate biologist type name to cosmid.N name */
if (wormIsGeneName(name))
{
syn = wormGeneToOrfNames(name);
if (syn != NULL)
name = syn->name;
}
if (strncmp(name, "g-", 2) == 0)
gdfCache = &wormGenieGdfCache;
else
gdfCache = &wormSangerGdfCache;
if ((g = wormGetSomeGdfGene(name, gdfCache)) == NULL)
return FALSE;
gdfGeneExtents(g, &lstart, &lend);
start = lstart;
end = lend;
/*wormClipRangeToChrom(chromIds[g->chromIx], &start, &end);*/
dnaSize = end-start;
dna = wormChromPart(chromIds[g->chromIx], start, dnaSize);
gdfOffsetGene(g, -start);
if (g->strand == '-')
{
reverseComplement(dna, dnaSize);
gdfRcGene(g, dnaSize);
}
pt = g->dataPoints;
for (i=0; i<g->dataCount; i += 2)
{
dyStringAppendN(dy, (dna+pt[i].start), (pt[i+1].start - pt[i].start));
}
*retDna = cloneString(dy->string);
dyStringFree(&dy);
gdfFreeGene(g);
return TRUE;
}
boolean getWormGeneDna(char *name, DNA **retDna, boolean upcExons)
/* Get the DNA associated with a gene. Optionally upper case exons. */
{
struct gdfGene *g;
struct slName *syn = NULL;
long lstart, lend;
int start, end;
int dnaSize;
DNA *dna;
struct wormGdfCache *gdfCache;
/* Translate biologist type name to cosmid.N name */
if (wormIsGeneName(name))
{
syn = wormGeneToOrfNames(name);
if (syn != NULL)
name = syn->name;
}
if (strncmp(name, "g-", 2) == 0)
gdfCache = &wormGenieGdfCache;
else
gdfCache = &wormSangerGdfCache;
if ((g = wormGetSomeGdfGene(name, gdfCache)) == NULL)
return FALSE;
gdfGeneExtents(g, &lstart, &lend);
start = lstart;
end = lend;
/* wormClipRangeToChrom(chromIds[g->chromIx], &start, &end); */
dnaSize = end-start;
*retDna = dna = wormChromPart(chromIds[g->chromIx], start, dnaSize);
gdfOffsetGene(g, -start);
if (g->strand == '-')
{
reverseComplement(dna, dnaSize);
gdfRcGene(g, dnaSize);
}
if (upcExons)
{
int i;
struct gdfDataPoint *pt = g->dataPoints;
for (i=0; i<g->dataCount; i += 2)
{
toUpperN(dna + pt[i].start, pt[i+1].start - pt[i].start);
}
}
gdfFreeGene(g);
return TRUE;
}
DNA *wormChromPartExonsUpper(char *chromId, int start, int size)
/* Return part of a worm chromosome with exons in upper case. */
{
DNA *dna = wormChromPart(chromId, start, size);
struct wormFeature *geneFeat = wormGenesInRange(chromId, start, start+size);
struct wormFeature *feat;
for (feat = geneFeat; feat != NULL; feat = feat->next)
{
char *name = feat->name;
if (!wormIsNamelessCluster(name))
{
struct gdfGene *gene = wormGetGdfGene(name);
gdfUpcExons(gene, feat->start, dna, size, start);
gdfFreeGene(gene);
}
}
slFreeList(&geneFeat);
return dna;
}
