Alpha-CENTAURI: assessing novel centromeric repeat sequence variation with long read sequencing

Abstract

Motivation: Long arrays of near-identical tandem repeats are a common feature of centromeric and subtelomeric regions in complex genomes. These sequences present a source of repeat structure diversity that is commonly ignored by standard genomic tools. Unlike reads shorter than the underlying repeat structure that rely on indirect inference methods, e.g. assembly, long reads allow direct inference of satellite higher order repeat structure. To automate characterization of local centromeric tandem repeat sequence variation we have designed Alpha-CENTAURI (ALPHA satellite CENTromeric AUtomated Repeat Identification), that takes advantage of Pacific Bioscience long-reads from whole-genome sequencing datasets. By operating on reads prior to assembly, our approach provides a more comprehensive set of repeat-structure variants and is not impacted by rearrangements or sequence underrepresentation due to misassembly.Results: We demonstrate the utility of Alpha-CENTAURI in characterizing repeat structure for alpha satellite containing reads in the hydatidiform mole (CHM1, haploid-like) genome. The pipeline is designed to report local repeat organization summaries for each read, thereby monitoring rearrangements in repeat units, shifts in repeat orientation and sites of array transition into non-satellite DNA, typically defined by transposable element insertion. We validate the method by showing consistency with existing centromere high order repeat references. Alpha-CENTAURI can, in principle, run on any sequence data, offering a method to generate a sequence repeat resolution that could be readily performed using consensus sequences available for other satellite families in genomes without high-quality reference assemblies.Availability and implementation: Documentation and source code for Alpha-CENTAURI are freely available at http://github.com/volkansevim/alpha-CENTAURI.Contact: ali.bashir@mssm.eduSupplementary information: Supplementary data are available at Bioinformatics online.

Higher-order organisation of extremely amplified, potentially functional and massively methylated 5S rDNA in European pikes (Esox sp.)

“…Large palindromes may thus transpose and seed 5S blocks to distal locations. Supporting this, a recent study of a human centromeric satellite (using long read PacBio sequencing) showed increased frequencies of inversions in acrocentric chromosomes compared to other chromosomes [71]. Perhaps, acrocentromeric positions of rDNA could be particularly vulnerable to such rearrangements and/or inversions.…”

centroFlye: Assembling Centromeres with Long Error-Prone Reads

“…Figure 1 shows a dot-plot of the consensus HOR sequence on the X centromere (referred to as DXZ1) formed by 12-monomers with length 2,055. HORs typically occupy multi megabase-sized segments that include rearrangements and transposon insertions (Sevim et al, 2016).…”

New Tools for Hop Cytogenomics: Identification of Tandem Repeat Families from Long-Read Sequences ofHumulus lupulus

“…We and others have successfully mined sequence data to identify tandem repeats that have been developed into FISH probes (Novak et al, 2013;Sevim et al, 2016;Novák et al, 2017;Easterling et al, 2018;Mlinarec et al, 2019). Here we carried out a thorough analysis of PacBio Single Molecule, Real-Time (SMRT) reads (n=1,037,871 reads), each consisting of sequences greater than 5,000 bp long.…”