comparison edger.xml @ 14:c5fa04118f83 draft

planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/edger commit 4b17306763bc8eb92c8967c7b4b57abcc514e670
author iuc
date Wed, 04 Sep 2024 15:49:53 +0000
parents 838b481dc6f9
children 119f41ec0ada
comparison
equal deleted inserted replaced
13:838b481dc6f9 14:c5fa04118f83
2 <description> 2 <description>
3 Perform differential expression of count data 3 Perform differential expression of count data
4 </description> 4 </description>
5 <macros> 5 <macros>
6 <token name="@TOOL_VERSION@">3.36.0</token> 6 <token name="@TOOL_VERSION@">3.36.0</token>
7 <token name="@VERSION_SUFFIX@">2</token> 7 <token name="@VERSION_SUFFIX@">5</token>
8 </macros> 8 </macros>
9 <edam_topics> 9 <edam_topics>
10 <edam_topic>topic_3308</edam_topic> 10 <edam_topic>topic_3308</edam_topic>
11 </edam_topics> 11 </edam_topics>
12 <edam_operations> 12 <edam_operations>
66 66
67 #if $anno.annoOpt=='yes': 67 #if $anno.annoOpt=='yes':
68 -a '$anno.geneanno' 68 -a '$anno.geneanno'
69 #end if 69 #end if
70 70
71 -C '${ ','.join( ['%s' % $x.contrast for x in $rep_contrast] ) }' 71 #if $formula:
72 -F '$formula'
73 #end if
74
75 #if $contrasts.contrastOpt == 'file':
76 -C '$contrasts.cinfo'
77 #else:
78 -C '${ ','.join( ['%s' % $x.contrast for x in $contrasts.rep_contrast] ) }'
79 #end if
72 80
73 #if $f.filt.filt_select == 'yes': 81 #if $f.filt.filt_select == 'yes':
74 #if $f.filt.cformat.format_select == 'cpm': 82 #if $f.filt.cformat.format_select == 'cpm':
75 -c '$f.filt.cformat.cpmReq' 83 -c '$f.filt.cformat.cpmReq'
76 -s '$f.filt.cformat.cpmSampleReq' 84 -s '$f.filt.cformat.cpmSampleReq'
120 <option value="matrix">Single Count Matrix</option> 128 <option value="matrix">Single Count Matrix</option>
121 </param> 129 </param>
122 <when value="files"> 130 <when value="files">
123 <repeat name="rep_factor" title="Factor" min="1"> 131 <repeat name="rep_factor" title="Factor" min="1">
124 <param name="factorName" type="text" label="Name" help="Name of experiment factor of interest (e.g. Genotype). One factor must be entered and there must be two or more groups per factor. Optional additional factors (e.g. Batch) can be entered using the Insert Factor button below, see Help section for more information. NOTE: Please only use letters, numbers or underscores, and the first character of each factor must be a letter"> 132 <param name="factorName" type="text" label="Name" help="Name of experiment factor of interest (e.g. Genotype). One factor must be entered and there must be two or more groups per factor. Optional additional factors (e.g. Batch) can be entered using the Insert Factor button below, see Help section for more information. NOTE: Please only use letters, numbers or underscores, and the first character of each factor must be a letter">
125 <sanitizer> 133 <validator type="empty_field"/>
126 <valid initial="string.letters,string.digits"> 134 <validator type="regex" message="Please only use letters, numbers or underscores">^[\w]+$</validator>
127 <add value="_"/>
128 </valid>
129 </sanitizer>
130 </param> 135 </param>
131 <repeat name="rep_group" title="Group" min="2" default="2"> 136 <repeat name="rep_group" title="Group" min="2" default="2">
132 <param name="groupName" type="text" label="Name" help="Name of group that the counts files belong to (e.g. WT or Mut). NOTE: Please only use letters, numbers or underscores (case sensitive), and the first character of each group must be a letter"> 137 <param name="groupName" type="text" label="Name" help="Name of group that the counts files belong to (e.g. WT or Mut). NOTE: Please only use letters, numbers or underscores (case sensitive), and the first character of each group must be a letter">
133 <sanitizer> 138 <sanitizer>
134 <valid initial="string.letters,string.digits"> 139 <valid initial="string.letters,string.digits">
174 <when value="yes"> 179 <when value="yes">
175 <param name="geneanno" type="data" format="tabular" label="Gene Annotations"/> 180 <param name="geneanno" type="data" format="tabular" label="Gene Annotations"/>
176 </when> 181 </when>
177 <when value="no"/> 182 <when value="no"/>
178 </conditional> 183 </conditional>
184 <!-- Optional formula -->
185 <param name="formula" type="text" optional="true" label="Formula for linear model" help="An optional formula for the EdgeR linear model, this will override the use of the fields in factors as a simple sum. The formula can only use elements available in the factors file. This needs to be exactly as EdgeR expect the formula, ie. `~ 0 + factor_A + factor_B:factor_C`. See EdgeR documentation for more details.">
186 <sanitizer invalid_char="">
187 <valid initial="string.letters,string.digits">
188 <add value="_"/>
189 <add value="-"/>
190 <add value="+"/>
191 <add value="*"/>
192 <add value="/"/>
193 <add value="^"/>
194 <add value=":"/>
195 <add value="."/>
196 <add value="~"/>
197 <add value=" "/>
198 <add value="("/>
199 <add value=")"/>
200 <add value="@"/>
201 <add value="$"/>
202 </valid>
203 </sanitizer>
204 </param>
179 <!-- Contrasts --> 205 <!-- Contrasts -->
180 <repeat name="rep_contrast" title="Contrast" min="1" default="1"> 206 <conditional name="contrasts">
181 <param name="contrast" type="text" label="Contrast of Interest" help="Names of two groups to compare separated by a hyphen e.g. Mut-WT. If the order is Mut-WT the fold changes in the results will be up/down in Mut relative to WT. If you have more than one contrast enter each separately using the Insert Contrast button below. For differences between contrasts use e.g. (MT.t1-MT.t0)-(WT.t1-WT.t0). For more info, see Chapter 8 in the limma User's guide: https://www.bioconductor.org/packages/release/bioc/vignettes/limma/inst/doc/usersguide.pdf or https://bioconductor.org/packages/release/bioc/vignettes/edgeR/inst/doc/edgeRUsersGuide.pdf page 36 for nested comparisons."> 207 <param name="contrastOpt" type="select" label="Input contrasts manually or through a file">
182 <validator type="empty_field"/> 208 <option value="manual">manually</option>
183 <validator type="regex" message="Please only use letters, numbers, parentheses or underscores">^[\w\-()]+$</validator> 209 <option value="file">file</option>
184 </param> 210 </param>
185 </repeat> 211 <when value="manual">
212 <repeat name="rep_contrast" title="Contrast" min="1" default="1">
213 <param name="contrast" type="text" label="Contrast of Interest" help="Names of two groups to compare separated by a hyphen e.g. Mut-WT. If the order is Mut-WT the fold changes in the results will be up/down in Mut relative to WT. If you have more than one contrast enter each separately using the Insert Contrast button below. For differences between contrasts use e.g. (MT.t1-MT.t0)-(WT.t1-WT.t0). For more info, see Chapter 8 in the limma User's guide: https://www.bioconductor.org/packages/release/bioc/vignettes/limma/inst/doc/usersguide.pdf or https://bioconductor.org/packages/release/bioc/vignettes/edgeR/inst/doc/edgeRUsersGuide.pdf page 36 for nested comparisons.">
214 <validator type="empty_field"/>
215 <validator type="regex" message="Please only use letters, numbers, parentheses or underscores">^[\w\-()]+$</validator>
216 </param>
217 </repeat>
218 </when>
219 <when value="file">
220 <param name="cinfo" optional="true" type="data" format="tabular" label="Contrasts File" help="Setting this file will ignore any manually added contrasts above, make sure to remove any contrast fields above pressing the trash bin icon, or the tool will fail. First line of the file must be a header, below that each separate contrast should be on a line. Contrast formulas need to be based on ther factors data and potentially the formula provided. See EdgeR documentation on contrasts for more details."/>
221 </when>
222 </conditional>
186 <!-- Filter Options --> 223 <!-- Filter Options -->
187 <section name="f" expanded="false" title="Filter Low Counts"> 224 <section name="f" expanded="false" title="Filter Low Counts">
188 <conditional name="filt"> 225 <conditional name="filt">
189 <param name="filt_select" type="select" label="Filter lowly expressed genes?" help="Treat genes with very low expression as unexpressed and filter out. See the Filter Low Counts section below for more information. Default: No"> 226 <param name="filt_select" type="select" label="Filter lowly expressed genes?" help="Treat genes with very low expression as unexpressed and filter out. See the Filter Low Counts section below for more information. Default: No">
190 <option value="no" selected="true">No</option> 227 <option value="no" selected="true">No</option>
262 <param name="counts" value="matrix.txt"/> 299 <param name="counts" value="matrix.txt"/>
263 <repeat name="rep_factor"> 300 <repeat name="rep_factor">
264 <param name="factorName" value="Genotype"/> 301 <param name="factorName" value="Genotype"/>
265 <param name="groupNames" value="Mut,Mut,Mut,WT,WT,WT"/> 302 <param name="groupNames" value="Mut,Mut,Mut,WT,WT,WT"/>
266 </repeat> 303 </repeat>
304 <param name="contrastOpt" value="manual"/>
267 <repeat name="rep_contrast"> 305 <repeat name="rep_contrast">
268 <param name="contrast" value="Mut-WT"/> 306 <param name="contrast" value="Mut-WT"/>
269 </repeat> 307 </repeat>
270 <repeat name="rep_contrast"> 308 <repeat name="rep_contrast">
271 <param name="contrast" value="WT-Mut"/> 309 <param name="contrast" value="WT-Mut"/>
300 <param name="counts" value="matrix-complex.txt"/> 338 <param name="counts" value="matrix-complex.txt"/>
301 <repeat name="rep_factor"> 339 <repeat name="rep_factor">
302 <param name="factorName" value="Genotype"/> 340 <param name="factorName" value="Genotype"/>
303 <param name="groupNames" value="MutA,MutA,MutA,MutB,MutB,MutB,WTA,WTA,WTA,WTB,WTB,WTB"/> 341 <param name="groupNames" value="MutA,MutA,MutA,MutB,MutB,MutB,WTA,WTA,WTA,WTB,WTB,WTB"/>
304 </repeat> 342 </repeat>
343 <param name="contrastOpt" value="manual"/>
305 <repeat name="rep_contrast"> 344 <repeat name="rep_contrast">
306 <param name="contrast" value="(MutA-MutB)-(WTA-WTB)"/> 345 <param name="contrast" value="(MutA-MutB)-(WTA-WTB)"/>
307 </repeat> 346 </repeat>
308 <param name="normalisationOption" value="TMM"/> 347 <param name="normalisationOption" value="TMM"/>
309 <output_collection name="outTables" count="1"> 348 <output_collection name="outTables" count="1">
331 <param name="counts" value="matrix.txt"/> 370 <param name="counts" value="matrix.txt"/>
332 <repeat name="rep_factor"> 371 <repeat name="rep_factor">
333 <param name="factorName" value="Genotype"/> 372 <param name="factorName" value="Genotype"/>
334 <param name="groupNames" value="Mut,Mut,Mut,WT,WT,WT"/> 373 <param name="groupNames" value="Mut,Mut,Mut,WT,WT,WT"/>
335 </repeat> 374 </repeat>
375 <param name="contrastOpt" value="manual"/>
336 <repeat name="rep_contrast"> 376 <repeat name="rep_contrast">
337 <param name="contrast" value="Mut-WT"/> 377 <param name="contrast" value="Mut-WT"/>
338 </repeat> 378 </repeat>
339 <param name="normalisationOption" value="TMM"/> 379 <param name="normalisationOption" value="TMM"/>
340 <output_collection name="outTables" count="1"> 380 <output_collection name="outTables" count="1">
354 <param name="counts" value="matrix.txt"/> 394 <param name="counts" value="matrix.txt"/>
355 <repeat name="rep_factor"> 395 <repeat name="rep_factor">
356 <param name="factorName" value="Genotype"/> 396 <param name="factorName" value="Genotype"/>
357 <param name="groupNames" value="Mut,Mut,Mut,WT,WT,WT"/> 397 <param name="groupNames" value="Mut,Mut,Mut,WT,WT,WT"/>
358 </repeat> 398 </repeat>
399 <param name="contrastOpt" value="manual"/>
359 <repeat name="rep_contrast"> 400 <repeat name="rep_contrast">
360 <param name="contrast" value="Mut-WT"/> 401 <param name="contrast" value="Mut-WT"/>
361 </repeat> 402 </repeat>
362 <param name="normalisationOption" value="TMM"/> 403 <param name="normalisationOption" value="TMM"/>
363 <output name="outReport"> 404 <output name="outReport">
381 </repeat> 422 </repeat>
382 <repeat name="rep_factor"> 423 <repeat name="rep_factor">
383 <param name="factorName" value="Batch"/> 424 <param name="factorName" value="Batch"/>
384 <param name="groupNames" value="b1,b2,b3,b1,b2,b3"/> 425 <param name="groupNames" value="b1,b2,b3,b1,b2,b3"/>
385 </repeat> 426 </repeat>
427 <param name="contrastOpt" value="manual"/>
386 <repeat name="rep_contrast"> 428 <repeat name="rep_contrast">
387 <param name="contrast" value="Mut-WT"/> 429 <param name="contrast" value="Mut-WT"/>
388 </repeat> 430 </repeat>
389 <param name="normalisationOption" value="TMM"/> 431 <param name="normalisationOption" value="TMM"/>
390 <output_collection name="outTables" count="1"> 432 <output_collection name="outTables" count="1">
400 <test expect_num_outputs="2"> 442 <test expect_num_outputs="2">
401 <param name="format" value="matrix"/> 443 <param name="format" value="matrix"/>
402 <param name="ffile" value="yes"/> 444 <param name="ffile" value="yes"/>
403 <param name="finfo" value="factorinfo.txt"/> 445 <param name="finfo" value="factorinfo.txt"/>
404 <param name="counts" value="matrix.txt"/> 446 <param name="counts" value="matrix.txt"/>
447 <param name="contrastOpt" value="manual"/>
405 <repeat name="rep_contrast"> 448 <repeat name="rep_contrast">
406 <param name="contrast" value="Mut-WT"/> 449 <param name="contrast" value="Mut-WT"/>
407 </repeat> 450 </repeat>
408 <param name="normalisationOption" value="TMM"/> 451 <param name="normalisationOption" value="TMM"/>
409 <output_collection name="outTables" count="1"> 452 <output_collection name="outTables" count="1">
422 <param name="counts" value="matrix.txt"/> 465 <param name="counts" value="matrix.txt"/>
423 <repeat name="rep_factor"> 466 <repeat name="rep_factor">
424 <param name="factorName" value="Genotype"/> 467 <param name="factorName" value="Genotype"/>
425 <param name="groupNames" value="Mut,Mut,Mut,WT,WT,WT"/> 468 <param name="groupNames" value="Mut,Mut,Mut,WT,WT,WT"/>
426 </repeat> 469 </repeat>
470 <param name="contrastOpt" value="manual"/>
427 <repeat name="rep_contrast"> 471 <repeat name="rep_contrast">
428 <param name="contrast" value="Mut-WT"/> 472 <param name="contrast" value="Mut-WT"/>
429 </repeat> 473 </repeat>
430 <param name="normalisationOption" value="TMM"/> 474 <param name="normalisationOption" value="TMM"/>
431 <output_collection name="outTables" count="2"> 475 <output_collection name="outTables" count="2">
449 <param name="counts" value="matrix.txt"/> 493 <param name="counts" value="matrix.txt"/>
450 <repeat name="rep_factor"> 494 <repeat name="rep_factor">
451 <param name="factorName" value="Genotype"/> 495 <param name="factorName" value="Genotype"/>
452 <param name="groupNames" value="Mut,Mut,Mut,WT,WT,WT"/> 496 <param name="groupNames" value="Mut,Mut,Mut,WT,WT,WT"/>
453 </repeat> 497 </repeat>
498 <param name="contrastOpt" value="manual"/>
454 <repeat name="rep_contrast"> 499 <repeat name="rep_contrast">
455 <param name="contrast" value="Mut-WT"/> 500 <param name="contrast" value="Mut-WT"/>
456 </repeat> 501 </repeat>
457 <param name="normalisationOption" value="TMM"/> 502 <param name="normalisationOption" value="TMM"/>
458 <param name="lrtOption" value="true"/> 503 <param name="lrtOption" value="true"/>
492 <param name="countsFile" value="WT3.counts,Mut3.counts"/> 537 <param name="countsFile" value="WT3.counts,Mut3.counts"/>
493 </repeat> 538 </repeat>
494 </repeat> 539 </repeat>
495 <param name="annoOpt" value="yes"/> 540 <param name="annoOpt" value="yes"/>
496 <param name="geneanno" value="anno.txt"/> 541 <param name="geneanno" value="anno.txt"/>
542 <param name="contrastOpt" value="manual"/>
497 <repeat name="rep_contrast"> 543 <repeat name="rep_contrast">
498 <param name="contrast" value="Mut-WT"/> 544 <param name="contrast" value="Mut-WT"/>
499 </repeat> 545 </repeat>
500 <repeat name="rep_contrast"> 546 <repeat name="rep_contrast">
501 <param name="contrast" value="WT-Mut"/> 547 <param name="contrast" value="WT-Mut"/>
528 <param name="counts" value="matrix.txt"/> 574 <param name="counts" value="matrix.txt"/>
529 <repeat name="rep_factor"> 575 <repeat name="rep_factor">
530 <param name="factorName" value="Genotype"/> 576 <param name="factorName" value="Genotype"/>
531 <param name="groupNames" value="Mut,Mut,Mut,WT,WT,WT"/> 577 <param name="groupNames" value="Mut,Mut,Mut,WT,WT,WT"/>
532 </repeat> 578 </repeat>
579 <param name="contrastOpt" value="manual"/>
533 <repeat name="rep_contrast"> 580 <repeat name="rep_contrast">
534 <param name="contrast" value="Mut-WT"/> 581 <param name="contrast" value="Mut-WT"/>
535 </repeat> 582 </repeat>
536 <param name="normalisationOption" value="TMM"/> 583 <param name="normalisationOption" value="TMM"/>
537 <param name="filt_select" value="yes"/> 584 <param name="filt_select" value="yes"/>
563 <param name="counts" value="matrix.txt"/> 610 <param name="counts" value="matrix.txt"/>
564 <repeat name="rep_factor"> 611 <repeat name="rep_factor">
565 <param name="factorName" value="Genotype"/> 612 <param name="factorName" value="Genotype"/>
566 <param name="groupNames" value="Mut,Mut,Mut,WT,WT,WT"/> 613 <param name="groupNames" value="Mut,Mut,Mut,WT,WT,WT"/>
567 </repeat> 614 </repeat>
615 <param name="contrastOpt" value="manual"/>
568 <repeat name="rep_contrast"> 616 <repeat name="rep_contrast">
569 <param name="contrast" value="Mut-WT"/> 617 <param name="contrast" value="Mut-WT"/>
570 </repeat> 618 </repeat>
571 <param name="normalisationOption" value="TMM"/> 619 <param name="normalisationOption" value="TMM"/>
572 <param name="filt_select" value="yes"/> 620 <param name="filt_select" value="yes"/>
597 <param name="counts" value="matrix.txt"/> 645 <param name="counts" value="matrix.txt"/>
598 <repeat name="rep_factor"> 646 <repeat name="rep_factor">
599 <param name="factorName" value="Genotype"/> 647 <param name="factorName" value="Genotype"/>
600 <param name="groupNames" value="Mut,Mut,Mut,WT,WT,WT"/> 648 <param name="groupNames" value="Mut,Mut,Mut,WT,WT,WT"/>
601 </repeat> 649 </repeat>
650 <param name="contrastOpt" value="manual"/>
602 <repeat name="rep_contrast"> 651 <repeat name="rep_contrast">
603 <param name="contrast" value="Mut-WT"/> 652 <param name="contrast" value="Mut-WT"/>
604 </repeat> 653 </repeat>
605 <param name="normalisationOption" value="TMM"/> 654 <param name="normalisationOption" value="TMM"/>
606 <param name="filt_select" value="yes"/> 655 <param name="filt_select" value="yes"/>
624 <not_has_text text="-0.0682"/> 673 <not_has_text text="-0.0682"/>
625 </assert_contents> 674 </assert_contents>
626 </element> 675 </element>
627 </output_collection> 676 </output_collection>
628 </test> 677 </test>
678 <!-- Ensure formula and contrast file work -->
679 <test expect_num_outputs="2">
680 <param name="format" value="matrix"/>
681 <param name="counts" value="matrix.txt"/>
682 <repeat name="rep_factor">
683 <param name="factorName" value="Genotype"/>
684 <param name="groupNames" value="Mut,Mut,Mut,WT,WT,WT"/>
685 </repeat>
686 <repeat name="rep_factor">
687 <param name="factorName" value="Batch"/>
688 <param name="groupNames" value="b1,b2,b3,b1,b2,b3"/>
689 </repeat>
690 <param name="contrastOpt" value="file"/>
691 <param name="cinfo" value="contrasts_file.txt"/>
692 <param name="formula" value="~ 0 + Genotype + Batch"/>
693 <param name="normalisationOption" value="TMM"/>
694 <output_collection name="outTables" count="3">
695 <element name="edgeR_Mut-WT" ftype="tabular">
696 <assert_contents>
697 <has_text_matching expression="GeneID.*logFC.*logCPM.*F.*PValue.*FDR"/>
698 <has_text_matching expression="11304.*0.4584"/>
699 </assert_contents>
700 </element>
701 <element name="edgeR_WT-Mut" ftype="tabular">
702 <assert_contents>
703 <has_text_matching expression="GeneID.*logFC.*logCPM.*F.*PValue.*FDR"/>
704 </assert_contents>
705 </element>
706 <element name="edgeR_(2*Mut_3*WT)-WT" ftype="tabular">
707 <assert_contents>
708 <has_text_matching expression="GeneID.*logFC.*logCPM.*F.*PValue.*FDR"/>
709 </assert_contents>
710 </element>
711 </output_collection>
712 </test>
629 </tests> 713 </tests>
630 <help><![CDATA[ 714 <help><![CDATA[
631 .. class:: infomark 715 .. class:: infomark
632 716
633 **What it does** 717 **What it does**
711 795
712 *Factor Name:* The name of the experimental factor being investigated e.g. Genotype, Treatment. One factor must be entered, the name should start with a letter and spaces must not be used. Optionally, additional factors can be included, these are variables that might influence your experiment e.g. Batch, Gender, Subject. If additional factors are entered, an additive linear model will be used. 796 *Factor Name:* The name of the experimental factor being investigated e.g. Genotype, Treatment. One factor must be entered, the name should start with a letter and spaces must not be used. Optionally, additional factors can be included, these are variables that might influence your experiment e.g. Batch, Gender, Subject. If additional factors are entered, an additive linear model will be used.
713 797
714 *Groups:* The names of the groups for the factor. The names should start with a letter, and only contain letters, numbers and underscores, other characters such as spaces and hyphens must not be used. If entered into the tool form above, the order must be the same as the samples (to which the groups correspond) are listed in the columns of the counts matrix, with the values separated by commas. 798 *Groups:* The names of the groups for the factor. The names should start with a letter, and only contain letters, numbers and underscores, other characters such as spaces and hyphens must not be used. If entered into the tool form above, the order must be the same as the samples (to which the groups correspond) are listed in the columns of the counts matrix, with the values separated by commas.
715 799
800 **Formula:**
801 By default the tool will construct a formula for modelling counts based on the contents of the factors files or the factors given.
802 This can be overriden by directly providing the EdgeR formula in section named Formula.
803
716 **Contrasts of Interest:** 804 **Contrasts of Interest:**
717 The contrasts you wish to make between levels. 805 The contrasts you wish to make between levels.
718 A common contrast would be a simple difference between two levels: "Mut-WT" 806 A common contrast would be a simple difference between two levels: "Mut-WT"
719 represents the difference between the mutant and wild type genotypes. 807 represents the difference between the mutant and wild type genotypes.
720 Multiple contrasts must be entered separately using the Insert Contrast button, spaces must not be used. 808 Multiple contrasts must be entered separately using the Insert Contrast button, spaces must not be used.
809
810 Alternatively, you can specify a file with contrasts. The file must contain a header (it's value is irrelevant)
811 and one contrast per line on the first column (other columns are ignored). If using this option, make sure to
812 remove any contrast section from the manual part, or the tool will fail.
721 813
722 **Filter Low Counts:** 814 **Filter Low Counts:**
723 Genes with very low counts across all libraries provide little evidence for differential expression. 815 Genes with very low counts across all libraries provide little evidence for differential expression.
724 In the biological point of view, a gene must be expressed at some minimal level before 816 In the biological point of view, a gene must be expressed at some minimal level before
725 it is likely to be translated into a protein or to be biologically important. In addition, the 817 it is likely to be translated into a protein or to be biologically important. In addition, the