83dbcadd73270e8b652fb4b5537235ca5727c7fa hiram Mon Dec 8 16:22:38 2025 -0800 fix the toga click code to use a simple local hideToggle class refs #35776 diff --git src/hg/hgc/togaClick.c src/hg/hgc/togaClick.c index 123f3246292..f446ec6bb5b 100644 --- src/hg/hgc/togaClick.c +++ src/hg/hgc/togaClick.c @@ -1,742 +1,734 @@ /* togaClick - click handling for TOGA tracks */ #include "common.h" #include "hgc.h" #include "togaClick.h" #include "string.h" #include "htmshell.h" #include "chromAlias.h" struct togaDataBB *togaDataBBLoad(char **row, bits16 fieldCount) /* Load a togaData from row fetched with select * from togaData * from database. Dispose of this with togaDataFree(). */ { struct togaDataBB *ret; AllocVar(ret); ret->projection = cloneString(row[0]); ret->ref_trans_id = cloneString(row[1]); ret->ref_region = cloneString(row[2]); ret->query_region = cloneString(row[3]); ret->chain_score = cloneString(row[4]); ret->chain_synteny = cloneString(row[5]); ret->chain_flank = cloneString(row[6]); ret->chain_gl_cds_fract = cloneString(row[7]); ret->chain_loc_cds_fract = cloneString(row[8]); ret->chain_exon_cov = cloneString(row[9]); ret->chain_intron_cov = cloneString(row[10]); ret->status = cloneString(row[11]); ret->perc_intact_ign_M = cloneString(row[12]); ret->perc_intact_int_M = cloneString(row[13]); ret->intact_codon_prop = cloneString(row[14]); ret->ouf_prop = cloneString(row[15]); ret->mid_intact = cloneString(row[16]); ret->mid_pres = cloneString(row[17]); ret->prot_alignment = cloneString(row[18]); ret->svg_line = cloneString(row[19]); ret->ref_link = cloneString(row[20]); ret->inact_mut_html_table = cloneString(row[21]); ret->exon_ali_html = cloneString(row[22]); /* read two optional new items, CDSseq and frame-corrected protein */ ret->CDSseq = NULL; if (fieldCount >= 35) /* 0-11 are bed core fields. This fct gets a pointer starting with element 12. 12+23 = 35 */ ret->CDSseq = cloneString(row[23]); ret->protseqFrameCorrected = NULL; if (fieldCount >= 36) /* 12+24 = 36 */ ret->protseqFrameCorrected = cloneString(row[24]); /* number and percent of mutated exons (additional data now shown for the gene loss status) */ ret->numExonsMutated = NULL; ret->percentExonsMutated = NULL; if (fieldCount >= 37) /* 12+25 = 37 */ ret->numExonsMutated = cloneString(row[25]); if (fieldCount >= 38) /* 12+26 = 38 */ ret->percentExonsMutated = cloneString(row[26]); return ret; } struct togaData *togaDataLoad(char **row) /* Load a togaData from row fetched with select * from togaData * from database. Dispose of this with togaDataFree(). */ { struct togaData *ret; AllocVar(ret); ret->projection = cloneString(row[0]); ret->ref_trans_id = cloneString(row[1]); ret->ref_region = cloneString(row[2]); ret->query_region = cloneString(row[3]); ret->chain_score = cloneString(row[4]); ret->chain_synteny = cloneString(row[5]); ret->chain_flank = cloneString(row[6]); ret->chain_gl_cds_fract = cloneString(row[7]); ret->chain_loc_cds_fract = cloneString(row[8]); ret->chain_exon_cov = cloneString(row[9]); ret->chain_intron_cov = cloneString(row[10]); ret->status = cloneString(row[11]); ret->perc_intact_ign_M = cloneString(row[12]); ret->perc_intact_int_M = cloneString(row[13]); ret->intact_codon_prop = cloneString(row[14]); ret->ouf_prop = cloneString(row[15]); ret->mid_intact = cloneString(row[16]); ret->mid_pres = cloneString(row[17]); ret->prot_alignment = cloneString(row[18]); ret->svg_line = cloneString(row[19]); return ret; } void togaDataBBFree(struct togaDataBB **pEl) /* Free a single dynamically allocated togaDatasuch as created * with togaDataLoad(). */ { struct togaDataBB *el; if ((el = *pEl) == NULL) return; freeMem(el->projection); freeMem(el->ref_trans_id); freeMem(el->ref_region); freeMem(el->query_region); freeMem(el->chain_score); freeMem(el->chain_synteny); freeMem(el->chain_flank); freeMem(el->chain_gl_cds_fract); freeMem(el->chain_loc_cds_fract); freeMem(el->chain_exon_cov); freeMem(el->chain_intron_cov); freeMem(el->status); freeMem(el->perc_intact_ign_M); freeMem(el->perc_intact_int_M); freeMem(el->intact_codon_prop); freeMem(el->ouf_prop); freeMem(el->mid_intact); freeMem(el->mid_pres); freeMem(el->prot_alignment); freeMem(el->svg_line); freeMem(el->ref_link); freeMem(el->inact_mut_html_table); freeMem(el->exon_ali_html); freez(pEl); } void togaDataFree(struct togaData **pEl) /* Free a single dynamically allocated togaDatasuch as created * with togaDataLoad(). */ { struct togaData *el; if ((el = *pEl) == NULL) return; freeMem(el->projection); freeMem(el->ref_trans_id); freeMem(el->ref_region); freeMem(el->query_region); freeMem(el->chain_score); freeMem(el->chain_synteny); freeMem(el->chain_flank); freeMem(el->chain_gl_cds_fract); freeMem(el->chain_loc_cds_fract); freeMem(el->chain_exon_cov); freeMem(el->chain_intron_cov); freeMem(el->status); freeMem(el->perc_intact_ign_M); freeMem(el->perc_intact_int_M); freeMem(el->intact_codon_prop); freeMem(el->ouf_prop); freeMem(el->mid_intact); freeMem(el->mid_pres); freeMem(el->prot_alignment); freeMem(el->svg_line); freez(pEl); } struct togaNucl *togaNuclLoad(char **row) /* Load a togaNucl from row fetched with select * from togaNucl * from database. Dispose of this with togaNuclFree(). */ { struct togaNucl *ret; AllocVar(ret); ret->transcript = cloneString(row[0]); ret->exon_num = cloneString(row[1]); ret->exon_region = cloneString(row[2]); ret->pid = cloneString(row[3]); ret->blosum = cloneString(row[4]); ret->gaps = cloneString(row[5]); ret->ali_class = cloneString(row[6]); ret->exp_region = cloneString(row[7]); ret->in_exp_region = cloneString(row[8]); ret->alignment = cloneString(row[9]); return ret; } void togaNuclFree(struct togaNucl **pEl) /* Free a single dynamically allocated togaNucl such as created * with togaNuclLoad(). */ { struct togaNucl *el; if ((el = *pEl) == NULL) return; freeMem(el->transcript); freeMem(el->exon_num); freeMem(el->exon_region); freeMem(el->pid); freeMem(el->blosum); freeMem(el->gaps); freeMem(el->ali_class); freeMem(el->exp_region); freeMem(el->in_exp_region); freeMem(el->alignment); freez(pEl); } struct togaInactMut *togaInactMutLoad(char **row) /* Load a togaInactMut from row fetched with select * from togaInactMut * from database. Dispose of this with togaInactMutFree(). */ { struct togaInactMut *ret; AllocVar(ret); ret->transcript = cloneString(row[0]); ret->exon_num = cloneString(row[1]); ret->position = cloneString(row[2]); ret->mut_class = cloneString(row[3]); ret->mutation = cloneString(row[4]); ret->is_inact = cloneString(row[5]); ret->mut_id = cloneString(row[6]); return ret; } void togaInactMutFree(struct togaInactMut **pEl) /* Free a single dynamically allocated togaInactMut such as created * with togaInactMutLoad(). */ { struct togaInactMut *el; if ((el = *pEl) == NULL) return; freeMem(el->transcript); freeMem(el->exon_num); freeMem(el->position); freeMem(el->mut_class); freeMem(el->mutation); freeMem(el->is_inact); freeMem(el->mut_id); freez(pEl); } void extractHLTOGAsuffix(char *suffix) /* Extract suffix from TOGA table name. Prefix must be HLTOGAannot */ { int suff_len = strlen(suffix); if (suff_len <= HLTOGA_BED_PREFIX_LEN) // we cannot chop first PREFIX_LEN characters { // TODO: NOT SURE IF IT WORKS; but this must not happen char empty[5] = { '\0' }; strcpy(suffix, empty); } else { // just start the string 11 characters upstream memmove(suffix, suffix + HLTOGA_BED_PREFIX_LEN, suff_len - HLTOGA_BED_PREFIX_LEN + 1); } } void HLprintQueryProtSeqForAli(char *proteinAlignment) { // take protein sequence alignment // print only the query sequence char *str = proteinAlignment; int printed_char_num = 0; while ((str = strstr(str, "que:")) != NULL) { str += 10; char ch; while ((ch = *str++) != '<') { if (ch != '-') { putchar(ch); ++printed_char_num; } if (printed_char_num == 80) { printed_char_num = 0; printf("
"); } } } } void print_with_newlines(const char *str) { int line_length = 80; // Number of characters per line int length = strlen(str); int i = 0; while (i < length) { /* Print up to 80 characters or the remainder of the string */ int chars_to_print = (length - i < line_length) ? (length - i) : line_length; printf("%.*s
", chars_to_print, &str[i]); i += chars_to_print; } } void doHillerLabTOGAGeneBig(char *database, struct trackDb *tdb, char *item, char *table_name) /* Put up TOGA Gene track info. */ // To think about -> put into a single bigBed // string: HTML formatted inact mut // string: HTML formatted exon ali section { int start = cartInt(cart, "o"); int end = cartInt(cart, "t"); char *chrom = cartString(cart, "c"); char *fileName = bbiNameFromSettingOrTable(tdb, NULL, tdb->table); struct bbiFile *bbi = bigBedFileOpenAlias(hReplaceGbdb(fileName), chromAliasFindAliases); struct lm *lm = lmInit(0); struct bigBedInterval *bbList = bigBedIntervalQuery(bbi, chrom, start, end, 0, lm); struct bigBedInterval *bb; char *fields[bbi->fieldCount]; for (bb = bbList; bb != NULL; bb = bb->next) { if (!(bb->start == start && bb->end == end)) continue; // our names are unique char *name = cloneFirstWordByDelimiterNoSkip(bb->rest, '\t'); boolean match = (isEmpty(name) && isEmpty(item)) || sameOk(name, item); if (!match) continue; char startBuf[16], endBuf[16]; bigBedIntervalToRow(bb, chrom, startBuf, endBuf, fields, bbi->fieldCount); break; } printf("

Projection v2 %s

\n", item); struct togaDataBB *info = togaDataBBLoad(&fields[11], bbi->fieldCount); // Bogdan: why 11? 0-11 are bed-like fields likely printf("Reference transcript: %s
", info->ref_link); printf("Genomic locus in reference: %s
\n", info->ref_region); printf("Genomic locus in query: %s
\n", info->query_region); printf("Projection classification: %s
\n", info->status); printf("Probability that query locus is orthologous: %s
\n", info->chain_score); // list of chain features (for orthology classification) -printf("Show features used for ortholog probability\n"); -printf("
\n"); -printf("

\n"); -printf("
\nFeature description:\n"); +printf("

\nFeature description:\n"); printf("For each projection (one reference transcript and one overlapping chain),\n"); printf("TOGA computes the following features by intersecting the reference coordinates of aligning\n"); printf("blocks in the chain with different gene parts (coding exons, UTR (untranslated region) exons, introns)\n"); -printf("and the respective intergenic regions.\n
\n"); +printf("and the respective intergenic regions.\n

\n"); -printf("We define the following variables:\n

\n"); +printf("

Using these variables, TOGA computes the following features:\n"); printf("

\n"); +printf("

\n"); -printf("\n
\n"); +printf("\n"); printf("
\n"); // show inact mut plot printf("

Visualization of inactivating mutations on exon-intron structure

\n"); printf("%s\n", info->svg_line); printf("
Exons shown in grey are missing (often overlap assembly gaps).\nExons shown in"); printf(" red or blue are deleted or do not align at all.\nRed indicates that the exon deletion "); printf("shifts the reading frame, while blue indicates that exon deletion(s) are framepreserving.
\n"); // GLP features -printf("Show features used for transcript classification\n"); -printf("
\n"); -printf("

\n
\n"); printf("
\n"); printf("

Predicted protein sequence

\n"); -printf("Show protein sequence of query\n"); -printf("
\n"); -// printf("{protein seq of the query without dashes or other things. Should end with *}\n"); -printf(""); +printf("Show protein sequence of query\n"); +printf("
\n"); +printf("\n

\n

\n"); if (info->protseqFrameCorrected != NULL) { - printf("Show frame-corrected protein sequence of query (potential frameshifts are masked)\n"); - printf("
\n"); - printf(""); + printf("Show frame-corrected protein sequence of query (potential frameshifts are masked)\n"); + printf("\n"); + printf("\n

\n
\n"); } if (info->CDSseq != NULL) { printf("
\n"); printf("

Predicted coding (DNA) sequence

\n"); - printf("Show coding sequence of query\n"); - printf("
\n"); - printf(""); + printf("Show coding sequence of query\n"); + printf("\n"); + printf("\n

\n
\n"); } // and show protein sequence printf("
\n"); printf("

Protein sequence alignment

\n"); -printf("Show alignment between reference and query\n"); -printf("
\n"); -printf("%s\n", info->prot_alignment); -printf("
\n"); +printf("Show alignment between reference and query\n"); +printf("\n"); // show inactivating mutations if required printf("
\n"); printf("

List of inactivating mutations

\n"); -printf("Show inactivating mutations\n"); -printf("
\n"); -printf("\n"); // init table +printf("Show inactivating mutations\n"); +printf("
\n"); // init table printf("\n"); printf("\n"); printf("%s\n", info->inact_mut_html_table); -printf("
Exon numberCodon numberMutation classMutationTreated as inactivatingMutation ID
\n"); +printf("

\n"); printf("
\n\n"); // show exons data printf("
\n"); printf("

Exon alignments

\n"); -printf("Show exon sequences and features\n"); -printf("
\n"); -// printf("%s\n", info->exon_ali_string); -printf("%s\n", info->exon_ali_html); +printf("Show exon sequences and features\n"); +printf("\n

\n"); -// TODO: check whether I need this -hPrintf(""); -hPrintf(""); -hPrintf(""); - - printTrackHtml(tdb); // and do I need this? } void doHillerLabTOGAGene(char *database, struct trackDb *tdb, char *item, char *table_name) /* Put up TOGA Gene track info. */ { //int start = cartInt(cart, "o"); char headerTitle[512]; char suffix[512]; strcpy(suffix, table_name); extractHLTOGAsuffix(suffix); safef(headerTitle, sizeof(headerTitle), "%s", item); genericHeader(tdb, headerTitle); printf("

TOGA gene annotation

\n"); // htmlHorizontalLine(); if (startsWith("bigBed", tdb->type)) { doHillerLabTOGAGeneBig(database, tdb, item, table_name); return; } struct sqlConnection *conn = hAllocConn(database); // define TOGA table names: initate with pre-defined prefixes char togaDataTableName[256]; char togaNuclTableName[256]; char togaInactMutTableName[256]; strcpy(togaDataTableName, HLTOGA_DATA_PREFIX); strcpy(togaNuclTableName, HLTOGA_NUCL_PREFIX); strcpy(togaInactMutTableName, HLTOGA_INACT_PREFIX); // add suffix strcat(togaDataTableName, suffix); strcat(togaNuclTableName, suffix); strcat(togaInactMutTableName, suffix); if (hTableExists(database, togaDataTableName)) { printf("

Projection %s


\n", item); char query[256]; struct sqlResult *sr = NULL; char **row; struct togaData *info = NULL; sqlSafef(query, sizeof(query), "select * from %s where transcript='%s'", togaDataTableName, item); sr = sqlGetResult(conn, query); if ((row = sqlNextRow(sr)) != NULL) { info = togaDataLoad(row); // parse sql output // fill HTML template: printf("Reference transcript: %s
", info->ref_trans_id, info->ref_trans_id); printf("Genomic locus in reference: %s
\n", info->ref_region); printf("Genomic locus in query: %s
\n", info->query_region); printf("Projection classification: %s
\n", info->status); printf("Probability that query locus is orthologous: %s
\n", info->chain_score); // list of chain features (for orthology classification) - printf("Show features used for ortholog probability\n"); - printf("
\n"); - printf("
    \n"); + printf("Show features used for ortholog probability\n"); + printf("

\n"); printf("
\nFeature description:\n"); printf("For each projection (one reference transcript and one overlapping chain),\n"); printf("TOGA computes the following features by intersecting the reference coordinates of aligning\n"); printf("blocks in the chain with different gene parts (coding exons, UTR (untranslated region) exons, introns)\n"); printf("and the respective intergenic regions.\n
\n"); printf("We define the following variables:\n
    \n"); printf("
  • c: number of reference bases in the intersection between chain blocks and coding exons of the gene under consideration.
    • \n"); printf("
  • C: number of reference bases in the intersection between chain blocks and coding exons of all genes.
    • \n"); printf("
  • a: number of reference bases in the intersection between chain blocks and coding exons and introns of the gene under consideration.
    • \n"); printf("
  • A: number of reference bases in the intersection between chain blocks and coding exons and introns of all genes and the intersection\n"); printf("between chain blocks and intergenic regions (excludes UTRs).
    • \n"); printf("
  • f: number of reference bases in chain blocks overlapping the 10 kb flanks of the gene under consideration.\n"); printf("Alignment blocks overlapping exons of another gene that is located in these 10 kb flanks are ignored.
    • \n"); printf("
  • i: number of reference bases in the intersection between chain blocks and introns of the gene under consideration.
    • \n"); printf("
  • CDS (coding sequence): length of the coding region of the gene under consideration.
    • \n"); printf("
  • I: sum of all intron lengths of the gene under consideration.
    • \n"); printf("
\n"); printf("Using these variables, TOGA computes the following features:\n"); printf("
    \n"); printf("
  • "global CDS fraction" as C / A. Chains with a high value have alignments that largely overlap coding exons,"); printf("which is a hallmark of paralogous or processed pseudogene chains. In contrast, chains with a low value also align many "); printf("intronic and intergenic regions, which is a hallmark of orthologous chains.
    • \n"); printf("
  • "local CDS fraction" as c / a. Orthologous chains tend to have a lower value, as intronic "); printf("regions partially align. This feature is not computed for single-exon genes.
    • \n"); printf("
  • "local intron fraction" as i / I. Orthologous chains tend to have a higher value."); printf("This feature is not computed for single-exon genes.
    • \n"); printf("
  • "flank fraction" as f / 20,000. Orthologous chains tend to have higher values,"); printf("as flanking intergenic regions partially align. This feature is important to detect orthologous loci of single-exon genes.
    • \n"); printf("
  • "synteny" as log10 of the number of genes, whose coding exons overlap by at least one base aligning"); printf("blocks of this chain. Orthologous chains tend to cover several genes located in a conserved order, resulting in higher synteny values.
    • \n"); printf("
  • "local CDS coverage" as c / CDS, which is only used for single-exon genes.
    • \n"); printf("
\n"); printf("\n
\n
\n"); htmlHorizontalLine(); // show inact mut plot printf("

Visualization of inactivating mutations on exon-intron structure

\n"); printf("%s
\n", info->svg_line); printf("
Exons shown in grey are missing (often overlap assembly gaps).\nExons shown in"); printf(" red or blue are deleted or do not align at all.\nRed indicates that the exon deletion "); printf("shifts the reading frame, while blue indicates that exon deletion(s) are framepreserving.
\n"); // GLP features - printf("Show features used for transcript classification\n"); - printf("
\n"); - printf("
    \n"); + printf("Show features used for transcript classification\n"); + printf("\n
    \n"); + printf("

\n
\n
\n"); printf("

Query protein sequence


"); - printf("Show protein sequence of query\n"); - printf("
\n"); - printf("{protein seq of the query without dashes or other things. Should end with *}\n"); - printf("
\n\n
\n"); + printf("Show protein sequence of query\n"); + printf("\n"); // and show protein sequence htmlHorizontalLine(); printf("

Protein sequence alignment


\n"); - printf("Show alignment between reference and query\n"); + printf("Show alignment between reference and query\n"); printf("
\n"); - printf("%s
\n", info->prot_alignment); + printf("

%s


\n", info->prot_alignment); printf("
\n

\n"); // do not forget to free toga data struct togaDataFree(&info); } else { // no data found, need to report this - printf("Not found data for %s\n", item); + printf("

No found data for %s

\n", item); } sqlFreeResult(&sr); } // show inactivating mutations if required printf("

List of inactivating mutations


\n"); if (hTableExists(database, togaInactMutTableName)) { char query[256]; struct sqlResult *sr = NULL; char **row; sqlSafef(query, sizeof(query), "select * from %s where transcript='%s'", togaInactMutTableName, item); sr = sqlGetResult(conn, query); printf("Show inactivating mutations\n"); printf("
\n"); printf("\n"); // init table printf("\n"); printf("\n"); while ((row = sqlNextRow(sr)) != NULL) { struct togaInactMut *info = NULL; info = togaInactMutLoad(row); printf("\n"); printf("\n", info->exon_num); printf("\n", info->position); printf("\n", info->mut_class); printf("\n", info->mutation); if (sameWord(info->is_inact, ONE_)){ printf("\n", YES_); } else { printf("\n", NO_); } printf("\n", info->mut_id); printf("\n"); togaInactMutFree(&info); } sqlFreeResult(&sr); printf("
Exon numberCodon numberMutation classMutationTreated as inactivatingMutation ID
%s%s%s%s%s%s%s
\n"); printf("
\n
\n"); } else { printf("Sorry, cannot find TOGAInactMut table.
\n"); } // show exons data htmlHorizontalLine(); printf("

Exon alignments


\n"); if (hTableExists(database, togaNuclTableName)) { char query[256]; struct sqlResult *sr = NULL; char **row; printf("Show exon sequences and features\n"); printf("
\n"); sqlSafef(query, sizeof(query), "select * from %s where transcript='%s'", togaNuclTableName, item); sr = sqlGetResult(conn, query); while ((row = sqlNextRow(sr)) != NULL) { struct togaNucl *info = NULL; info = togaNuclLoad(row); printf("
Exon number: %s

\n", info->exon_num); printf("Exon region: %s
\n", info->exon_region); printf("Nucleotide percent identity: %s | BLOSUM: %s
\n", info->pid, info->blosum); if (sameWord(info->gaps, ONE_)){ printf("Intersects assembly gaps: %s
\n", YES_); } else { printf("Intersects assembly gaps: %s
\n", NO_); } printf("Exon alignment class: %s
\n", info->ali_class); if (sameWord(info->in_exp_region, ONE_)){ printf("Detected within expected region (%s): %s
\n", info->exp_region, YES_); } else { printf("Detected within expected region (%s): %s
\n", info->exp_region, NO_); } // printf("Expected region: %s
\n", info->exp_region); printf("
\n"); printf("Sequence alignment between reference and query exon:
\n"); printf("%s
\n", info->alignment); togaNuclFree(&info); } sqlFreeResult(&sr); printf("
\n

\n"); } else { printf("Sorry, cannot find TOGANucl table.
\n"); } htmlHorizontalLine(); // TODO: check whether I need this printf("%s", hgTracksPathAndSettings()); hPrintf(""); hPrintf(""); hPrintf(""); printTrackHtml(tdb); // and do I need this? hFreeConn(&conn); }