Physical properties of galaxies at z>7 are of interest for understanding both the early phases of star formation and the process of cosmic reionization. Chemical abundance measurements offer valuable information on the integrated star formation history, and hence ionizing photon production, as well as the rapid gas accretion expected at such high redshifts. We use reported measurements of [O III] 88$\mu$m emission and star formation rate to estimate gas-phase oxygen abundances in five galaxies at z=7.1-9.1 using the direct T_e method. We find typical abundances 12+log(O/H) = 7.9 ($\sim$0.2 times the solar value) and an evolution of 0.9$\pm$0.5 dex in oxygen abundance at fixed stellar mass from z$\simeq$8 to 0. These results are compatible with theoretical predictions, albeit with large (conservative) uncertainties in both mass and metallicity. We assess both statistical and systematic uncertainties to identify promising means of improvement with the Atacama Large Millimeter Array (ALMA) and the James Webb Space Telescope (JWST). In particular we highlight [O III] 52$\mu$m as a valuable feature for robust metallicity measurements. Precision of 0.1-0.2 dex in T_e-based O/H abundance can be reasonably achieved for galaxies at z$\approx$5-8 by combining [O III] 52$\mu$m with rest-frame optical strong lines. It will also be possible to probe gas mixing and mergers via resolved T_e-based abundances on kpc scales. With ALMA and JWST, direct metallicity measurements will thus be remarkably accessible in the reionization epoch.

Mass outflow rates and loading factors are typically used to infer the quenching potential of galactic-scale outflows. However, these generally rely on observations of a single gas phase which can severely underestimate the total ejected gas mass. To address this, we use observations of high mass ($\geqslant$10$^{10}$ M$_{\odot}$), normal star-forming galaxies at $z\sim$0 from the MaNGA, xCOLD GASS, xGASS and ALFALFA surveys and a stacking of NaD, H$\alpha$, CO(1-0) and HI 21cm tracers with the aim of placing constraints on an average, total mass outflow rate and loading factor. We find detections of outflows in both neutral and ionised gas tracers, with no detections in stacks of molecular or atomic gas emission. Modelling of the outflow components reveals velocities of $|$v$_{\text{NaD}}|$=131 km s$^{-1}$ and $|$v$_{\text{H}\alpha}|$=439 km s$^{-1}$ and outflow rates of $\dot{M}_{\text{NaD}}$=7.55 M$_{\odot}$yr$^{-1}$ and $\dot{M}_{\text{H}\alpha}$=0.10 M$_{\odot}$yr$^{-1}$ for neutral and ionised gas, respectively. Assuming a molecular/atomic outflow velocity of 200 km s$^{-1}$, we derive upper limits of $\dot{M}_{\text{CO}}<$19.43 M$_{\odot}$yr$^{-1}$ and $\dot{M}_{\text{HI}}<$26.72 M$_{\odot}$yr$^{-1}$ for the molecular and atomic gas, respectively. Combining the detections and upper limits, we find average total outflow rates of $\dot{M}_{\text{tot}}\lesssim$27 M$_{\odot}$yr$^{-1}$ and a loading factor of $\eta_{\text{tot}}\lesssim$6.39, with molecular gas likely contributing $\lesssim$72% of the total mass outflow rate, and neutral and ionised gas contributing $\sim$28% and $<$1%, respectively. Our results suggest that, to first order, a degree of quenching via ejective feedback could occur in normal galaxies when considering all gas phases, even in the absence of an AGN.

Outflows form an integral component in regulating the gas cycling in and out of galaxies, although their impact on the galaxy hosts is still poorly understood. Here we present an analysis of 405 high mass (log M$_{*}$/M$_{\odot}\geqslant10$), star-forming galaxies (excluding AGN) with low inclinations at $z\sim$0, using stacking techniques of the NaD $\lambda\lambda$5889,5895 A neutral gas tracer in IFU observations from the MaNGA DR15 survey. We detect outflows in the central regions of 78/405 galaxies and determine their extent and power through the construction of stacked annuli. We find outflows are most powerful in central regions and extend out to $\sim$1R$_{e}$, with declining mass outflow rates and loading factors as a function of radius. The stacking of spaxels over key galaxy quantities reveals outflow detections in regions of high $\Sigma_{\text{SFR}}$ ($\gtrsim$0.01 M$_{\odot}$yr$^{-1}$kpc$^{-2}$) and $\Sigma_{M_{*}}$ ($\gtrsim$10$^{7}$ M$_{\odot}$kpc$^{-2}$) along the resolved main sequence. Clear correlations with $\Sigma_{\text{SFR}}$ suggest it is the main regulator of outflows, with a critical threshold of $\sim$0.01 M$_{\odot}$yr$^{-1}$kpc$^{-2}$ needed to escape the weight of the disk and launch them. Furthermore, measurements of the H$\delta$ and D$_{n}$4000 indices reveal virtually identical star formation histories between galaxies with outflows and those without. Finally, through stacking of HI 21 cm observations for a subset of our sample, we find outflow galaxies show reduced HI gas fractions at central velocities compared to their non-detection control counterparts, suggestive of some removal of HI gas, likely in the central regions of the galaxies, but not enough to completely quench the host.

Prior to the launch of JWST, Spitzer/IRAC photometry offers the only means of studying the rest-frame optical properties of z>7 galaxies. Many such high redshift galaxies display a red [3.6] - [4.5] micron colour, often referred to as the "IRAC excess", which has conventionally been interpreted as arising from intense [OIII]+Hbeta emission within the [4.5] micron bandpass. An appealing aspect of this interpretation is similarly intense line emission seen in star-forming galaxies at lower redshift as well as the redshift-dependent behaviour of the IRAC colours beyond z~7 modelled as the various nebular lines move through the two bandpasses. In this paper we demonstrate that, given the photometric uncertainties, established stellar populations with Balmer (4000 A, rest-frame) breaks, such as those inferred at z>9 where line emission does not contaminate the IRAC bands, can equally well explain the redshift-dependent behaviour of the IRAC colours in 7<z<9 galaxies. We discuss possible ways of distinguishing between the two hypotheses using ALMA measures of [OIII] 88 micron and dust continuum fluxes. Prior to further studies with JWST, we show that the distinction is important in determining the assembly history of galaxies in the first 500 Myr.

We present the Super Eight galaxies - a set of very luminous, high-redshift ($7.1<z<8.0$) galaxy candidates found in Brightest of Reionizing Galaxies (BoRG) Survey fields. The original sample includes eight galaxies that are $Y$-band dropout objects with $H$-band magnitudes of $m_H<25.5$. Four of these objects were originally reported in Calvi et al. 2016. Combining new Hubble Space Telescope (HST) WFC3/F814W imaging and $Spitzer$ IRAC data with archival imaging from BoRG and other surveys, we explore the properties of these galaxies. Photometric redshift fitting places six of these galaxies in the redshift range of $7.1<z<8.0$, resulting in three new high-redshift galaxies and confirming three of the four high-redshift galaxy candidates from Calvi et al. 2016. We calculate the half-light radii of the Super Eight galaxies using the HST F160W filter and find that the Super Eight sizes are in line with typical evolution of size with redshift. The Super Eights have a mean mass of log(M$_*$/M$_\odot$) $\sim10$, which is typical for sources in this luminosity range. Finally, we place our sample on the UV $z\sim8$ luminosity function and find that the Super Eight number density is consistent with other surveys in this magnitude and redshift range.

We report the results of an expanded search for z~9-10 candidates over the ~883 arcmin^2 CANDELS+ERS fields. This study adds 147 arcmin^2 to the search area we consider over the CANDELS COSMOS, UDS, and EGS fields, while expanding our selection to include sources with bluer J_{125}-H_{160} colors than our previous J_{125}-H_{160}>0.5 mag selection. In searching for new z~9-10 candidates, we make full use of all available HST, Spitzer/IRAC, and ground-based imaging data. As a result of our expanded search and use of broader color criteria, 3 new candidate z~9-10 galaxies are identified. We also find again the z=8.683 source previously confirmed by Zitrin+2015. This brings our sample of probable z~9-11 galaxy candidates over the CANDELS+ERS fields to 19 sources in total, equivalent to 1 candidate per 47 arcmin^2 (1 per 10 WFC3/IR fields). To be comprehensive, we also discuss 28 mostly lower likelihood z~9-10 candidates, including some sources that seem to be reliably at z>8 using the HST+IRAC data alone, but which the ground-based data show are much more likely at z<4. One case example is a bright z~9.4 candidate COS910-8 which seems instead to be at z~2. Based on this expanded sample, we obtain a more robust LF at z~9 and improved constraints on the volume density of bright z~9 and z~10 galaxies. Our improved z~9-10 results again reinforce previous findings for strong evolution in the UV LF at z>8, with a factor of ~10 evolution seen in the luminosity density from z~10 to z~8.

We present Keck/MOSFIRE $H$-band spectroscopy targeting C$\,$III] $\lambda$1907,1909 in a $z=7.5056$ galaxy previously identified via Ly$\alpha$ emission. We detect strong line emission at $1.621\pm0.002\,\mu$m with a line flux of ($2.63\pm0.52$)$\times10^{-18}$ erg s$^{-1}$ cm$^{-2}$. We tentatively identify this line as [C$\,$III] $\lambda$1907, but we are unable to detect C$\,$III] $\lambda$1909 owing to sky emission at the expected location. This gives a galaxy systemic redshift, $z_{sys}=7.5032\pm0.0003$, with a velocity offset to Ly$\alpha$ of $\Delta$v$_{Ly\alpha}$ = $88\pm27$ km s$^{-1}$. The ratio of combined C$\,$III]/Ly$\alpha$ is 0.30-0.45, one of the highest values measured for any $z>2$ galaxy. We do not detect Si$\,$III] $\lambda\lambda$1883, 1892, and place an upper limit on Si$\,$III]/C$\,$III] $<$ 0.35 ($2\sigma$). Comparing our results to photoionization models, the C$\,$III] equivalent width (W$_{\text{CIII]}} = 16.23\pm2.32\,$\AA), low Si$\,$III]/C$\,$III] ratio, and high implied [O$\,$III] equivalent width (from the $Spitzer$/IRAC [3.6]$-$[4.5]$\simeq$0.8 mag color) require sub-Solar metallicities ($Z\simeq0.1-0.2 Z_{\odot}$) and a high ionization parameter, log$\,$U $\gtrsim -1.5$. These results favor models that produce higher ionization, such as the BPASS models for the photospheres of high-mass stars, and that include both binary stellar populations and/or an IMF that extends to 300 M$_{\odot}$. The combined C$\,$III] equivalent width and [3.6]$-$[4.5] color are more consistent with ionization from young stars than AGN, however we cannot rule out ionization from a combination of an AGN and young stars. We make predictions for $James~Webb~Space~Telescope$ spectroscopy using these different models, which will ultimately test the nature of the ionizing radiation in this source.

We perform a stacking analysis of the neutral \nad\,$\lambda\lambda$5889,5895\,\AA\ ISM doublet using the SDSS DR7 spectroscopic data set to probe the prevalence and characteristics of cold (T\,$\lesssim$\,10$^{4}$\,K) galactic-scale gas flows in local (0.025$\leqslant z\leqslant$0.1) inactive and AGN-host galaxies across the SFR-M$_{*}$ plane. We find low-velocity outflows to be prevalent in regions of high SFRs and stellar masses (10 $\lesssim$log M$_{*}$/M$_{\odot}$ $\lesssim$ 11.5), however we do not find any detections in the low mass (log M$_{*}$/M$_{\odot}$ $\lesssim$ 10) regime. We also find tentative detections of inflowing gas in high mass galaxies across the star-forming population. We derive mass outflow rates in the range of 0.14-1.74\,M$_{\odot}$yr$^{-1}$ and upper limits on inflow rates <1\,M$_{\odot}$yr$^{-1}$, allowing us to place constraints on the mass loading factor ($\eta$=$\dot{M}_{\text{out}}$/SFR) for use in simulations of the local Universe. We discuss the fate of the outflows by comparing the force provided by the starburst to the critical force needed to push the outflow outward, and find the vast majority of the outflows unlikely to escape the host system. Finally, as outflow detection rates and central velocities do not vary strongly with the presence of a (weak) active supermassive black hole, we determine that star formation appears to be the primary driver of outflows at $z\sim$0.

A fundamental quest of modern astronomy is to locate the earliest galaxies and study how they influenced the intergalactic medium a few hundred million years after the Big Bang. The abundance of star-forming galaxies is known to decline from redshifts of about 6 to 10, but a key question is the extent of star formation at even earlier times, corresponding to the period when the first galaxies might have emerged. Here we present spectroscopic observations of MACS1149-JD1, a gravitationally lensed galaxy observed when the Universe was less than four per cent of its present age. We detect an emission line of doubly ionized oxygen at a redshift of $9.1096\pm0.0006$, with an uncertainty of one standard deviation. This precisely determined redshift indicates that the red rest-frame optical colour arises from a dominant stellar component that formed about 250 million years after the Big Bang, corresponding to a redshift of about 15. Our results indicate the it may be possible to detect such early episodes of star formation in similar galaxies with future telescopes.

We introduce xCOLD GASS, a legacy survey providing a census of molecular gas in the local Universe. Building upon the original COLD GASS survey, we present here the full sample of 532 galaxies with CO(1-0) measurements from the IRAM-30m telescope. The sample is mass-selected in the redshift interval $0.01<z<0.05$ from SDSS, and therefore representative of the local galaxy population with M$_{\ast}>10^9$M$_{\odot}$. The CO(1-0) flux measurements are complemented by observations of the CO(2-1) line with both the IRAM-30m and APEX telescopes, HI observations from Arecibo, and photometry from SDSS, WISE and GALEX. Combining the IRAM and APEX data, we find that the CO(2-1) to CO(1-0) luminosity ratio for integrated measurements is $r_{21}=0.79\pm0.03$, with no systematic variations across the sample. The CO(1-0) luminosity function is constructed and best fit with a Schechter function with parameters {$L_{\mathrm{CO}}^* = (7.77\pm2.11) \times 10^9\,\mathrm{K\,km\,s^{-1}\, pc^{2}}$, $\phi^{*} = (9.84\pm5.41) \times 10^{-4} \, \mathrm{Mpc^{-3}}$ and $\alpha = -1.19\pm0.05$}. With the sample now complete down to stellar masses of $10^9$M$_{\odot}$, we are able to extend our study of gas scaling relations and confirm that both molecular gas fraction and depletion timescale vary with specific star formation rate (or offset from the star-formation main sequence) much more strongly than they depend on stellar mass. Comparing the xCOLD GASS results with outputs from hydrodynamic and semi-analytic models, we highlight the constraining power of cold gas scaling relations on models of galaxy formation.

The ubiquity of Lyman alpha (Ly$\alpha$) emission in a sample of four bright [O III]-strong star-forming galaxies with redshifts above 7 has led to the suggestion that such luminous sources represent a distinct population compared to their fainter, more numerous, counterparts. The presence of Ly$\alpha$ emission within the reionization era could indicate that these sources created early ionized bubbles due to their unusually strong radiation, possibly due to the presence of active galactic nuclei. To test this hypothesis we have secured long integration spectra with XSHOOTER on the VLT for three $z\simeq$7 sources selected to have similar luminosities and prominent excess fluxes in the IRAC 3.6 or 4.5$\mu$m band, usually attributed to strong [O III] emission. We secured additional spectroscopy for one of these galaxies at $z$=7.15 using MOSFIRE at the Keck telescope. For this, the most well-studied source in our sample with the strongest IRAC excess, we detect prominent nebular emission from He II and NV indicative of a non-thermal source. For the other two sources at $z$=6.81 and $z$=6.85, for which no previous optical/near infrared spectroscopy was initially available, Ly$\alpha$ is seen in one and CIII] emission in the other. Although a modest sample, our results further support the hypothesis that the phenomenon of intense [O III] emission is associated preferentially with sources lying in early ionized bubbles. However, even though one of our sources at $z$=7.15 clearly indicates the presence of non-thermal radiation, such ionized bubbles may not uniquely arise in this manner. We discuss the unique advantages of extending such challenging diagnostic studies with JWST.

We present \ci\,(2--1) and multi-transition $^{12}$CO observations of a dusty star-forming galaxy, ACT\,J2029+0120, which we spectroscopically confirm to lie at $z$\,=\,2.64. We detect CO(3--2), CO(5--4), CO(7--6), CO(8--7), and \ci\,(2--1) at high significance, tentatively detect HCO$^{+}$(4--3), and place strong upper limits on the integrated strength of dense gas tracers (HCN(4--3) and CS(7--6)). Multi-transition CO observations and dense gas tracers can provide valuable constraints on the molecular gas content and excitation conditions in high-redshift galaxies. We therefore use this unique data set to construct a CO spectral line energy distribution (SLED) of the source, which is most consistent with that of a ULIRG/Seyfert or QSO host object in the taxonomy of the \textit{Herschel} Comprehensive ULIRG Emission Survey. We employ RADEX models to fit the peak of the CO SLED, inferring a temperature of T$\sim$117 K and $n_{\text{H}_2}\sim10^5$ cm$^{-3}$, most consistent with a ULIRG/QSO object and the presence of high density tracers. We also find that the velocity width of the \ci\ line is potentially larger than seen in all CO transitions for this object, and that the $L'_{\rm C\,I(2-1)}/L'_{\rm CO(3-2)}$ ratio is also larger than seen in other lensed and unlensed submillimeter galaxies and QSO hosts; if confirmed, this anomaly could be an effect of differential lensing of a shocked molecular outflow.

We report on the detailed analysis of a gravitationally-lensed Y-band dropout, A2744_YD4, selected from deep Hubble Space Telescope imaging in the Frontier Field cluster Abell 2744. Band 7 observations with the Atacama Large Millimeter Array (ALMA) indicate the proximate detection of a significant 1mm continuum flux suggesting the presence of dust for a star-forming galaxy with a photometric redshift of $z\simeq8$. Deep X-SHOOTER spectra confirms the high redshift identity of A2744_YD4 via the detection of Lyman $\alpha$ emission at a redshift $z$=8.38. The association with the ALMA detection is confirmed by the presence of [OIII] 88$\mu$m emission at the same redshift. Although both emission features are only significant at the 4 $\sigma$ level, we argue their joint detection and the positional coincidence with a high redshift dropout in the HST images confirms the physical association. Analysis of the available photometric data and the modest gravitational magnification ($\mu\simeq2$) indicates A2744_YD4 has a stellar mass of $\sim$ 2$\times$10$^9$ M$_{\odot}$, a star formation rate of $\sim20$ M$_{\odot}$/yr and a dust mass of $\sim$6$\times$10$^{6}$ M$_{\odot}$. We discuss the implications of the formation of such a dust mass only $\simeq$200 Myr after the onset of cosmic reionisation.

We identify 4 unusually bright (H < 25.5) galaxies from HST and Spitzer CANDELS data with probable redshifts z ~ 7-9. These identifications include the brightest-known galaxies to date at z > 7.5. As Y-band observations are not available over the full CANDELS program to perform a standard Lyman-break selection of z > 7 galaxies, we employ an alternate strategy using deep Spitzer/IRAC data. We identify z ~ 7.1 - 9.1 galaxies by selecting z >~ 6 galaxies from the HST CANDELS data that show quite red IRAC [3.6]-[4.5] colors, indicating strong [OIII]+Hbeta lines in the 4.5 micron band. This selection strategy was validated using a modest sample for which we have deep Y-band coverage, and subsequently used to select the brightest z > 7 sources. Applying the IRAC criteria to all HST-selected optical-dropout galaxies over the full ~900 arcmin**2 of the CANDELS survey revealed four unusually bright z ~ 7.1, 7.6, 7.9 and 8.6 candidates. The median [3.6]-[4.5] color of our selected z ~ 7.1-9.1 sample is consistent with rest-frame [OIII]+Hbeta EWs of ~1500A, in the [4.5] band. Keck/MOSFIRE spectroscopy has been independently reported for two of our selected sources, showing Ly-alpha at redshifts of 7.7302+/-0.0006 and 8.683^+0.001_-0.004, respectively. We present similar Keck/MOSFIRE spectroscopy for a third selected galaxy with a probable 4.7sigma Ly-alpha line at z_spec=7.4770+/-0.0008. All three have H-band magnitudes of ~25 mag and are ~0.5 mag more luminous (M(UV) ~ -22.0) than any previously discovered z ~ 8 galaxy, with important implications for the UV LF. Our 3 brightest, highest redshift z > 7 galaxies all lie within the CANDELS EGS field, providing a dramatic illustration of the potential impact of field-to-field variance.

In this paper, we present a unique data set of more than one year's worth of regular observations of comet C/2013 A1(Siding Spring) with TRAPPIST in Chile, along with low-resolution spectra obtained with the ESO/VLT FORS 2 instrument. The comet made a close approach to Mars on October 19, 2014 and was then observed by many space and ground-based telescopes. We followed the evolution of the OH, NH, CN, $\mathrm{C_3}$, and $\mathrm{C_2}$ production rates as well as the $Af\rho$ parameter as a proxy for the dust production. We detected an outburst two weeks after perihelion, with gas and dust production rates being multiplied by a factor five within a few days. By modelling the shape of the CN and $\mathrm{C_2}$ radial profiles, we determined that the outburst happened around on November 10 around 15:30 UT ($\pm$ 5h) and measured a gas ejection velocity of $1.1\pm0.2$ km/s. We used a thermal evolution model to reproduce the activity pattern and outburst. Our results are consistent with the progressive formation of a dust mantle explaining the shallow dependence of gas production rates, which may be partially blown off during the outburst. We studied the evolution of gas composition, using various ratios such as CN/OH, $\mathrm{C_2}$/OH, or $\mathrm{C_3}$/OH, which showed little or no variation with heliocentric distance including at the time of the outburst. This indicates a relative level of homogeneity of the nucleus composition.

We report the discovery of Lyman-alpha emission (Ly$\alpha$) in the bright galaxy EGSY-2008532660 (hereafter EGSY8p7) using the MOSFIRE spectrograph at the Keck Observatory. First reported by Roberts-Borsani et al. (2015), it was selected for spectroscopic observations because of its photometric redshift ($z_{phot}=8.57^{+0.22}_{-0.43}$), apparent brightness (H$_{160}=25.26\pm0.09$) and red Spitzer/IRAC [3.6]-[4.5] color indicative of contamination by strong oxygen emission in the [4.5] band. With a total integration of $\sim$4.3 hours, our data reveal an emission line at $\simeq$11776 {\AA} which we argue is likely Ly$\alpha$ at a redshift $z_{spec}=8.683^{+0.001}_{-0.004}$, in good agreement with the photometric estimate. The line was detected independently on two nights using different slit orientations and its detection significance is $\sim7.5\sigma$. An overlapping skyline contributes significantly to the uncertainty on the total line flux although the significance of the detected line is robust to a variety of skyline-masking procedures. By direct addition and a Gaussian fit, we estimate a 95\% confidence range of 1.0--2.5$\times10^{-17}$ erg s$^{-1}$ cm$^{-2}$, corresponding to a rest-frame equivalent width of 17--42 {\AA}. EGSY8p7 is the most distant galaxy confirmed spectroscopically to date, and the third luminous source in the EGS field beyond $z_{phot}\gtrsim7.5$ with detectable Ly$\alpha$ emission viewed at a time when the intergalactic medium is believed to be fairly neutral. Although the reionization process was probably patchy, we discuss whether luminous sources with prominent IRAC color excesses may harbor harder ionizing spectra than the dominant fainter population thereby creating earlier ionized bubbles. Further spectroscopic follow-up of such bright sources promises important insight into the early formation of galaxies.

We present a spectroscopic redshift measurement of a very bright Lyman break galaxy at z=7.7302+-0.0006 using Keck/MOSFIRE. The source was pre-selected photometrically in the EGS field as a robust z~8 candidate with H=25.0 mag based on optical non-detections and a very red Spitzer/IRAC [3.6]-[4.5] broad-band color driven by high equivalent width [OIII]+Hbeta line emission. The Lyalpha line is reliably detected at 6.1 sigma and shows an asymmetric profile as expected for a galaxy embedded in a relatively neutral inter-galactic medium near the Planck peak of cosmic reionization. The line has a rest-frame equivalent width of EW0=21+-4 A and is extended with V_FWHM=360+90-70 km/s. The source is perhaps the brightest and most massive z~8 Lyman break galaxy in the full CANDELS and BoRG/HIPPIES surveys, having assembled already 10^(9.9+-0.2) M_sol of stars at only 650 Myr after the Big Bang. The spectroscopic redshift measurement sets a new redshift record for galaxies. This enables reliable constraints on the stellar mass, star-formation rate, formation epoch, as well as combined [OIII]+Hbeta line equivalent widths. The redshift confirms that the IRAC [4.5] photometry is very likely dominated by line emission with EW0(OIII+Hbeta)= 720-150+180 A. This detection thus adds to the evidence that extreme rest-frame optical emission lines are a ubiquitous feature of early galaxies promising very efficient spectroscopic follow-up in the future with infrared spectroscopy using JWST and, later, ELTs.