Genome-wide rare variant analysis for thousands of phenotypes in over 70,000 exomes from two cohorts.
Ontology highlight
ABSTRACT:
SUBMITTER: Cirulli ET
PROVIDER: S-EPMC6987107 | biostudies-literature | 2020 Jan
REPOSITORIES: biostudies-literature
Publications
Genome-wide rare variant analysis for thousands of phenotypes in over 70,000 exomes from two cohorts.
Cirulli Elizabeth T ET White Simon S Read Robert W RW Elhanan Gai G Metcalf William J WJ Tanudjaja Francisco F Fath Donna M DM Sandoval Efren E Isaksson Magnus M Schlauch Karen A KA Grzymski Joseph J JJ Lu James T JT Washington Nicole L NL
Nature communications 20200128 1
Understanding the impact of rare variants is essential to understanding human health. We analyze rare (MAF < 0.1%) variants against 4264 phenotypes in 49,960 exome-sequenced individuals from the UK Biobank and 1934 phenotypes (1821 overlapping with UK Biobank) in 21,866 members of the Healthy Nevada Project (HNP) cohort who underwent Exome + sequencing at Helix. After using our rare-variant-tailored methodology to reduce test statistic inflation, we identify 64 statistically significant gene-bas ...[more]
Similar Datasets
Project description:Performance of a recently developed test for association between multivariate phenotypes and sets of genetic variants (MURAT) is demonstrated using measures of bone mineral density (BMD). By combining individual-level whole genome sequenced data from the UK10K study, and imputed genome-wide genetic data on individuals from the Study of Osteoporotic Fractures (SOF) and the Osteoporotic Fractures in Men Study (MrOS), a data set of 8810 individuals was assembled; tests of association were performed between autosomal gene-sets of genetic variants and BMD measured at lumbar spine and femoral neck. Distributions of p-values obtained from analyses of a single BMD phenotype are compared to those from the multivariate tests, across several region definitions and variant weightings. There is evidence of increased power with the multivariate test, although no new loci for BMD were identified. Among 17 genes highlighted either because there were significant p-values in region-based association tests or because they were in well-known BMD genes, 4 windows in 2 genes as well as 6 single SNPs in one of these genes showed association at genome-wide significant thresholds with the multivariate phenotype test but not with the single-phenotype test, Sequence Kernel Association Test (SKAT).
| S-EPMC5760616 | biostudies-literature
Project description:Genome-wide association studies have uncovered thousands of common variants associated with human disease, but the contribution of rare variants to common disease remains relatively unexplored. The UK Biobank contains detailed phenotypic data linked to medical records for approximately 500,000 participants, offering an unprecedented opportunity to evaluate the effect of rare variation on a broad collection of traits1,2. Here we study the relationships between rare protein-coding variants and 17,361 binary and 1,419 quantitative phenotypes using exome sequencing data from 269,171 UK Biobank participants of European ancestry. Gene-based collapsing analyses revealed 1,703 statistically significant gene-phenotype associations for binary traits, with a median odds ratio of 12.4. Furthermore, 83% of these associations were undetectable via single-variant association tests, emphasizing the power of gene-based collapsing analysis in the setting of high allelic heterogeneity. Gene-phenotype associations were also significantly enriched for loss-of-function-mediated traits and approved drug targets. Finally, we performed ancestry-specific and pan-ancestry collapsing analyses using exome sequencing data from 11,933 UK Biobank participants of African, East Asian or South Asian ancestry. Our results highlight a significant contribution of rare variants to common disease. Summary statistics are publicly available through an interactive portal ( http://azphewas.com/ ).
| S-EPMC8458098 | biostudies-literature
Project description:The recount2 resource is composed of over 70,000 uniformly processed human RNA-seq samples spanning TCGA and SRA, including GTEx. The processed data can be accessed via the recount2 website and the recount Bioconductor package. This workflow explains in detail how to use the recount package and how to integrate it with other Bioconductor packages for several analyses that can be carried out with the recount2 resource. In particular, we describe how the coverage count matrices were computed in recount2 as well as different ways of obtaining public metadata, which can facilitate downstream analyses. Step-by-step directions show how to do a gene-level differential expression analysis, visualize base-level genome coverage data, and perform an analyses at multiple feature levels. This workflow thus provides further information to understand the data in recount2 and a compendium of R code to use the data.
| S-EPMC5621122 | biostudies-literature
Project description:The utility of genotype imputation in genome-wide association studies is increasing as progressively larger reference panels are improved and expanded through whole-genome sequencing. Developing general guidelines for optimally cost-effective imputation, however, requires evaluation of performance issues that include the relative utility of study-specific compared with general/multipopulation reference panels; genotyping with various array scaffolds; effects of different ethnic backgrounds; and assessment of ranges of allele frequencies. Here we compared the effectiveness of study-specific reference panels to the commonly used 1000 Genomes Project (1000G) reference panels in the isolated Sardinian population and in cohorts of European ancestry including samples from Minnesota (USA). We also examined different combinations of genome-wide and custom arrays for baseline genotypes. In Sardinians, the study-specific reference panel provided better coverage and genotype imputation accuracy than the 1000G panels and other large European panels. In fact, even gene-centered custom arrays (interrogating ~200?000 variants) provided highly informative content across the entire genome. Gain in accuracy was also observed for Minnesotans using the study-specific reference panel, although the increase was smaller than in Sardinians, especially for rare variants. Notably, a combined panel including both study-specific and 1000G reference panels improved imputation accuracy only in the Minnesota sample, and only at rare sites. Finally, we found that when imputation is performed with a study-specific reference panel, cutoffs different from the standard thresholds of MACH-Rsq and IMPUTE-INFO metrics should be used to efficiently filter badly imputed rare variants. This study thus provides general guidelines for researchers planning large-scale genetic studies.
| S-EPMC4463504 | biostudies-literature
Project description:IntroductionEssential tremor (ET) is one of the most common movement disorders. Despite its high prevalence and heritability, its genetic etiology remains elusive with only a few susceptibility genes identified and poorly replicated. Our aim was to find novel candidate genes involved in ET predisposition through whole exome sequencing.MethodsWe studied eight multigenerational families (N = 40 individuals) with an autosomal-dominant inheritance using a comprehensive strategy combining whole exome sequencing followed by case-control association testing of prioritized variants in a separate cohort comprising 521 ET cases and 596 controls. We further performed gene-based burden analyses in an additional dataset comprising 789 ET patients and 770 healthy individuals to investigate whether there was an enrichment of rare deleterious variants within our candidate genes.ResultsFifteen variants co-segregated with disease status in at least one of the families, among which rs749875462 in CCDC183, rs535864157 in MMP10 and rs114285050 in GPR151 showed a nominal association with ET. However, we found no significant enrichment of rare variants within these genes in cases compared with controls. Interestingly, MMP10 protein is involved in the inflammatory response to neuronal damage and has been previously associated with other neurological disorders.ConclusionsWe prioritized a set of promising genes, especially MMP10, for further genetic and functional studies in ET. Our study suggests that rare deleterious coding variants that markedly increase susceptibility to ET are likely to be found in many genes. Future studies are needed to replicate and further infer biological mechanisms and potential disease causality for our identified genes.
| S-EPMC7856267 | biostudies-literature
Project description:Parkinson disease (PD) is a complex neurodegenerative disorder influenced by both environmental and genetic factors. While genome wide association studies have identified several susceptibility loci, many causal variants and genes underlying these associations remain undetermined. Identifying these is essential in order to gain mechanistic insight and identify biological pathways that may be targeted therapeutically. We hypothesized that gene-based enrichment of rare mutations is likely to be found within susceptibility loci for PD and may help identify causal genes. Whole-exome sequencing data from two independent cohorts were analyzed in tandem and by meta-analysis and a third cohort genotyped using the NeuroX-array was used for replication analysis. We employed collapsing methods (burden and the sequence kernel association test) to detect gene-based enrichment of rare, protein-altering variation within established PD susceptibility loci. Our analyses showed trends for three genes (GALC, PARP9 and SEC23IP), but none of these survived multiple testing correction. Our findings provide no evidence of rare mutation enrichment in genes within PD-associated loci, in our datasets. While not excluding that rare mutations in these genes may influence the risk of idiopathic PD, our results suggest that, if such effects exist, much larger sequencing datasets will be required for their detection.
| S-EPMC7529297 | biostudies-literature
Project description:OBJECTIVE:To determine whether a monogenic basis explains sudden infant death syndrome (SIDS) using an exome-wide focus. STUDY DESIGN:A cohort of 427 unrelated cases of SIDS (257 male; average age?=?2.7?±?1.9 months) underwent whole-exome sequencing. Exome-wide rare variant analyses were carried out with 278 SIDS cases of European ancestry (173 male; average age?=?2.7?±?1.98 months) and 973 ethnic-matched controls based on 6 genetic models. Ingenuity Pathway Analysis also was performed. The cohort was collected in collaboration with coroners, medical examiners, and pathologists by St George's University of London, United Kingdom, and Mayo Clinic, Rochester, Minnesota. Whole-exome sequencing was performed at the Genomic Laboratory, Kings College London, United Kingdom, or Mayo Clinic's Medical Genome Facility, Rochester, Minnesota. RESULTS:Although no exome-wide significant (P?<?2.5?×?10-6) difference in burden of ultra-rare variants was detected for any gene, 405 genes had a greater prevalence (P?<?.05) of ultra-rare nonsynonymous variants among cases with 17 genes at P?<?.005. Some of these potentially overrepresented genes may represent biologically plausible novel candidate genes for a monogenic basis for a portion of patients with SIDS. The top canonical pathway identified was glucocorticoid biosynthesis (P?=?.01). CONCLUSIONS:The lack of exome-wide significant genetic associations indicates an extreme heterogeneity of etiologies underlying SIDS. Our approach to understanding the genetic mechanisms of SIDS has far reaching implications for the SIDS research community as a whole and may catalyze new evidence-based SIDS research across multiple disciplines. Perturbations in glucocorticoid biosynthesis may represent a novel SIDS-associated biological pathway for future SIDS investigative research.
| S-EPMC6394853 | biostudies-literature
Project description:The mitochondrial genome is present at variable copy number between individuals. Mitochondria are vulnerable to oxidative stress, and their dysfunction may be associated with cardiovascular disease. The association of mitochondrial DNA copy number with cardiometabolic risk factors (lipids, glycaemic traits, inflammatory markers, anthropometry and blood pressure) was assessed in two independent cohorts of European origin women, one in whom outcomes were measured at mean (SD) age 30 (4.3) years (N=2278) and the second at 69.4 (5.5) years (N=2872). Mitochondrial DNA copy number was assayed by quantitative polymerase chain reaction. Associations were adjusted for smoking, sociodemographic status, laboratory factors and white cell traits. Out of a total of 12 outcomes assessed in both cohorts, mitochondrial DNA copy number showed little or no association with the majority (point estimates were close to zero and nearly all p-values were >0.01). The strongest evidence was for an inverse association in the older cohort with insulin (standardised beta [95%CI]: -0.06, [-0.098, -0.022], p=0.002), but this association did not replicate in the younger cohort. Our findings do not provide support for variation in mitochondrial DNA copy number having an important impact on cardio-metabolic risk factors in European origin women.
| S-EPMC5832987 | biostudies-literature
Project description:In genetic association analysis, a joint test of multiple distinct phenotypes can increase power to identify sets of trait-associated variants within genes or regions of interest. Existing multiphenotype tests for rare variants make specific assumptions about the patterns of association with underlying causal variants, and the violation of these assumptions can reduce power to detect association. Here, we develop a general framework for testing pleiotropic effects of rare variants on multiple continuous phenotypes using multivariate kernel regression (Multi-SKAT). Multi-SKAT models affect sizes of variants on the phenotypes through a kernel matrix and perform a variance component test of association. We show that many existing tests are equivalent to specific choices of kernel matrices with the Multi-SKAT framework. To increase power of detecting association across tests with different kernel matrices, we developed a fast and accurate approximation of the significance of the minimum observed P value across tests. To account for related individuals, our framework uses random effects for the kinship matrix. Using simulated data and amino acid and exome-array data from the METabolic Syndrome In Men (METSIM) study, we show that Multi-SKAT can improve power over single-phenotype SKAT-O test and existing multiple-phenotype tests, while maintaining Type I error rate.
| S-EPMC6330125 | biostudies-literature
Project description:OBJECTIVE:IgG4-related disease (IgG4-RD) is a heterogeneous, multiorgan condition of unclear aetiology that can cause organ failure. Difficulty recognising IgG4-RD contributes to diagnostic delays. We sought to identify key IgG4-RD phenotypes. METHODS:We used two cross-sectional studies assembled by an international, multispecialty network of IgG4-RD specialists who submitted 765 cases to derive and replicate phenotypic groups. Phenotype groups of disease manifestations and key covariate distributions across the identified groups were measured using latent class analysis. RESULTS:In the derivation cohort (n=493), we identified four groups with distinct manifestations: Group 1 (31%), Pancreato-Hepato-Biliary disease; Group 2 (24%), Retroperitoneal Fibrosis and/or Aortitis; Group 3 (24%), Head and Neck-Limited disease and Group 4 (22%), classic Mikulicz syndrome with systemic involvement. We replicated the identification of four phenotype groups in the replication cohort. Compared with cases in Groups 1, 2 and 4, respectively, cases in Group 3 were more likely to be female (OR 11.60 (95% CI 5.39 to 24.98), 10.35 (95% CI 4.63 to 23.15) and 9.24 (95% CI 3.53 to 24.20)) and Asian (OR 6.68 (95% CI 2.82 to 15.79), 7.43 (95% CI 2.97 to 18.56) and 6.27 (95% CI 2.27 to 17.29)). Cases in Group 4 had a higher median serum IgG4 concentration (1170?mg/dL) compared with groups 1-3 (316, 178 and 445?mg/dL, respectively, p<0.001). CONCLUSION:We identified four distinctive IgG4-RD phenotypes according to organ involvement. Being Asian or female may predispose individuals to head and neck-limited disease. These phenotypes serve as a framework for identifying IgG4-RD and studying its aetiology and optimal treatment.
| S-EPMC6996288 | biostudies-literature