Adjustment for covariates using summary statistics of genome-wide association studies

Identification of multiple novel genetic mechanisms that regulate chilling tolerance in Arabidopsis

INTRODUCTION

Cold stress adversely affects the growth and development of plants and limits the geographical distribution of many plant species. Accumulation of spontaneous mutations shapes the adaptation of plant species to diverse climatic conditions.

METHODS

The genome-wide association study of the phenotypic variation gathered by a newly designed phenomic platform with the over six millions single nucleotide polymorphic (SNP) loci distributed across the genomes of 417 Arabidopsis natural variants collected from various geographical regions revealed 33 candidate cold responsive genes.

RESULTS

Investigation of at least two independent insertion mutants for 29 genes identified 16 chilling tolerance genes governing diverse genetic mechanisms. Five of these genes encode novel leucine-rich repeat domain-containing proteins including three nucleotide-binding site-leucine-rich repeat (NBS-LRR) proteins. Among the 16 identified chilling tolerance genes, ADS2 and ACD6 are the only two chilling tolerance genes identified earlier.

DISCUSSION

The 12.5% overlap between the genes identified in this genome-wide association study (GWAS) of natural variants with those discovered previously through forward and reverse genetic approaches suggests that chilling tolerance is a complex physiological process governed by a large number of genetic mechanisms.

Whole Exome-Sequencing of Pooled Genomic DNA Samples to Detect Quantitative Trait Loci in Esotropia and Exotropia of Strabismus in Japanese

BACKGROUND

Esotropia and exotropia are two major phenotypes of comitant strabismus. It remains controversial whether esotropia and exotropia would share common genetic backgrounds. In this study, we used a quantitative trait locus (QTL)-sequencing pipeline for diploid plants to screen for susceptibility loci of strabismus in whole exome sequencing of pooled genomic DNAs of individuals.

METHODS

Pooled genomic DNA (2.5 ng each) of 20 individuals in three groups, Japanese patients with esotropia and exotropia, and normal members in the families, was sequenced twice after exome capture, and the first and second sets of data in each group were combined to increase the read depth. The SNP index, as the ratio of variant genotype reads to all reads, and Δ(SNP index) values, as the difference of SNP index between two groups, were calculated by sliding window analysis with a 4 Mb window size and 10 kb slide size. The rows of 200 "N"s were inserted as a putative 200-b spacer between every adjoining locus to depict Δ(SNP index) plots on each chromosome. SNP positions with depth < 20 as well as SNP positions with SNP index of <0.3 were excluded.

RESULTS

After the exclusion of SNPs, 12,242 SNPs in esotropia/normal group and 12,108 SNPs in exotropia/normal group remained. The patterns of the Δ(SNP index) plots on each chromosome appeared different between esotropia/normal group and exotropia/normal group. When the consecutive groups of SNPs on each chromosome were set at three patterns: SNPs in each cytogenetic band, 50 consecutive sliding SNPs, and SNPs in 4 Mb window size with 10 kb slide size, p values (Wilcoxon signed rank test) and Q values (false discovery rate) in a few loci as Manhattan plots showed significant differences in comparison between the Δ(SNP index) in the esotropia/normal group and exotropia/normal group.

CONCLUSIONS

The pooled DNA sequencing and QTL mapping approach for plants could provide overview of genetic background on each chromosome and would suggest different genetic backgrounds for two major phenotypes of comitant strabismus, esotropia and exotropia.

The promise and deceit of genomic selection component analyses

Selection component analyses (SCA) relate individual genotype to fitness components such as viability, fecundity and mating success. SCA are based on population genetic models and yield selection estimates directly in terms of predicted allele frequency change. This paper explores the statistical properties of gSCA: experiments that apply SCA to genome-wide scoring of SNPs in field sampled individuals. Computer simulations indicate that gSCA involving a few thousand genotyped samples can detect allele frequency changes of the magnitude that has been documented in field experiments on diverse taxa. To detect selection, imprecise genotyping from low-level sequencing of large samples of individuals provides much greater power than precise genotyping of smaller samples. The simulations also demonstrate the efficacy of 'haplotype matching', a method to combine information from a limited collection of whole genome sequence (the reference panel) with the much larger sample of field individuals that are measured for fitness. Pooled sequencing is demonstrated as another way to increase statistical power. Finally, I discuss the interpretation of selection estimates in relation to the Beavis effect, the overestimation of selection intensities at significant loci.

Approximate conditional phenotype analysis based on genome wide association summary statistics

Because single genetic variants may have pleiotropic effects, one trait can be a confounder in a genome-wide association study (GWAS) that aims to identify loci associated with another trait. A typical approach to address this issue is to perform an additional analysis adjusting for the confounder. However, obtaining conditional results can be time-consuming. We propose an approximate conditional phenotype analysis based on GWAS summary statistics, the covariance between outcome and confounder, and the variant minor allele frequency (MAF). GWAS summary statistics and MAF are taken from GWAS meta-analysis results while the traits covariance may be estimated by two strategies: (i) estimates from a subset of the phenotypic data; or (ii) estimates from published studies. We compare our two strategies with estimates using individual level data from the full GWAS sample (gold standard). A simulation study for both binary and continuous traits demonstrates that our approximate approach is accurate. We apply our method to the Framingham Heart Study (FHS) GWAS and to large-scale cardiometabolic GWAS results. We observed a high consistency of genetic effect size estimates between our method and individual level data analysis. Our approach leads to an efficient way to perform approximate conditional analysis using large-scale GWAS summary statistics.