Early evolutionary colocalization of the nuclear ribosomal 5S and 45S gene families in seed plants: evidence from the living fossil gymnosperm Ginkgo biloba

Intragenomic rDNA variation - the product of concerted evolution, mutation, or something in between?

The classical model of concerted evolution states that hundreds to thousands of ribosomal DNA (rDNA) units undergo homogenization, making the multiple copies of the individual units more uniform across the genome than would be expected given mutation frequencies and gene redundancy. While the universality of this over 50-year-old model has been confirmed in a range of organisms, advanced high throughput sequencing techniques have also revealed that rDNA homogenization in many organisms is partial and, in rare cases, even apparently failing. The potential underpinning processes leading to unexpected intragenomic variation have been discussed in a number of studies, but a comprehensive understanding remains to be determined. In this work, we summarize information on variation or polymorphisms in rDNAs across a wide range of taxa amongst animals, fungi, plants, and protists. We discuss the definition and description of concerted evolution and describe whether incomplete concerted evolution of rDNAs predominantly affects coding or non-coding regions of rDNA units and if it leads to the formation of pseudogenes or not. We also discuss the factors contributing to rDNA variation, such as interspecific hybridization, meiotic cycles, rDNA expression status, genome size, and the activity of effector genes involved in genetic recombination, epigenetic modifications, and DNA editing. Finally, we argue that a combination of approaches is needed to target genetic and epigenetic phenomena influencing incomplete concerted evolution, to give a comprehensive understanding of the evolution and functional consequences of intragenomic variation in rDNA.

Developing population identification tool based on polymorphism of rDNA for traditional Chinese medicine: Artemisia annua L

BACKGROUND

Artemisia annua, a well-known traditional Chinese medicine, is the main source for production of artemisinin, an anti-malaria drug. A. annua is distributed globally, with great diversity of morphological characteristics and artemisinin contents. Diverse traits among A. annua populations impeded the stable production of artemisinin, which needs an efficient tool to identify strains and assess population genetic homogeneity.

PURPOSE

In this study, ribosomal DNA (rDNA), were characterized for A. annua for strains identification and population genetic homogeneity assessment.

METHODS

The ribosomal RNA (rRNA) genes were identified using cmscan and assembled using rDNA unit of LQ-9 as a reference. rDNA among Asteraceae species were compared performing with 45S rDNA. The rDNA copy number was calculated based on sequencing depth. The polymorphisms of rDNA sequences were identified with bam-readcount, and confirmed by Sanger sequencing and restriction enzyme experiment. The ITS2 amplicon sequencing was used to verify the stability of ITS2 haplotype analysis.

RESULTS

Different from other Asteraceae species, 45S and 5S linked-type rDNA was only found in Artemisia genus. Rich polymorphisms of copy number and sequence of rDNA were identified in A. annua population. The haplotype composition of internal transcribed spacer 2 (ITS2) region which had moderate sequence polymorphism and relative short size was significantly different among A. annua strains. A population discrimination method was developed based on ITS2 haplotype analysis with high-throughput sequencing.

CONCLUSION

This study provides comprehensive characteristics of rDNA and suggests that ITS2 haplotype analysis is ideal tool for A. annua strain identification and population genetic homogeneity assessment.

Analysis of 5S rDNA Genomic Organization Through the RepeatExplorer2 Pipeline: A Simplified Protocol

The ribosomal RNA genes (rDNA) are universal genome components with a housekeeping function, given the crucial role of ribosomal RNA in the synthesis of ribosomes and thus for life-on-Earth. Therefore, their genomic organization is of considerable interest for biologists, in general. Ribosomal RNA genes have also been largely used to establish phylogenetic relationships, and to identify allopolyploid or homoploid hybridization.Here, we demonstrate how high-throughput sequencing data, through graph clustering implemented in RepeatExplorer2 pipeline ( https://repeatexplorer-elixir.cerit-sc.cz/galaxy/ ), can be helpful to decipher the genomic organization of 5S rRNA genes. We show that the linear shapes of cluster graphs are reminiscent to the linked organization of 5S and 35S rDNA (L-type arrangement) while the circular graphs correspond to their separate arrangement (S-type). We further present a simplified protocol based on the paper by (Garcia et al., Front Plant Sci 11:41, 2020) about the use of graph clustering of 5S rDNA homoeologs (S-type) to identify hybridization events in the species history. We found that the graph complexity (i.e., graph circularity in this case) is related to ploidy and genome complexity, with diploids typically showing circular-shaped graphs while allopolyploids and other interspecific hybrids display more complex graphs, with usually two or more interconnected loops representing intergenic spacers. When a three-genomic comparative clustering analysis from a given hybrid (homoploid/allopolyploid) and its putative progenitor species (diploids) is performed, it is possible to identify the corresponding homoeologous 5S rRNA gene families, and to elucidate the contribution of each putative parental genome to the 5S rDNA pool of the hybrid. Thus, the analysis of 5S rDNA cluster graphs by RepeatExplorer, together with information coming from other sources (e.g., morphology, cytogenetics) is a complementary approach for the determination of allopolyploid or homoploid hybridization and even ancient introgression events.

The Use of Ribosomal DNA for Comparative Cytogenetics

Fluorescence in situ hybridization (FISH) with ribosomal DNA (rDNA) sequences provides excellent chromosome markers for comparative cytogenetic analyses, especially in non-model plant species. The tandem repeat nature of a sequence and the presence of a highly conserved genic region make rDNA sequences relatively easy to isolate and clone. In this chapter, we describe the use of rDNA as markers for comparative cytogenetics studies. Traditionally, cloned probes labeled with Nick-translation have been used to detect rDNA loci. Recently, pre-labeled oligonucleotides are also employed quite frequently to detect both 35S and 5S rDNA loci. Ribosomal DNA sequences, together with other DNA probes in FISH/GISH or with fluorochromes such as CMA₃ banding or silver staining, are very useful tools in comparative analyses of plant karyotypes.

Revisiting the amalgaviral landscapes in plant transcriptomes expands the host range of plant amalgaviruses

Plant amalgaviruses are monopartite, double-stranded RNA viruses, capable of vertical transmission through seeds. An attempt to revisit plant transcriptome-assembled contigs for amalgaviral sequences identified 40 putative novel amalgaviruses in 35 plant species, nearly doubling the number of plant amalgaviruses. Of the 35 plant species, 33 are reported to host amalgaviruses for the first time, including a pteridophytic and two gymnospermic species. Coding-complete genomes of all identified viruses were recovered and the putative +1 programmed ribosomal frameshift (PRF) sites were determined. Genomes of 35 identified amalgaviruses contained the conserved +1 PRF motif 'UUU_CGN', while variant versions were predicted in five genomes. Phylogenetic analysis grouped pteridophyte- and gymnosperm-infecting amalgaviruses together in divergent sub-clades while few of the related angiosperm-infecting amalgaviruses infect members of the same plant family, reiterating the co-evolution of plant amalgaviruses and their hosts. The current study paves way for further studies on understanding biological properties of identified viruses.

The Divergence of Chromosome Structures and 45S Ribosomal DNA Organization in Cucumis debilis Inferred by Comparative Molecular Cytogenetic Mapping

Cucumis debilis W.J.de Wilde & Duyfjes is an annual and monoecious plant. This species is endemic to Southeast Asia, particularly Vietnam. However, C. debilis is rarely studied, and no detailed information is available regarding its basic chromosome number, 45S ribosomal DNA (rDNA) status, and divergence among other Cucumis species. In this study, we characterized the morphological characters and determined and investigated the basic chromosome number and chromosomal distribution of 45S rDNA of C. debilis using the fluorescent in situ hybridization (FISH) technique. A maximum likelihood tree was constructed by combining the chloroplast and internal transcribed spacer of 45S rDNAs to infer its relationship within Cucumis. C. debilis had an oval fruit shape, green fruit peel, and protrusion-like white spots during the immature fruit stage. FISH analysis using 45S rDNA probe showed three pairs of 45S rDNA loci located at the terminal region in C. debilis, similar to C. hystrix. Meanwhile, two, two, and five pairs of 45S rDNA loci were observed for C. melo, C. metuliferus, and C. sativus, respectively. One melon (P90) and cucumber accessions exhibited different chromosomal localizations compared with other members of Cucumis. The majority of Cucumis species showed the terminal location of 45S rDNA, but melon P90 and cucumber exhibited terminal–interstitial and all interstitial orientations of 45S rDNA loci. Based on molecular cytogenetics and phylogenetic evidence, C. debilis is more closely related to cucumber than melon. Therefore, C. debilis may serve as a potential parental accession for genetic improvement of cucumber through interspecific hybridization.

Genome structure and diversity among Cynanchum wilfordii accessions

BACKGROUND

Cynanchum wilfordii (Cw) and Cynanchum auriculatum (Ca) have long been used in traditional medicine and as functional food in Korea and China, respectively. They have diverse medicinal functions, and many studies have been conducted, including pharmaceutical efficiency and metabolites. Especially, Cw is regarded as the most famous medicinal herb in Korea due to its menopausal symptoms relieving effect. Despite the high demand for Cw in the market, both species are cultivated using wild resources with rare genomic information.

RESULTS

We collected 160 Cw germplasm from local areas of Korea and analyzed their morphological diversity. Five Cw and one Ca of them, which were morphologically diverse, were sequenced, and nuclear ribosomal DNA (nrDNA) and complete plastid genome (plastome) sequences were assembled and annotated. We investigated the genomic characteristics of Cw as well as the genetic diversity of plastomes and nrDNA of Cw and Ca. The Cw haploid nuclear genome was approximately 178 Mbp. Karyotyping revealed the juxtaposition of 45S and 5S nrDNA on one of 11 chromosomes. Plastome sequences revealed 1226 interspecies polymorphisms and 11 Cw intraspecies polymorphisms. The 160 Cw accessions were grouped into 21 haplotypes based on seven plastome markers and into 108 haplotypes based on seven nuclear markers. Nuclear genotypes did not coincide with plastome haplotypes that reflect the frequent natural outcrossing events.

CONCLUSIONS

Cw germplasm had a huge morphological diversity, and their wide range of genetic diversity was revealed through the investigation with 14 molecular markers. The morphological and genomic diversity, chromosome structure, and genome size provide fundamental genomic information for breeding of undomesticated Cw plants.

The Nuclear 35S rDNA World in Plant Systematics and Evolution: A Primer of Cautions and Common Misconceptions in Cytogenetic Studies

The ubiquitous presence of rRNA genes in nuclear, plastid, and mitochondrial genomes has provided an opportunity to use genomic markers to infer patterns of molecular and organismic evolution as well as to assess systematic issues throughout the tree of life. The number, size, location, and activity of the 35S rDNA cistrons in plant karyotypes have been used as conventional cytogenetic landmarks. Their scrutiny has been useful to infer patterns of chromosomal evolution and the data have been used as a proxy for assessing species discrimination, population differentiation and evolutionary relationships. The correct interpretation of rDNA markers in plant taxonomy and evolution is not free of drawbacks given the complexities derived from the lability of the genetic architecture, the diverse patterns of molecular change, and the fate and evolutionary dynamics of the rDNA units in hybrids and polyploid species. In addition, the terminology used by independent authors is somewhat vague, which often complicates comparisons. To date, no efforts have been reported addressing the potential problems and limitations involved in generating, utilizing, and interpreting the data from the 35S rDNA in cytogenetics. This review discusses the main technical and conceptual limitations of these rDNA markers obtained by cytological and karyological experimental work, in order to clarify biological and evolutionary inferences postulated in a systematic and phylogenetic context. Also, we provide clarification for some ambiguity and misconceptions in terminology usually found in published work that may help to improve the usage of the 35S ribosomal world in plant evolution.

Personal Perspectives on Plant Ribosomal RNA Genes Research: From Precursor-rRNA to Molecular Evolution

The history of rDNA research started almost 90 years ago when the geneticist, Barbara McClintock observed that in interphase nuclei of maize the nucleolus was formed in association with a specific region normally located near the end of a chromosome, which she called the nucleolar organizer region (NOR). Cytologists in the twentieth century recognized the nucleolus as a common structure in all eukaryotic cells, using both light and electron microscopy and biochemical and genetic studies identified ribosomes as the subcellular sites of protein synthesis. In the mid- to late 1960s, the synthesis of nuclear-encoded rRNA was the only system in multicellular organisms where transcripts of known function could be isolated, and their synthesis and processing could be studied. Cytogenetic observations of NOR regions with altered structure in plant interspecific hybrids and detailed knowledge of structure and function of rDNA were prerequisites for studies of nucleolar dominance, epistatic interactions of rDNA loci, and epigenetic silencing. In this article, we focus on the early rDNA research in plants, performed mainly at the dawn of molecular biology in the 60 to 80-ties of the last century which presented a prequel to the modern genomic era. We discuss - from a personal view - the topics such as synthesis of rRNA precursor (35S pre-rRNA in plants), processing, and the organization of 35S and 5S rDNA. Cloning and sequencing led to the observation that the transcribed and processed regions of the rRNA genes vary enormously, even between populations and species, in comparison with the more conserved regions coding for the mature rRNAs. Epigenetic phenomena and the impact of hybridization and allopolyploidy on rDNA expression and homogenization are discussed. This historical view of scientific progress and achievements sets the scene for the other articles highlighting the immense progress in rDNA research published in this special issue of Frontiers in Plant Science on "Molecular organization, evolution, and function of ribosomal DNA."

Gauging the trends of pseudogenes in plants

Pseudogenes, the debilitated parts of ancient genes, were previously scrapped off as junk or discarded genes with no functional significance. Pseudogenes have come under scrutiny for their functionality, since recent studies have unveiled their importance in the regulation of their corresponding parent genes and various biological mechanisms. Despite the enormous occurrence of pseudogenes in plants, the lack of experimental validation has contributed toward their unresolved roles in gene regulation. Contrarily, most of the studies associated with gene regulation have been mainly reported for humans, mice, and other mammalian genomes. Consequently, in order to present a cumulative report on plant-based pseudogenes research, an attempt has been made to assemble multiple studies presenting the pseudogene classification, the prediction and the determination of comparative accuracies of various computational pipelines, and recent trends in analyzing their biological functions, and regulatory mechanisms. This review represents the classical, as well as the recent advances on pseudogene identification and their potential roles in transcriptional regulation, which could possibly invigorate the quality of genome annotation, evolutionary analysis, and complexity surrounding the regulatory pathways in plants. Thus, when the ambiguous boundary girdling the pseudogenes eventually recedes on account of their explicit orchestration role, research in flora would no longer saunter compared to that on fauna.

Comparative Chromosomal Localization of 45S and 5S rDNA Sites in 76 Purple-Fleshed Sweet Potato Cultivars

In recent years, the purple-fleshed sweet potato has attracted more attention because of its high nutritional value. The cytogenetics of this crop is relatively unexplored, limiting our knowledge on its genetic diversity. Therefore, we conducted cytogenetic analysis of 76 purple-fleshed sweet potato cultivars to analyze the chromosome structure and distribution of 45S and 5S rDNA. We noted that only 62 cultivars had 90 chromosomes, and the others were aneuploid with 88, 89, 91, or 92 chromosomes. The number of 45S rDNA in the 76 cultivars varied from 16 to 21; these sites showed different signal sizes and intensities and were localized at the chromosomal termini or satellite. The number of 5S rDNA was relatively stable; 74 cultivars showed six sites located at the chromosomal sub-terminal or near the centromere. Only the ‘Quanzishu 96’ and ‘Yuzixiang 10’ showed seven and five 5S rDNA sites, respectively. Additionally, both parent cultivars of ‘Quanzishu 96’ showed 18 45S and six 5S rDNA sites. Overall, our results indicate a moderate diversity in the distribution pattern of rDNAs. Our findings provide comprehensive cytogenetic information for the identification of sweet potato chromosomes, which can be useful for developing a high-quality germplasm resource.

Advances in the cytogenetics of Annonaceae, the case of Annona cherimola L

Annonaceae represent the largest extant family among the early divergent angiosperms. Despite the long-standing interest in its evolutionary and taxonomic aspects, cytogenetic studies on this family remain extremely few even on economically important species. With this study, we realized a detailed characterization of the chromosomes of Annona cherimola (2n = 14) by a combination of in situ hybridization techniques, fluorochrome banding, and karyomorphological analysis. FISH revealed that 45S and 5S rDNA sites are co-localized in correspondence to the secondary constrictions of the SAT-chromosome pair. Some hypotheses on the organization of the linked 45S and 5S rDNA repeats have been made. FISH with Arabidopsis-type telomeric arrays demonstrated that the A. cherimola telomeres are constituted by TTTAGGG sequences and that they are exclusively localized at the extremities of chromosomes. An insight into the chromosome structure of A. cherimola was obtained by the self-GISH procedure which revealed highly repeated DNA sequences localized in the centromeric regions of all chromosomes. The correspondence of s-GISH signals with DAPI banding suggests that these sequences are the principal component of the centromeric heterochromatin of this species. The karyotype of A. cherimola here described is proposed as the basic reference karyotype for successive investigations in Annonaceae.

Cytogenetic analysis of Bienertia sinuspersici Akhani as the first step in genome sequencing

BACKGROUND

C₄ plants are efficient in suppressing photorespiration and enhancing carbon gain as compared to C₃ plants. Bienertia sinuspersici Akhani is one of the few species in the family Amaranthaceae that can perform C₄ photosynthesis within individual chlorenchyma cells, without the conventional Kranz anatomy in its leaf. This plant is salt-tolerant and is well-adapted to thrive in hot and humid climates. To date, there have been no reported cytogenetic analyses yet on this species.

OBJECTIVE

This study aims to provide a cytogenetic analysis of B. sinuspersici as the first step in genome sequencing.

METHODS

Fluorescence in situ hybridization (FISH) karyotype analysis was conducted using the metaphase chromosomes of B. sinuspersici probed with 5S and 45S rDNA and Arabidopsis-type telomeric repeats.

RESULTS

Results of the cytogenetic analysis confirmed that B. sinuspersici has 2n = 2x = 18 consisting of nine pairs of metacentric chromosomes. Two loci of 45S rDNA were found on the distal regions of the short arm of chromosome 7. Nine loci of 5S rDNA were found in the pericentromeric regions of chromosomes 1, 3, 4, 6, and 8, which also colocalized with Arabidopsis-type telomeric repeats; while four loci in the interstitial regions of chromosome 5 and 8 can be observed. The single locus of 5S rDNA that was found in chromosome 8 appears to be hemizygous.

CONCLUSION

The FISH karyotype analysis, based on the combination of rDNAs, telomeric tandem repeat markers and C₀t DNA chromosome landmarks, allowed efficient chromosome identification and provided useful information in characterizing the genome of B. sinuspersici.

In silico mining and FISH mapping of a chromosome-specific satellite DNA in Capsicum annuum L

BACKGROUND

A large proportion of eukaryote nuclear genomes is composed of repetitive DNA. Tracing the dynamics of repetitive elements in the genomes of related taxa can reveal important information about their phylogenic relationships as well as traits that have become distinct to a lineage.

OBJECTIVE

Study the genomic abundance and chromosomal location of repetitive DNA in Capsicum annuum L. to understand the repeat dynamics.

METHOD

We quantified repeated DNA content in the C. annuum genome using the RepeatExplorer pipeline.

RESULTS

About 42% of the C. annuum genome dataset comprised repetitive elements. Of these, 0.011, 0.98, 3.09, and 0.024% represented high and low confidence satellite repeats, putative long-terminal repeats (LTRs), and rDNA sequences, respectively. One novel high confidence 167-bp satellite repeat with a genomic proportion of 0.011%, Ca167TR, was identified. Furthermore, FISH with Ca167TR on metaphase chromosomes of C. annuum revealed signals in the subtelomeric regions of the short and long arms of chromosome 3 and 4, respectively.

CONCLUSION

Further understanding of the origin and associated functions of Ca167TR and other repeats in Capsicum will give us insights into the genomic relationships and functions of the genome.

Remarkable variation of ribosomal DNA organization and copy number in gnetophytes, a distinct lineage of gymnosperms

INTRODUCTION

Gnetophytes, comprising the genera Ephedra, Gnetum and Welwitschia, are an understudied, enigmatic lineage of gymnosperms with a controversial phylogenetic relationship to other seed plants. Here we examined the organization of ribosomal DNA (rDNA) across representative species.

METHODS

We applied high-throughput sequencing approaches to isolate and reconstruct rDNA units and to determine their intragenomic homogeneity. In addition, fluorescent in situ hybridization and Southern blot hybridization techniques were used to reveal the chromosome and genomic organization of rDNA.

KEY RESULTS

The 5S and 35S rRNA genes were separate (S-type) in Gnetum montanum, Gnetum gnemon and Welwitschia mirabilis and linked (L-type) in Ephedra altissima. There was considerable variability in 5S rDNA abundance, ranging from as few as ~4000 (W. mirabilis) to >100 000 (G. montanum) copies. A similar large variation was also observed in 5S rDNA locus numbers (two to 16 sites per diploid cell). 5S rRNA pseudogenes were interspersed between functional genes forming a single unit in E. altissima and G. montanum. Their copy number was comparable or even higher than that of functional 5S rRNA genes. In E. altissima internal transcribed spacers of 35S rDNA were long and intrinsically repetitive while in G. montanum and W. mirabilis they were short without the subrepeats.

CONCLUSIONS

Gnetophytes are distinct from other gymnosperms and angiosperms as they display surprisingly large variability in rDNA organization and rDNA copy and locus numbers between genera, with no relationship between copy numbers and genome sizes apparent. Concerted evolution of 5S rDNA units seems to have led to the amplification of 5S pseudogenes in both G. montanum and E. altissima. Evolutionary patterns of rDNA show both gymnosperm and angiosperm features underlining the diversity of the group.

Pre-labelled oligo probe-FISH karyotype analyses of four Araliaceae species using rDNA and telomeric repeat

BACKGROUND

The family Araliaceae contains many medicinal species including ginseng of which the whole genome sequencing analyses have been going on these days.

OBJECTIVE

To characterize the chromosomal distribution of 5S and 45S rDNAs and telomeric repeat in four ginseng related species of Aralia elata (Miq.) Seem., Dendropanax morbiferus H. Lév., Eleutherococcus sessiliflorus (Rupr. Et Maxim.) Seem., Kalopanax septemlobus (Thunb. ex A.Murr.) Koidz.

METHOD

Pre-labelled oligoprobe (PLOP)-fluorescence in situ hybridization (FISH) was carried out.

RESULTS

The chromosome number of A. elata was 2n = 24, whereas that of the other three species of D. morbiferus, E. sessiliflorus, and K. septemlobus was 2n = 48, corresponding to diploid and tetraploid, respectively, based on the basic chromosome number x = 12 in Araliaceae. In all four species, one pair of 5S signals were detected in the proximal regions of the short arms of chromosome 3, whereas in K. septemlobus, the 5S rDNA signals localized in the subtelomeric region of short arm of chromosome 3, while all the 45S rDNA signals localized at the paracentromeric region of the short arm of chromosome 1. And the telomeric repeat signals were detected at the telomeric region of both short and long arms of most chromosomes.

CONCLUSION

The PLOP-FISH was very effective compared with conventional FISH method. These results provide useful comparative cytogenetic information to better understand the genome structure of each species and will be useful to trace the history of ginseng genomic constitution.

Rapid and Efficient FISH using Pre-Labeled Oligomer Probes

Fluorescence in situ hybridization (FISH) is used to visualize the distribution of DNA elements within a genome. Conventional methods for FISH take 1-2 days. Here, we developed a simplified, rapid FISH technique using pre-labeled oligonucleotide probes (PLOPs) and tested the procedure using 18 PLOPs from 45S and 5S rDNA, Arabidopsis-type telomere, and newly-identified Panax ginseng-specific tandem repeats. The 16 developed rDNA PLOPs can be universally applied to plants and animals. The telomere PLOPs can be utilized in most plants with Arabidopsis-type telomeres. The ginseng-specific PLOP can be used to distinguish P. ginseng from related Panax species. Differential labeling of PLOPs allowed us to simultaneously visualize different target loci while reducing the FISH hybridization time from ~16 h to 5 min. PLOP-FISH is efficient, reliable, and rapid, making it ideal for routine analysis, especially of newly sequenced genomes using either universal or specific targets, such as novel tandem repeats identified from whole-genome sequencing data.

Cytogenetic features of rRNA genes across land plants: analysis of the Plant rDNA database

The online resource http://www.plantrdnadatabase.com/ stores information on the number, chromosomal locations and structure of the 5S and 18S-5.8S-26S (35S) ribosomal DNAs (rDNA) in plants. This resource was exploited to study relationships between rDNA locus number, distribution, the occurrence of linked (L-type) and separated (S-type) 5S and 35S rDNA units, chromosome number, genome size and ploidy level. The analyses presented summarise current knowledge on rDNA locus numbers and distribution in plants. We analysed 2949 karyotypes, from 1791 species and 86 plant families, and performed ancestral character state reconstructions. The ancestral karyotype (2n = 16) has two terminal 35S sites and two interstitial 5S sites, while the median (2n = 24) presents four terminal 35S sites and three interstitial 5S sites. Whilst 86.57% of karyotypes show S-type organisation (ancestral condition), the L-type arrangement has arisen independently several times during plant evolution. A non-terminal position of 35S rDNA was found in about 25% of single-locus karyotypes, suggesting that terminal locations are not essential for functionality and expression. Single-locus karyotypes are very common, even in polyploids. In this regard, polyploidy is followed by subsequent locus loss. This results in a decrease in locus number per monoploid genome, forming part of the diploidisation process returning polyploids to a diploid-like state over time.

A Highly Dense Genetic Map for Ginkgo biloba Constructed Using Sequence-Based Markers

Ginkgo biloba L. is a well-known living gymnosperm fossil that has medicinal and ornamental value. In this study, a high density genetic map was constructed with megagametophytes of 94 seeds from a single Ginkgo tree by employing the specific-locus amplified fragment (SLAF) sequencing technique. The average sequencing depth was 11.20×, which yielded 538,031 high-quality SLAFs. Among these SLAFs, 204,361 were heterozygous in the maternal tree and segregated in the progeny. The established map contained 12,263 SLAFs that were assigned to 12 linkage groups (LGs). The number of LGs on this map equaled the number of chromosomes in Ginkgo. The total map length was 1,671.77 cM, with an average distance of 0.89 cM between adjacent marker bins. Map evaluation based on the haplotype map and the heat map validated the high quality of the established map. Because Ginkgo is an economically and biologically important tree, strenuous efforts have been exerted to sequence its genome. This new map, built using sequence-based markers, will serve in the future as a fundamental platform for anchoring sequence assemblies along Ginkgo chromosomes. This map also provides a desirable platform for various genetic studies of Ginkgo, including gene/quantitative trait locus mapping and marker-aided selection.

A Portrait of Ribosomal DNA Contacts with Hi-C Reveals 5S and 45S rDNA Anchoring Points in the Folded Human Genome

Ribosomal RNAs (rRNAs) account for >60% of all RNAs in eukaryotic cells and are encoded in the ribosomal DNA (rDNA) arrays. The rRNAs are produced from two sets of loci: the 5S rDNA array resides exclusively on human chromosome 1, whereas the 45S rDNA array resides on the short arm of five human acrocentric chromosomes. The 45S rDNA gives origin to the nucleolus, the nuclear organelle that is the site of ribosome biogenesis. Intriguingly, 5S and 45S rDNA arrays exhibit correlated copy number variation in lymphoblastoid cells (LCLs). Here we examined the genomic architecture and repeat content of the 5S and 45S rDNA arrays in multiple human genome assemblies (including PacBio MHAP assembly) and ascertained contacts between the rDNA arrays and the rest of the genome using Hi-C datasets from two human cell lines (erythroleukemia K562 and lymphoblastoid cells). Our analyses revealed that 5S and 45S arrays each have thousands of contacts in the folded genome, with rDNA-associated regions and genes dispersed across all chromosomes. The rDNA contact map displayed conserved and disparate features between two cell lines, and pointed to specific chromosomes, genomic regions, and genes with evidence of spatial proximity to the rDNA arrays; the data also showed a lack of direct physical interaction between the 5S and 45S rDNA arrays. Finally, the analysis identified an intriguing organization in the 5S array with Alu and 5S elements adjacent to one another and organized in opposite orientation along the array. Portraits of genome folding centered on the ribosomal DNA array could help understand the emergence of concerted variation, the control of 5S and 45S expression, as well as provide insights into an organelle that contributes to the spatial localization of human chromosomes during interphase.

Astonishing 35S rDNA diversity in the gymnosperm species Cycas revoluta Thunb

In all eukaryotes, the highly repeated 35S ribosomal DNA (rDNA) sequences encoding 18S-5.8S-26S ribosomal RNA (rRNA) typically show high levels of intragenomic uniformity due to homogenisation processes, leading to concerted evolution of 35S rDNA repeats. Here, we compared 35S rDNA divergence in several seed plants using next generation sequencing and a range of molecular and cytogenetic approaches. Most species showed similar 35S rDNA homogeneity indicating concerted evolution. However, Cycas revoluta exhibits an extraordinary diversity of rDNA repeats (nucleotide sequence divergence of different copies averaging 12 %), influencing both the coding and non-coding rDNA regions nearly equally. In contrast, its rRNA transcriptome was highly homogeneous suggesting that only a minority of genes (<20 %) encode functional rRNA. The most common SNPs were C > T substitutions located in symmetrical CG and CHG contexts which were also highly methylated. Both functional genes and pseudogenes appear to cluster on chromosomes. The extraordinary high levels of 35S rDNA diversity in C. revoluta, and probably other species of cycads, indicate that the frequency of repeat homogenisation has been much lower in this lineage, compared with all other land plant lineages studied. This has led to the accumulation of methylation-driven mutations and pseudogenisation. Potentially, the reduced homology between paralogs prevented their elimination by homologous recombination, resulting in long-term retention of rDNA pseudogenes in the genome.

Complete chloroplast and ribosomal sequences for 30 accessions elucidate evolution of Oryza AA genome species

Cytoplasmic chloroplast (cp) genomes and nuclear ribosomal DNA (nR) are the primary sequences used to understand plant diversity and evolution. We introduce a high-throughput method to simultaneously obtain complete cp and nR sequences using Illumina platform whole-genome sequence. We applied the method to 30 rice specimens belonging to nine Oryza species. Concurrent phylogenomic analysis using cp and nR of several of specimens of the same Oryza AA genome species provides insight into the evolution and domestication of cultivated rice, clarifying three ambiguous but important issues in the evolution of wild Oryza species. First, cp-based trees clearly classify each lineage but can be biased by inter-subspecies cross-hybridization events during speciation. Second, O. glumaepatula, a South American wild rice, includes two cytoplasm types, one of which is derived from a recent interspecies hybridization with O. longistminata. Third, the Australian O. rufipogan-type rice is a perennial form of O. meridionalis.

Non-Random Distribution of 5S rDNA Sites and Its Association with 45S rDNA in Plant Chromosomes

5S and 45S rDNA sites are the best mapped chromosome regions in eukaryotic chromosomes. In this work, a database was built gathering information about the position and number of 5S rDNA sites in 784 plant species, aiming to identify patterns of distribution along the chromosomes and its correlation with the position of 45S rDNA sites. Data revealed that in most karyotypes (54.5%, including polyploids) two 5S rDNA sites (a single pair) are present, with 58.7% of all sites occurring in the short arm, mainly in the proximal region. In karyotypes of angiosperms with only 1 pair of sites (single sites) they are mostly found in the proximal region (52.0%), whereas in karyotypes with multiple sites the location varies according to the average chromosome size. Karyotypes with multiple sites and small chromosomes (<3 µm) often display proximal sites, while medium-sized (between 3 and 6 µm) and large chromosomes (>6 µm) more commonly show terminal or interstitial sites. In species with holokinetic chromosomes, the modal value of sites per karyotype was also 2, but they were found mainly in a terminal position. Adjacent 5S and 45S rDNA sites were often found in the short arm, reflecting the preferential distribution of both sites in this arm. The high frequency of genera with at least 1 species with adjacent 5S and 45S sites reveals that this association appeared several times during angiosperm evolution, but it has been maintained only rarely as the dominant array in plant genera.

Incomplete sequence homogenization in 45S rDNA multigene families: intermixed IGS heterogeneity within the single NOR locus of the polyploid species Medicago arborea (Fabaceae)

BACKGROUND AND AIMS

Ribosomal sequences have become the classical example of the genomic homogenization of nuclear multigene families. Despite theoretical advantages and modelling predictions that support concerted evolution of the 45S rDNA, several reports have found intragenomic polymorphisms. However, the origins and causes of these rDNA polymorphisms are difficult to assess because seed plants show a wide range of 45S rDNA loci number variation, especially in polyploids. Medicago arborea is a tetraploid species that has a single 45S rDNA locus. This feature makes this species a suitable case study to assess the fate of ribosomal IGS homogenization in polyploid species showing nucleolus organizer region (NOR) reduction.

METHODS

The intergenic spacer (IGS) region was amplified by long PCR and the fragments were cloned and sequenced by a primer-walking strategy. The physical mapping of the whole and partial IGS variants was assessed by fluorescent in situ hybridization (FISH) and fibre-FISH methods on mitotic chromosomes and extended DNA fibres, respectively.

KEY RESULTS

Two IGS fragments of 4·8 and 3·5 kb were obtained showing structural features of functional sequences. The shorter variant appears to be a truncated copy of the 4·8 kb fragment that lacks the duplication of the transcription initiation site region and the entire D region. The physical localization of the two IGS variants on metaphase chromosomes and extended DNA fibres using FISH corroborated their joint presence within the same locus. In addition, no spatial structure of the two variants was detected within the NOR.

CONCLUSIONS

The results suggest that full sequence homogenization is not operating within the NOR locus of M. arborea. The structure of the NOR locus reported here departs from the models of IGS heterogeneity present in plants and caution against assuming the widespread belief that intragenomic ribosomal heterogeneity is mainly due to sequence variation between paralogous loci.

Plant rDNA database: update and new features

The Plant rDNA database (www.plantrdnadatabase.com) is an open access online resource providing detailed information on numbers, structures and positions of 5S and 18S-5.8S-26S (35S) ribosomal DNA loci. The data have been obtained from >600 publications on plant molecular cytogenetics, mostly based on fluorescent in situ hybridization (FISH). This edition of the database contains information on 1609 species derived from 2839 records, which means an expansion of 55.76 and 94.45%, respectively. It holds the data for angiosperms, gymnosperms, bryophytes and pteridophytes available as of June 2013. Information from publications reporting data for a single rDNA (either 5S or 35S alone) and annotation regarding transcriptional activity of 35S loci now appears in the database. Preliminary analyses suggest greater variability in the number of rDNA loci in gymnosperms than in angiosperms. New applications provide ideograms of the species showing the positions of rDNA loci as well as a visual representation of their genome sizes. We have also introduced other features to boost the usability of the Web interface, such as an application for convenient data export and a new section with rDNA–FISH-related information (mostly detailing protocols and reagents). In addition, we upgraded and/or proofread tabs and links and modified the website for a more dynamic appearance. This manuscript provides a synopsis of these changes and developments.

Database URL: http://www.plantrdnadatabase.com

Dancing together and separate again: gymnosperms exhibit frequent changes of fundamental 5S and 35S rRNA gene (rDNA) organisation

In higher eukaryotes, the 5S rRNA genes occur in tandem units and are arranged either separately (S-type arrangement) or linked to other repeated genes, in most cases to rDNA locus encoding 18S–5.8S–26S genes (L-type arrangement). Here we used Southern blot hybridisation, PCR and sequencing approaches to analyse genomic organisation of rRNA genes in all large gymnosperm groups, including Coniferales, Ginkgoales, Gnetales and Cycadales. The data are provided for 27 species (21 genera). The 5S units linked to the 35S rDNA units occur in some but not all Gnetales, Coniferales and in Ginkgo (∼30% of the species analysed), while the remaining exhibit separate organisation. The linked 5S rRNA genes may occur as single-copy insertions or as short tandems embedded in the 26S–18S rDNA intergenic spacer (IGS). The 5S transcript may be encoded by the same (Ginkgo, Ephedra) or opposite (Podocarpus) DNA strand as the 18S–5.8S–26S genes. In addition, pseudogenised 5S copies were also found in some IGS types. Both L- and S-type units have been largely homogenised across the genomes. Phylogenetic relationships based on the comparison of 5S coding sequences suggest that the 5S genes independently inserted IGS at least three times in the course of gymnosperm evolution. Frequent transpositions and rearrangements of basic units indicate relatively relaxed selection pressures imposed on genomic organisation of 5S genes in plants.

New insights into the dynamics of plant cell nuclei and chromosomes

The plant lamin-like protein NMCP/AtLINC and orthologues of the SUN-KASH complex across the nuclear envelope (NE) show the universality of nuclear structure in eukaryotes. However, depletion of components in the connection complex of the NE in plants does not induce severe defects, unlike in animals. Appearance of the Rabl configuration is not dependent on genome size in plant species. Topoisomerase II and condensin II are not essential for plant chromosome condensation. Plant endoreduplication shares several common characteristics with animals, including involvement of cyclin-dependent kinases and E2F transcription factors. Recent finding regarding endomitosis regulator GIG1 shed light on the suppression mechanism of endomitosis in plants. The robustness of plants, compared with animals, is reflected in their genome redundancy. Spatiotemporal functional analyses using chromophore-assisted light inactivation, super-resolution microscopy, and 4D (3D plus time) imaging will reveal new insights into plant nuclear and chromosomal dynamics.

Expression of 5 S rRNA genes linked to 35 S rDNA in plants, their epigenetic modification and regulatory element divergence

BACKGROUND

In plants, the 5 S rRNA genes usually occur as separate tandems (S-type arrangement) or, less commonly, linked to 35 S rDNA units (L-type). The activity of linked genes remains unknown so far. We studied the homogeneity and expression of 5 S genes in several species from family Asteraceae known to contain linked 35 S-5 S units. Additionally, their methylation status was determined using bisulfite sequencing. Fluorescence in situ hybridization was applied to reveal the sub-nuclear positions of rDNA arrays.

RESULTS

We found that homogenization of L-type units went to completion in most (4/6) but not all species. Two species contained major L-type and minor S-type units (termed L(s)-type). The linked genes dominate 5 S rDNA expression while the separate tandems do not seem to be expressed. Members of tribe Anthemideae evolved functional variants of the polymerase III promoter in which a residing C-box element differs from the canonical angiosperm motif by as much as 30%. On this basis, a more relaxed consensus sequence of a plant C-box: (5'-RGSWTGGGTG-3') is proposed. The 5 S paralogs display heavy DNA methylation similarly as to their unlinked counterparts. FISH revealed the close association of 35 S-5 S arrays with nucleolar periphery indicating that transcription of 5 S genes may occur in this territory.

CONCLUSIONS

We show that the unusual linked arrangement of 5 S genes, occurring in several plant species, is fully compatible with their expression and functionality. This extraordinary 5 S gene dynamics is manifested at different levels, such as variation in intrachromosomal positions, unit structure, epigenetic modification and considerable divergence of regulatory motifs.