Search

Epidemiological studies report beneficial associations of higher educational attainment (EDU) with Alzheimer’s disease (AD). Prior genome-wide association studies (GWAS) also reported variants associated with AD and EDU separately. The analysis of pleiotropic associations with these phenotypes may shed light on EDU-related protection against AD. We performed pleiotropic meta-analyses using Fisher’s method and omnibus test applied to summary statistics for single nucleotide polymorphisms (SNPs) associated with AD and EDU in large-scale univariate GWAS at suggestive-effect (5 × 10−8 < p < 0.1) and genome-wide (p ≤ 5 × 10−8) significance levels. We report 53 SNPs that attained p ≤ 5 × 10−8 at least in one of the pleiotropic meta-analyses and were reported in the univariate GWAS at 5 × 10−8 < p < 0.1. Of them, there were 46 pleiotropic SNPs according to Fisher’s method. Additionally, Fisher’s method identified 25 of 206 SNPs with pleiotropic effects, which attained p ≤ 5 × 10−8 in the univariate GWAS. We showed that a large fraction of the pleiotropic associations was affected by a counterintuitive phenomenon of antagonistic genetic heterogeneity, which explains the increase, rather than decrease, of the significance of the pleiotropic associations in the omnibus test. Functional enrichment analysis showed that apart from cancers, gene set harboring the non-pleiotropic SNPs was characterized by late-onset AD and neurodevelopmental disorders. The pleiotropic gene set was characterized by a broad spectrum of progressive neurological and neuromuscular diseases and immune-mediated conditions, including progressive motor neuropathy, multiple sclerosis, Parkinson’s disease, and severe AD. Our results suggest that disentangling genes harboring variants with and without pleiotropic associations with AD and EDU is promising for dissecting heterogeneity in biological mechanisms of AD. Supplementary Information The online version contains supplementary material available at 10.1007/s11357-021-00484-1.

Summary The proportion of variation in complex traits that can be attributed to non-additive genetic effects has been a topic of intense debate. The availability of biobank-scale datasets of genotype and trait data from unrelated individuals opens up the possibility of obtaining precise estimates of the contribution of non-additive genetic effects. We present an efficient method to estimate the variation in a complex trait that can be attributed to additive (additive heritability) and dominance deviation (dominance heritability) effects across all genotyped SNPs in a large collection of unrelated individuals. Over a wide range of genetic architectures, our method yields unbiased estimates of additive and dominance heritability. We applied our method, in turn, to array genotypes as well as imputed genotypes (at common SNPs with minor allele frequency [MAF] > 1%) and 50 quantitative traits measured in 291,273 unrelated white British individuals in the UK Biobank. Averaged across these 50 traits, we find that additive heritability on array SNPs is 21.86% while dominance heritability is 0.13% (about 0.48% of the additive heritability) with qualitatively similar results for imputed genotypes. We find no statistically significant evidence for dominance heritability (p

Complex-trait genetics has advanced dramatically through methods to estimate the heritability tagged by SNPs, both genome-wide and in genomic regions of interest such as those defined by functional annotations. The models underlying many of these analyses are inadequate, and consequently many SNP-heritability results published to date are inaccurate. Here, we review the modelling issues, both for analyses based on individual genotype data and association test statistics, highlighting the role of a low-dimensional model for the heritability of each SNP. We use state-of-art models to present updated results about how heritability is distributed with respect to functional annotations in the human genome, and how it varies with allele frequency, which can reflect purifying selection. Our results give finer detail to the picture that has emerged in recent years of complex trait heritability widely dispersed across the genome. Confounding due to population structure remains a problem that summary statistic analyses cannot reliably overcome. Complex-trait genetics has advanced dramatically through methods to estimate the heritability tagged by SNPs, both genome-wide and in genomic regions of interest such as those defined by functional annotations. The models underlying many of these analyses are inadequate, and consequently many SNP-heritability results published to date are inaccurate. Here, we review the modelling issues, both for analyses based on individual genotype data and association test statistics, highlighting the role of a low-dimensional model for the heritability of each SNP. We use state-of-art models to present updated results about how heritability is distributed with respect to functional annotations in the human genome, and how it varies with allele frequency, which can reflect purifying selection. Our results give finer detail to the picture that has emerged in recent years of complex trait heritability widely dispersed across the genome. Confounding due to population structure remains a problem that summary statistic analyses cannot reliably overcome. Also see the video abstract here: https://youtu.be/WC2u03V65MQ.

Hypothesis-free Mendelian randomization studies provide a way to assess the causal relevance of a trait across the human phenome but can be limited by statistical power, sample overlap or complicated by horizontal pleiotropy. The recently described latent causal variable (LCV) approach provides an alternative method for causal inference which might be useful in hypothesis-free experiments across human phenome. We developed an automated pipeline for phenome-wide tests using the LCV approach including steps to estimate partial genetic causality, filter to a meaningful set of estimates, apply correction for multiple testing and then present the findings in a graphical summary termed causal architecture plot. We apply this pipeline to body mass index (BMI) and lipid traits as exemplars of traits where there is strong prior expectation for causal effects, and to dental caries and periodontitis as exemplars of traits where there is a need for causal inference. The results for lipids and BMI suggest that these traits are best viewed as contributing factors on a multitude of traits and conditions, thus providing additional evidence that supports viewing these traits as targets for interventions to improve health. On the other hand, caries and periodontitis are best viewed as a downstream consequence of other traits and diseases rather than a cause of ill health. The automated pipeline is implemented in the Complex-Traits Genetics Virtual Lab ( https://vl.genoma.io ) and results are available in https://view.genoma.io . We propose causal architecture plots based on phenome-wide partial genetic causality estimates as a new way visualizing the overall causal map of the human phenome.

Summary To facilitate scientific collaboration on polygenic risk scores (PRSs) research, we created an extensive PRS online repository for 35 common cancer traits integrating freely available genome-wide association studies (GWASs) summary statistics from three sources: published GWASs, the NHGRI-EBI GWAS Catalog, and UK Biobank-based GWASs. Our framework condenses these summary statistics into PRSs using various approaches such as linkage disequilibrium pruning/p value thresholding (fixed or data-adaptively optimized thresholds) and penalized, genome-wide effect size weighting. We evaluated the PRSs in two biobanks: the Michigan Genomics Initiative (MGI), a longitudinal biorepository effort at Michigan Medicine, and the population-based UK Biobank (UKB). For each PRS construct, we provide measures on predictive performance and discrimination. Besides PRS evaluation, the Cancer-PRSweb platform features construct downloads and phenome-wide PRS association study results (PRS-PheWAS) for predictive PRSs. We expect this integrated platform to accelerate PRS-related cancer research.

BackgroundAdolescent idiopathic scoliosis (AIS) is a genetically heterogeneous disease characterised by three-dimensional deformity of the spine in the absence of a congenital spinal anomaly or neurological musculoskeletal disorder. The clinical variability and incomplete penetrance of some genes linked with AIS indicate that this disease constitutes an oligogenic trait.ObjectiveWe aimed to explore the oligogenic nature of this disease and identify novel AIS genes.MethodsWe analysed rare damaging variants within AIS-associated genes by using exome sequencing in 40 AIS trios and 183 sporadic patients.ResultsMultiple variants within AIS-associated genes were identified in eight AIS trios, and five individuals harboured rare damaging variants in theFLNBgene. The patients showed more frequent oligogenicity than the controls. In the gene-based burden test, the top signal resided inFLNB. In functional studies, we found that the AIS-associatedFLNBvariants altered the protein’s conformation and subcellular localisation and its interaction with other proteins (TTC26 and OFD1) involved in AIS. The most compelling evidence of an oligogenic basis was that the number of rare damaging variants was recognised as an independent prognostic factor for curve progression in Cox regression analysis.ConclusionOur data indicate that AIS is an oligogenic disease and identifyFLNBas a susceptibility gene for AIS.

Mouse models have been instrumental in understanding human disease biology and proposing possible new treatments. The precise control of the environment and genetic composition of mice allows more rigorous observations, but limits the generalizability and translatability of the results into human applications. In the era of precision medicine, strategies using mouse models have to be revisited to effectively emulate human populations. Systems genetics is one promising paradigm that may promote the transition to novel precision medicine strategies. Here, we review the state-of-the-art resources and discuss how mouse systems genetics helps to understand human diseases and to advance the development of precision medicine, with an emphasis on the existing resources and strategies.

In transcriptome‐wide association studies (TWAS), gene expression values are predicted using genotype data and tested for association with a phenotype. The power of this approach to detect associations relies, at least in part, on the accuracy of the prediction. Here we compare the prediction accuracy of six different methods—LASSO, Ridge regression, Elastic net, Best Linear Unbiased Predictor, Bayesian Sparse Linear Mixed Model, and Random Forests—by performing cross‐validation using data from the Geuvadis Project. We also examine prediction accuracy (a) at different sample sizes, (b) when ancestry of the prediction model training and testing populations is different, and (c) when the tissue used to train the model is different from the tissue to be predicted. We find that, for most genes, the expression cannot be accurately predicted, but in general sparse statistical models tend to outperform polygenic models at prediction. Average prediction accuracy is reduced when the model training set size is reduced or when predicting across ancestries and is marginally reduced when predicting across tissues. We conclude that using sparse statistical models and the development of large reference panels across multiple ethnicities and tissues will lead to better prediction of gene expression, and thus may improve TWAS power.

BackgroundParkinson’s disease (PD) is a neurodegenerative disorder with complex genetic architecture. Besides rare mutations in high-risk genes related to monogenic familial forms of PD, multiple variants associated with sporadic PD were discovered via association studies.MethodsWe studied the whole-exome sequencing data of 340 PD cases and 146 ethnically matched controls from the Parkinson’s Progression Markers Initiative (PPMI) and performed burden analysis for different rare variant classes. Disease prediction models were built based on clinical, non-clinical and genetic features, including both common and rare variants, and two machine learning methods.ResultsWe observed a significant exome-wide burden of singleton loss-of-function variants (corrected p=0.037). Overall, no exome-wide burden of rare amino acid changing variants was detected. Finally, we built a disease prediction model combining singleton loss-of-function variants, a polygenic risk score based on common variants, and family history of PD as features and reached an area under the curve of 0.703 (95% CI 0.698 to 0.708). By incorporating a rare variant feature, our model increased the performance of the state-of-the-art classification model for the PPMI dataset, which reached an area under the curve of 0.639 based on common variants alone.ConclusionThe main finding of this study is to highlight the contribution of singleton loss-of-function variants to the complex genetics of PD and that disease risk prediction models combining singleton and common variants can improve models built solely on common variants.

To determine the impact of genetic variants on the brain, we used genetically informed brain atlases in genome-wide association studies of regional cortical surface area and thickness in 39,898 adults and 9136 children. We uncovered 440 genome-wide significant loci in the discovery cohort and 800 from a post hoc combined meta-analysis. Loci in adulthood were largely captured in childhood, showing signatures of negative selection, and were linked to early neurodevelopment and pathways associated with neuropsychiatric risk. Opposing gradations of decreased surface area and increased thickness were associated with common inversion polymorphisms. Inferior frontal regions, encompassing Broca’s area, which is important for speech, were enriched for human-specific genomic elements. Thus, a mixed genetic landscape of conserved and human-specific features is concordant with brain hierarchy and morphogenetic gradients.

The spatial distribution of genetic variants is jointly determined by geography, past demographic processes, natural selection, and its interplay with environmental variation. A fraction of these genetic variants are “causal alleles” that affect the manifestation of a complex trait. The effect exerted by these causal alleles on complex traits can be independent or dependent on the environment. Understanding the evolutionary processes that shape the spatial structure of causal alleles is key to comprehend the spatial distribution of complex traits. Natural selection, past population size changes, range expansions, consanguinity, assortative mating, archaic introgression, admixture, and the environment can alter the frequencies, effect sizes, and heterozygosities of causal alleles. This provides a genetic axis along which complex traits can vary. However, complex traits also vary along biogeographical and sociocultural axes which are often correlated with genetic axes in complex ways. The purpose of this review is to consider these genetic and environmental axes in concert and examine the ways they can help us decipher the variation in complex traits that is visible in humans today. This initiative necessarily implies a discussion of populations, traits, the ability to infer and interpret “genetic” components of complex traits, and how these have been impacted by adaptive events. In this review, we provide a history-aware discussion on these topics using both the recent and more distant past of our academic discipline and its relevant contexts.

Scientifically interesting as well as practically important phenotypes often belong to the realm of complex traits. To the extent that these traits are hereditary, they are usually ‘highly polygenic’. The study of such traits presents a challenge for researchers, as the complex genetic architecture of such traits makes it nearly impossible to utilise many of the usual methods of reverse genetics, which often focus on specific genes. In recent years, thousands of genome-wide association studies (GWAS) were undertaken to explore the relationships between complex traits and a large number of genetic factors, most of which are characterised by tiny effects. In this review, we aim to familiarise ‘wet biologists’ with approaches for the interpretation of GWAS results, to clarify some issues that may seem counterintuitive and to assess the possibility of using GWAS results in experiments on various complex traits.

The contemporary European genetic makeup formed in the last 8,000 years when local Western Hunter-Gatherers (WHGs) mixed with incoming Anatolian Neolithic farmers and Pontic Steppe pastoralists.1–3 This encounter combined genetic variants with distinct evolutionary histories and, together with new environmental challenges faced by the post-Neolithic Europeans, unlocked novel adaptations.4 Previous studies inferred phenotypes in these source populations, using either a few single loci5–7 or polygenic scores based on genome-wide association studies,8–10 and investigated the strength and timing of natural selection on lactase persistence or height, among others.6,11,12 However, how ancient populations contributed to present-day phenotypic variation is poorly understood. Here, we investigate how the unique tiling of genetic variants inherited from different ancestral components drives the complex traits landscape of contemporary Europeans and quantify selection patterns associated with these components. Using matching individual-level genotype and phenotype data for 27 traits in the Estonian biobank13 and genotype data directly from the ancient source populations, we quantify the contributions from each ancestry to present-day phenotypic variation in each complex trait. We find substantial differences in ancestry for eye and hair color, body mass index, waist/hip circumferences, and their ratio, height, cholesterol levels, caffeine intake, heart rate, and age at menarche. Furthermore, we find evidence for recent positive selection linked to four of these traits and, in addition, sleep patterns and blood pressure. Our results show that these ancient components were differentiated enough to contribute ancestry-specific signatures to the complex trait variability displayed by contemporary Europeans.

It is of great interest and potential to discover causal relationships between pairs of exposures and outcomes using genetic variants as instrumental variables (IVs) to deal with hidden confounding in observational studies. Two most popular approaches are Mendelian randomization (MR), which usually use independent genetic variants/SNPs across the genome, and transcriptome-wide association studies (TWAS) (or their generalizations) using cis-SNPs local to a gene (or some genome-wide and likely dependent SNPs), as IVs. In spite of their many promising applications, both approaches face a major challenge: the validity of their causal conclusions depends on three critical assumptions on valid IVs, and more generally on other modeling assumptions, which however may not hold in practice. The most likely as well as challenging situation is due to the wide-spread horizontal pleiotropy, leading to two of the three IV assumptions being violated and thus to biased statistical inference. More generally, we’d like to conduct a goodness-of-fit (GOF) test to check the model being used. Although some methods have been proposed as being robust to various degrees to the violation of some modeling assumptions, they often give different and even conflicting results due to their own modeling assumptions and possibly lower statistical efficiency, imposing difficulties to the practitioner in choosing and interpreting varying results across different methods. Hence, it would help to directly test whether any assumption is violated or not. In particular, there is a lack of such tests for TWAS. We propose a new and general GOF test, called TEDE (TEsting Direct Effects), applicable to both correlated and independent SNPs/IVs (as commonly used in TWAS and MR respectively). Through simulation studies and real data examples, we demonstrate high statistical power and advantages of our new method, while confirming the frequent violation of modeling (including valid IV) assumptions in practice and thus the importance of model checking by applying such a test in MR/TWAS analysis., Author summary With the increasing availability of large-scale GWAS summary data of various complex traits/diseases and software packages, it has become convenient and popular to apply Mendelian randomization (MR) and transcriptome-wide association studies (TWAS), using genetic variants as instrumental variables (IVs), to address fundamental and significant questions by unraveling causal relationships between complex or molecular traits such as gene expression and other complex traits. However, the validity of such causal conclusions critically depends on the validity of the model being used, including three key IV assumptions. In particular, with the wide-spread horizontal pleiotropy of genetic variants, two of the three IV assumptions may be violated, leading to biased inference from MR and TWAS. This issue may become more severe as more trait-associated genetic variants are used as IVs to increase the power of MR and TWAS. Although there are some methods to check the modeling assumptions for MR with independent genetic variants as IVs, there is barely any powerful one for TWAS (or more generally for MR and similar methods) with correlated SNPs as IVs. We propose such a powerful method applicable to both MR and TWAS with local or genome-wide, possibly correlated/dependent, SNPs as IVs, demonstrating its higher statistical power than several commonly used methods, while confirming the frequent violation of modeling/IV assumptions in TWAS with our example GWAS data of schizophrenia, Alzheimer’s disease and blood lipids. An important conclusion is that in practice it is necessary to conduct model checking in MR and TWAS, and our proposed method is expected to be useful for such a task.

Background Migraine is a complex neurological disorder that is considered the most common disabling brain disorder affecting 14 % of people worldwide. The present study sought to infer potential causal relationships between self-reported migraine and other complex traits, using genetic data and a hypothesis-free approach. Methods We leveraged available summary statistics from genome-wide association studies (GWAS) of 1,504 phenotypes and self-reported migraine and inferred pair-wise causal relationships using the latent causal variable (LCV) method. Results We identify 18 potential causal relationships between self-reported migraine and other complex traits. Hypertension and blood clot formations were causally associated with an increased migraine risk, possibly through vasoconstriction and platelet clumping. We observed that sources of abdominal pain and discomfort might influence a higher risk for migraine. Moreover, occupational and environmental factors such as working with paints, thinner or glues, and being exposed to diesel exhaust were causally associated with higher migraine risk. Psychiatric-related phenotypes, including stressful life events, increased migraine risk. In contrast, ever feeling unenthusiastic / disinterested for a whole week, a phenotype related to the psychological well-being of individuals, was a potential outcome of migraine. Conclusions Overall, our results suggest a potential vascular component to migraine, highlighting the role of vasoconstriction and platelet clumping. Stressful life events and occupational variables potentially influence a higher migraine risk. Additionally, a migraine could impact the psychological well-being of individuals. Our findings provide novel testable hypotheses for future studies that may inform the design of new interventions to prevent or reduce migraine risk and recurrence.

Genetic prediction of complex traits has great promise for disease prevention, monitoring, and treatment. The development of accurate risk prediction models is hindered by the wide diversity of genetic architecture across different traits, limited access to individual level data for training and parameter tuning, and the demand for computational resources. To overcome the limitations of the most existing methods that make explicit assumptions on the underlying genetic architecture and need a separate validation data set for parameter tuning, we develop a summary statistics-based nonparametric method that does not rely on validation datasets to tune parameters. In our implementation, we refine the commonly used likelihood assumption to deal with the discrepancy between summary statistics and external reference panel. We also leverage the block structure of the reference linkage disequilibrium matrix for implementation of a parallel algorithm. Through simulations and applications to twelve traits, we show that our method is adaptive to different genetic architectures, statistically robust, and computationally efficient. Our method is available at https://github.com/eldronzhou/SDPR., Author summary Recently there has been much interest in predicting an individual’s phenotype from genetic information, which has great promise for disease prevention, monitoring, and treatment. It has been found that there is great variation in the genetic architecture underlying different complex traits, including the number of genetic variants involved and the distribution of the effect sizes of genetic variants. How to model such genetic contribution is a key aspect for accurate prediction of complex traits. So far, most existing methods make specific assumptions about the shape of the genetic contribution. If these assumptions are not correct, the prediction accuracy might be compromised. Here we propose a method that learns the shape of the genetic contribution without making any explicit assumptions. We found that our method achieved robust performance when compared with other recently developed methods through simulation and real data analysis. Our method is also practically more feasible, since it supports the use of public summary statistics and consumes only small amount of computational resources.

The accuracy of polygenic risk scores (PRSs) to predict complex diseases increases with the training sample size. PRSs are generally derived based on summary statistics from large meta-analyses of multiple genome-wide association studies (GWAS). However, it is now common for researchers to have access to large individual-level data as well, such as the UK biobank data. To the best of our knowledge, it has not yet been explored how to best combine both types of data (summary statistics and individual-level data) to optimize polygenic prediction. The most widely used approach to combine data is the meta-analysis of GWAS summary statistics (Meta-GWAS), but we show that it does not always provide the most accurate PRS. Through simulations and using twelve real case-control and quantitative traits from both iPSYCH and UK Biobank along with external GWAS summary statistics, we compare Meta-GWAS with two alternative data-combining approaches, stacked clumping and thresholding (SCT) and Meta-PRS. We find that, when large individual-level data is available, the linear combination of PRSs (Meta-PRS) is both a simple alternative to Meta-GWAS and often more accurate.

[EN] Complex traits include, among many others, the evolution of eyes, wings, body forms, reproductive modes, human intelligence, social behaviour, diseases, and chromosome morphology. Dollo¿s law states that the evolution of complex traits is irreversible. However, potential exceptions have been proposed. Here, we investigated whether reticulation, a simple and elegant means by which complex characters may be reacquired, could account for suggested reversals in the evolution of complex characters using two datasets with sufficient genetic coverage and a total of five potential reversals. Our analyses uncovered a potential reversal in the evolution of parity mode and a potential reversal in the evolution of placentotrophy of fish (Cyprinodontiformes) as reticulation events. Moreover, in a reptile that exhibits a potential reversal from viviparity to oviparity (Zootoca vivipara), reticulation provided the most parsimonious explanation for sex chromosome evolution. Therefore, three of the five studied potential reversals were unravelled as reticulation events. This constitutes the first evidence that accounting for reticulation can fundamentally influence the interpretation of the evolution of complex traits, that testing for reticulation is crucial for obtaining robust phylogenies, and that complex ancestral characters may be reacquired through hybridization with a lineage that still exhibits the trait. Hybridization, rather than re-appearance of ancestral traits by means of small evolutionary steps, may thus account for suggested exceptions to Dollo¿s law. Consequently, ruling out reticulation is required to claim evolutionary reversal of complex characters and potential exceptions to Dollo's rule, EU Marie Curie‐Clarín CoFund, Grant/Award Number: ACA14–26; Swiss National Science Foundation, Grant/Award Numbers: PP00P3_152929/1, PPOOP3_128375; Spanish MINECO postdoc, Grant/Award Number: FPDI‐2013‐16116; Spanish Ministry of Education and Science, Grant/Award Number: CGL2008‐01522 and CGL2012‐32459

A role for hereditary influences in the susceptibility for chronic obstructive pulmonary disease (COPD) is widely recognized. Cytotoxic lymphocytes are implicated in COPD pathogenesis, and functions of these leukocytes are modulated by interactions between their killer-cell immunoglobulin-like receptors (KIR) and human leukocyte antigen (HLA)-Class I molecules on target cells. We hypothesized HLA-Class I and KIR inheritance affect risks for COPD. HLA-Class I alleles and KIR genotypes were defined by candidate gene analyses in multiple cohorts of COPD patients (total n=392) and control smokers with normal spirometry (total n=342). Compared to controls, COPD patients had over-representations of HLA-C*07 and activating KIR2DS1, with under-representations of HLA-C*12. Particular HLA-KIR permutations were synergistic; e.g. the presence of HLA-C*07 + KIR2DS1 + HLA-C12null vs. HLAC*07null + KIR2DS1null + HLA-C12 was associated with COPD, especially among HLA-C1 allotype homozygotes (OR=18.5, 95%CI=3.7-90.9, p

BackgroundHeight and other anthropometric measures are consistently found to associate with differential cancer risk. However, both genetic and mechanistic insights into these epidemiological associations are notably lacking. Conversely, inherited genetic variants in tumour suppressors and oncogenes increase cancer risk, but little is known about their influence on anthropometric traits.MethodsBy integrating inherited and somatic cancer genetic data from the Genome-Wide Association Study Catalog, expression Quantitative Trait Loci databases and the Cancer Gene Census, we identify SNPs that associate with different cancer types and differential gene expression in at least one tissue type, and explore the potential pleiotropic associations of these SNPs with anthropometric traits through SNP-wise association in a cohort of 500,000 individuals.ResultsWe identify three regulatory SNPs for three important cancer genes, FANCA, MAP3K1 and TP53 that associate with both anthropometric traits and cancer risk. Of particular interest, we identify a previously unrecognised strong association between the rs78378222[C] SNP in the 3' untranslated region (3'-UTR) of TP53 and both increased risk for developing non-melanomatous skin cancer (OR=1.36 (95% 1.31 to 1.41), adjusted p=7.62E−63), brain malignancy (OR=3.12 (2.22 to 4.37), adjusted p=1.43E−12) and increased standing height (adjusted p=2.18E−24, beta=0.073±0.007), lean body mass (adjusted p=8.34E−37, beta=0.073±0.005) and basal metabolic rate (adjusted p=1.13E−31, beta=0.076±0.006), thus offering a novel genetic link between these anthropometric traits and cancer risk.ConclusionOur results clearly demonstrate that heritable variants in key cancer genes can associate with both differential cancer risk and anthropometric traits in the general population, thereby lending support for a genetic basis for linking these human phenotypes.

ObjectiveTo determine the effects of HIV serostatus and disease severity on endothelial function in a large pooled cohort study of people living with HIV infection and HIV- controls. Approach and Results: We used participant-level data from 9 studies: 7 included people living with HIV (2 treatment-naïve) and 4 had HIV- controls. Brachial artery flow-mediated dilation (FMD) was measured using a standardized ultrasound imaging protocol with central reading. After data harmonization, multiple linear regression was used to examine the effects of HIV- serostatus, HIV disease severity measures, and cardiovascular disease risk factors on FMD. Of 2533 participants, 986 were people living with HIV (mean 44.4 [SD 11.8] years old) and 1547 were HIV- controls (42.9 [12.2] years old). The strongest and most consistent associates of FMD were brachial artery diameter, age, sex, and body mass index. The effect of HIV+ serostatus on FMD was strongly influenced by kidney function. In the highest tertile of creatinine (1.0 mg/dL), the effect of HIV+ serostatus was strong (β=-1.59% [95% CI, -2.58% to -0.60%], P=0.002), even after covariate adjustment (β=-1.36% [95% CI, -2.46% to -0.47%], P=0.003). In the lowest tertile (0.8 mg/dL), the effect of HIV+ serostatus was strong (β=-1.90% [95% CI, -2.58% to -1.21%], P

Although the cohort-level accuracy of polygenic risk scores (PRSs)-estimates of genetic value at the individual level-has been widely assessed, uncertainty in PRSs remains underexplored. In the present study, we show that Bayesian PRS methods can estimate the variance of an individual's PRS and can yield well-calibrated credible intervals via posterior sampling. For 13 real traits in the UK Biobank (n = 291,273 unrelated 'white British'), we observe large variances in individual PRS estimates which impact interpretation of PRS-based stratification; averaging across traits, only 0.8% (s.d. = 1.6%) of individuals with PRS point estimates in the top decile have corresponding 95% credible intervals fully contained in the top decile. We provide an analytical estimator for the expectation of individual PRS variance as a function of SNP heritability, number of causal SNPs and sample size. Our results showcase the importance of incorporating uncertainty in individual PRS estimates into subsequent analyses.

Significance Some of the genes responsible for the evolution of light skin pigmentation in Europeans show signals of positive selection in present-day populations. Recently, genome-wide association studies have highlighted the highly polygenic nature of skin pigmentation. It is unclear whether selection has operated on all of these genetic variants or just a subset. By studying variation in over a thousand ancient genomes from West Eurasia covering 40,000 y, we are able to study both the aggregate behavior of pigmentation-associated variants and the evolutionary history of individual variants. We find that the evolution of light skin pigmentation in Europeans was driven by frequency changes in a relatively small fraction of the genetic variants that are associated with variation in the trait today., Skin pigmentation is a classic example of a polygenic trait that has experienced directional selection in humans. Genome-wide association studies have identified well over a hundred pigmentation-associated loci, and genomic scans in present-day and ancient populations have identified selective sweeps for a small number of light pigmentation-associated alleles in Europeans. It is unclear whether selection has operated on all of the genetic variation associated with skin pigmentation as opposed to just a small number of large-effect variants. Here, we address this question using ancient DNA from 1,158 individuals from West Eurasia covering a period of 40,000 y combined with genome-wide association summary statistics from the UK Biobank. We find a robust signal of directional selection in ancient West Eurasians on 170 skin pigmentation-associated variants ascertained in the UK Biobank. However, we also show that this signal is driven by a limited number of large-effect variants. Consistent with this observation, we find that a polygenic selection test in present-day populations fails to detect selection with the full set of variants. Our data allow us to disentangle the effects of admixture and selection. Most notably, a large-effect variant at SLC24A5 was introduced to Western Europe by migrations of Neolithic farming populations but continued to be under selection post-admixture. This study shows that the response to selection for light skin pigmentation in West Eurasia was driven by a relatively small proportion of the variants that are associated with present-day phenotypic variation.

People in the Americas represent a diverse continuum of populations with varying degrees of admixture among African, European, and Amerindigenous ancestries. In the United States, populations with non-European ancestry remain understudied, and thus little is known about the genetic architecture of phenotypic variation in these populations. Using genotype data from the Hispanic Community Health Study/Study of Latinos, we find that Amerindigenous ancestry increased by an average of ~20% spanning 1940s-1990s in Mexican Americans. These patterns result from complex interactions between several population and cultural factors which shaped patterns of genetic variation and influenced the genetic architecture of complex traits in Mexican Americans. We show for height how polygenic risk scores based on summary statistics from a European-based genome-wide association study perform poorly in Mexican Americans. Our findings reveal temporal changes in population structure within Hispanics/Latinos that may influence biomedical traits, demonstrating a need to improve our understanding of admixed populations.

The use of genomic information to better understand and prevent common complex diseases has been an ongoing goal of genetic research. Over the past few years, research in this area has proliferated with several proposed methods of generating polygenic scores. This has been driven by the availability of larger data sets, primarily from genome-wide association studies and concomitant developments in statistical methodologies. Here we provide an overview of the methodological aspects of polygenic model construction. In addition, we consider the state of the field and implications for potential applications of polygenic scores for risk estimation within healthcare.

The selection pressures that have shaped the evolution of complex traits in humans remain largely unknown, and in some contexts highly contentious, perhaps above all where they concern mean trait differences among groups. To date, the discussion has focused on whether such group differences have any genetic basis, and if so, whether they are without fitness consequences and arose via random genetic drift, or whether they were driven by selection for different trait optima in different environments. Here, we highlight a plausible alternative: that many complex traits evolve under stabilizing selection in the face of shifting environmental effects. Under this scenario, there will be rapid evolution at the loci that contribute to trait variation, even when the trait optimum remains the same. These considerations underscore the strong assumptions about environmental effects that are required in ascribing trait differences among groups to genetic differences., This Essay discusses how complex traits may evolve under selection to maintain trait values in the face of changing environments; this scenario has implications for interpreting genetic differences between groups and signals of adaptation.

Phenotype prediction is a key goal for medical genetics. Unfortunately, most genome-wide association studies are done in European populations, which reduces the accuracy of predictions via polygenic scores in non-European populations. Here, we use population genetic models to show that human demographic history and negative selection on complex traits can result in population-specific genetic architectures. For traits where alleles with the largest effect on the trait are under the strongest negative selection, approximately half of the heritability can be accounted for by variants in Europe that are absent from Africa, leading to poor performance in phenotype prediction across these populations. Further, under such a model, individuals in the tails of the genetic risk distribution may not be identified via polygenic scores generated in another population. We empirically test these predictions by building a model to stratify heritability between European-specific and shared variants and applied it to 37 traits and diseases in the UK Biobank. Across these phenotypes, ∼30% of the heritability comes from European-specific variants. We conclude that genetic association studies need to include more diverse populations to enable the utility of phenotype prediction in all populations.

Despite more than two decades of extensive research focusing on nonalcoholic fatty liver disease (NAFLD), no approved therapy for steatohepatitis—the severe histological form of the disease—presently exists. More importantly, new drugs and small molecules with diverse molecular targets on the pathways of hepatocyte injury, inflammation, and fibrosis cannot achieve the primary efficacy endpoints. Precision medicine can potentially overcome this issue, as it is founded on extensive knowledge of the druggable genome/proteome. Hence, this review summarizes significant trends and developments in precision medicine with a particular focus on new potential therapeutic discoveries modeled via systems biology approaches. In addition, we computed and simulated the potential utility of the NAFLD polygenic risk score, which could be conceptually very advantageous not only for early disease detection but also for implementing actionable measures. Incomplete knowledge of the druggable NAFLD genome severely impedes the drug discovery process and limits the likelihood of identifying robust and safe drug candidates. Thus, we close this article with some insights into emerging disciplines, such as chemical genetics, that may accelerate accurate identification of the druggable NAFLD genome/proteome. Fil: Sookoian, Silvia Cristina. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Houssay. Instituto de Investigaciones Médicas. Universidad de Buenos Aires. Facultad de Medicina. Instituto de Investigaciones Médicas; Argentina Fil: Pirola, Carlos José. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Houssay. Instituto de Investigaciones Médicas. Universidad de Buenos Aires. Facultad de Medicina. Instituto de Investigaciones Médicas; Argentina

Genome-Wide Association Studies (GWAS) in large human cohorts have identified thousands of loci associated with complex traits and diseases. For identifying the genes and gene-associated variants that underlie complex traits in livestock, especially where sample sizes are limiting, it may help to integrate the results of GWAS for equivalent traits in humans as prior information. In this study, we sought to investigate the usefulness of results from a GWAS on human height as prior information for identifying the genes and gene-associated variants that affect stature in cattle, using GWAS summary data on samples sizes of 700,000 and 58,265 for humans and cattle, respectively. Using Fisher’s exact test, we observed a significant proportion of cattle stature-associated genes (30/77) that are also associated with human height (odds ratio = 5.1, p = 3.1e-10). Result of randomized sampling tests showed that cattle orthologs of human height-associated genes, hereafter referred to as candidate genes (C-genes), were more enriched for cattle stature GWAS signals than random samples of genes in the cattle genome (p = 0.01). Randomly sampled SNPs within the C-genes also tend to explain more genetic variance for cattle stature (up to 13.2%) than randomly sampled SNPs within random cattle genes (p = 0.09). The most significant SNPs from a cattle GWAS for stature within the C-genes did not explain more genetic variance for cattle stature than the most significant SNPs within random cattle genes (p = 0.87). Altogether, our findings support previous studies that suggest a similarity in the genetic regulation of height across mammalian species. However, with the availability of a powerful GWAS for stature that combined data from 8 cattle breeds, prior information from human-height GWAS does not seem to provide any additional benefit with respect to the identification of genes and gene-associated variants that affect stature in cattle., Author summary Provided that there is a similarity in the biological control of complex traits in different species, the results of genome-wide association study for a given trait in a model species can be utilized as prior information to identify genes and gene-associated variants that affect the same trait in a non-model species. Here we report that the number of genes found to be affecting height in both humans and cattle is more than can be expected by chance, suggesting a similarity in the biology of height in both species. In addition, genes that are associated with height in humans were found to be highly enriched with variants that are strongly associated with height in cattle. Using height as a model trait, we demonstrated that across-species prior information cannot compete with a within-species prior information, with respect to the identification of genes and the gene-associated variants that affect a complex trait. Across-species prior information can be useful, however, when within-species prior information is either unavailable or poor.

Many genetic mutations affecting phenotypes are presumed to do so via altering gene expression in particular cells or tissues, but identifying the specific genes involved has been challenging. A transcriptome-wide association study (TWAS) attempts to identify disease associated genes by first learning a predictive model on an eQTL dataset and then imputing gene expression levels into a larger genome-wide association study (GWAS). Finally, associations between predicted gene expressions and GWAS phenotype are identified.Here, we compared tree-based machine learning (ML) method of random forests (RF) with more widely used linear methods of lasso, ridge, and elastic net regression, for prediction of gene expression. We also developed a multi-task learning extension to RF which simultaneously makes use of information from multiple tissues (RF-MTL) and compared it to a multi-dataset version of lasso, the joint lasso, and to a single tissue RF. We found that for prediction of gene expression, RF, in general, outperformed linear approaches on our chosen eQTL dataset and that multi-tissue methods generally outperformed their single-tissue counterparts, with RF-MTL performing the best. Simulations showed that these benefits generally propagated to the next steps of the analysis, although highlighted that joint lasso had a tendency to erroneously identify genes in one tissue if there existed a disease signal for that gene in another.We tested all four methods on type 1 diabetes (T1D) GWAS and expression data for several immune cells and found that 46 genes were identified by at least one method, though only 7 by all methods. Joint lasso discovered the most T1D-associated genes, including 15 unique to that method, but this may reflect its higher false positive rate due to “overborrowing” information across tissues. RF-MTL found more unique associated genes than RF for 3 out 5 tissues. Compared to lasso-based analysis, the RF gene list was more likely to relate to T1D in an analysis of independent data types. We conclude that RF, both single- and multi-task version, is competitive and, for some cell types, superior to linear models conventionally used in the TWAS studies.Author summaryA transcriptome-wide association study (TWAS) is a way of integrating expression data and genome-wide association studies (GWAS), which allows for discovery of genes, rather than mutations, associated to traits of interest. In the TWAS framework, we first train predictive models on an eQTL dataset, then use these models to impute gene expression into a GWAS dataset. Finally, we look for significant associations between predicted gene expression and a GWAS trait. In this work, we compare non-linear method of random forests (RF) to linear models, customarily used in TWAS. Furthermore, we demonstrate that TWAS framework can naturally be extended to, and potentially benefit from, a multi-tissue setting, thereby taking advantage of the correlation between gene expression in different tissue types. We applied the RF, a selection of linear models, and the multi-tissue approaches to an eQTL dataset of monocytes and B cells and a large T1D GWAS. We found that RF outperform lasso in terms of predictive accuracy and the number of differentially expressed genes found, and that multi-dataset version of lasso discovered the most T1D-associated genes. Analysis of the gene lists produced for each method in independent data types (excluding genetic association data) showed all related to T1D, but that the RF methods ranked T1D higher in their lists than the linear methods. We conclude that RF is a useful addition to the TWAS tool box.