Repetitive DNA content in the maize genome is uncoupled from population stratification at SNP loci.

BACKGROUND: Repetitive DNA is a major component of plant genomes and is thought to be a driver of evolutionary novelty. Describing variation in repeat content among individuals and between populations is key to elucidating the evolutionary significance of repetitive DNA. However, the cost of producing references genomes has limited large-scale intraspecific comparisons to a handful of model organisms where multiple reference genomes are available. RESULTS: We examine repeat content variation in the genomes of 94 elite inbred maize lines using graph-based repeat clustering, a reference-free and rapid assay of repeat content. We examine population structure using genome-wide repeat profiles, and demonstrate the stiff-stalk and non-stiff-stalk heterotic populations are homogenous with regard to global repeat content. In contrast, and similar to previously reported results, the same individuals show clear differentiation, and aggregate into two populations when examining population structure using genome-wide SNPs. Additionally, we develop a novel kmer based technique to examine the chromosomal distribution of repeat clusters in silico and show a cluster dependent association with gene density. CONCLUSION: Our results indicate global repeat content variation in the heterotic populations of maize has not diverged, and is uncoupled from population stratification at SNP loci. We show that repeat families exhibit divergent patterns with regard to chromosomal distribution, some repeat clusters accumulate in regions of high gene density, whereas others aggregate in regions of low gene density.

Gene Fractionation and Function in the Ancient Subgenomes of Maize.

The maize genome experienced an ancient whole genome duplication ∼10 MYA and the duplicate subgenomes have since experienced reciprocal gene loss such that many genes have returned to single-copy status. This process has not affected the subgenomes equally; reduced gene expression in one of the subgenomes mitigates the consequences of mutations and gene deletions and is thought to drive higher rates of fractionation. Here, we use published data to show that, in accordance with predictions of this model, paralogs with greater expression contribute more to phenotypic variation compared with their lowly expressed counterparts. Furthermore, paralogous genes in the least-fractionated subgenome account for a greater degree of phenotypic diversity than those resident on the more-fractionated subgenome. Intriguingly, analysis of singleton genes reveals this difference persists even after fractionation is complete. Additionally, we show that the two subgenomes of maize may differ in their epigenetic profiles.

Persistence of subgenomes in paleopolyploid cotton after 60 my of evolution.

The importance of whole-genome multiplication (WGM) in plant evolution has long been recognized. In flowering plants, WGM is both ubiquitous and in many lineages cyclical, each round followed by substantial gene loss (fractionation). This process may be biased with respect to duplicated chromosomes, often with overexpression of genes in less fractionated relative to more fractionated regions. This bias is hypothesized to arise through downregulation of gene expression through silencing of local transposable elements (TEs). We assess differences in gene expression between duplicated regions of the paleopolyploid cotton genome and demonstrate that the rate of fractionation is negatively correlated with gene expression. We examine recent hypotheses regarding the source of fractionation bias and show that TE-mediated, positional downregulation is absent in the modern cotton genome, seemingly excluding this phenomenon as the primary driver of biased gene loss. Nevertheless, the paleo subgenomes of diploid cotton are still distinguishable with respect to TE content, targeting of 24-nt-small interfering RNAs and GC content, despite approximately 60 My of evolution. We propose that repeat content per se and differential recombination rates may drive biased fractionation following WGM. These data highlight the likely importance of ancient genomic fractionation biases in shaping modern crop genomes.

Analysis of the giant genomes of Fritillaria (Liliaceae) indicates that a lack of DNA removal characterizes extreme expansions in genome size.

Plants exhibit an extraordinary range of genome sizes, varying by > 2000-fold between the smallest and largest recorded values. In the absence of polyploidy, changes in the amount of repetitive DNA (transposable elements and tandem repeats) are primarily responsible for genome size differences between species. However, there is ongoing debate regarding the relative importance of amplification of repetitive DNA versus its deletion in governing genome size. Using data from 454 sequencing, we analysed the most repetitive fraction of some of the largest known genomes for diploid plant species, from members of Fritillaria. We revealed that genomic expansion has not resulted from the recent massive amplification of just a handful of repeat families, as shown in species with smaller genomes. Instead, the bulk of these immense genomes is composed of highly heterogeneous, relatively low-abundance repeat-derived DNA, supporting a scenario where amplified repeats continually accumulate due to infrequent DNA removal. Our results indicate that a lack of deletion and low turnover of repetitive DNA are major contributors to the evolution of extremely large genomes and show that their size cannot simply be accounted for by the activity of a small number of high-abundance repeat families.

Diploidization and genome size change in allopolyploids is associated with differential dynamics of low- and high-copy sequences.

Recent advances have highlighted the ubiquity of whole-genome duplication (polyploidy) in angiosperms, although subsequent genome size change and diploidization (returning to a diploid-like condition) are poorly understood. An excellent system to assess these processes is provided by Nicotiana section Repandae, which arose via allopolyploidy (approximately 5 million years ago) involving relatives of Nicotiana sylvestris and Nicotiana obtusifolia. Subsequent speciation in Repandae has resulted in allotetraploids with divergent genome sizes, including Nicotiana repanda and Nicotiana nudicaulis studied here, which have an estimated 23.6% genome expansion and 19.2% genome contraction from the early polyploid, respectively. Graph-based clustering of next-generation sequence data enabled assessment of the global genome composition of these allotetraploids and their diploid progenitors. Unexpectedly, in both allotetraploids, over 85% of sequence clusters (repetitive DNA families) had a lower abundance than predicted from their diploid relatives; a trend seen particularly in low-copy repeats. The loss of high-copy sequences predominantly accounts for the genome downsizing in N. nudicaulis. In contrast, N. repanda shows expansion of clusters already inherited in high copy number (mostly chromovirus-like Ty3/Gypsy retroelements and some low-complexity sequences), leading to much of the genome upsizing predicted. We suggest that the differential dynamics of low- and high-copy sequences reveal two genomic processes that occur subsequent to allopolyploidy. The loss of low-copy sequences, common to both allopolyploids, may reflect genome diploidization, a process that also involves loss of duplicate copies of genes and upstream regulators. In contrast, genome size divergence between allopolyploids is manifested through differential accumulation and/or deletion of high-copy-number sequences.

CenH3 evolution in diploids and polyploids of three angiosperm genera.

BACKGROUND: Centromeric DNA sequences alone are neither necessary nor sufficient for centromere specification. The centromere specific histone, CenH3, evolves rapidly in many species, perhaps as a coevolutionary response to rapidly evolving centromeric DNA. To gain insight into CenH3 evolution, we characterized patterns of nucleotide and protein diversity among diploids and allopolyploids within three diverse angiosperm genera, Brassica, Oryza, and Gossypium (cotton), with a focus on evidence for diversifying selection in the various domains of the CenH3 gene. In addition, we compare expression profiles and alternative splicing patterns for CenH3 in representatives of each genus. RESULTS: All three genera retain both duplicated CenH3 copies, while Brassica and Gossypium exhibit pronounced homoeologous expression level bias. Comparisons among genera reveal shared and unique aspects of CenH3 evolution, variable levels of diversifying selection in different CenH3 domains, and that alternative splicing contributes significantly to CenH3 diversity. CONCLUSIONS: Since the N terminus is subject to diversifying selection but the DNA binding domains do not appear to be, rapidly evolving centromere sequences are unlikely to be the primary driver of CenH3 sequence diversification. At present, the functional explanation for the diversity generated by both conventional protein evolution in the N terminal domain, as well as alternative splicing, remains unexplained.

Doubling down on genomes: polyploidy and crop plants.

Polyploidy, or whole genome multiplication, is ubiquitous among angiosperms. Many crop species are relatively recent allopolyploids, resulting from interspecific hybridization and polyploidy. Thus, an appreciation of the evolutionary consequences of (allo)polyploidy is central to our understanding of crop plant domestication, agricultural improvement, and the evolution of angiosperms in general. Indeed, many recent insights into plant biology have been gleaned from polyploid crops, including, but not limited to wheat, tobacco, sugarcane, apple, and cotton. A multitude of evolutionary processes affect polyploid genomes, including rapid and substantial genome reorganization, transgressive gene expression alterations, gene fractionation, gene conversion, genome downsizing, and sub- and neofunctionalization of duplicate genes. Often these genomic changes are accompanied by heterosis, robustness, and the improvement of crop yield, relative to closely related diploids. Historically, however, the genome-wide analysis of polyploid crops has lagged behind those of diploid crops and other model organisms. This lag is partly due to the difficulties in genome assembly, resulting from the genomic complexities induced by combining two or more evolutionarily diverged genomes into a single nucleus and by the significant size of polyploid genomes. In this review, we explore the role of polyploidy in angiosperm evolution, the domestication process and crop improvement. We focus on the potential of modern technologies, particularly next-generation sequencing, to inform us on the patterns and processes governing polyploid crop improvement and phenotypic change subsequent to domestication.

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

BACKGROUND: Tandemly arranged nuclear ribosomal DNA (rDNA), encoding 18S, 5.8S and 26S ribosomal RNA (rRNA), exhibit concerted evolution, a pattern thought to result from the homogenisation of rDNA arrays. However rDNA homogeneity at the single nucleotide polymorphism (SNP) level has not been detailed in organisms with more than a few hundred copies of the rDNA unit. Here we study rDNA complexity in species with arrays consisting of thousands of units. METHODS: We examined homogeneity of genic (18S) and non-coding internally transcribed spacer (ITS1) regions of rDNA using Roche 454 and/or Illumina platforms in four angiosperm species, Nicotiana sylvestris, N. tomentosiformis, N. otophora and N. kawakamii. We compared the data with Southern blot hybridisation revealing the structure of intergenic spacer (IGS) sequences and with the number and distribution of rDNA loci. RESULTS AND CONCLUSIONS: In all four species the intragenomic homogeneity of the 18S gene was high; a single ribotype makes up over 90% of the genes. However greater variation was observed in the ITS1 region, particularly in species with two or more rDNA loci, where >55% of rDNA units were a single ribotype, with the second most abundant variant accounted for >18% of units. IGS heterogeneity was high in all species. The increased number of ribotypes in ITS1 compared with 18S sequences may reflect rounds of incomplete homogenisation with strong selection for functional genic regions and relaxed selection on ITS1 variants. The relationship between the number of ITS1 ribotypes and the number of rDNA loci leads us to propose that rDNA evolution and complexity is influenced by locus number and/or amplification of orphaned rDNA units at new chromosomal locations.

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.

Next-generation sequencing and genome evolution in allopolyploids.

PREMISE OF THE STUDY: Hybridization and polyploidization (allopolyploidy) are ubiquitous in the evolution of plants, but tracing the origins and subsequent evolution of the constituent genomes of allopolyploids has been challenging. Genome doubling greatly complicates genetic analyses, and this has long hindered investigation in that most allopolyploid species are "nonmodel" organisms. However, recent advances in sequencing and genomics technologies now provide unprecedented opportunities to analyze numerous genetic markers in multiple individuals in any organism. METHODS: Here we review the application of next-generation sequencing technologies to the study of three aspects of allopolyploid genome evolution: duplicated gene loss and expression in two recently formed Tragopogon allopolyploids, intergenomic interactions and chromosomal evolution in Tragopogon miscellus, and repetitive DNA evolution in Nicotiana allopolyploids. KEY RESULTS: For the first time, we can explore on a genomic scale the evolutionary processes that are ongoing in natural allopolyploids and not be restricted to well-studied crops and genetic models. CONCLUSIONS: These approaches can be easily and inexpensively applied to many other plant species-making any evolutionarily provocative system a new "model" system.

Fall and rise of satellite repeats in allopolyploids of Nicotiana over c. 5 million years.

Allopolyploids represent natural experiments in which DNA sequences from different species are combined into a single nucleus and then coevolve, enabling us to follow the parental genomes, their interactions and evolution over time. Here, we examine the fate of satellite DNA over 5 million yr of divergence in plant genus Nicotiana (family Solanaceae). We isolated subtelomeric, tandemly repeated satellite DNA from Nicotiana diploid and allopolyploid species and analysed patterns of inheritance and divergence by sequence analysis, Southern blot hybridization and fluorescent in situ hybridization (FISH). We observed that parental satellite sequences redistribute around the genome in allopolyploids of Nicotiana section Polydicliae, formed c. 1 million yr ago (Mya), and that new satellite repeats evolved and amplified in section Repandae, which was formed c. 5 Mya. In some cases that process involved the complete replacement of parental satellite sequences. The rate of satellite repeat replacement is faster than theoretical predictions assuming the mechanism involved is unequal recombination and crossing-over. Instead we propose that this mechanism occurs with the deletion of large chromatin blocks and reamplification, perhaps via rolling circle replication.

Flow cytometry and GISH reveal mixed ploidy populations and Spartina nonaploids with genomes of S. alterniflora and S. maritima origin.

BACKGROUND: The genus Spartina exhibits extensive hybridization and includes classic examples of recent speciation by allopolyploidy. In the UK there are two hexaploid species, S. maritima and S. alterniflora, as well as the homoploid hybrid S. x townsendii (2n = 60) and a derived allododecaploid S. anglica (2n = 120, 122, 124); the latter two are considered to have originated in Hythe, southern England at the end of the 19th century. METHODS: Genomic in situ hybridization (GISH) and flow cytometry were used to characterize the genomic composition and distribution of these species and their ploidy levels at Eling Marchwood and Hythe, both near Southampton, southern England. KEY RESULTS: GISH identified approx. 60 chromosomes each of S. maritima and S. alterniflora origin in S. anglica and 62 chromosomes from S. alterniflora and 30 chromosomes from S. maritima in a nonaploid individual from Eling Marchwood, UK. GISH and flow cytometry also revealed that most (94 %) individuals examined at Hythe were hexaploid (the remaining two individuals (6 %) were dodedcaploid; n = 34), whereas hexaploid (approx. 36 % of plants), nonaploid (approx. 27 %) and dodecaploid (approx. 36 %) individuals were found at Eling Marchwood (n = 22). CONCLUSIONS: Nonaploid individuals indicate the potential for introgression between hexaploid and dodecaploid species, complicating the picture of polyploid-induced speciation within the genus. Despite the aggressive ecological habit of S. anglica, it has not out-competed S. x townsendii at Hythe (homoploid hybrids at a frequency of 94 %, n = 34), despite >100 years of coexistence. The success of GISH opens up the potential for future studies of polyploid-induced genome restructuring in this genus.

Independent Domestication of Two Old World Cotton Species.

Domesticated cotton species provide raw material for the majority of the world's textile industry. Two independent domestication events have been identified in allopolyploid cotton, one in Upland cotton (Gossypium hirsutum L.) and the other to Egyptian cotton (Gossypium barbadense L.). However, two diploid cotton species, Gossypium arboreum L. and Gossypium herbaceum L., have been cultivated for several millennia, but their status as independent domesticates has long been in question. Using genome resequencing data, we estimated the global abundance of various repetitive DNAs. We demonstrate that, despite negligible divergence in genome size, the two domesticated diploid cotton species contain different, but compensatory, repeat content and have thus experienced cryptic alterations in repeat abundance despite equivalence in genome size. Evidence of independent origin is bolstered by estimates of divergence times based on molecular evolutionary analysis of f7,000 orthologous genes, for which synonymous substitution rates suggest that G. arboreum and G. herbaceum last shared a common ancestor approximately 0.4-2.5 Ma. These data are incompatible with a shared domestication history during the emergence of agriculture and lead to the conclusion that G. arboreum and G. herbaceum were each domesticated independently.

A Cluster of Recently Inserted Transposable Elements Associated with siRNAs in Gossypium raimondii.

Stabilization of transposable element (TE) copy number involves the biosynthesis of short silencing RNAs (siRNAs) and further initialization of siRNA-mediated TE silencing. To gain insight into the relationship between the biosynthesis of siRNAs and their source TEs, we examined the co-evolutionary dynamics and expression of these two entities by characterizing the siRNA distribution across the genome of Gossypium raimondii Ulbr. We identified an unusual region at the 3' end of chromosome 1 with significantly enriched siRNA coverage. Analysis of the correlation pattern between uniquely mapped siRNAs and those mapping to multiple regions implicated active biogenesis of siRNAs from these potential young TEs. Furthermore, divergence estimates of TEs within this region confirmed that the majority of TEs are young. Active transcription of the source TEs and their positive correlation with expressed siRNAs indicates that sufficient expression of TEs may be necessary to generate siRNAs and maintain the silenced state of recently transposed TEs.

Ancient gene duplicates in Gossypium (cotton) exhibit near-complete expression divergence.

Whole genome duplication (WGD) is widespread in flowering plants and is a driving force in angiosperm diversification. The redundancy introduced by WGD allows the evolution of novel gene interactions and functions, although the patterns and processes of diversification are poorly understood. We identified ∼ 2,000 pairs of paralogous genes in Gossypium raimondii (cotton) resulting from an approximately 60 My old 5- to 6-fold ploidy increase. Gene expression analyses revealed that, in G. raimondii, 99.4% of the gene pairs exhibit differential expression in at least one of the three tissues (petal, leaf, and seed), with 93% to 94% exhibiting differential expression on a per-tissue basis. For 1,666 (85%) pairs, differential expression was observed in all tissues. These observations were mirrored in a time series of G. raimondii seed, and separately in leaf, petal, and seed of G. arboreum, indicating expression level diversification before species divergence. A generalized linear model revealed 92.4% of the paralog pairs exhibited expression divergence, with most exhibiting significant gene and tissue interactions indicating complementary expression patterns in different tissues. These data indicate massive, near-complete expression level neo- and/or subfunctionalization among ancient gene duplicates, suggesting these processes are essential in their maintenance over ∼ 60 Ma.

Independent, rapid and targeted loss of highly repetitive DNA in natural and synthetic allopolyploids of Nicotiana tabacum.

Allopolyploidy (interspecific hybridisation and polyploidy) has played a significant role in the evolutionary history of angiosperms and can result in genomic, epigenetic and transcriptomic perturbations. We examine the immediate effects of allopolyploidy on repetitive DNA by comparing the genomes of synthetic and natural Nicotiana tabacum with diploid progenitors N. tomentosiformis (paternal progenitor) and N. sylvestris (maternal progenitor). Using next generation sequencing, a recently developed graph-based repeat identification pipeline, Southern blot and fluorescence in situ hybridisation (FISH) we characterise two highly repetitive DNA sequences (NicCL3 and NicCL7/30). Analysis of two independent high-throughput DNA sequencing datasets indicates NicCL3 forms 1.6-1.9% of the genome in N. tomentosiformis, sequences that occur in multiple, discontinuous tandem arrays scattered over several chromosomes. Abundance estimates, based on sequencing depth, indicate NicCL3 is almost absent in N. sylvestris and has been dramatically reduced in copy number in the allopolyploid N. tabacum. Surprisingly elimination of NicCL3 is repeated in some synthetic lines of N. tabacum in their forth generation. The retroelement NicCL7/30, which occurs interspersed with NicCL3, is also under-represented but to a much lesser degree, revealing targeted elimination of the latter. Analysis of paired-end sequencing data indicates the tandem component of NicCL3 has been preferentially removed in natural N. tabacum, increasing the proportion of the dispersed component. This occurs across multiple blocks of discontinuous repeats and based on the distribution of nucleotide similarity among NicCL3 units, was concurrent with rounds of sequence homogenisation.