X chromosome effects and their interactions with mitochondrial effects

We report a simple and rapid method for detecting additive genetic variance due to X-linked loci in the absence of marker data for this chromosome. We examined the interaction of this method with an established method for detecting mitochondrial linkage (another source of sex-asymmetric genetic covariance). When applied to data from the Collaborative Study on the Genetics of Alcoholism, this method found evidence of X-chromosomal linkage for one continuous trait (ntth1) and one discrete trait (SPENT). Evidence of mitochondrial contribution was found for one discrete trait (CRAVING) and three continuous traits (ln(CIGPKYR), ecb21, and tth1). Results for ntth1 suggest that methods that do not also allow for male-female heterogeneity in environmental variance may be overly conservative in detection of X-chromosomal effects

Methylation of SOCS3 is inversely associated with metabolic syndrome in an epigenome-wide association study of obesity

Epigenetic mechanisms, including DNA methylation, mediate the interaction between gene and environment and may play an important role in the obesity epidemic. We assessed the relationship between DNA methylation and obesity in peripheral blood mononuclear cells (PBMCs) at 485,000 CpG sites across the genome in family members (8-90 y of age) using a discovery cohort (192 individuals) and a validation cohort (1,052 individuals) of Northern European ancestry. After Bonferroni-correction (Pα=0.05 = 1.31 × 10-7) for genome-wide significance, we identified 3 loci, cg18181703 (SOCS3), cg04502490 (ZNF771), and cg02988947 (LIMD2), where methylation status was associated with body mass index percentile (BMI%), a clinical index for obesity in children, adolescents, and adults. These sites were also associated with multiple metabolic syndrome (MetS) traits, including central obesity, fat depots, insulin responsiveness, and plasma lipids. The SOCS3 methylation locus was also associated with the clinical definition of MetS. In the validation cohort, SOCS3 methylation status was found to be inversely associated with BMI% (P = 1.75 × 10-6), waist to height ratio (P = 4.18 × 10-7), triglycerides (P = 4.01 × 10-4), and MetS (P = 4.01 × 10-7), and positively correlated with HDL-c (P = 4.57 × 10-8). Functional analysis in a sub cohort (333 individuals) demonstrated SOCS3 methylation and gene expression in PBMCs were inversely correlated (P = 2.93 × 10-4) and expression of SOCS3 was positively correlated with status of MetS (P = 0.012). We conclude that epigenetic modulation of SOCS3, a gene involved in leptin and insulin signaling, may play an important role in obesity and MetS

Genetic Analysis Workshop 17 mini-exome simulation

The data set simulated for Genetic Analysis Workshop 17 was designed to mimic a subset of data that might be produced in a full exome screen for a complex disorder and related risk factors in order to permit workshop participants to investigate issues of study design and statistical genetic analysis. Real sequence data from the 1000 Genomes Project formed the basis for simulating a common disease trait with a prevalence of 30% and three related quantitative risk factors in a sample of 697 unrelated individuals and a second sample of 697 individuals in large, extended pedigrees. Called genotypes for 24,487 autosomal markers assigned to 3,205 genes and simulated affection status, quantitative traits, age, sex, pedigree relationships, and cigarette smoking were provided to workshop participants. The simulating model included both common and rare variants with minor allele frequencies ranging from 0.07% to 25.8% and a wide range of effect sizes for these variants. Genotype-smoking interaction effects were included for variants in one gene. Functional variants were concentrated in genes selected from specific biological pathways and were selected on the basis of the predicted deleteriousness of the coding change. For each sample, unrelated individuals and family, 200 replicates of the phenotypes were simulated

Genetic basis of neurocognitive decline and reduced white-matter integrity in normal human brain aging

Identification of genes associated with brain aging should markedly improve our understanding of the biological processes that govern normal age-related decline. However, challenges to identifying genes that facilitate successful brain aging are considerable, including a lack of established phenotypes and difficulties in modeling the effects of aging per se, rather than genes that influence the underlying trait. In a large cohort of randomly selected pedigrees (n = 1,129 subjects), we documented profound aging effects from young adulthood to old age (18-83 y) on neurocognitive ability and diffusion-based white-matter measures. Despite significant phenotypic correlation between white-matter integrity and tests of processing speed, working memory, declarative memory, and intelligence, no evidence for pleiotropy between these classes of phenotypes was observed. Applying an advanced quantitative gene-by-environment interaction analysis where age is treated as an environmental factor, we demonstrate a heritable basis for neurocognitive deterioration as a function of age. Furthermore, by decomposing gene-by-aging (G × A) interactions, we infer that different genes influence some neurocognitive traits as a function of age, whereas other neurocognitive traits are influenced by the same genes, but to differential levels, from young adulthood to old age. In contrast, increasing white-matter incoherence with age appears to be nongenetic. These results clearly demonstrate that traits sensitive to the genetic influences on brain aging can be identified, a critical first step in delineating the biological mechanisms of successful aging

Genetic Variants Related to Cardiometabolic Traits Are Associated to B Cell Function, Insulin Resistance, and Diabetes Among AmeriCan Indians: The Strong Heart Family Study

Background: Genetic research may inform underlying mechanisms for disparities in the burden of type 2 diabetes mellitus among American Indians. Our objective was to assess the association of genetic variants in cardiometabolic candidate genes with B cell dysfunction via HOMA-B, insulin resistance via HOMA-IR, and type 2 diabetes mellitus in the Strong Heart Family Study (SHFS). Methods and Results: We examined the association of variants, previously associated with cardiometabolic traits (∼200,000 from Illumina Cardio MetaboChip), using mixed models of HOMA-B residuals corrected for HOMA-IR (cHOMA-B), log transformed HOMA-IR, and incident diabetes, adjusted for age, sex, population stratification, and familial relatedness. Center-specific estimates were combined using fixed effect meta-analyses. We used Bonferroni correction to account for multiple testing (P \u3c 4.13 × 10−7). We also assessed the association between variants in candidate diabetes genes with these metabolic traits. We explored the top SNPs in an independent, replication sample from Southwestern Arizona. We identified significant associations with cHOMA-B for common variants at 26 loci of which 8 were novel (PRSS7, FCRL5, PEL1, LRP12, IGLL1, ARHGEF10, PARVA, FLJ16686). The most significant variant association with cHOMA-B was observed on chromosome 5 for an intergenic variant near PARP8 (rs2961831, P = 6.39 × 10−9). In the replication study, we found a signal at rs4607517 near GCK/YKT6 (P = 0.01). Variants near candidate diabetes genes (especially GCK and KCNQ1) were also nominally associated with HOMA-IR and cHOMA-B. Conclusion: We identified variants at novel loci and confirmed those at known candidate diabetes loci associations for cHOMA-B. This study also provided evidence for association of variants at KCNQ2, CTNAA2, and KCNQ1with cHOMA-B among American Indians. Further studies are needed to account for the high heritability of diabetes among the American Indian participants of the SHFS cohort

Independent test assessment using the extreme value distribution theory

The new generation of whole genome sequencing platforms offers great possibilities and challenges for dissecting the genetic basis of complex traits. With a very high number of sequence variants, a naïve multiple hypothesis threshold correction hinders the identification of reliable associations by the overreduction of statistical power. In this report, we examine 2 alternative approaches to improve the statistical power of a whole genome association study to detect reliable genetic associations. The approaches were tested using the Genetic Analysis Workshop 19 (GAW19) whole genome sequencing data. The first tested method estimates the real number of effective independent tests actually being performed in whole genome association project by the use of an extreme value distribution and a set of phenotype simulations. Given the familiar nature of the GAW19 data and the finite number of pedigree founders in the sample, the number of correlations between genotypes is greater than in a set of unrelated samples. Using our procedure, we estimate that the effective number represents only 15 % of the total number of independent tests performed. However, even using this corrected significance threshold, no genome-wide significant association could be detected for systolic and diastolic blood pressure traits. The second approach implements a biological relevance-driven hypothesis tested by exploiting prior computational predictions on the effect of nonsynonymous genetic variants detected in a whole genome sequencing association study. This guided testing approach was able to identify 2 promising single-nucleotide polymorphisms (SNPs), 1 for each trait, targeting biologically relevant genes that could help shed light on the genesis of the human hypertension. The first gene, PFH14, associated with systolic blood pressure, interacts directly with genes involved in calcium-channel formation and the second gene, MAP4, encodes a microtubule-associated protein and had already been detected by previous genome-wide association study experiments conducted in an Asian population. Our results highlight the necessity of the development of alternative approached to improve the efficiency on the detection of reasonable candidate associations in whole genome sequencing studies

Do rare variant genotypes predict common variant genotypes?

The synthetic association hypothesis proposes that common genetic variants detectable in genome-wide association studies may reflect the net phenotypic effect of multiple rare polymorphisms distributed broadly within the focal gene rather than, as often assumed, the effect of common functional variants in high linkage disequilibrium with the focal marker. In a recent study, Dickson and colleagues demonstrated synthetic association in simulations and in two well-characterized, highly polymorphic human disease genes. The converse of this hypothesis is that rare variant genotypes must be correlated with common variant genotypes often enough to make the phenomenon of synthetic association possible. Here we used the exome genotype data provided for Genetic Analysis Workshop 17 to ask how often, how well, and under what conditions rare variant genotypes predict the genotypes of common variants within the same gene. We found nominal evidence of correlation between rare and common variants in 21-30% of cases examined for unrelated individuals; this rate increased to 38-44% for related individuals, underscoring the segregation that underlies synthetic association

Data for Genetic Analysis Workshop 18: human whole genome sequence, blood pressure, and simulated phenotypes in extended pedigrees

Genetic Analysis Workshop 18 (GAW18) focused on identification of genes and functional variants that influence complex phenotypes in human sequence data. Data for the workshop were donated by the T2D-GENES Consortium and included whole genome sequences for odd-numbered autosomes in 464 key individuals selected from 20 Mexican American families, a dense set of single-nucleotide polymorphisms in 959 individuals in these families, and longitudinal data on systolic and diastolic blood pressure measured at 1-4 examinations over a period of 20 years. Simulated phenotypes were generated based on the real sequence data and pedigree structures. In the design of the simulation model, gene expression measures from the San Antonio Family Heart Study (not distributed as part of the GAW18 data) were used to identify genes whose mRNA levels were correlated with blood pressure. Observed variants within these genes were designated as functional in the GAW18 simulation if they were nonsynonymous and predicted to have deleterious effects on protein function or if they were noncoding and associated with mRNA levels. Two simulated longitudinal phenotypes were modeled to have the same trait distributions as the real systolic and diastolic blood pressure data, with effects of age, sex, and medication use, including a genotype-medication interaction. For each phenotype, more than 1000 sequence variants in more than 200 genes present on the odd-numbered autosomes individually explained less than 0.01-2.78% of phenotypic variance. Cumulatively, variants in the most influential gene explained 7.79% of trait variance. An additional simulated phenotype, Q1, was designed to be correlated among family members but to not be associated with any sequence variants. Two hundred replicates of the phenotypes were simulated, with each including data for 849 individuals

Genetic signal maximization using environmental regression

Joint analyses of correlated phenotypes in genetic epidemiology studies are common. However, these analyses primarily focus on genetic correlation between traits and do not take into account environmental correlation. We describe a method that optimizes the genetic signal by accounting for stochastic environmental noise through joint analysis of a discrete trait and a correlated quantitative marker. We conducted bivariate analyses where heritability and the environmental correlation between the discrete and quantitative traits were calculated using Genetic Analysis Workshop 17 (GAW17) family data. The resulting inverse value of the environmental correlation between these traits was then used to determine a new β coefficient for each quantitative trait and was constrained in a univariate model. We conducted genetic association tests on 7,087 nonsynonymous SNPs in three GAW17 family replicates for Affected status with the β coefficient fixed for three quantitative phenotypes and compared these to an association model where the β coefficient was allowed to vary. Bivariate environmental correlations were 0.64 (± 0.09) for Q1, 0.798 (± 0.076) for Q2, and −0.169 (± 0.18) for Q4. Heritability of Affected status improved in each univariate model where a constrained β coefficient was used to account for stochastic environmental effects. No genome-wide significant associations were identified for either method but we demonstrated that constraining β for covariates slightly improved the genetic signal for Affected status. This environmental regression approach allows for increased heritability when the β coefficient for a highly correlated quantitative covariate is constrained and increases the genetic signal for the discrete trait