PRIMUS: improving pedigree reconstruction using mitochondrial and Y haplotypes.
Ontology highlight
ABSTRACT:
SUBMITTER: Staples J
PROVIDER: S-EPMC5963362 | biostudies-literature | 2016 Feb
REPOSITORIES: biostudies-literature
Publications
PRIMUS: improving pedigree reconstruction using mitochondrial and Y haplotypes.
Staples Jeffrey J Ekunwe Lynette L Lange Ethan E Wilson James G JG Nickerson Deborah A DA Below Jennifer E JE
Bioinformatics (Oxford, England) 20151029 4
<h4>Unlabelled</h4>PRIMUS is a pedigree reconstruction algorithm that uses estimates of genome-wide identity by descent to reconstruct pedigrees consistent with observed genetic data. However, when genetic data for individuals within a pedigree are missing, often multiple pedigrees can be reconstructed that fit the data. We report a major expansion of PRIMUS that uses mitochondrial (mtDNA) and non-recombining Y chromosome (NRY) haplotypes to eliminate many pedigree structures that are inconsiste ...[more]
Similar Datasets
Project description:Understanding and correctly utilizing relatedness among samples is essential for genetic analysis; however, managing sample records and pedigrees can often be error prone and incomplete. Data sets ascertained by random sampling often harbor cryptic relatedness that can be leveraged in genetic analyses for maximizing power. We have developed a method that uses genome-wide estimates of pairwise identity by descent to identify families and quickly reconstruct and score all possible pedigrees that fit the genetic data by using up to third-degree relatives, and we have included it in the software package PRIMUS (Pedigree Reconstruction and Identification of the Maximally Unrelated Set). Here, we validate its performance on simulated, clinical, and HapMap pedigrees. Among these samples, we demonstrate that PRIMUS can verify reported pedigree structures and identify cryptic relationships. Finally, we show that PRIMUS reconstructed pedigrees, all of which were previously unknown, for 203 families from a cohort collected in Starr County, TX (1,890 samples).
| S-EPMC4225580 | biostudies-literature
Project description:In humans, mitochondria contain their own DNA (mtDNA) that is inherited exclusively from the mother. The mitochondrial genome encodes 13 polypeptides that are components of oxidative phosphorylation to produce energy. Any disruption in these genes might interfere with energy production and thus contribute to metabolic derangement. Mitochondria also regulate several important cellular activities including cell death and calcium homeostasis. Aided by sharply declining costs of high-density genotyping, hundreds of mitochondrial variants will soon be available in several cohorts with pedigree structures. Association testing of mitochondrial variants with disease traits using pedigree data raises unique challenges because of the difficulty in separating the effects of nuclear and mitochondrial genomes, which display different modes of inheritance. Failing to correctly account for these effects might decrease power or inflate type I error in association tests. In this report, we sought to identify the best strategy for association testing of mitochondrial variants when genotype and phenotype data are available in pedigrees. We proposed several strategies to account for polygenic effects of the nuclear and mitochondrial genomes and we performed extensive simulation studies to evaluate type I error and power of these strategies. In addition, we proposed two permutation tests to obtain empirical P values for these strategies. Furthermore, we applied two of the analytical strategies to association analysis of 196 mitochondrial variants with blood pressure and fasting blood glucose in the pedigree rich, Framingham Heart Study. Finally, we discussed strategies for study design, genotyping, and data cleaning in association testing of mtDNA in pedigrees.
| S-EPMC4171957 | biostudies-literature
Project description:The aim of this study was to obtain unbiased estimates of the diversity parameters, the population history, and the degree of admixture in Cika cattle which represents the local admixed breeds at risk of extinction undergoing challenging conservation programs. Genetic analyses were performed on the genome-wide Single Nucleotide Polymorphism (SNP) Illumina Bovine SNP50 array data of 76 Cika animals and 531 animals from 14 reference populations. To obtain unbiased estimates we used short haplotypes spanning four markers instead of single SNPs to avoid an ascertainment bias of the BovineSNP50 array. Genome-wide haplotypes combined with partial pedigree and type trait classification show the potential to improve identification of purebred animals with a low degree of admixture. Phylogenetic analyses demonstrated unique genetic identity of Cika animals. Genetic distance matrix presented by rooted Neighbour-Net suggested long and broad phylogenetic connection between Cika and Pinzgauer. Unsupervised clustering performed by the admixture analysis and two-dimensional presentation of the genetic distances between individuals also suggest Cika is a distinct breed despite being similar in appearance to Pinzgauer. Animals identified as the most purebred could be used as a nucleus for a recovery of the native genetic background in the current admixed population. The results show that local well-adapted strains, which have never been intensively managed and differentiated into specific breeds, exhibit large haplotype diversity. They suggest a conservation and recovery approach that does not rely exclusively on the search for the original native genetic background but rather on the identification and removal of common introgressed haplotypes would be more powerful. Successful implementation of such an approach should be based on combining phenotype, pedigree, and genome-wide haplotype data of the breed of interest and a spectrum of reference breeds which potentially have had direct or indirect historical contribution to the genetic makeup of the breed of interest.
| S-EPMC4414476 | biostudies-literature
Project description:Typical bacterial strain differentiation methods are often challenged by high genetic similarity between strains. To address this problem, we introduce a novel in silico peptide fingerprinting method based on conventional wet-lab protocols that enables the identification of potential strain-specific peptides. These can be further investigated using in vitro approaches, laying a foundation for the development of biomarker detection and application-specific methods. This novel method aims at reducing large amounts of comparative peptide data to binary matrices while maintaining a high phylogenetic resolution. The underlying case study concerns the Bacillus cereus group, namely the differentiation of Bacillus thuringiensis, Bacillus anthracis and Bacillus cereus strains. Results show that trees based on cytoplasmic and extracellular peptidomes are only marginally in conflict with those based on whole proteomes, as inferred by the established Genome-BLAST Distance Phylogeny (GBDP) method. Hence, these results indicate that the two approaches can most likely be used complementarily even in other organismal groups. The obtained results confirm previous reports about the misclassification of many strains within the B. cereus group. Moreover, our method was able to separate the B. anthracis strains with high resolution, similarly to the GBDP results as benchmarked via Bayesian inference and both Maximum Likelihood and Maximum Parsimony. In addition to the presented phylogenomic applications, whole-peptide fingerprinting might also become a valuable complementary technique to digital DNA-DNA hybridization, notably for bacterial classification at the species and subspecies level in the future.
| S-EPMC5198984 | biostudies-literature
Project description:Current protocol of Amide Proton Transfer-weighted (APTw) imaging commonly starts with the acquisition of high-resolution T2-weighted (T2w) images followed by APTw imaging at particular geometry and locations (i.e. slice) determined by the acquired T2w images. Although many advanced MRI reconstruction methods have been proposed to accelerate MRI, existing methods for APTw MRI lacks the capability of taking advantage of structural information in the acquired T2w images for reconstruction. In this paper, we present a novel APTw image reconstruction framework that can accelerate APTw imaging by reconstructing APTw images directly from highly undersampled k-space data and corresponding T2w image at the same location. The proposed framework starts with a novel sparse representation-based slice matching algorithm that aims to find the matched T2w slice given only the undersampled APTw image. A Recurrent Feature Sharing Reconstruction network (RFS-Rec) is designed to utilize intermediate features extracted from the matched T2w image by a Convolutional Recurrent Neural Network (CRNN), so that the missing structural information can be incorporated into the undersampled APT raw image thus effectively improving the image quality of the reconstructed APTw image. We evaluate the proposed method on two real datasets consisting of brain data from rats and humans. Extensive experiments demonstrate that the proposed RFS-Rec approach can outperform the state-of-the-art methods.
| S-EPMC7577831 | biostudies-literature
Project description:BackgroundHaplotypes extracted from human DNA can be used for gene mapping and other analysis of genetic patterns within and across populations. A fundamental problem is, however, that current practical laboratory methods do not give haplotype information. Estimation of phased haplotypes of unrelated individuals given their unphased genotypes is known as the haplotype reconstruction or phasing problem.ResultsWe define three novel statistical models and give an efficient algorithm for haplotype reconstruction, jointly called HaploRec. HaploRec is based on exploiting local regularities conserved in haplotypes: it reconstructs haplotypes so that they have maximal local coherence. This approach--not assuming statistical dependence for remotely located markers--has two useful properties: it is well-suited for sparse marker maps, such as those used in gene mapping, and it can actually take advantage of long maps.ConclusionOur experimental results with simulated and real data show that HaploRec is a powerful method for the large scale haplotyping needed in association studies. With sample sizes large enough for gene mapping it appeared to be the best compared to all other tested methods (Phase, fastPhase, PL-EM, Snphap, Gerbil; simulated data), with small samples it was competitive with the best available methods (real data). HaploRec is several orders of magnitude faster than Phase and comparable to the other methods; the running times are roughly linear in the number of subjects and the number of markers. HaploRec is publicly available at http://www.cs.helsinki.fi/group/genetics/haplotyping.html.
| S-EPMC1766938 | biostudies-literature
Project description:Data on hundreds or thousands of single nucleotide polymorphisms (SNPs) provide detailed information about the relationships between individuals, but currently few tools can turn this information into a multigenerational pedigree. I present the r package sequoia, which assigns parents, clusters half-siblings sharing an unsampled parent and assigns grandparents to half-sibships. Assignments are made after consideration of the likelihoods of all possible first-, second- and third-degree relationships between the focal individuals, as well as the traditional alternative of being unrelated. This careful exploration of the local likelihood surface is implemented in a fast, heuristic hill-climbing algorithm. Distinction between the various categories of second-degree relatives is possible when likelihoods are calculated conditional on at least one parent of each focal individual. Performance was tested on simulated data sets with realistic genotyping error rate and missingness, based on three different large pedigrees (N = 1000-2000). This included a complex pedigree with overlapping generations, occasional close inbreeding and some unknown birth years. Parentage assignment was highly accurate down to about 100 independent SNPs (error rate <0.1%) and fast (<1 min) as most pairs can be excluded from being parent-offspring based on opposite homozygosity. For full pedigree reconstruction, 40% of parents were assumed nongenotyped. Reconstruction resulted in low error rates (<0.3%), high assignment rates (>99%) in limited computation time (typically <1 h) when at least 200 independent SNPs were used. In three empirical data sets, relatedness estimated from the inferred pedigree was strongly correlated to genomic relatedness.
| S-EPMC6849609 | biostudies-literature
Project description:Accurate and automated reconstruction of the in vivo human cerebral cortical surface from anatomical magnetic resonance (MR) images facilitates the quantitative analysis of cortical structure. Anatomical MR images with sub-millimeter isotropic spatial resolution improve the accuracy of cortical surface and thickness estimation compared to the standard 1-millimeter isotropic resolution. Nonetheless, sub-millimeter resolution acquisitions require averaging multiple repetitions to achieve sufficient signal-to-noise ratio and are therefore long and potentially vulnerable to subject motion. We address this challenge by synthesizing sub-millimeter resolution images from standard 1-millimeter isotropic resolution images using a data-driven supervised machine learning-based super-resolution approach achieved via a deep convolutional neural network. We systematically characterize our approach using a large-scale simulated dataset and demonstrate its efficacy in empirical data. The super-resolution data provide improved cortical surfaces similar to those obtained from native sub-millimeter resolution data. The whole-brain mean absolute discrepancy in cortical surface positioning and thickness estimation is below 100 μm at the single-subject level and below 50 μm at the group level for the simulated data, and below 200 μm at the single-subject level and below 100 μm at the group level for the empirical data, making the accuracy of cortical surfaces derived from super-resolution sufficient for most applications.
| S-EPMC7727379 | biostudies-literature
Project description:NMR spectroscopy is central to atomic resolution studies in biology and chemistry. Key to this approach are multidimensional experiments. Obtaining such experiments with sufficient resolution, however, is a slow process, in part since each time increment in every indirect dimension needs to be recorded twice, in quadrature. We introduce a modified compressed sensing (CS) algorithm enabling reconstruction of data acquired with random acquisition of quadrature components in gradient-selection NMR. We name this approach random quadrature detection (RQD). Gradient-selection experiments are essential to the success of modern NMR and with RQD, a 50 % reduction in the number of data points per indirect dimension is possible, by only acquiring one quadrature component per time point. Using our algorithm (CSRQD), high quality reconstructions are achieved. RQD is modular and combined with non-uniform sampling we show that this provides increased flexibility in designing sampling schedules leading to improved resolution with increasing benefits as dimensionality of experiments increases, with particular advantages for 4- and higher dimensional experiments.
| S-EPMC5504140 | biostudies-other
Project description:The ornamental crop species Hydrangea macrophylla exhibits diploid and triploid levels of ploidy and develops lacecap (wild type) or mophead inflorescences. In order to characterize a H. macrophylla germplasm collection, we determined the inflorescence type and the 2C DNA content of 120 plants representing 43 cultivars. We identified 78 putative diploid and 39 putative triploid plants by flow cytometry. In our collection 69 out of 98 flowering plants produced lacecap inflorescences, whereas 29 plants developed mophead inflorescences. Surprisingly, 12 cultivars included diploid as well as triploid plants, while 5 cultivars contained plants with different inflorescence types. We genotyped this germplasm collection using 12 SSR markers that detected 2-7 alleles per marker, and identified 51 different alleles in this collection. We detected 62 distinct fingerprints, revealing a higher genetic variation than the number of cultivars suggested. Only one genotype per cultivar is expected due to the vegetative propagation of Hydrangea cultivars; however we identified 25 cultivars containing 2-4 different genotypes. These different genotypes explained the variation in DNA content and inflorescence type. Diploid and triploid plants with the same cultivar name were exclusively mix-ups. We therefor assume, that 36% of the tested plants were mislabeled. Based on the "Wädenswil" pedigree, which includes 31 of the tested cultivars, we predicted cultivar-specific fingerprints and identified at least 21 out of 31 cultivars by SSR marker-based reconstruction of the "Wädenswil" pedigree. Furthermore, we detected 4 putative interploid crosses between diploid and triploid plants in this pedigree. These interploid crosses resulted in diploid or/and triploid offspring, suggesting that crosses with triploids were successfully applied in breeding of H. macrophylla.
| S-EPMC5915539 | biostudies-other