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

Uniparental silencing of 5S rRNA genes in plant allopolyploids - insights from Cardamine (Brassicaceae)

While the phenomenon of uniparental silencing of 35S rDNA in interspecific hybrids and allopolyploids is well documented, there is a notable absence of information regarding whether such silencing extends to the 5S RNA component of ribosomes. To address this gap in knowledge, we analyzed the 5S and 35S rDNA expression in Cardamine (Brassicaceae) allopolyploids, namely C. × insueta (2n = 3x = 24, genome composition RRA), C. flexuosa (2n = 4x = 32, AAHH), and C. scutata (2n = 4x = 32, PPAA) which share a common diploid ancestor (AA). We employed high-throughput sequencing of transcriptomes and genomes and phylogenetic analyses of 5S rRNA variants. The genomic organization of rDNA was further scrutinized through clustering and fluorescence in situ hybridization. In the C. × insueta allotriploid, we observed uniparental dominant expression of 5S and 35S rDNA loci. In the C. flexuosa and C. scutata allotetraploids, the expression pattern differed, with the 35S rDNA being expressed from the A subgenome, whereas the 5S rDNA was expressed from the partner subgenome. Both C. flexuosa and C. scutata but not C. × insueta showed copy and locus number changes. We conclude that in stabilized allopolyploids, transcription of ribosomal RNA components occurs from different subgenomes. This phenomenon appears to result in the formation of chimeric ribosomes comprising rRNA molecules derived from distinct parental origins. We speculate that the interplay of epigenetic silencing and rDNA rearrangements introduces an additional layer of variation in multimolecule ribosomal complexes, potentially contributing to the evolutionary success of allopolyploids.

Contrasting patterns of 5S rDNA repeats in European and Asian ecotypes of greater duckweed, Spirodela polyrhiza (Lemnaceae)

Ribosomal DNA (rDNA) contains highly conserved, specifically organized sequences encoding ribosomal RNAs (rRNAs) separated by variable non-transcribed intergenic spacers (NTSs) and is abundant in eukaryotic genomes. These characteristics make the rDNA an informative molecular target to study genome organization, molecular evolution, and phylogenetics. In this study, we characterized the 5S rDNA repeats in the greater duckweed Spiroldela polyrhiza, a species known for its small size, rapid growth, highly conserved genome organization, and low mutation rate. Sequence analysis of at least 12 individually cloned PCR fragments containing the 5S rDNA units for each of six ecotypes that originated from Europe (Ukraine) and Asia (China) revealed two distinct types of 5S rDNA repeats containing NTSs of different lengths and nucleotide compositions. The shorter 5S rDNA repeat units had a highly homogeneous 400-bp NTS, with few ecotype- or region-specific single-nucleotide polymorphisms (SNPs). The longer 5S rDNA units had NTSs of 1056-1084 bp with characteristic intra- and inter-genomic variants due to specific SNPs and insertions/deletions of 4-15-bp DNA elements. We also detected significant variability in the ratio of short/long 5S rDNA variants between ecotypes of S. polyrhiza. The contrasting dynamics of the two types of 5S rDNA units, combined with the unusually low repeat copy number (for plants) in S. polyrhiza (46-220 copies per genome), shows that this species could serve as an excellent model for examining the mechanisms of concerted evolution and functional significance of rDNA variability.

Telomere-to-telomere genome assembly of melon (Cucumis melo L. var. inodorus) provides a high-quality reference for meta-QTL analysis of important traits

Melon is an important horticultural crop with extensive diversity in many horticultural groups. To explore its genomic diversity, it is necessary to assemble more high-quality complete genomes from different melon accessions. Meanwhile, a large number of QTLs have been mapped in several studies. Integration of the published QTLs onto a complete genome can provide more accurate information for candidate gene cloning. To address these problems, a telomere-to-telomere (T2T) genome of the elite melon landrace Kuizilikjiz (Cucumis melo L. var. inodorus) was de novo assembled and all the published QTLs were projected onto it in this study. The results showed that a high-quality Kuizilikjiz genome with the size of 379.2 Mb and N50 of 31.7 Mb was de novo assembled using the combination of short reads, PacBio high-fidelity long reads, Hi-C data, and a high-density genetic map. Each chromosome contained the centromere and telomeres at both ends. A large number of structural variations were observed between Kuizilikjiz and the other published genomes. A total of 1294 QTLs published in 67 studies were collected and projected onto the T2T genome. Several clustered, co-localized, and overlapped QTLs were determined. Furthermore, 20 stable meta-QTLs were identified, which significantly reduced the mapping intervals of the initial QTLs and greatly facilitated identification of the candidate genes. Collectively, the T2T genome assembly together with the numerous projected QTLs will not only broaden the high-quality genome resources but also provide valuable and abundant QTL information for cloning the genes controlling important traits in melon.

Discrimination of Green Coffee (Coffea arabica and Coffea canephora) of Different Geographical Origin Based on Antioxidant Activity, High-Throughput Metabolomics, and DNA RFLP Fingerprinting

The genus Coffea is known for the two species C. arabica (CA) and C. canephora (CC), which are used to prepare the beverage coffee. Proper identification of green beans of coffee varieties is based on phenotypic and phytochemical/molecular characteristics. In this work, a combination of chemical (UV/Vis, HPLC-DAD-MS/MS, GC-MS, and GC-FID) and molecular (PCR-RFLP) fingerprinting was used to discriminate commercial green coffee accessions from different geographical origin. The highest content of polyphenols and flavonoids was always found in CC accessions, whereas CA showed lower values. ABTS and FRAP assays showed a significant correlation between phenolic content and antioxidant activity in most CC accessions. We identified 32 different compounds, including 28 flavonoids and four N-containing compounds. The highest contents of caffeine and melatonin were detected in CC accessions, whereas the highest levels of quercetin and kaempferol derivatives were found in CA accessions. Fatty acids of CC accessions were characterized by low levels of linoleic and cis octadecenoic acid and high amounts of elaidic acid and myristic acid. Discrimination of species according to their geographical origin was achieved using high-throughput data analysis, combining all measured parameters. Lastly, PCR-RFLP analysis was instrumental for the identification of recognition markers for the majority of accessions. Using the restriction enzyme AluI on the trnL-trnF region, we clearly discriminated C. canephora from C. arabica, whereas the cleavage performed by the restriction enzymes MseI and XholI on the 5S-rRNA-NTS region produced specific discrimination patterns useful for the correct identification of the different coffee accessions. This work extends our previous studies and provides new information on the complete flavonoid profile, combining high-throughput data with DNA fingerprinting to assess the geographical discrimination of green coffee.

Persistence of parental age effect on somatic mutation rates across generations in Arabidopsis

In the model plant Arabidopsis thaliana, parental age is known to affect somatic mutation rates in their immediate progeny and here we show that this age dependent effect persists across successive generations. Using a set of detector lines carrying the mutated uidA gene, we examined if a particular parental age maintained across five consecutive generations affected the rates of base substitution (BSR), intrachromosomal recombination (ICR), frameshift mutation (FS), and transposition. The frequency of functional GUS reversions were assessed in seedlings as a function of identical/different parental ages across generations. In the context of a fixed parental age, BSR/ICR rates were unaffected in the first three generations, then dropped significantly in the 4th and increased in most instances in the 5th generation (e.g. BSR (F1 38 = 0.9, F2 38 = 1.14, F3 38 = 1.02, F4 38 = 0.5, F5 38 = 0.76)). On the other hand, with advancing parental ages, BSR/ICR rates remained high in the first two/three generations, with a striking resemblance in the pattern of mutation rates (BSR (F1 38 = 0.9, F1 43 = 0.53, F1 48 = 0.79, F1 53 = 0.83 and F2 38 = 1.14, F2 43 = 0.57, F2 48 = 0.64, F2 53 = 0.94). We adopted a novel approach of identifying and tagging flowers pollinated on a particular day, thereby avoiding biases due to potential emasculation induced stress responses. Our results suggest a time component in counting the number of generations a plant has passed through self-fertilization at a particular age in determining the somatic mutation rates.

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.

Agro-Morphological and Cytogenetic Characterization of Colchicine-Induced Tetraploid Plants of Polemonium caeruleum L. (Polemoniaceae)

Polemonium caeruleum L. (Polemoniaceae) is a valuable medicinal herb with a wide spectrum of biological activities. Under natural conditions, the productivity of this species is rather low. In this study, colchicine-induced tetraploid plants (2n = 4x = 36) of P. caeruleum were obtained, and for the first time, their morphological and cytogenetic characterization was performed. In the tetraploid plants, raw material productivity and also the content of triterpene saponins were significantly higher than in the control diploids. The analysis of chromosome behavior at meiosis and FISH chromosome mapping of 45S and 5S rDNA generally demonstrated stability of both genomes in the tetraploid plants. Based on chromosome morphology and distribution patterns of the studied molecular cytogenetic markers, all chromosome pairs in karyotypes were identified, and chromosome karyograms and idiograms of P. caeruleum were constructed. The revealed specific microdiagnostic characteristics of P. caeruleum (strongly sinuous cells and anomocytic stomata of the leaf epidermis, and also glandular hairs along the veins) could be useful for raw material identification. In the obtained tetraploids, the predominance of large stomata on the lower leaf epidermis was determined. The studied tetraploids can be used in various breeding programs to obtain high-quality pharmaceutical raw materials of P. caeruleum.

Organization and epigenomic control of RNA polymerase III-transcribed genes in plants

The genetic information linearly scripted in chromosomes is wrapped in a ribonucleoprotein complex called chromatin. The adaptation of its compaction level and spatiotemporal organization refines gene expression in response to developmental and environmental cues. RNA polymerase III (RNAPIII) is responsible for the biogenesis of elementary non-coding RNAs. Their genes are subjected to high duplication and mutational rates, and invade nuclear genomes. Their insertion into different epigenomic environments raises the question of how chromatin packing affects their individual transcription. In this review, we provide a unique perspective to this issue in plants. In addition, we discuss how the genomic organization of RNAPIII-transcribed loci, combined with epigenetic differences, might participate to plant trait variations.

Optimization of Genome Knock-In Method: Search for the Most Efficient Genome Regions for Transgene Expression in Plants

Plant expression systems are currently regarded as promising alternative platforms for the production of recombinant proteins, including the proteins for biopharmaceutical purposes. However, the accumulation level of a target protein in plant expression systems is still rather low compared with the other existing systems, namely, mammalian, yeast, and E. coli cells. To solve this problem, numerous methods and approaches have been designed and developed. At the same time, the random nature of the distribution of transgenes over the genome can lead to gene silencing, variability in the accumulation of recombinant protein, and also to various insertional mutations. The current research study considered inserting target genes into pre-selected regions of the plant genome (genomic “safe harbors”) using the CRISPR/Cas system. Regions of genes expressed constitutively and at a high transcriptional level in plant cells (housekeeping genes) that are of interest as attractive targets for the delivery of target genes were characterized. The results of the first attempts to deliver target genes to the regions of housekeeping genes are discussed. The approach of “euchromatization” of the transgene integration region using the modified dCas9 associated with transcription factors is considered. A number of the specific features in the spatial chromatin organization allowing individual genes to efficiently transcribe are discussed.

5S Ribosomal DNA of Genus Solanum: Molecular Organization, Evolution, and Taxonomy

The Solanum genus, being one of the largest among high plants, is distributed worldwide and comprises about 1,200 species. The genus includes numerous agronomically important species such as Solanum tuberosum (potato), Solanum lycopersicum (tomato), and Solanum melongena (eggplant) as well as medical and ornamental plants. The huge Solanum genus is a convenient model for research in the field of molecular evolution and structural and functional genomics. Clear knowledge of evolutionary relationships in the Solanum genus is required to increase the effectiveness of breeding programs, but the phylogeny of the genus is still not fully understood. The rapidly evolving intergenic spacer region (IGS) of 5S rDNA has been successfully used for inferring interspecific relationships in several groups of angiosperms. Here, combining cloning and sequencing with bioinformatic analysis of genomic data available in the SRA database, we evaluate the molecular organization and diversity of IGS for 184 accessions, representing 137 species of the Solanum genus. It was found that the main mechanisms of IGS molecular evolution was step-wise accumulation of single base substitution or short indels, and that long indels and multiple base substitutions, which arose repeatedly during evolution, were mostly not conserved and eliminated. The reason for this negative selection seems to be association between indels/multiple base substitutions and pseudogenization of 5S rDNA. Comparison of IGS sequences allowed us to reconstruct the phylogeny of the Solanum genus. The obtained dendrograms are mainly congruent with published data: same major and minor clades were found. However, relationships between these clades and position of some species (S. cochoae, S. clivorum, S. macrocarpon, and S. spirale) were different from those of previous results and require further clarification. Our results show that 5S IGS represents a convenient molecular marker for phylogenetic studies on the Solanum genus. In particular, the simultaneous presence of several structural variants of rDNA in the genome enables the detection of reticular evolution, especially in the largest and economically most important sect. Petota. The origin of several polyploid species should be reconsidered.

A haplotype resolved chromosomal level avocado genome allows analysis of novel avocado genes

Avocado (Persea americana) is a member of the magnoliids, an early branching lineage of angiosperms that has high value globally with the fruit being highly nutritious. Here, we report a chromosome-level genome assembly for the commercial avocado cultivar Hass, which represents 80% of the world's avocado consumption. The DNA contigs produced from Pacific Biosciences HiFi reads were further assembled using a previously published version of the genome supported by a genetic map. The total assembly was 913 Mb with a contig N50 of 84 Mb. Contigs assigned to the 12 chromosomes represented 874 Mb and covered 98.8% of benchmarked single-copy genes from embryophytes. Annotation of protein coding sequences identified 48 915 avocado genes of which 39 207 could be ascribed functions. The genome contained 62.6% repeat elements. Specific biosynthetic pathways of interest in the genome were investigated. The analysis suggested that the predominant pathway of heptose biosynthesis in avocado may be through sedoheptulose 1,7 bisphosphate rather than via alternative routes. Endoglucanase genes were high in number, consistent with avocado using cellulase for fruit ripening. The avocado genome appeared to have a limited number of translocations between homeologous chromosomes, despite having undergone multiple genome duplication events. Proteome clustering with related species permitted identification of genes unique to avocado and other members of the Lauraceae family, as well as genes unique to species diverged near or prior to the divergence of monocots and eudicots. This genome provides a tool to support future advances in the development of elite avocado varieties with higher yields and fruit quality.

Genomic architecture of 5S rDNA cluster and its variations within and between species

BACKGROUND

Ribosomal DNAs (rDNAs) are arranged in purely tandem repeats, preventing them from being reliably assembled onto chromosomes during generation of genome assembly. The uncertainty of rDNA genomic structure presents a significant barrier for studying their function and evolution.

RESULTS

Here we generate ultra-long Oxford Nanopore Technologies (ONT) and short NGS reads to delineate the architecture and variation of the 5S rDNA cluster in the different strains of C. elegans and C. briggsae. We classify the individual rDNA's repeating units into 25 types based on the unique sequence variations in each unit of C. elegans (N2). We next perform assembly of the cluster by taking advantage of the long reads that carry these units, which led to an assembly of 5S rDNA cluster consisting of up to 167 consecutive 5S rDNA units in the N2 strain. The ordering and copy number of various rDNA units are consistent with the separation time between strains. Surprisingly, we observed a drastically reduced level of variation in the unit composition in the 5S rDNA cluster in the C. elegans CB4856 and C. briggsae AF16 strains than in the C. elegans N2 strain, suggesting that N2, a widely used reference strain, is likely to be defective in maintaining the 5S rDNA cluster stability compared with other wild isolates of C. elegans or C. briggsae.

CONCLUSIONS

The results demonstrate that Nanopore DNA sequencing reads are capable of generating assembly of highly repetitive sequences, and rDNA units are highly dynamic both within and between population(s) of the same species in terms of sequence and copy number. The detailed structure and variation of the 5S rDNA units within the rDNA cluster pave the way for functional and evolutionary studies.

Molecular and Cytogenetic Analysis of rDNA Evolution in Crepis Sensu Lato

Although Crepis was the first model plant group in which chromosomal changes were considered to play an important role in speciation, their chromosome structure and evolution have been barely investigated using molecular cytogenetic methods. The aim of the study was to provide a better understanding of the patterns and directions of Crepis chromosome evolution, using comparative analyses of rDNA loci number and localisation. The chromosome base number and chromosomal organisation of 5S and 35S rDNA loci were analysed in the phylogenetic background for 39 species of Crepis, which represent the evolutionary lineages of Crepis sensu stricto and Lagoseris, including Lapsana communis. The phylogenetic relationships among all the species were inferred from nrITS and newly obtained 5S rDNA NTS sequences. Despite high variations in rDNA loci chromosomal organisation, most species had a chromosome with both rDNA loci within the same (usually short) chromosomal arm. The comparative analyses revealed several independent rDNA loci number gains and loci repositioning that accompanied diversification and speciation in Crepis. Some of the changes in rDNA loci patterns were reconstructed for the same evolutionary lineages as descending dysploidy.

The Ribosomal DNA Loci of the Ancient Monocot Pistia stratiotes L. (Araceae) Contain Different Variants of the 35S and 5S Ribosomal RNA Gene Units

The freshwater plant water lettuce (Pistia stratiotes L.) grows in warm climatic zones and is used for phytoremediation and biomass production. P. stratiotes belongs to the Araceae, an ecologically and structurally diverse early monocot family, but the phylogenetic relationships among Araceae members are poorly understood. Ribosomal DNAs (rDNAs), including the 35S and 5S rDNA, encode the RNA components of ribosomes and are widely used in phylogenetic and evolutionary studies of various plant taxa. Here, we comprehensively characterized the chromosomal locations and molecular organization of 35S and 5S rDNA genes in water lettuce using karyological and molecular methods. Fluorescence in situ hybridization revealed a single location for the 35S and 5S rDNA loci, each on a different pair of the species' 28 chromosomes. Molecular cloning and nucleotide sequencing of 35S rDNA of P. stratiotes, the first representative Araceae sensu stricto in which such a study was performed, displayed typical structural characteristics. The full-length repeat showed high sequence conservation of the regions producing the 18S, 5.8S, and 25S rRNAs and divergence of the internal transcribed spacers ITS1 and ITS2 as well as the large intergenic spacer (IGS). Alignments of the deduced sequence of 18S rDNA with the sequences available for other Araceae and representatives of other clades were used for phylogenetic analysis. Examination of 11 IGS sequences revealed significant intra-genomic length variability due to variation in subrepeat number, with four types of units detected within the 35S rDNA locus of the P. stratiotes genome (estimated size 407 Mb/1C). Similarly, the 5S rDNA locus harbors gene units comprising a conserved 119-bp sequence encoding 5S rRNA and two types of non-transcribed spacer (NTS) sequences. Type I was classified into four subtypes, which apparently originated via progressive loss of subrepeats within the duplicated NTS region containing the 3' part of the 5S rRNA gene. The minor Type II NTS is shorter than Type I and differs in nucleotide composition. Some DNA clones containing two or three consecutive 5S rDNA repeats harbored 5S rDNA genes with different types of NTSs, confirming the mosaic composition of the 5S rDNA locus.

Sulfuriroseicoccus oceanibius gen. nov., sp. nov., a representative of the phylum Verrucomicrobia with a special cytoplasmic membrane

Here, we describe a novel bacterial strain, designated T37^T, which was isolated from the marine sediment of Xiaoshi Island, PR China. Growth of strain T37^T occurs at 15-40 °C (optimum 37 °C), pH 6.0-9.0 (optimum 7.5), and in the presence of 0.5-5.5% (w/v) NaCl (optimum 1.5%). Characteristic biochemical traits of the novel strain include MK-9 as the major menaquinone. The major fatty acids identified were iso-C_14:0 and C_16:1 ω9c (oleic acid). Phosphatidylethanolamine, phosphatidylglycerol, diphosphatidylglycerol, and phosphoglycolipids were the major cellular polar lipids. The G + C content of genomic DNA was 58.4 mol%. Unusual outer membrane features deduced from the analysis of cell morphology point towards the formation of an enlarged periplasmic space putatively used for the digestion of macromolecules. Phylogenetic analyses based on 16S rRNA genes and the genome indicated that strain T37^T represents a novel species and genus affiliated with a distinct family level lineage of the verrucomicrobial subdivision 1. Our polyphasic taxonomy approach places the novel strain in a new genus within the current family Verrucomicrobiaceae, order Verrucomicrobiales, class Verrucomicrobiae. Strain T37^T (= KCTC 72799 ^T = MCCC 1H00391^T) is the type strain of a novel species, for which the name Sulfuriroseicoccus oceanibius gen. nov., sp. nov. is proposed.

High-quality Arabidopsis thaliana Genome Assembly with Nanopore and HiFi Long Reads

Arabidopsis thaliana is an important and long-established model species for plant molecular biology, genetics, epigenetics, and genomics. However, the latest version of reference genome still contains significant number of missing segments. Here, we report a high-quality and almost complete Col-0 genome assembly with two gaps (Col-XJTU) using combination of Oxford Nanopore Technology ultra-long reads, PacBio high-fidelity long reads, and Hi-C data. The total genome assembly size is 133,725,193 bp, introducing 14.6 Mb of novel sequences compared to the TAIR10.1 reference genome. All five chromosomes of Col-XJTU assembly are highly accurate with consensus quality (QV) scores > 60 (ranging from 62 to 68), which are higher than those of TAIR10.1 reference (QV scores ranging from 45 to 52). We have completely resolved chromosome (Chr) 3 and Chr5 in a telomere-to-telomere manner. Chr4 has been completely resolved except the nucleolar organizing regions, which comprise long repetitive DNA fragments. The Chr1 centromere (CEN1), reportedly around 9 Mb in length, is particularly challenging to assemble due to the presence of tens of thousands of CEN180 satellite repeats. Using the cutting-edge sequencing data and novel computational approaches, we assembled about 4 Mb of sequence for CEN1 and a 3.5-Mb-long CEN2. We investigated the structure and epigenetics of centromeres. We detected four clusters of CEN180 monomers, and found that the centromere-specific histone H3-like protein (CENH3) exhibits a strong preference for CEN180 cluster 3. Moreover, we observed hypomethylation patterns in CENH3-enriched regions. We believe that this high-quality genome assembly, Col-XJTU, would serve as a valuable reference to better understand the global pattern of centromeric polymorphisms, as well as genetic and epigenetic features in plants.

Mosaic Arrangement of the 5S rDNA in the Aquatic Plant Landoltia punctata (Lemnaceae)

Duckweeds are a group of monocotyledonous aquatic plants in the Araceae superfamily, represented by 37 species divided into five genera. Duckweeds are the fastest growing flowering plants and are distributed around the globe; moreover, these plants have multiple applications, including biomass production, wastewater remediation, and making pharmaceutical proteins. Dotted duckweed (Landoltia punctata), the sole species in genus Landoltia, is one of the most resilient duckweed species. The ribosomal DNA (rDNA) encodes the RNA components of ribosomes and represents a significant part of plant genomes but has not been comprehensively studied in duckweeds. Here, we characterized the 5S rDNA genes in L. punctata by cloning and sequencing 25 PCR fragments containing the 5S rDNA repeats. No length variation was detected in the 5S rDNA gene sequence, whereas the nontranscribed spacer (NTS) varied from 151 to 524 bp. The NTS variants were grouped into two major classes, which differed both in nucleotide sequence and the type and arrangement of the spacer subrepeats. The dominant class I NTS, with a characteristic 12-bp TC-rich sequence present in 3-18 copies, was classified into four subclasses, whereas the minor class II NTS, with shorter, 9-bp nucleotide repeats, was represented by two identical sequences. In addition to these diverse subrepeats, class I and class II NTSs differed in their representation of cis-elements and the patterns of predicted G-quadruplex structures, which may influence the transcription of the 5S rDNA. Similar to related duckweed species in the genus Spirodela, L. punctata has a relatively low rDNA copy number, but in contrast to Spirodela and the majority of other plants, the arrangement of the 5S rDNA units demonstrated an unusual, heterogeneous pattern in L. punctata, as revealed by analyzing clones containing double 5S rDNA neighboring units. Our findings may further stimulate the research on the evolution of the plant rDNA and discussion of the molecular forces driving homogenization of rDNA repeats in concerted evolution.

Polycomb Repressive Complex 2 and KRYPTONITE regulate pathogen-induced programmed cell death in Arabidopsis

The Polycomb Repressive Complex 2 (PRC2) is well-known for its role in controlling developmental transitions by suppressing the premature expression of key developmental regulators. Previous work revealed that PRC2 also controls the onset of senescence, a form of developmental programmed cell death (PCD) in plants. Whether the induction of PCD in response to stress is similarly suppressed by the PRC2 remained largely unknown. In this study, we explored whether PCD triggered in response to immunity- and disease-promoting pathogen effectors is associated with changes in the distribution of the PRC2-mediated histone H3 lysine 27 trimethylation (H3K27me3) modification in Arabidopsis thaliana. We furthermore tested the distribution of the heterochromatic histone mark H3K9me2, which is established, to a large extent, by the H3K9 methyltransferase KRYPTONITE, and occupies chromatin regions generally not targeted by PRC2. We report that effector-induced PCD caused major changes in the distribution of both repressive epigenetic modifications and that both modifications have a regulatory role and impact on the onset of PCD during pathogen infection. Our work highlights that the transition to pathogen-induced PCD is epigenetically controlled, revealing striking similarities to developmental PCD.

Ancient Origin of Two 5S rDNA Families Dominating in the Genus Rosa and Their Behavior in the Canina-Type Meiosis

The genus Rosa comprises more than 100 woody species characterized by intensive hybridization, introgression, and an overall complex evolutionary history. Besides many diploid species (2n = 2x = 14) polyploids ranging from 3x to 10x are frequently found. Here we analyzed 5S ribosomal DNA in 19 species covering two subgenera and the major sections within subg. Rosa. In addition to diploids and polyploids with regular meiosis, we focused on 5x dogroses (Rosa sect. Caninae), which exhibit an asymmetric meiosis differentiating between bivalent- and univalent-forming chromosomes. Using genomic resources, we reconstructed 5S rDNA units to reveal their phylogenetic relationships. Additionally, we designed locus-specific probes derived from intergenic spacers (IGSs) and determined the position and number of 5S rDNA families on chromosomes. Two major 5S rDNA families (termed 5S_A and 5S_B, respectively) were found at variable ratios in both diploid and polyploid species including members of the early diverging subgenera, Rosa persica and Rosa minutifolia. Within subg. Rosa species of sect. Rosa amplified the 5S_A variant only, while taxa of other sections contained both variants at variable ratios. The 5S_B family was often co-localized with 35S rDNA at the nucleolar organizer regions (NOR) chromosomes, whereas the co-localization of the 5S_A family with NOR was only exceptionally observed. The allo-pentaploid dogroses showed a distinct distribution of 5S rDNA families between bivalent- and univalent-forming chromosomes. In conclusion, two divergent 5S rDNA families dominate rose genomes. Both gene families apparently arose in the early history of the genus, already 30 myrs ago, and apparently survived numerous speciation events thereafter. These observations are consistent with a relatively slow genome turnover in the Rosa genus.

Hybrid seed incompatibility in Capsella is connected to chromatin condensation defects in the endosperm

Hybridization of closely related plant species is frequently connected to endosperm arrest and seed failure, for reasons that remain to be identified. In this study, we investigated the molecular events accompanying seed failure in hybrids of the closely related species pair Capsella rubella and C. grandiflora. Mapping of QTL for the underlying cause of hybrid incompatibility in Capsella identified three QTL that were close to pericentromeric regions. We investigated whether there are specific changes in heterochromatin associated with interspecific hybridizations and found a strong reduction of chromatin condensation in the endosperm, connected with a strong loss of CHG and CHH methylation and random loss of a single chromosome. Consistent with reduced DNA methylation in the hybrid endosperm, we found a disproportionate deregulation of genes located close to pericentromeric regions, suggesting that reduced DNA methylation allows access of transcription factors to targets located in heterochromatic regions. Since the identified QTL were also associated with pericentromeric regions, we propose that relaxation of heterochromatin in response to interspecies hybridization exposes and activates loci leading to hybrid seed failure.

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.

Thousands of high-quality sequencing samples fail to show meaningful correlation between 5S and 45S ribosomal DNA arrays in humans

The ribosomal RNA genes (rDNA) are tandemly arrayed in most eukaryotes and exhibit vast copy number variation. There is growing interest in integrating this variation into genotype-phenotype associations. Here, we explored a possible association of rDNA copy number variation with autism spectrum disorder and found no difference between probands and unaffected siblings. Because short-read sequencing estimates of rDNA copy number are error prone, we sought to validate our 45S estimates. Previous studies reported tightly correlated, concerted copy number variation between the 45S and 5S arrays, which should enable the validation of 45S copy number estimates with pulsed-field gel-verified 5S copy numbers. Here, we show that the previously reported strong concerted copy number variation may be an artifact of variable data quality in the earlier published 1000 Genomes Project sequences. We failed to detect a meaningful correlation between 45S and 5S copy numbers in thousands of samples from the high-coverage Simons Simplex Collection dataset as well as in the recent high-coverage 1000 Genomes Project sequences. Our findings illustrate the challenge of genotyping repetitive DNA regions accurately and call into question the accuracy of recently published studies of rDNA copy number variation in cancer that relied on diverse publicly available resources for sequence data.

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.

Genome-wide scan for selection signatures reveals novel insights into the adaptive capacity in local North African cattle

Natural-driven selection is supposed to have left detectable signatures on the genome of North African cattle which are often characterized by the fixation of genetic variants associated with traits under selection pressure and/or an outstanding genetic differentiation with other populations at particular loci. Here, we investigate the population genetic structure and we provide a first outline of potential selection signatures in North African cattle using single nucleotide polymorphism genotyping data. After comparing our data to African, European and indicine cattle populations, we identified 36 genomic regions using three extended haplotype homozygosity statistics and 92 outlier markers based on Bayescan test. The 13 outlier windows detected by at least two approaches, harboured genes (e.g. GH1, ACE, ASIC3, HSPH1, MVD, BCL2, HIGD2A, CBFA2T3) that may be involved in physiological adaptations required to cope with environmental stressors that are typical of the North African area such as infectious diseases, extended drought periods, scarce food supply, oxygen scarcity in the mountainous areas and high-intensity solar radiation. Our data also point to candidate genes involved in transcriptional regulation suggesting that regulatory elements had also a prominent role in North African cattle response to environmental constraints. Our study yields novel insights into the unique adaptive capacity in these endangered populations emphasizing the need for the use of whole genome sequence data to gain a better understanding of the underlying molecular mechanisms.

Agro-Morphological, Microanatomical and Molecular Cytogenetic Characterization of the Medicinal Plant Chelidonium majus L

Chelidonium majus L. is a medicinal plant well-known as a valuable source of isoquinoline alkaloids, which has a variety of pharmacological properties including anti-viral and anti-bacterial effects. However, considerable intraspecific bio-morphological variability in C. majus complicates raw material identification and verification. For the first time, we have brought into cultivation five populations of C. majus subsp. majus originated from different regions, and performed their agro-morphological, microanatomical and molecular cytogenetic characterization. All examined populations produced high seed (18.6–19.9 kg/ha) and raw material (0.84–1.08 t/ha) yields; total alkaloid contents were within 0.30–0.38%. Nevertheless, significant differences in plant morphology and yield-contributing traits were observed. The performed microanatomical analysis of leaves and flowers in double- and normal-flowered plants revealed micro-diagnostic features (including tissue topography, types of stomata, laticifers, structure of leaf mesophyll, hairs, sepals and petals) important for identification of C. majus raw materials. The analysis of chromosome morphology, DAPI-banding patterns, FISH mapping of 45S and 5S rDNA and also chromosome behavior in meiosis allowed us to identify for the first time all chromosomes in karyotypes and confirm relative genotype stability of the studied plants. Our findings indicate that the examined C. majus populations can be used in further breeding programs.

Epigenetic silencing of clustered tRNA genes in Arabidopsis

Beyond their key role in translation, cytosolic transfer RNAs (tRNAs) are involved in a wide range of other biological processes. Nuclear tRNA genes (tDNAs) are transcribed by the RNA polymerase III (RNAP III) and cis-elements, trans-factors as well as genomic features are known to influence their expression. In Arabidopsis, besides a predominant population of dispersed tDNAs spread along the 5 chromosomes, some clustered tDNAs have been identified. Here, we demonstrate that these tDNA clusters are transcriptionally silent and that pathways involved in the maintenance of DNA methylation play a predominant role in their repression. Moreover, we show that clustered tDNAs exhibit repressive chromatin features whilst their dispersed counterparts contain permissive euchromatic marks. This work demonstrates that both genomic and epigenomic contexts are key players in the regulation of tDNAs transcription. The conservation of most of these regulatory processes suggests that this pioneering work in Arabidopsis can provide new insights into the regulation of RNA Pol III transcription in other organisms, including vertebrates.

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.

Field cress genome mapping: Integrating linkage and comparative maps with cytogenetic analysis for rDNA carrying chromosomes

Field cress (Lepidium campestre L.), despite its potential as a sustainable alternative oilseed plant, has been underutilized, and no prior attempts to characterize the genome at the genetic or molecular cytogenetic level have been conducted. Genetic maps are the foundation for anchoring and orienting annotated genome assemblies and positional cloning of candidate genes. Our principal goal was to construct a genetic map using integrated approaches of genetic, comparative and cytogenetic map analyses. In total, 503 F2 interspecific hybrid individuals were genotyped using 7,624 single nucleotide polymorphism markers. Comparative analysis demonstrated that ~57% of the sequenced loci in L. campestre were congruent with Arabidopsis thaliana (L.) genome and suggested a novel karyotype, which predates the ancestral crucifer karyotype. Aceto-orcein chromosome staining and fluorescence in situ hybridization (FISH) analyses confirmed that L. campestre, L. heterophyllum Benth. and their hybrids had a chromosome number of 2n = 2x = 16. Flow cytometric analysis revealed that both species possess 2C roughly 0.4 picogram DNA. Integrating linkage and comparative maps with cytogenetic map analyses assigned two linkage groups to their particular chromosomes. Future work could incorporate FISH utilizing A. thaliana mapped BAC clones to allow the chromosomes of field cress to be identified reliably.

Unique Epigenetic Features of Ribosomal RNA Genes (rDNA) in Early Diverging Plants (Bryophytes)

Introduction: In plants, the multicopy genes encoding ribosomal RNA (rDNA) typically exhibit heterochromatic features and high level of DNA methylation. Here, we explored rDNA methylation in early diverging land plants from Bryophyta (15 species, 14 families) and Marchantiophyta (4 species, 4 families). DNA methylation was investigated by methylation-sensitive Southern blot hybridization in all species. We also carried out whole genomic bisulfite sequencing in Polytrichum formosum (Polytrichaceae) and Dicranum scoparium (Dicranaceae) and used available model plant methyloms (Physcomitrella patents and Marchantia polymorpha) to determine rDNA unit-wide methylation patterns. Chromatin structure was analyzed using fluorescence in situ hybridization (FISH) and immunoprecipitation (CHIP) assays. Results: In contrast to seed plants, bryophyte rDNAs were efficiently digested with methylation-sensitive enzymes indicating no or low levels of CG and CHG methylation in these loci. The rDNA methylom analyses revealed variation between species ranging from negligible (<3%, P. formosum, P. patens) to moderate (7 and 17% in M. polymorpha and D. scoparium, respectively) methylation levels. There were no differences between coding and noncoding parts of rDNA units and between gametophyte and sporophyte tissues. However, major satellite repeat and transposable elements were heavily methylated in P. formosum and D. scoparium. In P. formosum rDNA, the euchromatic H3K4m3 and heterochromatic H3K9m2 histone marks were nearly balanced contrasting the angiosperms data where H3K9m2 typically dominates rDNA chromatin. In moss interphase nuclei, rDNA was localized at the nucleolar periphery and its condensation level was high. Conclusions: Unlike seed plants, the rRNA genes seem to escape global methylation machinery in bryophytes. Distinct epigenetic features may be related to rDNA expression and the physiology of these early diverging plants that exist in haploid state for most of their life cycles.

The multi-faceted regulation of nuclear tRNA gene transcription

Transfer RNAs are among the most ancient molecules of life on earth. Beyond their crucial role in protein synthesis as carriers of amino acids, they are also important players in a plethora of other biological processes. Many debates in term of biogenesis, regulation and function persist around these fascinating non-coding RNAs. Our review focuses on the first step of their biogenesis in eukaryotes, i.e. their transcription from nuclear genes. Numerous and complementary ways have emerged during evolution to regulate transfer RNA gene transcription. Here, we will summarize the different actors implicated in this process: cis-elements, trans-factors, genomic contexts, epigenetic environments and finally three-dimensional organization of nuclear genomes. © 2019 IUBMB Life, 2019 © 2019 IUBMB Life, 71(8):1099-1108, 2019.

An Arabidopsis FANCJ helicase homologue is required for DNA crosslink repair and rDNA repeat stability

Proteins of the Fanconi Anemia (FA) complementation group are required for crosslink (CL) repair in humans and their loss leads to severe pathological phenotypes. Here we characterize a homolog of the Fe-S cluster helicase FANCJ in the model plant Arabidopsis, AtFANCJB, and show that it is involved in interstrand CL repair. It acts at a presumably early step in concert with the nuclease FAN1 but independently of the nuclease AtMUS81, and is epistatic to both error-prone and error-free post-replicative repair in Arabidopsis. The simultaneous knock out of FANCJB and the Fe-S cluster helicase RTEL1 leads to induced cell death in root meristems, indicating an important role of the enzymes in replicative DNA repair. Surprisingly, we found that AtFANCJB is involved in safeguarding rDNA stability in plants. In the absence of AtRTEL1 and AtFANCJB, we detected a synergetic reduction to about one third of the original number of 45S rDNA copies. It is tempting to speculate that the detected rDNA instability might be due to deficiencies in G-quadruplex structure resolution and might thus contribute to pathological phenotypes of certain human genetic diseases.

Competition for Chiasma Formation Between Identical and Homologous (But Not Identical) Chromosomes in Synthetic Autotetraploids of Arabidopsis thaliana

Polyploid organisms provide additional opportunities to study meiosis in a more complex context since more than two potential homologous chromosomes are available. When the chromosome complement of a diploid individual is duplicated, each chromosome is accompanied by one identical and two homologous chromosomes within the same nucleus. In this situation, a competition in pairing/synapsis/chiasma formation between identical and homologous (but not necessarily identical) chromosomes can occur. Several studies have been conducted in different species to address whether there are preferences in crossover formation between identical rather than homologous chromosomes. In this study, multivalent and chiasma frequencies were cytologically analyzed in synthetic autotetraploids of Arabidopsis thaliana including the accessions Col, Ler, and the Col/Ler hybrid. Fluorescence in situ hybridization was conducted to identify each chromosome at metaphase I. The new Col and Ler tetraploids showed high multivalent frequencies, exceeding the theoretical 66.66% expected on a simple random end-pairing model, thus indicating that there are more than two autonomous synaptic sites per chromosome despite their small size. However, a significant excess of bivalent pairs was found in the Col/Ler hybrid, mainly due to the contribution of chromosomes 2 and 3. The mean chiasma frequencies of the three artificial autotetraploids were about twofold the corresponding mean cell chiasma frequencies of their diploid counterparts. The relative contribution of each chromosome to the total chiasma frequency was similar in the three genotypes, with the exception of a lower contribution of chromosome 3 in the hybrid. Preferences for chiasma formation between identical and homologous chromosomes were analyzed in Col/Ler 4x, taking advantage of the cytological differences between the accessions: variations in the size of the 45S rDNA region on the short arm of chromosome 2 and changes in the size and localization of the 5S rDNA region in chromosome 3. We observed a different behavior of chromosomes 2 and 3, i.e., random chiasma formation between identical and homologous chromosomes 2, and preferences for chiasma formation between homologous chromosomes 3. Hence, our results reveal the existence of chromosome-specific mechanisms responsible for these preferences.