A map of 75 human ribosomal protein genes

Human and mouse protein-noncoding snoRNA host genes with dissimilar nucleotide sequences show chromosomal synteny

snoRNAs are small protein-noncoding RNAs essential for pre-rRNA processing and ribosome biogenesis, and are encoded intronically in host genes (HGs) that are either protein coding or noncoding. mRNAs of protein-noncoding HGs differ in their nucleotide sequences among species. Although the reason for such sequential divergence has not been well explained, we present evidence here that such structurally different HGs have evolved from a common ancestral gene. We first identified two novel protein-noncoding HGs (mU50HG-a and mU50HG-b) that intronically encode a mouse ortholog of a human snoRNA, hU50. The sequences of mU50HG mRNA differed from that of hU50HG. However, a chromosome mapping study revealed that mU50HG is located at 9E3-1, the murine segment syntenic to human 6q15, where hU50HG is located. Synteny is a phenomenon whereby gene orthologs are arranged in the same order at equivalent chromosomal loci in different species; synteny between two species means it is highly likely that the genes have evolved from a common ancestral gene. We then extended this mapping study to other protein-noncoding snoRNA-HGs, and found again that they are syntenic, implying that they have evolved from genes of common ancestral species. Furthermore, on these syntenic segments, exons of adjacent protein-coding genes were found to be far better conserved than those of noncoding HGs, suggesting that the exons of protein-noncoding snoRNA-HGs have been much more fragile during evolution.

Global growth deficiencies in mice lacking the ribosomal protein HIP/RPL29

Because of their deleterious effects on developing organisms, ribosomal protein (RP) mutations have been poorly described in mammals, and only a few heterozygous mutations have been shown to be viable. This observation is believed to be due to the fact that each RP is an essential component in the assembly of a functional stable ribosome. Here, we created gene targeted mutant mice lacking HIP/RPL29, an RP associated with translationally active ribosomes in eukaryotes. In contrast to other RP mutants, HIP/RPL29 null mice are viable but are up to 50% smaller than their control littermates at weaning age. In null embryos, delayed global growth is first observed around mid-gestation, and postnatal lethality due to low birth weight results in distortion of the Mendelian ratio. Prenatal growth defects are not fully compensated for during adulthood, and null animals display proportionately smaller organs and stature, and reach sexual maturity considerably later when compared with their control siblings. Additionally, HIP/RPL29 null embryonic fibroblasts have decreased rates of proliferation and protein synthesis and exhibit reduced steady state levels of core RPs. Altogether, our findings provide conclusive genetic evidence that HIP/RPL29 functions as an important regulator of global growth by modulating the rate of protein synthesis.

Characterization of the Full-length cDNA, Chromosomal Localization, and Polymorphism of the Porcine RPLP0 Gene

RPLP0 gene encodes the acidic ribosomal phosphoprotein large P0 subunit, which is a component of the 60S subunit. The full-length cDNA sequence of porcine RPLP0 was obtained from skeletal muscle of fetal pig cDNA library and deposited in GenBank. The nucleotide sequence and the predicted protein sequence shared high sequence identity with other mammalian homologues. A C/A single nucleotide substitution in exon 5 was detected as Csp6?polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) shows allele frequency diversity among Tongcheng, Xiaomeishan, Yushan, Large White, Landrace, and Duroc breeds. Analyses of somatic cell hybrid panel (SCHP) and radiation hybrid (IMpRH) panel showed that the RPLP0 gene was mapped to SSC 14q22-q24 and was closely linked to locus SW1321 (25 cR, LOD = 14.54).

Novel Radiation Response Genes Identified in Gene-Trapped MCF10A Mammary Epithelial Cells

We have used a gene-trapping strategy to screen human mammary epithelial cells for radiation response genes. Relative mRNA expression levels of five candidate genes in MCF10A cells were analyzed, both with and without exposure to radiation. In all five cases, the trapped genes were significantly down-regulated after radiation treatment. Sequence analysis of the fusion transcripts identified the trapped genes: (1) the human androgen receptor, (2) the uncharacterized DREV1 gene, which has known homology to DNA methyltransferases, (3) the human creatine kinase gene, (4) the human eukaryotic translation elongation factor 1 beta 2, and (5) the human ribosomal protein L27. All five genes were down-regulated significantly after treatment with varying doses of ionizing radiation (0.10 to 4.0 Gy) and at varying times (2-30 h after treatment). The genes were also analyzed in human fibroblast and lymphoblastoid cell lines to determine whether the radiation response being observed was cell-type specific. The results verified that the observed radiation response was not a cell-type-specific phenomenon, suggesting that the genes play essential roles in the radiation damage control pathways. This study demonstrates the potential of the gene-trap approach for the identification and functional analysis of novel radiation response genes.

Ribosomal protein gene knockdown causes developmental defects in zebrafish

The ribosomal proteins (RPs) form the majority of cellular proteins and are mandatory for cellular growth. RP genes have been linked, either directly or indirectly, to various diseases in humans. Mutations in RP genes are also associated with tissue-specific phenotypes, suggesting a possible role in organ development during early embryogenesis. However, it is not yet known how mutations in a particular RP gene result in specific cellular changes, or how RP genes might contribute to human diseases. The development of animal models with defects in RP genes will be essential for studying these questions. In this study, we knocked down 21 RP genes in zebrafish by using morpholino antisense oligos to inhibit their translation. Of these 21, knockdown of 19 RPs resulted in the development of morphants with obvious deformities. Although mutations in RP genes, like other housekeeping genes, would be expected to result in nonspecific developmental defects with widespread phenotypes, we found that knockdown of some RP genes resulted in phenotypes specific to each gene, with varying degrees of abnormality in the brain, body trunk, eyes, and ears at about 25 hours post fertilization. We focused further on the organogenesis of the brain. Each knocked-down gene that affected the morphogenesis of the brain produced a different pattern of abnormality. Among the 7 RP genes whose knockdown produced severe brain phenotypes, 3 human orthologs are located within chromosomal regions that have been linked to brain-associated diseases, suggesting a possible involvement of RP genes in brain or neurological diseases. The RP gene knockdown system developed in this study could be a powerful tool for studying the roles of ribosomes in human diseases.

Development of a new set of reference genes for normalization of real-time RT-PCR data of porcine backfat and longissimus dorsi muscle, and evaluation with PPARGC1A

Background

An essential part of using real-time RT-PCR is that expression results have to be normalized before any conclusions can be drawn. This can be done by using one or multiple, validated reference genes, depending on the desired accuracy of the results. In the pig however, very little information is available on the expression stability of reference genes. The aim of this study was therefore to develop a new set of reference genes which can be used for normalization of mRNA expression data of genes expressed in porcine backfat and longissimus dorsi muscle, both representing an economically important part of a pig's carcass. Because of its multiple functions in fat metabolism and muscle fibre type composition, peroxisome proliferative activated receptor γ coactivator 1α (PPARGC1A) is a very interesting candidate gene for meat quality, and was an ideal gene to evaluate our developed set of reference genes for normalization of mRNA expression data of both tissue types.

Results

The mRNA expression stability of 10 reference genes was determined. The expression of RPL13A and SDHA appeared to be highly unstable. After normalization to the geometric mean of the three most stably expressed reference genes (ACTB, TBP and TOP2B), the results not only showed that the mRNA expression of PPARGC1A was significantly higher in each of the longissimus dorsi muscle samples than in backfat (P < 0.05), but also that the expression was significantly higher in the most cranial of the three muscle samples (P < 0.05).

Conclusion

This study provides a new set of reference genes (ACTB, TBP and TOP2B) suitable for normalization of real-time RT-PCR data of backfat and longissimus dorsi muscle in the pig. The obtained PPARGC1A expression results, after application of this set of reference genes, are a first step in unravelling the PPARGC1A expression pattern in the pig and provide a basis for possible selection towards improved meat quality while maintaining a lean carcass.

Rapid detection of trisomy 21 by gene dosage analysis using quantitative real-time polymerase chain reaction

AIM

Rapid detection of fetal aneuploidy helps inform a mother's choice about the course of her pregnancy. Obtaining results by fluorescent in situ hybridization (FISH) requires more than 24 h, and thus a more rapid method is needed.

METHODS

Conventional G-banding and FISH for chromosome 21 were performed for cultured amniocytes. Genomic DNA was extracted from uncultured amniocytes obtained from 23 patients. TaqMan polymerase chain reaction (PCR) primers were designed to amplify the potassium voltage gated channel gene on chromosome 21q22.12 and the ribosomal phosphoprotein gene on 18q21.1. Quantitative real-time PCR was performed for these two gene fragments and the differences of the threshold cycle (Ct) of the two genes (Ct 18-Ct 21) were calculated for each sample.

RESULTS

G-banding revealed that 19 patients had a normal karyotype and four had trisomy 21. FISH resulted in one case of a false positive. The Delta Ct values (Ct 18-Ct 21) of trisomy 21 patients were significantly higher than the values of individuals with normal karyotypes (P < 0.001) and there was no overlapping.

CONCLUSIONS

Fetal trisomy 21 is rapidly detectable by gene dosage analysis from amniocytes using quantitative real-time PCR.

Gene identification and analysis of transcripts differentially regulated in fracture healing by EST sequencing in the domestic sheep

BACKGROUND

The sheep is an important model animal for testing novel fracture treatments and other medical applications. Despite these medical uses and the well known economic and cultural importance of the sheep, relatively little research has been performed into sheep genetics, and DNA sequences are available for only a small number of sheep genes.

RESULTS

In this work we have sequenced over 47 thousand expressed sequence tags (ESTs) from libraries developed from healing bone in a sheep model of fracture healing. These ESTs were clustered with the previously available 10 thousand sheep ESTs to a total of 19087 contigs with an average length of 603 nucleotides. We used the newly identified sequences to develop RT-PCR assays for 78 sheep genes and measured differential expression during the course of fracture healing between days 7 and 42 postfracture. All genes showed significant shifts at one or more time points. 23 of the genes were differentially expressed between postfracture days 7 and 10, which could reflect an important role for these genes for the initiation of osteogenesis.

CONCLUSION

The sequences we have identified in this work are a valuable resource for future studies on musculoskeletal healing and regeneration using sheep and represent an important head-start for genomic sequencing projects for Ovis aries, with partial or complete sequences being made available for over 5,800 previously unsequenced sheep genes.

Reduced gene expression of clustered ribosomal proteins in Diamond-Blackfan anemia patients without RPS19 gene mutations

Diamond-Blackfan anemia (DBA) is a rare congenital pure red cell aplasia occasionally presenting physical anomalies. Ribosomal protein S19 gene (RPS19) is one of the causative genes for DBA; however, the pathologic mechanism of erythroblastopenia and abnormal morphology has not been clarified. To assess the pathophysiology of DBA, the gene expression profile of 2 representative patients carrying no RPS19 mutations was compared with that of aplastic anemia (AA) patients, assessed by the microarray analyses. The K-mean clustering analysis revealed the significant categorization of 28 ribosomal protein (RP) genes into a small set of group (994 genes) (P=2.39E-17), all of which were expressed at lower levels in DBA than in AA patients. RPS19 was categorized into the set of low expressing genes in DBA patients. No mutations were determined in the promoter and coding sequences of top 10 RP genes expressed at the levels over 1.2 of the AA/DBA ratio, in 3 DBA patients. These results indicated that the lower expression of RP gene group, even without the mutation, was a distinctive feature of DBA from AA, although the study number was small. The reduced RP gene expression, by itself, may suggest an underlying mechanism of the constitutional anemia.

Effects of Scutellariae Radix on gene expression in HEK 293 cells using cDNA microarray

The aim of the present study was to elucidate the molecular mechanisms underlying the anti-inflammatory effect of Scutellaria Radix (SR). The complementary DNA (cDNA) microarray method was used to survey the effects of SR on the changes of gene expression profile in HEK293 cells. Based on differential expression, 66 genes were selected for further analysis from 9,600 candidate genes in the microarray; 23 genes were validated by RT-PCR. The broad spectrum of the differentially expressed genes, including those associated with inflammation, immune response, energy metabolism, as well as others, such as ISGF3G, IL6ST, CD98, ATP5G2, PHKG2, YB-1 and SLC7A4, indicate overall cellular response to SR treatment. Our results suggest that the anti-inflammatory effect of SR may be related to IL6ST down-expression, and over-expression of CD98. Moreover, SR-related improvement in immune response may be related to the ISGF3G over-expression.

The complete set of Toxoplasma gondii ribosomal protein genes contains two conserved promoter elements

Recently we showed that de novo ribosome biosynthesis is transcriptionally regulated in Coccidia, depending on their life-cycle stage. Since the expression of ribosomal protein genes is likely coordinated, the transcriptional control of all Toxoplasma gondii ribosomal protein (RP) genes was analysed. Therefore, the complete set of all cytoplasmic RPs was defined, containing 79 different RPs in T. gondii. RP genes were randomly distributed over the genome, each with a unique upstream region with the exception of 8 RP genes which were paired in a head-to-head orientation. To study if the RP genes share conserved promoter elements, a database was made containing upstream sequences of all T. gondii RP genes. Promoter activity was confirmed for the upstream sequences of 8 RP genes, some of which are comparable in strength to the alpha-tubulin promoter. In the complete set of RP upstream sequences 2 novel and highly conserved elements were identified, named Toxoplasma Ribosomal Protein (TRP)-1 (consensus: TCGGCTTATATTCGG) and TRP-2 ([T/C]GCATGC[G/A]). TRP-1 and/or TRP-2 were present in 95% of all RP upstream sequences and moreover, were specifically localized in a small region near the presumptive transcriptional start site (10-330 bp upstream). Although TRP elements were mostly absent in known T. gondii promoters, they are present elsewhere in the T. gondii genome suggesting that they operate not only in RP genes but in a larger set of genes. The identification of TRP elements creates a basis to further study the underlying mechanism by which RP transcription is controlled in T. gondii.

Identification of differently expressed genes in human colorectal adenocarcinoma

AIM: To investigate the differently expressed genes in human colorectal adenocarcinoma.

METHODS: The integrated approach for gene expression profiling that couples suppression subtractive hybridization, high-throughput cDNA array, sequencing, bioinformatics analysis, and reverse transcriptase real-time quantitative polymerase chain reaction (PCR) was carried out. A set of cDNA clones including 1260 SSH inserts amplified by PCR was arrayed using robotic printing. The cDNA arrays were hybridized with florescent-labeled probes prepared from RNA of human colorectal adenocarcinoma (HCRAC) and normal colorectal tissues.

RESULTS: A total of 86 genes were identified, 16 unknown genes and 70 known genes. The transcription factor Sox9 influencing cell differentiation was downregulated. At the same time, Heat shock protein 10 KDis downregulated and Calmoulin is up-regulated.

CONCLUSION: Downregulation of heat shock protein 10 KD lost its inhibition of Ras, and then attenuated the Ras GTPase signaling pathway, increased cell proliferation and inhibited cell apoptosis. Down-regulated transcription factor So x 9 influences cell differentiation and cell-specific gene expression. Down-regulated So x 9 also decreases its binding to calmodulin, accumulates calmodulin as receptor-activated kinase and phosphorylase kinase due to the activation of PhK.

Methylthioadenosine phosphorylase gene expression is impaired in human liver cirrhosis and hepatocarcinoma

Methylthioadenosine phosphorylase (MTAP) is a key enzyme in the methionine and adenine salvage pathways. In mammals, the liver plays a central role in methionine metabolism, and this essential function is lost in the progression from liver cirrhosis to hepatocarcinoma. Deficient MTAP gene expression has been recognized in many transformed cell lines and tissues. In the present work, we have studied the expression of MTAP in human and experimental liver cirrhosis and hepatocarcinoma. We observe that MTAP gene expression is significantly reduced in human hepatocarcinoma tissues and cell lines. Interestingly, MTAP gene expression was also impaired in the liver of CCl4-cirrhotic rats and cirrhotic patients. We provide evidence indicating that epigenetic mechanisms, involving DNA methylation and histone deacetylation, may play a role in the silencing of MTAP gene expression in hepatocarcinoma. Given the recently proposed tumor suppressor activity of MTAP, our observations can be relevant to the elucidation of the molecular mechanisms of multistep hepatocarcinogenesis.

Nuclear factor-(kappa)B modulates the p53 response in neurons exposed to DNA damage

Previous studies have shown that DNA damage-evoked death of primary cortical neurons occurs in a p53 and cyclin-dependent kinase-dependent (CDK) manner. The manner by which these signals modulate death is unclear. Nuclear factor-kappaB (NF-kappaB) is a group of transcription factors that potentially interact with these pathways. Presently, we show that NF-kappaB is activated shortly after induction of DNA damage in a manner independent of the classic IkappaB kinase (IKK) activation pathway, CDKs, ATM, and p53. Acute inhibition of NF-kappaB via expression of a stable IkappaB mutant, downregulation of the p65 NF-kappaB subunit by RNA interference (RNAi), or pharmacological NF-kappaB inhibitors significantly protected against DNA damage-induced neuronal death. NF-kappaB inhibition also reduced p53 transcripts and p53 activity as measured by the p53-inducible messages, Puma and Noxa, implicating the p53 tumor suppressor in the mechanism of NF-kappaB-mediated neuronal death. Importantly, p53 expression still induces death in the presence of NF-kappaB inhibition, indicating that p53 acts downstream of NF-kappaB. Interestingly, neurons cultured from p65 or p50 NF-kappaB-deficient mice were not resistant to death and did not show diminished p53 activity, suggesting compensatory processes attributable to germline deficiencies, which allow p53 activation still to occur. In contrast to acute NF-kappaB inhibition, prolonged NF-kappaB inhibition caused neuronal death in the absence of DNA damage. These results uniquely define a signaling paradigm by which NF-kappaB serves both an acute p53-dependent pro-apoptotic function in the presence of DNA damage and an anti-apoptotic function in untreated normal neurons.

Expression Profiling of Ameloblastic Carcinoma

Ameloblastic carcinoma (AC) is a malignant epithelial odontogenic tumor that histologically retains the features of ameloblastic differentiation and exhibits cytological features of malignancy in the primary or recurrent tumor. It may develop within a preexisting ameloblastoma or arise de novo or from an odontogenic cyst. Expression profiling by DNA microarray is a new molecular technology that allows the analysis of cell and tissue gene expression. By using DNA microarrays containing 19,200 genes, several genes whose expression was significantly upregulated or downregulated were identified in a case of AC. The differentially expressed genes cover a broad range of functional activities: 1) transcription, 2) signaling transduction, 3) cell cycle regulation, 4) apoptosis control, and 5) differentiation. The data reported are, to our knowledge, the first genetic portrait of an AC. No final conclusion can be drawn; however, this portrait will be useful in investigating the biological behavior and in identifying possible gene targets for cancer therapy when more cases of this rare tumor are reported and compared.

X chromosomes, retrogenes and their role in male reproduction

Retrogenes originate from their progenitor genes by retroposition. Several retrogenes reported in recent studies are autosomal, originating from X-linked progenitor genes, and have evolved a testis-specific expression pattern. During male meiosis, sex chromosomes are segregated into a so-called 'XY' body and are silenced transcriptionally. It has been widely hypothesized that the silencing of the X chromosome during male meiosis is the driving force behind the retroposition of X-linked genes to autosomes during evolution. With the advent of sequenced genomes of many species, many retrogenes can be identified and characterized. The testis-specific retrogenes might be associated with human male infertility. My goal here is to integrate recent findings, highlight controversies in the field and identify areas for further study.

Annexin A2 Is a Novel RNA-binding Protein

Annexin A2 (ANXA2) is a Ca(2+)-binding protein that is up-regulated in virally transformed cell lines and in human tumors. Here, we show that ANXA2 binds directly to both ribonucleotide homopolymers and human c-myc RNA. ANXA2 was shown to bind specifically to poly(G) with high affinity (K(d) = 60 nM) and not to poly(A), poly(C), or poly(U). The binding of ANXA2 to poly(G) required Ca(2+) (A(50%) = 10 microM). The presence of RNA in the immunoprecipitates of ANXA2 isolated from HeLa cells established that ANXA2 formed a ribonucleoprotein complex in vivo. Sucrose gradient analysis showed that ANXA2 associates with ribonucleoprotein complexes and not with polyribosomes. Reverse transcriptase-PCR identified c-myc mRNA as a component of the ribonucleoprotein complex formed by ANXA2 in vivo, and binding studies confirmed a direct interaction between ANXA2 and c-myc mRNA. Transfection of LNCaP cells with the ANXA2 gene resulted in the up-regulation of c-Myc protein. These findings identify ANXA2 as a Ca(2+)-dependent RNA-binding protein that interacts with the mRNA of the nuclear oncogene, c-myc.

RPG: the Ribosomal Protein Gene database

RPG (http://ribosome.miyazaki-med.ac.jp/) is a new database that provides detailed information about ribosomal protein (RP) genes. It contains data from humans and other organisms, including Drosophila melanogaster, Caenorhabditis elegans, Saccharo myces cerevisiae, Methanococcus jannaschii and Escherichia coli. Users can search the database by gene name and organism. Each record includes sequences (genomic, cDNA and amino acid sequences), intron/exon structures, genomic locations and information about orthologs. In addition, users can view and compare the gene structures of the above organisms and make multiple amino acid sequence alignments. RPG also provides information on small nucleolar RNAs (snoRNAs) that are encoded in the introns of RP genes.

Millions of years of evolution preserved: a comprehensive catalog of the processed pseudogenes in the human genome

Processed pseudogenes were created by reverse-transcription of mRNAs; they provide snapshots of ancient genes existing millions of years ago in the genome. To find them in the present-day human, we developed a pipeline using features such as intron-absence, frame-disruption, polyadenylation, and truncation. This has enabled us to identify in recent genome drafts approximately 8000 processed pseudogenes (distributed from http://pseudogene.org). Overall, processed pseudogenes are very similar to their closest corresponding human gene, being 94% complete in coding regions, with sequence similarity of 75% for amino acids and 86% for nucleotides. Their chromosomal distribution appears random and dispersed, with the numbers on chromosomes proportional to length, suggesting sustained "bombardment" over evolution. However, it does vary with GC-content: Processed pseudogenes occur mostly in intermediate GC-content regions. This is similar to Alus but contrasts with functional genes and L1-repeats. Pseudogenes, moreover, have age profiles similar to Alus. The number of pseudogenes associated with a given gene follows a power-law relationship, with a few genes giving rise to many pseudogenes and most giving rise to few. The prevalence of processed pseudogenes agrees well with germ-line gene expression. Highly expressed ribosomal proteins account for approximately 20% of the total. Other notables include cyclophilin-A, keratin, GAPDH, and cytochrome c.

Heparan sulfate proteoglycans: Coordinators of multiple signaling pathways during chondrogenesis

Heparan sulfate proteoglycans are abundantly expressed in the pericellular matrix of both developing and mature cartilage. Increasing evidence indicates that the action of numerous chondroregulatory molecules depends on these proteoglycans. This review summarizes the current understanding of the interactions of heparan sulfate chains of cartilage proteoglycans with both soluble and nonsoluble ligands during the process of chondrogenesis. In addition, the consequences of mutating genes encoding heparan sulfate biosynthetic enzymes or heparan sulfate proteoglycan core proteins on cartilage development are discussed.

Identification and functional analysis of human transcriptional promoters

Genomic and full-length cDNA sequences provide opportunities for understanding human gene structure and transcriptional regulatory elements. The simplest regulatory elements to identify are promoters, as their positions are dictated by the location of transcription start sites. We aligned full-length cDNA clones from the Mammalian Gene Collection to the human genome rough draft sequence to estimate the start sites of more than 10,000 human transcripts. We selected genomic sequence just upstream from the 5' end of these cDNA sequences and designated these as putative promoters. We assayed the functions of 152 of these DNA fragments, chosen at random from the entire set, in a luciferase-based transfection assay in four human cultured cell types. Ninety-one percent of these DNA fragments showed significant transcriptional activity in at least one of the cell lines, whereas 89% showed activity in at least two of the lines. We analyzed the distributions of strengths of these promoter fragments in the different cell types and identified likely alternative promoters in a large fraction of the genes. These data indicate that this approach is an effective method for predicting human promoters and provide the first set of functional data collected in parallel for a large set of human promoters.

Ribosomal protein S19 expression during erythroid differentiation

The gene encoding ribosomal protein S19 (RPS19) has been shown to be mutated in 25% of the patients affected by Diamond-Blackfan anemia (DBA), a congenital erythroblastopenia. As the role of RPS19 in erythropoiesis is still to be defined, we performed studies on RPS19 expression during terminal erythroid differentiation. Comparative analysis of the genomic sequences of human and mouse RPS19 genes enabled the identification of 4 conserved sequence elements in the 5' region. Characterization of transcriptional elements allowed the identification of the promoter in the human RPS19 gene and the localization of a strong regulatory element in the third conserved sequence element. By Northern blot and Western blot analyses of murine splenic erythroblasts infected with the anemia-inducing strain Friend virus (FAV cells), RPS19 mRNA and protein expression were shown to decrease during terminal erythroid differentiation. We anticipate that these findings will contribute to further development of our understanding of the contribution of RPS19 to erythropoiesis.

Comparative analysis of ribosomal proteins in complete genomes: an example of reductive evolution at the domain scale

A comprehensive investigation of ribosomal genes in complete genomes from 66 different species allows us to address the distribution of r-proteins between and within the three primary domains. Thirty-four r-protein families are represented in all domains but 33 families are specific to Archaea and Eucarya, providing evidence for specialisation at an early stage of evolution between the bacterial lineage and the lineage leading to Archaea and Eukaryotes. With only one specific r-protein, the archaeal ribosome appears to be a small-scale model of the eukaryotic one in terms of protein composition. However, the mechanism of evolution of the protein component of the ribosome appears dramatically different in Archaea. In Bacteria and Eucarya, a restricted number of ribosomal genes can be lost with a bias toward losses in intracellular pathogens. In Archaea, losses implicate 15% of the ribosomal genes revealing an unexpected plasticity of the translation apparatus and the pattern of gene losses indicates a progressive elimination of ribosomal genes in the course of archaeal evolution. This first documented case of reductive evolution at the domain scale provides a new framework for discussing the shape of the universal tree of life and the selective forces directing the evolution of prokaryotes.

Functional second genes generated by retrotransposition of the X-linked ribosomal protein genes

We have identified a new class of ribosomal protein (RP) genes that appear to have been retrotransposed from X-linked RP genes. Mammalian ribosomes are composed of four RNA species and 79 different proteins. Unlike RNA constituents, each protein is typically encoded by a single intron- containing gene. Here we describe functional autosomal copies of the X-linked human RP genes, which we designated RPL10L (ribosomal protein L10-like gene), RPL36AL and RPL39L after their progenitors. Because these genes lack introns in their coding regions, they were likely retrotransposed from X-linked genes. The identities between the retrotransposed genes and the original X-linked genes are 89-95% in their nucleotide sequences and 92-99% in their amino acid sequences, respectively. Northern blot and PCR analyses revealed that RPL10L and RPL39L are expressed only in testis, whereas RPL36AL is ubiquitously expressed. Although the role of the autosomal RP genes remains unclear, they may have evolved to compensate for the reduced dosage of X-linked RP genes.

Identification and analysis of over 2000 ribosomal protein pseudogenes in the human genome

Mammals have 79 ribosomal proteins (RP). Using a systematic procedure based on sequence-homology, we have comprehensively identified pseudogenes of these proteins in the human genome. Our assignments are available at http://www.pseudogene.org or http://bioinfo.mbb.yale.edu/genome/pseudogene. In total, we found 2090 processed pseudogenes and 16 duplications of RP genes. In relation to the matching parent protein, each of the processed pseudogenes has an average relative sequence length of 97% and an average sequence identity of 76%. A small number (258) of them do not contain obvious disablements (stop codons or frameshifts) and, therefore, could be mistaken as functional genes, and 178 are disrupted by one or more repetitive elements. On average, processed pseudogenes have a longer truncation at the 5' end than the 3' end, consistent with the target-primed-reverse-transcription (TPRT) mechanism. Interestingly, on chromosome 16, an RPL26 processed pseudogene was found in the intron region of a functional RPS2 gene. The large-scale distribution of RP pseudogenes throughout the genome appears to result, chiefly, from random insertions with the numbers on each chromosome, consequently, proportional to its size. In contrast to RP genes, the RP pseudogenes have the highest density in GC-intermediate regions (41%-46%) of the genome, with the density pattern being between that of LINEs and Alus. This can be explained by a negative selection theory as we observed that GC-rich RP pseudogenes decay faster in GC-poor regions. Also, we observed a correlation between the number of processed pseudogenes and the GC content of the associated functional gene, i.e., relatively GC-poor RPs have more processed pseudogenes. This ranges from 145 pseudogenes for RPL21 down to 3 pseudogenes for RPL14. We were able to date the RP pseudogenes based on their sequence divergence from present-day RP genes, finding an age distribution similar to that for Alus. The distribution is consistent with a decline in retrotransposition activity in the hominid lineage during the last 40 Myr. We discuss the implications for retrotransposon stability and genome dynamics based on these new findings.

Eukaryotic mRNPs may represent posttranscriptional operons

Genomic array analysis of endogenous mammalian ribonucleoproteins has recently revealed three novel findings: (1) mRNA binding proteins are associated with unique subpopulations of messages, (2) the compositions of these mRNA subsets can vary with growth conditions, and (3) the same mRNA species can be found in multiple mRNP complexes. Based on these and other findings, we propose a model of posttranscriptional gene expression in which mRNA binding proteins regulate mRNAs as subpopulations during cell growth and development. This model predicts that functionally related genes are regulated posttranscriptionally as groups by specific mRNA binding proteins that recognize sequence elements in common among the mRNAs.

The human ribosomal protein genes: sequencing and comparative analysis of 73 genes

The ribosome, as a catalyst for protein synthesis, is universal and essential for all organisms. Here we describe the structure of the genes encoding human ribosomal proteins (RPs) and compare this class of genes among several eukaryotes. Using genomic and full-length cDNA sequences, we characterized 73 RP genes and found that (1) transcription starts at a C residue within a characteristic oligopyrimidine tract; (2) the promoter region is GC rich, but often has a TATA box or similar sequence element; (3) the genes are small (4.4 kb), but have as many as 5.6 exons on average; (4) the initiator ATG is in the first or second exon and is within plus minus 5 bp of the first intron boundaries in about half of cases; and (5) 5'- and 3'-UTRs are significantly smaller (42 bp and 56 bp, respectively) than the genome average. Comparison of RP genes from humans, Drosophila melanogaster, Caenorhabditis elegans, and Saccharomyces cerevisiae revealed the coding sequences to be highly conserved (63% homology on average), although gene size and the number of exons vary. The positions of the introns are also conserved among these species as follows: 44% of human introns are present at the same position in either D. melanogaster or C. elegans, suggesting RP genes are highly suitable for studying the evolution of introns.

Identification of differentially expressed genes in human uterine leiomyomas using differential display

In searching of differentially expressed genes in human uterine leiomyomas, differential display was used with twelve pairs of primers to compare human uterine leiomyomas with matched myometrium. False positives were eliminated by reverse Northern analysis. Positives were confirmed by Northern blot analysis.

RESULTS

Four of 69 cDNA fragments (3 up-regulated named L1, L2 and L3 and 1 down-regulated named M1 in leiomyoma) were confirmed by Northern analysis. Sequence comparison and Northern analysis proved that L1 is exactly the human ribosomal protein S19. It was present ubiquitously in 13 tissues tested but in various levels and even in different size. L1 was highly expressed in parotidean cystadenocarcinoma, pancreatic cancer and breast cancer examined. No mutations have been found in human uterine leiomyomas (n=6).

CONCLUSIONS

hRPS19 overexpression might be a universal signal in rapid cell growth tissues.

Molecular dissection of the human Y-chromosome

Human Y chromosome, earlier thought to be gene deficient, has attracted a great deal of attention owing to its supremacy in male sex determination and unique haplotype status in the genome. Studies on Y chromosome have shown the presence of different types of satellite DNA and several genes implicated with a variety of physical and physiological functions. The interaction of these repetitive DNA with genes in normal individuals and in patients with Y-chromosome-related genetic anomalies is still an unresolved issue and is actively being pursued. The fast changing scenario of the human genome project is likely to effect our overall understanding of the Y chromosome and Y-linked genetic anomalies in a big way. We provide a brief overview of the organization of Y chromosome with respect to several important loci encompassing both the arms and their likely involvement/modulation in genetic anomalies. The experimental approaches discussed here are envisaged to be of clinical relevance for the molecular diagnosis of the Y-linked disorders.

Comparative mapping of RPL3, a gene overexpressed in multiple obesity models

The ribosomal protein 3 gene is differentially expressed in hypothalamus and brown adipose tissue between mouse lines divergently selected for heat loss, and in skeletal muscle of the ob/ob mouse model. Unfortunately, multiple Rpl3-processed pseudogenes have hampered mapping of the functional gene copy in mammalian species. Using PCR amplification with intronic primer binding, we have mapped Rpl3 to MMU15, and have also localized RPL3 to BTA5 in cattle. Comparative mapping implicates a previously mapped copy of RPL3 on HSA22 as the functional copy of human RPL3, while predictive mapping places the porcine homologue on SSC5.

Mapping of 12 bovine ribosomal protein genes using a bovine radiation hybrid panel

Twelve bovine ribosomal protein genes, for which sequence data had been acquired from complementary deoxyribonucleic acid (cDNA) clones isolated from a cattle skin cDNA library, were mapped. As ribosomal protein genes are a group of highly conserved house keeping genes, specific primers were designed to span the intron-exon splice sites and to amplify intronic sequences, in order to obtain bovine-specific polymerase chain reaction (PCR) products. Two of 12 ribosomal protein genes were genotyped in this way and the remaining 10 were mapped using additional primers designed from within the intron. Eleven previously unmapped ribosomal protein genes were localized and one previously reported ribosomal protein gene localization was confirmed. The 12 ribosomal protein genes mapped in this study are spread over 10 chromosomes, including the X chromosome. The locations show conservation of comparative map position in cattle and human.

Screening of a HUVEC cDNA library with transplant-associated coronary artery disease sera identifies RPL7 as a candidate autoantigen associated with this disease

A HUVEC cDNA library was screened with sera from two patients who had developed transplant-associated coronary artery disease (TxCAD) following cardiac transplantation. A total of six positive clones were isolated from a primary screen of 40 000 genes. Subsequent DNA sequence analysis identified these to be lysyl tRNA synthetase, ribosomal protein L7, ribosomal protein L9, beta transducin and TANK. Another gene whose product could not be identified showed homology to a human cDNA clone (DKFZp566M063) derived from fetal kidney. Full-length constructs of selected genes were expressed as his-tag recombinant fusion proteins and used to screen a wider patient base by ELISA to determine prevalence and association with TxCAD. Of these ribosomal protein L7 showed the highest prevalence (55.6%) with TxCAD sera compared to 10% non-CAD.

The organization of cytoplasmic ribosomal protein genes in the Arabidopsis genome

Eukaryotic ribosomes are made of two components, four ribosomal RNAs, and approximately 80 ribosomal proteins (r-proteins). The exact number of r-proteins and r-protein genes in higher plants is not known. The strong conservation in eukaryotic r-protein primary sequence allowed us to use the well-characterized rat (Rattus norvegicus) r-protein set to identify orthologues on the five haploid chromosomes of Arabidopsis. By use of the numerous expressed sequence tag (EST) accessions and the complete genomic sequence of this species, we identified 249 genes (including some pseudogenes) corresponding to 80 (32 small subunit and 48 large subunit) cytoplasmic r-protein types. None of the r-protein genes are single copy and most are encoded by three or four expressed genes, indicative of the internal duplication of the Arabidopsis genome. The r-proteins are distributed throughout the genome. Inspection of genes in the vicinity of r-protein gene family members confirms extensive duplications of large chromosome fragments and sheds light on the evolutionary history of the Arabidopsis genome. Examination of large duplicated regions indicated that a significant fraction of the r-protein genes have been either lost from one of the duplicated fragments or inserted after the initial duplication event. Only 52 r-protein genes lack a matching EST accession, and 19 of these contain incomplete open reading frames, confirming that most genes are expressed. Assessment of cognate EST numbers suggests that r-protein gene family members are differentially expressed.

The human mitochondrial ribosomal protein genes: mapping of 54 genes to the chromosomes and implications for human disorders

Mitochondria possess their own translational machinery, which is composed of components distinct from their cytoplasmic counterparts. To investigate the possible involvement of mitochondrial ribosomal defects in human disease, we mapped nuclear genes that encode mitochondrial ribosomal proteins (MRPs). We generated sequence-tagged sites (STSs) of individual MRP genes that were able to be detected by PCR. They were placed on an STS content map of the human genome by typing of radiation hybrid panels. We located 54 MRP genes on the STS-content map and assigned these genes to cytogenetic bands of the human chromosomes. Although mitochondria are thought to have originated from bacteria, in which the genes encoding ribosomal proteins are clustered into operons, the mapped MRP genes are widely dispersed throughout the genome, suggesting that transfer of each MRP gene to the nuclear genome occurred individually. We compared the assigned positions with candidate regions for mendelian disorders and found certain genes that might be involved in particular diseases. This map provides a basis for studying possible roles of MRP defects in mitochondrial disorders.

Combined copy status of 18q21 genes in colorectal cancer shows frequent retention of SMAD7

Deletions of chromosome band 18q21 appear with very high frequency in a variety of carcinomas, especially in colorectal cancer. Potent tumor suppressor genes located in this region encode transforming growth factor beta (TGF-beta) signal transducers SMAD2 and SMAD4, and inactivation of either one leads to impaired TGF-beta-mediated cell growth/apoptosis. Following the assignment of SMAD7 to 18q21, we first refined the SMAD7 gene position within this region by genetically mapping SMAD7 between SMAD2 and SMAD4. Further, to compare the respective frequencies of genetic alterations of these three SMAD genes in colorectal cancer, we undertook a large-scale evaluation of the copy status of each of these genes on DNA samples from colorectal tumor biopsy material. Among a subset of 233 DNA samples for which data were available for all four genes, SMAD4, SMAD2, and the nearby gene DCC showed high deletion rates (66%, 64%, and 59%, respectively), whereas SMAD7 was deleted in only 48% of the tumors. Unexpectedly, we found some gene duplications; SMAD7 appears to be more frequently amplified (10%) than the three other genes (4-7%). Compiled data for SMAD genes in each tumor show that the most common combination (26% of all the tumors) consists of the simultaneous deletions of SMAD2 and SMAD4 associated with normal diploidy or even duplication of SMAD7. Since SMAD7 normally counteracts SMAD2 and SMAD4 in TGF-beta signaling, we hypothesize that the tumor might not benefit from simultaneous SMAD7 inactivation, thereby exerting selective pressure to retain or even to duplicate the SMAD7 gene.

Epigenetic silencing of PEG3 gene expression in human glioma cell lines

Genomic imprinting, the phenomenon in which alleles of genes are expressed differentially depending on their parental origins, has important consequences for mammalian development, and disturbance of normal imprinting leads to abnormal embryogenesis and some inherited diseases and is also associated with various cancers. In the context of screening for novel imprinted genes on human chromosome 19q13.4 with mouse A9 hybrids, we identified a maternal allele-specific methylated CpG island in exon 1 of paternally expressed imprinted gene 3 (PEG3), a gene that exhibits paternal allele-specific expression. Because PEG3 expression is downregulated in some gliomas and glioma cell lines, despite high-level expression in normal brain tissues, we investigated whether the loss of PEG3 expression is related to epigenetic modifications involving DNA methylation. We found monoallelic expression of PEG3 in all normal brain tissues examined and five of nine glioma cell lines that had both unmethylated and methylated alleles; the remaining four glioma cell lines exhibited gain of imprinting with hypermethylated alleles. In addition, treatment of glioma cell lines with the DNA demethylating agent 5-aza-2'-deoxycytidine reversed the silencing of PEG3 biallelically. In this article, we report that the epigenetic silencing of PEG3 expression in glioma cell lines depends on aberrant DNA methylation of an exonic CpG island, suggesting that PEG3 contributes to glioma carcinogenesis in certain cases.

A complete map of the human ribosomal protein genes: assignment of 80 genes to the cytogenetic map and implications for human disorders

Mapping of the human ribosomal protein (RP) genes has been completed, and all 80 different genes were placed on a cytogenetic map of the human genome. Because of the existence of processed pseudogenes, the localization of the RP genes was complicated, and five genes had remained to be mapped. Here we developed a novel strategy to identify sequence-tagged sites (STSs) at introns of the RP genes, and we localized RPL14, RPL22, RPL35, RPL36, and RPL39 within the chromosomes by radiation hybrid mapping. Unlike the case of eubacteria or archaebacteria, human RP genes are widely scattered about the genome. Together with the previous results, both sex chromosomes and 20 autosomes (all but chromosomes 7 and 21) were found to carry one or more RP genes. To explore the possible involvement of RP genes in human disorders, all 80 genes were assigned to cytogenetic bands according to a published cytogenetic BAC-STS map of the human genome. We compared the assigned positions with candidate regions for Mendelian disorders and found certain genes that might be involved in particular human disorders.

The sequence of the human genome

A 2.91-billion base pair (bp) consensus sequence of the euchromatic portion of the human genome was generated by the whole-genome shotgun sequencing method. The 14.8-billion bp DNA sequence was generated over 9 months from 27,271,853 high-quality sequence reads (5.11-fold coverage of the genome) from both ends of plasmid clones made from the DNA of five individuals. Two assembly strategies-a whole-genome assembly and a regional chromosome assembly-were used, each combining sequence data from Celera and the publicly funded genome effort. The public data were shredded into 550-bp segments to create a 2.9-fold coverage of those genome regions that had been sequenced, without including biases inherent in the cloning and assembly procedure used by the publicly funded group. This brought the effective coverage in the assemblies to eightfold, reducing the number and size of gaps in the final assembly over what would be obtained with 5.11-fold coverage. The two assembly strategies yielded very similar results that largely agree with independent mapping data. The assemblies effectively cover the euchromatic regions of the human chromosomes. More than 90% of the genome is in scaffold assemblies of 100,000 bp or more, and 25% of the genome is in scaffolds of 10 million bp or larger. Analysis of the genome sequence revealed 26,588 protein-encoding transcripts for which there was strong corroborating evidence and an additional approximately 12,000 computationally derived genes with mouse matches or other weak supporting evidence. Although gene-dense clusters are obvious, almost half the genes are dispersed in low G+C sequence separated by large tracts of apparently noncoding sequence. Only 1.1% of the genome is spanned by exons, whereas 24% is in introns, with 75% of the genome being intergenic DNA. Duplications of segmental blocks, ranging in size up to chromosomal lengths, are abundant throughout the genome and reveal a complex evolutionary history. Comparative genomic analysis indicates vertebrate expansions of genes associated with neuronal function, with tissue-specific developmental regulation, and with the hemostasis and immune systems. DNA sequence comparisons between the consensus sequence and publicly funded genome data provided locations of 2.1 million single-nucleotide polymorphisms (SNPs). A random pair of human haploid genomes differed at a rate of 1 bp per 1250 on average, but there was marked heterogeneity in the level of polymorphism across the genome. Less than 1% of all SNPs resulted in variation in proteins, but the task of determining which SNPs have functional consequences remains an open challenge.

The Mouse Brain Transcriptome by SAGE: Differences in Gene Expression between P30 Brains of the Partial Trisomy 16 Mouse Model of Down Syndrome (Ts65Dn) and Normals

Trisomy 21, or Down syndrome (DS), is the most common genetic cause of mental retardation. Changes in the neuropathology, neurochemistry, neurophysiology, and neuropharmacology of DS patients' brains indicate that there is probably abnormal development and maintenance of central nervous system structure and function. The segmental trisomy mouse (Ts65Dn) is a model of DS that shows analogous neurobehavioral defects. We have studied the global gene expression profiles of normal and Ts65Dn male and normal female mice brains (P30) using the serial analysis of gene expression (SAGE) technique. From the combined sample we collected a total of 152,791 RNA tags and observed 45,856 unique tags in the mouse brain transcriptome. There are 14 ribosomal protein genes (nine underexpressed) among the 330 statistically significant differences between normal male and Ts65Dn male brains, which possibly implies abnormal ribosomal biogenesis in the development and maintenance of DS phenotypes. This study contributes to the establishment of a mouse brain transcriptome and provides the first overall analysis of the differences in gene expression in aneuploid versus normal mammalian brain cells.

The mouse brain transcriptome by SAGE: differences in gene expression between P30 brains of the partial trisomy 16 mouse model of Down syndrome (Ts65Dn) and normals

Trisomy 21, or Down syndrome (DS), is the most common genetic cause of mental retardation. Changes in the neuropathology, neurochemistry, neurophysiology, and neuropharmacology of DS patients' brains indicate that there is probably abnormal development and maintenance of central nervous system structure and function. The segmental trisomy mouse (Ts65Dn) is a model of DS that shows analogous neurobehavioral defects. We have studied the global gene expression profiles of normal and Ts65Dn male and normal female mice brains (P30) using the serial analysis of gene expression (SAGE) technique. From the combined sample we collected a total of 152,791 RNA tags and observed 45,856 unique tags in the mouse brain transcriptome. There are 14 ribosomal protein genes (nine under expressed) among the 330 statistically significant differences between normal male and Ts65Dn male brains, which possibly implies abnormal ribosomal biogenesis in the development and maintenance of DS phenotypes. This study contributes to the establishment of a mouse brain transcriptome and provides the first overall analysis of the differences in gene expression in aneuploid versus normal mammalian brain cells.

Cloning, expression, and chromosome mapping of the murine Hip/Rpl29 gene

We previously have identified murine heparin/heparan sulfate-interacting protein (HIP) identical to mouse ribosomal protein L29 that is, like its human orthologue, distinctively expressed both on the cell surface and intracellularly in different adult tissues and cell types. In the present study, we show that mouse HIP/RPL29 is encoded by a single mRNA and that it is expressed to different extents in most of the tissues of the developing embryo without restriction to a specific cell type. We isolated the single-copy gene coding for murine Hip/Rpl29 among a large number of pseudogenes, established its structure, and assigned its location to distal chromosome 9. Similar to other ribosomal protein promoters, the promoter of Hip/Rpl29 is rich in polypyrimidine tracts, contains binding motifs for ubiquitously expressed transcription factors, and lacks a TATA box. Progressive 5' deletion analyses identified a strong enhancer region that includes CT-rich sequences and a potential consensus binding site for NF-kappaB. These data will provide valuable tools to progress the understanding of HIP/RPL29 function as a ribosomal protein and/or as a regulator of growth and cell adhesion through interaction with heparan sulfate proteoglycans.

Human and mouse homologues of the Drosophila melanogaster tweety (tty) gene: a novel gene family encoding predicted transmembrane proteins

We have cloned cDNA for TTYH1, a human homologue of the Drosophila melanogaster tweety (tty) gene. The 450-residue predicted protein shows 27% amino acid sequence identity (51% similarity) to the Drosophila protein, which contains an additional C-terminal repetitive region. A second Drosophila homologue exhibits 42% identity (65% similarity) to the tty protein. Mouse (Ttyh1), macaque, and Caenorhabditis elegans homologues were also identified, and the complete coding sequence for the mouse gene was determined. The mouse protein is 91% identical to the human protein. Hydrophobicity analysis of the tty-related proteins indicates that they represent a new family of membrane proteins with five potential membrane-spanning regions. The yeast FTR1 and FTH1 iron transporter proteins and the mammalian neurotensin receptors 1 and 2 have a similar hydrophobicity profile, although there is no detectable sequence homology to the tty-related proteins. This suggests that the tweety-related proteins could be involved in transport of iron or other divalent cations or alternatively that they may be membrane-bound receptors. TTYH1 was mapped to chromosome 19q13.4 by FISH and by radiation hybrid mapping using the Stanford G3 panel.

Structure of the human ubiquitin fusion gene Uba80 (RPS27a) and one of its pseudogenes

Ubiquitin is a highly conserved 76 amino acid protein that is generated in the cell by proteolysis of larger proteins containing either polyubiquitin chains or ubiquitin fused to carboxyl extension proteins (CEPs). In humans, the two human ubiquitin-CEP genes are Uba80 and Uba52, which code for ubiquitin fused to ribosomal protein S27a and L40, respectively. Working from a recently generated physical map of human chromosome 2p16, we determined the genetic and physical location and the genomic structure of the Uba80 gene in its entirety. A comparison of Uba80 to Uba52 revealed that the two genes share a conserved 5'-end structure, but that the structure of the ubiquitin coding regions was not conserved. Analysis of 400 bp of the promoter of Uba80 revealed strong similarity not only to the Uba52 promoter, but also to the other known human ribosomal gene promoters that have been identified to date. Homology searches also detected the presence of a pseudogene for Uba80, and the structure of this sequence feature is also reported.

Comparative organization of cattle chromosome 5 revealed by comparative mapping by annotation and sequence similarity and radiation hybrid mapping

A whole genome cattle-hamster radiation hybrid cell panel was used to construct a map of 54 markers located on bovine chromosome 5 (BTA5). Of the 54 markers, 34 are microsatellites selected from the cattle linkage map and 20 are genes. Among the 20 mapped genes, 10 are new assignments that were made by using the comparative mapping by annotation and sequence similarity strategy. A LOD-3 radiation hybrid framework map consisting of 21 markers was constructed. The relatively low retention frequency of markers on this chromosome (19%) prevented unambiguous ordering of the other 33 markers. The length of the map is 398.7 cR, corresponding to a ratio of approximately 2.8 cR(5,000)/cM. Type I genes were binned for comparison of gene order among cattle, humans, and mice. Multiple internal rearrangements within conserved syntenic groups were apparent upon comparison of gene order on BTA5 and HSA12 and HSA22. A similarly high number of rearrangements were observed between BTA5 and MMU6, MMU10, and MMU15. The detailed comparative map of BTA5 should facilitate identification of genes affecting economically important traits that have been mapped to this chromosome and should contribute to our understanding of mammalian chromosome evolution.

Ribosomal proteins in cell proliferation and apoptosis

Ribosomal proteins have the complex task of coordinating protein biosynthesis to maintain cell homeostasis and survival. Recent evidence suggests that a number of ribosomal proteins have secondary functions independent of their involvement in protein biosynthesis. A number of these proteins function as cell proliferation regulators and in some instances as inducers of cell death. Specifically, expression of human ribosomal protein L13a has been shown to induce apoptosis, presumably by arresting cell growth in the G2/M phase of the cell cycle. In addition, inhibition of expression of L13a induces apoptosis in target cells, suggesting that this protein is necessary for cell survival. Similar results have been obtained in the yeast Saccharomyces cerevisiae, where inactivation of the yeast homologues of L13a, rp22 and rp23, by homologous recombination results in severe growth retardation and death. In addition, a closely related ribosomal protein, L7, arrests cells in G1 and also induces apoptosis. Thus, it appears that a group of ribosomal proteins may function as cell cycle checkpoints and compose a new family of cell proliferation regulators.

Structure of the highly conserved HERC2 gene and of multiple partially duplicated paralogs in human

Recombination between chromosome-specific low-copy repeats (duplicons) is an underlying mechanism for several genetic disorders. Recently, a chromosome 15 duplicon was discovered in the common breakpoint regions of Prader-Willi and Angelman syndrome deletions. We identified previously the large HERC2 transcript as an ancestral gene in this duplicon, with approximately 11 HERC2-containing duplicons, and demonstrated that recessive mutations in mouse Herc2 lead to a developmental syndrome, juvenile development and fertility 2 (jdf2). We have now constructed and sequenced a genomic contig of HERC2, revealing a total of 93 exons spanning approximately 250 kb and a CpG island promoter. A processed ribosomal protein L41 pseudogene occurs in intron 2 of HERC2, and putative VNTRs occur in intron 70 (28 copies, approximately 76-bp repeat) and 3' exon 40 through intron 40 (6 copies, approximately 62-bp repeat). Sequence comparisons show that HERC2-containing duplicons have undergone several deletion, inversion, and dispersion events to form complex duplicons in 15q11, 15q13, and 16p11. To further understand the developmental role of HERC2, a highly conserved Drosophila ortholog was characterized, with 70% amino acid sequence identity to human HERC2 over the carboxy-terminal 743 residues. Combined, these studies provide significant insights into the structure of complex duplicons and into the evolutionary pathways of formation, dispersal, and genomic instability of duplicons. Our results establish that some genes not only have a protein coding function but can also play a structural role in the genome.

The intron-containing L3 ribosomal protein gene (RPL3): sequence analysis and identification of U43 and of two novel intronic small nucleolar RNAs

Isolation and sequencing of bovine and human intron-containing L3 ribosomal protein genes are here reported. They exhibit very similar organisation, both comprising 10 exons and nine introns. A polymorphic locus, involving a 19-bp deletion, was found in intron 6 of the human gene. The frequency of the two alleles has been estimated in 200 haploid genomes. In bovine and human genes intron sequences are rather different, except for limited regions, located in corresponding positions, which show a surprisingly high degree of identity. All these regions contain conserved features defining the box C/D class of small nucleolar RNAs. Demonstration is given that U43 small nucleolar RNA is encoded within the first intron of both bovine and human genes. Single nucleotide sequences, encoding two novel species of small nucleolar RNAs (U82, U83a and U83b), are located in introns 3, 5 and 7. Their expression has been investigated and a possible role of these molecules in 2'-O-ribose methylation of rRNAs is discussed.

Estimation of synteny conservation and genome compaction between pufferfish (Fugu) and human

BACKGROUND

Knowledge of the amount of gene order and synteny conservation between two species gives insights to the extent and mechanisms of divergence. The vertebrate Fugu rubripes (pufferfish) has a small genome with little repetitive sequence which makes it attractive as a model genome. Genome compaction and synteny conservation between human and Fugu were studied using data from public databases.

METHODS

Intron length and map positions of human and Fugu orthologues were compared to analyse relative genome compaction and synteny conservation respectively. The divergence of these two genomes by genome rearrangement was simulated and the results were compared to the real data.

RESULTS

Analysis of 199 introns in 22 orthologous genes showed an eight-fold average size reduction in Fugu, consistent with the ratio of total genome sizes. There was no consistent pattern relating the size reduction in individual introns or genes to gene base composition in either species. For genes that are neighbours in Fugu (genes from the same cosmid or GenBank entry), 40-50% have conserved synteny with a human chromosome. This figure may be underestimated by as much as two-fold, due to problems caused by incomplete human genome sequence data and the existence of dispersed gene families. Some genes that are neighbours in Fugu have human orthologues that are several megabases and tens of genes apart. This is probably caused by small inversions or other intrachromosomal rearrangements.

CONCLUSIONS

Comparison of observed data to computer simulations suggests that 4000-16 000 chromosomal rearrangements have occurred since Fugu and human shared a common ancestor, implying a faster rate of rearrangement than seen in human/mouse comparisons.