PediHaplotyper: software for consistent assignment of marker haplotypes in pedigrees.
Ontology highlight
ABSTRACT:
SUBMITTER: Voorrips RE
PROVIDER: S-EPMC4977329 | biostudies-literature | 2016
REPOSITORIES: biostudies-literature
Publications
PediHaplotyper: software for consistent assignment of marker haplotypes in pedigrees.
Voorrips Roeland E RE Bink Marco C A M MC Kruisselbrink Johannes W JW Koehorst-van Putten Herma J J HJ van de Weg W Eric WE
Molecular breeding : new strategies in plant improvement 20160808
In the study of large outbred pedigrees with many founders, individual bi-allelic markers, such as SNP markers, carry little information. After phasing the marker genotypes, multi-allelic loci consisting of groups of closely linked markers can be identified, which are called "haploblocks". Here, we describe PediHaplotyper, an R package capable of assigning consistent alleles to such haploblocks, allowing for missing and incorrect SNP data. These haploblock genotypes are much easier to interpret ...[more]
Similar Datasets
Project description:In ecology and forensics, some population assignment techniques use molecular markers to assign individuals to known groups. However, assigning individuals to known populations can be difficult if the level of genetic differentiation among populations is small. Most assignment studies handle independent markers, often by pruning markers in Linkage Disequilibrium (LD), ignoring the information contained in the correlation among markers due to LD.To improve the accuracy of population assignment, we present an algorithm, implemented in the HaploPOP software, that combines markers into haplotypes, without requiring independence. The algorithm is based on the Gain of Informativeness for Assignment that provides a measure to decide if a pair of markers should be combined into haplotypes, or not, in order to improve assignment. Because complete exploration of all possible solutions for constructing haplotypes is computationally prohibitive, our approach uses a greedy algorithm based on windows of fixed sizes. We evaluate the performance of HaploPOP to assign individuals to populations using a split-validation approach. We investigate both simulated SNPs data and dense genotype data from individuals from Spain and Portugal.Our results show that constructing haplotypes with HaploPOP can substantially reduce assignment error. The HaploPOP software is freely available as a command-line software at www.ieg.uu.se/Jakobsson/software/HaploPOP/.
| S-EPMC4521458 | biostudies-literature
Project description:Detection of genotyping errors is a necessary step to minimize false results in genetic analysis. This is especially important when the rate of genotyping errors is high, as has been reported for high-throughput sequence data. To detect genotyping errors in pedigrees, Mendelian inconsistent (MI) error checks exist, as do multi-point methods that flag Mendelian consistent (MC) errors for sparse multi-allelic markers. However, few methods exist for detecting MC genotyping errors, particularly for dense variants on large pedigrees. Here, we introduce an efficient method to detect MC errors even for very dense variants (e.g., SNPs and sequencing data) on pedigrees that may be large. Our method first samples inheritance vectors (IVs) using a moderately sparse but informative set of markers using a Markov chain Monte Carlo-based sampler. Using sampled IVs, we considered two test statistics to detect MC genotyping errors: the percentage of IVs inconsistent with observed genotypes (A1) or the posterior probability of error configurations (A2). Using simulations, we show that this method, even with the simpler A1 statistic, is effective for detecting MC genotyping errors in dense variants, with sensitivity almost as high as the theoretical best sensitivity possible. We also evaluate the effectiveness of this method as a function of parameters, when including the observed pattern for genotype, density of framework markers, error rate, allele frequencies, and number of sampled inheritance vectors. Our approach provides a line of defense against false findings based on the use of dense variants in pedigrees.
| S-EPMC4081466 | biostudies-literature
Project description:Mass spectrometry combined with enrichment strategies for phosphorylated peptides has been successfully employed for two decades to identify sites of phosphorylation. However, unambiguous phosphosite assignment is considered challenging. Given that site-specific phosphorylation events function as different molecular switches, validation of phosphorylation sites is of utmost importance. In our earlier study we developed a method based on simulated phosphopeptide spectral libraries, which enables highly sensitive and accurate phosphosite assignments. To promote more widespread use of this method, we here introduce a software implementation with improved usability and performance.We present SimPhospho, a fast and user-friendly tool for accurate simulation of phosphopeptide tandem mass spectra. Simulated phosphopeptide spectral libraries are used to validate and supplement database search results, with a goal to improve reliable phosphoproteome identification and reporting. The presented program can be easily used together with the Trans-Proteomic Pipeline and integrated in a phosphoproteomics data analysis workflow.SimPhospho is open source and it is available for Windows, Linux and Mac operating systems. The software and its user's manual with detailed description of data analysis as well as test data can be found at https://sourceforge.net/projects/simphospho/.Supplementary data are available at Bioinformatics online.
| S-EPMC6061695 | biostudies-other
Project description:Familial clustering without any prerequisite knowledge becomes often necessary in Behavioral Science, and forensic studies in case of great disasters like Tsunami and earthquake requiring body-identification without any usable information. However, there has been no well-established method for this purpose although conventional ones such as short tandem repeats (STR) and single nucleotide polymorphism (SNP), which might be applied with toil and moil to some extent. In this situation, we could find that the universal genome distance-measuring method genome profiling (GP), which is made up of three elemental techniques; random PCR, micro-temperature gradient gel electrophoresis (?TGGE), and computer processing for normalization, can do this purpose with ease when applied to mouse families. We also confirmed that the sequencing approach based on the ccgf (commonly conserved genetic fragment appearing in the genome profile) was not completely discriminative in this case. This is the first demonstration that the familial clustering can be attained without a priori sequence information to the level of discriminating strains and sibling relationships. This method can complement the conventional approaches in preliminary familial clustering.
| S-EPMC4629661 | biostudies-literature
Project description:The major histocompatibility complex (MHC) contains the most polymorphic genetic system in humans, the human leukocyte antigen (HLA) genes of the adaptive immune system. High allelic diversity in HLA is argued to be maintained by balancing selection, such as negative frequency-dependent selection or heterozygote advantage. Selective pressure against immune escape by pathogens can maintain appreciable frequencies of many different HLA alleles. The selection pressures operating on combinations of HLA alleles across loci, or haplotypes, have not been extensively evaluated since the high HLA polymorphism necessitates very large sample sizes, which have not been available until recently. We aimed to evaluate the effect of selection operating at the HLA haplotype level by analyzing HLA A~C~B~DRB1~DQB1 haplotype frequencies derived from over six million individuals genotyped by the National Marrow Donor Program registry. In contrast with alleles, HLA haplotype diversity patterns suggest purifying selection, as certain HLA allele combinations co-occur in high linkage disequilibrium. Linkage disequilibrium is positive (Dij'>0) among frequent haplotypes and negative (Dij'<0) among rare haplotypes. Fitting the haplotype frequency distribution to several population dynamics models, we found that the best fit was obtained when significant positive frequency-dependent selection (FDS) was incorporated. Finally, the Ewens-Watterson test of homozygosity showed excess homozygosity for 5-locus haplotypes within 23 US populations studied, with an average Fnd of 28.43. Haplotype diversity is most consistent with purifying selection for HLA Class I haplotypes (HLA-A, -B, -C), and was not inferred for HLA Class II haplotypes (-DRB1 and-DQB1). We discuss our empirical results in the context of evolutionary theory, exploring potential mechanisms of selection that maintain high linkage disequilibrium in MHC haplotype blocks.
| S-EPMC5590998 | biostudies-literature
Project description:Glycosaminoglycans (GAGs) covalently linked to proteoglycans (PGs) are characterized by repeating disaccharide units and variable sulfation patterns along the chain. GAG length and sulfation patterns impact disease etiology, cellular signaling, and structural support for cells. We and others have demonstrated the usefulness of tandem mass spectrometry (MS2) for assigning the structures of GAG saccharides; however, manual interpretation of tandem mass spectra is time-consuming, so computational methods must be employed. In the proteomics domain, the identification of monoisotopic peaks and charge states relies on algorithms that use averagine, or the average building block of the compound class being analyzed. Although these methods perform well for protein and peptide spectra, they perform poorly on GAG tandem mass spectra, because a single average building block does not characterize the variable sulfation of GAG disaccharide units. In addition, it is necessary to assign product ion isotope patterns to interpret the tandem mass spectra of GAG saccharides. To address these problems, we developed GAGfinder, the first tandem mass spectrum peak finding algorithm developed specifically for GAGs. We define peak finding as assigning experimental isotopic peaks directly to a given product ion composition, as opposed to deconvolution or peak picking, which are terms more accurately describing the existing methods previously mentioned. GAGfinder is a targeted, brute force approach to spectrum analysis that uses precursor composition information to generate all theoretical fragments. GAGfinder also performs peak isotope composition annotation, which is typically a subsequent step for averagine-based methods. Data are available via ProteomeXchange with identifier PXD009101.
| S-EPMC6030731 | biostudies-literature
Project description:We have developed a software analysis package, HapScope, which includes a comprehensive analysis pipeline and a sophisticated visualization tool for analyzing functionally annotated haplotypes. The HapScope analysis pipeline supports: (i) computational haplotype construction with an expectation-maximization or Bayesian statistical algorithm; (ii) SNP classification by protein coding change, homology to model organisms or putative regulatory regions; and (iii) minimum SNP subset selection by either a Brute Force Algorithm or a Greedy Partition Algorithm. The HapScope viewer displays genomic structure with haplotype information in an integrated environment, providing eight alternative views for assessing genetic and functional correlation. It has a user-friendly interface for: (i) haplotype block visualization; (ii) SNP subset selection; (iii) haplotype consolidation with subset SNP markers; (iv) incorporation of both experimentally determined haplotypes and computational results; and (v) data export for additional analysis. Comparison of haplotypes constructed by the statistical algorithms with those determined experimentally shows variation in haplotype prediction accuracies in genomic regions with different levels of nucleotide diversity. We have applied HapScope in analyzing haplotypes for candidate genes and genomic regions with extensive SNP and genotype data. We envision that the systematic approach of integrating functional genomic analysis with population haplotypes, supported by HapScope, will greatly facilitate current genetic disease research.
| S-EPMC137968 | biostudies-literature
Project description:Not all proteins are amenable to uniform isotopic labeling with 13C and 15N, something needed for the widely used, and largely deductive, triple resonance assignment process. Among them are proteins expressed in mammalian cell culture where native glycosylation can be maintained, and proper formation of disulfide bonds facilitated. Uniform labeling in mammalian cells is prohibitively expensive, but sparse labeling with one or a few isotopically enriched amino acid types is an option for these proteins. However, assignment then relies on accessing the best match between a variety of measured NMR parameters and predictions based on 3D structure, often from X-ray crystallography. Finding this match is a challenging process that has benefitted from many computational tools, including trained neural nets for chemical shift prediction, genetic algorithms for searches through a myriad of assignment possibilities, and now AI-based prediction of high-quality structures for protein targets. AssignSLP_GUI, a new version of a software package for assignment of resonances from sparsely-labeled proteins, uses many of these tools. These tools and new additions to the package are highlighted in an application to a sparsely-labeled domain from a glycoprotein, CEACAM1.
| S-EPMC9742323 | biostudies-literature
Project description:BACKGROUND:Identification of biomarkers associated with the prognosis of different cancer subtypes is critical to achieve better therapeutic assistance. In colorectal cancer (CRC) the discovery of stable and consistent survival markers remains a challenge due to the high heterogeneity of this class of tumors. In this work, we identified a new set of gene markers for CRC associated to prognosis and risk using a large unified cohort of patients with transcriptomic profiles and survival information. RESULTS:We built an integrated dataset with 1273 human colorectal samples, which provides a homogeneous robust framework to analyse genome-wide expression and survival data. Using this dataset we identified two sets of genes that are candidate prognostic markers for CRC in stages III and IV, showing either up-regulation correlated with poor prognosis or up-regulation correlated with good prognosis. The top 10 up-regulated genes found as survival markers of poor prognosis (i.e. low survival) were: DCBLD2, PTPN14, LAMP5, TM4SF1, NPR3, LEMD1, LCA5, CSGALNACT2, SLC2A3 and GADD45B. The stability and robustness of the gene survival markers was assessed by cross-validation, and the best-ranked genes were also validated with two external independent cohorts: one of microarrays with 482 samples; another of RNA-seq with 269 samples. Up-regulation of the top genes was also proved in a comparison with normal colorectal tissue samples. Finally, the set of top 100 genes that showed overexpression correlated with low survival was used to build a CRC risk predictor applying a multivariate Cox proportional hazards regression analysis. This risk predictor yielded an optimal separation of the individual patients of the cohort according to their survival, with a p-value of 8.25e-14 and Hazard Ratio 2.14 (95% CI: 1.75-2.61). CONCLUSIONS:The results presented in this work provide a solid rationale for the prognostic utility of a new set of genes in CRC, demonstrating their potential to predict colorectal tumor progression and evolution towards poor survival stages. Our study does not provide a fixed gene signature for prognosis and risk prediction, but instead proposes a robust set of genes ranked according to their predictive power that can be selected for additional tests with other CRC clinical cohorts.
| S-EPMC6288855 | biostudies-literature
Project description:Inbreeding (mating between relatives) can dramatically reduce the fitness of offspring by causing parts of the genome to be identical by descent. Thus, measuring individual inbreeding is crucial for ecology, evolution and conservation biology. We used computer simulations to test whether the realized proportion of the genome that is identical by descent (IBDG) is predicted better by the pedigree inbreeding coefficient (FP) or by genomic (marker-based) measures of inbreeding. Genomic estimators of IBDG included the increase in individual homozygosity relative to mean Hardy-Weinberg expected homozygosity (FH), and two measures (FROH and FE) that use mapped genetic markers to estimate IBDG. IBDG was more strongly correlated with FH, FE and FROH than with FP across a broad range of simulated scenarios when thousands of SNPs were used. For example, IBDG was more strongly correlated with FROH, FH and FE (estimated with ?10?000 SNPs) than with FP (estimated with 20 generations of complete pedigree) in populations with a recent reduction in the effective populations size (from Ne=500 to Ne=75). FROH, FH and FE generally explained >90% of the variance in IBDG (among individuals) when 35?K or more SNPs were used. FP explained <80% of the variation in IBDG on average in all simulated scenarios, even when pedigrees included 20 generations. Our results demonstrate that IBDG can be more precisely estimated with large numbers of genetic markers than with pedigrees. We encourage researchers to adopt genomic marker-based measures of IBDG as thousands of loci can now be genotyped in any species.
| S-EPMC4815495 | biostudies-literature