Research Article
Print
Research Article
Species limits in polymorphic mimetic Eniclases net-winged beetles from New Guinean mountains (Coleoptera, Lycidae)
expand article infoMatej Bocek, Ladislav Bocak
‡ Palacky University in Olomouc, Olomouc, Czech Republic
Open Access

Abstract

Species delimitation was compared in a group of closely related lineages of aposematically colored Eniclases (Coleoptera, Lycidae) using morphology, genetic distances, and Bayesian implementation of the Poisson Tree Processes model. A high diversity of net-winged beetles was found in previously unsampled regions of New Guinea and ten new species are described: Eniclases bicolor sp. n., E. bokondinensis sp. n., E. brancuccii sp. n., E. elelimensis sp. n., E. infuscatus sp. n., E. niger sp. n., E. pseudoapertus sp. n., E. pseudoluteolus sp. n., E. tikapurensis sp. n., and E. variabilis sp. n. Different levels of genetic and morphological diversification were identified in various sister-species pairs. As a result, both morphological and molecular analyses are used to delimit species. Sister-species with uncorrected pairwise genetic divergence as low as 0.45% were morphologically distinct not only in color pattern, but also in the relative size of eyes. Conversely, differences in color pattern regardless of their magnitude did not necessarily indicate genetic distance and intraspecific mimicry polymorphism was common. Additionally, genetic divergence without morphological differentiation was detected in one sister-species pair. Low dispersal propensity, diverse mimicry patterns, and mimetic polymorphism resulted in complex diversification of Eniclases and uncertain species delimitation in recently diversified lineages.

Keywords

Aposematism, bPTP model, Coleoptera , cox1 mtDNA, genetic distance, morphology, new species, species delimitation

Introduction

The Papuan beetle fauna remains one of the most under-investigated despite high species richness in the Oceanian region and a long history of taxonomic research in Papua (e.g. Kleine 1926). Recent studies show exceptional diversity (Riedel et al. 2013, Toussaint et al. 2014). Previous reports on Papuan net-winged beetles included chaotic classification and poor species delimitation resulting from formal morphological descriptions using a semaphoront without diagnostic characters. DNA data provide a tool to accelerate biodiversity research, but these data must be viewed as part of a process of cross-validation of hypotheses on species limits based on both morphology and genetic information (e.g. Meyer and Paulay 2005, Meier et al. 2006, Baselga et al. 2013, Jorger and Schrodl 2013). Here, we present a taxonomic study dealing with Eniclases Waterhouse, 1879 (Metriorrhynchini), a genus in which most species are aposematically colored. They belong to numerous mimetic complexes, mostly formed by net-winged beetles, which are unpalatable (Bocak and Yagi 2010). Additionally, the net-winged beetles have weak dispersal ability and consequently prone to rapid differentiation of local populations (Li et al. 2015).

Most Papuan net-winged beetles belong to genera known from Australia (Calder 1998, Sklenarova et al. 2014), but some including Eniclases, are endemic. These lineages diversified probably in New Guinea and adjacent islands (Sklenarova et al. 2013) and are diverse despite the supposedly short existence of New Guinea in its present form. New Guinea is a young landmass resulting from uplift of the northern margin of the Australian tectonic plate and accretion of oceanic islands about 5 million years ago (Hill and Hall 2002). Eniclases was revised by Bocak and Bocakova (1991), including valid 27 species. The morphology-based revision of Bocak and Bocakova (1991) primarily investigated material from the eastern part of New Guinea. The fauna of Western New Guinea remained poorly studied with few species reported from the lowlands at the northern coast (4 spp.), the Panai Lake region (3 spp.), the Fak Fak Peninsula (2 spp.), and the Star Mts. (Oksibil area, 2 spp.). The faunas of the Central Mts. region and the Bird’s Head Peninsula were unknown.

A taxonomic study based on material representing Eniclases from the western part of the island, mainly from the Central Mts, is presented. The aim was to compare species delimitations inferred from genetic distances (barcoding approach; Hebert et al. 2003, Meier et al. 2006) and the tree shape (the Bayesian Poisson Tree Processes method, bPTP; Zhang et al. 2013). The status of these putative species inferred from mtDNA data was tested by the presence of morphological characters. Specifically, we investigate intraspecific color polymorphism and morphological diagnostic characters. Descriptions and diagnoses of new species are presented using available evidence and possible scenarios for diversification of Eniclases are discussed.

Methods

Material and laboratory procedures

In total, 81 specimens of Eniclases from the western part of New Guinea: the Central, Cyclops and Arfak Mts. were kept in -20 °C until DNA extraction. Each specimen was given a voucher number (Figs 1–2 and Taxonomy section, where exact locality data are given), voucher specimens are deposited in the collection of the Laboratory of Molecular Systematics, Olomouc and all sequences were submitted to GeneBank (Accession Numbers KT256092–172). The cox1 + tRNA-Leu + cox2 mtDNA fragment (hereafter referred as cox1 only) was amplified using primers JerryM (CAACAYYTATTTTGRTTYTTTGG) and Marcy (TARTTCRTATGWRCAATAYCAYTGRTG) or JerryN (CAACAYYTATTYTGATTYTTYGG) and MarcyN (TTCRTAWGTTCARTATCATTGRTG). DNA extraction, PCR settings and cycle sequencing conditions follow Bocak and Bocakova (2008). The PCR products were purified using PCRμ96™ Plates (Millipore Inc.) and sequenced by an ABI 3130 automated sequencer using the BigDye® Terminator Cycle Sequencing Kit 1.1.

Figures 1–2. 

1 Phylogenetic tree of Eniclases inferred from the maximum likelihood optimality criterion; a basal part of outgroups omitted. The numbers at branches show bootstrap support greater than 50%, genetic divergence within respective putative species and posterior probabilities inferred from the bPTP model 2 Dated tree produced using Bayesian inference.

Phylogenetic analyses and species delimitation

Sequences were edited using the Sequencher 4.8 software package (Gene Codes Corp.) and combined with homologous sequences of 25 outgroup taxa representing Cautirina and Metriorrhynchina (all sequences taken from GenBank). The length invariable cox1 mtDNA fragment was aligned using ClustalX 2.1 (Thompson et al. 1997) under default parameters and the phylogenetic analysis was carried out under the maximum likelihood criterion using RAxML 7.3.1 (Stamatakis 2006) and the GTR + I + G model for all partitions identified by jmodeltest 2.1.7 (Darriba et al. 2012). All genes and codon positions in the protein coding fragments were partitioned and parameters unlinked for each partition. Bootstrap values (BS) were assessed by analyzing 1000 pseudoreplicates using the rapid bootstrap algorithm under the GTRCAT model (Stamatakis et al. 2008). Trees were edited and visualized in Dendroscope 3.1 (Huson et al. 2007). In order to estimate the origin of closely related lineages with different mimetic patterns, the tree was dated using Beast 1.8.1 (Drummond et al. 2012) with the same model as in the maximum likelihood analysis. As there is no fossil record of metriorrhynchine Lycidae, we relied on the molecular evolution rates used earlier and tested two rates, 0.0115 substitutions per lineage per million years (s/l/my, Brower 1994), which gave results in agreement with tectonics in earlier study on Metriorrhynchus (Bocak and Yagi 2010) and the rate 0.0177 s/l/my calculated by Papadopoulou et al. (2010) for Tenebrionidae. We used the uncorrected lognormal clock model, Speciation: Birth Death, and sampled every 1,000 generations from a 10 million generation MCMC chain. The first 4 million generations were discarded as burn-in after evaluation of likelihood values and the effective sample size >1000 in Tracer 1.6 (Rambaut et al. 2013).

The genetic intra- and interspecific differentiation among whole-fragment sequences was estimated using Kimura 2-parameter genetic distances in MEGA6 (Tamura et al. 2013) and clusters of similar sequences were identified in Species Identifier 1.7.7 using a 3% barcoding threshold (Meier et al. 2006). The densities for intra- and interspecific differences were visualized in R (http://r-project.org). Further, species level entities were defined using the maximum likelihood and Bayesian implementation of the Poisson Tree Processes (bPTP) model for species delimitation (the bPTP server at species.h-its.org; Zhang et al. 2013).

Morphological characters were observed for all sequenced specimens: measurements of the body size, maximum diameter of eyes in the lateral view, the minimum interocular distance in the frontal part of cranium, color patterns of the pronotum and elytra, shape of pronotum and antennae, and structure of the elytral costae. Photographs were taken using a binocular microscope Olympus SZX-16 and were assembled in Helicon Focus 6 (www.heliconsoft.com). Due to previously reported uniformity of the genitalia of both sexes (Bocak and Bocakova 1991), those of only a few species were illustrated.

Abbreviations

DEIM Deutsches Entomologisches Institute, Müncheberg, Germany

MHNP Muséum national d’Histoire naturelle, Paris, France

LMBC Voucher collection, Department of Zoology, UP Olomouc, Czech Republic

Results

Molecular phylogeny

The 1101 bp fragment of mtDNA was sequenced for 81 individuals of Eniclases from western New Guinea. The DNA fragment consisted of 782 base pairs (bp)cox1 mtDNA, 59 bp tRNA-Leu, and 260 bpcox2 mtDNA. We identified 38 unique haplotypes and the Kimura-2-parameter genetic distances among Eniclases sequences spanned 0.09–14.31%. The maximum likelihood analysis produced the tree in Fig. 1. Eniclases formed a monophylum with 100% bootstrap support and the fauna of the Central Mts. was represented by four deeply rooted clades of closely related species (Fig. 1). The dated tree inferred with 0.0115 s/l/my rate is shown in Fig. 2 and suggested that closely related species pairs originated within in the last million years. The higher rate of 1.77% produced even shallower splits for the terminal lineages (results not shown).

The genetic divergence between all haplotypes was evaluated and the highest densities were between 0.0–2.0% and 9.0–14% (Fig. 4). Further, morphological divergence was considered within genetically close terminals and three morphologically distinguishable species pairs were identified: E. elelimensis and E. bokondinensis (mean interspecific K-2-P distance 0.79%, maximum intraspecific distance 0.17%); E. infuscatus and E. bicolor (1.07% and 0.12%); E. niger and E. similis (0.59% and 0.16%). Relatively low divergence was found between E. variabilis and E. elelimensis (<1.59%), E. variabilis and E. bokondinensis (<1.72%) and between E. brancuccii and three species in its sister clade (3.90–4.59%). The clade of E. tikapurensis and Eniclases sp. A consisted of two sister-subclades with inter-clade genetic difference 1.12–1.28% (Fig. 1) and these candidate species were morphologically indistinguishable.

Figures 3–5. 

3 The positions, coordinates and elevations of sampled localities in the Central Mountains of New Guinea 4 Density plots of genetic distances of all Eniclases samples 5 Density plots of intra- and interspecific genetic distances between pairs of closely related species of Eniclases (E. infuscatus and E. bicolor; E. tikapurensis and Eniclases sp. A; E. variabilis, E. elelimensis and E. bokondinensis; E. niger and E. similis).

Furthermore, putative species were identified using genetic distance and the phylogenetic tree. The pairwise differences among cox1 mtDNA sequences merged haplotypes in 9 clusters when the threshold 3% was applied. The inferred clusters merged E. infuscatus and E. bicolor, E. tikapurensis + Eniclases sp. A, the clade E. elelimensis + E. bokondinensis + E. variabilis and the pair E. niger + E. similis (Fig. 1, the clades labeled with triangles and intra-clade genetic distances). Most species differed in coloration, often additionally in eye size and geographic origin (see taxonomy section for detailed morphological descriptions). The maximum likelihood PTP method proposed identical species limits as the distance approach (results not shown). Unlike these, the Bayesian implementation of the PTP model suggested 9–19 putative species, mean 11.97 species. The most supported partitions, considered further as putative species, are displayed in Figure 1 (12 putative species designated by squares and labels designating posterior probabilities, filled green squares designate partitions with the highest PPs, empty squares suboptimal, non-preferred partitions). The posterior delimitation probabilities (PP) were high only for E. apertus and E. brancuccii partitions (≥90%); other clades obtained moderate to very low PPs. On the other hand, the alternative species delimitations corresponding to those from distance analyses and delimitations based on morphological traits had even lower PPs (Fig. 1).

Fourteen species were identified in the sequenced material, including ten new to science and described in the Taxonomy section. The delimitation of species was based on morphological characters (the size of eyes, shape of male antennae and coloration), genetic distance, and delimitations of putative species inferred from the bPTP model.

Taxonomy

Eniclases Waterhouse, 1879

Eniclases Waterhouse, 1879: 66.

Type species

Lycus luteolus Waterhouse, 1878 by monotypy.

Diagnosis

Eniclases is similar to Trichalus Waterhouse, 1877 and they share a small to medium sized, dorso-ventrally flattened body; a characteristic shape of a pronotum with acutely projecting posterior angles and a hump in a posterior third of a pronotal margin; nine elytral costae in the humeral part of elytra (four robust primary costae and five weak secondary costae) and a shortened primary costa 1 (Figs 650). Unlike Trichalus, the median pronotal areola is absent in Eniclases and only two divergent longitudinal keels are present in the pronotum (Figs 629). Additionally, all Eniclases have a characteristic cap-shaped apex of the phallus (Figs 43–50) similar to those of Schizotrichalus Kleine, 1926 (Bocak 2002). The detailed redescription of Eniclases was published by Bocak and Bocakova (1991).

Figures 6–17. 

Habitus of Eniclases: 6 E. divaricatus, male 7 E. pseudoapertus sp. n., male 8 E. apertus, male 9 E. tikapurensis sp. n., male 10 E. bicolor sp. n., female 11 E. infuscatus sp. n., female 12 E. brancuccii sp. n., female 13, 14 E. similis, male 15, 16 E. niger sp. n., male 17 E. bokondinensis sp. n., female. Scale bars: 2 mm.

Figures 18–29. 

Habitus of Eniclases: 18 E. bokondinensis sp. n., female 19 E. elelimensis, female 20–29 E. variabilis sp. n. Scale bars: 2 mm.

Phylogenetic relationships

The shape of the pronotum, shortened elytral costa one, and shape of male genitalia (Figs 650) support close relationships of Eniclases and other trichaline genera. At present, the trichaline lineages form a subordinate clade within Metriorrhynchina and consist of Eniclases, Flabellotrichalus Pic, 1921, Microtrichalus Pic, 1921, Schizotrichalus Kleine, 1926, and Trichalus Waterhouse, 1877 (Bocak 2002, Sklenarova et al. 2014).

Figures 30–50. 

Basal antennomeres of (Figs 30–42): 30 E. apertus, male 31 E. bicolor sp. n., female 32 E. bokondinensis sp. n., female 33 E. brancuccii sp. n., male 34 E. divaricatus, male 35 E. infuscatus sp. n., female 36 E. pseudoluteolus, male 37 E. similis, male 38 E. variabilis sp. n., male 39 E. elelimensis, male 40 E. niger sp. n., female 41 E. tikapurensis, male 42 E. pseudoapertus sp. n., male. Male genitalia, ventral view (Figs 43–50): 43 E. pseudoapertus sp. n. 44 E. divaricatus 45 E. tikapurensis sp. n. 46 E. brancuccii sp. n. 47 E. variabilis sp. n. 48 E. pseudoluteolus sp. n. 49 E. niger sp. n. 50 E. similis. Scale bars: 0.5 mm.

Eniclases pseudoapertus sp. n.

Figs 7, 42

Material examined

Holotype. Male (Voucher number BM0080), Indonesia, Irian Jaya, Arfak Mts., Maibri village, 1570 m, Nov.–Dec. 1991 (LMBC).

Diagnosis

Eniclases pseudoapertus resembles E. apertus Pic, 1923 in the small body and black coloration (Fig. 7); E. apertus has smaller eyes than E. pseudoapertus (Tab. 1).

Table 1.

Measurements of Eniclases spp. (all measurements in mm, n.a. – not available).

Body length Width humeri Pronotum Eye diameter/distance
length width male female
E. pseudoapertus 6.3 1.6 0.75 1.2 1.4 n.a.
E. divaricatus 8.8–9.7 2.1–2.3 1.2–1.3 1.7–1.7 0.92–0.96 0.80–0.87
E. pseudoluteolus 9.3 2.3 1.15 1.6 0.9 n.a.
E. apertus 5.8–6.9 1.4–1.7 0.9 1.25 1.15–1.17
E. tikapurensis 9.5–11.1 2.0–2.5 1.1–1.3 1.4–1.7 1.11–1.40 0.92–1.06
E. bicolor 10.3 2.4 1.4 1.7 n.a. 0.71–0.74
E. infuscatus 12.1 2.5 1.25 1.6 n.a. 0.79–0.84
E. brancuccii 7.6–8.0 1.8–1.9 1.0–1.1 7.5–8.0 1.00 0.84–0.91
E. bokondinensis 9.2 2.05 1.0 1.35 n.a. 0.72–0.82
E. elelimensis 6.9–8.1 1.5–1.9 0.9–1.1 1.3–1.4 n.a. 0.78–0.89
E. variabilis 6.6–8.2 1.6–2.0 0.1–1.1 1.1–1.35 0.83–0.95 0.70–0.85
E. niger 9.2–11.6 2.2–2.8 1.3–1.6 9.0–11.5 1.17–1.28 0.89
E. similis 7.5–9.7 1.9–2.3 1.1–1.4 1.8–1.8 1.02–1.15 0.89

Description

Male. Body length 6.3 mm, uniformly dark colored, only trochanters and bases of femora light brown (Fig. 7). Head with large hemispherically prominent eyes, their maximum diameter 1.40 times minimal interocular distance. Antennae slender, serrate, almost parallel-sided, with short apical process of antennomere 3 (Fig. 42). Pronotum transverse, with almost straight sides, elytra with weak secondary costae and irregular cells.

Distribution

Indonesia, Arfak Mts.

Etymology

The specific name refers to similarity with E. apertus.

Eniclases divaricatus (Pic, 1921)

Figs 6, 34

Trichalus (Trichalolus) divaricatus Pic, 1921: 10.

Material examined

Lectotype. Female, New Guinea, Humboldt Bay, Doherty lgt., coll. Pic (MHNP). Other material examined. 4 males, 2 females (BM0001–2, 9, 15–17), Indonesia, Irian Jaya, Sentani, Cyclops Mts., 300 m, Nov.–Dec. 1991; female (BM0057), Indonesia, Irian Jaya, Elelim, path to Apalapsili, 600 m, Nov.–Dec. 1991 (LMBC).

Diagnosis

Eniclases divaricatus is the only Papuan species with the bicolored elytra and flabellate male antennae (Figs 6, 34). Additionally, this species has the characteristic pattern of bright humeri and dark colored suture, which is shared only with some specimens of E. similis from the same locality (Figs 6, 13). The similarly colored females cannot be distinguished as their relative size of eyes is similar (Table 1).

Redescription

Male. Body length 8.8–9.7 mm. Head, thorax, legs, and humeri yellow to orange, antennae except basal part, abdomen, and most of elytra dark colored. Head with moderately large, hemispherically prominent eyes, their diameter 0.92–0.96 times minimum interocular distance Antennae flabellate, lamella of antennomere 3 slightly longer than the body of antennomere, other lamellae similar (Fig. 34). Pronotum transverse, with almost straight, slightly elevated lateral margins, elytra with quite strong straight secondary costae and regular dense cells (Fig. 6).

Distribution

Central North New Guinea.

Eniclases infuscatus sp. n.

Figs 11, 35

Material examined

Holotype. Female (BM0050), Indonesia, Irian Jaya, Bokondini, 1300 m, Nov.–Dec. 1991 (LMBC). Paratype. Female (BM0062), Indonesia, Irian Jaya, Elelim, path to Apalapsili, 600 m, Nov.–Dec. 1991 (LMBC).

Diagnosis

Eniclases infuscatus has a unique color pattern among western Papuan Eniclases. The upper part of the body is yellow to orange and only tips of elytra and the posterior part of the lateral margins are infuscate. This species partly resembles in the coloration E. tikapurensis, which is slender and pale colored (Figs 9, 11).

Description

Female. Body length 12.1 mm, robust. Head brown, abdomen, meso- and metathorax dark colored, most of legs similarly colored, only trochanters and basal half of femora light brown; pronotum and elytra yellow to orange, only tips and posterior margins of elytra infuscate, transition between dark and bright parts of elytra gradual (Fig. 11). Head with small hemispherically prominent eyes, their diameter 0.79–0.84 times minimum interocular distance. Antennae serrate, antennomeres 3–4 triangular, antennomeres 5–10 almost parallel-sided (Fig. 35). Pronotum transverse, with almost straight lateral margins, elytra with weak, regular secondary costae and regular, often transverse, small cells.

Distribution

New Guinea, Bokondini region.

Etymology

The specific name refers to blackened edge of the apical part of elytra.

Eniclases bicolor sp. n.

Figs 10, 31

Material examined

Holotype. Female (BM0046), Indonesia, Irian Jaya, Elelim, path to Apalapsili, 600 m, Nov.–Dec. 1991 (LMBC). Paratypes. 2 females (BM0045, 47), same locality data as the holotype (LMBC).

Diagnosis

Eniclases bicolor resembles some forms of E. similis from the Cyclops Mts. and differs in a higher contrast between light colored costae and dark cells in a transitional area between the dark and light parts of their elytra (Fig. 10). Additionally, the females of E. bicolor have small eyes (Table 1).

Description

Female. Body length 10.3 mm, robust. Head, basal part of antennae, pro- and mesothorax, basal half of elytra and legs yellow to light orange, apical half of antennae, metathorax, half of elytra and abdomen dark, transition between dark and bright parts abrupt(Fig. 10). Head with small, hemispherically prominent eyes, their diameter 0.71–0.74 times minimum interocular distance, antennae slender, antennomeres 3–5 serrate, antennomeres 6–10 almost parallel-sided (Fig. 31). Pronotum transverse, with almost straight lateral margins and prominent posterior angles, elytra with weak, but straight secondary costae and regular, often transverse, small cells.

Distribution

New Guinea, Elelim region.

Etymology

The specific name refers to the coloration of elytra.

Eniclases tikapurensis sp. n.

Figs 9, 41, 45

Material examined

Holotype. Male (BM0039), Indonesia, Irian Jaya, Yiwika, N of Wamena, 2000 m, Nov.–Dec. 1991 (LMBC). Paratypes. 3 males, 2 females (BM0040–44), Indonesia, Irian Jaya, Tikapura village, 2200 m, Nov.–Dec. 1991 (LMBC).

Diagnosis

Eniclases tikapurensis resembles in general appearance E. papuensis Bocak & Bocakova, 1991 from the Panai Lake area approximately 250 km west of Yiwika. Both species are characteristic in pale hue of the yellow upper part (Fig. 9) but differ in the relative size of eyes when E. tikapurensis has much larger eyes (Table 1).

Description

Male. Body length 9.5–11.1 mm, slender (Table 1). Head, antennae, thorax, legs except trochanters and bases of femora bark brown to black, pronotum and elytra except posterior edge pale yellow (Fig. 9). Head with large, hemispherically prominent eyes, their diameter 1.11–1.40 times minimum interocular distance. Antennae serrate, antennomere 3 triangular, with pointed apical process, antennomere 4 parallel-sided in most of its length, its process shorter, shape of other antennomeres similar (Fig. 41). Pronotum transverse, lateral margins with weak bulge in basal third and prominent posterior angles, elytra with weak but well developed, straight secondary costae and regular, mostly quadrate, small elytral cells.

Distribution

New Guinea, Upper Baliem Valley.

Etymology

The specific name refers to the type locality, the village Tikapura, north of Tagime.

Remark

The genetically distant population from Bokondini is a sister to E. tikapurensis, but does not differ in any morphological character. As their delimitation would be based only on mtDNA sequence and further information on nuclear markers and geographical distribution is needed for these two putative cryptic species, we postpone the formal description of the population from Bokondini. The sequenced specimens representing the Bokondini population are designated as Eniclases sp. A. in Fig. 1.

Eniclases brancuccii sp. n.

Figs 12, 33, 46

Material examined

Holotype. Male (BM0006), Indonesia, Irian Jaya, Sentani, Cyclops Mts., 300 m, Nov.–Dec. 1991 (LMBC). Paratypes. 3 females (BM0005, 0007, 00010), the same data as the holotype (LMBC).

Diagnosis

Eniclases brancuccii resembles in color pattern three species: E. elelimensis, E. niger, and E. variabilis. The last two of them are polymorphic and only some individuals share the color pattern with E. brancuccii (Figs 12, 16, 24–25, 28–29). Eniclases niger differs in the large eyes and almost parallel-sided antennomeres 3–10 (Table 1, Figs 33–40). Eniclases elelimensis and E. variabilis have similar antennae (Figs 38–39) and slightly larger eyes compared to E. brancuccii, but reliable identification of these species can only be based on the DNA sequences. To the best of our knowledge, the similarly colored forms of these species do not occur together in a single locality.

Description

Male. Body length 7.6–8.0 mm, robust, head, antennae, thorax, elytra, and abdomen dark brown to black, pronotum, scutellum, and basal parts of femora yellow to orange, apical parts of femora, tibiae, and tarsi dark brown. Head with small, hemispherically prominent eyes, their diameter equals minimum interocular distance, antennae flat, slender, acutely serrate, apical process of antennomere 3 about half length of its body, further antennomeres similar in shape. Pronotum transverse, with apparent bulge in basal third, almost parallel-sided between bulge and posterior angles, elytra with weak but regular secondary costae and mostly regular, quadrate, small cells.

Distribution

New Guinea, Cyclops Mountains.

Etymology

The specific name ‘brancuccii’ is proposed in honor of the late Michel Brancucci, a specialist in Dytiscidae and Cantharidae.

Eniclases elelimensis sp. n.

Figs 19, 39

Material examined

Holotype. Female (BM0056), Indonesia, Irian Jaya, Elelim, path to Apalapsili, 600 m, Nov.–Dec. 1991 (LMBC). Paratypes. 2 females (BM0051–52), the same data as the holotype (LMBC).

Diagnosis

Eniclases elelimensis differs from E. brancuccii in slightly smaller eyes (Table 1), but no morphological difference has been found to discriminate this species from some similarly colored individuals of E. variabilis (Figs 24–25, 28–29) and morphology based identification can reliably assign similarly colored specimens only to the clade of E. variabilis and related species.

Description

Male. Body length 6.9–8.1 mm, robust, head, antennae, thorax, elytra, and abdomen dark brown to black, pronotum and femora yellow to light orange, apical part of femora, tibiae, and tarsi dark brown. Head with small, hemispherically prominent eyes, their diameter equals minimum interocular distance, antennae flat, slender, acutely serrate, apical process of antennomere 3 about half length of its body, further antennomeres similar in shape. Pronotum transverse, with apparent bulge in basal third, almost parallel-sided between bulge and posterior angles, elytra with weak but regular secondary costae and mostly regular, quadrate, small cells.

Distribution

New Guinea, Elelim region.

Etymology

The specific name refers to the type locality.

Remark

E. versicolor Kleine, 1926 from an unspecified locality in New Guinea was studied (Holotype, male, ‘Neuguinea, Coll. Kraatz’ deposited in DEIM). Eniclases versicolor is similar in general appearance, but differs in large male eyes. The female specimen of E. versicolor from the Jayapura district cited by Bocak and Bocakova (1991) might be conspecific with E. variabilis or E. elelimensis and these species might have allopatric distribution.

Eniclases bokondinensis sp. n.

Figs 17, 18, 32

Material examined

Holotype. Female (BM0095), Indonesia, Irian Jaya, Bokondini, 1900 m, Nov.–Dec. 1991 (LMBC). Paratypes. 2 females (BM0092, 94), the same data as the holotype (LMBC).

Diagnosis

Eniclases bokondinensis has a characteristic combination of the black pronotum and light yellow elytra with dark colored apex (Figs 1718). The shape of antennae and the relative size of eyes are similar to those of E. variabilis and E. elelimensis (Table 1, Figs 32, 38, 39). The co-mimics of E. bokondinensis are large-bodied and this species has a larger body than its closest relatives (Fig. 1, Table 1).

Description

Female. Body length 9.2 mm, head, antennae, thorax, and abdomen dark brown to black, elytra pale yellow in humeral half, gradually infuscate to apex (Figs 1718), trochanters and basal parts of femora brown, rest of legs black. Head with small, hemispherically prominent eyes, their diameter 0.72–0.82 interocular distance, antennae flat, slender, acutely serrate, apical process of antennomere 3 about half length of its body, further antennomeres similar in shape. Pronotum transverse, with weak bulge in basal third, elytra with weak but regular secondary costae and mostly regular, subquadrate cells (Figs 1718).

Distribution

New Guinea, Bokondini region.

Etymology

The specific name refers to the type locality.

Eniclases variabilis sp. n.

Figs 20–29, 38, 47

Material examined

Holotype. Male (BM0054), Indonesia, Irian Jaya, Elelim, path to Apalapsili, 600 m, Nov.–Dec. 1991 (LMBC). Paratypes. 1 male, 3 females (BM0048–49, 53, 55), same data as the holotype; 7 males, 8 females (BM0025–32, 34–36, 86, 88, 90–91), Indonesia, Irian Jaya, Bokondini, 1300 m, Nov.–Dec. 1991; 4 males, female (BM0063–67), Indonesia, Irian Jaya, Dombomi, Lower Pass valley, 1200 m; 2 males (BM0008, 12), Indonesia, Irian Jaya, Sentani, Cyclops Mts., 300 m, Nov.–Dec. 1991 (LMBC).

Diagnosis

Eniclases variabilis is a widespread, highly polymorphic species and resembles in general appearance several distinct, sympatric mimetic types. The color forms are illustrated in Figs 20–29. Similarly colored individuals of E. niger differ in the large eyes (1.17–1.28 times larger than eye distance) and acutely serrate antennae of E. variabilis (Figs 38, 40). The similarly colored individuals of E. elelimensis can only be identified using the DNA data.

Description

Male. Body length 6.6–8.2 mm. Head, antennae, thorax, and abdomen dark black, elytra orange yellow in humeral third, rest of elytra black, transitional zone between bright and dark part of elytra is narrow (Figs 20–29, 38), trochanters brown, rest of legs black. Head with small, hemispherically prominent eyes, their diameter 0.83–0.95 interocular distance, antennae flat, slender, acutely serrate, apical process of antennomere 3 about half length of its body, further antennomeres similar in shape. Pronotum transverse, with weak bulge in basal third, elytra with weak secondary costae and mostly regular, subquadrate cells (Fig. 20).

Distribution

Central North New Guinea.

Etymology

The specific name refers to exceptional variability in coloration.

Remark

Eniclases variabilis can be differently colored in various localities. Generally, the dark colored specimens (Fig. 27) were found in higher elevations in Bokondini (1300 m) and Dombomi (1200 m) and bright colored individuals (Figs 22–23, 26) were collected in lower elevations in Elelim (600 m) and Cyclops Mts. (300 m) (Fig. 3).

Eniclases pseudoluteolus sp. n.

Figs 36, 48

Material examined

Holotype. Male (BM0084), Indonesia, Irian Jaya, Maibri vill., Arfak Mts., 1600 m, Nov.–Dec. 1991 (LMBC).

Diagnosis

Eniclases pseudoluteolus belongs to the group of uniformly yellow species which additionally includes E. luteolus Waterhouse, 1878, E. nigriceps Bocak & Bocakova, 1991, E. fuscicornis Bocak & Bocakova, 1991, and E. pallidus Bocak & Bocakova, 1991. Most of them have large eyes (Bocak and Bocakova 1991). Two species, E. robustus Bocak & Bocakova, 1991 and E. proximus Bocak & Bocakova, 1991 have the similar size of eyes with E. pseudoluteolus. These species differ from E. pseudoluteolus in the slender antennomere 3 and light colored thorax (antennomere of E. pseudoluteolus as in Fig. 36). The similar species occur only in Eastern New Guinea in the vicinity of Wau and Mt. Hagen (Fig. 2).

Description

Male. Body length 9.3 mm, head, apical antennomeres, thorax, and abdomen dark brown to black, pronotum and elytra yellow, trochanters brown, rest of legs black. Head with small, hemispherically prominent eyes, their diameter 0.90 interocular distance, antennae flat, acutely serrate, antennomere 3 triangular, wide, its apical process short (Fig. 36), further antennomeres similar in shape, becoming slenderer to apex of antennae. Pronotum transverse, with weak bulge in basal third, elytra with weak secondary costae and mostly regular, subquadrate cells.

Distribution

New Guinea, Arfak Mts.

Etymology

The specific name refers to the similarity with E. luteolus.

Eniclases niger sp. n.

Figs 15, 16, 40, 49

Material examined

Holotype. Male (BM0059), Indonesia, Irian Jaya, Elelim, path to Apalapsili, 600 m (LMBC). Paratypes. 3 males (BM0058, 60–61), same data as the holotype; 2 males, 1 female (BM0033, 87, 89), Indonesia, Irian Jaya, Bokondini, 1300 m, Nov.–Dec. 1991 (LMBC).

Diagnosis

Eniclases niger is polymorphic in coloration and can be uniformly black or can have the brightly colored pronotum and scutellum (Figs 15–16). The latter specimens resemble some individuals of E. variabilis sp. n. when they occur in the same locality. These two species differ in the relative size of eyes (Table 1).

Description

Male. Body length 9.2–11.6 mm (Table 1), uniformly dark colored, only trochanters and bases of femora brown (Fig. 16). Head with large hemispherically prominent eyes, their diameter 1.17–1.28 times distance between eyes. Antennae serrate, with apical process about a third of antennomere stem (Fig. 40). Pronotum transverse, with almost straight lateral margins, lateral bulge inconspicuous, elytra with weak secondary costae and in some parts with irregular cells.

Distribution

New Guinea, Central Mountains.

Etymology

The specific name refers to body coloration.

Remark

Eniclases niger has two forms (Figs 15, 16). The specimens from Bokondini have the brightly colored pronotum and scutellum (Fig. 16), the individuals from other localities are uniformly dark colored (Fig. 15).

Eniclases similis Bocak & Bocakova, 1991

Figs 13, 14, 37, 50

Eniclases similis Bocak & Bocakova, 1991: 210.

Material examined

10 males, 2 females (BM0003–4, 11, 13–14, 19–24, 37), Indonesia, Irian Jaya, Sentani, Cyclops Mts., 300 m, Nov.–Dec. 1991 (LMBC).

Diagnosis

Eniclases similis differs from the similarly colored individuals of E. divaricatus in large eyes (Table 1) and from E. bicolor in the gradual transition between the dark and bright parts of elytra (Figs 10, 13–14). Despite distant relationships between E. similis and E. bicolor (Fig. 1) there is no clearly defined morphological character available for their identification.

Redescription

Male. Body length 7.5–9.7 mm, pronotum, humeral part of elytra, and legs yellow to orange, metathorax, abdomen, apical part of elytra, and sometimes antennae dark brown to black (Figs 13–14). Head with moderately large eyes, their diameter 1.17–1.28 times longer than interocular distance, antennae slender, acutely serrate, with apical process of antennomere 3 about third of antennomere stem. Pronotum transverse, with weak lateral bulge in posterior third, elytra with weak, but straight secondary costae and regular, subquadrate cells.

Distribution

New Guinea Highlands.

Remark

Eniclases similis has two extreme forms in the extent of the bright part of elytra (Figs 13, 14) and transitional forms between these. The similar color patterns are present in sympatrically occurring species of Trichalus. The species is identified as E. similis due to the absence of any morphological difference when compared to the holotype of E. similis. The type locality of E. similis is Karimui in Eastern New Guinea and identity of the species needs further investigation.

Discussion

The taxonomy of Eniclases has been based solely on the morphological species concept, which depends on the presence of identifiable diagnostic traits (Kleine 1926, Bocak and Bocakova 1991). Such species delimitation is difficult as these beetles are morphologically similar and their diagnostic characters are limited to the relative size of eyes, the shape of antennomeres and body coloration (Figs 642). Net-winged beetles are unpalatable and the advergence to the similar body shape and coloration has been shown in other metriorrhynchines (Bocak and Yagi 2010). Therefore, we suggest that Eniclases are similarly selected to resemble sympatrically occurring net-winged beetles. The high degree of resemblance can be demonstrated by similar color patterns of sympatrically occurring and distantly related E. similis and E. divaricatus (Figs 1, 6, 13). Already the morphology based revision of Eniclases by Bocak and Bocakova (1991) pointed to a high degree of uncertainty in species delimitation and suggested that some species, such as size and color variable E. luteolus, probably represent a group of morphologically similar species whose status cannot be resolved morphologically.

Morphological and DNA sequence diversification was investigated in a clade of 14 Eniclases net-winged beetles from the western part of New Guinea, delimited as separate species using mtDNA haplotypes, morphological characters and biogeography (DeSalle et al. 2005). We hypothesize that these clusters of individuals represent biological species and are reproductively isolated. Several species are unique in coloration, the shape of antennae, and size of eyes and they were inferred consensually as separate species by all DNA-based analyses: E. pseudoapertus (the closest interspecific match at 10.62%), E. divaricatus (11.02%), E. apertus (6.20%), and E. pseudoluteolus (9.26%). The genetic distances between these species and their closest relatives are higher than generally accepted intraspecific genetic variation (Hebert et al. 2003, Meier et al. 2006), the posterior delimitation probabilities for the branches leading to these species were all over 90% except the branch of E. divaricatus (69%).

A much higher degree of uncertainty was found in pairs of closely related terminals, which were refused as putative species by various methods (Fig. 1). We can delimit broadly defined species as suggested by the barcoding threshold (Fig. 1) and consider internal clades as forms without formal taxonomic delimitation or we can delimit each clade as a separate species. The first approach makes morphological definitions difficult, as at least in some cases, we would have to merge in a single species some individuals with clearly different morphology. Additionally, there are indications that the narrowly defined species are independently selected for different mimetic patterns in various localities and the color patterns supposedly support reproductive isolation (Bocak and Yagi 2010). Therefore, we prefer to split the closely related clades into separate species when they are supported by clear morphological difference (e.g. size of eyes as E. similis and E. niger or E. infuscatus and E. bicolor); they occur allopatrically, belong to different mimicry rings and simultaneously differ in the sequence of cox1 mtDNA (e.g., E. bokondinensis and E. elelimensis). In one case, the sister clades are defined by divergent DNA sequence and distribution, but no morphological divergence was identified (E. tikapurensis + Eniclases sp. A, Fig. 1). Below, we discuss reasons for species delimitation in detail.

The pair of E. infuscatus and E. bicolor represents sister species which differ in coloration (Figs 10–11) and the relative size of female eyes (Table 1). They were marginally inferred as separate species using the bPTP model (PP 0.47 and 0.52 versus 0.48 for the branch merging these two entities), but clustered together as a putative single species using the distance method (0.09–1.08% distance). Therefore, based on coloration, size of eyes, and genetic difference, we consider these subclades as separate species.

Similarly, E. bokondinensis and E. elelimensis represent genetically close species with different color patterns (Figs 1719). The bPTP model merged them in a single putative species and the distance method merged both of them with E. variabilis. Both species are known from the Central Mts.; E. bokondinensis from the mountain region north of Bokondini (~1900 m a. s. l.) and E. elelimensis from Elelim (600 m a. s. l., Fig. 2) about 80 km away. We suppose that in this case the memberships in different mimetic rings play a role in genetic differentiation between these two species (Bocak and Yagi 2010). Species with the E. elelimensis color pattern occur in low elevations and have never been collected in high mountains in the Bokondini area. This pattern was recorded as high as 1300 m a. s. l. in, but in a biotope different from mountain forests where E. bokondinensis occurs. We suppose that these species occur allopatrically despite the proximity of localities and additionally the colonization of high elevations could support their reproductive isolation (Toussaint et al. 2014).

Eniclases variabilis is a sister to the E. bokondinensis + E. elelimensis clade and was identified as a separate entity using the bPTP model. Other methods merged this species with its sister clade (genetic distance, Fig. 1) or a part of it (morphology, see taxonomy section). The reproductive isolation of E. variabilis is supported by the sympatric occurrence of genetically differentiated E. variabilis and E. elelimensis in Elelim. E. variabilis is extremely polymorphic (Figs 20–29) and resembles black colored E. niger, brightly colored E. similis and E. elelimensis types, and one color pattern is unique in Eniclases and resembles other net-winged beetles (Figs 20–29). No genetic differentiation has been found in their mtDNA and we consider all color forms as a single species.

The clade of E. similis + E. niger was merged into a single putative species by all DNA based analyses, but they can be identified by morphology and color patterns. Eniclases niger has large eyes (diameter/distance ratio 1.17–1.28) and E. similis has smaller eyes (1.02–1.15). We suppose that the daytime or evening, eventually night activity of respective species might be the reason for observed morphological differentiation. Additionally, these species belong to different mimetic complexes. E. niger is black colored (Fig. 15) or has pronotum bright and elytra completely black (Fig. 16) and E. similis is brightly colored. The bright patterns are similar to those of E. divaricatus and E. bicolor (Figs 6, 10, 13–14) and were recorded also in several Trichalus spp. in the respective localities (unpublished data). These sister-species, E. similis and E. niger, have not been collected sympatrically. E. niger occurs in lower mountain localities in the Central Mts. and E. similis in low elevations of the Cyclops Mts.

The clade of E. tikapurensis and Eniclases sp. A split in two subclades, which do not differ in morphology, but they are genetically distant. The levels of DNA distances between E. tikapurensis and Eniclases sp. A do not agree with the geographical distance of respective localities, when higher differentiation was found between populations from Bokondini and Tikapura (13 km apart) than between localities lying at the rim of the Baliem valley (Yiwika and Tikapura, 35 km apart). The observed genetic differentiation surpasses some cases when separate morphologically divergent sister species are delimited. Due to limited information we postpone formal description of the putative species from Bokondini.

Eniclases are variable in coloration (Figs 629), but the differences in coloration regardless of their depth do not necessarily mean that they can be used for delimitation of reproductively isolated lineages. We found intraspecific color variability in three species: E. similis, E. niger, and E. variabilis (Figs 13–14, 15–16, 20–29). We ascribe the polymorphism to the advergence to the most common models of other net-winged beetles as described by Bocak and Yagi (2010). Conversely, the color patterns can be used in several species for their delimitation, e.g. in E. bicolor and very characteristic E. bokondinensis. Molecular data or further morphological characters are generally needed to investigate color polymorphism and to support species limits. Uncorrelated morphological and genetic diversity has been reported in lyponiine lycids (Li et al. 2015) and other beetle families (e.g. Goldberg et al. 2012) and points to the necessary evaluation of all evidence when species are delimited (Jorger and Schrodl 2013). The cox1 mtDNA fragment alone cannot provide sufficient information on the diversification process, but even these limited data suggest that closely related lineages can develop distant mimetic patterns and start further morphological differentiation, for example in the relative size of eyes. The dated phylogeny suggests that all closely related sister species differentiated in the last million years (Fig. 2).

Acknowledgements

We are sincerely grateful to R. Bilkova for technical assistance and to colleagues who provided us with the material for this study. Timothy C. Bray critically read the MS prior to its submission.

This work was supported by the Czech Science Foundation (P506/11/1757) and the Palacky University (PrF2015-018).

References

  • Baselga A, Gomez-Rodriguez C, Novoa F, Vogler AP (2013) Rare failures of DNA barcodes to separate morphologically distinct species in a biodiversity survey of Iberian leaf beetles. PLoS ONE 8(9): e74854. doi: 10.1371/journal.pone.0074854
  • Bocak L (2002) Generic revision and phylogenetic analysis of the Metriorrhynchinae (Coleoptera, Lycidae). European Journal of Entomology 99: 315–351. doi: 10.14411/eje.2002.043
  • Bocak L, Bocakova M (1991) Revision of the genus Eniclases Waterhouse, 1879 (Coleoptera, Lycidae, Metriorrhynchinae). Mitteilungen der Münchener Entomologischen Gesellschaft 81: 203–226.
  • Bocak L, Bocakova M (2008) Phylogeny and classification of the family Lycidae (Insecta, Coleoptera). Annales Zoologici 58: 695–720. doi: 10.3161/000345408x396639
  • Bocak L, Yagi T (2010) Evolution of mimicry patterns in Metriorrhynchus (Coleoptera, Lycidae): The history of dispersal and speciation in Southeast Asia. Evolution 64: 39–52. doi: 10.1111/j.1558-5646.2009.00812.x
  • Brower AVZ (1994) Rapid morphological radiation and convergence among races of the butterfly Heliconius erato inferred from patterns of mitochondrial-DNA evolution. Proceedings of the National Academy of Sciences of the United States of America 91: 6491–6495. doi: 10.1073/pnas.91.14.6491
  • Calder A (1998) Zoological Catalogue of Australia – Coleoptera: Elateroidea. CSIRO Publishing, Melbourne.
  • Darriba D, Taboada GL, Doallo R, Posada D (2012) jModelTest 2: more models, new heuristics and parallel computing. Nature Methods 9(8): 772. doi: 10.1038/nmeth.2109
  • DeSalle R, Egan MG, Siddall M (2005) The unholy trinity: taxonomy, species delimitation and DNA barcoding. Philosophical Transactions of the Royal Society B-Biological Sciences 360: 1905–1916. doi: 10.1098/rstb.2005.1722
  • Drummond AJ, Suchard MA, Xie D, Rambaut A (2012) Bayesian Phylogenetics with BEAUti and the BEAST 1.7. Molecular Biology and Evolution 29: 1969–1973. doi: 10.1093/molbev/mss075
  • Goldberg CS, Tank DC, Uribe-Convers S, Bosworth WR, Marx HE, Waits LP (2012) Species designation of the Bruneau Dune tiger beetle (Cicindela waynei) is supported by phylogenetic analysis of mitochondrial DNA sequence data. Conservation Genetics 13: 373–380. doi: 10.1007/s10592-011-0295-9
  • Hebert PDN, Ratnasingham S, Dewaard JR (2003) Barcoding animal life: cytochrome c oxidase subunit 1 divergences among closely related species. Proceedings of the Royal Society B-Biological Sciences 270: S96–S99. doi: 10.1098/rsbl.2003.00
  • Hill KC, Hall R (2002) Mesozoic-Cainozoic evolution of Australia’s New Guinea margin in a West Pacific context. In: Hillis R, Müller RD (Eds) Defining Australia: The Australian Plate as Part of the Planet Earth.Geological Society of America Special Paper/Geological Society of Australia Special Paper 372: 265–290. doi: 10.1130/0-8137-2372-8.265
  • Huson DH, Richter DC, Rausch C, Dezulian T, Franz M, Rupp R (2007) Dendroscope: An interactive viewer for large phylogenetic trees. BMC Bioinformatics 8: 6. doi: 10.1186/1471-2105-8-460
  • Jorger KM, Schrodl M (2013) How to describe a cryptic species? Practical challenges of molecular taxonomy. Frontiers in Zoology 10: 27. doi: 10.1186/1742-9994-10-59
  • Kleine R (1926) ColeopteraLycidae. Nova Guinea 15: 91–195.
  • Li Y, Gunter N, Hong P, Bocak L (2015) DNA-based species delimitation separates highly divergent populations within morphologically coherent clades of poorly dispersing beetles. Zoological Journal of the Linnean Society 175: 59–72. doi: 10.1111/zoj.12262
  • Meier R, Shiyang K, Vaidya G, Ng PKL (2006) DNA barcoding and taxonomy in Diptera: a tale of high intraspecific variability and low identification success. Systematic Biology 55: 715–728. doi: 10.1080/10635150600969864
  • Meyer CP, Paulay G (2005) DNA barcoding: Error rates based on comprehensive sampling. PLoS Biology 3: 2229–2238. doi: 10.1371/journal.pbio.0030422
  • Papadopoulou A, Anastasiou I, Vogler AP (2010) Revisiting the Insect Mitochondrial Molecular Clock: The Mid-Aegean Trench Calibration. Molecular Biology and Evolution 27: 1659–1672. doi: 10.1093/molbev/msq051
  • Pic M (1921) Contribution à l’étude des Lycides. L’Echange 406: 9–12.
  • Riedel A, Sagata K, Surbakti S, Tanzler R, Balke M (2013) One hundred and one new species of Trigonopterus weevils from New Guinea. ZooKeys 280: 1–150. doi: 10.3897/zookeys.280.3906
  • Sklenarova K, Chesters D, Bocak L (2013) Phylogeography of poorly dispersing net-winged beetles: a role of drifting India in the origin of Afrotropical and Oriental fauna. PLoS ONE 8(6): e67957. doi: 10.1371/journal.pone.0067957
  • Sklenarova K, Kubecek V, Bocak L (2014) Subtribal classification of Metriorrhynchini (Insecta, Coleoptera, Lycidae): an integrative approach using molecular phylogeny and morphology of adults and larvae. Arthropod Systematics & Phylogeny 72: 37–54.
  • Stamatakis A (2006) RAxML-VI-HPC: Maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 22: 2688–2690. doi: 10.1093/bioinformatics/btl446
  • Stamatakis A, Hoover P, Rougemont J (2008) A Rapid Bootstrap Algorithm for the RAxML Web Servers. Systematic Biology 57: 758–771. doi: 10.1080/10635150802429642
  • Tamura K, Stecher G, Peterson D, Filipski A, Kumar S (2013) MEGA6: Molecular Evolutionary Genetics Analysis Version 6.0. Molecular Biology and Evolution 30: 2725–2729. doi: 10.1093/molbev/mst197
  • Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG (1997) The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Research 25: 4876–4882. doi: 10.1093/nar/25.24.4876
  • Toussaint EFA, Hall R, Monaghan MT, Sagata K, Ibalim S, Shaverdo HV, Vogler AP, Pons J, Balke M (2014) The towering orogeny of New Guinea as a trigger for arthropod megadiversity. Nature Communications 5: 10. doi: 10.1038/ncomms5001
  • Waterhouse CO (1879) Illustration of the Typical Specimens of Coleoptera in the Collection of the British Museum. Part I.-Lycidae. British Museum, London.
  • Zhang JJ, Kapli P, Pavlidis P, Stamatakis A (2013) A general species delimitation method with applications to phylogenetic placements. Bioinformatics 29: 2869–2876. doi: 10.1093/bioinformatics/btt499
login to comment