Complete genomic and epigenetic maps of human centromeres

The p-Arms of Human Acrocentric Chromosomes Play by a Different Set of Rules

The p-arms of the five human acrocentric chromosomes bear nucleolar organizer regions (NORs) comprising ribosomal gene (rDNA) repeats that are organized in a homogeneous tandem array and transcribed in a telomere-to-centromere direction. Precursor ribosomal RNA transcripts are processed and assembled into ribosomal subunits, the nucleolus being the physical manifestation of this process. I review current understanding of nucleolar chromosome biology and describe current exploration into a role for the NOR chromosomal context. Full DNA sequences for acrocentric p-arms are now emerging, aided by the current revolution in long-read sequencing and genome assembly. Acrocentric p-arms vary from 10.1 to 16.7 Mb, accounting for ∼2.2% of the genome. Bordering rDNA arrays, distal junctions, and proximal junctions are shared among the p-arms, with distal junctions showing evidence of functionality. The remaining p-arm sequences comprise multiple satellite DNA classes and segmental duplications that facilitate recombination between heterologous chromosomes, which is likely also involved in Robertsonian translocations.

Long-Read DNA Sequencing: Recent Advances and Remaining Challenges

DNA sequencing has revolutionized medicine over recent decades. However, analysis of large structural variation and repetitive DNA, a hallmark of human genomes, has been limited by short-read technology, with read lengths of 100-300 bp. Long-read sequencing (LRS) permits routine sequencing of human DNA fragments tens to hundreds of kilobase pairs in size, using both real-time sequencing by synthesis and nanopore-based direct electronic sequencing. LRS permits analysis of large structural variation and haplotypic phasing in human genomes and has enabled the discovery and characterization of rare pathogenic structural variants and repeat expansions. It has also recently enabled the assembly of a complete, gapless human genome that includes previously intractable regions, such as highly repetitive centromeres and homologous acrocentric short arms. With the addition of protocols for targeted enrichment, direct epigenetic DNA modification detection, and long-range chromatin profiling, LRS promises to launch a new era of understanding of genetic diversity and pathogenic mutations in human populations.

Alternative end-joining results in smaller deletions in heterochromatin relative to euchromatin

Pericentromeric heterochromatin is highly enriched for repetitive sequences prone to aberrant recombination. Previous studies showed that homologous recombination (HR) repair is uniquely regulated in this domain to enable 'safe' repair while preventing aberrant recombination. In Drosophila cells, DNA double-strand breaks (DSBs) relocalize to the nuclear periphery through nuclear actin-driven directed motions before recruiting the strand invasion protein Rad51 and completing HR repair. End-joining (EJ) repair also occurs with high frequency in heterochromatin of fly tissues, but how alternative EJ (alt-EJ) pathways operate in heterochromatin remains largely uncharacterized. Here, we induce DSBs in single euchromatic and heterochromatic sites using a new system that combines the DR- white reporter and I-SceI expression in spermatogonia of flies. Using this approach, we detect higher frequency of HR repair in heterochromatin, relative to euchromatin. Further, sequencing of mutagenic repair junctions reveals the preferential use of different EJ pathways across distinct euchromatic and heterochromatic sites. Interestingly, synthesis-dependent microhomology-mediated end joining (SD-MMEJ) appears differentially regulated in the two domains, with a preferential use of motifs close to the cut site in heterochromatin relative to euchromatin, resulting in smaller deletions. Together, these studies establish a new approach to study repair outcomes in fly tissues, and support the conclusion that heterochromatin uses more HR and less mutagenic EJ repair relative to euchromatin.

Genome assembly in the telomere-to-telomere era

De novo assembly is the process of reconstructing the genome sequence of an organism from sequencing reads. Genome sequences are essential to biology, and assembly 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 technological advances in long-read sequencing now enable near complete chromosome-level assembly, also known as telomere-to-telomere assembly, for many organisms. Here we review recent progress on assembly algorithms and protocols. We focus on how to derive near telomere-to-telomere assemblies and discuss potential future developments.

Global effects of identity and aging on the human sperm methylome

BACKGROUND

As the average age of fatherhood increases worldwide, so too does the need for understanding effects of aging in male germline cells. Molecular change, including epigenomic alterations, may impact offspring. Age-associated change to DNA cytosine methylation in the cytosine-guanine (CpG) context is a hallmark of aging tissues, including sperm. Prior studies have led to accurate models that predict a man's age based on specific methylation features in the DNA of sperm, but the relationship between aging and global DNA methylation in sperm remains opaque. Further clarification requires a more complete survey of the methylome with assessment of variability within and between individuals.

RESULTS

We collected sperm methylome data in a longitudinal study of ten healthy fertile men. We used whole-genome bisulfite sequencing of samples collected 10 to 18 years apart from each donor. We found that, overall, variability between donors far exceeds age-associated variation. After controlling for donor identity, we see significant age-dependent genome-wide change to the methylome. Notably, trends of change with age depend on genomic location or annotation, with contrasting signatures that correlate with gene density and proximity to centromeres and promoter regions.

CONCLUSIONS

We uncovered epigenetic signatures that reflect a stable process which begins in early adulthood, progressing steadily through most of the male lifespan, and warrants consideration in any future study of the aging sperm epigenome.

Centromere Plasticity With Evolutionary Conservation and Divergence Uncovered by Wheat 10+ Genomes

Centromeres (CEN) are the chromosomal regions that play a crucial role in maintaining genomic stability. The underlying highly repetitive DNA sequences can evolve quickly in most eukaryotes, and promote karyotype evolution. Despite their variability, it is not fully understood how these widely variable sequences ensure the homeostasis of centromere function. In this study, we investigated the genetics and epigenetics of CEN in a population of wheat lines from global breeding programs. We captured a high degree of sequences, positioning, and epigenetic variations in the large and complex wheat CEN. We found that most CENH3-associated repeats are Cereba element of retrotransposons and exhibit phylogenetic homogenization across different wheat lines, but the less-associated repeat sequences diverge on their own way in each wheat line, implying specific mechanisms for selecting certain repeat types as functional core CEN. Furthermore, we observed that CENH3 nucleosome structures display looser wrapping of DNA termini on complex centromeric repeats, including the repositioned CEN. We also found that strict CENH3 nucleosome positioning and intrinsic DNA features play a role in determining centromere identity among different lines. Specific non-B form DNAs were substantially associated with CENH3 nucleosomes for the repositioned centromeres. These findings suggest that multiple mechanisms were involved in the adaptation of CENH3 nucleosomes that can stabilize CEN. Ultimately, we proposed a remarkable epigenetic plasticity of centromere chromatin within the diverse genomic context, and the high robustness is crucial for maintaining centromere function and genome stability in wheat 10+ lines as a result of past breeding selections.

Conserved and divergent mechanisms of inner kinetochore assembly onto centromeric chromatin

Kinetochores are large protein complexes built on centromeric chromatin that mediate chromosome segregation. The inner kinetochore, or constitutive centromere-associated network (CCAN), assembles onto centromeres defined by centromere protein A (CENP-A) nucleosomes (CENP-ANuc), and acts as a platform for the regulated assembly of the microtubule-binding outer kinetochore. Recent cryo-EM work revealed structural conservation of CCAN, from the repeating human regional centromeres to the point centromere of budding yeast. Centromere recognition is determined mainly through engagement of duplex DNA proximal to the CENP-A nucleosome by a DNA-binding CENP-LN channel located at the core of CCAN. Additional DNA interactions formed by other CCAN modules create an enclosed DNA-binding chamber. This configuration explains how kinetochores maintain their tight grip on centromeric DNA to withstand the forces of chromosome segregation. Defining the higher-order architecture of complete kinetochore assemblies with implications for understanding the 3D organisation of regional centromeres and mechanisms of kinetochore dynamics, including how kinetochores sense and respond to tension, are important future directions.

Taming transposable elements in livestock and poultry: a review of their roles and applications

Livestock and poultry play a significant role in human nutrition by converting agricultural by-products into high-quality proteins. To meet the growing demand for safe animal protein, genetic improvement of livestock must be done sustainably while minimizing negative environmental impacts. Transposable elements (TE) are important components of livestock and poultry genomes, contributing to their genetic diversity, chromatin states, gene regulatory networks, and complex traits of economic value. However, compared to other species, research on TE in livestock and poultry is still in its early stages. In this review, we analyze 72 studies published in the past 20 years, summarize the TE composition in livestock and poultry genomes, and focus on their potential roles in functional genomics. We also discuss bioinformatic tools and strategies for integrating multi-omics data with TE, and explore future directions, feasibility, and challenges of TE research in livestock and poultry. In addition, we suggest strategies to apply TE in basic biological research and animal breeding. Our goal is to provide a new perspective on the importance of TE in livestock and poultry genomes.

DNA satellite and chromatin organization at house mouse centromeres and pericentromeres

Centromeres are essential for faithful chromosome segregation during mitosis and meiosis. However, the organization of satellite DNA and chromatin at mouse centromeres and pericentromeres is poorly understood due to the challenges of sequencing and assembling repetitive genomic regions. Using recently available PacBio long-read sequencing data from the C57BL/6 strain and chromatin profiling, we found that contrary to the previous reports of their highly homogeneous nature, centromeric and pericentromeric satellites display varied sequences and organization. We find that both centromeric minor satellites and pericentromeric major satellites exhibited sequence variations within and between arrays. While most arrays are continuous, a significant fraction is interspersed with non-satellite sequences, including transposable elements. Additionally, we investigated CENP-A and H3K9me3 chromatin organization at centromeres and pericentromeres using Chromatin immunoprecipitation sequencing (ChIP-seq). We found that the occupancy of CENP-A and H3K9me3 chromatin at centromeric and pericentric regions, respectively, is associated with increased sequence abundance and homogeneity at these regions. Furthermore, the transposable elements at centromeric regions are not part of functional centromeres as they lack CENP-A enrichment. Finally, we found that while H3K9me3 nucleosomes display a well-phased organization on major satellite arrays, CENP-A nucleosomes on minor satellite arrays lack phased organization. Interestingly, the homogeneous class of major satellites phase CENP-A and H3K27me3 nucleosomes as well, indicating that the nucleosome phasing is an inherent property of homogeneous major satellites. Overall, our findings reveal that house mouse centromeres and pericentromeres, which were previously thought to be highly homogenous, display significant diversity in satellite sequence, organization, and chromatin structure.

DNA 5-methylcytosine detection and methylation phasing using PacBio circular consensus sequencing

Long single-molecular sequencing technologies, such as PacBio circular consensus sequencing (CCS) and nanopore sequencing, are advantageous in detecting DNA 5-methylcytosine in CpGs (5mCpGs), especially in repetitive genomic regions. However, existing methods for detecting 5mCpGs using PacBio CCS are less accurate and robust. Here, we present ccsmeth, a deep-learning method to detect DNA 5mCpGs using CCS reads. We sequence polymerase-chain-reaction treated and M.SssI-methyltransferase treated DNA of one human sample using PacBio CCS for training ccsmeth. Using long (≥10 Kb) CCS reads, ccsmeth achieves 0.90 accuracy and 0.97 Area Under the Curve on 5mCpG detection at single-molecule resolution. At the genome-wide site level, ccsmeth achieves >0.90 correlations with bisulfite sequencing and nanopore sequencing using only 10× reads. Furthermore, we develop a Nextflow pipeline, ccsmethphase, to detect haplotype-aware methylation using CCS reads, and then sequence a Chinese family trio to validate it. ccsmeth and ccsmethphase can be robust and accurate tools for detecting DNA 5-methylcytosines.

Characterization of large-scale genomic differences in the first complete human genome

BACKGROUND

The first telomere-to-telomere (T2T) human genome assembly (T2T-CHM13) release is a milestone in human genomics. The T2T-CHM13 genome assembly extends our understanding of telomeres, centromeres, segmental duplication, and other complex regions. The current human genome reference (GRCh38) has been widely used in various human genomic studies. However, the large-scale genomic differences between these two important genome assemblies are not characterized in detail yet.

RESULTS

Here, in addition to the previously reported "non-syntenic" regions, we find 67 additional large-scale discrepant regions and precisely categorize them into four structural types with a newly developed website tool called SynPlotter. The discrepant regions (~ 21.6 Mbp) excluding telomeric and centromeric regions are highly structurally polymorphic in humans, where the deletions or duplications are likely associated with various human diseases, such as immune and neurodevelopmental disorders. The analyses of a newly identified discrepant region-the KLRC gene cluster-show that the depletion of KLRC2 by a single-deletion event is associated with natural killer cell differentiation in ~ 20% of humans. Meanwhile, the rapid amino acid replacements observed within KLRC3 are probably a result of natural selection in primate evolution.

CONCLUSION

Our study provides a foundation for understanding the large-scale structural genomic differences between the two crucial human reference genomes, and is thereby important for future human genomics studies.

Cancer aneuploidies are shaped primarily by effects on tumour fitness

Aneuploidies-whole-chromosome or whole-arm imbalances-are the most prevalent alteration in cancer genomes1,2. However, it is still debated whether their prevalence is due to selection or ease of generation as passenger events1,2. Here we developed a method, BISCUT, that identifies loci subject to fitness advantages or disadvantages by interrogating length distributions of telomere- or centromere-bounded copy-number events. These loci were significantly enriched for known cancer driver genes, including genes not detected through analysis of focal copy-number events, and were often lineage specific. BISCUT identified the helicase-encoding gene WRN as a haploinsufficient tumour-suppressor gene on chromosome 8p, which is supported by several lines of evidence. We also formally quantified the role of selection and mechanical biases in driving aneuploidy, finding that rates of arm-level copy-number alterations are most highly correlated with their effects on cellular fitness1,2. These results provide insight into the driving forces behind aneuploidy and its contribution to tumorigenesis.

Efficient Formation of Single-copy Human Artificial Chromosomes

Large DNA assembly methodologies underlie milestone achievements in synthetic prokaryotic and budding yeast chromosomes. While budding yeast control chromosome inheritance through ~125 bp DNA sequence-defined centromeres, mammals and many other eukaryotes use large, epigenetic centromeres. Harnessing centromere epigenetics permits human artificial chromosome (HAC) formation but is not sufficient to avoid rampant multimerization of the initial DNA molecule upon introduction to cells. Here, we describe an approach that efficiently forms single-copy HACs. It employs a ~750 kb construct that is sufficiently large to house the distinct chromatin types present at the inner and outer centromere, obviating the need to multimerize. Delivery to mammalian cells is streamlined by employing yeast spheroplast fusion. These developments permit faithful chromosome engineering in the context of metazoan cells.

Holocentromeres can consist of merely a few megabase-sized satellite arrays

The centromere is the chromosome region where microtubules attach during cell division. In contrast to monocentric chromosomes with one centromere, holocentric species usually distribute hundreds of centromere units along the entire chromatid. We assembled the chromosome-scale reference genome and analyzed the holocentromere and (epi)genome organization of the lilioid Chionographis japonica. Remarkably, each of its holocentric chromatids consists of only 7 to 11 evenly spaced megabase-sized centromere-specific histone H3-positive units. These units contain satellite arrays of 23 and 28 bp-long monomers capable of forming palindromic structures. Like monocentric species, C. japonica forms clustered centromeres in chromocenters at interphase. In addition, the large-scale eu- and heterochromatin arrangement differs between C. japonica and other known holocentric species. Finally, using polymer simulations, we model the formation of prometaphase line-like holocentromeres from interphase centromere clusters. Our findings broaden the knowledge about centromere diversity, showing that holocentricity is not restricted to species with numerous and small centromere units.

Long-read sequencing: An effective method for genetic analysis of CYP21A2 variation in congenital adrenal hyperplasia

BACKGROUND

The sequence similarity between CYP21A2 gene and its inactive pseudogene CYP21A1P, and copy number variation (CNV) caused by unequal crossover, make it challenging to characterize the CYP21A2 gene through traditional methods. This study aimed to evaluate the clinical utility of the long-read sequencing (LRS) method in carrier screening and genetic diagnosis of congenital adrenal hyperplasia (CAH) by comparing the efficiency of the LRS method with the conventional multiplex ligation-dependent probe amplification (MLPA) plus Sanger sequencing approaches in CYP21A2 analysis.

METHODS

In a retrospective study, full sequence analysis of the CYP21A2 and CYP21A1P was performed for three pedigrees through long-range locus-specific PCR followed by LRS based on the Pacific Biosciences (PacBio, California, USA) single-molecule real-time (SMRT) platform, and the results were compared with those obtained from next-generation sequencing (NGS)-based whole exome sequencing (WES) and the traditional methods of MLPA plus Sanger sequencing.

RESULTS

The LRS method successfully identified seven CYP21A2 variants , including three single nucleotide variants (NM_000500.9:c.1451G>C p.(Arg484Pro), c.293-13A/C>G (IVS2-13A/C>G), c.518T>A p.(Ile173Asn)), one 111-bp polynucleotide insertion, one set of 3'URT variants (NM_000500.9:c.*368T>C, c.*390A>G, c.*440C>T, c.*443T>C) and two types of chimeric genes and straightforwardly depicted the inheritance patterns of these variants within families. Moreover, the LRS method enabled us to determine the cis-trans configuration of multiple variants in one assay, without the need to analyze additional family samples. Compared with traditional methods, this LRS method can achieve a precise, comprehensive and intuitive result in the genetic diagnosis of 21-hydroxylase deficiency (21-OHD).

CONCLUSION

The LRS method is comprehensive in CYP21A2 analysis and intuitive in result presentation, which holds substantial promise in clinical application as a crucial tool for carrier screening and genetic diagnosis of CAH.

Rio1 downregulates centromeric RNA levels to promote the timely assembly of structurally fit kinetochores

Kinetochores assemble on centromeres via histone H3 variant CENP-A and low levels of centromere transcripts (cenRNAs). The latter are ensured by the downregulation of RNA polymerase II (RNAPII) activity, and cenRNA turnover by the nuclear exosome. Using S. cerevisiae, we now add protein kinase Rio1 to this scheme. Yeast cenRNAs are produced either as short (median lengths of 231 nt) or long (4458 nt) transcripts, in a 1:1 ratio. Rio1 limits their production by reducing RNAPII accessibility and promotes cenRNA degradation by the 5'-3'exoribonuclease Rat1. Rio1 similarly curtails the concentrations of noncoding pericenRNAs. These exist as short transcripts (225 nt) at levels that are minimally two orders of magnitude higher than the cenRNAs. In yeast depleted of Rio1, cen- and pericenRNAs accumulate, CEN nucleosomes and kinetochores misform, causing chromosome instability. The latter phenotypes are also observed with human cells lacking orthologue RioK1, suggesting that CEN regulation by Rio1/RioK1 is evolutionary conserved.

Molecular mechanisms protecting centromeres from self-sabotage and implications for cancer therapy

Centromeres play a crucial role in DNA segregation by mediating the cohesion and separation of sister chromatids during cell division. Centromere dysfunction, breakage or compromised centromeric integrity can generate aneuploidies and chromosomal instability, which are cellular features associated with cancer initiation and progression. Maintaining centromere integrity is thus essential for genome stability. However, the centromere itself is prone to DNA breaks, likely due to its intrinsically fragile nature. Centromeres are complex genomic loci that are composed of highly repetitive DNA sequences and secondary structures and require the recruitment and homeostasis of a centromere-associated protein network. The molecular mechanisms engaged to preserve centromere inherent structure and respond to centromeric damage are not fully understood and remain a subject of ongoing research. In this article, we provide a review of the currently known factors that contribute to centromeric dysfunction and the molecular mechanisms that mitigate the impact of centromere damage on genome stability. Finally, we discuss the potential therapeutic strategies that could arise from a deeper understanding of the mechanisms preserving centromere integrity.

The variation and evolution of complete human centromeres

We completely sequenced and assembled all centromeres from a second human genome and used two reference sets to benchmark genetic, epigenetic, and evolutionary variation within centromeres from a diversity panel of humans and apes. We find that centromere single-nucleotide variation can increase by up to 4.1-fold relative to other genomic regions, with the caveat that up to 45.8% of centromeric sequence, on average, cannot be reliably aligned with current methods due to the emergence of new α-satellite higher-order repeat (HOR) structures and two to threefold differences in the length of the centromeres. The extent to which this occurs differs depending on the chromosome and haplotype. Comparing the two sets of complete human centromeres, we find that eight harbor distinctly different α-satellite HOR array structures and four contain novel α-satellite HOR variants in high abundance. DNA methylation and CENP-A chromatin immunoprecipitation experiments show that 26% of the centromeres differ in their kinetochore position by at least 500 kbp-a property not readily associated with novel α-satellite HORs. To understand evolutionary change, we selected six chromosomes and sequenced and assembled 31 orthologous centromeres from the common chimpanzee, orangutan, and macaque genomes. Comparative analyses reveal nearly complete turnover of α-satellite HORs, but with idiosyncratic changes in structure characteristic to each species. Phylogenetic reconstruction of human haplotypes supports limited to no recombination between the p- and q-arms of human chromosomes and reveals that novel α-satellite HORs share a monophyletic origin, providing a strategy to estimate the rate of saltatory amplification and mutation of human centromeric DNA.

Fission yeast Srr1 and Skb1 promote isochromosome formation at the centromere

Rad51 maintains genome integrity, whereas Rad52 causes non-canonical homologous recombination leading to gross chromosomal rearrangements (GCRs). Here we find that fission yeast Srr1/Ber1 and Skb1/PRMT5 promote GCRs at centromeres. Genetic and physical analyses show that srr1 and skb1 mutations reduce isochromosome formation mediated by centromere inverted repeats. srr1 increases DNA damage sensitivity in rad51 cells but does not abolish checkpoint response, suggesting that Srr1 promotes Rad51-independent DNA repair. srr1 and rad52 additively, while skb1 and rad52 epistatically reduce GCRs. Unlike srr1 or rad52, skb1 does not increase damage sensitivity. Skb1 regulates cell morphology and cell cycle with Slf1 and Pom1, respectively, but neither Slf1 nor Pom1 causes GCRs. Mutating conserved residues in the arginine methyltransferase domain of Skb1 greatly reduces GCRs. These results suggest that, through arginine methylation, Skb1 forms aberrant DNA structures leading to Rad52-dependent GCRs. This study has uncovered roles for Srr1 and Skb1 in GCRs at centromeres.

Cycles of satellite and transposon evolution in Arabidopsis centromeres

Centromeres are critical for cell division, loading CENH3 or CENPA histone variant nucleosomes, directing kinetochore formation and allowing chromosome segregation^1,2. Despite their conserved function, centromere size and structure are diverse across species. To understand this centromere paradox^3,4, it is necessary to know how centromeric diversity is generated and whether it reflects ancient trans-species variation or, instead, rapid post-speciation divergence. To address these questions, we assembled 346 centromeres from 66 Arabidopsis thaliana and 2 Arabidopsis lyrata accessions, which exhibited a remarkable degree of intra- and inter-species diversity. A. thaliana centromere repeat arrays are embedded in linkage blocks, despite ongoing internal satellite turnover, consistent with roles for unidirectional gene conversion or unequal crossover between sister chromatids in sequence diversification. Additionally, centrophilic ATHILA transposons have recently invaded the satellite arrays. To counter ATHILA invasion, chromosome-specific bursts of satellite homogenization generate higher-order repeats and purge transposons, in line with cycles of repeat evolution. Centromeric sequence changes are even more extreme in comparison between A. thaliana and A. lyrata. Together, our findings identify rapid cycles of transposon invasion and purging through satellite homogenization, which drive centromere evolution and ultimately contribute to speciation.

Targeted assembly of ectopic kinetochores to induce chromosome-specific segmental aneuploidies

Cancer cells display persistent underlying chromosomal instability, with individual tumour types intriguingly exhibiting characteristic subsets of whole, and subchromosomal aneuploidies. Few methods to induce specific aneuploidies will exist, hampering investigation of functional consequences of recurrent aneuploidies, as well as the acute consequences of specific chromosome mis-segregation. We therefore investigated the possibility of sabotaging the mitotic segregation of specific chromosomes using nuclease-dead CRISPR-Cas9 (dCas9) as a cargo carrier to specific genomic loci. We recruited the kinetochore-nucleating domain of centromere protein CENP-T to assemble ectopic kinetochores either near the centromere of chromosome 9, or the telomere of chromosome 1. Ectopic kinetochore assembly led to increased chromosome instability and partial aneuploidy of the target chromosomes, providing the potential to induce specific chromosome mis-segregation events in a range of cell types. We also provide an analysis of putative endogenous repeats that could support ectopic kinetochore formation. Overall, our findings provide new insights into ectopic kinetochore biology and represent an important step towards investigating the role of specific aneuploidy and chromosome mis-segregation events in diseases associated with aneuploidy.

Centromere Innovations Within a Mouse Species

Mammalian centromeres direct faithful genetic inheritance and are typically characterized by regions of highly repetitive and rapidly evolving DNA. We focused on a mouse species, Mus pahari, that we found has evolved to house centromere-specifying CENP-A nucleosomes at the nexus of a satellite repeat that we identified and term π-satellite (π-sat), a small number of recruitment sites for CENP-B, and short stretches of perfect telomere repeats. One M. pahari chromosome, however, houses a radically divergent centromere harboring ~6 Mbp of a homogenized π-sat-related repeat, π-satB, that contains >20,000 functional CENP-B boxes. There, CENP-B abundance drives accumulation of microtubule-binding components of the kinetochore, as well as a microtubule-destabilizing kinesin of the inner centromere. The balance of pro- and anti-microtubule-binding by the new centromere permits it to segregate during cell division with high fidelity alongside the older ones whose sequence creates a markedly different molecular composition.

Defining a core configuration for human centromeres during mitosis

The biorientation of sister chromatids on the mitotic spindle, essential for accurate sister chromatid segregation, relies on critical centromere components including cohesin, the centromere-specific H3 variant CENP-A, and centromeric DNA. Centromeric DNA is highly variable between chromosomes yet must accomplish a similar function. Moreover, how the 50 nm cohesin ring, proposed to encircle sister chromatids, accommodates inter-sister centromeric distances of hundreds of nanometers on the metaphase spindle is a conundrum. Insight into the 3D organization of centromere components would help resolve how centromeres function on the mitotic spindle. We used ChIP-seq and super-resolution microscopy to examine the geometry of essential centromeric components on human chromosomes. ChIP-seq demonstrates that cohesin subunits are depleted in α-satellite arrays where CENP-A nucleosomes and kinetochores assemble. Cohesin is instead enriched at pericentromeric DNA. Structured illumination microscopy of sister centromeres is consistent, revealing a non-overlapping pattern of CENP-A and cohesin. We used single particle averaging of hundreds of mitotic sister chromatids to develop an average centromere model. CENP-A clusters on sister chromatids, connected by α-satellite, are separated by ~562 nm with a perpendicular intervening ~190 nM wide axis of cohesin. Two differently sized α-satellite arrays on chromosome 7 display similar inter-sister CENP-A cluster distance, demonstrating different sized arrays can achieve a common spacing. Our data suggest a working model for a common core configuration of essential centromeric components that includes CENP-A nucleosomes at the outer edge of extensible α-satellite DNA and pericentromeric cohesion. This configuration helps reconcile how centromeres function and serves as a foundation for future studies of additional components required for centromere function.

Disparity in the detection of chromosome 15 centromere in patients of African ancestry with a plasma cell neoplasm

Purpose

Fluorescence in situ hybridization (FISH) is the current gold standard assay that provides information related to risk stratification and therapeutic selection for individuals with plasma cell neoplasms. The differential hybridization of FISH probe sets in association with individuals' genetic ancestry has not been previously reported.

Methods

This retrospective study included 1224 bone marrow samples from individuals who had an abnormal plasma cell proliferative disorder FISH result and a concurrent conventional G-banded chromosome study. DNA from bone marrow samples obtained from the G-banded chromosome study was genotyped, and a biogeographical ancestry prediction was carried out.

Results

Using a cohort of individuals with a plasma cell neoplasm, we identified reduced hybridization of a chromosome 15 centromere FISH probe (D15Z4). Metaphase FISH studies of cells with 2 copies of chromosome 15 demonstrated a failure of the D15Z4 FISH probe to hybridize to one chromosome 15 centromere, revealing a false-positive monosomy 15 FISH result in some individuals. Surprisingly, individuals with a monosomy 15 FISH result had a median African ancestry of 77.2% (95% CI 74.1%-80.3%), compared with a median African ancestry of 2.2% (95% CI 2.0%-2.5%) in the non-monosomy 15 cohort (P value = 9.4 × 10-10). Thus, individuals with African ancestry had an 8.02-fold (95% CI 3.73-17.25) increased probability of having a false-positive monosomy 15 result (P value = 9.92 × 10-8).

Conclusion

This study emphasizes a concern regarding the reliability of diagnostic genomic tools and their application in interpreting genetic testing results in diverse patient populations. We discuss alternative methodologies to better represent different ancestry groups in clinical diagnostic testing.

De novo induction of a DNA-histone H3K9 methylation loop on synthetic human repetitive DNA in cultured tobacco cells

Genetic modifications in plants are crucial tools for fundamental and applied research. Transgene expression usually varies among independent lines or their progeny and is associated with the chromatin structure of the insertion site. Strategies based on understanding how to manipulate the epigenetic state of the inserted gene cassette would help to ensure transgene expression. Here, we report a strategy for chromatin manipulation by the artificial tethering of epigenetic effectors to a synthetic human centromeric repetitive DNA (alphoid DNA) platform in plant Bright-Yellow-2 (BY-2) culture cells. By tethering DNA-methyltransferase (Nicotiana tabacum DRM1), we effectively induced DNA methylation and histone methylation (H3K9me2) on the alphoid DNA platform. Tethering of the Arabidopsis SUVH9, which has been reported to lack histone methyltransferase activity, also induced a similar epigenetic state on the alphoid DNA in BY-2 cells, presumably by activating the RNA-dependent DNA methylation (RdDM) pathway. Our results emphasize that the interplay between DNA and histone methylation mechanisms is intrinsic to plant cells. We also found that once epigenetic modification states were induced by the tethering of either DRM1 or SUVH9, the modification was maintained even when the direct tethering of the effector was inhibited. Our system enables the analysis of more diverse epigenetic effectors and will help to elucidate the chromatin assembly mechanisms of plant cells.

A near-complete genome assembly of Brassica rapa provides new insights into the evolution of centromeres

Brassica rapa comprises many important cultivated vegetables and oil crops. However, Chiifu v3.0, the current B. rapa reference genome, still contains hundreds of gaps. Here, we presented a near-complete genome assembly of B. rapa Chiifu v4.0, which was 424.59 Mb with only two gaps, using Oxford Nanopore Technology (ONT) ultralong-read sequencing and Hi-C technologies. The new assembly contains 12 contigs, with a contig N50 of 38.26 Mb. Eight of the ten chromosomes were entirely reconstructed in a single contig from telomere to telomere. We found that the centromeres were mainly invaded by ALE and CRM long terminal repeats (LTRs). Moreover, there is a high divergence of centromere length and sequence among B. rapa genomes. We further found that centromeres are enriched for Copia invaded at 0.14 MYA on average, while pericentromeres are enriched for Gypsy LTRs invaded at 0.51 MYA on average. These results indicated the different invasion mechanisms of LTRs between the two structures. In addition, a novel repetitive sequence PCR630 was identified in the pericentromeres of B. rapa. Overall, the near-complete genome assembly, B. rapa Chiifu v4.0, offers valuable tools for genomic and genetic studies of Brassica species and provides new insights into the evolution of centromeres.

A draft human pangenome reference

Here the Human Pangenome Reference Consortium presents a first draft of the human pangenome reference. The pangenome contains 47 phased, diploid assemblies from a cohort of genetically diverse individuals1. These assemblies cover more than 99% of the expected sequence in each genome and are more than 99% accurate at the structural and base pair levels. Based on alignments of the assemblies, we generate a draft pangenome that captures known variants and haplotypes and reveals new alleles at structurally complex loci. We also add 119 million base pairs of euchromatic polymorphic sequences and 1,115 gene duplications relative to the existing reference GRCh38. Roughly 90 million of the additional base pairs are derived from structural variation. Using our draft pangenome to analyse short-read data reduced small variant discovery errors by 34% and increased the number of structural variants detected per haplotype by 104% compared with GRCh38-based workflows, which enabled the typing of the vast majority of structural variant alleles per sample.

Increased mutation and gene conversion within human segmental duplications

Single-nucleotide variants (SNVs) in segmental duplications (SDs) have not been systematically assessed because of the limitations of mapping short-read sequencing data1,2. Here we constructed 1:1 unambiguous alignments spanning high-identity SDs across 102 human haplotypes and compared the pattern of SNVs between unique and duplicated regions3,4. We find that human SNVs are elevated 60% in SDs compared to unique regions and estimate that at least 23% of this increase is due to interlocus gene conversion (IGC) with up to 4.3 megabase pairs of SD sequence converted on average per human haplotype. We develop a genome-wide map of IGC donors and acceptors, including 498 acceptor and 454 donor hotspots affecting the exons of about 800 protein-coding genes. These include 171 genes that have 'relocated' on average 1.61 megabase pairs in a subset of human haplotypes. Using a coalescent framework, we show that SD regions are slightly evolutionarily older when compared to unique sequences, probably owing to IGC. SNVs in SDs, however, show a distinct mutational spectrum: a 27.1% increase in transversions that convert cytosine to guanine or the reverse across all triplet contexts and a 7.6% reduction in the frequency of CpG-associated mutations when compared to unique DNA. We reason that these distinct mutational properties help to maintain an overall higher GC content of SD DNA compared to that of unique DNA, probably driven by GC-biased conversion between paralogous sequences5,6.

Inversion polymorphism in a complete human genome assembly

The telomere-to-telomere (T2T) complete human reference has significantly improved our ability to characterize genome structural variation. To understand its impact on inversion polymorphisms, we remapped data from 41 genomes against the T2T reference genome and compared it to the GRCh38 reference. We find a ~ 21% increase in sensitivity improving mapping of 63 inversions on the T2T reference. We identify 26 misorientations within GRCh38 and show that the T2T reference is three times more likely to represent the correct orientation of the major human allele. Analysis of 10 additional samples reveals novel rare inversions at chromosomes 15q25.2, 16p11.2, 16q22.1-23.1, and 22q11.21.

The Low-Copy-Number Satellite DNAs of the Model Beetle Tribolium castaneum

The red flour beetle Tribolium castaneum is an important pest of stored agricultural products and the first beetle whose genome was sequenced. So far, one high-copy-number and ten moderate-copy-number satellite DNAs (satDNAs) have been described in the assembled part of its genome. In this work, we aimed to catalog the entire collection of T. castaneum satDNAs. We resequenced the genome using Illumina technology and predicted potential satDNAs via graph-based sequence clustering. In this way, we discovered 46 novel satDNAs that occupied a total of 2.1% of the genome and were, therefore, considered low-copy-number satellites. Their repeat units, preferentially 140-180 bp and 300-340 bp long, showed a high A + T composition ranging from 59.2 to 80.1%. In the current assembly, we annotated the majority of the low-copy-number satDNAs on one or a few chromosomes, discovering mainly transposable elements in their vicinity. The current assembly also revealed that many of the in silico predicted satDNAs were organized into short arrays not much longer than five consecutive repeats, and some of them also had numerous repeat units scattered throughout the genome. Although 20% of the unassembled genome sequence masked the genuine state, the predominance of scattered repeats for some low-copy satDNAs raises the question of whether these are essentially interspersed repeats that occur in tandem only sporadically, with the potential to be satDNA "seeds".

KaryoCreate: A CRISPR-based technology to study chromosome-specific aneuploidy by targeting human centromeres

Aneuploidy, the presence of chromosome gains or losses, is a hallmark of cancer. Here, we describe KaryoCreate (karyotype CRISPR-engineered aneuploidy technology), a system that enables the generation of chromosome-specific aneuploidies by co-expression of an sgRNA targeting chromosome-specific CENPA-binding ɑ-satellite repeats together with dCas9 fused to mutant KNL1. We design unique and highly specific sgRNAs for 19 of the 24 chromosomes. Expression of these constructs leads to missegregation and induction of gains or losses of the targeted chromosome in cellular progeny, with an average efficiency of 8% for gains and 12% for losses (up to 20%) validated across 10 chromosomes. Using KaryoCreate in colon epithelial cells, we show that chromosome 18q loss, frequent in gastrointestinal cancers, promotes resistance to TGF-β, likely due to synergistic hemizygous deletion of multiple genes. Altogether, we describe an innovative technology to create and study chromosome missegregation and aneuploidy in the context of cancer and beyond.

De novo reconstruction of satellite repeat units from sequence data

Satellite DNA are long tandemly repeating sequences in a genome and may be organized as high-order repeats (HORs). They are enriched in centromeres and are challenging to assemble. Existing algorithms for identifying satellite repeats either require the complete assembly of satellites or only work for simple repeat structures without HORs. Here we describe Satellite Repeat Finder (SRF), a new algorithm for reconstructing satellite repeat units and HORs from accurate reads or assemblies without prior knowledge on repeat structures. Applying SRF to real sequence data, we showed that SRF could reconstruct known satellites in human and well-studied model organisms. We also found satellite repeats are pervasive in various other species, accounting for up to 12% of their genome contents but are often underrepresented in assemblies. With the rapid progress on genome sequencing, SRF will help the annotation of new genomes and the study of satellite DNA evolution even if such repeats are not fully assembled.

Variant calling and benchmarking in an era of complete human genome sequences

Genetic variant calling from DNA sequencing has enabled understanding of germline variation in hundreds of thousands of humans. Sequencing technologies and variant-calling methods have advanced rapidly, routinely providing reliable variant calls in most of the human genome. We describe how advances in long reads, deep learning, de novo assembly and pangenomes have expanded access to variant calls in increasingly challenging, repetitive genomic regions, including medically relevant regions, and how new benchmark sets and benchmarking methods illuminate their strengths and limitations. Finally, we explore the possible future of more complete characterization of human genome variation in light of the recent completion of a telomere-to-telomere human genome reference assembly and human pangenomes, and we consider the innovations needed to benchmark their newly accessible repetitive regions and complex variants.

Differences in the intraspecies copy number variation of Arabidopsis thaliana conserved and nonconserved miRNA genes

MicroRNAs (miRNAs) regulate gene expression by RNA interference mechanism. In plants, miRNA genes (MIRs) which are grouped into conserved families, i.e. they are present among the different plant taxa, are involved in the regulation of many developmental and physiological processes. The roles of the nonconserved MIRs-which are MIRs restricted to one plant family, genus, or even species-are less recognized; however, many of them participate in the responses to biotic and abiotic stresses. Both over- and underproduction of miRNAs may influence various biological processes. Consequently, maintaining intracellular miRNA homeostasis seems to be crucial for the organism. Deletions and duplications in the genomic sequence may alter gene dosage and/or activity. We evaluated the extent of copy number variations (CNVs) among Arabidopsis thaliana (Arabidopsis) MIRs in over 1000 natural accessions, using population-based analysis of the short-read sequencing data. We showed that the conserved MIRs were unlikely to display CNVs and their deletions were extremely rare, whereas nonconserved MIRs presented moderate variation. Transposon-derived MIRs displayed exceptionally high diversity. Conversely, MIRs involved in the epigenetic control of transposons reactivated during development were mostly invariable. MIR overlap with the protein-coding genes also limited their variability. At the expression level, a higher rate of nonvariable, nonconserved miRNAs was detectable in Col-0 leaves, inflorescence, and siliques compared to nonconserved variable miRNAs, although the expression of both groups was much lower than that of the conserved MIRs. Our data indicate that CNV rate of Arabidopsis MIRs is related with their age, function, and genomic localization.

Overexpression of Pericentromeric HSAT2 DNA Increases Expression of EMT Markers in Human Epithelial Cancer Cell Lines

Pericentromeric tandemly repeated DNA of human satellites 1, 2, and 3 (HS1, HS2, and HS3) is actively transcribed in some cells. However, the functionality of the transcription remains obscure. Studies in this area have been hampered by the absence of a gapless genome assembly. The aim of our study was to map a transcript that we have previously described as HS2/HS3 on chromosomes using a newly published gapless genome assembly T2T-CHM13, and create a plasmid overexpressing the transcript to assess the influence of HS2/HS3 transcription on cancer cells. We report here that the sequence of the transcript is tandemly repeated on nine chromosomes (1, 2, 7, 9, 10, 16, 17, 22, and Y). A detailed analysis of its genomic localization and annotation in the T2T-CHM13 assembly revealed that the sequence belonged to HSAT2 (HS2) but not to the HS3 family of tandemly repeated DNA. The transcript was found on both strands of HSAT2 arrays. The overexpression of the HSAT2 transcript increased the transcription of the genes encoding the proteins involved in the epithelial-to-mesenchymal transition, EMT (SNAI1, ZEB1, and SNAI2), and the genes that mark cancer-associated fibroblasts (VIM, COL1A1, COL11A1, and ACTA2) in cancer cell lines A549 and HeLa. Co-transfection of the overexpression plasmid and antisense nucleotides eliminated the transcription of EMT genes observed after HSAT2 overexpression. Antisense oligonucleotides also decreased transcription of the EMT genes induced by tumor growth factor beta 1 (TGFβ1). Thus, our study suggests HSAT2 lncRNA transcribed from the pericentromeric tandemly repeated DNA is involved in EMT regulation in cancer cells.

Gaps and complex structurally variant loci in phased genome assemblies

There has been tremendous progress in phased genome assembly production by combining long-read data with parental information or linked-read data. Nevertheless, a typical phased genome assembly generated by trio-hifiasm still generates more than 140 gaps. We perform a detailed analysis of gaps, assembly breaks, and misorientations from 182 haploid assemblies obtained from a diversity panel of 77 unique human samples. Although trio-based approaches using HiFi are the current gold standard, chromosome-wide phasing accuracy is comparable when using Strand-seq instead of parental data. Importantly, the majority of assembly gaps cluster near the largest and most identical repeats (including segmental duplications [35.4%], satellite DNA [22.3%], or regions enriched in GA/AT-rich DNA [27.4%]). Consequently, 1513 protein-coding genes overlap assembly gaps in at least one haplotype, and 231 are recurrently disrupted or missing from five or more haplotypes. Furthermore, we estimate that 6-7 Mbp of DNA are misorientated per haplotype irrespective of whether trio-free or trio-based approaches are used. Of these misorientations, 81% correspond to bona fide large inversion polymorphisms in the human species, most of which are flanked by large segmental duplications. We also identify large-scale alignment discontinuities consistent with 11.9 Mbp of deletions and 161.4 Mbp of insertions per haploid genome. Although 99% of this variation corresponds to satellite DNA, we identify 230 regions of euchromatic DNA with frequent expansions and contractions, nearly half of which overlap with 197 protein-coding genes. Such variable and incompletely assembled regions are important targets for future algorithmic development and pangenome representation.

HiCAT: a tool for automatic annotation of centromere structure

Significant improvements in long-read sequencing technologies have unlocked complex genomic areas, such as centromeres, in the genome and introduced the centromere annotation problem. Currently, centromeres are annotated in a semi-manual way. Here, we propose HiCAT, a generalizable automatic centromere annotation tool, based on hierarchical tandem repeat mining to facilitate decoding of centromere architecture. We apply HiCAT to simulated datasets, human CHM13-T2T and gapless Arabidopsis thaliana genomes. Our results are generally consistent with previous inferences but also greatly improve annotation continuity and reveal additional fine structures, demonstrating HiCAT's performance and general applicability.

Tigerfish designs oligonucleotide-based in situ hybridization probes targeting intervals of highly repetitive DNA at the scale of genomes

Fluorescent in situ hybridization (FISH) is a powerful method for the targeted visualization of nucleic acids in their native contexts. Recent technological advances have leveraged computationally designed oligonucleotide (oligo) probes to interrogate >100 distinct targets in the same sample, pushing the boundaries of FISH-based assays. However, even in the most highly multiplexed experiments, repetitive DNA regions are typically not included as targets, as the computational design of specific probes against such regions presents significant technical challenges. Consequently, many open questions remain about the organization and function of highly repetitive sequences. Here, we introduce Tigerfish, a software tool for the genome-scale design of oligo probes against repetitive DNA intervals. We showcase Tigerfish by designing a panel of 24 interval-specific repeat probes specific to each of the 24 human chromosomes and imaging this panel on metaphase spreads and in interphase nuclei. Tigerfish extends the powerful toolkit of oligo-based FISH to highly repetitive DNA.

Large haplotypes highlight a complex age structure within the maize pan-genome

The genomes of maize and other eukaryotes contain stable haplotypes in regions of low recombination. These regions, including centromeres, long heterochromatic blocks, and rDNA arrays, have been difficult to analyze with respect to their diversity and origin. Greatly improved genome assemblies are now available that enable comparative genomics over these and other nongenic spaces. Using 26 complete maize genomes, we developed methods to align intergenic sequences while excluding genes and regulatory regions. The centromere haplotypes (cenhaps) extend for megabases on either side of the functional centromere regions and appear as evolutionary strata, with haplotype divergence/coalescence times dating as far back as 450 thousand years ago (kya). Application of the same methods to other low recombination regions (heterochromatic knobs and rDNA) and all intergenic spaces revealed that deep coalescence times are ubiquitous across the maize pan-genome. Divergence estimates vary over a broad timescale with peaks at ∼16 and 300 kya, reflecting a complex history of gene flow among diverging populations and changes in population size associated with domestication. Cenhaps and other long haplotypes provide vivid displays of this ancient diversity.

Epigenetic centromere identity is precisely maintained through DNA replication but is uniquely specified among human cells

Centromere identity is defined and maintained epigenetically by the presence of the histone variant CENP-A. How centromeric CENP-A position is specified and precisely maintained through DNA replication is not fully understood. The recently released Telomere-to-Telomere (T2T) genome assembly containing the first complete human centromere sequences provides a new resource for examining CENP-A position. Mapping CENP-A position in clones of the same cell line to the T2T assembly identified highly similar CENP-A position after multiple cell divisions. In contrast, centromeric CENP-A epialleles were evident at several centromeres of different human cell lines, demonstrating the location of CENP-A enrichment and the site of kinetochore recruitment vary among human cells. Across the cell cycle, CENP-A molecules deposited in G1 phase are maintained in their precise position through DNA replication. Thus, despite CENP-A dilution during DNA replication, CENP-A is precisely reloaded onto the same sequences within the daughter centromeres, maintaining unique centromere identity among human cells.

Evolutionary analysis of a complete chicken genome

Microchromosomes are prevalent in nonmammalian vertebrates [P. D. Waters et al., Proc. Natl. Acad. Sci. U.S.A. 118 (2021)], but a few of them are missing in bird genome assemblies. Here, we present a new chicken reference genome containing all autosomes, a Z and a W chromosome, with all gaps closed except for the W. We identified ten small microchromosomes (termed dot chromosomes) with distinct sequence and epigenetic features, among which six were newly assembled. Those dot chromosomes exhibit extremely high GC content and a high level of DNA methylation and are enriched for housekeeping genes. The pericentromeric heterochromatin of dot chromosomes is disproportionately large and continues to expand with the proliferation of satellite DNA and testis-expressed genes. Our analyses revealed that the 41-bp CNM repeat frequently forms higher-order repeats (HORs) at the centromeres of acrocentric chromosomes. The centromere core regions where the kinetochore attaches often encompass telomeric sequence (TTAGGG)n, and in a one of the dot chromosomes, the centromere core recruits an endogenous retrovirus (ERV). We further demonstrate that the W chromosome shares some common features with dot chromosomes, having large arrays of hypermethylated tandem repeats. Finally, using the complete chicken chromosome models, we reconstructed a fine picture of chordate karyotype evolution, revealing frequent chromosomal fusions before and after vertebrate whole-genome duplications. Our sequence and epigenetic characterization of chicken chromosomes shed insights into the understanding of vertebrate genome evolution and chromosome biology.

Centromeres as universal hotspots of DNA breakage, driving RAD51-mediated recombination during quiescence

Centromeres are essential for chromosome segregation in most animals and plants yet are among the most rapidly evolving genome elements. The mechanisms underlying this paradoxical phenomenon remain enigmatic. Here, we report that human centromeres innately harbor a striking enrichment of DNA breaks within functionally active centromere regions. Establishing a single-cell imaging strategy that enables comparative assessment of DNA breaks at repetitive regions, we show that centromeric DNA breaks are induced not only during active cellular proliferation but also de novo during quiescence. Markedly, centromere DNA breaks in quiescent cells are resolved enzymatically by the evolutionarily conserved RAD51 recombinase, which in turn safeguards the specification of functional centromeres. This study highlights the innate fragility of centromeres, which may have been co-opted over time to reinforce centromere specification while driving rapid evolution. The findings also provide insights into how fragile centromeres are likely to contribute to human disease.

Human Satellite 1A analysis provides evidence of pericentromeric transcription

BACKGROUND

Pericentromeric regions of human chromosomes are composed of tandem-repeated and highly organized sequences named satellite DNAs. Human classical satellite DNAs are classified into three families named HSat1, HSat2, and HSat3, which have historically posed a challenge for the assembly of the human reference genome where they are misrepresented due to their repetitive nature. Although being known for a long time as the most AT-rich fraction of the human genome, classical satellite HSat1A has been disregarded in genomic and transcriptional studies, falling behind other human satellites in terms of functional knowledge. Here, we aim to characterize and provide an understanding on the biological relevance of HSat1A.

RESULTS

The path followed herein trails with HSat1A isolation and cloning, followed by in silico analysis. Monomer copy number and expression data was obtained in a wide variety of human cell lines, with greatly varying profiles in tumoral/non-tumoral samples. HSat1A was mapped in human chromosomes and applied in in situ transcriptional assays. Additionally, it was possible to observe the nuclear organization of HSat1A transcripts and further characterize them by 3' RACE-Seq. Size-varying polyadenylated HSat1A transcripts were detected, which possibly accounts for the intricate regulation of alternative polyadenylation.

CONCLUSION

As far as we know, this work pioneers HSat1A transcription studies. With the emergence of new human genome assemblies, acrocentric pericentromeres are becoming relevant characters in disease and other biological contexts. HSat1A sequences and associated noncoding RNAs will most certainly prove significant in the future of HSat research.

DeepConsensus improves the accuracy of sequences with a gap-aware sequence transformer

Circular consensus sequencing with Pacific Biosciences (PacBio) technology generates long (10-25 kilobases), accurate 'HiFi' reads by combining serial observations of a DNA molecule into a consensus sequence. The standard approach to consensus generation, pbccs, uses a hidden Markov model. We introduce DeepConsensus, which uses an alignment-based loss to train a gap-aware transformer-encoder for sequence correction. Compared to pbccs, DeepConsensus reduces read errors by 42%. This increases the yield of PacBio HiFi reads at Q20 by 9%, at Q30 by 27% and at Q40 by 90%. With two SMRT Cells of HG003, reads from DeepConsensus improve hifiasm assembly contiguity (NG50 4.9 megabases (Mb) to 17.2 Mb), increase gene completeness (94% to 97%), reduce the false gene duplication rate (1.1% to 0.5%), improve assembly base accuracy (Q43 to Q45) and reduce variant-calling errors by 24%. DeepConsensus models could be trained to the general problem of analyzing the alignment of other types of sequences, such as unique molecular identifiers or genome assemblies.

Transgenerational epigenetic control of constitutive heterochromatin, transposons, and centromeres

Epigenetic mechanisms are important not only for development but also for genome stability and chromosome dynamics. The latter types of epigenetic controls can often be transgenerational. Here, we review recent progress in two examples of transgenerational epigenetic control: i) the control of constitutive heterochromatin and transposable elements and ii) epigenetic mechanisms that regulate centromere specification and functions. We also discuss the biological significance of enigmatic associations among centromeres, transposons, and constitutive heterochromatin.

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

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

Haplotype-aware pantranscriptome analyses using spliced pangenome graphs

Pangenomics is emerging as a powerful computational paradigm in bioinformatics. This field uses population-level genome reference structures, typically consisting of a sequence graph, to mitigate reference bias and facilitate analyses that were challenging with previous reference-based methods. In this work, we extend these methods into transcriptomics to analyze sequencing data using the pantranscriptome: a population-level transcriptomic reference. Our toolchain, which consists of additions to the VG toolkit and a standalone tool, RPVG, can construct spliced pangenome graphs, map RNA sequencing data to these graphs, and perform haplotype-aware expression quantification of transcripts in a pantranscriptome. We show that this workflow improves accuracy over state-of-the-art RNA sequencing mapping methods, and that it can efficiently quantify haplotype-specific transcript expression without needing to characterize the haplotypes of a sample beforehand.

Innate Immunity in Cardiovascular Diseases-Identification of Novel Molecular Players and Targets

During the past few years, unexpected developments have driven studies in the field of clinical immunology. One driver of immense impact was the outbreak of a pandemic caused by the novel virus SARS-CoV-2. Excellent recent reviews address diverse aspects of immunological re-search into cardiovascular diseases. Here, we specifically focus on selected studies taking advantage of advanced state-of-the-art molecular genetic methods ranging from genome-wide epi/transcriptome mapping and variant scanning to optogenetics and chemogenetics. First, we discuss the emerging clinical relevance of advanced diagnostics for cardiovascular diseases, including those associated with COVID-19-with a focus on the role of inflammation in cardiomyopathies and arrhythmias. Second, we consider newly identified immunological interactions at organ and system levels which affect cardiovascular pathogenesis. Thus, studies into immune influences arising from the intestinal system are moving towards therapeutic exploitation. Further, powerful new research tools have enabled novel insight into brain-immune system interactions at unprecedented resolution. This latter line of investigation emphasizes the strength of influence of emotional stress-acting through defined brain regions-upon viral and cardiovascular disorders. Several challenges need to be overcome before the full impact of these far-reaching new findings will hit the clinical arena.

AlcoR: alignment-free simulation, mapping, and visualization of low-complexity regions in biological data

BACKGROUND

Low-complexity data analysis is the area that addresses the search and quantification of regions in sequences of elements that contain low-complexity or repetitive elements. For example, these can be tandem repeats, inverted repeats, homopolymer tails, GC-biased regions, similar genes, and hairpins, among many others. Identifying these regions is crucial because of their association with regulatory and structural characteristics. Moreover, their identification provides positional and quantity information where standard assembly methodologies face significant difficulties because of substantial higher depth coverage (mountains), ambiguous read mapping, or where sequencing or reconstruction defects may occur. However, the capability to distinguish low-complexity regions (LCRs) in genomic and proteomic sequences is a challenge that depends on the model's ability to find them automatically. Low-complexity patterns can be implicit through specific or combined sources, such as algorithmic or probabilistic, and recurring to different spatial distances-namely, local, medium, or distant associations.

FINDINGS

This article addresses the challenge of automatically modeling and distinguishing LCRs, providing a new method and tool (AlcoR) for efficient and accurate segmentation and visualization of these regions in genomic and proteomic sequences. The method enables the use of models with different memories, providing the ability to distinguish local from distant low-complexity patterns. The method is reference and alignment free, providing additional methodologies for testing, including a highly flexible simulation method for generating biological sequences (DNA or protein) with different complexity levels, sequence masking, and a visualization tool for automatic computation of the LCR maps into an ideogram style. We provide illustrative demonstrations using synthetic, nearly synthetic, and natural sequences showing the high efficiency and accuracy of AlcoR. As large-scale results, we use AlcoR to unprecedentedly provide a whole-chromosome low-complexity map of a recent complete human genome and the haplotype-resolved chromosome pairs of a heterozygous diploid African cassava cultivar.

CONCLUSIONS

The AlcoR method provides the ability of fast sequence characterization through data complexity analysis, ideally for scenarios entangling the presence of new or unknown sequences. AlcoR is implemented in C language using multithreading to increase the computational speed, is flexible for multiple applications, and does not contain external dependencies. The tool accepts any sequence in FASTA format. The source code is freely provided at https://github.com/cobilab/alcor.

The Dynamic Structure and Rapid Evolution of Human Centromeric Satellite DNA

The complete sequence of a human genome provided our first comprehensive view of the organization of satellite DNA associated with heterochromatin. We review how our understanding of the genetic architecture and epigenetic properties of human centromeric DNA have advanced as a result. Preliminary studies of human and nonhuman ape centromeres reveal complex, saltatory mutational changes organized around distinct evolutionary layers. Pockets of regional hypomethylation within higher-order α-satellite DNA, termed centromere dip regions, appear to define the site of kinetochore attachment in all human chromosomes, although such epigenetic features can vary even within the same chromosome. Sequence resolution of satellite DNA is providing new insights into centromeric function with potential implications for improving our understanding of human biology and health.