The genomic ecosystem of transposable elements in maize

HiTE: a fast and accurate dynamic boundary adjustment approach for full-length transposable element detection and annotation

Recent advancements in genome assembly have greatly improved the prospects for comprehensive annotation of Transposable Elements (TEs). However, existing methods for TE annotation using genome assemblies suffer from limited accuracy and robustness, requiring extensive manual editing. In addition, the currently available gold-standard TE databases are not comprehensive, even for extensively studied species, highlighting the critical need for an automated TE detection method to supplement existing repositories. In this study, we introduce HiTE, a fast and accurate dynamic boundary adjustment approach designed to detect full-length TEs. The experimental results demonstrate that HiTE outperforms RepeatModeler2, the state-of-the-art tool, across various species. Furthermore, HiTE has identified numerous novel transposons with well-defined structures containing protein-coding domains, some of which are directly inserted within crucial genes, leading to direct alterations in gene expression. A Nextflow version of HiTE is also available, with enhanced parallelism, reproducibility, and portability.

Natural polymorphisms in ZMET2 encoding a DNA methyltransferase modulate the number of husk layers in maize

DNA methylation affects agronomic traits and the environmental adaptability of crops, but the natural polymorphisms in DNA methylation-related genes and their contributions to phenotypic variation in maize (Zea mays) remain elusive. Here, we show that a polymorphic 10-bp insertion/deletion variant in the 3'UTR of Zea methyltransferase2 (ZMET2) alters its transcript level and accounts for variation in the number of maize husk layers. ZMET2 encodes a chromomethylase and is required for maintaining genome-wide DNA methylation in the CHG sequence context. Disruption of ZMET2 increased the number of husk layers and resulted in thousands of differentially methylated regions, a proportion of which were also distinguishable in natural ZMET2 alleles. Population genetic analyses indicated that ZMET2 was a target of selection and might play a role in the spread of maize from tropical to temperate regions. Our results provide important insights into the natural variation of ZMET2 that confers both global and locus-specific effects on DNA methylation, which contribute to phenotypic diversity in maize.

Genome Assembly of the Dyeing Poison Frog Provides Insights into the Dynamics of Transposable Element and Genome-Size Evolution

Genome size varies greatly across the tree of life and transposable elements are an important contributor to this variation. Among vertebrates, amphibians display the greatest variation in genome size, making them ideal models to explore the causes and consequences of genome size variation. However, high-quality genome assemblies for amphibians have, until recently, been rare. Here, we generate a high-quality genome assembly for the dyeing poison frog, Dendrobates tinctorius. We compare this assembly to publicly available frog genomes and find evidence for both large-scale conserved synteny and widespread rearrangements between frog lineages. Comparing conserved orthologs annotated in these genomes revealed a strong correlation between genome size and gene size. To explore the cause of gene-size variation, we quantified the location of transposable elements relative to gene features and find that the accumulation of transposable elements in introns has played an important role in the evolution of gene size in D. tinctorius, while estimates of insertion times suggest that many insertion events are recent and species-specific. Finally, we carry out population-scale mobile-element sequencing and show that the diversity and abundance of transposable elements in poison frog genomes can complicate genotyping from repetitive element sequence anchors. Our results show that transposable elements have clearly played an important role in the evolution of large genome size in D. tinctorius. Future studies are needed to fully understand the dynamics of transposable element evolution and to optimize primer or bait design for cost-effective population-level genotyping in species with large, repetitive genomes.

ACMGA: a reference-free multiple-genome alignment pipeline for plant species

BACKGROUND

The short-read whole-genome sequencing (WGS) approach has been widely applied to investigate the genomic variation in the natural populations of many plant species. With the rapid advancements in long-read sequencing and genome assembly technologies, high-quality genome sequences are available for a group of varieties for many plant species. These genome sequences are expected to help researchers comprehensively investigate any type of genomic variants that are missed by the WGS technology. However, multiple genome alignment (MGA) tools designed by the human genome research community might be unsuitable for plant genomes.

RESULTS

To fill this gap, we developed the AnchorWave-Cactus Multiple Genome Alignment (ACMGA) pipeline, which improved the alignment of repeat elements and could identify long (> 50 bp) deletions or insertions (INDELs). We conducted MGA using ACMGA and Cactus for 8 Arabidopsis (Arabidopsis thaliana) and 26 Maize (Zea mays) de novo assembled genome sequences and compared them with the previously published short-read variant calling results. MGA identified more single nucleotide variants (SNVs) and long INDELs than did previously published WGS variant callings. Additionally, ACMGA detected significantly more SNVs and long INDELs in repetitive regions and the whole genome than did Cactus. Compared with the results of Cactus, the results of ACMGA were more similar to the previously published variants called using short-read. These two MGA pipelines identified numerous multi-allelic variants that were missed by the WGS variant calling pipeline.

CONCLUSIONS

Aligning de novo assembled genome sequences could identify more SNVs and INDELs than mapping short-read. ACMGA combines the advantages of AnchorWave and Cactus and offers a practical solution for plant MGA by integrating global alignment, a 2-piece-affine-gap cost strategy, and the progressive MGA algorithm.

Genome-wide expression quantitative trait locus analysis reveals silk-preferential gene regulatory network in maize

Silk of maize (Zea mays L.) contains diverse metabolites with complicated structures and functions, making it a great challenge to explore the mechanisms of metabolic regulation. Genome-wide identification of silk-preferential genes and investigation of their expression regulation provide an opportunity to reveal the regulatory networks of metabolism. Here, we applied the expression quantitative trait locus (eQTL) mapping on a maize natural population to explore the regulation of gene expression in unpollinated silk of maize. We obtained 3,985 silk-preferential genes that were specifically or preferentially expressed in silk using our population. Silk-preferential genes showed more obvious expression variations compared with broadly expressed genes that were ubiquitously expressed in most tissues. We found that trans-eQTL regulation played a more important role for silk-preferential genes compared to the broadly expressed genes. The relationship between 38 transcription factors and 85 target genes, including silk-preferential genes, were detected. Finally, we constructed a transcriptional regulatory network around the silk-preferential gene Bx10, which was proposed to be associated with response to abiotic stress and biotic stress. Taken together, this study deepened our understanding of transcriptome variation in maize silk and the expression regulation of silk-preferential genes, enhancing the investigation of regulatory networks on metabolic pathways.

Dynamics of accessible chromatin regions and subgenome dominance in octoploid strawberry

Subgenome dominance has been reported in diverse allopolyploid species, where genes from one subgenome are preferentially retained and are more highly expressed than those from other subgenome(s). However, the molecular mechanisms responsible for subgenome dominance remain poorly understood. Here, we develop genome-wide map of accessible chromatin regions (ACRs) in cultivated strawberry (2n = 8x = 56, with A, B, C, D subgenomes). Each ACR is identified as an MNase hypersensitive site (MHS). We discover that the dominant subgenome A contains a greater number of total MHSs and MHS per gene than the submissive B/C/D subgenomes. Subgenome A suffers fewer losses of MHS-related DNA sequences and fewer MHS fragmentations caused by insertions of transposable elements. We also discover that genes and MHSs related to stress response have been preferentially retained in subgenome A. We conclude that preservation of genes and their cognate ACRs, especially those related to stress responses, play a major role in the establishment of subgenome dominance in octoploid strawberry.

Transposition of HOPPLA in siRNA-deficient plants suggests a limited effect of the environment on retrotransposon mobility in Brachypodium distachyon

Long terminal repeat retrotransposons (LTR-RTs) are powerful mutagens regarded as a major source of genetic novelty and important drivers of evolution. Yet, the uncontrolled and potentially selfish proliferation of LTR-RTs can lead to deleterious mutations and genome instability, with large fitness costs for their host. While population genomics data suggest that an ongoing LTR-RT mobility is common in many species, the understanding of their dual role in evolution is limited. Here, we harness the genetic diversity of 320 sequenced natural accessions of the Mediterranean grass Brachypodium distachyon to characterize how genetic and environmental factors influence plant LTR-RT dynamics in the wild. When combining a coverage-based approach to estimate global LTR-RT copy number variations with mobilome-sequencing of nine accessions exposed to eight different stresses, we find little evidence for a major role of environmental factors in LTR-RT accumulations in B. distachyon natural accessions. Instead, we show that loss of RNA polymerase IV (Pol IV), which mediates RNA-directed DNA methylation in plants, results in high transcriptional and transpositional activities of RLC_BdisC024 (HOPPLA) LTR-RT family elements, and that these effects are not stress-specific. This work supports findings indicating an ongoing mobility in B. distachyon and reveals that host RNA-directed DNA methylation rather than environmental factors controls their mobility in this wild grass model.

New whole-genome alignment tools are needed for tapping into plant diversity

Genome alignment is one of the most foundational methods for genome sequence studies. With rapid advances in sequencing and assembly technologies, these newly assembled genomes present challenges for alignment tools to meet the increased complexity and scale. Plant genome alignment is technologically challenging because of frequent whole-genome duplications (WGDs) as well as chromosome rearrangements and fractionation, high nucleotide diversity, widespread structural variation, and high transposable element (TE) activity causing large proportions of repeat elements. We summarize classical pairwise and multiple genome alignment (MGA) methods, and highlight techniques that are widely used or are being developed by the plant research community. We also outline the remaining challenges for precise genome alignment and the interpretation of alignment results in plants.

The good, the bad and the ugly of transposable elements annotation tools

Transposable elements are repetitive and mobile DNA segments that can be found in virtually all organisms investigated to date. Their complex structure and variable nature are particularly challenging from the genomic annotation point of view. Many softwares have been developed to automate and facilitate TEs annotation at the genomic level, but they are highly heterogeneous regarding documentation, usability and methods. In this review, we revisited the existing software for TE genomic annotation, concentrating on the most often used ones, the methodologies they apply, and usability. Building on the state of the art of TE annotation software we propose best practices and highlight the strengths and weaknesses from the available solutions.

Helenus and Ajax, Two Groups of Non-Autonomous LTR Retrotransposons, Represent a New Type of Small RNA Gene-Derived Mobile Elements

Terminal repeat retrotransposons in miniature (TRIMs) are short non-autonomous long terminal repeat (LTR) retrotransposons found from various eukaryotes. Cassandra is a unique TRIM lineage which contains a 5S rRNA-derived sequence in its LTRs. Here, two new groups of TRIMs, designated Helenus and Ajax, are reported based on bioinformatics analysis and the usage of Repbase. Helenus is found from fungi, animals, and plants, and its LTRs contain a tRNA-like sequence. It includes two LTRs and between them, a primer-binding site (PBS) and polypurine tract (PPT) exist. Fungal and plant Helenus generate 5 bp target site duplications (TSDs) upon integration, while animal Helenus generates 4 bp TSDs. Ajax includes a 5S rRNA-derived sequence in its LTR and is found from two nemertean genomes. Ajax generates 5 bp TSDs upon integration. These results suggest that despite their unique promoters, Helenus and Ajax are TRIMs whose transposition is dependent on autonomous LTR retrotransposon. These TRIMs can originate through an insertion of SINE in an LTR of TRIM. The discovery of Helenus and Ajax suggests the presence of TRIMs with a promoter for RNA polymerase III derived from a small RNA gene, which is here collectively termed TRIMp3.

Integrase-associated niche differentiation of endogenous large DNA viruses in crustaceans

IMPORTANCE

Crustacean genomes harbor sequences originating from a family of large DNA viruses called nimaviruses, but it is unclear why they are present. We show that endogenous nimaviruses selectively insert into repetitive sequences within the host genome, and this insertion specificity was correlated with different types of integrases, which are DNA recombination enzymes encoded by the nimaviruses themselves. This suggests that endogenous nimaviruses have colonized various genomic niches through the acquisition of integrases with different insertion specificities. Our results point to a novel survival strategy of endogenous large DNA viruses colonizing the host genomes. These findings may clarify the evolution and spread of nimaviruses in crustaceans and lead to measures to control and prevent the spread of pathogenic nimaviruses in aquaculture settings.

Quantifying the role of genome size and repeat content in adaptive variation and the architecture of flowering time in Amaranthus tuberculatus

Genome size variation, largely driven by repeat content, is poorly understood within and among populations, limiting our understanding of its significance for adaptation. Here we characterize intraspecific variation in genome size and repeat content across 186 individuals of Amaranthus tuberculatus, a ubiquitous native weed that shows flowering time adaptation to climate across its range and in response to agriculture. Sequence-based genome size estimates vary by up to 20% across individuals, consistent with the considerable variability in the abundance of transposable elements, unknown repeats, and rDNAs across individuals. The additive effect of this variation has important phenotypic consequences-individuals with more repeats, and thus larger genomes, show slower flowering times and growth rates. However, compared to newly-characterized gene copy number and polygenic nucleotide changes underlying variation in flowering time, we show that genome size is a marginal contributor. Differences in flowering time are reflected by genome size variation across sexes and marginally, habitats, while polygenic variation and a gene copy number variant within the ATP synthesis pathway show consistently stronger environmental clines than genome size. Repeat content nonetheless shows non-neutral distributions across the genome, and across latitudinal and environmental gradients, demonstrating the numerous governing processes that in turn influence quantitative genetic variation for phenotypes key to plant adaptation.

Repeat-based phylogenomics shed light on unclear relationships in the monocentric genus Juncus L. (Juncaceae)

The repetitive fraction (repeatome) of eukaryotic genomes is diverse and usually fast evolving, being an important tool for clarify plant systematics. The genus Juncus L. comprises 332 species, karyotypically recognized by having holocentric chromosomes. However, four species were recently described as monocentric, yet our understanding of their genome evolution is largely masked by unclear phylogenetic relationships. Here, we reassess the current Juncus systematics using low-coverage genome skimming data of 33 taxa to construct repeats, nuclear rDNA and plastome-based phylogenetic hypothesis. Furthermore, we characterize the repeatome and chromosomal distribution of Juncus-specific centromeric repeats/CENH3 protein to test the monocentricity reach in the genus. Repeat-base phylogenies revealed topologies congruent with the rDNA tree, but not with the plastome tree. The incongruence between nuclear and plastome chloroplast dataset suggest an ancient hybridization in the divergence of Juncotypus and Tenageia sections 40 Myr ago. The phylogenetic resolution at section level was better fitted with the rDNA/repeat-based approaches, with the recognition of two monophyletic sections (Stygiopsis and Tenageia). We found specific repeatome trends for the main lineages, such as the higher abundances of TEs in the Caespitosi and Iridifolii + Ozophyllum clades. CENH3 immunostaining confirmed the monocentricity of Juncus, which can be a generic synapomorphy for the genus. The heterogeneity of the repeatomes, with high phylogenetic informativeness, identified here may be correlated with their ancient origin (56 Mya) and reveals the potential of comparative genomic analyses for understanding plant systematics and evolution.

Combined analysis of transposable elements and structural variation in maize genomes reveals genome contraction outpaces expansion

Structural differences between genomes are a major source of genetic variation that contributes to phenotypic differences. Transposable elements, mobile genetic sequences capable of increasing their copy number and propagating themselves within genomes, can generate structural variation. However, their repetitive nature makes it difficult to characterize fine-scale differences in their presence at specific positions, limiting our understanding of their impact on genome variation. Domesticated maize is a particularly good system for exploring the impact of transposable element proliferation as over 70% of the genome is annotated as transposable elements. High-quality transposable element annotations were recently generated for de novo genome assemblies of 26 diverse inbred maize lines. We generated base-pair resolved pairwise alignments between the B73 maize reference genome and the remaining 25 inbred maize line assemblies. From this data, we classified transposable elements as either shared or polymorphic in a given pairwise comparison. Our analysis uncovered substantial structural variation between lines, representing both simple and complex connections between TEs and structural variants. Putative insertions in SNP depleted regions, which represent recently diverged identity by state blocks, suggest some TE families may still be active. However, our analysis reveals that within these recently diverged genomic regions, deletions of transposable elements likely account for more structural variation events and base pairs than insertions. These deletions are often large structural variants containing multiple transposable elements. Combined, our results highlight how transposable elements contribute to structural variation and demonstrate that deletion events are a major contributor to genomic differences.

Identification of transcriptionally active transposons in Barley

BACKGROUND

The genomes of many major crops including barley (Hordeum vulgare) consist of numerous transposons. Despite their important roles in crop genome evolution and morphological variations, most of these elements are silent or truncated and unable to be mobile in host genomes. Thus far, only a very limited number of active transposons were identified in plants.

RESULTS

We analyzed the barley full-length cDNA (FLcDNA) sequences and detected 71 unique FLcDNAs exhibiting significant sequence similarity to the extant transposase proteins. These FLcDNAs were then used to search against the genome of a malting barley cultivar 'Morex', seven new intact transposons were identified. Sequence alignments indicated that six intact transposons contained the entire FLcDNAs whereas another one served as 3' untranslated region (3' UTR) of a barley gene. Our reverse transcription-PCR (RT-PCR) experiment further confirmed the expression of these six transposons and revealed their differential expression. We conducted genome-wide transposon comparisons and detected polymorphisms of three transposon families between the genomes of 'Morex' and other three genotypes including the wild barley (Hordeum spontaneum, B1K-04-12) and two cultivated barley varieties, 'Golden Promise' and 'Lasa Goumang'. Lastly, we screened the transcripts of all annotated barley genes and found that some transposons may serve as the coding regions (CDSs) or UTRs of barley genes.

CONCLUSION

We identified six newly expressed transposons in the barley genome and revealed the recent mobility of three transposon families. Our efforts provide a valuable resource for understanding the effects of transposons on barley genome evolution and for developing novel molecular tools for barley genetic improvement and other research.

miRNA-like secondary structures in maize (Zea mays) genes and transposable elements correlate with small RNAs, methylation, and expression

RNA molecules carry information in their primary sequence and also their secondary structure. Secondary structure can confer important functional information, but it is also a signal for an RNAi-like host epigenetic response mediated by small RNAs (smRNAs). In this study, we used two bioinformatic methods to predict local secondary structures across features of the maize genome, focusing on small regions that had similar folding properties to pre-miRNA loci. We found miRNA-like secondary structures to be common in genes and most, but not all, superfamilies of RNA and DNA transposable elements (TEs). The miRNA-like regions map to a higher diversity of smRNAs than regions without miRNA-like structure, explaining up to 27% of variation in smRNA mapping for some TE superfamilies. This mapping bias is more pronounced among putatively autonomous TEs relative to nonautonomous TEs. Genome-wide, miRNA-like regions are also associated with elevated methylation levels, particularly in the CHH context. Among genes, those with miRNA-like secondary structure are 1.5-fold more highly expressed, on average, than other genes. However, these genes are also more variably expressed across the 26 nested association mapping founder lines, and this variability positively correlates with the number of mapping smRNAs. We conclude that local miRNA-like structures are a nearly ubiquitous feature of expressed regions of the maize genome, that they correlate with higher smRNA mapping and methylation, and that they may represent a trade-off between functional requirements and the potentially negative consequences of smRNA production.

Transposable element evolution in plant genome ecosystems

The relationship of transposable elements (TEs) with their host genomes has usually been seen as an arms race between TEs and their host genomes. Consequently, TEs are supposed to amplify by bursts of transposition, when the TE escapes host surveillance, followed by long periods of TE quiescence and efficient host control. Recent data obtained from an increasing number of assembled plant genomes and resequencing population datasets show that TE dynamics is more complex and varies among TE families and their host genomes. This variation ranges from large genomes that accommodate large TE populations to genomes that are very active in TE elimination, and from inconspicuous elements with very low activity to elements with high transposition and elimination rates. The dynamics of each TE family results from a long history of interaction with the host in a genome populated by many other TE families, very much like an evolving ecosystem.

Analyses of 600+ insect genomes reveal repetitive element dynamics and highlight biodiversity-scale repeat annotation challenges

Repetitive elements (REs) are integral to the composition, structure, and function of eukaryotic genomes, yet remain understudied in most taxonomic groups. We investigated REs across 601 insect species and report wide variation in RE dynamics across groups. Analysis of associations between REs and protein-coding genes revealed dynamic evolution at the interface between REs and coding regions across insects, including notably elevated RE-gene associations in lineages with abundant long interspersed nuclear elements (LINEs). We leveraged this large, empirical data set to quantify impacts of long-read technology on RE detection and investigate fundamental challenges to RE annotation in diverse groups. In long-read assemblies, we detected ∼36% more REs than short-read assemblies, with long terminal repeats (LTRs) showing 162% increased detection, whereas DNA transposons and LINEs showed less respective technology-related bias. In most insect lineages, 25%-85% of repetitive sequences were "unclassified" following automated annotation, compared with only ∼13% in Drosophila species. Although the diversity of available insect genomes has rapidly expanded, we show the rate of community contributions to RE databases has not kept pace, preventing efficient annotation and high-resolution study of REs in most groups. We highlight the tremendous opportunity and need for the biodiversity genomics field to embrace REs and suggest collective steps for making progress toward this goal.

All families of transposable elements were active in the recent wheat genome evolution and polyploidy had no impact on their activity

Bread wheat (Triticum aestivum L.) is a major crop and its genome is one of the largest ever assembled at reference-quality level. It is 15 Gb, hexaploid, with 85% of transposable elements (TEs). Wheat genetic diversity was mainly focused on genes and little is known about the extent of genomic variability affecting TEs, transposition rate, and the impact of polyploidy. Multiple chromosome-scale assemblies are now available for bread wheat and for its tetraploid and diploid wild relatives. In this study, we computed base pair-resolved, gene-anchored, whole genome alignments of A, B, and D lineages at different ploidy levels in order to estimate the variability that affects the TE space. We used assembled genomes of 13 T. aestivum cultivars (6x = AABBDD) and a single genome for Triticum durum (4x = AABB), Triticum dicoccoides (4x = AABB), Triticum urartu (2x = AA), and Aegilops tauschii (2x = DD). We show that 5%-34% of the TE fraction is variable, depending on the species divergence. Between 400 and 13,000 novel TE insertions per subgenome were detected. We found lineage-specific insertions for nearly all TE families in di-, tetra-, and hexaploids. No burst of transposition was observed and polyploidization did not trigger any boost of transposition. This study challenges the prevailing idea of wheat TE dynamics and is more in agreement with an equilibrium model of evolution.

Genome-wide characterization of long terminal repeat retrotransposons provides insights into trait evolution of four cucurbit species

Cucurbits are a diverse plant family that includes economically important crops, such as cucumber, watermelon, melon, and pumpkin. Knowledge of the roles that long terminal repeat retrotransposons (LTR-RTs) have played in diversification of cucurbit species is limited; to add to understanding of the roles of LTR-RTs, we assessed their distributions in four cucurbit species. We identified 381, 578, 1086, and 623 intact LTR-RTs in cucumber (Cucumis sativus L. var. sativus cv. Chinese Long), watermelon (Citrullus lanatus subsp. vulgaris cv. 97103), melon (Cucumis melo cv. DHL92), and Cucurbita (Cucurbita moschata var. Rifu), respectively. Among these LTR-RTs, the Ale clade of the Copia superfamily was the most abundant in all the four cucurbit species. Insertion time and copy number analysis revealed that an LTR-RT burst occurred approximately 2 million years ago in cucumber, watermelon, melon, and Cucurbita, and may have contributed to their genome size variation. Phylogenetic and nucleotide polymorphism analyses suggested that most LTR-RTs were formed after species diversification. Analysis of gene insertions by LTR-RTs revealed that the most frequent insertions were of Ale and Tekay and that genes related to dietary fiber synthesis were the most commonly affected by LTR-RTs in Cucurbita. These results increase our understanding of LTR-RTs and their roles in genome evolution and trait characterization in cucurbits.

Pol V produced RNA facilitates transposable element excision site repair in Arabidopsis

The plant-specific RNA Polymerase V (Pol V) plays a key role in gene silencing, but its role in repair of double stranded DNA breaks is unclear. Excision of the transposable element mPing creates double stranded breaks that are repaired by NHEJ. We measured mPing excision site repair in multiple DNA methylation mutants including pol V using an mPing : GFP reporter. Two independent mutant alleles of pol V showed less GFP expression, indicating that the Pol V protein plays a role in excision site repair. Sequence analysis of the pol V excision sites indicated an elevated rate of large deletions consistent with less efficient repair. These results clarify the role of Pol V, but not other RNA-directed DNA methylation proteins (Pol IV) or maintenance DNA methylation pathways ( MET1 ), in the repair of double-strand DNA breaks.

Recent transposable element bursts are associated with the proximity to genes in a fungal plant pathogen

The activity of transposable elements (TEs) contributes significantly to pathogen genome evolution. TEs often destabilize genome integrity but may also confer adaptive variation in pathogenicity or resistance traits. De-repression of epigenetically silenced TEs often initiates bursts of transposition activity that may be counteracted by purifying selection and genome defenses. However, how these forces interact to determine the expansion routes of TEs within a pathogen species remains largely unknown. Here, we analyzed a set of 19 telomere-to-telomere genomes of the fungal wheat pathogen Zymoseptoria tritici. Phylogenetic reconstruction and ancestral state estimates of individual TE families revealed that TEs have undergone distinct activation and repression periods resulting in highly uneven copy numbers between genomes of the same species. Most TEs are clustered in gene poor niches, indicating strong purifying selection against insertions near coding sequences, or as a consequence of insertion site preferences. TE families with high copy numbers have low sequence divergence and strong signatures of defense mechanisms (i.e., RIP). In contrast, small non-autonomous TEs (i.e., MITEs) are less impacted by defense mechanisms and are often located in close proximity to genes. Individual TE families have experienced multiple distinct burst events that generated many nearly identical copies. We found that a Copia element burst was initiated from recent copies inserted substantially closer to genes compared to older copies. Overall, TE bursts tended to initiate from copies in GC-rich niches that escaped inactivation by genomic defenses. Our work shows how specific genomic environments features provide triggers for TE proliferation in pathogen genomes.

Intracellular common gardens reveal niche differentiation in transposable element community during bacterial adaptive evolution

The distribution and abundance of transposable elements across the tree of life have significantly shaped the evolution of cellular organisms, but the underlying mechanisms shaping these ecological patterns remain elusive. Here we establish a "common garden" approach to study causal ecological interactions between a xenogeneic conditional lethal sacB gene and the community of transposable insertion sequences (ISs) in a multipartite prokaryote genome. Xenogeneic sacB of low, medium, or high GC content was individually inserted into three replicons of a model bacterium Sinorhizobium fredii, and exhibited replicon- and GC-dependent variation in genetic stability. This variation was largely attributable to multidimensional niche differentiation for IS community members. The transposition efficiency of major active ISs depended on the nucleoid-associated xenogeneic silencer MucR. Experimentally eliminating insertion activity of specific ISs by deleting MucR strongly demonstrated a dominant role of niche differentiation among ISs. This intracellular common garden approach in the experimental evolution context allows not only for evaluating genetic stability of natural and synthetic xenogeneic genes of different sequence signatures in host cells but also for tracking and testing causal relationships in unifying ecological principles in genome ecology.

The role of transposon inverted repeats in balancing drought tolerance and yield-related traits in maize

The genomic basis underlying the selection for environmental adaptation and yield-related traits in maize remains poorly understood. Here we carried out genome-wide profiling of the small RNA (sRNA) transcriptome (sRNAome) and transcriptome landscapes of a global maize diversity panel under dry and wet conditions and uncover dozens of environment-specific regulatory hotspots. Transgenic and molecular studies of Drought-Related Environment-specific Super eQTL Hotspot on chromosome 8 (DRESH8) and ZmMYBR38, a target of DRESH8-derived small interfering RNAs, revealed a transposable element-mediated inverted repeats (TE-IR)-derived sRNA- and gene-regulatory network that balances plant drought tolerance with yield-related traits. A genome-wide scan revealed that TE-IRs associate with drought response and yield-related traits that were positively selected and expanded during maize domestication. These results indicate that TE-IR-mediated posttranscriptional regulation is a key molecular mechanism underlying the tradeoff between crop environmental adaptation and yield-related traits, providing potential genomic targets for the breeding of crops with greater stress tolerance but uncompromised yield.

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.

ZmDREB2.9 Gene in Maize (Zea mays L.): Genome-Wide Identification, Characterization, Expression, and Stress Response

Dehydration-responsive element-binding (DREB) transcription factors of the A2 subfamily play key roles in plant stress responses. In this study, we identified and characterized a new A2-type DREB gene, ZmDREB2.9, in the Zea mays cv. B73 genome and compared its expression profile with those of the known A2-type maize genes ZmDREB2.1-2.8. ZmDREB2.9 was mapped to chromosome 8, contained 18 predicted hormone- and stress-responsive cis-elements in the promoter, and had two splice isoforms: short ZmDREB2.9-S preferentially expressed in the leaves, embryos, and endosperm and long ZmDREB2.9-L expressed mostly in the male flowers, stamens, and ovaries. Phylogenetically, ZmDREB2.9 was closer to A. thaliana DREB2A than the other ZmDREB2 factors. ZmDREB2.9-S, ZmDREB2.2, and ZmDREB2.1/2A were upregulated in response to cold, drought, and abscisic acid and may play redundant roles in maize stress resistance. ZmDREB2.3, ZmDREB2.4, and ZmDREB2.6 were not expressed in seedlings and could be pseudogenes. ZmDREB2.7 and ZmDREB2.8 showed similar transcript accumulation in response to cold and abscisic acid and could be functionally redundant. Our results provide new data on Z. mays DREB2 factors, which can be used for further functional studies as well as in breeding programs to improve maize stress tolerance.

Environmental response in gene expression and DNA methylation reveals factors influencing the adaptive potential of Arabidopsis lyrata

Understanding what factors influence plastic and genetic variation is valuable for predicting how organisms respond to changes in the selective environment. Here, using gene expression and DNA methylation as molecular phenotypes, we study environmentally induced variation among Arabidopsis lyrata plants grown at lowland and alpine field sites. Our results show that gene expression is highly plastic, as many more genes are differentially expressed between the field sites than between populations. These environmentally responsive genes evolve under strong selective constraint – the strength of purifying selection on the coding sequence is high, while the rate of adaptive evolution is low. We find, however, that positive selection on cis-regulatory variants has likely contributed to the maintenance of genetically variable environmental responses, but such variants segregate only between distantly related populations. In contrast to gene expression, DNA methylation at genic regions is largely insensitive to the environment, and plastic methylation changes are not associated with differential gene expression. Besides genes, we detect environmental effects at transposable elements (TEs): TEs at the high-altitude field site have higher expression and methylation levels, suggestive of a broad-scale TE activation. Compared to the lowland population, plants native to the alpine environment harbor an excess of recent TE insertions, and we observe that specific TE families are enriched within environmentally responsive genes. Our findings provide insight into selective forces shaping plastic and genetic variation. We also highlight how plastic responses at TEs can rapidly create novel heritable variation in stressful conditions.

A comprehensive map of preferentially located motifs reveals distinct proximal cis-regulatory sequences in plants

Identification of cis-regulatory sequences controlling gene expression is an arduous challenge that is being actively explored to discover key genetic factors responsible for traits of agronomic interest. Here, we used a genome-wide de novo approach to investigate preferentially located motifs (PLMs) in the proximal cis-regulatory landscape of Arabidopsis thaliana and Zea mays. We report three groups of PLMs in both the 5'- and 3'-gene-proximal regions and emphasize conserved PLMs in both species, particularly in the 3'-gene-proximal region. Comparison with resources from transcription factor and microRNA binding sites shows that 79% of the identified PLMs are unassigned, although some are supported by MNase-defined cistrome occupancy analysis. Enrichment analyses further reveal that unassigned PLMs provide functional predictions that differ from those derived from transcription factor and microRNA binding sites. Our study provides a comprehensive map of PLMs and demonstrates their potential utility for future characterization of orphan genes in plants.

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'.

Prediction of evolutionary constraint by genomic annotations improves functional prioritization of genomic variants in maize

Background

Crop improvement through cross-population genomic prediction and genome editing requires identification of causal variants at high resolution, within fewer than hundreds of base pairs. Most genetic mapping studies have generally lacked such resolution. In contrast, evolutionary approaches can detect genetic effects at high resolution, but they are limited by shifting selection, missing data, and low depth of multiple-sequence alignments. Here we use genomic annotations to accurately predict nucleotide conservation across angiosperms, as a proxy for fitness effect of mutations.

Results

Using only sequence analysis, we annotate nonsynonymous mutations in 25,824 maize gene models, with information from bioinformatics and deep learning. Our predictions are validated by experimental information: within-species conservation, chromatin accessibility, and gene expression. According to gene ontology and pathway enrichment analyses, predicted nucleotide conservation points to genes in central carbon metabolism. Importantly, it improves genomic prediction for fitness-related traits such as grain yield, in elite maize panels, by stringent prioritization of fewer than 1% of single-site variants.

Conclusions

Our results suggest that predicting nucleotide conservation across angiosperms may effectively prioritize sites most likely to impact fitness-related traits in crops, without being limited by shifting selection, missing data, and low depth of multiple-sequence alignments. Our approach—Prediction of mutation Impact by Calibrated Nucleotide Conservation (PICNC)—could be useful to select polymorphisms for accurate genomic prediction, and candidate mutations for efficient base editing. The trained PICNC models and predicted nucleotide conservation at protein-coding SNPs in maize are publicly available in CyVerse (10.25739/hybz-2957).

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02747-2.

The widespread nature of Pack-TYPE transposons reveals their importance for plant genome evolution

Pack-TYPE transposable elements (TEs) are a group of non-autonomous DNA transposons found in plants. These elements can efficiently capture and shuffle coding DNA across the host genome, accelerating the evolution of genes. Despite their relevance for plant genome plasticity, the detection and study of Pack-TYPE TEs are challenging due to the high similarity these elements have with genes. Here, we produced an automated annotation pipeline designed to study Pack-TYPE elements and used it to successfully annotate and analyse more than 10,000 new Pack-TYPE TEs in the rice and maize genomes. Our analysis indicates that Pack-TYPE TEs are an abundant and heterogeneous group of elements. We found that these elements are associated with all main superfamilies of Class II DNA transposons in plants and likely share a similar mechanism to capture new chromosomal DNA sequences. Furthermore, we report examples of the direct contribution of these TEs to coding genes, suggesting a generalised and extensive role of Pack-TYPE TEs in plant genome evolution.

Transposable Elements (TEs) are genetic DNA sequences able to move across the genome, and their transposition activity is associated with genome plasticity and gene evolution. However, most of these elements exhibit “selfish” behaviour, meaning that they mainly transpose their own DNA sequence and only exceptionally might rearrange the DNA of coding genes. Pack-TYPE TEs, found in plants, represent an important exception, and they can efficiently capture and shuffle DNA sequences captured from the genome, accelerating the evolution of genes. We provide here the first automatic pipeline designed explicitly for the annotation of Pack-TYPE TEs. We used our approach to systematically investigate Pack-TYPE TEs in the rice and maize reference genomes, and annotated thousands of new elements in these species. We demonstrate that Pack-TYPE elements are abundant in plants and we report several examples of coding genes originated as a consequence of the mobilization of these elements.

Zebrafish transposable elements show extensive diversification in age, genomic distribution, and developmental expression

There is considerable interest in understanding the effect of transposable elements (TEs) on embryonic development. Studies in humans and mice are limited by the difficulty of working with mammalian embryos and by the relative scarcity of active TEs in these organisms. The zebrafish is an outstanding model for the study of vertebrate development, and over half of its genome consists of diverse TEs. However, zebrafish TEs remain poorly characterized. Here we describe the demography and genomic distribution of zebrafish TEs and their expression throughout embryogenesis using bulk and single-cell RNA sequencing data. These results reveal a highly dynamic genomic ecosystem comprising nearly 2000 distinct TE families, which vary in copy number by four orders of magnitude and span a wide range of ages. Longer retroelements tend to be retained in intergenic regions, whereas short interspersed nuclear elements (SINEs) and DNA transposons are more frequently found nearby or within genes. Locus-specific mapping of TE expression reveals extensive TE transcription during development. Although two-thirds of TE transcripts are likely driven by nearby gene promoters, we still observe stage- and tissue-specific expression patterns in self-regulated TEs. Long terminal repeat (LTR) retroelements are most transcriptionally active immediately following zygotic genome activation, whereas DNA transposons are enriched among transcripts expressed in later stages of development. Single-cell analysis reveals several endogenous retroviruses expressed in specific somatic cell lineages. Overall, our study provides a valuable resource for using zebrafish as a model to study the impact of TEs on vertebrate development.

AnchorWave: Sensitive alignment of genomes with high sequence diversity, extensive structural polymorphism, and whole-genome duplication

Millions of species are currently being sequenced, and their genomes are being compared. Many of them have more complex genomes than model systems and raise novel challenges for genome alignment. Widely used local alignment strategies often produce limited or incongruous results when applied to genomes with dispersed repeats, long indels, and highly diverse sequences. Moreover, alignment using many-to-many or reciprocal best hit approaches conflicts with well-studied patterns between species with different rounds of whole-genome duplication. Here, we introduce Anchored Wavefront alignment (AnchorWave), which performs whole-genome duplication-informed collinear anchor identification between genomes and performs base pair-resolved global alignment for collinear blocks using a two-piece affine gap cost strategy. This strategy enables AnchorWave to precisely identify multikilobase indels generated by transposable element (TE) presence/absence variants (PAVs). When aligning two maize genomes, AnchorWave successfully recalled 87% of previously reported TE PAVs. By contrast, other genome alignment tools showed low power for TE PAV recall. AnchorWave precisely aligns up to three times more of the genome as position matches or indels than the closest competitive approach when comparing diverse genomes. Moreover, AnchorWave recalls transcription factor-binding sites at a rate of 1.05- to 74.85-fold higher than other tools with significantly lower false-positive alignments. AnchorWave complements available genome alignment tools by showing obvious improvement when applied to genomes with dispersed repeats, active TEs, high sequence diversity, and whole-genome duplication variation.

OUP accepted manuscript

Recognition of the important role of transposable elements (TEs) in eukaryotic genomes quickly led to a burgeoning literature modeling and estimating the effects of selection on TEs. Much of the empirical work on selection has focused on analyzing the site frequency spectrum (SFS) of TEs. But TE evolution differs from standard models in a number of ways that can impact the power and interpretation of the SFS. For example, rather than mutating under a clock-like model, transposition often occurs in bursts which can inflate particular frequency categories compared with expectations under a standard neutral model. If a TE burst has been recent, the excess of low-frequency polymorphisms can mimic the effect of purifying selection. Here, we investigate how transposition bursts affect the frequency distribution of TEs and the correlation between age and allele frequency. Using information on the TE age distribution, we propose an age-adjusted SFS to compare TEs and neutral polymorphisms to more effectively evaluate whether TEs are under selective constraints. We show that our approach can minimize instances of false inference of selective constraint, remains robust to simple demographic changes, and allows for a correct identification of even weak selection affecting TEs which experienced a transposition burst. The results presented here will help researchers working on TEs to more reliably identify the effects of selection on TEs without having to rely on the assumption of a constant transposition rate.

Stability of DNA methylation and chromatin accessibility in structurally diverse maize genomes

Accessible chromatin and unmethylated DNA are associated with many genes and cis-regulatory elements. Attempts to understand natural variation for accessible chromatin regions (ACRs) and unmethylated regions (UMRs) often rely upon alignments to a single reference genome. This limits the ability to assess regions that are absent in the reference genome assembly and monitor how nearby structural variants influence variation in chromatin state. In this study, de novo genome assemblies for four maize inbreds (B73, Mo17, Oh43, and W22) are utilized to assess chromatin accessibility and DNA methylation patterns in a pan-genome context. A more complete set of UMRs and ACRs can be identified when chromatin data are aligned to the matched genome rather than a single reference genome. While there are UMRs and ACRs present within genomic regions that are not shared between genotypes, these features are 6- to 12-fold enriched within regions between genomes. Characterization of UMRs present within shared genomic regions reveals that most UMRs maintain the unmethylated state in other genotypes with only ∼5% being polymorphic between genotypes. However, the majority (71%) of UMRs that are shared between genotypes only exhibit partial overlaps suggesting that the boundaries between methylated and unmethylated DNA are dynamic. This instability is not solely due to sequence variation as these partially overlapping UMRs are frequently found within genomic regions that lack sequence variation. The ability to compare chromatin properties among individuals with structural variation enables pan-epigenome analyses to study the sources of variation for accessible chromatin and unmethylated DNA.