Genome-wide association studies (GWAS) are a powerful means of identifying genes with disease-associated common variants, but they are not well-suited to detecting genes with disease-associated rare and low-frequency variants. In the current study of Behçet disease (BD), nonsynonymous variants (NSVs) identified by deep exonic resequencing of 10 genes found by GWAS (
IL10, IL23R
,
CCR1
,
STAT4
,
KLRK1
,
KLRC1, KLRC2, KLRC3, KLRC4
, and
ERAP1
) and 11 genes selected for their role in innate immunity
(IL1B, IL1R1, IL1RN, NLRP3, MEFV, TNFRSF1A, PSTPIP1, CASP1, PYCARD, NOD2, and TLR4
) were evaluated for BD association. A differential distribution of the rare and low-frequency NSVs of a gene in 2,461 BD cases compared with 2,458 controls indicated their collective association with disease. By stringent criteria requiring at least a single burden test with study-wide significance and a corroborating test with at least nominal significance, rare and low-frequency NSVs in one GWAS-identified gene,
IL23R
(
P
= 6.9 × 10
−5
), and one gene involved in innate immunity,
TLR4
(
P
= 8.0 × 10
−4
), were associated with BD. In addition, damaging or rare damaging
NOD2
variants were nominally significant across all three burden tests applied (
P
= 0.0063–0.045). Furthermore, carriage of the familial Mediterranean fever gene (
MEFV
) mutation Met694Val, which is known to cause recessively inherited familial Mediterranean fever, conferred BD risk in the Turkish population (OR, 2.65;
P
= 1.8 × 10
−12
). The disease-associated NSVs in
MEFV
and
TLR4
implicate innate immune and bacterial sensing mechanisms in BD pathogenesis.