Prerpints.org logo
Preprint
Article

Diversity of Unusual Ribosomal Genes and Ecological Origin of Rice

Altmetrics

Downloads

124

Views

60

Comments

0

A peer-reviewed article of this preprint also exists.

supplementary.pdf (958.81KB )

Submitted:

27 December 2023

Posted:

27 December 2023

You are already at the latest version

Alerts
Abstract
Tandemly organized rRNA genes are the typical example of a multigene family, where individual members evolve coordinately within but independently between species due to gene conversion and unequal crossing over. More frequently, in eukaryotic species with an inter-species hybrid origin, expression of unhomogenized rRNA genes from one progenitor is epigenetically silenced because of nucleolus dominance, and distinct rRNA genes may lose functionality and evolve faster. Interestingly, we obtained unusual ribosomal gene sequences from Oryza species that showed great diversity and did not appear in the present rice genomic sequence. The diversity of rDNA sequences indicated that the homogenization in rice is incomplete and explains the introgression of distinct rRNA gene families into ancestral rice genome before speciation and continent separation. The divergent large subunit (LSU) ribosomal genes are expressed, some of them differentially, depending on the N-fertilization of plants. Detection of differential transcripts of the rRNA genes suggested that rRNA gene families are not functionally equivalent. Phylogenetic analysis assigned Oryza species branching order to monocots, and monocot lineages probably have the same ecological origin by molecular clock calculation. Therefore, our results suggested that the geographical distances of continents separation cause barriers to the gene flow and homogenization among Oryza species which requires further explanation,
Keywords: 
Subject: Biology and Life Sciences  -   Agricultural Science and Agronomy

1. Introduction

Genera of the Gramineae, Oryza, does not only provide food for half of the world's population but is also used as a model plant by the scientific community. The genus Oryza contains about 27 species but is remarkable in the diverse ecological adaptations of its species, in total 11 genome types; the diploid six (AA, BB, CC, EE, FF, and GG), and tetraploid five (BBCC, CCDD, HHJJ, HHKK, and KKLL)[1,2]. For thousands of years, Rice has grown uniquely in sustainable high-output agroecosystems. In relation to genera containing other cereals, Oryza occupies a distinct phylogenetic position in a separate subfamily, the Ehrhartoideae [3]. Species in the genus Oryza and other genera closely related to Oryza have been extensively studied either because of their agronomically useful traits in rice genetic improvement (wild species in Oryza and Porteresia) or because of their economic value as part of cuisine (Zizania) and forage (Leersia)[4]. The origin and diversification of this tribe, in particular the origin of the genus Oryza and its divergence, remains largely unclear [5].
Protein production occurs on Ribosomes in all living cells. Ribosomes exhibit substantial differences across species, although universally conserved[6,7]. Some variations are observed in ribosomal RNAs (rRNAs) which play the role of central interface for hundreds of proteins. These proteins are ribosome assembly factors (AFs) and ribosomal proteins (RPs) which are highly variable in terms of length and sequence in various species [8,9,10]. Interestingly, eukaryotic RNAs have species-specific expansion segments (ESs) in comparison to prokaryotic rRNAs which are the hot spot of variations among different eukaryotic species[11,12]. Furthermore, in prokaryotes, ES/ES-like segments were also reported which are not common in eukaryotic species even heterogeneous rRNAs with differential expression of sequences have also been observed in various eukaryotic species, e.g., Plasmodium, zebrafish, mice, and Homo sapiens[10,13].
The duplication of single genes or of whole genomes followed by functional gene diversification are important processes that lead to evolutionary innovation[14]. The large proportion of genes with no recognizable homolog in the genome of the diploid Oryza sativa spp. indica suggests that the rice genome had been dupli cated one or several times during evolution [15] and that Oryza sativa as vertebrates [16] and most other angiosperms [17] may be degenerate, ancient polyploids. The contribution of hybridization, autopolyploidy, or allopolyploidy to the origin of the rice genome could be not estimated yet, because rice genes show unique gradients in GC content, codon, and amino acid usage which confound genome analyses based on comparisons of protein-encoding regions[18]. However, in some cases, the contribution of interspecies hybridization to the origin of a particular species could be also traced back by inference of the evolutionary history of its conserved rRNA genes. The highly conserved genes of 16S-18S and 23S-28S RNA genes are usually contained head to tail in repeated tandem units and evolve in a concerted manner. Within most species, they have nearly identical sequences among individual copies due to homogenization processes, which allowed them to treat them as single-copy genes in a phylogenetic context. Follo wing interbreeding between two different species, xenologous rRNA genes from one parent are usually epigenetically silenced because of nucleolar dominance and probably deleted later. Alternatively, xenologous rRNA genes may be also converted to one sequence type [19], or as in Quercus robur and Quercus petraea, parental rDNA families may be even maintained in the nucleus[20].
Genetic distance is the measurement of evolutionary changes between sequences from different organisms and is calculated for a pair of sequences by simply counting the number of nucleotides or amino acids that differ between them (Morgenstern et al. 2015). A small genetic distance between two sequences may suggest a recent common ancestor but is also consistent with a slower rate of sequence change and a more ancient common ancestor. Evolutionary rates depend on a combination of factors: generation time, population size, metabolic rate, the efficacy of DNA repair, and the degree to which mutations are beneficial or deleterious, all of which may vary among species. Often paleontology or biogeography can provide a date for one or more points in a phylogeny, which are then used to “calibrate” the timescale for the rest of the phylogeny[21]. Remarkably, sample sequencing at different times is the estimation of evolution which requires a faster evolutionary rate and/or negligible sampling times [22]. Therefore, evolution requires the usual phylogeny assuming that all its branches evolve at the same rate following a static molecular clock of sequence change [23] calculating the origin time and fossils existence of the organisms.
We are interested in understanding the evolution of rRNA genes due to the uneven distribution of crossing over which evolves coordinately within but independently between species [24]. Here, we have undertaken the functional characterization of unusual ribosomal gene sequences from Oryza species that showed great diversity but are absent in the present genomic sequence. We also estimated the divergence between the rRNA families in rice. Moreover, we observed that rice rDNA homogenization is not complete and provides an explanation to the introgression of distinct rRNA gene families into ancestral rice before speciation. Furthermore, our study showed that the large subunit (LSU) of ribosomal genes of all families are expressed depending on N fertilization of plants.
The introduction should briefly place the study in a broad context and highlight why it is important. It should define the purpose of the work and its significance. The current state of the research field should be carefully reviewed and key publications cited. Please highlight controversial and diverging hypotheses when necessary. Finally, briefly mention the main aim of the work and highlight the principal conclusions. As far as possible, please keep the introduction comprehensible to scientists outside your particular field of research. References should be numbered in order of appearance and indicated by a numeral or numerals in square brackets—e.g., [1] or [2,3], or [4,5,6]. See the end of the document for further details on references.

2. Materials and Methods

DNA and RNA extraction
Oryza officinalis and 17 Oryza species, O. alta (ACC105143, with IRGC accession number), O. australiensis (ACC100882), O. glaberrima (TOG5674), O. grandiglumis (ACC105669), O. granulata (ACC102119), O. latifolia (ACC100190), O. longiglumis (ACC105148), O. longistaminata, O. malampuzhaensis (ACC105329), O. minuta (ACC101141), O. nivara (ACC105763), O. officinalis (ACC101399), O. punctata (ACC105690), O. rhizomatis (ACC105432), O. ridleyi (ACC100821), O. rufipogon (ACC106423), O. sativa subsp. indica (IR36), were selected for further analysis. Plants including three-way hybrid Zea mays were grown in unfertilized rice field soil. For some experiments, plants were fertilized with 2 g NH4NO3 per kg soil. DNA and rRNA preparations from ribosomes were extracted from fresh leaves or roots kept in liquid nitrogen according to standard methods [25].
Gene sequencing
28S rRNA genes corresponding to the positions 184 to 1101 in the Saccharomyces cerevisiae 28S rRNA gene were amplified by PCR with forward primer 28f1 (5´-GAC CCC AGG TCA GGC GGG ACT ACC-3´) and reverse primer 28r1 (5´-GCT ATC CTG AGG GAA ACT TCG GAG G-3´) and sequenced with the ALF Express (Pharmacia, Uppsala, Sweden (SI-1). For specific amplification of a partial RDF2 rRNA tandem unit we used primers NS7 and 98rev. PCR conditions were as follows: initial denaturation at 95°C for 2 min, 30 cycles with denaturation at 95°C for 1 min, annealing at 60°C for 1 min, and extension at 72°C for 1 min, with a final denaturation at 72°C for 10 minutes. PCR products were cloned into the TOPO vector (Cat. No. K4500-40, Invitrogen Co.) according to the manufacturer‘s instructions and sequenced with standard primers. Reverse transcription PCR (RT-PCR) on total RNA or ribosomal RNA were done with RT-PCR beads (Amersham Pharmacia Biotech, Freiburg, Germany) following the manufacturer‘s protocol. Reverse transcription was carried out with primer 25r1. Direct sequencing of cDNA PCR products was carried out with forward primer 283fCy5 (5´-GCA (AG) CC CAA AT (CT) (AG) GG (ACT) G (AG) TAA AC-3´) and reverse primer 406rCy5 (5´-CA (AC) GCA CT (GCT) TTT GAC TCT CTT TTC-3´). For specific amplification of cDNA products of RC1, RC2 and RC3 rDNA families, the following primer pairs for nested PCR were used 98f2 (5´-CGC ACC GTT CGA ACT GTA GTC-3´) and 98r2 (5´-GTT ACA GCG TGG CAC CCC AAG G-3´), 914f (5´-GCC CAA CGT GAA AAT CGG GCA G-3´) and 914r (5´-GTA TCA CTT TGA GCC TCC ACC-3´), 931f (5´-CCC TAA TAA CCG AAT TGT AGT CTG G-3´) and 931r (5´-CTA GAT GGT TCG ATT GGT CTC ATC C-3´), respectively (SI-1). PCR conditions were as given above. Additionally, for O. officinalis and for O. longistaminata genomic libraries were constructed in the lambda ZAP express vector (Stratagene, San Diego, CA) and screened with 5’-digoxigenin labeled primer 28f1. The obtained 94 sequences were deposited in GenBank under accession numbers AF363835 to AF363923, AF375550, AF375551, AY097331, and AY097327 to AY097329.
Probes (RC and RDF) labeled and northern hybridization
Three specific probes RCp1, RCp2, and RCp3 digested from 28S rDNA clones (AF363841, AF363902, and AF363882 representing RC1, RC2, and RC3, respectively) with Ahd I or Bsm I and resulting 225, 261 and 276 bp fragments, respectively. Furthermore, RDF1, RDF2, and RDF3 were digested each with Ahd I or Bsm I and resulting 225 to 276 bp fragments corresponding to nucleotides 7-232, 7-268, or 7-283 of equivalent PCR products. These all were labeled with digoxigenin by random primed labeling. Labeled fragments correspond to variable regions B10-B18 of the secondary structure model of the Saccharomyces cerevisiae 28S rRNA gene. The RDF2 18S probe is a 3’ and 5’ with digoxigenin end labeled primer and corresponds in length to the sequence (SI-2B). Probes were hybridized at high stringency and visualized by CDP-Star (Roche, Mannheim, Germany) on an X-ray film according to the manual’s guide.
Terminated restriction fragment length polymorphism (T-RFLP)
PCR products were obtained from the cDNA templates with primer 98f2 with 5'-labeled Cy5. PCR products purified after agarose gel electrophoresis were digested with restriction enzymes and analyzed on an ALF Express automatic sequencer (Pharmacia, Uppsala, Sweden).
Sequences data analysis
DNA sequences were aligned by the Ribosomal Database II 28S (http://rdp.cme.msu.edu/cgis/seq_align.cgi?su=LSU) and the MultAlin alignment tools (http://prodes.toulouse.inra.fr/multalin/multalin.html). The alignment producing the shortest tree was used for (equal weighting) maximum parsimony, neighbor-joining, and maximum likelihood analyses with PAUP version 4.0b2a [5]. Parsimony jackknifing (1000 replicates) with TBR branch swapping, 10 random sequence additions, and gaps treated as missing data were used to provide estimates for internal support of clades. Neighbor-joining trees were inferred from maximum likelihood distances and 1000 bootstrap re-samplings. Maximum likelihood trees (1000 replicates) were estimated by a heuristic search procedure using the GTR+G+I and a four-category discrete gamma model of among-site rate variation. Average distances were estimated with the software package PHYLTEST, version 2.0 [26]. The gamma distribution shape parameter alpha was estimated from the data set using PAUP version 4.0 [27].

3. Results

3.1. Phylogeny Inferred from Concatenated Sequences of the Genus Oryza

Recombinant DNA technology has become an extremely important tool in biological sciences for manipulating DNA fragments to obtain desired characteristics. The rDNA families reported here were obtained in geographically isolated and widely distributed wild and cultivated species of the genus Oryza (Figure 1). We sequenced 94 clones obtained from PCR amplifications, containing partial sequences of 28S rRNA genes of 17 Oryza species. Sequence alignment, restriction enzyme sites, and used primers are listed in supplementary information 1 (SI-1). Phylogenetic analyses depict the descendent evolutionary lines of different genes, organisms, or species from a common ancestor. Therefore, our phylogenetic analyses are based on 864 alignment positions with 723-845 nucleotides and are consistent with previous results on the phylogeny of monocots based on 18S rDNA, rbcL, and atpB genes [28]. Further, we analyzed the genealogical relationships within rice cluster I (RCI) by networking approach allowing the assessment of a putative evolutionary origin of the ribosomal sequences through estimation of root probabilities. This analysis identified two networks (A and B) in RCI. Network A comprises two ribosomal sequences of O. glaberrima and O. punctata from Africa (AF363874, AF363875) which are located on a single branch in the phylogenetic tree (Figure 1). Network B is represented by all other sequences which are located on the well supported sister branch. In this network the highest root probability (0.114) was estimated for sequence AF363842 from O. sativa. For sequences AF363848, AF363856, AF363857, AF363864, and AF363868 root probabilities of 0.057 and for the rest of 0.028 were calculated. This result is consistent with an important role of O. sativa in the trans continental distribution of the rRNA families by interspecific hybridization and gene flow even to species of wild rice like O. longiglumis and O. australiensis, which are not known to exchange genes with O. sativa. The possibility that O. sativa plays an im por tant role in gene flow between rice species does not contradict our conclusion that the rRNA families in rice have been originally brought together by an ancient hybridiza tion event predating speciation. Because in the genome of O. sativa RC2 is missing but is present in the genome of all other rice species subjected by us to genomic Southern hybridization, this LSU rRNA family like the other ones was probably present in wild species of rice previous to interspecific hybridization with O. sativa.
Topologically congruent trees were obtained from all three phylogeny reconstruction methods. Three major rDNA families were assigned as RC1, RC2, and RC3 (Figure 1). Gene fragments obtained were highly divergent from the published O. sativa sequence (M11585). Interestingly, most of the obtained 28S rRNA gene sequences were low homology to that of the O. sativa sequence (M11585), even if the genome sequences of the indica and japonica subspecies have been assembled by whole-genome shotgun sequencing. The assembled sequence covers 93% of the 420-megabase genome. Homologs of 98% of the sequenced maize, wheat, and barley proteins are found in rice [29] providing the base for the development of cereals.
Preprints 94623 g001
Figure 1. Maximum likelihood (ML) and maximum parsimony (MP) consensus tree (for groups receiving >50% support) inferred from partial 28S rRNA gene sequences showing phylogenetic relationships of Oryza species and other angiosperms. Nodes are labeled with jackknife support values from MP, ML, or bootstrap values from neighbor-joining analyses (NJ). The lowest support value >55% in all three analyses is shown A to R. Three well-supported 28S rRNA gene clusters are assigned as RC1-3. Genome types are given for each Oryza species. Dicot sequences were used as outgroup. This is a figure. Schemes follow the same formatting.
Figure 1. Maximum likelihood (ML) and maximum parsimony (MP) consensus tree (for groups receiving >50% support) inferred from partial 28S rRNA gene sequences showing phylogenetic relationships of Oryza species and other angiosperms. Nodes are labeled with jackknife support values from MP, ML, or bootstrap values from neighbor-joining analyses (NJ). The lowest support value >55% in all three analyses is shown A to R. Three well-supported 28S rRNA gene clusters are assigned as RC1-3. Genome types are given for each Oryza species. Dicot sequences were used as outgroup. This is a figure. Schemes follow the same formatting.
Preprints 94623 g001

3.2. Validation of the PCR Sequence Artifacts through Southern Hybridization

To exclude the possibility of PCR artifacts for the diverged rDNA families, genomic southern hybridization with specific probes for 3 different rDNA families (RC1, RC2, and RC3) revealed that the different rDNA families were not PCR artifacts caused by sequence contamination or PCR recombination [30] (Figure 2). In a blot of genomic DNA, BseRI- and XcmI-digested total DNAs from five rice species were hybridized with the three specific probes. BseRI does not digest RC2 but RC1 and RC3 at nucleotides 60-63 (corresponding to the positions of Saccharomyces cerevisiae 28S rRNA gene). XcmI digests three rDNA families RC1, RC2, and RC3 at nucleotides 861-864, 929-931, and 914-918 (corresponding to the positions of Saccharomyces cerevisiae 28S rRNA gene) respectively. Accordingly, genomic DNA of RC1 and RC3 digested with BseRI and XcmI generated 800 bp and 860 bp fragments, respectively, were exhibited by hybridization (Figure 2). In contrast, in the specific RCp2 probed blot only 4.5 kb and 10 kb fragments hybridized due to absence of the BseRI restriction sites. The other hybridizing fragments (RC2, lanes 5 and 6; RC3, lane 6,) indicated that the rDNA families in rice are polymorphic.

3.3. Fragments Hybridization of RDF1-4

In control blots, RDF1-3 18S probes only hybridized with 18S gene fragments from the same but not from a different rDNA family (SI-2B-D). In a blot of genomic DNA, Bse RI-and Xcm I-digested total DNAs from five rice species were hybridized with the three RDF 28S probes. Bse RI does not restrict RDF2- but cuts RDF1- and RDF3 28S PCR fragments at nucleotides 68-71. Xcm I cuts all retrieved RDF1-, RDF2-, and RDF3 PCR fragments. Accordingly, in RDF1 and RDF3 probed lanes 800 bp and 860 bp fragments, respectively, were hybridizing (SI-2C).

3.4. Analyses of 28S Genes in Ribosome of O. officinalis

After conforming the polymorphic behavior of these rDNA fragments, we further analyzed 28S rDNA which are conserved universal DNA barcodes for living organisms including plants (Thompson and Christensen, 2011). Northern analyses from ribosome preparations showed that 28S genes belonging to rDNA families RC1, RC2, and RC3 are co-dominantly expressed and present in the ribosome of O. officinalis (Figure 3). A further indication that 28S rDNA family genes in rice might be functional came from expression analyses. This is surprising because their A + T content is increased (0.26% O. sativa; 0.24% O. longistaminata; 0.24% O. officinalis) in comparison to the much less diverged rice LSU rRNAs (0.21%; M11585, AY097330, AY097328, AY097327, respectively) within the same genome. An increased A + T content in alleles usually indicates a relaxation in selection which can follow duplication of genes [31]. This kind of rRNA gene transcription would indicate that rDNA family genes were not silenced and may have remained functional because of positive selection. Apart from the common 5´-cap-structure-dependent initiation of translation, also complementarity of rRNA within internal ribosome entry sides (IRES) of particular mRNAs [32], or binary interactions of ribosomal subunits with mRNA in general may be important for controlling the levels of particular proteins in eukaryotes [33]. In the ancestor(s) of cereal plants rRNA genes may have been selected for, which maintained this type of translational control for certain genes. This model predicts that mRNAs of particular genes are preferentially translated by rDNA family, or rDNA family hybrid ribosomes in rice. The high heterogeneity of rice ribosomes will provide a unique opportunity to test the ribosome filter hypothesis [33] which predicts that ribosomes will translate individual mRNAs with different relative efficiencies.

3.5. Conservation of RNA Secondary Structure

RNA secondary structure is similar to an alignment of protein or nucleic acid sequences, except that the RNA sequence folds back on itself by a “complementary bases'' pair, and it is also much more conserved than the sequence. The implications of conserved secondary structure changed 1 base of a pair, its partner is also changed so as to conserve that base pair. Ribosomal functions such as peptidyl transferase activity and interaction sites for translation factors and tRNAs have been mapped to this phylogenetically conserved structure. Parts of the LSU rRNA secondary structure of three copies of O. officinalis (AF363841, AF363902, and AF363882) are shown in Figure 4. The proposed secondary structure of one copy (AF363841) lost the B13_1 region (Figure 4a). The other changes also happened in B10, B11, B12, C1_1, and C1_3 regions. These changes may influence the interaction between rRNA and mRNA during the protein translation.
The consequences of changes in the region of the LSU rRNAs sequenced so far are unknown. However, a comparison of representative secondary structures with the variability map for eukaryotic LSU [34] gave no indica tion that LSU rRNA family genes are nonfunctional. Most conserved stem and loop structures in helixes B, C, and D are preserved with most of the differences in the variable parts. A further indication that LSU rRNA family genes in rice might be functional came from expression analyses.
Preprints 94623 g004
Figure 4. Proposed secondary structure of three copies of rRNA genes in O. officinalis. Arrows indicated the differences among the related loops or stems of secondary structure (B10, B11, B12, B13_1, C1_1 to C1_3). a, O. officinalis AF363841 represented cluster RC1; b, O. officinalis AF363902 represented cluster RC2; c, O. officinalis AF363882 represented cluster RC3.
Figure 4. Proposed secondary structure of three copies of rRNA genes in O. officinalis. Arrows indicated the differences among the related loops or stems of secondary structure (B10, B11, B12, B13_1, C1_1 to C1_3). a, O. officinalis AF363841 represented cluster RC1; b, O. officinalis AF363902 represented cluster RC2; c, O. officinalis AF363882 represented cluster RC3.
Preprints 94623 g004

3.6. Differential Expression of 28S ribosomal Genes with and without Fertilization

Applying a different fertilization regime with and without combined nitrogen on pot-grown O. officinalis plants revealed no differences in expression of RC1, RC2, and RC3 of 28S ribosomal genes (Figure 3). But further inspection of ribosomal gene expression after RC2-specific reverse transcription PCR followed by a terminated restriction fragments length polymorphism (T-RFLP) analyses of the products showed that dependent on N fertilization treatment, some RC2 28S rRNA gene were differentially expressed (Figure 5A,B). Differential expression was also confirmed by northern hybridization with specific 5’-digoxigenin labeled primers targeting sequences AF363902 and AF363911. 400 bp of readable sequence obtained from unfertilized plants are identical to AF363902 in RC2. Readable sequences from fertilized plants correspond to AF363902 and AF363911. This result suggests that switching of rRNA transcription units does not occur between rDNA families but within rDNA families in rice. In eukaryotes, switching of ribosome types has been first observed in Plasmodium [35], where several different rRNA genes are differentially transcribed during parasite development. However differential regulation of rRNA gene transcription within a repeat type has been not reported yet. It is unlikely that this mode of transcriptional regulation of rRNA genes happens by chance. Because nitrogen is the nutrient that is most limiting to plant growth in many ecosystems, a nitrogen-mediated regulation of rRNA expression may reflect an adaptation at the translational level to balance the supply of biologically available nitrogen under adverse conditions.

4. Discussion

The relationships of species within the genus Oryza which consists of 16 diploid and 7 allotetraploid species could not be resolved in 28S rRNA genes constructed phylogenetic tree. In the diploid species, O. longistaminata, O. officinalis, or O. australiensis from Africa, Asia, and Australia, respectively, three rDNA families RC1, RC2, and RC3 were detected within the same species. As in Quercus, the intra- and interspecific sequence differences within the clades were usually very low, which was explained for the two oak species by ongoing homogenization within but not between rDNA families [19]. The same observation is exhibited in rice. Some of the sequences on these branches namely from O. officinalis, O. sativa, and O. alta were identical to sequences from other species (O. longistaminata, O. longiglumis, and O. glaberrima respectively), which may be explained by species interbreeding, or very ancient speciation. In contrast to the Quercus species, which often hybridized with each other [36], natural crosses between wild rice species are unknown [37]. However, natural hybridization is common between wild rice and the worldwide cultivated O. sativa [38]. Species interbreeding as the reason for very similar or even identical 28S rRNA gene sequences in different species of Oryza can be excluded because comparative genomic mapping with AFLP markers had shown pre viously that all rice species have diverged independently from each other [39]. Recent speciation can be excluded because the phylogenetic analyses of the AFLP markers clearly indicated that the evolution of Oryza followed a polyphyletic pattern with eight monophyletic groups of species including those with identical sequences being separated from each other already early in evolution [40]. Furthermore, the finding of RC1 homologs in the draft genome sequences of O. sativa ssp. Indica [41] and japonica [42] confirmed the presence of RC1 in this species. Comparing 28S sequence AF363845 of RC1 rDNA family from O. sativa L. spp. indica IR36 to the draft genome sequence of the same subspecies by means of BLASTN gave a hit with 99.6% identity to the blasted 851 nucleotides in contig 21607 (http://btn.genomics.org.cn /blast/blast.php?name=rice) and a hit with 98% identity in contig HTC1232751301.1.1 when this sequence was blasted against the genome sequence of O. sativa L. ssp. japonica (http://www.tmri.org /index.html). No other high hits (>90% identity) were found, suggesting that RC1 genes are either not high copy genes or have not been sequenced or assembled yet. However, due to the high signal strength observed in Southern hybridization (Figure 2) at short exposure times (<5 min), the possibility that RC1-3 genes represent single copy genes can be excluded.
Sequence analyses of a partial PCR-amplified RDF2 tandem unit from O. officinalis (SI-2B) identified a unique truncated 18S rRNA gene in standard head to tail orientation upstream of the RDF2 18S rRNA gene. We found two distinct INDELs (insertion/deletion) in the P-site required for tRNA binding, one 48-base-pair (bp) deletion and one 14 bp insertion (SI-2A). Southern (SI-2C) and Northern (SI-2D) hybridization with a specific probe, revealed that this unique truncated 18S rDNA is widespread in rice and occurs even in maize. This result confirms that in the last common ancestor of maize, rice and all other cereals 50 Mya, this unique truncated 18S rRNA gene was already present, supporting the ancient origin of the rDNA families.
These results showed unambiguously, that in the diploid genomes of O. officinalis, O. longistaminata, O. glaberrima, O. sativa, and O. rufipogon all three rDNA families are present. Genomic southern hybridization also indicated that the rDNA families are more widespread in the genus Oryza than the initial PCR results suggested. The high diversity of 28S rRNA genes is consistent with a high variability in the positions of rDNA loci on the chromosomes in the genus Oryza [43]. This variability may be even higher than expected from fluorescence in situ hybridization, because for example in O. officinalis only three 28S rDNA clusters on chromosomes 4, 9 and 11 have been located [43] as opposed to the 4 highly diverged 28S rDNA sequences detected in our study.
Our results also showed that the rDNA families are ancient and predate the species split in the genus Oryza, because 28S gene sequences of several rice species were distributed among several rDNA families. The very little intraspecific sequence variability (Figure 1) within RC1-3 in combination with the high signal strength in southern hybridization (Figure 2) suggests high rates of concerted evolution within each rDNA family in spite of a large sequence divergence among families. In the rRNA multigene family of eukaryotes arrays of paralogous gene copies with low sequence divergence are typically homogenized by concerted evolution to nearly identical sequences within a species. Therefore, the presence of highly diverged copies within a multigene family would indicate that recombination between arrays is somehow suppressed, which may be attributed to a high sensitivity of the mismatch repair system to high sequence divergence [44] or to a non-telomere location on the chromosome [45]. This would suggest that diverged members of rRNA multigene families within a species are always xenologous, because a paralogous origin would be difficult to explain otherwise. Accordingly, the most plausible explanation for the evolutionary origin of the divergent rDNA families in rice as in Quercus [19] is that they were brought together by an ancient hybridization event.
In order to get a clue about the timing of this evolutionary event we tried to detect the xenologous rDNA families in another grass. The well-established phylogeny on the grass family and molecular-clock-based estimates indicate that maize and rice share their last common ancestor with the vast majority of the grass family and emerged 50 My ago, respectively [46]. These evolutionary relationships make maize an excellent candidate to evaluate the distribution of the rDNA families of rice within grasses. Rice species with several identical rDNA families in their genome, like O. officinalis, O. longistaminata, and O. australiensis, and possibly also O. grandiglumis are restricted to Asia, Africa, Australia, and South America ecological region, respectively [47]. Chang suggested that because of this allopatric distribution, speciation in the genus Oryza might have been driven by the successive breakup of Gondwanaland, which would place the origin of grasses into the Jurassic period (~140-210 My ago) (Table 1), much earlier than current estimates based on the first appearance of fossil grass pollen grains at the end of the Cretaceous (55-70 Mya) [48]. Monocots are probably not well preserved in the fossil record in the early Creteous (125-115 million years) [49], which raises the possibility that similar to vertebrate evolution [50], divergence times may significantly precede the appearance of relevant groups in the fossil record. The high statistical support for the branching nodes G within the rDNA families of maize-rice (Figure 1) allows us to infer maize and Oryza spices existed before the successive break-ups of Gondwanaland (SI-3). Therefore it was a possible Gondwanaland ancestry of rice based on our data.

5. Conclusions

Natural crosses between wild species of Oryza with different genomic types are unknown. However natural hybridization is common between wild species of rice and the worldwide cultivated O. sativa. Several evidence indicate that the distribution of the rRNA families in rice species was not a consequence of interspecific hybridization after speciation of the genus Oryza. (i) Comparative genomic fingerprints strongly support an independent evolution of Oryza species [51]. (ii) Before the onset of the worldwide distribution of O. sativa about 300 years ago, the geographic isolation of many rice species and short-lived pollen at least in Oryza sativa represented constraints making cross-con tinental interspecific hybridization within the genus Oryza by long-range transport of pollen. (iii) Chromosomal incompatibilities between rice species like O. sativa and O. officinalis or O. granulata either produce completely male sterile hybrids in sexual crosses with O. sativa or totally prevent species interbreeding. Nonetheless, the identity of partial LSU rDNA sequences affiliated with different rRNA families from different species is striking. If we consider the long span of genetic isolation after continent separation, the presence of identical rRNA gene sequences in different species appears to be unlikely but could be explained by the exchange of rRNA families through natural species interbreeding.

Supplementary Materials

The following supporting information can be downloaded at the website of this paper posted on Preprints.org. Supplementary 1 (Table S1: Used primers in this study); Supplementary 2; Figure SI-1. Codominant expression of RDF1-3 in rice. Presence of truncated RDF2 18S in ribosome preparations. (A) The portion of the aligned matrix from a unique truncated O. officinalis 18S ribosomal gene linked to RDF2 28S rRNA (AY097331) shows two INDELs in the decoding region that distinguishes it from other 18S rRNA genes. Dots indicate nucleotides identical to the top sequence, and dashes indicate gaps. The deleted and a putative alternative P-site are labeled. (B-D) Occurrence of the truncated 18S rRNA gene in rice. (B) Control Southern blot of EcoR I digested from rDNA clones with normal and truncated 18S rRNA genes (AY097331) hybridized with the RDF2 18S probe. Vector labeled with a star. (C). Presence of the truncated 18S rRNA in ribosome preparations from leaves of Zea mays and O. longistaminata (Oryza l.) grown in soil fertilized (+N) or not fertilized (-N) with combined nitrogen (D). Agarose gels (upper panel) and Northern blots (lower panel).; Figure S1: title

Author Contributions

“Conceptualization, Xiyu Tan and Guixiang Peng.; methodology, Muhammad Sajjid and Xiyu Tan.; software, Sidra Kaleem.; validation, Mehmood Jan., Raheel Munir and Arif Ali Khattak.; formal analysis, Abid Ali Abbas.; investigation, Muhammad Afzal and Xiyu Tan.; resources, Yihang Chen and Xiaolin Wang.; data curation, Muhammad Afzal.; writing—original draft preparation, Xiyu Tan.; writing—review and editing, Muhammad Afzal, Sidra Kaleem and Muhammad Sajjid.; visualization, Muhammad Afzal.; supervision, Muhammad Afzal and Zhiyuan Tan.; project administration, Muhammad Afzal and Zhiyuan Tan.; funding acquisition, Zhiyuan Tan. All authors have read and agreed to the published version of the manuscript.”

Funding

This work was supported by the National Natural Science Foundation of China (31970001) and National Key Research and Development Program (2022YFC3703100) and Science and Technology Project of Guangdong Province (2018A050506075).

Acknowledgments

Thanks to B. Reinhold-Hurek and T. Hurek for kindly instruction in the study.

Conflicts of Interest

“The authors declare no conflict of interest.”

References

  1. Wing, R.A.; Purugganan, M.D.; Zhang, Q. The rice genome revolution: from an ancient grain to Green Super Rice. Nature Reviews Genetics 2018, 19, 505–517. [Google Scholar] [CrossRef]
  2. Stein, J.C.; Yu, Y.; Copetti, D.; Zwickl, D.J.; Zhang, L.; Zhang, C.; Chougule, K.; Gao, D.; Iwata, A.; Goicoechea, J.L. Genomes of 13 domesticated and wild rice relatives highlight genetic conservation, turnover and innovation across the genus Oryza. Nature genetics 2018, 50, 285–296. [Google Scholar] [CrossRef]
  3. Tabassum, J.; Raza, Q.; Riaz, A.; Ahmad, S.; Rashid, M.A.R.; Javed, M.A.; Ali, Z.; Kang, F.; Khan, I.A.; Atif, R.M. Exploration of the genomic atlas of Dof transcription factor family across genus Oryza provides novel insights on rice breeding in changing climate. Frontiers in Plant Science 2022, 13, 1004359. [Google Scholar] [CrossRef]
  4. Gross, B.L.; Zhao, Z. Archaeological and genetic insights into the origins of domesticated rice. Proceedings of the National Academy of Sciences 2014, 111, 6190–6197. [Google Scholar] [CrossRef]
  5. Zhang, T.; Li, X.; Zhao, Z.; Wu, R.; Yang, Z.; He, G. Sequencing and Genomic Analysis of Sorghum DNA Introgression Variant Line R21 and Recipient Rice Jin Hui 1 Revealed Repetitive Element Variation. International Journal of Molecular Sciences 2022, 23, 11864. [Google Scholar] [CrossRef] [PubMed]
  6. Konikkat, S.; Woolford Jr, J.L. Principles of 60S ribosomal subunit assembly emerging from recent studies in yeast. Biochemical Journal 2017, 474, 195–214. [Google Scholar] [CrossRef] [PubMed]
  7. Mitterer, V.; Pertschy, B. RNA folding and functions of RNA helicases in ribosome biogenesis. RNA biology 2022, 19, 781–810. [Google Scholar] [CrossRef]
  8. Espinar-Marchena, F.J.; Fernández-Fernández, J.; Rodríguez-Galán, O.; Fernández-Pevida, A.; Babiano, R.; de la Cruz, J. Role of the yeast ribosomal protein L16 in ribosome biogenesis. FEBS J . 2016, 283, 2968–2985. [Google Scholar] [CrossRef] [PubMed]
  9. Wan, K.; Yabuki, Y.; Mizuta, K. Roles of Ebp2 and ribosomal protein L36 in ribosome biogenesis in Saccharomyces cerevisiae. Current genetics 2015, 61, 31–41. [Google Scholar] [CrossRef] [PubMed]
  10. Wang, X.; Yue, Z.; Xu, F.; Wang, S.; Hu, X.; Dai, J.; Zhao, G. Coevolution of ribosomal RNA expansion segment 7L and assembly factor Noc2p specializes the ribosome biogenesis pathway between Saccharomyces cerevisiae and Candida albicans. Nucleic Acids Research 2021, 49, 4655–4667. [Google Scholar] [CrossRef] [PubMed]
  11. Gómez Ramos, L.M.; Degtyareva, N.N.; Kovacs, N.A.; Holguin, S.Y.; Jiang, L.; Petrov, A.S.; Biesiada, M.; Hu, M.Y.; Purzycka, K.J.; Arya, D.P. Eukaryotic ribosomal expansion segments as antimicrobial targets. Biochemistry 2017, 56, 5288–5299. [Google Scholar] [CrossRef]
  12. Leppek, K.; Fujii, K.; Quade, N.; Susanto, T.T.; Boehringer, D.; Lenarčič, T.; Xue, S.; Genuth, N.R.; Ban, N.; Barna, M. Gene-and species-specific Hox mRNA translation by ribosome expansion segments. Molecular cell 2020, 80, 980–995. e913. [Google Scholar] [CrossRef]
  13. Paul, B.; Raj, K.K.; Murali, T.S.; Satyamoorthy, K. Species-specific genomic sequences for classification of bacteria. Computers in Biology and Medicine 2020, 123, 103874. [Google Scholar] [CrossRef]
  14. Raza, Q.; Rashid, M.A.R.; Waqas, M.; Ali, Z.; Rana, I.A.; Khan, S.H.; Khan, I.A.; Atif, R.M. Genomic diversity of aquaporins across genus Oryza provides a rich genetic resource for development of climate resilient rice cultivars. BMC Plant Biology 2023, 23, 172. [Google Scholar] [CrossRef]
  15. Chen, E.; Huang, X.; Tian, Z.; Wing, R.A.; Han, B. The genomics of Oryza species provides insights into rice domestication and heterosis. Annual review of plant biology 2019, 70, 639–665. [Google Scholar] [CrossRef] [PubMed]
  16. Wei, C.; Wang, Z.; Wang, J.; Teng, J.; Shen, S.; Xiao, Q.; Bao, S.; Feng, Y.; Zhang, Y.; Li, Y. Conversion between 100-million-year-old duplicated genes contributes to rice subspecies divergence. BMC genomics 2021, 22, 1–19. [Google Scholar] [CrossRef] [PubMed]
  17. He, L.; Hörandl, E. Does polyploidy inhibit sex chromosome evolution in angiosperms? Frontiers in Plant Science 2022, 13, 976765. [Google Scholar] [CrossRef] [PubMed]
  18. Bowers, J.E.; Tang, H.; Burke, J.M.; Paterson, A.H. GC content of plant genes is linked to past gene duplications. Plos one 2022, 17, e0261748. [Google Scholar] [CrossRef] [PubMed]
  19. Zhou, J.; Yang, Y.; Lv, Y.; Pu, Q.; Li, J.; Zhang, Y.; Deng, X.; Wang, M.; Wang, J.; Tao, D. Interspecific Hybridization Is an Important Driving Force for Origin and Diversification of Asian Cultivated Rice Oryza sativa L. Frontiers in Plant Science 2022, 13, 932737. [Google Scholar] [CrossRef] [PubMed]
  20. Palaniyappan, S.; Arunachalam, P.; Banumathy, S.; Mini, M.; Muthuramu, S. Genetic divergence and clustering studies in advanced breeding lines of rice (Oryza sativa L.). Electronic Journal of Plant Breeding 2020, 11, 499–504. [Google Scholar] [CrossRef]
  21. Cheng, L.; Han, Q.; Chen, F.; Li, M.; Balbuena, T.S.; Zhao, Y. Phylogenomics as an effective approach to untangle cross-species hybridization event: A case study in the family Nymphaeaceae. Frontiers in Genetics 2022, 13, 1031705. [Google Scholar] [CrossRef]
  22. Sun, K.; Li, D.; Xia, A.; Zhao, H.; Wen, Q.; Jia, S.; Wang, J.; Yang, G.; Zhou, D.; Huang, C. Targeted Identification of Rice Grain-Associated Gene Allelic Variation Through Mutation Induction, Targeted Sequencing, and Whole Genome Sequencing Combined with a Mixed-Samples Strategy. Rice 2022, 15, 1–15. [Google Scholar] [CrossRef] [PubMed]
  23. Zeng, H.; Chen, X.; Li, H.; Zhang, J.; Wei, Z.; Wang, Y. Interpopulation differences of retroduplication variations (RDVs) in rice retrogenes and their phenotypic correlations. Computational and Structural Biotechnology Journal 2021, 19, 600–611. [Google Scholar] [CrossRef] [PubMed]
  24. Parks, M.M.; Kurylo, C.M.; Batchelder, J.E.; Theresa Vincent, C.; Blanchard, S.C. Implications of sequence variation on the evolution of rRNA. Chromosome Research 2019, 27, 89–93. [Google Scholar] [CrossRef] [PubMed]
  25. Shi, C.; Li, W.; Zhang, Q.-J.; Zhang, Y.; Tong, Y.; Li, K.; Liu, Y.-L.; Gao, L.-Z. The draft genome sequence of an upland wild rice species, Oryza granulata. Scientific Data 2020, 7, 131. [Google Scholar] [CrossRef]
  26. Jagadeesh, D.; Kumar, M.P.; Chandrakanth, R.; Devaki, N. Molecular diversity of internal transcribed spacer among the monoconidial isolates of Magnaporthe oryzae isolated from rice in Southern Karnataka, India. Journal of Genetic Engineering and Biotechnology 2018, 16, 631–638. [Google Scholar] [CrossRef]
  27. Jayaswal, P.K.; Dogra, V.; Shanker, A.; Sharma, T.R.; Singh, N.K. A tree of life based on ninety-eight expressed genes conserved across diverse eukaryotic species. Plos one 2017, 12, e0184276. [Google Scholar] [CrossRef] [PubMed]
  28. Zou, X.-H.; Du, Y.-S.; Tang, L.; Xu, X.-W.; Doyle, J.J.; Sang, T.; Ge, S. Multiple origins of BBCC allopolyploid species in the rice genus (Oryza). Scientific reports 2015, 5, 14876. [Google Scholar] [CrossRef]
  29. Jain, R.; Jenkins, J.; Shu, S.; Chern, M.; Martin, J.A.; Copetti, D.; Duong, P.Q.; Pham, N.T.; Kudrna, D.A.; Talag, J. Genome sequence of the model rice variety KitaakeX. BMC genomics 2019, 20, 1–9. [Google Scholar] [CrossRef]
  30. Luo, S.; Peng, J.; Li, K.; Wang, M.; Kuang, H. Contrasting evolutionary patterns of the Rp1 resistance gene family in different species of Poaceae. Molecular biology and evolution 2011, 28, 313–325. [Google Scholar] [CrossRef]
  31. Piya, A.A.; DeGiorgio, M.; Assis, R. Predicting gene expression divergence between single-copy orthologs in two species. Genome Biology and Evolution 2023, 15, evad078. [Google Scholar] [CrossRef]
  32. Zhao, D.; Hamilton, J.P.; Hardigan, M.; Yin, D.; He, T.; Vaillancourt, B.; Reynoso, M.; Pauluzzi, G.; Funkhouser, S.; Cui, Y. Analysis of ribosome-associated mRNAs in rice reveals the importance of transcript size and GC content in translation. G3: Genes, Genomes, Genetics 2017, 7, 203–219. [Google Scholar] [CrossRef]
  33. Siwaszek, A.; Ukleja, M.; Dziembowski, A. Proteins involved in the degradation of cytoplasmic mRNA in the major eukaryotic model systems. RNA biology 2014, 11, 1122–1136. [Google Scholar] [CrossRef]
  34. Deng, H.; Cheema, J.; Zhang, H.; Woolfenden, H.; Norris, M.; Liu, Z.; Liu, Q.; Yang, X.; Yang, M.; Deng, X. Rice in vivo RNA structurome reveals RNA secondary structure conservation and divergence in plants. Molecular plant 2018, 11, 607–622. [Google Scholar] [CrossRef]
  35. McGee, J.P.; Armache, J.-P.; Lindner, S.E. Ribosome heterogeneity and specialization of Plasmodium parasites. PLoS pathogens 2023, 19, e1011267. [Google Scholar] [CrossRef] [PubMed]
  36. Burgarella, C.; Lorenzo, Z.; Jabbour-Zahab, R.; Lumaret, R.; Guichoux, E.; Petit, R.; Soto, A.; Gil, L. Detection of hybrids in nature: application to oaks (Quercus suber and Q. ilex). Heredity 2009, 102, 442–452. [Google Scholar] [CrossRef]
  37. He, N.; Wu, R.; Pan, X.; Peng, L.; Sun, K.; Zou, T.; Zhu, H.; Zeng, R.; Liu, Z.; Liu, G. Development and trait evaluation of chromosome single-segment substitution lines of O. meridionalis in the background of O. sativa. Euphytica 2017, 213, 1–15. [Google Scholar] [CrossRef]
  38. Kopecký, D.; Martín, A.; Smýkal, P. Interspecific hybridization and plant breeding: From historical retrospective through work of Mendel to current crops. Czech Journal of Genetics and Plant Breeding 2022. [Google Scholar] [CrossRef]
  39. Vemireddy, L.R.; Noor, S.; Satyavathi, V.; Srividhya, A.; Kaliappan, A.; Parimala, S.; Bharathi, P.M.; Deborah, D.A.; Rao, K.S.; Shobharani, N. Discovery and mapping of genomic regions governing economically important traits of Basmati rice. BMC plant biology 2015, 15, 1–19. [Google Scholar] [CrossRef] [PubMed]
  40. Siddiq, E.; Vemireddy, L.; Nagaraju, J. Basmati rices: genetics, breeding and trade. Agricultural Research 2012, 1, 25–36. [Google Scholar] [CrossRef]
  41. Reddy, M.M.; Ulaganathan, K. Draft genome sequence of Oryza sativa elite indica cultivar RP Bio-226. Frontiers in plant science 2015, 6, 896. [Google Scholar] [CrossRef] [PubMed]
  42. Wang, Y. A draft genome, resequencing, and metabolomes reveal the genetic background and molecular basis of the nutritional and medicinal properties of loquat (Eriobotrya japonica (Thunb.) Lindl). Horticulture research 2021, 8. [Google Scholar] [CrossRef] [PubMed]
  43. Dai, S.-f.; Zhu, X.-g.; Hutang, G.-r.; Li, J.-y.; Tian, J.-q.; Jiang, X.-h.; Zhang, D.; Gao, L.-z. Genome size variation and evolution driven by transposable elements in the genus oryza. Frontiers in Plant Science 2022, 13, 921937. [Google Scholar] [CrossRef] [PubMed]
  44. Guo, Y.L. Gene family evolution in green plants with emphasis on the origination and evolution of a rabidopsis thaliana genes. The Plant Journal 2013, 73, 941–951. [Google Scholar] [CrossRef] [PubMed]
  45. Ohmido, N.; Fukui, K.; Kinoshita, T. Recent advances in rice genome and chromosome structure research by fluorescence in situ hybridization (FISH). Proceedings of the Japan Academy, Series B 2010, 86, 103–116. [Google Scholar] [CrossRef]
  46. Huang, W.; Zhang, L.; Columbus, J.T.; Hu, Y.; Zhao, Y.; Tang, L.; Guo, Z.; Chen, W.; McKain, M.; Bartlett, M. A well-supported nuclear phylogeny of Poaceae and implications for the evolution of C4 photosynthesis. Molecular Plant 2022, 15, 755–777. [Google Scholar] [CrossRef]
  47. Gaikwad, K.B.; Singh, N.; Kaur, P.; Rani, S.; Babu H, P.; Singh, K. Deployment of wild relatives for genetic improvement in rice (Oryza sativa). Plant Breeding 2021, 140, 23–52. [Google Scholar] [CrossRef]
  48. Yang, X.; Fuller, D.Q.; Huan, X.; Perry, L.; Li, Q.; Li, Z.; Zhang, J.; Ma, Z.; Zhuang, Y.; Jiang, L. Barnyard grasses were processed with rice around 10000 years ago. Scientific Reports 2015, 5, 16251. [Google Scholar] [CrossRef]
  49. Coiffard, C.; Kardjilov, N.; Manke, I.; Bernardes-de-Oliveira, M.E. Fossil evidence of core monocots in the Early Cretaceous. Nature plants 2019, 5, 691–696. [Google Scholar] [CrossRef] [PubMed]
  50. Hertweck, K.L.; Kinney, M.S.; Stuart, S.A.; Maurin, O.; Mathews, S.; Chase, M.W.; Gandolfo, M.A.; Pires, J.C. Phylogenetics, divergence times and diversification from three genomic partitions in monocots. Botanical Journal of the Linnean Society 2015, 178, 375–393. [Google Scholar] [CrossRef]
  51. Gao, L.-Z.; Liu, Y.-L.; Zhang, D.; Li, W.; Gao, J.; Liu, Y.; Li, K.; Shi, C.; Zhao, Y.; Zhao, Y.-J. Evolution of Oryza chloroplast genomes promoted adaptation to diverse ecological habitats. Communications Biology 2019, 2, 278. [Google Scholar] [CrossRef] [PubMed]
Preprints 94623 g002
Figure 2. Southern hybridization showed the occurrence of the rDNA families RC1, RC2 and RC3 in rice genome with specific probes RCp1, RCp2 and RCp3, respectively; 1, 2, 3, 4, 5, 6 indicated genomic DNA from O. officinalis, O. sativa subsp. indica IR36, O. longistaminata, O. rufipogon from Nepal, O. rufipogon from IRRI and O. glaberrima digested with Bse RI and Xcm I, separately; agarose gel showed the electrophoresis of digested genomic DNA; M, lambda DNA/Pst I Marker, 29 fragments (size in bp): 11501, 5077, 4749, 4507, 2838, 2556, 2459, 2443, 2140, 1986, 1700, 1159, 1093, 805, 514, 468, 448, 339, 264, 247, 216, 211, 200, 164, 150, 94, 87, 72, 15.
Figure 2. Southern hybridization showed the occurrence of the rDNA families RC1, RC2 and RC3 in rice genome with specific probes RCp1, RCp2 and RCp3, respectively; 1, 2, 3, 4, 5, 6 indicated genomic DNA from O. officinalis, O. sativa subsp. indica IR36, O. longistaminata, O. rufipogon from Nepal, O. rufipogon from IRRI and O. glaberrima digested with Bse RI and Xcm I, separately; agarose gel showed the electrophoresis of digested genomic DNA; M, lambda DNA/Pst I Marker, 29 fragments (size in bp): 11501, 5077, 4749, 4507, 2838, 2556, 2459, 2443, 2140, 1986, 1700, 1159, 1093, 805, 514, 468, 448, 339, 264, 247, 216, 211, 200, 164, 150, 94, 87, 72, 15.
Preprints 94623 g002
Preprints 94623 g003
Figure 3. Northern hybridization showed the co-dominant expression of rDNA families RC1, RC2 and RC3 in rice plants. A, B, C, total RNA extracted from O. officinalis with specific probe RCp1, RCp2 and RCp3 hybridization, respectively; +N, -N indicated O. officinalis grown in soil fertilized and not fertilized. Agarose gels (upper panel) and Northern blots (lower panel).
Figure 3. Northern hybridization showed the co-dominant expression of rDNA families RC1, RC2 and RC3 in rice plants. A, B, C, total RNA extracted from O. officinalis with specific probe RCp1, RCp2 and RCp3 hybridization, respectively; +N, -N indicated O. officinalis grown in soil fertilized and not fertilized. Agarose gels (upper panel) and Northern blots (lower panel).
Preprints 94623 g003
Preprints 94623 g005
Figure 5. 28S rDNA-based cDNA sequence printouts (A) and corresponding terminal restriction fragment length polymorphism (T-RFLP) fingerprint patterns (B) showed the differential expression of 28S rRNA genes in rDNA family RC2. 28S rRNA genes extracted from roots and leaves of O. officinalis grown in soil fertilized (+N) and unfertilized (-N) with combined nitrogen. For T-RFLP analyses restriction sites of Pvu II, Bcl I and Bst XI were underlined in the sequence printout.
Figure 5. 28S rDNA-based cDNA sequence printouts (A) and corresponding terminal restriction fragment length polymorphism (T-RFLP) fingerprint patterns (B) showed the differential expression of 28S rRNA genes in rDNA family RC2. 28S rRNA genes extracted from roots and leaves of O. officinalis grown in soil fertilized (+N) and unfertilized (-N) with combined nitrogen. For T-RFLP analyses restriction sites of Pvu II, Bcl I and Bst XI were underlined in the sequence printout.
Preprints 94623 g005
Table 1. Dates of divergence of Oryza 28S rDNA lineages RC1, RC2 and RC3.
Table 1. Dates of divergence of Oryza 28S rDNA lineages RC1, RC2 and RC3.
Preprints 94623 i001
a Numbers in brackets give the number of sequences used for the analysis as represented in Figure 1. bA divergence time of 170-235 My ago was assumed for the monocot-dicot split (Yang et al., 1999).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

facebook logotwitter logolinkedin logo
MDPI Initiatives
Important Links
Subscribe

Disclaimer

Terms of Use

Privacy Policy

© 2024 MDPI (Basel, Switzerland) unless otherwise stated