Medical imaging is playing an increasingly crucial role in disease diagnosis. Numerous deep learning-based methods have been developed for imageguided automatic disease diagnosis. Most of the methods have harnessed conventional convolutional neural networks, which are directly applied in the regular image domain. However, some irregular spatial patterns revealed in medical images are also critical to disease diagnosis, since they can describe latent relations in different image regions of a subject (e.g., different focal lesions in an image) or between different groups (e.g., Alzheimer's disease and healthy control). Therefore, how to exploit and analyze irregular spatial patterns and their relations has become a research challenge in the field of imageguided disease diagnosis. To address this challenge, graph neural networks (GNNs) are proposed to perform the convolution operation on graphs. Graphs can naturally represent irregular spatial structures. Because of their ability to aggregate node features, edge features, and graph structure information to capture and learn hidden spatial patterns in irregular structures, GNN-based algorithms have achieved promising results in the detection of various diseases. In this paper, we introduce commonly used GNN-based algorithms and systematically review their applications to disease diagnosis. We summarize the workflow of GNN-based applications in disease diagnosis, ranging from localizing the regions of interest and edge construction to modeling. Furthermore, we discuss the limitations and outline potential research directions for GNNs in disease diagnosis.