<figure><img src="image_1.png" /></figure>

HAL
open science

# Soluble organic matter Molecular atlas of Ryugu reveals cold hydrothermalism on C-type asteroid parent body

Philippe Schmitt-Kopplin, Norbert Hertkorn, Mourad Harir, Franco Moritz,

Marianna Lucio, Lydie Bonal, Eric Quirico, Yoshinori Takano, Jason

Dworkin, Hiroshi Naraoka, et al.

► To cite this version:

Philippe Schmitt-Kopplin, Norbert Hertkorn, Mourad Harir, Franco Moritz, Marianna Lucio, et al. Soluble organic matter Molecular atlas of Ryugu reveals cold hydrothermalism on C-type asteroid parent body. Nature Communications, 2023, 14 (1), pp.6525. 10.1038/s41467-023-42075-y. hal-04420427

HAL Id: hal-04420427

https://hal.science/hal-04420427v1

Submitted on 30 Jan 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

# Soluble organic matter Molecular atlas of Ryugu reveals cold hydrothermalism on C-type asteroid parent body

Received: 15 February 2023

Accepted: 19 September 2023

Published online: 16 October 2023

<figure><img src="image_4.png" /></figure>

Check for updates

Philippe Schmitt-Kopplin Ⓔ<fnref n="1" />,²,³✉, Norbert Hertkorn Ⓔ², Mourad Harir², Franco Moritz², Marianna Lucio², Lydie Bonal⁴, Eric Quirico Ⓔ⁴, Yoshinori Takano Ⓔ⁵, Jason P. Dworkin Ⓔ⁶, Hiroshi Naraoka Ⓔ⁷, Shogo Tachibana Ⓔ⁸,⁹, Tomoki Nakamura¹⁰, Takaaki Noguchi Ⓔ¹¹, Ryuji Okazaki⁷, Hikaru Yabuta Ⓔ¹², Hisayoshi Yurimoto Ⓔ¹³, Kanako Sakamoto Ⓔ⁹, Toru Yada Ⓔ⁹, Masahiro Nishimura⁹, Aiko Nakato Ⓔ⁹, Akiko Miyazaki Ⓔ⁹, Kasumi Yogata Ⓔ⁹, Masanao Abe⁹, Tomohiro Usui Ⓔ⁹, Makoto Yoshikawa⁹, Takanao Saiki⁹, Satoshi Tanaka⁹, Fuyuto Terui⁹, Satoru Nakazawa Ⓔ⁹, Tatsuaki Okada Ⓔ⁹, Sei-ichiro Watanabe Ⓔ¹⁴, Yuichi Tsuda⁹ & Hayabusa2-initial-analysis SOM team*

The sample from the near-Earth carbonaceous asteroid (162173) Ryugu is analyzed in the context of carbonaceous meteorites soluble organic matter. The analysis of soluble molecules of samples collected by the Hayabusa2 spacecraft shines light on an extremely high molecular diversity on the C-type asteroid. Sequential solvent extracts of increasing polarity of Ryugu samples are analyzed using mass spectrometry with complementary ionization methods and structural information confirmed by nuclear magnetic resonance spectroscopy. Here we show a continuum in the molecular size and polarity, and no organomagnesium molecules are detected, reflecting a low temperature and water-rich environment on the parent body approving earlier mineralogical and chemical data. High abundance of sulfidic and nitrogen rich compounds as well as high abundance of ammonium ions confirm the water processing. Polycyclic aromatic hydrocarbons are also detected in a structural continuum of carbon saturations and oxidations, implying multiple origins of the observed organic complexity, thus involving generic processes such as earlier carbonization and serpentinization with successive low temperature aqueous alteration.

<footnote n="1">¹Technische Universität München, Analytische Lebensmittel Chemie, Maximus-von-Forum 2, 85354 Freising, Germany. ²Helmholtz Munich, Analytical Bio-
GeoChemistry, Ingolstaedter Landstraße 1, 85764 Neuherberg, Germany. ³Max Planck Institute for Extraterrestrial Physics, Gießebachstraße 1, 85748
Garching bei München, Germany. ⁴Université Grenoble Alpes, CNRS, CNES, IPAG, 38000 Grenoble, France. ⁵Biogeochemistry Research Center (BGC), Japan
Agency for Marine-Earth Science and Technology (JAMSTEC), 2-15 Natsushima, Yokosuka 237-0061, Japan. ⁶Solar System Exploration Division, NASA
Goddard Space Flight Center, Greenbelt, Maryland 20771, USA. ⁷Department of Earth and Planetary Sciences, Kyushu University, Motooka 744, Nishiku,
Fukuoka 819-0395, Japan. ⁸Tokyo Organization for Planetary and Space Science, University of Tokyo, Bunkyo-ku, Tokyo 113-0033, Japan. ⁹Institute of Space
and Astronautical Science, Japan Aerospace Exploration Agency (ISAS/JAXA), Sagamihara 252-5210, Japan. ¹⁰Department of Earth Material Science, Tohoku
University, Aoba-ku, Sendai 980-8578, Japan. ¹¹Division of Earth and Planetary Sciences, Kyoto University, Kyoto 606-8502, Japan. ¹²Department of Earth and
Planetary Sciences, Hiroshima University, Higashi-Hiroshima, Hiroshima 739-8526, Japan. ¹³Department of Earth and Planetary Sciences, Hokkaido University,
Kita-ku, Sapporo 060-0810, Japan. ¹⁴Graduate School of Environment Studies, Nagoya University, Nagoya 464-8601, Japan. *A list of authors and their
affiliations appears at the end of the paper. ✉e-mail: schmitt-kopplin@tum.de</footnote>
The soluble organic matter (SOM) in carbonaceous chondrites is chemically highly diverse and integrates in its chemistry much of the history on temperature events and water-rock interaction on the parent body¹⁻³. Two touch-down samplings by the Hayabusa2 spacecraft in February and July 2019 enabled the collection of surface (samples stored in Chamber A) and possibly sub-surface materials (samples stored in Chamber C) of the near-Earth C₉-Type asteroid (162173) Ryugu. Samples were returned to Earth on December 6, 2020 and with about 5.4 g material that constitutes the fourth returned sample from extraterrestrial bodies following Apollo, Luna and recently Chang'e 5 from the Moon, Stardust from comet 81P/Wild2 and Hayabusa from the near-Earth S-type asteroid Itokawa, respectively. These material from enabled first time the comparison of our knowledge on meteorites to a C-type asteroid⁴,⁵. Ryugu is considered as a primitive small Solar-System body and its low-albedo reflectance spectrum in the wavelength range of 1.8 to 3.0 µm⁶ and spectroscopic similarity⁷ made it from observation a good candidate in being a parent body of carbonaceous chondrites (CCs), rich in water and organics.

The analysis of the Ryugu material showed the dominance of hydrous silicate minerals including serpentine and saponite associated with dolomite, pyrrhotite and magnetite, reflecting episode of aqueous alteration⁸ and suggesting similarities with CI chondrites⁴,⁸,⁹. Element analysis confirmed Ryugu with the richest C, N, and H concentrations compared to various types of CCs with signatures comparable to the observed falls Ivuna and Orgueil CI chondrites⁴,¹⁰. The analyzed elements content of C, H, N, S, and pyrolyzable O (not including O of anhydrous silica) was found to be ~21.3 wt% including the carbonates, the sulfides as well as the soluble organic matter (SOM) and the macromolecular insoluble organic matter (IOM)⁵. Targeted organic analysis of Ryugu A0106 sample revealed ~15 amino acids, including many non-protein amino acids in an approximately racemic mixture indicating a non-biological origin. Their abundances were different to Orgueil CI, reflecting different chemosynthetic or alteration pathway conditions on the asteroid parent bodies⁴. High abundances of β-alanine, γ-amino-*n*-butyric acid, and δ-aminovaleric acid in the Ryugu samples indicated a higher peak temperature than the CI chondrites as these compounds are typical in CV and CO chondrites processed at elevated temperatures of up to ~300 °C⁵,¹¹,¹². Series of aliphatic amines (i.e., methylamine (CH₃NH₂), ethylamine (C₂H₅NH₂) and isopropylamine, *n*-propylamine (C₃H₇NH₂)) and short chained carboxylic acids were analyzed in Ryugu samples as well⁵. Further targets were polycyclic aromatic hydrocarbons (PAH) as they are ubiquitous in presolar synthesis in interstellar locations¹³ and the FTIR spectra of Ryugu samples were close to interstellar PAHs spectra, suggesting incorporation of interstellar PAHs to Ryugu during its accretion⁵. Four-ringed fluoranthene, pyrene, chrysene/triphenylene and methylated fluoranthene and pyrene as well as smaller PAHs as naphthalene, phenanthrene and anthracene were analyzed and their relative abundances were attributed to differential alteration survivals and solubilities in parent body fluids⁵. Alkylpyridines and alkylimidazoles (aromatic N-heterocycles) and alkylpiperidines (aliphatic N-heterocycles)⁵,¹⁴ were targeted in Ryugu samples and showed a different profile as found earlier in CM2 Murchison¹⁵,¹⁶ reflecting a different redox condition on the parent body.

The returned surface sample of Ryugu A0106 enable now a direct analysis of SOM from the surface of the parent body and the comparison with meteoritic SOM. Here, we report the molecular characteristics of Ryugu's SOM from non-targeted organic analysis. More specifically we compared the ultrahigh-resolution mass spectroscopy analysis profile of the methanol extracts to 36 CCs that experienced moderate temperatures and partial aqueous alteration to position Ryugu's organic diversity in the context of possible hydrothermalism.

# Results and discussion

## The molecular atlas of Ryugu SOM

The Ryugu A0106 sample was extracted sequentially with solvents of increasing polarity, starting with hexane, dichloromethane, methanol and water to generate/isolate chemical fractions of SOM with molecules of increasing polarity. We first optimized injection flow rates to generate maximum scan count and signal-to-noise response to fully use the minimum amount of 5 mg of available extracts for a maximal information output in mass spectrometry and NMR-spectroscopy. We used electrospray ionization in negative (ESI(-)) and positive (ESI(+) modes as well as atmospheric pressure photoionization in positive mode (APPI(+) in direct injection Fourier transform-ion cyclotron resonance mass spectrometry (DI-FT-ICR/MS) as previously tested and demonstrated powerful method on many meteorites earlier³,¹⁷. We also confirmed that the methanol extract of the sequential extraction procedure and the direct extraction of Murchison CM2 and Aguas Zarcas CM2 fragments converged in the same signal and compositional profiles, making these two methanol extracts directly comparable. One third of the small solution volume available was used for nuclear magnetic resonance (NMR) spectroscopy experiments.

All analyzed FTICR-MS spectra showed very high signal density in the mass range from m/z 120 to 700 very comparable when analyzing CCs, reflecting that rapidly collected CCs are not biased by terrestrial contamination and weathering. The methanol extract in ESI(-) showed regular intense signals of increased intensity that could be assigned to polythiols (disulfate) with a number of sulfur from S3 to S9 (Fig. S1). The less polar hexane and dichloromethane extracts further showed the presence of polythiols as monosulfonates and of polysulfanes with sulfur numbers S3-S8. From all CCs analyzed in this study (Table S1), such intense profiles were observed only in methanol extract of recently fallen Tarda C2-ung material¹⁸. Each nominal mass always showed multiple signals with a repeating pattern of mass spacing of 36.3845 mDa following a compositional substitution of CH₄ with respect to O reflecting overall structural regularity in chemical homologous series in the mixture¹. Moreover, the Gaussian distribution of the mass signals within each nominal mass represents the large number of possible isomers for each exact mass, changing from a higher oxygen degree and carbon unsaturation from lower to higher mass defect as reported¹⁹, making any detailed structural annotation nearly impossible. In ESI(+) and APPI(+) the same regularity pattern was observed as well, similarly to patterns earlier described with all meteorite SOMs³,¹⁷.

In total, with all ionization modes and extractions we distinguished a molecular richness with more than 200,000 signals at S/N = 3. We could reduce them with conservative filter rules to 23,100 elementary compositions in the CHNOS monoisotopic atomic space. These results agree with our previous studies for various classes of CCs from a small molecule mass spectrometry perspective²,³, demonstrating their huge chemical diversity compared to the biological space²⁰,²¹. Figure 1A–C shows the final counting in elementary compositions assigned from FTICR-MS exact mass data and the overlaps of the compositional space between the different ionization methods and extraction solvents. The overall most abundant chemical families were the N-containing molecules (16,950 formula), followed by CHO/Na (8843 formula), CHOS/Na (5647 formula), CH (1056 formula) and CHS (165 formula). Most abundant molecules were extracted in polar protic methanol, followed by dichloromethane. Highly polar and highly apolar compounds extracted with water and hexane, respectively, showed the most specific compounds in terms of carbon oxidation state, nitrogen or sulfur content (Figs. S2, S4, and S5). In APPI(+) the oxygenation degree of CHO, CHNO and CHOS was similar and independent of the solvent polarity (2-3,0-atoms per molecule). The number of nitrogen, however, increased with increasing polarity indicating the involvement of nitrogen in polar functional groups such as primary or secondary amines increasing molecular basicity and thus
<figure><img src="image_5.png" /><figcaption>**Fig. 1** | Molecular atlas of the Ryugu soluble organic matter. A counting of the number of elementary composition retrieved from exact mass analysis in ESI(-), ESI(+) and APPI(+) ionization modes of the sequential extractions in hexane, dichloromethane, methanol and water. B Venn-diagram of the ionization mode signatures, C abundances of the various chemical families in the various extracts and ionization modes, involving, CH, CHO, CHN, CHS, CHOS, CHNO, CHNOS,</figcaption></figure>

CHONa, CHNONa, CHOSNa and CHNOSNa. D Network representation of the methanol extract using the most frequent mass transitions and showing the strong structural connectivity and the regular gradients in O/C, H/C and *m/z*. E frequency histogram of the exact mass differences (i.e., CH₂, H₂, O, CO₂, H₂O, HCN, NH, NH₃, SO, S, SO₂, SO₃, and SO₄).

promoting APPI(+) ionization (Fig. S3). CH type of molecules detected in APPI(+) showed a continuum in abundance of aliphatic to highly aromatic compounds bearing multiple aromatic rings (Fig. S6). Similarly, in ESI(+), the oxygenation degree is not much affected in the various polarity extracts, but was higher than in APPI (3-4 oxygen atoms per molecule) and particularly characterized by a higher relative abundance of nitrogen-containing ions, especially in the methanol extract. Oxygen being involved in acidic functional groups such as carboxyls, hydroxyls, the profiles are directly impacted by solvent polarity in ESI(-) and show a higher number of oxygen-rich molecules when moving to polar extraction solvents (Figs. S4 and S5). In addition, organosulfurs are more oxygenated toward polar solvent reflecting the extractability of molecules with sulfonations (-SO₃) compared to less oxygen bearing compounds bearing sulfidic functional groups (-SH). Only the methanol extract showed composition with multiple heteroatoms.

Data mining of the molecular atlas showed that the most abundant mass differences were found to be equivalent to methylation (CH₂), hydration (H₂O), hydrogenation (H₂), oxygenation (O), sulfurizations (SO₄, SO₃, SO₂, SO, S), nitrogenations (HCN, NH, NH₃) as presented earlier with Murchison CM2 as well¹. The search for removal or addition of these chemical functionalities' exact masses between all pairs of 23,100 molecular formulas resulted in the reconstruction of a mass difference network (MDiN) covering 154,578 pairs (edges) from among 19,675 molecular formulas (nodes). The remaining molecular formulas

could not be connected using these mass transitions. The most abundant mass difference transformations (CH₂, H₂, O and H₂O) amounted to 47% of all detected mass differences, connecting 20,232, 19,349, 17,111, and 16,620 pairs of compositions, respectively. Nitrogen (HCN, NH, NH₃) and sulfur-based mass differences (S, SO, SO₂, SO₃, SO₄) accounted for 22.6% and 20.0% of all transformations found, respectively. These structural continuum leads to differential solvent solubility and ionization efficiencies in the various ionization modes in FTICR-MS (Fig. S7).

¹H NMR spectra of A0106 methanol extract provided quantitative assessment of key proton and carbon chemical environments, and showed distinct groups of narrow resonances from δH ~ 0.5-8.5 ppm in which alkyl branching and termination by methyl and carboxylic groups were the defining molecular features, followed by lower abundance of aliphatic OCH₃ units, olefinic (coarsely dH <6.5 ppm) and aromatic unsaturation (coarsely dH <7 ppm) (Fig. 2). Olefins were incorporated in alkyl systems whereas aromatic CₐrHₐ units likely comprised polycyclic aromatic rings with ill-constrained admixture of aromatic carboxylic acids and nitrogen heterocycles. Ammonium (¹⁴NH₄⁺) showed a distinct 1:1:1 triplet at δH ~ 4.5 ppm with ¹JNH = 43 Hz and 1.2% relative abundance, exceeding other meteorite extracts. Ratios of methyl (δH < 1 ppm), non-methyl alkyl (δH ~ 1.0-1.65 ppm), and protons proximate to carboxylic acids HOOC-CHα-CHβ· (CHβα: δH ~ 2.15 ± 0.15 ppm; CHβ: δH ~ 1.65 ± 0.5 ppm) defined connectivity networks in two dimensional NMR spectra (Fig. S9) and provided reconstruction of average features of aliphatic branching in the
methanol extract with minor contributions from aliphatic oxygenation and olefinic unsaturation in A0106 SOM.

The dominance of hydrous silicate minerals (i.e., saponite, ser-
pentine) and the absence of chondrules, already oriented the classifi-
cation of the return material close to CI-type meteorites like Orgueil⁹.
We compared the methanol extract of the sample A0106 with the
methanol extracts of Murchison CM2 and Orgueil CI, two well-studied
meteorites representative of their chondritic classes (Fig. 3). All
extracts show full mass patterns over a range up to m/z 600. The van

<figure><img src="image_6.png" /><figcaption>**Fig. 2 |** Nuclear magnetic resonance spectroscopy of A0106. ¹H NMR spectrum (800 MHz, CD₃OD) of mHDOM extract 4A0106, with ¹H NMR section integrals of main substructures (cf. Table S3). Apart from certain sharp NMR resonances which possibly denote individual aliphatic CHO molecules, mHDOM comprised relevant background of broad bulk ¹H NMR resonances that represent a huge diversity of aliphatic branched alkyl and remotely oxygenated aliphatic groups (δ<sub>H</sub> ~ 0.6-2.7 ppm) and aliphatic OCH<sub>3</sub> molecules (δ<sub>H</sub> ~ 3.2-4.2 ppm); the overall curvature suggests higher relative abundance of aliphatic alcohols and ethers (δ<sub>H</sub> < 3.8 ppm) in comparison with e.g. aliphatic esters (δ<sub>H</sub> > 3.8 ppm). Ammonium (¹⁴NH<sub>4</sub><sup>+</sup>) showed a distinct 1:1:1 triplet at δ<sub>H</sub> ~ 4.5 ppm with ¹NH = 42.7 Hz (1.2% of ¹H NMR integral; Fig. S8a-c).</figcaption></figure>

Krevelen diagrams of A0106 show the same regular pattern as observed with Murchison that correspond to homolog chemical series but with lower abundance towards higher *m/z*. Orgueil showed the lowest abundance in mass signals and covers the same range in oxygen-rich region as A0106. Mass difference analysis normalized on CH2-content (Fig. 1E) showed a closest similarity in functionalization of A0106 with Orgueil CI than with Murchison CM2. Both A0106 and Orgueil share the same profile in the abundance of oxygenated hydrocarbon and organosulfur compounds. Murchison and A0106 share the same profile in heteroatom abundance, Orgueil having more molecules with many sulfur atoms.

Comparing the ¹H NMR spectra of the methanol extracts of A0106 with Murchison CM2 meteorite¹⁹ revealed major distinction. ¹H NMR resonances were much narrower in A0106 than those in Murchison, indicating faster transverse (T₂) relaxation in Murchison extracts. This reflects the relative proportions of background bulk NMR resonances across the entire shift range δ<sub>H</sub> in Murchison exceeding that of A0106. This indicates a larger chemodiversity of low abundance atomic environments in Murchison. Considerable background ¹H NMR resonances of pure aliphatic units (CCCH; δ<sub>H</sub> ~ 0.7-2.5 ppm; Fig. S8B) probably arose from significant contributions of alicyclic rings in Murchison extracts.

To explore the high-dimensional dataset, we computed a princi-
pal component analysis (PCA) using SIMCA 13.0.3.0 (Umetrics) based
on the intensities of all annotated m/z of the CHNOS compositional
space unit-variance (UV) scaled and log transformed (Fig. 4A). We
divided the data in a work (25 meteorite samples) and prediction set (11
meteorite samples) as shown in Table S1. The 36 CCs were chosen as to
reflect a range of warm to low-temperatures history based on their
thermal index from previous studies or known to have experienced
strong water alteration and hydro-thermalism²². The work and pre-
diction set were chosen randomly out of the 36 to check for the validity
of the model. The variance absorbed by the first two valid
components²³ is R²X(cum)=0.51. Additionally, the provision ability of
the model was measured by the Q² index (Q²(cum)=0.39). The PCA
score scatter plot (Fig. 4A) differentiates the chemical space in two
main regions. In one side, the type-TI CCs having seen low-temperature

<figure><img src="image_7.png" /><figcaption>**Fig. 3 |** FT-ICR-MS data evaluation of sample A0106 compared to carbonaceous meteorites. Comparison of the mass edited H/C ratios and van Krevelen diagrams of sample A0106 with Murchison CM2 and Orgueil CI. The bubble size is proportional to the intensity in the original ESI(-)-FTICR-MS mass spectra and the color</figcaption></figure>

codes correspond to defined chemical classes of CHO(blue), CHNO(orange),
CHOS(green) and CHONS(red). Histograms on the right are the abundance profiles
of the heteroatoms as well as the oxygenation profiles of the CHO, CHNO and CHOS
chemical families.
<figure><img src="image_8.png" /><figcaption>**Fig. 4 | Hydrothermal prediction of sample A0106 based on meteoritic soluble organic matter data.** A Principal component analysis (PCA) and model positioning the sample A0106 within 36 carbonaceous chondrites of TI-type (we defined the TI-type in Quirico et al. as a combination of XRD/mineral observation and Raman/IR data). All TI meteorites experiences aqueous alteration at temperatures <300 °C. In green the work dataset of the samples taken for the model and in blue the</figcaption></figure>

simulated positions form predicted dataset. B The loading plot represents the m/z values that differentiates an increase in water alteration (left) and temperature (right) represented in mass edited H/C ratios, accordingly; the bubble size is proportional to the intensity in the original ESI(-)-FTICR-MS mass spectra and the color codes correspond to defined chemical classes of CHO(blue), CHNO(orange), CHOS(green) and CHONS(red).

in various water alterations and in the other side, the chondrites of CM1, CR1, C1/2-ung, CM1/2, CI types that saw much water in their parent body history (*hydro-thermalism*). The most relevant masses of these two groups could be extracted based on their highest loading values and a representation in van Krevelen diagrams Fig. 4B.

The sulfurization of nitrogen compounds is mostly abundant at slightly increased temperatures reached in the typical TI-type meteorites²⁴. Lower peak temperatures and higher water alteration corresponding to cold hydro-thermalism show higher abundance of saturated/mono-unsaturated long chained aliphatic carboxylic acid (Fig. 4B).

## Implication of chemical diversity on processes on the parent body

All solvent extracts of the A0106 surface sample from Ryugu showed very high molecular diversity in the mass range from 120 to 600

atomic mass units. Compounds containing nitrogen and sulfur were preferentially extracted in a solvent polarity gradient based on their chemical functionality (leading to continuum in polarity) and were detected as a function of their ionization specificity in ESI(-), ESI(+) or APPI(+) accordingly (Figs. S4 and S5). The molecular formula were observed in a continuum of oxygenation degree and relative abundances and corresponded to chemical families such as CH, CHO, CHN, CHS, CHOS, CHNO, CHNOS, CHONa, CHNONa, CHOSNa and CHNOSNa. The visualization of this chemical structural continuum of complex kerogen-like organic matter (OM) of the Ryugu surface sample was especially possible with Mass Difference Networks (MDiNs) that were reconstructed from mass spectral m/z values or theoretical masses after computation of molecular formulas (Figs. 2 and S7). The extreme molecular richness of these extracts reflects the high importance of chemical processes leading to this diversity involving the specific chemistry of nitrogen and sulfur²⁵. We recently
linked the nitrogen compounds content between the carbonaceous meteorite Murchison and ices of dense molecular clouds and showed that the post-aqueous products of these organic residues produced from interstellar ice analogs share compositional similarities with SOM of the Murchison meteorite²⁶. In particular, the nitrogen containing molecules in some CCs may originate from organic-rich ices, inherited from the dense molecular cloud of our solar system. Hydrothermal processing of this material may release ammonium ions as seen in the A0106 sample. A suite of amino acids and aliphatic amines were recently analyzed in the same A0106 samples²⁷ and showed different profiles than in CI, indicating multiple formation mechanisms on Ryugu's parent body. More recently, nucleobases were described in various Murchison, Murray and Tagish Lake meteorites showing that some of these nitrogen-rich photochemical derivatives produced in the interstellar medium could have been incorporated into asteroids during solar system formation²⁸.

We did not observe any organic compounds containing magnesium (i.e. CHOMg, CHOSMg), which were described in various CCs with increasing temperature history. The absence of organomagnesium compounds was recently shown in the recent CI-ung fall Flensburg²⁹ and in the C2-ung Tarda³⁰, two meteorites that also saw high degree of water alteration. These characteristics reflect the low temperature hydrothermal processing (<-150 °C) involved on the parent body²,³, as well as the severe aqueous alteration and hydrolysis of the whole system. These results on the SOM are in line with the reported studies on insoluble organic matter (IOM) showing that Ryugu's organic material was modified by aqueous alteration on the asteroid parent body³¹. Considering the abundance profile of thousands of oxygenated molecular formula of the methanol extract of 36 meteorites we were able to dress a predictive A0106 being highly water processed at low temperature such as observed in C2-ung Tarda, CI-ung, CM1/2, CM1, CRI or CI (Fig. 4).

Ryugu's molecular atlas is characterized by a continuum in oxidized heterocyclic organosulfur and organonitrogen compounds corresponding to a molecular array snapshot of ongoing abiotic reaction processes on the parent body, heavily impacted by water alteration. Abiotic molecular synthesis was described recently in the Martian meteorite Tissint and earlier in ALH 84001 as being the origin of the complex OM observed co-localized to evolving minerals³²,³³. Here as well, the chemical signature of the SOM may result from a complex history in carbonization and serpentinization, water-rock interactions generating and consuming water during olivine and pyroxenes dissolution and leaving behind serpentine and an array of oxygenated novel carbon species³⁴. Co-localization of carbon with mineral phases and coevolution within a redox homeostasis will oxidize newly synthesized hydrocarbons while mineral and/or organic electron donors in the mixtures will reduce³². These arrays of mixed and localized reactions will create further complex organic molecules in a structural continuum in mass, oxygenation, carbon saturation and heteroatom content. Our results from both ESI- and APPI-FTICR-MS show the abundance of PAHs and hydrocarbons (CH) in a wide molecular range as well as their oxygenation products (CHO). While a hydrogenation process may favor dearomatization i.e. of PAHs, water alteration and hydrolysis processes would rather lead to more polar oxygen-rich compounds to reach a thermodynamically stable steady state. Serpentinization processes are thermodynamically favored, exothermic and therefore can also generate local low heat of up to only 150–200 °C on small parent bodies, thereby activating an additional chemical process of carbonization to again generate oxidized carbon and water³⁵. Finding no organomagnesium compounds also confirms the high water alteration and hydrolysis processes affecting the organic matter in the system as essential processes in organic evolution on Ryugu. The two described geochemical processes may be thermodynamically coupled and involve closely the mineral and organic phases within their complex coevolution and both occur

during the water-rock interaction. In this way, sulfurization processes may also occur involving the reactive sulfur species generated during sulfur mineral evolution. These processes are described in terrestrial environments and the incorporation of reduced inorganic sulfur (S) into OM in anoxic environments is considered a vital pathway by which OM is preserved for millions of years and the iron sulfide precipitation is believed to compete with OM preservation by sulfurization³⁶. The sulfur chemistry continuum in CHOS and CHNOS observed in sample A0106 reflects similar abiotic processes on Ryugu parent body.

The high molecular diversity of SOM in the Ryugu sample A0106 was comparable to meteorite SOM having seen high water alteration processes at low temperatures. The methanol-extracted SOM profile are consistent with A0106 being closely related to CM1/2, CI or CM1 type of meteorites. The presence of homologous chemical series in nitrogen-containing organics may support an origin from organic ices, inherited from the dense molecular cloud of our solar system. The molecular atlas of Ryugu consists of a continuum in molecular size, oxygenations, and sulfurization of both hydrocarbons and nitrogen-containing organics and reflect a complex series of multiple and successive sources in abiotic organic synthesis possibly involving redox-based water-rock reaction, water dissociation (e.g. water radiolysis), with subsequent serpentinization carbonation and further organic solvolysis. Low amounts of organomagnesium confirmed low processing temperatures around 150 °C maximum in accordance to serpentinization processes. Thus, the SOM of Ryugu C-type asteroid surface is a source of organic diversity and complexity and may be regarded as a cradle of evolving prebiotic molecules in the Solar system, a source of molecular precursors of life itself.

## Methods

### Meteorite samples used for analysis and structural comparison

The list of the meteorites analyzed in the study and used for the model in Fig. 4 is given in Table S1.³⁷,³⁸

### Sample description

The Ryugu sample extracted was sample A0106, the surface samples stored in Chamber A. It was extracted sequentially with solvents of increasing polarity, starting with hexane, dichloromethane, methanol and water to generate/isolate chemical fractions of SOM with molecules of increasing polarity. This was done at JAXA and samples were distributed to various Labs for analysis as described in ref. 5.

### Fourier transform ion cyclotron resonance mass spectrometry

The experimental study was performed on a high-field FT-ICR mass spectrometer from Bruker Daltonics with a 12-T magnet from Magnex¹. A time domain transient with 4 MWords was obtained and Fourier-transformed into a frequency domain spectrum. The frequency domain was afterward converted to a mass spectrum by the SolarIX Control program of Bruker Daltonics. The ion excitations were generated in broadband mode (frequency sweep radial ion excitation) and 3000 scans were accumulated for each mass spectrum in a mass range of 147–1000 amu. Ions were accumulated for 300 ms before ICR ion detection. The pressure in the quadrupole/hexapole and ICR vacuum chamber was $3 \times 10^{-6}$ and $6 \times 10^{-10}$ mbar, respectively. For CID-MS/MS, ions were accumulated for 3 s.

The ESI source (Apollo II; Bruker Daltonics) was used in negative ionization mode and the APPI source in positive mode. The methanolic solutions were injected directly into the ionization source by means of a microliter pump at a flow rate of 120 µL h⁻¹. in ESI and at 500 µL h⁻¹ APPI. A source heating temperature of 200 °C was maintained and no nozzle-skimmer fragmentation was performed in the ionization source. The instrument was previously externally calibrated by using arginine negative cluster ions (5 mg L⁻¹ arginine in methanol).
FT-ICR mass spectra with *m/z* from 95 to 1000 amu were calibrated externally and internally to preclude alignment errors. Subsequently, the mass spectra were exported to peak lists at a signal-to-noise ratio ≥3. Elemental formulas were calculated combinatorically within a mass accuracy window of ±0.2 ppm for each peak in batch mode by an in-house software tool and validated via the senior-rule approach/cyclomatic number, assuming valence 2 for S and valence 4 (coordination number) for Mg²⁺.

FT-ICR-MS analysis enables highly resolved (*R* > 10⁶ at *m/z* 200) and accurate chemical mass analysis of electrospray generated ions within a 200 ppb error window over a wide mass range from *m/z* 100–1000. The weight of the ions is measured with a precision lower than the mass of an electron (Δ*m/z* = 0.0003 amu) and the specific signals can be differentiated with the same mass precision due to the ultrahigh resolution. These exact masses of the ions can routinely be converted into unique compositional formula bearing the light elements C-, H-, N-, O-, S-, Mg (or any other element in target), also taking account of their possible Cl and Na adducts and of their natural isotopic abundance.

The solvent extracts generated thousands of individual signals that were converted into elementary compositions (formula); these are all represented in van Krevelen type of diagrams (H/C vs. O/C) or related (H/C vs *m/z*) in which each formula is represented by a dot (the size of the dot is proportional to its abundance) as a projection of the relative oxygenation degree (O/C) and saturation degree (H/C) for various classes of compound types (CHO, CHNO, CHOS, CHNOS, CH, CHS, CHN, and selected Na⁺-adducts). Aromaticity equivalents Xc were calculated as reported earlier³⁹ and plotted as in Fig. S6.

## Nuclear magnetic resonance spectroscopy

A Bruker Avance III spectrometer and TopSpin 3.5/PL7 software were used to acquire nuclear magnetic resonance (NMR) spectra of aqueous extracts of Hayabusa methanolic extracts of which ~30 µg were exchanged with CD₃OD (99.96 % ²H, Aldrich) three times on a vacuum line¹⁹; after sealing in a 1.7 mm Bruker match tube, the sample was centrifuged without visible solid. Murchison extract A5¹⁹ was used here for comparison to sample A0106.

A cryogenic inverse geometry 5 mm z-gradient ¹H/¹³C/¹⁵N/³¹P QCI probe (*B*₀ = 18.8 T) was used for 1D ¹H NMR and proton-detected 2D NMR spectra. Transmitter pulses were at ~10 µs for ¹H and ¹³C. The one bond coupling constant ¹J(CH) used in 2D ¹H,¹³C DEPT-HSQC spectra (*hsqcedetgspsp2.2*) was set to 145 Hz; other conditions: ¹³C 90 deg decoupling pulse, GARP (70 µs); 50 kHz WURST 180 degree ¹³C inversion pulse (Wideband, Uniform, Rate, and Smooth Truncation; 1.2 ms); F2 (¹H): spectral width of 11160.7 Hz (13.95 ppm); 1.25 s relaxation delay; F1 (¹³C): SW = 36052 Hz (180 ppm). HSQC-derived NMR spectra were computed to a 8192 × 1024 matrix. Gradient (1 ms length, 450 µs recovery) and sensitivity enhanced sequences were used for all 2D NMR spectra. Absolute value COSY, and phase sensitive echo-antiecho TOCSY spectra (*cosygpmfppqf*, *dipsi2etgpsi*) used a spectral width of 9615.4 Hz and were computed to a 16,384 × 2048 matrix; other NMR acquisition conditions are given in Table S2.

## Data availability

All data from the mission are available at the DARTS archive www.darts.isax.jaxa.jp/planet/project/hayabusa2/ and on the Small Bodies Node of the NASA Planetary Data System https://pds-smallbodies.astro.umd.edu/data sb/missions/hayabusa2/. The samples of Ryugu are curated by the JAXA Astromaterials Science Research Group; distribution for analysis is through an Announcement of Opportunity available at https://jaxa-ryugu-sample-ao.net. The FTICR-MS raw data can be made available in contacting the corresponding author and can be adapted depending on further utilization.

## References

1. Ph. Schmitt-Kopplin, Z. et al. High molecular diversity of extra-terrestrial organic matter in Murchison meteorite revealed 40 years after its fall. *Proc. Natl Acad. Sci. USA* **107**, 2763–2768 (2010).

2. Ruf, A. et al. Previously unknown class of metalorganic compounds revealed in meteorites. *Proc. Natl Acad. Sci. USA* **114**, 2819–2824 (2017).

3. Matzka, M. et al. Thermal history of asteroid parent bodies is reflected in their metalorganic chemistry. *Astrophys. J. Lett.* **915**, L7 (2021).

4. Yokoyama, T. et al. Samples returned from the asteroid Ryugu are similar to Ivuna-type carbonaceous meteorites. *Science* **379**, eabn7850 (2023).

5. Naraoka, H. et al. Organic compounds from Ryugu – prebiotic molecules in a C-type asteroid. *Science* **379**, eabn9033 (2023).

6. Kitazato, K. et al. The surface composition of asteroid 162173 Ryugu from Hayabusa2 near-infrared spectroscopy. *Science* **364**, 272–275 (2019).

7. Hiroi, T., Pieters, C. M., Zolensky, M. E. & Lipschutz, M. E. Evidence of thermal metamorphism on the C, G, B, and F asteroids. *Science* **261**, 1016–1018 (1993).

8. Nakamura, T. et al. Formation and evolution of Cb-type asteroid Ryugu: direct evidence from returned samples. *Science* **379**, eabn8671 (2022).

9. Yada, T. et al. Preliminary analysis of the Hayabusa2 samples returned from C-type asteroid Ryugu. *Nat. Astron.* **6**, 214–220 (2022).

10. Kerridge, J. F. Carbon, hydrogen and nitrogen in carbonaceous chondrites: abundances and isotopic compositions in bulk samples. *Geochim. Cosmochim. Acta* **49**, 1707–1714 (1985).

11. Nazrul, I. M., Takeo, K. & Kensei, K. Reaction of amino acids in a supercritical water-flow reactor simulating submarine hydrothermal systems. *Bull. Chem. Soc. Jpn.* **76**, 1171–1178 (2003).

12. Li, J. & Brill, T. B. Spectroscopy of hydrothermal reactions, part 26: kinetics of decarboxylation of aliphatic amino acids and comparison with the rates of racemization. *Int. J. Chem. Kinet.* **35**, 602–610 (2003).

13. Allamandola, L. J., Tielens, A. G. G. M. & Baker, J. R. Polycyclic aromatic hydrocarbons and the unidentified infrared emission bands: auto exhaust along the Milky Way! *Astrophys. J.* **290**, L25–L28 (1985).

14. Naraoka, H. & Hashiguchi, M. Distinct distribution of soluble N-heterocyclic compounds between CM and CR chondrites. *Geochem. J.* **53**, 33–40 (2019).

15. Orthous-Daunay, F.-R. et al. Ultraviolet-photon fingerprints on chondritic large organic molecules. *Geochem. J.* **53**, 21–32 (2019).

16. Isa, J. et al. Aqueous alteration on asteroids simplifies soluble organic matter mixtures. *Astrophys. J. Lett.* **920**, L39 (2021).

17. Hertzog, J., Naraoka, H. & Schmitt-Kopplin, P. Profiling Murchison soluble organic matter for new organic compounds with APPI- and ESI-FT-ICR MS. *Life* **9**, 48 (2019).

18. Schmitt-Kopplin, P. et al. High molecular diversity and structural complexity revealed with ultrahigh resolution mass spectroscopy and nuclear magnetic resonance spectroscopy of Ryugu samples. *LPI Contrib.* **2678**, 1524 (2022).

19. Hertkorn, N., Harir, M. & Schmitt-Kopplin, P. Nontarget analysis of Murchison soluble organic matter by high-field NMR spectroscopy and FTICR mass spectrometry. *Magn. Reson. Chem.* **53**, 754–768 (2015).

20. Schmitt-Kopplin, P. et al. Systems chemical analytics: introduction to the challenges of chemical complexity analysis. *Faraday Discuss.* https://doi.org/10.1039/c9fd00078j (2019).

21. Meckenstock, R. U. et al. Water droplets in oil are microhabitats for microbial life. *Sci. Mag.* **345**, 673–676 (2014).
22. Quirico, E. et al. Prevalence and nature of heating processes in CM and C2-ungrouped chondrites as revealed by insoluble organic matter. *Geochim. Cosmochim. Acta* **241**, 17–37 (2018).

23. Bair, E., Hastie, T., Paul, D. & Tibshirani, R. Prediction by supervised principal components. *J. Am. Stat. Assoc.* **100**, 473 (2006).

24. Ruf, A. et al. Sulfur ion irradiation experiments simulating space weathering of Solar System body surfaces compound formation. *Astron. Astrophys.* **655**, A74 (2021).

25. Ruf, A. et al. Organosulfur compounds formed by sulfur ion bombardment of astrophysical ice analogs: implications for Moons, Comets, and Kuiper Belt Objects. *Astrophys. J. Lett.* **885**, 2G (2019).

26. Danger, G. et al. From molecular clouds to chondrites: exploring the link between molecular cloud ices and chondritic organics. *Nat. Commun.* **12**, 3538 (2021).

27. Parker, E. et al. Extraterrestrial amino acids and amines identified in asteroid Ryugu sample returned by the Hayabusa2 mission. *Geochem. Cosmochem. Acta* **347**, 42 (2023).

28. Oba, Y. et al. Identifying the wide diversity of extraterrestrial purine and pyrimidine nucleobases in carbonaceous meteorites. *Nat. Commun.* **13**, 2008 (2022).

29. Bischoff, A. C. et al. The old, unique C1 chondrite Flensburg - insight into the first processes of aqueous alteration, brecciation, and the diversity of water-bearing parent bodies and lithologies. *Geochim. Cosmochim. Acta* **293**, 142–186 (2021).

30. Aoudjehane, H. C. et al. Tarda an unusual carbonaceous chondrite meteorite fall from Morocco. In *52nd Lunar and Planetary Science Conference LPI Contrib. No. 2548 (NASA/ADS, 2021)*.

31. Yabuta, H. et al. Macromolecular organic matter in sample of the asteroid (162173) Ryugu. *Science* **379**, eabn9057 (2023).

32. Ph. Schmitt-Kopplin, M. et al. Complex carbonaceous matter in Tissint martian meteorites give insights into the diversity of organic geochemistry on Mars. *Sci. Adv.* **9**, eadd6439 (2023).

33. Steele, A. et al. Organic synthesis associated with serpentinization and carbonation on early Mars. *Science* **375**, 6577 (2022).

34. Schwarzenbach, E. M., Vogel, M., Früh-Green, G. L. & Boschi, C. Serpentinization, carbonation, and metasomatism of ultramafic sequences in the Northern Apennine Ophiolite (NW Italy). *AGU* **126**, 5 (2021).

35. Sforna, M. C. et al. Abiotic formation of condensed carbonaceous matter in the hydrating oceanic crust. *Nat. Commun.* **9**, 5049 (2018).

36. Aubakar, Y., Taylor, K. G., Coker, V., Vogelius, R. A. & van Dongen, B. E. Fundamental controls on organic matter preservation in organic- and sulfur-rich hydrocarbon source rocks. *Mar. Pet. Geol.* **141**, 105684 (2022).

37. Rubin, A. E. Petrography of refractory inclusions in CM2. 6 QUE 97990 and the origin of melilite-free spinel inclusions in CM chondrites. *Meteorit. Planet. Sci.* **42**, 1711 (2007).

38. Alexander, C. M. D., Howard, K. T., Bowden, R. & Fogel, M. L. The classification of CM and CR chondrites using bulk H, C and N abundances and isotopic compositions. *Geochim. Cosmochim. Acta* **123**, 244 (2013).

39. Yassine, M., Harir, M., Dabek-Zlotorzynska, E. & Schmitt-Kopplin, P. Structural characterization of organic aerosol using Fourier transform ion cyclotron resonance mass spectrometry: aromaticity equivalent approach. *Rapid Commun. Mass Spec.* **28**, 2445–2454 (2014).

## Acknowledgements

The Hayabusa2 project has been led by JAXA (Japan Aerospace Exploration Agency) in collaboration with DLR (German Space Center) and CNES (French Space Center), and supported by NASA and ASA

(Australian Space Agency). We thank all of the members of the Hayabusa2 project for their technical and scientific contributions. We thank Dr. Laurence Garvie from ASU (Arizona State University) and the Buseck Center Meteorite Studies for some selected fragments of meteorites used in this study. This research is partly supported by the Japan Society for the Promotion of Science (JSPS) under KAKENHI grant numbers JP20H00202, JP20H05846, JP20K20485, JP20K14549, JP21J00504, JP21H01203, JP21H04501, and JP21KK0062. J.P.D., J.C.A., E.T.P., D.P.G., H.L.M., J.E.E., and H.V.G. are grateful to NASA for support of the Consortium for Hayabusa2 Analysis of Organic Solubles. This research is also funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) – Project-ID 364653263 – TRR 235 (CRC 235).

## Author contributions

P.S.-K. designed this research. P.S.-K. and N.H. conducted experiments and analyzed data in cooperation with M.H., M.L., F.M., L.B., and E.Q. P.S.-K. wrote the paper. S.T., H. Yurimoto, T. Nakamura, T. Noguchi, R.O., H. Yabuta, and H.N. administered the initial analysis with J.D., Y. Takano, and M.A. T. Yada, M.N., K.Y., A.N., M.Y., A.M., and T.U. curated samples. K.S., T.O., S.N., F.T., S. Tanaka, and T.S. contributed to science operations of the spacecraft. S.-i.W. and Y. Tsuda administered the project. All authors discussed the results and commented on the manuscript.

## Funding

Open Access funding enabled and organized by Projekt DEAL.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41467-023-42075-y.

**Correspondence** and requests for materials should be addressed to Philippe Schmitt-Kopplin.

**Peer review information** *Nature Communications* thanks Marceau Lecasble and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. A peer review file is available.

**Reprints and permissions information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Open Access** This 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, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2023
Hayabusa2-initial-analysis SOM team

Philippe Schmitt-Kopplin✉️¹,²,³✉️, Norbert Hertkorn², Hiroshi Naraoka✉️⁷, Yoshinori Takano✉️⁵, Jason P. Dworkin✉️⁶, Kenji Hamase<fnref n="15" />, Aogu Furusho<fnref n="15" />, Minako Hashiguchi¹⁴, Kazuhiko Fukushima¹⁶, Dan Aoki¹⁶, José C. Aponte⁶, Eric T. Parker⁶, Daniel P. Glavin⁶, Hannah L. McLain⁶, Jamie E. Elsila⁶, Heather V. Graham⁶, John M. Eiler¹⁷, Alexander Ruf¹⁸, Francois-Regis Orthous-Daunay⁴, Junko Isa¹⁹, Véronique Vuitton⁴, Roland Thissen²⁰, Nanako O. Ogawa²¹, Saburo Sakai²¹, Toshihiro Yoshimura²¹, Toshiki Koga²¹, Haruna Sugahara⁹, Naohiko Ohkouchi²¹, Hajime Mita²², Yoshihiro Furukawa¹⁰, Yasuhiro Oba²³ & Shogo Tachibana✉️⁸,⁹

<footnote n="15">¹⁵Graduate School of Pharmaceutical Sciences, Kyushu University, Motooka 744, Nishi-ku, Fukuoka 819-0395, Japan. ¹⁶Graduate School of Bioagricultural Sciences, Nagoya University, Chigusa-ku, Nagoya 464-8601, Japan. ¹⁷Division of Geological and Planetary Sciences, California Institute of Technology, Pasadena, California 91125, USA. ¹⁸Excellence Cluster ORIGINS, Garching 85748, Germany. ¹⁹Earth-Life Science Institute (ELSI), Tokyo Institute of Technology, Meguro-ku, Tokyo 152-8550, Japan. ²⁰Université Paris-Saclay, CNRS, Institut de Chimie Physique, UMR8000, 91405 Orsay, France. ²¹Biogeochemistry Research Center (BGC), Japan Agency for Marine-Earth Science and Technology (JAMSTEC), 2-15 Natsushima, Yokosuka 237-0061, Japan. ²²Department of Life, Environment and Applied Chemistry, Fukuoka Institute of Technology, Higashi-ku, Fukuoka 811-0295, Japan. ²³Institute of Low Temperature Sciences (ILTS), Hokkaido University, Kita-ku, Sapporo 060-0810, Japan.</footnote>