ORIGINAL RESEARCH article

Front. Mar. Sci., 12 December 2017

Sec. Marine Biotechnology and Bioproducts

Volume 4 - 2017 | https://doi.org/10.3389/fmars.2017.00399

Sporochartines A–E, A New Family of Natural Products from the Marine Fungus Hypoxylon monticulosum Isolated from a Sphaerocladina Sponge

  • 1. Centre National de la Recherche Scientifique, Institut de Chimie des Substances Naturelles ICSN, Gif-sur-Yvette, France

  • 2. LEMAR, IRD, Centre National de la Recherche Scientifique, IFREMER, Université de Bretagne Occidentale, IUEM, Technopole Brest-Iroise, Plouzané, France

Abstract

Four new sporochartines B–E were isolated from the marine fungus Hypoxylon monticulosum CLL-205, isolated from a sponge belonging to the Sphaerocladina order and collected in Tahiti coast. Sporochartine A (1), the first representative of this family was previously isolated from the same fungus. The structures of sporochartines B–E were elucidated using 1D and 2D NMR, HRMS and IR data. Their configurations were established according to ROE correlations and comparison with the absolute configuration of sporochartine A (1) previously obtained from X-ray analysis. Sporochartines A–D (2–4) may be derived from endo Diels-Alderase type catalysis and sporochartine E (5) from an exo Diels-Alderase catalysis. The spatial conformation of sporochartines drastically influences the results of the cytotoxic bioassay against HCT-116, PC-3, and MCF-7 human cancer cell lines.

Introduction

The fungal Xylariaceae family includes more than 16 genera and 130 species (Sánchez-Ballesteros et al., ) and has been extensively investigated for the chemo diversity and biological activity of their metabolites (Stadler et al., , ). Among the 16 genera reported, Hypoxylon with 14 species is largely distributed in various marine and terrestrial habitats, and producing a large variety of bioactive compounds among which cohaerins (Quang et al., ; Surup et al., ), daldinins and daldinones (Quang et al., ; Gu et al., ), cytochalasin (Espada et al., ), fragiformin (Stadler et al., ), mitorubrinols (Quang et al., ), hypoxylonols (Fukai et al., ), hypoxylans (Kuhnert et al., ), hypoxyvermelhotins (Kuhnert et al., ), rickenyls (Kuhnert et al., ), rutilins (Quang et al., ), carneic acids (Quang et al., ), hymatoxins (Bodo et al., ; Borgschulte et al., ), malettinins (Angawi et al., ), hypoxysordarin (Daferner et al., ), lenormandins (Kuhnert et al., ), nodulisporic acids (Bills et al., ), schweinitzin A (Linh et al., ), truncatones (Sudarman et al., ), macrocyclic polyesters 15G256 family (Schlingmann et al., ), and sporothriolide (Krohn et al., ; Surup et al., ; Cao et al., ).

Sporothriolide belongs to the furofurandione family of natural compounds first published in 1994 (Krohn et al., ). This compound exhibits antifungal activity and benefits from substantial synthetic efforts (Sharma and Krishnnudu, ; Yu et al., ; Fernandes and Ingle, ; Ishihara et al., ). The name sporothriolide is related to Sporothrix sp. Hektoen and Perkins (strain 700), from which this compound was first isolated. Sporothrix genus belongs to a different ascomycete family, ophiostomataceae. The first report on sporothriolide in 1994 detailed both the structure and bioactivity of this product (Krohn et al., ). It shows that Sporothrix produces sporothriolide, dihydrosporothriolide, as well as various sporothriolide analogs with different side-chain length (canadensolide, discosiolide, avenaciolide, ethiosolide). The authors reported also the antifungal/herbicidal activities of these compounds (Krohn et al., ) (Figure 1).

Figure 1

). Middle: Monocyclic precursors reported in 2014 (Surup et al., ). Bottom: Deoxysporothric acid and sporochartine A reported in our previous work (Leman-Loubière et al., ).

Twenty years later, sporothriolide and dihydrosporothriolide were isolated from Hypoxylon monticulosum together with three monocyclic acid precursors: sporothric acid, isosporothric acid and dihydroisosporothric acid (Figure 1) (Surup et al., ). More recently, sporothriolide was isolated from Nodulisporium sp., an anamorph of Hypoxylon, and the herbicidal activity was confirmed (Cao et al., ).

In our previous contribution, we added to the scarce sporothriolide family two new compounds, deoxysporothric acid and a new complex architecture sporochartine A, combining sporothriolide and trienylfuranol A moieties. The trienylfuranol A was recently isolated from a different Hypoxylon submoniticulosum (Burgess et al., ).

In the present work we report four new sporochartines B to E. Their structures were elucidated using 1D and 2D NMR, HRMS, IR and comparison with sporochartine A data, for which the absolute configuration was previously established by X-Ray analysis. A Diels-Alderase type reaction is probably involved in the biosynthesis of the five isolated sporochartines, as discussed below.

The human cancer cell-lines cytotoxicity bioassay shows that the conformation of sporochartines has an impact on the biological activity.

Results and discussion

According to our previous work, sporochartine A (1) was obtained after 5 days cultivation of H. monticulosum CLL-205 in PDB broth (Leman-Loubière et al., ). By extending the cultivation of the same microorganism by a further 4 days, the ethyl actetate extract gives the HPLC chromatogram presented in Figure 2.

Figure 2

Sporochartines B–D were isolated as white powders by flash chromatography followed by semi-preparative HPLC. They had similar [M+H]+ HRESIMS molecular weights, molecular formula C24H34O6 and IR spectra compared with sporochartine A (1) (Table 2) (Leman-Loubière et al., ). The 1H and 13C NMR spectra of sporochartines B–D were similar to those of compound 1 (Tables 1, 2). Optical rotations , IR bands and HRESIMS are reported in Table 2.

Table 1

PositionSporochartine ASporochartine BSporochartine CSporochartine DSporochartine E
δH mult, (J in Hz)δH mult, (J in Hz)δH mult, (J in Hz)δH mult, (J in Hz)δH (J in Hz)
23.32, d (5.3)3.30, d (5.8)3.24, d (5.9)3.30, d (5.2)3.42, d (5.9)
55.10, dd (3.5, 5.3)5.13, dd (4.3, 5.9)5.01, dd (4.2, 5.8)5.09, dd, (3.6, 5.2)4.97, dd (4.1, 5.8)
64.45, m4.39, m4.40, m4.45, m4.44, m
71.78, m1.76, m1.76, m1.82, m1.81, m
1.91, m1.85, m1.85, m1.91, m1.89, m
81.45, m1.45, m1.45, m1.47, m1.44, m
91.37, m1.34, m1.34, m1.37, m1.36, m
101.30, m1.29, m1.29, m1.30, m1.29, m
111.30, m1.29, m1.29, m1.30, m1.28, m
120.89, t (6.8)0.88, t (6.9)0.89, t (6.9)0.89, t (7.0)0.88, t (6.8)
132.00, m2.04, m2.02, brdd (5.8, 14.0)1.96, m1.72, dd (9.5, 14.1)
2.13, m2.12, m2.26, dd (6.0, 14.6)
142.27, m2.25, m2.24, brd (21.4)2.26, brd (19.1)3.25, m
2.62, m2.79, m2.81, m2.60, m
155.94, brd (9.9)5.95, brd (10.9)5.95, brd (9.9)5.94, brd (9.9)6.16, d (10.2)
165.63, m5.54, brd (10.9)5.50, dq (2.0, 9.9)5.51, m5.56, ddd (2.2, 4.7, 5.5)
172.76, br t (5.8)3.23, br m3.19, br m2.80, br t (6.5)2.80, dd (5.5, 9.0)
185.65, dd (15.5, 7.2)5.82, ddd (1.5, 8.8, 15.4)5.66, dd (7.7, 15.4)5.65, dd (7.7, 15.6)5.65, dt (9.7, 16.9)
195.64, m5.76, dd (3.9, 15.4)5.67, dd (5.8, 15.4)5.57, dd (5.9, 15.6)5.27, d (10.0)
5.21, d (16.9)
204.16, t (4.8)4.19, m4.16, m4.20, t (5.3)3.22, dd (3.1, 9.7)
214.30, m4.27, m4.16, m4.07, quad (6.9)4.28, d (5.0)
221.55, m1.59, m1.60, m1.59, m1.53, dd (1.2, 5.9, 13.7)
2.39, m2.39, m2.39, q (6.5)2.40, dt (6.6, 12.6)2.39, ddd (6.6, 8.2, 14.0)
233.94, sext (6.3)4.07, m4.22, m4.26, sext (7.4)3.96, m
241.34, d (6.2)1.34, d (6.1)1.32, d (6.3)1.32, d (6.2)1.34, d (6.2)

1H NMR data for sporochartines A–E (Data acquired in CDCl3 at 500 MHz).

Table 2

PositionSporochartine ASporochartine BSporochartine CSporochartine DSporochartine E
δCδCδCδCδC
1171.8, C173.1, C172.8, C171.9, C172.7, C
250.6, CH47.2, CH47.2, CH50.9, CH49.7, CH
351.0, C51.0, C50.9, C49.2, C46.3, C
4177.0, C178.7, C178.6, C176.2, C175.6, C
578.2, CH78.7, CH78.7, CH78.0, CH76.4, CH
681.2 CH81.7 CH81.1 CH81.3 CH81.6, CH
728.9, CH228.9, CH229.0, CH229.0, CH229.0, CH2
825.5, CH225.3, CH225.6, CH225.5 CH225.5, CH2
929.1, CH229.1, CH229.1, CH229.2, CH229.1, CH2
1031.7, CH231.7, CH231.7, CH231.8, CH231.7, CH2
1122.7, CH222.9, CH222.9, CH222.7, CH222.7, CH2
1214.2, CH314.3, CH314.2, CH314.3, CH314.2, CH3
1320.9, CH226.9, CH227.1, CH220.9, CH224.1, CH2
1422.7, CH222.7, CH222.7, CH222.6, CH232.6, CH
15129.4, CH130.0, CH130.2, CH129.8, CH130.9, CH
16124.0, CH124.8, CH124.6, CH123.5 CH124.1, CH
1745.5, CH47.0, CH46.5, CH43.9, CH46.9, CH
18131.6, CH130.9, CH130.6, CH129.1, CH136.8, CH
19131.5, CH130.1, CH134.3, CH134.1, CH119.6, CH2
2083.3, CH82.8, CH84.8, CH84.9, CH87.0, CH
2174.8, CH73.9, CH77.3, CH77.1, CH72.5, CH
2242.6, CH242.4, CH242.3, CH241.2, CH243.4, CH2
2374.3, CH74.2, CH74.1, CH74.1, CH73.9, CH
2421.7 CH322.5 CH322.4 CH322.5 CH322.6, CH3
−57° (0.5, CHCl3)+72° (1.0, CHCl3)+93° (0.27, CHCl3)−152° (0.27, CHCl3)+51° (c 0.3, CHCl3)
[M+H]+ HRESIMS419.2423419.2423419.2431419.2429419.2433
IR3,521, 2,929, 2,859, 1,771, 1,452, 1,304, 1,175, 1,019 cm−13,468, 2,934, 2,863, 1,770, 1,303, 1,179, 1,071 cm−13,441, 2,957, 2,928, 2,858, 1,770, 1,454, 1,177, 1,017 cm−13,435, 2,932, 1,771, 1,178, 1,019 cm−13,501, 2,931, 1,766, 1,308, 1,187, 1,075 cm−1

13C NMR data for sporochartines A–E (Data acquired in CDCl3 at 125 MHz).

COSY and HMBC spectra, confirmed that sporochartines A–D had the same connectivities supporting similar planar scaffold (Figure 3). In addition, the common coupling constant of 15.4–15.6 Hz between H-18 and H-19 confirmed that the double bond C-18/C-19 is in E configuration.

Figure 3

Based on the previously reported absolute configuration of sporochartine A (1) and ROE correlations, we deduced the absolute configuration of sporochartines B–D (2–4) (Figure 4).

Figure 4

The common ROE correlations between H-2 and H-5 and between H-5 and H-6 requiring a cis orientation of these three protons was found in the sporochartine A–D. Therefore, the stereochemistry of the sporothriolide moiety was identical. Moreover, based on ROE correlations between H-20 and H-21 and between H-21 and H-23, the strerochemistry of the tetrahydrofurane moiety is also a common feature in sporochartines A–D.

For sporochartine B (2) (Figure 4), we did not observe ROE correlations between H-17 and H-2 and between H-17 and H-14b as in sporochartine A (1), while a new correlation is observed between H-17 and H-13. This data suggests that the carbon C-17 have opposite stereochemistry compared to 1 supporting a 3S,17R configuration of 2 (instead of 3S,17S in 1).

For sporochartine C (3) (Figure 4), the H-17/H-2 and H-17/H-14b correlations observed in sporochartine A (1) are absent. In addition, we observed a correlation between H-17 and H-13 and H-2 and H-14a in compound 3. Based on this data, we suggest that compound 3 has a 3R,17S configuration.

Sporochartine D (4) (Figure 4) conserved the correlations between H-17 and H-2 and between H-17 and H-14b reported for sporochartine A (1). Furthermore, the correlation between H-17 and H-13 is absent in both 4 and 1. In 4 we have an additional correlation between H-13 and H-2, absent in 1. These observations support the conclusion that 4 is the 3R,17R isomer of 1.

A new compound referred as sporochartine E (5) was also isolated as a white powder. Compound 5 has the same molecular formula C24H34O6 as compound 1, deduced from HRESIMS m/z [M+H]+ 419.2433. Here again we have eight degrees of unsaturation accounting for two γ-lactones, two double bonds, one six-membered cycle moiety and one tetrahydrofurane moiety.

Compound 5 has a terminal methylene group (at δC 119.6, δH 5.27 and δH 5.21) while the tetrahydrofuran moiety connected to C-19 in 1 is connected to C-14 in 5.

Based on COSY correlations (Figure 5), the sporothriolide moiety was the same in compound 5 as in 1. Moreover, COSY correlations from H-13 to H-19 through H-14 (δH 3.25), H-15 (δH 6.16), H-16 (δH 5.56), H-17 (δH 2.80) and H-18 (δH 5.65) together with HMBC correlation between H-13 and C-3 and H-17 and C-3 formed a cyclohexane fragment like in 1. Finally, by using the COSY correlations from H-24 (δH 1.34) to H-20 (δH 3.22), through H-23 (δH 3.96), H-22 (δH 1.53 and 2.39) and H-21 (δH 4.28) we deduce the tetrahydrofuran moiety. The HMBC correlations between H-20 and C-14 and C-15 allowed us to connect this tetrahydrofurane moiety to the sp3 methine C-14.

Figure 5

The absolute configuration of sporochartine E (5) was suggested using ROE correlations compared to the absolute configuration of sporochartine A (1) (Figure 6).

Figure 6

ROE correlations between H-2 and H-5/H-6 in the sporothriolide moiety and the ROE correlation between H-20 and H-21 and between H-23 and H-21 in the tetrahydrofuran moiety indicated a similar to that in 1.

Sporochartine E (5), showed a correlation between H-17 and H-2, like in compounds 1 and 4. H-2 also exhibited a correlation with H-13b but not with H-14. This suggests that C-3 and C-17 has the same relative configuration than 4. For C-14, we observed ROE correlations between H-14 and H-13a, H-13a, and H-17 and H-14 and H-24, suggesting a 3R, 14S, 17S configuration for compound 5.

Based on the structure of sporothriolide and the recently reported trienylfuranol A isolated from H. submoniticulosum, we suggested a hypothetic biosynthetic pathway of sporochartines, involving a “spiro” Diels-Alderase reaction as shown in Figure 7 (Klas et al., ; Byrne et al., ). The possibility of a non-enzymatic catalysis was excluded as reported previously (Leman-Loubière et al., ).

Figure 7

The cytotxicity of sporochartines was evaluated on three human cancer cell lines, HCT-116 (human colon carcinoma), PC-3 (prostate cancer cell lines) and MCF-7 (breast cancer cell line). The results presented in Table 3 are highly contrasting but nevertheless clearly indicate the impact of sporochartine conformation on the bioassay results.

Table 3

IC50in μM
HCT116PC3MCF7
Sporochartine A (1)>100>100>100
Sporochartine B (2)28.7 ± 2.5>1002.28 ± 0.1
Sporochartine C (3)7.2 ± 0.2113.4 ± 1.221.5 ± 0.3
Sporochartine D (4)>10015.2 ± 1.7>100
Sporochartine E (5)>1005.96 ± 0.28>100

IC50 values recorded for sporochartines A–E against three cancer cell lines, HCT-116 (human colon carcinoma), PC-3 (human prostate cancer cell lines), and MCF-7 (human breast cancer cell line).

Thus, sporochartine C (3) is toxic against the three cell lines with IC50 ranging from 7.2 to 21.5 μM. In contast, sporochartine A (1) is totally inactive at concentrations higher than 100 μM. This may be due to the substantial difference in the spatial conformation of compounds 1 and 3 (Figure 4).

The lower IC50 values were recorded for different sporochartines and against different cell lines, sporochartine B (2) for MCF-7 (2.28 μM), sporochartine C (3) for HCT-116 (7.2 μM) and sporochartine E (5) for PC-3 (5.96 μM).

Our future efforts will focus on the cytotoxic profile, biosynthesis and synthesis of sporochartines. The cytotoxicity profile reveals a non-cytotoxic sporochartine A (1), a large spectrum cytotoxic sporochartine C (3) and more cell line specific sporochartines B (2), D (4), and E (5). This finding merits future investigation on the mechanisms of action of these new scaffolds of cytotoxic compounds.

The biosynthesis of sporochartines, and the biosynthesis of its two moieties, sporothriolide and trienylfuranol A are still unknown. This opens new and promising opportunities for the discovery of novel biosynthetic microbial clusters.

Finally, having in hand hundreds of milligrams of sporothriolide, the hemi-synthesis of sporochartines is currently in progress based on a final Diels-Alder connection. The selectivity of the chemical catalysis and the proportion of different isomers will be compared to the microbial counterpart. According to our expertise in biocatalysis-based chemodiversification of natural compounds (Adelin et al., ; Martins et al., ), sporothriolide will be submitted to a panel of microorganisms in order to pursue the enrichment of sporothriolide related compounds.

Materials and methods

General experimental procedures

Optical rotations [α]D were measured using an Anton Paar MCP-300 polarimeter. IR spectra were obtained using a Perkin Elmer BX FT-IR spectrometer. NMR experiments were performed using a Bruker Avance 500 MHz in CDCl3 at room temperature. High-resolution mass spectra were obtained on a Waters LCT Premier XE spectrometer equipped with an ESI-TOF (electrospray-time of flight) by direct infusion of the purified compounds. Preparative HPLC was performed using Waters modules consisting of an autosampler 717, a pump 600, a photodiode array detector 2996 and an evaporative light-scattering detector, ELSD 2420. Prepacked silica gel Redisep columns were used for flash chromatography using a Combiflash-Companion chromatogram (Serlabo, France). All other chemicals and solvents were purchased from SDS (France).

Animal material

The Sphaerocladina sponge was collected on 17 December 2015 from the coast of Tahiti (9°45.421′S–139°08.275′W) at 20 m depth (Leman-Loubière et al., ).

Hypoxylon identification and cultivation

H. monticulosum CLL205 was isolated from the sponge Sphaerocladina and grown at 28°C on a PDB medium (Potatoes Dextrose Broth, DIFCO). The ITS rDNA gene amplification and sequencing were performed, and submitted to NCBI/BLAST database (GenBank). The primers used for PCR amplification were ITS1 F: CTT GGT CAT TTA GAG GAA GTA A (Tm: 55°C) and ITS4: TCC TCC GCT TAT TGA TATGC (Tm: 53°C). The GenBank accession number for H. monticulosum CLL205 sequence is SUB2477083 25758633.seq KY744359. H. monticulosum CLL205 was cultivated in a 2 L Erlenmeyer containing 1L of PDB medium (DIFCO) in a rotary shaker at 28°C and 130 rpm.

Compounds isolation

The culture broth was extracted with ethyl acetate (3 × 500 mL). The solvent was concentrated to dryness in vacuo to afford 430 mg of crude extract. 300 mg were submitted to flash chromatography on a Combiflash Companion using a Redisep 12 g silica column, eluted with a heptane-ethyl acetate mixture. After concentration in vacuo, we obtained sporothriolide (30 mg), compound 1 (9 mg), 2 (14 mg), 3 (4 mg), 4 (3 mg), 5 (1 mg).

Cytotoxicity assays

A tetrazolium dye [3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium-bromide; MTT]-based colorimetric assay was used to measure the inhibition on the proliferation of various human tumor cell lines HCT-116 (human colon carcinoma), PC-3 (prostate cancer cell lines) and MCF-7 (breast cancer cell line). The tested compounds were formulated in DMSO and added to the cells such that the final DMSO concentration ranged from 1 to 3%. Cells were grown in D-MEM medium supplemented with 10% fetal calf serum (Invitrogen), in the presence of penicillin, streptomycin, and fungizone, and plated in 96-well microplates. After 24 h of growth, cells were treated with target compounds from 100 μM to 10 nM. After 72 h, MTS reagent (Promega) was added, and the absorbance was monitored (490 nm) to measure the inhibition of cell proliferation compared to untreated cells. IC50 determination experiments were performed in separate duplicate experiments.

Isolated compounds

Sporochartine A (1)33: white needles, M.p. 86.5–87.9°C; −57 (c 0.5, CHCl3). See Tables 1, 2 for complete 1H, 13C NMR and IR data. HRESIMS m/z 419.2433 [M + H]+ (calcd for C24H35O6, 419.2433).

Sporochartine B (2): white powder; +72 (c 1.0, CHCl3). See Tables 1, 2 for complete 1H, 13C NMR and IR data. HRESIMS m/z 419.2419 [M + H]+ (calcd for C24H35O6, 419.2433).

Sporochartine C (3): white powder (4 mg); +93 (c 0.27, CHCl3). See Tables 1, 2 for complete 1H, 13C NMR and IR data. HRESIMS m/z 419.2433 [M + H]+ (calcd for C24H35O6, 419.2433).

Sporochartine D (4): white powder; −152 (c 0.27, CHCl3). See Tables 1, 2 for complete 1H, 13C NMR and IR data. HRESIMS m/z 419.2431 [M + H]+ (calcd for C24H35O6, 419.2433).

Sporochartine E (5): white powder; +51 (c 0.3, CHCl3). See Tables 1, 2 for complete 1H, 13C NMR and IR data. HRESIMS m/z 419.2425 [M + H]+ (calcd for C24H35O6, 419.2434).

Structural elucidation data are reported in the Supplementary Materials.

Associated content

Detailed 1D and 2DNMR, MS and IR spectra of sporochartines are available free of charge via the Internet at http://pubs.acs.org.

Statements

Author contributions

CL-L: microbiologie chemistry; GL: chemistry; CD: invetebrate investigation; JO: head of the team and science manager.

Acknowledgments

This work was supported by TASCMAR project funded by the European Union's Horizon 2020 research and innovation program under grant agreement n° 634674, and the Ph.D. support program of ICSN.

The authors are grateful to the CIBI Plateforme of CNRS-ICSN for cytotoxicity essays (Jérôme Bignon & Hélène Levaique).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2017.00399/full#supplementary-material

References

Summary

Keywords

Hypoxylon, Sphaerocladina, sporothriolide, sporochartines, cytotoxic compounds

Citation

Leman-Loubière C, Le Goff G, Debitus C and Ouazzani J (2017) Sporochartines A–E, A New Family of Natural Products from the Marine Fungus Hypoxylon monticulosum Isolated from a Sphaerocladina Sponge. Front. Mar. Sci. 4:399. doi: 10.3389/fmars.2017.00399

Received

17 October 2017

Accepted

24 November 2017

Published

12 December 2017

Volume

4 - 2017

Edited by

Antonio Trincone, Istituto di Chimica Biomolecolare (CNR), Italy

Reviewed by

Johannes F. Imhoff, GEOMAR Helmholtz Centre for Ocean Research Kiel (HZ), Germany; Usama Ramadan Abdelmohsen, University of Würzburg, Germany

Updates

Copyright

*Correspondence: Jamal Ouazzani

This article was submitted to Marine Biotechnology, a section of the journal Frontiers in Marine Science

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics