- Research article
- Open Access
The satellite DNA AflaSAT-1 in the A and B chromosomes of the grasshopper Abracris flavolineata
BMC Geneticsvolume 18, Article number: 81 (2017)
Satellite DNAs (satDNAs) are organized in repetitions directly contiguous to one another, forming long arrays and composing a large portion of eukaryote genomes. These sequences evolve according to the concerted evolution model, and homogenization of repeats is observed at the intragenomic level. Satellite DNAs are the primary component of heterochromatin, located primarily in centromeres and telomeres. Moreover, satDNA enrichment in specific chromosomes has been observed, such as in B chromosomes, that can provide clues about composition, origin and evolution of this chromosome. In this study, we isolated and characterized a satDNA in A and B chromosomes of Abracris flavolineata by integrating cytogenetic, molecular and genomics approaches at intra- and inter-population levels, with the aim to understand the evolution of satDNA and composition of B chromosomes.
AflaSAT-1 satDNA was shared with other species and in A. flavolineata, was associated with another satDNA, AflaSAT-2. Chromosomal mapping revealed centromeric blocks variable in size in almost all chromosomes (except pair 11) of A complement for both satDNAs, whereas for B chromosome, only a small centromeric signal occurred. In distinct populations, variable number of AflaSAT-1 chromosomal sites correlated with variability in copy number. Instead of such variability, low sequence diversity was observed in A complement, but monomers from B chromosome were more variable, presenting also exclusive mutations. AflaSAT-1 was transcribed in five tissues of adults in distinct life cycle phases.
The sharing of AflaSAT-1 with other species is consistent with the library hypothesis and indicates common origin in a common ancestor; however, AflaSAT-1 was highly amplified in the genome of A. flavolineata. At the population level, homogenization of repeats in distinct populations was documented, but dynamic expansion or elimination of repeats was also observed. Concerning the B chromosome, our data provided new information on the composition in A. flavolineata. Together with previous results, the sequences of heterochromatic nature were not likely highly amplified in the entire B chromosome. Finally, the constitutive transcriptional activity suggests a possible unknown functional role, which should be further investigated.
For more than half a century, researchers have attempted to understand the complexity and evolution of the repetitive DNA fraction of the genome. Repetitive DNAs are classified depending on the rate of repetitiveness and arrangement in the genome and include Transposable Elements (TEs) and minisatellite, microsatellite and satellite DNA (satDNA) [1, 2]. The satDNAs are a highly repetitive fraction of the genome that are organized by repetitions directly contiguous to one another, in tandem, forming long arrays with hundreds of copies generally composed of 100–1000 base pairs . These in tandem arrays compose most of the content of heterochromatin in eukaryotes, associated generally with the centromeres and telomeres [3,4,5,6]. These sequences evolve by the homogenization and fixation of different variants in a determined sexual population under the process of concerted evolution [7, 8]. Molecular mechanisms of DNA turnover, such as unequal crossing over and gene conversion, are primarily responsible for this homogenization pattern, leading to quantitative changes between species and generating specific satDNA subfamilies with differential arrangements and organization in a given genome [2, 9, 10]. For satDNA transcription, multiple functions are indicated to date , such as acting through RNAi as an epigenetic regulator of heterochromatin , playing a role as a structural element of centromeres, such as the alpha satellite in humans , and also regulating gene expression in yeast .
B chromosomes, which are extra elements to the standard complement (A complement) of some species, occur in approximately 15% of eukaryotes [15,16,17,18,19,20]. These chromosomes are known to be dispensable for normal development, for not recombining with standard A chromosomes and for their potential to present a drive mechanism and for accumulation . Very characteristically, B chromosomes frequently accumulate repetitive DNAs, including satDNAs. Because of the repetitive DNA content, for a long time, B chromosomes were considered genetically inert . However, studies revealed genes in B chromosomes [22,23,24,25,26,27] that could be transcriptionally active, such as in Crepis capillaris , Capreolus pygargus , Eyprepocnemis plorans , and Trichogramma kaykai . Moreover, sequences in B chromosomes can influence expression patterns from genes in A chromosomes [32,33,34]. These data suggest a putative biological role for B chromosomes in some species.
Although frequently reported in some groups, information on satDNA populating B chromosomes in insects is restricted to a few species, such as the grasshoppers Eyprepocnemis plorans [15, 35, 36] and Eumigus monticola , wasps [38, 39], Glossina  and Drosophila subsilvestris . In these B chromosomes, a direct association is observed between heterochromatin and accumulation of repetitive DNAs. In the organisms listed above, the isolated satDNAs helped in tracking the origin and evolution of B chromosomes.
In the grasshopper Abracris flavolineata with karyotype composed of 2n = 23, ×0 (males), one or two submetacentric B chromosomes were detected in individuals from a population sampled in Rio Claro/SP-Brazil . Some studies address B chromosome origin and evolution in this species using repetitive DNAs as markers [24, 42, 43], revealing clues about B chromosome composition and origin. Here, we used restriction enzymatic digestion to isolate the first satDNA in the species, named AflaSAT-1, with the aim to further understand satDNAs among grasshoppers at the population level and their contribution to B chromosome structure and evolution in A. flavolineata. Finally, AflaSAT-1 was studied using a combination of chromosomal and molecular analyses at the population level to uncover the genome organization and evolution in A and B chromosomes.
Isolation and characterization of satDNA sequences
The genomic DNA from individuals of the Rio Claro/SP population was digested using three different restriction endonucleases (RE), i.e., HindII, AluI and SmaI. A ladder pattern was revealed for HindII and SmaI enzymes with bands of approximately 150, 300, and 450 bp, corresponding to the monomer and multimers of a putative tandem sequence. The monomer generated by HindII was selected for further analysis (Fig. 1a). The monomer generated a sequence of 173 bp, after cloning and sequencing analysis. A convergent internal primer was designed and recovered a PCR fragment containing 137 bp, which formed a ladder pattern, typical for satDNAs (Fig. 1b, c). This sequence was named AflaSAT-1. Internal restriction sites for other enzymes were also recognized in AflaSAT-1 (Fig. 1c).
With PCR, AflaSAT-1 was detected in individuals in six populations that were sampled and in the μB-DNA from the Rio Claro/SP population. A total of 79 AflaSAT-1 sequences were recovered after cloning of individuals from Rio Claro/SP (15 clones), μB-DNA from Rio Claro/SP (10 clones), Cabo/PE (11 clones), Santa Bárbara do Pará/PA (10 clones), Paranaíta/MT (11 clones), Manaus/AM (12 clones), and Posadas/AR (10 clones). The AflaSAT-1 was similar in composition in all populations, with approximately 41.0% G + C base pairs. Different numbers of haplotypes were recognized depending on the population. Considering all sequences, 11 haplotypes were recognized. The number of mutations and variable sites was slightly higher in the sequences recovered from μB-DNA than in those from A chromosomes of the six populations. Nucleotide and haplotype diversity was also slightly variable, reaching the highest values in the Cabo/PE population (Table 1).
A minimum spanning tree revealed the relationships among the 11 haplotypes from all populations (Fig. 2). The tree of the 11 haplotypes was shared among at least five populations (Hap1, Hap3 and Hap5). The difference among the haplotypes was only one mutation (substitution), except for the haplotypes from exclusive μB-DNA that were differentiated by three or four mutations, Hap4 and Hap11. Finally, the μB-DNA presented a sequence grouped to Hap2, which occurred in Cabo/PE and Manaus/AM populations (Fig. 2).
Using the sequenced genome from one individual harboring one B chromosome, we identified one cluster that was similar to AflaSAT-1 in the assembly of RepeatExplorer. This cluster represented 1.04% of the genome, corresponding to the fourth most representative sequence. The cluster retrieved a total of 32 contigs, but only 19 entire sequences of AflaSAT-1 could be recovered, in addition to truncated sequences. Unexpectedly, in a few contigs and occurring with AflaSAT-1, a different sequence containing 242 bp was observed, which was named AflaSAT-2. For this repeat, we recovered ten entire sequences and some truncated repeats. Two arrangements were observed in the recovered contigs, only the AflaSAT-1 and AflaSAT-1 plus AflaSAT-2; the AflaSAT-2 was never recovered alone.
Chromosome distribution and copy number variation of AflaSAT-1 in populations of A. flavolineata
The FISH using AflaSAT-1 and AflaSAT-2 probes in individuals from Rio Claro/SP revealed large pericentromeric blocks in all chromosomes, with blocks varying in size, except for pair 11 that did not have visible signals (Fig. 3a-d). The B chromosome showed a small signal in the centromeric region (Fig. 3e, f). These two sequences were interglimed with one another, as observed in fiber-FISH experiments (Fig. 3g). Association of AflaSAT-1 and AflaSAT-2 forming a composed unit was also observed in sequences retrieved from sequenced genomes.
Analysis of chromosomal distribution of AflaSAT-1 in different populations revealed identical patterns for individuals from Rio Claro/SP and from Sta. Bárbara do Pará/PA (Fig. 4a). In the other populations, a slight variation was observed, i.e., in Cabo/PE, the pair 1 did not have a signal and the pair 11 had a small positive signal (Fig. 4b), whereas individuals from Posadas/AR had signals in all chromosomes (Fig. 4c).
Additionally, we observed copy number variation of AflaSAT-1 with qPCR. Males without and with one B chromosome from Rio Claro/SP showed differences, with the individuals harboring the B chromosome having a higher copy number, as expected. Females from Rio Claro/SP also had a higher copy number, which was attributed to two copies of the X chromosome. Compared with those of Rio Claro/SP, individuals from Cabo/PE had fewer copies, whereas individuals from Santa Bárbara do Pará/PA had similar copy numbers. Finally, individuals from Posadas/AR had the highest number of copies (Fig. 5). The high standard deviations, considering the absolute copy number quantified, indicated that the copy number of AflaSAT-1 is very unstable and variable individually (Additional file 1).
Occurrence of AflaSAT-1 in other grasshopper species
In the dot blot test for AflaSAT-1, among the six species tested, none revealed positive signals, and positive hybridization was restricted to genomic DNA from A. flavolineata. In the two species of the Ommatolampidinae, the occurrence of AflaSAT-1 was also tested using FISH, and no signals were observed (data not shown). However, with PCR amplification, the AflaSAT-1 monomer was detected in A. dilecta, V. rugulosa and R. bergi and further confirmed through sequencing of PCR products for each species (Additional file 2).
Transcriptional activity of AflaSAT-1
Based on RT-PCR analysis, AflaSAT-1 was transcriptionally active in all different organ tissues tested in male and female adults, with activity also in both embryos and nymphs. The different tissues included the head, saltatory leg, testis, ovariole and gastric caecum, in addition to embryos and first stage nymphs. These diverse samples were studied to test for constitutional transcription of AflaSAT-1, which could suggest putative function for the organism. Tissues from somatic and germ line lineages were also studied to reinforce putative constitutional function. The 1B male individuals also showed transcriptional activity for AflaSAT-1. For the analyses of all tissues and life stages, the most evident bands corresponded to monomers and tetramers of AflaSAT-1 DNA (Fig. 6), and the sequencing of monomers from randomly selected samples confirmed the specific transcriptional activity of AflaSAT-1.
AflaSAT-1 organization and evolution in a chromosomes of Abracris flavolineata
satDNA families have been described for only 12 grasshopper species, most members of Acrididae. The representation of AflaSAT-1 was low (1.04%) in the A. flavolineata genome but corresponded to the fourth most abundant repetitive DNA (Milani et al., unpublished data). In grasshoppers, different satDNA repeats commonly represent less than 1% of the genome, as observed L. migratoria and E. monticola [37, 44]. However, highly abundant sequences were observed in Schistocerca gregaria with different satDNA families representing approximately 18% (SG1) and 12% (SG2-alfa) . As expected, the general chromosomal distribution of AflaSAT-1 was coincident with pericentromeric heterochromatin, a common pattern among eukaryotes and observed extensively in grasshoppers [35, 37, 44,45,46]. The occurrence of another satDNA, i.e., AflaSAT-2, with the same distribution as that of AflaSAT-1 indicated a more complex structure, which was demonstrated by the two types of arrangement in the A. flavolineata genome. This arrangement is similar to that described in Drosophila buzzatii for pBuM-1 and pBuM-2 satDNA subfamilies . The similar distributions of non-homologous satDNAs observed in this study have also been reported in some Drosophila , rye  and rodents .
Based on interspecies analysis, AflaSAT-1 was not specific to A. flavolineata and was recovered from the genomes of three other species with PCR, although dot blot analysis did not reveal positive hybridization. Although not species-specific, the copy number of AflaSAT-1 in the genomes of A. dilecta, V. rugulosa and R. bergi is probably low, whereas AflaSAT-1 was highly amplified in the genome of A. flavolineata, generating a clustered organization. Divergence among sequences among the different species was discarded as an explanation for non-positive dot blot signals, because the sequencing of PCR products revealed almost no mutations among sequences. The sharing of AflaSAT-1 among different species is consistent with the “library hypothesis” by which related species share satDNA families that can be stochastically amplified or not in the diverging genomes . Considering R. bergi, as a member of Melanoplinae, and A. flavolineata, A. dilecta and V. rugulosa, as members of Ommatolampidinae and the sharing of AflaSAT-1, the origin of this satDNA likely occurred in a common ancestor before the divergence of the two subfamilies, which remits to approximately less than 73 mya, corresponding to the time of origin of Acrididae . Recently, the most comprehensive comparison between the satellitomes of two other grasshopper species, i.e., L. migratoria (Acrididae) and E. monticola (Pamphagidae), showed that the life span of a satDNA library is less than approximately 100 mya, corresponding to the most recent common ancestor between the two families . Our result was within this life span and showed that the divergence between satDNA libraries is variable in Acrididae, considering that the AflaSAT-1 was found in some species, depending on the subfamily, but not in others.
Low variability was observed for AflaSAT-1, and additionally, almost no distinct or population-specific profile was observed, which has also been documented for other species consistent with the concerted evolution model for homogenization of repeats at the intragenomic level [53,54,55,56,57]. By contrast, chromosomal variability was observed for cluster number, size, distribution and abundance of repeats, depending on the population. These data indicate that a random expansion or elimination of the satDNA may occur, generating different patterns in different populations. Variation related to satDNA mapping was also found among populations of the frog Physalaemus cuvieri  and in the grasshopper Eyprepocnemis plorans , suggesting a high rate of amplification or deletion mutational events. As expected, the population-dependent higher or lower number of chromosomal clusters for AflaSAT-1 was also consistent with the increase or decrease of AflaSAT-1 copy number, corroborating the dynamic of amplification and deletion of repeats. This change in copy number could be induced by several events, such as unequal crossing-over, slippage replication, extrachromosomal circular DNA and rolling-circle, which are common for repetitive DNAs [3, 60, 61].
AflaSAT-1 and its relationship with the B chromosome
To date, satDNA in B chromosomes has been reported in few species of grasshoppers, such as Eyprepocnemis plorans  and Eimugus monticola , which revealed important characteristics about the origin and evolution of B chromosomes. In A. flavolineata, AflaSAT-1 did not disclose clues concerning the specific origin of the B chromosome because of the occurrence in almost all centromeres of A complement. However, AflaSAT-1 helped to understand the composition and molecular evolution of the B chromosome. The C-positive heterochromatin blocks are absent in the B chromosome of A. flavolineata [24, 42], and in the chromosomal arms, different sequences are shared with the euchromatic regions of A complement. It is considered remarkable for mobile elements and anonymous sequences to be obtained from a microdissected B chromosome [43, 61]. The heterochromatin-enriched sequences of A complement that occur in the B chromosome are restricted to a small signal in the centromeric region, as described in this study for AflaSAT-1 and for the C 0 t-1 DNA . Together, these data suggest that the heterochromatic sequences were not highly amplified in the B chromosome of A. flavolineata, as is frequently observed in other species , and that this chromosome most likely bears unknown euchromatic sequences. Hence, genes and other satDNAs should be investigated for a clearer picture regarding this scenario.
Molecular sequence analysis of AflaSAT-1 recovered from the microdissected B chromosome revealed exclusive mutations, with a higher number of variable sites and mutations in this chromosome. The higher sequence variability of B chromosome than that of A complement is commonly reported and has been attributed to the higher tolerance of mutations because of the dispensable nature of chromosome B . Similar data are also observed for other types of repetitive DNAs, such as histone genes and 18S rRNA gene in the grasshopper L. migratoria  and in the fish Astyanax paranae . However, in E. plorans, two different sequences (satDNA and 45S rDNA), and in A. flavolineata, the U2 snDNA, revealed different homogenization patterns with high similarity between A and B chromosomes [61, 64]. The discrepancy in sequence variability observed for the two repetitive DNAs in the B chromosome of A. flavolineata, i.e., AflaSAT-1 and U2 snDNA, suggests that different sequences with variable roles have divergent evolutionary patterns.
AflaSAT-1 is transcriptionally active satellite sequence
SatDNA transcripts have been reported in different eukaryotes, and the evidence is accumulating that satellite transcription might be a common placement with putative structural or functional role for genomes . Included among the important roles of satDNAs are heterochromatin formation and regulation, involvement in centromere function, epigenetic chromatin silencing and modulation, and regulation of genes, among others .
Among insects, satDNA transcription has been described generally in Hymenoptera, Orthoptera, Diptera and Coleoptera  and in the lepidopteran species Cydia pomonella  and Plodia interpunctella . Differential transcription of satDNAs in this group is related to different developmental stages and tissues and also the sexes [54, 69, 70]. Considering this information and the constitutional transcription of AflaSAt-1, we suggest an unknown biological function for AflaSAT-1 satDNA, such as a regulatory element, with a structural or functional role, in the A. flavolineata genome. In future research, the quantitative differential transcription of AflaSAT-1 in A. flavolineata, depending on tissue, sex and presence/absence of B chromosome, must be determined to obtain more precise information for putative functional roles.
A total of 129 adult A. flavolineata were collected at five different sites in Brazil (BR) and one from Argentina (AR): Rio Claro/SP (São Paulo), 22°24′45″ S, 47°31′28″ W (30 individuals); Cabo/PE (Pernambuco), 8°17′15″ S, 35°2′7″ W (19 individuals); Sta. Bárbara do Pará/PA (Pará), 1°13′27″ S, 48°17′38″ W (60 individuals); Paranaíta/MT (Mato Grosso), 9°40′25″ S, 56°28′36″ W (5 individuals); Manaus/AM (Amazonas), 3°06′02″ S, 59°58′31″ W (4 individuals); and Posadas/Misiones/AR, 27o25’ S, 55o56’ W (14 individuals). Testes were fixed in 3:1 absolute ethanol:acetic acid and stored at −20 °C. Entire animals were immersed in absolute ethanol and stored at −20 °C for genomic DNA (gDNA) extraction. We obtained tissues for DNA extraction and for chromosomes from Rio Claro/SP, Cabo/PE, Sta Bárbara do Pará/PA and Posadas/AR populations, whereas from Manaus/AM and Paranaíta/MT populations, only tissues for DNA extraction were obtained. For each analysis, at least three individuals were used.
To obtain embryos, animals collected in Rio Claro/SP were placed in plastic boxes until oviposition. Embryo mitotic cells were obtained from embryo dissection approximately 15 days after deposition, following  method.
Restriction enzymatic digestion, cloning and primer design
The genomic DNA was extracted using the phenol-chloroform method  and stored at −20 °C until use. Genomic DNA of the specimens from the Rio Claro/SP population was fragmented by restriction enzymatic digestion using HindII, AluI and SmaI enzymes. The digested products were fractionated by electrophoresis in 1% agarose gel to verify a ladder pattern, typical for in tandem DNA sequences. Fragments with a ladder pattern were selected and purified using a Zymoclean™ Gel DNA Recovery Kit (Zymo Resarch Corp., The Epigenetics Company, USA) according to the manufacturer’s instructions.
The purified products were cloned using a pMOSBlue Blunt Ended Cloning Kit (GE Healthcare) with DH5α Escherichia coli as competent cells. Positive colonies were randomly chosen and screened by Polymerase Chain Reaction (PCR) using M13 primers set (F-5’GTAAAACGACGGCCAG and R-5’CAGGAAACAGCTATGAC) for DNA sequencing by Macrogen Inc. (Korea). Geneious v4.8.5 software  was used to check quality and exclude vector sequences. To amplify the satDNA sequence (AflaSAT-1) using PCR, the specific primer set (F-5’GACAGTTTTAAACACTTCCATTACAG and R-5’GACTGTGTTGATATCCAATAACA)was designated.
PCR amplification and sequence analysis
Using specific primers for AflaSAT-1, PCR amplification was conducted using the genomic DNA as template of animals from the six different populations and from the B chromosome previously microdissected (μB-DNA) from individuals from Rio Claro/SP . PCR products were visualized on a 1% agarose gel, and the bands were isolated and purified using a ZymocleanTM Gel DNA Recovery Kit (Zymo Research Corp., The Epigenetics Company, USA) according to the producer’s recommendations.
For cloning of the purified PCR products, pGEM-T easy vector (Promega, Madison, WI, USA) was used with DH5α Escherichia coli as competent cells. Positive clones were screened using the M13 primers and sequenced by Macrogen Inc. (Korea).
The monomer sequences were analyzed using Geneious v4.8.5 , and then, DNA polymorphism and haplotype recognition were checked using the DnaSP v.5.10.01 tool . For more accurate analysis, singleton sequences were discarded. A graphical haplotype network was constructed employing Network 184.108.40.206 software (http://www.fluxus-engineering.com). The correspondent haplotype sequences were deposited in NCBI database under the accession numbers MF752447-MF752457.
We searched the AflaSAT-1 repeats in the sequenced genome of A. flavolineata from an individual harboring one B chromosome. This genome was sequenced using Illumina Miseq paired-end (2 × 300), and the libraries were constructed using a Nextera DNA Library Preparation Kit and quantified by a KAPA Library Quantification Kit (Milani et al., in preparation). For these purposes, we applied the graph-based clustering and assembly using RepeatExplorer  and manually searched for sequences similar to AflaSAT-1 in the assembled contigs using Geneious v4.8.5 . For this search, we considered only clusters showing high graph density that are typical for satDNA . In the clusters containing the AflaSAT-1, we also observed another associated sequence (named AflaSAT-2), and a primer set to recover this sequence was designed, F-5’GGGTCTCGCGAAATGAGAC and R-5’GCTTTCTAAACGGAATCGAG.
Chromosome obtaining and fluorescent in situ hybridization (FISH)
Meiotic cells were obtained from testes, whereas mitotic cells were obtained from embryos. For conventional analyses used to check the general chromosomal structure of animals from different populations, the slides were prepared with tissue maceration and staining with Giemsa 5%.
The PCR products from satDNAs from Rio Claro/SP individuals were used as probes for FISH assays using chromosomes of individuals from the populations of Rio Claro/SP, Cabo/PE, Sta Bárbara do Pará/PA and Posadas/AR. Fragments were labeled by nick translation using digoxigenin-11-dUTP and detected by anti-digoxigenin rhodamine (Roche, Mannheim, Germany) or biotin-14-dUTP detected with streptavidin Alexa Fluor-488 conjugated (Invitrogen, San Diego, CA, USA). FISH experiments were conducted following , with adaptations by . Fiber-FISH was conducted according to . Slides were counterstained with DAPI (4′,6-Diamidine-2′-phenylindole) and mounted with VECTASHIELD (Vector, Burlingame, CA, USA). Pictures were captured using a DP70 cooled digital camera in gray scale coupled with an Olympus microscope BX51 equipped with a fluorescence lamp and appropriate filters. The images were pseudo-colored, merged and treated for brightness and contrast using Adobe Photoshop CS6.
Dot blot hybridization
To determine whether the AflaSAT-1 was shared with other species at the level of genus, subfamily and family, we conducted a dot blot experiment, following the descriptions from , using the sequence amplified from individuals from Rio Claro/SP as a probe for the other six species. The probe was tested against the genomic DNA of A. flavolineata as positive control, the congeneric species A. dilecta and one another species in the same subfamily (Ommatolampidinae), i.e., Vilerna rugulosa. Moreover, we used species from other subfamilies of Acrididae, including Melanoplinae (Ronderosia bergi), Cyrtacanthacridinae (Schistocerca pallens), Gomphocerinae (Amblytropidia robusta) and Leptysminae (Eumastusia koebelei koebelei). The presence of AflaSAT-1 was also tested using PCR in each species following the same conditions applied for A. flavolineata described above, and the results were checked by nucleotide sequencing.
Copy number detection by qPCR
To verify the copy number of AflaSAT-1 among different populations and in individuals harboring one B chromosome, we performed a relative quantification employing the gene dose ΔCt method . qPCR assays were performed using the gDNA from male samples of each population for which we obtained FISH results, i.e., Cabo/PE, Sta Bárbara do Pará/PA and Posadas/AR. For the Rio Claro/SP population, we tested males with and without one B chromosome and females without a B chromosome. The gene dosage ratios (GDR) were obtained following the same parameter of , using a Heat Shock Protein (Hsp70) as the reference gene (F-5’GGTGTGATGACCACTCTTATCAA and R-5’CACTTCAATTTGAGGCACACC), because no difference was detected in amplification ratio among males, females and 0B and 1B individuals. Our results support that the Hsp70 gene is placed in autosomes, considering no difference between males and females. The target and reference gene were analyzed simultaneously in duplicates of four independent samples, and the qPCR conditions were set at 95 °C for 10 min; 45 cycles of 95 °C for 15 s, and 60 °C for 1 min, performed in a StepOne Real-Time PCR System (Life Technologies, Carlsbad, CA). Specificity of the PCR products was confirmed by analysis of the dissociation curve.
RNA extraction and transcriptional analyses
The RNA was extracted from organ tissues of male and female adults, including head, saltatory leg, testis, ovariole and gastric caecum, and from embryos (with 20 days of development) and first stage nymphs (immediately after hatch) from the Rio Claro/SP population. We also tested these tissues from adult males with 1B. RNA was extracted using TRIzol® Reagent (Life Technologies) and then treated with Amplification Grade DNAse I (Sigma-Aldrich) to avoid DNA contaminants. Finally, the cDNA samples were obtained by reverse transcription (RT-PCR) using a High-Capacity cDNA Reverse Transcription Kit (Life Technologies) with the reaction conditions set to 25 °C for 10 min; 37 °C for 120 min, and 85 °C for 5 min. Afterward, the cDNA samples were used as templates for conventional PCR amplification of AflaSAT-1, with at least three biological replicates.
To validate the signals of transcription and to ensure that the samples were not contaminated with genomic DNA, the cDNA samples were used as source for amplification of Hsp70 as a control. The amplification of the Hsp70 gene in genomic DNA samples and uncontaminated cDNA samples revealed bands with different sizes of approximately 300 and 200 bp, respectively. An intron in the genomic DNA that was not present after RNA transcription and processing caused the difference. Some PCR products were sequenced to confirm the HSP70 sequence.
Fluorescence in situ hybridization
Polymerase chain reaction
Biémont C, Vieira C. Genetics: junk DNA as an evolutionary force. Nature. 2006;443:521–4.
Plohl M, Meštrović N, Mravinac B. Satellite DNA evolution. In: Garrido-Ramos MA, editor. Repetitive DNA. Genome dynamics. 7th ed. Basel: Karger; 2012. p. 126–52.
Charlesworth B, Sniegowski P, Stephan W. The evolutionary dynamics of repetitive DNA in eukaryotes. Nature. 1994;371:215–20.
John B, Miklos GL. Functional aspects of satellite DNA and heterochromatin. In: Bourne GH, Danielli JF, editors. International review of cytology, vol. 58. New York: Academic Press; 1979. p. 1–114.
Ugarković D, Plohl M. Variation in satellite DNA profiles-causes and effects. EMBO J. 2002;21:5955–9.
Palomeque T, Lorite P. Satellite DNA in insects: a review. Heredity. 2008;100:564–73.
Dover GA. Molecular drive in multigene families: how biological novelties arise, spread and are assimilated. Trends Genet. 1986;2:159–65.
Dover GA. Molecular drive. Trends Genet. 2002;18:587–9.
Pons J, Petitpierre E, Juan C. Evolutionary dynamics of satellite DNA family PIM357 in species of the genus Pimelia (Tenebrionidae, Coleoptera). Mol Biol Evol. 2002;19:1329–40.
Plohl M, Luchetti A, Meštrović N, Mantovani B. Satellite DNAs between selfishness and functionality: structure, genomics and evolution of tandem repeats in centromeric (hetero) chromatin. Gene. 2008;409:72–82.
Pezer Z, Brajković J, Feliciello I, Ugarković D. Satellite DNA-mediated effects on genome regulation. Genome Dyn. 2012;7:153–69.
Volpe TA, Kidner C, Hall IM, Teng G, Grewal SIS, Martienssen RA. Regulation of heterochromatic silencing and histone H3 lysine-9 methylation by RNAi. Science. 2002;297:1833–7.
Wong LH, Brettingham-Moore KH, Chan L, Quach JM, Anderson MA, Northrop EL, Hannan R, Saffery R, Shaw ML, Williams E, Choo KH. Centromere RNA is a key component for the assembly of nucleoproteins at the nucleolus and centromere. Genome Res. 2007;17:1146–60.
Vinces MD, Legendre M, Caldara M, Hagihara M, Verstrepen KJ. Unstable tandem repeats in promoters confer transcriptional evolvability. Science. 2009;324:1213–6.
López-León MD, Neves N, Schwarzacher T, Heslop-Harrison JS, Hewitt GM, Camacho JPM. Possible origin of a B chromosome deduced from its DNA composition using double FISH technique. Chromosom Res. 1994;2:87–92.
de Jesus MC, Galetti PM, Valentini SR, Moreira-Filho O. Molecular characterization and chromosomal localization of two families of satellite DNA in Prochilodus lineatus (Pisces, Prochilodontidae), a species with B chromosomes. Genetica. 2003;118:25–32.
Camacho JPM. B chromosomes. In: Gregory TR, editor. The evolution of the genome. Amsterdam: Elsevier; 2005. p. 223–86.
Houben A, Banaei-Moghaddam AM, Klemme S, Timmis JN. Evolution and biology of supernumerary B chromosomes. Cell and Mol Lif Sci. 2014;71:467–78.
Klemme S, Banaei-Moghaddam AM, Macas J, Wicker T, Novák P, Houben A. High-copy sequences reveal distinct evolution of the rye B chromosome. New Phytol. 2013;199:550–8.
Kumke K, Macas J, Fuchs J, Altschmied L, Kour J, Dhar MK, Houben A. Plantago lagopus B chromosome is enriched in 5S rDNA-derived satellite DNA. Cytogenet Genome Res. 2016;148:68–73.
Camacho JPM, Sharbel TF, Beukeboom LW. B chromosome evolution. Phil Trans R Soc Lond B. 2000;355:163–78.
Teruel M, Cabrero J, Perfectti F, Camacho JP. B chromosome ancestry revealed by histone genes in the migratory locust. Chromosoma. 2010;119:217–25.
Yoshida K, Terai Y, Mizoiri S, Aibara M, Nishihara H, Watanabe M, Kuroiwa A, Hirai H, Hirai Y, Matsuda Y, Okada N. B chromosomes have a functional effect on female sex determination in Lake Victoria cichlid fishes. PLoS Genet. 2011;7:8e1002203.
Bueno D, Palacios-Gimenez OM, Cabral-de-Mello DC. Chromosomal mapping of repetitive DNAs in Abracris flavolineata reveal possible ancestry for the B chromosome and surprisingly H3 histone spreading. PLoS One. 2013;8:e66532.
Valente GT, Conte MA, Fantinatti BE, Cabral-de-Mello DC, Carvalho RF, Vicari MR, Kocher TD, Martins C. Origin and evolution of B chromosomes in the cichlid fish Astatotilapia latifasciata based on integrated genomic analyses. Mol Biol Evol. 2014;31:2061–72.
Makunin AI, Dementyeva PV, Graphodatsky AS, Volobouev VT, Kukekova AV, Trifonov VA. Genes on B chromosomes of vertebrates. Mol Cytogenetics. 2014;7:99.
Makunin AI, Kichigin IG, Larkin DM, O'Brien PC, Ferguson-Smith MA, Yang F, Proskuryakova AA, Vorobieva NV, Chernyaeva EN, O'Brien SJ, Graphodatsky AS, Trifonov VA. Contrasting origin of B chromosomes in two cervids (Siberian roe deer and grey brocket deer) unravelled by chromosome-specific DNA sequencing. BMC Genomics. 2016;17:618.
Leach CR, Houben A, Field B, Pistrick K, Demidov D, Timmis JN. Molecular evidence for transcription of genes on a B chromosome in Crepis capillaris. Genetics. 2005;171:269–78.
Trifonov VA, Dementyeva PV, Larkin DM, O'Brien PC, Perelman PL, Yang F, Ferguson-Smith MA, Graphodatsky AS. Transcription of a protein-coding gene on B chromosomes of the Siberian roe deer (Capreolus pygargus). BMC Biol. 2013;11:90.
Ruiz-Estevez M, Lopez-Leon MD, Cabrero J, Camacho JPM. B-chromosome riboso- mal DNA is functional in the grasshopper Eyprepocnemis plorans. PLoS One. 2012;7:e36600.
Van Vugt JJFA, Nooijer S, Stouthamer R, Jong H. NOR activity and repeat sequences of the paternal sex ratio chromosome of the parasitoid wasp Trichogramma kaykai. Chromosoma. 2005;114:410–9.
Carchilan M, Kumke K, Mikolajewski S, Houben A. Rye B chromosomes are weakly transcribed and might alter the transcriptional activity of a chromosome sequences. Chromosoma. 2009;118:607–16.
Banaei-Moghaddam AM, Meier K, Karimi-Ashtiyani R, Houben A. Formation and expression of pseudogenes on the B chromosome of rye. Plant Cell. 2013;25:2536–44.
Huang W, Du Y, Zhao X, Jin W. B chromosome contains active genes and impacts the transcription of a chromosomes in maize (Zea mays L.). BMC Plant Biol. 2016;16:1.
López-León MD, Vázquez P, Hewitt GM, Camacho JPM. Cloning and sequence analysis of an extremely homogeneous tandemly repeated DNA in the grasshopper Eyprepocnemis plorans. Heredity. 1995;75:370–5.
Cabrero J, Bakkali M, Bugrov A, Warchalowska-Sliwa E, López-León MD, Perfectti F, Camacho JPM. Multiregional origin of B chromosomes in the grasshopper Eyprepocnemis plorans. Chromosoma. 2003;112:207–11.
Ruiz-Ruano FJ, Cabrero J, López-León MD, Camacho JPM. Satellite DNA content illuminates the ancestry of a supernumerary (B) chromosome. Chromosoma. 2017;126:487–500.
Nur U, Werren JH, Eickbush DG, Burke WD, Eickbush TH. A “selfish” B chromosome that enhances its transmission by eliminating the paternal genome. Science. 1988;240:512–4.
Eickbush DG, Eickbush TH, Werren JH. Molecular characterization of repetitive DNA sequences from a B chromosome. Chromosoma. 1992;101:575–83.
Amos A, Dover G. The distribution of repetitive DNAs between regular and supernumerary chromosomes in species of Glossina (tsetse): a two-step process in the origin of supernumeraries. Chromosoma. 1981;81:673–90.
Gutknecht J, Sperlich D, Bachmann L. A species specific satellite DNA family of Drosophila subsilvestris appearing predominantly in B chromosomes. Chromosoma. 1994;103:539–44.
Milani D, Cabral-de-Mello DC. Microsatellite organization in the grasshopper Abracris flavolineata (Orthoptera: Acrididae) revealed by FISH mapping: remarkable spreading in the a and B chromosomes. PLoS One. 2014;9:e97956.
Palacios-Gimenez OM, Bueno D, Cabral-de-Mello DC. Chromosomal mapping of two Mariner-like elements in the grasshopper Abracris flavolineata (orthoptera: Acrididae) reveals enrichment in euchromatin. Eur J Entomol. 2014;111:329–34.
Ruiz-Ruano FJ, López-León MD, Cabrero J, Camacho JPM. High-throughput analysis of the satellitome illuminates satellite DNA evolution. Scientific Reports. 2016;6:28333.
Camacho JM, Ruiz-Ruano FJ, Martín-Blázquez R, López-León MD, Cabrero J, Lorite P, Cabral-de-Mello DC, Bakkali M. A step to the gigantic genome of the desert locust: chromosome sizes and repeated DNAs. Chromosoma. 2015;124:263–75.
Rodríguez IE, Fernández-Calvín B, Capel J, La Vega D, García C. Equilocality and heterogeneity of constitutive heterochromatin: in situ localization of two families of highly repetitive DNA in Dociostaurus genei (Orthoptera). Heredity. 1996;76:70–6.
Kuhn GC, Sene FM, Moreira-Filho O, Schwarzacher T, Heslop-Harrison JS. Sequence analysis, chromosomal distribution and long-range organization show that rapid turnover of new and old pBuM satellite DNA repeats leads to different patterns of variation in seven species of the Drosophila buzzatii cluster. Chromosom Res. 2008;16:307–24.
Kuhn GCS, Teo CH, Schwarzacher T, Heslop-Harrison JS. Evolutionary dynamics and sites of illegitimate recombination revealed in the interspersion and sequence junctions of two nonhomologous satellite DNAs in cactophilic Drosophila species. Heredity. 2009;102:453–64.
Alkhimova OG, Mazurok NA, Potapova TA, Zakian SM, Heslop-Harrison JS, Vershinin AV. Diverse patterns of the tandem repeats organization in rye chromosomes. Chromosoma. 2004;113:42–52.
Paço A, Adega F, Meštrović N, Plohl M, Chaves R. The puzzling character of repetitive DNA in Phodopus genomes (Cricetidae, Rodentia). Chromosom Res. 2015;23:427–40.
Fry K, Salser W. Nucleotide sequences of HS-alfa satellite DNA from kangaroo rat Dipodomys ordii and characterisation of similar sequences in other rodents. Cell. 1977;12:1069–84.
Song H, Amédégnato C, Cigliano MM, Desutter-Grandcolas L, Heads SW, Huang Y, Otte D, Whiting MF. 300 million years of diversification: elucidating the patterns of orthopteran evolution based on comprehensive taxon and gene sampling. Cladistics. 2015;36:621–51.
Bachmann L, Venanzetti F, Sbordoni V. Characterization of a species-specific satellite DNA of Dolichopoda schiavazzi (Orthoptera, Rhaphidophoridae) cave crickets. J Mol Evol. 1994;39:274–81.
Lorite P, Carrillo JA, Tinaut A, Palomeque T. Comparative study of satellite DNA in ants of the Messor genus (hymenoptera, Formicidae). Gene. 2002;297:113–22.
Kuhn GCS, Sene FM. Characterization and interpopulation variability of a complex HpaI satellite DNA of Drosophila seriema (repleta group). Genetica. 2004;121:241–9.
Feliciello I, Chinali G, Ugarković D. Structure and evolutionary dynamics of the major satellite in the red flour beetle Tribolium castaneum. Genetica. 2011;139:999–1008.
Vieira P, Castagnone C, Mallez S, Espada M, Navas A, Mota M, Castagnone-Sereno P. Sequence variability of the MspI satellite DNA family of the pinewood nematode Bursaphelenchus xylophilus at different geographic scales. Mol Phylogenet Evol. 2014;70:120–9.
Vittorazzi SE, Lourenço LB, Recco-Pimentel SM. Long-time evolution and highly dynamic satellite DNA in leptodactylid and hylodid frogs. BMC Genet. 2014;15:1–11.
Cabrero J, López-León MD, Ruíz-Estévez M, Gómez R, Petitpierre E, Rufas JS, Massa B, Kamel Ben Halima M, JPM C. B1 was the ancestor B chromosome variant in the western Mediterranean area in the grasshopper Eyprepocnemis plorans. Cytogenet Genome Res. 2013;142:54–8.
Navrátilová A, Koblížková A, Macas J. Survey of extrachromosomal circular DNA derived from plant satellite repeats. BMC Plant Biol. 2008;8:90.
Menezes-de-Carvalho NZ, Palacios-Gimenez OM, Milani D, Cabral-de-Mello DC. High similarity of U2 snDNA sequence between a and B chromosomes in the grasshopper Abracris flavolineata. Mol Gen Genomics. 2015;290:1787–92.
Teruel M, Ruíz-Ruano FJ, Marchal JA, Sánchez A, Cabrero J, Camacho JPM, Perfectti F. Disparate molecular evolution of two types of repetitive DNAs in the genome of the grasshopper Eyprepocnemis plorans. Heredity. 2014;112:531–42.
Silva DMZA, Pansonato-Alves JC, Utsunomia R, Araya-Jaime C, Ruiz-Ruano FJ, Daniel SN, Hashimoto DT, Oliveira C, Camacho JPM, Porto-Foresti F, Foresti F. Delimiting the origin of a B chromosome by FISH mapping, chromosome painting and DNA sequence analysis in Astyanax paranae (Teleostei, Characiformes). PLoS One. 2014;9:e94896.
Ugarković D. Functional elements residing within satellite DNAs. EMBO Rep. 2005;6:1035–9.
Ferreira D, Meles S, Escudeiro A, Mendes-da-Silva A, Adega F, Chaves R. Satellite non-coding RNAs: the emerging players in cells, cellular pathways and cancer. Chromosom Res. 2015;23:479–93.
Pezer Z, Brajković J, Feliciello I, Ugarković D. Transcription of satellite DNAs in insects. Prog Mol Subcell Biol. 2011;51:161–78.
Věchtová P, Dalíková M, Sýkorová M, Žurovcová M, Füssy Z, Zrzavá M. CpSAT-1, a transcribed satellite sequence from the codling moth, Cydia pomonella. Genetica. 2016;144(4):385–95.
Dalíková D, Zrzavá M, Kubíčková S, Marec F. W-enriched satellite sequence in the Indian meal moth, Plodia interpunctella (Lepidoptera, Pyralidae). Chromosom Res. 2017; in press
Kuhn GC, Satellite DNA. Transcripts have diverse biological roles in Drosophila. Heredity (Edinb). 2015;115:1–2.
Rouleux-Bonnin F, Bigot S, Bigot Y. Structural and transcriptional features of Bombus terrestris satellite DNA and their potential involvement in the differentiation process. Genome. 2004;47:877–88.
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, Russell DW. Molecular cloning, a laboratory manual. 3rd edition. New York: Cold Spring Harbor Laboratory Press; 2001.
Drummond AJ, Ashton B, Cheung M, Heled J, Kearse M, Moir R, Stones-Havas S, Thierer T, Wilson A. Geneious v4.8.5. 2009; Available from. http://www.geneious.com.
Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;25:1451–2.
Novák P, Neumann P, Pech J, Steinhaisl J, Macas J. RepeatExplorer: a galaxy-based web server for genome-wide characterization of eukaryotic repetitive elements from next generation sequence reads. Bioinformatics. 2013;29:792–3.
Novák P, Neumann P, Macas J. Graph-based clustering and characterization of repetitive sequences in next-generation sequencing data. BMC Bioinformatis. 2010;11:378.
Pinkel D, Lanlegent J, Collins C, Fuscoe J, Segraves R, Lucas J, Gray J. Fluorescence in situ hybridization with human chromosome-specific libraries: detection of trisomy 21 and translocations of chromosome 4. Proc Natl Acad Sci U S A. 1988;85:9138–42.
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.
Barros AV, Sczepanski TS, Cabrero J, Camacho JPM, Vicari MR, Artoni RF. Fiber FISH reveals diferente patterns of high-resolution physical mapping for repetitive DNA in fish. Aquaculture. 2011;322-323:47–50.
Anjos A, Loreto V, Cabral-de-Mello DC. Organization of some repetitive DNAs and B chromosomes in the grasshopper Eumastusia koebelei koebelei (Rehn, 1909) (Orthoptera, Acrididae, Leptysminae). Comp Cytogenet. 2016;10:219–28.
Nguyen P, Sykorova M, Sichova J, Kůta V, Dalíková M, Čapková Frydrychová R, Neven LG, Sahara K, Marec F. Neo-sex chromosomes and adaptive potential in tortricid pests. Proc Natl Acad Sci U S A. 2013;110:6931–6.
This study was partially supported by the Fundação de Amparo a Pesquisa do Estado de São Paulo-FAPESP (process number 2011/19481-3, 2015/16661-1) and Coordenadoria de Aperfeiçoamento de Pessoal de Nível Superior-CAPES. DM was supported by FAPESP (process number 2015/05246-3). The authors are grateful to Dr. Thiago Gazoni for sampling in Manaus/AM, Paranaíta/MT and to the administration of Parque Estadual Edmundo Navarro de Andrade for collecting authorization. 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).
Ethics approval and consent to participate
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
qPCR data of gDNA used to calculate the sequence dose by Ct method of relative quantification (Nguyen et al. 2013). Gene dosage ratios (GDR) of the target genes were compared with autosomal gene Hsp70. (XLSX 90 kb)
(a) Dot blot analysis, (b) PCR electrophoresis of AflaSAT-1 and (c) Nucleotide alignment of positively amplified AflaSAT-1 fragments. Af = Abracris flavolineata, 1 = A. dilecta, 2 = Vilerna rugulosa, 3 = Ronderosia bergi, 4 = Schistocerca pallens, 5 = Amblytropidia robusta, 6 = Eumastusia koebelei koebelei. (JPEG 363 kb)