Search

Genome-wide association analyses (GWAS) have identified more than 300 single nucleotide polymorphisms at 163 independent loci associated with coronary artery disease (CAD). However, there is no full understanding about the causal genes for CAD and the mechanisms of their action. We aimed to perform a post GWAS analysis to identify genes whose polymorphism may influence the risk of CAD. Using the UK Biobank GWAS summary statistics, we performed a gene-based association analysis. We found 63 genes significantly associated with CAD due to their within-gene polymorphisms. Many of these genes are well known. Some known CAD genes such as FURIN and SORT1 did not show the gene-based association because their variants had low GWAS signals or gene-based association was inflated by the strong GWAS signal outside the gene. For several known CAD genes, we demonstrated that their effects could be explained not only or not at all by their own variants but by the variants within the neighboring genes controlling their expression. Using a number of bioinformatics techniques, we suggested potential mechanisms underlying gene-CAD associations. Three genes, CDK19, NCALD and ARHGEF12 were not previously associated with CAD. The role of these genes should be clarified in further studies.

Back pain is the leading cause of years lived with disability worldwide, yet surprisingly, little is known regarding the biology underlying this condition. The impact of genetics is known for chronic back pain: its heritability is estimated to be at least 40%. Large genome-wide association studies have shown that common variation may account for up to 35% of chronic back pain heritability; rare variants may explain a portion of the heritability not explained by common variants. In this study, we performed the first gene-based association analysis of chronic back pain using UK Biobank imputed data including rare variants with moderate imputation quality. We discovered 2 genes, SOX5 and PANX3, influencing chronic back pain. The SOX5 gene is a well-known back pain gene. The PANX3 gene has not previously been described as having a role in chronic back pain. We showed that the association of PANX3 with chronic back pain is driven by rare noncoding intronic polymorphisms. This result was replicated in an independent sample from UK Biobank and validated using a similar phenotype, dorsalgia, from FinnGen Biobank. We also found that the PANX3 gene is associated with intervertebral disk disorders. We can speculate that a possible mechanism of action of PANX3 on back pain is due to its effect on the intervertebral disks.

To date, more than 100 loci associated with coronary artery disease (CAD) have been detected in large-scale genome-wide studies. For some of the several hundreds of genes located in these loci, roles in the pathogenesis of the disease have been shown. However, the genetic mechanisms and specific genes controlling this disease are still not fully understood. This study is aimed atin silicosearch for new CAD genes. We performed a gene-based association analysis, where all polymorphic variants within a gene are analyzed simultaneously. The analysis was based on the results of the genome-wide association studies (GWAS) available from the open databases MICAD (120,575 people, 85,112 markers) and UK Biobank (337,199 people, 10,894,597 markers). We used the sumFREGAT package implementing a wide range of new methods for gene-based association analysis using summary statistics. We found 88 genes demonstrating significant gene-based associations. Forty-four of the identified genes were already known as CAD genes. Furthermore, we identified 28 additional genes in the known CAD loci. They can be considered as new candidate genes. Finally, we identified sixteen new genes (AGPAT4, ARHGEF12, BDP1, DHX58, EHBP1, FBF1, HSPB9, NPBWR2, PDLIM5, PLCB3, PLEKHM2, POU2F3, PRKD2, TMEM136, TTC29andUTP20) outside the known loci. Information about the functional role of these genes allows us to consider many of them as candidates for CAD. The 41 identified genes did not have significant GWAS signals and they were identified only due to simultaneous consideration of all variants within the gene in the framework of gene-based analysis. These results demonstrate that gene-based association analysis is a powerful tool for gene mapping. The method can utilize huge amounts of GWAS results accumulated in the world to map different traits and diseases. This type of studies is widely available, as it does not require additional material costs.

Gene-based association analysis is an effective gene-mapping tool. Many gene-based methods have been proposed recently. However, their power depends on the underlying genetic architecture, which is rarely known in complex traits, and so it is likely that a combination of such methods could serve as a universal approach. Several frameworks combining different gene-based methods have been developed. However, they all imply a fixed set of methods, weights and functional annotations. Moreover, most of them use individual phenotypes and genotypes as input data. Here, we introduce sumSTAAR, a framework for gene-based association analysis using summary statistics obtained from genome-wide association studies (GWAS). It is an extended and modified version of STAAR framework proposed by Li and colleagues in 2020. The sumSTAAR framework offers a wider range of gene-based methods to combine. It allows the user to arbitrarily define a set of these methods, weighting functions and probabilities of genetic variants being causal. The methods used in the framework were adapted to analyse genes with large number of SNPs to decrease the running time. The framework includes the polygene pruning procedure to guard against the influence of the strong GWAS signals outside the gene. We also present new improved matrices of correlations between the genotypes of variants within genes. These matrices estimated on a sample of 265,000 individuals are a state-of-the-art replacement of widely used matrices based on the 1000 Genomes Project data.

Motivation A huge number of genome-wide association studies (GWAS) summary statistics freely available in databases provide a new material for gene-based association analysis aimed at identifying rare genetic variants. Only a few of the many popular gene-based methods developed for individual genotype and phenotype data are adapted for the practical use of the GWAS summary statistics as input. Results We analytically prove and numerically illustrate that all popular powerful methods developed for gene-based association analysis of individual phenotype and genotype data can be modified to utilize GWAS summary statistics. We have modified and implemented all of the popular methods, including burden and kernel machine-based tests, multiple and functional linear regression, principal components analysis and others, in the R package sumFREGAT. Using real summary statistics for coronary artery disease, we show that the new package is able to detect genes not found by the existing packages. Availability and implementation The R package sumFREGAT is freely and publicly available at: https://CRAN.R-project.org/package=sumFREGAT. Supplementary information Supplementary data are available at Bioinformatics online.

Regional association analysis is a new statistical method which simultaneously considers all variants in a selected genome region. This method was created for the analysis of rare genetic variants, whose genotypes are determined by exome or genome sequencing. The gene is usually considered as a region. It was also proposed to use a regional analysis for testing of the association between a complex trait and a set of common variants genotyped by the panels developed for genome-wide association analysis. In this case, overlapping genome regions (sliding windows) are usually considered as a region. Since the size of such regions can be rather large, there is a risk of overestimation (inflation) of the test statistic and an increase in the type I error. In this work, the effect of the size of the region on the type I error was studied for traits with different heritability. The results of simulating experiments demonstrated that the physical size of the region but not the number of genetic variants in it is a limiting factor. The higher the trait heritability, the greater the type I error differs from the declared value. The analysis of a large number of real traits confirmed these conclusions. It is necessary to take into account these results during the interpretation of the results of regional association analysis conducted on large regions using common genetic variants.

Background Despite high heritability, little success was achieved in mapping genetic determinants of depression-related traits by means of genome-wide association studies. Methods To identify genes associated with depressive symptomology, we performed a gene-based association analysis of nonsynonymous variation captured using exome-sequencing and exome-chip genotyping in a genetically isolated population from the Netherlands ( n = 1999). Finally, we reproduced our significant findings in an independent population-based cohort ( n = 1604). Results We detected significant association of depressive symptoms with a gene NKPD1 ( p = 3.7 × 10 −08 ). Nonsynonymous variants in the gene explained 0.9% of sex- and age-adjusted variance of depressive symptoms in the discovery study, which is translated into 3.8% of the total estimated heritability ( h 2 = 0.24). Significant association of depressive symptoms with NKPD1 was also observed ( n = 1604; p = 1.5 × 10 −03 ) in the independent replication sample despite little overlap with the discovery cohort in the set of nonsynonymous genetic variants observed in the NKPD1 gene. Meta-analysis of the discovery and replication studies improved the association signal ( p = 1.0 × 10 −09 ). Conclusions Our study suggests that nonsynonymous variation in the gene NKPD1 affects depressive symptoms in the general population. NKPD1 is predicted to be involved in the de novo synthesis of sphingolipids, which have been implicated in the pathogenesis of depression.

Background Electrocardiographic measures of left ventricular hypertrophy (LVH) are used as predictors of cardiovascular risk. We combined linkage and association analyses to discover novel rare genetic variants involved in three such measures and two principal components derived from them. Methods The study was conducted among participants from the Erasmus Rucphen Family Study (ERF), a Dutch family-based sample from the southwestern Netherlands. Variance components linkage analyses were performed using Merlin. Regions of interest (LOD > 1.9) were fine-mapped using microarray and exome sequence data. Results We observed one significant LOD score for the second principal component on chromosome 15 (LOD score = 3.01) and 12 suggestive LOD scores. Several loci contained variants identified in GWAS for these traits; however, these did not explain the linkage peaks, nor did other common variants. Exome sequence data identified two associated variants after multiple testing corrections were applied. Conclusions We did not find common SNPs explaining these linkage signals. Exome sequencing uncovered a relatively rare variant in MAPK3K11 on chromosome 11 (MAF = 0.01) that helped account for the suggestive linkage peak observed for the first principal component. Conditional analysis revealed a drop in LOD from 2.01 to 0.88 for MAP3K11, suggesting that this variant may partially explain the linkage signal at this chromosomal location. MAP3K11 is related to the JNK pathway and is a pro-apoptotic kinase that plays an important role in the induction of cardiomyocyte apoptosis in various pathologies, including LVH. Electronic supplementary material The online version of this article (10.1186/s12920-018-0339-9) contains supplementary material, which is available to authorized users.

Summary We propose an automatic heuristic algorithm for splitting large pedigrees into fragments of no more than a user-specified bit size. The algorithm specifically aims to split large pedigrees where many close relatives are genotyped and to produce a set of sub-pedigrees for haplotype reconstruction, IBD computation or multipoint linkage analysis with the help of the Lander-Green-Kruglyak algorithm. We demonstrate that a set of overlapping pedigree fragments constructed with the help of our algorithm allows fast and effective haplotype reconstruction and detection of an allele’s parental origin. Moreover, we compared pedigree fragments constructed with the help of our algorithm and existing programs PedCut and Jenti for multipoint linkage analysis. Our algorithm demonstrated significantly higher linkage power than the algorithm of Jenti and significantly shorter running time than the algorithm of PedCut. The software package PedStr implementing our algorithms is available at http://mga.bionet.nsc.ru/soft/index.html.

For pedigrees with multiple loops, exact likelihoods could not be computed in an acceptable time frame and thus, approximate methods are used. Some of these methods are based on breaking loops and approximations of complex pedigree likelihoods using the exact likelihood of the corresponding zero-loop pedigree. Due to ignoring loops, this method results in a loss of genetic information and a decrease in the power to detect linkage. To minimize this loss, an optimal set of loop breakers has to be selected. In this paper, we present a graph theory based algorithm for automatic selection of an optimal set of loop breakers. We propose using a total relationship between measured pedigree members as a proxy to power. To minimize the loss of genetic information, we suggest selection of such breakers whose duplication in a pedigree would be accompanied by a minimal loss of total relationship between measured pedigree members. We show that our algorithm compares favorably with other existing loop-breaker selection algorithms in terms of conservation of genetic information, statistical power and CPU time of subsequent linkage analysis. We implemented our method in a software package LOOP_EDGE, which is available at http://mga.bionet.nsc.ru/nlru/.

Variation of serum protein allotypes serving as genetic markers of the blood has been analyzed in 29 populations of the domestic pig and subspecies of the wild boar. The population biodiversity and genetic structure have been estimated by two methods: by the frequencies of allotype combinations and with the use of a map constructed in the space of two principal components. The results obtained are the basis for determining the characteristics of the microevolution of wild boars and formation of the breeds of domestic pigs.

A method for reconstructing allele frequencies characteristic of an original ethnically homogeneous population before the start of migration processes is described. Information on both the ethnic group studied and offspring of interethnic marriages is used to estimate the allele frequencies. This makes it possible to increase the informativeness of the sample, which, in the case of ethnic heterogeneity, depends not only on allele frequencies and the total sample size, but also on the ethnic structure of the sample. The problem of estimating allele frequency in an ethnically heterogeneous sample has been solved analytically for diallelic loci. It has been demonstrated that, if offspring of interethnic marriages with the same degree of outbreeding is added to a sample of the ethnic group studied, the sample informativeness does not change. To utilize the information contained in the phenotypes of the offspring of interethnic marriages, representatives of the population from which migration occurs should be included into the sample. The size of the sample ensuring the preassigned accuracy of estimation is minimized at a certain ratio between the numbers of the offspring of interethnic marriages and the “immigrants.” To analyze polyallelic loci, a software package has been developed that allows estimating allele frequencies, determining the errors of these estimates, and planning the sample ensuring the preassigned accuracy of estimation. The package is available free at http://mga.bionet.nsc.ru/PopMixed/PopMixed.html.

Procedure is described to estimate allele frequencies in indigenous populations of Siberia using phenotype data not only for “pure-blood” representatives of the ethnic groups examined, but also for the descendants of mixed marriages. Implementation of the method requires reconstruction of the pedigree structure for the sample examined. Inclusion of the data on descendants of mixed marriages into the analysis increases the sample information content and decreases variance of the estimates obtained. The advantages of the method are illustrated using an example of Tundra Nentsy, for whom it was shown that variance of estimates at the analysis of the blood groups allele frequencies can be diminished approximately by a factor of 1.5.

Personality traits are complex phenotypes related to psychosomatic health. Individually, various gene finding methods have not achieved much success in finding genetic variants associated with personality traits. We performed a meta-analysis of four genome-wide linkage scans (N=6149 subjects) of five basic personality traits assessed with the NEO Five-Factor Inventory. We compared the significant regions from the meta-analysis of linkage scans with the results of a meta-analysis of genome-wide association studies (GWAS) (N∼17 000). We found significant evidence of linkage of neuroticism to chromosome 3p14 (rs1490265, LOD=4.67) and to chromosome 19q13 (rs628604, LOD=3.55); of extraversion to 14q32 (ATGG002, LOD=3.3); and of agreeableness to 3p25 (rs709160, LOD=3.67) and to two adjacent regions on chromosome 15, including 15q13 (rs970408, LOD=4.07) and 15q14 (rs1055356, LOD=3.52) in the individual scans. In the meta-analysis, we found strong evidence of linkage of extraversion to 4q34, 9q34, 10q24 and 11q22, openness to 2p25, 3q26, 9p21, 11q24, 15q26 and 19q13 and agreeableness to 4q34 and 19p13. Significant evidence of association in the GWAS was detected between openness and rs677035 at 11q24 (P-value=2.6 × 10 -06, KCNJ1). The findings of our linkage meta-analysis and those of the GWAS suggest that 11q24 is a susceptible locus for openness, with KCNJ1 as the possible candidate gene. © 2013 Macmillan Publishers Limited. All rights reserved.

The NEO-Five-Factor Inventory divides human personality traits into five dimensions: neuroticism, extraversion, openness, conscientiousness and agreeableness. In this study, we sought to identify regions harboring genes with large effects on the five NEO personality traits by performing genome-wide linkage analysis of individuals scoring in the extremes of these traits (>90th percentile). Affected-only linkage analysis was performed using an Illumina 6K linkage array in a family-based study, the Erasmus Rucphen Family study. We subsequently determined whether distinct, segregating haplotypes found with linkage analysis were associated with the trait of interest in the population. Finally, a dense single-nucleotide polymorphism genotyping array (Illumina 318K) was used to search for copy number variations (CNVs) in the associated regions. In the families with extreme phenotype scores, we found significant evidence of linkage for conscientiousness to 20p13 (rs1434789, log of odds (LOD) = 5.86) and suggestive evidence of linkage (LOD >2.8) for neuroticism to 19q, 21q and 22q, extraversion to 1p, 1q, 9p and 12q, openness to 12q and 19q, and agreeableness to 2p, 6q, 17q and 21q. Further analysis determined haplotypes in 21q22 for neuroticism (P-values = 0.009, 0.007), in 17q24 for agreeableness (marginal P-value = 0.018) and in 20p13 for conscientiousness (marginal P-values = 0.058, 0.038) segregating in families with large contributions to the LOD scores. No evidence for CNVs in any of the associated regions was found. Our findings imply that there may be genes with relatively large effects involved in personality traits, which may be identified with next-generation sequencing techniques. Molecular Psychiatry (2012) 17, 1031-1041; doi:10.1038/mp.2011.97; published online 9 August 2011

Phospho- and sphingolipids are crucial cellular and intracellular compounds. These lipids are required for active transport, a number of enzymatic processes, membrane formation, and cell signalling. Disruption of their metabolism leads to several diseases, with diverse neurological, psychiatric, and metabolic consequences. A large number of phospholipid and sphingolipid species can be detected and measured in human plasma. We conducted a meta-analysis of five European family-based genome-wide association studies (N = 4034) on plasma levels of 24 sphingomyelins (SPM), 9 ceramides (CER), 57 phosphatidylcholines (PC), 20 lysophosphatidylcholines (LPC), 27 phosphatidylethanolamines (PE), and 16 PE-based plasmalogens (PLPE), as well as their proportions in each major class. This effort yielded 25 genome-wide significant loci for phospholipids (smallest P-value = 9.88×10−204) and 10 loci for sphingolipids (smallest P-value = 3.10×10−57). After a correction for multiple comparisons (P-value, Author Summary Phospho- and sphingolipids are integral to membrane formation and are involved in crucial cellular functions such as signalling, membrane fluidity, membrane protein trafficking, neurotransmission, and receptor trafficking. In addition to severe monogenic diseases resulting from defective phospho- and sphingolipid function and metabolism, the evidence suggests that variations in these lipid levels at the population level are involved in the determination of cardiovascular and neurologic traits and subsequent disease. We took advantage of modern laboratory methods, including microarray-based genotyping and electrospray ionization tandem mass spectrometry, to hunt for genetic variation influencing the levels of more than 350 phospho- and sphingolipid phenotypes. We identified nine novel loci, in addition to confirming a number of previously described loci. Several other genetic regions provided substantial evidence of their involvement in these traits. All of these loci are strong candidates for further research in the field of lipid biology and are likely to yield considerable insights into the complex metabolic pathways underlying circulating phospho- and sphingolipid levels. Understanding these mechanisms might help to illuminate factors leading to the development of common cardiovascular and neurological diseases and might provide molecular targets for the development of new therapies.

Variation in the manifestation age is typical of many mitochondrial diseases. The estimation of penetrance of pathogenic mutations causing such diseases is usually conducted on samples of individuals whose age exceeds the maximum age of the disease manifestation. In the case of rare diseases, samples of sufficient size sometimes cannot be formed. In this study, we propose a method for estimating penetrance involving individuals of any age. The efficiency of the method is demonstrated using Leber hereditary optic neuropathy as an example. It is shown that the method provides an unbiased estimate of penetrance and considerably reduces the error of this estimate in comparison with a sample including individuals whose age exceeds the maximum age of disease manifestation.

Despite extensive research of genetic determinants of human adult height, the genes identified up until now allow to predict only a small proportion of the trait's variance. To identify new genes we analyzed 2,486 genotyped and phenotyped individuals in a large pedigree including 23,612 members in 18 generations. The pedigree was derived from a young genetically isolated Dutch population, where genetic heterogeneity is expected to be low and linkage disequilibrium has been shown to be increased. Complex segregation analysis confirmed high heritability of adult height, and suggested mixed model of height inheritance in this population. The estimates of the model parameters obtained from complex segregation analysis were used in parametric linkage analysis, which highlighted three genome-wide significant and additionally at least four suggestive loci involved in height. Significant peaks were located at the chromosomal regions 1p32 (LOD score = 3.35), 2p16 (LOD score = 3.29) and 16q24 (LOD score = 3.94). For the latter region, a strong association signal (FDR q < 0.05) was obtained for 19 SNPs, 17 of them were located in the CDH13 (cadherin 13) gene of which one (rs1035569) explained 1.5% of the total height variance.

Genes for height have gained interest for decades, but only recently have candidate genes started to be identified. We have performed linkage analysis and genome-wide association for height in approximately 4000 individuals from five European populations. A total of five chromosomal regions showed suggestive linkage and in one of these regions, two SNPs (rs849140 and rs1635852) were associated with height (nominal P = 7.0 x 10(-8) and P = 9.6 x 10(-7), respectively). In total, five SNPs across the genome showed an association with height that reached the threshold of genome-wide significance (nominal P < 1.6 x 10(-7)). The association with height was replicated for two SNPs (rs1635852 and rs849140) using three independent studies (n = 31 077, n=1268 and n = 5746) with overall meta P-values of 9.4 x 10(-10) and 5.3 x 10(-8). These SNPs are located in the JAZF1 gene, which has recently been associated with type II diabetes, prostate and endometrial cancer. JAZF1 is a transcriptional repressor of NR2C2, which results in low IGF1 serum concentrations, perinatal and early postnatal hypoglycemia and growth retardation when knocked out in mice. Both the linkage and association analyses independently identified the JAZF1 region affecting human height. We have demonstrated, through replication in additional independent populations, the consistency of the effect of the JAZF1 SNPs on height. Since this gene also has a key function in the metabolism of growth, JAZF1 represents one of the strongest candidates influencing human height identified so far.

Sphingolipids have essential roles as structural components of cell membranes and in cell signalling, and disruption of their metabolism causes several diseases, with diverse neurological, psychiatric, and metabolic consequences. Increasingly, variants within a few of the genes that encode enzymes involved in sphingolipid metabolism are being associated with complex disease phenotypes. Direct experimental evidence supports a role of specific sphingolipid species in several common complex chronic disease processes including atherosclerotic plaque formation, myocardial infarction (MI), cardiomyopathy, pancreatic β-cell failure, insulin resistance, and type 2 diabetes mellitus. Therefore, sphingolipids represent novel and important intermediate phenotypes for genetic analysis, yet little is known about the major genetic variants that influence their circulating levels in the general population. We performed a genome-wide association study (GWAS) between 318,237 single-nucleotide polymorphisms (SNPs) and levels of circulating sphingomyelin (SM), dihydrosphingomyelin (Dih-SM), ceramide (Cer), and glucosylceramide (GluCer) single lipid species (33 traits); and 43 matched metabolite ratios measured in 4,400 subjects from five diverse European populations. Associated variants (32) in five genomic regions were identified with genome-wide significant corrected p-values ranging down to 9.08×10−66. The strongest associations were observed in or near 7 genes functionally involved in ceramide biosynthesis and trafficking: SPTLC3, LASS4, SGPP1, ATP10D, and FADS1–3. Variants in 3 loci (ATP10D, FADS3, and SPTLC3) associate with MI in a series of three German MI studies. An additional 70 variants across 23 candidate genes involved in sphingolipid-metabolizing pathways also demonstrate association (p = 10−4 or less). Circulating concentrations of several key components in sphingolipid metabolism are thus under strong genetic control, and variants in these loci can be tested for a role in the development of common cardiovascular, metabolic, neurological, and psychiatric diseases., Author Summary Although several rare monogenic diseases are caused by defects in enzymes involved in sphingolipid biosynthesis and metabolism, little is known about the major variants that control the circulating levels of these important bioactive molecules. As well as being essential components of plasma membranes and endosomes, sphingolipids play critical roles in cell surface protection, protein and lipid transport and sorting, and cellular signalling cascades. Experimental evidence supports a role for sphingolipids in several common complex chronic metabolic, cardiovascular, or neurological disease processes. Therefore, sphingolipids represent novel and important intermediate phenotypes for genetic analysis, and discovering the genetic variants that influence their circulating concentrations is an important step towards understanding how the genetic control of sphingolipids might contribute to common human disease. We have identified 32 variants in 7 genes that have a strong effect on the circulating plasma levels of 33 distinct sphingolipids, and 43 matched metabolite ratios. In a series of 3 German MI studies, we see association with MI for variants in 3 of the genes tested. Further cardiovascular, metabolic, neurological, and psychiatric disease associations can be tested with the variants described here, which may identify additional disease risk and potentially useful therapeutic targets.