- Research article
- Open Access
Uncovering the evolutionary history of neo-XY sex chromosomes in the grasshopper Ronderosia bergii (Orthoptera, Melanoplinae) through satellite DNA analysis
BMC Evolutionary Biologyvolume 18, Article number: 2 (2018)
Neo-sex chromosome systems arose independently multiple times in evolution, presenting the remarkable characteristic of repetitive DNAs accumulation. Among grasshoppers, occurrence of neo-XY was repeatedly noticed in Melanoplinae. Here we analyzed the most abundant tandem repeats of R. bergii (2n = 22, neo-XY♂) using deep Illumina sequencing and graph-based clustering in order to address the neo-sex chromosomes evolution.
The analyses revealed ten families of satDNAs comprising about ~1% of the male genome, which occupied mainly C-positive regions of autosomes. Regarding the sex chromosomes, satDNAs were recorded within centromeric or interstitial regions of the neo-X chromosome and four satDNAs occurred in the neo-Y, two of them being exclusive (Rber248 and Rber299). Using a combination of probes we uncovered five well-defined cytological variants for neo-Y, originated by multiple paracentric inversions and satDNA amplification, besides fragmented neo-Y. These neo-Y variants were distinct in frequency between embryos and adult males.
The genomic data together with cytogenetic mapping enabled us to better understand the neo-sex chromosome dynamics in grasshoppers, reinforcing differentiation of neo-X and neo-Y and revealing the occurrence of multiple additional rearrangements involved in the neo-Y evolution of R. bergii. We discussed the possible causes that led to differences in frequency for the neo-Y variants between embryos and adults. Finally we hypothesize about the role of DNA satellites in R. bergii as well as putative historical events involved in the evolution of the R. bergii neo-XY.
Evolutionary evidence indicates that sex chromosomes have evolved independently from ordinary autosomal pairs (the proto-XY or the proto-WZ pairs) comprising one of the most dynamic chromosomal elements in the genome [1,2,3]. The sex chromosome systems are diverse with, for example, the XY system in mammals, the WZ system in Lepidoptera, birds and some reptiles. Moreover, the occurrence of multiple sex chromosomes (multiple Xs and Ys or Zs and Ws) is also a common phenomenon [4,5,6,7,8]. Common hallmarks of sex chromosome evolution include the partial or complete loss of recombination between them; abundant gene inactivation or loss (i.e., genetic degeneration); progressive accumulation/expansion of repetitive DNAs; and heterochromatinization of the Y or W chromosomes [6, 7, 9,10,11,12,13,14].
In Orthoptera the X0♂/XX♀ sex chromosome system is a conserved pattern found in most species. This modal system results from the loss of the Y chromosome in an ancestral XY♂/XX♀ species [15, 16]. However, derived variants, like neo-XY♂/neo-XX♀, neo-X1X2Y♂/neo-X1X1X2X2♀ or even X1X20♂/X1X1X2X2♀ evolved several times due to chromosomal rearrangements involving autosomes and sex chromosomes [15,16,17,18,19,20,21]. In animals, it is particularly evident that chromosomal rearrangements trigger recombination suppression because of the lack of homologous sequences in regions participating in the rearrangement [22,23,24,25]. The establishment of inversions, for example, is typically associated with recombination suppression and has been detected in the Y chromosome of the grasshopper Ronderosia bergii , mammals , the three-spine stickleback fish , and the plant Silene latifolia (Nicolas et al. ).
Ronderosia bergii is a Melanoplinae grasshopper with a 2n = 22, neo-XY♂ karyotype. The neo-X is metacentric and the neo-Y is acrocentric, with a larger small arm in comparison to autosomes, and rich in heterochromatin; they are the result of an X-A centric fusion followed by a large pericentric inversion involving more than 90% of the neo-Y length. The inversion restricted the segment of contact between the chromosomes to the distal regions during meiosis [20, 27, 28]. There is no published estimate about the origin time of Ronderosia genus to accurately estimate the age of the neo-XY system in R. bergii. However, the genus is endemic to South America  and members of the Melanoploid lineage have been present for at least 46 mya in this region . Classical cytogenetic and molecular analyses suggest that the R. bergii neo-Y has acquired derived features as well as retained ancestral features after its origin. For instance, the neo-Y chromosome accumulated repetitive DNA sequences while retaining the post-translational histone modifications of autosomes .
The heterochromatic nature of the R. bergii neo-Y chromosome suggests the accumulation of repetitive DNAs (e.g., satellite DNAs, satDNAs), as detected in the neo-Y of the plants Silene latifolia  and Rumex acetosa , the cervid Muntiacus muntjac [33, 34], among others. Here we performed a characterization of the most abundant satDNAs in the R. bergii genome. The analyses are mostly focused on the neo-XY sex chromosomes, and the expectation that information about sex chromosome origin and evolution could be obtained from satDNAs. We found that satDNAs predominate at centromeric heterochromatin of autosomes and neo-X chromosome, and throughout the neo-Y. Interestingly, FISH chromosomal localization of satDNAs in both embryonic neuroblasts and adult male metaphases revealed five neo-Y variants, explained by the occurrence of paracentric inversions involving two of the satDNAs identified here. In addition, amplification of two other satDNAs significantly increased the size of the neo-Y. This data brings new information about neo-XY sex chromosome evolution in R. bergii, revealing that the evolutionary history of these chromosomes is much more complex than previously thought.
Sampling, chromosome obtaining and DNA extraction
Adult animals (males and females) of Ronderosia bergii were collected in the Parque Estadual Edmundo Navarro de Andrade (Rio Claro, SP, Brazil) between May/2013 and March/2015 under authorization of COTEC (process number 341/2013). Testes were removed and fixed in Carnoy’s solution (3:1, 100% ethanol: absolute Acetic Acid). For embryo obtaining females were kept in captivity until oviposition. Mitotic chromosomes were obtained following the protocol proposed by Webb et al. . The phenol/chloroform-based procedure described in Sambrook and Russel  was used to obtain genomic DNA from males and females. For karyotypic analyses it was used the conventional staining with 5% Giemsa. White’s  terminology was used to recognize the arms of the neo-X chromosome at meiosis: the arm from the original X chromosome was referred to as XL and the arm that shares homology with the neo-Y was referred to as XR.
Illumina sequencing and graph-based clustering
Paired-ends libraries (2 × 300) were prepared as recommended by Illumina (Illumina Inc., San Diego, CA, USA), using Nextera DNA kit and genomic DNA isolated from a single male. Paired-end reads were preprocessed to check the quality of sequenced nucleotides with FASTQC  and quality trimmed using the FASTX-Toolkit suit . The trimmed paired-end reads were joined using the “fastq-join” software of the FASTX-Toolkit suit  with default options. To search for satDNAs in the R. bergii genome, 4,702,802 joined reads were used in the RepeatExplorer pipeline [39, 40]. We then searched for clusters showing higher graph density, which is a typical feature of satDNAs families analysis [39, 40] due to tandem repeats characteristic of the sequence.
satDNAs isolation and sequence analysis
The contigs from clusters displaying high graph density were submitted to Tandem Repeats Finder (TRF) algorithm ; TRF enabled the identification of DNA sequences that maximized the alignment scores between the different tandemly repeated monomers. We used 2, 3, 5 parameters for match, mismatch and indels, respectively in TRF analysis, and the score of 50 for minimum alignment was considered for reporting. To confirm the tandem organization and to identify monomers of the same family, we used the dotplot graphic alignment tool implemented in Dotlet . To identify a representative copy of a given satDNA family, the monomers with maximum length was taken. To check similarity with published sequences, monomers were used as the query in BLAST and Repbase searches using the National Center for Biotechnology Information (NCBI)  and RepeatMasking tools . In addition, these canonical monomers were BLASTed against the satellitome of Locusta migratoria  and Eumigus monticola . satDNAs copies were aligned using Muscle  implemented in MEGA5 . Nucleotide divergence (p divergence) and A + T content were estimated in MEGA5.
The presence of satDNAs families in both male and female were verified by polymerase chain reactions (PCR). The primers were designed based on consensus sequence of each satDNA family using the Primer3 software  or manually. PCR mix reactions contained 10× PCR Rxn Buffer, 0.2 mM MgCl2, 0.16 mM dNTPs, 2 mM of each primer, 1 U of Taq Platinum DNA Polymerase (Invitrogen, San Diego, CA, USA) and 50-100 ng/μl of template DNA. The PCR conditions were as follows: initial denaturation at 94 °C for 5 min and 30 cycles at 94 °C (30 s), 55 °C (30 s), and 72 °C (80 s), plus a final extension at 72 °C for 5 min. Monomeric bands were isolated and purified from agarose gel 1% using the Zymoclean™ Gel DNA Recovery Kit (Zymo Research Corp., The Epigenetics Company, USA) according to the recommendations of the fabricant. Afterwards, the monomers were used as template for PCR reamplification.
The purified PCR products were sequenced in both directions (Macrogen Inc., South Korea) in order to determine the isolation of the sequences of interest. The sequenced products were compared to the consensus sequences obtained by genome analysis. The monomer consensus sequences belonging to each of the ten satDNAs families obtained from genomic analysis were deposited into the NCBI database under the following accession numbers: MF765804-MF765813.
Probes and Fluorescence In situ Hybridization, measurement of sex chromosomes and distance between signals
The PCR products for each satDNA family with more than 50 bp were labeled by nick translation using biotin-14-dATP (Invitrogen) or digoxigenin-11-dUTP (Roche, Mannheim, Germany). SatDNAs with less than 50 bp were labeled directly at the 5′ end with biotin-14 dATP (Sigma-Aldrich, St Louis, MO, USA) during their synthesis. For single or two-color FISH the protocols proposed Pinkel et al.  with modifications  were followed using mitotic chromosome preparations. Probes labeled with digoxigenin-11-dUTP were detected using anti-digoxigenin rhodamine (Roche), while the probes labeled with biotin-14-dATP were detected using Streptavidin Alexa Fluor 488-conjugated (Invitrogen). The preparations were mounted in VECTASHIELD (Vector, Burlingame, CA, USA) with chromosomes counterstained using 4′,6-diamidine-2′-phenylindole (DAPI).
The sex chromosome sizes and the distance between some FISH signals present in neo-Y chromosome were measured using the software ImageTool version 3.0 (developed at the University of Texas Health Science Center at San Antonio, Texas and available from the Internet by anonymous FTP from maxrad6.uthscsa.edu).
SatDNAs identification and sequence analysis
Clustering analysis with the RepeatExplorer pipeline [39, 40] produced 1,030,168 clusters (containing 91,4% of the raw reads) with differences in size and genome composition. It included satDNAs and other non-characterized repetitive elements. A low copy number fraction of the genome is represented by singletons containing 8.6% of reads. For satDNAs searching the most abundant clusters (551 clusters) representing repetitive elements were examined.
The analysis of clusters generated by RepeatExplorer and using TRF algorithm, besides dotplots confirmed that 10 of the 551 clusters were resolved to be satDNAs sequences, with repeat units ranging from 5 to 285 bp long. Among them, six satDNAs sequences presented >50 bp length and were experimentally validated through PCR. The analyses revealed a ladder pattern in the gel electrophoresis that is characteristic for this kind of repetitive sequence. For the remaining four repeat families, with monomers <50 bp, we were not able to confirm them due to the very limited sequence window for PCR primer design and amplification. Features for satDNAs families are detailed in the Table 1. According to RepeatExplorer output all 10 satDNAs families together represented about 1.071% of the male genome, with the most abundant representing 0.551% and the less abundant 0.011% of the genome. Within satDNAs families the nucleotide divergence ranged from 6.2 to 17.9%. A + T content varied from 39.1 to 64.8% (Table 1). NCBI BLAST, internal BLAST, and Repbase searches performed using the consensus monomer sequence belonging to each repeat family as query did not reveal significant similarity with described sequences. The monomer consensus sequences of each satDNA family are shown in the Additional file 1. The satellites alignments are available upon request to the author.
Chromosomal localization of satDNAs
The ten satDNAs we identified were represented on both autosomal and sex chromosomes. As expected, satDNAs were localized primarily in segments of heterochromatin (i.e. pericentromeric areas), although non-heterochromatic segments were also labeled, including interstitial, proximal or terminal regions. Most of the 10 satDNAs occurred in several chromosomes but none was present in all chromosomes (Fig. 1, Table 1). Some of the satDNAs were localized on some autosomes (Rber1, Rber59, Rber61, Rber370, Rber491, Rber520), while other were restricted to few chromosome pairs, like Rber158, Rber185 (Fig. 1, Table 1). Two satDNAs were exclusively visualized through FISH in the neo-Y chromosome [Rber248 (Fig. 1c) and Rber299 (Fig. 1b)], with proximal segments, but not co-localized (Fig. 2). These two satDNAs were involved in paracentric inversions (see below). Besides these two satDNAs the neo-Y also harbored Rber1 and Rber158. Rber1 was mapped to a small proximal band while Rber158 was localized in three different segments in the long arm and one segment in the short arm of the neo-Y chromosome (Fig. 1a). Finally, seven satDNAs were present in the neo-X chromosome, six of them in the centromeric region and one (Rber158) in the interstitial position (Fig. 1, Table 1).
Variable neo-Y chromosomes with different satDNAs distribution
Intriguingly, we observed inter-individual variation in neo-Y chromosomes size and variable localization of the satDNAs that are exclusively placed in the neo-Y (Rber248 and Rber299). Regarding the size of the neo-Y two variants were observed in distinct males, one in which the neo-Y is almost half size of the neo-X chromosome and another one in which the neo-Y is about 0.78 of the neo-X chromosome observed in one individual (Fig. 3a; Additional file 2).
Concerning the location of Rber248 and Rber299, we studied 15 adult males and 22 embryos. The study of 22 embryos sampled randomly from distinct ooteca revealed four inter-individual distinct variants, named as I, II, III and IV. Pattern I was the most common; it occurred in 68% of the embryos (15 embryos) and showed that Rber248 was proximal to centromere while Rber299 was located just below this marker (Fig. 2a). In pattern II, Rber248 was proximal to the centromere while the Rber299 was located interstitially around the middle of the neo-Y chromosome; the distance between the Rber248 and Rber299 was similar to the distance between the Rber299 and the end of the long arm. This situation occurred in 22.7% of the embryos (5 embryo) (Fig. 2b; Additional file 3). In pattern III Rber248 was located in the same position while Rber299 was interstitial, like in the pattern II, but closer to the Rber248, occurring in only one embryo (4.54% of the sample). The distance between the two satDNAs was closer than between Rber299 and the end of the long arm (Fig. 2c; Additional file 3). For the three variants I, II and III the neo-Y had a similar size as shown in the Fig. 3a. Variant IV occurred in only one individual (4.54%) and in this variant the neo-Y was larger than in I, II and III. Moreover the signals for the satDNAs Rber248 and Rber299 were positioned in the middle of the chromosome, but in inverse order in relation to the centromere than the other variants, i.e., the Rber299 was more proximal than the Rber248 (Fig. 2d).
For variant IV, we also mapped three other satDNAs repeats that were observed in the neo-Y variant I (the most common in mitosis), i.e., Rber1, Rber59 and Rber158. These three satDNAs were found arranged in distinct position on the neo-Y, and also differing in clusters number in comparison to variant I. The Rber1 and Rber59 were located in centromere, extending for the short arm, besides a proximal signal in the long arm. For Rber1 additional signal with interstitial position was also observed, similar to those observed in variant I (Fig. 3b). The Rber158, as observed in the variant I presented three discrete signals in interstitial position, but the band that occurred in the short arm of variant I was not noticed (Fig. 3b).
In order to verify the neo-Y polymorphisms in the germline of adult males, we also performed analysis of meiotic cells at metaphase I. Among the four variants seen in embryos only three were observed in adults, i.e., I, II and III (Fig. 2a-c), and additionally another distinct variant (named variant V) not noticed in embryos was recognized, in which the Rber299 was placed terminally (Fig. 2e). Frequencies for the variants were different in comparison to those observed in embryos, i.e., variant I 7.15% (1 individual), variant II 57.1% (8 individuals), variant III 28.6% (4 individuals) and variant V 7.15% (1 individual). Additionally, in meiosis we found a neo-Y chromosome that was not possible to include in any of the described variants. This neo-Y was folding on itself at metaphase I, fragmented but joined by a thin chromatin in early anaphase I. The meiotic signals for Rber248 and Rber299 remained in the same position as those observed for mitotic variant I, II, II and V (Fig. 4).
General organization of satDNAs in R. bergii genome
In grasshoppers RepeatExplorer satDNAs analysis was successfully applied in three species, revealing clues about organization of satDNAs in chromosomes of Schistocerca gregaria , general patterns of evolution of satDNAs using in Locusta migratoria  and the origin of the B chromosome in Eumigus monticola . In comparison to the other species, the R. bergii satDNAs revealed differences concerning number of satDNA loci, clusterization for distinct families and specific chromosomal distribution, e.g., clustered in R. bergii and S. gregaria; clustered, dispersed and intermingled in E. monticola and L. migratoria, besides absence of band-like pattern detected through FISH [45, 46, 52]. In R. bergii, the absence of one satDNA occupying all centromeric or proximal areas of centromeres is remarkable, highlighting the differential centromere composition between chromosomes. In S. gregaria , L. migratoria  and E. monticola  at least one satDNA repeat was placed near to centromere and for E. monticola the pericentromeric location of EmoSat08-41 suggested involvement in centromeric function . Our data suggest that if the satDNAs of R. bergii play a centromeric function as in other organisms [53,54,55], distinct centromeres are governed by distinct satDNA repeats or a combination of them. On the other hand, it is possible that these satellites might not play a centromeric function.
A noticed similarity between the three grasshopper species with satDNAs studied until now is the predominance of satDNAs families rich in A + T base pairs, suggesting that heavy satDNAs could be a common feature of grasshopper genomes. Considering that species studied belongs to distantly phylogenetically related families or subfamilies the sharing of A + T-rich satDNAs could be an ancient, which was conserved after satDNA library divergence. This picture is completely different from that observed in the cricket Eneoptera surinamensis in which the analysis of satDNAs revealed predominance of G + C-rich satDNAs , suggesting a more complex history for satDNA library in the ancestor of Ensifera and Celifera orthopterans, which diverged about 300 mya . The most common A + T content of the R. bergii satDNAs is also a common hallmark for satDNAs embedded in the centromeric heterochromatin [53, 58]. A + T-rich centromeric satDNAs have been reported in other insect, such as Tribolium castaenum , T. brevicornis  and Drosophila species .
Our analysis revealed the occurrence of ten satDNAs that represent only about 1.072% of R. bergii genome. Even though grasshoppers present large genomes, the satDNAs represent a small fraction of their repetitive DNA content, e.g., L. migrtoria 2.39 and 2.74% in Southern and Northern lineages, respectively  and E. monticola about 1.79 and 1.91% in 0B and +B genomes, respectively . An exception is S. gregaria where the characterized satDNAs correspond to about 30.4% of the genome content . High amount of satDNAs were also reported in other insects, like the TCAST repeat that makes up 35% of the centromeric heterochromatin in Tribolium castaneum genome [61, 62]. These data suggest that in grasshoppers other repetitive elements could be more relevant for genome size increasing, like the Transposable Elements (TEs), as documented in L. migratoria with 60%, being ~24% and ~17% represented by DNA transposons and LINE retroelements, respectively .
Ronderosia bergii satDNAs lacks completely homology of satDNAs libraries from other grasshoppers, suggesting that satDNAs evolved repeatedly and independently diverging in each genome. Based on the divergence between satDNAs of L. migratoria and E. monticola Ruiz-Ruano et al.  proposed that in grasshoppers the satDNA lifespan is shorter than 100 mya, which according to Song et al.  corresponds to the most recent common ancestor for the two families to which these species belong, i.e. Acrididae and Pampharidae, respectively. Our data support this hypothesis considering that the diverging time between the Acrididae subfamilies Melanoplinae (R. bergii) and Oedipodinae (L. migratoria) is shorter than 100 mya .
SatDNAs content highlights cryptic neo-Y variants and illuminates the complex evolutionary history of the R.bergii neo-sex chromosomes
The mapping of satDNAs and their distribution in the neo-XY bring new information about sex-chromosome composition and evolution in R. bergii. In a previous paper Palacios-Gimenez et al.  described the occurrence of distinct repetitive DNAs in the R. bergii sex chromosomes, like microsatellite and multigene families, highlighting the accumulation of repetitive DNAs and for some of them distinct accumulation between the neo-X and neo-Y, also distinct heterochromatinization. Here it is documented that the satDNAs are also composing the neo-XY, and could influence their evolution and differentiation through amplification for the ten distinct families mapped. Our findings confirm our previous hypothesis that the neo-sex chromosomes are highly divergent in R. bergii, where homologous sequences are completely loss between them, at least with the mapped repetitive DNAs. This differential accumulation could be documented, for example, by Rber59, Rber61, Rber370, Rber491, Rber520 sequences that are accumulated in the centromere of neo-X, but they are absent in the neo-Y; as well as by accumulation of Rber158 in neo-Y and the exclusive observation of Rber248 and Rber299 in the neo-Y. Differential patterns concerning the accumulation of repetitive DNAs and heterochromatinization are common hallmark of many sex chromosome systems with suppressed recombination [9, 64,65,66,67,68], and were reported in several well studied insect species, such as in the Y of Drosophila melanogaster  and D. miranda  and in the W chromosomes of many lepidopteran species (e.g. ).
The mapping of the two neo-Y specific satDNAs, i.e., Rber248 and Rber299, revealed the occurrence of five cryptic variants for this chromosome, caused most parsimonially by paracentric inversions after neo-XY establishment and after pericentric inversion established in the neo-Y. In animals, it is evident that chromosomal inversions trigger recombination suppression because of the lack of homologous sequences in those regions involved in the rearrangement [22,23,24,25]. The establishment of inversions is typically associated with recombination suppression and has been detected in other species with Y chromosome, such as in mammals , the three-spine stickleback fish , and the plant Silene latifolia , among others.
Taken together, the former considerations and presented data lead us to suggest a model of evolution for the R. bergii neo-XY based on the satDNAs distribution patterns (Fig. 5). First the de novo origin of the neo-sex chromosomes involved an X-A centric fusion [27, 28] that is the most accepted hypothesis for the neo-sex chromosome origin in Orthoptera [15,16,17]. After the neo-XY arisen, recombining events are still able to occur between the XR arm, which corresponds to the ancestral autosome involved in the rearrangement, and the neo-Y, as documented in other Melanoplinae species by chiasmata analysis [17, 71]. After the neo-sex chromosomes have been established the C-shape disposition of neo-XY in meiotic metaphase led some authors to suggest the occurrence of a large pericentric inversion [15, 16] involving more than 90% of the neo-Y followed by its heterochromatinization and recombination suppression between the neo-X and the neo-Y [15, 16, 20]. These rearrangements and recombination suppression could have favored the accumulation of the satDNAs in the neo-Y, like Rber248 and Rber299 clustered near each other in the proximal region of long arm, nearby the centromere, as well as Rber1 in proximal position and Rber158 forming multiple clusters along the neo-Y. These characteristics were observed in the variant named I and specifically the Rber248 was closer to the centromere than Rber299. Afterward, it is more parsimonious from the variant I the occurrence of sets of chromosomal rearrangements generating the other variants, as follows: (i) three large and independent paracentric inversions involving Rber299 displacing this satDNAs towards interstitial or distal regions originating the variants II, III and V; (ii) a paracentric inversion involving both Rber248 and Rber299 and (iii) amplifications of Rber1 and Rber59 (not observed in the variant I, probably by low copy number and FISH resolution) generating the pattern IV, in which the neo-Y is enlarged. Our hypothesis as well as the occurrence of paracentric inversions are supported by the satDNA distribution observed in both the mitotic and meiotic cells (see bellow).
Some remarkable differences regarding the occurrence and frequency of variants were noticed between embryos and meiotic cells from adult males. For example, the variant V was observed exclusively in adults, while the variant IV was seen in embryos. In addition, a variant, that not fall in the other five ones was recognizable in meiosis. This neo-Y is folding itself at metaphase I and fragmented but joined by a thin chromatin at anaphase I. We suggested that this variant could be the result of multiple breakages involved in the inversions. Even though the neo-Y is largely non-recombining with the neo-X, this meiotic configuration could facilitate the propagation of tandem repeats by intra-chromosomal crossing over, which is of the mechanisms proposed for the propagation of tandem repeats along the genome [9, 72, 73], but it deserves more investigation.
The differences between frequency for the five variants observed in R. bergii embryos and adults, mainly for variants I and II, could be directly implicated with the evolutionary patterns at populational level involving selection and genetic drift, like proposed for other chromosomal rearrangments . We noticed that the variant I occur in 68.2 and 7.15% in mitosis and meiosis, respectively, while the variant II result in 22.7 and 57.1% in mitosis and meiosis, respectively. It is more parsimonious to explain these differences as results of genetic drift. The higher frequency of variant II in adults could be favored probably by selection; however, we are aware about this second explanation because we do not have enough evidences of selection acting on the neo-Y chromosome. In any case a larger sample and animals from other populations should be studied to confirm these statements. If inversions in R. bergii are evolving neutrally and variation in their frequency is only driven by genetic drift, it must be expected that similar proportion for the two variants would be maintained, both in embryos and adults. But our data suggests that variant II is more favored, being more frequent in adults. Other populations should be investigated to test if this pattern is a local adaptation or different geographical niches support different neo-Y variants. The other three variants should be less favored than the variant II, occurring in lower frequencies in adults.
Taken together, our results show accumulation and repositioning of satDNAs arrays resulting in an interesting cryptic polymorphisms for the R. bergii neo-Y. This suggests a complex evolutionary history of the neo-sex chromosomes than previously thought (see ), shedding light on understanding patterns of sex chromosome evolution in Orthoptera. There is no other case described until now in Orthoptera in which an established neo-Y chromosome is highly variable as in R. bergii. Because neo-sex chromosomes like those present in R. bergii could be an agent of species-level selection further analysis must be performed. Another interesting issue that allows to address these findings is at populational level. A large numbers of traits have been associated with inversion polymorphism in several taxa . Ronderosia bergii is characterized by the highly variable patterns of body colour along its distribution range and in some cases within the same locality . Until now, we do not know if in other distant populations of R. bergii, another neo-XY variant is favored by natural selection, and the frequencies found in meiosis and mitosis change following a polymorphism. If this occurs, it would be an interesting fact that could contribute to the knowledge of complex polymorphic populations presenting both the chromatic and morphological variants. Also, the neo-sex chromosomes may be involved in the evolution of complex population variants of R. bergii, which, in some cases, were erroneously described as different species.
FISH mapping using these markers in different populations, which included this phenotypic variations, will help us to understand the actual state of polymorphisms, and its possible influence on population isolation. In addition, controlled crosses using male with known neo-Y variant and analysis of offspring will give the possibility to test if inversions are or not occurring de novo.
Fluorescence in situ hybridization
Polymerase chain reaction
Ohno S. Sex chromosomes and sex linked genes. Berlin: Springer; 1967.
Bull JJ. Evolution of sex determining mechanisms. Menlo Park: Benjamin Cummings; 1983.
Rice WR. Evolution of the Y sex chromosome in animals. Bioscience. 1996;46:331–43.
McMillan D, Miethke P, Alsop AE, et al. Characterizing the chromosomes of the platypus (Ornithorhynchus anatinus). Chromosom Res. 2007;15:961–74.
Traut W, Sahara K, Marec F. Sex chromosomes and sex determination in Lepidoptera. Sex Dev. 2007;1:332–46.
Bergero R, Charlesworth D. The evolution of restricted recombination in sex chromosomes. Trends Ecol Evol. 2008;24:94–102.
Ming R, Bendahmane A, Renner SS. Sex chromosomes in land plants. Annu Rev Plant Biol. 2011;62:485–514.
Pokorná M, Giovannotti M, Kratochvíl L, et al. Strong conservation of the Z bird chromosome in reptilian genomes is revealed by comparative painting despite 275 million years divergence. Chromosoma. 2011;120:455–68.
Charlesworth B, Sniegowski P, Stephan W. The evolutionary dynamics of repetitive DNA in eukaryotes. Nature. 1994;371:215–20.
Skaletsky H, Kuroda-Kawaguchi T, Minx PJ, et al. The male- specific region of the human Y chromosome is a mosaic of discrete sequence classes. Nature. 2003;423:825–37.
Nicolas M, Marais G, Hykelova V, et al. A gradual process of recombination restriction in the evolutionary history of the sex chromosomes in dioecious plants. PLoS Biol. 2005;3:e4.
Matsubara K, Tarui H, Toriba M, et al. Evidence for different origin of sex chromosomes in snakes, birds, and mammals and step-wise differentiation of snake sex chromosomes. Proc Natl Acad Sci U S A. 2006;103:18190–5.
Bergero R, Qiu S, Charlesworth D. Gene loss from plant sex chromosome system. Curr Biol. 2015;25:1234–40.
Wright AE, Dean R, Zimmer F, Mank JE. How to make a sex chromosome. Nat Commun. 2016;7:12087.
White MJD. Animal cytology and evolution. Cambridge: Cambridge University Press; 1973.
Hewitt GM. Grasshoppers and crickets. Animal cytogenetics, vol. 3. Berlin: Insecta 1. Orthoptera. Gebrüder Borntraeger; 1979.
Castillo ER, Martí DA, Bidau CJ. Sex and neo-sex chromosomes in Orthoptera: a review. J Orthopt Res. 2010;19:213–31.
Palacios-Gimenez OM, Castillo ER, Martí DA, Cabral-de-Mello DC. Tracking the evolution of sex chromosome systems in Melanoplinae grasshoppers through chromosomal mapping of repetitive DNA sequences. BMC Evol Biol. 2013;3:167.
Palacios-Gimenez OM, Cabral-de-Mello DC. Repetitive DNA chromosomal organization in the cricket Cycloptiloides americanus: a case of the unusual X1X20 sex chromosome system in Orthoptera. Mol Gen Genomics. 2015;290:623–31.
Palacios-Gimenez OM, Marti DA, Cabral-de-Mello DC. Neo-sex chromosomes of Ronderosia bergii: insight into the evolution of sex chromosomes in grasshoppers. Chromosoma. 2015;124:353–65.
Palacios-Gimenez OM, Carvalho CR, Ferrari Soares FA, Cabral-de-Mello DC. Contrasting the chromosomal organization of repetitive DNAs in two Gryllidae crickets with highly divergent karyotypes. PLoS One. 2015;10:e0143540.
Jaarola M, Martin RH, Ashley T. Direct evidence for suppression of recombination within two pericentric inversions in humans: a new sperm-fish technique. Am J Hum Genet. 1998;63:218–24.
Lahn BT, Page DC. Four evolutionary strata on the human X chromosome. Nature. 1999;286:964–7.
Ross MT, Grafham DV, Coffey AJ, et al. The DNA sequence of the human X chromosome. Nature. 2005;434:325–37.
Kirkpatrick K. How and why chromosome inversions evolve. PLoS Biol. 2010;8:e1000501.
Ross J, Peichel CL. Molecular cytogenetics evidence of rearrangement of the Y chromosome of the threespine stickleback fish. Genetics. 2008;179:2173–82.
Díaz MO, Sáez FA. DNA synthesis in the neo-X neo-Y sex deter- mination system of Dichroplus bergii (Orthoptera: Acrididae). Chromosoma. 1968;24:10–6.
Cardoso H, Dutra A. The Neo-X Neo-Y sex pair in Acrididae, its structure and association. Chromosoma. 1979;70:323–36.
Cigliano MM. Ronderosia, a new genus of south American Melanoplinae (Orthoptera: Acrididae). J Orthoptera Res. 1999;6:1–19.
Chintauan-Marquier IC, Jordan S, Berthier P, et al. Evolutionary history and taxonomy of a short-horned grasshopper subfamily: the Melanoplinae (Orthoptera: Acrididae). Mol Phyl Evol. 2011;58:22–32.
Kejnovsky E, Hobza R, Cermak T, et al. The role of repetitive DNA in structure and evolution of sex chromosomes in plants. Heredity. 2009;102:533–41.
Steflova P, Tokan V, Vogel I, et al. Contrasting patterns of transposable element and satellite distribution on sex chromosomes (XY1Y2) in the dioecious plant Rumex acetosa. Genome Biol Evol. 2013;5:769–82.
Frönicke L, Scherthan H. Zoo-fluorescence in situ hybridization analysis of human and lndian muntjac karyotypes (Muntiacus muntjak vaginalis) reveals satellite DNA clusters at the margins of conserved syntenic segments. Chromosom Res. 1997;5:254–61.
Hartmann N, Scherthan H. Characterization of ancestral chromosome fusion points in the Indian muntjac deer. Chromosoma. 2004;112:213–20.
Webb GC, White MJD, Contreras N, Cheney J. Cytogenetics of the parthogenetic grasshopper Warramaba (formely Moraba) virgo and its bisexual relatives. IV. Chromosome banding studies. Chromosoma. 1978;67:309–39.
Sambrook J, Russel DW. Molecular cloning. A laboratory manual. 3rd ed. New York: Cold Spring Harbor Laboratory Press; 2001.
Andrews S. FastQC. A quality control tool for high throughput sequence data. 2012. Available at http://www.bioinformatics.babraham.ac.uk/projects/fastqc/.
Gordon A, Hannon GJ. Fastx-toolkit. FASTQ/A short-reads pre-processing tools. 2010. Available at http://hannonlab.cshl.edu/fastx_toolkit.
Novák P, Neumann P, Macas J. Graph-based clustering and characterization of repetitive sequences in next-generation sequencing data. BMC Bioinformatics. 2010;11:378.
Novák P, Neumann P, Pech J, et al. RepeatExplorer: a galaxy-based web server for genome-wide characterization of eukaryotic repetitive elements from next-generation sequence reads. Bioinformatics. 2013;29:792–3.
Benson G. Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 1999;27:573–80.
Junier T, Pagni M. Dotlet: diagonal plots in a web browser. Bioinformatics. 2000;16:178–9.
NCBI Resource Coordinators. Database resources of the National Center for biotechnology information. Nucleic Acids Res. 2016;44:D7–D19.
Kohany O, Gentles AJ, Hankus L, Jurka J. Annotation, submission and screening of repetitive elements in Repbase: RepbaseSubmitter and censor. BMC Bioinformatics. 2006;7:474.
Ruiz-Ruano FJ, López-León MD, Cabrero J, Camacho JPM. High-throughput analysis of the satellitome illuminates satellite DNA evolution. Sci Rep. 2016;6:28333.
Ruiz-Ruano FJ, Cabrero J, López-León MD, Camacho JPM. Satellite DNA content illuminates the ancestry of a supernumerary (B) chromosome. Chromosoma. 2016; https://doi.org/10.1007/s00412-016-0611-8.
Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792–7.
Tamura K, Peterson D, Stecher G, et al. MEGA: molecular evolutionary genetics using maximun likelihood, evolutionary distance, and maximun parsimony methods. Mol Biol Evol. 2011;28:2731–9.
Rozen S, Skaletsky H. Primer3 on the WWW for general users and for biologist programmers. In: Bioinformatics methods and protocols. New Jersey: Humana Press; 1999. p. 365-386.
Pinkel D, Straume T, Gray JW. Cytogenetic analysis using quantitative, high sensitivity, fluorescence hybridization. Proc Natl Acad Sci USA. 1986;83:2934–8.
Cabral-de-Mello DC, Moura RC, Martins C. Chromosomal mapping of repetitive DNAs in the beetle Dichotomius geminatus provides the first evidence for an association of 5S rRNA and histone H3 genes in insects, and repetitive DNA similarity between the B chromosome and A complement. Heredity. 2010;104:393–400.
Camacho JPM, Ruiz-Ruano FJ, Martín-Blazquéz R, et al. A step to the gigantic genome of the desert locust: chromosome sizes and repeated DNAs. Chromosoma. 2015;124:263–75.
Pezer Ž, Brajković J, Feliciello I, Ugarković Đ. Transcription of satDNAs in insects. In: Ugarković Ð, editor. Long non-coding RNAs, progress in molecular 51, and subcellular biology. Berlin-Heidelberg: Springer-Verlag; 2011. p. 161–78.
Pezer Ž, Ugarković Đ. Satellite DNA-associated siRNAs as mediators of heat shock response in insects. RNA Biol. 2012;9:587–95.
Plohl M, Meštrović N, Mravinac B. Centromere identity from DNA point of view. Chromosoma. 2014;123:313–25.
Palacios-Giemenez OM, Dias GB, Lima LG, Kuhn GCS, Ramos E, Martins C, Cabral-de-Mello DC. High-throughput analysis of the satellitome revealed enormous diversity of satellite DNAs in the neo-Y chromosome of the cricket Eneoptera surinamensis. Sci Rep. 2017;7:6422.
Song H, Amédégnato C, Cigliano MM, et al. 300 million years of diversification: elucidating the patterns of orthopteran evolution based on comprehensive taxon and gene sampling. Cladistics. 2015;31:621–51.
López-Flores I, Garrido-Ramos MA. The repetitive DNA content of eukaryotic genomes. In: Garrido-Ramos MA, editor. Repetitive DNA. Genome Dyn 7. Basel: Karger Publishers; 2012. p. 1–28.
Pavlek M, Gelfand Y, Plohl M, Meštrović N. Genome-wide analysis of tandem repeats in Tribolium castaneum genome reveals abundant and highly dynamic tandem repeat families with satellite DNA features in euchromatic chromosomal arms. DNA Res. 2015;22:387-01.
Mravinac B, Ugarković Đ, Franjević D, Plohl M. Long inversely oriented subunits form a complex monomer of Tribolium brevicornis satellite DNA. J Mol Evol. 2005;60:513–25.
Ugarković Đ, Podnar M, Plohl M. Satellite DNA of the red flour beetle Tribolium castaneum - comparative study of satellites from the genus Tribolium. Mol Biol Evol. 1996;13:1059–66.
Feliciello I, Chinali G, Ugarković Đ. Structure and evolutionary dynamics of the major satellite in the red flour beetle Tribolium castaneum. Genetica. 2011;139:999.
Wang X, Fang X, Yang P, et al. The locust genome provides insight into swarm formation and long-distance flight. Nat Commun. 2014;5:2957.
Lepesant JMJ, Cosseau C, Boissier J, et al. Chromatin structural changes around satellite repeats on the female sex chromosome in Schistosoma mansoni and their possible role in sex chromosome emergence. Genome Biol. 2012;13:R14.
Steinemann M, Steinemann S. The enigma of Y chromosome degeneration: TRAM, a novel retrotransposon is preferentially located on the neo-Y chromosome of Drosophila miranda. Genetics. 1997;145:261–6.
Steinemann S, Steinemann M. Retroelements: tools for sex chromosome evolution. Cytogenet Genome Res. 2005;110:134–43.
Matsunaga S. Junk DNA promotes sex chromosome evolution. Heredity. 2009;102:525–6.
Kaiser VB, Bachtrog D. Evolution of sex chromosome in insects. Annu Rev Genet. 2010;44:91–112.
Carvalho AB. Origin and evolution of the Drosophila Y chromosome. Curr Opin Genet Dev. 2002;12:664–8.
Vítková M, Fuková I, Kubíčková S, Marec F. Molecular divergence of the W chromosomes in pyralid moths (Lepidoptera). Chromosom Res. 2007;15:917–30.
Castillo ERD, Tafarel A, Martí DA. The early evolutionary history of neo-sex chromosomes in Neotropical grasshoppers, Boliviacris noroestensis (Orthoptera: Acrididae: Melanoplinae). Eur J Entomol. 2014;111:321–7.
Dover GA. Molecular drive: a cohesive mode of species evolution. Nature. 1982;299:111–7.
Walsh JB. Persistence of tandem arrays: implications for satellite and simple-sequence DNAs. Genetics. 1987;115:553Y567.
Hoffman AA, Rieseberg LH. Revisiting the impact of inversions in evolution: from population genetic markers to drivers of adaptive shifts and speciation. Annu Rev Ecol Evol Syst. 2008;39:21–42.
We thank “Parque Estadual Edmundo Navarro de Andrade” administration for sample collecting authorization and Antonio Sergio Pascon for excellent technical assistance in obtaining embryos.
This work was supported by FAPESP (process number 2014/11763-8) and Coordenadoria de Aperfeiçoamento de Pessoal de Nível Superior-CAPES. OMPG and DM were supported by Fundação de Amparo a Pesquisa do Estado de São Paulo-FAPESP (process numbers 2014/02038-8 and 2015/05246-3, respectively). DCCM was the recipient of a research productivity fellowship from the Conselho Nacional de Desenvolvimento Científico e Tecnológico-CNPq (process number 304758/2014-0). DAM and ERDC were supported by Argentina’s Consejo Nacional de Investigaciones Científicas y Técnicas - CONICET. The funders had no role in the study design, data collection and analysis, interpretation of data, or writing of the manuscript.
Availability of data and materials
DNA sequences have been deposited in GenBank under Accession numbers MF765804-MF765813. All remaining data supporting our findings is contained within the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
satDNAs monomer consensus sequences. (DOC 63 kb)
Table displaying the size ratio Y/X of the neo-Y variants in distinct mitotic metaphases. For measurement of sex chromosome sizes from ten metaphases it was used one individual for each variant. (DOCX 25 kb)
Distance between markers regarding the size of the neo-Y chromosome. In mitosis for variant II we studied five individuals and for variant III one individual. In meiosis for variant II we studied seven individuals, for variant III five individuals and for variant V one individual. For each individual five cells were analyzed. (DOCX 36 kb)