The genetic and epigenetic landscape of the Arabidopsis centromeres

Heterochromatin in plant meiosis

Heterochromatin is an organizational property of eukaryotic chromosomes, characterized by extensive DNA and histone modifications, that is associated with the silencing of transposable elements and repetitive sequences. Maintaining heterochromatin is crucial for ensuring genomic integrity and stability during the cell cycle. During meiosis, heterochromatin is important for homologous chromosome synapsis, recombination, and segregation, but our understanding of meiotic heterochromatin formation and condensation is limited. In this review, we focus on the dynamics and features of heterochromatin and how it condenses during meiosis in plants. We also discuss how meiotic heterochromatin influences the interaction and recombination of homologous chromosomes during prophase I.

The plant early recombinosome: a high security complex to break DNA during meiosis

KEY MESSAGE

The formacion of numerous unpredictable DNA Double Strand Breaks (DSBs) on chromosomes iniciates meiotic recombination. In this perspective, we propose a 'multi-key lock' model to secure the risky but necesary breaks as well as a 'one per pair of cromatids' model for the topoisomerase-like early recombinosome. During meiosis, homologous chromosomes recombine at few sites of crossing-overs (COs) to ensure correct segregation. The initiation of meiotic recombination involves the formation of DNA double strand breaks (DSBs) during prophase I. Too many DSBs are dangerous for genome integrity: if these DSBs are not properly repaired, it could potentially lead to chromosomal fragmentation. Too few DSBs are also problematic: if the obligate CO cannot form between bivalents, catastrophic unequal segregation of univalents lead to the formation of sterile aneuploid spores. Research on the regulation of the formation of these necessary but risky DSBs has recently advanced in yeast, mammals and plants. DNA DSBs are created by the enzymatic activity of the early recombinosome, a topoisomerase-like complex containing SPO11. This opinion paper reviews recent insights on the regulation of the SPO11 cofactors necessary for the introduction of temporally and spatially controlled DSBs. We propose that a 'multi-key-lock' model for each subunit of the early recombinosome complex is required to secure the formation of DSBs. We also discuss the hypothetical implications that the established topoisomerase-like nature of the SPO11 core-complex can have in creating DSB in only one of the two replicated chromatids of early prophase I meiotic chromosomes. This hypothetical 'one per pair of chromatids' DSB formation model could optimize the faithful repair of the self-inflicted DSBs. Each DSB could use three potential intact homologous DNA sequences as repair template: one from the sister chromatid and the two others from the homologous chromosomes.

The histone chaperones ASF1 and HIRA are required for telomere length and 45S rDNA copy number homeostasis

Genome stability is significantly influenced by the precise coordination of chromatin complexes that facilitate the loading and eviction of histones from chromatin during replication, transcription, and DNA repair processes. In this study, we investigate the role of the Arabidopsis H3 histone chaperones ANTI-SILENCING FUNCTION 1 (ASF1) and HISTONE REGULATOR A (HIRA) in the maintenance of telomeres and 45S rDNA loci, genomic sites that are particularly susceptible to changes in the chromatin structure. We find that both ASF1 and HIRA are essential for telomere length regulation, as telomeres are significantly shorter in asf1a1b and hira mutants. However, these shorter telomeres remain localized around the nucleolus and exhibit a comparable relative H3 occupancy to the wild type. In addition to regulating telomere length, ASF1 and HIRA contribute to silencing 45S rRNA genes and affect their copy number. Besides, ASF1 supports global heterochromatin maintenance. Our findings also indicate that ASF1 transiently binds to the TELOMERE REPEAT BINDING 1 protein and the N terminus of telomerase in vivo, suggesting a physical link between the ASF1 histone chaperone and the telomere maintenance machinery.

Nanopore ultra-long sequencing and adaptive sampling spur plant complete telomere-to-telomere genome assembly

The pursuit of complete telomere-to-telomere (T2T) genome assembly in plants, challenged by genomic complexity, has been advanced by Oxford Nanopore Technologies (ONT), which offers ultra-long, real-time sequencing. Despite its promise, sequencing length and gap filling remain significant challenges. This study optimized DNA extraction and library preparation, achieving DNA lengths exceeding 485 kb; average N50 read lengths of 80.57 kb, reaching up to 440 kb; and maximum reads of 5.83 Mb. Importantly, we demonstrated that combining ultra-long sequencing and adaptive sampling can effectively fill gaps during assembly, evidenced by successfully filling the remaining gaps of a near-complete Arabidopsis genome assembly and resolving the sequence of an unknown telomeric region in watermelon genome. Collectively, our strategies improve the feasibility of complete T2T genomic assemblies across various plant species, enhancing genome-based research in diverse fields.

Chromosome-level baobab genome illuminates its evolutionary trajectory and environmental adaptation

Baobab (Adansonia digitata) is a long-lived tree endemic to Africa with economic, ecological, and cultural importance, yet its genomic features are underexplored. Here, we report a chromosome-level reference genome anchored to 42 chromosomes for A. digitata, alongside draft assemblies for a sibling tree, two trees from distinct locations in Africa, and A. za from Madagascar. The baobab genome is uniquely rich in DNA transposons, which make up 33%, while LTR retrotransposons account for 10%. A. digitata experienced whole genome multiplication (WGM) around 30 million years ago (MYA), followed by a second WGM event 3-11 MYA, likely linked to autotetraploidy. Resequencing of 25 trees identify three subpopulations, with gene flow across West Africa distinct from East Africa. Gene enrichment and fixation index (Fst) analyses show baobab retained multiple circadian, flowering, and light-responsive genes, which likely support longevity through the UV RESISTANCE LOCUS 8 (UVR8) pathway. In sum, we provide genomic resources and insights for baobab breeding and conservation.

Alternating between even and odd ploidy levels switches on and off the recombination control, even near the centromeres

Meiotic recombination is a key biological process in plant evolution and breeding, as it generates genetic diversity in each generation through the formation of crossovers (COs). However, due to their importance in genome stability, COs are highly regulated in frequency and distribution. We previously demonstrated that this strict regulation of COs can be modified, both in terms of CO frequency and distribution, in allotriploid Brassica hybrids (2n = 3x = 29; AAC) resulting from a cross between Brassica napus (2n = 4x = 38; AACC) and Brassica rapa (2n = 2x = 20; AA). Using the recently updated B. napus genome now including pericentromeres, we demonstrated that COs occur in these cold regions in allotriploids, as close as 375 kb from the centromere. Reverse transcription quantitative PCR (RT-qPCR) of various meiotic genes indicated that Class I COs are likely involved in the increased recombination frequency observed in allotriploids. We also demonstrated that this modified recombination landscape can be maintained via successive generations of allotriploidy (odd ploidy level). This deregulated meiotic behavior reverts to strict regulation in allotetraploid (even ploidy level) progeny in the second generation. Overall, we provide an easy way to manipulate tight recombination control in a polyploid crop.

HEI10 coarsening, chromatin and sequence polymorphism shape the plant meiotic recombination landscape

Meiosis is a conserved eukaryotic cell division that produces spores required for sexual reproduction. During meiosis, chromosomes pair and undergo programmed DNA double-strand breaks, followed by homologous repair that can result in reciprocal crossovers. Crossover formation is highly regulated with typically few events per homolog pair. Crossovers additionally show wider spacing than expected from uniformly random placement - defining the phenomenon of interference. In plants, the conserved HEI10 E3 ligase is initially loaded along meiotic chromosomes, before maturing into a small number of foci, corresponding to crossover locations. We review the coarsening model that explains these dynamics as a diffusion and aggregation process, resulting in approximately evenly spaced HEI10 foci. We review how underlying chromatin states, and the presence of interhomolog polymorphisms, shape the meiotic recombination landscape, in light of the coarsening model. Finally, we consider future directions to understand the control of meiotic recombination in plant genomes.

Multi-omics analyzes of Rosa gigantea illuminate tea scent biosynthesis and release mechanisms

Rose is an important ornamental crop cultivated globally for perfume production. However, our understanding of the mechanisms underlying scent production and molecular breeding for fragrance is hindered by the lack of a reference genome for tea roses. We present the first complete telomere-to-telomere (T2T) genome of Rosa gigantea, with high quality (QV > 60), including detailed characterization of the structural features of repetitive regions. The expansion of genes associated with phenylpropanoid biosynthesis may account for the unique tea scent. We uncover the release rhythm of aromatic volatile organic compounds and their gene regulatory networks through comparative genomics and time-ordered gene co-expression networks. Analyzes of eugenol homologs demonstrate how plants attract pollinators using specialized phenylpropanoids in specific tissues. This study highlights the conservation and utilization of genetic diversity from wild endangered species through multi-omics approaches, providing a scientific foundation for enhancing rose fragrance via de novo domestication.

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.

Atlas of telomeric repeat diversity in Arabidopsis thaliana

BACKGROUND

Telomeric repeat arrays at the ends of chromosomes are highly dynamic in composition, but their repetitive nature and technological limitations have made it difficult to assess their true variation in genome diversity surveys.

RESULTS

We have comprehensively characterized the sequence variation immediately adjacent to the canonical telomeric repeat arrays at the very ends of chromosomes in 74 genetically diverse Arabidopsis thaliana accessions. We first describe several types of distinct telomeric repeat units and then identify evolutionary processes such as local homogenization and higher-order repeat formation that shape diversity of chromosome ends. By comparing largely isogenic samples, we also determine repeat number variation of the degenerate and variant telomeric repeat array at both the germline and somatic levels. Finally, our analysis of haplotype structure uncovers chromosome end-specific patterns in the distribution of variant telomeric repeats, and their linkage to the more proximal non-coding region.

CONCLUSIONS

Our findings illustrate the spectrum of telomeric repeat variation at multiple levels in A. thaliana-in germline and soma, across all chromosome ends, and across genetic groups-thereby expanding our knowledge of the evolution of chromosome ends.

Novel centromeric repetitive DNA elements reveal karyotype dynamics in polyploid sainfoin (Onobrychis viciifolia)

Polyploidy is a common feature in eukaryotes with one of paramount consequences leading to better environmental adaptation. Heterochromatin is often located at telomeres and centromeres and contains repetitive DNA sequences. Sainfoin (Onobrychis viciifolia) is an important perennial forage legume for sustainable agriculture. However, there are only a few studies on the sainfoin genome and chromosomes. In this study, novel tandem repetitive DNA sequences of the sainfoin genome (OnVi180, OnVi169, OnVi176 and OnVidimer) were characterized using bioinformatics, molecular and cytogenetic approaches. The OnVi180 and OnVi169 elements colocalized within functional centromeres. The OnVi176 and OnVidimer elements were localized in centromeric, subtelomeric and interstitial regions. We constructed a sainfoin karyotype that distinguishes the seven basic chromosome groups. Our study provides the first detailed description of heterochromatin and chromosome structure of sainfoin and proposes an origin of heterozygous ancestral genomes, possibly from the same ancestral diploid species, not necessarily from different species, or for chromosome rearrangements after polyploidy. Overall, we discuss our novel and complementary findings in a polyploid crop with unique and complex chromosomal features.

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.

Genome assembly in the telomere-to-telomere era

Genome sequences largely determine the biology and encode the history of an organism, and de novo assembly - the process of reconstructing the genome sequence of an organism from sequencing reads - has been a central problem in bioinformatics for four decades. Until recently, genomes were typically assembled into fragments of a few megabases at best, but now technological advances in long-read sequencing enable the near-complete assembly of each chromosome - also known as telomere-to-telomere assembly - for many organisms. Here, we review recent progress on assembly algorithms and protocols, with a focus on how to derive near-telomere-to-telomere assemblies. We also discuss the additional developments that will be required to resolve remaining assembly gaps and to assemble non-diploid genomes.

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.

Unveiling the distinctive traits of functional rye centromeres: minisatellites, retrotransposons, and R-loop formation

Centromeres play a vital role in cellular division by facilitating kinetochore assembly and spindle attachments. Despite their conserved functionality, centromeric DNA sequences exhibit rapid evolution, presenting diverse sizes and compositions across species. The functional significance of rye centromeric DNA sequences, particularly in centromere identity, remains unclear. In this study, we comprehensively characterized the sequence composition and organization of rye centromeres. Our findings revealed that these centromeres are primarily composed of long terminal repeat retrotransposons (LTR-RTs) and interspersed minisatellites. We systematically classified LTR-RTs into five categories, highlighting the prevalence of younger CRS1, CRS2, and CRS3 of CRSs (centromeric retrotransposons of Secale cereale) were primarily located in the core centromeres and exhibited a higher association with CENH3 nucleosomes. The minisatellites, mainly derived from retrotransposons, along with CRSs, played a pivotal role in establishing functional centromeres in rye. Additionally, we observed the formation of R-loops at specific regions of CRS1, CRS2, and CRS3, with both rye pericentromeres and centromeres exhibiting enrichment in R-loops. Notably, these R-loops selectively formed at binding regions of the CENH3 nucleosome in rye centromeres, suggesting a potential role in mediating the precise loading of CENH3 to centromeres and contributing to centromere specification. Our work provides insights into the DNA sequence composition, distribution, and potential function of R-loops in rye centromeres. This knowledge contributes valuable information to understanding the genetics and epigenetics of rye centromeres, offering implications for the development of synthetic centromeres in future plant modifications and beyond.

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.

Systematic analysis of bZIP gene family in Suaeda australis reveal their roles under salt stress

BACKGROUND

Suaeda australis is one of typical halophyte owing to high levels of salt tolerance. In addition, the bZIP gene family assumes pivotal functions in response to salt stress. However, there are little reports available regarding the bZIP gene family in S. australis.

RESULTS

In this study, we successfully screened 44 bZIP genes within S. australis genome. Subsequently, we conducted an extensive analysis, encompassing investigations into chromosome location, gene structure, phylogenetic relationship, promoter region, conserved motif, and gene expression profile. The 44 bZIP genes were categorized into 12 distinct groups, exhibiting an uneven distribution among the 9 chromosomes of S. australis chromosomes, but one member (Sau23745) was mapped on unanchored scaffolds. Examination of cis-regulatory elements revealed that bZIP promoters were closely related to anaerobic induction, transcription start, and light responsiveness. Comparative transcriptome analysis between ST1 and ST2 samples identified 2,434 DEGs, which were significantly enriched in some primary biological pathways related to salt response-regulating signaling based on GO and KEGG enrichment analysis. Expression patterns analyses clearly discovered the role of several differently expressed SabZIPs, including Sau08107, Sau08911, Sau11415, Sau16575, and Sau19276, which showed higher expression levels in higher salt concentration than low concentration and a response to salt stress. These expression patterns were corroborated through RT-qPCR analysis. The six differentially expressed SabZIP genes, all localized in the nucleus, exhibited positive regulation involved in the salt stress response. SabZIP14, SabZIP26, and SabZIP36 proteins could bind to the promoter region of downstream salt stress-related genes and activate their expressions.

CONCLUSIONS

Our findings offer valuable insights into the evolutionary trajectory of the bZIP gene family in S. australis and shed light on their roles in responding to salt stress. In addition to fundamental genomic information, these results would serve as a foundational framework for future investigations into the regulation of salt stress responses in S. australis.

Karyotypic and phenotypic condensation in allotetraploid wheats accompanied with reproductive strategy transformation: from natural evolution to domestication

MAIN CONCLUSION

Allotetraploid wheat reflects evolutionary divergence and domestication convergence in the karyotypic and phenotypic evolution, accompanied with the transformation from r- strategy to K- strategy in reproductive fitness. Allotetraploid wheat, the progenitor of hexaploidy bread wheat, has undergone 300,000 years of natural evolution and 10,000 years of domestication. The variations in karyotype and phenotype as well as fertility fitness have not been systematically linked. Here, by combining fluorescent in situ hybridization with the quantification of phenotypic and reproductive traits, we compared the karyotype, vegetative growth phenotype and reproductive fitness among synthesized, wild and domesticated accessions of allotetraploid wheat. We detected that the wild accessions showed dramatically high frequencies of homologous recombination and copy number variations of simple sequence repeats (SSR) comparing with synthetic and domesticated accessions. The phenotypic traits reflected significant differences among the populations shaped by distinct evolutionary processes. The diversity observed in wild accessions was significantly greater than that in domesticated ones, particularly in traits associated with vegetative growth and spike morphology. We found that the active pollen of domesticated accessions exhibited greater potential of germination, despite a lower rate of active pollen compared with the wild accessions, indicating a transformation in reproductive fitness strategy for pollen development in domesticated accessions compared to the wild accessions, from r-strategy to K-strategy. Our results demonstrate the condensation of karyotype and phenotype from natural wild accessions to domesticated accessions in allotetraploid wheats. Ecological strategy transformation should be seriously considered from evolution to domestication in polyploid plants, especially crops, which may provide a perspective on the adaptive evolution of polyploid plants.

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.

Multifunctional histone variants in genome function

Histones are integral components of eukaryotic chromatin that have a pivotal role in the organization and function of the genome. The dynamic regulation of chromatin involves the incorporation of histone variants, which can dramatically alter its structural and functional properties. Contrary to an earlier view that limited individual histone variants to specific genomic functions, new insights have revealed that histone variants exert multifaceted roles involving all aspects of genome function, from governing patterns of gene expression at precise genomic loci to participating in genome replication, repair and maintenance. This conceptual change has led to a new understanding of the intricate interplay between chromatin and DNA-dependent processes and how this connection translates into normal and abnormal cellular functions.

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.

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.

Beyond the Human Genome Project: The Age of Complete Human Genome Sequences and Pangenome References

The Human Genome Project was an enormous accomplishment, providing a foundation for countless explorations into the genetics and genomics of the human species. Yet for many years, the human genome reference sequence remained incomplete and lacked representation of human genetic diversity. Recently, two major advances have emerged to address these shortcomings: complete gap-free human genome sequences, such as the one developed by the Telomere-to-Telomere Consortium, and high-quality pangenomes, such as the one developed by the Human Pangenome Reference Consortium. Facilitated by advances in long-read DNA sequencing and genome assembly algorithms, complete human genome sequences resolve regions that have been historically difficult to sequence, including centromeres, telomeres, and segmental duplications. In parallel, pangenomes capture the extensive genetic diversity across populations worldwide. Together, these advances usher in a new era of genomics research, enhancing the accuracy of genomic analysis, paving the path for precision medicine, and contributing to deeper insights into human biology.

The Chaperone NASP Contributes to de Novo Deposition of the Centromeric Histone Variant CENH3 in Arabidopsis Early Embryogenesis

The centromere is an essential chromosome region where the kinetochore is formed to control equal chromosome distribution during cell division. The centromere-specific histone H3 variant CENH3 (also called CENP-A) is a prerequisite for the kinetochore formation. Since CENH3 evolves rapidly, associated factors, including histone chaperones mediating the deposition of CENH3 on the centromere, are thought to act through species-specific amino acid sequences. The functions and interaction networks of CENH3 and histone chaperons have been well-characterized in animals and yeasts. However, molecular mechanisms involved in recognition and deposition of CENH3 are still unclear in plants. Here, we used a swapping strategy between domains of CENH3 of Arabidopsis thaliana and the liverwort Marchantia polymorpha to identify specific regions of CENH3 involved in targeting the centromeres and interacting with the general histone H3 chaperone, nuclear autoantigenic sperm protein (NASP). CENH3's LoopN-α1 region was necessary and sufficient for the centromere targeting in cooperation with the α2 region and was involved in interaction with NASP in cooperation with αN, suggesting a species-specific CENH3 recognition. In addition, by generating an Arabidopsis nasp knock-out mutant in the background of a fully fertile GFP-CENH3/cenh3-1 line, we found that NASP was implicated for de novo CENH3 deposition after fertilization and thus for early embryo development. Our results imply that the NASP mediates the supply of CENH3 in the context of the rapidly evolving centromere identity in land plants.

Genomes of Meniocus linifolius and Tetracme quadricornis reveal the ancestral karyotype and genomic features of core Brassicaceae

Brassicaceae represents an important plant family from both a scientific and economic perspective. However, genomic features related to the early diversification of this family have not been fully characterized, especially upon the uplift of the Tibetan Plateau, which was followed by increasing aridity in the Asian interior, intensifying monsoons in Eastern Asia, and significantly fluctuating daily temperatures. Here, we reveal the genomic architecture that accompanied early Brassicaceae diversification by analyzing two high-quality chromosome-level genomes for Meniocus linifolius (Arabodae; clade D) and Tetracme quadricornis (Hesperodae; clade E), together with genomes representing all major Brassicaceae clades and the basal Aethionemeae. We reconstructed an ancestral core Brassicaceae karyotype (CBK) containing 9 pseudochromosomes with 65 conserved syntenic genomic blocks and identified 9702 conserved genes in Brassicaceae. We detected pervasive conflicting phylogenomic signals accompanied by widespread ancient hybridization events, which correlate well with the early divergence of core Brassicaceae. We identified a successive Brassicaceae-specific expansion of the class I TREHALOSE-6-PHOSPHATE SYNTHASE 1 (TPS1) gene family, which encodes enzymes with essential regulatory roles in flowering time and embryo development. The TPS1s were mainly randomly amplified, followed by expression divergence. Our results provide fresh insights into historical genomic features coupled with Brassicaceae evolution and offer a potential model for broad-scale studies of adaptive radiation under an ever-changing environment.

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

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

A telomere-to-telomere gap-free reference genome assembly of avocado provides useful resources for identifying genes related to fatty acid biosynthesis and disease resistance

Avocado (Persea americana Mill.) is an economically valuable plant because of the high fatty acid content and unique flavor of its fruits. Its fatty acid content, especially the relatively high unsaturated fatty acid content, provides significant health benefits. We herein present a telomere-to-telomere gapless genome assembly (841.6 Mb) of West Indian avocado. The genome contains 40 629 predicted protein-coding genes. Repeat sequences account for 57.9% of the genome. Notably, all telomeres, centromeres, and a nucleolar organizing region are included in this genome. Fragments from these three regions were observed via fluorescence in situ hybridization. We identified 376 potential disease resistance-related nucleotide-binding leucine-rich repeat genes. These genes, which are typically clustered on chromosomes, may be derived from gene duplication events. Five NLR genes (Pa11g0262, Pa02g4855, Pa07g3139, Pa07g0383, and Pa02g3196) were highly expressed in leaves, stems, and fruits, indicating they may be involved in avocado disease responses in multiple tissues. We also identified 128 genes associated with fatty acid biosynthesis and analyzed their expression patterns in leaves, stems, and fruits. Pa02g0113, which encodes one of 11 stearoyl-acyl carrier protein desaturases mediating C18 unsaturated fatty acid synthesis, was more highly expressed in the leaves than in the stems and fruits. These findings provide valuable insights that enhance our understanding of fatty acid biosynthesis in avocado.

Haplotype-resolved assembly of the mule duck genome using high-fidelity sequencing technology

Mule duck is vitally important to the production of global duck meat. Here, we present two high-quality haplotypes of a female mule duck (haplotype 1 (H1):1.28 Gb, haplotype 2 (H2): 1.40 Gb). The continuity (H1: contig N50 = 14.90 Mb, H2: contig N50 = 15.70 Mb) and completeness (BUSCO: H1 = 96.9%, H2 = 97.3%) are substantially better than those of other duck genomes. We detected the structural variations (SVs) in H1 and H2. We observed a positive correlation between autosome length and the number of SVs. Z chromosome was some deficient in deletions and insertions, but W chromosome was some excessive. A total of 1,451 genes were haplotype specific expression (HSEs). Among them, 737 specifically expressed in H1, and 714 specifically expressed in H2. We found that H1 and H2 HSEs tended to be involved in similar biological processes, such as myometrial relaxation and contraction pathways, muscle structure development and phosphorylation. Our haplotype-resolved genome assembly provides a powerful platform for future functional genomics, molecular breeding, and genome editing in mule duck.

Meiotic double-strand break repair DNA synthesis tracts in Arabidopsis thaliana

We report here the successful labelling of meiotic prophase I DNA synthesis in the flowering plant, Arabidopsis thaliana. Incorporation of the thymidine analogue, EdU, enables visualisation of the footprints of recombinational repair of programmed meiotic DNA double-strand breaks (DSB), with ~400 discrete, SPO11-dependent, EdU-labelled chromosomal foci clearly visible at pachytene and later stages of meiosis. This number equates well with previous estimations of 200-300 DNA double-strand breaks per meiosis in Arabidopsis, confirming the power of this approach to detect the repair of most or all SPO11-dependent meiotic DSB repair events. The chromosomal distribution of these DNA-synthesis foci accords with that of early recombination markers and MLH1, which marks Class I crossover sites. Approximately 10 inter-homologue cross-overs (CO) have been shown to occur in each Arabidopsis male meiosis and, athough very probably under-estimated, an equivalent number of inter-homologue gene conversions (GC) have been described. Thus, at least 90% of meiotic recombination events, and very probably more, have not previously been accessible for analysis. Visual examination of the patterns of the foci on the synapsed pachytene chromosomes corresponds well with expectations from the different mechanisms of meiotic recombination and notably, no evidence for long Break-Induced Replication DNA synthesis tracts was found. Labelling of meiotic prophase I, SPO11-dependent DNA synthesis holds great promise for further understanding of the molecular mechanisms of meiotic recombination, at the heart of reproduction and evolution of eukaryotes.

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.

Simultaneous profiling of chromatin accessibility and DNA methylation in complete plant genomes using long-read sequencing

Epigenetic regulations, including chromatin accessibility, nucleosome positioning and DNA methylation intricately shape genome function. However, current chromatin profiling techniques relying on short-read sequencing technologies fail to characterise highly repetitive genomic regions and cannot detect multiple chromatin features simultaneously. Here, we performed Simultaneous Accessibility and DNA Methylation Sequencing (SAM-seq) of purified plant nuclei. Thanks to the use of long-read nanopore sequencing, SAM-seq enables high-resolution profiling of m6A-tagged chromatin accessibility together with endogenous cytosine methylation in plants. Analysis of naked genomic DNA revealed significant sequence preference biases of m6A-MTases, controllable through a normalisation step. By applying SAM-seq to Arabidopsis and maize nuclei we obtained fine-grained accessibility and DNA methylation landscapes genome-wide. We uncovered crosstalk between chromatin accessibility and DNA methylation within nucleosomes of genes, TEs, and centromeric repeats. SAM-seq also detects DNA footprints over cis-regulatory regions. Furthermore, using the single-molecule information provided by SAM-seq we identified extensive cellular heterogeneity at chromatin domains with antagonistic chromatin marks, suggesting that bivalency reflects cell-specific regulations. SAM-seq is a powerful approach to simultaneously study multiple epigenetic features over unique and repetitive sequences, opening new opportunities for the investigation of epigenetic mechanisms.

Emerging roles of plant transcriptional gene silencing under heat

Plants continuously endure unpredictable environmental fluctuations that upset their physiology, with stressful conditions negatively impacting yield and survival. As a contemporary threat of rapid progression, global warming has become one of the most menacing ecological challenges. Thus, understanding how plants integrate and respond to elevated temperatures is crucial for ensuring future crop productivity and furthering our knowledge of historical environmental acclimation and adaptation. While the canonical heat-shock response and thermomorphogenesis have been extensively studied, evidence increasingly highlights the critical role of regulatory epigenetic mechanisms. Among these, the involvement under heat of heterochromatic suppression mediated by transcriptional gene silencing (TGS) remains the least understood. TGS refers to a multilayered metabolic machinery largely responsible for the epigenetic silencing of invasive parasitic nucleic acids and the maintenance of parental imprints. Its molecular effectors include DNA methylation, histone variants and their post-translational modifications, and chromatin packing and remodeling. This work focuses on both established and emerging insights into the contribution of TGS to the physiology of plants under stressful high temperatures. We summarized potential roles of constitutive and facultative heterochromatin as well as the most impactful regulatory genes, highlighting events where the loss of epigenetic suppression has not yet been associated with corresponding changes in epigenetic marks.

Centromere diversity: How different repeat-based holocentromeres may have evolved

In addition to monocentric eukaryotes, which have a single localized centromere on each chromosome, there are holocentric species, with extended repeat-based or repeat-less centromeres distributed over the entire chromosome length. At least two types of repeat-based holocentromeres exist, one composed of many small repeat-based centromere units (small unit-type), and another one characterized by a few large centromere units (large unit-type). We hypothesize that the transposable element-mediated dispersal of hundreds of short satellite arrays formed the small centromere unit-type holocentromere in Rhynchospora pubera. The large centromere unit-type of the plant Chionographis japonica is likely a product of simultaneous DNA double-strand breaks (DSBs), which initiated the de novo formation of repeat-based holocentromeres via insertion of satellite DNA, derived from extra-chromosomal circular DNAs (eccDNAs). The number of initial DSBs along the chromosomes must be higher than the number of centromere units since only a portion of the breaks will have incorporated eccDNA at an appropriate position to serve as future centromere unit sites. Subsequently, preferential incorporation of the centromeric histone H3 variant at these positions is assumed. The identification of repeat-based holocentromeres across lineages will unveil the centromere plasticity and elucidate the mechanisms underlying the diverse formation of holocentromeres.

How diverse a monocentric chromosome can be? Repeatome and centromeric organization of Juncus effusus (Juncaceae)

Juncus is the largest genus of Juncaceae and was considered holocentric for a long time. Recent findings, however, indicated that 11 species from different clades of the genus have monocentric chromosomes. Thus, the Juncus centromere organization and evolution need to be reassessed. We aimed to investigate the major repetitive DNA sequences of two accessions of Juncus effusus and its centromeric structure by employing whole-genome analyses, fluorescent in situ hybridization, CENH3 immunodetection, and chromatin immunoprecipitation sequencing. We showed that the repetitive fraction of the small J. effusus genome (~270 Mbp/1C) is mainly composed of Class I and Class II transposable elements (TEs) and satellite DNAs. Three identified satellite DNA families were mainly (peri)centromeric, with two being associated with the centromeric protein CENH3, but not strictly centromeric. Two types of centromere organization were discerned in J. effusus: type 1 was characterized by a single CENH3 domain enriched with JefSAT1-155 or JefSAT2-180, whereas type 2 showed multiple CENH3 domains interrupted by other satellites, TEs or genes. Furthermore, while type 1 centromeres showed a higher degree of satellite identity along the array, type 2 centromeres had less homogenized arrays along the multiple CENH3 domains per chromosome. Although the analyses confirmed the monocentric organization of J. effusus chromosomes, our data indicate a more dynamic arrangement of J. effusus centromeres than observed for other plant species, suggesting it may constitute a transient state between mono- and holocentricity.

Scalable telomere-to-telomere assembly for diploid and polyploid genomes with double graph

Despite advances in long-read sequencing technologies, constructing a near telomere-to-telomere assembly is still computationally demanding. Here we present hifiasm (UL), an efficient de novo assembly algorithm combining multiple sequencing technologies to scale up population-wide near telomere-to-telomere assemblies. Applied to 22 human and two plant genomes, our algorithm produces better diploid assemblies at a cost of an order of magnitude lower than existing methods, and it also works with polyploid genomes.

Complete genome assembly provides a high-quality skeleton for pan-NLRome construction in melon

Melon (Cucumis melo L.), being under intensive domestication and selective breeding, displays an abundant phenotypic diversity. Wild germplasm with tolerance to stress represents an untapped genetic resource for discovery of disease-resistance genes. To comprehensively characterize resistance genes in melon, we generate a telomere-to-telomere (T2T) and gap-free genome of wild melon accession PI511890 (C. melo var. chito) with a total length of 375.0 Mb and a contig N50 of 31.24 Mb. The complete genome allows us to dissect genome architecture and identify resistance gene analogs. We construct a pan-NLRome using seven melon genomes, which include 208 variable and 18 core nucleotide-binding leucine-rich repeat receptors (NLRs). Multiple disease-related transcriptome analyses indicate that most up-regulated NLRs induced by pathogens are shell or cloud NLRs. The T2T gap-free assembly and the pan-NLRome not only serve as essential resources for genomic studies and molecular breeding of melon but also provide insights into the genome architecture and NLR diversity.

H1 restricts euchromatin-associated methylation pathways from heterochromatic encroachment

Silencing pathways prevent transposable element (TE) proliferation and help to maintain genome integrity through cell division. Silenced genomic regions can be classified as either euchromatic or heterochromatic, and are targeted by genetically separable epigenetic pathways. In plants, the RNA-directed DNA methylation (RdDM) pathway targets mostly euchromatic regions, while CMT DNA methyltransferases are mainly associated with heterochromatin. However, many epigenetic features - including DNA methylation patterning - are largely indistinguishable between these regions, so how the functional separation is maintained is unclear. The linker histone H1 is preferentially localized to heterochromatin and has been proposed to restrict RdDM from encroachment. To test this hypothesis, we followed RdDM genomic localization in an h1 mutant by performing ChIP-seq on the largest subunit, NRPE1, of the central RdDM polymerase, Pol V. Loss of H1 resulted in NRPE1 enrichment predominantly in heterochromatic TEs. Increased NRPE1 binding was associated with increased chromatin accessibility in h1, suggesting that H1 restricts NRPE1 occupancy by compacting chromatin. However, RdDM occupancy did not impact H1 localization, demonstrating that H1 hierarchically restricts RdDM positioning. H1 mutants experience major symmetric (CG and CHG) DNA methylation gains, and by generating an h1/nrpe1 double mutant, we demonstrate these gains are largely independent of RdDM. However, loss of NRPE1 occupancy from a subset of euchromatic regions in h1 corresponded to the loss of methylation in all sequence contexts, while at ectopically bound heterochromatic loci, NRPE1 deposition correlated with increased methylation specifically in the CHH context. Additionally, we found that H1 similarly restricts the occupancy of the methylation reader, SUVH1, and polycomb-mediated H3K27me3. Together, the results support a model whereby H1 helps maintain the exclusivity of heterochromatin by preventing encroachment from other competing pathways.

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.

Combining DNA and protein alignments to improve genome annotation with LiftOn

As the number and variety of assembled genomes continues to grow, the number of annotated genomes is falling behind, particularly for eukaryotes. DNA-based mapping tools help to address this challenge, but they are only able to transfer annotation between closely-related species. Here we introduce LiftOn, a homology-based software tool that integrates DNA and protein alignments to enhance the accuracy of genome-scale annotation and to allow mapping between relatively distant species. LiftOn's protein-centric algorithm considers both types of alignments, chooses optimal open reading frames, resolves overlapping gene loci, and finds additional gene copies where they exist. LiftOn can reliably transfer annotation between genomes representing members of the same species, as we demonstrate on human, mouse, honey bee, rice, and Arabidopsis thaliana. It can further map annotation effectively across species pairs as far apart as mouse and rat or Drosophila melanogaster and D. erecta.

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.

The role of centromeric repeats and transcripts in kinetochore assembly and function

Centromeres are the chromosomal domains, where the kinetochore protein complex is formed, mediating proper segregation of chromosomes during cell division. Although the function of centromeres has remained conserved during evolution, centromeric DNA is highly variable, even in closely related species. In addition, the composition of the kinetochore complexes varies among organisms. Therefore, it is assumed that the centromeric position is determined epigenetically, and the centromeric histone H3 (CENH3) serves as an epigenetic marker. The loading of CENH3 onto centromeres depends on centromere-licensing factors, chaperones, and transcription of centromeric repeats. Several proteins that regulate CENH3 loading and kinetochore assembly interact with the centromeric transcripts and DNA in a sequence-independent manner. However, the functional aspects of these interactions are not fully understood. This review discusses the variability of centromeric sequences in different organisms and the regulation of their transcription through the RNA Pol II and RNAi machinery. The data suggest that the interaction of proteins involved in CENH3 loading and kinetochore assembly with centromeric DNA and transcripts plays a role in centromere, and possibly neocentromere, formation in a sequence-independent manner.

ModDotPlot-Rapid and interactive visualization of complex 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 Mbp genome of Arabidopsis thaliana in under 5 minutes 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.

Near telomere-to-telomere genome assemblies of two Chlorella species unveil the composition and evolution of centromeres in green algae

BACKGROUND

Centromeres play a crucial and conserved role in cell division, although their composition and evolutionary history in green algae, the evolutionary ancestors of land plants, remains largely unknown.

RESULTS

We constructed near telomere-to-telomere (T2T) assemblies for two Trebouxiophyceae species, Chlorella sorokiniana NS4-2 and Chlorella pyrenoidosa DBH, with chromosome numbers of 12 and 13, and genome sizes of 58.11 Mb and 53.41 Mb, respectively. We identified and validated their centromere sequences using CENH3 ChIP-seq and found that, similar to humans and higher plants, the centromeric CENH3 signals of green algae display a pattern of hypomethylation. Interestingly, the centromeres of both species largely comprised transposable elements, although they differed significantly in their composition. Species within the Chlorella genus display a more diverse centromere composition, with major constituents including members of the LTR/Copia, LINE/L1, and LINE/RTEX families. This is in contrast to green algae including Chlamydomonas reinhardtii, Coccomyxa subellipsoidea, and Chromochloris zofingiensis, in which centromere composition instead has a pronounced single-element composition. Moreover, we observed significant differences in the composition and structure of centromeres among chromosomes with strong collinearity within the Chlorella genus, suggesting that centromeric sequence evolves more rapidly than sequence in non-centromeric regions.

CONCLUSIONS

This study not only provides high-quality genome data for comparative genomics of green algae but gives insight into the composition and evolutionary history of centromeres in early plants, laying an important foundation for further research on their evolution.

Targeted A-to-G base editing in the organellar genomes of Arabidopsis with monomeric programmable deaminases

Plastids and mitochondria are 2 intracellular organelles containing DNA-encoding partial but essential components for their roles, photosynthesis, and respiration. Precise base editing in both plastid and mitochondrial genomes would benefit their gene functional analysis and crop breeding. Targeted base editing in organellar genomes relies on a protein-based genome-editing system that uses the TALE-DNA recognition motif with deaminases. This is because the efficient delivery of guide RNA for clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 systems into organelles is currently impossible. Since TALE-based base editors used in organellar genomes are usually dimeric types, in this study, we used targeted A-to-G base editing in Arabidopsis (Arabidopsis thaliana) plastid and mitochondrial genomes with monomeric TALE-based deaminase for easier assembling of vectors. As a result, inheritable targeted A-to-G base editing of adenosine triphosphatase subunit 6-2 (atp6-2) in plant mitochondrial genomes and of 16S ribosomal RNA (16S rRNA) in plastid genomes of Arabidopsis was successfully induced by monomeric TALE-based adenine deaminase (AD) without off-target mutations. The monomeric TALE-based adenine deaminases also demonstrated a preference for editing the 8th T on the same strand from the recognition end. Phenotypic analysis showed that A-to-G conversion at 1139A of plastid 16S rRNA conferred substantial spectinomycin resistance in Arabidopsis, but not the other 2 potential-resistant mutations at 1131T and 1137T, predicted from the previous bacterial data. Our study demonstrated the feasibility of monomeric TALE-based ADs in plant organelles and their potential contribution to the functional analyses of plant organelles with easier assembling.

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.

H1 restricts euchromatin-associated methylation pathways from heterochromatic encroachment

Silencing pathways prevent transposable element (TE) proliferation and help to maintain genome integrity through cell division. Silenced genomic regions can be classified as either euchromatic or heterochromatic, and are targeted by genetically separable epigenetic pathways. In plants, the RNA-directed DNA methylation (RdDM) pathway targets mostly euchromatic regions, while CMT DNA methyltransferases are mainly associated with heterochromatin. However, many epigenetic features - including DNA methylation patterning - are largely indistinguishable between these regions, so how the functional separation is maintained is unclear. The linker histone H1 is preferentially localized to heterochromatin and has been proposed to restrict RdDM from encroachment. To test this hypothesis, we followed RdDM genomic localization in an h1 mutant by performing ChIP-seq on the largest subunit, NRPE1, of the central RdDM polymerase, Pol V. Loss of H1 resulted in NRPE1 enrichment predominantly in heterochromatic TEs. Increased NRPE1 binding was associated with increased chromatin accessibility in h1 , suggesting that H1 restricts NRPE1 occupancy by compacting chromatin. However, RdDM occupancy did not impact H1 localization, demonstrating that H1 hierarchically restricts RdDM positioning. H1 mutants experience major symmetric (CG and CHG) DNA methylation gains, and by generating an h1/nrpe1 double mutant, we demonstrate these gains are largely independent of RdDM. However, loss of NRPE1 occupancy from a subset of euchromatic regions in h1 corresponded to loss of methylation in all sequence contexts, while at ectopically bound heterochromatic loci, NRPE1 deposition correlated with increased methylation specifically in the CHH context. Additionally, we found that H1 similarly restricts the occupancy of the methylation reader, SUVH1, and polycomb-mediated H3K27me3. Together, the results support a model whereby H1 helps maintain the exclusivity of heterochromatin by preventing encroachment from other competing pathways.

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.