Testing for association with rare variants in the coding and non-coding genome: RAVA-FIRST, a new approach based on CADD deleteriousness score.

Rare variant association tests (RVAT) have been developed to study the contribution of rare variants widely accessible through high-throughput sequencing technologies. RVAT require to aggregate rare variants in testing units and to filter variants to retain only the most likely causal ones. In the exome, genes are natural testing units and variants are usually filtered based on their functional consequences. However, when dealing with whole-genome sequence (WGS) data, both steps are challenging. No natural biological unit is available for aggregating rare variants. Sliding windows procedures have been proposed to circumvent this difficulty, however they are blind to biological information and result in a large number of tests. We propose a new strategy to perform RVAT on WGS data: "RAVA-FIRST" (RAre Variant Association using Functionally-InfoRmed STeps) comprising three steps. (1) New testing units are defined genome-wide based on functionally-adjusted Combined Annotation Dependent Depletion (CADD) scores of variants observed in the gnomAD populations, which are referred to as "CADD regions". (2) A region-dependent filtering of rare variants is applied in each CADD region. (3) A functionally-informed burden test is performed with sub-scores computed for each genomic category within each CADD region. Both on simulations and real data, RAVA-FIRST was found to outperform other WGS-based RVAT. Applied to a WGS dataset of venous thromboembolism patients, we identified an intergenic region on chromosome 18 enriched for rare variants in early-onset patients. This region that was missed by standard sliding windows procedures is included in a TAD region that contains a strong candidate gene. RAVA-FIRST enables new investigations of rare non-coding variants in complex diseases, facilitated by its implementation in the R package Ravages.

Polygenic risk of major depressive disorder as a risk factor for venous thromboembolism.

Major depressive disorder (MDD), bipolar disorder (BD), and schizophrenia (SCZ) are associated with an increased risk of cardiovascular diseases, including venous thromboembolism (VTE). The reasons for this are complex and include obesity, smoking, and use of hormones and psychotropic medications. Genetic studies have increasingly provided evidence of the shared genetic risk of psychiatric and cardiometabolic illnesses. This study aimed to determine whether a genetic predisposition to MDD, BD, or SCZ is associated with an increased risk of VTE. Genetic correlations using the largest genome-wide genetic meta-analyses summary statistics for MDD, BD, and SCZ (Psychiatric Genetics Consortium) and a recent genome-wide genetic meta-analysis of VTE (INVENT Consortium) demonstrated a positive association between VTE and MDD but not BD or SCZ. The same summary statistics were used to construct polygenic risk scores for MDD, BD, and SCZ in UK Biobank participants of self-reported White British ancestry. These were assessed for impact on self-reported VTE risk (10 786 cases, 285 124 controls), using logistic regression, in sex-specific and sex-combined analyses. We identified significant positive associations between polygenic risk for MDD and the risk of VTE in men, women, and sex-combined analyses, independent of the known risk factors. Secondary analyses demonstrated that this association was not driven by those with lifetime experience of mental illness. Meta-analyses of individual data from 6 additional independent cohorts replicated the sex-combined association. This report provides evidence for shared biological mechanisms leading to MDD and VTE and suggests that, in the absence of genetic data, a family history of MDD might be considered when assessing the risk of VTE.