The zebrafish (Danio rerio) is used extensively as a model species for studies on vertebrate development and for assessing chemical effects on reproduction. Despite this, the molecular mechanisms controlling zebrafish reproduction are poorly understood. We analyzed the transcriptomic profiles of the gonads of individual zebrafish, using a 17k oligonucleotide microarray, to define the molecular basis of sex and reproductive status in sexually mature fish. The gonadal transcriptome differed substantially between sexes. Among the genes overexpressed in females, 11 biological processes were overrepresented including mitochondrion organization and biogenesis, and cell growth and/or maintenance. Among the genes overexpressed in males, six biological processes were overrepresented including protein biosynthesis and protein metabolism. Analysis of the expression of gene families known to be involved in reproduction identified a number of genes differentially expressed between ovaries and testes including a number of sox genes and genes belonging to the insulin-like growth factor and the activin-inhibin pathways. Real-time quantitative PCR confirmed the expression profiles for nine of the most differentially expressed genes and indicated that many transcripts are likely to be switched off in one of the sexes in the gonads of adult fish. Significant differences were seen between the gonad transcriptomes of individual reproductively active females reflecting their stage of maturation, whereas the testis transcriptomes were remarkably similar between individuals. In summary, we have identified molecular processes associated with (gonadal) sex specificity in breeding zebrafish and established a strong relationship between individual ovarian transcriptomes and reproductive status in females.