Immune molecular profiling of whole blood drawn from a non-human primate cardiac xenograft model treated with anti-CD154 monoclonal antibodies.

Most studies of xenografts have been carried out with complex immunosuppressive regimens to prevent immune rejection; however, such treatments may be fatal owing to unknown causes. Here, we performed immune molecular profiling following anti-CD154 monoclonal antibody (mAb) treatment in heterotopic abdominal cardiac xenografts from α-1,3-galactosyltransferase-knockout pigs into cynomolgus monkeys to elucidate the mechanisms mediating the undesirable fatal side effects of immunosuppressive agents. Blood samples were collected from healthy monkeys as control and then at 2 days after xenograft transplantation and just before humane euthanasia; 94 genes related to the immune system were analyzed. The basic immunosuppressive regimen included cobra venom factor, anti-thymocyte globulin, and rituximab, with and without anti-CD154 mAbs. The maintenance therapy was followed with tacrolimus, MMF, and methylprednisolone. The number of upregulated genes was initially decreased on Day 2 (-/+ anti-CD154 mAb, 22/13) and then increased before euthanasia in recipients treated with anti-CD154 mAbs (-/+ anti-CD154 mAb, 30/37). The number of downregulated genes was not affected by anti-CD154 mAb treatment. Additionally, the number of upregulated genes increased over time for both groups. Interestingly, treatment with anti-CD154 mAbs upregulated coagulation inducers (CCL2/IL6) before euthanasia. In conclusion, immunosuppressive regimens used for cardiac xenografting affected upregulation of 6 inflammation genes (CXCL10, MPO, MYD88, NLRP3, TNFα, and TLR1) and downregulation of 8 genes (CCR4, CCR6, CD40, CXCR3, FOXP3, GATA3, STAT4, and TBX21).

The miR-302 cluster transcriptionally regulated by POUV, SOX and STAT5B controls the undifferentiated state through the post-transcriptional repression of PBX3 and E2F7 in early chick development.

Early chick development is a systematic process governed by the concerted action of multiple mechanisms that regulate transcription and post-transcriptional processes. Post-transcriptional microRNA-mediated regulation, with regard to lineage specification and differentiation in early chick development, requires further investigation. Here, we characterize the transcriptional and post-transcriptional regulation mechanisms in undifferentiated chick blastodermal cells. Expression of the miR-302 cluster, POUV, SOX2, and STAT5B decreased in a time-dependent manner in early chick development. We found that POUV, SOX2, and STAT5B regulate the transcription of the miR-302 cluster, as its 5'-flanking region contains binding elements for each transcription factor. Additionally, POUV, SOX2, and STAT5B maintain pluripotency by regulating genes containing the miR-302 cluster target sequence. For example, microRNAs from the miR-302 cluster can bind to PBX3 and E2F7 transcripts, thus acting as a post-transcriptional regulator that maintains the undifferentiated state of blastodermal cells by balancing the expression of genes related to pluripotency and differentiation. Based on these results, we suggest that both transcriptional and post-transcriptional regulation of the miR302 cluster is critical for intrinsically controlling the undifferentiated state of chick embryonic blastodermal cells. These findings may help our understanding of the cellular and molecular mechanisms that underlie developmental decisions during early chick development.

QTL and gene expression analyses identify genes affecting carcass weight and marbling on BTA14 in Hanwoo (Korean Cattle).

Causal mutations affecting quantitative trait variation can be good targets for marker-assisted selection for carcass traits in beef cattle. In this study, linkage and linkage disequilibrium analysis (LDLA) for four carcass traits was undertaken using 19 markers on bovine chromosome 14. The LDLA analysis detected quantitative trait loci (QTL) for carcass weight (CWT) and eye muscle area (EMA) at the same position at around 50 cM and surrounded by the markers FABP4SNP2774C>G and FABP4_µsat3237. The QTL for marbling (MAR) was identified at the midpoint of markers BMS4513 and RM137 in a 3.5-cM marker interval. The most likely position for a second QTL for CWT was found at the midpoint of tenth marker bracket (FABP4SNP2774C>G and FABP4_µsat3237). For this marker bracket, the total number of haplotypes was 34 with a most common frequency of 0.118. Effects of haplotypes on CWT varied from a -5-kg deviation for haplotype 6 to +8 kg for haplotype 23. To determine which genes contribute to the QTL effect, gene expression analysis was performed in muscle for a wide range of phenotypes. The results demonstrate that two genes, LOC781182 (p = 0.002) and TRPS1 (p = 0.006) were upregulated with increasing CWT and EMA, whereas only LOC614744 (p = 0.04) has a significant effect on intramuscular fat (IMF) content. Two genetic markers detected in FABP4 were the most likely QTL position in this QTL study, but FABP4 did not show a significant effect on both traits (CWT and EMA) in gene expression analysis. We conclude that three genes could be potential causal genes affecting carcass traits CWT, EMA, and IMF in Hanwoo.

Use of a bovine genome array to identify new biological pathways for beef marbling in Hanwoo (Korean Cattle).

BACKGROUND: Marbling (intramuscular fat) is a valuable trait that impacts on meat quality and an important factor determining price of beef in the Korean beef market. Animals that are destined for this high marbling market are fed a high concentrate ration for approximately 30 months in the Korean finishing farms. However, this feeding strategy leads to inefficiencies and excessive fat production. This study aimed to identify candidate genes and pathways associated with intramuscular fat deposition on highly divergent marbling phenotypes in adult Hanwoo cattle. RESULTS: Bovine genome array analysis was conducted to detect differentially expressed genes (DEGs) in m. longissimus with divergent marbling phenotype (marbling score 2 to 7). Three data-processing methods (MAS5.0, GCRMA and RMA) were used to test for differential expression (DE). Statistical analysis identified 21 significant transcripts from at least two data-processing methods (P < 0.01). All 21 differentially expressed genes were validated by real-time PCR. Results showed a high concordance in the gene expression fold change between the microarrays and the real time PCR data. Gene Ontology (GO) and pathway analysis demonstrated that some genes (ADAMTS4, CYP51A and SQLE) over expressed in high marbled animals are involved in a protein catabolic process and a cholesterol biosynthesis process. In addition, pathway analysis also revealed that ADAMTS4 is activated by three regulators (IL-17A, TNFα and TGFß1). QRT-PCR was used to investigate gene expression of these regulators in muscle with divergent intramuscular fat contents. The results demonstrate that ADAMTS4 and TGFß1 are associated with increasing marbling fat. An ADAMTS4/TGFß1 pathway seems to be associated with the phenotypic differences between high and low marbled groups. CONCLUSIONS: Marbling differences are possibly a function of complex signaling pathway interactions between muscle and fat. These results suggest that ADAMTS4, which is involved in connective tissue degradation, could play a role in an important biological pathway for building up marbling in cattle. Moreover, ADAMTS4 and TGFß1could potentially be used as an early biological marker for marbling fat content in the early stages of growth.

Developmental and Degenerative Characterization of Porcine Parthenogenetic Fetuses during Early Pregnancy.

The difference between early pregnancy and delivery rate is quite large in assisted reproduction techniques (ARTs), including animal cloning. However, it is not clear why the implanted fetuses aborted after the early pregnancy stage. In the present study, we tried to evaluate the developmental and morphological characteristics of porcine parthenogenetically activated (PA) embryos or fetuses by electric stimulation during the early pregnancy period. The implanted PA and artificially inseminated (AI) embryos and fetuses were collected at day 26 and 35 after embryo transfer, respectively. The developmental and morphological parameters in the PA embryos at day 26 were similar to the AI embryos. The size, weight, formation of major organs, and apoptotic cells were not statistically different in both embryos at day 26. However, the PA fetuses at day 35 showed ceased fetal development and degenerated with abnormal morphologies in their organs. The day 35 PA fetuses showed significantly higher apoptotic cells and lower methylation status in three differentially methylated regions of the H19 gene compared to their comparators. Therefore, the normal development of PA embryos and fetuses during early gestation could lead to these pregnancies being misinterpreted as normal and become one of the main reasons for the gap between early pregnancy and delivery rate.

Characterization of the replication-competent porcine endogenous retrovirus class B molecular clone originated from Korean domestic pig.

Xenotransplantation from pigs offers an opportunity to resolve the shortage of human organs. The porcine endogenous retrovirus (PERV) cannot be eliminated because of its presence in the germline DNA. Three subgroups of the replication-competent PERV (PERV-A, PERV-B, and PERV-C) have been identified in pigs. We constructed a molecular clone of PERV-B from a Korean domestic pig BAC clone containing PERV genomes, and its replication competency was characterized in human cells. The pol region of PERV-B was detected in the genomic DNA of 293T cells transfected with PERV-B (465D1) and in the genomic DNA of 293T cells infected with PERV-B (465D1) viruses. The 293T cells transfected with PERV-B (465D1) were maintained for 140 days. PERV-B (465D1) showed the stationary replication competence with no toxicity in the cell growth.

Characterization of beef transcripts correlated with tenderness and moisture.

To identify transcriptional markers for beef traits related to meat tenderness and moisture, we measured the transcriptome of the Longissimus dorsi skeletal muscle in 10 Korean native cattle (KNC). We analyzed the correlation between the beef transcriptome and measurements of four different beef traits, shear force (SF), water holding capacity (WHC), cooking loss (CL), and loin eye area (LEA). We obtained non-overlapping and unique panels of genes showing strong correlations (|r|>0.8) with SF, WHC, CL, and LEA, respectively. Functional studies of these genes indicated that SF was mainly related to energy metabolism, and LEA to rRNA processing. Interestingly, our data suggested that WHC is influenced by protein metabolism. Overall, the skeletal muscle transcriptome pointed to the importance of energy and protein metabolism in determining meat quality after the aging process. The panels of transcripts for beef traits may be useful for predicting meat tenderness and moisture.

Proteins in longissimus muscle of Korean native cattle and their relationship to meat quality.

Proteomic profiling by two-dimensional gel electrophoresis and mass spectrometry of longissimus dorsi muscle tissue from Korean native cattle identified seven proteins that are differentially expressed in animals producing low and high quality grade beef. The expression level of alpha actin is increased in high quality grade beef and the expression levels of T-complex protein 1 (TCP-1), heat shock protein beta-1 (HSP27), and inositol 1,4,5-triphosphate receptor type1 (IP3R1), a new protein to be associated with meat quality, are increased in low quality grade beef. In particular, the quantitation of HSP27 and IP3R1 by both silver staining and immunoblotting correlated well with intramuscular fat content, meat tenderness, and free calcium levels. The data suggest that HSP27 and IP3R1 are potential meat quality biomarkers and their identification provides new insight into the molecular mechanisms and pathways associated with overall beef quality.

Artificial selection increased body weight but induced increase of runs of homozygosity in Hanwoo cattle.

Artificial selection has been demonstrated to have a rapid and significant effect on the phenotype and genome of an organism. However, most previous studies on artificial selection have focused solely on genomic sequences modified by artificial selection or genomic sequences associated with a specific trait. In this study, we generated whole genome sequencing data of 126 cattle under artificial selection, and 24,973,862 single nucleotide variants to investigate the relationship among artificial selection, genomic sequences and trait. Using runs of homozygosity detected by the variants, we showed increase of inbreeding for decades, and at the same time demonstrated a little influence of recent inbreeding on body weight. Also, we could identify ~0.2 Mb runs of homozygosity segment which may be created by recent artificial selection. This approach may aid in development of genetic markers directly influenced by artificial selection, and provide insight into the process of artificial selection.

An accurate method for quantifying and analyzing copy number variation in porcine KIT by an oligonucleotide ligation assay.

BACKGROUND: Aside from single nucleotide polymorphisms, copy number variations (CNVs) are the most important factors in susceptibility to genetic disorders because they affect expression levels of genes. In previous studies, pyrosequencing, mini-sequencing, real-time PCR, invader assays and other techniques have been used to detect CNVs. However, the higher the copy number in a genome, the more difficult it is to resolve the copies, so a more accurate method for measuring CNVs and assigning genotype is needed. RESULTS: PCR followed by a quantitative oligonucleotide ligation assay (qOLA) was developed for quantifying CNVs. The accuracy and precision of the assay were evaluated for porcine KIT, which was selected as a model locus. Overall, the root mean squares of bias and standard deviation of qOLA were 2.09 and 0.45, respectively. These values are less than half of those in the published pyrosequencing assay for analyzing CNV in porcine KIT. Using a combined method of qOLA and another pyrosequencing for quantitative analysis of KIT copies with spliced forms, we confirmed the segregation of KIT alleles in 145 F1 animals with pedigree information and verified the correct assignment of genotypes. In a diagnostic test on 100 randomly sampled commercial pigs, there was perfect agreement between the genotypes obtained by grouping observations on a scatter plot and by clustering using the nearest centroid sorting method implemented in PROC FASTCLUS of the SAS package. In a test on 159 Large White pigs, there were only two discrepancies between genotypes assigned by the two clustering methods (98.7% agreement), confirming that the quantitative ligation assay established here makes genotyping possible through the accurate measurement of high KIT copy numbers (>4 per diploid genome). Moreover, the assay is sensitive enough for use on DNA from hair follicles, indicating that DNA from various sources could be used. CONCLUSION: We have established a high resolution quantification method using an oligonucleotide ligation assay to measure CNVs, and verified the reliability of genotype assignment for random animal samples using the nearest centroid sorting method. This new method will make it more practical to determine KIT CNV and to genotype the complicated Dominant White/KIT locus in pigs. This procedure could have wide applications for studying gene or segment CNVs in other species.

Identification of differentially expressed genes related to intramuscular fat development in the early and late fattening stages of hanwoo steers.

Marbling of cattle meat is dependent on the coordinated expression of multiple genes. Cattle dramatically increase their intramuscular fat content in the longissimus dorsi muscle between 12 and 27 months of age. We used the annealing control primer (ACP)-differential display RT-PCR method to identify differentially expressed genes (DEGs) that may participate in the development of intramuscular fat between early (12 months old) and late fattening stages (27 months old). Using 20 arbitrary ACP primers, we identified and sequenced 14 DEGs. BLAST searches revealed that expression of the MDH, PI4-K, ferritin, ICER, NID-2, WDNMI, telethonin, filamin, and desmin (DES) genes increased while that of GAPD, COP VII, ACTA1, CamK II, and nebulin decreased during the late fattening stage. The results of functional categorization using the Gene Ontology database for 14 known genes indicated that MDH, GAPD, and COP VII are involved in metabolic pathways such as glycolysis and the TCA cycle, whereas telethonin, filamin, nebulin, desmin, and ACTA1 contribute to the muscle contractile apparatus, and PI4-K, CamK II, and ICER have roles in signal transduction pathways regulated by growth factor or hormones. The final three genes, NID-2, WDNMI, and ferritin, are involved in iron transport and extracellular protein inhibition. The expression patterns were confirmed for seven genes (MDH, PI4-K, ferritin, ICER, nebulin, WDNMI, and telethonin) using real-time PCR. We found that the novel transcription repressor ICER gene was highly expressed in the late fattening stage and during bovine preadipocyte differentiation. This information may be helpful in selecting candidate genes that participate in intramuscular fat development in cattle.

Genome-wide analysis of the diversity and ancestry of Korean dogs.

There are various hypotheses on dog domestication based on archeological and genetic studies. Although many studies have been conducted on the origin of dogs, the existing literature about the ancestry, diversity, and population structure of Korean dogs is sparse. Therefore, this study is focused on the origin, diversity and population structure of Korean dogs. The study sample comprised four major categories, including non-dogs (coyotes and wolves), ancient, modern and Korean dogs. Selected samples were genotyped using an Illumina CanineHD array containing 173,662 single nucleotide polymorphisms. The genome-wide data were filtered using quality control parameters in PLINK 1.9. Only autosomal chromosomes were used for further analysis. The negative off-diagonal variance of the genetic relationship matrix analysis depicted, the variability of samples in each population. FIS (inbreeding rate within a population) values indicated, a low level of inbreeding within populations, and the patterns were in concordance with the results of Nei's genetic distance analysis. The lowest FST (inbreeding rate between populations) values among Korean and Chinese breeds, using a phylogenetic tree, multi-dimensional scaling, and a TreeMix likelihood tree showed Korean breeds are highly related to Chinese breeds. The Korean breeds possessed a unique and large diversity of admixtures compared with other breeds. The highest and lowest effective population sizes were observed in Korean Jindo Black (485) and Korean Donggyeong White (109), respectively. The historical effective population size of all Korean dogs showed declining trend from the past to present. It is important to take immediate action to protect the Korean dog population while conserving their diversity. Furthermore, this study suggests that Korean dogs have unique diversity and are one of the basal lineages of East Asian dogs, originating from China.

Production of transgenic pig as an Alzheimer's disease model using a multi-cistronic vector system.

Alzheimer's disease (AD) is a progressive neurodegenerative disease associated with memory loss and cognitive impairments. An AD transgenic (Tg) pig model would be useful for preclinical testing of therapeutic agents. We generated an AD Tg pig by somatic cell nuclear transfer (SCNT) using a multi-cistronic vector that harbored three AD-related genes with a total of six well-characterized mutations: hAPP (K670N/M671L, I716V, and V717I), hTau (P301L), and hPS1 (M146V and L286P). Four AD Tg cell lines were established from Jeju black pig ear fibroblasts (JB-PEFs); the resultant JB-PEFAD cells harbored transgene integration, expressed transgene mRNAs, and had normal karyotypes. Tg line #2-1, which expressed high levels of the transgenes, was used for SCNT; cleavage and blastocyst rates of embryos derived from this line were lower than those of Non-Tg. These embryos yielded three piglets (Jeju National University AD-Tg pigs, JNUPIGs) revealed by microsatellite testing to be genetically identical to JB-PEFAD. Transgenes were expressed in multiple tissues, and at especially high levels in brain, and Aß-40/42, total Tau, and GFAP levels were high in brains of the Tg animals. Five or more copies of transgenes were inserted into chromosome X. This is the first report of an AD Tg pig derived from a multi-cistronic vector.

Confirming single nucleotide polymorphisms from expressed sequence tag datasets derived from three cattle cDNA libraries.

Using the Phred/Phrap/Polyphred/Consed pipeline established in the National Livestock Research Institute of Korea, we predicted candidate coding single nucleotide polymorphisms (cSNPs) from 7,600 expressed sequence tags (ESTs) derived from three cDNA libraries (liver, M. longissimus dorsi, and intermuscular fat) of Hanwoo (Korean native cattle) steers. From the 7,600 ESTs, 829 contigs comprising more than two EST reads were assembled using the Phrap assembler. Based on the contig analysis, 201 candidate cSNPs were identified in 129 contigs, in which transitions (69%) outnumbered transversions (31%). To verify whether the predicted cSNPs are real, 17 SNPs involved in lipid and energy metabolism were selected from the ESTs. Twelve of these were confirmed to be real while five were identified as artifacts, possibly due to expressed sequence tag sequence error. Further analysis of the 12 verified cSNPs was performed using the program BLASTX. Five were identified as nonsynonymous cSNPs, five were synonymous cSNPs, and two SNPs were located in 3'-UTRs. Our data indicated that a relatively high SNP prediction rate (71%) from a large EST database could produce abundant cSNPs rapidly, which can be used as valuable genetic markers in cattle.

A large-insert porcine library with sevenfold genome coverage: a tool for positional cloning of candidate genes for major quantitative traits.

A porcine genomic bacterial artificial chromosome (BAC) library was constructed by cloning partial EcoRI-digested high-molecular-weight DNA from a Korean native boar into the EcoRI site of the pBACe3.6 vector. The library consists of about 165,000 clones with an average insert size of 125 kb, representing about seven genome equivalents of coverage. About 130,000 clones (corresponding to fivefold genome coverage) were arrayed in 14 superpools which were organized as four dimensional pools. The library was further characterized by PCR screening of 38 microsatellite probes. An average of 4.84 positive clones were selected per marker. This indicates that the library is unbiased and will be useful for initiating fine scale physical mapping of major QTL in pigs. The library is being used to isolate specific clones by screening with type I and type II marker clones located in the QTL region affecting intramuscular fat content on SSC6.

Heat shock protein B1 and its regulator genes are negatively correlated with intramuscular fat content in the longissimus thoracis muscle of Hanwoo (Korean cattle) steers.

In previous proteomic studies, heat shock protein ß 1 (HSPB1) was detected as a candidate protein related to meat quality in cattle. This study sought to determine if its gene expression was associated with intramuscular fat content in the longissimus thoracis muscle of Korean cattle (Hanwoo). Tissue from two groups of 10 steers each, low-marbling (mean intramuscular fat content, 7.4 ± 1.5%) and high-marbling (23.5 ± 2.8%), were used for immunoblotting, real-time PCR, and statistical analyses. HSPB1 expression in both mRNA and protein was shown to be negatively related to intramuscular fat content (P < 0.05). Pathway analysis found two genes, TNF receptor superfamily member 6 (FAS) and angiotensinogen (AGT), that were regulators of the HSPB1 gene. The expression of the two genes showed a negative correlation with intramuscular fat content (P < 0.05). These results suggest that HSPB1, FAS, and AGT may be good candidate genes associated with intramuscular fat content in the longissimus muscle of Korean cattle.

Neuronal genes for subcutaneous fat thickness in human and pig are identified by local genomic sequencing and combined SNP association study.

Obesity represents a major global public health problem that increases the risk for cardiovascular or metabolic disease. The pigs represent an exceptional biomedical model related to energy metabolism and obesity in humans. To pinpoint causal genetic factors for a common form of obesity, we conducted local genomic de novo sequencing, 18.2 Mb, of a porcine QTL region affecting fatness traits, and carried out SNP association studies for backfat thickness and intramuscular fat content in pigs. In order to relate the association studies in pigs to human obesity, we performed a targeted genome wide association study for subcutaneous fat thickness in a cohort population of 8,842 Korean individuals. These combined association studies in human and pig revealed a significant SNP located in a gene family with sequence similarity 73, member A (FAM73A) associated with subscapular skin-fold thickness in humans (rs4121165, GC-corrected p-value  = 0.0000175) and with backfat thickness in pigs (ASGA0029495, p-value  = 0.000031). Our combined association studies also suggest that eight neuronal genes are responsible for subcutaneous fat thickness: NEGR1, SLC44A5, PDE4B, LPHN2, ELTD1, ST6GALNAC3, ST6GALNAC5, and TTLL7. These results provide strong support for a major involvement of the CNS in the genetic predisposition to a common form of obesity.

Effects of cloned-cattle meat diet on reproductive parameters in pregnant rabbits.

In this paper, we report on the effects of a diet containing cloned-cattle meat on the reproductive parameters in pregnant rabbits. The artificially inseminated rabbits (gestation day 0) were fed a diet containing 5% or 10% of normal or cloned-cattle meat during the gestation period. Rabbits fed commercial pellet (no additional supplementations) were used as the control. Supplementation of cloned-cattle meat diets did not have any toxicologically significant effects on reproductive performance in dams (body weight, clinical signs, organ weight, and cesarean section analysis). And it also did not affect on fetal development (body and placental weight, and external, visceral and skeletal findings) compared to the controls. The only difference was a food consumption in the first week of gestation for all meat-based diet groups (p<0.05, 0.01, and 0.001, respectively). Our results collectively suggest that there are no obvious differences in reproductive parameters in pregnant rabbits fed cloned-cattle meat.

Analysis of the molecular and regulatory properties of active porcine endogenous retrovirus gamma-1 long terminal repeats in kidney tissues of the NIH-Miniature pig.

The pig genome contains the gamma 1 family of porcine endogenous retroviruses (PERVs), which are a major obstacle to the development of successful xenotransplantation from pig to human. Long terminal repeats (LTRs) found in PERVs are known to be essential elements for the control of the transcriptional activity of single virus by different transcription factors (TFs). To identify transcribed PERV LTR elements, RT-PCR and DNA sequencing analyses were performed. Twenty-nine actively transcribed LTR elements were identified in the kidney tissues of the NIH-Miniature pig. These elements were divided into two major groups (I and II), and four minor groups (I-1, I-2, I-3, and II-1), by the presence of insertion and deletion (INDEL) sequences. Group I elements showed strong transcriptional activity compared to group II elements. Four different LTR elements (PL1, PL2, PL3, and PL4) as representative of the groups were analyzed by using a transient transfection assay. The regulation of their promoter activity was investigated by treatment with M.SssI (CpG DNA methyltransferase) and garcinol (histone acetyltransferase inhibitor). The transcriptional activity of PERV LTR elements was significantly reduced by treatment with M.SssI. These data indicate that transcribed PERV LTR elements harbor sufficient promoter activity to regulate the transcription of a single virus, and the transcriptional activity of PERV LTRs may be controlled by DNA methylation events.

The Pig Genome Database (PiGenome): an integrated database for pig genome research.

We established the Pig Genome Database (PiGenome) for pig genome research. The PiGenome integrates and analyzes all publicly available genome-wide data on pigs, including UniGenes, sequence tagged sites (STS) markers, quantitative trait loci (QTLs) data, and bacterial artificial chromosome (BAC) contigs. In addition, we produced 69,545 expressed sequence tags (ESTs) from the full-length enriched cDNA libraries of six tissues and 182 BAC contig sequences, which are also included in the database. QTLs, genetic markers, and BAC end-sequencing information were collected from public databases. The full-length enriched EST data were clustered and assembled into unique sequences, contigs, and singletons. The PiGenome provides functional annotation, identification of transcripts, mapping of coding sequences, and SNP information. It also provides an advanced search interface, a disease browser, alternative-splicing events, and a comparative gene map of the pig. A graphical map view and genome browser can map ESTs, contigs, BAC contigs (from the National Institute of Animal Science), Sino-Danish Pig Genome Project transcripts, and UniGene onto pig genome sequences which include our 182 BAC contigs and publically available BAC sequences of the Wellcome Trust Sanger Institute. The PiGenome is accessible at http://pigenome.nabc.go.kr/ .