Dataset Information

Sparse Tensor Decomposition for Haplotype Assembly of Diploids and Polyploids.

ABSTRACT:

Background

Haplotype assembly is the task of reconstructing haplotypes of an individual from a mixture of sequenced chromosome fragments. Haplotype information enables studies of the effects of genetic variations on an organism's phenotype. Most of the mathematical formulations of haplotype assembly are known to be NP-hard and haplotype assembly becomes even more challenging as the sequencing technology advances and the length of the paired-end reads and inserts increases. Assembly of haplotypes polyploid organisms is considerably more difficult than in the case of diploids. Hence, scalable and accurate schemes with provable performance are desired for haplotype assembly of both diploid and polyploid organisms.

Results

We propose a framework that formulates haplotype assembly from sequencing data as a sparse tensor decomposition. We cast the problem as that of decomposing a tensor having special structural constraints and missing a large fraction of its entries into a product of two factors, U and [Formula: see text]; tensor [Formula: see text] reveals haplotype information while U is a sparse matrix encoding the origin of erroneous sequencing reads. An algorithm, AltHap, which reconstructs haplotypes of either diploid or polyploid organisms by iteratively solving this decomposition problem is proposed. The performance and convergence properties of AltHap are theoretically analyzed and, in doing so, guarantees on the achievable minimum error correction scores and correct phasing rate are established. The developed framework is applicable to diploid, biallelic and polyallelic polyploid species. The code for AltHap is freely available from https://github.com/realabolfazl/AltHap .

Conclusion

AltHap was tested in a number of different scenarios and was shown to compare favorably to state-of-the-art methods in applications to haplotype assembly of diploids, and significantly outperforms existing techniques when applied to haplotype assembly of polyploids.

SUBMITTER: Hashemi A

PROVIDER: S-EPMC5872563 | biostudies-literature | 2018 Mar

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Sparse Tensor Decomposition for Haplotype Assembly of Diploids and Polyploids.

Hashemi Abolfazl A Zhu Banghua B Vikalo Haris H

BMC genomics 20180321 Suppl 4

<h4>Background</h4>Haplotype assembly is the task of reconstructing haplotypes of an individual from a mixture of sequenced chromosome fragments. Haplotype information enables studies of the effects of genetic variations on an organism's phenotype. Most of the mathematical formulations of haplotype assembly are known to be NP-hard and haplotype assembly becomes even more challenging as the sequencing technology advances and the length of the paired-end reads and inserts increases. Assembly of ha ...[more]

PMID: 29589554

Similar Datasets

Project description:Functional magnetic resonance imaging (fMRI) is a powerful technique with the potential to estimate individual variations in behavioral and cognitive traits. Joint learning of multiple datasets can utilize their complementary information so as to improve learning performance, but it also gives rise to the challenge for data fusion to effectively integrate brain patterns elicited by multiple fMRI data. However, most of the current data fusion methods analyze each single dataset separately and further infer the relationship among them, which fail to utilize the multidimensional structure inherent across modalities and may ignore complex but important interactions. To address this issue, we propose a novel sparse tensor decomposition method to integrate multiple task-stimulus (paradigm) fMRI data. Seeing each paradigm fMRI as one modality, our proposed method considers the relationships across subjects and modalities simultaneously. In specific, a third-order tensor is first modeled by using the functional network connectivity (FNC) of subjects in multiple fMRI paradigms. A novel sparse tensor decomposition with the regularization terms is designed to factorize the tensor into a series of rank-one components, which can extract the shared components across modalities as the embedded features. The L2,1-norm regularizer (i.e., group sparsity) is enforced to select a few common features among multiple subjects. Validation of the proposed method is performed on realistic three paradigm fMRI datasets from the Philadelphia Neurodevelopmental Cohort (PNC) study, for the study of the relationship between the FNC and human cognitive abilities. Experimental results show our method outperforms several other competing methods in the prediction of individuals with different cognitive behaviors via the wide range achievement test (WRAT). Furthermore, our method discovers the FNC related to the cognitive behaviors, such as the connectivity associated with the default mode network (DMN) for three paradigms, and the connectivity between DMN and visual (VIS) domains within the emotion task.

Project description:Background and aimsApomixis or asexual seed reproduction is a key evolutionary mechanism in certain angiosperms providing them with reproductive assurance and isolation. Nevertheless, the frequency of apomixis is largely unknown, especially in groups with autonomous apomixis such as the diploid-polyploid genus Hieracium.MethodsUsing flow cytometric analyses, we determined the ploidy level and reproductive pathways (sexual vs. apomictic) for 7616 seeds originating from 946 plants belonging to >50 taxa sampled at 130 sites across Europe.Key resultsDiploid seeds produced by diploids were formed exclusively by the sexual pathway after double fertilization of reduced embryo sacs. An absolute majority of tri- and tetraploid seeds (99.6 %) produced by tri- and tetrapolyploid taxa were formed by autonomous apomixis. Only 20 polyploid seeds (0.4 %) were formed sexually. These seeds, which originated on seven polyploid accessions of four taxa, were formed after fertilization of either unreduced embryo sacs through a so-called triploid bridge or reduced embryo sacs, and frequently resulted in progeny with an increased ploidy. In addition, the formation of seedlings with increased ploidy (4x and 6x) was found in two triploid plants. This is the first firm evidence on functional facultative apomixis in polyploid members of Hieracium sensu stricto (s.s.).ConclusionsThe mode of reproduction in Hieracium s.s. is tightly associated with ploidy. While diploids produce seeds exclusively sexually, polyploids produce seeds by obligate or almost obligate apomixis. Strict apomixis can increase the reproductive assurance and this in turn can increase the colonization ability of apomicts. Nevertheless, our data clearly show that certain polyploid plants are still able to reproduce sexually and contribute to the formation of new cytotypes and genotypes. The finding of functional facultative apomicts is essential for future studies focused on evolution, inheritance and ecological significance of apomixis in this genus.

Dataset Information

Sparse Tensor Decomposition for Haplotype Assembly of Diploids and Polyploids.

Background

Results

Conclusion

Publications

Sparse Tensor Decomposition for Haplotype Assembly of Diploids and Polyploids.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets