Mobilization of retrotransposons in synthetic allotetraploid tobacco

The genome and population genomics of allopolyploid Coffea arabica reveal the diversification history of modern coffee cultivars

Coffea arabica, an allotetraploid hybrid of Coffea eugenioides and Coffea canephora, is the source of approximately 60% of coffee products worldwide, and its cultivated accessions have undergone several population bottlenecks. We present chromosome-level assemblies of a di-haploid C. arabica accession and modern representatives of its diploid progenitors, C. eugenioides and C. canephora. The three species exhibit largely conserved genome structures between diploid parents and descendant subgenomes, with no obvious global subgenome dominance. We find evidence for a founding polyploidy event 350,000-610,000 years ago, followed by several pre-domestication bottlenecks, resulting in narrow genetic variation. A split between wild accessions and cultivar progenitors occurred ~30.5 thousand years ago, followed by a period of migration between the two populations. Analysis of modern varieties, including lines historically introgressed with C. canephora, highlights their breeding histories and loci that may contribute to pathogen resistance, laying the groundwork for future genomics-based breeding of C. arabica.

Transposable elements: multifunctional players in the plant genome

Transposable elements (TEs) are indispensable components of eukaryotic genomes that play diverse roles in gene regulation, recombination, and environmental adaptation. Their ability to mobilize within the genome leads to gene expression and DNA structure changes. TEs serve as valuable markers for genetic and evolutionary studies and facilitate genetic mapping and phylogenetic analysis. They also provide insight into how organisms adapt to a changing environment by promoting gene rearrangements that lead to new gene combinations. These repetitive sequences significantly impact genome structure, function and evolution. This review takes a comprehensive look at TEs and their applications in biotechnology, particularly in the context of plant biology, where they are now considered "genomic gold" due to their extensive functionalities. The article addresses various aspects of TEs in plant development, including their structure, epigenetic regulation, evolutionary patterns, and their use in gene editing and plant molecular markers. The goal is to systematically understand TEs and shed light on their diverse roles in plant biology.

The Genomic Shock Hypothesis: Genetic and Epigenetic Alterations of Transposable Elements after Interspecific Hybridization in Plants

Transposable elements (TEs) are major components of plant genomes with the ability to change their position in the genome or to create new copies of themselves in other positions in the genome. These can cause gene disruption and large-scale genomic alterations, including inversions, deletions, and duplications. Host organisms have evolved a set of mechanisms to suppress TE activity and counter the threat that they pose to genome integrity. These includes the epigenetic silencing of TEs mediated by a process of RNA-directed DNA methylation (RdDM). In most cases, the silencing machinery is very efficient for the vast majority of TEs. However, there are specific circumstances in which TEs can evade such silencing mechanisms, for example, a variety of biotic and abiotic stresses or in vitro culture. Hybridization is also proposed as an inductor of TE proliferation. In fact, the discoverer of the transposons, Barbara McClintock, first hypothesized that interspecific hybridization provides a "genomic shock" that inhibits the TE control mechanisms leading to the mobilization of TEs. However, the studies carried out on this topic have yielded diverse results, showing in some cases a total absence of mobilization or being limited to only some TE families. Here, we review the current knowledge about the impact of interspecific hybridization on TEs in plants and the possible implications of changes in the epigenetic mechanisms.

TEMM: A Curated Data Resource for Transposon Element-Based Molecular Markers in Plants

Transposon elements (TEs) are mobile genetic elements that can insert themselves into new locations and modify the plant genome. In recent years, they have been used as molecular markers in plant breeding programs. TE-based molecular markers (TE-markers) are divided into two categories depending on the transcription mechanism of the TEs. The first category is retrotransposon-based molecular markers, which include RBIP, IRAP, REMAP, and iPBS. The second group is DNA-based-TE-markers, which include MITE, TE-junction, and CACTA TE-markers. These markers are a good tool for studying genetic diversity and can provide information on plants' phylogenetic and evolutionary history. They can help improve breeding programs to increase agronomic traits and develop new varieties. Overall, TE-markers play an important role in plant genetics and plant breeding and contribute to a better understanding of plant biology. Here, we present TEMM, a curated data resource for TE-markers in plants. Relevant research articles were screened to collect primer sequences and related information. Only articles containing primer sequences are added to the present data resource. TEMM contains 784 primers with their associated PCR reaction programs and their applications in various crops. These include 203 IPBS, 191 RBIP, 140 IRAP, 78 TE-junction, 76 IRAPS, 47 RBIP-IRAP, 16 IRAP-REMAP, 12 REMAP, 12 REMA-IRAP, 6 REMA, and 3 ISBP primers. The data resource is freely available at https://bioinformatics.um6p.ma/TEMM .

Measuring the Invisible: The Sequences Causal of Genome Size Differences in Eyebrights (Euphrasia) Revealed by k-mers

Genome size variation within plant taxa is due to presence/absence variation, which may affect low-copy sequences or genomic repeats of various frequency classes. However, identifying the sequences underpinning genome size variation is challenging because genome assemblies commonly contain collapsed representations of repetitive sequences and because genome skimming studies by design miss low-copy number sequences. Here, we take a novel approach based on k-mers, short sub-sequences of equal length k, generated from whole-genome sequencing data of diploid eyebrights (Euphrasia), a group of plants that have considerable genome size variation within a ploidy level. We compare k-mer inventories within and between closely related species, and quantify the contribution of different copy number classes to genome size differences. We further match high-copy number k-mers to specific repeat types as retrieved from the RepeatExplorer2 pipeline. We find genome size differences of up to 230Mbp, equivalent to more than 20% genome size variation. The largest contributions to these differences come from rDNA sequences, a 145-nt genomic satellite and a repeat associated with an Angela transposable element. We also find size differences in the low-copy number class (copy number ≤ 10×) of up to 27 Mbp, possibly indicating differences in gene space between our samples. We demonstrate that it is possible to pinpoint the sequences causing genome size variation within species without the use of a reference genome. Such sequences can serve as targets for future cytogenetic studies. We also show that studies of genome size variation should go beyond repeats if they aim to characterise the full range of genomic variants. To allow future work with other taxonomic groups, we share our k-mer analysis pipeline, which is straightforward to run, relying largely on standard GNU command line tools.

Global Patterns of Subgenome Evolution in Organelle-Targeted Genes of Six Allotetraploid Angiosperms

Whole-genome duplications (WGDs) are a prominent process of diversification in eukaryotes. The genetic and evolutionary forces that WGD imposes on cytoplasmic genomes are not well understood, despite the central role that cytonuclear interactions play in eukaryotic function and fitness. Cellular respiration and photosynthesis depend on successful interaction between the 3,000+ nuclear-encoded proteins destined for the mitochondria or plastids and the gene products of cytoplasmic genomes in multi-subunit complexes such as OXPHOS, organellar ribosomes, Photosystems I and II, and Rubisco. Allopolyploids are thus faced with the critical task of coordinating interactions between the nuclear and cytoplasmic genes that were inherited from different species. Because the cytoplasmic genomes share a more recent history of common descent with the maternal nuclear subgenome than the paternal subgenome, evolutionary "mismatches" between the paternal subgenome and the cytoplasmic genomes in allopolyploids might lead to the accelerated rates of evolution in the paternal homoeologs of allopolyploids, either through relaxed purifying selection or strong directional selection to rectify these mismatches. We report evidence from six independently formed allotetraploids that the subgenomes exhibit unequal rates of protein-sequence evolution, but we found no evidence that cytonuclear incompatibilities result in altered evolutionary trajectories of the paternal homoeologs of organelle-targeted genes. The analyses of gene content revealed mixed evidence for whether the organelle-targeted genes are lost more rapidly than the non-organelle-targeted genes. Together, these global analyses provide insights into the complex evolutionary dynamics of allopolyploids, showing that the allopolyploid subgenomes have separate evolutionary trajectories despite sharing the same nucleus, generation time, and ecological context.

An Eruption of LTR Retrotransposons in the Autopolyploid Genomes of Chrysanthemum nankingense (Asteraceae)

Whole genome duplication, associated with the induction of widespread genetic changes, has played an important role in the evolution of many plant taxa. All extant angiosperm species have undergone at least one polyploidization event, forming either an auto- or allopolyploid organism. Compared with allopolyploidization, however, few studies have examined autopolyploidization, and few studies have focused on the response of genetic changes to autopolyploidy. In the present study, newly synthesized C. nankingense autotetraploids (Asteraceae) were employed to characterize the genome shock following autopolyploidization. Available evidence suggested that the genetic changes primarily involved the loss of old fragments and the gain of novel fragments, and some novel sequences were potential long terminal repeat (LTR) retrotransposons. As Ty1-copia and Ty3-gypsy elements represent the two main superfamilies of LTR retrotransposons, the dynamics of Ty1-copia and Ty3-gypsy were evaluated using RT-PCR, transcriptome sequencing, and LTR retrotransposon-based molecular marker techniques. Additionally, fluorescence in situ hybridization(FISH)results suggest that autopolyploidization might also be accompanied by perturbations of LTR retrotransposons, and emergence retrotransposon insertions might show more rapid divergence, resulting in diploid-like behaviour, potentially accelerating the evolutionary process among progenies. Our results strongly suggest a need to expand the current evolutionary framework to include a genetic dimension when seeking to understand genomic shock following autopolyploidization in Asteraceae.

Genome-Wide Identification and Characterization of the RCI2 Gene Family in Allotetraploid Brassica napus Compared with Its Diploid Progenitors

Brassica napus and its diploid progenitors (B. rapa and B. oleracea) are suitable for studying the problems associated with polyploidization. As an important anti-stress protein, RCI2 proteins widely exist in various tissues of plants, and are crucial to plant growth, development, and stress response. In this study, the RCI2 gene family was comprehensively identified and analyzed, and 9, 9, and 24 RCI2 genes were identified in B. rapa, B. oleracea, and B. napus, respectively. Phylogenetic analysis showed that all of the identified RCI2 genes were divided into two groups, and further divided into three subgroups. Ka/Ks analysis showed that most of the identified RCI2 genes underwent a purifying selection after the duplication events. Moreover, gene structure analysis showed that the structure of RCI2 genes is largely conserved during polyploidization. The promoters of the RCI2 genes in B. napus contained more cis-acting elements, which were mainly involved in plant development and growth, plant hormone response, and stress responses. Thus, B. napus might have potential advantages in some biological aspects. In addition, the changes of RCI2 genes during polyploidization were also discussed from the aspects of gene number, gene structure, gene relative location, and gene expression, which can provide reference for future polyploidization analysis.

Genomic and Meiotic Changes Accompanying Polyploidization

Hybridization and polyploidy have been considered as significant evolutionary forces in adaptation and speciation, especially among plants. Interspecific gene flow generates novel genetic variants adaptable to different environments, but it is also a gene introgression mechanism in crops to increase their agronomical yield. An estimate of 9% of interspecific hybridization has been reported although the frequency varies among taxa. Homoploid hybrid speciation is rare compared to allopolyploidy. Chromosome doubling after hybridization is the result of cellular defects produced mainly during meiosis. Unreduced gametes, which are formed at an average frequency of 2.52% across species, are the result of altered spindle organization or orientation, disturbed kinetochore functioning, abnormal cytokinesis, or loss of any meiotic division. Meiotic changes and their genetic basis, leading to the cytological diploidization of allopolyploids, are just beginning to be understood especially in wheat. However, the nature and mode of action of homoeologous recombination suppressor genes are poorly understood in other allopolyploids. The merger of two independent genomes causes a deep modification of their architecture, gene expression, and molecular interactions leading to the phenotype. We provide an overview of genomic changes and transcriptomic modifications that particularly occur at the early stages of allopolyploid formation.

Genome downsizing after polyploidy: mechanisms, rates and selection pressures

An analysis of over 10 000 plant genome sizes (GSs) indicates that most species have smaller genomes than expected given the incidence of polyploidy in their ancestries, suggesting selection for genome downsizing. However, comparing ancestral GS with the incidence of ancestral polyploidy suggests that the rate of DNA loss following polyploidy is likely to have been very low (4-70 Mb/million years, 4-482 bp/generation). This poses a problem. How might such small DNA losses be visible to selection, overcome the power of genetic drift and drive genome downsizing? Here we explore that problem, focussing on the role that double-strand break (DSB) repair pathways (non-homologous end joining and homologous recombination) may have played. We also explore two hypotheses that could explain how selection might favour genome downsizing following polyploidy: to reduce (i) nitrogen (N) and phosphate (P) costs associated with nucleic acid synthesis in the nucleus and the transcriptome and (ii) the impact of scaling effects of GS on cell size, which influences CO₂ uptake and water loss. We explore the hypothesis that losses of DNA must be fastest in early polyploid generations. Alternatively, if DNA loss is a more continuous process over evolutionary time, then we propose it is a byproduct of selection elsewhere, such as limiting the damaging activity of repetitive DNA. If so, then the impact of GS on photosynthesis, water use efficiency and/or nutrient costs at the nucleus level may be emergent properties, which have advantages, but not ones that could have been selected for over generational timescales.

Comparative phylogenetics of repetitive elements in a diverse order of flowering plants (Brassicales)

Genome sizes of plants have long piqued the interest of researchers due to the vast differences among organisms. However, the mechanisms that drive size differences have yet to be fully understood. Two important contributing factors to genome size are expansions of repetitive elements, such as transposable elements (TEs), and whole-genome duplications (WGD). Although studies have found correlations between genome size and both TE abundance and polyploidy, these studies typically test for these patterns within a genus or species. The plant order Brassicales provides an excellent system to further test if genome size evolution patterns are consistent across larger time scales, as there are numerous WGDs. This order is also home to one of the smallest plant genomes, Arabidopsis thaliana-chosen as the model plant system for this reason-as well as to species with very large genomes. With new methods that allow for TE characterization from low-coverage genome shotgun data and 71 taxa across the Brassicales, we confirm the correlation between genome size and TE content, however, we are unable to reconstruct phylogenetic relationships and do not detect any shift in TE abundance associated with WGD.

Subgenome evolution in allotetraploid plants

Polyploidization is a well-known speciation and adaptation mechanism. Traces of former polyploidization events were discovered within many genomes, and especially in plants. Allopolyploidization by interspecific hybridization between two species is common. Among hybrid plants, many are domesticated species of agricultural interest and many of their genomes and of their presumptive parents have been sequenced. Hybrid genomes remain challenging to analyse because of the presence of multiple subgenomes. The genomes of hybrids often undergo rearrangement and degradation over time. Based on 10 hybrid plant genomes from six different genera, with hybridization dating from 10,000 to 5 million years ago, we assessed subgenome degradation, subgenomic intermixing and biased subgenome fractionation. The restructuring of hybrid genomes does not proceed proportionally with the age of the hybrid. The oldest hybrids in our data set display completely different fates: whereas the subgenomes of the tobacco plant Nicotiana benthamiana are in an advanced stage of degradation, the subgenomes of quinoa (Chenopodium quinoa) are exceptionally well conserved by structure and sequence. We observed statistically significant biased subgenome fractionation in seven out of 10 hybrids, which had different ages and subgenomic intermixing levels. Hence, we conclude that no correlation exists between biased fractionation and subgenome intermixing. Lastly, domestication may encourage or hinder subgenome intermixing, depending on the evolutionary context. In summary, comparative analysis of hybrid genomes and their presumptive parents allowed us to determine commonalities and differences between their evolutionary fates. In order to facilitate the future analysis of further hybrid genomes, we automated the analysis steps within manticore, which is publicly available at https://github.com/MatteoSchiavinato/manticore.git.

Gene and Transposable Element Expression Evolution Following Recent and Past Polyploidy Events in Spartina (Poaceae)

Gene expression dynamics is a key component of polyploid evolution, varying in nature, intensity, and temporal scales, most particularly in allopolyploids, where two or more sub-genomes from differentiated parental species and different repeat contents are merged. Here, we investigated transcriptome evolution at different evolutionary time scales among tetraploid, hexaploid, and neododecaploid Spartina species (Poaceae, Chloridoideae) that successively diverged in the last 6-10 my, at the origin of differential phenotypic and ecological traits. Of particular interest are the recent (19th century) hybridizations between the two hexaploids Spartina alterniflora (2n = 6x = 62) and S. maritima (2n = 6x = 60) that resulted in two sterile F1 hybrids: Spartina × townsendii (2n = 6x = 62) in England and Spartina × neyrautii (2n = 6x = 62) in France. Whole genome duplication of S. × townsendii gave rise to the invasive neo-allododecaploid species Spartina anglica (2n = 12x = 124). New transcriptome assemblies and annotations for tetraploids and the enrichment of previously published reference transcriptomes for hexaploids and the allododecaploid allowed identifying 42,423 clusters of orthologs and distinguishing 21 transcribed transposable element (TE) lineages across the seven investigated Spartina species. In 4x and 6x mesopolyploids, gene and TE expression changes were consistent with phylogenetic relationships and divergence, revealing weak expression differences in the tetraploid sister species Spartina bakeri and Spartina versicolor (<2 my divergence time) compared to marked transcriptome divergence between the hexaploids S. alterniflora and S. maritima that diverged 2-4 mya. Differentially expressed genes were involved in glycolysis, post-transcriptional protein modifications, epidermis development, biosynthesis of carotenoids. Most detected TE lineages (except SINE elements) were found more expressed in hexaploids than in tetraploids, in line with their abundance in the corresponding genomes. Comparatively, an astonishing (52%) expression repatterning and deviation from parental additivity were observed following recent reticulate evolution (involving the F1 hybrids and the neo-allododecaploid S. anglica), with various patterns of biased homoeologous gene expression, including genes involved in epigenetic regulation. Downregulation of TEs was observed in both hybrids and accentuated in the neo-allopolyploid. Our results reinforce the view that allopolyploidy represents springboards to new regulatory patterns, offering to worldwide invasive species, such as S. anglica, the opportunity to colonize stressful and fluctuating environments on saltmarshes.

Genomic imprinted genes in reciprocal hybrid endosperm of Brassica napus

BACKGROUND

Genomic imprinting results in the expression of parent-of-origin-specific alleles in the offspring. Brassica napus is an oil crop with research values in polyploidization. Identification of imprinted genes in B. napus will enrich the knowledge of genomic imprinting in dicotyledon plants.

RESULTS

In this study, we performed reciprocal crosses between B. napus L. cultivars Yangyou 6 (Y6) and Zhongshuang 11 (ZS11) to collect endosperm at 20 and 25 days after pollination (DAP) for RNA-seq. In total, we identified 297 imprinted genes, including 283 maternal expressed genes (MEGs) and 14 paternal expressed genes (PEGs) according to the SNPs between Y6 and ZS11. Only 36 genes (35 MEGs and 1 PEG) were continuously imprinted in 20 and 25 DAP endosperm. We found 15, 2, 5, 3, 10, and 25 imprinted genes in this study were also imprinted in Arabidopsis, rice, castor bean, maize, B. rapa, and other B. napus lines, respectively. Only 26 imprinted genes were specifically expressed in endosperm, while other genes were also expressed in root, stem, leaf and flower bud of B. napus. A total of 109 imprinted genes were clustered on rapeseed chromosomes. We found the LTR/Copia transposable elements (TEs) were most enriched in both upstream and downstream of the imprinted genes, and the TEs enriched around imprinted genes were more than non-imprinted genes. Moreover, the expression of 5 AGLs and 6 pectin-related genes in hybrid endosperm were significantly changed comparing with that in parent endosperm.

CONCLUSION

This research provided a comprehensive identification of imprinted genes in B. napus, and enriched the gene imprinting in dicotyledon plants, which would be useful in further researches on how gene imprinting regulates seed development.

Evolutionary dynamics of transposable elements and satellite DNAs in polyploid Spartina species

Repeated sequences and polyploidy play a central role in plant genome dynamics. Here, we analyze the evolutionary dynamics of repeats in tetraploid and hexaploid Spartina species that diverged during the last 10 million years within the Chloridoideae, one of the poorest investigated grass lineages. From high-throughput genome sequencing, we annotated Spartina repeats and determined what sequence types account for the genome size variation among species. We examined whether differential genome size evolution correlated with ploidy levels and phylogenetic relationships. We also examined the tempo of repeat sequence dynamics associated with allopatric speciation over the last 3-6 million years between hexaploid species that diverged on the American and European Atlantic coasts and tetraploid species from North and South America. The tetraploid S. spartinae, whose phylogenetic placement has been debated, exhibits a similar repeat content as hexaploid species, suggesting common ancestry. Genome expansion or contraction resulting from repeat dynamics seems to be explained mostly by the contrasting divergence times between species, rather than by genome changes triggered by ploidy level change per se. One 370 bp satellite may be exhibiting 'meiotic drive' and driving chromosome evolution in S. alterniflora. Our results provide crucial insights for investigating the genetic and epigenetic consequences of such differential repeat dynamics on the ecology and distribution of the meso- and neopolyploid Spartina species.

Characterization and Dynamics of Repeatomes in Closely Related Species of Hieracium (Asteraceae) and Their Synthetic and Apomictic Hybrids

The repetitive content of the plant genome (repeatome) often represents its largest fraction and is frequently correlated with its size. Transposable elements (TEs), the main component of the repeatome, are an important driver in the genome diversification due to their fast-evolving nature. Hybridization and polyploidization events are hypothesized to induce massive bursts of TEs resulting, among other effects, in an increase of copy number and genome size. Little is known about the repeatome dynamics following hybridization and polyploidization in plants that reproduce by apomixis (asexual reproduction via seeds). To address this, we analyzed the repeatomes of two diploid parental species, Hieracium intybaceum and H. prenanthoides (sexual), their diploid F1 synthetic and their natural triploid hybrids (H. pallidiflorum and H. picroides, apomictic). Using low-coverage next-generation sequencing (NGS) and a graph-based clustering approach, we detected high overall similarity across all major repeatome categories between the parental species, despite their large phylogenetic distance. Medium and highly abundant repetitive elements comprise ∼70% of Hieracium genomes; most prevalent were Ty3/Gypsy chromovirus Tekay and Ty1/Copia Maximus-SIRE elements. No TE bursts were detected, neither in synthetic nor in natural hybrids, as TE abundance generally followed theoretical expectations based on parental genome dosage. Slight over- and under-representation of TE cluster abundances reflected individual differences in genome size. However, in comparative analyses, apomicts displayed an overabundance of pararetrovirus clusters not observed in synthetic hybrids. Substantial deviations were detected in rDNAs and satellite repeats, but these patterns were sample specific. rDNA and satellite repeats (three of them were newly developed as cytogenetic markers) were localized on chromosomes by fluorescence in situ hybridization (FISH). In a few cases, low-abundant repeats (5S rDNA and certain satellites) showed some discrepancy between NGS data and FISH results, which is due partly to the bias of low-coverage sequencing and partly to low amounts of the satellite repeats or their sequence divergence. Overall, satellite DNA (including rDNA) was markedly affected by hybridization, but independent of the ploidy or reproductive mode of the progeny, whereas bursts of TEs did not play an important role in the evolutionary history of Hieracium.

Genomics of Evolutionary Novelty in Hybrids and Polyploids

It has long been recognized that hybridization and polyploidy are prominent processes in plant evolution. Although classically recognized as significant in speciation and adaptation, recognition of the importance of interspecific gene flow has dramatically increased during the genomics era, concomitant with an unending flood of empirical examples, with or without genome doubling. Interspecific gene flow is thus increasingly thought to lead to evolutionary innovation and diversification, via adaptive introgression, homoploid hybrid speciation and allopolyploid speciation. Less well understood, however, are the suite of genetic and genomic mechanisms set in motion by the merger of differentiated genomes, and the temporal scale over which recombinational complexity mediated by gene flow might be expressed and exposed to natural selection. We focus on these issues here, considering the types of molecular genetic and genomic processes that might be set in motion by the saltational event of genome merger between two diverged species, either with or without genome doubling, and how these various processes can contribute to novel phenotypes. Genetic mechanisms include the infusion of new alleles and the genesis of novel structural variation including translocations and inversions, homoeologous exchanges, transposable element mobilization and novel insertional effects, presence-absence variation and copy number variation. Polyploidy generates massive transcriptomic and regulatory alteration, presumably set in motion by disrupted stoichiometries of regulatory factors, small RNAs and other genome interactions that cascade from single-gene expression change up through entire networks of transformed regulatory modules. We highlight both these novel combinatorial possibilities and the range of temporal scales over which such complexity might be generated, and thus exposed to natural selection and drift.

Parental origin of the allotetraploid tobacco Nicotiana benthamiana

Nicotiana section Suaveolentes is an almost all-Australian clade of allopolyploid tobacco species including the important plant model Nicotiana benthamiana. The homology relationships of this clade and its formation are not completely understood. To address this gap, we assessed phylogenies of all individual genes of N. benthamiana and the well studied N. tabacum (section Nicotiana) and their homologues in six diploid Nicotiana species. We generated sets of 44 424 and 65 457 phylogenetic trees of N. benthamiana and N. tabacum genes, respectively, each collectively called a phylome. Members of Nicotiana sections Noctiflorae and Sylvestres were represented as the species closest to N. benthamiana in most of the gene trees. Analyzing the gene trees of the phylome we: (i) dated the hybridization event giving rise to N. benthamiana to 4-5 MyA, and (ii) separated the subgenomes. We assigned 1.42 Gbp of the genome sequence to section Noctiflorae and 0.97 Gbp to section Sylvestres based on phylome analysis. In contrast, read mapping of the donor species did not succeed in separating the subgenomes of N. benthamiana. We show that the maternal progenitor of N. benthamiana was a member of section Noctiflorae, and confirm a member of section Sylvestres as paternal subgenome donor. We also demonstrate that the advanced stage of long-term genome diploidization in N. benthamiana is reflected in its subgenome organization. Taken together, our results underscore the usefulness of phylome analysis for subgenome characterization in hybrid species.

Examination of Genomic and Transcriptomic Alterations in a Morphologically Stable Line, MU1, Generated by Intergeneric Pollination

Interspecific hybridization creates genetic variation useful for crop improvement. However, whether pollen from a different genus affects the genomic stability and/or transcriptome of the recipient species during intergeneric pollination has not been investigated. Here, we crossed japonica rice cv. Z12 with the maize accession B73 (pollen donor) and obtained a morphologically stable line, MU1, exhibiting moderate dwarfism, higher tiller number, and increased grain weight compared with Z12. To reveal the genetic basis of these morphological changes in MU1, we performed whole-genome resequencing of MU1 and Z12. Compared with Z12, MU1 showed 107,250 single nucleotide polymorphisms (SNPs) and 23,278 insertion/deletions (InDels). Additionally, 5'-upstream regulatory regions (5'UTRs) of 429 and 309 differentially expressed genes (DEGs) in MU1 contained SNPs and InDels, respectively, suggesting that a subset of these DEGs account for the variation in 5'UTRs. Transcriptome analysis revealed 2190 DEGs in MU1 compared with Z12. Genes up-regulated in MU1 were mainly involved in photosynthesis, generation of precursor metabolites, and energy and cellular biosynthetic processes; whereas those down-regulated in MU1 were involved in plant hormone signal transduction pathway and response to stimuli and stress processes. Quantitative PCR (qPCR) further identified the expression levels of the up- or down-regulated gene in plant hormone signal transduction pathway. The expression level changes of plant hormone signal transduction pathway may be significant for plant growth and development. These findings suggest that mutations caused by intergeneric pollination could be the important reason for changes of MU1 in agronomic traits.

Developing Transposable Element Marker System for Molecular Breeding

Transposable element (TE) marker system was developed considering the useful properties of the transposable elements such as their large number in the animal and plant genomes, high rate of insertion polymorphism, and ease of detection. Various methods have been employed for developing a large number of TE markers in several crop plants for genomics studies. Here we describe some of these methods including the recent whole genome search. We also review the application of TE markers in molecular breeding.

Copy Number Variation of Transposable Elements in Thinopyrum intermedium and Its Diploid Relative Species

Diploid and polyploid wild species of Triticeae have complex relationships, and the understanding of their evolution and speciation could help to increase the usability of them in wheat breeding as a source of genetic diversity. The diploid species Pseudoroegneria spicata (St), Thinopyrum bessarabicum (Jb), Dasypyrum villosum (V) derived from a hypothetical common ancestor are considered to be possible subgenome donors in hexaploid species Th. intermedium (JrJvsSt, where indices r, v, and s stand for the partial relation to the genomes of Secale, Dasypyrum, and Pseudoroegneria, respectively). We quantified 10 families of transposable elements (TEs) in P. spicata, Th. bessarabicum, D. villosum (per one genome), and Th. intermedium (per one average subgenome) using the quantitative real time PCR assay and compared their abundance within the studied genomes as well as between them. Sabrina was the most abundant among all studied elements in P. spicata, D. villosum, and Th. intermedium, and among Ty3/Gypsy elements in all studied species. Among Ty1/Copia elements, Angela-A and WIS-A showed the highest and close abundance with the exception of D. villosum, and comprised the majority of all studied elements in Th. bessarabicum. Sabrina, BAGY2, and Angela-A showed similar abundance among diploids and in Th. intermedium hexaploid; Latidu and Barbara demonstrated sharp differences between diploid genomes. The relationships between genomes of Triticeae species based on the studied TE abundance and the role of TEs in speciation and polyploidization in the light of the current phylogenetic models is discussed.

Relaxed purifying selection in autopolyploids drives transposable element over-accumulation which provides variants for local adaptation

Polyploidization is frequently associated with increased transposable element (TE) content. However, what drives TE dynamics following whole genome duplication (WGD) and the evolutionary implications remain unclear. Here, we leverage whole-genome resequencing data available for ~300 individuals of Arabidopsis arenosa, a well characterized natural diploid-autotetraploid plant species, to address these questions. Based on 43,176 TE insertions we detect in these genomes, we demonstrate that relaxed purifying selection rather than transposition bursts is the main driver of TE over-accumulation after WGD. Furthermore, the increased pool of TE insertions in tetraploids is especially enriched within or near environmentally responsive genes. Notably, we show that the major flowering-time repressor gene FLC is disrupted by a TE insertion specifically in the rapid-cycling tetraploid lineage that colonized mainland railways. Together, our findings indicate that tetrasomy leads to an enhanced accumulation of genic TE insertions, some of which likely contribute to local adaptation.

Why transposable elements (TEs) accumulate in polyploids and the evolutionary implications remain unclear. Here, the authors show that following whole genome duplication, relaxed purifying selection is the main driver of TE over-accumulation, which provides variants for rapid local adaptation.

Interspecific hybridisation and LTR-retrotransposon mobilisation-related structural variation in plants: A case study

The dynamics of long-terminal-repeat retrotransposons in two poplar species (Populus deltoides and P. nigra) and in an interspecific hybrid, recently synthesized, were investigated by analyzing the genomic abundance and transcription levels of a collection of 828 full-length retroelements identified in the genome sequence of P. trichocarpa, all occurring also in the genomes of P. deltoides and P. nigra. Overall, genomic abundance and transcription levels of many retrotransposons in the hybrid resulted higher or lower than expected by calculating the mean of the parental values. A bioinformatics procedure was established to ascertain the occurrence of the same retrotransposon loci in the three genotypes. The results indicated that retrotransposon abundance variations between the hybrid and the mean value of the parents were due to i) co-segregation of retrotransposon high- or low-abundant haplotypes; ii) new retroelement insertions; iii) retrotransposon loss. Concerning retrotransposon expression, this was generally low, with only 14/828 elements over- or under-expressed in the hybrid than expected by calculating the mean of the parents. It is concluded that interspecific hybridisation between the two poplar species determine quantitative variation and differential expression of some retrotransposons, with possible consequences for the genetic differentiation of the hybrid.

ChIP-cloning analysis uncovers centromere-specific retrotransposons in Brassica nigra and reveals their rapid diversification in Brassica allotetraploids

Centromeres are indispensable functional units of chromosomes. The evolutionary mechanisms underlying the rapid evolution of centromeric repeats, especially those following polyploidy, remain unknown. In this study, we isolated centromeric sequences of Brassica nigra, a model diploid progenitor (B genome) of the allopolyploid species B. juncea (AB genome) and B. carinata (BC genome) by chromatin immunoprecipitation of nucleosomes containing the centromere-specific histone CENH3. Sequence analysis detected no centromeric satellite DNAs, and most B. nigra centromeric repeats were found to originate from Tyl/copia-class retrotransposons. In cytological analyses, six of the seven analyzed repeat clusters had no FISH signals in A or C genomes of the related diploid species B. rapa and B. oleracea. Notably, five repeat clusters had FISH signals in both A and B subgenomes in the tetraploid B. juncea. In the tetraploid B. carinata, only CL23 displayed three pairs of signals in terminal or interstitial regions of the C-derived chromosome, and no evidence of colonization of CLs onto C-subgenome centromeres was found in B. carinata. This observation suggests that centromeric repeats spread and proliferated between genomes after polyploidization. CL3 and CRB are likely ancient centromeric sequences arising prior to the divergence of diploid Brassica which have detected signals across the genus. And in allotetraploids B. juncea and B. carinata, the FISH signal intensity of CL3 and CRB differed among subgenomes. We discussed possible mechanisms for centromeric repeat divergence during Brassica speciation and polyploid evolution, thus providing insights into centromeric repeat establishment and targeting.

Parental transposable element loads influence their dynamics in young Nicotiana hybrids and allotetraploids

The genomic shock hypothesis suggests that allopolyploidy is associated with genome changes driven by transposable elements, as a response to imbalances between parental insertion loads. To explore this hypothesis, we compared three allotetraploids, Nicotiana arentsii, N. rustica and N. tabacum, which arose over comparable time frames from hybridisation between increasingly divergent diploid species. We used sequence-specific amplification polymorphism (SSAP) to compare the dynamics of six transposable elements in these allopolyploids, their diploid progenitors and in corresponding synthetic hybrids. We show that element-specific dynamics in young Nicotiana allopolyploids reflect their dynamics in diploid progenitors. Transposable element mobilisation is not concomitant with immediate genome merger, but occurs within the first generations of allopolyploid formation. In natural allopolyploids, such mobilisations correlate with imbalances in the repeat profile of the parental species, which increases with their genetic divergence. Other restructuring leading to locus loss is immediate, nonrandom and targeted at specific subgenomes, independently of cross orientation. The correlation between transposable element mobilisation in allopolyploids and quantitative imbalances in parental transposable element loads supports the genome shock hypothesis proposed by McClintock.

Footprints of parasitism in the genome of the parasitic flowering plant Cuscuta campestris

A parasitic lifestyle, where plants procure some or all of their nutrients from other living plants, has evolved independently in many dicotyledonous plant families and is a major threat for agriculture globally. Nevertheless, no genome sequence of a parasitic plant has been reported to date. Here we describe the genome sequence of the parasitic field dodder, Cuscuta campestris. The genome contains signatures of a fairly recent whole-genome duplication and lacks genes for pathways superfluous to a parasitic lifestyle. Specifically, genes needed for high photosynthetic activity are lost, explaining the low photosynthesis rates displayed by the parasite. Moreover, several genes involved in nutrient uptake processes from the soil are lost. On the other hand, evidence for horizontal gene transfer by way of genomic DNA integration from the parasite's hosts is found. We conclude that the parasitic lifestyle has left characteristic footprints in the C. campestris genome.

The genome structure of Arachis hypogaea (Linnaeus, 1753) and an induced Arachis allotetraploid revealed by molecular cytogenetics

Peanut, Arachis hypogaea (Linnaeus, 1753) is an allotetraploid cultivated plant with two subgenomes derived from the hybridization between two diploid wild species, A. duranensis (Krapovickas & W. C. Gregory, 1994) and A. ipaensis (Krapovickas & W. C. Gregory, 1994), followed by spontaneous chromosomal duplication. To understand genome changes following polyploidy, the chromosomes of A. hypogaea, IpaDur1, an induced allotetraploid (A. ipaensis × A. duranensis)^4x and the diploid progenitor species were cytogenetically compared. The karyotypes of the allotetraploids share the number and general morphology of chromosomes; DAPI⁺ bands pattern and number of 5S rDNA loci. However, one 5S rDNA locus presents a heteromorphic FISH signal in both allotetraploids, relative to corresponding progenitor. Whilst for A. hypogaea the number of 45S rDNA loci was equivalent to the sum of those present in the diploid species, in IpaDur1, two loci have not been detected. Overall distribution of repetitive DNA sequences was similar in both allotetraploids, although A. hypogaea had additional CMA₃⁺ bands and few slight differences in the LTR-retrotransposons distribution compared to IpaDur1. GISH showed that the chromosomes of both allotetraploids had preferential hybridization to their corresponding diploid genomes. Nevertheless, at least one pair of IpaDur1 chromosomes had a clear mosaic hybridization pattern indicating recombination between the subgenomes, clear evidence that the genome of IpaDur1 shows some instability comparing to the genome of A. hypogaea that shows no mosaic of subgenomes, although both allotetraploids derive from the same progenitor species. For some reasons, the chromosome structure of A. hypogaea is inherently more stable, or, it has been at least, partially stabilized through genetic changes and selection.

Mixed-Ploidy Species: Progress and Opportunities in Polyploid Research

Mixed-ploidy species harbor a unique form of genomic and phenotypic variation that influences ecological interactions, facilitates genetic divergence, and offers insights into the mechanisms of polyploid evolution. However, there have been few attempts to synthesize this literature. We review here research on the cytotype distribution, diversity, and dynamics of intensively studied mixed-ploidy species and consider the implications for understanding mechanisms of polyploidization such as cytotype formation, establishment, coexistence, and post-polyploid divergence. In general, mixed-ploidy species are unevenly represented among families: they exhibit high cytotype diversity, often within populations, and frequently comprise rare and odd-numbered ploidies. Odd-ploidies often occur in association with asexuality. We highlight research hypotheses and opportunities that take advantage of the unique properties of ploidy variation.

Impact of transposable elements on polyploid plant genomes

BACKGROUND

The growing wealth of knowledge on whole-plant genome sequences is highlighting the key role of transposable elements (TEs) in plant evolution, as a driver of drastic changes in genome size and as a source of an important number of new coding and regulatory sequences. Together with polyploidization events, TEs should thus be considered the major players in evolution of plants.

SCOPE

This review outlines the major mechanisms by which TEs impact plant genome evolution and how polyploidy events can affect these impacts, and vice versa. These include direct effects on genes, by providing them with new coding or regulatory sequences, an effect on the epigenetic status of the chromatin close to genes, and more subtle effects by imposing diverse evolutionary constraints to different chromosomal regions. These effects are particularly relevant after polyploidization events. Polyploidization often induces bursts of transposition probably due to a relaxation in their epigenetic control, and, in the short term, this can increase the rate of gene mutations and changes in gene regulation due to the insertion of TEs next to or into genes. Over longer times, TE bursts may induce global changes in genome structure due to inter-element recombination including losses of large genome regions and chromosomal rearrangements that reduce the genome size and the chromosome number as part of a process called diploidization.

CONCLUSIONS

TEs play an essential role in genome and gene evolution, in particular after polyploidization events. Polyploidization can induce TE activity that may explain part of the new phenotypes observed. TEs may also play a role in the diploidization that follows polyploidization events. However, the extent to which TEs contribute to diploidization and fractionation bias remains unclear. Investigating the multiple factors controlling TE dynamics and the nature of ancient and recent polyploid genomes may shed light on these processes.

LTR-retrotransposons in plants: Engines of evolution

LTR retrotransposons are the most abundant group of transposable elements (TEs) in plants. These elements can fall inside or close to genes, and therefore influence their expression and evolution. This review aims to examine how LTR retrotransposons, especially Ty1-copia elements, mediate gene regulation and evolution. Various stimuli, including polyploidization and biotic and abiotic elicitors, result in the transcription and movement of these retrotransposons, and can facilitate adaptation. The presence of cis-regulatory motifs in the LTRs are central to their stress-mediated responses and are shared with host stress-responsive genes, showing a complex evolutionary history in which TEs provide new regulatory units to genes. The presence of retrotransposon remnants in genes that are necessary for normal gene function, demonstrates the importance of exaptation and co-option, and is also a consequence of the abundance of these elements in plant genomes. Furthermore, insertions of LTR retrotransposons in and around genes provide potential for alternative splicing, epigenetic control, transduction, duplication and recombination. These characteristics can become an active part of the evolution of gene families as in the case of resistance genes (R-genes). The character of TEs as exclusively selfish is now being re-evaluated. Since genome-wide reprogramming via TEs is a long evolutionary process, the changes we can examine are case-specific and their fitness advantage may not be evident until TE-derived motifs and domains have been completely co-opted and fixed. Nevertheless, the presence of LTR retrotransposons inside genes and as part of gene promoter regions is consistent with their roles as engines of plant genome evolution.

Isolation and application of P genome-specific DNA sequences of Agropyron Gaertn. in Triticeae

Different types of P genome sequences and markers were developed, which could be used to analyze the evolution of P genome in Triticeae and identify precisely wheat- A. cristatum introgression lines. P genome of Agropyron Gaertn. plays an important role in Triticeae and could provide many desirable genes conferring high yield, disease resistance, and stress tolerance for wheat genetic improvement. Therefore, it is significant to develop specific sequences and functional markers of P genome. In this study, 126 sequences were isolated from the degenerate oligonucleotide primed-polymerase chain reaction (DOP-PCR) products of microdissected chromosome 6PS. Forty-eight sequences were identified as P genome-specific sequences by dot-blot hybridization and DNA sequences analysis. Among these sequences, 22 displayed the characteristics of retrotransposons, nine and one displayed the characteristics of DNA transposons and tandem repetitive sequence, respectively. Fourteen of 48 sequences were determined to distribute on different regions of P genome chromosomes by fluorescence in situ hybridization, and the distributing regions were as following: all over P genome chromosomes, centromeres, pericentromeric regions, distal regions, and terminal regions. We compared the P genome sequences with other genome sequences of Triticeae and found that the similar sequences of the P genome sequences were widespread in Triticeae, but differentiation occurred to various extents. Additionally, thirty-four molecular markers were developed from the P genome sequences, which could be used for analyzing the evolutionary relationship among 16 genomes of 18 species in Triticeae and identifying P genome chromatin in wheat-A. cristatum introgression lines. These results will not only facilitate the study of structure and evolution of P genome chromosomes, but also provide a rapid detecting tool for effective utilization of desirable genes of P genome in wheat improvement.

Profiling Transposable Elements and Their Epigenetic Effects in Non-model Species

Taking transposable elements into consideration in surveys of genetic and epigenetic variation remains challenging in species lacking a high-quality reference genome. Here, molecular techniques reducing genome complexity and specifically targeting restructuring and methylation changes in TE genome fractions are described. In particular, methyl-sensitive transposon display (MSTD) uses isoschizomers and PCR amplifications to assess the methylation environment of TE insertions. MSTD offers reliable insights into genome-wide epigenetic changes associated with TEs, especially when used together with similar techniques tracking random sequences.

Rapid proliferation and nucleolar organizer targeting centromeric retrotransposons in cotton

Centromeric chromatin in most eukaryotes is composed of highly repetitive centromeric retrotransposons and satellite repeats that are highly variable even among closely related species. The evolutionary mechanisms that underlie the rapid evolution of centromeric repeats remain unknown. To obtain insight into the evolution of centromeric repeats following polyploidy, we studied a model diploid progenitor (Gossypium raimondii, D-genome) of the allopolyploid (AD-genome) cottons, G. hirsutum and G. barbadense. Sequence analysis of chromatin-immunoprecipitated DNA showed that the G. raimondii centromeric repeats originated from retrotransposon-related sequences. Comparative analysis showed that nine of the 10 analyzed centromeric repeats were absent from the centromeres in the A-genome and related diploid species (B-, F- and G-genomes), indicating that they colonized the centromeres of D-genome lineage after the divergence of the A- and D- ancestral species or that they were ancestrally retained prior to the origin of Gossypium. Notably, six of the nine repeats were present in both the A- and D-subgenomes in tetraploid G. hirsutum, and increased in abundance in both subgenomes. This finding suggests that centromeric repeats may spread and proliferate between genomes subsequent to polyploidization. Two repeats, Gr334 and Gr359 occurred in both the centromeres and nucleolar organizer regions (NORs) in D- and AD-genome species, yet localized to just the NORs in A-, B-, F-, and G-genome species. Contained within is a story of an established centromeric repeat that is eliminated and allopolyploidization provides an opportunity for reinvasion and reestablishment, which broadens our evolutionary understanding behind the cycles of centromeric repeat establishment and targeting.

Moving through the Stressed Genome: Emerging Regulatory Roles for Transposons in Plant Stress Response

The recognition of a positive correlation between organism genome size with its transposable element (TE) content, represents a key discovery of the field of genome biology. Considerable evidence accumulated since then suggests the involvement of TEs in genome structure, evolution and function. The global genome reorganization brought about by transposon activity might play an adaptive/regulatory role in the host response to environmental challenges, reminiscent of McClintock's original 'Controlling Element' hypothesis. This regulatory aspect of TEs is also garnering support in light of the recent evidences, which project TEs as "distributed genomic control modules." According to this view, TEs are capable of actively reprogramming host genes circuits and ultimately fine-tuning the host response to specific environmental stimuli. Moreover, the stress-induced changes in epigenetic status of TE activity may allow TEs to propagate their stress responsive elements to host genes; the resulting genome fluidity can permit phenotypic plasticity and adaptation to stress. Given their predominating presence in the plant genomes, nested organization in the genic regions and potential regulatory role in stress response, TEs hold unexplored potential for crop improvement programs. This review intends to present the current information about the roles played by TEs in plant genome organization, evolution, and function and highlight the regulatory mechanisms in plant stress responses. We will also briefly discuss the connection between TE activity, host epigenetic response and phenotypic plasticity as a critical link for traversing the translational bridge from a purely basic study of TEs, to the applied field of stress adaptation and crop improvement.

Small RNA changes in synthetic Brassica napus

Small RNAs and microRNAs were found to vary extensively in synthetic Brassica napus and subsequent generations, accompanied by the activation of transposable elements in response to hybridization and polyploidization. Resynthesizing B. napus by hybridization and chromosome doubling provides an approach to create novel polyploids and increases the usable genetic variability in oilseed rape. Although many studies have shown that small RNAs (sRNAs) act as important factor during hybridization and polyploidization in plants, much less is known on how sRNAs change in synthetic B. napus, particularly in subsequent generations after formation. We performed high-throughput sequencing of sRNAs in S1-S4 generations of synthetic B. napus and in the homozygous B. oleracea and B. rapa parent lines. We found that the number of small RNAs (sRNAs) and microRNAs (miRNAs) doubled in synthetic B. napus relative to the parents. The proportions of common sRNAs detected varied from the S1 to S4 generations, suggesting sRNAs are unstable in synthetic B. napus. The majority of miRNAs (67.2 %) were non-additively expressed in the synthesized Brassica allotetraploid, and 33.3 % of miRNAs were novel in the resynthesized B. napus. The percentage of miRNAs derived from transposable elements (TEs) also increased, indicating transposon activation and increased transposon-associated miRNA production in response to hybridization and polyploidization. The number of target genes for each miRNA in the synthesized Brassica allotetraploid was doubled relative to the parents, enhancing the complexity of gene expression regulation. The potential roles of miRNAs and their targets are discussed. Our data demonstrate generational changes in sRNAs and miRNAs in synthesized B. napus.

Polyploidy: Pitfalls and paths to a paradigm

Investigators have long searched for a polyploidy paradigm-rules or principles that might be common following polyploidization (whole-genome duplication, WGD). Here we attempt to integrate what is known across the more thoroughly investigated polyploid systems on topics ranging from genetics to ecology. We found that while certain rules may govern gene retention and loss, systems vary in the prevalence of gene silencing vs. homeolog loss, chromosomal change, the presence of a dominant genome (in allopolyploids), and the relative importance of hybridization vs. genome doubling per se. In some lineages, aspects of polyploidization are repeated across multiple origins, but in other species multiple origins behave more stochastically in terms of genetic and phenotypic change. Our investigation also reveals that the path to synthesis is hindered by numerous gaps in our knowledge of even the best-known systems. Particularly concerning is the absence of linkage between genotype and phenotype. Moreover, most recent studies have focused on the genetic and genomic attributes of polyploidy, but rarely is there an ecological or physiological context. To promote a path to a polyploidy paradigm (or paradigms), we propose a major community goal over the next 10-20 yr to fill the gaps in our knowledge of well-studied polyploids. Before a meaningful synthesis is possible, more complete data sets are needed for comparison-systems that include comparable genetic, genomic, chromosomal, proteomic, as well as morphological, physiological, and ecological data. Also needed are more natural evolutionary model systems, as most of what we know about polyploidy continues to come from a few crop and genetic models, systems that often lack the ecological context inherent in natural systems and necessary for understanding the drivers of biodiversity.

Transposable element evolution in the allotetraploid Capsella bursa-pastoris

PREMISE OF THE STUDY

Shifts in ploidy affect the evolutionary dynamics of genomes in a myriad of ways. Population genetic theory predicts that transposable element (TE) proliferation may follow because the genomewide efficacy of selection should be reduced and the increase in gene copies may mask the deleterious effects of TE insertions. Moreover, in allopolyploids, TEs may further accumulate because of hybrid breakdown of TE silencing. However, to date the evidence of TE proliferation following an increase in ploidy is mixed, and the relative importance of relaxed selection vs. silencing breakdown remains unclear.

METHODS

We used high-coverage whole-genome sequence data to evaluate the abundance, genomic distribution, and population frequencies of TEs in the self-fertilizing recent allotetraploid Capsella bursa-pastoris (Brassicaceae). We then compared the C. bursa-pastoris TE profile with that of its two parental diploid species, outcrossing C. grandiflora and self-fertilizing C. orientalis.

KEY RESULTS

We found no evidence that C. bursa-pastoris has experienced a large genomewide proliferation of TEs relative to its parental species. However, when centromeric regions are excluded, we found evidence of significantly higher abundance of retrotransposons in C. bursa-pastoris along the gene-rich chromosome arms compared with C. grandiflora and C. orientalis.

CONCLUSIONS

The lack of a genomewide effect of allopolyploidy on TE abundance, combined with the increases TE abundance in gene-rich regions, suggests that relaxed selection rather than hybrid breakdown of host silencing explains the TE accumulation in C. bursa-pastoris.

Direct vs. indirect effects of whole-genome duplication on prezygotic isolation in Chamerion angustifolium: Implications for rapid speciation

PREMISE OF THE STUDY

The depiction of polyploid speciation as instantaneous implies that strong prezygotic and postzygotic isolation form as a direct result of whole-genome duplication. However, the direct vs. indirect contributions of genome duplication to phenotypic divergence and prezygotic isolation are rarely quantified across multiple reproductive barriers.

METHODS

We compared the phenotypic differences between diploid and both naturally occurring and synthesized tetraploids (neotetraploids) of the plant Chamerion angustifolium. Using this information and additional published values for this species, we compared the magnitude of isolation (ecological, flowering, pollinator, and gametic) between diploids and natural-occurring tetraploids to that between diploids and neotetraploids.

KEY RESULTS

Differences among ploidy cytotypes were observed for eight of 12 vegetative and reproductive traits measured. Neotetraploids resembled diploids but differed from natural tetraploids with respect to four traits, including flowering time and plant height. Diploid-neotetraploid (2x-4xneo) experimental arrays exhibited lower pollinator fidelity to cytotype and seed set compared with 2x-4xnat arrays. Based on these results and published evidence, reproductive isolation between diploids and neotetraploids across all four life stages averaged 0.48 and deviated significantly from that between diploids and natural tetraploids (RI = 0.96).

CONCLUSIONS

Genome duplication causes phenotypic shifts and contributes directly to prezygotic isolation for some barriers (gametic isolation) but cannot account for the cumulative isolation from diploids observed in natural tetraploids. Therefore, conditions for species formation through genome duplication are not necessarily instantaneous and selection to strengthen prezygotic barriers in young polyploids is critical for the establishment of polyploid species in sympatry.

The origin of exon 3 skipping of paternal GLOBOSA pre-mRNA in some Nicotiana tabacum lines correlates with a point mutation of the very last nucleotide of the exon

In plants, genome duplication followed by genome diversification and selection is recognized as a major evolutionary process. Rapid epigenetic and genetic changes that affect the transcription of parental genes are frequently observed after polyploidization. The pattern of alternative splicing is also frequently altered, yet the related molecular processes remain largely unresolved. Here, we study the inheritance and expression of parental variants of three floral organ identity genes in allotetraploid tobacco. DEFICIENS and GLOBOSA are B-class genes, and AGAMOUS is a C-class gene. Parental variants of these genes were found to be maintained in the tobacco genome, and the respective mRNAs were present in flower buds in comparable amounts. However, among five tobacco cultivars, we identified two in which the majority of paternal GLOBOSA pre-mRNA transcripts undergo exon 3 skipping, producing an mRNA with a premature termination codon. At the DNA level, we identified a G-A transition at the very last position of exon 3 in both cultivars. Although alternative splicing resulted in a dramatic decrease in full-length paternal GLOBOSA mRNA, no phenotypic effect was observed. Our finding likely serves as an example of the initiation of homoeolog diversification in a relatively young polyploid genome.

Effects of genome structure variation, homeologous genes and repetitive DNA on polyploid crop research in the age of genomics

Compared to diploid species, allopolyploid crop species possess more complex genomes, higher productivity, and greater adaptability to changing environments. Next generation sequencing techniques have produced high-density genetic maps, whole genome sequences, transcriptomes and epigenomes for important polyploid crops. However, several problems interfere with the full application of next generation sequencing techniques to these crops. Firstly, different types of genomic variation affect sequence assembly and QTL mapping. Secondly, duplicated or homoeologous genes can diverge in function and then lead to emergence of many minor QTL, which increases difficulties in fine mapping, cloning and marker assisted selection. Thirdly, repetitive DNA sequences arising in polyploid crop genomes also impact sequence assembly, and are increasingly being shown to produce small RNAs to regulate gene expression and hence phenotypic traits. We propose that these three key features should be considered together when analyzing polyploid crop genomes. It is apparent that dissection of genomic structural variation, elucidation of the function and mechanism of interaction of homoeologous genes, and investigation of the de novo roles of repeat sequences in agronomic traits are necessary for genomics-based crop breeding in polyploids.

Genome reorganization in F1 hybrids uncovers the role of retrotransposons in reproductive isolation

Interspecific hybridization leads to new interactions among divergent genomes, revealing the nature of genetic incompatibilities having accumulated during and after the origin of species. Conflicts associated with misregulation of transposable elements (TEs) in hybrids expectedly result in their activation and genome-wide changes that may be key to species boundaries. Repetitive genomes of wild wheats have diverged under differential dynamics of specific long terminal repeat retrotransposons (LTR-RTs), offering unparalleled opportunities to address the underpinnings of plant genome reorganization by selfish sequences. Using reciprocal F1 hybrids between three Aegilops species, restructuring and epigenetic repatterning was assessed at random and LTR-RT sequences with amplified fragment length polymorphism and sequence-specific amplified polymorphisms as well as their methylation-sensitive counterparts, respectively. Asymmetrical reorganization of LTR-RT families predicted to cause conflicting interactions matched differential survival of F1 hybrids. Consistent with the genome shock model, increasing divergence of merged LTR-RTs yielded higher levels of changes in corresponding genome fractions and lead to repeated reorganization of LTR-RT sequences in F1 hybrids. Such non-random reorganization of hybrid genomes is coherent with the necessary repression of incompatible TE loci in support of hybrid viability and indicates that TE-driven genomic conflicts may represent an overlooked factor supporting reproductive isolation.

Transcriptionally active LTR retrotransposons in Eucalyptus genus are differentially expressed and insertionally polymorphic

BACKGROUND

In Eucalyptus genus, studies on genome composition and transposable elements (TEs) are particularly scarce. Nearly half of the recently released Eucalyptus grandis genome is composed by retrotransposons and this data provides an important opportunity to understand TE dynamics in Eucalyptus genome and transcriptome.

RESULTS

We characterized nine families of transcriptionally active LTR retrotransposons from Copia and Gypsy superfamilies in Eucalyptus grandis genome and we depicted genomic distribution and copy number in two Eucalyptus species. We also evaluated genomic polymorphism and transcriptional profile in three organs of five Eucalyptus species. We observed contrasting genomic and transcriptional behavior in the same family among different species. RLC_egMax_1 was the most prevalent family and RLC_egAngela_1 was the family with the lowest copy number. Most families of both superfamilies have their insertions occurring <3 million years, except one Copia family, RLC_egBianca_1. Protein theoretical models suggest different properties between Copia and Gypsy domains. IRAP and REMAP markers suggested genomic polymorphisms among Eucalyptus species. Using EST analysis and qRT-PCRs, we observed transcriptional activity in several tissues and in all evaluated species. In some families, osmotic stress increases transcript values.

CONCLUSION

Our strategy was successful in isolating transcriptionally active retrotransposons in Eucalyptus, and each family has a particular genomic and transcriptional pattern. Overall, our results show that retrotransposon activity have differentially affected genome and transcriptome among Eucalyptus species.

Short-term hybridisation activates Tnt1 and Tto1 Copia retrotransposons in wild tuber-bearing Solanum species

Interspecific hybridisation in tuber-bearing species of Solanum is a common phenomenon and represents an important source of variability, crucial for adaptation and speciation of potato species. In this regard, the effects of interspecific hybridisation on retrotransposon families present in the genomes, and their consequent effects on generation of genetic variability in wild tuber-bearing Solanum species, are poorly characterised. The aim of this study was to analyse the activity of retrotransposons in inter- and intraspecific hybrids between S. kurtzianum and S. microdontum, obtained by controlled crosses, and the effects on morphological, genetic and epigenetic variability. For genetic and epigenetic analysis, S-SAP (sequence-specific amplification polymorphism) and TMD (transposon methylation display) techniques were used, respectively, with specific primers for Tnt1 and Tto1 retrotransposon families (Order LTR, Superfamily Copia). The results indicate that at morphological level, interspecific hybrid genotypes differ from their parental species, whereas derived intraspecific hybrids do not. In both cases, we observed significant reductions in pollen grain viability, and a negative correlation with Tnt1 mobility. Both retrotransposons, Tto1 and Tnt1, were mobilised in the genotypes analysed, with mobility ranging from 0 to 7.8%. Furthermore, at the epigenetic level, demethylation was detected in the vicinity of Tnt1 and Tto1 in the hybrids compared with the parental genotypes. These patterns were positively correlated with the activity of the retrotransposons. The results suggest a possible mechanism through which hybridisation events generate genetic variability in tuber-bearing species of Solanum through retrotranposon activation.

Sequence-based ultra-dense genetic and physical maps reveal structural variations of allopolyploid cotton genomes

Background

SNPs are the most abundant polymorphism type, and have been explored in many crop genomic studies, including rice and maize. SNP discovery in allotetraploid cotton genomes has lagged behind that of other crops due to their complexity and polyploidy. In this study, genome-wide SNPs are detected systematically using next-generation sequencing and efficient SNP genotyping methods, and used to construct a linkage map and characterize the structural variations in polyploid cotton genomes.

Results

We construct an ultra-dense inter-specific genetic map comprising 4,999,048 SNP loci distributed unevenly in 26 allotetraploid cotton linkage groups and covering 4,042 cM. The map is used to order tetraploid cotton genome scaffolds for accurate assembly of G. hirsutum acc. TM-1. Recombination rates and hotspots are identified across the cotton genome by comparing the assembled draft sequence and the genetic map. Using this map, genome rearrangements and centromeric regions are identified in tetraploid cotton by combining information from the publicly-available G. raimondii genome with fluorescent in situ hybridization analysis.

Conclusions

We report the genotype-by-sequencing method used to identify millions of SNPs between G. hirsutum and G. barbadense. We construct and use an ultra-dense SNP map to correct sequence mis-assemblies, merge scaffolds into pseudomolecules corresponding to chromosomes, detect genome rearrangements, and identify centromeric regions in allotetraploid cottons. We find that the centromeric retro-element sequence of tetraploid cotton derived from the D subgenome progenitor might have invaded the A subgenome centromeres after allotetrapolyploid formation. This study serves as a valuable genomic resource for genetic research and breeding of cotton.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0678-1) contains supplementary material, which is available to authorized users.

Construction of a linkage map based on retrotransposon insertion polymorphisms in sweetpotato via high-throughput sequencing

Sweetpotato (Ipomoea batatas L.) is an outcrossing hexaploid species with a large number of chromosomes (2n = 6x = 90). Although sweetpotato is one of the world's most important crops, genetic analysis of the species has been hindered by its genetic complexity combined with the lack of a whole genome sequence. In the present study, we constructed a genetic linkage map based on retrotransposon insertion polymorphisms using a mapping population derived from a cross between 'Purple Sweet Lord' (PSL) and '90IDN-47' cultivars. High-throughput sequencing and subsequent data analyses identified many Rtsp-1 retrotransposon insertion sites, and their allele dosages (simplex, duplex, triplex, or double-simplex) were determined based on segregation ratios in the mapping population. Using a pseudo-testcross strategy, 43 and 47 linkage groups were generated for PSL and 90IDN-47, respectively. Interestingly, most of these insertions (~90%) were present in a simplex manner, indicating their utility for linkage map construction in polyploid species. Additionally, our approach led to savings of time and labor for genotyping. Although the number of markers herein was insufficient for map-based cloning, our trial analysis exhibited the utility of retrotransposon-based markers for linkage map construction in sweetpotato.

Polyploidy and novelty: Gottlieb's legacy

Nearly four decades ago, Roose & Gottlieb (Roose & Gottlieb 1976 Evolution 30, 818-830. (doi:10.2307/2407821)) showed that the recently derived allotetraploids Tragopogon mirus and T. miscellus combined the allozyme profiles of their diploid parents (T. dubius and T. porrifolius, and T. dubius and T. pratensis, respectively). This classic paper addressed the link between genotype and biochemical phenotype and documented enzyme additivity in allopolyploids. Perhaps more important than their model of additivity, however, was their demonstration of novelty at the biochemical level. Enzyme multiplicity-the production of novel enzyme forms in the allopolyploids-can provide an extensive array of polymorphism for a polyploid individual and may explain, for example, the expanded ranges of polyploids relative to their diploid progenitors. In this paper, we extend the concept of evolutionary novelty in allopolyploids to a range of genetic and ecological features. We observe that the dynamic nature of polyploid genomes-with alterations in gene content, gene number, gene arrangement, gene expression and transposon activity-may generate sufficient novelty that every individual in a polyploid population or species may be unique. Whereas certain combinations of these features will undoubtedly be maladaptive, some unique combinations of newly generated variation may provide tremendous evolutionary potential and adaptive capabilities.