How important are transposons for plant evolution?

Going broad and deep: sequencing-driven insights into plant physiology, evolution, and crop domestication

Deep sequencing is a term that has become embedded in the plant genomic literature in recent years and with good reason. A torrent of (largely) high-quality genomic and transcriptomic data has been collected and most of this has been publicly released. Indeed, almost 1000 plant genomes have been reported (www.plabipd.de) and the 2000 Plant Transcriptomes Project has long been completed. The EarthBioGenome project will dwarf even these milestones. That said, massive progress in understanding plant physiology, evolution, and crop domestication has been made by sequencing broadly (across a species) as well as deeply (within a single individual). We will outline the current state of the art in genome and transcriptome sequencing before we briefly review the most visible of these broad approaches, namely genome-wide association and transcriptome-wide association studies, as well as the compilation of pangenomes. This will include both (i) the most commonly used methods reliant on single nucleotide polymorphisms and short InDels and (ii) more recent examples which consider structural variants. We will subsequently present case studies exemplifying how their application has brought insight into either plant physiology or evolution and crop domestication. Finally, we will provide conclusions and an outlook as to the perspective for the extension of such approaches to different species, tissues, and biological processes.

Polymorphic inverted repeats near coding genes impact chromatin topology and phenotypic traits in Arabidopsis thaliana

Transposons are mobile elements that are commonly silenced to protect eukaryotic genome integrity. In plants, transposable element (TE)-derived inverted repeats (IRs) are commonly found near genes, where they affect host gene expression. However, the molecular mechanisms of such regulation are unclear in most cases. Expression of these IRs is associated with production of 24-nt small RNAs, methylation of the IRs, and drastic changes in local 3D chromatin organization. Notably, many of these IRs differ between Arabidopsis thaliana accessions, causing variation in short-range chromatin interactions and gene expression. CRISPR-Cas9-mediated disruption of two IRs leads to a switch in genome topology and gene expression with phenotypic consequences. Our data show that insertion of an IR near a gene provides an anchor point for chromatin interactions that profoundly impact the activity of neighboring loci. This turns IRs into powerful evolutionary agents that can contribute to rapid adaptation.

Pan-3D genome analysis reveals structural and functional differentiation of soybean genomes

BACKGROUND

High-order chromatin structure plays important roles in gene regulation. However, the diversity of the three-dimensional (3D) genome across plant accessions are seldom reported.

RESULTS

Here, we perform the pan-3D genome analysis using Hi-C sequencing data from 27 soybean accessions and comprehensively investigate the relationships between 3D genomic variations and structural variations (SVs) as well as gene expression. We find that intersection regions between A/B compartments largely contribute to compartment divergence. Topologically associating domain (TAD) boundaries in A compartments exhibit significantly higher density compared to those in B compartments. Pan-3D genome analysis shows that core TAD boundaries have the highest transcription start site (TSS) density and lowest GC content and repeat percentage. Further investigation shows that non-long terminal repeat (non-LTR) retrotransposons play important roles in maintaining TAD boundaries, while Gypsy elements and satellite repeats are associated with private TAD boundaries. Moreover, presence and absence variation (PAV) is found to be the major contributor to 3D genome variations. Nevertheless, approximately 55% of 3D genome variations are not associated with obvious genetic variations, and half of them affect the flanking gene expression. In addition, we find that the 3D genome may also undergo selection during soybean domestication.

CONCLUSION

Our study sheds light on the role of 3D genomes in plant genetic diversity and provides a valuable resource for studying gene regulation and genome evolution.

DDM1-Mediated TE Silencing in Plants

Epigenetic modifications are indispensable for regulating gene bodies and TE silencing. DECREASE IN DNA METHYLATION 1 (DDM1) is a chromatin remodeller involved in histone modifications and DNA methylation. Apart from maintaining the epigenome, DDM1 also maintains key plant traits such as flowering time and heterosis. The role of DDM1 in epigenetic regulation is best characterised in plants, especially arabidopsis, rice, maize and tomato. The epigenetic changes induced by DDM1 establish the stable inheritance of many plant traits for at least eight generations, yet DDM1 does not methylate protein-coding genes. The DDM1 TE silencing mechanism is distinct and has evolved independently of other silencing pathways. Unlike the RNA-directed DNA Methylation (RdDM) pathway, DDM1 does not depend on siRNAs to enforce the heterochromatic state of TEs. Here, we review DDM1 TE silencing activity in the RdDM and non-RdDM contexts. The DDM1 TE silencing machinery is strongly associated with the histone linker H1 and histone H2A.W. While the linker histone H1 excludes the RdDM factors from methylating the heterochromatin, the histone H2A.W variant prevents TE mobility. The DDM1-H2A.W strategy alone silences nearly all the mobile TEs in the arabidopsis genome. Thus, the DDM1-directed TE silencing essentially preserves heterochromatic features and abolishes mobile threats to genome stability.

Conserved noncoding sequences and de novo Mutator insertion alleles are imprinted in maize

Genomic imprinting is an epigenetic phenomenon in which differential allele expression occurs in a parent-of-origin-dependent manner. Imprinting in plants is tightly linked to transposable elements (TEs), and it has been hypothesized that genomic imprinting may be a consequence of demethylation of TEs. Here, we performed high-throughput sequencing of ribonucleic acids from four maize (Zea mays) endosperms that segregated newly silenced Mutator (Mu) transposons and identified 110 paternally expressed imprinted genes (PEGs) and 139 maternally expressed imprinted genes (MEGs). Additionally, two potentially novel paternally suppressed MEGs are associated with de novo Mu insertions. In addition, we find evidence for parent-of-origin effects on expression of 407 conserved noncoding sequences (CNSs) in maize endosperm. The imprinted CNSs are largely localized within genic regions and near genes, but the imprinting status of the CNSs are largely independent of their associated genes. Both imprinted CNSs and PEGs have been subject to relaxed selection. However, our data suggest that although MEGs were already subject to a higher mutation rate prior to their being imprinted, imprinting may be the cause of the relaxed selection of PEGs. In addition, although DNA methylation is lower in the maternal alleles of both the maternally and paternally expressed CNSs (mat and pat CNSs), the difference between the two alleles in H3K27me3 levels was only observed in pat CNSs. Together, our findings point to the importance of both transposons and CNSs in genomic imprinting in maize.

Plant pan-genomics and its applications

Plant genomes are so highly diverse that a substantial proportion of genomic sequences are not shared among individuals. The variable DNA sequences, along with the conserved core sequences, compose the more sophisticated pan-genome that represents the collection of all non-redundant DNA in a species. With rapid progress in genome sequencing technologies, pan-genome research in plants is now accelerating. Here we review recent advances in plant pan-genomics, including major driving forces of structural variations that constitute the variable sequences, methodological innovations for representing the pan-genome, and major successes in constructing plant pan-genomes. We also summarize recent efforts toward decoding the remaining dark matter in telomere-to-telomere or gapless plant genomes. These new genome resources, which have remarkable advantages over numerous previously assembled less-than-perfect genomes, are expected to become new references for genetic studies and plant breeding.

The role of LTR retrotransposons in plant genetic engineering: how to control their transposition in the genome

We briefly discuss that the similarity of LTR retrotransposons to retroviruses is a great opportunity for the development of a genetic engineering tool that exploits intragenic elements in the plant genome for plant genetic improvement. Long terminal repeat (LTR) retrotransposons are very similar to retroviruses but do not have the property of being infectious. While spreading between its host cells, a retrovirus inserts a DNA copy of its genome into the cells. The ability of retroviruses to cause infection with genome integration allows genes to be delivered to cells and tissues. Retrovirus vectors are, however, only specific to animals and insects, and, thus, are not relevant to plant genetic engineering. However, the similarity of LTR retrotransposons to retroviruses is an opportunity to explore the former as a tool for genetic engineering. Although recent long-read sequencing technologies have advanced the knowledge about transposable elements (TEs), the integration of TEs is still unable either to control them or to direct them to specific genomic locations. The use of existing intragenic elements to achieve the desired genome composition is better than using artificial constructs like vectors, but it is not yet clear how to control the process. Moreover, most LTR retrotransposons are inactive and unable to produce complete proteins. They are also highly mutable. In addition, it is impossible to find a full active copy of a LTR retrotransposon out of thousands of its own copies. Theoretically, if these elements were directly controlled and turned on or off using certain epigenetic mechanisms (inducing by stress or infection), LTR retrotransposons could be a great opportunity to develop a genetic engineering tool using intragenic elements in the plant genome. In this review, the recent developments in uncovering the nature of LTR retrotransposons and the possibility of using these intragenic elements as a tool for plant genetic engineering are briefly discussed.

Graph-based pangenomics maximizes genotyping density and reveals structural impacts on fungal resistance in melon

The genomic sequences segregating in experimental populations are often highly divergent from the community reference and from one another. Such divergence is problematic under various short-read-based genotyping strategies. In addition, large structural differences are often invisible despite being strong candidates for causal variation. These issues are exacerbated in specialty crop breeding programs with fewer, lower-quality sequence resources. Here, we examine the benefits of complete genomic information, based on long-read assemblies, in a biparental mapping experiment segregating at numerous disease resistance loci in the non-model crop, melon (Cucumis melo). We find that a graph-based approach, which uses both parental genomes, results in 19% more variants callable across the population and raw allele calls with a 2 to 3-fold error-rate reduction, even relative to single reference approaches using a parent genome. We show that structural variation has played a substantial role in shaping two Fusarium wilt resistance loci with known causal genes. We also report on the genetics of powdery mildew resistance, where copy number variation and local recombination suppression are directly interpretable via parental genome alignments. Benefits observed, even in this low-resolution biparental experiment, will inevitably be amplified in more complex populations.

Regulatory mechanism of a heat-activated retrotransposon by DDR complex in Arabidopsis thaliana

The RNA-directed DNA methylation (RdDM) pathway plays an essential role in the transposon silencing mechanism; the DDR complex, consisting of DRD1, DMS3, and RDM1, is an essential component of the RdDM pathway. ONSEN, identified in Arabidopsis, is a retrotransposon activated by heat stress at 37°C; however, studies on the regulation of ONSEN are limited. In this study, we analyzed the regulation of ONSEN activity by the DDR complex in Arabidopsis. We elucidated that loss of any component of the DDR complex increased ONSEN transcript levels. Transgenerational transposition of ONSEN was observed in the DDR-complex mutants treated with heat stress for 48 h. Furthermore, the DDR complex components DRD1, DMS3, and RDM1 played independent roles in suppressing ONSEN transcription and transposition. Moreover, we found that the duration of heat stress affects ONSEN activity. Therefore, the results of this study provide new insights into the retrotransposon regulatory mechanisms of the DDR complex in the RdDM pathway.

Transposable Element Tissue-Specific Response to Temperature Stress in the Stenothermal Fish Puntius tetrazona

Ray-finned fish represent a very interesting group of vertebrates comprising a variety of organisms living in different aquatic environments worldwide. In the case of stenothermal fish, thermal fluctuations are poorly tolerated, thus ambient temperature represents a critical factor. In this paper, we considered the tiger barb Puntius tetrazona, a freshwater fish belonging to the family Cyprinidae, living at 21-28 °C. We analyzed the available RNA-Seq data obtained from specimens exposed at 27 °C and 13 °C to investigate the transcriptional activity of transposable elements (TEs) and genes encoding for proteins involved in their silencing in the brain, gill, and liver. TEs are one of the tools generating genetic variability that underlies biological evolution, useful for organisms to adapt to environmental changes. Our findings highlighted a different response of TEs in the three analyzed tissues. While in the brain and gill, no variation in TE transcriptional activity was observed, a remarkable increase at 13 °C was recorded in the liver. Moreover, the transcriptional analysis of genes encoding proteins involved in TE silencing such as heterochromatin formation, the NuRD complex, and the RISC complex (e.g., AGO and GW182 proteins) highlighted their activity in the hepatic tissue. Overall, our findings suggested that this tissue is a target organ for this kind of stress, since TE activation might regulate the expression of stress-induced genes, leading to a better response of the organism to temperature changes. Therefore, this view corroborates once again the idea of a potential role of TEs in organism rapid adaptation, hence representing a promising molecular tool for species resilience.

Inpactor2: a software based on deep learning to identify and classify LTR-retrotransposons in plant genomes

LTR-retrotransposons are the most abundant repeat sequences in plant genomes and play an important role in evolution and biodiversity. Their characterization is of great importance to understand their dynamics. However, the identification and classification of these elements remains a challenge today. Moreover, current software can be relatively slow (from hours to days), sometimes involve a lot of manual work and do not reach satisfactory levels in terms of precision and sensitivity. Here we present Inpactor2, an accurate and fast application that creates LTR-retrotransposon reference libraries in a very short time. Inpactor2 takes an assembled genome as input and follows a hybrid approach (deep learning and structure-based) to detect elements, filter partial sequences and finally classify intact sequences into superfamilies and, as very few tools do, into lineages. This tool takes advantage of multi-core and GPU architectures to decrease execution times. Using the rice genome, Inpactor2 showed a run time of 5 minutes (faster than other tools) and has the best accuracy and F1-Score of the tools tested here, also having the second best accuracy and specificity only surpassed by EDTA, but achieving 28% higher sensitivity. For large genomes, Inpactor2 is up to seven times faster than other available bioinformatics tools.

Retrotransposons: How the continuous evolutionary front shapes plant genomes for response to heat stress

Long terminal repeat retrotransposons (LTR retrotransposons) are the most abundant group of mobile genetic elements in eukaryotic genomes and are essential in organizing genomic architecture and phenotypic variations. The diverse families of retrotransposons are related to retroviruses. As retrotransposable elements are dispersed and ubiquitous, their "copy-out and paste-in" life cycle of replicative transposition leads to new genome insertions without the excision of the original element. The overall structure of retrotransposons and the domains responsible for the various phases of their replication is highly conserved in all eukaryotes. The two major superfamilies of LTR retrotransposons, Ty1/Copia and Ty3/Gypsy, are distinguished and dispersed across the chromosomes of higher plants. Members of these superfamilies can increase in copy number and are often activated by various biotic and abiotic stresses due to retrotransposition bursts. LTR retrotransposons are important drivers of species diversity and exhibit great variety in structure, size, and mechanisms of transposition, making them important putative actors in genome evolution. Additionally, LTR retrotransposons influence the gene expression patterns of adjacent genes by modulating potential small interfering RNA (siRNA) and RNA-directed DNA methylation (RdDM) pathways. Furthermore, comparative and evolutionary analysis of the most important crop genome sequences and advanced technologies have elucidated the epigenetics and structural and functional modifications driven by LTR retrotransposon during speciation. However, mechanistic insights into LTR retrotransposons remain obscure in plant development due to a lack of advancement in high throughput technologies. In this review, we focus on the key role of LTR retrotransposons response in plants during heat stress, the role of centromeric LTR retrotransposons, and the role of LTR retrotransposon markers in genome expression and evolution.

A New Hope: A Hermaphroditic Nematode Enables Analysis of a Recent Whole Genome Duplication Event

Whole genome duplication (WGD) is often considered a major driver of evolution that leads to phenotypic novelties. However, the importance of WGD for evolution is still controversial because most documented WGD events occurred anciently and few experimental systems amenable to genetic analysis are available. Here, we report a recent WGD event in the hermaphroditic nematode Allodiplogaster sudhausi and present a comparison with a gonochoristic (male/female) sister species that did not undergo WGD. Self-fertilizing reproduction of A. sudhausi makes it amenable to functional analysis and an ideal system to study WGD events. We document WGD in A. sudhausi through karyotype analysis and whole genome sequencing, the latter of which allowed us to 1) identify functional bias in retention of protein domains and metabolic pathways, 2) show most duplicate genes are under evolutionary constraint, 3) show a link between sequence and expression divergence, and 4) characterize differentially expressed duplicates. We additionally show WGD is associated with increased body size and an abundance of repeat elements (36% of the genome), including a recent expansion of the DNA-hAT/Ac transposon family. Finally, we demonstrate the use of CRISPR/Cas9 to generate mutant knockouts, whereby two WGD-derived duplicate genes display functional redundancy in that they both need to be knocked out to generate a phenotype. Together, we present a novel experimental system that is convenient for examining and characterizing WGD-derived genes both computationally and functionally.

The Vinca minor genome highlights conserved evolutionary traits in monoterpene indole alkaloid synthesis

Vinca minor, also known as the lesser periwinkle, is a well-known species from the Apocynaceae, native to central and southern Europe. This plant synthesizes monoterpene indole alkaloids, which are a class of specialized metabolites displaying a wide range of bioactive- and pharmacologically important properties. Within the almost 50 monoterpene indole alkaloids it produces, V. minor mainly accumulates vincamine, which is commercially used as a nootropic. Using a combination of Oxford Nanopore Technologies long read- and Illumina short-read sequencing, a 679,098 Mb V. minor genome was assembled into 296 scaffolds with an N50 scaffold length of 6 Mb, and encoding 29,624 genes. These genes were functionally annotated and used in a comparative genomic analysis to establish gene families and to investigate gene family expansion and contraction across the phylogenetic tree. Furthermore, homology-based monoterpene indole alkaloid gene predictions together with a metabolic analysis across 4 different V. minor tissue types guided the identification of candidate monoterpene indole alkaloid genes. These candidates were finally used to identify monoterpene indole alkaloid gene clusters, which combined with synteny analysis allowed for the discovery of a functionally validated vincadifformine-16-hydroxylase, reinforcing the potential of this dataset for monoterpene indole alkaloids gene discovery. It is expected that access to these resources will facilitate the elucidation of unknown monoterpene indole alkaloid biosynthetic routes with the potential of transferring these pathways to heterologous expression systems for large-scale monoterpene indole alkaloid production.

Evolution of lmiRNAs and their targets from MITEs for rice adaptation

Twenty-four nucleotide long microRNAs (lmiRNAs) direct DNA methylation at target genes and regulate their transcription. The evolutionary origin of lmiRNAs and the range of lmiRNA-mediated regulation remain obscure. Here, we reannotated lmiRNAs and their targets in rice by applying stringent criteria. We found that the majority of lmiRNAs are derived from Miniature Inverted-repeat Transposable Elements (MITEs) and most sites targeted by MITE-derived lmiRNAs reside within MITEs, suggesting co-evolution of lmiRNAs and their targets through MITE amplification. lmiRNAs undergo dynamically changes under stress conditions and the genes targeted by lmiRNAs show an enrichment for stress-responsive genes, suggesting that lmiRNAs are widely involved in plant responses to stresses. We constructed the evolutionary histories of lmiRNAs and their targets. Nearly half of lmiRNAs emerged before or when the AA genome was diverged, while the emergence of lmiRNA targets coincided with or followed the emergence of lmiRNAs. Furthermore, we found that the sequences of a lmiRNA target site underwent variations, coincident with the divergence of rice accessions and the distribution of rice accessions in different geographical locations and climatic conditions. Our findings highlight MITEs as an important origin of lmiRNAs and suggest that the evolution of lmiRNA-target regulatory modules may contribute to rice adaptation to environmental changes.

Morphological, physiological, and transcriptional responses of the freshwater diatom Fragilaria crotonensis to elevated pH conditions

Harmful algal blooms (HABs) caused by the toxin-producing cyanobacteria Microcystis spp., can increase water column pH. While the effect(s) of these basified conditions on the bloom formers are a high research priority, how these pH shifts affect other biota remains understudied. Recently, it was shown these high pH levels decrease growth and Si deposition rates in the freshwater diatom Fragilaria crotonensis and natural Lake Erie (Canada-US) diatom populations. However, the physiological mechanisms and transcriptional responses of diatoms associated with these observations remain to be documented. Here, we examined F. crotonensis with a set of morphological, physiological, and transcriptomic tools to identify cellular responses to high pH. We suggest 2 potential mechanisms that may contribute to morphological and physiological pH effects observed in F. crotonensis. Moreover, we identified a significant upregulation of mobile genetic elements in the F. crotonensis genome which appear to be an extreme transcriptional response to this abiotic stress to enhance cellular evolution rates-a process we have termed "genomic roulette." We discuss the ecological and biogeochemical effects high pH conditions impose on fresh waters and suggest a means by which freshwater diatoms such as F. crotonensis may evade high pH stress to survive in a "basified" future.

De novo assembly and characterization of the draft genome of the cashew (Anacardium occidentale L.)

Cashew is the second most important tree nut crop in the global market. Cashew is a diploid and heterozygous species closely related to the mango and pistachio. Its improvement by conventional breeding is slow due to the long juvenile phase. Despite the economic importance, very little genomics/transcriptomics information is available for cashew. In this study, the Oxford nanopore reads and Illumina reads were used for de novo assembly of the cashew genome. The hybrid assembly yielded a 356.6 Mb genome corresponding to 85% of the estimated genome size (419 Mb). The BUSCO analysis showed 91.8% of genome completeness. Transcriptome mapping showed 92.75% transcripts aligned with the assembled genome. Gene predictions resulted in the identification of 31,263 genes coding for a total of 35,000 gene isoforms. About 46% (165 Mb) of the cashew genome comprised of repetitive sequences. Phylogenetic analyses of the cashew with nine species showed that it was closely related to Mangifera indica. Analysis of cashew genome revealed 3104 putative R-genes. The first draft assembly of the genome, transcriptome and R gene information generated in this study would be the foundation for understanding the molecular basis of economic traits and genomics-assisted breeding in cashew.

Effect of DNA methylation, modified by 5-azaC, on ecophysiological responses of a clonal plant to changing climate

Epigenetic regulation of gene expression is expected to be an important mechanism behind phenotypic plasticity. Whether epigenetic regulation affects species ecophysiological adaptations to changing climate remains largely unexplored. We compared ecophysiological traits between individuals treated with 5-azaC, assumed to lead to DNA demethylation, with control individuals of a clonal grass originating from and grown under different climates, simulating different directions and magnitudes of climate change. We linked the ecophysiological data to proxies of fitness. Main effects of plant origin and cultivating conditions predicted variation in plant traits, but 5-azaC did not. Effects of 5-azaC interacted with conditions of cultivation and plant origin. The direction of the 5-azaC effects suggests that DNA methylation does not reflect species long-term adaptations to climate of origin and species likely epigenetically adjusted to the conditions experienced during experiment set-up. Ecophysiology translated to proxies of fitness, but the intensity and direction of the relationships were context dependent and affected by 5-azaC. The study suggests that effects of DNA methylation depend on conditions of plant origin and current climate. Direction of 5-azaC effects suggests limited role of epigenetic modifications in long-term adaptation of plants. It rather facilitates fast adaptations to temporal fluctuations of the environment.

The evolution and function of transposons in epigenetic regulation in response to the environment

Transposable elements (TEs) make up a major proportion of plant genomes. Despite their prevalence genome-wide, TEs are often tossed aside as "junk DNA" since they rarely cause phenotypes, and epigenetic mechanisms silence TEs to prevent them from causing deleterious mutations through movement. While this bleak picture of TEs in genomes is true on average, a growing number of examples across many plant species point to TEs as drivers of phenotypic diversity and novel stress responses. Examples of TE-influenced phenotypes illustrate the many ways that novel transposition events can alter local gene expression and how this relates to potential variation in plant responses to environmental stress. Since TE families and insertions at the locus level lack evolutionary conservation, advancements in the field will require TE experts across diverse species to identify and utilize TE variation in their own systems as a means of crop improvement.

Experimentally heat-induced transposition increases drought tolerance in Arabidopsis thaliana

Eukaryotic genomes contain a vast diversity of transposable elements (TEs). Formerly often described as selfish and parasitic DNA sequences, TEs are now recognised as a source of genetic diversity and powerful drivers of evolution. However, because their mobility is tightly controlled by the host, studies experimentally assessing how fast TEs may mediate the emergence of adaptive traits are scarce. We exposed Arabidopsis thaliana high-copy TE lines (hcLines) with up to c. eight-fold increased copy numbers of the heat-responsive ONSEN TE to drought as a straightforward and ecologically highly relevant selection pressure. We provide evidence for increased drought tolerance in five out of the 23 tested hcLines and further pinpoint one of the causative mutations to an exonic insertion of ONSEN in the ribose-5-phosphate-isomerase 2 gene. The resulting loss-of-function mutation caused a decreased rate of photosynthesis, plant size and water consumption. Overall, we show that the heat-induced transposition of a low-copy TE increases phenotypic diversity and leads to the emergence of drought-tolerant individuals in A. thaliana. This is one of the rare empirical examples substantiating the adaptive potential of mobilised stress-responsive TEs in eukaryotes. Our work demonstrates the potential of TE-mediated loss-of-function mutations in stress adaptation.

The ecology of palm genomes: repeat-associated genome size expansion is constrained by aridity

Genome size varies 2400-fold across plants, influencing their evolution through changes in cell size and cell division rates which impact plants' environmental stress tolerance. Repetitive element expansion explains much genome size diversity, and the processes structuring repeat 'communities' are analogous to those structuring ecological communities. However, which environmental stressors influence repeat community dynamics has not yet been examined from an ecological perspective. We measured genome size and leveraged climatic data for 91% of genera within the ecologically diverse palm family (Arecaceae). We then generated genomic repeat profiles for 141 palm species, and analysed repeats using phylogenetically informed linear models to explore relationships between repeat dynamics and environmental factors. We show that palm genome size and repeat 'community' composition are best explained by aridity. Specifically, Ty3-gypsy and TIR elements were more abundant in palm species from wetter environments, which generally had larger genomes, suggesting amplification. By contrast, Ty1-copia and LINE elements were more abundant in drier environments. Our results suggest that water stress inhibits repeat expansion through selection on upper genome size limits. However, elements that may associate with stress-response genes (e.g. Ty1-copia) have amplified in arid-adapted palm species. Overall, we provide novel evidence of climate influencing the assembly of repeat 'communities'.

The relationship between transposable elements and ecological niches in the Greater Cape Floristic Region: A study on the genus Pteronia (Asteraceae)

Non-coding repetitive DNA (repeatome) is an active part of the nuclear genome, involved in its structure, evolution and function. It is dominated by transposable elements (TEs) and satellite DNA and is prone to the most rapid changes over time. The TEs activity presumably causes the global genome reorganization and may play an adaptive or regulatory role in response to environmental challenges. This assumption is applied here for the first time to plants from the Cape Floristic hotspot to determine whether changes in repetitive DNA are related to responses to a harsh, but extremely species-rich environment. The genus Pteronia (Asteraceae) serves as a suitable model group because it shows considerable variation in genome size at the diploid level and has high and nearly equal levels of endemism in the two main Cape biomes, Fynbos and Succulent Karoo. First, we constructed a phylogeny based on multiple low-copy genes that served as a phylogenetic framework for detecting quantitative and qualitative changes in the repeatome. Second, we performed a comparative analysis of the environments of two groups of Pteronia differing in their TEs bursts. Our results suggest that the environmental transition from the Succulent Karoo to the Fynbos is accompanied by TEs burst, which is likely also driving phylogenetic divergence. We thus hypothesize that analysis of rapidly evolving repeatome could serve as an important proxy for determining the molecular basis of lineage divergence in rapidly radiating groups.

Transposable Elements in the Revealing of Polymorphism-Based Differences in the Seeds of Flax Varieties Grown in Remediated Chernobyl Area

The nuclear reactor accident in Chernobyl, Ukraine, resulted in effects both locally and farther away. Most of the contaminated areas were the agricultural fields and forests. Experimental fields were established near Chernobyl—radioactively contaminated fields localized 5 km from Chernobyl Nuclear Power Plant as well as the remediated soil that is localized directly in the Chernobyl town. Two flax varieties growing under chronic exposition to ionizing radiation were used for this study—the local Ukrainian variety Kyivskyi and a commercial variety Bethune. The screening of the length polymorphism generated by transposable elements insertions were performed. All known types of common flax transposon, retrotransposons and iPBS approach were used. In the iPBS multiplex analyze, for the Kyivskyi variety, a unique addition was found in the seeds from the radioactive contaminated field and for the Bethune variety, a total of five amplicon additions were obtained and one deletion. For the TRIM Cassandra fingerprints, two amplicon additions were generated in the seeds from radioactive contaminated fields for the Bethune variety. In summary, the obtained data represent the genetic diversity between control and irradiated subgroups of flax seeds from Chernobyl area and the presence of activated transposable elements due to the irradiation stress.

Ovule siRNAs methylate protein-coding genes in trans

Twenty-four-nucleotide (nt) small interfering RNAs (siRNAs) maintain asymmetric DNA methylation at thousands of euchromatic transposable elements in plant genomes in a process called RNA-directed DNA methylation (RdDM). RdDM is dispensable for growth and development in Arabidopsis thaliana, but is required for reproduction in other plants, such as Brassica rapa. The 24-nt siRNAs are abundant in maternal reproductive tissue, due largely to overwhelming expression from a few loci in the ovule and developing seed coat, termed siren loci. A recent study showed that 24-nt siRNAs produced in the anther tapetal tissue can methylate male meiocyte genes in trans. Here we show that in B. rapa, a similar process takes place in female tissue. siRNAs are produced from gene fragments embedded in some siren loci, and these siRNAs can trigger methylation in trans at related protein-coding genes. This trans-methylation is associated with silencing of some target genes and may be responsible for seed abortion in RdDM mutants. Furthermore, we demonstrate that a consensus sequence in at least two families of DNA transposons is associated with abundant siren expression, most likely through recruitment of CLASSY3, a putative chromatin remodeler. This research describes a mechanism whereby RdDM influences gene expression and sheds light on the role of RdDM during plant reproduction.

Genome-wide association study of traits in sacred lotus uncovers MITE-associated variants underlying stamen petaloid and petal number variations

Understanding the genetic variants responsible for floral trait diversity is important for the molecular breeding of ornamental flowers. Widely used in water gardening for thousands of years, the sacred lotus exhibits a wide range of diversity in floral organs. Nevertheless, the genetic variations underlying various morphological characteristics in lotus remain largely unclear. Here, we performed a genome-wide association study of sacred lotus for 12 well-recorded ornamental traits. Given a moderate linkage disequilibrium level of 32.9 kb, we successfully identified 149 candidate genes responsible for seven flower traits and plant size variations, including many pleiotropic genes affecting multiple floral-organ-related traits, such as NnKUP2. Notably, we found a 2.75-kb presence-and-absence genomic fragment significantly associated with stamen petaloid and petal number variations, which was further confirmed by re-examining another independent population dataset with petal number records. Intriguingly, this fragment carries MITE transposons bound by siRNAs and is related to the expression differentiation of a nearby candidate gene between few-petalled and double-petalled lotuses. Overall, these genetic variations and candidate genes responsible for diverse lotus traits revealed by our GWAS highlight the role of transposon variations, particularly MITEs, in shaping floral trait diversity.

Transposable elements

Transposable elements are known by many names, including 'transposons', 'interspersed repeats', 'selfish genetic elements', 'jumping genes', and 'parasitic DNA', but here we will refer to them simply as transposable elements. Many biologists will have heard of transposable elements and their ability to transpose (change position) within the genome. But fewer may be aware of their varied influences on host biology, including contributions to the evolution of diverse host traits such as internal gestation, memory, colouration, and adaptive immunity. Transposable elements are a near ubiquitous feature of eukaryotic genomes, and they often comprise a substantial proportion of total genomic content. Consequently, transposable element genes are considered among the most abundant coding sequences in nature. Recent advances in genome sequencing have ushered in a golden age for transposable-element research, providing opportunities to greatly improve our understanding of the effects of transposable elements on host evolution and disease. However, our ability to detect and analyse transposable elements still faces significant challenges, impairing efforts to decipher their evolution, characterise their diversity, and elucidate their myriad host influences. Below, we summarise key aspects of transposable element biology in eukaryotes and discuss major outstanding research questions.

Loss of pod strings in common bean is associated with gene duplication, retrotransposon insertion and overexpression of PvIND

Fruit development has been central in the evolution and domestication of flowering plants. In common bean (Phaseolus vulgaris), the principal global grain legume staple, two main production categories are distinguished by fibre deposition in pods: dry beans, with fibrous, stringy pods; and stringless snap/green beans, with reduced fibre deposition, which frequently revert to the ancestral stringy state. Here, we identify genetic and developmental patterns associated with pod fibre deposition. Transcriptional, anatomical, epigenetic and genetic regulation of pod strings were explored through RNA-seq, RT-qPCR, fluorescence microscopy, bisulfite sequencing and whole-genome sequencing. Overexpression of the INDEHISCENT ('PvIND') orthologue was observed in stringless types compared with isogenic stringy lines, associated with overspecification of weak dehiscence-zone cells throughout the pod vascular sheath. No differences in DNA methylation were correlated with this phenotype. Nonstringy varieties showed a tandemly direct duplicated PvIND and a Ty1-copia retrotransposon inserted between the two repeats. These sequence features are lost during pod reversion and are predictive of pod phenotype in diverse materials, supporting their role in PvIND overexpression and reversible string phenotype. Our results give insight into reversible gain-of-function mutations and possible genetic solutions to the reversion problem, of considerable economic value for green bean production.

Transposable element polymorphisms improve prediction of complex agronomic traits in rice

KEY MESSAGE

Transposon insertion polymorphisms can improve prediction of complex agronomic traits in rice compared to using SNPs only, especially when accessions to be predicted are less related to the training set. Transposon insertion polymorphisms (TIPs) are significant sources of genetic variation. Previous work has shown that TIPs can improve detection of causative loci on agronomic traits in rice. Here, we quantify the fraction of variance explained by single nucleotide polymorphisms (SNPs) compared to TIPs, and we explore whether TIPs can improve prediction of traits when compared to using only SNPs. We used eleven traits of agronomic relevance from by five different rice population groups (Aus, Indica, Aromatic, Japonica, and Admixed), 738 accessions in total. We assess prediction by applying data split validation in two scenarios. In the within-population scenario, we predicted performance of improved Indica varieties using the rest of Indica accessions. In the across population scenario, we predicted all Aromatic and Admixed accessions using the rest of populations. In each scenario, Bayes C and a Bayesian reproducible kernel Hilbert space regression were compared. We find that TIPs can explain an important fraction of total genetic variance and that they also improve genomic prediction. In the across population prediction scenario, TIPs outperformed SNPs in nine out of the eleven traits analyzed. In some traits like leaf senescence or grain width, using TIPs increased predictive correlation by 30-50%. Our results evidence, for the first time, that TIPs genotyping can improve prediction on complex agronomic traits in rice, especially when accessions to be predicted are less related to training accessions.

Dynamic genome evolution in a model fern

The large size and complexity of most fern genomes have hampered efforts to elucidate fundamental aspects of fern biology and land plant evolution through genome-enabled research. Here we present a chromosomal genome assembly and associated methylome, transcriptome and metabolome analyses for the model fern species Ceratopteris richardii. The assembly reveals a history of remarkably dynamic genome evolution including rapid changes in genome content and structure following the most recent whole-genome duplication approximately 60 million years ago. These changes include massive gene loss, rampant tandem duplications and multiple horizontal gene transfers from bacteria, contributing to the diversification of defence-related gene families. The insertion of transposable elements into introns has led to the large size of the Ceratopteris genome and to exceptionally long genes relative to other plants. Gene family analyses indicate that genes directing seed development were co-opted from those controlling the development of fern sporangia, providing insights into seed plant evolution. Our findings and annotated genome assembly extend the utility of Ceratopteris as a model for investigating and teaching plant biology.

Ruminant-specific retrotransposons shape regulatory evolution of bovine immunity

Cattle are an important livestock species, and mapping the genomic architecture of agriculturally relevant traits such as disease susceptibility is a major challenge in the bovine research community. Lineage-specific transposable elements (TEs) are increasingly recognized to contribute to gene regulatory evolution and variation, but this possibility has been largely unexplored in ruminant genomes. We conducted epigenomic profiling of the type II interferon (IFN) response in bovine cells and found thousands of ruminant-specific TEs including MER41_BT and Bov-A2 elements predicted to act as IFN-inducible enhancer elements. CRISPR knockout experiments in bovine cells established that critical immune factors including IFNAR2 and IL2RB are transcriptionally regulated by TE-derived enhancers. Finally, population genomic analysis of 38 individuals revealed that a subset of polymorphic TE insertions may function as enhancers in modern cattle. Our study reveals that lineage-specific TEs have shaped the evolution of ruminant IFN responses and potentially continue to contribute to immune gene regulatory differences across modern breeds and individuals. Together with previous work in human cells, our findings demonstrate that lineage-specific TEs have been independently co-opted to regulate IFN-inducible gene expression in multiple species, supporting TE co-option as a recurrent mechanism driving the evolution of IFN-inducible transcriptional networks.

Environmental stress and transposons in plants

Transposons were once thought to be junk repetitive DNA in the genome. However, their importance gradually became apparent as it became clear that they regulate gene expression, which is essential for organisms to survive, and that they are important factors in the driving force of evolution. Since there are multiple transposons in the genomes of all organisms, transposons have likely been activated and increased in copy number throughout their long history. This review focuses on environmental stress as a factor in transposon activation, paying particular attention to transposons in plants that are activated by environmental stresses. It is now known that plants respond to environmental stress in various ways, and correspondingly, many transposons respond to stress. The relationship between environmental stress and transposons is reviewed, including the mechanisms of their activation and the effects of transposon activation on host plants.

Ongoing endeavors to detect mobilization of transposable elements

Transposable elements (TEs) are DNA sequences capable of mobilization from one location to another in the genome. Since the discovery of ‘Dissociation (Dc) locus’ by Barbara McClintock in maize (1), mounting evidence in the era of genomics indicates that a significant fraction of most eukaryotic genomes is composed of TE sequences, involving in various aspects of biological processes such as development, physiology, diseases and evolution. Although technical advances in genomics have discovered numerous functional impacts of TE across species, our understanding of TEs is still ongoing process due to challenges resulted from complexity and abundance of TEs in the genome. In this mini-review, we briefly summarize biology of TEs and their impacts on the host genome, emphasizing importance of understanding TE landscape in the genome. Then, we introduce recent endeavors especially in vivo retrotransposition assays and long read sequencing technology for identifying de novo insertions/TE polymorphism, which will broaden our knowledge of extraordinary relationship between genomic cohabitants and their host.

Chromosome-scale assembly with a phased sex-determining region resolves features of early Z and W chromosome differentiation in a wild octoploid strawberry

When sex chromosomes stop recombining, they start to accumulate differences. The sex-limited chromosome (Y or W) especially is expected to degenerate via the loss of nucleotide sequence and the accumulation of repetitive sequences. However, how early signs of degeneration can be detected in a new sex chromosome is still unclear. The sex-determining region of the octoploid strawberries is young, small, and dynamic. Using PacBio HiFi reads, we obtained a chromosome-scale assembly of a female (ZW) Fragaria chiloensis plant carrying the youngest and largest of the known sex-determining region on the W in strawberries. We fully characterized the previously incomplete sex-determining region, confirming its gene content, genomic location, and evolutionary history. Resolution of gaps in the previous characterization of the sex-determining region added 10 kb of sequence including a noncanonical long terminal repeat-retrotransposon; whereas the Z sequence revealed a Harbinger transposable element adjoining the sex-determining region insertion site. Limited genetic differentiation of the sex chromosomes coupled with structural variation may indicate an early stage of W degeneration. The sex chromosomes have a similar percentage of repeats but differ in their repeat distribution. Differences in the pattern of repeats (transposable element polymorphism) apparently precede sex chromosome differentiation, thus potentially contributing to recombination cessation as opposed to being a consequence of it.

Automatic curation of LTR retrotransposon libraries from plant genomes through machine learning

Transposable elements are mobile sequences that can move and insert themselves into chromosomes, activating under internal or external stimuli, giving the organism the ability to adapt to the environment. Annotating transposable elements in genomic data is currently considered a crucial task to understand key aspects of organisms such as phenotype variability, species evolution, and genome size, among others. Because of the way they replicate, LTR retrotransposons are the most common transposable elements in plants, accounting in some cases for up to 80% of all DNA information. To annotate these elements, a reference library is usually created, a curation process is performed, eliminating TE fragments and false positives and then annotated in the genome using the homology method. However, the curation process can take weeks, requires extensive manual work and the execution of multiple time-consuming bioinformatics software. Here, we propose a machine learning-based approach to perform this process automatically on plant genomes, obtaining up to 91.18% F1-score. This approach was tested with four plant species, obtaining up to 93.6% F1-score (Oryza granulata) in only 22.61 s, where bioinformatics methods took approximately 6 h. This acceleration demonstrates that the ML-based approach is efficient and could be used in massive sequencing projects.

Identification and Characterization of PTE-2, a Stowaway-like MITE Activated in Transgenic Chinese Cabbage Lines

Transposable elements (TEs) are DNA fragments that can be replicated or transposed within a genome. TEs make up a high proportion of the plant genome and contribute to genetic diversity and evolution, affecting genome structure or gene activity. Miniature inverted-repeat transposable elements (MITEs) are short, non-autonomous class II DNA transposable elements. MITEs have specific sequences, target site duplications(TSDs), and terminal inverted repeats(TIRs), which are characteristics of the classification of MITE families. In this study, a Stowaway-like MITE, PTE-2, was activated in transgenic Chinese cabbage lines. PTE-2 was revealed by in silico analysis as the putative activated element in transgenic Chinese cabbage lines. To verify the in silico analysis data, MITE insertion polymorphism (MIP) PCR was conducted and PTE-2 was confirmed to be activated in transgenic Chinese cabbage lines. The activation tendency of the copy elements of PTE-2 at different loci was also analyzed and only one more element was activated in the transgenic Chinese cabbage lines. Analyzing the sequence of MIP PCR products, the TSD sequence and TIR motif of PTE-2 were identified and matched to the characteristics of the Stowaway-like MITE family. In addition, the flanking region of PTE-2 was modified when PTE-2 was activated.

Ongoing endeavors to detect mobilization of transposable elements

Transposable elements (TEs) are DNA sequences capable of mobilization from one location to another in the genome. Since the discovery of 'Dissociation (Dc) locus' by Barbara McClintock in maize (1), mounting evidence in the era of genomics indicates that a significant fraction of most eukaryotic genomes is composed of TE sequences, involving in various aspects of biological processes such as development, physiology, diseases and evolution. Although technical advances in genomics have discovered numerous functional impacts of TE across species, our understanding of TEs is still ongoing process due to challenges resulted from complexity and abundance of TEs in the genome. In this mini-review, we briefly summarize biology of TEs and their impacts on the host genome, emphasizing importance of understanding TE landscape in the genome. Then, we introduce recent endeavors especially in vivo retrotransposition assays and long read sequencing technology for identifying de novo insertions/TE polymorphism, which will broaden our knowledge of extraordinary relationship between genomic cohabitants and their host. [BMB Reports 2022; 55(7): 305-315].

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.

Transcriptome Profiling of Transposon-Derived Long Non-coding RNAs Response to Hormone in Strawberry Fruit Development

Strawberry is an economically grown horticulture crop required for fruit consumption. The ripening of its fruit is a complex biological process regulated by various hormones. Abscisic acid (ABA) is a critical phytohormone involved in fruit ripening. However, little is known about the long non-coding RNAs (LncRNAs), especially transposon-derived LncRNA (TE-lncRNA), response to hormones during fruit ripening in octoploid strawberry. In the study, the transcriptome data of developing strawberry fruits treated with ABA and its inhibitor Nordihydroguaiaretic acid (NGDA) were analyzed to identify responsive LncRNAs and coding genes. A total of 14,552 LncRNAs were identified, including 8,617 transposon-derived LncRNAs (TE-LncRNAs), 412 LncRNAs (282 TE-LncRNAs), and 382 ABA-sensitive LncRNAs (231 TE-LncRNAs). Additionally, a weighted co-expression network analysis constructed 27 modules containing coding RNAs and LncRNAs. Seven modules, including "MEdarkorange" and "MElightyellow" were significantly correlated with ABA/NDGA treatments, resulting in 247 hub genes, including 21 transcription factors and 22 LncRNAs (15 TE-LncRNAs). Gene ontology enrichment analysis further revealed that ABA/NDGA-responsive modules, including LncRNAs, were associated with various metabolic pathways involved in strawberry fruit development and ripening, including lipid metabolism, organic acid metabolism, and phenylpropanoid metabolism. The current study identifies many high-confidence LncRNAs in strawberry, with a percentage of them being ABA pathway-specific and 22 hub-responsive LncRNAs, providing new insight into strawberry or other Rosaceae crop fruit ripening.

Genomes and demographic histories of the endangered Bretschneidera sinensis (Akaniaceae)

BACKGROUND

Bretschneidera sinensis is an endangered relic tree species in the Akaniaceae family and is sporadically distributed in eastern Asia. As opposed to its current narrow and rare distribution, the fossil pollen of B. sinensis has been found to be frequent and widespread in the Northern Hemisphere during the Late Miocene. B. sinensis is also a typical mycorrhizal plant, and its annual seedlings exhibit high mortality rates in absence of mycorrhizal development. The chromosome-level high-quality genome of B. sinensis will help us to more deeply understand the survival and demographic histories of this relic species.

RESULTS

A total of 25.39 Gb HiFi reads and 109.17 Gb Hi-C reads were used to construct the chromosome-level genome of B. sinensis, which is 1.21 Gb in length with the contig N50 of 64.13 Mb and chromosome N50 of 146.54 Mb. The identified transposable elements account for 55.21% of the genome. A total of 45,839 protein-coding genes were predicted in B. sinensis. A lineage-specific whole-genome duplication was detected, and 7,283 lineage-specific expanded gene families with functions related to the specialized endotrophic mycorrhizal adaptation were identified. The historical effective population size (Ne) of B. sinensis was found to oscillate greatly in response to Quaternary climatic changes. The Ne of B. sinensis has decreased rapidly in the recent past, making its extant Ne extremely lower. Our additional evolutionary genomic analyses suggested that the developed mycorrhizal adaption might have been repeatedly disrupted by environmental changes caused by Quaternary climatic oscillations. The environmental changes and an already decreased population size during the Holocene may have led to the current rarity of B. sinensis.

CONCLUSION

This is a detailed report of the genome sequences for the family Akaniaceae distributed in evergreen forests in eastern Asia. Such a high-quality genomic resource may provide critical clues for comparative genomics studies of this family in the future.

Gypsy retrotransposon-derived maize lncRNA GARR2 modulates gibberellin response

Long non-coding RNAs (lncRNAs) mediate diverse biological events mainly through the modulation of transcriptional hierarchy. The phytohormone gibberellin (GA) is essential for various aspects of plant growth and development. However, the roles of lncRNAs in the regulation of the GA response remain largely unknown. Through sequencing multiple strand-specific and ribosomal-depleted RNA libraries, we delineated the landscape of lncRNAs in maize (Zea mays). Out of identified lncRNAs, 445 GIBBERELLIN-RESPONSIVE lncRNAs (GARRs) were differentially expressed upon GA application. By the intersection of GARRs from normal-height and dwarf plants from an advanced backcross population, four shared GARRs (GARR1 to GARR4) were identified. Out of these four shared GARRs, GARR2 was derived from a Gypsy LTR retrotransposon. GA-responsive element P-boxes were identified upstream of GARR2. GARR2-edited lines exhibited a GA-induced phenotype. Editing of GARR2 resulted in changes in the transcriptional abundance of GA pathway components and endogenous GA contents. Besides GA, GARR2 affected the primary auxin response. An RNA pull-down assay revealed the HECT ubiquitin-protein ligase family member ZmUPL1 as a potential interaction target of GARR2. GARR2 influenced the abundance of ZmUPL1 in the GA response. Our study uncovers lncRNA players involved in the modulation of the GA response and guides the development of plant height ideotype driven by knowledge of the phytohormone GA.

Dynamically expressed small RNAs, substantially driven by genomic structural variants, contribute to transcriptomic changes during tomato domestication

Tomato has undergone extensive selections during domestication. Recent progress has shown that genomic structural variants (SVs) have contributed to gene expression dynamics during tomato domestication, resulting in changes of important traits. Here, we performed comprehensive analyses of small RNAs (sRNAs) from nine representative tomato accessions. We demonstrate that SVs substantially contribute to the dynamic expression of the three major classes of plant sRNAs: microRNAs (miRNAs), phased secondary short interfering RNAs (phasiRNAs), and 24-nucleotide heterochromatic siRNAs (hc-siRNAs). Changes in the abundance of phasiRNAs and 24-nucleotide hc-siRNAs likely contribute to the alteration of mRNA gene expression in cis during tomato domestication, particularly for genes associated with biotic and abiotic stress tolerance. We also observe that miRNA expression dynamics are associated with imprecise processing, alternative miRNA-miRNA* selections, and SVs. SVs mainly affect the expression of less-conserved miRNAs that do not have established regulatory functions or low abundant members in highly expressed miRNA families. Our data highlight different selection pressures on miRNAs compared to phasiRNAs and 24-nucleotide hc-siRNAs. Our findings provide insights into plant sRNA evolution as well as SV-based gene regulation during crop domestication. Furthermore, our dataset provides a rich resource for mining the sRNA regulatory network in tomato.

Local Adaptation and the Evolution of Genome Architecture in Threespine Stickleback

Theory predicts that local adaptation should favor the evolution of a concentrated genetic architecture, where the alleles driving adaptive divergence are tightly clustered on chromosomes. Adaptation to marine versus freshwater environments in threespine stickleback has resulted in an architecture that seems consistent with this prediction: divergence among populations is mainly driven by a few genomic regions harboring multiple quantitative trait loci for environmentally adapted traits, as well as candidate genes with well-established phenotypic effects. One theory for the evolution of these "genomic islands" is that rearrangements remodel the genome to bring causal loci into tight proximity, but this has not been studied explicitly. We tested this theory using synteny analysis to identify micro- and macro-rearrangements in the stickleback genome and assess their potential involvement in the evolution of genomic islands. To identify rearrangements, we conducted a de novo assembly of the closely related tubesnout (Aulorhyncus flavidus) genome and compared this to the genomes of threespine stickleback and two other closely related species. We found that small rearrangements, within-chromosome duplications, and lineage-specific genes (LSGs) were enriched around genomic islands, and that all three chromosomes harboring large genomic islands have experienced macro-rearrangements. We also found that duplicates and micro-rearrangements are 9.9× and 2.9× more likely to involve genes differentially expressed between marine and freshwater genotypes. While not conclusive, these results are consistent with the explanation that strong divergent selection on candidate genes drove the recruitment of rearrangements to yield clusters of locally adaptive loci.

Construction and characterization of a de novo draft genome of garden cress (Lepidium sativum L.)

Garden cress (Lepidium sativum L.) is a Brassicaceae crop recognized as a healthy vegetable and a medicinal plant. Lepidium is one of the largest genera in Brassicaceae, yet, the genus has not been a focus of extensive genomic research. In the present work, garden cress genome was sequenced using the long read high-fidelity sequencing technology. A de novo, draft genome assembly that spans 336.5 Mb was produced, corresponding to 88.6% of the estimated genome size and representing 90% of the evolutionarily expected orthologous gene content. Protein coding gene content was structurally predicted and functionally annotated, resulting in the identification of 25,668 putative genes. A total of 599 candidate disease resistance genes were identified by predicting resistance gene domains in gene structures, and 37 genes were detected as orthologs of heavy metal associated protein coding genes. In addition, 4289 genes were assigned as "transcription factor coding." Six different machine learning algorithms were trained and tested for their performance in classifying miRNA coding genomic sequences. Logistic regression proved the best performing trained algorithm, thus utilized for pre-miRNA coding loci identification in the assembly. Repetitive DNA analysis involved the characterization of transposable element and microsatellite contents. L. sativum chloroplast genome was also assembled and functionally annotated. Data produced in the present work is expected to constitute a foundation for genomic research in garden cress and contribute to genomics-assisted crop improvement and genome evolution studies in the Brassicaceae family.

Oil plant genomes: current state of the science

Vegetable oils are an indispensable nutritional component of the human diet as well as important raw materials for a variety of industrial applications such as pharmaceuticals, cosmetics, oleochemicals, and biofuels. Oil plant genomes are highly diverse, and their genetic variation leads to a diversity in oil biosynthesis and accumulation along with agronomic traits. This review discusses plant oil biosynthetic pathways, current state of genome assembly, polyploidy and asymmetric evolution of genomes of oil plants and their wild relatives, and research progress of pan-genomics in oil plants. The availability of complete high-resolution genomes and pan-genomes has enabled the identification of structural variations in the genomes that are associated with the diversity of agronomic and environment fitness traits. These and future genomes also provide powerful tools to understand crop evolution and to harvest the rich natural variations to improve oil crops for enhanced productivity, oil quality, and adaptability to changing environments.

SunUp and Sunset genomes revealed impact of particle bombardment mediated transformation and domestication history in papaya

Transgenic papaya is widely publicized for controlling papaya ringspot virus. However, the impact of particle bombardment on the genome remains unknown. The transgenic SunUp and its progenitor Sunset genomes were assembled into 351.5 and 350.3 Mb in nine chromosomes, respectively. We identified a 1.64 Mb insertion containing three transgenic insertions in SunUp chromosome 5, consisting of 52 nuclear-plastid, 21 nuclear-mitochondrial and 1 nuclear genomic fragments. A 591.9 kb fragment in chromosome 5 was translocated into the 1.64 Mb insertion. We assembled a gapless 9.8 Mb hermaphrodite-specific region of the Y^h chromosome and its 6.0 Mb X counterpart. Resequencing 86 genomes revealed three distinct groups, validating their geographic origin and breeding history. We identified 147 selective sweeps and defined the essential role of zeta-carotene desaturase in carotenoid accumulation during domestication. Our findings elucidated the impact of particle bombardment and improved our understanding of sex chromosomes and domestication to expedite papaya improvement.

Comprehensive analysis of both long and short read transcriptomes of a clonal and a seed-propagated model species reveal the prerequisites for transcriptional activation of autonomous and non-autonomous transposons in plants

Background

Transposable element (TE) transcription is a precursor to its mobilisation in host genomes. However, the characteristics of expressed TE loci, the identification of self-competent transposon loci contributing to new insertions, and the genomic conditions permitting their mobilisation remain largely unknown.

Results

Using Vitis vinifera embryogenic callus, we explored the impact of biotic stressors on transposon transcription through the exposure of the callus to live cultures of an endemic grapevine yeast, Hanseniaspora uvarum. We found that only 1.7–2.5% of total annotated TE loci were transcribed, of which 5–10% of these were full-length, and the expressed TE loci exhibited a strong location bias towards expressed genes. These trends in transposon transcription were also observed in RNA-seq data from Arabidopsis thaliana wild-type plants but not in epigenetically compromised Arabidopsis ddm1 mutants. Moreover, differentially expressed TE loci in the grapevine tended to share expression patterns with co-localised differentially expressed genes. Utilising nanopore cDNA sequencing, we found a strong correlation between the inclusion of intronic TEs in gene transcripts and the presence of premature termination codons in these transcripts. Finally, we identified low levels of full-length transcripts deriving from structurally intact TE loci in the grapevine model.

Conclusion

Our observations in two disparate plant models representing clonally and seed propagated plant species reveal a closely connected transcriptional relationship between TEs and co-localised genes, particularly when epigenetic silencing is not compromised. We found that the stress treatment alone was insufficient to induce large-scale full-length transcription from structurally intact TE loci, a necessity for non-autonomous and autonomous mobilisation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-022-00271-5.

Exploring transposable element-based markers to identify allelic variations underlying agronomic traits in rice

Transposable elements (TEs) are a major force in the production of new alleles during domestication; nevertheless, their use in association studies has been limited because of their complexity. We have developed a TE genotyping pipeline (TEmarker) and applied it to whole-genome genome-wide association study (GWAS) data from 176 Oryza sativa subsp. japonica accessions to identify genetic elements associated with specific agronomic traits. TE markers recovered a large proportion (69%) of single-nucleotide polymorphism (SNP)-based GWAS peaks, and these TE peaks retained ca. 25% of the SNPs. The use of TEs in GWASs may reduce false positives associated with linkage disequilibrium (LD) among SNP markers. A genome scan revealed positive selection on TEs associated with agronomic traits. We found several cases of insertion and deletion variants that potentially resulted from the direct action of TEs, including an allele of LOC_Os11g08410 associated with plant height and panicle length traits. Together, these findings reveal the utility of TE markers for connecting genotype to phenotype and suggest a potential role for TEs in influencing phenotypic variations in rice that impact agronomic traits.

Patterns of selection in the evolution of a transposable element

Transposable elements are a major component of most eukaryotic genomes. Here, we present a new approach which allows us to study patterns of natural selection in the evolution of transposable elements over short time scales. The method uses the alignment of all elements with intact gag/pol genes of a transposable element family from a single genome. We predict that the ratio of nonsynonymous to synonymous variants in the alignment should decrease as a function of the frequency of the variants, because elements with nonsynonymous variants that reduce transposition will have fewer progeny. We apply our method to Sirevirus long-terminal repeat retrotransposons that are abundant in maize and other plant species and show that nonsynonymous to synonymous variants declines as variant frequency increases, indicating that negative selection is acting strongly on the Sirevirus genome. The asymptotic value of nonsynonymous to synonymous variants suggests that at least 85% of all nonsynonymous mutations in the transposable element reduce transposition. Crucially, these patterns in nonsynonymous to synonymous variants are only predicted to occur if the gene products from a particular transposable element insertion preferentially promote the transposition of the same insertion. Overall, by using large numbers of intact elements, this study sheds new light on the selective processes that act on transposable elements.

Broad-spectrum fungal resistance in sorghum is conferred through the complex regulation of an immune receptor gene embedded in a natural antisense transcript

Sorghum (Sorghum bicolor), the fifth most widely grown cereal crop globally, provides food security for millions of people. Anthracnose caused by the fungus Colletotrichum sublineola is a major disease of sorghum worldwide. We discovered a major fungal resistance locus in sorghum composed of the nucleotide-binding leucine-rich repeat receptor gene ANTHRACNOSE RESISTANCE GENE1 (ARG1) that is completely nested in an intron of a cis-natural antisense transcript (NAT) gene designated CARRIER OF ARG1 (CARG). Susceptible genotypes express CARG and two alternatively spliced ARG1 transcripts encoding truncated proteins lacking the leucine-rich repeat domains. In resistant genotypes, elevated expression of an intact allele of ARG1, attributed to the loss of CARG transcription and the presence of miniature inverted-repeat transposable element sequences, resulted in broad-spectrum resistance to fungal pathogens with distinct virulence strategies. Increased ARG1 expression in resistant genotypes is also associated with higher histone H3K4 and H3K36 methylation. In susceptible genotypes, lower ARG1 expression is associated with reduced H3K4 and H3K36 methylation and increased expression of NATs of CARG. The repressive chromatin state associated with H3K9me2 is low in CARG-expressing genotypes within the CARG exon and higher in genotypes with low CARG expression. Thus, ARG1 is regulated by multiple mechanisms and confers broad-spectrum, strong resistance to fungal pathogens.

TE Density: a tool to investigate the biology of transposable elements

BACKGROUND

Transposable elements (TEs) are powerful creators of genotypic and phenotypic diversity due to their inherent mutagenic capabilities and in this way they serve as a deep reservoir of sequences for genomic variation. As agents of genetic disruption, a TE's potential to impact phenotype is partially a factor of its location in the genome. Previous research has shown TEs' ability to impact the expression of neighboring genes, however our understanding of this trend is hampered by the exceptional amount of diversity in the TE world, and a lack of publicly available computational methods that quantify the presence of TEs relative to genes.

RESULTS

Here, we have developed a tool to more easily quantify TE presence relative to genes through the use of only a gene and TE annotation, yielding a new metric we call TE Density. Briefly defined as the proportion of TE-occupied base-pairs relative to a window-size of the genome. This new pipeline reports TE density for each gene in the genome, for each type descriptor of TE (order and superfamily), and for multiple positions and distances relative to the gene (upstream, intragenic, and downstream) over sliding, user-defined windows. In this way, we overcome previous limitations to the study of TE-gene relationships by focusing on all TE types present in the genome, utilizing flexible genomic distances for measurement, and reporting a TE presence metric for every gene in the genome.

CONCLUSIONS

Together, this new tool opens up new avenues for studying TE-gene relationships, genome architecture, comparative genomics, and the tremendous diversity present of the TE world. TE Density is open-source and freely available at: https://github.com/sjteresi/TE_Density .