Introduction: A previous study of 200,000 exome-sequenced UK Biobank participants to test for association of rare coding variants with hypertension implicated two genes at exome-wide significance, DNMT3A and FES. A total of 42 genes had an uncorrected p value <0.001. These results were followed up in a larger sample of 470,000 exome-sequenced participants. Methods: Weighted burden analysis of rare coding variants in a new sample of 97,050 cases and 172,263 controls was carried out for these 42 genes. Those showing evidence for association were then analysed in the combined sample of 167,127 cases and 302,691 controls. Results: The association of DNMT3A and FES with hypertension was replicated in the new sample and they and the previously implicated gene NPR1, which codes for a membrane-bound guanylate cyclase, were all exome-wide significant in the combined sample. Also exome-wide significant as risk genes for hypertension were GUCY1A1, ASXL1, and SMAD6, while GUCY1B1 had a nominal p value of <0.0001. GUCY1A1 and GUCY1B1 code for subunits of a soluble guanylate cyclase. For two genes, DBH, which codes for dopamine beta hydroxylase, and INPPL1, rare coding variants predicted to impair gene function were protective against hypertension, again with exome-wide significance. Conclusion: The findings offer new insights into biological risk factors for hypertension which could be the subject of further investigation. In particular, genetic variants predicted to impair the function of either membrane-bound guanylate cyclase, activated by natriuretic peptides, or soluble guanylate cyclase, activated by nitric oxide, increase risk of hypertension. Conversely, variants impairing the function of dopamine beta hydroxylase, responsible for the synthesis of norepinephrine, reduce hypertension risk.