Plant centromeric retrotransposons: a structural and cytogenetic perspective

Comparative Analysis of Active LTR Retrotransposons in Sunflower (Helianthus annuus L.): From Extrachromosomal Circular DNA Detection to Protein Structure Prediction

Plant genomes possess numerous transposable element (TE) insertions that have occurred during evolution. Most TEs are silenced or diverged; therefore, they lose their ability to encode proteins and are transposed in the genome. Knowledge of active plant TEs and TE-encoded proteins essential for transposition and evasion of plant cell transposon silencing mechanisms remains limited. This study investigated active long terminal repeat (LTR) retrotransposons (RTEs) in sunflowers (Helianthus annuus), revealing heterogeneous and phylogenetically distinct RTEs triggered by epigenetic changes and heat stress. Many of these RTEs belong to three distinct groups within the Tekay clade, showing significant variations in chromosomal insertion distribution. Through protein analysis of these active RTEs, it was found that Athila RTEs and Tekay group 2 elements possess additional open reading frames (aORFs). The aORF-encoded proteins feature a transposase domain, a transmembrane domain, and nuclear localization signals. The aORF proteins of the Tekay subgroup exhibited remarkable conservation among over 500 Tekay members, suggesting their functional importance in RTE mobility. The predicted 3D structure of the sunflower Tekay aORF protein showed significant homology with Tekay proteins in rice, maize, and sorghum. Additionally, the structural features of aORF proteins resemble those of plant DRBM-containing proteins, suggesting their potential role in RNA-silencing modulation. These findings offer insights into the diversity and activity of sunflower RTEs, emphasizing the conservation and structural characteristics of aORF-encoded proteins.

Genomes of autotetraploid wild and cultivated Ziziphus mauritiana reveal polyploid evolution and crop domestication

Indian jujube (Ziziphus mauritiana) holds a prominent position in the global fruit and pharmaceutical markets. Here, we report the assemblies of haplotype-resolved, telomere-to-telomere genomes of autotetraploid wild and cultivated Indian jujube plants using a 2-stage assembly strategy. The generation of these genomes permitted in-depth investigations into the divergence and evolutionary history of this important fruit crop. Using a graph-based pan-genome constructed from 8 monoploid genomes, we identified structural variation (SV)-FST hotspots and SV hotspots. Gap-free genomes provide a means to obtain a global view of centromere structures. We identified presence-absence variation-related genes in 4 monoploid genomes (cI, cIII, wI, and wIII) and resequencing populations. We also present the population structure and domestication trajectory of the Indian jujube based on the resequencing of 73 wild and cultivated accessions. Metabolomic and transcriptomic analyses of mature fruits of wild and cultivated accessions unveiled the genetic basis underlying loss of fruit astringency during domestication of Indian jujube. This study reveals mechanisms underlying the divergence, evolution, and domestication of the autotetraploid Indian jujube and provides rich and reliable genetic resources for future research.

ZW sex chromosome structure in Amborella trichopoda

Sex chromosomes have evolved hundreds of times across the flowering plant tree of life; their recent origins in some members of this clade can shed light on the early consequences of suppressed recombination, a crucial step in sex chromosome evolution. Amborella trichopoda, the sole species of a lineage that is sister to all other extant flowering plants, is dioecious with a young ZW sex determination system. Here we present a haplotype-resolved genome assembly, including highly contiguous assemblies of the Z and W chromosomes. We identify a ~3-megabase sex-determination region (SDR) captured in two strata that includes a ~300-kilobase inversion that is enriched with repetitive sequences and contains a homologue of the Arabidopsis METHYLTHIOADENOSINE NUCLEOSIDASE (MTN1-2) genes, which are known to be involved in fertility. However, the remainder of the SDR does not show patterns typically found in non-recombining SDRs, such as repeat accumulation and gene loss. These findings are consistent with the hypothesis that dioecy is derived in Amborella and the sex chromosome pair has not significantly degenerated.

Chromosome-specific barcode system with centromeric repeat in cultivated soybean and wild progenitor

Wild soybean Glycine soja is the progenitor of cultivated soybean Glycine max Information on soybean functional centromeres is limited despite extensive genome analysis. These species are an ideal model for studying centromere dynamics for domestication and breeding. We performed a detailed chromatin immunoprecipitation analysis using centromere-specific histone H3 protein to delineate two distinct centromeric DNA sequences with unusual repeating units with monomer sizes of 90-92 bp (CentGm-1) and 413-bp (CentGm-4) shorter and longer than standard nucleosomes. These two unrelated DNA sequences with no sequence similarity are part of functional centromeres in both species. Our results provide a comparison of centromere properties between a cultivated and a wild species under the effect of the same kinetochore protein. Possible sequence homogenization specific to each chromosome could highlight the mechanism for evolutionary conservation of centromeric properties independent of domestication and breeding. Moreover, a unique barcode system to track each chromosome is developed using CentGm-4 units. Our results with a unifying centromere composition model using CentGm-1 and CentGm-4 superfamilies could have far-reaching implications for comparative and evolutionary genome research.

Allotetraploid nature of a wild potato species, Solanum stoloniferum Schlechtd. et Bché., as revealed by whole-genome sequencing

Mexican wild diploid potato species are reproductively isolated from A-genome species, including cultivated potatoes; thus, their genomic relationships remain unknown. Solanum stoloniferum Schlechtd. et Bché. (2n = 4x = 48, AABB) is a Mexican allotetraploid species frequently used in potato breeding. We constructed a chromosome-scale assembly of the S. stoloniferum genome using PacBio long-read sequencing and Hi-C scaffolding technologies. The final assembly consisted of 1742 Mb, among which 745 Mb and 713 Mb were anchored to the 12 A-genome and 12 B-genome chromosomes, respectively. Using the RNA-seq datasets, we detected 20 994 and 19 450 genes in the A and B genomes, respectively. Among these genes, 5138 and 3594 were specific to the A and B genomes, respectively, and 15 856 were homoeologous, of which 18.6-25.4% were biasedly expressed. Structural variations such as large pericentromeric inversions were frequently found between the A- and B-genome chromosomes. A comparison of the gene sequences from 38 diverse genomes of the related Solanum species revealed that the S. stoloniferum B genome and Mexican diploid species, with the exception of S. verrucosum, were monophyletically distinct from the S. stoloniferum A genome and the other A-genome species, indicating that the Mexican diploid species share the B genome. The content and divergence of transposable elements (TEs) revealed recent bursts and transpositions of TEs after polyploidization. Thus, the S. stoloniferum genome has undergone dynamic structural differentiation and TE mobilization and reorganization to stabilize the genomic imbalance. This study provides new insights into polyploid evolution and the efficient use of allotetraploid species in potato breeding.

3D chromatin maps of a brown alga reveal U/V sex chromosome spatial organization

Nuclear three dimensional (3D) folding of chromatin structure has been linked to gene expression regulation and correct developmental programs, but little is known about the 3D architecture of sex chromosomes within the nucleus, and how that impacts their role in sex determination. Here, we determine the sex-specific 3D organization of the model brown alga Ectocarpus chromosomes at 2 kb resolution, by mapping long-range chromosomal interactions using Hi-C coupled with Oxford Nanopore long reads. We report that Ectocarpus interphase chromatin exhibits a non-Rabl conformation, with strong contacts among telomeres and among centromeres, which feature centromere-specific LTR retrotransposons. The Ectocarpus chromosomes do not contain large local interactive domains that resemble TADs described in animals, but their 3D genome organization is largely shaped by post-translational modifications of histone proteins. We show that the sex determining region (SDR) within the U and V chromosomes are insulated and span the centromeres and we link sex-specific chromatin dynamics and gene expression levels to the 3D chromatin structure of the U and V chromosomes. Finally, we uncover the unique conformation of a large genomic region on chromosome 6 harboring an endogenous viral element, providing insights regarding the impact of a latent giant dsDNA virus on the host genome's 3D chromosomal folding.

Repeat-based holocentromeres of the woodrush Luzula sylvatica reveal insights into the evolutionary transition to holocentricity

In most studied eukaryotes, chromosomes are monocentric, with centromere activity confined to a single region. However, the rush family (Juncaceae) includes species with both monocentric (Juncus) and holocentric (Luzula) chromosomes, where centromere activity is distributed along the entire chromosome length. Here, we combine chromosome-scale genome assembly, epigenetic analysis, immuno-FISH and super-resolution microscopy to study the transition to holocentricity in Luzula sylvatica. We report repeat-based holocentromeres with an irregular distribution of features along the chromosomes. Luzula sylvatica holocentromeres are predominantly associated with two satellite DNA repeats (Lusy1 and Lusy2), while CENH3 also binds satellite-free gene-poor regions. Comparative repeat analysis suggests that Lusy1 plays a crucial role in centromere function across most Luzula species. Furthermore, synteny analysis between L. sylvatica (n = 6) and Juncus effusus (n = 21) suggests that holocentric chromosomes in Luzula could have arisen from chromosome fusions of ancestral monocentric chromosomes, accompanied by the expansion of CENH3-associated satellite repeats.

DANTE and DANTE_LTR: lineage-centric annotation pipelines for long terminal repeat retrotransposons in plant genomes

Long terminal repeat (LTR) retrotransposons constitute a predominant class of repetitive DNA elements in most plant genomes. With the increasing number of sequenced plant genomes, there is an ongoing demand for computational tools facilitating efficient annotation and classification of LTR retrotransposons in plant genome assemblies. Herein, we introduce DANTE, a computational pipeline for Domain-based ANnotation of Transposable Elements, designed for sensitive detection of these elements via their conserved protein domain sequences. The identified protein domains are subsequently inputted into the DANTE_LTR pipeline to annotate complete element sequences by detecting their structural features, such as LTRs, in adjacent genomic regions. Leveraging domain sequences allows for precise classification of elements into phylogenetic lineages, offering a more granular annotation compared with coarser conventional superfamily-based classification methods. The efficiency and accuracy of this approach were evidenced via annotation of LTR retrotransposons in 93 plant genomes. Results were benchmarked against several established pipelines, showing that DANTE_LTR is capable of identifying significantly more intact LTR retrotransposons. DANTE and DANTE_LTR are provided as user-friendly Galaxy tools accessible via a public server (https://repeatexplorer-elixir.cerit-sc.cz), installable on local Galaxy instances from the Galaxy tool shed or executable from the command line.

Evolution of Einkorn wheat centromeres is driven by the mutualistic interplay of two LTR retrotransposons

BACKGROUND

Centromere function is highly conserved across eukaryotes, but the underlying centromeric DNA sequences vary dramatically between species. Centromeres often contain a high proportion of repetitive DNA, such as tandem repeats and/or transposable elements (TEs). Einkorn wheat centromeres lack tandem repeat arrays and are instead composed mostly of the two long terminal repeat (LTR) retrotransposon families RLG_Cereba and RLG_Quinta which specifically insert in centromeres. However, it is poorly understood how these two TE families relate to each other and if and how they contribute to centromere function and evolution.

RESULTS

Based on conservation of diagnostic motifs (LTRs, integrase and primer binding site and polypurine-tract), we propose that RLG_Cereba and RLG_Quinta are a pair of autonomous and non-autonomous partners, in which the autonomous RLG_Cereba contributes all the proteins required for transposition, while the non-autonomous RLG_Quinta contributes GAG protein. Phylogenetic analysis of predicted GAG proteins showed that the RLG_Cereba lineage was present for at least 100 million years in monocotyledon plants. In contrast, RLG_Quinta evolved from RLG_Cereba between 28 and 35 million years ago in the common ancestor of oat and wheat. Interestingly, the integrase of RLG_Cereba is fused to a so-called CR-domain, which is hypothesized to guide the integrase to the functional centromere. Indeed, ChIP-seq data and TE population analysis show only the youngest subfamilies of RLG_Cereba and RLG_Quinta are found in the active centromeres. Importantly, the LTRs of RLG_Quinta and RLG_Cereba are strongly associated with the presence of the centromere-specific CENH3 histone variant. We hypothesize that the LTRs of RLG_Cereba and RLG_Quinta contribute to wheat centromere integrity by phasing and/or placing CENH3 nucleosomes, thus favoring their persistence in the competitive centromere-niche.

CONCLUSION

Our data show that RLG_Cereba cross-mobilizes the non-autonomous RLG_Quinta retrotransposons. New copies of both families are specifically integrated into functional centromeres presumably through direct binding of the integrase CR domain to CENH3 histone variants. The LTRs of newly inserted RLG_Cereba and RLG_Quinta elements, in turn, recruit and/or phase new CENH3 deposition. This mutualistic interplay between the two TE families and the plant host dynamically maintains wheat centromeres.

Celine, a long interspersed nuclear element retrotransposon, colonizes in the centromeres of poplar chromosomes

Centromeres in most multicellular eukaryotes are composed of long arrays of repetitive DNA sequences. Interestingly, several transposable elements, including the well-known long terminal repeat centromeric retrotransposon of maize (CRM), were found to be enriched in functional centromeres marked by the centromeric histone H3 (CENH3). Here, we report a centromeric long interspersed nuclear element (LINE), Celine, in Populus species. Celine has colonized preferentially in the CENH3-associated chromatin of every poplar chromosome, with 84% of the Celine elements localized in the CENH3-binding domains. In contrast, only 51% of the CRM elements were bound to CENH3 domains in Populus trichocarpa. These results suggest different centromere targeting mechanisms employed by Celine and CRM elements. Nevertheless, the high target specificity seems to be detrimental to further amplification of the Celine elements, leading to a shorter life span and patchy distribution among plant species compared with the CRM elements. Using a phylogenetically guided approach, we were able to identify Celine-like LINE elements in tea plant (Camellia sinensis) and green ash tree (Fraxinus pennsylvanica). The centromeric localization of these Celine-like LINEs was confirmed in both species. We demonstrate that the centromere targeting property of Celine-like LINEs is of primitive origin and has been conserved among distantly related plant species.

Insights into spinach domestication from genome sequences of two wild spinach progenitors, Spinacia turkestanica and Spinacia tetrandra

Cultivated spinach (Spinacia oleracea) is a dioecious species. We report high-quality genome sequences for its two closest wild relatives, Spinacia turkestanica and Spinacia tetrandra, which are also dioecious, and are used to study the genetics of spinach domestication. Using a combination of genomic approaches, we assembled genomes of both these species and analyzed them in comparison with the previously assembled S. oleracea genome. These species diverged c. 6.3 million years ago (Ma), while cultivated spinach split from S. turkestanica 0.8 Ma. In all three species, all six chromosomes include very large gene-poor, repeat-rich regions, which, in S. oleracea, are pericentromeric regions with very low recombination rates in both male and female genetic maps. We describe population genomic evidence that the similar regions in the wild species also recombine rarely. We characterized 282 structural variants (SVs) that have been selected during domestication. These regions include genes associated with leaf margin type and flowering time. We also describe evidence that the downy mildew resistance loci of cultivated spinach are derived from introgression from both wild spinach species. Collectively, this study reveals the genome architecture of spinach assemblies and highlights the importance of SVs during the domestication of cultivated spinach.

Non-B-form DNA is associated with centromere stability in newly-formed polyploid wheat

Non-B-form DNA differs from the classic B-DNA double helix structure and plays a crucial regulatory role in replication and transcription. However, the role of non-B-form DNA in centromeres, especially in polyploid wheat, remains elusive. Here, we systematically analyzed seven non-B-form DNA motif profiles (A-phased DNA repeat, direct repeat, G-quadruplex, inverted repeat, mirror repeat, short tandem repeat, and Z-DNA) in hexaploid wheat. We found that three of these non-B-form DNA motifs were enriched at centromeric regions, especially at the CENH3-binding sites, suggesting that non-B-form DNA may create a favorable loading environment for the CENH3 nucleosome. To investigate the dynamics of centromeric non-B form DNA during the alloploidization process, we analyzed DNA secondary structure using CENH3 ChIP-seq data from newly formed allotetraploid wheat and its two diploid ancestors. We found that newly formed allotetraploid wheat formed more non-B-form DNA in centromeric regions compared with their parents, suggesting that non-B-form DNA is related to the localization of the centromeric regions in newly formed wheat. Furthermore, non-B-form DNA enriched in the centromeric regions was found to preferentially form on young LTR retrotransposons, explaining CENH3's tendency to bind to younger LTR. Collectively, our study describes the landscape of non-B-form DNA in the wheat genome, and sheds light on its potential role in the evolution of polyploid centromeres.

Two telomere-to-telomere gapless genomes reveal insights into Capsicum evolution and capsaicinoid biosynthesis

Chili pepper (Capsicum) is known for its unique fruit pungency due to the presence of capsaicinoids. The evolutionary history of capsaicinoid biosynthesis and the mechanism of their tissue specificity remain obscure due to the lack of high-quality Capsicum genomes. Here, we report two telomere-to-telomere (T2T) gap-free genomes of C. annuum and its wild nonpungent relative C. rhomboideum to investigate the evolution of fruit pungency in chili peppers. We precisely delineate Capsicum centromeres, which lack high-copy tandem repeats but are extensively invaded by CRM retrotransposons. Through phylogenomic analyses, we estimate the evolutionary timing of capsaicinoid biosynthesis. We reveal disrupted coding and regulatory regions of key biosynthesis genes in nonpungent species. We also find conserved placenta-specific accessible chromatin regions, which likely allow for tissue-specific biosynthetic gene coregulation and capsaicinoid accumulation. These T2T genomic resources will accelerate chili pepper genetic improvement and help to understand Capsicum genome evolution.

The structure, function, and evolution of plant centromeres

Centromeres are essential regions of eukaryotic chromosomes responsible for the formation of kinetochore complexes, which connect to spindle microtubules during cell division. Notably, although centromeres maintain a conserved function in chromosome segregation, the underlying DNA sequences are diverse both within and between species and are predominantly repetitive in nature. The repeat content of centromeres includes high-copy tandem repeats (satellites), and/or specific families of transposons. The functional region of the centromere is defined by loading of a specific histone 3 variant (CENH3), which nucleates the kinetochore and shows dynamic regulation. In many plants, the centromeres are composed of satellite repeat arrays that are densely DNA methylated and invaded by centrophilic retrotransposons. In some cases, the retrotransposons become the sites of CENH3 loading. We review the structure of plant centromeres, including monocentric, holocentric, and metapolycentric architectures, which vary in the number and distribution of kinetochore attachment sites along chromosomes. We discuss how variation in CENH3 loading can drive genome elimination during early cell divisions of plant embryogenesis. We review how epigenetic state may influence centromere identity and discuss evolutionary models that seek to explain the paradoxically rapid change of centromere sequences observed across species, including the potential roles of recombination. We outline putative modes of selection that could act within the centromeres, as well as the role of repeats in driving cycles of centromere evolution. Although our primary focus is on plant genomes, we draw comparisons with animal and fungal centromeres to derive a eukaryote-wide perspective of centromere structure and function.

The Dynamic Interplay Between Ribosomal DNA and Transposable Elements: A Perspective From Genomics and Cytogenetics

Although both are salient features of genomes, at first glance ribosomal DNAs and transposable elements are genetic elements with not much in common: whereas ribosomal DNAs are mainly viewed as housekeeping genes that uphold all prime genome functions, transposable elements are generally portrayed as selfish and disruptive. These opposing characteristics are also mirrored in other attributes: organization in tandem (ribosomal DNAs) versus organization in a dispersed manner (transposable elements); evolution in a concerted manner (ribosomal DNAs) versus evolution by diversification (transposable elements); and activity that prolongs genomic stability (ribosomal DNAs) versus activity that shortens it (transposable elements). Re-visiting relevant instances in which ribosomal DNA-transposable element interactions have been reported, we note that both repeat types share at least four structural and functional hallmarks: (1) they are repetitive DNAs that shape genomes in evolutionary timescales, (2) they exchange structural motifs and can enter co-evolution processes, (3) they are tightly controlled genomic stress sensors playing key roles in senescence/aging, and (4) they share common epigenetic marks such as DNA methylation and histone modification. Here, we give an overview of the structural, functional, and evolutionary characteristics of both ribosomal DNAs and transposable elements, discuss their roles and interactions, and highlight trends and future directions as we move forward in understanding ribosomal DNA-transposable element associations.

High-quality Gossypium hirsutum and Gossypium barbadense genome assemblies reveal the landscape and evolution of centromeres

Centromere positioning and organization are crucial for genome evolution; however, research on centromere biology is largely influenced by the quality of available genome assemblies. Here, we combined Oxford Nanopore and Pacific Biosciences technologies to de novo assemble two high-quality reference genomes for Gossypium hirsutum (TM-1) and Gossypium barbadense (3-79). Compared with previously published reference genomes, our assemblies show substantial improvements, with the contig N50 improved by 4.6-fold and 5.6-fold, respectively, and thus represent the most complete cotton genomes to date. These high-quality reference genomes enable us to characterize 14 and 5 complete centromeric regions for G. hirsutum and G. barbadense, respectively. Our data revealed that the centromeres of allotetraploid cotton are occupied by members of the centromeric repeat for maize (CRM) and Tekay long terminal repeat families, and the CRM family reshapes the centromere structure of the At subgenome after polyploidization. These two intertwined families have driven the convergent evolution of centromeres between the two subgenomes, ensuring centromere function and genome stability. In addition, the repositioning and high sequence divergence of centromeres between G. hirsutum and G. barbadense have contributed to speciation and centromere diversity. This study sheds light on centromere evolution in a significant crop and provides an alternative approach for exploring the evolution of polyploid plants.

Evolving Together: Cassandra Retrotransposons Gradually Mirror Promoter Mutations of the 5S rRNA Genes

The 5S rRNA genes are among the most conserved nucleotide sequences across all species. Similar to the 5S preservation we observe the occurrence of 5S-related nonautonomous retrotransposons, so-called Cassandras. Cassandras harbor highly conserved 5S rDNA-related sequences within their long terminal repeats, advantageously providing them with the 5S internal promoter. However, the dynamics of Cassandra retrotransposon evolution in the context of 5S rRNA gene sequence information and structural arrangement are still unclear, especially: (1) do we observe repeated or gradual domestication of the highly conserved 5S promoter by Cassandras and (2) do changes in 5S organization such as in the linked 35S-5S rDNA arrangements impact Cassandra evolution? Here, we show evidence for gradual co-evolution of Cassandra sequences with their corresponding 5S rDNAs. To follow the impact of 5S rDNA variability on Cassandra TEs, we investigate the Asteraceae family where highly variable 5S rDNAs, including 5S promoter shifts and both linked and separated 35S-5S rDNA arrangements have been reported. Cassandras within the Asteraceae mirror 5S rDNA promoter mutations of their host genome, likely as an adaptation to the host's specific 5S transcription factors and hence compensating for evolutionary changes in the 5S rDNA sequence. Changes in the 5S rDNA sequence and in Cassandras seem uncorrelated with linked/separated rDNA arrangements. We place all these observations into the context of angiosperm 5S rDNA-Cassandra evolution, discuss Cassandra's origin hypotheses (single or multiple) and Cassandra's possible impact on rDNA and plant genome organization, giving new insights into the interplay of ribosomal genes and transposable elements.

Near telomere-to-telomere genome of the model plant Physcomitrium patens

The model plant Physcomitrium patens has played a pivotal role in enhancing our comprehension of plant evolution and development. However, the current genome harbours numerous regions that remain unfinished and erroneous. To address these issues, we generated an assembly using Oxford Nanopore reads and Hi-C mapping. The assembly incorporates telomeric and centromeric regions, thereby establishing it as a near telomere-to-telomere genome except a region in chromosome 1 that is not fully assembled due to its highly repetitive nature. This near telomere-to-telomere genome resolves the chromosome number at 26 and provides a gap-free genome assembly as well as updated gene models to aid future studies using this model organism.

Transposable elements: multifunctional players in the plant genome

Transposable elements (TEs) are indispensable components of eukaryotic genomes that play diverse roles in gene regulation, recombination, and environmental adaptation. Their ability to mobilize within the genome leads to gene expression and DNA structure changes. TEs serve as valuable markers for genetic and evolutionary studies and facilitate genetic mapping and phylogenetic analysis. They also provide insight into how organisms adapt to a changing environment by promoting gene rearrangements that lead to new gene combinations. These repetitive sequences significantly impact genome structure, function and evolution. This review takes a comprehensive look at TEs and their applications in biotechnology, particularly in the context of plant biology, where they are now considered "genomic gold" due to their extensive functionalities. The article addresses various aspects of TEs in plant development, including their structure, epigenetic regulation, evolutionary patterns, and their use in gene editing and plant molecular markers. The goal is to systematically understand TEs and shed light on their diverse roles in plant biology.

The telomere-to-telomere (T2T) genome provides insights into the evolution of specialized centromere sequences in sandalwood

BACKGROUND

Sandalwood, a prized hemiparasitic plant, is highly sought in the commercial market because of its aromatic core materia. The structure and stability of the genome are instrumental in the rapid adaptation of parasitic plants to their surroundings. However, there is a conspicuous lack of research on the genomic-level adaptive evolution of sandalwood.

RESULTS

In this study, we assembled a gap-free telomere-to-telomere (T2T) reference genome for Santalum album using PacBio HiFi, Hi-C, and ultra-long ONT data. The T2T reference genome (Sal_t2t) encompassed annotations of 24,171 genes and 25.34% repetitive sequences, in addition to all 10 centromeres and 20 telomeres across the 10 chromosomes. The results revealed that the 3 distinct parasitic species of Santalales had diverse centromeric compositions. The Copia-type long terminal repeat transposon emerged as the most significant in the S. album genome, constituting the primary sequence of the centromere and influencing gene expression. Third, in sandalwood, the presence of Copia affected the size of the centromeres and, consequently, the genome size. Identification of the sandalwood T2T genome in this study also enabled the identification of more precise organelle transfer fragments.

CONCLUSIONS

Our research provides a sandalwood T2T genome, laying the groundwork for future investigations on the evolution of energy organs in parasitic plants. Moreover, it offers novel insights into the function and evolution of centromeres, as well as the mechanisms of adaptation and parasitism.

Comparative analysis of repetitive DNA in dysploid and non-dysploid Phaseolus beans

Structural karyotype changes result from ectopic recombination events frequently associated with repetitive DNA. Although most Phaseolus species present relatively stable karyotypes with 2n = 22 chromosomes, the karyotypes of species of the Leptostachyus group show high rates of structural rearrangements, including a nested chromosome fusion that led to the dysploid chromosome number of the group (2n = 20). We examined the roles of repetitive landscapes in the rearrangements of species of the Leptostachyus group using genome-skimming data to characterize the repeatome in a range of Phaseolus species and compared them to species of that group (P. leptostachyus and P. macvaughii). LTR retrotransposons, especially the Ty3/gypsy lineage Chromovirus, were the most abundant elements in the genomes. Differences in the abundance of Tekay, Retand, and SIRE elements between P. macvaughii and P. leptostachyus were reflected in their total amounts of Ty3/gypsy and Ty1/copia. The satellite DNA fraction was the most divergent among the species, varying both in abundance and distribution, even between P. leptostachyus and P. macvaughii. The rapid turnover of repeats in the Leptostachyus group may be associated with the several rearrangements observed.

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.

Einkorn genomics sheds light on history of the oldest domesticated wheat

Einkorn (Triticum monococcum) was the first domesticated wheat species, and was central to the birth of agriculture and the Neolithic Revolution in the Fertile Crescent around 10,000 years ago1,2. Here we generate and analyse 5.2-Gb genome assemblies for wild and domesticated einkorn, including completely assembled centromeres. Einkorn centromeres are highly dynamic, showing evidence of ancient and recent centromere shifts caused by structural rearrangements. Whole-genome sequencing analysis of a diversity panel uncovered the population structure and evolutionary history of einkorn, revealing complex patterns of hybridizations and introgressions after the dispersal of domesticated einkorn from the Fertile Crescent. We also show that around 1% of the modern bread wheat (Triticum aestivum) A subgenome originates from einkorn. These resources and findings highlight the history of einkorn evolution and provide a basis to accelerate the genomics-assisted improvement of einkorn and bread wheat.

Characterization of centromeric DNA of Gossypium anomalum reveals sequence-independent enrichment dynamics of centromeric repeats

Centromeres in eukaryotes are composed of highly repetitive DNAs, which evolve rapidly and are thought to achieve a favorable structure in mature centromeres. However, how the centromeric repeat evolves into an adaptive structure is largely unknown. We characterized the centromeric sequences of Gossypium anomalum through chromatin immunoprecipitation against CENH3 antibodies. We revealed that the G. anomalum centromeres contained only retrotransposon-like repeats but were depleted in long arrays of satellites. These retrotransposon-like centromeric repeats were present in the African-Asian and Australian lineage species, suggesting that they might have arisen in the common ancestor of these diploid species. Intriguingly, we observed a substantial increase and decrease in copy numbers among African-Asian and Australian lineages, respectively, for the retrotransposon-derived centromeric repeats without apparent structure or sequence variation in cotton. This result indicates that the sequence content is not a decisive aspect of the adaptive evolution of centromeric repeats or at least retrotransposon-like centromeric repeats. In addition, two active genes with potential roles in gametogenesis or flowering were identified in CENH3 nucleosome-binding regions. Our results provide new insights into the constitution of centromeric repetitive DNA and the adaptive evolution of centromeric repeats in plants.

Assembly of the 81.6 Mb centromere of pea chromosome 6 elucidates the structure and evolution of metapolycentric chromosomes

Centromeres in the legume genera Pisum and Lathyrus exhibit unique morphological characteristics, including extended primary constrictions and multiple separate domains of centromeric chromatin. These so-called metapolycentromeres resemble an intermediate form between monocentric and holocentric types, and therefore provide a great opportunity for studying the transitions between different types of centromere organizations. However, because of the exceedingly large and highly repetitive nature of metapolycentromeres, highly contiguous assemblies needed for these studies are lacking. Here, we report on the assembly and analysis of a 177.6 Mb region of pea (Pisum sativum) chromosome 6, including the 81.6 Mb centromere region (CEN6) and adjacent chromosome arms. Genes, DNA methylation profiles, and most of the repeats were uniformly distributed within the centromere, and their densities in CEN6 and chromosome arms were similar. The exception was an accumulation of satellite DNA in CEN6, where it formed multiple arrays up to 2 Mb in length. Centromeric chromatin, characterized by the presence of the CENH3 protein, was predominantly associated with arrays of three different satellite repeats; however, five other satellites present in CEN6 lacked CENH3. The presence of CENH3 chromatin was found to determine the spatial distribution of the respective satellites during the cell cycle. Finally, oligo-FISH painting experiments, performed using probes specifically designed to label the genomic regions corresponding to CEN6 in Pisum, Lathyrus, and Vicia species, revealed that metapolycentromeres evolved via the expansion of centromeric chromatin into neighboring chromosomal regions and the accumulation of novel satellite repeats. However, in some of these species, centromere evolution also involved chromosomal translocations and centromere repositioning.

Gapless genome assembly of azalea and multi-omics investigation into divergence between two species with distinct flower color

The genus Rhododendron (Ericaceae), with more than 1000 species highly diverse in flower color, is providing distinct ornamental values and a model system for flower color studies. Here, we investigated the divergence between two parental species with different flower color widely used for azalea breeding. Gapless genome assembly was generated for the yellow-flowered azalea, Rhododendron molle. Comparative genomics found recent proliferation of long terminal repeat retrotransposons (LTR-RTs), especially Gypsy, has resulted in a 125 Mb (19%) genome size increase in species-specific regions, and a significant amount of dispersed gene duplicates (13 402) and pseudogenes (17 437). Metabolomic assessment revealed that yellow flower coloration is attributed to the dynamic changes of carotenoids/flavonols biosynthesis and chlorophyll degradation. Time-ordered gene co-expression networks (TO-GCNs) and the comparison confirmed the metabolome and uncovered the specific gene regulatory changes underpinning the distinct flower pigmentation. B3 and ERF TFs were found dominating the gene regulation of carotenoids/flavonols characterized pigmentation in R. molle, while WRKY, ERF, WD40, C2H2, and NAC TFs collectively regulated the anthocyanins characterized pigmentation in the red-flowered R simsii. This study employed a multi-omics strategy in disentangling the complex divergence between two important azaleas and provided references for further functional genetics and molecular breeding.

InpactorDB: A Plant LTR Retrotransposon Reference Library

LTR retrotransposons (LTR-RT) are major components of plant genomes. These transposable elements participate in the structure and evolution of genes and genomes through their mobility and their copy number amplification. For example, they are commonly used as evolutionary markers in genetic, genomic, and cytogenetic approaches. However, the plant research community is faced with the near absence of free availability of full-length, curated, and lineage-level classified LTR retrotransposon reference sequences. In this chapter, we will introduce InpactorDB, an LTR retrotransposon sequence database of 181 plant species representing 98 plant families for a total of 67,241 non-redundant elements. We will introduce how to use newly sequenced genomes to identify and classify LTR-RTs in a similar way with a standardized procedure using the Inpactor tool. InpactorDB is freely available at https://inpactordb.github.io .

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.

Impact of LTR-Retrotransposons on Genome Structure, Evolution, and Function in Curcurbitaceae Species

Long terminal repeat (LTR)-retrotransposons (LTR-RTs) comprise a major portion of many plant genomes and may exert a profound impact on genome structure, function, and evolution. Although many studies have focused on these elements in an individual species, their dynamics on a family level remains elusive. Here, we investigated the abundance, evolutionary dynamics, and impact on associated genes of LTR-RTs in 16 species in an economically important plant family, Cucurbitaceae. Results showed that full-length LTR-RT numbers and LTR-RT content varied greatly among different species, and they were highly correlated with genome size. Most of the full-length LTR-RTs were amplified after the speciation event, reflecting the ongoing rapid evolution of these genomes. LTR-RTs highly contributed to genome size variation via species-specific distinct proliferations. The Angela and Tekay lineages with a greater evolutionary age were amplified in Trichosanthes anguina, whereas a recent activity burst of Reina and another ancient round of Tekay activity burst were examined in Sechium edule. In addition, Tekay and Retand lineages belonging to the Gypsy superfamily underwent a recent burst in Gynostemma pentaphyllum. Detailed investigation of genes with intronic and promoter LTR-RT insertion showed diverse functions, but the term of metabolism was enriched in most species. Further gene expression analysis in G.pentaphyllum revealed that the LTR-RTs within introns suppress the corresponding gene expression, whereas the LTR-RTs within promoters exert a complex influence on the downstream gene expression, with the main function of promoting gene expression. This study provides novel insights into the organization, evolution, and function of LTR-RTs in Cucurbitaceae genomes.

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.

Diversity of the repetitive DNA fraction in Cestrum, the genus with the largest genomes within Solanaceae

BACKGROUND

Cestrum species present large genomes (2 C = ~ 24 pg), a high occurrence of B chromosomes and great diversity in heterochromatin bands. Despite this diversity, karyotypes maintain the chromosome number 2n = 16 (except when they present B chromosomes), and a relative similarity in chromosome morphology and symmetry. To deepen our knowledge of the Cestrum genome composition, low-coverage sequencing data of C. strigilatum and C. elegans were compared, including cytogenomic analyses of seven species.

METHODS AND RESULTS

Bioinformatics analyses showed retrotransposons comprising more than 70% of the repetitive fraction, followed by DNA transposons (~ 17%), but FISH assays using retrotransposon probes revealed inconspicuous and scattered signals. The four satellite DNA families here analyzed represented approximately 2.48% of the C. strigilatum dataset, and these sequences were used as probes in FISH assays. Hybridization signals were colocalized with all AT- and GC-rich sequences associated with heterochromatin, including AT-rich Cold-Sensitive Regions (CSRs). Although satellite probes hybridized in almost all tested species, a satDNA family named CsSat49 was highlighted because it predominates in centromeric regions.

CONCLUSIONS

Data suggest that the satDNA fraction is conserved in the genus, although there is variation in the number of FISH signals between karyotypes. Except to the absence of FISH signals with probes CsSat1 and CsSat72 in two species, the other satellites occurred in species of different phylogenetic clades. Some satDNA sequences have been detected in the B chromosomes, indicating that they are rich in preexisting sequences in the chromosomes of the A complement. This comparative study provides an important advance in the knowledge on genome organization and heterochromatin composition in Cestrum, especially on the distribution of satellite fractions between species and their importance for the B chromosome composition.

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

Genomes shed light on the evolution of Begonia, a mega-diverse genus

Clarifying the evolutionary processes underlying species diversification and adaptation is a key focus of evolutionary biology. Begonia (Begoniaceae) is one of the most species-rich angiosperm genera with c. 2000 species, most of which are shade-adapted. Here, we present chromosome-scale genome assemblies for four species of Begonia (B. loranthoides, B. masoniana, B. darthvaderiana and B. peltatifolia), and whole genome shotgun data for an additional 74 Begonia representatives to investigate lineage evolution and shade adaptation of the genus. The four genome assemblies range in size from 331.75 Mb (B. peltatifolia) to 799.83 Mb (B. masoniana), and harbor 22 059-23 444 protein-coding genes. Synteny analysis revealed a lineage-specific whole-genome duplication (WGD) that occurred just before the diversification of Begonia. Functional enrichment of gene families retained after WGD highlights the significance of modified carbohydrate metabolism and photosynthesis possibly linked to shade adaptation in the genus, which is further supported by expansions of gene families involved in light perception and harvesting. Phylogenomic reconstructions and genomics studies indicate that genomic introgression has also played a role in the evolution of Begonia. Overall, this study provides valuable genomic resources for Begonia and suggests potential drivers underlying the diversity and adaptive evolution of this mega-diverse clade.

Draft Sequencing Crested Wheatgrass Chromosomes Identified Evolutionary Structural Changes and Genes and Facilitated the Development of SSR Markers

Crested wheatgrass (Agropyron cristatum), a wild relative of wheat, is an attractive source of genes and alleles for their improvement. Its wider use is hampered by limited knowledge of its complex genome. In this work, individual chromosomes were purified by flow sorting, and DNA shotgun sequencing was performed. The annotation of chromosome-specific sequences characterized the DNA-repeat content and led to the identification of genic sequences. Among them, genic sequences homologous to genes conferring plant disease resistance and involved in plant tolerance to biotic and abiotic stress were identified. Genes belonging to the important groups for breeders involved in different functional categories were found. The analysis of the DNA-repeat content identified a new LTR element, Agrocen, which is enriched in centromeric regions. The colocalization of the element with the centromeric histone H3 variant CENH3 suggested its functional role in the grass centromere. Finally, 159 polymorphic simple-sequence-repeat (SSR) markers were identified, with 72 of them being chromosome- or chromosome-arm-specific, 16 mapping to more than one chromosome, and 71 mapping to all the Agropyron chromosomes. The markers were used to characterize orthologous relationships between A. cristatum and common wheat that will facilitate the introgression breeding of wheat using A. cristatum.

Transposable Element Populations Shed Light on the Evolutionary History of Wheat and the Complex Co-Evolution of Autonomous and Non-Autonomous Retrotransposons

Wheat has one of the largest and most repetitive genomes among major crop plants, containing over 85% transposable elements (TEs). TEs populate genomes much in the way that individuals populate ecosystems, diversifying into different lineages, sub-families and sub-populations. The recent availability of high-quality, chromosome-scale genome sequences from ten wheat lines enables a detailed analysis how TEs evolved in allohexaploid wheat, its diploids progenitors, and in various chromosomal haplotype segments. LTR retrotransposon families evolved into distinct sub-populations and sub-families that were active in waves lasting several hundred thousand years. Furthermore, It is shown that different retrotransposon sub-families were active in the three wheat sub-genomes, making them useful markers to study and date polyploidization events and chromosomal rearrangements. Additionally, haplotype-specific TE sub-families are used to characterize chromosomal introgressions in different wheat lines. Additionally, populations of non-autonomous TEs co-evolved over millions of years with their autonomous partners, leading to complex systems with multiple types of autonomous, semi-autonomous and non-autonomous elements. Phylogenetic and TE population analyses revealed the relationships between non-autonomous elements and their mobilizing autonomous partners. TE population analysis provided insights into genome evolution of allohexaploid wheat and genetic diversity of species, and may have implication for future crop breeding.

Separable roles for RNAi in regulation of transposable elements and viability in the fission yeast Schizosaccharomyces japonicus

RNA interference (RNAi) is a conserved mechanism of small RNA-mediated genome regulation commonly involved in suppression of transposable elements (TEs) through both post-transcriptional silencing, and transcriptional repression via heterochromatin assembly. The fission yeast Schizosaccharomyces pombe has been extensively utilised as a model for studying RNAi pathways. However, this species is somewhat atypical in that TEs are not major targets of RNAi, and instead small RNAs correspond primarily to non-coding pericentromeric repeat sequences, reflecting a specialised role for the pathway in promoting heterochromatin assembly in these regions. In contrast, in the related fission yeast Schizosaccharomyces japonicus, sequenced small RNAs correspond primarily to TEs. This suggests there may be fundamental differences in the operation of RNAi pathways in these two related species. To investigate these differences, we probed RNAi function in S. japonicus. Unexpectedly, and in contrast to S. pombe, we found that RNAi is essential in this species. Moreover, viability of RNAi mutants can be rescued by mutations implicated in enhancing RNAi-independent heterochromatin propagation. These rescued strains retain heterochromatic marks on TE sequences, but exhibit derepression of TEs at the post-transcriptional level. Our findings indicate that S. japonicus retains the ancestral role of RNAi in facilitating suppression of TEs via both post-transcriptional silencing and heterochromatin assembly, with specifically the heterochromatin pathway being essential for viability, likely due to a function in genome maintenance. The specialised role of RNAi in heterochromatin assembly in S. pombe appears to be a derived state that emerged after the divergence of S. japonicus.

The chromosomes of many species are populated by repetitive transposable elements that are able to “jump” throughout the genome. The consequences of these mobilisations can be catastrophic, resulting in disruption of genes or chromosomal rearrangements, thus organisms usually employ defence mechanisms to keep these elements inactivated. The most widespread of these systems is RNA interference, which utilises small RNA molecules to direct either packaging of transposable element DNA into repressive heterochromatin, or degradation of RNA transcripts. Many fundamental discoveries about RNAi function have been made in the model fission yeast Schizosaccharomyces pombe; however, this species is unusual as it does not generally employ RNAi to control its transposable elements. We found that in a lesser studied relative, Schizosaccharomyces japonicus, small RNAs are required to silence transposable elements, and that this silencing occurs via both formation of heterochromatin and degradation of transcripts. This dual function RNAi pathway targeting transposable elements that appear to cluster at centromeres is very similar to systems seen in complex multicellular organisms, thus our findings reveal S. japonicus to be an exciting emergent model in which to study RNAi and centromere function.

The Sister Chromatid Division of the Heteromorphic Sex Chromosomes in Silene Species and Their Transmissibility towards the Mitosis

Young sex chromosomes possess unique and ongoing dynamics that allow us to understand processes that have an impact on their evolution and divergence. The genus Silene includes species with evolutionarily young sex chromosomes, and two species of section Melandrium, namely Silene latifolia (24, XY) and Silene dioica (24, XY), are well-established models of sex chromosome evolution, Y chromosome degeneration, and sex determination. In both species, the X and Y chromosomes are strongly heteromorphic and differ in the genomic composition compared to the autosomes. It is generally accepted that for proper cell division, the longest chromosomal arm must not exceed half of the average length of the spindle axis at telophase. Yet, it is not clear what are the dynamics between males and females during mitosis and how the cell compensates for the presence of the large Y chromosome in one sex. Using hydroxyurea cell synchronization and 2D/3D microscopy, we determined the position of the sex chromosomes during the mitotic cell cycle and determined the upper limit for the expansion of sex chromosome non-recombining region. Using 3D specimen preparations, we found that the velocity of the large chromosomes is compensated by the distant positioning from the central interpolar axis, confirming previous mathematical modulations.

Advances in the Molecular Cytogenetics of Bananas, Family Musaceae

The banana is a staple food crop and represents an important trade commodity for millions of people living in tropical and subtropical countries. The most important edible banana clones originated from natural crosses between diploid Musa balbisiana and various subspecies of M. acuminata. It is worth mentioning that evolution and speciation in the Musaceae family were accompanied by large-scale chromosome structural changes, indicating possible reasons for lower fertility or complete sterility of these vegetatively propagated clones. Chromosomal changes, often accompanied by changes in genome size, are one of the driving forces underlying speciation in plants. They can clarify the genomic constitution of edible bananas and shed light on their origin and on diversification processes in members of the Musaceae family. This article reviews the development of molecular cytogenetic approaches, ranging from classical fluorescence in situ hybridization (FISH) using common cytogenetic markers to oligo painting FISH. We discuss differences in genome size and chromosome number across the Musaceae family in addition to the development of new chromosome-specific cytogenetic probes and their use in genome structure and comparative karyotype analysis. The impact of these methodological advances on our knowledge of Musa genome evolution at the chromosomal level is demonstrated. In addition to citing published results, we include our own new unpublished results and outline future applications of molecular cytogenetics in banana research.

The Highly Repeat-Diverse (Peri) Centromeres of White Lupin (Lupinus albus L.)

Plant genomes are known to be mainly composed of repetitive DNA sequences. Regardless of the non-genic function of these sequences, they are important for chromosome structure and stability during cell-cycle. Based on the recent available whole-genome assembly of white lupin (Lupinus albus L.; WL), we have in silico annotated and in situ mapped the main classes of DNA repeats identified with RepeatExplorer. A highly diverse and an abundance of satellite DNAs were found representing more than 10 families, where three of them were highly associated with CENH3-immunoprecipitated chromatin. Applying a strategy of several re-hybridization steps with different combinations of satDNA, rDNA, and LTR-RTs probes, we were able to construct a repeat-based chromosome map for the identification of most chromosome pairs. Two families of LTR retrotransposons, Ty1/copia SIRE and Ty3/gypsy Tekay, were highly abundant at pericentromeric regions, while the centromeric retrotransposon of WL (CRWL) from the CRM clade showed strong centromere-specific localization in most chromosomes and was also highly enriched with CENH3-immunoprecipitated chromatin. FISH mapping of repeat DNA showed some incongruences with the reference genome, which can be further used for improving the current version of the genome. Our results demonstrate that despite the relatively small genome of WL, a high diversity of pericentromeric repeats was found, emphasizing the rapid evolution of repeat sequences in plant genomes.

Centromere-Specific Retrotransposons and Very-Long-Chain Fatty Acid Biosynthesis in the Genome of Yellowhorn (Xanthoceras sorbifolium, Sapindaceae), an Oil-Producing Tree With Significant Drought Resistance

In-depth genome characterization is still lacking for most of biofuel crops, especially for centromeres, which play a fundamental role during nuclear division and in the maintenance of genome stability. This study applied long-read sequencing technologies to assemble a highly contiguous genome for yellowhorn (Xanthoceras sorbifolium), an oil-producing tree, and conducted extensive comparative analyses to understand centromere structure and evolution, and fatty acid biosynthesis. We produced a reference-level genome of yellowhorn, ∼470 Mb in length with ∼95% of contigs anchored onto 15 chromosomes. Genome annotation identified 22,049 protein-coding genes and 65.7% of the genome sequence as repetitive elements. Long terminal repeat retrotransposons (LTR-RTs) account for ∼30% of the yellowhorn genome, which is maintained by a moderate birth rate and a low removal rate. We identified the centromeric regions on each chromosome and found enrichment of centromere-specific retrotransposons of LINE1 and Gypsy in these regions, which have evolved recently (∼0.7 MYA). We compared the genomes of three cultivars and found frequent inversions. We analyzed the transcriptomes from different tissues and identified the candidate genes involved in very-long-chain fatty acid biosynthesis and their expression profiles. Collinear block analysis showed that yellowhorn shared the gamma (γ) hexaploidy event with Vitis vinifera but did not undergo any further whole-genome duplication. This study provides excellent genomic resources for understanding centromere structure and evolution and for functional studies in this important oil-producing plant.

An 82 bp tandem repeat family typical of 3' non-coding end of Gypsy/TAT LTR retrotransposons is conserved in Coffea spp. pericentromeres

Coffea spp. chromosomes are very small and accumulate a variety of repetitive DNA families around the centromeres. However, the proximal regions of Coffea chromosomes remain poorly understood, especially regarding the nature and organisation of the sequences. Taking advantage of the genome sequences of C. arabica (2n = 44), C. canephora, and C. eugenioides (C. arabica progenitors with 2n = 22) and good coverage genome sequencing of dozens of other wild Coffea spp., repetitive DNA sequences were identified, and the genomes were compared to decipher particularities of pericentromeric structures. The searches revealed a short tandem repeat (82 bp length) typical of Gypsy/TAT LTR retrotransposons, named Coffea_sat11. This repeat organises clusters with fragments of other transposable elements, comprising regions of non-coding RNA production. Cytogenomic analyses showed that Coffea_sat11 extends from the pericentromeres towards the middle of the chromosomal arms. This arrangement was observed in the allotetraploid C. arabica chromosomes, as well as in its progenitors. This study improves our understanding of the role of the Gypsy/TAT LTR retrotransposon lineage in the organisation of Coffea pericentromeres, as well as the conservation of Coffea_sat11 within the genus. The relationships between fragments of other transposable elements and the functional aspects of these sequences on the pericentromere chromatin were also evaluated. Highlights: A scattered short tandem repeat, typical of Gypsy/TAT LTR retrotransposons, associated with several fragments of other transposable elements, accumulates in the pericentromeres of Coffea chromosomes. This arrangement is preserved in all clades of the genus and appears to have a strong regulatory role in the organisation of chromatin around centromeres.

The genomic ecosystem of transposable elements in maize

Transposable elements (TEs) constitute the majority of flowering plant DNA, reflecting their tremendous success in subverting, avoiding, and surviving the defenses of their host genomes to ensure their selfish replication. More than 85% of the sequence of the maize genome can be ascribed to past transposition, providing a major contribution to the structure of the genome. Evidence from individual loci has informed our understanding of how transposition has shaped the genome, and a number of individual TE insertions have been causally linked to dramatic phenotypic changes. Genome-wide analyses in maize and other taxa have frequently represented TEs as a relatively homogeneous class of fragmentary relics of past transposition, obscuring their evolutionary history and interaction with their host genome. Using an updated annotation of structurally intact TEs in the maize reference genome, we investigate the family-level dynamics of TEs in maize. Integrating a variety of data, from descriptors of individual TEs like coding capacity, expression, and methylation, as well as similar features of the sequence they inserted into, we model the relationship between attributes of the genomic environment and the survival of TE copies and families. In contrast to the wholesale relegation of all TEs to a single category of junk DNA, these differences reveal a diversity of survival strategies of TE families. Together these generate a rich ecology of the genome, with each TE family representing the evolution of a distinct ecological niche. We conclude that while the impact of transposition is highly family- and context-dependent, a family-level understanding of the ecology of TEs in the genome can refine our ability to predict the role of TEs in generating genetic and phenotypic diversity.

The Genomics of Plant Satellite DNA

The twenty-first century began with a certain indifference to the research of satellite DNA (satDNA). Neither genome sequencing projects were able to accurately encompass the study of satDNA nor classic methodologies were able to go further in undertaking a better comprehensive study of the whole set of satDNA sequences of a genome. Nonetheless, knowledge of satDNA has progressively advanced during this century with the advent of new analytical techniques. The enormous advantages that genome-wide approaches have brought to its analysis have now stimulated a renewed interest in the study of satDNA. At this point, we can look back and try to assess more accurately many of the key questions that were left unsolved in the past about this enigmatic and important component of the genome. I review here the understanding gathered on plant satDNAs over the last few decades with an eye on the near future.

The nature of intraspecific and interspecific genome size variation in taxonomically complex eyebrights

BACKGROUND AND AIMS

Genome size varies considerably across the diversity of plant life. Although genome size is, by definition, affected by genetic presence/absence variants, which are ubiquitous in population sequencing studies, genome size is often treated as an intrinsic property of a species. Here, we studied intra- and interspecific genome size variation in taxonomically complex British eyebrights (Euphrasia, Orobanchaceae). Our aim is to document genome size diversity and investigate underlying evolutionary processes shaping variation between individuals, populations and species.

METHODS

We generated genome size data for 192 individuals of diploid and tetraploid Euphrasia and analysed genome size variation in relation to ploidy, taxonomy, population affiliation and geography. We further compared the genomic repeat content of 30 samples.

KEY RESULTS

We found considerable intraspecific genome size variation, and observed isolation-by-distance for genome size in outcrossing diploids. Tetraploid Euphrasia showed contrasting patterns, with genome size increasing with latitude in outcrossing Euphrasia arctica, but with little genome size variation in the highly selfing Euphrasia micrantha. Interspecific differences in genome size and the genomic proportions of repeat sequences were small.

CONCLUSIONS

We show the utility of treating genome size as the outcome of polygenic variation. Like other types of genetic variation, such as single nucleotide polymorphisms, genome size variation may be affected by ongoing hybridization and the extent of population subdivision. In addition to selection on associated traits, genome size is predicted to be affected indirectly by selection due to pleiotropy of the underlying presence/absence variants.

Telomere-to-telomere gapless chromosomes of banana using nanopore sequencing

Long-read technologies hold the promise to obtain more complete genome assemblies and to make them easier. Coupled with long-range technologies, they can reveal the architecture of complex regions, like centromeres or rDNA clusters. These technologies also make it possible to know the complete organization of chromosomes, which remained complicated before even when using genetic maps. However, generating a gapless and telomere-to-telomere assembly is still not trivial, and requires a combination of several technologies and the choice of suitable software. Here, we report a chromosome-scale assembly of a banana genome (Musa acuminata) generated using Oxford Nanopore long-reads. We generated a genome coverage of 177X from a single PromethION flowcell with near 17X with reads longer than 75 kbp. From the 11 chromosomes, 5 were entirely reconstructed in a single contig from telomere to telomere, revealing for the first time the content of complex regions like centromeres or clusters of paralogous genes.

Repetitive genomic regions and the inference of demographic history

Inference of demographic histories using whole-genome datasets has provided insights into diversification, adaptation, hybridization, and plant-pathogen interactions, and stimulated debate on the impact of anthropogenic interventions and past climate on species demography. However, the impact of repetitive genomic regions on these inferences has mostly been ignored by masking of repeats. We use the Populus trichocarpa genome (Pop_tri_v3) to show that masking of repeat regions leads to lower estimates of effective population size (Ne) in the distant past in contrast to an increase in Ne estimates in recent times. However, in human datasets, masking of repeats resulted in lower estimates of Ne at all time points. We demonstrate that repeats affect demographic inferences using diverse methods like PSMC, MSMC, SMC++, and the Stairway plot. Our genomic analysis revealed that the biases in Ne estimates were dependent on the repeat class type and its abundance in each atomic interval. Notably, we observed a weak, yet consistently significant negative correlation between the repeat abundance of an atomic interval and the Ne estimates for that interval, which potentially reflects the recombination rate variation within the genome. The rationale for the masking of repeats has been that variants identified within these regions are erroneous. We find that polymorphisms in some repeat classes occur in callable regions and reflect reliable coalescence histories (e.g., LTR Gypsy, LTR Copia). The current demography inference methods do not handle repeats explicitly, and hence the effect of individual repeat classes needs careful consideration in comparative analysis. Deciphering the repeat demographic histories might provide a clear understanding of the processes involved in repeat accumulation.

Full-length LTR retroelements in Capsicum annuum revealed a few species-specific family bursts with insertional preferences

Capsicum annuum is a species that has undergone an expansion of the size of its genome caused mainly by the amplification of repetitive DNA sequences, including mobile genetic elements. Based on information obtained from sequencing the genome of pepper, the estimated fraction of retroelements is approximately 81%, and previous results revealed an important contribution of lineages derived from Gypsy superfamily. However, the dynamics of the retroelements in the C. annuum genome is poorly understood. In this way, the present work seeks to investigate the phylogenetic diversity and genomic abundance of the families of autonomous (complete and intact) LTR retroelements from C. annuum and inspect their distribution along its chromosomes. In total, we identified 1151 structurally full-length retroelements (340 Copia; 811 Gypsy) grouped in 124 phylogenetic families in the base of their retrotranscriptase. All the evolutive lineages of LTR retroelements identified in plants were present in pepper; however, three of them comprise 83% of the entire LTR retroelements population, the lineages Athila, Del/Tekay, and Ale/Retrofit. From them, only three families represent 70.8% of the total number of the identified retroelements. A massive family-specific wave of amplification of two of them occurred in the last 0.5 Mya (GypsyCa_16; CopiaCa_01), whereas the third is more ancient and occurred 3.0 Mya (GypsyCa_13). Fluorescent in situ hybridization performed with family and lineage-specific probes revealed contrasting patterns of chromosomal affinity. Our results provide a database of the populations LTR retroelements specific to C. annuum genome. The most abundant families were analyzed according to chromosome insertional preferences, suppling useful tools to the design of retroelement-based markers specific to the species.

Sequence, Chromatin and Evolution of Satellite DNA

Satellite DNA consists of abundant tandem repeats that play important roles in cellular processes, including chromosome segregation, genome organization and chromosome end protection. Most satellite DNA repeat units are either of nucleosomal length or 5–10 bp long and occupy centromeric, pericentromeric or telomeric regions. Due to high repetitiveness, satellite DNA sequences have largely been absent from genome assemblies. Although few conserved satellite-specific sequence motifs have been identified, DNA curvature, dyad symmetries and inverted repeats are features of various satellite DNAs in several organisms. Satellite DNA sequences are either embedded in highly compact gene-poor heterochromatin or specialized chromatin that is distinct from euchromatin. Nevertheless, some satellite DNAs are transcribed into non-coding RNAs that may play important roles in satellite DNA function. Intriguingly, satellite DNAs are among the most rapidly evolving genomic elements, such that a large fraction is species-specific in most organisms. Here we describe the different classes of satellite DNA sequences, their satellite-specific chromatin features, and how these features may contribute to satellite DNA biology and evolution. We also discuss how the evolution of functional satellite DNA classes may contribute to speciation in plants and animals.

Large vs small genomes in Passiflora: the influence of the mobilome and the satellitome

While two lineages of retrotransposons were more abundant in larger Passiflora genomes, the satellitome was more diverse and abundant in the smallest genome analysed. Repetitive sequences are ubiquitous and fast-evolving elements responsible for size variation and large-scale organization of plant genomes. Within Passiflora genus, a tenfold variation in genome size, not attributed to polyploidy, is known. Here, we applied a combined in silico and cytological approach to study the organization and diversification of repetitive elements in three species of this genus representing its known range in genome size variation. Sequences were classified in terms of type and repetitiveness and the most abundant were mapped to chromosomes. We identified long terminal repeat (LTR) retrotransposons as the most abundant elements in the three genomes, showing a considerable variation among species. Satellite DNAs (satDNAs) were less representative, but highly diverse between subgenera. Our results clearly confirm that the largest genome species (Passiflora quadrangularis) presents a higher accumulation of repetitive DNA sequences, specially Angela and Tekay elements, making up most of its genome. Passiflora cincinnata, with intermediate genome and from the same subgenus, showed similarity with P. quadrangularis regarding the families of repetitive DNA sequences, but in different proportions. On the other hand, Passiflora organensis, the smallest genome, from a different subgenus, presented greater diversity and the highest proportion of satDNA. Altogether, our data indicates that while large genomes evolved by an accumulation of retrotransposons, the smallest genome known for the genus has evolved by diversification of different repeat types, particularly satDNAs.