Semimethylation is a feature of diffuse large B-cell lymphoma, and subgroups with poor prognosis are characterized by global hypomethylation and short telomere length.

BACKGROUND: Large B-cell lymphoma (LBCL) is the most common lymphoma and is known to be a biologically heterogeneous disease regarding genetic, phenotypic, and clinical features. Although the prognosis is good, one-third has a primary refractory or relapsing disease which underscores the importance of developing predictive biological markers capable of identifying high- and low-risk patients. DNA methylation (DNAm) and telomere maintenance alterations are hallmarks of cancer and aging. Both these alterations may contribute to the heterogeneity of the disease, and potentially influence the prognosis of LBCL. RESULTS: We studied the DNAm profiles (Infinium MethylationEPIC BeadChip) and relative telomere lengths (RTL) with qPCR of 93 LBCL cases: Diffuse large B-cell lymphoma not otherwise specified (DLBCL, n = 66), High-grade B-cell lymphoma (n = 7), Primary CNS lymphoma (n = 8), and transformation of indolent B-cell lymphoma (n = 12). There was a substantial methylation heterogeneity in DLBCL and other LBCL entities compared to normal cells and other B-cell neoplasms. LBCL cases had a particularly aberrant semimethylated pattern (0.15 ≤ ß ≤ 0.8) with large intertumor variation and overall low hypermethylation (ß > 0.8). DNAm patterns could not be used to distinguish between germinal center B-cell-like (GC) and non-GC DLBCL cases. In cases treated with R-CHOP-like regimens, a high percentage of global hypomethylation (ß < 0.15) was in multivariable analysis associated with worse disease-specific survival (DSS) (HR 6.920, 95% CI 1.499-31.943) and progression-free survival (PFS) (HR 4.923, 95% CI 1.286-18.849) in DLBCL and with worse DSS (HR 5.147, 95% CI 1.239-21.388) in LBCL. These cases with a high percentage of global hypomethylation also had a higher degree of CpG island methylation, including islands in promoter-associated regions, than the cases with less hypomethylation. Additionally, telomere length was heterogenous in LBCL, with a subset of the DLBCL-GC cases accounting for the longest RTL. Short RTL was independently associated with worse DSS (HR 6.011, 95% CI 1.319-27.397) and PFS (HR 4.689, 95% CI 1.102-19.963) in LBCL treated with R-CHOP-like regimens. CONCLUSION: We hypothesize that subclones with high global hypomethylation and hypermethylated CpG islands could have advantages in tumor progression, e.g. by inactivating tumor suppressor genes or promoting treatment resistance. Our findings suggest that cases with high global hypomethylation and thus poor prognosis could be candidates for alternative treatment regimens including hypomethylating drugs.

Exchange of subtelomeric regions between chromosomes 4q and 10q reverts the FSHD genotype and phenotype.

The most common form of facioscapulohumeral dystrophy (FSHD1) is caused by a partial loss of the D4Z4 macrosatellite repeat array in the subtelomeric region of chromosome 4. Patients with FSHD1 typically carry 1 to 10 D4Z4 repeats, whereas nonaffected individuals have 11 to 150 repeats. The ~150-kilobyte subtelomeric region of the chromosome 10q exhibits a ~99% sequence identity to the 4q, including the D4Z4 array. Nevertheless, contractions of the chr10 array do not cause FSHD or any known disease, as in most people D4Z4 array on chr10 is flanked by the nonfunctional polyadenylation signal, not permitting the DUX4 expression. Here, we attempted to correct the FSHD genotype by a CRISPR-Cas9-induced exchange of the chr4 and chr10 subtelomeric regions. We demonstrated that the induced t(4;10) translocation can generate recombinant genotypes translated into improved FSHD phenotype. FSHD myoblasts with the t(4;10) exhibited reduced expression of the DUX4 targets, restored PAX7 target expression, reduced sensitivity to oxidative stress, and improved differentiation capacity.

Annotations of four high-quality indigenous chicken genomes identify more than one thousand missing genes in subtelomeric regions and micro-chromosomes with high G/C contents.

BACKGROUND: Although multiple chicken genomes have been assembled and annotated, the numbers of protein-coding genes in chicken genomes and their variation among breeds are still uncertain due to the low quality of these genome assemblies and limited resources used in their gene annotations. To fill these gaps, we recently assembled genomes of four indigenous chicken breeds with distinct traits at chromosome-level. In this study, we annotated genes in each of these assembled genomes using a combination of RNA-seq- and homology-based approaches. RESULTS: We identified varying numbers (17,497-17,718) of protein-coding genes in the four indigenous chicken genomes, while recovering 51 of the 274 "missing" genes in birds in general, and 36 of the 174 "missing" genes in chickens in particular. Intriguingly, based on deeply sequenced RNA-seq data collected in multiple tissues in the four breeds, we found 571 ~ 627 protein-coding genes in each genome, which were missing in the annotations of the reference chicken genomes (GRCg6a and GRCg7b/w). After removing redundancy, we ended up with a total of 1,420 newly annotated genes (NAGs). The NAGs tend to be found in subtelomeric regions of macro-chromosomes (chr1 to chr5, plus chrZ) and middle chromosomes (chr6 to chr13, plus chrW), as well as in micro-chromosomes (chr14 to chr39) and unplaced contigs, where G/C contents are high. Moreover, the NAGs have elevated quadruplexes G frequencies, while both G/C contents and quadruplexes G frequencies in their surrounding regions are also high. The NAGs showed tissue-specific expression, and we were able to verify 39 (92.9%) of 42 randomly selected ones in various tissues of the four chicken breeds using RT-qPCR experiments. Most of the NAGs were also encoded in the reference chicken genomes, thus, these genomes might harbor more genes than previously thought. CONCLUSION: The NAGs are widely distributed in wild, indigenous and commercial chickens, and they might play critical roles in chicken physiology. Counting these new genes, chicken genomes harbor more genes than originally thought.

Development of a telomere vector-based approach to overcome limitations caused by lethal phenotypes in the study of essential genes in Magnaporthe oryzae.

Reverse genetic approaches are common tools in genomics for elucidating gene functions, involving techniques such as gene deletion followed by screening for aberrant phenotypes. If the generation of gene deletion mutants fails, the question arises whether the failure stems from technical issues or because the gene of interest (GOI) is essential, meaning that the deletion causes lethality. In this report, we introduce a novel method for assessing gene essentiality using the phytopathogenic ascomycete Magnaporthe oryzae. The method is based on the observation that telomere vectors are lost in transformants during cultivation without selection pressure. We tested the hypothesis that essential genes can be identified in deletion mutants co-transformed with a telomere vector. The M. oryzae gene MoPKC, described in literature as essential, was chosen as GOI. Using CRISPR/Cas9 technology transformants with deleted GOI were generated and backed up by a telomere vector carrying a copy of the GOI and conferring fenhexamid resistance. Transformants in which the GOI deletion in the genome was not successful lost the telomere vector on media without fenhexamid. In contrast, transformants with confirmed GOI deletion retained the telomere vector even in absence of fenhexamid selection. In the latter case, the maintenance of the telomere indicates that the GOI is essential for the surveillance of the fungi, as it would have been lost otherwise. The method presented here allows to test for essentiality of genes when no mutants can be obtained from gene deletion approaches, thereby expanding the toolbox for studying gene function in ascomycetes.

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 SUN-family protein Sad1 mediates heterochromatin spatial organization through interaction with histone H2A-H2B.

Heterochromatin is generally associated with the nuclear periphery, but how the spatial organization of heterochromatin is regulated to ensure epigenetic silencing remains unclear. Here we found that Sad1, an inner nuclear membrane SUN-family protein in fission yeast, interacts with histone H2A-H2B but not H3-H4. We solved the crystal structure of the histone binding motif (HBM) of Sad1 in complex with H2A-H2B, revealing the intimate contacts between Sad1HBM and H2A-H2B. Structure-based mutagenesis studies revealed that the H2A-H2B-binding activity of Sad1 is required for the dynamic distribution of Sad1 throughout the nuclear envelope (NE). The Sad1-H2A-H2B complex mediates tethering telomeres and the mating-type locus to the NE. This complex is also important for heterochromatin silencing. Mechanistically, H2A-H2B enhances the interaction between Sad1 and HDACs, including Clr3 and Sir2, to maintain epigenetic identity of heterochromatin. Interestingly, our results suggest that Sad1 exhibits the histone-enhanced liquid-liquid phase separation property, which helps recruit heterochromatin factors to the NE. Our results uncover an unexpected role of SUN-family proteins in heterochromatin regulation and suggest a nucleosome-independent role of H2A-H2B in regulating Sad1's functionality.

Deregulated DNA ADP-ribosylation impairs telomere replication.

The recognition that DNA can be ADP ribosylated provides an unexpected regulatory level of how ADP-ribosylation contributes to genome stability, epigenetics and immunity. Yet, it remains unknown whether DNA ADP-ribosylation (DNA-ADPr) promotes genome stability and how it is regulated. Here, we show that telomeres are subject to DNA-ADPr catalyzed by PARP1 and removed by TARG1. Mechanistically, we show that DNA-ADPr is coupled to lagging telomere DNA strand synthesis, forming at single-stranded DNA present at unligated Okazaki fragments and on the 3' single-stranded telomere overhang. Persistent DNA-linked ADPr, due to TARG1 deficiency, eventually leads to telomere shortening. Furthermore, using the bacterial DNA ADP-ribosyl-transferase toxin to modify DNA at telomeres directly, we demonstrate that unhydrolyzed DNA-linked ADP-ribose compromises telomere replication and telomere integrity. Thus, by identifying telomeres as chromosomal targets of PARP1 and TARG1-regulated DNA-ADPr, whose deregulation compromises telomere replication and integrity, our study highlights and establishes the critical importance of controlling DNA-ADPr turnover for sustained genome stability.

Association between telomere length and Plasmodium falciparum malaria endemicity in sub-Saharan Africans.

Leukocyte telomere length (LTL) varies significantly across human populations, with individuals of African ancestry having longer LTL than non-Africans. However, the genetic and environmental drivers of LTL variation in Africans remain largely unknown. We report here on the relationship between LTL, genetics, and a variety of environmental and climatic factors in ethnically diverse African adults (n = 1,818) originating from Botswana, Tanzania, Ethiopia, and Cameroon. We observe significant variation in LTL among populations, finding that the San hunter-gatherers from Botswana have the longest leukocyte telomeres and that the Fulani pastoralists from Cameroon have the shortest telomeres. Genetic factors explain ∼50% of LTL variation among individuals. Moreover, we observe a significant negative association between Plasmodium falciparum malaria endemicity and LTL while adjusting for age, sex, and genetics. Within Africa, adults from populations indigenous to areas with high malaria exposure have shorter LTL than those in populations indigenous to areas with low malaria exposure. Finally, we explore to what degree the genetic architecture underlying LTL in Africa covaries with malaria exposure.

Characterization of telomere variant repeats using long reads enables allele-specific telomere length estimation.

Telomeres are regions of repetitive DNA at the ends of linear chromosomes which protect chromosome ends from degradation. Telomere lengths have been extensively studied in the context of aging and disease, though most studies use average telomere lengths which are of limited utility. We present a method for identifying all 92 telomere alleles from long read sequencing data. Individual telomeres are identified using variant repeats proximal to telomere regions, which are unique across alleles. This high-throughput and high-resolution characterization of telomeres could be foundational to future studies investigating the roles of specific telomeres in aging and disease.

Genetic polymorphisms affecting telomere length and their association with cardiovascular disease in the Heinz-Nixdorf-Recall study.

Short telomeres are associated with cardiovascular disease (CVD). We aimed to investigate, if genetically determined telomere-length effects CVD-risk in the Heinz-Nixdorf-Recall study (HNRS) population. We selected 14 single-nucleotide polymorphisms (SNPs) associated with telomere-length (p<10-8) from the literature and after exclusion 9 SNPs were included in the analyses. Additionally, a genetic risk score (GRS) using these 9 SNPs was calculated. Incident CVD was defined as fatal and non-fatal myocardial infarction, stroke, and coronary death. We included 3874 HNRS participants with available genetic data and had no known history of CVD at baseline. Cox proportional-hazards regression was used to test the association between the SNPs/GRS and incident CVD-risk adjusting for common CVD risk-factors. The analyses were further stratified by CVD risk-factors. During follow-up (12.1±4.31 years), 466 participants experienced CVD-events. No association between SNPs/GRS and CVD was observed in the adjusted analyses. However, the GRS, rs10936599, rs2487999 and rs8105767 increase the CVD-risk in current smoker. Few SNPs (rs10936599, rs2487999, and rs7675998) showed an increased CVD-risk, whereas rs10936599, rs677228 and rs4387287 a decreased CVD-risk, in further strata. The results of our study suggest different effects of SNPs/GRS on CVD-risk depending on the CVD risk-factor strata, highlighting the importance of stratified analyses in CVD risk-factors.

Transcription stress at telomeres leads to cytosolic DNA release and paracrine senescence.

Transcription stress has been linked to DNA damage -driven aging, yet the underlying mechanism remains unclear. Here, we demonstrate that Tcea1-/- cells, which harbor a TFIIS defect in transcription elongation, exhibit RNAPII stalling at oxidative DNA damage sites, impaired transcription, accumulation of R-loops, telomere uncapping, chromatin bridges, and genome instability, ultimately resulting in cellular senescence. We found that R-loops at telomeres causally contribute to the release of telomeric DNA fragments in the cytoplasm of Tcea1-/- cells and primary cells derived from naturally aged animals triggering a viral-like immune response. TFIIS-defective cells release extracellular vesicles laden with telomeric DNA fragments that target neighboring cells, which consequently undergo cellular senescence. Thus, transcription stress elicits paracrine signals leading to cellular senescence, promoting aging.

Breaking down causes, consequences, and mediating effects of telomere length variation on human health.

BACKGROUND: Telomeres form repeated DNA sequences at the ends of chromosomes, which shorten with each cell division. Yet, factors modulating telomere attrition and the health consequences thereof are not fully understood. To address this, we leveraged data from 326,363 unrelated UK Biobank participants of European ancestry. RESULTS: Using linear regression and bidirectional univariable and multivariable Mendelian randomization (MR), we elucidate the relationships between leukocyte telomere length (LTL) and 142 complex traits, including diseases, biomarkers, and lifestyle factors. We confirm that telomeres shorten with age and show a stronger decline in males than in females, with these factors contributing to the majority of the 5.4% of LTL variance explained by the phenome. MR reveals 23 traits modulating LTL. Smoking cessation and high educational attainment associate with longer LTL, while weekly alcohol intake, body mass index, urate levels, and female reproductive events, such as childbirth, associate with shorter LTL. We also identify 24 traits affected by LTL, with risk for cardiovascular, pulmonary, and some autoimmune diseases being increased by short LTL, while longer LTL increased risk for other autoimmune conditions and cancers. Through multivariable MR, we show that LTL may partially mediate the impact of educational attainment, body mass index, and female age at childbirth on proxied lifespan. CONCLUSIONS: Our study sheds light on the modulators, consequences, and the mediatory role of telomeres, portraying an intricate relationship between LTL, diseases, lifestyle, and socio-economic factors.

Extrachromosomal telomere DNA derived from excessive strand displacements.

Alternative lengthening of telomeres (ALT) is a telomere maintenance mechanism mediated by break-induced replication, evident in approximately 15% of human cancers. A characteristic feature of ALT cancers is the presence of C-circles, circular single-stranded telomeric DNAs composed of C-rich sequences. Despite the fact that extrachromosomal C-rich single-stranded DNAs (ssDNAs), including C-circles, are unique to ALT cells, their generation process remains undefined. Here, we introduce a method to detect single-stranded telomeric DNA, called 4SET (Strand-Specific Southern-blot for Single-stranded Extrachromosomal Telomeres) assay. Utilizing 4SET, we are able to capture C-rich single-stranded DNAs that are near 200 to 1500 nucleotides in size. Both linear C-rich ssDNAs and C-circles are abundant in the fractions of cytoplasm and nucleoplasm, which supports the idea that linear and circular C-rich ssDNAs are generated concurrently. We also found that C-rich ssDNAs originate during Okazaki fragment processing during lagging strand DNA synthesis. The generation of C-rich ssDNA requires CST-PP (CTC1/STN1/TEN1-PRIMASE-Polymerase alpha) complex-mediated priming of the C-strand DNA synthesis and subsequent excessive strand displacement of the C-rich strand mediated by the DNA Polymerase delta and the BLM helicase. Our work proposes a model for the generation of C-rich ssDNAs and C-circles during ALT-mediated telomere elongation.

Telomere length is longer following diapause in two solitary bee species.

The mechanisms that underlie senescence are not well understood in insects. Telomeres are conserved repetitive sequences at chromosome ends that protect DNA during replication. In many vertebrates, telomeres shorten during cell division and in response to stress and are often used as a cellular marker of senescence. However, little is known about telomere dynamics across the lifespan in invertebrates. We measured telomere length in larvae, prepupae, pupae, and adults of two species of solitary bees, Osmia lignaria and Megachile rotundata. Contrary to our predictions, telomere length was longer in later developmental stages in both O. lignaria and M. rotundata. Longer telomeres occurred after emergence from diapause, which is a physiological state with increased tolerance to stress. In O. lignaria, telomeres were longer in adults when they emerged following diapause. In M. rotundata, telomeres were longer in the pupal stage and subsequent adult stage, which occurs after prepupal diapause. In both species, telomere length did not change during the 8 months of diapause. Telomere length did not differ by mass similarly across species or sex. We also did not see a difference in telomere length after adult O. lignaria were exposed to a nutritional stress, nor did length change during their adult lifespan. Taken together, these results suggest that telomere dynamics in solitary bees differ from what is commonly reported in vertebrates and suggest that insect diapause may influence telomere dynamics.

The impact of telomere length on prostate cancer aggressiveness, genomic instability and health disparities.

The telomere repetitive TTAGGG motif at the ends of chromosomes, serves to preserve genomic integrity and chromosomal stability. In turn, genomic instability is a hallmark of cancer-implicating telomere disturbance. Prostate cancer (PCa) shows significant ancestral disparities, with men of African ancestry at the greatest risk for aggressive disease and associated genomic instability. Yet, no study has explored the role of telomere length (TL) with respect to ancestrally driven PCa health disparities. Patient- and technically-matched tumour-blood whole genome sequencing data for 179 ancestrally defined treatment naïve PCa patients (117 African, 62 European), we assessed for TL (blood and tumour) associations. We found shortened tumour TL to be associated with aggressive PCa presentation and elevated genomic instabilities, including percentage of genome alteration and copy number gains, in men of African ancestry. For European patients, tumour TL showed significant associations with PCa driver genes PTEN, TP53, MSH2, SETBP1 and DDX11L1, while shorter blood TL (< 3200 base pairs) and tumour TL (< 2861 base pairs) were correlated with higher risk for biochemical recurrence. Concurring with previous studies linking TL to PCa diagnosis and/or prognosis, for the first time we correlated TL differences with patient ancestry with important implications for future treatments targeting telomere dysfunction.

PARP2 promotes Break Induced Replication-mediated telomere fragility in response to replication stress.

PARP2 is a DNA-dependent ADP-ribosyl transferase (ARTs) enzyme with Poly(ADP-ribosyl)ation activity that is triggered by DNA breaks. It plays a role in the Base Excision Repair pathway, where it has overlapping functions with PARP1. However, additional roles for PARP2 have emerged in the response of cells to replication stress. In this study, we demonstrate that PARP2 promotes replication stress-induced telomere fragility and prevents telomere loss following chronic induction of oxidative DNA lesions and BLM helicase depletion. Telomere fragility results from the activity of the break-induced replication pathway (BIR). During this process, PARP2 promotes DNA end resection, strand invasion and BIR-dependent mitotic DNA synthesis by orchestrating POLD3 recruitment and activity. Our study has identified a role for PARP2 in the response to replication stress. This finding may lead to the development of therapeutic approaches that target DNA-dependent ART enzymes, particularly in cancer cells with high levels of replication stress.

Reprogramming Chromosome Ends by Functional Histone Acetylation.

Cancers harness embryonic programs to evade aging and promote survival. Normally, sequences at chromosome ends called telomeres shorten with cell division, serving as a countdown clock to limit cell replication. Therefore, a crucial aspect of cancerous transformation is avoiding replicative aging by activation of telomere repair programs. Mouse embryonic stem cells (mESCs) activate a transient expression of the gene Zscan4, which correlates with chromatin de-condensation and telomere extension. Head and neck squamous cell carcinoma (HNSCC) cancers reactivate ZSCAN4, which in turn regulates the phenotype of cancer stem cells (CSCs). Our study reveals a new role for human ZSCAN4 in facilitating functional histone H3 acetylation at telomere chromatin. Next-generation sequencing indicates ZSCAN4 enrichment at telomere chromatin. These changes correlate with ZSCAN4-induced histone H3 acetylation and telomere elongation, while CRISPR/Cas9 knockout of ZSCAN4 leads to reduced H3 acetylation and telomere shortening. Our study elucidates the intricate involvement of ZSCAN4 and its significant contribution to telomere chromatin remodeling. These findings suggest that ZSCAN4 induction serves as a novel link between 'stemness' and telomere maintenance. Targeting ZSCAN4 may offer new therapeutic approaches to effectively limit or enhance the replicative lifespan of stem cells and cancer cells.

Study of association of leptin with leukocyte telomere length in a Chinese rural population.

BACKGROUND: Previous studies have demonstrated the relationship between adipocyte factors, insulin resistance, and other indicators with telomere length. However, these studies did not consider the influence of changes in different indicators on telomere length over time. Therefore, the aim of this study is to elucidate the impact of changes in adipocyte factors, HOMA-IR, and other indicators on the dynamic variation of telomere length. METHODS: The data were from a cohort study conducted in Ningxia, China. A total of 1624 subjects were analyzed. Adipokines and relative leukocyte telomere length (RLTL) were measured, and changes in Homeostatic Model Assessment for Insulin Resistance (HOMA-IR), Homeostatic Model Assessment for ß-Cell Function (HOMA-ß), and Quantitative Insulin Sensitivity Check Index (QUICKI) were calculated. Generalized linear models evaluated associations between changes in adipokines and RLTL changes. Furthermore, univariate analyses examined the effects of changes in adipokines and insulin resistance indicators on ΔRLTL. RESULTS: The research findings indicate that females generally have shorter telomeres compared to males. In comparison to the low-level group of Δleptin (LEP), the high-level group of ΔLEP shows a negative correlation with ΔRLTL (B=-1.32, 95% CI (-2.38, -0.27)). Even after multivariable adjustments, this relationship persists (B=-1.31, 95% CI (-2.24, -0.23)). Further analysis reveals that after adjusting for ΔHOMA-IR, ΔHOMA-ß, and ΔQUICKI, the high-level group of ΔLEP still exhibits a significant negative correlation with ΔRLTL (B=-1.37, 95% CI (-2.43, -0.31)). However, the interaction effects between ΔHOMA-IR, ΔHOMA-ß, ΔQUICKI, and ΔLEP do not affect ΔRLTL. CONCLUSIONS: Elevated levels of leptin were significantly correlated with shortened telomere length. This suggests that increased leptin levels may impact overall individual health by affecting telomere length, underscoring the importance of measures to reduce leptin levels to mitigate the onset and progression of related diseases.

First telomere-to-telomere gapless assembly of the rice blast fungus Pyricularia oryzae.

Rice blast caused by Pyricularia oryzae (syn., Magnaporthe oryzae) was one of the most destructive diseases of rice throughout the world. Genome assembly was fundamental to genetic variation identification and critically impacted the understanding of its ability to overcome host resistance. Here, we report a gapless genome assembly of rice blast fungus P. oryzae strain P131 using PacBio, Illumina and high throughput chromatin conformation capture (Hi-C) sequencing data. This assembly contained seven complete chromosomes (43,237,743 bp) and a circular mitochondrial genome (34,866 bp). Approximately 14.31% of this assembly carried repeat sequences, significantly greater than its previous assembled version. This assembly had a 99.9% complement in BUSCO evaluation. A total of 14,982 genes protein-coding genes were predicted. In summary, we assembled the first telomere-to-telomere gapless genome of P. oryzae, which would be a valuable genome resource for future research on the genome evolution and host adaptation.

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.