8urn:lsid:arphahub.com:pub:A71ED5FC-60ED-5DA3-AC8E-F6D2BB5B3573urn:lsid:zoobank.org:pub:C8FA3ADA-5585-4F26-9215-A520EE683979Comparative CytogeneticsCCG1993-07711993-078XPensoft Publishers10.3897/compcytogen.v7i3.53981796Research ArticleReptiliaSauriaMolecular biologySystematicsTaxonomyItalyKaryological evidence for diversification of Italian slow worm populations (Squamata, Anguidae)MezzasalmaMarcello1marcello.mezzasalma@unina.itGuarinoFabio Maria1ApreaGennaro1Agnese Petraccioli1CrottiniAngelica2https://orcid.org/0000-0002-8505-3050OdiernaGaetano1Dipartimento di Biologia, Università di Napoli Federico II, Via Cinthia I- 80126 Napoli, ItalyUniversità Federico IINaplesItalyCIBIO, Centro de Investigação em Biodiversidade e Recursos Genéticos, Campus Agrário de Vairão, R. Padre Armando Quintas, 4485-661 Vairão, Vila do Conde, PortugalTU BraunschweigItaly
2013239201373217227FFBB1644-BC5E-FFE6-DE56-FFA2FFD9FFA757542822420131872013Marcello Mezzasalma, Fabio Maria Guarino, Gennaro Aprea, Agnese Petraccioli, Angelica Crottini, Gaetano OdiernaThis is an open access article distributed under the terms of the Creative Commons Attribution License 3.0 (CC-BY), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
A karyological analysis on six Italian populations the slow worm (Anguis veronensis Pollini, 1818) was performed and their genetic differentiation at the mitochondrial 16S rRNA gene fragment from a Spanish sample has been assessed. The Italian populations were karyologically uniform, all showing 2n=44 elements, of which 20 were macrochromosomes and 24 microchromosomes. Comparison with literature data on Central European populations showed a difference on the morphology of the 10th chromosome pair: submetacentric in Italian populations and telocentric in the Central European ones. Our analysis showed the presence of a fragile site on chromosomes of this pair, suggesting its propensity for structural rearrangements. Analysis of the 16S rRNA gene fragment showed uniformity among Italian populations (uncorrected genetic distance of 0.4%), and their genetic distinctness from the Spanish individual (uncorrected genetic distance of 4.2%). Our results confirm the existence of two different Anguis fragilis Linnaeus, 1758 lineages, each one characterized by a different cytotype.
Until recentlythere were only two recognized species of the genus Anguis Linnaeus, 1758 in the Palaeartic region: Anguis cephallonica Werner, 1894 and Anguis fragilis Linnaeus, 1758, commonly known as slow worms. The first species was considered a Mediterranean endemic restricted to the Peloponnese and some Ionian islands, while the second (Anguis fragilis) was considered a widespread taxon distributed from Western Europe to NW Iran and from the Mediterranean coast to Scandinavia, with a broad altitudinal distribution, ranging from the sea level up to 2300 m above sea level (Gasc et al. 1997, Sindaco et al. 2006). Two recent studies (Gvoždik et al. 2010, 2013) based on the analysis of mitochondrial and nuclear genes and morphology redefined the distribution and the phylogenetic relationships of various populations of the genus Anguis and identified five main distinct evolutionary lineages: Anguis fragilis sensu stricto (present in Austria, Bosnia and Herzegovina, Bulgaria, Croatia, Czech Republic, France, Germany, Great Britain, Greece, Hungary, Italy, Macedonia, Montenegro, Norway, Poland, Portugal, Serbia, Slovakia, Slovenia, Spain, Sweden and Switzerland), Anguis graeca Bedriaga, 1881 (in Albania, Greece and Montenegro), Anguis cephallonica (limited to the Peloponnese), and Anguis colchica (Nordmann, 1840) (widely distributed from Eastern Europe to Iran) and subdivided into three main lineages (Anguis colchica colchica, Anguis colchica incerta and Anguis colchica orientalis) each with a geographically distinct distribution; and more recently Anguis veronensis Pollini, 1818 (present in Italian peninsula and some areas from south-eastern France). Both the mitochondrial (mtDNA) and the nuclear (nuDNA) analyses analyses identified Anguis cephallonica as the sister lineage of a wider clade comprising Anguis fragilis, Anguis colchica and Anguis graeca (Gvoždik et al. 2010), while the position of the Italian lineage still remains unresolved (Gvoždik et al. 2013). The proposed species overall parapatric distributions, excluding some apparently isolated populations of Anguis fragilis sensu stricto from Greece (Gvoždik et al. 2010). In a previous karyological study (Gigantino et al. 2002) on the herpetofauna from the Matese Regional Park (Campania, Southern Italy) showed that local slow worm populations had a different chromosome formula when compared with available karyological data of Anguis fragilis (Dalq 1921, Margot 1946). Those karyological data, indeed scarce and dated, referred to referred to Central European populations of Anguis fragilis and described a karyotype of 2n=44 elements, of which 20 were macrochromosomes (the first, the fourth and the fifth pairs shaped as metacentric and the rest as telocentric elements) and 24 microchromosomes, arm number (A.N.)=48. The populations of the Matese Regional Park differed in showing as submetacentric the elements of the tenth macrochromosome pair (A.N.=50). On the basis of the recent finding of Gvoždik et al. (2013) and our previous karyological data (Gigantino et al. 2002), we extended the chromosome analysis to other Italian populations. The results of this comparative karyological study are here presented together with the re-worked analysis on populations from Matese Regional Park. In addition, since in Balkan Peninsula two endemic Anguis species were found (Gvoždik et al. 2013), we also perfomed a molecular analysis of a fragment of the mitochondrial 16S rRNA gene in order to test the genetic uniformity of the studied Italian specimens.
Material and methods
Eleven individuals from six geographically distinct Italian localities were analysed in this work: Travacò Siccomario (Pavia, northern Italy), one male; Ancona (central Italy), two males; Valle Agricola (Caserta, southern Italy), one male and one female; Giffoni Valle Piana (Salerno, southern Italy), one male and one female; Ruvo del Monte (Potenza, southern Italy), two males and one females; Monte Cocuzzo (Cosenza, southern Italy), one males. As already successfully performed (Gigantino et al. 2002), adult individuals were sexed in the field by means of hemipenis extroversion before their release. For the samples of the Matese Regional Park we benefitted from stored chromosome suspensions obtained from organs (see Odierna et al. 2004 for details) kept in the herpetological collection of one of the co-author (G.F.M.). For all other samples, chromosomes were obtained by establishing blood cultures using a freshly collected blood sample (100–500 µl), drawn on the field from the caudal vein, incubated for 5 days in 3 ml of peripheral blood karyotyping medium (Biological Industries). Colcemid (100 ng/ml) was added to culture two hours before harvesting. After washing in
Phosphate Buffer Saline
(PBS) 1x, and incubation in hypotonic solution (KCl 0,075 M) per 30 min, cells were fixed in 3:1 methanol-acetic acid. This non-invasive method, useful for both karyologycal and molecular investigations, allows the avoidance of specimen sacrifice and animals were released at capture sites immediately after the sampling. In addition to standard staining method (5% Giemsa solution at pH 7) we performed various chromosome banding techniques: Ag-Nor banding (Howel and Black 1980), C-banding (Sumner 1972), C-banding+CMA3+DAPI staining (Odierna et al. 2004), Chromomycin A3-methyl green staining (Odierna et al. 2007); G banding performed on 10 days old slides by a brief incubation (10–30 sec) in a 0,05% trypsin solution (Odierna et al. 1994). For
Late Replicating banding pattern
(LR), bromedeoxyuridine (BrdU; 35 µg/ml) was added to blood cultures during the last six hours and differential staining was revealed by staining chromosomes with 4% Giemsa solution in 2% 4Na-EDTA for two minutes (Odierna et al. 2004).
70%-majority consensus tree derived from a Bayesian inference analysis of 321 bp of the mitochondrial 16S rRNA gene. Ophisaurus attenuatus was used as outgroup. Sequences retrieved from GenBank are marked with their accession numbers.
https://binary.pensoft.net/fig/12354Results and discussion
The alignment of the analysed 16S rRNA gene fragments showed a minimum of 14 nucleotidic substitutions (11 transitions and 3 transversions) and 2 insertion/deletions between samples from Italy and Spain, corresponding to an average uncorrected genetic distance of 4.2%, thus confirming the genetic distinctiveness of the Italian populations. On the contrary, the analyzed Italian specimens were genetically very uniform and showed an intraspecific uncorrected divergence of 0.4%, whereas their uncorrected genetic distance from the outgroup, Ophisaurus attenuatus, was 12.8%. According to Gvoždik et al. (2010) the average genetic divergence at the ND2 mitochondrial locus among the several Anguis species is about 7%, while the maximum intraspecific distance within the various Anguis colchica clades is 4.4%. Unfortunately, we did not use the same molecular marker, but the observed 4.2% distance at the more conserved 16S rRNA gene fragment among Spanish and Italian samples provide further evidenced for the differentiation of these two lineages. Even if the phylogenetic analyses resulted in a tree (Fig. 1) with largely unresolved relationships, we could recover a good support for the monophyly of the Italian specimens, as well as for a clear mithocondrial segregation between them and the extra-Italian Anguis samples analysed in this study. Metaphase plates suitable for chromosome analysis were obtained from four out the six investigated populations (Travacò Siccomario, Valle Agricola, Ruvo del Monte and Monte Cocuzzo). All individuals showed a karyotype of 2n = 44 elements, comprising ten pairs of macrochromosomes and twelve pairs of microchromosomes. The first macrochromosomes pair was metacentric, the fourth and fifth were submetacentric and the others were telocentric with the exception of the tenth pair (Fig. 2A; Table 1). Chromosomes morphology did not differ between sexes and in any of the various Italian localities studied. The Ag-NOR banding revealed the presence of NOR loci on three microchromosome pairs in all examined specimens (Fig. 2B, arrows). This condition is considered an apomorphic character in Squamata, derived from a plesiomorphic single NOR locus (Porter et al. 1994, Aprea et al. 1996). The C-banding + Giemsa and the C-banding+fluorochromes colorations rrevealed small pericentromeric C-bands on almost all macrochromosomes and on three pairs of microchromosomes, namely those NOR bearing ones (Fig. 2C, D, E). Both G and LR-banding did not reveal any difference among analysed populations or between sexes (Fig. 3A, B). However, LR banding highlighted the presence of a uncondensed trait in three out of 20 examined metaphase plates (15% on the short arms of one chromosome of the tenth pair). This uncondensed trait was observed only on metaphase plates from cultures with BrdU addition and probably is a fragile site. In fact, it is known that BrdU promotes the expression of fragile sites (e.g. Sutherland 1979, Stone and Stephens 1993) (Fig. 3B). Fragile sites are usually due to a DNA strand breakage and represent hot spots for recombination (Glover and Stein 1988). The presence of a fragile site on the tenth chromosome suggests the propensity for this chromosome for structural rearrangements. By chance, this chromosome pair karyologically differentiates Italian slow worm populations from Centro-European ones, this pair being shaped as submetacentric in the former populations and as telocentric in the Centro-European ones. A pericentric inversion can account for the different morphology of the tenth chromosome pair. Due to the absence of karyological data on other congeneric species, the polarity of chromosome rearrangements cannot be unambiguously assessed. The role of chromosome rearrangements in speciation is heavily debated (White 1978, Odierna et al. 1996, Simard et al. 2009), as it is difficult to ascertain if chromosome changes occur after species diversification or if they are involved in the speciation process in itself, acting as proximate cause of diversification. Chromosome polymorphism as well as interspecific differences characterized by a single chromosome inversion are known in lizards (Odierna et al. 2004, Kupriyanova et al. 2008, Aprea et al. 1996). Thus, the difference in the chromosome 10th pair might have occurred during or after the diversification between Italian and central European lineages. It should be also taken into consideration that the different morphology of the 10th pair may be a consequence of a centromere repositioning, so far observed only among different species using comparative genome analysis (Montefalcone et al. 1999, Cardone et al. 2007).
A Giemsa stained karyotype of an individual from Valle Agricola B Ag-NOR-banding stained metaphase plate of an individual from Ruvo del Monte C, D, E metaphase plate of an individual from Valle Agricola stained with C-banding (C) and C-Banding+ CMA3 (D)+ DAPI (E). Arrows in (B) point at NOR loci.
https://binary.pensoft.net/fig/12355
Karyotypes of a male from Travacò Siccomario (A) and of female from Ruvo del Monte (B) stained with G-banding and replication pattern, respectively. The arrow in (B) points at a fragile site.
https://binary.pensoft.net/fig/12356
Chromosome relative length (RL), centromeric index (CI) and chromosome shape (CS) of studied samples of Anguis veronensis. The values of RL and CI are expressed as mean ± standard deviation. Chromosome morphology was measured according to Levan et al. (1964). For the microchromosomes 11–22 only their complessive RL value is provided.
chrom.
RL
CI
CS
1
15,2 ± 3,5
0,45 ± 3,6
metacentric
2
11,7 ± 2,7
0,10 ± 3,8
telocentric
3
9,7 ± 3,3
0,08 ± 4,4
telocentric
4
7,6 ± 3,1
0,07 ± 3,5
telocentric
5
7,4 ± 3,4
0,42 ± 3,6
metacentric
6
7,1 ± 3,1
0,07 ± 4,1
telocentric
7
6,2 ± 3,1
0,10 ± 3,1
telocentric
8
5,7 ± 2,8
0,08 ± 2,9
telocentric
9
5,1 ± 2,7
0,07 ± 4,0
telocentric
10
4,6 ± 3,2
0,34 ± 3,8
submetacentric
11–22
19,7 ± 6,8
–
–
To conclude, our karyological study is consistent and provides further support for the recently discovered molecular and morphological differentiation of the Italian slow worm lineage (Gvoždik et al. 2013), who ascribe the Italian populations to the recently resurrected Anguis veronensis.
Acknowledgements
Animal sampling was conducted under following authorizations: n. SCN/2D/2000/9213 of 01/06/2000 by Ministero dell’Ambiente of Italy. We are grateful to A. Gentilli, S. Scali and V. Caputo Barucchi for providing blood aliquots of the samples from Pavia and Ancona. We thank Prof. Rakesh Rastogi for his criticisms and for language revision. The work of Angelica Crottini was supported by Fundação para a Ciência e a Tecnologia SFRH/ BPD/72908/2010. We acknowledge the project “Genomics and Evolutionary Biology” cofinanced by North Portugal Regional Operational Programme 2007/2013 (ON.2 – O Novo Norte), under the National Strategic Reference Framework (NSRF), through the European Regional Development Fund (ERDF).
ReferencesAlbertEMSan MauroDGarcia-ParisMRuberLZardoyaR (2009) Effect of taxon sampling on recovering the phylogeny of squamate reptiles based on complete mitochondrial genome and nuclear gene sequence data.441(1–2): 12-21. doi: 10.1016/j.gene.2008.05.014ApreaGOdiernaGCapriglioneTCaputoVMorescalchiAOlmoE (1996) Heterochromatin and NOR distribution in the chromosomes of six gekkonid species of the genus Phelsuma (Squamata: Gekkonidae).110: 343-349.CardoneMFLomientoMTetiMGMisceoDRobertoRCapozziOD’AddabboPVenturaMRocchiMArchidiaconoN (2007) Evolutionary history of chromosome 11 featuring four distinct centromere repositioning events in Catarrhini.90: 35-43. doi: 10.1016/j.ygeno.2007.01.007CastoeTAJiangZJGuWWangZOPollockDD (2008) Adaptive evolution and functional redesign of core metabolic proteins in snakes.3(5): e2201. doi: 10.1371/journal.pone.0002201DalqA (1921) Étude de la spermatogenèse chez l’orvet, Anguis fragilis.31: 347-452.GigantinoRApreaGCapriglioneTGuarinoFMOdiernaG (2002) Caratteristiche genomiche dei vertebrati ectotermi del Parco del Matese. I. Risultati dell’analisi cromosomica. In: OdiernaGGuarinoFM (Eds) I vertebrati ectotermi del Parco Regionale del Matese.Napoli, 29-55.GloverTWSteinCK (1988) Chromosome breakage and recombination at fragile sites.43(3): 265-273.GvoždikVJandzikDLymberakisPJablonskiDMoravecJ (2010) Slow worm, Anguis fragilis (Reptilia: Anguidae) as a species complex: genetic structure reveals deep divergences.55: 460-472. doi: 10.1016/j.ympev.2010.01.007GvoždikVBenkovskyNCrottiniABellatiAMoravecJRomanoASacchiRJandzikD (2013) An ancient lineage of slow worms, genus Anguis (Squamata: Anguidae), survived in the Italian Peninsula. doi: 10.1016/j.ympev.2013.05.004HallTA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT., 41: 95–98.HowellWMBlackDA (1980) Controlled silver staining of nucleolus organizer regions with a protective colloidal developer: 1-step method.36: 1014-1015. doi: 10.1007/BF01953855KupriyanovaLKuksinAOdiernaG (2008) Karyotype, chromosome structure, reproductive modalities of three Southern Eurasian populations of the common lacertid lizard, Zootoca vivipara (Jacquin, 1787).3: 99-106.MargotA (1946) Demonstration de l’absence d’heterochromosomes morphologiquement differencies chez deux éspeces de Sauriens: Anguis fragilis L. et Lacerta vivipara Jacquin.53: 555-596.MontefalconeGTempestaSRocchiMArchidiaconoN (1999) Centromere repositioning.9: 1184-1188. doi: 10.1101/gr.9.12.1184OdiernaGApreaGCapriglioneTCaputoVOlmoE (1994) Chromosomal diversification in the gekkonid Tarentola mauritanica (Reptilia, Gekkonidae).61: 325–330. doi: 10.1080/11250009409355902OdiernaGApreaGArribasOJCapriglioneTCaputoVOlmoE (1996) The karyology of the Iberian rock lizards.52: 542–550.OdiernaGApreaGCapriglioneTPukyM (2004) Chromosomal evidence for the double origin of viviparity in the European common lizard, Lacerta (Zootoca) vivipara.14: 157-160.OdiernaGApreaGAndreoneFBöhemeWVencesM (2007) Cytosystematics of hyperoliid frogs: Phylogeny of Heterixalus, low karyotypic variability in hyperoliines and separate phylogenetic position of Leptopelis.74: 71-81. doi: 10.1080/11250000601088564PalumbiSR (1996) Nucleic acids II: the polymerase chain reaction. In: HillisDMMoritzCMableBK (Eds) Molecular Systematics, 2nd Edition. Sunderland, USA, 205-247.PorterCAHaidukMWQueirozde K (1994) Evolution and phylogenetic significance of ribosomal gene location in chromosomes of squamate reptiles.1994: 302-315. doi: 10.2307/1446980PosadaD (2008) jModelTest: Phylogenetic Model Averaging.25(7): 1253-1256. doi: 10.1093/molbev/msn083RambautADrummondAJ (2007) http://beast.bio.ed.ac.uk/TracerRonquistFHuelsenbeckJP (2003) MRBAYES: Bayesian phylogenetic inference under mixed models.19: 1572-1574. doi: 10.1093/bioinformatics/btg180SambrookJFritschEFManiatisT (1989) Molecular cloning: a laboratory manual.2nd edition. Cold Spring Harbor Lab Press, New York, 1626 pp.SimardFAyalaDKamdemGCPombiMEtounaJOseKFotsingJMFontenilleDBesanskyNJCostantiniC (2009) Ecological niche partitioning between Anopheles gambiae molecular forms in Cameroon: the ecological side of speciation.9: 17-45. doi: 10.1186/1472-6785-9-17SindacoRDoriaGRazzettiEBerniniF (2006) Atlante degli Anfibi e dei Rettili d’Italia/Atlas of Italian Amphibians and Reptiles., 792 pp.StoneDMStephensKE (1993) Bromodeoxyuridine induces chromosomal fragile sites in the canine genome.46: 198-202. doi: 10.1002/ajmg.1320460220SumnerAT (1972) A simple technique for demonstrating centromeric heterochromatin.75: 304-306. doi: 10.1016/0014-4827(72)90558-7SutherlandGR (1979) Factors affecting expression in lymphocyte cultures.31: 125-135.TamuraKPetersonDPetersonNStecherGNeiMKumarS (2011) MEGA5: Molecular Evolutionary Genetics Analysis using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods.28: 2731-2739. doi: 10.1093/molbev/msr121WhiteMJD (1978) Modes of Speciation.WH Freeman. San Francisco, 455 pp.