Retrotransposon BARE-1 and Its Role in Genome Evolution in the Genus Hordeum

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.

Evolution of complex genome architecture in gymnosperms

Gymnosperms represent an ancient lineage that diverged from early spermatophytes during the Devonian. The long fossil records and low diversity in living species prove their complex evolutionary history, which included ancient radiations and massive extinctions. Due to their ultra-large genome size, the whole-genome assembly of gymnosperms has only generated in the past 10 years and is now being further expanded into more taxonomic representations. Here, we provide an overview of the publicly available gymnosperm genome resources and discuss their assembly quality and recent findings in large genome architectures. In particular, we describe the genomic features most related to changes affecting the whole genome. We also highlight new realizations relative to repetitive sequence dynamics, paleopolyploidy, and long introns. Based on the results of relevant genomic studies of gymnosperms, we suggest additional efforts should be made toward exploring the genomes of medium-sized (5–15 gigabases) species. Lastly, more comparative analyses among high-quality assemblies are needed to understand the genomic shifts and the early species diversification of seed plants.

Constitutive Heterochromatin in Eukaryotic Genomes: A Mine of Transposable Elements

Transposable elements (TEs) are abundant components of constitutive heterochromatin of the most diverse evolutionarily distant organisms. TEs enrichment in constitutive heterochromatin was originally described in the model organism Drosophila melanogaster, but it is now considered as a general feature of this peculiar portion of the genomes. The phenomenon of TE enrichment in constitutive heterochromatin has been proposed to be the consequence of a progressive accumulation of transposable elements caused by both reduced recombination and lack of functional genes in constitutive heterochromatin. However, this view does not take into account classical genetics studies and most recent evidence derived by genomic analyses of heterochromatin in Drosophila and other species. In particular, the lack of functional genes does not seem to be any more a general feature of heterochromatin. Sequencing and annotation of Drosophila melanogaster constitutive heterochromatin have shown that this peculiar genomic compartment contains hundreds of transcriptionally active genes, generally larger in size than that of euchromatic ones. Together, these genes occupy a significant fraction of the genomic territory of heterochromatin. Moreover, transposable elements have been suggested to drive the formation of heterochromatin by recruiting HP1 and repressive chromatin marks. In addition, there are several pieces of evidence that transposable elements accumulation in the heterochromatin might be important for centromere and telomere structure. Thus, there may be more complexity to the relationship between transposable elements and constitutive heterochromatin, in that different forces could drive the dynamic of this phenomenon. Among those forces, preferential transposition may be an important factor. In this article, we present an overview of experimental findings showing cases of transposon enrichment into the heterochromatin and their positive evolutionary interactions with an impact to host genomes.

The genome sequence provides insights into salt tolerance of Achnatherum splendens (Gramineae), a constructive species of alkaline grassland

Achnatherum splendens Trin. (Gramineae) is a constructive species of the arid grassland ecosystem in Northwest China and is a major forage grass. It has good tolerance of salt and drought stress in alkaline habitats. Here, we report its chromosome-level genome, determined through a combination of Illumina HiSeq sequencing, PacBio sequencing and Hi-C technology. The final assembly of the ~1.17 Gb genome sequence had a super-scaffold N50 of 40.3 Mb. A total of 57 374 protein-coding genes were annotated, of which 54 426 (94.5%) genes have functional protein annotations. Approximately 735 Mb (62.37%) of the assembly were identified as repetitive elements, and among these, LTRs (40.53%) constitute the highest proportion, having made a major contribution to the expansion of genome size in A. splendens. Phylogenetic analysis revealed that A. splendens diverged from the Brachypodium distachyon-Hordeum vulgare-Aegilops tauschii subclade around 37 million years ago (Ma) and that a clade comprising these four species diverged from the Phyllostachys edulis clade ~47 Ma. Genomic synteny indicates that A. splendens underwent an additional species-specific whole-genome duplication (WGD) 18-20 Ma, which further promoted an increase in copies of numerous saline-alkali-related gene families in the A. splendens genome. By transcriptomic analysis, we further found that many of these duplicated genes from this extra WGD exhibited distinct functional divergence in response to salt stress. This WGD, therefore, contributed to the strong resistance to salt stress and widespread arid adaptation of A. splendens.

Chromosome and Genome Diversity in the Genus Trifolium (Fabaceae)

Trifolium L. is an economically important genus that is characterized by variable karyotypes relating to its ploidy level and basic chromosome numbers. The advent of genomic resources combined with molecular cytogenetics provides an opportunity to develop our understanding of plant genomes in general. Here, we summarize the current state of knowledge on Trifolium genomes and chromosomes and review methodologies using molecular markers that have contributed to Trifolium research. We discuss possible future applications of cytogenetic methods in research on the Trifolium genome and chromosomes.

The Dynamism of Transposon Methylation for Plant Development and Stress Adaptation

Plant development processes are regulated by epigenetic alterations that shape nuclear structure, gene expression, and phenotypic plasticity; these alterations can provide the plant with protection from environmental stresses. During plant growth and development, these processes play a significant role in regulating gene expression to remodel chromatin structure. These epigenetic alterations are mainly regulated by transposable elements (TEs) whose abundance in plant genomes results in their interaction with genomes. Thus, TEs are the main source of epigenetic changes and form a substantial part of the plant genome. Furthermore, TEs can be activated under stress conditions, and activated elements cause mutagenic effects and substantial genetic variability. This introduces novel gene functions and structural variation in the insertion sites and primarily contributes to epigenetic modifications. Altogether, these modifications indirectly or directly provide the ability to withstand environmental stresses. In recent years, many studies have shown that TE methylation plays a major role in the evolution of the plant genome through epigenetic process that regulate gene imprinting, thereby upholding genome stability. The induced genetic rearrangements and insertions of mobile genetic elements in regions of active euchromatin contribute to genome alteration, leading to genomic stress. These TE-mediated epigenetic modifications lead to phenotypic diversity, genetic variation, and environmental stress tolerance. Thus, TE methylation is essential for plant evolution and stress adaptation, and TEs hold a relevant military position in the plant genome. High-throughput techniques have greatly advanced the understanding of TE-mediated gene expression and its associations with genome methylation and suggest that controlled mobilization of TEs could be used for crop breeding. However, development application in this area has been limited, and an integrated view of TE function and subsequent processes is lacking. In this review, we explore the enormous diversity and likely functions of the TE repertoire in adaptive evolution and discuss some recent examples of how TEs impact gene expression in plant development and stress adaptation.

The Evolutionary Patterns of Genome Size in Ensifera (Insecta: Orthoptera)

Genomic size variation has long been a focus for biologists. However, due to the lack of genome size data, the mechanisms behind this variation and the biological significance of insect genome size are rarely studied systematically. The detailed taxonomy and phylogeny of the Ensifera, as well as the extensive documentation concerning their morphological, ecological, behavioral, and distributional characteristics, make them a strong model for studying the important scientific problem of genome size variation. However, data on the genome size of Ensifera are rather sparse. In our study, we used flow cytometry to determine the genome size of 32 species of Ensifera, the smallest one being only 1C = 0.952 pg with the largest species up to 1C = 19.135 pg, representing a 20-fold range. This provides a broader blueprint for the genome size variation of Orthoptera than was previously available. We also completed the assembly of nine mitochondrial genomes and combined mitochondrial genome data from public databases to construct phylogenetic trees containing 32 species of Ensifera and three outgroups. Based on these inferred phylogenetic trees, we detected the phylogenetic signal of genome size variation in Ensifera and found that it was strong in both males and females. Phylogenetic comparative analyses revealed that there were no correlations between genome size and body size or flight ability in Tettigoniidae. Reconstruction of ancestral genome size revealed that the genome size of Ensifera evolved in a complex pattern, in which the genome size of the grylloid clade tended to decrease while that of the non-grylloid clade expanded significantly albeit with fluctuations. However, the evolutionary mechanisms underlying variation of genome size in Ensifera are still unknown.

Chromosome-scale genome assembly of areca palm (Areca catechu)

Areca palm (Areca catechu L.; family Arecaceae) is an important tropical medicinal crop and is also used for masticatory and religious purposes in Asia. Improvements to areca properties made by traditional breeding tools have been very slow, and further advances in its cultivation and practical use require genomic information, which is still unavailable. Here, we present a chromosome-scale reference genome assembly for areca by combining Illumina and PacBio data with Hi-C mapping technologies, covering the predicted A. catechu genome length (2.59 Gb, variety "Reyan#1") to an estimated 240× read depth. The assembly was 2.51 Gb in length with a scaffold N50 of 1.7Mb. The scaffolds were then further assembled into 16 pseudochromosomes, with an N50 of 172 Mb. Transposable elements comprised 80.37% of the areca genome, and 68.68% of them were long-terminal repeat retrotransposon elements. The areca palm genome was predicted to harbour 31,571 protein-coding genes and overall, 92.92% of genes were functionally annotated, including enriched and expanded families of genes responsible for biosynthesis of flavonoid, anthocyanin, monoterpenoid and their derivatives. Comparative analyses indicated that A. catechu probably diverged from its close relatives Elaeis guineensis and Cocos nucifera approximately 50.3 million years ago (Ma). Two whole genome duplication events in areca palm were found to be shared by palms and monocots, respectively. This genome assembly and associated resources represents an important addition to the palm genomics community and will be a valuable resource that will facilitate areca palm breeding and improve our understanding of areca palm biology and evolution.

Full-length LTR retroelements in Capsicum annuum revealed a few species-specific family bursts with insertional preferences

Capsicum annuum is a species that has undergone an expansion of the size of its genome caused mainly by the amplification of repetitive DNA sequences, including mobile genetic elements. Based on information obtained from sequencing the genome of pepper, the estimated fraction of retroelements is approximately 81%, and previous results revealed an important contribution of lineages derived from Gypsy superfamily. However, the dynamics of the retroelements in the C. annuum genome is poorly understood. In this way, the present work seeks to investigate the phylogenetic diversity and genomic abundance of the families of autonomous (complete and intact) LTR retroelements from C. annuum and inspect their distribution along its chromosomes. In total, we identified 1151 structurally full-length retroelements (340 Copia; 811 Gypsy) grouped in 124 phylogenetic families in the base of their retrotranscriptase. All the evolutive lineages of LTR retroelements identified in plants were present in pepper; however, three of them comprise 83% of the entire LTR retroelements population, the lineages Athila, Del/Tekay, and Ale/Retrofit. From them, only three families represent 70.8% of the total number of the identified retroelements. A massive family-specific wave of amplification of two of them occurred in the last 0.5 Mya (GypsyCa_16; CopiaCa_01), whereas the third is more ancient and occurred 3.0 Mya (GypsyCa_13). Fluorescent in situ hybridization performed with family and lineage-specific probes revealed contrasting patterns of chromosomal affinity. Our results provide a database of the populations LTR retroelements specific to C. annuum genome. The most abundant families were analyzed according to chromosome insertional preferences, suppling useful tools to the design of retroelement-based markers specific to the species.

Clustered and dispersed chromosomal distribution of the two classes of Revolver transposon family in rye (Secale cereale)

The chromosomal locations of a new class of Revolver transposon-like elements were analyzed by using FISH method on the metaphase chromosome in somatic cell division of the rye cultivar Petkus. First, the Revolver standard element probe λ2 was weakly hybridized throughout the rye chromosome, and comparatively large interstitial signals spotted with a dot shape were detected together with several telomeric regions. The dot shape interstitial signal was stably detected at one site on Chromosome (Chr) 1R (middle part of the interstitial region of the short arm), three sites on Chr 2R (distal part of the interstitial region and adjacent to the centromere on the short arm, middle part of the interstitial region of the long arm), and two sites on Chr 5R (middle part of the interstitial region and adjacent to the centromere on the long arm). The Revolver λ2 probe was effective for identification of 1R, 2R, and 5R chromosomes. On the other hand, Revolver nonautonomous element-specific L626-BARE-100 probe was strongly distributed throughout the rye chromosomes, and considerable numbers and diverse lengths of transcripts were detected by RT-PCR. Although the standard elements were found in localized clusters, the nonautonomous elements tended to be dispersed throughout the genome. Clustered nature of Revolver is a significantly rare case in genomics.

Coconut genome assembly enables evolutionary analysis of palms and highlights signaling pathways involved in salt tolerance

Coconut (Cocos nucifera) is the emblematic palm of tropical coastal areas all around the globe. It provides vital resources to millions of farmers. In an effort to better understand its evolutionary history and to develop genomic tools for its improvement, a sequence draft was recently released. Here, we present a dense linkage map (8402 SNPs) aiming to assemble the large genome of coconut (2.42 Gbp, 2n = 32) into 16 pseudomolecules. As a result, 47% of the sequences (representing 77% of the genes) were assigned to 16 linkage groups and ordered. We observed segregation distortion in chromosome Cn15, which is a signature of strong selection among pollen grains, favouring the maternal allele. Comparing our results with the genome of the oil palm Elaeis guineensis allowed us to identify major events in the evolutionary history of palms. We find that coconut underwent a massive transposable element invasion in the last million years, which could be related to the fluctuations of sea level during the glaciations at Pleistocene that would have triggered a population bottleneck. Finally, to better understand the facultative halophyte trait of coconut, we conducted an RNA-seq experiment on leaves to identify key players of signaling pathways involved in salt stress response. Altogether, our findings represent a valuable resource for the coconut breeding community.

Screening diversity and distribution of Copia retrotransposons reveals a specific amplification of BARE1 elements in genomes of the polyploid Hordeum murinum complex

We explored diversity, distribution and evolutionary dynamics of Ty1-Copia retrotransposons in the genomes of the Hordeum murinum polyploid complex and related taxa. Phylogenetic and fluorescent in situ hybridization (FISH) analyses of reverse transcriptase sequences identified four Copia families in these genomes: the predominant BARE1 (including three groups or subfamilies, A, B and C), and the less represented RIRE1, IKYA and TAR-1. Within the BARE1 family, BARE1-A elements and a subgroup of BARE1-B elements (named B1) have proliferated in the allopolyploid members of the H. murinum complex (H. murinum and H. leporinum), and in their extant diploid progenitor, subsp. glaucum. Moreover, we found a specific amplification of BARE1-B elements within each Hordeum species surveyed. The low occurrence of RIRE1, IKYA and TAR-1 elements in the allopolyploid cytotypes suggests that they are either weakly represented or highly degenerated in their diploid progenitors. The results demonstrate that BARE1-A and BARE1-B1 Copia elements are particularly well represented in the genomes of the H. murinum complex and constitute its genomic hallmark. No BARE1-A and -B1 homologs were detected in the reference barley genome. The similar distribution of RT-Copia probes across chromosomes of diploid, tetraploid and hexaploid taxa of the murinum complex shows no evidence of proliferation following polyploidization.

Low Long Terminal Repeat (LTR)-Retrotransposon Expression in Leaves of the Marine Phanerogam Posidonia Oceanica L

Seagrasses as Posidonia oceanica reproduce mostly by vegetative propagation, which can reduce genetic variability within populations. Since, in clonally propagated species, insurgence of genetic variability can be determined by the activity of transposable elements, we have estimated the activity of such repeat elements by measuring their expression level in the leaves of plants from a Mediterranean site, for which Illumina complementary DNA (cDNA) sequence reads (produced from RNAs isolated by leaves of plants from deep and shallow meadows) were publicly available. Firstly, we produced a collection of retrotransposon-related sequences and then mapped Illumina cDNA reads onto these sequences. With this approach, it was evident that Posidonia retrotransposons are, in general, barely expressed; only nine elements resulted transcribed at levels comparable with those of reference genes encoding tubulins and actins. Differences in transcript abundance were observed according to the superfamily and the lineage to which the retrotransposons belonged. Only small differences were observed between retrotransposon expression levels in leaves of shallow and deep Posidonia meadow stands, whereas one TAR/Tork element resulted differentially expressed in deep plants exposed to heat. It can be concluded that, in P. oceanica, the contribution of retrotransposon activity to genetic variability is reduced, although the nine specific active elements could actually produce new structural variations.

Genome assembly provides insights into the genome evolution and flowering regulation of orchardgrass

Orchardgrass (Dactylis glomerata L.) is an important forage grass for cultivating livestock worldwide. Here, we report an ~1.84-Gb chromosome-scale diploid genome assembly of orchardgrass, with a contig N50 of 0.93 Mb, a scaffold N50 of 6.08 Mb and a super-scaffold N50 of 252.52 Mb, which is the first chromosome-scale assembled genome of a cool-season forage grass. The genome includes 40 088 protein-coding genes, and 69% of the assembled sequences are transposable elements, with long terminal repeats (LTRs) being the most abundant. The LTRretrotransposons may have been activated and expanded in the grass genome in response to environmental changes during the Pleistocene between 0 and 1 million years ago. Phylogenetic analysis reveals that orchardgrass diverged after rice but before three Triticeae species, and evolutionarily conserved chromosomes were detected by analysing ancient chromosome rearrangements in these grass species. We also resequenced the whole genome of 76 orchardgrass accessions and found that germplasm from Northern Europe and East Asia clustered together, likely due to the exchange of plants along the 'Silk Road' or other ancient trade routes connecting the East and West. Last, a combined transcriptome, quantitative genetic and bulk segregant analysis provided insights into the genetic network regulating flowering time in orchardgrass and revealed four main candidate genes controlling this trait. This chromosome-scale genome and the online database of orchardgrass developed here will facilitate the discovery of genes controlling agronomically important traits, stimulate genetic improvement of and functional genetic research on orchardgrass and provide comparative genetic resources for other forage grasses.

A computational genome-wide analysis of long terminal repeats retrotransposon expression in sunflower roots (Helianthus annuus L.)

Long terminal repeats (LTR) retrotransposons have a major role in determining genome size, structure and function, thanks to their ability to transpose. We performed a meta-analysis of LTR-retrotransposon expression in roots of sunflower plantlets treated with different plant hormones, chemicals and NaCl. By using Illumina cDNA libraries, available from public repositories, we measured the number of reads matching the retrotranscriptase domains isolated from a whole genome library of retrotransposons. LTR-retrotransposons resulted in general barely expressed, except for 4 elements, all belonging to the AleII lineage, which showed high transcription levels in roots of both control and treated plants. The expression of retrotransposons in treated plants was slightly higher than in the control. Transcribed elements belonged to specific chromosomal loci and were not abundant in the genome. A few elements resulted differentially expressed depending on the treatment. Results suggest that, although most retrotransposons are not expressed, the transcription of such elements is related to their abundance, to their position in the chromosome and to their lineage.

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.

Reporting a Transcript from Iranian Viola Tricolor, Which May Encode a Novel Cyclotide-Like Precursor: Molecular and in silico Studies

The cyclotides are the largest known family of cyclic proteins, which are found in several plant families including Violaceae. They are circular bioactive peptides consisting of 28-37 amino acids, which possess a cyclic cystine knot (CCK) motif and could be useful in biotechnology and drug design as scaffolds for peptide-based drugs. This study describes our finding of a potentially novel gene transcript from the petals of the Iranian Viola tricolor (V. tricolor) flowers. This study is based on the cDNA screening method employed for isolation of cyclotide precursor genes and in silico analysis. Our study resulted in the finding of a novel cyclotide-like precursor from V. tricolor, which is documented in the NCBI by GenBank accession number: KP065812. The in silico analysis revealed that there are lots of similar sequences in many other plant families and they all exhibit some different features from previously discovered cyclotide precursors. The differences occur particularly in the main cyclotide domain that exists without the usual CCK structure. All of these hypothetical precursors have a conserved ER-signal sequence, a Cysteine (C)-rich sequence forming two zinc finger motifs and a cyclotide-like region containing several conserved elements including two highly conserved C residues. In conclusion, using the cDNA screening method we found a potentially new cyclotide-like precursor gene and in silico studies revealed its significant characteristics that may open up a new research line on the distribution and evolution of cyclotides.

Expansion of LINEs and species-specific DNA repeats drives genome expansion in Asian Gypsy Moths

Two subspecies of Asian gypsy moth (AGM), Lymantria dispar asiatica and L. dispar japonica, pose a serious alien invasive threat to North American forests. Despite decades of research on the ecology and biology of this pest, limited AGM-specific genomic resources are currently available. Here, we report on the genome sequences and functional content of these AGM subspecies. The genomes of L.d. asiatica and L.d. japonica are the largest lepidopteran genomes sequenced to date, totaling 921 and 999 megabases, respectively. Large genome size in these subspecies is driven by the accumulation of specific classes of repeats. Genome-wide metabolic pathway reconstructions suggest strong genomic signatures of energy-related pathways in both subspecies, dominated by metabolic functions related to thermogenesis. The genome sequences reported here will provide tools for probing the molecular mechanisms underlying phenotypic traits that are thought to enhance AGM invasiveness.

LTR-TEs abundance, timing and mobility in Solanum commersonii and S. tuberosum genomes following cold-stress conditions

Copia/Ale is the youngest lineage in both Solanum tuberosum and S. commersonii. Within it, we identified nightshade, a new LTR element active in the cultivated potato. From an evolutionary perspective, long-terminal repeat retrotransposons (LTR-RT) activity during stress may be viewed as a mean by which organisms can keep up rates of genetic adaptation to changing conditions. Potato is one of the most important crop consumed worldwide, but studies on LTR-RT characterization are still lacking. Here, we assessed the abundance, insertion time and activity of LTR-RTs in both cultivated Solanum tuberosum and its cold-tolerant wild relative S. commersonii genomes. Gypsy elements were more abundant than Copia ones, suggesting that the former was somehow more successful in colonizing potato genomes. However, Copia elements, and in particular, the Ale lineage, are younger than Gypsy ones, since their insertion time was in average ~ 2 Mya. Due to the ability of LTR-RTs to be circularized by the host DNA repair mechanisms, we identified via mobilome-seq a Copia/Ale element (called nightshade, informal name used for potato family) active in S. tuberosum genome. Our analyses represent a valuable resource for comparative genomics within the Solanaceae, transposon-tagging and for the design of cultivar-specific molecular markers in potato.

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.

Retrotransposons in Plant Genomes: Structure, Identification, and Classification through Bioinformatics and Machine Learning

Transposable elements (TEs) are genomic units able to move within the genome of virtually all organisms. Due to their natural repetitive numbers and their high structural diversity, the identification and classification of TEs remain a challenge in sequenced genomes. Although TEs were initially regarded as "junk DNA", it has been demonstrated that they play key roles in chromosome structures, gene expression, and regulation, as well as adaptation and evolution. A highly reliable annotation of these elements is, therefore, crucial to better understand genome functions and their evolution. To date, much bioinformatics software has been developed to address TE detection and classification processes, but many problematic aspects remain, such as the reliability, precision, and speed of the analyses. Machine learning and deep learning are algorithms that can make automatic predictions and decisions in a wide variety of scientific applications. They have been tested in bioinformatics and, more specifically for TEs, classification with encouraging results. In this review, we will discuss important aspects of TEs, such as their structure, importance in the evolution and architecture of the host, and their current classifications and nomenclatures. We will also address current methods and their limitations in identifying and classifying TEs.

LTR Retrotransposons Show Low Levels of Unequal Recombination and High Rates of Intraelement Gene Conversion in Large Plant Genomes

The accumulation and removal of transposable elements (TEs) is a major driver of genome size evolution in eukaryotes. In plants, long terminal repeat (LTR) retrotransposons (LTR-RTs) represent the majority of TEs and form most of the nuclear DNA in large genomes. Unequal recombination (UR) between LTRs leads to removal of intervening sequence and formation of solo-LTRs. UR is a major mechanism of LTR-RT removal in many angiosperms, but our understanding of LTR-RT-associated recombination within the large, LTR-RT-rich genomes of conifers is quite limited. We employ a novel read-based methodology to estimate the relative rates of LTR-RT-associated UR within the genomes of four conifer and seven angiosperm species. We found the lowest rates of UR in the largest genomes studied, conifers and the angiosperm maize. Recombination may also resolve as gene conversion, which does not remove sequence, so we analyzed LTR-RT-associated gene conversion events (GCEs) in Norway spruce and six angiosperms. Opposite the trend for UR, we found the highest rates of GCEs in Norway spruce and maize. Unlike previous work in angiosperms, we found no evidence that rates of UR correlate with retroelement structural features in the conifers, suggesting that another process is suppressing UR in these species. Recent results from diverse eukaryotes indicate that heterochromatin affects the resolution of recombination, by favoring gene conversion over crossing-over, similar to our observation of opposed rates of UR and GCEs. Control of LTR-RT proliferation via formation of heterochromatin would be a likely step toward large genomes in eukaryotes carrying high LTR-RT content.

A Mechanism for Genome Size Reduction Following Genomic Rearrangements

The factors behind genome size evolution have been of great interest, considering that eukaryotic genomes vary in size by more than three orders of magnitude. Using a model of two wild peanut relatives, Arachis duranensis and Arachis ipaensis, in which one genome experienced large rearrangements, we find that the main determinant in genome size reduction is a set of inversions that occurred in A. duranensis, and subsequent net sequence removal in the inverted regions. We observe a general pattern in which sequence is lost more rapidly at newly distal (telomeric) regions than it is gained at newly proximal (pericentromeric) regions - resulting in net sequence loss in the inverted regions. The major driver of this process is recombination, determined by the chromosomal location. Any type of genomic rearrangement that exposes proximal regions to higher recombination rates can cause genome size reduction by this mechanism. In comparisons between A. duranensis and A. ipaensis, we find that the inversions all occurred in A. duranensis. Sequence loss in those regions was primarily due to removal of transposable elements. Illegitimate recombination is likely the major mechanism responsible for the sequence removal, rather than unequal intrastrand recombination. We also measure the relative rate of genome size reduction in these two Arachis diploids. We also test our model in other plant species and find that it applies in all cases examined, suggesting our model is widely applicable.

Physiological and molecular genetic studies on two elicitors for improving the tolerance of six Egyptian soybean cultivars to cotton leaf worm

Cotton leaf worm (Spodoptera littoralis) is considered one of the most destructive agricultural pests in Egypt. Six soybean cultivars (Giza-21, Giza-22, Giza-35, Giza-82, Giza-83 and Giza-111) were grown under natural infection with cotton leaf worm. The effect of two elicitors, methyl jasmonate and sodium nitroprusside on enhancing the ability of susceptible cultivars to tolerate (Spodoptera littoralis) was studied. Giza-35 and Giza-111 showed tolerance performance under natural infection compared to Giza-22 and Giza-82 as sensitive ones, while Giza-83 and Giza-21 showed moderate tolerance. Both treatments positively affected seed yield and its components and fatty acid composition. Extracted fatty acids showed variable changes in treated plants compared with the untreated controls. Plants treated with the two elicitors showed an increase in Linoleic acid and Linolenic acid fatty acids and decrease in Palmitic acid and Palmitolic acid content. Treatment with methyl jasmonate was found to be more effective than sodium nitroprusside and enhanced resistance of the susceptible cultivars. Eight IRAP and iPBS retrotransposon-based markers were used to detect genetic differences among studied soybean cultivars and to develop molecular genetic markers for cotton leaf worm infestation. The technique successfully identified soybean genotypes in addition to nineteen molecular markers related to soybean tolerance.

Analysis of retrotransposon abundance, diversity and distribution in holocentric Eleocharis (Cyperaceae) genomes

BACKGROUND AND AIMS

Long terminal repeat-retrotransposons (LTR-RTs) comprise a large portion of plant genomes, with massive repeat blocks distributed across the chromosomes. Eleocharis species have holocentric chromosomes, and show a positive correlation between chromosome numbers and the amount of nuclear DNA. To evaluate the role of LTR-RTs in karyotype diversity in members of Eleocharis (subgenus Eleocharis), the occurrence and location of different members of the Copia and Gypsy superfamilies were compared, covering interspecific variations in ploidy levels (considering chromosome numbers), DNA C-values and chromosomal arrangements.

METHODS

The DNA C-value was estimated by flow cytometry. Genomes of Eleocharis elegans and E. geniculata were partially sequenced using Illumina MiSeq assemblies, which were a source for searching for conserved proteins of LTR-RTs. POL domains were used for recognition, comparing families and for probe production, considering different families of Copia and Gypsy superfamilies. Probes were obtained by PCR and used in fluorescence in situ hybridization (FISH) against chromosomes of seven Eleocharis species.

KEY RESULTS

A positive correlation between ploidy levels and the amount of nuclear DNA was observed, but with significant variations between samples with the same ploidy levels, associated with repetitive DNA fractions. LTR-RTs were abundant in E. elegans and E. geniculata genomes, with a predominance of Copia Sirevirus and Gypsy Athila/Tat clades. FISH using LTR-RT probes exhibited scattered and clustered signals, but with differences in the chromosomal locations of Copia and Gypsy. The diversity in LTR-RT locations suggests that there is no typical chromosomal distribution pattern for retrotransposons in holocentric chromosomes, except the CRM family with signals distributed along chromatids.

CONCLUSIONS

These data indicate independent fates for each LTR-RT family, including accumulation between and within chromosomes and genomes. Differential activity and small changes in LTR-RTs suggest a secondary role in nuclear DNA variation, when compared with ploidy changes.

The repetitive landscape of the 5100 Mbp barley genome

Background

While transposable elements (TEs) comprise the bulk of plant genomic DNA, how they contribute to genome structure and organization is still poorly understood. Especially in large genomes where TEs make the majority of genomic DNA, it is still unclear whether TEs target specific chromosomal regions or whether they simply accumulate where they are best tolerated.

Results

Here, we present an analysis of the repetitive fraction of the 5100 Mb barley genome, the largest angiosperm genome to have a near-complete sequence assembly. Genes make only about 2% of the genome, while over 80% is derived from TEs. The TE fraction is composed of at least 350 different families. However, 50% of the genome is comprised of only 15 high-copy TE families, while all other TE families are present in moderate or low copy numbers. We found that the barley genome is highly compartmentalized with different types of TEs occupying different chromosomal “niches”, such as distal, interstitial, or proximal regions of chromosome arms. Furthermore, gene space represents its own distinct genomic compartment that is enriched in small non-autonomous DNA transposons, suggesting that these TEs specifically target promoters and downstream regions. Furthermore, their presence in gene promoters is associated with decreased methylation levels.

Conclusions

Our data show that TEs are major determinants of overall chromosome structure. We hypothesize that many of the the various chromosomal distribution patterns are the result of TE families targeting specific niches, rather than them accumulating where they have the least deleterious effects.

Electronic supplementary material

The online version of this article (10.1186/s13100-017-0102-3) contains supplementary material, which is available to authorized users.

Rapid Increase in Genome Size as a Consequence of Transposable Element Hyperactivity in Wood-White (Leptidea) Butterflies

Characterizing and quantifying genome size variation among organisms and understanding if genome size evolves as a consequence of adaptive or stochastic processes have been long-standing goals in evolutionary biology. Here, we investigate genome size variation and association with transposable elements (TEs) across lepidopteran lineages using a novel genome assembly of the common wood-white (Leptidea sinapis) and population re-sequencing data from both L. sinapis and the closely related L. reali and L. juvernica together with 12 previously available lepidopteran genome assemblies. A phylogenetic analysis confirms established relationships among species, but identifies previously unknown intraspecific structure within Leptidea lineages. The genome assembly of L. sinapis is one of the largest of any lepidopteran taxon so far (643 Mb) and genome size is correlated with abundance of TEs, both in Lepidoptera in general and within Leptidea where L. juvernica from Kazakhstan has considerably larger genome size than any other Leptidea population. Specific TE subclasses have been active in different Lepidoptera lineages with a pronounced expansion of predominantly LINEs, DNA elements, and unclassified TEs in the Leptidea lineage after the split from other Pieridae. The rate of genome expansion in Leptidea in general has been in the range of four Mb/Million year (My), with an increase in a particular L. juvernica population to 72 Mb/My. The considerable differences in accumulation rates of specific TE classes in different lineages indicate that TE activity plays a major role in genome size evolution in butterflies and moths.

Different histories of two highly variable LTR retrotransposons in sunflower species

In the Helianthus genus, very large intra- and interspecific variability related to two specific retrotransposons of Helianthus annuus (Helicopia and SURE) exists. When comparing these two sequences to sunflower sequence databases recently produced by our lab, the Helicopia family was shown to belong to the Maximus/SIRE lineage of the Sirevirus genus of the Copia superfamily, whereas the SURE element (whose superfamily was not even previously identified) was classified as a Gypsy element of the Ogre/Tat lineage of the Metavirus genus. Bioinformatic analysis of the two retrotransposon families revealed their genomic abundance and relative proliferation timing. The genomic abundance of these families differed significantly among 12 Helianthus species. The ratio between the abundance of long terminal repeats and their reverse transcriptases suggested that the SURE family has relatively more solo long terminal repeats than does Helicopia. Pairwise comparisons of Illumina reads encoding the reverse transcriptase domain indicated that SURE amplification may have occurred more recently than that of Helicopia. Finally, the analysis of population structure based on the SURE and Helicopia polymorphisms of 32 Helianthus species evidenced two subpopulations, which roughly corresponded to species of the Helianthus and Divaricati/Ciliares sections. However, a number of species showed an admixed structure, confirming the importance of interspecific hybridisation in the evolution of this genus. In general, these two retrotransposon families differentially contributed to interspecific variability, emphasising the need to refer to specific families when studying genome evolution.

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.

Useful parasites: the evolutionary biology and biotechnology applications of transposable elements

Transposable elements usually comprise the most abundant nongenic fraction of eukaryotic genomes. Because of their capacity to selfreplicate and to induce a wide range of mutations, transposable elements have long been considered as 'parasitic' or 'selfish'. Today, we recognize that the findings about genomic changes affected by transposable elements have considerably altered our view of the ways in which genomes evolve and work. Numerous studies have provided evidences that mobile elements have the potential to act as agents of evolution by increasing, rearranging and diversifying the genetic repertoire of their hosts. With large-scale sequencing becoming increasingly available, more and more scientists come across transposable element sequences in their data. I will provide examples that transposable elements, although having signatures of 'selfish' DNA, play a significant biological role in the maintainance of genome integrity and providing novel regulatoty networks. These features, along with the transpositional and mutagenic capacity to produce a raw genetic diversity, make the genome mobile fraction, a key player in species adaptation and microevolution. The last but not least, transposable elements stand as informative DNA markers that may complement other conventional DNA markers. Altogether, transposable elements represent a promising, but still largely unexplored research niche and deserve to be included into the agenda of molecular ecologists, evolutionary geneticists, conservation biologists and plant breeders.

Patterns of Evolutionary Trajectories and Domestication History within the Genus Hordeum Assessed by REMAP Markers

The patterns of genetic diversity related to the taxonomy and domestication history of 85 accessions representing the main four species of the genus Hordeum were examined by retrotransposon-microsatellite amplified polymorphism (REMAP) markers based on the retrotransposon BARE-1. A substantial level of genetic polymorphisms at among- and within-species level was observed showing that this retrotransposon family and its adjacent genomic regions has been a target for genome dynamics during the evolution and domestication of barley. The obtained data are consistent with the current taxonomic status within the genus Hordeum. Similar level of genetic diversity was observed between the wild and the domesticated barley accessions suggesting that transposable elements` activity and accumulation may counteract the decrease of genome-wide diversity following domestication. In addition, eco-geographical sub-genome pools of the cultivated barley were identified in support to the theory of multiple origins of domestication within the genus Hordeum. We also provide conclusions about the relationship between accessions of different species and the putative routes of barley domestication. In conclusion, the retrotransposon BARE-1 stands as a reliable and perspective DNA marker for the assessment of the phylogenetic and domestication history in the genus Hordeum and other crop species.

Dynamics of genome size evolution in birds and mammals

Genome size in mammals and birds shows remarkably little interspecific variation compared with other taxa. However, genome sequencing has revealed that many mammal and bird lineages have experienced differential rates of transposable element (TE) accumulation, which would be predicted to cause substantial variation in genome size between species. Thus, we hypothesize that there has been covariation between the amount of DNA gained by transposition and lost by deletion during mammal and avian evolution, resulting in genome size equilibrium. To test this model, we develop computational methods to quantify the amount of DNA gained by TE expansion and lost by deletion over the last 100 My in the lineages of 10 species of eutherian mammals and 24 species of birds. The results reveal extensive variation in the amount of DNA gained via lineage-specific transposition, but that DNA loss counteracted this expansion to various extents across lineages. Our analysis of the rate and size spectrum of deletion events implies that DNA removal in both mammals and birds has proceeded mostly through large segmental deletions (>10 kb). These findings support a unified "accordion" model of genome size evolution in eukaryotes whereby DNA loss counteracting TE expansion is a major determinant of genome size. Furthermore, we propose that extensive DNA loss, and not necessarily a dearth of TE activity, has been the primary force maintaining the greater genomic compaction of flying birds and bats relative to their flightless relatives.

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.

Transposable elements and G-quadruplexes

A significant part of eukaryotic genomes is formed by transposable elements (TEs) containing not only genes but also regulatory sequences. Some of the regulatory sequences located within TEs can form secondary structures like hairpins or three-stranded (triplex DNA) and four-stranded (quadruplex DNA) conformations. This review focuses on recent evidence showing that G-quadruplex-forming sequences in particular are often present in specific parts of TEs in plants and humans. We discuss the potential role of these structures in the TE life cycle as well as the impact of G-quadruplexes on replication, transcription, translation, chromatin status, and recombination. The aim of this review is to emphasize that TEs may serve as vehicles for the genomic spread of G-quadruplexes. These non-canonical DNA structures and their conformational switches may constitute another regulatory system that, together with small and long non-coding RNA molecules and proteins, contribute to the complex cellular network resulting in the large diversity of eukaryotes.

The diversification of PHIS transposon superfamily in eukaryotes

BACKGROUND

PHIS transposon superfamily belongs to DNA transposons and includes PIF/Harbinger, ISL2EU, and Spy transposon groups. These three groups have similar DDE domain-containing transposases; however, their coding capacity, species distribution, and target site duplications (TSDs) are significantly different.

RESULTS

In this study, we systematically identified and analyzed PHIS transposons in 836 sequenced eukaryotic genomes using transposase homology search and structure approach. In total, 380 PHIS families were identified in 112 genomes and 168 of 380 families were firstly reported in this study. Besides previous identified PIF/Harbinger, ISL2EU, and Spy groups, three new types (called Pangu, NuwaI, and NuwaII) of PHIS superfamily were identified; each has its own distinctive characteristics, especially in TSDs. Pangu and NuwaII transposons are characterized by 5'-ANT-3' and 5'-C|TNA|G-3' TSDs, respectively. Both transposons are widely distributed in plants, fungi, and animals; the NuwaI transposons are characterized by 5'-CWG-3' TSDs and mainly distributed in animals.

CONCLUSIONS

Here, in total, 380 PHIS families were identified in eukaryotes. Among these 380 families, 168 were firstly reported in this study. Furthermore, three new types of PHIS superfamily were identified. Our results not only enrich the transposon diversity but also have extensive significance for improving genome sequence assembly and annotation of higher organisms.

The Ty1-copia LTR retroelement family PARTC is highly conserved in conifers over 200 MY of evolution

Long Terminal Repeat retroelements (LTR-RTs) are a major component of many plant genomes. Although well studied and described in angiosperms, their features and dynamics are poorly understood in gymnosperms. Representative complete copies of a Ty1-copia element isolate in Picea abies and named PARTC were identified in six other conifer species (Picea glauca, Pinus sylvestris, Pinus taeda, Abies sibirica, Taxus baccata and Juniperus communis) covering more than 200 million years of evolution. Here we characterized the structure of this element, assessed its abundance across conifers, studied the modes and timing of its amplification, and evaluated the degree of conservation of its extant copies at nucleotide level over distant species. We demonstrated that the element is ancient, abundant, widespread and its paralogous copies are present in the genera Picea, Pinus and Abies as an LTR-RT family. The amplification leading to the extant copies of PARTC occurred over long evolutionary times spanning 10s of MY and mostly took place after the speciation of the conifers analyzed. The level of conservation of PARTC is striking and may be explained by low substitution rates and limited removal mechanisms for LTR-RTs. These PARTC features and dynamics are representative of a more general scenario for LTR-RTs in gymnosperms quite different from that characterizing the vast majority of LTR-RT elements in angiosperms.

Transposable Element Insertion and Epigenetic Modification Cause the Multiallelic Variation in the Expression of FAE1 in Sinapis alba

Naturally occurring heritable variation provides a fundamental resource to reveal the genetic and molecular bases of traits in forward genetic studies. Here, we report the molecular basis of the differences in the four alleles E¹, E², E³, and e of the FATTY ACID ELONGATION1 (FAE1) gene controlling high, medium, low, and zero erucic content in yellow mustard (Sinapis alba). E¹ represents a fully functional allele with a coding DNA sequence (CDS) of 1521 bp and a promoter adjacent to the CDS. The null allele e resulted from an insertional disruption in the CDS by Sal-PIF, a 3100-bp PIF/Harbinger-like DNA transposon, whereas E² and E³ originated from the insertion of Sal-T1, a 4863-bp Copia-like retrotransposon, in the 5' untranslated region. E³ was identical to E² but showed cytosine methylation in the promoter region and was thus an epiallele having a further reduction in expression. The coding regions of E² and E³ also contained five single-nucleotide polymorphisms (SNPs) not present in E¹, but expression studies in Saccharomyces cerevisiae indicated that these SNPs did not affect enzyme functionality. These results demonstrate a comprehensive molecular framework for the interplay of transposon insertion, SNP/indel mutation, and epigenetic modification influencing the broad range of natural genetic variation in plants.

Environmental stress activation of plant long-terminal repeat retrotransposons

Genomic retrotransposons (RTs) are major components of most plant genomes. They spread throughout the genomes by a process termed retrotransposition, which consists of reverse transcription and reinsertion of the copied element into a new genomic location (a copy-and-paste system). Abiotic and biotic stresses activate long-terminal repeat (LTR) RTs in photosynthetic eukaryotes from algae to angiosperms. LTR RTs could represent a threat to the integrity of host genomes because of their activity and mutagenic potential by epigenetic regulation. Host genomes have developed mechanisms to control the activity of the retroelements and their mutagenic potential. Some LTR RTs escape these defense mechanisms, and maintain their ability to be activated and transpose as a result of biotic or abiotic stress stimuli. These stimuli include pathogen infection, mechanical damage, in vitro tissue culturing, heat, drought and salt stress, generation of doubled haploids, X-ray irradiation and many others. Reactivation of LTR RTs differs between different plant genomes. The expression levels of reactivated RTs are influenced by the transcriptional and post-transcriptional gene silencing mechanisms (e.g. DNA methylation, heterochromatin formation and RNA interference). Moreover, the insertion of RTs (e.g. Triticum aestivum L. Wis2-1A) into or next to coding regions of the host genome can generate changes in the expression of adjacent host genes of the host. In this paper, we review the ways that plant genomic LTR RTs are activated by environmental stimuli to affect restructuring and diversification of the host genome.

Contrasting evolutionary trajectories of multiple retrotransposons following independent allopolyploidy in wild wheats

Transposable elements (TEs) are expectedly central to genome evolution. To assess the impact of TEs in driving genome turnover, we used allopolyploid genomes, showing considerable deviation from the predicted additivity of their diploid progenitors and thus having undergone major restructuring. Genome survey sequencing was used to select 17 putatively active families of long terminal repeat retrotransposons. Genome-wide TE insertions were genotyped with sequence-specific amplified polymorphism (SSAP) in diploid progenitors and their derived polyploids, and compared with changes in random sequences to assess restructuring of four independent Aegilops allotetraploid genomes. Generally, TEs with different evolutionary trajectories from those of random sequences were identified. Thus, TEs presented family-specific and species-specific dynamics following polyploidy, as illustrated by Sabine showing proliferation in particular polyploids, but massive elimination in others. Contrasting with that, only a few families (BARE1 and Romani) showed proliferation in all polyploids. Overall, TE divergence between progenitors was strongly correlated with the degree of restructuring in polyploid TE fractions. TE families present evolutionary trajectories that are decoupled from genome-wide changes after allopolyploidy and have a pervasive impact on their restructuring.

Major repeat components covering one-third of the ginseng (Panax ginseng C.A. Meyer) genome and evidence for allotetraploidy

Ginseng (Panax ginseng) is a famous medicinal herb, but the composition and structure of its genome are largely unknown. Here we characterized the major repeat components and inspected their distribution in the ginseng genome. By analyzing three repeat-rich bacterial artificial chromosome (BAC) sequences from ginseng, we identified complex insertion patterns of 34 long terminal repeat retrotransposons (LTR-RTs) and 11 LTR-RT derivatives accounting for more than 80% of the BAC sequences. The LTR-RTs were classified into three Ty3/gypsy (PgDel, PgTat and PgAthila) and two Ty1/Copia (PgTork and PgOryco) families. Mapping of 30-Gbp Illumina whole-genome shotgun reads to the BAC sequences revealed that these five LTR-RT families occupy at least 34% of the ginseng genome. The Ty3/Gypsy families were predominant, comprising 74 and 33% of the BAC sequences and the genome, respectively. In particular, the PgDel family accounted for 29% of the genome and presumably played major roles in enlargement of the size of the ginseng genome. Fluorescence in situ hybridization (FISH) revealed that the PgDel1 elements are distributed throughout the chromosomes along dispersed heterochromatic regions except for ribosomal DNA blocks. The intensity of the PgDel2 FISH signals was biased toward 24 out of 48 chromosomes. Unique gene probes showed two pairs of signals with different locations, one pair in subtelomeric regions on PgDel2-rich chromosomes and the other in interstitial regions on PgDel2-poor chromosomes, demonstrating allotetraploidy in ginseng. Our findings promote understanding of the evolution of the ginseng genome and of that of related species in the Araliaceae.

Quadruplex-forming DNA sequences spread by retrotransposons may serve as genome regulators

Transposable elements (TEs) are ubiquitous genome inhabitants in eukaryotes. Increasing evidence shows that TEs are involved in regulatory networks of eukaryotic cells and contribute to genome evolution. Recently, we reported that many plant long-terminal repeat (LTR) retrotransposons contain DNA quadruplex-forming sequences at precise positions inside their LTRs and that quadruplexes are better preserved in evolutionary younger elements. As quadruplexes can modulate molecular processes, quadruplexes found at specific distances upstream and downstream from the endogenous TE promoter can affect transcription of the element. Moreover, quadruplexes found in solo LTRs, as well as in 3' ends of 5'-truncated copies of LINE-1 elements, can affect expression of neighboring genes. Here, we propose that this way retrotransposons can serve as vehicles for spread of DNA quadruplexes. Quadruplexes can thus fulfill a dual regulatory role-to influence both the retrotransposons carrying them and the neighboring host genes, e.g., by direct effect on transcription or by modifying the local chromatin state. Additionally, four-stranded DNA structures may serve as hotspots for recombination-based genome rearrangements.

Retrotransposon BARE displays strong tissue-specific differences in expression

The BARE retrotransposon comprises c. 10% of the barley (Hordeum vulgare) genome. It is actively transcribed, translated and forms virus-like particles (VLPs). For retrotransposons, the inheritance of new copies depends critically on where in the plant replication occurs. In order to shed light on the replication strategy of BARE in the plant, we have used immunolocalization and in situ hybridization to examine expression of the BARE capsid protein, Gag, at a tissue-specific level. Gag is expressed in provascular tissues and highly localized in companion cells surrounding the phloem sieve tubes in mature vascular tissues. BARE Gag and RNA was not seen in the shoot apical meristem of young seedlings, but appeared, following transition to flowering, in the developing floral spike. Moreover, Gag has a highly specific localization in pre-fertilization ovaries. The strong presence of Gag in the floral meristems suggests that newly replicated copies there will be passed to the next generation. BARE expression patterns are consistent with transcriptional regulation by predicted response elements in the BARE promoter, and in the ovary with release from epigenetic transcriptional silencing. To our knowledge, this is the first analysis of the expression of native retrotransposon proteins within a plant to be reported.

Transposon fingerprinting using low coverage whole genome shotgun sequencing in cacao (Theobroma cacao L.) and related species

Background

Transposable elements (TEs) and other repetitive elements are a large and dynamically evolving part of eukaryotic genomes, especially in plants where they can account for a significant proportion of genome size. Their dynamic nature gives them the potential for use in identifying and characterizing crop germplasm. However, their repetitive nature makes them challenging to study using conventional methods of molecular biology. Next generation sequencing and new computational tools have greatly facilitated the investigation of TE variation within species and among closely related species.

Results

(i) We generated low-coverage Illumina whole genome shotgun sequencing reads for multiple individuals of cacao (Theobroma cacao) and related species. These reads were analysed using both an alignment/mapping approach and a de novo (graph based clustering) approach. (ii) A standard set of ultra-conserved orthologous sequences (UCOS) standardized TE data between samples and provided phylogenetic information on the relatedness of samples. (iii) The mapping approach proved highly effective within the reference species but underestimated TE abundance in interspecific comparisons relative to the de novo methods. (iv) Individual T. cacao accessions have unique patterns of TE abundance indicating that the TE composition of the genome is evolving actively within this species. (v) LTR/Gypsy elements are the most abundant, comprising c.10% of the genome. (vi) Within T. cacao the retroelement families show an order of magnitude greater sequence variability than the DNA transposon families. (vii) Theobroma grandiflorum has a similar TE composition to T. cacao, but the related genus Herrania is rather different, with LTRs making up a lower proportion of the genome, perhaps because of a massive presence (c. 20%) of distinctive low complexity satellite-like repeats in this genome.

Conclusions

(i) Short read alignment/mapping to reference TE contigs provides a simple and effective method of investigating intraspecific differences in TE composition. It is not appropriate for comparing repetitive elements across the species boundaries, for which de novo methods are more appropriate. (ii) Individual T. cacao accessions have unique spectra of TE composition indicating active evolution of TE abundance within this species. TE patterns could potentially be used as a “fingerprint” to identify and characterize cacao accessions.