Brasilianoids A–F, New Meroterpenoids From the Sponge-Associated Fungus Penicillium brasilianum

3,5-Dimethylorsellinic acid (DMOA) derived meroterpenoids comprise an unique class of natural products with diverse scaffolds and with a broad spectrum of bioactivities. Bioinformatics analysis of the gene clusters in association with the qRT-PCR detection of the amplification of two key genes led to speculate that the sponge associated fungus Penicillium brasilianum WZXY-m122-9 is a potential producer of meroterpenoids. Chromatographic separation of the EtOAc extract of this fungal strain on a large-scale fermentation resulted in the isolation of six new DMOA-related meroterpenoids with trivial names of brasilianoids A–F (1-6), together with preaustinoid D and preaustinoid A2. The structures were determined by extensive analyses of spectroscopic data, including the X-ray diffraction and the ECD data for configurational assignment. Brasilianoids A and F showed an unprecedented skeleton with a γ-lactone in ring A, while brasilianoids B–C featured a 7/6/6/5/5 pentacyclic ring system finding in nature for the first time. The biosynthetic relationship among the isolated compounds was postulated. Compound 1 significantly stimulated the expression of filaggrin and caspase-14 in HaCaT cells in dose-dependent manner, while compounds 2 and 3 showed moderate inhibition against NO production in LPS-induced RAW 264.7 macrophages.


INTRODUCTION
The fungal meroterpenoids as the fascinating hybrid natural products are widely distributed in marine environments with diverse molecular architectures, that are assembled by terpene moieties with other precursors such as polyketide unit by various biosynthetic pathways (Iida et al., 2008;Silva et al., 2011;Guo et al., 2012;Qi et al., 2017). Among the fungus-derived meroterpenoids, a polyketide-terpenoid biosynthetic pathway that has C-alkylation of 3,5-dimethylorsellinic acid (DMOA) with farnesyl pyrophosphate (FPP) generated more than 100 secondary metabolites with a number of unique scaffolds (Geris and Simpson, 2009;Matsuda and Abe, 2016). The biogenetic pathways of these typical natural products have been extensively investigated, uncovering a series of synthetic gene clusters and functional enzymes Mori et al., 2017). The structure variation of the DMOA-based meroterpenoids was attributed to sequential cyclization, complex oxidative ring rearrangement, and recyclization. Based on the carbocyclic frameworks, the DMOA-FPP derived meroterpenoids can be classified into seven subtypes. Andrastins with a 6,6,6,5-tetra-carbocyclic skeleton are the potent inhibitors of RAS proteins, which are important for controlling cell division and the development of cancer (Nielsen et al., 2005). Terretonin-type congeners possessing a δ-lactone in ring D are originated from terrenoid (andrastintype) by D-ring expansion and unusual rearrangement of the methoxy group (Matsuda et al., 2015). Berkeleyonetype (or protoaustinoid-tyep) derivatives as the caspase-1 inhibitors are the meroterpenoids representing a series of unique and functionalized chemical scaffolds, which are characterized by the presence bicyclo[3.3,1]nonane or its rearranged bicyclo[3,2,1]octane unit in rings C and D. (Stierle et al., 2011), and are derived by the same intermediate as for andrastins with different rearrangement. Austinol and its analogues featured a pentacyclic scaffold with a spiro-δ-lactone in ring A and a γ -lactone in ring E, that were derived from protoaustinoid through oxidation and ring rearrangement (Matsuda et al., 2013). Chrysogenolides are a group of DMOAbased compounds with a unusual seven-numbered ring B, which exhibited the inhibition of nitric oxide production (Qi et al., 2017), while anditomin analogues featured the presence of a unique and highly oxygenated bridged-ring system (Matsuda et al., 2014). Fumigatonin and novofumigatonin are an additional subtype containing highly oxidized and complicatedly condensed ring system (Okuyama et al., 2008;Rank et al., 2008). These meroterpenoids have been reported to possess a range of biological activities (Rank et al., 2008;Zhang et al., 2012). Penicillium brasilianum is a remarkable microorganism with great potential to produce secondary metabolite with a variety of carbon skeletons and with interesting biological activities. Fungal strains of P. brasilianum are mostly isolated from terrestrial sources. They produced diverse metabolite scaffolds including meroterpenenes, polyketides, alkaloids, and cyclopeptides (Bazioli et al., 2017). In the course of our search for new terpenoids with pharmaceutical bioactivity, a sponge-associated fungal strain, P. brasilianum WZXY-m122-9, was selected due to the same fungal species has been reported to produce meroterpenoids with chemical diversity and novelty Bazioli et al., 2017). Anti-SMASH genome sequence analysis revealed that two gene clusters (clusters A and B) in WZXY-m122-9 are highly identical to the aus' clusters (Figure 1), that playing key role for the biosynthesis of DMOA-derived meroterpenoids in P. brasilianum MG11 . qRT-PCR detection showed that prenyltransferase gene (PM-122-9_1376) and terpene cyclase gene (PM-122-9_1374) in clusters A and B were highly expressed. These findings suggested WZXY-m122-9 strain to be a potential producer of terpenoids. Chromatographic separation of the EtOAc extract of the cultured fungus on a large scale resulted in the isolation of eight meroterpenoids, including six new compounds namely brasilianoids A-F (1-6) and known analogues preaustinoids D (7) and A2 (8) (Figure 2).

General Experimental Procedures
Optical rotations were measured with an Autopol III automatic polarimeter. IR spectra were measured with a Thermo Nicolet Nexus 470 FT-IR spectrometer. 1D and 2D NMR spectra were recorded on a Bruker Avance-400/600FT NMR spectrometer using TMS as internal standard. HRESIMS data were acquired on a Bruker APEX IV 70 eV FT-MS spectrometer. ESIMS data were obtained on a Finnigan MAT 95 mass spectrometer. The column chromatographic substrates included silica gel (200-300 mesh) and HF254 silica gel for TLC (Qingdao Marine Chemistry Co. Ltd.), Sephadex LH-20 (18-110 µm; Pharmacia Co., Ltd.); ODS (50 µm; YMC, Milford, MA). Semipreparative HPLC separation was performed with an Alltech instrument (426-HPLC pump) equipped with an UV detector at 210 nm and using a Prevail-C 18 column (Semipreparative, 5 µm). X-ray data were measured with a Bruker SMART APEX-II DUO instrument.

Fungal Strain and Identification
The fungus P. brasilianum WZXY-m122-9 was isolated from an unidentified marine sponge, which was collected in July 2016 from Weizhou Island in the South China Sea, and was identified by microscopic examination and 18S rDNA ITS sequence's BLAST in GenBank (GenBank accession number HM469396). A voucher specimen (WZXY-m122-9) was deposited at the State Key Laboratory of Natural and Biomimetic Drugs, Peking University, China.

Genome Sequencing and Analysis
Genome sequencing of P. brasilianum WZXY-m122-9 was performed by Sangon Biotech (Shanghai) Co., Ltd. (Shanghai, China) with an Illumina HiSeq 2000 system. Sequence assembly was performed with SPAdes version 3.5.0 (http:// cab.spbu.ru/software/spades/) to yield 1,367 contigs covering approximately 34.7 Mb. Gene prediction was then performed with Prokka (https://github.com/tseemann/prokka). Anti-SMASH (antibiotics and secondary metabolite analysis shell) analysis of genome sequences was performed to detect secondary metabolite gene clusters, and accurate gene cluster alignment was performed manually by comparisons with homologous genes found in the NCBI database. Anti-SMASH of WZXY-m122-9 genome analysis displayed that two gene clusters namely cluster A and cluster B, showed high identity compared with aus' clusters that had been reported in P. brasilianum MG11 .

Quantitative RT-PCR for Gene Expression
Degenerate primers were designed according to the conserved sequences of prenyltransferase (PM-122-9_1376 ′ ) and terpene cyclase (PM-122-9_1374) genes within clusters A and B, and the length of PCR products were 805 bp and 529 bp, respectively. The gene expression levels were detected by qRT-PCR. Total RNA of WZXY-m122-9 was obtained from the rice culture medium, FIGURE 1 | Comparison of the gene clusters between P. brasilianum WZXY-m122-9 and P. brasilianum MG11. (A) Comparison of the genes in cluster A of P. brasilianum WZXY-m122-9 with the aus's genes (cluster A) from P. brasilianum MG11. The homologous genes between the gene clusters were linked with lines. (B) Comparison of the genes in cluster B of P. brasilianum WZXY-m122-9 with the aus's genes (cluster B) from P. brasilianum MG11. The homologous genes between the gene clusters were linked with lines. PM-122-9-5653 and others mean the genes in P. brasilianum WZXY-m122-9. Frontiers in Chemistry | www.frontiersin.org and cDNA was synthesized from 1 µg of total RNA in a volume of 20 µL using TransScriptIIAll-in-One first-strand cDNA synthesis super mix (Transgene) according to the manufacturer's instructions. 0.4 µL cDNA, forward primer (10 µM), reverse primer (10 µM) and 10µL 2× TransStart Top Green qPCRSuperMix (Transgene) were used in subsequent RT-PCR reactions with supplement ddH 2 O to 20 µL. The specific primers included actin-F (5 ′ -ACCTGCTCTGCGACTACAAC-3 ′ ), actin- PM-122-9_1374-R (5 ′ -ACGGCGGAGTGTAGGAAGAA-3 ′ ). Each primer pair gave a single PCR product. The β-actin gene was treated as an internal reference gene, optimized PCR conditions were 94 • C for 30 s; 45 cycles of 94 • C for 5 s; 60 • C for 15 s; and 72 • C for 10 s; followed by dissociation stage. Three parallel tests for each reaction and recording their respective Ct values. According to amplification curves, every sample was different from no template control (NTC), which indicated that the genes were amplified from cDNA template. Under current culture conditions, the aus ′ gene cluster was expressed in the plasmids of WZXY-m122-9, and the expression level of prenyltransferase (PM-122-9_1376 ′ ) and terpene cyclase (PM-122-9_1374) genes was higher than β-actin due to Ct 1376 and Ct 1374 less than Ct actin .

Fermentation of the Fungus
The fermentation was carried out in 40 Fernbach flasks (500 mL), each containing 80 g of rice. Distilled H 2 O (100 mL) was added to each flask, and the contents were soaked overnight before autoclaving at 15 psi for 30 min. After cooling to room temperature, each flask was inoculated with 5.0 mL of the spore inoculum (10 7 /mL) and incubated at 25 • C for 25 days.

Assay for in Vitro Anti-HBV Effects
HepAD38 cells were cultured with DMEM in 48-well plate at 1 × 10 5 cells/well for 24 h, then treated with 10 µM of each compound for 72 h. After 4 days, the cells were washed with precooled PBS for 2 times. HBV progeny DNA of HepAD38 cells were extracted using QIAamp DNA Blood Mini kit (Biomiga) according to the manufacturer's instruction. Total DNA was reversely transcribed using PrimeScript RT reagent Kit (Takara, Dalian, China). The primers were designed and synthesized by Takara, and the sequences of the primers are indicated in Table  S9. PCR amplification was performed on an StepOne Plus real time PCR system (Applied Biosystems, Foster City, CA) using the SYBR Green Master Mix (Applied Biosystems, Foster City, CA). All experiments were performed in triplicate, and the relative levels of assayed HBV DNA were calculated with the deltadelta CT method using lamivudine as a positive control, and normalized to non-treated control.

Assay for Protective Effects on Skin Barrier Functions in in Vitro
HaCaT cells were cultured with DMEM in 6-well plate at 1 × 10 5 cells/well for 24 h, then treated with 20, 10, 5 µM of compounds 1-7 for 72 h. Total RNA of HaCaT cells were extracted using RNA Miniprep kit (Biomiga) according to the manufacturer's instruction. Then total RNA was reverse transcribed using PrimeScript RT reagent Kit (Takara, Dalian, China). The primers were designed and synthesized by Takara, the sequences of the primers are indicated in Table S10. PCR amplification was performed on an StepOne Plus real time PCR system (Applied Biosystems, Foster City, CA) using the SYBR Green Master Mix (Applied Biosystems, Foster City, CA). All experiments were performed in triplicate, and the relative levels of assayed mRNAs were calculated with the delta-delta CT method using ACTIN expressions as endogenous control, and normalized to non-treated control.
HaCaT cells were seeded at 5 × 10 5 cells per well in 6-well plates and incubated under starvation conditions for 24 h using serum-free MEM. HaCaT cells were exposed to UVB irradiation (UVB lamp: SUV-1000, Sigma-China) at a dose of 30 mJ/cm 2 . After UVB irradiation, DMEM medium containing 1 (20 µM) or epigallocatechin gallate (20 µM) was added and the cells were incubated for 24 h, and then the MTT assay was performed.
Assay for the Inhibition of Nitric Oxide Production in RAW264.7 Macrophages All compounds were assayed by measuring the inhibitory effects on NO production induced by LPS in mouse macrophage RAW 264.7 cells. The RAW 264.7 cells were seeded in each well of a 96-well plate at a concentration of 5 × 10 5 cells/ml in DMEM medium containing 10% fetal bovine serum, 2 mmol/l glutamine, 100 U/ml penicillin, and 100 µg/ml streptomycin. After incubation at 37 • C in a humidified atmosphere of 5% CO 2 air for 12 h, cells were treated with LPS (1 µg/ml) and compounds (50 µg/ml for each compound, and compound was dissolved in DMSO to dilute until the final DMSO concentration < 0.1%, v/v), then incubation for 24 h at 37 • C. NO production was determined by measuring the quantity of nitrite in the culture medium by the Griess reagent. Aminoguanidine was selected as a positive control. A 50 µl of the supernatant of culture medium was mixed with 100 µl of the Griess reagent. Using nitrite to generate a standard curve, nitrite production was measured by a microplate reader at 540 nm. Cell viability was examined by the MTT method (Wu et al., 2016). The compounds with inhibitory rate more than 50% at a dose of 50 µg/ml were further detected in the gradient concentrations to calculate IC 50 values.

Statistical Analysis
All quantitative values are given as means ± SEMs. Significant differences among the experimental groups were assessed by one-way ANOVA. Statistical significances were considered at the * P < 0.05 and * * P < 0.01 levels.

Crystal Structure Analysis
X-ray crystal data of 2-5 and 8 were acquired on an Eos CCD with a graphite monochromated Cu Kα radiation (λ = 1.5418 Å) (Supporting Information). These structures was solved with ShelXT-97 using direct methods and refined with ShelXL using full matrix least-squares on F 2 . Crystallographic data for 2-5 and 8

Identification of New Compounds
The molecular formula of brasilianoid A (1) Table 2), which were classified into six methyl, five methylene (one olefinic carbon), four methine, and eleven quaternary carbons involving two ketones and two ester carbonyl carbons. Analyses of the 1D and 2D NMR data revealed that 1 is a preaustinoid-related meroterpenoid, while the partial structure regarding rings C/D/E of 1 was corresponded to rings B/C/D of preaustinoid A2 (8)  , ascertained the location of a hydroxy group at C-1 (δ C 70.0). A methylene bridge across C-2 (δ C 48.0) and C-4 (δ C 86.8) was evident from the HMBC correlations from H 3 -14 (δ H 1.29, s) to C-4, C-5 (δ C 48.3), and C-15 (δ C 40.7), and from H 2 -15 to C-1, C-2, C-4, and C-5. Additional HMBC correlation from the carbonyl carbon C-3 (δ C 177.5) to H-1, H-2, and H 2 -15 in association with the degrees of molecular unsaturation confirmed the formation of a γ-lactone for ring A.
The relative configuration of 1 was established by the NOE data. As shown in Figure 3, the NOE correlations from H-15a (δ H 2.24) to OH-1 and H-5, and from H-9 (δ H 1.50) to OH-1, H-5, and H 3 -9 ′ , in association with the correlation between H 3 -12 (δ H 1.17) and H 3 -13 (δ H 0.76), clarified trans fusion between B/C and C/D rings. When H-1 was arbitrarily assigned to αorientation, the methylene bridge H 2 -15 was ascertained to be in β-orientation. The relative configuration in rings C to E was the same as that of preaustinoid A2, based on the similar NOE relationship. The spatial closure between H 3 -9 ′ and H-9 was recognized by their NOE interaction, while the NOE correlation between H 3 -12 and H 2 -7 ′ indicated that the exo-vinyl group across C-1 ′ and C-3 ′ was α-oriented. In biogenetic consideration, compound 1 was assumed to be derived from preaustinoid A2 by the formation of C-2/C-15 bond. Thus, both compounds were suggested to share the same absolute configuration in rings C-E. Since the absolute configuration of preaustinoid A2 was determined by the single crystal X-ray diffraction analysis using Cu K α radiation in present work, the stereogenic centers in rings C-E of 1 were supposed to be 8S, 9S, 10R, 3 ′ R, 5 ′ S, and 7 ′ R configurations (Supporting Information). Thus, the absolute configuration in rings A and B was in agreement with 1R, 2R, 4S, and 5S, based on the NOE relationships. This assignment was further supported by the ECD data (Stephens and Harada, 2010), of which the Cotton effects of experimental and calculated ECD  for the model molecule of 1 with 1R, 2R, 4S, 5S, 8S, 9S, 10R, 3 ′ R, 5 ′ S, and 7 ′ R were comparable (Figure 4). Brasilianoid B (2) was isolated as colorless monoclinic crystals. Its molecular formula was determined as C 25 H 32 O 6 based on the HRESIMS (m/z 473.2170 [M + HCOO] − ) and NMR data, containing ten degrees of unsaturation. The 1 H and 13 C NMR data of 2 (Tables 1, 2) regarding rings A-C resembled those of preaustinoid A2, while the distinction was attributed to rings D and E. An exocyclic olefinic group at C-2 ′ (δ C 147.6) in ring C was recognized by the HMBC correlations from the olefinic methylene protons H 2 -7 ′ (δ H 4.86, 5.09) to C-1 ′ (δ C 66.5) and C-3 ′ (δ C 55.3). Additional HMBC correlations from H 3 -8 ′ (δ H 1.20, s) to C-2 ′ , C-3 ′ , and a ketone carbon C-4 ′ (δ C 213.4), and from OH-6 ′ (δ H 6.95, s) to C-4 ′ , C-1 ′ , and C-6 ′ (δ C 90.8) established a cyclopentanone unit for ring D. Moreover, the HMBC correlations from H-5 ′ (δ H 4.30, q, J = 6.4 Hz) to C-1 ′ , C-4 ′ , C-6 ′ , and a carbonyl carbon C-10 ′ (δ C 172.6), in addition to the COSY coupling between H-5 ′ and H 3 -9 ′ (δ H 1.14, d, J = 6.4 Hz), deduced a γ-lactone unit for ring E, while a hydroxy group and a methyl group were located at C-6 ′ and C-5 ′ , respectively.
The relative configuration of 2 in rings A-C was the same as that of preaustinoid A2, based on the similar NOE relationships. The NOE correlations from OH-6 ′ to H-7a (δ H 2.23) and H-9 (δ H 1.99) demonstrated the spatial closure of these protons,  C 25 H 31 O 7 , 443.2070) and NMR data, requiring ten degrees of unsaturation. Analyses of 1D and 2D NMR data indicated that rings A to C are in accordance with those of 1, as evident from the proton spin systems from d,J = 4.5 Hz) to H 2 -15 (δ H 1.87, 2.54) and from H-5 (δ H 1.67) to H 2 -7 (δ H 1.63, 2.16) in the COSY spectrum, in addition to the key HMBC correlations observed from H 2 -15 to C-1, C-2, C-3, C-4 and C-5; H 3 -14 to C-4, C-5 and C-15; H 3 -12 to C-7, C-8, C-9, and C-1 ′ ; H 3 -13 to C-1, C-5, C-9, and C-10, and OH-1 to C-1, C-2 and C-10. Moreover, the HMBC and COSY correlations assigned the partial structure of rings D to E to be the same as that of 2. The similar NOE data, such as the key NOE correlations from OH-1 to H-15b (δ H 2.54), H-5 and H-9 and between H 3 -12 and H 3 -13 as observed in the NOESY spectrum of 1 suggested the same relative configuration in rings A-C of both 6 and 1. The NOE correlations from H-5 to H-15b and H-7b (δ H 2.16), H-7b to H-9 and OH-6 ′ , and between OH-6 ′ and H 3 -9 ′ (Figure 2) clarified these protons to be cofacial, indicating the same relative configuration in rings C-E of both 6 and 2. Additional NOE correlations from H-5 ′ to H 2 -7 ′ and H 3 -8 ′ supported the exo-vinyl bridge C-2 ′ to be in the same orientation as H-5 ′ . In order to assign the absolute configuration, the ECD data of 6 and its enantiomer were calculated at the B31YP/6-311++G(2d, 2p) level in the gas phase using the B3LYP/6-31G(d) optimized geometries after conformational searches via the MMFF94S force field (Ding et al., 2007). Comparison of the experimental ECD data with those calculated for the enantiomers indicated that 6 (Figure 4) led to the identification of 1R, 2S, 4R, 5S, 8S, 9S, 10R, 3 ′ R, 5 ′ R, 6 ′ R, and 7 ′ R configurations.
Analyses of 1D and 2D NMR data established the planar structure of 7 to be identical to preaustinoid D (Duan et al., 2016), while the similar NOE correlation and ECD data clarified both compounds to be identical. However, comparison of the experimental ECD data with that calculated for a model molecule of (5S, 8S, 9S, 10S, 3 ′ R, 5 ′ S, 6 ′ R, 7 ′ R)-7 (Figure 4) suggested that the structure of preaustinoid D should be revised to its enantiomer. Compound 8 was identical to preaustinoid A2 based on the comparison of its spectroscopic data with these reported in the literature (dos Santos and Rodrigues-Fo, 2003), while the absolute configuration was determined by the Flack parameters of the X-ray diffraction of the single crystals.

Biological Activity
Compounds 1-8 showed no cytotoxic activities toward tumor cell lines A549, U937 and SMMC7721 with IC 50 > 10 µM. However, compound 1 significantly stimulated the expression of filaggrin and caspase-14 in HaCaT cells with a dose dependent manner (Figure 6). Filaggrin is a key natural moisturizing factor that maintains the ability to regulate the skin moisture barrier (Rawlings and Harding, 2004), while the mutation of the filaggrin gene resulted in skin losing water to allow the entrance of bacteria, leading to allergies, irritation and infection (Nomura et al., 2007). Caspase-14 deficiency was associated with the accumulation of incompletely degraded filaggrin fragments within the stratum corneum, decreased stratum corneum hydration, increased transepidermal water loss, and sensitivity to UVB photodamage (Eckhart and Tschachler, 2011). To test the skin protective activity against UVB irradiation, the cytotoxicity of 1 against HaCaT cells was measured by the MTT assay. The cell viability was decreased to 70% compared to the normal group after exposure to UVB 30 mJ/cm 2 . Treatment of the damaged cells with compound 1 (20 µM) resulted in the cell viability increasing to 77%, while the positive control epigallocatechin gallate increased the cell viability to 75% at the same dose. This finding indicated that compound 1 is able reduce the UVB-induced cell damage.
Compounds 2 and 3 showed moderate inhibition toward NO production in LPS-induced RAW 264.7 macrophages ( Table 3). In addition, compounds 3-5 (10 µM) inhibited the DNA expression of the HBV virus in HepG2.2.15 cells with the inhibitory rates of 25, 15, and 10%, respectively, whereas lamivudine as a positive control exerted the inhibitory rate of 75% in the same dose (Supporting Information, Figure S8).

DISCUSSION
In present work, we reported P. brasilianum WZXY-m122-9 to be obtained from a marine sponge for the first time. It provides a new source to produce austin-related cyclization, Baeyer-Villiger oxidation to yield preaustinoid A1 (a product coexisted with 1-6 in the fungal strain) (Itoh et al., 2012). The latter compound can be converted to preaustinoid A2 by dehydrogenation. Compounds 1 and 6 presented as unique scaffold with γ-lactone unit in ring A. Both compounds were biogenetically presumed to be converted from preaustinoid A2 by the radical-based pathway, similar to the case for the biotransformation of preandiloid C to andiconin (Matsuda et al., 2014). During the formation of γ-lactone, Me-15 in preaustinoid A2 was catalyzed by a Fe(II)/α-KG-dependent dioxygenase to yield a radical intermediate, which subsequently formed a C-15/C-2 bond with the electron-rich olefin at C-2 to generate an intermediate with a free radical at C-1. The latter carbon was hydroxylated to yield compound 1. The formation of γlactone in 6 was suggested to follow the similar manner as for 1. Conversion of preaustinoid A2 to 2 and 3 was presumably mediated by the enzymes as the case of the biotransformation from preaustinoid A3 to austinol. Coexistence of C-5 ′ epimers (2, 3) with equal amount suggested that the stereogenic selection for the reduction of ketone at C-5 ′ to hydroxy group is unspecific. Compounds 4 and 5 are regarded to be derived from preaustinoid A1, which followed the similar rearrangement as that for the conversion of preaustinoid A2 to 2 and 3. We isolated a dioxygenase from the fungus, which was identical to AusE (Matsuda et al., 2013;Valiante et al., 2017). The dioxygenase can convert preaustinoid A2 to preaustinoid A3 in vitro experiment. However, The transformation of preaustinoid A2 to 1 or 6 was unsuccessful. This data suggested that the other Fe(II)/α-KG-dependent dioxygenase is required for the generation of the methylene bridge in 1 and 6. The methyl ester of 4 and 5 is derived by the hydrolysis of the lactone in ring-A to yield a carboxylic acid at the side chain, and then methyl esterification (Scheme 1). Although we found 4 and 5 in the HPLC chromatographic profile of EtOAc extract, whether both compounds to be derived by enzymatic step are uncertain.

CONCLUSION
In conclusion, compounds 1 and 6 possess unprecedented skeleton featuring the formation of C-2/C-15/C-4 bridge in ring A, while 2 and 3 with 7/6/6/5/5 pentacyclic skeletons of meroterpenoids are found from nature for the first time. The bioassay results indicated that the structural variation directly induced the bioactive spectrum. Compound 1 is the first example of natural products used to promote filaggrin and caspase-14 SCHEME 1 | Proposal of putative biosynthetic pathways for the compounds.
Frontiers in Chemistry | www.frontiersin.org expression for the protection of the UVB-induced cell damage, suggesting that it is a promising lead for the treatment of dermatological diseases.

AUTHOR CONTRIBUTIONS
JZ responses for the isolation of compounds. BY performed the biogenetic approach for synthetic gene screening. DL took the antiviral bioassay. SG assayed for the regulation of filaggrin and caspase-14. PP helped to revise the manuscript. WL elucidated the structures and edited the manuscript.

FUNDING
Financial supports from the MOST-973 program (2015CB755906) and the National Natural Science Foundation of China (81630089, 41376127, U1606403) are highly appreciated.