The absolute intensities of all 13 C NMR signals detected in all spectra were collected by automatic peak picking, after spectra processing (manually phased, baseline corrected, and referenced by setting the central resonance of DMSO-d6 at δ 39.80 ppm) using the TOPSPIN 4.1.3 software (Bruker, Rheinstetten, Germany). Each peak list was converted into a text file contain-ing peak positions and absolute peak intensities. Peaks positions from the whole set of spectra were then aligned by an algorithm written in the Python language. The principle was to divide the 13 C spectral width (from 0 to 240 ppm) into regular bins of 0.2 ppm width and to place the absolute intensity of each 13 C peak into the corresponding bin. The bins without any signal, regardless of fraction, were removed from the bin list. The resulting global table contains 12 columns, each corresponding a CPC fraction, and 403 rows, corresponding to the NMR spectral buckets for which at least one 13 C NMR peak was detected in at least one spectrum. This table was imported into the PermutMatrix software (version 1.9.3, LIRMM, Montpellier, France) for hierarchical clustering analysis (HCA). This operation reordered the rows of the global peak table so that similar rows, corresponding to similar chromatographic emergence profiles, were grouped together and lead to define clusters of chemical shifts values, with ideally one cluster defined per compound contained in the EtOAc extract of Larix decidua bark. The similarity between table rows was measured by the Euclidian distance, and data agglomeration was performed with the Ward's method. The resulting clusters of 13 C NMR chemical shifts were visualised as dendrograms on a heat map. Each 13 C NMR chemical shift cluster obtained from HCA was manually submitted to the VersaDB-generated database containing the structures and corresponding predicted NMR chemical shifts values (nmrshiftdb2) of 2790 natural metabolites found in the Pinaceae family (April 2022), via the structure search engine from ACD/NMR Workbook Suite 2012 (ACD/Labs, Ontario, Canada). This dereplication procedure was described in a previous article.(13) A 13 C NMR chemical shift difference between the predicted and experimental spectra was tolerated between 2 and 3 ppm, and the minimum number of 13 C query shift values to match was set at about 80% of the number of chemical shifts in the cluster. Finally, each structure proposal provided by the database query and its associated relative configurations was confirmed by interpretation of 1D and 2D NMR data ( 1 H, 13 C, HSQC, HMBC, COSY).