Characterization of a new family of tobacco highly repetitive DNA, GRS, specific for the Nicotiana tomentosiformis genomic component

The nature and organization of satellite DNAs in Petunia hybrida, related, and ancestral genomes

Introduction

The garden petunia, Petunia hybrida (Solanaceae) is a fertile, diploid, annual hybrid species (2n=14) originating from P. axillaris and P. inflata 200 years ago. To understand the recent evolution of the P. hybrida genome, we examined tandemly repeated or satellite sequences using bioinformatic and molecular cytogenetic analysis.

Methods

Raw reads from available genomic assemblies and survey sequences of P. axillaris N (PaxiN), P. inflata S6, (PinfS6), P. hybrida (PhybR27) and the here sequenced P. parodii S7 (PparS7) were used for graph and k-mer based cluster analysis of TAREAN and RepeatExplorer. Analysis of repeat specific monomer lengths and sequence heterogeneity of the major tandem repeat families with more than 0.01% genome proportion were complemented by fluorescent in situ hybridization (FISH) using consensus sequences as probes to chromosomes of all four species.

Results

Seven repeat families, PSAT1, PSAT3, PSAT4, PSAT5 PSAT6, PSAT7 and PSAT8, shared high consensus sequence similarity and organisation between the four genomes. Additionally, many degenerate copies were present. FISH in P. hybrida and in the three wild petunias confirmed the bioinformatics data and gave corresponding signals on all or some chromosomes. PSAT1 is located at the ends of all chromosomes except the 45S rDNA bearing short arms of chromosomes II and III, and we classify it as a telomere associated sequence (TAS). It is the most abundant satellite repeat with over 300,000 copies, 0.2% of the genomes. PSAT3 and the variant PSAT7 are located adjacent to the centromere or mid-arm of one to three chromosome pairs. PSAT5 has a strong signal at the end of the short arm of chromosome III in P. axillaris and P.inflata, while in P. hybrida additional interstitial sites were present. PSAT6 is located at the centromeres of chromosomes II and III. PSAT4 and PSAT8 were found with only short arrays.

Discussion

These results demonstrate that (i) repeat families occupy distinct niches within chromosomes, (ii) they differ in the copy number, cluster organization and homogenization events, and that (iii) the recent genome hybridization in breeding P. hybrida preserved the chromosomal position of repeats but affected the copy number of repetitive DNA.

The C4 protein encoded by tomato leaf curl Yunnan virus reverses transcriptional gene silencing by interacting with NbDRM2 and impairing its DNA-binding ability

In plants, cytosine DNA methylation is an efficient defense mechanism against geminiviruses, since methylation of the viral genome results in transcriptional gene silencing (TGS). As a counter-defense mechanism, geminiviruses encode viral proteins to suppress viral DNA methylation and TGS. However, the molecular mechanisms by which viral proteins contribute to TGS suppression remain incompletely understood. In this study, we found that the C4 protein encoded by tomato leaf curl Yunnan virus (TLCYnV) suppresses methylation of the viral genome through interacting with and impairing the DNA-binding ability of NbDRM2, a pivotal DNA methyltransferase in the methyl cycle. We show that NbDRM2 catalyzes the addition of methyl groups on specific cytosine sites of the viral genome, hence playing an important role in anti-viral defense. Underscoring the relevance of the C4-mediated suppression of NbDRM2 activity, plants infected by TLCYnV producing C4(S43A), a point mutant version of C4 unable to interact with NbDRM2, display milder symptoms and lower virus accumulation, concomitant with enhanced viral DNA methylation, than plants infected by wild-type TLCYnV. Expression of TLCYnV C4, but not of the NbDRM2-interaction compromised C4(S43A) mutant, in 16c-TGS Nicotiana benthamiana plants results in the recovery of GFP, a proxy for suppression of TGS. This study provides new insights into the molecular mechanisms by which geminiviruses suppress TGS, and uncovers a new viral strategy based on the inactivation of the methyltransferase NbDRM2.

A reference genome for Nicotiana tabacum enables map-based cloning of homeologous loci implicated in nitrogen utilization efficiency

BACKGROUND

Tobacco (Nicotiana tabacum) is an important plant model system that has played a key role in the early development of molecular plant biology. The tobacco genome is large and its characterisation challenging because it is an allotetraploid, likely arising from hybridisation between diploid N. sylvestris and N. tomentosiformis ancestors. A draft assembly was recently published for N. tabacum, but because of the aforementioned genome complexities it was of limited utility due to a high level of fragmentation.

RESULTS

Here we report an improved tobacco genome assembly, which, aided by the application of optical mapping, achieves an N50 size of 2.17 Mb and enables anchoring of 64% of the genome to pseudomolecules; a significant increase from the previous value of 19%. We use this assembly to identify two homeologous genes that explain the differentiation of the burley tobacco market class, with potential for greater understanding of Nitrogen Utilization Efficiency and Nitrogen Use Efficiency in plants; an important trait for future sustainability of agricultural production.

CONCLUSIONS

Development of an improved genome assembly for N. tabacum enables what we believe to be the first successful map-based gene discovery for the species, and demonstrates the value of an improved assembly for future research in this model and commercially-important species.

Complete Chloroplast Genome of Nicotiana otophora and its Comparison with Related Species

Nicotiana otophora is a wild parental species of Nicotiana tabacum, an interspecific hybrid of Nicotiana tomentosiformis and Nicotiana sylvestris. However, N. otophora is least understood as an alternative paternal donor. Here, we compared the fully assembled chloroplast (cp) genome of N. otophora and with those of closely related species. The analysis showed a cp genome size of 156,073 bp and exhibited a typical quadripartite structure, which contains a pair of inverted repeats separated by small and large single copies, containing 163 representative genes, with 165 microsatellites distributed unevenly throughout the genome. Comparative analysis of a gene with known function across Nicotiana species revealed 76 protein-coding sequences, 20 tRNA sequences, and 3 rRNA sequence shared between the cp genomes. The analysis revealed that N. otophora is a sister species to N. tomentosiformis within the Nicotiana genus, and Atropha belladonna and Datura stramonium are their closest relatives. These findings provide a valuable analysis of the complete N. otophora cp genome, which can identify species, elucidate taxonomy, and reconstruct the phylogeny of genus Nicotiana.

Evaluation of DNA bending models in their capacity to predict electrophoretic migration anomalies of satellite DNA sequences

DNA containing a sequence that generates a local curvature exhibits a pronounced retardation in electrophoretic mobility. Various theoretical models have been proposed to explain relationship between DNA structural features and migration anomaly. Here, we studied the capacity of 15 static wedge-bending models to predict electrophoretic behavior of 69 satellite monomers derived from four divergent families. All monomers exhibited retarded mobility in PAGE corresponding to retardation factors ranging 1.02-1.54. The curvature varied both within and across the groups and correlated with the number, position, and lengths of A-tracts. Two dinucleotide models provided strong correlation between gel mobility and curvature prediction; two trinucleotide models were satisfactory while remaining dinucleotide models provided intermediate results with reliable prediction for subsets of sequences only. In some cases, similarly shaped molecules exhibited relatively large differences in mobility and vice versa. Generally less accurate predictions were obtained in groups containing less homogeneous sequences possessing distinct structural features. In conclusion, relatively universal theoretical models were identified suitable for the analysis of natural sequences known to harbor relatively moderate curvature. These models could be potentially applied to genome wide studies. However, in silico predictions should be viewed in context of experimental measurement of intrinsic DNA curvature.

Tobacco karyotyping by accurate centromere identification and novel repetitive DNA localization

Tobacco (Nicotiana tabacum) is an amphidiploid species (2n = 4x = 48, genome constitution SSTT) derived from a natural hybrid between Nicotiana sylvestris (2n = 2x = 24, SS) and Nicotiana tomentosiformis (2n = 2x = 24, TT). Genomic in situ hybridization (GISH), using the genomic DNA from these ancestral species as probes, revealed the chromosomal origins (S or T) and the occurrence of intergenomic translocations in N. tabacum. Fluorescence in situ hybridization (FISH) was also used to distinguish between chromosomes. However, the use of repetitive DNA sequences as probes for FISH analysis is limited by an inability to identify all chromosomes. In addition to this limitation, the occurrence of chromosomal tertiary constrictions can easily lead to the misclassification of chromosomes. To overcome these issues, immunostaining with anti-N. tabacum centromere-specific histone H3 antibody was carried out to determine the centromere position of each chromosome, followed by FISH analysis with ten distinct repetitive DNA probes. This approach allowed us to identify 22 of the 24 chromosome pairs in N. tabacum and revealed novel intergenomic chromosome rearrangements and B-chromosome-like minichromosomes. Hence, the combination of immunostaining with FISH and GISH is critical to accurately karyotype tobacco.

Reference genomes and transcriptomes of Nicotiana sylvestris and Nicotiana tomentosiformis

BACKGROUND

Nicotiana sylvestris and Nicotiana tomentosiformis are members of the Solanaceae family that includes tomato, potato, eggplant and pepper. These two Nicotiana species originate from South America and exhibit different alkaloid and diterpenoid production. N. sylvestris is cultivated largely as an ornamental plant and it has been used as a diploid model system for studies of terpenoid production, plastid engineering, and resistance to biotic and abiotic stress. N. sylvestris and N. tomentosiformis are considered to be modern descendants of the maternal and paternal donors that formed Nicotiana tabacum about 200,000 years ago through interspecific hybridization. Here we report the first genome-wide analysis of these two Nicotiana species.

RESULTS

Draft genomes of N. sylvestris and N. tomentosiformis were assembled to 82.9% and 71.6% of their expected size respectively, with N50 sizes of about 80 kb. The repeat content was 72-75%, with a higher proportion of retrotransposons and copia-like long terminal repeats in N. tomentosiformis. The transcriptome assemblies showed that 44,000-53,000 transcripts were expressed in the roots, leaves or flowers. The key genes involved in terpenoid metabolism, alkaloid metabolism and heavy metal transport showed differential expression in the leaves, roots and flowers of N. sylvestris and N. tomentosiformis.

CONCLUSIONS

The reference genomes of N. sylvestris and N. tomentosiformis represent a significant contribution to the SOL100 initiative because, as members of the Nicotiana genus of Solanaceae, they strengthen the value of the already existing resources by providing additional comparative information, thereby helping to improve our understanding of plant metabolism and evolution.

Deciphering the complex leaf transcriptome of the allotetraploid species Nicotiana tabacum: a phylogenomic perspective

BACKGROUND

Polyploidization is an important mechanism in plant evolution. By analyzing the leaf transcriptomes taken from the allotetraploid Nicotiana tabacum (tobacco) and parental genome donors, N. sylvesteris (S-Genome) and N. tomentosiformis (T-Genome), a phylogenomic approach was taken to map the fate of homeologous gene pairs in this plant.

RESULTS

A comparison between the genes present in the leaf transcriptomes of N. tabacum and modern day representatives of its progenitor species demonstrated that only 33% of assembled transcripts could be distinguished based on their sequences. A large majority of the genes (83.6% of the non parent distinguishable and 87.2% of the phylogenetic topology analyzed clusters) expressed above background level (more than 5 reads) showed similar overall expression levels. Homeologous sequences could be identified for 968 gene clusters, and 90% (6% of all genes) of the set maintained expression of only one of the tobacco homeologs. When both homeologs were expressed, only 15% (0.5% of the total) showed evidence of differential expression, providing limited evidence of subfunctionalization. Comparing the rate of synonymous nucleotide substitution (Ks) and non-synonymous nucleotide substitution (Kn) provided limited evidence for positive selection during the evolution of tobacco since the polyploidization event took place.

CONCLUSIONS

Polyploidization is a powerful mechanism for plant speciation that can occur during one generation; however millions of generations may be necessary for duplicate genes to acquire a new function. Analysis of the tobacco leaf transcriptome reveals that polyploidization, even in a young tetraploid such as tobacco, can lead to complex changes in gene expression. Gene loss and gene silencing, or subfunctionalization may explain why both homeologs are not expressed by the associated genes. With Whole Genome Duplication (WGD) events, polyploid genomes usually maintain a high percentage of gene duplicates. The data provided little evidence of preferential maintenance of gene expression from either the T- or S-genome. Additionally there was little evidence of neofunctionalization in Nicotiana tabacum suggesting it occurs at a low frequency in young polyploidy.

Three TERT genes in Nicotiana tabacum

Telomerase is essential for proper functioning of telomeres in eukaryotes. We cloned and characterised genes for the protein subunit of telomerase (TERT) in the allotetraploid Nicotiana tabacum (tobacco) and its diploid progenitor species Nicotiana sylvestris and Nicotiana tomentosiformis with the aim of determining if allopolyploidy (hybridisation and genome duplication) influences TERT activity and divergence. Two of the three sequence variants present in the tobacco genome (NtTERT-C/s and NtTERT-D) revealed similarity to two sequence variants found in N. sylvestris and another variant (NtTERT-C/t) was similar to TERT of N. tomentosiformis. Variants of N. sylvestris origin showed less similarity to each other (80.5 % in the genomic region; 90.1 % in the coding sequence) than that between the NtTERT-C/s and NtTERT-C/t variants (93.6 and 97.2 %, respectively). The NtTERT-D variant was truncated at the 5' end, and indels indicated that it was a pseudogene. All tobacco variants were transcribed and alternatively spliced sequences were detected. Analysis of gene arrangements uncovered a novel exon in the N-terminal domain of TERT variants, a feature that is likely to be commonly found in Solanaceae species. In addition, species-specific duplications were observed within exon 5. The putative function, copy number and evolutionary origin of these NtTERT sequence variants are discussed.

Isolation of centromeric-tandem repetitive DNA sequences by chromatin affinity purification using a HaloTag7-fused centromere-specific histone H3 in tobacco

The centromere is a multi-functional complex comprising centromeric DNA and a number of proteins. To isolate unidentified centromeric DNA sequences, centromere-specific histone H3 variants (CENH3) and chromatin immunoprecipitation (ChIP) have been utilized in some plant species. However, anti-CENH3 antibody for ChIP must be raised in each species because of its species specificity. Production of the antibodies is time-consuming and costly, and it is not easy to produce ChIP-grade antibodies. In this study, we applied a HaloTag7-based chromatin affinity purification system to isolate centromeric DNA sequences in tobacco. This system required no specific antibody, and made it possible to apply a highly stringent wash to remove contaminated DNA. As a result, we succeeded in isolating five tandem repetitive DNA sequences in addition to the centromeric retrotransposons that were previously identified by ChIP. Three of the tandem repeats were centromere-specific sequences located on different chromosomes. These results confirm the validity of the HaloTag7-based chromatin affinity purification system as an alternative method to ChIP for isolating unknown centromeric DNA sequences. The discovery of more than two chromosome-specific centromeric DNA sequences indicates the mosaic structure of tobacco centromeres.

Analysis of two abundant, highly related satellites in the allotetraploid Nicotiana arentsii using double-strand conformation polymorphism analysis and sequencing

• Allopolyploidy, a driving force in plant evolution, can induce rapid structural changes in parental subgenomes. Here, we examined the fate of homologous subtelomeric satellites in intrasection allotetraploid Nicotiana arentsii formed from N. undulata and N. wigandioides progenitors < 200,000 yr ago. • We cloned and sequenced a number of monomers from progenitors and the allotetraploid. Structural features of both cloned and genomic monomers were studied using double-strand conformation polymorphism analysis. • Two homologous satellites were isolated from N. undulata (called NUNSSP) and N. wigandioides (NWISSP). While the NUNSSP monomers were highly homogeneous in nucleotide sequences, the NWISSP monomers formed two separate clades. Likewise, the genomic NUNSSP monomers showed less DNA conformation heterogeneity than NWISSP monomers, with distinct conformations. While both satellites predominantly occupy subtelomeric positions, a fraction of the NWISSP repeats was found in an intercalary location, supporting the hypothesis that dispersion prevents the repeats becoming homogeneous. Sequence, structural and chromosomal features of the parental satellites were faithfully inherited by N. arentsii. • Our study revealed that intergenomic homogenization of subtelomeric satellite repeats does not occur in N. arentsii allotetraploid. We propose that the sequence and structural divergence of subtelomeric satellites may render allopolyploid chromosomes less vulnerable to intergenomic exchanges.

Organization and evolution of subtelomeric satellite repeats in the potato genome

Subtelomeric domains immediately adjacent to telomeres represent one of the most dynamic and rapidly evolving regions in eukaryotic genomes. A common feature associated with subtelomeric regions in different eukaryotes is the presence of long arrays of tandemly repeated satellite sequences. However, studies on molecular organization and evolution of subtelomeric repeats are rare. We isolated two subtelomeric repeats, CL14 and CL34, from potato (Solanum tuberosum). The CL14 and CL34 repeats are organized as independent long arrays, up to 1-3 Mb, of 182 bp and 339 bp monomers, respectively. The CL14 and CL34 repeat arrays are directly connected with the telomeric repeats at some chromosomal ends. The CL14 repeat was detected at the subtelomeric regions among highly diverged Solanum species, including tomato (Solanum lycopersicum). In contrast, CL34 was only found in potato and its closely related species. Interestingly, the CL34 repeat array was always proximal to the telomeres when both CL14 and CL34 were found at the same chromosomal end. In addition, the CL34 repeat family showed more sequence variability among monomers compared with the CL14 repeat family. We conclude that the CL34 repeat family emerged recently from the subtelomeric regions of potato chromosomes and is rapidly evolving. These results provide further evidence that subtelomeric domains are among the most dynamic regions in eukaryotic genomes.

Comparative genomics and repetitive sequence divergence in the species of diploid Nicotiana section Alatae

Combining phylogenetic reconstructions of species relationships with comparative genomic approaches is a powerful way to decipher evolutionary events associated with genome divergence. Here, we reconstruct the history of karyotype and tandem repeat evolution in species of diploid Nicotiana section Alatae. By analysis of plastid DNA, we resolved two clades with high bootstrap support, one containing N. alata, N. langsdorffii, N. forgetiana and N. bonariensis (called the n = 9 group) and another containing N. plumbaginifolia and N. longiflora (called the n = 10 group). Despite little plastid DNA sequence divergence, we observed, via fluorescent in situ hybridization, substantial chromosomal repatterning, including altered chromosome numbers, structure and distribution of repeats. Effort was focussed on 35S and 5S nuclear ribosomal DNA (rDNA) and the HRS60 satellite family of tandem repeats comprising the elements HRS60, NP3R and NP4R. We compared divergence of these repeats in diploids and polyploids of Nicotiana. There are dramatic shifts in the distribution of the satellite repeats and complete replacement of intergenic spacers (IGSs) of 35S rDNA associated with divergence of the species in section Alatae. We suggest that sequence homogenization has replaced HRS60 family repeats at sub-telomeric regions, but that this process may not occur, or occurs more slowly, when the repeats are found at intercalary locations. Sequence homogenization acts more rapidly (at least two orders of magnitude) on 35S rDNA than 5S rDNA and sub-telomeric satellite sequences. This rapid rate of divergence is analogous to that found in polyploid species, and is therefore, in plants, not only associated with polyploidy.

Sequence homogenization and chromosomal localization of VicTR-B satellites differ between closely related Vicia species

Satellite sequences of the VicTR-B family are specific for the genus Vicia (Leguminosae), but their abundance varies among the species, being the highest in Vicia sativa and Vicia grandiflora. In this study, we have sequenced multiple randomly cloned VicTR-B fragments from these two species and analyzed their sequence variability, periodicity, and chromosomal localization. We have found that V. sativa VicTR-B sequences are homogeneous with respect to their nucleotide sequences and periodicity (monomers of 38 bp), whereas V. grandiflora repeats are considerably more variable, occurring in at least four distinct sequence subfamilies. Although the periodicity of 38 bp was conserved in most of the V. grandiflora sequences, one of the subfamilies was composed of higher-order repeats of 186 bp, which originated from a pentamer of the basic repeated unit. Individual VicTR-B subfamilies were preferentially located in either intercalary or subtelomeric regions of chromosomes. Interestingly, two V. grandiflora subfamilies with the highest similarity to V. sativa VicTR-B sequences were located in intercalary heterochromatic bands, showing similar chromosomal distribution as the majority of VicTR-B repeats in V. sativa. The other two V. grandiflora subfamilies showing a considerable divergence from V. sativa sequences were found to be accumulated at subtelomeric regions of V. grandiflora chromosomes.

Telomeres in evolution and evolution of telomeres

This paper examines telomeres from an evolutionary perspective. In the monocot plant order Asparagales two evolutionary switch-points in telomere sequence are known. The first occurred when the Arabidopsis-type telomere was replaced by a telomere based on a repeat motif more typical of vertebrates. The replacement is associated with telomerase activity, but the telomerase has low fidelity and this may have implications for the binding of telomeric proteins. At the second evolutionary switch-point, the telomere and its mode of synthesis are replaced by an unknown mechanism. Elsewhere in plants (Sessia, Vestia, Cestrum) and in arthropods, the telomere "typical" of the group is lost. Probably many other groups with "unusual" telomeres will be found. We question whether telomerase is indeed the original end-maintenance system and point to other candidate processes involving t-loops, t-circles, rolling circle replication and recombination. Possible evolutionary outcomes arising from the loss of telomerase activity in alternative lengthening of telomere (ALT) systems are discussed. We propose that elongation of minisatellite repeats using recombination/replication processes initially substitutes for the loss of telomerase function. Then in more established ALT groups, subtelomeric satellite repeats may replace the telomeric minisatellite repeat whilst maintaining the recombination/replication mechanisms for telomere elongation. Thereafter a retrotransposition-based end-maintenance system may become established. The influence of changing sequence motifs on the properties of the telomere cap is discussed. The DNA and protein components of telomeres should be regarded--as with any other chromosome elements--as evolving and co-evolving over time and responding to changes in the genome and to environmental stresses. We describe how telomere dysfunction, resulting in end-to-end chromosome fusions, can have a profound effect on chromosome evolution and perhaps even speciation.

Preferential elimination of repeated DNA sequences from the paternal, Nicotiana tomentosiformis genome donor of a synthetic, allotetraploid tobacco

Nicotiana tabacum (tobacco, 2n = 4x = 48) is a natural allotetraploid combining two ancestral genomes closely related to modern Nicotiana sylvestris and Nicotiana tomentosiformis. Here we examine the immediate consequences of allopolyploidy on genome evolution using 20 S4-generation plants derived from a single synthetic, S0 plant made by Burk in 1973 (Th37). Using molecular and cytogenetic methods we analysed 14 middle and highly repetitive sequences that together total approximately 4% of the genome. Two repeats related to endogenous geminiviruses (GRD5) and pararetroviruses (NtoEPRV), and two classes of satellite repeats (NTRS, A1/A2) were partially or completely eliminated at variable frequency (25-60%). These sequences are all from the N. tomentosiformis parent. Genomic in situ hybridization revealed additivity in chromosome numbers in two plants (2n = 48), while a third was aneuploid for an N. tomentosiformis-origin chromosome (2n = 49). Two plants had homozygous translocations between chromosomes of the S- and T-genomes. * The data demonstrate that genetic changes in synthetic tobacco were fast, targeted to the paternal N. tomentosiformis-donated genome, and some of the changes showed concordance with changes that presumably occurred during evolution of natural tobacco.

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.

Dynamic Changes in the Distribution of a Satellite Homologous to Intergenic 26-18S rDNA Spacer in the Evolution of Nicotiana

An ∼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 S3 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.

Molecular characterization of a family of tandemly repeated DNA sequences, TR-1, in heterochromatic knobs of maize and its relatives

Two families of tandem repeats, 180-bp and TR-1, have been found in the knobs of maize. In this study, we isolated 59 clones belonging to the TR-1 family from maize and teosinte. Southern hybridization and sequence analysis revealed that members of this family are composed of three basic sequences, A (67 bp); B (184 bp) or its variants B' (184 bp), 2/3B (115 bp), 2/3B' (115 bp); and C (108 bp), which are arranged in various combinations to produce repeat units that are multiples of approximately 180 bp. The molecular structure of TR-1 elements suggests that: (1) the B component may evolve from the 180-bp knob repeat as a result of mutations during evolution; (2) B' may originate from B through lateral amplification accompanied by base-pair changes; (3) C plus A may be a single sequence that is added to B and B', probably via nonhomologous recombination; and (4) 69 bp at the 3' end of B or B', and the entire sequence of C can be removed from the elements by an unknown mechanism. Sequence comparisons showed partial homologies between TR-1 elements and two centromeric sequences (B repeats) of the supernumerary B chromosome. This result, together with the finding of other investigators that the B repeat is also fragmentarily homologous to the 180-bp repeat, suggests that the B repeat is derived from knob repeats in A chromosomes, which subsequently become structurally modified. Fluorescence in situ hybridization localized the B repeat to the B centromere and the 180-bp and TR-1 repeats to the proximal heterochromatin knob on the B chromosome.

Molecular characterization and physical localization of highly repetitive DNA sequences from Brazilian Alstroemeria species

Highly repetitive DNA sequences were isolated from genomic DNA libraries of Alstroemeria psittacina and A. inodora. Among the repetitive sequences that were isolated, tandem repeats as well as dispersed repeats could be discerned. The tandem repeats belonged to a family of interlinked Sau3A subfragments with sizes varying from 68-127 bp, and constituted a larger HinfI repeat of approximately 400 bp. Southern hybridization showed a similar molecular organization of the tandem repeats in each of the Brazilian Alstroemeria species tested. None of the repeats hybridized with DNA from Chilean Alstroemeria species, which indicates that they are specific for the Brazilian species. In-situ localization studies revealed the tandem repeats to be localized in clusters on the chromosomes of A. inodora and A. psittacina: distal hybridization sites were found on chromosome arms 2PS, 6PL, 7PS, 7PL and 8PL, interstitial sites on chromosome arms 2PL, 3PL, 4PL and 5PL. The applicability of the tandem repeats for cytogenetic analysis of interspecific hybrids and their role in heterochromatin organization are discussed.

TAS49--a dispersed repetitive sequence isolated from subtelomeric regions of Nicotiana tomentosiformis chromosomes

We have isolated and characterized a new repetitive sequence, TAS49, from terminal restriction fragments of Nicotiana tomentosiformis genomic DNA by means of a modified vectorette approach. The TAS49 was found directly attached to telomeres of N. tabacum and one of its ancestors, N. tomentosiformis, and also at inner chromosome locations. No association with telomeres was detected neither in N. otophora nor in the second tobacco ancestor, N. sylvestris. PCR and Southern hybridization reveal similarities in the arrangement of TAS49 on the chromosomes of 9 species of the genus Nicotiana, implying its occurrence as a subunit of a conserved complex DNA repeat. TAS49 belongs to the family of dispersed repetitive sequences without features of transposons. The copy number of TAS49 varies widely in the genomes of 8 species analyzed being lowest in N. sylvestris, with 3300 copies per diploid genome. In N. tomentosiformis, TAS49 forms about 0.56% of the diploid genome, corresponding to 17400 copies. TAS49 units are about 460 bp long and show about 90% of mutual homology, but no significant homology to DNA sequences deposited in GenBank and EMBL. Although genomic clones of TAS49 contain an open reading frame encoding a proline-rich protein similar to plant extensins, no mRNA transcript was detected. TAS49 is extensively methylated at CpG and CpNpG sites and its chromatin forms nucleosomes phased with a 170 +/- 8 bp periodicity.

Comparative analysis of DNA methylation in tobacco heterochromatic sequences

Cytosine methylation levels and susceptibility to drug-induced hypomethylation have been studied in several Nicotiana tabacum (tobacco) DNA repetitive sequences. It has been shown using HapII, MspI, BamHI and Sau3AI methylation-sensitive restriction enzymes that the degree of 5'-mCmCG-3' methylation varied significantly between different repeats. There were almost saturation levels of 5-methylcytosine at the inner (3') cytosine position and variable degrees of methylation at the outer (5') cytosine at the enzyme recognition sites. The non-transcribed high copy satellite sequences (HRS60, GRS) displayed significant heterogeneity in methylation of their basic units while middle repetitive sequences (R8.1, GRD5, 5S rDNA) were more uniformly modified at both cytosine residues. Dihydroxypropyladenine (DHPA) treatment, which is thought to reduce DNA methyltransferase activity by increasing S-adenosylhomocysteine levels, resulted in extensive demethylation of the outer cytosine in all repeats, and the partial hypomethylation of cytosines at the inner positions in less densely methylated repeats such as HRS60 and GRS. The results suggest that hypomethylation of 5'-mCmCG-3' sites with DHPA is a gradual non-random process proceeding in the direction mCmCG-->CmCG-->CCG. The 18S-5.8S-25S rDNA was remarkably hypomethylated relative to the 5S rDNA at all restriction sites studied. Fluorescence in-situ hybridization showed that DNA decondensation within and between the 18S-5.8S-25S and 5S rDNA loci was variable in different nuclei. All nuclei had condensed and decondensed sequence. The chromatin of 18S-5.8S-25S rDNA was more readily digested with micrococcal nuclease than the 5S rDNA suggesting that the overall levels of decondensation were higher for 18S-5.8S-25S rDNA. Variable decondensation patterns within and between loci were also observed for GRS and HRS60. Cytosine methylation of the tobacco repeats is discussed with respect to transcription, overall levels of condensation and overall structure.

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.

Two tandemly repeated telomere-associated sequences in Nicotiana plumbaginifolia

Two tandemly repeated telomere-associated sequences, NP3R and NP4R, have been isolated from Nicotiana plumbaginifolia. The length of a repeating unit for NP3R and NP4R is 165 and 180 nucleotides respectively. The abundance of NP3R, NP4R and telomeric repeats is, respectively, 8.4 x 10(4), 6 x 10(3) and 1.5 x 10(6) copies per haploid genome of N. plumbaginifolia. Fluorescence in situ hybridization revealed that NP3R is located at the ends and/or in interstitial regions of all 10 chromosomes and NP4R on the terminal regions of three chromosomes in the haploid genome of N. plumbaginifolia. Sequence homology search revealed that not only are NP3R and NP4R homologous to HRS60 and GRS, respectively, two tandem repeats isolated from N. tabacum, but that NP3R and NP4R are also related to each other, suggesting that they originated from a common ancestral sequence. The role of these repeated sequences in chromosome healing is discussed based on the observation that two to three copies of a telomere-similar sequence were present in each repeating unit of NP3R and NP4R.

Evidence for a sequence-directed conformation periodicity in the genomic highly repetitive DNA detectable with single-strand-specific chemical probe potassium permanganate

A single-strand-specific chemical probe, potassium permanganate (KMnO4), was used to study the sequence-dependent conformation periodicity of tandem multicopy repetitive DNA sequences HRS60 and GRS (Nicotiana Species) at the level of single base pair and dinucleotide step. Local DNA structures, sensitive to KMnO4, revealed periodicity of 182 +/- 2 bp, equal to the length of repeat units. Permanganate-sensitive local structures were mapped to both DNA strands of genomic HRS60 sequences and were found to be linked to d(A)n tracts. These adenine tracts are located in the proximity of the intrinsically curved domains. Distamycin A increased reactivity of the DNA but decreased the specificity of DNA cleavage. Similar conformation periodicity has been detected also in the 'canrep' family of repeats (Brassica species). All studied repetitive sequences are predominantly located in the constitutive heterochromatin. We discuss the role of conformation periodicities in relation to a structural code for nucleosome phasing at tandem arrays of DNA repeats.

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.

T-DNA integrations in a new family of repetitive elements of Nicotiana tabacum

A number of T-DNA insertions in the genome of Nicotiana tabacum were characterized. One class of integrations was found to have occurred in a new family of highly repetitive sequences. Three genomic regions (ecoA, ecoB, and ecoC) were isolated, all of which contain basic units of 180 bp, organized in direct tandem repeats. Several of the 180-bp elements contain an EcoRI recognition site within the repeating unit and are therefore named "eco repeats." All members of this family are weakly homologous in sequence to a previously described class of repeat elements which contained a BamHI site (HRS60 repeat family), which suggests that both groups of sequences are of common evolutionary origin. The allotetraploid genome of N. tabacum is presumed to originate from the hybridization of two diploid genomes. The HRS60 elements previously described have been found exclusively in the genome of one of the ancestors, N. sylvestris, and in N. tabacum itself. Our DNA hybridization data suggest that the eco elements originate from the genome of the other ancestor, N. tomentosiformis. Whereas the HRS60 elements are transcriptionally silent, at least some eco elements appear to be transcribed.

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.