Transposable elements drive rapid phenotypic variation in Capsella rubella

The molecular basis of phenotypic evolution: beyond the usual suspects

It has been well documented that mutations in coding DNA or cis-regulatory elements underlie natural phenotypic variation in many organisms. However, the development of sophisticated functional tools in recent years in a wide range of traditionally non-model systems have revealed many 'unusual suspects' in the molecular bases of phenotypic evolution, including upstream open reading frames (uORFs), cryptic splice sites, and small RNAs. Furthermore, large-scale genome sequencing, especially long-read sequencing, has identified a cornucopia of structural variation underlying phenotypic divergence and elucidated the composition of supergenes that control complex multi-trait polymorphisms. In this review article we highlight recent studies that demonstrate this great diversity of molecular mechanisms producing adaptive genetic variation and the panoply of evolutionary paths leading to the 'grandeur of life'.

Transposable elements cause the loss of self-incompatibility in citrus

Self-incompatibility (SI) is a widespread prezygotic mechanism for flowering plants to avoid inbreeding depression and promote genetic diversity. Citrus has an S-RNase-based SI system, which was frequently lost during evolution. We previously identified a single nucleotide mutation in Sm-RNase, which is responsible for the loss of SI in mandarin and its hybrids. However, little is known about other mechanisms responsible for conversion of SI to self-compatibility (SC) and we identify a completely different mechanism widely utilized by citrus. Here, we found a 786-bp miniature inverted-repeat transposable element (MITE) insertion in the promoter region of the FhiS2-RNase in Fortunella hindsii Swingle (a model plant for citrus gene function), which does not contain the Sm-RNase allele but are still SC. We demonstrate that this MITE plays a pivotal role in the loss of SI in citrus, providing evidence that this MITE insertion prevents expression of the S-RNase; moreover, transgenic experiments show that deletion of this 786-bp MITE insertion recovers the expression of FhiS2-RNase and restores SI. This study identifies the first evidence for a role for MITEs at the S-locus affecting the SI phenotype. A family-wide survey of the S-locus revealed that MITE insertions occur frequently adjacent to S-RNase alleles in different citrus genera, but only certain MITEs appear to be responsible for the loss of SI. Our study provides evidence that insertion of MITEs into a promoter region can alter a breeding strategy and suggests that this phenomenon may be broadly responsible for SC in species with the S-RNase system.

The evolution of transposable elements in Brachypodium distachyon is governed by purifying selection, while neutral and adaptive processes play a minor role

Understanding how plants adapt to changing environments and the potential contribution of transposable elements (TEs) to this process is a key question in evolutionary genomics. While TEs have recently been put forward as active players in the context of adaptation, few studies have thoroughly investigated their precise role in plant evolution. Here, we used the wild Mediterranean grass Brachypodium distachyon as a model species to identify and quantify the forces acting on TEs during the adaptation of this species to various conditions, across its entire geographic range. Using sequencing data from more than 320 natural B. distachyon accessions and a suite of population genomics approaches, we reveal that putatively adaptive TE polymorphisms are rare in wild B. distachyon populations. After accounting for changes in past TE activity, we show that only a small proportion of TE polymorphisms evolved neutrally (<10%), while the vast majority of them are under moderate purifying selection regardless of their distance to genes. TE polymorphisms should not be ignored when conducting evolutionary studies, as they can be linked to adaptation. However, our study clearly shows that while they have a large potential to cause phenotypic variation in B. distachyon, they are not favored during evolution and adaptation over other types of mutations (such as point mutations) in this species.

Population genetic analysis based on the polymorphisms mediated by transposons in the genomes of pig

Transposable elements (TEs) mobility is capable of generating a large number of structural variants (SVs), which can have considerable potential as molecular markers for genetic analysis and molecular breeding in livestock. Our results showed that the pig genome contains mainly TE-SVs generated by short interspersed nuclear elements (51,873/76.49%), followed by long interspersed nuclear elements (11,131/16.41%), and more than 84% of the common TE-SVs (Minor allele frequency, MAF > 0.10) were validated to be polymorphic. Subsequently, we utilized the identified TE-SVs to gain insights into the population structure, resulting in clear differentiation among the three pig groups and facilitating the identification of relationships within Chinese local pig breeds. In addition, we investigated the frequencies of TEs in the gene coding regions of different pig groups and annotated the respective TE types, related genes, and functional pathways. Through genome-wide comparisons of Large White pigs and Chinese local pigs utilizing the Beijing Black pigs, we identified TE-mediated SVs associated with quantitative trait loci and observed that they were mainly involved in carcass traits and meat quality traits. Lastly, we present the first documented evidence of TE transduction in the pig genome.

Forces driving transposable element load variation during Arabidopsis range expansion

Genetic load refers to the accumulated and potentially life-threatening deleterious mutations in populations. Understanding the mechanisms underlying genetic load variation of transposable element (TE) insertion, a major large-effect mutation, during range expansion is an intriguing question in biology. Here, we used 1,115 global natural accessions of Arabidopsis (Arabidopsis thaliana) to study the driving forces of TE load variation during its range expansion. TE load increased with range expansion, especially in the recently established Yangtze River basin population. Effective population size, which explains 62.0% of the variance in TE load, high transposition rate, and selective sweeps contributed to TE accumulation in the expanded populations. We genetically mapped and identified multiple candidate causal genes and TEs, and revealed the genetic architecture of TE load variation. Overall, this study reveals the variation in TE genetic load during Arabidopsis expansion and highlights the causes of TE load variation from the perspectives of both population genetics and quantitative genetics.

Post-invasion selection acts on standing genetic variation despite a severe founding bottleneck

Invasive populations often have lower genetic diversity relative to the native-range populations from which they derive.1,2 Despite this, many biological invaders succeed in their new environments, in part due to rapid adaptation.3,4,5,6 Therefore, the role of genetic bottlenecks in constraining the adaptation of invaders is debated.7,8,9,10 Here, we use whole-genome resequencing of samples from a 10-year time-series dataset, representing the natural invasion of the Asian honey bee (Apis cerana) in Australia, to investigate natural selection occurring in the aftermath of a founding event. We find that Australia's A. cerana population was founded by as few as one colony, whose arrival was followed by a period of rapid population expansion associated with an increase of rare variants.11 The bottleneck resulted in a steep loss of overall genetic diversity, yet we nevertheless detected loci with signatures of positive selection during the first years post-invasion. When we investigated the origin of alleles under selection, we found that selection acted primarily on the variation introduced by founders and not on the variants that arose post-invasion by mutation. In all, our data highlight that selection on standing genetic variation can occur in the early years post-invasion, even where founding bottlenecks are severe.

Transposable Element Insertions Are Associated with Batesian Mimicry in the Pantropical Butterfly Hypolimnas misippus

Hypolimnas misippus is a Batesian mimic of the toxic African Queen butterfly (Danaus chrysippus). Female H. misippus butterflies use two major wing patterning loci (M and A) to imitate three color morphs of D. chrysippus found in different regions of Africa. In this study, we examine the evolution of the M locus and identify it as an example of adaptive atavism. This phenomenon involves a morphological reversion to an ancestral character that results in an adaptive phenotype. We show that H. misippus has re-evolved an ancestral wing pattern present in other Hypolimnas species, repurposing it for Batesian mimicry of a D. chrysippus morph. Using haplotagging, a linked-read sequencing technology, and our new analytical tool, Wrath, we discover two large transposable element insertions located at the M locus and establish that these insertions are present in the dominant allele responsible for producing mimetic phenotype. By conducting a comparative analysis involving additional Hypolimnas species, we demonstrate that the dominant allele is derived. This suggests that, in the derived allele, the transposable elements disrupt a cis-regulatory element, leading to the reversion to an ancestral phenotype that is then utilized for Batesian mimicry of a distinct model, a different morph of D. chrysippus. Our findings present a compelling instance of convergent evolution and adaptive atavism, in which the same pattern element has independently evolved multiple times in Hypolimnas butterflies, repeatedly playing a role in Batesian mimicry of diverse model species.

Genomic asymmetric epigenetic modification of transposable elements is involved in gene expression regulation of allopolyploid Brassica napus

Polyploids are common and have a wide geographical distribution and environmental adaptability. Allopolyploidy may lead to the activation of transposable elements (TE). However, the mechanism of epigenetic modification of TEs in the establishment and evolution of allopolyploids remains to be explored. We focused on the TEs of model allopolyploid Brassica napus (An An Cn Cn ), exploring the TE characteristics of the genome, epigenetic modifications of TEs during allopolyploidization, and regulation of gene expression by TE methylation. In B. napus, approximately 50% of the genome was composed of TEs. TEs increased with proximity to genes, especially DNA transposons. TE methylation levels were negatively correlated with gene expression, and changes in TE methylation levels were able to regulate the expression of neighboring genes related to responses to light intensity and stress, which promoted powerful adaptation of allopolyploids to new environments. TEs can be synergistically regulated by RNA-directed DNA methylation pathways and histone modifications. The epigenetic modification levels of TEs tended to be similar to those of the diploid parents during the genome evolution of B. napus. The TEs of the An subgenome were more likely to be modified, and the imbalance in TE number and epigenetic modification level in the An and Cn subgenomes may lead to the establishment of subgenome dominance. Our study analyzed the characteristics of TE location, DNA methylation, siRNA, and histone modification in B. napus and highlighted the importance of TE epigenetic modifications during the allopolyploidy process, providing support for revealing the mechanism of allopolyploid formation and evolution.

The chromosome-level genome assembly of Fraxinus americana provides insights into the evolution of Oleaceae plants

White ash (Fraxinus americana linn.) originates from the southeastern United States. It is a tall and fast-growing tree species with strong salt-alkali resistance and cold tolerance, making it an important reforestation species and widely planted worldwide. Here, we completed the chromosome-level reference genome assembly of F. americana based on Illumina, PacBio, and Hi-C reads, with a genome size of 878.98 Mb, an N50 of 3.27 Mb, and a heterozygosity rate of 0.3 %. Based on de novo prediction, transcriptome prediction, and homology-based protein prediction, we obtained 39,538 genes. Approximately 843.21 Mb of the assembly genome was composed of 37,928 annotated protein-coding genes, with a gene function annotation rate of 95.93 %. 99.94 % of the overlap clusters (877.44 Mb) were anchored to 23 chromosomes. Synteny analysis of F. americana and other Oleaceae plants showed that F. americana underwent frequent chromosome rearrangements. The amplification of the Ale transposons effectively promoted the genome size of F. americana. Compared with other Oleaceae plants, the Glutathione S-transferase (GST) gene family in the F. americana genome has undergone significant expansion, which may help F. americana cope with adverse natural environments. Furthermore, we found that key enzyme-coding gene families related to lignin biosynthesis were expanded and highly expressed in F. americana leaves. These key genes drive lignin synthesis and benefit F. americana in fast-growing, as well as resisting biotic and abiotic stress. Overall, the F. americana genome assembly provides insights into the evolution of Oleaceae plants and provides abundant resources for breeding and germplasm conservation of white ash.

Multi-integrated genomic data for Passiflora foetida provides insights into genome size evolution and floral development in Passiflora

Passiflora is a plant genus known for its extremely distinctive and colorful flowers and a wide range of genome size variation. However, how genome characteristics are related to flower traits among Passiflora species remains poorly understood. Here, we assembled a chromosome-scale genome of P. foetida, which belongs to the same subgenus as the commercial passionfruit P. edulis. The genome of P. foetida is smaller (424.16 Mb) and contains fewer copies of long terminal repeat retrotransposons (LTR-RTs). The disparity in LTR-RTs is one of the main contributors to the differences in genome sizes between these two species and possibly in floral traits. Additionally, we observed variation in insertion times and copy numbers of LTR-RTs across different transposable element (TE) lineages. Then, by integrating transcriptomic data from 33 samples (eight floral organs and flower buds at three developmental stages) with phylogenomic and metabolomic data, we conducted an in-depth analysis of the expression, phylogeny, and copy number of MIKC-type MADS-box genes and identified essential biosynthetic genes responsible for flower color and scent from glandular bracts and other floral organs. Our study pinpoints LRT-RTs as an important player in genome size variation in Passiflora species and provides insights into future genetic improvement.

Algorithms for the Uniqueness of the Longest Common Subsequence

Given several number sequences, determining the longest common subsequence is a classical problem in computer science. This problem has applications in bioinformatics, especially determining transposable genes. Nevertheless, related works only consider how to find one longest common subsequence. In this paper, we consider how to determine the uniqueness of the longest common subsequence. If there are multiple longest common subsequences, we also determine which number appears in all/some/none of the longest common subsequences. We focus on four scenarios: (1) linear sequences without duplicated numbers; (2) circular sequences without duplicated numbers; (3) linear sequences with duplicated numbers; (4) circular sequences with duplicated numbers. We develop corresponding algorithms and apply them to gene sequencing data.

The genetic basis of epigenetic variation and its consequences for adaptation

Recent population genomic studies in plants have shed new light on natural epigenetic variation by identifying key genetic determinants, "trans modifiers," that influence epigenetic states genome-wide and their interplay with environmental factors. Here, we review this progress by focusing on the epigenetic control of transposition and life-cycle transitions to highlight the ecological consequences of this genetic architecture and its evolutionary significance. This knowledge provides new opportunities to address long-standing questions about the establishment of environment-associated epigenetic variation and its relevance in adaptation.

The pig pangenome provides insights into the roles of coding structural variations in genetic diversity and adaptation

Structural variations have emerged as an important driving force for genome evolution and phenotypic variation in various organisms, yet their contributions to genetic diversity and adaptation in domesticated animals remain largely unknown. Here we constructed a pangenome based on 250 sequenced individuals from 32 pig breeds in Eurasia and systematically characterized coding sequence presence/absence variations (PAVs) within pigs. We identified 308.3-Mb nonreference sequences and 3438 novel genes absent from the current reference genome. Gene PAV analysis showed that 16.8% of the genes in the pangene catalog undergo PAV. A number of newly identified dispensable genes showed close associations with adaptation. For instance, several novel swine leukocyte antigen (SLA) genes discovered in nonreference sequences potentially participate in immune responses to productive and respiratory syndrome virus (PRRSV) infection. We delineated previously unidentified features of the pig mobilome that contained 490,480 transposable element insertion polymorphisms (TIPs) resulting from recent mobilization of 970 TE families, and investigated their population dynamics along with influences on population differentiation and gene expression. In addition, several candidate adaptive TE insertions were detected to be co-opted into genes responsible for responses to hypoxia, skeletal development, regulation of heart contraction, and neuronal cell development, likely contributing to local adaptation of Tibetan wild boars. These findings enhance our understanding on hidden layers of the genetic diversity in pigs and provide novel insights into the role of SVs in the evolutionary adaptation of mammals.

Bioinformatics Analysis of MSH1 Genes of Green Plants: Multiple Parallel Length Expansions, Intron Gains and Losses, Partial Gene Duplications, and Alternative Splicing

MutS homolog 1 (MSH1) is involved in the recombining and repairing of organelle genomes and is essential for maintaining their stability. Previous studies indicated that the length of the gene varied greatly among species and detected species-specific partial gene duplications in Physcomitrella patens. However, there are critical gaps in the understanding of the gene size expansion, and the extent of the partial gene duplication of MSH1 remains unclear. Here, we screened MSH1 genes in 85 selected species with genome sequences representing the main clades of green plants (Viridiplantae). We identified the MSH1 gene in all lineages of green plants, except for nine incomplete species, for bioinformatics analysis. The gene is a singleton gene in most of the selected species with conserved amino acids and protein domains. Gene length varies greatly among the species, ranging from 3234 bp in Ostreococcus tauri to 805,861 bp in Cycas panzhihuaensis. The expansion of MSH1 repeatedly occurred in multiple clades, especially in Gymnosperms, Orchidaceae, and Chloranthus spicatus. MSH1 has exceptionally long introns in certain species due to the gene length expansion, and the longest intron even reaches 101,025 bp. And the gene length is positively correlated with the proportion of the transposable elements (TEs) in the introns. In addition, gene structure analysis indicated that the MSH1 of green plants had undergone parallel intron gains and losses in all major lineages. However, the intron number of seed plants (gymnosperm and angiosperm) is relatively stable. All the selected gymnosperms contain 22 introns except for Gnetum montanum and Welwitschia mirabilis, while all the selected angiosperm species preserve 21 introns except for the ANA grade. Notably, the coding region of MSH1 in algae presents an exceptionally high GC content (47.7% to 75.5%). Moreover, over one-third of the selected species contain species-specific partial gene duplications of MSH1, except for the conserved mosses-specific partial gene duplication. Additionally, we found conserved alternatively spliced MSH1 transcripts in five species. The study of MSH1 sheds light on the evolution of the long genes of green plants.

Selection signatures and population dynamics of transposable elements in lima bean

The domestication process in lima bean (Phaseolus lunatus L.) involves two independent events, within the Mesoamerican and Andean gene pools. This makes lima bean an excellent model to understand convergent evolution. The mechanisms of adaptation followed by Mesoamerican and Andean landraces are largely unknown. Genes related to these adaptations can be selected by identification of selective sweeps within gene pools. Previous genetic analyses in lima bean have relied on Single Nucleotide Polymorphism (SNP) loci, and have ignored transposable elements (TEs). Here we show the analysis of whole-genome sequencing data from 61 lima bean accessions to characterize a genomic variation database including TEs and SNPs, to associate selective sweeps with variable TEs and to predict candidate domestication genes. A small percentage of genes under selection are shared among gene pools, suggesting that domestication followed different genetic avenues in both gene pools. About 75% of TEs are located close to genes, which shows their potential to affect gene functions. The genetic structure inferred from variable TEs is consistent with that obtained from SNP markers, suggesting that TE dynamics can be related to the demographic history of wild and domesticated lima bean and its adaptive processes, in particular selection processes during domestication.

Long-read direct RNA sequencing reveals epigenetic regulation of chimeric gene-transposon transcripts in Arabidopsis thaliana

Transposable elements (TEs) are accumulated in both intergenic and intragenic regions in plant genomes. Intragenic TEs often act as regulatory elements of associated genes and are also co-transcribed with genes, generating chimeric TE-gene transcripts. Despite the potential impact on mRNA regulation and gene function, the prevalence and transcriptional regulation of TE-gene transcripts are poorly understood. By long-read direct RNA sequencing and a dedicated bioinformatics pipeline, ParasiTE, we investigated the transcription and RNA processing of TE-gene transcripts in Arabidopsis thaliana. We identified a global production of TE-gene transcripts in thousands of A. thaliana gene loci, with TE sequences often being associated with alternative transcription start sites or transcription termination sites. The epigenetic state of intragenic TEs affects RNAPII elongation and usage of alternative poly(A) signals within TE sequences, regulating alternative TE-gene isoform production. Co-transcription and inclusion of TE-derived sequences into gene transcripts impact regulation of RNA stability and environmental responses of some loci. Our study provides insights into TE-gene interactions that contributes to mRNA regulation, transcriptome diversity, and environmental responses in plants.

Molecular mechanisms of adaptive evolution in wild animals and plants

Wild animals and plants have developed a variety of adaptive traits driven by adaptive evolution, an important strategy for species survival and persistence. Uncovering the molecular mechanisms of adaptive evolution is the key to understanding species diversification, phenotypic convergence, and inter-species interaction. As the genome sequences of more and more non-model organisms are becoming available, the focus of studies on molecular mechanisms of adaptive evolution has shifted from the candidate gene method to genetic mapping based on genome-wide scanning. In this study, we reviewed the latest research advances in wild animals and plants, focusing on adaptive traits, convergent evolution, and coevolution. Firstly, we focused on the adaptive evolution of morphological, behavioral, and physiological traits. Secondly, we reviewed the phenotypic convergences of life history traits and responding to environmental pressures, and the underlying molecular convergence mechanisms. Thirdly, we summarized the advances of coevolution, including the four main types: mutualism, parasitism, predation and competition. Overall, these latest advances greatly increase our understanding of the underlying molecular mechanisms for diverse adaptive traits and species interaction, demonstrating that the development of evolutionary biology has been greatly accelerated by multi-omics technologies. Finally, we highlighted the emerging trends and future prospects around the above three aspects of adaptive evolution.

Neutral and adaptive genetic diversity in plants: An overview

Genetic diversity is a prerequisite for evolutionary change in all kinds of organisms. It is generally acknowledged that populations lacking genetic variation are unable to evolve in response to new environmental conditions (e.g., climate change) and thus may face an increased risk of extinction. Although the importance of incorporating genetic diversity into the design of conservation measures is now well understood, less attention has been paid to the distinction between neutral (NGV) and adaptive (AGV) genetic variation. In this review, we first focus on the utility of NGV by examining the ways to quantify it, reviewing applications of NGV to infer ecological and evolutionary processes, and by exploring its utility in designing conservation measures for plant populations and species. Against this background, we then summarize the ways to identify and estimate AGV and discuss its potential use in plant conservation. After comparing NGV and AGV and considering their pros and cons in a conservation context, we conclude that there is an urgent need for a better understanding of AGV and its role in climate change adaptation. To date, however, there are only a few AGV studies on non-model plant species aimed at deciphering the genetic and genomic basis of complex trait variation. Therefore, conservation researchers and practitioners should keep utilizing NGV to develop relevant strategies for rare and endangered plant species until more estimates of AGV are available.

The nature and genomic landscape of repetitive DNA classes in Chrysanthemum nankingense shows recent genomic changes

BACKGROUND AND AIMS

Tandemly repeated DNA and transposable elements represent most of the DNA in higher plant genomes. High-throughput sequencing allows a survey of the DNA in a genome, but whole-genome assembly can miss a substantial fraction of highly repeated sequence motifs. Chrysanthemum nankingense (2n = 2x = 18; genome size = 3.07 Gb; Asteraceae), a diploid reference for the many auto- and allopolyploids in the genus, was considered as an ancestral species and serves as an ornamental plant and high-value food. We aimed to characterize the major repetitive DNA motifs, understand their structure and identify key features that are shaped by genome and sequence evolution.

METHODS

Graph-based clustering with RepeatExplorer was used to identify and classify repetitive motifs in 2.14 millions of 250-bp paired-end Illumina reads from total genomic DNA of C. nankingense. Independently, the frequency of all canonical motifs k-bases long was counted in the raw read data and abundant k-mers (16, 21, 32, 64 and 128) were extracted and assembled to generate longer contigs for repetitive motif identification. For comparison, long terminal repeat retrotransposons were checked in the published C. nankingense reference genome. Fluorescent in situ hybridization was performed to show the chromosomal distribution of the main types of repetitive motifs.

KEY RESULTS

Apart from rDNA (0.86 % of the total genome), a few microsatellites (0.16 %), and telomeric sequences, no highly abundant tandem repeats were identified. There were many transposable elements: 40 % of the genome had sequences with recognizable domains related to transposable elements. Long terminal repeat retrotransposons showed widespread distribution over chromosomes, although different sequence families had characteristic features such as abundance at or exclusion from centromeric or subtelomeric regions. Another group of very abundant repetitive motifs, including those most identified as low-complexity sequences (9.07 %) in the genome, showed no similarity to known sequence motifs or tandemly repeated elements.

CONCLUSIONS

The Chrysanthemum genome has an unusual structure with a very low proportion of tandemly repeated sequences (~1.02 %) in the genome, and a high proportion of low-complexity sequences, most likely degenerated remains of transposable elements. Identifying the presence, nature and genomic organization of major genome fractions enables inference of the evolutionary history of sequences, including degeneration and loss, critical to understanding biodiversity and diversification processes in the genomes of diploid and polyploid Chrysanthemum, Asteraceae and plants more widely.

Adaptive Evolution of Rhizobial Symbiosis beyond Horizontal Gene Transfer: From Genome Innovation to Regulation Reconstruction

There are ubiquitous variations in symbiotic performance of different rhizobial strains associated with the same legume host in agricultural practices. This is due to polymorphisms of symbiosis genes and/or largely unexplored variations in integration efficiency of symbiotic function. Here, we reviewed cumulative evidence on integration mechanisms of symbiosis genes. Experimental evolution, in concert with reverse genetic studies based on pangenomics, suggests that gain of the same circuit of key symbiosis genes through horizontal gene transfer is necessary but sometimes insufficient for bacteria to establish an effective symbiosis with legumes. An intact genomic background of the recipient may not support the proper expression or functioning of newly acquired key symbiosis genes. Further adaptive evolution, through genome innovation and reconstruction of regulation networks, may confer the recipient of nascent nodulation and nitrogen fixation ability. Other accessory genes, either co-transferred with key symbiosis genes or stochastically transferred, may provide the recipient with additional adaptability in ever-fluctuating host and soil niches. Successful integrations of these accessory genes with the rewired core network, regarding both symbiotic and edaphic fitness, can optimize symbiotic efficiency in various natural and agricultural ecosystems. This progress also sheds light on the development of elite rhizobial inoculants using synthetic biology procedures.

Local adaptation: Causal agents of selection and adaptive trait divergence

Divergent selection across the landscape can favor the evolution of local adaptation in populations experiencing contrasting conditions. Local adaptation is widely observed in a diversity of taxa, yet we have a surprisingly limited understanding of the mechanisms that give rise to it. For instance, few have experimentally confirmed the biotic and abiotic variables that promote local adaptation, and fewer yet have identified the phenotypic targets of selection that mediate local adaptation. Here, we highlight critical gaps in our understanding of the process of local adaptation and discuss insights emerging from in-depth investigations of the agents of selection that drive local adaptation, the phenotypes they target, and the genetic basis of these phenotypes. We review historical and contemporary methods for assessing local adaptation, explore whether local adaptation manifests differently across life history, and evaluate constraints on local adaptation.

The Boechera model system for evolutionary ecology

Model systems in biology expand the research capacity of individuals and the community. Closely related to Arabidopsis, the genus Boechera has emerged as an important ecological model owing to the ability to integrate across molecular, functional, and eco-evolutionary approaches. Boechera species are broadly distributed in relatively undisturbed habitats predominantly in western North America and provide one of the few experimental systems for identification of ecologically important genes through genome-wide association studies and investigations of selection with plants in their native habitats. The ecologically, evolutionarily, and agriculturally important trait of apomixis (asexual reproduction via seeds) is common in the genus, and field experiments suggest that abiotic and biotic environments shape the evolution of sex. To date, population genetic studies have focused on the widespread species B. stricta, detailing population divergence and demographic history. Molecular and ecological studies show that balancing selection maintains genetic variation in ~10% of the genome, and ecological trade-offs contribute to complex trait variation for herbivore resistance, flowering phenology, and drought tolerance. Microbiome analyses have shown that host genotypes influence leaf and root microbiome composition, and the soil microbiome influences flowering phenology and natural selection. Furthermore, Boechera offers numerous opportunities for investigating biological responses to global change. In B. stricta, climate change has induced a shift of >2 weeks in the timing of first flowering since the 1970s, altered patterns of natural selection, generated maladaptation in previously locally-adapted populations, and disrupted life history trade-offs. Here we review resources and results for this eco-evolutionary model system and discuss future research directions.

Polygenic adaptation contributes to the invasive success of the Colorado potato beetle

How invasive species cope with novel selective pressures with limited genetic variation is a fundamental question in molecular ecology. Several mechanisms have been proposed, but they can lack generality. Here, we addressed an alternative solution, polygenic adaptation, wherein traits that arise from multiple combinations of loci may be less sensitive to loss of variation during invasion. We tested the polygenic signal of environmental adaptation of Colorado potato beetle (CPB) introduced in Eurasia. Population genomic analyses showed declining genetic diversity in the eastward expansion of Eurasian populations, and weak population genetic structure (except for the invasion fronts in Asia). Demographic history showed that all populations shared a strong bottleneck about 100 years ago when CPB was introduced to Europe. Genome scans revealed a suite of genes involved in activity regulation functions that are plausibly related to cold stress, including some well-founded functions (e.g., the activity of phosphodiesterase, the G-protein regulator) and discrete functions. Such polygenic architecture supports the hypothesis that polygenic adaptation and potentially genetic redundancy can fuel the adaptation of CPB despite strong genetic depletion, thus representing a promising general mechanism for resolving the genetic paradox of invasion. More broadly, most complex traits based on polygenes may be less sensitive to invasive bottlenecks, thus ensuring the evolutionary success of invasive species in novel environments.

Pulmonate slug evolution is reflected in the de novo genome of Arion vulgaris Moquin-Tandon, 1855

Stylommatophoran pulmonate land slugs and snails successfully completed the water-to-land transition from an aquatic ancestor and flourished on land. Of the 30,000 estimated species, very few genomes have so far been published. Here, we assembled and characterized a chromosome-level genome of the “Spanish” slug, Arion vulgaris Moquin-Tandon, 1855, a notorious pest land slug in Europe. Using this reference genome, we conclude that a whole-genome duplication event occurred approximately 93–109 Mya at the base of Stylommatophora and might have promoted land invasion and adaptive radiation. Comparative genomic analyses reveal that genes related to the development of kidney, blood vessels, muscle, and nervous systems had expanded in the last common ancestor of land pulmonates, likely an evolutionary response to the terrestrial challenges of gravity and water loss. Analyses of A. vulgaris gene families and positively selected genes show the slug has evolved a stronger ability to counteract the greater threats of external damage, radiation, and water loss lacking a protective shell. Furthermore, a recent burst of long interspersed elements in the genome of A. vulgaris might affect gene regulation and contribute to rapid phenotype changes in A. vulgaris, which might be conducive to its rapid adaptation and invasiveness.

Function deficiency of GhOMT1 causes anthocyanidins over-accumulation and diversifies fibre colours in cotton (Gossypium hirsutum)

Naturally coloured cotton (NCC) fibres need little or no dyeing process in textile industry to low-carbon emission and are environment-friendly. Proanthocyanidins (PAs) and their derivatives were considered as the main components causing fibre coloration and made NCCs very popular and healthy, but the monotonous fibre colours greatly limit the wide application of NCCs. Here a G. hirsutum empurpled mutant (HS2) caused by T-DNA insertion is found to enhance the anthocyanidins biosynthesis and accumulate anthocyanidins in the whole plant. HPLC and LC/MS-ESI analysis confirmed the anthocyanidins methylation and peonidin, petunidin and malvidin formation are blocked. The deficiency of GhOMT1 in HS2 was associated with the activation of the anthocyanidin biosynthesis and the altered components of anthocyanidins. The transcripts of key genes in anthocyanidin biosynthesis pathway are significantly up-regulated in HS2, while transcripts of the genes for transport and decoration were at similar levels as in WT. To investigate the potential mechanism of GhOMT1 deficiency in cotton fibre coloration, HS2 mutant was crossed with NCCs. Surprisingly, offsprings of HS2 and NCCs enhanced PAs biosynthesis and increased PAs levels in their fibres from the accumulated anthocyanidins through up-regulated GhANR and GhLAR. As expected, multiple novel lines with improved fibre colours including orange red and navy blue were produced in their generations. Based on this work, a new strategy for breeding diversified NCCs was brought out by promoting PA biosynthesis. This work will help shed light on mechanisms of PA biosynthesis and bring out potential molecular breeding strategy to increase PA levels in NCCs.

Transposable element insertion: a hidden major source of domesticated phenotypic variation in Brassica rapa

Transposable element (TE) is prevalent in plant genomes. However, studies on their impact on phenotypic evolution in crop plants are relatively rare, because systematically identifying TE insertions within a species has been a challenge. Here, we present a novel approach for uncovering TE insertion polymorphisms (TIPs) using pan-genome analysis combined with population-scale resequencing, and we adopt this pipeline to retrieve TIPs in a Brassica rapa germplasm collection. We found that 23% of genes within the reference Chiifu-401-42 genome harbored TIPs. TIPs tended to have large transcriptional effects, including modifying gene expression levels and altering gene structure by introducing new introns. Among 524 diverse accessions, TIPs broadly influenced genes related to traits and acted a crucial role in the domestication of B. rapa morphotypes. As examples, four specific TIP-containing genes were found to be candidates that potentially involved in various climatic conditions, promoting the formation of diverse vegetable crops in B. rapa. Our work reveals the hitherto hidden TIPs implicated in agronomic traits and highlights their widespread utility in studies of crop domestication.

Chromosome-level genome assembly of the aquatic plant Nymphoides indica reveals transposable element bursts and NBS-LRR gene family expansion shedding light on its invasiveness

Nymphoides indica, an aquatic plant, is an invasive species that causes both ecological and economic damage in North America and elsewhere. However, the lack of genomic data of N. indica limits the in-depth analysis of this invasive species. Here, we report a chromosome-level genome assembly of nine pseudochromosomes of N. indica with a total size of ∼ 520 Mb. More than half of the N. indica genome consists of transposable elements (TEs), and a higher density of TEs around genes may play a significant role in response to an ever-changing environment by regulating the nearby gene. Additionally, our analysis revealed that N. indica only experienced a gamma (γ) whole-genome triplication event. Functional enrichment of the N. indica-specific and expanded gene families highlighted genes involved in the responses to hypoxia and plant–pathogen interactions, which may strengthen the ability to adapt to external challenges and improve ecological fitness. Furthermore, we identified 160 members of the nucleotide-binding site and leucine-rich repeat gene family, which may be linked to the defence response. Collectively, the high-quality N. indica genome reported here opens a novel avenue to understand the evolution and rapid invasion of Nymphoides spp.

Differential MC5R loss in whales and manatees reveals convergent evolution to the marine environment

Melanocortin 5 receptor (MC5R), which is expressed in the terminally differentiated sebaceous gland, is a G protein-coupled receptor (GPCR). MC5R exists mostly in mammals but is completely lost in whales; only the relic of MC5R can be detected in manatees, and phenotypically, they have lost sebaceous glands. Interestingly, whales and manatees are both aquatic mammals but have no immediate common ancestors. The loss of MC5R and sebaceous glands in whales and manatees is likely to be a result of convergent evolution. Here, we find that MC5R in whales and manatees are lost by two different mechanisms. Homologous recombination of MC5R in manatees and the insertion of reverse transcriptase in whales lead to the gene loss, respectively. On one hand, in manatees, there are two "TTATC" sequences flanking MC5R, and homologous recombination of the segments between the two "TTATC" sequences resulted in the partial loss of the sequence of MC5R. On the other hand, in whales, reverse transcriptase inserts between MC2R and RNMT on the chromosome led to the loss of MC5R. Based on these two different mechanisms for gene loss in whales and manatees, we finally concluded that MC5R loss might be the result of convergent evolution to the marine environment, and we explored the impact on biological function that is significant to environmental adaptation.

Rosa1, a Transposable Element-Like Insertion, Produces Red Petal Coloration in Rose Through Altering RcMYB114 Transcription

Rose (Rosa sp.) flowers have a rich diversity of colors resulting from the differential accumulation of anthocyanins, flavonols, and carotenoids. However, the genetic and molecular determinants of the red-petal trait in roses remains poorly understood. Here we report that a transposable element-like insertion (Rosa1) into RcMYB114, a R2R3-MYB transcription factor's promoter region causes its transcription, resulting in red petals. In red-petal varieties, RcMYB114 is expressed specifically in flower organs, but is absent from non-red varieties. Sequencing, yeast two-hybrid, transient transformation, and promoter activity assays of RcMYB114 independently confirmed the role of Rosa1 in altering RcMYB114's transcription and downstream effects on flower color. Genetic and molecular evidence confirmed that the Rosa1 transposable element-like insertion, which is a previously unknown DNA transposable element, is different from those in other plants and is a reliable molecular marker to screen red-petal roses.

Transposable Elements: Distribution, Polymorphism, and Climate Adaptation in Populus

Transposable elements (TEs) are a class of mobile genetic elements that make effects on shaping rapid phenotypic traits of adaptive significance. TE insertions are usually related to transcription changes of nearby genes, and thus may be subjected to purifying selection. Based on the available genome resources of Populus, we found that the composition of Helitron DNA family were highly variable and could directly influence the transcription of nearby gene expression, which are involving in stress-responsive, programmed cell death, and apoptosis pathway. Next, we indicated TEs are highly enriched in Populus trichocarpa compared with three other congeneric poplar species, especially located at untranslated regions (3'UTRs and 5'UTRs) and Helitron transposons, particularly 24-nt siRNA-targeted, are significantly associated with reduced gene expression. Additionally, we scanned a representative resequenced Populus tomentosa population, and identified 9,680 polymorphic TEs loci. More importantly, we identified a Helitron transposon located at the 3'UTR, which could reduce WRKY18 expression level. Our results highlight the importance of TE insertion events, which could regulate gene expression and drive adaptive phenotypic variation in Populus.

The important contribution of transposable elements to phenotypic variation and evolution

Transposable elements (TEs) are responsible for significant genomic variation in plants. Our understanding of the evolutionary forces shaping TE polymorphism has lagged behind other mutations because of the difficulty of accurately identifying TE polymorphism in short-read population genomic data. However, new approaches allow us to quantify TE polymorphisms in population datasets and address fundamental questions about the evolution of these polymorphisms. Here, we discuss how insertional biases shape where, when, and how often TEs insert throughout the genome. Next, we examine mechanisms by which TEs can affect phenotype. Finally, we evaluate current evidence for selection on TE polymorphisms. All together, it is clear that TEs are important, but underappreciated, contributors to intraspecific phenotypic variation, and that understanding the dynamics governing TE polymorphism is crucial for evolutionary biologists interested in the maintenance of variation.

A Wheat TaTOE1-B1 Transcript TaTOE1-B1-3 Can Delay the Flowering Time of Transgenic Arabidopsis

Flowering time is one of the most important agronomic traits in wheat production. A proper flowering time might contribute to the reduction or avoidance of biotic and abiotic stresses, adjust plant architecture, and affect the yield and quality of grain. In this study, TaTOE1-B1 in wheat produced three transcripts (TaTOE1-B1-1, TaTOE1-B1-2, and TaTOE1-B1-3) by alternative splicing. Compared to the longest transcript, TaTOE1-B1-1, TaTOE1-B1-3 has a deletion in the sixth exon (1219-1264 bp). Under long-day conditions, the heterologous overexpression of the TaTOE1-B1-3 gene delayed flowering, prolonged the vegetative growth time, and enlarged the vegetative body of Arabidopsis, but that of TaTOE1-B1-1 did not. As typical AP2 family members, TaTOE1-B1-1 and TaTOE1-B1-3 are mainly located in the nucleus and have transcriptional activation activities; the transcriptional activation region of TaTOE1-B1-3 is located in the C-terminal. In TaTOE1-B1-3 overexpression lines, the expression of flowering-related AtFT and AtSOC1 genes is significantly downregulated. In addition, this study confirms the protein-protein interaction between TaTOE1-B1-3 and TaPIFI, which may play an important role in flowering inhibition. These results provide a theoretical basis for the precise regulation of wheat flowering time.

High genetic and epigenetic variation of transposable elements: Potential drivers to rapid adaptive evolution for the noxious invasive weed Mikania micrantha

Why invasive species can rapidly adapt to novel environments is a puzzling question known as the genetic paradox of invasive species. This paradox is explainable in terms of transposable elements (TEs) activity, which are theorized to be powerful mutational forces to create genetic variation. Mikania micrantha, a noxious invasive weed, in this sense provides an excellent opportunity to test the explanation. The genetic and epigenetic variation of 21 invasive populations of M. micrantha in southern China have been examined by using transposon display (TD) and transposon methylation display (TMD) techniques to survey 12 TE superfamilies. Our results showed that M. micrantha populations maintained an almost equally high level of TE-based genetic and epigenetic variation and they have been differentiated into subpopulations genetically and epigenetically. A similar positive spatial genetic and epigenetic structure pattern was observed within 300 m. Six and seven TE superfamilies presented significant genetic and epigenetic isolation by distance (IBD) pattern. In total, 59 genetic and 86 epigenetic adaptive TE loci were identified. Of them, 51 genetic and 44 epigenetic loci were found to correlate with 25 environmental variables (including precipitation, temperature, vegetation coverage, and soil metals). Twenty-five transposon-inserted genes were sequenced and homology-based annotated, which are found to be involved in a variety of molecular and cellular functions. Our research consolidates the importance of TE-associated genetic and epigenetic variation in the rapid adaptation and invasion of M. micrantha.

The chromosome-scale genome of Magnolia officinalis provides insight into the evolutionary position of magnoliids

Magnolia officinalis, a representative tall aromatic tree of the Magnoliaceae family, is a medicinal plant that is widely used in diverse industries from medicine to cosmetics. We report a chromosome-scale draft genome of M. officinalis, in which ∼99.66% of the sequences were anchored onto 19 chromosomes with the scaffold N50 of 76.62 Mb. We found that a high proportion of repetitive sequences was a common feature of three Magnoliaceae with known genomic data. Magnoliids were a sister clade to eudicots-monocots, which provided more support for understanding the phylogenetic position among angiosperms. An ancient duplication event occurred in the genome of M. officinalis and was shared with Lauraceae. Based on RNA-seq analysis, we identified several key enzyme-coding gene families associated with the biosynthesis of lignans in the genome. The construction of the M. officinalis genome sequence will serve as a reference for further studies of Magnolia, as well as other Magnoliaceae.

Comparative Genome Analyses Highlight Transposon-Mediated Genome Expansion and the Evolutionary Architecture of 3D Genomic Folding in Cotton

Transposable element (TE) amplification has been recognized as a driving force mediating genome size expansion and evolution, but the consequences for shaping 3D genomic architecture remains largely unknown in plants. Here, we report reference-grade genome assemblies for three species of cotton ranging 3-fold in genome size, namely Gossypium rotundifolium (K2), G. arboreum (A2), and G. raimondii (D5), using Oxford Nanopore Technologies. Comparative genome analyses document the details of lineage-specific TE amplification contributing to the large genome size differences (K2, 2.44 Gb; A2, 1.62 Gb; D5, 750.19 Mb) and indicate relatively conserved gene content and synteny relationships among genomes. We found that approximately 17% of syntenic genes exhibit chromatin status change between active ("A") and inactive ("B") compartments, and TE amplification was associated with the increase of the proportion of A compartment in gene regions (∼7,000 genes) in K2 and A2 relative to D5. Only 42% of topologically associating domain (TAD) boundaries were conserved among the three genomes. Our data implicate recent amplification of TEs following the formation of lineage-specific TAD boundaries. This study sheds light on the role of transposon-mediated genome expansion in the evolution of higher-order chromatin structure in plants.

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 elements and introgression introduce genetic variation in the invasive ant Cardiocondyla obscurior

Introduced populations of invasive organisms have to cope with novel environmental challenges, while having reduced genetic variation caused by founder effects. The mechanisms associated with this "genetic paradox of invasive species" has received considerable attention, yet few studies have examined the genomic architecture of invasive species. Populations of the heart node ant Cardiocondyla obscurior belong to two distinct lineages, a New World lineage so far only found in Latin America and a more globally distributed Old World lineage. In the present study, we use population genomic approaches to compare populations of the two lineages with apparent divergent invasive potential. We find that the strong genetic differentiation of the two lineages began at least 40,000 generations ago and that activity of transposable elements (TEs) has contributed significantly to the divergence of both lineages, possibly linked to the very unusual genomic distribution of TEs in this species. Furthermore, we show that introgression from the Old World lineage is a dominant source of genetic diversity in the New World lineage, despite the lineages' strong genetic differentiation. Our study uncovers mechanisms underlying novel genetic variation in introduced populations of C. obscurior that could contribute to the species' adaptive potential.

The chromosome-level reference genome assembly for Dendrobium officinale and its utility of functional genomics research and molecular breeding study

Dendrobium officinale, an important medicinal plant of the genus Dendrobium in Orchidaceae family, has been used as traditional Chinese medicine (TCM) for nearly thousands of years. Here, we report the first chromosome-level reference genome of D. officinale, based on PacBio long-reads, Illumina short-reads and Hi-C data. The high-quality assembled genome is 1.23 Gb long, with contig N50 of 1.44 Mb. A total of 93.53% genome sequences were assembled into 19 pseudochromosomes with a super scaffold N50 of 63.07 Mb. Through comparative genomic analysis, we explored the expanded gene families of D. officinale, and also their impact on environmental adaptation and biosynthesis of secondary metabolites. We further performed detailed transcriptional analysis of D. officinale, and identified the candidate genes involved in the biosynthesis of three main active ingredients, including polysaccharides, alkaloids and flavonoids. In addition, the MODIFYING WALL LIGNIN-1 (MWL1) gene, which inferred from Genome-Wide Association Studies (GWAS) based on the resequencing date from D. officinale and five related species and their morphologic features, may contribute to the plant production (yield of stems) of D. officinale. Therefore, the high-quality reference genome reported in this study could benefits functional genomics research and molecular breeding of D. officinale.

The Structural, Functional and Evolutionary Impact of Transposable Elements in Eukaryotes

Transposable elements (TEs) are nearly ubiquitous in eukaryotes. The increase in genomic data, as well as progress in genome annotation and molecular biology techniques, have revealed the vast number of ways mobile elements have impacted the evolution of eukaryotes. In addition to being the main cause of difference in haploid genome size, TEs have affected the overall organization of genomes by accumulating preferentially in some genomic regions, by causing structural rearrangements or by modifying the recombination rate. Although the vast majority of insertions is neutral or deleterious, TEs have been an important source of evolutionary novelties and have played a determinant role in the evolution of fundamental biological processes. TEs have been recruited in the regulation of host genes and are implicated in the evolution of regulatory networks. They have also served as a source of protein-coding sequences or even entire genes. The impact of TEs on eukaryotic evolution is only now being fully appreciated and the role they may play in a number of biological processes, such as speciation and adaptation, remains to be deciphered.

Epigenetics and the success of invasive plants

Biological invasions impose ecological and economic problems on a global scale, but also provide extraordinary opportunities for studying contemporary evolution. It is critical to understand the evolutionary processes that underly invasion success in order to successfully manage existing invaders, and to prevent future invasions. As successful invasive species sometimes are suspected to rapidly adjust to their new environments in spite of very low genetic diversity, we are obliged to re-evaluate genomic-level processes that translate into phenotypic diversity. In this paper, we review work that supports the idea that trait variation, within and among invasive populations, can be created through epigenetic or other non-genetic processes, particularly in clonal invaders where somatic changes can persist indefinitely. We consider several processes that have been implicated as adaptive in invasion success, focusing on various forms of 'genomic shock' resulting from exposure to environmental stress, hybridization and whole-genome duplication (polyploidy), and leading to various patterns of gene expression re-programming and epigenetic changes that contribute to phenotypic variation or even novelty. These mechanisms can contribute to transgressive phenotypes, including hybrid vigour and novel traits, and may thus help to understand the huge successes of some plant invaders, especially those that are genetically impoverished. This article is part of the theme issue 'How does epigenetics influence the course of evolution?'

A Global Landscape of Miniature Inverted-Repeat Transposable Elements in the Carrot Genome

Miniature inverted-repeat transposable elements (MITEs) are the most abundant group of Class II mobile elements in plant genomes. Their presence in genic regions may alter gene structure and expression, providing a new source of functional diversity. Owing to their small size and lack of coding capacity, the identification of MITEs has been demanding. However, the increasing availability of reference genomes and bioinformatic tools provides better means for the genome-wide identification and analysis of MITEs and for the elucidation of their contribution to the evolution of plant genomes. We mined MITEs in the carrot reference genome DH1 using MITE-hunter and developed a curated carrot MITE repository comprising 428 families. Of the 31,025 MITE copies spanning 10.34 Mbp of the carrot genome, 54% were positioned in genic regions. Stowaways and Tourists were frequently present in the vicinity of genes, while Mutator-like MITEs were relatively more enriched in introns. hAT-like MITEs were relatively more frequently associated with transcribed regions, including untranslated regions (UTRs). Some carrot MITE families were shared with other Apiaceae species. We showed that hAT-like MITEs were involved in the formation of new splice variants of insertion-harboring genes. Thus, carrot MITEs contributed to the accretion of new diversity by altering transcripts and possibly affecting the regulation of many genes.

Jumpstarting evolution: How transposition can facilitate adaptation to rapid environmental changes

Because of their ability to replicate across genomes, transposable elements (TEs) represent major generators of large-effect mutations. As a result, chromatin-based mechanisms have evolved to control the mutational potential of TEs at multiple levels, from the epigenetic silencing of TE sequences, through the modulation of their integration space, up to the alleviation of the impact of new insertions. Although most TE insertions are highly deleterious, some can provide key adaptive variation. Together with their remarkable sensitivity to the environment and precise integration preferences, the unique characteristics of TEs place them as potent genomic engines of adaptive innovation. Herein, we review recent works exploring the regulation and impact of transposition in nature and discuss their implications for the evolutionary response of species to drastic environmental changes.

Genome-wide DNA mutations in Arabidopsis plants after multigenerational exposure to high temperatures

BACKGROUND

Elevated temperatures can cause physiological, biochemical, and molecular responses in plants that can greatly affect their growth and development. Mutations are the most fundamental force driving biological evolution. However, how long-term elevations in temperature influence the accumulation of mutations in plants remains unknown.

RESULTS

Multigenerational exposure of Arabidopsis MA (mutation accumulation) lines and MA populations to extreme heat and moderate warming results in significantly increased mutation rates in single-nucleotide variants (SNVs) and small indels. We observe distinctive mutational spectra under extreme and moderately elevated temperatures, with significant increases in transition and transversion frequencies. Mutation occurs more frequently in intergenic regions, coding regions, and transposable elements in plants grown under elevated temperatures. At elevated temperatures, more mutations accumulate in genes associated with defense responses, DNA repair, and signaling. Notably, the distribution patterns of mutations among all progeny differ between MA populations and MA lines, suggesting that stronger selection effects occurred in populations. Methylation is observed more frequently at mutation sites, indicating its contribution to the mutation process at elevated temperatures. Mutations occurring within the same genome under elevated temperatures are significantly biased toward low gene density regions, special trinucleotides, tandem repeats, and adjacent simple repeats. Additionally, mutations found in all progeny overlap significantly with genetic variations reported in 1001 Genomes, suggesting non-uniform distribution of de novo mutations through the genome.

CONCLUSION

Collectively, our results suggest that elevated temperatures can accelerate the accumulation, and alter the molecular profiles, of DNA mutations in plants, thus providing significant insight into how environmental temperatures fuel plant evolution.

The worldwide invasion of Drosophila suzukii is accompanied by a large increase of transposable element load and a small number of putatively adaptive insertions

Transposable Elements (TEs) are ubiquitous and mobile repeated sequences. They are major determinants of host fitness. Here, we characterized the TE content of the spotted wing fly Drosophila suzukii. Using a recently improved genome assembly, we reconstructed TE sequences de novo, and found that TEs occupy 47% of the genome and are mostly located in gene poor regions. The majority of TE insertions segregate at low frequencies, indicating a recent and probably ongoing TE activity. To explore TE dynamics in the context of biological invasions, we studied variation of TE abundance in genomic data from 16 invasive and six native populations of D. suzukii. We found a large increase of the TE load in invasive populations correlated with a reduced Watterson estimate of genetic diversity θ̂w a proxy of effective population size. We did not find any correlation between TE contents and bioclimatic variables, indicating a minor effect of environmentally induced TE activity. A genome-wide association study revealed that ca. 2,000 genomic regions are associated with TE abundance. We did not find, however, any evidence in such regions of an enrichment for genes known to interact with TE activity (e.g. transcription factor encoding genes or genes of the piRNA pathway). Finally, the study of TE insertion frequencies revealed 15 putatively adaptive TE insertions, six of them being likely associated with the recent invasion history of the species.

Comparative genome analyses highlight transposon-mediated genome expansion and the evolutionary architecture of 3D genomic folding in cotton

Transposable element (TE) amplification has been recognized as a driving force mediating genome size expansion and evolution, but the consequences for shaping 3D genomic architecture remains largely unknown in plants. Here, we report reference-grade genome assemblies for three species of cotton ranging 3-fold in genome size, namely Gossypium rotundifolium (K₂), G. arboreum (A₂), and G. raimondii (D₅), using Oxford Nanopore Technologies. Comparative genome analyses document the details of lineage-specific TE amplification contributing to the large genome size differences (K₂, 2.44 Gb; A₂, 1.62 Gb; D₅, 750.19 Mb) and indicate relatively conserved gene content and synteny relationships among genomes. We found that approximately 17% of syntenic genes exhibit chromatin status change between active (“A”) and inactive (“B”) compartments, and TE amplification was associated with the increase of the proportion of A compartment in gene regions (∼7,000 genes) in K₂ and A₂ relative to D₅. Only 42% of topologically associating domain (TAD) boundaries were conserved among the three genomes. Our data implicate recent amplification of TEs following the formation of lineage-specific TAD boundaries. This study sheds light on the role of transposon-mediated genome expansion in the evolution of higher-order chromatin structure in plants.

The transposable elements of the Drosophila serrata reference panel

Transposable elements (TEs) are an important component of the complex genomic ecosystem. Understanding the tempo and mode of TE proliferation, that is whether it is in maintained in transposition selection balance, or is induced periodically by environmental stress or other factors, is important for understanding the evolution of organismal genomes through time. Although TEs have been characterized in individuals or limited samples, a true understanding of the population genetics of TEs, and therefore the tempo and mode of transposition, is still lacking. Here, we characterize the TE landscape in an important model Drosophila, Drosophila serrata using the D. serrata reference panel, which is comprised of 102 sequenced inbred genotypes. We annotate the families of TEs in the D. serrata genome and investigate variation in TE copy number between genotypes. We find that many TEs have low copy number in the population, but this varies by family and includes a single TE making up to 50% of the genome content of TEs. We find that some TEs proliferate in particular genotypes compared with population levels. In addition, we characterize variation in each TE family allowing copy number to vary in each genotype and find that some TEs have diversified very little between individuals suggesting recent spread. TEs are important sources of spontaneous mutations in Drosophila, making up a large fraction of the total number of mutations in particular genotypes. Understanding the dynamics of TEs within populations will be an important step toward characterizing the origin of variation within and between species.

The URL1-ROC5-TPL2 transcriptional repressor complex represses the ACL1 gene to modulate leaf rolling in rice

Moderate leaf rolling is beneficial for leaf erectness and compact plant architecture. However, our understanding regarding the molecular mechanisms of leaf rolling is still limited. Here, we characterized a semi-dominant rice (Oryza sativa L.) mutant upward rolled leaf 1 (Url1) showing adaxially rolled leaves due to a decrease in the number and size of bulliform cells. Map-based cloning revealed that URL1 encodes the homeodomain-leucine zipper (HD-Zip) IV family member RICE OUTERMOST CELL-SPECIFIC 8 (ROC8). A single-base substitution in one of the two conserved complementary motifs unique to the 3'-untranslated region of this family enhanced URL1 mRNA stability and abundance in the Url1 mutant. URL1 (UPWARD ROLLED LEAF1) contains an ethylene-responsive element binding factor-associated amphiphilic repression motif and functions as a transcriptional repressor via interaction with the TOPLESS co-repressor OsTPL2. Rather than homodimerizing, URL1 heterodimerizes with another HD-ZIP IV member ROC5. URL1 could bind directly to the promoter and suppress the expression of abaxially curled leaf 1 (ACL1), a positive regulator of bulliform cell development. Knockout of OsTPL2 or ROC5 or overexpression of ACL1 in the Url1 mutant partially suppressed the leaf-rolling phenotype. Our results reveal a regulatory network whereby a transcriptional repression complex composed of URL1, ROC5, and the transcriptional corepressor TPL2 suppresses the expression of the ACL1 gene, thus modulating bulliform cell development and leaf rolling in rice.

RepeatProfiler: A pipeline for visualization and comparative analysis of repetitive DNA profiles

Study of repetitive DNA elements in model organisms highlights the role of repetitive elements (REs) in many processes that drive genome evolution and phenotypic change. Because REs are much more dynamic than single-copy DNA, repetitive sequences can reveal signals of evolutionary history over short time scales that may not be evident in sequences from slower-evolving genomic regions. Many tools for studying REs are directed toward organisms with existing genomic resources, including genome assemblies and repeat libraries. However, signals in repeat variation may prove especially valuable in disentangling evolutionary histories in diverse non-model groups, for which genomic resources are limited. Here, we introduce RepeatProfiler, a tool for generating, visualizing, and comparing repetitive element DNA profiles from low-coverage, short-read sequence data. RepeatProfiler automates the generation and visualization of RE coverage depth profiles (RE profiles) and allows for statistical comparison of profile shape across samples. In addition, RepeatProfiler facilitates comparison of profiles by extracting signal from sequence variants across profiles which can then be analysed as molecular morphological characters using phylogenetic analysis. We validate RepeatProfiler with data sets from ground beetles (Bembidion), flies (Drosophila), and tomatoes (Solanum). We highlight the potential of RE profiles as a high-resolution data source for studies in species delimitation, comparative genomics, and repeat biology.