Evolution and phylogenetic information content of the ribosomal DNA repeat unit in the Blattodea (Insecta)

Solving the 250-year-old mystery of the origin and global spread of the German cockroach, Blattella germanica

The origin of the German cockroach, Blattella germanica, is enigmatic, in part because it is ubiquitous worldwide in human-built structures but absent from any natural habitats. The first historical records of this species are from ca. 250 years ago (ya) from central Europe (hence its name). However, recent research suggests that the center of diversity of the genus is Asian, where its closest relatives are found. To solve this paradox, we sampled genome-wide markers of 281 cockroaches from 17 countries across six continents. We confirm that B. germanica evolved from the Asian cockroach Blattella asahinai approximately 2,100 ya, probably by adapting to human settlements in India or Myanmar. Our genomic analyses reconstructed two primary global spread routes, one older, westward route to the Middle East coinciding with various Islamic dynasties (~1,200 ya), and another younger eastward route coinciding with the European colonial period (~390 ya). While Europe was not central to the early domestication and spread of the German cockroach, European advances in long-distance transportation and temperature-controlled housing were likely important for the more recent global spread, increasing chances of successful dispersal to and establishment in new regions. The global genetic structure of German cockroaches further supports our model, as it generally aligns with geopolitical boundaries, suggesting regional bridgehead populations established following the advent of international commerce.

R2 and Non-Site-Specific R2-Like Retrotransposons of the German Cockroach, Blattella germanica

The structural and functional organization of the ribosomal RNA gene cluster and the full-length R2 non-LTR retrotransposon (integrated into a specific site of 28S ribosomal RNA genes) of the German cockroach, Blattella germanica, is described. A partial sequence of the R2 retrotransposon of the cockroach Rhyparobia maderae is also analyzed. The analysis of previously published next-generation sequencing data from the B. germanica genome reveals a new type of retrotransposon closely related to R2 retrotransposons but with a random distribution in the genome. Phylogenetic analysis reveals that these newly described retrotransposons form a separate clade. It is shown that proteins corresponding to the open reading frames of newly described retrotransposons exhibit unequal structural domains. Within these retrotransposons, a recombination event is described. New mechanism of transposition activity is discussed. The essential structural features of R2 retrotransposons are conserved in cockroaches and are typical of previously described R2 retrotransposons. However, the investigation of the number and frequency of 5′-truncated R2 retrotransposon insertion variants in eight B. germanica populations suggests recent mobile element activity. It is shown that the pattern of 5′-truncated R2 retrotransposon copies can be an informative molecular genetic marker for revealing genetic distances between insect populations.

Ribosomal DNA internal transcribed spacer 2 sequence analysis and phylogenetic comparison of seven cockroach species in northwestern Iran

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.

Intragenomic variations of multicopy ITS2 marker in Agrodiaetus blue butterflies (Lepidoptera, Lycaenidae)

The eukaryotic ribosomal DNA cluster consists of multiple copies of three genes, 18S, 5. 8S and 28S rRNAs, separated by multiple copies of two internal transcribed spacers, ITS1 and ITS2. It is an important, frequently used marker in both molecular cytogenetic and molecular phylogenetic studies. Despite this, little is known about intragenomic variations within the copies of eukaryotic ribosomal DNA genes and spacers. Here we present data on intraindividual variations of ITS2 spacer in three species of Agrodiaetus Hübner, 1822 blue butterflies revealed by cloning technique. We demonstrate that a distinctly different intragenomic ITS2 pattern exists for every individual analysed. ITS2 sequences of these species show significant intragenomic variation (up to 3.68% divergence), setting them apart from each other on inferred phylogenetic tree. This variation is enough to obscure phylogenetic relationships at the species level.

Phylogenetic information content of Copepoda ribosomal DNA repeat units: ITS1 and ITS2 impact

The utility of various regions of the ribosomal repeat unit for phylogenetic analysis was examined in 16 species representing four families, nine genera, and two orders of the subclass Copepoda (Crustacea). Fragments approximately 2000 bp in length containing the ribosomal DNA (rDNA) 18S and 28S gene fragments, the 5.8S gene, and the internal transcribed spacer regions I and II (ITS1 and ITS2) were amplified and analyzed. The DAMBE (Data Analysis in Molecular Biology and Evolution) software was used to analyze the saturation of nucleotide substitutions; this test revealed the suitability of both the 28S gene fragment and the ITS1/ITS2 rDNA regions for the reconstruction of phylogenetic trees. Distance (minimum evolution) and probabilistic (maximum likelihood, Bayesian) analyses of the data revealed that the 28S rDNA and the ITS1 and ITS2 regions are informative markers for inferring phylogenetic relationships among families of copepods and within the Cyclopidae family and associated genera. Split-graph analysis of concatenated ITS1/ITS2 rDNA regions of cyclopoid copepods suggested that the Mesocyclops, Thermocyclops, and Macrocyclops genera share complex evolutionary relationships. This study revealed that the ITS1 and ITS2 regions potentially represent different phylogenetic signals.

Molecular adaptation and resilience of the insect's nuclear receptor USP

Background

The maintenance of biological systems requires plasticity and robustness. The function of the ecdysone receptor, a heterodimer composed of the nuclear receptors ECR (NR1H1) and USP (NR2B4), was maintained in insects despite a dramatic divergence that occurred during the emergence of Mecopterida. This receptor is therefore a good model to study the evolution of plasticity. We tested the hypothesis that selection has shaped the Ligand-Binding Domain (LBD) of USP during evolution of Mecopterida.

Results

We isolated usp and cox1 in several species of Drosophilidae, Tenebrionidae and Blattaria and estimated non-synonymous/synonymous rate ratios using maximum-likelihood methods and codon-based substitution models. Although the usp sequences were mainly under negative selection, we detected relaxation at residues located on the surface of the LBD within Mecopterida families. Using branch-site models, we also detected changes in selective constraints along three successive branches of the Mecopterida evolution. Residues located at the bottom of the ligand-binding pocket (LBP) underwent strong positive selection during the emergence of Mecopterida. This change is correlated with the acquisition of a large LBP filled by phospholipids that probably allowed the stabilisation of the new Mecopterida structure. Later, when the two subgroups of Mecopterida (Amphiesmenoptera: Lepidoptera, Trichoptera; Antliophora: Diptera, Mecoptera, Siphonaptera) diverged, the same positions became under purifying selection. Similarly, several positions of the heterodimerisation interface experienced positive selection during the emergence of Mecopterida, rapidly followed by a phase of constrained evolution. An enlargement of the heterodimerisation surface is specific for Mecopterida and was associated with a reinforcement of the obligatory partnership between ECR and USP, at the expense of homodimerisation.

Conclusions

In order to explain the episodic mode of evolution of USP, we propose a model in which the molecular adaptation of this protein is seen as a process of resilience for the maintenance of the ecdysone receptor functionality.

Structure and molecular evolution of the ribosomal DNA external transcribed spacer in the cockroach genus Blattella

The ribosomal DNA (rDNA) cluster of insects contains several hundred repeating structural-functional units and, therefore, is a typical example of a multigene family. Eukaryotic ribosomal RNA (rRNA) genes (18S, 5.8S, and 28S like) are arranged in tandemly repeated clusters in the nucleolus organizers, separated by several spacers, namely the nontranscribed spacer, the external transcribed spacer (ETS), and the internal transcribed spacers. The nucleotide sequences of the ETS of the three closely related Blattella cockroach species, Blattella germanica (Linnaeus, 1767), Blattella asahinai (Mizukubo, 1981), and Blattella lituricollis (Walker, 1868), were determined and compared. The three species had relatively similar ETS lengths, and sequence differences among them could be explained by two types of rearrangements, namely deletions of subrepeats and nucleotide substitutions. Minor ETS variants in B. germanica differed from the major variant in the same way that the major ETS variants of the three Blattella species differed from each other. Concerted evolution and the birth-and-death models, which are often invoked to explain the diversity and evolution of the multigene families of rDNA clusters, are discussed in the light of our data. A new model is proposed to explain the evolutionary reorganization of the ETS region: evolution of rDNA by "magnification-and-fixation" is characterized by magnification of minor subrepeats, which become adaptive in a new rapidly changed environment, and subsequent fixation of this variant type as a major component of the multigene family of a new species.

Characterization of internal transcribed spacer (ITS1)-ITS2 region of ribosomal RNA gene from 25 species of Culicoides biting midges (Diptera: Ceratopogonidae) in Japan

We determined nucleotide sequences of the nuclear rDNA internal transcribed spacer (ITS)1-5.8S-ITS2a-2S-ITS2 region in 103 individuals of 25 Culicoides species (Diptera: Ceratopogonidae) from 11 locations in Japan. Ribosomal RNA genes, 5.8S and 2S rDNA, were highly conserved among the species with few variations. The ITS2a region showed length variation among species. Both ITS1 and ITS2 showed highly varied sequences among species. The noticeable indel regions among ITS1 sequences are present in some Culicoides species, separating species into two types having long or short ITS1 region. However, Culicoides cylindratus Kitaoka possesses both types of ITS1 in each individual; these results seem to indicate that the ITS1-long type was the prototype and the short type was produced through deletion in many Culicoides species. One species, belonging to subgenus Avaritia, possessed an Avaritia-specific sequence in ITS1 and phylogenetically formed a monophyletic group. Geographical genotypes in a species were not clear. Species-specific sequence features were observed, enabling molecular identification of Culicoides species.

Comparative analysis of sequences and secondary structures of the rRNA internal transcribed spacer 2 (ITS2) in pollen beetles of the subfamily Meligethinae (Coleoptera, Nitidulidae): potential use of slippage-derived sequences in molecular systematics

A comparative analysis of ITS2 sequences and secondary structures in 89 species of pollen beetles of the subfamily Meligethinae (Coleoptera, Nitidulidae) was performed. The ITS2 folding pattern was highly conserved and comparable with the general model proposed for eukaryotes. Simple sequence repeats (SSRs) were responsible for most of the observed nucleotide variability (approximately 1-3%) and length variation (359-459bp). When plotted on secondary structures, SSRs mapped in expansion segments positioned at the apices of three ITS2 helices ('A', 'B' and 'D1') and appeared to have evolved under mechanisms of compensatory slippage. Homologies among SSRs nucleotides could not be unambiguously assigned, and thus were not useful to resolve phylogeny. However, slippage-derived motifs provided some preliminary genetic support for newly proposed taxonomic arrangements of several genera and subgenera of Meligethinae, corroborating existing morphological and ecological datasets.

Intraspecific variation and population structure of the German cockroach, Blattella germanica, revealed with RFLP analysis of the non-transcribed spacer region of ribosomal DNA

Little information is available on genetic variation within and between populations of pest cockroaches. In this study, intraspecific HindIII polymorphism was investigated in the German cockroach, Blattella germanica (Linnaeus) (Dictyoptera, Blattaria: Blattellidae), using restriction fragment length polymorphisms (RFLP) of the non-transcribed spacer (NTS) region of ribosomal DNA (rDNA). Individual male insects were collected from infestations at three different pig farms. Each population was characterized by HindIII restriction fragment frequencies and haplotype (a particular X-chromosome pattern) frequencies. The inheritance of the X-chromosome HindIII rDNA patterns over 12 generations (3 years) follows Mendelian patterns, and the stability of this polymorphic marker indicates infrequent genetic recombination of variable sites. Although pairwise genetic distance measures were uncorrelated with geographical distance, the pattern of genetic differentiation of the three cockroach populations suggests that human-mediated transport of cockroaches is an important force in shaping the population genetic structure of cockroach infestations, at least at the regional scale of 10-100 km. Sequence variation in the ribosomal NTS is a useful marker, and RFLP of rDNA is a simple, robust and reproducible technique for differentiating recently diverged cockroach populations.

Intragenomic heterogeneity of internal transcribed spacer rDNA in neotropical malaria vector Anopheles aquasalis (Diptera: Culicidae)

Intragenomic heterogeneity of the internal transcribed spacer (ITS) array was investigated in Anopheles aquasalis Curry mosquitoes from two geographic locations in each of Brazil and Venezuela, and one in Suriname. Polymerase chain reaction-amplified copies of the ITS were cloned and sequenced. The length of the entire array ranged from 782 to 990 bp, with most variation due to microsatellite insertions in ITS1. We detected 40 different ITSL sequences and 15 different ITS2 sequences of the 71 to 72 clones examined. The sequence divergence within localities ranged from 0.002 to 0.043 for ITS1 and from 0 to 0.006 for ITS2. Point mutations were common to both spacer regions, but dinucleotide microsatellite repeats were restricted to ITS1. Sequences from neither ITS1 nor ITS2 had a diagnostic distribution or were informative in distinguishing these populations, providing additional support for the status of An. aquasalis as a single species.

Comparison of oligosaccharide compositions in male nuptial secretions of three cockroach species of the genus Blattella

Male cockroaches of the genus Blattella secrete courtship pheromones from abdominal tergal glands. Conspecific females mount courting males and feed on their tergal exudates. In the German cockroach, B. germanica, a complex mixture of oligosaccharides was previously identified as the major phagostimulants in the male tergal secretion. We compared the oligosaccharide compositions of three Blattella species from different habitats, B. germanica (house), B. lituricollis (field), and B. nipponica (forest), by gas chromatographic (GC) analysis of their trifluoroacetylated derivatives. The three possessed a series of maltooligosaccharides and oligoglucosyl trehaloses with alpha1-->4 and alpha1-->1 linkages in common. However, B. germanica was conspicuous for its accumulation of a oligoglucosyl trehalose with a alpha1-->6 linkage [O-alpha-D-glucopyranosyl-(1-->6)-alpha-D-glucopyranosyl alpha-D-glucopyranoside] as the most prominent component of the secretion. Such diversification of the male pheromonal compositions within the genus might be reflected by biosynthetic systems in males and chemosensory properties in females suited for their natural habitat.