Slow Rates of Evolution and Sequence Homogenization in an Ancient Satellite DNA Family of Sturgeons

Roberto de la Herrán*, Francesco Fontana, Massimo Lanfredi, Leonardo Congiu, Marilena Leis, Remigio Rossi, Carmelo Ruiz Rejón, Manuel Ruiz Rejón and Manuel A. Garrido-Ramos

*Departamento de Genética, Facultad de Ciencias, Universidad de Granada, Granada, Spain;
{dagger}Department of Biology, University of Ferrara, Ferrara, Italy

Satellite DNA sequences are noncoding, short, tandemly repeated sequences located mainly at centromeres and telomeres (Charlesworth, Sniegowski, and Stephan 1994Citation ). The main concerns of evolutionary patterns of satellite DNA sequences are their rapid change rate and their concerted mode of evolution. When a satellite DNA family is shared by several related species, this enables one to find specific diagnostic sites for use in phylogenetic studies (Arnason, Grettarsdottir, and Widegren 1992Citation ; Gretarsdottir and Arnason 1992Citation ). Here, we analyze an exceptional ancient satellite DNA family of sturgeons which does not fit the evolutionary features of most satellite DNAs. Specifically, we found that the evolution of this satellite DNA family is very slow, as reflected by low mutation and homogenization rates. Notwithstanding, our data offer some information on the tracing of mechanisms accounting for the formation and origin of this satellite DNA family and on the phylogeny of the sturgeons.

This satellite DNA family, the HindIII family, was previously found within the genome of Acipenser naccarii (Garrido-Ramos et al. 1997Citation ). Now, we extend this analysis in an attempt to determine the presence of this satellite DNA and its evolutionary pattern at the molecular level in the genomes of six species of Acipenser (the Adriatic sturgeon A. naccarii, the Siberian sturgeon Acipenser baerii, the Russian sturgeon Acipenser gueldenstaedtii, the sterlet Acipenser ruthenus, the white sturgeon Acipenser transmontanus and the Atlantic sturgeon Acipenser sturio) and one species of Huso (the great sturgeon Huso huso). Three individuals from each species were analyzed.

By conventional cloning or by PCR amplification (primers 5'-CGAACCTGTAAGCTT-3' and 5'-TGATCTTCAGAACTACC-3', designed from conserved parts of the HindIII satellite DNA) and posterior cloning (Garrido-Ramos et al. 1995, 1999Citation ), we obtained a total of 52 monomeric units of the HindIII satellite DNA: 11 from A. naccarii, 9 from A. baerii, 12 from A. gueldenstaedtii, 9 from A. transmontanus, 4 from A. ruthenus, and 7 from H. huso. The EMBL accession numbers for these sequences are AJ286564–AJ286615 and Z49941. Acipenser sturio proved not to have the HindIII satellite DNA sequences. This extreme has also been confirmed by in situ hybridization techniques (unpublished data).

According to biogeographic and molecular data on phylogenetic analyses of the sturgeon species (Birstein and deSalle 1998Citation ), our results indicate that we have analyzed an ancient satellite DNA more than 80 Myr old. As discussed below, since that time, the HindIII sequence pattern change has been very slow. This allows us to infer major abrupt changes that could have occurred in the process of formation and spreading of the HindIII satellite DNA. In all of the species, the length of the satellite monomers varied between 169 and 172 bp. However, in all of the species, the current monomer length appeared to be the consequence of several duplication events from smaller subunits (fig. 1A ). Analyzing the characteristics of this satellite DNA family, we propose the following events to explain the apparition of the current repeat unit of the HindIII satellite DNA family (fig. 1B ): (1) several partial duplications of the sequences AAATTAT and GGACC led to the amplification of a initial monomer of about 50 bp, generating a monomer of about 85 bp in length; (2) one entire duplication event led to the amplification of the monomer length from 85 bp to 170 bp; after that, (3) sequence divergence occurred between the two parts of the new monomer (these three events must have occurred early in the evolution of HindIII sequences, as they are shared by all species analyzed); finally, (4) a new event of partial duplication occurred, now of the sequence CTTTT in positions 1–5 and 11–15. This event must be recent (about 15 MYA according to Birstein and deSalle [1998Citation ]). This duplication could lead to the appearance of two HindIII subfamilies (the CTTTT subfamily and the AGCTC subfamily) in A. naccarii, A. baerii, and A. gueldenstaedtii, which currently coexist in the process of homogenization. The appearance of the CTTTT subfamily after the AGCTC subfamily is supported by the fact that the duplication event is observed only in the first part of the sequence and only in three species (apparently closely related), while it is absent in the other species, such as A. transmontanus, which has the sequence AGCTC (or variants) at positions 11–15. Although monomers having the CTTTT sequence in this region are constant (presumably because of their recent origin), the alternative subfamily is highly variable, with monomers having AGCTC, AGTTC, or AGAAC. This polymorphism is ancient, and it is shared by all the species having HindIII sequences, except A. naccarii, in which the AGCTC variant has apparently been fixed or almost fixed.



View larger version (29K):
[in this window]
[in a new window]
 
Fig. 1.—A, Sequence comparison between the two distinguishable parts of the HindIII monomers. Arrows indicate short tandem direct repeats within each part. See text for details. The sequence depicted in this figure is a consensus sequence obtained from the comparison of all 52 monomeric units analyzed. All of the species had similar HindIII sequences and the same subrepeat organization as that shown in this figure. B, Scheme indicating the possible origin and formation of HindIII sequences from an initial monomer of about 50 bp to a current monomer of 170 bp throughout several partial duplication events and one whole duplication event. See text for details.

 
The latter variability for A. naccarii, A. baerii, and A. gueldenstaedtii is the most significant detected in the HindIII satellite DNA of these species. Apart from the variable region, the intraspecific variabilities (nucleotide diversity; Nei 1987Citation ) for these three species were 2.27%, 8.17%, and 4.27%, respectively. Aside from the variable region between positions 11 and 15, intraspecific sequence variabilities in A. transmontanus, A. ruthenus, and H. huso were 6.02%, 5.71%, and 10.65%, respectively.

The 52 monomeric units sequenced were analyzed using the MEGA package (Kumar, Tamura, and Nei 1993Citation ). Genetic distances within and between species were calculated according to the Jukes-Cantor method. UPGMA and neighbor-joining trees showed similar topologies (fig. 2 ). Two main branches are visible: one leading to sequences of A. naccarii, A. baerii, A. gueldenstaedtii, and A. transmontanus, and another leading to sequences of A. ruthenus and H. huso, although some sequences of A. baerii are intermixed in this second branch. Two main features concerning the clustering are evident. First, the sequences appear intermixed without reflecting taxonomic affinity: this reflects the fact that many repeat variants are shared by different species and that genetic distances between monomers of different species are similar to or even shorter than genetic distances between monomers within species. Actually, by performing an analysis of molecular variance (AMOVA) (Excoffier, Smouse, and Quattro 1992Citation ) within and between species using ARLEQUIN software (Schneider et al. 1996Citation ), we found that 70% of the total genetic diversity could be ascribed to intraspecific variability, and the remaining 30% represented the amount of differentiation between species. Furthermore, within the first main branch, two groups of sequences corresponding to the two sequence subfamilies (the CTTTT and the AGCTC subfamilies, or variants) could be identified.



View larger version (21K):
[in this window]
[in a new window]
 
Fig. 2.—Neighbor-joining tree of the 52 monomeric HindIII sequences analyzed in this study (pHAn = clones of Acipenser naccarii; pGU = clones of Acipenser gueldenstaedtii; pBa = clones of Acipenser baerii; pTR = clones of Acipenser transmontanus; pRu = clones of Acipenser ruthenus; pHH = clones of Huso huso). Numbers are bootstrapping indices for the levels of support for individual nodes (5,000 replicates)

 
In general, satellite DNAs follow concerted evolution, with very few exceptions: the best known case is that of a satellite DNA in tsetse flies (Trick and Dover 1984Citation ), in which absence of concerted evolution was attributed to low rates of sequence homogenization. However, in this case, both the intraspecific and the interspecific sequence divergences were high. In contrast, the HindIII satellite DNA of sturgeons is characterized by low intraspecific variability and also by similar or even lower interspecific sequence divergence. Thus, the explanation for HindIII sequences may be related to the fact that not only the rate of homogenization, but also the rate of mutation is low. With respect to the rate of mutation, slow sequence evolution appears to be the norm for sturgeon DNA, either nuclear or mitochondrial (Kedrova, Vladychenskaia, and Antonov 1980Citation ; Birstein and deSalle 1998Citation ; Ludwig and Kirschbaum 1998Citation ; Tagliavini et al. 1999aCitation ), and for fish DNA in general (Martin, Naylor, and Palumbi 1992Citation ; Rico, Rico, and Hewitt 1996Citation ), but also for aquatic mammals and cold-blooded animals. To explain the unusually high conservation of sequences flanking microsatellites in whales, Scholötter, Amos, and Tautz (1991)Citation suggested that aquatic environments are less mutagenic than terrestrial ones. Martin, Naylor, and Palumbi (1992)Citation and Rico, Rico, and Hewitt (1996)Citation suggested that the observations that metabolic rates in poikilotherms are lower and that rates of DNA damage are proportional to specific metabolic rates (Schmidt-Nielsen, 1986Citation ; Shigenaga, Gimeno, and Ames 1989Citation ) could explain the lower rates of nucleotide substitution for mitochondrial DNA found in sharks and for microsatellite flanking sequences found in turtles and fish.

The homogenization rates should be interpreted in terms of parameters influencing homogenization in the chromosome pool, such as generation time and effective population size (Ohta and Dover 1983, 1984Citation ). Differences in homogenization rates between multigene families could be also related to functional requirements (Gojorobi and Nei 1984Citation ) and the structure of the repeat units (Zimmer et al. 1980Citation ). However, in the case of the HindIII sequences of sturgeons, we could explain the low homogenization rate among single species in terms of the primary rates of the homogenization process. That is, it is possible that the exchange between nonhomologous chromosomes having HindIII sequences is limited. The HindIII sequences are restricted to the centromeres of a small percentage of chromosomes in these species (unpublished data). However, these chromosomes are very heterogeneous in shape and size, and this could hamper interchromosomal exchange and homogenization (Jantsch et al. 1990Citation ). A fish family of more recent origin than sturgeons, the Sparidae (50–35 MYA; Tortonese 1975Citation ), shows a satellite DNA family which is localized in centromeres of all chromosomes in all species and follows concerted evolution (Garrido-Ramos et al. 1995, 1999Citation ). However, one difference between sparids and sturgeons is that the karyotypes of sparids are rather symmetrical, and there is thus no evident physical limitation to interchromosomal exchange. Actually, sparid species in which karyotypes are less symmetrical (for example, Diplodus bellotii) show higher intraspecific variability (Garrido-Ramos et al. 1999Citation ). However, we cannot disregard the possibility that events of interspecific hybridizations between sturgeon species followed by polyploidization (i.e., reticular evolution; Vasil'ev 1999Citation ) could preclude interspecific sequence divergence for HindIII satellite DNA.

Our data offer some clues concerning the controversial phylogenetic relationships of the species of the Acipenserinae subfamily of Acipenseridae. The HindIII satellite DNA is shared by six of the seven species analyzed here (A. naccarii, A. gueldenstaedtii, A. baerii, A. transmontanus, A. ruthenus, and Huso huso). If we assume that the sharing of a satellite DNA family is a strong signal of cladistic relationships, we find that H. huso is a species belonging to the Acipenser group, supporting the view of Birstein and deSalle (1998)Citation and the doubts (Artyukhin 1995Citation ) about the systematic position of the apparently artificial genus Huso. Furthermore, we demonstrated that the HindIII sequences are absent from the genome of A. sturio. These data support the hypothesis that A. sturio differs from the other species and most likely had an independent evolution (Birstein and deSalle 1998Citation ).

We found that the HindIII sequences were identical between A. gueldenstaedtii and A. naccarii and almost identical to those of A. baerii (although a significant proportion of its sequence variants are shared with species such as A. ruthenus and H. huso). Presumably, the first three species are closely related (Tagliavini et al. 1999aCitation ; Artyukhin 1995Citation ) and, together with A. transmontanus, must represent a group of species separated from the other group made up of A. ruthenus and H. huso (Tagliavini et al. 1999aCitation ). It is striking that H. huso and A. ruthenus are diploid species and the rest of species are tetraploid species (Fontana 1994Citation ; Fontana et al. 1998Citation ; Tagliavini et al. 1999bCitation ).

Acknowledgements

This work was supported by grants of the FEDER funds (1FD97-2289) and of Plan Andaluz de Investigación (Group CVI0200). This work was also supported by the Fourth Triennial Plan "Aquaculture in Seas and Lagoons" of the Italian Ministry of Agriculture, Food Resources and Forestry (4C14t), and by research grants from the Italian Ministry of University and Research. We are grateful for the valuable comments and suggestions on the manuscript given by two anonymous referees.

Footnotes

Pierre Capy, Reviewing Editor

1 Keywords: sturgeons satellite DNA concerted evolution slow mutation rates phylogeny Back

2 Address for correspondence and reprints: Roberto de la Herrán, Departamento de Genética, Facultad de Ciencias, Universidad de Granada, 18071 Granada, Spain. E-mail: rherran{at}ugr.es Back

literature cited

    Arnason, U., S. Grettarsdottir, and B. Widegren. 1992. Mysticete (baleen whale) relationships based upon the sequence of the common cetacean DNA satellite. Mol. Biol. Evol. 9:1018–1028.[Abstract]

    Artyukhin, E. N. 1995. On biogeography and relationships within the genus Acipenser. Sturg. Q. 3:6–8.

    Birstein, V. J., and R. deSalle. 1998. Molecular phylogeny of Acipenserinae. Mol. Phylogenet. Evol. 9:141–155.[ISI][Medline]

    Charlesworth, B., P. Sniegowski, and W. Stephan. 1994. The evolutionary dynamics of repetitive DNA in eukaryotes. Nature 371:215–220.

    Excoffier, L., P. E. Smouse, and J. M. Quattro. 1992. Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics 131:479–491.

    Fontana, F. 1994. Chromosomal nucleolar organizer regions in four sturgeon species as markers of karyotype evolution in Acipenseriformes (Pisces). Genome 37:888–892.

    Fontana, F., J. Tagliavini, L. Congiu, M. Lanfredi, M. Chicca, C. Laurente, and R. Rossi. 1998. Karyotypic characterization of the great sturgeon, Huso huso, by multiple staining techniques and fluorescent in situ hybridization. Mar. Biol. 132:495–501.[ISI]

    Garrido-Ramos, M. A., R. de la Herrán, M. Jamilena, R. Lozano, C. Ruiz Rejón, and M. Ruiz Rejón. 1999. Evolution of centromeric satellite DNA and its use in phylogenetic studies of the Sparidae family (Pisces, Perciformes). Mol. Phylogenet. Evol. 12:200–204.[ISI][Medline]

    Garrido-Ramos, M. A., M. Jamilena, R. Lozano, C. Ruiz Rejón, and M. Ruiz Rejón. 1995. The EcoRI centromeric satellite DNA of the Sparidae family (Pisces, Perciformes) contains a sequence motive common to other vertebrate centromeric satellite-DNA. Cytogenet. Cell Genet. 71:345–351.[ISI][Medline]

    Garrido-Ramos, M. A., M. C. Soriguer, R. de la Herrán, M. Jamilena, C. Ruiz Rejón, A. Domezain, J. A. Hernando, and M. Ruiz Rejón. 1997. Morphometric and genetic analysis as proof of the existence of two sturgeon species in the Guadalquivir river. Mar. Biol. 129:33–39.[ISI]

    Gojorobi, T., and M. Nei. 1984. Concerted evolution of the immunoglobulin VH gene family. Mol. Biol. Evol. 1:195–211.[Abstract]

    Gretarsdottir, G., and U. Arnason. 1992. Evolution of the common cetacean highly repetitive DNA component and the systematic position of Orcaella brevirostris. J. Mol. Evol. 34:201–208.[ISI][Medline]

    Jantsch, M., B. Hamilton, B. Mayr, and D. Schweizer. 1990. Meiotic chromosome behaviour reflects levels of sequence divergence in Sus scrofa domestica satellite DNA. Chromosoma 99:330–335.

    Kedrova, O. S., N. S. Vladychenskaia, and A. S. Antonov. 1980. Divergence of unique and repetitive sequences in the genomes of fish. Mol. Biol. (Moscow) 14:1001–1012.

    Kumar, S., K. Tamura, and M. Nei. 1993. MEGA: molecular evolutionary genetics analysis. Version 1.01. Pennsylvania State University, University Park.

    Ludwig, A., and F. Kirschbaum. 1998. Comparison of mitochondrial DNA sequences between the European and the Adriatic sturgeon. J. Fish Biol. 52:1289–1291.[ISI]

    Martin, A. P., G. J. P. Naylor, and S. R. Palumbi. 1992. Rates of mitochondrial DNA in sharks are slow compared with mammals. Nature 357:153–155.

    Nei, M. 1987. Molecular evolutionary genetics. Columbia University Press, New York.

    Ohta, T., and G. Dover. 1983. Population genetics of multigene families that are dispersed into two or more chromosomes. Proc. Natl. Acad. Sci. USA 80:4079–4083.

    ———. 1984. The cohesive population genetics of molecular drive. Genetics 108:501–521.

    Rico, C., I. Rico, and G. Hewitt. 1996. 470 million years of conservation of microsatellite loci among fish species. Proc. R. Soc. Lond. B Biol. Sci. 263:549–557.[ISI][Medline]

    Schlötterer, C., B. Amos, and D. Tautz. 1991. Conservation of polymorphic simple sequence loci in cetacean species. Nature 354:63–65.

    Schmidt-Nielsen, K. 1986. Scalling: why is animal size so important?. Cambridge University Press, New York.

    Schneider, S., J. M. Kueffer, D. Roessli, and L. Excoffier. 1996. Arlequin: a software environment for the analysis of population genetics data. Version 1.0. Genetics and Biometry Laboratory, University of Geneva, Switzerland.

    Shigenaga, M. K., C. J. Gimeno, and B. N. Ames. 1989. Urinary 8-hydroxy-2'-deoxyguanosine as a biological marker of in vivo oxidative DNA damage. Proc. Natl. Acad. Sci. USA 86:9697–9701.

    Tagliavini, J., G. Conterio, G. Gandolfi, and F. Fontana. 1999a. Mitochondrial DNA sequences of six sturgeon species and phylogenetic relationships within Acipenseridae. J. Appl. Ichthyol. 15:17–22.

    Tagliavini, J., P. Williot, L. Congiu, M. Chicca, M. Lanfredi, R. Rossi, and F. Fontana. 1999b. Molecular cytogenetic analysis of the karyotype of the European Atlantic sturgeon, Acipenser sturio. Heredity 83:520–525

    Tortonese, E. 1975. Fauna de Italia "Osteichthyes." Pesci Ossei. Vol. XI. Calderini, Bologna, Italy.

    Trick, M., and G. Dover. 1984. Unexpectedly slow homogenisation within a repetitive DNA family shared between two subspecies of Tsetse fly. J. Mol. Evol. 20:322–329.[ISI][Medline]

    Vasil'ev, V. P. 1999. Polyploidization by reticular speciation in Acipenseriform evolution: a working hypothesis. J. Appl. Ichthyol. 15:29–31.

    Zimmer, E. A., S. L. Martin, S. M. Beverley, Y. W. Kan, and A. C. Wilson. 1980. Rapid duplication and loss of genes coding for the {alpha} chains of hemoglobin. Proc. Natl. Acad. Sci. USA 77:2158–2162.

Accepted for publication September 26, 2000.