Next Article in Journal
Chromosomal Instability Characterizes Pediatric Medulloblastoma but Is Not Tolerated in the Developing Cerebellum
Previous Article in Journal
Unilateral Cervical Vagotomy Modulates Immune Cell Profiles and the Response to a Traumatic Brain Injury
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 17
10.3390/ijms23179847
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Microsatellite Variation in the Most Devastating Beetle Pests (Coleoptera: Curculionidae) of Agricultural and Forest Crops

by
Manee M. Manee
1,2,*,
Badr M. Al-Shomrani
1,2,
Musaad A. Altammami
1,3,
Hamadttu A. F. El-Shafie
4,
Atheer A. Alsayah
1,
Fahad M. Alhoshani
5 and
Fahad H. Alqahtani
1,2,*
1
National Center for Bioinformatics, King Abdulaziz City for Science and Technology, Riyadh 11442, Saudi Arabia
2
National Center for Agricultural Technology, King Abdulaziz City for Science and Technology, Riyadh 11442, Saudi Arabia
3
Department of Life Sciences, Imperial College London, London SW7 2AZ, UK
4
Date Palm Research Center of Excellence, King Faisal University, Al-Ahsa 31982, Saudi Arabia
5
National Center for Biotechnology, King Abdulaziz City for Science and Technology, Riyadh 11442, Saudi Arabia
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(17), 9847; https://doi.org/10.3390/ijms23179847
Submission received: 25 July 2022 / Revised: 19 August 2022 / Accepted: 24 August 2022 / Published: 30 August 2022
(This article belongs to the Section Molecular Biology)
Browse Figures
Versions Notes

Abstract

:
Weevils, classified in the family Curculionidae (true weevils), constitute a group of phytophagous insects of which many species are considered significant pests of crops. Within this family, the red palm weevil (RPW), Rhynchophorus ferrugineus, has an integral role in destroying crops and has invaded all countries of the Middle East and many in North Africa, Southern Europe, Southeast Asia, Oceania, and the Caribbean Islands. Simple sequence repeats (SSRs), also termed microsatellites, have become the DNA marker technology most applied to study population structure, evolution, and genetic diversity. Although these markers have been widely examined in many mammalian and plant species, and draft genome assemblies are available for many species of true weevils, very little is yet known about SSRs in weevil genomes. Here we carried out a comparative analysis examining and comparing the relative abundance, relative density, and GC content of SSRs in previously sequenced draft genomes of nine true weevils, with an emphasis on R. ferrugineus. We also used Illumina paired-end sequencing to generate draft sequence for adult female RPW and characterized it in terms of perfect SSRs with 1–6 bp nucleotide motifs. Among weevil genomes, mono- to trinucleotide SSRs were the most frequent, and mono-, di-, and hexanucleotide SSRs exhibited the highest GC content. In these draft genomes, SSR number and genome size were significantly correlated. This work will aid our understanding of the genome architecture and evolution of Curculionidae weevils and facilitate exploring SSR molecular marker development in these species.

1. Introduction

The family Curculionidae represents a highly diverse group of coleopteran insects that differ morphologically, ecologically, and behaviorally. Specifically, it comprises 17 subfamilies with over 50,000 described species [1,2]. Members of this family are generally called weevils (snout beetles), and most have a characteristic snout or beak, which is an elongation of the forepart of the head. Curculionidae includes the most damaging and devastating pests of horticultural, field, and forest crops in various ecosystems including rainforests, deserts, and grasslands; these species pose a real menace to global agricultural and forest produce [3,4,5]. For example, the rice weevil, Sitophilus oryzae, can cause 10–80% yield loss [6]. Meanwhile, the mountain pine beetle, Dendroctonus ponderosae, is considered the most important mortality agent for forest ecosystems in western North America and Europe. This weevil seriously influences deforestation and global carbon sequestration strategies [7,8]. Similarly, species of the genus Rhynchophorus, called palm weevils, cause substantial direct damage to several palms of economic importance, such as the edible date palm, oil palm, coconut palm, and the ornamental Canary Islands date palm [5]. They also damage palms indirectly through vectoring diseases or creating wounds that allow the entry of other pathogens [9,10]. Palm weevils also negatively affect the aesthetic value of palms used in urban landscape design [5].
Weevils also comprise extremely important invasive species that may present quarantine problems if they gain entry into new areas, which in modern times is more likely due to the global commercialization and movement of agricultural and forest products [11]. Moreover, it is not easy to detect these weevils during early stages of infestation, making them extremely difficult to control. Nevertheless, it is possible to manage weevils through combining cultural, biological, and chemical strategies in an integrated pest management program. When setting up such control strategies, proper identification and classification of the target beetles is essential to ensure their appropriateness [3]. Recently, [12] reported the first phylogenetic analysis of the subfamily Dryophthoridae within the family Curculionidae which is essential for proper identification and classification.
Microsatellites, also known as simple sequence repeats (SSRs), are 1–6 bp motifs present in both coding and non-coding regions of eukaryotic and prokaryotic genomes that have become the primary source of genetic markers for population analysis in insects due to their high levels of polymorphism [13]. It is well established that SSRs have high rates of mutation and thus have implications for genome organization and genetic variation [14,15]. In addition, SSRs play essential roles in genetic divergence and phenotypic diversity, aiding species in adapting to different environments [16]. Generation of SSR markers by using conventional methods has been challenging; however, in silico mining and analysis of SSRs has proven an effective approach.
To date, draft genome sequences have been released for nine species in the Curculionidae family: R. ferrugineus, Sitophilus oryzae, Hypothenemus hampei, D. ponderosae, Pissodes strobi, Elaeidobius kamerunicus, Ips nitidus, Listronotus oregonensis, and Listronotus bonariensis. This study aimed to identify and characterize microsatellites in the draft genomes of these major agricultural insect pests. The obtained data may contribute to ongoing efforts in managing this group of weevils.

2. Materials and Methods

2.1. Collection of Insect Samples

The female adult of the red palm weevil (RPW) R. ferrugineus used for this study was randomly selected from a colony reared at the insectary of the Date Palm Research Center of Excellence, King Faisal University, Saudi Arabia. The weevil was sexed based on the absence (female) of tuft hairs on the dorsal side of the rostrum [17]. The initial adult weevils used to start the colony were captured in pheromone-food baited traps deployed in an infested date palm plantation in Al-Ahsa, Saudi Arabia (Latitude: 25.268528 ∘ N, Longitude: 49.707218 ∘ E). The weevil colony has been kept for at least three generations, feeding on sugar cane and bolts of the popular “Khalas” date palm cultivar.

2.2. Sample Preparation and DNA Extraction

Tissue (20–30 mg) was obtained from adult female RPW for DNA extraction. Lysis buffer (600 μ L) consisting of 10 mM Tris-HCl, 400 mM NaCl, 100 mM EDTA, pH 8.0, 40 μ L 10% SDS, and 10 μ L Proteinase K (Qiagen, cat. no. 19131; Hilden, Germany) was added to the tissue and incubated overnight, after which the sample was centrifuged and the supernatant discarded. Pellets were resuspended in 1 mL PBS, then processed for DNA extraction and purification by using the KingFisherTM Flex Purification System (ThermoFisher Scientific, cat. no. 5400610; Waltham, MA, USA) and MagMAXTM DNA Multi-Sample Ultra 2.0 Kit (Applied Biosystems, cat. no. A36570; Waltham, MA, USA). The obtained DNA was quantified by using the Qubit dsDNA BR Assay Kit (Invitrogen, cat. no. Q32850; Waltham, MA, USA).

2.3. Next-Generation Sequencing and Genome Assembly

Whole-genome sequencing was outsourced to Macrogen (South Korea) and used paired-end sequencing with read length 151 nucleotides. Library preparation was carried out by using a TruSeq Nano DNA kit according to the sample library preparation protocol (Part # 15041110 Rev. D) on an Illumina NovaSeq 6000 System. De novo assembly was carried out by using SPAdes v3.13.1 with k-mer sizes of 21, 33, 55, and 77 [18]. QUAST v5.2.0 was used to assess the draft assembly metrics [19]. Draft genome completeness was evaluated with the Benchmarking Universal Single-Copy Orthologs (BUSCO) v4.0.6 [20] and the Arthropoda gene set (1013 genes).

2.4. Genome Sequences

The draft genome sequences of nine crop pests were selected for analysis of SSR distributions at genome level. These sequences were assembled at scaffold level according to the genomic resources of the NCBI. The genome sequences in FASTA format were obtained from the Genomes FTP site (ftp://ftp.ncbi.nlm.nih.gov/genomes/ (accessed on 16 May 2022)) and had the following accession numbers: GCA_012979105.1 (male RPW), GCA_014462685.1 (RPW larva), GCA_002938485.2 (S. oryzae), GCA_013372445.1 (H. hampei), GCA_020466585.1 (female D. ponderosae), GCA_020466635.1 (male D. ponderosae), GCA_016904865.1 (P. strobi), GCA_014849505.1 (E. kamerunicus), GCA_018691245.1 (I. nitidus), GCA_019359885.1 (L. oregonensis), and GCA_014170235.1 (L. bonariensis). Although unknown at the time of sequencing, the sex of the RPW larva sample was inferred to be female after analysis of male/female coverage ratios.
The completeness of the assemblies was assessed in relation to BUSCO v4.0.6 [20] based on the Arthropoda database (1013 genes). When investigating the distribution of SSRs in different genomic regions, only three draft genomes and corresponding GFF annotation files could be used: the R. ferrugineus larva and D. ponderosae male and female specimens. We also included the GFF file of Tribolium castaneum (red flour beetle, family Tenebrionidae) for comparison purposes.

2.5. Identification of Microsatellites

The software PERF v0.2.5 [21] was used to scan each entire genome and conduct genome-wide SSR mining. The following criteria were adopted to identify perfect SSRs: repeat lengths of 1 to 6 nucleotides and minimum repeat numbers of 12 repeats for mononucleotides, 7 repeats for dinucleotides, 5 repeats for trinucleotides, and 4 repeats for tetra-, penta- and hexanucleotides; these criteria are consistent with previous studies [22,23]. The remaining parameters were set as default. Repeats with unit patterns being circular permutations and/or reverse complements were deemed as a single type in this study [24,25]; for instance, depending on the reading frame and strand, the type “ACT” corresponds to ACT, CTA, TAC, ATG, GAT, and TGA. The relative frequency (number per Mb) and relative density (length in bp per Mb) of SSRs were utilized in comparing different types of SSR repeats or motifs.

2.6. Assigning Microsatellites to Genomic Regions

We determined exon sequences and gene coding sequences (CDSs) of the nine weevil genomes in this study according to the positions noted in genome annotation files in general feature format (GFF). Intergenic regions were defined as the interval sequences between two adjacent genes. Intronic regions were defined as interval sequences within genes that did not overlap any annotated exons. We identified the coordinates defining intergenic and intronic regions from GFF files by using the BEDtools subtract tool v2.30.0, and assigned the identified perfect SSRs to genomic compartments by using the BEDtools intersect tool v2.30.0 [26].

2.7. Statistical Analysis

All graphical and statistical analyses were carried out in the R programming environment (v4.0.4) (R Core Team, 2021). Pearson correlations determined by using the cor.test method were utilized to elucidate correlations between SSR data sets, including in terms of the number, relative frequency, relative density, and GC content of SSRs.

3. Results

3.1. Genome Assembly and Assessing of Draft Genome Completeness

The de novo assembly of female RPW was performed, generating a draft genome of 1121.36 Mb with a GC content of 43.96%. Contigs with lengths less than 200 bp were filtered out prior to the analysis. The final draft assembly resulted in 945,214 contigs that yielded the longest contig length of 720,101 bp with an N50 contig length of 7782 bp. To determine the completeness of each weevil genome assembly including our female RPW draft, we compared it against the BUSCO Arthropoda lineage dataset (arthropoda_odb10), which consisted of 1013 single-copy orthologs. This revealed that for eight of the sequenced species, 72.4–97.4% of those 1013 Arthropoda single-copy orthologs were completely present; the exception were R. ferrugineus adult male and E. kamerunicus, at 52.9% and 51%, respectively (Figure 1).

3.2. Identification and Characterization of Microsatellites in Beetle Genomes

Twelve draft genomes representing the insect species R. ferrugineus, S. oryzae, H. hampei, D. ponderosae, P. strobi, E. kamerunicus, I. nitidus, L. oregonensis, and L. bonariensis were scanned for perfect microsatellites by using PERF. We first carried out analyses to report all perfect SSRs in the RPW genomes without applying any search criteria (Supplementary Files S1–S3). All exhibited similar patterns of SSRs, as shown in Figure 2. When applying consistent search parameters, a total of 57,175, 50,723, and 67,261 perfect SSRs were identified with frequencies ranging from 50.99 to 114.11 SSRs/Mb in the adult female, adult male, and larval RPW genomes, respectively (Table 1). These perfect SSRs occupied about 0.13%, 0.14%, and 0.36% of the respective genome, had mean lengths of 25.91, 22.29, and 31.98 bp, and their relative densities ranged from 1320.92–3649.45 SSRs/Mb. The other true weevil genomes exhibited similar length proportions for their SSRs, ranging from 0.02% (E. kamerunicus) to 1.44% (L. oregonensis), as seen in Table 1. Number of SSRs was positively correlated with their relative frequency and density (Pearson r = 0.944, p < 0.01 and Pearson r = 0.937, p < 0.01, respectively). The genome size of these draft genomes was also significantly positively correlated with number of SSRs (Pearson r < 0.580, p < 0.05). In contrast, the GC content of SSRs was not significantly correlated with number of SSRs (Pearson r < −0.442, p = 0.150). The relative frequency and density of SSRs were also not significantly correlated with genome size (Pearson r < 0.370, p = 0.236 and Pearson r < 0.324, p = 0.305, respectively). For example, P. strobi has the largest genome (2025.02 Mb) among those surveyed, but was found to have lower SSR frequency (76.30 SSRs/Mb) compared to some other species with smaller genome sizes (Table 1).
Table 2 lists the respective number, length, relative frequency, relative density, and percentage of each of the six types of SSRs. The percentage and relative frequencies and densities of different SSR types were found to vary in the twelve draft genomes (Figure 3). Dinucleotide SSRs were the most frequent type in the R. ferrugineus adult male, adult female, and larva and in I. nitidus, with respective frequencies of 30.93, 35.52, 81.43, and 57.18 SSRs/Mb; these accounted for 60.66%, 54.77%, 71.35%, and 40.78% of SSRs in those draft genomes (Figure 3A,B). Meanwhile, mononucleotide SSRs were the most abundant type in S. oryzae, P. strobi, L. oregonensis, and L. bonariensis, with respective frequencies of 29.18, 23.87, 221.81, and 51.89 SSRs/Mb and comprising 26.64%, 31.29%, 53.71%, and 36.83% of all SSRs (Figure 3A). Trinucleotide SSRs were the most frequent type in H. hampei and in both female and male D. ponderosae, with frequencies of 24.78, 11.43, and 11.85 SSRs/Mb. Finally, tetranucleotide SSRs were the most abundant type in E. kamerunicus, with a frequency of 5.60 SSRs/Mb and accounting for 34.34% of SSRs.
Dinucleotide SSRs were found to have the highest densities, ranging from 956.40 to 10,152.94 bp/Mb in R. ferrugineus, S. oryzae, H. hampei, I. nitidus, L. oregonensis, and L. bonariensis (Figure 3C). Trinucleotide SSRs had the highest densities (198.78–561.66 bp/Mb) in D. ponderosae and P. strobi, whereas tetranucleotide SSRs had the highest density (91.25 bp/Mb) in E. kamerunicus (Figure 3C). Across the investigated genomes, hexanucleotide SSRs were the least abundant at frequencies below 1.93 SSRs/Mb, except in L. oregonensis, for which pentanucleotide SSRs were identified to be the least frequent (1.08 SSRs/Mb).
Next, GC content was investigated for the various types of SSRs (Figure 3D). The highest GC content was observed for hexanucleotide SSRs, which had values of 19.48–54.34%, except in P. strobi, for which genome mononucleotide SSRs exhibited the highest GC content at 43.26%. Meanwhile, the lowest levels of GC content were identified for dinucleotide SSRs in S. oryzae, R. ferrugineus, and L. oregonensis, at values of only 0.55–4.41%; for mononucleotide SSRs in H. hampei, E. kamerunicus, L. bonariensis, and D. ponderosae, at 0.01–12.24%; and for trinucleotide SSRs in P. strobi, at 8.82%.

3.3. Diversity of Microsatellite Motifs in Beetle Genomes

The microsatellites in the weevil genome assemblies examined here were found to be relatively AT-rich. To gain insight into this characteristic, we further analyzed the motif composition of SSRs. Motif abundance was found to vary across the draft genomes. More specifically, the investigated assemblies were identical in the degenerated number of repeat motifs for mono- to trinucleotide SSRs, at 2, 4, and 10 motifs respectively, but differed in the number of tetranucleotide, pentanucleotide, and hexanucleotide repeat motifs.
Among mononucleotide repeats, the predominate motif was (A)n, with total counts of 4385, 4951, 4228, 20603, 3798, 1305, 1215, 28851, 913, 7138, 283468, and 57712 SSRs in R. ferrugineus (F), R. ferrugineus (M), R. ferrugineus (L), H. hampei, D. ponderosae (F), D. ponderosae (M), P. strobi, E. kamerunicus, I. nitidus, L. oregonensis, and L. bonariensis respectively. This type accounted for 6.29–53.07% of all mononucleotide SSRs in the draft genomes (Figure 4). The frequency of the (A)n motif ranged from 3.39–219.19% SSRs/Mb, with the highest frequency observed in L. oregonensis and the lowest in E. kamerunicus. The (C)n motif type was far less abundant, accounting for just 0.01–12.62% of all mononucleotide SSRs in the twelve draft genomes.
Among dinucleotide SSRs, the most prominent type in ten draft genomes was the (AT)n motif, with frequencies ranging from 1.14 to 142.38 SSRs/Mb; the exceptions were H. hampei and I. nitidus, in which this motif comprised about 6.98–67.17% of dinucleotide SSRs (Figure 4). In H. hampei, the most frequent dinucleotide motif was the (AG)n repeat at 7.47 SSRs/Mb, accounting for 9.28% of all SSRs in that assembly. Meanwhile, in I. nitidus, the most prevalent dinucleotide motif was (AC)n with frequency 25.72 SSRs/Mb; this motif accounted for 18.35% of all dinucleotide SSRs in that genome. Notably, the (AG)n repeat was almost equally frequent in I. nitidus (24.67 SSRs/Mb). In all weevil assemblies, the least frequent dinucleotide SSR was the (CG)n motif.
For the trinucleotide repeat type, the (AAT)n repeat was the most frequent motif in eleven draft genomes, with frequencies ranging from 3.36 to 15.20 SSRs/Mb; these repeats accounted for 3.24–19.92% of all trinucleotide SSRs Figure 4). The exception was I. nitidus, in which the (AAC)n repeat was the most frequent trinucleotide motif, followed by the (AAT)n motif; these had frequencies of below 9 SSRs/Mb, and together accounted for 11.24% of all trinucleotide SSRs in that species.
Among tetranucleotide repeats, (AAAT)n was the most abundant in eleven assemblies with frequencies ranging from 1.69 to 10.65 SSRs/Mb and accounting for 1.81–19.06% of all tetranucleotide SSRs. The exception was again I. nitidus (Figure 4), in which the most frequent tetranucleotide motif was (AAAG)n, with frequency 2.66 SSRs/Mb and comprising about 1.89% of all tetranucleotide SSRs in that draft genome.
For pentanucleotide repeats, the most abundant motifs varied among species. (AAACC)n was the most abundant in the S. oryzae, with frequency of 2.68 SSRs/Mb and comprising about 2.45% of pentanucleotide SSRs in this draft genome. The (AATAT)n motif was the most frequent in the R. ferrugineus adult female, R. ferrugineus adult male and R. ferrugineus larva with frequencies of 0.19, 0.33, and 0.43 SSRs/Mb, respectively. Meanwhile, (AAATC)n and (AAATT)n motifs had similar frequencies of approximately 0.15 SSRs/Mb in the D. ponderosae adult female and male assemblies, accounting for 1.92% of all pentanucleotide SSRs. (AACCT)n repeats were the predominant pentanucleotide motif in H. hampei and I. nitidus, with frequencies below 3 SSRs/Mb. (ACGAG)n and (AATCT)n motif types were more abundant in the L. oregonensis, and L. bonariensis, with respective frequencies of 1.22 and 1.90 SSRs/Mb. Finally, P. strobi and E. kamerunicus were found to share their most frequent pentanucleotide motif, (AAATC)n, with a frequency below 0.6 SSRs/Mb.
Hexanucleotide motifs occurred at a far lower frequency in the examined weevil genomes than did other microsatellite repeat types. The (AAACCC)n motif was the most abundant hexanucleotide in the R. ferrugineus adult female and R. ferrugineus larva draft genomes, with frequencies of less than 0.07 SSRs/Mb, while the (ACATAT)n repeat was the most frequent in the R. ferrugineus adult male, with the frequency of 0.03 SSRs/Mb. The (AAATTC)n motif was the most frequent type in D. ponderosae, P. strobi, E. kamerunicus, and L. bonariensis, with frequencies below 0.4 SSRs/Mb. Meanwhile, (AAGAGG)n, (ACACAT)n, (AAAGAG)n, and (AAGACC)n motifs were the most abundant hexanucleotide repeats in S. oryzae, H. hampei, I. nitidus, and L. oregonensis, respectively.

3.4. Microsatellite Distribution and Motif Diversity According to Genomic Region

The distribution of SSRs across different genomic regions was investigated in four draft genomes representing three species (R. ferrugineus larva, female and male D. ponderosae, and T. castaneum) as described in the Methods. Specifically, microsatellite analysis was executed to examine the distribution of SSRs in exons, CDSs, and intronic and intergenic regions. The results revealed most mono- to hexanucleotide SSRs to have region-associated differences in terms of their relative abundance, density, and percentage, and those differences to vary between species; however, as expected, results in the female and male D. ponderosae were substantially similar. Overall, lower relative frequencies and densities of SSRs were observed in coding and noncoding regions than in intronic and intergenic regions (Figure 5). Microsatellites were most commonly identified in intergenic regions, followed in order by intronic regions, exons, and CDSs, with one exception: SSRs were found to be abundant in the intronic regions of T. castaneum (Figure 5B). In CDSs of the four assemblies, SSR frequency ranged from 0.95 to 4.97 SSRs/Mb; overall, coding regions contained 0.83–5.54% of SSRs. In exons, SSR frequency ranged from 0.95 to 3.90 SSRs/Mb except in T. castaneum, which had a frequency 8.03 SSRs/Mb; collectively, exonic regions accounted for 0.83–12.44% of SSRs in the four samples. In intronic regions of R. ferrugineus larva, female and male D. ponderosae, and T. castaneum, respectively, the observed SSR frequencies were 26.47, 8.99, 9.36, and 44.50 SSRs/Mb; in total, introns accounted for 22.98–27.98% of SSRs except in T. castaneum, where they comprised 48.71%. Finally, intergenic regions exhibited respective frequencies of 26.47, 8.99, 9.36, and 44.50 SSRs/Mb, and accounted for 37.05–75.37% of SSRs in the four assemblies. Overall, microsatellite densities were higher in noncoding regions than in coding regions: intronic regions had densities of 146.19–1063.02 bp/Mb, and intergenic regions of 289.79–2791.55 bp/Mb, while CDSs had densities of 15.58–119.29 bp/Mb and exons of 15.58–177.18 bp/Mb (Figure 5C).
Next, the GC content of microsatellites was examined according to genomic region (Figure 5C). Across the four assemblies, GC contents were mostly identical in coding regions (CDSs and exons), but were found to vary in noncoding regions (intronic and intergenic regions). The highest GC contents were observed for SSRs located in CDSs (48.88–52.85%), followed by those in exons (32.98–51.66%), whereas intronic regions had GC contents of 3.31–17.16% and intergenic regions of 3.51–18.61%.
Among CDSs and exons, trinucleotide SSRs were the most abundant type (0.77–5.35 SSRs/Mb) in all four genomes, while pentanucleotide SSRs were consistently the least frequent in the three curculionid assemblies (Figure 6A,B). For the tenebrionid T. castaneum, di- and hexanucleotide SSRs were the least abundant types in CDSs (0.07 SSRs/Mb) and exons (0.24 SSRs/Mb), respectively. In intronic and intergenic regions, trinucleotide SSRs were the most abundant type in D. ponderosae and T. castaneum, with frequencies of 17.94–2.99 SSRs/Mb, whereas dinucleotide SSRs were the most abundant type in R. ferrugineus (Figure 6C,D). Pentanucleotide SSRs were rare in intronic and intergenic regions, and hexanucleotide SSRs were the least abundant, with frequencies below 1.08 SSRs/Mb for all four genomes (Figure 6C,D).
Among the three beetle species examined here, motif types were found to vary quite obviously in different genomic regions (Figure 7). In coding regions of R. ferrugineus and T. castaneum, the predominant motifs were (AAG)n and (CCG)n, respectively, accounting for 15–22% of CDS and exonic SSRs (Figure 7A,B). Meanwhile, (AGC)n and (AAT)n respectively comprised the most abundant trinucleotide repeats in the CDSs and exonic regions of D. ponderosae. In noncoding regions of the R. ferrugineus genome, the (AT)n motif was the most abundant repeat, representing ∼67% of intronic and intergenic SSRs (Figure 7C,D). Meanwhile, intronic and intergenic regions of the T. castaneum assembly had (AAT)n as the most common repeat, with frequencies of approximately 16 SSRs/Mb. In D. ponderosae assemblies, (A)n and (AAT)n were the most abundant motifs in intronic regions and intergenic regions, with frequencies below 4 SSRs/Mb.

4. Discussion

The development of next-generation sequencing has allowed for the generation of a massive number of sequenced draft genomes, including those of non-model species. The availability of draft genomic sequences from Curculionidae weevils allowed us to investigate the distributions of microsatellites in members of this family. As far as we know, this is the first comprehensive report on the identification and analysis of SSRs 1–6 bp long in the entire draft genomes of nine curculionid beetles. We used computational techniques to search for microsatellites and compare the relative frequency, relative density, and GC content of SSRs in these beetles. Consistent search parameters were utilized so as to carry out the same analysis in each investigated draft genome. BUSCO results suggest these draft genomes are mostly comparable. Moreover, BUSCO indicated that our female RPW assembly is more complete than male RPW (GCA_012979105.1) [27] both in terms of complete single genes (92.2% versus 52.9%, respectively) and of missing genes (3.0% versus 15.13%, respectively). SSR repeat content differs between species, which might be a general phenomenon across taxa [33]. Previous studies reported SSRs to comprise 3% of the human genome [34], 0.04–0.44% of plant and fungal genomes [35,36,37], and 0.44–0.88% of primate genomes [22,38]. Here, our results showed that identified SSRs differ with the degree of coverage and comprise 0.02–1.44% of the draft genomes for these nine weevil species. Assemblies representing the same species exhibited similar proportions of SSRs, as seen in female and male D. ponderosae, whereas values differed between species. The observed variance in microsatellite proportion could result from differences in computational approaches utilized for SSR detection, incompleteness of genome assemblies, or actual variation in SSR content among these weevils [39]. Moreover, variation might even arise between closely related species [40,41].
Our findings suggest that in weevils, the number of SSRs is significantly positively correlated with genome size; this is inconsistent with the results reported in [35,42]. Nonetheless, a study reported that the number of SSRs was significantly associated with genome size in 136 insects [43], which agrees with our results. However, it is necessary to sequence more genomes of beetles from the Curculionidae family to solidify this conclusion. In this work, frequency and density of SSRs were not significantly correlated with genome size.
In all of the weevil species examined here, the six types of SSRs were not evenly distributed; rather, mono- to trinucleotide SSRs were the most prevalent. This finding is consistent with previous reports that mono- to trinucleotide SSR repeats are more frequent in 23 mosquito species [44] and six plant species [45]. Meanwhile, tetra- to hexanucleotide SSRs were the least frequent types in these draft genomes, an observation similar to what has been found in Palmae genomes [35] and Gossypium species [46]. More specifically, we observed dinucleotide SSRs to be the most frequent repeat type in R. ferrugineus and I. nitidus, consistent with dicotyledons [47] and Drosophila [14]. Mononucleotide SSRs were the dominant type in S. oryzae, P. strobi, L. oregonensis, and L. bonariensis, which is consistent with prior findings for Batocera horsfieldi [48] and eukaryotic genomes [39,49]. Finally, trinucleotide SSRs were the most abundant type in H. hampei, which is consistent with eukaryotes [50]. The higher abundance of SSRs with shorter motif lengths (mono-, di-, and trinucleotides) could be the result of a higher frequency of replication slippage over shorter repeat monomers. Additionally, repeat motifs may differ in the stability of secondary structures they form, which might also impact the evolutionary dynamics of their abundance and distribution [51]. However, no such analysis has been performed in weevils and the relative contributions of selection and the molecular mechanisms affecting the abundance of SSRs (e.g., slippage, rolling circle amplification, crossing over, gene conversion) is poorly understood in general [52].
We also observed SSR motifs within each microsatellite type to vary in abundance across the examined draft genomes. Among mononucleotide repeats, the most frequent motif was (A/T)n, occupying about 6.29–53.07% of mononucleotide SSRs in these genomes, similar to the trend previously reported across 100 insect species [43]. Of dinucleotide SSRs, the most abundant motifs were (AT)n and (AG)n, similar to palms [35], several insect species [43], and garden asparagus [53]. Regarding trinucleotide motifs, (AAT)n was the dominant motif in most weevil draft genomes, which is consistent with both mammals [22,49] and plants [35,54]. Of tetra-, penta-, and hexanucleotide SSRs, (AAAT)n, (AAAG)n, (AAACC)n, (AAATC)n, (AAATT)n, (AATAT)n, (AATCT)n, (ACATAT)n, and (AAATTC)n were the more frequent motifs. Overall, these findings are consistent with previous reports suggesting that AT-rich SSR motifs predominate [43,48]. The abundance of AT-rich SSRs might also reflect the overall base composition of insect genomes, which are often AT-rich themselves [43].
Strong evidence exists that the microsatellites are nonrandomly distributed across protein-coding regions, untranslated regions, and introns, and that they may play roles in gene expression and regulation [55,56]. Moreover, SSRs may play different functional roles in different genomic regions. We further investigated the distribution of SSRs in different genomic regions for four beetles from three species representing two beetle families (Curculionidae and Tenebrionidae). We found SSR abundance to differ among genomic regions in these genomes; moreover, the same genomic regions in different species showed notable similarity in SSR distribution, consistent with previous studies in mammals and plants [49,57]. SSRs were found to occur less frequently in coding regions than in noncoding regions, which aligns with previous reports [58,59]. Specifically, SSRs were greatly abundant in intergenic and intronic regions, less common in exons, and least abundant in CDSs. These results may suggest that SSRs in coding regions are subject to negative/purifying selection pressure [23].
Within CDSs and exons, trinucleotide SSRs were the most abundant repeat type, which echoes results from prior studies in mammals and plants [23,49]. The predominance of trinucleotide SSRs in coding regions may be due to frameshift mutations eliminating non-trimeric SSRs [60]. Inconsistent with previous reports in mosquitos, primates, mammals, and plants [35,44,49,58], we observed intronic and intergenic regions to feature trinucleotide SSRs as the most abundant repeat type in D. ponderosae and T. castaneum, but dinucleotide SSRs in R. ferrugineus.
Notably, SSRs exhibit bias toward a few specific nucleotide motifs according to the genomic region they occur in. In coding regions of the R. ferrugineus and T. castaneum genomes, (AAG)n and (CCG)n repeats predominated; meanwhile, there was a noticeable excess of the (AGC)n motif in the CDSs and exonic regions of D. ponderosae, similar to observations in Drosophila [14]. Consistent with previous reports [35,58,59], AT-rich motifs such as (AT)n, (AAT)n, and (A)n were the most abundant in the intronic and intergenic regions of the examined beetle genomes, which can be interpreted as confirming high AT content in the majority of the analyzed SSRs.
To evaluate the effects of nucleotide composition on SSR abundance, we examined GC content in relation to SSR type in the different genomic compartments of all nine weevil species. The results showed average GC content values (0.01–54.34%) to be much lower than AT content values, and moreover that the distribution of GC content was uneven; this is consistent with previous reports [35,49,58,59]. The greatest GC content values were mostly detected among hexanucleotide SSRs and the least for mono- and dinucleotide SSRs. In terms of genomic regions, CDSs demonstrated the most GC content, followed by exonic regions, then intergenic regions, and lastly intronic regions. These results suggest that high GC content is more frequently distributed in coding regions, consistent with results reported in [58]. The bias for high GC in coding regions has been suggested to increase the bendability of the double helix [61] and in turn contribute to maintain the higher transcriptional activity in these regions [62].
All told, this study performed the first comprehensive large-scale analysis of microsatellites in draft genomes of nine crop pests of the Curculionidae family, with a focus on common features of SSRs including their abundance patterns and variation characteristics. The findings of this work provide useful insights into the diversity and distributions of SSRs in these weevil species. The SSR number in these draft genomes was significantly correlated with genome size and but not significantly correlated with GC content. Mono- to trinucleotide SSRs were dominant in all examined species, but the occurrence, percentage, and density of each type of SSR varied between species. Overall, most SSRs were distributed in intronic and intergenic regions; within coding regions, trinucleotide SSRs predominated. Genomic microsatellite markers are widely used in population genetics and evolutionary studies because they are reliable, highly polymorphic, and easy to amplify [63]. Further refining our understanding of the characteristics of SSRs in weevil genomes will serve as a foundation for genetic research and the selection of SSR molecular markers in these beetles.

Supplementary Materials

The supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms23179847/s1.

Author Contributions

Conceptualization, M.M.M.; methodology, M.M.M.; software, M.M.M. and F.H.A.; validation, M.M.M. and F.H.A.; formal analysis, M.M.M., B.M.A.-S., M.A.A., H.A.F.E.-S., A.A.A., F.M.A., and F.H.A.; investigation, M.M.M.; resources, M.M.M. and F.H.A.; data curation, M.M.M., B.M.A.-S., and F.H.A.; writing—original draft preparation, M.M.M., M.A.A., H.A.F.E.-S., and F.H.A.; writing—review and editing, M.M.M.; visualization, M.M.M.; supervision, M.M.M.; project administration, M.M.M. and F.H.A.; funding acquisition, M.M.M. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Life Science and Environment Research Institute and the Center of Excellence for Genomics (grant 20-0078), King Abdulaziz City for Science and Technology, Saudi Arabia.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data generated and analyzed during this study are included in the published article, its Supplementary Files, and publicly available repositories. Raw reads from genome sequencing of the female R. ferrugineus have been deposited at NCBI Sequence Read Archive (SRA) under the BioProject accessions PRJNA848948. Draft genome of the female R. ferrugineus weevil can be found at https://doi.org/10.5281/zenodo.6878576 (accessed on 21 July 2022).

Acknowledgments

The authors would like to thank Guilherme Dias at the Department of Genetics and Institute of Bioinformatics, University of Georgia, for his valuable comments and suggestions. The authors would also thank Amer S. Alharthi at the General Directorate for Research and Innovation, King Abdulaziz City for Science and Technology, for his technical support.

Conflicts of Interest

The authors declare there are no competing interest.

References

  1. Bozdoğan, H.; Erbey, M.; Aksoy, H.A. Total amount of protein, lipid and carbohydrate of some adult species belong to curculionidae family (Coleoptera: Curculionidae). J. Entomol. Zool. Stud 2016, 4, 242–248. [Google Scholar]
  2. Bhatti, A.R.; Zia, A.; Mastoi, M.I.; Shehzad, M.I.A.; Iqbal, J. Tanymecus xanthuruschevrolat, 1880 (curculionidae: Entiminae), a new addition to curculionid fauna of pakistan. Pak. Entomol 2018, 40, 91–94. [Google Scholar]
  3. Rugman-Jones, P.F.; Hoddle, C.D.; Hoddle, M.S.; Stouthamer, R. The lesser of two weevils: Molecular-genetics of pest palm weevil populations confirm Rhynchophorus vulneratus (Panzer 1798) as a valid species distinct from R. ferrugineus (Olivier 1790), and reveal the global extent of both. PLoS ONE 2013, 8, e78379. [Google Scholar] [CrossRef]
  4. Aguirre, C.; Olivares, N.; Luppichini, P.; Hinrichsen, P. A PCR-based diagnostic system for differentiating two weevil species (Coleoptera: Curculionidae) of economic importance to the chilean citrus industry. J. Econ. Entomol. 2015, 108, 107–113. [Google Scholar] [CrossRef] [PubMed]
  5. Milosavljević, I.; El-Shafie, H.A.; Faleiro, J.R.; Hoddle, C.D.; Lewis, M.; Hoddle, M.S. Palmageddon: The wasting of ornamental palms by invasive palm weevils, Rhynchophorus spp. J. Pest Sci. 2019, 92, 143–156. [Google Scholar] [CrossRef]
  6. Chen, H.; Chen, Z.; Zhou, Y. Rice water weevil (Coleoptera: Curculionidae) in mainland China: Invasion, spread and control. Crop Prot. 2005, 24, 695–702. [Google Scholar] [CrossRef]
  7. Bentz, B.J.; Jönsson, A.M.; Schroeder, M.; Weed, A.; Wilcke, R.A.I.; Larsson, K. Ips typographus and Dendroctonus ponderosae models project thermal suitability for intra-and inter-continental establishment in a changing climate. Front. For. Glob. Chang. 2019, 2, 1. [Google Scholar]
  8. Hansen, E.M.; Amacher, M.C.; Van Miegroet, H.; Long, J.N.; Ryan, M.G. Carbon dynamics in central US Rockies lodgepole pine type after mountain pine beetle outbreaks. For. Sci. 2015, 61, 665–679. [Google Scholar] [CrossRef]
  9. Griffith, R.; Koshy, P. Chapter Il Nematode Parasites of Coconut and Other Paims. Plant Parasit. Nematodes Subtrop. Trop. Agric. 1990, 363. [Google Scholar] [CrossRef]
  10. Cruz, L.F.; Menocal, O.; Mantilla, J.; Ibarra-Juarez, L.A.; Carrillo, D. Xyleborus volvulus (Coleoptera: Curculionidae): Biology and fungal associates. Appl. Environ. Microbiol. 2019, 85, e01190-19. [Google Scholar]
  11. Faleiro, J. A review of the issues and management of the red palm weevil Rhynchophorus ferrugineus (Coleoptera: Rhynchophoridae) in coconut and date palm during the last one hundred years. Int. J. Trop. Insect Sci. 2006, 26, 135–154. [Google Scholar]
  12. Chamorro, M.L.; de Medeiros, B.A.; Farrell, B.D. First phylogenetic analysis of Dryophthorinae (Coleoptera, Curculionidae) based on structural alignment of ribosomal DNA reveals Cenozoic diversification. Ecol. Evol. 2021, 11, 1984–1998. [Google Scholar] [CrossRef] [PubMed]
  13. Ma, L.; Cao, L.J.; Hoffmann, A.A.; Gong, Y.J.; Chen, J.C.; Chen, H.S.; Wang, X.B.; Zeng, A.P.; Wei, S.J.; Zhou, Z.S. Rapid and strong population genetic differentiation and genomic signatures of climatic adaptation in an invasive mealybug. Divers. Distrib. 2020, 26, 610–622. [Google Scholar] [CrossRef]
  14. Katti, M.V.; Ranjekar, P.K.; Gupta, V.S. Differential distribution of simple sequence repeats in eukaryotic genome sequences. Mol. Biol. Evol. 2001, 18, 1161–1167. [Google Scholar] [CrossRef]
  15. Bagshaw, A.T. Functional mechanisms of microsatellite DNA in eukaryotic genomes. Genome Biol. Evol. 2017, 9, 2428–2443. [Google Scholar] [CrossRef] [Green Version]
  16. Kashi, Y.; King, D.G. Simple sequence repeats as advantageous mutators in evolution. Trends Genet. 2006, 22, 253–259. [Google Scholar] [CrossRef]
  17. Kaakeh, W. Longevity, fecundity, and fertility of the red palm weevil, Rynchophorus ferrugineus Olivier (Coleoptera: Curculionidae) on natural and artificial diets. Emir. J. Food Agric. 2005, 23–33. [Google Scholar] [CrossRef]
  18. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Pham, S.; Prjibelski, A.D.; et al. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 2012, 19, 455–477. [Google Scholar] [CrossRef]
  19. Gurevich, A.; Saveliev, V.; Vyahhi, N.; Tesler, G. QUAST: Quality assessment tool for genome assemblies. Bioinformatics 2013, 29, 1072–1075. [Google Scholar] [CrossRef]
  20. Simão, F.A.; Waterhouse, R.M.; Ioannidis, P.; Kriventseva, E.V.; Zdobnov, E.M. BUSCO: Assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 2015, 31, 3210–3212. [Google Scholar] [CrossRef]
  21. Avvaru, A.K.; Sowpati, D.T.; Mishra, R.K. PERF: An exhaustive algorithm for ultra-fast and efficient identification of microsatellites from large DNA sequences. Bioinformatics 2017, 27, 573. [Google Scholar] [CrossRef] [PubMed]
  22. Liu, S.; Hou, W.; Sun, T.; Xu, Y.; Li, P.; Yue, B.; Fan, Z.; Li, J. Genome-wide mining and comparative analysis of microsatellites in three macaque species. Mol. Genet. Genom. 2017, 292, 537–550. [Google Scholar] [CrossRef] [PubMed]
  23. Qi, W.H.; Jiang, X.M.; Yan, C.C.; Zhang, W.Q.; Xiao, G.S.; Yue, B.S.; Zhou, C.Q. Distribution patterns and variation analysis of simple sequence repeats in different genomic regions of bovid genomes. Sci. Rep. 2018, 8, 14407. [Google Scholar] [CrossRef] [PubMed]
  24. Jurka, J.; Pethiyagoda, C. Simple repetitive DNA sequences from primates: Compilation and analysis. J. Mol. Evol. 1995, 40, 120–126. [Google Scholar] [CrossRef]
  25. Li, C.Y.; Liu, L.; Yang, J.; Li, J.B.; Su, Y.; Zhang, Y.; Wang, Y.Y.; Zhu, Y.Y. Genome-wide analysis of microsatellite sequence in seven filamentous fungi. Interdiscip. Sci. Comput. Life Sci. 2009, 1, 141–150. [Google Scholar] [CrossRef] [PubMed]
  26. Quinlan, A.R.; Hall, I.M. BEDTools: A flexible suite of utilities for comparing genomic features. Bioinformatics 2010, 26, 841–842. [Google Scholar] [CrossRef] [PubMed]
  27. Hazzouri, K.M.; Sudalaimuthuasari, N.; Kundu, B.; Nelson, D.; Al-Deeb, M.A.; Le Mansour, A.; Spencer, J.J.; Desplan, C.; Amiri, K. The genome of pest Rhynchophorus ferrugineus reveals gene families important at the plant-beetle interface. Commun. Biol. 2020, 3, 323. [Google Scholar] [CrossRef]
  28. Dias, G.B.; Altammami, M.A.; El-Shafie, H.A.; Alhoshani, F.M.; Al-Fageeh, M.B.; Bergman, C.M.; Manee, M.M. Haplotype-resolved genome assembly enables gene discovery in the red palm weevil Rhynchophorus ferrugineus. Sci. Rep. 2021, 11, 9987. [Google Scholar] [CrossRef]
  29. Parisot, N.; Vargas-Chávez, C.; Goubert, C.; Baa-Puyoulet, P.; Balmand, S.; Beranger, L.; Blanc, C.; Bonnamour, A.; Boulesteix, M.; Burlet, N.; et al. The transposable element-rich genome of the cereal pest Sitophilus oryzae. BMC Biol. 2021, 19, 241. [Google Scholar] [CrossRef]
  30. Vega, F.E.; Brown, S.M.; Chen, H.; Shen, E.; Nair, M.B.; Ceja-Navarro, J.A.; Brodie, E.L.; Infante, F.; Dowd, P.F.; Pain, A. Draft genome of the most devastating insect pest of coffee worldwide: The coffee berry borer, Hypothenemus hampei. Sci. Rep. 2015, 5, 12525. [Google Scholar] [CrossRef]
  31. Apriyanto, A. Draft genome sequence, annotation, and SSR mining data of Elaeidobius kamerunicus Faust., an essential oil palm pollinating weevil. Data Brief 2021, 34, 106745. [Google Scholar] [CrossRef] [PubMed]
  32. Harrop, T.W.; Le Lec, M.F.; Jauregui, R.; Taylor, S.E.; Inwood, S.N.; van Stijn, T.; Henry, H.; Skelly, J.; Ganesh, S.; Ashby, R.L.; et al. Genetic diversity in invasive populations of argentine stem weevil associated with adaptation to biocontrol. Insects 2020, 11, 441. [Google Scholar] [CrossRef] [PubMed]
  33. Ellegren, H. Microsatellites: Simple sequences with complex evolution. Nat. Rev. Genet. 2004, 5, 435–445. [Google Scholar] [CrossRef] [PubMed]
  34. Subramanian, S.; Mishra, R.K.; Singh, L. Genome-wide analysis of microsatellite repeats in humans: Their abundance and density in specific genomic regions. Genome Biol. 2003, 4, R13. [Google Scholar] [CrossRef]
  35. Manee, M.M.; Al-Shomrani, B.M.; Al-Fageeh, M.B. Genome-wide characterization of simple sequence repeats in Palmae genomes. Genes Genom. 2020, 42, 597–608. [Google Scholar] [CrossRef] [Green Version]
  36. Qian, J.; Xu, H.; Song, J.; Xu, J.; Zhu, Y.; Chen, S. Genome-wide analysis of simple sequence repeats in the model medicinal mushroom Ganoderma lucidum. Gene 2013, 512, 331–336. [Google Scholar] [CrossRef]
  37. Karaoglu, H.; Lee, C.M.Y.; Meyer, W. Survey of simple sequence repeats in completed fungal genomes. Mol. Biol. Evol. 2005, 22, 639–649. [Google Scholar] [CrossRef]
  38. Xu, Y.; Li, W.; Hu, Z.; Zeng, T.; Shen, Y.; Liu, S.; Zhang, X.; Li, J.; Yue, B. Genome-wide mining of perfect microsatellites and tetranucleotide orthologous microsatellites estimates in six primate species. Gene 2018, 643, 124–132. [Google Scholar] [CrossRef]
  39. Sharma, P.C.; Grover, A.; Kahl, G. Mining microsatellites in eukaryotic genomes. Trends Biotechnol. 2007, 25, 490–498. [Google Scholar] [CrossRef]
  40. Webster, M.T.; Smith, N.G.; Ellegren, H. Microsatellite evolution inferred from human–chimpanzee genomic sequence alignments. Proc. Natl. Acad. Sci. USA 2002, 99, 8748–8753. [Google Scholar] [CrossRef]
  41. Pascual, M.; Schug, M.D.; Aquadro, C.F. High density of long dinucleotide microsatellites in Drosophila subobscura. Mol. Biol. Evol. 2000, 17, 1259–1267. [Google Scholar] [CrossRef] [PubMed]
  42. Chapman, M.A. Optimizing depth and type of high-throughput sequencing data for microsatellite discovery. Appl. Plant Sci. 2019, 7, e11298. [Google Scholar] [CrossRef] [PubMed]
  43. Ding, S.; Wang, S.; He, K.; Jiang, M.; Li, F. Large-scale analysis reveals that the genome features of simple sequence repeats are generally conserved at the family level in insects. BMC Genom. 2017, 18, 848. [Google Scholar] [CrossRef]
  44. Wang, X.T.; Zhang, Y.J.; Qiao, L.; Chen, B. Comparative analyses of simple sequence repeats (SSRs) in 23 mosquito species genomes: Identification, characterization and distribution (Diptera: Culicidae). Insect Sci. 2019, 26, 607–619. [Google Scholar] [CrossRef] [PubMed]
  45. Zhao, H.; Yang, L.; Peng, Z.; Sun, H.; Yue, X.; Lou, Y.; Dong, L.; Wang, L.; Gao, Z. Developing genome-wide microsatellite markers of bamboo and their applications on molecular marker assisted taxonomy for accessions in the genus Phyllostachys. Sci. Rep. 2015, 5, 8018. [Google Scholar] [CrossRef]
  46. Wang, Q.; Fang, L.; Chen, J.; Hu, Y.; Si, Z.; Wang, S.; Chang, L.; Guo, W.; Zhang, T. Genome-wide mining, characterization and development of microsatellite markers in Gossypium species. Sci. Rep. 2015, 5, 10638. [Google Scholar] [CrossRef]
  47. Kumpatla, S.P.; Mukhopadhyay, S. Mining and survey of simple sequence repeats in expressed sequence tags of dicotyledonous species. Genome 2005, 48, 985–998. [Google Scholar] [CrossRef]
  48. Peng, X.; Yang, Z.; Xu, L.; Wang, H.; Guo, C.; Hu, P. Genome Survey Sequencing and Identification of Genomic SSR Markers for Batocera Horsfieldi (Coleoptera: Cerambycidae). 2021. Available online: https://www.researchsquare.com/article/rs-498077/v1 (accessed on 21 July 2022).
  49. Manee, M.M.; Algarni, A.T.; Alharbi, S.N.; Al-Shomrani, B.M.; Ibrahim, M.A.; Binghadir, S.A.; Al-Fageeh, M.B. Genome-wide characterization and analysis of microsatellite sequences in camelid species. Mammal Res. 2020, 65, 359–373. [Google Scholar] [CrossRef]
  50. Kim, T.S.; Booth, J.G.; Gauch, H.G.; Sun, Q.; Park, J.; Lee, Y.H.; Lee, K. Simple sequence repeats in Neurospora crassa: Distribution, polymorphism and evolutionary inference. BMC Genom. 2008, 9, 31. [Google Scholar] [CrossRef]
  51. Bacolla, A.; Larson, J.E.; Collins, J.R.; Li, J.; Milosavljevic, A.; Stenson, P.D.; Cooper, D.N.; Wells, R.D. Abundance and length of simple repeats in vertebrate genomes are determined by their structural properties. Genome Res. 2008, 18, 1545–1553. [Google Scholar] [CrossRef]
  52. Charlesworth, B.; Sniegowski, P.; Stephan, W. The evolutionary dynamics of repetitive DNA in eukaryotes. Nature 1994, 371, 215–220. [Google Scholar] [CrossRef] [PubMed]
  53. Li, S.; Zhang, G.; Li, X.; Wang, L.; Yuan, J.; Deng, C.; Gao, W. Genome-wide identification and validation of simple sequence repeats (SSRs) from Asparagus officinalis. Mol. Cell. Probes 2016, 30, 153–160. [Google Scholar] [CrossRef] [PubMed]
  54. Xiao, J.; Zhao, J.; Liu, M.; Liu, P.; Dai, L.; Zhao, Z. Genome-wide characterization of simple sequence repeat (SSR) loci in Chinese jujube and jujube SSR primer transferability. PLoS ONE 2015, 10, e0127812. [Google Scholar]
  55. Sureshkumar, S.; Todesco, M.; Schneeberger, K.; Harilal, R.; Balasubramanian, S.; Weigel, D. A genetic defect caused by a triplet repeat expansion in Arabidopsis thaliana. Science 2009, 323, 1060–1063. [Google Scholar] [CrossRef] [PubMed]
  56. Li, Y.C.; Korol, A.B.; Fahima, T.; Nevo, E. Microsatellites within genes: Structure, function, and evolution. Mol. Biol. Evol. 2004, 21, 991–1007. [Google Scholar] [CrossRef]
  57. Morgante, M.; Hanafey, M.; Powell, W. Microsatellites are preferentially associated with nonrepetitive DNA in plant genomes. Nat. Genet. 2002, 30, 194–200. [Google Scholar] [CrossRef]
  58. Qi, W.H.; Yan, C.c.; Li, W.J.; Jiang, X.M.; Li, G.Z.; Zhang, X.Y.; Hu, T.Z.; Li, J.; Yue, B.S. Distinct patterns of simple sequence repeats and GC distribution in intragenic and intergenic regions of primate genomes. Aging (Albany NY) 2016, 8, 2635. [Google Scholar] [CrossRef]
  59. Hong, C.P.; Piao, Z.Y.; Kang, T.W.; Batley, J.; Yang, T.; Hur, Y.; Bhak, J.; Park, B.; Edwards, D.; Lim, Y.P.; et al. Genomic distribution of simple sequence repeats in Brassica rapa. Mol. Cells 2007, 23, 349. [Google Scholar]
  60. Metzgar, D.; Bytof, J.; Wills, C. Selection against frameshift mutations limits microsatellite expansion in coding DNA. Genome Res. 2000, 10, 72–80. [Google Scholar]
  61. Vinogradov, A.E. DNA helix: The importance of being GC-rich. Nucleic Acids Res. 2003, 31, 1838–1844. [Google Scholar] [CrossRef]
  62. Kudla, G.; Lipinski, L.; Caffin, F.; Helwak, A.; Zylicz, M. High guanine and cytosine content increases mRNA levels in mammalian cells. PLoS Biol. 2006, 4, e180. [Google Scholar] [CrossRef]
  63. Zhang, Z.Q. Animal Biodiversity: An Outline of Higher-Level Classification and Survey of Taxonomic Richness; Magnolia Press: Auckland, New Zealand, 2011. [Google Scholar]
Ijms 23 09847 g001 550
Figure 1. Assessment of the assembled weevil genomes using the arthropoda_odb10 BUSCO dataset (1013 single-copy orthologs).
Figure 1. Assessment of the assembled weevil genomes using the arthropoda_odb10 BUSCO dataset (1013 single-copy orthologs).
Ijms 23 09847 g001
Ijms 23 09847 g002 550
Figure 2. Type composition of SSRs identified in the RPW genome assemblies without applying any search criteria. (A–C) represent adult female, adult male, and larval samples, respectively.
Figure 2. Type composition of SSRs identified in the RPW genome assemblies without applying any search criteria. (A–C) represent adult female, adult male, and larval samples, respectively.
Ijms 23 09847 g002
Ijms 23 09847 g003 550
Figure 3. Comparison of SSR types according to their percentage, frequency, density, and GC content in the weevil genome assemblies. Percentages were determined by dividing the number of SSRs of a given type by the total number for that species. (A–D) represent the percentage, frequency, density, and GC content of SSRs, respectively.
Figure 3. Comparison of SSR types according to their percentage, frequency, density, and GC content in the weevil genome assemblies. Percentages were determined by dividing the number of SSRs of a given type by the total number for that species. (A–D) represent the percentage, frequency, density, and GC content of SSRs, respectively.
Ijms 23 09847 g003
Ijms 23 09847 g004 550
Figure 4. Percentages of top SSR motif types in the weevil draft genomes. Percentages were calculated by dividing the number of SSRs of each motif type by the total number of SSRs in each assembly. (A–L) represent the weevil draft genomes R. ferrugineus (F), R. ferrugineus (M), R. ferrugineus (L), S. oryzae, H. hampei, D. ponderosae (F), D. ponderosae (M), P. strobi, E. kamerunicus, I. nitidus, L. oregonensis, and L. bonariensis.
Figure 4. Percentages of top SSR motif types in the weevil draft genomes. Percentages were calculated by dividing the number of SSRs of each motif type by the total number of SSRs in each assembly. (A–L) represent the weevil draft genomes R. ferrugineus (F), R. ferrugineus (M), R. ferrugineus (L), S. oryzae, H. hampei, D. ponderosae (F), D. ponderosae (M), P. strobi, E. kamerunicus, I. nitidus, L. oregonensis, and L. bonariensis.
Ijms 23 09847 g004
Ijms 23 09847 g005 550
Figure 5. Comparison of SSR percentage, frequency, density, and GC content among different genomic regions in four beetle draft genomes. (A–D) represent the percentage, frequency, density, and GC content of SSRs, respectively.
Figure 5. Comparison of SSR percentage, frequency, density, and GC content among different genomic regions in four beetle draft genomes. (A–D) represent the percentage, frequency, density, and GC content of SSRs, respectively.
Ijms 23 09847 g005
Ijms 23 09847 g006 550
Figure 6. Relative frequency of mono- to hexanucleotide SSRs in different genomic regions of four draft beetle genomes. (A–D) represent CDSs, exons, intronic regions, and intergenic regions, respectively.
Figure 6. Relative frequency of mono- to hexanucleotide SSRs in different genomic regions of four draft beetle genomes. (A–D) represent CDSs, exons, intronic regions, and intergenic regions, respectively.
Ijms 23 09847 g006
Ijms 23 09847 g007 550
Figure 7. The most frequent SSR motif types in different genomic regions of four draft beetle genomes. (A–D) represent CDSs, exons, intronic regions, and intergenic regions, respectively.
Figure 7. The most frequent SSR motif types in different genomic regions of four draft beetle genomes. (A–D) represent CDSs, exons, intronic regions, and intergenic regions, respectively.
Ijms 23 09847 g007
Table
Table 1. Overview of draft genomes for the nine weevil species.
Table 1. Overview of draft genomes for the nine weevil species.
Insect NameCommon NameGenome Size (Mb)Number of SSRsFrequency (SSR/Mb)Density (bp/Mb)SSRs Content (%)Reference
R. ferrugineus (F)Female red palm weevil1121.3657,17550.991320.920.13This study
R. ferrugineus (M)Male red palm weevil782.1050,72364.861445.930.14[27]
R. ferrugineus (L)Red palm weevil larva589.4067,261114.113649.450.36[28]
S. oryzaeRice weevil770.5784,391109.523287.110.33[29]
H. hampeiCoffee berry borer162.5713,09280.533260.240.33[30]
D. ponderosae (F)Female mountain pine beetle223.74650529.07481.680.05Unpublished
D. ponderosae (M)Male mountain pine beetle224.79680330.26511.440.05Unpublished
P. strobiWhite pine weevil2025.02154,51176.301516.870.15Unpublished
E. kamerunicusAfrican oil palm weevil269.64439716.31249.980.02[31]
I. nitidusQinghai spruce bark beetle345.0048372140.213127.270.31Unpublished
L. oregonensisCarrot weevil1293.28534,123412.9914,406.261.44Unpublished
L. bonariensisArgentine stem weevil1112.44156,716140.883976.450.40[32]
Table
Table 2. Number, length, frequency, and density of SSRs by type (mono- to hexanucleotide repeats) in the investigated weevil genomes.
Table 2. Number, length, frequency, and density of SSRs by type (mono- to hexanucleotide repeats) in the investigated weevil genomes.
Repeat TypeParameterR. ferrugineus (F)R. ferrugineus (M)R. ferrugineus (L)S. oryzaeH. hampeiD. ponderosae (F)D. ponderosae (M)P. strobiE. kamerunicusI. nitidusL. oregonensisL. bonariensis
Mono-Number of SSRs63935748490622,48439241415139548,34593912,459286,85957,721
Total length (bp)98,92979,57769,520429,33170,50918,62518,802725,18111,752178,8284,316,826809,526
Average length (bp)15.4713.8414.1719.0917.9713.1613.4815.0012.5214.3515.0514.02
Frequency (SSR/Mb)5.707.358.3229.1824.146.326.2123.873.4836.11221.8151.89
Density (bp/Mb)88.22101.75117.95557.16433.7183.2483.64358.1143.58518.333337.89727.70
Di-Number of SSRs34,68127,78147,993196,963077944112818,61251719,728190,06446,373
Total length (bp)1,084,748748,0021,823,9361,157,184335,41615,83419,970360,8067412493,79813,130,6002,463,404
Average length (bp)31.2826.9238.0058.75109.0016.7717.7019.3914.3425.0369.0953.12
Frequency (SSR/Mb)30.9335.5281.4325.5618.934.225.029.191.9257.18146.9641.69
Density (bp/Mb)967.35956.403094.551501.732063.2070.7788.84178.1727.491431.2810,152.942214.42
Tri-Number of SSRs85178897740715,78640282558266345,395105010,35723,75216,301
Total length (bp)142,764149,838126,207284,63179,14344,47546,9411,137,37815,882196,965457,914304,446
Average length (bp)16.7616.8417.0418.0319.6482117.3917.6325.0615.1319.0219.2818.68
Frequency (SSR/Mb)7.6011.3812.5720.4924.7811.4311.8522.423.8930.0218.3714.65
Density (bp/Mb)127.31191.58214.13369.38486.82198.78208.82561.6658.90570.90354.07273.67
Tetra-Number of SSRs53686700551519,69816271310132332,6981510391521,08529,736
Total length (bp)94,688116,57297,284422,34030,38022,68822,752621,16024,60495,912406,368654,160
Average length (bp)17.6417.4017.6421.4418.6717.3217.2018.9916.2924.5019.2721.99
Frequency (SSR/Mb)4.798.579.3625.5610.015.865.8916.155.6011.3516.3026.73
Density (bp/Mb)84.44149.05165.06548.09186.87101.40101.21306.7491.25278.00314.21588.04
Penta-Number of SSRs13181324115153753822302396630348128957865784
Total length (bp)30,55528,99525,330135,84513,17549005145147,940696070,845136,620158,590
Average length (bp)23.1821.9022.0125.2734.4921.3021.5222.3120.0054.9623.6127.42
Frequency (SSR/Mb)1.181.691.956.982.351.031.063.271.293.744.475.20
Density (bp/Mb)27.2537.0742.98176.2981.0421.9022.8973.0625.81205.35105.64142.56
Hexa-Number of SSRs89827328913525448552831336246577801
Total length (bp)29,53878728718103,61413981248135679,23079242,576183,01233,432
Average length (bp)32.8928.8430.1776.6425.8926.0024.6527.9924.0068.2327.8341.74
Frequency (SSR/Mb)0.800.350.491.750.330.210.241.400.121.815.090.72
Density (bp/Mb)26.3410.0714.79134.478.605.586.0339.132.94123.41141.5130.05
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Manee, M.M.; Al-Shomrani, B.M.; Altammami, M.A.; El-Shafie, H.A.F.; Alsayah, A.A.; Alhoshani, F.M.; Alqahtani, F.H. Microsatellite Variation in the Most Devastating Beetle Pests (Coleoptera: Curculionidae) of Agricultural and Forest Crops. Int. J. Mol. Sci. 2022, 23, 9847. https://doi.org/10.3390/ijms23179847

AMA Style

Manee MM, Al-Shomrani BM, Altammami MA, El-Shafie HAF, Alsayah AA, Alhoshani FM, Alqahtani FH. Microsatellite Variation in the Most Devastating Beetle Pests (Coleoptera: Curculionidae) of Agricultural and Forest Crops. International Journal of Molecular Sciences. 2022; 23(17):9847. https://doi.org/10.3390/ijms23179847

Chicago/Turabian Style

Manee, Manee M., Badr M. Al-Shomrani, Musaad A. Altammami, Hamadttu A. F. El-Shafie, Atheer A. Alsayah, Fahad M. Alhoshani, and Fahad H. Alqahtani. 2022. "Microsatellite Variation in the Most Devastating Beetle Pests (Coleoptera: Curculionidae) of Agricultural and Forest Crops" International Journal of Molecular Sciences 23, no. 17: 9847. https://doi.org/10.3390/ijms23179847

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop