The contributions of transposable elements to the structure, function, and evolution of plant genomes

Repeat competition and ecological shifts drive the evolution of the mobilome in Rhynchospora Vahl (Cyperaceae), the holocentric beaksedges

BACKGROUND AND AIMS

Genomic changes triggered by polyploidy, chromosomal rearrangements and/ or environmental stress are among factors that affect the activity of mobile elements, particularly long terminal repeat retrotransposons (LTR-RTs) and DNA transposons. Because these elements can proliferate and move throughout host genomes, altering the genetic, epigenetic and nucleotypic landscape, they have been recognized as a relevant evolutionary force. Beaksedges (Rhynchospora) stand out for their wide cosmopolitan distribution, high diversity (~400 species) and holocentric chromosomes related to high karyotypic diversity and a centromere-specific satDNA, Tyba. This makes the genus an interesting model to investigate the interactions between repetitive elements, phylogenetic relationships and ecological variables.

METHODS

Here we used comparative phylogenetic methods to investigate the forces driving the evolution of the entire set of mobile elements (mobilome) in the holocentric genus Rhynchospora. We statistically tested the impact of phylogenetic relationships, abundance of holocentromeric satDNA Tyba, diversity of repeatome composition, ecological variables and chromosome number in mobile element diversification.

KEY RESULTS

Tyba abundance was found to be inversely correlated with LTR-RT content. Decrease of LTR abundance and diversity was also related to increase in chromosome number (likely due to fission events) and colonization of dry environments in the northern hemisphere. In contrast, we found constant LTR insertions throughout time in species with lower chromosome numbers in rainier environments in South America. A multivariate model showed that different traits drive LTR abundance, especially repeat diversity and Tyba abundance. Other mobile elements, such as non-LTR RTs and DNA transposons, had insufficient abundance to be included in our models.

CONCLUSIONS

Our findings suggest that LTR evolution is strongly impacted by the holocentric characteristics of Rhynchospora chromosomes, correlating with species diversification and biome shifts, and supporting a holokinetic drive model of evolution and a competitive scenario with Tyba. Altogether, our results present evidence of multi-trait influence on LTR-RT dynamics and provide a broader understanding of transposable element evolution in a macroevolutionary context.

De novo, high-quality assembly and annotation of the halophyte grass Aeluropus littoralis draft genome and identification of A20/AN1 zinc finger protein family

BACKGROUND

Aeluropus littoralis is considered a valuable natural forage plant for ruminant livestock and is highly tolerant to extreme abiotic stresses, especially salinity, drought, and heat. It is a monocotyledonous halophyte, has salt glands, performs C4-type photosynthesis and has a close genetic relationship with cereal crops. Moreover, previous studies have shown its huge potential as a reservoir of genes and promoters to understand and improve abiotic stress tolerance in crops.

RESULTS

The sequencing and hybrid assembly of the A. littoralis genome (2n = 2X = 20) using short and long reads from the BGISeq-500 and PacBio high-fidelity (HiFi) sequencing platforms, respectively. Using the k-mer analysis method, the haploid genome size of A. littoralis was estimated to be 360 Mb (with a heterozygosity rate of 1.88%). The hybrid assembled genome included 4,078 contigs with a GC content of 44% and covered 348 Mb. The longest contig and the N50 values were 5.1 Mb and 133.77 kb, respectively. The Benchmarking Universal Single Copy Ortholog (BUSCO) value was 91.1%, indicating good integrity of the assembled genome. The discovered repetitive elements accounted for 90.6 Mb, representing 26.03% of the total genome, and included a significant component of transposable elements (11.48%, ~40 Mb). Using a homology-based approach, 35,147 genes were predicted from the genome assembly. We next focused our analysis on the zinc-finger A20/AN1 gene family, a member of which (AlSAP) was previously shown to confer increased tolerance to osmotic and salt stresses when it was over-expressed in tobacco, wheat, and rice. Here, we identified the complete set of members of this family in the Aeluropus littoralis genome, thereby laying the foundation for their future functional analysis in cereal crops. In addition, the expression patterns of four novel genes from this family were analyzed by qPCR.

CONCLUSION

This resource and our findings will contribute to improve the current understanding of salinity tolerance in halophytes while providing useful genes and allelic variation to improve salinity and drought tolerance in cereals through genetic engineering and gene editing, respectively.

Construction of a comprehensive library of repeated sequences for the annotation of Citrus genomes

BACKGROUND

The comprehensive annotation of repeated sequences in genomes is an essential prerequisite for studying the dynamics of these sequences over time and their involvement in gene regulation. Currently, the diversity of repeated sequences in Citrus genomes is only partially characterized because the annotations have been performed using heterogeneous bioinformatics tools, each with its specificity and dedicated only to the annotation of transposable elements.

RESULTS

We combined complementary repeat-finding programs including REPET, CAULIFINDER, and TAREAN, to enable the identification of all types of repetitive sequences found in plant genomes, including transposable elements, endogenous caulimovirids, and satellite DNAs. A fine-grained annotation method was first developed to create a consensus sequence library of repeated sequences identified in the genome assemblies of C. medica, C. micrantha, C. reticulata, and C. maxima, the four ancestral parental species involved in the formation of economically valuable cultivated Citrus varieties. A second, faster annotation method was developed to enrich the dataset by adding new repeated sequences retrieved from genome assemblies of other Citrus species and closely related species belonging to the Aurantioideae subfamily. The final reference library contains 3,091 consensus sequences, of which 94.5% are transposable elements. The diversity of endogenous caulimovirids was characterized for the first time within the genus Citrus, contributing 160 consensus sequences to the final reference library. Finally, 10 satellite DNAs were also identified.

CONCLUSION

Combining multiple repeat detection methods enables the comprehensive annotation of all repeated sequences in Citrus genomes. Using the final reference library reported in this work will improve our understanding of the dynamics of repeated sequences during Citrus speciation, particularly following the genome duplication and hybridization events that led to modern cultivars. The exploration of repeat position insertions along chromosomes using the developed web interface, RepeatLoc Citrus, will also make it possible to further investigate the role of transposable elements and endogenous caulimovirids in genome structure and gene regulation in Citrus species.

Phenotypic and genomic signatures across wild Rosa species open new horizons for modern rose breeding

The cultivation and domestication of roses reflects cultural exchanges and shifts in aesthetics that have resulted in today's most popular ornamental plant group. However, the narrow genetic foundation of cultivated roses limits their further improvement. Wild Rosa species harbour vast genetic diversity, yet their utilization is impeded by taxonomic confusion. Here we generated a phased and gap-free reference genome of Rosa persica for phylogenetic and population genomic analyses of a large collection of Rosa samples. The robust nuclear and plastid phylogenies support most of the morphology-based traditional taxonomy of Rosa. Population genomic analyses disclosed potential genetic exchanges among sections, indicating the northwest and southwest of China as two independent centres of diversity for Rosa. Analyses of domestication traits provide insights into selection processes related to flower colour, fragrance, double flower and resistance. This study provides a comprehensive understanding of rose domestication and lays a solid foundation for future re-domestication and innovative breeding efforts using wild resources.

Widespread impact of transposable elements on the evolution of post-transcriptional regulation in the cotton genus Gossypium

BACKGROUND

Transposable element (TE) expansion has long been known to mediate genome evolution and phenotypic diversity in organisms, but its impact on the evolution of post-transcriptional regulation following species divergence remains unclear.

RESULTS

To address this issue, we perform long-read direct RNA sequencing, polysome profiling sequencing, and small RNA sequencing in the cotton genus Gossypium, the species of which range more than three folds in genome size. We find that TE expansion contributes to the turnover of transcription splicing sites and regulatory sequences, leading to changes in alternative splicing patterns and the expression levels of orthologous genes. We also find that TE-derived upstream open reading frames and microRNAs serve as regulatory elements mediating differences in the translation levels of orthologous genes. We further identify genes that exhibit lineage-specific divergence at the transcriptional, splicing, and translational levels, and showcase the high flexibility of gene expression regulation in the evolutionary process.

CONCLUSIONS

Our work highlights the significant role of TE in driving post-transcriptional regulation divergence in the cotton genus. It offers insights for deciphering the evolutionary mechanisms of cotton species and the formation of biological diversity.

Tandem LTR-retrotransposon structures are common and highly polymorphic in plant genomes

BACKGROUND

LTR-retrotransposons (LTR-RT) are a major component of plant genomes and important drivers of genome evolution. Most LTR-RT copies in plant genomes are defective elements found as truncated copies, nested insertions or as part of more complex structures. The recent availability of highly contiguous plant genome assemblies based on long-read sequences now allows to perform detailed characterization of these complex structures and to evaluate their importance for plant genome evolution.

RESULTS

The detailed analysis of two rice loci containing complex LTR-RT structures showed that they consist of tandem arrays of LTR copies sharing internal LTRs. Our analyses suggests that these LTR-RT tandems are the result of a single insertion and not of the recombination of two independent LTR-RT elements. Our results also suggest that gypsy elements may be more prone to form these structures. We show that these structures are highly polymorphic in rice and therefore have the potential to generate genetic variability. We have developed a computational pipeline (IDENTAM) that scans genome sequences and identifies tandem LTR-RT candidates. Using this tool, we have detected 266 tandems in a pangenome built from the genomes of 76 accessions of cultivated and wild rice, showing that tandem LTR-RT structures are frequent and highly polymorphic in rice. Running IDENTAM in the Arabidopsis, almond and cotton genomes showed that LTR-RT tandems are frequent in plant genomes of different size, complexity and ploidy level. The complexity of differentiating intra-element variations at the nucleotide level among haplotypes is very high, and we found that graph-based pangenomic methodologies are appropriate to resolve these structures.

CONCLUSIONS

Our results show that LTR-RT elements can form tandem arrays. These structures are relatively abundant and highly polymorphic in rice and are widespread in the plant kingdom. Future studies will contribute to understanding how these structures originate and whether the variability that they generate has a functional impact.

The calcium-dependent protein kinase CmaCPK4 regulates sex determination in pumpkin (Cucurbita maxima D.)

Pumpkin (Cucurbita maxima D.) is typically monoecious with individual male and female flowers, and its yield is associated with the degree of femaleness, i.e. the ratio of female to male flowers produced by the plant. Subgynoecy represents a sex form with a high degree of femaleness, but the regulatory mechanisms in pumpkin remain poorly understood. In this study, using the F2 population crossed from the subgynoecious line 2013-12 and the monoecious line 9-6, we initially identified a recessive locus to control the subgynoecious trait and named it sg1. After bulked segregant analysis with whole-genome resequencing and molecular marker linkage analysis, the sg1 locus was mapped to pumpkin Chromosome 2. Genetic sequence analysis found a pumpkin calcium-dependent protein kinase (CPK) gene, CmaCPK4, in the mapping interval as the candidate gene. A retrotransposon insertion identified within the promoter elevated CmaCPK4 expression in 2013-12. Morphological characterization of near-isogenic lines containing the sg1 allele showed increases in the ratio of female flowers and high ethylene contents in terminal buds compared with the receptor parent. Heterologous overexpression of CmaCPK4 significantly increased the ratio of female flowers in cucumber (Cucumis sativus). Furthermore, CmaCPK4 directly interacts with and phosphorylates 1-aminocyclopropane-1-carboxylate synthase 5 (CmaACS5) and 1-aminocyclopropane-1-carboxylate synthase 7 (CmaACS7), resulting in increased ethylene content in 2013-12, which affected pumpkin sex determination. These findings provide insights into the role of the CmaCPK4-CmaACS5/CmaACS7 module in ethylene-induced sex determination in pumpkin.

Genome-wide identification and expression analysis of phytochrome gene family in Aikang58 wheat (Triticum aestivum L.)

Phytochromes are essential photoreceptors in plants that sense red and far-red light, playing a vital role in regulating plant growth and development through light signal transduction. Despite extensive research on phytochromes in model plants like Arabidopsis and rice, they have received relatively little attention in wheat. In this study, we employed bioinformatics methods to identify eight TaAkPHY genes in the Aikang58 wheat variety. Based on gene structure, conserved domains, and phylogenetic relationships, the TaAkPHY gene family exhibits a high degree of conservation. Synteny analysis revealed the evolutionary history of the PHY genes in Aikang58 and Chinese Spring wheat (Triticum aestivum L.), barley (Hordeum vulgare L.), rice (Oryza sativa L.), maize (Zea mays L.), quinoa (Chenopodium quinoa Willd.), soybean [Glycine max (L.) Merr.], and Arabidopsis [Arabidopsis thaliana (L.) Heynh.]. Among these species, wheat is most closely related to barley, followed by rice and maize. The cis-acting element analysis indicates that the promoter regions of TaAkPHY genes contain a large number of CAT-box, CGTCA-motif, GC-motif, etc., which are mainly involved in plant development, hormone response, and stress response. Gene expression profiling demonstrated that TaAkPHY genes exhibit varying expression levels across different tissues and are induced by various stress conditions and plant hormone treatments. Co-expression network analysis suggested that TaAkPHY genes may specifically regulate downstream genes associated with stress responses, chloroplast development, and circadian rhythms. Additionally, the least absolute shrinkage and selection operator (LASSO) regression algorithm in machine learning was used to screen transcription factors such as bHLH, WRKY, and MYB that influenced the expression of TaAkPHY genes. This method helps to quickly extract key influencing factors from a large amount of complex data. Overall, these findings provide new insights into the role of phytochromes in wheat growth, development, and stress responses, laying a foundation for future research on phytochromes in wheat.

Miniature-inverted-repeat transposable elements contribute to phenotypic variation regulation of rice induced by space environment

Introduction

Rice samples exposed to the space environment have generated diverse phenotypic variations. Miniature-inverted-repeat transposable elements (MITEs), often found adjacent to genes, play a significant role in regulating the plant genome. Herein, the contribution of MITEs in regulating space-mutagenic phenotypes was explored.

Methods

The space-mutagenic phenotype changes in the F3 to F5 generations of three space-mutagenic lines from the rice varieties Dongnong423 (DN423) and Dongnong (DN416) were meticulously traced. Rice leaves samples at the heading stage from three space-mutagenic lines were subjected to high coverage whole-genome bisulfite sequencing and whole-genome sequencing. These analyses were conducted to investigate the effects of MITEs related epigenetic and genetic variations on space-mutagenic phenotypes.

Results and discussion

Studies have indicated that MITEs within gene regulatory regions might contribute to the formation and differentiation of space-mutagenic phenotypes. The space environment has been shown to induce the transposable elements insertion polymorphisms of MITEs (MITEs-TIPs), with a notable preference for insertion near genes involved in stress response and phenotype regulation. The space-induced MITEs-TIPs contributed to the formation of space-mutagenic phenotype by modulating the expression of gene near the insertion site. This study underscored the pivotal role of MITEs in modulating plant phenotypic variation induced by the space environment, as well as the transgenerational stability of these phenotypic variants.

The smallest angiosperm genomes may be the price for effective traps of bladderworts

BACKGROUND

Species of the carnivorous family Lentibulariaceae exhibit the smallest genomes in flowering plants. We explored the hypothesis that their minute genomes result from the unique mitochondrial cytochrome c oxidase (COX) mutation. The mutation may boost mitochondrial efficiency, which is especially useful for suction-bladder traps of Utricularia, but also increase DNA-damaging reactive oxygen species, leading to genome shrinkage through deletion-biased DNA repair. We aimed to explore the impact of this mutation on genome size, providing insights into genetic mutation roles in plant genome evolution under environmental pressures.

METHODS

We compiled and measured genome and mean chromosome sizes for 127 and 67 species, respectively, representing all three genera (Genlisea, Pinguicula and Utricularia) of Lentibulariaceae. We also isolated and analysed COX sequences to detect the mutation. Through phylogenetic regressions and Ornstein-Uhlenbeck models of trait evolution, we assessed the impact of the COX mutation on the genome and chromosome sizes across the family.

RESULTS

Our findings reveal significant correlations between the COX mutation and smaller genome and chromosome sizes. Specifically, species carrying the ancestral COX sequence exhibited larger genomes and chromosomes than those with the novel mutation. This evidence supports the notion that the COX mutation contributes to genome downsizing, with statistical analyses confirming a directional evolution towards smaller genomes in species harbouring these mutations.

CONCLUSIONS

Our study confirms that the COX mutation in Lentibulariaceae is associated with genome downsizing, probably driven by increased reactive oxygen species production and subsequent DNA damage requiring deletion-biased repair mechanisms. While boosting mitochondrial energy output, this genetic mutation compromises genome integrity and may potentially affect recombination rates, illustrating a complex trade-off between evolutionary advantages and disadvantages. Our results highlight the intricate processes by which genetic mutations and environmental pressures shape genome size evolution in carnivorous plants.

Centromere drive may propel the evolution of chromosome and genome size in plants

BACKGROUND

Genome size is influenced by natural selection and genetic drift acting on variations from polyploidy and repetitive DNA sequences. We hypothesized that centromere drive, where centromeres compete for inclusion in the functional gamete during meiosis, may also affect genome and chromosome size. This competition occurs in asymmetric meiosis, where only one of the four meiotic products becomes a gamete. If centromere drive influences chromosome size evolution, it may also impact post-polyploid diploidization, where a polyploid genome is restructured to function more like a diploid through chromosomal rearrangements, including fusions. We tested if plant lineages with asymmetric meiosis exhibit faster chromosome size evolution compared to those with only symmetric meiosis, which lack centromere drive as all four meiotic products become gametes. We also examined if positive selection on centromeric histone H3 (CENH3), a protein that can suppress centromere drive, is more frequent in these asymmetric lineages.

METHODS

We analysed plant groups with different meiotic modes: asymmetric in gymnosperms and angiosperms, and symmetric in bryophytes, lycophytes and ferns. We selected species based on available CENH3 gene sequences and chromosome size data. Using Ornstein-Uhlenbeck evolutionary models and phylogenetic regressions, we assessed the rates of chromosome size evolution and the frequency of positive selection on CENH3 in these clades.

RESULTS

Our analyses showed that clades with asymmetric meiosis have a higher frequency of positive selection on CENH3 and increased rates of chromosome size evolution compared to symmetric clades.

CONCLUSIONS

Our findings support the hypothesis that centromere drive accelerates chromosome and genome size evolution, potentially also influencing the process of post-polyploid diploidization. We propose a model which in a single framework helps explain the stability of chromosome size in symmetric lineages (bryophytes, lycophytes and ferns) and its variability in asymmetric lineages (gymnosperms and angiosperms), providing a foundation for future research in plant genome evolution.

The chromatin remodeling factor OsINO80 promotes H3K27me3 and H3K9me2 deposition and maintains TE silencing in rice

The INO80 chromatin remodeling complex plays a critical role in shaping the dynamic chromatin environment. The diverse functions of the evolutionarily conserved INO80 complex have been widely reported. However, the role of INO80 in modulating the histone variant H2A.Z is controversial. Moreover, whether INO80 helps regulate heterochromatin remains unknown. Here, we characterize the regulatory effects of OsINO80 on protein-coding genes and transposable elements (TEs) in rice. Upon OsINO80 overexpression in rice, we found three types of OsINO80-occupied regions with different chromatin signatures: type I (enriched with H2A.Z), type II (enriched with H3K9me2), and type III (deficient in H2A.Z/H3K9me2). Loss of OsINO80 results in a decrease in H3K27me3, but not H2A.Z, at type I regions as well as a decrease in H3K9me2 at type II regions, which correlates with TE activation and transposition. Our findings reveal that OsINO80 facilitates H3K27me3 establishment, promotes H3K9me2 deposition, and maintains TE silencing.

Genome Studies in Amaranthus cruentus L. and A. hypochondriacus L. Based on Repeatomic and Cytogenetic Data

Amaranthus cruentus L. and Amaranthus hypochondriacus L. are valuable and promising food crops for multi-purpose use that are distributed worldwide in temperate, subtropical, and tropical zones. However, their karyotypes and genomic relationships still remain insufficiently studied. For the first time, a comparative repeatome analysis of A. cruentus and A. hypochondriacus was performed based on the available NGS data; bioinformatic analyses using RepeatExplorer/TAREAN pipelines; and chromosome FISH mapping of 45S rDNA, 5S rDNA, and the most abundant satellite DNAs. In the repeatomes of these species, interspecific variations in the amount of Ty3/Gypsy and Ty1/Copia retroelements, DNA transposons, ribosomal, and satellite DNA were detected. In the repeatomes of both species, shared satDNAs with high sequence similarity were identified. The chromosome distribution patterns of four effective molecular markers, 45S rDNA, 5S rDNA, AmC4, and AmC9, allowed us to identify all chromosome pairs in the species karyotypes, construct unique karyograms of A. cruentus and A. hypochondriacus, and confirm the close relationship between their genomes. These results are important for comparative karyotypic studies within the genus Amaranthus. Our findings demonstrated that cytogenomic analyses might provide important data on genomic relationships within Amaranthus and increase knowledge on genome organization in these valuable crops.

Detection and classification of long terminal repeat sequences in plant LTR-retrotransposons and their analysis using explainable machine learning

BACKGROUND

Long terminal repeats (LTRs) represent important parts of LTR retrotransposons and retroviruses found in high copy numbers in a majority of eukaryotic genomes. LTRs contain regulatory sequences essential for the life cycle of the retrotransposon. Previous experimental and sequence studies have provided only limited information about LTR structure and composition, mostly from model systems. To enhance our understanding of these key sequence modules, we focused on the contrasts between LTRs of various retrotransposon families and other genomic regions. Furthermore, this approach can be utilized for the classification and prediction of LTRs.

RESULTS

We used machine learning methods suitable for DNA sequence classification and applied them to a large dataset of plant LTR retrotransposon sequences. We trained three machine learning models using (i) traditional model ensembles (Gradient Boosting), (ii) hybrid convolutional/long and short memory network models, and (iii) a DNA pre-trained transformer-based model using k-mer sequence representation. All three approaches were successful in classifying and isolating LTRs in this data, as well as providing valuable insights into LTR sequence composition. The best classification (expressed as F1 score) achieved for LTR detection was 0.85 using the hybrid network model. The most accurate classification task was superfamily classification (F1=0.89) while the least accurate was family classification (F1=0.74). The trained models were subjected to explainability analysis. Positional analysis identified a mixture of interesting features, many of which had a preferred absolute position within the LTR and/or were biologically relevant, such as a centrally positioned TATA-box regulatory sequence, and TG..CA nucleotide patterns around both LTR edges.

CONCLUSIONS

Our results show that the models used here recognized biologically relevant motifs, such as core promoter elements in the LTR detection task, and a development and stress-related subclass of transcription factor binding sites in the family classification task. Explainability analysis also highlighted the importance of 5'- and 3'- edges in LTR identity and revealed need to analyze more than just dinucleotides at these ends. Our work shows the applicability of machine learning models to regulatory sequence analysis and classification, and demonstrates the important role of the identified motifs in LTR detection.

Activation characteristics of Ty3-retrotransposons after spaceflight and genetic stability of insertion sites in rice progeny

Introduction

The space environment is mutagenic and may induce genomic and phenotypic variations. Exploring the changes in transposon activity in the rice genome under space radiation is of great significance.

Methods

To analyze the activation characteristics of Ty3-retrotransposons and genetic stability of insertion sites in rice progeny after spaceflight, seeds of Nipponbare, DN416, and DN423 were exposed on board the SJ-10 recoverable satellite for 12.5 days. The differential methylation and transcription levels of Ty3-retrotransposons in the genome of Nipponbare's F0 generation after spaceflight, as well as the genetic stability of Ty3-retrotransposon insertion sites in DN416 and DN423 from F3 to F5 generations, was analyzed.

Results

The study found that the retrotransposons of ancient and young transposon families underwent demethylation from the tillering to heading stages of Nipponbare plants, which were F0 generation of space-exposed seeds, when the Nipponbare seeds were hit by single space high charge and energy (HZE) particles with LET ≥ 100 keV/μm. the transcription levels significantly increased in ancient transposon families (osr30, osr40, and rire10) and young transposon families (dagul, rn215-125, osr37, RLG_15, osr34, rire8, rire3, rire2, and hopi) (p ≤ 0.05) when LET > 100 keV/μm. Furthermore, the young Ty3-retrotransposons, which included the hopi, squiq, dasheng, rire2, rire3, rire8, osr34, rn_215-125, dagul, and RLG_15 families, underwent 1 to 8 transpositions in the F3 to F5 of DN416 and DN423 mutants, and some of these transposon insertion sites were stably inherited.

Discussion

The research holds great significance for understanding the activation characteristics of Ty3-retrotransposons in the rice genome induced by space radiation and the genetic characteristics of transposon insertion sites in its progeny.

A land plant phylogenetic framework for GLABROUS INFLORESCENCE STEMS (GIS), SUPERMAN, JAGGED and allies plus their TOPLESS co-repressor

Members of the plant specific family of C1-1i zincfinger transcriptionfactors (ZF-TFs), such as SUPERMAN, JAGGED, KNUCKLES or GIS,regulatediversedevelopmental processes including sexual reproduction. C1-1is consist of one zinc-finger and one to two EAR domains, connected by large intrinsically disordered regions (IDR). While the role of C1-i1 ZF-TFs in development processes is well known for some genes in Arabidopsis, rice or tomatoa comprehensive and broadphylogenetic background is lacking, yet knowledge of orthology is a requirement for a better understanding of C1-1i-Zf-TFs diverse roles in plants. Here, we provide a fine-grained and land plant wide classification of C1-1i sub-families and their known co-repressors TOPLESS and TOPLESS RELATED. Our work combines the identification of orthologous groups with Maximum-Likelihood phylogeny reconstructions and digital gene expression analyses mining high quality land plant genomes and transcriptomes to generate a comprehensive framework of C1-1i ZF-TF evolution. We show that C1-1i's are low to moderate copy genesand that orthologous genesonly partiallyhaveconserved sub-family and life cycle stage dependent expression pattern across land plants while others are highly diverged. Our workprovides the phylogenetic framework for C1-1i ZF-TFs, s and strengthen C1-1 ZF-TFs as a potential model for IDR-research in plants.

PTGS is dispensable for the initiation of epigenetic silencing of an active transposon in Arabidopsis

Transposable elements (TEs) are repressed in plants through transcriptional gene silencing (TGS), maintained epigenetic silencing marks such as DNA methylation. However, the mechanisms by which silencing is first installed remain poorly understood in plants. Small interfering (si)RNAs and post-transcriptional gene silencing (PTGS) are believed to mediate the initiation of TGS by guiding the first deposition of DNA methylation. To determine how this silencing installation works, we took advantage of ÉVADÉ (EVD), an endogenous retroelement in Arabidopsis, able to recapitulate true de novo silencing with a sequence of PTGS followed by a TGS. To test whether PTGS is required for TGS, we introduce active EVD into RNA-DEPENDENT-RNA-POLYMERASE-6 (RDR6) mutants, an essential PTGS component. EVD activity and silencing are monitored across several generations. In the absence of PTGS, silencing of EVD is still achieved through installation of RNA-directed DNA methylation (RdDM). Our study shows that PTGS is dispensable for de novo EVD silencing. Although we cannot rule out that PTGS might facilitate TGS, or control TE activity, initiation of epigenetic silencing can take place in its absence.

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.

Repeatome Analysis and Satellite DNA Chromosome Patterns in Hedysarum Species

The cosmopolitan genus Hedysarum L. (Fabaceae) is divided into sections Hedysarum, Stracheya, and Multicaulia. This genus includes many valuable medicinal, melliferous, and forage species. The species taxonomy and genome relationships within the sections are still unclear. We examined intra- and interspecific diversity in the section (sect.) Hedysarum based on repeatome analyses using NGS data, bioinformatic technologies, and chromosome FISH mapping of 35S rDNA, 5S rDNA, and the identified satellite DNA families (satDNAs). A comparison of repeatomes of H. alpinum, H. theinum, and H. flavescens revealed differences in their composition. However, similarity in sequences of most satDNAs indicated a close relationship between genomes within sect. Hedysarum. New effective satDNA chromosomal markers were detected, which is important for karyotype analyses within Hedysarum. Intra- and interspecific variability in the chromosomal distribution patterns of the studied markers were revealed, and species karyograms were constructed. These results provided new insight into the karyotype structures and genomic diversity within sect. Hedysarum, clarified the systematic position of H. sachalinense and H. arcticum, and confirmed the distant genomic relationships between species from sections Hedysarum and Multicaulia. Our findings are important for further comparative genome studies within the genus Hedysarum.

Analysis of pericentromere composition and structure elucidated the history of the Hieracium alpinum L. genome, revealing waves of transposable elements insertions

BACKGROUND

The centromere is one of the key regions of the eukaryotic chromosome. While maintaining its function, centromeric DNA may differ among closely related species. Here, we explored the composition and structure of the pericentromeres (a chromosomal region including a functional centromere) of Hieracium alpinum (Asteraceae), a member of one of the most diverse genera in the plant kingdom. Previously, we identified a pericentromere-specific tandem repeat that made it possible to distinguish reads within the Oxford Nanopore library attributed to the pericentromeres, separating them into a discrete subset and allowing comparison of the repeatome composition of this subset with the remaining genome.

RESULTS

We found that the main satellite DNA (satDNA) monomer forms long arrays of linear and block types in the pericentromeric heterochromatin of H. alpinum, and very often, single reads contain forward and reverse arrays and mirror each other. Beside the major, two new minor satDNA families were discovered. In addition to satDNAs, high amounts of LTR retrotransposons (TEs) with dominant of Tekay lineage, were detected in the pericentromeres. We were able to reconstruct four main TEs of the Ty3-gypsy and Ty1-copia superfamilies and compare their relative positions with satDNAs. The latter showed that the conserved domains (CDs) of the TE proteins are located between the newly discovered satDNAs, which appear to be parts of ancient Tekay LTRs that we were able to reconstruct. The dominant satDNA monomer shows a certain similarity to the GAG CD of the Angela retrotransposon.

CONCLUSIONS

The species-specific pericentromeric arrays of the H. alpinum genome are heterogeneous, exhibiting both linear and block type structures. High amounts of forward and reverse arrays of the main satDNA monomer point to multiple microinversions that could be the main mechanism for rapid structural evolution stochastically creating the uniqueness of an individual pericentromeric structure. The traces of TEs insertion waves remain in pericentromeres for a long time, thus "keeping memories" of past genomic events. We counted at least four waves of TEs insertions. In pericentromeres, TEs particles can be transformed into satDNA, which constitutes a background pool of minor families that, under certain conditions, can replace the dominant one(s).

Recurrent excision of a hAT-like transposable element in CmAPRR2 leads to the "shooting star" melon phenotype

The external appearance of fruit commodities is an essential trait that has profound effects on consumer preferences. A natural melon variety, characterized by an uneven and patchy arrangement of dark green streaks and spots on the white-skinned rind, resembles shooting stars streaking across the sky; thus, this variety is called "Shooting Star" (SS). To investigate the mechanism underlying the SS melon rind pattern, we initially discovered that the variegated dark green color results from chlorophyll accumulation on the white skin. We then constructed a segregation population by crossing a SS inbred line with a white rind (WR) inbred line and used bulk segregant analysis (BSA) revealed that the SS phenotype is controlled by a single dominant gene, CmAPRR2, which has been previously confirmed to determine dark green coloration. Further genomic analysis revealed a hAT-like transposable element (TE) inserted in CmAPRR2. This TE in CmAPRR2 is recurrently excised from rind tissues, activating the expression of CmAPRR2. This activation promotes the accumulation of chlorophyll, leading to the variegated dark green color on the rind, and ultimately resulting in the SS rind phenotype. Therefore, we propose that the SS phenotype results from the recurrent excision of the hAT-like TE in CmAPRR2.

The genome of Eleocharis vivipara elucidates the genetics of C3-C4 photosynthetic plasticity and karyotype evolution in the Cyperaceae

Eleocharis vivipara, an amphibious sedge in the Cyperaceae family, has several remarkable properties, most notably its alternate use of C3 photosynthesis underwater and C4 photosynthesis on land. However, the absence of genomic data has hindered its utility for evolutionary and genetic research. Here, we present a high-quality genome for E. vivipara, representing the first chromosome-level genome for the Eleocharis genus, with an approximate size of 965.22 Mb mainly distributed across 10 chromosomes. Its Hi-C pattern, chromosome clustering results, and one-to-one genome synteny across two subgroups indicates a tetraploid structure with chromosome count 2n = 4x = 20. Phylogenetic analysis suggests that E. vivipara diverged from Cyperus esculentus approximately 32.96 million years ago (Mya), and underwent a whole-genome duplication (WGD) about 3.5 Mya. Numerous fusion and fission events were identified between the chromosomes of E. vivipara and its close relatives. We demonstrate that E. vivipara has holocentromeres, a chromosomal feature which can maintain the stability of such chromosomal rearrangements. Experimental transplantation and cross-section studies showed its terrestrial culms developed C4 Kranz anatomy with increased number of chloroplasts in the bundle sheath (BS) cells. Gene expression and weighted gene co-expression network analysis (WGCNA) showed overall elevated expression of core genes associated with the C4 pathway, and significant enrichment of genes related to modified culm anatomy and photosynthesis efficiency. We found evidence of mixed nicotinamide adenine dinucleotide - malic enzyme and phosphoenolpyruvate carboxykinase type C4 photosynthesis in E. vivipara, and hypothesize that the evolution of C4 photosynthesis predates the WGD event. The mixed type is dominated by subgenome A and supplemented by subgenome B. Collectively, our findings not only shed light on the evolution of E. vivipara and karyotype within the Cyperaceae family, but also provide valuable insights into the transition between C3 and C4 photosynthesis, offering promising avenues for crop improvement and breeding.

The physiology of plants in the context of space exploration

The stress that the space environment can induce on plant physiology is of both abiotic and biotic nature. The abiotic space environment is characterized by ionizing radiation and altered gravity, geomagnetic field (GMF), pressure, and light conditions. Biotic interactions include both pathogenic and beneficial interactions. Here, we provide an overall picture of the effects of abiotic and biotic space-related factors on plant physiology. The knowledge required for the success of future space missions will lead to a better understanding of fundamental aspects of plant physiological responses, thus providing useful tools for plant breeding and agricultural practices on Earth.

Unveiling the Genome-Wide Consequences of Range Expansion and Mating System Transitions in Primula vulgaris

Genetic diversity is heterogeneously distributed among populations of the same species, due to the joint effects of multiple demographic processes, including range contractions and expansions, and mating systems shifts. Here, we ask how both processes shape genomic diversity in space and time in the classical Primula vulgaris model. This perennial herb originated in the Caucasus region and was hypothesized to have expanded westward following glacial retreat in the Quaternary. Moreover, this species is a long-standing model for mating system transitions, exemplified by shifts from heterostyly to homostyly. Leveraging a high-quality reference genome of the closely related Primula veris and whole-genome resequencing data from both heterostylous and homostylous individuals from populations encompassing a wide distribution of P. vulgaris, we reconstructed the demographic history of P. vulgaris. Results are compatible with the previously proposed hypothesis of range expansion from the Caucasus region approximately 79,000 years ago and suggest later shifts to homostyly following rather than preceding postglacial colonization of England. Furthermore, in accordance with population genetic theoretical predictions, both processes are associated with reduced genetic diversity, increased linkage disequilibrium, and reduced efficacy of purifying selection. A novel result concerns the contrasting effects of range expansion versus shift to homostyly on transposable elements, for the former, process is associated with changes in transposable element genomic content, while the latter is not. Jointly, our results elucidate how the interactions among range expansion, transitions to selfing, and Quaternary climatic oscillations shape plant evolution.

Telomere-to-Telomere Haplotype-Resolved Genomes of Agrocybe chaxingu Reveals Unique Genetic Features and Developmental Insights

Agrocybe chaxingu is a widely cultivated edible fungus in China, which is rich in nutrients and medicinal compounds. However, the lack of a high-quality genome hinders further research. In this study, we assembled the telomere-to-telomere genomes of two sexually compatible monokaryons (CchA and CchB) derived from a primarily cultivated strain AS-5. The genomes of CchA and CchB were 50.60 Mb and 51.66 Mb with contig N50 values of 3.95 Mb and 3.97 Mb, respectively. Each contained 13 complete chromosomes with telomeres at both ends. The high mapping rate, uniform genome coverage, high LAI score, all BUSCOs with 98.5%, and all base accuracy exceeding 99.999% indicated the high level of integrity and quality of these two assembled genomes. Comparison of the two genomes revealed that approximately 30% of the nucleotide sequences between homologous chromosomes were non-syntenic, including 19 translocations, 36 inversions, and 15 duplications. An additional gene CchA_000467 was identified at the Mat A locus of CchA, which was observed exclusively in the Cyclocybe cylindracea species complex. A total of 613 (4.26%) and 483 (3.4%) unique genes were identified in CchA and CchB, respectively, with over 80% of these being hypothetical proteins. Transcriptomic analysis revealed that the expression levels of unique genes in CchB were significantly higher than those in CchA, and both CchA and CchB had unique genes specifically expressed at stages of mycelium and fruiting body. It was indicated that the growth and development of the A. chaxingu strain AS-5 required the coordinated action of two different nuclei, with CchB potentially playing a more significant role. These findings contributed to a more profound comprehension of the growth and developmental processes of basidiomycetes.

Species of the Sections Hedysarum and Multicaulia of the Genus Hedysarum (Fabaceae): Taxonomy, Distribution, Chromosomes, Genomes, and Phylogeny

The genus Hedysarum L. (Fabaceae) includes about 200 species of annual and perennial herbs distributed in Asia, Europe, North Africa, and North America. Many species of this genus are valuable medicinal, melliferous, and forage resources. In this review, we consider the taxonomic history of the genus Hedysarum, the chromosomal organization of the species from the sections Hedysarum and Multicaulia, as well as phylogenetic relationships between these sections. According to morphological, genetic, and phylogenetic data, the genus Hedysarum is divided into three main sections: Hedysarum (= syn. Gamotion), Multicaulia, and Stracheya. In species of this genus, two basic chromosome numbers, x = 7 (section Hedysarum) and x = 8 (sections Multicaulia and Stracheya), were determined. The systematic positions of some species within the sections are still uncertain due to their morphological similarities. The patterns of distribution of molecular chromosomal markers (45S rDNA, 5S rDNA, and different satellite DNAs) in karyotypes of various Hedysarum species made it possible to determine their ploidy status and also specify genomic relationships within the sections Hedysarum and Multicaulia. Recent molecular phylogenetic studies clarified significantly the taxonomy and evolutionary development of the genus Hedysarum.

The super-pangenome of Populus unveils genomic facets for its adaptation and diversification in widespread forest trees

Understanding the underlying mechanisms and links between genome evolution and adaptive innovations stands as a key goal in evolutionary studies. Poplars, among the world's most widely distributed and cultivated trees, exhibit extensive phenotypic diversity and environmental adaptability. In this study, we present a genus-level super-pangenome comprising 19 Populus genomes, revealing the likely pivotal role of private genes in facilitating local environmental and climate adaptation. Through the integration of pangenomes with transcriptomes, methylomes, and chromatin accessibility mapping, we unveil that the evolutionary trajectories of pangenes and duplicated genes are closely linked to local genomic landscapes of regulatory and epigenetic architectures, notably CG methylation in gene-body regions. Further comparative genomic analyses have enabled the identification of 142 202 structural variants across species that intersect with a significant number of genes and contribute substantially to both phenotypic and adaptive divergence. We have experimentally validated a ∼180-bp presence/absence variant affecting the expression of the CUC2 gene, crucial for leaf serration formation. Finally, we developed a user-friendly web-based tool encompassing the multi-omics resources associated with the Populus super-pangenome (http://www.populus-superpangenome.com). Together, the present pioneering super-pangenome resource in forest trees not only aids in the advancement of breeding efforts of this globally important tree genus but also offers valuable insights into potential avenues for comprehending tree biology.

Transposable elements cause the loss of self-incompatibility in citrus

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

Ancient hybridization and repetitive element proliferation in the evolutionary history of the monocot genus Amomum (Zingiberaceae)

Genome size variation is a crucial aspect of plant evolution, influenced by a complex interplay of factors. Repetitive elements, which are fundamental components of genomic architecture, often play a role in genome expansion by selectively amplifying specific repeat motifs. This study focuses on Amomum, a genus in the ginger family (Zingiberaceae), known for its 4.4-fold variation in genome size. Using a robust methodology involving PhyloNet reconstruction, RepeatExplorer clustering, and repeat similarity-based phylogenetic network construction, we investigated the repeatome composition, analyzed repeat dynamics, and identified potential hybridization events within the genus. Our analysis confirmed the presence of four major infrageneric clades (A-D) within Amomum, with clades A-C exclusively comprising diploid species (2n = 48) and clade D encompassing both diploid and tetraploid species (2n = 48 and 96). We observed an increase in the repeat content within the genus, ranging from 84% to 89%, compared to outgroup species with 75% of the repeatome. The SIRE lineage of the Ty1-Copia repeat superfamily was prevalent in most analyzed ingroup genomes. We identified significant difference in repeatome structure between the basal Amomum clades (A, B, C) and the most diverged clade D. Our investigation revealed evidence of ancient hybridization events within Amomum, coinciding with a substantial proliferation of multiple repeat groups. This finding supports the hypothesis that ancient hybridization is a driving force in the genomic evolution of Amomum. Furthermore, we contextualize our findings within the broader context of genome size variations and repeatome dynamics observed across major monocot lineages. This study enhances our understanding of evolutionary processes within monocots by highlighting the crucial roles of repetitive elements in shaping genome size and suggesting the mechanisms that drive these changes.

Look4LTRs: a Long terminal repeat retrotransposon detection tool capable of cross species studies and discovering recently nested repeats

Plant genomes include large numbers of transposable elements. One particular type of these elements is flanked by two Long Terminal Repeats (LTRs) and can translocate using RNA. Such elements are known as LTR-retrotransposons; they are the most abundant type of transposons in plant genomes. They have many important functions involving gene regulation and the rise of new genes and pseudo genes in response to severe stress. Additionally, LTR-retrotransposons have several applications in biotechnology. Due to the abundance and the importance of LTR-retrotransposons, multiple computational tools have been developed for their detection. However, none of these tools take advantages of the availability of related genomes; they process one chromosome at a time. Further, recently nested LTR-retrotransposons (multiple elements of the same family are inserted into each other) cannot be annotated accurately - or cannot be annotated at all - by the currently available tools. Motivated to overcome these two limitations, we built Look4LTRs, which can annotate LTR-retrotransposons in multiple related genomes simultaneously and discover recently nested elements. The methodology of Look4LTRs depends on techniques imported from the signal-processing field, graph algorithms, and machine learning with a minimal use of alignment algorithms. Four plant genomes were used in developing Look4LTRs and eight plant genomes for evaluating it in contrast to three related tools. Look4LTRs is the fastest while maintaining better or comparable F1 scores (the harmonic average of recall and precision) to those obtained by the other tools. Our results demonstrate the added benefit of annotating LTR-retrotransposons in multiple related genomes simultaneously and the ability to discover recently nested elements. Expert human manual examination of six elements - not included in the ground truth - revealed that three elements belong to known families and two elements are likely from new families. With respect to examining recently nested LTR-retrotransposons, three out of five were confirmed to be valid elements. Look4LTRs - with its speed, accuracy, and novel features - represents a true advancement in the annotation of LTR-retrotransposons, opening the door to many studies focused on understanding their functions in plants.

Karyotype and LTR-RTs analysis provide insights into oak genomic evolution

BACKGROUND

Whole-genome duplication and long terminal repeat retrotransposons (LTR-RTs) amplification in organisms are essential factors that affect speciation, local adaptation, and diversification of organisms. Understanding the karyotype projection and LTR-RTs amplification could contribute to untangling evolutionary history. This study compared the karyotype and LTR-RTs evolution in the genomes of eight oaks, a dominant lineage in Northern Hemisphere forests.

RESULTS

Karyotype projections showed that chromosomal evolution was relatively conservative in oaks, especially on chromosomes 1 and 7. Modern oak chromosomes formed through multiple fusions, fissions, and rearrangements after an ancestral triplication event. Species-specific chromosomal rearrangements revealed fragments preserved through natural selection and adaptive evolution. A total of 441,449 full-length LTR-RTs were identified from eight oak genomes, and the number of LTR-RTs for oaks from section Cyclobalanopsis was larger than in other sections. Recent amplification of the species-specific LTR-RTs lineages resulted in significant variation in the abundance and composition of LTR-RTs among oaks. The LTR-RTs insertion suppresses gene expression, and the suppressed intensity in gene regions was larger than in promoter regions. Some centromere and rearrangement regions indicated high-density peaks of LTR/Copia and LTR/Gypsy. Different centromeric regional repeat units (32, 78, 79 bp) were detected on different Q. glauca chromosomes.

CONCLUSION

Chromosome fusions and arm exchanges contribute to the formation of oak karyotypes. The composition and abundance of LTR-RTs are affected by its recent amplification. LTR-RTs random retrotransposition suppresses gene expression and is enriched in centromere and chromosomal rearrangement regions. This study provides novel insights into the evolutionary history of oak karyotypes and the organization, amplification, and function of LTR-RTs.

Repeat turnover meets stable chromosomes: repetitive DNA sequences mark speciation and gene pool boundaries in sugar beet and wild beets

Sugar beet and its wild relatives share a base chromosome number of nine and similar chromosome morphologies. Yet, interspecific breeding is impeded by chromosome and sequence divergence that is still not fully understood. Since repetitive DNAs are among the fastest evolving parts of the genome, we investigated, if repeatome innovations and losses are linked to chromosomal differentiation and speciation. We traced genome and chromosome-wide evolution across 13 beet species comprising all sections of the genera Beta and Patellifolia. For this, we combined short and long read sequencing, flow cytometry, and cytogenetics to build a comprehensive framework that spans the complete scale from DNA to chromosome to genome. Genome sizes and repeat profiles reflect the separation into three gene pools with contrasting evolutionary patterns. Among all repeats, satellite DNAs harbor most genomic variability, leading to fundamentally different centromere architectures, ranging from chromosomal uniformity in Beta and Patellifolia to the formation of patchwork chromosomes in Corollinae/Nanae. We show that repetitive DNAs are causal for the genome expansions and contractions across the beet genera, providing insights into the genomic underpinnings of beet speciation. Satellite DNAs in particular vary considerably between beet genomes, leading to the evolution of distinct chromosomal setups in the three gene pools, likely contributing to the barriers in beet breeding. Thus, with their isokaryotypic chromosome sets, beet genomes present an ideal system for studying the link between repeats, genomic variability, and chromosomal differentiation and provide a theoretical fundament for understanding barriers in any crop breeding effort.

The global distribution of angiosperm genome size is shaped by climate

Angiosperms, which inhabit diverse environments across all continents, exhibit significant variation in genome sizes, making them an excellent model system for examining hypotheses about the global distribution of genome size. These include the previously proposed large genome constraint, mutational hazard, polyploidy-mediated, and climate-mediated hypotheses. We compiled the largest genome size dataset to date, encompassing 16 017 (> 5% of known) angiosperm species, and analyzed genome size distribution using a comprehensive geographic distribution dataset for all angiosperms. We observed that angiosperms with large range sizes generally had small genomes, supporting the large genome constraint hypothesis. Climate was shown to exert a strong influence on genome size distribution along the global latitudinal gradient, while the frequency of polyploidy and the type of growth form had negligible effects. In contrast to the unimodal patterns along the global latitudinal gradient shown by plant size traits and polyploid proportions, the increase in angiosperm genome size from the equator to 40-50°N/S is probably mediated by different (mostly climatic) mechanisms than the decrease in genome sizes observed from 40 to 50°N northward. Our analysis suggests that the global distribution of genome sizes in angiosperms is mainly shaped by climatically mediated purifying selection, genetic drift, relaxed selection, and environmental filtering.

Genomic evidence for rediploidization and adaptive evolution following the whole-genome triplication

Whole-genome duplication (WGD), or polyploidy, events are widespread and significant in the evolutionary history of angiosperms. However, empirical evidence for rediploidization, the major process where polyploids give rise to diploid descendants, is still lacking at the genomic level. Here we present chromosome-scale genomes of the mangrove tree Sonneratia alba and the related inland plant Lagerstroemia speciosa. Their common ancestor has experienced a whole-genome triplication (WGT) approximately 64 million years ago coinciding with a period of dramatic global climate change. Sonneratia, adapting mangrove habitats, experienced extensive chromosome rearrangements post-WGT. We observe the WGT retentions display sequence and expression divergence, suggesting potential neo- and sub-functionalization. Strong selection acting on three-copy retentions indicates adaptive value in response to new environments. To elucidate the role of ploidy changes in genome evolution, we improve a model of the polyploidization-rediploidization process based on genomic evidence, contributing to the understanding of adaptive evolution during climate change.

Regulatory mechanism of heat-active retrotransposons by the SET domain protein SUVH2

New transposon insertions are deleterious to genome stability. The RNA-directed DNA methylation (RdDM) pathway evolved to regulate transposon activity via DNA methylation. However, current studies have not yet clearly described the transposition regulation. ONSEN is a heat-activated retrotransposon that is activated at 37°C. The plant-specific SUPPRESSOR OF VARIEGATION 3-9 HOMOLOG (SUVH) family proteins function downstream of the RdDM pathway. The SUVH protein families are linked to TE silencing by two pathways, one through DNA methylation and the other through chromatin remodeling. In this study, we analyzed the regulation of ONSEN activity by SUVH2. We observed that ONSEN transcripts were increased; however, there was no transpositional activity in Arabidopsis suvh2 mutant. The suvh2 mutant produced siRNAs from the ONSEN locus under heat stress, suggesting that siRNAs are involved in suppressing transposition. These results provide new insights into the regulatory mechanisms of retrotransposons that involve siRNA in the RdDM pathway.

Duck pan-genome reveals two transposon insertions caused bodyweight enlarging and white plumage phenotype formation during evolution

Structural variations (SVs) are a major source of domestication and improvement traits. We present the first duck pan-genome constructed using five genome assemblies capturing ∼40.98 Mb new sequences. This pan-genome together with high-depth sequencing data (∼46.5×) identified 101,041 SVs, of which substantial proportions were derived from transposable element (TE) activity. Many TE-derived SVs anchoring in a gene body or regulatory region are linked to duck's domestication and improvement. By combining quantitative genetics with molecular experiments, we, for the first time, unraveled a 6945 bp Gypsy insertion as a functional mutation of the major gene IGF2BP1 associated with duck bodyweight. This Gypsy insertion, to our knowledge, explains the largest effect on bodyweight among avian species (27.61% of phenotypic variation). In addition, we also examined another 6634 bp Gypsy insertion in MITF intron, which triggers a novel transcript of MITF, thereby contributing to the development of white plumage. Our findings highlight the importance of using a pan-genome as a reference in genomics studies and illuminate the impact of transposons in trait formation and livestock breeding.

Complete mitochondrial genome of Syzygium samarangense reveals genomic recombination, gene transfer, and RNA editing events

Wax apple (Syzygium samarangense) is a commercial fruit that belongs to one of the most species-rich tree genera in the world. We report here the first complete S. samarangense mitogenome obtained using a hybrid assembly strategy. The mitogenome was a 530,242 bp circular molecule encoding 61 unique genes accounting for 7.99% of the full-length genome. Additionally, 167 simple sequence repeats, 19 tandem repeats, and 529 pairs of interspersed repeats were identified. Long read mapping and Sanger sequencing revealed the involvement of two forward repeats (35,843 bp and 22,925 bp) in mediating recombination. Thirteen homologous fragments in the chloroplast genome were identified, accounting for 1.53% of the mitogenome, and the longest fragment was 2,432 bp. An evolutionary analysis showed that S. samarangense underwent multiple genomic reorganization events and lost at least four protein-coding genes (PCGs) (rps2, rps7, rps11, and rps19). A total of 591 RNA editing sites were predicted in 37 PCGs, of which nad1-2, nad4L-2, and rps10-2 led to the gain of new start codons, while atp6-1156, ccmFC-1315 and rps10-331 created new stop codons. This study reveals the genetic features of the S. samarangense mitogenome and provides a scientific basis for further studies of traits with an epistatic basis and for germplasm identification.

High-quality genome assemblies for two Australimusa bananas (Musa spp.) and insights into regulatory mechanisms of superior fiber properties

Bananas (Musa spp.) are monocotyledonous plants with high genetic diversity in the Musaceae family that are cultivated mainly in tropical and subtropical countries. The fruits are a popular food, and the plants themselves have diverse uses. Four genetic groups (genomes) are thought to have contributed to current banana cultivars: Musa acuminata (A genome), Musa balbisiana (B genome), Musa schizocarpa (S genome), and species of the Australimusa section (T genome). However, the T genome has not been effectively explored. Here, we present the high-quality TT genomes of two representative accessions, Abaca (Musa textilis), with high-quality natural fiber, and Utafun (Musa troglodytarum, Fe'i group), with abundant β-carotene. Both the Abaca and Utafun assemblies comprise 10 pseudochromosomes, and their total genome sizes are 613 Mb and 619 Mb, respectively. Comparative genome analysis revealed that the larger size of the T genome is likely attributable to rapid expansion and slow removal of transposons. Compared with those of Musa AA or BB accessions or sisal (Agava sisalana), Abaca fibers exhibit superior mechanical properties, mainly because of their thicker cell walls with a higher content of cellulose, lignin, and hemicellulose. Expression of MusaCesA cellulose synthesis genes peaks earlier in Abaca than in AA or BB accessions during plant development, potentially leading to earlier cellulose accumulation during secondary cell wall formation. The Abaca-specific expressed gene MusaMYB26, which is directly regulated by MusaMYB61, may be an important regulator that promotes precocious expression of secondary cell wall MusaCesAs. Furthermore, MusaWRKY2 and MusaNAC68, which appear to be involved in regulating expression of MusaLAC and MusaCAD, may at least partially explain the high accumulation of lignin in Abaca. This work contributes to a better understanding of banana domestication and the diverse genetic resources in the Musaceae family, thus providing resources for Musa genetic improvement.

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.

Improved chromosome-level genome assembly of Indian sandalwood (Santalum album)

Santalum album is a well-known aromatic and medicinal plant that is highly valued for the essential oil (EO) extracted from its heartwood. In this study, we present a high-quality chromosome-level genome assembly of S. album after integrating PacBio Sequel, Illumina HiSeq paired-end and high-throughput chromosome conformation capture sequencing technologies. The assembled genome size is 207.39 M with a contig N50 of 7.33 M and scaffold N50 size of 18.31 M. Compared with three previously published sandalwood genomes, the N50 length of the genome assembly was longer. In total, 94.26% of the assembly was assigned to 10 pseudo-chromosomes, and the anchor rate far exceeded that of a recently released value. BUSCO analysis yielded a completeness score of 94.91%. In addition, we predicted 23,283 protein-coding genes, 89.68% of which were functionally annotated. This high-quality genome will provide a foundation for sandalwood functional genomics studies, and also for elucidating the genetic basis of EO biosynthesis in S. album.

A super-pangenome of the North American wild grape species

BACKGROUND

Capturing the genetic diversity of wild relatives is crucial for improving crops because wild species are valuable sources of agronomic traits that are essential to enhance the sustainability and adaptability of domesticated cultivars. Genetic diversity across a genus can be captured in super-pangenomes, which provide a framework for interpreting genomic variations.

RESULTS

Here we report the sequencing, assembly, and annotation of nine wild North American grape genomes, which are phased and scaffolded at chromosome scale. We generate a reference-unbiased super-pangenome using pairwise whole-genome alignment methods, revealing the extent of the genomic diversity among wild grape species from sequence to gene level. The pangenome graph captures genomic variation between haplotypes within a species and across the different species, and it accurately assesses the similarity of hybrids to their parents. The species selected to build the pangenome are a great representation of the genus, as illustrated by capturing known allelic variants in the sex-determining region and for Pierce's disease resistance loci. Using pangenome-wide association analysis, we demonstrate the utility of the super-pangenome by effectively mapping short reads from genus-wide samples and identifying loci associated with salt tolerance in natural populations of grapes.

CONCLUSIONS

This study highlights how a reference-unbiased super-pangenome can reveal the genetic basis of adaptive traits from wild relatives and accelerate crop breeding research.

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.

Comparative Analysis of Transposable Elements in Strawberry Genomes of Different Ploidy Levels

Transposable elements (TEs) make up a large portion of plant genomes and play a vital role in genome structure, function, and evolution. Cultivated strawberry (Fragaria x ananassa) is one of the most important fruit crops, and its octoploid genome was formed through several rounds of genome duplications from diploid ancestors. Here, we built a pan-genome TE library for the Fragaria genus using ten published strawberry genomes at different ploidy levels, including seven diploids, one tetraploid, and two octoploids, and performed comparative analysis of TE content in these genomes. The TEs comprise 51.83% (F. viridis) to 60.07% (F. nilgerrensis) of the genomes. Long terminal repeat retrotransposons (LTR-RTs) are the predominant TE type in the Fragaria genomes (20.16% to 34.94%), particularly in F. iinumae (34.94%). Estimating TE content and LTR-RT insertion times revealed that species-specific TEs have shaped each strawberry genome. Additionally, the copy number of different LTR-RT families inserted in the last one million years reflects the genetic distance between Fragaria species. Comparing cultivated strawberry subgenomes to extant diploid ancestors showed that F. vesca and F. iinumae are likely the diploid ancestors of the cultivated strawberry, but not F. viridis. These findings provide new insights into the TE variations in the strawberry genomes and their roles in strawberry genome evolution.

Introduction of Plant Transposon Annotation for Beginners

Transposons are mobile DNA sequences that contribute large fractions of many plant genomes. They provide exclusive resources for tracking gene and genome evolution and for developing molecular tools for basic and applied research. Despite extensive efforts, it is still challenging to accurately annotate transposons, especially for beginners, as transposon prediction requires necessary expertise in both transposon biology and bioinformatics. Moreover, the complexity of plant genomes and the dynamic evolution of transposons also bring difficulties for genome-wide transposon discovery. This review summarizes the three major strategies for transposon detection including repeat-based, structure-based, and homology-based annotation, and introduces the transposon superfamilies identified in plants thus far, and some related bioinformatics resources for detecting plant transposons. Furthermore, it describes transposon classification and explains why the terms 'autonomous' and 'non-autonomous' cannot be used to classify the superfamilies of transposons. Lastly, this review also discusses how to identify misannotated transposons and improve the quality of the transposon database. This review provides helpful information about plant transposons and a beginner's guide on annotating these repetitive sequences.

Complete mitogenome assembly of Selenicereus monacanthus revealed its molecular features, genome evolution, and phylogenetic implications

BACKGROUND

Mitochondria are the powerhouse of the cell and are critical for plant growth and development. Pitaya (Selenicereus or Hylocereus) is the most important economic crop in the family Cactaceae and is grown worldwide, however its mitogenome is unreported.

RESULTS

This study assembled the complete mitogenome of the red skin and flesh of pitaya (Selenicereus monacanthus). It is a full-length, 2,290,019 bp circular molecule encoding 59 unique genes that only occupy 2.17% of the entire length. In addition, 4,459 pairs of dispersed repeats (≥ 50 bp) were identified, accounting for 84.78% of the total length, and three repeats (394,588, 124,827, and 13,437 bp) mediating genomic recombination were identified by long read mapping and Sanger sequencing. RNA editing events were identified in all 32 protein-coding genes (PCGs), among which four sites (nad1-2, nad4L-2, atp9-copy3-223, and ccmFC-1309) were associated with the initiation or termination of PCGs. Seventy-eight homologous fragments of the chloroplast genome were identified in the mitogenome, the longest having 4,523 bp. In addition, evolutionary analyses suggest that S. monacanthus may have undergone multiple genomic reorganization events during evolution, with the loss of at least nine PCGs (rpl2, rpl10, rps2, rps3, rps10, rps11, rps14, rps19, and sdh3).

CONCLUSIONS

This study revealed the genetic basis of the S. monacanthus mitogenome, and provided a scientific basis for further research on phenotypic traits and germplasm resource development.

Genome-wide analysis of transposable elements and satellite DNA in Humulus scandens, a dioecious plant with XX/XY1Y2 chromosomes

Transposable elements (TEs) and satellite DNAs, two major categories of repetitive sequences, are expected to accumulate in non-recombining genome regions, including sex-linked regions, and contribute to sex chromosome evolution. The dioecious plant, Humulus scandens, can be used for studying the evolution of the XX/XY1Y2 sex chromosomes. In this study, we thoroughly examined the repetitive components of male and female H. scandens using next-generation sequencing data followed by bioinformatics analysis and florescence in situ hybridization (FISH). The H. scandens genome has a high overall repetitive sequence composition, 68.30% in the female and 66.78% in the male genome, with abundant long terminal repeat (LTR) retrotransposons (RTs), including more Ty3/Gypsy than Ty1/Copia elements, particularly two Ty3/Gypsy lineages, Tekay and Retand. Most LTR-RT lineages were found dispersed across the chromosomes, though CRM and Athila elements were predominately found within the centromeres and the pericentromeric regions. The Athila elements also showed clearly higher FISH signal intensities in the Y1 and Y2 chromosomes than in the X or autosomes. Three novel satellite DNAs were specifically distributed in the centromeric and/or telomeric regions, with markedly different distributions on the X, Y1, and Y2 chromosomes. Combined with FISH using satellite DNAs to stain chromosomes during meiotic diakinesis, we determined the synapsis pattern and distinguish pseudoautosomal regions (PARs). The results indicate that the XY1Y2 sex chromosomes of H. scandens might have originated from a centric fission event. This study improves our understanding of the repetitive sequence organization of H. scandens genome and provides a basis for further analysis of their chromosome evolution process.

MegaLTR: a web server and standalone pipeline for detecting and annotating LTR-retrotransposons in plant genomes

LTR-retrotransposons (LTR-RTs) are a class of RNA-replicating transposon elements (TEs) that can alter genome structure and function by moving positions, repositioning genes, shifting exons, and causing chromosomal rearrangements. LTR-RTs are widespread in many plant genomes and constitute a significant portion of the genome. Their movement and activity in eukaryotic genomes can provide insight into genome evolution and gene function, especially when LTR-RTs are located near or within genes. Building the redundant and non-redundant LTR-RTs libraries and their annotations for species lacking this resource requires extensive bioinformatics pipelines and expensive computing power to analyze large amounts of genomic data. This increases the need for online services that provide computational resources with minimal overhead and maximum efficiency. Here, we present MegaLTR as a web server and standalone pipeline that detects intact LTR-RTs at the whole-genome level and integrates multiple tools for structure-based, homologybased, and de novo identification, classification, annotation, insertion time determination, and LTR-RT gene chimera analysis. MegaLTR also provides statistical analysis and visualization with multiple tools and can be used to accelerate plant species discovery and assist breeding programs in their efforts to improve genomic resources. We hope that the development of online services such as MegaLTR, which can analyze large amounts of genomic data, will become increasingly important for the automated detection and annotation of LTR-RT elements.

Vitis labrusca genome assembly reveals diversification between wild and cultivated grapevine genomes

Wild grapevines are important genetic resources in breeding programs to confer adaptive fitness traits and unique fruit characteristics, but the genetics underlying these traits, and their evolutionary origins, are largely unknown. To determine the factors that contributed to grapevine genome diversification, we performed comprehensive intragenomic and intergenomic analyses with three cultivated European (including the PN40024 reference genome) and two wild North American grapevine genomes, including our newly released Vitis labrusca genome. We found the heterozygosity of the cultivated grapevine genomes was twice as high as the wild grapevine genomes studied. Approximately 30% of V. labrusca and 48% of V. vinifera Chardonnay genes were heterozygous or hemizygous and a considerable number of collinear genes between Chardonnay and V. labrusca had different gene zygosity. Our study revealed evidence that supports gene gain-loss events in parental genomes resulted in the inheritance of hemizygous genes in the Chardonnay genome. Thousands of segmental duplications supplied source material for genome-specific genes, further driving diversification of the genomes studied. We found an enrichment of recently duplicated, adaptive genes in similar functional pathways, but differential retention of environment-specific adaptive genes within each genome. For example, large expansions of NLR genes were discovered in the two wild grapevine genomes studied. Our findings support variation in transposable elements contributed to unique traits in grapevines. Our work revealed gene zygosity, segmental duplications, gene gain-and-loss variations, and transposable element polymorphisms can be key driving forces for grapevine genome diversification.

Abundance and Diversification of Repetitive Elements in Decapoda Genomes

Repetitive elements are a major component of DNA sequences due to their ability to propagate through the genome. Characterization of Metazoan repetitive profiles is improving; however, current pipelines fail to identify a significant proportion of divergent repeats in non-model organisms. The Decapoda order, for which repeat content analyses are largely lacking, is characterized by extremely variable genome sizes that suggest an important presence of repetitive elements. Here, we developed a new standardized pipeline to annotate repetitive elements in non-model organisms, which we applied to twenty Decapoda and six other Crustacea genomes. Using this new tool, we identified 10% more repetitive elements than standard pipelines. Repetitive elements were more abundant in Decapoda species than in other Crustacea, with a very large number of highly repeated satellite DNA families. Moreover, we demonstrated a high correlation between assembly size and transposable elements and different repeat dynamics between Dendrobranchiata and Reptantia. The patterns of repetitive elements largely reflect the phylogenetic relationships of Decapoda and the distinct evolutionary trajectories within Crustacea. In summary, our results highlight the impact of repetitive elements on genome evolution in Decapoda and the value of our novel annotation pipeline, which will provide a baseline for future comparative analyses.