Nitric oxide (NO) is essential to many physiological functions and operates in several signaling pathways. It is not understood how and when the different isoforms of nitric oxide synthase (NOS), the enzyme responsible for NO production, evolved in metazoans. This study investigates the number and structure of metazoan NOS enzymes by genome data mining and direct cloning of Nos genes from the lamprey. In total, 181 NOS proteins are analyzed from 33 invertebrate and 63 vertebrate species. Comparisons among protein and gene structures, combined with phylogenetic and syntenic studies, provide novel insights into how NOS isoforms arose and diverged. Protein domains and gene organization--that is, intron positions and phases--of animal NOS are remarkably conserved across all lineages, even in fast-evolving species. Phylogenetic and syntenic analyses support the view that a proto-NOS isoform was recurrently duplicated in different lineages, acquiring new structural configurations through gains and losses of protein motifs. We propose that in vertebrates a first duplication took place after the agnathan-gnathostome split followed by a paralog loss. A second duplication occurred during early tetrapod evolution, giving rise to the three isoforms--I, II, and III--in current mammals. Overall, NOS family evolution was the result of multiple gene and genome duplication events together with changes in protein architecture.