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.v11i2.1028210282Short CommunicationInsectaOrthopteraBiogeographyEvolutionary biologyPopulation geneticsBrazilChromosome mapping in Abracrisflavolineata (De Geer, 1773) (Orthoptera) from the Iguaçu National Park – Foz do Iguaçu, Paraná, BrazilPineMariana Bozina1GalloRaquel Bozini1da SilvaCarlos Roberto Maximiano1PezentiLarissa Forim1De DomenicoFernando Campos2LoretoVilma3da RosaRenatarenata-darosa@uel.brhttps://orcid.org/0000-0003-4258-72441Departamento de Biologia Geral, CCB, Universidade Estadual de Londrina (UEL), Londrina, Paraná, BrazilUniversidade Estadual de LondrinaLondrinaBrazilMuseu de Zoologia, Instituto de Biociências, Universidade de São Paulo (USP), São Paulo, São Paulo, BrazilUniversidade de São PauloSão PauloBrazilDepartamento de Genética, CCB, Universidade Federal de Pernambuco, Recife, BrazilCCB, Universidade Federal de PernambucoRecifeBrazil
Corresponding author: Renata da Rosa (renata-darosa@uel.br)
Academic editor: A. Bugrov
201731032017112203212FF93A671-FFA4-FFE9-D346-AE319170FFB0C0BB769F-53A6-4B31-8226-CC4DBF3890B55756122408201601022017Mariana Bozina Pine, Raquel Bozini Gallo, Carlos Roberto Maximiano da Silva, Larissa Forim Pezenti, Fernando Campos De Domenico, Vilma Loreto, Renata da RosaThis 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/C0BB769F-53A6-4B31-8226-CC4DBF3890B5
In this paper, we present the cytomolecular analysis of a population of Abracrisflavolineata collected in the largest fragment of the Brazilian Atlantic forest, the Iguaçu National Park. The diploid number in males was 23 (22+X0), with two large pairs (1–2), 7 medium (3–9), 2 small (10–11) and the X chromosome of medium size. Heterochromatic blocks were evident in the pericentromeric regions of all chromosomes. Heterogeneity in the distribution of heterochromatin was observed, with a predominance of DAPI+ blocks. However, some chromosomes showed CMA3+ blocks and other DAPI+/CMA3+ blocks. The 18S rDNA sites were distributed on the short arms of 5 pairs. In two of these pairs, such sites were in the same chromosome bearing 5S rDNA, and one of the bivalents, they were co-located. Histone H3 genes were found on one bivalent. The results added to the existing cytogenetic studies provided evidence of great karyotypic plasticity in the species. This pliancy may be the result of vicariant events related to the geographical distribution of different populations of A.flavolineata.
AcrididaeBrazilian Atlantic forestchromosome bandingfluorescence in situ hybridizationgrasshopperCitation
Bozina Pine MB, Gallo RB, da Silva CRM, Pezenti LF, De Domenico FC, Loreto V, Rosa R (2017) Chromosome mapping in Abracris flavolineata (De Geer, 1773) (Orthoptera) from the Iguaçu National Park – Foz do Iguaçu, Paraná, Brazil. Comparative Cytogenetics 11(2): 203–212. https://doi.org/10.3897/CompCytogen.v11i2.10282
Introduction
The family Acrididae is one of the most speciose, heterogeneous and conceivably the most recent in the Acridoidea group (Song et al. 2015). Among its subfamilies, Ommatolampidinae comprises 9 tribes and over 50 genera, with geographical distribution in North, Central and South America (Amédégnato 1974, Carbonell 1977). This subfamily is associated with secondary growth of vegetation in dry forest areas (Amédégnato and Decamps 1980, Rowell 1987), or even in clearings within Tropical rainforests (Braker 1991, Sperber 1996). The genus Abracris Walker, 1870 belongs to this subfamily and is distributed in the Neotropical region, from Mexico to Argentina (Roberts and Carbonell 1981).
Although there are several cytogenetic studies of Neotropical Acrididae species, especially of the subfamily Ommatolampidinae, there are few studies that elucidate the molecular structure of chromosomes. Specific C-banding techniques for the identification and classification of chromosomes have been applied more often in Abracrisflavolineata species (De Geer, 1773) (Cella and Ferreira 1991, Rocha et al. 2011).
According to Cella and Ferreira (1991), Abracrisflavolineata shows 2n=24/23 (females/males), with the XX/X0 sex chromosome system. Seven subtelocentric, two metacentric and two submetacentric pairs, and the subtelocentric chromosome X make up the karyotype of A.flavolineata. In addition, this species displayed B chromosomes (Cella and Ferreira 1991, Bueno et al. 2013). Molecular markers were applied to understand their molecular composition and mechanisms of evolution (Bueno et al. 2013, Milani and Cabral-de-Mello 2014, Palacios-Gimenez et al. 2014).
Knowledge of Ommatolampidinae is limited, considering the diversity of the subfamily. Few cytogenetic tools have been applied, and the only samples that have been analyzed were from southeastern Brazil (Rio Claro/São Paulo state). To improve the knowledge about this species, in particular about this group, it was necessary, not only to use more cytogenetic tools but also study samples from other regions of Brazil. Thus, this work aimed to study specimens of A.flavolineata collected in the Iguaçu National Park (Southern Brazil), a paramount area of the Brazilian Atlantic, using different cytogenetic markers in order to understand the evolutionary mechanisms present in this group of insects.
Material and methods
The current study used twenty male specimens of Abracrisflavolineata collected from the Iguaçu National Park, Foz do Iguaçu, Paraná State, Brazil – 25°37'40.67"S; 54°27'45.29"W (DDM). The individuals were identified and deposited in the Museu de Zoologia da Universidade de São Paulo (MZUSP) and collected with the permission of Instituto Chico Mendes de Conservação da Biodiversidade – ICMBio, protocol number 31946-2. The insects were anesthetized and dissected before fixing their testes in methanol: acetic acid (3:1). Chromosome preparations for mitotic and meiotic analyses were made through cell suspension by maceration in one drop of 45% acetic acid. Heterochromatin distribution was analyzed by Giemsa C-banding (Sumner 1972). The GC- and AT-rich bands were detected with chromomycin A3 (CMA3) and 4’-6-diamino-2-phenylindole (DAPI), respectively (Schweizer et al. 1983).
In addition to the karyotype studies, genomic DNA was extracted from the muscle tissue of a male specimen using the phenol/chloroform procedure described by Sambrook and Russel (2001). Unlabelled 18S rDNA, 5S rDNA and Histone H3 gene probes were generated by polymerase chain reaction (PCR) using the primers: 18S rDNAF 5‘-CCTG AGAAACGGCTACCACATC-3’ and 18S rDNAR 5‘-GAGTCTCGTTCGTTATCGGA-3’(Whiting 2002); 5SrDNAF 5’-AACGACCATACCACGCTGAA-3’ and 5SrDNAR 5’-AA GCGGTCCCCCATCTAAGT-3’ (Loreto et al. 2008a); H3F 5’-ATATCCTTRGGCATRAT RGTGAC-3’ and H3R 5’-ATGGCTCGTACCAAGCAGACVGC-3’ (Colgan et al. 1998). The probes isolated by PCR were labeled with digoxigenin-11-dUTP and biotin-11-dATP by PCR. The in situ hybridization procedure was performed according to Pinkel et al. (1986).
Results and discussion
All samples of Abracrisflavolineata collected in the Iguaçu National Park presented 2n=23 (♂) and an XX/X0 sex-determination system (Fig. 1). The chromosomes were classified into two large pairs (1–2), 7 medium pairs (3–9) and 2 small pairs (10–11). The X chromosome is medium-sized. Pairs 1 to 7 and the X chromosomes are subtelo-acrocentric; pairs 8 and 9 are submetacentric, and pairs 10 and 11 are metacentric (Fig. 1a). The X chromosome presented itself univalent at diakinesis (Fig. 1b).
F3CD3727-9D31-52F4-9B7E-25ADEF98C1BD
Mitotic and meiotic stages of Abracrisflavolineata by Giemsa conventional staining (a–c), C-banding (d–f) and fluorochromes staining (g–i): a mitotic metaphase b diakinesis c metaphase II d pachytene e diakinesis f metaphase II g mitotic metaphase by DAPI staining h mitotic metaphase by CMA3 staining i overlapping DAPI/CMA3. The numbers correspond to autosomes. X corresponds to the sex chromosome. The arrows indicate the discrete heterochromatic bands. Bar= 5µm.
https://binary.pensoft.net/fig/128846
The karyotype structure observed in this study is conserved regarding the diploid number and morphology of the A chromosomes, corroborating the data already established for the species by Cella and Ferreira (1991) and by Bueno et al. (2013) studying specimens from Rio Claro (southeastern Brazil). However, these authors reported the occurrence of supernumerary chromosomes, which were not found in samples from National Iguaçu Park. The sex-determination pattern considered the ancestor of grasshoppers was present in specimens of Abracrisflavolineata, with males with the X0 and females with the XX sex chromosome system. In Acrididae, some species, such as Rhammatocerusbrasiliensis (Bruner, 1904), R.brunneri (Giglio-Tos, 1895), R.palustris Carbonell, 1988, R.pictus (Bruner, 1900), Orthoscapheusrufipes (Thunberg, 1824), among others, show a karyotype similar to that of Abracrisflavolineata, with 2n = 23,X0 and 24,XX, being considered the ancestral pattern of the family (White 1973; Loreto et al. 2008b; Rocha et al. 2011).
Heterochromatic blocks were shown in the pericentromeric regions of all chromosomes (Fig. 1d–i). Chromosome 4 showed a slight heterochromatic block in the proximal region of the long arm (Fig. 1f). Nevertheless, this block was also viewed in most meiotic stages. It was possible note discrete heterochromatic bands in some chromosomes (Fig. 1d–f). These bands appear in some cells and not others, and this may be the result of differences in the condensation of chromosomes in different cells or a technical artifact. The position of heterochromatin at pericentromeric regions is a common feature in Acrididae, being observed in several species (Loreto and Souza 2000, Rocha et al. 2004). The heterochromatin plays a very important role in the evolution of karyotypes. The pericentromeric preferential distribution may be related to equilocal transfer between non-homologous chromosomes of similar size by positioning in the nucleus by bouquet configuration (Schweizer and Loidl 1987).
Staining with fluorochromes CMA3/DAPI showed heterogeneity in the distribution of heterochromatin on autosomes: (i) DAPI+ bands located at pericentromeric regions in most chromosomes; (ii) pericentromeric DAPI+/CMA3+ bands on 2 and 6 chromosomes; (iii) a discrete CMA3+ band in the terminal region on chromosome 4; and (iv) DAPI+ adjacent to CMA3+ blocks on three pairs (5, 7 and 9). The X chromosome showed adjacency of DAPI+-CMA3+-DAPI+ bands (Figure 1g-i). The distribution of these blocks can be best demonstrated in the ideogram (Figure 3).
There are three patterns of distribution of GC-rich blocks in grasshoppers: (i) CMA3+ bands related to NORs location (Souza et al. 1998, Loreto and Souza 2000, Rocha et al. 2012); (ii) CMA3+ in all chromosomes (Souza et al. 1998, Pereira and Souza 2000); and, (iii) GC-rich blocks in some chromosomes (Loreto and Souza 2000, Souza et al. 2003, Loreto et al. 2005, Souza and Melo 2007, Rocha et al. 2012). AT-rich heterochromatin is rarely encountered in Acrididae, as is the case Arcypterafusca (Pallas, 1773), A.tornosi Bolívar, 1884 (Bella and Gosálvez 1991) and Dociostaurusgenei (Ocskay, 1832) (Rodríguez Iñigo et al. 1993). In Abracrisflavolineata, studies show that most heterochromatin is GC-rich and there are no AT-rich blocks (Bueno et al. 2013). Thus, the pattern observed in specimens of A.flavolineata in this study shows a considerable karyotypic differentiation, where heterochromatinization processes, such as heterochromatin spreading, can be responsible for the differentiation of the Iguaçu National Park population, as in the grasshopper of the genus Tropidacris Scudder, 1869 (Rocha et al. 2015). Moreover, this differentiation of heterochromatin may be related to action of transposable elements (Grewal and Jia 2007).
Although the position of rDNA sites was similar to that previously reported for other populations of this species, on terminal region, the number of sites was different from that observed in previous studies. The 18S rDNA sites were distributed on terminal regions of the short arms of all analyzed specimens of 5 pairs (1, 2, 4, 7 and 10), with no variations. Two of these sites were on the same chromosome as 5S rDNA sites. These sites were co-located on one of the bivalents (Fig. 2a and b). Bueno et al. (2013) reported a variation from 5 to 9 chromosomes bearing 18S rDNA in the species of the southeastern population. Conversely, the 5S rDNA sites were distributed in the terminal portion of the bivalent 1, in the interstitial region of bivalent 2 and pericentromeric region of 7. This pattern was more conservative where two sites (1 and 2) were similar to those observed by Bueno et al. (2013). A difference in the present study was the co-location of 18S rDNA and 5S sites on pair 7, where these sequences are located in heterochromatic regions. Cabral-de-Mello et al. (2011a) analyzed the distribution of 45S and 5S rDNA in 29 species of family Acrididae, and noted the predominance of the location of sites on separate chromosomes (80.3% of the clusters). Considering this study, the co-location of DNA sites may be associated with recent evolutionary processes in the species from the Iguaçu National Park.
D5719A82-3BEA-5F80-95F5-C59C3A287E5E
Fluorescent in situ hybridization of A.flavolineata in meiocytes: a, b 18S rDNA (green) and 5S rDNA (red) probes c, d histone H3 gene probe: a pachytene b–d two cells in metaphase II. The numbers correspond to autosomes. X corresponds to the sex chromosome. The arrows indicate the histone H3 gene sites. Bar= 5µm.
https://binary.pensoft.net/fig/128847
Histone H3 genes were located on a corresponding bivalent of pair 8 (Fig. 2c and d). This location is similar in number to that found in most species of grasshoppers, distributed on a single pair of chromosomes (Cabrero et al. 2009, Cabral-de-Mello 2011a, 2011b, Neto et al. 2013). However, our data differ from those reported for the population of A.flavolineata of southeastern Brazil, in which these genes were observed on all chromosomes, including the X chromosome, but not on B chromosomes (Bueno et al. 2013). As in our results, the distribution of the histone gene was adjacent to the CMA3+ and DAPI+ heterochromatic blocks (Fig. 3), these genes are probably intercalated to heterochromatin blocks. Thus, the repetitive DNA present in the heterochromatin may play an important role in the dispersion of this sequence leading to differences in distribution between different histone genes of A.flavolineata populations studied herein.
27E40BC7-BF12-57A8-8A36-6192CAD4ECB9
Idiogram showing the mapping of different sequences studied in Abracrisflavolineata from the Iguaçu National Park.
https://binary.pensoft.net/fig/128848
Although the chromosome number and karyotype formula of the specimens studied here were the same as those described in the southeastern population (Cella and Ferreira 1991, Bueno et al. 2013), a considerable variation in the distribution of AT and GC-rich bands and specific gene sequences (18S, 5S and H3) was found. This variation refers to significant chromosomal changes, as spreading of heterochromatin can be involved in the evolution of this species. The populations studied by Cella and Ferreira (1991) and Bueno et al. (2013) inhabit Atlantic rainforest fragments, which once formed a continuous biome. The forest fragments mentioned above are approximately 790 km away from the population studied here (Iguaçu National Park). Such distance, coupled with the fact that these species probably underwent vicariance events that led to the disruption of these populations, caused reproductive isolation that led to chromosomal differences. The data obtained, together with the existing cytogenetic studies, allow us to suggest that Abracrisflavolineata has an extensive karyotype plasticity.
Acknowledgements
The authors are grateful to the National Iguaçu Park staff members for their technical assistance; to Edson Mendes Francisco for his help with the sample collection. This work was supported by Universidade Estadual de Londrina. The researchers received permission from Instituto Chico Mendes de Conservação da Biodiversidade – ICMBio to collect grasshoper specimens.
ReferencesAmédégnatoC (1974) Les genes d’acridiens néotropicaux, leur classification par families, sousfamilies et tribus.3: 193–204.AmédégnatoCDecampsM (1980) Évolution dês populations d’orthoptéres d’amazonie du nord-ouest dans les cultures traditionneles et hes formations secondaries d’origine anthropique.9: 1–33.BellaJLGosálvezJ (1991) C-banding with specific fluorescent DNA-ligands: a new approach to constitutive heterochromatin heterogeneity.1: 44–52. https://doi.org/10.3109/10520299109110549BrakerE (1991) Natural history of a neotropical gap-inhabiting grasshopper.23: 41–50. https://doi.org/10.2307/2388686BuenoDPalacios-GimenezOMCabral-de-MelloDC (2013) Chromosomal mapping of repetitive DNAs in the grasshopper Abracrisflavolineata reveal possible ancestry of the B chromosome and H3 histone spreading8: e66532. https://doi.org/10.1371/journal.pone.0066532Cabral-de-MelloDCCabreroJLópez-LeónMDCamachoJPM (2011a) Evolutionary dynamics of 5S rDNA location in acridid grasshoppers and its relationship with H3 histone gene and 45S rDNA location.139: 921–931. https://doi.org/10.1007/s10709-011-9596-7Cabral-de-MelloDCMartinsCSouzaMJMouraRC (2011b) Cytogenetic mapping of 5S and 18S rRNAs and H3 histone genes in 4 ancient Proscopiidae grasshopper species: contribution to understanding the evolutionary dynamics of multigene families.132: 89–93. https://doi.org/10.1159/000317476CabreroJLópez-LeónMDTeruelMCamachoJPM (2009) Chromosome mapping of H3 and H4 histone gene clusters in 35 species of acridid grasshoppers.17: 397–404. https://doi.org/10.1007/s10577-009-9030-5CarbonellCS (1977) Origin, evolution and distribuition of the neotropical acridomorph fauna (Orthoptera): A preliminary hypothesis.36: 153–175.CellaDMFerreiraA (1991) The cytogenetics of Abracrisflavolineata (Orthoptera, Caelifera, Ommatolampinae, Abracrini).14: 315–329. http://www.gmb.org.br/Revistas/V14/v14a23.pdfColganDJMcLauchlanAWilsonGDFLivingstonSPEdgecombeGDMacaranasJCassisGGrayMR (1998) Histone H3 and U2 snRNA DNA sequences and arthropod molecular evolution.46: 419–437. https://doi.org/10.1071/ZO98048GrewalSISJiaS (2007) Heterochromatin revisited.8: 35–46. https://doi.org/10.1038/nrg2008LoretoVSouzaMJ (2000) Karyotype, constitutive heterochromatin and nucleolar organizer regions (NORs) in Belosacriscoccineipes (Acrididae-Leptysminae).23: 575–579. https://doi.org/10.1590/S1415-47572000000300013LoretoVStadtlerEde MeloNFSouzaMJ (2005) A comparative cytogenetic analysis between the grasshopper species Chromacrisnuptialis and C.speciosa (Romaleidae): constitutive heterochromatin variability and rDNA sites.125: 253–260. https://doi.org/10.1007/s10709-005-0368-0LoretoVCabreroJLópez-LeónMDCamachoJPMSouzaMJ (2008a) Possible autosomal origin of macro B chromosomes in two grasshopper species.16: 233–241. https://doi.org/10.1007/s10577-007-1188-0LoretoVCabreroJLópez-LeónMDCamachoJPMSouzaMJ (2008b) Comparative analysis of rDNA location in five Neotropical gomphocerine grasshopper species.132: 95–101. https://doi.org/10.1007/s10709-007-9152-7MilaniDCabral-de-MelloDC (2014) . Microsatellite organization in the grasshopper Abracrisflavolineata (Orthoptera: Acrididae) revealed by FISH mapping: remarkable spreading in the A and B Chromosomes. 9: e97956. https://doi.org/10.1371/journal.pone.0097956NetoMSRSouzaMJLoretoV (2013) Chromosomal evolution of rDNA and H3 histone genes in representative Romaleidae grasshoppers from northeast Brazil. 6: 41. https://doi.org/10.1186/1755-8166-6-41Palacios-GimenezOMBuenoDCabral-de-MelloDC (2014) Chromosomal mapping of two Mariner-like elements in the grasshopper Abracrisflavolineata (Orthoptera: Acrididae) reveals enrichment in euchromatin.111(3): 329–334. https://doi.org/10.14411/eje.2014.052PereiraLGSouzaMJ (2000) Nature and distribution of constitutive heterochromatin and NOR location in the grasshopper Phaeopariamegacephala (Romaleidae: Orthoptera).103: 111–119.PinkelDStraumeTGrayJW (1986) Cytogenetic analysis using quantitative, high-sensitivity, fluorescence hybridization.83: 2934–2938. https://doi.org/10.1073/pnas.83.9.2934RobertsHRCarbonellCS (1981) A revision of the neotropical genus Abracris and related genera (Orthoptera, Acrididae, Ommatolampinae).133: 1–14. http://www.jstor.org/stable/4064766RochaMFSouzaMJMouraRC (2004) Karyotype analysis, constitutive heterochromatin and NOR distribution in five grasshopper species of the subfamily Leptysminae (Acrididae).57: 107–116. https://doi.org/10.1080/00087114.2004.10589378RochaMFMeloNFSouzaMJ (2011) Comparative cytogenetic analysis of two grasshopper species of the tribe Abracrini (Ommatolampinae, Acrididae).34: 214–219. https://doi.org/10.1590/S1415-47572011000200008RochaMFSouzaMJMouraRC (2012) Karyotypic analysis, constitutive heterochromatin and NOR distribution in five grasshopper species of the subfamily Leptysminae (Acrididae).57: 107–116. https://doi.org/10.1080/00087114.2004.10589378RochaMFPineMBOliveiraEFASLoretoVGalloRBDa SilvaCRMDe DomenicoFCDa RosaR (2015) Spreading of heterochromatin and karyotype differentiation in two Tropidacris Scudder, 1869 species (Orthoptera, Romaleidae).9(3): 435–450. https://doi.org/10.3897/CompCytogen.v9i3.5160Rodriguez-IñigoEBellaJLLaVega CG (1993) Heterochromatin differentiation between two species of the genus Dociostaurus (Orthoptera: Acrididae).70: 458–465. https://doi.org/10.1038/hdy.1993.67RowellCHF (1987) The biogeography of Costa Rican acridid grasshoppers, in relation to their phylogenetic origins and ecology. In: BaccettiBM (Ed.) Evolutionary biology of orthopteroid insects., 470–482.SambrookJRussellDW (2001) Cold Spring Harbor, New York, 1659 pp.SchweizerDLoidlJ (1987) A model for heterochromatin dispersion and the evolution of C-bands patterns. In: StahlALucianiJMVagner-CapodanoAM (Eds) Chromosomes Today., 61–74. https://doi.org/10.1007/978-94-010-9166-4_7SchweizerDMendelakMWhiteMJDContrerasN (1983) Cytogenetics of the parthenogenetic grasshopper Warramabavirgo and its bisexual relatives. X. Pattern of florescent banding.88: 227–236. https://doi.org/10.1007/BF00285625SongHAmédégnatoCCiglianoMMDesutter-GrandcolasLHeadsSHHuangYOtteDWhitingMF (2015) 300 million years of diversification: elucidating the patterns of orthopteran evolution based on comprehensive taxon and gene sampling. . 1–31. https://doi.org/10.1111/cla.12116SouzaMJRufasJSOrellanaJ (1998) Constitutive heterochromatin, NOR location and FISH in the grasshopper Xyleusangulatus (Romaleidae).51: 73–80. https://doi.org/10.1080/00087114.1998.10589122SouzaMJHaverPROMeloNF (2003) Karyotype, C- and florescence banding patterns,NOR location and FISH in the grasshopper Xestotrachelusrobustus (Romaleidae).56: 261–267. https://doi.org/10.1080/00087114.2003.10589333SouzaMJMeloNF (2007) Chromosome study in Schistocerca (Orthoptera-Acrididae-Cyrtacanthacridinae): Karyotypes and distribution patterns of constitutive heterochromatin and nucleolus organizer regions (NORs).30: 54–59. https://doi.org/10.1590/S1415-47572007000100011SperberCF (1996) Field diet of the grasshopper Abracrisdilecta Walker (Orthoptera, Acrididae).13: 127–135. https://doi.org/10.1590/S0101-81751996000100012SumnerAT (1972) A simple techinique for demonstrating centromeric heterochromatin.75: 304–306. https://doi.org/10.1016/0014-4827(72)90558-7WhiteMJD (1973) Animal cytology and evolution. 3rd ed. Cambridge University Press, London, 455p.WhitingMF (2002) Mecoptera is paraphyletic: multiple genes and phylogeny of Mecoptera and Siphonaptera.31: 93–104. https://doi.org/10.1046/j.0300- 3256.2001.00095.x