Not yet reviewed by Pith; the record is open.

Extreme emission line galaxies (EELGs) are typically characterized by high equivalent widths (EWs) which are driven by elevated specific star formation rates (sSFR) in low-mass galaxies with subsolar metallicities and little dust. Such extreme systems are rare in the local universe, but the number density of EELGs increases with redshift. Such starburst galaxies are currently presumed to be the main drivers of hydrogen reionization over 5.5<z<15, which serves to motivate many of the searches for high-z EELGs. We aim to characterize the physical properties of a sample of ~730 EELGs at 4<z<9 photometrically selected from the CEERS survey using JWST/NIRCam. We validate our method and demonstrate the main physical properties of a subset of EELGs using NIRSpec spectra. We create synthetic NIRCam observations of EELGs using empirical templates based on ~2000 local metal-poor starbursts to select EELGs based on color-color criteria. We study their properties based on SED fitting and flux excess from emission lines in the photometric filters. Our sample has a mean stellar mass of $10^{7.84}$Msun with high sSFRs with a mean value of $10^{-7.03}$ yr$^{-1}$. We consider a delayed-$\tau$ model for the star formation history and find our sample of EELGs are young with a mean value of the time after the onset of star formation of 45Myr. We find that they have similar line ratios to local metal-poor starbursts with high log([OIII]/H$\beta$)>0.4-1 which indicates that star formation may be the dominant source of ionization. Based on the photometric fluxes, we find an increase of EW([OIII]+H$\beta$) with sSFR and $\Sigma_{SFR}$, and a decrease with age and stellar mass. The sample of EELGs can reach $\Sigma_{SFR}>$10Msun yr$^{-1}$kpc$^{-2}$ which indicate they are strong candidates of LyC leakers. Another indirect indicator is the high values of O32>5 that can be reached for some galaxies in the sample.