![[Feedback]](/icons/banner/feedbackACT.gif){kind=icon}
![[For Subscribers]](/icons/banner/subscriberACT.gif){kind=icon}
![[Archive]](/icons/banner/archiveACT.gif){kind=icon}
![[Search]](/icons/banner/searchACT.gif){kind=icon}
![[Contents]](/icons/banner/tocACT.gif){kind=icon}
|
|
|
|
|
||
|
Published online before print
November 19, 2007, 10.1101/gr.6858507 Genome Res. 18:13-18, 2008 ©2008 by Cold Spring Harbor Laboratory Press; ISSN 1088-9051/08 $5.00
Letter Genomic architecture and inheritance of human ribosomal RNA gene clusters1 Graduate Center for Toxicology, University of Kentucky, Lexington, Kentucky 40515, USA; 2 Department of Microbiology, Immunology and Molecular Genetics, University of Kentucky, Lexington, Kentucky 40515, USA; 3 Department of Internal Medicine, Markey Cancer Center, University of Kentucky, Lexington, Kentucky 40515, USA; 4 Department of Microbiology, Immunology and Molecular Genetics, Markey Cancer Center, University of Kentucky, Lexington, Kentucky 40515, USA
The finishing of the Human Genome Project largely completed the detailing of human euchromatic sequences; however, the most highly repetitive regions of the genome still could not be assembled. The 12 gene clusters producing the structural RNA components of the ribosome are critically important for cellular viability, yet fall into this unassembled region of the Human Genome Project. To determine the extent of human variation in ribosomal RNA gene content (rDNA) and patterns of rDNA cluster inheritance, we have determined the physical lengths of the rDNA clusters in peripheral blood white cells of healthy human volunteers. The cluster lengths exhibit striking variability between and within human individuals, ranging from 50 kb to >6 Mb, manifest essentially complete heterozygosity, and provide each person with their own unique rDNA electrophoretic karyotype. Analysis of these rDNA fingerprints in multigenerational human families demonstrates that the rDNA clusters are subject to meiotic rearrangement at a frequency >10% per cluster, per meiosis. With this high intrinsic recombinational instability, the rDNA clusters may serve as a unique paradigm of potential human genomic plasticity.
Each human ribosome contains four RNA molecules (rRNA) that play indispensable roles in all protein translation. In order to produce sufficient rRNA for the highly abundant ribosomes, the genes encoding the rRNA are represented genomically in multiple copies. The 5S rRNA molecule (OMIM 180420 [OMIM] ) is produced from a cluster of repeated 2.2-kb genes (Little and Braaten 1989  ; Sorensen and Frederiksen 1991) on maternal and paternal Chromosome 1q42 (Sorensen et al. 1991). Similarly, the remaining three rRNA molecules—18S, 5.8S, and 28S—are produced by post-transcriptional processing of a 45S precursor transcript expressed from clusters of repeated 43-kb genes (Gonzalez and Sylvester 1995) on maternal and paternal 13p12, 14p12, 15p12, 21p12, and 22p12 (OMIM 180450
[OMIM]
–180454) (Henderson et al. 1972). Overall, the rDNA clusters consist of repeated individual rRNA genes abutting each other in a largely head-to-tail orientation, devoid of intervening non-rDNA sequences (Little and Braaten 1989; Sorensen and Frederiksen 1991; Shiels et al. 1997; Schofer et al. 1998; Caburet et al. 2005).
The repetitive, clustered nature of the rDNA was noted in early quantitative hybridization studies (Schmickel and Knoller 1977
Our experimental approach is shown in Figure 1. Very high molecular weight genomic DNA is prepared from human peripheral blood (Birren and Lai 1993 ), and intact rDNA clusters are released from the bulk genomic DNA by restriction digestion. Since the individual rRNA genes are very highly self-similar (Sorensen et al. 1990; Gonzalez and Sylvester 2001), liberating intact rDNA gene clusters is feasible, as a restriction enzyme without a recognition site in a single rDNA gene will typically lack a recognition site throughout the entire cluster. The intact clusters are separated by size in a pulsed field gel and detected by hybridization with either a 5S or 45S rDNA-specific radiolabeled probe.
5S rDNA clusters (Fig. 2A) isolated from anonymous human blood donors demonstrate tremendous allelic diversity, with most of the donors showing two distinct 5S rDNA cluster lengths that clearly indicate differentially sized 5S rDNA clusters between maternal and paternal 1q42. Indeed, we have yet to observe a person homozygous for 5S rDNA cluster lengths. Occasionally, as seen in donor #10, we resolve three 5S rDNA bands; the presence of an extra band is enzyme-specific, and we infer that a chance polymorphism in one of the 5S repeats has coincidentally created this specific restriction site within the gene cluster.
The Human Genome Project (build 36.2) (International Human Genome Sequencing Consortium 2004 ) shows 17.2 5S rDNA repeats on 1q42 with an overall sequence conservation of 99.6% (Clamp et al. 2004). The Genome Project and other sources (Little and Braaten 1989; Sorensen and Frederiksen 1991; Timofeeva et al. 1993) give a unit 5S rDNA repeat length of 2236 bp. Using our 5S rDNA cluster length data and the lengths of the unit 5S rDNA repeat and non-rDNA flanking sequences (International Human Genome Sequencing Consortium 2004), we directly calculated the number of human 5S rDNA repeats (Table 1) in a total of 27 human donors (18 shown in Fig. 2A and nine additional unrelated individuals from Fig. 3A). Our average 5S rDNA cluster length of 98 repeats is consistent with earlier hybridization studies (Sorensen and Frederiksen 1991); however, the Human Genome Project’s total of 17.2 repeats is >2.3 standard deviations outside this value, likely as a result of the difficulty the Genome Project has in correctly assembling very highly self-similar sequences (She et al. 2004).
45S rDNA cluster lengths from the same anonymous donors (Fig. 2B) also exhibit striking person-to-person variability. Consistent with 45S clusters on the five pairs of human acrocentric chromosomes, our pulsed field gel analysis revealed eight to 10 bands per person. We assume that individuals showing fewer than 10 bands have yet additional bands >10 Mb that are incapable of entering a pulsed field gel (Birren and Lai 1993 ), although complete loss of rDNA sequences on specific acrocentric chromosomes is an alternative possibility. The 45S rDNA cluster size variation is consistent with our 5S rDNA cluster data: the 45S rDNA clusters display distinct and independent size variation on parentally homologous chromosomes. The 45S rDNA clusters shown here range in length from  70 kb to >6 Mb. For a unit 45S rDNA repeat size of 43 kb (Gonzalez and Sylvester 1995), these data indicate cluster sizes from a single 45S rDNA unit alone up to >140 repeats. While the lack of resolution above 6 Mb precludes an exact determination of human total diploid 45S rDNA content, our results are consistent with the upper range of earlier hybridization experiments of 600 copies per diploid human cell (Schmickel 1973). In contrast to the 5S rDNA clusters on 1q42, there are currently no data from the Human Genome Project concerning 45S rDNA cluster lengths. Beyond variation in cluster length between individuals, we observed other intriguing features of 5S and 45S rDNA architecture, including evidence of significant meiotic rearrangement. In our analysis of the 5S rDNA clusters from four different human families (Fig. 3A), it can be seen that the paternal 5S rDNA cluster of PYFI-02B was inherited with a meiotic rearrangement. We verified the presence of this de novo 5S rDNA cluster length, relative to the 5S clusters in PYFI-02B’s parents PYFI-10 and PYFI-20, with a panel of three different restriction enzymes (Fig. 3B). Commercial short tandem repeat (STR) testing verified that PYFI-20 is, indeed, the biological father of PYFI-02B (see Supplemental Material).
The higher resolution of the 5S rDNA clusters in Figure 3B also revealed an additional faint band in PYFI-02B of intermediate length at
The 45S rDNA clusters display an even larger degree of both meiotic rearrangement and somatic mosaicism (Fig. 4A). There are clear examples of inherited rearrangement in five of the offspring (ARDO-01, PYFI-02B, VAGR-02A, VAGR-02B, and THAE-02), as indicated by the bands enclosed in dotted circles. For example, VAGR-02B possesses a cluster of 4.6 Mb that is 1.6 Mb larger than the largest cluster from either his mother (VAGR-10) or his father (VAGR-20). A similar expansion is observed in the 6-Mb cluster of PYFI-02B that is at least 0.5 Mb larger than the largest rDNA cluster in either of his parents. Since clusters of repeated genes are potentially susceptible to recombinational erosion by single-strand annealing (Stark et al. 2004
A high-resolution gel for the 45S rDNA clusters <1 Mb in size (Fig. 4B) provides a superior view of this meiotic rearrangement. VAGR-02A has clusters of 625 kb and 715 kb (dotted circles) that are not found in either his mother (VAGR-10) or his father (VAGR-20). Likewise, the 950-kb cluster found in THAE-02 is not found in either of his parents (THAE-10 and THAE-20). Collectively, we observe 11 clear examples of 5S or 45S rDNA cluster meiotic restructuring in five of the eight parent–child trios analyzed.
Similar to the 5S rDNA cluster analysis of PYFI-02B (Fig. 3B), minor bands arising from 45S rDNA cluster somatic mosaicism are also apparent (Fig. 4B). PYFI-10 shows the strongest minor banding with a 910-kb length cluster at
Collectively, we observed one occurrence of 5S and 10 occurrences of 45S rDNA gene cluster meiotic restructuring in five of eight parent–child trios analyzed, for a total of 11 clearly identifiable events. Given a potential total of 96 different rDNA clusters (eight children, each with two 5S and up to 10 45S rDNA clusters), this yields an average recombination rate of 11% per generation per gene cluster (95% confidence interval: 6.4%–20%) (Agresti and Coull 1998 ). This rapid meiotic restructuring likely accounts for the present-day interperson variability in both rDNA quantity and genomic organization. We cannot currently ascertain the relative degree to which this diversity is generated by recombination between rDNA clusters on nonhomologous chromosomes, parental homologs, or sister chromatids, although familial analysis of flow-sorted acrocentric chromosomes may be able to answer this question.
Molecular combing experiments indicate that the rDNA clusters largely consist of repeats oriented in a head-to-tail configuration, but also with a substantial fraction of individual repeats within any given rDNA cluster in the inverted or palindromic orientation (Caburet et al. 2005
The human 5S and 45S rDNA gene clusters are clearly capable of very rapid copy number evolution. Although there may be considerable flexibility in rRNA expression at the transcriptional level via epigenetic control (Kawasaki et al. 1992 We anticipate that further study of the rDNA gene clusters, which are fundamentally required for cellular metabolism yet are highly variable in copy number and exhibit significant instability, will provide further insight into the role of copy number variation in human genomic architecture and stability.
High molecular weight genomic DNA was prepared in agarose from peripheral blood via ammonium chloride erythrocyte lysis (generally according to methods described in Birren and Lai 1993 ) with a final concentration of 1 x 107 white blood cells per milliliter of 0.8% low-melting temperature agarose. Agarose/DNA slices of 15 µL volume were digested with 50 units of indicated restriction enzymes (New England Biolabs) overnight at 37°C and subsequently sealed with low melting temperature agarose into gels for electrophoresis. All electrophoresis was performed using a CHEF Mapper (Bio-Rad) system. Gels were subsequently processed according to Lueders and Fewell (1994) and Ehtesham and Hasnain (1991) with minor modifications.
Hybridization probes
45S rDNA PCR products were subsequently digested with XhoI and XbaI, subcloned individually into pBluescript II SK– (Stratagene), and verified by DNA sequencing. XhoI/XbaI fragments from these plasmids were 32P radiolabeled using a Prime-It II kit (Stratagene).
Pulsed field gels The conditions in Figures 2B and 4A were genomic DNA digested with EcoRV; 1% SeaKem Gold agarose (Cambrex); 1x TAE (40 mM Tris, 40 mM acetate, 1 mM EDTA at pH 8.0); 14°C; 2 V/cm; 106° included angle; 5–40-min field switching with linear ramp; and 92 h of total electrophoresis. The conditions in Figure 3A were genomic DNA digested with HindIII; 1% PFC agarose (Bio-Rad); 0.5x TBE; 14°C; 6 V/cm; 120° included angle; 10–40-sec field switching with linear ramp; and 20 h of total electrophoresis.
The conditions in Figure 3B were 1% PFC agarose; 0.5x TBE; 14°C; 6 V/cm; sequence of electric field vectors relative to 0° as "forward": –60°, 180°, +60°, –60°, +60°, 180°, –60°, +60°; field durations: ±60° field vectors: 10 sec; 180° field vectors: 5 sec, 42 h of total electrophoresis. Pulse field parameters were adapted from Clark et al. (1988). The conditions in Figure 4B were as in Figure 2A, except we used a "ramp factor" (Bio-Rad) of 0.357.
Gel post-electrophoresis processing and hybridization
Commercial STR allele testing
Human subjects
We thank our human families and the Kentucky Blood Center for providing blood. Funding for these studies was provided by the Markey Cancer Center and aided by an Institutional Research Grant from the American Cancer Society (A.J.P.). D.M.S. performed the experiments in Figures 2A and 3. M.W.K. helped prepare genomic DNA from blood. H.H.P. recruited and consented human subjects. A.J.P. performed the experiments in Figures 2B and 4. D.M.S. and A.J.P. wrote the paper. All authors discussed the results and commented on the manuscript.
5 Corresponding author.
E-mail andrew.pierce{at}uky.edu; fax (859) 257-8940. [Supplemental material is available online at www.genome.org.] Article published online before print. Article and publication date are at http://www.genome.org/cgi/doi/10.1101/gr.6858507
Agresti, A. and Coull, B.A. 1998. Approximate is better than "exact" for interval estimation of binomial proportions. Am. Stat. 2: 119–125. Antonarakis, S.E., Adelsberger, P.A., Petersen, M.B., Binkert, F., and Schinzel, A.A. 1990. Analysis of DNA polymorphisms suggests that most de novo dup(21q) chromosomes in patients with Down syndrome are isochromosomes and not translocations. Am. J. Hum. Genet. 47: 968–972.[Medline] Arnheim, N., Krystal, M., Schmickel, R., Wilson, G., Ryder, O., and Zimmer, E. 1980. Molecular evidence for genetic exchanges among ribosomal genes on nonhomologous chromosomes in man and apes. Proc. Natl. Acad. Sci. 77: 7323–7327. Birren, B. and Lai, E., eds. 1993. Pulsed field gel electrophoresis—A practical guide. Academic Press, San Diego, CA. Boyd, L.J., Livingston, J.S., Brown, M.G., Lawce, H.J., Gilhooly, J.T., Wildin, R.S., Linck, L.M., Magenis, R.E., and Pillers, D.A. 2005. Meiotic exchange event within the stalk region of an inverted chromosome 22 results in a recombinant chromosome with duplication of the distal long arm. Am. J. Med. Genet. A 138: 355–360.[Medline] Caburet, S., Conti, C., Schurra, C., Lebofsky, R., Edelstein, S.J., and Bensimon, A. 2005. Human ribosomal RNA gene arrays display a broad range of palindromic structures. Genome Res. 15: 1079–1085. Clamp, M., Cuff, J., Searle, S.M., and Barton, G.J. 2004. The Jalview Java alignment editor. Bioinformatics 20: 426–427. Clark, S.M., Lai, E., Birren, B.W., and Hood, L. 1988. A novel instrument for separating large DNA molecules with pulsed homogeneous electric fields. Science 241: 1203–1205. Delany, M.E., Muscarella, D.E., and Bloom, S.E. 1994. Effects of rRNA gene copy number and nucleolar variation on early development: Inhibition of gastrulation in rDNA-deficient chick embryos. J. Hered. 85: 211–217. Delany, M.E., Taylor, R.L., and Bloom, S.E. 1995. Teratogenic development in chicken embryos associated with a major deletion in the rRNA gene cluster. Dev. Growth Differ. 37: 403–412.[CrossRef] Draptchinskaia, N., Gustavsson, P., Andersson, B., Pettersson, M., Willig, T.N., Dianzani, I., Ball, S., Tchernia, G., Klar, J., Matsson, H., et al. 1999. The gene encoding ribosomal protein S19 is mutated in Diamond-Blackfan anaemia. Nat. Genet. 21: 169–175.[CrossRef][Medline] Ehtesham, N.Z. and Hasnain, S.E. 1991. Direct in-gel hybridization, without blotting, using nick-translated cloned DNA probe. Biotechniques 11: 718–722.[Medline] Ellis, S.R. and Massey, A.T. 2006. Diamond Blackfan anemia: A paradigm for a ribosome-based disease. Med. Hypotheses 66: 643–648.[CrossRef][Medline] Gazda, H.T., Grabowska, A., Merida-Long, L.B., Latawiec, E., Schneider, H.E., Lipton, J.M., Vlachos, A., Atsidaftos, E., Ball, S.E., Orfali, K.A., et al. 2006. Ribosomal protein S24 gene is mutated in Diamond-Blackfan anemia. Am. J. Hum. Genet. 79: 1110–1118.[CrossRef][Medline] Gonzalez, I.L. and Sylvester, J.E. 1995. Complete sequence of the 43-kb human ribosomal DNA repeat: Analysis of the intergenic spacer. Genomics 27: 320–328.[CrossRef][Medline] Gonzalez, I.L. and Sylvester, J.E. 2001. Human rDNA: Evolutionary patterns within the genes and tandem arrays derived from multiple chromosomes. Genomics 73: 255–263.[CrossRef][Medline] Guanti, G. and Petrinelli, P. 1974. rDNA and acrocentric chromosomes in man. Cell Differ. 2: 319–324.[CrossRef][Medline] Henderson, A.S., Warburton, D., and Atwood, K.C. 1972. Location of ribosomal DNA in the human chromosome complement. Proc. Natl. Acad. Sci. 69: 3394–3398. International Human Genome Sequencing Consortium. 2004. Finishing the euchromatic sequence of the human genome. Nature 431: 931–945.[CrossRef][Medline] Kawasaki, K., Minoshima, S., Kudoh, J., Fukuyama, R., and Shimizu, N. 1992. Methylation status of ribosomal RNA gene clusters in the flow-sorted human acrocentric chromosomes. Mamm. Genome 3: 173–178.[CrossRef][Medline] Kuick, R., Asakawa, J., Neel, J.V., Kodaira, M., Satoh, C., Thoraval, D., Gonzalez, I.L., and Hanash, S.M. 1996. Studies of the inheritance of human ribosomal DNA variants detected in two-dimensional separations of genomic restriction fragments. Genetics 144: 307– 316.[Abstract] Little, R.D. and Braaten, D.C. 1989. Genomic organization of human 5 S rDNA and sequence of one tandem repeat. Genomics 4: 376–383.[CrossRef][Medline] Lueders, K.K. and Fewell, J.W. 1994. Hybridization of DNA in dried gels provides increased sensitivity compared with hybridization to blots. Biotechniques 16: 66–67.[Medline] Markovic, V.D., Worton, R.G., and Berg, J.M. 1978. Evidence for the inheritance of silver-stained nucleolus organizer regions. Hum. Genet. 41: 181–187.[Medline] Mayer, C., Schmitz, K.M., Li, J., Grummt, I., and Santoro, R. 2006. Intergenic transcripts regulate the epigenetic state of rRNA genes. Mol. Cell 22: 351–361.[CrossRef][Medline] Miller, D.A., Dev, V.G., Tantravahi, R., and Miller, O.J. 1976. Suppression of human nucleolus organizer activity in mouse–human somatic hybrid cells. Exp. Cell Res. 101: 235–243.[CrossRef][Medline] Page, S.L. and Shaffer, L.G. 1998. Chromosome stability is maintained by short intercentromeric distance in functionally dicentric human Robertsonian translocations. Chromosome Res. 6: 115–122.[CrossRef][Medline] Page, S.L., Shin, J.C., Han, J.Y., Choo, K.H., and Shaffer, L.G. 1996. Breakpoint diversity illustrates distinct mechanisms for Robertsonian translocation formation. Hum. Mol. Genet. 5: 1279–1288. Schmickel, R.D. 1973. Quantitation of human ribosomal DNA: Hybridization of human DNA with ribosomal RNA for quantitation and fractionation. Pediatr. Res. 7: 5–12.[Medline] Schmickel, R.D. and Knoller, M. 1977. Characterization and localization of the human genes for ribosomal ribonucleic acid. Pediatr. Res. 11: 929–935.[Medline] Schofer, C., Weipoltshammer, K., Almeder, M., and Wachtler, F. 1998. Arrangement of individual human ribosomal DNA fragments on stretched DNA fibers. Histochem. Cell Biol. 110: 201–205.[CrossRef][Medline] She, X., Jiang, Z., Clark, R.A., Liu, G., Cheng, Z., Tuzun, E., Church, D.M., Sutton, G., Halpern, A.L., and Eichler, E.E. 2004. Shotgun sequence assembly and recent segmental duplications within the human genome. Nature 431: 927–930.[CrossRef][Medline] Shiels, C., Coutelle, C., and Huxley, C. 1997. Analysis of ribosomal and alphoid repetitive DNA by fiber-FISH. Cytogenet. Cell Genet. 76: 20–22.[Medline] Sorensen, P.D. and Frederiksen, S. 1991. Characterization of human 5S rRNA genes. Nucleic Acids Res. 19: 4147–4151. Sorensen, P.D., Simonsen, H., and Frederiksen, S. 1990. Nucleotide sequence of a human 5S rRNA gene. Nucleic Acids Res. 18: 3060. Sorensen, P.D., Lomholt, B., Frederiksen, S., and Tommerup, N. 1991. Fine mapping of human 5S rRNA genes to chromosome 1q42.11–q42.13. Cytogenet. Cell Genet. 57: 26–29.[Medline] Stark, J.M., Pierce, A.J., Oh, J., Pastink, A., and Jasin, M. 2004. Genetic steps of mammalian homologous repair with distinct mutagenic consequences. Mol. Cell. Biol. 24: 9305–9316. Sullivan, B.A., Jenkins, L.S., Karson, E.M., Leana-Cox, J., and Schwartz, S. 1996. Evidence for structural heterogeneity from molecular cytogenetic analysis of dicentric Robertsonian translocations. Am. J. Hum. Genet. 59: 167–175.[Medline] Taylor, E.F. and Martin-DeLeon, P.A. 1981. Familial silver staining patterns of human nucleolus organizer regions (NORs). Am. J. Hum. Genet. 33: 67–76.[Medline] Timofeeva, M., Kost, M.V., Tikhomirova, T.P., Sarafanov, A.G., Elbert, B.L., Alimov, A.A., Kvitsiniia, S.L., Beritashvili, D.R., and Zelenin, A.V. 1993. Organization of a 5S ribosomal RNA gene cluster in the human genome. Mol. Biol. 27: 861–868. Worton, R.G., Sutherland, J., Sylvester, J.E., Willard, H.F., Bodrug, S., Dube, I., Duff, C., Kean, V., Ray, P.N., and Schmickel, R.D. 1988. Human ribosomal RNA genes: Orientation of the tandem array and conservation of the 5' end. Science 239: 64–68.
Received June 28, 2007; accepted in revised format September 17, 2007.  CiteULike  Connotea  Del.icio.us  Digg  Reddit  Technorati What's this?
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Top
Abstract
Results
3 = 145,506 bp, 
