Analysis of a full-length cDNA library constructed from swine olfactory bulb for elucidation of expressed genes and their transcription initiation sites

Microduplication of Xq24 and Hartsfield syndrome with holoprosencephaly, ectrodactyly, and clefting

The combination of holoprosencephaly and ectrodactyly, also known as Hartsfield syndrome, represents a unique genetic entity. An X-linked recessive mode of transmission has been suggested for this condition based on the observation that male patients have preferentially been affected. Thus far, no candidate genes have been suggested on the X chromosome. We report a male patient with a full-blown Hartsfield syndrome phenotype who had microduplication at Xq24 involving four genes. He presented with bilateral ectrodactyly of the hands (both hands had four fingers with a deep gap between the 2nd and 3rd digits), cleft lip and palate, and a depressed nasal bridge. Magnetic resonance imaging of the brain revealed lobar holoprosencephaly. His G-banded karyotype was normal. Array comparative genomic hybridization (CGH) using the Agilent 244K Whole Human Genome CGH array revealed a microduplication at Xq24 of 210 kb. Parental testing revealed that the deletion was derived from the asymptomatic mother. Of the genes on the duplicated interval, the duplications of SLC25A43 and SLC25A5 appeared to be the most likely to explain the patient's phenotype. From a clinical standpoint, it is important to point out that the propositus, who performs relatively well with holoprosencephaly and has a developmental quotient around 70, has survived multiple life-threatening episodes of hypernatremia. Awareness of the risk of hypernatremia is of great importance for the anticipatory management of patients with ectrodactyly and an oral cleft, even in the absence of overt hypotelorism.

Molecular characterization and tissue expression of ovine PSAM6 gene from muscle full-length cDNA library of black-boned sheep

An ovine PSMA6 gene was obtained from muscle full-length cDNA library of black-boned sheep. The sequences for the PSAM6 gene of Romney sheep and Yunling black goat were also generated in this study. Sequence analysis revealed that nucleotide sequence of this gene was not homologous to any of the known sheep genes, and its open reading frame encodes a protein that contains the putative conserved domain of proteasome subunit alpha type 6 (PSAM6). The nucleotide sequence had higher identity with other animals. However, one mutation of A to G at the site of 383 bp, leading to an amino acid mutation of Asn to Ser, was found only in the black-boned sheep. Tissue expression analysis indicated that this gene was generally expressed in most tissues and differently expressed in tissues of black-boned sheep. This the first report of the ovine PSAM6 gene.

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.

Reporting 678 putative cSNPs from full-length enriched cDNA sequences of the Korean native pig

Sequences from the clones of full-length enriched cDNA libraries serve as valuable resources for functional genomic studies. We have analysed 1970 high-quality chromatograms (Phred value >or= 30) that were obtained from sequencing the 5' ends of brainstem, liver, neocortex and spleen clones derived from full-length enriched cDNA libraries from Korean native pigs. In addition, 50,000 pig expressed sequence tag (EST) sequence trace files were obtained from Genbank and combined with our sequencing information to facilitate SNP identification in silico. The process generated 8118 contigs, of which 239 included minimum one sequence from Korean native pig and contained 678 putative coding single nucleotide polymorphisms (cSNPs). Of these, 33 putative cSNPs were randomly selected for confirmatory analysis and validated using 20 pigs from four different breeds (Duroc, Landrace, Yorkshire, Korean native pig). Of the 33 putative cSNPs, 20 were confirmed (61%), which was similar to the frequency reported in other studies. We also identified 15 new cSNPs from the validation process, which were not detected by our in silico analysis. Our study shows that analysing genetically diverse pig breeds including the Korean native pig could serve as a useful strategy for generating a large number of cSNPs.

An integrated RH map of porcine chromosome 10

BACKGROUND

Whole genome radiation hybrid (WG-RH) maps serve as "scaffolds" to significantly improve the orientation of small bacterial artificial chromosome (BAC) contigs, order genes within the contigs and assist assembly of a sequence-ready map for virtually any species. Here, we report the construction of a porcine: human comparative map for pig (Sus scrofa) chromosome 10 (SSC10) using the IMNpRH2(12,000-rad) porcine WG-RH panel, integrated with the IMpRH(7000-rad) WG-RH, genetic and BAC fingerprinted contig (FPC) maps.

RESULTS

Map vectors from the IMNpRH2(12,000-rad) and IMpRH(7,000-rad) panels were merged to construct parallel framework (FW) maps, within which FW markers common to both panels have an identical order. This strategy reduced map discrepancies between the two panels and significantly improved map accuracy. A total of 216 markers, including 50 microsatellites (MSs), 97 genes and ESTs, and 69 BAC end sequences (BESs), were ordered within two linkage groups at two point (2 pt) LOD score of 8. One linkage group covers SSC10p with accumulated map distances of 738.2 cR(7,000) and 1814.5 cR(12,000), respectively. The second group covers SSC10q at map distances of 1336.9 cR(7,000) and 3353.6 cR(12,000), yielding an overall average map resolution of 16.4 kb/cR(12,000) or 393.5 kb per marker on SSC10. This represents an approximately 2.5-fold increase in map resolution over the IMpRH(7,000-rad) panel. Based on 127 porcine markers that have homologous sequences in the human genome, a detailed comparative map between SSC10 and human (Homo sapiens) chromosome (HSA) 1, 9 and 10 was built.

CONCLUSION

This initial comparative RH map of SSC10 refines the syntenic regions between SSC10 and HSA1, 9 and 10. It integrates the IMNpRH2(12,000-rad) and IMpRH(7,000-rad), genetic and BAC FPC maps and provides a scaffold to close potential gaps between contigs prior to genome sequencing and assembly. This map is also useful in fine mapping of QTLs on SSC10.

Gene expression profiles of rat olfactory bulb at developmental stage

Microarray analysis may be a useful tool to identify some candidate genes related to the development of olfactory bulbs. In the present study, gene expression profiles of olfactory bulbs from postnatal day 1 (P1) rats and postnatal day 35 (P35) rats were analyzed by oligonucleotide-microarray and expression levels of some selected genes were also confirmed by RT-PCR and in situ hybridization. 9146 genes were commonly identified in six microarray chips. Among these genes, 76 were up-regulated and 130 were down-regulated three-folds or more at P1 olfactory bulbs. Out of these 76 up-regulated genes, 24 genes were annotated based on the NCBI database of reference sequences and expression levels of these 24 genes were confirmed by RT-PCR. Among them, 2 interesting genes (neurogenic differentiation 1 and retinoid acid receptor alpha) were localized in the P1 olfactory bulb by the use of in situ hybridization technique. Our results may provide basic information to identify genes associated with functional growth of olfactory bulbs.

High-resolution comprehensive radiation hybrid maps of the porcine chromosomes 2p and 9p compared with the human chromosome 11

We are constructing high-resolution, chromosomal 'test' maps for the entire pig genome using a 12,000-rad WG-RH panel (IMNpRH2(12,000-rad))to provide a scaffold for the rapid assembly of the porcine genome sequence. Here we present an initial, comparative map of human chromosome (HSA) 11 with pig chromosomes (SSC) 2p and 9p. Two sets of RH mapping vectors were used to construct the RH framework (FW) maps for SSC2p and SSC9p. One set of 590 markers, including 131 microsatellites (MSs), 364 genes/ESTs, and 95 BAC end sequences (BESs) was typed on the IMNpRH2(12,000-rad) panel. A second set of 271 markers (28 MSs, 138 genes/ESTs, and 105 BESs) was typed on the IMpRH(7,000-rad) panel. The two data sets were merged into a single data-set of 655 markers of which 206 markers were typed on both panels. Two large linkage groups of 72 and 194 markers were assigned to SSC2p, and two linkage groups of 84 and 168 markers to SSC9p at a two-point LOD score of 10. A total of 126 and 114 FW markers were ordered with a likelihood ratio of 1000:1 to the SSC2p and SSC9p RH(12,000-rad) FW maps, respectively, with an accumulated map distance of 4046.5 cR(12,000 )and 1355.2 cR(7,000 )for SSC2p, and 4244.1 cR(12,000) and 1802.5 cR(7,000) for SSC9p. The kb/cR ratio in the IMNpRH2(12,000-rad) FW maps was 15.8 for SSC2p, and 15.4 for SSC9p, while the ratio in the IMpRH(7,000-rad) FW maps was 47.1 and 36.3, respectively, or an approximately 3.0-fold increase in map resolution in the IMNpRH(12,000-rad) panel over the IMpRH(7,000-rad) panel. The integrated IMNpRH(12,000-rad) andIMpRH(7,000-rad) maps as well as the genetic and BAC FPC maps provide an inclusive comparative map between SSC2p, SSC9p and HSA11 to close potential gaps between contigs prior to sequencing, and to identify regions where potential problems may arise in sequence assembly.

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.

Assignment of 117 genes from HSA5 to the porcine IMpRH map and generation of a dense human?pig comparative map

Twenty-two and eight significant quantitative trait loci for economically important traits have been located on porcine chromosomes (SSC) 2q and SSC16 respectively, both of which have been shown to correspond to human chromosome 5 (HSA5) by chromosome painting. To provide a comprehensive comparative map for efficient selection of candidate genes, we assigned 117 genes from HSA5 using a porcine radiation hybrid (IMpRH) panel. Sixty-six genes were assigned to SSC2 and 48 to SSC16. One gene was suggested to link to SSC2 markers and another to SSC6. One gene did not link to any gene, expressed sequence tag or marker in the map, including those in the present investigation. This study demonstrated the following: (1) SSC2q21-q28 corresponds to the region ranging from 74.0 to 148.2 Mb on HSA5q13-q32 and the region from 176.0 to 179.3 Mb on HSA5q35; (2) SSC16 corresponds to the region from 1.4 to 68.7 Mb on HSA5p-q13 and to the region from 150.4 to 169.1 Mb on HSA5q32-q35 and (3) the conserved synteny between HSA5 and SSC2q21-q28 is interrupted by at least two sites and the synteny between HSA5 and SSC16 is also interrupted by at least two sites.

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.

Assignment of 101 genes localized in HSA10 to a swine RH (IMpRH) map to generate a dense human-swine comparative map

Economically important traits such as growth and backfat in pigs have been shown to be influenced by genes in swine chromosome (SSC) 10q12-->qter corresponding to human chromosome (HSA) 10p. However, since gene information in the swine chromosomal region was limited, we attempted to generate a dense comparative map between SSC10 and HSA10 by mapping the 115 genes of HSA10 to a swine RH map (IMpRH map). In the mapping ten genes were assigned to SSC10, 88 to SSC14, and one to SSC3. One gene was suggested to link to SSC3, and another to SSC9. The correspondences between HSA10 and SSC10 and between HSA10 and SSC14 were essentially consistent with the observations obtained from bi/uni-directional chromosome painting or other results. This study further indicated that a large number of intrachromosomal rearrangements occurred in the synteny-conserved regions following species separation.

Assignment of 204 genes localized on HSA17 to a porcine RH (IMpRH) map to generate a dense comparative map between pig and human/mouse

Bi- and uni-directional chromosome painting (ZOO-FISH) and gene mapping have revealed correspondences between human chromosome (HSA) 17 and porcine chromosome (SSC) 12 harboring economically important quantitative trait loci. In the present study, we have assigned 204 genes localized on HSA17 to SSC12 to generate a comprehensive comparative map between HSA17 and SSC12. Two hundred fifty-five primer pairs were designed using porcine sequences orthologous with human genes. Of the 255 primer pairs, 208 (81.6%) were used to assign the corresponding genes to porcine chromosomes using the INRA-Minnesota 7000-rad porcine x Chinese hamster whole genome radiation hybrid (IMpRH) panel. Two hundred three genes were integrated into the SSC12 IMpRH linkage maps; and one gene, PPARBP, was found to link to THRA1 located in SSC12 but not incorporated into the linkage maps. Three genes (GIT1, SLC25A11, and HT008) were suggested to link to SSC12 markers, and the remaining gene (RPL26) did not link to any genes/expressed sequence tags/markers registered, including those in the present study. A comparison of the gene orders among SSC12, HSA17, and mouse chromosome 11 indicates that intra-chromosomal rearrangements occurred frequently in this ancestral mammalian chromosome during speciation.

A 12,000-rad porcine radiation hybrid (IMNpRH2) panel refines the conserved synteny between SSC12 and HSA17

Reverse or bidirectional Zoo-FISH suggests that synteny between porcine chromosome 12 (SSC12) and human chromosome 17 (HSA17) is completely conserved. The construction of a high-resolution radiation hybrid (RH) map for SSC12 provides a unique opportunity to determine whether chromosomal synteny is reflected at the molecular level by comparative gene mapping of SSC12 and HSA17. We report an initial, high-resolution RH map of SSC12 on the 12,000-rad IMNpRH2 panel using CarthaGene software. This map contains a total of 320 markers, including 20 microsatellites and 300 ESTs/genes, covering approximately 4836.9 cR12,000. The markers were ordered in 16 linkage groups at LOD 6.0 using framework markers previously mapped on the IMpRH7000-rad SSC12 and porcine genetic maps. Ten linkage groups ordered more than 10 markers, with the largest containing 101 STSs. The resolution of the current RH map is approximately 15.3 kb/cR on SSC12, a significant improvement over the second-generation EST SSC12 RH7000-rad map of 103 ESTs and 15 framework markers covering approximately 2287.2 cR7000. Compared to HSA17, six distinct segments were identified, revealing macro-rearrangements within the apparently complete synteny between SSC12 and HSA17. Further analysis of the order of 245 genes (ESTs) on HSA17 and SSC12 also revealed several micro-rearrangements within a synteny segment. A high-resolution SSC12 RH12,000-rad map will be useful in fine-mapping QTL and as a scaffold for sequencing this chromosome.

A dense comparative gene map between human chromosome 19q13.3-->q13.4 and a homologous segment of swine chromosome 6

The human chromosome (HSA)19q region has been shown to correspond to swine chromosome (SSC) 6q11-->q21 by bi-directional chromosomal painting and gene mapping. However, since the precise correspondence has not been determined, 26 genes localized in HSA19q13.3-->q13.4 were assigned to the SSC6 region mainly by radiation hybrid (RH) mapping, and additionally, by somatic cell hybrid panel (SCHP) mapping, and fluorescent in situ hybridization (FISH). Out of the 26 genes, 24 were assigned to a swine RH map with LOD scores greater than 6 (threshold of significance). The most likely order of the 24 genes along SSC6 was calculated by CarthaGene, revealing that the order is essentially the same as that in HSA19q13.3-->q13.4. For AURKC and RPS5 giving LOD scores not greater than 6, SCHP mapping and FISH were additionally performed; SCHP mapping assigned AURKC and RPS5 to SSC6q22-->q23 and SSC6q21, respectively, which is consistent with the observation of FISH. Consequently, all the genes (26 genes) examined in the present study were shown to localize in SSC6q12-->q23, and the order of the genes along the chromosomes was shown to be essentially the same in swine and human, though several intrachromosomal rearrangements were observed between the species.