Abstract

In most organisms, superoxide dismutases (SODs) are among the most effective antioxidant enzymes that regulate the reactive oxygen species (ROS) generated by oxidative energy metabolism. ROS are considered main proximate causes of aging. However, it remains unclear if SOD activities are associated with organismal longevity. The queens of eusocial insects, such as termites, ants, and honeybees, exhibit extraordinary longevity in comparison with the nonreproductive castes, such as workers. Therefore, the queens are promising candidates to study the underlying mechanisms of aging. Here, we found that queens have higher Cu/Zn-SOD activity than nonreproductive individuals of the termite Reticulitermes speratus. We identified three Cu/Zn-SOD sequences and one Mn-SOD sequence by RNA sequencing in R. speratus. Although the queens showed higher Cu/Zn-SOD activity than the nonreproductive individuals, there were no differences in their expression levels of the Cu/Zn-SOD genes RsSOD1 and RsSOD3A. Copper (Cu2+ and Cu+) is an essential cofactor for Cu/Zn-SOD enzyme activity, and the queens had higher concentrations of copper than the workers. These results suggest that the high Cu/Zn-SOD activity of termite queens is related to their high levels of the cofactor rather than gene expression. This study highlights that Cu/Zn-SOD activity contributes to extraordinary longevity in termites.

1. Introduction

Queens of termites, ants, and honeybees are representative eusocial insects that live up to 10 times longer than nonreproductive workers, although they possess the same genome [13]. Generally, longevity negatively correlates with reproduction among species [4]. Most animal species show a gradual decline in reproduction with age [5]. Indeed, sterile female flies have longer lifespans than fertile females [6], and germline ablation extends worm longevity [7]. However, the queens of eusocial insects are believed to live for long periods while also laying large numbers of eggs [8]. Because of their unique characteristics, eusocial insects are expected to facilitate the discovery of the mechanisms involved in longevity determination [9].

In most organisms, superoxide dismutases (SODs) are the first line of antioxidant defense against the toxic reactive oxygen species (ROS) that are generated by aerobic metabolism. SODs convert the toxic superoxide anion radicals into hydrogen peroxide, which is subsequently detoxified by catalase, glutathione peroxidase, and peroxiredoxin. At present, two types of Cu/Zn-SODs, which are localized to either cytoplasmic compartments (SOD1) or extracellular elements (SOD3), and Mn-SOD (SOD2), which is exclusively located in mitochondria, have been identified in a variety of organisms [10, 11]. The accumulation of ROS-mediated damage is associated with aging and negative effects on longevity [12, 13]. Therefore, SODs contribute to stress resistance associated with the lifespan of an organism. Indeed, many studies have reported that SODs contributed to the longevity of the fruit fly Drosophila melanogaster [1416], the yeast Saccharomyces cerevisiae [17], and the nematode Caenorhabditis elegans [18]. However, there are also many conflicting reports [1922]. The impact of high SOD activity levels on longevity remains to be determined.

Unlike short-lived organisms, such as fruit flies, yeasts, and nematodes, termites sustain their maturity for a long time and sexually reproduce each year, which leads to an extended period of strong selection that should promote long life for the reproductives (queens and kings). Therefore, termite reproductives may be suitable model organisms to investigate if SODs are proximate causes of longevity. Reticulitermes speratus is one of the most studied subterranean termites, with regard to its reproductive system [23], symbiotic system [24], pheromone communication [25, 26], and antioxidant system [27, 28]. In particular, a previous study revealed that R. speratus produced numerous neotenic queens by an asexual queen succession system [29], which results in adequate sample numbers for several experiments. Hence, we selected the termite R. speratus as a model organism.

In this study, we investigated if high expression of SODs is associated with the extreme lifespan of the queens of R. speratus. First, we identified three Cu/Zn-SOD sequences and one Mn-SOD sequence from the RNA sequencing (RNA-seq) data in R. speratus. Then, we tested if the levels of SODs are higher in termite queens than in nonreproductive individuals.

2. Materials and Methods

2.1. Termites

Animal ethics committee approval was not required for this study, which used insect species. Eight colonies of the termite R. speratus (workers, soldiers, nymphs, and mature neotenic queens) were collected from the experimental forest of Yamaguchi University, which is part of Mt. Himeyama in Yamaguchi, western Japan. We used pooled samples from different colonies for each experiment, as described in Table S1. These insect samples were preserved at −80°C until use.

2.2. Molecular Analysis of R. speratus SODs

The whole transcriptome of R. speratus was examined using next-generation RNA-seq technology in a previous study [26]. We obtained mRNA sequences of SOD genes from the transcriptome data by performing a BLAST search with the amino acid sequences of translated SOD genes in the termite Zootermopsis nevadensis and other model organisms (Table S2). The presence and location of signal peptide cleavage sites in SOD proteins were predicted using the SignalP 4.1 server [30] and TargetP 1.1 server [31, 32]. We performed multiple amino acid sequence alignments with CLUSTALW and conducted phylogenetic analyses using the molecular evolutionary genetics analysis software MEGA7 [33]. Gene evolutionary history was inferred using the maximum likelihood method based on the Whelan and Goldman model [34], which is the best model based on the Bayesian information criterion.

2.3. Protein Extraction

Whole insect body samples that had been stored at −80°C were ground to powder in liquid nitrogen, then homogenised by sonication in tubes with 20 mM Tris–HCl containing 2% protease inhibitor cocktail (v/v), and followed by centrifugation at 17,000 ×g for 10 min at 4°C. The supernatants containing proteins were transferred to new tubes and used as samples for analyses. For each sample, protein concentration was measured with a bicinchoninic acid assay kit before extraction. The protein samples were preserved at −80°C until use in antioxidant activity assays.

2.4. SOD Activity Assays

The activities of antioxidant enzymes were determined as in a previous report [35]. Briefly, we quantified SOD activity by using 2-(-4-iodophenyl)-3-(4-nitrophenyl)-5-(2,4-disulfophenyl)-2H-tetrazolium (WST-1; Dojindo) to detect superoxide anion radicals. The reaction mixture contained diluted xanthine oxidase (approximately 0.2 units), 0.1 mM xanthine, 25 μM WST-1, 0.1 mM EDTA, and 50 mM Na2CO3 (pH 10.2) in a total volume of 3 mL. The absorbance at 438 nm was monitored at 30°C for 1 min. One unit was defined as the amount of enzyme required to inhibit 50% of an absorbance change of 0.06 per minute, which was equivalent to 0.8 units as determined by the standard procedure using cytochrome c. Mn-SOD activities were defined as 1 mM NaCN-resistant activity. We performed 6–9 biological replicates for queens and nonreproductive individuals of R. speratus.

2.5. Measurement of Copper Concentrations

Copper (Cu2+ and Cu+), which is a cofactor of Cu/Zn-SOD, was extracted from termite queens and workers with 0.1 M HCl and measured using Metallo Assay Copper LS (Metallogenics) according to the accompanying manual. The measurements were based on the DiBr-PAESA method. The concentrations of copper were calculated as the molecular weights of copper relative to the protein weights. Sixteen biological replicates were performed, each with a queen and four workers.

2.6. Quantitative Real-Time PCR

We designed primer pairs for each the SOD genes using Primer3 (version 1.1.4; [36]; Table S3). Using ISOGEN reagent (Nippon Gene), total RNA was extracted individually from the whole bodies of termite workers, soldiers, nymphs, and queens that had been frozen with liquid nitrogen and stored at −80°C until extraction. Immediately following extraction, cDNA was synthesized from the RNA using a PrimeScript™ RT reagent kit (Takara) and preserved at −20°C. Quantitative real-time PCR (qRT-PCR) was performed using a LightCycler® (Roche) with QuantiTect® SYBR® Green PCR reagents (Qiagen). All procedures were performed in accordance with each manufacturer’s protocol. GAPDH was selected as the reference gene. Relative expression levels were calculated using a typical ∆∆Ct method. We performed 9–12 biological replicates for queens and nonreproductive individuals of R. speratus.

2.7. Statistical Analysis

R software package (version 3.2.2) was used for all statistical analyses. We performed unpaired t-tests followed by value corrections using Holm’s method for multiple comparisons. All data in graphs are presented as the mean ± standard error of the mean (SEM), and all calculated values are provided in the figure legends. Differences were considered significant when the value was <0.05.

3. Results

3.1. Identification of Cu/Zn-SOD and Mn-SOD Sequences in R. speratus

We previously examined the whole transcriptome of R. speratus using next-generation RNA-seq technology [26]. Via a subsequent homology search of the known amino acid sequence data from the termite Z. nevadensis and several model insects, we inferred the existence of three Cu/Zn-SOD sequences and one Mn-SOD sequence in the R. speratus transcriptome (Table S2). We determined, by sequence alignments with the well-characterized SOD1 and SOD3 of Homo sapiens, Mus musculus, C. elegans, and D. melanogaster, that the R. speratus Cu/Zn-SODs had conserved all of the active sites that are essential for metal binding and disulphide formation [37], suggesting that they are active Cu/Zn-SODs (Figure 1). We also confirmed that the R. speratus Mn-SOD has a conserved manganese-binding site (Figure S1). Previous reports concluded that insects and other organisms have cytoplasmic and extracellular SODs [10, 38]. Our phylogenetic analysis revealed that one of the R. speratus sequences clustered with the cytoplasmic Cu/Zn-SODs of other organisms, whereas the other two sequences clustered with extracellular SODs (Figure S2). Furthermore, we determined that the predicted extracellular sequences have signal peptide cleavage sites (Figure 1 and Table 1). Although these phylogenetic analyses indicated that R. speratus has three Cu/Zn-SODs, we discarded RsSOD3B (FX985482) because of low reliability, based on the high E-value even in comparison to Z. nevadensis (Table S2). Thus, we treated RsSOD1 (FX985484) as the cytoplasmic Cu/Zn-SOD gene, RsSOD3A (FX985481) as the extracellular Cu/Zn-SOD gene, and RsSOD2 (FX985483) as the Mn-SOD gene in our analysis.

3.2. Termite Queens Have High Cu/Zn-SOD Activity in Comparison with Nonreproductive Individuals

Previously, we determined that R. speratus queens showed markedly lower levels of oxidative damage in comparison with nonreproductive workers, which is partly due to high catalase and peroxiredoxin expression and activity [27]. Therefore, we investigated if termite queens had higher SOD activities than nonreproductive individuals. We showed that Cu/Zn-SOD activity in R. speratus queens was higher than in nonreproductive individuals (Figure 2(a)). However, we found no differences in the expression of RsSOD1 and RsSOD3A, which encode cytoplasmic Cu/Zn-SOD and extracellular Cu/Zn-SOD, respectively, between the queens and nonreproductive individuals (Figures 2(b) and 2(c)). We next examined if the copper concentrations corresponded to the Cu/Zn-SOD activity in termites. We found that termite queens have higher copper concentrations than nonreproductive workers (Figure 3). On the other hand, the Mn-SOD activity in R. speratus queens differed only from that in soldiers (Figure 4(a)). In addition, the queens and nonreproductive individuals had similar levels of RsSOD2 expression (Figure 4(b)).

4. Discussion

To date, the determining factors of organismal longevity are among the most fascinating problems for many scientists. Eusocial insects have drawn attention as models for the study of the roles of specific genes in the aging process [9]. In the present study, we found that termite queens have higher Cu/Zn-SOD activity than nonreproductive individuals in the subterranean termite R. speratus. Given the long lifespan of termite queens, these findings suggest that high Cu/Zn-SOD activity may contribute to longevity. Our findings contract with those of earlier studies of eusocial Hymenoptera, which suggested that queen longevity was not associated with cytoplasmic Cu/Zn-SOD activity [39, 40]. In Hymenoptera, Cu/Zn-SOD is reportedly highly expressed in venom glands [41] and secreted in venom [38]. This unique trait may obscure a relationship between SOD levels and longevity. Although eusocial insects have attracted considerable attention for aging research, termites have not yet been investigated. This study highlights the importance of studying eusocial insects like termites, in addition to Hymenoptera, in the search for the mechanisms that contribute to the extraordinarily long lifespan of insect queens.

Oral fluid exchange (trophallaxis) is important for nutritional dynamics and communication in termites because queens depend on nestmate workers for most nutrients. Interestingly, we found that termite queens have higher copper concentrations than workers. Subterranean termites acquire several metal ions, including copper, from the soil and other foods [42]. This implies that higher copper concentrations in termite queens occur by trophic accumulation in castes. High copper concentrations are associated with increased Cu/Zn-SOD activity [43]. In the present study, we demonstrated that there were no significant differences in RsSOD1 and RsSOD3A expression among the termite castes although the queens had higher Cu/Zn-SOD activity than nonreproductive individuals. There are often discrepancies between enzyme activity and gene expression levels. A previous report indicated that copper is critical for Cu/Zn-SOD activity and modulates enzymatic activity in the absence of changes in gene expression [43]. Therefore, the difference in copper concentrations between the queens and workers may be critical for their disparate levels of Cu/Zn-SOD activity. These results indicate that trophallaxis (from workers to queens) not only reduces the energy cost of foraging in termite queens but also provides nutritional benefits to the queens. Thus, trophallaxis may play an important role in nutrient accumulation contributing to termite queen longevity.

It remains unclear if increases in Cu/Zn-SOD activity are localized to specific tissues in termite queens. In mammals, ROS and SODs are thought to play important roles in several aspects of reproductive physiology because superoxide is generated and SOD is expressed in the ovary [44]. Additionally, several studies have reported a relationship between reproduction and antioxidant systems in insects [4547]. Therefore, Cu/Zn-SOD may have an important role to play in reproductive processes, as well as in longevity, in termites. These findings suggest that Cu/Zn-SOD activity may be higher in reproductive tissues than in other somatic tissues. To determine if Cu/Zn-SOD is associated with reproduction in termites, future studies featuring tissue-specific analyses of Cu/Zn-SOD are required. Moreover, it is estimated that copper tends to accumulate in the reproductive tissues in termite queens. It is also necessary to determine the mechanisms by which copper transfer and accumulation occur in termite.

Conflicts of Interest

The authors declare no competing financial interests.

Authors’ Contributions

Eisuke Tasaki and Yoshihito Iuchi conceived and designed the study. Eisuke Tasaki, Kazuya Kobayashi, and Kenji Matsuura collected the biological material. Kazuya Kobayashi performed RNA-seq analysis and Eisuke Tasaki constructed the phylogenetic tree. Eisuke Tasaki performed the enzymatic activity assays, qRT-PCR, and other biochemical analyses. Eisuke Tasaki drafted the manuscript with contributions from Kazuya Kobayashi, Kenji Matsuura, and Yoshihito Iuchi.

Acknowledgments

This work was supported by Japan Society for the Promotion of Science Grant numbers 26660113 to Yoshihito Iuchi and 25221206 to Kenji Matsuura.

Supplementary Materials

Supplementary 1. Figure S1: alignment and structure conservation of Mn-SOD sequences. Multiple sequence alignment of Homo sapiens, Mus musculus, Caenorhabditis elegans, Drosophila melanogaster, Zootermopsis nevadensis, and Reticulitermes speratus. Conserved functional residues are highlighted in bold characters in grey (circled M = manganese-binding site).

Supplementary 2. Figure S2: maximum likelihood molecular phylogenetic tree of Cu/Zn-SOD sequences. We inferred the evolutionary history of the Cu/Zn-SODs using the maximum likelihood method based on the Whelan and Goldman model. The tree with the highest log likelihood (−2592.76) is shown. The percentage of trees in which the associated taxa clustered together is shown above the branches. Initial tree(s) for the heuristic search were obtained automatically by applying Neighbor-Joining and BioNJ algorithms to a matrix of pairwise distances estimated using a JTT model and then selecting the topology with a superior log likelihood value. A discrete gamma distribution was used to model the evolutionary rate differences among the sites (5 categories [+G, parameter = 1.8410]). The rate variation model allowed for some sites to be evolutionarily invariable ([+/], 21.40% sites). The tree is drawn to scale, with branch lengths measured in the number of substitutions per site. The analysis involved 14 amino acid sequences. All positions containing gaps and missing data were eliminated. There were a total of 137 positions in the final dataset. Evolutionary analyses were conducted in MEGA7.

Supplementary 3. Table S1: termite samples.

Supplementary 4. Table S2: target gene details.

Supplementary 5. Table S3: primer sequences.