A widespread metabolic gene cluster family in metazoans

A widespread metabolic gene cluster family in metazoans | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Biological Sciences - Article A widespread metabolic gene cluster family in metazoans Bradley Moore, Natalie Grayson, Paul Scesa, Malia Moore, Jean-Baptiste Ledoux, and 5 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-4859447/v1 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 13 Jun, 2025 Read the published version in Nature Chemical Biology → Version 1 posted You are reading this latest preprint version Abstract Octocorals are unique among metazoans in their prolific production of bioactive terpenoid natural products that rival the chemical diversity of plants and microbes. We recently established that these cnidarians uniformly express terpene cyclases and that their encoding genes often reside within putative biosynthetic gene clusters (BGCs), a feature uncommon in animal genomes. In this work, we report the discovery and characterization of a widespread gene cluster family for the biosynthesis of briarane diterpenoids that number over 700 molecules specific to the Scleralcyonaceans, one of the two octocoral orders. We sequenced five genomes from evolutionarily distinct families of briarane-producing octocorals to complement three publicly available briarane-producing coral genomes, enabling the discovery of a conserved five-gene cluster composed of a terpene cyclase, three cytochrome P450s, and a short-chain dehydrogenase. Using Escherichia coli and Saccharomyces cerevisiae as hosts and homologous briarane biosynthesis genes from eight corals, we reconstituted the biosynthesis of cembrene B γ-lactone, which contains the γ-lactone structural feature distinctive of briarane diterpenoids. The discovery of the genomic basis of briarane biosynthesis establishes that animals, like microbes and plants, employ gene cluster families to produce specialized metabolites. Further, the presence of BGCs in octocoral proves that the formation and maintenance of BGCs over evolutionary time is a more widespread phenomenon in specialized metabolite biosynthesis than previously realized. Biological sciences/Biochemistry/Enzymes/Multienzyme complexes Earth and environmental sciences/Ocean sciences/Marine biology Biological sciences/Biotechnology/Metabolic engineering Biological sciences/Biological techniques/Sequencing/DNA sequencing Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Introduction Specialized metabolites, or natural products, are small molecules that while not vital for primary growth, enable a producing organism to survive in an ecological niche by mediating time-, tissue-, or stimulus-dependent functions. 1 Although most of these compounds are known from plants and microbes, they are widespread in all kingdoms of life, including animals, where they too serve diverse purposes including communication, development, and defense. 2 Examples include volatile sex-pheromones from frogs 3 and spiders, 4 polyketide pigments in parrots 5 and sea urchins, 6 sugar-derived hormones in nematodes, 7 hydroquinones employed by rove beetles to fend off predators, 8 and peptide toxins used by cone snails to stun their prey. 9 The overwhelming majority of research on specialized metabolite biosynthesis has been conducted in bacteria and fungi, where genes coding for a biosynthetic pathway are typically arranged in physically co-localized, biosynthetic gene clusters (BGCs). 10,11 In contrast, the few characterized biosynthetic pathways in animals are not encoded in BGCs, but have their encoding genes dispersed throughout the genome. This is currently regarded as the general rule for biosynthetic pathways in non-fungal eukaryotes. 12 However, improved genomic resources in other eukaryotes like plants and algae has revealed the colocalization of non-paralogous genes into BGCs, where they code for specialized metabolic pathways (Fig. 1 ). 1,13,14 Due to their abundance, eukaryotic BGCs were first recognized in fungi where they are now regarded as the norm in specialized metabolite biosynthesis. A prominent example is the lovastatin ( 1 ) BGC, which codes for nine non-homologous genes within a compact genomic region of 40 kilobase pairs (kb) and encodes the complete set of enzymes needed to produce 1 , a resistance gene, a transporter, and a transcription factor (Fig. 1 ). 15 The more recently discovered BGCs in plant genomes show clear differences from fungal BGCs. As generally the case in plant genomes, they show larger intergenic regions and often do not code for the complete biosynthetic pathway. An example is the tomatine ( 2 ) BGC, which codes for six steps in the pathway starting from cholesterol. 16 However, the first enzymatic step of the pathway is encoded by a gene residing over 7,000 kb away from this cluster, while two additional genes are coded as a gene pair on a different chromosome (Fig. 1 ). 17 There are only a few examples of small molecule BGCs in animals that have since been identified. However, their origins link to microbial gene sets acquired via horizontal gene transfer (HGT) as in the case of mycosporine-like amino acids such as gadusol ( 3 ) that widely occur in vertebrates 18 (Fig. 1 ) and a penicillin biosynthesis gene pair from the springtail Folsomia candida that shows high homology to fungal sequences. 19,20 In contrast, the genes involved in the biosynthesis of the hybrid polyketide/ non-ribosomal peptide hormone nemamide ( 4 ) from the nematode Caenorhabditis elegans do not show clear sequence homology with microbial genes and thus may have evolved de novo in nematodes. 21 However, none of these genes are colocalized in a cluster, but dispersed over large genomic regions across four separate chromosomes. Although uncommon, a few gene clusters in animal genomes that are not the result of horizontal gene transfer are known, such as the Hox 22 and b-globin 23 clusters in humans, as well as a laccase gene cluster recently discovered in rove beetles. 24 However, these gene clusters are made up of paralogous genes and have arisen from successive tandem duplications from one ancestor gene, and only the last example is involved in small molecule biosynthesis. Generally, the few characterized biosynthetic pathways in animals, including the extraordinarily well investigated examples of steroid hormones 25 or melatonin ( 5 ) 26 in humans, are not encoded in BGCs, but have their encoding genes dispersed throughout the genome. It was therefore surprising when we recently identified potential terpenoid BGCs in octocorals, consisting of up to 13 co-localized genes spanning up to 500 kb that did not result from HGT or pure tandem duplications. 27,28 Instead, over evolutionary time, coral gene paralogs from different protein families appear to have been recruited to single genomic scaffolds, forming potential BGCs. We observed related putative gene clusters from Renilla muelleri and Erythropodium caribeorum , corals that belong to different families (Pennatuloidea vs. Erythropodiidae) and differ in their lifestyle and morphology. Yet, both corals are known sources of briarane-type diterpenoids, a Scleralcyonacea-specific structural family of diterpenoids, numbering over 700 representatives and showing a wide array of potent biological properties. 29 Here, we test the hypothesis that briarane diterpenoids are encoded by broadly distributed and ancient biosynthetic gene clusters, which is unprecedented across metazoans. We assembled and analyzed high-quality genome sequences from six key briarane producing octocorals from different taxonomic families, revealing a conserved gene cluster family. Using a combination of protein biochemistry and pathway engineering, we reconstituted the biosynthesis of cembrene B γ-lactone, which contains the γ-lactone structural feature distinctive of all briarane diterpenoids. These results provide compelling evidence that corals use BGCs to synthesize their abundant bioactive terpenoids, opening the possibility that other metazoans may have similarly evolved BGCs in their own specialized metabolism. Results Chemical Validation of Briarane Production We previously characterized homologous cembrene B synthases identified from the genomes of R. muelleri and E. caribaeorum. 27,28 Both encoding genes were colocalized with cytochromes P450 genes, which also shared homology with their counterparts on the two distinct genomes. The tantalizing possibility that coral animals evolved gene cluster families for the biosynthesis of structurally related specialized metabolites motivated us to evaluate whether diverse corals known to contain briarane diterpenoids also encode cembrene B synthases in their genomes. As such, we first set out to confirm literature reports of briarane producing octocoral families. Previously, Pennatuloidea, Ellisellidae, Coralliidae, Briareidae, and Erythropodiidae were established as briarane containing octocorals (Fig. 2 a). 30–36 As each of these five coral families belong to the Scleralcyonacea order, we further examined the remaining 14 families within this taxonomic order and found that five others reported terpenoid chemistry, yet not with briarane scaffolds (Fig. 2 b, SI Table 1 ) . 37 For the remaining nine families, no terpenoid literature was found. We thus collected coral tissue from representative species across five known briarane producing families, including Renilla koellikeri (family Pennatuloidea), Stylatula elongata (family Pennatuloidea), Dichotella gemmacea (family Ellisellidae), the branching and encrusting morphologies of Briareum asbestinum (family Briareidae), E. caribaeorum (family Erythropodiidae), and Corallium rubrum (family Coralliidae). R. koellikeri collected, by SCUBA in San Diego, is not a reported source for briaranes, but is related to the Atlantic R. muelleri that produces briaranes 38 and whose genome encodes cembrene B synthase as part of a putative BGC. 27,28 We thus analyzed the tissue of R. koellikeri and isolated 11-hydroxyptilosarcenone ( 6 ) whose chemical structure was confirmed by NMR spectroscopy ( NMR supplementary note, SI Table 2, SI Fig. 1 ). This compound was previously identified from the sea pen Ptilosarcus gurneyi , but not yet from a Renilla sea pansy. 39,40 Limitations in accessible biomass for the other octocorals targeted in this study led us to use molecular networking analyses of crude extracts and semi-pure fractions to confirm briarane production. We used 11-hydroxyptilosarcenone ( 6 ) as a standard in the molecular network to identify a cluster containing 54 nodes representing unique masses of briarane compounds (Fig. 2 , SI Fig. 2 , SI Table 3 ). This analysis not only expanded the ptilosarcenone network from R. koellikeri (Pennatuloidea) but also established related briaranes in D. gemmacea (Ellisellidae), E. caribeorum (Erythropodiidae), and B. asbestinum (Briareidae). Two nodes corresponding to compounds in these corals had exact masses consistent with literature values: gemmacolide R ( 7 ) from D. gemmacea and erythrolide D ( 8 ) from E. caribaeorum (Fig. 2 c, SI Table 1 ). 41,42 However, we did not detect briarane nodes from extracts of C. rubrum (Coralliidae) or S. elongata (Pennatuloidea). While there is literature precedent of S. elongata and corals in the family Coralliidae, although not specifically C. rubrum , making briaranes, we included these corals in our study assuming that they should have the genomic capacity to produce briaranes even though the samples we analyzed were devoid at the time of collection. 33,34,43 Genome assembly and syntenic analysis of the briarane gene cluster family in Scleralcyonacea corals At the start of this study, a few briarane containing coral draft genomes with very low continuity (contig N50) were available: E. caribaeorum (Contig N50 2 kb) 28 , R. muelleri (N50 70.5 kb) 44 , and Pteroeides caledonicum (N50 4 kb) ( SI Table 4 for assembly accessions and statistics).To capture intact BGCs, we proceeded to generate higher quality assemblies of corals across four briarane producing families, obtaining long read sequencing data for R. koellikeri (family Pennatuloidea), S. elongata (family Pennatuloidea), D. gemmacea (family Ellisellidae), the branching and encrusting morphologies of B. asbestinum (family Briareidae), and E. caribaeorum (family Erythropodiidae). We assembled 26 to 78 gigabase pairs (Gb) of sequencing reads into genomes ranging in size from 231 megabase pairs (Mb) ( R. koellikeri ) to 1,237 Mb ( B. asbestinum ) ( Table 1 ). To assess the completeness of our assemblies, we estimated genome sizes with K-mer based analyses using Illumina short reads ( SI Fig. 3 ), producing genomes ranging from 144 Mb ( R. koellikeri ) to 1,155 Mb ( B. asbestinum ) ( Table 1 ). Illumina short read coverage was too shallow for E. caribaeorum to allow for K-mer based analysis ( SI Table 5 for all sequencing data used). Additionally, we included our recently reported chromosomally resolved C. rubrum (family Coralliidae) genome at 475 Mb (20 chromosomes, 951 scaffolds) 45,46 in our analysis to enable BGC exploration across all five known briarane producing coral families. The genome assemblies are now some of the most contiguous for octocorals with contig N50 values above 196 kb ( Table 1 ,Fig. 3 a, SI Table 6 ). The chromosomal resolution assembly of C. rubrum had a scaffold N50 of 16,290 kb, on par with the current best scaffolded octocoral assembly from Xenia sp. with a scaffold N50 of 14,832 kb. 47 Benchmarking universal single-copy ortholog (BUSCO) completeness scores ranged from 86.4–88% ( Table 1 ). BUSCO scores are highly dependent on the database used for the analysis and the low scores could reflect the use of the metazoan database, which was the closest available. The lower-than-expected completeness is on par with the BUSCO for the Xenia assembly, which was 88.1% complete (Fig. 3 a). 47 We further assembled transcriptomes to guide our genomic analyses by identifying actively transcribed genes and determining intron bounds ( SI Table 6 ). Using the putative briarane BGC gene sequences from R. muelleri and E. caribeorum as queries, we analyzed the five newly established genomes of D. gemmacea, B. asbestinum, R. koellikeri, C. rubrum , and S. elongata and the public P. caledonicum for homologous genes. In total, the genomes from these eight octocoral species span five of the 20 octocoral families in the order Scleralcyonacea. To determine if non-briarane producing octocoral families also have syntenic BGCs, we also evaluated the publicly available Heliopora coerulea (family Helioporidae) genome. 48 We identified homologs of the cembrene B terpene cyclase (cbTC) in all eight octocorals shown or suspected to produce briaranes, whereas the H. coerulea genome did not contain a homologous sequence, although it did encode nine other TCs. In the eight octocoral species containing a homologous cbTC gene, seven contained iterations of the briarane BGC (Fig. 3 b). Given that the P. caledonicum genome had poor contiguity, we were not surprised to find the cbTC homolog on a 3 kb contig and expect that it too is part of a BGC. We observed one conserved short-chain dehydrogenase (cbSDH) and up to three cytochrome P450s (cbCYPa, cbCYPb, and cbCYPc) per briarane BGC (Fig. 3 b, SI Fig. 4 ). In the eight octocorals with briarane BGCs, four maintained the full five-gene briarane BGC spanning 22 kb in R. muelleri up to 40 kb in C. rubrum , albeit with different gene organization. Notably, in the cases of S. elongata, E. caribeorum , and B. asbestinum , the orthologous genes were split between two BGC loci ( SI Table 7 ). For S. elongata and B. asbestinum , their respective cbTC and cbCYPc genes reside on one contig, while cbCYPa, cbCYPb, and at least two cbSDHs are on a second contig. In the case of E. caribeorum , cbCYPc resides outside the BGC contig. Also, in the case of B. asbestinum , the orthologous genes have much longer gaps resulting in a significantly larger BGCs spanning ~ 150 kb. S. elongata and B. asbestinum further have multiple copies of the SDH gene at three (85%, 89% and 90% sequence similarity) and two (46% sequence similarity) copies, respectively ( SI Fig. 5 ). The duplications of the SDHs and gene rearrangements suggest lineage specific independent evolution across octocoral species. Since we identified the same set of homologous cytochromes P450s in a wide range of octocoral species, we conducted a phylogenetic analysis of these sequences in the background of broader cnidarian cytochromes P450s to shed light on their evolutionary origin. Our analysis included CYPs from four non-octocoral cnidarians as a reference set: Hydra vulgaris (24 CYPs), Acropora digitifera (24 CYPs), Aurelia aurita (37 CYPs), and Nematostella vectensis (70 CYPs). 49 From octocorals we included CYPs from six briarane producing species D. gemmacea, B. asbestinum, R. koellikeri, S. elongata, C. rubrum , and E. caribaeorum , as well as the non-briarane producers Xenia sp. and H. coerulea , yielding 318 octocoral CYP sequences at ~ 40 per genome ( SI Table 8 ). We combined the annotated reference cnidarian CYPs with the annotated octocoral CYPs to build a phylogenetic tree (Fig. 3 c, SI Fig. 6 ). The CYPs from the four reference cnidarians were previously associated with accepted CYP clans and families, 49 thereby enabling us to examine the evolutionary context of the octocoral BGC CYPs. Notably, some octocoral CYPs fell into the clades associated with the cnidarian CYP clans, suggesting that they fulfill primary, widespread functions generally present in cnidarians. Others, including the BGC-associated CYPs, formed new monophyletic clades that were octocoral specific. The three briarane CYPs, cbCYPa-c, form monophyletic clades consisting of cbCYPa, cbCYPb, and cbCYPc homologs from different species (Fig. 3 c). While cbCYPb and cbCYPc share a common ancestor that may suggest related enzymatic functions, cbCYPa is distantly related and instead lies closer to clan 46 with predicted functions associated with steroid biosynthesis (Fig. 3 c). Interestingly, H. coerulea , while also from the order Scleralcyonacea has CYPs closely related to the cbCYPb and cbCYPc clades, but these genes are not co-localized hinting at the loss of the briarane BGC in this lineage. To test whether our phylogenetic approach could be used for genome mining in octocorals with lower quality assemblies, we next included CYP sequences from the poor contiguity genome assembly of P. caledonicum. We found corresponding sequences present in two of the three clades, thereby identifying cbCYPa and cbCYPb orthologs in the genome (Fig. 3 d). Table 1. Assembly statistics and BUSCO assessment of five octocoral genomes Organism Contigs (#) Assembly size (Mb) K-mer Based Genome size (Mb) Contig N50 (kb) Longest (kb) BUSCO (%C/ %D) Briareum asbestinum 22,526 1,254 1,155 332 4,126 87.3/4.6 Dichotella gemmacea 8,171 550 401 196 2,036 87.5/15.6 Renilla koellikeri 557 231 144 4,895 22,942 86.9/2.6 Stylatula elongata 392 360 269 5,139 39,131 88.0/1.0 Erythropodium caribeorum 8,289 300 NA 233 1,576 86.4/0.6 Corallium rubrum 45 20 chromosomes, 951 scaffolds 545 540 1,600 (scaffold N50 18,521kb) 115,000 88.5/1.2 Biochemical validation of the briarane BGC To provide support that the briarane BGC is involved in the early stages of briarane biosynthesis, we set out to provide biochemical validation of the clustered biosynthesis genes ( SI Table 9 ). We first tested whether all the homologous cbTC terpene synthases functioned as cembrene B cyclases as we originally established in R. muelleri and E. caribaeorum . 27,28 We expressed each, 6 in total, in Saccharomyces cerevisiae and evaluated their production profiles by GC-MS analysis (Fig. 4 , SI Fig. 7 ). In all cases, we confirmed that GGPP ( 9 ) was cyclized into cembrene B ( 10 ) as the sole product. For the remaining enzyme assays, we tested the functions of genes from the R. koellikeri BGC and additional representative homologs identified by synteny. All tested enzyme homologs showed identical reactivity. We evaluated the functions of the three co-localized CYP subtypes in pairwise fashion by co-expression with cognate cbTC genes using previously established methods (Fig. 4 ). 50–52 We consistently measured the production of 19-hydroxycembrene B ( 11 ) in yeast upon co-expression with cbCYPb homologs. Its formation proceeded in a regiospecific manner with no other detectable monooxygenated products ( SI Fig. 8 ). Scale-up fermentation and chromatographic purification provided 11 in pure form, which we characterized by NMR spectroscopic analyses ( SI Tables 10 and 11, NMR supplementary note, SI Fig. 9 ). When we instead co-expressed cbCYPc homologs with cbTC genes, we observed formation of 17-hydroxycembrene ( S1 ) which could be detected by LCMS ( SI Fig. 10 ) and characterized by NMR ( SI Tables 11 and 12, NMR supplementary note, SI Fig. 9 ). Because cbCYPb and cbCYPc were found to oxidize complimentary positions on 10 with respect to the briarane pathway, we tested the cbCYPb product 11 by feeding it to yeast expressing cbCYPc homologs ( SI Fig. 11 ). We clearly observed a single oxidized product that upon purification and combined spectroscopic ( SI Tables 10 and 11, NMR supplementary note, SI Fig. 9 ) and x-ray crystallographic characterization ( SI Fig. 12, SI Table 12 ) was established as (7 S )-7,19-dihydroxycembrene B ( 12 ). Likewise, when yeast harboring cbCYPb genes were fed S1 , we similarly obtained 12 , albeit with lower rates of conversion (see Methods). The cooperative and related functions of the cbCYPb and cbCYPc enzymes correlates well with their shared evolution (Fig. 3 c). Despite numerous attempts with different substrate and gene pairs, we have yet to observe a function for the third cytochrome P450, cbCYPa. As a variety of cembrene double bond isomers are produced by corals, we considered the possibility that the described CYP450 enzymes could act non-specifically on other cembrene isomers. We thus co-expressed cbCYPb homologs with terpene synthases that selectively produce either cembrene A or C. However, no oxidized products were produced ( SI Fig. 13 ), showing high specificity of cbCYPb for cembrene B. We next evaluated the co-localized short-chain dehydrogenase cbSDH homologs for their possible involvement in further transforming the diol 12 . We expressed cbSDHs in Escherichia coli to conduct a series of in vitro incubation experiments. Incubation of purified cbSDH enzymes with NADP and 11 produced aldehyde ( S2 ), while incubation with the diol 12 produced cembrene B g-lactone ( 13 ), a structure feature diagnostic of all briarane diterpenoids ( SI Figs. 14 and 15 ). Two consecutive dehydrogenation steps were detected by LCMS chromatography via the formation of a m/z 301.2183 product ion ([M + H] + , calcd. m/z 301.2163 for C 20 H 29 O 2 + ). We scaled-up the in vitro reaction of 12 and NADP with the cbSDH to provide lactone 13 , which was purified and characterized spectroscopically. Scale up experiments were conducted with representative characterized genes to unambiguously determine all structures. The structure was supported by NMR ( SI Tables 11 and 12, NMR supplementary note, SI Fig. 9 ) and the diagnostic γ-lactone carbonyl stretch in the IR spectrum ( v max 1757 cm − 1 ). It is possible that 13 forms from 12 via a five-membered hemiacetal intermediate ( SI Fig. 16 ). Discussion We discovered a broadly conserved BGC family in corals and biochemically validated its function in the early biosynthetic steps toward briarane diterpenoids. The distinctive family of bioactive molecules numbers over 700 characterized members that are exclusively found in corals. Our work established the biosynthetic logic of briarane’s diagnostic γ-lactone core structure by four of the co-clustered gene products, cbTC, cbCYPb, cbCYPc, and cbSDH, confirming the clustered function of the briarane locus. The enzymatic transformations begin with cembrene B ( 10 ) formation by cbTC, then proceed with subsequent 19- and 7-hydroxylations carried out by cbCYPb and cbCYPc, respectively. Finally, in a unique mechanism of lactone biosynthesis, cbSDH transforms diol 12 directly into lactone 13 (Fig. 5 A), likely via aldehyde formation at the 19-position, hemiacetal formation and a second dehydrogenation to produce the diagnostic lactone ( SI Fig. 16 ). We have not yet characterized the enzyme(s) catalyzing bond formation between carbons 1 and 10. 29 A noncanonical class 2 TC acting on lactone 13 could introduce this feature to produce 14 (Fig. 5 A). 53 Alternatively, a Δ 1,14 epoxide intermediate, potentially installed by a P450, might precede a transannular cyclization reaction to 18 (Fig. 5 B). In animals, previous reports indicate that genes can be found in clusters when they result from successive tandem duplications or HGT. However, such factors are not at play here in the octocorals. Homologs of the briarane TC and CYP genes are ancient, found throughout octocorallia (> 500 Mya). The briarane BGC itself is likely evolutionarily ancient and evolved soon after the division of the octocoral orders Scleralcyonacea and Malacalcyonacea. It is only found in briarane-containing corals from Sclearcyonacea and has similar synteny in distantly related corals such as Renilla koellikeri (Family Pennatoloidea) and Erythropodium caribeorum (Family Erythropodiidae) (Fig. 5 C). 54 A key evolutionary event in the formation of the briarane BGC was the duplication of the related cbCYPb and -c genes early in the evolution of Scleralcyonacea. Supporting this timing, cbCYPb / -c pair is found across the Scleralcyonacea coral families Pennatuloidea, Ellisellidae, Coralliidae, Briareidae, and Erythropodiidae, but not elsewhere in sequenced octocorals. Subsequently, we speculate that cbTC, cbSDH, and cbCYPa were recruited and incorporated adjacent to the duplicated cbCYPb–c gene pair, also before speciation in the Scleralcyonacea, via chromosomal rearrangements to form the BGC — a phenomenon not previously reported in metazoans. Following the formation of the briarane BGC, all or parts of it were lost in a group of octocorals, as no briaranes have been reported from several Scleralcyonacea families (Fig. 5 C). Supporting the gene loss scenario, H. coerulea (family Helioporidae), does not harbor a cembrene B synthase in its genome but retains cbCYPb and -c relatives, suggesting that these genes were once part of a briarane BGC. More complex evolutionary scenarios are also observable in sequenced genomes, such as the BGCs from S. elongata and B. asbestinum , which contain additional SDH gene copies. Their appearance coincides with a split of the gene cluster into two parts. The observation of metazoan BGC formation and decay suggests potential mechanisms of chromosomal evolution not previously observed in animals. In contrast to a well preserved central briarane BGC, genes for specialized, late-stage tailoring processes are not colocalized. The core briarane skeleton is modified in many ways to create the > 700 known briarane derivatives (Fig. 5 A). For example, diterpenoids such as 11-hydroxyptilosarcenone ( 6 ), gemmacolide R ( 7 ), bathyptilone A ( 15 ), and erythrolide B ( 16 ), found in different coral species, have been modified by diverse enzymatic reactions including oxidation, chlorination, and acylation, by genes found in different chromosomal locations. This may represent a general rule in metazoan BGCs, in which the required core genes are found in the BGC, while those for variable and species-specific tailoring steps are distributed through the genomes. While the innate functions of coral terpenoids are largely unknown outside of a few examples having feeding deterrent properties, 55 we suspect that the conservation of their BGCs indicates central, unrealized biological roles pertaining to some of the more common coral terpenoid classes. 56 With the characterization of a coral terpenoid BGC, future genetic experiments are positioned to address their function once genome editing methods have been established in octocorals as recently accomplished in the hexacoral Acropora millepora. 57 The stage is set for genome mining and bioengineering studies to delineate biosynthetic cascade reactions to rare, bioactive coral molecules for their biosustainable production and application as therapeutic compounds. The briarane gene cluster family is likely the first of many to be described. We previously reported that some corals contain 70 kb and 500 kb putative BGCs with upward of 13 terpenoid biosynthesis genes. 27,28 With over 4,000 terpenoids characterized from corals (determined by searching the CMNPD database), 58 gene clustering is likely common in this ancient phylum. Looking ahead, we have yet to see if the active assembly of biosynthetic gene clusters by evolutionary forces is an exceedingly rare event that is limited to octocorals or if this is a more common phenomenon in animals and eukaryotes in general that has been overlooked so far due to limited genomic information in the eukaryotic kingdom. Declarations Acknowledgments We thank P. Zerofski (Scripps Institution of Oceanography) for collecting R. koellikeri and S. elongata , J. Sprung (Two Little Fishes Inc) for the D. gemmacea tissue, and J. Garrabou (CIIMAR) for the C. rubrum tissue. For her help in species identification, we would like to recognize C. S. McFadden (Harvey Mudd College). For coral images, we thank C. S. McFadden- D. gemmacea , F. Zuberer (CNRS)- C. rubrum , P. Webster @underwaterpat- S. elongata , and J. Simpson (Herbert Wertheim UF Scripps Institute for Biomedical Innovation & Technology)- B. asbestinum and E. caribaeorum . We thank T. Damiani (IOCB Prague) and A. Rodriguez (UC San Diego) for their assistance with feature based molecular networking. Chemical isolation and elucidation of the R. koellikeri briarane diterpene was aided by A. Bogdanov (Scripps Institution of Oceanography). This work was supported by the National Institutes of Health (R01-GM146224 to B.S.M., R35-GM148283 to E.W.S, and K99-GM148783 to P.D.S.), a Margaret A. Davidson Graduate Fellowship to N.E.G. (NERRS NA22NOS4200050), an NSF Graduate Research Fellowship to M.L.M., Tang Genomics Fund to T.P.M,, and Fundaço para a Ciência e a Tecnologia funds (UIDB/04423/2020, UIDP/04423/2020 and 2021.00855.CEECIND) to J.B.L. X-ray diffraction research reported in this publication was supported by the Office of the Director, NIH under award S10-OD030326. NMR data obtained at UC San Diego (the Scripps Institution of Oceanography and Biomolecular NMR Facilities) was assisted by B. Duggan. Additional NMR data was obtained at the University of Utah Health Sciences NMR Core. GCMS data was obtained at the University of Utah Health Sciences Proteomics Core by Q. Pearce. Author Contributions: Conceptualization: PDS, IB, EWS, BSM Data creation: NEG, PDS, MLM, IB Formal analysis and validation: NEG, PDS, IB, JGG, TA Funding acquisition: PDS, JBL, TPM, EWS, BSM Investigation and methodology: NEG, PDS, MLM, IB Project administration: PDS, IB, EWS, BSM Resources: JBL, TPM, EWS, BSM Writing – original draft: NEG, PDS, IB, BSM Writing – review and editing: All authors Competing interests: Authors declare that they have no competing interests. Data and materials availability: The feature based molecular network parameters used and raw data are available through the MASSIVE dataset ID: MSV000094792. Accessions numbers for the genomes generated in this paper are as follows; R. koellikeri (SAMN40621396) , B. asbestinum (SAMN40621398) , D. gemmacea (SAMN40621399) , and S. elongata (SAMN40621397). References Medema, M.H., de Rond, T., and Moore, B.S. (2021). Mining genomes to illuminate the specialized chemistry of life. Nat. Rev. Genet. 22 , 553–571. https://doi.org/10.1038/s41576-021-00363-7. Torres, J.P., and Schmidt, E.W. (2019). The biosynthetic diversity of the animal world. J. Biol. Chem. 294 , 17684–17692. https://doi.org/10.1074/jbc.REV119.006130. Poth, D., Wollenberg, K.C., Vences, M., and Schulz, S. (2012). Volatile Amphibian Pheromones: Macrolides from Mantellid Frogs from Madagascar. Angew. Chem. Int. Ed. 51 , 2187–2190. https://doi.org/10.1002/anie.201106592. Chinta, S.P., Goller, S., Lux, J., Funke, S., Uhl, G., and Schulz, S. (2010). The Sex Pheromone of the Wasp Spider Argiope bruennichi. Angew. Chem. Int. Ed. 49 , 2033–2036. https://doi.org/10.1002/anie.200906311. Masello, J.F., Lubjuhn, T., and Quillfeldt, P. (2008). Is the structural and psittacofulvin-based coloration of wild burrowing parrots Cyanoliseus patagonus condition dependent? J. Avian Biol. 39 , 653–662. https://doi.org/10.1111/j.1600-048X.2008.04417.x. Anderson, H.A., Mathieson, J.W., and Thomson, R.H. (1969). Distribution of spinochrome pigments in echinoids. Comp. Biochem. Physiol. 28 , 333–345. https://doi.org/10.1016/0010-406X(69)91347-4. Butcher, R.A., Fujita, M., Schroeder, F.C., and Clardy, J. (2007). Small-molecule pheromones that control dauer development in Caenorhabditis elegans. Nat. Chem. Biol. 3 , 420–422. https://doi.org/10.1038/nchembio.2007.3. Dettner, K., and Schwinger, G. (1982). Defensive secretions of three oxytelinae rove beetles (Coleoptera: Staphylinidae). J. Chem. Ecol. 8 , 1411–1420. https://doi.org/10.1007/BF01403104. Olivera, B.M., Gray, W.R., Zeikus, R., McIntosh, J.M., Varga, J., Rivier, J., de Santos, V., and Cruz, L.J. (1985). Peptide Neurotoxins from Fish-Hunting Cone Snails. Science 230 , 1338–1343. https://doi.org/10.1126/science.4071055. Chevrette, M.G., Gutiérrez-García, K., Selem-Mojica, N., Aguilar-Martínez, C., Yañez-Olvera, A., Ramos-Aboites, H.E., Hoskisson, P.A., and Barona-Gómez, F. (2020). Evolutionary dynamics of natural product biosynthesis in bacteria. Nat. Prod. Rep. 37 , 566–599. https://doi.org/10.1039/C9NP00048H. Rokas, A., Mead, M.E., Steenwyk, J.L., Raja, H.A., and Oberlies, N.H. (2020). Biosynthetic gene clusters and the evolution of fungal chemodiversity. Nat. Prod. Rep. 37 , 868–878. https://doi.org/10.1039/C9NP00045C. Nützmann, H.-W., Scazzocchio, C., and Osbourn, A. (2018). Metabolic Gene Clusters in Eukaryotes. Annu. Rev. Genet. 52 , 159–183. https://doi.org/10.1146/annurev-genet-120417-031237. Steele, T.S., Brunson, J.K., Maeno, Y., Terada, R., Allen, A.E., Yotsu-Yamashita, M., Chekan, J.R., and Moore, B.S. (2022). Domoic acid biosynthesis in the red alga Chondria armata suggests a complex evolutionary history for toxin production. Proc. Natl. Acad. Sci. 119 , e2117407119. https://doi.org/10.1073/pnas.2117407119. Smit, S.J., and Lichman, B.R. (2022). Plant biosynthetic gene clusters in the context of metabolic evolution. Nat. Prod. Rep. 39 , 1465–1482. https://doi.org/10.1039/D2NP00005A. Ryngajłło, M., Boruta, T., and Bizukojć, M. (2021). Complete genome sequence of lovastatin producer Aspergillus terreus ATCC 20542 and evaluation of genomic diversity among A. terreus strains. Appl. Microbiol. Biotechnol. 105 , 1615–1627. https://doi.org/10.1007/s00253-021-11133-0. Itkin, M., Heinig, U., Tzfadia, O., Bhide, A.J., Shinde, B., Cardenas, P.D., Bocobza, S.E., Unger, T., Malitsky, S., Finkers, R., et al. (2013). Biosynthesis of Antinutritional Alkaloids in Solanaceous Crops Is Mediated by Clustered Genes. Science 341 , 175–179. https://doi.org/10.1126/science.1240230. Nützmann, H.-W., Huang, A., and Osbourn, A. (2016). Plant metabolic clusters – from genetics to genomics. New Phytol. 211 , 771–789. https://doi.org/10.1111/nph.13981. Osborn, A.R., Almabruk, K.H., Holzwarth, G., Asamizu, S., LaDu, J., Kean, K.M., Karplus, P.A., Tanguay, R.L., Bakalinsky, A.T., and Mahmud, T. (2015). De novo synthesis of a sunscreen compound in vertebrates. eLife 4 , e05919. https://doi.org/10.7554/eLife.05919. Roelofs, D., Timmermans, M.J.T.N., Hensbergen, P., van Leeuwen, H., Koopman, J., Faddeeva, A., Suring, W., de Boer, T.E., Mariën, J., Boer, R., et al. (2013). A Functional Isopenicillin N Synthase in an Animal Genome. Mol. Biol. Evol. 30 , 541–548. https://doi.org/10.1093/molbev/mss269. Faddeeva-Vakhrusheva, A., Kraaijeveld, K., Derks, M.F.L., Anvar, S.Y., Agamennone, V., Suring, W., Kampfraath, A.A., Ellers, J., Le Ngoc, G., van Gestel, C.A.M., et al. (2017). Coping with living in the soil: the genome of the parthenogenetic springtail Folsomia candida. BMC Genomics 18 , 493. https://doi.org/10.1186/s12864-017-3852-x. Feng, L., Gordon, M.T., Liu, Y., Basso, K.B., and Butcher, R.A. (2021). Mapping the biosynthetic pathway of a hybrid polyketide-nonribosomal peptide in a metazoan. Nat. Commun. 12 , 4912. https://doi.org/10.1038/s41467-021-24682-9. de Rosa, R., Grenier, J.K., Andreeva, T., Cook, C.E., Adoutte, A., Akam, M., Carroll, S.B., and Balavoine, G. (1999). Hox genes in brachiopods and priapulids and protostome evolution. Nature 399 , 772–776. https://doi.org/10.1038/21631. Storz, J.F., Opazo, J.C., and Hoffmann, F.G. (2013). Gene duplication, genome duplication, and the functional diversification of vertebrate globins. Mol. Phylogenet. Evol. 66 , 469–478. https://doi.org/10.1016/j.ympev.2012.07.013. Kitchen, S.A., Naragon, T.H., Brückner, A., Ladinsky, M.S., Quinodoz, S.A., Badroos, J.M., Viliunas, J.W., Kishi, Y., Wagner, J.M., Miller, D.R., et al. (2024). The genomic and cellular basis of biosynthetic innovation in rove beetles. Cell 0 . https://doi.org/10.1016/j.cell.2024.05.012. Turcu, A.F., Rege, J., Auchus, R.J., and Rainey, W.E. (2020). 11-Oxygenated androgens in health and disease. Nat. Rev. Endocrinol. 16 , 284–296. https://doi.org/10.1038/s41574-020-0336-x. Tan, D.-X., Hardeland, R., Back, K., Manchester, L.C., Alatorre-Jimenez, M.A., and Reiter, R.J. (2016). On the significance of an alternate pathway of melatonin synthesis via 5-methoxytryptamine: comparisons across species. J. Pineal Res. 61 , 27–40. https://doi.org/10.1111/jpi.12336. Burkhardt, I., de Rond, T., Chen, P.Y.-T., and Moore, B.S. (2022). Ancient plant-like terpene biosynthesis in corals. Nat. Chem. Biol. 18 , 664–669. https://doi.org/10.1038/s41589-022-01026-2. Scesa, P.D., Lin, Z., and Schmidt, E.W. (2022). Ancient defensive terpene biosynthetic gene clusters in the soft corals. Nat. Chem. Biol. 18 , 659–663. https://doi.org/10.1038/s41589-022-01027-1. Moon, N.G., and Harned, A.M. (2018). Synthetic explorations of the briarane jungle: progress in developing a synthetic route to a common family of diterpenoid natural products. R. Soc. Open Sci. 5 , 172280. https://doi.org/10.1098/rsos.172280. La, M.-P., Li, J., Li, C., Tang, H., Liu, B.-S., Sun, P., Zhuang, C.-L., Li, T.-J., and Zhang, W. (2014). Briarane Diterpenoids from the Gorgonian Dichotella gemmacea. Mar. Drugs 12 , 6178–6189. https://doi.org/10.3390/md12126178. Barsby, T., and Kubanek, J. (2005). Isolation and Structure Elucidation of Feeding Deterrent Diterpenoids from the Sea Pansy, Renilla reniformis. J. Nat. Prod. 68 , 511–516. https://doi.org/10.1021/np049609u. Fenical, W., and Pawlik, J.R. (1991). Defensive properties of secondary metabolites from the Caribbean gorgonian coral Erythropodium caribaeorum. Mar. Ecol. Prog. Ser. 75 , 1–8. Harvell, C.D., Fenical, W., Roussis, V., Ruesink, J.L., Griggs, C.C., and Greene, C.H. (1993). Local and geographic variation in the defensive chemistry of a West Indian gorgonian coral (Briareum asbestinum). Mar. Ecol. Prog. Ser. 93 , 165–173. Ksebati, M.B., and Schmitz, F.J. (1986). Diterpenes from a Soft Coral, Minabea SP., from Truk Lagoon. Bull. Sociétés Chim. Belg. 95 , 835–851. https://doi.org/10.1002/bscb.19860950910. Wratten, S.J., and Faulkner, D.J. (1979). Some diterpenes from the sea pen stylatula SP. Tetrahedron 35 , 1907–1912. https://doi.org/10.1016/0040-4020(79)85047-4. Yang, J., Zhang, S., Qia, S., Pan, J., Qiu, Y., Tao, S., Yin, H., and Li, Q. (2007). Briarane-type diterpenoids from the China gorgonian coral Subergorgia reticulata . Biochem. Syst. Ecol. 35 , 770–773. https://doi.org/10.1016/j.bse.2007.06.003. MarinLit - A database of the marine natural products literature https://marinlit.rsc.org/. Keifer, P.A., Rinehart, K.L., and Hooper, I.R. (1986). Renillafoulins, antifouling diterpenes from the sea pansy Renilla reniformis (Octocorallia). J. Org. Chem. 51 , 4450–4454. https://doi.org/10.1021/jo00373a020. Hendrickson, R.L., and Cardellina, J.H. (1986). Structure and stereochemistry of insecticidal diterpenes from the sea penXXX. Tetrahedron 42 , 6565–6570. https://doi.org/10.1016/S0040-4020(01)88120-5. Nurco, D.J., Conklin, D.E., Shapiro, N.S., and Tran, E. (2010). The absolute structure of ptilosarcenone 2.5-hydrate, a diterpenoid briarane from the orange sea pen Ptilosarcus gurneyi (Gray). Acta Crystallogr. Sect. E Struct. Rep. Online 67 , o181–o182. https://doi.org/10.1107/S1600536810050142. Li, C., La, M.-P., Sun, P., Kurtan, T., Mandi, A., Tang, H., Liu, B.-S., Yi, Y.-H., Li, L., and Zhang, W. (2011). Bioactive (3Z,5E)-11,20-Epoxybriara-3,5-dien-7,18-olide Diterpenoids from the South China Sea Gorgonian Dichotella gemmacea. Mar. Drugs 9 , 1403–1418. https://doi.org/10.3390/md9081403. Pordesimo, E.O., Schmitz, F.J., Ciereszko, L.S., Hossain, M.B., and Van der Helm, D. (1991). New briarein diterpenes from the Caribbean gorgonians Erythropodium caribaeorum and Briareum sp. J. Org. Chem. 56 , 2344–2357. https://doi.org/10.1021/jo00007a020. Molina, S.L., Forero, A.M., Ayala, F.I., Puyana, M., Zea, S., Castellanos, L., Muñoz, D., Arboleda, G., Sandoval-Hernández, A.G., and Ramos, F.A. (2020). Metabolic Profiling of the Soft Coral Erythropodium caribaeorum (Alcyonacea: Anthothelidae) from the Colombian Caribbean Reveals Different Chemotypes. Mar. Drugs 18 , 4. https://doi.org/10.3390/md18010004. Jiang, J.B., Quattrini, A.M., Francis, W.R., Ryan, J.F., Rodríguez, E., and McFadden, C.S. (2019). A hybrid de novo assembly of the sea pansy (Renilla muelleri) genome. GigaScience 8 , giz026. https://doi.org/10.1093/gigascience/giz026. Ledoux, J.-B., Gomez-Garrido, J., Cruz, F., Ferreira, F.C., Matos, A., Sarropoulou, X., Ramirez-Calero, S., Aurelle, D., Lopez-Sendino, P., Grayson, N., et al. (2024). Chromosome-level genome assembly and annotation of Corallium rubrum: a Mediterranean coral threatened by overharvesting and climate change. Preprint at bioRxiv, https://doi.org/10.1101/2024.07.13.603384 https://doi.org/10.1101/2024.07.13.603384. Corominas, M., Marquès-Bonet, T., Arnedo, M.A., Bayés, M., Belmonte, J., Escrivà, H., Fernández, R., Gabaldón, T., Garnatje, T., Germain, J., et al. (2024). The Catalan initiative for the Earth BioGenome Project: contributing local data to global biodiversity genomics. NAR Genomics Bioinforma. 6 , lqae075. https://doi.org/10.1093/nargab/lqae075. Hu, M., Zheng, X., Fan, C.-M., and Zheng, Y. (2020). Lineage dynamics of the endosymbiotic cell type in the soft coral Xenia. Nature 582 , 534–538. https://doi.org/10.1038/s41586-020-2385-7. Ip, J.C.-H., Ho, M.-H., Chan, B.K.K., and Qiu, J.-W. (2023). A draft genome assembly of reef-building octocoral Heliopora coerulea. Sci. Data 10 , 381. https://doi.org/10.1038/s41597-023-02291-z. Pankov, K.V., McArthur, A.G., Gold, D.A., Nelson, D.R., Goldstone, J.V., and Wilson, J.Y. (2021). The cytochrome P450 (CYP) superfamily in cnidarians. Sci. Rep. 11 , 9834. https://doi.org/10.1038/s41598-021-88700-y. DeJong, J.M., Liu, Y., Bollon, A.P., Long, R.M., Jennewein, S., Williams, D., and Croteau, R.B. (2006). Genetic engineering of taxol biosynthetic genes in Saccharomyces cerevisiae. Biotechnol. Bioeng. 93 , 212–224. https://doi.org/10.1002/bit.20694. Pompon, D., Louerat, B., Bronine, A., and Urban, P. (1996). [6] Yeast expression of animal and plant P450s in optimized redox environments. In Methods in Enzymology, E. F. Johnson and M. R. Waterman, eds. (Academic Press), pp. 51–64. https://doi.org/10.1016/S0076-6879(96)72008-6. Ro, D.-K., Paradise, E.M., Ouellet, M., Fisher, K.J., Newman, K.L., Ndungu, J.M., Ho, K.A., Eachus, R.A., Ham, T.S., Kirby, J., et al. (2006). Production of the antimalarial drug precursor artemisinic acid in engineered yeast. Nature 440 , 940–943. https://doi.org/10.1038/nature04640. Rudolf, J.D., and Chang, C.-Y. (2020). Terpene synthases in disguise: enzymology, structure, and opportunities of non-canonical terpene synthases. Nat. Prod. Rep. 37 , 425–463. https://doi.org/10.1039/C9NP00051H. McFadden, C.S., Ofwegen, L.P. van, and Quattrini, A.M. (2022). Revisionary systematics of Octocorallia (Cnidaria: Anthozoa) guided by phylogenomics. Bull. Soc. Syst. Biol. 1 . https://doi.org/10.18061/bssb.v1i3.8735. Barsby, T., and Kubanek, J. (2005). Isolation and structure elucidation of feeding deterrent diterpenoids from the sea pansy, Renilla reniformis. J Nat Prod 68 , 511–516. https://doi.org/10.1021/np049609u. He, J., Xin, P., Ma, X., Chu, J., and Wang, G. (2020). Gibberellin Metabolism in Flowering Plants: An Update and Perspectives. Front. Plant Sci. 11 . https://doi.org/10.3389/fpls.2020.00532. Cleves, P.A., Strader, M.E., Bay, L.K., Pringle, J.R., and Matz, M.V. (2018). CRISPR/Cas9-mediated genome editing in a reef-building coral. Proc. Natl. Acad. Sci. 115 , 5235–5240. https://doi.org/10.1073/pnas.1722151115. Lyu, C., Chen, T., Qiang, B., Liu, N., Wang, H., Zhang, L., and Liu, Z. (2021). CMNPD: a comprehensive marine natural products database towards facilitating drug discovery from the ocean. Nucleic Acids Res. 49 , D509–D515. https://doi.org/10.1093/nar/gkaa763. Additional Declarations There is NO Competing Interest. Supplementary Files BriaraneTerpenoidSI080324.pdf Supporting Information Cite Share Download PDF Status: Published Journal Publication published 13 Jun, 2025 Read the published version in Nature Chemical Biology → Version 1 posted You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-4859447","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Biological Sciences - Article","associatedPublications":[],"authors":[{"id":337754915,"identity":"935a9a20-3ea2-48fa-a924-7942caee875c","order_by":0,"name":"Bradley Moore","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAAsklEQVRIiWNgGAWjYBACAygtByIkGBiYiddiTLqWxAaitZjzn0588HGPTfqG48cf3mCosAbrxQssZ+RuNpzxLC13w5kcYwuGM+mEtRjc4N0mzXPgcO6GGzxsEoxth4nQcv4sSMv/dIMb7M8kGP8Ro+VALkjLgQSDGwxmEowNRGiB+OVAsuFMkF8SjqUbE9Rizn9244MPB+zk+UAh9qHGWpagFlSQQJryUTAKRsEoGAW4AABOOUEisHSNjQAAAABJRU5ErkJggg==","orcid":"https://orcid.org/0000-0002-4652-1253","institution":"University of California, San Diego","correspondingAuthor":true,"prefix":"","firstName":"Bradley","middleName":"","lastName":"Moore","suffix":""},{"id":337754916,"identity":"944a1fa9-d502-4eec-8559-f935e94708a0","order_by":1,"name":"Natalie Grayson","email":"","orcid":"","institution":"University of California, San Diego","correspondingAuthor":false,"prefix":"","firstName":"Natalie","middleName":"","lastName":"Grayson","suffix":""},{"id":337754918,"identity":"86416d74-e5c5-476d-af84-3ce34590438e","order_by":2,"name":"Paul Scesa","email":"","orcid":"","institution":"University of Utah","correspondingAuthor":false,"prefix":"","firstName":"Paul","middleName":"","lastName":"Scesa","suffix":""},{"id":337754920,"identity":"e1035f52-cc0f-44ee-97e6-404ebbbfbf77","order_by":3,"name":"Malia Moore","email":"","orcid":"","institution":"University of California, San Diego","correspondingAuthor":false,"prefix":"","firstName":"Malia","middleName":"","lastName":"Moore","suffix":""},{"id":337754922,"identity":"32ef6185-2d5f-411d-839f-98487c75ab56","order_by":4,"name":"Jean-Baptiste Ledoux","email":"","orcid":"https://orcid.org/0000-0001-8796-6163","institution":"CIIMAR-Interdisciplinary centre of marine and environmental research - University of Porto. Terminal de Cruzeiros de Leixões, Av. General Norton de Matos S/N 4450-208 Matosinhos, Portugal","correspondingAuthor":false,"prefix":"","firstName":"Jean-Baptiste","middleName":"","lastName":"Ledoux","suffix":""},{"id":337754923,"identity":"0f72e57f-d2e3-4939-9bb9-af5c4d8a7d82","order_by":5,"name":"Jessica Gomez-Garrido","email":"","orcid":"","institution":"Centro Nacional de Análisis Genómico","correspondingAuthor":false,"prefix":"","firstName":"Jessica","middleName":"","lastName":"Gomez-Garrido","suffix":""},{"id":337754925,"identity":"bd559655-f141-4da2-be6a-0421009c086a","order_by":6,"name":"Tyler Alioto","email":"","orcid":"","institution":"Centro Nacional de Análisis Genómico","correspondingAuthor":false,"prefix":"","firstName":"Tyler","middleName":"","lastName":"Alioto","suffix":""},{"id":337754926,"identity":"53b4e0a3-e01c-401d-9869-b36000bfce29","order_by":7,"name":"Todd Michael","email":"","orcid":"https://orcid.org/0000-0001-6272-2875","institution":"Salk Institute for Biological Studies","correspondingAuthor":false,"prefix":"","firstName":"Todd","middleName":"","lastName":"Michael","suffix":""},{"id":337754927,"identity":"9722f61e-40f2-4405-8528-c3555e0e4b6d","order_by":8,"name":"Immo Burkhardt","email":"","orcid":"https://orcid.org/0000-0001-9515-4042","institution":"University of California, San Diego","correspondingAuthor":false,"prefix":"","firstName":"Immo","middleName":"","lastName":"Burkhardt","suffix":""},{"id":337754929,"identity":"5bffccb3-d531-46cb-8d19-4cb197ad27ba","order_by":9,"name":"Eric Schmidt","email":"","orcid":"https://orcid.org/0000-0001-5839-694X","institution":"Univeristy of Utah","correspondingAuthor":false,"prefix":"","firstName":"Eric","middleName":"","lastName":"Schmidt","suffix":""}],"badges":[],"createdAt":"2024-08-05 05:50:19","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-4859447/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-4859447/v1","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1038/s41589-025-01927-y","type":"published","date":"2025-06-13T04:00:00+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":62927132,"identity":"16e9ec7c-7dfe-43c8-b01a-59fa31defb1c","added_by":"auto","created_at":"2024-08-21 06:57:36","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":196390,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eGene clustering trends across eukaryotes\u003c/strong\u003e. Fungal gene clusters are often relatively compact and encode a complete biosynthetic pathway as exemplified in the lovastatin (\u003cstrong\u003e1\u003c/strong\u003e) BGC from \u003cem\u003eA. terreus\u003c/em\u003e. Plant BGCs often encode only partial pathways and are more extended across the chromosomes compared to fungal BGCs, as shown for the two partial tomatine (\u003cstrong\u003e2\u003c/strong\u003e) BGCs from \u003cem\u003eS. lycopersicum\u003c/em\u003e. Genes coding for specialized metabolite biosynthesis in animals sometimes occur in clusters like the genes for gadusol (\u003cstrong\u003e3\u003c/strong\u003e) production, but they show high similarity and synteny to microbial genes and are therefore likely a result of horizontal gene transfer. Other animal genes that are not HGT products are usually not colocalized like for nemamide (\u003cstrong\u003e4\u003c/strong\u003e) or melatonin (\u003cstrong\u003e5\u003c/strong\u003e) production in \u003cem\u003eC. elegans\u003c/em\u003e and \u003cem\u003eH. sapiens\u003c/em\u003e, respectively. The investigation and characterization of a relatively tight BGC occurring in numerous octocorals like \u003cem\u003eR. koellikeri\u003c/em\u003e, putatively encoding enzymes to produce briaranes like\u003cstrong\u003e 6\u003c/strong\u003e is the focus of this work.\u003c/p\u003e","description":"","filename":"image1.png","url":"https://assets-eu.researchsquare.com/files/rs-4859447/v1/ff50b9c4f70c73c9f770c8fd.png"},{"id":62927136,"identity":"6760a28e-27db-4efd-994b-ec1abb5be3f1","added_by":"auto","created_at":"2024-08-21 06:57:36","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":447090,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eDistribution of briarane diterpenoids in the phylum octocorallia. A:\u003c/strong\u003e Phylogeny of octocoral rooted at Scleralcyonacea. Families containing known briarane producer species are denoted with a star. Families with reported other terpenoid chemistry are marked with triangles, while families with no reported chemical data are marked with squares. \u003cstrong\u003eB: \u003c/strong\u003eOrganisms used in the molecular network. \u003cstrong\u003eC: \u003c/strong\u003eThe briarane terpenoid-containing molecular network. Highlighted nodes include 11-hydroxyptilosarcenone (\u003cstrong\u003e6\u003c/strong\u003e) from \u003cem\u003eR. koellikeri, \u003c/em\u003eerythrolide D (\u003cstrong\u003e8\u003c/strong\u003e) from \u003cem\u003eE. caribaeorum\u003c/em\u003e, and gemmacolide R (\u003cstrong\u003e7\u003c/strong\u003e) from \u003cem\u003eD. gemmacea\u003c/em\u003e. The full molecular network is shown in \u003cstrong\u003eSI Figure 1\u003c/strong\u003e.\u003c/p\u003e","description":"","filename":"image2.png","url":"https://assets-eu.researchsquare.com/files/rs-4859447/v1/fcaafebc731e1c5ada213940.png"},{"id":62927135,"identity":"646c9449-8774-49e3-92b7-b0a218011710","added_by":"auto","created_at":"2024-08-21 06:57:36","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":103878,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eGenome guided discovery of a conserved octocoral BGC family. A:\u003c/strong\u003e comparison of BUSCO and contig N50 data from 17 octocoral genomes. \u003cstrong\u003eB: \u003c/strong\u003eSynteny of the briarane BGCs in Scleralcyonacea corals. Intron sequences have been removed for ease of syntenic visualization; for full BGC to scale depiction, see \u003cstrong\u003eSI Figure 4\u003c/strong\u003e. \u003cstrong\u003eC: \u003c/strong\u003ePhylogeny of cnidarian CYPs. The tree includes all CYPs from eight octocorals (\u003cem\u003eH. coerulea, D. gemmacea, B. asbestinum, R. koellikeri, S. elongata, C. rubrum, E. caribaeorum, \u003c/em\u003eand \u003cem\u003eXenia \u003c/em\u003esp.)\u003cem\u003e \u003c/em\u003eand from four non-octocoral cnidarians (\u003cem\u003eH. vulgaris, A. digitifera, A. aurita, N. vectensis\u003c/em\u003e) as reference. Previously established CYP clans\u003csup\u003e49\u003c/sup\u003e\u0026nbsp;are marked with blue and contain octocoral and other cnidarian CYPs. Unmarked clades contain only octocoral CYP sequences. \u003cstrong\u003eD: \u003c/strong\u003eZoomed-in branches of the briarane BGC cbCYPa-c’s distinct clades, highlighted coral (cbCYPa), lavender (cbCYPb), and mustard (cbCYPc), labeled accordingly with octocoral species initials followed by enzyme type.\u003c/p\u003e","description":"","filename":"image3.png","url":"https://assets-eu.researchsquare.com/files/rs-4859447/v1/47b8b2b077f08f20c74cec58.png"},{"id":62927131,"identity":"4b45cf93-19da-44fc-94ef-ee895fa20a0c","added_by":"auto","created_at":"2024-08-21 06:57:36","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":234459,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eBiochemical validation of clustered biosynthesis genes across syntenic homologs.\u003c/strong\u003e \u003cstrong\u003eA\u003c/strong\u003e Annotation of the \u003cem\u003eR. koelliker\u003c/em\u003ei briarane BGC. This prototypical example includes the terpene cyclase cbTC, the three different CYP450 subtypes cbCYPa, cbCYPb and cbCYPc, and the short-chain dehydrogenase cbSDH. \u003cstrong\u003eB\u003c/strong\u003e Scheme depicting the overall pathway to cembrene B γ-lactone \u003cstrong\u003e13\u003c/strong\u003e. This pathway is mediated stepwise by initial cyclization of geranylgeranyl diphosphate (GGPP-\u003cstrong\u003e9\u003c/strong\u003e) to \u003cstrong\u003e10\u003c/strong\u003e by cbTC, followed by mono-oxygenation to \u003cstrong\u003e11\u003c/strong\u003e by cbCYPb. Mono-oxygenation of \u003cstrong\u003e11\u003c/strong\u003e to diol \u003cstrong\u003e12\u003c/strong\u003e followed by two consecutive dehydrogenations by cbSDH to γ-lactone \u003cstrong\u003e13\u003c/strong\u003e. \u003cstrong\u003eC\u003c/strong\u003e Chromatographic data for enzymatic assays demonstrating the depicted reactions. From left to right: GCMS EIC trace at 272.1 \u003cem\u003em/z\u003c/em\u003e showing the formation of \u003cstrong\u003e10\u003c/strong\u003e; LCMS EIC trace at 271.2 \u003cem\u003em/z\u003c/em\u003e ([M-OH]\u003csup\u003e+\u003c/sup\u003e ion) showing \u003cstrong\u003e11; \u003c/strong\u003eLCMS EIC trace at 287.2 \u003cem\u003em/z\u003c/em\u003e ([M-OH]\u003csup\u003e+\u003c/sup\u003e ion) showing \u003cstrong\u003e12\u003c/strong\u003e; and LCMS EIC trace at 301.2 \u003cem\u003em/z\u003c/em\u003e ([M+H]\u003csup\u003e+\u003c/sup\u003e ion) showing the formation of \u003cstrong\u003e13. \u003c/strong\u003eAll traces were compared to the relevant purified reference compound (see\u003cstrong\u003e SI Figures 7,8, 11, 15\u003c/strong\u003e) and negative control (bottom trace). The negative control for the SDH assay was reaction mixture without enzyme; for all others, the negative control was the crude extract of yeast harboring empty vector. The full four gene pathways from \u003cem\u003eR. koellikeri\u003c/em\u003e and \u003cem\u003eD. gemmacea\u003c/em\u003e were characterized as well as additional representative homologs identified by synteny. All tested enzyme homologs showed identical reactivity; n.t. = not tested.\u003c/p\u003e","description":"","filename":"image4.png","url":"https://assets-eu.researchsquare.com/files/rs-4859447/v1/70cf5d59f672703b5c14691d.png"},{"id":62927613,"identity":"42cb097d-db5d-4378-bb9a-e681658e3b2a","added_by":"auto","created_at":"2024-08-21 07:05:36","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":209173,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eBiosynthetic pathway for briarane production and estimated timepoints for the evolution of the briarane BGC.\u003c/strong\u003e \u003cstrong\u003eA\u003c/strong\u003eClustered enzymes characterized in this work produce the lactone \u003cstrong\u003e13\u003c/strong\u003e, a central intermediate in briarane biosynthesis by subsequent terpene cyclization, dual hydroxylation and two dehydrogenations. A missing central transformation is a 1,10-cyclization, the last fundamental step to obtain a compound containing the complete briarane scaffold (\u003cstrong\u003e14\u003c/strong\u003e). These steps must be conserved in all briarane producers (blue background). Other modifications of the briarane core are most likely introduced later and are lineage and species specific (purple background). These biosynthetic steps are evolutionarily younger and might not be present in gene clusters. \u003cstrong\u003eB\u003c/strong\u003e A putative cytochrome P450 mediated 1,10-ring closing reaction, initiated by epoxidation of \u003cstrong\u003e13\u003c/strong\u003e to \u003cstrong\u003e17\u003c/strong\u003e, which could be opened to produce alcohol \u003cstrong\u003e18\u003c/strong\u003e bearing the briarane scaffold. The resulting oxidized 14 position occurs as either oxygenated or olefinic position in most briaranes. \u003cstrong\u003eC\u003c/strong\u003e Simplified octocoral phylogeny, rooted at Scleralcyonacea. Estimated timescale of BGC evolution. The circular mark in purple denotes the latest acquisition time of the briarane BGC, given that Ellisellidae, Pennatuloidea, as well as Briareidae, and Erythropodiidae (members of the top 9 families) are shown to contain the briarane BGC and are known briarane producers. The purple-colored X-mark denotes a likely BGC deletion event, in which the whole gene cluster or parts of it were lost. None of the 8 families in this clade are reported to produce briaranes and the only genome-sequenced species, \u003cem\u003eH. coerulea\u003c/em\u003e, does not harbor a cembrene B synthase.\u003c/p\u003e","description":"","filename":"image5.png","url":"https://assets-eu.researchsquare.com/files/rs-4859447/v1/36284d6989458d957e4b48b1.png"},{"id":84597026,"identity":"cc522b63-6d2b-4341-bb46-9a397b7de7ef","added_by":"auto","created_at":"2025-06-14 07:07:08","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":2223170,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-4859447/v1/30322fcf-6cfb-4237-9d6c-f6f9ce3ce3b0.pdf"},{"id":62927137,"identity":"9487ecff-c370-4304-8b97-703391575860","added_by":"auto","created_at":"2024-08-21 06:57:36","extension":"pdf","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":7459844,"visible":true,"origin":"","legend":"Supporting Information","description":"","filename":"BriaraneTerpenoidSI080324.pdf","url":"https://assets-eu.researchsquare.com/files/rs-4859447/v1/c3f70431bb0b2416d8c1f06b.pdf"}],"financialInterests":"There is \u003cb\u003eNO\u003c/b\u003e Competing Interest.","formattedTitle":"A widespread metabolic gene cluster family in metazoans","fulltext":[{"header":"Introduction","content":"\u003cp\u003eSpecialized metabolites, or natural products, are small molecules that while not vital for primary growth, enable a producing organism to survive in an ecological niche by mediating time-, tissue-, or stimulus-dependent functions.\u003csup\u003e1\u003c/sup\u003e Although most of these compounds are known from plants and microbes, they are widespread in all kingdoms of life, including animals, where they too serve diverse purposes including communication, development, and defense.\u003csup\u003e2\u003c/sup\u003e Examples include volatile sex-pheromones from frogs\u003csup\u003e3\u003c/sup\u003e and spiders,\u003csup\u003e4\u003c/sup\u003e polyketide pigments in parrots\u003csup\u003e5\u003c/sup\u003e and sea urchins,\u003csup\u003e6\u003c/sup\u003e sugar-derived hormones in nematodes,\u003csup\u003e7\u003c/sup\u003e hydroquinones employed by rove beetles to fend off predators,\u003csup\u003e8\u003c/sup\u003e and peptide toxins used by cone snails to stun their prey.\u003csup\u003e9\u003c/sup\u003e\u003c/p\u003e \u003cp\u003eThe overwhelming majority of research on specialized metabolite biosynthesis has been conducted in bacteria and fungi, where genes coding for a biosynthetic pathway are typically arranged in physically co-localized, biosynthetic gene clusters (BGCs).\u003csup\u003e10,11\u003c/sup\u003e In contrast, the few characterized biosynthetic pathways in animals are not encoded in BGCs, but have their encoding genes dispersed throughout the genome. This is currently regarded as the general rule for biosynthetic pathways in non-fungal eukaryotes.\u003csup\u003e12\u003c/sup\u003e However, improved genomic resources in other eukaryotes like plants and algae has revealed the colocalization of non-paralogous genes into BGCs, where they code for specialized metabolic pathways (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e).\u003csup\u003e1,13,14\u003c/sup\u003e\u003c/p\u003e \u003cp\u003eDue to their abundance, eukaryotic BGCs were first recognized in fungi where they are now regarded as the norm in specialized metabolite biosynthesis. A prominent example is the lovastatin (\u003cb\u003e1\u003c/b\u003e) BGC, which codes for nine non-homologous genes within a compact genomic region of 40 kilobase pairs (kb) and encodes the complete set of enzymes needed to produce \u003cb\u003e1\u003c/b\u003e, a resistance gene, a transporter, and a transcription factor (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e).\u003csup\u003e15\u003c/sup\u003e The more recently discovered BGCs in plant genomes show clear differences from fungal BGCs. As generally the case in plant genomes, they show larger intergenic regions and often do not code for the complete biosynthetic pathway. An example is the tomatine (\u003cb\u003e2\u003c/b\u003e) BGC, which codes for six steps in the pathway starting from cholesterol.\u003csup\u003e16\u003c/sup\u003e However, the first enzymatic step of the pathway is encoded by a gene residing over 7,000 kb away from this cluster, while two additional genes are coded as a gene pair on a different chromosome (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e).\u003csup\u003e17\u003c/sup\u003e There are only a few examples of small molecule BGCs in animals that have since been identified. However, their origins link to microbial gene sets acquired via horizontal gene transfer (HGT) as in the case of mycosporine-like amino acids such as gadusol (\u003cb\u003e3\u003c/b\u003e) that widely occur in vertebrates\u003csup\u003e18\u003c/sup\u003e (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e) and a penicillin biosynthesis gene pair from the springtail \u003cem\u003eFolsomia candida\u003c/em\u003e that shows high homology to fungal sequences.\u003csup\u003e19,20\u003c/sup\u003e In contrast, the genes involved in the biosynthesis of the hybrid polyketide/ non-ribosomal peptide hormone nemamide (\u003cb\u003e4\u003c/b\u003e) from the nematode \u003cem\u003eCaenorhabditis elegans\u003c/em\u003e do not show clear sequence homology with microbial genes and thus may have evolved de novo in nematodes.\u003csup\u003e21\u003c/sup\u003e However, none of these genes are colocalized in a cluster, but dispersed over large genomic regions across four separate chromosomes.\u003c/p\u003e \u003cp\u003eAlthough uncommon, a few gene clusters in animal genomes that are not the result of horizontal gene transfer are known, such as the Hox \u003csup\u003e22\u003c/sup\u003eand b-globin\u003csup\u003e23\u003c/sup\u003e clusters in humans, as well as a laccase gene cluster recently discovered in rove beetles.\u003csup\u003e24\u003c/sup\u003e However, these gene clusters are made up of paralogous genes and have arisen from successive tandem duplications from one ancestor gene, and only the last example is involved in small molecule biosynthesis. Generally, the few characterized biosynthetic pathways in animals, including the extraordinarily well investigated examples of steroid hormones\u003csup\u003e25\u003c/sup\u003e or melatonin (\u003cb\u003e5\u003c/b\u003e)\u003csup\u003e26\u003c/sup\u003e in humans, are not encoded in BGCs, but have their encoding genes dispersed throughout the genome.\u003c/p\u003e \u003cp\u003eIt was therefore surprising when we recently identified potential terpenoid BGCs in octocorals, consisting of up to 13 co-localized genes spanning up to 500 kb that did not result from HGT or pure tandem duplications.\u003csup\u003e27,28\u003c/sup\u003e Instead, over evolutionary time, coral gene paralogs from different protein families appear to have been recruited to single genomic scaffolds, forming potential BGCs. We observed related putative gene clusters from \u003cem\u003eRenilla muelleri\u003c/em\u003e and \u003cem\u003eErythropodium caribeorum\u003c/em\u003e, corals that belong to different families (Pennatuloidea vs. Erythropodiidae) and differ in their lifestyle and morphology. Yet, both corals are known sources of briarane-type diterpenoids, a Scleralcyonacea-specific structural family of diterpenoids, numbering over 700 representatives and showing a wide array of potent biological properties.\u003csup\u003e29\u003c/sup\u003e\u003c/p\u003e \u003cp\u003eHere, we test the hypothesis that briarane diterpenoids are encoded by broadly distributed and ancient biosynthetic gene clusters, which is unprecedented across metazoans. We assembled and analyzed high-quality genome sequences from six key briarane producing octocorals from different taxonomic families, revealing a conserved gene cluster family. Using a combination of protein biochemistry and pathway engineering, we reconstituted the biosynthesis of cembrene B γ-lactone, which contains the γ-lactone structural feature distinctive of all briarane diterpenoids. These results provide compelling evidence that corals use BGCs to synthesize their abundant bioactive terpenoids, opening the possibility that other metazoans may have similarly evolved BGCs in their own specialized metabolism.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e"},{"header":"Results","content":"\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e \u003ch2\u003eChemical Validation of Briarane Production\u003c/h2\u003e \u003cp\u003eWe previously characterized homologous cembrene B synthases identified from the genomes of \u003cem\u003eR. muelleri\u003c/em\u003e and \u003cem\u003eE. caribaeorum.\u003c/em\u003e\u003csup\u003e27,28\u003c/sup\u003e Both encoding genes were colocalized with cytochromes P450 genes, which also shared homology with their counterparts on the two distinct genomes. The tantalizing possibility that coral animals evolved gene cluster families for the biosynthesis of structurally related specialized metabolites motivated us to evaluate whether diverse corals known to contain briarane diterpenoids also encode cembrene B synthases in their genomes. As such, we first set out to confirm literature reports of briarane producing octocoral families. Previously, Pennatuloidea, Ellisellidae, Coralliidae, Briareidae, and Erythropodiidae were established as briarane containing octocorals (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ea).\u003csup\u003e30\u0026ndash;36\u003c/sup\u003e As each of these five coral families belong to the Scleralcyonacea order, we further examined the remaining 14 families within this taxonomic order and found that five others reported terpenoid chemistry, yet not with briarane scaffolds (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eb, \u003cb\u003eSI Table\u0026nbsp;1\u003c/b\u003e) .\u003csup\u003e37\u003c/sup\u003e For the remaining nine families, no terpenoid literature was found.\u003c/p\u003e \u003cp\u003eWe thus collected coral tissue from representative species across five known briarane producing families, including \u003cem\u003eRenilla koellikeri\u003c/em\u003e (family Pennatuloidea), \u003cem\u003eStylatula elongata\u003c/em\u003e (family Pennatuloidea), \u003cem\u003eDichotella gemmacea\u003c/em\u003e (family Ellisellidae), the branching and encrusting morphologies of \u003cem\u003eBriareum asbestinum\u003c/em\u003e (family Briareidae), \u003cem\u003eE. caribaeorum\u003c/em\u003e (family Erythropodiidae), and \u003cem\u003eCorallium rubrum\u003c/em\u003e (family Coralliidae). \u003cem\u003eR. koellikeri\u003c/em\u003e collected, by SCUBA in San Diego, is not a reported source for briaranes, but is related to the Atlantic \u003cem\u003eR. muelleri\u003c/em\u003e that produces briaranes\u003csup\u003e38\u003c/sup\u003e and whose genome encodes cembrene B synthase as part of a putative BGC.\u003csup\u003e27,28\u003c/sup\u003e We thus analyzed the tissue of \u003cem\u003eR. koellikeri\u003c/em\u003e and isolated 11-hydroxyptilosarcenone (\u003cb\u003e6\u003c/b\u003e) whose chemical structure was confirmed by NMR spectroscopy (\u003cb\u003eNMR supplementary note, SI Table\u0026nbsp;2, SI\u003c/b\u003e Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e). This compound was previously identified from the sea pen \u003cem\u003ePtilosarcus gurneyi\u003c/em\u003e, but not yet from a \u003cem\u003eRenilla\u003c/em\u003e sea pansy.\u003csup\u003e39,40\u003c/sup\u003e\u003c/p\u003e \u003cp\u003eLimitations in accessible biomass for the other octocorals targeted in this study led us to use molecular networking analyses of crude extracts and semi-pure fractions to confirm briarane production. We used 11-hydroxyptilosarcenone (\u003cb\u003e6\u003c/b\u003e) as a standard in the molecular network to identify a cluster containing 54 nodes representing unique masses of briarane compounds (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e, \u003cb\u003eSI\u003c/b\u003e Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e, \u003cb\u003eSI Table\u0026nbsp;3\u003c/b\u003e). This analysis not only expanded the ptilosarcenone network from \u003cem\u003eR. koellikeri\u003c/em\u003e (Pennatuloidea) but also established related briaranes in \u003cem\u003eD. gemmacea\u003c/em\u003e (Ellisellidae), \u003cem\u003eE. caribeorum\u003c/em\u003e (Erythropodiidae), and \u003cem\u003eB. asbestinum\u003c/em\u003e (Briareidae). Two nodes corresponding to compounds in these corals had exact masses consistent with literature values: gemmacolide R (\u003cb\u003e7\u003c/b\u003e) from \u003cem\u003eD. gemmacea\u003c/em\u003e and erythrolide D (\u003cb\u003e8\u003c/b\u003e) from \u003cem\u003eE. caribaeorum\u003c/em\u003e (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ec, \u003cb\u003eSI Table\u0026nbsp;1\u003c/b\u003e).\u003csup\u003e41,42\u003c/sup\u003e However, we did not detect briarane nodes from extracts of \u003cem\u003eC. rubrum\u003c/em\u003e (Coralliidae) or \u003cem\u003eS. elongata\u003c/em\u003e (Pennatuloidea). While there is literature precedent of \u003cem\u003eS. elongata\u003c/em\u003e and corals in the family Coralliidae, although not specifically \u003cem\u003eC. rubrum\u003c/em\u003e, making briaranes, we included these corals in our study assuming that they should have the genomic capacity to produce briaranes even though the samples we analyzed were devoid at the time of collection.\u003csup\u003e33,34,43\u003c/sup\u003e\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec4\" class=\"Section2\"\u003e \u003ch2\u003eGenome assembly and syntenic analysis of the briarane gene cluster family in Scleralcyonacea corals\u003c/h2\u003e \u003cp\u003eAt the start of this study, a few briarane containing coral draft genomes with very low continuity (contig N50) were available: \u003cem\u003eE. caribaeorum\u003c/em\u003e (Contig N50 2 kb)\u003csup\u003e28\u003c/sup\u003e, \u003cem\u003eR. muelleri\u003c/em\u003e (N50 70.5 kb)\u003csup\u003e44\u003c/sup\u003e, and \u003cem\u003ePteroeides caledonicum\u003c/em\u003e (N50 4 kb) (\u003cb\u003eSI Table\u0026nbsp;4\u003c/b\u003e for assembly accessions and statistics).To capture intact BGCs, we proceeded to generate higher quality assemblies of corals across four briarane producing families, obtaining long read sequencing data for \u003cem\u003eR. koellikeri\u003c/em\u003e (family Pennatuloidea), \u003cem\u003eS. elongata\u003c/em\u003e (family Pennatuloidea), \u003cem\u003eD. gemmacea\u003c/em\u003e (family Ellisellidae), the branching and encrusting morphologies of \u003cem\u003eB. asbestinum\u003c/em\u003e (family Briareidae), and \u003cem\u003eE. caribaeorum\u003c/em\u003e (family Erythropodiidae). We assembled 26 to 78 gigabase pairs (Gb) of sequencing reads into genomes ranging in size from 231 megabase pairs (Mb) (\u003cem\u003eR. koellikeri\u003c/em\u003e) to 1,237 Mb (\u003cem\u003eB. asbestinum\u003c/em\u003e) (\u003cb\u003eTable\u0026nbsp;1\u003c/b\u003e). To assess the completeness of our assemblies, we estimated genome sizes with K-mer based analyses using Illumina short reads (\u003cb\u003eSI\u003c/b\u003e Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e), producing genomes ranging from 144 Mb (\u003cem\u003eR. koellikeri\u003c/em\u003e) to 1,155 Mb (\u003cem\u003eB. asbestinum\u003c/em\u003e) (\u003cb\u003eTable\u0026nbsp;1\u003c/b\u003e). Illumina short read coverage was too shallow for \u003cem\u003eE. caribaeorum\u003c/em\u003e to allow for K-mer based analysis (\u003cb\u003eSI Table\u0026nbsp;5\u003c/b\u003e for all sequencing data used). Additionally, we included our recently reported chromosomally resolved \u003cem\u003eC. rubrum\u003c/em\u003e (family Coralliidae) genome at 475 Mb (20 chromosomes, 951 scaffolds)\u003csup\u003e45,46\u003c/sup\u003e in our analysis to enable BGC exploration across all five known briarane producing coral families. The genome assemblies are now some of the most contiguous for octocorals with contig N50 values above 196 kb (\u003cb\u003eTable\u0026nbsp;1\u003c/b\u003e,Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003ea, \u003cb\u003eSI Table\u0026nbsp;6\u003c/b\u003e). The chromosomal resolution assembly of \u003cem\u003eC. rubrum\u003c/em\u003e had a scaffold N50 of 16,290 kb, on par with the current best scaffolded octocoral assembly from \u003cem\u003eXenia\u003c/em\u003e sp. with a scaffold N50 of 14,832 kb.\u003csup\u003e47\u003c/sup\u003e Benchmarking universal single-copy ortholog (BUSCO) completeness scores ranged from 86.4\u0026ndash;88% (\u003cb\u003eTable\u0026nbsp;1\u003c/b\u003e). BUSCO scores are highly dependent on the database used for the analysis and the low scores could reflect the use of the metazoan database, which was the closest available. The lower-than-expected completeness is on par with the BUSCO for the \u003cem\u003eXenia\u003c/em\u003e assembly, which was 88.1% complete (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003ea).\u003csup\u003e47\u003c/sup\u003e We further assembled transcriptomes to guide our genomic analyses by identifying actively transcribed genes and determining intron bounds (\u003cb\u003eSI Table\u0026nbsp;6\u003c/b\u003e).\u003c/p\u003e \u003cp\u003eUsing the putative briarane BGC gene sequences from \u003cem\u003eR. muelleri\u003c/em\u003e and \u003cem\u003eE. caribeorum\u003c/em\u003e as queries, we analyzed the five newly established genomes of \u003cem\u003eD. gemmacea, B. asbestinum, R. koellikeri, C. rubrum\u003c/em\u003e, and \u003cem\u003eS. elongata\u003c/em\u003e and the public \u003cem\u003eP. caledonicum\u003c/em\u003e for homologous genes. In total, the genomes from these eight octocoral species span five of the 20 octocoral families in the order Scleralcyonacea. To determine if non-briarane producing octocoral families also have syntenic BGCs, we also evaluated the publicly available \u003cem\u003eHeliopora coerulea\u003c/em\u003e (family Helioporidae) genome.\u003csup\u003e48\u003c/sup\u003e We identified homologs of the cembrene B terpene cyclase (cbTC) in all eight octocorals shown or suspected to produce briaranes, whereas the \u003cem\u003eH. coerulea\u003c/em\u003e genome did not contain a homologous sequence, although it did encode nine other TCs. In the eight octocoral species containing a homologous cbTC gene, seven contained iterations of the briarane BGC (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eb). Given that the \u003cem\u003eP. caledonicum\u003c/em\u003e genome had poor contiguity, we were not surprised to find the cbTC homolog on a 3 kb contig and expect that it too is part of a BGC.\u003c/p\u003e \u003cp\u003eWe observed one conserved short-chain dehydrogenase (cbSDH) and up to three cytochrome P450s (cbCYPa, cbCYPb, and cbCYPc) per briarane BGC (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eb, \u003cb\u003eSI\u003c/b\u003e Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e). In the eight octocorals with briarane BGCs, four maintained the full five-gene briarane BGC spanning 22 kb in \u003cem\u003eR. muelleri\u003c/em\u003e up to 40 kb in \u003cem\u003eC. rubrum\u003c/em\u003e, albeit with different gene organization. Notably, in the cases of \u003cem\u003eS. elongata, E. caribeorum\u003c/em\u003e, and \u003cem\u003eB. asbestinum\u003c/em\u003e, the orthologous genes were split between two BGC loci (\u003cb\u003eSI Table\u0026nbsp;7\u003c/b\u003e). For \u003cem\u003eS. elongata\u003c/em\u003e and \u003cem\u003eB. asbestinum\u003c/em\u003e, their respective cbTC and cbCYPc genes reside on one contig, while cbCYPa, cbCYPb, and at least two cbSDHs are on a second contig. In the case of \u003cem\u003eE. caribeorum\u003c/em\u003e, cbCYPc resides outside the BGC contig. Also, in the case of \u003cem\u003eB. asbestinum\u003c/em\u003e, the orthologous genes have much longer gaps resulting in a significantly larger BGCs spanning\u0026thinsp;~\u0026thinsp;150 kb. \u003cem\u003eS. elongata\u003c/em\u003e and \u003cem\u003eB. asbestinum\u003c/em\u003e further have multiple copies of the SDH gene at three (85%, 89% and 90% sequence similarity) and two (46% sequence similarity) copies, respectively (\u003cb\u003eSI\u003c/b\u003e Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e). The duplications of the SDHs and gene rearrangements suggest lineage specific independent evolution across octocoral species.\u003c/p\u003e \u003cp\u003eSince we identified the same set of homologous cytochromes P450s in a wide range of octocoral species, we conducted a phylogenetic analysis of these sequences in the background of broader cnidarian cytochromes P450s to shed light on their evolutionary origin. Our analysis included CYPs from four non-octocoral cnidarians as a reference set: \u003cem\u003eHydra vulgaris\u003c/em\u003e (24 CYPs), \u003cem\u003eAcropora digitifera\u003c/em\u003e (24 CYPs), \u003cem\u003eAurelia aurita\u003c/em\u003e (37 CYPs), and \u003cem\u003eNematostella vectensis\u003c/em\u003e (70 CYPs).\u003csup\u003e49\u003c/sup\u003e From octocorals we included CYPs from six briarane producing species \u003cem\u003eD. gemmacea, B. asbestinum, R. koellikeri, S. elongata, C. rubrum\u003c/em\u003e, and \u003cem\u003eE. caribaeorum\u003c/em\u003e, as well as the non-briarane producers \u003cem\u003eXenia\u003c/em\u003e sp. and \u003cem\u003eH. coerulea\u003c/em\u003e, yielding 318 octocoral CYP sequences at ~\u0026thinsp;40 per genome (\u003cb\u003eSI Table\u0026nbsp;8\u003c/b\u003e). We combined the annotated reference cnidarian CYPs with the annotated octocoral CYPs to build a phylogenetic tree (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003ec, \u003cb\u003eSI Fig.\u0026nbsp;6\u003c/b\u003e). The CYPs from the four reference cnidarians were previously associated with accepted CYP clans and families,\u003csup\u003e49\u003c/sup\u003e thereby enabling us to examine the evolutionary context of the octocoral BGC CYPs. Notably, some octocoral CYPs fell into the clades associated with the cnidarian CYP clans, suggesting that they fulfill primary, widespread functions generally present in cnidarians. Others, including the BGC-associated CYPs, formed new monophyletic clades that were octocoral specific. The three briarane CYPs, cbCYPa-c, form monophyletic clades consisting of cbCYPa, cbCYPb, and cbCYPc homologs from different species (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003ec). While cbCYPb and cbCYPc share a common ancestor that may suggest related enzymatic functions, cbCYPa is distantly related and instead lies closer to clan 46 with predicted functions associated with steroid biosynthesis (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003ec). Interestingly, \u003cem\u003eH. coerulea\u003c/em\u003e, while also from the order Scleralcyonacea has CYPs closely related to the cbCYPb and cbCYPc clades, but these genes are not co-localized hinting at the loss of the briarane BGC in this lineage.\u003c/p\u003e \u003cp\u003eTo test whether our phylogenetic approach could be used for genome mining in octocorals with lower quality assemblies, we next included CYP sequences from the poor contiguity genome assembly of \u003cem\u003eP. caledonicum.\u003c/em\u003e We found corresponding sequences present in two of the three clades, thereby identifying cbCYPa and cbCYPb orthologs in the genome (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003ed).\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"No\" id=\"Taba\" border=\"1\"\u003e \u003ccolgroup cols=\"7\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c6\" colnum=\"6\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c7\" colnum=\"7\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/th\u003e \u003cth align=\"left\" colspan=\"6\" nameend=\"c7\" namest=\"c2\"\u003e \u003cp\u003eTable\u0026nbsp;1. Assembly statistics and BUSCO assessment of five octocoral genomes\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eOrganism\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eContigs (#)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eAssembly size (Mb)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eK-mer Based Genome size (Mb)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eContig N50 (kb)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eLongest (kb)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003eBUSCO (%C/ %D)\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eBriareum asbestinum\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e22,526\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e1,254\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e1,155\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e332\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e4,126\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e87.3/4.6\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eDichotella gemmacea\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e8,171\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e550\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e401\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e196\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e2,036\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e87.5/15.6\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eRenilla koellikeri\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e557\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e231\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e144\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e4,895\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e22,942\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e86.9/2.6\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eStylatula elongata\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e392\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e360\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e269\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e5,139\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e39,131\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e88.0/1.0\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eErythropodium caribeorum\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e8,289\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e300\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eNA\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e233\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e1,576\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e86.4/0.6\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eCorallium rubrum\u003c/em\u003e\u003csup\u003e45\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e20 chromosomes, 951 scaffolds\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e545\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e540\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e1,600\u003c/p\u003e \u003cp\u003e(scaffold N50 18,521kb)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e115,000\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e88.5/1.2\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec5\" class=\"Section2\"\u003e \u003ch2\u003eBiochemical validation of the briarane BGC\u003c/h2\u003e \u003cp\u003eTo provide support that the briarane BGC is involved in the early stages of briarane biosynthesis, we set out to provide biochemical validation of the clustered biosynthesis genes (\u003cb\u003eSI Table\u0026nbsp;9\u003c/b\u003e). We first tested whether all the homologous cbTC terpene synthases functioned as cembrene B cyclases as we originally established in \u003cem\u003eR. muelleri\u003c/em\u003e and \u003cem\u003eE. caribaeorum\u003c/em\u003e.\u003csup\u003e27,28\u003c/sup\u003e We expressed each, 6 in total, in \u003cem\u003eSaccharomyces cerevisiae\u003c/em\u003e and evaluated their production profiles by GC-MS analysis (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e, \u003cb\u003eSI Fig.\u0026nbsp;7\u003c/b\u003e). In all cases, we confirmed that GGPP (\u003cb\u003e9\u003c/b\u003e) was cyclized into cembrene B (\u003cb\u003e10\u003c/b\u003e) as the sole product.\u003c/p\u003e \u003cp\u003eFor the remaining enzyme assays, we tested the functions of genes from the \u003cem\u003eR. koellikeri\u003c/em\u003e BGC and additional representative homologs identified by synteny. All tested enzyme homologs showed identical reactivity. We evaluated the functions of the three co-localized CYP subtypes in pairwise fashion by co-expression with cognate cbTC genes using previously established methods (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e).\u003csup\u003e50\u0026ndash;52\u003c/sup\u003e We consistently measured the production of 19-hydroxycembrene B (\u003cb\u003e11\u003c/b\u003e) in yeast upon co-expression with cbCYPb homologs. Its formation proceeded in a regiospecific manner with no other detectable monooxygenated products (\u003cb\u003eSI Fig.\u0026nbsp;8\u003c/b\u003e). Scale-up fermentation and chromatographic purification provided \u003cb\u003e11\u003c/b\u003e in pure form, which we characterized by NMR spectroscopic analyses (\u003cb\u003eSI Tables\u0026nbsp;10 and 11, NMR supplementary note, SI Fig.\u0026nbsp;9\u003c/b\u003e). When we instead co-expressed cbCYPc homologs with cbTC genes, we observed formation of 17-hydroxycembrene (\u003cb\u003eS1\u003c/b\u003e) which could be detected by LCMS (\u003cb\u003eSI Fig.\u0026nbsp;10\u003c/b\u003e) and characterized by NMR (\u003cb\u003eSI Tables\u0026nbsp;11 and 12, NMR supplementary note, SI Fig.\u0026nbsp;9\u003c/b\u003e). Because cbCYPb and cbCYPc were found to oxidize complimentary positions on \u003cb\u003e10\u003c/b\u003e with respect to the briarane pathway, we tested the cbCYPb product \u003cb\u003e11\u003c/b\u003e by feeding it to yeast expressing cbCYPc homologs (\u003cb\u003eSI Fig.\u0026nbsp;11\u003c/b\u003e). We clearly observed a single oxidized product that upon purification and combined spectroscopic (\u003cb\u003eSI Tables\u0026nbsp;10 and 11, NMR supplementary note, SI Fig.\u0026nbsp;9\u003c/b\u003e) and x-ray crystallographic characterization (\u003cb\u003eSI Fig.\u0026nbsp;12, SI Table\u0026nbsp;12\u003c/b\u003e) was established as (7\u003cem\u003eS\u003c/em\u003e)-7,19-dihydroxycembrene B (\u003cb\u003e12\u003c/b\u003e). Likewise, when yeast harboring cbCYPb genes were fed \u003cb\u003eS1\u003c/b\u003e, we similarly obtained \u003cb\u003e12\u003c/b\u003e, albeit with lower rates of conversion (see Methods). The cooperative and related functions of the cbCYPb and cbCYPc enzymes correlates well with their shared evolution (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003ec). Despite numerous attempts with different substrate and gene pairs, we have yet to observe a function for the third cytochrome P450, cbCYPa.\u003c/p\u003e \u003cp\u003eAs a variety of cembrene double bond isomers are produced by corals, we considered the possibility that the described CYP450 enzymes could act non-specifically on other cembrene isomers. We thus co-expressed cbCYPb homologs with terpene synthases that selectively produce either cembrene A or C. However, no oxidized products were produced (\u003cb\u003eSI Fig.\u0026nbsp;13\u003c/b\u003e), showing high specificity of cbCYPb for cembrene B.\u003c/p\u003e \u003cp\u003eWe next evaluated the co-localized short-chain dehydrogenase cbSDH homologs for their possible involvement in further transforming the diol \u003cb\u003e12\u003c/b\u003e. We expressed cbSDHs in \u003cem\u003eEscherichia coli\u003c/em\u003e to conduct a series of in vitro incubation experiments. Incubation of purified cbSDH enzymes with NADP and \u003cb\u003e11\u003c/b\u003e produced aldehyde (\u003cb\u003eS2\u003c/b\u003e), while incubation with the diol \u003cb\u003e12\u003c/b\u003e produced cembrene B g-lactone (\u003cb\u003e13\u003c/b\u003e), a structure feature diagnostic of all briarane diterpenoids (\u003cb\u003eSI Figs.\u0026nbsp;14 and 15\u003c/b\u003e). Two consecutive dehydrogenation steps were detected by LCMS chromatography via the formation of a \u003cem\u003em/z\u003c/em\u003e 301.2183 product ion ([M\u0026thinsp;+\u0026thinsp;H]\u003csup\u003e+\u003c/sup\u003e, calcd. \u003cem\u003em/z\u003c/em\u003e 301.2163 for C\u003csub\u003e20\u003c/sub\u003eH\u003csub\u003e29\u003c/sub\u003eO\u003csub\u003e2\u003c/sub\u003e\u003csup\u003e+\u003c/sup\u003e). We scaled-up the in vitro reaction of \u003cb\u003e12\u003c/b\u003e and NADP with the cbSDH to provide lactone \u003cb\u003e13\u003c/b\u003e, which was purified and characterized spectroscopically. Scale up experiments were conducted with representative characterized genes to unambiguously determine all structures. The structure was supported by NMR (\u003cb\u003eSI Tables\u0026nbsp;11 and 12, NMR supplementary note, SI Fig.\u0026nbsp;9\u003c/b\u003e) and the diagnostic γ-lactone carbonyl stretch in the IR spectrum (\u003cem\u003ev\u003c/em\u003e\u003csub\u003emax\u003c/sub\u003e1757 cm\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e). It is possible that \u003cb\u003e13\u003c/b\u003e forms from \u003cb\u003e12\u003c/b\u003e via a five-membered hemiacetal intermediate (\u003cb\u003eSI Fig.\u0026nbsp;16\u003c/b\u003e).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e"},{"header":"Discussion","content":"\u003cp\u003eWe discovered a broadly conserved BGC family in corals and biochemically validated its function in the early biosynthetic steps toward briarane diterpenoids. The distinctive family of bioactive molecules numbers over 700 characterized members that are exclusively found in corals. Our work established the biosynthetic logic of briarane\u0026rsquo;s diagnostic γ-lactone core structure by four of the co-clustered gene products, cbTC, cbCYPb, cbCYPc, and cbSDH, confirming the clustered function of the briarane locus. The enzymatic transformations begin with cembrene B (\u003cb\u003e10\u003c/b\u003e) formation by cbTC, then proceed with subsequent 19- and 7-hydroxylations carried out by cbCYPb and cbCYPc, respectively. Finally, in a unique mechanism of lactone biosynthesis, cbSDH transforms diol \u003cb\u003e12\u003c/b\u003e directly into lactone \u003cb\u003e13\u003c/b\u003e (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eA), likely via aldehyde formation at the 19-position, hemiacetal formation and a second dehydrogenation to produce the diagnostic lactone (\u003cb\u003eSI Fig.\u0026nbsp;16\u003c/b\u003e). We have not yet characterized the enzyme(s) catalyzing bond formation between carbons 1 and 10.\u003csup\u003e29\u003c/sup\u003e A noncanonical class 2 TC acting on lactone \u003cb\u003e13\u003c/b\u003e could introduce this feature to produce \u003cb\u003e14\u003c/b\u003e (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eA).\u003csup\u003e53\u003c/sup\u003e Alternatively, a Δ\u003csup\u003e1,14\u003c/sup\u003e epoxide intermediate, potentially installed by a P450, might precede a transannular cyclization reaction to \u003cb\u003e18\u003c/b\u003e (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eB).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eIn animals, previous reports indicate that genes can be found in clusters when they result from successive tandem duplications or HGT. However, such factors are not at play here in the octocorals. Homologs of the briarane TC and CYP genes are ancient, found throughout octocorallia (\u0026gt;\u0026thinsp;500 Mya). The briarane BGC itself is likely evolutionarily ancient and evolved soon after the division of the octocoral orders Scleralcyonacea and Malacalcyonacea. It is only found in briarane-containing corals from Sclearcyonacea and has similar synteny in distantly related corals such as \u003cem\u003eRenilla koellikeri\u003c/em\u003e (Family Pennatoloidea) and \u003cem\u003eErythropodium caribeorum\u003c/em\u003e (Family Erythropodiidae) (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eC).\u003csup\u003e54\u003c/sup\u003e A key evolutionary event in the formation of the briarane BGC was the duplication of the related cbCYPb and -c genes early in the evolution of Scleralcyonacea. Supporting this timing, cbCYPb / -c pair is found across the Scleralcyonacea coral families Pennatuloidea, Ellisellidae, Coralliidae, Briareidae, and Erythropodiidae, but not elsewhere in sequenced octocorals. Subsequently, we speculate that cbTC, cbSDH, and cbCYPa were recruited and incorporated adjacent to the duplicated cbCYPb\u0026ndash;c gene pair, also before speciation in the Scleralcyonacea, via chromosomal rearrangements to form the BGC \u0026mdash; a phenomenon not previously reported in metazoans. Following the formation of the briarane BGC, all or parts of it were lost in a group of octocorals, as no briaranes have been reported from several Scleralcyonacea families (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eC). Supporting the gene loss scenario, \u003cem\u003eH. coerulea\u003c/em\u003e (family Helioporidae), does not harbor a cembrene B synthase in its genome but retains cbCYPb and -c relatives, suggesting that these genes were once part of a briarane BGC. More complex evolutionary scenarios are also observable in sequenced genomes, such as the BGCs from \u003cem\u003eS. elongata\u003c/em\u003e and \u003cem\u003eB. asbestinum\u003c/em\u003e, which contain additional SDH gene copies. Their appearance coincides with a split of the gene cluster into two parts. The observation of metazoan BGC formation and decay suggests potential mechanisms of chromosomal evolution not previously observed in animals.\u003c/p\u003e \u003cp\u003eIn contrast to a well preserved central briarane BGC, genes for specialized, late-stage tailoring processes are not colocalized. The core briarane skeleton is modified in many ways to create the \u0026gt;\u0026thinsp;700 known briarane derivatives (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eA). For example, diterpenoids such as 11-hydroxyptilosarcenone (\u003cb\u003e6\u003c/b\u003e), gemmacolide R (\u003cb\u003e7\u003c/b\u003e), bathyptilone A (\u003cb\u003e15\u003c/b\u003e), and erythrolide B (\u003cb\u003e16\u003c/b\u003e), found in different coral species, have been modified by diverse enzymatic reactions including oxidation, chlorination, and acylation, by genes found in different chromosomal locations. This may represent a general rule in metazoan BGCs, in which the required core genes are found in the BGC, while those for variable and species-specific tailoring steps are distributed through the genomes.\u003c/p\u003e \u003cp\u003eWhile the innate functions of coral terpenoids are largely unknown outside of a few examples having feeding deterrent properties,\u003csup\u003e55\u003c/sup\u003e we suspect that the conservation of their BGCs indicates central, unrealized biological roles pertaining to some of the more common coral terpenoid classes.\u003csup\u003e56\u003c/sup\u003e With the characterization of a coral terpenoid BGC, future genetic experiments are positioned to address their function once genome editing methods have been established in octocorals as recently accomplished in the hexacoral \u003cem\u003eAcropora millepora.\u003c/em\u003e\u003csup\u003e57\u003c/sup\u003e\u003c/p\u003e \u003cp\u003eThe stage is set for genome mining and bioengineering studies to delineate biosynthetic cascade reactions to rare, bioactive coral molecules for their biosustainable production and application as therapeutic compounds. The briarane gene cluster family is likely the first of many to be described. We previously reported that some corals contain 70 kb and 500 kb putative BGCs with upward of 13 terpenoid biosynthesis genes.\u003csup\u003e27,28\u003c/sup\u003e With over 4,000 terpenoids characterized from corals (determined by searching the CMNPD database),\u003csup\u003e58\u003c/sup\u003e gene clustering is likely common in this ancient phylum. Looking ahead, we have yet to see if the active assembly of biosynthetic gene clusters by evolutionary forces is an exceedingly rare event that is limited to octocorals or if this is a more common phenomenon in animals and eukaryotes in general that has been overlooked so far due to limited genomic information in the eukaryotic kingdom.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eAcknowledgments\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eWe thank P. Zerofski (Scripps Institution of Oceanography) for collecting \u003cem\u003eR. koellikeri\u003c/em\u003e and \u003cem\u003eS. elongata\u003c/em\u003e, J. Sprung (Two Little Fishes Inc) for the \u003cem\u003eD. gemmacea\u0026nbsp;\u003c/em\u003etissue, and J. Garrabou (CIIMAR) for the \u003cem\u003eC. rubrum\u003c/em\u003e tissue. For her help in species identification, we would like to recognize C. S. McFadden (Harvey Mudd College). For coral images, we thank C. S. McFadden- \u003cem\u003eD. gemmacea\u003c/em\u003e, F. Zuberer (CNRS)- \u003cem\u003eC. rubrum\u003c/em\u003e, P. Webster @underwaterpat- \u003cem\u003eS. elongata\u003c/em\u003e, and J. Simpson (Herbert Wertheim UF Scripps Institute for Biomedical Innovation \u0026amp; Technology)- \u003cem\u003eB. asbestinum\u0026nbsp;\u003c/em\u003eand \u003cem\u003eE. caribaeorum\u003c/em\u003e. We thank T. Damiani (IOCB Prague) and A. Rodriguez (UC San Diego) for their assistance with feature based molecular networking. Chemical isolation and elucidation of the \u003cem\u003eR. koellikeri\u003c/em\u003e briarane diterpene was aided by A. Bogdanov (Scripps Institution of Oceanography). This work was supported by the National Institutes of Health (R01-GM146224 to B.S.M., R35-GM148283 to E.W.S, and K99-GM148783 to P.D.S.), a Margaret A. Davidson Graduate Fellowship to N.E.G. (NERRS NA22NOS4200050), an NSF Graduate Research Fellowship to M.L.M., Tang Genomics Fund to T.P.M,, and Fundaço para a Ciência e a Tecnologia funds (UIDB/04423/2020, UIDP/04423/2020 and 2021.00855.CEECIND) to J.B.L. X-ray diffraction research reported in this publication was supported by the Office of the Director, NIH under award S10-OD030326. NMR data obtained at UC San Diego (the Scripps Institution of Oceanography and Biomolecular NMR Facilities) was assisted by B. Duggan. Additional NMR data was obtained at the University of Utah Health Sciences NMR Core. GCMS data was obtained at the University of Utah Health Sciences Proteomics Core by Q. Pearce.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthor Contributions:\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eConceptualization: PDS, IB, EWS, BSM\u003c/p\u003e\n\u003cp\u003eData creation: NEG, PDS, MLM, IB\u003c/p\u003e\n\u003cp\u003eFormal analysis and validation: NEG, PDS, IB, JGG, TA\u003c/p\u003e\n\u003cp\u003eFunding acquisition: PDS, JBL, TPM, EWS, BSM\u003c/p\u003e\n\u003cp\u003eInvestigation and methodology: NEG, PDS, MLM, IB\u003c/p\u003e\n\u003cp\u003eProject administration: PDS, IB, EWS, BSM\u003c/p\u003e\n\u003cp\u003eResources: JBL, TPM, EWS, BSM\u003c/p\u003e\n\u003cp\u003eWriting \u0026ndash; original draft: NEG, PDS, IB, BSM\u003c/p\u003e\n\u003cp\u003eWriting \u0026ndash; review and editing: All authors\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCompeting interests:\u003c/strong\u003e\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eAuthors declare that they have no competing interests.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eData and materials availability:\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe feature based molecular network parameters used and raw data are available through the MASSIVE dataset ID: MSV000094792.\u003c/p\u003e\n\u003cp\u003eAccessions numbers for the genomes generated in this paper are as follows; \u003cem\u003eR. koellikeri\u0026nbsp;\u003c/em\u003e(SAMN40621396)\u003cem\u003e, B. asbestinum\u0026nbsp;\u003c/em\u003e(SAMN40621398)\u003cem\u003e, D. gemmacea\u0026nbsp;\u003c/em\u003e(SAMN40621399)\u003cem\u003e,\u003c/em\u003e and\u003cem\u003e\u0026nbsp;S. elongata\u0026nbsp;\u003c/em\u003e(SAMN40621397).\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n\u003cli\u003eMedema, M.H., de Rond, T., and Moore, B.S. (2021). Mining genomes to illuminate the specialized chemistry of life. Nat. Rev. Genet. \u003cem\u003e22\u003c/em\u003e, 553\u0026ndash;571. https://doi.org/10.1038/s41576-021-00363-7.\u003c/li\u003e\n\u003cli\u003eTorres, J.P., and Schmidt, E.W. (2019). The biosynthetic diversity of the animal world. J. Biol. Chem. \u003cem\u003e294\u003c/em\u003e, 17684\u0026ndash;17692. https://doi.org/10.1074/jbc.REV119.006130.\u003c/li\u003e\n\u003cli\u003ePoth, D., Wollenberg, K.C., Vences, M., and Schulz, S. (2012). Volatile Amphibian Pheromones: Macrolides from Mantellid Frogs from Madagascar. Angew. Chem. Int. Ed. \u003cem\u003e51\u003c/em\u003e, 2187\u0026ndash;2190. https://doi.org/10.1002/anie.201106592.\u003c/li\u003e\n\u003cli\u003eChinta, S.P., Goller, S., Lux, J., Funke, S., Uhl, G., and Schulz, S. (2010). The Sex Pheromone of the Wasp Spider Argiope bruennichi. Angew. Chem. Int. Ed. \u003cem\u003e49\u003c/em\u003e, 2033\u0026ndash;2036. https://doi.org/10.1002/anie.200906311.\u003c/li\u003e\n\u003cli\u003eMasello, J.F., Lubjuhn, T., and Quillfeldt, P. (2008). Is the structural and psittacofulvin-based coloration of wild burrowing parrots Cyanoliseus patagonus condition dependent? J. Avian Biol. \u003cem\u003e39\u003c/em\u003e, 653\u0026ndash;662. https://doi.org/10.1111/j.1600-048X.2008.04417.x.\u003c/li\u003e\n\u003cli\u003eAnderson, H.A., Mathieson, J.W., and Thomson, R.H. (1969). Distribution of spinochrome pigments in echinoids. Comp. Biochem. Physiol. \u003cem\u003e28\u003c/em\u003e, 333\u0026ndash;345. https://doi.org/10.1016/0010-406X(69)91347-4.\u003c/li\u003e\n\u003cli\u003eButcher, R.A., Fujita, M., Schroeder, F.C., and Clardy, J. (2007). Small-molecule pheromones that control dauer development in Caenorhabditis elegans. Nat. Chem. Biol. \u003cem\u003e3\u003c/em\u003e, 420\u0026ndash;422. https://doi.org/10.1038/nchembio.2007.3.\u003c/li\u003e\n\u003cli\u003eDettner, K., and Schwinger, G. (1982). Defensive secretions of three oxytelinae rove beetles (Coleoptera: Staphylinidae). J. Chem. Ecol. \u003cem\u003e8\u003c/em\u003e, 1411\u0026ndash;1420. https://doi.org/10.1007/BF01403104.\u003c/li\u003e\n\u003cli\u003eOlivera, B.M., Gray, W.R., Zeikus, R., McIntosh, J.M., Varga, J., Rivier, J., de Santos, V., and Cruz, L.J. (1985). Peptide Neurotoxins from Fish-Hunting Cone Snails. Science \u003cem\u003e230\u003c/em\u003e, 1338\u0026ndash;1343. https://doi.org/10.1126/science.4071055.\u003c/li\u003e\n\u003cli\u003eChevrette, M.G., Guti\u0026eacute;rrez-Garc\u0026iacute;a, K., Selem-Mojica, N., Aguilar-Mart\u0026iacute;nez, C., Ya\u0026ntilde;ez-Olvera, A., Ramos-Aboites, H.E., Hoskisson, P.A., and Barona-G\u0026oacute;mez, F. (2020). Evolutionary dynamics of natural product biosynthesis in bacteria. Nat. Prod. Rep. \u003cem\u003e37\u003c/em\u003e, 566\u0026ndash;599. https://doi.org/10.1039/C9NP00048H.\u003c/li\u003e\n\u003cli\u003eRokas, A., Mead, M.E., Steenwyk, J.L., Raja, H.A., and Oberlies, N.H. (2020). Biosynthetic gene clusters and the evolution of fungal chemodiversity. Nat. Prod. Rep. \u003cem\u003e37\u003c/em\u003e, 868\u0026ndash;878. https://doi.org/10.1039/C9NP00045C.\u003c/li\u003e\n\u003cli\u003eN\u0026uuml;tzmann, H.-W., Scazzocchio, C., and Osbourn, A. (2018). Metabolic Gene Clusters in Eukaryotes. Annu. Rev. Genet. \u003cem\u003e52\u003c/em\u003e, 159\u0026ndash;183. https://doi.org/10.1146/annurev-genet-120417-031237.\u003c/li\u003e\n\u003cli\u003eSteele, T.S., Brunson, J.K., Maeno, Y., Terada, R., Allen, A.E., Yotsu-Yamashita, M., Chekan, J.R., and Moore, B.S. (2022). Domoic acid biosynthesis in the red alga Chondria armata suggests a complex evolutionary history for toxin production. Proc. Natl. Acad. Sci. \u003cem\u003e119\u003c/em\u003e, e2117407119. https://doi.org/10.1073/pnas.2117407119.\u003c/li\u003e\n\u003cli\u003eSmit, S.J., and Lichman, B.R. (2022). Plant biosynthetic gene clusters in the context of metabolic evolution. Nat. Prod. Rep. \u003cem\u003e39\u003c/em\u003e, 1465\u0026ndash;1482. https://doi.org/10.1039/D2NP00005A.\u003c/li\u003e\n\u003cli\u003eRyngajłło, M., Boruta, T., and Bizukojć, M. (2021). Complete genome sequence of lovastatin producer Aspergillus terreus ATCC 20542 and evaluation of genomic diversity among A. terreus strains. Appl. Microbiol. Biotechnol. \u003cem\u003e105\u003c/em\u003e, 1615\u0026ndash;1627. https://doi.org/10.1007/s00253-021-11133-0.\u003c/li\u003e\n\u003cli\u003eItkin, M., Heinig, U., Tzfadia, O., Bhide, A.J., Shinde, B., Cardenas, P.D., Bocobza, S.E., Unger, T., Malitsky, S., Finkers, R., et al. (2013). Biosynthesis of Antinutritional Alkaloids in Solanaceous Crops Is Mediated by Clustered Genes. Science \u003cem\u003e341\u003c/em\u003e, 175\u0026ndash;179. https://doi.org/10.1126/science.1240230.\u003c/li\u003e\n\u003cli\u003eN\u0026uuml;tzmann, H.-W., Huang, A., and Osbourn, A. (2016). Plant metabolic clusters \u0026ndash; from genetics to genomics. New Phytol. \u003cem\u003e211\u003c/em\u003e, 771\u0026ndash;789. https://doi.org/10.1111/nph.13981.\u003c/li\u003e\n\u003cli\u003eOsborn, A.R., Almabruk, K.H., Holzwarth, G., Asamizu, S., LaDu, J., Kean, K.M., Karplus, P.A., Tanguay, R.L., Bakalinsky, A.T., and Mahmud, T. (2015). De novo synthesis of a sunscreen compound in vertebrates. eLife \u003cem\u003e4\u003c/em\u003e, e05919. https://doi.org/10.7554/eLife.05919.\u003c/li\u003e\n\u003cli\u003eRoelofs, D., Timmermans, M.J.T.N., Hensbergen, P., van Leeuwen, H., Koopman, J., Faddeeva, A., Suring, W., de Boer, T.E., Mari\u0026euml;n, J., Boer, R., et al. (2013). A Functional Isopenicillin N Synthase in an Animal Genome. Mol. Biol. Evol. \u003cem\u003e30\u003c/em\u003e, 541\u0026ndash;548. https://doi.org/10.1093/molbev/mss269.\u003c/li\u003e\n\u003cli\u003eFaddeeva-Vakhrusheva, A., Kraaijeveld, K., Derks, M.F.L., Anvar, S.Y., Agamennone, V., Suring, W., Kampfraath, A.A., Ellers, J., Le Ngoc, G., van Gestel, C.A.M., et al. (2017). Coping with living in the soil: the genome of the parthenogenetic springtail Folsomia candida. BMC Genomics \u003cem\u003e18\u003c/em\u003e, 493. https://doi.org/10.1186/s12864-017-3852-x.\u003c/li\u003e\n\u003cli\u003eFeng, L., Gordon, M.T., Liu, Y., Basso, K.B., and Butcher, R.A. (2021). Mapping the biosynthetic pathway of a hybrid polyketide-nonribosomal peptide in a metazoan. Nat. Commun. \u003cem\u003e12\u003c/em\u003e, 4912. https://doi.org/10.1038/s41467-021-24682-9.\u003c/li\u003e\n\u003cli\u003ede Rosa, R., Grenier, J.K., Andreeva, T., Cook, C.E., Adoutte, A., Akam, M., Carroll, S.B., and Balavoine, G. (1999). Hox genes in brachiopods and priapulids and protostome evolution. Nature \u003cem\u003e399\u003c/em\u003e, 772\u0026ndash;776. https://doi.org/10.1038/21631.\u003c/li\u003e\n\u003cli\u003eStorz, J.F., Opazo, J.C., and Hoffmann, F.G. (2013). Gene duplication, genome duplication, and the functional diversification of vertebrate globins. Mol. Phylogenet. Evol. \u003cem\u003e66\u003c/em\u003e, 469\u0026ndash;478. https://doi.org/10.1016/j.ympev.2012.07.013.\u003c/li\u003e\n\u003cli\u003eKitchen, S.A., Naragon, T.H., Br\u0026uuml;ckner, A., Ladinsky, M.S., Quinodoz, S.A., Badroos, J.M., Viliunas, J.W., Kishi, Y., Wagner, J.M., Miller, D.R., et al. (2024). The genomic and cellular basis of biosynthetic innovation in rove beetles. Cell \u003cem\u003e0\u003c/em\u003e. https://doi.org/10.1016/j.cell.2024.05.012.\u003c/li\u003e\n\u003cli\u003eTurcu, A.F., Rege, J., Auchus, R.J., and Rainey, W.E. (2020). 11-Oxygenated androgens in health and disease. Nat. Rev. Endocrinol. \u003cem\u003e16\u003c/em\u003e, 284\u0026ndash;296. https://doi.org/10.1038/s41574-020-0336-x.\u003c/li\u003e\n\u003cli\u003eTan, D.-X., Hardeland, R., Back, K., Manchester, L.C., Alatorre-Jimenez, M.A., and Reiter, R.J. (2016). On the significance of an alternate pathway of melatonin synthesis via 5-methoxytryptamine: comparisons across species. J. Pineal Res. \u003cem\u003e61\u003c/em\u003e, 27\u0026ndash;40. https://doi.org/10.1111/jpi.12336.\u003c/li\u003e\n\u003cli\u003eBurkhardt, I., de Rond, T., Chen, P.Y.-T., and Moore, B.S. (2022). Ancient plant-like terpene biosynthesis in corals. Nat. Chem. Biol. \u003cem\u003e18\u003c/em\u003e, 664\u0026ndash;669. https://doi.org/10.1038/s41589-022-01026-2.\u003c/li\u003e\n\u003cli\u003eScesa, P.D., Lin, Z., and Schmidt, E.W. (2022). Ancient defensive terpene biosynthetic gene clusters in the soft corals. Nat. Chem. Biol. \u003cem\u003e18\u003c/em\u003e, 659\u0026ndash;663. https://doi.org/10.1038/s41589-022-01027-1.\u003c/li\u003e\n\u003cli\u003eMoon, N.G., and Harned, A.M. (2018). Synthetic explorations of the briarane jungle: progress in developing a synthetic route to a common family of diterpenoid natural products. R. Soc. Open Sci. \u003cem\u003e5\u003c/em\u003e, 172280. https://doi.org/10.1098/rsos.172280.\u003c/li\u003e\n\u003cli\u003eLa, M.-P., Li, J., Li, C., Tang, H., Liu, B.-S., Sun, P., Zhuang, C.-L., Li, T.-J., and Zhang, W. (2014). Briarane Diterpenoids from the Gorgonian Dichotella gemmacea. Mar. Drugs \u003cem\u003e12\u003c/em\u003e, 6178\u0026ndash;6189. https://doi.org/10.3390/md12126178.\u003c/li\u003e\n\u003cli\u003eBarsby, T., and Kubanek, J. (2005). Isolation and Structure Elucidation of Feeding Deterrent Diterpenoids from the Sea Pansy, Renilla reniformis. J. Nat. Prod. \u003cem\u003e68\u003c/em\u003e, 511\u0026ndash;516. https://doi.org/10.1021/np049609u.\u003c/li\u003e\n\u003cli\u003eFenical, W., and Pawlik, J.R. (1991). Defensive properties of secondary metabolites from the Caribbean gorgonian coral Erythropodium caribaeorum. Mar. Ecol. Prog. Ser. \u003cem\u003e75\u003c/em\u003e, 1\u0026ndash;8.\u003c/li\u003e\n\u003cli\u003eHarvell, C.D., Fenical, W., Roussis, V., Ruesink, J.L., Griggs, C.C., and Greene, C.H. (1993). Local and geographic variation in the defensive chemistry of a West Indian gorgonian coral (Briareum asbestinum). Mar. Ecol. Prog. Ser. \u003cem\u003e93\u003c/em\u003e, 165\u0026ndash;173.\u003c/li\u003e\n\u003cli\u003eKsebati, M.B., and Schmitz, F.J. (1986). Diterpenes from a Soft Coral, Minabea SP., from Truk Lagoon. Bull. Soci\u0026eacute;t\u0026eacute;s Chim. Belg. \u003cem\u003e95\u003c/em\u003e, 835\u0026ndash;851. https://doi.org/10.1002/bscb.19860950910.\u003c/li\u003e\n\u003cli\u003eWratten, S.J., and Faulkner, D.J. (1979). Some diterpenes from the sea pen stylatula SP. Tetrahedron \u003cem\u003e35\u003c/em\u003e, 1907\u0026ndash;1912. https://doi.org/10.1016/0040-4020(79)85047-4.\u003c/li\u003e\n\u003cli\u003eYang, J., Zhang, S., Qia, S., Pan, J., Qiu, Y., Tao, S., Yin, H., and Li, Q. (2007). Briarane-type diterpenoids from the China gorgonian coral \u003cem\u003eSubergorgia reticulata\u003c/em\u003e. Biochem. Syst. Ecol. \u003cem\u003e35\u003c/em\u003e, 770\u0026ndash;773. https://doi.org/10.1016/j.bse.2007.06.003.\u003c/li\u003e\n\u003cli\u003eMarinLit - A database of the marine natural products literature https://marinlit.rsc.org/.\u003c/li\u003e\n\u003cli\u003eKeifer, P.A., Rinehart, K.L., and Hooper, I.R. (1986). Renillafoulins, antifouling diterpenes from the sea pansy Renilla reniformis (Octocorallia). J. Org. Chem. \u003cem\u003e51\u003c/em\u003e, 4450\u0026ndash;4454. https://doi.org/10.1021/jo00373a020.\u003c/li\u003e\n\u003cli\u003eHendrickson, R.L., and Cardellina, J.H. (1986). Structure and stereochemistry of insecticidal diterpenes from the sea penXXX. Tetrahedron \u003cem\u003e42\u003c/em\u003e, 6565\u0026ndash;6570. https://doi.org/10.1016/S0040-4020(01)88120-5.\u003c/li\u003e\n\u003cli\u003eNurco, D.J., Conklin, D.E., Shapiro, N.S., and Tran, E. (2010). The absolute structure of ptilosarcenone 2.5-hydrate, a diterpenoid briarane from the orange sea pen Ptilosarcus gurneyi (Gray). Acta Crystallogr. Sect. E Struct. Rep. Online \u003cem\u003e67\u003c/em\u003e, o181\u0026ndash;o182. https://doi.org/10.1107/S1600536810050142.\u003c/li\u003e\n\u003cli\u003eLi, C., La, M.-P., Sun, P., Kurtan, T., Mandi, A., Tang, H., Liu, B.-S., Yi, Y.-H., Li, L., and Zhang, W. (2011). Bioactive (3Z,5E)-11,20-Epoxybriara-3,5-dien-7,18-olide Diterpenoids from the South China Sea Gorgonian Dichotella gemmacea. Mar. Drugs \u003cem\u003e9\u003c/em\u003e, 1403\u0026ndash;1418. https://doi.org/10.3390/md9081403.\u003c/li\u003e\n\u003cli\u003ePordesimo, E.O., Schmitz, F.J., Ciereszko, L.S., Hossain, M.B., and Van der Helm, D. (1991). New briarein diterpenes from the Caribbean gorgonians Erythropodium caribaeorum and Briareum sp. J. Org. Chem. \u003cem\u003e56\u003c/em\u003e, 2344\u0026ndash;2357. https://doi.org/10.1021/jo00007a020.\u003c/li\u003e\n\u003cli\u003eMolina, S.L., Forero, A.M., Ayala, F.I., Puyana, M., Zea, S., Castellanos, L., Mu\u0026ntilde;oz, D., Arboleda, G., Sandoval-Hern\u0026aacute;ndez, A.G., and Ramos, F.A. (2020). Metabolic Profiling of the Soft Coral Erythropodium caribaeorum (Alcyonacea: Anthothelidae) from the Colombian Caribbean Reveals Different Chemotypes. Mar. Drugs \u003cem\u003e18\u003c/em\u003e, 4. https://doi.org/10.3390/md18010004.\u003c/li\u003e\n\u003cli\u003eJiang, J.B., Quattrini, A.M., Francis, W.R., Ryan, J.F., Rodr\u0026iacute;guez, E., and McFadden, C.S. (2019). A hybrid de novo assembly of the sea pansy (Renilla muelleri) genome. GigaScience \u003cem\u003e8\u003c/em\u003e, giz026. https://doi.org/10.1093/gigascience/giz026.\u003c/li\u003e\n\u003cli\u003eLedoux, J.-B., Gomez-Garrido, J., Cruz, F., Ferreira, F.C., Matos, A., Sarropoulou, X., Ramirez-Calero, S., Aurelle, D., Lopez-Sendino, P., Grayson, N., et al. (2024). Chromosome-level genome assembly and annotation of Corallium rubrum: a Mediterranean coral threatened by overharvesting and climate change. Preprint at bioRxiv, https://doi.org/10.1101/2024.07.13.603384 https://doi.org/10.1101/2024.07.13.603384.\u003c/li\u003e\n\u003cli\u003eCorominas, M., Marqu\u0026egrave;s-Bonet, T., Arnedo, M.A., Bay\u0026eacute;s, M., Belmonte, J., Escriv\u0026agrave;, H., Fern\u0026aacute;ndez, R., Gabald\u0026oacute;n, T., Garnatje, T., Germain, J., et al. (2024). The Catalan initiative for the Earth BioGenome Project: contributing local data to global biodiversity genomics. NAR Genomics Bioinforma. \u003cem\u003e6\u003c/em\u003e, lqae075. https://doi.org/10.1093/nargab/lqae075.\u003c/li\u003e\n\u003cli\u003eHu, M., Zheng, X., Fan, C.-M., and Zheng, Y. (2020). Lineage dynamics of the endosymbiotic cell type in the soft coral Xenia. Nature \u003cem\u003e582\u003c/em\u003e, 534\u0026ndash;538. https://doi.org/10.1038/s41586-020-2385-7.\u003c/li\u003e\n\u003cli\u003eIp, J.C.-H., Ho, M.-H., Chan, B.K.K., and Qiu, J.-W. (2023). A draft genome assembly of reef-building octocoral Heliopora coerulea. Sci. Data \u003cem\u003e10\u003c/em\u003e, 381. https://doi.org/10.1038/s41597-023-02291-z.\u003c/li\u003e\n\u003cli\u003ePankov, K.V., McArthur, A.G., Gold, D.A., Nelson, D.R., Goldstone, J.V., and Wilson, J.Y. (2021). The cytochrome P450 (CYP) superfamily in cnidarians. Sci. Rep. \u003cem\u003e11\u003c/em\u003e, 9834. https://doi.org/10.1038/s41598-021-88700-y.\u003c/li\u003e\n\u003cli\u003eDeJong, J.M., Liu, Y., Bollon, A.P., Long, R.M., Jennewein, S., Williams, D., and Croteau, R.B. (2006). Genetic engineering of taxol biosynthetic genes in Saccharomyces cerevisiae. Biotechnol. Bioeng. \u003cem\u003e93\u003c/em\u003e, 212\u0026ndash;224. https://doi.org/10.1002/bit.20694.\u003c/li\u003e\n\u003cli\u003ePompon, D., Louerat, B., Bronine, A., and Urban, P. (1996). [6] Yeast expression of animal and plant P450s in optimized redox environments. In Methods in Enzymology, E. F. Johnson and M. R. Waterman, eds. (Academic Press), pp. 51\u0026ndash;64. https://doi.org/10.1016/S0076-6879(96)72008-6.\u003c/li\u003e\n\u003cli\u003eRo, D.-K., Paradise, E.M., Ouellet, M., Fisher, K.J., Newman, K.L., Ndungu, J.M., Ho, K.A., Eachus, R.A., Ham, T.S., Kirby, J., et al. (2006). Production of the antimalarial drug precursor artemisinic acid in engineered yeast. Nature \u003cem\u003e440\u003c/em\u003e, 940\u0026ndash;943. https://doi.org/10.1038/nature04640.\u003c/li\u003e\n\u003cli\u003eRudolf, J.D., and Chang, C.-Y. (2020). Terpene synthases in disguise: enzymology, structure, and opportunities of non-canonical terpene synthases. Nat. Prod. Rep. \u003cem\u003e37\u003c/em\u003e, 425\u0026ndash;463. https://doi.org/10.1039/C9NP00051H.\u003c/li\u003e\n\u003cli\u003eMcFadden, C.S., Ofwegen, L.P. van, and Quattrini, A.M. (2022). Revisionary systematics of Octocorallia (Cnidaria: Anthozoa) guided by phylogenomics. Bull. Soc. Syst. Biol. \u003cem\u003e1\u003c/em\u003e. https://doi.org/10.18061/bssb.v1i3.8735.\u003c/li\u003e\n\u003cli\u003eBarsby, T., and Kubanek, J. (2005). Isolation and structure elucidation of feeding deterrent diterpenoids from the sea pansy, Renilla reniformis. J Nat Prod \u003cem\u003e68\u003c/em\u003e, 511\u0026ndash;516. https://doi.org/10.1021/np049609u.\u003c/li\u003e\n\u003cli\u003eHe, J., Xin, P., Ma, X., Chu, J., and Wang, G. (2020). Gibberellin Metabolism in Flowering Plants: An Update and Perspectives. Front. Plant Sci. \u003cem\u003e11\u003c/em\u003e. https://doi.org/10.3389/fpls.2020.00532.\u003c/li\u003e\n\u003cli\u003eCleves, P.A., Strader, M.E., Bay, L.K., Pringle, J.R., and Matz, M.V. (2018). CRISPR/Cas9-mediated genome editing in a reef-building coral. Proc. Natl. Acad. Sci. \u003cem\u003e115\u003c/em\u003e, 5235\u0026ndash;5240. https://doi.org/10.1073/pnas.1722151115.\u003c/li\u003e\n\u003cli\u003eLyu, C., Chen, T., Qiang, B., Liu, N., Wang, H., Zhang, L., and Liu, Z. (2021). CMNPD: a comprehensive marine natural products database towards facilitating drug discovery from the ocean. Nucleic Acids Res. \u003cem\u003e49\u003c/em\u003e, D509\u0026ndash;D515. https://doi.org/10.1093/nar/gkaa763.\u003c/li\u003e\n\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":true,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"nature-portfolio","isNatureJournal":true,"hasQc":false,"allowDirectSubmit":false,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"","title":"Nature Portfolio","twitterHandle":"","acdcEnabled":false,"dfaEnabled":false,"editorialSystem":"ejp","reportingPortfolio":"","inReviewEnabled":true,"inReviewRevisionsEnabled":false},"keywords":"","lastPublishedDoi":"10.21203/rs.3.rs-4859447/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-4859447/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eOctocorals are unique among metazoans in their prolific production of bioactive terpenoid natural products that rival the chemical diversity of plants and microbes. We recently established that these cnidarians uniformly express terpene cyclases and that their encoding genes often reside within putative biosynthetic gene clusters (BGCs), a feature uncommon in animal genomes. In this work, we report the discovery and characterization of a widespread gene cluster family for the biosynthesis of briarane diterpenoids that number over 700 molecules specific to the Scleralcyonaceans, one of the two octocoral orders. We sequenced five genomes from evolutionarily distinct families of briarane-producing octocorals to complement three publicly available briarane-producing coral genomes, enabling the discovery of a conserved five-gene cluster composed of a terpene cyclase, three cytochrome P450s, and a short-chain dehydrogenase. Using \u003cem\u003eEscherichia coli\u003c/em\u003e and \u003cem\u003eSaccharomyces cerevisiae\u003c/em\u003e as hosts and homologous briarane biosynthesis genes from eight corals, we reconstituted the biosynthesis of cembrene B γ-lactone, which contains the γ-lactone structural feature distinctive of briarane diterpenoids. The discovery of the genomic basis of briarane biosynthesis establishes that animals, like microbes and plants, employ gene cluster families to produce specialized metabolites. Further, the presence of BGCs in octocoral proves that the formation and maintenance of BGCs over evolutionary time is a more widespread phenomenon in specialized metabolite biosynthesis than previously realized.\u003c/p\u003e","manuscriptTitle":"A widespread metabolic gene cluster family in metazoans","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2024-08-21 06:57:31","doi":"10.21203/rs.3.rs-4859447/v1","editorialEvents":[],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"nature-chemical-biology","isNatureJournal":true,"hasQc":false,"allowDirectSubmit":false,"externalIdentity":"nchembio","sideBox":"Learn more about [Nature Chemical Biology](http://www.nature.com/nchembio/)","snPcode":"","submissionUrl":"","title":"Nature Chemical Biology","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"ejp","reportingPortfolio":"Nature Research","inReviewEnabled":true,"inReviewRevisionsEnabled":false}}],"origin":"","ownerIdentity":"0414e94b-5580-4de4-b1f6-905e9cddd044","owner":[],"postedDate":"August 21st, 2024","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[{"id":35778244,"name":"Biological sciences/Biochemistry/Enzymes/Multienzyme complexes"},{"id":35778245,"name":"Earth and environmental sciences/Ocean sciences/Marine biology"},{"id":35778246,"name":"Biological sciences/Biotechnology/Metabolic engineering"},{"id":35778247,"name":"Biological sciences/Biological techniques/Sequencing/DNA sequencing"}],"tags":[],"updatedAt":"2025-06-14T07:07:02+00:00","versionOfRecord":{"articleIdentity":"rs-4859447","link":"https://doi.org/10.1038/s41589-025-01927-y","journal":{"identity":"nature-chemical-biology","isVorOnly":false,"title":"Nature Chemical Biology"},"publishedOn":"2025-06-13 04:00:00","publishedOnDateReadable":"June 13th, 2025"},"versionCreatedAt":"2024-08-21 06:57:31","video":"","vorDoi":"10.1038/s41589-025-01927-y","vorDoiUrl":"https://doi.org/10.1038/s41589-025-01927-y","workflowStages":[]},"version":"v1","identity":"rs-4859447","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-4859447","identity":"rs-4859447","version":["v1"]},"buildId":"qtupq5eGEP_6zYnWcrvyt","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}