Survey Sequencing and Comparative Analysis of the Elephant Shark (Callorhinchus milii) Genome

Owing to their phylogenetic position, cartilaginous fishes (sharks, rays, skates, and chimaeras) provide a critical reference for our understanding of vertebrate genome evolution. The relatively small genome of the elephant shark, Callorhinchus milii, a chimaera, makes it an attractive model cartilaginous fish genome for whole-genome sequencing and comparative analysis. Here, the authors describe survey sequencing (1.4× coverage) and comparative analysis of the elephant shark genome, one of the first cartilaginous fish genomes to be sequenced to this depth. Repetitive sequences, represented mainly by a novel family of short interspersed element–like and long interspersed element–like sequences, account for about 28% of the elephant shark genome. Fragments of approximately 15,000 elephant shark genes reveal specific examples of genes that have been lost differentially during the evolution of tetrapod and teleost fish lineages. Interestingly, the degree of conserved synteny and conserved sequences between the human and elephant shark genomes are higher than that between human and teleost fish genomes. Elephant shark contains putative four Hox clusters indicating that, unlike teleost fish genomes, the elephant shark genome has not experienced an additional whole-genome duplication. These findings underscore the importance of the elephant shark as a critical reference vertebrate genome for comparative analysis of the human and other vertebrate genomes. This study also demonstrates that a survey-sequencing approach can be applied productively for comparative analysis of distantly related vertebrate genomes.     

Survey sequencing and comparative analysis of the elephant shark (Callorhinchus milii) genome

Owing to their phylogenetic position, cartilaginous fishes (sharks, rays, skates, and chimaeras) provide a critical reference for our understanding of vertebrate genome evolution. The relatively small genome of the elephant shark, Callorhinchus milii, a chimaera, makes it an attractive model cartilaginous fish genome for whole-genome sequencing and comparative analysis. Here, the authors describe survey sequencing (1.4× coverage) and comparative analysis of the elephant shark genome, one of the first cartilaginous fish genomes to be sequenced to this depth. Repetitive sequences, represented mainly by a novel family of short interspersed element–like and long interspersed element–like sequences, account for about 28% of the elephant shark genome. Fragments of approximately 15,000 elephant shark genes reveal specific examples of genes that have been lost differentially during the evolution of tetrapod and teleost fish lineages. Interestingly, the degree of conserved synteny and conserved sequences between the human and elephant shark genomes are higher than that between human and teleost fish genomes. Elephant shark contains putative four Hox clusters indicating that, unlike teleost fish genomes, the elephant shark genome has not experienced an additional whole-genome duplication. These findings underscore the importance of the elephant shark as a critical reference vertebrate genome for comparative analysis of the human and other vertebrate genomes. This study also demonstrates that a survey-sequencing approach can be applied productively for comparative analysis of distantly related vertebrate genomes.     

Allosteric water and phosphate effects in Hoplosternum littorale hemoglobins.

This paper reports the results obtained using the osmotic stress method applied to the purified cathodic and anodic hemoglobins (Hbs) from the catfish Hoplosternum littorale, a species that displays facultative accessorial air oxygenation. We demonstrate that water potential affects the oxygen affinity of H. littorale Hbs in the presence of an inert solute (sucrose). Oxygen affinity increases when water activity increases, indicating that water molecules stabilize the high-affinity state of the Hb. This effect is the same as that observed in tetrameric vertebrate Hbs. We show that both anodic and cathodic Hbs show conformational substrates similar to other vertebrate Hbs. For both Hbs, addition of anionic effectors, especially chloride, strongly increases the number of water molecules bound, although anodic Hb did not exhibit sensitivity to saturating levels of ATP. Accordingly, for both Hbs, we propose that the deoxy conformations coexist in at least two anion-dependent allosteric states, T(o) and T(x), as occurs for human Hb. We found a single phosphate binding site for the cathodic Hb.     

Crystal structures of deoxy- and carbonmonoxyhemoglobin F1 from the hagfish Eptatretus burgeri

Hagfish are extremely primitive jawless fish of disputed ancestry. Although generally classed with lampreys as cyclostomes ("round mouths"), it is clear that they diverged from them several hundred million years ago. The crystal structures of the deoxy and CO forms of hemoglobin from a hagfish (Eptatretus burgeri) have been solved at 1.6 and 2.1 A, respectively. The deoxy crystal contains one dimer and two monomers in a unit cell, with the dimer being similar to that found in lamprey deoxy-Hb, but with a larger interface and different relative orientation of the partner chains. Ile(E11) and Gln(E7) obstruct ligand binding in the deoxy form and make room for ligands in the CO form, but no interaction path between the two hemes could be identified. The BGH core structure, which forms the alpha1beta1 interface of all vertebrate alpha2beta2 tetrameric Hbs, is conserved in hagfish and lamprey Hbs. It was shown previously that human and cartilaginous fish Hbs have independently evolved stereochemical mechanisms other than the movement of the proximal histidine to regulate ligand binding at the hemes. Our results therefore suggest that the formation of the alpha2beta2 tetramer using the BGH core and the mechanism of quaternary structure change evolved between the branching points of hagfish and lampreys from other vertebrates.     

Molecular cloning of C4 gene and identification of the class III complement region in the shark MHC

To clarify the evolutionary origin of the linkage of the MHC class III complement genes with the MHC class I and II genes, we isolated C4 cDNA from the banded hound shark (Triakis scyllium). Upon phylogenetic tree analysis, shark C4 formed a well-supported cluster with C4 of higher vertebrates, indicating that the C3/C4 gene duplication predated the divergence of cartilaginous fish from the main line of vertebrate evolution. The deduced amino acid sequence predicted the typical C4 three-subunits chain structure, but without the histidine residue catalytic for the thioester bond, suggesting the human C4A-like specificity. The linkage analysis of the complement genes, one C4 and two factor B (Bf) genes, to the shark MHC was performed using 56 siblings from two typing panels of T. scyllium and Ginglymostoma cirratum. The C4 and one of two Bf genes showed a perfect cosegregation with the class I and II genes, whereas two recombinants were identified for the other Bf gene. These results indicate that the linkage between the complement C4 and Bf genes, as well as the linkage between these complement genes and the MHC class I and II genes were established before the emergence of cartilaginous fish >460 million years ago.     

Effects of substitutions of lysine and aspartic acid for asparagine at beta 108 and of tryptophan for valine at alpha 96 on the structural and functional properties of human normal adult hemoglobin: roles of alpha 1 beta 1 and alpha 1 beta 2 subunit interfaces in the cooperative oxygenation process.

Using our Escherichia coli expression system, we have produced five mutant recombinant (r) hemoglobins (Hbs): r Hb (alpha V96 W), r Hb Presbyterian (beta N108K), r Hb Yoshizuka (beta N108D), r Hb (alpha V96W, beta N108K), and r Hb (alpha V96W, beta N108D). These r Hbs allow us to investigate the effect on the structure-function relationship of Hb of replacing beta 108Asn by either a positively charged Lys or a negatively charged Asp as well as the effect of replacing alpha 96Val by a bulky, nonpolar Trp. We have conducted oxygen-binding studies to investigate the effect of several allosteric effectors on the oxygenation properties and the Bohr effects of these r Hbs. The oxygen affinity of these mutants is lower than that of human normal adult hemoglobin (Hb A) under various experimental conditions. The oxygen affinity of r Hb Yoshizuka is insensitive to changes in chloride concentration, whereas the oxygen affinity of r Hb Presbyterian exhibits a pronounced chloride effect. r Hb Presbyterian has the largest Bohr effect, followed by Hb A, r Hb (alpha V96W), and r Hb Yoshizuka. Thus, the amino acid substitution in the central cavity that increases the net positive charge enhances the Bohr effect. Proton nuclear magnetic resonance studies demonstrate that these r Hbs can switch from the R quaternary structure to the T quaternary structure without changing their ligation states upon the addition of an allosteric effector, inositol hexaphosphate, and/or by reducing the temperature. r Hb (alpha V96W, beta N108K), which has the lowest oxygen affinity among the hemoglobins studied, has the greatest tendency to switch to the T quaternary structure. The following conclusions can be derived from our results: First, if we can stabilize the deoxy (T) quaternary structure of a hemoglobin molecule without perturbing its oxy (R) quaternary structure, we will have a hemoglobin with low oxygen affinity and high cooperativity. Second, an alteration of the charge distribution by amino acid substitutions in the alpha 1 beta 1 subunit interface and in the central cavity of the hemoglobin molecule can influence the Bohr effect. Third, an amino acid substitution in the alpha 1 beta 1 subunit interface can affect both the oxygen affinity and cooperativity of the oxygenation process. There is communication between the alpha 1 beta 1 and alpha 1 beta 2 subunit interfaces during the oxygenation process. Fourth, there is considerable cooperativity in the oxygenation process in the T-state of the hemoglobin molecule.     

Structural analysis of fish versus mammalian hemoglobins: effect of the heme pocket environment on autooxidation and hemin loss

The underlying stereochemical mechanisms for the dramatic differences in autooxidation and hemin loss rates of fish versus mammalian hemoglobins (Hb) have been examined by determining the crystal structures of perch, trout IV, and bovine Hb at high and low pH. The fish Hbs autooxidize and release hemin approximately 50- to 100-fold more rapidly than bovine Hb. Five specific amino acid replacements in the CD corner and along the E helix appear to cause the increased susceptibility of fish Hbs to oxidative degradation compared with mammalian Hbs. Ile is present at the E11 helical position in most fish Hb chains whereas a smaller Val residue is present in all mammalian alpha and beta chains. The larger IleE11 side chain sterically hinders bound O(2) and facilitates dissociation of the neutral superoxide radical, enhancing autooxidation. Lys(E10) is found in most mammalian Hb and forms favorable electrostatic and hydrogen bonding interactions with the heme-7-propionate. In contrast, Thr(E10) is present in most fish Hbs and is too short to stabilize bound heme, and causes increased rates of hemin dissociation. Especially high rates of hemin loss in perch Hb are also due to a lack of electrostatic interaction between His(CE3) and the heme-6 propionate in alpha subunits whereas this interaction does occur in trout IV and bovine Hb. There is also a larger gap for solvent entry into the heme crevice near beta CD3 in the perch Hb (approximately 8 A) compared with trout IV Hb (approximately 6 A) which in turn is significantly higher than that in bovine Hb (approximately 4 A) at low pH. The amino acids at CD4 and E14 differ between bovine and the fish Hbs and have the potential to modulate oxidative degradation by altering the orientation of the distal histidine and the stability of the E-helix. Generally rapid rates of lipid oxidation in fish muscle can be partly attributed to the fact that fish Hbs are highly susceptible to oxidative degradation.     

Contribution of alpha140Tyr and beta37Trp to the near-UV CD spectra on quaternary structure transition of human hemoglobin A

The CD band of human adult hemoglobin (Hb A) at 280 approximately 290 nm shows a pronounced change from a small positive band to a definite negative band on the oxy (R) to deoxy (T) structural transition. This change has been suggested to be due to environmental alteration of Tyrs (alpha42, alpha140, and beta145) or beta37 Trp residues located at the alpha1beta2 subunit interface by deoxygenation. In order to evaluate contributions of alpha140Tyr and beta37Trp to this change of CD band, we compared the CD spectra of two mutant Hbs, Hb Rouen (alpha140Tyr-->His) and Hb Hirose (beta37Trp-->Ser) with those of Hb A. Both mutant Hbs gave a distinct, but smaller negative CD band at 287nm in the deoxy form than that of deoxyHb A. Contributions of alpha140Tyr and beta37Trp to the negative band at the 280 approximately 290 nm region were estimated from difference spectra to be 30% and 26%, respectively. These results indicate that the other aromatic amino acid residues, alpha42Tyr and beta145Tyr, at the alpha1beta2 interface, are also responsible for the change of the CD band upon the R-->T transition of Hb A.     

Discrimination between adaptive and neutral amino acid substitutions in vertebrate hemoglobins

Summary A discriminant analysis on the basis of the physicochemical properties of amino acid residues is developed to investigate the accumulation pattern of amino acid substitutions in a family of proteins. The application of this analysis to vertebrate hemoglobins reveals the following new results. (1) The major components of teleost fish and amphibian hemoglobins showing the Root effect are sharply discriminated from mammalian hemoglobins in several regions of the and chains, whereas shark, minor components of teleost fish and amphibian, reptile, and bird hemoglobins showing no Root effect exhibit a gradual change to mammalian hemoglobin in a straightforward way. This result suggests at least two lines of molecular evolution in vertebrate hemoglobins. (2) The nonadult hemoglobin chains are allocated to the latter line, i.e., tadpole, , and chains are similar to shark and trout I chains, and and chains are similar to some of the reptile chains. (3) In any case, most of the amino acid residues causing the discrimination are located near the sites that carry the amino acid residues conserved well throughout all classes of vertebrates, suggesting that modifications adapting to the respective living conditions or respiratory organs have taken place effectively near the amino acid residues essential for the manifestation of cooperative oxygen binding. (4) The amino acid residues at other sites are changed from one to another species even within the same class, showing a constant substitution rate as a whole. These amino acid substitutions may be nearly neutral, being under a weak functional constraint. The number of sites allowing such neutral substitutions is rather small, less than one-half of all the sites in the adult hemoglobins of bird and mammal, whereas it amounts to two-thirds in teleost fish hemoglobins.     

Hemoglobin structure/function and globin-gene evolution in the Arctic fish Liparis tunicatus

The importance of the Arctic, in contributing to the knowledge of the overall ensemble of adaptive processes influencing the evolution of marine organisms, calls for investigations on molecular adaptations in Arctic fish. Unlike the vast majority of Antarctic Notothenioidei, several Arctic species display high hemoglobin multiplicity. The blood of four species, the spotted wolffish of the family Anarhichadidae and three Gadidae, contains three functionally distinct major components. Similar to many Antarctic notothenioids, Arctic Liparis tunicatus (suborder Cottoidei, family Liparidae) has one major hemoglobin (Hb 1) accompanied by a minor component (Hb 2). This paper reports the structural and functional characterisation of Hb 1 of L. tunicatus. This hemoglobin shows low oxygen affinity, and pronounced Bohr and Root effects. The amino-acid sequence of the beta chain displays an unusual substitution in NA2 (beta2) at the phosphate-binding site, and the replacement of Val E11 (beta67) with Ile. Similar to some Antarctic fish Hbs, electron paramagnetic resonance spectra reveal the formation of a ferric penta-coordinated species even at physiological pH. The amino-acid sequences have also been used to gain insight into the evolutionary history of globins of polar fish. L. tunicatus globins appear close to the notothenioid clades as predicted by teleostean phylogenies. Close phylogenetic relationships between Cottoidei and Notothenioidei, together with their life style, seem to be the main factor driving the globin-sequence evolution.     

Site-directed mutagenesis in hemoglobin: test of functional homology of the F9 amino acid residues of hemoglobin alpha and beta chains

The cysteine residue at F9(93) of the human hemoglobin (Hb A) beta chain, conserved in mammalian and avian hemoglobins, is located near the functionally important alpha1-beta2 interface and C-terminal region of the beta chain and is reactive to sulfhydryl reagents. The functional roles of this residue are still unclear, although regulation of local blood flow through allosteric S-nitrosylation of this residue is proposed. To clarify the role of this residue and its functional homology to F9(88) of the alpha chain, we measured oxygen equilibrium curves, UV-region derivative spectra, Soret-band absorption spectra, the number of titratable -SH groups with p-mercuribenzoate and the rate of reaction of these groups with 4, 4'-dipyridine disulfide for three recombinant mutant Hbs with single amino acid substitutions: Ala-->Cys at 88alpha (rHb A88alphaC), Cys-->Ala at 93beta (rHb C93betaA) and Cys-->Thr at 93beta (rHb C93betaT). These Hbs showed increased oxygen affinities and impaired allosteric effects. The spectral data indicated that the R to T transition upon deoxygenation was partially restricted in these Hbs. The number of titratable -SH groups of liganded form was 3.2-3.5 for rHb A88alphaC compared with 2.2 for Hb A, whereas those for rHb C93betaA and rHb C93betaT were negligibly small. The reduction of rate of reaction with 4,4'-dipyridine disulfide upon deoxygenation in rHb A88alphaC was smaller than that in Hb A. Our experimental data have shown that the residues at 88alpha and 93beta have definite roles but they have no functional homology. Structure-function relationships in our mutant Hbs are discussed.     

Ancient vertebrate conserved noncoding elements have been evolving rapidly in teleost fishes

Vertebrate genomes contain thousands of conserved noncoding elements (CNEs) that often function as tissue-specific enhancers. In this study, we have identified CNEs in human, dog, chicken, Xenopus, and four teleost fishes (zebrafish, stickleback, medaka, and fugu) using elephant shark, a cartilaginous vertebrate, as the base genome and investigated the evolution of these ancient vertebrate CNEs (aCNEs) in bony vertebrate lineages. Our analysis shows that aCNEs have been evolving at different rates in different bony vertebrate lineages. Although 78-83% of CNEs have diverged beyond recognition ("lost") in different teleost fishes, only 24% and 40% have been lost in the chicken and mammalian lineages, respectively. Relative rate tests of substitution rates in CNEs revealed that the teleost fish CNEs have been evolving at a significantly higher rate than those in other bony vertebrates. In the ray-finned fish lineage, 68% of aCNEs were lost before the divergence of the four teleosts. This implicates the "fish-specific" whole-genome duplication in the accelerated evolution and the loss of a large number of both copies of duplicated CNEs in teleost fishes. The aCNEs are rich in tissue-specific enhancers and thus many of them are likely to be evolutionarily constrained cis-regulatory elements. The rapid evolution of aCNEs might have affected the expression patterns driven by them. Transgenic zebrafish assay of some human CNE enhancers that have been lost in teleosts has indicated instances of conservation or changes in trans-acting factors between mammals and fishes.     

Functionally distinct haemoglobins of the cryopelagic Antarctic teleost Pagothenia borchgrevinki

Pagothenia borchgrevinki , has a higher haemoglobin concentration than other Antarctic notothenioids and the high oxygen capacity may correlate with the relatively active mode of life of this fish. The fish has five haemoglobins (Hb C, Hb 0, Hb 1, Hb 2 and Hb 3) with Hb 1 accounting for 70–80% of the total, and Hb C being present in trace amounts. Hb 1 and Hb 2 are functionally similar in terms of Bohr and Root effects. Hb 3 has a weaker Bohr effect than Hb 1 and Hb 2, and the Root effect is similar to that of Hb 1. Hb 0 has a strong Bohr effect and the Root effect is enhanced to a larger extent by the physiological effectors chlorides and phosphates than that of the other components with the exception of Hb C. The heats of oxygenation are lower than those of temperate fish haemoglobins. Temperature variations may have a different effect on the functional properties of each haemoglobin, and chloride and phosphates may play an important role in the conformational change between the oxy and deoxy structures. The complete amino acid sequences of Hb 1 and Hb 0, as well as partial N-terminal or internal sequences of the other haemoglobins, have been established. The high multiplicity of functionally distinct haemoglobins indicates that P. borchgrevinki , has a specialized haemoglobin system.     

Enhanced inhibition of polymerization of sickle cell hemoglobin in the presence of recombinant mutants of human fetal hemoglobin with substitutions at position 43 in the gamma-chain

Four recombinant mutants of human fetal hemoglobin [Hb F (alpha2gamma2)] with amino acid substitutions at the position 43 of the gamma-chain, rHb (gammaD43L), rHb (gammaD43E), rHb (gammaD43W), and rHb (gammaD43R), have been expressed in our Escherichia coli expression system and used to investigate their inhibitory effect on the polymerization of deoxygenated sickle cell hemoglobin (Hb S). Oxygen-binding studies show that rHb (gammaD43E), rHb (gammaD43W), and rHb (gammaD43R) exhibit higher oxygen affinity than human normal adult hemoglobin (Hb A), Hb F, or rHb (gammaD43L), and all four rHbs are cooperative in binding O2. Proton nuclear magnetic resonance (NMR) studies of these four rHbs indicate that the quaternary and tertiary structures around the heme pockets are similar to those of Hb F in both deoxy (T) and liganded (R) states. Solution light-scattering experiments indicate that these mutants remain mostly tetrameric in the liganded (R) state. In equimolar mixtures of Hb S and each of the four rHb mutants (gammaD43L, gammaD43E, gammaD43R, and gammaD43W), the solubility (Csat) of each of the pairs of Hbs is higher than that of a similar mixture of Hb S and Hb A, as measured by dextran-Csat experiments. Furthermore, the Csat values for Hb S/rHb (gammaD43L), Hb S/rHb (gammaD43E), and Hb S/rHb (gammaD43R) mixtures are substantially higher than that for Hb S/Hb F. The results suggest that these three mutants of Hb F are more effective than Hb F in inhibiting the polymerization of deoxy-Hb S in equimolar mixtures.     

A study on the quaternary structure change of hemoglobin in the ligation process

In order to inquire into the molecular mechanism underlying the cooperative ligand binding to hemoglobin (Hb), conformational interaction at the interfaces between subunits are investigated on the basis of the atomic coordinates of human deoxy and human carbonmonoxy Hbs. Hypothetical intermediate structures are used, each of which is obtained from the procedure where one or more subunits in deoxy Hb are replaced by the corresponding CO-liganded subunits in carbonmonoxy Hb using the method of superimposition of two sets of atomic coordinates. When either alpha or beta subunit is substituted with the corresponding subunit in carbonmonoxy Hb, serious steric hindrances are produced between alpha 1FG4(92)Arg and beta 2C3(37)Trp or between alpha 1C6(41)Thr and beta 2FG4(97)His, all of which belong to the allosteric core affected directly by ligand binding. These steric hindrances become more serious when both alpha 1(alpha 2) and beta 2(beta 1) subunits are substituted. Therefore the change in the relative distance between iron atom and porphyrin by ligation results in strain in the C-terminal residues as an effect of the steric hindrance between the FG and C segments. However, no steric hindrance can be seen between subunits when the subunits in carbonmonoxy Hb are substituted with the corresponding subunits in deoxy Hb. The nature of the quaternary structural change from liganded to deoxy Hb seems to be different from that from deoxy to liganded Hb.     

Molecular evolution of shark and other vertebrate DNases I.

We purified pancreatic deoxyribonuclease I (DNase I) from the shark Heterodontus japonicus using three-step column chromatography. Although its enzymatic properties resembled those of other vertebrate DNases I, shark DNase I was unique in being a basic protein. Full-length cDNAs encoding the DNases I of two shark species, H. japonicus and Triakis scyllia, were constructed from their total pancreatic RNAs using RACE. Nucleotide sequence analyses revealed two structural alterations unique to shark enzymes: substitution of two Cys residues at positions 101 and 104 (which are well conserved in all other vertebrate DNases I) and insertion of an additional Thr or Asn residue into an essential Ca(2+)-binding site. Site-directed mutagenesis of shark DNase I indicated that both of these alterations reduced the stability of the enzyme. When the signal sequence region of human DNase I (which has a high alpha-helical structure content) was replaced with its amphibian, fish and shark counterparts (which have low alpha-helical structure contents), the activity expressed by the chimeric mutant constructs in transfected mammalian cells was approximately half that of the wild-type enzyme. In contrast, substitution of the human signal sequence region into the amphibian, fish and shark enzymes produced higher activity compared with the wild-types. The vertebrate DNase I family may have acquired high stability and effective expression of the enzyme protein through structural alterations in both the mature protein and its signal sequence regions during molecular evolution.     

Biochemical Characterization of Hemoglobins from Caspian Sea Sturgeons (<Emphasis Type="Italic">Acipenser persicus and Acipenser stellatus</Emphasis>)

Hemoglobin (Hb) variability is a commonly used index of phylogenetic differentiation and molecular adaptation in fish enabling them to adapt to different ecological conditions. In this study, the characteristics of Hbs from two Sturgeon species of the Southern Caspian Sea Basin were investigated. After extraction and separation of hemoglobin from whole blood, the polyacrylamide gel electrophoresis (SDS-PAGE), cellulose acetate electrophoresis, and isoelectric focusing (IEF) were used to confirm Hb variabilities in these fishes. We showed that although both species have variable Hbs with different isoelectric points, their dominant Hbs can be identified. Ion exchange on CM-cellulose chromatography was used for purification of the dominant Hbs from these fishes. The accuracy of the methods was confirmed by IEF and SDS-PAGE. Spectral studies using fluorescence spectrophotometery indicated that although the Hbs from these fishes had similar properties they exhibited clear differences with human Hb. A comparative study of Hbs alpha-helix secondary substructures was performed by circular dichroism spectropolarimetry (CD) analysis. UV–vis spectrophotometery was also utilized to measure oxygen affinity of Hbs by sodium dithionite. Oxygen affinities of these Hbs were compared using Hb–oxygen dissociation curves. Together, these results demonstrate a significant relationship between oxygen affinity of fish hemoglobins and environmental partial pressure of oxygen.     

The evolution of multiple isotypic IgM heavy chain genes in the shark

The IgM H chain gene organization of cartilaginous fishes consists of 15-200 miniloci, each with a few gene segments (V(H)-D1-D2-J(H)) and one C gene. This is a gene arrangement ancestral to the complex IgH locus that exists in all other vertebrate classes. To understand the molecular evolution of this system, we studied the nurse shark, which has relatively fewer loci, and characterized the IgH isotypes for organization, functionality, and the somatic diversification mechanisms that act upon them. Gene numbers differ slightly between individuals ( approximately 15), but five active IgM subclasses are always present. Each gene undergoes rearrangement that is strictly confined within the minilocus; in B cells there is no interaction between adjacent loci located > or =120 kb apart. Without combinatorial events, the shark IgM H chain repertoire is based on junctional diversity and, subsequently, somatic hypermutation. We suggest that the significant contribution by junctional diversification reflects the selected novelty introduced by RAG in the early vertebrate ancestor, whereas combinatorial diversity coevolved with the complex translocon organization. Moreover, unlike other cartilaginous fishes, there are no germline-joined VDJ at any nurse shark mu locus, and we suggest that such genes, when functional, are species-specific and may have specialized roles. With an entire complement of IgM genes available for the first time, phylogenetic analyses were performed to examine how the multiple Ig loci evolved. We found that all domains changed at comparable rates, but V(H) appears to be under strong positive selection for increased amino acid sequence diversity, and surprisingly, so does Cmicro2.