Abstract
Mitochondrial protein-coding genes (mt genes) encode subunits forming complexes of crucial cellular pathways, including those involved in the vital process of oxidative phosphorylation (OXPHOS). Despite the vital role of the mitochondrial genome (mt genome) in the survival of organisms, little is known with respect to its adaptive implications within marine invertebrates. The molluscan Class Cephalopoda is represented by a marine group of species known to occupy contrasting environments ranging from the intertidal to the deep sea, having distinct metabolic requirements, varied body shapes and highly advanced visual and nervous systems that make them highly competitive and successful worldwide predators. Thus, cephalopods are valuable models for testing natural selection acting on their mitochondrial subunits (mt subunits). Here, we used concatenated mt genes from 17 fully sequenced mt genomes of diverse cephalopod species to generate a robust mitochondrial phylogeny for the Class Cephalopoda. We followed an integrative approach considering several branches of interest–covering cephalopods with distinct morphologies, metabolic rates and habitats–to identify sites under positive selection and localize them in the respective protein alignment and/or tridimensional structure of the mt subunits. Our results revealed significant adaptive variation in several mt subunits involved in the energy production pathway of cephalopods: ND5 and ND6 from Complex I, CYTB from Complex III, COX2 and COX3 from Complex IV, and in ATP8 from Complex V. Furthermore, we identified relevant sites involved in protein-interactions, lining proton translocation channels, as well as disease/deficiencies related sites in the aforementioned complexes. A particular case, revealed by this study, is the involvement of some positively selected sites, found in Octopoda lineage in lining proton translocation channels (site 74 from ND5) and in interactions between subunits (site 507 from ND5) of Complex I.
Highlights
The mitochondrial subunits of the respiratory chain complexes–cytochrome c oxidase subunit 1–3 (COX1, COX2, COX3), cytochrome b (CYTB), NADH dehydrogenase subunit 1–6 (ND1, ND2, ND3, ND4, ND5, ND6), NADH dehydrogenase subunit 4L (ND4L), ATPase F0 subunit 6 (ATP6) and ATPase F0 subunit 8 (ATP8)–encoded by 13 mitochondrial genes in most metazoans, are involved in several key evolutionary processes of eukaryotes having a major role in the production of energy [1]
Since some Cephalopoda mt genomes have duplicated protein-coding genes (S2 Table)–a feature only found in Oegopsida and Bathyteuthoidea species (S2 Table) [24, 25]: with cox1, cox2, atp6 and atp8 genes duplicated in Dosidicus gigas (NC_009734), Sthenoteuthis oualaniensis (EU660577), Todarodes pacificus (NC_006354), Watasenia scintillans (NC_007893), Architeuthis dux (FJ429092) and Bathyteuthis abyssicola (AP012225) and cox3 that is duplicated in all the above species except for Bathyteuthis abyssicola (S2 Table)–the codon based CDS alignment of the aforementioned genes allowed the verification that the duplicated genes are highly conserved; some are exact copies and others present slight nucleotide substitutions (S2 Table)
To complement our analyses with the spatial position of these sites in a tridimensional plane, and since the X-ray crystal structures of the mitochondrial encoded proteins from cephalopods are not available on the Protein Data Bank (PDB) [59], we modeled its 3D structures for the Octopus vulgaris species, from ND5 (YP_112440.1), ND6 (YP_112444.1), CYTB (YP_112443.1), COX2 (YP_112437.1), COX3 (YP_112433.1) and ATP8 (YP_112438.1), using the I-TASSER server with default parameters [60]
Summary
The mitochondrial subunits of the respiratory chain complexes (mt subunits)–cytochrome c oxidase subunit 1–3 (COX1, COX2, COX3), cytochrome b (CYTB), NADH dehydrogenase subunit 1–6 (ND1, ND2, ND3, ND4, ND5, ND6), NADH dehydrogenase subunit 4L (ND4L), ATPase F0 subunit 6 (ATP6) and ATPase F0 subunit 8 (ATP8)–encoded by 13 mitochondrial genes (mt genes) in most metazoans, are involved in several key evolutionary processes of eukaryotes having a major role in the production of energy [1]. The mt subunits interact with nuclear-encoded proteins establishing four of the five multi-subunit enzyme complexes–Complex I (ND1, ND2, ND3, ND4, ND4L, ND5 and ND6), III (CYTB), IV (COX1, COX2 and COX3) and V (ATP6 and ATP8)–involved in the respiratory chain of aerobic cells. The effect of any amino acid substitution can vary greatly depending on its location in the protein structure. The amino acid substitutions can improve or decrease aerobic capacity and be linked with life history traits and environmental adaptation [6,7,8]
Published Version (Free)
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have