Extreme thrombocytosis with an aggressive evolution harboring a novel variant of calreticulin (CALR) in exon 3.

We describe the case of a patient with extreme thrombocytosis whose evolution was rapidly fatal. No cause of secondary thrombocytosis was found. There was no sign of myelofibrosis but the megakaryocytes were small and dysplastic. The patient presented a calreticulin (CALR) variant in exon 3 (C105S), as well as concomitant mutations of ASXL1, U2AF1, and EZH2. This variant of CALR has never been described before, and after sorting, all identified mutations were found in myeloid cells but not in lymphoid cells. Therefore, the diagnosis of a frontier case of myelodysplastic syndrome/myeloproliferative neoplasm (MDS/MPN) was made. A treatment with hydroxycarbamide was started because of a high risk of thrombosis. Upon worsening of the hematological status two new mutations appeared, SETBP1 and ETV6, and the CALR mutation was still detectable, as well as the three other mutations found in the chronic stage. Our results show that this variant could contribute to MDS/MPN pathogenesis in that patient.

SYNPHONI: scale-free and phylogeny-aware reconstruction of synteny conservation and transformation across animal genomes.

SUMMARY: Current approaches detect conserved genomic order either at chromosomal (macrosynteny) or at subchromosomal scales (microsynteny). The latter generally requires collinearity and hard thresholds on syntenic region size, thus excluding a major proportion of syntenies with recent expansions or minor rearrangements. 'SYNPHONI' bridges the gap between micro- and macrosynteny detection, providing detailed information on both synteny conservation and transformation throughout the evolutionary history of animal genomes. AVAILABILITY AND IMPLEMENTATION: Source code is freely available at https://github.com/nsmro/SYNPHONI, implemented in Python 3.9. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Emergence of distinct syntenic density regimes is associated with early metazoan genomic transitions.

BACKGROUND: Animal genomes are strikingly conserved in terms of local gene order (microsynteny). While some of these microsyntenies have been shown to be coregulated or to form gene regulatory blocks, the diversity of their genomic and regulatory properties across the metazoan tree of life remains largely unknown. RESULTS: Our comparative analyses of 49 animal genomes reveal that the largest gains of synteny occurred in the last common ancestor of bilaterians and cnidarians and in that of bilaterians. Depending on their node of emergence, we further show that novel syntenic blocks are characterized by distinct functional compositions (Gene Ontology terms enrichment) and gene density properties, such as high, average and low gene density regimes. This is particularly pronounced among bilaterian novel microsyntenies, most of which fall into high gene density regime associated with higher gene coexpression levels. Conversely, a majority of vertebrate novel microsyntenies display a low gene density regime associated with lower gene coexpression levels. CONCLUSIONS: Our study provides first evidence for evolutionary transitions between different modes of microsyntenic block regulation that coincide with key events of metazoan evolution. Moreover, the microsyntenic profiling strategy and interactive online application (Syntenic Density Browser, available at: http://synteny.csb.univie.ac.at/ ) we present here can be used to explore regulatory properties of microsyntenic blocks and predict their coexpression in a wide-range of animal genomes.

Kinome expression profiling to target new therapeutic avenues in multiple myeloma.

Multiple myeloma (MM) account for approximately 10% of hematological malignancies and is the second most common hematological disorder. Kinases inhibitors are widely used and their efficiency for the treatment of cancers has been demonstrated. Here, in order to identify kinases of potential therapeutic interest for the treatment of MM, we investigated the prognostic impact of the kinome expression profile in large cohorts of patients. We identified 36 kinome-related genes significantly linked with a prognostic value to MM, and built a kinome index based on their expression. The Kinome Index (KI) is linked to prognosis, proliferation, differentiation, and relapse in MM. We then tested inhibitors targeting seven of the identified protein kinas-es (PBK, SRPK1, CDC7-DBF4, MELK, CHK1, PLK4, MPS1/TTK) in human myeloma cell lines. All tested inhibitors significantly reduced the viability of myeloma cell lines, and we confirmed the potential clinical interest of three of them on primary myeloma cells from patients. In addition, we demonstrated their ability to potentialize the toxicity of conventional treatments, including Melphalan and Lenalidomide. This highlights their potential beneficial effect in myeloma therapy. Three kinases inhibitors (CHK1i, MELKi and PBKi) overcome resistance to Lenalidomide, while CHK1, PBK and DBF4 inhibitors re-sensitize Melphalan resistant cell line to this conventional therapeutic agent. Altogether, we demonstrate that kinase inhibitors could be of therapeutic interest especially in high-risk myeloma patients defined by the KI. CHEK1, MELK, PLK4, SRPK1, CDC7-DBF4, MPS1/TTK and PBK inhibitors could represent new treatment options either alone or in combination with Melphalan or IMiD for refractory/relapsing myeloma patients.

The hydroxymethylome of multiple myeloma identifies FAM72D as a 1q21 marker linked to proliferation.

Cell identity relies on the cross-talk between genetics and epigenetics and their impact on gene expression. Oxidation of 5-methylcytosine (5mC) into 5-hydroxymethylcytosine (5hmC) is the first step of an active DNA demethylation process occurring mainly at enhancers and gene bodies and, as such, participates in processes governing cell identity in normal and pathological conditions. Although genetic alterations are well documented in multiple myeloma (MM), epigenetic alterations associated with this disease have not yet been thoroughly analyzed. To gain insight into the biology of MM, genome-wide 5hmC profiles were obtained and showed that regions enriched in this modified base overlap with MM enhancers and super enhancers and are close to highly expressed genes. Through the definition of a MM-specific 5hmC signature, we identified FAM72D as a poor prognostic gene located on 1q21, a region amplified in high risk myeloma. We further uncovered that FAM72D functions as part of the FOXM1 transcription factor network controlling cell proliferation and survival and we evidenced an increased sensitivity of cells expressing high levels of FOXM1 and FAM72 to epigenetic drugs targeting histone deacetylases and DNA methyltransferases.

A wnt2 ortholog in the sea urchin Paracentrotus lividus.

Members of the wnt gene family encode secreted glycoproteins that mediate critical intercellular communications in metazoans. Large-scale genome and transcriptome analyses have shown that this family is composed of 13 distinct subfamilies. These analyses have further established that the number of wnt genes per subfamily varies significantly between metazoan phyla, highlighting that gene duplication and gene loss events have shaped the complements of wnt genes during evolution. In sea urchins, for example, previous work reported the absence of representatives of both the WNT2 and WNT11 subfamilies in two different species, Paracentrotus lividus and Strongylocentrotus purpuratus. Recently, however, we identified a gene encoding a WNT2 ortholog in P. lividus and, based on that finding, we also reanalyzed the genome of S. purpuratus. Yet, we found no evidence of a bona fide wnt2 gene in S. purpuratus. Furthermore, we established that the P. lividus wnt2 gene is selectively expressed in vegetal tissues during embryogenesis, in a pattern that is similar, although not identical, to that of other P. lividus wnt genes. Taken together, this study amends previous work on the P. lividus wnt complement and reveals an unexpected variation in the number of wnt genes between closely related sea urchin species.

Retinoic acid signaling and neurogenic niche regulation in the developing peripheral nervous system of the cephalochordate amphioxus.

The retinoic acid (RA) signaling pathway regulates axial patterning and neurogenesis in the developing central nervous system (CNS) of chordates, but little is known about its roles during peripheral nervous system (PNS) formation and about how these roles might have evolved. This study assesses the requirement of RA signaling for establishing a functional PNS in the cephalochordate amphioxus, the best available stand-in for the ancestral chordate condition. Pharmacological manipulation of RA signaling levels during embryogenesis reduces the ability of amphioxus larvae to respond to sensory stimulation and alters the number and distribution of ectodermal sensory neurons (ESNs) in a stage- and context-dependent manner. Using gene expression assays combined with immunohistochemistry, we show that this is because RA signaling specifically acts on a small population of soxb1c-expressing ESN progenitors, which form a neurogenic niche in the trunk ectoderm, to modulate ESN production during elongation of the larval body. Our findings reveal an important role for RA signaling in regulating neurogenic niche activity in the larval amphioxus PNS. Although only few studies have addressed this issue so far, comparable RA signaling functions have been reported for neurogenic niches in the CNS and in certain neurogenic placode derivatives of vertebrates. Accordingly, the here-described mechanism is likely a conserved feature of chordate embryonic and adult neural development.

Global miRNA expression analysis identifies novel key regulators of plasma cell differentiation and malignant plasma cell.

MicroRNAs (miRNAs) are small noncoding RNAs that attenuate expression of their mRNA targets. Here, we developed a new method and an R package, to easily infer candidate miRNA-mRNA target interactions that could be functional during a given biological process. Using this method, we described, for the first time, a comprehensive integrated analysis of miRNAs and mRNAs during human normal plasma cell differentiation (PCD). Our results reveal 63 miRNAs with significant temporal changes in their expression during normal PCD. We derived a high-confidence network of 295 target relationships comprising 47 miRNAs and 141 targets. These relationships include new examples of miRNAs that appear to coordinately regulate multiple members of critical pathways associated with PCD. Consistent with this, we have experimentally validated a role for the miRNA-30b/c/d-mediated regulation of key PCD factors (IRF4, PRDM1, ELL2 and ARID3A). Furthermore, we found that 24 PCD stage-specific miRNAs are aberrantly overexpressed in multiple myeloma (MM) tumor plasma cells compared to their normal counterpart, suggesting that MM cells frequently acquired expression changes in miRNAs already undergoing dynamic expression modulation during normal PCD. Altogether, our analysis identifies candidate novel key miRNAs regulating networks of significance for normal PCD and malignant plasma cell biology.

DNMTi/HDACi combined epigenetic targeted treatment induces reprogramming of myeloma cells in the direction of normal plasma cells.

BACKGROUND: Multiple myeloma (MM) is the second most common hematologic malignancy. Aberrant epigenetic modifications have been reported in MM and could be promising therapeutic targets. As response rates are overall limited but deep responses occur, it is important to identify those patients who could indeed benefit from epigenetic-targeted therapy. METHODS: Since HDACi and DNMTi combination have potential therapeutic value in MM, we aimed to build a GEP-based score that could be useful to design future epigenetic-targeted combination trials. In addition, we investigated the changes in GEP upon HDACi/DNMTi treatment. RESULTS: We report a new gene expression-based score to predict MM cell sensitivity to the combination of DNMTi/HDACi. A high Combo score in MM patients identified a group with a worse overall survival but a higher sensitivity of their MM cells to DNMTi/HDACi therapy compared to a low Combo score. In addition, treatment with DNMTi/HDACi downregulated IRF4 and MYC expression and appeared to induce a mature BMPC plasma cell gene expression profile in myeloma cell lines. CONCLUSION: In conclusion, we developed a score for the prediction of primary MM cell sensitivity to DNMTi/HDACi and found that this combination could be beneficial in high-risk patients by targeting proliferation and inducing maturation.

A comprehensive survey of wnt and frizzled expression in the sea urchin Paracentrotus lividus.

WNT signaling is, in all multicellular animals, an essential intercellular communication pathway that is critical for shaping the embryo. At the molecular level, WNT signals can be transmitted by several transduction cascades, all activated chiefly by the binding of WNT ligands to receptors of the FRIZZLED family. The first step in assessing the biological functions of WNT signaling during embryogenesis is thus the establishment of the spatiotemporal expression profiles of wnt and frizzled genes in the course of embryonic development. To this end, using quantitative polymerase chain reaction, Northern blot, and in situ hybridization assays, we report here the comprehensive expression patterns of all 11 wnt and 4 frizzled genes present in the genome of the sea urchin Paracentrotus lividus during its embryogenesis. Our findings indicate that the expression of these wnt ligands and frizzled receptors is highly dynamic in both time and space. We further establish that all wnt genes are chiefly transcribed in the vegetal hemisphere of the embryo, whereas expression of the frizzled genes is distributed more widely across the embryonic territories. Thus, in P. lividus, WNT ligands might act both as short- and long-range signaling molecules that may operate in all cell lineages and tissues to control various developmental processes during embryogenesis.

Influence of atmospheric pressure on infrarenal abdominal aortic aneurysm rupture.

BACKGROUND: Meteorologic conditions have a significant impact on the occurrence of cardiovascular events. Previous studies have shown that abdominal aortic aneurysm rupture (AAAR) may be associated with atmospheric pressure, with conflicting results. Therefore, we aimed to further investigate the nature of the correlation between atmospheric pressure variations and AAAR. METHOD: Hospital admissions related to AAAR between 2005-2009 were assessed in 19 districts of metropolitan France and correlated with geographically and date-matched mean atmospheric pressures. In parallel and from 2005-2009, all fatal AAARs as reported by death certificates were assessed nationwide and correlated to local atmospheric pressures at the time of aortic rupture. RESULTS: Four hundred ninety-four hospital admissions related to AAAR and 6,358 deaths nationwide by AAAR were identified between 2005-2009. Both in-hospital ruptures and aneurysm-related mortality had seasonal variations, with peak/trough incidences in January and June, respectively. Atmospheric pressure peaks occurred during winter. Univariate analysis revealed a significant association (P < 0.001) of high mean atmospheric pressure values and AAAR. After multivariate analysis, mean maximum 1-month prerupture atmospheric pressure had a persistent correlation with both in-hospital relative risk (1.05 [95% confidence interval: 1.03-1.06]; P < 0.0001) and aneurysm rupture-related mortality relative risk (1.02 [95% confidence interval: 1.01-1.03]; P < 0.0001). CONCLUSION: The annual incidence of AAAR is nonhomogeneous with a peak incidence in winter, and is independently associated with mean maximum 1-month prerupture atmospheric pressure.

Efficacy and safety of ultrasound-guided distal blocks for analgesia without motor blockade after ambulatory hand surgery.

PURPOSE: To assess the suitability of ultrasound-guided (USG), single-injection distal block(s) for pain management after outpatient hand and wrist bone surgery. METHODS: We conducted a retrospective review of 125 of 198 consecutive ambulatory surgery patients who underwent hand and wrist bone surgery between June 2010 and January 2012. All patients received a USG axillary block using a short-acting local anesthetic (lidocaine) and secondary 1, 2, or 3 (median, radial, or ulnar) USG distal analgesic block(s) using a long-acting local anesthetic (ropivacaine). All patients were contacted by phone on the first postoperative day. All patients received a concomitant prescription of acetaminophen and nonsteroidal anti-inflammatory drugs with opioids as a rescue treatment. Effectiveness and duration of the distal nerve blocks, compliance with analgesic treatment and rescue opioids requirement, opioid-related side effects, prolonged upper limb motor block, quality of sleep on first postoperative night, and patient satisfaction were evaluated. RESULTS: Most distal analgesic blocks were effective (120 of 125; 96%), with an average duration of nearly 12 hours On the first day after surgery, 28 patients (23%) had a numeric verbal scale greater than 3, although 14 of them had taken the rescue opioids. No patient reported prolonged motor blockade or insensate limb. Opioid-related side effects occurred in 23% of patients. CONCLUSIONS: After hand and wrist bone surgery, USG selective distal blocks using a long-acting local anesthetic, combined with oral analgesics, were effective in a large majority of patients. However, pain control was suboptimal for some especially painful procedures such as wrist surgery, trapeziometacarpal arthrodesis, and finger amputation. TYPE OF STUDY/LEVEL OF EVIDENCE: Therapeutic IV.

Characterization of a transitional preplasmablast population in the process of human B cell to plasma cell differentiation.

The early steps of differentiation of human B cells into plasma cells are poorly known. We report a transitional population of CD20(low/-)CD38(-) preplasmablasts along differentiation of human memory B cells into plasma cells in vitro. Preplasmablasts lack documented B cell or plasma cell (CD20, CD38, and CD138) markers, express CD30 and IL-6R, and secrete Igs at a weaker level than do plasmablasts or plasma cells. These preplasmablasts further differentiate into CD20(-)CD38(high)CD138(-) plasmablasts and then CD20(-)CD38(high)CD138(+) plasma cells. Preplasmablasts were fully characterized in terms of whole genome transcriptome profiling and phenotype. Preplasmablasts coexpress B and plasma cell transcription factors, but at a reduced level compared with B cells, plasmablasts, or plasma cells. They express the unspliced form of XBP1 mRNA mainly, whereas plasmablasts and plasma cells express essentially the spliced form. An in vivo counterpart (CD19(+)CD20(low/-)CD38(-)IL-6R(+) cells) of in vitro-generated preplasmablasts could be detected in human lymph nodes (0.06% of CD19(+) cells) and tonsils (0.05% of CD19(+) cells). An open access "B to Plasma Cell Atlas," which makes it possible to interrogate gene expression in the process of B cell to plasma cell differentiation, is provided. Taken together, our findings show the existence of a transitional preplasmablast population using an in vitro model of plasma cell generation and of its in vivo counterpart in various lymphoid tissues.

Robust 3D modeling reveals spatiosyntenic properties of animal genomes.

Animal genomes are organized into chromosomes that are remarkably conserved in their gene content, forming distinct evolutionary units (synteny). Using versatile chromosomal modeling, we infer three-dimensional topology of genomes from representative clades spanning the earliest animal diversification. We apply a partitioning approach using interaction spheres to compensate for varying quality of topological data. Using comparative genomics approaches, we test whether syntenic signal at gene pair, local, and whole chromosomal scale is reflected in the reconstructed spatial organization. We identify evolutionarily conserved three-dimensional networks at all syntenic scales revealing novel evolutionarily maintained interactors associated with known conserved local gene linkages (such as hox). We thus present evidence for evolutionary constraints that are associated with three-, rather than just two-, dimensional animal genome organization, which we term spatiosynteny. As more accurate topological data become available, together with validation approaches, spatiosynteny may become relevant in understanding the functionality behind the observed conservation of animal chromosomes.

Combined inhibition of Wee1 and Chk1 as a therapeutic strategy in multiple myeloma.

Multiple myeloma (MM) is a hematological malignancy characterized by an abnormal clonal proliferation of malignant plasma cells. Despite the introduction of novel agents that have significantly improved clinical outcome, most patients relapse and develop drug resistance. MM is characterized by genomic instability and a high level of replicative stress. In response to replicative and DNA damage stress, MM cells activate various DNA damage signaling pathways. In this study, we reported that high CHK1 and WEE1 expression is associated with poor outcome in independent cohorts of MM patients treated with high dose melphalan chemotherapy or anti-CD38 immunotherapy. Combined targeting of Chk1 and Wee1 demonstrates synergistic toxicities on MM cells and was associated with higher DNA double-strand break induction, as evidenced by an increased percentage of γH2AX positive cells subsequently leading to apoptosis. The therapeutic interest of Chk1/Wee1 inhibitors' combination was validated on primary MM cells of patients. The toxicity was specific of MM cells since normal bone marrow cells were not significantly affected. Using deconvolution approach, MM patients with high CHK1 expression exhibited a significant lower percentage of NK cells whereas patients with high WEE1 expression displayed a significant higher percentage of regulatory T cells in the bone marrow. These data emphasize that MM cell adaptation to replicative stress through Wee1 and Chk1 upregulation may decrease the activation of the cell-intrinsic innate immune response. Our study suggests that association of Chk1 and Wee1 inhibitors may represent a promising therapeutic approach in high-risk MM patients characterized by high CHK1 and WEE1 expression.

Inhibition of the Protein Arginine Methyltransferase PRMT5 in High-Risk Multiple Myeloma as a Novel Treatment Approach.

Multiple myeloma (MM) is an incurable clonal plasma cell malignancy. Subsets of patients have high-risk features linked with dismal outcome. Therefore, the need for effective therapeutic options remains high. Here, we used bio-informatic tools to identify novel targets involved in DNA repair and epigenetics and which are associated with high-risk myeloma. The prognostic significance of the target genes was analyzed using publicly available gene expression data of MM patients (TT2/3 and HM cohorts). Hence, protein arginine methyltransferase 5 (PRMT5) was identified as a promising target. Druggability was assessed in OPM2, JJN3, AMO1 and XG7 human myeloma cell lines using the PRMT5-inhibitor EPZ015938. EPZ015938 strongly reduced the total symmetric-dimethyl arginine levels in all cell lines and lead to decreased cellular growth, supported by cell line dependent changes in cell cycle distribution. At later time points, apoptosis occurred, as evidenced by increased AnnexinV-positivity and cleavage of PARP and caspases. Transcriptome analysis revealed a role for PRMT5 in regulating alternative splicing, nonsense-mediated decay, DNA repair and PI3K/mTOR-signaling, irrespective of the cell line type. PRMT5 inhibition reduced the expression of upstream DNA repair kinases ATM and ATR, which may in part explain our observation that EPZ015938 and the DNA-alkylating agent, melphalan, have combinatory effects. Of interest, using a low-dose of mTOR-inhibitor, we observed that cell viability was partially rescued from the effects of EPZ015938, indicating a role for mTOR-related pathways in the anti-myeloma activity of EPZ015938. Moreover, PRMT5 was shown to be involved in splicing regulation of MMSET and SLAMF7, known genes of importance in MM disease. As such, we broaden the understanding of the exact role of PRMT5 in MM disease and further underline its use as a possible therapeutic target.

Comprehensive characterization of the epigenetic landscape in Multiple Myeloma.

Background: Human multiple myeloma (MM) cell lines (HMCLs) have been widely used to understand the molecular processes that drive MM biology. Epigenetic modifications are involved in MM development, progression, and drug resistance. A comprehensive characterization of the epigenetic landscape of MM would advance our understanding of MM pathophysiology and may attempt to identify new therapeutic targets. Methods: We performed chromatin immunoprecipitation sequencing to analyze histone mark changes (H3K4me1, H3K4me3, H3K9me3, H3K27ac, H3K27me3 and H3K36me3) on 16 HMCLs. Results: Differential analysis of histone modification profiles highlighted links between histone modifications and cytogenetic abnormalities or recurrent mutations. Using histone modifications associated to enhancer regions, we identified super-enhancers (SE) associated with genes involved in MM biology. We also identified promoters of genes enriched in H3K9me3 and H3K27me3 repressive marks associated to potential tumor suppressor functions. The prognostic value of genes associated with repressive domains and SE was used to build two distinct scores identifying high-risk MM patients in two independent cohorts (CoMMpass cohort; n = 674 and Montpellier cohort; n = 69). Finally, we explored H3K4me3 marks comparing drug-resistant and -sensitive HMCLs to identify regions involved in drug resistance. From these data, we developed epigenetic biomarkers based on the H3K4me3 modification predicting MM cell response to lenalidomide and histone deacetylase inhibitors (HDACi). Conclusions: The epigenetic landscape of MM cells represents a unique resource for future biological studies. Furthermore, risk-scores based on SE and repressive regions together with epigenetic biomarkers of drug response could represent new tools for precision medicine in MM.

The BLM helicase is a new therapeutic target in multiple myeloma involved in replication stress survival and drug resistance.

Multiple myeloma (MM) is a hematologic cancer characterized by accumulation of malignant plasma cells in the bone marrow. To date, no definitive cure exists for MM and resistance to current treatments is one of the major challenges of this disease. The DNA helicase BLM, whose depletion or mutation causes the cancer-prone Bloom's syndrome (BS), is a central factor of DNA damage repair by homologous recombination (HR) and genomic stability maintenance. Using independent cohorts of MM patients, we identified that high expression of BLM is associated with a poor outcome with a significant enrichment in replication stress signature. We provide evidence that chemical inhibition of BLM by the small molecule ML216 in HMCLs (human myeloma cell lines) leads to cell cycle arrest and increases apoptosis, likely by accumulation of DNA damage. BLM inhibition synergizes with the alkylating agent melphalan to efficiently inhibit growth and promote cell death in HMCLs. Moreover, ML216 treatment re-sensitizes melphalan-resistant cell lines to this conventional therapeutic agent. Altogether, these data suggest that inhibition of BLM in combination with DNA damaging agents could be of therapeutic interest in the treatment of MM, especially in those patients with high BLM expression and/or resistance to melphalan.

Ancestral Role of Ecdysis-Related Neuropeptides in Animal Life Cycle Transitions.

Ecdysis or molting evolved ∼535 mya in Ecdysozoa, the most diverse and species-rich animal superphylum.1 A cascade of ecdysis-related neuropeptides (ERNs) controls the innate behavioral programs required for cuticle shedding in some ecdysozoan lineages (e.g., arthropods)2-12 but is lacking in others (e.g., nematodes).13 We recently reported on the surprisingly ancient bilaterian origin of key ERNs, such as eclosion hormone (EH), crustacean cardioactive neuropeptide (CCAP), myoinhibitory peptide (MIP), bursicon alpha (Bursα), and bursicon beta (Bursß).13,14 Thus, ERNs far predate the emergence of ecdysis, but the question as to their ancestral functions remains unresolved. Here, we compare the ERN toolkits and temporal expression profiles of six ecdysozoans (tardigrades, crustaceans, and insects), eight lophotrochozoans (planarians, annelids, and mollusks), and five deuterostomes (crinoids, sea urchins, and hemichordates). Our results show that the major, coordinated upregulation of ERNs always coincides with a transition between key life history stages, such as hatching in direct developers and metamorphosis in indirect developers. This implies that ERNs already played an ancestral role in the switch from embryonic or larval ontogeny to juvenile maturation in the last common ancestor of Nephrozoa. Consequently, the transcriptional signature of invertebrate life cycle transitions presented here was already in place in the Precambrian and was only secondarily co-opted into regulating the molting process at the dawn of Ecdysozoa.

Ecdysis-related neuropeptide expression and metamorphosis in a non-ecdysozoan bilaterian.

Ecdysis-related neuropeptides (ERNs), including eclosion hormone, crustacean cardioactive peptide, myoinhibitory peptide, bursicon alpha, and bursicon beta regulate molting in insects and crustaceans. Recent evidence further revealed that ERNs likely play an ancestral role in invertebrate life cycle transitions, but their tempo-spatial expression patterns have not been investigated outside Arthropoda. Using RNA-seq and in situ hybridization, we show that ERNs are broadly expressed in the developing nervous system of a mollusk, the polyplacophoran Acanthochitona fascicularis. While some ERN-expressing neurons persist from larval to juvenile stages, others are only present during settlement and metamorphosis. These transient neurons belong to the "ampullary system," a polyplacophoran-specific larval sensory structure. Surprisingly, however, ERN expression is absent from the apical organ, another larval sensory structure that degenerates before settlement is completed in A. fascicularis. Our findings thus support a role of ERNs in A. fascicularis metamorphosis but contradict the common notion that the apical organ-like structures shared by various aquatic invertebrates (i.e., cnidarians, annelids, mollusks, echinoderms) are of general importance for this process.