Evidence for a sub-Jovian planet in the young TWA 7 disk 

Planets are thought to form from dust and gas in protoplanetary disks, with debris disks being the remnants of planet formation. Aged a few million up to a few billion years, debris disks have lost their primordial gas, and their dust is produced by steady-state collisions between larger, rocky bodies1,2. Tens of debris disks, with sizes of tens, sometimes hundreds, of astronomical units have been resolved with high-spatial-resolution, high-contrast imagers at optical and near-infrared or (sub)millimetre interferometers3,4. They commonly show cavities, ring-like structures and gaps, which are often regarded as indirect signatures of the presence of planets that gravitationally interact with unseen planetesimals2,5. However, no planet responsible for these features has been detected yet, probably because of the limited sensitivity (typically 2–10MJ) of high-contrast imaging instruments (see, for example, refs.6,7,8,9) before the James Webb Space Telescope. Here we have used the unprecedented sensitivity of the James Webb Space Telescope’s Mid-Infrared Instrument10,11in the thermal infrared to search for such planets in the disk of the approximately 6.4-Myr-old star TWA 7. With its pole-on orientation, this three-ring debris disk is indeed ideally suited for such a detection. We unambiguously detected a source 1.5 arcsec from the star, which is best interpreted as a cold, sub-Jupiter-mass planet. Its estimated mass (about 0.3MJ) and position (about 52au, de-projected) can thoroughly account for the main disk structures.

The disk around TWA 7 is one of the youngest (6.4 ± 1 Myr old; ref.12) debris disks known to date. TWA 7 is a close (about 34 pc; ref.13), low-mass (0.46M☉; ref.14) member of the young TW Hydra association, sometimes classified as a weak-line, non-accreting T-Tauri star15. The disk, resolved by the Near Infrared Camera and Multi-Object Spectrometer (NICMOS) on the Hubble Space Telescope (HST)16, is one of the rare ones resolved around M stars. It is seen almost pole-on17,18, a most favourable configuration to precisely estimate its radial distribution and to look for planets. The most recent modelling of the disk surface density deduced from polarimetric data from the Spectro-Polarimetic High Contrast Imager for Exoplanets Research (SPHERE) on the Very Large Telescope (VLT)18includes a ring peaking at 28auand extending out to more than 100au(R1), a narrow (less than 7aufull-width at half-maximum) ring at 52au(R2) and a broader (more than 40aufull-width at half-maximum) structure (93au; R3; Figs. 1 and 4 from ref.18). No planet has been detected so far, with detection limits roughly estimated to 0.5–1 MJbeyond 50au(Extended Data Fig.6andSupplementary Information).

The coronagraphic images of TWA 7 taken with the F1140C filter (central wavelength = 11.3 μm, bandwidth = 0.8 μm) of the Mid-Infrared Instrument (MIRI) on the James Webb Space Telescope (JWST)19were obtained on 21 June 2024 during cycle 2 (ID 3662; principal investigator, A.-M.L.; Extended Data Table1). The details on the data reduction procedure are described in theMethods. The main critical step is the subtraction of the residual diffracted light leaking from the MIRI coronagraph using a reference star observed with the same set-up. This process is necessary to bring the contrast ratio with respect to the star to a level of 10−5–10−4beyond angular separations of about 0.5 arcsec. The final image, presented in Fig.1, reveals three sources within 10 arcsec from TWA 7, the properties of which are listed in Extended Data Table2. One source at about 4.7 arcsec from TWA 7 (position angle 107°) was classified as a stellar background source, already detected in ancillary optical data from the Space Telescope Imaging Spectrograph (STIS) on the HST as well as near-infrared data from the NICMOS on the HST and the SPHERE on the VLT. The second one, located about 6.7 arcsec east of TWA 7, and spatially resolved, has no counterpart in the data from the SPHERE on the VLT, or in the data from the NICMOS or the STIS on the HST (taking into account TWA 7 proper motion (−118.751 ± 0.023 mas yr−1, −19.648 ± 0.026 mas yr−1; ref.20)). Its location in the MIRI image is consistent with that of a bright source in Atacama Large Millimeter Array (ALMA) band 7 (346 GHz) data from 2016 (ref.21), given the proper motion of the TWA 7 system between the ALMA and MIRI observations. This object therefore has the characteristics of a highly reddened background source. It is reminiscent of the JWST observations of the HR 8799 multi-planetary system, for which az≈ 1 galaxy was detected in the MIRI data taken at both 10 and 15 μm, and in ALMA band 7 data as well, but never identified in the near-infrared22. The third source, located about 1.5 arcsec northwest of TWA 7 (about 51au, projected separation) is unique to these MIRI observations. Hence, this source (hereafter, CC#1) is extremely red, and is not compatible with any background or foreground star.

North is up, and east is left. The status of three identified sources is indicated. Note that the faint signal north of the background galaxy is an artefact. bgd, background.

No data are available to test whether this third source shares a common proper motion with the star. In this context, in the following, we discuss the possible nature of this object. The first origin that one can consider is a Solar System object. Yet, most Solar System objects have proper motion between 5 and 40 arcsec h−1(ref.23). Even very remote, low-proper-motion small Solar System objects such as the dwarf planets Eris (semi-major axis of about 68au) and Sedna (semi-major axis of about 510au) showed proper motion of 1.4 arcsec h−1(ref.24) and 1.7  arcsec h−1(ref.25), respectively, at the time of their discovery. No trail was observed during the 2-h-long exposures, and no apparent motion was observed between the two images recorded during the sequences taken 2 h apart, indicating that CC#1 has a proper motion of less of than about 0.05  arcsec h−1. A Solar System object with such a low proper motion would be located at more than 200au. For such a cold object, reflected light would dominate in the MIRI F1140C filter and would require a Neptune-like size to fit the flux measured for CC#1 (for a geometric albedo of 0.1–0.3). It was checked whether such a scenario could be compatible with the hypothetical Planet Nine26. On the basis of the constraints from the planetary ephemeris27and the predictions of the orbit of Planet Nine28, a Solar System origin for CC#1 can definitely be excluded.

The second possible origin is a background galaxy. Like the background galaxy seen east of TWA 7, CC#1 has no reported counterpart at optical and near-infrared wavelengths. However, in contrast to this galaxy, it has no detected counterpart in the ALMA band 7 data (see details inSupplementary Information). The detection of this unresolved source at 11.3 μm, its non-detection in ALMA band 7 and the measured upper limits at 1.6 μm in the data from the SPHERE on VLT (Supplementary Information) could still be compatible with intermediate-redshift star-forming galaxies. Using such galaxy templates at various redshifts and published galaxy counts in JWST fields of view, we estimated that the probability of finding one such galaxy in a region of 1.5 arcsec radius centred on TWA 7 is about 0.34% (Methods). This probability is low, albeit non-zero. The location of the source with respect to the disk structure, right in a circumstellar ring gap (see below), makes the galaxy hypothesis even more unlikely.

The third and last possible origin is a planet. A forward modelling approach is used to constrain the properties of this planet. Using the HADES model29, it is possible to find fits for the JWST photometric data point while accounting for the 5σupper limits provided by the high-contrast images at 1.59 μm (H2) and 1.67 μm (H3) from the SPHERE on the VLT (Fig.2). HADES considers the thermal evolution of the planet; it assumes that the planet and the stars are coeval. Atmospheric fits incorporating water clouds indicate a narrow range, regardless of other parameters, for the effective temperature between 305 and 335 K (Extended Data Fig.3and Extended Data Table3), and a mass of about 0.3MJ. A metallicity range above solar is required. Additional data will be necessary to further constrain this parameter.

Modelled photometry is shown with crosses. The 1–15-μm spectra correspond to representative solutions that fit the observed flux of TWA 7b in the JWST F1140C filter (blue), with respect to the 5σupper limits from the VLT SPHERE data in the H2 and H3 filters (red), and are consistent with an age of 6.4 ± 1 Myr. The bandwidths of H2 and H3 are 0.052 μm and 0.054 μm, respectively. A zoomed-in view around H2 and H3 highlights their bandwidths and the integrated model spectrum points below the upper limit. met, metallicity;Teff, effective temperature;fsed, sedimentation rate.

As a comparison and extra mass estimation, evolutionary models of cold, low-mass planets30combined with our estimation of effective temperature are used; they lead to a comparable mass of close to 0.3MJ(Extended Data Fig.4) for metallicities less than 2.5, as available in their framework. These two consistent results indicate that the planet mass is substantially below 1MJ. The current best estimate is around 0.3MJ, and is only weakly dependent on the underlying details of the two models used. It, however, depends on the age of the planet, assumed here to be coeval with its parent star. A younger planet would lead to a smaller mass.

The observed source is located right on the R2 narrow ring, and, moreover, within a region identified by ref.18as underdense compared to the rest of the ring (Fig.3a). This is very reminiscent of simulations of resonant rings predicted by early works31,32for closer and less massive planets, which led to the proposal of such a possible situation for the TWA 7 system18.

a, Polarimetric image (in log scale) of the disk composed of the sum of three epochs (26 April 2016 presented in ref.17, 20 March 2017 presented in ref.18, and a new epoch, 8 February 2022, reduced as in ref.18) from the SPHERE Infrared Dual-Band Imager and Spectrograph (IRDIS), with the MIRI image (resampled to the SPHERE pixel size) as an overlay with contours. The log of these data is provided in Supplementary Table2. The peak densities of the rings are also indicated. The central hatched disk is a numerical mask to hide the stellar residuals.b, Disk simulations. Top view of a disk of massless planetesimals perturbed by a 0.34-MJplanet at 52au, on a circular orbit, after 6 Myr (see details in text). The orbit of the perturbing planet is sketched in green and the location of the planet on its orbit is shown in red.

DedicatedN-body simulations were conducted for a planet with a mass of 0.34MJ, located at 52auaround the 0.46M☉central star. This value is consistent with the measured projected separation, assuming that the planet and the ≈13°-inclined disk are coplanar. The simulation also included a disk of 200,000 planetesimals, distributed between 20 and 130au. These parameters were selected to roughly match the boundaries of the observed disk. The planetesimal disk was assumed to be coplanar with the planet’s orbit, with an initial surface density proportional to 1/r. The eccentricities of the planetesimals were chosen to range between 0 and 0.01, and the planet was assumed to evolve on a circular orbit. The simulations, which spanned 10 Myr, were performed using the symplecticN-body code swift_rmv3 (ref.33), which provides a first-order treatment of close encounters.

A top view of the resulting distribution of planetesimals after 6 Myr is shown in Fig.3b. Planetary perturbations efficiently carve the disk over about 30au, but leave a narrow ring at 52au, as well as a relative void (underdensity) around the planet. The latter structure is characteristic of a ring of co-orbital planetesimals with the planets, trapped in a 1:1 resonance with it. The similarity between the TWA 7 disk image and the simulation (Fig.3a,b) is remarkable. In this context, the observed R2 ring would correspond to the ring of co-orbiters with the planet. We nevertheless expect some differences between both distributions owing to the radiation pressure acting on the grains. Additional information on the grain size distribution is needed to compute the effects of radiation pressure and refine the dynamical modelling.

The low likelihood of a background galaxy, the successful fit of the MIRI flux and SPHERE upper limits by a 0.3-MJplanet spectrum and the fact that an approximately 0.3-MJplanet at the observed position would naturally explain the structure of the R2 ring, its underdensity at the planet’s position and the gaps provide compelling evidence supporting a planetary origin for the observed source. Like the planet β Pictoris b, which is responsible for an inner warp in a well-resolved—from optical to millimetre wavelengths—debris disk34, TWA 7b is very well suited for further detailed dynamical modelling of disk–planet interactions. To do so, deep disk images at short and millimetre wavelengths are needed to constrain the disk properties (grain sizes and so on). Refining the planet mass determination can be carried out with additional JWST photometry and possibly spectroscopy. Measuring the orbital parameters (eccentricity, in particular) is more challenging given the long orbital period (about 550 yr) of the planet. Yet, one notes that a planet on an eccentric orbit would rapidly destroy R2.

As it is angularly well resolved from the star, TWA 7b is suited for direct spectroscopic investigations, providing the opportunity to study the interior and the atmosphere of a non-irradiated sub-Jupiter-mass, cold (about 320 K) exoplanet, and start comparative studies with our much older and cooler Solar System giants, as well as with the recently imaged cold (about 275 K) but more massive (6MJ) planet eps Ind Ab35. Improved estimations of its metallicity and temperature will further constrain its mass.

The present results show that the JWST MIRI has opened up a new window in the study of sub-Jupiter-mass planets using direct imaging. Indeed, TWA 7b (about 100M⊕) is at least ten times lighter than the exoplanets directly imaged so far, and planets as light as 25–30M⊕could have been detected if present at 1.5 arcsec from the star or beyond.

Coronagraphic observations were performed with the 4QPM_1140 coronagraph paired with the F1140C filter. The details of the observations are given in Extended Data Table1. We obtained two roll angles (difference of 7.835°) to mitigate the attenuation of the coronagraph in the field of view, in case an object falls close to one phase transition of the 4QPM. Each coronagraphic observation was 2 h long, hence a total of 4 h on the science target. Background observations were observed immediately after the science exposures in a two-point dithering mode, for a total of 4 h.

A reference star, CD-23-9765, was observed back-to-back with the target in the same configuration with the aim to subtract the starlight diffraction after the coronagraph. The reference shares similar brightness and spectral type with the target, and is angularly close. It was observed with nine-point small-grid dithering (SGD) to apply post-processing algorithms such as principal component analysis (PCA)36. In total, the reference star was observed for about 1 h and comes with dedicated background observations.

Using comparison with simulated coronagraphic images, as in ref.37, we were able to estimate a pointing accuracy on the 4QPM of about 2 mas per axis, significantly lower than the 10-mas step of the SGD. We also confirmed the detector coordinates of the 4QPM_1140 mask (119.758, 112.158 as provided in the JWST Calibration Reference Data System).

The data reduction follows the steps described in refs.21,38. Level 1 data are retrieved from the Mikulski Archive for Space Telescopes (MAST), processed with v1.14.0 of the pipeline together with Calibration Reference Data System file 1241. Images are registered to the coronagraph centre. Calibrated files (‘cal’ files) are produced in-house with the JWST pipeline for each roll, by subtracting the background and converting the photometric units (data numbers per second to mega-janskys per steradian). The background is built from the minimum per pixel of the two dithers. We skipped the flat-field correction, which is not appropriate for the MIRI coronagraph22.

We took advantage of the diversity brought by the SGD mode to build a reference frame to subtract the stellar diffraction. We tested various algorithms and retained a linear combination (which uses the downhill simplex minimization) of the nine SGD reference star images, as well as the PCA, as the two algorithms providing the best detection of CC#1. To mitigate the over-subtraction effect, a numerical masking is implemented to ignore some parts of the image. We obtained the best compromise by selecting the annular region between 0.5 arcsec and 3 arcsec, and the three sources were masked with a 1λ/Dpatch, withλrepresenting the observing wavelength andDthe telescope diameter (we checked that the CC#1 flux measurements are similar with a different region: 2–3 arcsec). We proceeded similarly for the PCA, using eight components to build the final image that is subtracted to the data, in an annular region from 0.5 arcsec to 5 arcsec (point sources not masked). Despite the bad-pixel correction applied on the raw coronagraphic images, the subtracted images with the reference star are still affected by a few bad pixels, both with the linear combination and PCA. We further apply aσ-clipping function to correct for these remaining bad pixels. Finally, the images are rotated to align north up, considering the aperture position angles: 121.45° for roll 1 and 129.27° for roll 2, as well as the V3 axis orientation on the detector (4.835°).

Next, extracting the flux and position of the CC#1 requires modelling its point spread function (PSF), which can vary spatially and nonlinearly owing to the attenuation of the coronagraph for which the phase transitions extend across the whole field of view. We used both the diffraction code developed in refs.39,40and WebbPSF41to simulate the MIRI PSF, taking into account the coronagraph, considering the configuration of mask, stop and filter (‘FQPM1140’, ‘MASKFQPM’ and ‘F1140C’, according to the WebbPSF terminology). The position of CC#1 is approximated with a Gaussian fit, passing the sky coordinates to the former diffraction code and detector coordinates to WebbPSF to calculate the PSF of CC#1 accounting for the coronagraph attenuation. We measured the coronagraph transmission with both PSF estimates. The flux of the PSF model at the position of CC#1 is integrated in a 1.5λ/Daperture (to match the aperture used for photometric measurements), and compared to that at 10 arcsec (far away from the coronagraph influence). The two approaches give similar transmissions: 0.66 and 0.62 in roll 1 and 0.31 and 0.28 in roll 2. Therefore, CC#1 is significantly closer to one 4QPM transition in roll 2 than in roll 1, so both its astrometry and photometry can be affected (Extended Data Fig.1). We measured a signal-to-noise ratio of 30 and 18, respectively, for roll 1 and roll 2, using the linear combination method. In the roll 2 image, the PSF is more asymmetrical as the planet is closer to the quadrant edge in comparison to roll 1, so we decided to consider only roll 1 data for the photometric analysis.

On the basis of these CC#1 PSF models, we extracted the flux and the photometry of the object by minimizing the residuals between the reduced JWST data and a PSF model in a 1.5λ/Darea with three free parameters (positions and flux) and using either a downhill simplex algorithm or the Nelder–Mead algorithm42. For comparison, we also used aperture photometry, but this required implementation of an aperture correction based on simulated PSF (ratio of the total flux in the PSF to the flux integrated in the 1.5λ/Daperture).

The flux extraction is applied both on the photometrically calibrated files (.cal), which directly provides CC#1 flux in mega-janskys per steradian, and on the uncalibrated files (.rate), which requires to measure a contrast with respect to the non-coronagraphic image of the star. As detailed previously22, the contrast measurement relies on commissioning data either on target acquisition images that come with the telescope pointing procedure but are obtained with a neutral density filter, or from images obtained on and off the coronagraph on another star. The method using target acquisition shows some net discrepancies, probably because the target acquisition filter is very broad (about 8–18 μm) and the targets have different spectral types (M3 for TWA 7 and K0 or K5 for the commissioning targets). Besides, the emission of the TWA 7 disk is expected to become significant beyond 15 μm. As a result, we did not use the target acquisition method in the following. The final photometric values are based on the linear combination and PCA technique to suppress the starlight. We averaged the values of the different methods (calibrated files and contrast, aperture and PSF model for the photometric extraction) and the error bar is built from the extreme values ((max − min)/2) for being conservative. We measured a flux density of 5.60 ± 0.97 × 10−19W m−2μm−1for CC#1. The fluxes of the other sources are given in Extended Data Table2. Note that we also tested the typical ‘injection–recovery’ method directly in the raw data, but did not notice any significant differences for the extracted flux of the planet.

To estimate the probability that the source labelled CC#1 is a galaxy, we have taken into account the three constraints on the fluxes or upper limits at 1.6, 11 and 870 μm. The source has a flux of 22 μJy at 11 μm, and is not detected with ALMA at 870 μm, but with a tapered resolution of 2 arcsec. The measured 3σupper limit for an unresolved source at the position of the MIRI source is 96 μJy (ref.21). Combining all ALMA observations (see Extended Data Fig.5), the 3σupper limit is 76 μJy, in a beam of 0.29 arcsec × 0.24 arcsec (see a detailed analysis and estimation in theSupplementary Information). The third constraint is from the non-detection at 1.6 μm with VLT SPHERE, with an upper limit of approximately 0.6 μJy for a point source (about 60 mas in size). However, we have to take into account the fact that a low-zgalaxy could be extended in the calculation of its maximum flux at 2 μm. Indeed, the size of a galaxy is expected to vary from one wavelength to the other: at 1.6 μm, the disk of old stars dominates, so the source is more extended, of the order of 10 kpc, whereas at 11 μm and 870 μm, the nuclear star-forming region dominates, and will appear more concentrated, of the order of 100–500 pc. We therefore considered the flux limit of an extended (0.6 arcsec) source in the SPHERE data (that is, 60 μJy).

Two types of galaxy might comply with these three constraints: star-forming galaxies or active nucleus (AGN) or a combination of both, with a redshift betweenz= 0.1 and 1. For such galaxies, at lower redshifts, where 1 arcsec is smaller than 2 kpc, the source would probably look extended (up to 5 arcsec) for wavelengths between 1 and 4 μm (ref.43). At higher redshifts, the peak of the emission usually located at 100 μm will enter the ALMA domain at about 1 mm, and it should have been detected by ALMA44.

As the density of the cosmic star formation rate increases considerably with redshift betweenz= 0 and 1, starbursts atz= 0 are good templates for star-forming galaxies atz= 0.1 to 1. The starburst is in general nuclear. It can be highly peaked in the centre, or distributed in a ring of about 100-pc radius (like in M82), but the emitting size at 100-μm rest frame will be no more than typically 1 kpc. This corresponds to an angular size of about 0.55 arcsec atz= 0.1, about 0.16 arcsec atz= 0.5, and about 0.12 arcsec atz= 1. We also considered AGN templates. The nucleus is then a point source at long wavelengths, while the galaxy host is extended (about 10 kpc) at the shortest wavelengths.

The probability to find az= 0.1 to 1 star-forming galaxy with a flux at 11 μm of 22 ± 3.8 μJy (hence, within a range of 7.6 μJy) in a field of view of 10 × 10 arcsec can be estimated readily through the source counts, from refs.45,46,47,48. A probability of 50% is found in the redshift rangez= 0.1 to 1.

To take into account the other constraints at 1.6 and 870 µm, we used the Spitzer Wide-Area Infrared Extragalactic Survey (SWIRE) templates, from ref.49,http://www.iasf-milano.inaf.it/~polletta/templates/swire_templates.htmlwhich include 14 templates of starbursts and AGN, representative of star-forming and active intermediate-redshift galaxies. These templates were computed for a half-dozen redshifts, and calibrated to have a flux of 22 μJy at 11 μm. An illustrative plot of their spectral energy distribution (SED) is shown in Extended Data Fig.2, in which the flux constraints at 1.6 and 870 µm are also indicated by black symbols (square and triangle).

We used the rest-frame 8-μm luminosity function of galaxies at redshifts betweenz= 0.1 andz= 1 in the Great Observatories Origins Deep Survey (GOODS) fields50to estimate the abundance of these star-forming galaxies. The 8-μm luminosity was deduced from the templates. This led to a probability of 12% to find such a galaxy in the right redshift range (z= 0.1 to 1) in the 10 × 10 arcsec field of view. We estimated that there are about 80% star-forming galaxies and 20% active nuclei, including low-luminosity ones such as Seyfert galaxies51,52. As can be seen in Extended Data Fig.2, all templates comply with the three constraints belowz= 0.5, and some starbursts are found brighter at 870 μm above this redshift. We therefore reduced the probability accordingly forz> 0.5, attributing each star-forming template equal probability, as shown by the observations53. The final probability to find such a galaxy in the 10 × 10 arcsec field of view is about 5%, with significant uncertainties: +3/−2%. Hence, in a radius of 1.5 arcsec around TWA 7, the probability to have such a galaxy is about 0.34% (+0.22/−0.14%). We note that the distribution of galaxies on the sky might be clustered in some rare places, and our error bars should be increased, but no more than 30% given the typical errors on luminosity function of galaxies45,47,50. Altogether, the probability of having a galaxy satisfying the various constraints within 1.5 arcsec would be 0.34% (+0.29/−0.18%). Note that the current calculation applies to these specific observational constraints. The probability for CC#1 to be a rare type of background galaxy compatible with the JWST, ALMA and SPHERE observational constraints is therefore low. Such a low-probability event can admittedly arise would a large number of independent datasets be analysed. However, this is not the case for this study. Indeed, only two datasets among about 30 JWST MIRI datasets dedicated to exoplanet searches are analysed. Moreover, most of the other datasets are not as deep as the present ones, and there are no corresponding 1.6- and 870-μm data available. Hence, the likelihood that any research group would have found such a peculiar background galaxy contaminant within the whole corpus of relevant JWST data is much lower than the simple 10.2% (30 × 0.34%) one could derive.

All observational data in this work are available through public data archives. The JWST data were collected through General Observer programme number 3662 (principal investigator A.-M.L.) and will be available through the MAST database (https://mast.stsci.edu) from 21 June 2025. The data described here are available via the MAST archive athttps://doi.org/10.17909/4qvz-sw62. The present study made use of ALMA data collected during the 2015.1.01015.S programme (principal investigator A. Bayo). The present study made use of detection limits derived from SPHERE data collected during the programme 198.C-0209, and of SPHERE data collected during the programmes 097.C-0319(A), 105.209E.001 and 198.C-0209(F). The data are available athttps://archive.eso.org.

JWST MIRI data reduction and analysis codes are available via GitHub athttps://github.com/mathildemalin/ExoCAT. Data reduction was carried out using the publicly available JWST pipelinehttps://jwst-pipeline.readthedocs.io/en/latest/, and the stpf packagehttps://stpsf.readthedocs.io/en/latest/to simulate PSFs for JWST. We used various functions of the following software packages to perform the analysis and create the figures: numpy, astropy, scipy, matplotlib and photutils.

Backman, D. E. & Paresce, F. inProtostars and Planets IV639 (1993).

Wyatt, M. C. Evolution of debris disks.Annu. Rev. Astron. Astrophys.46, 339–383 (2008).

Matrà, L. et al. REsolved ALMA and SMA Observations of Nearby Stars (REASONS): A population of 74 resolved planetesimal belts at millimetre wavelengths.Astron. Astrophys.693, A151 (2025).

Pawellek, N. & Krivov, A.Resolved Debris Discs(Astrophysikalisches Institut und Universitäts-Sternwarte);https://www.astro.uni-jena.de/index.php/theory/catalog-of-resolved-debris-disks.html.

Shannon, A. et al. The unseen planets of double belt debris disc systems.Mon. Not. R. Astron. Soc.462, L116–L120 (2016).

Su, K. Y. L. et al. Asteroid belts in debris disk twins: Vega and Fomalhaut.Astrophys. J.763, 118 (2013).

Milli, J. et al. Near-infrared scattered light properties of the HR 4796 A dust ring. A measured scattering phase.Astron. Astrophys.599, A108 (2017).

Pearce, T. et al. Planet populations inferred from debris discs. Insights from 178 debris systems in the ISPY, LEECH, and LIStEN planet-hunting surveys.Astron. Astrophys.659, A135 (2022).

Faramaz, V. et al. A detailed characterization of HR 8799’s debris disk with ALMA in band 7.Astron. J.161, 271 (2021).

Gardner, J. P. et al. The James Webb Space Telescope Mission.Publ. Astron. Soc. Pac.135, 068001 (2023).

Wright, et al. The Mid-Infrared Instrument for the James Webb Space Telescope, II: design and build.Publ. Astron. Soc. Pac.125, 595 (2015).

Binks, A. S. et al. A kinematically hot population of young stars in the solar neighbourhood Mon.Not. R. Astron. Soc.494, 2429–2439 (2020).

Gaia Collaboration et al. Gaia data release 3. Summary of the content and survey properties.Astron. Astrophys.674, A1 (2023).

Stassun, K. G. et al. The TESS input catalog and candidate target list.Astron. J.156, 102 (2018).

Nicholson, B. A. et al. The surface magnetic activity of the weak-line T Tauri stars TWA 7 and TWA 25.Mon. Not. R. Astron. Soc.504, 2461–2473 (2021).

Choquet, E. et al. First images of debris disks around TWA 7, TWA 25, HD 35650, and HD 377.Astrophys. J.817, L2 (2016).

Olofsson, J. et al. Resolving faint structures in the debris disk around TWA 7. Tentative detections of an outer belt, a spiral arm, and a dusty cloud.Astron. Astrophys.617, A109 (2018).

Ren, B. et al. A layered debris disk around M Star TWA 7 in scattered light.Astrophys. J.914, 95 (2021).

Rieke, et al. The Mid-Infrared Instrument for the James Webb Space Telescope, VII: the MIRI detectors.Publ. Astron. Soc. Pac.127, 665 (2015).

VizieR Online Data Catalog: Gaia EDR3(Gaia Collaboration, 2020).

Bayo, A. et al. Sub-millimetre non-contaminated detection of the disc around TWA 7 by ALMA.Mon. Not. R. Astron. Soc.486, 5552–5557 (2019).

Boccaletti, A. et al. Imaging detection of the inner dust belt and the four exoplanets in the HR 8799 system with JWST’s MIRI coronagraph.Astron. Astrophys.686, A33 (2024).

Mahlke, M. et al. Mining the Kilo-Degree Survey for solar system objects.Astron. Astrophys.610, A21 (2018).

Brown, M. E. et al. Discovery of a planetary-sized object in the Scattered Kuiper Belt.Astrophys. J.635, L97–L100 (2005).

Brown, M. E. et al. Discovery of a candidate inner Oort Cloud planetoid.Astrophys. J.617, 645 (2004).

Batygin, K. et al. Evidence for a distant giant planet in the Solar System.Astron. J.151, 22 (2016).

Fienga, A. et al. New constraints on the location of P9 obtained with the INPOP19a planetary ephemeris.Astron. Astrophys.640, A6 (2020).

Brown, M. E. & Batygin, K. The orbit of Planet Nine.Astron. J.162, 219 (2021).

Wilkinson, C. Breaking degeneracies in exoplanetary parameters through self-consistent atmosphere-interior modelling.Astron. Astrophys.692, A113 (2024).

Linder, E. et al. Evolutionary models of cold and low-mass planets: cooling curves, magnitudes, and detectability.Astron. Astrophys.623, A85 (2019).

Ozernoy, L. M. et al. Signatures of extrasolar planets in dust debris disks.Astrophys. J.537, L147 (2000).

Stark, C. C. & Kuchner, M. J. The detectability of exo-Earths and super-Earths via resonant signatures in exozodiacal clouds.Astrophys. J.686, 637 (2008).

Levison, H. F. & Duncan, M. J. The long-term dynamical behavior of short-period comets.Icarus108, L18 (1994).

Lagrange, A.-M. et al. A giant planet imaged in the disk of the young star beta Pictoris.Science329, 57–59 (2010).

Matthews, E. C. et al. A temperate super-Jupiter imaged with JWST in the mid-infrared.Nature633, 789–792 (2024).

ArticleCASPubMedPubMed CentralGoogle Scholar

LaJoie, C.-P. et al. Small-grid dithers for the JWST coronagraphs.SPIE Conf. Ser.9904, 99045K (2016).

Boccaletti, A. et al. JWST/MIRI coronagraphic performances as measured on-sky.Astron. Astrophys.667, A165 (2022).

Mâlin, M. et al. Unveiling the HD 95086 system at mid-infrared wavelengths with JWST/MIRI.Astron. Astrophys.690, 316.

Boccaletti, A. et al. The four-quadrant phase mask coronagraph. IV. First light at the Very Large Telescope.Astron. Astrophys.116, 825 (2005).

Boccaletti, A. et al. The Mid-Infrared Instrument for the James Webb Space Telescope, V: predicted performance of the MIRI coronagraphs.Publ. Astron. Soc. Pac.127, 633 (2015).

Perrin, M. D. et al. Updated point spread function simulations for JWST with WebbPS.Proc. SPIE9143, 91433X (2014).

Nelder, J. A. & Mead, R. A simplex method for function minimization.Comput. J.7, 308–313 (1995).

ArticleMathSciNetGoogle Scholar

Salo, H. et al. The Spitzer Survey of Stellar Structure in Galaxies (S4G): multi-component decomposition strategies and data release.Astrophys. J. Supp. Ser.219, 4 (2015).

Combes, F. Molecular gas in distant galaxies from ALMA studies.Astron. Astrophys. Rev.26, 5 (2018).

Ling, C.-T. et al. Galaxy source counts at 7.7, 10, and 15 μm with the James Webb Space Telescope.Mon. Not. R. Astron. Soc.517, 85 (2022).

Wu, C. K. W. et al. Source counts at 7.7-21 μm in CEERS field with JWST.Mon. Not. R. Astron. Soc.523, 5187–5197 (2023).

Stone, M. A. et al. 5-25 μm galaxy number counts from deep JWST data.Astrophys. J.972, 62 (2024).

Sajkov, L. et al. Halfway to the peak: The JWST MIRI 5.6 μm number counts and source population.Astrophys. J.977, 115 (2024).

Polletta, M. et al. Spectral energy distributions of hard X-ray selected active galactic nuclei in the XMM-Newton Medium Deep Survey.Astrophys. J.663, 81–102 (2007).

Caputi, K. et al. The infrared luminosity function of galaxies at redshifts z=1 and z~2 in the GOODS fields.Astrophys. J.660, 97 (2007).

Jarrett, T. H. et al. The Spitzer-WISE Survey of the ecliptic poles.Astrophys. J.735, 112 (2011).

Dale, D. et al. A two-parameter model for the infrared/submillimeter/radio spectral energy distributions of galaxies and active galactic nuclei.Astrophys. J.784, 83 (2014).

Bernard-Salas, J. et al. A Spitzer high-resolution mid-infrared spectral atlas of starburst galaxies.Astrophys. J. Supp. Ser.184, 230 (2009).

This work is based on observations made with the JWST of NASA (National Aeronautics and Space Administration), the European Space Agency and the Canadian Space Agency. The data were obtained from the Mikulski Archive for Space Telescopes at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-03127 for JWST. These observations are associated with programme number 3662. This paper makes use of ALMA data. ALMA is a partnership of the European Organisation for Astronomical Research in the Southern Hemisphere (representing its member states), the National Science Foundation (USA) and the National Institutes of Natural Sciences (Japan), together with the National Research Council of Canada (Canada), the National Science Council and the Academia Sinica Institute of Astronomy and Astrophysics (Taiwan) and the Korea Astronomy and Space Science Institute (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by the European Organisation for Astronomical Research in the Southern Hemisphere, the National Radio Astronomy Observatory operated by Associated Universities, Inc., and the National Astronomical Observatory of Japan. This paper makes use of observations collected at the European Organisation for Astronomical Research in the Southern Hemisphere. This project has received funding from the European Research Council under the European Union’s Horizon 2020 research and innovation programme (COBREX; grant agreement number 885593). L.M. acknowledges funding by the European Union through the E-BEANS ERC project (grant number 100117693), and by the Irish Research Council under grant number IRCLA-2022-3788. Views and opinions expressed are those of the authors alone and do not necessarily reflect those of the European Union or the ERC Executive Agency. Neither the European Union nor the granting authority can be held responsible for them.

LIRA, Observatoire de Paris, Université PSL, Sorbonne Université, Université Paris Cité, CY Cergy Paris Université, CNRS, Meudon, France

A.-M. Lagrange, C. Wilkinson, A. Boccaletti, C. Perrot, D. Rouan, A. Chomez, B. Charnay, Q. Kral, P. Thebault, F. Kiefer, A. Radcliffe, J. Mazoyer & S. Stasevic

Université Grenoble Alpes, CNRS, IPAG, Grenoble, France

A.-M. Lagrange, H. Beust, A. Chomez, J. Milli & P. Delorme

Department of Physics & Astronomy, Johns Hopkins University, Baltimore, MD, USA

Space Telescope Science Institute, Baltimore, MD, USA

M. Mâlin, A. Carter & K. A. Crotts

School of Physics, Trinity College Dublin, University of Dublin, Dublin, Ireland

Observatoire de Paris, LUX, PSL University, Collège de France, Sorbonne University, CNRS, Paris, France

Observatoire de la Côte d’Azur, Université Côte d’Azur, Nice, France

Centre de Recherche Astrophysique de Lyon, Université Claude Bernard Lyon 1, CNRS, ENS, Saint-Genis-Laval, France

European Southern Observatory, Garching bei Muenchen, Germany

J. Olofsson, A. Bayo & M. Langlois

Département d’Informatique de l’École Normale Supérieure (ENS-PSL, CNRS, Inria), Paris, France

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

A.-M.L. led the proposal, the interpretation of the various data and the writing of the paper. A. Boccaletti, M.M. and C.P. designed the JWST MIRI observing sequence, and reduced the data and participated in their interpretation. C.W. conducted the SED fitting. B.C. participated in the vetting of the Solar System origin of the source. A.R., B.C. and S.M. participated in discussions on the SED fitting. A. Chomez provided the SPHERE detection limit maps used in this paper. A. Chomez, O.F., M.L., P.D. and G.C. participated in their interpretation. S.S. and T.B. tried to improve the detection limits using alternative algorithms. J. Milli participated in the building of the proposal, produced the SPHERE polarimetric images and participated the interpretation of the JWST MIRI data. A. Bayo and J.O. provided their analysis of the ALMA data, and participated in their interpretation. L.M. reanalysed the ALMA data and led the interpretation of these data. D.R. and F.C. conducted the estimation of the galaxy background probability. H.B. performed theN-body dynamical simulations. Q.K. and P.T. participated in the writing of the proposal and in the discussion of the disk–planet interactions. J. Mazoyer participated in the discussion of the disk properties. F.K. participated in the analysis of Gaia data, and general discussion. A. Carter and K.C. participated in the discussion on the astrometry of the source. P.R. participated in the proposal writing and the analysis of the data. All authors commented on the manuscript.

Correspondence toA.-M. Lagrange.

The authors declare no competing interests.

Naturethanks the anonymous reviewers for their contribution to the peer review of this work.

Publisher’s noteSpringer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Orientation of the 4QPM and its phase transitions in the sky plane for the two telescope rolls. North is up, East is left.

The curves are distinguished by their colors, and the corresponding labels, valid for all panels. All curves are calibrated to have a flux of 22 mJy at 11 μm. The two additional constraints are marked by black symbols (square and triangle): the flux should be lower than 60 mJy at 1.6 μm, and 96 mJy at 870 μm. The SED which are not consistent with the limits are plotted in dashed lines.

Corner plot of cloudy forward modeling using SPHERE upper limits and age constraints of 6.4+/− 1 Myr. The associated priors are listed in Extended Data Table3. The mass found is 0.34+/− 0.06 MJ(considering errors from the MCMC only). Core given in Earth mass, Tintcorresponds to the intrinsic temperature and fsedthe sedimentation rate of the considered clouds.

Thermal evolution curves from ref.30, showing the effective temperature evolution of cloudy planets with an [Fe/H] of 0.4 dex. The estimated age and effective temperature of TWA 7b supposing the planet and star are coeval is represented by the blue point. The effective temperature used is 316+19/−23K (derived from the forward modeling) and the age 6.4+/− 1 Myr.

Combined ALMA 0.88 mm image of the TWA 7 system obtained with Briggs 0.5 weighting, centered at the phase center of the April compact configuration observations. The background galaxy is clearly detected East of the expected stellar position, at a position consistent with the East source detected by MIRI (black and white circle). The stellar location at each of the 2016 ALMA epochs is shown by the green star (positions largely overlapping), whereas the stellar location at the 2024 MIRI observation is shown by the orange star. The position of the CC#1 source at the epoch of the 2024 MIRI observation is shown by the cyan circle. The image (not primary beam-corrected) has a resolution of 0.19” x 0.18” (shown as the circle in the bottom left of the image) and an RMS noise level of 23 μJy/beam.

Contrast 5-sigma confidence level curves of the SPHERE observation used to compute the upper limits on the candidate companion flux.

Supplementary Notes 1 and 2, Tables 1 and 2 and References.

Open AccessThis article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visithttp://creativecommons.org/licenses/by-nc-nd/4.0/.

Lagrange, AM., Wilkinson, C., Mâlin, M.et al.Evidence for a sub-Jovian planet in the young TWA 7 disk .Nature(2025). https://doi.org/10.1038/s41586-025-09150-4

DOI:https://doi.org/10.1038/s41586-025-09150-4

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

Nerve-to-cancer transfer of mitochondria during cancer metastasis

The nervous system has a pivotal role in cancer biology, and pathological investigations have linked intratumoural nerve density to metastasis1. However, the precise impact of cancer-associated neurons and the communication channels at the nerve–cancer interface remain poorly understood. Previous cancer denervation models in rodents and humans have highlighted robust cancer dependency on nerves, but the underlying mechanisms that drive nerve-mediated cancer aggressivity remain unknown2,3. Here we show that cancer-associated neurons enhance cancer metabolic plasticity by transferring mitochondria to cancer cells. Breast cancer denervation and nerve–cancer coculture models confirmed that neurons significantly improve tumour energetics. Neurons cocultured with cancer cells undergo metabolic reprogramming, resulting in increased mitochondrial mass and subsequent transfer of mitochondria to adjacent cancer cells. To precisely track the fate of recipient cells, we developed MitoTRACER, a reporter of cell-to-cell mitochondrial transfer that permanently labels recipient cancer cells and their progeny. Lineage tracing and fate mapping of cancer cells acquiring neuronal mitochondria in primary tumours revealed their selective enrichment at metastatic sites following dissemination. Collectively, our data highlight the enhanced metastatic capabilities of cancer cells that receive mitochondria from neurons in primary tumours, shedding new light on how the nervous system supports cancer metabolism and metastatic dissemination.

Cancer plasticity has a pivotal role in cancer heterogeneity and the emergence of adaptive phenotypes4. Although cancer cells often rely on glycolysis for energy, numerous studies indicate that they adapt their metabolism, including oxidative phosphorylation (OXPHOS), to changing environments. The metabolic plasticity of cancer cells is becoming increasingly recognized as essential for cancer progression by enabling cancer cells to efficiently navigate through the metastatic cascade5,6. Studies on cancer metabolic plasticity have primarily focused on cell-autonomous mechanisms, in which cells modify their metabolic program and perform metabolic rewiring by regulating metabolic enzymes and adapting alternative substrates. The non-autonomous mechanisms of metabolic plasticity are less well understood and can involve complex interactions within the cancer microenvironment7. For instance, stromal cells, including fibroblasts and immune cells, can share a variety of metabolites, growth factors and cytokines, contributing to the metabolic reprogramming of cancer8. Overall, non-cell-autonomous mechanisms of the cancer cells’ metabolic plasticity are still unclear and may represent promising targets for preventing metastatic dissemination.

Many cancers are innervated, and in recent years, cancer neurobiology has garnered considerable attention1. Studies scrutinizing the impact of the nervous system in various cancers have established a critical role of cancer-infiltrating neurons in controlling cancer aggressivity, including supporting cancer cell metabolism through direct metabolic interactions1,9. Pathological analyses in patients with cancer consistently associated cancer innervation with negative clinical outcomes, and targeted ablation of intratumoural nerves can suppress cancer growth in a tissue-specific manner1,10,11,12. The neuronal component of the cancer microenvironment can arise from the recruitment of pre-existing nerves through perineural invasion or the de novo generation of nerves within the cancer stroma through cancer neurogenesis13,14. Nerve withdrawal in human prostate cancer demonstrated impaired tumour growth and metabolism, and metabolomic analysis of denervated cancers across rodent and human species revealed disrupted metabolic efficiency that is characterized by their downregulated mitochondrial metabolism and transition of the cancer cells into a more glycolysis-dependent production of energy2. Overall, converging evidence indicates that denervated cancer actively suffers from nerve withdrawal, with a predominant effect on their energetic metabolism1,15. These studies demonstrated the cancer metabolic dependency on nerves and suggested the existence of metabolic support mechanisms acting at the nerve–cancer interface. However, their nature remains unknown, and characterizing these mechanisms, which represent a critical gap in the understanding of the metabolic support provided by nerves to cancer cells, is the focus of this study.

Here we demonstrate the metabolic dependency of breast cancer cells on nerves. Using orthogonal in vitro and in vivo methods, we observed the prevalence of transfers of mitochondria between neurons and cancer cells. To trace the fate of the recipient cells, we developed MitoTRACER, a genetic reporter that permanently marks cancer cells receiving mitochondria from donor cells. Our approach revealed that neuron-derived mitochondria enhance cancer cell metabolic capacity, stemness and resistance to metastatic stressors. Fate mapping of the recipient cells in vivo has shown their increased metastatic capacities. In the human context, multispectral imaging with machine-learning deconvolution showed increased mitochondrial mass in metastatic cells, and perineural invasion was associated with higher mitochondrial content in cancer cells near nerves. Finally, cancer chemically denervated with botulinum neurotoxin type A (BoNT/A) confirmed reduced mitochondrial load in cancer cells. These findings highlight nerve-driven metabolic support as essential to cancer metabolic plasticity and metastatic potential.

Previous studies on the effects of BoNT/A-mediated prostate denervation conducted in humans and rats uncovered the cancer dependency on nerves and the metabolic reprogramming of prostate cancer cells after nerve withdrawal2. Similarly, breast cancer denervation resulted in a marked decrease in cancer growth15. However, the impact of cancer-infiltrating nerves in breast cancer metabolism has not been explored. We investigated the consequences of nerve withdrawal in breast cancer using two complementary mouse models. We first adapted a BoNT/A-mediated denervation protocol to the 4T1 model of triple-negative breast cancer (TNBC; Fig.1a,b) and to a human ductal carcinoma in situ xenograft model16(Extended Data Fig.1a–d). In both settings, transcriptomic profiling of cancer cells from control and denervated tumours revealed distinct expression signatures following denervation (Fig.1band Extended Data Fig.1a,b), and Gene Ontology analyses have underscored the prevalent downregulation of metabolic processes (Extended Data Fig.1b–dand Supplementary Tables1–3). Gene set enrichment analysis of the ductal carcinoma in situ model identified the tricarboxylic acid cycle as the most suppressed pathway overall (Extended Data Fig.1cand Supplementary Table2). Histopathological examination of this model has shown a reduction in the incidence of invasive lesions from 55% in control mice to 12% in denervated mice (Extended Data Fig.1e), underscoring the functional importance of innervation for breast cancer progression.

a,b, Pre-denervation breast cancer model using BoNT/A injections in BALB/c mice followed by implantation of 4T1mCherrycells. Transcriptomic analysis of 4T1mCherrycells revealed a distinct transcriptomic signature and downregulation of metabolic processes. SSC-H, side scatter height; RNA-seq, RNA sequencing.c, Confocal micrograph of SVZ-NSCsGFPmixed with 4T1mCherrycells. White arrows show the establishment of neuron–cancer contacts. Scale bar, 50 μm.d, 4T1 cells FACS-isolated from coculture (left) showed increased OXPHOS capacities (Seahorse assay; mean ± s.e.m; representative profile (n= 3); right). O, oligomycin; F, carbonyl cyanide-p-trifluoromethoxyphenylhydrazone; R + A, rotenone and antimycin A; OCR, oxygen consumption rate.e,f, Time-lapse confocal microscopy (e) and flow cytometry (f) highlighted mitochondrial transfer from SVZ-NSCsCCO-GFPto 4T1mCherrycells (n= 6 independent cocultures). Scale bar, 10 μm. The white arrows show transferred mitochondria.g,h, 3D reconstruction shows transfer through tunnelling nanotubes between peripheral nervous system-derived 50B11-DRGCCO-GFPcells and 4T1mCherrycells. Red arrow shows the tunnelling nanotube structure, and white arrows show transferred mitochondria. Scale bars, 20 μm.i, Quantification of direct cell–cell contact (23.04%) versus distant (0.59%) transfers using Transwell inserts. Normalized transfer rate is calculated as the percentage of 4T1mCherry+/GFP+cells among the eGFP+cells in the coculture. Mean ± s.d., Student’s two-tailed unpairedt-test; ****P< 0.0001 (n= 6 independent cocultures). Ctrl, control.j, Cytochalasin B (Cyto B) reduces mitochondrial transfer. DMSO, dimethylsulfoxide. Mean ± s.d., Student’s two-tailed unpairedt-test, ***P= 0.001 (n= 3 independent cocultures).k, Mitochondrial transfer rates vary with donor cells from different origins. MEF, mouse embryonic fibroblast. Mean ± s.d. (n= 5 independent cocultures).l, 4T1mCherrycells rendered devoid of mtDNA (ρ0) were cocultured with ρ+SVZ-NSCsGFPand isolated by FACS at various times to monitor the transfers. mUNG1, Y147A mutant of uracil-N-glycosylase.m, PCR of mtDNA content in FACS-sorted ρ04T1mCherrycells showed progressive reacquisition of mtDNA. Expression of GFP and mCherry confirmed the purity of the 4T1 cells. Nuclear DNA (nDNA) was used as the loading control.n, MitoTracker microscopy imaging of ρ+4T1, ρ04T1 and ρ04T1 cells FACS-sorted from coculture shows rescue in their mitochondrial morphology.o, Reacquisition of SVZ-derived mtDNA in ρ04T1 cells restores uridine-independent growth.p,q, FACS-isolated ρ04T1 cells rescued by mitochondrial transfer regained OXPHOS (p; Seahorse assay, mean ± s.e.m.; representative profile (n= 3)) and proliferative capacities (q; direct cell counting, mean ± s.d.,n= 6 independent cultures); two-way analysis of variance (ANOVA), ****P< 0.0001; NS, not significant.a, Created in BioRender. S. Grelet (2025)https://biorender.com/0j8zovf.d, Created in BioRender. S. Grelet (2025)https://biorender.com/oxxilqq.i, Created in BioRender. S. Grelet (2025)https://biorender.com/8tdz09x.l, Created in BioRender. S. Grelet (2025)https://biorender.com/pm5yh64.

To explore in more detail the mechanisms associated with cancer metabolic dependencies on nerves in breast cancer progression, we developed an in vitro nerve–cancer coculture. We mixed the 4T1 aggressive mouse breast carcinoma cells with neuronal stem cells (NSCs) from the mouse subventricular zone (SVZ) (Extended Data Fig.2a–d), as they have been proved to be a source of cancer innervation14. After stimulation by 4T1 cancer cells, SVZ-NSCs rapidly undergo differentiation, evidenced by their morphological transition from round, neural progenitor cells to more elongated cells with neurite extensions and forming close contact with the cancer cell (Fig.1cand Extended Data Fig.2e). SVZ-NSCs can differentiate into neurons, astrocytes or oligodendrocyte glial cells. Using our coculture model, we confirmed the neuronal commitment of NSCs, as shown by microscopy through the expression of the neuronal marker tubulin β3 (TUBB3) (Extended Data Fig.2f). Absolute quantification by flow cytometry showed that more than 90% of progenitor cells differentiated into TUBB3+and MAP2+neurons but remained negative for the glial cell markers O4 and ALDH1L1 (Extended Data Fig.2g–i). SVZ-NSCs also exhibited neuronal functional characteristics, including calcium pulsatile activity (Extended Data Fig.2jand Supplementary Video2) and the ability to generate action potentials in response to depolarizing current injection, with an action potential threshold averaging −46 mV ± 5 mV as observed through whole-cell electrophysiology recordings (Extended Data Fig.2k,l).

After establishing the biological relevance of the nerve–cancer coculture in vitro, we next confirmed that the cancer’s metabolic dependency on nerves, observed in mouse xenografts, could also be replicated in a more simplistic in vitro model. Cancer cells in monoculture or cultivated in the presence of neurons were isolated before their metabolic analysis. In comparison to those cultivated alone, the cancer cells isolated from the nerve–cancer coculture exhibited upregulated mitochondrial respiration (Fig.1d), characterized by significant increases in their basal and maximal mitochondrial respiration as well as enhanced spare respiratory capacities (Extended Data Fig.3a). This reflects the establishment of the nerve–cancer metabolic dependencies in the coculture model in vitro.

As previously described, cancer-induced neuronal progenitor differentiation is integral to the nerve–cancer cross-talk established during cancer innervation1. For instance, we and others have previously demonstrated how the cancer-driven expression of axon guidance molecules such as semaphorin 4F controls the cancer-induced neuronal precursor differentiation into the cancer stroma and increases intratumoural nerve density and cancer aggressivity17,18,19. Physiologically, differentiating neuronal progenitors must undergo a metabolic shift to increase their mitochondrial metabolism and meet their specific energetic demands during physiological differentiation20(Extended Data Fig.3b,c). We examined whether similar reprogramming also occurs during cancer-driven neuronal differentiation and confirmed a robust increase in the mitochondrial mass of SVZ-NSCs exposed to 4T1 cancer cells (Extended Data Fig.3d), with absolute quantification of mitochondrial DNA (mtDNA) showing an increasing mtDNA load from about 16 to 226 mtDNA/nuclear DNA copies per neuron after their cancer-induced differentiation (Extended Data Fig.3e). Finally, genetic fluorescent labelling of SVZ-NSC mitochondria confirmed morphological changes of the SVZ-NSCs exposed to the cancer cells and the development of an extended mitochondrial network across the nerve–cancer coculture, transitioning from globular structures in NSC monoculture to thin and elongated tubular structures extending throughout the NSCs in coculture with breast cancer cells (Extended Data Fig.3f). Such changes represent typical hallmarks of metabolic reprogramming, in which maturing neurons transition from glycolytic to mitochondrial oxidative metabolism20,21.

The establishment of close nerve–cancer cross-talk in our in vitro cocultures, the significant increase in neuronal mitochondria abundance and the enhancement of mitochondrial metabolism in isolated cancer cells following exposure to neuronal cells in coculture suggested that metabolism-related collaboration mechanisms initiated at the nerve–cancer interface. The cell–cell transfer of mitochondria has recently become a subject of investigation as it has crucial functions in health and disease22,23,24, and was recently shown to have profound effects on cancer progression25,26. Moreover, research has demonstrated that astrocytes can provide metabolic support to glioblastoma cells by mitochondria transfer in the central nervous system27.

We therefore investigated whether mitochondrial transfer could occur within the peripheral nervous system and assessed whether breast cancer cells exposed to neurons could acquire neuron-derived mitochondria. To test this hypothesis, we established nerve–cancer cocultures by combining central nervous system-derived neurons (SVZ-NSCs) or dorsal root ganglia-derived neurons (50B11-DRG) with 4T1 breast cancer cells expressing the mCherry fluorophore (4T1mCherry+). Neurons were genetically modified to express enhanced green fluorescent protein (eGFP)-labelled mitochondria (SVZ-NSCCCO-GFPand 50B11-DRGCCO-GFP) (Fig.1e–hand Supplementary Video1). Confocal microscopy and flow cytometry analysis of the coculture confirmed mitochondrial transfer from neurons to cancer cells (Fig.1e–hand Extended Data Fig.4a) and the formation of tunnelling nanotube-like structures, facilitating organelle transfer (Fig.1g, red arrow), as confirmed by three-dimensional (3D) reconstruction (Fig.1h, white arrows). Flow cytometry (Fig.1fand Extended Data Fig.4a) showed acquisition of a double-positive 4T1GFP+/mCherry+subpopulation of cells in the cocultures that reflects the acquisition of eGFP-labelled mitochondria into the recipient 4T1mCherry+cells. Although this population accounts for an average of 0.96% of the total cells in the coculture (Fig.1f), it represents only a snapshot of the current mitochondria transfers, and the SVZ-NSCs providing the mitochondria represent, on average, only 3.06% of the coculture. When normalized to the population of eGFP-labelled cells in the coculture, the double-positive (4T1mCherry+/GFP+) fraction accounted for an equivalent of 31.4% of the donor population. This normalization method was used to standardize the transfer rate from donor to recipient cells for different experimental conditions. Beyond the direct cell–cell contact-mediated transfer of mitochondria, we investigated whether distant mechanisms, such as microvesicles, could be involved. Distant coculture using Transwell inserts confirmed that cell–cell contact is the primary route of mitochondrial transfer, although distant transfers also occurred (Fig.1i). Inhibition of tunnelling nanotube formation further validated the role of these structures in the transfer process without a significant impact on cell viability (Fig.1jand Extended Data Fig.4b), as previously shown27. Using different cell lines, we next tested the capacities of cells with diverse origins (50B11-DRG neurons, SVZ neurons, PC12 pheochromocytoma, HT-22 hippocampal, Neuro2A neuroblastoma, 3T3-L1 pre-adipocytes, mouse embryonic fibroblasts, NMuMG normal mouse mammary gland and 4T1 breast cancer cells) in transferring their mitochondria to cancer cells (Fig.1k). Every tested donor significantly transferred mitochondria to the recipient 4T1 cells. Notably, cell lines of neuronal origin exhibited higher mitochondrial transfer rates. In addition, we confirmed mitochondria transfer in nerve–cancer cocultures using human cancer cell lines (Extended Data Fig.4c).

We subsequently validated the transfer of neuronal mitochondria to cancer cells through an additional orthogonal validation approach. We generated rho-zero (ρ0) 4T1 cancer cells lacking mtDNA and cocultured them with rho-plus (ρ+) SVZ-NSCs28(Fig.1l). PCR amplification of mtDNA from the ρ04T1mCherrycells confirmed a complete loss of mtDNA and revealed their gradual reacquisition of mtDNA through coculture with SVZ-NSCs (Fig.1mand Extended Data Fig.4d). The mitochondrial morphology was altered in ρ0cells, with the mitochondria exhibiting a globular and fragmented appearance as described previously29. This morphology was restored following the mtDNA restoration through nerve–cancer coculture (Fig.1n). ρ0cells are also phenotypically characterized by the complete loss of OXPHOS capacities and auxotrophy for uridine30,31. Therefore, ρ0cells require uridine complementation to grow in culture in vitro. We tested whether transferred mitochondria are functional in recipient cells. We grew ρ04T1 cells in the presence of uridine either alone or in coculture with ρ+SVZ-NSCs. After 5 days of coculture, 4T1mCherry+cells were sorted by fluorescence-activated cell sorting (FACS) and subcultured without uridine. In the absence of uridine, ρ04T1 cells did not form colonies, but a subset of the ρ04T1 cells enriched from coculture with ρ+neurons formed viable colonies (Fig.1o), reflecting rescue of uridine synthesis, therefore demonstrating the functionality of the transferred mitochondria. The SVZ-NSC-derived mitochondria transferred to the ρ04T1 cells also rescued their mitochondrial respiration and proliferative capacities as observed by Seahorse extracellular flux metabolic analysis (Fig.1p, ρ0+ SVZ) and proliferation assays, respectively (Fig.1q, ρ0+ SVZ).

We next tested the occurrence of nerve–cancer transfers of mitochondria in vivo. Clinically, human prostate cancer samples showed an increased mitochondrial load in cancer cells associated with perineural invasion (Fig.2a). High-throughput multispectral imaging quantification based on machine learning32highlighted that cancer cells closer to nerves had a significantly higher mitochondrial load than cells farther from nerves (Fig.2band Supplementary Tables4and5). To assess whether this increased mitochondrial load may result from mitochondrial transfer, we analysed prostate tissue from a clinical trial (NCT01520441), in which human prostate cancer was chemically denervated. Multispectral imaging revealed that cancer cells of the BoNT/A-denervated side of the prostate had a lower mitochondrial load than the saline-injected side, supporting the hypothesis that nerves promote mitochondrial load in cancer cells, with mitochondrial transfer as a contributing factor (Fig.2cand Supplementary Tables6and7).

a, Histopathology of human prostate cancer with perineural invasion shows increased mitochondrial content near nerves (mitochondria are visualized by periodic acid–Schiff staining (magenta); nerves are visualized by diaminobenzidine staining (brown)). Representative profile (n= 72 patients). Scale bar, 180 μm.b, Multispectral imaging combined with machine learning-based image deconvolution, spatial correlation, and quantification indicate significantly higher mitochondrial loads in prostate cancer cells near nerves (perineural;n= 72, 40,007 cells) compared to distant cancer cells (distant;n= 58, 20,766 cells). The line shows the median, the box boundaries show the 25th and 75th percentiles, and the whiskers show the minimum and maximum values. Two-sided Welch’st-test, ***P< 0.001.c, BoNT/A-mediated denervation reduces mitochondrial load in human prostate cancer cells (paired analysis: saline versus BoNT/A; saline:n= 10,918 cells, BoNT/A:n= 14,186 cells). Two-sided Welch’st-test, ****P= 1.463 × 10−113. Clinical trial (NCT01520441).d, Mouse DRG neurons innervating the mammary gland were labelled with lentivirus (LV) to tag neuronal mitochondria before injection of 4T1mCherrycells into mammary fat pads (MFPs). Cancer cells were isolated post-tumour growth for mitochondrial transfer analysis to detect host-derived mitochondrial transfer.e, Neuronal mitochondria were labelled using lentiviruses expressing either nuclear-localized (GFP-NLS, non-transferable) or mitochondria-localized (GFP-OMP25, transferable) eGFP under the synaptin1 (Syn1) promoter. Flow cytometry of cancer cells identified eGFP+subpopulations, indicating neuronal mitochondrial transfer between mouse host neurons and cancer xenografts.f, Sanger sequencing enabled the detection of mtDNA polymorphisms between host BALB/c cell and 4T1 cancer cell mtDNA.g, Oxford Nanopore sequencing analysis of mtDNA heteroplasmy in FACS-isolated cancer cells demonstrated host-to-cancer mitochondrial transfer. BoNT/A-mediated pre-denervation at the xenograft site significantly reduced mitochondrial transfer to cancer cells (salinen= 8, BoNT/An= 9; median values indicated). The line shows the median, the box boundaries show the 25th and 75th percentiles, and the whiskers show the minimum and maximum values. One-tailed unpaired Student’st-test *P= 0.0316 (n= 8 saline,n= 9 BoNT/A mice).c, Created in BioRender. S. Grelet (2025)https://biorender.com/mprd95w.d, Created in BioRender. S. Grelet (2025)https://biorender.com/98cu18d.g, Created in BioRender. S. Grelet (2025)https://biorender.com/xapzg43.

Using a BALB/c mouse xenograft model, we evaluated mitochondrial transfer from mouse host neurons to 4T1 breast cancer cells in vivo (Fig.2d). We designed a lentiviral construct encoding a neuron-specific, Syn1-GFP-OMP25, mitochondria-anchored eGFP reporter to label the mitochondria of host mammary fat pad neurons (Fig.2eand Extended Data Fig.5a–c). After genetic modification of the mouse DRG innervating the lower mammary fat pads, 4T1mCherry+cancer cells were injected into the corresponding fat pads, and flow cytometry analysis of the primary cancer showed a subpopulation of cancer cells exhibiting the green signal of eGFP, consistent with mitochondrial transfer from mouse neurons to cancer cells in situ (Fig.2e). We developed a similar lentiviral construct as a control, in which eGFP was targeted to the nucleus (GFP-NLS); this construct was non-transferable and showed no signal transfer between mouse nerves and 4T1mCherry+cells (Fig.2e).

We further assessed mitochondrial transfer by analysing mtDNA heteroplasmy arising from a mtDNA polymorphism identified between the mouse host and cancer cells, as identified by Sanger sequencing (Fig.2f), with this analysis serving as an additional independent approach. 4T1mCherry+cells from cancer xenografts were FACS-sorted and analysed for their mtDNA content through Nanopore sequencing, which confirmed the overall acquisition of mouse-derived mtDNA by cancer cells (Fig.2g). The neuronal origin of these transfers was further validated through BoNT/A-mediated chemical denervation, which demonstrated that neurons accounted for approximately 35% of the total mitochondrial transfers between host neurons and cancer cells (Fig.2g).

As the previous results established the biological relevance of nerve-to-cancer mitochondrial transfer in vitro and in vivo, and suggested its clinical relevance, it was important to investigate the effects of mitochondrial transfer on cancer cell biology. However, the ρ04T1 cell model (Fig.1l–q) was limited in biological relevance for studying the functional consequences of mitochondrial transfer in recipient cells. In addition, the mitochondrial-bound methods used in vitro (Fig.1e–k) and in vivo (Fig.2d,e) were limited to mitochondrial transfer at the time of analysis and could not evaluate previous transfers. Furthermore, after an eGFP-labelled mitochondrion entered a recipient cell, the eGFP signal faded quickly because the recipient cells did not express the mitochondrial genetic reporter. Therefore, these methods do not allow for distinguishing between cells in the same culture that received mitochondria versus those that did not, nor do they enable lineage tracing of recipient cells and their progeny. To address this limitation, we designed a new genetic reporter strategy known as MitoTRACER to permanently label the recipient cancer cells after mitochondrial transfer from neurons and distinguish them from cells that did not receive mitochondria.

In the MitoTRACER method, the recipient cells constitutively express a red (DsRed-Express2) fluorophore (hereafter referred to as red cells) until they receive mitochondria from the donor cells, which triggers the removal of the red fluorescence expression and activates the permanent expression of the eGFP green fluorophore in these cells (hereafter referred to as green cells; Fig.3a). In brief, the SVZ-NSCMitoTRACERdonor cells were genetically labelled with a mitochondrial-anchored Cre recombinase, and recipient 4T1 cancer cells were equipped with a loxP-DsRed-Express2-Stop-loxP-eGFP switch (Fig.3b). The addition of a specific tobacco etch virus protease (TEVp) proteolytic cleavage site to the donor construct was required for the release of the mitochondria-bound Cre recombinase and its nuclear translocation into the recipient cells (Fig.3cand Extended Data Fig.6a). In the MitoTRACER coculture, following the transfer of nerve-derived mitochondria to the recipient cancer cell, the Cre recombinase triggers the red-to-green switch of the recipient cells, as observed by using time-lapse fluorescence microscopy and flow cytometry analyses of the nerve–cancer coculture (Fig.3d,e). Time-lapse microscopy of the MitoTRACER coculture revealed the development of a green signal in recipient cells following the establishment of tunnelling nanotube connections, further supporting the role of these structures in mitochondrial transfer (Extended Data Fig.6band Supplementary Video3). The use of various donor cell types, such as SVZ (Supplementary Video4), mouse embryonic fibroblasts (Supplementary Video5) and 3T3-L1 pre-adipocytes (Supplementary Video6), demonstrated successful mitochondrial tracking through time-lapse imaging.

a,b, MitoTRACER strategy. Donor neurons express mitochondria-targeted Cre recombinase (iCre) with an SV40 nuclear localization signal (NLS-iCre), linked to the OMP25 mitochondrial outer membrane domain. Recipient cells express both a loxP-DsRed-Express2-Stop-loxP-eGFP switch and the TEVp. After transfer, TEVp cleaves NLS-iCre, enabling nuclear localization and excision of DsRed-Express2, resulting in a permanent change from DsRed (red) to eGFP (green) expression. LTR, long terminal repeat.c, 4T1loxP-DsRed-Express2-Stop-loxP-eGFPco-expressing both MitoTRACER and TEVp shows efficient NLS-iCre cleavage and eGFP expression activation. No unintended cleavage was detected in the absence of TEVp expression. The same sample extracts were loaded in different gels. Representative experiment (n= 3). FL, full length.d,e, Confocal microscopy (d) and flow cytometry (e) of SVZ neuron–4T1 coculture confirmed red-to-green conversion, confirming the transfers and suitability of the approach for high-throughput analysis and collection of recipient cells. WT, wild type. Scale bar, 50 μm.f, Dose-dependent increase of mitochondrial transfer with donor-to-recipient ratios (1:1 to 4:1). The centre line shows the median, the box boundaries show the 25th and 75th percentiles, and the whiskers show the minimum and maximum values. Student’s two-tailed unpairedt-test, ****P< 0.0001 (n= 5 independent cocultures).g, Time-dependent and cumulative increase of mitochondrial transfer from day 1 to day 3. The centre line shows the median, the box boundaries show the 25th and 75th percentiles, and the whiskers show the minimum and maximum values. Student’s two-tailed unpairedt-test, ****P< 0.0001 (n= 6 independent cocultures).h,i, Western blot (h) and densitometry analysis (i) of the MitoTRACER subcellular localization through the HA tag expression on the construct confirmed its mitochondrial localization. Expression of COXIV and HSP90 validated the subcellular fraction purity. The same sample extracts were loaded in different gels. Representative experiment (n= 3) (mean ± s.d.;n= 4 independent experiments). C, cytoplasmic fraction; M, mitochondrial fraction.j, Coculture using recipient cells lacking TEVp confirmed the signal’s specificity. The line shows the median, the box boundaries show the 25th–75th percentiles, and the whiskers show the minimum and maximum values. Student’s two-tailed unpairedt-test,P= 2.16463 × 10−10, ****P< 0.0001 (n= 5 independent cocultures).a, Created in BioRender. S. Grelet (2025)https://biorender.com/aa3gfx0.b, Created in BioRender. S. Grelet (2025)https://biorender.com/ytn18rx.

No unintended signal activation was observed in recipient cells until donor cells were introduced into the coculture. Dose–response (Fig.3f) and time-course (Fig.3g) analyses of the coculture demonstrated the sensitivity of the MitoTRACER approach in capturing the dynamics of mitochondrial transfer and its cumulative capability. Subcellular fractionation of MitoTRACER expression confirmed that nearly all of the construct is anchored to mitochondrial organelles (Fig.3h,i). Finally, MitoTRACER coculture with recipient cells lacking TEV protease showed no signal, indicating that no unintended self-cleavage and further transfer through non-mitochondrial routes, such as secretory pathways, occurred (Fig.3j). Together, these findings demonstrate the efficacy of this system in allowing real-time observation of mitochondrial transfer to recipient cells and its permanent nature, enabling lineage tracing of recipient cells.

We used MitoTRACER to examine the biological impact of the nerve–cancer transfer of mitochondria in the recipient cancer cells. We also examined their fate during cancer progression, both in vitro and in vivo. From the MitoTRACER coculture, we FACS-isolated cells that have received mitochondria from the neurons (4T1MitoTRACERGreen) and those that have not (4T1MitoTRACERRed). Separate subcultures of the red versus green cells obtained from the MitoTRACER coculture revealed distinct growth patterns. Green cells having received mitochondria had a higher propensity for anchorage-independent growth patterns and exhibited the development of spheres throughout the culture (Fig.4a). Anchorage-independent growth capacities of cancer cells are associated with their stemness potential and are usually tied to specific metabolic profiles driven by mitochondrial metabolism33,34,35. We confirmed the increased stemness of the recipient cancer cells by a mammosphere formation assay (Fig.4b). Metabolic profiling of mitochondrial metabolism of the recipient cells indicated enhanced respiratory capacities of the green recipient cancer cells compared to red cells exposed to the neurons but not having received mitochondria and to parental 4T1 cells that were not exposed to neurons (Fig.4c,d). Energetic mapping of the red versus green cells indicated a shift towards a more energetic status, corresponding to increases in both the oxygen consumption rate and the extracellular acidification rate (Fig.4e), and this metabolic shift was associated with a significant increase in basal and maximal respiration capacities of the cancer cells, increased coupling and increased ATP production (Fig.4d,f). Analysis of the functional outcomes of mitochondrial transfer revealed improved redox balance in the recipient (green) cancer cells, as indicated by their higher levels of reduced glutathione (GSH) (Fig.4g). This increase translated into an improved GSH to oxidized glutathione (GSSG) ratio (Fig.4h), a key marker of cellular redox status36. We next confirmed that enhanced redox balance was associated with an increased capacity of recipient cancer cells to withstand oxidative stress (Fig.4i) at physiological doses (9–75 µM) and greater resistance to shear stress (Fig.4j).

a, 4T1 recipient (green) cells sorted after MitoTRACER coculture showed spontaneous sphere formation capacities. Scale bar, 500 μm.b, Mammosphere formation assay confirmed increased stemness potential in green cells (***P= 0.0005,n= 8 independent cultures). Student’s unpaired two-tailedt-test.c,d, Recipient green cells show enhanced mitochondrial OXPHOS capacities. CC, coculture. Representative profile (n= 3); mean ± s.e.m., parental:n= 11, red:n= 13, green:n= 12 cell cultures; Student’s unpaired two-tailedt-test, *P= 0.037, ***P< 0.001, ****P< 0.0001.e, Energy map shows a more aerobic and energetic phenotype in green cells compared to the red counterpart (mean ± s.d.n= 24). ECAR, extracellular acidification rate.f, Luminescence-based assay of total cellular ATP content showed significantly higher levels in green cells (**P= 0.0044,n= 3 independent cultures). RLU, relative light units normalized per cell. Mean ± s.d., Student’s two-tailed pairedt-test.g,h, Green cells exhibited increased GSH (***P= 0.0002,n= 6 independent cultures) and overall improved GSH/GSSG ratios (**P= 0.0028 (g), **P< 0.01 (h),n= 6 independent cultures). Mean ± s.d., Student’s two-tailed pairedt-test. FC, fold change.i,j, Green cells exhibited higher tolerance to H2O2-induced oxidative stress (i; mean ± s.d., Student’s two-tailed pairedt-test, NS, not significant; two-way ANOVA,P= 0.0007) and greater resistance to shear stress (j; two-way ANOVA,P< 0.001). Representative profile (n= 3).k, Modified Boyden chamber assay revealed no change in intrinsic invasion capacity between green and red cells in vitro (NS, not significant; Student’s two-tailedt-test,n= 3 independent cultures).l, In vivo, metastatic progression was enhanced in green versus red cells in mouse mammary fat pad xenografts, as observed by increased liver metastasis (mean ± s.d.; Student’s two-tailed unpairedt-test, *P= 0.01997,n= 8 mice).m, Haematoxylin and eosin liver sections showed metastatic lesions, and Ki67 immunostaining confirmed the cancerous character of the lesions. Scale bar, 4 mm.n,o, Multispectral imaging in human breast cancer and matched metastatic sites revealed increased mitochondrial content at metastatic versus primary sites. Representative image from patient breast cancer samples (n) and matching metastasis and mitochondrial score quantification curve in the samples set (o). The line shows the median, the box boundaries show the 25th and 75th percentiles, and the whiskers show the minimum and maximum values (two-sided Welch’st-test, ***P< 0.001;n= 8 patients).k,l, Created in BioRender. S. Grelet (2025)https://biorender.com/35rzyib.

Together, metabolic plasticity, improved redox balance and increased oxidative and shear stress resistance are hallmarks of metastatic cancer cells37,38,39. We therefore investigated whether mitochondrial transfer could enhance metastatic behaviours in cancer cells. In vitro analysis of the recipient cancer cells’ invasive potential revealed no increase in intrinsic invasiveness (Fig.4k). However, in vivo xenografts with recipient cells demonstrated significantly higher metastatic potential than their ‘non-recipient’ counterparts (Fig.4l,m), suggesting that mitochondrial transfer may contribute to the metastatic cascade beyond the invasion process.

Pathological analysis of human breast cancer samples further underscored the role of mitochondria in tumour dissemination, with metastatic cells exhibiting a significant increase in mitochondrial load (Fig.4n,oand Supplementary Tables8and9). Together, our findings clearly suggest that mitochondrial transfer from neurons to cancer cells may enhance metastatic behaviour by strengthening the resilience of cells against metastatic stressors, such as oxidative and shear stress, thereby augmenting their metastatic potential through adaptive mechanisms.

The observed differences in mitochondrial load between primary and metastatic cancer suggested that mitochondria-recipient cancer cells in primary tumours may possess enhanced metastatic potential. The metastatic cascade is inefficient, with cancer cells encountering multiple stressors that impede their successful dissemination and growth at secondary sites40. Metabolic reprogramming and plasticity have emerged as crucial adaptive mechanisms for the successful metastatic dissemination of cancer cells5,41,42,43. Thus, we reasoned that neuron-derived mitochondria confer enhanced metabolic adaptability and resilience, enabling recipient cancer cells to better spread and ultimately survive and proliferate at distant sites. To test this hypothesis, we performed lineage tracing to follow the fate of primary cancer cells receiving mitochondria from nerves.

We developed a preclinical model of the nerve–cancer transfer of mitochondria (Fig.5). We combined the MitoTRACER coculture approach with the 4T1 mammary fat pad xenograft model that can fully recapitulate the breast cancer progression steps and metastasis of TNBC44,45. Mixed-cell spheroids of the MitoTRACER coculture (Fig.5a) were transplanted into mammary fat pads. When the xenograft reached appropriate size, cancer cells were collected from the primary tumour and both lung and brain tissues, which are prevalent metastatic sites in TNBC. The cells were then analysed by flow cytometry to probe the ratio between 4T1Red+versus 4T1Green+cells (Fig.5b). We observed the development of an average of about 5.4% 4T1GFP+cells within the primary tumour. This proportion was significantly enriched in both lung and brain tissues to reach 27.3% and 46.0% of the total cancer cells in those sites, respectively (Fig.5c). This demonstrates that cells that acquired mitochondria from the SVZ-NSCs within the primary tumour, or their progeny, are more likely to form distant metastases successfully. To extend the observations with implanted mixed-cell spheroids to a model of host-mediated neuronal mitochondrial transfer, we genetically modified mouse DRG in vivo using a lentivirus encoding a MitoTRACER construct driven by the synaptin1 promoter (LV-CRE-OMP25) or its non-transferable, nuclear-targeted variant (LV-CRE-NLS; Fig.5d). At 10 days after lentiviral transduction, SWITCH-TEVp-expressing recipient cells were injected into the mammary fat pad, and tumours were subsequently excised to assess eGFP fluorophore expression, indicative of mouse-derived neuronal mitochondrial transfer into cancer cells (Fig.5e). We observed only 1.6% eGFP+cells in the primary tumour, which was less than the 5.4% observed with the mixed-cell spheroids, probably owing to the lower nerve density in the tumour in situ, when compared to the spheroid model. Lineage tracing reproduced by labelling the endogenous nerve of the house mouse using lentiviruses injected into the DRG area also demonstrated a marked enrichment of eGFP+recipient cells in metastatic sites, with significant increases in the brain and liver compared to the primary tumour (Fig.5f).

a, 3D spheroids of SVZ-NSCsMitoTRACERmixed with 4T1 recipient cancer cells confirmed mitochondria transfer, evidenced by red-to-green fluorescence conversion in 4T1 cells (white arrows).b,c, MitoTRACER spheroids were transplanted into mammary fat pads of BALB/c mice (b); after cancer progression, flow cytometry of cells isolated from primary tumours, lungs and brains revealed selective enrichment of green fluorescent cells (eGFP+) in lung and brain metastases compared to the primary tumour (c). Mean ± s.d.; Student’s two-tailed pairedt-test, lung: *P= 0.018, brain: ***P= 1.2275 × 10−6;n= 9 mice; ANOVA,P=0.0005.d, Schematic of the procedure for lineage tracing of mitochondrial transfer between mouse host mammary neurons and cancer cells in vivo. Lentivirus expressing Syn1-MitoTRACER construct was injected into the DRG area, innervating the lower mammary fat pads. After 10 days, the SWITCH-TEVp-expressing 4T1 recipient cells were injected. Following tumour growth, tissues (primary, lung, brain and liver) were analysed for mitochondrial transfer by flow cytometry.e, Flow cytometry of primary tumours showed eGFP+cells, confirming mitochondrial transfer. Control lentivirus expressing Syn1-driven nuclear-localized Cre (LV-CRE-NLS) showed no green signal, validating the specificity of our approach.f, Ex vivo quantification of mitochondrial transfer in the primary tumour and lineage tracing of metastatic development in the lung, brain and liver showed significant enrichment of eGFP+cells in metastatic sites compared to the primary tumour (ANOVA,P< 0.0001), with significant enrichment in brain and liver metastases (mean ± s.d., Student’s one-tailed pairedt-test values, brain: **P=0.0096, liver: *P= 0.0215;n= 5 mice). Top panel shows distribution of red and green cells.g, Validation using syngeneic B16-F1 melanoma cells co-injected with SVZ-NSCsMitoTRACERas mixed-cell spheroids in C57BL/6 mice showed significant mitochondrial transfer enrichment in brain metastases (mean ± s.d.,  ANOVA,P= 0.0009; Student’s one-tailed pairedt-test, brain: **P= 0.0017;n= 8 mice). Top panel shows distribution of red and green cells.b, Created in BioRender. S. Grelet (2025)https://biorender.com/culfgzj.d, Created in BioRender. S. Grelet (2025)https://biorender.com/2w7o7rr.

We tested whether nerve-mediated mitochondrial transfer may affect distant metastasis in another cancer progression model using the B16-F10 melanoma xenograft. The observations about transferred mitochondria in primary and metastatic tumours were similar in the melanoma and breast cancer models, but the rate of mitochondrial transfer was lower in the melanoma than in the 4T1 breast cancer model. Furthermore, the melanoma model showed no significant enrichment of green eGFP+cells in the lungs or liver but marked enrichment in the brain, corroborating findings from the breast cancer model and providing evidence that neuronal mitochondrial transfer may promote brain metastases with various cancer types (Fig.5g).

The nerve–cancer interplay was initially discovered in prostate cancer, in which cancer cells can promote cancer innervation by expressing axon guidance molecules, such as semaphorins, promoting neuronal progenitor differentiation in the cancer stroma and the establishment of the nerve–cancer interface17,46,47,48. Initial observations in prostate cancer have shown how the aberrant expression of semaphorin 4F by cancer cells promotes increased prostate nerve density and cancer aggressiveness18,49. Later studies supported the relevance of semaphorin 4F in other biological contexts such as gastric cancer19and breast cancer17, in which cancer cell plasticity triggers the expression of semaphorin 4F to promote cancer innervation and metastasis50. Our previous studies demonstrated how cancer innervation and cancer cell plasticity are intimately linked17,50, and although progress has been made in understanding the mechanisms leading to cancer-mediated neuronal differentiation and the establishment of the nerve–cancer interface, the mechanisms and functional effects of neurons on breast cancer progression and metastatic potential remain incompletely understood. Our study demonstrates how cancer-induced neuron differentiation leads to marked neuronal metabolic reprogramming with clear functional consequences. We show that cancer-associated neurons are a significant source of mitochondria transferred to the cancer cell to induce their metabolic reprogramming and increase stemness potential.

Pathological analysis and animal denervation models, including the breast carcinoma denervation model presented herein, have consistently associated cancer nerve density with cancer invasion and metastasis51,52,53,54. Metabolic plasticity and stemness potential are essential hallmarks of cancer metastasis55. To test whether the nerve–cancer transfers of mitochondria relate to the development of distant metastasis, we developed the MitoTRACER genetic reporter capable of permanently marking recipient cells. We used this approach to create a preclinical model of nerve–cancer transfer of mitochondria in vivo. The fate mapping experiment revealed enrichment of mitochondria-recipient cancer cells or their progeny at metastatic tumour sites relative to the primary tumour, reflecting their increased capacities to achieve the metastatic colonization steps successfully. Although the exact molecular mechanisms remain unclear, our study shows that cancer cells acquiring mitochondria from neurons gain adaptive advantages, enabling resilience against metastatic stressors such as oxidative39and shear stress37,38,56, which are well-known key barriers to metastasis. From a probabilistic point of view, cancer cells that energetically outperform others have a greater probability of moving and seeding. It is probably a selection process by which the incorporation of these neuronal mitochondria into the cancer cells provides an increased capacity to survive the metastatic process.

We noted selective enrichment of mitochondria-recipient cancer cells in brain metastases in both breast cancer and melanoma models. This suggests that mitochondrial acquisition from neurons may prime cancer cells to adapt more effectively to the brain’s unique microenvironment. Previous research supports this idea, as brain-metastatic cancer cells must metabolically adapt to survive in the brain’s nutrient-poor environment57due to the high energy demands of neurons58.

Neuronal mitochondria are renowned for their highly efficient metabolic potential, which contributes to the superior metabolic efficiency of neurons compared to epithelial cells59. Therefore, cancer cell acquisition of neuron-derived mitochondria may provide them with an adequate metabolic arsenal to survive and successfully proliferate in the brain environment. Our findings also suggest that neurons may be particularly effective donors of mitochondria, probably owing to their abundant mitochondrial content and the establishment of strong contacts at the nerve–cancer interface. Further studies are needed to clarify whether the impact of neuronal mitochondria is primarily due to their intrinsic metabolic efficiency or simply their highest transfer frequency.

Collectively, our findings provide a compelling metabolic explanation for the observed dependency between cancer cells and nerves, potentially extending to broader contexts. These results advocate for more in-depth studies into underlying mechanisms and therapeutic strategies targeting nerve–cancer mitochondrial transfers to prevent metastatic disease.

The 4T1 (CRL-2539), Neuro2A (CCL-131) and NMuMG (CRL-1636) cell lines were purchased from ATCC. HT-22 (ESA111) and 3T3-L1 (EF3001) were purchased from Kerafast. The 50B11 rat DRG neuronal cells were obtained from the GRCF Biorepository & Cell Center (Johns Hopkins University). MEFs were provided by S. Lloyd. Primary SVZ-NSCs were isolated from fresh BALB/c mouse brains as described previously60(Extended Data Fig.2a–d). The cell lines were routinely checked for mycoplasma contamination by using the Mycostrip Test Kit (Invivogen), and through Hoechst-33342 staining. The cell lines used on this study has been authenticated by genomic profiling. Cancer cells were cultured using Dulbecco’s modified Eagle’s medium (DMEM; Life Technologies) supplemented with 5% fetal bovine serum (FBS; Fisher Scientific), 5% bovine calf serum (Fisher Scientific), antimicrobial solution (penicillin, streptomycin and amphotericin B; 1% antibiotic–antimycotic, Life Technologies) and antimycoplasma solution (1.25 mg l−1plasmocin prophylactic, InvivoGen). SVZ-NSCs were cultured in specialized medium (Neurobasal, Life Technologies) supplemented with serum-free supplement (2% B-27, Life Technologies), 1 mg l−1heparin (STEMCELL Technologies),l-alanyl-l-glutamine dipeptide supplement (1% GlutaMAX, Life Technologies), antimicrobial solution (1% antibiotic–antimycotic, Life Technologies), 10 ng ml−1recombinant human epidermal growth factor (Bio-Techne) and 10 ng ml−1of recombinant human fibroblast growth factor basic (Bio-Techne). 50B11 cells were cultured in specialized medium (Neurobasal, Life Technologies) supplemented with both 10% serum and the serum-free supplement (2% B-27, Life Technologies),l-alanyl-l-glutamine dipeptide supplement (0.27% GlutaMAX, Life Technologies), antimicrobial solution (1% antibiotic–antimycotic, Life Technologies) and 0.22% glucose. SVZ-NSC adherence in vitro was aided by adding 0.25 µl ml−1ECMatrix-511 Silk E8 Laminin Substrate (MilliporeSigma). For serum-free experiments, FBS and bovine calf serum were substituted with a serum-free medium (10% Knockout Serum Replacement (KSR), Life Technologies). Dialysed FBS (Thermo Fisher Scientific) was substituted for FBS in uridine depletion assays. ρ0cells were cultured in medium supplemented with 50 µg ml−1uridine (Sigma-Aldrich) and 1 mM sodium pyruvate (Fisher Scientific). Antibiotic selection pressure was applied with combinations of 10–20 µg ml−1puromycin (Fisher Scientific), 100–200 µg ml−1phleomycin D1 (Zeocin, InvivoGen) and 10–20 µg ml−1Blasticidin (InvivoGen) to maintain transgene expression. Cells were cultured in an incubator at 37 °C and 5% carbon dioxide and were tested routinely for mycoplasma contamination.

In short-term nerve–cancer coculture experiments (≤5 days), 4T1 cells were initially seeded in DMEM and allowed to adhere. After 1 day, the medium was changed to KSR medium, and SVZ-NSCs were added the following day. For long-term coculture experiments, 4T1 cells were maintained in DMEM, with SVZ-NSCs added weekly and 4T1 cells passaged twice a week. SVZ-NSCs were introduced into the 4T1 cell cultures at 2–4% of the total cell population, matching the percentage observed in primary cancers formed from 4T1 cell xenografts61. For the 3D coculture and reconstruction imaging experiment, a total of 1.2 × 10550B11CCO-GFPcells and 1.2 × 1054T1mCherrycells (both in 10 µl PBS) were mixed with 60 µl of Growth Factor Reduced Matrigel (Corning, catalogue number 356231) and seeded onto the glass area of uncoated glass-bottomed 35-mm dishes (MatTek). The dishes were incubated at 37 °C in a humidified atmosphere of 5% CO2in air for 1 h. Following incubation, 2 ml of Neurobasal complete medium containing 75 µM forskolin (MilliporeSigma) was added. Confocal images were captured 2 days later using a Nikon A1 confocal laser microscope.

Lenti-X 293T cells (Takara Bio USA) were used for lentivirus generation, and Phoenix-AMPHO cells (ATCC) were used for retrovirus production. Cells were plated in complete DMEM 1 day before transfection using a transfection reagent (FuGENE, 3 µl µg−1of plasmid; Promega).

For lentivirus production, Lenti-X 293T cells were transfected by lipofection with the following components: 1,000 ng of transfer plasmid, 750 ng of packaging plasmid (psPAX2; Addgene) and 250 ng of vesicular stomatitis virus G envelope-expressing plasmid (pMD2.G; Addgene). Both psPAX2 and pMD2.G were gifts from D. Trono, École Polytechnique Fédérale de Lausanne (Addgene plasmid numbers 12260 and 12259). For in vivo work, the lentivirus was further concentrated in a concentrator solution (40% PEG-8000; 1.2 M NaCl in PBS) added at a 1:3 ratio of solution/supernatant, gently shaken for 2 days at 4 °C, collected through 1,600gcentrifugation and resuspended in a 1:100 volume of PBS and stored at −80 °C.

For retrovirus production, Phoenix-AMPHO cells were transfected with 1,000 ng of plasmid DNA. The medium was replaced with KSR-containing medium 1 day after transfection, and the viral supernatant was collected over a period of 2 days. For viral transduction, target cells were plated at 60% confluence and transduced in the presence of 10 µg ml−1Polybrene (MilliporeSigma) in KSR medium. SVZ-NSCs were cultured in Neurobasal medium without Polybrene to maintain cell viability. The pLenti-CMV-GFP-puro plasmid was provided by E. Campeau (Addgene plasmid number 17448). The pLVX-EF1a-CCO-IRES-puromycin plasmid was provided by D. Andrews (Addgene plasmid number 134861). The pFUGW mCherry-KASH plasmid was a gift from H. MacGillavry (Addgene plasmid number 131505). The CMV-loxP-DsRed-loxP-eGFP plasmid was provided by D. Gilkes (Addgene plasmid number 141148). The pLenti-CMV-puro-2A-TEVp plasmid was a gift from M. Tripodi (Addgene plasmid number 99610).

ρ04T1 cell lines were developed as described previously28. The complete and permanent loss of mtDNA in 4T1 cells was achieved by transient co-overexpression of a mUNG1 and UL12.5M185 herpesvirus protein. 4T1mCherrycells were cultured on 6-cm plates. At 1 day after seeding, cells were lipofected using 15 µl of FuGENE transfection reagent and 2,500 ng each of pMA4008 and pMA3790 plasmids provided by M. Alexeyev62. Lipofection was performed on medium with uridine and sodium pyruvate and without antibiotics, and medium was changed after 1 day to fresh medium with uridine, sodium pyruvate and antibiotics. Single-cell sorting based on strong positive eGFP expression (top 3%) was performed after 3 days into 96-well plates supplemented with uridine, sodium pyruvate and 20% FBS. The medium was changed weekly and twice weekly when colonies were observed.

The MitoTRACER donor cells were prepared using retroviral or lentiviral transduction, followed by selection with 20 µg ml−1Blasticidin. To further purify cells with the highest MitoTRACER expression, we incorporated a GFP-11 strand into the MitoTRACER construct, allowing for additional selection by FACS. If necessary, the MitoTRACER cells were further purified to ensure consistent and robust expression by transiently transfecting them with 5,000 ng of GFP1-10 plasmid and sorting based on GFP signal 3 days later. The plasmid pQCXIP-GFP1-10 was provided by Y. Hata (Addgene plasmid number 68715). Loss of GFP signal following transient selection was confirmed by flow cytometry before utilizing the cells in MitoTRACER coculture experiments.

DNA was extracted from cells using a kit (GeneJET Genomic DNA Purification Kit, Thermo Fisher Scientific). Reaction mixes were used for end-point PCR (DreamTaq Hot Start Green PCR Master Mix, Thermo Fisher Scientific) and quantitative PCR (PowerTrack SYBR Green Master Mix, Thermo Fisher Scientific). The quantitative PCR was performed with a thermocycler (QuantStudio 7 Pro, Thermo Fisher Scientific). Sample concentration was adjusted to 10 ng µl−1, and PCR was performed with 20 ng of DNA template. The end-point PCR products were separated on acrylamide/bis-acrylamide (29:1) 10% gels (Thermo Fisher Scientific) and stained with ethidium bromide (Thermo Fisher Scientific). Mouse mtDNA amplification was performed as described previously63using the mMitoF1 5′-CTAGAAACCCCGAAACCAAA-3′ and mMitoR1 5′-CCAGCTATCACCAAGCTCGT-3′ primers. The mtDNA copy numbers were normalized to the nuclear copy numbers using the B2MF1 5′-ATGGGAAGCCGAACATACTG-3′ and B2M-R1 5′-CAGTCTCAGTGGGGGTGAAT-3′ primers. Uncropped gels are provided in Supplementary Fig.1.

Western blot was performed by sodium dodecyl sulfate polyacrylamide gel electrophoresis. Whole-cell lysates were prepared from 2–5 × 106cells in 300 µl of lysis buffer (20 mM Tris, pH 7.4, 1% Triton X-100, 10% glycerol, 137 mM sodium chloride, 2 mM ethylene diamine tetraacetic acid, 1 mM sodium orthovanadate) and protease inhibitors (Thermo Fisher Scientific). Subcellular fractionation was performed using the Mitochondria Isolation Kit for Cultured Cells (Life Technologies) according to the manufacturer’s instructions. With mitochondrial-bound proteins, sonication (30 pulses, each 20 s) was performed instead of centrifugation of excess material to ensure that mitochondrial proteins remained in suspension. Lysates were clarified by centrifugation at 4 °C for 30 min at 17,000g. Typically, whole-cell lysates (5–20 µg) were separated on 10% or 12% acrylamide minigels and transferred to a membrane (Immuno-Blot, Bio-Rad).

Samples were separated with electrophoresis (Any kD gels and Mini-PROTEAN apparatus; Bio-Rad). The membrane was blocked for 30 min in wash buffer (0.1% Tween 20 in PBS) containing 5% nonfat dry milk and incubated overnight with primary antibody (Cre recombinase (1:1,000, number 15036, lot 2), mCherry (1:1,000, number 43590, lot 2) or α-tubulin (1:1,000 number 2144, lot 5), Cell Signaling Technology; or GFP (1:200, number sc-9996, lot E2521) or HA tag (1:1,000, number sc-7392 lot I0992), Santa Cruz Biotechnology) that was diluted in the same buffer. After extensive washing, the blot was incubated with secondary antibody (1:1,000, Thermo Fisher PI31430 or PI31460) for 30 min in blocking buffer, washed and processed with a western blot detection system (ChemiDoc MP Imaging System, Bio-Rad). Uncropped blots are provided in Supplementary Fig.1.

Cell imaging was performed with a fluorescence (BZ-X810, Keyence) or inverted confocal microscope (A1R, Nikon). Cell sorting was performed using FACS (FACSAria II (BD Biosciences) with analysis of GFP (laser excitation, 488 nm; detection, 525/30 band-pass filter width (BP)) and DsRed-Express2 and mCherry Red fluorophores (561 nm; 610/20 BP). The gating strategy is provided in Supplementary Fig.2.

Dye staining (MitoTracker, Invitrogen) was performed according to the protocol from the manufacturer. Sorting was performed with 1 cell per well in 96-well plates (single-cell sorting) or 5,000 cells per well in 6-well plates containing 2 ml of DMEM with 10% FBS. Immunostaining and cytometry were performed as follows: cells were collected by trypsinization, fixed in 4% paraformaldehyde for 15 min at room temperature, and permeabilized with 0.1% Triton X-100 diluted in 0.5% bovine serum albumin. Fixed cells were incubated for 1 h at room temperature with primary antibodies, washed in PBS and incubated for 30 min at room temperature with an Alexa Fluor 647 secondary antibody (1:10,000, Thermo Fisher A21235, lot 2836809; A21244, lot 2674387). For animal xenograft samples, deparaffinized slides were microwave-boiled (heat-induced epitope retrieval) with 10 mM sodium citrate, blocked in 1% bovine serum albumin, incubated in primary antibodies overnight, washed, incubated in secondary antibodies, Alexa 647 and Alexa 568, for 40–60 min (1:250, Thermo Fisher A78952; lot 3034155, A11004, lot 2198584) and counterstained with 4′,6-diamidino-2-phenylindole (DAPI). The following primary antibodies were used: O4 (1:10,000, R&D Systems MAB1326, lot HWW143051), ALDH1L1 (1:10,000, OriGene TA501868S, lot A01), Map2 (1:1,000 Proteintech 17490-I-AP, lot 00126726), GFP (1:200, Invitrogen A10262, lot 2738237) and β3 Tubb3 (1:400, Santa Cruz sc80005, lot A1821).

Whole-cell current clamp recordings were performed on neurons differentiated from cocultures of SVZ-neural progenitor cell mixed with 4T1 cancer cells using an Axopatch 200B amplifier and Digidata 1322A, with data acquisition through pClamp 8 software (Molecular Devices). A Zeiss Axiovert microscope with epifluorescence was used to identify and record neural cells expressing eGFP in the nerve–cancer coculture. Cells were clamped at −80 to −90 mV, and depolarizing currents were injected using 10 current ramp increments from 20 to 200 pm (800 ms duration), interleaved by 5-s intervals. Membrane voltages were sampled at 5 kHz and filtered at 2 kHz. The action potential threshold was determined from the dV/dtderivative calculated using IgorPro 6 (WaveMetrics)64,65.

The external bath solution consisted of 120 mM NaCl, 6 mM KCl, 1 mM CaCl2, 1 mM MgCl2, 10 mM HEPES and 20 mM glucose (pH 7.4). Patch pipettes were filled with an internal solution containing 125 mM K+gluconate, 4 mM KCl, 4 mM NaCl, 1 mM MgCl2, 10 mM HEPES, 4 mM MgATP, 0.3 mM Na2GTP and 10 mM phosphocreatine (pH 7.26). Osmolarity of all solutions was verified with an Osmette III osmometer (Precision Systems). Series resistance and membrane capacitance were not compensated but were measured at the beginning and end of each recording to ensure quality. Electrophysiology figures were prepared using IgorPro and GraphPad Prism 9.

The ρ04T1mCherry+cells were cocultured with SVZ-NSCsCCO-GFPand sorted on the basis of mCherry and eGFP expression using flow cytometry. The sorted cells were seeded on DMEM with dialysed FBS, with or without uridine and sodium pyruvate, and stained with 2% crystal violet in 20% methanol. The rescued cells were cultured in the presence of antibiotic (Zeocin), associated with the cancer cell expression of the mCherry fluorophore. 4T1 cells were checked for the presence of mCherry and the absence of eGFP expression to confirm the 4T1 nature of the collected cells and exclude contamination with SVZ-NSCs during FACS.

Cells were plated at 4 × 104cells per well (Seahorse XFe96 plate, Agilent) on the day before assay, and the assay tools and calibrant solution (Seahorse FluxPak, Agilent) were acclimated in an incubator without carbon dioxide on the night before assay according to the protocol from the manufacturer. An assay kit was used (Seahorse XF Mito Stress Test Kit, Agilent) with oligomycin (1.5 µM), FCCP (1.0 µM), and rotenone and antimycin A (0.5 µM). On the assay day, the medium was changed (Seahorse XF DMEM, Agilent) and supplemented with 10 mM glucose (Agilent), 1 mM pyruvate (Agilent) and 2 mM glutamine (Agilent). Uridine and sodium pyruvate were added for assays with ρ0cells. Assays were performed according to the protocol from the manufacturer, and wells were normalized by direct cell counting to determine the absolute number of cells per well. Cells were stained with Hoechst dye and counted directly with an automatic microscope reader (Celigo S, Nexcelom Bioscience). The absence of cross-contamination following the sorting of coculture was confirmed by the analysis of the presence of eGFP+cells in the mCherry-FACS-sorted cancer cell samples.

In cases in which inserts were used to physically block contact between cell types, 100,000 4T1mCherrycells were seeded in the bottom well of a 6-well Transwell system (Corning, 3412) or 6-well plate, and the following day, the medium was changed to 10% KSR-containing DMEM, and 200,000 SVZ-NSCsCCO-GFPwere added either to the insert or directly to the well. For the dose–response and time–response assays, 4T1mCherrycells were seeded at 100,000 cells per well and changed to KSR-containing medium the next day. For the dose–response assay, between 25,000 and 400,000 SVZ-NSCsMitoTRACERwere then added to the culture. For the time–response assay, 200,000 SVZ-NSCsMitoTRACERwere added daily until the time of FACS analysis. FACS analysis was carried out on the third day after adding SVZ-NSCsMitoTRACER.

The MitoTRACER construct was cloned into the pMXs-IRES-Blasticidin retroviral vector RTV-016. We used the HA-MITO plasmid pMXs-3XHA-EGFP-OMP25 as a template. pMXs-3XHA-EGFP-OMP25 was a gift from D. Sabatini (Addgene plasmid number 83356). The OMP25 mitochondrial tagging element and attached 3′ linker (Pro-Arg-His) were conserved. The HA-MITO plasmid was digested to insert our construct design consisting of the assembly of GFP-11–HA tag–TEVp cleavage recognition site–linker–iCre recombinase–SV40-NLS. The GFP-11 fragment was added between the OMP25 and TEVp cleavage site for convenience to generate cells having the desired amount of reporter using GFP1-10 fluorescence complementation and FACS. The construct also included a Blasticidin S deaminase gene separated by an internal ribosome entry site (IRES), enabling the selection of cells that were transduced successfully by the retroviruses. All constructs were validated by Nanopore sequencing. For in situ labelling, the MitoTRACER construct was transferred into a lentiviral construct with theSyn1promoter (Addgene plasmid number 71427). The CRE-NLS control was generated through site-directed mutagenesis (Q5 Site-Directed Mutagenesis E0554S, New England Biolabs) by removing the OMP25 section of MitoTRACER. Similar constructs were generated by replacing the MitoTRACER with the eGFP fluorophore instead to generate the lentiviral Syn1-GFP-OMP25 and Syn1-GFP-NLS constructs.

The recipient cell construct for the MitoTRACER method was a DsRED and eGFPloxPswitch with TEVp expression under the control of an SV40 promoter and tied to a BleoR gene by a 2A self-cleaving peptide. This construct was derived from the CMV-loxP-DsRed-loxP-eGFP plasmid (Addgene number 141148) that was digested with DraIII and SwaI restriction enzymes to insert the F2A-TEVp at the Ct end of BleoR with a DNA assembly kit (NEBuilder HiFi DNA Assembly, New England Biolabs).

Laser-capture microdissection, RNA extraction and cDNA microarray analysis were conducted as previously detailed2,66. Cancer cells from BoNT/A and saline groups were procured through laser-capture microdissection from 8-μm frozen sections cut using a cryostat from the cancerous areas of fresh, unfixed xenograft specimens. The sections were mounted on non-charged glass slides and stored at −80 °C until analysis. Following staining and dehydration as per the manufacturer’s protocol, a laser-capture microdissection was performed before cDNA preparation, and microarray analyses as previously described2. After performing quantile normalization, we carried out principal component analysis as it was noted that samples from the BoNT/A and saline groups have clear distinct expression patterns. A total of 116 genes (135 probes) were identified as differently regulated by an empirical Bayes test with a false discovery rate cutoff <0.30 and a twofold change as an exploratory approach. To identify the biological processes and pathways associated with the differentially expressed genes, we performed a GATHER analysis67and used gene set enrichment analysis68to score association of BoNT/A treatment with Kyoto Encyclopedia of Genes and Genomes (KEGG)69pathways.

Reads were trimmed for quality control to remove low-quality bases and adaptor sequences. Samples failing quality assurance and quality control thresholds (based on read quality, length or count) were excluded. Quality-trimmed reads were aligned to theMus musculusmm39 genome using Bowtie2, in paired-end mode with a 500-bp fragment length. SAM and BAM files were processed with SAMtools and Sambamba, and alignment quality was assessed using SAMtools idxstats, retaining only reads aligning to mm39. Read counts were normalized with DESeq2’s median-of-ratios method, and genes with expression below 1.0 (geometric mean) were filtered out. Gene set ANOVA was performed on gene sets after excluding low expression levels (geometric mean <1.0) on sets containing between 2 and 50 genes. ANOVA was applied using a log-normal model with the criterion Akaike information criterion (AIC), followed by false discovery rate (FDR) correction for multiple comparisons. KEGG pathways were analysed to identify differentially regulated gene sets.

To analyse mtDNA transfer in situ, cancer cells were FACS-sorted to isolate mCherry+4T1 cancer cells from tumour tissues. DNA was extracted from collected cells using the GeneJET Genomic DNA Purification Kit (Thermo Fisher Scientific). The mtDNA region containing the mutation site was amplified using PCR with Q5 High-Fidelity DNA Polymerase (New England Biolabs) with primers designed to amplify a 420-base-pair product around the variation point (forward: 5′-CTAGAAACCCCGAAACCAAA-3′, reverse: 5′-TCATACTAACAGTGTTGCATC-3′).

The PCR reaction was column-purified and analysed using Oxford Nanopore Technology sequencing (Plasmidsaurus). A custom Python script (available via GitHub athttps://github.com/GreletLab/mtDNA-heteroplasmy) was developed specifically for this project to process the obtained raw reads, align them to the target sequence, classify each read as wild type (mouse host-derived) or mutated (cancer cell-derived) and generate statistical output. The script was validated using mtDNA sequence from pure 4T1 cell extract or BALB/c mouse tissue, as well as through the manual counting of some samples to ensure the accuracy of the output.

Intraductal transplantation of human ductal carcinoma in situ cells was performed as previously described13. In brief, before transplantation into 6-week-old virgin female SCID-beige mice, DCIS.COM cells were resuspended as single cells in PBS and counted. A 30-gauge Hamilton syringe, 50-μl capacity, with a blunt-ended 1/2-inch needle, was used to deliver the cells. The mice were anaesthetized, and a Y-incision was made on the abdomen to allow the skin covering the inguinal mammary fat pads to be peeled back to expose the inguinal gland. The nipple of the inguinal gland was snipped so that the needle could be directly inserted through the nipple. Two microlitres of cell culture medium (with 0.1% trypan blue) containing cells at a concentration of 2,500 to 5,000 cells per microlitre was injected; the injected liquid can be visually detected in the duct. The skin flaps were repositioned normally and held together with wound clips.

BoNT/A denervation was performed as previously described8, adapting the procedure through intraductal delivery of the toxin 4 weeks after the intraductal transplantation of the cells, and the experiment was concluded at 10 weeks. As previously described, the loss of nerve function, whether through physical or chemical denervation procedures, resulted in axonal atrophy and decreased intratumoural nerve density with an effect akin to Wallerian degeneration6,16,70,71. BoNT/A was reconstituted in 0.1 ml of 0.9% saline without preservatives to achieve a concentration of 1 U μl−1and was utilized within 2 h of reconstitution. For the intraductal model, 15 U kg−1of BoNT/A per xenograft was carefully injected using a 30-gauge needle and a precision glass syringe, with meticulous attention to avoid extra ductal leakage. For the mammary fat pad injection model, 0.3 units of BoNT/A in 100 µl PBS was injected into the mammary fat pad 1 week before the injection. One million 4T1 cancer cells (100 µl of a suspension of 10 million cells per ml) were injected into the same mammary fat pad. Using calipers, the tumour’s height and width were subsequently estimated.

Cancer cells isolated from the MitoTRACER coculture and established as long-term subcultures were seeded at 5,000 cells per ml in 96-well ultralow-attachment plates. Mammospheres were collected on the seventh day of culture, washed with PBS, and sedimented for 30 min. All mammospheres with diameters greater than 50 µm were counted with a fluorescence microscope (BZ-X810, Keyence).

Cell invasion was assessed using the Corning BioCoat Matrigel Invasion Chamber (Fisher Scientific), following the manufacturer’s instructions. In brief, the Matrigel-coated chambers were rehydrated with serum-free medium for 2 h at 37 °C in a humidified incubator with 5% CO2. After rehydration, 50,000 cells were seeded in the upper chamber in a serum-free DMEM medium; the lower chamber contained DMEM supplemented with 20% FBS as a chemoattractant. The plates were incubated for 24 h at 37 °C. Non-invaded cells and Matrigel were removed from the upper membrane surface using a cotton swab. The invaded cells on the bottom side of the membrane were fixed and stained with DAPI. Fluorescence microscopy was used to image the stained cells, and a custom image analysis script (available via GitHub athttps://github.com/GreletLab/DAPI-count) quantified the number of cells that migrated through the membrane. The script was previously validated by traditional manual counting to compare to the output data and ensure its accuracy. All assays were performed in independent triplicate.

The shear stress assay was performed as previously described38. In brief, cancer cells were resuspended at a concentration of 1 × 105cells per ml and loaded into the injection pump using a 30-ml syringe. The suspension was dispensed through a 26-gauge, 150-mm blunt tip Luer lock needle at a controlled flow rate of 0.25 ml s−1. The cells were subjected to 10 shear cycles, after which cell viability and quantity were assessed using a Celigo image cytometer. Cell viability was determined with Hoechst–propidium iodide staining to evaluate cell quantity and cell death.

To induce oxidative stress, cancer cells were cultured in vitro and treated with increasing concentrations of hydrogen peroxide. Cells were seeded at a density of 50 × 104cells per well in a 96-well plate and allowed to adhere overnight. The following day, cells were exposed to a range of hydrogen peroxide concentrations from 100 mM to 9.375 µM in serum-free medium for 24 h. After treatment, cells were washed with PBS and analysed for recovery before Hoechst–propidium iodide staining to evaluate cell quantity and cell death using the Celigo image cytometer.

Nerve–cancer cocultures were grown in Nunc Lab-Tek II chambered coverglass slides (Roskilde) at 37 °C with 5% CO2. Cells were stained with Fluo-4AM calcium indicator dye and Hoechst 33342 (Invitrogen) just before imaging, following the manufacturer’s instructions. Time-lapse videos were acquired using a Nikon A1r confocal microscope equipped with a stage-top incubator for temperature and gas control, using a 20× objective lens (numerical aperture 0.8). Images were captured every 15 s over a duration of 10 min. Calcium fluorescence intensity was quantified by defining a region of interest around the cell to be measured and calculating the mean fluorescence intensity of the green channel within the region of interest for each time point.

Mixed-cell MitoTRACER spheroids were injected into the mammary fat pads of 6-week-old female BALB/c mice. In accordance with the Institutional Animal Care and Use Committee (IACUC) protocol for this study, tumour volumes were monitored every 2–3 days using digital calipers, and humane end points were applied if tumours exceeded 2 cm3or signs of distress or ulceration appeared. Tumour volume was estimated through (short side × short side × long side)/2. The maximal size set by IACUC was 2 cm3and was not exceeded in any experiment.

At end point, primary tumours, lungs and brains were collected, rinsed in PBS, and processed for dissociation, cell culture and flow cytometry to quantify red-to-green cell ratios in both primary and metastatic sites. Dissociated cell suspensions were filtered through 40-μm strainers, centrifuged at 300gfor 5 min, and seeded in DMEM for further analysis. To confirm metastasis at distant sites, a subset of cells was cultured in vitro with Blasticidin selection for 2–4 weeks before re-examination of red and green fluorescence signals.

Tumour dissociation was performed using the Tumor Dissociation Kit (Miltenyi Biotec 130-096-730); tumours or organs weighing between 0.04 and 1 g were cut into 2–4-mm pieces and placed in 1–2.5 ml of enzyme medium. Tumours weighing more than 0.2 g were placed in C-tubes and processed using the gentleMACS Dissociator both before and after incubation. Samples not placed in C-tubes were manually separated after incubation using a filter. All samples were then strained through a 40-µm or 70-µm filter into a 50-ml tube and rinsed with 10–20 ml of RPMI or DMEM. The samples were centrifuged at 500gfor 7 min, after which the supernatant was discarded, and flow buffer was added.

The human breast cancer tissue microarray was purchased from Bio-Techne. The prostate cancer perineural array and controls were obtained at Baylor College of Medicine. Tissue sections from four patients with prostate cancer were obtained from theNCT01520441phase 1/2 trial designed as a proof of principle that nerves affect the biology of cancer in humans. Patients served as their own controls by receiving BoNT/A injections into the right peripheral and transition zones and sham saline injections into the left peripheral and transition zones2. A paired treated cancer and control was available for only one patient and was used for this study.

Tissue microarray slides and clinical trial slides were dual-stained with antibodies to mitochondria (NeoBiotechnologies, catalogue number MSM2-740-P1ABX) and PGP 9.5 (Invitrogen, PA5-29012) using the Dako Omnis instrument with Envision Flex HRP and high-pH reagents (Agilent, Dako Omnis, GV800). Antigen retrieval was performed using EnV FLEX TRS at high pH (97 °C for 30 min). The PGP 9.5 antibody was used at a dilution of 1:1,500 with a 30-min room-temperature incubation, and the mitochondria antibody was used at 1:5,000 with a 30-min room-temperature incubation. Diaminobenzidine and magenta were used as chromogens for PGP 9.5 and mitochondria, respectively.

Slides with dual immunohistochemical staining were converted to high-resolution digital images through multispectral imaging using a Nuance Multispectral Microscope (PerkinElmer). The images were analysed using inForm analysis software (Akoya Biosciences), as previously described32. The system was trained to recognize and segment cancer, stroma and empty spaces. After segmentation of image compartments, each cell within the identified compartments was analysed separately for nuclear and cytoplasmic regions. Mitochondrial load was quantified by measuring the optical density of magenta staining specifically in the cytoplasm of cancer cells.

Statistical analyses were conducted using Prism 10 for macOS (GraphPad) and Excel for MacOS (Microsoft). Unless otherwise specified in the figure legends, all in vitro data analyses were performed using unpairedt-tests for pairwise comparisons. For clinical data interpretation, high-throughput analysis of mitochondrial load in cancer cells from clinical histology samples was conducted. The equality of variances was first assessed using Levene’s test. A standard Student’st-test was used for comparisons when variances were equal, and Welch’st-test was applied when variances were unequal. For animal studies, unless specified otherwise, data were analysed using ANOVA with Tukey’s multiple comparison test. Two-tailedt-tests were used for exploratory experiments, and one-tailedt-tests were used when a directional effect was supported by prior evidence presented earlier in the study. The number of independent animals used is provided in the figure legends. Shear stress and Seahorse metabolic assays were conducted independently three times, with multiple culture replicates per assay. For Seahorse measurements, 8 to 16 technical replicates were included per assay, depending on the number of biological conditions tested in the 96-well plate. The experiments were analysed using the Wave Controller software (Agilent Technologies) and Prism 10 for MacOS, and graphs shown represent a typical experiment. Microscopy, western blot and PCR experiments were performed at least in triplicate, and representative images are provided. Uncropped western blots and PCR gels are provided in the Supplementary Information.Pvalues are reported as exact values or symbolically as follows: *P< 0.05, **P< 0.01, ***P< 0.001 and ****P< 0.0001. APvalue of less than 0.05 was considered to indicate nominal statistical significance.

Six-week-old virgin female SCID-beige mice were used for the intraductal breast cancer model. The mammary fat pad injection model was performed on BALB/c female mice (age, 5 weeks) and flank injection of B16-F10 in C57BL/6 male mice (age, 6 weeks). Mice were housed under a 12 h/12 h light/dark cycle in a barrier vivarium. Room temperature (22.8 °C ± 1.7 °C) and relative humidity (30–70%) were continuously monitored and maintained. Animals were housed in cages of four, and each experimental condition was assigned to one cage. For each experiment, one cage was randomly selected from a pool of pre-assigned condition-specific cages. Although allocation was not randomized at the individual animal level, random selection at the cage level ensured unbiased assignment. All mice were age- and sex-matched, and housed under identical environmental and handling conditions to minimize potential covariates.

The study was conducted according to the guidelines of the Declaration of Helsinki. All experiments and procedures were conducted in accordance with the guidelines described in the Guide for the Care and Use of Laboratory Animals (National Institutes of Health). Approvals for animal work were obtained from the Baylor College of Medicine Animal Care and Use and Human Subjects Committee and the University of South Alabama IACUC. Informed consent was obtained from all human participants involved in the study.

Further information on research design is available in theNature Portfolio Reporting Summarylinked to this article.

Plasmids generated for this study have been deposited in the Addgene database (greletlab;https://www.addgene.org/Simon_Grelet/). The bioinformatic datasets for the RNA-seq and RNA microarray experiments are available in the Sequence Read Archive (SRA) of the National Center for Biotechnology Information (NCBI) under the accession numbersGSE288069andGSE292157, respectively.Source dataare provided with this paper.

Custom code developed for this study is available via GitHub athttps://github.com/GreletLab/mtDNA-heteroplasmyandhttps://github.com/GreletLab/DAPI-count.

Ayala, G. Neuroepithelial interactions in cancer.Annu. Rev. of Pathol.18, 493–514 (2023).

Coarfa, C. et al. Influence of the neural microenvironment on prostate cancer.Prostate78, 128–139 (2018).

Magnon, C. & Hondermarck, H. The neural addiction of cancer.Nat. Rev. Cancer23, 317–334 (2023).

Torborg, S. R., Li, Z., Chan, J. E. & Tammela, T. Cellular and molecular mechanisms of plasticity in cancer.Trends Cancer8, 735–746 (2022).

ArticleCASPubMedPubMed CentralGoogle Scholar

Faubert, B., Solmonson, A. & DeBerardinis, R. J. Metabolic reprogramming and cancer progression.Science368, eaaw5473 (2020).

ArticleCASPubMedPubMed CentralGoogle Scholar

Mosier, J. A., Schwager, S. C., Boyajian, D. A. & Reinhart-King, C. A. Cancer cell metabolic plasticity in migration and metastasis.Clin. Exp. Metastasis38, 343–359 (2021).

DeBerardinis, R. J. & Chandel, N. S. Fundamentals of cancer metabolism.Sci. Adv.2, e1600200 (2016).

ArticlePubMedPubMed CentralGoogle Scholar

Xing, Y., Zhao, S., Zhou, B. P. & Mi, J. Metabolic reprogramming of the tumour microenvironment.FEBS J.282, 3892–3898 (2015).

Silverman, D. A. et al. Cancer-associated neurogenesis and nerve-cancer crosstalk.Cancer Res.81, 1431–1440 (2020).

Reavis, H. D., Chen, H. I. & Drapkin, R. Tumor innervation: cancer has some nerve.Trends Cancer6, 1059–1067 (2020).

ArticleCASPubMedPubMed CentralGoogle Scholar

Ali, S. R., Jordan, M., Nagarajan, P. & Amit, M. Nerve density and neuronal biomarkers in cancer.Cancers14, 4817 (2022).

ArticleCASPubMedPubMed CentralGoogle Scholar

Li, D. et al. High nerve density in breast cancer is associated with poor patient outcome.FASEB Bioadv.4, 391–401 (2022).

ArticleCASPubMedPubMed CentralGoogle Scholar

Liebig, C., Ayala, G., Wilks, J. A., Berger, D. H. & Albo, D. Perineural invasion in cancer.Cancer115, 3379–3391 (2009).

Mauffrey, P. et al. Progenitors from the central nervous system drive neurogenesis in cancer.Nature569, 672–678 (2019).

Kappos, E. A. et al. Denervation leads to volume regression in breast cancer.J. Plast. Reconstr. Aesthet. Surg.71, 833–839 (2018).

Behbod, F. et al. An intraductal human-in-mouse transplantation model mimics the subtypes of ductal carcinoma in situ.Breast Cancer Res.11, R66 (2009).

ArticlePubMedPubMed CentralGoogle Scholar

Grelet, S. et al. TGFβ-induced expression of long noncoding lincRNA Platr18 controls breast cancer axonogenesis.Life Sci. Alliance5, e202101261 (2022).

Ding, Y. et al. Semaphorin 4F as a critical regulator of neuroepithelial interactions and a biomarker of aggressive prostate cancer.Clin. Cancer Res.19, 6101–6111 (2013).

ArticleCASPubMedPubMed CentralGoogle Scholar

Wang, H. et al. Semaphorin4F is a potential biomarker for clinical progression and prognosis in gastric cancer.Int. J. Clin. Exp. Pathol.16, 210–224 (2023).

CASPubMedPubMed CentralGoogle Scholar

Agostini, M. et al. Metabolic reprogramming during neuronal differentiation.Cell Death Differ.23, 1502–1514 (2016).

ArticleCASPubMedPubMed CentralGoogle Scholar

Zheng, X. et al. Metabolic reprogramming during neuronal differentiation from aerobic glycolysis to neuronal oxidative phosphorylation.eLife5, e13374 (2016).

Borcherding, N. & Brestoff, J. R. The power and potential of mitochondria transfer.Nature623, 283–291 (2023).

ArticleCASPubMedPubMed CentralGoogle Scholar

Nakai, R. et al. Mitochondria transfer-based therapies reduce the morbidity and mortality of Leigh syndrome.Nat. Metab.6, 1886–1896 (2024).

Brestoff, J. R. et al. Intercellular mitochondria transfer to macrophages regulates white adipose tissue homeostasis and is impaired in obesity.Cell Metab.33, 270–282 (2021).

Kidwell, C. U. et al. Transferred mitochondria accumulate reactive oxygen species, promoting proliferation.Elife12, e85494 (2023).

ArticleCASPubMedPubMed CentralGoogle Scholar

Goliwas, K. F. et al. Mitochondrial transfer from cancer-associated fibroblasts increases migration in aggressive breast cancer.J. Cell Sci.136, jcs260419 (2023).

ArticleCASPubMedPubMed CentralGoogle Scholar

Watson, D. C. et al. GAP43-dependent mitochondria transfer from astrocytes enhances glioblastoma tumorigenicity.Nat. Cancer4, 648–664 (2023).

Khozhukhar, N., Spadafora, D., Rodriguez, Y. & Alexeyev, M. Elimination of mitochondrial DNA from mammalian cells.Curr. Protoc. Cell Biol.78, 20.11.1–20.11.14 (2018).

Kujoth, G. C. et al. Mitochondrial DNA mutations, oxidative stress, and apoptosis in mammalian aging.Science309, 481–484 (2005).

King, M. P. & Attardi, G. Human cells lacking mtDNA: repopulation with exogenous mitochondria by complementation.Science246, 500–503 (1989).

Spees, J. L., Olson, S. D., Whitney, M. J. & Prockop, D. J. Mitochondrial transfer between cells can rescue aerobic respiration.Proc. Natl Acad. Sci. USA103, 1283–1288 (2006).

ArticleCASPubMedPubMed CentralGoogle Scholar

Ruder, S. et al. Development and validation of a quantitative reactive stroma biomarker (qRS) for prostate cancer prognosis.Hum. Pathol.122, 84–91 (2022).

ArticleCASPubMedPubMed CentralGoogle Scholar

Sancho, P., Barneda, D. & Heeschen, C. Hallmarks of cancer stem cell metabolism.Br. J. Cancer114, 1305–1312 (2016).

ArticleCASPubMedPubMed CentralGoogle Scholar

Schafer, Z. T. et al. Antioxidant and oncogene rescue of metabolic defects caused by loss of matrix attachment.Nature461, 109–113 (2009).

ArticleCASPubMedPubMed CentralGoogle Scholar

Fiorillo, M., Sotgia, F. & Lisanti, M. P. “Energetic” cancer stem cells (e-CSCs): a new hyper-metabolic and proliferative tumor cell phenotype, driven by mitochondrial energy.Front. Oncol.8, 677 (2019).

ArticlePubMedPubMed CentralGoogle Scholar

Migaud, M. E., Ziegler, M. & Baur, J. A. Regulation of and challenges in targeting NAD+metabolism.Nat. Rev. Mol. Cell Biol.25, 822–840 (2024).

Huang, Q. et al. Fluid shear stress and tumor metastasis.Am. J. Cancer Res.8, 763–777 (2018).

CASPubMedPubMed CentralGoogle Scholar

Barnes, J. M., Nauseef, J. T. & Henry, M. D. Resistance to fluid shear stress is a conserved biophysical property of malignant cells.PLoS ONE7, e50973 (2012).

ArticlePubMedPubMed CentralGoogle Scholar

Piskounova, E. et al. Oxidative stress inhibits distant metastasis by human melanoma cells.Nature527, 186–191 (2015).

ArticleCASPubMedPubMed CentralGoogle Scholar

Dongre, A. & Weinberg, R. A. New insights into the mechanisms of epithelial-mesenchymal transition and implications for cancer.Nat. Rev. Mol. Cell Biol.20, 69–84 (2019).

Lehuédé, C., Dupuy, F., Rabinovitch, R., Jones, R. G. & Siegel, P. M. Metabolic plasticity as a determinant of tumor growth and metastasis.Cancer Res.76, 5201–5208 (2016).

Faubert, B. & DeBerardinis, R. J. Analyzing tumor metabolism in vivo.Annu. Rev. Cancer Biol.1, 99–117 (2017).

Valastyan, S. & Weinberg, R. A. Tumor metastasis: molecular insights and evolving paradigms.Cell147, 275–292 (2011).

ArticleCASPubMedPubMed CentralGoogle Scholar

Bailey-Downs, L. C. et al. Development and characterization of a preclinical model of breast cancer lung micrometastatic to macrometastatic progression.PLoS ONE9, e98624 (2014).

ArticlePubMedPubMed CentralGoogle Scholar

Yang, J. et al. Twist, a master regulator of morphogenesis, plays an essential role in tumor metastasis.Cell117, 927–939 (2004).

Ayala, G. E. et al. Cancer-related axonogenesis and neurogenesis in prostate cancer.Clin. Cancer Res.14, 7593–7603 (2008).

Allen, J. K. et al. Sustained adrenergic signaling promotes intratumoral innervation through BDNF induction.Cancer Res.78, 3233–3242 (2018).

ArticleCASPubMedPubMed CentralGoogle Scholar

Jurcak, N. R. et al. Axon guidance molecules promote perineural invasion and metastasis of orthotopic pancreatic tumors in mice.Gastroenterology157, 838–850 (2019).

Ayala, G. E. et al. In vitro dorsal root ganglia and human prostate cell line interaction: redefining perineural invasion in prostate cancer.Prostate49, 213–223 (2001).

Galappaththi, S. L., Katz, B., Howze, P. H., Hoover, G. & Grelet, S. A CRISPR/Cas9-based assay for high-throughput studies of cancer-induced innervation.Cancers15, 2026 (2023).

ArticleCASPubMedPubMed CentralGoogle Scholar

Magnon, C. et al. Autonomic nerve development contributes to prostate cancer progression.Science341, 1236361 (2013).

Shao, J. et al. Autonomic nervous infiltration positively correlates with pathological risk grading and poor prognosis in patients with lung adenocarcinoma.Thorac. Cancer7, 588–598 (2016).

ArticleCASPubMedPubMed CentralGoogle Scholar

Huang, D. et al. Nerve fibers in breast cancer tissues indicate aggressive tumor progression.Medicine93, e172 (2014).

ArticlePubMedPubMed CentralGoogle Scholar

Albo, D. et al. Neurogenesis in colorectal cancer is a marker of aggressive tumor behavior and poor outcomes.Cancer117, 4834–4845 (2011).

Krstic, J., Trivanovic, D., Jaukovic, A., Santibanez, J. F. & Bugarski, D. Metabolic plasticity of stem cells and macrophages in cancer.Front. Immunol.8, 939 (2017).

ArticlePubMedPubMed CentralGoogle Scholar

Ildiz, E. S., Gvozdenovic, A., Kovacs, W. J. & Aceto, N. Travelling under pressure – hypoxia and shear stress in the metastatic journey.Clin. Exp. Metastasis40, 375–394 (2023).

ArticlePubMedPubMed CentralGoogle Scholar

Fischer, G. M. et al. Molecular profiling reveals unique immune and metabolic features of melanoma brain metastases.Cancer Discov.9, 628–645 (2019).

ArticleCASPubMedPubMed CentralGoogle Scholar

Srinivasan, E. S., Deshpande, K., Neman, J., Winkler, F. & Khasraw, M. The microenvironment of brain metastases from solid tumors.Neuro-oncol. Adv.3, v121–v132 (2021).

Bray, N. Many makes of mitochondria.Nat. Rev. Neurosci.20, 645–645 (2019).

Walker, T. L. & Kempermann, G. One mouse, two cultures: isolation and culture of adult neural stem cells from the two neurogenic zones of individual mice.J. Vis. Exp.84, e51225 (2014).

Kaduri, M. et al. Targeting neurons in the tumor microenvironment with bupivacaine nanoparticles reduces breast cancer progression and metastases.Sci. Adv.7, eabj5435 (2021).

ArticleCASPubMedPubMed CentralGoogle Scholar

Spadafora, D., Kozhukhar, N., Chouljenko, V. N., Kousoulas, K. G. & Alexeyev, M. F. Methods for efficient elimination of mitochondrial DNA from cultured cells.PLoS ONE11, e0154684 (2016).

ArticlePubMedPubMed CentralGoogle Scholar

Malik, A. N., Czajka, A. & Cunningham, P. Accurate quantification of mouse mitochondrial DNA without co-amplification of nuclear mitochondrial insertion sequences.Mitochondrion29, 59–64 (2016).

Raman, I. M. & Bean, B. P. Resurgent sodium current and action potential formation in dissociated cerebellar Purkinje neurons.J. Neurosci.17, 4517–4526 (1997).

ArticleCASPubMedPubMed CentralGoogle Scholar

Gutierrez, C., Cox, C. L., Rinzel, J. & Sherman, S. M. Dynamics of low-threshold spike activation in relay neurons of the cat lateral geniculate nucleus.J. Neurosci.21, 1022–1032 (2001).

ArticleCASPubMedPubMed CentralGoogle Scholar

Farach, A. et al. Neuronal trans-differentiation in prostate cancer cells.Prostate76, 1312–1325 (2016).

ArticleCASPubMedPubMed CentralGoogle Scholar

Chang, J. T. & Nevins, J. R. GATHER: a systems approach to interpreting genomic signatures.Bioinformatics22, 2926–2933 (2006).

Subramanian, A. et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles.Proc. Natl Acad. Sci. USA102, 15545–15550 (2005).

ArticleCASPubMedPubMed CentralGoogle Scholar

Kanehisa, M., Goto, S., Furumichi, M., Tanabe, M. & Hirakawa, M. KEGG for representation and analysis of molecular networks involving diseases and drugs.Nucleic Acids Res.38, D355–D360 (2010).

Davletov, B., Bajohrs, M. & Binz, T. Beyond BOTOX: advantages and limitations of individual botulinum neurotoxins.Trends Neurosci.28, 446–452 (2005).

Conforti, L., Gilley, J. & Coleman, M. P. Wallerian degeneration: an emerging axon death pathway linking injury and disease.Nat. Rev. Neurosci.15, 394–409 (2014).

Sayols, S. rrvgo: a Bioconductor package for interpreting lists of Gene Ontology terms.microPubl. Biol.https://doi.org/10.17912/micropub.biology.000811(2023).

This study was made possible by the support and funding provided by the Frederick P. Whiddon College of Medicine at the University of South Alabama and its Department of Biochemistry and Molecular Biology, and the Mitchell Cancer Institute. We thank P. Howze, B. Moore, A. Cioroch and A. Parton for their technical assistance; A. Jones for assistance with the Illumina runs; the College of Medicine leadership and E. Trepman for editorial support; S. Lloyd for providing the MEF cells; the Flow Cytometry Shared Resource, and S. McLellan, for support in flow cytometry analysis; M. Schuler, D. Miller and the staff of the Department of Comparative Medicine and the Frederick P. Whiddon College of Medicine at the University of South Alabama, for their assistance with animal work, as well as J. Audia and members of the IACUC committee for support in animal care and guidance with the animal models; A. Prakash, M. Migaud, C. Davies and R. Honkanen for discussions about the project; and the Department of Biochemistry and Molecular Biology at Frederick P. Whiddon College of Medicine, University of South Alabama for support in completing this project. Fluorescence microscopy was performed at the Bioimaging Core Research Facilities at the Frederick P. Whiddon College of Medicine, University of South Alabama and at the Center for Advanced Microscopy, a Nikon Center of Excellence, McGovern Medical School, UTHealth Houston (RRID: SCR_025962). Illumina sequencing runs were performed at the Vanderbilt Technologies for Advanced Genomics core facility, supported by CTSA Grant (5UL1 RR024975-03), the Vanderbilt Ingram Cancer Center (P30 CA68485), the Vanderbilt Vision Center (P30 EY08126), and NIH/NCRR (G20 RR030956). The authors are grateful to the patients and their families who donated tissue for research. This research was financed by the National Institutes of Health (NIH) R01 HL140182 to M.T.L., CCTS Pilot Award UL1TR003096-05 to S. Grelet, the Breast Cancer Research Foundation of Alabama (BCRFA2024-2156) to S. Grelet, the Mitchell Endowment and the Patricia Cobb Rodgers Endowment from the Mitchell Cancer Institute to S. Grelet, and startup funds from the Department of Biochemistry and Molecular Biology, Frederick P. Whiddon College of Medicine, University of South Alabama and the Mitchell Cancer Institute to S. Grelet.

Department of Biochemistry and Molecular Biology, Frederick P. Whiddon College of Medicine, University of South Alabama, Mobile, AL, USA

Gregory Hoover, Shila Gilbert, Olivia Curley, William Hixson, Joel F. Andrews & Simon Grelet

Mitchell Cancer Institute, University of South Alabama, Mobile, AL, USA

Gregory Hoover, Shila Gilbert, Olivia Curley, William Hixson, Terry W. Pierce, Joel F. Andrews & Simon Grelet

Department of Microbiology and Immunology, Frederick P. Whiddon College of Medicine, University of South Alabama, Mobile, AL, USA

Department of Physiology and Cell Biology, Frederick P. Whiddon College of Medicine, University of South Alabama, Mobile, AL, USA

Mike T. Lin & Mikhail F. Alexeyev

Department of Pathology and Laboratory Medicine, University of Texas Health Science Center at Houston, McGovern School of Medicine, Houston, TX, USA

Yi Ding, Ping Bu & Gustavo Ayala

Department of Pathology and Laboratory Medicine, University of Kansas Medical Center, Kansas City, KS, USA

Department of Molecular and Cellular Biology, Baylor College of Medicine, Houston, TX, USA

Department of Integrative Biology & Pharmacology, University of Texas Health Science Center at Houston, Houston, TX, USA

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

S. Grelet supervised the study, designed the study and conceived methodological development, and acquired funding. G.A. collaborated in study design and methodological development and supervised pathological interpretations; G.H., S. Gilbert and O.C. assisted in developing the cell–genetic models and other methodological development throughout the study; S. Gilbert, O.C., C.O. and S. Grelet performed the TNBC and melanoma animal work; C.O. and S. Grelet developed the SVZ-NSC cultures and the nerve–cancer coculture model; G.H., Y.D., P.B., G.A. and S. Grelet developed the DRG-based nerve–cancer coculture model; J.F.A., G.A. and S. Grelet assisted with microscopy experiments; G.H., M.F.A. and S. Grelet assisted in the generation of ρ0cells; S. Grelet performed RNA-seq and Nanopore bioinformatics interpretation and wrote the code presented in the study; Y.D. and P.B. assisted in the processing and interpretation of pathology samples; J.T.C. designed and performed microarray bioinformatic and statistical gene-expression analyses in the Botox denervation model of ductal carcinoma in situ; F.B. and D.M. developed the intraductal breast carcinoma model; G.H., M.T.L., G.A. and S. Grelet performed and interpreted electrophysiology experiments and assisted in manuscript editing; and W.H. and T.W.P. assisted in experimental development and troubleshooting. S. Grelet wrote the manuscript, and S. Grelet and G.A. validated it. All authors reviewed the manuscript and approved the final version for publication.

Correspondence toGustavo AyalaorSimon Grelet.

The MitoTRACER technology developed in this manuscript is covered under a pending patent application titled Methods and applications for monitoring mitochondrial transfer between donor and receiver cells, US Patent and Trademark Office, application number 18/885,864. S. Grelet and G.H. are listed as inventors of the technology.

Naturethanks Jonathan Brestoff, Yuan Pan and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.Peer reviewer reportsare available.

Publisher’s noteSpringer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

(A)Ductal carcinoma in situ innervation model to evaluate nerve-cancer dependencies in breast cancer. At 4 weeks after intraductal transplant of ductal carcinoma in situ cells into mice, primary cancers were chemically denervated by BoNT/A treatment or treated with saline control (saline, n = 20 mice; BoNT/A, n = 25 mice).(B)Heatmap of differentially expressed genes between saline and BoNT/A-treated groups. Transcriptomic analysis of microdissected primary cancer cells showed 116 genes were significantly altered: 43 upregulated and 73 downregulated in the BoNT/A-treated group.(C)Gene Set Enrichment Analysis (GSEA) revealed significant downregulation of metabolic pathways, with the tricarboxylic acid (TCA) cycle being the most impacted in cancer cells derived from denervated tumors.(D)Treemap visualization of reduced Gene Ontology (GO) terms generated using the rrvgo72interpretation platform confirmed the regulation of biological processes, predominantly linked to cancer cell metabolism in denervated tumors. Each square of color represents distinct clusters of GO terms grouped by semantic similarity. The size of the squares corresponds to the score of the representative GO term, with higher scores indicating greater significance in the analysis. The list of reduced GO terms obtained through rrvGO was provided in Supplementary TableS3.(E)Quantification of invasive lesions at 10 weeks through microscopy-based pathological assessment by a certified pathologist of tumors that extend beyond the originating duct, into the host stroma, after breaching the basement membrane. A significant reduction in the mean number of invasive lesions was observed in the BoNT/A-treated group compared to the saline group, indicating decreased invasiveness of cancer cells following denervation (mean ± SEM; unpaired t-test; p = 0.0014).a, Created in BioRender. S. Grelet (2025)https://biorender.com/q4jviiq.

(A)Subventricular zone neuronal stem cells (SVZ-NSCs) were isolated from a fresh 6-week-old female BALB/c mouse brain. The dorsal view of the dissected mouse brain is shown (scale bar, 1 cm).(B)Lateral view of the mouse brain showing the translucent SVZ region. The diagram shows the area (orange) from which SVZ-NSCs were extracted.(C)Isolated SVZ-NSCs were cultured on laminin-coated plates. The cells differentiated into mature neurons using a differentiation medium.(D)The SVZ-NSCs were maintained as neurospheres on uncoated plates, confirming stemness (self-renewal).(E)Fluorescence microscopy of SVZ-NSCs transduced with GFP-expressing lentiviruses and cultured independently (SVZGFPAlone) showed moderate neurite development, evidenced by the absence of neurite extensions. SVZ-NSCs cocultured with 4T1 cells (SVZ-NSCGFP+ 4T1) showed a marked increase in SVZ-NSC differentiation, evidenced by extensive neurite outgrowth traversing the coculture. Scale: 100μm.(F)Immunohistochemistry of SVZ-NSCGFPcocultured with 4T1 cells shows SVZ-NSCs expressing the neuronal marker Tubulin β-3 (Tubb3), confirming neuronal differentiation. Scale: 100 μm.(G)Flow cytometry analysis of nerve-cancer cocultures indicates the commitment of neuronal precursors to neurons. Nearly all neurons expressed MAP2 and Tubulin β−3 (TUBB3) but were negative for O4 and ALDH1L1, confirming that these cells were not oligodendrocytes or astrocytes. (Mean ± S.D, Student’s two-tailed unpaired t-test, n = 3).(H)Histogram plot from flow cytometry illustrating the neuronal commitment of SVZ-NSCs in the nerve-cancer coculture. EGFP-expressing SVZ-NSCs, introduced into the coculture and comprising approximately 1.6% of the total cell population in the series, were gated based on eGFP expression to assess SVZ marker expression using antibodies against O4, ALDH1L1, TUBB3, and MAP2. The omission of the primary antibodies in co-culture consisting of eGFP-expressing SVZ-NSCs mixed with 4T1 cancer cells was used as control.(I)Validation of antibodies used to quantify the neuronal commitment of SVZ-NSCs in the nerve-cancer coculture. SVZ-NSCs were treated with triiodothyronine (T3) to induce oligodendrocyte commitment and stained with the O4 antibody. To induce astrocyte differentiation, cells were treated with fetal bovine serum (FBS) and stained with the ALDH1L1 antibody. For neuronal differentiation, cells were treated with differentiation medium and stained with Tubulin β3 (TUBB3) or Microtubule-Associated Protein 2 (MAP2) antibodies. The omission of primary antibodies on SVZ-NSCs culture was used as control.(J)Calcium flux analysis in nerve-cancer coculture using Fluo-4 AM staining. (Top) Time-lapse imaging showed pulsatile calcium activity in neurons. (Bottom) Signal dynamics with 10 μM nifedipine treatment demonstrated a blockade of calcium flux, abolishing pulsatile activity.(K)Current-clamp recording of neurons in nerve-cancer coculture. Whole-cell recordings were used to record the membrane potentials of neurons. Cells clamped at −80 mV, were depolarized with incremental current ramp injections ranging from 20–200 pA, as shown in the inset. Cocultured neurons showed action potential activity.(L)Distribution of action potential threshold measurements between neurons analyzed in the nerve-cancer coculture. Action potentials were elicited in all recorded neurons (mean ± SEM n = 9 neuron cells acquisition; 100% response).

(A)Related to Fig.1d- Oxygen consumption rate (OCR) analysis comparing basal, maximal, and spare respiratory capacities of 4T1 cells alone versus coculture with SVZ-NSCs. Coculture led to a significant increase in all three parameters, suggesting enhanced mitochondrial function (Representative Profile (n = 3); Mean ± SEM; Student’s two-tailed unpaired t-test ****P< 0.0001).(B)MitoSOX mitochondrial indicator staining (red) of SVZ-NSCs before (SVZ-NSC) and after differentiation (NSC Diff.). There was a significant increase in mitochondrial activity observed after differentiation, evidenced by an increase in MitoSOX staining.(C)Cancer-driven neuronal differentiation markedly upregulated mitochondrial load in SVZ-NSCs. PCR of mtDNA content in SVZ-NSCs under varied concentrations of differentiation medium (0% to 80%) showed a progressive increase in mtDNA content throughout differentiation (nDNA, nuclear DNA loading control).(D)FACS analysis following mitochondrial staining (MitoTracker) shows SVZ-NSCs in three conditions: undifferentiated and cultured alone (blue, SVZ-NSCs alone), differentiated with nerve growth factor (red, SVZ-NSCs + NGF), and differentiated by coculture with 4T1 cancer cells (green, SVZ-NSCs in coculture). Representative Profile (n = 3 idenpendent cell cultures).(E)FACS and quantitative polymerase chain reaction (qPCR) of mitochondrial DNA (mtDNA) to nuclear DNA (nDNA) ratio in SVZ-NSCs alone (Alone) or in coculture with 4T1mCherrycells (With 4T1) (mean ± S.D.). There was increased mitochondrial mass in SVZ-NSCs after exposure to 4T1 cells.(F)Fluorescence microscopy of SVZ-NSCs having eGFP-labeled mitochondria (SVZ-NSCCCO-GFP). The mitochondria had globular morphology and were localized around the nucleus. In contrast, SVZ-NSCCCO-GFPcocultured with 4T1mCherrycells had mitochondria with elongated morphology and were dispersed throughout the cell cytoplasm. Scale: 10 μm.

(A)Quantification of mitochondrial transfer. Flow cytometry confirmed strong eGFP fluorescence in 50B11 DRG-derived cells (50B11CCO-GFP). When cocultured with 4T1mCherry+cells, mitochondrial transfer was indicated by the acquisition of eGFP fluorescence in the 4T1mCherry+cells. The gating strategy was established using 50B11CCO-GFPand 4T1-mCherry cells cultured independently.(B)Cell viability of cultures treated with cytochalasin B (Cyto B) for 24 h. The percentage of live cells was calculated using Hoechst/PI staining and automated counting with a fluorescence cytometer (Celigo). No significant change in cell viability was observed after treatment in both SVZ-NSC and 4T1 cancer cells.(C)Quantification of mitochondrial transfer from (Left) SVZ-NSCCCO-GFPor (Right) 50B11 mitochondrial donor cells to cancer cell lines derived from human breast (MDA-MB-231), lung (A549), and prostate (PC3) carcinomas. The normalized transfer rate was calculated as the percentage of recipient cells within the eGFP+ population in the coculture, standardized across biological conditions. Mean ± S.D. (n = 6).(D)Related to Fig.1m. Quantitative PCR analysis of mtDNA acquisition in ρ04T1 cancer cells exposed to SVZ-NSCs and FACS isolated from the co-culture. Mean ± S.D. (n = 3).

(A)Microscopy validation of the Syn1-GFP-NLS and Syn1-GFP-OMP25 lentiviral constructs in PC12 neuronal cells and 4T1 breast cancer cells. Only neurons displayed expression of the synapsin 1-driven constructs, with eGFP localized to the nucleus (localization signal, NLS) or mitochondria (OMP25). DAPI was used as a counterstain for the nuclei. Scale: 50μm.(B)Flow cytometry validated the Syn1-GFP-OMP25 construct and confirmed neuron-specific expression by showing strong expression in PC12 neuronal cells and no activity in 4T1 cancer cells. This analysis also determined the titer of the production batch.(C)Immunohistochemical staining of mouse mammary fat pad xenograft tissue nerves using tubulin beta-3 (Tubb3) and labeling of the eGFP fluorophore confirms the expression of the transgene in the mouse host nerves. 4X magnification of the xenograft area. Zoom: Nerve fiber at 40X magnification. Magenta: tubulin beta-3, Green: eGFP, Blue: DAPI. Scale: 500 μm.

(A)Validation of the cloning constructs. We tested the plasmid constructs by co-transfection of MitoTRACER with increasing amounts of a separated TEVp-encoding plasmid into 4T1 cells stably expressing the LoxP-DsRed Express2-STOP-LoxP-eGFP recipient genetic switch but lacking the TEVp construct. Fluorescence microscopy showed activation of the Red-to-Green fluorophore switch when all components of the system were expressed in the cells by co-transfection with the TEVp. This confirmed the Red-to-Green switch and the requirement of TEVp to ensure feasibility. Scale: 100 μm.(B)Time-lapse fluorescence microscopy over 12 h of SVZ-MitoTRACER mitochondrial donor cells cocultured with 4T1 SWITCH-TEVp recipient cells, illustrating the dynamics of mitochondrial transfer. The process included the establishment of close cell-cell contacts, the development of tunneling nanotube structures, and subsequent red-to-green fluorescence transition in recipient cells. Arrowheads: blue recipient cells; yellow, tunneling nanotube structure. Captures from Video S3.

Gating strategies in flow cytometry experiments.a, Flow cytometry gating strategy used in the study for the mitochondria transfers using the CCO-GFP or OMP25-GFP mitochondrial reporters.b, Flow cytometry gating strategy for transfer experiments using the MitoTRACER approach to determine the red-to-green ratio.

Time-lapse microscopy of 4T1mCherrycells cocultured with SVZ-NSCsCCO-GFP. Confocal imaging of the coculture showing the transfer of green mitochondria from SVZ-NSCCCO-GFPdonor cells to 4T1mCherrycancer cells. Images were acquired using confocal microscopy, withz-stacking adjusted to focus on the transferred mitochondria in the recipient cell as presented in Fig. 1e.

Time-lapse microscopy showing spontaneous calcium pulsatile activity in the neuron of the coculture. Related to Extended Data Fig. 2j.

Time-lapse microscopy of mitochondrial transfer dynamics. Time-lapse microscopy of mitochondrial transfer dynamics in a nerve–cancer coculture consisting of SVZ-NSCMitoTRACERmitochondrial donor cells and 4T1SWITCH-TEVprecipient cells. Blue arrows indicate the recipient cells, and yellow arrows highlight the development of tunnelling nanotube structures originating from the donor cells, followed by the rapid red-to-green fluorescence transition in the recipient cells. This underscores the important function of cell–cell contact in facilitating mitochondrial transfer.

MitoTRACER analysis of nerve–cancer mitochondrial transfer. Time-lapse microscopy over 48 h showing the dynamics of mitochondrial transfer, highlighted by the red-to-green fluorescence transition in the coculture of neuron-derived mitochondrial donor cells (SVZ-NSCsMitoTRACER) and 4T1SWITCH-TEVprecipient cells.

MitoTRACER analysis of fibroblast-to-cancer mitochondrial transfer. Time-lapse microscopy over 48 h showing the dynamics of mitochondrial transfer, marked by the red-to-green fluorescence transition in the coculture of MEF-derived mitochondrial donor cells (MEFMitoTRACER) and 4T1SWITCH-TEVprecipient cells.

MitoTRACER analysis of adipocyte-to-cancer mitochondrial transfer. Time-lapse microscopy over 48 h showing the dynamics of mitochondrial transfer, highlighted by the red-to-green fluorescence transition in the coculture of pre-adipocyte-derived mitochondrial donor cells (3T3-L1MitoTRACER) and 4T1SWITCH-TEVprecipient cells. After stimulation by cancer cells, 3T3-L1 cells exhibited changes in morphology, including increased size and granularity, suggesting differentiation to adipocytes in the coculture.

Open AccessThis article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visithttp://creativecommons.org/licenses/by-nc-nd/4.0/.

Hoover, G., Gilbert, S., Curley, O.et al.Nerve-to-cancer transfer of mitochondria during cancer metastasis.Nature(2025). https://doi.org/10.1038/s41586-025-09176-8

DOI:https://doi.org/10.1038/s41586-025-09176-8

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

Spin-qubit control with a milli-kelvin CMOS chip

A key virtue of spin qubits is their sub-micron footprint, enabling a single silicon chip to host the millions of qubits required to execute useful quantum algorithms with error correction1,2,3. However, with each physical qubit needing multiple control lines, a fundamental barrier to scale is the extreme density of connections that bridge quantum devices to their external control and readout hardware4,5,6. A promising solution is to co-locate the control system proximal to the qubit platform at milli-kelvin temperatures, wired up by miniaturized interconnects7,8,9,10. Even so, heat and crosstalk from closely integrated control have the potential to degrade qubit performance, particularly for two-qubit entangling gates based on exchange coupling that are sensitive to electrical noise11,12. Here we benchmark silicon metal-oxide-semiconductor (MOS)-style electron spin qubits controlled by heterogeneously integrated cryo-complementary metal-oxide-semiconductor (cryo-CMOS) circuits with a power density sufficiently low to enable scale-up. Demonstrating that cryo-CMOS can efficiently perform universal logic operations for spin qubits, we go on to show that milli-kelvin control has little impact on the performance of single- and two-qubit gates. Given the complexity of our sub-kelvin CMOS platform, with about 100,000 transistors, these results open the prospect of scalable control based on the tight packaging of spin qubits with a ‘chiplet-style’ control architecture.

Utility-scale quantum computing probably requires millions of physical qubits, operated by auxiliary classical systems that generate more than a trillion control signals per second13,14. In realizing this vast and complex platform, silicon qubits present advantages with their small footprint15, long coherence times16and inherent compatibility with VLSI (very large scale integrated) control circuits. Although the potential for integrated control has been a key motivator for the progress in silicon-based qubits over the past two decades, so far this aspect has remained largely undeveloped.

Despite the advantages of integrated control4,5, a serious concern arises from the heat and crosstalk generated by modern complementary metal-oxide-semiconductor (CMOS) circuits. In relation to heat, this problem is eased by recent work showing that spin qubits continue to function at elevated temperatures17,18,19. But two-qubit entangling gates remain sensitive to electrical noise11,12, arising, for instance, from volt-scale, sub-nanosecond switching of proximal CMOS transistors. One way of partially mitigating these adverse effects is to separate the control system to 4 K, connecting to milli-kelvin qubits using long cables8,20. Cable connectivity poses an additional barrier to scaling up the control interface4, given the extreme density of interconnects required to operate even modest numbers of qubits.

Here we demonstrate the control of MOS-style silicon spin qubits using a heterogeneously integrated cryo-CMOS chip operating at milli-kelvin temperatures, as shown in Fig.1a. Heterogeneous, ‘chiplet-style’ integration, as opposed to monolithic circuits, decouples the hot and noisy control system from sensitive qubits and retains the potential for dense, lithographically defined chip-to-chip interconnects needed to manage the wiring challenge inherent to spin qubits. We demonstrate that this chiplet architecture supports a control scheme that leverages a global resonance field to enable complete universal control of spin qubits using the baseband pulses that can be generated efficiently with proximal, low-power cryo-CMOS.

a, The cryo-CMOS and the qubit chip are mounted on the same circuit board at milli-kelvin temperatures and wire-bonded together. All other control systems are at room temperature.b, An electron micrograph of a nominally identical silicon device to that measured here. There are 32 CLFG cells on the chip, one of which is connected to gate J that controls the coupling between the two quantum dots on the qubit chip. The other cell is connected to gate B, which acts as an additional barrier gate. All other fast-pulse gates (P1, P2 and SET) and the d.c. barrier gates are connected to room-temperature electronics.c, Schematic of a cell and its electrical connection to gate J of the silicon device. Pulsing on this gate acts to modulate the tunnel coupling between the two quantum dots needed for single- and two-qubit controls. Scale bar, 100 nm (b).

The details of the CMOS control chip have been reported previously with an early conceptual demonstration using GaAs quantum dot structures7. The effect of milli-kelvin CMOS on qubit performance, however, has remained an open question until the present work. As the spin degree of freedom is decoupled from electrical noise, integrated CMOS is expected to have only a minor impact on single-qubit operations. By contrast, coupling spins by Heisenberg exchange creates the most sensitive probe of voltage noise known11,12because the exchange energy can depend exponentially on gate voltage. Countering this intuition, we show that even for noise-sensitive two-qubit gates, our chiplet architecture, comprising some 100,000 transistors, does not lead to a measurable reduction in coherence time.

An electron micrograph of a silicon-metal-oxide-semiconductor qubit device is shown in Fig.1b. The device is fabricated on an isotopically purified28Si epilayer with a residual29Si concentration of 800 ppm (ref.21) and SiO2isolating layer with metal gates patterned in aluminium. Quantum dots hosting single spin qubits are formed under the plunger gates (P1 and P2) at the Si/SiO2interface, and an exchange gate (J) modulates the tunnel coupling between the two dots, essential for two-qubit operations. A radiofrequency single-electron transistor (RF-SET)22detects the charge state of the quantum dots on microsecond timescales by leveraging an off-chip LC resonator operating near 400 MHz (ref.23), and a proximal microwave antenna generates an oscillating magnetic field for spin resonance control (Fig.1b).

The exchange gate J and a barrier gate B are wire-bonded to the cryo-CMOS control chip7, which is implemented in 28 nm fully depleted silicon-on-insulator technology (FDSOI). The chip contains a serial peripheral interface for handling digital input instructions and a finite state machine (FSM) for on-chip digital logic. The FSM configures 32 analogue ‘charge-lock fast-gate’ (CLFG) circuit blocks, each of which can be used to control a gate electrode on the quantum device (Fig.1c). In this configuration, the charge is periodically stored and shuffled between small capacitors, leveraging the low leakage of transistors at cryogenic temperatures that maintain the potential during quantum operations. The cryo-CMOS chip also incorporates a ring oscillator and configurable register designed as a programmable internal trigger (see Extended Data Fig.3afor oscillator schematics). Here, for convenience, we opt for external triggering.

The core functionality of a CLFG cell is to lock a static voltage bias and enable a fast pulse between two voltage levels, as outlined in Fig.1c. For instance, targeting gate J, the CLFG cell is programmed to first bring the gate to a potentialVout, equal to the potentialVholdof an external source. Opening the switchGlockunder the control of the FSM ‘charge locks’ this potential on the gate capacitor. Although this floating capacitor is now galvanically disconnected from the source, a pulse can be induced by toggling the potential on the top plate of this capacitance betweenVhighandVlow, as shown in Fig.1c. This toggling is produced autonomously by the programmed on-chip FSM, leading to a modified outputVoutby ΔVpulse= (Cpulse/CP+Cpulse) × (Vhigh−Vlow), whereCPis the parasitic capacitance. This mechanism has previously been shown to produce pulse amplitudes of 100 mV at a power of about 20 nW MHz−1(ref.7). Below, we demonstrate how this architecture can be used to efficiently control spin qubits.

To evaluate single-qubit gates, we first establish a baseline using all room-temperature electronics for control. Following the usual protocol for two-spin manipulation1,24,25, the singlet state is first prepared in the (1, 3) charge configuration using pulses applied to detuning gates P1 and P2, with (n,m) labelling the number of electrons in each dot under P1 and P2, respectively. A pulse applied to the J gate, connected toVhold, then increases the barrier, separating the two electrons into each dot, in which they are independently addressed using the microwave antenna through their unique resonance frequency (fESR= 13.9 GHz for a fieldB0= 0.5 T). Free-induction decay (FID) of the target spin is produced by applying microwave power to the on-chip electron spin resonance (ESR) line. Finally, a second pulse of the J gate returns the spins to the readout configuration, in which Pauli spin blockade enables spin-to-charge conversion26,27and measurement by the RF-SET. The shot-averaged readout signal as a function of microwave pulse time and frequency is shown in Fig.2a. Beyond FID, we further establish our room temperature baseline by performing pulse sequences implementing Hahn echo (to measure coherence timeT2; Extended Data Fig.4) and randomized benchmarking (to measure qubit control fidelity)28,29. We have intentionally limited the number of qubit gates bonded to the CMOS control chip to facilitate direct comparison with room temperature control in the same cooldown. The present prototype CMOS chip contains 32 CLFG cells to drive 32 qubit gate electrodes.

a, Single-qubit Rabi oscillations (Q1) as a function of microwave (MW) frequencyfMWand pulse timetMW, performed with room temperature (RT) control.b, Single-qubit randomized benchmarking (Q1) under various cryo-CMOS conditions. Traces are offset for clarity. Each data point is the average of 300 randomized sequences at 100 shots each.c, Single-qubit\({T}_{2}^{* }\)coherence time as a function of select cryo-CMOS parameters. Unless otherwise indicated, all data use cryo-CMOS control. Each data point is the average of 100 shots with 4 repeats for a total of 400 single shots.d, Mixing chamber temperature with cryo-CMOS power. Error bars represent the 95% confidence level. Osc., oscillator.

In this single-qubit measurement, the function of the J-gate pulse is to separate the two-spin system for controlled rotation by spin resonance. As such, electrical noise, coupled through the J gate or other means, is unlikely to affect qubit fidelity in the limit that the pulse amplitude and duration are sufficiently large to fully separate the spins. Even so, we now evaluate the impact of cryo-CMOS control on single-qubit performance by performing the same protocol outlined above, but now with charge-locking applied to the J gate and the pulse produced using a CLFG cell under control of the FSM. Again, we generate FID data and quantitatively compare the CMOS and room temperature control using randomized benchmarking protocols, as shown in Fig.2b. A slight degradation in qubit fidelity is observed (0.07%), probably because of unmitigated heat from the CMOS. We discuss heating in detail below.

Although electrical noise at the J gate does not directly couple with single spins, heat and drift in gate potential over longer timescales can affect qubit performance. Gate noise can also produce d.c. Stark shift of the qubit frequency in certain regimes (discussed further below). To investigate these mechanisms, we extract the time-ensemble average coherence time\({T}_{2}^{* }\)for each qubit, repeatedly measured as each circuit block of the CMOS chip is powered up. Comparing the data in Fig.2cagain shows a small impact with respect to our room temperature baseline, correlating with a slight rise in the base temperature of the refrigerator (see Fig.2dand its caption for details).

A drawback of the control scheme outlined above is its reliance on qubit-specific microwave pulses, which, despite cryo-CMOS gate control, require additional room-temperature microwave generators (and cables) for each qubit. We next demonstrate an alternate control approach that leverages a continuous wave global microwave field, sourced from room temperature but common to all qubits15,30. Key to this scheme is the ability to tune the spin resonance frequency of a qubit using a gate voltage31,32. This d.c. Stark shift enables a gate pulse to bring a qubit into resonance with the global field for a controlled amount of time to produce a rotation in the qubit state vector. With a single microwave tone from room temperature, cryo-CMOS produces the ‘baseband’ gate pulses that independently bring each qubit into and out of resonance with the global field.

The global microwave scheme requires a calibration of the unique d.c.-Stark shift produced by the J gate on, for example, qubitQ2, as shown in Fig.3a. The sizable shift (about 1 MHz per 10 mV) is well-matched to the voltage pulses that can be efficiently generated with proximal low-power CMOS. Here, the spins are initialized to a mixed orT−(|↓↓⟩) state off-resonance by a detuning pulse to a relaxation hot-spot33,34. The J-gate pulse produces the time-controlled Stark shift as shown in Fig.3b. Repeating this sequence as a function of pulse length yields the coherent oscillations (Fig.3c), and coherence metrics can be extracted (Fig.3d). For this measurement, the width of the pulse is set by the timing of the trigger fed to the CMOS from room temperature. Conceptually, this reliance on room temperature triggering may seem to be a limitation of our CMOS circuits. However, we note that fine time resolution is needed only to map out coherent oscillations. By contrast, once calibrated, logic gates require a fixed time pulse, for instance, 1.034 μs to produce a π/2 rotation. As such, these fixed time pulses are straightforward to implement with our CMOS architecture. Furthermore, we note that high fidelity control is possible with fixed time width pulses by precise tuning of the pulse amplitude and bias, shifting the requirement for high-resolution timing to the resolution of the voltage source.

a, ESR spectra as a function of microwave (MW) frequency andJgate voltage.Q2exhibits a significant Stark shift, which allows for on-resonance single-qubit operations (X,Z; triangle), and off-resonance loading and measurement (M; star).b, Schematic of the pulsing sequence when using effective global control. Qubit rotations are determined byJpulse time, with the microwave tone fixed to the maximumtpulsetime, extending into load and read stages.c, Rabi chevron ofQ2.d, Ramsey coherence time, using global control and cryo-CMOS pulsing atB0= 0.5 T. Error margin represents the 95% confidence level. a.u., arbitrary units.

Finally, we turn to evaluate two-qubit logic gates, which provide the most stringent test for crosstalk or electrical noise stemming from proximal milli-kelvin CMOS control. Here, the target qubit is rotated about thez-axis depending on the state of the control qubit, with the gate-tunable exchange interaction modulating the coupling between the two electrons1,24. As a baseline, we first perform a decoupled controlled phase gate (DCZ) using all room temperature instrumentation. The DCZ gate incorporates a spin-echo sequence with coherent rotation about thez-axis of angleϕ=J(ϵ)texchangeħ, enabled by turning on exchange for a controlled time with J-gate pulse of durationtexchange(Extended Data Fig.2d). The resulting readout probability with J-gate pulse width is shown in Fig.4a. The DCZ gate is sensitive to high-frequency electrical noise arising either directly from exchange-gate voltage fluctuations or indirectly from noise in the (gradient) magnetic field or variation ing-factors from gate-induced movement in the position of the electron wavefunction.

a, DCZ oscillations as a function of exchange timetexchangeandVJ, performed using room-temperature (RT) control.b, Room-temperature- and cryo-CMOS-enabled DCZ oscillations at a setVJ, indicated ina. Set level is determined byVhold, which can be tuned for stronger qubit interaction.c, Comparison of visibility between room-temperature and cryo-CMOS control using identical pulsing methods on all gates. Increased state preparation and measurement error is present because of two-level-only approach to pulsingVJ.d, CZ coherence times of our two-qubit control conditions under various cryo-CMOS parameters.Jpulses are generated by cryo-CMOS unless otherwise indicated. Each data point is the average of 100 shots, with 4 repeats for a total of 400 single shots. Error bars represent the 95% confidence level. Osc., oscillator.

Comparing cryo-CMOS control with our room-temperature baseline, we observe that the coherent exchange oscillations show similar behaviour (Fig.4b). These results immediately confirm the utility of proximal milli-kelvin CMOS for controlling two-qubit logic gates. Close inspection perhaps suggests a suppression in visibility for the CMOS data, which probably stems from the limited two-state resolution of the voltage pulses used to tune the readout and preparation state, which for this qubit device require significantly different tunnel rates than those used in two-qubit control. Although it is not uncommon to find tunnel rates that are very similar for qubit control as state preparation and measurement (SPAM), in the present device, the differing tunnel rates precluded more quantitative measures of SPAM error and two-qubit fidelity using cryo-CMOS35. Future improvements in qubit tunability and fidelity will also enhance sensitivity to new noise sources, including further assessment of the control platform.

As a noise diagnostic tool, it is also worth noting that the DCZ gate is limited because the spin-echo sequence decouples the spin dynamics from low-frequency noise. Removing the echo pulses then opens the bandwidth to now include all of the low-frequency components down to quasi-d.c., offering a better measure of the total aggregate noise inherent in the system. A comparison of room-temperature control and cryo-CMOS is made in the data shown in Fig.4c, now without spin-echo. These datasets constitute a measure of the ensemble average coherence time associated with the exchange gate,\({T}_{2,{\rm{CZ}}}^{* }\). Finally, using this parameter as a wideband measure of noise, Fig.4dcompares\({T}_{2,{\rm{CZ}}}^{* }\)for room-temperature control, cryo-CMOS with a single charge-locked cell, all 32 cells locked (mirroringJ-gate pulses), and as a function of CMOS oscillator frequency. A slight reduction in\({T}_{2,{\rm{CZ}}}^{* }\)(around 20%) is observed at the highest clock frequencies, which, given the slight corresponding increase in refrigerator temperature, can be explained as arising from parasitic heating (for more detailed two-qubit performance data see Extended Data Fig.5). We note that no additional (electrical) noise is observed beyond the thermal noise contribution associated with the small increase in temperature from CMOS power dissipation (for further discussion, see Extended Data Fig.4).

Our cryo-CMOS control chip consists of complex mixed-signal circuits realized using more than 100,000 transistors. Most of these transistors are used in the digital sub-systems and related circuit blocks, accounting for a fixed overhead power of tens of microwatts. On top of this constant offset power from the digital blocks, the CLFG analogue cells each contribute approximately 20 nW MHz−1when generating 100 mV amplitude pulses, enabling many thousands of cells (and thus gate pulses) to fit within the cooling budget of a commercial dilution refrigerator (around 1 mW at 100 mK) (ref.7). Apart from the cooling limits of the refrigerator, however, a challenge arises in the thermal management of hot control systems to ensure the routing of heat bypasses proximal, cold quantum devices. Here, we have made no attempt to mitigate this parasitic heating, simply wire-bonding the chips together in a standard package. This arrangement can lead to elevated electron temperatures in the quantum device (Extended Data Fig.3) even when the refrigerator remains cold (Fig.2d) and is the likely explanation for the small impact we observe in qubit fidelity when the largest CMOS circuits are powered up at the highest clock rates. As such, we emphasize that there is a notable opportunity to suppress parasitic heating by using separate parallel cooling pathways for the CMOS chip and quantum plane7. The use of heterogeneous, rather than monolithic, integration opens new thermal configuration options in this regard.

Beyond direct heating, the close presence of 100,000 transistors, with volt-scale biasing and sub-nanosecond rise and fall times, can create an exceedingly noisy environment in which to operate electrically sensitive qubits. It is surprising that the CMOS chip has only a small impact on qubit performance relative to previous experiments with room-temperature control17,19. Furthermore, the small degradation in fidelity is probably explained entirely from parasitic heating, rather than from electrical noise. Certainly, our use of CMOS design rules that minimize external crosstalk are important; however, beyond these, we suggest three additional aspects that probably reduce electrical noise. First, as the physical temperature of the CMOS die is a few hundred milli-kelvin, thermal noise contributions are substantially suppressed. Second, the chip-to-chip interconnect probably has a relatively low bandwidth, filtering noise above a few gigahertz. Last, we note that the action of the CLFG circuits effectively decouples the CMOS from the quantum device when in charge-lock mode, except for a very small coupling capacitor. Taken collectively, these aspects further underscore the utility of heterogeneous over monolithic integration for mitigating crosstalk and heating. Apart from addressing the challenges posed by scaling up qubits, cryo-CMOS using a chiplet architecture may also prove useful in generating ultrafast, low thermal noise control pulses that probe fundamental physics in mesoscale quantum devices36.

In conclusion, the results presented here demonstrate the viability of heterogeneous, milli-kelvin CMOS for generating the volt-scale biases and milli-volt pulses needed to control spin qubits at scale. Beyond addressing the interconnect bottleneck posed by cryogenic qubit platforms, these results show that degradation in qubit performance from milli-kelvin CMOS is very limited. Although our focus here has been controlling spin qubits based on single electrons, we draw attention to the inherent compatibility of our control architecture with other flavours of spin qubits, for instance, exchange-only qubits that leverage square voltage pulses exclusively37. Pairing cryo-CMOS-based control with highly compatible radiofrequency readout approaches that exploit dense frequency multiplexing38enables a highly integrated and scalable spin qubit platform.

Measurements are performed in a Bluefors LD400 dilution refrigerator with a base temperature of 7 mK. The qubit chip and cryo-CMOS chip are packaged on the same FR4 printed circuit board (PCB)39, separated by 3 mm and wire-bonded together. The daughterboard PCB is placed on a custom motherboard and electrically connected by an interposer. The motherboard manages signals from room temperature, and the ESR line connects directly to the daughterboard through a miniSMP. The PCB setup is mounted in a magnetic field, on a cold finger. The external magnetic field is supplied by an Oxford MercuryiPS-M magnet power supply. The d.c. voltages are generated by an in-house custom-made digital-to-analogue converter (DAC). A Quantum Machines OPX+ generates room-temperature dynamic pulses, the 400 MHz signal for radiofrequency readout, I/Q and pulse modulation for the microwave source as well as trigger lines to the cryo-CMOS chip and microwave source (Extended Data Fig.1). Dynamic and d.c. voltage sources are combined at room temperature using custom-made voltage combiners. The OPX has a sampling rate of 1 GS s−1and a clock rate of 250 MHz. The microwave tone is generated by a Keysight PSG E8267D Vector Signal Generator with a signal spanning 250 kHz to 31.8 GHz. Single-qubit gates are operated at a microwave frequency of 13.9 GHz.

We use an RF-SET and radiofrequency reflectometry readout, comprising a Low Noise Factory LNF-LNC0.3-14A amplifier at the 4 K stage of the refrigerator and a Minicircuits ZX60-P103LN+ at room temperature for signal amplification. The directional coupler in Extended Data Fig.1is a Minicircuits ZEDC-10-182-S+ 10–1,800 MHz.

All measurements are performed in the same cooldown to enable careful comparison between control methods without device variation that can occur from thermal cycling. For room-temperature operation, a pass-through switchGlockis closed and room-temperature-sourced exchange gate pulses are delivered throughVhold(see Extended Data Fig.1for instrument connection details).

Our cryo-CMOS receives programming instructions, power, d.c. bias and dynamic voltage levels from an in-house room-temperature DAC40and is configured to receive an external trigger from the OPX+, as well as room-temperature dynamic pulses for pass-through room-temperature control. Some of the programming instructions as well as the external trigger, are received by the same input line. Appropriate input is handled by a Minicircuits RC-4SPDT-A18 DC-18 GHz radiofrequency switch, which takes inputs from both the DAC and OPX+. The radiofrequency switch is programmed to work in unison with charge-locking commands sent to the cryo-CMOS, switching inputs from the DAC to the OPX+ once programming is complete.

Owing to the low leakage rate of our cryo-CMOS transistors, there are no charge-lock refresh cycles of our charge-locked gate when performing experiments. Maximum experiment times are approximately 1 h, usually single-qubit RBM or PSD measurements. Cryo-CMOS circuit blocks were measured to be functional within the bias range DD1P0 = 0.65V, N/PMOS backgate = 1.5/−1.5 V, VDD1p8 = 1.2 V–1.6 V.

The datasets generated and/or analysed relevant to this study are available from the corresponding author upon request and further available at Zenodo (https://zenodo.org/records/15080656)41.

Analysis codes supporting the findings of this study are available from the corresponding authors upon request and further available at Zenodo (https://zenodo.org/records/15080656)41.

Loss, D., & DiVincenzo, D. P. Quantum computation with quantum dots.Phys. Rev. A57, 120–126 (1998).

Maurand, R. et al. A CMOS silicon spin qubit.Nat. Commun.7, 13575 (2016).

Veldhorst, M., Eenink, H. G. J., Yang, C. H., & Dzurak, A. S. Silicon CMOS architecture for a spin-based quantum computer.Nat. Commun.8, 1766 (2017).

Reilly, D. J. Challenges in scaling-up the control interface of a quantum computer. In2019 IEEE International Electron Devices Meeting (IEDM), 31.7.1–6 (2019).

Vandersypen, L. M. K. et al. Interfacing spin qubits in quantum dots and donors—hot, dense, and coherent.npj Quantum Inf.3, 34 (2017).

Degenhardt, C., Geck, L., Kruth, A., Vliex, P. & van Waasen, S. CMOS based scalable cryogenic control electronics for qubits. In2017 IEEE International Conference on Rebooting Computing (ICRC), 1–4 (IEEE, 2017).

Pauka, S. J. et al. A cryogenic CMOS chip for generating control signals for multiple qubits.Nat. Electron.4, 64–70 (2021).

Hornibrook, J. M. et al. Cryogenic control architecture for large-scale quantum computing.Phys. Rev. Appl.3, 024010 (2015).

Patra, B. et al. Cryo-CMOS Circuits and Systems for Quantum Computing Applications.IEEE J. Solid-State Circuits53, 309–321 (2018).

Geck, L., Kruth, A., Bluhm, H., van Waasen, S. & Heinen, S. Control electronics for semiconductor spin qubits.Quantum Sci. Technol.5, 015004 (2019).

Dial, O. E. et al. Charge noise spectroscopy using coherent exchange oscillations in a singlet-triplet qubit.Phys. Rev. Lett.110, 146804 (2013).

Thiney, V. et al. In-flight detection of few electrons using a singlet-triplet spin qubit.Phys. Rev. Res.4, 043116 (2022).

Beverland, M. E. et al. Assessing requirements to scale to practical quantum advantage (2022). Preprint atarxiv.org/abs/2211.07629(2022).

Reilly, D. J. Engineering the quantum-classical interface of solid-state qubits.npj Quantum Inf.1, 15011 (2015).

Kane, B. E. A silicon-based nuclear spin quantum computer.Nature393, 133–137 (1998).

Veldhorst, M. et al. An addressable quantum dot qubit with fault-tolerant control-fidelity.Nat. Nanotechnol.9, 981–985 (2014).

Yang, C. H. et al. Operation of a silicon quantum processor unit cell above one kelvin.Nature580, 350–354 (2020).

Petit, L. et al. Universal quantum logic in hot silicon qubits.Nature580, 355–359 (2020).

Huang, J. Y. et al. High-fidelity spin qubit operation and algorithmic initialization above 1 K.Nature627, 772–777 (2024).

Xue, X. et al. CMOS-based cryogenic control of silicon quantum circuits.Nature593, 205–210 (2021).

Itoh, K. M., & Watanabe, H. Isotope engineering of silicon and diamond for quantum computing and sensing applications.MRS Commun.4, 143–157 (2014).

Angus, S. J., Ferguson, A. J., Dzurak, A. S. & Clark, R. G. A silicon radio-frequency single electron transistor.Appl. Phys. Lett.92, 112103 (2008).

Reilly, D. J., Marcus, C. M., Hanson, M. P. & Gossard, A. C. Fast single-charge sensing with a rf quantum point contact.Appl. Phys. Lett.91, 162101 (2007).

Petta, J. R. et al. Coherent manipulation of coupled electron spins in semiconductor quantum dots.Science309, 2180–2184 (2005).

Martins, F. et al. Noise suppression using symmetric exchange gates in spin qubits.Phys. Rev. Lett.116, 116801 (2016).

ArticleMathSciNetGoogle Scholar

Ono, K., Austing, D. G., Tokura, Y. & Tarucha, S. Current rectification by Pauli exclusion in a weakly coupled double quantum dot system.Science297, 1313–1317 (2002).

Johnson, A. C., Petta, J. R., Marcus, C. M., Hanson, M. P. & Gossard, A. C. Singlet-triplet spin blockade and charge sensing in a few-electron double quantum dot.Phys. Rev. B72, 165308 (2005).

Knill, E. et al. Randomized benchmarking of quantum gates.Phys. Rev. A77, 012307 (2008).

Magesan, E., Gambetta, J. M. & Emerson, J. Characterizing quantum gates via randomized benchmarking.Phys. Rev. A85, 042311 (2012).

Vahapoglu, E. et al. Coherent control of electron spin qubits in silicon using a global field.npj Quantum Inf.8, 126 (2022).

Laucht, A. et al. Electrically controlling single-spin qubits in a continuous microwave field.Sci. Adv.1, e1500022 (2015).

Laucht, A. et al. A dressed spin qubit in silicon.Nat. Nanotechnol.12, 61–66 (2017).

Watson, T. F. et al. A programmable two-qubit quantum processor in silicon.Nature555, 633–637 (2018).

Huang, W. et al. Fidelity benchmarks for two-qubit gates in silicon.Nature569, 532–536 (2019).

Nielsen, E. et al. Gate set tomography.Quantum5, 557 (2021).

Dubois, J. et al. Minimal-excitation states for electron quantum optics using levitons.Nature502, 659–663 (2013).

Weinstein, A. J. et al. Universal logic with encoded spin qubits in silicon.Nature615, 817–822 (2023).

Hornibrook, J. M. et al. Frequency multiplexing for readout of spin qubits.Appl. Phys. Lett.104, 103108 (2014).

Colless, J. I. & Reilly, D. J. Modular cryogenic interconnects for multi-qubit devices.Rev. Sci. Instrum.85, 114706 (2014).

Waddy, S.Quantum Readout and Control: Scalable Hardware and Techniques for Solid State Quantum Systems. PhD thesis. The Univ. Sydney (2023).

Bartee, S. Data used in “Spin-qubit control with a milli-kelvin CMOS chip”.Zenodohttps://doi.org/10.5281/zenodo.15080655(2025).

Cywiński, Ł., Lutchyn, R. M., Nave, C. P. & Sarma, S. D. How to enhance dephasing time in superconducting qubits.Phys. Rev. B77, 174509 (2008).

Álvarez, G. A., & Suter, D. Measuring the spectrum of colored noise by dynamical decoupling.Phys. Rev. Lett.107, 230501 (2011).

Medford, J. et al. Scaling of dynamical decoupling for spin qubits.Phys. Rev. Lett.108, 086802 (2012).

Philips, S. G. J. et al. Universal control of a six-qubit quantum processor in silicon.Nature609, 919–924 (2022).

Undseth, B. et al. Hotter is easier: unexpected temperature dependence of spin qubit frequencies.Phys. Rev. X13, 041015 (2023).

We thank Y. Yang and R. Kalra for technical contributions and discussions. This research was supported by Microsoft Corporation (CPD 1-4) and by the Australian Research Council Centre of Excellence for Engineered Quantum Systems (EQUS, CE170100009). We acknowledge support from the Australian Research Council (FL190100167 and CE170100012), the US Army Research Office (W911NF-23-10092), the US Air Force Office of Scientific Research (FA2386-22-1-4070) and the NSW Node of the Australian National Fabrication Facility as well as the Research and Prototype Foundry Core Research Facility at the University of Sydney, also part of the Australian National Fabrication Facility. The views and conclusions in this document are those of the authors and should not be interpreted as representing the official policies, either expressed or implied, of Microsoft Corporation, the Army Research Office, the US Air Force or the US government. The US government is authorized to reproduce and distribute reprints for government purposes, notwithstanding any copyright notation herein. R.Y.S. and S.S. acknowledge support from the Sydney Quantum Academy.

Open access funding provided by the University of Sydney.

Present address: Diraq, Sydney, New South Wales, Australia

ARC Centre of Excellence for Engineered Quantum Systems, School of Physics, The University of Sydney, Sydney, New South Wales, Australia

Samuel K. Bartee, Kun Zuo, Kushal Das, Sebastian J. Pauka & David J. Reilly

Diraq, Sydney, New South Wales, Australia

Will Gilbert, Tuomo Tanttu, Chih Hwan Yang, Nard Dumoulin Stuyck, Wee Han Lim, Santiago Serrano, Christopher C. Escott, Fay E. Hudson, Arne Laucht & Andrew S. Dzurak

School of Electrical Engineering and Telecommunications, University of New South Wales, Sydney, New South Wales, Australia

Will Gilbert, Tuomo Tanttu, Chih Hwan Yang, Nard Dumoulin Stuyck, Rocky Y. Su, Wee Han Lim, Santiago Serrano, Christopher C. Escott, Fay E. Hudson, Arne Laucht & Andrew S. Dzurak

Emergence Quantum, Sydney, New South Wales, Australia

Kushal Das, Sebastian J. Pauka & David J. Reilly

School of Fundamental Science and Technology, Keio University, Yokohama, Japan

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

S.K.B., W.G., K.Z., T.T., C.H.Y., N.D.S., R.Y.S., S.J.P. and D.J.R. designed the experiments. K.D., S.J.P. and D.J.R. designed and characterized the CMOS chip. S.K.B. and K.Z. performed experiments with input from W.G., K.D., T.T., S.S., C.C.E. and A.L., under the supervision of A.S.D. and D.J.R.; W.H.L. and F.E.H. fabricated the device on enriched28Si wafers supplied by K.M.I.; S.K.B., K.Z. and D.J.R. wrote the paper with input from all co-authors.

Correspondence toDavid J. Reilly.

A.S.D. is the CEO and director of Diraq. S.K.B., W.G., T.T., C.H.Y., N.D.S., W.H.L., C.C.E., F.E.H. and A.L. declare equity interest in Diraq. D.J.R. is the CEO and Director of Emergence Quantum. K.D. and S.J.P. declare equity interests in Emergence Quantum. Other authors declare no competing interests.

Naturethanks the anonymous reviewer(s) for their contribution to the peer review of this work.Peer reviewer reportsare available.

Publisher’s noteSpringer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

For further hardware information, seeMethods.

a, the barrier to the SET reservoir is lowered for a set number of electrons to occupy the dots (i). This barrier is raised to isolate the dots (ii), and the electron states are initialized via tuning of the plunger gates P1 and P2 (iii).b, Single-qubit pulsing scheme for ESR measurements. A diabatic detuning ramp (1 → 2) from the (0,4) to (1,3) state initialises the double quantum dot into a singlet state. The qubit is then pulsed using theJ-gate to the control point (2 → 3) at which point a MW pulse is applied that rotates the target qubit in resonance with the MW frequency. The qubit is pulsed back to ♦(3 → 4), then pulsed into the Pauli Spin Blockade regime (4 → 5) for readout. Detuningϵpulses use gatesP1 andP2 and are always RT-operated.J-gate pulses are generated either at RT or by cryo-CMOS. The approximate operation points of single qubit (X, Z), two qubit (CZ), and readout points (M1Q,M2Q) are indicated by a triangle (▴), square (▪), star (⋆) and diamond (♦) respectively.c(inset of b), basicJ-gate pulsing scheme for ESR measurements when using cryo-CMOS. Here, the gateJis in a “charge-locked” state held atVHOLD. An external trigger is used to pulse between the twoJ-gate levels, whose separation is determined byΔVJ,PULSE.d, DCZ sequence focusing onJ-gate pulses. The qubits are initialised into aT−state, beforeXandX/2 gates are performed on the control and target qubits respectively. A spin-echo sequence is inserted between the exchange gates, which precedes the final projection for measurement of the two-qubit spin state.

a, By forward biasing an ESD protection diode on-chip, the power draw can be programmed to mimic that of any oscillator trim or division value upon oscillator activation. Artificial power draw also allows for higher resolution investigation versus other parameters, as well as extrapolation beyond what the maximal cryo-CMOS power draw possible with it’s feature set.b, Schematic of the cryo-CMOS ring oscillator. 81 in-series inverters are connected to a three-bit multiplexer, controlling the tap-off point. The frequency is tunable by programmable inputs into the oscillator trim. This frequency is further divided by the oscillator divider, which is eight-bit tunable. The ultimate output frequency is then passed on to the FSM.c, Hahn echo coherence time and Ramsey coherence time of Q1, and measured electron temperature as a function of mixing chamber temperature. Base effective electron temperature is approximately 850 mK, which only starts to increase once the mixing chamber exceeds 700 mK. Electron temperature and mixing chamber temperature equalize at 1 K. Each data point consists of 100 shots with 4 repeats for a total of 400 single shots. Error bars represent the 95 % confidence level.d(1-3)-(0,4) charge occupation probability. Solid lines are Fermi distribution fits, allowing for extraction of effective electron temperature. Fitting to distributions with high mixing chamber temperature (at which point becoming equal to that of the sample) allows for determination of the lever arm.

a,\({T}_{2}^{Hahn}\)coherence of both qubits as a function of divided frequency from the cryo-CMOS oscillator. The output clock frequency (here set to 30 MHz) passed to the fast state machine is divided by an integer between 1 and 255, see Extended Data Fig.2for schematic details. The inset indicates whether RT or cryo-CMOS control is used; all RT pulsing schemes are replicable by cryo-CMOS. Activating the oscillator immediately causes a drop in coherence due to the extra thermal dissipation, and lowering the divider value increases this dissipation slightly.b, Qubit coherence while directly changing oscillator frequency. The divider is set to a constant value of 255, and increasing oscillator frequency leads to an increase in power draw of the cryo-CMOS chip, reflected in the mixing chamber temperature.c, Similar to Fig.2(c),\({T}_{2}^{Hahn}\)coherence is also explored as a function of various parameters under cryo-CMOS control conditions.d, e,\({T}_{2}^{* }\)of both qubits while changing the divided and oscillator frequencies, similar to a and b. For all data points, each data point consists of 100 shots with 4 repeats for a total of 400 single shots.f, The noise power spectral density (PSD) is explored as a function of oscillator frequency. Noise spectroscopy, based on the Carr-Purcell-Meiboom-Gill (CPMG) protocol42,43,44, uses a single qubit as a noise probe. A slight increase in the overall white noise level is observed over the detectable frequency range when the oscillator is at its maximum frequency. This increase in white noise is consistent with parasitic heating leading to an increase in temperature and thermal noise contribution to the cumulative noise background. We note that we do not observe additional electrical frequency spurs or artifacts beyond the thermal noise contribution. At higher frequencies, we observe an increase in PSD, likely due to high-power driving or mis-calibration of microwave pulses19,45,46. Error bars represent the 95 % confidence level.

a, ESR spectrum as a function ofVJ. Exchange starts to open at aroundVJ= 1.37 V. Measurements are done with room temperatureJpulses; RT or cryo-CMOS control is indicated in the lower right hand side of each figure.b, c, CZ oscillations as a function of exchange timetexchangeandVJ. MultipleJlevels, not replicable by cryo-CMOS control are used here for optimal performance. The initial state of the target control qubits is indicated in the bottom right hand corner. In (b), theJgate drives two-qubit exchange at a sensitivity of 25dec\V.d, CZ oscillations using cryo-CMOS control, showing the difference in visibility compared tob-caveraged over 100 shots. Traces are offset for clarity. For all further figures, pulses are two-level and if generated at RT, are replicable by cryo-CMOS.e, similar to Extended Data Fig.4(a), the two qubit CZ coherence time\({T}_{2}^{exchange}\)is explored as a function of divided frequency. The oscillator is set to 30 MHz for all divisions, and pulses are generated from room temperature.f, CZ coherence time\({T}_{2}^{exchange}\)is now explored as a function of oscillator frequency. The divider is set to 255 for all measurements.g, h, Measurements ine, fare replicated here, however now using a DCZ pulsing protocol, allowing for longer coherence times.i, DCZ coherence is explored using cryo-CMOS control under various parameters, similar to Fig.4(d). For all individual data points, each data point consists of 100 shots with 4 repeats for a total of 400 single shots. Error bars represent the 95 % confidence level.

Open AccessThis article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visithttp://creativecommons.org/licenses/by/4.0/.

Bartee, S.K., Gilbert, W., Zuo, K.et al.Spin-qubit control with a milli-kelvin CMOS chip.Nature(2025). https://doi.org/10.1038/s41586-025-09157-x

DOI:https://doi.org/10.1038/s41586-025-09157-x

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

Coherent bunching of anyons and dissociation in an interference experiment

Aharonov–Bohm interference of fractional quasiparticles in the quantum Hall effect generally reveals their elementary charge (e*)1,2,3,4,5,6,7,8,9,10,11,12,13,14,15. Recently, our interferometry experiments with several ‘particle states’ reported flux periods of ΔΦ= (e/e*)Φ0(withΦ0the flux quantum) at moderate temperatures16. Here we report interference measurements of ‘particle–hole conjugated’ states at filling factorsν= 2/3, 3/5 and 4/7, which revealed unexpected flux periodicities of ΔΦ=ν−1Φ0. The measured shot-noise Fano factor (F) of the partitioned quasiparticles in each of the quantum point contacts of the interferometer wasF=ν(ref.17) rather than that of the elementary chargeF=e*/e(refs.18,19). These observations indicate that the interference of bunched (clustered) elementary quasiparticles occurred for coherent pairs, triples and quadruplets, respectively. A small metallic gate (top gate), deposited in the centre of the interferometer bulk, formed an antidot (or a dot) when charged, thus introducing local quasiparticles at the perimeter of the (anti)dot. Surprisingly, such charging led to a dissociation of the ‘bunched quasiparticles’ and, thus, recovered the conventional flux periodicity set by the elementary charge of the quasiparticles. However, the shot-noise Fano factor (of each quantum point contact) consistently remained atF=ν, possibly due to the neutral modes accompanying the conjugated states. The two observations—bunching and debunching (or dissociation)—were not expected by current theories. Similar effects may arise in Jain’s ‘particle states’ (at lower temperatures) and at even denominator fractional quantum Hall states20.

This is a preview of subscription content,access via your institution

Access Nature and 54 other Nature Portfolio journals

Get Nature+, our best-value online-access subscription

Receive 51 print issues and online access

Prices may be subject to local taxes which are calculated during checkout

The data that support the plots within this paper and other findings of this study are publicly available athttps://doi.org/10.5281/zenodo.15395283(ref.39).Source dataare provided with this paper.

Zhang, Y. M. et al. Distinct signatures for Coulomb blockade and Aharonov-Bohm interference in electronic Fabry-Perot interferometers.Phys. Rev. B79, 241304 (2009).

Ofek, N. et al. Role of interactions in an electronic Fabry–Perot interferometer operating in the quantum Hall effect regime.Proc. Natl Acad. Sci. USA107, 5276–5281 (2010).

ArticleCASPubMedPubMed CentralGoogle Scholar

Ronen, Y. et al. Aharonov-Bohm effect in graphene-based Fabry-Pérot quantum Hall interferometers.Nat. Nanotechnol.16, 563–569 (2021).

Nakamura, J., Liang, S., Gardner, G. C. & Manfra, M. J. Direct observation of anyonic braiding statistics.Nat. Phys.16, 931–936 (2020).

Kim, J. et al. Aharonov–Bohm interference and statistical phase-jump evolution in fractional quantum Hall states in bilayer graphene.Nat. Nanotechnol.19, 1619–1626 (2024).

Nakamura, J. et al. Aharonov–Bohm interference of fractional quantum Hall edge modes.Nat. Phys.15, 563–569 (2019).

Déprez, C. et al. A tunable Fabry-Pérot quantum Hall interferometer in graphene.Nat. Nanotechnol.16, 555–562 (2021).

ArticlePubMedPubMed CentralGoogle Scholar

Samuelson, N. L. et al. Anyonic statistics and slow quasiparticle dynamics in a graphene fractional quantum Hall interferometer. Preprint athttps://arxiv.org/abs/2403.19628(2024).

Werkmeister, T. et al. Anyon braiding and telegraph noise in a graphene interferometer.Science388, 730–735 (2024).

Chamon, C. D. C., Freed, D. E., Kivelson, S. A., Sondhi, S. L. & Wen, X. G. Two point-contact interferometer for quantum Hall systems.Phys. Rev. B55, 2331–2343 (1997).

Willett, R. L. et al. Interference measurements of non-abeliane/4 & abeliane/2 quasiparticle braiding.Phys. Rev.13, 011028 (2023).

Yang, W. et al. Evidence for correlated electron pairs and triplets in quantum Hall interferometers.Nat. Commun.15, 10064 (2024).

ArticleCASPubMedPubMed CentralGoogle Scholar

Ji, Y. et al. An electronic Mach-Zehnder interferometer.Nature422, 415–418 (2003).

Neder, I., Heiblum, M., Levinson, Y., Mahalu, D. & Umansky, V. Unexpected behavior in a two-path electron interferometer.Phys. Rev. Lett.96, 016804 (2006).

Kundu, H. K., Biswas, S., Ofek, N., Umansky, V. & Heiblum, M. Anyonic interference and braiding phase in a Mach-Zehnder interferometer.Nat. Phys.19, 515–521 (2023).

Ghosh, B. et al. Anyonic braiding in a chiral Mach–Zehnder interferometer.Nat. Phys.https://doi.org/10.1038/s41567-025-02960-3(in the press).

Bhattacharyya, R., Banerjee, M., Heiblum, M., Mahalu, D. & Umansky, V. Melting of interference in the fractional quantum Hall effect: appearance of neutral modes.Phys. Rev. Lett.122, 246801 (2019).

de-Picciotto, R. et al. Direct observation of a fractional charge.Phys. B: Condens. Matter249–251, 395–400 (1998).

Reznikov, M., de Picciotto, R., Griffiths, T. G., Heiblum, M. & Umansky, V. Observation of quasiparticles with one-fifth of an electron’s charge.Nature399, 238–241 (1999).

Chung, Y. C., Heiblum, M. & Umansky, V. Scattering of bunched fractionally charged quasiparticles.Phys. Rev. Lett.91, 216804 (2003).

Schuster, R. et al. Phase measurement in a quantum dot via a double-slit interference experiment.Nature385, 417–420 (1997).

Neder, I. et al. Interference between two indistinguishable electrons from independent sources.Nature448, 333–337 (2007).

dePicciotto, R. et al. Direct observation of a fractional charge.Nature389, 162–164 (1997).

Saminadayar, L., Glattli, D. C., Jin, Y. & Etienne, B. Observation of thee/3 fractionally charged Laughlin quasiparticle.Phys. Rev. Lett.79, 2526–2529 (1997).

Simmons, J. A. et al. Resistance fluctuations in the integral-quantum-Hall-effect and fractional-quantum-Hall-effect regimes.Phys. Rev. B44, 12933–12944 (1991).

Lee, J.-Y. M. et al. Partitioning of diluted anyons reveals their braiding statistics.Nature617, 277–281 (2023).

Glidic, P. et al. Cross-correlation investigation of anyon statistics in the ν = 1/3 and 2/5 fractional quantum Hall states.Phys. Rev.13, 011030 (2023).

Ruelle, M. et al. Comparing fractional quantum Hall Laughlin and Jain topological orders with the anyon collider.Phys. Rev.13, 011031 (2023).

Bartolomei, H. et al. Fractional statistics in anyon collisions.Science368, 173–177 (2020).

ArticleMathSciNetCASPubMedGoogle Scholar

Biswas, S. et al. Shot noise does not always provide the quasiparticle charge.Nat. Phys.18, 1476–1481 (2022).

Deviatov, E. V., Egorov, S. V., Biasiol, G. & Sorba, L. Quantum Hall Mach-Zehnder interferometer at fractional filling factors.Europhys. Lett.100, 67009 (2012).

Batra, N., Wei, Z., Vishveshwara, S. & Feldman, D. E. Anyonic Mach-Zehnder interferometer on a single edge of a two-dimensional electron gas.Phys. Rev. B108, L241302 (2023).

Bid, A., Ofek, N., Heiblum, M., Umansky, V. & Mahalu, D. Shot noise and charge at the 2/3 composite fractional quantum Hall state.Phys. Rev. Lett.103, 236802 (2009).

Biswas, S., Kundu, H. K., Umansky, V. & Heiblum, M. Electron pairing of interfering interface-based edge modes.Phys. Rev. Lett.131, 096302 (2023).

Choi, H. K. et al. Robust electron pairing in the integer quantum Hall effect regime.Nat. Commun.6, 7435 (2015).

Frigeri, G. A., Scherer, D. D. & Rosenow, B. Sub-periods and apparent pairing in integer quantum Hall interferometers.Europhys. Lett.126, 67007 (2019).

Kim, J. et al. Aharonov-Bohm interference in even-denominator fractional quantum Hall states. Preprint athttps://arxiv.org/abs/2412.19886(2024).

Moore, G. & Read, N. Nonabelions in the fractional quantum Hall effect.Nucl. Phys. B360, 362–396 (1991).

ArticleMathSciNetGoogle Scholar

Ghosh, B. Coherent bunching of anyons and their dissociation in interference experiments.Zenodohttps://doi.org/10.5281/zenodo.15395283(2025).

M.H. thanks Y. Ronen for fruitful discussions. D.F.M. acknowledges many illuminating conversations on quantum Hall interferometry with Y. Ronen. B.G. thanks A. K. Paul for his helpful comment that improved the device. M.L. thanks the Ariane de Rothschild Women Doctoral Program for their support. D.F.M. acknowledges the support of the Israel Science Foundation (Grant No. 2572/21) and the Deutsche Forschungsgemeinschaft (DFG) within the CRC network TR 183 (Project Grant No. 277101999). M.H. acknowledges the support of the European Research Council under the European Union’s Horizon 2020 research and innovation programme (Grant Agreement No. 833078) and the support of the Israel Science Foundation (Grant No. 1510/22).

Braun Center for Sub-Micron Research, Department of Condensed Matter Physics, Weizmann Institute of Science, Rehovot, Israel

Bikash Ghosh, Maria Labendik, Vladimir Umansky & Moty Heiblum

Department of Condensed Matter Physics, Weizmann Institute of Science, Rehovot, Israel

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

B.G. fabricated the devices. B.G. and M.L. performed the measurements and analysed the data with the input from M.H. M.H. supervised the design, execution and data analysis in the experiment. D.F.M. worked on the theoretical aspects and data analysis. V.U. grew the GaAs heterostructures. All authors contributed to the writing of the manuscript.

The authors declare no competing interests.

Naturethanks the anonymous reviewers for their contribution to the peer review of this work.Peer reviewer reportsare available.

Publisher’s noteSpringer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Notes, methods, Figs. 1–11 and references.

Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.

Ghosh, B., Labendik, M., Umansky, V.et al.Coherent bunching of anyons and dissociation in an interference experiment.Nature(2025). https://doi.org/10.1038/s41586-025-09143-3

DOI:https://doi.org/10.1038/s41586-025-09143-3

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

These are the top companies and countries for industry research

Search author on:PubMedGoogle Scholar

More than 1 in 20 of all research moves within the United States for the period 2020–24 were from academia to companies, according to data from Sydney-based analytics company League of Scholars. Moves were calculated based on verified Google Scholar accounts affiliated with company e-mail domains, which means moves in countries where Google Scholar use is relatively low among corporate researchers, such as China, might be underrepresented.

US companies lead as top destinations for academics moving to corporate research roles, based on 68,210 verified global cross-sector movements for 2020–24. Nearly 4% of these transitions were to Amazon and Google, with just over 2.5% specifically to Amazon. Samsung in South Korea and Huawei in China are the only non-US companies in this list.

Of the corporate-research output in the Nature Index, these are the five fastest-rising natural-sciences fields — measured by absolute change in article Count from 2020 to 2024. All fields sit within the top 10 overall.

The top natural-sciences* fields in corporate-authored research in the Nature Index, as measured by article count, is shown for 2020–24. Physical chemistry is the fastest-rising field overall, increasing its Count by 171.29% for the period.

Gain free access to this article, as well as selected content from this journal and more onnature.com

doi: https://doi.org/10.1038/d41586-025-01929-9

This article is part ofNature Index 2025 Science Inc., an editorially independent supplement. Advertisers have no influence over the content. For more information about Nature Index,see the homepage.

Can industry fill the gap left by US research funding cuts?

What incentives do companies need to publish research?

Here’s why China’s science and innovation model is thriving

Three science-led start-ups making their mark

Partner content: How AI is combing through real-world data to transform healthcare

Partner content: Swarms of AI robots could extend human capabilities

Partner content: Quantum leaps promise to reshape the internet

Partner content: How exploring the unknown keeps the innovations coming

Partner content: Giving plastics a whole lease of life

Partner content: Advancing precision nutrition for better health

Can industry fill the gap left by US research funding cuts?

Three science-led start-ups making their mark

What incentives do companies need to publish research?

NIH-funded science must now be free to read instantly: what you should know

Can industry fill the gap left by US research funding cuts?

AI, peer review and the human activity of science

Can industry fill the gap left by US research funding cuts?

Three science-led start-ups making their mark

Founded in 2018, Westlake University is a new type of non-profit research-oriented university in Hangzhou, China, supported by public a…

The University of Toronto now recruiting for the Eric and Wendy Schmidt AI in Science Postdoctoral Fellowship. Valued at $85,000 CDN per year.

The Integrated Science Lab (IceLab) offers three postdoctoral scholarships that will be affiliated with one of six possible multidisciplinary projects

Aim to a global leader in materials science research society, and create a diverse, innovative, and collaborative research environment

Department of Frontier Materials, Suzhou Laboratory

The Max Planck Institute of Molecular Cell Biology and Genetics (MPI-CBG) in Dresden, Germany, is seeking an outstanding individual to fill the po…

Max Planck Institute of Molecular Cell Biology and Genetics

Is AI powering Big Brother? Surveillance research is on the rise

Search author on:PubMedGoogle Scholar

You have full access to this article via your institution.

Computer-vision research powers surveillance technology

Computer-vision research powers surveillance technology

AI research on computer vision is almost always used in ways that enable surveillance, according to a newNaturestudy.

Computer vision is what allows driverless cars to navigate, automatic tagging of people's faces on social media, and the detection of cancerous cells, but many people have argued it could also power mass surveillance. The new study trawled through over 40,000 documents to discover that the vast majority of research in computer vision powers surveillance.

Such technology could be used to secure borders and catch criminals, but researchers are concerned that this could impact democratic freedoms by allowing easy identification and tracking of people. The study could help empower those who want to protect people from being surveilled, by highlighting how this research is being conducted and used.

doi: https://doi.org/10.1038/d41586-025-01936-w

Meet the MIT engineer who invented an AI-powered way to restore art

How AI is revealing the language of the birds

Watch this bird-inspired drone leap into the air

ChatGPT has a language problem — but science can fix it

What ChatGPT is and what it's not: a three minute guide

Can AI build a virtual cell? Scientists race to model life’s smallest unit

AI, peer review and the human activity of science

Don’t sleepwalk from computer-vision research into surveillance

Don’t sleepwalk from computer-vision research into surveillance

Wake up call for AI: computer-vision research increasingly used for surveillance

Brain-reading devices raise ethical dilemmas — researchers propose protections

Don’t sleepwalk from computer-vision research into surveillance

Computer-vision research is hiding its role in creating ‘Big Brother’ technologies

Does using ChatGPT change your brain activity? Study sparks debate

Founded in 2018, Westlake University is a new type of non-profit research-oriented university in Hangzhou, China, supported by public a…

The University of Toronto now recruiting for the Eric and Wendy Schmidt AI in Science Postdoctoral Fellowship. Valued at $85,000 CDN per year.

The Integrated Science Lab (IceLab) offers three postdoctoral scholarships that will be affiliated with one of six possible multidisciplinary projects

Aim to a global leader in materials science research society, and create a diverse, innovative, and collaborative research environment

Department of Frontier Materials, Suzhou Laboratory

The Max Planck Institute of Molecular Cell Biology and Genetics (MPI-CBG) in Dresden, Germany, is seeking an outstanding individual to fill the po…

Max Planck Institute of Molecular Cell Biology and Genetics

Hormone helps mice live longer and avoid obesity

Raising the level of a hormone in fat tissue (pictured, fat cells; artificially coloured) could promote healthy ageing. Credit: Steve Gschmeissner/SPL

A genetic modification that causes mice to overexpress a particular hormone prevents them from gaining weight and allows them to live longer than unmodified mice1.

Access Nature and 54 other Nature Portfolio journals

Get Nature+, our best-value online-access subscription

Receive 51 print issues and online access

Prices may be subject to local taxes which are calculated during checkout

doi: https://doi.org/10.1038/d41586-025-01896-1

Gliniak, C. M.et al. Cell Metab. https://doi.org/10.1016/j.cmet.2025.05.011 (2025).

Can a pill replace exercise? Swigging this molecule gives mice benefits of working out

Is there a link between gut microbes and ageing?

How your brain controls ageing — and why zombie cells could be key

Founded in 2018, Westlake University is a new type of non-profit research-oriented university in Hangzhou, China, supported by public a…

The University of Toronto now recruiting for the Eric and Wendy Schmidt AI in Science Postdoctoral Fellowship. Valued at $85,000 CDN per year.

The Integrated Science Lab (IceLab) offers three postdoctoral scholarships that will be affiliated with one of six possible multidisciplinary projects

Aim to a global leader in materials science research society, and create a diverse, innovative, and collaborative research environment

Department of Frontier Materials, Suzhou Laboratory

The Max Planck Institute of Molecular Cell Biology and Genetics (MPI-CBG) in Dresden, Germany, is seeking an outstanding individual to fill the po…

Max Planck Institute of Molecular Cell Biology and Genetics

Is AI watching you? The hidden links between research and surveillance

Search author on:PubMedGoogle Scholar

Search author on:PubMedGoogle Scholar

You have full access to this article via your institution.

Download theNature Podcast25 June 2025

We’d like to learn more about our listeners, please help us out byfilling in this short survey.

A significant amount of research in the AI field of computer vision is being used to analyse humans in ways that support the development of surveillance technologies, according to new research. By analysing the contents of thousands of research papers, the team behind the work showed that 90% of studies, and 86% of patents resulting from them, involved data relating to imaging humans. Although there are many positive applications for this technology, such as in medical diagnostics, this work shows evidence of a pipeline from computer-vision research to surveillance.

Research article:Kalluri et al.

News and Views:Computer-vision research is hiding its role in creating ‘Big Brother’ technologies

Video:Is AI powering Big Brother? Surveillance research is on the rise

News:Wake up call for AI: computer-vision research increasingly used for surveillance

Editorial:Don’t sleepwalk from computer-vision research into surveillance

A minuscule robot that can manipulate liquid droplets, and the sensors that can identify hydrothermal explosions at Yellowstone National Park.

Research Highlight:This tiny robot moves mini-droplets with ease

Research Highlight:Sensors pinpoint the exact time of a Yellowstone explosion

This week, the Vera C. Rubin Observatory in Chile has unveiled the first images from its 3,200 megapixel digital camera. We discuss the images, and how the facility will help researchers learn more about the Universe.

News:First images from world’s largest digital camera leave astronomers in awe

How scientists created hexanitrogen, a new molecule made of six nitrogen atoms, and the why researchers are excited about the first confirmed skull of an extinct Denisovan.

Chemistry World:Most energetic molecule ever made is stable – in liquid nitrogen

Nature:First ever skull from ‘Denisovan’ reveals what ancient people looked like

Subscribe to Nature Briefing, an unmissable daily round-up of science news, opinion and analysis free in your inbox every weekday.

Never miss an episode. Subscribe to the Nature Podcast onApple Podcasts,Spotify,YouTube Musicor your favourite podcast app.An RSS feed for the Nature Podcastis available too.

doi: https://doi.org/10.1038/d41586-025-02004-z

Is AI powering Big Brother? Surveillance research is on the rise

Flight simulator for moths reveals they navigate by starlight

Mum’s a Neanderthal, Dad’s a Denisovan: First discovery of an ancient-human hybrid

Biggest Denisovan fossil yet spills ancient human’s secrets

Can industry fill the gap left by US research funding cuts?

These are the top companies and countries for industry research

Three science-led start-ups making their mark

Morphodynamics of human early brain organoid development

An open AI model could help medical experts to interpret chest X-rays

How the brain separates real images from those it imagines

Are these the world’s oldest rocks? Fresh data confirm contested claim

Sensors pinpoint the exact time of a Yellowstone explosion

One of the world’s richest lithium deposits began inside a mega-volcano

Founded in 2018, Westlake University is a new type of non-profit research-oriented university in Hangzhou, China, supported by public a…

The University of Toronto now recruiting for the Eric and Wendy Schmidt AI in Science Postdoctoral Fellowship. Valued at $85,000 CDN per year.

The Integrated Science Lab (IceLab) offers three postdoctoral scholarships that will be affiliated with one of six possible multidisciplinary projects

Aim to a global leader in materials science research society, and create a diverse, innovative, and collaborative research environment

Department of Frontier Materials, Suzhou Laboratory

The Max Planck Institute of Molecular Cell Biology and Genetics (MPI-CBG) in Dresden, Germany, is seeking an outstanding individual to fill the po…

Max Planck Institute of Molecular Cell Biology and Genetics

Pioneering but overlooked 1938 fusion experiment is recreated at last

Physicist Arthur Ruhlig was the first to fuse deuterium and tritium (top and bottom left; artist’s impression), creating helium (top right) and releasing a large amount of energy. Credit: Science Photo Library

Physicists have recreated the first experiment to fuse two specific versions of hydrogen – a result that was overlooked for 85 years1.

Access Nature and 54 other Nature Portfolio journals

Get Nature+, our best-value online-access subscription

Receive 51 print issues and online access

Prices may be subject to local taxes which are calculated during checkout

doi: https://doi.org/10.1038/d41586-025-01991-3

Tornow, W.et al.Phys. Rev. C.111, 064618 (2025).

Did the US wipe out Iran's nuclear programme? What researchers know

Richard Garwin obituary: designer of the hydrogen bomb who advocated for arms control

An asymmetric fission island driven by shell effects in light fragments

Founded in 2018, Westlake University is a new type of non-profit research-oriented university in Hangzhou, China, supported by public a…

The University of Toronto now recruiting for the Eric and Wendy Schmidt AI in Science Postdoctoral Fellowship. Valued at $85,000 CDN per year.

The Integrated Science Lab (IceLab) offers three postdoctoral scholarships that will be affiliated with one of six possible multidisciplinary projects

Aim to a global leader in materials science research society, and create a diverse, innovative, and collaborative research environment

Department of Frontier Materials, Suzhou Laboratory

The Max Planck Institute of Molecular Cell Biology and Genetics (MPI-CBG) in Dresden, Germany, is seeking an outstanding individual to fill the po…

Max Planck Institute of Molecular Cell Biology and Genetics

Can industry fill the gap left by US research funding cuts?

Search author on:PubMedGoogle Scholar

Nicholas Flanders and Etosha Cave, co-founders of carbon-conversion start-up Twelve in Berkeley, California. Credit: Mark Leongfor Nature Index

As the United Statescuts billions of dollars in federal funds for research, the role of the private sector in driving and supporting science and innovation has been put under the spotlight. Could US pharmaceutical giants support drug- or vaccine-development projects that have lost their funding? Will tech and engineering firms use their clout to advance defundedresearch in green energy? And if companies do fill the funding gap, how might the priorities of research projects shift to meet business goals and pressures?

Gain free access to this article, as well as selected content from this journal and more onnature.com

doi: https://doi.org/10.1038/d41586-025-01925-z

This article is part ofNature Index 2025 Science Inc., an editorially independent supplement. Advertisers have no influence over the content. For more information about Nature Index,see the homepage.

What incentives do companies need to publish research?

Here’s why China’s science and innovation model is thriving

Three science-led start-ups making their mark

These are the top companies and countries for industry research

Partner content: How AI is combing through real-world data to transform healthcare

Partner content: Swarms of AI robots could extend human capabilities

Partner content: Quantum leaps promise to reshape the internet

Partner content: How exploring the unknown keeps the innovations coming

Partner content: Giving plastics a whole lease of life

Partner content: Advancing precision nutrition for better health

These are the top companies and countries for industry research

Three science-led start-ups making their mark

What incentives do companies need to publish research?

NIH-funded science must now be free to read instantly: what you should know

These are the top companies and countries for industry research

AI, peer review and the human activity of science

These are the top companies and countries for industry research

Three science-led start-ups making their mark

Founded in 2018, Westlake University is a new type of non-profit research-oriented university in Hangzhou, China, supported by public a…

The University of Toronto now recruiting for the Eric and Wendy Schmidt AI in Science Postdoctoral Fellowship. Valued at $85,000 CDN per year.

The Integrated Science Lab (IceLab) offers three postdoctoral scholarships that will be affiliated with one of six possible multidisciplinary projects

Aim to a global leader in materials science research society, and create a diverse, innovative, and collaborative research environment

Department of Frontier Materials, Suzhou Laboratory

The Max Planck Institute of Molecular Cell Biology and Genetics (MPI-CBG) in Dresden, Germany, is seeking an outstanding individual to fill the po…

Max Planck Institute of Molecular Cell Biology and Genetics