Drosophila selenophosphate synthetase 1 regulates vitamin B6 metabolism: prediction and confirmation
- Kwang Hee Lee†1, 2,
- Myoung Sup Shim†1,
- Jin Young Kim1,
- Hee Kyoung Jung1,
- Eunji Lee2, 3,
- Bradley A Carlson4,
- Xue-Ming Xu4,
- Jin Mo Park5,
- Dolph L Hatfield4,
- Taesung Park2, 3 and
- Byeong Jae Lee1, 2Email author
© Lee et al; licensee BioMed Central Ltd. 2011
Received: 21 March 2011
Accepted: 24 August 2011
Published: 24 August 2011
There are two selenophosphate synthetases (SPSs) in higher eukaryotes, SPS1 and SPS2. Of these two isotypes, only SPS2 catalyzes selenophosphate synthesis. Although SPS1 does not contain selenophosphate synthesis activity, it was found to be essential for cell growth and embryogenesis in Drosophila. The function of SPS1, however, has not been elucidated.
Differentially expressed genes in Drosophila SL2 cells were identified using two-way analysis of variance methods and clustered according to their temporal expression pattern. Gene ontology analysis was performed against differentially expressed genes and gene ontology terms related to vitamin B6 biosynthesis were found to be significantly affected at the early stage at which megamitochondria were not formed (day 3) after SPS1 knockdown. Interestingly, genes related to defense and amino acid metabolism were affected at a later stage (day 5) following knockdown. Levels of pyridoxal phosphate, an active form of vitamin B6, were decreased by SPS1 knockdown. Treatment of SL2 cells with an inhibitor of pyridoxal phosphate synthesis resulted in both a similar pattern of expression as that found by SPS1 knockdown and the formation of megamitochondria, the major phenotypic change observed by SPS1 knockdown.
These results indicate that SPS1 regulates vitamin B6 synthesis, which in turn impacts various cellular systems such as amino acid metabolism, defense and other important metabolic activities.
Selenium has been reported to provide many health benefits in animals, including humans, when obtained from the diet in adequate amounts. For example, selenium has been known to play roles in cancer prevention, aging retardation, immune augmentation, prevention of heart diseases, muscle development and development [[1–4] and references therein]. Many of the health benefits of selenium are mediated by selenoproteins, which contain selenocysteine (Sec) as a selenium containing amino acid .
Selenophosphate synthetase (SPS) synthesizes selenophosphate (SeP), the active selenium donor in Sec biosynthesis, using selenide and ATP as substrates . SeP serves as a selenium donor during Sec biosynthesis . Sec is contained in all selenoproteins . SPS was first isolated from Escherichia coli as one of the enzymes involved in selenoprotein synthesis and was designated SelD . Only one type of SPS, SelD, exists in lower eukaryotes and eubacteria, however, there are two isoforms of SPS, SPS1 and SPS2, that occur in higher eukaryotes . One of the major differences in the sequences between SPS1 and SPS2 is that SPS1 has an arginine at the position corresponding to Sec in SPS2 .
Although it is not clear why there are two SPSs in higher eukaryotes, recent studies have shown that SPS2 synthesizes SeP from selenide and ATP in vitro, while SPS1 does not have this activity . Loss of function in NIH3T3 cells using RNA interference technology showed that SPS2 is required for selenoprotein biosynthesis, while SPS1 does not affect the biosynthesis of this protein class . While some insects such as the red beetle and silkworm have lost the selenoprotein synthesizing machinery including SPS2, SPS1 is still encoded in the genome of these insects, suggesting SPS1 is required for a function other than SeP synthesis .
Although SPS1 does not catalyze SeP biosynthesis, it plays essential roles in the cell. When the gene encoding SPS1 (SPS1, also designated patufet) was deleted in Drosophila, the embryo showed lethality during development , and reactive oxygen species (ROS) levels increased . The haploinsufficiency of genes involved in the Ras-regulated signaling pathway was also suppressed by SPS1 knockout in Drosophila. From the finding that the SelD (E. coli SPS) mutant of E. coli can be complemented by human SPS1 only when L-Sec is supplemented in the medium, it was suggested that SPS1 is involved in the recycling of Sec . However, the means by which SPS1 may be involved in Sec recycling has not been determined. Recently, it was found that the targeted depletion of SPS1 by RNA interference in Drosophila SL2 cells causes growth inhibition, ROS induction and megamitochondrial formation by increasing intracellular glutamine levels . Interestingly, human SPS1 was found to interact with the soluble liver antigen, which was recently identified as eukaryotic Sec synthase (SecS), and the binding reaction was enhanced by Sec tRNA methylase designated SECp43 [18, 19]. It should be noted that SecS is a pyridoxal phosphate (PLP)-dependent enzyme and, therefore, the uptake and/or activation of vitamin B6 may be related to selenium metabolism [20, 21].
Vitamin B6 is a water-soluble compound that contains a pyridine ring. Vitamin B6 is present in nature as several different forms such as pyridoxal (PL), pyridoxine (PN), pyridoxamine (PM) and their 5'-phosphorylated forms . Before use, these vitamers are converted to PLP, which is the metabolically active form. PLP is used as a cofactor for PLP-dependent enzymes, where the pyridine ring acts as an electron sink during enzymatic reactions. Since animals, including humans, cannot synthesize vitamin B6, they must obtain it from their diet . PLP can be synthesized through several different pathways, and two types of enzymes, kinases and oxidases, participate in these pathways. For PM to be converted to PLP, it is first phosphorylated by a kinase (PL/PM/PN kinase) to form pyridoxamine phosphate (PMP), and then the PMP is oxidized to form PLP using an oxidase (PMP/PNP oxidase). PN can also be converted to PLP using the same kinase and oxidase used for PM. In this case, the phosphorylated intermediate is pyridoxine phosphate (PNP). However, PL can be directly converted to PLP by phosphorylation using a kinase . Therefore, kinases and oxidases are important enzymes for PLP synthesis.
There are more than 100 PLP-dependent enzymes in a cell that perform essential roles in various metabolic pathways including amino acid metabolism (such as amino acid synthesis and degradation), fatty acid metabolism (such as the synthesis of polyunsaturated fatty acids), and carbohydrate metabolism (such as the breakdown of glycogen) [ and references therein]. The PLP-dependent enzymes that participate in amino acid metabolism can be classified into 4 categories: transaminase, racemase, decarboxylase and α,β-eliminase . Interestingly, the biosynthesis of Sec can be mediated by cystathionine β-synthase (CBS) using serine as a precursor and it can also be synthesized by cystathionine γ-lyase (CGL) from selenocystathionine [26, 27]. Both CBS and CGL are PLP-dependent enzymes . In addition, enzymes that are involved in the degradation of Sec, such as selenocysteine lyase (SCL), D-selenocystine α, and β-lyase, use PLP as a cofactor . Recently, it was found that SCL can interact with SPS1 . Therefore, it seems that vitamin B6 participates in the metabolism of Sec, i.e., in the biosynthesis and/or decomposition of Sec.
In the present study, we found that the knockdown of SPS1 led to the down regulation of genes involved in PLP biosynthesis, which, in turn, induced the formation of megamitochondria and the expression of genes responsible for innate immunity. Our findings suggest that SPS1 primarily regulates PLP biosynthesis, and the intracellular PLP level affects various biological processes such as amino acid metabolism, megamitochondrial formation and innate immune response.
Identification and temporal clustering of differentially expressed genes
To analyze the expression pattern of DEGs generated by SPS1 knockdown, clustering of DEGs was performed according to their temporal expression using by self-organizing map (SOM) algorithms . As a result, the DEGs were classified into 6 clusters (Figure 1B). Genes belonging to cluster 1 (33 genes) showed continuous increase in their expression by SPS1 knockdown, and most of them showed more than 4-fold increase on day 5. The expression patterns of genes in cluster 2 (77 genes) were similar to those of cluster 1, but the average expression level was lower than that of cluster 1. Genes in cluster 3 (9 genes) showed down-up patterns of expression. The expression of cluster 4 genes (12 genes) was decreased until day 3, and the expression level was maintained afterward. The expression pattern of genes in cluster 5 (27 genes) was a down-down type. Genes in cluster 6 (80 genes) showed an expression pattern similar to that of cluster 5 genes. However, the average level of expression of cluster 5 genes was much lower than that of cluster 6 genes.
Using six clusters resulted from above, the expression ratios of DEGs composing a cluster were drawn as a box plot according to their sampling date (days 1, 3, and 5). As shown in Figure 1C, the median values (Q2s) of all clusters were close to zero on day 1. However, Q2s of clusters 3, 4 and 5 on day 3 were significantly decreased. On day 5, Q2s of clusters 1 and 2 were significantly increased, while those of clusters 4, 5 and 6 decreased. The interquartile ranges (IQRs) of each cluster were compared to select cluster(s) whose IQRs were significantly deviated. Clusters 3, 4 and 5 revealed significant down regulation compared to the other clusters on day 3. The IQRs of those clusters on day 3 were lower than -0.75. Therefore, the threshold to select clusters whose expression was significantly changed at a specific sampling date was set to the absolute value of 0.75 (see the dashed lines in Figure 1C). A gene pool composing the selected clusters that showed the same expression pattern at the same sampling date was used as a gene-set for gene ontology analysis. As shown in Figure 1C, there is no cluster showing that their IQRs were located at the outside of the threshold range (-0.75~ +0.75) on day 1; thus, no gene was selected for GO analysis from day 1 samples. However, on day 3, the IQRs of clusters 3, 4 and 5 were lower than the lower threshold (-0.75), and the genes in these clusters were defined as the early/down gene-set because their expressions were decreased. Clusters 1 and 2 showed a significant increase in their expression on day 5, and the genes in those clusters were defined as the late/up gene-set. On the other hand, genes in clusters 4, 5 and 6 showed significant down-regulation in their expression, and they were defined as the late/down gene-set (the dotted boxes in Figure 1C; Additional File 2 for the list of genes in these gene-sets).
Identification of statistically overrepresented biological processes by gene ontology analysis
List of biological processes selected from gene ontology analysis
Represented biological process
Max. corrected p-value
Vitamin B6 biosynthesis
AttB, AttD, CecB, DptB, Dro, Drs, Mtk, egr, pirk, PGRP-LF, PGRP-SD, W, Cyp6a8, Cyp12a4, Toll-7
Carboxylic acid metabolism
(Amino acid metabolism)
arg, CG8745, Gs1, Oat, Pepck, yellow-f
Validation of expression by quantitative PCR
Intracellular pyridoxal phosphate level was decreased by SPS1 knockdown
Inhibition of PLP biosynthesis and SPS1 knockdown showed similar expression patterns
Because intracellular PLP levels were significantly reduced after SPS1 knockdown, it can be assumed that PLP biosynthesis is the primary target of SPS1, and the inhibition of PLP synthesis by treating cells with inhibitors will cause similar gene expression patterns as those resulting from SPS1 knockdown. To test this hypothesis, Drosophila cells were treated with 4-DPN for 5 days, and the expression level of genes selected by GO analysis was measured with RT-PCR. As shown in Figure 3B, the level of expression of the early/down genes (CG31472 and CG11899) was not changed by 4-DPN treatment. Because 4-DPN inhibits only the function of proteins that participate in PLP synthesis and does not affect the expression of genes encoding those proteins, it is reasonable that 4-DPN does not affect the expression of CG31472 and CG11899. However, the treatment of 4-DPN affected the expression of genes comprising the late/up and late/down gene-sets. Of the 17 genes tested, 14 genes (82%) showed expression patterns similar to those observed by microarray analysis. It should be noted that the late gene-sets include genes responsible for defense response and amino acid metabolism. These results strongly suggest that PLP synthesis is the primary target of SPS1 and that intracellular PLP levels regulate other important biological processes such as defense system and amino acid metabolism.
The reduction of intracellular PLP level inhibits cell growth and induces megamitochondrial formation
We assumed that the genes whose expression was changed at the early stage after knockdown are involved in the primary target process regulated by SPS1. To identify the primary target, DEGs were isolated after microarray analysis and classified according to their temporal expression pattern; GO terms of early changed DEGs were analyzed using BinGO software. It is interesting that only PLP biosynthesis was predicted from the early/down gene set, even though the parameters were changed. As shown in Table 1, the DEGs in the early/down gene set that are involved in vitamin B6 synthesis are CG31472 and CG11899. CG31472 is an ortholog of mammalian pyridoxine phosphate oxidase (PNPO), which catalyzes PLP production from PMP and PNP and PL production from PN or PM by oxidizing the substrates . The function of CG11899 was not determined experimentally. However, it has high homology with mammalian phosphoserine aminotransferase and PdxC of E. coli, which are responsible for producing 4-phospho-hydroxy threonine, a precursor of the pyridoxine ring . Therefore, it seems that CG11899 plays a role in producing precursors of vitamin B6. Interestingly, intracellular PLP levels were decreased even though only two genes among four genes that are involved in the PLP biosynthesis pathway in Drosophila cell were down-regulated (see Additional File 5). This result suggests that these two genes are involved in an essential step of PLP biosynthesis, or SPS1 may also regulate the other proteins involved in PLP biosynthesis post-transcriptionally.
Because PLP is used as a cofactor for various enzymes that are important for many metabolic pathways, including amino acid metabolism, the inhibition of PLP biosynthesis will lead to the inhibition of cell growth. The inhibition of cell growth induced by SPS1 knockdown seems to be mediated by a decrease in intracellular PLP levels. Specific inhibition of PLP synthesis by 4-DPN treatment led to growth inhibition (Figure 4A), suggesting the growth inhibition by SPS1 knockdown is caused by down-regulation of PLP synthesis.
As described in the Results, down-regulation of genes responsible for PLP synthesis stimulated the expression of DEGs that participate in the defense response. In addition, most of the late gene-sets showed the same pattern of expression as that seen when cells were treated with 4-DPN (Figure 3B). The relationship between vitamin B6 and cellular defense, however, has not been demonstrated before this study. Previously, it was reported that the knockdown of SPS1 induced diphthericin expression in Drosophila SL2 cell when a genome-wide knockdown was performed . The inhibition of PLP synthesis also induced the expression of various AMPs, including dipththericin. Therefore, SPS1 plays a key role in innate immune responses, including AMP production, by regulating PLP level in the cell. The mechanism by which vitamin B6 regulates the innate immune system remains to be elucidated.
The fact that the treatment of 4-DPN, like SPS1 knockdown, induced megamitochondrial formation indicates that intracellular glutamine levels increased with the inhibition of PLP synthesis. Because PLP is used as a cofactor for enzymes that have transaminase activity, it is reasonable to assume that low levels of PLP will lead to the inhibition of synthesis of amino acids such as glutamate or glutamine. However, the inhibition of PLP biosynthesis induced the expression of Gs1 and l(2)01810 (Figure 4C). These two genes are involved in the increase of intracellular glutamine levels . These results suggest that the lack of PLP in the cell provides a signal for compensatory induction of some genes responsible for amino acid metabolism. PLP regulation of the expression of Gs1 and l(2)01810 has not been elucidated.
Although SPS1 was found to regulate the biosynthesis of vitamin B6, the mechanism or signal pathway to which SPS1 is related has not been determined. Because SPS1 is localized to both plasma and nuclear membranes , it can be speculated that SPS1 regulates signal transduction by transducing signals on the plasma membrane or by transporting messengers or transcription factors through the nuclear membrane. The treatment of cell with 4-DPN or SPS1 knockdown induced the expression of PGRP-SD and Toll-7, which are involved in the Toll signaling pathway, and PGRP-LF, which is an activator of the IMD pathway (Figure 5). In addition, Tamo, which is a negative regulator for nuclear import of Dorsal, was found to be one of the down-regulated DEGs. These results strongly suggest that PLP, which is regulated by SPS1, participates in both the Toll and the IMD pathways.
Interestingly, SPS1 knockdown induced down-regulation of CG1753, which encodes cystathionine β-synthase (see Table 1). Cystathionine β-synthase catalyzes both L-cystathionine and L-selenocysteine synthesis . Therefore, it seems that SPS1 regulates the synthesis of Sec indirectly by regulating the expression of Sec synthesizing enzymes.
In this study, we predicted that vitamin B6 biosynthesis is the primary target of SPS1 by employing bioinformatics methods such as microarray and GO analyses and confirmed the prediction experimentally by showing that PLP levels were decreased by SPS1 knockdown and that the inhibition of PLP biosynthesis caused the same phenotypes as SPS1 knockdown.
Materials were purchased from the following sources: Drosophila Schneider cell line 2 (SL2) was purchased from Invitrogen, HyQ SFX-Insect medium from Hyclone, T3 Megascript kit from Ambion, RNeasy mini kit from Qiagen, GeneChip Drosophila genome 2.0 array from Affymetrix, SYBR Green mix from Applied Biosystems, TRIzol reagent from Invitrogen, Moloney murine leukemia virus reverse transcriptase from Super-Bio, 4-deoxypyridoxine hydrochloride from TCI, 5',6,6'-tetrachloro-1,1',3,3'-tetraethylbenzimidazolyl-carbocyanine iodide (JC-1) from Molecular Probes, and oligonucleotides from Cosmo Genetech. The sequences of oligos used for RT-PCR are listed in Additional File 6.
SL2 cell culture and RNA interference
SL2 cell culture and preparation of double-stranded RNAs were carried out as described . Briefly, for RNA interference, 0.25 × 106 cells were plated on a 24-well plate containing 0.5 ml of HyQ SFX-Insect medium. Four micrograms of dsRNAs were added directly to the medium and incubated for 48 hr and cells were split into appropriate culture dishes for further incubation and other experiments.
Microarray experiments were performed using the GeneChip Drosophila genome 2.0 array. After the addition of double stranded RNAs targeting SPS1 to the culture medium, total RNA was extracted from SL2 cells treated with or without SPS1 dsRNA on day 1, 3 and 5 after treatment using the RNeasy mini kit according to the manufacturer's instructions. The cells that were not treated with any dsRNA were used as controls. The RNA quality was checked using Experion (Applied Biosystems) according to the manufacturer's instructions. Five micrograms of total RNAs were reverse transcribed with oligo-dT primer containing a T7 RNA polymerase promoter (TAATACGACTCACTATAGGG). Biotin-labeled cRNAs were generated from the cDNA sample by in vitro transcription with T7 RNA polymerase. The labeled cRNAs were fragmented to an average size of 35-200 bases by mild alkaline treatment at 94°C for 40 min. Fragmented cRNAs were hybridized with probes that are on GeneChip Drosophila genome 2.0 array, and the chips were washed and stained in the Affymetrix Fluidics Station 450 by following the procedures established by Affymetrix (Affymetrix GeneChip R Expression Analysis Technical Manual). The signals were scanned using the GeneChip Scanner 3000 7G (Affymetrix).
Analysis of microarray data
The raw data were imported into Acuity 4.0 software (Molecular Devices, Inc.), and a background adjustment and normalization were performed using robust multichip average (RMA) and quantile methods, respectively, implemented in Acuity 4.0 software [39, 40]. To identify differentially expressed genes (DEGs), a two-way analysis of variance (ANOVA) model was used and fitted using the R software http://www.r-project.org, as described by Park et al. Two models were considered to identify DEGs. Model 1 contains group and time effects as well as their interactions. Model 1 allows the expression level of genes to change over time (days 1, 3 and 5) and these change patterns to differ between groups (control and knockdown). Model 2 includes only group and time effects assuming that the expression level of genes changes over time but these change patterns are the same between groups. From Model 1, DEGs were identified by the genes with significant interaction effects, while from Model 2 DEGs were identified by the genes with significant group effects. The p-values were adjusted by Westfall and Young's method . The genes with adjusted p-values less than 0.1 were identified.
To classify DEGs according to their temporal expression pattern, DEGs were clustered using a self-organizing map (SOM) algorithm implemented in Acuity 4.0 . The ratios of normalized log2 values of DEGs between SPS1 knockdown cells and control cells were used as input data and the SOM map size was set to 3 × 2. The ranges of expression ratios of DEGs within each cluster at each sampling date were displayed by box plot using R software. The interquartile ranges (IQRs) of each cluster were compared to select cluster(s) whose IQRs were significantly deviated. The criterion for determining clusters within which gene expressions were changed significantly was set to 0.75, i.e., when the interquartile range (IQR) of a cluster was larger than +0.75 or smaller than -0.75, the cluster was selected as significantly changed. This is because 0.75 is the threshold value to isolate clusters on day 3 (see Results for more details). The genes composing a cluster selected at the early stage (day 3) were defined as an early responding gene-set and those composing a cluster selected at the late stage (day 5) were defined as a late responding gene-set.
GO analysis was performed by BiNGO version 2.3 , which is plugged in Cytoscape . Gene symbols of each gene-set were used as input data. The parameters were set as follows: assessment was set to overrepresentation, statistical test to binomial test, multiple testing correction to FWER correction, significance level to 0.05. Among GO evidence codes, inferred from electronically annotated (IEA) were discarded. The most significant pathway was predicted by considering the selected GO terms and visualized output.
RT-PCR and quantitative real time RT-PCR
RT- PCR and real time PCR were carried out as described . Briefly, total RNA was isolated from the cells using the TRIzol reagent. cDNAs were synthesized from total RNAs with Moloney murine leukemia virus reverse transcriptase and oligo (dT) primers according to the manufacturer's protocols. RT-PCR was performed with 0.1 μg of template total RNA and specific primers (Additional File 6). RT-PCR products were electrophoresed on a 2% agarose gel and visualized by ethidium bromide. For the measurement of relative mRNA levels of each gene, real time PCR was carried out using an ABI 7300 real time PCR system (Applied Biosystems) as follows. cDNAs were amplified using SYBR Green mix and specific primers for 40 cycles [initial incubation at 50°C for 2 min and then at 95°C for 10 min, and 40 cycles (95°C for 15 sec, 55°C for 1 min and 72°C for 1 min)]. Output data were obtained as Ct values using Sequence Detection Software (SDS) version 1.3 (7300 System, Applied Biosystems) and the differential mRNA expression of each gene between control and knockdown cell was calculated using the comparative Ct method . RP49 mRNA, an internal control, was amplified along with the target genes, and the Ct value of RP49 used to normalize the expression of target genes.
Measurement of intracellular PLP concentration
Cellular PLP levels were determined using the method previously described  with minor modifications. At day 5 after treatment with dsRNA or 4-DPN, cells were washed with phosphate buffered saline and harvested. Cells (6 × 107) were lysed by resuspension in 600 μl of distilled water. Cell extracts were induced to produce the semicarbazon derivative of PLP as follows: 40 μl of 250 mg/ml of both semicarbazide and glycine were added into 500 μl of cell extracts or PLP standard. The mixture was vortexed and incubated at room temperature in the dark for 30 min. Proteins were then precipitated by adding 50 μl of 60% HClO4 into the mixture, and the solution was thoroughly mixed for 1 min. The solution was clarified by centrifugation for 10 min at 15,000 × g, and 30-50 μl of a 25% NaOH solution was added to the supernatant to achieve a pH between 3.0 and 5.0. HPLC was performed using a ZORBAX SB-C18 column (4.6 mm × 25 cm, PN 880975902) and an isocratic mobile phase consisting of 60 mM sodium phosphate (pH 6.5), 400 mg/l EDTA and 9.5% methanol at a flow-rate of 1 ml/min, and the derivatized PLP was quantified using a Waters™ 474 scanning fluorescence detector by setting excitation and emission wavelengths to 380 and 450 nm, respectively.
Mitochondrial staining and confocal microscopy
Mitochondrial staining and confocal microscopy were carried out as described . Briefly, SL2 cells (0.5 × 106) were plated onto a chambered coverglass one day before staining. Cells were incubated with 1 μg/ml JC-1 for 30 min at 25°C, washed three times with HyQ-SFX-Insect medium and observed with a LSM510 confocal microscope (Carl Zeiss) at 512 × 512 pixel resolution through an X63 C-Apochromat objective. Excitation wavelengths for JC-1 aggregate and JC-1 monomer were 543 and 488 nm, respectively.
This work was supported by the Priority Research Centers Program and Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (Grant Nos 2009-0094020 and 2011-0012947 to BJL) and in part by the Intramural Research Program of the National Institutes of Health, National Cancer Institute, Center for Cancer Research. MSS, JYK, KHL and HKJ were supported by Brain Korea 21 Research Fellowship from the Korea Ministry of Education and Human Resources Development.
- Boosalis MG: The Role of Selenium in Chronic Disease. Nutr Clin Pract. 2008, 23 (2): 152-160. 10.1177/0884533608314532.PubMedView Article
- Flohé L: Selenium in mammalian spermiogenesis. J Biol Chem. 2007, 388 (10): 987-995. 10.1515/BC.2007.112.
- Hatfield DL, Berry MJ, Gladyshev VN: Selenium: Its molecular biology and role in human health. 2006, New York: Springer-Verlag Inc., 2View Article
- Tamura T, Yamamoto S, Takahata M, Sakaguchi H, Tanaka H, Stadtman TC, Inagaki K: Selenophosphate synthetase genes from lung adenocarcinoma cells: Sps1 for recycling L-selenocysteine and Sps2 for selenite assimilation. Proc Natl Acad Sci USA. 2004, 101 (46): 16162-16167. 10.1073/pnas.0406313101.PubMed CentralPubMedView Article
- Ehrenreich A, Forchhammer K, Tormay P, Veprek B, Böck A: Selenoprotein synthesis in E. coli. Purification and characterisation of the enzyme catalysing selenium activation. Eur J Biochem. 1992, 206 (3): 767-773. 10.1111/j.1432-1033.1992.tb16983.x.PubMedView Article
- Glass RS, Singh WP, Jung W, Veres Z, Scholz TD, Stadtman TC: Monoselenophosphate: synthesis, characterization, and identity with the prokaryotic biological selenium donor, compound SePX. Biochemistry. 1993, 32 (47): 12555-12559. 10.1021/bi00210a001.PubMedView Article
- Lu J, Holmgren A: Selenoproteins. J Biol Chem. 2009, 284 (2): 723-727.PubMedView Article
- Leinfelder W, Forchhammer K, Veprek B, Zehelein E, Böck A: In vitro synthesis of selenocysteinyl-tRNAUCA from seryl-tRNAUCA: involvement and characterization of the selD gene product. Proc Natl Acad Sci USA. 1990, 87 (2): 543-547. 10.1073/pnas.87.2.543.PubMed CentralPubMedView Article
- Guimaraes MJ, Peterson D, Vicari A, Cock BG, Copeland NG, Gilbert DJ, Jenkins NA, Ferrick DA, Kastelein RA, Bazan JF, Zlotnik A: Identification of a novel selD homolog from eukaryotes, bacteria, and archaea: Is there an autoregulatory mechanism in selenocysteine metabolism?. Proc Natl Acad Sci USA. 1996, 93 (26): 15086-15091. 10.1073/pnas.93.26.15086.PubMed CentralPubMedView Article
- Low SC, Harney JW, Berry MJ: Cloning and functional characterization of human selenophosphate synthetase, an essential component of selenoprotein synthesis. J Biol Chem. 1995, 270 (37): 21659-21664. 10.1074/jbc.270.37.21659.PubMedView Article
- Xu XM, Carlson BA, Mix H, Zhang Y, Saira K, Glass RS, Berry MJ, Gladyshev VN, Hatfield DL: Biosynthesis of selenocysteine on its tRNA in eukaryotes. PLoS Biol. 2007, 5 (1): e4-10.1371/journal.pbio.0050004.PubMed CentralPubMedView Article
- Xu XM, Carlson BA, Irons R, Mix H, Zhong N, Gladyshev VN, Hatfield DL: Selenophosphate synthetase 2 is essential for selenoprotein biosynthesis. Biochem J. 2007, 404 (1): 115-120. 10.1042/BJ20070165.PubMed CentralPubMedView Article
- Lobanov AV, Hatfield DL, Gladyshev VN: Selenoproteinless animals: selenophosphate synthetase SPS1 functions in a pathway unrelated to selenocysteine biosynthesis. Protein Sci. 2008, 17 (1): 176-182.PubMed CentralPubMedView Article
- Alsina B, Serras F, Baguna J, Corominas M: Patufet, the gene encoding the Drosophila melanogaster homologue of selenophosphate synthetase, is involved in imaginal disc morphogenesis. Mol Gen Genet. 1998, 257 (2): 113-123. 10.1007/s004380050630.PubMedView Article
- Morey M, Corominas M, Serras F: DIAP1 suppresses ROS-induced apoptosis caused by impairment of the selD/sps1 homolog in Drosophila. J Cell Sci. 2003, 116 (22): 4597-4604. 10.1242/jcs.00783.PubMedView Article
- Morey M, Serras F, Baguñà J, Hafen E, Corominas M: Modulation of the Ras/MAPK signalling pathway by the redox function of selenoproteins in Drosophila melanogaster. Dev Biol. 2001, 238 (1): 145-156. 10.1006/dbio.2001.0389.PubMedView Article
- Shim MS, Kim JY, Jung HK, Lee KH, Xu XM, Carlson BA, Kim KW, Kim IY, Hatfield DL, Lee BJ: Elevation of glutamine level by selenophosphate synthetase 1 knockdown induces megamitochondrial formation in Drosophila cells. J Biol Chem. 2009, 284 (47): 32881-32894. 10.1074/jbc.M109.026492.PubMed CentralPubMedView Article
- Small-Howard A, Morozova N, Stoytcheva Z, Forry EP, Mansell JB, Harney JW, Carlson BA, Xu XM, Hatfield DL, Berry MJ: Supramolecular complexes mediate selenocysteine incorporation in vivo. Mol Cell Biol. 2006, 26 (6): 2337-2346. 10.1128/MCB.26.6.2337-2346.2006.PubMed CentralPubMedView Article
- Xu XM, Carlson BA, Zhang Y, Mix H, Kryukov GV, Glass RS, Berry MJ, Gladyshev VN, Hatfield DL: New developments in selenium biochemistry: selenocysteine biosynthesis in eukaryotes and archaea. Biol Trace Elem Res. 2007, 119 (3): 234-241. 10.1007/s12011-007-8003-9.PubMedView Article
- Forchhammer K, Böck A: Selenocysteine synthase from Escherichia coli. Analysis of the reaction sequence. J Biol Chem. 1991, 266 (10): 6324-6328.PubMed
- Ganichkin OM, Xu XM, Carlson BA, Mix H, Hatfield DL, Gladyshev VN, Wahl MC: Structure and catalytic mechanism of eukaryotic selenocysteine synthase. J Biol Chem. 2008, 283 (9): 5849-5865.PubMedView Article
- Fitzpatrick TB, Amrhein N, Kappes B, Macheroux P, Tews I, Raschle T: Two independent routes of de novo vitamin B6 biosynthesis: not that different after all. Biochem J. 2007, 407 (1): 1-13. 10.1042/BJ20070765.PubMedView Article
- Mooney S, Leuendorf JE, Hendrickson C, Hellmann H: Vitamin B6: A long known compound of surprising complexity. Molecules. 2009, 14 (1): 329-351. 10.3390/molecules14010329.PubMedView Article
- Gonza'lez E, Danehower D, Daub ME: Vitamer levels, stress response, enzyme activity, and gene regulation of Arabidopsis lines mutant in the pyridoxine/pyridoxamine 5'-phosphate oxidase (PDX3) and the pyridoxal kinase (SOS4) genes involved in the vitamin B6 salvage pathway1[W][OA]. Plant Physiol. 2007, 145 (3): 985-996. 10.1104/pp.107.105189.View Article
- Yoshimura T, Goto M: D-amino acids in the brain: structure and function of pyridoxal phosphate-dependent racemases. FEBS J. 2008, 275 (14): 3527-3537. 10.1111/j.1742-4658.2008.06516.x.PubMedView Article
- Esaki N, Nakamura T, Tanaka H, Suzuki T, Morino Y, Soda K: Enzymatic synthesis of selenocysteine in rat liver. Biochemistry. 1981, 20 (15): 4492-4496. 10.1021/bi00518a039.PubMedView Article
- Meier M, Janosik M, Kery V, Kraus JP, Burkhard P: Structure of human cystathionine β-synthase: a unique pyridoxal 5'-phophate-dependent heme protein. EMBO J. 2001, 20 (15): 3910-3916. 10.1093/emboj/20.15.3910.PubMed CentralPubMedView Article
- Aitken SM, Kirsch JF: The enzymology of cystathionine biosynthesis: strategies for the control of substrate and reaction specificity. Arch Biochem Biophys. 2005, 433 (1): 166-175. 10.1016/j.abb.2004.08.024.PubMedView Article
- Soda K, Oikawa T, Esaki N: Vitamin B6 enzymes participating in selenium amino acid metabolism. BioFactors. 1999, 10 (2): 257-262. 10.1002/biof.5520100225.PubMedView Article
- Tobe R, Mihara H, Kurihara T, Esaki N: Identification of proteins interacting with selenocysteine lyase. Biosci Biotechnol Biochem. 2009, 73 (5): 1230-1232. 10.1271/bbb.90065.PubMedView Article
- Tamayo P, Slonim D, Mesirov J, Zhu Q, Kitareewan S, Dmitrovsky E, Lander ES, Golub TR: Interpreting patterns of gene expression with self-organizing maps: Methods and application to hematopoietic differentiation. Proc Natl Acad Sci USA. 1999, 96 (6): 2907-2912. 10.1073/pnas.96.6.2907.PubMed CentralPubMedView Article
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardsom E, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology. Nat Genet. 2000, 25 (1): 25-29. 10.1038/75556.PubMed CentralPubMedView Article
- Maere S, Heymans K, Kuiper M: BiNGO: a cytoscape plugin to assess overrepresentation of gene ontology categories in biological networks. Bioinformatics. 2005, 21 (16): 3448-3449. 10.1093/bioinformatics/bti551.PubMedView Article
- Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Roy Stat Soc B. 1995, 57: 289-300.
- Musayev FN, Di Salvo ML, Ko TP, Schirch V, Safo MK: Structure and properties of recombinant human pyridoxine 5'-phosphate oxidase. Protein Sci. 2003, 12 (7): 1455-1463. 10.1110/ps.0356203.PubMed CentralPubMedView Article
- Drewke C, Klein M, Clade D, Arenz A, Müller R, Leistner E: 4-O-Phosphoryl-L-threonine, a substrate of the pdxC(serC) gene product involved in vitamin B6 biosynthesis. FEBS Lett. 1996, 390 (2): 179-182. 10.1016/0014-5793(96)00652-7.PubMedView Article
- Foley E, O'Farrell PH: Functional dissection of an innate immune response by a genome-wide RNAi screen. PLoS Biol. 2004, 2 (8): e203-10.1371/journal.pbio.0020203.PubMed CentralPubMedView Article
- Kim JY, Lee KH, Shim MS, Shin H, Xu XM, Carlson BA, Hatfield DL, Lee BJ: Human selenophosphate synthetase 1 has five splice variants with unique interactions, subcellular localizations and expression patterns. Biochem Biophys Res Commun. 2010, 397 (1): 53-58. 10.1016/j.bbrc.2010.05.055.PubMed CentralPubMedView Article
- Bolstad BM, Irizarry RA, Astrand M, Speed TP: A Comparison of normalization methods for high density oligonucleotide array data based on bias and variance. Bioinformatics. 2003, 19 (2): 185-193. 10.1093/bioinformatics/19.2.185.PubMedView Article
- Irizarry RA, Bolstad BM, Collin F, Cope LM, Hobbs B, Speed TP: Summaries of Affymetrix GeneChip probe level data. Nucleic Acids Res. 2003, 31 (4): e15-10.1093/nar/gng015.PubMed CentralPubMedView Article
- Park T, Yi SG, Lee S, Lee S, Yoo DH, Ahn JI, Lee YS: Statistical tests for identifying differentially expressed genes in time course microarray experiments. Bioinformatics. 2003, 19 (6): 694-703. 10.1093/bioinformatics/btg068.PubMedView Article
- Westfall PH, Young SS: Resampling-based multiple testing: examples and methods for p-value adjustment. 1993, New York, Wiley
- Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13 (11): 2498-2504. 10.1101/gr.1239303.PubMed CentralPubMedView Article
- Schmittgen TD, Livak KJ: Analyzing real-time PCR data by the comparative CT method. Nat Protoc. 2008, 3 (6): 1101-1108. 10.1038/nprot.2008.73.PubMedView Article
- Perry C, Yu S, Chen J, Matharu KS, Stover PJ: Effect of vitamin B6 availability on serine hydroxymethyltransferase in MCF-7 cells. Arch Biochem Biophys. 2007, 462 (1): 21-27. 10.1016/j.abb.2007.04.005.PubMed CentralPubMedView Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.