Genomics of homoploid hybrid speciation: diversity and transcriptional activity of long terminal repeat retrotransposons in hybrid sunflowers

The genomic landscape of transposable elements in yeast hybrids is shaped by structural variation and genotype-specific modulation of transposition rate

Transposable elements (TEs) are major contributors to structural genomic variation by creating interspersed duplications of themselves. In return, structural variants (SVs) can affect the genomic distribution of TE copies and shape their load. One long-standing hypothesis states that hybridization could trigger TE mobilization and thus increase TE load in hybrids. We previously tested this hypothesis (Hénault et al., 2020) by performing a large-scale evolution experiment by mutation accumulation (MA) on multiple hybrid genotypes within and between wild populations of the yeasts Saccharomyces paradoxus and Saccharomyces cerevisiae. Using aggregate measures of TE load with short-read sequencing, we found no evidence for TE load increase in hybrid MA lines. Here, we resolve the genomes of the hybrid MA lines with long-read phasing and assembly to precisely characterize the role of SVs in shaping the TE landscape. Highly contiguous phased assemblies of 127 MA lines revealed that SV types like polyploidy, aneuploidy, and loss of heterozygosity have large impacts on the TE load. We characterized 18 de novo TE insertions, indicating that transposition only has a minor role in shaping the TE landscape in MA lines. Because the scarcity of TE mobilization in MA lines provided insufficient resolution to confidently dissect transposition rate variation in hybrids, we adapted an in vivo assay to measure transposition rates in various S. paradoxus hybrid backgrounds. We found that transposition rates are not increased by hybridization, but are modulated by many genotype-specific factors including initial TE load, TE sequence variants, and mitochondrial DNA inheritance. Our results show the multiple scales at which TE load is shaped in hybrid genomes, being highly impacted by SV dynamics and finely modulated by genotype-specific variation in transposition rates.

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.

Mutation load in sunflower inversions is negatively correlated with inversion heterozygosity

Recombination is critical both for accelerating adaptation and purging deleterious mutations. Chromosomal inversions can act as recombination modifiers that suppress local recombination in heterozygotes and thus, under some conditions, are predicted to accumulate such mutations. In this study, we investigated patterns of recombination, transposable element abundance and coding sequence evolution across the genomes of 1,445 individuals from three sunflower species, as well as within nine inversions segregating within species. We also analyzed the effects of inversion genotypes on 87 phenotypic traits to test for overdominance. We found significant negative correlations of long terminal repeat retrotransposon abundance and deleterious mutations with recombination rates across the genome in all three species. However, we failed to detect an increase in these features in the inversions, except for a modest increase in the proportion of stop codon mutations in several very large or rare inversions. Consistent with this finding, there was little evidence of overdominance of inversions in phenotypes that may relate to fitness. On the other hand, significantly greater load was observed for inversions in populations polymorphic for a given inversion compared to populations monomorphic for one of the arrangements, suggesting that the local state of inversion polymorphism affects deleterious load. These seemingly contradictory results can be explained by the low frequency of inversion heterozygotes in wild sunflower populations, apparently due to divergent selection and associated geographic structure. Inversions contributing to local adaptation represent ideal recombination modifiers, acting to facilitate adaptive divergence with gene flow, while largely escaping the accumulation of deleterious mutations.

OUP accepted manuscript

Interspecific hybridization is often seen as a genomic stress that may lead to new gene expression patterns and deregulation of transposable elements (TEs). The understanding of expression changes in hybrids compared with parental species is essential to disentangle their putative role in speciation processes. However, to date we ignore the detailed mechanisms involved in genomic deregulation in hybrids. We studied the ovarian transcriptome and epigenome of the Drosophila buzzatii and Drosophila koepferae species together with their F₁ hybrid females. We found a trend toward underexpression of genes and TE families in hybrids. The epigenome in hybrids was highly similar to the parental epigenomes and showed intermediate histone enrichments between parental species in most cases. Differential gene expression in hybrids was often associated only with changes in H3K4me3 enrichments, whereas differential TE family expression in hybrids may be associated with changes in H3K4me3, H3K9me3, or H3K27me3 enrichments. We identified specific genes and TE families, which their differential expression in comparison with the parental species was explained by their differential chromatin mark combination enrichment. Finally, cis–trans compensatory regulation could also contribute in some way to the hybrid deregulation. This work provides the first study of histone content in Drosophila interspecific hybrids and their effect on gene and TE expression deregulation.

New Insights into Interspecific Hybridization in Lemna L. Sect. Lemna (Lemnaceae Martinov)

Duckweeds have been increasingly studied in recent years, both as model plants and in view of their potential applications as a new crop in a circular bioeconomy perspective. In order to select species and clones with the desired attributes, the correct identification of the species is fundamental. Molecular methods have recently provided a more solid base for taxonomy and yielded a consensus phylogenetic tree, although some points remain to be elucidated. The duckweed genus Lemna L. comprises twelve species, grouped in four sections, which include very similar sister species. The least taxonomically resolved is sect. Lemna, presenting difficulties in species delimitation using morphological and even barcoding molecular markers. Ambiguous species boundaries between Lemna minor L. and Lemna japonica Landolt have been clarified by Tubulin Based Polymorphism (TBP), with the discovery of interspecific hybrids. In the present work, we extended TBP profiling to a larger number of clones in sect. Lemna, previously classified using only morphological features, in order to test that classification, and to investigate the possible existence of other hybrids in this section. The analysis revealed several misidentifications of clones, in particular among the species L. minor, L. japonica and Lemna gibba L., and identified six putative ‘L. gibba’ clones as interspecific hybrids between L. minor and L. gibba.

Reactivation of transposable elements following hybridization in fission yeast

Hybridization is thought to reactivate transposable elements (TEs) that were efficiently suppressed in the genomes of the parental hosts. Here, we provide evidence for this “genomic shock hypothesis” in the fission yeast Schizosaccharomyces pombe. In this species, two divergent lineages (Sp and Sk) have experienced recent, likely human-induced, hybridization. We used long-read sequencing data to assemble genomes of 37 samples derived from 31 S. pombe strains spanning a wide range of ancestral admixture proportions. A comprehensive TE inventory revealed exclusive presence of long terminal repeat (LTR) retrotransposons. Sequence analysis of active full-length elements, as well as solo LTRs, revealed a complex history of homologous recombination. Population genetic analyses of syntenic sequences placed insertion of many solo LTRs before the split of the Sp and Sk lineages. Most full-length elements were inserted more recently, after hybridization. With the exception of a single full-length element with signs of positive selection, both solo LTRs and, in particular, full-length elements carry signatures of purifying selection indicating effective removal by the host. Consistent with reactivation upon hybridization, the number of full-length LTR retrotransposons, varying extensively from zero to 87 among strains, significantly increases with the degree of genomic admixture. This study gives a detailed account of global TE diversity in S. pombe, documents complex recombination histories within TE elements, and provides evidence for the “genomic shock hypothesis.”

Testing the Genomic Shock Hypothesis Using Transposable Element Expression in Yeast Hybrids

Transposable element (TE) insertions are a source of structural variation and can cause genetic instability and gene expression changes. A host can limit the spread of TEs with various repression mechanisms. Many examples of plant and animal interspecific hybrids show disrupted TE repression leading to TE propagation. Recent studies in yeast did not find any increase in transposition rate in hybrids. However, this does not rule out the possibility that the transcriptional or translational activity of TEs increases following hybridization because of a disruption of the host TE control mechanisms. Thus, whether total expression of a TE family is higher in hybrids than in their parental species remains to be examined. We leveraged publically available RNA-seq and ribosomal profiling data on yeast artificial hybrids of the Saccharomyces genus and performed differential expression analysis of their LTR retrotransposons (Ty elements). Our analyses of total mRNA levels show that Ty elements are generally not differentially expressed in hybrids, even when the hybrids are exposed to a low temperature stress condition. Overall, only 2/26 Ty families show significantly higher expression in the S. cerevisiae × S. uvarum hybrids while there are 3/26 showing significantly lower expression in the S. cerevisiae x S. paradoxus hybrids. Our analysis of ribosome profiling data of S. cerevisiae × S. paradoxus hybrids shows similar translation efficiency of Ty in both parents and hybrids, except for Ty1_cer showing higher translation efficiency. Overall, our results do not support the hypothesis that hybridization could act as a systematic trigger of TE expression in yeast and suggest that the impact of hybridization on TE activity is strain and TE specific.

Transposable Element Mobilization in Interspecific Yeast Hybrids

Barbara McClintock first hypothesized that interspecific hybridization could provide a “genomic shock” that leads to the mobilization of transposable elements (TEs). This hypothesis is based on the idea that regulation of TE movement is potentially disrupted in hybrids. However, the handful of studies testing this hypothesis have yielded mixed results. Here, we set out to identify if hybridization can increase transposition rate and facilitate colonization of TEs in Saccharomyces cerevisiae × Saccharomyces uvarum interspecific yeast hybrids. Saccharomyces cerevisiae have a small number of active long terminal repeat retrotransposons (Ty elements), whereas their distant relative S. uvarum have lost the Ty elements active in S. cerevisiae. Although the regulation system of Ty elements is known in S. cerevisiae, it is unclear how Ty elements are regulated in other Saccharomyces species, and what mechanisms contributed to the loss of most classes of Ty elements in S. uvarum. Therefore, we first assessed whether TEs could insert in the S. uvarum sub-genome of a S. cerevisiae × S. uvarum hybrid. We induced transposition to occur in these hybrids and developed a sequencing technique to show that Ty elements insert readily and nonrandomly in the S. uvarum genome. We then used an in vivo reporter construct to directly measure transposition rate in hybrids, demonstrating that hybridization itself does not alter rate of mobilization. However, we surprisingly show that species-specific mitochondrial inheritance can change transposition rate by an order of magnitude. Overall, our results provide evidence that hybridization can potentially facilitate the introduction of TEs across species boundaries and alter transposition via mitochondrial transmission, but that this does not lead to unrestrained proliferation of TEs suggested by the genomic shock theory.

Evolutionary and genomic comparisons of hybrid uninucleate and nonhybrid Rhizoctonia fungi

The basidiomycetous fungal genus, Rhizoctonia, can cause severe damage to many plants and is composed of multinucleate, binucleate, and uninucleate species differing in pathogenicity. Here we generated chromosome-scale genome assemblies of the three nuclear types of Rhizoctonia isolates. The genomic comparisons revealed that the uninucleate JN strain likely arose by somatic hybridization of two binucleate isolates, and maintained a diploid nucleus. Homeolog gene pairs in the JN genome have experienced both decelerated or accelerated evolution. Homeolog expression dominance occurred between JN subgenomes, in which differentially expressed genes show potentially less evolutionary constraint than the genes without. Analysis of mating-type genes suggested that Rhizoctonia maintains the ancestral tetrapolarity of the Basidiomycota. Long terminal repeat-retrotransposons displayed a reciprocal correlation with the chromosomal GC content in the three chromosome-scale genomes. The more aggressive multinucleate XN strain had more genes encoding enzymes for host cell wall decomposition. These findings demonstrate some evolutionary changes of a recently derived hybrid and in multiple nuclear types of Rhizoctonia.

Hybridization History and Repetitive Element Content in the Genome of a Homoploid Hybrid, Yucca gloriosa (Asparagaceae)

Hybridization in plants results in phenotypic and genotypic perturbations that can have dramatic effects on hybrid physiology, ecology, and overall fitness. Hybridization can also perturb epigenetic control of transposable elements, resulting in their proliferation. Understanding the mechanisms that maintain genomic integrity after hybridization is often confounded by changes in ploidy that occur in hybrid plant species. Homoploid hybrid species, which have no change in chromosome number relative to their parents, offer an opportunity to study the genomic consequences of hybridization in the absence of change in ploidy. Yucca gloriosa (Asparagaceae) is a young homoploid hybrid species, resulting from a cross between Yucca aloifolia and Yucca filamentosa. Previous analyses of ∼11 kb of the chloroplast genome and nuclear-encoded microsatellites implicated a single Y. aloifolia genotype as the maternal parent of Y. gloriosa. Using whole genome resequencing, we assembled chloroplast genomes from 41 accessions of all three species to re-assess the hybrid origins of Y. gloriosa. We further used re-sequencing data to annotate transposon abundance in the three species and mRNA-seq to analyze transcription of transposons. The chloroplast phylogeny and haplotype analysis suggest multiple hybridization events contributing to the origin of Y. gloriosa, with both parental species acting as the maternal donor. Transposon abundance at the superfamily level was significantly different between the three species; the hybrid was frequently intermediate to the parental species in TE superfamily abundance or appeared more similar to one or the other parent. In only one case-Copia LTR transposons-did Y. gloriosa have a significantly higher abundance relative to either parent. Expression patterns across the three species showed little increased transcriptional activity of transposons, suggesting that either no transposon release occurred in Y. gloriosa upon hybridization, or that any transposons that were activated via hybridization were rapidly silenced. The identification and quantification of transposon families paired with expression evidence paves the way for additional work seeking to link epigenetics with the important trait variation seen in this homoploid hybrid system.

Genomic and Transcriptomic Survey Provides New Insight into the Organization and Transposition Activity of Highly Expressed LTR Retrotransposons of Sunflower (Helianthus annuus L.)

LTR retrotransposons (RTEs) play a crucial role in plant genome evolution and adaptation. Although RTEs are generally silenced in somatic plant tissues under non-stressed conditions, some expressed RTEs (exRTEs) escape genome defense mechanisms. As our understanding of exRTE organization in plants is rudimentary, we systematically surveyed the genomic and transcriptomic organization and mobilome (transposition) activity of sunflower (Helianthus annuus L.) exRTEs. We identified 44 transcribed RTEs in the sunflower genome and demonstrated their distinct genomic features: more recent insertion time, longer open reading frame (ORF) length, and smaller distance to neighboring genes. We showed that GAG-encoding ORFs are present at significantly higher frequencies in exRTEs, compared with non-expressed RTEs. Most exRTEs exhibit variation in copy number among sunflower cultivars and one exRTE Gagarin produces extrachromosomal circular DNA in seedling, demonstrating recent and ongoing transposition activity. Nanopore direct RNA sequencing of full-length RTE RNA revealed complex patterns of alternative splicing in RTE RNAs, resulting in isoforms that carry ORFs for distinct RTE proteins. Together, our study demonstrates that tens of expressed sunflower RTEs with specific genomic organization shape the hidden layer of the transcriptome, pointing to the evolution of specific strategies that circumvent existing genome defense mechanisms.

The effect of hybridization on transposable element accumulation in an undomesticated fungal species

Transposable elements (TEs) are mobile genetic elements that can profoundly impact the evolution of genomes and species. A long-standing hypothesis suggests that hybridization could deregulate TEs and trigger their accumulation, although it received mixed support from studies mostly in plants and animals. Here, we tested this hypothesis in fungi using incipient species of the undomesticated yeast Saccharomyces paradoxus. Population genomic data revealed no signature of higher transposition in natural hybrids. As we could not rule out the elimination of past transposition increase signatures by natural selection, we performed a laboratory evolution experiment on a panel of artificial hybrids to measure TE accumulation in the near absence of selection. Changes in TE copy numbers were not predicted by the level of evolutionary divergence between the parents of a hybrid genotype. Rather, they were highly dependent on the individual hybrid genotypes, showing that strong genotype-specific deterministic factors govern TE accumulation in yeast hybrids.

Drosophila Interspecific Hybridization Causes A Deregulation of the piRNA Pathway Genes

Almost all eukaryotes have transposable elements (TEs) against which they have developed defense mechanisms. In the Drosophila germline, the main transposable element (TE) regulation pathway is mediated by specific Piwi-interacting small RNAs (piRNAs). Nonetheless, for unknown reasons, TEs sometimes escape cellular control during interspecific hybridization processes. Because the piRNA pathway genes are involved in piRNA biogenesis and TE control, we sequenced and characterized nine key genes from this pathway in Drosophila buzzatii and Drosophila koepferae species and studied their expression pattern in ovaries of both species and their F1 hybrids. We found that gene structure is, in general, maintained between both species and that two genes—armitage and aubergine—are under positive selection. Three genes—krimper, methyltransferase 2, and zucchini—displayed higher expression values in hybrids than both parental species, while others had RNA levels similar to the parental species with the highest expression. This suggests that the overexpression of some piRNA pathway genes can be a primary response to hybrid stress. Therefore, these results reinforce the hypothesis that TE deregulation may be due to the protein incompatibility caused by the rapid evolution of these genes, leading to a TE silencing failure, rather than to an underexpression of piRNA pathway genes.

Hybrid dysgenesis in Drosophila virilis results in clusters of mitotic recombination and loss-of-heterozygosity but leaves meiotic recombination unaltered

BACKGROUND

Transposable elements (TEs) are endogenous mutagens and their harmful effects are especially evident in syndromes of hybrid dysgenesis. In Drosophila virilis, hybrid dysgenesis is a syndrome of incomplete gonadal atrophy that occurs when males with multiple active TE families fertilize females that lack active copies of the same families. This has been demonstrated to cause the transposition of paternally inherited TE families, with gonadal atrophy driven by the death of germline stem cells. Because there are abundant, active TEs in the male inducer genome, that are not present in the female reactive genome, the D. virilis syndrome serves as an excellent model for understanding the effects of hybridization between individuals with asymmetric TE profiles.

RESULTS

Using the D. virilis syndrome of hybrid dysgenesis as a model, we sought to determine how the landscape of germline recombination is affected by parental TE asymmetry. Using a genotyping-by-sequencing approach, we generated a high-resolution genetic map of D. virilis and show that recombination rate and TE density are negatively correlated in this species. We then contrast recombination events in the germline of dysgenic versus non-dysgenic F1 females to show that the landscape of meiotic recombination is hardly perturbed during hybrid dysgenesis. In contrast, hybrid dysgenesis in the female germline increases transmission of chromosomes with mitotic recombination. Using a de novo PacBio assembly of the D. virilis inducer genome we show that clusters of mitotic recombination events in dysgenic females are associated with genomic regions with transposons implicated in hybrid dysgenesis.

CONCLUSIONS

Overall, we conclude that increased mitotic recombination is likely the result of early TE activation in dysgenic progeny, but a stable landscape of meiotic recombination indicates that either transposition is ameliorated in the adult female germline or that regulation of meiotic recombination is robust to ongoing transposition. These results indicate that the effects of parental TE asymmetry on recombination are likely sensitive to the timing of transposition.

Population Genomics in Rhamdia quelen (Heptapteridae, Siluriformes) Reveals Deep Divergence and Adaptation in the Neotropical Region

Rhamdia quelen, a Neotropical fish with hybridization between highly divergent mitochondrial DNA (mtDNA) lineages, represents an interesting evolutionary model. Previous studies suggested that there might be demographic differences between coastal lagoons and riverine environments, as well as divergent populations that could be reproductively isolated. Here, we investigated the genetic diversity pattern of this taxon in the Southern Neotropical Basin system that includes the La Plata Basin, Patos-Merin lagoon basin and the coastal lagoons draining to the SW Atlantic Ocean, through a population genomics approach using 2b-RAD-sequencing-derived single nucleotide polymorphisms (SNPs). The genomic scan identified selection footprints associated with divergence and suggested local adaptation environmental drivers. Two major genomic clusters latitudinally distributed in the Northern and Southern basins were identified, along with consistent signatures of divergent selection between them. Population structure based on the whole set of loci and on the presumptive neutral vs. adaptive loci showed deep genomic divergence between the two major clusters. Annotation of the most consistent SNPs under divergent selection revealed some interesting candidate genes for further functional studies. Moreover, signals of adaptation to a coastal lagoon environment mediated by purifying selection were found. These new insights provide a better understanding of the complex evolutionary history of R. quelen in the southernmost basin of the Neotropical region.

Importance of Pollinator-Mediated Interspecific Pollen Transfer for Angiosperm Evolution

Understanding how pollen moves between species is critical to understanding speciation, diversification, and evolution of flowering plants. For co-flowering species that share pollinators, competition through interspecific pollen transfer (IPT) can profoundly impact floral evolution, decreasing female fitness via heterospecific pollen deposition on stigmas and male fitness via pollen misplacement during visits to heterospecific flowers. The pollination literature demonstrates that such reproductive interference frequently selects for reproductive character displacement in floral traits linked to pollinator attraction, pollen placement, and mating systems and has also revealed that IPT between given pairs of species is typically asymmetric. More recent work is starting to elucidate its importance to the speciation process, clarifying the link between IPT and current and historical patterns of hybridization, the evolution of phenotypic novelty through adaptive introgression, and the rise of reproductive isolation. Our review aims to stimulate further research on IPT as a ubiquitous mechanism that plays a central role in angiosperm diversification.

Jumping genes: Genomic ballast or powerhouse of biological diversification

Studying hybridization has the potential to elucidate challenging questions in evolutionary biology such as the nature of adaptive genetic variation and reproductive isolation. A growing body of work highlights that the merging of divergent genomes goes beyond the reshuffling of standing variation from related species and promotes mutations (Abbott et al., ). However, to what extent such genome instability generates evolutionary significant variation remains largely elusive. In this issue of Molecular Ecology, Dennenmoser et al. () report considerable dynamics of transposable elements (TEs) in a recent invasive fish species of hybrid origin (Cottus; Figure ). It adds to the recent examples from plants to support TE-specific genome variation following hybridization. Insights from early, as well as established, hybrids are largely coherent with increased TE activity, and this fish system thus represents an inspiring opportunity to further address the possible association between genome dynamics and "rapid evolution of hybrid species." This work based on genome (re)sequencing contrasts with prior transcriptomics or PCR-based studies of TEs and illustrates how unprecedented amount of information promises a better understanding of the multiple patterns of variation across eukaryotic genomes; provided that we get the better of methodological advances. As discussed here, unbiased assessment of TE variation from genome surveys indeed remains a challenge precluding firm conclusions to be reached about the evolutionary significance of TEs. Despite methodological and conceptual developments that appear necessary to unambiguously uncover the unexplored iceberg below the known tip, the role of coding genes vs. TEs in promoting adaptation and speciation might be clarified in a not so remote future.

Hybrid Genome Evolution by Transposition: An Update

Contrary to the view that hybrids are lineages devoid of evolutionary value, a number of case studies that have been lately reported show how hybrids are at the origin of many species. Some well-documented cases demonstrate that bursts of transposition often follow hybridization, generating new genetic variability. Studies in hybrid transposition strongly suggest that epigenetic changes and divergence in piRNA pathways drive deregulation in TE landscapes. Here, I have focused on mechanisms acting in Drosophila hybrids between two cactophilic species. The results reported here show that while hybrid instability by transposition is a genome-wide event, deregulation by TE overexpression in hybrid ovaries is not a general rule. When piRNA pools of ovaries are studied, results show that TEs with parental differences higher than 2-fold in their piRNA amounts are not more commonly deregulated in hybrids than TEs with similar levels, partially discrediting the generality of the maternal cytotype hypothesis. Some promising results on the piRNA pathway global failure hypothesis, which states that accumulated divergence of piRNA effector proteins is responsible for hybrid TE deregulation, have also been obtained. Altogether, these results suggest that TE deregulation might be driven by several interacting mechanisms. A natural scenario is proposed in which genome instability by transposition leads to hybrid genome reorganization. Small hybrid populations, subjected to natural selection helped by genetic drift, evolve new adaptations adapted to novel environments. The final step is either introgression or even a new hybrid species.

Intraspecific and intraorganismal copy number dynamics of retrotransposons and tandem repeat in Aegilops speltoides Tausch (Poaceae, Triticeae)

Transposable elements (TE) and tandem repeats (TR) compose the largest fraction of the plant genome. The abundance and repatterning of repetitive DNA underlie intrapopulation polymorphisms and intraspecific diversification; however, the dynamics of repetitive elements in ontogenesis is not fully understood. Here, we addressed the genotype-specific and tissue-specific abundances and dynamics of the Ty1-copia, Ty3-gypsy, and LINE retrotransposons and species-specific Spelt1 tandem repeat in wild diploid goatgrass, Aegilops speltoides Tausch. Copy numbers of TEs and TR were estimated by real-time quantitative PCR in vegetative and generative tissues in original plants from contrasting allopatric populations and artificial intraspecific hybrids. The results showed that between leaves and somatic spike tissues as well as in progressive microsporogenesis of individual genotypes, the copy numbers of three TEs correlatively oscillated between 2- to 4-fold and the TR copy numbers fluctuated by 18- to 440-fold. Inter-individual and intraorganismal TEs and TR copy number dynamics demonstrate large-scale parallelism with extensive chromosomal rearrangements that were detected using fluorescent in situ hybridization in parental and hybrid genotypes. The data obtained indicate that tissue-specific differences in the abundance and pattern of repetitive sequences emerge during cell proliferation and differentiation in ontogenesis and reflect the reorganization of individual genomes in changing environments, especially in small peripheral population(s) under the influence of rapid climatic changes.

Transposable Element Misregulation Is Linked to the Divergence between Parental piRNA Pathways in Drosophila Hybrids

Interspecific hybridization is a genomic stress condition that leads to the activation of transposable elements (TEs) in both animals and plants. In hybrids between Drosophila buzzatii and Drosophila koepferae, mobilization of at least 28 TEs has been described. However, the molecular mechanisms underlying this TE release remain poorly understood. To give insight on the causes of this TE activation, we performed a TE transcriptomic analysis in ovaries (notorious for playing a major role in TE silencing) of parental species and their F1 and backcrossed (BC) hybrids. We find that 15.2% and 10.6% of the expressed TEs are deregulated in F1 and BC1 ovaries, respectively, with a bias toward overexpression in both cases. Although differences between parental piRNA (Piwi-interacting RNA) populations explain only partially these results, we demonstrate that piRNA pathway proteins have divergent sequences and are differentially expressed between parental species. Thus, a functional divergence of the piRNA pathway between parental species, together with some differences between their piRNA pools, might be at the origin of hybrid instabilities and ultimately cause TE misregulation in ovaries. These analyses were complemented with the study of F1 testes, where TEs tend to be less expressed than in D. buzzatii. This can be explained by an increase in piRNA production, which probably acts as a defence mechanism against TE instability in the male germline. Hence, we describe a differential impact of interspecific hybridization in testes and ovaries, which reveals that TE expression and regulation are sex-biased.

Transcriptomic analysis reveals overdominance playing a critical role in nicotine heterosis in Nicotiana tabacum L

BACKGROUND

As a unique biological phenomenon, heterosis has been concerned with the superior performance of the heterosis than either parents. Despite several F1 hybrids, containing supernal nicotine content, had been discovered and applied to heterosis utilization in Nicotiana tabacum L., nevertheless, the potential molecular mechanism revealing nicotine heterosis has not been illustrated clearly.

RESULT

Phenotypically, the F1 hybrids (Vall6 × Basma) show prominent heterosis in nicotine content by 3 years of field experiments. Transcriptome analysis revealed that genes participating in nicotine anabolism (ADC, PMT, MPO, QPT, AO, QS, QPT, A622, BBLs) and nicotine transport (JAT2, MATE1 and 2, NUP1 and 2) showed an upregulated expression in the hybrid, a majority of which demonstrated an overdominant performance. RT-PCR confirmed that nicotine anabolism was induced in the hybrid.

CONCLUSIONS

These findings strongly suggest that nicotine synthesis and transport efficiency improved in hybrid and overdominance at gene-expression level played a critical role in heterosis of nicotine metabolism.

Genomic abundance and transcriptional activity of diverse gypsy and copia long terminal repeat retrotransposons in three wild sunflower species

BACKGROUND

Long terminal repeat (LTR) retrotransposons are highly abundant in plant genomes and require transcriptional activity for their proliferative mode of replication. These sequences exist in plant genomes as diverse sublineages within the main element superfamilies (i.e., gypsy and copia). While transcriptional activity of these elements is increasingly recognized as a regular attribute of plant transcriptomes, it is currently unknown the extent to which different sublineages of these elements are transcriptionally active both within and across species. In the current report, we utilize next generation sequencing methods to examine genomic copy number abundance of diverse LTR retrotransposon sublineages and their corresponding levels of transcriptional activity in three diploid wild sunflower species, Helianthus agrestis, H. carnosus and H. porteri.

RESULTS

The diploid sunflower species under investigation differ in genome size 2.75-fold, with 2C values of 22.93 for H. agrestis, 12.31 for H. carnosus and 8.33 for H. porteri. The same diverse gypsy and copia sublineages of LTR retrotransposons were identified across species, but with gypsy sequences consistently more abundant than copia and with global gypsy sequence abundance positively correlated with nuclear genome size. Transcriptional activity was detected for multiple copia and gypsy sequences, with significantly higher activity levels detected for copia versus gypsy. Interestingly, of 11 elements identified as transcriptionally active, 5 exhibited detectable expression in all three species and 3 exhibited detectable expression in two species.

CONCLUSIONS

Combined analyses of LTR retrotransposon genomic abundance and transcriptional activity across three sunflower species provides novel insights into genome size evolution and transposable element dynamics in this group. Despite considerable variation in nuclear genome size among species, relatively conserved patterns of LTR retrotransposon transcriptional activity were observed, with a highly overlapping set of copia and gypsy sequences observed to be transcriptionally active across species. A higher proportion of copia versus gypsy elements were found to be transcriptionally active and these sequences also were expressed at higher levels.

Major structural genomic alterations can be associated with hybrid speciation in Aegilops markgrafii (Triticeae)

During evolutionary history many grasses from the tribe Triticeae have undergone interspecific hybridization, resulting in allopolyploidy; whereas homoploid hybrid speciation was found only in rye. Homoeologous chromosomes within the Triticeae preserved cross-species macrocolinearity, except for a few species with rearranged genomes. Aegilops markgrafii, a diploid wild relative of wheat (2n = 2x = 14), has a highly asymmetrical karyotype that is indicative of chromosome rearrangements. Molecular cytogenetics and next-generation sequencing were used to explore the genome organization. Fluorescence in situ hybridization with a set of wheat cDNAs allowed the macrostructure and cross-genome homoeology of the Ae. markgrafii chromosomes to be established. Two chromosomes maintained colinearity, whereas the remaining were highly rearranged as a result of inversions and inter- and intrachromosomal translocations. We used sets of barley and wheat orthologous gene sequences to compare discrete parts of the Ae. markgrafii genome involved in the rearrangements. Analysis of sequence identity profiles and phylogenic relationships grouped chromosome blocks into two distinct clusters. Chromosome painting revealed the distribution of transposable elements and differentiated chromosome blocks into two groups consistent with the sequence analyses. These data suggest that introgressive hybridization accompanied by gross chromosome rearrangements might have had an impact on karyotype evolution and homoploid speciation in Ae. markgrafii.

Copy number increases of transposable elements and protein-coding genes in an invasive fish of hybrid origin

Evolutionary dynamics of structural genetic variation in lineages of hybrid origin is not well explored, although structural mutations may increase in controlled hybrid crosses. We therefore tested whether structural variants accumulate in a fish of recent hybrid origin, invasive Cottus, relative to both parental species Cottus rhenanus and Cottus perifretum. Copy-number variation in exons of 10,979 genes was assessed using comparative genome hybridization arrays. Twelve genes showed significantly higher copy numbers in invasive Cottus compared to both parents. This coincided with increased expression for three genes related to vision, detoxification and muscle development, suggesting possible gene dosage effects. Copy number increases of putative transposons were assessed by comparative mapping of genomic DNA reads against a de novo assembly of 1,005 repetitive elements. In contrast to exons, copy number increases of repetitive elements were common (20.7%) in invasive Cottus, whereas decrease was very rare (0.01%). Among the increased repetitive elements, 53.8% occurred at higher numbers in C. perifretum compared to C. rhenanus, while only 1.4% were more abundant in C. rhenanus. This implies a biased mutational process that amplifies genetic material from one ancestor. To assess the frequency of de novo mutations through hybridization, we screened 64 laboratory-bred F2 offspring between the parental species for copy-number changes at five candidate loci. We found no evidence for new structural variants, indicating that they are too rare to be detected given our sampling scheme. Instead, they must have accumulated over more generations than we observed in a controlled cross.

Long Terminal Repeat Retrotransposon Content in Eight Diploid Sunflower Species Inferred from Next-Generation Sequence Data

The most abundant transposable elements (TEs) in plant genomes are Class I long terminal repeat (LTR) retrotransposons represented by superfamilies gypsy and copia. Amplification of these superfamilies directly impacts genome structure and contributes to differential patterns of genome size evolution among plant lineages. Utilizing short-read Illumina data and sequence information from a panel of Helianthus annuus (sunflower) full-length gypsy and copia elements, we explore the contribution of these sequences to genome size variation among eight diploid Helianthus species and an outgroup taxon, Phoebanthus tenuifolius. We also explore transcriptional dynamics of these elements in both leaf and bud tissue via RT-PCR. We demonstrate that most LTR retrotransposon sublineages (i.e., families) display patterns of similar genomic abundance across species. A small number of LTR retrotransposon sublineages exhibit lineage-specific amplification, particularly in the genomes of species with larger estimated nuclear DNA content. RT-PCR assays reveal that some LTR retrotransposon sublineages are transcriptionally active across all species and tissue types, whereas others display species-specific and tissue-specific expression. The species with the largest estimated genome size, H. agrestis, has experienced amplification of LTR retrotransposon sublineages, some of which have proliferated independently in other lineages in the Helianthus phylogeny.

Drosophila Females Undergo Genome Expansion after Interspecific Hybridization

Genome size (or C-value) can present a wide range of values among eukaryotes. This variation has been attributed to differences in the amplification and deletion of different noncoding repetitive sequences, particularly transposable elements (TEs). TEs can be activated under different stress conditions such as interspecific hybridization events, as described for several species of animals and plants. These massive transposition episodes can lead to considerable genome expansions that could ultimately be involved in hybrid speciation processes. Here, we describe the effects of hybridization and introgression on genome size of Drosophila hybrids. We measured the genome size of two close Drosophila species, Drosophila buzzatii and Drosophila koepferae, their F1 offspring and the offspring from three generations of backcrossed hybrids; where mobilization of up to 28 different TEs was previously detected. We show that hybrid females indeed present a genome expansion, especially in the first backcross, which could likely be explained by transposition events. Hybrid males, which exhibit more variable C-values among individuals of the same generation, do not present an increased genome size. Thus, we demonstrate that the impact of hybridization on genome size can be detected through flow cytometry and is sex-dependent.

Evidence for extensive parallelism but divergent genomic architecture of adaptation along altitudinal and latitudinal gradients in Populus trichocarpa

Adaptation to climate across latitude and altitude reflects shared climatic constraints, which may lead to parallel adaptation. However, theory predicts that higher gene flow should favor more concentrated genomic architectures, which would lead to fewer locally maladapted recombinants. We used exome capture to resequence the gene space along a latitudinal and two altitudinal transects in the model tree Populus trichocapra. Adaptive trait phenotyping was coupled with FST outlier tests and sliding window analysis to assess the degree of parallel adaptation as well as the genomic distribution of outlier loci. Up to 51% of outlier loci overlapped between transect pairs and up to 15% of these loci overlapped among all three transects. Genomic clustering of adaptive loci was more pronounced for altitudinal than latitudinal transects. In both altitudinal transects, there was a larger number of these 'islands of divergence', which were on average longer and included several of exceptional physical length. Our results suggest that recapitulation of genetic clines over latitude and altitude involves extensive parallelism, but that steep altitudinal clines generate islands of divergence. This suggests that physical proximity of genes in coadapted complexes may buffer against the movement of maladapted alleles from geographically proximal but climatically distinct populations.

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

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

Does hybridization with a widespread congener threaten the long-term persistence of the Eastern Alpine rare local endemic Knautia carinthiaca?

Interspecific hybridization, especially when regularly followed by backcrossing (i.e., introgressive hybridization), conveys a substantial risk for many endangered organisms. This is particularly true for narrow endemics occurring within distributional ranges of widespread congeners. An excellent example is provided by the plant genus Knautia (Caprifoliaceae): Locally endemic K. carinthiaca is reported from two isolated populations in southern Austria situated within an area predominantly occupied by widespread K. arvensis. While K. carinthiaca usually inhabits low-competition communities on rocky outcrops, K. arvensis occurs mainly in dry to mesic managed grasslands, yet both species can coexist in marginal environments and were suspected to hybridize. Flow cytometry revealed that diploid K. carinthiaca only occurs at its locus classicus, whereas the second locality is inhabited by the morphologically similar but tetraploid K. norica. In the, therefore, single population of K. carinthiaca, flow cytometry and AFLP fingerprinting showed signs of introgressive hybridization with diploid K. arvensis. Hybridization patterns were also reflected in intermediate habitat preferences and morphology of the hybrids. Environmental barriers to gene flow seem to prevent genetic erosion of K. carinthiaca individuals from the core ecological niches, restricting most introgressed individuals to peripheral habitats. Efficient conservation of K. carinthiaca will require strict protection of its habitat and ban on forest clear cuts in a buffer zone to prevent invasion of K. arvensis. We demonstrate the large potential of multidisciplinary approaches combining molecular, cytometric, and ecological tools for a reliable inventory and threat assessment of rare species.

Process and pattern in cichlid radiations - inferences for understanding unusually high rates of evolutionary diversification

The cichlid fish radiations in the African Great Lakes differ from all other known cases of rapid speciation in vertebrates by their spectacular trophic diversity and richness of sympatric species, comparable to the most rapid angiosperm radiations. I review factors that may have facilitated these radiations and compare these with insights from recent work on plant radiations. Work to date suggests that it was a coincidence of ecological opportunity, intrinsic ecological versatility and genomic flexibility, rapidly evolving behavioral mate choice and large amounts of standing genetic variation that permitted these spectacular fish radiations. I propose that spatially orthogonal gradients in the fit of phenotypes to the environment facilitate speciation because they allow colonization of alternative fitness peaks during clinal speciation despite local disruptive selection. Such gradients are manifold in lakes because of the interaction of water depth as an omnipresent third spatial dimension with other fitness-relevant variables. I introduce a conceptual model of adaptive radiation that integrates these elements and discuss its applicability to, and predictions for, plant radiations.