Power Analysis for Genetic Association Test (PAGEANT) provides insights to challenges for rare variant association studies

Ravages: An R package for the simulation and analysis of rare variants in multicategory phenotypes

Current software packages for the analysis and the simulations of rare variants are only available for binary and continuous traits. Ravages provides solutions in a single R package to perform rare variant association tests for multicategory, binary and continuous phenotypes, to simulate datasets under different scenarios and to compute statistical power. Association tests can be run in the whole genome thanks to C++ implementation of most of the functions, using either RAVA-FIRST, a recently developed strategy to filter and analyse genome-wide rare variants, or user-defined candidate regions. Ravages also includes a simulation module that generates genetic data for cases who can be stratified into several subgroups and for controls. Through comparisons with existing programmes, we show that Ravages complements existing tools and will be useful to study the genetic architecture of complex diseases. Ravages is available on the CRAN at https://cran.r-project.org/web/packages/Ravages/ and maintained on Github at https://github.com/genostats/Ravages.

Rare genetic variants explain missing heritability in smoking

Common genetic variants explain less variation in complex phenotypes than inferred from family-based studies, and there is a debate on the source of this 'missing heritability'. We investigated the contribution of rare genetic variants to tobacco use with whole-genome sequences from up to 26,257 unrelated individuals of European ancestries and 11,743 individuals of African ancestries. Across four smoking traits, single-nucleotide-polymorphism-based heritability ([Formula: see text]) was estimated from 0.13 to 0.28 (s.e., 0.10-0.13) in European ancestries, with 35-74% of it attributable to rare variants with minor allele frequencies between 0.01% and 1%. These heritability estimates are 1.5-4 times higher than past estimates based on common variants alone and accounted for 60% to 100% of our pedigree-based estimates of narrow-sense heritability ([Formula: see text], 0.18-0.34). In the African ancestry samples, [Formula: see text] was estimated from 0.03 to 0.33 (s.e., 0.09-0.14) across the four smoking traits. These results suggest that rare variants are important contributors to the heritability of smoking.

Whole Exome Sequencing Study Identifies Novel Rare Risk Variants for Habitual Coffee Consumption Involved in Olfactory Receptor and Hyperphagia

Habitual coffee consumption is an addictive behavior with unknown genetic variations and has raised public health issues about its potential health-related outcomes. We performed exome-wide association studies to identify rare risk variants contributing to habitual coffee consumption utilizing the newly released UK Biobank exome dataset (n = 200,643). A total of 34,761 qualifying variants were imported into SKAT to conduct gene-based burden and robust tests with minor allele frequency <0.01, adjusting the polygenic risk scores (PRS) of coffee intake to exclude the effect of common coffee-related polygenic risk. The gene-based burden and robust test of the exonic variants found seven exome-wide significant associations, such as OR2G2 (PSKAT = 1.88 × 10−9, PSKAT-Robust = 2.91 × 10−17), VEZT1 (PSKAT = 3.72 × 10−7, PSKAT-Robust = 1.41 × 10−7), and IRGC (PSKAT = 2.92 × 10−5, PSKAT-Robust = 1.07 × 10−7). These candidate genes were verified in the GWAS summary data of coffee intake, such as rs12737801 (p = 0.002) in OR2G2, and rs34439296 (p = 0.008) in IRGC. This study could help to extend genetic insights into the pathogenesis of coffee addiction, and may point to molecular mechanisms underlying health effects of habitual coffee consumption.

The impact of rare germline variants on human somatic mutation processes

Somatic mutations are an inevitable component of ageing and the most important cause of cancer. The rates and types of somatic mutation vary across individuals, but relatively few inherited influences on mutation processes are known. We perform a gene-based rare variant association study with diverse mutational processes, using human cancer genomes from over 11,000 individuals of European ancestry. By combining burden and variance tests, we identify 207 associations involving 15 somatic mutational phenotypes and 42 genes that replicated in an independent data set at a false discovery rate of 1%. We associate rare inherited deleterious variants in genes such as MSH3, EXO1, SETD2, and MTOR with two phenotypically different forms of DNA mismatch repair deficiency, and variants in genes such as EXO1, PAXIP1, RIF1, and WRN with deficiency in homologous recombination repair. In addition, we identify associations with other mutational processes, such as APEX1 with APOBEC-signature mutagenesis. Many of the genes interact with each other and with known mutator genes within cellular sub-networks. Considered collectively, damaging variants in the identified genes are prevalent in the population. We suggest that rare germline variation in diverse genes commonly impacts mutational processes in somatic cells.

Association of FAAH p.Pro129Thr and COMT p.Ala72Ser with schizophrenia and comorbid substance use through next-generation sequencing: an exploratory analysis

OBJECTIVE

Individuals with schizophrenia and substance use disorders have a poor prognosis and increased psychiatric symptoms. The present study aimed to explore the association of 106 genes in individuals with schizophrenia and comorbid substance use through a next-generation sequencing (NGS) analysis and different in silico algorithms.

METHODS

We included 105 individuals diagnosed with schizophrenia and a family history of schizophrenia, of whom 49 (46.67%) presented comorbid substance use. Using NGS, we sequenced 106 genes previously associated with schizophrenia. Logistic regression models were used to assess differences in allele frequencies, and a generalized gene-set analysis was performed at the gene level. Functional annotations were performed using different algorithms and databases.

RESULTS

We identified a total of 3,109 variants, of which 25 were associated with schizophrenia and comorbid substance use and were located in regulatory and coding regions. We found low-frequency variants in COMT p.Ala72Ser, independently of p.Val158Met, that were associated with substance use. The endocannabinoid functional variant FAAH p.Pro129Thr was also associated with substance use.

CONCLUSIONS

Genetic variants of genes related to dopaminergic and cannabinoid neurotransmitter systems were associated with comorbid substance use in schizophrenia. Nevertheless, more studies with larger sample sizes are needed to confirm our findings.

Identifying therapeutic drug targets using bidirectional effect genes

Prioritizing genes for translation to therapeutics for common diseases has been challenging. Here, we propose an approach to identify drug targets with high probability of success by focusing on genes with both gain of function (GoF) and loss of function (LoF) mutations associated with opposing effects on phenotype (Bidirectional Effect Selected Targets, BEST). We find 98 BEST genes for a variety of indications. Drugs targeting those genes are 3.8-fold more likely to be approved than non-BEST genes. We focus on five genes (IGF1R, NPPC, NPR2, FGFR3, and SHOX) with evidence for bidirectional effects on stature. Rare protein-altering variants in those genes result in significantly increased risk for idiopathic short stature (ISS) (OR = 2.75, p = 3.99 × 10-8). Finally, using functional experiments, we demonstrate that adding an exogenous CNP analog (encoded by NPPC) rescues the phenotype, thus validating its potential as a therapeutic treatment for ISS. Our results show the value of looking for bidirectional effects to identify and validate drug targets.

Exome Sequencing Identifies Abnormalities in Glycosylation and ANKRD36C in Patients with Immune-Mediated Thrombotic Thrombocytopenic Purpura

BACKGROUND

 Immune-mediated thrombotic thrombocytopenic purpura (iTTP) is a potentially fatal blood disorder, resulting from autoantibodies against ADAMTS13 (a disintegrin and metalloproteinase with a thrombospondin type 1 motif, member 13). However, the mechanism underlying anti-ADAMTS13 autoantibody formation is not known, nor it is known how genetic aberrations contribute to the pathogenesis of iTTP.

METHODS

 Here we performed whole exome sequencing (WES) of DNA samples from 40 adult patients with iTTP and 15 local healthy subjects with no history of iTTP and other hematological disorders.

RESULTS

 WES revealed variations in the genes involved in protein glycosylation, including O-linked glycosylation, to be a major pathway affected in patients with iTTP. Moreover, variations in the ANKRD gene family, particularly ANKRD36C and its paralogs, were also more prevalent in patients with iTTP than in the healthy controls. The ANKRD36 family of proteins have been implicated in inflammation. Mass spectrometry revealed a dramatic alternation in plasma glycoprotein profile in patients with iTTP compared with the healthy controls.

CONCLUSION

 Altered glycosylation may affect the disease onset and progression in various ways: it may predispose patients to produce ADAMTS13 autoantibodies or affect their binding properties; it may also alter clearance kinetics of hemostatic and inflammatory proteins. Together, our findings provide novel insights into plausible mechanisms underlying the pathogenesis of iTTP.

Dynamic incorporation of multiple in silico functional annotations empowers rare variant association analysis of large whole-genome sequencing studies at scale

Large-scale whole-genome sequencing studies have enabled the analysis of rare variants (RVs) associated with complex phenotypes. Commonly used RV association tests have limited scope to leverage variant functions. We propose STAAR (variant-set test for association using annotation information), a scalable and powerful RV association test method that effectively incorporates both variant categories and multiple complementary annotations using a dynamic weighting scheme. For the latter, we introduce 'annotation principal components', multidimensional summaries of in silico variant annotations. STAAR accounts for population structure and relatedness and is scalable for analyzing very large cohort and biobank whole-genome sequencing studies of continuous and dichotomous traits. We applied STAAR to identify RVs associated with four lipid traits in 12,316 discovery and 17,822 replication samples from the Trans-Omics for Precision Medicine Program. We discovered and replicated new RV associations, including disruptive missense RVs of NPC1L1 and an intergenic region near APOC1P1 associated with low-density lipoprotein cholesterol.

The contribution of rare genetic variants to the pathogenesis of polycystic ovary syndrome

Polycystic ovary syndrome (PCOS) is a highly heritable disorder, but only a small proportion of the heritability can be accounted for by common genetic risk variants identified to date. It is possible that variants with lower allele frequencies that cannot be detected using genome-wide association study arrays contribute to PCOS. Here, we discuss the challenges inherent to studying rare genetic variants in complex disease and review several recent studies that have used DNA sequencing techniques to investigate whether rare variants play a role in PCOS pathogenesis. We evaluate these findings in the context of the latest literature in PCOS and complex disease genetics.

U-PASS: unified power analysis and forensics for qualitative traits in genetic association studies

Summary

Despite the availability of existing calculators for statistical power analysis in genetic association studies, there has not been a model-invariant and test-independent tool that allows for both planning of prospective studies and systematic review of reported findings. In this work, we develop a web-based application U-PASS (Unified Power analysis of ASsociation Studies), implementing a unified framework for the analysis of common association tests for binary qualitative traits. The application quantifies the shared asymptotic power limits of the common association tests, and visualizes the fundamental statistical trade-off between risk allele frequency and odds ratio. The application also addresses the applicability of asymptotics-based power calculations in finite samples, and provides guidelines for single-SNP-based association tests. In addition to designing prospective studies, U-PASS enables researchers to retrospectively assess the statistical validity of previously reported associations.

Availability and implementation

U-PASS is an open-source R Shiny application. A live instance is hosted at https://power.stat.lsa.umich.edu. Source is available on https://github.com/Pill-GZ/U-PASS.

Supplementary information

Supplementary data are available at Bioinformatics online.

Pulling the covers in electronic health records for an association study with self-reported sleep behaviors

The electronic health record (EHR) contains rich histories of clinical care, but has not traditionally been mined for information related to sleep habits. Here, we performed a retrospective EHR study based on a cohort of 3,652 individuals with self-reported sleep behaviors documented from visits to the sleep clinic. These individuals were obese (mean body mass index 33.6 kg/m²) and had a high prevalence of sleep apnea (60.5%), however we found sleep behaviors largely concordant with prior prospective cohort studies. In our cohort, average wake time was 1 hour later and average sleep duration was 40 minutes longer on weekends than on weekdays (p < 10^-12). Sleep duration varied considerably as a function of age and tended to be longer in females and in whites. Additionally, through phenome-wide association analyses, we found an association of long weekend sleep with depression, and an unexpectedly large number of associations of long weekday sleep with mental health and neurological disorders (q < 0.05). We then sought to replicate previously published genetic associations with morning/evening preference on a subset of our cohort with extant genotyping data (n = 555). While those findings did not replicate in our cohort, a polymorphism (rs3754214) in high linkage disequilibrium with a previously published polymorphism near TARS2 was associated with long sleep duration (p < 0.01). Collectively, our results highlight the potential of the EHR for uncovering the correlates of human sleep in real-world populations.