A novel snoRNA (U73) is encoded within the introns of the human and mouse ribosomal protein S3a genes

The mouse ribosomal protein S3a-encoding gene (mRPS3a) was cloned and sequenced in this study. mRPS3a shares identical exon/intron structure with its human counterpart. Both genes are split to six exons and exhibit remarkable conservation of the promoter region (68.8% identity in the 250 bp upstream of cap site) and coding region (the proteins differ in two amino acids). mRPS3a displays many features common to other r-protein genes, including the CpG-island at 5′-end of the gene, cap site within an oligopyrimidine tract and no consensus TATA or CAAT boxes. However, mRPS3a represents a rare subclass of r-protein genes that possess a long coding sequence in the first exon. Comparison of human and mouse S3a genes revealed sequence fragments with striking similarity within introns 3 and 4. Here we demonstrate that these sequences encode for a novel small nucleolar RNA (snoRNA) designated U73. U73 contains C, D and D′ boxes and a 12-nucleotide antisense complementarity to the 28S ribosomal RNA. These features place U73 into the family of intron-encoded antisense snoRNAs that guide site-specific 2′-O-ribose methylation of pre-rRNA. We propose that U73 is involved in methylation of the G1739 residue of the human 28S rRNA. In addition, we present the mapping of human ribosomal protein S3a gene (hRPS3a) and internally nested U73 gene to the human chromosome 4q31.2–3.     

Distributions of exons and introns in the human genome

The human genome is revisited using exon and intron distribution profiles. The 26,564 annotated genes in the human genome (build October, 2003) contain 233,785 exons and 207,344 introns. On average, there are 8.8 exons and 7.8 introns per gene. About 80% of the exons on each chromosome are < 200 bp in length. < 0.01% of the introns are < 20 bp in length and < 10% of introns are more than 11,000 bp in length. These results suggest constraints on the splicing machinery to splice out very long or very short introns and provide insight to optimal intron length selection. Interestingly, the total length in introns and intergenic DNA on each chromosome is significantly correlated to the determined chromosome size with a coefficient of correlation r = 0.95 and r = 0.97, respectively. These results suggest their implication in genome design.     

Unique features of Chinese hamster S13 gene relative to its human and Xenopus analogs.

We have cloned and sequenced the ribosomal protein S13 gene from the Chinese hamster fibroblast HA-1 cells. The predicted protein encoded by this gene is identical to the human ribosomal protein S13, except for one amino acid substitution at residue 29, which is an alanine in the hamster protein and a threonine in that of humans. The physical organization of the six exons and five introns in the hamster S13 gene is also identical to that found in the human and Xenopus genes with respect to the amino acid codes, even though there are small differences in the lengths of the introns. The striking feature is that unlike its human and Xenopus counterparts, which encode two U14 snoRNAs in two separate introns, the hamster S13 gene encodes no U14 snoRNA. Instead, the hamster gene has a pseudo-U14 coding sequence in its third intron. Our data support the idea that the single copy of the hsc70/U14 gene, which we had previously characterized, is the only source for the production of both U14 snoRNA and hsc70 mRNA species in hamster HA-1 cells.     

The human cytochrome P450 3A locus. Gene evolution by capture of downstream exons

Using a bacterial artificial chromosome (BAC) clone, we have mapped the human cytochrome P450 3A (CYP3A) locus containing the genes encoding for CYP3A4, CYP3A5 and CYP3A7. The genes lie in a head-to-tail orientation in the order of 3A4, 3A7 and 3A5. In both intergenic regions (3A4–3A7 and 3A7–3A5), we have detected several additional cytochrome P450 3A exons, forming two CYP3A pseudogenes. These pseudogenes have the same orientation as the CYP3A genes. To our surprise, a 3A7 mRNA species has been detected in which the exons 2 and 13 of one of the pseudogenes (the one that is downstream of 3A7) are spliced after the 3A7 terminal exon. This results in an mRNA molecule that consists of the 13 3A7 exons and two additional exons at the 3′ end. The additional two exons originating from the pseudogene are in an altered reading frame and consequently have the capability to code a completely different amino acid sequence than the canonical CYP3A exons 2 and 13. These findings may represent a generalized evolutionary process with genes having the potential to capture neighboring sequences and use them as functional exons.     

Genomic organization of mouse Dishevelled genes

We report the isolation and characterization of two new genomic loci corresponding to the mouse Dishevelled (Dvl) genes Dvl2 and Dvl3. The Dvl genes are homologs of the Drosophila dsh segment polarity gene, and are involved in the Wnt/wingless signal transduction pathway. Dvl2 and Dvl3 genomic clones were isolated from a mouse 129 strain λFIXII genomic library and have identical exon/intron organization to Dvll. All three Dishevelled genes span 15 exons and 14 introns and have a number of conserved splice junction sites.     

In contrast to the nematode and fruit fly all 9 intron positions of the sea anemone lamin gene are conserved in human lamin genes

We identified the single gene for nuclear lamin in the genome draft of the sea anemone Nematostella vectensis, a member of the cnidaria, a very old metazoan phylum. The gene consists of 10 exons and 9 introns. Strikingly all 9 intron positions are conserved in the human lamin B genes, which have only 1 (lamin B1) or 2 (lamin B2) additional introns. Using the information on neighboring genes we propose that the human lamin B1 gene on chromosome 5 is the true homolog of the Nematostella lamin gene, while the lamin B2 gene on chromosome 19 arose during vertebrate evolution. In marked contrast to this conservation of gene structure are the results in the rapidly evolving genomes of Drosophila and Caenorhabditis elegans. Here the lamin genes have much fewer introns and these occur often at novel positions. In the single nematode lamin gene and the Drosophila lamin Dmo gene no intron position coincides with an intron in the sea anemone lamin gene.     

Complete genomic structure of human rpS3: identification of functional U15b snoRNA in the fifth intron

Analysis of the complete genomic structure of the human ribosomal protein S3 (rpS3) gene revealed the presence of a functional U15b snoRNA gene in its intron. Human ribosomal protein S3 (rpS3) gene of 6115 bp long has been identified to contain six introns and seven exons in this study. The first and fifth introns of human S3 gene contain functional U15 snoRNA genes. Although Xenopus and Fugu counterparts also have six introns and seven exons, S3 gene of Fugu contains two functional U15 snoRNAs in the fourth and sixth introns and two pseudo genes for U15 snoRNAs in the first and fifth introns. In Xenopus S1 gene encoding ribosomal protein S3, however, three of its six introns contain U15 snoRNA gene sequence. Sequence comparison of the U15 genes from Xenopus, Fugu and human revealed that the regions involved in binding to 28S rRNA and the consensus sequence (C, D and D' boxes) for snoRNAs are highly conserved among those genes from these three species. Human U15a and U15b RNAs which are derived from the first and the fifth introns, respectively, have been identified to be functional by microinjection of human U15a and U15b snoRNAs into Xenopus oocyte. Northern blot and primer extension analyses confirm that human U15b snoRNA is expressed in vivo.     

Computer identification of snoRNA genes using a Mammalian Orthologous Intron Database

Based on comparative genomics, we created a bioinformatic package for computer prediction of small nucleolar RNA (snoRNA) genes in mammalian introns. The core of our approach was the use of the Mammalian Orthologous Intron Database (MOID), which contains all known introns within the human, mouse and rat genomes. Introns from orthologous genes from these three species, that have the same position relative to the reading frame, are grouped in a special orthologous intron table. Our program SNO.pl searches for conserved snoRNA motifs within MOID and reports all cases when characteristic snoRNA-like structures are present in all three orthologous introns of human, mouse and rat sequences. Here we report an example of the SNO.pl usage for searching a particular pattern of conserved C/D-box snoRNA motifs (canonical C- and D-boxes and the 6 nt long terminal stem). In this computer analysis, we detected 57 triplets of snoRNA-like structures in three mammals. Among them were 15 triplets that represented known C/D-box snoRNA genes. Six triplets represented snoRNA genes that had only been partially characterized in the mouse genome. One case represented a novel snoRNA gene, and another three cases, putative snoRNAs. Our programs are publicly available and can be easily adapted and/or modified for searching any conserved motifs within mammalian introns.     

Intronic U14 snoRNAs of Xenopus laevis are located in two different parent genes and can be processed from their introns during early oogenesis.

U14 is a member of the rapidly growing family of intronic small nucleolar RNAs (snoRNAs) that are involved in pre-rRNA processing and ribosome biogenesis. These snoRNA species are encoded within introns of eukaryotic protein coding genes and are synthesized via an intron processing pathway. Characterization of Xenopus laevis U14 snoRNA genes has revealed that in addition to the anticipated location of U14 within introns of the amphibian hsc70 gene (introns 4, 5 and 7), additional intronic U14 snoRNAs are also found in the ribosomal protein S13 gene (introns 3 and 4). U14 is thus far a unique intronic snoRNA in that it is encoded within two different parent genes of a single organism. Northern blot analysis revealed that U14 snoRNAs accumulate during early oocyte development and are rapidly expressed after the mid-blastula transition of developing embryos. Microinjection of hsc70 pre-mRNAs into developing oocytes demonstrated that oocytes as early as stages II and III are capable of processing U14 snoRNA from the pre-mRNA precursor. The ability of immature oocytes to process intronic snoRNAs is consistent with the observed accumulation of U14 during oocyte maturation and the developmentally regulated synthesis of rRNA during oogenesis.     

The organization of the gene for the human cerebroside sulfate activator protein

The organization of 14 exons covering 97% of the cDNA sequence of human cerebroside sulfate activator protein precursor has been determined from two overlapping EMBL-4 human genomic clones extending over 17kb. All exons and exon/intron splice junctions and five introns were sequenced. Exon 8 consists of only 9 bp and is involved in alternative splicing which generates three different mRNAs of cerebroside sulfate activator precursor.     

鳜鱼胰蛋白酶和淀粉酶与胃蛋白酶原基因的克隆与序列分析

采用RT-PCR及RACE法,克隆得到鳜鱼(Siniperca chuatsi)肝胰脏胰蛋白酶(trypsin, Try)、淀粉酶(amylase, Amy)基因 cDNA全序列.结果表明,鳜鱼Try基因cDNA全长为896 bp,其中开放阅读框 (open reading frame,ORF)为744 bp,编码247个氨基酸. 序列同源性分析发现,鳜鱼Try与 斑马鱼(Danio rerio)、非洲爪蟾(Xenopus laevis)、 小鼠Try和人TRY氨基酸序列同源性分别为81.4%、75.3%、74.5%和71.4%.鳜鱼Amy 基因cDNA全长为1 647 bp,其中ORF为1 539 bp,编码512个氨基酸.鳜鱼Amy与斑马鱼 、非洲爪蟾、小鼠Amy和人AMY氨基酸序列同源性分别为79.7%、75.4%、71.9%和70.9%. 同时对鳜鱼基因组进行PCR,获得鳜鱼Try、Amy与胃蛋白酶原(pepsinogen, Pep)全基因组DNA序列.序列分析表明,鳜鱼Try基因由4个内含子和5个外显子组成,全长1 362 bp;鳜鱼Amy基因由8个内含子和9个外显子组成,全长4 267 bp;鳜鱼Pep基因由8个内含子和9个外显子组成,全长 4 032 bp,与其它脊椎动物基因结构相似.应用Genome walker方法在鳜鱼克隆得到长度分别为1 189 bp、413 bp和527 bp的Try、Amy和Pep基因的5′侧翼区序列以及1段长为704 bp的Pep 基因3′侧翼区序列,并利用相关软件预测其中具有多个可调节其表达的调控元件.鳜鱼Try、Am y和Pep基因组全序列的克隆及其序列、结构分析和分子系统进化等的研究,为鱼类消化代谢相关基因的生理功能及表达调控机理进一步研究提供依据.     

Structure of the rat p53 tumor suppressor gene.

Aberration within the p53 tumor suppressor gene is the most frequently identified genetic damage in human cancer. Regulatory functions proposed for the p53 protein include modulation of the cell cycle, cellular differentiation, signal transduction, and gene expression. Additionally, the p53 gene product may guard the genome against incorporation of damaged DNA. To facilitate study of its role in carcinogenesis using a common animal model, we determined the structure of the rat p53 gene. We identified 18 splice sites and defined 25 bases of the intervening sequences adjacent to these sites. We also discovered an allelic polymorphism that occurs within intron 5 of the gene. The rat gene approximates the mouse ortholog. It is 12 kb in length with the non-coding exon 1 separated from exon 2 by 6.2 kb of intervening sequence. The location and size of all rat gene introns approximate those of the mouse. Whereas the mouse and human genes each contain 11 exons, the rat p53 gene is composed of only 10. No intervening sequence occurs between the region of the rat gene corresponding to exons 6 and 7 of the mouse and human p53 genes. This implies intron 6 may be functionally insignificant for species in which it is retained. To extrapolate to p53 involvement in human tumorigenesis, we suggest that mutational events within intron 6 may not be of pathological significance unless splicing is hindered.     

Genomic organization of four β-1,4-endoglucanase genes in plant-parasitic cyst nematodes and its evolutionary implications

The genomic organization of genes encoding β-1,4-endoglucanases (cellulases) from the plant-parasitic cyst nematodes Heterodera glycines and Globodera rostochiensis (HG-eng1, Hg-eng2, GR-eng1, and GR-eng2) was investigated. HG-eng1 and GR-eng1 both contained eight introns and structural domains of 2151 and 2492 bp, respectively. HG-eng2 and GR-eng2 both contained seven introns and structural domains of 2324 and 2388 bp, respectively. No significant similarity in intron sequence or size was observed between HG-eng1 and HG-eng2, whereas the opposite was true between GR-eng1 and GR-eng2. Intron positions among all four cyst nematode cellulase genes were conserved identically in relation to the predicted amino acid sequence. HG-eng1, GR-eng1, and GR-eng2 had several introns demarcated by 5′-GC…AG-3′ in the splice sites, and all four nematode cellulase genes had the polyadenylation and cleavage signal sequence 5′-GAUAAA-3′—both rare occurences in eukaryotic genes. The 5′- flanking regions of each nematode cellulase gene, however, had signature sequences typical of eukaryotic promoter regions, including a TATA box, bHLH-type binding sites, and putative silencer, repressor, and enhancer elements. Database searches and subsequent phylogenetic comparison of the catalytic domain of the nematode cellulases placed the nematode genes in one group, with Family 5, subfamily 2, glycosyl hydrolases from Scotobacteria and Bacilliaceae as the most homologous groups. The overall amino acid sequence identity among the four nematode cellulases was from 71 to 83%, and the amino acid sequence identity to bacterial Family 5 cellulases ranged from 33 to 44%. The eukaryotic organization of the four cyst nematode cellulases suggests that they share a common ancestor, and their strong homology to prokaryotic glycosyl hydrolases may be indicative of an ancient horizontal gene transfer.     

The yeast Hansenula wingei U3 snoRNA gene contains an intron and its coding sequence co-evolved with the 5' ETS region of the pre-ribosomal RNA.

The 5' external transcribed spacer (ETS) region of the pre-rRNA in Saccharomyces cerevisiae contains a sequence with 10 bp of perfect complementarity to the U3 snoRNA. Base pairing between these sequences has been shown to be required for 18S rRNA synthesis, although interaction over the full 10 bp of complementarity is not required. We have identified the homologous sequence in the 5' ETS from the evolutionarily distant yeast Hansenula wingei; unexpectedly, this shows two sequence changes in the region predicted to base pair to U3. By PCR amplification and direct RNA sequencing, a single type of U3 snoRNA coding sequence was identified in H. wingei. As in the S. cerevisiae U3 snoRNA genes, it is interrupted by an intron with features characteristic of introns spliced in a spliceosome. Consequently, this unusual property is not restricted to the yeast genus Saccharomyces. The introns of the H. wingei and S. cerevisiae U3 genes show strong differences in length and sequence, but are located at the same position in the U3 sequence, immediately upstream of the phylogenetically conserved Box A region. The 3' domains of the H. wingei and S. cerevisiae U3 snoRNAs diverge strongly in primary sequence, but have very similar predicted secondary structures. The 5' domains, expected to play a direct role in pre-ribosomal RNA maturation, are more conserved. The sequence predicted to base pair to the pre-rRNA contains two nucleotide substitutions in H. wingei that restore 10 bp of perfect complementarity to the 5' ETS. This is a strong phylogenetic evidence for the importance of the U3/pre-rRNA interaction.     

Large-scale analysis of the genetic and epigenetic alterations in hepatocellular carcinoma from Southeast China

Our knowledge about molecular alterations during hepatocarcinogenesis is still fragmentary, due to lack of comprehensive genetic and epigenetic analyses in the same set of hepatocellular carcinomas (HCCs). In this study, we conducted a large-scale analysis, including mutation screening in 50 genes and methylation assays in three genes in 54 pairs of HCCs and their neighboring non-cancerous tissues. All samples were collected from the residents in Southeast China. We found HBV infection and chronic hepatitis/cirrhosis in 83.3% and 98.1% of the cases, respectively. Mutations were identified in 18 out of 54 (33.3%) samples, with p53 alterations in 14 cases and β-catenin mutations in four tumors. No mutations were identified in the neighboring tissues. Interestingly, 9 out of 14 (64.3%) tumors carrying p53 mutations displayed substitution of serine by arginine at codon 249, a characteristic change believed to be induced by aflatoxin-B1. Furthermore, p53 mutation was significantly associated with shorter recurrence-free survival (P = 0.004). The results also revealed aberrant methylation in two or more genes in as high as 90% of tumors and 40% of adjacent tissues. The frequency of RASSF1A hypermethylation was much higher than that of p16INK4a and HAI2 in both HCC and neighboring tissues, indicating that deregulation of RASSF1A may precede the other two genes. These data suggest that aberrant methylation occurs before mutation and is an early event in the development of this set of HCC. Our findings highlight p53 as a prognostic factor of HCC and RASSF1A as a potential target in preventing malignant transformation of hepatocytes.