The construction of the International Thermonuclear Experimental Reactor (ITER) has started and the Broader Approach (BA) activities aiming at early realization of fusion energy are in progress.  In order to accurately predict the behavior of burning fusion plasmas and to develop reliable schemes controlling them, the development of integrated simulation codes for burning plasmas is urgently needed.  Fusion reactions and external plasma heating and control modify the velocity distribution functions of electrons and ions and affect transport phenomena and various instabilities as well as the heating and current-drive efficiencies. The aim of this research project is to develop a kinetic integrated simulation code based on the time evolution of the velocity distribution functions and, through the time-dependent integrated analysis of toroidal plasmas including non-axisymmetric effects, to conduct the performance prediction of burning plasmas and the development of control schemes. This project will also contribute to the international collaboration for the development of ITER integrated modeling codes and the domestic collaboration for the code development intended for computer simulations in the BA activities.