Is Radon Emission in Caves Causing Deletions in Satellite DNA Sequences of Cave-Dwelling Crickets?

Abstract

The most stable isotope of radon, 222Rn, represents the major source of natural radioactivity in confined environments such as mines, caves and houses. In this study, we explored the possible radon-related effects on the genome of Dolichopoda cave crickets (Orthoptera, Rhaphidophoridae) sampled in caves with different concentrations of radon. We analyzed specimens from ten populations belonging to two genetically closely related species, D. geniculata and D. laetitiae, and explored the possible association between the radioactivity dose and the level of genetic polymorphism in a specific family of satellite DNA (pDo500 satDNA). Radon concentration in the analyzed caves ranged from 221 to 26000 Bq/m3. Specimens coming from caves with the highest radon concentration showed also the highest variability estimates in both species, and the increased sequence heterogeneity at pDo500 satDNA level can be explained as an effect of the mutation pressure induced by radon in cave. We discovered a specific category of nuclear DNA, the highly repetitive satellite DNA, where the effects of the exposure at high levels of radon-related ionizing radiation are detectable, suggesting that the satDNA sequences might be a valuable tool to disclose harmful effects also in other organisms exposed to high levels of radon concentration.

About the Authors:

Giuliana Allegrucci

* E-mail: [email protected]

Affiliation: Department of Biology, University of Rome Tor Vergata, Rome, Italy

Valerio Sbordoni

Affiliation: Department of Biology, University of Rome Tor Vergata, Rome, Italy

Donatella Cesaroni

Affiliation: Department of Biology, University of Rome Tor Vergata, Rome, Italy

Introduction

Radon is a radioactive gas occurring naturally. It is part of the normal radioactive chain of uranium and represents the decay product of radium. It is a rare gas and usually migrates freely through faults and fragmented soils and may accumulate in caves and / or water. The most stable isotope of radon, 222Rn, has a half-life of about 4 days and due to this characteristic, its concentration decreases with increasing distance from the production area. Ground water has generally higher concentrations of 222Rn than surface water because the radon is continuously produced by the radium present in the rocks. 222Rn can be significantly high in hot sulfur spring waters [1]. Due to these characteristics, 222Rn represents the major source of natural radioactivity in confined environments such as mines, caves and houses. Typical domestic exposures are about 100 Becquerel per cubic meter (Bq/m3) indoors and 10–20 Bq/m3 outdoors [2]. Concentration limits of radon for domestic areas are variable and depend on the organization; the European Union established two threshold values, one for the old houses (400 Bq/m3) and one for the new ones (200 Bq/m3), while the US-EPA (2007) put the limit at concentration of 74 Bq/m3. Studies have demonstrated a significant and dose-related excess of lung cancer in radon-exposed miners (National Research Council 1988) and several ecologic studies have found increased rates of leukaemia in regions with elevated levels of radon in homes [3,4, 5, 6,7].

In caves, radon concentration is known to vary within an extremely wide range [8, 9]. Natural caves of volcanic origin can be characterized by exceedingly high levels of radon because of the presence of uranium and therefore of the decay chain products of uranium series [10, 11]. Artificial caves as cellars, Etruscan graves, and Roman cisterns are often built with tuff, a type of rock consisting of consolidated volcanic ash ejected during a volcanic eruption. In such environments, radon concentration may be very high.

The occurrence of a wide spectrum of radon concentration in Italian caves, and the possibility to find some of these caves constantly inhabited by Dolichopoda cave crickets (Orthoptera, Rhaphidophoridae), led us to evaluate these insects as a suitable model to study the effects of radon on cave life.

Dolichopoda cave crickets are strictly dependent upon caves and several populations inhabit cave-like habitats, such as rock crevices and ravines, cellars, catacombs, aqueducts, Etruscan tombs and other similar man-made hypogean environments. They have long been studied in our laboratory from a wide array of genetic and ecological aspects addressed to understand their evolution and phylogeny [12, 13, 14, 15, 16, 17, 18, 19, 20, 21]. A preliminary study, carried out through the Comet assay, suggested a statistically significant dose-effect increase of DNA damage in specimens of Dolichopoda from radon-polluted caves, especially for the brain cells [22].

Dolichopoda populations and species have also been investigated for processes of molecular evolution of satellite DNA (satDNA), [23, 24, 25]. SatDNA is a class of non-coding DNA typically organized in large homogeneous arrays of tandemly arranged repetition units. These units are usually located in the heterochromatic parts of the chromosomes in the regions close to the centromeres and telomeres. Repeat size can vary largely within and between species from only a few base pairs up to several thousand base pairs [25 and references therein]. Three specific satDNA families have been characterized for Dolichopoda species, two of them being species-specific (pDo102 and pDsPv400) and one (pDo500) occurring in all Dolichopoda species [23, 24]. A potential hammerhead (HH) ribozyme is embedded within the pDo500 tandemly repeat satDNA [26, 27].

In the present study, we explored whether increasing level of radon-related ionizing radiation could induce an increased mutation pressure at the level of the pDo500 tandemly repeat satDNA. We considered two species, D. laetitiae (Menozzi 1920) and D. geniculata (Costa, 1860) that, as demonstrated in previous studies [13, 14, 15, 18, 19], are genetically closely related. Population samples were collected in caves showing various amount of radon concentration, in order to investigate the possible association between radon-related ionizing radiation and the level of polymorphism in the pDo500 tandemly repeat satDNA

Materials and Methods

This study was formerly approved by Regione Lazio in Italy (Dipartimento del Territorio), as a contribution to the knowledge of radon effects on insects constantly subjected to radioactivity. None of the field surveys in the present study involved endangered or protected species and no permission was necessary for the studied areas. Specimens for the DNA analysis were collected in seven different caves showing different radon concentration (Table 1, Fig 1). Sequences of pDo500 tandemly repeats satDNA were derived from samples belonging to two populations of D. laetitiae and to five populations of D. geniculata. We also retrieved pDo500 satDNA sequences from GenBank for population samples coming from other three caves whose radon’s concentration values were known: two population samples of D. geniculata, [25, 28] and one of D. laetitiae [25, 29]. See Table 1 for details.

[Figure omitted. See PDF.]

Fig 1. Sampling sites of Dolichopoda populations considered in this study (see also Table 1).

Different colors refer to different species: pink to D. laetitiae and blue to D. geniculata.

https://doi.org/10.1371/journal.pone.0122456.g001

[Figure omitted. See PDF.]

Table 1. Dolichopoda population samples included in this study.

https://doi.org/10.1371/journal.pone.0122456.t001

Radon measurements

The presence of radon in caves was detected by the Alfa track detector LR115. This detector has a particular film capable of measuring 222Rn concentration. Its working is based on the principle that the radon's alpha particles leave traces on a film coated with a thin layer of gelatin. It has to be placed in a stable and dry location for an adequate time. In this study, the Alfa track detector LR115 was located at the center of the cave in all considered sites except for the Pastena cave (PAS). The latter is a large cave subdivided in two distinct rooms and the specialized equipment was located at the center of each room. Following the manufacturer instructions, the Alfa track detector was left in the caves for a month, to obtain the measure of the average concentration of radon based on traces left by alpha particles (certificate numbers from 13843 to 13849, in accordance with U.S EPA National Radon Proficiency Program EPA—CFA Recommended Test Report Format).

Laboratory procedures

Genomic DNA was extracted from leg muscles using the Sigma-Aldrich GenElute Mammalian genomic DNA Miniprep Kit, following the instructions.

PDo500 satDNA sequences were amplified with the following primers, 5'-GTTTTACACGTTCACTGCAG-3' and 5' GACACATTGATGAGACTGCAG-3' [24]. The PCR conditions were as follow: 95°C for 3 minutes, followed by 30 cycles of denaturation at 95°C for 30 seconds, annealing at 50°C for 30 seconds, elongation at 72°C for 30 seconds and one final elongation step at 72°C for 2 minutes. The obtained PCR products were cloned using the pGEM-T Easy Vector kit (Promega). Positive clones were selected through PCR amplification using the reverse and forward M13 primers. The obtained PCR products were purified using the enzymatic digestion (ExoSAP-IT, Affymetrix, U.K.) and sequenced using the ABI-3730 Genetic Analyzer. Alignment was carried out using Clustal X 1.81 [30].

Data analyses

SatDNA repeat polymorphism, considering the estimates of nucleotide diversity (π) and the average number of nucleotide differences (K), was investigated using DNAsp software [31] for each sampled population. DNAsp was also used to perform sliding window analysis in order to detect regions of high sequence conservation. We carried out this analysis by considering sites with alignment gap in the window length. The window size was set to 30 with step size of 5. The analysis was performed on both the complete alignment for each population and on consensus sequences for each species.

Insertion-deletion polymorphism was also analyzed, using the multiallelic option in DNAsp. The total number of indel events, the average indel length per event, the number of indel haplotypes and the indel haplotype diversity were calculated for each population.

To investigate possible relations among the radon concentration in caves, the amount of satDNA polymorphism, and possibly the taxonomic status of each population, two type of multivariate analyses were carried out. In particular, multivariate ordination of Dolichopoda population samples based on polymorphism’s measures was studied by Factorial Correspondence Analysis (FCA), [32], using xlstat 2014. The radon concentration in caves and the taxonomic status of each population were considered as supplementary variables. In this way, these two variables were not taken into account for the computation of the representation space and their coordinates were computed a posteriori.

A multivariate multiple regression analysis (manova) was carried out using past software [33] to compare radon concentration in caves and the taxonomic status of each population with the measures of satDNA polymorphism. In particular, the frequency of polymorphic sites (PSF) and indel sites, the average length of indel and the haplotype diversity per population were log transformed and considered in this analysis. In order to exclude the possibility that our results were influenced by the unbalanced sample size, we carried out a Linear Mixed Model, using xlstat 2014, considering the environmental radioactivity as explanatory variable and the length of indel, calculated for each individual, as the response variable. Radioactivity measures (Bq/m3) were log transformed and considered as fixed effect. To test the non-independency of data from the same cave, each cave population was considered as a random effect. Likelihood ratio tests (LRT) were computed using the restricted maximum likelihood (REML) method, as implemented in xlstat 2014. The significance of the fixed effect was tested by a LRT between the full model and a null model comprising only the intercept and the random effects.

Results

Radon measurements

Table 1 shows the mean values of natural radioactivity for each of ten measured caves. The highest levels of radioactivity were reported in MTR1 and PTV caves, showing values of 25,997 and 13,200 Bq/m3, respectively. All other caves, except for CLP and PSC, showed 222Rn concentration measures ranging between 982 and 2,700 Bq/m3, which are much higher than the lowest threshold value established by European Union for the radon concentration in dwellings (400 Bq/m3). Only CLP and PSC caves exhibited radon concentration levels close to the threshold established for human housing.

Satellite DNA

We analyzed 163 satDNA repeats of the pDo500 family from seven sampled populations, each represented with 4–95 sequences. The length of the pDo500 sequences ranged between 458 and 481 bp. The total alignment consisted of 500 positions. The average nucleotide composition was T = 35.2%, C = 23.9%, A = 21.6% and G = 19.3%. The estimated transition/transversion bias (R) was 0.79.

Results from sliding window analyses, carried out separately for each population and each species, indicated that population samples coming from caves with the highest radioactivity showed also the highest nucleotide diversity in the peaks of local maxima, regardless of the species (MTR1 in D.laetitiae and PTV in D.geniculata, Fig 2A and 2B). Results from sliding window analyses carried out on consensus sequences for each species (Fig 2C) indicated higher nucleotide diversity in D. geniculata (π ranged from 0 to 0.1) than in D. laetitiae (π ranged from 0 to 0.07).

[Figure omitted. See PDF.]

Fig 2. Sliding window analyses of the pDo500 satellite DNA in populations of D. laetitiae and D. geniculata included in this study.

The analysis was performed for each population of D. laetitiae (A), D. geniculata (B) and on consensus sequences for each species (C). The value of nucleotide diversity (π) was obtained by a sliding window size of 30 with step size 5.

https://doi.org/10.1371/journal.pone.0122456.g002

Table 2 shows the estimates of pDo500 satDNA polymorphism within population samples. Frequency of polymorphic sites (PSF), frequency of indel sites (SIF), and average length of indel (ALI), are remarkably higher in population samples coming from the three caves, with the highest radioactivity measures: MTR1 and MTR2 housing D. laetitiae and PTV housing D. geniculata. One specimen from MTR1, five from MTR2 and five from PTV showed also clones with very short sequences ranging from 102 to 378 bp. These short sequences were characterized by a large gap in the middle of the pDo500 satDNA sequence. In particular, deletion covered the region between the first 10 bases of pDo500 and the position 296 in the alignment and did not overlap the region of the potential hammerhead (HH) ribozyme embedded within the pDo500 satDNA [25; 26]. These short sequences were taken into account only when indel polymorphism analysis was carried out.

[Figure omitted. See PDF.]

Table 2. Polymorphism estimates in sampled populations of Dolichopoda cave crickets.

https://doi.org/10.1371/journal.pone.0122456.t002

Multiple regression analysis (manova) revealed a significant correlation between the polymorphism estimates and the radon concentration in cave (F = 4.02, Wilk’s lambda = 0.040; P = 0.033). In particular, the frequency of indel sites (SIF, p = 0.008) and the average length of indel (ALI, p = 0.004) were statistically significant correlated with the levels of radioactivity in cave. On the other hand, the nucleotide diversity (K, p = 0.057) and the average number of nucleotide differences (π, p = 0.078) showed a high tendency to be dependent from the taxonomic status of each population.

Results obtained from LMM analysis, carried out for each individual, showed a significant regression line (P = 0.05; Table 3), indicating an increase of the length of indel related to the radioactivity levels measured in the caves.

Fig 3 reports results from FCA with the first two axes explaining together 95.12% of the total variance. The first axis clearly separates population samples subjected to high radioactivity from all the others.

[Figure omitted. See PDF.]

Fig 3. Factorial Correspondence Analysis on polymorphism measures of Dolichopoda populations analyzed in this study.

Ordination of populations on the plane is described by the first two axes. Population samples coming from caves with high radon concentration are displaced into different portions of the ordination plane.

https://doi.org/10.1371/journal.pone.0122456.g003

[Figure omitted. See PDF.]

Table 3. Estimated values, standard errors (SE), t-values and significance (P-values) obtained by Likelihood Ratio Test (LRT) of full vs. reduced models, for the coefficients of fixed effects in Linear Mixed Models (LMMs) for the average length of indel (ALI, see text for details).

https://doi.org/10.1371/journal.pone.0122456.t003

Discussion

In recent years, the impact of chemicals and physical pollutants on the functionality of DNA has been investigated in many animal species [34]. Most of the studies evaluated the biological response to the agents considering gene mutation, chromosome aberration, sister chromatid exchanges, DNA damage by Comet assay, micronuclei [22, 35, 36, 37]. In this study, the possible biological response to the environmental radioactivity was investigated by considering a specific category of nuclear DNA, the satDNA, and two genetically closely related species. We used two species, in order to verify if, regardless of the specific genetic variability, they could have the same biological response to the environmental contaminant.

The analysis of variability showed that D. geniculata is more polymorphic than D. laetitiae, as expected [25]; (Table 2, Fig 2). However, we found significant correlations between some polymorphism estimates and radon concentration in caves in both species, regardless of the degree of variability expressed by each one. In particular, both manova and FCA analyses (Fig 3) revealed that the indel polymorphism (SIF, ALI; Table 2) is significantly correlated with radon concentration (Bq/m3; Table 1); while the nucleotide diversity and the average number of nucleotide differences appear to be species dependent. In particular, the two localities MTR1, hosting D. laetitiae, and PTV, hosting D. geniculata, showed 222Rn levels higher of one to two orders of magnitude than the other caves (Table 1). This very high radioactivity can be explained by the presence of sulfur springs in MTR1 and by both the very low circulation of air and the presence of numerous faults and fractures in PTV. Samples coming from these two sites showed an average length of indels (ALI; Table 2) greater than that observed in other caves. High levels of pDo500 polymorphism (SIF, ALI, PSF; Table 2) were detected also in Dolichopoda samples from MTR2 site, although the latter showed a radioactivity level of one order of magnitude lower than in MTR1 and PTV caves (Table 1). This result might be explained by assuming that samples from MTR1 and MTR2 belong to the same population. Indeed, these localities are very close and are located in the Natural Reserve of Monterano, MTR1 is an old sulfur mine, MTR2 is a cellar in an old ruined house built with tuff and is located at 100 meters away in front of MTR1, being separated by the Mignone river. The two caves are surrounded by woods, an optimal environment for Dolichopoda that, at night, exit the cave to forage and move around. Results from mtDNA Cytochrome Oxidase I sequences showed that individuals coming from the one or the other site are very similar, showing at most a singleton difference (not published data) and Dolichopoda cave crickets show generally a certain degree of mtDNA variability both between and within populations [19, 20, 21]. Therefore, it is reasonable to consider these two sites as hosting a single population and to expect that individuals of Dolichopoda transfer from one cave to another, being definitely subjected to the same radon dose.

The LMM analysis, carried out for each individual to attempt to correct for our unequal sample size, confirmed these results (Table 3), suggesting, again, that the radioactivity levels measured in the caves appear to be responsible for the gaps’ length observed in our samples.

The biological significance of satDNA has been the object of several discussion and generally, based on the diversity of satDNA in nucleotide sequences, length of repeats, genomic abundance, a specific function has not been yet assigned to this genomic region. However, a number of possible functions have been hypothesized [38] and most of them are related to heterochromatin and/or centromere formation and function.

Samples from populations in hypogean environments with the highest radioactivity showed also the highest frequency of indel sites and clones with pDo500 repeat sequences shorter than the standard, but that are integer at level of HH ribozyme region. Previous studies suggested that the HH region of the pDo500 sequence family has a functional role in Dolichopoda cave crickets although its function is unclear and remains to be investigated [27]. Our data seem to support the hypothesis that the HH ribozyme region could have an important role, since this region is never affected by the events of insertion / deletion.

Finally, our results are consistent with those from [22] where six out of the ten populations here studied were analyzed for DNA primary damage through Comet assay. Both haemocytes and brain cells taken from individuals from radon-polluted caves were tested and compared to a control group of cave crickets reared in absence of radon. Results indicated a statistically significant dose-effect increase of DNA damage in all caves, especially for the brain cells. In conclusion, we can infer from present data that the increased sequence heterogeneity at pDo500 satDNA level can be explained as an effect of the mutation pressure induced by radon in cave. Furthermore, we discovered a specific category of nuclear DNA, the highly repetitive satDNA sequences, where the effects of the exposure at high levels of radon-related ionizing radiation are detectable. Future researches could be addressed to evaluate and investigate if satDNA might be a valuable tool to reveal the effects of radon in other organisms.

Acknowledgments

We are indebted to Mauro Rampini and Giorgio Pintus who helped one of us (GA) in collecting samples in the different caves. Sara Cocchi helped us in the lab working. We also thank Gabriele Gentile for his suggestions in cloning techniques and Paolo Gratton for his help in LMM statistical analysis. Stefano De Felici helped us in the implementation of Fig 1.

Author Contributions

Conceived and designed the experiments: GA VS DC. Performed the experiments: GA. Analyzed the data: GA. Contributed reagents/materials/analysis tools: GA. Wrote the paper: GA DC VS.

Citation: Allegrucci G, Sbordoni V, Cesaroni D (2015) Is Radon Emission in Caves Causing Deletions in Satellite DNA Sequences of Cave-Dwelling Crickets? PLoS ONE 10(3): e0122456. https://doi.org/10.1371/journal.pone.0122456

References

1. Iyengar MAR. The environmental behavior of Radium. Tech Rep series, 1990: 310

2. Sperrin M, Gillmore G, Denman T Radon concentration variations in a Mendip cave cluster. Environ Manage Health 2001; 12: 476–482.

3. Albering HJ, Hageman G J, Kleinjans JCS, Engelen JJM, Koulishcer L, Herens C. Indoor radon exposure and cytogenetic damage. Lancet 1992; 340: 739. pmid:1355848

4. Khan MA, Cross FT, Buschbom RL, Brooks AL. Inhaled radon-induced genotoxicity in Wistar rat, Syrian Hamster, and Chinese hamster deep-lung fibroblasts in vivo. Mutat Res. 1995; 334: 131–137. pmid:7885364

5. Jostes RF. Genetic, cytogenetic, and carcinogenic effects of radon: a review. Mutat Res. 1996; 340: 125–139. pmid:8692177

6. Bilban M, Vaupotic J. Chromosome aberrations study of pupils in high radon level elementary school. Health Phys. 2001; 80: 157–163. pmid:11197464

7. Cristaldi M, Ieradi LA, Udroiu I, Zilli R. Comparative evaluation of background micronucleous frequencies in domestic mammals. Mutat Res. 2004; 559: 1–9. pmid:15066568

8. Hakl J, Hunyadi I, Varhegyi A. Radon monitoring in caves. In: Durrani SA, Ilic R, editors. Radon measurements by etched track detectors. World Scientific; 1997. pp. 261–283.

9. Cigna AA.The distribution of radon concentration in caves. Int J Speleol. 2003; 32: 113–115.

10. Baubron JC, Allard P, Sabroux JC, Tedesco D, Toutain JP. Soil gas emanations as precursory indicators of volcanic eruptions. J geol Soc Lond. 1991; 148: 571–576.

11. Aytekin H, Baldik R, Celebi N, Ataksor B, Tasdelen M, Kopuz G. Radon measurements in the caves of Zonguldak (Turkey). Radiat Prot Dosimetry 2006; 118: 117–121. pmid:16120690

12. Sbordoni V, De Matthaeis E, Cobolli Sbordoni M. Phosphoglucomutase polymorfism and natural selection in populations of the cave cricket Dolichopoda geniculata. J. Zoolog. Sist. Evol. Res. 1976; 14: 292–299.

13. Sbordoni V, Allegrucci G, Cesaroni D, Sbordoni Cobolli M, De Matthaeis E. Genetic structure of populations and species of Dolichopoda cave crickets: Evidence of peripatric divergence. Boll. Zool. 1985; 52: 139–156.

14. Sbordoni V, Allegrucci G, Caccone A, Carchini G, Cesaroni D. Microevolutionary studies in Dolichopodinae cave crickets. In: Baccetti B editor. Evolutionary Biology of Orthopteroid Insects. Horwood Ltd. Publ, Chichester, U.K. 1987. pp. 514–540

15. Allegrucci G, Cesaroni D, Sbordoni V. Adaptation and speciation of Dolichopoda cave crickets (Orth. Rhaph.): geographic variation of morphometric indices and allozyme frequencies. Biol J Linn Soc Lond. 1987; 31: 151–160.

16. De Pasquale L, Cesaroni D, Di Russo C, Sbordoni V. Trophic niche, age structure and seasonality in Dolichopoda cave crickets. Ecography 1995; 18: 217–224.

17. Allegrucci G, Minasi MG, Sbordoni V. Patterns of gene flow and genetic structure in cave-dwelling crickets of the Tuscan endemic, Dolichopoda schiavazzii (Orthoptera, Rhaphidophoridae). Heredity 1997; 78: 665–673.

18. Cesaroni D, Matarazzo P, Allegrucci G, Sbordoni V. Comparing patterns of geographic variation in cave crickets by combining geostatistic methods and Mantel tests. J Biogeogr. 1997; 24: 419–431.

19. Allegrucci G, Todisco V, Sbordoni V. Molecular phylogeography of Dolichopoda cave crickets (Orthoptera, Rhaphidophoridae): a scenario suggested by mitochondrial DNA. Mol Phylogenet Evol. 2005; 37: 153–164. pmid:15964214

20. Allegrucci G, Rampini M, Gratton P, Todisco V, Sbordoni V. Testing phylogenetic hypotheses for reconstructing the evolutionary history of Dolichopoda cave crickets in the eastern Mediterranean. J Biogeogr. 2009; 36: 1785–1797.

21. Allegrucci G, Trucchi E, Sbordoni V. Tempo and mode of species diversification in Dolichopoda cave crickets (Orthoptera, Rhaphidophoridae). Mol Phylogenet Evol. 2011; 60: 108–121. pmid:21514393

22. Gustavino B, Meschini R, Franzetti G, Gratton P, Allegrucci G, Sbordoni V. Genotoxicity testing for radon exposure: Dolichopoda (Orthoptera, Rhaphidophoridae) as potential bio-indicator of confined environments. Curr Zool. 2014; 60: 299–307.

23. Bachmann L, Venanzetti F, Sbordoni V. Characterization of a species-specific satellite DNA family of Dolichopoda schiavazzii (Orthoptera, Rhaphidophoridae) cave crickets. J Mol Evol. 1994; 39: 274–281. pmid:7932789

24. Bachmann L, Venanzetti F, Sbordoni V. Tandemly repeated satellite DNA of D. schiavazzii: a test for models on the evolution of highly repetitive DNA. J Mol Evol. 1996; 43:135–144. pmid:8660438

25. Martinsen L, Venanzetti F, Johnsen A, Sbordoni V, Bachmann L. Molecular evolution of the pDo500 family in Dolichopoda cave crickets (Orthoptera, Rhaphidophoridae). BMC Evol Biol. 2009; 9: 301–314. pmid:20038292

26. Rojas AA, Vazquez-Tello A, Ferbeyre G, Venanzetti F, Bachmann L, Paquin B, et al. Hammerhead-mediated processing of satellite pDo500 family transcripts from Dolichopoda cave crickets. Nucleic Acids Res. 2000; 28:4037–4043. pmid:11024185

27. Martinsen L, Johnsen A, Venanzetti F, Bachmann L. Phylogenetic footprinting of non-coding RNA: hammerhead ribozyme sequences in a satellite DNA family of Dolichopoda cave crickets (Orthoptera, Rhaphidophoridae). BMC Evol Biol. 2010; 10:3–11. pmid:20047671

28. De Pasquale L. Misure della concentrazione di radon e variabilità genetica di popolazioni cavernicole di Dolichopoda (Orthoptera, Rhaphidophoridae). Tesi di dottorato in Scienze Ambientali, “Ambiente e Uomo in Appennino”, Università degli studi dell’Aquila. 1998.

29. Bellocchi E, Boifava F, Dal Molin L, Marchetto G. L’indice di langlier e la radioattività ambientale indoor in due grotte dell’altopiano carsico del Faedo-Casaron (Vicenza). Speleologia Veneta 2009; 17: 77–90.

30. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG. The ClustalX windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997; 24: 4876–4882. pmid:9396791

31. Librado P, Rozas J. DNA sp v.5. A software for comprehensive analysis of DNA polymorphism data. Bioinformatics 2009; 25: 1451–1452. pmid:19346325

32. Benzécri JP. L’analyse des données. II. L’analyse des correspondenses. Dunod, Paris; 1973

33. Hammer Ø, Harper DAT, Ryan PD. PAST: Paleontological statistics software package for education and data analysis. Palaeontologia Electronica. 2001; 4(1): 9pp. Available: http://palaeo-electronica.org/2001_1/past/issue1_01.htm

34. Angeletti D, Carere C. Comparative ecogenotoxicology: Monitoring the DNA of wildlife. Curr Zool. 2014; 60: 252–25

35. Dixon DR, Pruski AM, Dixon LRJ, Jha AN. Marine invertebrate eco-genotoxicology: A methodological overview. Mutagenesis 2002; 17: 495–507. pmid:12435847

36. Lewis C, Galloway T. Genotoxic damage in polychaetes: A study of species and cell-type sensitivities. Mutat Res. 2008; 654: 69–75. pmid:18579434

37. Lacaze E, Geffard O, Bony S, Devaux A. Genotoxicity assessment in the amphipod Gammarus fossarum by use of the alkaline Comet assay. Mutat Res. 2010; 700: 32–38. pmid:20451657

38. Palomeque T, Lorite P. Satellite DNA in insects: a review. Heredity 2008; 100:564–573. pmid:18414505