Apoptosis-Inducing Factor 2 (AIF-2) Mediates a Caspase-Independent Apoptotic Pathway in the Tropical Sea Cucumber (Holothuria leucospilota)

Apoptosis, also known as programmed cell death, is a biological process that is critical for embryonic development, organic differentiation, and tissue homeostasis of organisms. As an essential mitochondrial flavoprotein, the apoptosis-inducing factor (AIF) can directly mediate the caspase-independent mitochondrial apoptotic pathway. In this study, we identified and characterized a novel AIF-2 (HlAIF-2) from the tropical sea cucumber Holothuria leucospilota. HlAIF-2 contains a conserved Pyr_redox_2 domain and a putative C-terminal nuclear localization sequence (NLS) but lacks an N-terminal mitochondrial localization sequence (MLS). In addition, both NADH- and FAD-binding domains for oxidoreductase function are conserved in HlAIF-2. HlAIF-2 mRNA was ubiquitously detected in all tissues and increased significantly during larval development. The transcript expression of HlAIF-2 was significantly upregulated after treatment with CdCl2, but not the pathogen-associated molecular patterns (PAMPs) in primary coelomocytes. In HEK293T cells, HlAIF-2 protein was located in the cytoplasm and nucleus, and tended to transfer into the nucleus by CdCl2 incubation. Moreover, there was an overexpression of HlAIF-2-induced apoptosis in HEK293T cells. As a whole, this study provides the first evidence for heavy metal-induced apoptosis mediated by AIF-2 in sea cucumbers, and it may contribute to increasing the basic knowledge of the caspase-independent apoptotic pathway in ancient echinoderm species.


Introduction
Apoptosis, also known as programmed cell death, is an essential biological process that plays critical roles in embryonic development, organic differentiation, and normal tissue homeostasis in metazoans [1]. During immune responses, apoptosis is also important in the cell-mediated killing mechanism for target cells that are invaded by pathogens [2]. The core functions of apoptosis are considered to be mediated by the classical intrinsic and extrinsic pathways with initiator and executioner caspases [3]. In addition, several caspaseindependent apoptotic pathways play important roles in the immune system, including the release of apoptosis-inducing factor (AIF) from mitochondria, which can induce caspaseindependent peripheral chromatin condensation and large-scale DNA fragmentation in the nucleus [4].
AIFs are essential mitochondrial flavoproteins with multiple cellular functions, including the maintenance of electron transport chain function, the regulation of reactive oxygen

Molecular Cloning and Sequence Analysis of HlAIF-2
By using 3′-/5′-RACE approaches, the full-length cDNA sequence of HlAIF-2 was obtained from H. leucospilota and deposited in GenBank under the accession number OM417064. The open reading frame (ORF) of HlAIF-2 cDNA is 1119 bp in length and is predicted to encode a protein of 372 amino acids (a.a.) ( Figure S1). The calculated molecular weight of HlAIF-2 is 40.94 kDa, and the estimated isoelectric point is 5.50. HlAIF-2 contains a putative C-terminal NLS (residues 291-322) but lacks an N-terminal MLS (Figure S1).
Based on the SMART program, a conserved Pyr_redox_2 domain (residues 11-300) was predicted in the HlAIF-2 a.a. sequence ( Figure 1A). In addition, a casein kinase II phosphorylation site and a protein kinase C phosphorylation site were further indicated in the HlAIF-2 a.a. sequence by the ScanProsite program ( Figure 1A).

Phylogenetic, Homology and Structural Analysis
The results of the phylogenetic analysis revealed that AIFs from multiple animal species were classified into three branches: namely, AIF-1, AIF-2, and AIF-3 ( Figure 1B). The branch of AIF-2s was further separated into two clades: vertebrate AIF-2s and invertebrate AIF-2s. Our newly identified HlAIF-2 was found in the clade of invertebrate AIF-2s and shared a close evolutionary distance with the A. japonicus AIF-2 ( Figure 1B). Multiple alignments of a.a. sequences showed that AIF-2s from different species in echinoderms and vertebrates shared considerably conserved sequences ( Figure 2A). Most, if not all, a.a. which were supposed to interact with FAD and NADH were strongly conserved in AIFs, as precisely mapped in Figure 2A. Additionally, the core consensus for the typical motif GXGXXG was found at two distinct regions of the HlAIF-2 a.a sequence (residues 17-22 and 148-153, Figure 2A). Three-dimensional (3-D) modeling was performed for the vertebrate AIF-2 from humans (Homo sapiens) and the echinoderm AIF-2s from the sea cucumber A. japonicus and H. leucospilota. As shown in Figure 2B-D, the probably conservative binding sites of NAD/FAD were highly comparable based on their 3-D structures. The results of the phylogenetic analysis revealed that AIFs from multiple animal species were classified into three branches: namely, AIF-1, AIF-2, and AIF-3 ( Figure 1B). The branch of AIF-2s was further separated into two clades: vertebrate AIF-2s and invertebrate AIF-2s. Our newly identified HlAIF-2 was found in the clade of invertebrate AIF-2s and shared a close evolutionary distance with the A. japonicus AIF-2 ( Figure 1B). Multiple alignments of a.a. sequences showed that AIF-2s from different species in echinoderms and vertebrates shared considerably conserved sequences ( Figure 2A). Most, if not all, a.a. which were supposed to interact with FAD and NADH were strongly conserved in AIFs, as precisely mapped in Figure 2A. Additionally, the core consensus for the typical motif GXGXXG was found at two distinct regions of the HlAIF-2 a.a sequence (residues 17-22 and 148-153, Figure 2A). Three-dimensional (3-D) modeling was performed for the vertebrate AIF-2 from humans (Homo sapiens) and the echinoderm AIF-2s from the sea cucumber A. japonicus and H. leucospilota. As shown in Figure 2B-D, the probably conservative binding sites of NAD/FAD were highly comparable based on their 3-D structures. (A) A.a. sequence alignment of AIF-2 in nine Deuterostomia species. The conserved a.a. residues of H. leucospilota are presented in WebLogo format, and the conserved and similar a.a. residues between different species are labeled in dark blue and light blue, respectively. Residues that interact with FAD or NAD (in H. leucospilota) are marked as "F" or "N", respectively. Two core consensus sequences of the typical motif "GXGXXG" are boxed in red lines. (B-D) Comparison of the 3-D protein of AIF-2 among human H. sapiens and the sea cucumber A. japonicus and H. leucospilota. Space-filling symbols indicate the most conserved binding sites of FAD (orange) and NAD (blue).

Expression Patterns of HlAIF-2 among Different Tissues
The mRNA expression pattern of HlAIF-2 was analyzed in various tissues and different embryonic developmental stages by qPCR. As shown in Figure 3A, HlAIF-2 mRNA was ubiquitously expressed in all the examined tissues, and the strongest expression was found in the intestine, followed by the transverse vessel, rete mirabile, Cuvierian tubules, esophagus, respiratory tree, body wall, coelomocytes, muscle, gonads, and polian vesicle ( Figure 3A). However, the expression of HlAIF-2 mRNA in the polian vesicle reached half of that in the transverse vessel ( Figure 3A).
x FOR PEER REVIEW 5 of 13

Expression Patterns of HlAIF-2 among Different Tissues
The mRNA expression pattern of HlAIF-2 was analyzed in various tissues and different embryonic developmental stages by qPCR. As shown in Figure 3A, HlAIF-2 mRNA was ubiquitously expressed in all the examined tissues, and the strongest expression was found in the intestine, followed by the transverse vessel, rete mirabile, Cuvierian tubules, esophagus, respiratory tree, body wall, coelomocytes, muscle, gonads, and polian vesicle ( Figure 3A). However, the expression of HlAIF-2 mRNA in the polian vesicle reached half of that in the transverse vessel ( Figure 3A).

Expression Patterns of HlAIF-2 during Embryonic and Larval Development
As shown in Figure 3B, HlAIF-2 mRNA was constitutively expressed in all detected embryonic and larval developmental stages, and the highest expression level was observed at the auricularia stage. After that, the expression level of HlAIF-2 mRNA decreased sharply, reached its bottom at the pentactula stage, and increased again at the juvenile stage. Generally, HlAIF-2 mRNA remained expressed at low levels in the embryonic stages but changed significantly in the larval stages.

HlAIF-2 Expression in Response to Challenges of CdCl2, LPS, and Poly (I:C)
Temporal expression of HlAIF-2 mRNA in the coelomocytes was detected after challenge with cadmium chloride (CdCl2, 20 μM) as an oxidative stress ( Figure 4A). After exposure to CdCl2, the expression of HlAIF-2 was first upregulated with a 13.28-fold change (p < 0.001) at 12 h, followed by a 22.68-fold change (p < 0.001) at 24 h. In a parallel experiment, treatments with lipopolysaccharides (LPS) or polyriboinosinic polyribocytidylic acid [poly (I:C)] did not alter the expression level of HlAIF-2 ( Figure 4B,C).

Expression Patterns of HlAIF-2 during Embryonic and Larval Development
As shown in Figure 3B, HlAIF-2 mRNA was constitutively expressed in all detected embryonic and larval developmental stages, and the highest expression level was observed at the auricularia stage. After that, the expression level of HlAIF-2 mRNA decreased sharply, reached its bottom at the pentactula stage, and increased again at the juvenile stage. Generally, HlAIF-2 mRNA remained expressed at low levels in the embryonic stages but changed significantly in the larval stages.

HlAIF-2 Expression in Response to Challenges of CdCl 2 , LPS, and Poly (I:C)
Temporal expression of HlAIF-2 mRNA in the coelomocytes was detected after challenge with cadmium chloride (CdCl 2 , 20 µM) as an oxidative stress ( Figure 4A). After exposure to CdCl 2 , the expression of HlAIF-2 was first upregulated with a 13.28-fold change (p < 0.001) at 12 h, followed by a 22.68-fold change (p < 0.001) at 24 h. In a parallel experiment, treatments with lipopolysaccharides (LPS) or polyriboinosinic polyribocytidylic acid [poly (I:C)] did not alter the expression level of HlAIF-2 ( Figure 4B,C).

Subcellular Localization of HlAIF-2 in HEK293T Cells
The subcellular location of HlAIF-2 was determined by transfection into HEK293T cells in the presence and absence of Cd 2+ . In the control group without Cd 2+ treatment, HlAIF-2 was located in both the cytoplasm and nucleus of the HEK293T cells. After incubation with CdCl 2 for 24 h, part of the cytoplasmic HlAIF-2 translocated into the nucleus ( Figure 5), indicating that HlAIF-2 tended to transfer into the nucleus during Cd 2+ -induced apoptosis. lenge with cadmium chloride (CdCl2, 20 μM) as an oxidative stress ( Figure 4A). After exposure to CdCl2, the expression of HlAIF-2 was first upregulated with a 13.28-fold change (p < 0.001) at 12 h, followed by a 22.68-fold change (p < 0.001) at 24 h. In a parallel experiment, treatments with lipopolysaccharides (LPS) or polyriboinosinic polyribocytidylic acid [poly (I:C)] did not alter the expression level of HlAIF-2 ( Figure 4B,C).  (I:C) (10 μg/mL). Data are presented as mean ± SE (n = 3), and significant differences are analyzed using one-way ANOVA, and shown as n.s. p > 0.05 and *** p < 0.001.

Subcellular Localization of HlAIF-2 in HEK293T Cells
The subcellular location of HlAIF-2 was determined by transfection into HEK293T cells in the presence and absence of Cd 2+ . In the control group without Cd 2+ treatment, HlAIF-2 was located in both the cytoplasm and nucleus of the HEK293T cells. After incubation with CdCl2 for 24 h, part of the cytoplasmic HlAIF-2 translocated into the nucleus ( Figure 5), indicating that HlAIF-2 tended to transfer into the nucleus during Cd 2+ -induced apoptosis. Figure 5. Subcellular localization of HlAIF-2 in HEK293T cells. "DAPI" represents the DAPI-stained cell nuclei; "EGFP-AIF" represents the EGFP-labeled HlAIF-2 protein; "Merge" represents the combination of cell nuclei and HlAIF-2 protein. Cd 2+ treatment (12 h) could trigger HlAIF-2 nuclear translocation, compared with the "Control" group. The arrows indicate the typical cells located in cytoplasm or translocated into nuclei.

Effects of HlAIF-2 Overexpression on Cell Apoptosis
The function of HlAIF-2 in the mediation of apoptosis was validated by transfection with the pcDNA3.1/HA/HlAIF-2 plasmid in HEK293T cells. As shown in Figure 6, the apoptosis rate was detected by TUNEL assay. The results demonstrated that HlAIF-2 overexpression could significantly induce apoptosis with DNA fragmentation in cell nu- Figure 5. Subcellular localization of HlAIF-2 in HEK293T cells. "DAPI" represents the DAPIstained cell nuclei; "EGFP-AIF" represents the EGFP-labeled HlAIF-2 protein; "Merge" represents the combination of cell nuclei and HlAIF-2 protein. Cd 2+ treatment (12 h) could trigger HlAIF-2 nuclear translocation, compared with the "Control" group. The arrows indicate the typical cells located in cytoplasm or translocated into nuclei.

Effects of HlAIF-2 Overexpression on Cell Apoptosis
The function of HlAIF-2 in the mediation of apoptosis was validated by transfection with the pcDNA3.1/HA/HlAIF-2 plasmid in HEK293T cells. As shown in Figure 6, the apoptosis rate was detected by TUNEL assay. The results demonstrated that HlAIF-2 overexpression could significantly induce apoptosis with DNA fragmentation in cell nuclei ( Figure 6A), and the percentage of apoptosis was 17.05%, 31.54%, and 55.48% in the blank group, control group, and experimental group, respectively ( Figure 6B).

Discussion
It is generally known that invertebrates lack adaptive immunity and so, as an alternative, innate immunity becomes a vital part of their immune system against invading pathogens and environmental stresses [29,30]. Apoptosis is a highly regulated and controlled process that confers advantages for organisms, and AIF is an ancient and conserved apoptotic executor that mediates apoptosis via a caspase-independent mitochondrial pathway [5,6].
Based on their widespread existence in various species, from invertebrates to humans [18,19,31,32], AIFs have been proposed to have an ancient and conserved pyridine nucleotide-disulfide oxidoreductase domain (Pyr_redox domain), which could generate superoxide rather than exhibit antioxidant activity [5]. Our current study found that HlAIF-2 contained a Pyr_redox_2 domain and a deduced C-terminal NLS but lacked a recognizable MLS ( Figure 1A). AIF precursors are usually synthesized in the cytoplasm and then imported into mitochondria through N-terminal MLS. However, in the presence of FAD, the MLS of AIF may be removed by proteolysis to produce mature AIF protein [4]. In contrast, AIF lacking MLS can spontaneously bind FAD and refold into mature AIF with a potential apoptosis-promoting function [4,8]. Normally, AIF performs mitochondrial functions and translocates to the nucleus only under the induction of apoptotic signals [33]. Once DNA damage occurs, the permeability of the mitochondrial outer membrane changes, and AIF is released from mitochondria [4,33]. However, subcellular localization experiments indicated that most HlAIF-2 was localized in both the cytoplasm and nucleus of HEK293T cells without Cd 2+ treatment ( Figure 5). When Cd 2+ was added, HlAIF-2 in the cytoplasm tended to concentrate into the nucleus, indicating that overexpression of HlAIF-2 may increase basic apoptosis. In addition, the similar phenomenon that Cd 2+ treatment (12 h) could trigger AIF nuclear translocation dose-dependently has been observed in rat cells, previously [15,34].
Phylogenetic analysis showed that vertebrate AIF-2s first grouped with AIF-1s and AIF-3s, and then clustered into a branch of invertebrate AIF-2s ( Figure 2B), illustrating that HlAIF-2 is a fairly ancient gene with a conserved structure. The a.a. sequence of HlAIF-2 shared high similarity with AIF-2s in other species (Figure 2A). AIFs are known for their oxidoreductase function, which is endowed by their NADH-and FAD-binding domains [35]. The Pyr_redox_2 domain is actually a smaller NADH-binding domain apoptosis by TUNEL assay. "DAPI" represents the DAPI-stained cell nuclei; "FITC" represents the FITC-stained fractured DNA fragments (marker for apoptosis); "Merged" represents the combination of cell nuclei and fractured DNA fragments. (B) Comparison of apoptosis rates for HEK293T cells in different groups. "Blank" represents the blank group (untransfected HEK293T cells); "Control" represents the control group (HEK293T cells transfected with pcDNA3.1/HA); "AIF" represents the experimental group (HEK293T cells transfected with pcDNA3.1/HA/HlAIF-2). The values are expressed as mean ± SE (n = 3), and significant difference was analyzed by the Student's t-test and shown as *** p < 0.001.

Discussion
It is generally known that invertebrates lack adaptive immunity and so, as an alternative, innate immunity becomes a vital part of their immune system against invading pathogens and environmental stresses [29,30]. Apoptosis is a highly regulated and controlled process that confers advantages for organisms, and AIF is an ancient and conserved apoptotic executor that mediates apoptosis via a caspase-independent mitochondrial pathway [5,6].
Based on their widespread existence in various species, from invertebrates to humans [18,19,31,32], AIFs have been proposed to have an ancient and conserved pyridine nucleotide-disulfide oxidoreductase domain (Pyr_redox domain), which could generate superoxide rather than exhibit antioxidant activity [5]. Our current study found that HlAIF-2 contained a Pyr_redox_2 domain and a deduced C-terminal NLS but lacked a recognizable MLS ( Figure 1A). AIF precursors are usually synthesized in the cytoplasm and then imported into mitochondria through N-terminal MLS. However, in the presence of FAD, the MLS of AIF may be removed by proteolysis to produce mature AIF protein [4]. In contrast, AIF lacking MLS can spontaneously bind FAD and refold into mature AIF with a potential apoptosis-promoting function [4,8]. Normally, AIF performs mitochondrial functions and translocates to the nucleus only under the induction of apoptotic signals [33]. Once DNA damage occurs, the permeability of the mitochondrial outer membrane changes, and AIF is released from mitochondria [4,33]. However, subcellular localization experiments indicated that most HlAIF-2 was localized in both the cytoplasm and nucleus of HEK293T cells without Cd 2+ treatment ( Figure 5). When Cd 2+ was added, HlAIF-2 in the cytoplasm tended to concentrate into the nucleus, indicating that overexpression of HlAIF-2 may increase basic apoptosis. In addition, the similar phenomenon that Cd 2+ treatment (12 h) could trigger AIF nuclear translocation dose-dependently has been observed in rat cells, previously [15,34].
Phylogenetic analysis showed that vertebrate AIF-2s first grouped with AIF-1s and AIF-3s, and then clustered into a branch of invertebrate AIF-2s ( Figure 2B), illustrating that HlAIF-2 is a fairly ancient gene with a conserved structure. The a.a. sequence of HlAIF-2 shared high similarity with AIF-2s in other species (Figure 2A). AIFs are known for their oxidoreductase function, which is endowed by their NADH-and FAD-binding domains [35]. The Pyr_redox_2 domain is actually a smaller NADH-binding domain within a larger FAD-binding domain [9]. The binding sites of NADH and FAD in AIF-2s are highly conserved among different species (Figure 2A). Consistently, the 3-D structure of the HlAIF-2 protein was highly comparable with those of AIF-2s from H. sapiens and A. japonicus ( Figure 2B), including NADH-and FAD-binding domains.
Studies have shown that AIFs are widely distributed in various tissues in mammals [36,37]. In the present study, the transcripts of HlAIF-2 were detected in all the tested tissues, with the highest expression level in the intestine ( Figure 3A). Similarly, the intestine was the tissue with the highest expression level of AIF mRNA in Pacific white shrimps [18]. Previous studies showed that HSP70 could inhibit the nuclear translocation of AIFM1 during hibernation and thermal stimulation in A. japonicus, indicating a potential antiapoptotic response in the intestinal cells of sea cucumbers [19]. On the other hand, AIF-2 is reported to be involved in neural differentiation during embryonic development in vertebrates [37,38]. However, our present study showed that HlAIF-2 expression continuously remained at a low level during embryonic development but increased significantly at the larval stages ( Figure 3B). It is possible that the nerves of the sea cucumber initially formed at the larval stages, as was reported in A. japonicus, with the formation of five radial symmetrical nerve structures at the base of the oral tentacle [39,40]. In addition, the intestine of sea cucumbers gradually matures at the larval stages [41,42], and they need to accumulate nutrition for the transformation of planktonic to benthic lifestyles [21,43]. Hence, the roles of HlAIF-2 in the embryonic and larval stages are speculated to be related to neurogenesis and intestinal development.
With the increasing global attention to marine's sustainable development, an increasing number of studies are focused on the responses of marine animals to polluted marine environments. Among them, heavy metal (e.g., Cd 2+ ) pollution will cause irreversible damage to plants, animals, aquatic life, and humans [30]. Cd 2+ can induce mitochondrial oxidative stress and endoplasmic reticulum stress, consequently leading to apoptosis. Mitochondria release pro-apoptotic proteins via both the caspase-dependent pathways and the caspase-independent pathways, and the mechanism of Cd 2+ -induced apoptosis is believed to be complex [15]. In mammalian cells, oxidative stress caused by Cd 2+ may induce the release of cytochrome c from mitochondria, followed by the activation of intracellular procaspase-9 protein, to form apoptosomes and produce active caspase-9 and caspase-3 proteins in a cascade, which ultimately induce apoptosis [44]. On the other hand, AIFs are released directly from the mitochondrial membrane and translocated to the nucleus to undergo apoptosis when mitochondria suffer from oxidative stress [45]. The mitochondrial apoptotic pathway has also been reported in the invertebrate Pacific oyster Crassostrea gigas [46]. In this study, we explored the caspase-independent mitochondrial apoptotic pathway that was induced by heavy metals in sea cucumbers. The expression of HlAIF-2 in coelomocytes was significantly upregulated by challenge with Cd 2+ (Figure 4A). Combined with subcellular localization experiments, showing that HlAIF-2 concentrated from the cytoplasm to the nucleus during Cd 2+ challenge, it is speculated that a conserved AIF-2-mediated apoptotic pathway could be induced by heavy metal stress in echinoderms.
To date, the roles of AIF-2 overexpression in the induction of apoptosis are still controversial. A previous study showed that overexpressed human AMID (AIF-2) could induce apoptosis in HEK293T cells in a dose-dependent manner, resulting in the condensation of chromatin and the formation of apoptotic bodies [10], while another study found that neither chromatin fragmentation nor protein translocation was observed with AMID overexpression [47]. Recent studies have shown that FSP1 (AIF-2) is an unrecognized antiferroptotic gene [11] that can inhibit the proliferation of lipid peroxides and prevent lipid damage and, consequently, ferroptosis [48]. The present study showed that overexpressed HlAIF-2 could induce apoptosis after transfection into HEK293T cells (Figure 6), suggesting a conversed AIF-2-mediated apoptotic pathway in an ancient echinoderm model.
In conclusion, a novel AIF (HlAIF-2) was identified from the sea cucumber H. leucospilota in this study. With its conserved functional domains and NLS, HlAIF-2 can actively respond to Cd 2+ -induced oxidative stress in coelomocytes by translocation from the mitochondrial membrane into the nucleus. Moreover, the overexpressed HlAIF-2 could induce apoptosis in HEK293T cells, with characteristics of DNA fragmentation in cell nuclei. These results collectively suggested that HlAIF-2 was a conserved apoptotic executor that participated in the mitochondrial apoptotic pathway, that could be induced by heavy metal stresses. In addition, HlAIF-2 mRNA was significantly upregulated during sea cucumber larval development and was expressed in the intestine with the highest expression level, suggesting its role in neurogenesis and intestinal development. Given that apoptosis is a complicated mechanism that can be mediated by multiple pathways, the conservation and specificities of apoptosis in echinoderms still need to be further investigated.

Animals and Tissue Collection
Healthy tropical sea cucumbers (H. leucospilota) weighing 100 ± 10 g were obtained from Daya Bay (Shenzhen, China) and temporarily reared in an aquarium with filtrated and aerated seawater (salinity of 35‰ and temperature of 30 • C) for a week before experiments. Sea cucumbers were dissected on ice, and the tissues were collected, frozen rapidly in liquid nitrogen, and stored at −80 • C until RNA extraction. The coelomic fluids were centrifuged at 1000× g for 10 min at 4 • C to harvest the coelomocytes, which were then kept in 1 mL of TRIzol reagent (Invitrogen, Carlsbad, CA, USA) at −80 • C. All animal experiments were conducted in accordance with the guidelines of the South China Sea Institute of Oceanology, Chinese Academy of Sciences, and this research does not contain any studies with human participants.

Molecular Cloning of HlAIF-2 Full-Length cDNA
Total RNA from the intestine of H. leucospilota was extracted using TRIzol reagent (Invitrogen) and reverse-transcribed to synthesize first-strand cDNA, using the PrimeScript™ II 1st Strand cDNA Synthesis Kit (Takara, Kusatsu, Japan). To obtain corresponding full-length cDNA sequences, 3 -and 5 -rapid amplification of cDNA ends (RACE) was performed using the 3 Full Race Core Set Ver. 2.0 and 5 Full Race Kit (Takara) with gene-specific primers (3 RACE1/3 RACE2 and 5 RACE1/5 RACE2, Table S1), respectively, which were designed based on a partial sequence for the HlAIF-2 homolog from a transcriptomic library of H. leucospilota coelomocytes previously constructed by our lab [49].

Detection of HlAIF-2 Transcript by Real-Time PCR
Total RNA was extracted with TRIzol reagent, digested with gDNA Eraser (Takara), and reverse-transcribed using the PrimeScript™ RT Reagent Kit (Takara) for quantitative PCR (qPCR). Specific primers QHlAIF-2-F and QHlAIF-2-R (Table S1) were designed based on the obtained HlAIF-2 cDNA sequences. qPCRs were performed using SYBR Premix Ex Taq™ II (Takara) in a final volume of 20 µL, with the conditions of 40 cycles of 95 • C for 5 s and 60 • C for 30 s. In this experiment, Hlβ-actin was used as an internal control to verify qPCR results.

Plasmid Construction, Cell Line Culture and Transfection
The coding region of HlAIF-2 was amplified by PCR, using the gene-specific primers PHlAIF-2-F and PHlAIF-2-R (Table S1) and subcloned into the expression vectors pEGFP-N1 (Promega, Madison, WI, USA) and pcDNA3.1/HA (Invitrogen) by homologous recombination using a Hieff Clone Plus One Step Cloning Kit (Yeasen, Shanghai, China). All the plasmids used for transfection were extracted from overnight bacterial cultures using a Plasmid MiniPrep DNA Kit (Axygen, Union City, CA, USA), and all constructed recombinant plasmids were subsequently verified by DNA sequencing.

Subcellular Localization and Translocation
After transfection for 24 h, HEK293 cells transfected with pEGFP-N1/HlAIF-2 plasmid were, then, treated with CdCl 2 (20 µM) for 12 h. As a control, HEK293 cells were transfected with pEGFP-N1/HlAIF-2 plasmid and cultured for 36 h without CdCl 2 treatment. Subsequently, HEK293 cells were rinsed with PBS for 5 min, fixed with precooled 4% paraformaldehyde for 10 min, rinsed again with PBS for 5 min, treated with 0.5% Triton X-100 for 10 min, and stained with DAPI (1 mg/mL) for 10 min in the dark. Finally, HEK293 cells transfected with fluorescent vectors were directly observed by a confocal fluorescence microscope (Leica, Wetzlar, Germany).

Detection of Apoptosis
After transfection for 4 h, the untransfected HEK293 cells (blank group) and HEK293 cells that were transfected with pcDNA3.1/HA (control group) or pcDNA3.1/HA/HlAIF-2 (experimental group) were cultured in 10 mL of DMEM containing 10 % FCS at 37 • C with 5 % CO 2 for 48 h. Then, the apoptotic cells in the three groups were observed by a terminal deoxynucleotidyl transferase (TdT)-mediated dUTP nick-end labeling (TUNEL) assay, as described previously [23].

Data Transformation and Statistical Analysis
All data are presented as the mean ± standard error (SEM). Statistical analysis was performed using one-way ANOVA followed by Fisher's least significant difference (LSD) test with SPSS 22.0 (IBM Software, Armonk, NY, USA), and statistical significance was determined at n.s. p > 0.05, * p < 0.05, ** p < 0.01, and *** p < 0.001.

Supplementary Materials:
The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/ijms23063008/s1. Author Contributions: T.C. and C.R. conceived and designed the experiments; X.L., X.W. and Z.E. performed the experiments; X.L., T.C., X.J. and P.L. analyzed the data; T.C., X.J., P.L., C.H. and C.R. contributed reagents/materials/analysis tools; X.L., T.C. and C.R. wrote the paper. All authors have read and agreed to the published version of the manuscript.