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.v12i1.2174421744Research ArticlePhalangiidaeGeneticsKaryosystematicsCenozoicSouthern AfricaMolecular technique reveals high variability of 18S rDNA distribution in harvestmen (Opiliones, Phalangiidae) from South AfricaŠťáhlavskýFrantišekstahlf@natur.cuni.cz1OpatovaVera2JustPavel1LotzLeon N.3HaddadCharles R.https://orcid.org/0000-0002-2317-77604Department of Zoology, Faculty of Science, Charles University, Viničná 7, CZ-12844 Praha, Czech RepublicCharles UniversityPrahaCzech RepublicDepartment of Biological Sciences and Auburn University Museum of Natural History, Auburn University, Auburn, AL 36849, USAAuburn UniversityAuburnUnited States of AmericaDepartment of Arachnology, National Museum, P.O. Box 266, Bloemfontein 9300, South AfricaNational MuseumBloemfonteinSouth AfricaDepartment of Zoology and Entomology, University of the Free State, P.O. Box 339, Bloemfontein 9300, South AfricaUniversity of the Free StateBloemfonteinSouth Africa
Corresponding author: František Šťáhlavský (stahlf@natur.cuni.cz)
Academic editor: M. Schneider
20181302201812141593D21320C-FFFB-1E7C-7A28-9309FFD79D6BBD37B31A-62D0-48B6-B7AD-FCFF7050FD0812225291810201727122017František Šťáhlavský, Vera Opatova, Pavel Just, Leon N. Lotz, Charles R. HaddadThis is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.http://zoobank.org/BD37B31A-62D0-48B6-B7AD-FCFF7050FD08
The knowledge of cytogenetics in the harvestmen family Phalangiidae has been based on taxa from the Northern Hemisphere. We performed cytogenetic analysis on Guruiaafricana (Karsch, 1878) (2n=24) and four species of the genus Rhampsinitus Simon, 1879 (2n=24, 26, 34) from South Africa. Fluorescence in situ hybridization with an 18S rDNA probe was used to analyze the number and the distribution of this cluster in the family Phalangiidae for the first time. The results support the cytogenetic characteristics typical for the majority of harvestmen taxa, i.e. the predominance of small biarmed chromosomes and the absence of morphologically well-differentiated sex chromosomes as an ancestral state. We identified the number of 18S rDNA sites ranging from two in R.qachasneki Kauri, 1962 to seven in one population of R.leighi Pocock, 1903. Moreover, we found differences in the number and localization of 18S rDNA sites in R.leighi between populations from two localities and between sexes of R.capensis (Loman, 1898). The heterozygous states of the 18S rDNA sites in these species may indicate the presence of XX/XY and ZZ/ZW sex chromosomes, and the possible existence of these systems in harvestmen is discussed. The variability of the 18S rDNA sites indicates intensive chromosomal changes during the differentiation of the karyotypes, which is in contrast to the usual uniformity in chromosomal morphology known from harvestmen so far.
Šťáhlavský F, Opatova V, Just P, Lotz LN, Haddad CR (2018) Molecular technique reveals high variability of 18S rDNA distribution in harvestmen (Opiliones, Phalangiidae) from South Africa. Comparative Cytogenetics 12(1): 41–59. https://doi.org/10.3897/CompCytogen.v12i1.21744
Introduction
The harvestmen (Opiliones) represent one of the oldest terrestrial arachnid lineages (Dunlop et al. 2004). They are traditionally divided into four suborders (Cyphophthalmi, Eupnoi, Dyspnoi and Laniatores) (e.g. Pinto-da-Rocha et al. 2007), altogether comprising more than 6500 described species (Kury 2017). Harvestmen are present on all continents (except Antarctica) (Giribet and Kury 2007) and have a wide range of ecological functions within the ecosystem (e.g. predators, detritivores, scavengers). Typically, harvestmen have limited dispersal capabilities and show tendencies to local endemicity (Giribet and Kury 2007). Despite the species richness and ancient geographical isolation of many groups, harvestmen possess similar cytogenetic characteristics within each suborder.
So far, the chromosomes of 90 species of harvestmen have been examined, which ranks them as the third best cytogenetically explored arachnid order (Tsurusaki et al. 2017). All groups usually have biarmed monocentric chromosomes (Tsurusaki 2007). The highest numbers of chromosomes are known from Laniatores (2n=25–109, median 81) (Schneider et al. 2008) and Cyphopthalmi (2n=24–52, median 30) (Svojanovská et al. 2016). Both Dyspnoi (2n=10–28, median 16) and Eupnoi (2n=10–36, median 22) typically possess a lower number of chromosomes (Tsurusaki et al. 2017). In some species, intraspecific variability in chromosome numbers has been reported. Given the low vagility of harvestmen, the variability could be caused by a fixation of local karyotype changes in isolated populations. Geographical variation among populations is particularly known from sabaconids (Dypnoi) and sclerosomatids (Eupnoi) (Tsurusaki 2007), where narrow hybrid zones exist between the populations with different chromosome races (e.g. Tsurusaki et al. 1991, Gorlov and Tsurusaki 2000a). However, the intraspecific variability in the chromosome number may also be attributed to the presence of B chromosomes, i.e. the supernumerary chromosomes that remain univalent during meiosis and are usually heterochromatic and non-coding. The number of B chromosomes can vary even within one individual (White 1973). In opilionids, the presence of B chromosomes has been documented in the sclerosomatid harvestman Psathyropustenuipes L. Koch, 1878 (Eupnoi), where some individuals possessed up to 18 B chromosomes (Tsurusaki 1993). Finally, the intraspecific variability could also be related to the presence of the heteromorphic XX/XY or ZZ/ZW sex chromosome systems. The XX/XY system evolved independently in both Dyspnoi and Eupnoi (Tsurusaki 2007), while the ZZ/ZW system is documented uniquely in Mitopusmorio (Fabricius, 1779) (Eupnoi: Phalangiidae) from Rishiri island in Japan (Tsurusaki and Cokendolpher 1990).
The majority of the available chromosome data on harvestmen originates from the suborder Eupnoi (58 species examined) (Tsurusaki et al. 2017). However, all the Eupnoi taxa examined so far, with exception of Holmbergianaweyenberghii (Holmberg, 1876) (Sclerosomatidae) from Argentina (Rodríguez Gil and Mola 2010), occur in the Northern Hemisphere. Despite our relatively good knowledge of the cytogenetics of this group, it is likely that the variability within Eupnoi could still be underestimated, due to the obvious geographical bias of the available data. The suborder Cyphophthalmi could serve as an example of a similar situation; it was traditionally perceived as cytogenetically uniform, but a surprising amount of diversity was revealed when analyses encompassed material from a broad geographical range (Svojanovská et al. 2016).
The suborder Eupnoi is divided into six families comprising 1822 species (Pinto-da-Rocha et al. 2007, Kury 2013). The family Phalangiidae (394 described species), with the centre of diversity in the Northern Hemisphere (Giribet and Kury 2007), represents the second most diverse group of Eupnoi (Kury 2017). The presence of Eupnoi in South Africa is traditionally explained by the group’s active dispersion from the Mediterranean region (Staręga 1984). Presently, only one species, Guruiaafricana (Karsch, 1878), and 30 species of the genus Rhampsinitus Simon, 1879 belonging to the family Phalangiidae are known from South Africa (Lotz 2009). In this study, we analyse the karyotypes of South African representatives of the family Phalangiidae and present the first data on the cytogenetics of this family of harvestmen from the Afrotropical Region. Moreover, we also use fluorescence in situ hybridization (FISH) to identify the number and position of the 18S rRNA genes, which is the first implementation of this method in the family Phalangiidae. The marker for 18S rRNA is frequently used to reveal concealed karyotype differences in groups with similar chromosomes (e.g. Cabral de Mello et al. 2011, Mattos et al. 2014, Sember et al. 2015) and to help inferring specific chromosomal changes along the course of the karyotype evolution (e.g. Nguyen et al. 2010, Grzywacz et al. 2011, Panzera et al. 2012).
Material and methods
The material examined in this study is deposited in the National Museum, Bloemfontein, South Africa (NMBA). We analyzed five species of harvestmen belonging to the family Phalangiidae from different localities in South Africa:
Rhampsinitusqachasneki Kauri, 1962: KwaZulu-Natal: Royal Natal National Park (28.7101°S, 28.9336°E), 1 male, 1 female (NMBAO01025, NMBAO01026).
The specimens were individually hand collected and kept alive until the gonad dissection. Chromosome preparation follows the “plate spreading” method (Traut 1976) widely used in arachnids (e.g. Šťáhlavský and Král 2004). After the dissection, the gonad tissue was hypotonized in 0.075 M KCl for 20 min, subsequently fixed in methanol: acetic acid (3:1) solution for 20 min, and finally dissociated in a drop of 60% acetic acid directly on the surface of the microscope slide. After the dissociation, the slide with the suspension was moved onto a histological plate (45 °C) to allow the evaporation of the liquid component of the suspension. The chromosomes were stained in a 5% Giemsa solution in Sörensen phosphate buffer for 20 min. Chromosomes were photographed using an ORCA-AG monochromatic camera (Hamamatsu) in an Olympus IX81 microscope. The karyotypes were analyzed from photographs using the LEVAN plugin (Sacamato and Zacaro 2009) for IMAGEJ 1.47 program (http://imagej.nih.gov/ij/). We used the classification of chromosomes of Levan et al. (1964) and measured the relative length of the chromosomes for the diploid set (Suppl. material 1).
FISH detection of 18S rDNA
The probe for 18S rDNA was prepared from Euscorpiussicanus (Koch, 1837). The whole genomic DNA was extracted using Genomic DNA Minikit - Tissue (Geneaid), following the manufacturer’s guidelines. The 18S rDNA fragment (GenBank accession number MG761815) were amplified with the following primer combination 18S-Gal forward: 5‘-CGAGCGCTTTTATTAGACCA-3‘ and 18S-Gal reverse: 5‘-GGTTCACCTACGGAAACCTT-3‘ (Fuková et al. 2005), and 50 ng DNA template. The PCR protocol was as follows: 95°C for 3 min, 30 cycles at 94°C for 30 s, 55 °C for 30 s and 72 °C for 2 min, final extension at 72 °C for 3 min. PCR products were purified with GenElute PCR Clean-Up Kit (Sigma-Aldrich). The probe was labelled using PCR with biotin-14-dUTP (Roche) using Nick Translation Kit (Abbott Molecular) following the manufacturer’s guidelines.
The FISH protocol was performed following Forman et al. (2013). Briefly, chromosome preparations were treated with RNase A (200 µg/ml in 2× SSC) for 60 min and then washed twice in 2× SSC for 5 min. Chromosomes were denaturized at 68 °C for 3 min 30 s in 70 % formamide in 2× SSC. The probe mixture, comprising 20 ng of the probe and 2.5 µl of salmon sperm DNA (Sigma-Aldrich) for each slide, was applied. The slides were left to hybridize overnight in a dark chamber at 37 °C. The following day, the preparations were treated with Cy3-conjugated streptavidin, followed by application of biotinylated antistreptavidin and another dose of Cy3-conjugated streptavidin. Chromosome preparations were counterstained with DAPI (4',6-diamidino-2-phenylindole), contained in FluoroshieldTM (Sigma-Aldrich) and observed in an Olympus IX81 microscope equipped with an ORCA-AG monochromatic charge-coupled device camera (Hamamatsu). The images were pseudocolored (red for Cy3 and blue for DAPI) and superimposed with Cell^R software (Olympus Soft Imaging Solutions GmbH).
ResultsGuruiaafricana (Karsch, 1878)
The diploid number of chromosomes in the four males analysed was 24 (Fig. 1a). The karyotype of this species comprised five pairs of metacentric (pairs No. 1, 3, 5, 6, 9), seven pairs of submetacentric and one pair of subtelocentric (pair No. 12) chromosomes (Suppl. material 1). The chromosomes gradually decreased in length from 6.49 % to 2.74 % of the diploid set (Suppl. material 1). In this species, we observed visible modification of condensation of the chromosomes, such as positive heteropycnosis during mitosis. In the chromosomes, there were visible positive heteropycnotic blocks at the position of the centromere of mitotic metaphases (pairs No. 2, 5, 7, 10) (Fig. 1a); however, in some chromosomes, the heteropycnotic parts expanded to a large proportion of the arms (pairs No. 1, 3, 4, 8, 11), or even across the whole length of the chromosome (pairs No. 6, 9). Large positive heteropycnotic blocks were visible also during prophase I of meiosis, mainly during zygotene and pachytene (Fig. 1b). During the diffuse stage and the diplotene the heteropycnosis was only moderate (Fig. 1c). Later, during diakinesis (Fig. 1d), metaphase I and metaphase II (Fig. 1e), the whole chromosomes seemed to be isopycnotic (Fig. 1d, e). We did not clearly detect visible heteromorphic bivalents during diakinesis and metaphase I (Fig. 1d).
66AAA55B-A5E2-557F-B7A6-3CB0CA988C41
Chromosomes of Guruiaafricana (2n=24) males after Giemsa staining (A–E) and after FISH with 18S rDNA (A partly, F–I). A karyotype based on mitotic metaphase B pachytene C diplotene D diakinesis E one sister cell of metaphase II F–I are the same cells as B–E after FISH with 18S rDNA. White arrows indicate the position of 18S rDNA and black arrows indicate the same position after Giemsa staining. Bar = 10 μm.
https://binary.pensoft.net/fig/185714
Two pairs of 18S rDNA clusters were detected by FISH in this species. The 18S rDNA probe signals were localized interstitially at approximately one third of the short arm of chromosome pair No. 1 and approximately at one third of the long arm of chromosome pair No. 3 (Fig. 1a). The signals on chromosome pair No. 3 were more intensive than on the longest chromosomes. During pachytene we identified the position of the 18S rDNA clusters on the positive heteropycnotic areas of the chromosomes, and the clusters were usually associated in this phase (Fig. 1f). Also during the diffuse stage and the diplotene the chromosomes bearing the NORs were close to each other (Fig. 1g). The condensation of the chromosomes caused the visualization of the 18S rDNA in the terminal part of the chromosomes during diakinesis and metaphase II (Fig. 1h, i).
Rhampsinituscapensis (Loman, 1898)
The diploid number of chromosomes in the two males and two females analysed was 34 (Fig. 2a, b). The karyotype of males comprised two pairs of metacentric (pairs No. 4, 16), three pairs of submetacentric (pairs No. 2, 5, 6, 14), eight pairs of subtelocentric, and four pairs of acrocentric (pairs No. 8, 11, 12, 15) chromosomes (Suppl. material 1). The first pair of chromosomes was slightly longer (5.49 % of the diploid set) than the following chromosome pair (4.50 %), and the rest of the chromosomes gradually decreased in length to 1.49 % of the diploid set (Suppl. material 1). In some chromosomes, positive heteropycnotic blocks were visible at the position of the centromere during mitotic metaphase in both sexes (Figs 2a, b), and also during pachytene in females (Fig. 2d) (we did not observe meiosis in males). During this phase we only observed one pair with the positive heteropycnotic areas expanded to the larger proportion of the arms (Fig. 2d). In pachytene of females we detected one heteromorphic bivalent (Fig. 2d). Three signals of 18S rDNA clusters were detected by FISH in males of this species (Fig. 2a). The 18S rDNA probe signals were localized in the terminal position of the short arms of chromosome pair No. 14 and at one chromosome of pair No. 16. In females, we detected five signals of 18S rDNA clusters. During the mitotic metaphase they seemed to be in a terminal position (Figs 2b, c); however, during the pachytene it was clear that one pair was in a subterminal position (Fig. 2e). Interestingly, one cluster of 18S rDNA was very large (Fig. 2c) and covered nearly the whole arm of the large metacentric chromosome from the heteromorphic bivalent (Fig. 2e).
5B3C268B-3648-5043-99DE-97FCA4899460
Chromosomes of Rhampsinituscapensis (2n=34) (A-E) and Rhampsinitusdiscolor (F–K) (2n=24) after Giemsa staining (A, B, D, F, G, J) and after FISH with 18S rDNA (A and F partly, C, E, H, I, K). A male karyotype based on mitotic metaphase B, C female mitotic metaphase D, E female pachytene, bottom left inset shows a reconstruction of the heteromorphic bivalent; crosses indicate the subterminal position of 18S rDNA; arrowhead indicates bivalent with the expanded positive heteropycnotic area F male karyotype based on two sister cells of metaphase II G, H male mitotic metaphase I diakinesis J, K two sister cells of metaphase II. White arrows indicate the position of 18S rDNA and black arrows indicate the position after Giemsa staining. Asterisks indicate large signal of heteromorphic bivalent. Bar = 10 μm.
The diploid number of chromosomes in the two males analysed was 24 (Fig. 2f). The karyotype comprised eight pairs of metacentric, two pairs of submetacentric (pairs No. 6, 9), and two pairs of subtelocentric (pairs No. 8, 12) chromosomes (Suppl. material 1). The chromosomes gradually decreased in length from 5.96 % to 2.69 % of the diploid set (Suppl. material 1). In this species, we observed modification in spiralization, visible as positive heteropycnosis of large parts of the arms of almost all chromosomes during mitotic prophase. These large blocks of positive heteropycnotic regions were visible only in a few chromosomes later during mitotic metaphase (Fig. 2g). Large positive heteropycnotic blocks were visible also during meiosis, from leptotene. These positive heteropycnotic parts usually associated during pachytene, and they were still clearly visible during diplotene. Some chromosomes still showed the positive heteropycnosis during diakinesis and metaphase II (Fig. 2j). During meiosis, we did not detect clearly visible heteromorphic bivalents (Fig. 2i). Two pairs of 18S rDNA clusters were detected by FISH in this species (Fig. 2h, i, k). The 18S rDNA probe signals were localized at the terminal parts of the first and third pairs of chromosomes (Fig. 2f).
Rhampsinitusleighi Pocock, 1903
We identified two cytotypes from two different localities in this species. The diploid numbers of chromosomes in both of them were 26 (Fig. 3a, b), but the morphology of the chromosomes and the number and position of 18S rDNA clusters differed.
The karyotype of cytotype I from Vernon Crookes comprised six pairs of metacentric, five pairs of submetacentric (pairs No. 2, 4, 8, 9, 13), and one pair of subtelocentric (pair No. 5) chromosomes. Moreover, the karyotype of cytotype I included one heteromorphic bivalent (pair No. 11) composed of one metacentric and one submetacentric chromosome. The chromosomes of the first two pairs were distinctly longer (6.18 % and 5.67 % of the diploid set) than the next chromosomes, which gradually decrease in length from 4.29 % to 2.42 % of the diploid set (Suppl. material 1). In some chromosomes, positive heteropycnotic blocks were visible at the position of the centromere in mitotic metaphases, and in some pairs heteropycnotic parts expanded to a large proportion of the arms (Fig. 3a). Large heteropycnotic blocks were also visible during the early meiotic phases such as zygotene, pachytene (Fig. 3c) and diplotene (Fig. 3e). Later, during metaphase I and metaphase II the whole chromosome seemed to be isopycnotic. We only found bivalents with one chiasma in both males from Vernon Crookes during the diplotene (Fig. 3e, Suppl. material 2). Altogether seven signals of 18S rDNA clusters were detected by FISH in cytotype I (Fig. 3f). All of them were localized in the terminal position of the chromosomes. Interestingly, we identified a heterozygous state in two pairs of chromosomes (Fig. 3a). We identified a small signal in the terminal position of one chromosome of pair No. 5 (Fig. 3a), which is also visible during pachytene (Fig. 3d) and diplotene (Fig. 3f). In pair No. 11 we identified an expanded part of one chromosome that covered 18S rDNA clusters, which makes the bivalent heteromorphic (compare Fig. 3c, d, e, f).
E78BB791-4AB6-546B-801F-D946C0812BB6
Chromosomes of Rhampsinitusleighi (2n=26) from Vernon Crookes (A, C–F) and from Silaka (B, G–M) after Giemsa staining (A–C, E, G, I, J, L) and after FISH with 18S rDNA (A and B partly, D, F, H, K, M). A, B karyotypes of males based on mitotic metaphase C, D pachytene with heterozygous bivalent (asterisk), bottom right inset shows a reconstruction of the heteromorphic bivalent E, F diplotene with heterozygous bivalent (asterisk) G, H mitotic prophase I diplotene, arrowheads indicate bivalents with two chiasmata. J, K diakinesis L, M one sister cell of metaphase II. White arrows indicate the position of 18S rDNA and black arrows indicate the same position after Giemsa staining. Asterisks indicate heteromorphic bivalent with large signal of rDNA. Bar = 10 μm.
https://binary.pensoft.net/fig/185716
The karyotype of cytotype II from Silaka (Fig. 3b) comprised seven pairs of metacentric and five pairs of submetacentric (pairs No. 3, 7, 9, 10, 12) chromosomes. Similar to cytotype I, cytotype II also included one heteromorphic bivalent (pair No. 11) composed of one metacentric and one submetacentric chromosome (Suppl. material 1). The chromosomes of the first three pairs were slightly longer (6.04 %, 5.66 %, and 5.10 % of the diploid set) and the following chromosomes gradually decreased in length from 4.42 % to 2.21 % of the diploid set (Suppl. material 1). Also in cytotype II, positive heteropycnotic blocks were visible at the position of the centromere of the chromosomes during mitotic prometaphase (Fig. 3g) and metaphase, and in some pairs the heteropycnotic parts expanded to a large portion of the arms (Fig. 3b). Large heteropycnotic blocks were also visible during the early meiotic phases until diplotene (Fig. 3i). Later, during diakinesis, metaphase I and metaphase II, the whole chromosome seemed to be isopycnotic (Figs 3j, l). We found from four to eight bivalents with two chiasmata in the male from Silaka during diplotene (Fig. 3i, Suppl. material 2). Six signals of 18S rDNA clusters were detected by FISH in cytotype II (Figs 3b, h, k). All of them were localized in the terminal position of the chromosomes (Figs 3b, h, k, m). We detected a heterozygous state in one pair of chromosomes (Fig. 3b). In pair No. 11 we identified an expanded part of one chromosome that covered 18S rDNA clusters, which makes the bivalent heteromorphic (compare Fig. 3j, k).
Rhampsinitusqachasneki Kauri, 1962
The diploid number of chromosomes in this species was 24 (Fig. 4a, b). The karyotype of the male comprised six pairs of metacentric, four pairs of submetacentric (pair No. 2, 6, 9, 11), and two pairs of subtelocentric (pair No. 8, 12) chromosomes (Suppl. material 1). The chromosomes gradually decreased in length from 5.74 % to 2.87 % of the diploid set (Suppl. material 1). In some chromosomes, positive heteropycnotic blocks were visible at the position of the centromere of mitotic metaphases, but these blocks did not expand across the whole arms of the chromosomes. Large heteropycnotic blocks were also visible during the early phases of meiosis such as zygotene, pachytene and diplotene; however, later during the metaphase I and metaphase II the entire chromosomes seemed to be isopycnotic (Fig. 4e). One pair of 18S rDNA clusters was detected by FISH in the terminal part of one pair of chromosomes in this species (pair No. 1) (Fig. 4a, c, d, f).
C88352B3-3C8E-57C7-9EC2-0280D754C65B
Chromosomes of Rhampsinitusqachasneki (2n=24) males after Giemsa staining (A, B, E) and after FISH with 18S rDNA (A partly, C, D, F). A male karyotype based on two sister cells of metaphase II B, C mitotic metaphase D diakinesis E, F two sister cells of metaphase II. White arrows indicate the position of 18S rDNA and black arrows indicate the same position after Giemsa staining. Bar = 10 μm.
https://binary.pensoft.net/fig/185717Discussion
Phalangiid harvestmen have a centre of diversity in the Northern Hemisphere (Giribet and Kury 2007), and their presence in South Africa is traditionally explained by dispersal from the Mediterranean region (Staręga 1984). In this study, we focus on the genera Guruia and Rhampsinitus that are endemic to the Afrotropical biogeographical region (Staręga 1984), in order to provide further information about the cytogenetics of the family Phalangiidae from previously understudied parts of the World.
The results of the cytogenetic analyses of Guruiaafricana (2n=24) and four species of the genus Rhampsinitus (2n=24, 26, 34) show that the African taxa share the basic karyotype characteristics with the Northern Hemisphere phalangiids analysed so far (13 species). However, detailed information about the karyotypes is only available in three species, whereas the records for the remaining species only comprise chromosome numbers (see Tsurusaki et al. 2017). In Phalangiidae the chromosome number ranges from 2n=16 in Oligolophustridens (C.L. Koch, 1836) (Tsurusaki 2007) to 2n=36 in Rilaenatriangularis (Herbst, 1799) (Sokolow 1930), with the karyotypes predominated by biarmed chromosomes. The most frequent number of chromosomes in this group is 2n=24 (8 species) and 2n=32 (5 species) (Tsurusaki et al. 2017, present study). The hypothesised ancestral number of chromosomes in the order Opiliones is around 30 (Svojanovská et al. 2016). The low numbers of the chromosomes in Eupnoi may be explained in two different ways. The first hypothesis assumes independent reductions of the proposed ancestral number of chromosomes shared by the whole order in each family of Phalangioidea. The family Phalangiidae still harbours some species with an ancestral number close to 30 (see Tsurusaki et al. 2017). Therefore the 2n=34, found in Rhampsinituscapensis (the highest number of chromosomes), would represent a karyotype derived by centric fissions, since one-armed chromosomes (subtelocentrics and acrocentrics) are still common in this species’ karyotype. The rest of Rhampsinitus species (2n=24 and 26) and Guruiaafricana (2n=24) would thus possess karyotypes derived by chromosome rearrangements, like fusions and inversions. This hypothesis would imply an independent reduction of the diploid number in the families Protolophidae (2n=18–22) and Sclerosomatidae (2n=10–24) (see Tsurusaki et al. 2017). The diploid number 2n=30 is also known from the family Caddidae (Tsurusaki and Cokendolpher 1990), a sister group to the superfamily Phalangioidea (Giribet et al. 2010), which would further indicate that the diploid number close to 30 probably represents the ancestral state in Eupnoi.
Alternatively, the second hypothesis presumes the ancestral number of chromosomes of the whole suborder Eupnoi being between 20 and 24 (the most frequent number in this suborder; see Tsurusaki et al. 2017). In this case, all higher numbers of chromosomes would represent a derived state driven by chromosomal fissions. However, the direction of karyotype evolution in Eupnoi can only be resolved by interpreting the cytogenetic data within a wider phylogenetic framework. Unfortunately, a robust phylogenetic scheme for Eupnoi is not currently available.
The data concerning the number and position of the 18S rDNA clusters using the FISH technique in arachnids are still limited (see Svojanovská et al. 2016), and most of them come from scorpions (e.g. Mattos et al. 2014, Sadílek et al. 2015). The results suggest that the position of 18S rDNA clusters can be highly conserved among different genera/species of scorpions, e.g. in a terminal position of the genera Rhopalurus Thorell, 1876 and Tityus C. L. Koch, 1836 (e.g. Mattos et al. 2014) from South America, or in an interstitial position in Androctonus Ehrenberg, 1828 from Africa and Asia (Sadílek et al. 2015).
Contrastingly, the variability of the number of 18S rDNA sites may be high in some arachnids (e.g. Forman et al. 2013). However, the interpretation of the differences in the number of 18S rDNA clusters requires a precise knowledge of both the ancestral number of these clusters and the phylogeny of the group, otherwise the course of the karyotype evolution cannot be specified (e.g. Svojanovská et al. 2016), as in case of our results. There is no information available directly concerning the 18S rDNA sites of phalangiid harvestmen in the literature that would allow a comparison to our results. The only record available for the entire suborder Eupnoi belongs to the sclerosomatid harvestman Psathyropustenuipes (L. Koch, 1878), which possesses one NOR (a transcription unit that encodes 28S, 5.8S and 18S rRNAs) in the diploid set (Gorlov and Tsurusaki 2000b). However, the NOR was visualised by silver staining, which only detects the NORs that are active during the preceding interphase (Miller et al. 1976), and thus the real number may be underestimated (Forman et al. 2013). Generally, it is hypothesized that one pair of NORs, as well as 18S rDNA sites, represents the ancestral state in the class Arachnida (Forman et al. 2013). We detected one pair of 18S rDNA clusters only in one species included in the analyses (Rhampsinitusqachasneki), while the number of clusters varied between 2 in R.qachasneki and 7 in R.leighi. The high variability in the number, position and even size of rDNA clusters between and also within the analysed species suggests intensive chromosomal changes and rapid evolutionary dynamics of 18S rDNA clusters on chromosomes of South African phalangiids. Rhampsinitus has a terminal position of 18S rDNA clusters, similarly to basal harvestmen of the suborder Cyphophthalmi (Svojanovská et al. 2016), and other arachnid groups such as spiders (e.g. Král et al. 2013), amblypygids (Paula-Neto et al. 2013), and South American scorpions (e.g. Mattos et al. 2014). The proximity to the telomeric region facilitates the efficiency of ectopic recombination (Goldman and Lichten 1996) and the preferential replacement of rDNA genes into a new subtelomeric position (Nguyen et al. 2010), as detected in all Rhampsinitus species. The effect of ectopic recombination may also be increased due to the association of positive heteropycnotic areas of different chromosomes during the meiotic prophase (e.g. Figs 1b, 3c). However, Guruiaafricana has 18S rDNA clusters placed in an interstitial position of one arm of biarmed chromosomes, which may be a consequence of another chromosome rearrangements, transposable element insertions or ectopic recombination (see Cabrero and Camacho 2008). Some 18S rDNA clusters in Rhampsinitus species were longer (see Fig. 3a, b), likely due to the presence of a higher number of copies of the 18S rDNA elements. This observation can be also explained by the duplication of the chromosome part including this cluster, or alternatively as a consequence of an insertion of transposable elements, which potentially accelerates the genomic reorganization after a speciation event (e.g. Symonová et al. 2013).
We also detected intraspecific variability and a heterozygous state in R.capensis and R.leighi that may correspond to population polymorphism. The low dispersal capability, typical for most harvestmen (Giribet and Kury 2007), can lead to a genetic differentiation among the populations on a relatively small geographical scale in many arachnid groups (e.g. Kotrbová et al. 2016, Opatova et al. 2016). The intraspecific variability in the number of 18S rDNA clusters between the populations from the two localities of R.leighi could thus be explained as an independent accumulation of the chromosomal changes at each locality, as is assumed for the wolf spider Wadicosafidelis (O. Pickard-Cambridge, 1872) (Forman et al. 2013). However, the harvestmen are also an understudied group, so the differences in 18S rDNA clusters and karyotypes between both populations could indicate a possible existence of cryptic species that are common in this order (e.g. Murienne et al. 2010). Given this scenario, the presence of a heterozygous state could be a result of crossbreeding.
Interestingly, we also identified heteromorphic bivalents in both populations of R.leighi caused by distinctly enlarged 18S rDNA clusters on one chromosome. Similarly, the heteromorphic bivalents were also detected in R.capensis, but the enlarged rDNA sites were only identified in females. The heterozygous states in the size of 18S rDNA clusters could indicate the existence of XY sex chromosomes in R.leighi and ZW sex chromosomes in R.capensis. The XY sex chromosome system is known in the harvestmen Sabaconmakinoi Suzuki, 1949 (Dyspnoi, Sabaconidae) (Tsurusaki 1989) and from 13 species of Sclerosomatidae (Eupnoi) (see Tsurusaki et al. 2017). In other arachnid orders the XY system is common in ticks (Oliver 1977), but it was also identified in few species from different spider families (Král et al. 2006, Řezáč et al. 2006), and some species of the pseudoscorpion family Neobisiidae (Troiano 1990, 1997, Šťáhlavský and Zaragoza 2008). The ZW system is completely exceptional in arachnids (Tsurusaki 2007). It is described only from Mitopusmorio and it is hypothesized in Olioglophusaspersus (Karsch, 1881) (both Phalangiidae) (Tsurusaki and Cokendolpher 1990). However, the ZW system in M.morio was described from a low quality single standard stained mitotic metaphase of one female (see Tsurusaki and Cokendolpher 1990). In this study, heteromorphic bivalents were detected in both females of R.capensis (Fig. 2b, d); however, we did not obtain any information about meiosis and chromosome pairing in males. The number and position of 18S rDNA clusters in mitotic metaphases of the analysed males does not indicate the presence of homomorphic chromosomes that would correspond to the heteromorphic state in females. Similarly, the evidence of a XY system in R.leighi is unclear, since we have no information concerning the organization of 18S rDNA clusters in females.
The size variation of rDNA clusters on homologous chromosomes is frequently considered a result of unequal crossing-over or gene duplication that does not represent sex-specific differences (e.g. Foresti et al. 1981). However, when the sex locus is located close to the rDNA, the recombination in this region can be limited and thus accelerate the evolution of the sex chromosome(s) (Reed and Phillips 1997). For that reason, we cannot rule out that the heteromorphic bivalents in R.leighi and R.capensis may represent the initial stages of sex differentiation from the ancestral nondifferentiated state typical for harvestmen (Tsurusaki 2007). Our analysis indicates that the differences in morphology and size of the chromosomes, considered in some harvestmen to be the sex chromosomes (e.g. Tsurusaki and Cokendolpher 1990), may in fact represent the variability in rDNA clusters induced by ectopic recombination, unequal crossing-over or gene duplication mentioned above. It would be necessary to analyse the species with heteromorphic sex chromosomes in order to know whether this heteromorphism is also associated with the rDNA. These findings bring a new light to our understanding of the presence of heteromorphic bivalents in harvestmen, and they highlight the necessity for the implementation of molecular cytogenetic techniques such as CGH in order to assess the presence of sex chromosomes in this group.
Acknowledgment
We are grateful to J.A. Neethling for technical support during the field work and to J. Štundlová for technical assistance during probe preparation. The project was supported by the Czech Republic Ministry of Education, Youth and Sports grant SVV 260 434 / 2017 and by grants from the National Research Foundation of South Africa in the Knowledge Interchange and Collaboration programme (No.’s 97495 and 105318) and Competitive Funding for Rated Researchers programme (No. 95569) to the last author. The specimens were collected under the collecting permits OP 4072/2016 from Ezemvelo KZN Wildlife (KwaZulu-Natal Province) and RA-0201 from the Eastern Cape Parks and Tourism Agency (Eastern Cape Province). Microscopy was performed in the Laboratory of Confocal and Fluorescence Microscopy, co-financed by the European Regional Development Fund and the state budget of the Czech Republic. Project no. CZ.1.05/4.1.00/16.0347 and CZ.2.16/3.1.00/21515.
ReferencesCabral de MelloDCMouraRCMartinsC (2011) Cytogenetic mapping of rRNAs and histone H3 genes in 14 species of Dichotomius (Coleoptera, Scarabaeidae, Scarabaeinae) beetles.134: 127–135. https://doi.org/10.1159/000326803CabreroJCamachoJPM (2008) Location and expression of ribosomal RNA genes in grasshoppers: abundance of silent and cryptic loci.16: 595–607. https://doi.org/10.1007/s10577-008-1214-xDunlopJAAndersonLIKerpHHassH (2004) A harvestman (Arachnida, Opiliones) from the Early Devonian Rhynie cherts, Aberdeenshire, Scotland. Transactions of the Royal Society of Edinburgh: Earth Sciences.94: 341–354. https://doi.org/10.1017/S0263593300000730ForestiFAlmeidaToledo LFToledoSA (1981) Polymorphic nature of nucleolus organizer regions in fishes.31: 137–144. https://doi.org/10.1159/000131639FormanMNguyenPHulaVKrálJ (2013) Sex chromosome pairing and extensive NOR polymorphism in Wadicosafidelis (Araneae: Lycosidae).141: 43–49. https://doi.org/10.1159/000351041FukováINguyenPMarecF (2005) Codling moth cytogenetics: karyotype, chromosomal location of rDNA, and molecular differentiation of sex chromosomes.48: 1083–1092. https://doi.org/10.1139/G05-063GiribetGKuryAB (2007) Phylogeny and Biogeography. In: Pinto-da-RochaRMachadoGGiribetG (Eds) Harvestmen: The Biology of Opiliones., 62–87.GiribetGVogtLGonzálezAPSharmaPKuryAB (2010) A multilocus approach to harvestman (Arachnida: Opiliones) phylogeny with emphasis on biogeography and the systematics of Laniatores.26: 408–437. https://doi.org/10.1111/j.1096-0031.2009.00296.xGoldmanASLichtenM (1996) The efficiency of meiotic recombination between dispersed sequences in Saccharomycescerevisiae depends upon their chromosomal location.144: 43–55. http://www.genetics.org/content/144/1/43.shortGorlovIPTsurusakiN (2000a) Staggered clines in a hybrid zone between two chromosome races of the harvestman Gagrellopsisnodulifera (Arachnida: Opiliones).54: 176–190. https://doi.org/10.1111/j.0014-3820.2000.tb00018.xGorlovIPTsurusakiN (2000b) Morphology and meiotic/mitotic behavior of B chromosomes in a Japanese harvestman, Metagagrellatenuipes (Arachnida: Opiliones): No evidence for B accumulation mechanisms.17: 349–355. https://doi.org/10.2108/jsz.17.349GrzywaczBMaryanska-NadachowskaAChobanovDPKaramyshevaTWarchalowska-SliwaE (2011) Comparative analysis of the location of rDNA in the Palaearctic bushcricket genus Isophya (Orthoptera: Tettigoniidae: Phaneropterinae).108: 509–517. https://doi.org/10.14411/eje.2011.066KotrbováJOpatovaVGardiniGŠťáhlavskýF (2016) Karyotype diversity of pseudoscorpions of the genus Chthonius (Pseudoscorpiones, Chthoniidae) in the Alps.10: 325–345. https://doi.org/10.3897/CompCytogen.v10i3.8906KrálJKořínkováTKrkavcováLMusilováJFormanMHerreraIMAHaddadCRVítkováMHenriquesSVargasJGPHedinM (2013) Evolution of karyotype, sex chromosomes, and meiosis in mygalomorph spiders (Araneae: Mygalomorphae).109: 377–408. https://doi.org/10.1111/bij.12056KrálJMusilováJŠťáhlavskýFŘezáčMAkanZEdwardsRLCoyleFAAlmerjeCR (2006) Evolution of the karyotype and sex chromosome systems in basal clades of araneomorph spiders.14: 859–880. https://doi.org/10.1007/s10577-006-1095-9KuryAB (2013) Order Opiliones Sundevall, 1833. In: ZhangZ-Q (Ed.) Animal Biodiversity: An Outline of Higher-level Classification and Survey of Taxonomic Richness (Addenda 2013).3703: 27–33. https://doi.org/10.11646/zootaxa.3703.1.7KuryAB (2017) Classification of Opiliones. Museu Nacional/UFRJ website. Online at: http://www.museunacional.ufrj.br/mndi/Aracnologia/opiliones.htmlLevanAFredgaKSandbergAA (1964) Nomenclature for centromeric position of chromosomes.52: 201–220. https://doi.org/10.1111/j.1601-5223.1964.tb01953.xLotzLN (2009) Harvestmen (Arachnida: Opiliones) in southern Africa - an annotated catalogue with notes on distribution.25: 1–45.MattosVFCarvalhoLSCellaDMSchneiderMC (2014) Location of 45S ribosomal genes in mitotic and meiotic chromosomes of buthid scorpions.31: 603–607. https://doi.org/10.2108/zs140005MillerDADewVTantravahiRMillerO (1976) Suppression of human nucleolar organizer activity in mouse human somatic hybrid cells.101: 235–243. https://doi.org/10.1016/0014-4827(76)90373-6MurienneJKaramanIGiribetG (2010) Explosive evolution of an ancient group of Cyphophthalmi (Arachnida: Opiliones) in the Balkan Peninsula.37: 90–102. https://doi.org/10.1111/j.1365-2699.2009.02180.xNguyenPSaharaKYoshidoAMarecF (2010) Evolutionary dynamics of rDNA clusters on chromosomes of moths and butterflies (Lepidoptera).138: 343–354. https://doi.org/10.1007/s10709-009-9424-5OliverJH (1977) Cytogenetics of mites and ticks.22: 407–429. https://doi.org/10.1146/annurev.en.22.010177.002203OpatovaVBondJEArnedoMA (2016) Uncovering the role of the Western Mediterranean tectonics in shaping the diversity and distribution of the trap-door spider genus Ummidia (Araneae, Ctenizidae).43: 1955–1966. https://doi.org/10.1111/jbi.12838PanzeraYPitaSFerreiroMJFerrandisILagesCPérezRSilvaAEGuerraMPanzeraF (2012) High dynamics of rDNA cluster location in kissing bug holocentric chromosomes (Triatominae, Heteroptera).138: 56–67. https://doi.org/10.1159/000341888Paula-NetoEAraujoDCarvalhoLSCellaDMSchneiderMC (2013) Chromosomal characteristics of a Brazilian whip spider (Amblypygi) and evolutionary relationships with other arachnid orders.12: 3726–3734. https://doi.org/10.4238/2013.September.19.3Pinto-da-RochaRGiribetG (2007) Taxonomy. In: Pinto-da-RochaRMachadoGGiribetG (Eds) Harvestmen: The Biology of Opiliones., 88–246.ReedKMPhillipsRB (1997) Polymorphism of the nucleolus organizer region (NOR) on the putative sex chromosomes of Arctic char (Salvelinusalpinus) is not sex related.5: 221–227. https://doi.org/10.1023/A:1018411417816Rodríguez-GilSGMolaLM (2010) Chromosome complement and meiosis of Holmbergianaweyenberghii (Opiliones: Sclerosomatidae: Gagrellinae) from Argentina.69: 164–170. https://www.biotaxa.org/RSEA/article/view/28588ŘezáčMKrálJMusilováJPekárS (2006) Unusual karyotype diversity in the European spiders of the genus Atypus (Araneae: Atypidae).143: 123–129. https://doi.org/10.1111/j.2006.0018-0661.01949.xSadílekDNguyenPKoçHKovaříkFYağmurEAŠťáhlavskýF (2015) Molecular cytogenetics of Androctonus scorpions: an oasis of calm in the turbulent karyotype evolution of the diverse family Buthidae.115: 69–75. https://doi.org/10.1159/000445863SakamotoYZacaroAA (2009) LEVAN, an ImageJ plugin for morphological cytogenetic analysis of mitotic and meiotic chromosomes. Initial version. http://rsbweb.nih.gov/ij/SchneiderMCZacaroAAPinto-da-RochaRCandidoDMCellaDM (2008) Conventional and ultrastructural analyses of the chromosomes of Discocyrtuspectinifemur (Opiliones, Laniatores, Gonyleptidae).47: 203–207. https://doi.org/10.1111/j.1439-0469.2008.00490.x.SemberABohlenJŠlechtováVAltmanováMSymonováRRábP (2015) Karyotype differentiation in 19 species of river loach fishes (Nemacheilidae, Teleostei): extensive variability associated with rDNA and heterochromatin distribution and its phylogenetic and ecological interpretation. BMC Evolutionary Biology 15: 251. https://doi.org/10.1186/s12862-015-0532-9SokolowI (1930) Untersuchungen über die Spermatogenese bei den Arachniden. IV. Über die spermatogenese der Phalangiden (Opiliones).10: 164–194. https://doi.org/10.1007/BF02450640StaręgaW (1984) Revision der Phalangiidae (Opiliones), III. Die afrikanischen Gattungen der Phalangiinae, nebst Katalog aller afrikanischen Arten der Familie.38: 1–79. http://rcin.org.pl/Content/58500/WA058_10835_P255-T38_Annal-Zool-Nr-1.pdfSvojanovskáHNguyenPHiřmanMTufIHWahabRAHaddadCRŠťáhlavskýF (2016) Karyotype evolution in harvestmen of the suborder Cyphophthalmi (Opiliones).148: 227–236. https://doi.org/10.1159/000445863SymonováRMajtánováZSemberAStaaksGBOBohlenJFreyhofJRábováMRábP (2013) Genome differentiation in a species pair of coregonine fishes: an extremely rapid speciation driven by stress-activated retrotransposons mediating extensive ribosomal DNA multiplications. BMC Evolutionary Biology 13: 42. https://doi.org/10.1186/1471-2148-13-42ŠťáhlavskýFKrálJ (2004) Karyotype analysis and achiasmatic meiosis in pseudoscorpions of the family Chthoniidae (Arachnida: Pseudoscorpiones).140: 49–60. https://doi.org/10.1111/j.1601–5223.2004.01783.xTrautW (1976) Pachytene mapping in the female silkworm Bombyxmori L. (Lepidoptera).58: 275–284. https://doi.org/10.1007/BF00292094TroianoG (1990) Karyotype and male meiosis of four species of Roncus L. Koch, 1873 (Pseudoscorpionida, Neobisiidae).57: 1–9. https://doi.org/10.1080/11250009009355666TroianoG (1997) Further studies on karyology of the pseudoscorpions of the gen. Roncus: the karyotype of Roncusgestroi and Roncusbelluatii.50: 271–279. https://doi.org/10.1080/00087114.1997.10797401TsurusakiN (1989) Geographic variation of chromosomes in Sabaconmakinoi Suzuki (Arachnida, Opiliones, Sabaconidae).44: 111–116.TsurusakiN (1993) Geographic variation of the number of B-chromosomes in Metagagrellatenuipes (Opiliones, Phalangiidae, Gagrellinae).33: 659–665.TsurusakiN (2007) Cytogenetics. In: Pinto-da-RochaRMachadoGGiribetG (Eds) Harvestmen: The Biology of Opiliones., 266–279.TsurusakiNCokendolpherJ (1990) Chromosomes of sixteen species of harvestmen (Arachnida, Opiliones, Caddidae and Phalangididae).18: 151–166. http://www.jstor.org/stable/3705833TsurusakiNMurakamiMShimokawaK (1991) Geographic variation of chromosomes in the Japanese harvestman, Gagrellopsisnodulifera with special reference to hybrid zone in western Honshu.8: 265–275.TsurusakiNSvojanovskáHSchöenhoferAŠťáhlavskýF (2017) The harvestmen cytogenetic database. Version 4.0. Online at: http://www.arthropodacytogenetics.bio.br/harvestmendatabase/index.htmlWhiteMJD (1973) Animal cytology and evolution, 3rd edn., Cambridge University Press.ZaragozaJAŠťáhlavskýF (2008) A new Roncus species (Pseudoscorpiones: Neobisiidae) from Montseny Natural Park (Catalonia, Spain), with remarks on karyology.1693: 27–40.Supplementary materials10.3897/compcytogen.v12i1.21744.suppl112261384B717422-13C0-5A50-A8FD-1AFFCA99C48D
Table S1
Microsoft Excel Worksheet (.xls).
Measurements of the relative chromosome length (% RCL) and arm ratio of South African harvestmen from family Phalangiidae (± standard deviation). Numbers of measured metaphases are given in brackets following the species names.
https://binary.pensoft.net/file/185718This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited.František Šťáhlavský, Vera Opatova, Pavel Just, Leon N. Lotz, Charles R. Haddad10.3897/compcytogen.v12i1.21744.suppl21226140E514BC49-A11A-5202-A1F5-5FDB0946314D
Table S2
Microsoft Excel Worksheet (.xls).
Numbers of nuclei with different numbers of chiasmata in diplotene of Rhampsinitusleighi (2n = 26) from Vernon Crookes (cytotype I) and Silaka (cytotype II).
https://binary.pensoft.net/file/185719This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited.František Šťáhlavský, Vera Opatova, Pavel Just, Leon N. Lotz, Charles R. Haddad