Dataset Information

Extreme sampling design in genetic association mapping of quantitative trait loci using balanced and unbalanced case-control samples.

ABSTRACT: It is extremely expensive to conduct large sample size array- or sequencing based genome scale association studies. For a quantitative trait, an extreme case-control study design may improve the power and reduce the cost of variant calling. We investigated the performance of extreme study design when various proportions of samples are selected from the tails of phenotype distribution. Using simulations, we show that when risk genotypes become rare in the population and effect size is relatively small, it is beneficial to carry out an extreme sampling study. In particular, the number of selected cases and controls can even be unbalanced such that power is further increased, compared with a balanced selection. Our application to two data sets: methadone dose data and yearling weight data, demonstrated that similar results for full data analysis can be obtained using extreme sampling with only a fraction of the data. Using power analysis with simulated data and an experimental data application, we conclude that when full data is unavailable due to restricted budget, it is rewarding to employ an extreme sampling design in the sense that there can be immense cost reductions and qualitatively similar power as in the full data analysis.

SUBMITTER: Li Y

PROVIDER: S-EPMC6820758 | biostudies-literature | 2019 Oct

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Extreme sampling design in genetic association mapping of quantitative trait loci using balanced and unbalanced case-control samples.

Li Yi Y Levran Orna O Kim JongJoo J Zhang Tiejun T Chen Xingdong X Suo Chen C

Scientific reports 20191029 1

It is extremely expensive to conduct large sample size array- or sequencing based genome scale association studies. For a quantitative trait, an extreme case-control study design may improve the power and reduce the cost of variant calling. We investigated the performance of extreme study design when various proportions of samples are selected from the tails of phenotype distribution. Using simulations, we show that when risk genotypes become rare in the population and effect size is relatively ...[more]

PMID: 31664079

Similar Datasets

Project description:BackgroundSingle nucleotide polymorphisms (SNPs) which capture a significant impact on a trait can be identified with genome-wide association studies. High linkage disequilibrium (LD) among SNPs makes it difficult to identify causative variants correctly. Thus, often target regions instead of single SNPs are reported. Sample size has not only a crucial impact on the precision of parameter estimates, it also ensures that a desired level of statistical power can be reached. We study the design of experiments for fine-mapping of signals of a quantitative trait locus in such a target region.MethodsA multi-locus model allows to identify causative variants simultaneously, to state their positions more precisely and to account for existing dependencies. Based on the commonly applied SNP-BLUP approach, we determine the z-score statistic for locally testing non-zero SNP effects and investigate its distribution under the alternative hypothesis. This quantity employs the theoretical instead of observed dependence between SNPs; it can be set up as a function of paternal and maternal LD for any given population structure.ResultsWe simulated multiple paternal half-sib families and considered a target region of 1 Mbp. A bimodal distribution of estimated sample size was observed, particularly if more than two causative variants were assumed. The median of estimates constituted the final proposal of optimal sample size; it was consistently less than sample size estimated from single-SNP investigation which was used as a baseline approach. The second mode pointed to inflated sample sizes and could be explained by blocks of varying linkage phases leading to negative correlations between SNPs. Optimal sample size increased almost linearly with number of signals to be identified but depended much stronger on the assumption on heritability. For instance, three times as many samples were required if heritability was 0.1 compared to 0.3. An R package is provided that comprises all required tools.ConclusionsOur approach incorporates information about the population structure into the design of experiments. Compared to a conventional method, this leads to a reduced estimate of sample size enabling the resource-saving design of future experiments for fine-mapping of candidate variants.

Project description:A new mixed-model method was developed for mapping quantitative trait loci (QTL) by incorporating multiple polygenic covariance structures. First, we used genome-wide markers to calculate six different kinship matrices. We then partitioned the total genetic variance into six variance components, one corresponding to each kinship matrix, including the additive, dominance, additive × additive, dominance × dominance, additive × dominance, and dominance × additive variances. The six different kinship matrices along with the six estimated polygenic variances were used to control the genetic background of a QTL mapping model. Simulation studies showed that incorporating epistatic polygenic covariance structure can improve QTL mapping resolution. The method was applied to yield component traits of rice. We analyzed four traits (yield, tiller number, grain number, and grain weight) using 278 immortal F2 crosses (crosses between recombinant inbred lines) and 1619 markers. We found that the relative importance of each type of genetic variance varies across different traits. The total genetic variance of yield is contributed by additive × additive (18%), dominance × dominance (14%), additive × dominance (48%), and dominance × additive (15%) variances. Tiller number is contributed by additive (17%), additive × additive (22%), and dominance × additive (43%) variances. Grain number is mainly contributed by additive (42%), additive × additive (19%), and additive × dominance (31%) variances. Grain weight is almost exclusively contributed by the additive (73%) variance plus a small contribution from the additive × additive (10%) variance. Using the estimated genetic variance components to capture the polygenic covariance structure, we detected 39 effects for yield, 39 effects for tiller number, 24 for grain number, and 15 for grain weight. The new method can be directly applied to polygenic-effect-adjusted genome-wide association studies (GWAS) in human and other species.

Dataset Information

Extreme sampling design in genetic association mapping of quantitative trait loci using balanced and unbalanced case-control samples.

Publications

Extreme sampling design in genetic association mapping of quantitative trait loci using balanced and unbalanced case-control samples.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets