Incomplete homogenization of 18 S ribosomal DNA coding regions in Arabidopsis thaliana

Cytogenomic analysis unveils mixed molecular evolution and recurrent chromosomal rearrangements shaping the multigene families on Schistocerca grasshopper genomes

Multigene families are essential components of eukaryotic genomes and play key roles either structurally and functionally. Their modes of evolution remain elusive even in the era of genomics, because multiple multigene family sequences coexist in genomes, particularly in large repetitive genomes. Here, we investigate how the multigene families 18S rDNA, U2 snDNA, and H3 histone evolved in 10 species of Schistocerca grasshoppers with very large and repeat-enriched genomes. Using sequenced genomes and fluorescence in situ hybridization mapping, we find substantial differences between species, including the number of chromosomal clusters, changes in sequence abundance and nucleotide composition, pseudogenization, and association with transposable elements (TEs). The intragenomic analysis of Schistocerca gregaria using long-read sequencing and genome assembly unveils conservation for H3 histone and recurrent pseudogenization for 18S rDNA and U2 snDNA, likely promoted by association with TEs and sequence truncation. Remarkably, TEs were frequently associated with truncated copies, were also among the most abundant in the genome, and revealed signatures of recent activity. Our findings suggest a combined effect of concerted and birth-and-death models driving the evolution of multigene families in Schistocerca over the last 8 million years, and the occurrence of intra- and interchromosomal rearrangements shaping their chromosomal distribution. Despite the conserved karyotype in Schistocerca, our analysis highlights the extensive reorganization of repetitive DNAs in Schistocerca, contributing to the advance of comparative genomics for this important grasshopper genus.

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.

Age-Dependent and Tissue-Specific Alterations in the rDNA Clusters of the Panax ginseng C. A. Meyer Cultivated Cell Lines

Long-term cultivation of Panax ginseng cell lines leads to a decreasing synthesis of the biologically active substances used in traditional medicine. To gain insight into the cellular mechanisms which may influence this process, we analyzed variations within the rDNA cluster of the Oriental ginseng cell lines. The cell lines were cultivated for 6 and 24 years; the number of nucleoli and chromosomes was analyzed. The complete 18S rDNA sequences were cloned and sequenced. The nucleotide polymorphism and phylogenetic relations of the sequences were analyzed, and the secondary structures for separate 18S rRNA regions were modeled. The 18S rDNA accumulated mutations during cell cultivation that correlate well with an increase in the number of chromosomes and nucleoli. The patterns of nucleotide diversity are culture-specific and the increasing polymorphism associates with cytosine methylation sites. The secondary structures of some 18S rRNA regions and their interaction can alter during cultivation. The phylogenetic tree topologies are particular for each cell line.The observed alterations in rDNA clusters are associated with a somaclonal variation, leading to changes in the pattern of intracellular synthesis during cell cultivation. The identified divergent rRNAs could provide additional gene expression regulation in P. ginseng cells by forming heterogeneous ribosomes.

Nucleolar Dominance in a Tetraploidy Hybrid Lineage Derived From Carassius auratus red var. () × Megalobrama amblycephala ()

Nucleolar dominance is related to the expression of 45S rRNA genes inherited from one progenitor due to the silencing of the other progenitor’s rRNA genes. To investigate nucleolar dominance associated with tetraploidization, we analyzed the changes regarding the genetic traits and expression of 45S rRNA genes in tetraploidy hybrid lineage including F₁ allotetraploids (4n = 148) and F₂ autotetraploids (4n = 200) derived from the distant hybridization of Carassius auratus red var. (2n = 100) () ×Megalobrama amblycephala (2n = 48) (). Results showed that nucleolar dominance from the females was established in F₁ hybrids and it was inherited in F₂ hybrids, suggesting that tetraploidization can lead to rapid establishment of nucleolar dominance in the hybrid origin’s tetraploid lineage. These results extend the knowledge of nucleolar dominance in polyploidy hybrid animals, which are of significance for the evolution of hybrids in vertebrates.

Evidence of non-tandemly repeated rDNAs and their intragenomic heterogeneity in Rhizophagus irregularis

Arbuscular mycorrhizal fungus (AMF) species are some of the most widespread symbionts of land plants. Our much improved reference genome assembly of a model AMF, Rhizophagus irregularis DAOM-181602 (total contigs = 210), facilitated a discovery of repetitive elements with unusual characteristics. R. irregularis has only ten or 11 copies of complete 45S rDNAs, whereas the general eukaryotic genome has tens to thousands of rDNA copies. R. irregularis rDNAs are highly heterogeneous and lack a tandem repeat structure. These findings provide evidence for the hypothesis that rDNA heterogeneity depends on the lack of tandem repeat structures. RNA-Seq analysis confirmed that all rDNA variants are actively transcribed. Observed rDNA/rRNA polymorphisms may modulate translation by using different ribosomes depending on biotic and abiotic interactions. The non-tandem repeat structure and intragenomic heterogeneity of AMF rDNA/rRNA may facilitate successful adaptation to various environmental conditions, increasing host compatibility of these symbiotic fungi.

Evolutionary trends in animal ribosomal DNA loci: introduction to a new online database

Ribosomal DNA (rDNA) loci encoding 5S and 45S (18S-5.8S-28S) rRNAs are important components of eukaryotic chromosomes. Here, we set up the animal rDNA database containing cytogenetic information about these loci in 1343 animal species (264 families) collected from 542 publications. The data are based on in situ hybridisation studies (both radioactive and fluorescent) carried out in major groups of vertebrates (fish, reptiles, amphibians, birds, and mammals) and invertebrates (mostly insects and mollusks). The database is accessible online at www.animalrdnadatabase.com . The median number of 45S and 5S sites was close to two per diploid chromosome set for both rDNAs despite large variation (1-74 for 5S and 1-54 for 45S sites). No significant correlation between the number of 5S and 45S rDNA loci was observed, suggesting that their distribution and amplification across the chromosomes follow independent evolutionary trajectories. Each group, irrespective of taxonomic classification, contained rDNA sites at any chromosome location. However, the distal and pericentromeric positions were the most prevalent (> 75% karyotypes) for 45S loci, while the position of 5S loci was more variable. We also examined potential relationships between molecular attributes of rDNA (homogenisation and expression) and cytogenetic parameters such as rDNA positions, chromosome number, and morphology.

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

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

Haplotype Detection from Next-Generation Sequencing in High-Ploidy-Level Species: 45S rDNA Gene Copies in the Hexaploid Spartina maritima

Gene and whole-genome duplications are widespread in plant nuclear genomes, resulting in sequence heterogeneity. Identification of duplicated genes may be particularly challenging in highly redundant genomes, especially when there are no diploid parents as a reference. Here, we developed a pipeline to detect the different copies in the ribosomal RNA gene family in the hexaploid grass Spartina maritima from next-generation sequencing (Roche-454) reads. The heterogeneity of the different domains of the highly repeated 45S unit was explored by identifying single nucleotide polymorphisms (SNPs) and assembling reads based on shared polymorphisms. SNPs were validated using comparisons with Illumina sequence data sets and by cloning and Sanger (re)sequencing. Using this approach, 29 validated polymorphisms and 11 validated haplotypes were reported (out of 34 and 20, respectively, that were initially predicted by our program). The rDNA domains of S. maritima have similar lengths as those found in other Poaceae, apart from the 5′-ETS, which is approximately two-times longer in S. maritima. Sequence homogeneity was encountered in coding regions and both internal transcribed spacers (ITS), whereas high intragenomic variability was detected in the intergenic spacer (IGS) and the external transcribed spacer (ETS). Molecular cytogenetic analysis by fluorescent in situ hybridization (FISH) revealed the presence of one pair of 45S rDNA signals on the chromosomes of S. maritima instead of three expected pairs for a hexaploid genome, indicating loss of duplicated homeologous loci through the diploidization process. The procedure developed here may be used at any ploidy level and using different sequencing technologies.

The translation initiation factor eIF3i up-regulates vascular endothelial growth factor A, accelerates cell proliferation, and promotes angiogenesis in embryonic development and tumorigenesis

Vascular endothelial growth factor A (VEGFA) is a critical proangiogenic factor that is activated by hypoxia at both the transcriptional and post-transcriptional levels. In hypoxia conditions, stabilized hypoxia-inducible factor 1α (HIF1A) is the key regulator for transcriptional activation of VEGFA. However, the post-transcriptional control of VEGFA expression remains poorly understood. Here, we report that the eukaryotic translation initiation factor 3i (eIF3i) is required for VEGFA protein expression in both normal embryonic and tumorigenic angiogenesis. eIF3i is dynamically expressed in the early stages of zebrafish embryogenesis and in human hepatocellular carcinoma tissues. eIF3i homozygous mutant zebrafish embryos show severe angiogenesis defects and human hepatocellular cancer cells with depletion of eIF3i to induce less angiogenesis in tumor models. Under hypoxia, the HIF1A protein can interact with its binding sequence in the eIF3i promoter and activate eIF3i transcription. The expression of VEGFA, which should rise in hypoxia, is significantly inhibited by eIF3i siRNA treatment. Moreover, eIF3i knockdown did not cause a general translation repression but specifically reduced the translation efficiency of the VEGFA mRNAs. Taken together, our results suggest that eIF3i is induced by HIF1A under hypoxia and controls normal and tumorigenic angiogenesis through regulating VEGFA protein translation.

RAPD, SCAR and conserved 18S rDNA markers for a red-listed and endemic medicinal plant species, Knema andamanica (Myristicaceae)

Knema andamanica is a red-listed endemic medicinal species of Myristicaceae restricted to Andaman and Nicobar (A&N) Islands, India. This species is used in tribal medicines and has immense bioprospective potential. With a view to generate suitable genomic markers for classification and identification, we have generated RAPD, SCAR and conserved 18S rDNA markers from K. andamanica. A unique 585 bp fragment, that distinguished it from seven other related species of Myristicaceae was first amplified using the random primer OPE 06 and converted to SCAR marker (GenBank accession # JN228256). The conserved sequences of 18S rDNA loci from K. andamanica were also amplified and sequenced (GenBank accession #JN228265). The sequence revealed deviations including 18 variable regions and 15 indels that were unique to K. andamanica. These markers can help in definite identification of K. andamanica even at the juvenile stages.

Next generation sequencing analysis reveals a relationship between rDNA unit diversity and locus number in Nicotiana diploids

BACKGROUND

Tandemly arranged nuclear ribosomal DNA (rDNA), encoding 18S, 5.8S and 26S ribosomal RNA (rRNA), exhibit concerted evolution, a pattern thought to result from the homogenisation of rDNA arrays. However rDNA homogeneity at the single nucleotide polymorphism (SNP) level has not been detailed in organisms with more than a few hundred copies of the rDNA unit. Here we study rDNA complexity in species with arrays consisting of thousands of units.

METHODS

We examined homogeneity of genic (18S) and non-coding internally transcribed spacer (ITS1) regions of rDNA using Roche 454 and/or Illumina platforms in four angiosperm species, Nicotiana sylvestris, N. tomentosiformis, N. otophora and N. kawakamii. We compared the data with Southern blot hybridisation revealing the structure of intergenic spacer (IGS) sequences and with the number and distribution of rDNA loci.

RESULTS AND CONCLUSIONS

In all four species the intragenomic homogeneity of the 18S gene was high; a single ribotype makes up over 90% of the genes. However greater variation was observed in the ITS1 region, particularly in species with two or more rDNA loci, where >55% of rDNA units were a single ribotype, with the second most abundant variant accounted for >18% of units. IGS heterogeneity was high in all species. The increased number of ribotypes in ITS1 compared with 18S sequences may reflect rounds of incomplete homogenisation with strong selection for functional genic regions and relaxed selection on ITS1 variants. The relationship between the number of ITS1 ribotypes and the number of rDNA loci leads us to propose that rDNA evolution and complexity is influenced by locus number and/or amplification of orphaned rDNA units at new chromosomal locations.

Highly conserved low-copy nuclear genes as effective markers for phylogenetic analyses in angiosperms

Organismal phylogeny provides a crucial evolutionary framework for many studies and the angiosperm phylogeny has been greatly improved recently, largely using organellar and rDNA genes. However, low-copy protein-coding nuclear genes have not been widely used on a large scale in spite of the advantages of their biparental inheritance and vast number of choices. Here, we identified 1083 highly conserved low-copy nuclear genes by genome comparison. Furthermore, we demonstrated the use of five nuclear genes in 91 angiosperms representing 46 orders (73% of orders) and three gymnosperms as outgroups for a highly resolved phylogeny. These nuclear genes are easy to clone and align, and more phylogenetically informative than widely used organellar genes. The angiosperm phylogeny reconstructed using these genes was largely congruent with previous ones mainly inferred from organellar genes. Intriguingly, several new placements were uncovered for some groups, including those among the rosids, the asterids, and between the eudicots and several basal angiosperm groups. These conserved universal nuclear genes have several inherent qualities enabling them to be good markers for reconstructing angiosperm phylogeny, even eukaryotic relationships, further providing new insights into the evolutionary history of angiosperms.