Ancient origins of arthropod moulting pathway components
Abstract
Ecdysis (moulting) is the defining character of Ecdysoza (arthropods, nematodes and related phyla). Despite superficial similarities, the signalling cascade underlying moulting differs between Panarthropoda and the remaining ecdysozoans. Here, we reconstruct the evolution of major components of the ecdysis pathway. Its key elements evolved much earlier than previously thought and are present in non-moulting lophotrochozoans and deuterostomes. Eclosion hormone (EH) and bursicon originated prior to the cnidarian-bilaterian split, whereas ecdysis-triggering hormone (ETH) and crustacean cardioactive peptide (CCAP) evolved in the bilaterian last common ancestor (LCA). Identification of EH, CCAP and bursicon in Onychophora and EH, ETH and CCAP in Tardigrada suggests that the pathway was present in the panarthropod LCA. Trunk, an ancient extracellular signalling molecule and a well-established paralog of the insect peptide prothoracicotropic hormone (PTTH), is present in the non-bilaterian ctenophore Mnemiopsis leidyi. This constitutes the first case of a ctenophore signalling peptide with homology to a neuropeptide.
https://doi.org/10.7554/eLife.46113.001eLife digest
Animals such as insects, crabs and spiders belong to one of the most species-rich animal groups, called the arthropods. These animals have exoskeletons, which are hard, external coverings that support their bodies. Arthropods shed their exoskeletons as they grow, a process called ecdysis or moulting, and this behaviour is controlled by a set of hormones and small protein-like molecules called neuropeptides that allow communication between neurons.
Other animals, such as roundworms, also moult; and together with arthropods they are classified into a group called the Ecdysozoa. Since moulting is a common behaviour in ecdysozoans, it was previously assumed that its signalling components had evolved in the common ancestor of roundworms and arthropods, although differences in the moulting machinery between both groups exist.
Here, De Oliveira et al. investigate the evolutionary origins of the arthropod moulting machinery and find that some of the hormones and neuropeptides involved appeared long before the arthropods themselves.
Database searches showed that important hormones and neuropeptides involved in arthropod moulting can be found in diverse animal groups, such as jellyfish, molluscs and starfish, confirming that these molecules evolved before the last common ancestor of roundworms and arthropods. These animals must therefore use the hormones and neuropeptides in many processes unrelated to moulting. De Oliveira et al. also found that roundworms have lost most of these molecules, and that moulting in these animals must be driven by a different complement of hormones and neuropeptides.
These results invite research into the role of moulting hormones and neuropeptides in animals outside the Ecdysozoa. They also show that signalling pathways and the processes they regulate are highly adaptable: two animals can use the same hormone in entirely different processes, but conversely, the same behaviour may be regulated by different molecules depending on the animal. This means that the evolution of a process and the evolution of its regulation can be decoupled, a finding that has important implications for the study of signalling pathways and their evolution.
https://doi.org/10.7554/eLife.46113.002Introduction
Ecdysis or moulting, which describes the process of shedding the outer integument, the cuticle, is a defining feature of Ecdysozoa (arthropods, tardigrades, onychophorans, nematodes and related phyla) (Aguinaldo et al., 1997; Schmidt-Rhaesa et al., 1998; de Rosa et al., 1999; Dunn et al., 2008; Telford et al., 2008). Despite superficial similarities of the ‘moulting behaviour’ within Ecdysozoa, the neuroendocrine components underlying this process remain elusive for the majority of the ecdysozoans outside of Arthropoda. This includes well-established model organisms such as the nematode Caenorhabditis elegans, for which the gene regulatory network responsible for ecdysis remains to be fully resolved (Frand et al., 2005; reviewed by Page et al., 2014 and Lažetić and Fay, 2017).
In arthropods, ecdysis can be divided into three distinct stages, pre-ecdysis, ecdysis and post-ecdysis. Each of these stages correlates with major behavioural, molecular and cellular changes and encompasses a series of specific muscular contractions controlled by a cascade of hormones and neuropeptides (Truman, 2005). Studies in insects have revealed that the major components of this peptidergic signalling pathway are ecdysis-triggering hormone (ETH), eclosion hormone (EH), crustacean cardioactive peptide (CCAP) and bursicon (Gammie and Truman, 1997a; Gammie and Truman, 1997b; Zitnan et al., 1999; Clark et al., 2004; Kim et al., 2006a; Kim et al., 2006b; Arakane et al., 2008; Lee et al., 2013). The process begins with the release of prothoracicotropic hormone (PTTH) from neurohemal organs. PTTH initiates a signalling cascade that results in the biosynthesis of ecdysteroids (i.e., steroid hormones synthesised from ingested cholesterol), including ecdysone (E) and 20-hydroxyecdysone (20E) (Figure 1). The decline of the ecdysone titre due to the ecdysone-inactivating enzyme cytochrome P450 protein Cyp18a1 (Guittard et al., 2011; reviewed by Rewitz et al., 2013) triggers the release of ETH that, in turn, causes the release of EH. These two hormones mutually enhance one another in a positive feedback loop to control and regulate pre-ecdysis behaviour (Figure 1). With the ensuing release of CCAP, caused by EH, pre-ecdysis ceases and the ecdysis motor program is initiated. Finally, bursicon responds to the increasing levels of CCAP and initiates post-ecdysis behaviour and cuticle tanning (Figure 1).
Comparative biochemical, genomic and transcriptomic analyses revealed that ecdysteroids and the required genes responsible for their biosynthesis are present outside of Ecdysozoa, showing that some key molecular players of moulting predate the origin of Ecdysozoa (Mendis et al., 1984; Nolte et al., 1986; Garcia et al., 1989; Barker et al., 1990; Schumann et al., 2018). Such integrative and comparative analyses have so far not been conducted on the major components of the peptidergic signalling system underlying moulting. To fill this gap in knowledge, we explored the distribution of PTTH, ETH, EH, CCAP and bursicon ligand-receptor pairs across Metazoa.
Results and discussion
PTTH is a neurohormone with a proposed origin at the base of Arthropoda that is believed to have evolved from the duplication of the ancient and widely distributed bilaterian signalling molecule-encoding gene trunk (Rewitz et al., 2009; Jékely, 2013). By screening 39 metazoan genomes and 57 transcriptomes (Supplementary file 1), we found that the PTTH peptide is present in Drosophila and Tribolium but absent in the house spider Parasteatoda tepidariorum and the crustacean Parhyale hawaiensis, suggesting that PTTH is an insect innovation (Figure 2A,B, Figure 2—figure supplement 1A, Figure 2—source data 1). The ablation of PTTH-producing neurons in Drosophila generates an imbalance in ecdysone biosynthesis, causing developmental delay, prolonged duration of feeding and larger individuals with reduced fecundity (McBrayer et al., 2007). These findings indicate that PTTH, at least in Drosophila, regulates developmental timing and body size, but is not essential for moulting. Trunk, the paralog of ptth, has previously been identified in arthropods, annelids, mollusks and the cephalochordate Branchiostoma floridae (Jékely, 2013). Our study expands the phyletic distribution of trunk to onychophorans, tardigrades, gastrotrichs, brachiopods, nemerteans, ectoprocts, phoronids and hemichordates (Figure 2A,B, Figure 2—figure supplement 1A, Figure 2—source data 1). Although we did not find a trunk ortholog in the genome of the sea anemone Nematostella vectensis, our similarity searches against the NCBI protein database led to the identification of this protein in other anthozoans, namely Stylophora pistillata (PFX31008.1) and Orbicella faveolata (XP_020630744.1 and XP_020630745.1). More importantly, our multi-species screen also recovered a trunk-like peptide in the ctenophore Mnemiopsis leidyi with high similarity (p-value < 1e-05) to lophotrochozoan, deuterostome and ecdysozoan trunk sequences (Figure 2—figure supplement 1A; Figure 3A,B). By similarity-based clustering we were able to demonstrate homology of the ctenophore trunk-like peptide with the insect trunk paralog, ptth (Figure 3A, Figure 3—source data 1; see also Rewitz et al., 2009; Jékely, 2013). This extends the phyletic distribution of trunk to the ctenophores (see, e.g., Halanych, 2004; Dunn et al., 2008; Moroz et al., 2014; Jékely et al., 2015; Pisani et al., 2015 for discussion).
-
Figure 2—source data 1
- https://doi.org/10.7554/eLife.46113.011
-
Figure 2—source data 2
- https://doi.org/10.7554/eLife.46113.012
-
Figure 2—source data 3
- https://doi.org/10.7554/eLife.46113.013
-
Figure 2 - source data 4
- https://doi.org/10.7554/eLife.46113.014
-
Figure 2—source data 5
- https://doi.org/10.7554/eLife.46113.015
-
Figure 3—source data 1
- https://doi.org/10.7554/eLife.46113.017
PTTH and trunk share a common receptor, the tyrosine kinase torso (Rewitz et al., 2009). Similar to its ligands, torso (Rewitz et al., 2009) proved also to be much more ancient than commonly assumed. We identified torso sequences in deuterostomes, lophotrochozoans, cnidarians and ecdysozoans (Figure 2—figure supplement 2), indicating that the trunk-torso neuropeptide signalling pathway dates back at least as far as the last common ancestor of Cnidaria, Ctenophora and Bilateria and is thus not restricted to Bilateria as suggested previously (e.g., Jékely, 2013) (Figures 2A,B and 4A).
In insects, the first hormone released in response to decreasing ecdysone levels is usually ETH (Zitnan et al., 1996; Zitnan et al., 1999) although in the lepidopteran Manduca sexta the neuropeptide corozanin acts as the trigger for the release of ETH from the epitracheal glands (Kim et al., 2004). Knockdown of the eth gene in Drosophila (Park et al., 2002) and of eth and its receptors in Tribolium and Schistocerca (Arakane et al., 2008; Lenaerts et al., 2017) lead to lethality at the expected onset of ecdysis, demonstrating the essential role of the ETH peptidergic signalling system in moulting (Park et al., 2002; Arakane et al., 2008; Lenaerts et al., 2017; Shi et al., 2017). Our screening and phylogenetic analyses confirmed the presence of ETH and its receptor in tardigrades and arthropods, thus corroborating previous studies (Figure 2B, Figure 2—figure supplement 3, Figure 2—source data 2) (Zitnan et al., 1996; Zitnan et al., 1999; Park et al., 2002; Arakane et al., 2008; Veenstra et al., 2012; Lenaerts et al., 2017; Koziol, 2018; Zhu et al., 2019). For Arthropoda, the ETH ligand was only found in insects (Drosophila and Tribolium), but was lacking in the crustacean Parhyale hawaiensis and the arachnid Parasteatoda tepidariorum. However, studies on the two mites Panonychus citri and Tetranychus urticae as well as several decapods have shown the presence of the ETH ligand in chelicerates and crustaceans (in which the homology was reconfirmed by our clustering analysis) (Veenstra et al., 2012; Veenstra, 2016; Zhu et al., 2019). Surprisingly, we did not find the entire ETH signalling pathway in the two onychophoran genomes analysed herein (Figures 2B and 4B, Figure 2—figure supplement 3, Figure 2—source data 2). Due to the fragmented nature of the onychophoran genomes a final statement whether this signalling system was indeed lost in this lineage cannot be made at present.
No eth ortholog was identified outside the panarthropods, including the nematode C. elegans, suggesting that this gene originated in the last common ancestor of Panarthropoda (Figure 2B). Our findings on the distribution of the ETH receptor are in agreement with the results of previous studies and demonstrate the presence of this receptor in arthropods (insects, crustaceans and arachnids), mollusks, nemerteans, brachiopods, echinoderms, cephalochordates (in which we found a substantial expansion of eth-receptor homologs in the genomes of Branchiostoma floridae and B. belcheri), vertebrates and acoels (Park et al., 2003; Roller et al., 2010; Veenstra et al., 2012; Mirabeau and Joly, 2013; Lenaerts et al., 2017; Thiel et al., 2018; Zhu et al., 2019) (Figures 2B and 4A, Figure 2—figure supplement 3, Figure 2—source data 2). This provides evidence for the presence of individual components of the ETH signalling system already at the base of Bilateria (Figures 2A and 4A).
Eclosion hormone (EH) was first identified as a blood-borne factor (Truman and Riddiford, 1970) in three lepidopteran species, Hyalophora cecropia, Antheraea polyphemus and Antheraea pernyi. A positive feedback loop between EH and ETH was found in Manduca and suggested in Tribolium (Ewer et al., 1997; Arakane et al., 2008), whereas in Drosophila, EH has been described either acting downstream of ETH (Kim et al., 2006a; Kim et al., 2006b) or in a positive endocrine feedback loop with ETH (Krüger et al., 2015). Despite EH being a key regulator of ecdysis in insects, eh knockout in Drosophila melanogaster did not abolish ecdysis, but instead produced flies with discrete behavioural deficits such as slow and uncoordinated eclosion. This shows that EH is involved, but does not play an essential role, in moulting in these insects (McNabb et al., 1997). This result, however, has been contested in a more recent study using eh null mutants in Drosophila (Krüger et al., 2015), showing that the lack of eh function is lethal during larval fruit fly ecdysis.
Traditionally considered to be confined to arthropods, recent studies showed the presence of EH and its receptor, a guanylyl cyclase, in echinoderms and tardigrades (Zandawala et al., 2017; Koziol, 2018). Our study corroborates these findings but considerably expands the presence of the EH ligand to cnidarians, acoels, hemichordates, lophotrochozoans (mollusks, annelids, nemerteans and phoronids) and onychophorans (Figure 2B, Figure 2—figure supplement 1B, Figure 2—source data 3). All EH ligand orthologs harbour the six cysteine conserved residues (Zitnan et al., 2007) except for Cnidaria, in which only five are present. The identification of the EH receptor in ambulacrarians, mollusks, annelids, nemerteans and phoronids suggests co-evolution of this ligand-receptor pair throughout Metazoa (Figure 2B, Figure 2—figure supplement 4, Figure 2—source data 3). Although the eh-receptor gene was not found in Cnidaria, Xenacoelomorpha and Onychophora, its distribution includes the Brachiopoda, Ectoprocta and Cephalochordata lineages. Consequently, our findings shift the ancestry of this peptidergic pathway back to the cnidarian-bilaterian split (Figures 2A and 4A).
First isolated from the shore crab Carcinus maenas, crustacean cardioactive peptide (CCAP) is a highly conserved amidated neuropeptide that increases heart rate in crustaceans and insects (Stangier et al., 1987; Cheung et al., 1992; Lehman et al., 1993; Suggs et al., 2016). CCAP has multiple functions in addition to its cardioacceleratory activity, such as accelerating the frequency and amplitude of oviduct contractions in the locust Locusta migratoria (Donini et al., 2001) and regulating the release of digestive enzymes in the cockroach Periplaneta americana (Sakai et al., 2006). CCAP is important for ecdysis in crustaceans and insects where it initiates the stereotyped sequence of behaviours that mark the end of the pre-ecdysis stage (Gammie and Truman, 1997a; Gammie and Truman, 1997b; Phlippen et al., 2000; Arakane et al., 2008; Lee et al., 2013). However, transgenic Drosophila larvae lacking CCAP neurons moult normally and only exhibit a prolonged pre-ecdysis behaviour (Clark et al., 2004).
Only three studies focusing on the CCAP signalling pathway components are available outside of Arthropoda. In the snail Lymnaea stagnalis (Vehovszky et al., 2005), immunostaining revealed a dense network of CCAP-positive fibres that likely function to regulate parts of the feeding behaviour. In the oyster Saccostrea glomerata (In et al., 2016) and in the cuttlefish Sepia officinalis (Endress et al., 2018), in vivo bioassays using synthesised neuropeptides and immunohistochemistry suggested that the CCAP signalling pathway is involved in reproduction (e.g., spawning, oocyte transport, egg-laying). Additionally, Sepia CCAP has been shown to increase the tonus of the vena cava, demonstrating its role in the regulation of hemolymph circulation (Endress et al., 2018). These results indicate that in both, mollusks and arthropods, CCAP functions in feeding, reproduction and regulation of hemolymph circulation, suggesting that these may have been its ancestral roles. In arthropods, co-option of CCAP into the ecdysis pathway expanded this set of functions to include moulting.
The CCAP receptor is a G protein-coupled receptor (GPCR) that was first described from the Drosophila genome and subsequently identified in many other insects (Cazzamali et al., 2003; Arakane et al., 2008; Vogel et al., 2013). We confirm here the presence of a CCAP ligand in mollusks, annelids, arthropods and tardigrades, as stated earlier (Veenstra, 2010; Jékely, 2013; Mirabeau and Joly, 2013; Conzelmann et al., 2013; Stewart et al., 2014; Ahn et al., 2017; Zhang et al., 2018; Koziol, 2018), and extend the distribution of the CCAP ligand to three additional lophotrochozoan phyla (Nemertea, Platyhelminthes and Rotifera) as well as to the remaining panarthropod phylum Onychophora (Figures 2B and 4B, Figure 2—figure supplement 1C, Figure 2—source data 4). Interestingly, the CCAP ligand is absent from all investigated deuterostome genomes analysed (Figure 2—figure supplement 1C, Figure 2B, Figure 2—source data 4). The ccap-receptor ortholog was found in acoels, lophotrochozoans and panarthropods except Onychophora (Figure 2—figure supplement 5). Surprisingly, we found the receptor also in all deuterostome phyla (Echinodermata, Hemichordata, Vertebrata, Cephalochordata) except Tunicata (Figure 2—figure supplement 5, Figure 2—source data 4). These results reinforce the suggested origin of the ligand-receptor pair at the base of Bilateria and points towards a possible loss of the CCAP ligand in the Deuterostomia lineage (Figures 2B and 4A, Figure 2—figure supplement 5, Figure 2—source data 4).
Bursicon was identified as a neurohormone responsible for cuticle sclerotization and melanisation (tanning) during post-ecdysis (Cottrell, 1962a; Cottrell, 1962b; Fraenkel and Hsiao, 1965). Recent studies have shown that bursicon also has a mild effect on the regulation of pre-ecdysis and is important for the proper execution of post-ecdysis in Manduca, Drosophila and Tribolium as well as for the development of wings and other integumentary structures (Baker and Truman, 2002; Dewey et al., 2004; Arakane et al., 2008; Bai and Palli, 2010). Together with the ecdydis-triggering hormone signalling system, bursicon is an indispensable component of the moulting behaviour in insects (Arakane et al., 2008). Previous studies show that bursicon is present outside Ecdysozoa, for example in the anthozoan Nematostella vectensis, the echinoderm Strongylocentrotus purpuratus as well as in annelids and mollusks (Jékely, 2013; Conzelmann et al., 2013; Stewart et al., 2014; Ahn et al., 2017; Zhang et al., 2018). Our work confirms the presence of the complete bursicon signalling system in all arthropod genomes analysed here and extends its distribution (receptor and/or ligand) to the hemichordate, nemertean, phoronid, rotifer and onychophoran phyla (Figure 2B, Figure 2—figure supplement 1D, Figure 2—source data 5). Interestingly, bursicon and its receptor rickets are absent in tardigrades, suggesting the loss of the bursicon peptidergic signalling in this lineage (Figures 2B and 4B, Figure 2—figure supplement 6, Figure 2—source data 5). The latter findings are corroborated by independent proneuropeptide and peptide prohormone surveys in the tardigrade genomic and EST data that also failed to detect this ligand-receptor pair in different tardigrade species (Christie et al., 2011; Koziol, 2018). We did not identify the receptor in any deuterostome or non-arthropod ecdysozoan lineage (Figure 2B).
The PTTH, ETH, EH, CCAP and bursicon peptide signalling systems are lacking in the nematode Caenorhabditis elegans (cf. our study and, e.g., Page et al., 2014; Lažetić and Fay, 2017; Figure 2B, Figure 2—figure supplement 1). Additionally, other key moulting components, such as the ecdysteroid ecdysone (E), 20-hydroxyecdysone (20E), and various halloween gene products have also been reported absent from the C. elegans genome (Frand et al., 2005; Schumann et al., 2018). An extensive body of research on moulting in C. elegans suggests an entirely different molecular machinery controlling this behaviour in this free-living nematode (Russel et al., 2011; for review Lažetić and Fay, 2017).
Interestingly, however, E and 20E were identified in parasitic nematodes (Cleator et al., 1987; Shea et al., 2004) and, outside Ecdysozoa, in the platyhelminth Monieza expansa, the gastropod mollusks Lymnaea stagnalis and Helix pomatia as well as in the hirudinean annelid Hirudo medicinalis (Mendis et al., 1984; Nolte et al., 1986; Garcia et al., 1989; Barker et al., 1990).
Conclusion
We show that key peptidergic components of the arthropod ecdysis pathway emerged prior to the protostome-deuterostome split, and thus considerably earlier than commonly assumed. EH, CCAP and the bursicon signalling systems are more widespread among non-moulting animals than previously appreciated. The presence of the eth-receptor ortholog in ecdysozoans, lophotrochozoans and deuterostomes, in combination with the restriction of its known ligand to insects, arachnids and tardigrades, suggests a scenario in which promiscuous ligand/receptor relationships can lead to novel signalling interactions that provide new opportunities for natural selection to generate novel functions (Figure 2B). The identification of the near complete suite of the peptidergic arthropod ecdysis pathway components in Onychophora and Tardigrada strongly suggests that the entire pathway was at least functional in the last common ancestor of Panarthropoda and maybe as early as in the ur-ecdysozoan (Figures 2B and 4B). However, considering the crucial role of the ETH and bursicon signalling systems in insect moulting, together with the apparent secondary loss of ETH in Onychophora and bursicon in Tardigrada (Figure 4B), the consequences of harbouring only the partial set of the ecdysis signalling genes should be the focus of future assessments. Independent recruitment of novel peptidergic components into insect ecdysis has been shown (cf. Kim et al., 2004; Kim et al., 2006a; Kim et al., 2006b; extensively reviewed by Zitnan et al., 1996), illustrating the evolutionary plasticity of this signalling pathway and calling for more detailed functional investigations into the role of individual components during moulting of the various ecdysozoan lineages.
Materials and methods
Data collection, filtering, sequence reconstruction and proteome prediction
Request a detailed protocolTo obtain a comprehensive sampling across Metazoa, ecdysozoan, deuterostome and non-bilaterian protein-coding sequence (CDS) databases were downloaded from publicly available sites and combined with previous lophotrochozoan transcriptomes (see De Oliveira et al., 2019). The acoel transcriptomic data were pre-processed and assembled as described in De Oliveira et al. (2019). The databases include representatives from the following phyla: Porifera, Ctenophora, Cnidaria, Placozoa, Xenacoelomorpha, Echinodermata, Hemichordata, Chordata, Annelida, Brachiopoda, Ectoprocta, Entoprocta, Gastrotricha, Mollusca, Nemertea, Phoronida, Platyhelminthes, Rotifera, Arthropoda, Tardigrada, Onychophora and Nematoda. The choanoflagellate Monosiga brevicollis was used as outgroup. Supplementary file 1 summarises the databases and the publicly available repositories from which they were obtained. Sequence read archive (SRA) accession numbers for xenacoelomorph databases are also shown.
Sensitive similarity searches with jackhmmer
Request a detailed protocolSensitive probabilistic iterative similarity searches based on profile hidden Markov models (HMMs) were performed with jackhmmer (Johnson et al., 2010) against the respective metazoan and choanoflagellate databases. Insect eh, eth ccap, ptth and bursicon orthologs were retrieved from NCBI (National Center for Biotechnology Information) and their respective receptors from Vogel et al. (2013). These sequences were used as queries in the similarity searches. The searches were performed under the default parameters using varying e-value thresholds (1 to 1e-06) controlled by the options –E and –domE, as defined in jackhmmer. The best hits found in the metazoan and choanoflagellate databases were stored in fasta format and used in the subsequent analyses.
Clustering and phylogenetic analyses
Request a detailed protocolEH, ETH, CCAP, PTTH and bursicon ligand candidates retrieved from the metazoan and choanoflagellate databases were used as input, together with their respective insect orthologs, in the program clans (Frickey and Lupas, 2004) under different e-value thresholds (0.1 to 1e-06) and blast programs, that is blastp or psiblast (Camacho et al., 2009). Singleton sequences (isolated unconnected sequences) were excluded from the map. To further improve the orthology assessment, multiple sequence alignments were performed with mafft (Katoh and Standley, 2013) and the presence of shared conserved amino acid regions and residues were investigated with aliview (Larsson, 2014). The final 3D maps were collapsed into 2D after the clustering for easier visualisation.
Putative EH, ETH, CCAP, PTTH and bursicon receptor candidates retrieved from the metazoan and choanoflagellate databases were aligned with mafft together with their respective orthologs, when found, and subsequently trimmed with BMGE software under the following parameters: –h 1 –b 1 –m BLOSUM30 –t AA (Criscuolo and Gribaldo, 2010). Outgroups for the phylogenetic analyses were defined according to Vogel et al. (2013).
Phylogenetic analyses were performed using RAxML (Stamatakis, 2014), PhyML (Guindon et al., 2010) and mrbayes (Ronquist et al., 2012) softwares using the appropriate best-fit model of amino acid substitution. RaxML was executed under default parameters and rapid bootstrap. PhyML was executed under the default parameters and an optimised starting tree (-o tlr option). The number of bootstrap values was set to 1.000 in RaxML and PhyML and the number of generations used in mrbayes was determined using a convergence diagnostic. All runs in mrbayes were performed with the samplefreq and relative burn-in defined as 1000 and 25%, respectively. The three final phylogenetic trees obtained for each of the four different receptors were visualised and combined with TreeGraph2 (Stöver and Müller, 2010).
Data availability
Request a detailed protocolAll data generated in the course of this study are included in this article (Figure 2—source datas 1–5 and Figure 3—source data 1). The accession numbers for the publicly available datasets used in this work are available in Supplementary file 1. The 3D cluster peptide maps can be visualised and manipulated using the program clans (Frickey and Lupas, 2004); see ftp://ftp.tuebingen.mpg.de/pub/protevo/CLANS/). The multiple sequence alignment files can be viewed with aliview (Larsson, 2014). The phylogenetic tree files can be viewed using Figtree (http://tree.bio.ed.ac.uk/software/figtree/) or TreeGraph2 (Stöver and Müller, 2010).
Data availability
All data generated or analysed during this study are included in the manuscript and supporting files. Source data files have also been provided. The molecular databases analysed in this study are publicly available and novel annotated sequences are included in the Source data 1-5. The list of investigated species and their respective links to a direct download are presented in the Supplementary File 1 (Table S1). The 3D proneuropeptide/prohormone maps as well as all the multiple sequence alignments and the phylogenetic trees generated in this study are available in the Source Data 1-5 enclosed in the original submission. The 3D maps in .rtf format can be visualised and inspected with the software clans (ftp://ftp.tuebingen.mpg.de/pub/protevo/CLANS/). The multiple sequence alignments used in the phylogenetic inferences can be graphically visualised using aliview (http://www.ormbunkar.se/aliview/#DOWNLOAD). The phylogenetic tree files can be viewed using an appropriate phylogetic tree viewer such as Figtree (http://tree.bio.ed.ac.uk/software/figtree/).
References
-
Mutations in the Drosophila glycoprotein hormone receptor, rickets, eliminate neuropeptide-induced tanning and selectively block a stereotyped behavioral programThe Journal of Experimental Biology 205:2555–2565.
-
Ecdysteroids in helminths and annelidsInvertebrate Reproduction & Development 18:1–11.https://doi.org/10.1080/07924259.1990.9672124
-
BLAST+: architecture and applicationsBMC Bioinformatics 10:421.https://doi.org/10.1186/1471-2105-10-421
-
Molecular identification of a Drosophila G protein-coupled receptor specific for crustacean cardioactive peptideBiochemical and Biophysical Research Communications 303:146–152.https://doi.org/10.1016/S0006-291X(03)00302-4
-
Bioinformatic prediction of arthropod/nematode-like peptides in non-arthropod, non-nematode members of the ecdysozoaGeneral and Comparative Endocrinology 170:480–486.https://doi.org/10.1016/j.ygcen.2010.11.002
-
Identity and tissue localization of free and conjugated ecdysteroids in adults of Dirofilaria immitis and Ascaris suumMolecular and Biochemical Parasitology 25:93–105.https://doi.org/10.1016/0166-6851(87)90022-3
-
The imaginal ecdysis of blowflies. The control of cuticuar hardening and darkeningThe Journal of Experimental Biology 39:395–412.
-
The imaginal ecdysis of blowflies. Detection of the blood-borne darkening factor and determination of some of its propertiesThe Journal of Experimental Biology 39:413–430.
-
Crustacean cardioactive peptide is a modulator of oviduct contractions in Locusta migratoriaJournal of Insect Physiology 47:277–285.https://doi.org/10.1016/S0022-1910(00)00112-8
-
Animal phylogeny and its evolutionary implicationsAnnual Review of Ecology, Evolution, and Systematics 45:371–395.https://doi.org/10.1146/annurev-ecolsys-120213-091627
-
Control of insect ecdysis by a positive-feedback endocrine system: roles of eclosion hormone and ecdysis triggering hormoneThe Journal of Experimental Biology 200:869–881.
-
Bursicon, a hormone which mediates tanning of the cuticle in the adult fly and other insectsJournal of Insect Physiology 11:513–556.https://doi.org/10.1016/0022-1910(65)90137-X
-
An endogenous elevation of cGMP increases the excitability of identified insect neurosecretory cellsJournal of Comparative Physiology A: Sensory, Neural, and Behavioral Physiology 180:329–337.https://doi.org/10.1007/s003590050052
-
Neuropeptide Hierarchies and the Activation of Sequential Motor Behaviors in the Hawkmoth, Manduca sextaThe Journal of Neuroscience 17:4389–4397.https://doi.org/10.1523/JNEUROSCI.17-11-04389.1997
-
Ecdysteroid metabolism in leechesInvertebrate Reproduction & Development 15:57–68.https://doi.org/10.1080/07924259.1989.9672022
-
The new view of animal phylogenyAnnual Review of Ecology, Evolution, and Systematics 35:229–256.https://doi.org/10.1146/annurev.ecolsys.35.112202.130124
-
MAFFT multiple sequence alignment software version 7: improvements in performance and usabilityMolecular Biology and Evolution 30:772–780.https://doi.org/10.1093/molbev/mst010
-
Precursors of neuropeptides and peptide hormones in the genomes of tardigradesGeneral and Comparative Endocrinology 267:116–127.https://doi.org/10.1016/j.ygcen.2018.06.012
-
The occurrence of ecdysteroids in the cestode, Moniezia expansaMolecular and Biochemical Parasitology 10:123–138.https://doi.org/10.1016/0166-6851(84)90001-X
-
Ecdysteroids in the dorsal bodies of pulmonates (gastropoda): Synthesis and release of ecdysoneComparative Biochemistry and Physiology Part A: Physiology 84:777–782.https://doi.org/10.1016/0300-9629(86)90405-6
-
Enzymology of the nematode cuticle: a potential drug target?International Journal for Parasitology: Drugs and Drug Resistance 4:133–141.https://doi.org/10.1016/j.ijpddr.2014.05.003
-
Deletion of the ecdysis-triggering hormone gene leads to lethal ecdysis deficiencyDevelopment 129:493–503.
-
Two subtypes of ecdysis triggering hormone receptor in Drosophila MelanogasterThe Journal of Biological Chemistry.
-
Ecdysis of decapod crustaceans is associated with a dramatic release of crustacean cardioactive peptide into the haemolymphThe Journal of Experimental Biology 203:521–536.
-
Developmental checkpoints and feedback circuits time insect maturationCurrent Topics in Developmental Biology 103:1–33.https://doi.org/10.1016/B978-0-12-385979-2.00001-0
-
Regulation of the C. elegans molt by pqn-47Developmental Biology 360:297–309.https://doi.org/10.1016/j.ydbio.2011.09.025
-
Halloween genes in Panarthropods and the evolution of the early moulting pathway in ecdysozoaRoyal Society Open Science 5:180888.https://doi.org/10.1098/rsos.180888
-
CCAP and FMRFamide-like peptides accelerate the contraction rate of the antennal accessory pulsatile organs (auxiliary hearts) of mosquitoesThe Journal of Experimental Biology 219:2388–2395.https://doi.org/10.1242/jeb.141655
-
The evolution of the ecdysozoaPhilosophical Transactions of the Royal Society B: Biological Sciences 363:1529–1537.https://doi.org/10.1098/rstb.2007.2243
-
Xenacoelomorph neuropeptidomes reveal a major expansion of neuropeptide systems during early bilaterian evolutionMolecular Biology and Evolution 35:2528–2543.https://doi.org/10.1093/molbev/msy160
-
Hormonal control of insect ecdysis: endocrine cascades for coordinating behavior with physiologyVitamins & Hormones 73:1–30.
-
Neurohormones and neuropeptides encoded by the genome of Lottia gigantea, with reference to other mollusks and insectsGeneral and Comparative Endocrinology 167:86–103.https://doi.org/10.1016/j.ygcen.2010.02.010
-
In silico cloning of genes encoding neuropeptides, neurohormones and their putative G-protein coupled receptors in a spider miteInsect Biochemistry and Molecular Biology 42:277–295.https://doi.org/10.1016/j.ibmb.2011.12.009
-
Molecular characterization of ecdysis triggering hormone and its receptor in Citrus red mite (Panonychus citri)Comparative Biochemistry and Physiology Part A: Molecular & Integrative Physiology 230:100–105.https://doi.org/10.1016/j.cbpa.2019.01.003
-
Complex steroid-peptide-receptor cascade controls insect ecdysisGeneral and Comparative Endocrinology 153:88–96.https://doi.org/10.1016/j.ygcen.2007.04.002
-
Insect Endocrinology253–309, Neuroendocrine regulation of ecdysis, Insect Endocrinology, Academic Press.
Article and author information
Author details
Funding
Austrian Science Fund (P29455-B29)
- Andreas Wanninger
Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (6090/13-3)
- André Luiz de Oliveira
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Copyright
© 2019, de Oliveira et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 5,185
- views
-
- 630
- downloads
-
- 37
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Evolutionary Biology
- Genetics and Genomics
Evolutionary arms races can arise at the contact surfaces between host and viral proteins, producing dynamic spaces in which genetic variants are continually pursued. However, the sampling of genetic variation must be balanced with the need to maintain protein function. A striking case is given by protein kinase R (PKR), a member of the mammalian innate immune system. PKR detects viral replication within the host cell and halts protein synthesis to prevent viral replication by phosphorylating eIF2α, a component of the translation initiation machinery. PKR is targeted by many viral antagonists, including poxvirus pseudosubstrate antagonists that mimic the natural substrate, eIF2α, and inhibit PKR activity. Remarkably, PKR has several rapidly evolving residues at this interface, suggesting it is engaging in an evolutionary arms race, despite the surface’s critical role in phosphorylating eIF2α. To systematically explore the evolutionary opportunities available at this dynamic interface, we generated and characterized a library of 426 SNP-accessible nonsynonymous variants of human PKR for their ability to escape inhibition by the model pseudosubstrate inhibitor K3, encoded by the vaccinia virus gene K3L. We identified key sites in the PKR kinase domain that harbor K3-resistant variants, as well as critical sites where variation leads to loss of function. We find K3-resistant variants are readily available throughout the interface and are enriched at sites under positive selection. Moreover, variants beneficial against K3 were also beneficial against an enhanced variant of K3, indicating resilience to viral adaptation. Overall, we find that the eIF2α-binding surface of PKR is highly malleable, potentiating its evolutionary ability to combat viral inhibition.
-
- Evolutionary Biology
- Genetics and Genomics
It is well established that several Homo sapiens populations experienced admixture with extinct human species during their evolutionary history. Sometimes, such a gene flow could have played a role in modulating their capability to cope with a variety of selective pressures, thus resulting in archaic adaptive introgression events. A paradigmatic example of this evolutionary mechanism is offered by the EPAS1 gene, whose most frequent haplotype in Himalayan highlanders was proved to reduce their susceptibility to chronic mountain sickness and to be introduced in the gene pool of their ancestors by admixture with Denisovans. In this study, we aimed at further expanding the investigation of the impact of archaic introgression on more complex adaptive responses to hypobaric hypoxia evolved by populations of Tibetan/Sherpa ancestry, which have been plausibly mediated by soft selective sweeps and/or polygenic adaptations rather than by hard selective sweeps. For this purpose, we used a combination of composite-likelihood and gene network-based methods to detect adaptive loci in introgressed chromosomal segments from Tibetan WGS data and to shortlist those presenting Denisovan-like derived alleles that participate to the same functional pathways and are absent in populations of African ancestry, which are supposed to do not have experienced Denisovan admixture. According to this approach, we identified multiple genes putatively involved in archaic introgression events and that, especially as regards TBC1D1, RASGRF2, PRKAG2, and KRAS, have plausibly contributed to shape the adaptive modulation of angiogenesis and of certain cardiovascular traits in high-altitude Himalayan peoples. These findings provided unprecedented evidence about the complexity of the adaptive phenotype evolved by these human groups to cope with challenges imposed by hypobaric hypoxia, offering new insights into the tangled interplay of genetic determinants that mediates the physiological adjustments crucial for human adaptation to the high-altitude environment.