The effect of polyploidy and hybridization on the evolution of floral colour in Nicotiana (Solanaceae)

Background and Aims Speciation in angiosperms can be accompanied by changes in floral colour that may influence pollinator preference and reproductive isolation. This study investigates whether changes in floral colour can accompany polyploid and homoploid hybridization, important processes in angiosperm evolution.Methods Spectral reflectance of corolla tissue was examined for 60 Nicotiana (Solanaceae) accessions (41 taxa) based on spectral shape (corresponding to pigmentation) as well as bee and hummingbird colour perception in order to assess patterns of floral colour evolution. Polyploid and homoploid hybrid spectra were compared with those of their progenitors to evaluate whether hybridization has resulted in floral colour shifts.Key Results Floral colour categories in Nicotiana seem to have arisen multiple times independently during the evolution of the genus. Most younger polyploids displayed an unexpected floral colour, considering those of their progenitors, in the colour perception of at least one pollinator type, whereas older polyploids tended to resemble one or both of their progenitors.Conclusions Floral colour evolution in Nicotiana is weakly constrained by phylogeny, and colour shifts do occur in association with both polyploid and homoploid hybrid divergence. Transgressive floral colour in N. tabacum has arisen by inheritance of anthocyanin pigmentation from its paternal progenitor while having a plastid phenotype like its maternal progenitor. Potentially, floral colour evolution has been driven by, or resulted in, pollinator shifts. However, those polyploids that are not sympatric (on a regional scale) with their progenitor lineages are typically not divergent in floral colour from them, perhaps because of a lack of competition for pollinators.     

Pseudogenes contribute to the extreme diversity of nuclear ribosomal DNA in the hard coral Acropora

One characteristic of Indo-Pacific Acropora spp. is extremely high diversity in rDNA sequences at both the species and individual levels. In order to test the hypothesis that pseudogenes may contribute to this diversity, three kinds of analyses were conducted. First, for A. millepora (the species containing the most diverse suite of rDNA types), RT-PCR was used to determine which 5.8S rDNA types are expressed. Second, as previous studies have indicated that interspecific hybridization has occurred in the genus Acropora and silencing of rDNA loci via nucleolar dominance has been shown in some cases to involve methylation, patterns of variation were examined at methylation-susceptible sites. Third, patterns of substitution at conserved sites (including those that are likely to contribute to secondary structure in rRNA) in the 5.8S rDNA were examined. These analyses consistently indicated that one rDNA sequence type present in a broad range of Indo-Pacific Acropora species is likely to consist predominantly of pseudogenes. Patterns of variation also suggest that species may differ with respect to which rDNA sequence types have been silenced and which are active. These pseudogenes are likely to have arisen as a consequence of the introduction of highly divergent rDNA types into single genomes by interspecific hybridization events, and we attribute the extreme rDNA diversity characteristic of many Acropora species to both the independent evolution of these silenced rDNA types and to the suppressive effects of high sequence diversity on homogenization processes acting on functional loci.     

Kinetics of cold inactivation and heat reactivation of the ribulose diphosphate carboxylase activity of crystalline Fraction I proteins isolated from different species and hybrids of Nicotiana

The maximum specific ribulose diphosphate carboxylase activityof crystalline Fraction I proteins isolated from 8 species,and 12 reciprocal, interspecific, F₁ hybrids of Nicotiana changedto minimum specific activity between 16 and 20 hr after theenzymes at 25?C were placed in ice baths. The minimum statewas preserved for more than one year without further changewhen crystals of enzymes were suspended in buffer at 0?C. Atany time during cold storage, the maximum state of specificactivity could be regained by dissolving the crystals in buffercontaining NaCl and heating the solution at temperatures upto 50?C. Maximum activity was attained by heating at 50?C for20 min whereas more than 100 hr was required at 25?C. An equilibrium was demonstrated between maximally active andminimally active enzyme molecules. The equilibrium constantwas a function of the storage temperatures between 0? and 25?C.The molar enthalpy, H, separating maximally active from minimallyactive enzyme molecules was 3.75?0.1 kcal/mole, an energy differenceconsistent with previous findings that the change in specificactivity did not result from dissociation of subunits or depolymerizationof aggregates. The effect of cold storage of Nicotiana enzyme would seem tooffer a molecular basis for explaining the effect of cold temperatureson photosynthetic CO₂ fixation by winter wheat leaves whichwas shown by Sawada and Miyachi (5) to produce an accumulationof ribulose 1–5 diphosphate and diversion of the normalpath of photosynthetic CO₂ fixation away from phosphoglycericacid to the appearance of fixed carbon in serine and glycineby the glycolate path. (Received October 22, 1973; )     

Species-specific evolution of telomeric and rDNA repeats in the tobacco composite genome

In order to investigate possible interactions between parental genomes in the composite genome of Nicotiana tabacum we have analyzed the organization of telomeric (TTTAGGG)_n and ribosomal gene (rDNA) repeats in the progenitor genomes Nicotiana sylvestris and Nicotiana tomentosiformis or Nicotiana otophora. Telomeric arrays in the Nicotiana species tested are heterogeneous in length ranging from 20 to 200 kb in N. sylvestris, from 20 to 50 kb in N. tomentosiformis, from 15 to 100kb in N. otophora, and from 40 to 160kb in N. tabacum. The patterns of rDNA repeats (18S, 5.8S, 25S RNA) appeared to be highly homogeneous and speciesspecific; no parental rDNA units corresponding to N. sylvestris, N. tomentosiformis or N. otophora were found in the genome of N. tabacum by Southern hybridization. The results provide evidence for a species-specific evolution of telomeric and ribosomal repeats in the tobacco composite genome.     

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

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

Karyotype diversification and evolution in diploid and polyploid South American Hypochaeris (Asteraceae) inferred from rDNA localization and genetic fingerprint data

Background and Aims: Changes in chromosome structure and number play an importantrole in plant evolution. A system well-suited to studying differentmodes of chromosome evolution is the genus Hypochaeris (Asteraceae)with its centre of species' diversity in South America. AllSouth American species uniformly have a chromosome base numberof x = 4 combined with variation in rDNA number and distribution,and a high frequency of polyploidy. The aim of this paper isto assess directions and mechanisms of karyotype evolution inSouth American species by interpreting both newly obtained andprevious data concerning rDNA localization in a phylogeneticcontext. Methods: Eleven Hypochaeris species from 18 populations were studiedusing fluorescence in situ hybridization (FISH) with 35S and5S rDNA probes. A phylogenetic framework was established fromneighbour-net analysis of amplified fragment length polymorphism(AFLP) fingerprint data. Key Results: A single 5S rDNA locus is invariably found on the short armof chromosome 2. Using 35S rDNA loci, based on number (one ortwo) and localization (interstitial on the long arm of chromosome2, but sometimes lacking, and terminal or interstitial on theshort arm of chromosome 3, only very rarely lacking), sevenkaryotype groups can be distinguished; five of these includepolyploids. Karyotype groups with more than one species do notform monophyletic groups. Conclusions: Early evolution of Hypochaeris in South America was characterizedby considerable karyotype differentiation resulting from independentderivations from an ancestral karyotype. There was marked diversificationwith respect to the position and evolution of the 35S rDNA locuson chromosome 3, probably involving inversions and/or transpositions,and on chromosome 2 (rarely 3) concerning inactivation and loss.Among these different karyotype assemblages, the apargioidesgroup and its derivatives constitute by far the majority ofspecies.     

Intragenic Hill-Robertson interference influences selection intensity on synonymous mutations in Drosophila

Natural selection influences synonymous mutations and synonymouscodon usage in many eukaryotes to improve the efficiency oftranslation in highly expressed genes. Recent studies of genecomposition in eukaryotes have shown that codon usage also variesindependently of expression levels, both among genes and atthe intragenic level. Here, we investigate rates of evolution(Ks) and intensity of selection (_s) on synonymous mutationsin two groups of genes that differ greatly in the length oftheir exons, but with equivalent levels of gene expression andrates of crossing-over in Drosophila melanogaster. We estimate_s using patterns of divergence and polymorphism in 50 Drosophilagenes (100 kb of coding sequence) to take into account possiblevariation in mutation trends across the genome, among genesor among codons. We show that genes with long exons exhibithigher Ks and reduced _s compared to genes with short exons.We also show that Ks and _s vary significantly across long exons,with higher Ks and reduced _s in the central region comparedto flanking regions of the same exons, hence indicating thatthe difference between genes with short and long exons can bemostly attributed to the central region of these long exons.Although amino acid composition can also play a significantrole when estimating Ks and _s, our analyses show that the differencesin Ks and _s between genes with short and long exons and acrosslong exons cannot be explained by differences in protein composition.All these results are consistent with the Interference Selection(IS) model that proposes that the Hill-Robertson (HR) effectcaused by many weakly selected mutations has detectable evolutionaryconsequences at the intragenic level in genomes with recombination.Under the IS model, exon size and exon-intron structure influencethe effectiveness of selection, with long exons showing reducedeffectiveness of selection when compared to small exons andthe central region of long exons showing reduced intensity ofselection compared to flanking coding regions. Finally, ourresults further stress the need to consider selection on synonymousmutations and its variation—among and across genes andexons—in studies of protein evolution.     

The phylogenetic potential of entire 26S rDNA sequences in plants

18S ribosomal RNA genes are the most widely used nuclear sequences for phylogeny reconstruction at higher taxonomic levels in plants. However, due to a conservative rate of evolution, 18S rDNA alone sometimes provides too few phylogenetically informative characters to resolve relationships adequately. Previous studies using partial sequences have suggested the potential of 26S or large-subunit (LSU) rDNA for phylogeny retrieval at taxonomic levels comparable to those investigated with 18S rDNA. Here we explore the patterns of molecular evolution of entire 26S rDNA sequences and their impact on phylogeny retrieval. We present a protocol for PCR amplification and sequencing of entire (approximately 3.4 kb) 26S rDNA sequences as single amplicons, as well as primers that can be used for amplification and sequencing. These primers proved useful in angiosperms and Gnetales and likely have broader applicability. With these protocols and primers, entire 26S rDNA sequences were generated for a diverse array of 15 seed plants, including basal eudicots, monocots, and higher eudicots, plus two representatives of Gnetales. Comparisons of sequence dissimilarity indicate that expansion segments (or divergence domains) evolve 6.4 to 10.2 times as fast as conserved core regions of 26S rDNA sequences in plants. Additional comparisons indicate that 26S rDNA evolves 1.6 to 2.2 times as fast as and provides 3.3 times as many phylogenetically informative characters as 18S rDNA; compared to the chloroplast gene rbcL, 26S rDNA evolves at 0.44 to 1.0 times its rate and provides 2.0 times as many phylogenetically informative characters. Expansion segment sequences analyzed here evolve 1.2 to 3.0 times faster than rbcL, providing 1.5 times the number of informative characters. Plant expansion segments have a pattern of evolution distinct from that found in animals, exhibiting less cryptic sequence simplicity, a lower frequency of insertion and deletion, and greater phylogenetic potential.      

PAML 4: phylogenetic analysis by maximum likelihood

PAML, currently in version 4, is a package of programs for phylogeneticanalyses of DNA and protein sequences using maximum likelihood(ML). The programs may be used to compare and test phylogenetictrees, but their main strengths lie in the rich repertoire ofevolutionary models implemented, which can be used to estimateparameters in models of sequence evolution and to test interestingbiological hypotheses. Uses of the programs include estimationof synonymous and nonsynonymous rates (d_N and d_S) between twoprotein-coding DNA sequences, inference of positive Darwinianselection through phylogenetic comparison of protein-codinggenes, reconstruction of ancestral genes and proteins for molecularrestoration studies of extinct life forms, combined analysisof heterogeneous data sets from multiple gene loci, and estimationof species divergence times incorporating uncertainties in fossilcalibrations. This note discusses some of the major applicationsof the package, which includes example data sets to demonstratetheir use. The package is written in ANSI C, and runs underWindows, Mac OSX, and UNIX systems. It is available at http://abacus.gene.ucl.ac.uk/software/paml.html.     

The Close Relationship Between the A and B Genomes in Avena L. (Poaceae) Determined by Molecular Cytogenetic Analysis of Total Genomic, Tandemly and Dispersed Repetitive DNA Sequences

The genusAvena L. (Poaceae) consists of diploid, tetraploid,and hexaploid species, with the B genome known only in tetraploidspecies and the D genome in the hexaploid species. DNA:DNAinsitu hybridization, using total genomic DNA from diploidAvenastrigosa Schreb. (A_sgenome) as a probe, labelled all 28 chromosomesof the AB tetraploidAvena vaviloviana (Malz.) Mordv. stronglyand uniformly, revealing the close relationship between thesetwo genomes. Comparison of patterns of size-separated DNA restrictionfragments between the diploidA. strigosa and the tetraploidA.vaviloviana , using 32 different restriction enzymes, revealedno differences. Southern hybridization using total AB genomicDNA as a probe also gave no differences in banding patternsbetween the two genomes, even when a large excess of A genomicDNA was used as a block. From anA. vaviloviana genomic library,1800 colonies were blotted and probed sequentially with A andAB genomic DNA, but no colony was identified to be B genomespecific. DNA digests of AB genome tetraploids with restrictionenzymeHae III gave a strong band at 4.2 kb. Clone pAbKB3, derivedfrom the 4.2 kb band, was found to be part of aTy1-copia -likeretrotransposon present in A and B genome chromosomes. ClonedrRNA genes were used forin situ hybridization and showed thatdiploidA. strigosa has four major sites for 18S-25S rDNA andtwo pairs of sites for 5S rDNA (pairs on the same satellitedchromosome, on different chromosome arms), while 4xA. vavilovianahas eight major sites for 18S-25S rDNA and four pairs of sitesfor 5S rDNA (pairs on the same satellited chromosome, on differentchromosome arms). A repetitive sequence from rye pSc119.2, showeddispersed hybridization, while the telomeric sequence in clonepLT11 hybridized to telomeres. Again no discrimination was possiblebetween A and B genome chromosomes. The molecular similaritiesbetween the diploidA. strigosa and thebarbata group tetraploidsclearly indicate that thebarbata group of tetraploids arosefrom A_sdiploids through autotetraploidization. Avena ; evolution; repetitive sequences; in situ hybridization; retrotransposons; genome organization     

Mapping of ribosomal DNA and (TTAGGG)n telomeric sequence by FISH in the bivalve Patinopecten yessoensis (Jay, 1857)

The chromosome set of Patinopecten yessoensis (Jay, 1857) wascharacterized using Giemsa staining, DAPI staining and fluorescencein situ hybridization (FISH) with three repetitive DNA probes[18S–28S rDNA, 5S rDNA and telomeric (TTAGGG)_n]. DAPIstaining showed that AT-rich regions were located on the centromereof almost all chromosomes and interstitial banding was not observed.FISH showed that 18S–28S rDNA spread over the short armsof two subtelocentric chromosome pairs and 5S rDNA was locatedon the long arm of one subtelocentric chromosome pair. SequentialFISH demonstrated that 18S–28S and 5S rDNA were locatedon different chromosomes. FISH also showed that the vertebratetelomeric sequence (TTAGGG)_n was located on both ends of eachchromosome and no interstitial signals were detected. Sequential18S–28S rDNA and (TTAGGG)_n FISH indicated that repeatedunits of the two multicopy families were closely associatedon the same chromosome pair. (Received 4 January 2007; accepted 1 September 2007)     

Limitation to Photosynthesis in Water-stressed Leaves: Stomata vs. Metabolism and the Role of ATP

Decreasing relative water content (RWC) of leaves progressivelydecreases stomatal conductance (g_s), slowing CO₂ assimilation(A) which eventually stops, after which CO₂ is evolved. In somestudies, photosynthetic potential (A_pot), measured under saturatingCO₂, is unaffected by a small loss of RWC but becomes progressivelymore inhibited, and less stimulated by elevated CO₂, below athreshold RWC (Type 1 response). In other studies, A_pot andthe stimulation of A by elevated CO₂ decreases progressivelyas RWC falls (Type 2 response). Decreased A_pot is caused byimpaired metabolism. Consequently, as RWC declines, the relativelimitation of A by g_s decreases, and metabolic limitation increases.Causes of decreased A_pot are considered. Limitation of ribulosebisphosphate (RuBP) synthesis is the likely cause of decreasedA_pot at low RWC, not inhibition or loss of photosynthetic carbonreduction cycle enzymes, including RuBP carboxylase/oxygenase(Rubisco). Limitation of RuBP synthesis is probably caused byinhibition of ATP synthesis, due to progressive inactivationor loss of Coupling Factor resulting from increasing ionic (Mg²⁺)concentration, not to reduced capacity for electron or protontransport, or inadequate trans-thylakoid proton gradient (pH).Inhibition of A_pot by accumulation of assimilates or inadequateinorganic phosphate is not considered significant. DecreasedATP content and imbalance with reductant status affect cellmetabolism substantially: possible consequences are discussedwith reference to accumulation of amino acids and alterationsin protein complement under water stress.     

Diverse evolutionary pathways shaped 5S rDNA of species of tribe Anemoneae (Ranunculaceae) and reveal phylogenetic signal

Anemone sensu lato (including Pulsatilla and Hepatica), tribe Anemoneae (Ranunculaceae), is arranged into two subgenera, Anemone and Anemonidium, with basic chromosome numbers x = 8 and x = 7, respectively. We elucidated the level of divergence of 5S rDNA unit arrays between the subgenera, determined intra‐individual and interspecific sequence variation and tested 5S rDNA phylogenetic signal in revealing the origin of polyploid species. High intra‐individual nucleotide diversity and the presence of 5S rDNA unit array length variants and pseudogenes indicate that weak homogenization forces have shaped 5S rDNA in the investigated species. Our results show that 5S rDNA evolved through two major changes: diversification of 5S rDNA into two lineages, one with long (subgenus Anemone) and one with short 5S rDNA unit arrays (subgenus Anemonidium); and subsequent contraction and expansion of 5S rDNA unit arrays. Phylogenetic analysis based on 5S rDNA supports the hypothesis that A. parviflora could be a parental species and donor of the subgenome D to the allopolyploids A. multifida (BBDD) and A. baldensis (AABBDD). In A. baldensis interlocus exchange possibly occurred, followed by subsequent replacement of the 5S rDNA from subgenome D with those from subgenome B. Here we present evidence that both models, concerted and birth‐and‐death evolution, were probably involved in the evolution of the 5S rDNA multigene family in subgenera Anemone and Anemonidium.     

Concerted evolution of 18–5.8–26S rDNA repeats in Nicotiana allotetraploids

We review and extend data showing concerted evolution of parental 18–5.8–26S nuclear ribosomal DNA (18–26S rDNA) gene families in three natural Nicotiana allotetraploids ( N. tabacum , N. rustica and N. arentsii , each 2 n  = 4 x  = 48) and one synthetic N. tabacum line (Th37, ♀ N. sylvestris (2 n  = 24) × ♂ N. tomentosiformis (2 n  = 24)). The origin of the gene families was analysed by sequence polymorphisms in the intergenic spacer (IGS) region and the number of chromosomal loci by fluorescence in situ hybridization (FISH). FISH revealed that the number and locations of 18–26S rDNA in the natural allopolyploids was the sum of those found in the diploid progenitors. However, the rDNA restriction patterns showed polymorphisms in the IGS that were not additive, suggesting that parental rDNA clusters were partially ( N. tabacum, N. rustica ) or completely ( N. arentsii ) overwritten by hybrid-specific units. Thus the Nicotiana allotetraploids show evidence of concerted evolution, including both intralocus and interlocus gene conversion. A feral N. tabacum collected in Bolivia had a higher proportion of unconverted parental rDNA units than cultivated tobacco varieties, suggesting either that rDNA homogenization is accelerated by inbreeding or multiple origins of tobacco. There is no evidence for the elimination of N. sylvestris- derived rDNA units in the synthetic Th37 tobacco line as occurred in natural tobacco, although several novel rDNA unit variants were found in most but not all the hybrid plants. Factors that may control the occurrence and extent of rDNA homogenization are discussed for allopolyploids in Nicotiana and other taxa.  © 2004 The Linnean Society of London, Biological Journal of the Linnean Society , 2004, 82 , 615–625.     

Recent and recurrent polyploidy in Tragopogon (Asteraceae): cytogenetic, genomic and genetic comparisons

Tragopogon mirus Ownbey and T. miscellus Ownbey are allopolyploids that formed repeatedly during the past 80 years following the introduction of three diploids (T. dubius Scop., T. pratensis L. and T. porrifolius L.) from Europe to western North America. These polyploid species of known parentage are useful for studying the consequences of recent and recurrent polyploidization. We summarize recent analyses of the cytogenetic, genomic and genetic consequences of polyploidy in Tragopogon. Analyses of rDNA ITS (internal transcribed spacer) + ETS (external transcribed spacer) sequence data indicate that the parental diploids are phylogenetically well separated within Tragopogon (a genus of perhaps 150 species), in agreement with isozymic and cpDNA data. Using Southern blot and cloning experiments on tissue from early herbarium collections of T. mirus and T. miscellus (from 1949) to represent the rDNA repeat condition closer to the time of polyploidization than samples collected today, we have demonstrated concerted evolution of rDNA. Concerted evolution is ongoing, but has not proceeded to completion in any polyploid population examined; rDNA repeats of the diploid T. dubius are typically lost or converted in both allopolyploids, including populations of independent origin. Molecular cytogenetic studies employing rDNA probes, as well as centromeric and subtelomeric repeats isolated from Tragopogon, distinguished all chromosomes among the diploid progenitors (2n = 12). The diploid chromosome complements are additive in both allopolyploids (2n = 24); there is no evidence of major chromosomal rearrangements in populations of either T. mirus or T. miscellus. cDNA‐AFLP display revealed differences in gene expression between T. miscellus and its diploid parents, as well as between populations of T. miscellus of reciprocal origin. Approximately 5% of the genes examined in the allopolyploid populations have been silenced, and an additional 4% exhibit novel gene expression relative to their diploid parents. Some of the differences in gene expression represent maternal or paternal effects. Multiple origins of a polyploid species not only affect patterns of genetic variation in natural populations, but also contribute to differential patterns of gene expression and may therefore play a major role in the long‐term evolution of polyploids. © 2004 The Linnean Society of London, Biological Journal of the Linnean Society, 2004, 82 , 485–501.     

Phylogenetic analysis of Hordeum (Poaceae) as inferred by nuclear rDNA ITS sequences

Hordeum (Poaceae, Triticeae) occurs with 31 species worldwide in temperate regions, with the exception of Australasia. About 50% of the species are polyploids (4x, 6x) or occur as di- and polyploid forms. To analyze the phylogenetic relationships among diploid and polyploid taxa of the genus the nuclear rDNA internal transcribed spacer region (ITS) was analyzed for 91 accessions, representing all Hordeum species, together with 10 outgroup species. PCR products were either directly sequenced (outgroups) or cloned and eight clones per individual were analyzed. Phylogenetic analysis revealed four major clades that concur with the four genome groups in Hordeum (H, I, Xa, and Xu). Allopolyploids, putative autopolyploids, and species groups within the closely related H-genome clade could be identified. The ITS data indicate times of independent evolution of paralogous rDNA clusters on different chromosomes intermitted by sweeps of homogenization among these clusters and bi-directional homogenization of the clusters in diploids. Penalized likelihood analysis revealed an age of about 12 million years (my) for the genus and indicated the start of a rapid radiation in the H-genome group about 2.5 my ago in South America and Asia.     

Silenced rRNA genes are activated and substitute for partially eliminated active homeologs in the recently formed allotetraploid,Tragopogon mirus (Asteraceae)

To study the relationship between uniparental rDNA (encoding 18S, 5.8S and 26S ribosomal RNA) silencing (nucleolar dominance) and rRNA gene dosage, we studied a recently emerged (within the last 80 years) allotetraploid Tragopogon mirus (2n=24), formed from the diploid progenitors T. dubius (2n=12, D-genome donor) and T. porrifolius (2n=12, P-genome donor). Here, we used molecular, cytogenetic and genomic approaches to analyse rRNA gene activity in two sibling T. mirus plants (33A and 33B) with widely different rRNA gene dosages. Plant 33B had ~400 rRNA genes at the D-genome locus, which is typical for T. mirus, accounting for ~25% of total rDNA. We observed characteristic expression dominance of T. dubius-origin genes in all organs. Its sister plant 33A harboured a homozygous macrodeletion that reduced the number of T. dubius-origin genes to about 70 copies (~4% of total rDNA). It showed biparental rDNA expression in root, flower and callus, but not in leaf where D-genome rDNA dominance was maintained. There was upregulation of minor rDNA variants in some tissues. The RNA polymerase I promoters of reactivated T. porrifolius-origin rRNA genes showed reduced DNA methylation, mainly at symmetrical CG and CHG nucleotide motifs. We hypothesise that active, decondensed rDNA units are most likely to be deleted via recombination. The silenced homeologs could be used as a ‘first reserve'' to ameliorate mutational damage and contribute to evolutionary success of polyploids. Deletion and reactivation cycles may lead to bidirectional homogenisation of rRNA arrays in the long term.