Construction of a bacterial artificial chromosome library for the Rongchang pig breed and its use for the identification of genes involved in intramuscular fat deposition

In a search for genes affecting intramuscular fat deposition, we constructed a bacterial artificial chromosome (BAC) library for the whole genome of Rongchang pig, a domestic Chinese swine breed. The library consisted of approximately 192,000 clones, with an averaged insert size of 116 kb. Frequency of non-insert clone of the BAC library was no higher than 1.8%, based on estimation of 220 BAC clones randomly selected. We estimated the coverage of the library to be more than seven porcine genome equivalents. Subsequent screening of the BAC library with a three-step PCR procedure resulted in identification of seven candidate genes that were potentially involved in intramuscular fat deposition. The number of positive BAC clones ranged from 2 to 4 for each of the seven genes. One positive clone, containing the lipin1 gene, was fully sequenced by shotgun method to generate 118,041 bp porcine genomic sequences. The BAC clone contained complete DNA sequence of porcine lipin1 gene including all the exons and introns. Our results indicate that this BAC library is a useful tool for gene identification and help to serve as an important resource for future porcine genomic study.     

Comparative human-mouse-rat sequence analysis of the ICAM gene cluster on HSA 19p13.2 and a 185-kb porcine region from SSC 2q

The human intercellular adhesion molecule gene (ICAM) cluster is located in a GC-rich and gene-rich region on HSA 19p13.2. We determined the complete DNA sequence of a 185-kb porcine bacterial artificial chromosome (BAC) clone containing parts of the ICAM gene cluster. We used the porcine sequence for a detailed comparative analysis between human, pig, mouse and rat. The 185 kb of porcine sequence covered 220 kb of homologous sequence in the human genome, which adds to the growing evidence that the porcine genome is somewhat smaller than the human genome. The genomic sequences of the four species showed a high level of conserved synteny and no rearrangements in gene order were observed. During evolution, the ICAM3 gene was inactivated by mutation in the mouse and rat genome, whereas it is still present in the human and pig genome. The loss of Icam3 in rodent genomes might be relevant for rodent-specific properties of the T-cell-mediated immune response. All the other investigated genes are conserved across all four investigated sequences.     

Improving the comparative map of SSC2p-q13 by sample sequencing of BAC clones

To improve the comparative map for pig chromosome 2 and increase the gene density on this chromosome, a porcine bacterial artificial chromosome (BAC) library was screened with 17 microsatellite markers and 18 genes previously assigned to pig chromosome 2. Fifty-one BAC clones located in the region of a maternally imprinted quantitative trait locus for backfat thickness (BFT) were identified. From these BACs 372 kb were sample sequenced. The average read length of a subclone was 442 basepair (bp). Contig assembly analysis showed that every bp was sequenced 1.28 times. Subsequently, sequences were compared with sequences in the nucleotide databases to identify homology with other mammalian sequences. Sequence identity was observed with sequences derived from 35 BACs. The average percentage identity with human sequences was 87.6%, with an average length of 143 bp. In total, sample sequencing of all BACs resulted in sequence identity with 29 human genes, 13 human expressed sequence tags (ESTs), 17 human genomic clones, one rat gene, one porcine gene and nine porcine ESTs. Eighteen genes located on human chromosome 11 and 19, and seven genes from other human locations, one rat gene and one porcine gene were assigned to pig chromosome 2 for the first time. The new genes were added to the radiation hybrid map at the same position as the locus from which the BAC that was sequenced was derived. In total 57 genes were placed on the radiation hybrid map of SSC2p-q13.     

Multiple isoforms of porcine aromatase are encoded by three distinct genes

Cytochrome P450 aromatase, a product of the CYP 19 gene and the terminal enzyme in the estrogen biosynthetic pathway, is synthesized by the ovary, endometrium, placenta, and peri-implantation embryos in the pig and other mammals, albeit to varying levels, implying its functional role(s) in pregnancy events. The aromatase produced by the pig tissues exists as three distinct isoforms (type I - ovary, type II - placenta, and type III - embryo), with presumed differences in substrate specificities, expression levels, activity, and mode of regulation. In order to delineate the molecular mechanisms whereby estrogen synthesis is regulated in these diverse tissues, the present study examined if these aromatase isoforms represent products of multiple genes or of a single gene via complex splicing mechanisms. Porcine genomic DNA from a single animal was used as a template in the polymerase chain reaction (PCR) to amplify isoform-specific sequences corresponding to exons 4 and 7, respectively. Nucleotide sequence analysis of the generated fragments revealed the presence of only clones corresponding to the three known aromatase types. Screening a porcine Bacterial Artificial Chromosome (BAC) library for aromatase gene by PCR yielded a single clone approximately 80 kb in length. Southern blot analysis, using probes specific for exons 1A-1B, 2-3, 4-9, and 10 sequences indicated that the BAC genomic clone contains the entirety of the coding exons as well as the proximal promoter region. Sequence analysis of the fragment generated with exon 4 primers determined that this BAC clone contains only the type II gene. The presence and relative orientation of the untranslated 5'- exons 1A and 1B, previously demonstrated for the type III isoform were evaluated in the BAC clone and genomic DNA by PCR. The 265 bp fragment generated from both PCR reactions was confirmed by sequence analysis to contain exons 1A and 1B that are located contiguous to each other and separated by only three bp. A diagnostic procedure for typing aromatase isoforms was developed, based on the presence of specific restriction sites within isoform-specific exons. The use of this protocol confirmed the existence of only three aromatase isoforms in the porcine genome and indicated changes in aromatase types expressed by the uterine endometrium as a function of pregnancy stage. The presence of distinct genes encoding each of the aromatase isoform predicts important differences in the mechanisms underlying the molecular evolution and regulation of porcine aromatase, unique from those of other mammals, and suggests a critical role for P450 aromatase steroidal products in uterine functions related to pregnancy events.     

Genetic and physical mapping, expression analysis and partial sequence of porcine PER1

The porcine PER1 gene was mapped to chromosome 12q1.4-->q1.5 using fluorescence in situ hybridisation. A polymorphic microsatellite marker (S0601) was isolated from a BAC clone shown to contain the PER1 gene. Linkage analysis assigned S0601 distal to ALOX12 on SSC12, providing further evidence for the conservation of synteny between HSA17 and SSC12. RT-PCR analysis demonstrated the expression of PER1 in all 11 tissues tested, consistent with the data from other mammalian species. Part of the PER1 gene was sequenced, homologous to exons 2-14 of the human gene and encoding the N-terminus of porcine PER1. The predicted amino acid sequence of the partial pig PER1 protein shares over 96% identity with its human orthologue.     

Analysis of the genomic structure of the porcine CD1 gene cluster

CD1 is an MHC class I-like protein that presents lipid antigens to T cell receptors. We determined 470,187 bp of the genomic sequence encompassing the region encoding porcine CD1 genes. We identified 16 genes in this region and newly identified CD1A2, CD1B, CD1C, CD1D, and CD1E. Porcine CD1 genes were located in clusters between KIRREL and olfactory receptor (OR) genes, as observed in humans, although they were divided into two regions by a region encoding OR genes. Comparison of the genomic sequences of CD1 gene loci in pigs with other mammals showed that separation of the CD1 gene cluster by ORs was observed only in pigs. CD1A duplication in the porcine genome was estimated to have occurred after the divergence of the human and porcine. This analysis of the genomic sequence of the porcine CD1 family will contribute to our understanding of the evolution of mammalian CD1 genes.     

Exploiting dinucleotide microsatellites conserved among mammalian species

Dinucleotide microsatellites are useful for gene mapping projects. Depending upon definition of conservation, published estimates of dinucleotide microsatellite conservation levels vary dramatically (30% to 100%). This study focused on well-characterized genes that contain microsatellites in the human genome. The objective was to examine the feasibility of developing microsatellite markers within genes on the basis of the assumption of microsatellite conservation across distantly related species. Eight genes (Gamma-actin, carcinoembryonic antigen, apolipoprotein A-II, cardiac beta myosin heavy chain, laminin B2 chain, MHC class I CD8 alpha chain, c-reactive protein, and retinoblastoma susceptibility protein) containing large dinucleotide repeat units (N ≥ 15), complete genomic structure information, and homologous gene sequences in a second species were selected. Heterologous primers were designed from conserved exon sequences flanking a microsatellite motif. PCR products from bovine and porcine genomic DNA were tested for the presence of microsatellite sequences by Southern blot hybridization with biotin-labeled (CA)₁₂ oligonucleotides. Fragments containing microsatellites were cloned and sequenced. Homology was verified by sequence comparisons between human and corresponding bovine or porcine fragments. Four of sixteen (25%) cross-amplified PCR products contained dinucleotide repetitive sequences with repeat unit lengths of 5 to 23. Two dinucleotide repetitive sequences showed microsatellite length polymorphism, and an additional sequence displayed single-strand conformational polymorphism. Results from this study suggest that exploitation of conserved microsatellite sequences is a useful approach for developing specific genetic markers for comparative mapping purposes. Received: 7 July 1995 / Accepted: 28 September 1995     

Genomic structure and expression of the soluble guanylyl cyclase alpha2 subunit gene in the medaka fish Oryzias latipes

A cDNA clone encoding the soluble guanylyl cyclase alpha2 subunit was isolated from medaka fish (Oryzias latipes) and designated as OlGCS-alpha2. The OlGCS-alpha2 cDNA was 3,192 bp in length and the open reading frame (ORF) encodes a protein of 805 amino acids. The deduced amino acid sequence has high similarity to that of the mammalian alpha2 subunit gene except for the N-terminal regulatory domain. The C-terminal 5 amino acids, "RETSL", which have been reported to interact with the post synaptic density protein (PSD)-95 were conserved. An RNase protection assay with adult fish organs showed that OlGCS-alpha2 was expressed mainly in the brain and testis. The complete nucleotide sequence (about 41 kbp) of the OlGCS-alpha2 genomic DNA clone isolated from a medaka fish BAC library indicated that the OlGCS-alpha2 gene consisted of 9 exons and 8 introns. The 5'-flanking region and larger introns, such as introns 1, 4, and 7, contained the several fragments conserved in the nucleotide sequences of Rex6 (non-long terminal repeat retrotransposon), MHC class I genomic region, and OlGC1, the medaka fish homolog of the mammalian guanylyl cyclase B gene. Linkage analysis on the medaka fish chromosome demonstrated that the OlGCS-alpha2 gene was mapped to LG13; this mapping position was different from those for the OlGCS-alpha1 and OlGCS-beta1 genes (LG1).     

Targeted analysis of orthologous phytochrome A regions of the sorghum,maize, and rice genomes using comparative gene-island sequencing

A "gene-island" sequencing strategy has been developed that expedites the targeted acquisition of orthologous gene sequences from related species for comparative genome analysis. A 152-kb bacterial artificial chromosome (BAC) clone from sorghum (Sorghum bicolor) encoding phytochrome A (PHYA) was fully sequenced, revealing 16 open reading frames with a gene density similar to many regions of the rice (Oryza sativa) genome. The sequences of genes in the orthologous region of the maize (Zea mays) and rice genomes were obtained using the gene-island sequencing method. BAC clones containing the orthologous maize and rice PHYA genes were identified, sheared, subcloned, and probed with the sorghum PHYA-containing BAC DNA. Sequence analysis revealed that approximately 75% of the cross-hybridizing subclones contained sequences orthologous to those within the sorghum PHYA BAC and less than 25% contained repetitive and/or BAC vector DNA sequences. The complete sequence of four genes, including up to 1 kb of their promoter regions, was identified in the maize PHYA BAC. Nine orthologous gene sequences were identified in the rice PHYA BAC. Sequence comparison of the orthologous sorghum and maize genes aided in the identification of exons and conserved regulatory sequences flanking each open reading frame. Within genomic regions where micro-colinearity of genes is absolutely conserved, gene-island sequencing is a particularly useful tool for comparative analysis of genomes between related species.     

Novel gene exon homologous to pancreatic phospholipase A2: sequence and chromosomal mapping of both human genes

We described previously the cloning and DNA sequence of the human gene encoding pancreatic phospholipase A2 [DNA 5, 519]. When pancreatic phospholipase A2 (PLA2) cDNA was used to screen a human genomic library, two classes of clones were obtained. One class encoded the pancreatic enzyme, and a second class encoded one exon of an apparently related PLA2. No additional PLA2 gene exons displayed sufficient homology to be detected by the probe. A homologous sequence in both rat and porcine genomic DNA was detected by DNA blot hybridization, and the corresponding gene fragments were cloned and sequenced. Within the deduced amino acid sequences, the presence of known functional residues along with the high degree of interspecies conservation suggests the genes encode a functional PLA2 enzyme form. The encoded sequence lacks Cys11, as do the "type II" viperid venom and other nonpancreatic mammalian PLA2 enzymes. The sequence is distinct from porcine intestinal PLA2 and appears not to be a direct homolog of the recently published rabbit ascites and rat platelet enzymes. Hybridization of DNA probes containing sequences from these genes to genomic DNA blots of mouse/human somatic cell hybrids permitted chromosomal assignment for both. The pancreatic gene mapped to human chromosome 12, and the homologous gene mapped to chromosome 1.     

水稻中一个NBS-LRR抗病同源基因家族的克隆和分析

利用克隆的抗病基因同源序列RS13作为探针,从水稻IR64的BAC文库中筛选到4个阳性克隆,其中一个克隆14E19能够覆盖其余3个克隆。对14E19进行全序列测定和分析,获得了73kb的全长DNA序列,基因预测显示其上有4个编码NBS-LRR结构域的基因(NL),分别命名为NL-A,B,C和D。对具有相同基因组背景的IRBB56同一染色体位置上跨度更大的BAC克隆106P13进行分析,发现其上有10个NL同源拷贝,其中4个同14E19上的NL一样。搜索日本晴、93—11、广陆矮4号的序列,发现三者有类似的同源序列。但与已知的NBS-LRR抗病基因同源性较低,说明NL是一个至少由10个成员(分别命名为NL-A至J)组成的新基因家族。对NL家族进行RT-PCR和cDNA库筛选分析,发现NL-B基因能够在抗白叶枯病品系IRBB4中表达,暗示该基因参与了抗病反应。     

Cloning, transcription and chromosomal localization of the porcine whey acidic protein gene and its expression in HC11 cell line

The whey acidic protein (WAP) is the major whey protein of rodent, rabbit and camel. Recently, it was identified in the milk of swine (Simpson et al., 1998. J. Mol. Endocrinol. 20, 27-35). In this paper, the cloning of the pig WAP cDNA and of bacterial artificial chromosome (BAC) construct containing the entire porcine WAP gene is reported. The comparison of the coding sequence of the pig WAP gene to rodent or lagomorph WAP sequence already published demonstrated that only exon sequences are partially conserved. The porcine WAP gene was localized on the subtelomeric region of the chromosome 18. The estimation of the expression of the swine WAP gene in the mammary gland from lactating animals revealed a high level of expression. In order to compare the expression level of the porcine WAP gene from the large genomic fragment which contained 70 kb downstream and 50 kb upstream the pig WAP gene or the smaller one (1 kb downstream and 2.4 kb upstream), these two genomic fragments were transfected in HC11 cell line. The BAC construct was expressed 15 times higher than the plasmid when reported to the integrated copy number. This report suggests that the HC11 cell line is a useful tool to identify the regulatory sequences of milk protein genes.     

Gene identification in a complex chromosomal continuum by local genomic cross-referencing

Most higher plants have complex genomes containing large quantities of repetitive DNA interspersed with low-copy-number sequences. Many of these repetitive DNAs are mobile and have homology to RNAs in various cell types. This can make it difficult to identify the genes in a long chromosomal continuum. It was decided to use genic sequence conservation and grass genome co-linearity as tools for gene identification. A bacterial artificial chromosome (BAC) clone containing sorghum genomic DNA was selected using a maize Adh1 probe. The 165 kb sorghum BAC was tested for hybridization to a set of clones representing the contiguous 280 kb of DNA flanking maize Adh1. None of the repetitive maize DNAs hybridized, but most of the low-copy-number sequences did. A low-copy-number sequence that did cross-hybridize was found to be a gene, while one that did not was found to be a low-copy-number retrotransposon that was named Reina. Regions of cross-hybridization were co-linear between the two genomes, but closer together in the smaller sorghum genome. These results indicate that local genomic cross-referencing by hybridization of orthologous clones can be an efficient and rapid technique for gene identification and studies of genome organization.     

Comparative analysis of a Brassica BAC clone containing several major aliphatic glucosinolate genes with its corresponding Arabidopsis sequence.

We compared the sequence of a 101-kb-long bacterial artificial chromosome (BAC) clone (B21H13) from Brassica oleracea with its homologous region in Arabidopsis thaliana. This clone contains a gene family involved in the synthesis of aliphatic glucosinolates. The A. thaliana homologs for this gene family are located on chromosome IV and correspond to three 2-oxoglutarate-dependent dioxygenase (AOP) genes. We found that B21H13 harbors 23 genes, whereas the equivalent region in Arabidopsis contains 37 genes. All 23 common genes have the same order and orientation in both Brassica and Arabidopsis. The 16 missing genes in the broccoli BAC clone were arranged in two major blocks of 5 and 7 contiguous genes, two singletons, and a twosome. The 118 exons comprising these 23 genes have high conservation between the two species. The arrangement of the AOP gene family in A. thaliana is as follows: AOP3 (GS-OHP) - AOP2 (GS-ALK) - pseudogene - AOP1. In contrast, in B. oleracea (broccoli and collard), two of the genes are duplicated and the third, AOP3, is missing. The remaining genes are arranged as follows: Bo-AOP2.1 (BoGSL-ALKa) - pseudogene - AOP2.2 (BoGSL-ALKb) - AOP1.1 - AOP1.2. When the survey was expanded to other Brassica accessions, we found variation in copy number and sequence for the Brassica AOP2 homologs. This study confirms that extensive rearrangements have taken place during the evolution of the Brassicacea at both gene and chromosomal levels.     

Molecular cloning and expression of the CRISP family of proteins in the boar

The family of mammalian cysteine-rich secretory proteins (CRISP) have been well characterized in the rat, mouse, and human. Here we report the molecular cloning and expression analysis of CRISP1, CRISP2, and CRISP3 in the boar. A partial sequence published in the National Center for Biotechnology Information (NCBI) database was used to derive the full-length sequences for CRISP1 and CRISP2 using rapid amplification of cDNA ends. RT-PCR confirmed the expression of these mRNAs in the boar reproductive tract, and real time RT-PCR showed CRISP1 to be highly expressed throughout the epididymis, with CRISP2 highly expressed in the testis. A search of the porcine genomic sequence in the NCBI database identified a BAC (CH242-199E6) encoding the CRISP1 gene. This BAC is derived from porcine Chromosome 7 and is syntenic with the regions of the mouse, rat, and human genomes encoding the CRISP gene family. This BAC was found to encode a third CRISP protein with a predicted amino acid sequence of high similarity to human CRISP3. Using RT-PCR we show that CRISP3 expression in the boar reproductive tract is confined to the prostate. Recombinant porcine (rp) CRISP2 protein was produced and purified. When incubated with capacitated boar sperm, rpCRISP2 induced an acrosome reaction, consistent with its demonstrated ability to alter the activity of calcium channels.     

Complete amino acid sequence of human cartilage link protein (CRTL1) deduced from cDNA clones and chromosomal assignment of the gene

Little is known about the primary amino acid structure of human cartilage link protein (CRTL1). We screened a human genomic library with a cDNA encoding the 3' untranslated region and the adjoining B1 domain of chicken link protein. One clone was isolated and characterized. A 3.5-kb EcoRI-KpnI fragment from this genomic clone that contains the human B1 exon was used to map the gene to chromosome 5q13----q14.1. The same fragment was used to screen a cDNA library prepared from mRNA of Caco-2, a human colon tumor cell line. Two overlapping clones were isolated and shown to encode all of CRTL1. The deduced amino acid sequence is 354 residues long. The amino acid sequence shows a striking degree of identity to the porcine (96%), rat (96%), and chicken (85%) link protein sequences. Furthermore, there is greater than 86% homology between the 3' untranslated region of the genes encoding human and porcine link proteins. These results indicate that there has been strong evolutionary pressure against changes in the coding and 3' untranslated regions of the gene encoding cartilage link protein.