During the plastic deformation of hexagonal metals, it is easy to generate the point defect clusters with complex shapes and configurations due to their anisotropic properties. The interactions among these clusters and between these clusters and moving dislocations significantly influence the physical and mechanical properties of hexagonal materials. However, none of these issues in particular concerning the evolutions of vacancy clusters, the formation of microvoids, and the crack nucleation and propagation, is comprehensively understood on an atomic scale. In the present work, we first employ the activation-relaxation technique, in combination with ab initio and interatomic potential calculations, to systematically investigate vacancy cluster configurations in titanium and the transformation between these clusters. The results indicate the stable and metastable configurations of vacancy clusters at various sizes and activation energies of their dissociation, combination and migration. It is found that the formation and migration energies decrease with the size of vacancy cluster increasing. Small vacancy clusters stabilize at configurations with special symmetry, while large clusters transform into microvoids or microcracks. High-throughput molecular dynamics simulations are subsequently employed to investigate the influences of these clusters on plastic deformation under tensile loading. The clusters are found to facilitate the crack nucleation by providing lower critical stress, which decreases with the size of the vacancy clusters increasing. Under tensile loading, cracks are first nucleated at small clusters and then grow up, while large clusters form microvoids and cracks directly grow up.