Estimation of pleiotropy between complex diseases using single-nucleotide polymorphism-derived genomic relationships and restricted maximum likelihood

Application of linear mixed models to study genetic stability of height and body mass index across countries and time

Background: It is now possible to estimate genetic correlations between two independent samples when there is no overlapping phenotypic information. We applied the latest bivariate genomic methods to children in the UK and older adults in Sweden to ask two questions. Are the same variants driving individual differences in anthropometric traits in these two populations, and are these variants as important in childhood as they are later in life?

Methods: A sample of 3152 11-year-old children in the UK was compared with a sample of 6813 adults with an average age of 65 in Sweden. Genotypes were imputed from 1000 genomes with combined 9 767 136 single nucleotide polymorphisms meeting quality control criteria in both samples. Two cross-sample GCTA-GREML analyses and linkage disequilibrium (LD) score regressions were conducted to assess genetic correlations across more than 50 years: child versus adult height and child versus adult body mass index (BMI). Consistency of effects was tested using the recently proposed polygenic scoring method.

Results: For height, GCTA-GREML and LD score indicated strong genetic stability between children and adults, 0.58 (0.16) and 1.335 (1.09), respectively. For BMI, both methods produced similarly strong estimates of genetic stability 0.75 (0.26) and 0.855 (0.49), respectively. In height, adult polygenic score explained 60% of genetic variance in childhood and 10% of variance in BMI.

Conclusions: Here we replicated and extended previous findings of longitudinal genetic stability in anthropometric traits to cross-cultural dimensions, and showed that for height but not BMI these variants are as important in childhood as they are in adulthood.

Genetic determination of height-mediated mate choice

BACKGROUND

Numerous studies have reported positive correlations among couples for height. This suggests that humans find individuals of similar height attractive. However, the answer to whether the choice of a mate with a similar phenotype is genetically or environmentally determined has been elusive.

RESULTS

Here we provide an estimate of the genetic contribution to height choice in mates in 13,068 genotyped couples. Using a mixed linear model we show that 4.1% of the variation in the mate height choice is determined by a person's own genotype, as expected in a model where one's height determines the choice of mate height. Furthermore, the genotype of an individual predicts their partners' height in an independent dataset of 15,437 individuals with 13% accuracy, which is 64% of the theoretical maximum achievable with a heritability of 0.041. Theoretical predictions suggest that approximately 5% of the heritability of height is due to the positive covariance between allelic effects at different loci, which is caused by assortative mating. Hence, the coupling of alleles with similar effects could substantially contribute to the missing heritability of height.

CONCLUSIONS

These estimates provide new insight into the mechanisms that govern mate choice in humans and warrant the search for the genetic causes of choice of mate height. They have important methodological implications and contribute to the missing heritability debate.

Genome-Wide Estimates of Heritability for Social Demographic Outcomes

An increasing number of studies that are widely used in the demographic research community have collected genome-wide data from their respondents. It is therefore important that demographers have a proper understanding of some of the methodological tools needed to analyze such data. This article details the underlying methodology behind one of the most common techniques for analyzing genome-wide data, genome-wide complex trait analysis (GCTA). GCTA models provide heritability estimates for health, health behaviors, or indicators of attainment using data from unrelated persons. Our goal was to describe this model, highlight the utility of the model for biodemographic research, and demonstrate the performance of this approach under modifications to the underlying assumptions. The first set of modifications involved changing the nature of the genetic data used to compute genetic similarities between individuals (the genetic relationship matrix). We then explored the sensitivity of the model to heteroscedastic errors. In general, GCTA estimates are found to be robust to the modifications proposed here, but we also highlight potential limitations of GCTA estimates.

The contribution of rare variation to prostate cancer heritability

We report targeted sequencing of 63 known prostate cancer risk regions in a multi-ancestry study of 9,237 men and use the data to explore the contribution of low-frequency variation to disease risk. We show that SNPs with minor allele frequencies (MAFs) of 0.1-1% explain a substantial fraction of prostate cancer risk in men of African ancestry. We estimate that these SNPs account for 0.12 (standard error (s.e.) = 0.05) of variance in risk (∼42% of the variance contributed by SNPs with MAF of 0.1-50%). This contribution is much larger than the fraction of neutral variation due to SNPs in this class, implying that natural selection has driven down the frequency of many prostate cancer risk alleles; we estimate the coupling between selection and allelic effects at 0.48 (95% confidence interval [0.19, 0.78]) under the Eyre-Walker model. Our results indicate that rare variants make a disproportionate contribution to genetic risk for prostate cancer and suggest the possibility that rare variants may also have an outsize effect on other common traits.

A large genome-wide association study of age-related macular degeneration highlights contributions of rare and common variants

Advanced age-related macular degeneration (AMD) is the leading cause of blindness in the elderly with limited therapeutic options. Here, we report on a study of >12 million variants including 163,714 directly genotyped, most rare, protein-altering variant. Analyzing 16,144 patients and 17,832 controls, we identify 52 independently associated common and rare variants (P < 5×10^–8) distributed across 34 loci. While wet and dry AMD subtypes exhibit predominantly shared genetics, we identify the first signal specific to wet AMD, near MMP9 (difference-P = 4.1×10^–10). Very rare coding variants (frequency < 0.1%) in CFH, CFI, and TIMP3 suggest causal roles for these genes, as does a splice variant in SLC16A8. Our results support the hypothesis that rare coding variants can pinpoint causal genes within known genetic loci and illustrate that applying the approach systematically to detect new loci requires extremely large sample sizes.

Genetic Control over mtDNA and Its Relationship to Major Depressive Disorder

Control over the number of mtDNA molecules per cell appears to be tightly regulated, but the mechanisms involved are largely unknown. Reversible alterations in the amount of mtDNA occur in response to stress suggesting that control over the amount of mtDNA is involved in stress-related diseases including major depressive disorder (MDD). Using low-coverage sequence data from 10,442 Chinese women to compute the normalized numbers of reads mapping to the mitochondrial genome as a proxy for the amount of mtDNA, we identified two loci that contribute to mtDNA levels: one within the TFAM gene on chromosome 10 (rs11006126, p value = 8.73 × 10(-28), variance explained = 1.90%) and one over the CDK6 gene on chromosome 7 (rs445, p value = 6.03 × 10(-16), variance explained = 0.50%). Both loci replicated in an independent cohort. CDK6 is thus a new molecule involved in the control of mtDNA. We identify increased rates of heteroplasmy in women with MDD, and show from an experimental paradigm using mice that the increase is likely due to stress. Furthermore, at least one heteroplasmic variant is significantly associated with changes in the amount of mtDNA (position 513, p value = 3.27 × 10(-9), variance explained = 0.48%) suggesting site-specific heteroplasmy as a possible link between stress and increase in amount of mtDNA. These findings indicate the involvement of mitochondrial genome copy number and sequence in an organism's response to stress.

Quantifying the heritability of glioma using genome-wide complex trait analysis

Genome-wide association studies (GWAS) have successfully identified a number of common single-nucleotide polymorphisms (SNPs) influencing glioma risk. While these SNPs only explain a small proportion of the genetic risk it is unclear how much is left to be detected by other, yet to be identified, common SNPs. Therefore, we applied Genome-Wide Complex Trait Analysis (GCTA) to three GWAS datasets totalling 3,373 cases and 4,571 controls and performed a meta-analysis to estimate the heritability of glioma. Our results identify heritability estimates of 25% (95% CI: 20-31%, P = 1.15 × 10(-17)) for all forms of glioma - 26% (95% CI: 17-35%, P = 1.05 × 10(-8)) for glioblastoma multiforme (GBM) and 25% (95% CI: 17-32%, P = 1.26 × 10(-10)) for non-GBM tumors. This is a substantial increase from the genetic variance identified by the currently identified GWAS risk loci (~6% of common heritability), indicating that most of the heritable risk attributable to common genetic variants remains to be identified.

Contrasting genetic architectures of schizophrenia and other complex diseases using fast variance-components analysis

Heritability analyses of genome-wide association study (GWAS) cohorts have yielded important insights into complex disease architecture, and increasing sample sizes hold the promise of further discoveries. Here we analyze the genetic architectures of schizophrenia in 49,806 samples from the PGC and nine complex diseases in 54,734 samples from the GERA cohort. For schizophrenia, we infer an overwhelmingly polygenic disease architecture in which ≥71% of 1-Mb genomic regions harbor ≥1 variant influencing schizophrenia risk. We also observe significant enrichment of heritability in GC-rich regions and in higher-frequency SNPs for both schizophrenia and GERA diseases. In bivariate analyses, we observe significant genetic correlations (ranging from 0.18 to 0.85) for several pairs of GERA diseases; genetic correlations were on average 1.3 tunes stronger than the correlations of overall disease liabilities. To accomplish these analyses, we developed a fast algorithm for multicomponent, multi-trait variance-components analysis that overcomes prior computational barriers that made such analyses intractable at this scale.

Uncovering the Genetic Architectures of Quantitative Traits

The aim of a genome-wide association study (GWAS) is to identify loci in the human genome affecting a phenotype of interest. This review summarizes some recent work on conceptual and methodological aspects of GWAS. The average effect of gene substitution at a given causal site in the genome is the key estimand in GWAS, and we argue for its fundamental importance. Implicit in the definition of average effect is a linear model relating genotype to phenotype. The fraction of the phenotypic variance ascribable to polymorphic sites with nonzero average effects in this linear model is called the heritability, and we describe methods for estimating this quantity from GWAS data. Finally, we show that the theory of compressed sensing can be used to provide a sharp estimate of the sample size required to identify essentially all sites contributing to the heritability of a given phenotype.

PAPA: a flexible tool for identifying pleiotropic pathways using genome-wide association study summaries

: Pleiotropy is common in the genetic architectures of complex diseases. To the best of our knowledge, no analysis tool has been developed for identifying pleiotropic pathways using multiple genome-wide association study (GWAS) summaries by now. Here, we present PAPA, a flexible tool for pleiotropic pathway analysis utilizing GWAS summary results. The performance of PAPA was validated using publicly available GWAS summaries of body mass index and waist-hip ratio of the GIANT datasets. PAPA identified a set of pleiotropic pathways, which have been demonstrated to be involved in the development of obesity.

AVAILABILITY AND IMPLEMENTATION

PAPA program, document and illustrative example are available at http://sourceforge.net/projects/papav1/files/

CONTACT

: fzhxjtu@mail.xjtu.edu.cn

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genome-wide analysis of genetic correlation in dementia with Lewy bodies, Parkinson's and Alzheimer's diseases

The similarities between dementia with Lewy bodies (DLB) and both Parkinson's disease (PD) and Alzheimer's disease (AD) are many and range from clinical presentation, to neuropathological characteristics, to more recently identified, genetic determinants of risk. Because of these overlapping features, diagnosing DLB is challenging and has clinical implications since some therapeutic agents that are applicable in other diseases have adverse effects in DLB. Having shown that DLB shares some genetic risk with PD and AD, we have now quantified the amount of sharing through the application of genetic correlation estimates, and show that, from a purely genetic perspective, and excluding the strong association at the APOE locus, DLB is equally correlated to AD and PD.

An atlas of genetic correlations across human diseases and traits

Identifying genetic correlations between complex traits and diseases can provide useful etiological insights and help prioritize likely causal relationships. The major challenges preventing estimation of genetic correlation from genome-wide association study (GWAS) data with current methods are the lack of availability of individual-level genotype data and widespread sample overlap among meta-analyses. We circumvent these difficulties by introducing a technique-cross-trait LD Score regression-for estimating genetic correlation that requires only GWAS summary statistics and is not biased by sample overlap. We use this method to estimate 276 genetic correlations among 24 traits. The results include genetic correlations between anorexia nervosa and schizophrenia, anorexia and obesity, and educational attainment and several diseases. These results highlight the power of genome-wide analyses, as there currently are no significantly associated SNPs for anorexia nervosa and only three for educational attainment.

Identification of Promising Mutants Associated with Egg Production Traits Revealed by Genome-Wide Association Study

Egg number (EN), egg laying rate (LR) and age at first egg (AFE) are important production traits related to egg production in poultry industry. To better understand the knowledge of genetic architecture of dynamic EN during the whole laying cycle and provide the precise positions of associated variants for EN, LR and AFE, laying records from 21 to 72 weeks of age were collected individually for 1,534 F2 hens produced by reciprocal crosses between White Leghorn and Dongxiang Blue-shelled chicken, and their genotypes were assayed by chicken 600 K Affymetrix high density genotyping arrays. Subsequently, pedigree and SNP-based genetic parameters were estimated and a genome-wide association study (GWAS) was conducted on EN, LR and AFE. The heritability estimates were similar between pedigree and SNP-based estimates varying from 0.17 to 0.36. In the GWA analysis, we identified nine genome-wide significant loci associated with EN of the laying periods from 21 to 26 weeks, 27 to 36 weeks and 37 to 72 weeks. Analysis of GTF2A1 and CLSPN suggested that they influenced the function of ovary and uterus, and may be considered as relevant candidates. The identified SNP rs314448799 for accumulative EN from 21 to 40 weeks on chromosome 5 created phenotypic differences of 6.86 eggs between two homozygous genotypes, which could be potentially applied to the molecular breeding for EN selection. Moreover, our finding showed that LR was a moderate polygenic trait. The suggestive significant region on chromosome 16 for AFE suggested the relationship between sex maturity and immune in the current population. The present study comprehensively evaluates the role of genetic variants in the development of egg laying. The findings will be helpful to investigation of causative genes function and future marker-assisted selection and genomic selection in chickens.

Genome-wide genetic homogeneity between sexes and populations for human height and body mass index

Sex-specific genetic effects have been proposed to be an important source of variation for human complex traits. Here we use two distinct genome-wide methods to estimate the autosomal genetic correlation (rg) between men and women for human height and body mass index (BMI), using individual-level (n = ∼44 000) and summary-level (n = ∼133 000) data from genome-wide association studies. Results are consistent and show that the between-sex genetic correlation is not significantly different from unity for both traits. In contrast, we find evidence of genetic heterogeneity between sexes for waist-hip ratio (rg = ∼0.7) and between populations for BMI (rg = ∼0.9 between Europe and the USA) but not for height. The lack of evidence for substantial genetic heterogeneity for body size is consistent with empirical findings across traits and species.

Genome-wide association studies for feed intake and efficiency in two laying periods of chickens

Background

Feed contributes to over 60 % of the total production costs in the poultry industry. Increasing feed costs prompt geneticists to include feed intake and efficiency as selection goals in breeding programs. In the present study, we used an F₂ chicken population in a genome-wide association study (GWAS) to detect potential genetic variants and candidate genes associated with daily feed intake (FI) and feed efficiency, including residual feed intake (RFI) and feed conversion ratio (FCR).

Methods

A total of 1534 F₂ hens from a White Leghorn and Dongxiang reciprocal cross were phenotyped for feed intake and efficiency between 37 and 40 weeks (FI1, RFI1, and FCR1) and between 57 and 60 weeks (FI2, RFI2, and FCR2), and genotyped using the chicken 600 K single nucleotide polymorphism (SNP) genotyping array. Univariate, bivariate, and conditional genome-wide association studies (GWAS) were performed with GEMMA, a genome-wide efficient mixed model association algorithm. The statistical significance threshold for association was inferred by the simpleM method.

Results

We identified eight genomic regions that each contained at least one genetic variant that showed a significant association with FI. Genomic regions on Gallus gallus (GGA) chromosome 4 coincided with known quantitative trait loci (QTL) that affect feed intake of layers. Of particular interest, eight SNPs on GGA1 in the region between 169.23 and 171.55 Mb were consistently associated with FI in both univariate and bivariate GWAS, which explained 3.72 and 2.57 % of the phenotypic variance of FI1 and FI2, respectively. The CAB39L gene can be considered as a promising candidate for FI1. For RFI, a haplotype block on GGA27 harbored a significant SNP associated with RFI2. The major allele of rs315135692 was favorable for a lower RFI, with a phenotypic difference of 3.35 g/day between opposite homozygous genotypes. Strong signals on GGA1 were detected in the bivariate GWAS for FCR.

Conclusions

The results demonstrated the polygenic nature of feed intake. GWAS identified novel variants and confirmed a QTL that was previously reported for feed intake in chickens. Genetic variants associated with feed efficiency may be used in genomic breeding programs to select more efficient layers.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0161-1) contains supplementary material, which is available to authorized users.

Shared genetics underlying epidemiological association between endometriosis and ovarian cancer

Epidemiological studies have demonstrated associations between endometriosis and certain histotypes of ovarian cancer, including clear cell, low-grade serous and endometrioid carcinomas. We aimed to determine whether the observed associations might be due to shared genetic aetiology. To address this, we used two endometriosis datasets genotyped on common arrays with full-genome coverage (3194 cases and 7060 controls) and a large ovarian cancer dataset genotyped on the customized Illumina Infinium iSelect (iCOGS) arrays (10 065 cases and 21 663 controls). Previous work has suggested that a large number of genetic variants contribute to endometriosis and ovarian cancer (all histotypes combined) susceptibility. Here, using the iCOGS data, we confirmed polygenic architecture for most histotypes of ovarian cancer. This led us to evaluate if the polygenic effects are shared across diseases. We found evidence for shared genetic risks between endometriosis and all histotypes of ovarian cancer, except for the intestinal mucinous type. Clear cell carcinoma showed the strongest genetic correlation with endometriosis (0.51, 95% CI = 0.18-0.84). Endometrioid and low-grade serous carcinomas had similar correlation coefficients (0.48, 95% CI = 0.07-0.89 and 0.40, 95% CI = 0.05-0.75, respectively). High-grade serous carcinoma, which often arises from the fallopian tubes, showed a weaker genetic correlation with endometriosis (0.25, 95% CI = 0.11-0.39), despite the absence of a known epidemiological association. These results suggest that the epidemiological association between endometriosis and ovarian adenocarcinoma may be attributable to shared genetic susceptibility loci.

Genetic sharing and heritability of paediatric age of onset autoimmune diseases

Autoimmune diseases (AIDs) are polygenic diseases affecting 7-10% of the population in the Western Hemisphere with few effective therapies. Here, we quantify the heritability of paediatric AIDs (pAIDs), including JIA, SLE, CEL, T1D, UC, CD, PS, SPA and CVID, attributable to common genomic variations (SNP-h(2)). SNP-h(2) estimates are most significant for T1D (0.863±s.e. 0.07) and JIA (0.727±s.e. 0.037), more modest for UC (0.386±s.e. 0.04) and CD (0.454±0.025), largely consistent with population estimates and are generally greater than that previously reported by adult GWAS. On pairwise analysis, we observed that the diseases UC-CD (0.69±s.e. 0.07) and JIA-CVID (0.343±s.e. 0.13) are the most strongly correlated. Variations across the MHC strongly contribute to SNP-h(2) in T1D and JIA, but does not significantly contribute to the pairwise rG. Together, our results partition contributions of shared versus disease-specific genomic variations to pAID heritability, identifying pAIDs with unexpected risk sharing, while recapitulating known associations between autoimmune diseases previously reported in adult cohorts.

Differences in Common Genetic Predisposition to Ischemic Stroke by Age and Sex

Supplemental Digital Content is available in the text.

Evidence from epidemiological studies points to differences in factors predisposing to stroke by age and sex. Whether these arise because of different genetic influences remained untested. Here, we use data from 4 genome-wide association data sets to study the relationship between genetic influence on stroke with both age and sex.

Genome-wide association study dissects genetic architecture underlying longitudinal egg weights in chickens

BACKGROUND

As a major economic trait in chickens, egg weight (EW) receives widespread interests in breeding, production and consumption. However, limited information is available for underlying genetic architecture of longitudinal trend in EW. Herein, we measured EWs at nine time points from onset of laying to 60 week of age, and conducted comprehensive genome-wide association studies (GWAS) in 1,534 F2 hens derived from reciprocal crosses between White Leghorn and Dongxiang chickens.

RESULTS

Egg weights at all ages except the first egg weight (FEW) exhibited high SNP-based heritability estimates (0.47~0.60). Strong pair-wise genetic correlations (0.77~1.00) were found among all EWs. Nine separate univariate genome-wide screens suggested 73 signals showing significant associations with longitudinal EWs. After multivariate and conditional analyses, four variants on three chromosomes remained independent contributions. The minor alleles at two loci exerted consistent and positive substitution effects on EWs, and other two were negative. The four loci together accounted for 3.84 % of the phenotypic variance for FEW and 7.29~11.06 % for EWs from 32 to 60 week of age. We obtained five candidate genes, of which NCAPG harbors a non-synonymous SNP (rs14491030) causing a valine-to-alanine amino-acid substitution. Genome partitioning analysis indicated a strong linear correlation between the variance explained by each chromosome and its length, which provided evidence that EW follows a highly polygenic nature of inheritance.

CONCLUSIONS

Identification of significant genetic causes that together implicate EWs at different ages will greatly advance our understanding of the genetic basis behind longitudinal EWs, and would be helpful to illuminate the future breeding direction on how to select desired egg size.

Genetic variance estimation with imputed variants finds negligible missing heritability for human height and body mass index

We propose a method (GREML-LDMS) to estimate heritability for human complex traits in unrelated individuals using whole-genome sequencing data. We demonstrate using simulations based on whole-genome sequencing data that ∼97% and ∼68% of variation at common and rare variants, respectively, can be captured by imputation. Using the GREML-LDMS method, we estimate from 44,126 unrelated individuals that all ∼17 million imputed variants explain 56% (standard error (s.e.) = 2.3%) of variance for height and 27% (s.e. = 2.5%) of variance for body mass index (BMI), and we find evidence that height- and BMI-associated variants have been under natural selection. Considering the imperfect tagging of imputation and potential overestimation of heritability from previous family-based studies, heritability is likely to be 60-70% for height and 30-40% for BMI. Therefore, the missing heritability is small for both traits. For further discovery of genes associated with complex traits, a study design with SNP arrays followed by imputation is more cost-effective than whole-genome sequencing at current prices.

Epigenomic mapping and effect sizes of noncoding variants associated with psychotropic drug response

AIM

To provide insight into potential regulatory mechanisms of gene expression underlying addiction, analgesia, psychotropic drug response and adverse drug events, genome-wide association studies searching for variants associated with these phenotypes has been undertaken with limited success. We undertook analysis of these results with the aim of applying epigenetic knowledge to aid variant discovery and interpretation.

METHODS

We applied conditional imputation to results from 26 genome-wide association studies and three candidate gene-association studies. The analysis workflow included data from chromatin conformation capture, chromatin state annotation, DNase I hypersensitivity, hypomethylation, anatomical localization and biochronicity. We also made use of chromatin state data from the epigenome roadmap, transcription factor-binding data, spatial maps from published Hi-C datasets and 'guilt by association' methods.

RESULTS

We identified 31 pharmacoepigenomic SNPs from a total of 2024 variants in linkage disequilibrium with lead SNPs, of which only 6% were coding variants. Interrogation of chromatin state using our workflow and the epigenome roadmap showed agreement on 34 of 35 tissue assignments to regulatory elements including enhancers and promoters. Loop boundary domains were inferred by association with CTCF (CCCTC-binding factor) and cohesin, suggesting proximity to topologically associating domain boundaries and enhancer clusters. Spatial interactions between enhancer-promoter pairs detected both known and previously unknown mechanisms. Addiction and analgesia SNPs were common in relevant populations and exhibited large effect sizes, whereas a SNP located in the promoter of the SLC1A2 gene exhibited a moderate effect size for lithium response in bipolar disorder in patients of European ancestry. SNPs associated with drug-induced organ injury were rare but exhibited the largest effect sizes, consistent with the published literature.

CONCLUSION

This work demonstrates that an in silico bioinformatics-based approach using integrative analysis of a diversity of molecular and morphological data types can discover pharmacoepigenomic variants that are suitable candidates for further validation in cell lines, animal models and human clinical trials.

Association analyses identify 38 susceptibility loci for inflammatory bowel disease and highlight shared genetic risk across populations

Ulcerative colitis and Crohn's disease are the two main forms of inflammatory bowel disease (IBD). Here we report the first trans-ancestry association study of IBD, with genome-wide or Immunochip genotype data from an extended cohort of 86,640 European individuals and Immunochip data from 9,846 individuals of East Asian, Indian or Iranian descent. We implicate 38 loci in IBD risk for the first time. For the majority of the IBD risk loci, the direction and magnitude of effect are consistent in European and non-European cohorts. Nevertheless, we observe genetic heterogeneity between divergent populations at several established risk loci driven by differences in allele frequency (NOD2) or effect size (TNFSF15 and ATG16L1) or a combination of these factors (IL23R and IRGM). Our results provide biological insights into the pathogenesis of IBD and demonstrate the usefulness of trans-ancestry association studies for mapping loci associated with complex diseases and understanding genetic architecture across diverse populations.

Do Molecular Markers Inform About Pleiotropy?

The availability of dense panels of common single-nucleotide polymorphisms and sequence variants has facilitated the study of statistical features of the genetic architecture of complex traits and diseases via whole-genome regressions (WGRs). At the onset, traits were analyzed trait by trait, but recently, WGRs have been extended for analysis of several traits jointly. The expectation is that such an approach would offer insight into mechanisms that cause trait associations, such as pleiotropy. We demonstrate that correlation parameters inferred using markers can give a distorted picture of the genetic correlation between traits. In the absence of knowledge of linkage disequilibrium relationships between quantitative or disease trait loci and markers, speculating about genetic correlation and its causes (e.g., pleiotropy) using genomic data is conjectural.

Familial aggregation and childhood blood pressure

There is growing concern about elevated blood pressure (BP) in children. The evidence for familial aggregation of childhood BP is substantial. Twin studies have shown that a large part of the familial aggregation of childhood BP is due to genes. The first part of this review provides the latest progress in gene finding for childhood BP, focusing on the combined effects of multiple loci identified from the genome-wide association studies on adult BP. We further review the evidence on the contribution of the genetic components of other family risk factors to the familial aggregation of childhood BP including obesity, birth weight, sleep quality, sodium intake, parental smoking, and socioeconomic status. At the end, we emphasize the promise of using genomic-relatedness-matrix restricted maximum likelihood (GREML) analysis, a method that uses genome-wide data from unrelated individuals, in answering a number of unsolved questions in the familial aggregation of childhood BP.

Markers of Psychological Differences and Social and Health Inequalities: Possible Genetic and Phenotypic Overlaps

Associations between markers of ostensible psychological characteristics and social and health inequalities are pervasive but difficult to explain. In some cases, there may be causal influence flowing from social and health inequalities to psychological differences, whereas sometimes it may be the other way around. Here, we focus on the possibility that some markers that we often consider as indexing different domains of individual differences may in fact reflect at least partially overlapping genetic and/or phenotypic bases. For example, individual differences in cognitive abilities and educational attainment appear to reflect largely overlapping genetic influences, whereas cognitive abilities and health literacy may be almost identical phenomena at the phenotypic, never mind genetic, level. We make the case for employing molecular genetic data and quantitative genetic techniques to better understand the associations of psychological individual differences with social and health inequalities. We illustrate these arguments by using published findings from the Lothian Birth Cohort and the Generation Scotland studies. We also present novel findings pertaining to longitudinal stability and change in older age personality traits and some correlates of the change, molecular genetic data-based heritability estimates of Neuroticism and Extraversion, and the genetic correlations of these personality traits with markers of social and health inequalities.

A Fast Method that Uses Polygenic Scores to Estimate the Variance Explained by Genome-wide Marker Panels and the Proportion of Variants Affecting a Trait

Several methods have been proposed to estimate the variance in disease liability explained by large sets of genetic markers. However, current methods do not scale up well to large sample sizes. Linear mixed models require solving high-dimensional matrix equations, and methods that use polygenic scores are very computationally intensive. Here we propose a fast analytic method that uses polygenic scores, based on the formula for the non-centrality parameter of the association test of the score. We estimate model parameters from the results of multiple polygenic score tests based on markers with p values in different intervals. We estimate parameters by maximum likelihood and use profile likelihood to compute confidence intervals. We compare various options for constructing polygenic scores, based on nested or disjoint intervals of p values, weighted or unweighted effect sizes, and different numbers of intervals, in estimating the variance explained by a set of markers, the proportion of markers with effects, and the genetic covariance between a pair of traits. Our method provides nearly unbiased estimates and confidence intervals with good coverage, although estimation of the variance is less reliable when jointly estimated with the covariance. We find that disjoint p value intervals perform better than nested intervals, but the weighting did not affect our results. A particular advantage of our method is that it can be applied to summary statistics from single markers, and so can be quickly applied to large consortium datasets. Our method, named AVENGEME (Additive Variance Explained and Number of Genetic Effects Method of Estimation), is implemented in R software.

Evaluation of Penalized and Nonpenalized Methods for Disease Prediction with Large-Scale Genetic Data

Owing to recent improvement of genotyping technology, large-scale genetic data can be utilized to identify disease susceptibility loci and this successful finding has substantially improved our understanding of complex diseases. However, in spite of these successes, most of the genetic effects for many complex diseases were found to be very small, which have been a big hurdle to build disease prediction model. Recently, many statistical methods based on penalized regressions have been proposed to tackle the so-called “large P and small N” problem. Penalized regressions including least absolute selection and shrinkage operator (LASSO) and ridge regression limit the space of parameters, and this constraint enables the estimation of effects for very large number of SNPs. Various extensions have been suggested, and, in this report, we compare their accuracy by applying them to several complex diseases. Our results show that penalized regressions are usually robust and provide better accuracy than the existing methods for at least diseases under consideration.

Lifetime Socioeconomic Status, Historical Context, and Genetic Inheritance in Shaping Body Mass in Middle and Late Adulthood

This study demonstrates body mass in middle and late adulthood as a consequence of the complex interplay among individuals' genes, lifetime socioeconomic experiences, and the historical context in which they live. Drawing on approximately 9,000 genetic samples from the Health and Retirement Study, we first investigate how socioeconomic status (SES) over the life course moderates the impact of 32 established obesity-related genetic variants on body mass index (BMI) in middle and late adulthood. Further, we consider differences across birth cohorts in the genetic influence on BMI and cohort variations in the moderating effects of life-course SES on the genetic influence. Our analyses suggest that persistently low SES over the life course or downward mobility (e.g., high SES in childhood but low SES in adulthood) amplified the genetic influence on BMI, while persistently high SES or upward mobility (e.g., low SES in childhood but high SES in adulthood) compensated for such influence. For more recent birth cohorts, while the genetic influence on BMI became stronger, the moderating effects of lifetime SES on the genetic influence were weaker compared to earlier cohorts. We discuss these findings in light of social changes during the obesity epidemic in the United States.

Genomic architecture of human neuroanatomical diversity

Human brain anatomy is strikingly diverse and highly inheritable: genetic factors may explain up to 80% of its variability. Prior studies have tried to detect genetic variants with a large effect on neuroanatomical diversity, but those currently identified account for <5% of the variance. Here, based on our analyses of neuroimaging and whole-genome genotyping data from 1765 subjects, we show that up to 54% of this heritability is captured by large numbers of single-nucleotide polymorphisms of small-effect spread throughout the genome, especially within genes and close regulatory regions. The genetic bases of neuroanatomical diversity appear to be relatively independent of those of body size (height), but shared with those of verbal intelligence scores. The study of this genomic architecture should help us better understand brain evolution and disease.

Pleiotropy across academic subjects at the end of compulsory education

Research has shown that genes play an important role in educational achievement. A key question is the extent to which the same genes affect different academic subjects before and after controlling for general intelligence. The present study investigated genetic and environmental influences on, and links between, the various subjects of the age-16 UK-wide standardized GCSE (General Certificate of Secondary Education) examination results for 12,632 twins. Using the twin method that compares identical and non-identical twins, we found that all GCSE subjects were substantially heritable, and that various academic subjects correlated substantially both phenotypically and genetically, even after controlling for intelligence. Further evidence for pleiotropy in academic achievement was found using a method based directly on DNA from unrelated individuals. We conclude that performance differences for all subjects are highly heritable at the end of compulsory education and that many of the same genes affect different subjects independent of intelligence.

Iron and hepcidin as risk factors in atherosclerosis: what do the genes say?

Background

Previous reports suggested a role for iron and hepcidin in atherosclerosis. Here, we evaluated the causality of these associations from a genetic perspective via (i) a Mendelian randomization (MR) approach, (ii) study of association of atherosclerosis-related single nucleotide polymorphisms (SNPs) with iron and hepcidin, and (iii) estimation of genomic correlations between hepcidin, iron and atherosclerosis.

Results

Analyses were performed in a general population sample. Iron parameters (serum iron, serum ferritin, total iron-binding capacity and transferrin saturation), serum hepcidin and genome-wide SNP data were available for N = 1,819; non-invasive measurements of atherosclerosis (NIMA), i.e., presence of plaque, intima media thickness and ankle-brachial index (ABI), for N = 549. For the MR, we used 12 iron-related SNPs that were previously identified in a genome-wide association meta-analysis on iron status, and assessed associations of individual SNPs and quartiles of a multi-SNP score with NIMA. Quartile 4 versus quartile 1 of the multi-SNP score showed directionally consistent associations with the hypothesized direction of effect for all NIMA in women, indicating that increased body iron status is a risk factor for atherosclerosis in women. We observed no single SNP associations that fit the hypothesized directions of effect between iron and NIMA, except for rs651007, associated with decreased ferritin concentration and decreased atherosclerosis risk. Two of six NIMA-related SNPs showed association with the ratio hepcidin/ferritin, suggesting that an increased hepcidin/ferritin ratio increases atherosclerosis risk. Genomic correlations were close to zero, except for hepcidin and ferritin with ABI at rest [−0.27 (SE 0.34) and −0.22 (SE 0.35), respectively] and ABI after exercise [−0.29 (SE 0.34) and −0.30 (0.35), respectively]. The negative sign indicates an increased atherosclerosis risk with increased hepcidin and ferritin concentrations.

Conclusions

Our results suggest a potential causal role for hepcidin and ferritin in atherosclerosis, and may indicate that iron status is causally related to atherosclerosis in women.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-015-0246-4) contains supplementary material, which is available to authorized users.

Genetic Factors Influencing Coagulation Factor XIII B-Subunit Contribute to Risk of Ischemic Stroke

Background and Purpose—

Abnormal coagulation has been implicated in the pathogenesis of ischemic stroke, but how this association is mediated and whether it differs between ischemic stroke subtypes is unknown. We determined the shared genetic risk between 14 coagulation factors and ischemic stroke and its subtypes.

Methods—

Using genome-wide association study results for 14 coagulation factors from the population-based TwinsUK sample (N≈2000 for each factor), meta-analysis results from the METASTROKE consortium ischemic stroke genome-wide association study (12 389 cases, 62 004 controls), and genotype data for 9520 individuals from the WTCCC2 ischemic stroke study (3548 cases, 5972 controls—the largest METASTROKE subsample), we explored shared genetic risk for coagulation and stroke. We performed three analyses: (1) a test for excess concordance (or discordance) in single nucleotide polymorphism effect direction across coagulation and stroke, (2) an estimation of the joint effect of multiple coagulation-associated single nucleotide polymorphisms in stroke, and (3) an evaluation of common genetic risk between coagulation and stroke.

Results—

One coagulation factor, factor XIII subunit B (FXIIIB), showed consistent effects in the concordance analysis, the estimation of polygenic risk, and the validation with genotype data, with associations specific to the cardioembolic stroke subtype. Effect directions for FXIIIB-associated single nucleotide polymorphisms were significantly discordant with cardioembolic disease (smallest P=5.7×10⁻⁰⁴); the joint effect of FXIIIB-associated single nucleotide polymorphisms was significantly predictive of ischemic stroke (smallest P=1.8×10⁻⁰⁴) and the cardioembolic subtype (smallest P=1.7×10⁻⁰⁴). We found substantial negative genetic covariation between FXIIIB and ischemic stroke (rG=−0.71, P=0.01) and the cardioembolic subtype (rG=−0.80, P=0.03).

Conclusions—

Genetic markers associated with low FXIIIB levels increase risk of ischemic stroke cardioembolic subtype.

Genes associated with RSV lower respiratory tract infection and asthma: the application of genetic epidemiological methods to understand causality

Infants with respiratory syncytial virus (RSV) lower respiratory tract infections (LRIs) are at increased risk for childhood asthma. The objectives of this article are to review the genes associated with both RSV LRI and asthma, review analytic approaches to assessing shared genetic risk and propose a future perspective on how these approaches can help us to understand the role of infant RSV infection as both an important risk factor for asthma and marker of shared genetic etiology between the two conditions. The review of shared genes and thus pathways associated with severity of response to RSV infection and asthma risk can help us to understand mechanisms of disease and ultimately propose new and novel targets for primary prevention of both diseases.

Familiality and SNP heritability of age at onset and episodicity in major depressive disorder

BACKGROUND

Strategies to dissect phenotypic and genetic heterogeneity of major depressive disorder (MDD) have mainly relied on subphenotypes, such as age at onset (AAO) and recurrence/episodicity. Yet, evidence on whether these subphenotypes are familial or heritable is scarce. The aims of this study are to investigate the familiality of AAO and episode frequency in MDD and to assess the proportion of their variance explained by common single nucleotide polymorphisms (SNP heritability).

METHOD

For investigating familiality, we used 691 families with 2-5 full siblings with recurrent MDD from the DeNt study. We fitted (square root) AAO and episode count in a linear and a negative binomial mixed model, respectively, with family as random effect and adjusting for sex, age and center. The strength of familiality was assessed with intraclass correlation coefficients (ICC). For estimating SNP heritabilities, we used 3468 unrelated MDD cases from the RADIANT and GSK Munich studies. After similarly adjusting for covariates, derived residuals were used with the GREML method in GCTA (genome-wide complex trait analysis) software.

RESULTS

Significant familial clustering was found for both AAO (ICC = 0.28) and episodicity (ICC = 0.07). We calculated from respective ICC estimates the maximal additive heritability of AAO (0.56) and episodicity (0.15). SNP heritability of AAO was 0.17 (p = 0.04); analysis was underpowered for calculating SNP heritability of episodicity.

CONCLUSIONS

AAO and episodicity aggregate in families to a moderate and small degree, respectively. AAO is under stronger additive genetic control than episodicity. Larger samples are needed to calculate the SNP heritability of episodicity. The described statistical framework could be useful in future analyses.

Implications of pleiotropy: challenges and opportunities for mining Big Data in biomedicine

Pleiotropy arises when a locus influences multiple traits. Rich GWAS findings of various traits in the past decade reveal many examples of this phenomenon, suggesting the wide existence of pleiotropic effects. What underlies this phenomenon is the biological connection among seemingly unrelated traits/diseases. Characterizing the molecular mechanisms of pleiotropy not only helps to explain the relationship between diseases, but may also contribute to novel insights concerning the pathological mechanism of each specific disease, leading to better disease prevention, diagnosis and treatment. However, most pleiotropic effects remain elusive because their functional roles have not been systematically examined. A systematic investigation requires availability of qualified measurements at multilayered biological processes (e.g., transcription and translation). The rise of Big Data in biomedicine, such as high-quality multi-omics data, biomedical imaging data and electronic medical records of patients, offers us an unprecedented opportunity to investigate pleiotropy. There will be a great need of computationally efficient and statistically rigorous methods for integrative analysis of these Big Data in biomedicine. In this review, we outline many opportunities and challenges in methodology developments for systematic analysis of pleiotropy, and highlight its implications on disease prevention, diagnosis and treatment.

Integrating Multiple Correlated Phenotypes for Genetic Association Analysis by Maximizing Heritability

Many correlated disease variables are analyzed jointly in genetic studies in the hope of increasing power to detect causal genetic variants. One approach involves assessing the relationship between each phenotype and each SNP individually and using a Bonferroni correction for the effective number of tests conducted. Alternatively, one can apply a multivariate regression or a dimension reduction technique, such as principal component analysis, and test for the association with the principal components of the phenotypes rather than the individual phenotypes. Inspired by the previous approaches of combining phenotypes to maximize heritability at individual SNPs, in this paper, we propose to construct a maximally heritable (MaxH) phenotype by taking advantage of the estimated total heritability and co-heritability. The heritability and co-heritability only need to be estimated once; therefore, our method is applicable to genome-wide scans. The MaxH phenotype is a linear combination of the individual phenotypes with increased heritability and power over the phenotypes being combined. Simulations show that the heritability and power achieved agree well with the theory for large samples and two phenotypes. We compare our approach with commonly used methods and assess both the heritability and the power of the MaxH phenotype. Moreover, we provide suggestions for how to choose the phenotypes for combination. An application of our approach to a GWAS on chronic obstructive pulmonary disease shows its practical relevance.

Estimation of Genetic Relationships Between Individuals Across Cohorts and Platforms: Application to Childhood Height

Combining genotype data across cohorts increases power to estimate the heritability due to common single nucleotide polymorphisms (SNPs), based on analyzing a Genetic Relationship Matrix (GRM). However, the combination of SNP data across multiple cohorts may lead to stratification, when for example, different genotyping platforms are used. In the current study, we address issues of combining SNP data from different cohorts, the Netherlands Twin Register (NTR) and the Generation R (GENR) study. Both cohorts include children of Northern European Dutch background (N = 3102 + 2826, respectively) who were genotyped on different platforms. We explore imputation and phasing as a tool and compare three GRM-building strategies, when data from two cohorts are (1) just combined, (2) pre-combined and cross-platform imputed and (3) cross-platform imputed and post-combined. We test these three strategies with data on childhood height for unrelated individuals (N = 3124, average age 6.7 years) to explore their effect on SNP-heritability estimates and compare results to those obtained from the independent studies. All combination strategies result in SNP-heritability estimates with a standard error smaller than those of the independent studies. We did not observe significant difference in estimates of SNP-heritability based on various cross-platform imputed GRMs. SNP-heritability of childhood height was on average estimated as 0.50 (SE = 0.10). Introducing cohort as a covariate resulted in ≈2 % drop. Principal components (PCs) adjustment resulted in SNP-heritability estimates of about 0.39 (SE = 0.11). Strikingly, we did not find significant difference between cross-platform imputed and combined GRMs. All estimates were significant regardless the use of PCs adjustment. Based on these analyses we conclude that imputation with a reference set helps to increase power to estimate SNP-heritability by combining cohorts of the same ethnicity genotyped on different platforms. However, important factors should be taken into account such as remaining cohort stratification after imputation and/or phenotypic heterogeneity between and within cohorts. Whether one should use imputation, or just combine the genotype data, depends on the number of overlapping SNPs in relation to the total number of genotyped SNPs for both cohorts, and their ability to tag all the genetic variance related to the specific trait of interest.

Human fertility, molecular genetics, and natural selection in modern societies

Research on genetic influences on human fertility outcomes such as number of children ever born (NEB) or the age at first childbirth (AFB) has been solely based on twin and family-designs that suffer from problematic assumptions and practical limitations. The current study exploits recent advances in the field of molecular genetics by applying the genomic-relationship-matrix based restricted maximum likelihood (GREML) methods to quantify for the first time the extent to which common genetic variants influence the NEB and the AFB of women. Using data from the UK and the Netherlands (N = 6,758), results show significant additive genetic effects on both traits explaining 10% (SE = 5) of the variance in the NEB and 15% (SE = 4) in the AFB. We further find a significant negative genetic correlation between AFB and NEB in the pooled sample of –0.62 (SE = 0.27, p-value = 0.02). This finding implies that individuals with genetic predispositions for an earlier AFB had a reproductive advantage and that natural selection operated not only in historical, but also in contemporary populations. The observed postponement in the AFB across the past century in Europe contrasts with these findings, suggesting an evolutionary override by environmental effects and underscoring that evolutionary predictions in modern human societies are not straight forward. It emphasizes the necessity for an integrative research design from the fields of genetics and social sciences in order to understand and predict fertility outcomes. Finally, our results suggest that we may be able to find genetic variants associated with human fertility when conducting GWAS-meta analyses with sufficient sample size.

The association between lower educational attainment and depression owing to shared genetic effects? Results in ~25,000 subjects

An association between lower educational attainment (EA) and an increased risk for depression has been confirmed in various western countries. This study examines whether pleiotropic genetic effects contribute to this association. Therefore, data were analyzed from a total of 9662 major depressive disorder (MDD) cases and 14,949 controls (with no lifetime MDD diagnosis) from the Psychiatric Genomics Consortium with additional Dutch and Estonian data. The association of EA and MDD was assessed with logistic regression in 15,138 individuals indicating a significantly negative association in our sample with an odds ratio for MDD 0.78 (0.75-0.82) per standard deviation increase in EA. With data of 884,105 autosomal common single-nucleotide polymorphisms (SNPs), three methods were applied to test for pleiotropy between MDD and EA: (i) genetic profile risk scores (GPRS) derived from training data for EA (independent meta-analysis on ~120,000 subjects) and MDD (using a 10-fold leave-one-out procedure in the current sample), (ii) bivariate genomic-relationship-matrix restricted maximum likelihood (GREML) and (iii) SNP effect concordance analysis (SECA). With these methods, we found (i) that the EA-GPRS did not predict MDD status, and MDD-GPRS did not predict EA, (ii) a weak negative genetic correlation with bivariate GREML analyses, but this correlation was not consistently significant, (iii) no evidence for concordance of MDD and EA SNP effects with SECA analysis. To conclude, our study confirms an association of lower EA and MDD risk, but this association was not because of measurable pleiotropic genetic effects, which suggests that environmental factors could be involved, for example, socioeconomic status.

Examining the role of common genetic variants on alcohol, tobacco, cannabis and illicit drug dependence: genetics of vulnerability to drug dependence

BACKGROUND AND AIMS

Twin and family studies suggest that genetic influences are shared across substances of abuse. However, despite evidence of heritability, genome-wide association and candidate gene studies have indicated numerous markers of limited effects, suggesting that much of the heritability remains missing. We estimated (1) the aggregate effect of common single nucleotide polymorphisms (SNPs) on multiple indicators of comorbid drug problems that are typically employed across community and population-based samples, and (2) the genetic covariance across these measures.

PARTICIPANTS

A total of 2596 unrelated subjects from the Study of Addiction: Genetics and Environment provided information on alcohol, tobacco, cocaine, cannabis and other illicit substance dependence. Phenotypic measures included: (1) a factor score based on DSM-IV drug dependence diagnoses (DD), (2) a factor score based on problem use (PU; i.e. 1+ DSM-IV symptoms) and (3) dependence vulnerability (DV; a ratio of DSM-IV symptoms to the number of substances used).

FINDINGS

Univariate and bivariate genome-wide complex trait analyses of this selected sample indicated that common SNPs explained 25-36% of the variance across measures, with DD and DV having the largest effects [h(2) SNP (standard error) = 0.36 (0.13) and 0.33 (0.13), respectively; PU = 0.25 (0.13)]. Genetic effects were shared across the three phenotypic measures of comorbid drug problems [rDD-PU = 0.92 (0.08), rDD-DV = 0.97 (0.08) and rPU-DV = 0.96 (0.07)].

CONCLUSION

At least 20% of the variance in the generalized vulnerability to substance dependence is attributable to common single nucleotide polymorphisms. The additive effect of common single nucleotide polymorphisms is shared across important indicators of comorbid drug problems.

Efficient multiple-trait association and estimation of genetic correlation using the matrix-variate linear mixed model

Multiple-trait association mapping, in which multiple traits are used simultaneously in the identification of genetic variants affecting those traits, has recently attracted interest. One class of approaches for this problem builds on classical variance component methodology, utilizing a multitrait version of a linear mixed model. These approaches both increase power and provide insights into the genetic architecture of multiple traits. In particular, it is possible to estimate the genetic correlation, which is a measure of the portion of the total correlation between traits that is due to additive genetic effects. Unfortunately, the practical utility of these methods is limited since they are computationally intractable for large sample sizes. In this article, we introduce a reformulation of the multiple-trait association mapping approach by defining the matrix-variate linear mixed model. Our approach reduces the computational time necessary to perform maximum-likelihood inference in a multiple-trait model by utilizing a data transformation. By utilizing a well-studied human cohort, we show that our approach provides more than a 10-fold speedup, making multiple-trait association feasible in a large population cohort on the genome-wide scale. We take advantage of the efficiency of our approach to analyze gene expression data. By decomposing gene coexpression into a genetic and environmental component, we show that our method provides fundamental insights into the nature of coexpressed genes. An implementation of this method is available at http://genetics.cs.ucla.edu/mvLMM.

Age- and gender-dependent heterogeneous proportion of variation explained by SNPs in quantitative traits reflecting human health

Age-related effects are often included as covariates in the analytical model for genome-wide association analysis of quantitative traits reflecting human health. Nevertheless, previous studies have hardly examined the effects of age on the proportion of variation explained by single nucleotide polymorphisms (PVSNP) in these traits. In this study, the PVSNP estimates of body mass index (BMI), waist-to-hip ratio, pulse pressure, high-density lipoprotein cholesterol level, triglyceride level (TG), low-density lipoprotein cholesterol level, and glucose level were obtained from Korean consortium metadata partitioned by gender or by age. Restricted maximum likelihood estimates of the PVSNP were obtained in a mixed model framework. Previous studies using pedigree data suggested possible differential heritability of certain traits with regard to gender, which we observed in our current study (BMI and TG; P < 0.05). However, the PVSNP analysis based on age revealed that, with respect to every trait tested, individuals aged 40 to 49 exhibited significantly lower PVSNP estimates than individuals aged 50 to 59 or 60 to 69 (P < 0.05). The consistent heterogeneous PVSNP with respect to age may be due to degenerated genetic functions in individuals between the ages of 50 and 69. Our results suggest the genetic mechanism of age- and gender-dependent PVSNP of quantitative traits related to human health should be further examined.

Joint analysis of psychiatric disorders increases accuracy of risk prediction for schizophrenia, bipolar disorder, and major depressive disorder

Genetic risk prediction has several potential applications in medical research and clinical practice and could be used, for example, to stratify a heterogeneous population of patients by their predicted genetic risk. However, for polygenic traits, such as psychiatric disorders, the accuracy of risk prediction is low. Here we use a multivariate linear mixed model and apply multi-trait genomic best linear unbiased prediction for genetic risk prediction. This method exploits correlations between disorders and simultaneously evaluates individual risk for each disorder. We show that the multivariate approach significantly increases the prediction accuracy for schizophrenia, bipolar disorder, and major depressive disorder in the discovery as well as in independent validation datasets. By grouping SNPs based on genome annotation and fitting multiple random effects, we show that the prediction accuracy could be further improved. The gain in prediction accuracy of the multivariate approach is equivalent to an increase in sample size of 34% for schizophrenia, 68% for bipolar disorder, and 76% for major depressive disorders using single trait models. Because our approach can be readily applied to any number of GWAS datasets of correlated traits, it is a flexible and powerful tool to maximize prediction accuracy. With current sample size, risk predictors are not useful in a clinical setting but already are a valuable research tool, for example in experimental designs comparing cases with high and low polygenic risk.

Genetics and intelligence differences: five special findings

Intelligence is a core construct in differential psychology and behavioural genetics, and should be so in cognitive neuroscience. It is one of the best predictors of important life outcomes such as education, occupation, mental and physical health and illness, and mortality. Intelligence is one of the most heritable behavioural traits. Here, we highlight five genetic findings that are special to intelligence differences and that have important implications for its genetic architecture and for gene-hunting expeditions. (i) The heritability of intelligence increases from about 20% in infancy to perhaps 80% in later adulthood. (ii) Intelligence captures genetic effects on diverse cognitive and learning abilities, which correlate phenotypically about 0.30 on average but correlate genetically about 0.60 or higher. (iii) Assortative mating is greater for intelligence (spouse correlations ~0.40) than for other behavioural traits such as personality and psychopathology (~0.10) or physical traits such as height and weight (~0.20). Assortative mating pumps additive genetic variance into the population every generation, contributing to the high narrow heritability (additive genetic variance) of intelligence. (iv) Unlike psychiatric disorders, intelligence is normally distributed with a positive end of exceptional performance that is a model for 'positive genetics'. (v) Intelligence is associated with education and social class and broadens the causal perspectives on how these three inter-correlated variables contribute to social mobility, and health, illness and mortality differences. These five findings arose primarily from twin studies. They are being confirmed by the first new quantitative genetic technique in a century-Genome-wide Complex Trait Analysis (GCTA)-which estimates genetic influence using genome-wide genotypes in large samples of unrelated individuals. Comparing GCTA results to the results of twin studies reveals important insights into the genetic architecture of intelligence that are relevant to attempts to narrow the 'missing heritability' gap.

Genetic architecture of white matter hyperintensities differs in hypertensive and nonhypertensive ischemic stroke

Supplemental Digital Content is available in the text.

Background and Purpose—

Epidemiological studies suggest that white matter hyperintensities (WMH) are extremely heritable, but the underlying genetic variants are largely unknown. Pathophysiological heterogeneity is known to reduce the power of genome-wide association studies (GWAS). Hypertensive and nonhypertensive individuals with WMH might have different underlying pathologies. We used GWAS data to calculate the variance in WMH volume (WMHV) explained by common single nucleotide polymorphisms (SNPs) as a measure of heritability (SNP heritability [H_SNP]) and tested the hypothesis that WMH heritability differs between hypertensive and nonhypertensive individuals.

Methods—

WMHV was measured on MRI in the stroke-free cerebral hemisphere of 2336 ischemic stroke cases with GWAS data. After adjustment for age and intracranial volume, we determined which cardiovascular risk factors were independent predictors of WMHV. Using the genome-wide complex trait analysis tool to estimate H_SNP for WMHV overall and within subgroups stratified by risk factors found to be significant in multivariate analyses.

Results—

A significant proportion of the variance of WMHV was attributable to common SNPs after adjustment for significant risk factors (H_SNP=0.23; P=0.0026). H_SNP estimates were higher among hypertensive individuals (H_SNP=0.45; P=7.99×10⁻⁵); this increase was greater than expected by chance (P=0.012). In contrast, estimates were lower, and nonsignificant, in nonhypertensive individuals (H_SNP=0.13; P=0.13).

Conclusions—

A quarter of variance is attributable to common SNPs, but this estimate was greater in hypertensive individuals. These findings suggest that the genetic architecture of WMH in ischemic stroke differs between hypertensives and nonhypertensives. Future WMHV GWAS studies may gain power by accounting for this interaction.

Heritability and genome-wide analyses of problematic peer relationships during childhood and adolescence

Peer behaviour plays an important role in the development of social adjustment, though little is known about its genetic architecture. We conducted a twin study combined with a genome-wide complex trait analysis (GCTA) and a genome-wide screen to characterise genetic influences on problematic peer behaviour during childhood and adolescence. This included a series of longitudinal measures (parent-reported Strengths-and-Difficulties Questionnaire) from a UK population-based birth-cohort (ALSPAC, 4-17 years), and a UK twin sample (TEDS, 4-11 years). Longitudinal twin analysis (TEDS; N ≤ 7,366 twin pairs) showed that peer problems in childhood are heritable (4-11 years, 0.60 < twin-h(2) ≤ 0.71) but genetically heterogeneous from age to age (4-11 years, twin-r(g) = 0.30). GCTA (ALSPAC: N ≤ 5,608, TEDS: N ≤ 2,691) provided furthermore little support for the contribution of measured common genetic variants during childhood (4-12 years, 0.02 < GCTA-h(2)(Meta) ≤ 0.11) though these influences become stronger in adolescence (13-17 years, 0.14 < GCTA-h (2)(ALSPAC) ≤ 0.27). A subsequent cross-sectional genome-wide screen in ALSPAC (N ≤ 6,000) focussed on peer problems with the highest GCTA-heritability (10, 13 and 17 years, 0.0002 < GCTA-P ≤ 0.03). Single variant signals (P ≤ 10(-5)) were followed up in TEDS (N ≤ 2835, 9 and 11 years) and, in search for autism quantitative trait loci, explored within two autism samples (AGRE: N Pedigrees = 793; ACC: N Cases = 1,453/N Controls = 7,070). There was, however, no evidence for association in TEDS and little evidence for an overlap with the autistic continuum. In summary, our findings suggest that problematic peer relationships are heritable but genetically complex and heterogeneous from age to age, with an increase in common measurable genetic variation during adolescence.

Heritability of pulmonary function estimated from genome-wide SNPs in healthy Japanese adults

BACKGROUND

Pulmonary function is a heritable trait, and recent genome-wide association studies (GWASs) have identified a number of loci influencing the trait. Genome-wide Complex Trait Analysis (GCTA) is a novel method provided by a software package that estimates the total additive genetic influence caused by common single nucleotide polymorphisms (SNPs) on whole-genome arrays. We conducted a GWAS and assessed the heritability of pulmonary function in an adult Japanese population using this approach.

METHODS

We initially conducted a GWAS on %forced vital capacity (FVC), %forced expiratory volume (FEV1) and FEV1/FVC in healthy Japanese adults (N=967). We then examined the heritability of these traits using GCTA with a total of 480,026 SNPs. We also estimated the genetic impact of the 24 genes identified as susceptibility genes to FEV1/FVC in six previous GWASs on the heritability of FEV1/FVC in the Japanese population.

RESULTS

The heritabilities for %FVC, %FEV1, and FEV1/FVC were 71.2%, 51.9% and 41.6%, respectively. These results corresponded to previous heritability estimates for pulmonary function obtained by GCTA or by twin studies. The 24 previously reported pulmonary function genes accounted for 4.3-12.0% of the entire estimated heritability of FEV1/FVC.

CONCLUSIONS

This study demonstrated that the heritability of pulmonary function traits can be explained by the additive effects of multiple common SNPs in healthy Japanese adults. The pulmonary function genes reported in previous GWASs of non-Japanese populations showed a definite impact of the genes on FEV1/FVC, thus indicating the presence of common pathways related to this trait beyond ethnicity.

GPA: a statistical approach to prioritizing GWAS results by integrating pleiotropy and annotation

Results from Genome-Wide Association Studies (GWAS) have shown that complex diseases are often affected by many genetic variants with small or moderate effects. Identifications of these risk variants remain a very challenging problem. There is a need to develop more powerful statistical methods to leverage available information to improve upon traditional approaches that focus on a single GWAS dataset without incorporating additional data. In this paper, we propose a novel statistical approach, GPA (Genetic analysis incorporating Pleiotropy and Annotation), to increase statistical power to identify risk variants through joint analysis of multiple GWAS data sets and annotation information because: (1) accumulating evidence suggests that different complex diseases share common risk bases, i.e., pleiotropy; and (2) functionally annotated variants have been consistently demonstrated to be enriched among GWAS hits. GPA can integrate multiple GWAS datasets and functional annotations to seek association signals, and it can also perform hypothesis testing to test the presence of pleiotropy and enrichment of functional annotation. Statistical inference of the model parameters and SNP ranking is achieved through an EM algorithm that can handle genome-wide markers efficiently. When we applied GPA to jointly analyze five psychiatric disorders with annotation information, not only did GPA identify many weak signals missed by the traditional single phenotype analysis, but it also revealed relationships in the genetic architecture of these disorders. Using our hypothesis testing framework, statistically significant pleiotropic effects were detected among these psychiatric disorders, and the markers annotated in the central nervous system genes and eQTLs from the Genotype-Tissue Expression (GTEx) database were significantly enriched. We also applied GPA to a bladder cancer GWAS data set with the ENCODE DNase-seq data from 125 cell lines. GPA was able to detect cell lines that are biologically more relevant to bladder cancer. The R implementation of GPA is currently available at http://dongjunchung.github.io/GPA/.