Source author record

Nikita Alexeev

Nikita Alexeev appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

8works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

8 published item(s)

preprint2019arXiv

Orienting Ordered Scaffolds: Complexity and Algorithms

Despite the recent progress in genome sequencing and assembly, many of the currently available assembled genomes come in a draft form. Such draft genomes consist of a large number of genomic fragments (scaffolds), whose order and/or orientation (i.e., strand) in the genome are unknown. There exist various scaffold assembly methods, which attempt to determine the order and orientation of scaffolds along the genome chromosomes. Some of these methods (e.g., based on FISH physical mapping, chromatin conformation capture, etc.) can infer the order of scaffolds, but not necessarily their orientation. This leads to a special case of the scaffold orientation problem (i.e., deducing the orientation of each scaffold) with a known order of the scaffolds. We address the problem of orientating ordered scaffolds as an optimization problem based on given weighted orientations of scaffolds and their pairs (e.g., coming from pair-end sequencing reads, long reads, or homologous relations). We formalize this problem using notion of a scaffold graph (i.e., a graph, where vertices correspond to the assembled contigs or scaffolds and edges represent connections between them). We prove that this problem is NP-hard, and present a polynomial-time algorithm for solving its special case, where orientation of each scaffold is imposed relatively to at most two other scaffolds. We further develop an FPT algorithm for the general case of the OOS problem.

preprint2016arXiv

Combinatorial Scoring of Phylogenetic Networks

Construction of phylogenetic trees and networks for extant species from their characters represents one of the key problems in phylogenomics. While solution to this problem is not always uniquely defined and there exist multiple methods for tree/network construction, it becomes important to measure how well the constructed networks capture the given character relationship across the species. In the current study, we propose a novel method for measuring the specificity of a given phylogenetic network in terms of the total number of distributions of character states at the leaves that the network may impose. While for binary phylogenetic trees, this number has an exact formula and depends only on the number of leaves and character states but not on the tree topology, the situation is much more complicated for non-binary trees or networks. Nevertheless, we develop an algorithm for combinatorial enumeration of such distributions, which is applicable for arbitrary trees and networks under some reasonable assumptions.

preprint2016arXiv

Singular Values Distribution of Squares of Elliptic Random Matrices and Type B Narayana Polynomials

We consider Gaussian elliptic random matrices $X$ of a size $N \times N$ with parameter $ρ$, i.e., matrices whose pairs of entries $(X_{ij}, X_{ji})$ are mutually independent Gaussian vectors, $E X_{ij} = 0$, $E X^2_{ij} = 1$ and $E X_{ij} X_{ji} = ρ$. We are interested in the asymptotic distribution of eigenvalues of the matrix $W =\frac{1}{N^2} X^2 X^{*2}$. We have shown that this distribution is defined by its moments and we provide a recurrent relation for these moments. We have proven that the (symmetrized) asymptotic distribution is determined by its free cumulants, which are Narayana polynomials of type B: $$c_{2n} = \sum_{k=0}^n \binom{n}{k}^2 ρ^{2k}.$$

preprint2015arXiv

A Computational Method for the Rate Estimation of Evolutionary Transpositions

Genome rearrangements are evolutionary events that shuffle genomic architectures. Most frequent genome rearrangements are reversals, translocations, fusions, and fissions. While there are some more complex genome rearrangements such as transpositions, they are rarely observed and believed to constitute only a small fraction of genome rearrangements happening in the course of evolution. The analysis of transpositions is further obfuscated by intractability of the underlying computational problems. We propose a computational method for estimating the rate of transpositions in evolutionary scenarios between genomes. We applied our method to a set of mammalian genomes and estimated the transpositions rate in mammalian evolution to be around 0.26.

preprint2011arXiv

Hultman numbers, polygon gluings and matrix integrals

The Hultman numbers enumerate permutations whose cycle graph has a given number of alternating cycles (they are relevant to the Bafna-Pevzner approach to genome comparison and genome rearrangements). We give two new interpretations of the Hultman numbers: in terms of polygon gluings and as integrals over the space of complex matrices, and derive some properties of their generating functions.

preprint2011arXiv

On the asymptotic distribution of singular values of products of large rectangular random matrices

We consider products of independent large random rectangular matrices with independent entries. The limit distribution of the expected empirical distribution of singular values of such products is computed. The distribution function is described by its Stieltjes transform, which satisfies some algebraic equation. In the particular case of square matrices we get a well-known distribution which moments are Fuss-Catalan numbers.

preprint2010arXiv

Asymptotic distribution of singular values of powers of random matrices

Let $x$ be a complex random variable such that ${\E {x}=0}$, ${\E |x|^2=1}$, ${\E |x|^{4} < \infty}$. Let $x_{ij}$, $i,j \in \{1,2,...\}$ be independet copies of $x$. Let ${\Xb=(N^{-1/2}x_{ij})}$, $1\leq i,j \leq N$ be a random matrix. Writing $\Xb^*$ for the adjoint matrix of $\Xb$, consider the product $\Xb^m{\Xb^*}^m$ with some $m \in \{1,2,...\}$. The matrix $\Xb^m{\Xb^*}^m$ is Hermitian positive semi-definite. Let $λ_1,λ_2,...,λ_N$ be eigenvalues of $\Xb^m{\Xb^*}^m$ (or squared singular values of the matrix $\Xb^m$). In this paper we find the asymptotic distribution function \[ G^{(m)}(x)=\lim_{N\to\infty}\E{F_N^{(m)}(x)} \] of the empirical distribution function \[ {F_N^{(m)}(x)} = N^{-1} \sum_{k=1}^N {\mathbb{I}{\{λ_k \leq x\}}}, \] where $\mathbb{I} \{A\}$ stands for the indicator function of event $A$. The moments of $G^{(m)}$ satisfy \[ M^{(m)}_p=\int_{\mathbb{R}}{x^p dG^{(m)}(x)}=\frac{1}{mp+1}\binom{mp+p}{p}. \] In Free Probability Theory $M^{(m)}_p$ are known as Fuss--Catalan numbers. With $m=1$ our result turns to a well known result of Marchenko--Pastur 1967.

preprint2010arXiv

On the asymptotic distribution of the singular values of powers of random matrices

We consider powers of random matrices with independent entries. Let $X_{ij}, i,j\ge 1$, be independent complex random variables with $\E X_{ij}=0$ and $\E |X_{ij}|^2=1$ and let $\mathbf X$ denote an $n\times n$ matrix with $[\mathbf X]_{ij}=X_{ij}$, for $1\le i, j\le n$. Denote by $s_1^{(m)}\ge...\ge s_n^{(m)}$ the singular values of the random matrix $\mathbf W:={n^{-\frac m2}} \mathbf X^m$ and define the empirical distribution of the squared singular values by $$ \mathcal F_n^{(m)}(x)=\frac1n\sum_{k=1}^nI_{\{s_k^{(m)}^2\le x\}}, $$ where $I_{\{B\}}$ denotes the indicator of an event $B$. We prove that under a Lindeberg condition for the fourth moment that the expected spectral distribution $F_n^{(m)}(x)=\E \mathcal F_n^{(m)}(x)$ converges to the distribution function $G^{(m)}(x)$ defined by its moments $$ α_k(m):=\int_{\mathbb R}x^k\,d\,G(x)=\frac {1}{mk+1}\binom{km+k}{k}. $$