Identification and Positive Selection Analysis of Metallothionein fromDeep Sea Clams in Vesicomyidae
KONG Xue1,2,3, CHEN Sunan1, ZHU Xinhui1, WU Ke1, CAI Yuefeng1,2, ZHI Ying1, LI Yanan4, ZHANG Haibin5, SHEN Xin1,2
1. School of Marine Science and Fisheries, Jiangsu Ocean University, Lianyungang 222000, China; 2. Jiangsu Key Laboratory of Marine Bioresources and Environment/Jiangsu Key Laboratory of MarineBiotechnology, Jiangsu Ocean University, Lianyungang 222000, China; 3. Key Laboratory of MarineBiogenetic Resources, Ministry of Natural Resources, Xiamen 361000, China; 4. College of Animal Science andTehnology, Zhongkai University of Agriculture and Engineering, Guangzhou 510225, China; 5. Institute ofDeep-sea Science and Engineering, Chinese Academy of Sciences, Sanya 572000, China
Abstract:Deep-sea cold seep and hydrothermal ecosystems are characterized by extreme physical and chemical environments, including high hydrostatic pressure, darkness, high concentrations of heavy metals and nutrient limitation. As the dominant taxa in this habitat, clams in Vesicomyidae provide an ideal model for analyzing the adaptive evolution of deep-sea organisms. In this study, we conducted phylogenetic analysis and positive selection pressure detection of metallothioneins (MTs) in deep/shallow-sea bivalves with the aim of elucidating the molecular basis of their adaptation to the deep-sea environment. Five MTs coding sequences with molecular weight of 6—7 ku and isoelectric point of 7.2—8.5 were identified from the omics data of four species in Vesicomyidae. Multiple sequence comparisons revealed that MT-10 isoforms contained 19 conserved cysteine residues in both deep- and shallow-sea species, forming nine CXC conserved motifs, except for CmagMT-10Ⅱ. The phylogenetic evolutionary tree showed that the species (AmMT-10, ApMT-10, CmagMT-10Ⅰ, PoMT-10) in Vesicomyidae were clustered into a single branch and clustered into a bigger evolutionary branch with shallow-water clams MT-10s. Then, a branch-site model for selection pressure analysis, with the deep-sea clam MT-10 as the foreground branch and the shallow sea shellfish MT-10 as the background branch, showed that ω (non-synonymous substitution rate/synonymous substitution rate) was 30.705 2 (P<0.05). The assay identified two strong selection sites (30E, 46T) with amino acid substitution patterns (30P/S/A → E; 46V/A → T) leading to a shift in charge properties. Protein stability prediction by DUET and DDGun showed that these substitutions resulted in a negative shift in ΔΔG values (maximum -1.1 kcal/mol), suggesting that increased structural flexibility may enhance substrate/ligand binding efficiency.It is hypothesized that under the selective pressure of the deep-sea environment, the clam MT-10s in Vesicomyidae achieves functional adaptive evolution through charge modification at key sites. The finding provides new molecular evidence for understanding of the extreme environmental adaptation of deep-sea organisms.
[1] BLINDAUER C A, LESZCZYSZYN O I. Metallothioneins:unparalleled diversity in structures and functions for metal ion homeostasis and more[J]. Natural Product Reports,2010,27(5):720-741. [2] BINZ P A, KÄGI J H R. Metallothionein:Molecular Evolution and Classification[M]//KLAASSEN C D. Metallothionein IV. Basel:Birkhäuser Basel,1999:7-13. [3] CHIAVERINI N, DE LEY M. Protective effect of metallothionein on oxidative stress-induced DNA damage[J]. Free Radical Research,2010,44(6):605-613. [4] DUTSCH-WICHEREK M, SIKORA J, TOMASZEWSKA R. The possible biological role of metallothionein in apoptosis[J]. Frontiers in Bioscience,2008,13:4029-4038. [5] CAPDEVILA M, BOFILL R, PALACIOS Ò, et al. State-of-the-art of metallothioneins at the beginning of the 21st century[J]. Coordination Chemistry Reviews,2012,256(1/2):46-62. [6] 赵汉取,施沁璇,沈萍萍,等.低浓度Cd2+胁迫对青鱼组织SOD活性和MT诱导的影响[J].水生态学杂志,2014,35(2):90-94. [7] FELIX-PORTILLO M, MARTINEZ-QUINTANA J A, PEREGRINO-URIARTE A B, et al. The metallothionein gene from the white shrimp Litopenaeus vannamei:characterization and expression in response to hypoxia[J]. Marine Environmental Research,2014,101:91-100. [8] ZHANG J Q, WANG J, GUI T S, et al. A copper-induced metallothionein gene from Exopalaemon carinicauda and its response to heavy metal ions[J]. International Journal of Biological Macromolecules,2014,70:246-250. [9] REN F, JIANG H, SUN J L, et al. Cloning, characterization, expression, and copper sensitivity of the metallothionein-1 gene in the Chinese mitten crab, Eriocheir sinensis[J]. Molecular Biology Reports,2011,38(4):2383-2393. [10] LI L, SHEN Y C, LIANG J R, et al. Accumulation and depuration of Cd and its effect on the expressions of metallothionein and apoptotic genes in Litopenaeus vannamei[J]. Bulletin of Environmental Contamination and Toxicology,2021,106(3):501-506. [11] MACKAY E A, OVERNELL J, DUNBAR B, et al. Complete amino acid sequences of five dimeric and four monomeric forms of metallothionein from the edible mussel Mytilus edulis[J]. European Journal of Biochemistry,1993,218(1):183-194. [12] PARK J S, CHUNG S, PARK I S, et al. Purification and characterization of metallothionein-like cadmium binding protein from Asian periwinkle Littorina brevicula[J]. Comparative Biochemistry and Physiology Part C:Toxicology & Pharmacology,2002,131(4):425-431. [13] LEE S Y, NAM Y K. Transcriptional responses of metallothionein gene to different stress factors in Pacific abalone (Haliotis discushannai)[J]. Fish & Shellfish Immunology,2016,58:530-541. [14] 周湖明.近江牡蛎富集和排出Cd、Pb及其与金属硫蛋白含量相关性的研究[D].湛江:广东海洋大学,2013. [15] ADKINS J F, MCINTYRE K, SCHRAG D P. The salinity, temperature, and δ18O of the glacial deep ocean[J]. Science,2002,298(5599):1769-1773. [16] MARTÍNEZ-MARTÍNEZ I, NAVARRO-FERNÁNDEZ J, GARCÍA-CARMONA F, et al. Characterization and structural modeling of a novel thermostable glycine oxidase from Geobacillus kaustophilus HTA426[J]. Proteins:Structure, Function, and Bioinformatics,2008,70(4):1429-1441. [17] FENG D, QIU J W, HU Y, et al. Cold seep systems in the South China Sea:an overview[J]. Journal of Asian Earth Sciences,2018,168:3-16. [18] VRIJENHOEK R C. Genetic diversity and connectivity of deep-sea hydrothermal vent metapopulations[J]. Molecular Ecology, 2010,19(20):4391-4411. [19] SISCAR R,KOENIG S, TORREBLANCA A, et al. The role of metallothionein and selenium in metal detoxification in the liver of deep-sea fish from the NW Mediterranean Sea[J]. The Science of the Total Environment,2014,466/467:898-905. [20] GONZALEZ-REY M, SERAFIM A, COMPANY R, et al. Detoxification mechanisms in shrimp:comparative approach between hydrothermal vent fields and estuarine environments[J]. Marine Environmental Research,2008,66(1):35-37. [21] WONG Y H, SUN J, HE L S, et al. High-throughput transcriptome sequencing of the cold seep mussel Bathymodiolus platifrons[J]. Scientific Reports,2015,5:16597. [22] MA L, WANG W X. Subcellular metal distribution in two deep-sea mollusks:insight of metal adaptation and detoxification near hydrothermal vents[J]. Environmental Pollution,2020,266:115303. [23] 李俊杰,于朋,李昀,等.花鲈Claudins基因家族的鉴定、进化分析及对环境盐度的表达响应[J].海洋科学,2022,46(8):65-78. [24] 周慧,吴海霞,谢伟,等.温度和光照对克氏原螯虾生物钟基因PcClk与PcCry1表达的影响[J].广西科学,2022,29(2):393-404. [25] 吴海霞,时昕晔,周开亚,等.中华绒螯蟹翻译起始因子4G基因在幼体发育阶段的mRNA表达模式与适应性进化分析[J].南京师大学报(自然科学版),2022,45(3):96-107. [26] OLU K, SIBUET M, HARMEGNIES F, et al. Spatial distribution of diverse cold seep communities living on various diapiric structures of the southern Barbados prism[J]. Progress in Oceanography,1996,38(4):347-376. [27] JOHNSON S B, KRYLOVA E M, AUDZIJONYTE A, et al. Phylogeny and origins of chemosynthetic vesicomyid clams[J]. Systematics and Biodiversity,2017,15(4):346-360. [28] HONGO Y, IKUTA T, TAKAKI Y, et al. Expression of genes involved in the uptake of inorganic carbon in the gill of a deep-sea vesicomyid clam harboring intracellular thioautotrophic bacteria[J]. Gene,2016,585(2):228-240. [29] LAN Y, SUN J, ZHANG W P, et al. Host-symbiont interactions in deep-sea chemosymbiotic vesicomyid clams:insights from transcriptome sequencing[J]. Frontiers in Marine Science,2019,6:680. [30] IP J C, XU T, SUN J, et al. Host-endosymbiont genome integration in a deep-sea chemosymbiotic clam[J]. Molecular Biology and Evolution,2021,38(2):502-518. [31] COMBOSCH D J, COLLINS T M, GLOVER E A, et al. A family-level tree of life for bivalves based on a Sanger-sequencing approach[J]. Molecular Phylogenetics and Evolution,2017,107:191-208. [32] VITHANA M V P, XU M, ZHAO X, et al. Geological and geophysical signatures of the east Pacific Rise 8°~10°N[J]. Solid Earth Sciences,2019,4(2):66-83. [33] LETUNIC I, BORK P. 20 years of the SMART protein domain annotation resource[J]. Nucleic Acids Research,2018,46(D1):D493-D496. [34] LARKIN M A, BLACKSHIELDS G, BROWN N P, et al. Clustal W and clustal X version 2.0[J]. Bioinformatics,2007,23(21):2947-2948. [35] TAMURA K, PETERSON D, PETERSON N, et al. MEGA5:molecular evolutionary genetics analysis using maximum likelihood,evolutionary distance, and maximum parsimony methods[J]. Molecular Biology and Evolution,2011,28(10):2731-2739. [36] GAO F L, CHEN C J, ARAB D A, et al. EasyCodeML:a visual tool for analysis of selection using CodeML[J]. Ecology and Evolution,2019,9(7):3891-3898. [37] PIRES D E V, ASCHER D B, BLUNDELL T L. DUET:a server for predicting effects of mutations on protein stability using an integrated computational approach[J]. Nucleic Acids Research,2014,42(Web Server issue):W314-W319. [38] MONTANUCCI L, CAPRIOTTI E, FRANK Y, et al. DDGun:an untrained method for the prediction of protein stability changes upon single and multiple point variations[J]. BMC Bioinformatics,2019,20(Suppl 14):335. [39] 王诚远.池蝶蚌金属硫蛋白对重金属诱导响应及其转录调控机制[D].南昌:南昌大学,2018. [40] SOAZIG L, MARC L. Potential use of the levels of the mRNA of a specific metallothionein isoform (MT-20) in mussel (Mytilus edulis) as a biomarker of cadmium contamination[J]. Marine Pollution Bulletin,2003,46(11):1450-1455. [41] HARDIVILLIER Y, DENIS F, DEMATTEI M V, et al. Metal influence on metallothionein synthesis in the hydrothermal vent mussel Bathymodiolus thermophilus[J]. Comparative Biochemistry and Physiology Part C:Toxicology & Pharmacology,2006,143(3):321-332. [42] DI GIULIO M. A comparison of proteins from Pyrococcus furiosus and Pyrococcus abyssi:barophily in the physicochemical properties of amino acids and in the genetic code[J]. Gene,2005,346:1-6. [43] REEVES E P, SEEWALD J S, SACCOCIA P, et al. Geochemistry of hydrothermal fluids from the PACMANUS, Northeast Pual and Vienna Woods hydrothermal fields, Manus Basin, Papua New Guinea[J]. Geochimica et Cosmochimica Acta,2011,75(4):1088-1123. [44] DEMINA L L, GALKIN S V. On the role of abiogenic factors in the bioaccumulation of heavy metals by the hydrothermal fauna of the Mid-Atlantic Ridge[J]. Oceanology,2008,48(6):784-797. [45] SCHMIDT C, BEYERSMANN D. Transient peaks in zinc and metallothionein levels during differentiation of 3T3L1 cells[J]. Archives of Biochemistry and Biophysics,1999,364(1):91-98. [46] 吕新芳,毛伟腾,滑朝阳,等.海洋无脊椎动物金属硫蛋白研究进展[J].海洋通报,2015,34(3):241-246. [47] GERET F, ROUSSE N, RISO R, et al. Metal compartmentalization and metallothionein isoforms in mussels from the Mid-Atlantic Ridge;preliminary approach to the fluid-organism relationship[J]. Cahiers De Biologie Marine,1998,39:291-293. [48] WANG Q, ZHANG L B, YANG D L, et al. Molecular diversity and evolution of defensins in the Manila clam Ruditapes philippinarum[J]. Fish & Shellfish Immunology,2015,47(1):302-312. [49] TONG C, LIN Y Q, ZHANG C F, et al. Transcriptome-wide identification, molecular evolution and expression analysis of Toll-like receptor family in a Tibet fish, Gymnocypris przewalskii[J]. Fish & Shellfish Immunology,2015,46(2):334-345. [50] LAIMER J, HOFER H, FRITZ M, et al. MAESTRO—multi agent stability prediction upon point mutations[J]. BMC Bioinformatics,2015,16:116-129. [51] LENGYEL C S E, WILLIS L J, MANN P, et al. Mutations designed to destabilize the receptor-bound conformation increase MICA-NKG2D association rate and affinity[J]. Journal of Biological Chemistry,2007,282(42):30658-30666. [52] SHERRY D, PANDIAN R, SAYED Y. Non-active site mutations in the HIV protease:diminished drug binding affinity is achieved through modulating the hydrophobic sliding mechanism[J]. International Journal of Biological Macromolecules,2022,217:27-41. [53] LI Y N, KONG X, CHEN J W, et al. Characteristics of the copper, zinc superoxide dismutase of a hadal sea cucumber (Paelopatides sp. ) from the Mariana trench[J]. Marine Drugs,2018,16(5):169-184. [54] WANG S L, LI Y L, HAN Z, et al. Molecular characterization of a novel N-acetylneuraminate lyase from a deep-sea symbiotic Mycoplasma[J]. Marine Drugs,2018,16(3):80-93. [55] LI Y N, CHEN Z F, ZHANG H B, et al. Analysis of environmental selection pressure of superoxide dismutase in deep-sea sea cucumber[J]. Journal of Oceanology and Limnology,2024,42(3):893-904.