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: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.

Project description:BackgroundThe choice of probe set algorithms for expression summary in a GeneChip study has a great impact on subsequent gene expression data analysis. Spiked-in cRNAs with known concentration are often used to assess the relative performance of probe set algorithms. Given the fact that the spiked-in cRNAs do not represent endogenously expressed genes in experiments, it becomes increasingly important to have methods to study whether a particular probe set algorithm is more appropriate for a specific dataset, without using such external reference data.ResultsWe propose the use of the probe set redundancy feature for evaluating the performance of probe set algorithms, and have presented three approaches for analyzing data variance and result bias using two sample t-test statistics from redundant probe sets. These approaches are as follows: 1) analyzing redundant probe set variance based on t-statistic rank order, 2) computing correlation of t-statistics between redundant probe sets, and 3) analyzing the co-occurrence of replicate redundant probe sets representing differentially expressed genes. We applied these approaches to expression summary data generated from three datasets utilizing individual probe set algorithms of MAS5.0, dChip, or RMA. We also utilized combinations of options from the three probe set algorithms. We found that results from the three approaches were similar within each individual expression summary dataset, and were also in good agreement with previously reported findings by others. We also demonstrate the validity of our findings by independent experimental methods.ConclusionAll three proposed approaches allowed us to assess the performance of probe set algorithms using the probe set redundancy feature. The analyses of redundant probe set variance based on t-statistic rank order and correlation of t-statistics between redundant probe sets provide useful tools for data variance analysis, and the co-occurrence of replicate redundant probe sets representing differentially expressed genes allows estimation of result bias. The results also suggest that individual probe set algorithms have dataset-specific performance.

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