Abstract
Tandem repeats are important parts of eukaryotic genomes being crucial e.g., for centromere and telomere function and chromatin modulation. In Lepidoptera, knowledge of tandem repeats is very limited despite the growing number of sequenced genomes. Here we introduce seven new satellite DNAs (satDNAs), which more than doubles the number of currently known lepidopteran satDNAs. The satDNAs were identified in genomes of three species of Crambidae moths, namely Ostrinia nubilalis, Cydalima perspectalis, and Diatraea postlineella, using graph-based computational pipeline RepeatExplorer. These repeats varied in their abundance and showed high variability within and between species, although some degree of conservation was noted. The satDNAs showed a scattered distribution, often on both autosomes and sex chromosomes, with the exception of both satellites in D. postlineella, in which the satDNAs were located at a single autosomal locus. Three satDNAs were abundant on the W chromosomes of O. nubilalis and C. perspectalis, thus contributing to their differentiation from the Z chromosomes. To provide background for the in situ localization of the satDNAs, we performed a detailed cytogenetic analysis of the karyotypes of all three species. This comparative analysis revealed differences in chromosome number, number and location of rDNA clusters, and molecular differentiation of sex chromosomes.
Introduction
Sex chromosomes are extremely dynamic components of the genomes with high inter- and intraspecific variability. They evolved de novo multiple times from an ordinary pair of autosomes across various Eukaryote taxa (; ; ; ) via acquisition of a master sex-determining locus in one of the two homologs of an autosomal pair (; ; ; ). This event is followed by suppression of recombination, which leads to accumulation of repetitive DNAs and gene pseudogenization. Ultimately, this gradual genetic erosion may result in the disappearance of Y or W chromosomes (; ).
Accumulation of repeats leads to heterochromatinization of Y or W in many species (; ). Repetitive DNAs that accumulate on sex chromosomes often include satellite DNAs (satDNAs). These non-coding sequences occur in the genomes in hundreds to thousands of copies arranged in tandem in a head-to-tail manner. In general, satDNAs are highly dynamic and their nucleotide sequence, copy number, length of monomers, and chromosomal location can change quickly (). Empirical data obtained from various animal taxa, such as orthopterans (; ), mammals (; ), fishes (; ), reptiles (), and anurans (), have demonstrated the role of satDNAs in sex chromosome differentiation.
Lepidoptera (moths and butterflies) with about 160,000 described species () is the largest animal group with female heterogamety (). Most lepidopterans harbor a simple WZ sex chromosome system in females, although Z0 or multiple sex chromosomes have been reported in some species (; ; ). Extensive efforts to understand the sex chromosome composition and evolution in Lepidoptera have been made using variable approaches, from classical cytogenetics to Z-linked gene mapping, comparative genomic hybridization (CGH) of male and female genomic DNAs (gDNAs), genomic in situ hybridization (GISH) of female gDNA, W-chromosome painting, and sequencing of male and female genomes (; , ; , ; ; ; ; ; Zrzavá et al., 2018). These studies documented that the Z is a gene-rich, autosome-like chromosome with conserved synteny (; ). In contrast, the W chromosome is largely heterochromatic, composed of repetitive sequences and mostly lacking protein-coding genes (; ). The heterochromatin nature of the W chromosomes leads to the formation of sex chromatin, a roundish body in interphase nuclei of females (; ).
Due to the prevalence of repeats, attempts to sequence the lepidopteran W chromosome are scarce (; ; ). These efforts, however, have shown that the W consists mainly of mobile elements. Other studies that focused directly on satDNA detected a few satellites located on the W chromosome in several species, namely in Plodia interpunctella (), Mamestra brassicae (), and Spodoptera frugiperda (). However, low amounts of satDNA appear to be a general feature of lepidopteran genomes, which are generally rather small, ranging from 0.29 pg in Danaus plexippus (Danaidae) to 1.94 pg in Euchlaena irraria (Geometridae) (, Animal Genome Size Database). So far, only a total of five satDNAs have been identified in all lepidopteran species investigated (; ; ; ; ). They showed variable patterns of chromosomal distribution including W specific satDNA (), satDNA shared exclusively by Z and W chromosomes (), and satDNAs spread on multiple chromosomes (; ).
The small sizes of lepidopteran genomes and their repetitive composition, mainly consisting of transposable elements (TEs), may have hindered the isolation of satDNAs using classical methods, i.e., restriction endonuclease digestion of gDNA, which depends on high abundance of a particular repeat (). Only recently, this difficulty has been overcome by bioinformatic analysis of low-coverage sequenced genomes using RepeatExplorer (). This approach allowed the characterization of multiple satDNAs and generated valuable chromosomal and genomic information, for example in crickets (), Drosophila species (), fishes (; ), amphibians (), and some plants (; Zwyrtková et al., 2020). However, this approach has not yet been applied to Lepidoptera.
The lepidopteran family Crambidae has more than 10,000 species described worldwide (), some of which are serious pest of agricultural crops such as sugarcane, maize, rice, and sorghum, and are economically important (; ; ). Despite the overall significance of crambid moths, very little is known about their genome architecture, because cytogenetics of the Crambidae was poorly explored. In a few crambid species, only chromosome numbers are known (), including the reduced n = 17 in the sugarcane borer Diatraea saccharalis () and the ancestral n = 31 in Ostrinia representatives (; ; ). Sex chromosomes have only been studied in O. scapulalis () and O. nubilalis (). In O. nubilalis, a more detailed chromosomal analysis was performed with the assignment of the 31 chromosomes by gene-based fluorescence in situ hybridization (FISH) mapping (). To expand our knowledge of this group, we analyzed karyotypes and genomes of three representatives of Crambidae, namely the European corn borer (Ostrinia nubilalis), one of the most important pests causing economical losses to corn growers (), the box tree moth (Cydalima perspectalis), an Asian species that has been introduced in Europe, causing defoliation of the ornamental shrub Buxus spp. (), and the Guatemalan sugarcane borer (Diatraea postlineella), a pest which causes significant damage to sugarcane in Guatemala (; ). Our aim was to investigate the role of tandem repeats in general architecture of genomes and in sex chromosome differentiation in studied species, as little is known about these sequences in Lepidoptera, which are otherwise important components of eukaryotic genomes. To accomplish this, we used for the first time in Lepidoptera clustering analysis performed by RepeatExplorer on male and female genomic sequences of all three species to identify and map satDNAs. Further, we characterized the chromosomes of the studied species in order to provide a cytogenetic background not only to our FISH experiments but also to existing or upcoming sequencing projects. In particular, we focused on sex chromosomes, knowledge of which is important for potential pest control, such as the sterile insect technique, as O. nubilalis is one of the most important pests causing economical losses to corn growers ().
Materials and Methods
Animals
We studied three species of the family Crambidae, the box tree moth (C. perspectalis), the Guatemalan sugarcane borer (D. postlineella), and the European corn borer (O. nubilalis). A small laboratory colony of C. perspectalis was established from the last instar larvae collected from infested shrubs of Buxus spp. in Valtice (Czechia) in May 2018. The colony was kept for several generations on fresh leaves of Buxus spp. Eggs of D. postlineella were obtained from a mass-reared colony at the Santa Ana sugarcane farm, Escuintla, Guatemala City, Guatemala. Larvae of this species were reared on artificial diet prepared according to the supplier’s recipe. O. nubilalis (Z and E strains) was obtained from a laboratory colony kept at the Department of Entomology, Max Planck Institute of Chemical Ecology, Jena, Germany. For genomic analysis, we used sequenced genomes from both O. nubilalis strains, but only Z-strain was used for further chromosomal studies. Larvae were reared on an artificial wheat germ-based diet () with some modifications, and adults fed a 10% solution of honey in water. Cultures of all three species were maintained at 20–22°C and 12-h light/12-h dark regime.
Polyploid Nuclei Preparation and Microdissection of the W Chromatin
Polyploid interphase nuclei were prepared from Malpighian tubules of male and female larvae of the last instar (see ). Malpighian tubules were dissected in physiological solution, fixed in Carnoy’s fixative (ethanol, chloroform, acetic acid, 6:3:1) on a slide for about 1 min, stained with 1.25% lactic acetic orcein for 3–5 min, and mounted in the staining solution. The slides were then examined under a light microscope for the presence of W chromatin in females and its absence in males (reviewed in ).
For laser microdissection of W chromatin bodies, we followed the procedure described in with some modifications. Malpighian tubules were dissected from last instar female larvae, swollen for 10 min in a hypotonic solution (75 mM KCl) and fixed in Carnoy’s fixative for 15 min. The tubules were then transferred into a drop of 60% acetic acid on a glass slide coated with a polyethylene naphthalate membrane (Goodfellow, Huntingdon, United Kingdom) and torn to pieces using tungsten needles. The cell suspension was spread at 40°C using a heating plate and stained with 5% Giemsa for 7 min. Microdissection of W chromatin bodies was performed using the PALM MicroLaser System (Carl Zeiss MicroImaging GmbH, Munich, Germany) as described in .
Chromosome Preparations
Spread chromosome preparations were made as described previously (; ). Mitotic chromosomes were obtained from wing imaginal disks of the last instar larvae of both sexes and meiotic chromosomes in the pachytene stage from testes of the penultimate and last instar larvae and ovaries of the last instar larvae and early pupae. Wing imaginal disks and testes were dissected in a physiological solution, hypotonized for 10 min in 75 mM KCl and fixed in Carnoy’s solution for at least 15 min, whereas the ovaries were transferred to Carnoy’s solution for 15 min immediately after dissection. The fixed tissues were macerated on a slide in a drop of 60% acetic acid and then spread at 45°C using a heating plate. The slides were inspected under a phase contrast microscope and preparations of sufficient quality were passed through an ethanol series (70, 80, and 100%, 30 s each) and stored at –20°C until use.
Genome Sequencing, Satellite DNA Identification and Analysis
Genomic DNA was extracted from one male larva and one female larva of each species using the CTAB (hexadecyltrimethylammonium bromide) DNA isolation procedure according to with the following modifications. The concentration of EDTA in the extraction buffer was doubled, i.e., 40 mM, the homogenized material was incubated overnight, the centrifugation steps were performed at 14,000 g and prior to isopropanol precipitation, the samples were treated with 62.25 μg/mL RNase A for 30 min at 37°C.
Paired-end sequencing (2 × 150 bp) using the Illumina HiSeq 4000 system was performed by Novogene (HK) Co., Ltd., (Hong Kong, China). We checked the quality of reads with FastQC () and processed the reads though Toolkit suit using default options () to be used as input for RepeatExplorer (available at1) clustering. We made separate identification of satDNAs for each three species. A comparative analysis of satDNAs between the sexes for D. postlineella and C. perspectalis was performed using 250,000 reads of each sex as input. For O. nubilalis, genomes of both sexes from the two strains (E and Z) were used, being 125,000 reads per genome. After clustering by RepeatExplorer (, ) using default options, the satDNAs identified as satDNAs by TAREAN tool () were considered for subsequent analysis. The tandem arrangement was checked by dotplot and Tandem Repeats Finder, TRF ().
For similarity analysis, we performed all-against-all comparison of monomers of the recovered satDNAs using RepeatMasker () and the “rm_homology.py” script2. We used RepeatMasker to calculate divergence and abundance of each satDNA in the female and the male genomes at intraspecific level. For this purpose, we randomly selected 7.5 million read pairs per library obtained by seqtk tool3 and aligned them against dimers of consensus satDNA sequences. To estimate the average Kimura 2-parameter distances (K2P) for each satDNA family, we used the calcDivergenceFromAlign.pl script from the RepeatMasker utility tool. Genomic abundance for each satDNA family was calculated according to the proportion of nucleotides aligned with the reference consensus sequence divided by the library size. Finally, we compared the divergence of satDNAs between sexes generating repeat landscapes showing the relative abundance of repeat elements on the Y-axis and 1% intervals of K2P distance from the consensus on the X-axis. In addition, we also performed the RepeatMasker analysis at interspecific level to check the occurrence of different satDNAs in the species studied and to check their abundance. Similarity of the satDNAs identified with previously characterized sequences was checked by blast search in NCBI and Repbase.
DNA Probes and Fluorescence in situ Hybridization
For each studied species, we prepared specific gDNA probes, W-chromosome painting probes, and satDNA probes. We also prepared non-species specific probes to localize conserved sequences including telomeric repeats and 18S rDNA. These probes were used for different FISH experiments following different protocols, as specified below.
Female and male gDNAs extracted by the CTAB procedure were labeled by nick translation with Cy3-dUTP for female gDNA or fluorescein-12-dUTP (both Jena Bioscience, Jena, Germany) for male gDNA, or vice versa. The nick translation mixture was composed of 500 ng of gDNA, 25 μM of each dATP, dCTP, and dGTP, 9 μM of dTTP, 16 μM of labeled nucleotides, nick translation buffer (50 mM Tris–HCl, pH 7.5, 5 mM MgCl2, 0.005% BSA), 10 mM β-mercaptoethanol, 20 U DNA Polymerase I (ThermoFisher Scientific, Waltham, MA, United States), and 0.005 U DNase I (ThermoFisher Scientific). The mixture was incubated for 2 h at 15°C. These labeled gDNA probes were used for CGH following the protocol described in with modifications proposed by .
W-chromosome painting probes were prepared as described in . Briefly, samples of microdissected W chromatin were amplified by PCR using a GenomePlex Single Cell Whole Genome Amplification Kit, WGA4 (Sigma-Aldrich, St. Louis, MO, United States). The amplified product was re-amplified and labeled with either Cy3-dUTP or fluorescein-12-dUTP by PCR using a GenomePlex WGA Reamplification Kit, WGA3 (Sigma-Aldrich). These probes were used in W-chromosome painting experiments at intraspecific level and for cross-species W-painting following the CGH protocol (). For intraspecific experiments, we combined the W-painting probe with a (TTAGG)n telomeric probe obtained by non-template PCR according to and using primers (TTAGG)5 and (CCTAA)5. The telomeric probe was labeled with Cy3-dUTP by using the improved nick translation procedure of with some modifications (see ).
18S rDNA and satDNA probes were amplified by PCR from female gDNA and labeled by PCR or nick translation with biotin-16-dUTP (Roche Diagnostics, Mannheim, Germany). The 18S rDNA was obtained from codling moth (Cydia pomonella) gDNA using the primers described by . This probe was labeled by nick translation. For satDNA amplification we designed primers using primer3 () or manually (Supplementary Table 1). PCR was carried out in 20-μl reaction volume containing 1.5 μM of each primer, 1 × Ex Taq buffer (TaKaRa), 0.5 U Ex Taq polymerase (TaKaRa), 200 μM of each nucleotide and either 50–100 ng gDNA. The thermal cycling profile consisted of initial denaturation at 95°C for 5 min, 30 cycles of denaturation at 94°C for 30 s, annealing at variable temperature (see Supplementary Table 1) for 30 s, elongation at 72°C for 1 min and final elongation 72°C for 3 min. After PCR, the fragments of satDNAs were visualized on a 1.5% agarose gel and monomers were extracted from the gel and re-amplified by PCR. The re-amplified products were sequenced by the Sanger method using the service of SEQme (Dobříš, Czechia) to verify their reliability. For satDNAs probe labeling, a second PCR was done with incorporation of biotin-16-dUTP. For FISH with these probes, we followed the protocol described in . The probes were detected by Cy3-conjugated streptavidin (Jackson Immuno Res. Labs. Inc., West Grove, PA, United States).
In experiments in which the W-chromosome painting probe and satDNAs were mapped on the same slide, we applied two rounds of FISH. We first performed W-chromosome painting using the CGH protocol (), then FISH with satDNA probe (). In all FISH experiments, chromosomes were counterstained with DAPI (4′,6-diamidino-2-phenylindole, dihydrochloride; Sigma-Aldrich) and slides mounted with antifade based on DABCO (1,4-diazabicyclo(2.2.2)-octane; Sigma-Aldrich).
Results
Karyotypes and Chromosomal Location of 18S rDNA
Diploid chromosome numbers were determined by analysis of mitotic metaphase from wing imaginal disks of females and males stained with DAPI. In C. perspectalis and O. nubilalis we observed 2n = 62 chromosomes [Supplementary Figure 1A; cf. for chromosome number in O. nubilalis]. However, in D. postlineella we found a reduced number of chromosomes to 2n = 42 (Supplementary Figure 1B). We also counted bivalents in pachytene of females or males, confirming the number of chromosomes in all three species (Supplementary Figures 2A–C). In D. postlineella, meiotic bivalents often formed clumps at the pachytene stage (Supplementary Figure 2B), probably due to their longer length compared to the other two species. None of the studied species showed large heterochromatin blocks.
A single sex chromatin body, deeply stained with orcein and of regular spherical shape, was observed in each polyploid nucleus of Malpighian tubules from females, indicating the presence of the W chromosome in all three species (Supplementary Figures 3A,C,E). As expected, sex chromatin was absent in the polyploid nuclei of Malphigian tubules from males (Supplementary Figures 3B,D,E). This finding suggests that the three crambid species have a WZ/ZZ sex chromosome system, which was confirmed by further analysis.
Fluorescence in situ hybridization with 18S rDNA probe applied to male pachytene nuclei revealed significant differences between the species studied in the number and location of rDNA clusters (Figures 1a–c). In C. perspectalis, two clusters were found, each located in the terminal position of a different bivalent (Figure 1a), whereas in D. postlineella only one interstitial cluster was observed (Figure 1b). In contrast, four terminal clusters in four different bivalents were found in O. nubilalis (Figure 1c).
FIGURE 1
Differentiation of Sex Chromosomes by CGH and W-Painting Probes
To investigate the gross molecular differentiation of W and Z chromosomes, we performed CGH. In all three species, the W chromosome showed stronger binding of probes derived from both female and male gDNAs compared to the other chromosomes. In addition, all species showed more intense labeling of female-derived probes in comparison with male-derived probes on the W chromosome, indicating that this chromosome consists mainly of female-specific or female-enriched sequences (Figures 2a–l).
FIGURE 2
The W-painting probe specifically highlighted the entire W chromosome at the intraspecies level, allowing easy identification of the sex chromosome bivalent and confirming considerable molecular differentiation of the W chromosome from other chromosomes in all three species. The telomeric probe hybridized only to the end of the chromosomes, which also confirmed that the sex chromosome bivalent is composed of only two elements (WZ) (Figures 3a,d,g). Using cross-species W-painting, we investigated the molecular differentiation of W chromosomes between the three species. Cross-hybridization of W-painting probes resulted in scattered signals on the W chromosome of another species in all cases, but in varying numbers and intensities, depending on the species (Figures 3b,c,e,f,h,i). In general, after FISH with a W-probe from another species, W chromosomes were decorated with several to multiple clusters of hybridization signals along the entire length. The highest intensity of labeling was recorded in the W chromosome of D. postlineella using the W-painting probe of O. nubilalis (Figure 3e). In contrast, the lowest signal intensity in cross-species W-painting experiments was observed using the W-probe from C. perspectalis in the W chromosome of O. nubilalis (Figure 3h). Only weak hybridizations signals were observed on some autosomes (results not shown).
FIGURE 3
satDNAs Genomic Characterization and Chromosomal Mapping
By clustering analysis, separately for each species, in RepeatExplorer using the TAREAN report, we identified a total of seven putative satDNAs in the three species studied, some of them with high confidence and others with low confidence (a classification provided by TAREAN) (Supplementary Figure 4). One satDNA was identified in C. perspectalis (Cper-Sat01), two in D. postlineella (Dpos-Sat01 and Dpos-Sat02), and four in O. nubilalis (Onub-Sat01, Onub-Sat02, Onub-Sat03, and Onub-Sat04). Subsequent analyses were performed for these seven putative satDNAs, demonstrating their tandem arrangement by dotplot and TRF and also by PCR, which showed a typical ladder-like pattern. Monomers of these satDNAs were highly variable in size, ranging from 123 to 2,244 bp. These satDNAs were A + T-rich with G + C content ranging from 36.44 to 49.11%. They represented only a small portion of the male and female genomes in each species, ranging from 0.01909 to 0.15464%, showing similar abundance between the sexes. The divergences for each satDNA family ranged from 1.32 to 11.34% and were similar between the sexes. These data are summarized in Table 1 and the sequence logos are shown in Supplementary Figure 4. SatDNA landscapes, which were generated to compare the repeats between sexes, aiming to check the possible amplification of variants with distinct degrees of divergence, did not reveal any sex-specific differences (Figure 4). We also searched for all satDNA families at the interspecific level by RepeatMasker and found that they are present in all three genomes, but with a lower proportion compared to the species in which they were identified (Supplementary Table 2). Finally, our Repbase and NCBI searching did not reveal relevant similarity with any described sequences. The sequences were deposited in GenBank under the accession numbers MW369067–MW369073.
TABLE 1
| Species satDNA family | Monomer size (bp) | G + C content (%) | Abundance (%) | Divergence (%) | ||
| Female | Male | Female | Male | |||
| Cydalima perspectalis | ||||||
| Cper-Sat01 | 2,244 | 43.36 | 0.14476 | 0.12624 | 7.39 | 7.14 |
| Diatraea postlineella | ||||||
| Dpos-Sat01 | 1,391 | 36.44 | 0.04980 | 0.04895 | 10.1 | 11.23 |
| Dpos-Sat02 | 1,380 | 40.57 | 0.01929 | 0.02431 | 1.32 | 1.37 |
| Ostrinia nubilalis | ||||||
| Onub-Sat01 | 528 | 39.77 | 0.15464 | 0.11740 | 9.38 | 9.55 |
| Onub-Sat02 | 453 | 40.39 | 0.05552 | 0.07105 | 11.34 | 11.31 |
| Onub-Sat03 | 957 | 49.11 | 0.02807 | 0.02993 | 2.05 | 2.08 |
| Onub-Sat04 | 123 | 37.39 | 0.02839 | 0.02604 | 8.86 | 9.28 |
Main characteristics of SatDNA families identified in Crambidae genomes.
FIGURE 4
Fluorescence in situ hybridization mapping of satDNAs revealed a variable distribution pattern with clusters on autosomes and sex chromosomes (Figures 5–7). Hybridization signals of Cper-Sat01 were scattered in all chromosomes of C. perspectalis, with no apparent enrichment in specific chromosomal regions (Figure 5a). However, analysis of interphase nuclei revealed enrichment of this satDNA on W chromatin (Figure 5b). At the pachytene stage, in combined analysis with the W-probe, clusters of Cper-Sat01 were enriched and distributed along the entire length of the W chromosome, but were almost absent on the Z chromosome (Figure 5c). In D. postlineella, two clusters of Dpos-Sat01 were located interstitially in a pair of autosomes carrying the nucleolar organizer region (NOR) formed by the major rDNA (Figure 6a). FISH with Dpos-Sat02 probe identified a single cluster that was located on another autosomal pair in the interstitial position (Figure 6b).
FIGURE 5
FIGURE 6
FIGURE 7
Fluorescence in situ hybridization mapping of four satDNAs identified in the genome of O. nubilalis showed scattered hybridization signals on most chromosomes. However, these satDNAs were distributed differently in various chromosomal regions. Onub-Sat02 and Onub-Sat03 were virtually scattered and formed clusters of small dots along the entire length of most chromosomes, while Onub-Sat01 and Onub-Sat04 were enriched primarily in the terminal regions of the chromosomes (Figures 7a,b). Concerning the WZ sex chromosomes, Onub-Sat01 (Figure 7c), Onub-Sat02 (Figure 7d), and Onub-Sat03 (Figure 7e) showed a similar distribution in both sex chromosomes. Remarkably, large terminal clusters occurred in Onub-Sat01, in addition to several smaller clusters located mainly on the W chromosome (Figure 7c). In contrast, Onub-Sat04 was located exclusively at one termini of the Z chromosome (Figure 7f).
Discussion
Karyotypes and Molecular Differentiation of Sex Chromosomes in Crambidae
In this study, we made progress in understanding the basic features of karyotypes and molecular differentiation of sex chromosomes in the economically important family of Lepidoptera, Crambidae. Although the number of Crambidae species karyotyped so far is low, a high variability of haploid numbers, ranging from n = 10 to n = 41 has been reported (; ; ; ; ). These variations in the number of chromosomes below and above the ancestral number of n = 31 (; ) suggest a dynamic karyotype evolution by chromosome fusions and fissions, respectively. Fusions probably also played a major role in reducing the number of chromosomes to n = 21 found in this study in D. postlineella. The occurrence of large chromosomes in D. postlineella, which were difficult to individualize in pachytene compared to O. nubilalis and C. perspectalis with n = 31, also supports the fusion hypothesis. The fact that a related species, the sugarcane borer D. saccharalis, also has a reduced chromosome number (n = 17; ), suggests that this is an evolutionary trend in the genus Diatraea. In lepidopteran species, a reduction in the chromosome number by fusion is generally more common compared to an increase by fission, which could be attributed to more deleterious effects of fissions in meiotic products (; ).
The karyotypes of crambid species differed greatly in the number and chromosomal position of rDNA clusters. A single rDNA cluster in the interstitial position, as found in D. postlineella, is a common pattern for Lepidoptera (). Two terminal rDNA clusters found in the karyotype of C. perspectalis were also observed in the Mediterranean flour moth, Ephestia kuehniella, from the sister family Pyralidae (). Interestingly, the Z-strain of O. nubilalis showed four terminal rDNA clusters (this study), while specimens from another population showed five terminal clusters per haploid genome (see ), suggesting intra-specific variability in the number of rDNA clusters. Besides Crambidae, rDNA variability between species has also been reported in the lepidopteran family Tortricidae, but was mainly manifested by different chromosomal positions of rDNA clusters (). The data obtained here for Crambidae add to the growing evidence of high rDNA dynamics in Lepidoptera, which is mainly attributed to ectopic recombination not only in Lepidoptera but also in other insects (; ; ).
We have shown that the studied crambid species have a WZ/ZZ (female/male) sex chromosome system, which is common in Lepidoptera (). In all three species, the W chromosome is highly molecularly differentiated from the Z chromosome and consists mainly of female-specific and/or female-enriched repetitive sequences, as in some other Lepidoptera (; ,; Zrzavá et al., 2018). Despite this apparent similarity in gross molecular composition, cross-hybridization experiments revealed that W chromosomes exhibit distinct molecular divergence between species. Nevertheless, some degree of homology has been demonstrated, suggesting partial conservation of repeats on the W chromosomes. However, this degree of homology is not consistent with the phylogeny of Crambidae subfamilies (; Zhu et al., 2018; ), since the W probe from O. nubilalis (Pyraustinae) showed the highest number of signals in the W chromosome of distant D. postlineella (Crambinae), but a lower number of signals in the W chromosome of more closely related C. perspectalis (Spilomelinae). These results support the hypothesis of fast independent molecular divergence of W chromosomes with occasional conservation of some repetitive sequences. Partial homology between W chromosomes was also observed in species of the sister family Pyralidae, but the extent of homology was consistent with phylogenetical relationships (). In contrast, a high molecular divergence was found between the W chromosomes of two congeneric moths of the genus Abraxas, Geometridae (Zrzavá et al., 2018).
satDNAs in Crambidae
The interspecies occurrence of the seven satDNAs identified herein is consistent with the library hypothesis (). A common feature of these newly identified satDNAs is the enrichment in A + T base pairs (light satDNAs), similar to the other five satDNAs previously described in Lepidoptera (; ; ; ; ). Interestingly, the monomer length of some of the satDNAs identified herein is high, reaching 2,244 bp for Cper-Sat01 in C. perspectalis and more than 1,000 bp for two satDNAs in D. postlineella. This is an unusual feature in insects, which generally tend to have satDNA monomer sizes less than 600 bp (). Exceptions have been found, for example, in the ant Monomorium subopacum with a monomer size of 2.5 kb (), the kissing bug Triatoma infestans with about 1 kb repeat unit (), and in the beetles Misolampus goutodii and Hippodamia variegata with 1.2 and 2 kb repeat units, respectively (; ). It is worth noting that the use of different methods for satDNA prospection could affect the capacity of discovering satDNA families with greater or smaller monomer lengths.
The absence of large heterochromatic blocks in most studied lepidopterans suggests the low abundance of tandem repeats forming long arrays in the chromosomes. This assumption is supported by the extremely low abundance of satDNAs in the three crambid species, in which satDNAs reached a maximum of 0.255% of the genome (mean value between the male and female genomes), as found here in O. nubilalis. Such a small proportion of satDNAs is in line with the generally small amount of heterochromatin in Lepidoptera, which is mostly observed only in the W chromosome () or associated with NOR regions (; ). Only in exceptional cases were numerous heterochromatin blocks observed in other chromosomes (, ; Zrzavá et al., 2018). This contrasts with other insects, such as hemipterans and grasshoppers, in which a high abundance of satDNAs was revealed by bioinformatics tools (; ). This may be due to the large genome size of the latter insect groups compared to the generally small size of lepidopteran genomes (, Animal Genome Size Database) and also due to predominance of interspersed repeats in Lepidoptera ().
Five of the seven satDNAs identified in Crambidae show a scattered distribution on chromosomes. This appears to be a recurrent pattern in Lepidoptera, as it was also observed in two other satDNAs mapped on chromosomes of moth species from distant families, Saturniidae and Tortricidae (; ). The remaining two mapped satDNAs, MBSAT1 in M. brassicae (Noctuidae) and PiSAT1 in P. interpunctella (Pyralidae), were limited to sex chromosomes only (; ). The predominant scattered distribution is in contrast to patterns observed in organism with large heterochromatin-rich genomes (e.g., ; ). This might be caused by the fact that the formation of long arrays of tandem repeats requires a special environment, which is usually lacking in Lepidoptera. recently studied satDNAs of the legume plant Lathyrus sativus using ultra-long nanopore reads. They discovered that besides the long arrays of tandem repeats located in centromeric and subtelomeric heterochromatin blocks, satDNAs also formed short arrays in other parts of the genome associated with and probably originated from LTR transposons. The authors speculated that the formation of a long array of satDNA requires restriction of recombination, a process which might reduce the length of the array via homologous recombination. This condition is fulfilled, for example, in centromeres and sex chromosomes (; ). In holokinetic chromosomes of Lepidoptera, which lack a primary constriction (i.e., centromere), there are usually no suitable regions for long arrays of tandem repeats other than the W chromosome, which might be the reason why satDNAs often show a scattered distribution of small clusters. The occurrence of satDNAs in euchromatin, as found in Lepidoptera, has been reported in some insects (; ; ). Such satDNAs could be located near genes and modulate their expression (; ). This aspect deserves more research in Lepidoptera, due to the frequent location of satDNAs in euchromatin and transcription of some satDNAs (; ).
Distribution of satDNAs on Sex Chromosomes
Although it is well documented that the lepidopteran W chromosome is mostly composed of heterochromatin rich in repetitive DNA sequences, especially TEs, its specific molecular composition is largely unknown (; , ; ). Our study contributed to this knowledge and demonstrated the occurrence of satDNAs on the W chromosomes of C. perspectalis and O. nubilalis, but not in D. postlineella. It should be noted that, unlike Lepidoptera, satDNAs are widely reported in the Y chromosomes of various animals (; ; ; ) and plants (, ). Among the tandem sequences identified on the lepidopteran W chromosomes are rDNA in few species (; ; ; Zrzavá et al., 2018) and two satDNAs, PiSAT1 in P. interpunctella () and MBSAT1 in a M. brassicae cell line (). However, these satDNAs form conspicuous blocks on the W chromosomes, differing from the satDNAs of C. perspectalis and O. nubilalis, which are arranged as scattered clusters on the W chromosomes, similar to autosomes. Interestingly, AT-1 of C. pomonella is underrepresented on the W chromosome () and two satDNAs identified in D. postlineella are located exclusively on autosomes (this study). This demonstrates the great plasticity of satDNA arrangement in lepidopteran genomes, from absence to abundance on the W chromosomes and from a single cluster to a high number of scattered clusters.
The results of our bioinformatic analysis of satDNAs in both sexes of the crambid species studied suggest that in female genomes, ultimately in the W chromosome, there was no extensive amplification of these repeats, as is often reported in the Y or W chromosome in other animals (; ; ). However, our complementary chromosomal analysis revealed enrichment of some satDNAs on the W (C. perspectalis and O. nubilalis) or Z chromosome (O. nubilalis). Because satDNAs are highly dynamic in copy number (), the difference in abundance could be attributed to interindividual variability of autosomal clusters, masking the difference in abundance between sex chromosomes.
As observed for euchromatic autosomes, the absence of long arrays surprisingly also applies to the heterochromatic W chromosome in the crambid species studied. Of the seven satDNAs analyzed, four were present on the W chromosome, three of which were abundant there. The patterns of the W-signals of the respective satDNA resembled those observed on the other chromosomes, indicating that the clusters were relatively small and scattered. The only satDNA that formed slightly larger subtelomeric clusters was the Onub-Sat01 (Figure 7c). Of all the lepidopteran satDNAs known today, the only one which forms long arrays in the genome is MBSAT1 in the cabbage moth, M. brassicae, which has been shown to be vastly abundant on both the W and Z chromosomes, occupying ca 1.9% of the genome (). However, it would be worth to reanalyze this case, as these results have been obtained from a cell line that showed a drastically reduced chromosome number (n = 11) compared to the number reported for the wild population of this species (n = 31) (; ). This suggests that the cell line may have undergone genome reconstruction accompanied by amplification of this satDNA, which may not occur in wild populations. Thus, the fact that there are no long arrays of satDNA on the W chromosome that is non-recombining because lepidopteran females lack recombination on all of their chromosomes remains unexplained.
Our CGH results suggest a high abundance of W-enriched and/or W-specific sequences on the W chromosome of all three crambid species. Combining CGH with satDNA mapping helped us to better understand nature of these sequences. In D. postlineella, the W-enriched and/or W-specific sequences are not satellites because the only two satDNAs we detected are low-copy autosomal repeats. Hence, the hybridization signals of the female gDNA probe, which highlighted the W chromosome, were generated by other types of sequences, such as mobile elements or microsatellites. In the case of C. perspectalis, the only identified Cper-Sat01 is abundant on both autosomes and the W chromosome, while it is underrepresented or absent on the Z chromosome. Thus, it is enriched in females and contributes to the W-highlighting using the female gDNA probe. In O. nubilalis, satDNAs mapped to autosomes as well as both sex chromosomes, except Onub-Sat04, which was absent on the W and formed a small cluster at one Z-chromosome end. Hence, the W-enriched and/or W-specific sequences are in fact other repeats than the satDNAs studied.
Conclusion
Our study significantly contributed to the understanding of karyotype diversification, genome architecture, and sex chromosome evolution in Lepidoptera. The karyotypes of the studied species of the family Crambidae differ in some aspects, namely in the gross architecture (i.e., diploid numbers), probably due to chromosomal fusions, and in fine structures, as evidenced by variability in the distribution of repetitive DNAs. We revealed low abundance and high variability of satDNAs in these Crambidae species, which also contributed to the plasticity of sex chromosomes. The W chromosomes are highly differentiated between the three species due to independent evolution, although some degree of random conservation (not consistent with species phylogeny) of some anonymous repeats has been found. Finally, the combination of genomic and chromosomal data allowed the number of satDNAs identified in Lepidoptera to be doubled, opening new avenues for the study of this genome fraction.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author contributions
FM and DC-D-M contributed to conception and design of the study. DC-D-M and MZ performed the experiments and analyzed the data. DC-D-M, MZ, and FM interpreted the data and wrote the first draft of the manuscript. SK prepared the W chromosome probes. PR obtained samples of Diatraea postlineella and ensured the correct determination of this species. All authors contributed to manuscript revision, read and approved the submitted version.
Funding
This study was supported by grants 17-13713S and 20-13784S of the Czech Science Foundation, CNPq-Conselho Nacional de Pesquisa e Desenvolvimento, and FAPESP-Fundação de Amparo à Pesquisa do Estado de São Paulo (BE.PQ, process number 2017/097319, Regular Research Grant, process number 2019/19069-7). SK acknowledges support from grant RO 0520 of the Ministry of Education, Youth, and Sports of the Czechia. This work was conducted as part of the IAEA National project Using the Sterile Insect Technique (SIT) to Assess it as an Alternative for the Control of the Sugarcane Borer in Pilot Areas (GUA5017).
Acknowledgments
We would like to thank David G. Heckel and Melanie Unbehend for samples of Ostrinia nubilalis and Luis C. Arroyo from the Santa Ana sugarcane farm in Escuintla, Guatemala, for preparing Diatraea postlineella samples from the rearing facility. We would also like to thank Alma Solis (ARS, USDA) for access to an unpublished manuscript about D. postlineella and relevant pyraloid references. Computational resources were supplied by the project “e-Infrastruktura CZ” (e-INFRA LM2018140) provided within the program Projects of Large Research, Development and Innovations Infrastructures.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The reviewers DS and PR declared a shared affiliation with several of the authors, DC-D-M and MZ respectively, to the handling editor at the time of review.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.661417/full#supplementary-material
Supplementary Figure 1Mitotic chromosome spreads obtained from wing imaginal disks of female larvae of (A)Cydalima perspectalis (prometaphase) and (B)Diatraea postlineella (early anaphase) and stained with DAPI. Bar = 10 μm.
Supplementary Figure 2DAPI-stained pachytene bivalents obtained from testes (top images). (A)Cydalima perspectalis (n = 31), (B)Diatraea postlineella (n = 21) and (C)Ostrinia nubilalis (n = 31). A schematic drawing of the same pachytene nucleus is shown at the bottom of each panel; this way of schematic representation was used to determine the number of bivalents. Bar = 10 μm.
Supplementary Figure 3Orcein-stained polyploid nuclei obtained from Malpighian tubule cells of female (A,C,E) and male (B,D,F) larvae showing sex chromatin status in (A,B)Cydalima perspectalis, (C,D)Diatraea postlineella, and (E,F)Ostrinia nubilalis. Note the presence of a deeply stained W chromatin body (arrows) in females, while it is absent in males. Bar = 10 μm.
Supplementary Figure 4Cluster layouts and sequence logos for satDNAs. Graphical layouts were obtained from RepeatExplorer output. Sequence logos were generated from full-length consensus sequences from the most frequent k-mers for each satDNA.
Supplementary Table 1DNA sequences of primers used for PCR amplification of seven satDNAs and 18S rDNA in Crambidae species. Cper, Cydalima perspectalis; Dpos, Diatraea postlineella; Onub, Ostrinia nubilalis.
Supplementary Table 2Analysis of abundance of identified satDNAs at the interspecific level in the genomes of Crambidae species.
References
1
AbbottJ. K.NordénA. K.HanssonB. (2017). Sex chromosome evolution: historical insights and future perspectives.Proc. Biol. Sci.284:20162806. 10.1098/rspb.2016.2806
2
AbeH.MitaK.YasukochiY.OshikiT.ShimadaT. (2005). Retrotransposable elements on the W chromosome of the silkworm, Bombyx mori.Cytogenet. Genome Res.110144–151. 10.1159/000084946
3
AcostaM. J.MarchalJ. A.MartínezS.PuermaE.BullejosM.Díaz de la GuardiaR.et al (2007). Characterization of the satellite DNA Msat-160 from the species Chionomys nivalis (Rodentia, Arvicolinae).Genetica13043–51. 10.1007/s10709-006-0018-1
4
AholaV.LehtonenR.SomervuoP.SalmelaL.KoskinenP.RastasP.et al (2014). The Glanville fritillary genome retains an ancient karyotype and reveals selective chromosomal fusions in Lepidoptera.Nat. Commun.5:4737. 10.1038/ncomms5737
5
AndrewsS. (2010). FastQC: a Quality Control Tool for High Throughput Sequence Data. Available online at: http://www.bioinformatics.babrahamac.uk/projects/fastqc/. (Accessed March 18, 2019).
6
BardellaV. B.MilaniD.Cabral-de-MelloD. C. (2020). Analysis of Holhymenia histrio genome provides insight into the satDNA evolution in an insect with holocentric chromosomes.Chromosome Res.28369–380. 10.1007/s10577-020-09642-1
7
BensonG. (1999). Tandem repeats finder: a program to analyze DNA sequences.Nucleic Acids Res.27573–580. 10.1093/nar/27.2.573
8
BrajkovićJ.FelicielloI.Bruvo-MađarićB.UgarkovićÐ (2012). Satellite DNA-like elements associated with genes within euchromatin of the beetle Tribolium castaneum.G3 (Bethesda)2931–941. 10.1534/g3.112.003467
9
BullJ. J. (1983). Evolution of Sex Determining Mechanisms.San Francisco, CA: The Benjamin/Cummings Publishing Company, Inc.
10
Cabral-de-MelloD. C.MarecF. (2021). Universal fluorescence in situ hybridization (FISH) protocol for mapping repetitive DNAs in insects and other arthropods.Mol. Genet. Genom.10.1007/s00438-021-01765-2Online ahead of print.
11
Cabral-de-MelloD. C.MouraR. C.MartinsC. (2011). Cytogenetic mapping of rRNAs and histone H3 genes in 14 species of Dichotomius (Coleoptera, Scarabaeidae, Scarabaeinae) beetles.Cytogenet. Genome Res.134127–135. 10.1159/000326803
12
CamachoJ. P.Ruiz-RuanoF. J.Martín-BlázquezR.López-LeónM. D.CabreroJ.LoriteP.et al (2015). A step to the gigantic genome of the desert locust: chromosome sizes and repeated DNAs.Chromosoma124263–275. 10.1007/s00412-014-0499-0
13
CharlesworthB. (1991). The evolution of sex chromosomes.Science2511030–1033. 10.1126/science.1998119
14
CharlesworthD.CharlesworthB.MaraisG. (2005). Steps in the evolution of heteromorphic sex chromosomes.Heredity95118–128. 10.1038/sj.hdy.6800697
15
CrepaldiC.Parise-MaltempiP. P. (2020). Heteromorphic sex chromosomes and their DNA content in fish: an insight through satellite DNA accumulation in Megaleporinus elongatus.Cytogenet. Genome Res.16038–46. 10.1159/000506265
16
da SilvaM. J.Fogarin DestroR.GazoniT.NarimatsuH.Pereira, dos SantosP. S.et al (2020). Great abundance of satellite DNA in Proceratophrys (Anura, Odontophrynidae) revealed by genome sequencing.Cytogenet. Genome Res.160141–147. 10.1159/000506531
17
DalíkováM.ZrzaváM.HladováI.NguyenP.ŠonskýI.FlegrováM.et al (2017a). New insights into the evolution of the W chromosome in Lepidoptera.J. Hered.108709–719. 10.1093/jhered/esx063
18
DalíkováM.ZrzaváM.KubíčkováS.MarecF. (2017b). W-enriched satellite sequence in the Indian meal moth, Plodia interpunctella (Lepidoptera, Pyralidae).Chromosome Res.25241–252. 10.1007/s10577-017-9558-8
19
de VosJ. M.AugustijnenH.BätscherL.LucekK. (2020). Speciation through chromosomal fusion and fission in Lepidoptera.Philos. Trans. R. Soc. Lond. B Biol. Sci.37520190539. 10.1098/rstb.2019.0539
20
EscudeiroA.AdegaF.RobinsonT. J.Heslop-HarrisonJ. S.ChavesR. (2019). Conservation, divergence, and functions of centromeric satellite DNA families in the Bovidae.Gen. Biol. Evol.111152–1165. 10.1093/gbe/evz061
21
FelicielloI.AkrapI.UgarkovićÐ (2015). Satellite DNA modulates gene expression in the beetle Tribolium castaneum after heat stress.PLoS Genet.11:e1005466. 10.1371/journal.pgen.1005466
22
FerrettiA. B. S. M.MilaniD.Palacios-GimenezO. M.Ruiz-RuanoF. J.Cabral-de-MelloD. C. (2020). High dynamism for neo-sex chromosomes: satellite DNAs reveal complex evolution in a grasshopper.Heredity125124–137. 10.1038/s41437-020-0327-7
23
FerrettiA. B. S. M.Ruiz-RuanoF. J.MilaniD.LoretoV.MartíD. A.RamosE.et al (2019). How dynamic could be the 45S rDNA cistron? An intriguing variability in a grasshopper species revealed by integration of chromosomal and genomic data.Chromosoma128165–175. 10.1007/s00412-019-00706-8
24
FraïsseC.PicardM. A. L.VicosoB. (2017). The deep conservation of the Lepidoptera Z chromosome suggests a non-canonical origin of the W.Nat. Commun.8:1486. 10.1038/s41467-017-01663-5
25
FryK.SalserW. (1977). Nucleotide sequences of HS-α satellite DNA from kangaroo rat Dipodomys ordii and characterization of similar sequences in other rodents.Cell121069–1084. 10.1016/0092-8674(77)90170-2
26
FukováI.NguyenP.MarecF. (2005). Codling moth cytogenetics: karyotype, chromosomal location of rDNA, and molecular differentiation of sex chromosomes.Genome481083–1092. 10.1139/g05-063
27
FukováI.TrautW.VítkováM.KubíčkováS.MarecF. (2007). Probing the W chromosome of the codling moth, Cydia pomonella, with sequences from microdissected sex chromatin.Chromosoma116135–145. 10.1007/s00412-006-0086-0
28
FurmanB. L. S.MetzgerD. C. H.DaroltiI.WrightA. E.SandkamB. A.AlmeidaP.et al (2020). Sex chromosome evolution: so many exceptions to the rules.Genome Biol. Evol.12750–763. 10.1093/gbe/evaa081
29
Garrido-RamosM. A. (2017). Satellite DNA: an evolving topic.Genes8:230. 10.3390/genes8090230
30
GattoK. P.MattosJ. V.SegerK. R.LourençoL. B. (2018). Sex chromosome differentiation in the frog genus Pseudis involves satellite DNA and chromosome rearrangements.Front. Genet.9:301. 10.3389/fgene.2018.00301
31
GiovannottiM.CerioniP. N.RojoV.OlmoE.SlimaniT.SplendianiA.et al (2018). Characterization of a satellite DNA in the genera Lacerta and Timon (Reptilia, Lacertidae) and its role in the differentiation of the W chromosome.J. Exp. Zool. Part B Mol. Develop. Evol.33083–95. 10.1002/jez.b.22790
32
GordonA.HannonG. J. (2010). Fastx-toolkit. FASTQ/A Short-reads Pre-processing Tools. Available online at: http://hannonlab.cshl.edu/fastx_toolkit. (Accessed April 8, 2020).
33
GregoryT. R. (2020). Animal Genome Size Database. Available online at: http://www.genomesize.com. (Accessed November 26, 2020).
34
GuthrieW. D.DollingerE. J.StetsonJ. F. (1965). Chromosome studies of the European corn borer, smartweed borer, and lotus borer (Pyralidae).Ann. Entomol. Soc. Amer.58100–105. 10.1093/aesa/58.1.100
35
HejníčkováM.KouteckýP.PotockýP.ProvazníkováI.VoleníkováA.DalíkováM.et al (2019). Absence of W chromosome in Psychidae moths and implications for the theory of sex chromosome evolution in Lepidoptera.Genes10:1016. 10.3390/genes10121016
36
HobzaR.CeganR.JesionekW.KejnovskyE.VyskotB.KubatZ. (2017). Impact of repetitive elements on the Y chromosome formation in plants.Genes8:302. 10.3390/genes8110302
37
HobzaR.KubatZ.CeganR.JesionekW.VyskotB.KejnovskyE. (2015). Impact of repetitive DNA on sex chromosome evolution in plants.Chromosome Res.23561–570. 10.1007/s10577-015-9496-2
38
IjdoJ. W.WellsR. A.BaldiniA.ReedersS. T. (1991). Improved telomere detection using a telomere repeat probe (TTAGGG)n generated by PCR.Nucleic Acids Res.19:4780. 10.1093/nar/19.17.4780
39
KageymaD.TrautW. (2004). Opposite sex-specific effects of Wolbachia and interference with the sex determination of its host Ostrinia scapulalis.Proc. Biol. Sci.271251–258. 10.1098/rspb.2003.2604
40
KatoA.AlbertP. S.VejaJ. M.BirchlerJ. A. (2006). Sensitive fluorescence in situ hybridization signal detection in maize using directly labeled probes produced by high concentration DNA polymerase nick translation.Biotechnic Histoch.8171–78. 10.1080/10520290600643677
41
KubickovaS.CernohorskaH.MusilovaP.RubesJ. (2002). The use of laser microdissection for the preparation of chromosome-specific painting probes in farm animals.Chromosome Res.10571–577. 10.1023/A:1020914702767
42
LégerT.MallyR.NeinhuisC.NussM. (2021). Refining the phylogeny of Crambidae with complete sampling of subfamilies (Lepidoptera, Pyraloidea).Zool. Scripta.5084–99. 10.1111/zsc.12452
43
LewisL. C.LynchR. E. (1969). Rearing the European corn borer, Ostrinia nubilalis (Hubner), on diets containing corn leaf and wheat germ.Iowa State J. Sci.449–14.
44
LoriteP.CarrilloJ. A.AguilarJ. A.PalomequeT. (2004). Isolation and characterization of two families of satellite DNA with repetitive units of 135 bp and 2.5 kb in the ant Monomorium subopacum (Hymenoptera, Formicidae).Cytogenet. Genome Res.10583–92. 10.1159/000078013
45
LuY. J.KochertG. D.IsenhourD. J.AdangM. J. (1994). Molecular characterization of a strain-specific repeated DNA sequence in the fall armyworm Spodoptera frugiperda (Lepidoptera: Noctuidae).Insect Mol. Biol.3123–130. 10.1111/j.1365-2583.1994.tb00159.x
46
LynchR. E. (1980). European corn borer: yield losses in relation to hybrid and stage of corn development.J. Econ. Entomol.73159–164. 10.1093/jee/73.1.159
47
MahendranB.AcharyaC.DashR.GhoshS. K.KunduS. C. (2006). Repetitive DNA in tropical tasar silkworm Antheraea mylitta.Gene37051–57. 10.1016/j.gene.2005.11.010
48
MandrioliM.ManicardiG. C.MarecF. (2003). Cytogenetic and molecular characterization of the MBSAT1 satellite DNA in holokinetic chromosomes of the cabbage moth, Mamestra brassicae (Lepidoptera).Chromosome Res.1151–56. 10.1023/A:1022058032217
49
MarecF.TrautW. (1994). Sex chromosome pairing and sex chromatin bodies in W-Z translocation strains of Ephestia kuehniella (Lepidoptera).Genome37426–435. 10.1139/g94-060
50
Mata-SucreY.SaderM.Van-LumeB.GagnonE.Pedrosa-HarandA.LeitchI. J.et al (2020). How diverse is heterochromatin in the Caesalpinia group? Cytogenomic characterization of Erythrostemon hughesii Gagnon & G.P. Lewis (Leguminosae: Caesalpinioideae).Planta252:49. 10.1007/s00425-020-03453-8
51
MediouniJ.FukováI.FrydrychováR.DhouibiM. H.MarecF. (2004). Karyotype, sex chromatin and sex chromosome differentiation in the carob moth, Ectomyelois ceratoniae (Lepidoptera: Pyralidae).Caryologia57184–194. 10.1080/00087114.2004.10589391
52
MeissleM.MouronP.MusaT.BiglerF.PonsX.VasileiadisV. P.et al (2010). Pests, pesticide use and alternative options in European maize production: current status and future prospects.J. Appl. Entomol.134357–375. 10.1111/j.1439-0418.2009.01491.x
53
MoraP.VelaJ.Ruiz-RuanoF. J.Ruiz-MenaA.MontielE. E.PalomequeT.et al (2020). Satellitome analysis in the ladybird beetle Hippodamia variegata (Coleoptera, Coccinellidae). Genes11:783. 10.3390/genes11070783
54
MunroeE.SolisM. A. (1999). “The Pyraloidea,” in Handbook of Zoology IV: Lepidoptera, Moths and Butterflies, Vol. 1, Arthropoda, Insect, Vol.4, Part 35, ed.KristensenN. (Berlin: Walter de Gruyter & Co.), 491.
55
NacamboS.LeuthardtF. L. G.WanH.LiH.HayeT.BaurB.et al (2013). Development characteristics of the box-tree moth Cydalima perspectalis and its potential distribution in Europe.J. Appl. Entomol.13814–26. 10.1111/jen.12078
56
NavrátilováA.KoblížkováA.MacasJ. (2008). Survey of extrachromosomal circular DNA derived from plant satellite repeats.BMC Plant Biol.8:90. 10.1186/1471-2229-8-90
57
NguyenP.SaharaK.YoshidoA.MarecF. (2010). Evolutionary dynamics of rDNA clusters on chromosomes of moths and butterflies (Lepidoptera).Genetica138343–354. 10.1007/s10709-009-9424-5
58
NguyenP.SýkorováM.ŠíchováJ.KůtaV.DalíkováM.Čapková FrydrychováR.et al (2013). Neo-sex chromosomes and adaptive potential in tortricid pests.Proc. Natl. Acad. Sci. U S A1106931–6936. 10.1073/pnas.1220372110
59
NovákP.NeumannP.MacasJ. (2010). Graph-based clustering and characterization of repetitive sequences in next-generation sequencing data.BMC Bioinform.11:378. 10.1186/1471-2105-11-378
60
NovákP.NeumannP.PechJ.SteinhaislJ.MacasJ. (2013). RepeatExplorer: a Galaxy-based web server for genome-wide characterization of eukaryotic repetitive elements from next-generation sequence reads.Bioinformatics29792–793. 10.1093/bioinformatics/btt054
61
NovákP.RobledilloL. A.KoblížkováA.VrbováI.NeumannP.MacasJ. (2017). TAREAN: a computational tool for identification and characterization of satellite DNA from unassembled short reads.Nucleic Acids Res.45:e111. 10.1093/nar/gkx257
62
NussM.LandryB.MallyR.VeglianteF.TranknerA.BauerF.et al (2003–2020). Global informationsystem on Pyraloidea. Available online at: http://www.pyral oidea.org. (Accessed December 16, 2020).
63
Palacios-GimenezO. M.DiasG. B.LimaL. G.KuhnG. C. S.RamosE.MartinsC.et al (2017). High-throughput analysis of the satellitome revealed enormous diversity of satellite DNAs in the neo-Y chromosome of the cricket Eneoptera surinamensis.Sci. Rep.7:6422. 10.1038/s41598-017-06822-6828
64
Palacios-GimenezO. M.MilaniD.SongH.MartíD. A.López-LeónM. D.Ruiz-RuanoF. J.et al (2020). Eight million years of satellite DNA evolution in grasshoppers of the genus Schistocerca illuminate the ins and outs of the library hypothesis.Genome Biol. Evol.1288–102. 10.1093/gbe/evaa018
65
PalomequeT.LoriteP. (2008). Satellite DNA in insects: a review.Heredity100564–573. 10.1038/hdy.2008.24
66
PitaS.PanzeraF.MoraP.VelaJ.CuadradoA.SánchesA.et al (2017). Comparative repeatome analysis on Triatoma infestans Andean and non-Andean lineages, main vector of Chagas disease.PLoS One12:e0181635. 10.1371/journal.pone.0181635
67
PonsJ.PetitpierreE.JuanC. (1993). Characterization of the heterochromatin of the darkling beetle Misolampus goudoti: cloning of two satellite DNA families and digestion of chromosomes with restriction enzymes.Hereditas119179–185. 10.1111/j.1601-5223.1993.00179.x
68
PringleE. G.BaxterS. W.WebsterC. L.PapanicolaouA.LeeS. F.JigginsC. D. (2007). Synteny and chromosome evolution in the Lepidoptera: evidence from mapping in Heliconius melpomene.Genetics177417–426. 10.1534/genetics.107.073122
69
RegierJ. C.MitterC.SolisM. A.HaydenJ. E.LandryB.NussM.et al (2012). A molecular phylogeny for the pyraloid moths (Lepidoptera: Pyraloidea) and its implications for higher-level classification.Syst. Entomol.37635–656. 10.1111/j.1365-3113.2012.00641.x
70
RobinsonR. (1971). Lepidoptera Genetics.Oxford: Pergamon.
71
SaharaK.MarecF.EickhoffU.TrautW. (2003). Moth sex chromatin probed by comparative genomic hybridization (CGH).Genome46339–342. 10.1139/G03-003
72
SaharaK.YoshidoA.ShibataF.Fujikawa-KojimaN.OkabeT.Tanaka-OkuyamaM.et al (2013). FISH identification of Helicoverpa armigera and Mamestra brassicae chromosomes by BAC and fosmid probes.Insect Biochem. Mol. Biol.43644–653. 10.1016/j.ibmb.2013.04.003
73
SaharaK.YoshidoA.TrautW. (2012). Sex chromosome evolution in moths and butterflies.Chromosome Res.2083–94. 10.1007/s10577-011-9262-z
74
SaitohK. (1959). The chromosome numbers of some species of moths.Jap. J. Genet.3484–87.
75
Serrano-FreitasE. A.SilvaD. M. Z. A.Ruiz-RuanoF. J.UtsunomiaR.Araya-JaimeC.OliveiraC.et al (2020). Satellite DNA content of B chromosomes in the characid fish Characidium gomesi supports their origin from sex chromosomes.Mol. Genet. Genom.295195–207. 10.1007/s00438-019-01615-2
76
ŠíchováJ.NguyenP.DalíkováM.MarecF. (2013). Chromosomal evolution in tortricid moths: conserved karyotypes with diverged features.PLoS One8:e64520. 10.1371/journal.pone.0064520
77
ŠíchováJ.OhnoM.DincăV.WatanabeM.SaharaK.MarecF. (2016). Fissions, fusions, and translocations shaped the karyotype and multiple sex chromosome constitution of the northeast-Asian wood white butterfly, Leptidea amurensis.Biol. J. Linnean Soc.118457–471. 10.1111/bij.12756
78
ŠíchováJ.VoleníkováA.DincăV.NguyenP.VilaR.SaharaK.et al (2015). Dynamic karyotype evolution and unique sex determination systems in Leptidea wood white butterflies.BMC Evol. Biol.15:89. 10.1186/s12862-015-0375-4
79
SilvaB. S. M. L.HeringerP.DiasG. B.SvartmanM.KuhnG. C. S. (2019). De novo identification of satellite DNAs in the sequenced genomes of Drosophila virilis and D. americana using the RepeatExplorer and TAREAN pipelines.PLoS One14:e0223466. 10.1371/journal.pone.0223466
80
SilvaD. M. Z. A.UtsunomiaR.Ruiz-RuanoF. J.DanielS. N.Porto-ForestiF.HashimotoD. T.et al (2017). High-throughput analysis unveils a highly shared satellite DNA library among three species of fish genus Astyanax.Sci. Rep.7:12726. 10.1038/s41598-017-12939-7
81
SmitA. F. A.HubleyR.GreenP. (2013–2015). RepeatMasker Open-4.0. Available online at: http://www.repeatmasker.org. (Accessed Dec 14, 2018).
82
SmithS. G. (1945). Heteropycnosis as a means of diagnosing sex.J. Hered.36195–196. 10.1093/oxfordjournals.jhered.a105498
83
SolisM. A. (1997). “Snout moths: unraveling the taxonomic diversity of a speciose group in the neotropics,” in Biodiversity II: Understanding and Protecting our Biological Resources, edsReaka-KudlaM. L.WilsonD.WilsonE. O. (Washington, D. C: Joseph Henry Press), 231–242.
84
SolisM. A.MetzM. (2016). An illustrated guide to the identification of the known species of Diatraea Guilding (Lepidoptera: Crambidae: Crambinae) based on genitalia.Zookeys.56573–121. 10.3897/zookeys.565.6797
85
SolisM. A.SchefferS. J.LewisM. L.RendonP. (accepted). Diatraea postlineella Schaus (Lepidoptera: Crambidae) from Guatemala: molecular identity and host plant.Proc. Entomol. Soc. Wash.
86
SproulJ. S.KhostD. E.EickbushD. G.NegmS.WeiX.WongI.et al (2020). Dynamic evolution of euchromatic satellites on the X chromosome in Drosophila melanogaster and the simulans clade.Mol. Biol. Evol.372241–2256. 10.1093/molbev/msaa078
87
TallaV.SuhA.KalsoomF.DincăV.VilaR.FribergM.et al (2017). Rapid increase in genome size as a consequence of transposable element hyperactivity in wood-white (Leptidea) butterflies.Genome Biol. Evol.92491–2505. 10.1093/gbe/evx163
88
ThakurR.GautamD. C. (2013). Chromosome studies on four species of moths.Cytologia78327–331. 10.1508/cytologia.78.327
89
TrautW.MarecF. (1996). Sex chromatin in Lepidoptera.Quar. Rev. Biol.71239–256. 10.1086/419371
90
TrautW.MarecF. (1997). Sex chromosome differentiation in some species of Lepidoptera (Insecta).Chromosome Res.5283–291. 10.1023/B:CHRO.0000038758.08263.c3
91
TrautW.SaharaK.MarecF. (2007). Sex chromosomes and sex determination in Lepidoptera.Sex. Dev.1332–346. 10.1159/000111765
92
TrautW.SaharaK.OttoT. D.MarecF. (1999). Molecular differentiation of sex chromosomes probed by comparative genomic hybridization.Chromosoma108173–180. 10.1007/s004120050366
93
TrautW.VogelH.GlöcknerG.HartmannE.HeckelD. G. (2013). High-throughput sequencing of a single chromosome: a moth W chromosome.Chromosome Res.110491–505. 10.1007/s10577-013-9376-6
94
UntergasserA.CutcutacheI.KoressaarT.YeJ.FairclothB. C.RemmM.et al (2012). Primer3-new capabilities and interfaces.Nucleic Acids Res.40:e115. 10.1093/nar/gks596
95
UtsunomiaR.SilvaD. M. Z. A.Ruiz-RuanoF. J.GoesC. A. G.MeloS.RamosL. P.et al (2019). Satellitome landscape analysis of Megaleporinus macrocephalus (Teleostei, Anostomidae) reveals intense accumulation of satellite sequences on the heteromorphic sex chromosome.Sci. Rep.9:5856. 10.1038/s41598-019-42383-8
96
Van NieukerkenE. J.KailaL.KitchingI. J.KristensenN. P.LeesD. C.MinetJ.et al (2011). “Order Lepidoptera Linnaeus, 1758,” in Animal Biodiversity: An Outline of Higher-Level Classification and Survey of Taxonomic Richness, ed.ZhangZ. Q. (Auckland: Magnolia Press).
97
Van’t HofA. E.MarecF.SaccheriI. J.BrakefieldP. M.ZwaanB. J. (2008). Cytogenetic characterization and AFLP-based genetic linkage mapping for the butterfly Bicyclus anynana, covering all 28 karyotyped chromosomes.PLoS One3:e3882. 10.1371/journal.pone.0003882
98
Van’t HofA. E.NguyenP.DalíkováM.EdmondsN.MarecF.SaccheriI. J. (2013). Linkage map of the peppered moth, Biston betularia (Lepidoptera, Geometridae): a model of industrial melanism.Heredity110283–295. 10.1038/hdy.2012.84
99
VěchtováP.DalíkováM.SýkorováM.ŽurovcováM.FüssyZ.ZrzaváM. (2016). CpSAT-1, a transcribed satellite sequence from the codling moth, Cydia pomonella.Genetica144385–395. 10.1007/s10709-016-9907-0
100
VirkkiN. (1963). Gametogenesis in the sugarcane borer moth, Diatraea saccharalis (F.) (Crambidae).J. Agricult. University Puerto Rico47102–137. 10.46429/jaupr.v47i2.12944
101
VítkováM.FukováI.KubíčkováS.MarecF. (2007). Molecular divergence of the W chromosomes in pyralid moths (Lepidoptera).Chromosome Res.15917–930. 10.1007/s10577-007-1173-7
102
VondrakT.RobledilloL. A.NovákP.KoblížkováA.NeumannP.MacasJ. (2020). Characterization of repeat arrays in ultra-long nanopore reads reveals frequent origin of satellite DNA from retrotransposon-derived tandem repeats.Plant J.101484–500. 10.1111/tpj.14546
103
WeiK. H. C.BarbashD. A. (2015). Never settling down: frequent changes in sex chromosomes.PLoS Biol.13:e1002077. 10.1371/journal.pbio.1002077
104
WinnepenninckxB.BackeljauT.de WachterR. (1993). Extraction of high molecular weight DNA from molluscs.Trends Genet.12:407. 10.1016/0168-9525(93)90102-n
105
WrightA. E.DeanR.ZimmerF.MankJ. E. (2016). How to make a sex chromosome.Nat. Commun.7:12087. 10.1038/ncomms12087
106
YasukochiY.OhnoM.ShibataF.JourakuA.NakanoR.IshikawaY.et al (2016). A FISH-based chromosome map for the European corn borer yields insights into ancient chromosomal fusions in the silkworm.Heredity11675–83. 10.1038/hdy.2015.72
107
YoshidoA.MarecF.SaharaK. (2005). Resolution of sex chromosome constitution by genomic in situ hybridization and fluorescence in situ hybridization with (TTAGG)n telomeric probe in some species of Lepidoptera.Chromosoma114193–202. 10.1007/s00412-005-0013-9
108
YoshidoA.ŠíchováJ.PospíšilováK.NguyenP.VoleníkováA.ŠafářJ.et al (2020). Evolution of multiple sex-chromosomes associated with dynamic genome reshuffling in Leptidea wood-white butterflies.Heredity125138–154. 10.1038/s41437-020-0325-9
109
ZhuW.YanJ.SongJ.YouP. (2018). The first mitochondrial genomes for Pyralinae (Pyralidae) and Glaphyriinae (Crambidae), with phylogenetic implications of Pyraloidea.PLoS One13:e0194672. 10.1371/journal.pone.0194672
110
ZrzaváM.HladováI.DalíkováM.ŠíchováJ.ÕunapE.KubíčkováS.et al (2018). Sex chromosomes of the iconic moth Abraxas grossulariata (Lepidoptera, Geometridae) and its congener A. sylvata.Genes9:279. 10.3390/genes9060279
111
ZwyrtkováJ.NěmečkováA.ČížkováJ.HolušováK.KapustováV.SvačinaR.et al (2020). Comparative analyses of DNA repeats and identification of a novel Fesreba centromeric element in fescues and ryegrasses.BMC Plant Biol.20:280. 10.1186/s12870-020-02495-0
Summary
Keywords
Lepidoptera, repetitive DNAs, holocentric chromosomes, tandem repeat, W chromatin
Citation
Cabral-de-Mello DC, Zrzavá M, Kubíčková S, Rendón P and Marec F (2021) The Role of Satellite DNAs in Genome Architecture and Sex Chromosome Evolution in Crambidae Moths. Front. Genet. 12:661417. doi: 10.3389/fgene.2021.661417
Received
30 January 2021
Accepted
04 March 2021
Published
30 March 2021
Volume
12 - 2021
Edited by
Ricardo Utsunomia, Federal Rural University of Rio de Janeiro, Brazil
Reviewed by
André Marques, Max Planck Institute for Plant Breeding Research, Germany; Pablo Mora Ruiz, University of South Bohemia in České Budějovice, Czechia; Duílio M. Z. A. Silva, São Paulo State University, Brazil
Updates
Copyright
© 2021 Cabral-de-Mello, Zrzavá, Kubíčková, Rendón and Marec.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Diogo C. Cabral-de-Mello, cabral.mello@unesp.brFrantišek Marec, marec@entu.cas.cz
This article was submitted to Evolutionary and Population Genetics, a section of the journal Frontiers in Genetics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.