EST-based gene discovery in pig: virtual expression patterns and comparative mapping to human

Genome-wide profiling of Sus scrofa circular RNAs across nine organs and three developmental stages

The spatio-temporal expression patterns of Circular RNA (circRNA) across organs and developmental stages are critical for its function and evolution analysis. However, they remain largely unclear in mammals. Here, we comprehensively analysed circRNAs in nine organs and three skeletal muscles of Guizhou miniature pig (S. scrofa), a widely used biomedical model animal. We identified 5,934 circRNAs and analysed their molecular properties, sequence conservation, spatio-temporal expression pattern, potential function, and interaction with miRNAs. S. scrofa circRNAs show modest sequence conservation with human and mouse circRNAs, are flanked by long introns, exhibit low abundance, and are expressed dynamically in a spatio-temporally specific manner. S. scrofa circRNAs show the greatest abundance and complexity in the testis. Notably, 31% of circRNAs harbour well-conserved canonical miRNA seed matches, suggesting that some circRNAs act as miRNAs sponges. We identified 149 circRNAs potentially associated with muscle growth and found that their host genes were significantly involved in muscle development, contraction, chromatin modification, cation homeostasis, and ATP hydrolysis-coupled proton transport; moreover, this set of genes was markedly enriched in genes involved in tight junctions and the calcium signalling pathway. Finally, we constructed the first public S. scrofa circRNA database, allowing researchers to query comprehensive annotation, expression, and regulatory networks of circRNAs.

Genital Chlamydia trachomatis: understanding the roles of innate and adaptive immunity in vaccine research

Chlamydia trachomatis is the leading cause of bacterial sexually transmitted disease worldwide, and despite significant advances in chlamydial research, a prophylactic vaccine has yet to be developed. This Gram-negative obligate intracellular bacterium, which often causes asymptomatic infection, may cause pelvic inflammatory disease (PID), ectopic pregnancies, scarring of the fallopian tubes, miscarriage, and infertility when left untreated. In the genital tract, Chlamydia trachomatis infects primarily epithelial cells and requires Th1 immunity for optimal clearance. This review first focuses on the immune cells important in a chlamydial infection. Second, we summarize the research and challenges associated with developing a chlamydial vaccine that elicits a protective Th1-mediated immune response without inducing adverse immunopathologies.

Activation of the transcription factor nuclear factor-kappa B in uterine luminal epithelial cells by interleukin 1 Beta 2: a novel interleukin 1 expressed by the elongating pig conceptus

Conceptus mortality is greatest in mammals during the peri-implantation period, a time when conceptuses appose and attach to the uterine surface epithelium while releasing proinflammatory molecules. Interleukin 1 beta (IL1B), a master proinflammatory cytokine, is released by the primate, rodent, and pig blastocyst during the peri-implantation period and is believed to be essential for establishment of pregnancy. The gene encoding IL1B has duplicated in the pig, resulting in a novel gene. Preliminary observations indicate that the novel IL1B is specifically expressed by pig conceptuses during the peri-implantation period. To verify this, IL1B was cloned from mRNA isolated from Day 12 pig conceptuses and compared with IL1B cloned from mRNA isolated from pig peripheral blood leukocytes (PBLs). The pig conceptuses, but not the PBLs, expressed a novel IL1B, referred to here as interleukin 1 beta 2 (IL1B2). Porcine endometrium was treated with recombinant porcine interleukin 1 beta 1 (IL1B1), the prototypical cytokine, and IL1B2 proteins. Immunohistochemistry and real-time RT-PCR were used to measure activation of nuclear factor-kappa B (NFKB) and NFKB-regulated transcripts, respectively, within the endometrium. Both IL1B1 and IL1B2 activated NFKB in the uterine luminal epithelium within 4 h. The NFKB activation and related gene expression, however, were lower in endometrium treated with IL1B2, suggesting that the conceptus-derived cytokine may have reduced activity within the uterus. In conclusion, the peri-implantation pig conceptus expresses a novel IL1B that can activate NFKB within the uterine surface epithelium, likely creating a proinflammatory microenvironment during establishment of pregnancy in the pig.

Genomic profiling to improve embryogenesis in the pig

Over the past decade the technology to characterize transcription during embryogenesis has progressed from estimating a single transcript to a reliable description of the entire transcriptome. Northern blots were followed by sequencing ESTs, quantitative real time PCR, cDNA arrays, custom oligo arrays, and more recently, deep sequencing. The amount of information that can be generated is overwhelming. The challenge now is how to glean information from these vast data sets that can be used to understand development and to improve methods for creating and culturing embryos in vitro, and for reducing reproductive loss. The use of ESTs permitted the identification of SPP1 as an oviductal component that could reduce polyspermy. Microarrays identified LDL and NMDA as components to replace BSA in embryo culture media. Deep sequencing implicated arginine, glycine, and folate as components that should be adjusted in our current culture system, and identified a characteristic of embryo metabolism that is similar to cancer and stem cells. Not only will these characterizations aid in improving in vitro production of embryos, but will also be useful for identifying, or creating conditions for donor cells that will be more likely to result in normal development of cloned embryos. The easily found targets have been identified, and now more sophisticated methods are being employed to advance our understanding of embryogenesis. Here the technology to study the global transcriptome is reviewed followed by specific examples of how the technology has been used to understand and improve porcine embryogenesis both in vitro and in vivo.

Endometrial gene expression profiling in pregnant Meishan and Yorkshire pigs on day 12 of gestation

Background

Litter size in pigs is a major factor affecting the profitability in the pig industry. The peri-implantation window in pigs is characterized by the coordinated interactions between the maternal uterine endometrium and the rapidly elongating conceptuses and represents a period of time during which a large percentage of the developing conceptuses are lost. However, the gene expression and regulatory networks in the endometrium contributing to the establishment of the maternal: placental interface remain poorly understood.

Results

We characterized the endometrial gene expression profile during the peri-implantation stage of development by comparing two breeds that demonstrate very different reproductive efficiencies. We employed the porcine Affymetrix GeneChip® to assay the transcriptomic profiles of genes expressed in the uterine endometrium obtained from Meishan and Yorkshire gilts (n = 4 for each breed) on day 12 of gestation (M12 and Y12, respectively). Total of 17,076 probesets were identified as "present" in at least two arrays. A mixed model-based statistical analysis predicted a total of 2,656 (q < 0.1) transcripts as differentially expressed between Meishan and Yorkshire pigs. Eighteen differentially expressed transcripts of interest were validated by quantitative real-time PCR. Gene ontology (GO) annotation revealed that the known functions of the differentially expressed genes were involved in a series of important biological processes relevant to early pregnancy establishment in the pig.

Conclusions

The results identified endometrial gene expression profiles of two breeds differing in litter size and identified candidate genes that are related to known physiological pathways related to reproductive prolificacy. These findings provide a deeper understanding of molecular pathways differing between two breeds at the critical peri-implantation stage of pregnancy, which can be utilized to better understand the events contributing to pregnancy establishment in the pig.

Animal models for studying female genital tract infection with Chlamydia trachomatis

Chlamydia trachomatis is a Gram-negative obligate intracellular bacterial pathogen. It is the leading cause of bacterial sexually transmitted disease in the world, with more than 100 million new cases of genital tract infections with C. trachomatis occurring each year. Animal models are indispensable for the study of C. trachomatis infections and the development and evaluation of candidate vaccines. In this paper, the most commonly used animal models to study female genital tract infections with C. trachomatis will be reviewed, namely, the mouse, guinea pig, and nonhuman primate models. Additionally, we will focus on the more recently developed pig model.

A gene expression atlas of the domestic pig

BACKGROUND

This work describes the first genome-wide analysis of the transcriptional landscape of the pig. A new porcine Affymetrix expression array was designed in order to provide comprehensive coverage of the known pig transcriptome. The new array was used to generate a genome-wide expression atlas of pig tissues derived from 62 tissue/cell types. These data were subjected to network correlation analysis and clustering.

RESULTS

The analysis presented here provides a detailed functional clustering of the pig transcriptome where transcripts are grouped according to their expression pattern, so one can infer the function of an uncharacterized gene from the company it keeps and the locations in which it is expressed. We describe the overall transcriptional signatures present in the tissue atlas, where possible assigning those signatures to specific cell populations or pathways. In particular, we discuss the expression signatures associated with the gastrointestinal tract, an organ that was sampled at 15 sites along its length and whose biology in the pig is similar to human. We identify sets of genes that define specialized cellular compartments and region-specific digestive functions. Finally, we performed a network analysis of the transcription factors expressed in the gastrointestinal tract and demonstrate how they sub-divide into functional groups that may control cellular gastrointestinal development.

CONCLUSIONS

As an important livestock animal with a physiology that is more similar than mouse to man, we provide a major new resource for understanding gene expression with respect to the known physiology of mammalian tissues and cells. The data and analyses are available on the websites http://biogps.org and http://www.macrophages.com/pig-atlas.

The minipig as nonrodent species in toxicology--where are we now?

Over the past 3 decades minipigs have moved from being an obscure alternative to dogs and nonhuman primates to being a standard animal model in regulatory toxicity studies. This article covers the use of minipigs as a model in the context of nonclinical drug safety and provides an overview of the minipig's developmental history and relates minipigs to other animal species commonly used in toxicology; and the minipig's translational power is supported by 43 case studies of marketed drug products covered. Special focus is given to criteria for selecting minipigs in nonclinical programs supporting the development of new medicines; the use of swine in the assessment of food additives, agrochemicals, and pesticides; as well as a regulatory perspective on the use of minipigs in Food and Drug Administration (FDA)-regulated products. This article presents the main points conveyed at a symposium held at the 2010 American College of Toxicology meeting in Baltimore, Maryland.

Validation of the Chlamydia trachomatis genital challenge pig model for testing recombinant protein vaccines

Chlamydia trachomatis is a Gram-negative obligate intracellular bacterial pathogen that is the leading cause of bacterial sexually transmitted disease in humans in developing countries. A vaccination programme is considered to be the best approach to reduce the prevalence of C. trachomatis infections. However, there are still no commercial C. trachomatis vaccines. In order to develop effective C. trachomatis vaccines, it is important to identify those antigens that elicit a protective immune response, and to develop new and adequate methods and adjuvants for effective vaccine delivery, as conventional methods have failed to induce protective immunity. In order to test different vaccine candidates, animal models are needed. Former studies have used non-primate monkeys, mice or guinea pig infection models. The present study used a pig model for testing recombinant protein vaccines. Two recombinant proteins, polymorphic membrane protein G (PmpG), and secretion and cellular translocation protein C (SctC), were tested for their ability to create protection in a pig C. trachomatis challenge model. The vaccines were administered subcutaneously with GNE adjuvant. Six weeks later, animals were challenged intravaginally with C. trachomatis serovar E. After a further 4 weeks, the pigs were euthanized. PmpG-immunized pigs were better protected than pigs immunized with the less promising SctC candidate vaccine antigen. Interestingly, significant protection was apparently not correlated with a strong humoral immune response upon subcutaneous immunization. In conclusion, the pig model is useful for studying the efficacy of vaccine candidates against genital human C. trachomatis infection.

Activation method does not alter abnormal placental gene expression and development in cloned pigs

Nuclear transfer efficiency is low and is thought to be caused by inadequate placental development. The objective of this study was to identify differentially expressed transcripts in pig placentas derived from in vivo fertilization, in vitro fertilization or nuclear transfer at Day 30 of gestation. Three activation methods were compared: electrical fusion/activation, electrical fusion/activation followed by treatment with reversible proteasomal inhibitor, MG132 or electrical fusion followed by activation with Thimerosal/DTT. Extraembryonic membranes were collected 30 days after artificial insemination (IVV) or embryo transfer (IVF and NT). Extraembryonic membrane cDNAs labeled with Cy5 and a reference cDNA labeled with Cy3 were hybridized to a pig reproductive tissue-specific 19,968 spot cDNA microarray. Images acquired and assessed by using Genepix Pro 4.0 were analyzed by Genespring 7.3.1. ANOVA (P < 0.05) identified 227 differentially expressed transcripts between the five treatments and 0 between the three activation methods. The nuclear transfer groups were pooled and compared to in vivo samples, identifying 34 up- and 19 down-regulated transcripts (>2-fold change, P < 0.05). Ten transcripts were validated by real-time PCR. UPTI, PAG2, and GLUD1 protein was quantified by Western blot and densitometry verified that UPTI and PAG2 proteins had an expression pattern that mirrored mRNA abundance (P < 0.05). Localization patterns were also determined for UPTI, PAG2, GLUD2 and 14-3-3 gamma in Day 35 extraembryonic membranes. Observed differences in gene and protein expression in nuclear transfer extraembryonic membranes indicate that an impaired fetal-maternal interface, and not the activation method, may be causing defects observed in cloned pigs.

Full-length enriched cDNA library construction from tissues related to energy metabolism in pigs

Genome sequencing of the pig is being accelerated because of its importance as an evolutionary and biomedical model animal as well as a major livestock animal. However, information on expressed porcine genes is insufficient to allow annotation and use of the genomic information. A series of expressed sequence tags of 5' ends of five full-length enriched cDNA libraries (SUSFLECKs) were functionally characterized. SUSFLECKs were constructed from porcine abdominal fat, induced fat cells, loin muscle, liver, and pituitary gland, and were composed of non-normalized and normalized libraries. A total of 55,658 ESTs that were sequenced once from the 5' ends of clones were produced and assembled into 17,684 unique sequences with 7,736 contigs and 9,948 singletons. In Gene Ontology analysis, two significant biological process leaf nodes were found: gluconeogenesis and translation elongation. In functional domain analysis based on the Pfam database, the beta transducin repeat domain of WD40 protein was the most frequently occurring domain. Twelve genes, including SLC25A6, EEF1G, EEF1A1, COX1, ACTA1, SLA, and ANXA2, were significantly more abundant in fat tissues than in loin muscle, liver, and pituitary gland in the SUSFLECKs. These characteristics of SUSFLECKs determined by EST analysis can provide important insight to discover the functional pathways in gene networks and to expand our understanding of energy metabolism in the pig.

A robust linkage map of the porcine autosomes based on gene-associated SNPs

BACKGROUND

Genetic linkage maps are necessary for mapping of mendelian traits and quantitative trait loci (QTLs). To identify the actual genes, which control these traits, a map based on gene-associated single nucleotide polymorphism (SNP) markers is highly valuable. In this study, the SNPs were genotyped in a large family material comprising more than 5,000 piglets derived from 12 Duroc boars crossed with 236 Danish Landrace/Danish Large White sows. The SNPs were identified in sequence alignments of 4,600 different amplicons obtained from the 12 boars and containing coding regions of genes derived from expressed sequence tags (ESTs) and genomic shotgun sequences.

RESULTS

Linkage maps of all 18 porcine autosomes were constructed based on 456 gene-associated and six porcine EST-based SNPs. The total length of the averaged-sex whole porcine autosome was estimated to 1,711.8 cM resulting in an average SNP spacing of 3.94 cM. The female and male maps were estimated to 2,336.1 and 1,441.5 cM, respectively. The gene order was validated through comparisons to the cytogenetic and/or physical location of 203 genes, linkage to evenly spaced microsatellite markers as well as previously reported conserved synteny. A total of 330 previously unmapped genes and ESTs were mapped to the porcine autosome while ten genes were mapped to unexpected locations.

CONCLUSION

The linkage map presented here shows high accuracy in gene order. The pedigree family network as well as the large amount of meiotic events provide good reliability and make this map suitable for QTL and association studies. In addition, the linkage to the RH-map of microsatellites makes it suitable for comparison to other QTL studies.

Chlamydia trachomatis infection: host immune responses and potential vaccines

Chlamydia trachomatis causes genital tract infections that affect men, women, and children on a global scale. This review focuses on innate and adaptive immune responses in the female reproductive tract (FRT) to genital tract infections with C. trachomatis. It covers C. trachomatis infections and highlights our current knowledge of genital tract infections, serovar distribution, infectious load, and clinical manifestations of these infections in women. The unique features of the immune system of the FRT will be discussed and will include a review of our current knowledge of innate and adaptive immunity to chlamydial infections at this mucosal site. The use of animal models to study the pathogenesis of, and immunity to, Chlamydia infection of the female genital tract will also be discussed and a review of recent immunization and challenge experiments in the murine model of chlamydial FRT infection will be presented.

A pig multi-tissue normalised cDNA library: large-scale sequencing, cluster analysis and 9K micro-array resource generation

Background

Domestic animal breeding and product quality improvement require the control of reproduction, nutrition, health and welfare in these animals. It is thus necessary to improve our knowledge of the major physiological functions and their interactions. This would be greatly enhanced by the availability of expressed gene sequences in the databases and by cDNA arrays allowing the transcriptome analysis of any function.

The objective within the AGENAE French program was to initiate a high-throughput cDNA sequencing program of a 38-tissue normalised library and generate a diverse microarray for transcriptome analysis in pig species.

Results

We constructed a multi-tissue cDNA library, which was normalised and subtracted to reduce the redundancy of the clones. Expressed Sequence Tags were produced and 24449 high-quality sequences were released in EMBL database. The assembly of all the public ESTs (available through SIGENAE website) resulted in 40786 contigs and 54653 singletons. At least one Agenae sequence is present in 11969 contigs (12.5%) and in 9291 of the deeper-than-one-contigs (22.8%). Sequence analysis showed that both normalisation and subtraction processes were successful and that the initial tissue complexity was maintained in the final libraries. A 9K nylon cDNA microarray was produced and is available through CRB-GADIE. It will allow high sensitivity transcriptome analyses in pigs.

Conclusion

In the present work, a pig multi-tissue cDNA library was constructed and a 9K cDNA microarray designed. It contributes to the Expressed Sequence Tags pig data, and offers a valuable tool for transcriptome analysis.

Advances in swine transcriptomics

The past five years have seen a tremendous rise in porcine transcriptomic data. Available porcine Expressed Sequence Tags (ESTs) have expanded greatly, with over 623,000 ESTs deposited in Genbank. ESTs have been used to expand the pig-human comparative maps, but such data has also been used in many ways to understand pig gene expression. Several methods have been used to identify genes differentially expressed (DE) in specific tissues or cell types under different treatments. These include open screening methods such as suppression subtractive hybridization, differential display, serial analysis of gene expression, and EST sequence frequency, as well as closed methods that measure expression of a defined set of sequences such as hybridization to membrane arrays and microarrays. The use of microarrays to begin large-scale transcriptome analysis has been recently reported, using either specialized or broad-coverage arrays. This review covers published results using the above techniques in the pig, as well as unpublished data provided by the research community, and reports on unpublished Affymetrix data from our group. Published and unpublished bioinformatics efforts are discussed, including recent work by our group to integrate two broad-coverage microarray platforms. We conclude by predicting experiments that will become possible with new anticipated tools and data, including the porcine genome sequence. We emphasize that the need for bioinformatics infrastructure to efficiently store and analyze the expanding amounts of gene expression data is critical, and that this deficit has emerged as a limiting factor for acceleration of genomic understanding in the pig.

Radiation hybrid mapping of 18 positional and physiological candidate genes for arthrogryposis multiplex congenita on porcine chromosome 5

We report the chromosomal assignment of 18 porcine genes to human homologues using the INRA-Minnesota swine radiation hybrid panel (IMpRH). These genes (CACNA1C, COL2A1, CPNE8, C3F, C12ORF4, DDX11, GDF11, HOXC8, KCNA1, MDS028, TMEM106C, NR4A1, PHB2, PRICKLE1, Q6ZUQ4, SCN8A, TUBA8 and USP18) are located on porcine chromosome 5 (SSC5) and represent positional and functional candidates for arthrogryposis multiplex congenita (AMC), which maps to SSC5. CPNE8, PRICKLE1, Q6ZUQ4 and TUBA8 were mapped to the interval for pig AMC between microsatellites SW152 and SW904. Three SNPs in TUBA8 co-segregated with the AMC phenotype in 230 pigs of our research population without recombination and could be used as a genetic marker test for AMC. In addition, we provide evidence that a small chromosomal region of HSA22q11.2 evolutionarily corresponds to SSC5q12-q22 (and contains the human homologues of porcine SW152, Q6ZUQ4, TUBA8 and USP18), while the regions flanking HSA22q11.2 on SSC5 correspond to HSA12p13 and HSA12q12. We identified seven distinct chromosomal blocks, further supporting extensive rearrangements between genes on HSA12 and HSA22 in the AMC region on SSC5.

A wing expressed sequence tag resource for Bicyclus anynana butterflies, an evo-devo model

BACKGROUND

Butterfly wing color patterns are a key model for integrating evolutionary developmental biology and the study of adaptive morphological evolution. Yet, despite the biological, economical and educational value of butterflies they are still relatively under-represented in terms of available genomic resources. Here, we describe an Expression Sequence Tag (EST) project for Bicyclus anynana that has identified the largest available collection to date of expressed genes for any butterfly.

RESULTS

By targeting cDNAs from developing wings at the stages when pattern is specified, we biased gene discovery towards genes potentially involved in pattern formation. Assembly of 9,903 ESTs from a subtracted library allowed us to identify 4,251 genes of which 2,461 were annotated based on BLAST analyses against relevant gene collections. Gene prediction software identified 2,202 peptides, of which 215 longer than 100 amino acids had no homology to any known proteins and, thus, potentially represent novel or highly diverged butterfly genes. We combined gene and Single Nucleotide Polymorphism (SNP) identification by constructing cDNA libraries from pools of outbred individuals, and by sequencing clones from the 3' end to maximize alignment depth. Alignments of multi-member contigs allowed us to identify over 14,000 putative SNPs, with 316 genes having at least one high confidence double-hit SNP. We furthermore identified 320 microsatellites in transcribed genes that can potentially be used as genetic markers.

CONCLUSION

Our project was designed to combine gene and sequence polymorphism discovery and has generated the largest gene collection available for any butterfly and many potential markers in expressed genes. These resources will be invaluable for exploring the potential of B. anynana in particular, and butterflies in general, as models in ecological, evolutionary, and developmental genetics.

Generation and analysis of large-scale expressed sequence tags (ESTs) from a full-length enriched cDNA library of porcine backfat tissue

Background

Genome research in farm animals will expand our basic knowledge of the genetic control of complex traits, and the results will be applied in the livestock industry to improve meat quality and productivity, as well as to reduce the incidence of disease. A combination of quantitative trait locus mapping and microarray analysis is a useful approach to reduce the overall effort needed to identify genes associated with quantitative traits of interest.

Results

We constructed a full-length enriched cDNA library from porcine backfat tissue. The estimated average size of the cDNA inserts was 1.7 kb, and the cDNA fullness ratio was 70%. In total, we deposited 16,110 high-quality sequences in the dbEST division of GenBank (accession numbers: DT319652-DT335761). For all the expressed sequence tags (ESTs), approximately 10.9 Mb of porcine sequence were generated with an average length of 674 bp per EST (range: 200–952 bp). Clustering and assembly of these ESTs resulted in a total of 5,008 unique sequences with 1,776 contigs (35.46%) and 3,232 singleton (65.54%) ESTs. From a total of 5,008 unique sequences, 3,154 (62.98%) were similar to other sequences, and 1,854 (37.02%) were identified as having no hit or low identity (<95%) and 60% coverage in The Institute for Genomic Research (TIGR) gene index of Sus scrofa. Gene ontology (GO) annotation of unique sequences showed that approximately 31.7, 32.3, and 30.8% were assigned molecular function, biological process, and cellular component GO terms, respectively. A total of 1,854 putative novel transcripts resulted after comparison and filtering with the TIGR SsGI; these included a large percentage of singletons (80.64%) and a small proportion of contigs (13.36%).

Conclusion

The sequence data generated in this study will provide valuable information for studying expression profiles using EST-based microarrays and assist in the condensation of current pig TCs into clusters representing longer stretches of cDNA sequences. The isolation of genes expressed in backfat tissue is the first step toward a better understanding of backfat tissue on a genomic basis.

Specific-pathogen-free pigs as an animal model for studying Chlamydia trachomatis genital infection

The purpose of the present study was to evaluate pigs as a large-animal model for female genital infection with two Chlamydia trachomatis human serovar E strains. Sixteen-week-old specific-pathogen-free female pigs (gilts) were intravaginally infected with the trachoma type E reference strain Bour or the urogenital serovar E strain 468. Several conclusions can be drawn from our findings on the pathogenicity of a primary C. trachomatis genital infection in gilts. First of all, we demonstrated that the serovar E strains Bour and 468 could ascend in the genital tract of gilts. The serovar E strains could replicate in the superficial columnar cervical epithelium and in the superficial epithelial layer of the uterus, which are known to be the specific target sites for a C. trachomatis genital infection in women. Second, inflammation and pathology occurred at the replication sites. Third, the organisms could trigger a humoral immune response, as demonstrated by the presence of immunoglobulin M (IgM), IgG, and IgA in both serum and genital secretion samples. Our findings imply that the pig model might be useful for studying the pathology, pathogenesis, and immune response to a C. trachomatis infection of the genital system.

Advances in livestock genomics: Opening the barn door

Genome research in animals used in agriculture has progressed rapidly in recent years, moving from rudimentary genome maps to trait maps to gene discovery. These advances are the result of animal genome projects following closely in the footsteps of the Human Genome Project, which has opened the door to genome research in farm animals. In return, genome research in livestock species is contributing to our understanding of chromosome evolution and to informing the human genome. Enhancement of these contributions plus the much anticipated application of DNA-based tools to animal health and production can be expected as livestock genomics enters its sequencing era.

Chorionic mRNA expression and N-glycodiversity of pregnancy-associated glycoprotein family (PAG) of the European bison (Bison bonasus)

Placental PAG mRNA expression and N-glycodiversity of multiple PAG proteins secreted in vitro by trophectoderm (chorion epithelium) of wild pecoran Bovidae taxons was not examined previously. The study on European bison (Eb) aimed: (1) to determine placental PAG mRNA expression by in situ hybridisation; (2) to identify a profile of pecoran PAG protein family secreted in vitro by cotyledonary (CT) explants; (3) to examine N-glycodiversity of the PAG proteins in this wild taxon. In addition, we compared (4) a profile and N-glycodiversity of the PAG protein family secreted in vitro by CT and interCT-trophectoderm (intCT-TRD) explants of domestic ruminants. Cotyledonary sections of the Eb were used for in situ hybridisation (ISH) with (35)S-labelled probes produced with porcine PAG cDNA as templates. Various CT and intCT-TRD explants were long-term cultured in vitro. Chorionic proteins were isolated from media, ultra-filtrated (>10 kDa MWCO) and analysed by PAGE-Western blotting with various polyclonal anti-PAG sera. Protein samples with or without enzymatic deglycosylation were examined after different times of explant cultures. Released chorionic proteins were deglycosylated by N-glycanase F (PNGase F+) and compared to glycosylated forms (PNGase F-). This is the first paper demonstrating the PAG-like mRNA transcript expression (by ISH) and N-glycodiversity of immuno-reactive PAG-like proteins (produced in vitro by chorionic explants) of European bison. Various PAG proteins of Eb (EbPAG) were secreted by CT explants during long-term in vitro studies. Major approximately 78, approximately 67 and approximately 65 kDa EbPAG-like proteins were reduced by enzymatic deglycosylation (at least by 10 kDa). Considerably smaller amounts of approximately 45 kDa EbPAG-like proteins were also observed. In addition, we have found that various PAG proteins (30-73 kDa) were secreted by bovine CT explants, during long-term in vitro cultures. Corresponding amounts of PAG proteins, similar in M(r), were also secreted by intCT-TRD explants, whose tissues were not utilised for PAG protein extraction during other scientists' previous studies. It seems that the M(r)-heterogeneity and N-glycodiversity of the PAG protein family can play very important role during feto-placental interactions in Bovidae species.

Farm animal genomics and informatics: an update

Farm animal genomics is of interest to a wide audience of researchers because of the utility derived from understanding how genomics and proteomics function in various organisms. Applications such as xenotransplantation, increased livestock productivity, bioengineering new materials, products and even fabrics are several reasons for thriving farm animal genome activity. Currently mined in rapidly growing data warehouses, completed genomes of chicken, fish and cows are available but are largely stored in decentralized data repositories. In this paper, we provide an informatics primer on farm animal bioinformatics and genome project resources which drive attention to the most recent advances in the field. We hope to provide individuals in biotechnology and in the farming industry with information on resources and updates concerning farm animal genome projects.

Physical mapping of genes in the porcine ovarian transcriptome

Reproductive efficiency and associated traits are of major economic importance to the swine industry and have been more difficult to improve genetically than other production traits. Integration of phenotypical data with gene mapping and expression studies provides a powerful approach for dissection of the genetic basis regulating complex traits. We developed a total of 101 polymerase chain reaction-based markers, representing 91 unique genes, for expressed sequence tags previously reported to be putatively differentially expressed in the porcine ovarian transcriptome of a swine line selected on an index of high ovulation rate and embryonic survival. These were subsequently used in physical mapping experiments with a porcine radiation hybrid and somatic cell hybrid panels. Our results increased the information content of the porcine physical map useful for comparative mapping by c. 10%. Moreover, the mapped genes are likely to be biologically relevant to the molecular mechanisms that control ovulation rate in the pig. A total of 12 differentially expressed genes were mapped to regions previously reported to contain quantitative trait loci affecting swine ovulation rate.

EST-based analysis of gene expression in the porcine brain

Since pig is an important livestock species worldwide, its gene expression has been investigated intensively, but rarely in brain. In order to study gene expression profiles in the pig central nervous system, we sequenced and analyzed 43,122 high-quality 5’ end expressed sequence tags (ESTs) from porcine cerebellum, cortex cerebrum, and brain stem cDNA libraries, involving several different prenatal and postnatal developmental stages. The initial ESTs were assembled into 16,101 clusters and compared to protein and nucleic acid databases in GenBank. Of these sequences, 30.6% clusters matched protein databases and represented function known sequences; 75.1% had significant hits to nucleic acid databases and partial represented known function; 73.3% matched known porcine ESTs; and 21.5% had no matches to any known sequences in GenBank. We used the categories defined by the Gene Ontology to survey gene expression in the porcine brain.

Mapping of 443 porcine EST improves the comparative maps for SSC1 and SSC7 with the human genome

Numerous mapping studies of complex traits in the pig have resulted in quantitative trait loci (QTL) intervals of 10-20 cM. To improve the chances to identify the genes located in such intervals, increased expressed sequence tags (EST)-based marker density, coupled with comparative mapping with species whose genomes have been sequenced such as human and mouse, is the most efficient tool. In this study, we mapped 443 porcine EST with a radiation hybrid (RH) panel (384 had LOD > 6.0) and a somatic cell hybrid panel. Requiring no discrepancy between two-point and multipoint RH data allowed robust assignment of 309 EST, of which most were located on porcine chromosomes (SSC) 1, 4, 7, 8 and X. Moreover, we built framework maps for two chromosomes, SSC1 and SSC7, with mapped QTL in regions with known rearrangement between pig and human genomes. Using the Blast tool, we found orthologies between 407 of the 443 pig cDNA sequences and human genes, or to existing pig genes. Our porcine/human comparative mapping results reveal possible new homologies for SSC1, SSC3, SSC5, SSC6, SSC12 and SSC14 and add markers in synteny breakpoints for chromosome 7.

Transcriptional Profiling of Pig Embryogenesis by Using a 15-K Member Unigene Set Specific for Pig Reproductive Tissues and Embryos1

Differential mRNA expression patterns were evaluated between germinal vesicle oocytes (pgvo), four-cell (p4civv), blastocyst (pblivv), and in vitro-produced four-cell (p4civp) and in vitro-produced blastocyst (pblivp) stage embryos to determine key transcripts responsible for early embryonic development in the pig. Five comparisons were made: pgvo to p4civv, p4civv to pblivv, pgvo to pblivv, p4civv to p4civp, and pblivv to pblivp. ANOVA (P < 0.05) was performed with the Benjamini and Hochberg false-discovery-rate multiple correction test on each comparison. A comparison of pgvo to p4civv, p4civv to pblivv, and pgvo to pblivv resulted in 3214, 1989, and 4528 differentially detected cDNAs, respectively. Real-time PCR analysis on seven transcripts showed an identical pattern of changes in expression as observed on the microarrays, while one transcript deviated at a single cell stage. There were 1409 and 1696 differentially detected cDNAs between the in vitro- and in vivo-produced embryos at the four-cell and blastocyst stages, respectively, without the Benjamini and Hochberg false-discovery-rate multiple correction test. Real-time polymerase chain reaction (PCR) analysis on four genes at the four-cell stage showed an identical pattern of gene expression as found on the microarrays. Real-time PCR analysis on four of five genes at the blastocyst stage showed an identical pattern of gene expression as found on the microarrays. Thus, only 1 of the 39 comparisons of the pattern of gene expression exhibited a major deviation between the microarray and the real-time PCR. These results illustrate the complex mechanisms involved in pig early embryonic development.

Transcriptional profiling using a novel cDNA array identifies differential gene expression during porcine embryo elongation

A novel porcine cDNA array, containing 1,015 PCR products selected for embryonic expression, was used for transcriptional profiling of conceptuses at four stages of peri-implantation development. Total conceptus RNA from small spherical, large spherical, tubular, and filamentous stages was amplified, converted to cDNA, and hybridized to membranes. Initially, normalized signal intensities obtained using cDNA from total RNA or from amplified RNA were compared. Uniform distribution of P-values associated with t-tests conducted for each gene indicated no evidence that amplification introduced bias. Analysis of data obtained by using amplified targets and the novel array identified genes differentially expressed across stages. Such genes were identified by testing for significant stage effects in gene-specific mixed models. A total of nine genes were declared differentially expressed. Six of the nine genes had P-values less than 0.001, and a false discovery rate of approximately 17% was associated with this significance threshold. Two out of six genes were significant when using the Bonferroni method to control the probability of one or more false positives. The other three genes had P-values between 0.001 and 0.01 and exhibited differences greater than twofold between stages. All four genes selected for confirmation (steroidogenic acute regulatory protein, interleukin 1 beta, transforming growth factor beta 3, and thymosin beta 10) were shown to be differentially expressed by using quantitative real time RT-PCR. Our study shows that RNA amplification is useful for transcriptional profiling with limiting porcine embryonic RNA, and that this novel targeted array can detect differential gene expression during trophoblastic elongation. Finally, our results contribute to an increased understanding of the temporal patterns of expression of known genes controlling conceptus development, as well as identify novel genes also differentially regulated during implantation.

Large-Scale Generation and Analysis of Expressed Sequence Tags from Porcine Ovary1

One method to identify the factors that control ovarian function is to characterize the genes that are expressed in ovary. In the present study, cDNA libraries from fetal, neonatal, and prepubertal porcine ovaries, pubertal ovaries on different days of the estrous cycle (Days 0 [follicle], 5, and 12 [follicle and corpus luteum]), and follicles isolated from weaned sows (diameter, 2, 4, 6, and 8 mm) were constructed and sequenced. A total of 22 176 cDNAs were sequenced, of which 15 613 were of sufficient quality for clustering. Clustering of cDNAs resulted in 8507 contigs, 6294 (74%) of which were comprised of a single sequence. Sixty-eight percent of the contigs had consensus sequences that were homologous to existing Tentative Consensus (TC) sequences or mature transcripts (ET) in The Institute for Genomic Research Porcine Gene Index. The consensus sequences were classified according to the Gene Ontology Index. Most cDNA-encoded proteins were components of the nucleus, ribosome, or mitochondrion. The proteins primarily functioned in binding, catalysis, and transport. Nearly 75% of the proteins were involved in metabolism and cell growth and/or maintenance. Analysis of the cDNA frequency across different libraries demonstrated differential gene expression within different-size follicles, between follicles and corpora lutea, and across developmental time-points. The expression of selected genes (analyzed by ribonuclease protection assay and Northern blotting) was consistent with the frequency of their respective cDNA in the individual libraries. This porcine ovary unigene set will be useful for identifying factors and mechanisms controlling ovarian follicular development in a variety of species.

Comparative mapping of human chromosome 10 to pig chromosomes 10 and 14

Identification of predictive markers in QTL regions that impact production traits in commercial populations of swine is dependent on construction of dense comparative maps with human and mouse genomes. Chromosomal painting in swine suggests that large genomic blocks are conserved between pig and human, while mapping of individual genes reveals that gene order can be quite divergent. High-resolution comparative maps in regions affecting traits of interest are necessary for selection of positional candidate genes to evaluate nucleotide variation causing phenotypic differences. The objective of this study was to construct an ordered comparative map of human chromosome 10 and pig chromosomes 10 and 14. As a large portion of both pig chromosomes are represented by HSA10, genes at regularly spaced intervals along this chromosome were targeted for placement in the porcine genome. A total of 29 genes from human chromosome 10 were mapped to porcine chromosomes 10 (SSC10) and 14 (SSC14) averaging about 5 Mb distance of human DNA per marker. Eighteen genes were assigned by linkage in the MARC mapping population, five genes were physically assigned with the IMpRH mapping panel and seven genes were assigned on both maps. Seventeen genes from human 10p mapped to SSC10, and 12 genes from human 10q mapped to SSC14. Comparative maps of mammalian species indicate that chromosomal segments are conserved across several species and represent syntenic blocks with distinct breakpoints. Development of comparative maps containing several species should reveal conserved syntenic blocks that will allow us to better define QTL regions in livestock.

High-resolution comparative mapping of pig Chromosome 4, emphasizing the FAT1 region

The first quantitative trait locus (QTL) in pigs, FAT1, was found on Chromosome 4 (SSC4) using a Wild Boar intercross. Further mapping has refined the FAT1 QTL to a region with conserved synteny to both human Chromosomes 1 and 8. To both improve the comparative map of the entire SSC4 and to define the specific human chromosome region with conserved synteny to FAT1, we have now mapped 103 loci to pig Chromosome 4 using a combination of radiation hybrid and linkage mapping. The physical data and linkage analysis results are in very good agreement. Comparative analysis revealed that gene order is very well conserved across SSC4 compared to both HSA1 and HSA8. The breakpoint in conserved synteny was refined to an area of about 23 cR on the q arm of SSC4 corresponding to a genetic distance of less than 0.5 cM. Localizations of the centromeres do not seem to have been conserved between the two species. No remnants of the HSA1 centromere were detected on the corresponding region on SSC4 and traces from the centromeric region of SSC4 cannot clearly be revealed on the homologous region on HSA8. This refined SSC4 map and the comparative analysis will be a great aid in the search for the genes underlying the FAT1 locus.

Gene expression profiling of the bovine gastrointestinal tract

Basal gene expression levels across the bovine gastrointestinal tract (GI) were examined in an attempt to formulate genetic explanations for the differences in function that are known or thought to exist between the various regions. Gene expression along the tract was studied through the random sequencing of a total of 16 412 clones from seven tissue-specific cDNA libraries spanning its length. The expressed sequence tags (ESTs) within each library were clustered to reduce clone redundancy and obtain longer consensus sequences. BLASTN and BLASTX searches against the NCBI human RefSeq databases were used to find putative matches for the bovine sequences and gene ontology assignments were made. Notable similarities and differences in gene expression were observed among the various compartments of the GI tract of the bovine. Many of the prominent transcripts have yet to be reliably identified and the prominence of others may be worthy of further examination. This collection of ESTs represents an important resource for the future construction of a GI tract specific microarray for further gene expression studies.Key words: gene expression, EST, gastrointestinal tract, cattle.

PLET1 (C11orf34), a highly expressed and processed novel gene in pig and mouse placenta, is transcribed but poorly spliced in human

Sequencing of porcine cDNAs identified a novel EST with high frequency in placenta tissue. Full-length PLET1 (placenta-expressed transcript 1, also called C11orf34) matched a mouse cDNA and many bovine and mouse ESTs but no human transcripts or ESTs. However, the porcine cDNA matched several putative exons within a human genomic DNA fragment on chromosome 11. This human locus is in a region of conserved synteny with pig chromosome 9, to which the porcine gene was subsequently mapped. RNA blot hybridization showed that this gene had high expression in porcine and mouse conceptus and throughout placenta development. In situ hybridization using mouse placenta showed PLET1 expression in trophoblast cells of the labyrinth, as well as in spongiotrophoblast and glycogen trophoblast cells. However, no expression of PLET1 was detected by RNA blot analysis of human placenta, although RT-PCR analysis detected very small amounts of partially spliced RNA that were significantly less abundant than the RNA levels in mouse placenta. Donor and acceptor splicing site sequences in the exons of the human gene are poorly conserved and may be the cause of inefficient splicing found specifically in human tissue. Our data correct GenomeScan annotation of this region of the human genome and describe functional gene discovery in mammals not recognized in human EST projects.

Embryogenomics of pre-implantation mammalian development: current status

Pre-implantation development is marked by many critical molecular events, including the maternal to zygotic transition and the first differentiation of cells. Understanding such events is important, for both basic reproductive biology and practical applications, including regenerative medicine and livestock production. Scarcity of materials has hampered the progress of the field, but systematic genomics approaches are beginning to be applied to the study of pre-implantation development, resulting in unprecedented amounts of data about the pre-implantation process. The first step in embryogenomics is to collect and sequence cDNAs (expressed sequence tags (ESTs)) for genes that are expressed and function in these early embryos. Mouse work is the most advanced, with 140111 ESTs derived from all stages of pre-implantation development currently available in the public sequence database. For other mammals, at present only approximately 1000 ESTs can be found in the public database, but efforts by several groups are generating cDNA libraries and ESTs. In the present review, the current status of the implementation of these investigative tools for mammalian pre-implantation embryos is discussed.