Deciphering the sexually dimorphic endocrine regulatory network of Pacific white shrimp (Litopenaeus vannamei) after unilateral eyestalk ablation

Abstract

Eyestalk ablation (ESA) is widely used in decapod aquaculture to induce maturation by altering endocrine function. While pronounced sexual dimorphism in growth and development is common in decapods, the endocrine regulatory mechanisms underlying sex-specific responses to ESA remain unclear. Here, we comprehensively investigated sexual dimorphic responses to unilateral ESA in the most important farming decapod—Litopenaeus vannamei—via multi-tissue transcriptomics. Integrated mRNA and sRNA profiling of eyestalk and hepatopancreas tissues revealed pronounced sex-biased differential expression. Only a limited set of differentially expressed genes (DEGs) and transcripts (DETs) were common to both males and females (101 DEGs and 16 DETs in eyestalk; 478 DEGs and 231 DETs in hepatopancreas), with the majority being sex specific. We found unilateral ESA could trigger responses in hormonal pathways and metabolic processes related to carbohydrates, lipids, and proteins, while sex differences in metabolic (lipase and amylase activities), immune responses, and hormone signaling (ecdysteroids and sesquiterpenoids) were revealed. Co-expression network analysis identified sex-specific responsive modules and hub genes (e.g., Trypsin and Chit1). In addition, we found sexually dimorphic post-transcriptional regulation, including widespread alternative splicing, polyadenylation-induced 3′UTR length variation (longer in females and shortened in males), and microRNA–target interactions (e.g., miR-750, miR-2788, and miR-263a/b). Our study systematically reveals sexually dimorphic responses of L. vannamei to unilateral ESA from physiological, gene regulatory network, and post-transcriptional regulation aspects, advancing the understanding of the endocrine and molecular basis of sex dimorphism in decapods and offering insights for shrimp aquaculture improvement.

Keywords

Eyestalk ablation / Sex bias / Hormonal regulation / Alternative splicing / Polyadenylation / MicroRNA

References

[1] Asmat-Ullah M, Rozaimi R, Fazhan H, Shu-Chien AC, Wang Y, Waiho K. Eyestalk ablation to increase ovarian maturation in mud crabs. J vis Exp, 2023, 193e65039

[2] Amankwah BK, Wang C, Zhou T, Liu J, Shi L, Wang W, Chan S-M. Eyestalk ablation, a prerequisite for crustacean reproduction: a review. Isr J Aquacult-Bamid, 2019, 71: 1580

[3] Avram O, Kigel A, Vaisman-Mentesh A, Kligsberg S, Rosenstein S, Dror Y, Pupko T, Wine Y. PASA: proteomic analysis of serum antibodies web server. PLoS Comput Biol, 2021, 17 e1008607

[4] Bai H, Qiao H, Li F, Fu H, Sun S, Zhang W, Jin S, Gong Y, Jiang S, Xiong Y. Molecular characterization and developmental expression of vitellogenin in the oriental river prawn Macrobrachium nipponense and the effects of RNA interference and eyestalk ablation on ovarian maturation. Gene, 2015, 562: 22-31

[5] Betel D, Wilson M, Gabow A, Marks DS, Sander C. The microRNA.org resource: targets and expression. Nucleic Acids Res, 2008, 36: D149-D153

[6] Bicknell AA, Cenik C, Chua HN, Roth FP, Moore MJ. Introns in UTRs: why we should stop ignoring them. BioEssays, 2012, 34: 1025-1034

[7] Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics, 2014, 30: 2114-2120

[8] Campbell DLM, Coman GJ, Lee C, Schmoelzl S, Rombenso AN. The shrimp welfare dilemma: a scientific perspective on eyestalk ablation. Rev Aquacult, 2025, 17 e70082

[9] Cantalapiedra CP, Hernández-Plaza A, Letunic I, Bork P, Huerta-Cepas J. eggNOG-mapper v2: functional annotation, orthology assignments, and domain prediction at the metagenomic scale. Mol Biol Evol, 2021, 38: 5825-5829

[10] Chan S, Wang W, Shi L, Liu J, Zhou T, Wang C, Amankwah BK. Eyestalk ablation, a prerequisite for crustacean reproduction: a review. Isr J Aquacult-Bamid, 2019, 71: 14

[11] Chaves da Silva PGC, Corrêa CL, de Carvalho SL, Allodi S. The crustacean central nervous system in focus: subacute neurodegeneration induces a specific innate immune response. PLoS ONE, 2013, 8 e80896

[12] Chekulaeva M, Hentze MW, Ephrussi A. Bruno acts as a dual repressor of oskar translation, promoting mRNA oligomerization and formation of silencing particles. Cell, 2006, 124: 521-533

[13] Chen H-Y, Toullec J-Y, Lee C-Y. The crustacean hyperglycemic hormone superfamily: progress made in the past decade. Front Endocrinol, 2020, 11 578958

[14] Cho K-H, Daubnerová I, Park Y, Zitnan D, Adams ME. Secretory competence in a gateway endocrine cell conferred by the nuclear receptor βFTZ-F1 enables stage-specific ecdysone responses throughout development in Drosophila. Dev Biol, 2014, 385: 253-262

[15] Ding F, Cui P, Wang Z, Zhang S, Ali S, Xiong L. Genome-wide analysis of alternative splicing of pre-mRNA under salt stress in Arabidopsis. BMC Genom, 2014, 15 431

[16] Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, Batut P, Chaisson M, Gingeras TR. STAR: ultrafast universal RNA-seq aligner. Bioinformatics, 2013, 29: 15-21

[17] Famutimi OG, Adebiyi VG, Akinmolu BG, Dada OV, Adewale IO. Trypsin, chymotrypsin and elastase in health and disease. Future J Pharm Sci, 2024, 10: 126

[18] Fansler MM, Mitschka S, Mayr C. Quantifying 3'UTR length from scRNA-seq data reveals changes independent of gene expression. Nat Commun, 2024, 15: 4050

[19] FAO (2024) The State of World Fisheries and Aquaculture 2024—Blue Transformation in action. Rome

[20] Fingerman M. Crustacean endocrinology: a retrospective, prospective, and introspective analysis. Physiol Zool, 1997, 70: 257-269

[21] Gitterle T, Salte R, Gjerde B, Cock J, Johansen H, Salazar M, Lozano C, Rye M. Genetic (co)variation in resistance to White Spot Syndrome Virus (WSSV) and harvest weight in Penaeus (Litopenaeus) vannamei. Aquaculture, 2005, 246: 139-149

[22] Gruber AR, Martin G, Müller P, Schmidt A, Gruber AJ, Gumienny R, Mittal N, Jayachandran R, Pieters J, Keller W, van Nimwegen E, Zavolan M. Global 3' UTR shortening has a limited effect on protein abundance in proliferating T cells. Nat Commun, 2014, 5: 5465

[23] Guo S, Lin S. mRNA alternative polyadenylation (APA) in regulation of gene expression and diseases. Genes Dis, 2023, 10: 165-174

[24] Guo Z, Zhu L, Cheng Z, Dong L, Guo L, Bai Y, Wu Q, Wang S, Yang X, Xie W, Crickmore N, Zhou X, Lafont R, Zhang Y. A midgut transcriptional regulatory loop favors an insect host to withstand a bacterial pathogen. Innovation, 2024, 5100675

[25] Hari Deva Muthu B, Guru A, Sudhakaran G, Murugan R, Arshad A, Arockiaraj J. Double-edged sword role of shrimp miRNA explains an evolutionary language between shrimp-pathogen interactions that unties the knot of shrimp infection. Rev Aquacult, 2022, 14: 578-593

[26] Huang DW, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res, 2009, 37: 1-13

[27] Huang B, Zhang L, Tang X, Zhang G, Li L. Genome-wide analysis of alternative splicing provides insights into stress adaptation of the Pacific oyster. Mar Biotechnol, 2016, 18: 598-609

[28] Irion U, St Johnston D. Bicoid RNA localization requires specific binding of an endosomal sorting complex. Nature, 2007, 445: 554-558

[29] Jan CH, Friedman RC, Ruby JG, Bartel DP. Formation, regulation and evolution of Caenorhabditis elegans 3'UTRs. Nature, 2011, 469: 97-101

[30] Ji Z, Lee JY, Pan Z, Jiang B, Tian B. Progressive lengthening of 3' untranslated regions of mRNAs by alternative polyadenylation during mouse embryonic development. Proc Natl Acad Sci U S A, 2009, 106: 7028-7033

[31] Keegan LP, Gallo A, O'Connell MA. The many roles of an RNA editor. Nat Rev Genet, 2001, 2: 869-878

[32] Kingan TG, Gray WR, Žitňan D, Adams ME. Regulation of ecdysis-triggering hormone release by eclosion hormone. J Exp Biol, 1997, 200(Pt 24): 3245-3256

[33] Krüger J, Rehmsmeier M. RNAhybrid: microRNA target prediction easy, fast and flexible. Nucleic Acids Res, 2006, 34: W451-454

[34] Lai EC. MicroRNAs are complementary to 3' UTR sequence motifs that mediate negative post-transcriptional regulation. Nat Genet, 2002, 30: 363-374

[35] Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinform, 2008, 9: 1-13

[36] Lee J-H, Suryaningtyas IT, Yoon T-H, Shim JM, Park H, Kim H-W. Transcriptomic analysis of the hepatopancreas induced by eyestalk ablation in shrimp, Litopenaeus vannamei. Comp Biochem Physiol D Genom Proteom, 2017, 24: 99-110

[37] Lee S, Chen YC, Gillen AE, Taliaferro JM, Deplancke B, Li H, Lai EC. Diverse cell-specific patterns of alternative polyadenylation in Drosophila. Nat Commun, 2022, 13: 5372

[38] Li B, Ruotti V, Stewart RM, Thomson JA, Dewey CN. RNA-Seq gene expression estimation with read mapping uncertainty. Bioinformatics, 2010, 26: 493-500

[39] Lianoglou S, Garg V, Yang JL, Leslie CS, Mayr C. Ubiquitously transcribed genes use alternative polyadenylation to achieve tissue-specific expression. Gene Dev, 2013, 27: 2380-2396

[40] Liu X, Xie H, Liu W, Zuo J, Li S, Tian Y, Zhao J, Bai M, Li J, Bao L, Han J, Zhang ZC. Dynamic regulation of alternative polyadenylation by PQBP1 during neurogenesis. Cell Rep, 2024, 43 114525

[41] Luo X, Chen T, Zhong M, Jiang X, Zhang L, Ren C, Hu C. Differential regulation of hepatopancreatic vitellogenin (VTG) gene expression by two putative molt-inhibiting hormones (MIH1/2) in Pacific white shrimp (Litopenaeus vannamei). Peptides, 2015, 68: 58-63

[42] Macdonald PM, Struhl G. Cis-acting sequences responsible for anterior localization of bicoid mRNA in Drosophila embryos. Nature, 1988, 336: 595-598

[43] Maniatis T, Tasic B. Alternative pre-mRNA splicing and proteome expansion in metazoans. Nature, 2002, 418: 236-243

[44] Marini F, Scherzinger D, Danckwardt S. TREND-DB-a transcriptome-wide atlas of the dynamic landscape of alternative polyadenylation. Nucleic Acids Res, 2021, 49: D243-d253

[45] Masamha CP, Xia Z, Yang J, Albrecht TR, Li M, Shyu A-B, Li W, Wagner EJ. CFIm25 links alternative polyadenylation to glioblastoma tumour suppression. Nature, 2014, 510: 412-416

[46] Mastrangelo AM, Marone D, Laidò G, De Leonardis AM, de Vita P. Alternative splicing: enhancing ability to cope with stress via transcriptome plasticity. Plant Sci, 2012, 185–186: 40-49

[47] Matlin AJ, Clark F, Smith CW. Understanding alternative splicing: towards a cellular code. Nat Rev Mol Cell Biol, 2005, 6: 386-398

[48] Mayr C. What are 3' UTRs doing?. CSH Perspect Biol, 2019, 11a034728

[49] McGlincy NJ, Smith CW. Alternative splicing resulting in nonsense-mediated mRNA decay: what is the meaning of nonsense?. Trends Biochem Sci, 2008, 33: 385-393

[50] Meng X, Zhang M, Gao B, Lv J, Li J, Liu P. Integrative proteomic and microRNA analysis: insights into mechanisms of eyestalk ablation-induced ovarian maturation in the swimming crab Portunus trituberculatus. Front Endocrinol, 2020, 11: 533

[51] Miura P, Shenker S, Andreu-Agullo C, Westholm JO, Lai EC. Widespread and extensive lengthening of 3' UTRs in the mammalian brain. Genome Res, 2013, 23: 812-825

[52] Mofayezi A, Jadaliha M, Zangeneh FZ, Khoddami V. Poly(A) tale: from A to A; RNA polyadenylation in prokaryotes and eukaryotes. Wiley Interdiscip Rev RNA, 2024, 15 e1837

[53] Mykles DL. Signaling pathways that regulate the crustacean molting gland. Front Endocrinol, 2021, 12 674711

[54] Nagaraju GPC. Reproductive regulators in decapod crustaceans: an overview. J Exp Biol, 2011, 214: 3-16

[55] Ozsolak F, Kapranov P, Foissac S, Kim SW, Fishilevich E, Monaghan AP, John B, Milos PM. Comprehensive polyadenylation site maps in yeast and human reveal pervasive alternative polyadenylation. Cell, 2010, 143: 1018-1029

[56] Palacios E, Carreño D, Rodrfguez-Jaramillo M, Racotta I. Effect of eyestalk ablation on maturation, larval performance, and biochemistry of white Pacific shrimp, Penaeus vannamei, broodstock. J Appl Aquacult, 1999, 9: 1-23

[57] Pamuru RR. Narayan EJ. Endocrinology of reproduction in crustaceans. Comparative endocrinology of animals, 2019, London, IntechOpen116

[58] Patro R, Duggal G, Love MI, Irizarry RA, Kingsford C. Salmon provides fast and bias-aware quantification of transcript expression. Nat Methods, 2017, 14: 417-419

[59] Perazzolo LM, Gargioni R, Ogliari P, Barracco MAA. Evaluation of some hemato-immunological parameters in the shrimp Farfantepenaeus paulensis submitted to environmental and physiological stress. Aquaculture, 2002, 214: 19-33

[60] Pérez-Rostro CI, Ibarra AM. Heritabilities and genetic correlations of size traits at harvest size in sexually dimorphic Pacific white shrimp (Litopenaeus vannamei) grown in two environments. Aquac Res, 2003, 34: 1079-1085

[61] Pertea G, Pertea M. GFF utilities: GffRead and GffCompare. Res, 2020, 9: 304

[62] Pertea M, Pertea GM, Antonescu CM, Chang TC, Mendell JT, Salzberg SL. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat Biotechnol, 2015, 33: 290-295

[63] Qu Z, Bendena WG, Nong W, Siggens KW, Noriega FG, Kai Z-p, Zang Y-y, Koon AC, Chan HYE, Chan TF, Chu KH, Lam HM, Akam M, Tobe SS, Lam Hui JH. MicroRNAs regulate the sesquiterpenoid hormonal pathway in Drosophila and other arthropods. Proc R Soc B, 2017, 284: 20171827

[64] Qu Z, Bendena WG, Tobe SS, Hui JHL. Juvenile hormone and sesquiterpenoids in arthropods: biosynthesis, signaling, and role of microRNA. J Steroid Biochem, 2018, 184: 69-76

[65] Qu Z, Yiu WC, Yip HY, Nong W, Yu CWC, Lee IHT, Wong AYP, Wong NWY, Cheung FKM, Chan TF, Lau KF, Zhong S, Chu KH, Tobe SS, Ferrier DEK, Bendena WG, Hui JHL. Micro-RNA clusters integrate evolutionary constraints on expression and target affinities: the miR-6/5/4/286/3/309 cluster in Drosophila. Mol Biol Evol, 2020, 37: 2955-2965

[66] Rao X, Huang X, Zhou Z, Lin X. An improvement of the 2ˆ(-delta delta CT) method for quantitative real-time polymerase chain reaction data analysis. Biostat Bioinforma Biomath, 2013, 3: 71-85

[67] Ray D, Kazan H, Cook KB, Weirauch MT, Najafabadi HS, Li X, Gueroussov S, Albu M, Zheng H, Yang A, Na H, Irimia M, Matzat LH, Dale RK, Smith SA, Yarosh CA, Kelly SM, Nabet B, Mecenas D, Li W, et al.. A compendium of RNA-binding motifs for decoding gene regulation. Nature, 2013, 499: 172-177

[68] Rosenzweig R, Nillegoda NB, Mayer MP, Bukau B. The Hsp70 chaperone network. Nat Rev Mol Cell Biol, 2019, 20: 665-680

[69] Sainz Hernández JC, Cordova Murueta JH. Activity of trypsin from Litopenaeus vannamei. Aquaculture, 2009, 290: 190-195

[70] Sainz-Hernández JC, Racotta IS, Dumas S, Hernández-López J. Effect of unilateral and bilateral eyestalk ablation in Litopenaeus vannamei male and female on several metabolic and immunologic variables. Aquaculture, 2008, 283: 188-193

[71] Sandberg R, Neilson JR, Sarma A, Sharp PA, Burge CB. Proliferating cells express mRNAs with shortened 3' untranslated regions and fewer microRNA target sites. Science, 2008, 320: 1643-1647

[72] Sanfilippo P, Wen J, Lai EC. Landscape and evolution of tissue-specific alternative polyadenylation across Drosophila species. Genome Biol, 2017, 18: 229

[73] Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res, 2003, 13: 2498-2504

[74] Sittikankaew K, Pootakham W, Sonthirod C, Sangsrakru D, Yoocha T, Khudet J, Nookaew I, Uawisetwathana U, Rungrassamee W, Karoonuthaisiri N. Transcriptome analyses reveal the synergistic effects of feeding and eyestalk ablation on ovarian maturation in black tiger shrimp. Sci Rep, 2020, 10: 3239

[75] Spies N, Burge CB, Bartel DP. 3' UTR-isoform choice has limited influence on the stability and translational efficiency of most mRNAs in mouse fibroblasts. Genome Res, 2013, 23: 2078-2090

[76] Sun Y, Zhang Y, Liu Y, Xue S, Geng X, Hao T, Sun J. Changes in the organics metabolism in the hepatopancreas induced by eyestalk ablation of the Chinese mitten crab Eriocheir sinensis determined via transcriptome and DGE analysis. PLoS ONE, 2014, 9 e95827

[77] Tang L, Huang F, You W, Poetsch A, Nóbrega RH, Power DM, Zhu T, Liu K, Wang HY, Wang Q, Xu X, Feng B, Schartl M, Shao C. ceRNA crosstalk mediated by ncRNAs is a novel regulatory mechanism in fish sex determination and differentiation. Genome Res, 2022, 32: 1502-1515

[78] Tian J, Wu W, Li J, Wan X, Zhao Z, Xi R, Hu X, Pan M, Xue Y, Yu W. Development dilemma of Litopenaeus vannamei industry in China, current countermeasures taken and its implications for the World Shrimp Aquaculture Industry. Isr J Aquacult-Bamid, 2024, 76: 106-116

[79] Tong R, Pan L, Zhang X, Li Y. Neuroendocrine-immune regulation mechanism in crustaceans: a review. Rev Aquacult, 2022, 14: 378-398

[80] Treerattrakool S, Udomkit A, Eurwilaichitr L, Sonthayanon B, Panyim S. Expression of biologically active crustacean hyperglycemic hormone (CHH) of Penaeus monodon in Pichia pastoris. Mar Biotechnol, 2003, 5: 373-379

[81] Truman JW, Riddiford LM. Endocrine insights into the evolution of metamorphosis in insects. Annu Rev Entomol, 2002, 47: 467-500

[82] Uawisetwathana U, Leelatanawit R, Klanchui A, Prommoon J, Klinbunga S, Karoonuthaisiri N. Insights into eyestalk ablation mechanism to induce ovarian maturation in the black tiger shrimp. PLoS ONE, 2011, 6 e24427

[83] Venkat S, Tisdale AA, Schwarz JR, Alahmari AA, Maurer HC, Olive KP, Eng KH, Feigin ME. Alternative polyadenylation drives oncogenic gene expression in pancreatic ductal adenocarcinoma. Genome Res, 2020, 30: 347-360

[84] Waiho K, Ikhwanuddin M, Afiqah-Aleng N, Shu-Chien AC, Wang Y, Ma H, Fazhan H. Transcriptomics in advancing portunid aquaculture: a systematic review. Rev Aquacult, 2022, 14: 2064-2088

[85] Wainwright G, Webster SG, Wilkinson MC, Chung JS, Rees HH. Structure and significance of mandibular organ-inhibiting hormone in the crab, Cancer pagurus. Involvement in multihormonal regulation of growth and reproduction. J Biol Chem, 1996, 271: 12749-12754

[86] Wang R, Zheng D, Yehia G, Tian B. A compendium of conserved cleavage and polyadenylation events in mammalian genes. Genome Res, 2018, 28: 1427-1441

[87] Wang Y, Xie Z, Kutschera E, Adams JI, Kadash-Edmondson KE, Xing Y. rMATS-turbo: an efficient and flexible computational tool for alternative splicing analysis of large-scale RNA-seq data. Nat Protoc, 2024, 19: 1083-1104

[88] Webster SG, Keller R. Purification, characterisation and amino acid composition of the putative moult-inhibiting hormone (MIH) of Carcinus maenas (Crustacea, Decapoda). J Comp Physiol B, 1986, 156: 617-624

[89] Wu YS, Chang CH, Nan FH. Steroid hormone “cortisone” and “20-hydroxyecdysone” involved in the non-specific immune responses of white shrimp (Litopenaeus vannamei). Fish Shellfish Immunol, 2016, 56: 272-277

[90] Yang Z, Yang X, Du J, Wei C, Liu P, Hu J, Bao Z, Qu Z. Comparative transcriptome analysis of hepatopancreas reveals sexual dimorphic response to methyl farnesoate injection in Litopenaeus vannamei. Int J Mol Sci, 2024, 25: 8152

[91] Yao L, Wang H, Song Y, Sui G. BioQueue: a novel pipeline framework to accelerate bioinformatics analysis. Bioinformatics, 2017, 33: 3286-3288

[92] Zhang L, Pan L, Xu L, Si L. Effects of ammonia-N exposure on the concentrations of neurotransmitters, hemocyte intracellular signaling pathways and immune responses in white shrimp Litopenaeus vannamei. Fish Shellfish Immunol, 2018, 75: 48-57

[93] Zhang X, Yuan J, Sun Y, Li S, Gao Y, Yu Y, Liu C, Wang Q, Lv X, Zhang X, Ma KY, Wang X, Lin W, Wang L, Zhu X, Zhang C, Zhang J, Jin S, Yu K, Kong J, et al.. Penaeid shrimp genome provides insights into benthic adaptation and frequent molting. Nat Commun, 2019, 10: 356

[94] Zhang X, Zhang X, Yuan J, Li F. The responses of alternative splicing during heat stress in the Pacific white shrimp Litopenaeus vannamei. Genes, 2023, 14: 1473

[95] Zhang Y, Huang Z, Lu W, Liu Z. Alternative cleavage and polyadenylation: key regulatory mechanisms in health and disease. RNA Biol, 2025, 22: 1-33

[96] Zheng D, Wang R, Ding Q, Wang T, Xie B, Wei L, Zhong Z, Tian B. Cellular stress alters 3′UTR landscape through alternative polyadenylation and isoform-specific degradation. Nat Commun, 2018, 9: 2268

[97] Žitňan D, Kim YJ, Žitňanová I, Roller L, Adams ME. Complex steroid–peptide–receptor cascade controls insect ecdysis. Gen Comp Endocrinol, 2007, 153: 88-96

[98] Zmora N, Trant J, Zohar Y, Chung JS. Molt-inhibiting hormone stimulates vitellogenesis at advanced ovarian developmental stages in the female blue crab, Callinectes sapidus 1: an ovarian stage dependent involvement. Saline Syst, 2009, 5 7

RIGHTS & PERMISSIONS

The Author(s)

Just Accepted

This article has successfully passed peer review and final editorial review, and will soon enter typesetting, proofreading and other publishing processes. The currently displayed version is the accepted final manuscript. The officially published version will be updated with format, DOI and citation information upon launch. We recommend that you pay attention to subsequent journal notifications and preferentially cite the officially published version. Thank you for your support and cooperation.