How to run SQL Queries against the database?

As most of the databases in the world, our database uses a special language named SQL 
(Structured Query Language) to execute queries. 

This language allows different operations, but the most common is to retrieve the data. For 
this purpose SELECT statement is used. 

How SELECT works? 

Well, it is quite intuitive once you know the basic things. 

First it is good to know: 

 The name of the table that you want to query. 
 The structure of the table (fields, data types…) 

Let’s start with an example: 

 We want to retrieve all the data included in the table: all_info. The list of the tables 
available is shown in the drop‐down menu Table Info:  

 Select the option "all_info", and the description of the table will appear: 

 
… and scrolling: 

 
This basically means that: 

The table “all_info” has 13 fields (columns) with the names:  

cluster_id    unique cluster identifier (ranges from 1 .. 29608) 
cluster_name     ong name of the cluster (rA_c02_1 ‐ re‐assembly number 02 cluster 1)
cluster_seq  DNA sequence of the cluster 
cluster_length    length of the sequence (in nucleotides) 
P_reads  number of reads in POLYP stage 
S_reads   number of reads in STROBILA stage 
E_reads  number of reads in EPHYRA stage 
pep_seq  predicted peptide sequence (or No Prediction tag) 
pep_length  length of the peptide (in amino acids) 
score  peptide prediction score (given by EstScan) 
start_nuc  position of the first nucleotide of the ORF 
stop_nuc  position of the last nucleotide of the ORF 

The table also indicates the type of the data: numeric (int or smallint) or text (varchar, 
mediumtext). 

Now we are going to construct and run our first query (retrieve all the data from a given 
table). 

Just type the following on the query field:  

SELECT * 
FROM all_info 

Which means: Retrieve (SELECT) all the fields (*) FROM table named: all_info. 

 
NOTE: The syntax is case insensitive, so “SELECT” is the same as “select” or “SElecT”.  
The same happens with the column names and symbols. 

 Press the button   and a table with the results should appear (If not, 
check that you have written the query correctly): 

IMPORTANT: The results are limited to 1000 rows by default. Depending of the query, the 
results can be very large and memory consuming. Type "0" in the field for using no limit. 

Saving results as CSV 

In order to export the result to a file, the button   can be used. This will save 
the results and the filters applied to the table in a CSV file that one can download to the 
computer and open in Ms Excel, for example. 

 
Hiding Columns 
 
Question: How can I retrieve only data from columns cluster_id, P_reads, S_reads, E_reads 
and total_reads, for example?  
 
Answer: Just change the * for the field names separated by commas (,) except the last. 

SELECT cluster_id, P_reads, S_reads, E_reads, total_reads  
FROM all_info 

Ordering Data 

Question: How do I order the data ascendant or descendant? 
Answer: There are two possibilities, but with a slightly difference: 

 Click on the column’s header. This will order the values of this column ascendant or 
descendant, but only affects to the values displayed. Here we sorted all the rows 
according to the expression level in strobila (see small arrowhead near "S_reads"). 
  

 Use the clause ORDER BY in your query and use ASC or DESC to indicate the order. 
This affects to all the values in the database table. 

SELECT cluster_id, P_reads, S_reads, E_reads, total_reads  
FROM all_info 
ORDER BY S_reads DESC 

Filtering Data 

Question: How do I filter the data? 
Answer: There are two possibilities: 

 Use the filter fields on the column header (text or numeric data). But, once time 
more, it only affects to the values displayed. If your result is greater than the limit 
specified (1000 by default) it could be some values on the results table that are not 
showed. 

For numeric values the following operators are available:  

equal to: = N 
greater than: > N 
lesser than: < N 
lesser or equal: <= N 
greater or equal: >= N 
range of values: N1 .. N2 

  
  Use the clause WHERE and then the condition of your filter (it allows more complex 
searches). This is the recommended way: 

SELECT cluster_id, P_reads, S_reads, E_reads, total_reads  
FROM all_info 
WHERE S_reads>1000 and E_reads<1000 

 Example with arithmetical and logical operations: 

SELECT cluster_id, P_reads, S_reads, E_reads, total_reads 
FROM all_info 
WHERE (S_reads+1) / (P_reads+1) >= 500 

 
 SELECT cluster_id, P_reads, S_reads, E_reads, total_reads 
FROM all_info 
WHERE (S_reads+1) / (P_reads+1) >= 500 and total_reads>2000 

 Adding an ORDER BY: 

SELECT cluster_id, P_reads, S_reads, E_reads, total_reads 
FROM all_info 
WHERE (S_reads+1) / (P_reads+1) >= 500 and total_reads>2000 
ORDER BY E_reads DESC  

 
Selecting sub‐sets of genes 

Question: How can I retrieve the list of strobila‐specific genes? For example, we want to see 
only the clusters where more than 80% of all reads originate from the strobila stage. 
Moreover, we want to retrieve only the clusters where the total number of reads is >=10. 

Answer: We have to add additional conditions to our previous query. 80% threshold means 
that dividing the number of reads in strobila (S_reads) by the total number of reads 
(total_reads) we should get values >= 0.8. Total read number (total_reads value) should be 
>=10. So, we need to use a clause WHERE with two conditions: 

SELECT cluster_id, P_reads, S_reads, E_reads, total_reads  
FROM all_info 
WHERE S_reads / total_reads >= 0.8 and total_reads >=10  

As a result we will get a list of 345 strobila‐specific clusters (see Fuchs et al. Fig. 2B): 

And now let us get the list of the polyp‐specific genes. We need to change just one 
parameter (P_reads) and our new query will be: 

SELECT cluster_id, P_reads, S_reads, E_reads, total_reads  
FROM all_info 
WHERE P_reads / total_reads >= 0.8 and total_reads >=10  

As a result we will get a list of 336 polyp‐specific clusters (see Fuchs et al. Fig2B): 
  

Getting the list of the ephyra‐specific genes is easy now. Here is the corresponding query: 

SELECT cluster_id, P_reads, S_reads, E_reads, total_reads  
FROM all_info 
WHERE E_reads / total_reads >= 0.8 and total_reads >=10  

 
Working with the microarray data 

The advantage of a relational database (like MySQL) is that it allows to work with large data 
sets and gives absolute flexibility in "asking" question of any level of complexity. It is easy to 
link different data types together, for example, sequence data with the corresponding 
expression values, peptide prediction, images and so on. For analysing the data one needs to 
describe the "question" as a set of mathematical and logical operators (in a similar way like 
in R, MatLab and similar programs). 

In the following examples we will use the table "array_normalized".  

This table contains mean signal values from independent experiments 
((replicate_1+replicate_2+replicate_3)/3). Mean signal values across the stages (polyp, 14 
days 5‐Aza‐Cytidine, 14 days control, ..., Ephyra) have been normalized based on the 
expression level of elongation factor‐1‐alpha (EF1α). Values in the table has not been 
subjected to logarithmic transformation. Log2 or Log10 transformation is important for 
presenting data in a form of a heat map, but for comparing expression that operation is not 
necessary.  

The table "array_normalized" contains 12 field (columns): 

id_entry   unique entry identifier (primary key) 
id_oligo  unique oligonucleotide name 
P_signal  expression in POLYP (24h at 10°C) 
AZA14_signal  expression in POLYP (14 days at 10°C, incubated in 5‐Aza‐cytidine) 
CON14_signal  expression in POLYP (14 days at 10°C, DMSO control) 
AZA16_signal  expression in POLYP (16 days at 10°C, incubated in 5‐Aza‐cytidine) 
CON16_signal  expression in POLYP (16 days at 10°C, DMSO control) 
ES_signal  expression in STROBILA with 1 segment 
LS_signal  expression in STROBILA with 5 segments 
E_signal  expression in EPHYRA (freshly detached) 
cl_name  long name of the cluster (1‐RA_1 ‐ cluster 1, rc_8‐RA_8 ‐ cluster 8) 
cl_id  unique cluster identifier (ranges from 1 .. 29608) 

1) To view all the values from the table type: 

SELECT * FROM array_normalized 
  

The table with results should appear (if not, please check that the query has been correctly 
written): 

IMPORTANT: The results are limited to 1000 rows by default. Depending of the query, the 
results can be very large and memory consuming. Type "0" in the field for using no limit. 

 
2) To find all the genes where expression in early strobila is 100 times stronger than in a 
polyp type: 
SELECT * FROM array_normalized 
WHERE ES_signal / P_signal >= 100 
ORDER by cl_id ASC 

IMPORTANT: Results will be sorted according to the cluster identifiers in ascending order 
(ORDER by cl_id ASC). 

You can also sort the results by clicking on the column's headers. 

3) To identify genes which are up‐regulated during the temperature induction and might 
function as a strobilation inducer one will need a bit more complex query with many 
conditions (now we will describe the hypothetical model in Fig.4A): 

SELECT * FROM array_normalized 
where 
(P_signal+AZA14_signal+CON14_signal+AZA16_signal+CON16_signal+ES_signal+LS_signal+E
_signal>=100) 
and P_signal<50 
and CON14_signal>AZA14_signal 
and CON16_signal>AZA16_signal 
and ES_signal / P_signal>=5 
and LS_signal / P_signal>=10 
and LS_signal >= 1000 
order by cl_id ASC 

 
Here is the short explanation of the query: 

1) We want to select genes which are expressed not extremely weak ‐ cumulative expression 
must be >= 100: 

(P_signal+AZA14_signal+CON14_signal+AZA16_signal+CON16_signal+ES_signal+LS_signal+E
_signal>=100) 

2) and the expression in the polyp stage must be weak: 

and P_signal<50 

3) now we check that the genes are 5‐AZA‐cytidin sensitive and the expression increases at 
cold temperature: 

and CON14_signal>AZA14_signal 
and CON16_signal>AZA16_signal 
and ES_signal / P_signal>=5 
and LS_signal / P_signal>=10 

4) expression in late strobila must be relatively high: 

and LS_signal >= 1000 

5) ordering according to the cluster IDs (gene idenifiers): 

order by cl_id ASC 

As a result we get a list of potential strobilation inducers (27 clusters) represented by the 
heat map in Fig.4B in Fuchs et al. Expression dynamics of these genes follows the model 
represented in the Fig.4A. 

(See screenshot in the next page) 
  

Thus, by using a simple set of commands one can extract a lot of information with nearly 
unlimited flexibility. It is also possible to combine data from several tables (JOIN statement).  

There is more information about databases and SQL at, for example: 
http://www.w3schools.com/sql/ 

One can also consult MySQL Reference Manual at: 
http://dev.mysql.com/doc/refman/5.5/en/index.html