We present the discovery and statistical analysis of $12\;660$ spotted variable stars toward and inside the Galactic bulge from over two-decade-long Optical Gravitational Lensing Experiment (OGLE) data. We devise a new method of dereddening of individual stars toward the Galactic bulge where strong and highly nonuniform extinction is present. In effect, $11\;812$ stars were classified as giants and $848$ as dwarfs. Well defined correlations between the luminosity, variability amplitude and rotation period were found for the giants. Rapidly rotating dwarfs with periods $P \leq 2$ d show I-band amplitudes lower than 0.2 mag which is substantially less than the amplitudes of up to 0.8 mag observed in giants and slowly rotating dwarfs. We also notice that amplitudes of stars brighter than $I_0 \approx 16$ mag do not exceed 0.3-0.4 mag. We divide the stars into three groups characterized by correlation between light and color variations. The positive correlation is characteristic for stars that are cooler when fainter, which results from the variable coverage of the stellar surface with spots similar to the sunspots. The variability of stars that are cooler when brighter (negative correlation) can be characterized by chemical spots with overabundance of heavy elements inside and variable line-blanketing effect, which is observed in chemically peculiar stars. The null correlation may results from very high level of the magnetic activity with rapidly variable magnetic fields. This division is readily visible on the color-magnitude diagram (CMD), which suggests that it may depend on the radius of the stars. We detect 79 flaring objects and discuss briefly their properties. Among others, we find that relative brightening during flares is correlated with brightness amplitude.

We present the results of a study of the type II Cepheid ($P_{puls} = 4.974 d$) in the eclipsing binary system OGLE-LMC-T2CEP-098 ($P_{orb} = 397.2 d$). The Cepheid belongs to the peculiar W Vir group, for which the evolutionary status is virtually unknown. It is the first single-lined system with a pulsating component analyzed using the method developed by Pilecki et al. (2013). We show that the presence of a pulsator makes it possible to derive accurate physical parameters of the stars even if radial velocities can be measured for only one of the components. We have used four different methods to limit and estimate the physical parameters, eventually obtaining precise results by combining pulsation theory with the spectroscopic and photometric solutions. The Cepheid radius, mass and temperature are $25.3 \pm 0.2 R_\odot$, $1.51 \pm 0.09 M_\odot$ and $5300 \pm 100 K$, respectively, while its companion has similar size ($26.3 R_\odot$), but is more massive ($6.8 M_\odot$) and hotter ($9500 K$). Our best estimate for the p-factor of the Cepheid is $1.30 \pm 0.03$. The mass, position on the period-luminosity diagram, and pulsation amplitude indicate that the pulsating component is very similar to the Anomalous Cepheids, although it has a much longer period and is redder in color. The very unusual combination of the components suggest that the system has passed through a mass transfer phase in its evolution. More complicated internal structure would then explain its peculiarity.

Single star evolution does not allow extremely low-mass stars to cross the classical instability strip (IS) during the Hubble time. However, within binary evolution framework low-mass stars can appear inside the IS once the mass transfer (MT) is taken into account. Triggered by a discovery of low-mass 0.26 Msun RR Lyrae-like variable in a binary system, OGLE-BLG-RRLYR-02792, we investigate the occurrence of similar binary components in the IS, which set up a new class of low-mass pulsators. They are referred to as Binary Evolution Pulsators (BEPs) to underline the interaction between components, which is crucial for substantial mass loss prior to the IS entrance. We simulate a population of 500 000 metal-rich binaries and report that 28 143 components of binary systems experience severe MT (loosing up to 90% of mass), followed by at least one IS crossing in luminosity range of RR Lyrae (RRL) or Cepheid variables. A half of these systems enter the IS before the age of 4 Gyr. BEPs display a variety of physical and orbital parameters, with the most important being the BEP mass in range 0.2-0.8 Msun, and the orbital period in range 10-2500 d. Based on the light curve only, BEPs can be misclassified as genuine classical pulsators, and as such they would contaminate genuine RRL and classical Cepheid variables at levels of 0.8% and 5%, respectively. We state that the majority of BEPs will remain undetected and we discuss relevant detection limitations.