Molecular epidemiology of human parainfluenza viruses type 1 (HPIV1) was investigated. Samples were collected from patients hospitalized in Croatia during the three consecutive epidemic seasons (2011-2014). Results indicated co-circulation of two major genetic clusters of HPIV1. Samples from the current study refer to clades II and III in a phylogenetic tree of haemagglutinin-neuraminidase (HN) gene. Additional phylogenetic trees of fusion (F) and phosphoprotein (P) genes confirmed the topology. Analysis of nucleotide diversity of entire P, F and HN genes demonstrated similar values: 0.0255, 0.0236 and 0.0237, respectively. However, amino acid diversity showed F protein to be the most conserved, while P protein was the most tolerant to mutations. Potential N- and O-glycosylation sites suggested that HPIV1 HN protein is abundantly glycosylated, and a specific N-glycosylation pattern could distinguish between clades II and III. Analysis of potential O-glycosylation sites in F protein indicated that samples from this study have two potential O-glycosylation sites, while publicly available sequences have five potential sites. This study provides data on the molecular characterization and epidemic pattern of HPIV1 in Croatia.