From $5.5$ months of dual-band optical photometric monitoring at the $1$ mmag level, BRITE-Constellation has revealed two simultaneous types of variability in the O4I(n)fp star $\zeta$ Puppis: one single periodic non-sinusoidal component superimposed on a stochastic component. The monoperiodic component is the $1.78$ d signal previously detected by Coriolis/SMEI, but this time along with a prominent first harmonic. The shape of this signal changes over time, a behaviour that is incompatible with stellar oscillations but consistent with rotational modulation arising from evolving bright surface inhomogeneities. By means of a constrained non-linear light curve inversion algorithm we mapped the locations of the bright surface spots and traced their evolution. Our simultaneous ground-based multi-site spectroscopic monitoring of the star unveiled cyclical modulation of its He II $\lambda4686$ wind emission line with the $1.78$-day rotation period, showing signatures of Corotating Interaction Regions (CIRs) that turn out to be driven by the bright photospheric spots observed by BRITE. Traces of wind clumps are also observed in the He II $\lambda4686$ line and are correlated with the amplitudes of the stochastic component of the light variations probed by BRITE at the photosphere, suggesting that the BRITE observations additionally unveiled the photospheric drivers of wind clumps in $\zeta$ Pup and that the clumping phenomenon starts at the very base of the wind. The origins of both the bright surface inhomogeneities and the stochastic light variations remain unknown, but a subsurface convective zone might play an important role in the generation of these two types of photospheric variability.

This paper aims to precisely determine the masses and detect pulsation modes in the two massive components of Beta Cen with BRITE-Constellation photometry. In addition, seismic models for the components are considered and the effects of fast rotation are discussed. This is done to test the limitations of seismic modeling for this very difficult case. A simultaneous fit of visual and spectroscopic orbits is used to self-consistently derive the orbital parameters, and subsequently the masses, of the components. The derived masses are equal to 12.02 +/- 0.13 and 10.58 +/- 0.18 M_Sun. The parameters of the wider, A - B system, presently approaching periastron passage, are constrained. Analysis of the combined blue- and red-filter BRITE-Constellation photometric data of the system revealed the presence of 19 periodic terms, of which eight are likely g modes, nine are p modes, and the remaining two are combination terms. It cannot be excluded that one or two low-frequency terms are rotational frequencies. It is possible that both components of Beta Cen are Beta Cep/SPB hybrids. An attempt to use the apparent changes of frequency to distinguish which modes originate in which component did not succeed, but there is potential for using this method when more BRITE data become available. Agena seems to be one of very few rapidly rotating massive objects with rich p- and g-mode spectra, and precisely known masses. It can therefore be used to gain a better understanding of the excitation of pulsations in relatively rapidly rotating stars and their seismic modeling. Finally, this case illustrates the potential of BRITE-Constellation data for the detection of rich-frequency spectra of small-amplitude modes in massive pulsating stars.