Cycles of satellite and transposon evolution in Arabidopsis centromeres

A selfish supergene causes meiotic drive through both sexes in Drosophila

Meiotic drivers are selfish genetic elements that bias their own transmission during meiosis or gamete formation. Due to the fundamental differences between male and female meiosis in animals and plants, meiotic drivers operate through distinct mechanisms in the two sexes: In females, they exploit the asymmetry of meiosis to ensure their inclusion in the egg, whereas in males, they eliminate competing gametes after symmetric meiosis. Meiotic drive is commonly reported in males, where it strongly influences the evolution of spermatogenesis, while the few known cases in females have highlighted its crucial role in centromere evolution. Despite a growing number of examples in a wide range of organisms, meiotic drive has so far only been observed in one sex or the other since its discovery nearly 100 y ago. Here, we show that a selfish X chromosome known to cause meiotic drive in male Drosophila testacea flies also causes meiotic drive in females. We find that this X chromosome has supergene architecture, harboring extensive structural rearrangements that suppress recombination between the two X chromosomes. This has contributed to a substantial expansion of its size compared to the wild-type chromosome, partly due to the accumulation of species-specific repetitive elements. Our findings suggest that female meiotic drive may play an important role in the evolutionary dynamics of polymorphic structural variants that suppress recombination, including inversions, translocations, and supergenes.

A telomere-to-telomere genome assembly of cotton provides insights into centromere evolution and short-season adaptation

Cotton (Gossypium hirsutum L.) is a key allopolyploid crop with global economic importance. Here we present a telomere-to-telomere assembly of the elite variety Zhongmian 113. Leveraging technologies including PacBio HiFi, Oxford Nanopore Technology (ONT) ultralong-read sequencing and Hi-C, our assembly surpasses previous genomes in contiguity and completeness, resolving 26 centromeric and 52 telomeric regions, 5S rDNA clusters and nucleolar organizer regions. A phylogenetically recent centromere repositioning on chromosome D08 was discovered specific to G. hirsutum, involving deactivation of an ancestral centromere and the formation of a unique, satellite repeat-based centromere. Genomic analyses evaluated favorable allele aggregation for key agronomic traits and uncovered an early-maturing haplotype derived from an 11 Mb pericentric inversion that evolved early during G. hirsutum domestication. Our study sheds light on the genomic origins of short-season adaptation, potentially involving introgression of an inversion from primitively domesticated forms, followed by subsequent haplotype differentiation in modern breeding programs.

Post-polyploidization centromere evolution in cotton

Upland cotton (Gossypium hirsutum) accounts for more than 90% of the world's cotton production and, as an allotetraploid, is a model plant for polyploid crop domestication. In the present study, we reported a complete telomere-to-telomere (T2T) genome assembly of Upland cotton accession Texas Marker-1 (T2T-TM-1), which has a total size of 2,299.6 Mb, and annotated 79,642 genes. Based on T2T-TM-1, interspecific centromere divergence was detected between the A- and D-subgenomes and their corresponding diploid progenitors. Centromere-associated repetitive sequences (CRCs) were found to be enriched for Gypsy-like retroelements. Centromere size expansion, repositioning and structure variations occurred post-polyploidization. It is interesting that CRC homologs were transferred from the diploid D-genome progenitor to the D-subgenome, invaded the A-subgenome and then underwent post-tetraploidization proliferation. This suggests an evolutionary advantage for the CRCs of the D-genome progenitor, presents a D-genome-adopted inheritance of centromere repeats after polyploidization and shapes the dynamic centromeric landscape during polyploidization in polyploid species.

Analysis of 30 chromosome-level Drosophila genome assemblies reveals dynamic evolution of centromeric satellite repeats

BACKGROUND

The Drosophila genus is ideal for studying genome evolution due to its relatively simple chromosome structure and small genome size, with rearrangements mainly restricted to within chromosome arms, such as Muller elements. However, work on the rapidly evolving repetitive genomic regions, composed of transposons and tandem repeats, have been hampered by the lack of genus-wide chromosome-level assemblies.

RESULTS

Integrating long-read genomic sequencing and chromosome capture technology, here we produce and annotate 30 chromosome-level genome assemblies within the Drosophila genus. Based on this dataset, we reveal the evolutionary dynamics of genome rearrangements across the Drosophila phylogeny, including the identification of genomic regions that show comparatively high structural stability throughout evolution. Moreover, within the ananassae subgroup, we uncover the emergence of new chromosome conformations and the rapid expansion of novel satellite DNA sequence families, which form large and continuous pericentromeric domains with higher-order repeat structures that are reminiscent of those observed in the human and Arabidopsis genomes.

CONCLUSIONS

These chromosome-level genome assemblies present a valuable resource for future research, the power of which is demonstrated by our analysis of genome rearrangements and chromosome evolution. In addition, based on our findings, we propose the ananassae subgroup as an ideal model system for studying the evolution of centromere structure.

A multidisciplinary and integrative review of the structural genome and epigenome of Capsicum L. species

MAIN CONCLUSION

We revised and integrated the genomic and epigenomic data into a comparative Capsicum ideogram, evidencing the advances and future perspectives. Capsicum L. (Solanaceae) genome has been characterized concerning karyotype, nuclear and chromosomal genome size, genome sequencing and physical mapping. In addition, the epigenome has been investigated, showing chromosomal distribution of epimarks in histone amino acids. Genetic and epigenetic discoveries have given light to understanding the structure and organization of the Capsicum "omics". In addition, interspecific and intraspecific similarities and diversities have been identified, characterized and compared in taxonomic and evolutive scenarios. The journey through Capsicum studies allows us to know the 2n = 2x = 24 and 2n = 2x = 26 chromosome numbers, as well as the relatively homomorphic karyotype, and the 1C chromosomal DNA content. In addition, Capsicum "omics" diversity has mainly been evidenced from the nuclear 1C value, as well as from repeatome composition and mapping. Like this, Capsicum provides several opportunities for "omics", ecological, agronomic and conservation approaches, as well as subjects that can be used at different levels of education. In this context, we revisit and integrate Capsicum data about the genome size, karyotype, sequencing and cytogenomics, pointing out the progress and impact of this knowledge in taxonomic, evolutive and agronomic contexts. We also noticed gaps, which can be a focus of further studies. From this multidisciplinary and integrative review, we intend to show the beauty and intrigue of the Capsicum genome and epigenome, as well as the outcomes of these similarities and differences.

Convergent Patterns of Karyotype Evolution Underlying Karyotype Uniformity in Conifers

Karyotype diversity plays an important role in speciation and diversification. However, gymnosperms, particularly conifers, exhibit remarkable karyotype uniformity. To explore the evolutionary processes shaping karyotypes in gymnosperms, the karyotype evolutionary history is reconstructed through comparative genomic analyses. Synteny analysis confirms the absence of ancient polyploidy in conifers and its rarity across the gymnosperms as a whole. Further analysis reveals convergent patterns of reciprocal translocations between nonhomologous chromosomes in conifer genomes. Centromeric-centromeric reciprocal translocations (CRTs) have been identified as the primary mechanism of karyotype evolution in conifers, while telomeric-centromeric reciprocal translocations (TRTs) significantly contributed to descending dysploidy within Cupressales. A graph-based method is utilized to infer the detailed evolutionary pathways from the proto-gymnosperm karyotype (n = 12) to modern conifer karyotypes (n = 11-12). In conclusion, the scarcity of both polyploidy and dysploidy contributes to the karyotype uniformity of gymnosperms and potentially also to their lower species richness compared to angiosperms. However, the pervasive CRTs and occasional TRTs underlie this "apparent uniformity", supporting the "karyotype orthoselection" hypothesis. This study provides new insights into the mechanisms maintaining karyotype uniformity in conifers and the role of karyotype evolution in their diversification.

High-quality genome of black wolfberry (Lycium ruthenicum Murr.) provides insights into the genetics of anthocyanin biosynthesis regulation

Black wolfberry (Lycium ruthenicum Murr.) is an important plant for ecological preservation. In addition, its fruits are rich in anthocyanins and have important edible and medicinal value. However, a high-quality chromosome-level genome for this species is not yet available, and the regulatory mechanisms involved in the biosynthesis of anthocyanins are unclear. In this study, haploid material was used to assemble a high-quality chromosome-level reference genome of Lycium ruthenicum, resulting in a genome size of 2272 Mb with contig N50 of 92.64 Mb, and 38 993 annotated gene models. In addition, the evolution of this genome and large-scale variations compared with the Ningxia wolfberry Lycium barbarum were determined. Importantly, homology annotation identified 86 genes involved in the regulatory pathway of anthocyanin biosynthesis, five of which [LrCHS1 (evm.TU.Chr05.295), LrCHS2 (evm.TU.Chr09.488), LrAOMT (evm.TU.Chr09.809), LrF3'5'H (evm.TU.Chr06.177), and LrAN2.1 (evm.TU.Chr05.2618)] were screened by differential expression analysis and correlation analysis using a combination of transcriptome and metabolome testing. Overexpression of these genes could significantly up- or downregulate anthocyanin-related metabolites. These results will help accelerate the functional genomic research of L. ruthenicum, and the elucidation of the genes involved in anthocyanin synthesis will be beneficial for breeding new varieties and further exploring its ecological conservation potential.

Extensive genome evolution distinguishes maize within a stable tribe of grasses

Over the last 20 million years, the Andropogoneae tribe of grasses has evolved to dominate 17% of global land area. Domestication of these grasses in the last 10,000 years has yielded our most productive crops, including maize, sugarcane, and sorghum. The majority of Andropogoneae species, including maize, show a history of polyploidy - a condition that, while offering the evolutionary advantage of multiple gene copies, poses challenges to basic cellular processes, gene expression, and epigenetic regulation. Genomic studies of polyploidy have been limited by sparse sampling of taxa in groups with multiple polyploidy events. Here, we present 33 genome assemblies from 27 species, including chromosome-scale assemblies of maize relatives Zea and Tripsacum. In maize, the after-effects of polyploidy have been widely studied, showing reduced chromosome number, biased fractionation of duplicate genes, and transposable element (TE) expansions. While we observe these patterns within the genus Zea, 12 other polyploidy events deviate significantly. Those tetraploids and hexaploids retain elevated chromosome number, maintain nearly complete complements of duplicate genes, and have only stochastic TE amplifications. These genomes reveal variable outcomes of polyploidy, challenging simple predictions and providing a foundation for understanding its evolutionary implications in an ecologically and economically important clade.

Unravelling genomic drivers of speciation in Musa through genome assemblies of wild banana ancestors

Hybridization between wild Musa species and subspecies from Southeast Asia is at the origin of cultivated bananas. The genomes of these cultivars are complex mosaics involving nine genetic groups, including two previously unknown contributors. This study provides continuous genome assemblies for six wild genetic groups, one of which represents one of the unknown ancestor, identified as M.acuminata ssp. halabanensis. The second unknown ancestor partially present in a seventh assembly appears related to M. a. ssp. zebrina. These assemblies provide key resources for banana genetics and for improving cultivar assemblies, including that of the emblematic triploid Cavendish. Comparative and phylogenetic analyses reveal an ongoing speciation process within Musa, characterised by large chromosome rearrangements and centromere differentiation through the integration of different types of repeated sequences, including rDNA tandem repeats. This speciation process may have been favoured by reproductive isolation related to the particular context of climate and land connectivity fluctuations in the Southeast Asian region.

TIPPo: A User-Friendly Tool for De Novo Assembly of Organellar Genomes with High-Fidelity Data

Plant cells have two major organelles with their own genomes: chloroplasts and mitochondria. While chloroplast genomes tend to be structurally conserved, the mitochondrial genomes of plants, which are much larger than those of animals, are characterized by complex structural variation. We introduce TIPPo, a user-friendly, reference-free assembly tool that uses PacBio high-fidelity long-read data and that does not rely on genomes from related species or nuclear genome information for the assembly of organellar genomes. TIPPo employs a deep learning model for initial read classification and leverages k-mer counting for further refinement, significantly reducing the impact of nuclear insertions of organellar DNA on the assembly process. We used TIPPo to completely assemble a set of 54 complete chloroplast genomes. No other tool was able to completely assemble this set. TIPPo is comparable with PMAT in assembling mitochondrial genomes from most species but does achieve even higher completeness for several species. We also used the assembled organelle genomes to identify instances of nuclear plastid DNA (NUPTs) and nuclear mitochondrial DNA (NUMTs) insertions. The cumulative length of NUPTs/NUMTs positively correlates with the size of the nuclear genome, suggesting that insertions occur stochastically. NUPTs/NUMTs show predominantly C:G to T:A changes, with the mutated cytosines typically found in CG and CHG contexts, suggesting that degradation of NUPT and NUMT sequences is driven by the known elevated mutation rate of methylated cytosines. Small interfering RNA loci are enriched in NUPTs and NUMTs, consistent with the RdDM pathway mediating DNA methylation in these sequences.

Transposon-triggered epigenetic chromatin dynamics modulate EFR-related pathogen response

Infectious diseases drive wild plant evolution and impact crop yield. Plants, like animals, sense biotic threats through pattern recognition receptors (PRRs). Overly robust immune responses can harm plants; thus, understanding the tuning of defense response mechanisms is crucial for developing pathogen-resistant crops. In this study, we found that an inverted-repeat transposon (EFR-associated IR, Ea-IR) located between the loci encoding PRRs ELONGATION FACTOR-TU RECEPTOR (EFR) and myosin XI-k (XI-k) in Arabidopsis affects chromatin organization, promoting the formation of a repressive chromatin loop. Upon pathogen infection, chromatin changes around EFR and XI-k correlate with increased EFR transcription. Pathogen-induced chromatin opening causes RNA polymerase II readthrough, producing a longer, Ea-IR-containing XI-k transcript, processed by Dicer-like enzymes into small RNAs, which reset chromatin to a repressive state attenuating the immune response after infection. Arabidopsis accessions lacking Ea-IR have higher basal EFR levels and resistance to pathogens. We show a scenario in which a transposon, chromatin organization and gene expression interact to fine-tune immune responses, during both the course of infection and the course of evolution.

Centrophilic retrotransposon integration via CENH3 chromatin in Arabidopsis

In organisms ranging from vertebrates to plants, major components of centromeres are rapidly evolving repeat sequences, such as tandem repeats (TRs) and transposable elements (TEs), which harbour centromere-specific histone H3 (CENH3)1,2. Complete centromere structures recently determined in human and Arabidopsis suggest frequent integration and purging of retrotransposons within the TR regions of centromeres3-5. Despite the high impact of 'centrophilic' retrotransposons on the paradox of rapid centromere evolution, the mechanisms involved in centromere targeting remain poorly understood in any organism. Here we show that both Ty3 and Ty1 long terminal repeat retrotransposons rapidly turnover within the centromeric TRs of Arabidopsis species. We demonstrate that the Ty1/Copia element Tal1 (Transposon of Arabidopsis lyrata 1) integrates de novo into regions occupied by CENH3 in Arabidopsis thaliana, and that ectopic expansion of the CENH3 region results in spread of Tal1 integration regions. The integration spectra of chimeric TEs reveal the key structural variations responsible for contrasting chromatin-targeting specificities to centromeres versus gene-rich regions, which have recurrently converted during the evolution of these TEs. Our findings show the impact of centromeric chromatin on TE-mediated rapid centromere evolution, with relevance across eukaryotic genomes.

Identification of novel genes responsible for a pollen killer present in local natural populations of Arabidopsis thaliana

Gamete killers are genetic loci that distort segregation in the progeny of hybrids because the killer allele promotes the elimination of the gametes that carry the sensitive allele. They are widely distributed in eukaryotes and are important for understanding genome evolution and speciation. We had previously identified a pollen killer in hybrids between two distant natural accessions of Arabidopsis thaliana. This pollen killer involves three genetically linked genes, and we previously reported the identification of the gene encoding the antidote that protects pollen grains from the killer activity. In this study, we identified the two other genes of the pollen killer by using CRISPR-Cas9 induced mutants. These two genes are necessary for the killer activity that we demonstrated to be specific to pollen. The cellular localization of the pollen killer encoded proteins suggests that the pollen killer activity involves the mitochondria. Sequence analyses reveal predicted domains from the same families in the killer proteins. In addition, the C-terminal half of one of the killer proteins is identical to the antidote, and one amino acid, crucial for the antidote activity, is also essential for the killer function. Investigating more than 700 worldwide accessions of A. thaliana, we confirmed that the locus is subject to important structural rearrangements and copy number variation. By exploiting available de novo genomic sequences, we propose a scenario for the emergence of this pollen killer in A. thaliana. Furthermore, we report the co-occurrence and behavior of killer and sensitive genotypes in several local populations, a prerequisite for studying gamete killer evolution in the wild. This highlights the potential of the Arabidopsis model not only for functional studies of gamete killers but also for investigating their evolutionary trajectories at complementary geographical scales.

Gretl-variation GRaph Evaluation TooLkit

MOTIVATION

As genome graphs are powerful data structures for representing the genetic diversity within populations, they can help identify genomic variations that traditional linear references miss, but their complexity and size makes the analysis of genome graphs challenging. We sought to develop a genome graph analysis tool that helps these analyses to become more accessible by addressing the limitations of existing tools. Specifically, we improve scalability and user-friendliness, and we provide many new statistics tailored to variation graphs for graph evaluation, including sample-specific features.

RESULTS

We developed an efficient, comprehensive, and integrated tool, gretl, to analyze genome graphs and gain insights into their structure and composition by providing a wide range of statistics. gretl can be utilized to evaluate different graphs, compare the output of graph construction pipelines with different parameters, as well as perform an in-depth analysis of individual graphs, including sample-specific analysis. With the assistance of gretl, novel patterns of genetic variation and potential regions of interest can be identified, for later, more detailed inspection. We demonstrate that gretl outperforms other tools in terms of speed, particularly for larger genome graphs.

AVAILABILITY AND IMPLEMENTATION

Commented Rust source code and documentation is available under MIT license at https://github.com/MoinSebi/gretl together with Python scripts and step-by-step usage examples. The package is available at Bioconda for easy installation.

The complete genome assembly of Nicotiana benthamiana reveals the genetic and epigenetic landscape of centromeres

Nicotiana benthamiana is a model organism widely adopted in plant biology. Its complete assembly remains unavailable despite several recent improvements. To further improve its usefulness, we generate and phase the complete 2.85 Gb genome assembly of allotetraploid N. benthamiana. We find that although Solanaceae centromeres are widely dominated by Ty3/Gypsy retrotransposons, satellite-based centromeres are surprisingly common in N. benthamiana, with 11 of 19 centromeres featured by megabase-scale satellite arrays. Interestingly, the satellite-enriched and satellite-free centromeres are extensively invaded by distinct Gypsy retrotransposons which CENH3 protein more preferentially occupies, suggestive of their crucial roles in centromere function. We demonstrate that ribosomal DNA is a major origin of centromeric satellites, and mitochondrial DNA could be employed as a core component of the centromere. Subgenome analysis indicates that the emergence of satellite arrays probably drives new centromere formation. Altogether, we propose that N. benthamiana centromeres evolved via neocentromere formation, satellite expansion, retrotransposon enrichment and mtDNA integration.

Bridging the gap: unravelling plant centromeres in the telomere-to-telomere era

Centromeres are specific regions of the chromosomes that play a pivotal role in the segregation of chromosomes, by facilitating the loading of the kinetochore, which forms the link between the chromosomes to the spindle fibres during cell division. In plants and animals, these regions often form megabase-scale loci of tandemly repeated DNA sequences, which have presented a challenge to genomic studies even in model species. The functional designation of centromeres is determined epigenetically by the incorporation of a centromere-specific variant of histone H3. Recent developments in long-read sequencing technology have allowed the assembly of these regions for the first time and have prompted a reassessment of fidelity of centromere function and the evolutionary dynamics of these regions.

Estimating dispersal rates and locating genetic ancestors with genome-wide genealogies

Spatial patterns in genetic diversity are shaped by individuals dispersing from their parents and larger-scale population movements. It has long been appreciated that these patterns of movement shape the underlying genealogies along the genome leading to geographic patterns of isolation-by-distance in contemporary population genetic data. However, extracting the enormous amount of information contained in genealogies along recombining sequences has, until recently, not been computationally feasible. Here, we capitalize on important recent advances in genome-wide gene-genealogy reconstruction and develop methods to use thousands of trees to estimate per-generation dispersal rates and to locate the genetic ancestors of a sample back through time. We take a likelihood approach in continuous space using a simple approximate model (branching Brownian motion) as our prior distribution of spatial genealogies. After testing our method with simulations we apply it to Arabidopsis thaliana. We estimate a dispersal rate of roughly 60 km2/generation, slightly higher across latitude than across longitude, potentially reflecting a northward post-glacial expansion. Locating ancestors allows us to visualize major geographic movements, alternative geographic histories, and admixture. Our method highlights the huge amount of information about past dispersal events and population movements contained in genome-wide genealogies.

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

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.

Turnover of retroelements and satellite DNA drives centromere reorganization over short evolutionary timescales in Drosophila

Centromeres reside in rapidly evolving, repeat-rich genomic regions, despite their essential function in chromosome segregation. Across organisms, centromeres are rich in selfish genetic elements such as transposable elements and satellite DNAs that can bias their transmission through meiosis. However, these elements still need to cooperate at some level and contribute to, or avoid interfering with, centromere function. To gain insight into the balance between conflict and cooperation at centromeric DNA, we take advantage of the close evolutionary relationships within the Drosophila simulans clade-D. simulans, D. sechellia, and D. mauritiana-and their relative, D. melanogaster. Using chromatin profiling combined with high-resolution fluorescence in situ hybridization on stretched chromatin fibers, we characterize all centromeres across these species. We discovered dramatic centromere reorganization involving recurrent shifts between retroelements and satellite DNAs over short evolutionary timescales. We also reveal the recent origin (<240 Kya) of telocentric chromosomes in D. sechellia, where the X and fourth centromeres now sit on telomere-specific retroelements. Finally, the Y chromosome centromeres, which are the only chromosomes that do not experience female meiosis, do not show dynamic cycling between satDNA and TEs. The patterns of rapid centromere turnover in these species are consistent with genetic conflicts in the female germline and have implications for centromeric DNA function and karyotype evolution. Regardless of the evolutionary forces driving this turnover, the rapid reorganization of centromeric sequences over short evolutionary timescales highlights their potential as hotspots for evolutionary innovation.

Retrotransposon-driven environmental regulation of FLC leads to adaptive response to herbicide

The mobilization of transposable elements is a potent source of mutations. In plants, several stransposable elements respond to external cues, fuelling the hypothesis that natural transposition can create environmentally sensitive alleles for adaptation. Here we report on the detailed characterization of a retrotransposon insertion within the first intron of the Arabidopsis floral-repressor gene FLOWERING LOCUS C (FLC) and the discovery of its role for adaptation. The insertion mutation augments the environmental sensitivity of FLC by affecting the balance between coding and non-coding transcripts in response to stress, thus expediting flowering. This balance is modulated by DNA methylation and orchestrated by IBM2, a factor involved in the processing of intronic heterochromatic sequences. The stress-sensitive allele of FLC has spread across populations subjected to recurrent chemical weeding, and we show that retrotransposon-driven acceleration of the life cycle represents a rapid response to herbicide application. Our work provides a compelling example of a transposable element-driven environmentally sensitive allele that confers an adaptive response in nature.

Exploring the complexity of genome size reduction in angiosperms

The genome sizes of angiosperms decreased significantly more than the genome sizes of their ancestors (pteridophytes and gymnosperms). Decreases in genome size involve a highly complex process, with remnants of the genome size reduction scattered across the genome and not directly linked to specific genomic structures. This is because the associated mechanisms operate on a much smaller scale than the mechanisms mediating increases in genome size. This review thoroughly summarizes the available literature regarding the molecular mechanisms underlying genome size reductions and introduces Utricularia gibba and Arabidopsis thaliana as model species for the examination of the effects of these molecular mechanisms. Additionally, we propose that phosphorus deficiency and drought stress are the major external factors contributing to decreases in genome size. Considering these factors affect almost all land plants, angiosperms likely gained the mechanisms for genome size reductions. These environmental factors may affect the retention rates of deletions, while also influencing the mutation rates of deletions via the functional diversification of the proteins facilitating double-strand break repair. The biased retention and mutation rates of deletions may have synergistic effects that enhance deletions in intergenic regions, introns, transposable elements, duplicates, and repeats, leading to a rapid decrease in genome size. We suggest that these selection pressures and associated molecular mechanisms may drive key changes in angiosperms during recurrent cycles of genome size decreases and increases.

Plant kinetochore complex: composition, function, and regulation

The kinetochore complex, an important protein assembly situated on the centromere, plays a pivotal role in chromosome segregation during cell division. Like in animals and fungi, the plant kinetochore complex is important for maintaining chromosome stability, regulating microtubule attachment, executing error correction mechanisms, and participating in signaling pathways to ensure accurate chromosome segregation. This review summarizes the composition, function, and regulation of the plant kinetochore complex, emphasizing the interactions of kinetochore proteins with centromeric DNAs (cenDNAs) and RNAs (cenRNAs). Additionally, the applications of the centromeric histone H3 variant (the core kinetochore protein CENH3, first identified as CENP-A in mammals) in the generation of ploidy-variable plants and synthesis of plant artificial chromosomes (PACs) are discussed. The review serves as a comprehensive roadmap for researchers delving into plant kinetochore exploration, highlighting the potential of kinetochore proteins in driving technological innovations in synthetic genomics and plant biotechnology.

The loach haplotype-resolved genome and the identification of Mex3a involved in fish air breathing

Fish air breathing is crucial for the transition of vertebrates from water to land. So far, the genes involved in fish air breathing have not been well identified. Here, we performed gene enrichment analysis of positively selected genes (PSGs) in loach (Misgurnus anguillicaudatus, an air-breathing fish) in comparison to Triplophysa tibetana (a non-air-breathing fish), haplotype-resolved genome assembly of the loach, and gene evolutionary analysis of air-breathing and non-air-breathing fishes and found that the PSG mex3a originated from ancient air-breathing fish species. Deletion of Mex3a impaired loach air-breathing capacity by inhibiting angiogenesis through its interaction with T-box transcription factor 20. Mex3a overexpression significantly promoted angiogenesis. Structural analysis and point mutation revealed the critical role of the 201st amino acid in loach Mex3a for angiogenesis. Our findings innovatively indicate that the ancient mex3a is a fish air-breathing gene, which holds significance for understanding fish air breathing and provides a valuable resource for cultivating hypoxia-tolerant fish varieties.

Unlocking plant genetics with telomere-to-telomere genome assemblies

Contiguous genome sequence assemblies will help us to realize the full potential of crop translational genomics. Recent advances in sequencing technologies, especially long-read sequencing strategies, have made it possible to construct gapless telomere-to-telomere (T2T) assemblies, thus offering novel insights into genome organization and function. Plant genomes pose unique challenges, such as a continuum of ancient to recent polyploidy and abundant highly similar and long repetitive elements. Owing to progress in sequencing approaches, for most crop plants, chromosome-scale reference genome assemblies are available, but T2T assembly construction remains challenging. Here we describe methods for haplotype-resolved, gapless T2T assembly construction in plants, including various crop species. We outline the impact of T2T assemblies in elucidating the roles of repetitive elements in gene regulation, as well as in pangenomics, functional genomics, genome-assisted breeding and targeted genome manipulation. In conjunction with sequence-enriched germplasm repositories, T2T assemblies thus hold great promise for basic and applied plant sciences.

A telomere-to-telomere cotton genome assembly reveals centromere evolution and a Mutator transposon-linked module regulating embryo development

Assembly of complete genomes can reveal functional genetic elements missing from draft sequences. Here we present the near-complete telomere-to-telomere and contiguous genome of the cotton species Gossypium raimondii. Our assembly identified gaps and misoriented or misassembled regions in previous assemblies and produced 13 centromeres, with 25 chromosomal ends having telomeres. In contrast to satellite-rich Arabidopsis and rice centromeres, cotton centromeres lack phased CENH3 nucleosome positioning patterns and probably evolved by invasion from long terminal repeat retrotransposons. In-depth expression profiling of transposable elements revealed a previously unannotated DNA transposon (MuTC01) that interacts with miR2947 to produce trans-acting small interfering RNAs (siRNAs), one of which targets the newly evolved LEC2 (LEC2b) to produce phased siRNAs. Systematic genome editing experiments revealed that this tripartite module, miR2947-MuTC01-LEC2b, controls the morphogenesis of complex folded embryos characteristic of Gossypium and its close relatives in the cotton tribe. Our study reveals a trans-acting siRNA-based tripartite regulatory pathway for embryo development in higher plants.

Retrotransposon addiction promotes centromere function via epigenetically activated small RNAs

Retrotransposons have invaded eukaryotic centromeres in cycles of repeat expansion and purging, but the function of centromeric retrotransposons has remained unclear. In Arabidopsis, centromeric ATHILA retrotransposons give rise to epigenetically activated short interfering RNAs in mutants in DECREASE IN DNA METHYLATION1 (DDM1). Here we show that mutants that lose both DDM1 and RNA-dependent RNA polymerase have pleiotropic developmental defects and mis-segregate chromosome 5 during mitosis. Fertility and segregation defects are epigenetically inherited with centromere 5, and can be rescued by directing artificial small RNAs to ATHILA5 retrotransposons that interrupt tandem satellite repeats. Epigenetically activated short interfering RNAs promote pericentromeric condensation, chromosome cohesion and chromosome segregation in mitosis. We propose that insertion of ATHILA silences centromeric transcription, while simultaneously making centromere function dependent on retrotransposon small RNAs in the absence of DDM1. Parallels are made with the fission yeast Schizosaccharomyces pombe, where chromosome cohesion depends on RNA interference, and with humans, where chromosome segregation depends on both RNA interference and HELLSDDM1.

The gap-free genome of Forsythia suspensa illuminates the intricate landscape of centromeres

Forsythia suspensa, commonly known as weeping forsythia, holds significance in traditional medicine and horticulture. Despite its ecological and cultural importance, the existing reference genome presents challenges with duplications and gaps, hindering in-depth genomic analyses. Here, we present a Telomere-to-Telomere (T2T) assembly of the F. suspensa genome, integrating Oxford Nanopore Technologies (ONT) ultra-long, Hi-C datasets, and high-fidelity (HiFi) sequencing data. The T2T reference genome (Fsus-CHAU) consists of 14 chromosomes, totaling 688.79 Mb, and encompasses 33 932 predicted protein-coding genes. Additionally, we characterize functional centromeres in the F. suspensa genome by developing a specific CENH3 antibody. We demonstrate that centromeric regions in F. suspensa exhibit a diverse array of satellites, showcasing distinctive types with unconventional lengths across various chromosomes. This discovery offers implications for the adaptability of CENH3 and the potential influence on centromere dynamics. Furthermore, after assessing the insertion time of full-length LTRs within centromeric regions, we found that they are older compared to those across the entire genome, contrasting with observations in other species where centromeric retrotransposons are typically young. We hypothesize that asexual reproduction may impact retrotransposon dynamics, influencing centromere evolution. In conclusion, our T2T assembly of the F. suspensa genome, accompanied by detailed genomic annotations and centromere analysis, significantly enhances F. suspensa potential as a subject of study in fields ranging from ecology and horticulture to traditional medicine.

The centromere landscapes of four karyotypically diverse Papaver species provide insights into chromosome evolution and speciation

Understanding the roles played by centromeres in chromosome evolution and speciation is complicated by the fact that centromeres comprise large arrays of tandemly repeated satellite DNA, which hinders high-quality assembly. Here, we used long-read sequencing to generate nearly complete genome assemblies for four karyotypically diverse Papaver species, P. setigerum (2n = 44), P. somniferum (2n = 22), P. rhoeas (2n = 14), and P. bracteatum (2n = 14), collectively representing 45 gapless centromeres. We identified four centromere satellite (cenSat) families and experimentally validated two representatives. For the two allopolyploid genomes (P. somniferum and P. setigerum), we characterized the subgenomic distribution of each satellite and identified a "homogenizing" phase of centromere evolution in the aftermath of hybridization. An interspecies comparison of the peri-centromeric regions further revealed extensive centromere-mediated chromosome rearrangements. Taking these results together, we propose a model for studying cenSat competition after hybridization and shed further light on the complex role of the centromere in speciation.

Chromosome-Scale Assembly and Annotation of Eight Arabidopsis thaliana Ecotypes

The plant Arabidopsis thaliana is a model system used by researchers through much of plant research. Recent efforts have focused on discovering the genomic variation found in naturally occurring ecotypes isolated from around the world. These ecotypes have come from diverse climates and therefore have faced and adapted to a variety of abiotic and biotic stressors. The sequencing and comparative analysis of these genomes can offer insight into the adaptive strategies of plants. While there are a large number of ecotype genome sequences available, the majority were created using short-read technology. Mapping of short-reads containing structural variation to a reference genome bereft of that variation leads to incorrect mapping of those reads, resulting in a loss of genetic information and introduction of false heterozygosity. For this reason, long-read de novo sequencing of genomes is required to resolve structural variation events. In this article, we sequenced the genomes of eight natural variants of A. thaliana using nanopore sequencing. This resulted in highly contiguous assemblies with >95% of the genome contained within five contigs. The sequencing results from this study include five ecotypes from relict and African populations, an area of untapped genetic diversity. With this study, we increase the knowledge of diversity we have across A. thaliana ecotypes and contribute to ongoing production of an A. thaliana pan-genome.

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.

ModDotPlot-rapid and interactive visualization of tandem repeats

MOTIVATION

A common method for analyzing genomic repeats is to produce a sequence similarity matrix visualized via a dot plot. Innovative approaches such as StainedGlass have improved upon this classic visualization by rendering dot plots as a heatmap of sequence identity, enabling researchers to better visualize multi-megabase tandem repeat arrays within centromeres and other heterochromatic regions of the genome. However, computing the similarity estimates for heatmaps requires high computational overhead and can suffer from decreasing accuracy.

RESULTS

In this work, we introduce ModDotPlot, an interactive and alignment-free dot plot viewer. By approximating average nucleotide identity via a k-mer-based containment index, ModDotPlot produces accurate plots orders of magnitude faster than StainedGlass. We accomplish this through the use of a hierarchical modimizer scheme that can visualize the full 128 Mb genome of Arabidopsis thaliana in under 5 min on a laptop. ModDotPlot is bundled with a graphical user interface supporting real-time interactive navigation of entire chromosomes.

AVAILABILITY AND IMPLEMENTATION

ModDotPlot is available at https://github.com/marbl/ModDotPlot.

Causes and Consequences of Varying Transposable Element Activity: An Evolutionary Perspective

Transposable elements (TEs) are genomic parasites found in nearly all eukaryotes, including humans. This evolutionary success of TEs is due to their replicative activity, involving insertion into new genomic locations. TE activity varies at multiple levels, from between taxa to within individuals. The rapidly accumulating evidence of the influence of TE activity on human health, as well as the rapid growth of new tools to study it, motivated an evaluation of what we know about TE activity thus far. Here, we discuss why TE activity varies, and the consequences of this variation, from an evolutionary perspective. By studying TE activity in nonhuman organisms in the context of evolutionary theories, we can shed light on the factors that affect TE activity. While the consequences of TE activity are usually deleterious, some have lasting evolutionary impacts by conferring benefits on the host or affecting other evolutionary processes.

Characterization of the CsCENH3 protein and centromeric DNA profiles reveal the structures of centromeres in cucumber

Centromeres in eukaryotes mediate the accurate segregation of chromosomes during cell division. They serve as essential functional units of chromosomes and play a core role in the process of genome evolution. Centromeres are composed of satellite repeats and highly repetitive centromeric retrotransposons (CRs), which vary greatly even among closely related species. Cucumber (Cucumis sativus) is a globally cultivated and economically important vegetable and the only species in the Cucumis genus with seven pairs of chromosomes. Therefore, studying the centromeres of the Cucumis subgenus may yield valuable insights into its genome structure and evolution. Using chromatin immunoprecipitation (ChIP) techniques, we isolated centromeric DNA from cucumber reference line 9930. Our investigation into cucumber centromeres uncovered the centromeric satellite sequence, designated as CentCs, and the prevalence of Ty1/Copia long terminal repeat retrotransposons. In addition, active genes were identified in the CsCENH3 nucleosome regions with low transcription levels. To the best of our knowledge, this is the first time that characterization of centromeres has been achieved in cucumber. Meanwhile, our results on the distribution of CentCs and CsCRs in the subgenus Cucumis indicate that the content of centromeric repeats in the wild variants was significantly reduced compared with the cultivated cucumber. The results provide evidence for centromeric DNA amplification that occurred during the domestication process from wild to cultivated cucumber. Furthermore, these findings may offer new information for enhancing our understanding of phylogenetic relationships in the Cucumis genus.

Cotton transposon-related variome reveals roles of transposon-related variations in modern cotton cultivation

INTRODUCTION

Transposon plays a vital role in cotton genome evolution, contributing to the expansion and divergence of genomes within the Gossypium genus. However, knowledge of transposon activity in modern cotton cultivation is limited.

OBJECTIVES

In this study, we aimed to construct transposon-related variome within Gossypium genus and reveal role of transposon-related variations during cotton cultivation. In addition, we try to identify valuable transposon-related variations for cotton breeding.

METHODS

We utilized graphical genome construction to build up the graphical transposon-related variome. Based on the graphical variome, we integrated t-test, eQTL analysis and Mendelian Randomization (MR) to identify valuable transposon activities and elite genes. In addition, a convolutional neural network (CNN) model was constructed to evaluate epigenomic effects of transposon-related variations.

RESULTS

We identified 35,980 transposon activities among 10 cotton genomes, and the diversity of genomic and epigenomic features was observed among 21 transposon categories. The graphical cotton transposon-related variome was constructed, and 9,614 transposon-related variations with plasticity in the modern cotton cohort were used for eQTL, phenotypic t-test and Mendelian Randomization. 128 genes were identified as gene resources improving fiber length and strength simultaneously. 4 genes were selected from 128 genes to construct the elite gene panel whose utility has been validated in a natural cotton cohort and 2 accessions with phenotypic divergence. Based on the eQTL analysis results, we identified transposon-related variations involved in cotton's environmental adaption and human domestication, providing evidence of their role in cotton's adaption-domestication cooperation.

CONCLUSIONS

The cotton transposon-related variome revealed the role of transposon-related variations in modern cotton cultivation, providing genomic resources for cotton molecular breeding.

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 variation and evolution of complete human centromeres

Human centromeres have been traditionally very difficult to sequence and assemble owing to their repetitive nature and large size1. As a result, patterns of human centromeric variation and models for their evolution and function remain incomplete, despite centromeres being among the most rapidly mutating regions2,3. Here, using long-read sequencing, we completely sequenced and assembled all centromeres from a second human genome and compared it to the finished reference genome4,5. We find that the two sets of centromeres show at least a 4.1-fold increase in single-nucleotide variation when compared with their unique flanks and vary up to 3-fold in size. Moreover, we find that 45.8% of centromeric sequence cannot be reliably aligned using standard methods owing to the emergence of new α-satellite higher-order repeats (HORs). DNA methylation and CENP-A chromatin immunoprecipitation experiments show that 26% of the centromeres differ in their kinetochore position by >500 kb. To understand evolutionary change, we selected six chromosomes and sequenced and assembled 31 orthologous centromeres from the common chimpanzee, orangutan and macaque genomes. Comparative analyses reveal a nearly complete turnover of α-satellite HORs, with characteristic idiosyncratic changes in α-satellite HORs for each species. Phylogenetic reconstruction of human haplotypes supports limited to no recombination between the short (p) and long (q) arms across centromeres and reveals that novel α-satellite HORs share a monophyletic origin, providing a strategy to estimate the rate of saltatory amplification and mutation of human centromeric DNA.

A pan-genome of 69 Arabidopsis thaliana accessions reveals a conserved genome structure throughout the global species range

Although originally primarily a system for functional biology, Arabidopsis thaliana has, owing to its broad geographical distribution and adaptation to diverse environments, developed into a powerful model in population genomics. Here we present chromosome-level genome assemblies of 69 accessions from a global species range. We found that genomic colinearity is very conserved, even among geographically and genetically distant accessions. Along chromosome arms, megabase-scale rearrangements are rare and typically present only in a single accession. This indicates that the karyotype is quasi-fixed and that rearrangements in chromosome arms are counter-selected. Centromeric regions display higher structural dynamics, and divergences in core centromeres account for most of the genome size variations. Pan-genome analyses uncovered 32,986 distinct gene families, 60% being present in all accessions and 40% appearing to be dispensable, including 18% private to a single accession, indicating unexplored genic diversity. These 69 new Arabidopsis thaliana genome assemblies will empower future genetic research.

Forces driving transposable element load variation during Arabidopsis range expansion

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

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.

Three near-complete genome assemblies reveal substantial centromere dynamics from diploid to tetraploid in Brachypodium genus

BACKGROUND

Centromeres are critical for maintaining genomic stability in eukaryotes, and their turnover shapes genome architectures and drives karyotype evolution. However, the co-evolution of centromeres from different species in allopolyploids over millions of years remains largely unknown.

RESULTS

Here, we generate three near-complete genome assemblies, a tetraploid Brachypodium hybridum and its two diploid ancestors, Brachypodium distachyon and Brachypodium stacei. We detect high degrees of sequence, structural, and epigenetic variations of centromeres at base-pair resolution between closely related Brachypodium genomes, indicating the appearance and accumulation of species-specific centromere repeats from a common origin during evolution. We also find that centromere homogenization is accompanied by local satellite repeats bursting and retrotransposon purging, and the frequency of retrotransposon invasions drives the degree of interspecies centromere diversification. We further investigate the dynamics of centromeres during alloploidization process, and find that dramatic genetics and epigenetics architecture variations are associated with the turnover of centromeres between homologous chromosomal pairs from diploid to tetraploid. Additionally, our pangenomes analysis reveals the ongoing variations of satellite repeats and stable evolutionary homeostasis within centromeres among individuals of each Brachypodium genome with different polyploidy levels.

CONCLUSIONS

Our results provide unprecedented information on the genomic, epigenomic, and functional diversity of highly repetitive DNA between closely related species and their allopolyploid genomes at both coarse and fine scale.

A centromere map based on super pan-genome highlights the structure and function of rice centromeres

Rice (Oryza sativa) is a significant crop worldwide with a genome shaped by various evolutionary factors. Rice centromeres are crucial for chromosome segregation, and contain some unreported genes. Due to the diverse and complex centromere region, a comprehensive understanding of rice centromere structure and function at the population level is needed. We constructed a high-quality centromere map based on the rice super pan-genome consisting of a 251-accession panel comprising both cultivated and wild species of Asian and African rice. We showed that rice centromeres have diverse satellite repeat CentO, which vary across chromosomes and subpopulations, reflecting their distinct evolutionary patterns. We also revealed that long terminal repeats (LTRs), especially young Gypsy-type LTRs, are abundant in the peripheral CentO-enriched regions and drive rice centromere expansion and evolution. Furthermore, high-quality genome assembly and complete telomere-to-telomere (T2T) reference genome enable us to obtain more centromeric genome information despite mapping and cloning of centromere genes being challenging. We investigated the association between structural variations and gene expression in the rice centromere. A centromere gene, OsMAB, which positively regulates rice tiller number, was further confirmed by expression quantitative trait loci, haplotype analysis and clustered regularly interspaced palindromic repeats (CRISPR)/CRISPR-associated protein 9 methods. By revealing the new insights into the evolutionary patterns and biological roles of rice centromeres, our finding will facilitate future research on centromere biology and crop improvement.

Genomic and cytogenetic analyses reveal satellite repeat signature in allotetraploid okra (Abelmoschus esculentus)

BACKGROUND

Satellite repeats are one of the most rapidly evolving components in eukaryotic genomes and play vital roles in genome regulation, genome evolution, and speciation. As a consequence, the composition, abundance and chromosome distribution of satellite repeats often exhibit variability across various species, genome, and even individual chromosomes. However, we know little about the satellite repeat evolution in allopolyploid genomes.

RESULTS

In this study, we investigated the satellite repeat signature in five okra (Abelmoschus esculentus) accessions using genomic and cytogenetic methods. In each of the five accessions, we identified eight satellite repeats, which exhibited a significant level of intraspecific conservation. Through fluorescence in situ hybridization (FISH) experiments, we observed that the satellite repeats generated multiple signals and exhibited variations in copy number across chromosomes. Intriguingly, we found that five satellite repeats were interspersed with centromeric retrotransposons, signifying their involvement in centromeric satellite repeat identity. We confirmed subgenome-biased amplification patterns of these satellite repeats through existing genome assemblies or dual-color FISH, indicating their distinct dynamic evolution in the allotetraploid okra subgenome. Moreover, we observed the presence of multiple chromosomes harboring the 35 S rDNA loci, alongside another chromosomal pair carrying the 5 S rDNA loci in okra using FISH assay. Remarkably, the intensity of 35 S rDNA hybridization signals varied among chromosomes, with the signals predominantly localized within regions of relatively weak DAPI staining, associated with GC-rich heterochromatin regions. Finally, we observed a similar localization pattern between 35 S rDNA and three satellite repeats with high GC content and confirmed their origin in the intergenic spacer region of the 35 S rDNA.

CONCLUSIONS

Our findings uncover a unique satellite repeat signature in the allotetraploid okra, contributing to our understanding of the composition, abundance, and chromosomal distribution of satellite repeats in allopolyploid genomes, further enriching our understanding of their evolutionary dynamics in complex allopolyploid genomes.

Structural variation and DNA methylation shape the centromere-proximal meiotic crossover landscape in Arabidopsis

BACKGROUND

Centromeres load kinetochore complexes onto chromosomes, which mediate spindle attachment and allow segregation during cell division. Although centromeres perform a conserved cellular function, their underlying DNA sequences are highly divergent within and between species. Despite variability in DNA sequence, centromeres are also universally suppressed for meiotic crossover recombination, across eukaryotes. However, the genetic and epigenetic factors responsible for suppression of centromeric crossovers remain to be completely defined.

RESULTS

To explore the centromere-proximal meiotic recombination landscape, we map 14,397 crossovers against fully assembled Arabidopsis thaliana (A. thaliana) genomes. A. thaliana centromeres comprise megabase satellite repeat arrays that load nucleosomes containing the CENH3 histone variant. Each chromosome contains a structurally polymorphic region of ~3-4 megabases, which lack crossovers and include the satellite arrays. This polymorphic region is flanked by ~1-2 megabase low-recombination zones. These recombination-suppressed regions are enriched for Gypsy/Ty3 retrotransposons, and additionally contain expressed genes with high genetic diversity that initiate meiotic recombination, yet do not crossover. We map crossovers at high-resolution in proximity to CEN3, which resolves punctate centromere-proximal hotspots that overlap gene islands embedded in heterochromatin. Centromeres are densely DNA methylated and the recombination landscape is remodelled in DNA methylation mutants. We observe that the centromeric low-recombining zones decrease and increase crossovers in CG (met1) and non-CG (cmt3) mutants, respectively, whereas the core non-recombining zones remain suppressed.

CONCLUSION

Our work relates the genetic and epigenetic organization of A. thaliana centromeres and flanking pericentromeric heterochromatin to the zones of crossover suppression that surround the CENH3-occupied satellite repeat arrays.

A near-complete assembly of asparagus bean provides insights into anthocyanin accumulation in pods

Asparagus bean (Vigna unguiculata ssp. sesquipedialis), a subspecies of V. unguiculata, is a vital legume crop widely cultivated in Asia for its tender pods consumed as vegetables. However, the existing asparagus bean assemblies still contain numerous gaps and unanchored sequences, which presents challenges to functional genomics research. Here, we present an improved reference genome sequence of an elite asparagus bean variety, Fengchan 6, achieved through the integration of nanopore ultra-long reads, PacBio high-fidelity reads, and Hi-C technology. The improved assembly is 521.3 Mb in length and demonstrates several enhancements, including a higher N50 length (46.4 Mb), an anchor ratio of 99.8%, and the presence of only one gap. Furthermore, we successfully assembled 14 telomeres and all 11 centromeres, including four telomere-to-telomere chromosomes. Remarkably, the centromeric regions cover a total length of 38.1 Mb, providing valuable insights into the complex architecture of centromeres. Among the 30 594 predicted protein-coding genes, we identified 2356 genes that are tandemly duplicated in segmental duplication regions. These findings have implications for defence responses and may contribute to evolutionary processes. By utilizing the reference genome, we were able to effectively identify the presence of the gene VuMYB114, which regulates the accumulation of anthocyanins, thereby controlling the purple coloration of the pods. This discovery holds significant implications for understanding the underlying mechanisms of color determination and the breeding process. Overall, the highly improved reference genome serves as crucial resource and lays a solid foundation for asparagus bean genomic studies and genetic improvement efforts.

Pan-centromere reveals widespread centromere repositioning of soybean genomes

Centromere repositioning refers to a de novo centromere formation at another chromosomal position without sequence rearrangement. This phenomenon was frequently encountered in both mammalian and plant species and has been implicated in genome evolution and speciation. To understand the dynamic of centromeres on soybean genome, we performed the pan-centromere analysis using CENH3-ChIP-seq data from 27 soybean accessions, including 3 wild soybeans, 9 landraces, and 15 cultivars. Building upon the previous discovery of three centromere satellites in soybean, we have identified two additional centromere satellites that specifically associate with chromosome 1. These satellites reveal significant rearrangements in the centromere structures of chromosome 1 across different accessions, consequently impacting the localization of CENH3. By comparative analysis, we reported a high frequency of centromere repositioning on 14 out of 20 chromosomes. Most newly emerging centromeres formed in close proximity to the native centromeres and some newly emerging centromeres were apparently shared in distantly related accessions, suggesting their emergence is independent. Furthermore, we crossed two accessions with mismatched centromeres to investigate how centromere positions would be influenced in hybrid genetic backgrounds. We found that a significant proportion of centromeres in the S9 generation undergo changes in size and position compared to their parental counterparts. Centromeres preferred to locate at satellites to maintain a stable state, highlighting a significant role of centromere satellites in centromere organization. Taken together, these results revealed extensive centromere repositioning in soybean genome and highlighted how important centromere satellites are in constraining centromere positions and supporting centromere function.

Female meiotic drive in plants: mechanisms and dynamics

Female meiosis is fundamentally asymmetric, creating an arena for genetic elements to compete for inclusion in the egg to maximize their transmission. Centromeres, as mediators of chromosomal segregation, are prime candidates to evolve via 'female meiotic drive'. According to the centromere-drive model, the asymmetry of female meiosis ignites a coevolutionary arms race between selfish centromeres and kinetochore proteins, the by-product of which is accelerated sequence divergence. Here, I describe and compare plant models that have been instrumental in uncovering the mechanistic basis of female meiotic drive (maize) and the dynamics of active selfish centromeres in nature (monkeyflowers). Then, I speculate on the mechanistic basis of drive in monkeyflowers, discuss how centromere strength influences chromosomal segregation in plants, and describe new insights into the evolution of plant centromeres.