Abstract
The Subolesin/Akirin constitutes a good model for the study of functional evolution because these proteins have been conserved throughout the metazoan and play a role in the regulation of different biological processes. Here, we investigated the evolutionary history of Subolesin/Akirin with recent results on their structure, protein-protein interactions and function in different species to provide insights into the functional evolution of these regulatory proteins, and their potential as vaccine antigens for the control of ectoparasite infestations and pathogen infection. The results suggest that Subolesin/Akirin evolved conserving not only its sequence and structure, but also its function and role in cell interactome and regulome in response to pathogen infection and other biological processes. This functional conservation provides a platform for further characterization of the function of these regulatory proteins, and how their evolution can meet species-specific demands. Furthermore, the conserved functional evolution of Subolesin/Akirin correlates with the protective capacity shown by these proteins in vaccine formulations for the control of different arthropod and pathogen species. These results encourage further research to characterize the structure and function of these proteins, and to develop new vaccine formulations by combining Subolesin/Akirin with interacting proteins for the control of multiple ectoparasite infestations and pathogen infection.
Introduction
Akirin, from the Japanese “akiraka ni suru" meaning “making things clear," was first identified by Goto et al. (2008) as a key component of the immune deficiency (IMD) and Tumor necrosis factor (TNF)/Toll-like receptor (TLR)-nuclear factor-kappa B (NF-kB) (TNF/TLR) signaling pathways in Drosophila melanogaster and Mus musculus, respectively. However, previous reports identified akirin as a gene involved in developmental processes in flies (Peña-Rangel et al., 2002; DasGupta et al., 2005). Subolesin, from the Latin “suboles" meaning “progeny," was first reported in 2003 with its discovery as the candidate protective antigen 4D8 by expression library immunization in Ixodes scapularis (Almazán et al., 2003). Gene orthology is a key concept in functional evolution (Koonin, 2005). Orthologs genes, defined as derived from a single ancestral gene that diverged during speciation, usually perform equivalent or identical functions, while paralogs that originated after gene duplication are considered to have more divergent functions (Koonin, 2005; Adipietro et al., 2012). Studies at genome level have identified many orthologs genes between divergent species, but the functional equivalency of the proteins encoded by these genes has not been fully characterized (Koonin, 2005).
Subolesin/Akirin are encoded by orthologs evolutionarily conserved throughout the metazoan that play a role in the regulation of different biological processes including immune response (de la Fuente et al., 2006a; Beutler and Moresco, 2008; Goto et al., 2008; Galindo et al., 2009; Macqueen, 2009; Macqueen and Johnston, 2009; Wan and Lenardo, 2010). Two akirin paralogs encoding Akirin1 and Akirin2 have been identified in vertebrates, but the functional homolog for invertebrate Subolesin/Akirin appears to be Akirin2 ortholog (Beutler and Moresco, 2008; Goto et al., 2008; Galindo et al., 2009; Macqueen and Johnston, 2009).
Understanding the function of the cell interactome (protein-protein physical and functional interactions) and regulome (transcription factors-target genes interactions) in response to infection is critical toward a better understanding of host-pathogen interactions and the identification of potential targets for new interventions for the prevention and control of tick infestations and tick-borne diseases (Rioualen et al., 2017; de la Fuente, 2018). Subolesin/Akirin are involved in both cell interactome and regulome, and constitute a good model for the study of the functional evolution of these processes in response to infection. In this review, we integrated the evolutionary history of Subolesin/Akirin with recent results on their structure and function in different species to provide insights into the functional evolution of these regulatory proteins, and their potential as vaccine antigens for the control of ectoparasite infestations and pathogen infection.
Evolution of Subolesin/Akirin
The phylogenetic analysis of subolesin/akirin coding sequences using an updated sequence database (Figure 1A and Supplementary Figure S1) expanded the information on the evolution of these genes, and supported the results reported previously by Macqueen and Johnston (2009) that akirin1 and akirin2 are vertebrate-specific paralogs that form a separate clade from invertebrate subolesin/akirin. In some vertebrate species, Akirins constitute a family of paralog proteins that probably originated as a result of whole-genome duplications (Macqueen and Johnston, 2009; Macqueen et al., 2010a,b; Liu et al., 2015). After akirin duplication, akirin1 evolved faster than akirin2, the ortholog of tick subolesin (Figure 1A; Macqueen and Johnston, 2009). Furthermore, the loss of some akirin paralogs may have also occurred after genome duplications (Macqueen et al., 2010b; Liu et al., 2015). For example, the subolesin/akirin gene family consists of a single member in invertebrates (subolesin/akirin), birds and reptiles (subolesin/akirin2), two members in amphibians and mammals (akirin1 and akirin2), two to three members in teleosts, and more than three members in Salmonidae (Macqueen et al., 2010a,b; Liu et al., 2015).
FIGURE 1
The phylogenetic analysis of subolesin gene sequences was performed in 42 species belonging to 6 and 1 genera of hard (family Ixodidae) and soft (family Argasidae) ticks, respectively. The analysis corroborated previous results showing a reductive evolution in protein length (de la Fuente et al., 2006a; Figure 1B). The Subolesin amino acid (aa) sequence evolved from 173 to 184 aa in Ornithodoros, Ixodes and Haemaphysalis spp. to 161 aa in Rhipicephalus spp. (Figure 1B). It is generally accepted that evolution proceeds toward greater complexity at both the organismal and genomic levels. However, numerous examples of reductive evolution of parasites and symbionts have been described to challenge this notion (Wolf and Koonin, 2013). Wolf and Koonin (2013) proposed that quantitatively, the evolution of genomes appears to be dominated by reduction and simplification, punctuated by episodes of complexification. The reductive evolution process has been particularly documented in genomes that replicate within the domain of a host genome (Andersson and Kurland, 1998; Driscoll et al., 2017), but it has also been proposed to be involved in the origin of bacteria from eukaryotes (Staley, 2017). In arthropods, reductive evolution has been implicated in the evolutionary origin of other proteins such as type IV classical cadherins (Sasaki et al., 2017).
The protein length is subjected to systematic variation that relates to the cellular context in which it functions (Wang et al., 2011). For growth rate-optimized cells, the reduction in protein length constitutes an advantage by increasing their mass-normalized kinetic efficiencies (Ehrenberg and Kurland, 1984; Kurland et al., 2007; Wang et al., 2011). Consequently, shorter proteins that retain maximum functional rates are expected to support faster cell growth rates than longer proteins with similar kinetic characteristics. Wang et al. (2011) proposed the use of the terms “domain” and “linker” to refer to protein folded domains and nondomain regions, respectively. Proteins with nondomain sequences are proteins intrinsically unstructured or natively unfolded that lack a stable tertiary structure but have a dynamic range of conformations (Orengo and Thornton, 2005). These proteins appear to be more abundant in eukaryotes when compared to prokaryotes and are usually involved in binding and molecular recognition (Brown et al., 2011). Subolesin/Akirin were characterized as a linker with three predicted disordered nondomain regions that resulted in unstructured proteins (Prudencio et al., 2010; see also below). These results agreed with the findings of Wang et al. (2011) that the evolutionary reductive constraints on protein lengths are preferentially expressed in linker sequences.
It is difficult to establish a comprehensive record of ticks developmental rates because most of the experiments in previously published papers have been done at different regimes of temperature, relative humidity and photoperiod, all factors affecting the time in which ticks complete each developmental stage. However, data compiled by Hoogstraal (1956), Morel (2003) and Horak et al. (2018) under similar conditions established that in the range of 24–28°C, ticks of the genera Hyalomma and Rhipicephalus complete their life cycle in about 33% less time than ticks of the genera Ixodes or Amblyomma. The Hyalomma and Rhipicephalus spp. are considered the two most recent genera of ticks, while Ixodes spp. and Amblyomma are among the most ancient splits of tick lineages (Mans et al., 2016; Figure 1B). Therefore, it is possible that the reductive evolution of Subolesin is associated with faster developmental rates in Rhipicephalus and Hyalomma spp. when compared to more ancient tick species even if they are sympatric. The faster developmental rate in recently evolved tick species may be associated with increasing cell growth rates that have been associated with reductive evolution (Wang et al., 2011). However, the complete association between existing data about developmental rates and evolutionary features of ticks requires further research.
At the genome level, subolesin/akirin exon-intron architecture shows a clear evolutionary pattern (Figure 2). As shown for subolesin/akirin coding sequences (Figure 1A), the vertebrate-specific paralogs form a separate clade from invertebrate genes (Figure 2). The subolesin/akirin genes evolved from 4 exons in invertebrates to 5 exons in vertebrates. The exon sizes showed a pattern of larger to shorter for vertebrate akirin1 exons I to V and mosquito akirin exons I to IV (Figure 2). However, for tick subolesin and vertebrate akirin2 genes, the size of exons II and III, respectively, were larger or similar than that of preceding exon (Figure 2). Furthermore, while the length of 11 nucleotides (nt) in the last exon was conserved across evolution, the length of the penultimate exon evolved from 157 nt in mosquitoes to 72 nt in ticks and vertebrates (Figure 2). However, the length of the exon-intron sequence increased from 3,992 nt in mosquitoes to 22,829 nt in ticks and then decreased to 8,862 nt in fish to increase back again to 26,088 nt in humans (Figure 2). These results do not correlate with genome sizes of these organisms (Gregory, 2005), and may originated from still unknown evolutionary events.
FIGURE 2
Subolesin/Akirin Structure and Its Interactions With Dna and Transcription Factors
Akirins are involved in cellular processes that are regulated by specific domains and binding sites. The rat Akirin2 (or 14-3-3β Interactant 1) described by Komiya et al. (2008) is 48% identical to the I. scapularis Subolesin and was previously used to identify the conserved Akirin nuclear localization signal (NLS) domains and its binding sites (Macqueen and Johnston, 2009). The pairwise sequence alignment shows that both NLS domains of the rat Akirin2 are conserved in Subolesin, but with variations in a few binding sites (Figure 3A). Subolesin binding sites 1 and 4 are similar to Akirin2 with Subolesin binding sites 2 and 3 each possessing a single substitution. The Subolesin binding site 5, however, has a double Leu insertion when compared to Akirin2 (Figure 3A). Currently there are no resolved Subolesin/Akirin structures and standard sequence-based bioinformatics methods lack parameters for locating structural homologs in the Protein Databank (PDB). Therefore, several logical steps were taken to correctly model the I. scapularis Subolesin. The Subolesin sequence (Figure 3A) was initially submitted to I-TASSER (Zhang, 2008), a protein multiple threading algorithm that is considered a top competitor in the Critical Assessment of Structure Prediction1. The I-TASSER algorithm resulted in five distinct Subolesin models that were then individually submitted to the DALI server (Holm and Laakso, 2016) for identifying PDB structural homologs of similar length with minimalα-carbon backbone deviations between the two global structures. Since Subolesin/Akirin are effectors of the IMD/TNF/TLR Relish/NF-kB signaling pathways (Goto et al., 2008; Naranjo et al., 2013), a match with homologous structures that potentially participate within these pathways was also a criterion in selecting an appropriate Subolesin model. This logical approach resulted in a Subolesin model homologous in structure, with only 7% residue conservation, to the genetically engineered catalytic domain of the transposase Sleeping Beauty (Ivics et al., 1997; Zanesi et al., 2013; Figure 3B).
FIGURE 3
As part of the Tc1/mariner transposon superfamily, the Sleeping Beauty transposase has been engineered for genetic screening studies, leading to the discovery of several genes activated by Sleeping Beauty transposon insertions that participate in the NF-kB signaling pathway (Zanesi et al., 2013). The Sleeping Beauty transposase sequence is composed of an N-terminus paired-like domain with a leucine zipper (∼90 residues long) and the C-terminus folds as the catalytic domain (Ivics et al., 1997). The catalytic domain of Sleeping Beauty was resolved with a DNA transposon end and modeled with a target DNA revealing the mechanism of hyperactive Sleeping Beauty mutation screening studies while discovering novel variants for future screenings (Voigt et al., 2016). The Sleeping Beauty crystal structure details that its catalytic domain has a global homology to Ribonuclease H (RNase H) (Voigt et al., 2016). The RNase H-like protein fold forms a catalytic triad (Asp-Asp-Glu) that coordinates metal ions involved in excision and insertions of DNA (Voigt et al., 2016). By resolving the Sleeping Beauty catalytic domain, Voigt et al. (2016) also discovered that the Gly-rich box (located on the clamp loop) is involved in protein-protein interactions, specifically with partnering monomers in the DNA complex. The conserved positions of the catalytic triad and the Gly-rich box, however, are not present in Subolesin/Akirin sequences.
Prior to acting in the Tc1/mariner transposon system, Sleeping Beauty must enter the nucleus. Passage to the nucleus is controlled by NLS domains that have a strong affinity to karyopherin/importin receptors, proteins responsible for transporting NLS-tagged “cargo” in and out of the nucleus via nuclear pores (Leung et al., 2003). The N-terminus of the Sleeping Beauty catalytic domain contains a NLS domain that is quite long (17 residues) and is actually a bipartite NLS (Ivics et al., 1997). A monopartite NLS domain has the formulation Lys-Lys/Arg-[X]-Lys-Lys/Arg, where [X] is any other (∼2) amino acids. A bipartite NLS domain has a linker sequence, where [X] is ∼10 residues long (Makkerh et al., 1996). Mutations of these upstream NLS residues prior to the linker sequence [X] has shown to reduce protein entry into the nucleus (Dingwall et al., 1988), and inhibit NLS binding to karyopherin/importin receptors (Leung et al., 2003). These upstream, positively charged residue pairs are in the aligned sequences of Figure 3A (orange-labeled residues), indicating that Subolesin/Akirin2 NLS1 is bipartite. The two positively charged residues highlighted near NLS2 (Figure 3A) indicate that the NLS2 of Subolesin/Akirin is actually a longer monopartite domain. The modeled Subolesin structure has its NLS1 positioned on the clamp loop, which differs from Sleeping Beauty that is coordinated at the N-terminus α-helix (I) (Figure 3B). The Subolesin NLS2 domain, not present in the catalytic domain of Sleeping Beauty, is located on an α-helix (II) outside the central core of the protein. The absence of the catalytic triad of Sleeping Beauty (Voigt et al., 2016) and metal binding sites in Subolesin/Akirin support that these proteins do not act as a transposase. Additionally, the Subolesin/Akirin binding sites have long been recognized by mutation studies (Komiya et al., 2008), and Subolesin/Akirin RNA interference (RNAi) experiments have shown to disrupt the Relish signaling pathway (Goto et al., 2008; Naranjo et al., 2013). However, as discussed bellow, the Subolesin/Akirin interactome has not been fully characterized, and whether Subolesin/Akirin binding partners are only proteins or also include nucleic acids.
The superposed structures of Subolesin and Sleeping Beauty depict a global homology with a low α-carbon backbone deviation of 0.3 nm (Figure 3B). There are, however, missing and disordered secondary structures. The Sleeping Beauty catalytic domain has five β-sheets surrounded by five α-helices. As previously mentioned, the N-terminus α-helix I of Sleeping Beauty that contains its NLS domain is missing in the Subolesin model (encircled in Figure 3B), thereby shifting the Subolesin NLS1 domain to the clamp loop. Moreover, the β-sheets of Subolesin are highly disordered. Future experiments should resolve the stacking and conformations of the disordered Subolesin β-sheets by X-ray crystallography. Nevertheless, the remaining four α-helices (II-V) are structurally conserved, and the Subolesin/Akirin binding sites 2-5 are positioned on or approximating α-helix III (Figure 3B). As in the primary sequence (Figure 3A), the Subolesin/Akirin binding site 1 is structurally distant from the other sites, located on the N-terminus clamp loop near the position of NLS1 (Figure 3B). As a transposase, the clamp loop of the resolved Sleeping Beauty catalytic domain is not in its DNA-bound conformation. Therefore, Voigt et al. (2016) modeled the clamp loop after the DNA-bound transposase, Mos1, from Drosophila mauritiana (Richardson et al., 2009). The Mos1 also has poor sequence identity to Subolesin (<5%), but are structurally homologous with α-carbon backbone deviation of 0.34 nm. This led to a subsequent I-TASSER simulation using the template Mos1 as conducted by Voigt et al. (2016). The resulting model has Subolesin bound to a DNA duplex with an adequate clamp loop conformation that extends downstream the duplex (Figure 3C).
The Subolesin-DNA complex show several residues on α-helices IV and V that approximate the DNA 5′-end of the sense strand and the 3′-end of the antisense strand (Figure 3C). Four of the seven residues positioned on α-helix IV mainly interact with the phosphate backbone of the antisense strand. The remaining three residues that interact with both strands are after the pivotal kink of α-helix V (encircled in Figure 3C). The alignment in Figure 3A highlights these novel DNA-binding residues and indicates the α-helix on which they are positioned. The binding sites 2–5 and NLS2 are distal to the DNA interacting site, while the clamp loop containing NLS1 and binding site one is downstream the DNA duplex. Given the structural coordination of Subolesin bound to DNA (Figure 3C), binding sites 2–5 and NLS2 may interact with nucleotides upstream the DNA or with additional co-transcription factors (CTF) (Figure 3D). The transcription factors (TF) Relish/NF-kB form DNA-protein complexes with CTFs, and Subolesin is hypothesized to act as a CTF of Relish via an intermediate CTF (Figure 3D and described in the next section). Furthermore, the extended DNA downstream position of the Subolesin clamp loop with binding site one will hypothetically coordinate the CTF and possibly the TF (Figure 3D). If Subolesin is a CTF, how will it conform while the mRNA is being transcribed? Positively charged residues, specifically Lys, recognize RNA strands via electrostatic interactions (Law et al., 2006). Given conformational flexibility of the Subolesin clamp loop and the fact that it does not contain the Sleeping Beauty Gly-rich box, the proximity of positively charged Subolesin NLS1 domain residues downstream the DNA duplex may guide transcribing mRNAs for post-transcriptional processing.
Function of Subolesin/Akirin
The innate immune response acts as the first line of defense against pathogen infection in all metazoans, and constitutes the only immune response in invertebrates (Medzhitov et al., 1997; Shaw et al., 2017). It has been shown that invertebrate Subolesin/Akirin and vertebrate Akirin2 act in concert with Relish/NF-kB to induce the expression of a subset of downstream pathway elements in the IMD and TNF/TLR signaling pathways involved in the immune response to pathogen infection. This function of Subolesin/Akirin has been documented in ticks (de la Fuente et al., 2008, 2017a; Naranjo et al., 2013; Gulia-Nuss et al., 2016; Shaw et al., 2017), fruit fly (Beutler and Moresco, 2008; Goto et al., 2008; Wan and Lenardo, 2010; Bonnay et al., 2014), human (Goto et al., 2008), salmon (Macqueen et al., 2010b), shrimp (Hou et al., 2013), Japanese flounder (Yang et al., 2013), amphioxus (Yan et al., 2013), rock bream (Kasthuri et al., 2013), Chinese loach (Xue et al., 2014), Hong Kong oyster (Qu et al., 2014), turbot (Yang et al., 2011), sea louse (Valenzuela-Muñoz and Gallardo-Escárate, 2014), mouse (Beutler and Moresco, 2008; Goto et al., 2008; Wan and Lenardo, 2010; Tartey et al., 2015), croaker (Liu et al., 2015); shrimp (Liu Y. et al., 2016; Liu et al., 2016), and seahorse (Pavithiran et al., 2018; Figures 4, 5).
FIGURE 4
FIGURE 5
Subolesin/Akirin are also involved in the Relish/NF-kB independent gene regulation (Figure 4), thus playing a role in various biological processes in addition to the immune response (Figure 5). These processes include animal reproduction and development, causing lethal embryonic or reduced growth phenotypes in knockout mice, fruit fly, ticks, and nematodes (Maeda et al., 2001; de la Fuente et al., 2006a; Goto et al., 2008; Carpio et al., 2013; Qu et al., 2014), metazoan myogenesis (Marshall et al., 2008; Salerno et al., 2009; Macqueen et al., 2010a; Mobley et al., 2014; Sun et al., 2016), Xenopus neural development (Liu et al., 2017), meiosis/carcinogenesis (Komiya et al., 2008; Macqueen et al., 2010b; Clemons et al., 2013), tick stress response, feeding, growth and reproduction (Almazán et al., 2003, 2005; de la Fuente et al., 2006a, 2008; Smith et al., 2009; Busby et al., 2012; Rahman et al., 2018), pathogen infection and transmission in ticks (de la Fuente et al., 2006b, 2016, 2017a; Zivkovic et al., 2010a,b; Busby et al., 2012; Hajdušek et al., 2013) and turbot (Yang et al., 2011), human glioblastoma cell apoptosis (Krossa et al., 2015), cattle marbling (Sasaki et al., 2009; Watanabe et al., 2011; Kim et al., 2013), and mouse mitogen-induced B cell cycle progression and humoral immune responses (Tartey et al., 2015). For example, as previously reported in I. scapularis and other tick species (de la Fuente and Kocan, 2006; de la Fuente et al., 2006a,b, 2011, 2013; Merino et al., 2013a; de la Fuente and Contreras, 2015), Subolesin appears to function in multiple biological processes such as tick response to infection, feeding, reproduction, development and stress response (Figure 6).
FIGURE 6
Akirin1 and Akirin2 have also different functions in vertebrates, which is illustrated by the role of Akirin1 in myogenesis while Akirin2 promotes meiosis/carcinogenesis (Macqueen and Johnston, 2009; Macqueen et al., 2010a,b; Figure 4). These different functions may be related to the Akirin subcellular localization. While Akirin1 is found in the nucleus, Subolesin/Akirin2 is located in both cytoplasm and nucleus (de la Fuente et al., 2011; Antunes et al., 2014; Krossa et al., 2015; Pavithiran et al., 2018; Figure 7). The subcellular localization of Subolesin/Akirin2 is probably related to its structure, which as discussed above contains NLS domains that are involved in protein transport in and out of the nucleus via nuclear pores (Leung et al., 2003).
FIGURE 7
In summary and based on current information, Subolesin/Akirin evolved with similar functions in both invertebrates and vertebrates (Figure 5). The annotation of some biological processes described in certain taxa only may be due to the presence of species-specific functions or more likely a consequence of the still incomplete characterization of Subolesin/Akirin function in the different species.
Subolesin/Akirin Role in Cell Interactome and Regulome
Subolesin/Akirin are proteins without catalytic or DNA-binding capacity. How Subolesin/Akirin regulate gene expression is unknown but likely involve interactions with proteins with DNA-binding or chromatin-remodeling capacity (de la Fuente et al., 2008; Beutler and Moresco, 2008; Nowak and Baylies, 2012; Nowak et al., 2012; Naranjo et al., 2013). In this way, Subolesin/Akirin link the activities of transcription factors with those of chromatin remodeling complexes to influence gene expression in a context-dependent manner (Nowak and Baylies, 2012; Nowak et al., 2012). For example, Subolesin/Akirin and Akirin2 as key components of the innate immune response can directly or indirectly interact with other regulatory proteins such as “14-3-3” proteins, DNA methyltransferase-associated protein 1 (DMAP1) and the basic helix–loop–helix transcription factor (Twist) to up- or down-regulate transcription (Gonzalez and Baylies, 2005; Goto et al., 2008, 2014; Komiya et al., 2008; de la Fuente et al., 2008; Beutler and Moresco, 2008; Nowak et al., 2012; Naranjo et al., 2013; Tartey et al., 2014; Gulia-Nuss et al., 2016; Shaw et al., 2017).
The Subolesin/Akirin role in the cell interactome and regulome in response to different stimuli has not been characterized. Recently, we proposed a method based on the graph theory for the analysis of human and tick cell proteome in response to A. phagocytophilum infection (Estrada-Peña et al., 2018). This approach resulted in a network of interacting proteins and cell processes clustered in biological pathways, and ranked with indexes representing the topology of the proteome influenced by features of the interactome and regulome. The results evidenced differences in the response to A. phagocytophilum infection between human and tick cells, and supported that human neutrophils but not tick cells limit pathogen infection through differential representation of ras-related proteins (Estrada-Peña et al., 2018). Herein, this method was applied to predict the position of Subolesin in the regulome of tick cells and in response to A. phagocytophilum infection (Figure 8). The results showed that Subolesin is deeply involved in the core of transcription processes but also in other secondary processes such as transcription from RNA polymerase II promoter, DNA repair, and chromatin remodeling (Figure 8). Furthermore, other processes that change in infected cells when compared to uninfected cells (i.e., signal transduction, regulation of transcription, and response to heat) are deeply linked to the central transcription process. The putative Subolesin role in these processes varied between infected and uninfected cells (the width of the lines is proportional to the importance of the link between proteins and processes; Figure 8). For example, it appears that particularly in infected cells other proteins but Subolesin have a more prominent role in the strong protein link with transcription and transcription from RNA polymerase II promoter processes (Figure 8). These results predict the role that Subolesin plays in the regulation of different biological processes, and its differential role in response to A. phagocytophilum infection in tick cells. However, the lack of a more prominent role for Subolesin may respond to the fact that this protein does not bind directly to DNA but interacts with other proteins to exert its regulatory function.
FIGURE 8
In an attempt to provide information on the Subolesin/Akirin interactome, the information on Subolesin/Akirin-protein physical and functional interactions was compiled from the String protein-protein interactions database2 (Figure 9 and Supplementary Dataset S1). Based on the analysis of protein-protein interactions, the results did not allow establishing an evolutionary signature of the Subolesin/Akirin2 interactome (Figure 9), probably due to the limited information available. Nevertheless, similar Subolesin/Akirin2 interacting proteins were described in fly and mouse (B7PRT9, Brahma/SWI2-related protein BRG-1) and in fish and rat (B7P8Y4, Arginyl-tRNA synthetase) suggesting possible evolutionarily conserved protein-protein interactions (Figure 9). To further gain insight into the evolution of the Subolesin/Akirin2 interactome, instead of looking only at protein-protein interactions, the interacting proteins were annotated according to the biological processes in which they are involved (Figure 9). The results of this analysis showed that the biological processes affected by the Subolesin/Akirin interactome are evolutionarily conserved, with metabolic process (MP), cellular process (CP) and biological regulation (BR) being among the most represented processes in all organisms (Figure 9).
FIGURE 9
Protective Capacity of Subolesin/Akirin for the Control of Ectoparasite Vector Infestations and Pathogen Infection
Subolesin was discovered and characterized as a tick protective antigen for the control of I. scapularis infestations (Almazán et al., 2003, 2005; Sonenshine et al., 2006). Since then, Subolesin/Akirin showed a protective capacity in vaccines for the control of infestations by different arthropod species and pathogen infection and transmission (reviewed by de la Fuente et al., 2006a, 2011, 2013; de la Fuente and Kocan, 2006, 2014; Merino et al., 2013a,b; de la Fuente and Contreras, 2015). The putative mechanism for Subolesin vaccine protection was described by de la Fuente et al. (2011). They showed that by still unknown mechanisms anti-Subolesin antibodies could enter into tick cells where they can interact with cytosolic Subolesin to prevent its translocation to the nucleus and therefore the possibility to exert it regulatory functions.
The development of vaccines for the control of multiple arthropod ectoparasites constitutes a priority for targeting the various species infesting the same host (de la Fuente et al., 2011; Contreras and de la Fuente, 2016b). In this context, Subolesin/Akirin appears as a promising vaccine protective antigen due to its conservation in sequence and function through evolution (Figures 1, 5, 9). In fact, recent results support the potential of Subolesin/Akirin as a vaccine protective antigen for the control of multiple ectoparasite vector species and transmitted pathogens (Figure 10). The arthropod ectoparasite species in which Subolesin/Akirin have shown protective capacity by affecting different phases of their life cycles include the genera Aedes, Anopheles, Phlebotomus, Caligus, Dermanyssus, Ornithodoros, Ixodes, Haemaphysalis, Amblyomma, Dermacentor, Hyalomma and Rhipicephalus (Almazán et al., 2003, 2005, 2010, 2012; de la Fuente et al., 2006a, 2010, 2011, 2013; Canales et al., 2009; Harrington et al., 2009; Zivkovic et al., 2010a; Moreno-Cid et al., 2011, 2013; Carpio et al., 2011; Bensaci et al., 2012; Merino et al., 2011a,b, 2013b; Carreón et al., 2012; Manzano-Román et al., 2012, 2015; de la Fuente and Kocan, 2014; da Costa et al., 2014; Shakya et al., 2014; Torina et al., 2014; Contreras et al., 2015; de la Fuente and Contreras, 2015; Lu et al., 2016; Olds et al., 2016; Contreras and de la Fuente, 2016a; Kumar et al., 2017; Villar et al., 2017; Figure 10). The reduction in pathogen infection has been shown for tick-borne pathogens, Anaplasma marginale, A. phagocytophilum, Borrelia burgdorferi s.s. and Babesia bovis, and the mosquito-borne pathogen, Plasmodium berghei (de la Fuente et al., 2011; Merino et al., 2011b, 2013a; Bensaci et al., 2012; de la Fuente and Kocan, 2014; da Costa et al., 2014; Figure 10). Furthermore, protective epitopes were mapped in Subolesin/Akirin (Prudencio et al., 2010) and chimeric antigens were designed showing protective capacity in vaccinated hosts against tick infestations (Merino et al., 2013b; Contreras and de la Fuente, 2016a).
FIGURE 10
Considering the protective capacity shown by Subolesin/Akirin on different vector and pathogen species, future research directions will include the characterization of its protective capacity in other arthropod ectoparasite species, and the combination with other vector-derived and pathogen-derived antigens to increase vaccine efficacy for the control of both vector infestations and pathogen infection (Schetters et al., 2016; de la Fuente et al., 2017b; de la Fuente, 2018).
Conclusion and Future Directions
Significant advances have been made recently toward understanding the evolution and function of Subolesin/Akirin. Our results suggest that Subolesin/Akirin evolved conserving not only its sequence and structure, but also its function and role in cell interactome and regulome in response to pathogen infection and other biological processes. However, major challenges remain in fully characterizing the interactome and function of these proteins, their role in the cell regulome in response to different stimuli, and how their evolution can meet species-specific demands. Furthermore, the structure of Subolesin/Akirin and interacting molecules should be resolved by X-ray crystallography to better understand their function. Finally, the conserved functional evolution of Subolesin/Akirin correlates with the protective capacity shown by these proteins in vaccine formulations for the control of different arthropod and pathogen species, and encourage further research to develop new vaccine formulations by combining Subolesin/Akirin with interacting proteins for the control of multiple ectoparasite infestations and pathogen infection.
Statements
Author contributions
JdlF conceived the paper. SA-J, MV, AC-C, JV, and AE-P performed the data analyses. PA and SA-J performed the microscopy studies. JF, SA-J, JV, and AC-C wrote the manuscript. All authors approved and contributed to the final version of the manuscript.
Funding
This work was financially supported by the Ministerio de Economía, Industria y Competitividad (Spain) grant BFU2016-79892-P. MV was funded by the Universidad de Castilla La Mancha, Spain. JV was supported by Project FIT (Pharmacology, Immunotherapy, nanoToxicology), funded by the European Regional Development Fund. The publication fee was partially supported by the CSIC Open Access Publication Support Initiative through its Unit of Information Resources for Research (URICI).
Acknowledgments
We acknowledge a grant for the development of a research organization RVO: RO0516 from the Veterinary Research Institute, Brno, Czechia.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2018.01612/full#supplementary-material
FIGURE S1Phylogenetic tree of akirin and subolesin nucleotide sequences. The figure displays a Neighbor Joining (NJ) phylogenetic tree of 361 nucleotide sequences belonging to 152 families, 73 orders and 15 classes. GenBank accession numbers and species names are shown. Tree reconstruction method was as described in Figure 1.
DATASET S1Annotation of the Subolesin/Akirin interactome. The information on Subolesin/Akirin-protein interactions was compiled from the literature and String protein-protein interactions database v.10.5 (https://string-db.org). Only predictions with medium (or better) confidence (0.4) limited to the top 10 interactions with protein-protein interaction (PPI) enrichment p-value ≤ 0.5 were considered. To compare the different species, protein annotations were standardize by identity to I. scapularis/I. ricinus-D. melanogaster-H. sapiens order of priority.
References
1
AdipietroK. A.MainlandJ. D.MatsunamiH. (2012). Functional evolution of mammalian odorant receptors.PLoS Genet.8:e1002821. 10.1371/journal.pgen.1002821
2
AlmazánC.Blas-MachadoU.KocanK. M.YoshiokaJ. H.BlouinE. F.MangoldA. J.et al (2005). Characterization of three Ixodes scapularis cDNAs protective against tick infestations.Vaccine234403–4416. 10.1016/j.vaccine.2005.04.012
3
AlmazánC.KocanK. M.BergmanD. K.Garcia-GarciaJ. C.BlouinE. F.de la FuenteJ. (2003). Identification of protective antigens for the control of Ixodes scapularis infestations using cDNA expression library immunization.Vaccine211492–1501. 10.1016/S0264-410X(02)00683-7
4
AlmazánC.LagunesR.VillarM.CanalesM.Rosario-CruzR.JongejanF.et al (2010). Identification and characterization of Rhipicephalus (Boophilus) microplus candidate protective antigens for the control of cattle tick infestations.Parasitol. Res.106471–479. 10.1007/s00436-009-1689-1
5
AlmazánC.Moreno-CantúO.Moreno-CidJ. A.GalindoR. C.CanalesM.VillarM.et al (2012). Control of tick infestations in cattle vaccinated with bacterial membranes containing surface-exposed tick protective antigens.Vaccine30265–272. 10.1016/j.vaccine.2011.10.102
6
AnderssonS. G.KurlandC. G. (1998). Reductive evolution of resident genomes.Trends Microbiol.6263–268. 10.1016/S0966-842X(98)01312-2
7
AntunesS.MerinoO.MosquedaJ.Moreno-CidJ. A.Bell-SakyiL.FragkoudisR.et al (2014). Tick capillary feeding for the study of proteins involved in tick-pathogen interactions as potential antigens for the control of tick infestation and pathogen infection.Parasit. Vectors7:42. 10.1186/1756-3305-7-42
8
AyllónN.VillarM.GalindoR. C.KocanK. M.ŠímaR.LópezJ. A.et al (2015). Systems biology of tissue-specific response to Anaplasma phagocytophilum reveals differentiated apoptosis in the tick vector Ixodes scapularis.PLoS Genet.11:e1005120. 10.1371/journal.pgen.1005120
9
BensaciM.BhattacharyaD.ClarkR.HuL. T. (2012). Oral vaccination with Vaccinia virus expressing the tick antigen subolesin inhibits tick feeding and transmission of Borrelia burgdorferi.Vaccine306040–6046. 10.1016/j.vaccine.2012.07.053
10
BeutlerB.MorescoE. M. (2008). Akirins versus infection.Nat. Immunol.97–9. 10.1038/ni0108-7
11
BonnayF.NguyenX. H.Cohen-BerrosE.TroxlerL.BatscheE.CamonisJ.et al (2014). Akirin specifies NF-κB selectivity of Drosophila innate immune response via chromatin remodeling.EMBO J.332349–2362. 10.15252/embj.201488456
12
BrownC. J.JohnsonA. K.DunkerA. K.DaughdrillG. W. (2011). Evolution and disorder.Curr. Opin. Struct. Biol.21441–446. 10.1016/j.sbi.2011.02.005
13
BusbyA. T.AyllónN.KocanK. M.BlouinE. F.de la FuenteG.GalindoR. C.et al (2012). Expression of heat-shock proteins and subolesin affects stress responses, Anaplasma phagocytophilum infection and questing behavior in the tick, Ixodes scapularis.Med. Vet. Entomol.2692–102. 10.1111/j.1365-2915.2011.00973.x
14
CanalesM.NaranjoV.AlmazánC.MolinaR.TsurutaS. A.SzabóM. P. J.et al (2009). Conservation and immunogenicity of the mosquito ortholog of the tick protective antigen, subolesin.Parasitol. Res.10597–111. 10.1007/s00436-009-1368-2
15
CarpioY.BasabeL.AcostaJ.RodríguezA.MendozaA.LispergerA.et al (2011). Novel gene isolated from Caligus rogercresseyi: a promising target for vaccine development against sea lice.Vaccine292810–2820. 10.1016/j.vaccine.2011.01.109
16
CarpioY.GarcíaC.PonsT.HaussmannD.Rodríguez-RamosT.BasabeL.et al (2013). Akirins in sea lice: first steps toward a deeper understanding.Exp. Parasitol.135188–199. 10.1016/j.exppara.2013.06.018
17
CarreónD.Pérez de la LastraJ. M.AlmazánC.CanalesM.Ruiz-FonsF.BoadellaM.et al (2012). Vaccination with BM86, subolesin and akirin protective antigens for the control of tick infestations in white tailed deer and red deer.Vaccine30273–279. 10.1016/j.vaccine.2011.10.099
18
ClemonsA. M.BrockwayH. M.YinY.KasinathanB.ButterfieldY. S.JonesS. J.et al (2013). Akirin is required for diakinesis bivalent structure and synaptonemal complex disassembly at meiotic prophase I.Mol. Biol. Cell.241053–1067. 10.1091/mbc.E12-11-0841
19
ContrerasM.de la FuenteJ. (2016a). Control of Ixodes ricinus and Dermacentor reticulatus tick infestations in rabbits vaccinated with the Q38 Subolesin/Akirin chimera.Vaccine343010–3013. 10.1016/j.vaccine.2016.04.092
20
ContrerasM.de la FuenteJ. (2016b). Vaccinomics approach to the development of vaccines for the control of multiple ectoparasite infestations.Nova Acta Leopold.411185–200. 10.1186/1471-2164-12-105
21
ContrerasM.Moreno-CidJ. A.DomingosA.CanalesM.Díez-DelgadoI.Pérez de la LastraJ. M.et al (2015). Bacterial membranes enhance the immunogenicity and protective capacity of the surface exposed tick Subolesin-Anaplasma marginale MSP1a chimeric antigen.Ticks Tick Borne Dis.6820–828. 10.1016/j.ttbdis.2015.07.010
22
da CostaM.Pinheiro-SilvaR.AntunesS.Moreno-CidJ. A.VillarM.de la FuenteJ.et al (2014). Mosquito Akirin as a potential antigen for malaria control.Malar. J.13:470. 10.1186/1475-2875-13-470
23
DasGuptaR.KaykasA.MoonR. T.PerrimonN. (2005). Functional genomic analysis of the Wnt-wingless signaling pathway.Science308826–833. 10.1126/science.1109374
24
de la FuenteJ. (2018). Controlling ticks and tick-borne diseases…looking forward.Ticks Tick Borne Dis.91354–1357. 10.1016/j.ttbdis.2018.04.001
25
de la FuenteJ.AlmazaìnC.Blas-MachadoU.NaranjoV.MangoldA. J.BlouinE. F.et al (2006a). The tick protective antigen, 4D8, is a conserved protein involved in modulation of tick blood ingestion and reproduction.Vaccine244082–4095.
26
de la FuenteJ.AlmazánC.BlouinE. F.NaranjoV.KocanK. M. (2006b). Reduction of tick infections with Anaplasma marginale and A. phagocytophilum by targeting the tick protective antigen subolesin.Parasitol. Res.10085–91.
27
de la FuenteJ.AntunesS.BonnetS.Cabezas-CruzA.DomingosA.Estrada-PeñaA.et al (2017a). Tick-pathogen interactions and vector competence: identification of molecular drivers for tick-borne diseases.Front. Cell. Infect. Microbiol.7:114. 10.3389/fcimb.2017.00114
28
de la FuenteJ.ContrerasM.Estrada-PeñaA.Cabezas-CruzA. (2017b). Targeting a global health problem: vaccine design and challenges for the control of tick-borne diseases.Vaccine355089–5094. 10.1016/j.vaccine.2017.07.097
29
de la FuenteJ.ContrerasM. (2015). Tick vaccines: current status and future directions.Expert Rev. Vaccines141367–1376. 10.1586/14760584.2015.1076339
30
de la FuenteJ.KocanK. M. (2006). Strategies for development of vaccines for control of ixodid tick species.Parasite Immunol.28275–283. 10.1111/j.1365-3024.2006.00828.x
31
de la FuenteJ.KocanK. M. (2014). “Development of vaccines for control of tick infestations and interruption of pathogen transmission,” inBiology of Ticks, 2nd Edn, edsSonenshineD. E.RoeR. M. (New York, NY: Oxford University Press), 333–352.
32
de la FuenteJ.Manzano-RomanR.NaranjoV.KocanK. M.ZivkovicZ.BlouinE. F.et al (2010). Identification of protective antigens by RNA interference for control of the lone star tick, Amblyomma americanum.Vaccine281786–1795. 10.1016/j.vaccine.2009.12.007
33
de la FuenteJ.Maritz-OlivierC.NaranjoV.AyoubiP.NijhofA. M.AlmazánC.et al (2008). Evidence of the role of tick subolesin in gene expression.BMC Genomics9:372. 10.1186/1471-2164-9-372
34
de la FuenteJ.Moreno-CidJ. A.CanalesM.VillarM.Pérez de la LastraJ. M.KocanK. M.et al (2011). Targeting arthropod subolesin/akirin for the development of a universal vaccine for control of vector infestations and pathogen transmission.Vet. Parasitol.18117–22. 10.1016/j.vetpar.2011.04.018
35
de la FuenteJ.Moreno-CidJ. A.GalindoR. C.AlmazánC.KocanK. M.MerinoO.et al (2013). Subolesin/Akirin vaccines for the control of arthropod vectors and vector-borne pathogens.Transbound. Emerg. Dis. 60(Suppl. 2), 172–178. 10.1111/tbed.12146
36
de la FuenteJ.VillarM.Cabezas-CruzA.Estrada-PeñaA.AyllónN.AlberdiP. (2016). Tick-host-pathogen interactions: conflict and cooperation.PLoS Pathog.12:e1005488. 10.1371/journal.ppat.1005488
37
DingwallC.RobbinsJ.DilworthS. M.RobertsB.RichardsonW. D. (1988). The nucleoplasmin nuclear location sequence is larger and more complex than that of SV-40 large T antigen.J. Cell Biol.107841–849. 10.1083/jcb.107.3.841
38
DriscollT. P.VerhoeveV. I.GuillotteM. L.LehmanS. S.RennollS. A.Beier-SextonM.et al (2017). Wholly Rickettsia! Reconstructed metabolic profile of the quintessential bacterial parasite of eukaryotic cells.mBio8:e859-17. 10.1128/mBio.00859-17
39
EhrenbergM.KurlandC. G. (1984). Costs of accuracy determined by a maximal growth rate constraint.Q. Rev. Biophys.1745–82. 10.1017/S0033583500005254
40
Estrada-PeñaA.VillarM.Artigas-JerónimoS.LópezV.AlberdiP.Cabezas-CruzA.et al (2018). Use of graph theory to characterize human and arthropod vector cell protein response to infection.Front. Cell. Infect. Microbiol.8:265. 10.3389/fcimb.2018.00265
41
FelsensteinJ. (1985). Confidence limits on phylogenies: an approach using the bootstrap.Evolution39783–791. 10.1111/j.1558-5646.1985.tb00420.x
42
GalindoR. C.Doncel-PérezE.ZivkovicZ.NaranjoV.GortazarC.MangoldA. J.et al (2009). Tick subolesin is an ortholog of the akirins described in insects and vertebrates.Dev. Comp. Immunol.33612–617. 10.1016/j.dci.2008.11.002
43
Giraldo-CalderónG. I.EmrichS. J.MacCallumR. M.MaslenG.DialynasE.TopalisP.et al (2015). VectorBase: an updated bioinformatics resource for invertebrate vectors and other organisms related with human diseases.Nucleic Acids Res.43D707–D713. 10.1093/nar/gku1117
44
GonzalezK.BayliesM. (2005). “bhringi: A novel Twist co-regulator”, inProceedings of the Program and Abstracts 46th Annual Drosophila Research Conference, San Diego, CA, 46:320B.
45
GotoA.FukuyamaH.ImlerJ. L.HoffmannJ. A. (2014). The chromatin regulator DMAP1 modulates activity of the nuclear factor B (NF-B) transcription factor Relish in the Drosophila innate immune response.J. Biol. Chem.28920470–20476. 10.1074/jbc.C114.553719
46
GotoA.MatsushitaK.GesellchenV.El ChamyL.KuttenkeulerD.TakeuchiO.et al (2008). Akirins are highly conserved nuclear proteins required for NF-kappaB-dependent gene expression in Drosophila and mice.Nat. Immunol.997–104. 10.1038/ni1543
47
GregoryT. R. (2005). Animal Genome Size Database. Avilable at: http://www.genomesize.com
48
Gulia-NussM.NussA. B.MeyerJ. M.SonenshineD. E.RoeR. M.WaterhouseR. M.et al (2016). Genomic insights into the Ixodes scapularis tick vector of Lyme disease.Nat. Commun.7:10507. 10.1038/ncomms10507
49
HajdušekO.ŠímaR.AyllónN.JaloveckáM.PernerJ.de la FuenteJ.et al (2013). Interaction of the tick immune system with transmitted pathogens.Front. Cell. Infect. Microbiol.3:26. 10.3389/fcimb.2013.00026
50
HarringtonD.CanalesM.de la FuenteJ.de LunaC.RobinsonK.GuyJ.et al (2009). Immunisation with recombinant proteins subolesin and Bm86 for the control of Dermanyssus gallinae in poultry.Vaccine274056–4063. 10.1016/j.vaccine.2009.04.014
51
HolmL.LaaksoL. M. (2016). Dali server update.Nucleic Acids Res.44W351–W355. 10.1093/nar/gkw357
52
HoogstraalH. (1956). African Ixodoidea. I. Ticks of the Sudan (with Special Reference to Equatoria Province and with Preliminary Reviews of the Genera Boophilus, Margaropus and Hyalomma). Washington, DC: Department of the Navy, 1101.
53
HorakI. G.HeyneH.WilliamsR.GallivanG. J.SpickettA. M.BezuidenhoutJ. D.et al (2018). The Ixodid Ticks (Acari: Ixodidae) of Southern Africa. Hamburg: Springer, 675. 10.1007/978-3-319-70642-9
54
HouF.WangX.QianZ.LiuQ.LiuY.HeS.et al (2013). Identification and functional studies of Akirin, a potential positive nuclear factor of NF-κB signaling pathways in the Pacific white shrimp, Litopenaeus vannamei.Dev. Comp. Immunol.41703–714. 10.1016/j.dci.2013.08.005
55
IvicsZ.HackettP. B.PlasterkR. H.IzsvakZ. (1997). Molecular reconstruction of Sleeping Beauty, a Tc1-like transposon from fish, and its transposition in human cells.Cell91501–510. 10.1016/S0092-8674(00)80436-5
56
KasthuriS. R.UmasuthanN.WhangI.WanQ.LimB. S.JungH. B.et al (2013). Akirin2 homologues from rock bream, Oplegnathus fasciatus: genomic and molecular characterization and transcriptional expression analysis.Fish Shellfish Immunol.35740–747. 10.1016/j.fsi.2013.06.006
57
KatohK.RozewickiJ.YamadaK. D. (2017). MAFFT online service: multiple sequence alignment, interactive sequence choice and visualization.Brief. Bioinform.10.1093/bib/bbx108 [Epub ahead of print].
58
KatohK.StandleyD. M. (2013). MAFFT multiple sequence alignment software version 7: improvements in performance and usability.Mol. Biol. Evol.30772–780. 10.1093/molbev/mst010
59
KimH.LeeS. K.HongM. W.ParkS. R.LeeY. S.KimJ. W.et al (2013). Association of a single nucleotide polymorphism in the akirin 2 gene with economically important traits in Korean native cattle.Anim. Genet.44750–753. 10.1111/age.12055
60
KomiyaY.KurabeN.KatagiriK.OgawaM.SugiyamaA.KawasakiY.et al (2008). A novel binding factor of 14-3-3beta functions as a transcriptional repressor and promotes anchorage-independent growth, tumorigenicity, and metastasis.J. Biol. Chem.28318753–18764. 10.1074/jbc.M802530200
61
KooninE. V. (2005). Orthologs, paralogs, and evolutionary genomics.Annu. Rev. Genet.39309–338. 10.1371/journal.pgen.1002821
62
KrossaS.SchmittA. D.HattermannK.FritschJ.ScheidigA. J.MehdornH. M.et al (2015). Down regulation of Akirin-2 increases chemosensitivity in human glioblastomas more efficiently than Twist-1.Oncotarget621029–21045. 10.18632/oncotarget.3763
63
KumarB.ManjunathacharH. V.NagarG.RavikumarG.de la FuenteJ.SaravananB. C.et al (2017). Functional characterization of candidate antigens of Hyalomma anatolicum and evaluation of its cross-protective efficacy against Rhipicephalus microplus.Vaccine355682–5692. 10.1016/j.vaccine.2017.08.049
64
KumarS.StecherG.TamuraK. (2016). MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets.Mol. Biol. Evol.331870–1874. 10.1093/molbev/msw054
65
KurlandC. G.CanbäckB.BergO. G. (2007). The origins of modern proteomes.Biochimie891454–1463. 10.1016/j.biochi.2007.09.004
66
LawM. J.LindeM. E.ChambersE. J.OubridgeC.KatsambaP. S.NilssonL.et al (2006). The role of positively charged amino acids and electrostatic interactions in the complex of U1A protein and U1 hairpin II RNA.Nucleic Acids Res.34275–285. 10.1093/nar/gkj436
67
LeungS. W.HarremanM. T.HodelM. R.HodelA. E.CorbettA. H. (2003). Dissection of the karyopherin alpha nuclear localization signal (NLS)-binding groove: functional requirements for NLS binding.J. Biol. Chem.27841947–41953. 10.1074/jbc.M307162200
68
LiuN.WangX. W.SunJ. J.WangL.ZhangH. W.ZhaoX. F.et al (2016). Akirin interacts with Bap60 and 14-3-3 proteins to regulate the expression of antimicrobial peptides in the kuruma shrimp (Marsupenaeus japonicus).Dev. Comp. Immunol.5580–89. 10.1016/j.dci.2015.10.015
69
LiuT.GaoY.XuT. (2015). Evolution of akirin family in gene and genome levels and coexpressed patterns among family members and rel gene in croaker.Dev. Comp. Immunol.5217–25. 10.1016/j.dci.2015.04.010
70
LiuX.XiaY.TangJ.MaL.LiC.MaP.et al (2017). Dual roles of Akirin2 protein during Xenopus neural development.J. Biol. Chem.2925676–5684. 10.1074/jbc.M117.777110
71
LiuY.SongL.SunY.LiuT.HouF.LiuX. (2016). Comparison of immune response in Pacific white shrimp, Litopenaeus vannamei, after knock down of Toll and IMD gene in vivo.Dev. Comp. Immunol.6041–52. 10.1016/j.dci.2016.02.004
72
LuP.ZhouY.YuY.CaoJ.ZhangH.GongH.et al (2016). RNA interference and the vaccine effect of a subolesin homolog from the tick Rhipicephalus haemaphysaloides.Exp. Appl. Acarol.68113–126. 10.1007/s10493-015-9987-z
73
MacqueenD. (2009). Commentary on Galindo et al. [Dev. Comp. Immunol. 33(4) (2009) 612-617].Dev. Comp. Immunol.33877; authorrely 878–879. 10.1016/j.dci.2009.02.011
74
MacqueenD. J.JohnstonI. A. (2009). Evolution of the multifaceted eukaryotic akirin gene family.BMC Evol. Biol.9:34. 10.1186/1471-2148-9-34
75
MacqueenD. J.BowerN. I.JohnstonI. A. (2010a). Positioning the expanded akirin gene family of Atlantic salmon within the transcriptional networks of myogenesis.Biochem. Biophys. Res. Commun.200599–605. 10.1016/j.bbrc.2010.08.110
76
MacqueenD. J.KristjánssonB. K.JohnstonI. A. (2010b). Salmonid genomes have a remarkably expanded akirin family, coexpressed with genes from conserved pathways governing skeletal muscle growth and catabolism.Physiol. Genomics42134–148. 10.1152/physiolgenomics.00045.2010
77
MaedaI.KoharaY.YamamotoM.SugimotoA. (2001). Large-scale analysis of gene function in Caenorhabditis elegans by high-throughput RNAi.Curr. Biol.11171–176. 10.1016/S0960-9822(01)00052-5
78
MakkerhJ. P.DingwallC.LaskeyR. A. (1996). Comparative mutagenesis of nuclear localization signals reveals the importance of neutral and acidic amino acids.Curr. Biol.61025–1027. 10.1016/S0960-9822(02)00648-6
79
MansB. J.De CastroM. H.PienaarR.De KlerkD.GavenP.GenuS.et al (2016). Ancestral reconstruction of tick lineages.Ticks Tick Borne Dis.7509–535. 10.1016/j.ttbdis.2016.02.002
80
Manzano-RománR.Díaz-MartínV.OleagaA.Pérez-SánchezR. (2015). Identification of protective linear B-cell epitopes on the subolesin/akirin orthologues of Ornithodoros spp. soft ticks.Vaccine331046–1055. 10.1016/j.vaccine.2015.01.015
81
Manzano-RománR.Díaz-MartínV.OleagaA.Siles-LucasM.Pérez-SánchezR. (2012). Subolesin/akirin orthologs from Ornithodoros spp. soft ticks: cloning, RNAi gene silencing and protective effect of the recombinant proteins.Vet. Parasitol.185248–259. 10.1016/j.vetpar.2011.10.032
82
MarshallA.SalernoM. S.ThomasM.DaviesT.BerryC.DyerK.et al (2008). Mighty is a novel promyogenic factor in skeletal myogenesis.Exp. Cell Res.3141013–1029. 10.1016/j.yexcr.2008.01.004
83
MedzhitovR.JanewayC. A.Jr. (1997). Innate immunity: impact on the adaptive immune response.Curr. Opin. Immunol.94–9. 10.1016/S0952-7915(97)80152-5
84
MerinoO.AlberdiP.Pérez de la LastraJ. M.de la FuenteJ. (2013a). Tick vaccines and the control of tick-borne pathogens.Front. Cell. Infect. Microbiol.3:30. 10.3389/fcimb.2013.00030
85
MerinoO.AntunesS.MosquedaJ.Moreno-CidJ. A.Pérez de la LastraJ. M.Rosario-CruzR.et al (2013b). Vaccination with proteins involved in tick-pathogen interactions reduces vector infestations and pathogen infection.Vaccine315889–5896. 10.1016/j.vaccine.2013.09.037
86
MerinoO.AlmazánC.CanalesM.VillarM.Moreno-CidJ. A.Estrada-PeñaA.et al (2011a). Control of Rhipicephalus (Boophilus) microplus infestations by the combination of subolesin vaccination and tick autocidal control after subolesin gene knockdown in ticks fed on cattle.Vaccine292248–2254. 10.1016/j.vaccine.2011.01.050
87
MerinoO.AlmazánC.CanalesM.VillarM.Moreno-CidJ. A.GalindoR. C.et al (2011b). Targeting the tick protective antigen subolesin reduces vector infestations and pathogen infection by Anaplasma marginale and Babesia bigemina.Vaccine298575–8579. 10.1016/j.vaccine.2011.09.023
88
MobleyC. B.FoxC. D.FergusonB. S.AminR. H.DalboV. J.BaierS.et al (2014). L-leucine, beta-hydroxy-beta-methylbutyric acid (HMB) and creatine monohydrate prevent myostatin-induced Akirin-1/Mighty mRNA down-regulation and myotube atrophy.J. Int. Soc. Sports Nutr.11:38. 10.1186/1550-2783-11-38
89
MorelP. C. (2003). Les Tiques d’Afrique et du Bassin Meìditerraneìen (1965–1995). Can Tho: CIRAD-EMVT, 1342.
90
Moreno-CidJ. A.JiménezM.CornelieS.MolinaR.AlarcónP.LacroixM.-N.et al (2011). Characterization of Aedes albopictus akirin for the control of mosquito and sand fly infestations.Vaccine2977–82. 10.1016/j.vaccine.2010.10.011
91
Moreno-CidJ. A.Pérez de la LastraJ. M.VillarM.JiménezM.PinalR.Estrada-PeñaA.et al (2013). Control of multiple arthropod vector infestations with subolesin/akirin vaccines.Vaccine311187–1196. 10.1016/j.vaccine.2012.12.073
92
NaranjoN.AyllónN.Pérez de la LastraJ. M.GalindoR. C.KocanK. M.BlouinE. F.et al (2013). Reciprocal regulation of NF-kB (Relish) and Subolesin in the tick vector, Ixodes scapularis.PLoS One8:e65915. 10.1371/journal.pone.0065915
93
NowakS. J.AiharaH.GonzalezK.NibuY.BayliesM. K. (2012). Akirin links twist-regulated transcription with the Brahma chromatin remodeling complex during embryogenesis.PLoS Genet.8:e1002547. 10.1371/journal.pgen.1002547
94
NowakS. J.BayliesM. K. (2012). Akirin: a context-dependent link between transcription and chromatin remodeling.Bioarchitecture2209–213. 10.4161/bioa.22907
95
OldsC. L.MwauraS.OdongoD. O.ScolesG. A.BishopR.DaubenbergerC. (2016). Induction of humoral immune response to multiple recombinant Rhipicephalus appendiculatus antigens and their effect on tick feeding success and pathogen transmission.Parasit. Vectors9:484. 10.1186/s13071-016-1774-0
96
OrengoC. A.ThorntonJ. M. (2005). Protein families and their evolution-a structural perspective.Annu. Rev. Biochem.74867–900. 10.1146/annurev.biochem.74.082803.133029
97
PavithiranA.BathigeS. D. N. K.KugapreethanR.PriyathilakaT. T.YangH.KimM. J.et al (2018). A comparative study of three akirin genes from big belly seahorse Hippocampus abdominalis: molecular, transcriptional and functional characterization.Fish Shellfish Immunol.74584–592. 10.1016/j.fsi.2018.01.025
98
Peña-RangelM. T.RodriguezI.Riesgo-EscovarJ. R. (2002). A misexpression study examining dorsal thorax formation in Drosophila melanogaster.Genetics1601035–1050.
99
PrudencioC. R.Pérez de la LastraJ. M.CanalesM.VillarM.de la FuenteJ. (2010). Mapping protective epitopes in the tick and mosquito subolesin ortholog proteins.Vaccine285398–5406. 10.1016/j.vaccine.2010.06.021
100
QuF.XiangZ.ZhangY.LiJ.ZhangY.YuZ. (2014). The identification of the first molluscan Akirin2 with immune defense function in the Hong Kong oyster Crassostrea hongkongensis.Fish Shellfish Immunol.41455–465. 10.1016/j.fsi.2014.09.029
101
RahmanM. K.Saiful IslamM.YouM. (2018). Impact of subolesin and cystatin knockdown by RNA interference in adult female Haemaphysalis longicornis (Acari: Ixodidae) on blood engorgement and reproduction.Insects9:E39. 10.3390/insects9020039
102
RichardsonJ. M.CollomsS. D.FinneganD. J.WalkinshawM. D. (2009). Molecular architecture of the Mos1 paired-end complex: the structural basis of DNA transposition in a eukaryote.Cell1381096–1108. 10.1016/j.cell.2009.07.012
103
RioualenC.Da CostaQ.ChetritB.Charafe-JauffretE.GinestierC.BidautG. (2017). HTS-Net: an integrated regulome-interactome approach for establishing network regulation models in high-throughput screenings.PLoS One12:e0185400. 10.1371/journal.pone.0185400
104
SalernoM. S.DyerK.BracegirdleJ.PlattL.ThomasM.SiriettV.et al (2009). Akirin1 (Mighty), a novel promyogenic factor regulates muscle regeneration and cell chemotaxis.Exp. Cell Res.3152012–2021. 10.1016/j.yexcr.2009.04.014
105
SasakiM.Akiyama-OdaY.OdaH. (2017). Evolutionary origin of type IV classical cadherins in arthropods.BMC Evol. Biol.17:142. 10.1186/s12862-017-0991-2
106
SasakiS.YamadaT.SukegawaS.MiyakeT.FujitaT.MoritaM.et al (2009). Association of a single nucleotide polymorphism in akirin 2 gene with marbling in Japanese Black beef cattle.BMC Res. Notes2:131. 10.1186/1756-0500-2-131
107
SchettersT.BishopR.CramptonM.KopáčkP.Lew-TaborA.Maritz-OlivierC.et al (2016). Cattle tick vaccine researchers join forces in CATVAC.Parasit. Vectors9:105. 10.1186/s13071-016-1386-8
108
ShakyaM.KumarB.NagarG.de la FuenteJ.GhoshS. (2014). Subolesin: a candidate vaccine antigen for the control of cattle tick infestations in Indian situation.Vaccine323488–3494. 10.1016/j.vaccine.2014.04.053
109
ShawD. K.WangX.BrownL. J.Oliva ChávezA. S.ReifK. E.SmithA. A.et al (2017). Infection-derived lipids elicit an immune deficiency circuit in arthropods.Nat. Commun.8:14401. 10.1038/ncomms14401
110
SmithA.GuoX.de la FuenteJ.NaranjoN.KocanK. M.KaufmanW. R. (2009). The impact of RNA interference of the subolesin and voraxin genes in male Amblyomma hebraeum (Acari: Ixodidae) on female engorgement and oviposition.Exp. Appl. Acarol.4771–86. 10.1007/s10493-008-9195-1
111
SonenshineD. E.KocanK. M.de la FuenteJ. (2006). Tick control: further thoughts on a research agenda.Trends Parasitol.22550–551. 10.1016/j.pt.2006.09.003
112
StaleyJ. T. (2017). Domain Cell Theory supports the independent evolution of the Eukarya, Bacteria and Archaea and the Nuclear Compartment Commonality hypothesis.Open Biol.7:170041. 10.1098/rsob.170041
113
SunW.HuangH.MaS.GanX.ZhuM.LiuH.et al (2016). Akirin2 could promote the proliferation but not the differentiation of duck myoblasts via the activation of the mTOR/p70S6K signaling pathway.Int. J. Biochem. Cell. Biol.79298–307. 10.1016/j.biocel.2016.08.032
114
TarteyS.MatsushitaK.ImamuraT.WakabayashiA.OriD.MinoT.et al (2015). Essential function for the nuclear protein Akirin2 in B cell activation and humoral immune responses.J. Immunol.195519–527. 10.4049/jimmunol.1500373
115
TarteyS.MatsushitaK.VandenbonA.OriD.ImamuraT.MinoT.et al (2014). Akirin2 is critical for inducing inflammatory genes by bridging IκB-ζ and the SWI/SNF complex.EMBO J.332332–2348. 10.15252/embj.201488447
116
TorinaA.Moreno-CidJ. A.BlandaV.Fernández de MeraI. G.Pérez de la LastraJ. M.ScimecaS.et al (2014). Control of tick infestations and pathogen prevalence in cattle and sheep farms vaccinated with the recombinant Subolesin-Major Surface Protein 1a chimeric antigen.Parasit. Vectors7:10. 10.1186/1756-3305-7-10
117
Valenzuela-MuñozV.Gallardo-EscárateC. (2014). TLR and IMD signaling pathways from Caligus rogercresseyi (Crustacea: Copepoda): in silico gene expression and SNPs discovery.Fish Shellfish Immunol.36428–434. 10.1016/j.fsi.2013.12.019
118
VillarM.MarinaA.de la FuenteJ. (2017). Applying proteomics to tick vaccine development: where are we?Expert Rev. Proteomics14211–221. 10.1080/14789450.2017.1284590
119
VoigtF.WiedemannL.ZulianiC.QuerquesI.SebeA.MatesL.et al (2016). Sleeping Beauty transposase structure allows rational design of hyperactive variants for genetic engineering.Nat. Commun.7:11126. 10.1038/ncomms11126
120
WanF.LenardoM. J. (2010). The nuclear signaling of NF-kappaB: current knowledge, new insights, and future perspectives.Cell Res.2024–33. 10.1038/cr.2009.137
121
WangM.KurlandC. G.Caetano-AnollésG. (2011). Reductive evolution of proteomes and protein structures.Proc. Natl. Acad. Sci U.S.A.10811954–11958. 10.1073/pnas.1017361108
122
WatanabeN.SatohY.FujitaT.OhtaT.KoseH.MuramatsuY.et al (2011). Distribution of allele frequencies at TTN g.231054C > T, RPL27A g.3109537C > T and AKIRIN2 c.∗188G > A between Japanese Black and four other cattle breeds with differing historical selection for marbling.BMC Res. Notes4:10. 10.1186/1756-0500-4-10
123
WolfY. I.KooninE. V. (2013). Genome reduction as the dominant mode of evolution.Bioessays35829–837. 10.1002/bies.201300037
124
XueX.WangL.ChenY.ZhangX.LuoH.LiZ.et al (2014). Identification and molecular characterization of an Akirin2 homolog in Chinese loach (Paramisgurnus dabryanus).Fish Shellfish Immunol.36435–443. 10.1016/j.fsi.2013.12.021
125
YanJ.DongX.KongY.ZhangY.JingR.FengL. (2013). Identification and primary immune characteristics of an amphioxus akirin homolog.Fish Shellfish Immunol.35564–571. 10.1016/j.fsi.2013.05.020
126
YangC. G.WangX. L.WangL.ZhangB.ChenS. L. (2011). A new Akirin1 gene in turbot (Scophthalmus maximus): molecular cloning, characterization and expression analysis in response to bacterial and viral immunological challenge.Fish Shellfish Immunol301031–1041. 10.1016/j.fsi.2011.01.028
127
YangC. G.WangX. L.ZhangB.SunB.LiuS. S.ChenS. L. (2013). Screening and analysis of PoAkirin1 and two related genes in response to immunological stimulants in the Japanese flounder (Paralichthys olivaceus).BMC Mol. Biol.14:10. 10.1186/1471-2199-14-10
128
ZanesiN.BalattiV.RiordanJ.BurchA.RizzottoL.PalamarchukA.et al (2013). A Sleeping Beauty screen reveals NF-kB activation in CLL mouse model.Blood1214355–4358. 10.1182/blood-2013-02-486035
129
ZhangY. (2008). I-TASSER server for protein 3D structure prediction.BMC Bioinformatics9:40. 10.1186/1471-2105-9-40
130
ZivkovicZ.EstevesE.AlmazánC.DaffreS.NijhofA. M.KocanK. M.et al (2010a). Differential expression of genes in salivary glands of male Rhipicephalus (Boophilus) microplus in response to infection with Anaplasma marginale.BMC Genomics11:186. 10.1186/1471-2164-11-186
131
ZivkovicZ.TorinaA.MitraR.AlongiA.ScimecaS.KocanK. M.et al (2010b). Subolesin expression in response to pathogen infection in ticks.BMC Immunol.11:7. 10.1186/1471-2172-11-7
Summary
Keywords
immune response, vaccine, interactome, regulome, phylogeny, tick, Anaplasma phagocytophilum
Citation
Artigas-Jerónimo S, Villar M, Cabezas-Cruz A, Valdés JJ, Estrada-Peña A, Alberdi P and de la Fuente J (2018) Functional Evolution of Subolesin/Akirin. Front. Physiol. 9:1612. doi: 10.3389/fphys.2018.01612
Received
03 July 2018
Accepted
25 October 2018
Published
13 November 2018
Volume
9 - 2018
Edited by
Itabajara Da Silva Vaz Jr., Universidade Federal do Rio Grande do Sul (UFRGS), Brazil
Reviewed by
Sukanya Narasimhan, Yale University, United States; Carlos Alexandre Sanchez Ferreira, Pontifícia Universidade Católica do Rio Grande do Sul, Brazil
Updates
Copyright
© 2018 Artigas-Jerónimo, Villar, Cabezas-Cruz, Valdés, Estrada-Peña, Alberdi and de la Fuente.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: José de la Fuente, jose_delafuente@yahoo.com
†These authors have contributed equally to this work
This article was submitted to Invertebrate Physiology, a section of the journal Frontiers in Physiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.