Next generation sequencing reveals genome downsizing in allotetraploid Nicotiana tabacum, predominantly through the elimination of paternally derived repetitive DNAs

We used next generation sequencing to characterize and compare the genomes of the recently derived allotetraploid, Nicotiana tabacum (<200,000 years old), with its diploid progenitors, Nicotiana sylvestris (maternal, S-genome donor), and Nicotiana tomentosiformis (paternal, T-genome donor). Analysis of 14,634 repetitive DNA sequences in the genomes of the progenitor species and N. tabacum reveal all major types of retroelements found in angiosperms (genome proportions range between 17-22.5% and 2.3-3.5% for Ty3-gypsy elements and Ty1-copia elements, respectively). The diploid N. sylvestris genome exhibits evidence of recent bursts of sequence amplification and/or homogenization, whereas the genome of N. tomentosiformis lacks this signature and has considerably fewer homogenous repeats. In the derived allotetraploid N. tabacum, there is evidence of genome downsizing and sequences loss across most repeat types. This is particularly evident amongst the Ty3-gypsy retroelements in which all families identified are underrepresented in N. tabacum, as is 35S ribosomal DNA. Analysis of all repetitive DNA sequences indicates the T-genome of N. tabacum has experienced greater sequence loss than the S-genome, revealing preferential loss of paternally derived repetitive DNAs at a genome-wide level. Thus, the three genomes of N. sylvestris, N. tomentosiformis, and N. tabacum have experienced different evolutionary trajectories, with genomes that are dynamic, stable, and downsized, respectively.     

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.     

Analysis of multiple copies of geminiviral DNA in the genome of four closely related Nicotiana species suggest a unique integration event

Previously, we discovered multiple direct repeats of geminivirus-related DNA (GRD) sequences clustered at a single chromosomal position in Nicotiana tabacum (tobacco). Here we show that, in addition to tobacco, multiple copies of these elements occur in the genomes of three related Nicotiana species, all in the section Tomentosae: N. tomentosiformis, N. tomentosa and N. kawakamii, but not in 9 other more distantly related Nicotiana species, nor in various other solanaceous and non-solanaecous plants. DNA sequence analysis of 18 GRD copies reveal 4 distinct, but highly related, sub-families: GRD5, GRD3 and GRD53 in tobacco; GRD5 in N. tomentosiformis and N. kawakamii; and GRD2 in N. tomentosa. In addition to novel sequences, all elements share significant but varying lengths of DNA sequence similarity with the geminiviral replication origin plus the adjacent rep gene. There is extended sequence similarity to REP protein at the deduced amino acid sequence level, including motifs associated with other rolling circle replication proteins. Our data suggest that all GRD elements descend from a unique geminiviral integration event, most likely in a common ancestor of these Tomentosae species.     

Elimination and rearrangement of parental rDNA in the allotetraploid Nicotiana tabacum.

Origin and rearrangement of ribosomal DNA repeats in natural allotetraploid Nicotiana tabacum are described. Comparative sequence analysis of the intergenic spacer (IGS) regions of Nicotiana tomentosiformis (the paternal diploid progenitor) and Nicotiana sylvestris (the maternal diploid progenitor) showed species-specific molecular features. These markers allowed us to trace the molecular evolution of parental rDNA in the allopolyploid genome of N. tabacum; at least the majority of tobacco rDNA repeats originated from N. tomentosiformis, which endured reconstruction of subrepeated regions in the IGS. We infer that after hybridization of the parental diploid species, rDNA with a longer IGS, donated by N. tomentosiformis, dominated over the rDNA with a shorter IGS from N. sylvestris; the latter was then eliminated from the allopolyploid genome. Thus, repeated sequences in allopolyploid genomes are targets for molecular rearrangement, demonstrating the dynamic nature of allopolyploid genomes.     

The origin and evolution of geminivirus-related DNA sequences in Nicotiana

A horizontal transmission of a geminiviral DNA sequence, into the germ line of an ancestral Nicotiana, gave rise to multiple repeats of geminivirus-related DNA, GRD, in the genome. We follow GRD evolution in Nicotiana tabacum (tobacco), an allotetraploid, and its diploid relatives, and show GRDs are derived from begomoviruses. GRDs occur in two families: the GRD5 family's ancestor integrated into the common ancestor of three diploid species, Nicotiana kawakamii, Nicotiana tomentosa and Nicotiana tomentosiformis, on homeologous group 4 chromosomes. The GRD3 family was acquired more recently on chromosome 2 in a lineage of N. tomentosiformis, the paternal ancestor of tobacco. Both GRD families include individual members that are methylated and diverged. Using relative rates of synonymous and nonsynonymous nucleotide substitutions, we tested for evidence of selection on GRD units and found none within the GRD3 and GRD5 families. However, the substitutions between GRD3 and GRD5 do show a significant excess of synonymous changes, suggesting purifying selection and hence a period of autonomous evolution between GRD3 and GRD5 integration. We observe in the GRD3 family, features of Helitrons, a major new class of putative rolling-circle replicating eukaryotic transposon, not found in the GRD5 family or geminiviruses. We speculate that the second integration event, resulting in the GRD3 family, involved a free-living geminivirus, a Helitron and perhaps also GRD5. Thus our data point towards recurrent dynamic interplay between geminivirus and plant DNA in evolution.     

Evolution and structure of 5S rDNA loci in allotetraploid Nicotiana tabacum and its putative parental species

Nicotiana tabacum (tobacco) is an allotetraploid derived from ancestors of the modern diploids, N. sylvestris and N. tomentosiformis. We identified and characterized two distinct families of 5S ribosomal DNA (rDNA) in N. tabacum; one family had an average 431 bp unit length and the other a 646 bp unit length. In the diploid species, N. sylvestris and N. tomentosiformis, the 5S rDNA unit lengths are 431 bp and 644 bp respectively. The non-coding spacer sequence of the short unit in tobacco had high sequence homology to the spacer of N. sylvestris5S rDNA, while the longer spacer of tobacco had high homology with the 5S spacer of N. tomentosiformis. This suggests that the two 5S families in tobacco have their origin in the diploid ancestors. The longer spacer sequence had a GC rich sub-region (called the T-genome sub-region) that was absent in the short spacer. Pulsed field gel analysis and fluorescent in situ hybridization to tobacco metaphase chromosomes showed that the two families of 5S rDNA units are spatially separate at two chromosomal loci, on chromosomes S8 (short family) and T8 (long family). The repeat copy number at each chromosomal locus showed heterogeneity between different tobacco cultivars, with a tendency for a decrease in the copy number of one family to be compensated by an increase in the copy number of the second family. Sequence analysis reveals there is as much diversity in 5S family units within the diploid species as there is within the T and S-genome 5S family units respectively, suggesting 5S diversification within each family had occurred before tobacco speciation. There is no evidence of interlocus homogenization of the two 5S families in tobacco. This is therefore substantially different to 18-26S rDNA where interlocus gene conversion has substantially influenced most sequences of S and T genome origin; possible reasons are discussed.     

Molecular cytogenetic analyses and phylogenetic studies in the Nicotiana section Tomentosae

Phylogenetic schemes based on changing DNA sequence have made a major impact on our understanding of evolutionary relationships and significantly built on knowledge gained by morphological and anatomical studies. Here we present another approach to phylogeny, using fluorescent in situ hybridisation. The phylogenetic scheme presented is likely to be robust since it is derived from the chromosomal distribution of ten repetitive sequences with different functions and evolutionary constraints [GRS, HRS60, NTRS, the Arabidopsis-type telomere repeat (TTTAGGG)_n, 18S-5.8S-26S ribosomal DNA (rDNA), 5S rDNA, and four classes of geminiviral-related DNA (GRD)]. The basic karyotypes of all the plant species investigated Nicotiana tomentosiformis, N. kawakamii, N. tomentosa, N. otophora, N. setchellii, N. glutinosa (all section Tomentosae), and N. tabacum (tobacco, section Genuinae) are similar (x=12) but the distribution of genic and non-genic repeats is quite variable, making the karyotypes distinct. We found sequence dispersal, and locus gain, amplification and loss, all within the regular framework of the basic genomic structure. We predict that the GRD classes of sequence integrated into an ancestral genome only once in the evolution of section Tomentosae and thereafter spread by vertical transmission and speciation into four species. Since GRD is similar to a transgenic construct that was inserted into the N. tabacum genome, its fate over evolutionary time is interesting in the context of the debate on genetically modified organisms and the escape of genes into the wild. Nicotiana tabacum is thought to be an allotetraploid between presumed progenitors of N. sylvestris (maternal, S-genome donor) and a member of section Tomentosae (T-genome donor). Of section Tomentosae, N. tomentosiformis has the most similar genome to the T genome of tobacco and is therefore the most likely paternal genome donor. It is known for N. tabacum that gene conversion has converted most 18S-5.8S-26S rDNA units of N. sylvestris origin into units of an N. tomentosiformis type. Clearly if such a phenomenon were widespread across the genome, genomic in situ hybridisation (GISH) to distinguish the S and T genomes would probably not work since conversion would tend to homogenise the genomes. The fact that GISH does work suggests a limited role for gene conversion in the evolution of N. tabacum. Received: 8 November 1999; in revised form: 23 February 2000 / Accepted: 1 March 2000     

Origin of Nicotiana tabacum detected by primary structure of fraction I protein.

Nicotiana tabacum is believed to have arisen after hybridization of Nicotiana sylvestris with a species in the Tomentosae section of the genus Nicotiana. Recent biochemical experiments have confirmed the conclusions from previous cytogenetic studies that N. sylvestris was the maternal parent and have indicated that Nicotiana tomentosiformis was the paternal parent. However, these studies did not take into account the possibility that a new species of Nicotiana, called K-12, discovered in South America in 1968, could also have been one of the parents. Fraction I proteins were purified from N. tabacum and its putative progenitors, and separated into large and small subunits. Chymotryptic peptides of each subunit were analyzed by ion exchange column chromatography with a gradient elution system. Among 38 resolved peaks of the large subunits, 2 peaks were found to be different among the putative species. Since only N. sylvestris showed an identical chromatogram with N. tabacum, N. sylvestris was concluded to be the maternal progenitor, as the genetic information for the large subunit of Fraction I protein was known to be inherited by the cytoplasmic mode. On the other hand, the small subunit of Fraction I protein is inherited by the Mendelian mode and therefore N. tabacum, an allopolyploid, could be expected to contain two types of small subunits, one derived from N. sylvestris and the other from a paternal progenitor. N. sylvestris lacks two of the 25 chymotryptic peptides of the small subunit of N. tabacum. Among 3 putative paternal progenitors, these two peaks appeared only in N. tomentosiformis, but not in Nicotiana otophora or K-12. Thus, N. tomentosiformis was concluded to be a paternal progenitor of N. tabacum. The conclusion was verified by comparing chymotryptic peptides of small subunits from three amphidiploids of N. sylvestris crossed with N. tomentosiformis, N. sylvestris crossed with N. otophora snd N. sylvestris crossed with K-12. The analytical results showed that only the progeny of N. sylvestris crossed with N. tomentosiformis contained the same small subunits as N. tabacium.     

A distinct endogenous pararetrovirus family in Nicotiana tomentosiformis, a diploid progenitor of polyploid tobacco

A distinct endogenous pararetrovirus (EPRV) family corresponding to a previously unknown virus has been identified in the genome of Nicotiana tomentosiformis, a diploid ancestor of allotetraploid tobacco (Nicotiana tabacum). The putative virus giving rise to N. tomentosiformis EPRVs (NtoEPRVs) is most similar to tobacco vein clearing virus, an episomal form of a normally silent EPRV family in Nicotiana glutinosa; it is also related to a putative virus giving rise to the NsEPRV family in Nicotiana sylvestris (the second diploid progenitor of tobacco) and in the N. sylvestris fraction of the tobacco genome. The copy number of NtoEPRVs is significantly higher in N. tomentosiformis than in tobacco. This suggests that after the polyploidization event, many copies were lost from the polyploid genome or were accumulated specifically in the diploid genome. By contrast, the copy number of NsEPRVs has remained constant in N. sylvestris and tobacco, indicating that changes have occurred preferentially in the NtoEPRV family during evolution of the three Nicotiana species. NtoEPRVs are often flanked by Gypsy retrotransposon-containing plant DNA. Although the mechanisms of NtoEPRV integration, accumulation, and/or elimination are unknown, these processes are possibly linked to retrotransposon activity.     

Distribution of the Tnt1 retrotransposon family in the amphidiploid tobacco (Nicotiana tabacum) and its wild Nicotiana relatives

Transposable elements can generate considerable genetic diversity. Here we examine the distribution of the Tnt1 retrotransposon family in representative species of the genus Nicotiana . We show that multiple Tnt1 insertions are found in all Nicotiana species. However, Tnt1 insertions are too polymorphic to reveal species relationships. This indicates that Tnt1 has amplified rapidly and independently after Nicotiana speciation. We compare patterns of Tnt1 insertion in allotetraploid tobacco ( N. tabacum ) with those in the diploid species that are most closely related to the progenitors of tobacco, N. sylvestris (S-genome donor) and N. tomentosiformis (T-genome donor). We found no evidence for Tnt1 insertion sites of N. otophora origin in tobacco. Nicotiana sylvestris has a higher Tnt1 content than N. tomentosiformis and the elements are distributed more uniformly across the genome. This is reflected in tobacco where there is a higher Tnt1 content in S-genome chromosomes. However, the total Tnt1 content of tobacco is not the sum of the two modern-day parental species. We also observed tobacco-specific Tnt1 insertions and an absence of tobacco Tnt1 insertion sites in the diploid relatives. These data indicate Tnt1 evolution subsequent to allopolyploidy. We explore the possibility that fast evolution of Tnt1 is associated with 'genomic-shock' arising out of interspecific hybridization and allopolyploidy.  © 2004 The Linnean Society of London, Biological Journal of the Linnean Society , 2004, 82 , 639–649.     

Endogenous pararetroviruses of allotetraploid Nicotiana tabacum and its diploid progenitors, N. sylvestris and N. tomentosiformis

Endogenous pararetroviruses (EPRVs) represent a new class of dispersed repetitive DNA in plants. The genomes of many Nicotiana species and other solanaceous plants are rich in EPRVs. Distinct EPRV families are present in N. sylvestris ( Ns ) and in N. tomentosiformis ( Nto ), the two diploid progenitors of allotetraploid N. tabacum . Nicotiana EPRVs represent an interesting type of repetitive sequence to analyse in polyploids because of their potential impact on plant fitness and the epigenetic architecture of plant genomes. The Ns EPRV family appears identical in N. sylvestris and N. tabacum , indicating little change has occurred in either species since polyploid formation. By contrast, the Nto EPRV family is larger in N. tomentosiformis than in N. tabacum , suggesting either preferential elimination from the polyploid genome or specific accumulation in the diploid genome following polyploidization. The lability of Nto EPRVs might be enhanced by a frequent association with gypsy retrotransposons. Although some EPRVs are probably benign, others are potentially pathogenic or, conversely, determinants of virus resistance. Normally quiescent EPRVs can be reactivated and cause symptoms of infection in hybrids of species that differ in their EPRV content. EPRVs that furnish immunity to the free virus exemplify the selective value of so-called 'junk' DNA. Variation in the abundance and distribution of EPRVs among related species can be useful in taxonomic and evolutionary studies.  © 2004 The Linnean Society of London, Biological Journal of the Linnean Society , 2004, 82 , 627–638.     

Origin of sequence heterogeneity of the small subunit of fraction 1 protein from Nicotiana tabacum

The partial amino acid sequences of the small subunit of Fraction 1 protein from N. sylvestris, N. tomentosiformis and N. tabacum were determined. The sequence of N. sylvestris is NH2. Gln-Val-Trp-Pro-Pro-Ile-Asn-----Tyr COOH. In the sequence up to the 7th amino acid and C-terminus, differences were only found at the 6th position in the three species, where N. sylvestris and N. tomentosiformis show Ile and Tyr, respectively. N. tabacum show both Ile and Tyr in almost equal amount at this position. These results confirmed a previous hypothesis that N. tabacum had been evolved through the hybridization of N. sylvestris and N. tomentosiformis.     

AFLP analysis of genetic polymorphism and evolutionary relationships among cultivated and wild Nicotiana species.

Amplified fragment length polymorphism (AFLP) analysis was used to determine the degree of intra- and inter-specific genetic variation in the genus Nicotiana. Forty-six lines of cultivated tobacco (Nicotiana tabacum L.) and seven wild Nicotiana species, including N. sylvestris, N. tomentosiformis, N. otophora, N. glutinosa, N. suaveolens, N. rustica, and N. longiflora, were analyzed, using at least eight different oligonucleotide primer combinations capable of detecting a minimum of 50 polymorphic bands per primer pair. The amount of genetic polymorphism present among cultivated tobacco lines (N. tabacum) was limited, as evidenced by the high degree of similarity in the AFLP profiles of cultivars collected worldwide. Six major clusters were found within cultivated tobacco that were primarily based upon geographic origin and manufacturing quality traits. A greater amount of genetic polymorphism exists among wild species of Nicotiana than among cultivated forms. Pairwise comparisons of the AFLP profiles of wild and cultivated Nicotiana species show that polymorphic bands present in N. tabacum can be found in at least one of three proposed wild progenitor species (i.e., N. sylvestris, N. tomentosiformis, and N. otophora). This observation provides additional support for these species contributing to the origin of N. tabacum.     

Six active phage-type RNA polymerase genes in Nicotiana tabacum

In higher plants, a small nuclear gene family encodes mitochondrial as well as chloroplast RNA polymerases (RNAP) homologous to the bacteriophage T7-enzyme. The Arabidopsis genome contains three such RpoT genes, while in monocotyledonous plants only two copies have been found. Analysis of Nicotiana tabacum, a natural allotetraploid, identified six different RpoT sequences. The study of the progenitor species of tobacco, N. sylvestris and N. tomentosiformis, uncovered that the sequences represent two orthologous sets each of three RpoT genes (RpoT1, RpoT2 and RpoT3). Interestingly, while the organelles are inherited exclusively from the N. sylvestris maternal parent, all six RpoT genes are expressed in N. tabacum. GFP-fusions of Nicotiana RpoT1 revealed mitochondrial targeting properties. Constructs containing the amino-terminus of RpoT2 were imported into mitochondria as well as into plastids. Thus, the dual-targeting feature, first described for Arabidopsis RpoT;2, appears to be conserved among eudicotyledonous plants. Tobacco RpoT3 is targeted to chloroplasts and the RNA is differentially expressed in plants lacking the plastid-encoded RNAP. Remarkably, translation of RpoT3 mRNA has to be initiated at a CUG codon to generate a functional plastid transit peptide. Thus, besides AGAMOUS in Arabidopsis, Nicotiana RpoT3 provides a second example for a non-viral plant mRNA that is exclusively translated from a non-AUG codon.     

Rapid evolution of parental rDNA in a synthetic tobacco allotetraploid line

Unidirectional gene conversion of rDNA units has occurred in the evolution of natural tobacco (Nicotiana tabacum). In this paper we report the use of the synthetic tobacco line Th37, 4n (N. sylvestris × N. tomentosiformis), to study early rDNA evolution associated with allopolyploidy. At least three classes of newly amplified rDNA unit variants were identified (17/20 plants). Their presence was often accompanied by near-complete elimination of N. tomentosiformis-donated rDNA units (15/20 plants). Novel rDNA units were of N. tomentosiformis-type and contained rearranged subrepeats in the intergenic spacer. The maternal N. sylvestris-derived units were unchanged, except for some alteration in the ratio of individual gene family members. A cytogenetic analysis revealed rDNA sites on N. sylvestris-derived chromosomes S10, S11, and S12 and N. tomentosiformis-derived chromosomes T3 and in some cases T4. An rDNA locus does not occur on N. tomentosiformis chromosome 4. The locus on chromosome T4 of some hybrids correlates with the occurrence of the novel units that probably amplified at the locus. Combined with an analysis of tobacco cultivars, the data indicate that an initial burst of rDNA evolution associated with allopolyploidy was followed by a slower process that led towards reduced complexity and a decreased number of rDNA variants.