Ultraconserved elements from transcriptome and genome data provide insight into the phylogenomics of Sternorrhyncha (Insecta: Hemiptera).

Sternorrhyncha, one of the four major suborders of Hemiptera, is a phytophagous taxon inclusive of nearly 18 000 described species. The phylogenetic relationships within the taxon and the earliest-branching lineage of its infraorders remain incompletely understood. This study attempted to illuminate the phylogenetic relationships within Sternorrhyncha through the use of maximum likelihood, Bayesian inference and maximum parsimony analyses, employing ultraconserved element (UCE) data from 39 genomic and 62 transcriptomic datasets and thereby representing most families within the taxon. The probe set Hemiptera 2.7Kv1 was used to recover a total of 2731 UCE loci: from 547 to 1699 (with an average of 1084) across all genomic datasets and from 108 to 849 (with an average of 329) across all transcriptomic datasets. All three types of phylogenetic analyses employed in this study produced robust statistical support for Sternorrhyncha being a monophyletic group. The different methods of phylogenetic analysis produced inconsistent descriptions of topological structure at the infraorder level: while maximum likelihood and Bayesian inference analyses produced strong statistical evidence (100%) indicating the clade Psylloidea + Aleyrodoidea to be a sister of the clade Aphidoidea (Aphidomorpha) + Coccoidea (Coccomorpha), the maximum parsimony analysis failed to recover a similar result. Our results also provide detail on the phylogenetic relationships within each infraorder. This study presents the first use of UCE data to investigate the phylogeny of Sternorrhyncha. It also shows the viability of amalgamating genomic and transcriptomic data in studies of phylogenetic relationships, potentially highlighting a resource-efficient approach for future inquiries into diverse taxa through the integration of varied data sources.

Role of dietary nutrients and metabolism in colorectal cancer.

Colorectal cancer (CRC) is one of the most common malignancies and the leading causes of cancer related deaths worldwide. The development of CRC is driven by a combination of genetic and environmental factors. There is growing evidence that changes in dietary nutrition may modulate the CRC risk, and protective effects on the risk of developing CRC have been advocated for specific nutrients such as glucose, amino acids, lipid, vitamins, micronutrients and prebiotics. Metabolic crosstalk between tumor cells, tumor microenvironment components and intestinal flora further promote proliferation, invasion and metastasis of CRC cells and leads to treatment resistance. This review summarizes the research progress on CRC prevention, pathogenesis, and treatment by dietary supplementation or deficiency of glucose, amino acids, lipids, vitamins, micronutri-ents, and prebiotics, respectively. The roles played by different nutrients and dietary crosstalk in the tumor microenvironment and metabolism are discussed, and nutritional modulation is inspired to be beneficial in the prevention and treatment of CRC.

Advances in pelvic imaging parameters predicting surgical difficulty in rectal cancer.

Due to the fixed bony structure of the pelvis, the pelvic operation space is limited, complicating the surgical operation of rectal cancer, especially middle and low rectal cancer. The closer the tumor is to the anal verge, the smaller the operative field and operating space, the longer the operative time, and the greater the incidence of intraoperative side injuries and postoperative complications. To date, there is still no clear definition of a difficult pelvis that affects the surgical operation of rectal cancer. Few related research reports exist in the literature, and views on this aspect are not the same between countries. Therefore, it is particularly important to predict the difficulty of rectal cancer surgery in a certain way before surgery and to select the surgical method most suitable for each case during the treatment of rectal cancer.

KIF15 facilitates gastric cancer via enhancing proliferation, inhibiting apoptosis, and predict poor prognosis.

BACKGROUND: Kinesin superfamily proteins (KIFs) can transport membranous organelles and protein complexes in an ATP-dependent manner. Kinesin family member 15 (KIF15) is overexpressed in various cancers. However, the function of KIF15 in gastric cancer (GC) is still unclear. METHODS: GC patients' data from The Cancer Genome Atlas (TCGA) were analyzed by bioinformatics methods. The expression of KIF15 was examined in GC and paracarcinoma tissues from 41 patients to verify the analysis results. The relationship between KIF15 expression and clinical characteristics were also observed by bioinformatics methods. Kaplan-Meier survival analysis of 122 GC patients in our hospital was performed to explore the relationship between KIF15 expression levels and GC patients' prognosis. KIF15 was downregulated in GC cell lines AGS and SGC-7901 by transfecting a lentivirus-mediated shRNA plasmid targeting KIF15. In vitro, GC cell proliferation and apoptosis were detected by MTT assay, colony formation assay, and Annexin V-APC staining. In vivo, xenograft experiments were used to verify the in vitro results. Furthermore, Human Apoptosis Antibody Array kit was used to screen possible targets of KIF15 in GC cell lines. RESULTS: The bioinformatics results showed that KIF15 expression levels were higher in GC tissues than in normal tissues. IHC showed same results. High expression of KIF15 was statistical correlated with high age and early histologic stage. Kaplan-Meier curves indicated that high KIF15 expression predict poor prognosis in patients with GC. MTT assay and colony formation assay showed that KIF15 promote GC cell proliferation. Annexin V-APC staining found that KIF15 can inhibit GC cell apoptosis. Xenograft experiments reveal that downregulating KIF15 can inhibit GC tumor growth and promote GC apoptosis. Through detection of 43 anti-apoptotic proteins by the Human Apoptosis Antibody Array kit, it was confirmed that knocking down KIF15 can reduce seven anti-apoptotic proteins expression. CONCLUSIONS: Taken together, our study revealed a critical role for KIF15 to inhibit GC cell apoptosis and promote GC cell proliferation. KIF15 may decrease anti-apoptotic proteins expression by regulating apoptosis pathways. High expression of KIF15 predicts a poor prognosis in patients with GC. KIF15 might be a novel prognostic biomarker and a therapeutic target for GC.

Circular RNA circ-DONSON facilitates gastric cancer growth and invasion via NURF complex dependent activation of transcription factor SOX4.

BACKGROUND: Circular RNAs (circRNAs) are a novel type of noncoding RNAs and play important roles in tumorigenesis, including gastric cancer (GC). However, the functions of most circRNAs remain poorly understood. In our study, we aimed to investigate the functions of a new circRNA circ-DONSON in GC progression. METHODS: The expression of circ-DONSON in gastric cancer tissues and adjacent normal tissues was analyzed by bioinformatics method, qRT-PCR, Northern blotting and in situ hybridization (ISH). The effects of circ-DONSON on GC cell proliferation, apoptosis, migration and invasion were measured by using CCK8, colony formation, EdU, immunofluorescence (IF), FACS and Transwell assays. qRT-PCR and Western blotting were utilized to validate how circ-DONSON regulates SOX4 expression. ChIP, DNA fluorescence in situ hybridization (DNA-FISH) and DNA accessibility assays were used to investigate how circ-DONSON regulates SOX4 transcription. The interaction between circ-DONSON and NURF complex was evaluated by mass spectrum, RNA immunoprecipitation (RIP), pulldown and EMSA assays. Xenograft mouse model was used to analyze the effect of circ-DONSON on GC growth in vivo. RESULTS: Elevated expression of circ-DONSON was observed in GC tissues and positively associated with advanced TNM stage and unfavorable prognosis. Silencing of circ-DONSON significantly suppressed the proliferation, migration and invasion of GC cells while promoting apoptosis. circ-DONSON was localized in the nucleus, recruited the NURF complex to SOX4 promoter and initiated its transcription. Silencing of the NURF complex subunit SNF2L, BPTF or RBBP4 similarly attenuated GC cell growth and increased apoptosis. circ-DONSON knockdown inhibited GC growth in vivo. CONCLUSION: circ-DONSON promotes GC progression through recruiting the NURF complex to initiate SOX4 expression.

Comparative Mitogenomic Analysis of the Eurydema Genus in the Context of Representative Pentatomidae (Hemiptera: Heteroptera) Taxa.

The family Pentatomidae, the largest within the superfamily Pentatomoidae, comprises about 5,000 species; many of which are economically important pests. Although the phylogeny of Pentatomidae species has been studied using various molecular markers, their phylogenetic relationships remain controversial. Recently, mitochondrial genomes (mitogenomes) have been extensively employed to examine the phylogenetics and evolution of different insects, and in this study, we sequenced complete/near-complete mitochondrial genomes from five shield bug species of Eurydema to gain a better understanding of phylogenetic relationships in the Pentatomidae. The five mitogenomes ranged in length from 15,500 to 16,752 bp and comprised 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), 2 ribosomal RNAs (rRNAs), and a control region. We compared mitogenomic characteristics of the Pentatomidae and constructed phylogenetic trees using Bayesian inference and maximum likelihood methods. Our results showed that gene arrangements, base composition, start/stop codons, gene overlaps, and RNA structures were conserved within the Pentatomidae and that congeneric species shared more characteristics. Saturation and heterogeneity analyses revealed that our PCGs and PCGRNA datasets were valid for phylogenetic analysis. Phylogenetic analyses showed consistent topologies based on BI and ML methods. These analyses strongly supported that Eurydema species belong to the tribe Strachiini, and formed a sister group with Pentatomini. The relationships among Eurydema species were shown to be consistent with their morphological features. (Strachiini + Pentatomini) was found to be a stable sibling of the clade comprising Cappaeini, Graphosomini, and Carpocorini. Furthermore, our results indicated that Graphosoma rubrolineatum (Heteroptera: Pentatomidae) belongs to the Pentatominae and not the Podopinae.

Complete Mitochondrial Genome of Dinorhynchus dybowskyi (Hemiptera: Pentatomidae: Asopinae) and Phylogenetic Analysis of Pentatomomorpha Species.

Dinorhynchus dybowskyi (Hemiptera: Pentatomidae: Asopinae) is used as a biological control agent against various insect pests for its predatory. In the present study, the complete mitochondrial genome (mitogenome) of the species was sequenced using the next-generation sequencing technology. The results showed that the mitogenome is 15,952 bp long, including 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs), and a control region. Furthermore, the gene order and orientation of this mitogenome are identical to those of most heteropterans. There are 21 intergenic spacers (of length 1-28 bp) and 13 overlapping regions (of length 1-23 bp) throughout the genome. The control region is 1,291 bp long. The start codon of the PCGs is ATN, except cox1 (TTG), and stop codon is TAA, except nad1 (TAG). The 22 tRNAs exhibit a typical cloverleaf secondary structure, except trnS1, which lacks a dihydrouridine (DHU) arm and trnV, where the DHU arm forms a simple loop. The analyses based on nucleotide sequences of the 13 PCGs by Bayesian Inference and maximum likelihood methods. The results support the monophyly of five superfamilies Aradoidea, Pentatomoidea, Pyrrhocoroidea, Lygaeoidea, and Coreoidea. Within Pentatomoidea, the relationship observed is as follows: (Plataspidae + Urostylididae) + (Pentatomidae + (Acanthosomatidae + (Cydnidae + (Scutelleridae + (Dinidoridae + Tessaratomidae))))), and D. dybowskyi was placed in Pentatomidae and close to Eurydema gebleri.

Overexpression of Wnt7α protein predicts poor survival in patients with colorectal carcinoma.

Wnt7α (wingless-type MMTV integration site family, member 7A) is a secreted glycoprotein that plays a critical role in tumorigenesis and development by controlling cell proliferation and differentiation. Whether Wnt7α has the properties of an oncogene or not is an interesting issue because of its diverse expression in different tumors. In the present study, Wnt7α protein expression was evaluated through immunohistochemistry and Western blot analysis. Univariate and multivariate analyses were applied to explore the associations between Wnt7α staining score and various clinical parameters, including overall survival (OS) and disease-free survival (DFS), and a total of 212 patients with colorectal cancer (CRC) were surveyed. Wnt7α was strongly expressed in most CRC tissues but weakly expressed in adjacent normal mucosa, colorectal adenomas, and colonic polyps. High levels of Wnt7α expression were strongly associated with tumor size (P = 0.006), lymph node involvement (P < 0.001), and the international tumor-node-metastasis (TNM) stage (P = 0.005). Patients with strong Wnt7α expression showed significantly poorer OS and DFS than patients with weak Wnt7α expression (P < 0.0001, both). Multivariate Cox analysis confirmed that Wnt7α protein expression and TNM stage are independent factors of adverse OS and DFS in CRC patients. Taken together, our results present evidence that Wnt7α overexpression is associated with an unfavorable prognosis and that positive Wnt7α, in addition to TNM stage, may be an independent prognosis factor influencing OS and DFS prediction in CRC patients.

Improving DNA barcoding library of armored scale insects (Hemiptera: Diaspididae) in China.

DNA barcoding is used to identify cryptic species, survey environmental samples, and estimate phyletic and genetic diversity. Armored scale insects are phytophagous insects and are the most species-rich taxa in the Coccoidea superfamily. This study developed a DNA barcode library for armored scale insect species collected from southern China during 2021-2022. We sequenced a total of 239 specimens, recognized as 50 morphological species, representing two subfamilies and 21 genera. Sequencing analysis revealed that the average G + C content of the cytochrome oxidase subunit I (COI) gene sequence was very low (~18.06%) and that the average interspecific divergence was 10.07% while intraspecific divergence was 3.20%. The intraspecific divergence value was inflated by the high intraspecific divergence in ten taxa, which may indicate novel species overlooked by current taxonomic treatments. All the Automated Barcode Gap Discovery, Assemble Species by Automatic Partitioning, Taxon DNA analysis and Bayesian Poisson Tree Process methods yielded largely consistent results, indicating a robust and credible species delimitation. Based on these results, an intergeneric distance threshold of ≤ 5% was deemed appropriate for the differentiation of armored scale insect species in China. This study establishes a comprehensive barcode library for the identification of armored scale insects, future research, and application.

Global Potential Distribution of Invasive Species Pseudococcus viburni (Hemiptera: Pseudococcidae) under Climate Change.

The potential distribution range and management strategies for P. viburni are poorly understood. Based on historical distribution data and environmental factors, the present study predicted the potentially suitable areas for P. viburni spread under different climate change scenarios using MaxEnt (maximum entropy). The results showed that precipitation of the coldest quarter (Bio19), precipitation seasonality (Bio15), and mean temperature of the wettest quarter (Bio8) were the most important environmental factors determining the distribution of P. viburni. Under the current climate conditions, its potential suitable areas are southern China, the whole of Japan, North America (especially the eastern part of the United States), the southwestern part of South America, the Mediterranean coast and most of Europe, the central part of Africa, i.e., the south of the Sahara Desert, and most of the southern coast of Australia. The total area of habitats suitable for this insect pest is predicted to be increased in the future. In order to prevent P. viburni transmission and spread, there is a need to strengthen the monitoring and quarantine measures against this pest at the Southern ports.

Novel insights into hotspots of insect vectors of GLRaV-3: Dynamics and global distribution.

Grapevine leafroll-associated virus 3 (GLRaV-3) is the most prevalent and economically damaging virus in grapevines and is found on nearly all continents, except Antarctica. Ten mealybugs act as vector insects transmitting the GLRaV-3. Understanding the potential distribution range of vector insects under climate change is crucial for preventing and managing vector insects and controlling and delaying the spread of GLRaV-3. This study investigated the potential geographical range of insect vectors of GLRaV-3 worldwide using MaxEnt (maximum entropy) based on occurrence data under environmental variables. The potential distributions of these insects were projected for the 2030s, 2050s, 2070s, and 2090s under the three climate change scenarios. The results showed that the potential distribution range of most vector insects is concentrated in Southeastern North America, Europe, Asia, and Southeast Australia. Most vector insects contract their potential distribution ranges under climate-change conditions. The stacked model suggested that potential distribution hotspots of vector insects were present in Southeastern North America, Europe, Southeast Asia, and Southeast Australia. The potential distribution range of hotspots would shrink with climate change. These results provide important information for governmental decision-makers and farmers in developing control and management strategies against vector insects of GLRaV-3. They can also serve as references for studies on other insect vectors.

Potential Global Distribution of Paracoccus marginatus, under Climate Change Conditions, Using MaxEnt.

The papaya mealybug, Paracoccus marginatus, is an invasive pest species found all over the world. It is native to Mexico and Central America, but is now present in more than 50 countries and regions, seriously threatening the economic viability of the agricultural and forestry industry. In the current study, the global potential distribution of P. marginatus was predicted under current and future climatic conditions using MaxEnt. The results of the model assessment indicated that the area under the curve of the receiver operating characteristic ( ROC-AUC) was 0.949, while the TSS value was 0.820. The results also showed that the three variables with the greatest impact on the model were min temperature of coldest month (bio6), precipitation of wettest month (bio13), and precipitation of coldest quarter (bio19), with corresponding contributions of 46.8%, 31.1%, and 13.1%, respectively. The results indicated that the highly suitable areas were mainly located in tropical and subtropical regions, including South America, southern North America, Central America, Central Africa, Australia, the Indian subcontinent, and Southeast Asia. Under four climate scenarios in the 2050s and 2070s, the area of suitability will change very little. Moreover, the results showed that the area of suitable areas in 2070s increased under all four climate scenarios compared to the current climate. In contrast, the area of suitable habitat increases from the current to the 2050s under the SSP370 and SSP585 climate scenarios. The current study could provide a reference framework for the future control and management of papaya mealybug and other invasive species.

Two new species of Fiorinia Targioni-Tozzetti (Hemiptera: Coccoidea: Diaspididae) from China.

The adult females of two new species of armored scale, Fiorinia tianshuiensis Wei & Feng sp. n. and F. dinghuensis Wei & Feng sp. n., are described and illustrated from specimens collected from China. A key is provided for the all Fiorilim species known from China.

A new species of Aulacaspis and a revived combination of Diaspididae (Hemiptera, Coccomorpha) from China.

A new species of armored scale insect, Aulacaspisfanjingshanensissp. nov. is described and illustrated based on adult female specimens collected on Rosaceae plants in China. A key to the Aulacaspis species known from Guizhou Province of China is provided. Our molecular study suggests that Aulacaspisschizosoma (Takagi, 1970) is not a true member of the genus Aulacaspis; the genus Superturmaspis Chen, 1983 is revived and A.schizosoma is transferred to it as Superturmaspisschizosoma (Takagi, 1970), revived combination, based on a molecular phylogeny.

Comparative mitogenomics and phylogenetic analyses of the genus Menida (Hemiptera, Heteroptera, Pentatomidae).

In order to explore the genetic diversity and phylogenetic relationship of the genus Menida Motschulsky, 1861 and reveal the molecular evolution of the family Pentatomidae, subfamily Pentatominae, complete mitochondrial genomes of three species of Menida were sequenced, and the phylogenetic relationships of tribes within the subfamily Pentatominae were studied based on these results. The mitochondrial genomes of Menidamusiva (Jakovlev, 1876), M.lata Yang, 1934, and M.metallica Hsiao & Cheng, 1977 were 16,663 bp, 16,463 bp, and 16,418 bp, respectively, encoding 37 genes and including 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes, and a control region. The mitochondrial genome characteristics of Menida were compared and analyzed, and the phylogenetic tree of the Pentatominae was constructed based on the mitochondrial genome datasets using Bayesian inference (BI) and maximum likelihood (MI) methods. The results showed that gene arrangements, nucleotide composition, codon preference, gene overlaps, and RNA secondary structures were highly conserved within the Menida and had more similar characteristics in Pentatominae. The phylogenetic analysis shows a highly consistent topological structure based on BI and ML methods, which supported that the genus Menida belongs to the Pentatominae and is closely related to Hoplistoderini. The examined East Asian species of Menida form a monophyletic group with the internal relationships: (M.musiva + (M.lata + (M.violacea + M.metallica))). In addition, these results support the monophyly of Eysarcorini and Strachiini. Placosternum and Cappaeini are stable sister groups in the evolutionary branch of Pentatominae. The results of this study enrich the mitochondrial genome databases of Pentatominae and have significance for further elucidation of the phylogenetic relationships within the Pentatominae.

Invasion risks presented by Gonopsis affinis and the use of Trissolcus mitsukurii as a biological control agent under present and future climate conditions.

BACKGROUND: Gonopsis affinis (Uhler) is a stinkbug that represents a significant threat to the production of rice (Oryza sativa L.), sugarcane (Saccharum officinarum L.) and eulalia (Miscanthus sinensis (Andersson)), and has been listed as a sugarcane pest in Japan. Trissolcus mitsukurii Ashmead is an egg parasitoid of G. affinis. To determine the potential of T. mitsukurii to be a biological control agent for G. affinis, we aim to predict the current and future areas of suitable habitat for these two species and their overlap with areas of present crop production. We developed MaxEnt models using two different variable selection methods and compared the two for T. mitsukurii with a CLIMEX model. RESULTS: The results showed extensive suitable areas for G. affinis under current climate conditions in East Asia, West Africa, Madagascar, and South America. These ranges overlap with areas currently being used for the production of the three crops in question. More than half overlap with areas of suitable habitat for T. mitsukurii. The most critical environmental variable determining habitat suitability for G. affinis was showed to be precipitation of warmest quarter, whilst for T. mitsukurii it was minimum temperature of the coldest month. CONCLUSION: Based on our assessment we recommend the immediate implementation of monitoring and invasion prevention measures for G. affinis in southwest China, the Malay Archipelago and West Africa. We suggest that T. mitsukurii be considered for use as a biological control agent in East Asia, Madagascar, Florida and Brazil in the case of future invasions by G. affinis. © 2023 Society of Chemical Industry.

The complete mitochondrial genome of Aeschrocoristuberculatus and A.ceylonicus (Hemiptera, Pentatomidae) and its phylogenetic implications.

Aeschrocoristuberculatus and A.ceylonicus (Hemiptera, Pentatomidae, Pentatominae) are mainly distributed in southern China, India, Myanmar, and Sri Lanka. Both species are also common agricultural pests. However, only the morphology of the genus Aeschrocoris has previously been studied, and molecular data have been lacking. In this study, the whole mitochondrial genomes of A.tuberculatus and A.ceylonicus are and annotated. The lengths of the complete mitochondrial genomes of the two species are 16,134 bp and 16,142 bp, respectively, and both contain 37 typical genes, including 13 protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and a control region. The mitochondrial genome structure, gene order, nucleotide composition, and codon usage of A.tuberculatus and A.ceylonicus are consistent with those of typical Pentatomidae. Most PCGs of both species use ATN as the start codon, except atp8, nad1, and cox1, which use TTG as the start codon. cox1, cox2, and atp6 use a single T, and nad1 use TAG as the stop codon; the remaining PCGs have TAA as the stop codon. The A+T contents of the two species are 73.86% and 74.08%, respectively. All tRNAs have a typical cloverleaf structure, with the exception of trnS1, which lacks a dihydrouridine arm. The phylogenetic tree is reconstructed using the maximum-likelihood method based on the newly obtained mitochondrial genome sequences and 87 existing mitochondrial genomes of Pentatomoidea from the NCBI database and two species of Lygaeoidea as outgroups. The phylogenetic trees strongly support the following relationships: (Urostylididae + ((Acanthosomatidae + ((Cydnidae + (Dinidoridae + Tessaratomidae)) + (Scutelleridae + Plataspidae))) + Pentatomidae). This study enriches the mitochondrial genome database of Pentatomoidea and provides a reference for further phylogenetic studies.

Two Species Delimitation of Pseudaulacaspis (Hemiptera: Diaspididae) Based on Morphology, Molecular Clustering, and Niche Differentiation.

Pseudaucalaspis pentagona and P. prunicola are notorious pests and commonly feed on various ornamental plants and fruit trees worldwide. The two species share many host-plant species, and are similar in morphological characteristics and life cycle, making it difficult to distinguish to distinguish between them. In this study, morphological characteristics, molecular evidence, and ecological niches were used to define these species. We performed PCA analysis on 22 morphological characteristics that allowed the delineation of the species. We then sequenced the COI gene of both species revealing five populations of P. pentagona and one population of P. prunicola, and the higher support rate could distinguish the two species. We also identified the potential distribution area of the two species based on the MaxEnt niche model, which showed that the degree of niche overlap was high, but that they occupied different niches. Ultimately, we combined three lines of evidence to show that the two species are distinctly different. This study supports species definition using combined morphology, genetics, and ecology and provides a theoretical basis for the effective control of these two pests in the future.

First complete mitochondrial genome of the tribe Coccini (Hemiptera, Coccomorpha, Coccidae) and its phylogenetic implications.

Soft scale insects (Hemiptera, Coccidae) are important pests of various agricultural and horticultural crops and ornamental plants. They have negative impacts on agriculture and forestry. The tribe Coccini represents one of the most ancient evolutionary lineages of soft scale insects. However, no complete Coccini mitochondrial genome (mitogenome) is available in public databases. Here, we described the complete mitogenome of Coccushesperidum L., 1758. The 15,566 bp mitogenome of C.hesperidum had a high A+T content (83.4%) and contained a typical set of 37 genes, with 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs) and two ribosomal RNA genes (rRNAs). Only seven tRNAs had the typical clover-leaf secondary structure and the remaining tRNAs lacked the DHU arm, TψC arm or both. Moreover, a comparative analysis of all reported scale insect mitogenomes from GenBank database was performed. The mitogenomes of scale insects showed high similarities in base composition and A+T content. Additionally, our phylogenetic analysis confirmed the monophyly of Coccomorpha and revealed that the archaeococcoids were the most basal lineage within Coccomorpha, while Ericeruspela and Didesmococcuskoreanus, belonging to Coccidae, were often mixed with Aclerdidae, making Coccidae a paraphyletic group. These findings expand the mitogenome database of scale insects and provide new insights on mitogenome evolution for future studies across different insect groups.

Comparative analysis and phylogeny of mitochondrial genomes of Pentatomidae (Hemiptera: Pentatomoidea).

The Phyllocephalini is a group of herbivorous insects in Pentatomidae, which lack distinctive morphological characteristics and systematic studies. Up to now, there are only two complete mitochondrial genomes of Phyllocephalini have been reported. In this study, we sequenced and analyzed the complete mitochondrial genomes of three Phyllocephalini species, Gonopsis coccinea, Gonopsimorpha nigrosignata, and Chalcopis glandulosus, which were 16,534, 16,531, and 16,534 bp in length, respectively. The mitochondrial genomes contained 37 genes, including 13 protein-coding genes, two rRNA genes, 22 tRNA genes, and a control region. The gene arrangement was consistent with that of the putative ancestral insect, with no rearrangement. The cox1 gene of Pentatomidae showed the lowest evolutionary rate among the protein-coding genes, the mean genetic distance of species, genera, and subfamilies of Pentatomidae increased hierarchically based on cox1 gene. The 16S rRNA of Pentatomidae was more conserved than 12S rRNA in sequence and secondary structure. All tRNAs could be folded into a typical cloverleaf structure except trnS1. The stem region was more conserved than the loop region in the secondary structure of tRNAs within Pentatomidae. Gonopsis coccinea and Gonopsimorpha nigrosignata had one type of tandem repetition unit in the control region, while C. glandulosus had two types. The heterogeneity analysis of Pentatomidae showed that Phyllocephalinae was the most heterogeneous. Phylogenetic trees based on the newly obtain mitochondrial genomes along with other 50 mitochondrial genomes of Pentatomidae using Bayesian Inference and Maximum Likelihood strongly supported the following three relationships: (((Anaxilaus + (Plautia + Glaucias)) + (Nezara + Palomena)) + (Eysarcorini + Carpocorini)), (Hoplistoderini + (Menidini + Asopinae)), and ((Sephelini + Halyini) + (Caystrini + (Cappaeini + (Placosternum + Phyllocephalini)))). The relationships within Phyllocephalini were (Chalcopis + (Dalsira + (Gonopsimorpha + Gonopsis))). Our results provide valuable molecular data for further phylogenetic analyses of Pentatomidae.