Complex behaviors are often observed at a spectrum in the population, and psychiatric disorders represent extremes of such behavioral spectra. While grasping the underlying cellular and molecular basis of these disorders represents a major challenge, it is believed that studies of complex behaviors in model organisms, where genotyping and phenotyping can be more conveniently carried out and cause-effect relationships can be further discerned, will help address this challenge. Here we report the characterization of a natural dark aversion behavior in larval zebrafish, which is previously shown to be fear or anxiety-associated. Phenotyping ∼200 individuals using a light/dark choice assay uncovered that, while a majority of individuals displayed medium level of dark aversion (mda), a small number of individuals exhibited strong dark aversion (sda), and a third small cohort showed variable dark aversion (vda). Through selective breeding and phenotyping of the next generation, we demonstrated that both the sda and vda traits are heritable, with sda being invariable while vda being highly variable across multiple trials. Additionally, sda appears to be recessive and vda appears to be dominant over the common allele(s) in the population. Moreover, compared to vda, sda showed increased thigmotaxis (preference for the walls in an open field), another measure of anxiety. Together, these findings reveal a naturally heritable variation of anxiety-like behavior in a tractable model organism, thereby laying foundation for future dissection of the underlying molecular and cellular mechanisms.