Quantitative trait loci, G×E and G×G for glycemic traits: response to metformin and placebo in the Diabetes Prevention Program (DPP).

The complex genetic architecture of type-2-diabetes (T2D) includes gene-by-environment (G×E) and gene-by-gene (G×G) interactions. To identify G×E and G×G, we screened markers for patterns indicative of interactions (relationship loci [rQTL] and variance heterogeneity loci [vQTL]). rQTL exist when the correlation between multiple traits varies by genotype and vQTL occur when the variance of a trait differs by genotype (potentially flagging G×G and G×E). In the metformin and placebo arms of the DPP (n = 1762) we screened 280,965 exomic and intergenic SNPs, for rQTL and vQTL patterns in association with year one changes from baseline in glycemia and related traits (insulinogenic index [IGI], insulin sensitivity index [ISI], fasting glucose and fasting insulin). Significant (p < 1.8 × 10-7) rQTL and vQTL generated a priori hypotheses of individual G×E tests for a SNP × metformin treatment interaction and secondarily for G×G screens. Several rQTL and vQTL identified led to 6 nominally significant (p < 0.05) metformin treatment × SNP interactions (4 for IGI, one insulin, and one glucose) and 12G×G interactions (all IGI) that exceeded experiment-wide significance (p < 4.1 × 10-9). Some loci are directly associated with incident diabetes, and others are rQTL and modify a trait's relationship with diabetes (2 diabetes/glucose, 2 diabetes/insulin, 1 diabetes/IGI). rs3197999, an ISI/insulin rQTL, is a possible gene damaging missense mutation in MST1, is associated with ulcerative colitis, sclerosing cholangitis, Crohn's disease, BMI and coronary artery disease. This study demonstrates evidence for context-dependent effects (G×G & G×E) and the complexity of these T2D-related traits.

A versatile omnibus test for detecting mean and variance heterogeneity.

Recent research has revealed loci that display variance heterogeneity through various means such as biological disruption, linkage disequilibrium (LD), gene-by-gene (G × G), or gene-by-environment interaction. We propose a versatile likelihood ratio test that allows joint testing for mean and variance heterogeneity (LRT(MV)) or either effect alone (LRT(M) or LRT(V)) in the presence of covariates. Using extensive simulations for our method and others, we found that all parametric tests were sensitive to nonnormality regardless of any trait transformations. Coupling our test with the parametric bootstrap solves this issue. Using simulations and empirical data from a known mean-only functional variant, we demonstrate how LD can produce variance-heterogeneity loci (vQTL) in a predictable fashion based on differential allele frequencies, high D', and relatively low r² values. We propose that a joint test for mean and variance heterogeneity is more powerful than a variance-only test for detecting vQTL. This takes advantage of loci that also have mean effects without sacrificing much power to detect variance only effects. We discuss using vQTL as an approach to detect G × G interactions and also how vQTL are related to relationship loci, and how both can create prior hypothesis for each other and reveal the relationships between traits and possibly between components of a composite trait.

A family-based joint test for mean and variance heterogeneity for quantitative traits.

Traditional quantitative trait locus (QTL) analysis focuses on identifying loci associated with mean heterogeneity. Recent research has discovered loci associated with phenotype variance heterogeneity (vQTL), which is important in studying genetic association with complex traits, especially for identifying gene-gene and gene-environment interactions. While several tests have been proposed to detect vQTL for unrelated individuals, there are no tests for related individuals, commonly seen in family-based genetic studies. Here we introduce a likelihood ratio test (LRT) for identifying mean and variance heterogeneity simultaneously or for either effect alone, adjusting for covariates and family relatedness using a linear mixed effect model approach. The LRT test statistic for normally distributed quantitative traits approximately follows χ(2)-distributions. To correct for inflated Type I error for non-normally distributed quantitative traits, we propose a parametric bootstrap-based LRT that removes the best linear unbiased prediction (BLUP) of family random effect. Simulation studies show that our family-based test controls Type I error and has good power, while Type I error inflation is observed when family relatedness is ignored. We demonstrate the utility and efficiency gains of the proposed method using data from the Framingham Heart Study to detect loci associated with body mass index (BMI) variability.

Mitochondrial genomic variation associated with higher mitochondrial copy number: the Cache County Study on Memory Health and Aging.

BACKGROUND: The mitochondria are essential organelles and are the location of cellular respiration, which is responsible for the majority of ATP production. Each cell contains multiple mitochondria, and each mitochondrion contains multiple copies of its own circular genome. The ratio of mitochondrial genomes to nuclear genomes is referred to as mitochondrial copy number. Decreases in mitochondrial copy number are known to occur in many tissues as people age, and in certain diseases. The regulation of mitochondrial copy number by nuclear genes has been studied extensively. While mitochondrial variation has been associated with longevity and some of the diseases known to have reduced mitochondrial copy number, the role that the mitochondrial genome itself has in regulating mitochondrial copy number remains poorly understood. RESULTS: We analyzed the complete mitochondrial genomes from 1007 individuals randomly selected from the Cache County Study on Memory Health and Aging utilizing the inferred evolutionary history of the mitochondrial haplotypes present in our dataset to identify sequence variation and mitochondrial haplotypes associated with changes in mitochondrial copy number. Three variants belonging to mitochondrial haplogroups U5A1 and T2 were significantly associated with higher mitochondrial copy number in our dataset. CONCLUSIONS: We identified three variants associated with higher mitochondrial copy number and suggest several hypotheses for how these variants influence mitochondrial copy number by interacting with known regulators of mitochondrial copy number. Our results are the first to report sequence variation in the mitochondrial genome that causes changes in mitochondrial copy number. The identification of these variants that increase mtDNA copy number has important implications in understanding the pathological processes that underlie these phenotypes.

Estimating population genetic parameters and comparing model goodness-of-fit using DNA sequences with error.

It is known that sequencing error can bias estimation of evolutionary or population genetic parameters. This problem is more prominent in deep resequencing studies because of their large sample size n, and a higher probability of error at each nucleotide site. We propose a new method based on the composite likelihood of the observed SNP configurations to infer population mutation rate theta = 4N(e)micro, population exponential growth rate R, and error rate epsilon, simultaneously. Using simulation, we show the combined effects of the parameters, theta, n, epsilon, and R on the accuracy of parameter estimation. We compared our maximum composite likelihood estimator (MCLE) of theta with other theta estimators that take into account the error. The results show the MCLE performs well when the sample size is large or the error rate is high. Using parametric bootstrap, composite likelihood can also be used as a statistic for testing the model goodness-of-fit of the observed DNA sequences. The MCLE method is applied to sequence data on the ANGPTL4 gene in 1832 African American and 1045 European American individuals.

The importance of context to the genetic architecture of diabetes-related traits is revealed in a genome-wide scan of a LG/J × SM/J murine model.

Variations in diabetic phenotypes are caused by complex interactions of genetic effects, environmental factors, and the interplay between the two. We tease apart these complex interactions by examining genome-wide genetic and epigenetic effects on diabetes-related traits among different sex, diet, and sex-by-diet cohorts in a Mus musculus model. We conducted a genome-wide scan for quantitative trait loci that affect serum glucose and insulin levels and response to glucose stress in an F(16) Advanced Intercross Line of the LG/J and SM/J intercross (Wustl:LG,SM-G16). Half of each sibship was fed a high-fat diet and half was fed a relatively low-fat diet. Context-dependent genetic (additive and dominance) and epigenetic (parent-of-origin imprinting) effects were characterized by partitioning animals into sex, diet, and sex-by-diet cohorts. We found that different cohorts often have unique genetic effects at the same loci, and that genetic signals can be masked or erroneously assigned to specific cohorts if they are not considered individually. Our data demonstrate that the effects of genes on complex trait variation are highly context-dependent and that the same genomic sequence can affect traits differently depending on an individual's sex and/or dietary environment. Our results have important implications for studies of complex traits in humans.

Genetic, epigenetic, and gene-by-diet interaction effects underlie variation in serum lipids in a LG/JxSM/J murine model.

Variation in serum cholesterol, free-fatty acids, and triglycerides is associated with cardiovascular disease (CVD) risk factors. There is great interest in characterizing the underlying genetic architecture of these risk factors, because they vary greatly within and among human populations and between the sexes. We present results of a genome-wide scan for quantitative trait loci (QTL) affecting serum cholesterol, free-fatty acids, and triglycerides in an F(16) advanced intercross line of LG/J and SM/J (Wustl:LG,SM-G16). Half of the population was fed a high-fat diet and half was fed a relatively low-fat diet. Context-dependent genetic (additive and dominance) and epigenetic (imprinting) effects were characterized by partitioning animals into sex, diet, and sex-by-diet cohorts. Here we examine genetic, environmental, and genetic-by-environmental interactions of QTL overlapping previously identified loci associated with CVD risk factors, and we add to the serum lipid QTL landscape by identifying new loci.

Inferring population mutation rate and sequencing error rate using the SNP frequency spectrum in a sample of DNA sequences.

One challenge of analyzing samples of DNA sequences is to account for the nonnegligible polymorphisms produced by error when the sequencing error rate is high or the sample size is large. Specifically, those artificial sequence variations will bias the observed single nucleotide polymorphism (SNP) frequency spectrum, which in turn may further bias the estimators of the population mutation rate theta =4N mu for diploids. In this paper, we propose a new approach based on the generalized least squares (GLS) method to estimate theta, given a SNP frequency spectrum in a random sample of DNA sequences from a population. With this approach, error rate epsilon can be either known or unknown. In the latter case, epsilon can be estimated given an estimation of theta. Using coalescent simulation, we compared our estimators with other estimators of theta. The results showed that the GLS estimators are more efficient than other theta estimators with error, and the estimation of epsilon is usable in practice when the theta per bp is small. We demonstrate the application of the estimators with 10-kb noncoding region sequence sampled from a human population and provide suggestions for choosing theta estimators with error.

Quantitative trait loci associated with susceptibility to therapy-related acute murine promyelocytic leukemia in hCG-PML/RARA transgenic mice.

Therapy-related acute myelogenous leukemia (t-AML) is an important late adverse effect of alkylator chemotherapy. Susceptibility to t-AML has a genetic component, yet specific genetic variants that influence susceptibility are poorly understood. We analyzed an F(2) intercross (n = 282 mice) between mouse strains resistant or susceptible to t-AML induced by the alkylator ethyl-N-nitrosourea (ENU) to identify genes that regulate t-AML susceptibility. Each mouse carried the hCG-PML/RARA transgene, a well-characterized initiator of myeloid leukemia. In the absence of ENU treatment, transgenic F(2) mice developed leukemia with higher incidence (79.4% vs 12.5%) and at earlier time points (108 days vs 234 days) than mice in the resistant background. ENU treatment of F(2) mice further increased incidence (90.4%) and shortened median survival (171 vs 254 days). We genotyped F(2) mice at 384 informative single nucleotide polymorphisms across the genome and performed quantitative trait locus (QTL) analysis. Thirteen QTLs significantly associated with leukemia-free survival, spleen weight, or white blood cell count were identified on 8 chromosomes. These results suggest that susceptibility to ENU-induced leukemia in mice is a complex trait governed by genes at multiple loci. Improved understanding of genetic risk factors should lead to tailored treatment regimens that reduce risk for patients predisposed to t-AML.

Phylogenetics applied to genotype/phenotype association and selection analyses with sequence data from angptl4 in humans.

Genotype/phenotype association analyses (Treescan) with plasma lipid levels and functional site prediction methods (TreeSAAP and PolyPhen) were performed using sequence data for ANGPTL4 from 3,551 patients in the Dallas Heart Study. Biological assays of rare variants in phenotypic tails and results from a Treescan analysis were used as "known" variants to assess the site prediction abilities of PolyPhen and TreeSAAP. The E40K variant in European Americans and the R278Q variant in African Americans were significantly associated with multiple lipid phenotypes. Combining TreeSAAP and PolyPhen performed well to predict "known" functional variants while reducing noise from false positives.

Replication of long-bone length QTL in the F9-F10 LG,SM advanced intercross.

Quantitative trait locus (QTL) mapping techniques are frequently used to identify genomic regions associated with variation in phenotypes of interest. However, the F(2) intercross and congenic strain populations usually employed have limited genetic resolution resulting in relatively large confidence intervals that greatly inhibit functional confirmation of statistical results. Here we use the increased resolution of the combined F(9) and F(10) generations (n = 1455) of the LG,SM advanced intercross to fine-map previously identified QTL associated with the lengths of the humerus, ulna, femur, and tibia. We detected 81 QTL affecting long-bone lengths. Of these, 49 were previously identified in the combined F(2)-F(3) population of this intercross, while 32 represent novel contributors to trait variance. Pleiotropy analysis suggests that most QTL affect three to four long bones or serially homologous limb segments. We also identified 72 epistatic interactions involving 38 QTL and 88 novel regions. This analysis shows that using later generations of an advanced intercross greatly facilitates fine-mapping of confidence intervals, resolving three F(2)-F(3) QTL into multiple linked loci and narrowing confidence intervals of other loci, as well as allowing identification of additional QTL. Further characterization of the biological bases of these QTL will help provide a better understanding of the genetics of small variations in long-bone length.

Irinotecan pharmacogenetics: influence of pharmacodynamic genes.

PURPOSE: Irinotecan is an important drug for the treatment of solid tumors. Although genes involved in irinotecan pharmacokinetics have been shown to influence toxicity, there are no data on pharmacodynamic genes. CDC45L, NFKB1, PARP1, TDP1, and XRCC1 have been shown to influence the cytotoxic action of camptothecins, including irinotecan. Polymorphisms in the drug target of camptothecins, topoisomerase I (TOP1), and downstream effectors may influence patient outcomes to irinotecan therapy. We undertook a retrospective candidate gene haplotype association study to investigate this hypothesis. EXPERIMENTAL DESIGN: Haplotype compositions of six candidate genes were constructed in European (n = 93), East Asian (n = 94), and West African (n = 95) populations. Haplotype-tagging single nucleotide polymorphisms (htSNP) were selected based on genealogic relationships between haplotypes. DNA samples from 107 European, advanced colorectal cancer patients treated with irinotecan-based regimens were genotyped for htSNPs as well as three coding region SNPs. Associations between genetic variants and toxicity (grade 3/4 diarrhea and neutropenia) or efficacy (objective response) were assessed. RESULTS: TOP1 and TDP1 htSNPs were related to grade 3/4 neutropenia (P = 0.04) and response (P = 0.04), respectively. Patients homozygous for an XRCC1 haplotype (GGCC-G) were more likely to show an objective response to therapy than other patients (83% versus 30%; P = 0.02). This effect was also seen in a multivariate analysis (odds ratio, 11.9; P = 0.04). No genetic variants were associated with diarrhea. CONCLUSIONS: This is the first comprehensive pharmacogenetic investigation of irinotecan pharmacodynamic factors, and our findings suggest that genetic variation in the pharmacodynamic genes may influence the efficacy of irinotecan-containing therapies in advanced colorectal cancer patients.

Genome-wide association study for variants that modulate relationships between cerebrospinal fluid amyloid-beta 42, tau, and p-tau levels.

BACKGROUND: A relationship quantitative trait locus exists when the correlation between multiple traits varies by genotype for that locus. Relationship quantitative trait loci (rQTL) are often involved in gene-by-gene (G×G) interactions or gene-by-environmental interactions, making them a powerful tool for detecting G×G. METHODS: We performed genome-wide association studies to identify rQTL between tau and Aß42 and ptau and Aß42 with over 3000 individuals using age, gender, series, APOE ε2, APOE ε4, and two principal components for population structure as covariates. Each significant rQTL was separately screened for interactions with other loci for each trait in the rQTL model. Parametric bootstrapping was used to assess significance. RESULTS: We found four significant tau/Aß42 rQTL from three unique locations and six ptau/Aß42 rQTL from five unique locations. G×G screens with these rQTL produced four significant G×G interactions (one Aß42, two ptau, and one tau) with four rQTL where each second locus was from a unique location. On follow-up, rs1036819 and rs74025622 were associated with Alzheimer's disease (AD) case/control status; rs15205 and rs79099429 were associated with rate of decline. CONCLUSIONS: The two most significant rQTL (rs8027714 and rs1036819) for ptau/Aß42 are on different chromosomes and both are strong hits for pelvic organ prolapse. While diseases of the nervous system can cause pelvic organ prolapse, it is unlikely related to the ptau/Aß42 relationship but may suggest that these two loci share a pathway. In addition to a ptau/Aß42 rQTL and association with AD case/control status, rs1036819 is a strong rQTL for case/control status/Aß42 and for tau/Aß42. It resides in the ZFAT gene, which is related to autoimmune thyroid disease. For tau, rs9817620 interacts with the tau/Aß42 rQTL rs74025622. It is in the CHL1 gene, which is a neural cell adhesion molecule and may be involved in signal transduction pathways. CHL1 is related to BACE1, which is a ß-secretase enzyme that initiates production of the ß-amyloid peptide involved in AD and is a primary drug target. Overall, there are numerous loci that affect the relationship between these important AD endophenotypes and some are due to interactions with other loci. Some affect the risk of AD and/or rate of progression.

Characterization of HIV diversity, phylodynamics and drug resistance in Washington, DC.

BACKGROUND: Washington DC has a high burden of HIV with a 2.0% HIV prevalence. The city is a national and international hub potentially containing a broad diversity of HIV variants; yet few sequences from DC are available on GenBank to assess the evolutionary history of HIV in the US capital. Towards this general goal, here we analyze extensive sequence data and investigate HIV diversity, phylodynamics, and drug resistant mutations (DRM) in DC. METHODS: Molecular HIV-1 sequences were collected from participants infected through 2015 as part of the DC Cohort, a longitudinal observational study of HIV+ patients receiving care at 13 DC clinics. Sequences were paired with Cohort demographic, risk, and clinical data and analyzed using maximum likelihood, Bayesian and coalescent approaches of phylogenetic, network and population genetic inference. We analyzed 601 sequences from 223 participants for int (~864 bp) and 2,810 sequences from 1,659 participants for PR/RT (~1497 bp). RESULTS: Ninety-nine and 94% of the int and PR/RT sequences, respectively, were identified as subtype B, with 14 non-B subtypes also detected. Phylodynamic analyses of US born infected individuals showed that HIV population size varied little over time with no significant decline in diversity. Phylogenetic analyses grouped 13.5% of the int sequences into 14 clusters of 2 or 3 sequences, and 39.0% of the PR/RT sequences into 203 clusters of 2-32 sequences. Network analyses grouped 3.6% of the int sequences into 4 clusters of 2 sequences, and 10.6% of the PR/RT sequences into 76 clusters of 2-7 sequences. All network clusters were detected in our phylogenetic analyses. Higher proportions of clustered sequences were found in zip codes where HIV prevalence is highest (r = 0.607; P<0.00001). We detected a high prevalence of DRM for both int (17.1%) and PR/RT (39.1%), but only 8 int and 12 PR/RT amino acids were identified as under adaptive selection. We observed a significant (P<0.0001) association between main risk factors (men who have sex with men and heterosexuals) and genotypes in the five well-supported clusters with sufficient sample size for testing. DISCUSSION: Pairing molecular data with clinical and demographic data provided novel insights into HIV population dynamics in Washington, DC. Identification of populations and geographic locations where clustering occurs can inform and complement active surveillance efforts to interrupt HIV transmission.

Tree scanning: a method for using haplotype trees in phenotype/genotype association studies.

We use evolutionary trees of haplotypes to study phenotypic associations by exhaustively examining all possible biallelic partitions of the tree, a technique we call tree scanning. If the first scan detects significant associations, additional rounds of tree scanning are used to partition the tree into three or more allelic classes. Two worked examples are presented. The first is a reanalysis of associations between haplotypes at the Alcohol Dehydrogenase locus in Drosophila melanogaster that was previously analyzed using a nested clade analysis, a more complicated technique for using haplotype trees to detect phenotypic associations. Tree scanning and the nested clade analysis yield the same inferences when permutation testing is used with both approaches. The second example is an analysis of associations between variation in various lipid traits and genetic variation at the Apolipoprotein E (APOE) gene in three human populations. Tree scanning successfully identified phenotypic associations expected from previous analyses. Tree scanning for the most part detected more associations and provided a better biological interpretative framework than single SNP analyses. We also show how prior information can be incorporated into the tree scan by starting with the traditional three electrophoretic alleles at APOE. Tree scanning detected genetically determined phenotypic heterogeneity within all three electrophoretic allelic classes. Overall, tree scanning is a simple, powerful, and flexible method for using haplotype trees to detect phenotype/genotype associations at candidate loci.

Beta-2 adrenergic receptor genotypes and haplotypes in different ethnic groups.

The human beta-2 adrenergic receptor (beta2AR) is responsible for the binding of endogenous catecholamines and their exogenously administered agonists and antagonists. Three functional polymorphisms in codons 16, 27 and 164 have been described which have clinical importance for several diseases, including asthma, hypertension, heart failure, cystic fibrosis and obesity, as well as response to beta-agonist therapy. These were evaluated in 726 individuals from 8 distinct ethnic populations (Chinese, Filipino, Southwest Asian, Saudi, Ghanaian, Kenyan, Sudanese, and European from Scotland). The results show that most haplotypes are shared among all populations, yet there are marked differences in their frequency distributions geographically. The genetic distance tree is different from standard human population distance trees, implying a different mode of evolution for this locus than that for human population gene-flow history. The multilocus frequency differences between the observed clusters of populations correspond to historical haplotype groupings that have been found to be functionally different with respect to multiple medically related phenotypes. Further studies are needed to see if functional relationships are the same across populations.

Mitochondrial haplotypes associated with biomarkers for Alzheimer's disease.

Various studies have suggested that the mitochondrial genome plays a role in late-onset Alzheimer's disease, although results are mixed. We used an endophenotype-based approach to further characterize mitochondrial genetic variation and its relationship to risk markers for Alzheimer's disease. We analyzed longitudinal data from non-demented, mild cognitive impairment, and late-onset Alzheimer's disease participants in the Alzheimer's Disease Neuroimaging Initiative with genetic, brain imaging, and behavioral data. We assessed the relationship of structural MRI and cognitive biomarkers with mitochondrial genome variation using TreeScanning, a haplotype-based approach that concentrates statistical power by analyzing evolutionarily meaningful groups (or clades) of haplotypes together for association with a phenotype. Four clades were associated with three different endophenotypes: whole brain volume, percent change in temporal pole thickness, and left hippocampal atrophy over two years. This is the first study of its kind to identify mitochondrial variation associated with brain imaging endophenotypes of Alzheimer's disease. Our results provide additional evidence that the mitochondrial genome plays a role in risk for Alzheimer's disease.

APOE modulates the correlation between triglycerides, cholesterol, and CHD through pleiotropy, and gene-by-gene interactions.

Relationship loci (rQTL) exist when the correlation between multiple traits varies by genotype. rQTL often occur due to gene-by-gene (G × G) or gene-by-environmental interactions, making them a powerful tool for detecting G × G. Here we present an empirical analysis of apolipoprotein E (APOE) with respect to lipid traits and incident CHD leading to the discovery of loci that interact with APOE to affect these traits. We found that the relationship between total cholesterol (TC) and triglycerides (ln TG) varies by APOE isoform genotype in African-American (AA) and European-American (EA) populations. The e2 allele is associated with strong correlation between ln TG and TC while the e4 allele leads to little or no correlation. This led to a priori hypotheses that APOE genotypes affect the relationship of TC and/or ln TG with incident CHD. We found that APOE*TC was significant (P = 0.016) for AA but not EA while APOE*ln TG was significant for EA (P = 0.027) but not AA. In both cases, e2e2 and e2e3 had strong relationships between TC and ln TG with CHD while e2e4 and e4e4 results in little or no relationship between TC and ln TG with CHD. Using ARIC GWAS data, scans for loci that significantly interact with APOE produced four loci for African Americans (one CHD, one TC, and two HDL). These interactions contribute to the rQTL pattern. rQTL are a powerful tool to identify loci that modify the relationship between risk factors and disease and substantially increase statistical power for detecting G × G.