Transcriptomic Analyses Reveal Insights into the Shared Regulatory Network of Phenolic Compounds and Steviol Glycosides in Stevia rebaudiana.

Stevia rebaudiana (Bertoni) is a highly valuable crop for the steviol glycoside content in its leaves, which are no-calorie sweeteners hundreds of times more potent than sucrose. The presence of health-promoting phenolic compounds, particularly flavonoids, in the leaf of S. rebaudiana adds further nutritional value to this crop. Although all these secondary metabolites are highly desirable in S. rebaudiana leaves, the genes regulating the biosynthesis of phenolic compounds and the shared gene network between the regulation of biosynthesis of steviol glycosides and phenolic compounds still need to be investigated in this species. To identify putative candidate genes involved in the synergistic regulation of steviol glycosides and phenolic compounds, four genotypes with different contents of these compounds were selected for a pairwise comparison RNA-seq analysis, yielding 1136 differentially expressed genes. Genes that highly correlate with both steviol glycosides and phenolic compound accumulation in the four genotypes of S. rebaudiana were identified using the weighted gene co-expression network analysis. The presence of UDP-glycosyltransferases 76G1, 76H1, 85C1, and 91A1, and several genes associated with the phenylpropanoid pathway, including peroxidase, caffeoyl-CoA O-methyltransferase, and malonyl-coenzyme A:anthocyanin 3-O-glucoside-6″-O-malonyltransferase, along with 21 transcription factors like SCL3, WRK11, and MYB111, implied an extensive and synergistic regulatory network involved in enhancing the production of such compounds in S. rebaudiana leaves. In conclusion, this work identified a variety of putative candidate genes involved in the biosynthesis and regulation of particular steviol glycosides and phenolic compounds that will be useful in gene editing strategies for increasing and steering the production of such compounds in S. rebaudiana as well as in other species.

Future-Proofing Agriculture: De Novo Domestication for Sustainable and Resilient Crops.

The worldwide agricultural system confronts a significant challenge represented by the increasing demand for food in the face of a growing global population. This challenge is exacerbated by a reduction in cultivable land and the adverse effects of climate change on crop yield quantity and quality. Breeders actively embrace cutting-edge omics technologies to pursue resilient genotypes in response to these pressing issues. In this global context, new breeding techniques (NBTs) are emerging as the future of agriculture, offering a solution to introduce resilient crops that can ensure food security, particularly against challenging climate events. Indeed, the search for domestication genes as well as the genetic modification of these loci in wild species using genome editing tools are crucial steps in carrying out de novo domestication of wild plants without compromising their genetic background. Current knowledge allows us to take different paths from those taken by early Neolithic farmers, where crop domestication has opposed natural selection. In this process traits and alleles negatively correlated with high resource environment performance are probably eradicated through artificial selection, while others may have been lost randomly due to domestication and genetic bottlenecks. Thus, domestication led to highly productive plants with little genetic diversity, owing to the loss of valuable alleles that had evolved to tolerate biotic and abiotic stresses. Recent technological advances have increased the feasibility of de novo domestication of wild plants as a promising approach for crafting optimal crops while ensuring food security and using a more sustainable, low-input agriculture. Here, we explore what crucial domestication genes are, coupled with the advancement of technologies enabling the precise manipulation of target sequences, pointing out de novo domestication as a promising application for future crop development.

Discovering the Repeatome of Five Species Belonging to the Asteraceae Family: A Computational Study.

Genome divergence by repeat proliferation and/or loss is a process that plays a crucial role in species evolution. Nevertheless, knowledge of the variability related to repeat proliferation among species of the same family is still limited. Considering the importance of the Asteraceae family, here we present a first contribution towards the metarepeatome of five Asteraceae species. A comprehensive picture of the repetitive components of all genomes was obtained by genome skimming with Illumina sequence reads and by analyzing a pool of full-length long terminal repeat retrotransposons (LTR-REs). Genome skimming allowed us to estimate the abundance and variability of repetitive components. The structure of the metagenome of the selected species was composed of 67% repetitive sequences, of which LTR-REs represented the bulk of annotated clusters. The species essentially shared ribosomal DNA sequences, whereas the other classes of repetitive DNA were highly variable among species. The pool of full-length LTR-REs was retrieved from all the species and their age of insertion was established, showing several lineage-specific proliferation peaks over the last 15-million years. Overall, a large variability of repeat abundance at superfamily, lineage, and sublineage levels was observed, indicating that repeats within individual genomes followed different evolutionary and temporal dynamics, and that different events of amplification or loss of these sequences may have occurred after species differentiation.

Genome-Wide Analysis of WOX Multigene Family in Sunflower (Helianthus annuus L.).

The WUSCHEL-related homeobox (WOX) is a family of specific transcription factors involved in plant development and response to stress, characterized by the presence of a homeodomain. This study represents the first comprehensive characterization of the WOX family in a member of the Asteraceae family, the sunflower (H. annuus L.). Overall, we identified 18 putative HaWOX genes divided by phylogenetic analysis in three major clades (i.e., ancient, intermediate, and WUS). These genes showed conserved structural and functional motifs. Moreover, HaWOX has homogeneously distributed on H. annuus chromosomes. In particular, 10 genes originated after whole segment duplication events, underpinning a possible evolution of this family along with the sunflower genome. In addition, gene expression analysis evidenced a specific pattern of regulation of the putative 18 HaWOX during embryo growth and in ovule and inflorescence meristem differentiation, suggesting a pivotal role for this multigenic family in sunflower development. The results obtained in this work improved the understanding of the WOX multigenic family, providing a resource for future study on functional analysis in an economically valuable species such as sunflower.

Genome-wide identification and characterization of exapted transposable elements in the large genome of sunflower (Helianthus annuus L.).

Transposable elements (TEs) are an important source of genome variability, playing many roles in the evolution of eukaryotic species. Besides well-known phenomena, TEs may undergo the exaptation process and generate the so-called exapted transposable element genes (ETEs). Here we present a genome-wide survey of ETEs in the large genome of sunflower (Helianthus annuus L.), in which the massive amount of TEs, provides a significant source for exaptation. A library of sunflower TEs was used to build TE-specific Hidden Markov Model profiles, to search for all available sunflower gene products. In doing so, 20 016 putative ETEs were identified and further investigated for the characteristics that distinguish TEs from genes, leading to the validation of 3530 ETEs. The analysis of ETEs transcription patterns under different stress conditions showed a differential regulation triggered by treatments mimicking biotic and abiotic stress; furthermore, the distribution of functional domains of differentially regulated ETEs revealed a relevant presence of domains involved in many aspects of cellular functions. A comparative genomic investigation was performed including species representative of Asterids and appropriate outgroups: the bulk of ETEs that resulted were specific to the sunflower, while few ETEs presented orthologues in the genome of all analyzed species, making the hypothesis of a conserved function. This study highlights the crucial role played by exaptation, actively contributing to species evolution.

A chromosome-length genome assembly and annotation of blackberry (Rubus argutus, cv. "Hillquist").

Blackberries (Rubus spp.) are the fourth most economically important berry crop worldwide. Genome assemblies and annotations have been developed for Rubus species in subgenus Idaeobatus, including black raspberry (R. occidentalis), red raspberry (R. idaeus), and R. chingii, but very few genomic resources exist for blackberries and their relatives in subgenus Rubus. Here we present a chromosome-length assembly and annotation of the diploid blackberry germplasm accession "Hillquist" (R. argutus). "Hillquist" is the only known source of primocane-fruiting (annual-fruiting) in tetraploid fresh-market blackberry breeding programs and is represented in the pedigree of many important cultivars worldwide. The "Hillquist" assembly, generated using Pacific Biosciences long reads scaffolded with high-throughput chromosome conformation capture sequencing, consisted of 298 Mb, of which 270 Mb (90%) was placed on 7 chromosome-length scaffolds with an average length of 38.6 Mb. Approximately 52.8% of the genome was composed of repetitive elements. The genome sequence was highly collinear with a novel maternal haplotype-resolved linkage map of the tetraploid blackberry selection A-2551TN and genome assemblies of R. chingii and red raspberry. A total of 38,503 protein-coding genes were predicted, of which 72% were functionally annotated. Eighteen flowering gene homologs within a previously mapped locus aligning to an 11.2 Mb region on chromosome Ra02 were identified as potential candidate genes for primocane-fruiting. The utility of the "Hillquist" genome has been demonstrated here by the development of the first genotyping-by-sequencing-based linkage map of tetraploid blackberry and the identification of possible candidate genes for primocane-fruiting. This chromosome-length assembly will facilitate future studies in Rubus biology, genetics, and genomics and strengthen applied breeding programs.

Characterisation of LTR-Retrotransposons of Stevia rebaudiana and Their Use for the Analysis of Genetic Variability.

Stevia rebaudiana is one of the most important crops belonging to the Asteraceae family. Stevia is cultivated all over the world as it represents a valid natural alternative to artificial sweeteners thanks to its leaves, which produce steviol glycosides that have high sweetening power and reduced caloric value. In this work, the stevia genome sequence was used to isolate and characterise full-length long-terminal repeat retrotransposons (LTR-REs), which account for more than half of the genome. The Gypsy retrotransposons were twice as abundant as the Copia ones. A disproportionate abundance of elements belonging to the Chromovirus/Tekay lineage was observed among the Gypsy elements. Only the SIRE and Angela lineages represented significant portions of the genome among the Copia elements. The dynamics with which LTR-REs colonised the stevia genome were also estimated; all isolated full-length elements turned out to be relatively young, with a proliferation peak around 1-2 million years ago. However, a different analysis conducted by comparing sequences encoding retrotranscriptase showed the occurrence of an older period in which there was a lot of LTR-RE proliferation. Finally, a group of isolated full-length elements belonging to the lineage Angela was used to analyse the genetic variability in 25 accessions of S. rebaudiana using the Inter-Retrotransposon Amplified Polymorphism (IRAP) protocol. The obtained fingerprints highlighted a high degree of genetic variability and were used to study the genomic structures of the different accessions. It was hypothesised that there are four ancestral subpopulations at the root of the analysed accessions, which all turned out to be admixed. Overall, these data may be useful for genome sequence annotations and for evaluating genetic variability in this species, which may be useful in stevia breeding.

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.

In Silico Genome-Wide Characterisation of the Lipid Transfer Protein Multigenic Family in Sunflower (H. annuus L.).

The sunflower (Helianthus annuus L.) is among the most widely cultivated crops in the world due to the oilseed production. Lipid transfer proteins (LTPs) are low molecular mass proteins encoded by a broad multigenic family in higher plants, showing a vast range of functions; these proteins have not been characterised in sunflower at the genomic level. In this work, we exploited the reliable genome sequence of sunflower to identify and characterise the LTP multigenic family in H. annuus. Overall, 101 sunflower putative LTP genes were identified using a homology search and the HMM algorithm. The selected sequences were characterised through phylogenetic analysis, exon-intron organisation, and protein structural motifs. Sunflower LTPs were subdivided into four clades, reflecting their genomic and structural organisation. This gene family was further investigated by analysing the possible duplication origin of genes, which showed the prevalence of tandem and whole genome duplication events, a result that is in line with polyploidisation events that occurred during sunflower genome evolution. Furthermore, LTP gene expression was evaluated on cDNA libraries constructed on six sunflower tissues (leaf, root, ligule, seed, stamen, and pistil) and from roots treated with stimuli mimicking biotic and abiotic stress. Genes encoding LTPs belonging to three out of four clades responded specifically to external stimuli, especially to abscisic acid, auxin, and the saline environment. Interestingly, genes encoding proteins belonging to one clade were expressed exclusively in sunflower seeds. This work is a first attempt of genome-wide identification and characterisation of the LTP multigenic family in a plant species.

Impact of transposable elements on the evolution of complex living systems and their epigenetic control.

Transposable elements (TEs) contribute to genomic innovations, as well as genome instability, across a wide variety of species. Popular designations such as 'selfish DNA' and 'junk DNA,' common in the 1980s, may be either inaccurate or misleading, while a more enlightened view of the TE-host relationship covers a range from parasitism to mutualism. Both plant and animal hosts have evolved epigenetic mechanisms to reduce the impact of TEs, both by directly silencing them and by reducing their ability to transpose in the genome. However, TEs have also been co-opted by both plant and animal genomes to perform a variety of physiological functions, ranging from TE-derived proteins acting directly in normal biological functions to innovations in transcription factor activity and also influencing gene expression. Their presence, in fact, can affect a range of features at genome, phenotype, and population levels. The impact TEs have had on evolution is multifaceted, and many aspects still remain unexplored. In this review, the epigenetic control of TEs is contextualized according to the evolution of complex living systems.

Moderate Salinity Stress Affects Expression of Main Sugar Metabolism and Transport Genes and Soluble Carbohydrate Content in Ripe Fig Fruits (Ficus carica L. cv. Dottato).

Fig trees (Ficus carica L.) are commonly grown in the Mediterranean area, where salinity is an increasing problem in coastal areas. Young, fruiting plants of cv. Dottato were subjected to moderate salt stress (100 mM NaCl added to irrigation water) for 48 days before fruit sampling. To clarify the effect of salinity stress, we investigated changes in the transcription of the main sugar metabolism-related genes involved in the synthesis, accumulation and transport of soluble carbohydrates in ripe fruits by quantitative real-time PCR as well as the content of soluble sugars by quantitative 1H nuclear magnetic resonance spectroscopy. A general increase in the transcript levels of genes involved in the transport of soluble carbohydrates was observed. Alkaline-neutral and Acid Invertases transcripts, related to the synthesis of glucose and fructose, were up-regulated in ripe fruits of NaCl-stressed plants without a change in the content of D-glucose and D-fructose. The increases in sucrose and D-sorbitol contents were likely the result of the up-regulation of the transcription of Sucrose-Synthase- and Sorbitol-Dehydrogenase-encoding genes.

LTR-retrotransposon dynamics in common fig (Ficus carica L.) genome.

BACKGROUND: Long Terminal Repeat retrotransposons (LTR-REs) are repetitive DNA sequences that constitute a large part of the genome. The improvement of sequencing technologies and sequence assembling strategies has achieved genome sequences with much greater reliability than those of the past, especially in relation to repetitive DNA sequences. RESULTS: In this study, we analysed the genome of Ficus carica L., obtained using third generation sequencing technologies and recently released, to characterise the complete complement of full-length LTR-REs to study their dynamics during fig genome evolution. A total of 1867 full-length elements were identified. Those belonging to the Gypsy superfamily were the most abundant; among these, the Chromovirus/Tekay lineage was the most represented. For the Copia superfamily, Ale was the most abundant lineage. Measuring the estimated insertion time of each element showed that, on average, Ivana and Chromovirus/Tekay were the youngest lineages of Copia and Gypsy superfamilies, respectively. Most elements were inactive in transcription, both constitutively and in leaves of plants exposed to an abiotic stress, except for some elements, mostly belonging to the Copia/Ale lineage. A relationship between the inactivity of an element and inactivity of genes lying in close proximity to it was established. CONCLUSIONS: The data reported in this study provide one of the first sets of information on the genomic dynamics related to LTR-REs in a plant species with highly reliable genome sequence. Fig LTR-REs are highly heterogeneous in abundance and estimated insertion time, and only a few elements are transcriptionally active. In general, the data suggested a direct relationship between estimated insertion time and abundance of an element and an inverse relationship between insertion time (or abundance) and transcription, at least for Copia LTR-REs.

DNA Modification Patterns within the Transposable Elements of the Fig (Ficus carica L.) Genome.

Transposable element activity can be harmful to the host's genome integrity, but it can also provide selective advantages. One strategy to cope with transposons is epigenetic control through DNA base modifications. We report the non-canonic DNA modification dynamics of fig (Ficus carica L.) by exploiting high-quality genome reference and related N4-methylcytosine (4mC) and N6-methyladenine (6mA) data. Overall, 1.49% of transposon nucleotides showed either 4mC or 6mA modifications: the 4mC/6mA ratio was similar in Class I and Class II transposons, with a prevalence of 4mC, which is comparable to coding genes. Different percentages of 4mC or 6mA were observed among LTR-retrotransposon lineages and sub-lineages. Furthermore, both the Copia and Gypsy retroelements showed higher modification rates in the LTR and coding regions compared with their neighbour regions. Finally, the unconventional methylation of retrotransposons is unrelated to the number of close genes, suggesting that the 4mC and 6mA frequency in LTR-retrotransposons should not be related to transcriptional repression in the adjacency of the element. In conclusion, this study highlighted unconventional DNA modification patterns in fig transposable elements. Further investigations will focus on functional implications, in regards to how modified retroelements affect the expression of neighbouring genes, and whether these epigenetic markers can spread from repeats to genes, shaping the plant phenotype.

Structural characterization and duplication modes of pseudogenes in plants.

We identified and characterized the pseudogene complements of five plant species: four dicots (Arabidopsis thaliana, Vitis vinifera, Populus trichocarpa and Phaseolus vulgaris) and one monocot (Oryza sativa). Retroposition was considered of modest importance for pseudogene formation in all investigated species except V. vinifera, which showed an unusually high number of retro-pseudogenes in non coding genic regions. By using a pipeline for the classification of sequence duplicates in plant genomes, we compared the relative importance of whole genome, tandem, proximal, transposed and dispersed duplication modes in the pseudo and functional gene complements. Pseudogenes showed higher tendencies than functional genes to genomic dispersion. Dispersed pseudogenes were prevalently fragmented and showed high sequence divergence at flanking regions. On the contrary, those deriving from whole genome duplication were proportionally less than expected based on observations on functional loci and showed higher levels of flanking sequence conservation than dispersed pseudogenes. Pseudogenes deriving from tandem and proximal duplications were in excess compared to functional loci, probably reflecting the high evolutionary rate associated with these duplication modes in plant genomes. These data are compatible with high rates of sequence turnover at neutral sites and double strand break repairs mediated duplication mechanisms.

The plastic genome: The impact of transposable elements on gene functionality and genomic structural variations.

Transposable elements (TEs) are DNA sequences that can change their position within genomes. TEs are present in most organisms and can be an important genomic component. Their activities are manifold: restructuring of genome size, chromosomal rearrangements, induction of gene mutations, and alteration of gene activity by insertion near or within promoters, intronic regions, or enhancer. There are several examples of mutations and other genetic variations determined by the activity of TEs, associated with the evolution of prokaryotic and eukaryotic organisms and the domestication of plants. Generally, TE mobilization occurs when the organism is subjected to stress, which can include both biotic and abiotic stresses, polyploidy conditions, and interspecific hybridizations, very common events in plants. TEs are widely distributed among organisms. TEs also play essential roles in evolution, but most of them are either dormant or inactive. This is mainly determined by epigenetic silencing mechanisms, regulatory systems, and control systems that aim to limit its proliferation. Furthermore, the host has recruited many genes originated from TEs as transcriptional regulators, especially in defense against pathogens and invasive genetic elements; this phenomenon is called molecular domestication. Therefore, TEs are responsible for horizontal gene transfer and the movement of genetic material between organisms, even phylogenetically distant, with a consequent remixing of their gene pools.

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.

On the Trail of Tetu1: Genome-Wide Discovery of CACTA Transposable Elements in Sunflower Genome.

Much has been said about sunflower (Helianthus annuus L.) retrotransposons, representing the majority of the sunflower's repetitive component. By contrast, class II transposons remained poorly described within this species, as they present low sequence conservation and are mostly lacking coding domains, making the identification and characterization of these transposable elements difficult. The transposable element Tetu1, is a non-autonomous CACTA-like element that has been detected in the coding region of a CYCLOIDEA (CYC) gene of a sunflower mutant, tubular ray flower (turf). Based on our knowledge of Tetu1, the publicly available genome of sunflower was fully scanned. A combination of bioinformatics analyses led to the discovery of 707 putative CACTA sequences: 84 elements with complete ends and 623 truncated elements. A detailed characterization of the identified elements allowed further classification into three subgroups of 347 elements on the base of their terminal repeat sequences. Only 39 encode a protein similar to known transposases (TPase), with 10 TPase sequences showing signals of activation. Finally, an analysis of the proximity of CACTA transposons to sunflower genes showed that the majority of CACTA elements are close to the nearest gene, whereas a relevant fraction resides within gene-encoding sequences, likely interfering with sunflower genome functionality and organization.

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.

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.