IL17RB genetic variants are associated with acamprosate treatment response in patients with alcohol use disorder: A proteomics-informed genomics study.

Acamprosate is a Food and Drug Administration (FDA) approved medication for the treatment of alcohol use disorder (AUD). However, only a subset of patients achieves optimal treatment outcomes. Currently, no biological measures are utilized to predict response to acamprosate treatment. We applied our established pharmaco-omics informed genomics strategy to identify potential biomarkers associated with acamprosate treatment response. Specifically, our previous open-label acamprosate clinical trial recruited 442 patients with AUD who were treated with acamprosate for three months. We first performed proteomics using baseline plasma samples to identify potential biomarkers associated with acamprosate treatment outcomes. Next, we applied our established "proteomics-informed genome-wide association study (GWAS)" research strategy, and identified 12 proteins, including interleukin-17 receptor B (IL17RB), associated with acamprosate treatment response.​ A GWAS for IL17RB concentrations identified several genome-wide significant signals. Specifically, the top hit single nucleotide polymorphism (SNP) rs6801605 with a minor allele frequency is 38% in the European American population mapped 4 kilobase (Kb) upstream of IL17RB, and intron 1 of the choline dehydrogenase (CHDH) gene on chromosome 3 (p: 4.8E-20). The variant genotype (AA) for the SNP rs6801605 was associated with lower IL17RB protein expression. In addition, we identified a series of genetic variants in IL17RB that were associated with acamprosate treatment outcomes. Furthermore, the variantgenotypes for all of those IL17RB SNPs were protective for alcohol relapse. Finally, we demonstrated that the basal level of mRNA expression of IL17RB was inversely correlated with those of nuclear factor-κB (NF-κB) subunits, and a significantly higher expression of NF-κB subunits was observed in AUD patients who relapsed to alcohol use. In summary, this study illustrates that IL17RB genetic variants might contribute to acamprosate treatment outcomes. This series of studies represents an important step toward generating functional hypotheses that could be tested to gain insight into mechanisms underlying acamprosate treatment response phenotypes. (The ClinicalTrials.gov Identifier: NCT00662571).

Genome-wide association analyses identify 95 risk loci and provide insights into the neurobiology of post-traumatic stress disorder.

Post-traumatic stress disorder (PTSD) genetics are characterized by lower discoverability than most other psychiatric disorders. The contribution to biological understanding from previous genetic studies has thus been limited. We performed a multi-ancestry meta-analysis of genome-wide association studies across 1,222,882 individuals of European ancestry (137,136 cases) and 58,051 admixed individuals with African and Native American ancestry (13,624 cases). We identified 95 genome-wide significant loci (80 new). Convergent multi-omic approaches identified 43 potential causal genes, broadly classified as neurotransmitter and ion channel synaptic modulators (for example, GRIA1, GRM8 and CACNA1E), developmental, axon guidance and transcription factors (for example, FOXP2, EFNA5 and DCC), synaptic structure and function genes (for example, PCLO, NCAM1 and PDE4B) and endocrine or immune regulators (for example, ESR1, TRAF3 and TANK). Additional top genes influence stress, immune, fear and threat-related processes, previously hypothesized to underlie PTSD neurobiology. These findings strengthen our understanding of neurobiological systems relevant to PTSD pathophysiology, while also opening new areas for investigation.

Multi-ancestry meta-analysis of tobacco use disorder identifies 461 potential risk genes and reveals associations with multiple health outcomes.

Tobacco use disorder (TUD) is the most prevalent substance use disorder in the world. Genetic factors influence smoking behaviours and although strides have been made using genome-wide association studies to identify risk variants, most variants identified have been for nicotine consumption, rather than TUD. Here we leveraged four US biobanks to perform a multi-ancestral meta-analysis of TUD (derived via electronic health records) in 653,790 individuals (495,005 European, 114,420 African American and 44,365 Latin American) and data from UK Biobank (ncombined = 898,680). We identified 88 independent risk loci; integration with functional genomic tools uncovered 461 potential risk genes, primarily expressed in the brain. TUD was genetically correlated with smoking and psychiatric traits from traditionally ascertained cohorts, externalizing behaviours in children and hundreds of medical outcomes, including HIV infection, heart disease and pain. This work furthers our biological understanding of TUD and establishes electronic health records as a source of phenotypic information for studying the genetics of TUD.

Differential Serum Levels of CACNA1C, Circadian Rhythm and Stress Response Molecules in Subjects with Bipolar Disorder: Associations with Genetic and Clinical Factors.

Background: Many patients with bipolar disorder (BD) do not respond to or have difficulties tolerating lithium and/or other mood stabilizing agents. There is a need for personalized treatments based on biomarkers in guiding treatment options. The calcium voltage-gated channel CACNA1C is a promising candidate for developing personalized treatments. CACNA1C is implicated in BD by genome-wide association studies and several lines of evidence suggest that targeting L-type calcium channels could be an effective treatment strategy. However, before such individualized treatments can be pursued, biomarkers predicting treatment response need to be developed. Methods: As a first step in testing the hypothesis that CACNA1C genotype is associated with serum levels of CACNA1C, we conducted ELISA measures on serum samples from 100 subjects with BD and 100 control subjects. Results: We observed significantly higher CACNA1C (p<0.01) protein levels in subjects with BD. The risk SNP (rs11062170) showed functional significance as subjects homozygous for the risk allele (CC) had significantly greater CACNA1C protein levels compared to subjects with one (p=0.013) or no copies (p=0.009). We observed higher somatostatin (SST) (p<0.003) protein levels and lower levels of the clock protein ARTNL (p<0.03) and stress signaling factor corticotrophin releasing hormone (CRH) (p<0.001) in BD. SST and PER2 protein levels were associated with both alcohol dependence and lithium response. Conclusions: Our findings represent the first evidence for increased serum levels of CACNA1C in BD. Along with altered levels of SST, ARNTL, and CRH our findings suggest CACNA1C is associated with circadian rhythm and stress response disturbances in BD.

Molecular mechanisms involved in alcohol craving, IRF3, and endoplasmic reticulum stress: a multi-omics study.

Alcohol use disorder (AUD) is the most prevalent substance use disorder worldwide. Acamprosate and naltrexone are anti-craving drugs used in AUD pharmacotherapy. However, molecular mechanisms underlying their anti-craving effect remain unclear. This study utilized a patient-derived induced pluripotent stem cell (iPSC)-based model system and anti-craving drugs that are used to treat AUD as "molecular probes" to identify possible mechanisms associated with alcohol craving. We examined the pathophysiology of craving and anti-craving drugs by performing functional genomics studies using iPSC-derived astrocytes and next-generation sequencing. Specifically, RNA sequencing performed using peripheral blood mononuclear cells from AUD patients with extreme values for alcohol craving intensity prior to treatment showed that inflammation-related pathways were highly associated with alcohol cravings. We then performed a genome-wide assessment of chromatin accessibility and gene expression profiles of induced iPSC-derived astrocytes in response to ethanol or anti-craving drugs. Those experiments identified drug-dependent epigenomic signatures, with IRF3 as the most significantly enriched motif in chromatin accessible regions. Furthermore, the activation of IRF3 was associated with ethanol-induced endoplasmic reticulum (ER) stress which could be attenuated by anti-craving drugs, suggesting that ER stress attenuation might be a target for anti-craving agents. In conclusion, we found that craving intensity was associated with alcohol consumption and treatment outcomes. Our functional genomic studies suggest possible relationships among craving, ER stress, IRF3 and the actions of anti-craving drugs.

Clinical and genetic contributions to medical comorbidity in bipolar disorder: a study using electronic health records-linked biobank data.

Bipolar disorder is a chronic and complex polygenic disease with high rates of comorbidity. However, the independent contribution of either diagnosis or genetic risk of bipolar disorder to the medical comorbidity profile of individuals with the disease remains unresolved. Here, we conducted a multi-step phenome-wide association study (PheWAS) of bipolar disorder using phenomes derived from the electronic health records of participants enrolled in the Mayo Clinic Biobank and the Mayo Clinic Bipolar Disorder Biobank. First, we explored the conditions associated with a diagnosis of bipolar disorder by conducting a phenotype-based PheWAS followed by LASSO-penalized regression to account for correlations within the phenome. Then, we explored the conditions associated with bipolar disorder polygenic risk score (BD-PRS) using a PRS-based PheWAS with a sequential exclusion approach to account for the possibility that diagnosis, instead of genetic risk, may drive such associations. 53,386 participants (58.7% women) with a mean age at analysis of 67.8 years (SD = 15.6) were included. A bipolar disorder diagnosis (n = 1479) was associated with higher rates of psychiatric conditions, injuries and poisonings, endocrine/metabolic and neurological conditions, viral hepatitis C, and asthma. BD-PRS was associated with psychiatric comorbidities but, in contrast, had no positive associations with general medical conditions. While our findings warrant confirmation with longitudinal-prospective studies, the limited associations between bipolar disorder genetics and medical conditions suggest that shared environmental effects or environmental consequences of diagnosis may have a greater impact on the general medical comorbidity profile of individuals with bipolar disorder than its genetic risk.

Polygenic Prediction of Cellular Immune Responses to Mumps Vaccine.

In this report, we provide a follow-up analysis of a previously published genome-wide association study of host genetic variants associated with inter-individual variations in cellular immune responses to mumps vaccine. Here we report the results of a polygenic score (PGS) analysis showing how common variants can predict mumps vaccine response. We found higher PGS for IFNγ, IL-2, and TNFα were predictive of higher post-vaccine IFNγ (p-value = 2e-6), IL-2 (p = 2e-7), and TNFα (p = 0.004) levels, respectively. Control of immune responses after vaccination is complex and polygenic in nature. Our results suggest that the PGS-based approach enables better capture of the combined genetic effects that contribute to mumps vaccine-induced immunity, potentially offering a more comprehensive understanding than traditional single-variant GWAS. This approach will likely have broad utility in studying genetic control of immune responses to other vaccines and to infectious diseases.

Pharmacogenomic overlap between antidepressant treatment response in major depression & antidepressant associated treatment emergent mania in bipolar disorder.

There is increasing interest in individualizing treatment selection for more than 25 regulatory approved treatments for major depressive disorder (MDD). Despite an inconclusive efficacy evidence base, antidepressants (ADs) are prescribed for the depressive phase of bipolar disorder (BD) with oftentimes, an inadequate treatment response and or clinical concern for mood destabilization. This study explored the relationship between antidepressant response in MDD and antidepressant-associated treatment emergent mania (TEM) in BD. We conducted a genome-wide association study (GWAS) and polygenic score analysis of TEM and tested its association in a subset of BD-type I patients treated with SSRIs or SNRIs. Our results did not identify any genome-wide significant variants although, we found that a higher polygenic score (PGS) for antidepressant response in MDD was associated with higher odds of TEM in BD. Future studies with larger transdiagnostic depressed cohorts treated with antidepressants are encouraged to identify a neurobiological mechanism associated with a spectrum of depression improvement from response to emergent mania.

Prevalence and associations of multiple hypnotic prescriptions in a clinical sample.

STUDY OBJECTIVES: We examined the prevalence of multiple hypnotic prescriptions and its association with clinical and demographic characteristics from the electronic health record (EHR) in the Mayo Clinic Biobank. METHODS: Adult participants enrolled in the Mayo Clinic Biobank with an EHR number of ≥ 1 year were included (n = 52,940). Clinical and demographic characteristics were compared between participants who were and were not prescribed any hypnotic approved for insomnia by the US Food and Drug Administration and/or trazodone and in those prescribed a single vs multiple (≥ 2) hypnotics. A phenotype-based, phenome-wide association study (PheWAS) examining associations between hypnotic prescriptions and diagnoses across the EHR was performed adjusting for demographic and other confounders. RESULTS: A total of 17,662 (33%) participants were prescribed at least 1 hypnotic and 5,331 (10%) received ≥ 2 hypnotics. Participants who were prescribed a hypnotic were more likely to be older, female, White, with a longer EHR, and a greater number of diagnostic codes (all P < .001). Those with multiple hypnotic prescriptions were more likely to be younger, female, with a longer EHR, and a greater number of diagnostic codes (all P < .001) compared with those prescribed a single hypnotic. The PheWAS revealed that participants with multiple hypnotic prescriptions had higher rates of mood disorders, anxiety disorders, suicidal ideation, restless legs syndrome, and chronic pain (all P < 1 e-10). CONCLUSIONS: Receiving multiple hypnotic prescriptions is common and associated with a greater prevalence of psychiatric, chronic pain, and sleep-related movement disorders. Future studies should examine potential genetic associations with multiple hypnotic prescriptions to personalize treatments for chronic insomnia. CITATION: Kolla BP, Mansukhani MP, Chakravorty S, Frank JA, Coombes BJ. Prevalence and associations of multiple hypnotic prescriptions in a clinical sample. J Clin Sleep Med. 2024;20(5):793-800.

Patterns of pharmacotherapy for bipolar disorder: A GBC survey.

OBJECTIVES: To understand treatment practices for bipolar disorders (BD), this study leveraged the Global Bipolar Cohort collaborative network to investigate pharmacotherapeutic treatment patterns in multiple cohorts of well-characterized individuals with BD in North America, Europe, and Australia. METHODS: Data on pharmacotherapy, demographics, diagnostic subtypes, and comorbidities were provided from each participating cohort. Individual site and regional pooled proportional meta-analyses with generalized linear mixed methods were conducted to identify prescription patterns. RESULTS: This study included 10,351 individuals from North America (n = 3985), Europe (n = 3822), and Australia (n = 2544). Overall, participants were predominantly female (60%) with BD-I (60%; vs. BD-II = 33%). Cross-sectionally, mood-stabilizing anticonvulsants (44%), second-generation antipsychotics (42%), and antidepressants (38%) were the most prescribed medications. Lithium was prescribed in 29% of patients, primarily in the Australian (31%) and European (36%) cohorts. First-generation antipsychotics were prescribed in 24% of the European versus 1% in the North American cohort. Antidepressant prescription rates were higher in BD-II (47%) compared to BD-I (35%). Major limitations were significant differences among cohorts based on inclusion/exclusion criteria, data source, and time/year of enrollment into cohort. CONCLUSIONS: Mood-stabilizing anticonvulsants, second-generation antipsychotics, and antidepressants were the most prescribed medications suggesting prescription patterns that are not necessarily guideline concordant. Significant differences exist in the prescription practices across different geographic regions, especially the underutilization of lithium in the North American cohorts and the higher utilization of first-generation antipsychotics in the European cohorts. There is a need to conduct future longitudinal studies to further explore these differences and their impact on outcomes, and to inform and implement evidence-based guidelines to help improve treatment practices in BD.

The genetic contribution to the comorbidity of depression and anxiety: a multi-site electronic health records study of almost 178 000 people.

BACKGROUND: Depression and anxiety are common and highly comorbid, and their comorbidity is associated with poorer outcomes posing clinical and public health concerns. We evaluated the polygenic contribution to comorbid depression and anxiety, and to each in isolation. METHODS: Diagnostic codes were extracted from electronic health records for four biobanks [N = 177 865 including 138 632 European (77.9%), 25 612 African (14.4%), and 13 621 Hispanic (7.7%) ancestry participants]. The outcome was a four-level variable representing the depression/anxiety diagnosis group: neither, depression-only, anxiety-only, and comorbid. Multinomial regression was used to test for association of depression and anxiety polygenic risk scores (PRSs) with the outcome while adjusting for principal components of ancestry. RESULTS: In total, 132 960 patients had neither diagnosis (74.8%), 16 092 depression-only (9.0%), 13 098 anxiety-only (7.4%), and 16 584 comorbid (9.3%). In the European meta-analysis across biobanks, both PRSs were higher in each diagnosis group compared to controls. Notably, depression-PRS (OR 1.20 per s.d. increase in PRS; 95% CI 1.18-1.23) and anxiety-PRS (OR 1.07; 95% CI 1.05-1.09) had the largest effect when the comorbid group was compared with controls. Furthermore, the depression-PRS was significantly higher in the comorbid group than the depression-only group (OR 1.09; 95% CI 1.06-1.12) and the anxiety-only group (OR 1.15; 95% CI 1.11-1.19) and was significantly higher in the depression-only group than the anxiety-only group (OR 1.06; 95% CI 1.02-1.09), showing a genetic risk gradient across the conditions and the comorbidity. CONCLUSIONS: This study suggests that depression and anxiety have partially independent genetic liabilities and the genetic vulnerabilities to depression and anxiety make distinct contributions to comorbid depression and anxiety.

Factors associated with suicide attempts in the antecedent illness trajectory of bipolar disorder and schizophrenia.

BACKGROUND: Factors associated with suicide attempts during the antecedent illness trajectory of bipolar disorder (BD) and schizophrenia (SZ) are poorly understood. METHODS: Utilizing the Rochester Epidemiology Project, individuals born after 1985 in Olmsted County, MN, presented with first episode mania (FEM) or psychosis (FEP), subsequently diagnosed with BD or SZ were identified. Patient demographics, suicidal ideation with plan, self-harm, suicide attempts, psychiatric hospitalizations, substance use, and childhood adversities were quantified using the electronic health record. Analyses pooled BD and SZ groups with a transdiagnostic approach given the two diseases were not yet differentiated. Factors associated with suicide attempts were examined using bivariate methods and multivariable logistic regression modeling. RESULTS: A total of 205 individuals with FEM or FEP (BD = 74, SZ = 131) were included. Suicide attempts were identified in 39 (19%) patients. Those with suicide attempts during antecedent illness trajectory were more likely to be female, victims of domestic violence or bullying behavior, and have higher rates of psychiatric hospitalizations, suicidal ideation with plan and/or self-harm, as well as alcohol, drug, and nicotine use before FEM/FEP onset. Based on multivariable logistic regression, three factors remained independently associated with suicidal attempts: psychiatric hospitalization (OR = 5.84, 95% CI 2.09-16.33, p < 0.001), self-harm (OR = 3.46, 95% CI 1.29-9.30, p = 0.014), and nicotine use (OR = 3.02, 95% CI 1.17-7.76, p = 0.022). CONCLUSION: Suicidal attempts were prevalent during the antecedents of BD and SZ and were associated with several risk factors before FEM/FEP. Their clinical recognition could contribute to improve early prediction and prevention of suicide during the antecedent illness trajectory of BD and SZ.

Antidepressant-Associated Treatment Emergent Mania: A Meta-Analysis to Guide Risk Modeling Pharmacogenomic Targets of Potential Clinical Value.

BACKGROUND: The purpose of this study was to review the association between the SLC6A4 5-HTTLPR polymorphism and antidepressant (AD)-associated treatment emergent mania (TEM) in bipolar disorder alongside starting a discussion on the merits of developing risk stratification models to guide when not to provide AD treatment for bipolar depression. METHODS: Studies that examined the association between clinical and genetic risk factors, specifically monoaminergic transporter genetic variation, and TEM were identified. A meta-analysis was performed using the odds ratio to estimate the effect size under the Der-Simonian and Laird model. RESULTS: Seven studies, referencing the SLC6A4 5-HTTLPR polymorphism and TEM (total N = 1578; TEM+ =594, TEM- = 984), of 142 identified articles were included. The time duration between the start of the AD to emergence of TEM ranged from 4 to 12 weeks. There was a nominally significant association between the s allele of the 5-HTTLPR polymorphism and TEM (odds ratio, 1.434; 95% confidence interval, 1.001-2.055; P = 0.0493; I2 = 52%). No studies have investigated norepinephrine or dopamine transporters. CONCLUSION: Although the serotonin transporter genetic variation is commercially available in pharmacogenomic decision support tools, greater efforts, more broadly, should focus on complete genome-wide approaches to determine genetic variants that may contribute to TEM. Moreover, these data are exemplary to the merits of developing risk stratification models, which include both clinical and biological risk factors, to guide when not to use ADs in bipolar disorder. Future studies will need to validate new risk models that best inform the development of personalized medicine best practices treating bipolar depression.

Racial differences in pathways to care preceding first episode mania or psychosis: a historical cohort prodromal study.

Background: There is evidence suggesting racial disparities in diagnosis and treatment in bipolar disorder (BD) and schizophrenia (SZ). The purpose of this study is to compare psychiatric diagnoses and psychotropic use preceding a first episode of mania (FEM) or psychosis (FEP) in racially diverse patients. Methods: Using a comprehensive medical records linkage system (Rochester Epidemiology Project, REP), we retrospectively identified individuals diagnosed with BD or SZ and a documented first episode of mania or psychosis. Illness trajectory before FEP/FEM were characterized as the time from first visit for a mental health complaint to incident case. Pathways to care and clinical events preceding FEP/FEM were compared based on subsequent incident case diagnosis (BD or SZ) and self-reported race (White vs. non-White). Results: A total of 205 (FEM = 74; FEP = 131) incident cases were identified in the REP. Duration of psychiatric antecedents was significantly shorter in non-White patients, compared to White patients (2.2 ± 4.3 vs. 7.4 ± 6.6 years; p < 0.001) with an older age at time of first visit for a mental health complaint (15.7 ± 6.3 vs. 11.1 ± 6.0 years; p = 0.005). There were no significant differences by race in FEM pathway to care or age of first seeking mental health. Overall non-White patients had lower rates of psychotropic use. Conclusion: These data are unable to ascertain reasons for shorter duration of psychiatric antecedents and later age of seeking care, and more broadly first age of initial symptom presentation. If symptoms are confirmed to be earlier than first time seeking care in both groups, it would be important to identify barriers that racial minorities face to access timely psychiatric care and optimize early intervention strategies.

Discovery of 95 PTSD loci provides insight into genetic architecture and neurobiology of trauma and stress-related disorders.

Posttraumatic stress disorder (PTSD) genetics are characterized by lower discoverability than most other psychiatric disorders. The contribution to biological understanding from previous genetic studies has thus been limited. We performed a multi-ancestry meta-analysis of genome-wide association studies across 1,222,882 individuals of European ancestry (137,136 cases) and 58,051 admixed individuals with African and Native American ancestry (13,624 cases). We identified 95 genome-wide significant loci (80 novel). Convergent multi-omic approaches identified 43 potential causal genes, broadly classified as neurotransmitter and ion channel synaptic modulators (e.g., GRIA1, GRM8, CACNA1E ), developmental, axon guidance, and transcription factors (e.g., FOXP2, EFNA5, DCC ), synaptic structure and function genes (e.g., PCLO, NCAM1, PDE4B ), and endocrine or immune regulators (e.g., ESR1, TRAF3, TANK ). Additional top genes influence stress, immune, fear, and threat-related processes, previously hypothesized to underlie PTSD neurobiology. These findings strengthen our understanding of neurobiological systems relevant to PTSD pathophysiology, while also opening new areas for investigation.

A genome-wide association study of antidepressant-induced mania.

Antidepressant-induced mania (AIM) is a side effect of antidepressant treatment that is characterized by mania or hypomania after the start of medication. It is likely polygenic, but its genetic component remains largely unexplored. We aim to conduct the first genome-wide association study of AIM in 814 bipolar disorder patients of European ancestry. We report no significant findings from our single-marker or gene-based analyses. Our polygenic risk score analyses also did not yield significant results with bipolar disorder, antidepressant response, or lithium response. Our suggestive findings on the hypothalamic-pituitary-adrenal axis and the opioid system in AIM require independent replications.

Multi-ancestry meta-analysis of tobacco use disorder prioritizes novel candidate risk genes and reveals associations with numerous health outcomes.

Tobacco use disorder (TUD) is the most prevalent substance use disorder in the world. Genetic factors influence smoking behaviors, and although strides have been made using genome-wide association studies (GWAS) to identify risk variants, the majority of variants identified have been for nicotine consumption, rather than TUD. We leveraged five biobanks to perform a multi-ancestral meta-analysis of TUD (derived via electronic health records, EHR) in 898,680 individuals (739,895 European, 114,420 African American, 44,365 Latin American). We identified 88 independent risk loci; integration with functional genomic tools uncovered 461 potential risk genes, primarily expressed in the brain. TUD was genetically correlated with smoking and psychiatric traits from traditionally ascertained cohorts, externalizing behaviors in children, and hundreds of medical outcomes, including HIV infection, heart disease, and pain. This work furthers our biological understanding of TUD and establishes EHR as a source of phenotypic information for studying the genetics of TUD.

Whole Genome Analysis of Venous Thromboembolism: the Trans-Omics for Precision Medicine Program.

BACKGROUND: Risk for venous thromboembolism has a strong genetic component. Whole genome sequencing from the TOPMed program (Trans-Omics for Precision Medicine) allowed us to look for new associations, particularly rare variants missed by standard genome-wide association studies. METHODS: The 3793 cases and 7834 controls (11.6% of cases were individuals of African, Hispanic/Latino, or Asian ancestry) were analyzed using a single variant approach and an aggregate gene-based approach using our primary filter (included only loss-of-function and missense variants predicted to be deleterious) and our secondary filter (included all missense variants). RESULTS: Single variant analyses identified associations at 5 known loci. Aggregate gene-based analyses identified only PROC (odds ratio, 6.2 for carriers of rare variants; P=7.4×10-14) when using our primary filter. Employing our secondary variant filter led to a smaller effect size at PROC (odds ratio, 3.8; P=1.6×10-14), while excluding variants found only in rare isoforms led to a larger one (odds ratio, 7.5). Different filtering strategies improved the signal for 2 other known genes: PROS1 became significant (minimum P=1.8×10-6 with the secondary filter), while SERPINC1 did not (minimum P=4.4×10-5 with minor allele frequency <0.0005). Results were largely the same when restricting the analyses to include only unprovoked cases; however, one novel gene, MS4A1, became significant (P=4.4×10-7 using all missense variants with minor allele frequency <0.0005). CONCLUSIONS: Here, we have demonstrated the importance of using multiple variant filtering strategies, as we detected additional genes when filtering variants based on their predicted deleteriousness, frequency, and presence on the most expressed isoforms. Our primary analyses did not identify new candidate loci; thus larger follow-up studies are needed to replicate the novel MS4A1 locus and to identify additional rare variation associated with venous thromboembolism.