Ani Manichaikul

Robust relationship inference in genome-wide association studies

MOTIVATION Genome-wide association studies (GWASs) have been widely used to map loci contributing to variation in complex traits and risk of diseases in humans. Accurate specification of familial relationships is crucial for family-based GWAS, as well as in population-based GWAS with unknown (or unrecognized) family structure. The family structure in a GWAS should be routinely investigated using the SNP data prior to the analysis of population structure or phenotype. Existing algorithms for relationship inference have a major weakness of estimating allele frequencies at each SNP from the entire sample, under a strong assumption of homogeneous population structure. This assumption is often untenable. RESULTS Here, we present a rapid algorithm for relationship inference using high-throughput genotype data typical of GWAS that allows the presence of unknown population substructure. The relationship of any pair of individuals can be precisely inferred by robust estimation of their kinship coefficient, independent of sample composition or population structure (sample invariance). We present simulation experiments to demonstrate that the algorithm has sufficient power to provide reliable inference on millions of unrelated pairs and thousands of relative pairs (up to 3rd-degree relationships). Application of our robust algorithm to HapMap and GWAS datasets demonstrates that it performs properly even under extreme population stratification, while algorithms assuming a homogeneous population give systematically biased results. Our extremely efficient implementation performs relationship inference on millions of pairs of individuals in a matter of minutes, dozens of times faster than the most efficient existing algorithm known to us. AVAILABILITY Our robust relationship inference algorithm is implemented in a freely available software package, KING, available for download at http://people.virginia.edu/∼wc9c/KING.

Metabolic network reconstruction of Chlamydomonas offers insight into light-driven algal metabolism

Metabolic network reconstruction encompasses existing knowledge about an organisms metabolism and genome annotation, providing a platform for omics data analysis and phenotype prediction. The model alga Chlamydomonas reinhardtii is employed to study diverse biological processes from photosynthesis to phototaxis. Recent heightened interest in this species results from an international movement to develop algal biofuels. Integrating biological and optical data, we reconstructed a genome‐scale metabolic network for this alga and devised a novel light‐modeling approach that enables quantitative growth prediction for a given light source, resolving wavelength and photon flux. We experimentally verified transcripts accounted for in the network and physiologically validated model function through simulation and generation of new experimental growth data, providing high confidence in network contents and predictive applications. The network offers insight into algal metabolism and potential for genetic engineering and efficient light source design, a pioneering resource for studying light‐driven metabolism and quantitative systems biology.

Genetic loci associated with plasma phospholipid N-3 fatty acids: A Meta-Analysis of Genome-Wide association studies from the charge consortium

Long-chain n-3 polyunsaturated fatty acids (PUFAs) can derive from diet or from α-linolenic acid (ALA) by elongation and desaturation. We investigated the association of common genetic variation with plasma phospholipid levels of the four major n-3 PUFAs by performing genome-wide association studies in five population-based cohorts comprising 8,866 subjects of European ancestry. Minor alleles of SNPs in FADS1 and FADS2 (desaturases) were associated with higher levels of ALA (p = 3×10−64) and lower levels of eicosapentaenoic acid (EPA, p = 5×10−58) and docosapentaenoic acid (DPA, p = 4×10−154). Minor alleles of SNPs in ELOVL2 (elongase) were associated with higher EPA (p = 2×10−12) and DPA (p = 1×10−43) and lower docosahexaenoic acid (DHA, p = 1×10−15). In addition to genes in the n-3 pathway, we identified a novel association of DPA with several SNPs in GCKR (glucokinase regulator, p = 1×10−8). We observed a weaker association between ALA and EPA among carriers of the minor allele of a representative SNP in FADS2 (rs1535), suggesting a lower rate of ALA-to-EPA conversion in these subjects. In samples of African, Chinese, and Hispanic ancestry, associations of n-3 PUFAs were similar with a representative SNP in FADS1 but less consistent with a representative SNP in ELOVL2. Our findings show that common variation in n-3 metabolic pathway genes and in GCKR influences plasma phospholipid levels of n-3 PUFAs in populations of European ancestry and, for FADS1, in other ancestries.

Poor Performance of Bootstrap Confidence Intervals for the Location of a Quantitative Trait Locus

The aim of many genetic studies is to locate the genomic regions (called quantitative trait loci, QTL) that contribute to variation in a quantitative trait (such as body weight). Confidence intervals for the locations of QTL are particularly important for the design of further experiments to identify the gene or genes responsible for the effect. Likelihood support intervals are the most widely used method to obtain confidence intervals for QTL location, but the nonparametric bootstrap has also been recommended. Through extensive computer simulation, we show that bootstrap confidence intervals behave poorly and so should not be used in this context. The profile likelihood (or LOD curve) for QTL location has a tendency to peak at genetic markers, and so the distribution of the maximum-likelihood estimate (MLE) of QTL location has the unusual feature of point masses at genetic markers; this contributes to the poor behavior of the bootstrap. Likelihood support intervals and approximate Bayes credible intervals, on the other hand, are shown to behave appropriately.

Allelic Polymorphism within the TAS1R3 Promoter Is Associated with Human Taste Sensitivity to Sucrose

Human sweet taste perception is mediated by the heterodimeric G protein-coupled receptor encoded by the TAS1R2 and TAS1R3 genes. Variation in these genes has been characterized, but the functional consequences of such variation for sweet perception are unknown. We found that two C/T single-nucleotide polymorphisms (SNPs) located at positions -1572 (rs307355) and -1266 (rs35744813) upstream of the TAS1R3 coding sequence strongly correlate with human taste sensitivity to sucrose and explain 16% of population variability in perception. By using a luciferase reporter assay, we demonstrated that the T allele of each SNP results in reduced promoter activity in comparison to the C alleles, consistent with the phenotype observed in humans carrying T alleles. We also found that the distal region of the TAS1R3 promoter harbors a composite cis-acting element that has a strong silencing effect on promoter activity. We conclude that the rs307355 and rs35744813 SNPs affect gene transcription by altering the function of this regulatory element. A worldwide population survey reveals that the T alleles of rs307355 and rs35744813 occur at lowest frequencies in European populations. We propose that inherited differences in TAS1R3 transcription account for a substantial fraction of worldwide differences in human sweet taste perception.

Genome-Wide Association Studies Identify CHRNA5/3 and HTR4 in the Development of Airflow Obstruction

RATIONALE Genome-wide association studies (GWAS) have identified loci influencing lung function, but fewer genes influencing chronic obstructive pulmonary disease (COPD) are known. OBJECTIVES Perform meta-analyses of GWAS for airflow obstruction, a key pathophysiologic characteristic of COPD assessed by spirometry, in population-based cohorts examining all participants, ever smokers, never smokers, asthma-free participants, and more severe cases. METHODS Fifteen cohorts were studied for discovery (3,368 affected; 29,507 unaffected), and a population-based family study and a meta-analysis of case-control studies were used for replication and regional follow-up (3,837 cases; 4,479 control subjects). Airflow obstruction was defined as FEV(1) and its ratio to FVC (FEV(1)/FVC) both less than their respective lower limits of normal as determined by published reference equations. MEASUREMENTS AND MAIN RESULTS The discovery meta-analyses identified one region on chromosome 15q25.1 meeting genome-wide significance in ever smokers that includes AGPHD1, IREB2, and CHRNA5/CHRNA3 genes. The region was also modestly associated among never smokers. Gene expression studies confirmed the presence of CHRNA5/3 in lung, airway smooth muscle, and bronchial epithelial cells. A single-nucleotide polymorphism in HTR4, a gene previously related to FEV(1)/FVC, achieved genome-wide statistical significance in combined meta-analysis. Top single-nucleotide polymorphisms in ADAM19, RARB, PPAP2B, and ADAMTS19 were nominally replicated in the COPD meta-analysis. CONCLUSIONS These results suggest an important role for the CHRNA5/3 region as a genetic risk factor for airflow obstruction that may be independent of smoking and implicate the HTR4 gene in the etiology of airflow obstruction.

Genome-wide meta-analysis of observational studies shows common genetic variants associated with macronutrient intake

Background: Macronutrient intake varies substantially between individuals, and there is evidence that this variation is partly accounted for by genetic variants. Objective: The objective of the study was to identify common genetic variants that are associated with macronutrient intake. Design: We performed 2-stage genome-wide association (GWA) meta-analysis of macronutrient intake in populations of European descent. Macronutrients were assessed by using food-frequency questionnaires and analyzed as percentages of total energy consumption from total fat, protein, and carbohydrate. From the discovery GWA (n = 38,360), 35 independent loci associated with macronutrient intake at P < 5 × 10−6 were identified and taken forward to replication in 3 additional cohorts (n = 33,533) from the DietGen Consortium. For one locus, fat mass obesity-associated protein (FTO), cohorts with Illumina MetaboChip genotype data (n = 7724) provided additional replication data. Results: A variant in the chromosome 19 locus (rs838145) was associated with higher carbohydrate (β ± SE: 0.25 ± 0.04%; P = 1.68 × 10−8) and lower fat (β ± SE: −0.21 ± 0.04%; P = 1.57 × 10−9) consumption. A candidate gene in this region, fibroblast growth factor 21 (FGF21), encodes a fibroblast growth factor involved in glucose and lipid metabolism. The variants in this locus were associated with circulating FGF21 protein concentrations (P < 0.05) but not mRNA concentrations in blood or brain. The body mass index (BMI)–increasing allele of the FTO variant (rs1421085) was associated with higher protein intake (β ± SE: 0.10 ± 0.02%; P = 9.96 × 10−10), independent of BMI (after adjustment for BMI, β ± SE: 0.08 ± 0.02%; P = 3.15 × 10−7). Conclusion: Our results indicate that variants in genes involved in nutrient metabolism and obesity are associated with macronutrient consumption in humans. Trials related to this study were registered at clinicaltrials.gov as NCT00005131 (Atherosclerosis Risk in Communities), NCT00005133 (Cardiovascular Health Study), NCT00005136 (Family Heart Study), NCT00005121 (Framingham Heart Study), NCT00083369 (Genetic and Environmental Determinants of Triglycerides), NCT01331512 (InCHIANTI Study), and NCT00005487 (Multi-Ethnic Study of Atherosclerosis).

A Model Selection Approach for the Identification of Quantitative Trait Loci in Experimental Crosses, Allowing Epistasis

The identification of quantitative trait loci (QTL) and their interactions is a crucial step toward the discovery of genes responsible for variation in experimental crosses. The problem is best viewed as one of model selection, and the most important aspect of the problem is the comparison of models of different sizes. We present a penalized likelihood approach, with penalties on QTL and pairwise interactions chosen to control false positive rates. This extends the work of Broman and Speed to allow for pairwise interactions among QTL. A conservative version of our penalized LOD score provides strict control over the rate of extraneous QTL and interactions; a more liberal criterion is more lenient on interactions but seeks to maintain control over the rate of inclusion of false loci. The key advance is that one needs only to specify a target false positive rate rather than a prior on the number of QTL and interactions. We illustrate the use of our model selection criteria as exploratory tools; simulation studies demonstrate reasonable power to detect QTL. Our liberal criterion is comparable in power to two Bayesian approaches.

Metabolic network analysis integrated with transcript verification for sequenced genomes

With sequencing of thousands of organisms completed or in progress, there is a growing need to integrate gene prediction with metabolic network analysis. Using Chlamydomonas reinhardtii as a model, we describe a systems-level methodology bridging metabolic network reconstruction with experimental verification of enzyme encoding open reading frames. Our quantitative and predictive metabolic model and its associated cloned open reading frames provide useful resources for metabolic engineering.

The X chromosome in quantitative trait locus mapping.

The X chromosome requires special treatment in the mapping of quantitative trait loci (QTL). However, most QTL mapping methods, and most computer programs for QTL mapping, have focused exclusively on autosomal loci. We describe a method for appropriate treatment of the X chromosome for QTL mapping in experimental crosses. We address the important issue of formulating the null hypothesis of no linkage appropriately. If the X chromosome is treated like an autosome, a sex difference in the phenotype can lead to spurious linkage on the X chromosome. Further, the number of degrees of freedom for the linkage test may be different for the X chromosome than for autosomes, and so an X chromosome-specific significance threshold is required. To address this issue, we propose a general procedure to obtain chromosome-specific significance thresholds that controls the genomewide false positive rate at the desired level. We apply our methods to data on gut length in a large intercross of mice carrying the Sox10Dom mutation, a model of Hirschsprung disease. We identified QTL contributing to variation in gut length on chromosomes 5 and 18. We found suggestive evidence of linkage to the X chromosome, which would be viewed as strong evidence of linkage if the X chromosome was treated as an autosome. Our methods have been implemented in the package R/qtl.