- Research note
- Open Access
- Published:
Chromosome length genome assembly of the redbanded stink bug, Piezodorus guildinii (Westwood)
BMC Research Notes volume 15, Article number: 115 (2022)
Abstract
Objective
The redbanded stink bug (RBSB), Piezodorus guildinii (Hemiptera: Pentatomidae), is native to the Caribbean Basin and is currently considered an invasive pest in Florida, Louisiana, Mississippi, and Texas in the southern United States. Although RBSB is an economically important invasive pest in the USA, relatively few studies have been conducted to understand molecular mechanisms, population genetic structure, and the genetic basis of resistance to insecticides. The objective of this work was to obtain a high-quality genome assembly to develop genomic resources to conduct population genetic, genomic, and physiological studies of the RBSB.
Results
The genome of RBSB was sequenced with Pacific Biosciences technology followed by two rounds of scaffolding using Chicago libraries and HiC proximity ligation to obtain a high-quality assembly. The genome assembly contained 800 scaffolds larger than 1 kbp and the N50 was 170.84 Mbp. The largest scaffold was 222.22 Mbp and 90% of the genome was included in the 7 scaffolds larger than 118 Mbp. The number of megabase scaffolds also matched the number of chromosomes in this insect. The genome sequence will facilitate the development of resources to conduct studies on genetics, transcriptomics, and physiology of RBSB.
Introduction
The redbanded stink bug (RBSB), Piezodorus guildinii (Westwood) (Hemiptera: Pentatomidae), is native to the Caribbean Basin and is currently considered an invasive pest of soybeans and several other commercially grown crops in Florida, Louisiana, Mississippi, and Texas in the southern United States [1,2,3]. Uncontrolled outbreaks of RBSB can cause significant economic damage to soybeans from early seed development stages to mature seeds [4]. Although RBSB is an economically important invasive pest in the USA, relatively few studies have been conducted to understand genetic, population genetic structure, and genetic basis of resistance to insecticides. Biology, ecology, host plants, and pest status of this insect has been previously studied [2, 5,6,7,8,9,10,11,12]. Resistance to insecticides in RBSB have been documented [13], but the genetic basis of the insecticide resistance in this insect is not well understood. So far, only one population genetic study has been carried out on this species using 1,337 SNP markers that identified the presence of genetic structure separating populations in USA and Brazil [14]. In order to develop genetic resources to conduct functional genomic, population genetic, and physiological studies, we sequenced the genome of RBSB (Additional file 1: Fig. S1) with Pacific Biosciences long read technology using Chicago libraries and then we assembled the draft input assembly that was used for scaffolding Illumina short reads from HiC proximity ligation libraries to obtain a high-quality assembly (Fig. 1, Additional file 1: Figs. S2 and S3).
HiC linkage density histogram of Piezodorus guildinii genome assembly. The x and y axes in this histogram show the mapping positions of the first and second read in a read pair, respectively, grouped into bins. The intensity of color of each square represents the number of read pairs within that bin. Scaffolds less than 1 Mb were excluded from this histogram
Main text
Methods
Genomic DNA libraries for PacBio sequencing and initial assemblies were prepared by Dovetail Genomics (Scottes Valley, CA, USA). A Qiagen tissue kit and a tip-20 mini column (Qiagen, Germantown, MD) was used to isolate high molecular weight DNA from a field collected female RBSB. DNA was quantified using Qubit 2.0 fluorometer (Life Technologies, Carlsbad, CA) and a PacBio SMRTbell library with approximately 20 kbp was constructed using SMRTbell Express Template Prep Kit 2.0 (PacBio, Menlo Park, CA) following the manufacturer’s protocol. DNA sequencing was performed on a PacBio Sequel II sequencer using Sequel II 8 M SMRT cells generating 108 Gb of continuous long reads (CLR). Initial assembly was performed using the Wtdbg2 v.2.5. assembler [15] with the following parameters: --genome_size 1.0 g --read_type sq --min_read_len 20,000 --min_aln_len 8192. Blob Tools v1.1.1 [16] with default parameters ([-t HITS…] [-x TAXRULE…] [-m 0.0] [-d 0.0] [–tax_collision_random]) was used to identify potential contamination based on BLAST (v2.9) results of the assembly against the NT database. The assembly was filtered to remove potential haplotypic duplications using Purge Dups v1.1.2 [17] (parameters: -2 -T cutoffs) to obtain a purged draft assembly used for scaffolding. Scaffolding was performed using HiRise v2.1.5 pipeline [18] with default settings. Proximity ligation libraries were prepared using Dovetail Omni-C library protocol by digesting formaldehyde fixed chromatin with a DNAse I, repairing chromatin ends, and biotinylated bridge adapter ligation followed by proximity ligation of adapter containing ends. Then crosslinks were reversed, DNA was purified and was treated to remove biotin that was not internal to ligated fragments. Biotin-labelled DNA fragments were isolated using streptavidin beads and enriched by PCR. The library was sequenced on an Illumina HiSeqX platform to produce approximately 30 × sequence coverage. The draft de novo assembly from PacBio reads and Dovetail Omni-C proximity ligated library reads were used as input data for HiRise v2.1.5 pipeline. Dovetail Omni-C library sequences were aligned to draft input assembly using BWA-mem v0.7.17-r1188 [19, 20] using parameters 5SP -T0. The separations of Dovetail OmniC read pairs mapped within draft scaffolds were analyzed by HiRise to produce a likelihood model for genomic distance between read pairs, and the model was used to identify and break putative mis-joins, to score prospective joins, and make joins above a threshold.
The genome completeness evaluation was performed with BUSCO version 5.2.2 [21] using the Arthropoda (arthropoda_odb10.2019-11-20 1013) and Hemiptera (hemiptera_odb10.2019-11-20 2510) databases with 1,013 and 2,510 BUSCOs, respectively. The repeat annotation was done with RepeatModeler using the Dfam TE tools docker container version 1.4 (https://github.com/Dfam-consortium/TETools). The repeat classification was performed with RepeatClassifier Version 2.0.2 and RepeatMasker version 4.1.2-p1 [22] to identify the types of repeats in the RBSB genome. RepeatMasker was run in sensitive mode with rmblastn version 2.11.0 + . The database of repeats used for classification was Dfam 3.4 [23].
To identify reads derived from mitochondrial DNA (mtDNA), PacBio sequence reads from the shotgun library were mapped to the Nezara viridula (L.) mtDNA genome (Accession: EF208087.1) with relaxed parameters (length fraction 0.3 and similarity fraction 0.75). Mapped sequence reads were extracted, and de novo assembled at 85% similarity and 75% length fraction. A consensus of a 20,764 nt RBSB mtDNA contig with 84% identity to the stink bug Eurydema ventralis (Kolenati) mtDNA (Accession:MG584837.1) was selected for use as the reference to map the PacBio reads with a higher stringency at 85% similarity fraction and 75% length fraction. A total of 2,099 out of 9,113,332 reads were mapped to the RBSB mtDNA reference with > 150-fold sequence coverage across the coding regions and > 75-fold sequence coverage in the control region. Consensus of this contig was used to obtain 18,889 bp mitochondrial genome of RBSB.
Results and discussion
Published and unpublished genome assemblies of other pentatomids available in databases indicate that the RBSB genome assembly presented here is comparable to chromosome length assembly of Aelia acuminata (N50 = 172.2 Mbp, largest scaffold = 235.2 Mbp; [24]) and superior to the assemblies of Euschistus heros (N50 = 2.46 Mbp; PRJNA489772), and Halyomorpha halys (N50 = 802 Kbp; [25]). The assembled size of the RBSB genome was 1.205 Gbp with an N50 of 170.835 Mbp, N90 of 118.462 Mbp, L50 of 4, and L90 of 7. The final genome assembly contained 800 scaffolds larger than 1 Kbp and the smallest and the largest scaffold were 118.462 and 222.218 Mbp, respectively (Table 1). Karyotyping of RBSB identified six pairs of autosomes and a pair of sex chromosomes designated X and Y (2n = 14) [26]. The seven largest scaffolds in the RBSB genome assembly matched the haploid chromosome number in RBSB and other pentatomid bugs.
The genome assembly of RBSB was also found to be highly complete for single-copy markers conserved within the Arthropoda and Hemiptera clades with 96.5% and 96.2% completeness values (Table 2). The low duplication and fragmentation coupled with the completeness underlines the integrity of this genome assembly.
We performed two rounds of repeat annotations for the RBSB genome. It was analyzed for known repeat families in Insecta present in the DFAM 2.4 [27] (Additional file 1: Table S1). LINE (11.64%) were the predominant retroelements found followed by SINEs (1.21%). The Tc1-IS630-Pogo transposon was the predominant DNA transposon repeat family (4.51%). A de novo repeat annotation using RepeatModeler [28] identified 2338 RepeatScout/RECON families and 181 LTR repeat families. All annotations are available to the research community at the AgriVectors portal [29].
The complete sequence of the mitochondrial genome, which is often used in population genetics and molecular identification of insects is not currently available for this species. Availability of the complete mitochondrial genome assembled in this study (BioSample SAMN23701154) will provide additional DNA markers to conduct population studies that require relatively highly variable mitochondrial genes such as NADH dehydrogenases (ND) and cytochrome B.
Additionally, having an official gene set predicted by the NCBI eukaryotic annotation pipeline (https://www.ncbi.nlm.nih.gov/genome/annotation_euk/process/) based on this high-quality genome will aid in conducting expression profiling experiments intended to elucidate physiological responses to various host plant species and insecticides. Large scale insect genomics projects like the i5k [30] and more recently Ag100Pest [31] have also highlighted the long-term benefits of building open access databases and genomics resources for the community.
A limited number of molecular resources are currently available for RBSB in public databases, which include 107 microsatellite sequences, six partial cytochrome oxidase I (COI) subunit sequences [32], a 3.18 Mbp partial assembly with 1932 genomic contigs (N50 = 1494 bp; PRJNA263369), and 17 Gbp of genotype by sequencing (GBS) data in the NCBI Sequence Read Archive [14]. Transcriptome, genome, or proteome data are not available for this economically important pest species. Besides the high-quality genome assembly, we have generated triplicate RNASeq libraries for all life stages from eggs, first instar to fifth instar nymphs, and adult males and females that will be sequenced to obtain a minimum coverage of 25 million (2 × 150 paired end) reads for each library. Oxford Nanopore long reads will also be obtained by pooling mRNA from all life stages of RBSB. Long Oxford Nanopore reads will facilitate annotation of full-length gene models in the genome as well as qualitative identification of isoforms from each life stage. The NCBI structural annotation will be followed by functional annotation to identify gene ontology terms and pathways[33] which will be made available on the AgriVectors portal [29]. In addition, we plan to sequence the RBSB genomic DNA with Oxford Nanopore reads to identify epigenetic modifications such as methylation in the genome [34, 35].
Having a well characterized high quality genome assembly will provide a more robust foundation for developing genetic markers for population genetic studies, linkage mapping, and identifying genomic regions associated with regulation of gene expression, host selection, and insecticide resistance by comparative genomic studies using published genomes and transcriptomes of other insects and pentatomid stink bugs [24, 25, 36,37,38].
Limitations
Genomic DNA library preparation and assembly were performed with proprietary methods developed by a service provider with PacBio and Illumina sequencing versions currently available. Library construction and sequencing methods, assembly software, and other data used for comparative analysis may be updated in the future. PacBio continuous long reads (CLR) may contain insertion and deletion errors, some of which may have escaped correction during the assembly process.
Availability of data and materials
All raw sequencing data and assemblies have been submitted to NCBI BioProject PRJNA686660.
References
Bundy CS, Esquivel JF, Panizzi AR, Eger JE, Davis JA, Jones WA. Piezodorus guildinii (Westwood). In: McPherson JE, editor. Invasive stink bugs and related species (Pentatomoidea). Boca Raton: CRC Press; 2018. p. 425–52.
Panizzi A, Slansky F Jr. Legume host impact on performance of adult Piezodorus guildinii (Westwood)(Hemiptera: Pentatomidae). Environ Entomol. 1985;14:237–42.
Panizzi A, Slansky F Jr. New host plant records for the stink bug Piezodorus guildinii in Florida (Hemiptera: Pentatomidae). Fla Entomol. 1985;68:215–6.
Correaferreira B, Moscardi F. Seasonal occurrence and host spectrum of egg parasitoids associated with soybean stink bugs. Biol Control. 1995;5:196–202.
Panizzi AR, Smith JG. Biology of Piezodorus guildinii: oviposition, development time, adult sex ratio, and longevity. Ann Entomol Soc Am. 1977;70:35–9.
Panizzi A. Performance of Piezodorus guildinii on four species of Indigofera legumes. Entomol Exp Appl. 1992;63:221–8.
Link D, Panichi J, Concatto L. Oviposition by Piezodorus guildinii (Westwood, 1837) on bean plants. Rev Cent Cienc Rur. 1980;10:271–6.
Cividanes FJ, Parra JR. Ecological zoning of Nezara viridula (L.), Piezodorus guildinii (West.) and Euschistus heros (Fabr.)(Heteroptera: Pentatomidae) in four soyabean-producing states of Brazil. Anais Soc Entomol Brasil. 1994;23:219–26.
Oliveira ÉD, Panizzi AR. Performance of nymphs and adults of Piezodorus guildinii (Westwood)(Hemiptera: Pentatomidae) on soybean pods at different developmental stages. Braz Arch Biol Tech. 2003;46:187–92.
Bastola A. Ecology and Biology of the Redbanded Stink Bug, Piezodorus guildinii (Westwood) in Louisiana [Dissertation]. Baton Rouge: Louisiana State University; 2017.
Zerbino M, Miguel L, Altier N, Panizzi A. Overwintering of Piezodorus guildinii (Heteroptera, Pentatomidae) populations. Neotrop Entomol. 2020;49:179–90.
Gomez VA, Gaona EF, Arias OR, De Lopez MB, Ocampos OE. Biological aspects of Piezodorus guildinii (Westwood)(Hemiptera: Pentatomidae) reared in laboratory with different insects diets. Rev Soc Entomol Argentina. 2013;72:27–34.
Baur M, Sosa-Gomez D, Ottea J, Leonard B, Corso I, Da Silva J, et al. Susceptibility to insecticides used for control of Piezodorus guildinii (Heteroptera: Pentatomidae) in the United States and Brazil. J Econ Entomol. 2010;103:869–76.
Zucchi MI, Cordeiro EM, Allen C, Novello M, Viana JPG, Brown PJ, et al. Patterns of genome-wide variation, population differentiation and SNP discovery of the red banded stink bug (Piezodorus guildinii). Sci Rep. 2019;9:1–11.
Ruan J, Li H. Fast and accurate long-read assembly with wtdbg2. Nat Methods. 2020;17:155–8.
Laetsch DR, Blaxter ML. BlobTools: Interrogation of genome assemblies. F1000Res. 2017;6:1287.
Guan D, McCarthy SA, Wood J, Howe K, Wang Y, Durbin R. Identifying and removing haplotypic duplication in primary genome assemblies. Bioinformatics. 2020;36:2896–8.
Putnam NH, O’Connell BL, Stites JC, Rice BJ, Blanchette M, Calef R, et al. Chromosome-scale shotgun assembly using an in vitro method for long-range linkage. Genome Res. 2016;26:342–50.
Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25:1754–60.
Lieberman-Aiden E, Van Berkum NL, Williams L, Imakaev M, Ragoczy T, Telling A, et al. Comprehensive mapping of long-range interactions reveals folding principles of the human genome. Science. 2009;326:289–93.
Manni M, Berkeley MR, Seppey M, Simao FA, Zdobnov EM. BUSCO update: novel and streamlined workflows along with broader and deeper phylogenetic coverage for scoring of eukaryotic, prokaryotic, and viral genomes. arXiv preprint arXiv:210611799. 2021.
Smit AFA, Hubley R, Green P. RepeatMasker Open-4.0. 2015. http://www.repeatmasker.org Accessed 15 Aug 2021.
Storer J, Hubley R, Rosen J, Wheeler TJ, Smit AF. The Dfam community resource of transposable element families, sequence models, and genome annotations. Mob DNA. 2021;12:2.
Crowley L, Barclay MVL. University of Oxford and Wytham Woods Genome Acquisition Lab et al. The genome sequence of the bishop’s mitre shieldbug, Aelia acuminata (Linnaeus, 1758). Wellcome Open Res. 2021;6:320.
Sparks ME, Rhoades JH, Nelson DR, Kuhar D, Lancaster J, Lehner B, et al. A Transcriptome Survey Spanning Life Stages and Sexes of the Harlequin Bug Murgantia histrionica. Insects. 2017. https://doi.org/10.3390/insects8020055.
Rebagliati PJ, Mola LM, Papeschi AG. Karyotype and meiotic behaviour of the holokinetic chromosomes of six Argentine species of Pentatomidae (Heteroptera). Caryologia. 2001;54:339–47.
Hubley R, Finn RD, Clements J, Eddy SR, Jones TA, Bao W, et al. The Dfam database of repetitive DNA families. Nucleic Acids Res. 2016;44:D81–9.
Flynn JM, Hubley R, Goubert C, Rosen J, Clark AG, Feschotte C, et al. RepeatModeler2 for automated genomic discovery of transposable element families. Proc Nat Acad Sci USA. 2020;117:9451–7.
Saha S, Cooper WR, Hunter W, Mueller L, Consortium A. an open access resource portal for arthropod vectors and agricultural pathosystems: AgriVectors. org. The 1st International Electronic Conference on Entomology 2021; https://doi.org/10.3390/IECE-10576.
i5K-Consortium. The i5K Initiative: advancing arthropod genomics for knowledge, human health, agriculture, and the environment. J Hered. 2013;104:595–600.
Childers AK, Geib SM, Sim SB, Poelchau MF, Coates BS, Simmonds TJ, et al. The USDA-ARS Ag100Pest initiative: high-quality genome assemblies for agricultural pest arthropod research. Insects. 2021;12:626.
Greenstone M, Tillman P, Hu J. Predation of the newly invasive pest Megacopta cribraria (Hemiptera: Plataspidae) in soybean habitats adjacent to cotton by a complex of predators. J Econ Entomol. 2014;107:947–54.
Saha S, Cooksey AM, Childers AK, Poelchau MF, McCarthy FM. Workflows for rapid functional annotation of diverse arthropod genomes. BioRxiv. 2021. https://doi.org/10.1101/2021.06.12.448177.
Rand AC, Jain M, Eizenga JM, Musselman-Brown A, Olsen HE, Akeson M, et al. Mapping DNA methylation with high-throughput nanopore sequencing. Nat Methods. 2017;14:411–3.
Yuen ZWS, Srivastava A, Daniel R, McNevin D, Jack C, Eyras E. Systematic benchmarking of tools for CpG methylation detection from nanopore sequencing. Nat Commun. 2021;12:3438.
Sparks ME, Bansal R, Benoit JB, Blackburn MB, Chao H, Chen M, et al. Brown marmorated stink bug, Halyomorpha halys (Stål), genome: putative underpinnings of polyphagy, insecticide resistance potential and biology of a top worldwide pest. BMC Genomics. 2020;21:1–26.
Sparks ME, Nelson DR, Haber AI, Weber DC, Harrison RL. Transcriptome sequencing of the striped cucumber beetle, Acalymma vittatum (F.), reveals numerous sex-specific transcripts and xenobiotic detoxification genes. BioTech. 2020;9:21.
Cao C, Sun L, Wen R, Shang Q, Ma L, Wang Z. Characterization of the transcriptome of the Asian gypsy moth Lymantria dispar identifies numerous transcripts associated with insecticide resistance. Pestic Biochem Physiol. 2015;119:54–61.
Acknowledgements
We thank Calvin Pierce (USDA ARS SIMRU) for assistance with RBSB colony maintenance and sample preparation and Dr. Nathan Little for critically reading an earlier version of this manuscript. The use or mention of a trademark or proprietary product does not constitute an endorsement, guarantee, or warranty of the product and does not imply its approval to the exclusion of other suitable products by the U.S. Department of Agriculture, an equal opportunity employer.
Funding
This research was funded by USDA Agricultural Research Service in house research project 6066-22000-091-00D of Southern Insect Management Research Unit. Partial funding was provided by the ARS/State Potato Partnership Program.
Author information
Authors and Affiliations
Contributions
Conceptualization, OPP; methodology, KCA, OPP, SS; formal analysis, KCA, OPP, SS; investigation, KCA, OPP, SS, resources, OPP, GVPR, LAM; data curation, OPP, SS; writing-original draft preparation, OPP, SS; writing-review and editing, OPP, SS, KCA, LAM and GVPR; project administration, OPP. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
Authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1:
Figure S1. Female (top left) and male (top right) of Piezodorus guildinii. Figure S2. Comparison of the cumulative lengths of the purged input assembly and the final HiRise scaffolds. Figure S3. Cumulative insert size distribution of HiC paired-end reads mapped within a chromosome. Table S1. The number of repeat units, the total length and the percentage of different repeat families identified in Piezodorus guildinii genome.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Saha, S., Allen, K.C., Mueller, L.A. et al. Chromosome length genome assembly of the redbanded stink bug, Piezodorus guildinii (Westwood). BMC Res Notes 15, 115 (2022). https://doi.org/10.1186/s13104-022-05924-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13104-022-05924-5
Keywords
- Redbanded stink bug
- Piezodorus guildinii
- Genome assembly
- Pentatomidae
- Soybean pest