Search

Microbial symbioses have had profound impacts on the evolution of animals. Conversely, changes in host biology may impact the evolutionary trajectory of symbionts themselves. Blattabacterium cuenoti is present in almost all cockroach species and enables hosts to subsist on a nutrient-poor diet. To investigate if host biology has impacted Blattabacterium at the genomic level, we sequenced and analyzed 25 genomes from Australian soil-burrowing cockroaches (Blaberidae: Panesthiinae), which have undergone at least seven separate subterranean, subsocial transitions from above-ground, wood-feeding ancestors. We find at least three independent instances of genome erosion have occurred in Blattabacterium strains exclusive to Australian soil-burrowing cockroaches. These shrinkages have involved the repeated inactivation of genes involved in amino acid biosynthesis and nitrogen recycling, the core role of Blattabacterium in the host-symbiont relationship. The most drastic of these erosions have occurred in hosts thought to have transitioned underground the earliest relative to other lineages, further suggestive of a link between gene loss in Blattabacterium and the burrowing behavior of hosts. As Blattabacterium is unable to fulfill its core function in certain host lineages, these findings suggest soil-burrowing cockroaches must acquire these nutrients from novel sources. Our study represents one of the first cases, to our knowledge, of parallel host adaptations leading to concomitant parallelism in their mutualistic symbionts, further underscoring the intimate relationship between these two partners., (© 2024 The Authors. Insect Science published by John Wiley & Sons Australia, Ltd on behalf of Institute of Zoology, Chinese Academy of Sciences.)

Despite the many advances of the genomic era, there is a persistent problem in assessing the uncertainty of phylogenomic hypotheses. We see this in the recent history of phylogenetics for cockroaches and termites (Blattodea), where huge advances have been made, but there are still major inconsistencies between studies. To address this, we present a phylogenetic analysis of Blattodea that emphasizes identification and quantification of uncertainty. We analyze 1183 gene domains using three methods (multi-species coalescent inference, concatenation, and a supermatrix-supertree hybrid approach) and assess support for controversial relationships while considering data quality. The hybrid approach-here dubbed "tiered phylogenetic inference"-incorporates information about data quality into an incremental tree building framework. Leveraging this method, we are able to identify cases of low or misleading support that would not be possible otherwise, and explore them more thoroughly with follow-up tests. In particular, quality annotations pointed towards nodes with high bootstrap support that later turned out to have large ambiguities, sometimes resulting from low-quality data. We also clarify issues related to some recalcitrant nodes: Anaplectidae's placement lacks unbiased signal, Ectobiidae s.s. and Anaplectoideini need greater taxon sampling, the deepest relationships among most Blaberidae lack signal. As a result, several previous phylogenetic uncertainties are now closer to being resolved (e.g., African and Malagasy "Rhabdoblatta" spp. are the sister to all other Blaberidae, and Oxyhaloinae is sister to the remaining Blaberidae). Overall, we argue for more approaches to quantifying support that take data quality into account to uncover the nature of recalcitrant nodes., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Published by Elsevier Inc.)

Cockroaches are well-known pests and quarantined organisms worldwide. Due to morphological diversity and a lack of molecular data, their classification and identification are facing challenges. This study performed classification, identification, and DNA barcoding for cockroaches collected from China. Seventy-six samples were morphologically identified as seven species of two superfamilies that included Blattella germanica, Eublaberus posticus and Blaptica dubia belonging to the superfamily Blaberoidea, and Periplaneta americana, Periplaneta lateralis, Periplaneta fuliginosa and Periplaneta australasiae belonging to the superfamily Blattoidea. Based on sequence alignments of nine ribosomal and mitochondrial genes across the order Blattaria retrieved from GenBank, rDNA ITS2-517 bp and mtDNA 16S-327 bp were screened as candidates for molecular identification. Universal primers were designed for PCR amplification, cloning, and sequencing of the 37 representative samples. Sequence alignments and phylogeny analysis showed that both ITS2 and 16S confirmed samples 1-9, 20-24, and 25-29 as B. germanica, P. americana, and P. lateralis, respectively; only 16S (not ITS2) confirmed samples 10-14, 15-19, 30-34, and 35-37 as E. posticus, Blap. dubia, P. fuliginosa, and P. australasiae, respectively, indicating that 16S was a better target than ITS2 for molecular identification of cockroaches. Conservative motif and divergence analysis further revealed that ITS2 sequences vary significantly among different taxa, whereas 16S sequences are relatively conserved. There is an obvious barcoding gap between maximum intraspecific divergence and minimum interspecific divergence (2.57 % vs. 5.62 %) for ITS2, but not for 16S (6.15 % vs. 2.63 %). Therefore, it was confirmed that ITS2 is an ideal DNA barcode for molecular identification of cockroaches at lower category., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier B.V. All rights reserved.)

Rates of nucleotide substitution vary substantially across the Tree of Life, with potentially confounding effects on phylogenetic and evolutionary analyses. A large acceleration in mitochondrial substitution rate occurs in the cockroach family Nocticolidae, which predominantly inhabit subterranean environments. To evaluate the impacts of this among-lineage rate heterogeneity on estimates of phylogenetic relationships and evolutionary timescales, we analyzed nuclear ultraconserved elements (UCEs) and mitochondrial genomes from nocticolids and other cockroaches. Substitution rates were substantially elevated in nocticolid lineages compared with other cockroaches, especially in mitochondrial protein-coding genes. This disparity in evolutionary rates is likely to have led to different evolutionary relationships being supported by phylogenetic analyses of mitochondrial genomes and UCE loci. Furthermore, Bayesian dating analyses using relaxed-clock models inferred much deeper divergence times compared with a flexible local clock. Our phylogenetic analysis of UCEs, which is the first genome-scale study to include all 13 major cockroach families, unites Corydiidae and Nocticolidae and places Anaplectidae as the sister lineage to the rest of Blattoidea. We uncover an extraordinary level of genetic divergence in Nocticolidae, including two highly distinct clades that separated ~115 million years ago despite both containing representatives of the genus Nocticola. The results of our study highlight the potential impacts of high among-lineage rate variation on estimates of phylogenetic relationships and evolutionary timescales., (© The Author(s) 2024. Published by Oxford University Press on behalf of the Society of Systematic Biologists.)

Gene duplication is a fundamental evolutionary process which provides opportunities to acquire new gene functions. In the case of the insulin receptors (InRs) in cockroaches and close-related insects, two successive duplications determined the occurrence of three InR genes: InR2, InR1 and InR3, the last two forming a sister cluster to InR2. The biological role of each of the gene duplicates and whether they resulted from neofunctionalization or subfunctionalization is still unclear. The analysis of the sequences from different lineages did not detect positive selection as driving the divergence of InR1 and InR3, discarding neofunctionalization, and suggesting that there is no functional divergence between both gene copies. Using the cockroach Blattella germanica as a model, we have determined that BgInR2 is the gene with the highest expression levels in all the tissues analyzed, both in adult females and males, as well as in nymphs and embryos. BgInR3 is second in expression levels while BgInR1 is expressed at lower levels and only in some tissues. The selective depletion by RNAi of each of the three InRs, analyzed in terms of phenotype and fat body transcriptomic profiles, resulted in essentially redundant effects, with a magnitude approximately proportional to the level of expression of the respective InR. Therefore, the results indicate that the InR duplicates likely experienced a subfunctionalization process, by which the three InRs maintained similar functions but contributing to those functions proportionally to their expression levels., (Copyright © 2024 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Genetic changes that enabled the evolution of eusociality have long captivated biologists. More recently, attention has focussed on the consequences of eusociality on genome evolution. Studies have reported higher molecular evolutionary rates in eusocial hymenopteran insects compared with their solitary relatives. To investigate the genomic consequences of eusociality in termites, we analysed nine genomes, including newly sequenced genomes from three non-eusocial cockroaches. Using a phylogenomic approach, we found that termite genomes have experienced lower rates of synonymous substitutions than those of cockroaches, possibly as a result of longer generation times. We identified higher rates of non-synonymous substitutions in termite genomes than in cockroach genomes, and identified pervasive relaxed selection in the former (24-31% of the genes analysed) compared with the latter (2-4%). We infer that this is due to reductions in effective population size, rather than gene-specific effects (e.g. indirect selection of caste-biased genes). We found no obvious signature of increased genetic load in termites, and postulate efficient purging of deleterious alleles at the colony level. Additionally, we identified genomic adaptations that may underpin caste differentiation, such as genes involved in post-translational modifications. Our results provide insights into the evolution of termites and the genomic consequences of eusociality more broadly.

The mitogenome sequence data have been widely used in inferring the phylogeny of insects. In this study, we determined the complete mitogenome for Macrotermes sp. (Termitidae, Macrotermitinae) using next-generation sequencing. Macrotermes sp. possesses a typical insect mitogenome, displaying an identical gene order and gene content to other existing termite mitogenomes. We present the first prediction of the secondary structure of ribosomal RNA genes in termites. The rRNA secondary structures of Macrotermes sp. exhibit similarities to closely related insects and also feature distinctive characteristics in their helical structures. Together with 321 published mitogenomes of termites as ingroups and 8 cockroach mitogenomes as outgroups, we compiled the most comprehensive mitogenome sequence matrix for Termitoidae to date. Phylogenetic analyses were conducted using datasets employing different data coding strategies and various inference methods. Robust relationships were recovered at the family or subfamily level, demonstrating the utility of comprehensive mitogenome sampling in resolving termite phylogenies. The results supported the monophyly of Termitoidae, and consistent relationships within this group were observed across different analyses. Mastotermitidae was consistently recovered as the sister group to all other termite families. The families Hodotermitidae, Stolotermitidae, and Archotermopsidae formed the second diverging clade, followed by the Kalotermitidae. The Neoisoptera was consistently supported with strong node support, with Stylotermitidae being sister to the remaining families. Rhinotermitidae was found to be non-monophyletic, and Serritermitidae nested within the basal clades of Rhinotermitidae and was sister to Psammotermitinae. Overall, our phylogenetic results are largely consistent with earlier mitogenome studies., (© The Author(s) 2024. Published by Oxford University Press on behalf of Entomological Society of America.)

Subterranean termites in the genus Reticulitermes Holmgren 1913 are among the most economically important wood-destroying pests in the western United States. Yet, there remains uncertainty regarding the taxonomy and biology of the species in this genus. The 2 species described as having distributions in this region are the western subterranean termite, Reticulitermes hesperus Banks, and the arid land subterranean termite, Reticulitermes tibialis Banks. Taxonomic studies utilizing cuticular hydrocarbon (CHC) profiles, agonistic behavior, flight phenology, and mitochondrial DNA (mtDNA) suggested that R. hesperus is a species complex comprised of 2 or more sympatric, yet reproductively isolated species. To further delineate these taxa, we examined multiple genes from samples of Reticulitermes collected in the western United States. Alates collected after recent spring and fall mating flights, as well as previously collected workers, were subjected to CHC phenotyping and DNA sequence analyses that targeted mitochondrial cytochrome oxidase subunit II (COII), mitochondrial 16S rRNA, and nuclear Internal Transcribed Spacer 1 and 2 (ITS1 and 2). Phylogenetic analyses conducted also included published sequences of other putative western Reticulitermes species. Results suggest that at least 5 species of Reticulitermes may be present in California and that Reticulitermes in Arizona consistently group into multiple clades, including samples previously identified as R. tibialis in a sister clade. These analyses further support the species status of qualitatively different CHC phenotypes and that alates swarming in spring vs. fall are reproductively isolated species., (© The Author(s) 2023. Published by Oxford University Press on behalf of Entomological Society of America. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

(1) Unravelling the molecular basis underlying major evolutionary transitions can shed light on how complex phenotypes arise. The evolution of eusociality, a major evolutionary transition, has been demonstrated to be accompanied by enhanced gene regulation. Numerous pieces of evidence suggest the major impact of transposon insertion on gene regulation and its role in adaptive evolution. Transposons have been shown to be play a role in gene duplication involved in the eusocial transition in termites. However, evidence of the molecular basis underlying the eusocial transition in Blattodea remains scarce. Could transposons have facilitated the eusocial transition in termites through shifts of gene expression? (2) Using available cockroach and termite genomes and transcriptomes, we investigated if transposons insert more frequently in genes with differential expression in queens and workers and if those genes could be linked to specific functions essential for eusocial transition. (3) The insertion rate of transposons differs among differentially expressed genes and displays opposite trends between termites and cockroaches. The functions of termite transposon-rich queen- and worker-biased genes are related to reproduction and ageing and behaviour and gene expression, respectively. (4) Our study provides further evidence on the role of transposons in the evolution of eusociality, potentially through shifts in gene expression.

Present-day pterygote insects have two pairs of wings, one in the mesothorax (T2), the other in the metathorax (T3), and both have diverged in structure and function in different groups. Studies in endopterygote and paraneopteran species have shown that the gene Ultrabithorax ( Ubx ) specifies the identity and wing structure in T3, whereas the gene apterous ( ap ) significantly contributes to forming modified T2 wings. We wondered whether these Ubx and ap mechanisms operate in the lineage of polyneopterans. To explore this possibility, we used the cockroach Blattella germanica (Polyneoptera and Blattodea), in which the T2 wings are sclerotized (tegmina), whereas those of the T3 are membranous. We found that Ubx determines the structure of T3 and the membranous wing, while ap significantly contributes to form the sclerotized T2 tegmina. These results along with the studies carried out on the beetle Tribolium castaneum by Tomoyasu and collaborators suggest that ap plays an important role in the sclerotization and melanization of the T2 wings in neopteran groups that have sclerotized forewings. In turn, the sclerotizing properties of ap demonstrated in beetles and cockroaches suggest that the origin of this function goes back to the emergence of Neoptera, in the mid Devonian.

German cockroaches (Blattella germanica L.) harbor and disperse medically important pathogens and are a source of allergens that impact human health and wellbeing. Management of this pest requires an understanding of their distribution and dispersal. In this study, we collected German cockroaches from three apartment buildings in New Jersey, USA. We identified single-nucleotide polymorphisms (SNPs) from DNA extractions using next generation sequencing. We analyzed the SNPs and characterized cockroach population genetic structure using Fst, principal component, phylogenetic, and STRUCTURE analyses. We found significant differences in German cockroach population structure among the buildings. Within buildings, we found variable population structure that may be evidence for multiple colonization events. This study shows that SNPs derived from next generation sequencing provide a powerful tool for analyzing the genetic population structure of these medically important pests., (© The Author(s) 2022. Published by Oxford University Press on behalf of Entomological Society of America. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

In this paper, we report genome size (GS) values for nine cockroaches (order Blattodea, families Blattidae, Blaberidae and Ectobiidae, ex Blattelidae,), three of which are original additions to the ten already present in the GS database: the death's head roach (Blaberus craniifer), the Surinam cockroach (Pycnoscelus surinamensis) and the Madeira cockroach (Leucophaea maderae). Regarding the American cockroach (Periplaneta americana), the GS database contains two contrasting values (2.72 vs 3.41 pg); likely, the 2.72 pg value is the correct one as it is strikingly similar to our sperm DNA content evaluation (2.80 ± 0.11 pg). Also, we suggest halving the published GS of the Argentine cockroach Blaptica dubia and the spotted cockroach (the gray cockroach) Nauphoeta cinerea discussing (i) the occurrence of a correlation between increasing 2N chromosome number and GS within the order Blattodea; and (ii) the possible occurrence of a polyploidization phenomenon doubling a basic GS of 0.58 pg of some termite families (superfamily Blattoidea, epifamily Termitoidae).

Cockroaches, a sanitary pest, are essential species in insect developmental and metamorphic studies due to their easy feeding and hemimetabolous characteristics. Altogether with well-annotated genome sequences, these advantages have made American cockroach, Periplaneta americana, an important hemimetabolous insect model. Limited by the shortage of knockout strategy, effective RNA interference (RNAi)-based gene knockdown becomes an indispensable technique in functional gene research of P. americana. The present protocol describes the RNAi operation techniques in P. americana. The protocol includes (1) selection of the P. americana at proper developmental stages, (2) preparation for the injection setting, (3) dsRNA injection, and (4) gene knockdown efficiency detection. RNAi is a powerful reverse genetic tool in P. americana. The majority of P. americana tissues are sensitive to extracellular dsRNA. Its simplicity allows researchers to quickly obtain dysfunctional phenotypes under one or multiple targeting dsRNA injections, enabling researchers to better use the P. americana for developmental and metamorphic studies.

Subsocial Cryptocercus cockroaches are the sister group to termites and considered to be socially monogamous. Because genetic monogamy is a suggested requirement for evolution of cooperative breeding/eusociality, particularly in hymenopterans, clarification of the mating biology of Cryptocercus would help illuminate evolutionary trends in eusocial insects. To investigate possible extra-pair paternity in C. punctulatus, microsatellite markers were used to analyse offspring parentage, the stored sperm in females and results of experimental manipulation of sperm competition. Extra-pair paternity was common in field-collected families, but a lack of maternal alleles in several nymphs suggests sampling error or adoption. Isolating prereproductive pairs and assaying subsequently produced nymphs confirmed that nymphs lacked alleles from the pair male in 40% of families, with extra-pair male(s) siring 27%-77% of nymphs. Sperm of extra-pair males was detected in the spermatheca of 51% of paired prereproductive females. Mate switching and surgical manipulation of male mating ability indicated a tendency towards last male sperm precedence. Overall, the results demonstrate that about half of young females are serially monogamous during their maturational year, but bond, overwinter and produce their only set of offspring in company of the last mated male (=pair male). Repeated mating by the pair male increases the number of nymphs sired, but because many females use stored sperm of previous copulatory partners to fertilize eggs, pair males extend parental care to unrelated nymphs. The results suggest that genetic monogamy either developed in the termite ancestor after splitting from the Cryptocercus lineage, or that genetic monogamy may not be a strict prerequisite for the evolution of termite eusociality., (© 2021 John Wiley & Sons Ltd.)

Praying mantises are distributed all over the world. Though some Mantodea mitogenomes have been reported, an evolutionary genomic and phylogenetic analysis study lacks the latest taxonomic system. In the present study, four new mitogenomes were sequenced and annotated. Deroplatys truncate, D. lobate, Amorphoscelis chinensis and Macromantis sp. belong to Deroplatyidae, Amorphoscelidae and Photinaidae family, respectively. Our results indicated that the ATP8 gene may be lost in D. truncate and D. lobata mt genome, and four tRNA genes have not been found in D. truncate, D. lobata and Macromantis sp. A dN/dS pair analysis was conducted and it was found that all genes have evolved under purifying selection. Furthermore, we tested the phylogenetic relationships between the eight families of the Mantodea, including 35 species of praying Mantis. Based on the complete mitochondrial genome data, it was also suggested as sister to Deroplatyidae + Mantidae, Metallyticus sp., the only representative of Metallyticidae, is sister to the remaining mantises. Our results support the taxonomic system of Schwarz and Roy and are consistent with previous studies., Competing Interests: Co-Author Guy Smagghe is editor for the PLoS one journal and this does not alter the authors’ adherence to all PLoS one Policies on sharing data and materials.

Insects display an immense diversity in melanin pigmentation, which is generated by the interplay between the regulatory genes (that provide general patterning information) and effector genes (that provide coloration of the pattern). However, recent studies encompassing several different orders (Hemiptera, Blattodea, Coleoptera, and Lepidoptera) have shown that knockdowns of melanin producing genes alone can generate distinct region-specific patterns. This review surveys the most recent studies to further document the regional patterning of effector genes, and highlights the new advances and their implications for future research., (Copyright © 2021. Published by Elsevier Ltd.)

Two new species of the cockroach genus Eucorydia Hebard, 1929 from the Nansei Islands in Southwest Japan were compared to two closely related congeners, Eucorydia yasumatsui Asahina, 1971 and Eucorydia dasytoides (Walker, 1868). Eucorydia donanensis Yanagisawa, Sakamaki, and Shimano sp. nov. from Yonaguni-jima Island was characterized by an overall length of 12.5-14.5 mm in males. The dorsal side of the male abdomen was entirely dark purple and there was an obscure orange band running down the middle of the tegmen. Eucorydia tokaraensis Yanagisawa, Sakamaki, and Shimano sp. nov. was characterized by an overall length of 12.0-13.0 mm in males and a distinct orange band running down the middle of the tegmen. Eucorydia yasumatsui , E. donanensis , E. tokaraensis and the zonata population of E. dasytoides were divided into four lineages in a maximum-likelihood tree generated from a dataset concatenated from five (two nuclear, 28S rRNA, histone H3, and three mitochondrial, COII, 12S rRNA, 16S rRNA) genes. We recognized the three Japanese lineages E. yasumatsui , E. donanensis , and E. tokaraensis as distinct species, which were also supported by the pairwise genetic distances (5.4-7.8%, K2P) of the COI sequences. Morphometric analysis was performed on the genitalia. A principal component analysis plot revealed that the sizes of the genitalia in the three Japanese species were similar to each other and smaller than that of the zonata population of E. dasytoides . The analysis also revealed that the three Japanese species were distinguished from each other by combinations of the sizes of L3 and L7 sclerites and the shape of R2 sclerite, with some overlapping exceptions.

The development of new models to study diabetes in invertebrates is important to ensure adherence to the 3R's principle and to expedite knowledge of the complex molecular events underlying glucose toxicity. Streptozotocin (STZ)-an alkylating and highly toxic agent that has tropism to mammalian beta cells-is used as a model of type 1 diabetes in rodents, but little is known about STZ effects in insects. Here, the cockroach; Nauphoeta cinerea was used to determine the acute toxicity of 74 and 740 nmol of STZ injection per cockroach. STZ increased the glucose content, mRNA expression of glucose transporter 1 (GLUT1) and markers of oxidative stress in the head. Fat body glycogen, insect survival, acetylcholinesterase activity, triglyceride content and viable cells in head homogenate were reduced, which may indicate a disruption in glucose utilization by the head and fat body of insects after injection of 74 and 740 nmol STZ per nymph. The glutathione S-transferase (GST) activity and reduced glutathione levels (GSH) were increased, possibly via activation of nuclear factor erythroid 2 related factor as a compensatory response against the increase in reactive oxygen species. Our data present the potential for metabolic disruption in N. cinerea by glucose analogues and opens paths for the study of brain energy metabolism in insects. We further phylogenetically demonstrated conservation between N. cinerea glucose transporter 1 and the GLUT of other insects in the Neoptera infra-class.

Macrotermes barneyi, widely distributed in southern China, is the major fungus-growing termite in the subfamily Macrotermitinae. It has no flagellated protists in the guts. Here, we report occurrence of gregarine, a protozoan parasite in the digestive tract of M. barneyi. The general morphology and ultrastructure of the gregarine gamonts and syzygies by light micrograph and scanning electron micrograph are presented. SSU rDNA sequence analysis showed that the termite gregarine has the highest identity (90.10%) to that of Gregarina blattarum from cockroaches. Phylogenetic analysis based on the SSU rDNA sequences from diverse insect eugregarines indicated that the gregarine from M. barneyi is phylogenetically close to G. blattarus, L. erratica and G. tropica from Gregarinidae and Leidyanidae families, and may represent a novel species. This study expands our knowledge about the diversity of terrestrial eugregarines parasitizing in termites.

The subsocial life style and wood-feeding capability of Cryptocercus gives us an evolutionary key to unlock some outstanding questions in biology. With the advent of the Genomics Era, there is an unprecedented opportunity to address the evolution of eusociality and the acquisition of lignocellulases at the genetic level. However, to quantify gene expression, an appropriate normalization strategy is warranted to control for the non-specific variations among samples across different experimental conditions. To search for the internal references, 10 housekeeping genes from a gut transcriptome of a wood-feeding cockroach, Cryptocercus punctulatus, were selected as the candidates for the RT-qPCR analysis. The expression profiles of these candidates, including ACT, EF1α, GAPDH, HSP60, HSP70, αTUB, UBC, RPS18, ATPase and GST, were analyzed using a panel of analytical tools, including geNorm, NormFinder, BestKeeper, and comparative ΔC T method. RefFinder, a comprehensive ranking system integrating all four above-mentioned algorithms, rated ACT as the most stable reference gene for different developmental stages and tissues. Expression analysis of the target genes, Hex-1 and Cell-1, using the most or the least appropriate reference genes and a single or multiple normalizers signified this research. Our finding is the first step toward establishing a standardized RT-qPCR analysis in Cryptocercus.

Circadian clocks control rhythms in physiology and behavior entrained to 24 h light-dark cycles. Despite of conserved general schemes, molecular circadian clockworks differ between insect species. With RNA interference (RNAi) we examined an ancient circadian clockwork in a basic insect, the hemimetabolous Madeira cockroach Rhyparobia maderae. With injections of double-stranded RNA (dsRNA) of cockroach period (Rm´per), timeless 1 (Rm´tim1), or cryptochrome 2 (Rm´cry2) we searched for essential components of the clock´s core negative feedback loop. Single injections of dsRNA of each clock gene into adult cockroaches successfully and permanently knocked down respective mRNA levels within ~two weeks deleting daytime-dependent mRNA rhythms for Rm´per and Rm´cry2. Rm´perRNAi or Rm´cry2RNAi affected total mRNA levels of both genes, while Rm´tim1 transcription was independent of both, also keeping rhythmic expression. Unexpectedly, circadian locomotor activity of most cockroaches remained rhythmic for each clock gene knockdown employed. It expressed weakened rhythms and unchanged periods for Rm´perRNAi and shorter periods for Rm´tim1RNAi and Rm´cry2RNAi.As a hypothesis of the cockroach´s molecular clockwork, a basic network of switched differential equations was developed to model the oscillatory behavior of clock cells expressing respective clock genes. Data were consistent with two synchronized main groups of coupled oscillator cells, a leading (morning) oscillator, or a lagging (evening) oscillator that couple via mutual inhibition. The morning oscillators express shorter, the evening oscillators longer endogenous periods based on core feedback loops with either PER, TIM1, or CRY2/PER complexes as dominant negative feedback of the clockwork. We hypothesize that dominant morning oscillator cells with shorter periods express PER, but not CRY2, or TIM1 as suppressor of clock gene expression, while two groups of evening oscillator cells with longer periods either comprise TIM1 or CRY2/PER suppressing complexes. Modelling suggests that there is an additional negative feedback next to Rm´PER in cockroach morning oscillator cells., Competing Interests: The authors have declared that no competing interests exist.

Variations of speciation and extinction rates determine the fate of clades through time. Periods of high diversification and extinction (possibly mass-extinction events) can punctuate the evolutionary history of various clades, but they remain loosely defined for many biological groups, especially nonmarine invertebrates like insects. Here, we examine whether the cockroaches, mantises and termites (altogether included in Dictyoptera) have experienced episodic pulses of speciation or extinction and how these pulses may be associated with environmental fluctuations or mass extinctions. We relied on molecular phylogeny and fossil data to shed light on the times and rates at which dictyopterans diversified. The diversification of Dictyoptera has alternated between (i) periods of high diversification in the late Carboniferous, Early-Middle Triassic, Early Cretaceous and middle Palaeogene, and (ii) periods of high extinction rates particularly at the Permian-Triassic boundary, but not necessarily correlated with the major global biodiversity crises as in the mid-Cretaceous. This study advocates the importance of analyzing, when possible, both molecular phylogeny and fossil data to unveil diversification and extinction periods for a given group. The causes and consequences of extinction must be studied beyond mass-extinction events alone to gain a broader understanding of how clades wax and wane., (© The Willi Hennig Society 2020.)

Cold tolerance is considered an important factor determining the geographic distribution of insects. We have previously shown that despite its tropical origin, the cockroach Gromphadorinha coquereliana is capable of surviving exposures to cold. However, the freezing tolerance of this species had not yet been examined. Low temperature is known to alter membrane integrity in insects, but whether chilling or freezing compromises DNA integrity remains a matter of speculation. In the present study, we subjected the G. coquereliana adults to freezing to determine their supercooling point (SCP) and evaluated whether the cockroaches were capable of surviving partial and complete freezing. Next, we conducted single cell gel electrophoresis (SCGE) assays to determine whether heat, cold and freezing altered hemocyte DNA integrity. The SCP of this species was high and around -4.76°C, which is within the typical range of freezing-tolerant species. Most cockroaches survived to 1 day after partial ice formation (20% mortality), but died progressively in the next few days after cold stress (70% mortality after 4 days). One day after complete freezing, most insects died (70% mortality), and after 4 days, 90% of them had succumbed. The SCGE assays showed substantial levels of DNA damage in hemocytes. When cockroaches were heat-stressed, the level of DNA damage was similar to that observed in the freezing treatment, though all heat-stressed insects survived. The present study shows that G. coquereliana can be considered as moderately freeze-tolerant, and that extreme low temperature stress can affect DNA integrity, suggesting that this cockroach may possess an efficient DNA repair system., Competing Interests: Competing interestsThe authors declare that there are no conflicts of interest, financial or otherwise. The funder (National Science Center, Poland) had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript., (© 2019. Published by The Company of Biologists Ltd.)

The American cockroach ( Periplaneta americana ) is a medicinal insect. Its extract is used clinically to promote wound healing and tissue regeneration, but the effective medicinal components and mechanisms are not yet clear. It has been reported that human thymosin beta 4 (Tβ4) may accelerate skin wound healing, however, the role of P. americana thymosin (Pa-THYs) is still poorly understood. In the present study, we identify and analyze the DNA sequences of Pa-THYs by bioinformatics analysis. Then we clone, express, and purify the Pa-THYs proteins and evaluate the activity of recombinant Pa-THYs proteins by cell migration and proliferation assays in NIH/3T3 cells. To elucidate the role of Pa-THYs in wound healing, a mouse model is established, and we evaluate wound contraction, histopathological parameters, and the expressions of several key growth factors after Pa-THYs treatment. Our results showed that three THY variants were formed by skipping splicing of exons. Pa-THYs could promote fibroblast migration, but have no effect on fibroblast proliferation. In wound repair, Pa-THYs proteins could effectively promote wound healing through stimulating dermal tissue regeneration, angiogenesis, and collagen deposition. On the molecular mechanism, Pa-THYs also stimulated the expression of several key growth factors to promote wound healing. The data suggest that Pa-THYs could be a potential drug for promoting wound repair.

Insects are the only known animals in which sexual differentiation is controlled by sex-specific splicing. The doublesex transcription factor produces distinct male and female isoforms, which are both essential for sex-specific development. dsx splicing depends on transformer , which is also alternatively spliced such that functional Tra is only present in females. This pathway has evolved from an ancestral mechanism where dsx was independent of tra and expressed and required only in males. To reconstruct this transition, we examined three basal, hemimetabolous insect orders: Hemiptera, Phthiraptera, and Blattodea. We show that tra and dsx have distinct functions in these insects, reflecting different stages in the changeover from a transcription-based to a splicing-based mode of sexual differentiation. We propose that the canonical insect tra-dsx pathway evolved via merger between expanding dsx function (from males to both sexes) and narrowing tra function (from a general splicing factor to dedicated regulator of dsx )., Competing Interests: JW, ED, XB, CS, AW, MA, AK No competing interests declared, (© 2019, Wexler et al.)

In the endopterygote Drosophila melanogaster, Zelda is an activator of the zygotic genome during the maternal-to-zygotic transition (MZT). Zelda binds cis-regulatory elements (TAGteam heptamers), making chromatin accessible for gene transcription. Zelda has been studied in other endopterygotes: Apis mellifera and Tribolium castaneum, and the paraneopteran Rhodnius prolixus. We studied Zelda in the cockroach Blattella germanica, a hemimetabolan, short germ-band, and polyneopteran species. B. germanica Zelda has the complete set of functional domains, which is typical of species displaying ancestral features concerning embryogenesis. Interestingly, we found D. melanogaster TAGteam heptamers in the B. germanica genome. The canonical one, CAGGTAG, is present at a similar proportion in the genome of these two species and in the genome of other insects, suggesting that the genome admits as many CAGGTAG motifs as its length allows. Zelda-depleted embryos of B. germanica show defects involving blastoderm formation and abdomen development, and genes contributing to these processes are down-regulated. We conclude that in B. germanica, Zelda strictly activates the zygotic genome, within the MZT, a role conserved in more derived endopterygote insects. In B. germanica, zelda is expressed during MZT, whereas in D. melanogaster and T. castaneum it is expressed beyond this transition. In these species and A. mellifera, Zelda has functions even in postembryonic development. The expansion of zelda expression beyond the MZT in endopterygotes might be related with the evolutionary innovation of holometabolan metamorphosis. DATABASES: The RNA-seq datasets of B. germanica, D. melanogaster, and T. castaneum are accessible at the GEO databases GSE99785, GSE18068, GSE63770, and GSE84253. In addition, the RNA-seq library from T. castaneum adult females is available at SRA: SRX021963. The B. germanica reference genome is available as BioProject PRJNA203136., (© 2019 Federation of European Biochemical Societies.)

Termites are the principal decomposers in tropical and subtropical ecosystems around the world. Time-calibrated molecular phylogenies show that some lineages of Neoisoptera diversified during the Oligocene and Miocene, and acquired their pantropical distribution through transoceanic dispersal events, probably by rafting in wood. In this paper, we intend to resolve the historical biogeography of one of the earliest branching lineages of Neoisoptera, the Rhinotermitinae. We used the mitochondrial genomes of 27 species of Rhinotermitinae to build two robust time-calibrated phylogenetic trees that we used to reconstruct the ancestral distribution of the group. Our analyses support the monophyly of Rhinotermitinae and all genera of Rhinotermitinae. Our molecular clock trees provided time estimations that diverged by up to 15.6 million years depending on whether or not 3rd codon positions were included. Rhinotermitinae arose 50.4-64.6 Ma (41.7-74.5 Ma 95% HPD). We detected four disjunctions among biogeographic realms, the earliest of which occurred 41.0-56.6 Ma (33.0-65.8 Ma 95% HPD), and the latest of which occurred 20.3-34.2 Ma (15.9-40.4 Ma 95% HPD). These results show that the Rhinotermitinae acquired their distribution through a combination of transoceanic dispersals and dispersals across land bridges., (Copyright © 2018 Elsevier Inc. All rights reserved.)

The digestive system of selected phytophagous insects has been examined as a potential prospecting resource for identification of novel cellulolytic enzymes with potential industrial applications. In contrast to other model species, however, limited detailed information is available that characterizes cellulolytic activity and systems in basal hexapod groups. As part of a screening effort to identify insects with highly active cellulolytic systems, we have for the first time, identified species of Zygentoma that displayed the highest relative cellulase activity levels when compared to all other tested insect groups under the experimental conditions, including model species for cellulolytic systems such as termite and cockroach species in Rhinotermitidae (formerly Isoptera) and Cryptocercidae (formerly Blattodea). The goal of the present study was to provide a morphohistological characterization of cellulose digestion and to identify highly active cellulase enzymes present in digestive fluids of Zygentoma species. Morphohistological characterization supported no relevant differences in the digestive system of firebrat (Thermobia domestica) and the gray silverfish (Ctenolepisma longicaudata). Quantitative and qualitative cellulase assays identified the foregut as the region with the highest levels of cellulase activity in both T. domestica and C. longicaudata. However, T. domestica was found to have higher endoglucanase, xylanase and pectinase activities compared to C. longicaudata. Using nano liquid chromatography coupled to tandem mass spectrometry (nanoLC/MS/MS) and a custom gut transcriptome we identified cellulolytic enzymes from digestive fluids of T. domestica. Among the identified enzymes we report putative endoglucanases matching to insect or arthropod enzymes and glucan endo-1,6-β-glucosidases matching bacterial enzymes. These findings support combined activities of endogenous and symbiont-derived plant cell wall degrading enzymes in lignocellulose digestion in Zygentoma and advance our understanding of cellulose digestion in a primitive insect group., Competing Interests: The authors have declared that no competing interests exist.

Phylogenetic relationships among subgroups of cockroaches and termites are still matters of debate. Their divergence times and major phenotypic transitions during evolution are also not yet settled. We addressed these points by combining the first nuclear phylogenomic study of termites and cockroaches with a thorough approach to divergence time analysis, identification of endosymbionts, and reconstruction of ancestral morphological traits and behaviour. Analyses of the phylogenetic relationships within Blattodea robustly confirm previously uncertain hypotheses such as the sister-group relationship between Blaberoidea and remaining Blattodea, and Lamproblatta being the closest relative to the social and wood-feeding Cryptocercus and termites. Consequently, we propose new names for various clades in Blattodea: Cryptocercus + termites = Tutricablattae; Lamproblattidae + Tutricablattae = Kittrickea; and Blattoidea + Corydioidea = Solumblattodea. Our inferred divergence times contradict previous studies by showing that most subgroups of Blattodea evolved in the Cretaceous, reducing the gap between molecular estimates of divergence times and the fossil record. On a phenotypic level, the blattodean ground-plan is for egg packages to be laid directly in a hole while other forms of oviposition, including ovovivipary and vivipary, arose later. Finally, other changes in egg care strategy may have allowed for the adaptation of nest building and other novelties.

Objectives: The current study was conducted to identify cockroach species (Blattodea) of northwestern Iran in public places using morphological characteristics and ribosomal DNA internal transcribed spacer 2 (rDNA-ITS2). Sequences were analyzed with Basic Local Alignment Search Tool (BLAST) searches, Neighbor-Joining methods based on and Tamura-Nei phylogenetic analyses. In addition, eight cockroach rDNA-ITS2 sequences from China, India, Iran and the United States obtained from GenBank were compared to those obtained in this study., Results: Specimens collected in Iran were identified as Periplaneta americana (L.), Shelfordella lateralis (Walker), Blatta orientalis (L.) (Blattodea: Blattidae), Blattella germanica (L.), Supella longipalpa (F.) (Blattodea: Ectobiidae), Polyphaga aegyptiaca (L.), and Polyphaga saussurei (Dohrn) (Blattodea: Corydiidae). rDNA-ITS2 nucleotide sequence analysis showed 100% similarity between P. aegyptiaca and P. saussurei species collected from Iran despite morphological differences. However, ITS2 sequence of P. americana submitted from China showed 30.49-31.71% difference to P. americana sequences from Iran and the United States. The results highlight the importance of morphological identification of cockroach species before conducting molecular techniques.

Assessing support for molecular phylogenies is difficult because the data is heterogeneous in quality and overwhelming in quantity. Traditionally, node support values (bootstrap frequency, Bayesian posterior probability) are used to assess confidence in tree topologies. Other analyses to assess the quality of phylogenetic data (e.g. Lento plots, saturation plots, trait consistency) and the resulting phylogenetic trees (e.g. internode certainty, parameter permutation tests, topological tests) exist but are rarely applied. Here we argue that a single qualitative analysis is insufficient to assess support of a phylogenetic hypothesis and relate data quality to tree quality. We use six molecular markers to infer the phylogeny of Blattodea and apply various tests to assess relationship support, locus quality, and the relationship between the two. We use internode-certainty calculations in conjunction with bootstrap scores, alignment permutations, and an approximately unbiased (AU) test to assess if the molecular data unambiguously support the phylogenetic relationships found. Our results show higher support for the position of Lamproblattidae, high support for the termite phylogeny, and low support for the position of Anaplectidae, Corydioidea and phylogeny of Blaberoidea. We use Lento plots in conjunction with mutation-saturation plots, calculations of locus homoplasy to assess locus quality, identify long branch attraction, and decide if the tree's relationships are the result of data biases. We conclude that multiple tests and metrics need to be taken into account to assess tree support and data robustness., (Copyright © 2018 Elsevier Inc. All rights reserved.)

Sterile castes are a defining criterion of eusociality; investigating their evolutionary origins can critically advance theory. In termites, the soldier caste is regarded as the first acquired permanently sterile caste. Previous studies showed that juvenile hormone (JH) is the primary factor inducing soldier differentiation, and treatment of workers with artificial JH can generate presoldier differentiation. It follows that a shift from a typical hemimetabolous JH response might be required for soldier formation during the course of termite evolution within the cockroach clade. To address this possibility, analysis of the role of JH and its signaling pathway was performed in the termite Zootermopsis nevadensis and compared with the wood roach Cryptocercus punctulatus , a member of the sister group of termites. Treatment with a JH analog (JHA) induced a nymphal molt in C. punctulatus RNA interference (RNAi) of JH receptor Methoprene tolerant ( Met ) was then performed, and it inhibited the presoldier molt in Z. nevadensis and the nymphal molt in C. punctulatus Knockdown of Met in both species inhibited expression of 20-hydroxyecdysone (20E; the active form of ecdysone) synthesis genes. However, in Z. nevadensis , several 20E signaling genes were specifically inhibited by Met RNAi. Consequently, RNAi of these genes were performed in JHA-treated termite individuals. Knockdown of 20E signaling and nuclear receptor gene, Hormone receptor 39 ( HR39 / FTZ-F1β ) resulted in newly molted individuals with normal worker phenotypes. This is the first report of the JH-Met signaling feature in termites and Cryptocercus JH-dependent molting activation is shared by both taxa and mediation between JH receptor and 20E signalings for soldier morphogenesis is specific to termites., (Copyright © 2018 by the Genetics Society of America.)

Social insects show an extreme degree of phenotypic plasticity. In highly eusocial species, this manifests in the generation of distinct castes with extreme differences in both morphology and life span. The molecular basis of these differences is highly entangled and not fully understood, but several recent studies demonstrated that insulin/insulin-like growth factor signaling (IIS) is one of the key pathways. Here, we investigate the molecular evolution of insect insulin receptors (InRs), which are membrane-bound dimers that enable IIS by relaying extracellular signals to intracellular signaling cascades. Classic models of invertebrate IIS include only one InR gene, but some recent studies on less commonly studied insects have found two InRs, which act in an antagonistic manner to facilitate polyphenism in at least one documented case. We search 22 arthropod genomes and identify several InR copies and their evolutionary origin that were lacking from previous annotations. Phylogenetic analysis shows that the two insect InR genes date back at least 400 million years to a common ancestor of winged insects. Most notably, we also identified the evolutionary origin of a third InR copy that is unique to the clade of Blattodea, just before therein the eusocial termites evolved. One of the InR paralogs consistently shows caste-biased expression in all three termites, which strongly suggests a role in caste differentiation. These results have important ramifications for past and future InR inhibition/InR knockdown experiments in insects and they provide a set of key genes regulating life span and morphology in termite castes., (© 2018 Wiley Periodicals, Inc.)

Almost all examined cockroaches harbor an obligate intracellular endosymbiont, Blattabacterium cuenoti. On the basis of genome content, Blattabacterium has been inferred to recycle nitrogen wastes and provide amino acids and cofactors for its hosts. Most Blattabacterium strains sequenced to date harbor a genome of ∼630 kbp, with the exception of the termite Mastotermes darwiniensis (∼590 kbp) and Cryptocercus punctulatus (∼614 kbp), a representative of the sister group of termites. Such genome reduction may have led to the ultimate loss of Blattabacterium in all termites other than Mastotermes. In this study, we sequenced 11 new Blattabacterium genomes from three species of Cryptocercus in order to shed light on the genomic evolution of Blattabacterium in termites and Cryptocercus. All genomes of Cryptocercus-derived Blattabacterium genomes were reduced (∼614 kbp), except for that associated with Cryptocercus kyebangensis, which comprised 637 kbp. Phylogenetic analysis of these genomes and their content indicates that Blattabacterium experienced parallel genome reduction in Mastotermes and Cryptocercus, possibly due to similar selective forces. We found evidence of ongoing genome reduction in Blattabacterium from three lineages of the C. punctulatus species complex, which independently lost one cysteine biosynthetic gene. We also sequenced the genome of the Blattabacterium associated with Salganea taiwanensis, a subsocial xylophagous cockroach that does not vertically transmit gut symbionts via proctodeal trophallaxis. This genome was 632 kbp, typical of that of nonsubsocial cockroaches. Overall, our results show that genome reduction occurred on multiple occasions in Blattabacterium, and is still ongoing, possibly because of new associations with gut symbionts in some lineages.

We acquired 21 complete mitogenomes and 6 nearly complete mitogenomes of mantises belonging to 8 families (Hymenopodidae, Iridopterygidae, Mantidae, Metallyticidae, Sibyllidae, Tarachodidae, Thespidae, Toxoderidae) using 14 pairs of mantid specific primer sets and found that 5 species of mantises have duplicate copies (2-4) of trnR: Ambivia undata, Creobroter jiangxiensis, Creobroter urbanus, Phyllothelys sp1. and Theopropus elegans while two novel gene arrangements CR-I-NCR-I-NCR-I-NCR-I-NCR-I-NCR*-Q-M and COII-K*-D*-K-D*-K-D*-K-D were found in Schizocephala bicornis and Stenotoxodera porioni, respectively. The multiple copies of trnR are caused by independent duplications. The gene arrangements in Stenotoxodera porioni with three identical copies of trnK can be explained as mid-way through the TDRL process while the form of gene arrangement in Schizocephala bicornis is unclear. In the phylogeny at the family level, the monophyly of Liturgusidae and Iridopterygidae was supported, whereas the monophyly of Hymenopodidae, Mantidae and Tarachodidae wasn't. The features of mantis mitochondrial genomes including high duplication rates of trnR, trnK and trnI indicate that Mantodea mitochondrial genomes maybe a useful model system for studying gene duplication. However, derived gene arrangements may not be appropriate for phylogenetic inference in Mantodea as they aren't synapomorphy and aren't shared by close relatives., (Copyright © 2018 Elsevier B.V. All rights reserved.)

RNA interference (RNAi) has been widely applied for uncovering the biological functions of numerous genes, and has been envisaged as a pest control tool operating by disruption of essential gene expression. Although different methods, such as injection, feeding, and soaking, have been reported for successful delivery of double-stranded RNA (dsRNA), the efficiency of RNAi through oral delivery of dsRNA is highly variable among different insect groups. The German cockroach, Blattella germanica, is highly sensitive to the injection of dsRNA, as shown by many studies published previously. The present study describes a method to demonstrate that the dsRNA encapsulated with liposome carriers is sufficient to retard the degradation of dsRNA by midgut juice. Notably, the continuous feeding of dsRNA encapsulated by liposomes significantly reduces the tubulin expression in the midgut, and led to the death of cockroaches. In conclusion, the formulation and utilization of dsRNA lipoplexes, which protect dsRNA against nucleases, could be a practical use of RNAi for insect pest control in the future.

Following the acceptance of plate tectonics theory in the latter half of the 20th century, vicariance became the dominant explanation for the distributions of many plant and animal groups. In recent years, however, molecular-clock analyses have challenged a number of well-accepted hypotheses of vicariance. As a widespread group of insects with a fossil record dating back 300 My, cockroaches provide an ideal model for testing hypotheses of vicariance through plate tectonics versus transoceanic dispersal. However, their evolutionary history remains poorly understood, in part due to unresolved relationships among the nine recognized families. Here, we present a phylogenetic estimate of all extant cockroach families, as well as a timescale for their evolution, based on the complete mitochondrial genomes of 119 cockroach species. Divergence dating analyses indicated that the last common ancestor of all extant cockroaches appeared ∼235 Ma, ∼95 My prior to the appearance of fossils that can be assigned to extant families, and before the breakup of Pangaea began. We reconstructed the geographic ranges of ancestral cockroaches and found tentative support for vicariance through plate tectonics within and between several major lineages. We also found evidence of transoceanic dispersal in lineages found across the Australian, Indo-Malayan, African, and Madagascan regions. Our analyses provide evidence that both vicariance and dispersal have played important roles in shaping the distribution and diversity of these insects.

In this study, the complete mitochondrial genome of Cryptocercus meridianus was sequenced. The circular mitochondrial genome is 15,322 bp in size and contains 13 protein-coding genes, two ribosomal RNA genes (12S rRNA and 16S rRNA), 22 transfer RNA genes, and one D-loop region. We compare the mitogenome of C. meridianus with that of C. relictus and C. kyebangensis . The base composition of the whole genome was 45.20%, 9.74%, 16.06%, and 29.00% for A, G, C, and T, respectively; it shows a high AT content (74.2%), similar to the mitogenomes of C. relictus and C. kyebangensis . The protein-coding genes are initiated with typical mitochondrial start codons except for cox1 with TTG. The gene order of the C. meridianus mitogenome differs from the typical insect pattern for the translocation of tRNA-Ser AGN , while the mitogenomes of the other two Cryptocercus species, C. relictus and C. kyebangensis , are consistent with the typical insect pattern. There are two very long non-coding intergenic regions lying on both sides of the rearranged gene tRNA-Ser AGN . The phylogenetic relationships were constructed based on the nucleotide sequence of 13 protein-coding genes and two ribosomal RNA genes. The mitogenome of C. meridianus is the first representative of the order Blattodea that demonstrates rearrangement, and it will contribute to the further study of the phylogeny and evolution of the genus Cryptocercus and related taxa., Competing Interests: The authors declare no conflict of interest.

Our current taxonomic perspective on Entamoeba is largely based on small-subunit ribosomal RNA genes (SSU rDNA) from Entamoeba species identified in vertebrate hosts with minor exceptions such as E. moshkovskii from sewage water and E. marina from marine sediment. Other Entamoeba species have also been morphologically identified and described from non-vertebrate species such as insects; however, their genetic diversity remains unknown. In order to further disclose the diversity of the genus, we investigated Entamoeba spp. in the intestines of three cockroach species: Periplaneta americana, Blaptica dubia, and Gromphadorhina oblongonota. We obtained 134 Entamoeba SSU rDNA sequences from 186 cockroaches by direct nested PCR using the DNA extracts of intestines from cockroaches, followed by scrutinized BLASTn screening and phylogenetic analyses. All the sequences identified in this study were distinct from those reported from known Entamoeba species, and considered as novel Entamoeba ribosomal lineages. Furthermore, they were positioned at the base of the clade of known Entamoeba species and displayed remarkable degree of genetic diversity comprising nine major groups in the three cockroach species. This is the first report of the diversity of SSU rDNA sequences from Entamoeba in non-vertebrate host species, and should help to understand the genetic diversity of the genus Entamoeba.

Cockroach allergens can lead to serious allergy and asthma symptoms. Termites are evolutionarily related to cockroaches, cohabitate in human dwellings, and represent an increasing pest problem in the United States. The Formosan subterranean termite (Coptotermes formosanus) is one of the most common species in the southern United States. Several assays were used to determine if C. formosanus termite proteins cross-react with cockroach allergens. Expressed sequence tag and genomic sequencing results were searched for homology to cockroach allergens using BLAST 2.2.21 software. Whole termite extracts were analyzed by mass-spectrometry, immunoassay with IgG and scFv antibodies to cockroach allergens, and human IgE from serum samples of cockroach allergic patients. Expressed sequence tag and genomic sequencing results indicate greater than 60% similarity between predicted termite proteins and German and American cockroach allergens, including Bla g 2/Per a 2, Bla g 3/Per a 3, Bla g 5, Bla g 6/Per a 6, Bla g 7/Per a 7, Bla g 8, Per a 9, and Per a 10. Peptides from whole termite extract were matched to those of the tropomyosin (Bla g 7), arginine kinase (Per a 9), and myosin (Bla g 8) cockroach allergens by mass-spectrometry. Immunoblot and ELISA testing revealed cross-reaction between several proteins with IgG and IgE antibodies to cockroach allergens. Several termite proteins, including the hemocyanin and tropomyosin orthologs of Blag 3 and Bla g 7, were shown to crossreact with cockroach allergens. This work presents support for the hypothesis that termite proteins may act as allergens and the findings could be applied to future allergen characterization, epitope analysis, and clinical studies.

Complete mitochondrial genomes (mitogenomes) of two cockroach species, Periplaneta australasiae and Neostylopyga rhombifolia, 15,605 bp and 15,711 bp in length, respectively, were determined. As reported for other cockroach mitogenomes, the two mitogenomes possessed typical ancestral insect mitogenome gene composition and arrangement. Only several small intergenic spacers were found: one, which was common in all sequenced cockroach mitogenomes except for the genus Cryptocercus, was between tRNA-Ser (UCN) and ND1 and contained a 7bp highly conserved motif (WACTTAA). Three different types of short tandem repeats in the N. rhombifolia control region (CR) were observed. The homologous alignments of these tandem repeats with other six cockroach mitogenome CRs revealed a low similarity. Three conserved sequence blocks (CSB) were detected in both cockroach mitochondrial CRs. CSB1 was specific for blattinine mitogenomes and was highly conserved with 95% similarity, speculating that this block was a possible molecular synapomorphy for this subfamily. CSB3 located nearby downstream of CSB1 and has more variations within blattinine mitogenomes compared with CSB1. The CSB3 was capable of forming stable stem-loop structure with a small T-stretch in the loop portion. We assessed the influence of four datasets and two inference methods on topology within Orthopteroidea. All genes excluding the third codon positions of PCGs could generate more stable topology, and higher posterior probabilities than bootstrap values were presented at some branch nodes. The phylogenetic analysis with different datasets and analytical methods supported the monophyly of Dictyoptera and supported strongly the proposal that Isoptera should be classified as a family (Termitidae) of the Blattaria. Specifically, Shelfordella lateralis was inserted in the clade Periplaneta. Considering the K2P genetic distance, morphological characters, and the phylogenetic trees, we suggested that S. lateralis should be placed in the genus Periplaneta.

Group-level properties, such as collective movements or decisions, can be considered an outcome of the interplay between individual behavior and social interactions. However, the respective influences of individual preferences and social interactions are not evident. In this research, we study the implications of behavioral variability on the migration dynamics of a group of gregarious insects (Periplaneta americana) subjected to two different disturbance regimes (one without disturbances and another one with high frequency of disturbances). The results indicate that individuals presented consistent behavior during the nighttime (active phase of cockroaches) in both conditions. Moreover, we used a modeling approach to test the role of personality during the migration process. The model considers identical individuals (no personality) without memory and no direct inter-attraction between individuals. The agreement between theoretical and experimental results shows that behavioral variability play a secondary role during migration dynamics. Our results showing individual personality during the nighttime (spontaneous decision to forage) but not during the emigration process (induced by environmental disturbances) highlight the plasticity of personality traits.

During the period between 1999 and 2006, wood-feeding cockroaches in the Cryptocercus punctulatus Scudder species complex were collected throughout Great Smoky Mountains National Park, USA. The chromosome numbers of insects from 59 sites were determined, and phylogenetic analyses were performed based on mitochondrial COII and nuclear ITS2 DNA. The distribution of the three male karyotypes found in the park (2n = 37, 39, and 45) is mapped and discussed in relation to recent disturbances and glacial history. Clades of the three karyotype groups meet near the ridgeline separating North Carolina from Tennessee in the center of the park, suggesting that these may have originated from separate lower elevation refugia after the last glacial maximum. The timing of divergence and a significant correlation between elevation difference and genetic distance in two of the clades supports this hypothesis. The ecological role of the cockroaches in the park is discussed., (© The Authors 2017. Published by Oxford University Press on behalf of Entomological Society of America.)

A β-glucosidase (BG), PaBG1b, from the xylophagous cockroach Panesthia angustipennis spadica was heterologously expressed in the methylotrophic yeast Pichia pastoris, purified, and biochemically characterized. Post-translational modification and N-terminal sequencing analysis demonstrated that the expression product was comprised of two polypeptides with different N-terminal sequences, presumably due to the presence of lysine-arginine (KR) sequence in the putative mature region. Substrate specificity analysis showed that PaBG1b hydrolyzed a broad range of substrates including cellohexaose, with the preference for aryl β-d-fucosyl linkage and laminaribiose. Although the glucose tolerance of PaBG1b was moderate (K i =200.3±1.1mM), PaBG1b demonstrated high specific activity and catalytic efficiency towards cellobiose with V max and k cat /K m values of 436.7±6.3U/mg and 109.8mM -1 s -1 , respectively. In addition, PaBG1b was not inhibited by cellobiose up to the highest concentration tested (100mM). Collectively, our work demonstrates that PaBG1b is a potentially valuable BG for commercial bioethanol production from cellulose., (Copyright © 2016 Elsevier Inc. All rights reserved.)

The first complete mitochondrial genome in Blaberidae has been reported in this research. Opisthoplatia orientalis (Blaberidae, Epilamprinae) known as ground cockroach with golden edge is distributed in East and South Asia and widely used for thrombolytic therapy in China. Meanwhile, Blaptica dubia (Blaberidae, Blaberinae) has been adopted as feeder insect for various kinds of pets all over the world. In the present study, we investigated the complete mitochondrial genome of O. orientalis and B. dubia, and the mitogenome is 18 724  and 17 340 bp in length, respectively. The circular molecule consists of 13 protein coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and a non-coding control region, with an AT content of 75.7% for O. orientalis and 72.8% for B. dubia. A preliminary phylogenetic analysis has been carried out with 11 related species and the status of these two species is further confirmed.

We identified the antifungal gene termicin in three species of Cryptocercus woodroaches. Cryptocercus represents the closest living cockroach lineage of termites, which suggests that the antifungal role of termicin evolved prior to the divergence of termites from other cockroaches. An analysis of Cryptocercus termicin and two β-1,3-glucanase genes (GNBP1 and GNBP2), which appear to work synergistically with termicin in termites, revealed evidence of selection in these proteins. We identified the signature of past selective sweeps within GNBP2 from Cryptocercus punctulatus and Cryptocercus wrighti. The signature of past selective sweeps was also found within termicin from Cryptocercus punctulatus and Cryptocercus darwini. Our analysis further suggests a phenotypically identical variant of GNBP2 was maintained within Cryptocercus punctulatus, Cryptocercus wrighti, and Cryptocercus darwini while synonymous sites diverged. Cryptocercus termicin and GNBP2 appear to have experienced similar selective pressure to that of their termite orthologues in Reticulitermes. This selective pressure may be a result of ubiquitous entomopathogenic fungal pathogens such as Metarhizium. This study further reveals the similarities between Cryptocercus woodroaches and termites.