GraphChainer: Co-linear Chaining for Accurate Alignment of Long Reads to Variation Graphs

Aligning reads to a variation graph is a standard task in pangenomics, with downstream applications in e.g., improving variant calling. While the vg toolkit (Garrison et al., Nature Biotechnology, 2018) is a popular aligner of short reads, GraphAligner (Rautiainen and Marschall, Genome Biology, 2020) is the state-of-the-art aligner of long reads. GraphAligner works by finding candidate read occurrences based on individually extending the best seeds of the read in the variation graph. However, a more principled approach recognized in the community is to co-linearly chain multiple seeds. We present a new algorithm to co-linearly chain a set of seeds in an acyclic variation graph, together with the first efficient implementation of such a co-linear chaining algorithm into a new aligner of long reads to variation graphs, GraphChainer. Compared to GraphAligner, at a normalized edit distance threshold of 40%, it aligns 9% to 12% more reads, and 15% to 19% more total read length, on real PacBio reads from human chromosomes 1 and 22. On both simulated and real data, GraphChainer aligns between 97% and 99% of all reads, and of total read length. At the more stringent normalized edit distance threshold of 30%, GraphChainer aligns up to 29% more total real read length than GraphAligner. GraphChainer is freely available at https://github.com/algbio/GraphChainer

Theory and simulations of condensin mediated loop extrusion in DNA

Condensation of hundreds of mega-base-pair-long human chromosomes in a small nuclear volume is a spectacular biological phenomenon. This process is driven by the formation of chromosome loops. The ATP consuming motor, condensin, interacts with chromatin segments to actively extrude loops. Motivated by real-time imaging of loop extrusion (LE), we created an analytically solvable model, predicting the LE velocity and step size distribution as a function of external load. The theory fits the available experimental data quantitatively, and suggests that condensin must undergo a large conformational change, induced by ATP binding, bringing distant parts of the motor to proximity. Simulations using a simple model confirm that the motor transitions between an open and a closed state in order to extrude loops by a scrunching mechanism, similar to that proposed in DNA bubble formation during bacterial transcription. Changes in the orientation of the motor domains are transmitted over ~50 nm, connecting the motor head and the hinge, thus providing an allosteric basis for LE.

Imputation-free reconstructions of three-dimensional chromosome architectures in human diploid single-cells using allele-specified contacts

The sparseness of chromosomal contact information and the presence of homologous chromosomes with very similar nucleotide sequences make Hi-C analysis difficult. We propose a new algorithm using allele-specific single-nucleotide variations (SNVs) to reconstruct the three-dimensional (3D) chromosomal architectures from the Hi-C dataset of single diploid cells. Our algorithm has a function to discriminate SNVs specifically found between homologous chromosomes to our “recurrence plot”-based algorithm to estimate the 3D chromosome structure, which does not require imputation for ambiguous segment information. The new algorithm can efficiently reconstruct 3D chromosomal structures in single human diploid cells by employing only Hi-C segment pairs containing allele-specific SNVs. The datasets of the remaining pairs of segments without allele-specific SNVs are used to validate the estimated chromosome structure. This approach was used to reconstruct the 3D structures of human chromosomes in single diploid cells at a 1-Mb resolution. Introducing a subsequent mathematical measure further improved the resolution to 40-kb or 100-kb. The reconstruction data reveals that human chromosomes form chromosomal territories and take fractal structures where the mean dimension is a non-integer value. We also validate our approach by estimating 3D protein/polymer structures.

About the origin of the acrocentric part of non-acrocentric satellited chromosomes in humans

Variants in size of the acrocentric short arms (acro-ps) are normally not reported and considered as chromosomal heteromorphisms (CHMs) without any influence on the carrier’s phenotype. However, if acro-ps are translocated to ends of A-chromosomes (i.e. human chromosomes 1-22 and X or Y), those rearrangements are studied in more detail. The aim of the study: Here we characterized 11 healthy carriers of a non-acrocentric satellited chromosomes der(A)t(A;acro)(pter or qter;p1?1.2) to determine the frequency of chromosome 15p and 22p in such rearrangements. Materials and methods: 11 carriers of one (10 cases) or two (1 case) der(A)t(A;acro) were identified during routine cytogenetic analyses. They were originally referred due to infertility or due to a mentally retarded child with otherwise abnormal karyotype. Here derivative chromosomes were studied by fluorescence in situ hybridization applying probes D15Z1 (specific for 15p11.2) and D22Z4 (specific for 22p11.2). As there are no DNA-sequences available for 13p11.2, 14p11.2 and 21p11.2 these regions could not be tested. Results: D15Z1 sequences were identified in 1 out of 12 derivatives der(A)t(A;acro). D22Z1 could not be detected in any of the 11 remainder derivatives. However, only 3 of the 12 der(A)t(A;acro) had acro-ps large enough to potentially comprise sub-band p11.2. Conclusion: In contrast to der(Y)t(Y;acro)(q12;p1?1.2), where in at least 65% of the cases the acro-p part contains D15Z1 sequences, here it could be shown that in der(A)t(A;acro) 15p involvement can be substantiated much less frequently. Also, in none of the two groups D22Z4-sequences were detected in acro-p-parts yet. Besides, breakpoint of acro-pparts in der(A)t(A;acro) seem to be in ~75% of the cases distal from p11.2.

MONI: A Pangenomics Index for Finding MEMs

Recently, Gagie et al. proposed a version of the FM-index, called the r-index, that can store thousands of human genomes on a commodity computer. Then Kuhnle et al. showed how to build the r-index efficiently via a technique called prefix-free parsing (PFP) and demonstrated its effectiveness for exact pattern matching. Exact pattern matching can be leveraged to support approximate pattern matching but the r-index itself cannot support efficiently popular and important queries such as finding maximal exact matches (MEMs). To address this shortcoming, Bannai et al. introduced the concept of thresholds, and showed that storing them together with the r-index enables efficient MEM finding --- but they did not say how to find those thresholds. We present a novel algorithm that applies PFP to build the r-index and find the thresholds simultaneously and in linear time and space with respect to the size of the prefix-free parse. Our implementation called MONI can rapidly find MEMs between reads and large sequence collections of highly repetitive sequences. Compared to other read aligners -- PuffAligner, Bowtie2, BWA-MEM, and CHIC -- MONI used 2--11 times less memory and was 2--32 times faster for index construction. Moreover, MONI was less than one thousandth the size of competing indexes for large collections of human chromosomes. Thus, MONI represents a major advance in our ability to perform MEM finding against very large collections of related references. Availability: MONI is publicly available at https://github.com/maxrossi91/moni.

TCAB1 is necessary for telomerase assembly

The ribonucleoprotein telomerase counteracts telomere shortening by adding repetitive sequences to the ends of human chromosomes. Telomerase is composed of the reverse transcriptase TERT, the telomerase RNA, and several auxiliary proteins that associate with the telomerase RNA, including TCAB1. TCAB1 is necessary for telomere maintenance in human cells and has been proposed to play a role in telomerase trafficking to Cajal bodies and telomeres, and in telomerase RNA folding. Here we show that, contrary to previous findings, TCAB1 is essential for telomerase assembly. We demonstrate that in the absence of TCAB1, the telomerase RNA is trapped in the nucleolus, a phase separated nuclear organelle, while TERT localizes to the nucleoplasm and is excluded from the nucleolus. Thus, nucleolar phase separation constitutes a barrier that counteracts telomerase assembly and TCAB1 is required to extract the telomerase RNA from the nucleolus, providing a molecular mechanism for the essential role of TCAB1 in telomerase function.