Dataset Information

How many marker loci are necessary? Analysis of dominant marker data sets using two popular population genetic algorithms.

ABSTRACT: The number of marker loci required to answer a given research question satisfactorily is especially important for dominant markers since they have a lower information content than co-dominant marker systems. In this study, we used simulated dominant marker data sets to determine the number of dominant marker loci needed to obtain satisfactory results from two popular population genetic analyses: STRUCTURE and AMOVA (analysis of molecular variance). Factors such as migration, level of population differentiation, and unequal sampling were varied in the data sets to mirror a range of realistic research scenarios. AMOVA performed well under all scenarios with a modest quantity of markers while STRUCTURE required a greater number, especially when populations were closely related. The popular ?K method of determining the number of genetically distinct groups worked well when sampling was balanced, but underestimated the true number of groups with unbalanced sampling. These results provide a window through which to interpret previous work with dominant markers and we provide a protocol for determining the number of markers needed for future dominant marker studies.

SUBMITTER: Nelson MF

PROVIDER: S-EPMC3797491 | biostudies-literature | 2013 Sep

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

How many marker loci are necessary? Analysis of dominant marker data sets using two popular population genetic algorithms.

Nelson Michael F MF Anderson Neil O NO

Ecology and evolution 20130828 10

The number of marker loci required to answer a given research question satisfactorily is especially important for dominant markers since they have a lower information content than co-dominant marker systems. In this study, we used simulated dominant marker data sets to determine the number of dominant marker loci needed to obtain satisfactory results from two popular population genetic analyses: STRUCTURE and AMOVA (analysis of molecular variance). Factors such as migration, level of population ...[more]

PMID: 24223282

Similar Datasets

Project description:Background and purposeThe measurement of cortical shrinkage is a candidate marker of disease progression in Alzheimer's. This study evaluated the performance of two pipelines: Civet-CLASP (v1.1.9) and Freesurfer (v5.3.0).MethodsImages from 185 ADNI1 cases (69 elderly controls (CTR), 37 stable MCI (sMCI), 27 progressive MCI (pMCI), and 52 Alzheimer (AD) patients) scanned at baseline, month 12, and month 24 were processed using the two pipelines and two interconnected e-infrastructures: neuGRID (https://neugrid4you.eu) and VIP (http://vip.creatis.insa-lyon.fr). The vertex-by-vertex cross-algorithm comparison was made possible applying the 3D gradient vector flow (GVF) and closest point search (CPS) techniques.ResultsThe cortical thickness measured with Freesurfer was systematically lower by one third if compared to Civet's. Cross-sectionally, Freesurfer's effect size was significantly different in the posterior division of the temporal fusiform cortex. Both pipelines were weakly or mildly correlated with the Mini Mental State Examination score (MMSE) and the hippocampal volumetry. Civet differed significantly from Freesurfer in large frontal, parietal, temporal and occipital regions (p<0.05). In a discriminant analysis with cortical ROIs having effect size larger than 0.8, both pipelines gave no significant differences in area under the curve (AUC). Longitudinally, effect sizes were not significantly different in any of the 28 ROIs tested. Both pipelines weakly correlated with MMSE decay, showing no significant differences. Freesurfer mildly correlated with hippocampal thinning rate and differed in the supramarginal gyrus, temporal gyrus, and in the lateral occipital cortex compared to Civet (p<0.05). In a discriminant analysis with ROIs having effect size larger than 0.6, both pipelines yielded no significant differences in the AUC.ConclusionsCivet appears slightly more sensitive to the typical AD atrophic pattern at the MCI stage, but both pipelines can accurately characterize the topography of cortical thinning at the dementia stage.

Project description:BackgroundDue to the high cost and low reproducibility of many microarray experiments, it is not surprising to find a limited number of patient samples in each study, and very few common identified marker genes among different studies involving patients with the same disease. Therefore, it is of great interest and challenge to merge data sets from multiple studies to increase the sample size, which may in turn increase the power of statistical inferences. In this study, we combined two lung cancer studies using microarray GeneChip, employed two gene shaving methods and a two-step survival test to identify genes with expression patterns that can distinguish diseased from normal samples, and to indicate patient survival, respectively.ResultsIn addition to common data transformation and normalization procedures, we applied a distribution transformation method to integrate the two data sets. Gene shaving (GS) methods based on Random Forests (RF) and Fisher's Linear Discrimination (FLD) were then applied separately to the joint data set for cancer gene selection. The two methods discovered 13 and 10 marker genes (5 in common), respectively, with expression patterns differentiating diseased from normal samples. Among these marker genes, 8 and 7 were found to be cancer-related in other published reports. Furthermore, based on these marker genes, the classifiers we built from one data set predicted the other data set with more than 98% accuracy. Using the univariate Cox proportional hazard regression model, the expression patterns of 36 genes were found to be significantly correlated with patient survival (p < 0.05). Twenty-six of these 36 genes were reported as survival-related genes from the literature, including 7 known tumor-suppressor genes and 9 oncogenes. Additional principal component regression analysis further reduced the gene list from 36 to 16.ConclusionThis study provided a valuable method of integrating microarray data sets with different origins, and new methods of selecting a minimum number of marker genes to aid in cancer diagnosis. After careful data integration, the classification method developed from one data set can be applied to the other with high prediction accuracy.

Project description:MotivationAdvancements in high-throughput genomic sequencing are delivering genomic pathogen data at an unprecedented rate, positioning statistical phylogenetics as a critical tool to monitor infectious diseases globally. This rapid growth spurs the need for efficient inference techniques, such as Hamiltonian Monte Carlo (HMC) in a Bayesian framework, to estimate parameters of these phylogenetic models where the dimensions of the parameters increase with the number of sequences N. HMC requires repeated calculation of the gradient of the data log-likelihood with respect to (wrt) all branch-length-specific (BLS) parameters that traditionally takes O(N2) operations using the standard pruning algorithm. A recent study proposes an approach to calculate this gradient in O(N), enabling researchers to take advantage of gradient-based samplers such as HMC. The CPU implementation of this approach makes the calculation of the gradient computationally tractable for nucleotide-based models but falls short in performance for larger state-space size models, such as Markov-modulated and codon models. Here, we describe novel massively parallel algorithms to calculate the gradient of the log-likelihood wrt all BLS parameters that take advantage of graphics processing units (GPUs) and result in many fold higher speedups over previous CPU implementations.ResultsWe benchmark these GPU algorithms on three computing systems using three evolutionary inference examples exploring complete genomes from 997 dengue viruses, 62 carnivore mitochondria and 49 yeasts, and observe a >128-fold speedup over the CPU implementation for codon-based models and >8-fold speedup for nucleotide-based models. As a practical demonstration, we also estimate the timing of the first introduction of West Nile virus into the continental Unites States under a codon model with a relaxed molecular clock from 104 full viral genomes, an inference task previously intractable.Availability and implementationWe provide an implementation of our GPU algorithms in BEAGLE v4.0.0 (https://github.com/beagle-dev/beagle-lib), an open-source library for statistical phylogenetics that enables parallel calculations on multi-core CPUs and GPUs. We employ a BEAGLE-implementation using the Bayesian phylogenetics framework BEAST (https://github.com/beast-dev/beast-mcmc).

Project description:Rapid characterization of novel NB-LRR-associated resistance to Phomopsis cane spot on grapevine using high-throughput sampling and low-coverage sequencing for genotyping, locus mapping and transcriptome analysis provides insights into genetic resistance to a hemibiotrophic fungus. Phomopsis cane and leaf spot, caused by the hemibiotrophic fungus Diaporthe ampelina (syn = Phomopsis viticola), reduces the productivity in grapevines. Host resistance was studied on three F1 families derived from crosses involving resistant genotypes 'Horizon', Illinois 547-1, Vitis cinerea B9 and V. vinifera 'Chardonnay'. All families had progeny with extremely susceptible phenotypes, developing lesions on both dormant canes and maturing fruit clusters. Segregation of symptoms was observed under natural levels of inoculum in the field, while phenotypes on green shoots were confirmed under controlled inoculations in greenhouse. High-density genetic maps were used to localize novel qualitative resistance loci named Rda1 and Rda2 from V. cinerea B9 and 'Horizon', respectively. Co-linearity between reference genetic and physical maps allowed localization of Rda2 locus between 1.5 and 2.4 Mbp on chromosome 7, and Rda1 locus between 19.3 and 19.6 Mbp of chromosome 15, which spans a cluster of five NB-LRR genes. Further dissection of this locus was obtained by QTL mapping of gene expression values 14 h after inoculation across a subset of the 'Chardonnay' × V. cinerea B9 progeny. This provided evidence for the association between transcript levels of two of these NB-LRR genes with Rda1, with increased NB-LRR expression among susceptible progeny. In resistant parent V. cinerea B9, inoculation with D. ampelina was characterized by up-regulation of SA-associated genes and down-regulation of ethylene pathways, suggesting an R-gene-mediated response. With dominant effects associated with disease-free berries and minimal symptoms on canes, Rda1 and Rda2 are promising loci for grapevine genetic improvement.

Dataset Information

How many marker loci are necessary? Analysis of dominant marker data sets using two popular population genetic algorithms.

Publications

How many marker loci are necessary? Analysis of dominant marker data sets using two popular population genetic algorithms.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets