Abstract
Smoking is the leading cause of preventable death worldwide. Accordingly, effort has been devoted to determining the genetic variants that contribute to smoking risk. Genome-wide association studies have identified several variants in nicotinic acetylcholine receptor genes that contribute to nicotine dependence risk. We previously undertook pooled sequencing of the coding regions and flanking sequence of the CHRNA5, CHRNA3, CHRNB4, CHRNA6 and CHRNB3 genes and found that rare missense variants at conserved residues in CHRNB4 are associated with reduced risk of nicotine dependence among African Americans. We identified 10 low frequency (<5%) non-synonymous variants in CHRNB4 and investigated functional effects by co-expression with normal α3 or α4 subunits in human embryonic kidney cells. Voltage-clamp was used to obtain acetylcholine and nicotine concentration-response curves and qRT-PCR, western blots and cell-surface ELISAs were performed to assess expression levels. These results were used to functionally weight genetic variants in a gene-based association test. We find that there is a highly significant correlation between carrier status weighted by either acetylcholine EC50 (β = -0.67, r2 = 0.017, P = 2 × 10(-4)) or by response to low nicotine (β = -0.29, r2 = 0.02, P = 6 × 10(-5)) when variants are expressed with the α3 subunit. In contrast, there is no significant association when carrier status is unweighted (β = -0.04, r2 = 0.0009, P = 0.54). These results highlight the value of functional analysis of variants and the advantages to integrating such data into genetic studies. They also suggest that an increased sensitivity to low concentrations of nicotine is protective from the risk of developing nicotine dependence.