We carried out $^{12}$CO($J$ = 1-0) observations of the Galactic gamma-ray supernova remnant (SNR) Kesteven 79 using the Nobeyama Radio Observatory 45 m radio telescope, which has an angular resolution of $\sim20$ arcsec. We identified molecular and atomic gas interacting with Kesteven 79 whose radial velocity is $\sim80$ km s$^{-1}$. The interacting molecular and atomic gases show good spatial correspondence with the X-ray and radio shells, which have an expanding motion with an expanding velocity of $\sim4$ km s$^{-1}$. The molecular gas associated with the radio and X-ray peaks also exhibits a high-intensity ratio of CO 3-2/1-0 $>$ 0.8, suggesting a kinematic temperature of $\sim24$ K, owing to heating by the supernova shock. We determined the kinematic distance to the SNR to be $\sim5.5$ kpc and the radius of the SNR to be $\sim8$ pc. The average interstellar proton density inside of the SNR is $\sim360$ cm$^{-3}$, of which atomic protons comprise only $\sim10$ $\%$. Assuming a hadronic origin for the gamma-ray emission, the total cosmic-ray proton energy above 1 GeV is estimated to be $\sim5 \times 10^{48}$ erg.

M16, the Eagle Nebula, is an outstanding HII region where extensive high-mass star formation is taking place in the Sagittarius Arm, and hosts the remarkable "pillars" observed with HST. We made new CO observations of the region in the 12CO J=1--0 and J=2--1 transitions with NANTEN2. These observations revealed for the first time that a giant molecular cloud of $\sim 1.3 \times 10^5$ \Msun \ is associated with M16, which is elongated vertically to the Galactic plane over 35 pc at a distance of 1.8 kpc. We found a cavity of the molecular gas of $\sim 10$ pc diameter toward the heart of M16 at \lbeq (16.95\degree, 0.85\degree), where more than 10 O-type stars and $\sim 400$ stars are associated, in addition to a close-by molecular cavity toward a Spitzer bubble N19 at \lbeq (17.06\degree, 1.0\degree). We found three velocity components which show spatially complementary distribution in the entire M16 giant molecular cloud (GMC) including NGC6611 and N19, suggesting collisional interaction between them. Based on the above results we frame a hypothesis that collision between the red-shifted and blue-shifted components at a relative of $\sim 10$ \kms \ triggered formation of the O-type stars in the M16 GMC in the last 1-2 Myr. The collision is two fold in the sense that one of the collisional interactions is major toward the M16 cluster and the other toward N19 with a RCW120 type, the former triggered most of the O star formation with almost full ionization of the parent gas, and the latter an O star formation in N19.

We analyzed the 2.6-mm CO and 21-cm HI lines toward the Magellanic superbubble 30 Doradus C, in order to reveal the associated molecular and atomic gas. We uncovered five molecular clouds in a velocity range from 251 to 276 km s$^{-1}$ toward the western shell. The non-thermal X-rays are clearly enhanced around the molecular clouds on a pc scale, suggesting possible evidence for magnetic field amplification via shock-cloud interaction. The thermal X-rays are brighter in the eastern shell, where there are no dense molecular or atomic clouds, opposite to the western shell. The TeV $\gamma$-ray distribution may spatially match the total interstellar proton column density as well as the non-thermal X-rays. If the hadronic $\gamma$-ray is dominant, the total energy of the cosmic-ray protons is at least $\sim1.2 \times 10^{50}$ erg with the estimated mean interstellar proton density $\sim60$ cm$^{-3}$. In addition the $\gamma$-ray flux associated with the molecular cloud (e.g., MC3) could be detected and resolved by the Cherenkov Telescope Array (CTA). This should permit CTA to probe the diffusion of cosmic-rays into the associated dense ISM.