Genetic Basis of Hydrothermal Vent Adaptation in Bythograeidae Crabs: Insights from Adaptive Evolution of Mitochondrial Protein Coding Genes

Zhengfei Wang*, Dan Tang, Xuejia Shi, Huayun Guo, Xiuping Chen, Daizheng Zhang and Boping Tang*

Jiangsu Key Laboratory for Bioresources of Saline Soils, Jiangsu Synthetic Innovation Center for Coastal Bio-Agriculture, Jiangsu Provincial Key Laboratory of Coastal Wetland Bioresources and Environmental Protection, School of Ocean and Biological Engineering, Yancheng Teachers University, Yancheng 224001, Jiangsu Province, China

ABSTRACT

Deep-sea hydrothermal vents are very extreme environments. Bythograeidae crabs are considered the most common and abundant species and predators at the top of food chain in the hydrothermal vent ecosystem. However, the genetic basis for the adaptation of hydrothermal vents crabs to the harsh environment remains poorly explored. The objectives of this study were to increase our understanding of the mechanisms of hypoxia adaptability in vent crabs and to confirm if there is a correlation between mitochondrial protein coding genes (PCGs) and adaptation to the extreme hydrothermal vent environment. Thirteen PCGs from mitochondrial genomes of 48 Brachyura species and one Diogenidae species were examined. Each of the genes was investigated and compared to orthologous sequences using PAML, Datamonkey, and TreeSAAP. Nine mitochondrial PCGs (ATP6, ATP8, COX1, COX3, CYTB, ND1, ND2, ND4, and ND5) were validated to have undergone positive selection (i.e., directional selection) in vent crabs by at least two methods. A series of putatively selected codons was localized in or close to the important functional regions (protein binding region and helical transmembrane region) in the mitochondrial protein structure. These results help explain why Bythograeidae crabs are capable of living in the hydrothermal vents and suggest that these crabs might have acquired an enhanced capacity for energy metabolism in an extreme hypoxic environment. These findings highlight the critical role of PCGs in the evolution of extreme environmental tolerance by Bythograeidae crabs.

Article Information

Received 07 July 2018

Revised 12 November2018

Accepted 20 December 2018

Available online 19 June 2019

Authors’ Contribution

ZFW, HYG and ZDZ designed and conceived the experiment. ZFW, SXJ, DT, XPC and BPT performed the data analysis and drafted the manuscript.

Key words

Adaptive evolution, Bythograeidae crabs, Hydrothermal vents, Mitochondrial PCGs, Positive selection.

DOI: http://dx.doi.org/10.17582/journal.pjz/2019.51.5.1721.1731

* Corresponding authors: wangzf@yctu.edu.cn;

boptang@163.com

0030-9923/2019/0005-1721 $ 9.00/0

Copyright 2019 Zoological Society of Pakistan

Introduction

Deep-sea hydrothermal vents are very unusual environments that have attracted the attention of scientists since their discovery on the Galapagos Ridge in 1976 (Corliss, 1977; Jannasch and Mottl, 1985; Bettencourt et al., 2010; Sun et al., 2017). Hydrothermal vents are considered short-lived habitats because their temperature is much warmer than that of the surrounding deep-sea (Childress and Fisher, 1992). Hydrothermal vents are characterized by darkness, lack of photosynthetically derived nutrients, low oxygen concentrations, high hydrostatic pressure, and high concentrations of heavy metals and other toxic substances (Hourdez and Lallier, 2007; Martin et al., 2008; Sun et al., 2017). Hydrothermal vent communities illustrate an apparent paradox of displaying a very high biomass despite very challenging environmental conditions when compared to the surrounding deep-sea fauna (Hourdez and Lallier, 2007). Hydrothermal vent ecosystems support dense populations of macrobenthos which are fueled by simple reduced molecules, such as methane and hydrogen sulfide, with the help of chemoautotrophic endosymbionts (Hourdez and Lallier, 2007; Sun et al., 2017). Thus, the organisms living in hydrothermal vents have developed adaptive mechanisms to survive in this hostile environment (Min et al., 2017).

Brachyuran crabs (Bythograeidae) are biomass dominant and many are known to inhabit deep-sea hydrothermal vents (Min et al., 2017). These crabs are considered the most common and abundant species of predator at the top of food chain in the vent ecosystem (Kim et al., 2013). Survival in such extreme conditions requires unique anatomical and physiological adaptations (Bettencourt et al., 2010). These crabs are characterized with specific physiological and morphological traits including reduced eyestalks with vestigial cornea and whole white body (Min et al., 2017). However, studies on the genetic basis for adaptation by vent crabs to the harsh hydrothermal vent environment remain insufficient.

Mitochondria consume the greatest amount (85–90%) of oxygen in cells for oxidative phosphorylation (OXPHOS), which is the primary metabolic pathway for ATP production (Solaini et al., 2010). Five protein complexes in mitochondria are prerequisite for OXPHOS, including nicotinamide adenine dinucleotide (NADH) dehydrogenase (complex I), succinate dehydrogenase (complex II), cytochrome bc1 complex (complex III), cytochrome c oxidase (complex IV), and ATP synthase (complex V) (Eubel et al., 2004; Carroll et al., 2009). The enzymes involved in these complexes are encoded by both nuclear and mitochondrial genes (Wang et al., 2017). Among them, 13 enzymes are encoded by mitochondrial DNA (mtDNA), including seven NADPH reductase subunits (ND1, ND2, ND3, ND4, ND4L, ND5, and ND6), three cytochrome oxidase subunits (COX1, COX2, and COX3), two ATP synthetic subunits (ATP6 and ATP8), and cytochrome b (Cytb) (Kulawiec et al., 2008). Hypoxia is a major threat to the OXPHOS pathway, and if oxygen levels are very low, insufficient ATP availability will result in cell death (Santore et al., 2002). Therefore, mitochondrial protein coding genes are typically chosen to explore the molecular basis of energy metabolism and tolerance to hypoxia (Ning et al., 2010; Shen et al., 2010; Yang et al., 2014; Shi et al., 2018; Wang et al., 2017; Guo et al., 2018).

Thus, the objectives of this study were to increase our understanding of the mechanisms of hypoxia adaptability in vent crabs and to confirm if there was a correlation between mitochondrial protein coding genes (PCGs) and adaptation to the extreme hydrothermal vent environment by comparing the sequences of these genes in Bythograeidae with those in other crabs. Our results demonstrate that a series of positive selection events occurred specifically on the mitochondrial PCGs and on the ancestral and interior lineages of Bythograeidae, suggesting that adaptive evolution of these genes was necessary for Bythograeidae crabs to adapt to the hydrothermal vent environment.

 

Materials and methods

Mitochondrial PCGs and primary treatments

The analysis was comprised of 48 Brachyura species from 19 families (i.e., Homolidae, Leucosiidae, Portunidae, Mithracidae, Majidae, Xanthidae, Menippidae, Bythograeidae, Macrophthalmidae, Varunidae, Grapsidae, Ocypodidae, Dotillidae, Xenograpsidae, Sesarmidae, Mictyridae, Parathelphusidae, Potamidae, and Raninidae), and one Diogenidae species (Clibanarius infraspinatus) as the outgroup (Supplementary Table I). We downloaded the full-length coding sequence (CDS) of 13 mitochondrial PCGs from the National Center for Biotechnology Information (NCBI: https://www.ncbi.nlm.nih.gov/). We used two alignment methods (i.e. CLUSTAL and MUSCLE) as implemented in MEGA 7.0 (Kumar et al., 2016) to align the nucleotide sequences of each mitochondrial PCG and verified them by visual inspection.

Phylogenetic reconstruction

A phylogenetic tree of 48 Brachyura species and one outgroup C. infraspinatus (Supplementary Table I) was reconstructed from concatenated 13 mitochondrial PCGs using the Bayesian inference (BI) and maximum likelihood (ML) algorithms. The BI and ML trees were reconstructed using MrBayes 3.2.6 (Huelsenbeck and Ronquist, 2001; Ronquist et al., 2012) and RaxML (Stamatakis, 2014), respectively. jModeltest (Darriba et al., 2012) was used to select the best model, and the SYM+I+G model was optimized to analyze nucleotide alignment. MtArt+I+G was the appropriate model for the amino acid sequence dataset according to ProtTest 3.4 (Abascal et al., 2005; Darriba et al., 2011) based on Akaike’s information criterion. Two simultaneous runs of 10,000,000 generations were conducted for the matrix in the BI analysis. We made two simultaneous runs, sampling trees every 1,000 generations, with three heated and one cold chain to encourage swapping among the Markov-chain Monte Carlo chains. Convergence of the sampled parameters and potential autocorrelation (effective sampling size/ESS for all parameters > 200) were investigated in Tracer 1.6 (http://tree.bio.ed.ac.uk/software/tracer/). Additionally, the average standard deviation of split frequencies between both runs was checked (<0.01). The Bayesian posterior probabilities were obtained from the 50% majority rule consensus of the post-burn-in trees sampled at stationarity, after removing the first 25% of trees as a “burn-in” stage. The resulting phylogenetic trees were visualized in FigTree v1.4.3 (http://tree.bio.ed.ac.uk/software/figtree/).

Molecular evolutionary analysis

The non-synonymous (dN) to synonymous (dS) rate ratio ω is a measure of selective pressure, with values of ω = 1, > 1, and < 1 indicating neutral selection, positive selection, and purifying selection, respectively (Ohta, 1992). The PAML package (Yang, 2007) was used to determine whether adaptive evolution might have occurred in the mitochondrial PCGs of vent crabs. Alignment and consensus trees were used for the posterior molecular evolutionary analysis (Fig. 1).

 

We used a pair of site models for positive selection at individual codons in vent crabs for each gene, i.e., M8 and M8a (Yang and Nielsen, 2000). The M8a model only allows codons to evolve neutrally or under purifying selection (ω < 1), whereas the M8 model includes a class of sites with ω > 1 (Swanson et al., 2003). We used the free-ratio model (Yang, 1998; Yang and Nielsen, 1998) and branch-site model implemented in CODEML to evaluate whether positive selection was restricted to specific Bythograeidae lineages (Yang, 2007). The positive selection analysis was restricted to the branch of interest (branch leading to the most recently reconstructed ancestor of Bythograeidae, marked a in Figure 1; or the internal branches of Bythograeidae, marked b–i in Fig. 1). The improved branch-site model A (test 2) was performed for every gene in each foreground lineage, which facilitated analysis of the datasets, including all Brachyura species (branch a in Fig. 1) and Bythograeidae only (branches b–i in Fig. 1). The likelihood ratio test (LRT) statistic (2ΔL) approximates a chi-square distribution and was used to compare nested likelihood models. The Bayes empirical Bayes (BEB) approach (Yang et al., 2005) was used to identify amino acids under selection for PAML. We considered candidates sites with posterior probability ≥ 80%.

Positive selected sites were further employed in a series of ML methods implemented in the Datamonkey web server (http://www.datamonkey.org), which has the advantage of improving the dN/dS ratio estimate by incorporating variation in the rate of synonymous substitution (Pond and Frost, 2005; Poon et al., 2009). The fixed-effect likelihood (FEL) calculates site-by-site dN/dS without assuming a prior distribution, and random-effect likelihood (REL) assumes a prior distribution across sites (Poon et al., 2009). Sites with p-values < 0.1 for FEL, and Bayes factor > 50 for REL were considered candidates under positive selection. We then detected selective pressure using the TreeSAAP program (Woolley et al., 2003), which further supports PAML and Datamonkey at the protein physicochemical level. TreeSAAP detected selection based on 31 physicochemical amino acid properties, which were all magnitude category 6–8 changes, and p-values ≤ 0.01 were used as an index for the degree of radical amino acid substitution and positive selection. Finally, to provide further insight into the underlying effects of these positively selected sites, we mapped them onto the protein secondary and three-dimensional structures.

 

Results

Phylogenetic analyses

In this study, the 13 PCGs from the mitochondrial genomes of 48 Brachyura species and one outgroup (C. infraspinatus) were examined. We constructed phylogenetic trees using two different methods (BI and ML) with the concatenated 13 mitochondrial PCGs. The relationships of the gene trees obtained from BI and ML were similar to those previously estimated with morphological and molecular data (Dixon et al., 2003; Shen et al., 2013). We used the consensus tree (Fig. 1) that included all species employed in the present study as the working topology in subsequent analyses.

 

Table I.- Free-ratio (M1 vs M0) analyses of selective pattern on the mitochondrial protein-encoding genes in Brachyura.

Gene	Model	np	ln	2lnL	p value	Parameter estimates
ATP6	M1	193	-14326.21	189.40513	3.03281E-08	ω variation for each branch
	M0	98	-14420.92			ω = 0.02257
ATP8	M1	193	-4967.05	123.717862	0.025502809	ω variation for each branch
	M0	98	-5028.91			0.11106
COX1	M1	193	-25807.19	467.874576	0	ω variation for each branch
	M0	98	-26041.13			0.00895
COX2	M1	193	-13365.38	262.154572	0	ω variation for each branch
	M0	98	-13496.46			0.01912
COX3	M1	193	-14709.41	233.876422	0.0169	ω variation for each branch
	M0	98	-14826.35			0.0169
CYTB	M1	193	-22808.85	539.729226	0	ω variation for each branch
	M0	98	-23078.72			0.01812
ND1	M1	193	-19321.50	299.057304	0	ω variation for each branch
	M0	98	-19471.02			0.01929
ND2	M1	193	-30772.24	361.761826	0	ω variation for each branch
	M0	98	-30953.12			0.03881
ND3	M1	193	-8037.91	186.08815	7.12E-08	ω variation for each branch
	M0	98	-8130.95			0.03067
ND4	M1	193	-31544.00	352.666296	0	ω variation for each branch
	M0	98	-31720.33			0.02634
ND4L	M1	193	-6555.20	145.703226	0.000640494	ω variation for each branch
	M0	98	-6628.05			0.02176
ND5	M1	193	-42402.26	473.378558	0	ω variation for each branch
	M0	98	-42638.95			0.03335
ND6	M1	193	-15377.73	233.091418	1.31E-13	ω variation for each branch
	M0	98	-15494.28			0.03779

 

Molecular evolution of the mitochondrial PCGs in Bythograeidae crabs

We first used the M0 (one-ratio) model that only allows a single ω ratio for all Brachyura branches. The ω values obtained for all 13 PCGs ranged from 0.00895 to 0.11106 and were significantly < 1 (Table I), suggesting that strong purifying selection played a central role in the evolution of Brachyura mitochondrial PCGs to maintain their important functions in OXPHOS. However, the M1 (free-ratio) model of PAML estimated independent ω along all branches of the phylogeny, which was significantly better than the M0 model (p < 0.05, Table I) for all 13 mitochondrial PCGs, suggesting heterogeneous selective pressures on the different lineages.

Several models were compared to further test whether the evidence for positive selection was restricted to the common ancestor and internal branches of Bythograeidae. In the entire Brachyura dataset, the branch-site models were used to test for positive selection in individual codons for the common ancestor of the Bythograeidae (branch a in Fig. 1). Only the ND1 gene was under positive selection in this lineage, and three codons were identified to be under selection by the BEB approach with posterior probabilities ≥ 0.80 (Table II). Next, to evaluate the diversified selection pressures on various lineages within Bythograeidae crabs, branch-site models were used to test for positive selection in individual codons in each lineage across the Bythograeidae phylogeny (branches b–i in Fig. 1). The results show evidence of positive selection in three Bythograeidae-specific lineages; that is, the lineage leading to the last common ancestor of Gandalfus and Austinograea (branch b), the terminal branch of the Segonzacia mesatlantica (branch c) for COX3, and the common ancestor of Gandalfus (branch d) for ND2 and ND5 (Table III; Supplementary Table II).

 

Table II.- Selective pressure analyses (Branch-site model) of mitochondrial PCGs (13 genes) in all crabs dataset and evidence of positive selection of ND1 in the ancestral of Bythograeidae.

Gene	Models	lnL	2lnL	p value	Parameters	Positive selected sites*
ATP6 (branch a)	ma	-14379.6935			ω0 = 0.02 ω1 = 1.0 ω2 = 1.0
	ma0	-14379.6935	0	1.000	ω0 = 0.02 ω1 = 1.0 ω2 = 1.0
ATP8 (branch a)	ma	-4833.41437			ω0 = 0.081 ω1 = 1.0 ω2 = 1.0
	ma0	-4833.41436	-2.8E-05	1.000	ω0 = 0.081 ω1 = 1.0 ω2 = 1.0
COX1 (branch a)	ma	-26017.2764			ω0 = 0.009 ω1 = 1.0 ω2 = 1.0
	ma0	-26017.2764	0	1.000	ω0 = 0.009 ω1 = 1.0 ω2 = 1.0
COX2 (branch a)	ma	-13432.8195			ω0 = 0.018 ω1 = 1.0 ω2 = 1.0
	ma0	-13432.8195	6E-06	0.998	ω0 = 0.018 ω1 = 1.0 ω2 = 1.0
COX3 (branch a)	ma	-14754.578			ω0 = 0.016 ω1 = 1.0 ω2 = 53.483
	ma0	-14754.8316	0.50716	0.476	ω0 = 0.016 ω1 = 1.0 ω2 = 1.0
CYTB (branch a)	ma	-23042.8517			ω0 = 0.018 ω1 = 1.0 ω2 = 1.0
	ma0	-23042.8517	0	1.000	ω0 = 0.018 ω1 = 1.0 ω2 = 1.0
ND1 (branch a)	ma	-19441.9204			ω0 = 0.019 ω1 = 1.0 ω2 = 259.504	251 0.826; 252 0.990; 311 0.992
	ma0	-19445.9633	8.08589	0.004	ω0 = 0.019 ω1 = 1.0 ω2 = 1.0
ND2 (branch a)	ma	-30879.6932			ω0 = 0.039 ω1 = 1.0 ω2 = 19.642
	ma0	-30880.2728	1.15927	0.282	ω0 = 0.039 ω1 = 1.0 ω2 = 1.0
ND3 (branch a)	ma	-8086.63184			ω0 = 0.028 ω1 = 1.0 ω2 = 1.0
	ma0	-8086.63183	-1E-05	1.000	ω0 = 0.028 ω1 = 1.0 ω2 = 1.0
ND4 (branch a)	ma	-31469.8153			ω0 = 0.03 ω1 = 1.0 ω2 = 1.0
	ma0	-31469.8153	0	1.000	ω0 = 0.03 ω1 = 1.0 ω2 = 1.0
ND4L (branch a)	ma	-6625.28512			ω0 = 0.021 ω1 = 1.0 ω2 = 999.0
	ma0	-6625.90273	1.23523	0.266	ω0 = 0.021 ω1 = 1.0 ω2 = 1.0
ND5 (branch a)	ma	-42271.4485			ω0 = 0.034 ω1 = 1.0 ω2 = 8.124
	ma0	-42271.7596	0.6223	0.430	ω0 = 0.034 ω1 = 1.0 ω2 = 1.0
ND6 (branch a)	ma	-15304.8926			ω0 = 0.037 ω1 = 1.0 ω2 = 999.0
	ma0	-15305.5764	1.3676	0.242	ω0 = 0.037 ω1 = 1.0 ω2 = 1.0

*Position of the site in the amino acid sequence (Integer); followed by the posterior probabilities for site (Decimals).

 

Table III.- Evidence of positive selection of mitochondrial PCGs in the interior branchs of Bythograeidae.

Gene	Models	lnL	2lnL	p value	Parameters	Positive selected sites*
COX3	branch b
	ma	-1983.9203			ω0 = 0.006 ω1 = 1.0 ω2 = 52.197	232 0.931
	ma0	-1985.9746	4.108574	0.042666323	ω0 = 0.005 ω1 = 1.0 ω2 = 1.0
	branch c
	ma	-1983.9202			ω0 = 0.006 ω1 = 1.0 ω2 = 50.628	232 0.929
	ma0	-1985.9746	4.10866	0.042664153	ω0 = 0.005 ω1 = 1.0 ω2 = 1.0
ND2	branch d
	ma	-3304.0371			ω0 = 0.034 ω1 = 1.0 ω2 = 9.771	81 0.828; 224 0.989
	ma0	-3306.0675	4.060794	0.043889858	ω0 = 0.034 ω1 = 1.0 ω2 = 1.0
ND5	branch d
	ma	-4972.0567			ω0 = 0.021 ω1 = 1.0 ω2 = 34.114	247 0.932
	ma0	-4974.2993	4.485046	0.034192626	ω0 = 0.021 ω1 = 1.0 ω2 = 1.0

*Position of the site in the amino acid sequence (Integer); followed by the posterior probabilities for site (Decimals).

 

The Bythograeidae-only datasets were also used to test whether specific codons in the mitochondrial PCGs had been subjected to positive selection. The site model (M8 vs. M8a) results detected no significant differences in the estimates using the model incorporating selection (M8) and the neutral model (M8a) for the mitochondrial PCGs (Supplementary Table III). The FEL and REL implemented in Datamonkey which have the advantage of improving the estimate of the ω ratio incorporated variation in the rate of synonymous substitution (Pond and Frost, 2005). A total of 12 sites (FEL: 1 in ND1, 1 in ND4, 2 in ND5; REL: 4 in ATP6, 1 in ATP8, 1 in COX1, 1 in CYTB, 1 in ND4) were identified at seven genes in the Bythograeidae crabs. In addition, 43.75% (14/32) of the putative positively selected sites from seven genes (three in ATP6, one in COX1, three in COX3, one in ND1, two in ND2, one in ND4, and three in ND5) were undergoing radical changes using the protein-level approach performed in TreeSAAP (Supplementary Table IV), which provided additional evidence for positive selection on Bythograeidae crabs.

A positively selected gene is usually more reliable if it can be supported by two or more different methods (Wang et al., 2015). In the present study, COX3 and ND2 were detected using PAML and Datamonkey, whereas five genes (ATP6, ATP8, COX1, CYTB, and ND4) underwent positive selection according to Datamonkey and TreeSAAP. ND1 and ND5 were detected to have undergone strong positive selection according to the three methods (i.e., PAML, Datamonkey, and TreeSAAP). Thus, nine mitochondrial PCGs (ATP6, ATP8, COX1, COX3, CYTB, ND1, ND2, ND4, and ND5) were validated by at least two methods (Supplementary Table IV) and were regarded as positively selected genes.

Structural links to protein function

We mapped all radical amino acid sites onto secondary and tertiary structures to gain insight into the functional significance of the putatively selected sites. Twenty-seven sites (three in ATP6, two in ATP8, one in COX1, three in COX3, 3 in CYTB, three in ND1, three in ND2, two in ND4, and seven in ND5) in nine positively selected genes were situated in the protein binding region (Supplementary Table IV). In addition, a large number of radical changes in amino acids of the positively selected genes were located primarily in the helical transmembrane region (Supplementary Table IV). These results show that most positively selected and radical change sites were located within or close to the functional domain (Supplementary Fig. S1).

 

Discussion

Deep-sea vent crabs (Crustacea, Decapoda, Bythograeidae) are a good model to study physiological adaptation to extreme physical and chemical conditions (Kadar and Powell, 2006). Vent crabs are characterized by specific morphological traits, such as a white body and reduced eyestalks with vestigial corneas to survive in the extreme environment (Min et al., 2017). According to previous studies, organisms living in hydrothermal vents have developed enhanced oxygen extraction capacity and specific ways to process sulfide (Williams, 1980; Sanders et al., 1988; Segonzac et al., 1993; Hourdez and Lallier, 2007). However, the genetic basis of these adaptations to the hydrothermal vent environment remains poorly explored. Therefore, the present study is the first systematic survey of all 13 mitochondrial PCGs in vent crabs to assess signs of positive selection and to determine whether environmental tolerance has influenced the evolution of these genes.

The ω values that we obtained from the one-ratio model were significantly less than 1, providing support for the expected presence of purifying selection acting on all mitochondrial genes to maintain important functions during energy metabolism (Mamirova et al., 2007; Popadin et al., 2007; Stewart et al., 2008). The signal for positive selection is usually swamped by continuous negative selection that occurs on most sites in a gene sequence because positive selection mainly acts on only a few sites and for a short period of evolutionary time (Zhang et al., 2005; Shen et al., 2010). These reasons may partly explain why we did not detect positive selection using the site models.

 

Table IV.- Positive selection in mitochondrial PCGs based on the analysis of the three methods (PAML, DATAMONKEY, or TreeSAAP).

Gene	PAML*	DATAMONKEY*	TreeSAAP*
ATP6	—	+	+
ATP8	—	+	+
COX1	—	+	+
COX3	+	—	+
CYTB	—	+	+
ND1	+	+	+
ND2	+	—	+
ND4	—	+	+
ND5	+	+	+

+, positive selection was detected. —, no positive selection was detected. *results are listed in Supplementary Table IV.

 

Branch-site models are powerful for distinguishing positive selection from the relaxation of purifying selection (Zhang et al., 2005). Thus, four genes (COX3, ND1, ND2, and ND5) showed evidence for significant positive selection according to the LRTs in the branch-site model, and a series of codons were identified as candidate sites that had undergone positive selection (posterior probability ≥ 80%; Tables II, III). According to physiological and morphological researches, Segonzacia, Gandalfus and Austinograea species haemolymph possess higher buffer capacity than littoral species (Tsuchida and Fujikura, 2000; Chausson et al., 2004; Hamasaki et al., 2010). Positive selection of these genes in three Bythograeidae-specific lineages (Table III; Supplementary Table II) might therefore play an important role in promoting their adaptation to hydrothermal vent environment. In addition, adaptive evolution was also supported by evidence that the positively selected sites were detected in Datamonkey and TreeSAAP. Nine genes were determined to have undergone positive selection (Table IV). Particularly, the high proportion of putatively selected codons was localized on or near to the important functional regions (protein binding region and helical transmembrane region) in the mitochondrial protein structure. These nucleotide- and protein-level results indicate extensive adaptive evolution of the mitochondrial PCGs, and suggest that positive selection may be the major driving force for the evolution of mitochondrial PCGs in vent crabs.

Mitochondrial OXPHOS sustains organelle function and plays a central role in cellular energy metabolism (Koopman et al., 2013). The OXPHOS system consists of five multisubunit complexes (Complex I–V), All OXPHOS complexes except complex II, which is exclusively derived from the nuclear DNA, contain subunits that are encoded by both nuclear and mitochondrial genes (Distelmaier et al., 2009; Koopman et al., 2010, 2013). Complex I is the largest OXPHOS enzyme and is central to energy transformation in many prokaryotes and most eukaryotes (Brandt, 2006; Treberg and Brand, 2011). Complex I requires a set of 14 evolutionary conserved “core subunits”, consisting of seven mtDNA-encoded ND subunits and seven nDNA-encoded subunits to oxidize NADH to NAD+ and donates the released electrons to the electron carrier coenzyme Q10 (Koopman et al., 2010; Hirst, 2011). Therefore, the observed positive selection in the ND1, ND2, ND3, and ND5 genes in vent crabs suggests enhanced capability for catalysis using complex I. Complex III contains 11 subunits, one of which is encoded by the mtDNA (CYTB) (Koopman et al., 2013). Complex IV consists of 14 subunits, three of which are mtDNA-encoded (COX1, COX2, and COX3) and catalyze electrons donated to molecular oxygen (O2) to form water (Koopman et al., 2013). About 95% of the O2 organisms breathe is consumed by this complex (Fergusonmiller et al., 2012). The positive selection findings at these genes (CYTB, COX1, and COX3) suggest that vent crabs may have adaptively enhanced oxygen utilization during the OXPHOS process. Complex V (ATP synthase) couples proton flow from the intermembrane space back to the matrix by converting ADP and inorganic phosphate to ATP (Koopman et al., 2013). Complex V is comprised of 19 subunits, two of which are encoded by mtDNA (ATP6 and ATP8). A series of positively selected sites observed in the ATP6 and ATP8 genes was localized in or near the functional regions on the protein structure (Supplementary Table IV), indicating that ATP synthesis was advanced to a certain degree. In summary, these results explain why Bythograeidae crabs can live in hydrothermal vents, and suggest that these crabs might have acquired an enhanced capacity for energy metabolism in an extreme hypoxic environment. Of course, further functional assays are needed to confirm these associations in the future.

Sulfide combines with cytochrome c to suffocate animals to death at a low concentration. Thus, the high concentration of sulfide in the hydrothermal vents is another challenge for organisms. According to previous studies, release of cytochrome c from the mitochondria is a key initial step in the apoptotic process, and cytochrome c oxidase (COX) of the respiratory chain plays a key role in reducing toxicity (Ott et al., 2002; Min and Xu, 2007). Considering the function of COX1 and COX3 in reducing toxicity, positive selection of these genes suggests that vent crabs have evolved an enhanced capacity for inhibiting damage from toxic substances. Similar results were also found in a study in which the COX1 and COX2 mutations helped alvinocaridid shrimp resist sulfide deposits and toxic substances (Wang et al., 2017). Therefore, parallel/convergent studies of these toxic-related genes are necessary in these two species groups in the future.

 

Conclusions

The present study comprehensively investigated 13 mitochondrial PCGs in Bythograeidae crabs. Significant positive selection was examined at nine PCGs, specific for ND1 and ND5, and strong positive selection was detected by the three methods (i.e., PAML, Datamonkey, and TreeSAAP). These results are well matched with the complex adaptation of Bythograeidae crabs to the hydrothermal vent environment. We anticipate that future studies will elucidate the physiology and biochemistry of these PCGs and seek whether a similar evolutionary mechanism underlies animals inhabiting extreme environments.

 

Acknowledgments

This study was funded by the National Natural Science Foundation of China (Grant Number 31702014), and Doctoral Scientific Research Foundation of Yancheng Teachers University to ZFW, and Open Foundation of Jiangsu Key Laboratory for Bioresources of Saline Soils (Grant Number JKLBS2016007).

 

Supplementary material

There is supplementary material associated with this article. Access the material online at: http://dx.doi.org/10.17582/journal.pjz/2019.51.5.1721.1731

Statement of conflict of interest

The authors declare that there is no conflict of interests regarding the publication of this article.

References

Abascal, F., Zardoya, R. and Posada, D., 2005. ProtTest: Selection of best-fit models of protein evolution. Bioinformatics, 21: 2104-2105. https://doi.org/10.1093/bioinformatics/bti263

Bettencourt, R., Pinheiro, M., Egas, C., Gomes, P., Afonso, M., Timothy, S. and Santos, R.S., 2010. High-throughput sequencing and analysis of the gill tissue transcriptome from the deep-sea hydrothermal vent mussel Bathymodiolus azoricus. BMC Genom., 11: 559. https://doi.org/10.1186/1471-2164-11-559

Brandt, U., 2006. Energy converting NADH: Quinone oxidoreductase (complex I). Annu. Rev. Biochem., 75: 69-92. https://doi.org/10.1146/annurev.biochem.75.103004.142539

Carroll, J., Fearnley, I., Wang, Q. and Walker, J., 2009. Measurement of the molecular masses of hydrophilic and hydrophobic subunits of ATP synthase and complex I in a single experiment. Anal. Biochem., 395: 249-255. https://doi.org/10.1016/j.ab.2009.08.006

Chausson, F., Sanglier, S., Leize, E., Hagège, A., Bridges, C.R., Sarradin, P.M., Shillito, B., Lallier, F.H. and Zal, F. 2004. Respiratory adaptations to the deep-sea hydrothermal vent environment: The case of Segonzacia mesatlantica, a crab from the Mid-Atlantic Ridge. Micron, 35: 31-41. https://doi.org/10.1016/j.micron.2003.10.010

Childress, J. and Fisher, C., 1992. The biology of hydrothermal vent animals: Physiology, biochemistry, and autotrophic symbioses. Oceanogr. Mar. Biol., 30: 337-441.

Corliss, J., 1977. Oases of life in the cold abyss. Natl. Geogr. Mag., 152: 441-453.

Darriba, D., Taboada, G., Doallo, R. and Posada, D., 2012 jModelTest 2: More models, new heuristics and parallel computing. Nat. Methods, 9: 772. https://doi.org/10.1038/nmeth.2109

Darriba, D., Taboada, G., Doallo, R. and Posada, D., 2011 ProtTest 3: Fast selection of best-fit models of protein evolution. Bioinformatics, 27: 1164-1165. https://doi.org/10.1093/bioinformatics/btr088

Distelmaier, F., Koopman, W., Rodenburg, J., Mayatepek, E., Willems, P. and Smeitink, J., 2009. Mitochondrial complex I deficiency: from organelle dysfunction to clinical disease. Brain, 132: 833-842. https://doi.org/10.1093/brain/awp058

Dixon, C., Ahyong, S. and Schram, F., 2003. A new hypothesis of decapod phylogeny. Crustaceana, 76: 935-975. https://doi.org/10.1163/156854003771997846

Eubel, H. and Braun, H., 2004. Identification and characterization of respirasomes in potato mitochondria. Pl. Physiol., 134: 1450-1459. https://doi.org/10.1104/pp.103.038018

Fergusonmiller, S., Hiser, C. and Liu, J., 2012. Gating and regulation of the cytochrome c oxidase proton pump. Biochim. biophys. Acta, 1817: 489-494. https://doi.org/10.1016/j.bbabio.2011.11.018

Guo, H., Yang, H., Tao, Y., Tang, D., Wu, Q., Wang, Z. and Tang, B., 2018. Mitochondrial OXPHOS genes provides insights into genetics basis of hypoxia adaptation in anchialine cave shrimps. Genes Genom., 40: 1169-1180. https://doi.org/10.1007/s13258-018-0674-4

Hamasaki, K., Nakajima, K., Tsuchida, S., Kado, R. and Kitada, S., 2010. Number and duration of zoeal stages of the hydrothermal vent crab Gandalfus yunohana from laboratory reared specimens. J. Crust. Biol., 30: 236-240. https://doi.org/10.1651/09-3199.1

Hirst, J., 2011. Why does mitochondrial complex I have so many subunits? Biochem. J., 437: 1-3. https://doi.org/10.1042/BJ20110918

Hourdez, S. and Lallier, F., 2007. Adaptations to hypoxia in hydrothermal-vent and cold-seep invertebrates. Rev. environ. Sci. Biotechnol., 6: 143-159. https://doi.org/10.1007/s11157-006-9110-3

Huelsenbeck, J. and Ronquist, F., 2001. MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics, 17: 754-755. https://doi.org/10.1093/bioinformatics/17.8.754

Jannasch, H. and Mottl, M., 1985. Geomicrobiology of deep-sea hydrothermal vents. Science, 229: 717. https://doi.org/10.1126/science.229.4715.717

Kadar, E. and Powell, J., 2006. Post-capture investigations of hydrothermal vent macro-invertebrates to study adaptations to extreme environments. Rev. environ. Sci. Biotechnol., 5: 193-201. https://doi.org/10.1007/s11157-006-0006-z

Kim, S., Lee, K. and Ju, S., 2013. Nuclear mitochondrial pseudogenes in Austinograea alayseae hydrothermal vent crabs (Crustacea: Bythograeidae): effects on DNA barcoding. Mol. Ecol. Resour., 13: 781-787. https://doi.org/10.1111/1755-0998.12119

Koopman, W., Distelmaier, F., Smeitink, J. and Willems, P., 2013. OXPHOS mutations and neurodegeneration. EMBO J., 32: 9-29. https://doi.org/10.1038/emboj.2012.300

Koopman, W., Nijtmans, L., Dieteren, C., Roestenberg, P., Valsecchi, F., Smeitink, J.A. and Willems, P.H., 2010. Mammalian mitochondrial complex I: Biogenesis, regulation, and reactive oxygen species generation. Antioxid. Redox Signal., 12: 1431-1470. https://doi.org/10.1089/ars.2009.2743

Kulawiec, M., Safina, A., Desouki, M., Still, I., Matsui, S. and Bakin, A. and Singh, K.K., 2008. Tumorigenic transformation of human breast epithelial cells induced by mitochondrial DNA depletion. Cancer Biol. Ther., 7: 1732-1743. https://doi.org/10.4161/cbt.7.11.6729

Kumar, S., Stecher, G. and Tamura, K., 2016. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol., 33: 1870-1874. https://doi.org/10.1093/molbev/msw054

Mamirova, L., Popadin, K. and Gelfand, M., 2007. Purifying selection in mitochondria, free-living and obligate intracellular proteobacteria. BMC Evol. Biol., 7: 1-12. https://doi.org/10.1186/1471-2148-7-1

Martin, W., Baross, J., Kelley, D. and Russell, M., 2008. Hydrothermal vents and the origin of life. Nat. Rev. Microbiol., 6: 805-814. https://doi.org/10.1038/nrmicro1991

Min, H., Song, C., Yuan, L., Li, C. and Cui, Z., 2017. Exploring the molecular basis of adaptive evolution in hydrothermal vent crab Austinograea alayseae by transcriptome analysis. PLoS One, 12: e0178417. https://doi.org/10.1371/journal.pone.0178417

Min, L. and Xu, J., 2007. Detoxifying function of cytochrome c against oxygen toxicity. Mitochondrion, 7: 13-16. https://doi.org/10.1016/j.mito.2006.11.011

Ning, T., Xiao, H., Li, J., Hua, S. and Zhang, Y., 2010 Adaptive evolution of the mitochondrial ND6 gene in the domestic horse. Genet. Mol. Res., 9: 144-150. https://doi.org/10.4238/vol9-1gmr705

Ohta, T., 1992. The nearly neutral theory of molecular evolution. Annu. Rev. Ecol. Syst., 23: 263-286. https://doi.org/10.1146/annurev.es.23.110192.001403

Ott, M., Robertson, J., Gogvadze, V., Zhivotovsky, B. and Orrenius, S., 2002. Cytochrome c release from mitochondria proceeds by a two-step process. Proc. natl. Acad. Sci., 99: 1259-1263. https://doi.org/10.1073/pnas.241655498

Pond, S. and Frost, S., 2005. Datamonkey: Rapid detection of selective pressure on individual sites of codon alignments. Bioinformatics, 21: 2531-2533. https://doi.org/10.1093/bioinformatics/bti320

Poon, A., Frost, D. and Pond, S., 2009. Detecting signatures of selection from DNA sequences using Datamonkey. Bioinform. DNA, 537: 163-183. https://doi.org/10.1007/978-1-59745-251-9_8

Popadin, K., Polishchuk, L., Mamirova, L, Knorre, D. and Gunbin, K., 2007. Accumulation of slightly deleterious mutations in mitochondrial protein-coding genes of large versus small mammals. Proc. natl. Acad. Sci., 104: 13390-13395. https://doi.org/10.1073/pnas.0701256104

Ronquist, F., Teslenko, M., van der Mark, P., Ayres, D., Darling, A., Höhna, S., Larget, B., Liu, L., Suchard, M.A. and Huelsenbeck, J.P., 2012. MrBayes 3.2: Efficient Bayesian phylogenetic inference and model choice across a large model space. Syst. Biol., 61: 539-542. https://doi.org/10.1093/sysbio/sys029

Sanders, N., Arp, A. and Childress, J., 1988. Oxygen binding characteristics of the hemocyanins of two deep-sea hydrothermal vent crustaceans. Respir. Physiol., 71: 57-67. https://doi.org/10.1016/0034-5687(88)90115-6

Santore, M., Mcclintock, D., Lee, V., Budinger, G. and Chandel, N., 2002. Anoxia-induced apoptosis occurs through a mitochondria-dependent pathway in lung epithelial cells. Am. J. Physiol. Lung Cell Mol. Physiol., 282: 727-734. https://doi.org/10.1152/ajplung.00281.2001

Segonzac, M., Desaintlaurent, M. and Casanova, B., 1993. Enigma of the trophic adaptation of the shrimp Alvinocarididae in hydrothermal areas along the Mid-Atlantic Ridge. Cah. Debiol. Marine, 34: 535-571.

Shen, H., Braband, A. and Scholtz, G., 2013. Mitogenomic analysis of decapod crustacean phylogeny corroborates traditional views on their relationships. Mol. Physiol. Evol., 66: 776-789. https://doi.org/10.1016/j.ympev.2012.11.002

Shen, Y., Liang, L., Zhu, Z., Zhou, W., Irwin, D. and Zhang, Y., 2010. Adaptive evolution of energy metabolism genes and the origin of flight in bats. Proc. natl. Acad. Sci., 107: 8666-8671. https://doi.org/10.1073/pnas.0912613107

Shi, Y., Hu, Y., Wang, J., Elzo, M., Yang, X. and Lai, S., 2018. Genetic diversities of MT-ND1 and MT-ND2 genes are associated with high-altitude adaptation in yak. Mitochond. DNA Part A, 29: 485-494.

Solaini, G., Baracca, A., Lenaz, G. and Sgarbi, G., 2010. Hypoxia and mitochondrial oxidative metabolism. Acta Biochim. biophys. Sin., 1797: 1171. https://doi.org/10.1016/j.bbabio.2010.02.011

Stamatakis, A., 2014. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics, 30: 1312-1313. https://doi.org/10.1093/bioinformatics/btu033

Stewart, J., Freyer, C., Elson, J., Wredenberg, A., Cansu, Z., Trifunovic, A. and Larsson, N.G., 2008. Strong purifying selection in transmission of mammalian mitochondrial DNA. PLoS Biol., 6: e10. https://doi.org/10.1371/journal.pbio.0060010

Sun, J., Zhang, Y., Xu, T., Zhang, Y., Mu, H. and Zhang, Y., 2017. Adaptation to deep-sea chemosynthetic environments as revealed by mussel genomes. Nat. Ecol. Evol., 1: 121. https://doi.org/10.1038/s41559-017-0121

Swanson, W., Nielsen, J. and Yang, Q., 2003. Pervasive adaptive evolution in mammalian fertilization proteins. Mol. Evol., 20: 18-20. https://doi.org/10.1093/oxfordjournals.molbev.a004233

Treberg, J. and Brand, M., 2011. A model of the proton translocation mechanism of complex I. J. Biol., 286: 17579-17584. https://doi.org/10.1074/jbc.M111.227751

Tsuchida, S. and Fujikura, K., 2000. Heterochely, relative growth, and gonopod morphology in the bythograeid crab, Austinograea williamsi (Decapoda, Brachyura). J. Crust. Biol., 20: 407-414. https://doi.org/10.1651/0278-0372(2000)020[0407:HRGAGM]2.0.CO;2

Wang, Z., Zhuo, C., Xu, S., Ren, W., Zhou, K. and Yang, G., 2015. Obesity is healthy for cetaceans? Evidence from pervasive positive selection in genes related to triacylglycerol metabolism. Scient. Rep., 5: 14187. https://doi.org/10.1038/srep14187

Wang, Z., Shi, X., Sun, L., Bai, Y., Zhang, D. and Tang, B., 2017. Evolution of mitochondrial energy metabolism genes associated with hydrothermal vent adaption of Alvinocaridid shrimps. Genes Genom., 39: 1-10. https://doi.org/10.1007/s13258-017-0600-1

Williams, A., 1980. A new crab family from the vicinity of submarine thermal vents on the Galapagos rift (Crustacea: Decapoda: Brachyura). Proc. biol. Soc. Wash., 93: 443-472.

Woolley, S., Johnson, J., Smith, M., Crandall, K. and McClellan, D., 2003. TreeSAAP: Selection on amino acid properties using phylogenetic trees. Bioinformatics, 19: 671-672. https://doi.org/10.1093/bioinformatics/btg043

Yang, Y., Xu, S., Xu, J., Guo, Y. and Yang, G., 2014. Adaptive evolution of mitochondrial energy metabolism genes associated with increased energy demand in flying insects. PLoS One, 9: e99120. https://doi.org/10.1371/journal.pone.0099120

Yang, Z., 2007. PAML 4: Phylogenetic analysis by maximum likelihood. Mol. Biol. Evol., 24: 1586-1591. https://doi.org/10.1093/molbev/msm088

Yang, Z., 1998. Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution. Mol. Biol. Evol., 15: 568-573. https://doi.org/10.1093/oxfordjournals.molbev.a025957

Yang, Z., Wong, S. and Nielsen, R., 2005. Bayes empirical bayes inference of amino acid sites under positive selection. Mol. Biol. Evol., 22: 1107-1118. https://doi.org/10.1093/molbev/msi097

Yang, Z. and Nielsen, R., 2000. Estimating synonymous and nonsynonymous substitution rates under realistic evolutionary models. Mol. Biol. Evol., 17: 32-43. https://doi.org/10.1093/oxfordjournals.molbev.a026236

Yang, Z. and Nielsen, R., 1998. Synonymous and nonsynonymous rate variation in nuclear genes of mammals. J. Mol. Evol., 46: 409-418. https://doi.org/10.1007/PL00006320

Zhang, J., Nielsen, R. and Yang, Z., 2005. Evaluation of an improved branch-site likelihood method for detecting positive selection at the molecular level. Mol. Biol. Evol., 22: 2472-2479. https://doi.org/10.1093/molbev/msi237