A series of in situ synchrotron X-ray diffraction (XRD) measurements were carried out, combined with first-principles calculations, to study structural phase transitions of selenium at high pressures and room temperature. Several phase transitions were observed, among which an isostructural phase transition was found at around 120 GPa for the first time. Evolved from the rhombohedral (space group R
3
m) structure (Se-V), the new phase (Se-V′) exhibited an interesting increase of lattice parameter a at pressures from 120 to 148 GPa, known as negative linear compressibility (NLC). The discovery of NLC behavior observed in this work is mainly attributed to the accuracy and fine steps controlled by the membrane system for in situ XRD data collected with an exposure time of 0.5 s. After 140 GPa, a body-centered cubic (b.c.c.) structure Se-VI (space group Im
3
m) was formed, which remains stable up to 210 GPa, the highest pressure achieved in this study. The bulk moduli of phases Se-V, Se-V′ and Se-VI were estimated to be 83 ± 2, 321 ± 2 and 266 ± 7 GPa, respectively, according to the P–V curve fit by the third-order Birch–Murnaghan equation of state. The Se-V′ phase shows a bulk modulus almost 4 times larger than that of the Se-V phase, which is mainly due to the effect of its NLC. NLC in a higher pressure range is always more significant in terms of fundamental mechanism and new materials discovery, yet it has barely been reported at pressures above 100 GPa. This will hopefully inspire future studies on potential NLC behaviors in other materials at ultra-high pressure.