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.

Dancing together and separate again: gymnosperms exhibit frequent changes of fundamental 5S and 35S rRNA gene (rDNA) organisation.

In higher eukaryotes, the 5S rRNA genes occur in tandem units and are arranged either separately (S-type arrangement) or linked to other repeated genes, in most cases to rDNA locus encoding 18S-5.8S-26S genes (L-type arrangement). Here we used Southern blot hybridisation, PCR and sequencing approaches to analyse genomic organisation of rRNA genes in all large gymnosperm groups, including Coniferales, Ginkgoales, Gnetales and Cycadales. The data are provided for 27 species (21 genera). The 5S units linked to the 35S rDNA units occur in some but not all Gnetales, Coniferales and in Ginkgo (∼30% of the species analysed), while the remaining exhibit separate organisation. The linked 5S rRNA genes may occur as single-copy insertions or as short tandems embedded in the 26S-18S rDNA intergenic spacer (IGS). The 5S transcript may be encoded by the same (Ginkgo, Ephedra) or opposite (Podocarpus) DNA strand as the 18S-5.8S-26S genes. In addition, pseudogenised 5S copies were also found in some IGS types. Both L- and S-type units have been largely homogenised across the genomes. Phylogenetic relationships based on the comparison of 5S coding sequences suggest that the 5S genes independently inserted IGS at least three times in the course of gymnosperm evolution. Frequent transpositions and rearrangements of basic units indicate relatively relaxed selection pressures imposed on genomic organisation of 5S genes in plants.

Making a functional diploid: from polysomic to disomic inheritance.

One little understood feature of polyploid speciation is the transition from polysomic to disomic inheritance, and much recent attention has focused on the role of pairing genes in this process. Using computer simulations we studied the effects of mutations, chromosomal inversions, chiasma, neofunctionalization, subfunctionalization and selection on the evolution of disomic inheritance in a polyploid over 10 000 generations. We show that: the evolution of pairing genes is not essential for the establishment of disomic inheritance, as genetic drift, coupled with a threshold for homologue pairing fidelity, is sufficient to explain the transition from polysomic to disomic inheritance; high rates of recombination increase the number of generations required for disomic inheritance to become established; both neofunctionalization and subfunctionalization speed up the transition to disomic inheritance. The data suggest that during polyploid species establishment, selection will favour reduced chiasma number and/or more focused distribution. The data also suggest a new role for subfunctionalization in that it can drive disomic inheritance. The evolution of subfunctionalization in genes across the genome will then act to maintain genes in syntenic blocks and may explain why such regions are so highly conserved.

Mobilization of retrotransposons in synthetic allotetraploid tobacco.

Allopolyploidy is a major driving force in plant evolution and can induce rapid structural changes in the hybrid genome. As major components of plant genomes, transposable elements are involved in these changes. In a previous work, we observed turnover of retrotransposon insertions in natural allotretraploid tobacco (Nicotiana tabacum). Here, we studied the early stages of allopolyploid formation by monitoring changes at retrotransposon insertion sites in the Th37 synthetic tobacco. We used sequence-specific amplification polymorphism (SSAP) to study insertion patterns of two populations of the Tnt1 retrotransposon in Th37 S4 generation plants, and characterized the nature of polymorphic insertion sites. We observed significant amplification of young Tnt1 populations. Newly transposed copies were amplified from maternal elements and were highly similar to Tnt1A tobacco copies amplified in response to microbial factors. A high proportion of paternal SSAP bands were not transmitted to the hybrid, corresponding to various rearrangements at paternal insertion sites, including indels or the complete loss of the Tnt1/flanking junction. These data indicate that major changes, such as retrotransposon amplification and molecular restructuring in or around insertion sites, occur rapidly in response to allopolyploidy.

Frequent silencing of rDNA loci on the univalent-forming genomes contrasts with their stable expression on the bivalent-forming genomes in polyploid dogroses (Rosa sect. Caninae).

The polyploid species in Rosa section Caninae (2n=21, 28 or 35) are characterized by an unusual reproductive system known as odd (or asymmetric) meiosis. Only two chromosome sets form bivalents in meiosis, whereas the remaining chromosomes are transmitted as univalents through the female germline. Evolution of ribosomal rRNA genes (rDNA) does not seem to be significantly affected by interlocus homogenization in dogroses. As a consequence, most species contain several rDNA families falling into two main clades (beta and gamma) thought to be differentially distributed between bivalent and univalent chromosomes, respectively. Here, we have investigated expression of rRNA gene families in five pentaploid species (R. canina, R. rubiginosa, R. dumalis, R. sherardii and R. caesia, 2n=35) and in one tetraploid (R. mollis, 2n=28). Using extensive sequencing of ITS clones and cleaved amplified polymorphism sequence (CAPS) analysis, we found that the beta-family was constitutively expressed in all species. However, there was large variation in the expression patterns of families constituting the gamma-clade. In addition, a single family can be active in one species, whereas silenced in another. The data show that the families on bivalent-forming chromosomes dominate rDNA expression in all dogrose species. We hypothesize that genes on bivalent genomes are stably expressed, whereas those on univalent genomes undergo variable levels of epigenetic silencing. Nonetheless, mosaic expression of univalent genomes could contribute to phenotypic variation between the species.

Reduced inducibility of SOCS3 by interferon gamma associates with higher resistance of human breast cancer lines as compared to normal mammary epithelial cells.

The resistance to interferons (IFNs) limits their anticancer therapeutic efficacy. Here we studied the antiproliferative effect of interferon gamma in relation to SOCS3 expression in a panel of breast cancer cell lines and normal mammary epithelial cells. Compared to normal cells most breast cancer lines (7/8) were highly resistant to IFN-gamma. Using Northern blot and real time RT-PCR we investigated transcription of SOCS3 genes. All normal epithelial cells (4/4) showed SOCS3 induction (2-14 fold) while most breast cancer lines did not or weakly activated SOCS3 after the interferon gamma treatment. Among the cancer lines, the MDA-MB-468 cells showed increased sensitivity to IFN-gamma and relatively high level of SOCS3 induction (2-3 fold). Together, there was a good correlation

The asymmetric meiosis in pentaploid dogroses (Rosa sect. Caninae) is associated with a skewed distribution of rRNA gene families in the gametes.

In pentaploid dogroses, Rosa section Caninae (2n=5x=35), the pollen transmits one basic genome (x=7) derived from the seven segregating bivalents, whereas the egg transmits four basic genomes (4x=28) one set derived from the segregation of seven bivalents and three sets of univalent-forming chromosomes. Chromosomes from all five genomes carry 18-5.8-26S nuclear ribosomal DNA (rDNA) sites. This mode of sexual reproduction, known as permanent odd polyploidy, can potentially lead to the independent evolution of rDNA on bivalent- and univalent-forming chromosomes. To test this hypothesis, we analyzed rRNA gene families in pollen and somatic leaf tissue of R. canina, R. rubiginosa and R. dumalis. Six major rRNA gene families (alpha, beta, beta' gamma, delta and epsilon) were identified based on several highly polymorphic sites in the internal transcribed spacers (ITSs). At least two of the major rRNA gene families were found in each species indicating that rDNAs have not been homogenized across subgenomes. A comparison of ITS1 sequences from leaf and pollen showed differences: the shared beta rRNA gene family was more abundant among pollen clones compared to leaf clones and must constitute a major part of the rDNA loci on bivalent-forming chromosomes. The gamma and delta families were underrepresented in pollen genomes and are probably located predominantly (or solely) on the univalents. The results support the hypothesis that pentaploid dogroses inherited a bivalent-forming genome from a common proto-canina ancestor, a likely donor of the beta rDNA family. Allopolyploidy with distantly related species is likely to have driven evolution of Rosa section Caninae.

The ups and downs of genome size evolution in polyploid species of Nicotiana (Solanaceae).

BACKGROUND: In studies looking at individual polyploid species, the most common patterns of genomic change are that either genome size in the polyploid is additive (i.e. the sum of parental genome donors) or there is evidence of genome downsizing. Reports showing an increase in genome size are rare. In a large-scale analysis of 3008 species, genome downsizing was shown to be a widespread biological response to polyploidy. Polyploidy in the genus Nicotiana (Solanaceae) is common with approx. 40 % of the approx. 75 species being allotetraploid. Recent advances in understanding phylogenetic relationships of Nicotiana species and dating polyploid formation enable a temporal dimension to be added to the analysis of genome size evolution in these polyploids. METHODS: Genome sizes were measured in 18 species of Nicotiana (nine diploids and nine polyploids) ranging in age from <200,000 years to approx. 4.5 Myr old, to determine the direction and extent of genome size change following polyploidy. These data were combined with data from genomic in situ hybridization and increasing amounts of information on sequence composition in Nicotiana to provide insights into the molecular basis of genome size changes. KEY RESULTS AND CONCLUSIONS: By comparing the expected genome size of the polyploid (based on summing the genome size of species identified as either a parent or most closely related to the diploid progenitors) with the observed genome size, four polyploids showed genome downsizing and five showed increases. There was no discernable pattern in the direction of genome size change with age of polyploids, although with increasing age the amount of genome size change increased. In older polyploids (approx. 4.5 million years old) the increase in genome size was associated with loss of detectable genomic in situ hybridization signal, whereas some hybridization signal was still detected in species exhibiting genome downsizing. The possible significance of these results is discussed.

Rapid concerted evolution of nuclear ribosomal DNA in two Tragopogon allopolyploids of recent and recurrent origin.

We investigated concerted evolution of rRNA genes in multiple populations of Tragopogon mirus and T. miscellus, two allotetraploids that formed recurrently within the last 80 years following the introduction of three diploids (T. dubius, T. pratensis, and T. porrifolius) from Europe to North America. Using the earliest herbarium specimens of the allotetraploids (1949 and 1953) to represent the genomic condition near the time of polyploidization, we found that the parental rDNA repeats were inherited in roughly equal numbers. In contrast, in most present-day populations of both tetraploids, the rDNA of T. dubius origin is reduced and may occupy as little as 5% of total rDNA in some individuals. However, in two populations of T. mirus the repeats of T. dubius origin outnumber the repeats of the second diploid parent (T. porrifolius), indicating bidirectional concerted evolution within a single species. In plants of T. miscellus having a low rDNA contribution from T. dubius, the rDNA of T. dubius was nonetheless expressed. We have apparently caught homogenization of rDNA repeats (concerted evolution) in the act, although it has not proceeded to completion in any allopolyploid population yet examined.

Molecular cytogenetics and tandem repeat sequence evolution in the allopolyploid Nicotiana rustica compared with diploid progenitors N. paniculata and N. undulata.

Nicotiana rustica (2n = 4x = 48) is a natural allotetraploid composed of P and U genomes which are closely related to genomes of diploid species N. paniculata and N. undulata. Genomic in situ hybridization (GISH) also confirms that the diploid parents, or close relatives, are the ancestors of N. rustica. In order to study genetic interactions between ancestral genomes in the allotetraploid, we isolated three families of repetitive sequences, two from N. paniculata (NPAMBE and NPAMBO) and one from N. undulata (NUNSSP). Southern blot hybridization revealed that the sequences are digested with a range of restriction enzymes into regular ladder patterns indicating a tandem arrangement of high copy repeats possessing monomeric units of about 180 bp. The three-tandem sequences belong to a larger Nicotiana tandem repeat family called here the HRS-60 family. Members of this family are found in all Nicotiana species studied. Fluorescence in situ hybridization (FISH) analysis localized the satellite repeats to subtelomeric regions of most chromosomes of N. paniculata and N. undulata. The pattern of sequence distribution on the P- and U-genomes of N. rustica was similar to the putative parents N. paniculata and N. undulata respectively. However, NPAMBO repeats appear to be reduced and rearranged in N. rustica that may suggest evolution within the P genome. GISH and FISH with the tandem repeat probes failed to reveal intergenomic translocations as might be predicted from the nucleocytoplasmic interaction hypothesis.

Dynamic changes in the distribution of a satellite homologous to intergenic 26-18S rDNA spacer in the evolution of Nicotiana.

An approximately 135-bp sequence called the A1/A2 repeat was isolated from the transcribed region of the 26-18S rDNA intergenic spacer (IGS) of Nicotiana tomentosiformis. Fluorescence in situ hybridization (FISH) and Southern blot analysis revealed its occurrence as an independent satellite (termed an A1/A2 satellite) outside of rDNA loci in species of Nicotiana section Tomentosae. The chromosomal location, patterns of genomic dispersion, and copy numbers of its tandemly arranged units varied between the species. In more distantly related Nicotiana species the A1/A2 repeats were found only at the nucleolar organizer regions (NOR). There was a trend toward the elimination of the A1/A2 satellite in N. tabacum (tobacco), an allotetraploid with parents closely related to the diploids N. sylvestris and N. tomentosiformis. This process may have already commenced in an S(3) generation of synthetic tobacco. Cytosine residues in the IGS were significantly hypomethylated compared with the A1/A2 satellite. There was no clear separation between the IGS and satellite fractions in sequence analysis of individual clones and we found no evidence for CG suppression. Taken together the data indicate a dynamic nature of the A1/A2 repeats in Nicotiana genomes, with evidence for recurrent integration, copy number expansions, and contractions.

Variability in CpNpG methylation in higher plant genomes.

The methylation status of ribosomal gene (rRNA) clusters have been investigated in a large variety of angiosperm species. Here we have analysed methylation in ribosomal gene (rRNA) clusters using MspI, HpaII, BstNI, EcoRII and CfoI restriction enzymes in combination with Southern hybridization to the 25S rDNA probe. It was shown that cytosine methylation at CpG dinucleotides and CpNpG trinucleotides occurred in all plant genomes examined. Methylation of rDNA units at CpG dinucleotides (studied with CfoI) was high in all species tested with approx. 40-70% of units being completely or nearly completely methylated. In contrast, the extent of the CpNpG methylation (studied with MspI and EcoRII) varied significantly between species; the percentage of the rDNA fraction entirely methylated at CpNpG trinucleotides ranged from less than 1% to almost 90% depending on the genome studied. Larger interspecies than within species variation was also observed among several non-transcribing repetitive sequences. In a small genome of A. thaliana, the CpNpG methylation appeared to be highly compartmentalized into the repetitive fraction. The methylation of trinucleotides was abundant in large A+T-rich genomes and it is proposed that the CpA(T)pG trinucleotides may help to maintain a high density of methylatable targets in plant repeated sequences.

Telomerase activity in plant cells.

Telomerase is a ribonucleoprotein enzyme which elongates the G-rich strand of telomeric DNA to compensate for the progressive reduction in its length due to incomplete replication of chromosome ends, which in human somatic cells leads to cell cycle arrest upon shortening of telomeres to a critical length. To examine the possible involvement of telomerase in metabolism of plant genetic material, we used cells of Nicotiana tabacum strain TBY-2, a stable long-term culture which has kept a constant pattern of restriction fragments from chromosome termini during its 6 month period of cultivation in our laboratory. In a direct assay for telomerase, a 5' end-labeled plant telomeric oligonucleotide 5' (TTTAGGG)(3')6 was elongated in a TBY-2 cell extract, showing a pausing pattern which is a characteristic feature of telomerases from other organisms. The elongation was inhibited by RNase A pretreatment of the extract. We conclude that plant cells possess telomerase which is used for maintenance of their telomeres.

The telomeric sequence is directly attached to the HRS60 subtelomeric tandem repeat in tobacco chromosomes.

PCR and primers derived from the telomeric repeat (CCCTAAA)n and from the tobacco subtelomeric tandemly repetitive sequence HRS60 (EMBL X12489) were used to amplify the region linking the two loci. A 131 bp PCR product was obtained both from total tobacco DNA and from the DNA fraction enriched for telomeres. Its sequence only consists of the telomeric primer and the attached region of the HRS60 repetitive unit up to the end of the sequence complementary to the HRS60 primer. The site of direct continuity between the two sequences is formed by a (dA)7 tract.

Chromatin fragmentation associated with apoptotic changes in tobacco cells exposed to cold stress.

Programmed cell death (PCD) may be triggered by a variety of environmental stimuli. In this report we show that low temperature treatment of tobacco BY-2 cells results in specific chromatin changes. The early stage was characterised by chromatin condensation associated with specific endonucleolytic cleavage of the genome into fragments of 50-100 kbp in size. Later, after 2 weeks of the cold treatment, a ladder of nucleosomal units (178 bp) and their multiples occurred. Chromatin changes were accompanied by a general decrease in cell viability. However, the cell culture retained about 11% of living cells even after prolonged incubation in the cold suggesting the presence of a cold-resistant population of cells. The results support the view that PCD was activated by the cold stress. We suggest that cold-stressed tobacco BY-2 culture might be a useful system for investigation of PCD in plant cells.

Sequence-specific hypomethylation of the tobacco genome induced with dihydroxypropyladenine, ethionine and 5-azacytidine.

Higher plant DNA is methylated at CG and CNG targets. In this study we have investigated the tobacco methylation system in tissue culture using the methylation inhibitors 5-azacytidine (5-azaC), dihydroxypropyladenine (DHPA) and ethionine (Ethi), and methylation-sensitive restriction endonucleases HpaII, MspI, HhaI, EcoRII, ScrFI, and Fnu4HI. Surprisingly, CAG/CTG sequences, contrary to CG doublets and CCG/CGG triplets, appeared to be refractory to the inhibitory effect of 5-azaC. Thus 5-azaC cannot be considered a general inhibitor of DNA methylation in tobacco cells. On the other hand, DHPA, the inhibitor of S-adenosylhomocysteine (SAH) hydrolase, and Ethi caused hypomethylation of both CAG/CTG and CCG/CGG triplets but not of the CG doublets. The sensitivity of triplet-specific methylation to the inhibition of SAH hydrolase suggests the possibility that plant-specific DNA methylation at CNG targets might be modulated by alterations of the SAH/S-Adenosylmethionine ratio in plant cells.

Cross-talk between posttranscriptionally silenced neomycin phosphotransferase II transgenes.

Tobacco plants containing a transgene locus with two chimeric neomycin phosphotransferase II (nptII) genes in tail-to-tail orientation (locus 1) show posttranscriptional gene silencing. The silenced nptII transgenes of locus 1 can downregulate the expression of homologous nptII transgenes in hybrid plants. The 3' region of the silenced nptII genes located in the center of the inverted repeat locus 1 is extensively methylated. Moreover, 3' segments of in trans-inactivated transgenes also become methylated, suggesting cross-talk between homologous posttranscriptionally silenced genes. Our results are in accordance with the hypothesis that this cross-talk can be mediated by specially featured RNAs.

Drug-induced hypomethylation of a posttranscriptionally silenced transgene locus of tobacco leads to partial release of silencing.

The effect of DNA methylation upon posttranscriptional gene silencing (PTGS) has been investigated in transgenic tobacco lines showing PTGS and methylation of the neomycin phosphotransferase II (nptII) reporter genes. Application of the hypomethylation drugs dihydroxypropyladenine or 5-azacytidine resulted in approximately 30% reduced methylation of cytosines located in a non-symmetrical context in the 3' untranslated region of the nptII transgenes. The hypomethylation was accompanied by up to 12-fold increase in NPTII protein levels, suggesting that methylation of non-symmetrical motifs may account for an increased degree of PTGS. Models for the possible role of DNA methylation in PTGS are discussed.

Hypomethylation of CNG targets induced with dihydroxypropyladenine is rapidly reversed in the course of mitotic cell division in tobacco.

We followed the mitotic transmission of an experimentally induced hypomethylated state of several tobacco repetitive sequences in callus culture and plants. The initial hypomethylation was induced by a hypomethylation drug, dihydroxypropyladenine (DHPA), the competitive inhibitor of cellular S-adenosylhomocysteine hydrolase, which is known to preferentially inhibit methylation at CNG and non-symmetrical motifs while having a negligible effect on methylation at CG motifs. The deprivation of this drug resulted in an almost immediate remethylation of cytosines at CNG motifs ( MspI and EcoRII sites) leading us to conclude that, the hypomethylation effect of dihydroxypropyladenine is rather transient and differs from that of 5-azacytidine which often induces heritable changes in methylation patterns. The results suggest that de novo methylation of CNG motifs is a rapid and meiotically independent process on DNA sequences with pre-existing CG methylation.

DNA curvature of the tobacco GRS repetitive sequence family and its relation to nucleosome positioning.

Recently, a highly repetitive DNA sequence family (GRS) from tobacco was described in our laboratory. These sequences were found to be localized predominantly in the pericentromeric heterochromatin of tobacco chromosomes. To test the hypothesis that these sequences play an important role in the formation of heterochromatin, we investigated the DNA curvature of the GRS sequences and its possible impact to the chromatin structure at these loci. Application of the nearest-neighbour wedge model of intrinsic DNA curvature for the GRS1 family member predicted two loci of curvature: a major bend at the 5' end of the sequence and a minor bend of opposite direction at the centre of the GRS1. The presence of the major and the minor loci of DNA curvature was studied experimentally using permutation analysis and site-directed mutagenesis. The experimental results were consistent with the computer predictions. We gave evidence that the described DNA curvature is also present in the entire GRS family. Genomic statistical sequencing showed the conservation of the major bend sequence determinants in the members of the GRS family. To investigate the chromatin structure at the GRS sequences, we determined the nucleosome positioning in vivo at these sequences using thermal cycle primer extension. A relation between the curvature pattern and the histone octamer position was observed: the major bend is excluded from the nucleosome surface to the linker region, while the minor bend is distributed along the core DNA. The suggestion is made that the sequences in the minor locus of curvature define the rotational setting of the nucleosome, and a possible role of the major bend as a factor, which defines the translational setting, is discussed.