Human bocavirus (HBoV) is a recently discovered parvovirus associated with respiratory and gastroenteric infections in children. To date, four distinct subtypes have been identified worldwide. HBoV1 is the most frequently detected bocavirus in clinical samples derived from the respiratory tract. HBoV has a single-stranded DNA genome, which encodes two nonstructural proteins, NS1 and NP1, and two structural proteins, VP1 and VP2. Despite a large number of available HBoV sequences, the molecular evolution of this virus remains enigmatic. Here, we applied bioinformatic methods to measure the codon usage bias in 156 HBoV genomes and analyzed the factors responsible for preferential use of various synonymous codons. The effective number of codons (ENC) indicates a highly conserved, gene-specific codon usage bias in the HBoV genome. The structural genes exhibit a higher degree of codon usage bias than the non-structural genes. Natural selection emerged as dominant factor influencing the codon usage bias in the HBoV genome. Other factors that influence the codon usage include mutational pressure, gene length, protein properties, and the relative abundance of dinucleotides. The results presented in this study provide important insight into the molecular evolution of HBoV and may serve as a primer for HBoV gene expression studies and development of safe and effective vaccines to prevent infection.