Complete Sequence Construction of the Highly Repetitive Ribosomal RNA Gene Repeats in Eukaryotes Using Whole Genome Sequence Data

A new method for determining ribosomal DNA copy number shows differences between Saccharomyces cerevisiae populations

Ribosomal DNA genes (rDNA) encode the major ribosomal RNAs and in eukaryotes typically form tandem repeat arrays. Species have characteristic rDNA copy numbers, but there is substantial intra-species variation in copy number that results from frequent rDNA recombination. Copy number differences can have phenotypic consequences, however difficulties in quantifying copy number mean we lack a comprehensive understanding of how copy number evolves and the consequences. Here we present a genomic sequence read approach to estimate rDNA copy number based on modal coverage to help overcome limitations with existing mean coverage-based approaches. We validated our method using Saccharomyces cerevisiae strains with known rDNA copy numbers. Application of our pipeline to a global sample of S. cerevisiae isolates showed that different populations have different rDNA copy numbers. Our results demonstrate the utility of the modal coverage method, and highlight the high level of rDNA copy number variation within and between populations.

Variability of Human rDNA

In human cells, ribosomal DNA (rDNA) is arranged in ten clusters of multiple tandem repeats. Each repeat is usually described as consisting of two parts: the 13 kb long ribosomal part, containing three genes coding for 18S, 5.8S and 28S RNAs of the ribosomal particles, and the 30 kb long intergenic spacer (IGS). However, this standard scheme is, amazingly, often altered as a result of the peculiar instability of the locus, so that the sequence of each repeat and the number of the repeats in each cluster are highly variable. In the present review, we discuss the causes and types of human rDNA instability, the methods of its detection, its distribution within the locus, the ways in which it is prevented or reversed, and its biological significance. The data of the literature suggest that the variability of the rDNA is not only a potential cause of pathology, but also an important, though still poorly understood, aspect of the normal cell physiology.

Structure of the intergenic spacers in chicken ribosomal DNA

BACKGROUND

Ribosomal DNA (rDNA) repeats are situated in the nucleolus organizer regions (NOR) of chromosomes and transcribed into rRNA for ribosome biogenesis. Thus, they are an essential component of eukaryotic genomes. rDNA repeat units consist of rRNA gene clusters that are transcribed into single pre-rRNA molecules, each separated by intergenic spacers (IGS) that contain regulatory elements for rRNA gene cluster transcription. Because of their high repeat content, rDNA sequences are usually absent from genome assemblies. In this work, we used the long-read sequencing technology to describe the chicken IGS and fill the knowledge gap on rDNA sequences of one of the key domesticated animals.

METHODS

We used the long-read PacBio RSII technique to sequence the BAC clone WAG137G04 (Wageningen BAC library) known to contain chicken NOR elements and the HGAP workflow software suit to assemble the PacBio RSII reads. Whole-genome sequence contigs homologous to the chicken rDNA repetitive unit were identified based on the Gallus_gallus-5.0 assembly with BLAST. We used the Geneious 9.0.5 and Mega software, maximum likelihood method and Chickspress project for sequence evolution analysis, phylogenetic tree construction and analysis of the raw transcriptome data.

RESULTS

Three complete IGS sequences in the White Leghorn chicken genome and one IGS sequence in the red junglefowl contig AADN04001305.1 (Gallus_gallus-5.0) were detected. They had various lengths and contained three groups of tandem repeats (some of them being very GC rich) that form highly organized arrays. Initiation and termination sites of rDNA transcription were located within small and large unique regions (SUR and LUR), respectively. No functionally significant sites were detected within the tandem repeat sequences.

CONCLUSIONS

Due to the highly organized GC-rich repeats, the structure of the chicken IGS differs from that of IGS in human, apes, Xenopus or fish rDNA. However, the chicken IGS shares some molecular organization features with that of the turtles, which are other representatives of the Sauropsida clade that includes birds and reptiles. Our current results on the structure of chicken IGS together with the previously reported ribosomal gene cluster sequence provide sufficient data to consider that the complete chicken rDNA sequence is assembled with confidence in terms of molecular DNA organization.

Rapid and Sensitive Detection of Repetitive Nucleic Acid Sequences Using Magnetically Modulated Biosensors

Repetitive DNA sequences are abundant in the genome of most biological species. These sequences are naturally "preamplified", which makes them a preferential target for a variety of biological assays. Current methods to detect specific DNA sequences are based on the quantitative polymerase chain reaction (PCR), which relies on target amplification by Taq polymerase and uses a fluorescent resonance energy transfer (FRET)-based probe. Here, to rapidly detect a repetitive DNA sequence, we combine a highly sensitive magnetic modulation biosensing (MMB) system and a modified double-quenched FRET-based probe. The high numbers of copies of the female-specific XhoI sequence of the domestic chicken (Gallus gallus), combined with the low background fluorescence levels of the modified double-quenched probe and the high sensitivity of the MMB system, allow us to determine the chick sex in ovo within 13 min, with 100% sensitivity and specificity. Compared to quantitative PCR, the presented assay is 4-9 times faster. More broadly, by specifically tailoring the primers and probe, the proposed assay can detect any target DNA sequence, either repetitive or nonrepetitive.

Ribosomal DNA-connecting ribosome biogenesis and chromosome biology

Ribosomal DNA, the topic of this special issue, has long fascinated biologists. The RNA products of the ribosomal DNA are the ribosomal RNAs that are part of the ribosome. In this special issue, we focus on the sequence, molecular organization, repair, stability, copy number, and peculiar genetics of this region of the genome. The locus can impact not only the translational capability of cells, but also genome organization, stability and integrity, providing a link between translation and chromosome biology.

The conservation landscape of the human ribosomal RNA gene repeats

Ribosomal RNA gene repeats (rDNA) encode ribosomal RNA, a major component of ribosomes. Ribosome biogenesis is central to cellular metabolic regulation, and several diseases are associated with rDNA dysfunction, notably cancer, However, its highly repetitive nature has severely limited characterization of the elements responsible for rDNA function. Here we make use of phylogenetic footprinting to provide a comprehensive list of novel, potentially functional elements in the human rDNA. Complete rDNA sequences for six non-human primate species were constructed using de novo whole genome assemblies. These new sequences were used to determine the conservation profile of the human rDNA, revealing 49 conserved regions in the rDNA intergenic spacer (IGS). To provide insights into the potential roles of these conserved regions, the conservation profile was integrated with functional genomics datasets. We find two major zones that contain conserved elements characterised by enrichment of transcription-associated chromatin factors, and transcription. Conservation of some IGS transcripts in the apes underpins the potential functional significance of these transcripts and the elements controlling their expression. Our results characterize the conservation landscape of the human IGS and suggest that noncoding transcription and chromatin elements are conserved and important features of this unique genomic region.

Dosage effects of human ribosomal genes (rDNA) in health and disease

Human ribosomal RNA genes encoding a pre-transcript of the three major ribosomal RNA (18S, 5.8S, and 28S rRNA) are tandemly repeated in human genome. Their total copy number varies from 250 to 670 per diploid genome with a mean of approximately 420 copies, but only a fraction of them is transcriptionally active. The functional consequences of human ribosomal RNA gene dosage are not widely known and often assumed to be negligible. Here, we review the facts of rRNA gene dosage effects on normal growth and aging, stress resistance of healthy individuals, and survivability of patients with chromosomal abnormalities, as well as on the risk and severity of some multifactorial diseases with proven genetic predisposition. An original hypothesis that rRNA gene dosage can be a modulating factor involved in the pathogenesis of schizophrenia and rheumatoid arthritis is put forward.

Variation in human chromosome 21 ribosomal RNA genes characterized by TAR cloning and long-read sequencing

Despite the key role of the human ribosome in protein biosynthesis, little is known about the extent of sequence variation in ribosomal DNA (rDNA) or its pre-rRNA and rRNA products. We recovered ribosomal DNA segments from a single human chromosome 21 using transformation-associated recombination (TAR) cloning in yeast. Accurate long-read sequencing of 13 isolates covering ∼0.82 Mb of the chromosome 21 rDNA complement revealed substantial variation among tandem repeat rDNA copies, several palindromic structures and potential errors in the previous reference sequence. These clones revealed 101 variant positions in the 45S transcription unit and 235 in the intergenic spacer sequence. Approximately 60% of the 45S variants were confirmed in independent whole-genome or RNA-seq data, with 47 of these further observed in mature 18S/28S rRNA sequences. TAR cloning and long-read sequencing enabled the accurate reconstruction of multiple rDNA units and a new, high-quality 44 838 bp rDNA reference sequence, which we have annotated with variants detected from chromosome 21 of a single individual. The large number of variants observed reveal heterogeneity in human rDNA, opening up the possibility of corresponding variations in ribosome dynamics.

Epichloë hybrida, sp. nov., an emerging model system for investigating fungal allopolyploidy

Endophytes of the genus Epichloë (Clavicipitaceae, Ascomycota) frequently occur within cool-season grasses and form interactions with their hosts that range from mutualistic to antagonistic. Many Epichloë species have arisen via interspecific hybridization, resulting in species with two or three subgenomes that retain all or nearly all of their original parental genomes, a process termed allopolyploidization. Here, we characterize Epichloë hybrida, sp. nov., a mutualistic species that has increasingly become a model system for investigating allopolyploidy in fungi. The Epichloë species so far identified as the closest known relatives of the two progenitors of E. hybrida are E. festucae var. lolii and E. typhina. We confirm that the nuclear genome of E. hybrida contains two homeologs of most protein-coding genes from E. festucae and E. typhina, with genome-wide gene expression analysis indicating a slight bias in overall gene expression from the E. typhina subgenome. Mitochondrial DNA is detectable only from E. festucae, whereas ribosomal DNA is detectable only from E. typhina. Inheriting ribosomal DNA from just one parent might be expected to preferentially favor interactions with ribosomal proteins from the same parent, but we find that ribosomal protein genes from both parental subgenomes are nearly all expressed equally in E. hybrida. Finally, we provide a comprehensive set of resources for this model system that are intended to facilitate further study of fungal hybridization by other researchers.