Transcription factors (TFs) bind to specific regions of DNA known as transcription factor binding sites (TFBSs) and modulate gene expression by interacting with the transcriptional machinery. TFBSs are typically located upstream of target genes, within a few thousand base pairs of the transcription start site. The binding of TFs to TFBSs influences the recruitment of the transcriptional machinery, thereby regulating gene transcription in a precise and specific manner. This chapter provides practical examples and case studies demonstrating the extraction of upstream gene regions from the genome, identification of TFBSs using PWMEnrich R/Bioconductor package, interpretation of results, and preparation of publication-ready figures and tables. The EOMES promoter is used as a case study for single DNA sequence analysis, revealing potential regulation by the LHX9-FOXP1 complex during embryonic development. Additionally, an example is presented on how to investigate TFBSs in the upstream regions of a group of genes, using a case study of differentially expressed genes in response to human parainfluenza virus type 1 (HPIV1) infection and interferon-beta. Key regulators identified in this context include the STAT1:STAT2 heterodimer and interferon regulatory factor family proteins. The presented protocol is designed to be accessible to individuals with basic computer literacy. Understanding the interactions between TFs and TFBSs provides insights into the complex transcriptional regulatory networks that govern gene expression, with broad implications for several fields such as developmental biology, immunology, and disease research.