LTR retrotransposon dynamics in the evolution of the olive (Olea europaea) genome

PlantLTRdb: An interactive database for 195 plant species LTR-retrotransposons

LTR-retrotransposons (LTR-RTs) are a large group of transposable elements that replicate through an RNA intermediate and alter genome structure. The activities of LTR-RTs in plant genomes provide helpful information about genome evolution and gene function. LTR-RTs near or within genes can directly alter gene function. This work introduces PlantLTRdb, an intact LTR-RT database for 195 plant species. Using homology- and de novo structure-based methods, a total of 150.18 Gbp representing 3,079,469 pseudomolecules/scaffolds were analyzed to identify, characterize, annotate LTR-RTs, estimate insertion ages, detect LTR-RT-gene chimeras, and determine nearby genes. Accordingly, 520,194 intact LTR-RTs were discovered, including 29,462 autonomous and 490,732 nonautonomous LTR-RTs. The autonomous LTR-RTs included 10,286 Gypsy and 19,176 Copia, while the nonautonomous were divided into 224,906 Gypsy, 218,414 Copia, 1,768 BARE-2, 3,147 TR-GAG and 4,2497 unknown. Analysis of the identified LTR-RTs located within genes showed that a total of 36,236 LTR-RTs were LTR-RT-gene chimeras and 11,619 LTR-RTs were within pseudo-genes. In addition, 50,026 genes are within 1 kbp of LTR-RTs, and 250,587 had a distance of 1 to 10 kbp from LTR-RTs. PlantLTRdb allows researchers to search, visualize, BLAST and analyze plant LTR-RTs. PlantLTRdb can contribute to the understanding of structural variations, genome organization, functional genomics, and the development of LTR-RT target markers for molecular plant breeding. PlantLTRdb is available at https://bioinformatics.um6p.ma/PlantLTRdb.

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.

Understanding the Role of PIN Auxin Carrier Genes under Biotic and Abiotic Stresses in Olea europaea L

The PIN-FORMED (PIN) proteins represent the most important polar auxin transporters in plants. Here, we characterized the PIN gene family in two olive genotypes, the Olea europaea subsp. europaea var. sylvestris and the var. europaea (cv. ‘Farga’). Twelve and 17 PIN genes were identified for vars. sylvestris and europaea, respectively, being distributed across 6 subfamilies. Genes encoding canonical OePINs consist of six exons, while genes encoding non-canonical OePINs are composed of five exons, with implications at protein specificities and functionality. A copia-LTR retrotransposon located in intron 4 of OePIN2b of var. europaea and the exaptation of partial sequences of that element as exons of the OePIN2b of var. sylvestris reveals such kind of event as a driving force in the olive PIN evolution. RNA-seq data showed that members from the subfamilies 1, 2, and 3 responded to abiotic and biotic stress factors. Co-expression of OePINs with genes involved in stress signaling and oxidative stress homeostasis were identified. This study highlights the importance of PIN genes on stress responses, contributing for a holistic understanding of the role of auxins in plants.

The Singular Evolution of Olea Genome Structure

The current view of plant genome evolution proposes that genome size has mainly been determined by polyploidisation and amplification/loss of transposons, with a minor role played by other repeated sequences, such as tandem repeats. In cultivated olive (Olea europaea subsp. europaea var. europaea), available data suggest a singular model of genome evolution, in which a massive expansion of tandem-repeated sequences accompanied changes in nuclear architecture. This peculiar scenario highlights the importance of focusing on Olea genus evolution, to shed light on mechanisms that led to its present genomic structure. Next-generation sequencing technologies, bioinformatics and in situ hybridisation were applied to study the genomic structure of five related Olea taxa, which originated at different times from their last common ancestor. On average, repetitive DNA in the Olea taxa ranged from ~59% to ~73% of the total genome, showing remarkable differences in terms of composition. Among repeats, we identified 11 major families of tandem repeats, with different abundances in the analysed taxa, five of which were novel discoveries. Interestingly, overall tandem repeat abundance was inversely correlated to that of retrotransposons. This trend might imply a competition in the proliferation of these repeat classes. Indeed, O. paniculata, the species closest to the Olea common ancestor, showed very few tandem-repeated sequences, while it was rich in long terminal repeat retrotransposons, suggesting that the amplification of tandem repeats occurred after its divergence from the Olea ancestor. Furthermore, some tandem repeats were physically localised in closely related O. europaea subspecies (i.e., cultivated olive and O. europaea subsp. cuspidata), which showed a significant difference in tandem repeats abundance. For 4 tandem repeats families, a similar number of hybridisation signals were observed in both subspecies, apparently indicating that, after their dissemination throughout the olive genome, these tandem repeats families differentially amplified maintaining the same positions in each genome. Overall, our research identified the temporal dynamics shaping genome structure during Olea speciation, which represented a singular model of genome evolution in higher plants.

Genome-Wide Identification and Development of LTR Retrotransposon-Based Molecular Markers for the Melilotus Genus

Melilotus is an important genus of legumes with industrial and medicinal value, partly due to the production of coumarin. To explore the genetic diversity and population structure of Melilotus, 40 accessions were analyzed using long terminal repeat (LTR) retrotransposon-based markers. A total of 585,894,349 bp of LTR retrotransposon sequences, accounting for 55.28% of the Melilotus genome, were identified using bioinformatics tools. A total of 181,040 LTR retrotransposons were identified and classified as Gypsy, Copia, or another type. A total of 350 pairs of primers were designed for assessing polymorphisms in 15 Melilotus albus accessions. Overall, 47 polymorphic primer pairs were screened for their availability and transferability in 18 Melilotus species. All the primer pairs were transferable, and 292 alleles were detected at 47 LTR retrotransposon loci. The average polymorphism information content (PIC) value was 0.66, which indicated that these markers were highly informative. Based on unweighted pair group method with arithmetic mean (UPGMA) dendrogram cluster analysis, the 18 Melilotus species were classified into three clusters. This study provides important data for future breeding programs and for implementing genetic improvements in the Melilotus genus.

Genome-wide analysis of HSP90 gene family in the Mediterranean olive (Olea europaea subsp. europaea) provides insight into structural patterns, evolution and functional diversity

Plants regularly experience multiple abiotic and biotic pressures affecting their normal development. The 90-kDa heat shock protein (HSP90) plays a dynamic role in countering abiotic and biotic stresses via a plethora of functional mechanisms. The HSP90 has been investigated in many plant species. However, there is little information available about this gene family in the cultivated Mediterranean olive tree, Olea europaea subsp. europaea var. europaea. In the current study, we systematically performed genome-wide identification and characterization of the HSP90 gene family in O. europaea var. europaea (OeHSP90s). Twelve regular OeHSP90s were identified, which were phylogenetically grouped into two major clusters and four sub-clusters, showing five paralogous gene pairs evolving under purifying selection. All of the 12 proteins contained a Histidine kinase-like ATPase (HATPase_c) domain, justifying the role played by HSP90 proteins in ATP binding and hydrolysis. The predicted 3D structure of OeHSP90 proteins provided information to understand their functions at the biochemical level. Consistent with their phylogenetic relationships, OeHSP90 members were predicted to be localized in different cellular compartments, suggesting their involvement in various subcellular processes. In consonance with their spatial organization, olive HSP90 family members were found to share similar motif arrangements and similar number of exons. We found that OeHSP90 promoters contained various cis-acting elements associated with light responsiveness, hormone signaling pathways and reaction to various stress conditions. In addition, expression sequence tags (ESTs) analysis offered a view of OeHSP90 tissue- and developmental stage specific pattern of expression. Proteins interacting with OeHSP90s were predicted and their potential roles were discussed. Overall, our results offer premises for further investigation of the implication of HSP90 genes in the physiological processes of the olive and its adaptation to stresses.

Genome wide annotation and characterization of young, intact long terminal repeat retrotransposons (In-LTR-RTs) of seven legume species

Availability of genome sequence of different legume species has provided an opportunity to characterize the abundance, distribution, and divergence of canonical intact long terminal retrotransposons (In-LTR-RT) superfamilies. Among seven legume species, Arachis ipaensis (Aip) showed the highest number of full-length canonical In-LTR-RTs (3325), followed by Glycine max (Gma, 2328), Vigna angularis (Van, 1625), Arachis durensis (Adu, 1348), Lotus japonicus (Lja, 1294), Medicago truncatula (Mtr, 788), and Circer arietinum (Car, 124). Divergence time analysis demonstrated that the amplification timeframe of LTR-RTs dramatically varied in different families. The average insertion time of Copia element varied from 0.51 (Van) to 1.37 million years ago (Mya) (Adu, and Aip), whereas that of Gypsy was between 0.22 (Mtr) and 1.82 Mya (Adu). Bayesian phylogenetic tree analysis suggested that the 1397 and 1917 reverse transcriptase (RT) domains of Copia and Gypsy families of the seven legume species were clustered into 7 and 14 major groups, respectively. The highest proportion (approximately 94.79-100%) of transposable element (TE)-associated genes assigned to pathways was mapped to metabolism-related pathways in all species. The results enabled the structural understanding of full-length In-LTR-RTs and will be valuable resource for the further study of the impact of TEs on gene structure and expression in legume species.

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.

Physiological, epigenetic and genetic regulation in some olive cultivars under salt stress

Cultivated olive, a typical fruit crop species of the semi-arid regions, could successfully face the new scenarios driven by the climate change through the selection of tolerant varieties to salt and drought stresses. In the present work, multidisciplinary approaches, including physiological, epigenetic and genetic studies, have been applied to clarify the salt tolerance mechanisms in olive. Four varieties (Koroneiki, Royal de Cazorla, Arbequina and Picual) and a related form (O. europaea subsp. cuspidata) were grown in a hydroponic system under different salt concentrations from zero to 200 mM. In order to verify the plant response under salt stress, photosynthesis, gas exchange and relative water content were measured at different time points, whereas chlorophyll and leaf concentration of Na+, K+ and Ca2+ ions, were quantified at 43 and 60 days after treatment, when stress symptoms became prominent. Methylation sensitive amplification polymorphism (MSAP) technique was used to assess the effects of salt stress on plant DNA methylation. Several fragments resulted differentially methylated among genotypes, treatments and time points. Real time quantitative PCR (RT-qPCR) analysis revealed significant expression changes related to plant response to salinity. Four genes (OePIP1.1, OePetD, OePI4Kg4 and OeXyla) were identified, as well as multiple retrotransposon elements usually targeted by methylation under stress conditions.

The landscape of transposable elements and satellite DNAs in the genome of a dioecious plant spinach (Spinacia oleracea L.)

Background

Repetitive sequences, including transposable elements (TEs) and satellite DNAs, occupy a considerable portion of plant genomes. Analysis of the repeat fraction benefits the understanding of genome structure and evolution. Spinach (Spinacia oleracea L.), an important vegetable crop, is also a model dioecious plant species for studying sex determination and sex chromosome evolution. However, the repetitive sequences of the spinach genome have not been fully investigated.

Results

We extensively analyzed the repetitive components of draft spinach genome, especially TEs and satellites, by different strategies. A total of 16,002 full-length TEs were identified. Among the most abundant long terminal repeat (LTR) retrotransposons (REs), Copia elements were overrepresented compared with Gypsy ones. Angela was the most dominating Copia lineage; Ogre/Tat was the most abundant Gypsy lineage. The mean insertion age of LTR-REs was 1.42 million years; approximately 83.7% of these elements were retrotransposed during the last two million years. RepeatMasker totally masked about 64.05% of the spinach genome, with LTR-REs, non-LTR-REs, and DNA transposons occupying 49.2, 2.4, and 5.6%, respectively. Fluorescence in situ hybridization (FISH) analysis showed that most LTR-REs dispersed all over the chromosomes, by contrast, elements of CRM lineage were distributed at the centromeric region of all chromosomes. In addition, Ogre/Tat lineage mainly accumulated on sex chromosomes, and satellites Spsat2 and Spsat3 were exclusively located at the telomeric region of the short arm of sex chromosomes.

Conclusions

We reliably annotated the TE fraction of the draft genome of spinach. FISH analysis indicates that Ogre/Tat lineage and the sex chromosome-specific satellites DNAs might participate in sex chromosome formation and evolution. Based on FISH signals of microsatellites, together with 45S rDNA, a fine karyotype of spinach was established. This study improves our knowledge of repetitive sequence organization in spinach genome and aids in accurate spinach karyotype construction.

Electronic supplementary material

The online version of this article (10.1186/s13100-019-0147-6) contains supplementary material, which is available to authorized users.

Specific LTR-Retrotransposons Show Copy Number Variations between Wild and Cultivated Sunflowers

The relationship between variation of the repetitive component of the genome and domestication in plant species is not fully understood. In previous work, variations in the abundance and proximity to genes of long terminal repeats (LTR)-retrotransposons of sunflower (Helianthus annuus L.) were investigated by Illumina DNA sequencingtocompare cultivars and wild accessions. In this study, we annotated and characterized 22 specific retrotransposon families whose abundance varies between domesticated and wild genotypes. These families mostly belonged to the Chromovirus lineage of the Gypsy superfamily and were distributed overall chromosomes. They were also analyzed in respect to their proximity to genes. Genes close to retrotransposon were classified according to biochemical pathways, and differences between domesticated and wild genotypes are shown. These data suggest that structural variations related to retrotransposons might have occurred to produce phenotypic variation between wild and domesticated genotypes, possibly by affecting the expression of genes that lie close to inserted or deleted retrotransposons and belong to specific biochemical pathways as those involved in plant stress responses.

Intact DNA purified from flow-sorted nuclei unlocks the potential of next-generation genome mapping and assembly in Solanum species

Next-generation genome mapping through nanochannels (Bionano optical mapping) of plant genomes brings genome assemblies to the ‘nearly-finished’ level for reliable and detailed gene annotations and assessment of structural variations. Despite the recent progress in its development, researchers face the technical challenges of obtaining sufficient high molecular weight (HMW) nuclear DNA due to cell walls which are difficult to disrupt and to the presence of cytoplasmic polyphenols and polysaccharides that co-precipitate or are covalently bound to DNA and might cause oxidation and/or affect the access of nicking enzymes to DNA, preventing downstream applications. Here we describe important improvements for obtaining HMW DNA that we tested on Solanum crops and wild relatives. The methods that we further elaborated and refined focus on

•

Improving flexibility of using different tissues as source materials, like fast-growing root tips and young leaves from seedlings or in vitro plantlets.

•

Obtaining nuclei suspensions through either lab homogenizers or by chopping.

•

Increasing flow sorting efficiency using DAPI (4′,6-diamidino-2-phenylindole) and PI (propidium iodide) DNA stains, with different lasers (UV or 488 nm) and sorting platforms such as the FACSAria and FACSVantage flow sorters, thus making it appropriate for more laboratories working on plant genomics.

The obtained nuclei are embedded into agarose plugs for processing and isolating uncontaminated HMW DNA, which is a prerequisite for nanochannel-based next-generation optical mapping strategies.

AOX1-Subfamily Gene Members in Olea europaea cv. "Galega Vulgar"-Gene Characterization and Expression of Transcripts during IBA-Induced in Vitro Adventitious Rooting

Propagation of some Olea europaea L. cultivars is strongly limited due to recalcitrant behavior in adventitious root formation by semi-hardwood cuttings. One example is the cultivar ”Galega vulgar”. The formation of adventitious roots is considered a morphological response to stress. Alternative oxidase (AOX) is the terminal oxidase of the alternative pathway of the plant mitochondrial electron transport chain. This enzyme is well known to be induced in response to several biotic and abiotic stress situations. This work aimed to characterize the alternative oxidase 1 (AOX1)-subfamily in olive and to analyze the expression of transcripts during the indole-3-butyric acid (IBA)-induced in vitro adventitious rooting (AR) process. OeAOX1a (acc. no. MF410318) and OeAOX1d (acc. no. MF410319) were identified, as well as different transcript variants for both genes which resulted from alternative polyadenylation events. A correlation between transcript accumulation of both OeAOX1a and OeAOX1d transcripts and the three distinct phases (induction, initiation, and expression) of the AR process in olive was observed. Olive AOX1 genes seem to be associated with the induction and development of adventitious roots in IBA-treated explants. A better understanding of the molecular mechanisms underlying the stimulus needed for the induction of adventitious roots may help to develop more targeted and effective rooting induction protocols in order to improve the rooting ability of difficult-to-root cultivars.

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.

Recent developments in olive (Olea europaea L.) genetics and genomics: applications in taxonomy, varietal identification, traceability and breeding

The latest results in DNA markers application and genomic studies in olive. Olive (Olea europaea L.) is among the most ancient tree crops worldwide and the source of oil beneficial for human health. Despite this, few data on olive genetics are available in comparison with other cultivated plant species. Molecular information is mainly linked to molecular markers and their application to the study of DNA variation in the Olea europaea complex. In terms of genomic research, efforts have been made in sequencing, heralding the era of olive genomic. The present paper represents an update of a previous review work published in this journal in 2011. The review is again mainly focused on DNA markers, whose application still constitutes a relevant percentage of the most recently published researches. Since the olive genomic era has recently started, the latest results in this field are also being discussed.

Genome-wide analysis of LTR-retrotransposon diversity and its impact on the evolution of the genus Helianthus (L.)

Background

Genome divergence by mobile elements activity and recombination is a continuous process that plays a key role in the evolution of species. Nevertheless, knowledge on retrotransposon-related variability among species belonging to the same genus is still limited. Considering the importance of the genus Helianthus, a model system for studying the ecological genetics of speciation and adaptation, we performed a comparative analysis of the repetitive genome fraction across ten species and one subspecies of sunflower, focusing on long terminal repeat retrotransposons at superfamily, lineage and sublineage levels.

Results

After determining the relative genome size of each species, genomic DNA was isolated and subjected to Illumina sequencing. Then, different assembling and clustering approaches allowed exploring the repetitive component of all genomes. On average, repetitive DNA in Helianthus species represented more than 75% of the genome, being composed mostly by long terminal repeat retrotransposons. Also, the prevalence of Gypsy over Copia superfamily was observed and, among lineages, Chromovirus was by far the most represented. Although nearly all the same sublineages are present in all species, we found considerable variability in the abundance of diverse retrotransposon lineages and sublineages, especially between annual and perennial species.

Conclusions

This large variability should indicate that different events of amplification or loss related to these elements occurred following species separation and should have been involved in species differentiation. Our data allowed us inferring on the extent of interspecific repetitive DNA variation related to LTR-RE abundance, investigating the relationship between changes of LTR-RE abundance and the evolution of the genus, and determining the degree of coevolution of different LTR-RE lineages or sublineages between and within species. Moreover, the data suggested that LTR-RE abundance in a species was affected by the annual or perennial habit of that species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-4050-6) contains supplementary material, which is available to authorized users.

Recessive male sterility in cabbage (Brassica oleracea var. capitata) caused by loss of function of BoCYP704B1 due to the insertion of a LTR-retrotransposon

The LTR-retrotransposon insertion in BoCYP704B1 is proved to be the primary cause of the male sterility in cabbage. Effective allele-specific markers were developed for marker-assisted selection of male sterile gene. 83121A is a spontaneous male sterile mutant identified from cabbage. Genetic analysis indicated that male sterility is controlled by a single recessive gene. Pollen wall formation in the 83121A mutant was severely defective, with a lack of sporopollenin or exine. To understand the mechanisms of male sterility in 83121A, transcription analysis using RNA-Seq was carried out in the buds of the male sterile line 83121A and the male fertile line 83121B, which are near-isogenic lines differing only in the fertility trait. Via expression analysis of differentially expressed genes involved in pollen exine development before the bicellular pollen stage, BoCYP704B1 was identified as a candidate gene, which was approximately downregulated 30-fold in 83121A. BoCYP704B1 is a member of the evolutionarily conserved CYP704B family, which is essential for sporopollenin formation. The BoCYP704B1 transcript is specifically detected in the developing anthers of wild-type cabbage. Further sequence analysis revealed that a 5424-bp long terminal repeat-retrotransposon (LTR-RT) was inserted into the first exon of BoCYP704B1 in 83121A, which is not found in wild-type plants. The insertion of LTR-RT not only reduced the expression of BoCYP704B1 but also altered structure of protein encoded by BoCYP704B1. Moreover, linkage analysis showed that the homozygotic mutational BoCYP704B1 always cosegregated with male sterility. These data suggest that the LTR-RT insertion in BoCYP704B1 hinders sporopollenin formation in 83121A leading to male sterility. The allele-specific markers developed in this study were effective for marker-assisted selection of the male sterile gene.

Identification and characterisation of Short Interspersed Nuclear Elements in the olive tree (Olea europaea L.) genome

Short Interspersed Nuclear Elements (SINEs) are nonautonomous retrotransposons in the genome of most eukaryotic species. While SINEs have been intensively investigated in humans and other animal systems, SINE identification has been carried out only in a limited number of plant species. This lack of information is apparent especially in non-model plants whose genome has not been sequenced yet. The aim of this work was to produce a specific bioinformatics pipeline for analysing second generation sequence reads of a non-model species and identifying SINEs. We have identified, for the first time, 227 putative SINEs of the olive tree (Olea europaea), that constitute one of the few sets of such sequences in dicotyledonous species. The identified SINEs ranged from 140 to 362 bp in length and were characterised with regard to the occurrence of the tRNA domain in their sequence. The majority of identified elements resulted in single copy or very lowly repeated, often in association with genic sequences. Analysis of sequence similarity allowed us to identify two major groups of SINEs showing different abundances in the olive tree genome, the former with sequence similarity to SINEs of Scrophulariaceae and Solanaceae and the latter to SINEs of Salicaceae. A comparison of sequence conservation between olive SINEs and LTR retrotransposon families suggested that SINE expansion in the genome occurred especially in very ancient times, before LTR retrotransposon expansion, and presumably before the separation of the rosids (to which Oleaceae belong) from the Asterids. Besides providing data on olive SINEs, our results demonstrate the suitability of the pipeline employed for SINE identification. Applying this pipeline will favour further structural and functional analyses on these relatively unknown elements to be performed also in other plant species, even in the absence of a reference genome, and will allow establishing general evolutionary patterns for this kind of repeats in plants.

Construction of BAC Libraries from Flow-Sorted Chromosomes

Cloned DNA libraries in bacterial artificial chromosome (BAC) are the most widely used form of large-insert DNA libraries. BAC libraries are typically represented by ordered clones derived from genomic DNA of a particular organism. In the case of large eukaryotic genomes, whole-genome libraries consist of a hundred thousand to a million clones, which make their handling and screening a daunting task. The labor and cost of working with whole-genome libraries can be greatly reduced by constructing a library derived from a smaller part of the genome. Here we describe construction of BAC libraries from mitotic chromosomes purified by flow cytometric sorting. Chromosome-specific BAC libraries facilitate positional gene cloning, physical mapping, and sequencing in complex plant genomes.

Repetitive DNA and Plant Domestication: Variation in Copy Number and Proximity to Genes of LTR-Retrotransposons among Wild and Cultivated Sunflower (Helianthus annuus) Genotypes

The sunflower (Helianthus annuus) genome contains a very large proportion of transposable elements, especially long terminal repeat retrotransposons. However, knowledge on the retrotransposon-related variability within this species is still limited. We used next-generation sequencing (NGS) technologies to perform a quantitative and qualitative survey of intraspecific variation of the retrotransposon fraction of the genome across 15 genotypes—7 wild accessions and 8 cultivars—of H. annuus. By mapping the Illumina reads of the 15 genotypes onto a library of sunflower long terminal repeat retrotransposons, we observed considerable variability in redundancy among genotypes, at both superfamily and family levels. In another analysis, we mapped Illumina paired reads to two sets of sequences, that is, long terminal repeat retrotransposons and protein-encoding sequences, and evaluated the extent of retrotransposon proximity to genes in the sunflower genome by counting the number of paired reads in which one read mapped to a retrotransposon and the other to a gene. Large variability among genotypes was also ascertained for retrotransposon proximity to genes. Both long terminal repeat retrotransposon redundancy and proximity to genes varied among retrotransposon families and also between cultivated and wild genotypes. Such differences are discussed in relation to the possible role of long terminal repeat retrotransposons in the domestication of sunflower.

Analysis of the repetitive component and retrotransposon population in the genome of a marine angiosperm, Posidonia oceanica (L.) Delile

Posidonia oceanica is a monocotyledonous marine plant that plays a crucial role in maintaining the Mediterranean environment. Despite its ecological importance, basic knowledge of the functional and structural genomics of this species is still limited, as it is for the other seagrasses. Here, for the first time, we report data on the repetitive component of the genome of this seagrass using a low coverage of Illumina sequences and different assembly approaches. A dataset of 19,760 assembled sequences, mostly belonging to the repetitive fraction of the genome, was produced and annotated. Based on mapping Illumina reads onto this dataset, the genome structure of P. oceanica and its repetitive component was inferred. A very large proportion of the genome is represented by long-terminal-repeat (LTR) retrotransposons of both the Copia and Gypsy superfamilies. Posidonia LTR-retrotransposons were classified and their sequences analysed. Gypsy elements belong to three main lineages, while Copia ones belong to seven lineages. Gypsy elements were more represented than Copia ones in the set of assembled sequences and in the genome. Analysis of sequence variability indicated that Gypsy lineages have experienced amplification in more recent times compared to Copia ones.

Genome-wide analysis of LTR-retrotransposons in oil palm

BACKGROUND

The oil palm (Elaeis guineensis Jacq.) is a major cultivated crop and the world's largest source of edible vegetable oil. The genus Elaeis comprises two species E. guineensis, the commercial African oil palm and E. oleifera, which is used in oil palm genetic breeding. The recent publication of both the African oil palm genome assembly and the first draft sequence of its Latin American relative now allows us to tackle the challenge of understanding the genome composition, structure and evolution of these palm genomes through the annotation of their repeated sequences.

METHODS

In this study, we identified, annotated and compared Transposable Elements (TE) from the African and Latin American oil palms. In a first step, Transposable Element databases were built through de novo detection in both genome sequences then the TE content of both genomes was estimated. Then putative full-length retrotransposons with Long Terminal Repeats (LTRs) were further identified in the E. guineensis genome for characterization of their structural diversity, copy number and chromosomal distribution. Finally, their relative expression in several tissues was determined through in silico analysis of publicly available transcriptome data.

RESULTS

Our results reveal a congruence in the transpositional history of LTR retrotransposons between E. oleifera and E. guineensis, especially the Sto-4 family. Also, we have identified and described 583 full-length LTR-retrotransposons in the Elaeis guineensis genome. Our work shows that these elements are most likely no longer mobile and that no recent insertion event has occurred. Moreover, the analysis of chromosomal distribution suggests a preferential insertion of Copia elements in gene-rich regions, whereas Gypsy elements appear to be evenly distributed throughout the genome.

CONCLUSIONS

Considering the high proportion of LTR retrotransposon in the oil palm genome, our work will contribute to a greater understanding of their impact on genome organization and evolution. Moreover, the knowledge gained from this study constitutes a valuable resource for both the improvement of genome annotation and the investigation of the evolutionary history of palms.