Transcriptome analysis of the endangered dung beetle Copris tripartitus (Coleoptera: Scarabaeidae) and characterization of genes associated to immunity, growth, and reproduction.

BACKGROUND: Dung beetles recycle organic matter through the decomposition of feces and support ecological balance. However, these insects are threatened by the indiscriminate use of agrochemicals and habitat destruction. Copris tripartitus Waterhouse (Coleoptera: Scarabaeidae), a dung beetle, is listed as a class-II Korean endangered species. Although the genetic diversity of C. tripartitus populations has been investigated through analysis of mitochondrial genes, genomic resources for this species remain limited. In this study, we analyzed the transcriptome of C. tripartitus to elucidate functions related to growth, immunity and reproduction for the purpose of informed conservation planning. RESULTS: The transcriptome of C. tripartitus was generated using next-generation Illumina sequencing and assembled de novo using a Trinity-based platform. In total, 98.59% of the raw sequence reads were processed as clean reads. These reads were assembled into 151,177 contigs, 101,352 transcripts, and 25,106 unigenes. A total of 23,450 unigenes (93.40%) were annotated to at least one database. The largest proportion of unigenes (92.76%) were annotated to the locally curated PANM-DB. A maximum of 5,512 unigenes had homologous sequences in Tribolium castaneum. Gene Ontology (GO) analysis revealed a maximum of 5,174 unigenes in the Molecular function category. Further, in Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, a total of 462 enzymes were associated with established biological pathways. Based on sequence homology to known proteins in PANM-DB, representative immunity, growth, and reproduction-related genes were screened. Potential immunity-related genes were categorized into pattern recognition receptors (PRRs), the Toll-like receptor signaling pathway, the MyD88- dependent pathway, endogenous ligands, immune effectors, antimicrobial peptides, apoptosis, and adaptation-related transcripts. Among PRRs, we conducted detailed in silico characterization of TLR-2, CTL, and PGRP_SC2-like. Repetitive elements such as long terminal repeats, short interspersed nuclear elements, long interspersed nuclear elements and DNA elements were enriched in the unigene sequences. A total of 1,493 SSRs were identified among all unigenes of C. tripartitus. CONCLUSIONS: This study provides a comprehensive resource for analysis of the genomic topography of the beetle C. tripartitus. The data presented here clarify the fitness phenotypes of this species in the wild and provide insight to support informed conservation planning.

Transcriptome Characterization for Non-Model Endangered Lycaenids, Protantigius superans and Spindasis takanosis, Using Illumina HiSeq 2500 Sequencing.

The Lycaenidae butterflies, Protantigius superans and Spindasis takanosis, are endangered insects in Korea known for their symbiotic association with ants. However, necessary genomic and transcriptomics data are lacking in these species, limiting conservation efforts. In this study, the P. superans and S. takanosis transcriptomes were deciphered using Illumina HiSeq 2500 sequencing. The P. superans and S. takanosis transcriptome data included a total of 254,340,693 and 245,110,582 clean reads assembled into 159,074 and 170,449 contigs and 107,950 and 121,140 unigenes, respectively. BLASTX hits (E-value of 1.0 × 10(-5)) against the known protein databases annotated a total of 46,754 and 51,908 transcripts for P. superans and S. takanosis. Approximately 41.25% and 38.68% of the unigenes for P. superans and S. takanosis found homologous sequences in Protostome DB (PANM-DB). BLAST2GO analysis confirmed 18,611 unigenes representing Gene Ontology (GO) terms and a total of 5259 unigenes assigned to 116 pathways for P. superans. For S. takanosis, a total of 6697 unigenes were assigned to 119 pathways using the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database. Additionally, 382,164 and 390,516 Simple Sequence Repeats (SSRs) were compiled from the unigenes of P. superans and S. takanosis, respectively. This is the first report to record new genes and their utilization for conservation of lycaenid species population and as a reference information for closely related species.

Molecular and immunohistochemical characterization of the chitinase gene from Pieris rapae granulovirus.

The chitinase gene of baculoviruses is expressed in the late phase of virus replication in insects and possesses high exo- and endochitinase activity, which can hydrolyze chitin in the body of the insect, thus promoting terminal host liquefaction. Alphabaculovirus viral chitinases (vChitA) have been well analyzed, but information regarding viral chitinases from betabaculoviruses is limited. Whole-genome sequencing of a Korean isolate of Pieris rapae GV (PiraGV-K) predicted a putative chitinase gene corresponding to ORF10. The PiraGV-K chitinase gene had a coding sequence of 1,761 bp, encoding a protein of 586 amino acid (aa) residues, including an 18-aa putative signal peptide. Time course induction pattern observed by SDS-PAGE and subsequent Western blot with anti-PiraGV-K chitinase antibody revealed the cleavage of the signal peptide from the intact chitinase. Edman sequencing analysis was further conducted to confirm the exact nature of the mature chitinase, and the N-terminal amino acid sequence (KPGAP) exactly matched the sequence following the signal peptide sequence. The transcriptomics of PiraGV-K chitinase in infected P. rapae larvae, examined by real-time PCR, revealed a significant 75-fold increase after four days of feeding with PiraGV-K-treated leaves, with a subsequent decline at the later stages of infection. Confocal microscopic analysis showed that PiraGV-K chitinase possibly exists as a secreted protein, with strong chitinase-specific signals in fat body cells and integument at four days postinfection. Furthermore, immunogold labeling and electron microscopy studies localized the PiraGV-K chitinase in the cytoplasm and sparsely within vacuolar structures in the fat body apart from the extensive aggregation in the cuticular lining of the integument.

Molecular and immunohistochemical characterization of granulin gene encoded in Pieris rapae granulovirus genome.

Pieris rapae granulovirus (PiraGV) is highly pathogenic to the cabbage butterfly (P. rapae), an important pest of cultivated cabbages and mustard crops. It therefore holds significant promise towards exploitation as a potent bio-control agent in the field controlling the pest population. Whole-genome elucidation of the Korean isolate of the granulovirus (PiraGV-K), reported the presence of a granulin gene corresponding to ORF 1 in its genome. Comprehensive studies towards functional characterization of the gene, established that it is composed of 744 nucleotides and encodes a peptide of 247 amino acid residues. It possessed significant homology with AoGV and ClanGV with 87% identity at amino acid level. Multiple alignment data suggests that the C-terminus region of the gene had three different conserved regions. Time-course studies conducted in PiraGV-K infected P. rapae larvae revealed a significant upsurge of the transcript (134-fold) at 4 days post infection followed by a significant decline at the most advanced stages of infection. Anti-PiraGV-K granulin antibody was produced and western blot conducted with the infected larvae further confirmed the induction pattern with a protein of 30 kDa. Immunofluorescent staining showed a granulin-specific signal in fat body and integument of the infected larvae. Granulin-specific signals were noticed 2 days post infection with the eventual systemic spread of infection to the associated tracheal matrix witnessed at 4 days post infection. Immunogold labeling and electron microscopic studies further proved the cytopathological effects as the presence of numerous membrane-bound vesicles with nucleocapsids and abruption of intercellular junctions in fat body and hypertrophied cells in the integument.

Reconstructing intraordinal relationships in Lepidoptera using mitochondrial genome data with the description of two newly sequenced lycaenids, Spindasis takanonis and Protantigius superans (Lepidoptera: Lycaenidae).

Lepidoptera is one of the largest insect orders, but the phylogenetic relationships within this order, have yet to be adequately described. Among these unresolved relationships include those regarding the monophyly of the Macrolepidoptera and interfamilial relationships of the true butterflies superfamily Papilionoidea. We present two new mitochondrial genomes (mitogenomes) belonging to the butterfly family Lycaenidae to explore the phylogenetic relationships existing among lepidopteran superfamilies and true butterfly families from a mitogenome perspective, and to evaluate the characteristics of the lepidopteran mitogenomes. Our consensus phylogeny of the Lepidoptera largely supported the superfamilial relationships (((((Bombycoidea+Geometroidea)+Noctuoidea)+Pyraloidea)+Papilionoidea)+Tortricoidea), signifying a lack of support for a traditionally defined Macrolepidoptera. The familial relationships of the true butterflies concordantly recovered the previously proposed phylogenetic hypothesis (((Lycaenidae+Nymphalidae)+Pieridae)+Papilionidae). The test for the effect of optimization schemes (exclusion and inclusion of third codon position of PCGs and two rRNA genes, with and without partitions) on the resolution and relationships within the Lepidoptera have demonstrated that the majority of analyses did not substantially alter the relevant topology and node support, possibly as the result of relatively strong signal in mitogenomes for intraordinal relationships in Lepidoptera.

Mitochondrial DNA sequence variation of the swallowtail butterfly, Papilio xuthus, and the cabbage butterfly, Pieris rapae.

We analyzed a portion of mitochondrial COI gene sequences (658 bp) to investigate the genetic diversity and geographic variation of the swallowtail butterfly, Papilio xuthus L. (Lepidoptera: Papilionidae), and the cabbage butterfly, Pieris rapae L. (Lepidoptera: Pieridae). Papilio xuthus showed a moderate level of sequence divergence (0.91% at maximum) in 15 haplotypes, whereas Pi. rapae showed a moderate to high level of sequence divergence (1.67% at maximum) in 30 haplotypes, compared with other relevant studies. Analyses of population genetic structure showed that most populations are not genetically differentiated in both species. The distribution pattern of both species appears to be consistent with category IV of the phylogeographic pattern sensu Avise: a phylogenetic continuity, an absence of regional isolation of mtDNA clones, and extensive distribution of close clones. The observed pattern of genetic diversity and geographic variation of the two butterfly species seem to reflect the abundant habitats, abundant host plants, and flying abilities in connection with the lack of historical biogeographic barriers.

Transcriptome Profile of the Asian Giant Hornet (Vespa mandarinia) Using Illumina HiSeq 4000 Sequencing: De Novo Assembly, Functional Annotation, and Discovery of SSR Markers.

Vespa mandarinia found in the forests of East Asia, including Korea, occupies the highest rank in the arthropod food web within its geographical range. It serves as a source of nutrition in the form of Vespa amino acid mixture and is listed as a threatened species, although no conservation measures have been implemented. Here, we performed de novo assembly of the V. mandarinia transcriptome by Illumina HiSeq 4000 sequencing. Over 60 million raw reads and 59,184,811 clean reads were obtained. After assembly, a total of 66,837 unigenes were clustered, 40,887, 44,455, and 22,390 of which showed homologous matches against the PANM, Unigene, and KOG databases, respectively. A total of 15,675 unigenes were assigned to Gene Ontology terms, and 5,132 unigenes were mapped to 115 KEGG pathways. The zinc finger domain (C2H2-like), serine/threonine/dual specificity protein kinase domain, and RNA recognition motif domain were among the top InterProScan domains predicted for V. mandarinia sequences. Among the unigenes, we identified 534,922 cDNA simple sequence repeats as potential markers. This is the first transcriptomic analysis of the wasp V. mandarinia using Illumina HiSeq 4000. The obtained datasets should promote the search for new genes to understand the physiological attributes of this wasp.

The complete mitochondrial genome of the mountainous duskywing, Erynnis montanus (Lepidoptera: Hesperiidae): a new gene arrangement in Lepidoptera.

The mountainous duskywing, Erynnis montanus, belongs to a lepidopteran family Hesperiidae. The 15,530-bp long complete mitochondrial genome (mitogenome) of the species has the typical gene content of animals (13 protein-coding genes, two rRNA genes, 22 tRNA genes and one major non-coding A+T-rich region). As typical in lepidopteran mitogenome E. montanus mitogenome also contained a high A/T content in the whole genome (81.7%) and the CGA (arginine) as the start codon for the COI gene. Unlike other lepidopteran species, including two sequenced skippers, the E. montanus mitogenome has a unique arrangement tRNA(Ser)-tRNA(Asn), instead of the tRNA(Asn)-tRNA(Ser) found unanimously in other lepidopteran species, providing a new gene arrangement in Lepidoptera. Such rearrangement probably was likely caused by duplication of gene block tRNA(Ser)-tRNA(Asn) and subsequent random loss of tRNA(Asn) in the first copy and tRNA(Ser) in the second copy, resulting in the arrangement tRNA(Ser)-tRNA(Asn).

Analysis of the genome of a Korean isolate of the Pieris rapae granulovirus enabled by its separation from total host genomic DNA by pulse-field electrophoresis.

BACKGROUND: Most traditional genome sequencing projects involving viruses include the culture and purification of the virus particles. However, purification of virions may yield insufficient material for traditional sequencing. The electrophoretic method described here provides a strategy whereby the genomic DNA of the Korean isolate of Pieris rapae granulovirus (PiraGV-K) could be recovered in sufficient amounts for sequencing by purifying it directly from total host DNA by pulse-field gel electrophoresis (PFGE). METHODOLOGY/PRINCIPAL FINDINGS: The total genomic DNA of infected P. rapae was embedded in agarose plugs, treated with restriction nuclease and methylase, and then PFGE was used to separate PiraGV-K DNA from the DNA of P. rapae, followed by mapping of fosmid clones of the purified viral DNA. The double-stranded circular genome of PiraGV-K was found to encode 120 open reading frames (ORFs), which covered 92% of the sequence. BLAST and ORF arrangement showed the presence of 78 homologs to other genes in the database. The mean overall amino acid identity of PiraGV-K ORFs was highest with the Chinese isolate of PiraGV (~99%), followed up with Choristoneura occidentalis ORFs at 58%. PiraGV-K ORFs were grouped, according to function, into 10 genes involved in transcription, 11 involved in replication, 25 structural protein genes, and 15 auxiliary genes. Genes for Chitinase (ORF 10) and cathepsin (ORF 11), involved in the liquefaction of the host, were found in the genome. CONCLUSIONS/SIGNIFICANCE: The recovery of PiraGV-K DNA genome by pulse-field electrophoretic separation from host genomic DNA had several advantages, compared with its isolation from particles harvested as virions or inclusions from the P. rapae host. We have sequenced and analyzed the 108,658 bp PiraGV-K genome purified by the electrophoretic method. The method appears to be generally applicable to the analysis of genomes of large viruses.

Complete mitochondrial genome of the nerippe fritillary butterfly, Argynnis nerippe (Lepidoptera: Nymphalidae).

The complete mitochondrial genome sequence of the nerippe fritillary butterfly, Argynnis nerippe, which is listed as an endangered species in Korea, is described with an emphasis on the A+T-rich region. The 15,140-bp long circular molecule consisted of 13 protein-coding genes, two rRNA genes, 22 tRNA genes and 1 control region, known in insect as the A+T-rich region, as found in typical metazoans. The 329-bp long A+T-rich region located between srRNA and tRNA(Met) possessed the highest A/T content (95.7%) than any other region of the genome. Along with the several conserved sequences found typically in the lepidopteran insects the genome contained one tRNA(Met)-like and tRNA(Leu)(UUR)-like sequence in the A+T-rich region.

Phylogenetic relationships of true butterflies (Lepidoptera: Papilionoidea) inferred from COI, 16S rRNA and EF-1α sequences.

The molecular phylogenetic relationships among true butterfly families (superfamily Papilionoidea) have been a matter of substantial controversy; this debate has led to several competing hypotheses. Two of the most compelling of those hypotheses involve the relationships of (Nymphalidae + Lycaenidae) + (Pieridae + Papilionidae) and (((Nymphalidae + Lycaenidae) + Pieridae) + Papilionidae). In this study, approximately 3,500 nucleotide sequences from cytochrome oxidase subunit I (COI), 16S ribosomal RNA (16S rRNA), and elongation factor-1 alpha (EF-1α) were sequenced from 83 species belonging to four true butterfly families, along with those of three outgroup species belonging to three lepidopteran superfamilies. These sequences were subjected to phylogenetic reconstruction via Bayesian Inference (BI), Maximum Likelihood (ML), and Maximum Parsimony (MP) algorithms. The monophyletic Pieridae and monophyletic Papilionidae evidenced good recovery in all analyses, but in some analyses, the monophylies of the Lycaenidae and Nymphalidae were hampered by the inclusion of single species of the lycaenid subfamily Miletinae and the nymphalid subfamily Danainae. Excluding those singletons, all phylogenetic analyses among the four true butterfly families clearly identified the Nymphalidae as the sister to the Lycaenidae and identified this group as a sister to the Pieridae, with the Papilionidae identified as the most basal linage to the true butterfly, thus supporting the hypothesis: (Papilionidae + (Pieridae + (Nymphalidae + Lycaenidae))).

Complete mitogenome sequence of the jewel beetle, Chrysochroa fulgidissima (Coleoptera: Buprestidae).

In this study, we determined the complete mitochondrial genome of the jewel beetle, Chrysochroa fulgidissima (Coleoptera: Buprestidae), from four overlapping fragments. The 15,592 bp C. fulgidissima mitogenome exhibits a gene arrangement and content identical to the most common type in insects. The start codon of the C. fulgidissima COI gene is unusual, in that no typical ATN codon is available. The 875 bp A+T-rich region is the shortest among the coleopteran mitogenomes that have thus far been sequenced in their entirety. The most unusual feature of the genome is the presence of three tRNA-like sequences within the A+T-rich region: two tRNA(Leu)(UUR)-like sequences and one tRNA(Asn)-like sequence. These sequence stretches have the proper anticodon sequence and the potential to form secondary structures, but also harbor many mismatches in the stems. Phylogenetic analysis using a concatenation of 13 amino acid sequences of protein-coding genes among the available sequenced species of coleopteran superfamilies (Buprestoidea and Elateroidea belonging to the infraorder Elateriformia, and Chrysomeloidea and Tenebrioroidea belonging to the infraorder Cucujiformia) by Bayesian inference, maximum-likelihood analysis, and maximum-parsimony analysis revealed a lack of strong support for monophyletic Elateriformia.

Complete nucleotide sequence and organization of the mitogenome of the red-spotted apollo butterfly, Parnassius bremeri (Lepidoptera: Papilionidae) and comparison with other lepidopteran insects.

The 15,389-bp long complete mitogenome of the endangered red-spotted apollo butterfly, Parnassius bremeri (Lepidoptera: Papilionidae) was determined in this study. The start codon for the COI gene in insects has been extensively discussed, and has long remained a matter of some controversy. Herein, we propose that the CGA (arginine) sequence functions as the start codon for the COI gene in lepidopteran insects, on the basis of complete mitogenome sequences of lepidopteran insects, including P. bremeri, as well as additional sequences of the COI start region from a diverse taxonomic range of lepidopteran species (a total of 53 species from 15 families). In our extensive search for a tRNA-like structure in the A+T-rich region, one tRNA(Trp)-like sequence and one tRNA(Leu) (UUR)-like sequence were detected in the P. bremeri A+T-rich region, and one or more tRNA-like structures were detected in the A+T-rich region of the majority of other sequenced lepidopteran insects, thereby indicating that such features occur frequently in the lepidopteran mitogenomes. Phylogenetic analysis using the concatenated 13 amino acid sequences and nucleotide sequences of PCGs of the four macrolepidopteran superfamilies together with the Tortricoidea and Pyraloidea resulted in the successful recovery of a monophyly of Papilionoidea and a monophyly of Bombycoidea. However, the Geometroidea were unexpectedly identified as a sister group of the Bombycoidea, rather than the Papilionoidea.