JESAM: CORBA software components to create and publish EST alignments and clusters

MOTIVATION: Expressed Sequence Tags (ESTs) are cheap, easy and quick to obtain relative to full genomic sequencing and currently sample more eukaryotic genes than any other data source. They are particularly useful for developing Sequence Tag Sites (STSs for mapping), polymorphism discovery, disease gene hunting, mass spectrometer proteomics, and most ironically for finding genes and predicting gene structure after the great effort of genomic sequencing. However, ESTs have many problems and the public EST databases contain all the errors and high redundancy intrinsic to the submitted data so it is often found that derived database views, which reduce both errors and redundancy, are more effective starting points for research than the original raw submissions. Existing derived views such as EST cluster databases and consensus databases have never published supporting evidence or intermediary results leading to difficulties trusting, correcting, and customizing the final published database. These difficulties have lead many groups to wastefully repeat the complex intermediary work of others in order to offer slightly different final views. A better approach might be to discover the most expensive common calculations used by all the approaches and then publish all intermediary results. Given a globally accessible database with a suitable component interface, like the JESAM software described in this paper, the creation of customized EST-derived databases could be achieved with minimum effort. RESULTS: Databases of EST and full-length mRNA sequences for four model organisms have been self-compared by searching for overlaps consistent with contiguity. The sequence comparisons are performed in parallel using a PVM process farm and previous results are stored to allow incremental updates with minimal effort. The overlap databases have been published with CORBA interfaces to enable flexible global access as demonstrated by example Java applet browsers. Simple cDNA supercluster databases built as alignment database clients are themselves published via CORBA interfaces browsable with prototypical applets. A comparison with UniGene Mouse and Rat databases revealed undesirable features in both and the advantages of contrasting perspectives on complex data. AVAILABILITY: The software is packaged as two Jar files available from: URL: http://corba.ebi.ac.uk/EST/jesam/jesam. html. One jar contains all the Java source code, and the other contains all the C, C++ and IDL code. Links to working examples of the alignment and cluster viewers (if remote firewall permits) can be found at http://corba.ebi.ac.uk/EST. All the Washington University mouse EST traces are available for browsing at the same URL.     

Cloning and characterization of a novel Krüppel-like zinc finger gene, ZNF268, expressed in early human embryo

With the aim of identifying genes involved in early human embryonic development, we have isolated a cDNA clone representing a novel human zinc finger gene ZNF268 from 3 week old human embryo cDNA library using a differential hybridization strategy. The complete cDNA sequence of ZNF268 contained an open reading frame of 2841 nucleotides that encodes a 947 amino acid protein with an N-terminal Krüppel-associated box (KRAB) domain and 24 C-terminal zinc fingers. Northern blot analysis showed that ZNF268 mRNA is mainly expressed in 3-5 week old human embryos suggesting it could play certain roles in the embryogenesis. The gene consists of six exons spanning about 22 kb of genomic DNA. According to the genomic sequence from the HTGS database, the ZNF268 gene is assigned to human chromosome 5.     

Physical and transcriptional map of a 3-Mb region of mouse chromosome 1 containing the gene for the neural tube defect mutant loop-tail (Lp)

The Lp mouse mutant provides a model for the severe human neural tube defect (NTD), cranio-rachischisis. To identify the Lp gene, a positional cloning approach has been adopted. Previously, linkage analysis in a large intraspecific backcross was used to map the Lp locus to distal mouse chromosome 1. Here we report a detailed physical map of this region. The interval surrounding Lp has been cloned in a yeast artificial chromosome (YAC) contig consisting of 63 clones spanning approximately 3.2 Mb. Fifty sequence tagged sites (STSs) have been used to construct the contig and establish marker order across the interval. Based on the high level of conserved synteny between distal mouse chromosome 1 and human 1q21-q24, many of these STSs were designed from expressed sequences identified by cross-screening human and mouse databases of expressed sequence tags. Added to other known genes in the region, a total of 29 genes were located and ordered within the contig. Seven novel polymorphisms were identified within the region, allowing refinement of the genetic map and a reduction in the size of the physical interval containing the Lp gene. The Lp interval, between D1Mit113 and Tagln2, can be spanned by two nonchimeric overlapping YACs that define a physical distance of approximately 1 Mb. Within this region, 10 potential candidate genes have been mapped. The materials and genes described here will provide a resource for the identification and further study of the mutated Lp gene that causes this severe neural tube defect and will provide candidates for other defects known to map to the homologous region on human chromosome 1q.     

C lambda 2 and C lambda 4 immunoglobulin light chain genes in a wild-derived inbred mouse strain

A genomic clone containing the lambda constant (C) region genes (C lambda 2S and C lambda 4S) has been isolated from a genomic library from the mouse strain SPE. SPE is an inbred strain derived from progenitors trapped near Grenada in Spain and has been classified as mus 3 or mus spretus. The sequence of the C lambda 2S gene is virtually identical to that of BALB/c both in the coding region and in flanking sequences, suggesting that it is an expressed gene in the SPE strain. By contrast, the C lambda 4S gene on the same cluster has diverged in sequence from that of BALB/c and contains a large deletion that precludes its normal expression. Whereas BALB/c J lambda 4 region contains substitutions that probably preclude its usage, the SPE J lambda 4 gene includes all sequences required for a functional J gene. Comparison of the C lambda 2S and C lambda 4S gene sequences with those available for BALB/c C lambda 3 and C lambda 1 confirms the close relationship between the C lambda 1-C lambda 4 and C lambda 2-C lambda 3 gene pairs. The C lambda 3 gene of BALB/c is more closely related to C lambda 2S than is C lambda 1 of BALB/c to C lambda 4S. If it is assumed that C lambda 1 and C lambda 2 are respective duplicates of C lambda 4 and C lambda 3 and that these duplications occurred at the same time, then the C lambda 2 gene has been under stronger selective pressure than C lambda 4.     

Evidence for bacterial origin of heat shock RNA-1

The heat shock RNA-1 (HSR1) is a noncoding RNA (ncRNA) reported to be involved in mammalian heat shock response. HSR1 was shown to significantly stimulate the heat-shock factor 1 (HSF1) trimerization and DNA binding. The hamster HSR1 sequence was reported to consist of 604 nucleotides (nt) plus a poly(A) tail and to have only a 4-nt difference with the human HSR1. In this study, we present highly convincing evidence for bacterial origin of the HSR1. No HSR1 sequence was found by exhaustive sequence similarity searches of the publicly available eukaryotic nucleotide sequence databases at the NCBI, including the expressed sequence tags, genome survey sequences, and high-throughput genomic sequences divisions of GenBank, as well as the Trace Archive database of whole genome shotgun sequences, and genome assemblies. Instead, a putative open reading frame (ORF) of HSR1 revealed strong similarity to the amino-terminal region of bacterial chloride channel proteins. Furthermore, the 5′ flanking region of the putative HSR1 ORF showed similarity to the 5′ upstream regions of the bacterial protein genes. We propose that the HSR1 was derived from a bacterial genome fragment either by horizontal gene transfer or by bacterial infection of the cells. The most probable source organism of the HSR1 is a species belonging to the order Burkholderiales.     

Isolation, in silico characterization and chromosomal localization of a group of cDNAs from ciliated epithelial cells after in vitro ciliogenesis

Background

Immotile cilia syndrome (ICS) or primary ciliary dyskinesia (PCD) is an autosomal recessive disorder in humans in which the beating of cilia and sperm flagella is impaired. Ciliated epithelial cell linings are present in many tissues. To understand ciliary assembly and motility, it is important to isolate those genes involved in the process.

Results

Total RNA was isolated from cultured ciliated nasal epithelial cells after in vitro ciliogenesis and expressed sequenced tags (ESTs) were generated. The functions and locations of 63 of these ESTs were derived by BLAST from two public databases. These ESTs are grouped into various classes. One group has high homology not only with the mitochondrial genome but also with one or more chromosomal DNAs, suggesting that very similar genes, or genes with very similar domains, are expressed from both mitochondrial and nuclear DNA. A second class comprises genes with complete homology with part of a known gene, suggesting that they are the same genes. A third group has partial homology with domains of known genes. A fourth group, constituting 33% of the ESTs characterized, has no significant homology with any gene or EST in the database.

Conclusions

We have shown that sufficient information about the location of ESTs could be derived electronically from the recently completed human genome sequences. This strategy of EST localization should be significantly useful for mapping and identification of new genes in the forthcoming human genome sequences with the vast number of ESTs in the dbEST database.     

A second-generation integrated map of the silkworm reveals synteny and conserved gene order between lepidopteran insects

A second-generation linkage map was constructed for the silkworm, Bombyx mori, focusing on mapping Bombyx sequences appearing in public nucleotide databases and bacterial artificial chromosome (BAC) contigs. A total of 874 BAC contigs containing 5067 clones (22% of the library) were constructed by PCR-based screening with sequence-tagged sites (STSs) derived from whole-genome shotgun (WGS) sequences. A total of 523 BAC contigs, including 342 independent genes registered in public databases and 85 expressed sequence tags (ESTs), were placed onto the linkage map. We found significant synteny and conserved gene order between B. mori and a nymphalid butterfly, Heliconius melpomene, in four linkage groups (LGs), strongly suggesting that using B. mori as a reference for comparative genomics in Lepidotera is highly feasible.     

Molecular cloning and preliminary functional analysis of two novel human KRAB zinc finger proteins, HKr18 and HKr19.

cDNA clones encoding two novel human KRAB zinc finger proteins, HKr18 and HKr19, were isolated from a human testis cDNA library. Their corresponding genes were later identified in sequences originating from chromosomes 19 and 7, respectively. On the basis of the collected information from gene and cDNA sequences, Hkr18 was found to be a protein of 94 kDa with 20 zinc finger motifs in its C terminus. The HKr19 is a smaller protein, with a molecular weight of 56 kDa containing 11 zinc finger motifs. Both HKr18 and HKr19 contained a KRAB A as well as a KRAB B domain in their N termini. Northern blot analysis showed expression of HKr18 in all human tissues tested, indicating a ubiquitous expression pattern. In contrast, HKr19 showed a more restricted tissue distribution, with detectable expression primarily in testis and fetal tissues. The HKr19 protein is a member of the large ZNF91 subfamily of KRAB zinc finger genes. A PCR-based analysis of the expression of HKr19 and other closely related genes showed that lymphoid, myeloid, and nonhematopoietic cells expressed different sets of these genes. This latter finding indicates that some members of the ZNF91 family may be involved in regulating lineage commitment during hematopoietic development. Transfection of various parts of HKr19 into human embryonic kidney cells (HEK 293 cells) showed that the entire protein and its zinc finger region were toxic to these cells when expressed at high levels. In contrast, the KRAB domain and the linker region seemed to be well tolerated.     

Duplicated KOX zinc finger gene clusters flank the centromere of human chromosome 10: evidence for a pericentric inversion during primate evolution.

Two related zinc finger (ZNF) gene clusters from the pericentromeric region of human chromosome 10, defined by cDNAs of the KOX series, have been cloned in yeast artificial chromosomes (YACs). The two clusters evolved by duplication of an ancestral gene cluster before the divergence of the human and great ape lineages. Included in cluster A are the ZNF gene sequences ZNF11A, ZNF33A, and ZNF37A, while cluster B comprises the related sequences ZNF11B, ZNF33B and ZNF37B. Genes from both clusters are expressed: cDNAs KOX2, KOX31 and KOX21 derive from ZNF11B, ZNF33A and ZNF37A, respectively. Further YACs have been isolated which link ZNF11A and ZNF33A to another gene, ZNF25, defined by cDNA clone KOX19. Therefore ZNF25 also forms part of cluster A, but has no counterpart in cluster B. Surprisingly, the KOX ZNF gene clusters are located on opposite sides of the centromere: cluster A maps to 10p11.2, while cluster B is in 10q11.2. This suggests the occurrence during primate evolution of a previously undescribed pericentric inversion subsequent to the cluster duplication. The evolution of this subset of KOX ZNF genes has therefore involved three types of genetic event: local gene duplication, gene cluster duplication, and chromosome rearrangement.     

Generation of sequence-tagged sites from Xp22.3 by isolating commonAlu-PCR products of radiation hybrids retaining overlapping human X chromosome fragments

Several human diseases have been mapped to Xp22.3 on the distal short arm of the human X chromosome, and many genes in this area have been found to be expressed from the inactive X chromosome. To facilitate physical mapping and characterization of this interesting region, we have constructed a battery of radiation hybrids containing human X chromosomal fragments, and isolated two hybrid clones with overlapping fragments of Xp22.3. Alu-PCR on these hybrids and identification of sequences common to both hybrids allowed the isolation of six sequence-tagged sites (STSs) from Xp22.3. Five of the STSs were mapped to individual YACs comprising a recently constructed contig of this region. These novel STSs are useful markers for further physical characterization of this part of the genome. Received: 4 May 1995 / Revised: 27 September 1995     

Construction of an Integrated Physical and Gene Map of Human Chromosome 20p12 Providing Candidate Genes for Alagille Syndrome

Physical mapping and localization of eSTS markers were used to generate an integrated physical and gene map covering a ∼10-Mb region of human chromosome 20p12 containing the Alagille syndrome (AGS) locus. Seventy-four STSs, 28 of which were derived from cDNA sequences, mapped with an average resolution of 135 kb. The 28 eSTS markers define 20 genes. Six known genes, namely CHGB, BMP2, PLCB1, PLCB4, SNAP, and HJ1, were precisely mapped. Among the genes identified, one maps in the smallest region of overlap of the deletions associated with AGS and could therefore be regarded as a candidate gene for Alagille syndrome.     

A comparative map of chicken chromosome 24 and human chromosome 11

To improve the physical and comparative map of chicken chromosome 24 (GGA24; former linkage group E49C20W21) bacterial artificial chromosome (BAC) contigs were constructed around loci previously mapped on this chromosome by linkage analysis. The BAC clones were used for both sample sequencing and BAC end sequencing. Sequence tagged site (STS) markers derived from the BAC end sequences were used for chromosome walking. In total 191 BAC clones were isolated, covering almost 30% of GGA24, and 76 STS were developed (65 STS derived from BAC end sequences and 11 STS derived within genes). The partial sequences of the chicken BAC clones were compared with sequences present in the EMBL/GenBank databases, and revealed matches to 19 genes, expressed sequence tags (ESTs) and genomic clones located on human chromosome 11q22-q24 and mouse chromosome 9. Furthermore, 11 chicken orthologues of human genes located on HSA11q22-q24 were directly mapped within BAC contigs of GGA24. These results provide a better alignment of GGA24 with the corresponding regions in human and mouse and identify several intrachromosomal rearrangements between chicken and mammals.