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

Transcriptional repression of GTL1 under water-deficit stress promotes anthocyanin biosynthesis to enhance drought tolerance

The transcription factor GT2-LIKE 1 (GTL1) has been implicated in orchestrating a transcriptional network of diverse physiological, biochemical, and developmental processes. In response to water-limiting conditions, GTL1 is a negative regulator of stomatal development, but its potential rolein other water-deficit responses is unknown. We hypothesized that GTL1 regulates transcriptome changes associated with drought tolerance over leaf developmental stages. To test the hypothesis, gene expression was profiled by RNA-seq analysis in emerging and expanding leaves of wild-type and a drought-tolerant gtl1-4 knockout mutant under well-watered and water-deficit conditions. Our comparative analysis of genotype-treatment combinations within leaf developmental age identified 459 and 1073 differentially expressed genes in emerging and expanding leaves, respectively, as water-deficit responsive GTL1-regulated genes. Transcriptional profiling identified a potential role of GTL1 in two important pathways previously linked to drought tolerance: flavonoid and polyamine biosynthesis. In expanding leaves, negative regulation of GTL1 under water-deficit conditions promotes biosynthesis of flavonoids and anthocyanins that may contribute to drought tolerance. Quantification of polyamines did not support a role for GTL1 in these drought-responsive pathways, but this is likely due to the complex nature of polyamine synthesis and turnover. Our global transcriptome analysis suggests that transcriptional repression of GTL1 by water deficit allows plants to activate diverse pathways that collectively contribute to drought tolerance.

Monoallelically expressed noncoding RNAs form nucleolar territories on NOR-containing chromosomes and regulate rRNA expression

Out of the several hundred copies of rRNA genes arranged in the nucleolar organizing regions (NOR) of the five human acrocentric chromosomes, ~50% remain transcriptionally inactive. NOR-associated sequences and epigenetic modifications contribute to the differential expression of rRNAs. However, the mechanism(s) controlling the dosage of active versus inactive rRNA genes within each NOR in mammals is yet to be determined. We have discovered a family of ncRNAs, SNULs (Single NUcleolus Localized RNA), which form constrained sub-nucleolar territories on individual NORs and influence rRNA expression. Individual members of the SNULs monoallelically associate with specific NOR-containing chromosomes. SNULs share sequence similarity to pre-rRNA and localize in the sub-nucleolar compartment with pre-rRNA. Finally, SNULs control rRNA expression by influencing pre-rRNA sorting to the DFC compartment and pre-rRNA processing. Our study discovered a novel class of ncRNAs influencing rRNA expression by forming constrained nucleolar territories on individual NORs.

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.

The chromatin source-sink hypothesis: a shared mode of chromatin-mediated regulations

We propose that several chromatin-mediated regulatory processes are dominated by source-sink relationships in which factors operate as 'sources' to produce or provide a resource and compete with each other to occupy separate 'sinks'. In this model, large portions of genomic DNA operate as 'sinks', which are filled by 'sources', such as available histone variants, covalent modifications to histones, the readers of these modifications and non-coding RNAs. Competing occupation for the sinks by different sources leads to distinct states of genomic equilibrium in differentiated cells. During dynamic developmental events, such as sexual reproduction, we propose that dramatic and rapid reconfiguration of source-sink relationships modifies chromatin states. We envision that re-routing of sources could occur by altering the dimensions of the sink, by reconfiguration of existing sink occupation or by varying the size of the source, providing a central mechanism to explain a plethora of epigenetic phenomena, which contribute to phenotypic variegation, zygotic genome activation and nucleolar dominance.

45S rDNA Diversity In Natura as One Step towards Ribosomal Heterogeneity in Arabidopsis thaliana

The keystone of ribosome biogenesis is the transcription of 45S rDNA. The Arabidopsis thaliana genome contains hundreds of 45S rDNA units; however, they are not all transcribed. Notably, 45S rDNA units contain insertions/deletions revealing the existence of heterogeneous rRNA genes and, likely, heterogeneous ribosomes for rRNAs. In order to obtain an overall picture of 45S rDNA diversity sustaining the synthesis of rRNAs and, subsequently, of ribosomes in natura, we took advantage of 320 new occurrences of Arabidopsis thaliana as a metapopulation named At66, sampled from 0 to 1900 m of altitude in the eastern Pyrenees in France. We found that the 45S rDNA copy number is very dynamic in natura and identified new genotypes for both 5' and 3' External Transcribed Spacers (ETS). Interestingly, the highest 5'ETS genotype diversity is found in altitude while the highest 3'ETS genotype diversity is found at sea level. Structural analysis of 45S rDNA also shows conservation in natura of specific 5'ETS and 3'ETS sequences/features required to control rDNA expression and the processing of rRNAs. In conclusion, At66 is a worthwhile natural laboratory, and unraveled 45S rDNA diversity represents an interesting starting material to select subsets for rDNA transcription and alter the rRNA composition of ribosomes both intra- and inter-site.

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.

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.

Substantial rDNA copy number reductions alter timing of development and produce variable tissue-specific phenotypes in C. elegans

The genes that encode ribosomal RNAs are present in several hundred copies in most eukaryotes. These vast arrays of repetitive ribosomal DNA (rDNA) have been implicated not just in ribosome biogenesis, but also aging, cancer, genome stability, and global gene expression. rDNA copy number is highly variable among and within species; this variability is thought to associate with traits relevant to human health and disease. Here we investigate the phenotypic consequences of multicellular life at the lower bounds of rDNA copy number. We use the model Caenorhabditis elegans, which has previously been found to complete embryogenesis using only maternally provided ribosomes. We find that individuals with rDNA copy number reduced to ∼5% of wild type are capable of further development with variable penetrance. Such individuals are sterile and exhibit severe morphological defects, particularly in post-embryonically dividing tissues such as germline and vulva. Developmental completion and fertility are supported by an rDNA copy number ∼10% of wild type, with substantially delayed development. Worms with rDNA copy number reduced to ∼33% of wild type display a subtle developmental timing defect that was absent in worms with higher copy numbers. Our results support the hypothesis that rDNA requirements vary across tissues and indicate that the minimum rDNA copy number for fertile adulthood is substantially less than the lowest naturally observed total copy number. The phenotype of individuals with severely reduced rDNA copy number is highly variable in penetrance and presentation, highlighting the need for continued investigation into the biological consequences of rDNA copy number variation.

Pushing the limits of HiFi assemblies reveals centromere diversity between two Arabidopsis thaliana genomes

Although long-read sequencing can often enable chromosome-level reconstruction of genomes, it is still unclear how one can routinely obtain gapless assemblies. In the model plant Arabidopsis thaliana, other than the reference accession Col-0, all other accessions de novo assembled with long-reads until now have used PacBio continuous long reads (CLR). Although these assemblies sometimes achieved chromosome-arm level contigs, they inevitably broke near the centromeres, excluding megabases of DNA from analysis in pan-genome projects. Since PacBio high-fidelity (HiFi) reads circumvent the high error rate of CLR technologies, albeit at the expense of read length, we compared a CLR assembly of accession Eyach15-2 to HiFi assemblies of the same sample. The use of five different assemblers starting from subsampled data allowed us to evaluate the impact of coverage and read length. We found that centromeres and rDNA clusters are responsible for 71% of contig breaks in the CLR scaffolds, while relatively short stretches of GA/TC repeats are at the core of >85% of the unfilled gaps in our best HiFi assemblies. Since the HiFi technology consistently enabled us to reconstruct gapless centromeres and 5S rDNA clusters, we demonstrate the value of the approach by comparing these previously inaccessible regions of the genome between the Eyach15-2 accession and the reference accession Col-0.

Impaired function of rDNA transcription initiation machinery leads to derepression of ribosomal genes with insertions of R2 retrotransposon

Eukaryotic genomes harbor hundreds of rRNA genes, many of which are transcriptionally silent. However, little is known about selective regulation of individual rDNA units. In Drosophila melanogaster, some rDNA repeats contain insertions of the R2 retrotransposon, which is capable to be transcribed only as part of pre-rRNA molecules. rDNA units with R2 insertions are usually inactivated, although R2 expression may be beneficial in cells with decreased rDNA copy number. Here we found that R2-inserted rDNA units are enriched with HP1a and H3K9me3 repressive mark, whereas disruption of the heterochromatin components slightly affects their silencing in ovarian germ cells. Surprisingly, we observed a dramatic upregulation of R2-inserted rRNA genes in ovaries lacking Udd (Under-developed) or other subunits (TAF1b and TAF1c-like) of the SL1-like complex, which is homologues to mammalian Selective factor 1 (SL1) involved in rDNA transcription initiation. Derepression of rRNA genes with R2 insertions was accompanied by a reduction of H3K9me3 and HP1a enrichment. We suggest that the impairment of the SL1-like complex affects a mechanism of selective activation of intact rDNA units which competes with heterochromatin formation. We also propose that R2 derepression may serve as an adaptive response to compromised rRNA synthesis.

Characterization of the Different Levels of Variation in 45S rRNA Genes

Although commonly regarded as nearly identical in sequence, 45S ribosomal RNA genes represent a massive source of genetic variation at different levels. Unfortunately, due to their repetitive nature and the difficulty to assemble their multiple copies in tandem, these important genomic elements remain largely unexplored in sequencing projects. Here, I describe how to exploit next generation sequencing data to estimate their copy number in an organism and detect true polymorphic sites within and among individuals. Furthermore, for species that carry multiple 45S ribosomal RNA gene clusters, I show how to make use of experimental populations to assign some of these variants to their cluster of origin.

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

The epigenetic origin of life history transitions in plants and algae

Plants and algae have a complex life history that transitions between distinct life forms called the sporophyte and the gametophyte. This phenomenon-called the alternation of generations-has fascinated botanists and phycologists for over 170 years. Despite the mesmerizing array of life histories described in plants and algae, we are only now beginning to learn about the molecular mechanisms controlling them and how they evolved. Epigenetic silencing plays an essential role in regulating gene expression during multicellular development in eukaryotes, raising questions about its impact on the life history strategy of plants and algae. Here, we trace the origin and function of epigenetic mechanisms across the plant kingdom, from unicellular green algae through to angiosperms, and attempt to reconstruct the evolutionary steps that influenced life history transitions during plant evolution. Central to this evolutionary scenario is the adaption of epigenetic silencing from a mechanism of genome defense to the repression and control of alternating generations. We extend our discussion beyond the green lineage and highlight the peculiar case of the brown algae. Unlike their unicellular diatom relatives, brown algae lack epigenetic silencing pathways common to animals and plants yet display complex life histories, hinting at the emergence of novel life history controls during stramenopile evolution.

Licorice Germplasm Resources Identification Using DNA Barcodes Inner-Variants

Based on the gradual transformation from wild growth to artificial cultivation, the accurate authentication of licorice seeds contributes to the first committed step of its quality control and is pivotal to ensure the clinical efficacy of licorice. However, it is still challenging to obtain genetically stable licorice germplasm resources due to the multi-source, multi-heterozygous, polyploid, and hybrid characteristics of licorice seeds. Here, a new method for determining the heterozygosity of licorice seed mixture, based on the various sites, and finding the composition characteristics of licorice seed is preliminarily designed and proposed. Namely, high-throughput full-length multiple DNA barcodes(HFMD), based on ITS multi-copy variation exist, the full-length amplicons of ITS2, psbA-trnH and ITS are directly sequenced by rDNA through the next-generation sequence(NGS) and single-molecule real-time (SMRT) technologies. By comparing the three sequencing methods, our results proved that SMRT sequencing successfully identified the complete gradients of complex mixed samples with the best performance. Meanwhile, HFMD is a brilliant and feasible method for evaluating the heterozygosity of licorice seeds. It shows a perfect interpretation of DNA barcoding and can be applied in multi-base multi-heterozygous and polyploid species.

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.

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.

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.

It Is Just a Matter of Time: Balancing Homologous Recombination and Non-homologous End Joining at the rDNA Locus During Meiosis

Ribosomal RNA genes (rDNAs) are located in large domains of hundreds of rDNA units organized in a head-to-tail manner. The proper and stable inheritance of rDNA clusters is of paramount importance for survival. Yet, these highly repetitive elements pose a potential risk to the genome since they can undergo non-allelic exchanges. Here, we review the current knowledge of the organization of the rDNA clusters in Arabidopsis thaliana and their stability during meiosis. Recent findings suggest that during meiosis, all rDNA loci are embedded within the nucleolus favoring non-homologous end joining (NHEJ) as a repair mechanism, while DNA repair via homologous recombination (HR) appears to be a rare event. We propose a model where (1) frequent meiotic NHEJ events generate abundant single nucleotide polymorphisms and insertions/deletions within the rDNA, resulting in a heterogeneous population of rDNA units and (2) rare HR events dynamically change rDNA unit numbers, only to be observed in large populations over many generations. Based on the latest efforts to delineate the entire rDNA sequence in A. thaliana, we discuss evidence supporting this model. The results compiled so far draw a surprising picture of rDNA sequence heterogeneity between individual units. Furthermore, rDNA cluster sizes have been recognized as relatively stable when observing less than 10 generations, yet emerged as major determinant of genome size variation between different A. thaliana ecotypes. The sequencing efforts also revealed that transcripts from the diverse rDNA units yield heterogenous ribosome populations with potential functional implications. These findings strongly motivate further research to understand the mechanisms that maintain the metastable state of rDNA loci.

Gradual evolution of allopolyploidy in Arabidopsis suecica

Most diploid organisms have polyploid ancestors. The evolutionary process of polyploidization is poorly understood but has frequently been conjectured to involve some form of 'genome shock', such as genome reorganization and subgenome expression dominance. Here we study polyploidization in Arabidopsis suecica, a post-glacial allopolyploid species formed via hybridization of Arabidopsis thaliana and Arabidopsis arenosa. We generated a chromosome-level genome assembly of A. suecica and complemented it with polymorphism and transcriptome data from all species. Despite a divergence around 6 million years ago (Ma) between the ancestral species and differences in their genome composition, we see no evidence of a genome shock: the A. suecica genome is colinear with the ancestral genomes; there is no subgenome dominance in expression; and transposon dynamics appear stable. However, we find changes suggesting gradual adaptation to polyploidy. In particular, the A. thaliana subgenome shows upregulation of meiosis-related genes, possibly to prevent aneuploidy and undesirable homeologous exchanges that are observed in synthetic A. suecica, and the A. arenosa subgenome shows upregulation of cyto-nuclear processes, possibly in response to the new cytoplasmic environment of A. suecica, with plastids maternally inherited from A. thaliana. These changes are not seen in synthetic hybrids, and thus are likely to represent subsequent evolution.

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.

ASY1 acts as a dosage-dependent antagonist of telomere-led recombination and mediates crossover interference in Arabidopsis

During meiosis, interhomolog recombination produces crossovers and noncrossovers to create genetic diversity. Meiotic recombination frequency varies at multiple scales, with high subtelomeric recombination and suppressed centromeric recombination typical in many eukaryotes. During recombination, sister chromatids are tethered as loops to a polymerized chromosome axis, which, in plants, includes the ASY1 HORMA domain protein and REC8-cohesin complexes. Using chromatin immunoprecipitation, we show an ascending telomere-to-centromere gradient of ASY1 enrichment, which correlates strongly with REC8-cohesin ChIP-seq data. We mapped crossovers genome-wide in the absence of ASY1 and observe that telomere-led recombination becomes dominant. Surprisingly, asy1/+ heterozygotes also remodel crossovers toward subtelomeric regions at the expense of the pericentromeres. Telomeric recombination increases in asy1/+ occur in distal regions where ASY1 and REC8 ChIP enrichment are lowest in wild type. In wild type, the majority of crossovers show interference, meaning that they are more widely spaced along the chromosomes than expected by chance. To measure interference, we analyzed double crossover distances, MLH1 foci, and fluorescent pollen tetrads. Interestingly, while crossover interference is normal in asy1/+, it is undetectable in asy1 mutants, indicating that ASY1 is required to mediate crossover interference. Together, this is consistent with ASY1 antagonizing telomere-led recombination and promoting spaced crossover formation along the chromosomes via interference. These findings provide insight into the role of the meiotic axis in patterning recombination frequency within plant genomes.

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.

Chromatin Organization in Early Land Plants Reveals an Ancestral Association between H3K27me3, Transposons, and Constitutive Heterochromatin

Genome packaging by nucleosomes is a hallmark of eukaryotes. Histones and the pathways that deposit, remove, and read histone modifications are deeply conserved. Yet, we lack information regarding chromatin landscapes in extant representatives of ancestors of the main groups of eukaryotes, and our knowledge of the evolution of chromatin-related processes is limited. We used the bryophyte Marchantia polymorpha, which diverged from vascular plants circa 400 mya, to obtain a whole chromosome genome assembly and explore the chromatin landscape and three-dimensional genome organization in an early diverging land plant lineage. Based on genomic profiles of ten chromatin marks, we conclude that the relationship between active marks and gene expression is conserved across land plants. In contrast, we observed distinctive features of transposons and other repetitive sequences in Marchantia compared with flowering plants. Silenced transposons and repeats did not accumulate around centromeres. Although a large fraction of constitutive heterochromatin was marked by H3K9 methylation as in flowering plants, a significant proportion of transposons were marked by H3K27me3, which is otherwise dedicated to the transcriptional repression of protein-coding genes in flowering plants. Chromatin compartmentalization analyses of Hi-C data revealed that repressed B compartments were densely decorated with H3K27me3 but not H3K9 or DNA methylation as reported in flowering plants. We conclude that, in early plants, H3K27me3 played an essential role in heterochromatin function, suggesting an ancestral role of this mark in transposon silencing.

Challenges and Approaches to Genotyping Repetitive DNA

Individuals within a species can exhibit vast variation in copy number of repetitive DNA elements. This variation may contribute to complex traits such as lifespan and disease, yet it is only infrequently considered in genotype-phenotype associations. Although the possible importance of copy number variation is widely recognized, accurate copy number quantification remains challenging. Here, we assess the technical reproducibility of several major methods for copy number estimation as they apply to the large repetitive ribosomal DNA array (rDNA). rDNA encodes the ribosomal RNAs and exists as a tandem gene array in all eukaryotes. Repeat units of rDNA are kilobases in size, often with several hundred units comprising the array, making rDNA particularly intractable to common quantification techniques. We evaluate pulsed-field gel electrophoresis, droplet digital PCR, and Nextera-based whole genome sequencing as approaches to copy number estimation, comparing techniques across model organisms and spanning wide ranges of copy numbers. Nextera-based whole genome sequencing, though commonly used in recent literature, produced high error. We explore possible causes for this error and provide recommendations for best practices in rDNA copy number estimation. We present a resource of high-confidence rDNA copy number estimates for a set of S. cerevisiae and C. elegans strains for future use. We furthermore explore the possibility for FISH-based copy number estimation, an alternative that could potentially characterize copy number on a cellular level.

The effects of Arabidopsis genome duplication on the chromatin organization and transcriptional regulation

Autopolyploidy is widespread in higher plants and important for agricultural yield and quality. However, the effects of genome duplication on the chromatin organization and transcriptional regulation are largely unknown in plants. Using High-throughput Chromosome Conformation Capture (Hi-C), we showed that autotetraploid Arabidopsis presented more inter-chromosomal interactions and fewer short-range chromatin interactions compared with its diploid progenitor. In addition, genome duplication contributed to the switching of some loose and compact structure domains with altered H3K4me3 and H3K27me3 histone modification status. 539 genes were identified with altered transcriptions and chromatin interactions in autotetraploid Arabidopsis. Especially, we found that genome duplication changed chromatin looping and H3K27me3 histone modification in Flowering Locus C. We propose that genome doubling modulates the transcription genome-wide by changed chromatin interactions and at the specific locus by altered chromatin loops and histone modifications.

Genetic and epigenetic variation in 5S ribosomal RNA genes reveals genome dynamics in Arabidopsis thaliana

Organized in tandem repeat arrays in most eukaryotes and transcribed by RNA polymerase III, expression of 5S rRNA genes is under epigenetic control. To unveil mechanisms of transcriptional regulation, we obtained here in depth sequence information on 5S rRNA genes from the Arabidopsis thaliana genome and identified differential enrichment in epigenetic marks between the three 5S rDNA loci situated on chromosomes 3, 4 and 5. We reveal the chromosome 5 locus as the major source of an atypical, long 5S rRNA transcript characteristic of an open chromatin structure. 5S rRNA genes from this locus translocated in the Landsberg erecta ecotype as shown by linkage mapping and chromosome-specific FISH analysis. These variations in 5S rDNA locus organization cause changes in the spatial arrangement of chromosomes in the nucleus. Furthermore, 5S rRNA gene arrangements are highly dynamic with alterations in chromosomal positions through translocations in certain mutants of the RNA-directed DNA methylation pathway and important copy number variations among ecotypes. Finally, variations in 5S rRNA gene sequence, chromatin organization and transcripts indicate differential usage of 5S rDNA loci in distinct ecotypes. We suggest that both the usage of existing and new 5S rDNA loci resulting from translocations may impact neighboring chromatin organization.

The fate of ribosomal RNA genes in spontaneous polyploid dogrose hybrids [Rosa L. sect. Caninae (DC.) Ser.] exhibiting non-symmetrical meiosis

Dogroses represent an exceptional system for studying the effects of genome doubling and hybridization: their asymmetrical meiosis enables recombination in bi-parentally inherited chromosomes but prevents it in maternally inherited ones. We employed fluorescent in situ hybridization, genome skimming, amplicon sequencing of genomic and cDNA as well as conventional cloning of nuclear ribosomal DNA in two phylogenetically distinct pentaploid (2n = 5x = 35) species, Rosa canina and Rosa inodora, and their naturally occurring reciprocal hybrids, Rosa dumalis (5x) and Rosa agrestis (5x, 6x). Both progenitor species differed in composition, meiotic behaviour and expression of rDNA loci: R. canina (five 18S and 5-8 5S loci) was dominated by the Canina ribotypes, but R. inodora (four 18S loci and 7-8 5S loci) by the Rubiginosa ribotype. The co-localized 5S/18S loci occurred on either bivalent-forming (R. canina) or univalent-forming (R. inodora) chromosomes. Ribosomal DNA loci were additively inherited; however, the Canina ribotypes were dominantly expressed, even in genotypes with relatively low copy number of these genes. Moreover, we observed rDNA homogenization towards the paternally transmitted Canina ribotype in 6x R. agrestis. The here-observed variation in arrangement and composition of rDNA types between R. canina and R. inodora suggests the involvement of different genomes in bivalent formation. This results supports the hypothesis that the asymmetrical meiosis arose at least twice by independent ancient hybridization events.