cDNA microarray analysis of differential gene expression in boar testes during the prepubertal period.
Growth and development of animals depend on a number of factors that regulate gene expression and cell behavior. Transient changes of gene expression in various tissues and organs during the early postnatal period have been amply demonstrated (Chen et al., 2004; Ponsuksili et al., 2007; Wu et al., 2008). Such differential expression of genes at this critical time in early development is required for proper differentiation and proliferation of cells in a certain tissue. Thus, it is important to determine the temporal profile of gene expression during the post-natal development period.
The boar has a number of unique features concerning the synthesis of uncommon steroid hormones that include unusually high serum levels of anabolic steroids such as androstenedione and 19-nortestosterone (Debruychere and Van Peteghem, 1991; Schwarzenberger et al., 1993). The boar is one of only a few males that have a high estrogen concentration in the serum (Schwarzenberger et al., 1993). It is noteworthy that production of these anabolic steroids and estrogens is dramatically increased between post-natal weeks 2~4, followed by a transient decrease to basal levels at 2 months of age until puberty (Schwarzenberger et al., 1993). Significant production of anabolic steroids in the boar results from the action of cytochrome P450 aromatase (Choi et al., 1997; Kao et al., 2000). Our previous study has shown that such transient increase of 19-nortestosterone concentration in the serum during early pre-natal development is closely related with expression of steroidogenic enzymes and estrogen receptors in boar testis (Choi et al., 2007). High concentrations of 19-nortestosterone and 19-norandrostenedione have been detected between 1~2 weeks of age in the testis among various boar tissues examined (Choi et al., 2007). In addition, expression levels of cytochrome P450 aromatase, 17[alpha]-hydroxylase, and 3[beta]-hydroxysteroid dehydrogenase genes are significantly increased at 2 weeks of age, compared with 0, 1, 3, 5, and 6 weeks of age (Choi et al., 2007). Moreover, the latter study has found that expression of estrogen receptor [alpha] and [beta] in the boar testis is greatly increased at 2 weeks of age. Thus, up-regulation of gene expression of steroidogenic enzymes in the boar testis during the neonatal period could be responsible for significant production of anabolic steroids. Indeed, our recent study has demonstrated differential expression of a number of steroidogenic enzymes in boar testis during the neonatal and pre-pubertal periods (Choi et al., 2009). In this study, we further have documented that mRNA and protein levels of enzymes involved in the synthesis of testosterone and estrogen are significantly increased during neonatal development at 2-3 weeks of age compared with levels at pre-puberty (4 months of age). Along with the findings from a previous study (Choi et al., 2007), these observations imply that boar testis around 2 weeks of age during the post-natal period is metabolically active and undergo significant morphological and endocrinological changes.
Complementary DNA (cDNA) microarray analysis is a powerful biological technique commonly used to simultaneously detect differentially expressing genes. A number of previous researches have utilized cDNA microarray analyses to reveal differential expression of genes in various pig tissues (Kim et al., 2005; Guo et al., 2008). In the present study, we first constructed a normalized cDNA library from 2-weeks-old boar testis. Then we attempted to determine the differential gene expression profile in boar testis at post-natal ages of 3 days, 2 weeks, 7 weeks, and 13 weeks. These results were used in an attempt to characterize unknown genes found to be highly expressed at 2 weeks of age.
MATERIALS AND METHODS
Animals and tissue collection
The male reproductive tract from boars (Landrace x Yorksire x Duroc) was obtained from local farms in Kyungbuk province, Republic of Korea, by surgical castration. The testes were separated from the rest of the male reproductive tract in ice-cold phosphate buffered saline (PBS) and were rapidly frozen in liquid nitrogen. The testes were stored at -80[degrees]C until required for further studies. Testes from a total of twenty-one 10~14-days-old piglets were used for construction of a normalized cDNA library. In addition, testes were obtained from boars at 3 days (n = 4), 2 weeks (n = 3), 7 weeks (n = 2), and 13 weeks (n = 6) of age for comparison of gene expression using DNA microarray analysis.
Construction of normalized cDNA library
Total RNA was isolated from testes using TRIzol[TM] (Invitrogen, Carlsbad, CA, USA) and stored in diethylpyrocarbonate (DEPC)-treated distilled water as described previously (Choi et al., 2007). The quantities and qualities of total RNA were determined using a ND-1000 spectrophotometer (Nanodrop Technologies, Wilmington, DE, USA) and agarose gel electrophoresis, respectively. Reverse transcription reactions were carried out to synthesize first-strand cDNA from total RNA in a mixture of 1 [micro]g of total RNAs, oligo (dT)18 primer, 10 mM of dNTP, 4 [micro]l of 5 x buffer, 1 of RNase out (RNase inhibitor, Takara, Shiga, Japan), and 1 [micro]l of Superscript II reverse transcriptase (Invitrogen). Newly synthesized cDNAs were stored at -20[degrees]C until used later. To construct a normalized cDNA library, total RNA collected from the testes at 10~14 days-of-age were pooled and mRNAs were isolated using the Absolutely mRNA[TM] purification kit (Stratagene, La Jolla, CA, USA). Highly purified mRNAs were reverse-transcribed to generate first-stranded DNAs using oligo-dT primer and reverse transcriptase, and second-stranded DNAs were generated using first-stranded DNA templates and DNA polymerase. These purified double-strand (ds)DNAs were then hybridized with the reverse-transcribed cDNAs described above and treated with duplex-specific nuclease at 65[degrees]C for 2 min and at 68[degrees]C for 5 h to remove DNA-DNA or DNA-RNA duplexes. Isolated cDNAs were size-fractionated and collected by a gel filtration column. Quality and the size of the collected cDNAs were verified using a 5% non-denaturing acrylamide gel. To construct the cDNA library, 5'-ends of newly synthesized dscDNAs were ligated with EcoRI adapter DNA by DNA ligase, and the 3'-ends of the cDNA were trimmed by XhoI. These cDNAs were size-fractionated through a gel filtration column and separated by agarose gel electrophoresis. cDNAs of 600 bps or larger in size were collected and isolated by phenol-chloroform extraction and ethanol precipitation. Isolated cDNAs were inserted into a pBK-CMV vector and packed using a bacteriophage packaging extract (Gigapack III Gold Packaging Extract; Stratagene). After titration, a total of 300,000 independent clones were confirmed; cDNA clones were aliquoted in 384 well plates and stored at -80[degrees]C.
DNA sequencing and data analysis
Plasmid DNA from cDNAs clones used to construct the cDNA library were purified by the standard alkaline lysis method and were amplified by the dye-determinator cycling method (BigDye v3.1) using 5' universal primers in the vector (Stratagene). Then, sequences of plasmid DNAs were determined using an ABI3730 XL capillary sequencer (Applied Biosystems Inc, Foster City, CA, USA). The chromatogram data were initially submitted to Phred (Ewing et al., 1998) for base calling and quality assignment. The trace files were trimmed with trim-alt 0.12 and the sequence below 100 bps was removed. In addition, vector trimming was conducted with a cross-match software (http://www.phrap.org), and poly-A was removed using a trimmest of Emboss package (Rice et al., 2000). The prepared multifasta format data were clustered and assembled using the TGICL package (Pertea et al., 2003). Finally, the contigs and singletons were analyzed with NCBI local BLAST (Altschul et al., 1990). Sequences were searched against the NCBI NR databases (2005. 05). Matches with E-value of 1 x [e.sup.-10] for BLASTX were considered as putative hits.
cDNA microarray fabrication
A cDNA microarray was prepared from the cDNA library. Polymerase chain reaction (PCR) was conducted according to a standard protocol. The T7 (5'-TAATACGACTCACTATAGGG-3') and T3 (5'-ATTAACCCTCACTAAAGGGA-3') promoter primers were used for amplification in a 50 [micro]l reaction mix consisting of Premix taqTM (TaKaRa Bio Inc., Shiga, Japan), 0.4 [micro]M of each primer and 1 [micro]l of DNA extract. Thermocycling conditions consisted of an initial denaturation step at 94[degrees]C for 3 min, then 35 cycles of 94[degrees]C for 1 min, annealing at 58[degrees]C for 1.5 min and extension at 72 [degrees]C for 1 min, with final extension at 72[degrees]C for 3 min. PCR products were visualized under ultraviolet (UV) light on 1.0% agarose gel stained with ethidium bromide and reaction mixtures were subjected to 35 cycles of amplification. The size and amount of the PCR products were verified using 1% agarose gel. The PCR products were purified by ethanol precipitation followed by re-suspension in 15 [micro]l of spotting solution (GenoCheck, Korea) and spotting onto a CMT-GAPS II glass slide (Corning; Corning, NY, USA) with a pixsys 5500 arrayer (Cartesian Technologies, Irvine, CA, USA) using Stealth Micro spotting pins. For normalization of signals, a housekeeping gene ([beta]-actin) and a positive control of Arabidopsis thaliana genes were included on each slide. The printed slides were processed according to CMT-GAPS II slide protocol. Briefly, the spots were rehydrated with 1 x SSC for 1 min and then DNAs were linked using a UV cross-linker (Stratagene). The slides were soaked in succinic anhydride/sodium borate solution for 15 min with gentle agitation and then transferred to a 95[degrees]C water bath for 2 min. The slides were quickly transferred to 95% ethanol for 1 min and then dried by centrifugation at 3,000 rpm for 20 s.
DNA chip hybridization and data analysis
Total RNAs from boar testes at 3 days, 2 weeks, 7 weeks, and 13 weeks of age were isolated using TRIzol[TM] (Invitrogen). Fluorescence-labeled cDNA probes were prepared from 30 [micro]g of total RNA by oligo [(dT).sub.18]-primed polymerization using SuperScript II reverse transcriptase (Invitrogen) in a total volume of 30 [micro]l including 25 mM each of dATP, dTTP, and dGTP, 10 mM of dCTP, and 3 mM of Cy3 or Cy5 labeled dCTP (NEN Life Science Product Inc., Boston, MA, USA). After reverse transcription, the labeled cDNAs were concentrated using ethanol precipitation followed by re-suspension of Cy3 and Cy5 labeled cDNAs in 10 [micro]l of hybridization solution (GenoCheck, Seoul, Korea). Two cDNAs labeled with different dyes were mixed in an equal volume and denatured at 95[degrees]C for 2 min, followed by incubation at 45[degrees]C for 20 min. The cDNA mixture was then placed on a spotted slide and covered by a cover slip. The slides were hybridized in a hybridization chamber at 62[degrees]C for 12 h. The slides were sequentially washed in 2 x SSC, 0.1% sodium dodecyl sulfate (SDS) for 2 min, 1 x SSC for 3 min, and 0.2 x SSC for 2 min at room temperature. The slides were centrifuged at 3,000 rpm for 20 s to dry. After the hybridization reaction, slides were scanned with GenePix 4000B scanner (Axon Instruments, Union City, CA, USA). The scanned images were analyzed with GenePix Pro 5.1 software (Axon Instruments) and GeneSpring GX 7.3.1 software (Agilent Technologies, Santa Clara, CA, USA). To allow for an algorithm-mediated elimination of all bad spots, no data points were eliminated by visual inspection from the initial GenePix image. To filter out unreliable data, spots with signal-to-noise ratio below 100 were not included in the data. Data were normalized by global, lowness, print-tip, and scaled normalization for data reliability. Data were clustered in groups of genes that behaved similarly across time course experiments using GeneSpring GX 7.3.1 (Agilent Technologies). We used an algorithm based on the Pearson correlation to separate the genes of similar patterns. The distance cut-off was statistically considered as significant when there were two-fold differences between time course experiments. The correlation cut-off was 0.95, and the accuracy of microarray analyses in the present study was confirmed by real-time PCR analysis.
Leg muscle satellite cells and Leydig cells of the testis were isolated from boars (110 kg body weight) at a local slaughtering house. Skeletal muscle was washed in PBS. Then, each muscle was minced using sterilized scissors and digested by trypsin-EDTA (GIBCO-BRL, Valencia, CA, USA) for 2 h. Digested tissue was centrifuged at 90 g for 3 min and the upper phase was filtered using 40 [micro]m pore size cell strainer. The filtrate was centrifuged at 2,500 rpm at room temperature for 20 min. The collected cell pellet was washed three times using Dulbecco's Modified Eagle's Medium (DMEM; HyClone Laboratories, Logan, UT, USA) containing 1% penicillin/ streptomycin (GIBCO-BRL). The muscle satellite cells were cultured in DMEM supplemented with 10% fetal bovine serum (FBS; GIBCO-BRL), 1% penicillin/ streptomycin, and 0.1% Fungi Zone (GIBCO-BRL) in a 5% C[O.sub.2] incubator at 37[degrees]C. Porcine testes were digested by trypsin-EDTA and the digested cells were cultured in DMEM/Ham's F-12 (1/1, v/v) (GIBCO-BRL) supplemented with 10% FBS. After 3 h, the medium was removed and replaced with fresh medium. Purified Leydig cells were cultured for 4 days.
Real-time PCR analysis and statistical analysis
The first strand cDNA was synthesized using 1 [micro]g of total RNA template isolated from boar testes, as described earlier. Real-time PCR analysis was performed on an ABI7500 apparatus (Applied Biosystems) in Power SYBR[R] green PCR Master Mix[TM] (Applied Biosystems) according to the manufacturer's instructions. Primer sets for real-time PCR were designed by Primer 3 software (http://frodo.wi.mit.edu) from gene sequences (Table 1). The cDNA product was amplified by PCR with Taq DNA polymerase (Bioneer, Daejeon, Korea). A melting curve was analyzed to check the absence of mispriming. The sizes of the PCR products were confirmed by agarose gel electrophoresis. The fold difference of each gene expression was determined by [2.sup.-[DELTA][DELTA]Ct] formula (Lay et al., 2002). All values are presented as mean [+ or -] SEM. Data were analyzed by ANOVA according to the general linear model procedure, followed by Duncan's new multiple range test to determine significant differences at p<0.05. All statistical procedures were performed with the SAS[R] ver. 9.12 software package (SAS, Cary, NC, USA).
RESULTS AND DISCUSSION
A normalized cDNA library of 10~14 days old boar testis was generated in the present study. Titration of the normalized cDNA library resulted in approximately 3 x [10.sup.5] independent clones. Among these, 2,016 clones were randomly selected and the 5' end sequence of each was analyzed using the universal primer in the vector. After eliminating the lower quality DNA sequence (Phred score >13) and vector sequence, the total size of the nucleotide sequences was 803,532 bps, in which 198 clones were smaller than 100 bps in size and 1,818 clones were 100-899 bps in size (Table 2). The sequencing success rate was 90.2% and the average read length was 437 bps. In addition, clustering followed by assembly of the nucleotide sequences of the 1,818 clones yielded 423 contigs and 403 singletons. This data indicated that there were 826 unique DNA sequences. Cluster 1-68 contained 5-22 sequences, and other clusters had only two sequences (Figure 1).
Of these, 203 contigs and 326 singletons showed significant matches with known-function genes, whereas the other 220 contigs and 77 singletons were unknown-function genes. The clones of 423 contigs and 403 singletons were initially utilized to fabricate a cDNA microarray. Sixty-two out of 826 genes were excluded from a final construction of the microarray after checking with PCR and we constructed a DNA microarray for the genes above 100 bp. The DNA microarray comprised a total of 719 genes, which contained PCR-amplified cDNAs of 709 genes, two negative control genes, and eight positive control genes. Using this fabricated cDNA microarray, analysis was carried out to determine differentially expressed genes in boar testes among different age groups (3 days, 2 weeks, 7 weeks, and 13 weeks of age). The design of the cDNA hybridization for DNA chip analysis is shown in Figure 2A, and representative images of the hybridization are shown in Figure 2B. To detect the differential gene expression between different developmental periods, hierarchical clustering of the genes expressed in testes at different stages were shown considering each stage as controlled at a time. Day 3, and week 2, 7 and 13 samples were considered as control at a time to observe up-, and down-regulated genes in the other samples (Figure 2C). Results from the microarray analysis showed that 19 genes were up-regulated and 10 genes were down-regulated in the boar testis at 2 weeks of age, compared to those at 3 days of age (Table 3). Similarly, nine genes were down-regulated and 12 genes were up-regulated at 7 weeks of age. In addition, there were nine up-regulated genes and eight down-regulated genes at 13 weeks of age. Among the 19 genes highly up-regulated at 2 weeks of age, 12 genes had known functions and seven genes had unknown functions. These known functions included bovine osteonectin, pig guanosine diphosphate dissociation inhibitor 2, pig muscle creatine kinase, and pig skeletal alpha actin gene. In addition, further analysis revealed that expression levels of seven unknown-function genes at 2 weeks of age were 2-6-fold higher than those in the other age groups (Table 3 and Figure 3).
[FIGURE 1 OMITTED]
To verify results of the cDNA microarray analysis, real-time PCR analyses were performed with the pooled RNA samples used for microarray analysis, as well as with the individual RNA samples used to make the pooled RNA samples. With pooled RNA samples, three no-hit genes (6B06, 8E11, 5B11) showed high expression levels at 2 weeks of age, compared with those at other ages (Figure 4A). However, with individual RNA samples, six of seven unknown function genes did not show significantly increased expression at 2 weeks of age (Figure 4B). Expression of a no-hit gene, 5B11, was significantly increased at 2 weeks of age. The full-length nucleotide sequencings and the blast analyses of the 8E11 and 5B11 clones showed that 8E11 and 5B11 genes had high homology with human hemoglobin alpha subunit, and Homo sapiens hypothetical LOC643669 transcript variant 1, respectively (Figure 5). Furthermore, we examined tissue expression of the presence of 5B11 and 5G07 genes in pig Leydig cells, muscle cells, and adipocyte cells (Figure 6). Real-time PCR analyses showed that both 5B11 and 5G07 genes were present in pig Leydig cells and adipocyte cells (Figure 6A and B). However, expression of only the 5B11 gene was evident in pig muscle cells (Figure 6A).
A significant production of 19-nortestosterone in the boar during early post-natal development is detected at 2 weeks of age (Schwarzenberger et al., 1993). However, the precise physiological function of 19-nortestosterone in the boar is unclear. In the present study, cDNA microarray analysis identified genes differentially expressed in boar testes during early development. A normalized cDNA library was constructed using the mRNA isolated from the boar testes at 10~14-days-of-age, at which time a significant amount of 19-nortestosterone is secreted from the testis (Choi et al., 2007). An advantage of a normalized cDNA library over other cDNA library methods is the reduced redundancy of the clones in the library (Sasaki et al., 1994). The library clones presently sequenced included a number of genes coding for enzymes involved in the synthesis or metabolism of steroid hormones such as cytochrome P450 17[alpha]-hydroxylase, cytochrome P450 side chain cleavage enzyme, cytochrome b, 17[beta]-estradiol dehydrogenase, 3[beta]-hydroxysteroid dehydrogenase, and catechol O-methyltransferase. Interestingly, the result from clustering of the clones based upon their sequence homology revealed that the largest contig was assembled by 22 clones and contained coding sequence for carbonyl reductase 1 (CBR1). CBR1 is a member of the short-chain dehydrogenase/ reductase superfamily, which is responsible for catalyzed reduction of ketones on androgens and progesterone in porcines (Tanaka et al., 1992; Hoffmann and Master, 2007). Our recent study showed that the expression of CBR1 mRNA and protein in pig testis increases according to age during the neonatal period and is exclusively localized in Leydig cells (Choi et al., 2009). Thus, CBR1 may be involved in steroid hormone synthesis in porcine Leydig cells. In this regard, catechol-O-methyltransferase is an enzyme involved in the breakdown of dopamine and catechol-estrogens (Boulton and Eisenhofer, 1998).
[FIGURE 2 OMITTED]
[FIGURE 3 OMITTED]
In the present study, cDNA microarray analyses showed that a number of genes were expressed at high levels in boar testis at 2 weeks of age (Table 3). Among these genes, the seven no-hit genes were found to be up-regulated at 2 weeks of age. However, only three genes (6B06, 8E11, 5B11) were verified to be true when the pooled RNA samples were used for real time PCR analysis. In contrast, the 6B06 and 5B11 genes could be confirmed when individual RNA samples were used. The reason for the disagreement in gene expression pattern between the microarray analysis and real time PCR data is due to the individual variations of gene expression in the testes caused by the difference of development stage (10~14-days-of-age). Pigs display a dramatic change of testis development during the neonatal period (Franca et al., 2000; Moran et al., 2002). Besides, in this experiment normalized cDNA microarray fabrication was performed by using the PCR product from the constructed cDNA library. For the construction of the cDNA library, total RNA was extracted from the whole tissue of testis. Therefore, there is high probability of presence of RNA from untargeted cells and tissues. However, for the confirmation of cDNA microarray data total RNA was extracted from Leydig cell culture. This may be the reason for being able to confirm only 1 gene expression by real time RTPCR among the 7 up-regulated genes by cDNA microarray analysis. The 5B11 and 8E11 genes, initially considered as no-hit genes due to the lack of match with the gene sequence in the database, turned out to be hemoglobin and hypothetical protein, respectively. Because these genes are expressed not only in Leydig cells, but also in myogenic satellite cells and adipocyte cells, it would be of interest to learn the physiological role of these genes in boar testes during early development.
[FIGURE 4 OMITTED]
In conclusion, we have developed a useful and efficient research tool and materials which can be used to investigate differential gene expression in boar testes during early development. Our interest in the distinct molecular character of the pig testis at 2 weeks was due to detection of high concentrations of 19-nortestosterone and 19-norandrostenedione at this age. With this objective, this research showed up-regulation of a higher number of genes (19 genes) in 2 weeks old pig testes compared to the other ages studied. Among 19 genes, 7 no hit genes were detected by the cDNA library and two of them were confirmed by real time RT PCR. We speculate their involvement in the high concentration of 19-nortestosterone and 19-norandrostenedione in the 2 week-old boar at this point. Further studies of the known and unknown genes identified in the study will lead to a better understanding of the role of 19-nortestosterone in this species.
[FIGURE 5. i OMITTED]
[FIGURE 5. ii OMITTED]
[FIGURE 6 OMITTED]
The authors would like to thank all members of the molecular biology laboratory for their assistance. This project was supported by a grant (20050401-034-712) from the BioGreen 21 Program, Rural Development Administration, Republic of Korea.
Altschul, S. F., W. Gish, W. Miller, E. W. Meyers and D. J. Lipman. 1990. Basic local alignment search tool. J. Mol. Biol. 215:403-410.
Boulton, A. A. and G. Eisenhofer. 1998. Catecholamine metabolism. From molecular understanding to clinical diagnosis and treatment. Adv. Pharmacol. 42:273-292.
Chen, H. W., S. L. Yu, W. J. Chen, P. C. Yang, C. T. Chien, H. Y. Chou, H. N. Li, K. Peck, C. H. Huang, F. Y. Lin, J. J. Chen and Y. T. Lee. 2004. Dynamic changes of gene expression profiles during postnatal development of the heart in mice. Heart 90:927-934.
Choi, I., J.-Y. Kim, E. J. Lee, Y. Y. Kim, C. S. Chung, J. Chang, N.-J. Choi, H.-J. Chung and K.-H. Lee. 2009. Ontogeny of expression and localization of steroidogenic enzymes in the neonatal and prepubertal pig testes. J. Androl. 30:57-74.
Choi, I., D. L. Troyer, D. L. Cornwell, K. R. Kirby-Dobbels, W. R. Collante and F. A. Simmen. 1997. Closely related genes encode developmental and tissue isoforms of porcine cytochrome P450 aromatase. DNA Cell Biol. 16:769-777.
Choi, N. J., J. H. Hyun, J. M. Choi, E. J. Lee, K. H. Cho, Y. Kim, J. Chang, I. B. Chung, C. S. Chung and I. Choi. 2007. Testicular expression steroidogenic enzyme genes is related to a transient increase in serum 19-nortestosterone during neonatal development in pigs. Asian-Aust. J. Anim. Sci. 20:1832-1842.
Debruyckere, G. and C. Van Peteghem. 1991. Detection of 19-nortestosterine and its urinary metabolites in miniature pigs by gas chromatography-mass spectrometry. J. Chromatogr. 564:393-403.
Ewing, B., L. Hillier, M. C. Wendl and P. Green. 1998. Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res. 8:175-185.
Franca, L. R., V. A. Silva, H. Chiarini-Garcia, S. K. Garcia and L. Debeljuk. 2000. Cell proliferation and hormonal changes during postnatal development of the testis in the pig. Biol. Reprod. 63:1629-1636.
Guo, W., S. H. Wang, H. J. Cao, K. Xu, J. Zhang, Z. L. Du, W. Lu, J. D. Feng, N. Li, C. H. Wu and L. Zhang. 2008. Gene microarray analysis for porcine adipose tissue: comparison of gene expression between Chinese Xiang pig and large white. Asian-Aust. J. Anim Sci. 21:11-18.
Hoffmann, F. and E. Maser. 2007. Carbonyl reductase and pluripotent hydroxysteroid dehydrogenase of the short-chain dehydrogenase/reductase superfamily. Drug Metab. Rev. 39:87-144.
Kao, Y. C., T. Higashiyama, X. Sun, T. Okubo, C. Yarborough, I. Choi, Y. Osawa, F. A. Simmen and S. Chen. 2000. Catalytic differences between porcine blastocyst and placental aromatase isozymes. Eur. J. Biochem. 267:6134-6139.
Kim, C. W., K. T. Chang, Y. H. Hong, W. Y. Jung, E. J. Kwon, K. K. Cho, K. H. Chung, B. W. Kim, J. G. Lee, J. S. Yeo, Y. S. Kang and Y. K. Joo. 2005. cDNA microarray analysis of the gene expression profile of swine muscle. Asian-Aust. J. Anim. Sci. 18:1080-1087.
Lay, S. L., C. Lefrere, C. Trautwein, I. Dugail and S. Krief. 2002. Insulin and sterol-regulatory element-binding protein-1c (SREBP-1C) regulation of gene expression in 3T3-L1 adipocytes. J. Biol. Chem. 277:35625-35634.
Pertea, G., X. Huang, F. Liang, V. Antonescu, R. Sultana, S. Karamycheva, Y. Lee, J. White, F. Cheung, B. Parvizi, J. Tsai and J. Quackenbush. 2003. TIGR gene indices clustering tools (TGICL): a software system for fast clustering of large EST datasets. Bioinformatics 19:651-652.
Ponsuksili, S., E. Murani, C. Walz, M. Schwerin and K. Wimmers. 2007. Pre- and postnatal hepatic gene expression profiles of two pig breeds differing in body composition: insight into pathways of metabolic regulation. Physiol. Genomics 29:267-279.
Rice, P., I. Longden and A. Bleasby. 2000. EMBOSS: The European molecular biology open software suite. Trends Genet. 16:276-277.
Sasaki, Y. F., D. Avusawa and M. Oishi. 1994. Construction of a normalized cDNA library by introduction of a semi-solid mRNA--cDNA hybridization system. Nucleic Acids Res. 22:987-992.
Schwarzenberger, F., G. S. Toole, H. L. Christie and J. I. Raeside. 1993. Plasma levels of several androgens and estrogens from birth to puberty in male domestic pigs. Acta Endocrinol. (Copenh) 128:173-177.
Stanton, H. C. and S. K. Woo. 1978. Development of adrenal medullary function in swine. Am. J. Physiol. 234:E137-E145
Tanaka, M., S. Ohno, S. Adachi, S. Nakajin, M. Shinoda and Y. Nagahama. 1992. Pig testicular 20 beta-hydroxysteroid dehydrogenase exhibits carbonyl reductase-like structure and activity. cDNA cloning of pig testicular 20 beta-hydroxysteroid dehydrogenase. J. Biol. Chem. 267:13451-13455.
Wu, F., S. Lee, M. Schumacher, A. Jun and S. Chakravarti. 2008. Differential gene expression patterns of the developing and adult mouse cornea compared to the les and tendon. Exp. Eye Res. 87:214-225.
Dong-Mok Lee (1, a), Ki-Ho Lee (2, a), Jin Ho Choi (1), Jin Hee Hyun (3), Eun Ju Lee (3), Prati Bajracharya (3), Yong Seok Lee (4), Jongsoo Chang (5), Chung Soo Chung (6) and Inho Choi (1, 3) *
(1) School of Biotechnology, Yeungnam University, Gyeongsan, Korea
(2) Department of Biochemistry and Molecular Biology, and Medical Sciences Research Institute, Eulji University, Daejeon, Korea.
(3) Department of Biotechnology, Yeungnam University, Gyeongsan, Korea.
(4) Department of Parasitology, College of Medicine and Frontier Inje Research for Science and Technology, Inje University, Busan, Korea
(5) Department of Agricultural Sciences, Korea National Open University, Seoul, Korea.
(6) Department of Animal Science, Chungbuk National University, Cheongju, Korea.
(a) These two authors contributed equally to this work.
* Corresponding Author: Inho Choi. Tel: +82-53-810-3024, Fax: +82-53-810-4769, E-mail: firstname.lastname@example.org
Received January 2, 2009; Accepted April 2, 2009
Table 1. Primer sequences for real-time RT PCR analysis Molecule Forward primer sequence Reverse primer sequence (5' [right arrow] 3') (5' [right arrow] 3') GAPDH GTCCACTGGTGTCTTCACGA GTCATGAGTCCCTCCACGAT HPRT GGCAAAACAATGCAAACCTT GGGACTCCAGATGTTTCCAA PPIA AGCACTGGGGAGAAAGGATT GCCATCCAACCACTCAGTCT 2-F-06 AACCCAAACACTCCAAAACG CAGTTCGGTCGTCTGGATTT 3-B-04 AAGCAGGTGGACGTGTAACC AAGCAGGTGAGTTGGCTAGG 4-H-05 CGGACAGAGGTTGGAGAGAG GGTGTCAGGTGCTTCAGGAT 5-B-11 CACGCCTCTCTGGACAAGTT GCCCACTCAGGCTTTATTCA 5-G-07 CCAAATGCACAGATGTCCTG CACCTCCGAGATGACCTGTT 6-B-06 CGAGGAGCACTGCACTAGGT CTTCCTAGGACCAGCCCTTC 8-E-11 CAAGCCGTGGAGTAGGTTGT TGCTCGGCCCTAGTTTACTC PPIA = Cyclophilin A, GAPDH = Glyceraldehyde-3-phsphate dehydrogenase, HPRT = Hypoxanthine phosphoribosyltransferase. Table 2. The nucleotide sequence results of normalized cDN library clones Read length Read number Success 800-899 5 (phred score >13 700-799 231 and sequence length 600-699 262 over 100 bp after 500-599 241 vector trimming) 400-499 277 300-399 222 200-299 302 100-199 278 Subtotal 1,818 Fail 198 Total 2,016 Table 3. List of differentially expressed genes in boar testes at 2 weeks, 7 weeks, and 13 weeks compared to day 3 along with their annotation D3 2 Wks 7 Wks 13 Wks Annotation 1.00 2.08 0.81 0.99 Bovine osteonectin 1.00 3.96 1.36 0.66 Canis familiaris type I procollagen pro-alpha 2 chain (COL1A2) 1.00 2.32 0.84 1.02 Homo sapiens full open reading frame cDNA clone RZPDo834D0918D for gene SPARC 1.00 3.67 1.47 1.14 hypothetical protein [Pongo pygmaeus] 1.00 2.81 1.12 0.97 Sus scrofa guanosine diphosphate dissociation inhibitor 2 1.00 5.96 1.33 1.00 Sus scrofa skeletal alpha actin gene 1.00 2.66 0.40 0.34 Unnamed protein product [Mus musculus] 1.00 4.50 1.29 1.03 PREDICTED: Canis familiaris LOC488984 (LOC488984) 1.00 4.61 0.86 0.46 Homo sapiens clone FLC0562 PRO2841 1.00 2.24 0.58 0.57 PREDICTED: Canis familiaris similar to pregnancy-related serine protease (LOC488789) 1.00 2.36 0.68 0.40 Human mRNA for pro-alpha-1 type 3 collagen 1.00 2.60 0.73 0.63 Sus scrofa muscle creatine kinase (CKM) 1.00 2.43 0.81 0.78 No hit 1.00 2.87 0.47 0.39 No hit 1.00 2.37 0.79 0.84 No hit 1.00 6.67 1.02 0.70 No hit 1.00 3.44 0.43 0.37 No hit 1.00 2.92 0.47 0.48 No hit 1.00 3.07 1.34 1.22 No hit 1.00 1.19 11.05 2.06 ovulatory protein-2 precursor [Salvelinus fontinalis] 1.00 1.62 3.52 1.31 Pig MHC class II (SLA-DR-alpha) 1.00 1.32 3.07 1.31 Sus scrofa MHC class I antigen (SLA-1) 1.00 1.59 2.72 1.15 Bos taurus BAC CH240-60O13 1.00 1.43 2.36 1.03 Sus scrofa CD74 antigen (CD74) 1.00 0.86 2.31 0.93 Cow glutathione peroxidase (GPx) plasma isoform 1.00 1.16 2.04 0.89 Sus scrofa MHC class II antigen (SLA-DRB) 1.00 0.84 2.86 0.92 No hit 1.00 1.22 2.52 1.10 No hit 1.00 6.43 15.49 11.24 Sus scrofa Ig gamma 4 constant region (IgG4) 1.00 1.20 1.98 2.99 Gekko japonicus GekBS180P 1.00 1.10 1.86 2.93 Equus caballus procollagen alpha-1 type III precursor (COL3A1) 1.00 1.70 1.44 2.57 Sus scrofa clone Clu_30933.scr.msk.p1.Contig5 1.00 0.89 1.85 2.21 Full-length cDNA clone CS0DI087YM04 of Placenta Cot 25-normalized of Homo sapiens 1.00 1.00 1.52 2.10 Sus scrofa JAK1 mRNA for Janus kinase 1 1.00 1.09 1.28 2.08 Sus scrofa clone Clu_36446 1.00 1.29 1.84 2.30 No hit 1.00 1.23 1.93 2.14 No hit 1.00 0.48 0.86 1.23 Homo sapiens mRNA for KIAA1746 protein 1.00 0.46 0.82 0.87 Bos taurus endoplasmic reticulum thioredoxin superfamily member, 18 kDa (TLP19) 1.00 6.01 2.33 2.70 Bos taurus organic anion transporting polypeptide 2b1 (SLCO2B1) 1.00 0.41 0.80 1.19 Full-length cDNA clone CS0DC001YF17 of Neuroblastoma Cot 25-normalized of Homo sapiens 1.00 0.43 0.79 0.83 PREDICTED: Canis familiaris similar to SSP411 protein (LOC491082) 1.00 0.61 2.03 2.51 PREDICTED: Pan troglodytes similar to alpha 2 macroglobulin (LOC465372) 1.00 0.40 0.88 1.05 S.scrofa (BleyI-L) mRNA for Leydig insulin-like hormone 1.00 0.47 0.93 0.96 No hit 1.00 0.41 0.72 0.73 No hit 1.00 0.40 0.76 0.77 No hit 1.00 4.61 0.86 0.46 Homo sapiens clone FLC0562 PRO2841 1.00 2.02 0.51 0.61 Homo sapiens CDC5 cell division cycle 5-like (S. pombe) 1.00 1.96 0.72 0.49 hypothetical protein KIAA0603 1.00 2.24 0.58 0.57 PREDICTED: Canis familiaris similar to pregnancy-related serine protease (LOC488789) 1.00 2.66 0.40 0.34 Unnamed protein product [Mus musculus] 1.00 2.36 0.68 0.40 Human mRNA for pro-alpha-1 type 3 collagen 1.00 2.60 0.73 0.63 Sus scrofa muscle creatine kinase (CKM) 1.00 3.44 0.43 0.37 No hit 1.00 2.92 0.47 0.48 No hit 1.00 2.87 0.47 0.39 No hit 1.00 1.93 0.56 0.36 No hit 1.00 1.76 0.58 0.34 No hit 1.00 0.67 1.34 0.50 TPA: Sus scrofa mRNA for putative ISG12(a) protein (ISG12(a) gene) 1.00 0.65 0.54 0.48 Sus scrofa gastrin-binding protein 1.00 0.70 0.53 0.47 Sus scrofa membrane bound cytochrome b5 (CYB5) 1.00 0.64 0.55 0.46 Sus scrofa clone Clu_3292 1.00 0.62 0.54 0.46 PREDICTED: Pan troglodytes similar to 40S ribosomal protein S2 (LOC468578) 1.00 0.50 0.78 0.45 Bos taurus 7-dehydrocholesterol reductase (DHCR7) 1.00 0.78 0.57 0.48 No hit 1.00 0.77 0.81 0.48 No hit
|Printer friendly Cite/link Email Feedback|
|Author:||Lee, Dong-Mok; Lee, Ki-Ho; Choi, Jin Ho; Hyun, Jin Hee; Lee, Eun Ju; Bajracharya, Prati; Lee, Yong S|
|Publication:||Asian - Australasian Journal of Animal Sciences|
|Date:||Aug 1, 2009|
|Previous Article:||Genetic relationships among different breeds of Chinese gamecocks revealed by mtDNA variation.|
|Next Article:||Microsatellite-based genetic diversity and evolutionary relationships of six dog breeds.|