Defining the role of common variation in the genomic and biological architecture of adult human height

Temporal trends, regional variation and socio-economic differences in height, BMI and body proportions among German conscripts, 1956-2010

OBJECTIVE

We analyse temporal trends and regional variation among the most recent available anthropometric data from German conscription in the years 2008-2010 and their historical contextualization since 1956. Design/setting/subjects The overall sample included German conscripts (N 13 857 313) from 1956 to 2010.

RESULTS

German conscripts changed from growing in height to growing in breadth. Over the analysed 54 years, average height of 19-year-old conscripts increased by 6·5 cm from 173·5 cm in 1956 (birth year 1937) to 180·0 cm in 2010 (birth year 1991). This increase plateaued since the 1990s (1970s birth years). The increase in average weight, however, did not lessen during the last two decades but increased in two steps: at the end of the 1980s and after 1999. The weight and BMI distributions became increasingly right-skewed, the prevalence of overweight and obesity increased from 11·6 % and 2·1 % in 1984 to 19·9 % and 8·5 % in 2010, respectively. The north-south gradient in height (north = taller) persisted during our observations. Height and weight of conscripts from East Germany matched the German average between the early 1990s and 2009. Between the 1980s and the early 1990s, the average chest circumference increased, the average difference between chest circumference when inhaling and exhaling decreased, as did leg length relative to trunk length.

CONCLUSIONS

Measuring anthropometric data for military conscripts yielded year-by-year monitoring of the health status of young men at a proscribed age. Such findings contribute to a more precise identification of groups at risk and thus help with further studies and to target interventions.

Association Study of GWAS-Derived Loci with Height in Brazilian Children: Importance of MAP3K3, MMP24 and IGF1R Polymorphisms for Height Variation

BACKGROUND/AIM

The single nucleotide polymorphisms (SNPs) rs2282978 (CDK6), rs2425019 (MMP24), rs8081612 (MAP3K3), rs2871865 (IGF1R) and rs3782415 (SOCS2) were among the SNPs most strongly associated with height in a meta-analysis of 47 genome-wide association studies (GWAS) involving 114,223 adults from six ethnic groups. The present study aimed to examine associations between these SNPs and height in Brazilian children.

METHODS

Cross-sectional heights of 1,008 healthy unrelated 4.4- to 9.7-year-old children were evaluated. All genotypes were determined by allele-specific polymerase chain reactions. Height standard deviation scores (SDS) were generated for this population and regressed on allele counts. Linear regressions were performed to estimate the effect of individual SNPs or a polygenic allelic score on height.

RESULTS

The T allele of rs8081612 (MAP3K3), the C allele of rs2871865 (IGF1R) and the G allele of rs2425019 (MMP24) were significantly associated with a 0.091-SDS greater height (95% CI 0.089-0.093, p = 0.001) by polygenic analysis. The mean height SDS difference between children with 2 'tall' alleles and children with 4 'tall' alleles was 0.24 SDS (95% CI 0.05-0.43, p = 0.01). The observed allelic effect is consistent with that found in previous GWAS.

CONCLUSIONS

Polymorphisms in MAP3K3, MMP24 and IGF1R act additively on height in children of an admixed population. These results demonstrate the importance of these loci for children's height.

Bipartite Community Structure of eQTLs

Genome Wide Association Studies (GWAS) and expression quantitative trait locus (eQTL) analyses have identified genetic associations with a wide range of human phenotypes. However, many of these variants have weak effects and understanding their combined effect remains a challenge. One hypothesis is that multiple SNPs interact in complex networks to influence functional processes that ultimately lead to complex phenotypes, including disease states. Here we present CONDOR, a method that represents both cis- and trans-acting SNPs and the genes with which they are associated as a bipartite graph and then uses the modular structure of that graph to place SNPs into a functional context. In applying CONDOR to eQTLs in chronic obstructive pulmonary disease (COPD), we found the global network “hub” SNPs were devoid of disease associations through GWAS. However, the network was organized into 52 communities of SNPs and genes, many of which were enriched for genes in specific functional classes. We identified local hubs within each community (“core SNPs”) and these were enriched for GWAS SNPs for COPD and many other diseases. These results speak to our intuition: rather than single SNPs influencing single genes, we see groups of SNPs associated with the expression of families of functionally related genes and that disease SNPs are associated with the perturbation of those functions. These methods are not limited in their application to COPD and can be used in the analysis of a wide variety of disease processes and other phenotypic traits.

Large-scale studies have identified thousands of genetic variants associated with different phenotypes without explaining their function. Expression quantitative trait locus analysis associates the compendium of genetic variants with expression levels of individual genes, providing the opportunity to link those variants to functions. But the complexity of those associations has caused most analyses to focus solely on genetic variants immediately adjacent to the genes they may influence. We describe a method that embraces the complexity, representing all variant-gene associations as a bipartite graph. The graph contains highly modular, functional communities in which disease-associated variants emerge as those likely to perturb the structure of the network and the function of the genes in these communities.

Functional Analysis of Mouse G6pc1 Mutations Using a Novel In Situ Assay for Glucose-6-Phosphatase Activity and the Effect of Mutations in Conserved Human G6PC1/G6PC2 Amino Acids on G6PC2 Protein Expression

Elevated fasting blood glucose (FBG) has been associated with increased risk for development of type 2 diabetes. Single nucleotide polymorphisms (SNPs) in G6PC2 are the most important common determinants of variations in FBG in humans. Studies using G6pc2 knockout mice suggest that G6pc2 regulates the glucose sensitivity of insulin secretion. G6PC2 and the related G6PC1 and G6PC3 genes encode glucose-6-phosphatase catalytic subunits. This study describes a functional analysis of 22 non-synonymous G6PC2 SNPs, that alter amino acids that are conserved in human G6PC1, mouse G6pc1 and mouse G6pc2, with the goal of identifying variants that potentially affect G6PC2 activity/expression. Published data suggest strong conservation of catalytically important amino acids between all four proteins and the related G6PC3 isoform. Because human G6PC2 has very low glucose-6-phosphatase activity we used an indirect approach, examining the effect of these SNPs on mouse G6pc1 activity. Using a novel in situ functional assay for glucose-6-phosphatase activity we demonstrate that the amino acid changes associated with the human G6PC2 rs144254880 (Arg79Gln), rs149663725 (Gly114Arg) and rs2232326 (Ser324Pro) SNPs reduce mouse G6pc1 enzyme activity without affecting protein expression. The Arg79Gln variant alters an amino acid mutation of which, in G6PC1, has previously been shown to cause glycogen storage disease type 1a. We also demonstrate that the rs368382511 (Gly8Glu), rs138726309 (His177Tyr), rs2232323 (Tyr207Ser) rs374055555 (Arg293Trp), rs2232326 (Ser324Pro), rs137857125 (Pro313Leu) and rs2232327 (Pro340Leu) SNPs confer decreased G6PC2 protein expression. In summary, these studies identify multiple G6PC2 variants that have the potential to be associated with altered FBG in humans.

Fast set-based association analysis using summary data from GWAS identifies novel gene loci for human complex traits

We propose a method (fastBAT) that performs a fast set-based association analysis for human complex traits using summary-level data from genome-wide association studies (GWAS) and linkage disequilibrium (LD) data from a reference sample with individual-level genotypes. We demonstrate using simulations and analyses of real datasets that fastBAT is more accurate and orders of magnitude faster than the prevailing methods. Using fastBAT, we analyze summary data from the latest meta-analyses of GWAS on 150,064-339,224 individuals for height, body mass index (BMI), and schizophrenia. We identify 6 novel gene loci for height, 2 for BMI, and 3 for schizophrenia at PfastBAT < 5 × 10(-8). The gain of power is due to multiple small independent association signals at these loci (e.g. the THRB and FOXP1 loci for schizophrenia). The method is general and can be applied to GWAS data for all complex traits and diseases in humans and to such data in other species.

Association between Adult Height and Risk of Colorectal, Lung, and Prostate Cancer: Results from Meta-analyses of Prospective Studies and Mendelian Randomization Analyses

BACKGROUND

Observational studies examining associations between adult height and risk of colorectal, prostate, and lung cancers have generated mixed results. We conducted meta-analyses using data from prospective cohort studies and further carried out Mendelian randomization analyses, using height-associated genetic variants identified in a genome-wide association study (GWAS), to evaluate the association of adult height with these cancers.

METHODS AND FINDINGS

A systematic review of prospective studies was conducted using the PubMed, Embase, and Web of Science databases. Using meta-analyses, results obtained from 62 studies were summarized for the association of a 10-cm increase in height with cancer risk. Mendelian randomization analyses were conducted using summary statistics obtained for 423 genetic variants identified from a recent GWAS of adult height and from a cancer genetics consortium study of multiple cancers that included 47,800 cases and 81,353 controls. For a 10-cm increase in height, the summary relative risks derived from the meta-analyses of prospective studies were 1.12 (95% CI 1.10, 1.15), 1.07 (95% CI 1.05, 1.10), and 1.06 (95% CI 1.02, 1.11) for colorectal, prostate, and lung cancers, respectively. Mendelian randomization analyses showed increased risks of colorectal (odds ratio [OR] = 1.58, 95% CI 1.14, 2.18) and lung cancer (OR = 1.10, 95% CI 1.00, 1.22) associated with each 10-cm increase in genetically predicted height. No association was observed for prostate cancer (OR = 1.03, 95% CI 0.92, 1.15). Our meta-analysis was limited to published studies. The sample size for the Mendelian randomization analysis of colorectal cancer was relatively small, thus affecting the precision of the point estimate.

CONCLUSIONS

Our study provides evidence for a potential causal association of adult height with the risk of colorectal and lung cancers and suggests that certain genetic factors and biological pathways affecting adult height may also affect the risk of these cancers.

The Longitudinal Aging Study Amsterdam: cohort update 2016 and major findings

The Longitudinal Aging Study Amsterdam (LASA) is an ongoing longitudinal study of older adults in the Netherlands, which started in 1992. LASA is focused on the determinants, trajectories and consequences of physical, cognitive, emotional and social functioning. The study is based on a nationally representative sample of older adults aged 55 years and over. The findings of the LASA study have been reported in over 450 publications so far (see www.lasa-vu.nl ). In this article we describe the background and the design of the LASA study, and provide an update of the methods. In addition, we provide a summary of the major findings from the period 2011-2015.

Adiposity-Dependent Regulatory Effects on Multi-tissue Transcriptomes

Obesity is a global epidemic that is causally associated with a range of diseases, including type 2 diabetes and cardiovascular disease, at the population-level. However, there is marked heterogeneity in obesity-related outcomes among individuals. This might reflect genotype-dependent responses to adiposity. Given that adiposity, measured by BMI, is associated with widespread changes in gene expression and regulatory variants mediate the majority of known complex trait loci, we sought to identify gene-by-BMI (G × BMI) interactions on the regulation of gene expression in a multi-tissue RNA-sequencing (RNA-seq) dataset from the TwinsUK cohort (n = 856). At a false discovery rate of 5%, we identified 16 cis G × BMI interactions (top cis interaction: CHURC1, rs7143432, p = 2.0 × 10(-12)) and one variant regulating 53 genes in trans (top trans interaction: ZNF423, rs3851570, p = 8.2 × 10(-13)), all in adipose tissue. The interactions were adipose-specific and enriched for variants overlapping adipocyte enhancers, and regulated genes were enriched for metabolic and inflammatory processes. We replicated a subset of the interactions in an independent adipose RNA-seq dataset (deCODE genetics, n = 754). We also confirmed the interactions with an alternate measure of obesity, dual-energy X-ray absorptiometry (DXA)-derived visceral-fat-volume measurements, in a subset of TwinsUK individuals (n = 682). The identified G × BMI regulatory effects demonstrate the dynamic nature of gene regulation and reveal a functional mechanism underlying the heterogeneous response to obesity. Additionally, we have provided a web browser allowing interactive exploration of the dataset, including of association between expression, BMI, and G × BMI regulatory effects in four tissues.

Opportunities and challenges of big data for the social sciences: The case of genomic data

In this paper, we draw attention to one unique and valuable source of big data, genomic data, by demonstrating the opportunities they provide to social scientists. We discuss different types of large-scale genomic data and recent advances in statistical methods and computational infrastructure used to address challenges in managing and analyzing such data. We highlight how these data and methods can be used to benefit social science research.

Phenome-wide analysis of genome-wide polygenic scores

Genome-wide polygenic scores (GPS), which aggregate the effects of thousands of DNA variants from genome-wide association studies (GWAS), have the potential to make genetic predictions for individuals. We conducted a systematic investigation of associations between GPS and many behavioral traits, the behavioral phenome. For 3152 unrelated 16-year-old individuals representative of the United Kingdom, we created 13 GPS from the largest GWAS for psychiatric disorders (for example, schizophrenia, depression and dementia) and cognitive traits (for example, intelligence, educational attainment and intracranial volume). The behavioral phenome included 50 traits from the domains of psychopathology, personality, cognitive abilities and educational achievement. We examined phenome-wide profiles of associations for the entire distribution of each GPS and for the extremes of the GPS distributions. The cognitive GPS yielded stronger predictive power than the psychiatric GPS in our UK-representative sample of adolescents. For example, education GPS explained variation in adolescents' behavior problems (~0.6%) and in educational achievement (~2%) but psychiatric GPS were associated with neither. Despite the modest effect sizes of current GPS, quantile analyses illustrate the ability to stratify individuals by GPS and opportunities for research. For example, the highest and lowest septiles for the education GPS yielded a 0.5 s.d. difference in mean math grade and a 0.25 s.d. difference in mean behavior problems. We discuss the usefulness and limitations of GPS based on adult GWAS to predict genetic propensities earlier in development.

Meta-Analyses of Genome-Wide Association Data Hold New Promise for Addiction Genetics

Meta-analyses of genome-wide association study data have begun to lead to promising new discoveries for behavioral and psychiatrically relevant phenotypes (e.g., schizophrenia, educational attainment). We outline how this methodology can similarly lead to novel discoveries in genomic studies of substance use disorders, and discuss challenges that will need to be overcome to accomplish this goal. We illustrate our approach with the work of the newly established Substance Use Disorders workgroup of the Psychiatric Genomics Consortium.

Across-cohort QC analyses of GWAS summary statistics from complex traits

Genome-wide association studies (GWASs) have been successful in discovering SNP trait associations for many quantitative traits and common diseases. Typically, the effect sizes of SNP alleles are very small and this requires large genome-wide association meta-analyses (GWAMAs) to maximize statistical power. A trend towards ever-larger GWAMA is likely to continue, yet dealing with summary statistics from hundreds of cohorts increases logistical and quality control problems, including unknown sample overlap, and these can lead to both false positive and false negative findings. In this study, we propose four metrics and visualization tools for GWAMA, using summary statistics from cohort-level GWASs. We propose methods to examine the concordance between demographic information, and summary statistics and methods to investigate sample overlap. (I) We use the population genetics F_st statistic to verify the genetic origin of each cohort and their geographic location, and demonstrate using GWAMA data from the GIANT Consortium that geographic locations of cohorts can be recovered and outlier cohorts can be detected. (II) We conduct principal component analysis based on reported allele frequencies, and are able to recover the ancestral information for each cohort. (III) We propose a new statistic that uses the reported allelic effect sizes and their standard errors to identify significant sample overlap or heterogeneity between pairs of cohorts. (IV) To quantify unknown sample overlap across all pairs of cohorts, we propose a method that uses randomly generated genetic predictors that does not require the sharing of individual-level genotype data and does not breach individual privacy.

Evidence of Polygenic Adaptation in the Systems Genetics of Anthropometric Traits

Many signals of natural selection have been identified in the human genome. However, except for some single-locus mechanisms, most molecular processes generating these adaptation signals are still unknown. We developed an approach that integrates datasets related to genome-wide association studies (GWAS) with information about systems biology and genetic signatures of natural selection to identify evidence of polygenic adaptation. Specifically, we focused on five anthropometric measurements: body mass index (BMI), height, waist-to-hip ratio adjusted for BMI (WHR), and waist circumference adjusted for BMI (WC), and sex differences for WHR and WC. We performed an enrichment analysis for signals of natural selection in protein interaction networks associated with anthropometric traits in European populations. The adaptation signals-enriched gene networks associated highlighted epistatic interactions in the context of polygenic selection for the investigated traits. These polygenic mechanisms indicated intriguing selective mechanisms related to the anthropometric traits: adult locomotory behavior for BMI, infection resistance for height, interplay between lipid transport and immune systems for WHR, and female-specific polygenic adaptation for WHR and WC. In conclusion, we observed evidence of polygenic adaptation in the context of systems genetics of anthropometric traits that indicates polygenic mechanisms related to the natural selection in European populations.

Nationwide Genomic Study in Denmark Reveals Remarkable Population Homogeneity

Denmark has played a substantial role in the history of Northern Europe. Through a nationwide scientific outreach initiative, we collected genetic and anthropometrical data from ∼800 high school students and used them to elucidate the genetic makeup of the Danish population, as well as to assess polygenic predictions of phenotypic traits in adolescents. We observed remarkable homogeneity across different geographic regions, although we could still detect weak signals of genetic structure reflecting the history of the country. Denmark presented genomic affinity with primarily neighboring countries with overall resemblance of decreasing weight from Britain, Sweden, Norway, Germany, and France. A Polish admixture signal was detected in Zealand and Funen, and our date estimates coincided with historical evidence of Wend settlements in the south of Denmark. We also observed considerably diverse demographic histories among Scandinavian countries, with Denmark having the smallest current effective population size compared to Norway and Sweden. Finally, we found that polygenic prediction of self-reported adolescent height in the population was remarkably accurate (R² = 0.639 ± 0.015). The high homogeneity of the Danish population could render population structure a lesser concern for the upcoming large-scale gene-mapping studies in the country.

Gene set analysis for interpreting genetic studies

Interpretation of genome-wide association study (GWAS) results is lacking behind the discovery of new genetic associations. Consequently, there is an urgent need for data-driven methods for interpreting genetic association studies. Gene set analysis (GSA) can identify aetiologic pathways and functional annotations and may hence point towards novel biological insights. However, despite the growing availability of GSA tools, the sizeable amount of variants identified for a vast number of complex traits, and many irrefutably trait-associated gene sets, the gap between discovery and interpretation remains. More efficient interpretation requires more complete and consistent gene set representations of biological pathways, phenotypes and functional annotations. In this review, I examine different types of gene sets, discuss how inconsistencies in gene set definitions impact GSA, describe how GSA has helped to elucidate biology and outline potential future directions.

Common Genetic Polymorphisms Influence Blood Biomarker Measurements in COPD

Implementing precision medicine for complex diseases such as chronic obstructive lung disease (COPD) will require extensive use of biomarkers and an in-depth understanding of how genetic, epigenetic, and environmental variations contribute to phenotypic diversity and disease progression. A meta-analysis from two large cohorts of current and former smokers with and without COPD [SPIROMICS (N = 750); COPDGene (N = 590)] was used to identify single nucleotide polymorphisms (SNPs) associated with measurement of 88 blood proteins (protein quantitative trait loci; pQTLs). PQTLs consistently replicated between the two cohorts. Features of pQTLs were compared to previously reported expression QTLs (eQTLs). Inference of causal relations of pQTL genotypes, biomarker measurements, and four clinical COPD phenotypes (airflow obstruction, emphysema, exacerbation history, and chronic bronchitis) were explored using conditional independence tests. We identified 527 highly significant (p < 8 X 10-10) pQTLs in 38 (43%) of blood proteins tested. Most pQTL SNPs were novel with low overlap to eQTL SNPs. The pQTL SNPs explained >10% of measured variation in 13 protein biomarkers, with a single SNP (rs7041; p = 10-392) explaining 71%-75% of the measured variation in vitamin D binding protein (gene = GC). Some of these pQTLs [e.g., pQTLs for VDBP, sRAGE (gene = AGER), surfactant protein D (gene = SFTPD), and TNFRSF10C] have been previously associated with COPD phenotypes. Most pQTLs were local (cis), but distant (trans) pQTL SNPs in the ABO blood group locus were the top pQTL SNPs for five proteins. The inclusion of pQTL SNPs improved the clinical predictive value for the established association of sRAGE and emphysema, and the explanation of variance (R2) for emphysema improved from 0.3 to 0.4 when the pQTL SNP was included in the model along with clinical covariates. Causal modeling provided insight into specific pQTL-disease relationships for airflow obstruction and emphysema. In conclusion, given the frequency of highly significant local pQTLs, the large amount of variance potentially explained by pQTL, and the differences observed between pQTLs and eQTLs SNPs, we recommend that protein biomarker-disease association studies take into account the potential effect of common local SNPs and that pQTLs be integrated along with eQTLs to uncover disease mechanisms. Large-scale blood biomarker studies would also benefit from close attention to the ABO blood group.

A Track Record on SHOX: From Basic Research to Complex Models and Therapy

SHOX deficiency is the most frequent genetic growth disorder associated with isolated and syndromic forms of short stature. Caused by mutations in the homeobox gene SHOX, its varied clinical manifestations include isolated short stature, Léri-Weill dyschondrosteosis, and Langer mesomelic dysplasia. In addition, SHOX deficiency contributes to the skeletal features in Turner syndrome. Causative SHOX mutations have allowed downstream pathology to be linked to defined molecular lesions. Expression levels of SHOX are tightly regulated, and almost half of the pathogenic mutations have affected enhancers. Clinical severity of SHOX deficiency varies between genders and ranges from normal stature to profound mesomelic skeletal dysplasia. Treatment options for children with SHOX deficiency are available. Two decades of research support the concept of SHOX as a transcription factor that integrates diverse aspects of bone development, growth plate biology, and apoptosis. Due to its absence in mouse, the animal models of choice have become chicken and zebrafish. These models, therefore, together with micromass cultures and primary cell lines, have been used to address SHOX function. Pathway and network analyses have identified interactors, target genes, and regulators. Here, we summarize recent data and give insight into the critical molecular and cellular functions of SHOX in the etiopathogenesis of short stature and limb development.

Endogenous androgen exposures and ischemic heart disease, a separate sample Mendelian randomization study

BACKGROUND

Evolutionary biology suggests growth and reproduction trade-off against longevity. Correspondingly estrogen supplementation failed to increase lifespan. Testosterone supplementation is widely used by older men, although regulators have warned of its cardiovascular risk. No large trial of testosterone exists. We examined how genetic determinants of up-regulation (follicle-stimulating hormone (FSH)) and down-regulation (anti-Müllerian hormone (AMH) and testicular dysgenesis syndrome (TDS)) of mainly the male reproductive system are associated with ischemic heart disease (IHD).

METHODS

Separate sample instrumental variable analysis with genetic instruments, i.e., Mendelian randomization, was used to obtain unconfounded estimates using large case-control studies of coronary artery disease/myocardial infarction (CAD/MI) with extensive genotyping, i.e., CARDIoGRAMplusC4D (64,374 CAD/MI cases, 130,681controls), or CARDIoGRAMplusC4D 1000 Genomes (60,801 cases, 123,504 controls).

RESULTS

Genetically predicted FSH was positively associated with CAD/MI (odds ratio (OR) 1.08, 95% confidence interval (CI) 1.03 to 1.13 per mIU/mL FSH). Genetically predicted AMH and TDS were inversely associated with CAD/MI (OR 0.93, 95% CI 0.87 to 0.998 per ng/mL log AMH and OR 0.89, 95% CI 0.81 to 0.98 per log OR higher risk of TDS).

CONCLUSIONS

As expected from evolutionary biology, genetically predicted FSH, related to higher androgens in men and women, was positively associated with IHD, while genetically predicted AMH and TDS, related to lower androgens in men, were inversely associated with IHD. Androgens might be a modifiable causal factor underlying men's greater vulnerability to IHD, with corresponding implications for use of testosterone supplementation as well as for prevention and treatment of IHD.

Changing Polygenic Penetrance on Phenotypes in the 20(th) Century Among Adults in the US Population

This study evaluates changes in genetic penetrance—defined as the association between an additive polygenic score and its associated phenotype—across birth cohorts. Situating our analysis within recent historical trends in the U.S., we show that, while height and BMI show increasing genotypic penetrance over the course of 20^th Century, education and heart disease show declining genotypic effects. Meanwhile, we find genotypic penetrance to be historically stable with respect to depression. Our findings help inform our understanding of how the genetic and environmental landscape of American society has changed over the past century, and have implications for research which models gene-environment (GxE) interactions, as well as polygenic score calculations in consortia studies that include multiple birth cohorts.

Genotypic Context and Epistasis in Individuals and Populations

Genes encode components of coevolved and interconnected networks. The effect of genotype on phenotype therefore depends on genotypic context through gene interactions known as epistasis. Epistasis is important in predicting phenotype from genotype for an individual. It is also examined in population studies to identify genetic risk factors in complex traits and to predict evolution under selection. Paradoxically, the effects of genotypic context in individuals and populations are distinct and sometimes contradictory. We argue that predicting genotype from phenotype for individuals based on population studies is difficult and, especially in human genetics, likely to result in underestimating the effects of genotypic context.

Genetic evidence for natural selection in humans in the contemporary United States

Recent findings from molecular genetics now make it possible to test directly for natural selection by analyzing whether genetic variants associated with various phenotypes have been under selection. I leverage these findings to construct polygenic scores that use individuals' genotypes to predict their body mass index, educational attainment (EA), glucose concentration, height, schizophrenia, total cholesterol, and (in females) age at menarche. I then examine associations between these scores and fitness to test whether natural selection has been occurring. My study sample includes individuals of European ancestry born between 1931 and 1953 who participated in the Health and Retirement Study, a representative study of the US population. My results imply that natural selection has been slowly favoring lower EA in both females and males, and are suggestive that natural selection may have favored a higher age at menarche in females. For EA, my estimates imply a rate of selection of about -1.5 mo of education per generation (which pales in comparison with the increases in EA observed in contemporary times). Although they cannot be projected over more than one generation, my results provide additional evidence that humans are still evolving-albeit slowly, especially compared with the rapid changes that have occurred over the past few generations due to cultural and environmental factors.

New quality measure for SNP array based CNV detection

MOTIVATION

Only a few large systematic studies have evaluated the impact of copy number variants (CNVs) on common diseases. Several million individuals have been genotyped on single nucleotide variation arrays, which could be used for genome-wide CNVs association studies. However, CNV calls remain prone to false positives and only empirical filtering strategies exist in the literature. To overcome this issue, we defined a new quality score (QS) estimating the probability of a CNV called by PennCNV to be confirmed by other software.

RESULTS

Out-of-sample comparison showed that the correlation between the consensus CNV status and the QS is twice as high as it is for any previously proposed CNV filters. ROC curves displayed an AUC higher than 0.8 and simulations showed an increase up to 20% in statistical power when using QS in comparison to other filtering strategies. Superior performance was confirmed also for alternative consensus CNV definition and through improving known CNV-trait associations.

AVAILABILITY AND IMPLEMENTATION

http://goo.gl/T6yuFM CONTACT: zoltan.kutalik@unil.ch or aurelien@mace@unil.chSupplementary information: Supplementary data are available at Bioinformatics online.

Hunter-gatherer genomics: evolutionary insights and ethical considerations

Hunting and gathering societies currently comprise only a small proportion of all human populations. However, the geographic and environmental diversity of modern hunter-gatherer groups, their inherent dependence on ecological resources, and their connection to patterns of behavior and subsistence that represent the vast majority of human history provide opportunities for scientific research to deliver major insights into the evolutionary history of our species. We review recent evolutionary genomic studies of hunter-gatherers, focusing especially on those that identify and functionally characterize phenotypic adaptations to local environments. We also call attention to specific ethical issues that scientists conducting hunter-gatherer genomics research ought to consider, including potential social and economic tensions between traditionally mobile hunter-gatherers and the land ownership-based nation-states by which they are governed, and the implications of genomic-based evidence of long-term evolutionary associations with particular habitats.

Half dozen of one, six billion of the other: What can small- and large-scale molecular systems biology learn from one another?

Small-scale molecular systems biology, by which we mean the understanding of a how a few parts work together to control a particular biological process, is predicated on the assumption that cellular regulation is arranged in a circuit-like structure. Results from the omics revolution have upset this vision to varying degrees by revealing a high degree of interconnectivity, making it difficult to develop a simple, circuit-like understanding of regulatory processes. We here outline the limitations of the small-scale systems biology approach with examples from research into genetic algorithms, genetics, transcriptional network analysis, and genomics. We also discuss the difficulties associated with deriving true understanding from the analysis of large data sets and propose that the development of new, intelligent, computational tools may point to a way forward. Throughout, we intentionally oversimplify and talk about things in which we have little expertise, and it is likely that many of our arguments are wrong on one level or another. We do believe, however, that developing a true understanding via molecular systems biology will require a fundamental rethinking of our approach, and our goal is to provoke thought along these lines.

Contrasting the Genetic Architecture of 30 Complex Traits from Summary Association Data

Variance-component methods that estimate the aggregate contribution of large sets of variants to the heritability of complex traits have yielded important insights into the genetic architecture of common diseases. Here, we introduce methods that estimate the total trait variance explained by the typed variants at a single locus in the genome (local SNP heritability) from genome-wide association study (GWAS) summary data while accounting for linkage disequilibrium among variants. We applied our estimator to ultra-large-scale GWAS summary data of 30 common traits and diseases to gain insights into their local genetic architecture. First, we found that common SNPs have a high contribution to the heritability of all studied traits. Second, we identified traits for which the majority of the SNP heritability can be confined to a small percentage of the genome. Third, we identified GWAS risk loci where the entire locus explains significantly more variance in the trait than the GWAS reported variants. Finally, we identified loci that explain a significant amount of heritability across multiple traits.

Assessing the genetic architecture of epithelial ovarian cancer histological subtypes

Epithelial ovarian cancer (EOC) is one of the deadliest common cancers. The five most common types of disease are high-grade and low-grade serous, endometrioid, mucinous and clear cell carcinoma. Each of these subtypes present distinct molecular pathogeneses and sensitivities to treatments. Recent studies show that certain genetic variants confer susceptibility to all subtypes while other variants are subtype-specific. Here, we perform an extensive analysis of the genetic architecture of EOC subtypes. To this end, we used data of 10,014 invasive EOC patients and 21,233 controls from the Ovarian Cancer Association Consortium genotyped in the iCOGS array (211,155 SNPs). We estimate the array heritability (attributable to variants tagged on arrays) of each subtype and their genetic correlations. We also look for genetic overlaps with factors such as obesity, smoking behaviors, diabetes, age at menarche and height. We estimated the array heritabilities of high-grade serous disease ([Formula: see text] = 8.8 ± 1.1 %), endometrioid ([Formula: see text] = 3.2 ± 1.6 %), clear cell ([Formula: see text] = 6.7 ± 3.3 %) and all EOC ([Formula: see text] = 5.6 ± 0.6 %). Known associated loci contributed approximately 40 % of the total array heritability for each subtype. The contribution of each chromosome to the total heritability was not proportional to chromosome size. Through bivariate and cross-trait LD score regression, we found evidence of shared genetic backgrounds between the three high-grade subtypes: serous, endometrioid and undifferentiated. Finally, we found significant genetic correlations of all EOC with diabetes and obesity using a polygenic prediction approach.

Detection and interpretation of shared genetic influences on 42 human traits

We performed a scan for genetic variants associated with multiple phenotypes by comparing large genome-wide association studies (GWAS) of 42 traits or diseases. We identified 341 loci (at a false discovery rate of 10%) associated with multiple traits. Several loci are associated with multiple phenotypes; for example, a nonsynonymous variant in the zinc transporter SLC39A8 influences seven of the traits, including risk of schizophrenia (rs13107325: log-transformed odds ratio (log OR) = 0.15, P = 2 × 10(-12)) and Parkinson disease (log OR = -0.15, P = 1.6 × 10(-7)), among others. Second, we used these loci to identify traits that have multiple genetic causes in common. For example, variants associated with increased risk of schizophrenia also tended to be associated with increased risk of inflammatory bowel disease. Finally, we developed a method to identify pairs of traits that show evidence of a causal relationship. For example, we show evidence that increased body mass index causally increases triglyceride levels.

Ovarian Physiology and GWAS: Biobanks, Biology, and Beyond

Ovarian function is central to female fertility, and several genome-wide association studies (GWAS) have been carried out to elucidate the genetic background of traits and disorders that reflect and affect ovarian physiology. While GWAS have been successful in reporting numerous genetic associations and highlighting involved pathways relevant to reproductive aging, for ovarian disorders, such as premature ovarian insufficiency and polycystic ovary syndrome, research has lagged behind due to insufficient study sample size. Novel approaches to study design and analysis methods that help to fit GWAS findings into biological context will improve our knowledge about genetics governing ovarian function in fertility and disease, and provide input for clinical tools and better patient management.

Making sense of big data in health research: Towards an EU action plan

Medicine and healthcare are undergoing profound changes. Whole-genome sequencing and high-resolution imaging technologies are key drivers of this rapid and crucial transformation. Technological innovation combined with automation and miniaturization has triggered an explosion in data production that will soon reach exabyte proportions. How are we going to deal with this exponential increase in data production? The potential of "big data" for improving health is enormous but, at the same time, we face a wide range of challenges to overcome urgently. Europe is very proud of its cultural diversity; however, exploitation of the data made available through advances in genomic medicine, imaging, and a wide range of mobile health applications or connected devices is hampered by numerous historical, technical, legal, and political barriers. European health systems and databases are diverse and fragmented. There is a lack of harmonization of data formats, processing, analysis, and data transfer, which leads to incompatibilities and lost opportunities. Legal frameworks for data sharing are evolving. Clinicians, researchers, and citizens need improved methods, tools, and training to generate, analyze, and query data effectively. Addressing these barriers will contribute to creating the European Single Market for health, which will improve health and healthcare for all Europeans.

Genetic and environmental influences on height from infancy to early adulthood: An individual-based pooled analysis of 45 twin cohorts

Height variation is known to be determined by both genetic and environmental factors, but a systematic description of how their influences differ by sex, age and global regions is lacking. We conducted an individual-based pooled analysis of 45 twin cohorts from 20 countries, including 180,520 paired measurements at ages 1-19 years. The proportion of height variation explained by shared environmental factors was greatest in early childhood, but these effects remained present until early adulthood. Accordingly, the relative genetic contribution increased with age and was greatest in adolescence (up to 0.83 in boys and 0.76 in girls). Comparing geographic-cultural regions (Europe, North-America and Australia, and East-Asia), genetic variance was greatest in North-America and Australia and lowest in East-Asia, but the relative proportion of genetic variation was roughly similar across these regions. Our findings provide further insights into height variation during childhood and adolescence in populations representing different ethnicities and exposed to different environments.

Meta-analysis of 375,000 individuals identifies 38 susceptibility loci for migraine

Migraine is a debilitating neurological disorder affecting around one in seven people worldwide, but its molecular mechanisms remain poorly understood. There is some debate about whether migraine is a disease of vascular dysfunction or a result of neuronal dysfunction with secondary vascular changes. Genome-wide association (GWA) studies have thus far identified 13 independent loci associated with migraine. To identify new susceptibility loci, we carried out a genetic study of migraine on 59,674 affected subjects and 316,078 controls from 22 GWA studies. We identified 44 independent single-nucleotide polymorphisms (SNPs) significantly associated with migraine risk (P < 5 × 10(-8)) that mapped to 38 distinct genomic loci, including 28 loci not previously reported and a locus that to our knowledge is the first to be identified on chromosome X. In subsequent computational analyses, the identified loci showed enrichment for genes expressed in vascular and smooth muscle tissues, consistent with a predominant theory of migraine that highlights vascular etiologies.

Fetal haemoglobin in sickle-cell disease: from genetic epidemiology to new therapeutic strategies

Sickle-cell disease affects millions of individuals worldwide, but the global incidence is concentrated in Africa. The burden of sickle-cell disease is expected to continue to rise over the coming decades, adding to stress on the health infrastructures of many countries. Although the molecular cause of sickle-cell disease has been known for more than half a century, treatment options remain greatly limited. Allogeneic haemopoietic stem-cell transplantation is the only existing cure but is limited to specialised clinical centres and remains inaccessible for most patients. Induction of fetal haemoglobin production is a promising strategy for the treatment of sickle-cell disease. In this Series paper, we review scientific breakthroughs in epidemiology, genetics, and molecular biology that have brought reactivation of fetal haemoglobin to the forefront of sickle-cell disease research. Improved knowledge of the regulation of fetal haemoglobin production in human beings and the development of genome editing technology now support the design of innovative therapies for sickle-cell disease that are based on fetal haemoglobin.

Bias Characterization in Probabilistic Genotype Data and Improved Signal Detection with Multiple Imputation

Missing data are an unavoidable component of modern statistical genetics. Different array or sequencing technologies cover different single nucleotide polymorphisms (SNPs), leading to a complicated mosaic pattern of missingness where both individual genotypes and entire SNPs are sporadically absent. Such missing data patterns cannot be ignored without introducing bias, yet cannot be inferred exclusively from nonmissing data. In genome-wide association studies, the accepted solution to missingness is to impute missing data using external reference haplotypes. The resulting probabilistic genotypes may be analyzed in the place of genotype calls. A general-purpose paradigm, called Multiple Imputation (MI), is known to model uncertainty in many contexts, yet it is not widely used in association studies. Here, we undertake a systematic evaluation of existing imputed data analysis methods and MI. We characterize biases related to uncertainty in association studies, and find that bias is introduced both at the imputation level, when imputation algorithms generate inconsistent genotype probabilities, and at the association level, when analysis methods inadequately model genotype uncertainty. We find that MI performs at least as well as existing methods or in some cases much better, and provides a straightforward paradigm for adapting existing genotype association methods to uncertain data.

Genetic research has been focused at analysis of datapoints that are assumed to be deterministically known. However, the majority of current, high throughput data is only probabilistically known, and proper methods for handing such uncertain genotypes are limited. Here, we build on existing theory from the field of statistics to introduce a general framework for handling probabilistic genotype data obtained through genotype imputation. This framework, called Multiple Imputation, matches or improves upon existing methods for handling uncertainty in basic analysis of genetic association. As opposed to such methods, our work furthermore extends to more advanced analysis, such as mixed-effects models, with no additional complication. Importantly, it generates posterior probabilities of association that are intrinsically weighted by the certainty of the underlying data, a feature unmatched by other existing methods. Multiple Imputation is also fully compatible with meta-analysis. Finally, our analysis of probabilistic genotype data brings into focus the accuracy and unreliability of imputation’s estimated probabilities. Taken together, these results substantially increase the utility of imputed genotypes in statistical genetics, and may have strong implications for analysis of sequencing data moving forward.

Cryptic Genetic Variation in Evolutionary Developmental Genetics

Evolutionary developmental genetics has traditionally been conducted by two groups: Molecular evolutionists who emphasize divergence between species or higher taxa, and quantitative geneticists who study variation within species. Neither approach really comes to grips with the complexities of evolutionary transitions, particularly in light of the realization from genome-wide association studies that most complex traits fit an infinitesimal architecture, being influenced by thousands of loci. This paper discusses robustness, plasticity and lability, phenomena that we argue potentiate major evolutionary changes and provide a bridge between the conceptual treatments of macro- and micro-evolution. We offer cryptic genetic variation and conditional neutrality as mechanisms by which standing genetic variation can lead to developmental system drift and, sheltered within canalized processes, may facilitate developmental transitions and the evolution of novelty. Synthesis of the two dominant perspectives will require recognition that adaptation, divergence, drift and stability all depend on similar underlying quantitative genetic processes-processes that cannot be fully observed in continuously varying visible traits.

Intergenic GWAS SNPs are key components of the spatial and regulatory network for human growth

Meta-analysis of genome-wide association studies has resulted in the identification of hundreds of genetic variants associated with growth and stature. Determining how these genetic variants influence growth is important, but most are non-coding, and there is little understanding of how these variants contribute to adult height. To determine the mechanisms by which human variation contributes to growth, we combined spatial genomic connectivity (high-throughput conformation capture) with functional (gene expression, expression Quantitative Trait Loci) data to determine how non-genic loci associated with infant length, pubertal and adult height and contribute to gene regulatory networks. This approach identified intergenic single-nucleotide polymorphisms (SNPs) ∼85 kb upstream of FBXW11 that spatially connect with distant loci. These regulatory connections are reinforced by evidence of SNP-enhancer effects and altered expression in genes influencing the action of human growth hormone. Functional assays provided evidence for enhancer activity of the intergenic region near FBXW11 that harbors SNP rs12153391, which is associated with an expression Quantitative Trait Loci. Our results suggest that variants in this locus have genome-wide effects as key modifiers of growth (both overgrowth and short stature) acting through a regulatory network. We believe that the genes and pathways connected with this regulatory network are potential targets that could be investigated for diagnostic, prenatal and carrier testing for growth disorders. Finally, the regulatory networks we generated illustrate the power of using existing datasets to interrogate the contribution of intergenic SNPs to common syndromes/diseases.

A method to estimate the contribution of regional genetic associations to complex traits from summary association statistics

Despite considerable efforts, known genetic associations only explain a small fraction of predicted heritability. Regional associations combine information from multiple contiguous genetic variants and can improve variance explained at established association loci. However, regional associations are not easily amenable to estimation using summary association statistics because of sensitivity to linkage disequilibrium (LD). We now propose a novel method, LD Adjusted Regional Genetic Variance (LARGV), to estimate phenotypic variance explained by regional associations using summary statistics while accounting for LD. Our method is asymptotically equivalent to a multiple linear regression model when no interaction or haplotype effects are present. It has several applications, such as ranking of genetic regions according to variance explained or comparison of variance explained by two or more regions. Using height and BMI data from the Health Retirement Study (N = 7,776), we show that most genetic variance lies in a small proportion of the genome and that previously identified linkage peaks have higher than expected regional variance.

Genetic Techniques in the Evaluation of Short Stature

Normal growth is a complex dynamic process dependent on the coordination of multiple factors including genetics, nutrition and hormones that are all working in balance. This chapter will review selected features of commonly utilized genetic techniques such as chromosomal analysis, microarray analysis, targeted gene screening and whole exome sequencing that are being used to identify genes influencing growth. As genetic technologies continue to improve and become more accessible many of these techniques will help to provide a better understanding of mechanisms underlying abnormal growth and will eventually lead to novel management approaches for abnormal growth.

Genetic variants associated with subjective well-being, depressive symptoms, and neuroticism identified through genome-wide analyses

Very few genetic variants have been associated with depression and neuroticism, likely because of limitations on sample size in previous studies. Subjective well-being, a phenotype that is genetically correlated with both of these traits, has not yet been studied with genome-wide data. We conducted genome-wide association studies of three phenotypes: subjective well-being (n = 298,420), depressive symptoms (n = 161,460), and neuroticism (n = 170,911). We identify 3 variants associated with subjective well-being, 2 variants associated with depressive symptoms, and 11 variants associated with neuroticism, including 2 inversion polymorphisms. The two loci associated with depressive symptoms replicate in an independent depression sample. Joint analyses that exploit the high genetic correlations between the phenotypes (|ρ^| ≈ 0.8) strengthen the overall credibility of the findings and allow us to identify additional variants. Across our phenotypes, loci regulating expression in central nervous system and adrenal or pancreas tissues are strongly enriched for association.

Genome-wide association study of cognitive functions and educational attainment in UK Biobank (N=112 151)

People's differences in cognitive functions are partly heritable and are associated with important life outcomes. Previous genome-wide association (GWA) studies of cognitive functions have found evidence for polygenic effects yet, to date, there are few replicated genetic associations. Here we use data from the UK Biobank sample to investigate the genetic contributions to variation in tests of three cognitive functions and in educational attainment. GWA analyses were performed for verbal-numerical reasoning (N=36 035), memory (N=112 067), reaction time (N=111 483) and for the attainment of a college or a university degree (N=111 114). We report genome-wide significant single-nucleotide polymorphism (SNP)-based associations in 20 genomic regions, and significant gene-based findings in 46 regions. These include findings in the ATXN2, CYP2DG, APBA1 and CADM2 genes. We report replication of these hits in published GWA studies of cognitive function, educational attainment and childhood intelligence. There is also replication, in UK Biobank, of SNP hits reported previously in GWA studies of educational attainment and cognitive function. GCTA-GREML analyses, using common SNPs (minor allele frequency>0.01), indicated significant SNP-based heritabilities of 31% (s.e.m.=1.8%) for verbal-numerical reasoning, 5% (s.e.m.=0.6%) for memory, 11% (s.e.m.=0.6%) for reaction time and 21% (s.e.m.=0.6%) for educational attainment. Polygenic score analyses indicate that up to 5% of the variance in cognitive test scores can be predicted in an independent cohort. The genomic regions identified include several novel loci, some of which have been associated with intracranial volume, neurodegeneration, Alzheimer's disease and schizophrenia.

Rejection odds and rejection ratios: A proposal for statistical practice in testing hypotheses

Much of science is (rightly or wrongly) driven by hypothesis testing. Even in situations where the hypothesis testing paradigm is correct, the common practice of basing inferences solely on p-values has been under intense criticism for over 50 years. We propose, as an alternative, the use of the odds of a correct rejection of the null hypothesis to incorrect rejection. Both pre-experimental versions (involving the power and Type I error) and post-experimental versions (depending on the actual data) are considered. Implementations are provided that range from depending only on the p-value to consideration of full Bayesian analysis. A surprise is that all implementations - even the full Bayesian analysis - have complete frequentist justification. Versions of our proposal can be implemented that require only minor modifications to existing practices yet overcome some of their most severe shortcomings.

Obesity and Multiple Sclerosis: A Mendelian Randomization Study

BACKGROUND

Observational studies have reported an association between obesity, as measured by elevated body mass index (BMI), in early adulthood and risk of multiple sclerosis (MS). However, bias potentially introduced by confounding and reverse causation may have influenced these findings. Therefore, we elected to perform Mendelian randomization (MR) analyses to evaluate whether genetically increased BMI is associated with an increased risk of MS.

METHODS AND FINDINGS

Employing a two-sample MR approach, we used summary statistics from the Genetic Investigation of Anthropometric Traits (GIANT) consortium and the International MS Genetics Consortium (IMSGC), the largest genome-wide association studies for BMI and MS, respectively (GIANT: n = 322,105; IMSGC: n = 14,498 cases and 24,091 controls). Seventy single nucleotide polymorphisms (SNPs) were genome-wide significant (p < 5 x 10-8) for BMI in GIANT (n = 322,105) and were investigated for their association with MS risk in the IMSGC. The effect of each SNP on MS was weighted by its effect on BMI, and estimates were pooled to provide a summary measure for the effect of increased BMI upon risk of MS. Our results suggest that increased BMI influences MS susceptibility, where a 1 standard deviation increase in genetically determined BMI (kg/m2) increased odds of MS by 41% (odds ratio [OR]: 1.41, 95% CI 1.20-1.66, p = 2.7 x 10-5, I2 = 0%, 95% CI 0-29). Sensitivity analyses, including MR-Egger regression, and the weighted median approach provided no evidence of pleiotropic effects. The main study limitations are that, while these sensitivity analyses reduce the possibility that pleiotropy influenced our results, residual pleiotropy is difficult to exclude entirely.

CONCLUSION

Genetically elevated BMI is associated with risk of MS, providing evidence for a causal role for obesity in MS etiology. While obesity has been associated with many late-life outcomes, these findings suggest an important consequence of childhood and/or early adulthood obesity.

The Genetics of Success: How Single-Nucleotide Polymorphisms Associated With Educational Attainment Relate to Life-Course Development

A previous genome-wide association study (GWAS) of more than 100,000 individuals identified molecular-genetic predictors of educational attainment. We undertook in-depth life-course investigation of the polygenic score derived from this GWAS using the four-decade Dunedin Study (N = 918). There were five main findings. First, polygenic scores predicted adult economic outcomes even after accounting for educational attainments. Second, genes and environments were correlated: Children with higher polygenic scores were born into better-off homes. Third, children's polygenic scores predicted their adult outcomes even when analyses accounted for their social-class origins; social-mobility analysis showed that children with higher polygenic scores were more upwardly mobile than children with lower scores. Fourth, polygenic scores predicted behavior across the life course, from early acquisition of speech and reading skills through geographic mobility and mate choice and on to financial planning for retirement. Fifth, polygenic-score associations were mediated by psychological characteristics, including intelligence, self-control, and interpersonal skill. Effect sizes were small. Factors connecting DNA sequence with life outcomes may provide targets for interventions to promote population-wide positive development.

Assortative mating and differential fertility by phenotype and genotype across the 20th century

This study asks two related questions about the shifting landscape of marriage and reproduction in US society over the course of the last century with respect to a range of health and behavioral phenotypes and their associated genetic architecture: (i) Has assortment on measured genetic factors influencing reproductive and social fitness traits changed over the course of the 20th century? (ii) Has the genetic covariance between fitness (as measured by total fertility) and other traits changed over time? The answers to these questions inform our understanding of how the genetic landscape of American society has changed over the past century and have implications for population trends. We show that husbands and wives carry similar loadings for genetic factors related to education and height. However, the magnitude of this similarity is modest and has been fairly consistent over the course of the 20th century. This consistency is particularly notable in the case of education, for which phenotypic similarity among spouses has increased in recent years. Likewise, changing patterns of the number of children ever born by phenotype are not matched by shifts in genotype-fertility relationships over time. Taken together, these trends provide no evidence that social sorting is becoming increasingly genetic in nature or that dysgenic dynamics have accelerated.

Robust analysis of secondary phenotypes in case-control genetic association studies

The case-control study is a common design for assessing the association between genetic exposures and a disease phenotype. Though association with a given (case-control) phenotype is always of primary interest, there is often considerable interest in assessing relationships between genetic exposures and other (secondary) phenotypes. However, the case-control sample represents a biased sample from the general population. As a result, if this sampling framework is not correctly taken into account, analyses estimating the effect of exposures on secondary phenotypes can be biased leading to incorrect inference. In this paper, we address this problem and propose a general approach for estimating and testing the population effect of a genetic variant on a secondary phenotype. Our approach is based on inverse probability weighted estimating equations, where the weights depend on genotype and the secondary phenotype. We show that, though slightly less efficient than a full likelihood-based analysis when the likelihood is correctly specified, it is substantially more robust to model misspecification, and can out-perform likelihood-based analysis, both in terms of validity and power, when the model is misspecified. We illustrate our approach with an application to a case-control study extracted from the Framingham Heart Study. Copyright © 2016 John Wiley & Sons, Ltd.

Risk of psychiatric illness from advanced paternal age is not predominantly from de novo mutations

The offspring of older fathers have higher risk of psychiatric disorders such as schizophrenia and autism. Paternal-age-related de novo mutations are widely assumed to be the underlying causal mechanism, and, although such mutations must logically make some contribution, there are alternative explanations (for example, elevated liability to psychiatric illness may delay fatherhood). We used population genetic models based on empirical observations of key parameters (for example, mutation rate, prevalence, and heritability) to assess the genetic relationship between paternal age and risk of psychiatric illness. These models suggest that age-related mutations are unlikely to explain much of the increased risk of psychiatric disorders in children of older fathers. Conversely, a model incorporating a weak correlation between age at first child and liability to psychiatric illness matched epidemiological observations. Our results suggest that genetic risk factors shared by older fathers and their offspring are a credible alternative explanation to de novo mutations for risk to children of older fathers.

FINEMAP: efficient variable selection using summary data from genome-wide association studies

MOTIVATION

The goal of fine-mapping in genomic regions associated with complex diseases and traits is to identify causal variants that point to molecular mechanisms behind the associations. Recent fine-mapping methods using summary data from genome-wide association studies rely on exhaustive search through all possible causal configurations, which is computationally expensive.

RESULTS

We introduce FINEMAP, a software package to efficiently explore a set of the most important causal configurations of the region via a shotgun stochastic search algorithm. We show that FINEMAP produces accurate results in a fraction of processing time of existing approaches and is therefore a promising tool for analyzing growing amounts of data produced in genome-wide association studies and emerging sequencing projects.

AVAILABILITY AND IMPLEMENTATION

FINEMAP v1.0 is freely available for Mac OS X and Linux at http://www.christianbenner.com

CONTACT

: christian.benner@helsinki.fi or matti.pirinen@helsinki.fi.

Genome-wide association study of serum coenzyme Q10 levels identifies susceptibility loci linked to neuronal diseases

Coenzyme Q₁₀ (CoQ₁₀) is a lipophilic redox molecule that is present in membranes of almost all cells in human tissues. CoQ₁₀ is, amongst other functions, essential for the respiratory transport chain and is a modulator of inflammatory processes and gene expression. Rare monogenetic CoQ₁₀ deficiencies show noticeable symptoms in tissues (e.g. kidney) and cell types (e.g. neurons) with a high energy demand. To identify common genetic variants influencing serum CoQ₁₀ levels, we performed a fixed effects meta-analysis in two independent cross-sectional Northern German cohorts comprising 1300 individuals in total. We identified two genome-wide significant susceptibility loci. The best associated single nucleotide polymorphism (SNP) was rs9952641 (P value = 1.31 × 10 ^-8, β = 0.063, CI_0.95 [0.041, 0.085]) within the COLEC12 gene on chromosome 18. The SNP rs933585 within the NRXN-1 gene on chromosome 2 also showed genome wide significance (P value = 3.64 × 10 ^-8, β = -0.034, CI_0.95 [-0.046, -0.022]). Both genes have been previously linked to neuronal diseases like Alzheimer's disease, autism and schizophrenia. Among our 'top-10' associated variants, four additional loci with known neuronal connections showed suggestive associations with CoQ₁₀ levels. In summary, this study demonstrates that serum CoQ₁₀ levels are associated with common genetic loci that are linked to neuronal diseases.