A closed-loop supply chain structure organises material and information flows from origin points to consumption points, including production, recycling, disposal, and other reverse logistic activities. Some integration problems arise with this structure including production, inventory, location, routing, distribution, collection, recycling, and routing. The integration problems that are facing scientific researchers include inventory routing, location routing, and location inventory. This study considers the integration problem of a closed-loop supply chain for the production, distribution, collection, and recycling quantities, along with the distribution and collection routes for each time period of a finite planning horizon. We refer to this problem as the "Closed-Loop Supply Chain Integrated Production-Inventory-Distribution-Routing Problem" (CLSC-PRP). A mathematical model is proposed that is the first to determine both quantities and routes for the CLSC-PRP simultaneously. As the problem is known to be NP-hard in terms of computational complexity, a simulated annealing-based decomposition heuristic is developed for solving large-scale CLSC-PRP instances. The results of the proposed mathematical model for the CLSC-PRP are compared with the results of the developed heuristic and two separate models that manage forward and backward production routing problems. An extensive comparative study indicated the following: (i) the proposed model was able to reduce the cost required for operating the total supply chain by an average of 12%, along with providing a positive impact on the environment and (ii) the proposed heuristic is able to generate solutions that are close to optimal in most cases.