Characterization of the bHLH family of transcriptional regulators in the acoel S. roscoffensis and their putative role in neurogenesis
© The Author(s) 2018
Received: 25 January 2018
Accepted: 22 March 2018
Published: 29 March 2018
The basic helix-loop-helix (bHLH) family of transcription factors is one of the largest superfamilies of regulatory transcription factors and is widely used in eukaryotic organisms. They play an essential role in a range of metabolic, physiological, and developmental processes, including the development of the nervous system (NS). These transcription factors have been studied in many metazoans, especially in vertebrates but also in early branching metazoan clades such as the cnidarians and sponges. However, currently very little is known about their expression in the most basally branching bilaterian group, the xenacoelomorphs. Recently, our laboratory has characterized the full complement of bHLH in the genome of two members of the Xenacoelomorpha, the xenoturbellid Xenoturbella bocki and the acoel Symsagittifera roscoffensis. Understanding the patterns of bHLH gene expression in members of this phylum (in space and time) provides critical new insights into the conserved roles of the bHLH and their putative specificities in this group. Our focus is on deciphering the specific roles that these genes have in the process of neurogenesis.
Here, we analyze the developmental expression of the whole complement of bHLH genes identified in the acoel S. roscoffensis. Based on their expression patterns, several members of bHLH class A appear to have specific conserved roles in neurogenesis, while other class A genes (as well as members of other classes) have likely taken on more generalized functions. All gene expression patterns are described in embryos and early juveniles.
Our results suggest that the main roles of the bHLH genes of S. roscoffensis are evolutionarily conserved, with a specific subset dedicated to patterning the nervous system: SrAscA, SrAscB, SrHes/Hey, SrNscl, SrSrebp, SrE12/E47 and SrOlig.
Xenacoelomorpha is a phylum constituted by small, mostly marine, benthic worms that share a relatively “simple” morphology (reviewed in [1, 2]). They are bilaterally symmetrical, and their bodies are covered by a ciliated epithelium with a mouth being the only digestive opening to the exterior. Common features found in other metazoans such as the presence of circulatory or excretory systems, anus and coelom are completely absent in this group [3–5]. Xenacoelomorpha is divided into three clades: Xenoturbellida, Nemertodermatida and Acoela . The latter two are grouped into the clade Acoelomorpha, the sister group of Xenoturbellida. For a long time, the phylogenetic affinities of these clades have been a matter of intense debate; see [6–9]. However, the latest phylogenetic analysis carried out by Canon et al.  seemed to resolve this conflict, proposing, based on analysis with strong support, a monophyletic Xenacoelomorpha as the sister group to all remaining bilaterians.
There are several reasons why xenacoelomorphs are an interesting set of biological systems in which to carry out comparative molecular and developmental studies. They belong to a monophyletic group (sharing a common ancestor) which members show a high diversity in the complexity of many regulatory families and in the organization of anatomical architectures (i.e. nervous system, arrangement of musculature, position of the mouth, the morphology of copulatory apparatus). This fact allows us to compare different ways of “constructing” and patterning organ systems within a set of interrelated animals. We should suggest that, in addition, a significant practical advantage of studying the xenacoelomorphs could be the fact that they seem to possess fewer cell types and organs than most bilaterian animals (lacking proper through-gut, nephridia, complex glands), making their system more amenable to our future (comprehensive) research efforts. In addition, a better understanding of the genetic control of developmental processes in xenacoelomorphs will provide key insights into the origins and diversification of the bilaterians [10–12].
One of the key innovations linked to the emergence of Bilateria is the origin of centralized nervous systems. How these compact brains are assembled from simpler nerve nets remains a matter of debate that is mostly grounded in the lack of knowledge that we have on the molecular mechanisms that differentially control the development and assembly of nerve nets, cords and compact brains (see, for instance, [8, 12–15]. The comparative approach should provide an answer. In this context the use of members of the Xenacoelomorpha is ideal, since they have different nervous system morphologies (with variable degrees of condensation) all derived from a single, common, ancestor. Briefly, in the group of xenoturbellids, a unique basiepithelial nerve net surrounds the animal body while some acoels have, in addition to a nerve net, an anteriorly concentrated nervous system [1, 16] (for a comprehensive review, see: ). The nervous system architecture of the most divergent class of acoels (Crucimusculata) represents, most probably, one of the first instances of the acquisition of a compact brain in bilaterian evolution (nephrozoans would have acquired a compact brain independently). Acoel embryos possess a unique early development program that is known as “duet spiral” cleavage . The first three micromere duets give rise to the ectodermal layer including epidermal and neural progenitors. The formation of the organs’ anlage starts in mid-embryonic stages and the symmetrical brain primordium can be observed at the anterior pole, subepidermally, at early stages .
One of the most studied acoel species is Symsagittifera roscoffensis. The nervous system of this acoel is arranged in an anterior domain forming a compact brain with neural cell bodies surrounding a neuropil, divided into two lobes and connected by three commissures [17, 20, 21]. The ventral part of the neuropil projects anteriorly to a commissural ring that surrounds the frontal organ. In addition, three pairs of cords arise from the brain and run along the anterior–posterior body axis, in a specific dorso-ventral distribution. Two dorsal cords arise from the posterior part of the brain (specifically from the third commissure), while the remaining four nerve cords lie more ventrally (a ventral central pair and a latero-ventral pair). The most prominent sensory organ is the acoel statocyst, located in the anterior part of the body and surrounded by the brain neuropil. Anterior to the statocyst are a pair of ocelli consisting of several sensory cells and a pigment cell [12, 16, 20, 21]. Juveniles of S. roscoffensis are about 220 μm after hatching and the brain occupies more than a third of their body length, a striking difference when compared to the adult specimens, where the brain occupies only a small anterior region (approximately an eighth of the animal’s length) [12, 16, 20–22].
In this context, we carried out a systematic characterization of a well-known superfamily of transcriptional regulators, the basic helix-loop-helix (bHLH) proteins. They are widely present in eukaryotes and play an important role in metabolic, physiological and developmental processes [23–26]. Some of them are involved in the development and patterning of the nervous system in many bilaterians (but also in cnidarians (see for instance: [27–29] also reviewed in [26, 30]).
Members of the bHLH superfamily encode for proteins containing a characteristic 60 amino acids long bHLH domain that includes a N-terminal DNA binding basic (b) region followed by two α-helices connected by a loop region (HLH) of variable length. The HLH domain promotes dimerization, allowing the formation of homo- or heterodimeric complexes between different family members. Some bHLHs also include additional domains involved in protein–protein interactions such as the “leucine zipper,” PAS (Per–Arnt–Sim) and the “orange” domains [24, 31].
In the past 10 years, since the pioneering study of the origin and the diversification of the bHLH carried out by , new full sets of bHLH have been identified via the thorough analysis of many metazoan genomes [25, 32–35] (also, in non-metazoans ). The complex group of metazoan bHLH transcription factors have been classified, using molecular phylogenetic analysis, into 48 orthologous families (45 different families sensu Simionato et al.  and three new ones sensu Gyoja et al. ). Based on phylogenetic affinities and general biochemical properties, the orthologous families are organized into six “higher-order” groups named A to F [24, 25, 31, 37, 38]. Group A is especially relevant here since it includes most of the bHLH genes with neurogenic functions in other bilaterians.
These genes were cloned and analyzed in the context of a thorough characterization of the bHLH in the xenacoelomorph genomes. In order to characterize the expression patterns of bHLH genes in xenacoelomorphs, we started by conducting an analysis of this group of genes in one (experimentally amenable) member of the Acoela, the species S. roscoffensis, using both colored in situ hybridization (ISH) in embryos and juveniles and double fluorescence in situ hybridization (FISH) in juveniles. Our data are also discussed in the context of the evolution of neural system patterning.
bHLH sequence identification
Symsagittifera roscoffensis DNA sequences were extracted from the genome assembly and the embryo (mixed stages) transcriptome and were published in . In the mentioned analysis, sequences were extracted using HMMER (v.3.0)  and classified by phylogenetic methods.
Adult specimens of S. roscoffensis were collected in Carantec (Brittany, France) during two of their reproductive periods (in April 2016 and 2017). They were cleaned and placed in petri dishes with fresh seawater until they spawned the cocoons in the media. Every cocoon contains several embryos, mostly synchronized (between ten and twenty). We collected the cocoons at different time intervals, with the time when we detected the fresh spawn being the “zero” time. The other times corresponded to the number of hours elapsed between the detection of fresh spawn and fixing the cocoons for in situ analysis. For instance, a sample from 0 to 6 h contained embryos that developed within the first 6 h. Embryos were separated from the adults by filtering and then treated using 0.01% Pronase (Sigma) and 0.1% thioglycolate (Sigma) in seawater to permeabilize the egg shells. After cleaning them with seawater, embryos were fixed in 4% formaldehyde (methanol free) overnight at 4 °C. A large proportion of cocoons were not fixed, with the aim of obtaining hatchlings later on. These hatchlings were relaxed using 7% magnesium chloride and fixed, as were the embryos. Some specimens were allowed to develop to older juveniles. After fixation all samples were cleaned three times in 1 × PBS and dehydrated progressively with a methanol series (25–50–75% methanol in PBS).
Gene cloning and in situ hybridization
DIG-labeled and fluorescein-labeled RNA probes were synthesized using the DIG-RNA labeling kit from Sigma, following the manufacturer’s instructions (Id_number11277073910, Sigma-Aldrich). After precipitation, the riboprobes were diluted in hybridization buffer to a final working concentration of 2–1 ng/μl (depending of the probe). In situ hybridization on whole embryos and juveniles was performed following the protocol published by  with a few modifications: (1) the main solvents, PB-Tween and PB-Triton were replaced with TNT (0.1 M TRIS–HCl, pH7.5 + 0.15 M NaCl + 0.05% Tween-20 detergent in RNAse-free water); (2) the proteinase K and glycine steps were suppressed, in order to reduce the damage produced to the samples; (3) we included a step in which we increased the temperature to 80 °C during one of the two washes with hybridization buffer (HB), to reduce the background. After this step, we proceed to pre-hybridize overnight. We incubated the sample with the corresponding probe for at least 3 days. The hybridization temperatures were between 55 and 61 °C depending on the probe. The specimens were mounted in 70% glycerol and analyzed using a Zeiss Axiophot microscope (Carl Zeiss MicroImaging GmbH) equipped with a Leica DFC 300FX camera.
Colorimetric in situ protocols highlight the domains of highest gene expression. For detailed aspects of the patterns, we rely on the more sensitive fluorescence in situ alternatives.
For all clones, sense probes were synthesized that were used as negative controls for hybridization.
Double fluorescence in situ hybridization
Juveniles of S. roscoffensis used in FISH analysis needed a photo-bleaching treatment step after the re-hydration. In order to reduce background due to auto-fluorescence, as much as possible, we immersed the specimens in a solution of 1.5% hydrogen peroxide, 5% formamide and 0.5% 0.5 × SSC in water (RNAse free), during 15 min and under a white light. Samples were washed twice in PBT (1 × PBS + 0.1% TritonX-100). The following protocol is a short version of the ISH protocol that it is described in detail in . For FISH all the hybridization probes were diluted to 1 ng/μl. Afterward, the samples were incubated in anti-DIG-POD 1:500, overnight at 4 °C (Sigma-Aldrich, Id_number11207733910). More than four washes, over 2–3 h, in MAB-TritonX-100 0.1% were used to eliminate the rest of the antibody. The signal was developed in TSA red 1:300, in the so-called TSA Buffer (solution of 2 M NaCl + 100 mM Borate buffer, pH 8.5), over 2–6 h. To stop the development of the signal, samples were washed in PBT. The antibody quenching was made using 1% H2O2 in PBT 0.1% for 45 min at RT. After washing 2 times with PBT, a 2nd quenching step was done with 2 × SSC, 50% Formamide and 0.1% TritonX-100 for 10 min at 56 °C. The samples were washed twice in PBT and blocked again, previous to the incubation with the second antibody: anti-DNP-HRP 1:200. After the antibody wash, with MAB-TritonX-100 0.1%, the signal was developed in TSA green 1:300 in TSA Buffer. After the double or single FISH protocols, and in the cases required, we proceeded to combine this procedure with the immunostaining of the samples (using the species-specific anti-synaptotagmin antibody; as a reference, pan-neuronal marker), as explained in the following section.
It is important to note that the general absence of well-defined tissues and the presence of nuclear intermingling in most of the “parenchymal” (external-mesodermal; internal-digestive) tissues of the Acoela makes especially difficult to perform in situ hybridization and interpret detailed patterns (as has been noted in other acoel papers).
As in the previous section, for all clones, sense probes were synthesized that were used as controls for hybridization.
Immunostaining was performed using the protocols outlined in . S. roscoffensis specimens were incubated in primary anti-synaptotagmin (dilution 1:500) antibodies (previously pre-absorbed) and reacted with the secondary antibody [Alexa Fluor goat anti-rabbit 532 (Molecular Probes, Eugene, OR)]. The anti-synaptotagmin antibody was raised in our laboratory using the specific S. roscoffensis sequence from a transcriptome analysis (see ). Preimune serum was used as control for all immunochemical experiments.
Embryo cell counting
Given the difficulty of staging the embryos in the laboratory, we relied, as a good approximation for developmental time, on the number of cells. Embryo cell counting was performed by incubation of all samples for 10 min with Dapi. A sample of 12 individuals from a pool of 12–24 h’ post-fertilization embryos were scanned completely using a Leica SP2 confocal laser microscope and their stacks were processed using the software Fiji  with the plugin “cell counter” (author: Kurt de Vos; see https://imagej.net/Cell_Counter). This procedure allowed us to obtain a precise count of the total number of cells per embryo when performing in situ procedures.
Results and discussion
bHLH genes’ developmental expression patterns
Some S. roscoffensis bHLH genes seem to share similar expression domains with their homologs in bilaterians, or metazoans, hence pointing to the possibility of a conservation of roles over evolutionary time. In the next section, we analyzed these patterns by following the grouping established for the bHLH genes. We describe here the expression patterns of 13 genes as detected by ISH, with a subsequent, more detailed, focus on those that could be involved in the nervous system development, which were thus analyzed by (the more sensitive) FISH Our descriptions are always made in the context of what is known for other metazoan members of the same groups.
Group A genes families: Achaete-Scute, E12/E47, Twist, Net, Nscl and Olig
Nine genes of the whole S. roscoffensis bHLH gene complement were classified as members of the group A: SrAscA, SrAscB, SrAsc_like, SrE12/E47, SrNet, SrNscl, SrOlig/Beta (finally classified as an Olig ortholog) and SrTwist (see  for general classification of bHLH genes or Figs. 1, 2).
The metazoan Achaete-Scute family transcription factor is divided in two subfamilies: Achaete-Scute A and B. The genes belonging to family A (also named in other clades as “proneural genes”) possess a highly conserved neurogenic role across a wide range of metazoans. It is well known that the Achaete-Scute complex members provide critical proneural function during embryogenesis and the development of adult sense organs in Drosophila melanogaster . This role is preserved in other metazoans such as the beetle Tribolium castaneun, where TcAsh is also necessary for the formation of the neural precursor , and in mice, where the gene Mash-1 is essential for the generation of autonomic and olfactory neurons. In the cnidarian Nematostella vectensis have shown that the homologous gene NvashA is specifically expressed in a differentiating subset of neural cell types of the embryonic ectoderm [47–49].
Different functions have been associated with the Achaete-Scute gene orthologs belonging to subfamily B. These genes play diverse roles, and in most cases, they are not directly involved in the development of the nervous system. In fact, most of them regulate the expression of downstream genes in different tissues; for example, the mammalian bHLH Hash-2/Mash-2/Ascl-2 are crucial for development of the placenta , the intestinal stem cell fate  or the specification of cell types in the immune system . Nevertheless, some patterns seem to be clearly associated with the neural tissue, for instance, the expression of Mash2 (a mammalian ortholog) in Schwann cells or that of the neuronal-specific planarian (S. polychroa) ortholog Spol-AscB . In our study, we found that the ISH expression pattern of SrAscB was restricted to a circular area during embryonic development (Fig. 3), while in the juveniles it was localized specifically in the half anterior region of the animal, in two bilaterally symmetric clusters of cells that were connected by a track crossing the anterior–posterior body axis (Fig. 3). Analyzing this pattern by the more sensitive FISH technique, the cell populations that express SrAscB mRNA clearly form part of the anterior neural cords (Fig. 5; Additional file 1: Fig. S1) though, in addition, they are expressed in some commissures constituting part of the brain. At the posterior end of the expression domain the pattern is shaped as a transversal band of cells, crossing the anterior–posterior body axis at the level of the mouth (also observed by FISH in the Additional file 1: Fig. S1). The domains that express SrAscB have a location slightly ventral (Additional file 1: Fig. S1). The nature of these structures (brain, cords and commissures of the nervous system) in the anterior region of the acoel juvenile suggests that SrAscB is expressed specifically in a neuronal population.
Among metazoans, the E proteins play critical roles in cell growth, specification and differentiation, including the neurons. The acoel bHLH Sr_E12/E47 gene showed a high level of expression in two circular domains of the embryo that later on became two anterior-medial regions located immediately posterior to the statocyst (Fig. 3). At first sight, our results seem to differ from those obtained for the E12/E47 gene orthologs in, for instance, H. sapiens and Drosophila (Daughterless), in which expression is found in most tissues. However, we should point out that different studies have found an increment in the mRNA expression of E12/E47 (and other E2A mRNAs) in some areas of rapid cell proliferation and differentiation in several tissues, including neural tissue. These levels decrease progressively during neurogenesis and become almost undetectable in the adult nervous system . This is consistent with our observations, with the strongest expression pattern in embryos at around 24 hpf, and lower levels in the anterior body of the juvenile (Fig. 3). The expression domain in juveniles is in an area partially overlapping the location of the brain (Fig. 3 ISH). The pattern is consistent with an early role in the development of neurons described for other organisms , putatively in populations of non-terminally differentiated neurons or neuronal progenitors. E12/E47 genes form heterodimers with other group A bHLH factors [24, 57]. In fact, some authors have suggested a specific role of E proteins in early neural differentiation  and, in the same vein, recent studies have confirmed that E proteins orchestrate neural stem cell lineage progression . bHLH genes encoding other E proteins with a similar role are seen in, for instance, the planarians, in which the gene e22/23 is expressed in the CNS . What seems clear is that both genes Sr_AscB and Sr_E12/E17 are most probably associated with the development of the nervous system.
In contrast to the mRNA expression patterns found in the above members of group A, the other members of this group (SrTwist, SrNet, SrNscl and SrOlig) showed diverse expression patterns mostly in the middle body region. SrTwist was expressed in a discrete spatial domain well delimited in the embryo’s animal pole and organized as a pair of bilateral bands on both sides of juvenile specimens. This finding is in agreement with a previous studies from our laboratory that analyzed bilaterian mesodermal gene expression have already revealed the expression pattern of twist orthologs in embryo, juvenile and adult stages of S. roscoffensis  and adults of the acoel Isodiametra pulchra , suggesting its expression in part of the gonads, the male copulatory organ (only in I. pulchra) and neoblasts, all of which are mesodermal derivatives. Furthermore, the mesodermal role of twist homologs, for instance, in Drosophila, is well known .
The expression patterns of the remaining three genes, SrNet, SrNscl and SrOlig, showed a scattered distribution throughout the body though slightly different from each other. The expression domain of SrNet gene is a bit larger than that of SrOlig and SrNscl (Fig. 3). In the case of S. roscoffensis SrNet was expressed in a region located in the middle of the embryo. In juveniles, the SrNet expression signal covered the statocyst, forming different lines in the anterior territory of the body, with a high expression in the posterior half of the body (Fig. 3). In agreement with this SrNet expression pattern, and taking into account the roles described for Net homologs in other animal models (for instance in the jellyfish Podocoryne carnea, Net is expressed in the entocodon, a mesoderm-like structure that gives rise to the striated and smooth muscle of the bell  and in Drosophila DmNet is required to maintain the inter-vein regions during development ), we suggest that the acoel Net homolog might have a mesodermal role, without any significant participation in neurogenesis, although clearly further studies are required.
Among metazoans, the contribution of Nscl gene orthologs to the development of sensory cells and neurons has been well-described [69, 70]. The S. roscoffensis bHLH gene SrNscl showed, in the juvenile, a very low expression, in a domain that seemed to be entirely posterior to the statocyst (Fig. 3). Interestingly this gene showed, however, what seemed to be a strong expression in an extensive region in our embryos (compared with the expression levels, at similar times and probe concentrations, of other bHLH). Nscl mouse’s orthologs show strong expression during the post gastrulation period, most likely due to their role in early neurogenesis . Functional analysis of the Nscl-1 mouse ortholog demonstrated the need for this gene for correct neural cell differentiation, and in combination with Nscl-2, it is required to control the migration of neuronal precursor cells [72, 73]. With the data obtained, and due to the low expression levels in juveniles, we cannot conclude with certainty that the expression of SrNscl occurs within the nervous system.
Group B gene families: Max and Srebp
The bHLHs genes belonging to group B in S. roscoffensis are SrMax and SrSrebp. The metazoan orthologs of Max are involved in cellular proliferation, development and differentiation. Several studies have shown that the MAX protein forms heterodimers with other transcription factors such as the bHLHs MYC and MAD . MYC-MAX heterodimers are involved in the transcriptional activation of different target genes [75, 76]. In S. roscoffensis, MAX is the only known member of these complexes, since MYC, MAD and other related bHLH transcription factors such as MLX were not found in the genomes or transcriptomes . This is interesting in itself since it suggests a developmental function for Max that is independent of heterodimerization. Expression of SrMax mRNA in juveniles is completely ubiquitous; however, this was not the case in embryos, where even though the signal covered a large part of the embryo an enriched expression domain was still observed (Fig. 3).
The other member of group B in our species is SREBP. Srebp orthologs are involved in animal homeostasis, in the regulation of sterol metabolism. They regulate the gene expression of most of the enzymes involved in cholesterol biosynthesis. SrSrebp expression in S. roscoffensis is widespread, at all stages analyzed (Fig. 3) consistent with having a similar (generalized) role in metabolic regulation, perhaps including supporting glial cells (see  for evidence of glial cells in S. roscoffensis).
In the case of mice and rats, the Srebp gene orthologs are also expressed in several cell types, including astrocytes, oligodendrocytes and Schwann cells, those that are very active in lipid metabolism . Moreover, in the adult planarian S. mediterranea this gene shows intense expression over the whole body, where it is expressed as well in different cell types .
Group C gene families: ARNT
The only representative of group C in our acoel species is ARNT (aryl hydrocarbon receptor nuclear translocator). ARNT can form heterodimers with several bHLH proteins, and its function depends on its dimerization partners. For this particular reason, members of the ARNT family tend to be widely expressed . Specifically, the AHR/ARNT system controls processes such as oxidation/anti-oxidation, epidermal barrier function, photo-induced response, melanogenesis, and innate immunity [33, 42] . In our genomic/transcriptomic analysis, we have found only the gene encoding for ARNT transcription factor, and not the gene for its partner AHR, contrary to what has been described in most metazoans [3, 5, 16]. SrArnt mRNA was expressed in the entire embryo with small areas inside the domain showing higher levels of expression (Fig. 3). In juveniles, probably due to the low level of signal, we detected a very faint expression, widespread in the body (Fig. 3), suggesting that ARNT in S. roscoffensis is probably expressed in a variety of cell populations.
Group E gene families: HES/HEY
Summarizing the developmental expression of neural bHLH transcription factors in S. roscoffensis
In the acoel nervous system, a similar system of neural specification seems to be in place, where a subset of the bHLH genes are used in different (and overlapping) neural domains. We are aware, though, that this is a preliminary characterization of the whole superfamily of genes and that assigning specific functions to genes or combinations of them is still premature, in absence of more detailed and/or complementary data, such as Q-PCR and in situ hybridizations on histological sections. However, this study provides us with the cartography of neural expression domains, a roadmap to further investigations.
A final caveat: The bHLH complement in the Acoela
The complement of bHLH detected in the acoel S. roscoffensis is reduced in comparison with other metazoan organisms studied. Most of the described roles of the bHLH orthologs in other organisms could not be performed with the small complement of bHLH present in S. roscoffensis. There are several different possible reasons for the detection of a reduced set of bHLHs in this clade. First, there is the obvious possibility that some bHLH appear to be missing because their sequences are too divergent and therefore, difficult to identify (see Figs. 1; 2, for a schematic representation through the different phyla). This is not unlikely in acoels, as their genomes have clearly changed a lot because of the high rate of sequence evolution. This has generated some clear genomic modifications, among which genetic losses are common. A reduction in the number of protein complements has been reported in several families previously [12, 16]. Other possible factors are the highly divergent sequences of members of other families, for instance the sequences belonging to the Wnt family of ligands identified in S. roscoffensis and the acoel Hofstenia miamia. They are clearly derived, and this has made it impossible to classify them into the well-known metazoan families [12, 83]. An alternative possibility for the scarce number of bHLH relatives in the acoel genome is that we are using a newly sequenced (draft) genome, which we assume is almost complete, but could still be missing some fragments of the genome. In this context, we should mention that the parallel use of transcriptomes, from adults and embryos, has not provided any new sequences that were not present in the genome.
The genome of S. roscoffensis contains a relatively low number of bHLH transcription factors. However, a careful phylogenetic analysis of the sequences has revealed that this clade has bHLHs from various subfamilies. Acoels possess, at least, a bHLH belonging to five of the six classical high-order groups. The conservation of the bHLH sequences analyzed is also paralleled by the conservation of the different gene expression domains. Clear correlations can be made between the expression domains in acoels and those reported for other metazoans, suggesting conservation of roles over evolutionary time. Our study detected the presence of a pool of bHLH genes (SrAscA, SrAscB, SrE12/E47 SrOlig and SrHes/Hey) with expression patterns specific to territories within the nervous system, most probably in different cell populations of neurons and/or neural precursors. Their expression, which starts in most cases during embryogenesis, suggests that they are involved in the early specification of neural precursors and the later formation of the nervous system (or subdomains thereof). This analysis, together with that of other genes studied previously in our laboratory, such as the Hox family members Cdx and SoxB, constitutes a first step toward a clear understanding of how the nervous system is assembled in acoels. The study of many other regulatory genes in our laboratory (unpublished) suggests a complex gene network controlling the development of the acoel nervous system. The further characterization and detailed knowledge of the genes involved in neural development and patterning will help to address specific developmental issues in the future; issues that would be impossible to tackle in absence of this cartography of expression domains. Having access to the sequence of the S. roscoffensis genome has been the key factor in this study.
EP-A, PM and SGS conceived the experiments. E P-A and PM wrote the paper. SGS corrected the manuscript and figures. E P-A did the experimental work, as part of her Ph.D. thesis. All authors read and approved the final manuscript.
We would like to acknowledge the support offered by Drs. Cristian Cañestro (UB; Barcelona) and Anna Aragay (IBMB-CSIC; Barcelona), which was critical in order to finish this research project. E. Perea has been a recipient of an APIF (UB) doctoral fellowship. We should also acknowledge the continuous help provided by Brenda Gavilán, Beatriz Albuixech and Carlos Herrera during most of the duration of this project. We greatly appreciate the permanent support of the members of the CLSM facility of the University of Barcelona. PM would like to acknowledge Prof. William McGinnis (UCSD), in which laboratory he learned the use of basic tools for multiplex in situ hybridization. We would like to thank the referees for their thorough, constructive, comments to this paper. They have helped us improving the manuscript.
The authors declare that none of them have any competing interests.
Availability of data and materials
Sequences and alignments are publically available through the NCBI.
Consent for publication
All authors have approved the modified version of this paper.
Ethics approval and consent to participate
We are following the rules imposed by the Universitat de Barcelona to work with animals. Our species are not vertebrates, farm or wild animals of endangered status.
Does not apply. This project was done in the absence of state or regional funding. However, we are indebted to the University of Barcelona and the support of two members of the department for their generous financial help.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Achatz JG, Martínez P. The nervous system of Isodiametra pulchra (Acoela) with a discussion on the neuroanatomy of the Xenacoelomorpha and its evolutionary implications. Front Zool. 2012;9:27. https://doi.org/10.1186/1742-9994-9-27.View ArticlePubMed CentralGoogle Scholar
- Haszprunar G. Review of data for a morphological look on Xenacoelomorpha (Bilateria incertae sedis). Org Divers Evol. 2016;16:363–89.View ArticleGoogle Scholar
- Achatz JG, Chiodin M, Salvenmoser W, Tyler S, Martinez P. The Acoela: on their kind and kinships, especially with nemertodermatids and xenoturbellids (Bilateria incertae sedis). Org Divers Evol. 2013;13:267–86.View ArticlePubMedGoogle Scholar
- Nakano H. What is Xenoturbella? Zool Lett. 2015;1:22. https://doi.org/10.1186/s40851-015-0018-z.View ArticleGoogle Scholar
- Hejnol A, Pang K. Xenacoelomorpha’s significance for understanding bilaterian evolution. Curr Opin Genet Dev. 2016;39:48–54.View ArticlePubMedGoogle Scholar
- Philippe H, Brinkmann H, Copley RR, Moroz LL, Nakano H, Poustka AJ, et al. Acoelomorph flatworms are deuterostomes related to Xenoturbella. Nature. 2011;470:255–8.View ArticlePubMed CentralPubMedGoogle Scholar
- Hejnol A, Schnabel R. What a couple of dimensions can do for you: comparative developmental studies using 4D microscopy—examples from tardigrade development. Integr Comp Biol. 2006;46:151–61.View ArticlePubMedGoogle Scholar
- Hejnol A, Lowe CJ. Embracing the comparative approach: how robust phylogenies and broader developmental sampling impacts the understanding of nervous system evolution. Philos Trans R Soc B. 2015;370:20150045. https://doi.org/10.1098/rstb.2015.0045.View ArticleGoogle Scholar
- Ruiz-Trillo I, Paps J. Acoelomorpha: earliest branching bilaterians or deuterostomes? Org Divers Evol. 2016;16:391–9.View ArticleGoogle Scholar
- Cannon JT, Vellutini BC, Smith J, Ronquist F, Jondelius U, Hejnol A. Xenacoelomorpha is the sister group to Nephrozoa. Nature. 2016;530:89–93. https://doi.org/10.1038/nature16520.View ArticleGoogle Scholar
- Perea-Atienza E, Gavilan B, Chiodin M, Abril JF, Hoff KJ, Poustka AJ, et al. The nervous system of Xenacoelomorpha: a genomic perspective. J Exp Biol. 2015;218:618–28. https://doi.org/10.1242/jeb.110379.View ArticleGoogle Scholar
- Gavilán B, Perea-Atienza E, Martínez P. Xenacoelomorpha: a case of independent nervous system centralization? Philos Trans R Soc B Biol Sci. 2016;371(1685):20150039.View ArticleGoogle Scholar
- Parker GH. The elementary nervous system. Philadelphia: JB Lippincott; 1919.View ArticleGoogle Scholar
- Northcutt RG. Evolution of centralized nervous systems: two schools of evolutionary thought. Proc Natl Acad Sci. 2012;09:10626–33.View ArticleGoogle Scholar
- Martín-Durán JM, Pang K, Børve A, Lê HS, Furu A, Cannon JT, et al. Convergent evolution of bilaterian nerve cords. Nature. 2017;553:45.View ArticlePubMedGoogle Scholar
- Perea-Atienza E, Gavilán B, Chiodin M, Abril JF, Hoff KJ, Poustka AJ, et al. The nervous system of Xenacoelomorpha: a genomic perspective. J Exp Biol. 2015;218:618–28.View ArticlePubMedGoogle Scholar
- Martinez P, Hartenstein V, Sprecher S. Xenacoelomorpha nervous systems. In: Murray Sherman S, editor. Oxford research encyclopedia of neuroscience. New York: Oxford University Press; 2017.Google Scholar
- Henry JQ, Martindale MQ, Boyer BC. The unique developmental program of the acoel flatworm, Neochildia fusca. Dev Biol. 2000;220:285–95. https://doi.org/10.1006/dbio.2000.9628.View ArticleGoogle Scholar
- Ramachandra NB, Gates RD, Ladurner P, Jacobs DK, Hartenstein V. Embryonic development in the primitive bilaterian Neochildia fusca: normal morphogenesis and isolation of POU genes Brn-1 and Brn-3. Dev Genes Evol. 2002;212:55–69.View ArticlePubMedGoogle Scholar
- Bery A, Cardona A, Martinez P, Hartenstein V. Structure of the central nervous system of a juvenile acoel, Symsagittifera roscoffensis. Dev Genes Evol. 2010;220:61–76.View ArticlePubMed CentralPubMedGoogle Scholar
- Semmler H, Chiodin M, Bailly X, Martinez P, Wanninger A. Steps towards a centralized nervous system in basal bilaterians: insights from neurogenesis of the acoel Symsagittifera roscoffensis. Dev Growth Differ. 2010;52:701–13.View ArticlePubMedGoogle Scholar
- Bery A, Martínez P. Acetylcholinesterase activity in the developing and regenerating nervous system of the acoel Symsagittifera roscoffensis. Acta Zool. 2011;92:383–92.View ArticleGoogle Scholar
- Ledent V, Vervoort M. The basic helix-loop-helix protein family: comparative genomics and phylogenetic analysis. Genome Res. 2001;11:754–70. https://doi.org/10.1101/gr.177001.View ArticlePubMed CentralGoogle Scholar
- Massari ME, Murre C. MINIREVIEW helix-loop-helix proteins: regulators of transcription in eucaryotic organisms. Mol Cell Biol. 2000;20:429–40.View ArticlePubMed CentralPubMedGoogle Scholar
- Simionato E, Ledent V, Richards G, Thomas-Chollier M, Kerner P, Coornaert D, et al. Origin and diversification of the basic helix-loop-helix gene family in metazoans: insights from comparative genomics. BMC Evol Biol. 2007;7ii:33.View ArticleGoogle Scholar
- Gyoja F. Basic helix-loop-helix transcription factors in evolution: roles in development of mesoderm and neural tissues. Genesis. 2017;55:e23051.View ArticleGoogle Scholar
- Jarman AP, Grau Y, Jan LY, Jan YN. atonal is a proneural gene that directs chordotonal organ formation in the Drosophila peripheral nervous system. Cell. 1993;73:1307–21.View ArticlePubMedGoogle Scholar
- Guillemot F, Lo L-C, Johnson JE, Auerbach A, Anderson DJ, Joyner AL. Mammalian achaete-scute homolog 1 is required for the early development of olfactory and autonomic neurons. Cell. 1993;75:463–76. https://doi.org/10.1016/0092-8674(93)90381-Y.View ArticleGoogle Scholar
- Sun Y, Nadal-Vicens M, Misono S, Lin MZ, Zubiaga A, Hua X, et al. Neurogenin promotes neurogenesis and inhibits glial differentiation by independent mechanisms. Cell. 2001;104:365–76.View ArticlePubMedGoogle Scholar
- Bertrand N, Castro DS, Guillemot F. Proneural genes and the specification of neural cell types. Nat Rev Neurosci. 2002;3:517–30. https://doi.org/10.1038/nrn874.View ArticleGoogle Scholar
- Vervoort M, Ledent V. The evolution of the neural basic helix-loop-helix proteins. Sci World J. 2001;1:396–426.View ArticleGoogle Scholar
- Gyoja F, Kawashima T, Satoh N. A genomewide survey of bHLH transcription factors in the coral Acropora digitifera identifies three novel orthologous families, pearl, amber, and peridot. Dev Genes Evol. 2012;222:63–76.View ArticlePubMedGoogle Scholar
- Cowles MW, Brown DDR, Nisperos SV, Stanley BN, Pearson BJ, Zayas RM. Genome-wide analysis of the bHLH gene family in planarians identifies factors required for adult neurogenesis and neuronal regeneration. Development. 2013;140:4691–702. https://doi.org/10.1242/dev.098616.View ArticleGoogle Scholar
- Gyoja F. A genome-wide survey of bHLH transcription factors in the Placozoan Trichoplax adhaerens reveals the ancient repertoire of this gene family in metazoan. Gene. 2014;542:29–37. https://doi.org/10.1016/j.gene.2014.03.024.View ArticleGoogle Scholar
- Fortunato SAV, Vervoort M, Adamski M, Adamska M. Conservation and divergence of bHLH genes in the calcisponge Sycon ciliatum. EvoDevo. 2016;7:1–12.View ArticleGoogle Scholar
- Sebé-Pedrós A, De Mendoza A, Lang BF, Degnan BM, Ruiz-Trillo I. Unexpected repertoire of metazoan transcription factors in the unicellular holozoan capsaspora owczarzaki. Mol Biol Evol. 2011;28:1241–54.View ArticlePubMedGoogle Scholar
- Atchley WR, Fitch WM. A natural classification of the basic helix-loop-helix class of transcription factors. Evol Contrib Walter M Fitch. 1997;94:5172–6.Google Scholar
- Jones S. An overview of the basic helix-loop-helix proteins. Genome Biol. 2004;5:226. https://doi.org/10.1186/gb-2004-5-6-226.View ArticlePubMed CentralGoogle Scholar
- Eddy SR. Accelerated profile HMM searches. PLoS Comput Biol. 2011;7:e1002195.View ArticlePubMed CentralPubMedGoogle Scholar
- Perea-Atienza E, Gavilán B, Martínez P. Immunostaining and in situ hybridization of the developing Acoel nervous system. Methods in Brain Development. Methods in Molecular Biology. Springer; 2018 (in press).Google Scholar
- Schneider CA, Rasband WS, Eliceiri KW. NIH Image to ImageJ: 25 years of image analysis. Nat Methods. 2012;9:671–5.View ArticlePubMed CentralPubMedGoogle Scholar
- Reitzel AM, Tarrant AM, Levy O. Circadian clocks in the cnidaria: environmental entrainment, molecular regulation, and organismal outputs. Integr Comp Biol. 2013;53:118–30.View ArticlePubMed CentralPubMedGoogle Scholar
- Li K-L, Lu T-M, Yu J-K. Genome-wide survey and expression analysis of the bHLH-PAS genes in the amphioxus Branchiostoma floridae reveal both conserved and diverged expression patterns between cephalochordates and vertebrates. EvoDevo. 2014;5:20. https://doi.org/10.1186/2041-9139-5-20.View ArticlePubMed CentralGoogle Scholar
- Richards GS, Rentzsch F. Transgenic analysis of a SoxB gene reveals neural progenitor cells in the cnidarian Nematostella vectensis. Development. 2014;141:4681–9.View ArticlePubMedGoogle Scholar
- Dambly-Chaudiere C, Ghysen A. Independent subpatterns of sense organs require independent genes of the achaete-scute complex in Drosophila larvae. Genes Dev. 1987;1:297–306.View ArticleGoogle Scholar
- Wheeler SR, Carrico ML, Wilson BA, Brown SJ, Skeath JB. The expression and function of the achaete-scute genes in Tribolium castaneum reveals conservation and variation in neural pattern formation and cell fate specification. Development. 2003;130:4373–81.View ArticlePubMedGoogle Scholar
- Layden MJ, Boekhout M, Martindale MQ. Nematostella vectensis achaete-scute homolog NvashA regulates embryonic ectodermal neurogenesis and represents an ancient component of the metazoan neural specification pathway. Development. 2012;139:1013–22.View ArticlePubMed CentralPubMedGoogle Scholar
- Layden MJ, Martindale MQ. Non-canonical Notch signaling represents an ancestral mechanism to regulate neural differentiation. Evodevo. 2014;5:30.View ArticlePubMed CentralPubMedGoogle Scholar
- Richards GS, Rentzsch F. Regulation of Nematostella neural progenitors by SoxB, Notch and bHLH genes. Development. 2015;142:3332–42. https://doi.org/10.1242/dev.123745.View ArticlePubMed CentralGoogle Scholar
- García-Bellido A, de Celis JF. The complex tale of the achaete-scute complex: a paradigmatic case in the analysis of gene organization and function during development. Genetics. 2009;182:631–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Koide S, Yoshida I, Tsuji A, Matsuda Y. The expression of proprotein convertase PACE4 is highly regulated by Hash-2 in placenta: possible role of placenta-specific basic helix-loop-helix transcription factor, human achaete-scute homologue-2. J Biochem. 2003;134(3):433–40. https://doi.org/10.1093/jb/mvg161.View ArticleGoogle Scholar
- Van der Flier LG, van Gijn ME, Hatzis P, Kujala P, Haegebarth A, Stange DE, et al. Transcription factor achaete scute-like 2 controls intestinal stem cell fate. Cell. 2009;136:903–12. https://doi.org/10.1016/j.cell.2009.01.031.View ArticleGoogle Scholar
- Liu X, Chen X, Zhong B, Wang A, Wang X, Chu F, et al. Transcription factor achaete-scute homologue 2 initiates follicular T-helper-cell development. Nature. 2014;507:513–8. https://doi.org/10.1038/nature12910.View ArticlePubMed CentralGoogle Scholar
- Monjo F, Romero R. Embryonic development of the nervous system in the planarian Schmidtea polychroa. Dev Biol. 2015;397:305–19. https://doi.org/10.1016/j.ydbio.2014.10.021.View ArticleGoogle Scholar
- Roberts VJ, Steenbergen R, Murre C. Localization of E2A mRNA expression in developing and adult rat tissues. Proc Natl Acad Sci USA. 1993;90:7583–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Fischer B, Azim K, Hurtado-Chong A, Ramelli S, Fernández M, Raineteau O. E-proteins orchestrate the progression of neural stem cell differentiation in the postnatal forebrain. Neural Dev. 2014;9:23. https://doi.org/10.1186/1749-8104-9-23.View ArticlePubMed CentralGoogle Scholar
- Lassar AB, Davis RL, Wright WE, Kadesch T, Murre C, Voronova A, et al. Functional activity of myogenic HLH proteins requires hetero-oligomerization with E12/E47-like proteins in vivo. Cell. 1991;66:305–15. https://doi.org/10.1016/0092-8674(91)90620-E.View ArticleGoogle Scholar
- Heng JIT, Tan SS. The role of class I HLH genes in neural development—Have they been overlooked? BioEssays. 2003;25:709–16.View ArticleGoogle Scholar
- Chiodin M. The evolution of bilaterian body-plan: perspectives from the developmental genetics of the Acoela. PhD Thesis. MARTA CHIODIN. Universitat de Barcelona, Spain; 2013.Google Scholar
- Chiodin M, Børve A, Berezikov E, Ladurner P, Martinez P, Hejnol A. Mesodermal gene expression in the acoel Isodiametra pulchra indicates a low number of mesodermal cell types and the endomesodermal origin of the gonads. PLoS ONE. 2013;8:2.View ArticleGoogle Scholar
- Baylies MK, Bate M. Twist: a myogenic switch in Drosophila. Science. 1996; 272:1481–84. http://science.sciencemag.org/content/272/5267/1481.abstract.
- Müller P, Seipel K, Yanze N, Reber-Müller S, Streitwolf-Engel R, Stierwald M, et al. Evolutionary aspects of developmentally regulated helix-loop-helix transcription factors in striated muscle of jellyfish. Dev Biol. 2003;255:216–29.View ArticlePubMedGoogle Scholar
- Brentrup D, Lerch H, Jackle H, Noll M. Regulation of Drosophila wing vein patterning: net encodes a bHLH protein repressing rhomboid and is repressed by rhomboid-dependent Egfr signalling. Development. 2000;127:4729–41. http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=PubMed&dopt=Citation&list_uids=11023875.
- Mizuguchi R, Sugimori M, Takebayashi H, Kosako H, Nagao M, Yoshida S, et al. Combinatorial roles of Olig2 and Neurogenin2 in the coordinated induction of pan-neuronal and subtype-specific properties of motoneurons. Neuron. 2001;31:757–71.View ArticlePubMedGoogle Scholar
- Novitch BG, Chen AI, Jessell TM. Coordinate regulation of motor neuron subtype identity and pan-neuronal properties by the bHLH repressor Olig2. Neuron. 2001;31:773–89.View ArticlePubMedGoogle Scholar
- Lu QR, Sun T, Zhu Z, Ma N, Garcia M, Stiles CD, et al. Common developmental requirement for Olig function indicates a motor neuron/oligodendrocyte connection. Cell. 2002;109:75–86.View ArticlePubMedGoogle Scholar
- Zhou Q, Anderson DJ. The bHLH transcription factors OLIG2 and OLIG1 couple neuronal and glial subtype specification. Cell. 2002;109:61–73.View ArticlePubMedGoogle Scholar
- Sugimori M, Nagao M, Bertrand N, Parras CM, Guillemot F, Nakafuku M. Combinatorial actions of patterning and HLH transcription factors in the spatiotemporal control of neurogenesis and gliogenesis in the developing spinal cord. Development. 2007;134:1617–29.View ArticlePubMedGoogle Scholar
- Begley CG, Lipkowitz S, Göbel V, Mahon KA, Bertness V, Green AR, et al. Molecular characterization of NSCL, a gene encoding a helix-loop-helix protein expressed in the developing nervous system. Proc Natl Acad Sci USA. 1992;89:38–42.View ArticlePubMed CentralPubMedGoogle Scholar
- Kim WY. NeuroD1 is an upstream regulator of NSCL1. Biochem Biophys Res Commun. 2012;419:27–31. https://doi.org/10.1016/j.bbrc.2012.01.100.View ArticleGoogle Scholar
- Gobel V, Lipkowitz S, Kozak CA, Kirsch IR. Nscl-2—a basic domain helix-loop-helix gene expressed in early neurogenesis. Cell Growth Differ. 1992;3:143–8.PubMedGoogle Scholar
- Krüger M, Ruschke K, Braun T. NSCL-1 and NSCL-2 synergistically determine the fate of GnRH-1 neurons and control necdin gene expression. EMBO J. 2004;23:4353–64.View ArticlePubMed CentralPubMedGoogle Scholar
- Schmid T, Kruger M, Braun T. NSCL-1 and -2 control the formation of precerebellar nuclei by orchestrating the migration of neuronal precursor cells. J Neurochem. 2007;102:2061–72.View ArticlePubMedGoogle Scholar
- Gallant P. Myc/Max/Mad in invertebrates: the evolution of the max network. Curr Top Microbiol Immunol. 2006;302:235–53.PubMedGoogle Scholar
- Eilers M, Eisenman RN. Myc’s broad reach. Genes Dev. 2008;22:2755–66. https://doi.org/10.1101/gad.1712408.View ArticlePubMed CentralGoogle Scholar
- Hartl M, Mitterstiller A-M, Valovka T, Breuker K, Hobmayer B, Bister K. Stem cell-specific activation of an ancestral myc protooncogene with conserved basic functions in the early metazoan Hydra. Proc Natl Acad Sci USA. 2010;107:4051–6.View ArticlePubMed CentralPubMedGoogle Scholar
- Camargo N, Smit AB, Verheijen MHG. SREBPs: SREBP function in glia-neuron interactions. FEBS J. 2009;276:628–36.View ArticlePubMedGoogle Scholar
- McIntosh BE, Hogenesch JB, Bradfield CA. Mammalian per-arnt-sim proteins in environmental adaptation. Ann Rev Physiol. 2010;72:625–45.View ArticleGoogle Scholar
- Denison MS, Soshilov AA, He G, Degroot DE, Zhao B. Exactly the same but different: promiscuity and diversity in the molecular mechanisms of action of the aryl hydrocarbon (dioxin) receptor. Toxicol Sci. 2011;124:1–22.View ArticlePubMed CentralPubMedGoogle Scholar
- Guillemot F. Spatial and temporal specification of neural fates by transcription factor codes. Development. 2007;134:3771–80. https://doi.org/10.1242/dev.006379.View ArticleGoogle Scholar
- Ishibashi M, Ang S, Shiota K, Nakanishi S, Kageyama R, Guillemot F. Targeted disruption of mammalian hairy and Enhancer of split homolog-1 (HES-1) leads to up regulation of neural hehx-loop-hellx factors, premature neurogenesls, and severe neural tube defects. Genes Dev. 1995;1:3136–48.View ArticleGoogle Scholar
- Kageyama R, Ohtsuka T, Kobayashi T. Roles of Hes genes in neural development. Dev Growth Differ. 2008;50(SUPPL):1.Google Scholar
- Srivastava M, Mazza-Curll KL, Van Wolfswinkel JC, Reddien PW. Whole-body acoel regeneration is controlled by Wnt and Bmp-Admp signaling. Curr Biol. 2014;24:1107–13. https://doi.org/10.1016/j.cub.2014.03.042.View ArticleGoogle Scholar