Selective nucleolus organizer inactivation in Arabidopsis is a chromosome position-effect phenomenon

Cyto-molecular characterization of rDNA and chromatin composition in the NOR-associated satellite in Chestnut (Castanea spp.)

The American chestnut (Castanea dentata, 2n = 2x = 24), once known as the "King of the Appalachian Forest", was decimated by chestnut blight during the first half of the twentieth century by an invasive fungus (Cryphonectria parasitica). The Chinese chestnut (C. mollissima, 2n = 2x = 24), in contrast to American chestnut, is resistant to this blight. Efforts are being made to transfer this resistance to American chestnut through backcross breeding and genetic engineering. Both chestnut genomes have been genetically mapped and recently sequenced to facilitate gene discovery efforts aimed at assisting molecular breeding and genetic engineering. To complement and extend this genomic work, we analyzed the distribution and organization of their ribosomal DNAs (35S and 5S rDNA), and the chromatin composition of the nucleolus organizing region (NOR)-associated satellites. Using fluorescent in situ hybridization (FISH), we have identified two 35S (one major and one minor) and one 5S rDNA sites. The major 35S rDNA sites are terminal and sub-terminal in American and Chinese chestnuts, respectively, originating at the end of the short arm of the chromosome, extending through the secondary constriction and into the satellites. An additional 5S locus was identified in certain Chinese chestnut accessions, and it was linked distally to the major 35S site. The NOR-associated satellite in Chinese chestnut was found to comprise a proximal region packed with 35S rDNA and a distinct distal heterochromatic region. In contrast, the American chestnut satellite was relatively small and devoid of the distal heterochromatic region.

Low dose ribosomal DNA P-loop mutation affects development and enforces autophagy in Arabidopsis

Arabidopsis contains hundreds of ribosomal DNA copies organized within the nucleolar organizing regions (NORs) in chromosomes 2 and 4. There are four major types of variants of rDNA, VAR1-4, based on the polymorphisms of 3' external transcribed sequences. The variants are known to be differentially expressed during plant development. We created a mutant by the CRISPR-Cas9-mediated excision of ~ 25 nt from predominantly NOR4 ribosomal DNA copies, obtaining mosaic mutational events on ~ 5% of all rDNA copies. The excised region consists of P-loop and Helix-82 segments of 25S rRNA. The mutation led to allelic, dosage-dependent defects marked by lateral root inhibition, reduced size, and pointy leaves, all previously observed for defective ribosomal function. The mutation in NOR4 led to dosage compensation from the NOR2 copies by elevated expression of VAR1 in mutants and further associated single-nucleotide variants, thus, resulting in altered rRNA sub-population. Furthermore, the mutants exhibited rRNA maturation defects specifically in the minor pathway typified by 32S pre-rRNA accumulation. Density-gradient fractionation and subsequent RT-PCR of rRNA analyses revealed that mutated copies were not incorporated into the translating ribosomes. The mutants in addition displayed an elevated autophagic flux as shown by the autophagic marker GFP-ATG8e, likely related to ribophagy.

Downregulation of ribosomal RNA (rRNA) genes in human head and neck squamous cell carcinoma (HNSCC) cells correlates with rDNA promoter hypermethylation

Eukaryotes carry hundreds of ribosomal RNA (rRNA) genes as tandem arrays, which generate rRNA for protein synthesis. Humans carry ∼ 400 rRNA gene copies and their expression is epigenetically regulated. Dysregulation of rRNA synthesis and ribosome biogenesis are characteristic features of cancers. Targeting aberrant rRNA expression for cancer therapy is being explored. Head and neck squamous cell carcinoma (HNSCC) is among the most prevalent cancers globally. Using quantitative PCR and bisulfite sequencing, we show that rRNA genes are downregulated and their promoters are hypermethylated in HNSCC cell lines. These findings may have relevance for prognosis and diagnosis of HNSCC.

Sequence and epigenetic landscapes of active and silent nucleolus organizer regions in Arabidopsis

Arabidopsis thaliana has two ribosomal RNA (rRNA) gene loci, nucleolus organizer regions NOR2 and NOR4, whose complete sequences are missing in current genome assemblies. Ultralong DNA sequences assembled using an unconventional approach yielded ~5.5- and 3.9-Mbp sequences for NOR2 and NOR4 in the reference strain, Col-0. The distinct rRNA gene subtype compositions of the NORs enabled the positional mapping of their active and inactive regions, using RNA sequencing to identify subtype-specific transcripts and DNA sequencing to identify subtypes associated with flow-sorted nucleoli. Comparisons of wild-type and silencing-defective plants revealed that most rRNA gene activity occurs in the central region of NOR4, whereas most, but not all, genes of NOR2 are epigenetically silenced. Intervals of low CG and CHG methylation overlap regions where gene activity and gene subtype homogenization are high. Collectively, the data reveal the genetic and epigenetic landscapes underlying nucleolar dominance (differential NOR activity) and implicate transcription as a driver of rRNA gene concerted evolution.

Uniparental expression of ribosomal RNA in ×Festulolium grasses: a link between the genome and nucleolar dominance

Genome or genomic dominance (GD) is a phenomenon observed in hybrids when one parental genome becomes dominant over the other. It is manifested by the replacement of chromatin of the submissive genome by that of the dominant genome and by biased gene expression. Nucleolar dominance (ND) - the functional expression of only one parental set of ribosomal genes in hybrids - is another example of an intragenomic competitive process which, however, concerns ribosomal DNA only. Although GD and ND are relatively well understood, the nature and extent of their potential interdependence is mostly unknown. Here, we ask whether hybrids showing GD also exhibit ND and, if so, whether the dominant genome is the same. To test this, we used hybrids between Festuca and Lolium grasses (Festulolium), and between two Festuca species in which GD has been observed (with Lolium as the dominant genome in Festulolium and F. pratensis in interspecific Festuca hybrids). Using amplicon sequencing of ITS1 and ITS2 of the 45S ribosomal DNA (rDNA) cluster and molecular cytogenetics, we studied the organization and expression of rDNA in leaf tissue in five hybrid combinations, four generations and 31 genotypes [F. pratensis × L. multiflorum (F1, F2, F3, BC1), L. multiflorum × F. pratensis (F1), L. multiflorum × F. glaucescens (F2), L. perenne × F. pratensis (F1), F. glaucescens × F. pratensis (F1)]. We have found that instant ND occurs in Festulolium, where expression of Lolium-type rDNA reached nearly 100% in all F1 hybrids and was maintained through subsequent generations. Therefore, ND and GD in Festulolium are manifested by the same dominant genome (Lolium). We also confirmed the concordance between GD and ND in an interspecific cross between two Festuca species.

Reaching for the off switch in nucleolar dominance

Nucleolus organizer regions (NORs) are eukaryotic chromosomal loci where ribosomal RNA (rRNA) genes are clustered, typically in hundreds to thousands of copies. Transcription of these rRNA genes by RNA polymerase I and processing of their transcripts results in the formation of the nucleolus, the sub-nuclear domain in which ribosomes are assembled. Approximately 90 years ago, cytogenetic observations revealed that NORs inherited from the different parents of an interspecific hybrid sometimes differ in morphology at metaphase. Fifty years ago, those chromosomal differences were found to correlate with differences in rRNA gene transcription and the phenomenon became known as nucleolar dominance. Studies of the past 30 years have revealed that nucleolar dominance results from selective rRNA gene silencing, involving repressive chromatin modifications, and occurs in pure species as well as hybrids. Recent evidence also indicates that silencing depends on the NOR in which an rRNA gene is located, and not on the gene's sequence. In this perspective, we discuss how our thinking about nucleolar dominance has shifted over time from the kilobase scale of individual genes to the megabase scale of NORs and chromosomes and questions that remain unanswered in the search for a genetic and biochemical understanding of the off switch.

NICOTIANAMINE SYNTHASE activity affects nucleolar iron accumulation and impacts rDNA silencing and RNA methylation in Arabidopsis

In plant cells, a large pool of iron (Fe) is contained in the nucleolus, as well as in chloroplasts and mitochondria. A central determinant for intracellular distribution of Fe is nicotianamine (NA) generated by NICOTIANAMINE SYNTHASE (NAS). Here, we used Arabidopsis thaliana plants with disrupted NAS genes to study the accumulation of nucleolar iron and understand its role in nucleolar functions and more specifically in rRNA gene expression. We found that nas124 triple mutant plants, which contained lower quantities of the iron ligand NA, also contained less iron in the nucleolus. This was concurrent with the expression of normally silenced rRNA genes from nucleolar organizer regions 2 (NOR2). Notably, in nas234 triple mutant plants, which also contained lower quantities of NA, nucleolar iron and rDNA expression were not affected. In contrast, in both nas124 and nas234, specific RNA modifications were differentially regulated in a genotype dependent manner. Taken together, our results highlight the impact of specific NAS activities in RNA gene expression. We discuss the interplay between NA and nucleolar iron with rDNA functional organization and RNA methylation.

Switch them off or not: selective rRNA gene repression in grasses

Nucleolar dominance (ND) is selective epigenetic silencing of 35-48S rDNA loci. In allopolyploids, it is frequently manifested at the cytogenetic level by the inactivation of nucleolar organiser region(s) (NORs) inherited from one or several evolutionary ancestors. Grasses are ecologically and economically one of the most important land plant groups, which have frequently evolved through hybridisation and polyploidisation events. Here we review common and unique features of ND phenomena in this monocot family from cytogenetic, molecular, and genomic perspectives. We highlight recent advances achieved by using an allotetraploid model grass, Brachypodium hybridum, where ND commonly occurs at a population level, and we cover modern genomic approaches that decipher structural features of core arrays of NORs.

Structural variation among assembled genomes facilitates development of rapid and low-cost NOR-linked markers and NOR-telomere junction mapping in Arabidopsis

KEY MESSAGE

Genome-wide structural variants we identified and new NOR-linked markers we developed would be useful for future genome-wide association studies (GWAS), and for new gene/trait mapping purposes. Bioinformatic alignment of the assembled genomes of Col-0 and Sha ecotypes of Arabidopsis thaliana revealed ~ 13,000 genome-wide structural variants involving simple insertions or deletions and repeat contractions or expansions. Using some of these structural variants, we developed new, rapid, and low-cost PCR-based molecular markers that are genetically linked to the nucleolus organizer regions (NORs). A. thaliana has two NORs, one each on chromosome 2 (NOR2) and chromosome 4 (NOR4). Both NORs are ~ 4 Mb each, and hundreds of 45S ribosomal RNA (rRNA) genes are tandemly arrayed at these loci. Using previously characterized recombinant inbred lines (RILs) derived from Sha x Col-0  crosses, we validated the utility of the newly developed NOR-linked markers in genetically mapping rRNA genes and the associated telomeres to either NOR2 or NOR4. Lastly, we sequenced Sha genome using Oxford Nanopore Technology (ONT) and used the data to obtain sequences of NOR-telomere junctions, and with the help of RILs, we mapped them as new genetic markers to their respective NORs (NOR2-TEL2N and NOR4-TEL4N). The structural variants obtained from this study would serve as valuable data for genome-wide association studies (GWAS), and to rapidly design more genome-wide genetic (molecular) markers for new gene/trait mapping purposes.

Under the magnifying glass: The ups and downs of rDNA copy number

The ribosomal DNA (rDNA) in Drosophila is found as two additive clusters of individual 35 S cistrons. The multiplicity of rDNA is essential to assure proper translational demands, but the nature of the tandem arrays expose them to copy number variation within and between populations. Here, we discuss means by which a cell responds to insufficient rDNA copy number, including a historical view of rDNA magnification whose mechanism was inferred some 35 years ago. Recent work has revealed that multiple conditions may also result in rDNA loss, in response to which rDNA magnification may have evolved. We discuss potential models for the mechanism of magnification, and evaluate possible consequences of rDNA copy number variation.

Structure and Methylation of 35S rDNA in Allopolyploids Anemone multifida (2n = 4x = 32, BBDD) and Anemone baldensis (2n = 6x = 48, AABBDD) and Their Parental Species Show Evidence of Nucleolar Dominance

Transcriptional silencing of 35S rDNA loci inherited from one parental species is occurring relatively frequently in allopolyploids. However, molecular mechanisms by which it is selected for transcriptional silencing remain unclear. We applied NGS, silver staining and bisulfite sequencing to study the structure, expression and methylation landscape of 35S rDNA in two allopolyploids of common origin, allotetraploid Anemone multifida (2n = 4x = 32, genome composition BBDD) and allohexaploid A. baldensis (2n = 6x = 48, AABBDD), and their genome donors, A. sylvestris (2n = 16, AA), A. cylindrica (2n = 16, BB) and A. parviflora (2n = 16, DD). The size of the recovered 35S rDNA units varied from 10,489 bp in A. cylindrica to 12,084 bp in A. sylvestris. Anemone showed an organization typical of most ribosomal 35S rDNA composed of NTS, ETS, rRNA genes, TTS and TIS with structural features of plant IGS sequences and all functional elements needed for rRNA gene activity. The NTS was more variable than the ETS and consisted of SRs which are highly variable among Anemone. Five to six CpG-rich islands were found within the ETS. CpG island located adjacent to the transcription initiation site (TIS) was highly variable regarding the sequence size and methylation level and exhibited in most of the species lower levels of methylation than CpG islands located adjacent to the 18S rRNA gene. Our results uncover hypomethylation of A. sylvestris- and A. parviflora-derived 35S rDNA units in allopolyploids A. multifida and A. baldensis. Hypomethylation of A. parviflora-derived 35S rDNA was more prominent in A. baldensis than in A. multifida. We showed that A. baldensis underwent coupled A. sylvestris-derived 35S rDNA array expansion and A. parviflora-derived 35S rDNA copy number decrease that was accompanied by lower methylation level of A. sylvestris-derived 35S rDNA units in comparison to A. parviflora-derived 35S rDNA units. These observations suggest that in A. baldensis nucleolar dominance is directed toward A. sylvestris-derived chromosomes. This work broadens our current knowledge of the 35S rDNA organization in Anemone and provides evidence of the progenitor-specific 35S rDNA methylation in nucleolar dominance.

Tracing the Evolution of the Angiosperm Genome from the Cytogenetic Point of View

Cytogenetics constitutes a branch of genetics that is focused on the cellular components, especially chromosomes, in relation to heredity and genome structure, function and evolution. The use of modern cytogenetic approaches and the latest microscopes with image acquisition and processing systems enables the simultaneous two- or three-dimensional, multicolour visualisation of both single-copy and highly-repetitive sequences in the plant genome. The data that is gathered using the cytogenetic methods in the phylogenetic background enable tracing the evolution of the plant genome that involve changes in: (i) genome sizes; (ii) chromosome numbers and morphology; (iii) the content of repetitive sequences and (iv) ploidy level. Modern cytogenetic approaches such as FISH using chromosome- and genome-specific probes have been widely used in studies of the evolution of diploids and the consequences of polyploidy. Nowadays, modern cytogenetics complements analyses in other fields of cell biology and constitutes the linkage between genetics, molecular biology and genomics.

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."

Modeling the 3D genome of plants

Chromosomes are the carriers of inheritable traits and define cell function and development. This is not only based on the linear DNA sequence of chromosomes but also on the additional molecular information they are associated with, including the transcription machinery, histone modifications, and their three-dimensional folding. The synergistic application of experimental approaches and computer simulations has helped to unveil how these organizational layers of the genome interplay in various organisms. However, such multidisciplinary approaches are still rarely explored in the plant kingdom. Here, we provide an overview of our current knowledge on plant 3D genome organization and review recent efforts to integrate cutting-edge experiments from microscopy and next-generation sequencing approaches with theoretical models. Building on these recent approaches, we propose possible avenues to extend the application of theoretical modeling in the characterization of the 3D genome organization in plants.

Genetic, Epigenetic, Genomic and Microbial Approaches to Enhance Salt Tolerance of Plants: A Comprehensive Review

Simple Summary

Globally, soil salinity, which refers to salt-affected soils, is increasing due to various environmental factors and human activities. Soil salinity poses one of the most serious challenges in the field of agriculture as it significantly reduces the growth and yield of crop plants, both quantitatively and qualitatively. Over the last few decades, several studies have been carried out to understand plant biology in response to soil salinity stress with a major emphasis on genetic and other hereditary components. Based on the outcome of these studies, several approaches are being followed to enhance plants’ ability to tolerate salt stress while still maintaining reasonable levels of crop yields. In this manuscript, we comprehensively list and discuss various biological approaches being followed and, based on the recent advances in the field of molecular biology, we propose some new approaches to improve salinity tolerance of crop plants. The global scientific community can make use of this information for the betterment of crop plants. This review also highlights the importance of maintaining global soil health to prevent several crop plant losses.

Abstract

Globally, soil salinity has been on the rise owing to various factors that are both human and environmental. The abiotic stress caused by soil salinity has become one of the most damaging abiotic stresses faced by crop plants, resulting in significant yield losses. Salt stress induces physiological and morphological modifications in plants as a result of significant changes in gene expression patterns and signal transduction cascades. In this comprehensive review, with a major focus on recent advances in the field of plant molecular biology, we discuss several approaches to enhance salinity tolerance in plants comprising various classical and advanced genetic and genetic engineering approaches, genomics and genome editing technologies, and plant growth-promoting rhizobacteria (PGPR)-based approaches. Furthermore, based on recent advances in the field of epigenetics, we propose novel approaches to create and exploit heritable genome-wide epigenetic variation in crop plants to enhance salinity tolerance. Specifically, we describe the concepts and the underlying principles of epigenetic recombinant inbred lines (epiRILs) and other epigenetic variants and methods to generate them. The proposed epigenetic approaches also have the potential to create additional genetic variation by modulating meiotic crossover frequency.

Analyses of the Updated "Animal rDNA Loci Database" with an Emphasis on Its New Features

We report on a major update to the animal rDNA loci database, which now contains cytogenetic information for 45S and 5S rDNA loci in more than 2600 and 1000 species, respectively. The data analyses show the following: (i) A high variability in 5S and 45S loci numbers, with both showing 50-fold or higher variability. However, karyotypes with an extremely high number of loci were rare, and medians generally converged to two 5S sites and two 45S rDNA sites per diploid genome. No relationship was observed between the number of 5S and 45S loci. (ii) The position of 45S rDNA on sex chromosomes was relatively frequent in some groups, particularly in arthropods (14% of karyotypes). Furthermore, 45S rDNA was almost exclusively located in microchromosomes when these were present (in birds and reptiles). (iii) The proportion of active NORs (positively stained with silver staining methods) progressively decreased with an increasing number of 45S rDNA loci, and karyotypes with more than 12 loci showed, on average, less than 40% of active loci. In conclusion, the updated version of the database provides some new insights into the organization of rRNA genes in chromosomes. We expect that its updated content will be useful for taxonomists, comparative cytogeneticists, and evolutionary biologists.

The integrative biology of genetic dominance

Dominance is a basic property of inheritance systems describing the link between a diploid genotype at a single locus and the resulting phenotype. Models for the evolution of dominance have long been framed as an opposition between the irreconcilable views of Fisher in 1928 supporting the role of largely elusive dominance modifiers and Wright in 1929, who viewed dominance as an emerging property of the structure of enzymatic pathways. Recent theoretical and empirical advances however suggest that these opposing views can be reconciled, notably using models investigating the regulation of gene expression and developmental processes. In this more comprehensive framework, phenotypic dominance emerges from departures from linearity between any levels of integration in the genotype‐to‐phenotype map. Here, we review how these different models illuminate the emergence and evolution of dominance. We then detail recent empirical studies shedding new light on the diversity of molecular and physiological mechanisms underlying dominance and its evolution. By reconciling population genetics and functional biology, we hope our review will facilitate cross‐talk among research fields in the integrative study of dominance evolution.

Gene dosage compensation of rRNA transcript levels in Arabidopsis thaliana lines with reduced ribosomal gene copy number

The 45S rRNA genes (rDNA) are among the largest repetitive elements in eukaryotic genomes. rDNA consists of tandem arrays of rRNA genes, many of which are transcriptionally silenced. Silent rDNA repeats may act as 'back-up' copies for ribosome biogenesis and have nuclear organization roles. Through Cas9-mediated genome editing in the Arabidopsis thaliana female gametophyte, we reduced 45S rDNA copy number (CN) to a plateau of ∼10%. Two independent lines had rDNA CNs reduced by up to 90% at the T7 generation, named low copy number (LCN) lines. Despite drastic reduction of rDNA copies, rRNA transcriptional rates, and steady-state levels remained the same as wild-type plants. Gene dosage compensation of rRNA transcript levels was associated with reduction of silencing histone marks at rDNA loci and altered Nucleolar Organiser Region 2 organization. Although overall genome integrity of LCN lines appears unaffected, a chromosome segmental duplication occurred in one of the lines. Transcriptome analysis of LCN seedlings identified several shared dysregulated genes and pathways in both independent lines. Cas9 genome editing of rRNA repeats to generate LCN lines provides a powerful technique to elucidate rDNA dosage compensation mechanisms and impacts of low rDNA CN on genome stability, development, and cellular processes.

Targeted Enrichment of rRNA Gene Tandem Arrays for Ultra-Long Sequencing by Selective Restriction Endonuclease Digestion

Large regions of nearly identical repeats, such as the 45S ribosomal RNA (rRNA) genes of Nucleolus Organizer Regions (NORs), can account for major gaps in sequenced genomes. To assemble these regions, ultra-long sequencing reads that span multiple repeats have the potential to reveal sets of repeats that collectively have sufficient sequence variation to unambiguously define that interval and recognize overlapping reads. Because individual repetitive loci typically represent a small proportion of the genome, methods to enrich for the regions of interest are desirable. Here we describe a simple method that achieves greater than tenfold enrichment of Arabidopsis thaliana 45S rRNA gene sequences among ultra-long Oxford Nanopore Technology sequencing reads. This method employs agarose-embedded genomic DNA that is subjected to restriction endonucleases digestion using a cocktail of enzymes predicted to be non-cutters of rRNA genes. Most of the genome is digested into small fragments that diffuse out of the agar plugs, whereas rRNA gene arrays are retained. In principle, the approach can also be adapted for sequencing other repetitive loci for which gaps exist in a reference genome.

Nucleolar rDNA folds into condensed foci with a specific combination of epigenetic marks

Arabidopsis thaliana 45S ribosomal genes (rDNA) are located in tandem arrays called nucleolus organizing regions on the termini of chromosomes 2 and 4 (NOR2 and NOR4) and encode rRNA, a crucial structural element of the ribosome. The current model of rDNA organization suggests that inactive rRNA genes accumulate in the condensed chromocenters in the nucleus and at the nucleolar periphery, while the nucleolus delineates active genes. We challenge the perspective that all intranucleolar rDNA is active by showing that a subset of nucleolar rDNA assembles into condensed foci marked by H3.1 and H3.3 histones that also contain the repressive H3K9me2 histone mark. By using plant lines containing a low number of rDNA copies, we further found that the condensed foci relate to the folding of rDNA, which appears to be a common mechanism of rDNA regulation inside the nucleolus. The H3K9me2 histone mark found in condensed foci represents a typical modification of bulk inactive rDNA, as we show by genome-wide approaches, similar to the H2A.W histone variant. The euchromatin histone marks H3K27me3 and H3K4me3, in contrast, do not colocalize with nucleolar foci and their overall levels in the nucleolus are very low. We further demonstrate that the rDNA promoter is an important regulatory region of the rDNA, where the distribution of histone variants and histone modifications are modulated in response to rDNA activity.

The rDNA Loci-Intersections of Replication, Transcription, and Repair Pathways

Genes encoding ribosomal RNA (rDNA) are essential for cell survival and are particularly sensitive to factors leading to genomic instability. Their repetitive character makes them prone to inappropriate recombinational events arising from collision of transcriptional and replication machineries, resulting in unstable rDNA copy numbers. In this review, we summarize current knowledge on the structure and organization of rDNA, its role in sensing changes in the genome, and its linkage to aging. We also review recent findings on the main factors involved in chromatin assembly and DNA repair in the maintenance of rDNA stability in the model plants Arabidopsis thaliana and the moss Physcomitrella patens, providing a view across the plant evolutionary tree.

Sequencing of the Arabidopsis NOR2 reveals its distinct organization and tissue-specific rRNA ribosomal variants

Despite vast differences between organisms, some characteristics of their genomes are conserved, such as the nucleolus organizing region (NOR). The NOR is constituted of multiple, highly repetitive rDNA genes, encoding the catalytic ribosomal core RNAs which are transcribed from 45S rDNA units. Their precise sequence information and organization remain uncharacterized. Here, using a combination of long- and short-read sequencing technologies we assemble contigs of the Arabidopsis NOR2 rDNA domain. We identify several expressed rRNA gene variants which are integrated into translating ribosomes in a tissue-specific manner. These findings support the concept of tissue specific ribosome subpopulations that differ in their rRNA composition and provide insights into the higher order organization of NOR2.

The nucleolus organizing region (NOR) consists of multiple, highly repetitive rDNA genes. Here Sims et al. use both long- and short-read sequencing to determine the organization and sequence of Arabidopsis NOR2 rDNA and show that different rRNA gene variants are integrated into translating ribosomes in a tissue-specific manner.

To Be or Not to Be Expressed: The First Evidence of a Nucleolar Dominance Tissue-Specificity in Brachypodium hybridum

Nucleolar dominance (ND) is an epigenetic, developmentally regulated phenomenon that describes the selective inactivation of 35S rDNA loci derived from one progenitor of a hybrid or allopolyploid. The presence of ND was documented in an allotetraploid grass, Brachypodium hybridum (genome composition DDSS), which is a polyphyletic species that arose from crosses between two putative ancestors that resembled the modern B. distachyon (DD) and B. stacei (SS). In this work, we investigated the developmental stability of ND in B. hybridum genotype 3-7-2 and compared it with the reference genotype ABR113. We addressed the question of whether the ND is established in generative tissues such as pollen mother cells (PMC). We examined condensation of rDNA chromatin by fluorescence in situ hybridization employing state-of-art confocal microscopy. The transcription of rDNA homeologs was determined by reverse-transcription cleaved amplified polymorphic sequence analysis. In ABR113, the ND was stable in all tissues analyzed (primary and adventitious root, leaf, and spikes). In contrast, the 3-7-2 individuals showed a strong upregulation of the S-genome units in adventitious roots but not in other tissues. Microscopic analysis of the 3-7-2 PMCs revealed extensive decondensation of the D-genome loci and their association with the nucleolus in meiosis. As opposed, the S-genome loci were always highly condensed and localized outside the nucleolus. These results indicate that genotype-specific loss of ND in B. hybridum occurs probably after fertilization during developmental processes. This finding supports our view that B. hybridum is an attractive model to study ND in grasses.

G4 Structures in Control of Replication and Transcription of rRNA Genes

Genes encoding 45S ribosomal RNA (rDNA) are known for their abundance within eukaryotic genomes and for their unstable copy numbers in response to changes in various genetic and epigenetic factors. Commonly, we understand as epigenetic factors (affecting gene expression without a change in DNA sequence), namely DNA methylation, histone posttranslational modifications, histone variants, RNA interference, nucleosome remodeling and assembly, and chromosome position effect. All these were actually shown to affect activity and stability of rDNA. Here, we focus on another phenomenon - the potential of DNA containing shortly spaced oligo-guanine tracts to form quadruplex structures (G4). Interestingly, sites with a high propensity to form G4 were described in yeast, animal, and plant rDNAs, in addition to G4 at telomeres, some gene promoters, and transposons, suggesting the evolutionary ancient origin of G4 as a regulatory module. Here, we present examples of rDNA promoter regions with extremely high potential to form G4 in two model plants, Arabidopsis thaliana and Physcomitrella patens. The high G4 potential is balanced by the activity of G4-resolving enzymes. The ability of rDNA to undergo these "structural gymnastics" thus represents another layer of the rich repertoire of epigenetic regulations, which is pronounced in rDNA due to its highly repetitive character.

Roles of ASYMMETRIC LEAVES2 (AS2) and Nucleolar Proteins in the Adaxial-Abaxial Polarity Specification at the Perinucleolar Region in Arabidopsis

Leaves of Arabidopsis develop from a shoot apical meristem grow along three (proximal-distal, adaxial-abaxial, and medial-lateral) axes and form a flat symmetric architecture. ASYMMETRIC LEAVES2 (AS2), a key regulator for leaf adaxial-abaxial partitioning, encodes a plant-specific nuclear protein and directly represses the abaxial-determining gene ETTIN/AUXIN RESPONSE FACTOR3 (ETT/ARF3). How AS2 could act as a critical regulator, however, has yet to be demonstrated, although it might play an epigenetic role. Here, we summarize the current understandings of the genetic, molecular, and cellular functions of AS2. A characteristic genetic feature of AS2 is the presence of a number of (about 60) modifier genes, mutations of which enhance the leaf abnormalities of as2. Although genes for proteins that are involved in diverse cellular processes are known as modifiers, it has recently become clear that many modifier proteins, such as NUCLEOLIN1 (NUC1) and RNA HELICASE10 (RH10), are localized in the nucleolus. Some modifiers including ribosomal proteins are also members of the small subunit processome (SSUP). In addition, AS2 forms perinucleolar bodies partially colocalizing with chromocenters that include the condensed inactive 45S ribosomal RNA genes. AS2 participates in maintaining CpG methylation in specific exons of ETT/ARF3. NUC1 and RH10 genes are also involved in maintaining the CpG methylation levels and repressing ETT/ARF3 transcript levels. AS2 and nucleolus-localizing modifiers might cooperatively repress ETT/ARF3 to develop symmetric flat leaves. These results raise the possibility of a nucleolus-related epigenetic repression system operating for developmental genes unique to plants and predict that AS2 could be a molecule with novel functions that cannot be explained by the conventional concept of transcription factors.

The fate of 35S rRNA genes in the allotetraploid grass Brachypodium hybridum

Nucleolar dominance (ND) consists of the reversible silencing of 35S/45S rDNA loci inherited from one of the ancestors of an allopolyploid. The molecular mechanisms by which one ancestral rDNA set is selected for silencing remain unclear. We applied a combination of molecular (Southern blot hybridization and reverse-transcription cleaved amplified polymorphic sequence analysis), genomic (analysis of variants) and cytogenetic (fluorescence in situ hybridization) approaches to study the structure, expression and epigenetic landscape of 35S rDNA in an allotetraploid grass that exhibits ND, Brachypodium hybridum (genome composition DDSS), and its putative progenitors, Brachypodium distachyon (DD) and Brachypodium stacei (SS). In progenitor genomes, B. stacei showed a higher intragenomic heterogeneity of rDNA compared with B. distachyon. In all studied accessions of B. hybridum, there was a reduction in the copy number of S homoeologues, which was accompanied by their inactive transcriptional status. The involvement of DNA methylation in CG and CHG contexts in the silencing of the S-genome rDNA loci was revealed. In the B. hybridum allotetraploid, ND is stabilized towards the D-genome units, irrespective of the polyphyletic origin of the species, and does not seem to be influenced by homoeologous 35S rDNA ratios and developmental stage.

Systematic Review of Plant Ribosome Heterogeneity and Specialization

Plants dedicate a high amount of energy and resources to the production of ribosomes. Historically, these multi-protein ribosome complexes have been considered static protein synthesis machines that are not subject to extensive regulation but only read mRNA and produce polypeptides accordingly. New and increasing evidence across various model organisms demonstrated the heterogeneous nature of ribosomes. This heterogeneity can constitute specialized ribosomes that regulate mRNA translation and control protein synthesis. A prominent example of ribosome heterogeneity is seen in the model plant, Arabidopsis thaliana, which, due to genome duplications, has multiple paralogs of each ribosomal protein (RP) gene. We support the notion of plant evolution directing high RP paralog divergence toward functional heterogeneity, underpinned in part by a vast resource of ribosome mutants that suggest specialization extends beyond the pleiotropic effects of single structural RPs or RP paralogs. Thus, Arabidopsis is a highly suitable model to study this phenomenon. Arabidopsis enables reverse genetics approaches that could provide evidence of ribosome specialization. In this review, we critically assess evidence of plant ribosome specialization and highlight steps along ribosome biogenesis in which heterogeneity may arise, filling the knowledge gaps in plant science by providing advanced insights from the human or yeast fields. We propose a data analysis pipeline that infers the heterogeneity of ribosome complexes and deviations from canonical structural compositions linked to stress events. This analysis pipeline can be extrapolated and enhanced by combination with other high-throughput methodologies, such as proteomics. Technologies, such as kinetic mass spectrometry and ribosome profiling, will be necessary to resolve the temporal and spatial aspects of translational regulation while the functional features of ribosomal subpopulations will become clear with the combination of reverse genetics and systems biology approaches.

Does the chromosomal position of 35S rDNA sites influence their transcription? A survey on Nothoscordum species (Amaryllidaceae)

35S ribosomal DNA (rDNA) sites are the regions where the ribosomal genes 18S, 5.8S and 25S, responsible for the formation of the nucleoli, are found. The fact that rDNA sites have non-random distribution on chromosomes suggests that their positions may influence their transcription. To identify if the preferentially transcribed rDNA sites occupy specific position, six species (nine cytotypes) of the genus Nothoscordum were analyzed using two different techniques to impregnate the nucleolar organizer regions (NORs) with silver nitrate. Both techniques strongly stained NORs, but one of them also stained the proximal region of all chromosomes, suggesting the existence of another group of argentophilic proteins in this region. In species with rDNA sites in acrocentric and metacentric chromosomes, sites located on the short arms of the acrocentric chromosomes were preferentially activated. On the other hand, in species with rDNA sites restricted to the short arms of the acrocentrics, all of them were activated, whereas in those species with sites restricted to the terminal region of metacentric chromosomes, the frequency of active sites was always lower than expected. This indicate that, at least in Nothoscordum, the transcription of an rDNA site is influenced by its chromosomal position, and may explain, at least partially, the strongly non-random distribution of these sites in plant and animal chromosomes.

Regulation of Nucleolar Dominance in Drosophila melanogaster

In eukaryotic genomes, ribosomal RNA (rRNA) genes exist as tandemly repeated clusters, forming ribosomal DNA (rDNA) loci. Each rDNA locus typically contains hundreds of rRNA genes to meet the high demand of ribosome biogenesis. Nucleolar dominance is a phenomenon whereby individual rDNA loci are entirely silenced or transcribed, and is believed to be a mechanism to control rRNA dosage. Nucleolar dominance was originally noted to occur in interspecies hybrids, and has been shown to occur within a species (i.e., nonhybrid context). However, studying nucleolar dominance within a species has been challenging due to the highly homogenous sequence across rDNA loci. By utilizing single nucleotide polymorphisms between X rDNA and Y rDNA loci in males, as well as sequence variations between two X rDNA loci in females, we conducted a thorough characterization of nucleolar dominance throughout development of Drosophila melanogaster. We demonstrate that nucleolar dominance is a developmentally regulated program that occurs in nonhybrid, wild-type D. melanogaster, where Y rDNA dominance is established during male embryogenesis, whereas females normally do not exhibit dominance between two X rDNA loci. By utilizing various chromosomal complements (e.g., X/Y, X/X, X/X/Y) and a chromosome rearrangement, we show that the short arm of the Y chromosome including the Y rDNA likely contains information that instructs the state of nucleolar dominance. Our study begins to reveal the mechanisms underlying the selection of rDNA loci for activation/silencing in nucleolar dominance in the context of nonhybrid D. melanogaster.

Competition of Parental Genomes in Plant Hybrids

Interspecific hybridization represents one of the main mechanisms of plant speciation. Merging of two genomes from different subspecies, species, or even genera is frequently accompanied by whole-genome duplication (WGD). Besides its evolutionary role, interspecific hybridization has also been successfully implemented in multiple breeding programs. Interspecific hybrids combine agronomic traits of two crop species or can be used to introgress specific loci of interests, such as those for resistance against abiotic or biotic stresses. The genomes of newly established interspecific hybrids (both allopolyploids and homoploids) undergo dramatic changes, including chromosome rearrangements, amplifications of tandem repeats, activation of mobile repetitive elements, and gene expression modifications. To ensure genome stability and proper transmission of chromosomes from both parental genomes into subsequent generations, allopolyploids often evolve mechanisms regulating chromosome pairing. Such regulatory systems allow only pairing of homologous chromosomes and hamper pairing of homoeologs. Despite such regulatory systems, several hybrid examples with frequent homoeologous chromosome pairing have been reported. These reports open a way for the replacement of one parental genome by the other. In this review, we provide an overview of the current knowledge of genomic changes in interspecific homoploid and allopolyploid hybrids, with strictly homologous pairing and with relaxed pairing of homoeologs.

Ribosome Biogenesis in Plants: From Functional 45S Ribosomal DNA Organization to Ribosome Assembly Factors

The transcription of 18S, 5.8S, and 18S rRNA genes (45S rDNA), cotranscriptional processing of pre-rRNA, and assembly of mature rRNA with ribosomal proteins are the linchpins of ribosome biogenesis. In yeast (Saccharomyces cerevisiae) and animal cells, hundreds of pre-rRNA processing factors have been identified and their involvement in ribosome assembly determined. These studies, together with structural analyses, have yielded comprehensive models of the pre-40S and pre-60S ribosome subunits as well as the largest cotranscriptionally assembled preribosome particle: the 90S/small subunit processome. Here, we present the current knowledge of the functional organization of 45S rDNA, pre-rRNA transcription, rRNA processing activities, and ribosome assembly factors in plants, focusing on data from Arabidopsis (Arabidopsis thaliana). Based on yeast and mammalian cell studies, we describe the ribonucleoprotein complexes and RNA-associated activities and discuss how they might specifically affect the production of 40S and 60S subunits. Finally, we review recent findings concerning pre-rRNA processing pathways and a novel mechanism involved in a ribosome stress response in plants.

Meiotic DNA Repair in the Nucleolus Employs a Nonhomologous End-Joining Mechanism

Ribosomal RNA genes are arranged in large arrays with hundreds of rDNA units in tandem. These highly repetitive DNA elements pose a risk to genome stability since they can undergo nonallelic exchanges. During meiosis, DNA double-strand breaks (DSBs) are induced as part of the regular program to generate gametes. Meiotic DSBs initiate homologous recombination (HR), which subsequently ensures genetic exchange and chromosome disjunction. In Arabidopsis (Arabidopsis thaliana), we demonstrate that all 45S rDNA arrays become transcriptionally active and are recruited into the nucleolus early in meiosis. This shields the rDNA from acquiring canonical meiotic chromatin modifications and meiotic cohesin and allows only very limited meiosis-specific DSB formation. DNA lesions within the rDNA arrays are repaired in an RAD51-independent but LIG4-dependent manner, establishing that nonhomologous end-joining maintains rDNA integrity during meiosis. Utilizing ectopically integrated rDNA repeats, we validate our findings and demonstrate that the rDNA constitutes an HR-refractory genome environment.

Ribosomal RNA genes shape chromatin domains associating with the nucleolus

Genomic interactions can occur in addition to those within chromosome territories and can be organized around nuclear bodies. Several studies revealed how the nucleolus anchors higher order chromatin structures of specific chromosome regions displaying heterochromatic features. In this review, we comment on advances in this emerging field, with a particular focus on a recent study published by Quinodoz et al., that developed a new method to characterize simultaneous genomic interactions in the same cell. Highlighting studies conducted in animal and plant cells, we then discuss the establishment of inactive chromatin at nucleolus organizer region (NOR)-bearing chromosomes.

Ribosomal DNA loci derived from Brachypodium stacei are switched off for major parts of the life cycle of Brachypodium hybridum

Nucleolar dominance is an epigenetic phenomenon that occurs in some plant and animal allopolyploids and hybrids, whereby only one ancestral set of 35S rRNA genes retains the ability to form the nucleolus while the rDNA loci derived from the other progenitor are transcriptionally silenced. There is substantial evidence that nucleolar dominance is regulated developmentally. This study focuses upon the establishment and/or maintenance of nucleolar dominance during different stages of development in the model grass allotetraploid Brachypodium hybridum. Fluorescence in situ hybridization with a 25S rDNA probe to cells in three-dimensional cytogenetic preparations showed that nucleolar dominance is present not only in root meristematic and differentiated cells of this species, but also in male meiocytes at prophase I, tetrads of microspores, and different embryonic tissues. The inactive state of Brachypodium stacei-originated rDNA loci was confirmed by silver staining. Only B. distachyon-derived 35S rDNA loci formed nucleoli in the aforementioned tissues, whereas B. stacei-like loci remained highly condensed and thus transcriptionally suppressed. The establishment of nucleolar dominance during earlier stages of B. hybridum embryo development cannot be ruled out. However, we propose that gradual pseudogenization of B. stacei-like loci in the evolution of the allotetraploid seems to be more likely.

Two jatropha karyotypes constructed from meiotic pachytene chromosomes: Pericentric distribution of heterochromatin and variation in repetitive DNAs

Jatropha (Jatropha curcas) is an oil-bearing plant used for biodiesel production. Construction of its standard karyotype and identification of the euchromatin/heterochromatin distribution associated with gene expression and meiotic recombination are essential to fully characterize its genome. Here, we developed a J. curcas karyotype based on meiotic pachytene chromosomes. In addition, a karyotype of J. integerrima, a useful species for jatropha breeding, was also constructed. Five out of eleven J. curcas chromosomes were metacentric, but only two were metacentric in J. integerrima. Almost all of the heterochromatin was distributed around the pericentric regions. The interstitial and distal regions were euchromatic without heterochromatic knobs, except for small heterochromatin regions associated with the subtelomeric repeat sequence JcSat1. These pericentric heterochromatin distribution patterns, together with chromosome structure data and the results of FISH probing with rDNA and JcSat1, allowed us to classify all chromosomes of both species. The two species had two 35S rDNA loci and one 5S rDNA locus; one 35S rDNA locus in J. integerrima was located on the interstitial region of the short arms. In addition, JcSat1 was found at only the heterochromatic ends of the J. curcas chromosome, not the J. integerrima chromosome. Despite the same chromosome number, the two pachytene chromosome-based karyotypes suggest variation in chromosome structure and distribution of repetitive DNAs in these two species.

Structural features of two major nucleolar organizer regions (NORs), Nor-B1 and Nor-B2, and chromosome-specific rRNA gene expression in wheat

The reference genome sequence of wheat 'Chinese Spring' (CS) is now available (IWGSC RefSeq v1.0), but the core sequences defining the nucleolar organizer regions (NORs) have not been characterized. We estimated that the total copy number of the rDNA units in the wheat genome is 11 160, of which 30.5%, 60.9% and 8.6% are located on Nor-B1 (1B), Nor-B2 (6B) and other NORs, respectively. The total length of the NORs is estimated to be 100 Mb, corresponding to approximately 10% of the unassembled portion of the genome not represented in RefSeq v1.0. Four subtypes (S1-S4) of the rDNA units were identified based on differences within the 3' external transcribed spacer regions in Nor-B1 and Nor-B2, and quantitative PCR indicated locus-specific variation in rDNA subtype contents. Expression analyses of rDNA subtypes revealed that S1 was predominantly expressed and S2 weakly expressed, in contrast to the relative abundance of rDNA subtypes in the wheat genome. These results suggest a regulation mechanism of differential rDNA expression based on sequence differences. S3 expression increased in the ditelosomic lines Dt1BL and Dt6BL, suggesting that S3 is subjected to chromosome-mediated silencing. Structural differences were detected in the regions surrounding the NOR among homoeologous chromosomes of groups 1 and 6. The adjacent regions distal to the major NORs were expanded compared with their homoeologous counterparts, and the gene density of these expanded regions was relatively low. We provide evidence that these regions are likely to be important for autoregulation of the associated major NORs as well as silencing of minor NORs.

Canonical and Noncanonical Actions of Arabidopsis Histone Deacetylases in Ribosomal RNA Processing

Ribosome biogenesis is a fundamental process required for all cellular activities. Histone deacetylases play critical roles in many biological processes including transcriptional repression and rDNA silencing. However, their function in pre-rRNA processing remains poorly understood. Here, we discovered a previously uncharacterized role of Arabidopsis thaliana histone deacetylase HD2C in pre-rRNA processing via both canonical and noncanonical manners. HD2C interacts with another histone deacetylase HD2B and forms homo- and/or hetero-oligomers in the nucleolus. Depletion of HD2C and HD2B induces a ribosome-biogenesis deficient phenotype and aberrant accumulation of 18S pre-rRNA intermediates. Our genome-wide analysis revealed that HD2C binds and represses the expression of key genes involved in ribosome biogenesis. Using RNA immunoprecipitation and sequencing, we further uncovered a noncanonical mechanism of HD2C directly associating with pre-rRNA and small nucleolar RNAs to regulate rRNA methylation. Together, this study reveals a multifaceted role of HD2C in ribosome biogenesis and provides mechanistic insights into how histone deacetylases modulate rRNA maturation at the transcriptional and posttranscriptional levels.

Dynamic Changes in Plant Nuclear Organization in Response to Environmental and Developmental Signals

The functional organization of the plant nuclear pore, nuclear envelope, and nucleoplasm marks dynamically changing environmental cues and developmental programs.

Arabidopsis ATRX Modulates H3.3 Occupancy and Fine-Tunes Gene Expression

Histones are essential components of the nucleosome, the major chromatin subunit that structures linear DNA molecules and regulates access of other proteins to DNA. Specific histone chaperone complexes control the correct deposition of canonical histones and their variants to modulate nucleosome structure and stability. In this study, we characterize the Arabidopsis thaliana Alpha Thalassemia-mental Retardation X-linked (ATRX) ortholog and show that ATRX is involved in histone H3 deposition. Arabidopsis ATRX mutant alleles are viable, but show developmental defects and reduced fertility. Their combination with mutants of the histone H3.3 chaperone HIRA (Histone Regulator A) results in impaired plant survival, suggesting that HIRA and ATRX function in complementary histone deposition pathways. Indeed, ATRX loss of function alters cellular histone H3.3 pools and in consequence modulates the H3.1/H3.3 balance in the cell. H3.3 levels are affected especially at genes characterized by elevated H3.3 occupancy, including the 45S ribosomal DNA (45S rDNA) loci, where loss of ATRX results in altered expression of specific 45S rDNA sequence variants. At the genome-wide scale, our data indicate that ATRX modifies gene expression concomitantly to H3.3 deposition at a set of genes characterized both by elevated H3.3 occupancy and high expression. Together, our results show that ATRX is involved in H3.3 deposition and emphasize the role of histone chaperones in adjusting genome expression.

Epistatic and allelic interactions control expression of ribosomal RNA gene clusters in Arabidopsis thaliana

BACKGROUND

Ribosomal RNA (rRNA) accounts for the majority of the RNA in eukaryotic cells, and is encoded by hundreds to thousands of nearly identical gene copies, only a subset of which are active at any given time. In Arabidopsis thaliana, 45S rRNA genes are found in two large ribosomal DNA (rDNA) clusters and little is known about the contribution of each to the overall transcription pattern in the species.

RESULTS

By taking advantage of genome sequencing data from the 1001 Genomes Consortium, we characterize rRNA gene sequence variation within and among accessions. Notably, variation is not restricted to the pre-rRNA sequences removed during processing, but it is also present within the highly conserved ribosomal subunits. Through linkage mapping we assign these variants to a particular rDNA cluster unambiguously and use them as reporters of rDNA cluster-specific expression. We demonstrate that rDNA cluster-usage varies greatly among accessions and that rDNA cluster-specific expression and silencing is controlled via genetic interactions between entire rDNA cluster haplotypes (alleles).

CONCLUSIONS

We show that rRNA gene cluster expression is controlled via complex epistatic and allelic interactions between rDNA haplotypes that apparently regulate the entire rRNA gene cluster. Furthermore, the sequence polymorphism we discovered implies that the pool of rRNA in a cell may be heterogeneous, which could have functional consequences.