BACKGROUND: Rare variants are believed to play a substantial role in the genetic architecture of mental disorders, particularly in coding regions. However, limited evidence supports the impact of rare variants on anxiety. METHODS: Using whole-exome sequencing data from 200,643 participants in the UK Biobank, we investigated the contribution of rare variants to anxiety. Firstly, we computed genetic risk score (GRS) of anxiety utilizing genotype data and summary data from a genome-wide association study (GWAS) on anxiety disorder. Subsequently, we identified individuals within the lowest 50% GRS, a subgroup more likely to carry pathogenic rare variants. Within this subgroup, we classified individuals with the highest 10% 7-item Generalized Anxiety Disorder scale (GAD-7) score as cases (N = 1869), and those with the lowest 10% GAD-7 score were designated as controls (N = 1869). Finally, we conducted gene-based burden tests and single-variant association analyses to assess the relationship between rare variants and anxiety. RESULTS: Totally, 47,800 variants with MAF ≤0.01 were annotated as non-benign coding variants, consisting of 42,698 nonsynonymous SNVs, 489 nonframeshift substitution, 236 frameshift substitution, 617 stop-gain and 40 stop-loss variants. After variation aggregation, 5066 genes were included in gene-based association analysis. Totally, 11 candidate genes were detected in burden test, such as RNF123 (P = 3.40 × 10), MOAP1(P = 4.35 × 10), CCDC110 (P = 5.83 × 10). Single-variant test detected 9 rare variants, such as rs35726701(RNF123)(P = 3.16 × 10) and rs16942615(CAMTA2) (P = 4.04 × 10). Notably, RNF123, CCDC110, DNAH2, and CSKMT gene were identified in both tests. CONCLUSIONS: Our study identified novel candidate genes for anxiety in protein-coding regions, revealing the contribution of rare variants to anxiety.