The discovery of Iberobaeniidae (Coleoptera: Elateroidea): a new family of beetles from Spain, with immatures detected by environmental DNA sequencing

The ongoing exploration of biodiversity and the implementation of new molecular tools continue to unveil hitherto unknown lineages. Here, we report the discovery of three species of neotenic beetles for which we propose the new family Iberobaeniidae. Complete mitochondrial genomes and rRNA genes recovered Iberobaeniidae as a deep branch in Elateroidea, as sister to Lycidae (net-winged beetles). Two species of the new genus Iberobaenia, Iberobaenia minuta sp. nov. and Iberobaenia lencinai sp. nov. were found in the adult stage. In a separate incidence, a related sequence was identified in bulk samples of soil invertebrates subjected to shotgun sequencing and mitogenome assembly, which was traced to a larval voucher specimen of a third species of Iberobaenia Iberobaenia shows characters shared with other elateroid neotenic lineages, including soft-bodiedness, the hypognathous head, reduced mouthparts with reduced labial palpomeres, and extremely small-bodied males without strengthening structures due to miniaturization. Molecular dating shows that Iberobaeniidae represents an ancient relict lineage originating in the Lower Jurassic, which possibly indicates a long history of neoteny, usually considered to be evolutionarily short-lived. The apparent endemism of Iberobaeniidae in the Mediterranean region highlights the importance of this biodiversity hotspot and the need for further species exploration even in the well-studied European continent.

Independent evolution of sexual dimorphism and female-limited mimicry in swallowtail butterflies (Papilio dardanus and Papilio phorcas)

Several species of swallowtail butterflies (genus Papilio) are Batesian mimics that express multiple mimetic female forms, while the males are monomorphic and nonmimetic. The evolution of such sex-limited mimicry may involve sexual dimorphism arising first and mimicry subsequently. Such a stepwise scenario through a nonmimetic, sexually dimorphic stage has been proposed for two closely related sexually dimorphic species: Papilio phorcas, a nonmimetic species with two female forms, and Papilio dardanus, a female-limited polymorphic mimetic species. Their close relationship indicates that female-limited polymorphism could be a shared derived character of the two species. Here, we present a phylogenomic analysis of the dardanus group using 3964 nuclear loci and whole mitochondrial genomes, showing that they are not sister species and thus that the sexually dimorphic state has arisen independently in the two species. Nonhomology of the female polymorphism in both species is supported by population genetic analysis of engrailed, the presumed mimicry switch locus in P. dardanus. McDonald-Kreitman tests performed on SNPs in engrailed showed the signature of balancing selection in a polymorphic population of P. dardanus, but not in monomorphic populations, nor in the nonmimetic P. phorcas. Hence, the wing polymorphism does not balance polymorphisms in engrailed in P. phorcas. Equally, unlike in P. dardanus, none of the SNPs in P. phorcas engrailed were associated with either female morph. We conclude that sexual dimorphism due to female polymorphism evolved independently in both species from monomorphic, nonmimetic states. While sexual selection may drive male-female dimorphism in nonmimetic species, in mimetic Papilios, natural selection for protection from predators in females is an alternative route to sexual dimorphism.

The mitogenome of Hydropsyche pellucidula (Hydropsychidae): first gene arrangement in the insect order Trichoptera

We describe the mitochondrial genome of Hydropsyche pellucidula Curtis 1834, which is first described for the suborder Annulipalpia and the first in the order Trichoptera to show a non-canonical gene order. The mitogenome was obtained by de novo assembly of shotgun sequenced total genomic DNA using Illumina Miseq technology, which produced an average coverage of 115× and a minimum coverage of 48×. The mitochondrial genome includes 13 protein-coding genes, 2 rRNAs and 22 tRNAs. The genome is characterized by a rearrangement in the relative position of protein-coding and ribosomal genes. This mitogenome sequence will be useful for studying the family Hydropsychidae, which is commonly used for freshwater pollution biomonitoring.

Why barcode? High-throughput multiplex sequencing of mitochondrial genomes for molecular systematics

Mitochondrial genome sequences are important markers for phylogenetics but taxon sampling remains sporadic because of the great effort and cost required to acquire full-length sequences. Here, we demonstrate a simple, cost-effective way to sequence the full complement of protein coding mitochondrial genes from pooled samples using the 454/Roche platform. Multiplexing was achieved without the need for expensive indexing tags ('barcodes'). The method was trialled with a set of long-range polymerase chain reaction (PCR) fragments from 30 species of Coleoptera (beetles) sequenced in a 1/16th sector of a sequencing plate. Long contigs were produced from the pooled sequences with sequencing depths ranging from ∼10 to 100× per contig. Species identity of individual contigs was established via three 'bait' sequences matching disparate parts of the mitochondrial genome obtained by conventional PCR and Sanger sequencing. This proved that assembly of contigs from the sequencing pool was correct. Our study produced sequences for 21 nearly complete and seven partial sets of protein coding mitochondrial genes. Combined with existing sequences for 25 taxa, an improved estimate of basal relationships in Coleoptera was obtained. The procedure could be employed routinely for mitochondrial genome sequencing at the species level, to provide improved species 'barcodes' that currently use the cox1 gene only.

Testing Introgressive Hybridization Hypotheses Using Statistical Network Analysis of Nuclear and Cytoplasmic Haplotypes in the Leaf Beetle Timarcha goettingensis Species Complex

Previous studies of leaf beetles (Chrysomelidae) in the Timarcha goettingensis species complex using mitochondrial (cox2) and nuclear (ITS-2 rRNA) markers revealed two main clades confined to the Iberian Peninsula and the rest of Europe but showing incongruent distributions indicative of gene exchange between both groups. Because of the anastomosing nature of hybridization, which disrupts the cladistic structure of character variation, phylogenetic trees might be inappropriate to represent and study this process. Here we test for evidence of hybridization in the T. goettingensis complex by analyzing the extra homoplasy arising in hybrid genomes from the simultaneous analysis of genetically independent markers. Haplotype networks obtained by Templeton's statistical parsimony analysis were generated for combined (concatenated) cox2 and ITS-2 sequences from 167 individuals of the T. goettingensis complex. Networks were used to detect runs of homoplasious characters physically clustered along a nucleotide sequence, as evidence for recombination between both gene partitions. A hypergeometric tail probability for the chance occurrence of physically clustered character changes on the connections linking networks of genotypes was applied. The test recognized two instances of statistically significant clustering, indicating the presence of cox2-ITS-2 mosaic genotypes and reticulation of both main T. goettingensis clades, supporting the reticulate origin of samples of T. maritima in southwestern France and T. sinuatocollis/T. monserratensis in the eastern Pyrenees. Although the assessment of reticulation in DNA sequences does not provide direct proof for hybridization, the geographical distribution of mosaic genotypes in the vicinity of "pure" genotypes supports the effect of gene flow between the two divergent lineages. The study demonstrates the utility of statistical parsimony networks for the detection of hybrids in the growing number of phylogeographic studies based on multiple gene markers.

DNA taxonomy, phylogeny and Pleistocene diversification of the Cicindela hybrida species group (Coleoptera: Cicindelidae)

Species delimitation is complicated where morphological variation is continuous or poorly subdivided, but for taxonomic convenience it is common practice to separate and name geographical groups to capture this variation. DNA-based approaches may be used to test if these groups in fact represent historically divided, discrete species entities. The Cicindela hybrida complex (Coleoptera: Cicindelidae) is an assemblage of up to seven morphologically recognized species and 15 subspecies with wide distribution in the Palaearctic region. We sequenced a discontinuous segment of 1899 bp of mtDNA including three regions (coxI, rrnL+trnL2+nad1, cob) for a total of 99 specimens from 36 sampling localities across Europe, revealing 48 haplotypes. Four major clades could be identified corresponding to geographical groups from central Iberia, Ukraine, central Europe, and a band from the Atlantic Iberian coast to northern Europe. Taking into account further subdivisions within these clades, four of the six named species included in the analysis were recognizable by applying various procedures for species delimitation. Age estimates from calibrated molecular clocks date the diversification of the hybrida group within the past 2 million years (Myr), and the separation of the northern clade within 0.4 Myr. Nested clade analysis revealed the rapid range expansion of the northern group consistent with postglacial dispersal, but we did not find support for specific source population(s) in the postulated southern refugia. The evolutionary framework based on mtDNA sequences is shown to identify species entities as discrete clusters of closely related sequences and provides an objective system for delineating and recognizing hierarchically structured groups. In the case of the C. hybrida complex, these groups largely coincided with those established from morphology. The study adds further support to the utility of mtDNA-based sequence profiles (the 'DNA taxonomy') as a rapid and objective synthesis of evolutionary diversity and as reference system for communication.

Incongruent nuclear and mitochondrial phylogeographic patterns in the Timarcha goettingensis species complex (Coleoptera, Chrysomelidae)

Phylogeographic analyses have mostly been based on single-gene genealogies but it is unclear how conclusions from such studies depend on the choice of gene markers. We conducted a nested geographical clade analysis [A.R. Templeton, E. Routman and C.A. Phillips (1995) Genetics 140: 767-782] based on nuclear rDNA internal transcribed spacer region 2 (ITS2) sequences in the Timarcha goettingensis species complex (Coleoptera, Chrysomelidae), and compared the inferences with an updated version of previously published results using mitochondrial cytochrome oxidase II (COIl) sequences. Inferences from ITS2 suggest that patterns of marker distribution are mostly explained by restricted gene flow with isolation by distance. In contrast, COII revealed a history of geographical structure resulting from episodic population contiguous-range expansions. Both markers also show different genealogical patterns, which are associated to the effects of genetic introgression in a putative hybrid zone between two major lineages in the complex. Altogether, these differences are attributed to distinct population and/or evolutionary dynamics of the markers, and offer a more accurate phylogeographic description for the T. goettingensis complex.

Mitochondrial DNA phylogeography and population history of Meladema diving beetles on the Atlantic Islands and in the Mediterranean basin (Coleoptera, Dytiscidae)

The phylogeny and population history of Meladema diving beetles (Coleoptera, Dytiscidae) were examined using mitochondrial DNA sequence from 16S ribosomal RNA and cytochrome oxidase I genes in 51 individuals from 22 populations of the three extant species (M. imbricata endemic to the western Canary Islands, M. lanio endemic to Madeira and M. coriacea widespread in the Western Mediterranean and on the western Canaries), using a combination of phylogenetic and nested clade analyses. Four main lineages are observed within Meladema, representing the three recognized species plus Corsican populations of M. coriacea. Phylogenetic analyses demonstrate the sister relationship of the two Atlantic Island taxa, and suggest the possible paraphyly of M. coriacea. A molecular clock approach reveals that speciation within the genus occurred in the Early Pleistocene, indicating that the Atlantic Island endemics are not Tertiary relict taxa as had been proposed previously. Our results point to past population bottlenecks in all four lineages, with recent (Late-Middle Pleistocene) range expansion in non-Corsican M. coriacea and M. imbricata. Within the Canary Islands, M. imbricata seems to have independently colonized La Gomera and La Palma from Tenerife (although a colonization of La Palma from La Gomera cannot be discarded), and M. coriacea has independently colonized Tenerife and Gran Canaria from separate mainland lineages. In the Mediterranean basin this species apparently colonized Corsica on a single occasion, relatively early in its evolutionary history (Early Pleistocene), and has colonized Mallorca recently on multiple occasions. On the only island where M. coriacea and M. imbricata are broadly sympatric (Tenerife), we report evidence of bidirectional hybridization between the two species.

Discovery of Aspidytidae, a new family of aquatic Coleoptera

The six extant aquatic families of Hydradephaga (Coleoptera) known so far represent a diverse group of beetles morphologically highly modified for life in the water. We report the discovery of a new genus with two species from South Africa and China, which differ greatly from all extant families, but resemble the Jurassic-Cretaceous dagger Liadytidae (the dagger symbol indicates that the taxa are known only as fossils). Based on a combined phylogenetic analysis of molecular and morphological data we erect a new family, Aspidytidae, which is the sister group of Dytiscidae plus Hygrobiidae. We propose a new scenario for the evolution of swimming behaviour in adephagan beetles, in which the transition into the aquatic environment is followed by complex and repeated changes in lifestyles, including the secondary complete loss of swimming ability in Aspidytidae.

Comparison of EST libraries from seven beetle species: towards a framework for phylogenomics of the Coleoptera

Relatively little is known about Coleoptera genes and genomes and how these compare in different taxa. We describe here the construction, DNA sequencing and sequence comparisons of cDNA libraries from seven beetle species. A total of 6717 bacterial colonies were screened for cDNA insert containing plasmids and 2784 size selected clones were 5'- and 3'-end sequenced to produce 1620 assembled sequences. Similarity comparisons with existing protein sequence databases revealed that 65.1% had matches (E < 10(-4)) in other organisms, with greater numbers of matches in Drosophila melanogaster than Caenorhabditis elegans and Saccharomyces cerevisiae databases. tBlastX comparisons also revealed numerous similarity hits (E < 10(-20)) in intra- and interlibrary comparisons. These results show the potential of small cDNA libraries for discovery and comparative analysis of genes useful for phylogenomic and functional studies.

Extensive trans-species mitochondrial polymorphisms in the carabid beetles Carabus subgenus Ohomopterus caused by repeated introgressive hybridization

To study the potential importance of introgressive hybridization to the evolutionary diversification of a carabid beetle lineage, we studied intraspecific and trans-species polymorphisms in the mitochondrial NADH dehydrogenase subunit 5 (ND5) gene sequence (1083 bp) in four species of the subgenus Ohomopterus (genus Carabus) in central and eastern Honshu, Japan. Of the four species, C. insulicola is parapatric with the other three, and can hybridize naturally with at least two. This species possesses two haplotypes of remote lineages. We classified ND5 haplotypes using polymerase chain reaction-restriction fragment length polymorphism with TaqI endonuclease for 524 specimens, and sequenced 143 samples. Analysis revealed that each species was polyphyletic in its mitochondrial DNA phylogeny, representing a marked case of trans-species polymorphism. Recent one-way introgression of mitochondria from C. arrowianus nakamurai to C. insulicola, and from C. insulicola to C. esakii, was inferred from the frequency of identical sequences between these species and from direct evidence of hybridization in their contact zones. Other intraspecific polymorphisms in the four species may be due to undetected introgressive hybridization (e.g. C. insulicola to C. maiyasanus) or from stochastic lineage sorting of ancestral polymorphisms. This beetle group has a genital lock-and-key system, with species-specific or subspecies-specific genital morphology that may act as a barrier to hybridization. However, our results demonstrate that introgressive hybridization has occurred multiple times, at least for mitochondria, despite differences among, and stability within, morphological characters that distinguish local populations. Thus, hybridization and introgression could have been key processes in the evolutionary diversification of Ohomopterus.

Exploring data interaction and nucleotide alignment in a multiple gene analysis of Ips (Coleoptera: Scolytinae)

The possibility of gene tree incongruence in a species-level phylogenetic analysis of the genus Ips (Coleoptera: Scolytidae) was investigated based on mitochondrial 16S rRNA (16S) and nuclear elongation factor-1 alpha (EF-1 alpha) sequences, and existing cytochrome oxidase I (COI) and nonmolecular data sets. Separate cladistic analyses of the data partitions resulted in partially discordant most-parsimonious trees but revealed only low conflict of the phylogenetic signal. Interactions among data partitions, which differed in the extent of sequence divergence (COI > 16S > EF-1 alpha), base composition, and homoplasy, revealed that much of the branch support emerges only in the simultaneous analysis, particularly for deeper nodes in the tree, which are almost entirely supported through "hidden support" (sensu Gatesy et al., Cladistics 15:271-313, 1999). Apparent incongruence between data partitions is in part due to suboptimal alignments and bias of character transformations, but little evidence supports invoking incongruent phylogenetic histories of genetic loci. There is also no justification for eliminating or downweighting gene partitions on the basis of their apparent homoplasy or incongruence with other partitions, because the signal emerges only in the interaction of all data. In comparison with traditional taxonomy, the pini, plastographus, and perturbatus groups are polyphyletic, whereas the grandicollis group is monophyletic except for inclusion of the (monophyletic) calligraphus group. The latidens group and some European species are distantly related and closer to other genera within Ipini. Our robust cladogram was used to revise the classification of Ips. We provide new diagnoses for Ips and four subgeneric taxa.

Sequence alignment of 18S ribosomal RNA and the basal relationships of Adephagan beetles: evidence for monophyly of aquatic families and the placement of Trachypachidae

Current hypotheses regarding family relationships in the suborder Adephaga (Coleoptera) are conflicting. Here we report full-length 18S ribosomal RNA sequences of 39 adephagans and 13 outgroup taxa. Data analysis focused on the impact of sequence alignment on tree topology, using two principally different approaches. Tree alignments, which seek to minimize indels and substitutions on the tree in a single step, as implemented in an approximate procedure by the computer program POY, were contrasted with a more traditional procedure based on alignments followed by phylogenetic inference based on parsimony, likelihood, and distance analyses. Despite substantial differences between the procedures, phylogenetic conclusions regarding basal relationships within Adephaga and relationships between the four suborders of Coleoptera were broadly similar. The analysis weakly supports monophyly of Adephaga, with Polyphaga usually as its sister, and the two small suborders Myxophaga and Archostemata basal to them. In some analyses, however, Polyphaga was reconstructed as having arisen from within Hydradephaga. Adephaga generally split into two monophyletic groups, corresponding to the terrestrial Geadephaga and the aquatic Hydradephaga, as initially proposed by Crowson in 1955, consistent with a single colonization of the aquatic environment by adephagan ancestors and contradicting the recent proposition of three independent invasions. A monophyletic Hydradephaga is consistently, though not strongly, supported under most analyses, and a parametric bootstrapping test significantly rejects an hypothesis of nonmonophyly. The enigmatic Trachypachidae, which exhibit many similarities to aquatic forms but whose species are entirely terrestrial, were usually recovered as a basal lineage within Geadephaga. Strong evidence opposes the view that terrestrial trachypachids are related to the dytiscoid water beetles.

The effect of habitat type on speciation rates and range movements in aquatic beetles: inferences from species-level phylogenies

Most aquatic beetles in the family Dytiscidae are tightly associated either with running (lotic) or stagnant (lentic) water bodies. The range size of lotic species is known to be, on average, much smaller than that of lentic species, presumably as a result of differences in dispersal strategies in each habitat type. We explored possible effects of these differences on clade evolution and speciation rates by comparing species-level phylogenies based on cytochrome oxidase I (COI) and 16S rRNA mitochondrial genes for two genera, the lentic Ilybius and the lotic Deronectes. The expectation that species turnover is higher in lotic lineages due to their lower dispersal propensity compared to lentic species was not strongly supported. Deronectes displays a higher frequency of recent splits than Ilybius, consistent with the hypothesis, but the difference was not significant compared to expected patterns under a constant speciation rate null model. Similarly, when the degree of sympatry was plotted against relative node age, more allopatric splits were evident in the lentic Deronectes, suggesting a slower rate of range movement since speciation, but the differences were not significant. We discuss two explanations for our failure to detect differences between the two clades. First, current methods for analysing species-level phylogenies may be sensitive to taxonomic and sampling artefacts. Second, lentic and lotic clades may indeed display similar levels of species turnover despite occupying very different habitats at different spatial scales. More work is needed to investigate the effects of population level processes and spatial scale on macroevolutionary dynamics.

Incongruence of mitochondrial and nuclear gene trees in the Carabid beetles Ohomopterus

We studied the molecular phylogeny of the carabid subgenus Ohomopterus (genus Carabus), using two mitochondrial (mt) DNA regions (16SrRNA and NADH dehydrogenase subunit 5) and three nuclear DNA regions (wingless, phosphoenolpyruvate carboxykinase, and an anonymous locus). We revisited the previously reported incongruence between the distribution of mtDNA markers and morphologically defined species (Su et al., 1996; J. Mol. Evol. 43:662-671), which those authors attributed to "type switching", a concerted change in many morphological characters that results in the repeated evolution of a particular morphological type. Our mtDNA gene tree obtained from 44 individuals representing all 15 currently recognized species of Ohomopterus revealed that haplotypes isolated from individuals of a single "species" were frequently separated into distant clades, confirming the previous report. The three nuclear markers generally conformed better-with the morphologically defined species than did the mitochondrial markers. The phylogenetic signal in mtDNA and nuclear DNA data differed strongly, and these two partitions were significantly incongruent with each other according to the incongruence length difference test of Farris et al. (1994; Cladistics 10:315-320), although the three nuclear partitions were not homogeneous either. Our results did not support the type-switching hypothesis that had been proposed to fit the morphological data to the mitochondrial gene tree: The incongruence of the mtDNA tree with other nuclear markers indicates that the mtDNA-based tree does not reflect species history any better than the morphological data do. Incongruence of gene trees in Ohomopterus may have been promoted by the complex processes of geographic isolation and hybridization in the Japanese Archipelago that have led to occasional gene flow and recombination between separated entities. The occurrence of reticulate patterns in this group is intriguing, because species of Ohomopterus exhibit extremely divergent genitalic structures that represent a highly efficient reproductive isolation mechanism.

Polymorphic microsatellite markers identified in individual Plasmodium falciparum oocysts from wild-caught Anopheles mosquitoes

The oocysts of Plasmodiom falciparam that form on the midgut wall of anopheline vectors represent the meiotic products of the sexual cycle of the parasite and are useful for genetic analyses. Recognizing the limitations of current markers for the analysis of genetic variation, we describe a protocol for the DNA isolation and PCR amplification of microsatellite loci from individual oocysts. We present a simple method for the storage and isolation of individual Plasmodium oocysts from wild-caught mosquitoes and demonstrate that it is possible to extract sufficient oocyst DNA for approximately 50 amplifications. Genotyping was successful for 10 microsatellite loci and, using Genescan technology, we found extensive microsatellite polymorphism in an analysis of a dozen oocysts from each of 2 Anopheles gambiae female specimens collected within 6 months at the same location in northern Tanzania. All of the loci exhibited at least 2 alleles, with a maximum of 5. Homozygous and heterozygous oocysts could easily be discerned, but most of the oocysts were homozygous. The procedure has great potential for shedding light on genetic mechanisms operating during meiosis in mosquito vectors and the influence this may have on the genetic structure of natural populations of P. falciparum.

How slippage-derived sequences are incorporated into rRNA variable-region secondary structure: implications for phylogeny reconstruction

We analyzed the type and frequency of mutational changes in hypervariable rRNA regions, using the highly length-variable region V4 of the small subunit rRNA locus of tiger beetles (Cicindelidae) as an example. Phylogenetic analysis of indels in closely related species showed that (1) most indels are single nucleotides (usually A or T and sometimes G) or di-nucleotides of A and T. These occur at numerous foci, and they exhibit a strong bias for duplication of 5' single and di-nucleotide motifs but not 3' motifs. (2) Insertions/deletions in stem-forming regions affected paired and unpaired bases with about equal frequency but they did not disrupt the secondary structure. (3) Recurring mutations involving short repeats of the same bases caused parallel evolution of similar sequence motifs in the rRNA of different lineages. The observed types of change are consistent with the propostion that slippage is the main mutational mechanism. Slippage-derived sequences tend to be self-complementary, and therefore the stem-loop structure could be self-organizing as a consequence of the underlying mutational mechanism. Thus, the secondary structure in the cicindelid V4 region may be conserved due to the dynamics of the mutational mechanism rather than to functional constraints. These processes may also have a tendency to produce similar primary sequences irrespective of phylogenetic associations. The findings have implications for sequence alignment in phylogenetic analysis and should caution against the use of secondary structure to improve the determination of positional homology in hypervariable regions.

Species richness: does flower power explain beetle-mania?

The huge species richness of beetles has been attributed to their colonisation of flowering plants, but a vegetarian diet may not be the sole secret of the beetles' success.

Modelling the secondary structures of slippage-prone hypervariable RNA regions: the example of the tiger beetle 18S rRNA variable region V4

Variable regions within ribosomal RNAs frequently vary in length as a result of incorporating products of slippage. This makes constructing secondary structure models problematic because base homology is difficult or impossible to establish between species. Here, we model such a region by comparing the results of the MFOLD suboptimal folding algorithm for different species to identify conserved structures. Based on the reconstruction of base change on a phylogenetic tree of the species and comparison against null models of character change, we devise a statistical analysis to assess support of these structures from compensatory and semi-compensatory (i.e. G.C to G.U or A.U to G.U) mutations. As a model system we have used variable region V4 from cicindelid (tiger beetle) small subunit ribosomal RNAs (SSU rRNAs). This consists of a mixture of conserved and highly variable subregions and has been subject to extensive comparative analysis in the past. The model that results is similar to a previously described model of this variable region derived from a different set of species and contains a novel structure in the central, highly variable part. The method we describe may be useful in modelling other RNA regions that are subject to slippage.

Phylogeny of North American Cicindela tiger beetles inferred from multiple mitochondrial DNA sequences

Tiger beetles in the genus Cicindela (Coleoptera: Cicindelidae) have been used as a model system for studies in ecology and conservation biology. Work on this group will greatly benefit from the availability of a phylogenetic hypothesis. We selected a representative sample of 23 North American Cicindela and 6 outgroups to reconstruct a phylogeny based on 1896 nucleotide positions from three mitochondrial genes (Cytochrome b, Cytochrome oxidase III, and 16S rRNA). Cladistic analysis of these three data sets yielded widely different tree topologies, but character conflict between them appears to be relatively low. The combined analysis of all data resulted in three similar shortest trees of 3453 steps. One of these was also recovered after successive weighting and was considered the best estimate of relationships. The most basal taxa of North American Cicindela (s.l.) were in the cosmopolitan subgenus Cylindera. The derived taxa were in the subgenus Cicindela (s. str.), a group dominating at higher latitudes in the Nearctic and Palearctic Region. The molecular analysis was essentially in agreement with the traditional classification which has been worked out based on male genitalic structures by E. Rivalier (1954, Rev. Entomol. Française 21:249-268). In the molecular analysis, Rivalier's species groups and subgenera were mostly found to be composed of closely related taxa but several of them were not monophyletic. Implicit in the traditional classification is a sequence from basal to derived groups which we found to be essentially reversed in the molecular analysis. We also discuss the conceptual differences in the establishment of the traditional classification by Rivalier (1954) and the cladistic analysis presented in this study.

Phylogenetic analysis of slippage-like sequence variation in the V4 rRNA expansion segment in tiger beetles (Cicindelidae)

Sequence variation in the middle part of the small-subunit rRNA was studied for representatives of the major groups in the family Cicindelidae (Coleoptera). All taxa exhibited a much expanded segment in variable region V4 compared to D. melanogaster. This expanded segment was not found in other groups of beetles, including three taxa in the closely related Carabidae. Secondary structure predictions indicate that the expanded segment folds into a single stem-loop structure in all taxa. Despite its structural conservation, the fragment differs strongly in primary sequence, even between closely related sister taxa. Several features of these sequences are consistent with slippage replication as the mechanism that has generated this sequence variation: the level of internal sequence repetition as measured by the relative simplicity factor (RSF), its variation in length between close relatives, and the strong nucleotide bias compared to the remainder of the gene. With few exceptions, there was also a correlation between sequence length and the level of sequence repetition, frequently interpreted as the result of slippage. Phylogenies inferred from the expansion segment were not consistent with existing hypotheses from other molecular data for the group. This indicates that DNA sequences in this region are not homologous throughout the entire Cicindelidae, but it leaves open the possibility that this expansion segment can be used for phylogeny reconstruction within subgroups. The implications of a phylogenetic approach to the understanding of slippage-like evolution are discussed.

A molecular phylogeny of the tiger beetles (Cicindelidae): congruence of mitochondrial and nuclear rDNA data sets

As a basis for comparative studies in the family Cicindelidae (tiger beetles), a phylogenetic hypothesis for major groups was established based on molecular data, focusing on the taxa of the New World. Two markers, mitochondrial 16S rRNA and nuclear 18S rRNA, were sequenced for a total of 912 nucleotides and 30 taxa. All Cicindelidae were inferred to be a monophyletic group by the presence of a segment in the 5' region of the 18S rDNA which is predicted to fold into a single stem-loop structure not present in the outgroup. The total number of variable and potentially informative positions was smaller in the 18S data set, but their higher internal consistency resulted in a strong phylogenetic signal. Inferred relationships were stable over a wide range of alignment parameters, but the inclusion of alignment-sensitive positions and presumptive gap characters added substantially to the resolution and overall support of the tree. Both data sets were in general congruence. The combined evidence suggests that relationships in Cicindelidae as currently recognized have to be reconsidered. No convincing evidence could be found for the traditional subdivision into two subfamilies, Cicindelinae and Collyrinae. The latter is embedded within the Cicindelinae and clearly not the sister to all other cicindelids. Most of the current tribes and subtribes, particularly in the Megacephalini, are paraphyletic groups, whereas monophyly of the large tribe Cicindelini is well supported by the molecular data.

Homologues of the engrailed gene from five molluscan classes

We used the polymerase chain reaction (PCR) to amplify, clone, and sequence 10 engrailed homeodomains from 8 species in the five major molluscan classes, including the serially organized chiton (Polyplacophora) lineage. The Drosophila melanogaster gene engrailed (en) is one of several genes involved in embryonic segment polarity determination. Studies of engrailed sequence and expression in molluscs are of interest due to questions regarding the evolution and homology of segmentation in these taxa. Nucleotide and deduced amino acid sequence comparisons reflect evolutionary conservation within helices of the en homeodomain and ancient divergences in the region 3' to the homeodomain.

Evolution and phylogenetic information content of the ITS-1 region in the tiger beetle Cicindela dorsalis

Sequence divergence in the internal transcribed spacer region 1 (ITS-1) of the ribosomal DNA locus was assessed in subspecies of the coastal North American tiger beetle, Cicindela dorsalis. The spacer region was amplified using the polymerase chain reaction and cloned for sequencing. Of a total of 50 clones obtained from 12 specimens, 42 clones were different in at least one nucleotide position. In a parsimony analysis of these sequences, the main phylogenetic distinction was found to separate sequences from the Gulf of Mexico and the Atlantic Ocean. Within these two assemblages phylogenetic resolution was low, and the variation within individuals was almost as high as the variation within the entire lineage. The pattern of sequence variation suggests the existence of two forms of the ITS-1 that are maintained on different chromosomes. Polymorphisms of limited geographical distribution could be detected, and 41 additional clones were partly sequenced, to assess the geographic distribution of these polymorphisms in more detail. In a population aggregation analysis, the geographic pattern of ITS-1 distribution was basically congruent with that obtained in earlier studies from mitochondrial DNA in the same C. dorsalis populations.

Using molecular and ecological data to diagnose endangered populations of the puritan tiger beetle Cicindela puritana

Populations of the puritan tiger beetle Cicindela puritana in the eastern United States were found to be highly threatened at the Connecticut River, whereas several large populations on the western shore and newly discovered populations on the eastern shore of the Chesapeake Bay appeared to be less endangered. We assessed if the disjunct C. puritana subgroups are genetically distinct and therefore should be treated as separate units for conservation purposes. A total of 13 individuals from the Connecticut River and 27 individuals from the Chesapeake Bay were each analysed by sequencing of up to 837 base pairs of mitochondrial DNA per individual. Five different haplotypes could be distinguished. In a phylogenetic analysis of these DNA sequences that included four related Cicindela species as out-groups, haplotypes from the Chesapeake Bay represent a distinct clade. The conservation status of these populations was evaluated using a phylogenetic approach based on cladistic analysis and the framework of the phylogenetic species concept. According to this analysis, beetles from the Connecticut River and the Chesapeake Bay have to be considered as independent units. Populations from the eastern and western shore of Chesapeake Bay are not split in more than one unit using the same criteria, although they exhibited some degree of genetic subdivision. The results from the mtDNA analysis were corroborated by ecological parameters in that the Chesapeake Bay populations can be distinguished from all congeners by their different habitat association.

Comparison of the sequences of the nagE operons from Klebsiella pneumoniae and Escherichia coli K12: enhanced variability of the enzyme IIN-acetylglucosamine in regions connecting functional domains

The nagE operon, encoding the enzyme II specific for N-acetylglucosamine (EIINag), and adjacent DNA from the chromosome of Klebsiella pneumoniae were sequenced and compared with the corresponding sequence from Escherichia coli K12. The deduced EIINag sequences differ in 72 out of 651 amino acids, the K. pneumoniae sequence being three residues longer. The amino acid differences were distributed unevenly, and were most frequent in regions connecting the three functional domains of the protein. In the nagE-nagB intergenic region, two promoter, two operator, and one CAP consensus sequence with regulatory functions were highly conserved. The nag structural genes from both species were very similar (83% DNA similarity; 89% amino acid similarity) except for frequent AT to GC exchanges in the wobble base of codons in K. pneumoniae DNA relative to the E. coli DNA.

Salmonella typhimurium mutants defective in flagellar filament regrowth and sequence similarity of FliI to F0F1, vacuolar, and archaebacterial ATPase subunits

Many flagellar proteins are exported by a flagellum-specific export pathway. In an initial attempt to characterize the apparatus responsible for the process, we designed a simple assay to screen for mutants with export defects. Temperature-sensitive flagellar mutants of Salmonella typhimurium were grown at the permissive temperature (30 degrees C), shifted to the restrictive temperature (42 degrees C), and inspected in a light microscope. With the exception of switch mutants, they were fully motile. Next, cells grown at the permissive temperature had their flagellar filaments removed by shearing before the cells were shifted to the restrictive temperature. Most mutants were able to regrow filaments. However, flhA, fliH, fliI, and fliN mutants showed no or greatly reduced regrowth, suggesting that the corresponding gene products are involved in the process of flagellum-specific export. We describe here the sequences of fliH, fliI, and the adjacent gene, fliJ; they encode proteins with deduced molecular masses of 25,782, 49,208, and 17,302 Da, respectively. The deduced sequence of FliI shows significant similarity to the catalytic beta subunit of the bacterial F0F1 ATPase and to the catalytic subunits of vacuolar and archaebacterial ATPases; except for limited similarity in the motifs that constitute the nucleotide-binding or catalytic site, it appears unrelated to the E1E2 class of ATPases, to other proteins that mediate protein export, or to a variety of other ATP-utilizing enzymes. We hypothesize that FliI is either the catalytic subunit of a protein translocase for flagellum-specific export or a proton translocase involved in local circuits at the flagellum.

Involvement of the histidine protein (HPr) of the phosphotransferase system in chemotactic signaling of Escherichia coli K-12

It is known that in mutants of Escherichia coli lacking the histidine protein (HPr) of the carbohydrate: phosphotransferase system, all substrates of the system can be taken up in the presence of the fructose-regulated HPr-like protein FPr (gene fruF). Although this protein fully substituted for HPr in transport and phosphorylation, we found that it was not able to complement efficiently for HPr in mediating chemotaxis toward phosphotransferase system substrates. Furthermore, transport activity and chemotaxis could be genetically dissected by the exchange of single amino acids in HPr. The results suggest a specific role of HPr in chemotactic signaling. We propose a possible link of signal transduction pathways for phosphotransferase system- and methyl chemotaxis protein-dependent substrates via HPr.

Structures and homologies of carbohydrate: phosphotransferase system (PTS) proteins

The bacterial phosphotransferase system (PTS) is the major transport system for many carbohydrates that are phosphorylated concomitantly with the translocation step through the membrane (group translocation). It consists of two general proteins, enzyme I and histidine protein (HPr), and a series of more than 15 substrate-specific enzymes II (EII). The sequences of several of these derived from Gram-positive and Gram-negative bacteria were compared, which allowed the possible identification of the following functional domains: membrane-bound pore, substrate-binding site, linker domains, transphosphorylation domain and primary phosphorylation site. Several EIIs have been analysed in the meantime, also by topological tests, by sequential deletion of the corresponding structural genes, and by construction of intergenic hybrids between different domains of several EIIs. These data suggest evolutionary relationships between different EIIs; they also enable a general model to be constructed of EIIs as carbohydrate transport systems, phosphotransferases, chemoreceptors in chemotaxis and as part of a global regulatory network.

Alternative route for biosynthesis of amino sugars in Escherichia coli K-12 mutants by means of a catabolic isomerase

By inserting a lambda placMu bacteriophage into gene glmS encoding glucosamine 6-phosphate synthetase (GlmS), the key enzyme of amino sugar biosynthesis, a nonreverting mutant of Escherichia coli K-12 that was strictly dependent on exogenous N-acetyl-D-glucosamine or D-glucosamine was generated. Analysis of suppressor mutations rendering the mutant independent of amino sugar supply revealed that the catabolic enzyme D-glucosamine-6-phosphate isomerase (deaminase), encoded by gene nagB of the nag operon, was able to fulfill anabolic functions in amino sugar biosynthesis. The suppressor mutants invariably expressed the isomerase constitutively as a result of mutations in nagR, the locus for the repressor of the nag regulon. Suppression was also possible by transformation of glmS mutants with high-copy-number plasmids expressing the gene nagB. Efficient suppression of the glmS lesion, however, required mutations in a second locus, termed glmX, which has been localized to 26.8 min on the standard E. coli K-12 map. Its possible function in nitrogen or cell wall metabolism is discussed.

Analysis of the nag regulon from Escherichia coli K12 and Klebsiella pneumoniae and of its regulation

Four genes, nagR, A, B and E, clustered in the nag locus of Escherichia coli K12 and Klebsiella pneumoniae, were cloned and physically mapped, and the corresponding gene products involved in amino sugar metabolism identified. Expression of the nag genes was also analysed using a series of lacZ fusions. In both bacteria, the genes are arranged in two divergent operons and controlled by a common NagR repressor. The corresponding gene nagR was found to map in the first operon together with the promoter proximal gene nagB, encoding the enzyme D-glucosamine isomerase (deaminase) (NagB) and the middle gene nagA, coding for N-acetyl-glucosamine deacetylase (NagA). Polar mutations in nagB and nagA prevent the efficient expression of nagR and cause constitutive expression of all nag genes. This includes the gene nagE encoding Enzyme IINag of the phosphoenolpyruvate-dependent carbohydrate phosphotransferase system (PTS), encoded in the second divergently transcribed operon. No further gene is found in this operon which in both organisms is directly adjacent to the gene glnS. It is interesting that the NagR repressor also affects the mannose PTS (genes manX, Y, Z), the second transport system involved in amino sugar uptake and phosphorylation.

Suppression of IIIGlc-defects by enzymes IINag and IIBgl of the PEP:carbohydrate phosphotransferase system

The Enzymes II of the PEP:carbohydrate phosphotransferase system (PTS) specific for N-acetylglucosamine (IINag) and beta-glucosides (IIBgl) contain C-terminal domains that show homology with Enzyme IIIGlc of the PTS. We investigated whether one or both of the Enzymes II could substitute functionally for IIIGlc. The following results were obtained: (i) Enzyme IINag, synthesized from either a chromosomal or a plasmid-encoded nagE+ gene could replace IIIGlc in glucose, methyl alpha-glucoside and sucrose transport via the corresponding Enzymes II. An Enzyme IINag with a large deletion in the N-terminal domain but with an intact C-terminal domain could also replace IIIGlc in IIGlc-dependent glucose transport. (ii) After decryptification of the Escherichia coli bgl operon, Enzyme IIBgl could substitute for IIIGlc. (iii) Phospho-HPr-dependent phosphorylation of methyl alpha-glucoside via IINag/IIGlc is inhibited by antiserum against IIIGlc as is N-acetylglucosamine phosphorylation via IINag. (iv) In strains that contained the plasmid which coded for IINag, a protein band with a molecular weight of 62,000 D could be detected with antiserum against IIIGlc. We conclude from these results that the IIIGlc-like domain of Enzyme IINag and IIBgl can replace IIIGlc in IIIGlc-dependent carbohydrate transport and phosphorylation.

Complementation of a truncated membrane-bound Enzyme IINag from Klebsiella pneumoniae with a soluble Enzyme III in Escherichia coli K12

Cloning and analysis of the gene nagE encoding Enzyme IINag (EIINag) from Klebsiella pneumoniae revealed strong similarities with the corresponding gene from Escherichia coli K12. Truncated EIINag proteins were generated by inserting a series of Tn1725 transposons into the structural gene; the positions of the insertions were mapped by restriction enzyme analysis, and the activity of the polypeptides determined by in vitro and in vivo tests. Insertions in the region encoding the amino-terminal half of the protein invariably abolished transport and phosphorylation activity, while truncated proteins lacking a C-terminal domain homologous to the soluble Enzyme III (crr gene) could be complemented by this molecule to nearly wild-type activity.

Indirect role of adenylate cyclase and cyclic AMP in chemotaxis to phosphotransferase system carbohydrates in Escherichia coli K-12

Most strains of Escherichia coli K-12 lacking the enzyme adenylate cyclase showed normal chemotaxis toward carbohydrates taken up and phosphorylated by the phosphoenolpyruvate-dependent carbohydrate: phosphotransferase system. The normal reaction was observed even in the absence of externally added cyclic adenosine 3',5'-phosphate, provided that the enzyme II chemoreceptors and the flagella were synthesized. In the CA8306 series of strains, however, the cya-854 deletion abolished chemotaxis toward phosphotransferase system carbohydrates even though growth on and transport of these carbohydrates were not affected. This abnormal phenotype was due to the presence of a specific mutation in strain CA8306 which mapped in or close to the crp locus and apparently prevented expression of a hitherto unidentified molecule involved in enzyme II-mediated signal transduction. This molecule is neither a pts protein nor a cyclic adenosine 3',5'-phosphate-binding protein.