Sawyer and Malik 1010.1073/pnas.0605468103
Fig. 4. Multiple alignments of NEJ1 and XRS2 used for PAML analysis. Gapped regions have been pruned from the alignments.
References
1. Kellis M, Patterson N, Endrizzi M, Birren, B, Lander ES (2003) Nature 423:241-254.
2. Boyle EI, Weng S, Gollub J, Jin H, Botstein D, Cherry JM, Sherlock G (2004) Bioinformatics 20:3710-3715.
3. Bustamante CD, Fledel-Alon A, Williamson S, Nielsen R, Hubisz MT, Glanowski S, Tanenbaum DM, White TJ, Sninsky JJ, Hernandez RD, et al. (2005) Nature 437:1153-1157.
4. Clark AG, Glanowski S, Nielsen R, Thomas PD, Kejariwal A, Todd MA, Tanenbaum DM, Civello D, Lu F, Murphy B, et al. (2003) Science 302:1960-1963.
5. Nielsen R, Bustamante C, Clark AG, Glanowski S, Sackton TB, Hubisz MJ, Fledel-Alon A, Tanenbaum DM, Civello D, White TJ, et al. (2005) PLoS Biol 3:e170.
6. Voight BF, Kudaravalli S, Wen, X, Pritchard JK (2006) PLoS Biol 4:e72.
7. Wang ET, Kodama G, Baldi, P, Moyzis RK (2006) Proc Natl Acad Sci USA 103:135-140.
Table 2. Reexamination of genes with whole-gene values of dN/dS > 1 (based on the analysis of Kellis et al.)
|
ORF |
dN/dSa |
Ra |
Sa |
L bp |
Overlapsb |
Feature typec |
|
7.21 |
26 |
1 |
441 |
YDR432W* |
Dubious |
|
|
5.28 |
33 |
6 |
351 |
YLR281C*, YLR279W |
Dubious |
|
|
YKL111C |
4.31 |
25 |
3 |
336 |
YKL112W |
Dubious |
|
3.82 |
29 |
6 |
327 |
YLR262C* |
Dubious |
|
|
YAR068W |
2.33 |
NAd |
486 |
None |
Uncharacterized |
|
|
1.80 |
22 |
6 |
327 |
YJL008C |
Dubious |
|
|
1.17 |
33 |
12 |
360 |
YGL073W |
Dubious |
|
|
1.10 |
NAd |
NAd |
306 |
None |
Dubious |
|
|
1.05 |
49 |
18 |
345 |
None |
||
|
1.00 |
4 |
0 |
1011 |
None |
Verified |
|
|
1.00 |
0 |
20 |
723 |
None |
Verified |
|
|
1.00 |
0 |
48 |
609 |
None |
Verified |
|
|
1.00 |
0 |
53 |
915 |
None |
Verified |
|
|
1.00 |
0 |
18 |
285 |
None |
Verified |
|
|
1.00 |
0 |
5 |
366 |
None |
Verified |
|
|
1.00 |
0 |
34 |
525 |
None |
Verified |
|
|
1.00 |
0 |
11 |
327 |
None |
Verified |
|
|
1.00 |
0 |
9 |
483 |
None |
Verified |
|
|
1.00 |
2 |
1 |
171 |
YDL229W |
Dubious |
|
|
1.00 |
0 |
0 |
360 |
None |
Verified |
|
|
1.00 |
0 |
3 |
303 |
None |
Verified |
|
|
1.00 |
0 |
3 |
303 |
None |
Verified |
|
|
1.00 |
0 |
21 |
411 |
None |
Verified |
|
|
1.00 |
0 |
13 |
312 |
None |
Verified |
Previous analysis (1) of the S. cerevisiae-S. paradoxus genomes found nine genes with whole-gene values of dN/dS greater than one. We find that eight of these genes are now annotated as "dubious" ORFs on the Saccharomyces Genome Database (SGD), and all but three overlap significantly with another ORF on the same or opposite strand. Our reanalysis suggests that there are no S. cerevisiae-S. paradoxus orthologous gene pairs with a whole-gene value of dN/dS > 1. Furthermore, most of the genes identified as having dN/dS = 1 actually have no replacement changes, with one exception being YNL326C.
* Overlaps a gene on the same strand, but the two ORFs are in different reading frames.
a dN/dS values reported are as calculated by Kellis et al. [Kellis M, Patterson N, Endrizzi M, Birren B, Lander ES (2003) Nature 423:241-254]. "R" and "S" are actual numbers of replacement and synonymous changes calculated in the present study. This indicates a consistent problem in the calculation of dN/dS ratios when there are no R changes.
b This column reports instances where the gene in question overlaps another ORF, either on the same or opposite strand. "C" and "W" at the end of gene names indicate Watson or Crick strand.
c Listed by the Saccharomyces Genome Database (SGD) at Stanford. SGD defines a dubious ORF as follows: ORFs which are not conserved in other Saccharomyces species and for which there is no experimental evidence that a gene product is produced in S. cerevisiae. Many ORFs classified as "Dubious" are small and overlap a larger ORF of the class "Verified" or "Uncharacterized"; however, overlap with another ORF does not mandate that an ORF be classified as "Dubious."
d No unambiguous S. paradoxus orthologs have been identified by the Fungal Alignment tool at SGD.
e YNL326C is also called PFA3. The SGD description for this gene is: "Palmitoyltransferase for Vac8p, required for vacuolar membrane fusion; contains an Asp-His-His-Cys-cysteine rich (DHHC-CRD) domain; autoacylates; required for vacuolar integrity under stress conditions"
f There are multiple S. paradoxus orthologs for this gene identified by the Fungal Alignment tool at SGD.