Evolution of the sulfide-binding function within the globin multigenic family of the deep-sea hydrothermal vent...

The giant extracellular hexagonal bilayer hemoglobin (HBL-Hb) of the deep-sea hydrothermal vent tube worm Riftia pachyptila is able to transport simultaneously O(2) and H(2)S in the blood from the gills to a specific organ: the trophosome that harbors sulfide-oxidizing endosymbionts. This vascular HBL-Hb is made of 144 globins from which four globin types (A1, A2, B1, and B2) coevolve. The H(2)S is bound at a specific location (not on the heme site) onto two of these globin types. In order to understand how such a function emerged and evolved in vestimentiferans and other related annelids, six partial cDNAs corresponding to the six globins known to compose the multigenic family of R. pachyptila have been identified and sequenced. These partial sequences (ca. 120 amino acids, i.e., 80% of the entire protein) were used to reconstruct molecular phylogenies in order to trace duplication events that have led to the family organization of these globins and to locate the position of the free cysteine residues known to bind H(2)S. From these sequences, only two free cysteine residues have been found to occur, at positions Cys + 1 (i.e., 1 a.a. from the well-conserved distal histidine) and Cys + 11 (i.e., 11 a.a. from the same histidine) in globins B2 and A2, respectively. These two positions are well conserved in annelids, vestimentiferans, and pogonophorans, which live in sulfidic environments. The structural comparison of the hydrophobic environment that surrounds these cysteine residues (the sulfide-binding domain) using hydrophobic cluster analysis plots, together with the cysteine positions in paralogous strains, suggests that the sulfide-binding function might have emerged before the annelid radiation in order to detoxify this toxic compound. Moreover, globin evolutionary rates are highly different between paralogous strains. This suggests that either the two globin subfamilies involved in the sulfide-binding function (A2 and B2) have evolved under strong directional selective constraints (negative selection) and that the two other globins (A1 and B1) have accumulated more substitutions through positive selection or have evolved neutrally after a relaxation of selection pressures. A likely scenario on the evolution of this multigenic family is proposed and discussed from this data set.