Phylogenetic Signals from Nepomorpha (Insecta: Hemiptera: Heteroptera) Mouthparts: Stylets Bundle, Sense Organs, and Labial Segments

The present study is a cladistic analysis of morphological characters focusing on the file of the mandible, the apices of the maxillae, the rupturing device on the maxillae, the internal structures of the mouthparts, and the external morphology of the labial segments as well as the distribution of labial sensilla in true water bugs (Hemiptera: Heteroptera, infraorder Nepomorpha). The study is based on data referring to sixty-two species representing all nepomorphan families (Heteroptera), together with one outgroup species representing the infraorders Gerromorpha (Mesoveliidae). The morphological data matrix consists of forty-eight characters. The present hypothesis supports the monophyly of the Nepomorpha and the monophyly of all families. The new modification in the systematic classification has been proposed: ((Nepidae + Belostomatidae), (Diaprepocoridae + Corixidae + Micronectidae), (Ochteridae + Gelastocoridae), Aphelocheiridae, Potamocoridae, Naucoridae, Notonectidae, and (Pleidae + Helotrephidae)).

As far as the approach to cladistic relationships among infraorders of the Heteroptera is concerned, several various hypotheses have been proposed with respect to the systematic position of the Nepomorpha as well as other infraorders (generally without considering the Aradimorpha) and discussed in phylogenetic studies.
According to Schuh [6] the Enicocephalomorpha (first branch) is a basal, sister group to all remaining infraorders of the Heteroptera. In the following branches the Dipsocoromorpha, Gerromorpha, Nepomorpha, Leptopodomorpha, Cimicomorpha, and Pentatomomorpha have been placed. In this arrangement the Gerromorpha is the sister group to the Nepomorpha.
Similarly, on the basis of morphological evidence,Štys [7,8] placed the Enicocephalomorpha in the most basal clade of the Heteroptera. In Zrzavy's [9] system of relationships it was indicated that the Enicocephalomorpha was the sister group of Dipsocoromorpha + Gerromorpha and together formed a basal heteropteran clade in relation to the unresolved relationships (polytomy) among the Nepomorpha, Leptopodomorpha, and Cimicomorpha + Pentatomomorpha.
Furthermore, on the basis of morphological characters, Mahner [10] proposed a hypothesis that the Nepomorpha (Cryptocerata) should be placed as the basal, sister taxon to the remaining Heteroptera which also coincided with the hypothesis of Shcherbakov and Popov [11], based on fossil morphological evidence, although in both cases unresolved relationships were notated among the remaining infraorders.
Wheeler et al. [12] generally reached a substantial congruence between the molecular data and most of the morphological data used by Schuh [6] in the system of classification of heteropteran infraorders, even though a distinct result 2 The Scientific World Journal seemed to be the establishing of the sister clade Nepomorpha + Leptopodomorpha (as in Figure 6 in [12]). In the classifications of infraorders based on characters of male genitalia Yang [13] pointed out that the Enicocephalomorpha was the sister group to the remaining Heteroptera; however, in the arrangement such as the Enicocephalomorpha + (Leptopodomorpha + Cimicomorpha + Pentatomomorpha) + (Dipsocoromorpha + (Nepomorpha + Gerromorpha)), the above mentioned groups were presented as three unresolved branches. Recently, in the infraordinal relationships based on whole sequences of 18S rDNA whose alignment was modified by the secondary structure of rRNA, Xie et al. [14] obtained results featuring single branches for the Enicocephalomorpha, Nepomorpha, Leptopodomorpha, and two clades: Gerromorpha + Dipsocoromorpha and Cimicomorpha + Pentatomomorpha. Besides, using 64 morphological characters and DNA sequence data from the mitochondrial genes encoding COI+II and 16S rRNA and the nuclear gene encoding 28S rRNA, Damgaard [15] established the relationship of (Enicocephalomorpha + (Dipsocoromorpha + (Gerromorpha + Nepomorpha))). On the basis of multiple genes in many species of the heteropteran infraorders, Li et al. [16] revealed that the Nepomorpha was the most basal group.
Several different systematic position of the Nepomorpha within the Heteroptera have been indicated and essentially, only in three studies, the Gerromorpha has been estimated as the outgroup (sister clade) to the nepomorphans [6,12,15].
In historical views, various hypotheses have been proposed with respect to the relationships within the taxa of the Nepomorpha. In the first evaluations of the relationships among the true water bugs [20], the Corixidae were considered to be primitive and treated as a sister group of the remaining families. Later, China [21] proposed a scheme of the relationships among nepomorphan families in which the Ochteridae were treated as relatively the most primitive group based on the possession of ocelli and a respiratory system typical of terrestrial bugs. Comparative studies of the mouthparts [19,22,23], the egg structures, and other characters of embryology of the Heteroptera [24] supported the hypothesis that had been proposed by China [21]. Furthermore, according to several authors such as Popov [25], Rieger [26], Mahner [10], and Hebsgaard et al. [27], the Belostomatidae and Nepidae (Nepoidea) were introduced at the basal position and estimated as a sister group of the remaining nepomorphan families. However, their scenario of relationships was essentially different with respect to other nepomorphan families. Popov [25], Mahner [10], and Hebsgaard et al. [27] placed the Corixidae (Corixoidea) as the second group (second branch) in the arrangement of relationships system. In Rieger's [26] system the clade Ochteridae + Gelastocoridae (Ochteroidea) was shown in the second branch, while the Corixidae were presented as the third branch. A similar concept of the relationships among the Potamocoridae, Naucoridae, and Aphelocheiridae could be seen in systems developed by Popov [25] and Mahner [10]. Rieger [26] indicated the clade Naucoridae + Potamocoridae; however, Hebsgaard et al. [27] indicated the clade Aphelocheiridae + Potamocoridae and ranked it as a new superfamily, Aphelocheiroidea; furthermore, they placed the Ochteridae and Gelastocoridae (Ochteroidea) in a new position, that is, as a branch under the Naucoridae (Naucoroidea). Popov [25], Rieger [26], Mahner [10], and Hebsgaard et al. [27] generally agreed in the classification and relationships of the Notonectidae, Pleidae, and Helotrephidae. The new relationships of some families of the Nepomorpha postulated by Hua et al. [28] are interesting due to the specific location of the Pleidae. According to these authors, the Pleidae derive from the Nepomorpha as a new heteropteran infraorder; the Plemorpha and the monophyletic infraorder Nepomorpha consist of five superfamilies with the following relationships: (Corixoidea + ((Naucoroidea + Notonectoidea) (Ochteroidea + Nepoidea)). Nevertheless, not all families which were recognized within the Nepomorpha in the study by Hua et al. [28] have been analysed so far. Moreover, other relationships of super(families) of the Nepomorpha based on four Hox genes have been indicated by Li et al. [29]. According to their study, the most basal lineage is the Ochteroidea, whereas the Notonectoidea include only the Notonectidae and form a new sister relationships with (Pleoidea + Naucoroidea) and (Nepoidea + Corixoidea).
Members of the true bug group of the Nepomorpha have attracted the attention of researchers by displaying a variety of body structure forms and lifestyles. All of these subjects have been widely discussed in literature. However, the number of papers dealing with nepomorphan (sub)family-level relationships still remains relatively small and the phylogenetic affinities of some family groups still require clarification.
The water bugs classified as the Nepomorpha include about 2000 species worldwide [30], and as for the composition of this infraorder, there are 13 families. In the classification ofŠtys and Jansson [31] 11 families of the Nepomorpha were distinguished. Two subfamilies (Diaprepocorinae and Micronectinae) that belonged to the Corixidae were elevated to the rank family level of the Micronectidae and Diaprepocoridae by Nieser [32] and have been accepted at that position by other researchers [33][34][35][36][37].
In order to meet the requirements of the many ways of life adopted by members of the Nepomorpha, the morphology of the species in various families displays a great variety of modifications [25,38]. The representatives of most families live in water (aquatic bugs), except for the gelastocorids and ochterids, which occupy habitats at the water's edge [3,25,30,[39][40][41][42], like the Saldidae of the Leptopodomorpha. Those that remain submerged include fast swimmers inhabiting the open water, including the corixids and notonectids [38,[43][44][45] and slow-moving benthic species that breathe through long respiratory siphons, such as many nepids and belostomatids [25]. Morphologically, the group is characterized by the shortness of the antenna that is typically concealed, either partly or entirely, by the eyes [25,30,46] and families of most species can be identified immediately on the basis of size and general body shape. The general morphology is usually similar throughout the larval stages, and the family affiliation of the first instar larvae can immediately be recognized after only a cursory examination in almost all cases [38,47]. However, the general morphology of corixids differs in several ways from that of other groups of true bugs [25,48].
The Scientific World Journal 3 Some of the Nepomorpha species are mostly predators whereas most corixids are plant feeding; however, there are several species which prefer feeding on animals or a mixed type of nutrition [43][44][45][49][50][51].
So far, a significant range of various studies have been conducted with regard to the nepomorphans. The most comprehensive papers on the subject have been written by China [21] on general biology and morphology of water bugs; Popham [52] on the respiration of aquatic bugs; Cobben [19,24] on embryology and eggs, male genitalia, and mouthpart structures of the Heteroptera; Parsons [22,23,48] on triturating device, salivary pump, thorax, and labial skeleton; Popov [25] on general morphology and fossils study; Rieger [26] on the structures of the head and prothorax of Ochterus; Cassis and Silveira [53] on morphology and interrelationships in the Gelastocoridae (Nerthrinae). Also, the principal work of Mahner [10] included a number of comprehensive useful data on morphology within the Nepomorpha and their phylogeny. The first studies combining the morphological and genetic data in the Nepomorpha were conducted by Hebsgaard et al. [27] and by Hua et al. [28], who studied phylogenetic relationships based on the genomes. Generally, many other studies focusing on various fields of morphology and biology of the nepomorphans have been conducted by researchers.
Comprehensive studies of characters of the labium in the nepomorphans and in several individual species of the Corixoidea were conducted by Parsons [22,48] and previously by Griffith [54], Bentwitz [55], and Puchkova [56]. The papers of the above mentioned authors generally presented researches focusing strictly on the labium.
The studies by , 2014 in press) have also provided a number of useful new observations on the morphology of maxillae and mandibles, labial sensilla and labial segments within the Nepomorpha, which can be considered from the phylogenetical perspective of this infraorder.
The results achieved in these studies allowed establishing many new features in comparison to previous researches (the present compilation of data is meant to provide a summarizing description of characters). The use of the available data in the present analysis of the characters is justified by the fact that in the last combined phylogenetic analysis (using morphological and molecular data) conducted by Hebsgaard et al. [27] only four characters of the labium were incorporated. Moreover, the labial sensilla as well as the maxillary and mandibular structures were not previously estimated in the cladistic analysis of the Nepomorpha. The evaluation of all these characters was conducted only on the basal ground plan of assumptions ( [57][58][59]). For this reason, the recent new descriptions of new characters of the mouthparts have provided an opportunity for reassessing the phylogenetic relationships within nepomorphan groups through conducting new analyses based on the available new data.
In doing so, an opportunity has been taken to reexamine the relationships within the Nepomorpha and to make a comparison between the concept presented by Hebsgaard et al. [27], based on relationships achieved from the morphological data and also from the final combined data (morphological and genetic) and the currently available new strict consensus on the phylogenetical tree. The principal approach of the present cladistic analyses with respect to the concept of Hebsgaard et al. [27] has been possible due to the use of the same methodology in the analyses focusing on a great number of the same species or genera. Furthermore, in the studies of the Nepomorpha conducted by Hua et al. [28] as well as Li et al. [29] different relationships of super(families) were presented based on various molecular data using a cladistic analysis, so the comparison of their results with the present morphological data in this area is also possible. In addition, the present detailed study of these characters in the Nepomorpha also provides an opportunity for comparison with other morphological hypotheses regarding the phylogenetic relationships, proposed by China [21], Popov [25], Rieger [26], and Mahner [10], even though their analyses were not algorithmic analyses.
The main goal of this paper is to clarify the significance of the characters of mouthparts structures (labial segments and sensilla, external and internal structures of maxillae and mandibles) in the relationships of the nepomorphan families based on cladistic analyses. Table 1 and used for the purpose of the study came from the collections of the Natural History Museum in Vienna, Zoological Museum of the State Moscow University, and the Paleontological Institute of the Russian Academy of Sciences in Moscow. The new characters presented in this paper have been described based on SEM images of the mouthpart structures. The SEM photographs were taken with a Hitachi scanning electron microscope.

Range of Characters.
A preliminary estimation of the characters of the maxillae, mandible structure, and labial sensilla with respect to their phylogenetic value based on the ground plan was compared with the basic model within the group (i.e., the basal taxa of the Nepidae and Belostomatidae) and with the more diverse forms of these structures in more evolutionarily advanced groups (i.e., Ochteridae, Gelastocoridae, Aphelocheiridae, Naucoridae, Pleidae, Helotrephidae, Notonectidae, Diaprepocoridae, Corixidae, and Micronectidae). All these data were taken from papers authored by Brożek [57,58]. The construction of the labial segments of most nepomorphans was presented by Brożek [59]; moreover, details of morphological characters of the labium in the Corixoidea (Nepomorpha) were also investigated by Brożek (2014 in press). Presently, the total of all characters proposed by , 2014 in press) provides an improved characteristics of these features which can be combined for the purpose of a more precise coding (Tables 2 and 3) with respect to the outgroup and analyzed from the phylogenetic perspective.

Outgroup Selection.
In the present study, the outgroup of the Gerromorpha was accepted with respect to the Nepomorpha according to the hypothesis proposed by 4 The Scientific World Journal  Wheeler et al. [12]. Different variants of phylogenetic relationships among infraorders of the Heteroptera have been established as mentioned in the Introduction. Additionally, the nepomorphan characters are polarized with respect to the Mesoveliidae, because this family is the most plesiomorphic one within the Gerromorpha [17]. The presently selected species of Mesovelia is congruent with the outgroup used in the study by Hebsgaard et al. [27]. Choosing the same outgroup as in the study by Hebsgaard et al. [27] provides a chance for the comparison of morphological characters of mouthparts which is methodologically correct, that is, by identical direction of polarization in the analysis.

Type of Coding.
Several characters of the outgroup used for the purpose of this paper originated from the description of the elements of the rostrum by Andersen [17] and Cobben [19] while others were based on the materials prepared presently (description and Figures 1(a)-3(d)) by Brożek. Characters and states selected as being of interest are marked as (Kn (state number K0-K47)). All of them are presented in Tables 2 and 3 Table 2. Character states were written into a standard character by taxon character state matrix (Table 3) with unknown characters coded as a question mark (?).
The studied species were coded as having individual characters (Table 3) to provide a more accurate reflection of the observed morphology, rather than trying to achieve uniformity of coding within the (sub)families. The characters used for analysis were based exclusively on adult structures.
The morphological characters presented in Tables 2 and  3 were shown according to the ground plan characters and the outgroup was shown as in the previous studies by Brożek [57][58][59]. In the present cladistic analysis all these characters were regarded as nonadditive and equally weighted in order to avoid regarding them in an a priori manner and to conduct estimation through algorithms adopted by the cladistic software. Additionally, such an estimation was necessary due to the presence of characters coded as unknown (?) in the Potamocoridae.

Programs Used for Cladistic Analysis.
Morphological data ( Table 3) were analysed using the parsimony programs NONA [60] and Winclada (BETA) ver. 0.9.9 [61] with equal weight characters and heuristic search with TBR transformation option. However, Goloboff [62,63] presented convincing justification for using implied weights in cladistic analysis and his method has since been widely used, with some authors preferring it to equal weights.
The nonhomoplasies and homoplasies on cladograms were searched using unambiguous and slow optimization in order to evaluate how the character data on the cladogram changed. Winclada, in particular, apply itself to investigating synapomorphies supporting nodes as it allowed for the mapping of all characters and states simultaneously. Additionally, characters were also analysed using the heuristic search option of PAUP * 4.0 [64]. All characters were used as nonordered, of equal weight with ACCTRAN transformation option, and character polarities were determined in the context of the phylogenetic analysis. The topology of trees and the arrangement of terminal taxa as well as a length, consistency, and retention index obtained in PAUP and NONA program were similar. NONA [60] and PIWE [65] were also used for the calculation of Bremer support values (decay index) for branches [66]. Runs were conducted using the following commands: Mult * 10; Max * ; and subsequently HOLD 1000; SUB 1; FIND * ; HOLD 2000; SUB 3; FIND * ; HOLD 4000; SUB 5; FIND * ; HOLD 5000; SUB 15; FIND * ; BSUPPORT. Bremer support values, shown in Figure 7, were calculated as measures of branch support up to 15 steps away from the most parsimonious solution. These values were also checked in the TNT program [67].
Bootstrap support implemented in Winclada [68] of 1000 resampling replicates was used to study the level of character 6 The Scientific World Journal         [19]. Presently, the original photographical documentation of Mesovelia (Figures 1(a)-1(h)) indicated that mandibles were evenly serrated apically and equipped with seven short spines (Figure 1(a)). On the basis of such appearance of the mandibular file they were included in the categories of evenly serrated (short spines) and medium length (K0 (0)). The observation of maxillary stylets (Figures  1(b) and 1(c)) showed that the apices were symmetrical (both apices straight and slightly narrow and flat (K1 (0)). The maxillary spines were stiff, long, and forming regular and dense external (brdex and brvex) and internal rows (brvin) along the edges of the maxillae, that is, exposed a rupturing device (Figure 1(a)) (K2 (0)).
In the cross-section (Figure 1(d)), the locked maxillae (Rmx and Lmx) appeared to be pentagonal in shape with the dorsal side distinctly tapered and wider than the ventral side. On the dorsal and ventral sides both had one pair of protuberant external lobe processes (depr, depl, vepr, and vepl) (K3 (0)). The mandibles (Rmd and Lmd) (K4 (0)) were placed on the lateral suboval walls of the maxillae between the dorsal and ventral lobes.
In Mesovelia there is substantial variation in the set of labial sensilla in comparison to the representatives of the Nepomorpha. Several short chaetica sensilla (CH3) (mechanosensilla) were found on the dorsal and ventral side of the III and IV segments (Figures 1(e) and 1(f)) (K5 (0)), whereas slightly longer chaetica sensilla (CH2) (K6 (0)) and long chaetica sensilla (CH1) (K7 (0)) were observed on the I and II segments (Figure 1(g)). Essentially, only one pair of the proprioceptive hairs (mechanosensilla) was situated on the ventral and dorsal sides of the II segment (K8 (0)) and one pair on the ventral side (K9 (0)). In Mesovelia, the characters mentioned in Table 2 from K10 to K22, K24, and K26 (0) were estimated as absent characters. Near the labial tip on the dorsal and ventral sides one pair of trichoid sensilla (TRS) was found (probably bimodal sensilla: mechanoreceptors and gustatory) (K23 (0) and K25 (0)). The characteristic type of sensilla was a plate-like, elongated sensillum (Wpples) (Figure 1(h)) present in the Mesoveliidae (K27 (0)) and 14 The Scientific World Journal  the Hebridae, while in the Nepomorpha it was absent [57].
In Mesovelia (Figure 1(h)) the peg-in-pit sensilla (poreless coeloconic sensilla) were not identified on the labial tip (K28 (0)); however, four peg sensilla (contact-chemoreceptive sensilla, mechano-and chemoreceptors) were observed centrally on the labial tip (K29 (0) and K30 (0)). These sensilla were inserted in the socket on the smooth surface of the tip (K31 (0)). In this species the sensilla on the labial segments were essentially less numerous and not very distinguished as well as unevenly arranged (K32 (0)). Generally, all these sensilla were classified on the basis of characters distinguished in many publications referring to this subject [58,[69][70][71][72]. The labium in Mesovelia showed a substantial similarity in structure to the representatives of the Nepomorpha (except for the Corixoidea). The labial apex on the ventral side was equipped with one oval plate (ap) (Figure 2(a)) (K33 (0)), which was similar to some species of the Nepomorpha. The Mesoveliidae (and Hebridae) appeared to be the only group with four large intercalary sclerites (is-dr, is-dl, is-vl, and is-vr (is-vr is invisible only in Figure 2(b)). These sclerites were situated on the distal edge of the third segment, and they surrounded the dorsal and ventral sides of the labium (K34 (0)). The edges of the dorsal surface of the labium were not in contact medially and the stylet groove was open (K35 (0)). The first labial segment was short ventrally and longer dorsally, generally ring shaped (K36 (0)) ( Figure 2(c)). The dorsal surface of this segment was covered by the labrum.
The second segment on the dorsal side was smooth (i.e., in that part the segment was not divided) (K37 (0)) ( Figure 2(c)) and the dorsal edges of the segment were in contact, so that the stylet groove was closed (K38 (0)). Also the lateral surface of the segment was uniform (without no incision) (K39 (0)). The stylet groove of the two basal segments of the labium was covered by an epipharyngeal projection. The second segment was usually the smallest of the four labial segments (Figure 2(c)) (K40 (0)). The third labial segment was by far the longest (K41 (0)). Typically, it was swollen proximally and tapered distally. The fourth segment was distally shorter than the preceding segment and tapering towards the pointed apex (K42 (0)). Ventrally, the distal edge of the first segment was hidden and the midventral condyle was estimated as lack of data (K43 (?)); on the third segment the midventral condyle was putatively present (K44 (0)). Evidently, the condyle was not observed on the proximal edge of the fourth segment (K45 (0)) ( Figure 2(d)). The second segment was connected with the third segment dorsally by a wider band of membrane (K46 (0)), dorsal articulation (cd)) ( Figure 2(c)). The labium was four segmented and tubular shaped (K47 (0)) ( Figure 2(e)).

Nepomorpha:
Potamocoridae. The set of characters required for the present analysis was selected on the basis of data from literature. The Potamocoridae are basically a group which has been scarcely investigated with respect to their morphology. Several studies focused only on the general morphology of the body [3,18,19,[73][74][75][76]; however anatomical details of particular elements of the body parts are not known. In the studies of the Nepomorpha conducted by Brożek [57][58][59] the Potamocoridae were not analyzed, as the material of those families was unavailable. Due to this, in the present study only several characters were analyzed which had been described previously by several authors. In the Potamocoridae most of characters (K) mentioned in the Table 2 were coded as unknown (lack of data). According to Cobben [19] maxillary stylets of Potamocoris sp. (Figure 13(B), pp [36][37] are structurally entirely different from the typical naucorid stylets. On the basis of the review of many maxillary stylets of nepomorphan taxa [57] it is possible to compare the maxillary stylets of Potamocoris sp. with other nepomorphan species. On the basis of their appearance, maxillary stylets of Potamocoris sp. (Figure 3(a)) were classified as stylets with asymmetrical apices (the right one (Rmx) was straight and narrow; the left one (Lmx) was wide and curved) (K2 (4) [19]. The right apex (ax) is wider than the left one; the right apex is slightly curved while the left one is straight and narrow. The inner system of spines on the maxillae edges from the reduced and hidden rupturing device (brdex and brvin). On the right maxilla seven short spines are present on the internal ventral side (brvin). On the dorsal side the spines are strongly reduced (brdex). On the left maxilla, internally, one tuft of short spines (brvin) can be observed. I: first segment, II: second segment, III: third segment, IV: fourth segment, ax: apex of maxillae, cd: dorsal condyle (the articulation between the second and third segments on the dorsal surface), gr: labial groove, tp: triangular plate of the second segment, Lr: labrum, Rmd: right maxilla, Lmx: left maxilla, Rmd: right mandible, Lmd: left mandible. the right maxilla (Rmx) (brdex and brvin) and one tuft with minor spines (brvin) on the left maxilla (Lmx)). When the maxillae were locked the spines were externally invisible. Such a system of spines was evaluated as the rupturing device almost reduced and hidden (K3 (9)), like the Corixoidea, Pleidae, and Helotrephidae [57]. The coded characters from K4 (?) to K34 (?) were treated mainly as a lack of data and referred to labial sensilla.
According to van Doesburg [18] the labium of Potamocoris nieseri was broad at the base, tapering to the tip of its third segment. The last segment was slightly shorter than the second one. On the basis of the drawing by van Doesburg [18] (Figure 2(a) (Potamocoris nieseri), pp. 22) it was possible to estimate that there were no intercalary sclerites (K35 (8)) ( Figure 3(b)). A similar conclusion was drawn by Cobben [19]. Even though the drawing was based on the picture from the light microscope, certain structures of the labium could be recognized and compared to the SEM images of the labium of other nepomorphans. According to my experience, the drawing of Potamocoris nieseri showed the fourth segment of the labium and its appearance was similar to the labium of naucorids and pleids.
I took the liberty of describing the drawing made by van Doesburg [18] according to categories, which were used by Brożek [59], and to introduce these characters to the present analysis. The first segment (I) and partly the second (II) segment were covered by the triangular labrum (Lr). There was no certainty as to the type of the stylet groove of the first segment, so that the K36 (?) was estimated as a lack of data. The first segment was rather narrow (K37 (1)); the lateral sides were visible and reaching to the base of the labrum (Figures  3(b) and 3(c)). In Figure 3(b), the laterally and dorsally visible elements belonging to the second segment (II) corresponded to the elements (tp and cp) putatively marked in Figure 3(c). The dorsal surface of the second segment in this species could be divided into a triangular, flat plate (tp) and a second plate (cp) (K38 (1)). They were placed symmetrically on the left and right side of the stylet groove. In the second segment the stylet groove was usually open up to the half-length of the segment (K39 (2)) ( Figure 3(c)). For the Potamocoris it was assumed that the lateral surface of the second segment was smooth (K40 (0)) as in most nepomorphans. The first and second segments were short (K41 (0). The third segment (III) (mentioned as the second one by van Doesburg [18]) was long (K42 (1)) in comparison to the first and second segment, and the fourth one was shorter than the third (K43 (0)). The midventral condyle (K44 (?), K45 (?), and K46 (?)), either present or absent in various nepomorphans, was estimated as a lack of data for Potamocoris nieseri. Dorsally, the third and second segment had two points of articulation (band shaped (K47(0)). On the basis of Figure 3(b) it could be suggested that the labium was four segmented and tubular shaped (K48 (0).

Characters (K0-4): The Shapes of Maxillae and Mandibular Stylets of the Nepomorpha according to
Brożek [57] ( Table 2). General stylet structures were used in a prior cladistic analysis of relationships within the Heteroptera by Cobben [19]. Presently, the condition of stylets found in the Nepomorpha indicated a substantial variation in stylet structure within the group. Ten morphologically distinct types of files (K0) were identified on the mandibular tip in individual species, as well as eight distinct types of maxillary endings (K1) in individual species and ten distinct types of rupturing devices (K2) of the maxillae. The features of the internal maxillary (K3) and mandibular (K4) structures shared a common connection model, differing only by virtue of specific appendages in different subfamilies.

Characters (K5-32): Labial Sensilla Types and Distribution Patterns of Sensilla in the Nepomorpha according to
Brożek [58] (Table 2). Twenty-one morphologically distinct types of the mechanosensilla as well as two types of the trichoid sensilla (contact-chemoreceptive sensilla) were identified on all labial segments in representatives of the subfamilies. The chaetica sensilla (CH3, CH2, and CH1) were present in various layouts on the segments (K5, K6, and K7). The proprioceptive sensilla were positioned on the dorsal side (K8) (either one pair or more pairs) and on the ventral side (K9) (one pair) on the second segment of the labium. Several variously shaped mechanosensilla were specific for individual species (K10-K22). Variously shaped trichoid sensilla are placed on the IV segment (K23, K24, and K25) and on the III segment (K26). Near to the labial tip, subapically, the elongated plate sensillum was present in representatives of the outgroup (K27); however, it was absent in the Nepomorpha. On the labial tip of the nepomorphans, three morphologically distinct types of chemosensilla were identified: one type of the peg-in-pit sensilla (K28) and two types of papillae sensilla (K29), as well as various types of their distribution. In addition, these sensilla were present in various numbers, from a few to a dozen (K30). The sensilla were inserted in the labial tip, either smooth or folded (K31). The mechanosensilla were present and placed in groups or rows distributed along the labium near the labial groove on the dorsal side; the sensilla were also unevenly scattered over the ventral surface of that segments (K32). Brożek [59] (Table 2). Within the thirteen families, six morphologically distinct forms of the apical plate (K33) of the labium and several intercalary sclerites (K34) were identified. Although in most investigated taxa of the nepomorphans subsequent segments of the labium (I, II, III, and IV) were shaped similarly, individual characters in some (sub)families differed (K35-42). The presence of the midventral condyle on the distal edge of the first segment (K43) and the third segment (K44) was observed, but not in all species. A new position of the midventral condyle on the proximal edge of the fourth labial segment (K45) was distinguished in several groups. Additionally, three types of articulation (K46) on the dorsal side between the third and second segments were interpreted as the new characters in relation to previous studies of this area.

Characters (K33-47): Shape of the Labial Segments of the Nepomorpha according to
The labium showed a substantial variation in the structure and segmental development between the Corixoidea and the remaining nepomorphans. The Corixoidea appeared to be the only group in which the first and second segment were completely lost on the dorsal side. Generally, the labium is triangular-shaped and short (K47)); however there had been evidence that the third and fourth segments were conspicuously present ( [77]; Brożek, 2014 in press).

Morphological Characters Mapped on the Parsimonious
Tree. Character analysis (complete data matrix presented in Tables 2 and 3).
The heuristic search strategy yielded 100 parsimonious trees, 199-98 steps long and with the consistency index = 72 and retention index = 92. Two of the shortest trees (198 steps long) (Figures 4 and 5) and consensus tree ( Figure 6, 221 steps long) with the complete mapping of all morphological characters as nonhomoplasious = syn(apomorphies) and homoplasious represent the hypothesis with reference to the relationship within the Nepomorpha given below. The most parsimonious tree with branch support (bs = 1 for 15 in individual branches) values [66] is shown in Figure 7. The bootstrap analysis of morphological characters is also shown in Figure 7.
The monophyly of the lineage that includes the Notonectidae and the Pleidae + Helotrephidae is supported by one synapomorphy . Two synapomorphies, that is, (5-7) and , have been found for the Notonectidae (Notonectoidea). The monophyly of the superfamily Pleoidea (Pleidae and Helotrephidae) is supported by two unambiguous synapomorphies: (37-4) and (46-3). Each family also brings an individual synapomorphy: Pleidae) and Helotrephidae).
In this tree topology, the most (super/sub)families are found to be monophyletic; on the basis of the present data only the family of Potamocoridae (Potamocoris nieseri) is problematic, as no autapomorphy has been found.
The second equally parsimonious tree ( Figure 5) hypothesizes the monophyly of the Nepomorpha and also finds sister relationships among most taxa in a similar way as in Figure 4. A major difference with respect to the previous tree (Figure 4), obtained also under equal weights, is the position of the Corixoidea. They are placed as a basal taxon instead of the Nepoidea ( Figure 5). The most synapomorphies and autapomorphies marked in black box are the same as in the reconstruction discussed above (Figure 4).
The ambiguity in relationships among nepomorphan taxa are illustrated in the consensus tree ( Figure 6). The unresolved relationships among the some species are pointed within the Corixoidea and Belostomatidae. Also the unresolved relationships are visible among the Ochteridae, Nerthrinae, and Gelastocorinae and among subfamilies of the Naucoridae (Laccocorinae, Limnocorinae, Cryphocricinae, and Naucorinae) and Potamocoridae. The polytomies are also visible among the species of Notonectidae and Helotrephidae.
Generally, the bootstrap analysis is seldom used for morphological analyses; however, its use for the purpose of the present study seems to be necessary. In 100 parsimonious trees with the same parameters (L, CI, and RI), the nodes change within the range of the analyses, mainly regarding the positions of the Corixoidea and the Nepoidea. The bootstrap analysis (Figure 7) shows that the character of dataset is robust with regard to the hypothesis of the monophyly of the Nepomorpha (i.e., the clade is found in 100% of the trees). The basal group of the Nepoidea is indicated in 87% of the trees, while the Nepidae is found in 98% and Belostomatidae in 96%, respectively. Those high rates of support are also maintained for the subfamilies Ranatrinae: 84%, Nepinae: 96%, and Belostomatinae: 69%. The Corixoidea, placed as a sister groups with respect to the remaining nepomorphans, is evaluated in 61% of these trees. The clade Corixoidea has received 100% support in these trees, although individual families are found in 50% of the trees.
Other taxa, except for the Nepoidea and the Corixoidea, are visible in 64% of these trees. Nonetheless, there is a weak bootstrap support (slightly above 50% of the bootstrap trees) in this dataset for a sister group relationship between the Ochteridae and the Gelastocoridae. The clade Nerthrinae + Gelastocorinae is hypothesized to be monophyletic in 44% of the trees, a fairly low bootstrap value that reflects the unstable position of the Nerthrinae in the equally parsimonious trees.        the Aphelocheiridae is found in 73% of the trees, and a relatively high bootstrap value of 77% also supports this relationship for the Potamocoridae. However, lower bootstrap support values (53-57-55% in three nodes) have been received for the Naucoridae. Nevertheless, the clade Cheirochelinae is found in 63% of the trees, the Laccocorinae in 50%, and three other of subfamilies in 67% of these trees. A sister group relationship between the Helotrephidae + Pleidae and the Notonectidae is found in 64% of the trees. A sister group relationship between the Helotrephidae + Pleidae is found in 94% of the trees and reflects the unambiguous phylogenetic position of this clade in the equally parsimonious trees. Moreover, the relationships within the Notonectidae evaluated in 83%-50% of the bootstrap trees are robustly supported in this dataset.
There is no bootstrap support (attainment of the 50% bootstrap level) for the hypothesized relationships between the representatives of Macrocoris, Limnocoris, Ambrysus, Cryphocricos, Naucoris, Neomacrocoris, and Namtokocoris as well as the helotrephid species in this dataset, even though the relationships among them are consistent in all the shortest trees. This lack of bootstrap support, in contrast to the consistent placement of these taxa in the shortest trees, reflects the fact that relatively few, but highly consistent, characters support the nodes.

The Main Phylogenetic Hypothesis of Relationships within the Nepomorpha (to Be Presented in a Planned Discussion).
Previously, the phylogeny of the Nepomorpha was discussed on the basis of various morphological criteria used by China [21], Popov, [25], Rieger [26], and Mahner [10] and the hypotheses proposed by them brought about several different solutions.
Recent hypotheses regarding the relationships among taxa within the Nepomorpha based on rigorous cladistic assumptions and on molecular and morphological studies have been proposed by Hebsgaard et al. [27], and hypotheses based on molecular studies have been proposed by Hua et al. [28] and Li et al. [29]. Phylogenetic analysis of Hebsgaard et al. [27] was generally congruent with the traditional classification of Mahner [10]; however, a new superfamily of the Aphelocheiroidea (Aphelocheiridae + Potamocoridae) was evaluated, and moreover the Naucoroidea were restricted to only one family (Naucoridae) and the Ochteroidea received a new position (Ochteridae + Gelastocoridae). According to Hebsgaard et al. [27], the system of classification of the Nepomorpha included seven monophyletic superfamilies, namely, (Nepoidea, Corixoidea, Aphelocheiroidea, Naucoroidea, Ochteroidea, Notonectoidea, and Pleoidea (Pleidae + Helotrephidae)). A revised (or suggested) classification of the Nepomorpha by Hebsgaard et al. [27] based on a molecular dataset (genome) found support for just five superfamilies in the new distribution of these taxa (Corixoidea + ((Naucoroidea + Notonectoidea) + (Ochteroidea + Nepoidea). Two superfamilies from Hebsgaard et al. [27] system were lost; the Pleoidea (Pleidae + Helotrephidae) was placed in the new infraorder of the Plemorpha, while the Aphelocheiroidea sensu Hebsgaard's et al. were included into the Naucoroidea sensu Hua et al. [29].
The above presented phylogenetic analyses resulted in totally different hypotheses regarding the Nepomorpha; therefore future studies in this field seem necessary.

How the New Values for Phylogeny and Classification System of the Nepomorpha Represent the Dataset Concerning
Mouthparts and Labial Sensilla Structures? Presently, the monophyly of the Nepomorpha is supported by an unambiguous synapomorphy (lack of the elongated plate sensillum on the labium (27-1)). This hypothesis is concordant with that of Popov [25], Mahner [10], Hebsgaard et al. [27], and Li et al. [29] but contradictory to the views of Hua et al. [28], who treated the Nepomorpha as a monophyletic group, excluding the Pleoidea.
Problems with the relationships of families (or superfamilies) represented by the 62 species in the present analysis seem far more extensive and complicated. The two most parsimonious trees shown in Figures 4 and 5 provide distinctly different solutions.

4.3.
Relationships of (Super)families. The present system of relationships among taxa demonstrated in the first tree ( Figure 4) alludes to most of the previous hypotheses [10,[25][26][27] with respect to the basal position of the superfamily Nepoidea (Belostomatidae and Nepidae). Presently, the superfamily is supported by three unambiguously optimized synapomorphic characters and each family is also strongly evaluated through optimized characters (i.e., the Nepidae by three characters and the Belostomatidae by four characters). In addition, the subfamilies (Nepinae and Ranatrinae) are regarded as monophyletic groups in contrast to the estimations provided by Hebsgaard et al. [27] and Mahner [10], which indicated the paraphyletic characters of the Nepinae. The next position in the presented system of relationships treats the Corixoidea as a sister group of the remaining taxa of the nepomorphans. For them, three unambiguously optimized synapomorphic characters have been found. Three families are listed within this taxon and each of them is characterized by one autapomorphy. Such an arrangement of the Corixoidea (with one family Corixidae) finds support in the studies of Popov [25], Mahner [10], and Hebsgaard et al. [27]. Then, the superfamily Ochteroidea (Ochteridae and Gelastocoridae) is positioned bellow the Corixoidea, however, in a different position from the one it had in the cladogram developed by Hebsgaard et al. [27]. In several other studies, the Ochteroidea have been placed in various positions across the system of classification. The present result confirms the hypothesis proposed by Popov [25] and Mahner [10]. The Aphelocheiridae and Potamocoridae are positioned in a similar way to the final tree ( Figure 23; Hebsgaard et al. [27]); however, the superfamily Aphelocheiroidea (sensu Hebsgaard et al. [27]) presently is not formed. The Aphelocheiridae are a sister group to the Potamocoridae as a separate branch. At next step, the Potamocoridae are located separately as a sister group to the Naucoridae. Previous studies, mainly by Rieger [26], indicated a close relationship between the Potamocoridae and the Naucoridae; however, Popov [25] and Mahner [10] found a relationship between the Potamocoridae and a clade Naucoridae + Aphelocheiridae. Presently, the Naucoridae are regarded as a monophyletic family (one synapomorphy has been found) also confirmed by the studies of Hebsgaard et al. [27]. As for the remaining groups in the tree (Figure 4), namely, the Notonectidae, Helotrephidae, and Pleidae, their relationships are reminiscent of the relationships indicated by Hebsgaard et al. [27]. Actually, the clade Helotrephidae + Pleidae (Pleoidea) is also a sister group to the Notonectidae (Notonectoidea).
The classification system and the relationships among super(families) of the Nepomorpha with the Ochteroidea as the basal lineage and the configuration of the clade Pleoidea + Naucoroidea as well as the clade Nepoidea + Corixoidea presented by  are totally different from the current data (Figures 4, 5, and 6) and the previous study by Hua et al. [28], Hebsgaard et al. [27], Manher [10], and Rieger [26]. It ought to be emphasized that the Ochteroidea as a basal group has been indicated only by China [21]. Nevertheless, Popov [25] suggested that the Nepomorpha could derive from ochterid-like ancestors but that they rather derived from saldid-like stock. However, on the basis of detailed studies of the comparative morphology of the families of the Nepomorpha Parsons [22,48,77] suggested that the Ochteridae and Gelastocoridae were more specialized.

Relationship Nepomorpha in the Groups of Taxa.
The consideration of the relationships presented in the cladogram (Figures 4 and 5) in wider range of comparisons among the families yields interesting results. Essentially, in the three cladograms there is a visible group of families (Ochteridae, Gelastocoridae, Aphelocheiridae, Potamocoridae, Naucoridae, Notonectidae, Helotrephidae, and Pleidae) supported by one synapomorphy (the pit chemoreceptive sensillum is present in the mentioned taxa) which correspond to the group Tripartita previously indicated by Manher [10]. This group was found in the morphological analysis an in the simultaneous morphological and molecular analysis conducted by Hebsgaard et al. [27] but it was not supported by the same molecular data. The present study also strongly supports a group composed of the families Aphelocheiridae, Potamocoridae, Naucoridae, Notonectidae, and Pleidae + Helotrephidae based on three synapomorphies. These families correspond to the group Cibariopectinata distinguished by Mahner [10] as well as to the Cibariopectinata composed of a polytomy of the clades Potamocoridae, Aphelocheiridae + Naucoridae, and Notonectidae + Pleoidea (Pleidae + Helotrephidae) presented only in the morphological analysis by Hebsgaard et al. [27]. Presently obtained data with respect to the Notonectidae as a sister group to the Pleidae + Helotrephidae (Pleoidea) are congruent with previously obtained data reported by Hebsgaard et al. [27], Manher [10], Rieger [26], Popov [25], and China [21] except for the concepts proposed by Hua et al. [28] and Li et al. [29].

Taxonomic and Phylogenetic Placement of Corixoidea.
With regard to the tree in Figure 5, there is substantial congruence among the results of the present analyses and some hypotheses of relationships proposed by Hua et al. [29]. These include the essential points in the phylogenetical estimation of the Nepomorpha. The basal position of the Corixoidea ( Figure 5) diagnosed according to the present analyses resembles the results of analyses obtained by Hua et al. [28]. Nonetheless, two points of obvious ambiguity distinguish these analyses: according to Hua et al. [28] the Notonectidae are placed as the sister group to the Aphelocheiridae + Naucoridae, whereas in the present tree ( Figure 5) the distribution of the remaining taxa corresponds to the tree in Figure 4 (Aphelocheiridae, Potamocoridae, Naucoridae, Notonectidae + (Pleidae + Helotrephidae)). Moreover, as has already been mentioned above, according to Hua et al. [28] the Pleoidea do not belong to the Nepomorpha.
The present placement of the Corixoidea is different from previous several hypotheses. As for Hua et al. [28], they stated that in their results the Corixoidea was always the most basal taxon within the Nepomorpha, whereas presently the Corixoidea in the basal position is estimated only in 40% of the trees.
The aberrant morphology of the Corixidae (Corixoidea) has puzzled phylogeneticists and hence several different hypotheses have been developed about the place of the Corixidae in the systematics. Börner [78] proposed a separate division of the Sandaliorrhyncha family. However, it is now well established that the Corixidae belongs to the Nepomorpha. Both Parsons [22] and Popov [25] indicated a divergence of the Corixidae after the Nepoidea in their phylogenetic dendrograms and they stated that that group was very advanced developmentally and represented many apomorphic states. The evidence pointing to derived characters of corixoids is significant and several examples can be cited. The triangular-shaped labium is an evolutionary novelty in this group; however, it derived from the tubular four-segmented labium of other ancestral nepomorphans. In turn, the mandibles of corixid bugs shared a common pattern with other water bugs [79], but Brożek [57] pointed out similarities in the mandibles of the corixids only with respect to the mandibles of the Gelastocoridae. Moreover, the structure of the maxillary stylet in corixids was their specific characteristics, not encountered elsewhere. Characteristics of the internal structure of the mouthparts show a similar type across the Nepomorpha, indicating that the Corixoidea belongs to this infraorder. In most nepomorphan taxa the sensilla are placed along the long axis of the labium, The Scientific World Journal 27 while in the Corixoidea (except for the Cymatiainae) these sensilla are placed in the transverse bands on the labium. Different types of contrast in the distribution of sensilla on the surface of the labium can be noticed between the Corixoidea and the remaining nepomorphan families. There exists a transverse pattern of distribution of the sensilla and an autapomorphy in the case of the Corixoidea (except for the Cymatiainae). Many other modified body structures of the Corixoidea have also reached a new level of adaptation among nepomorphan taxa, and therefore they represent an advanced systematic position contrary to the suggestion of Hua et al. [28].

Unresolved Ancestral Nodes in the Consensus
Tree. The consensus tree ( Figure 6) formed on the basis of 100 parsimonious trees shows poorly resolved ancestral nodes leading to the divergence into respective taxa. This indicates that there is a substantial degree of disagreement among the trees regarding individual parameters, although the characters/characters states have been weighed. This especially refers to two groups: the Nepoidea and the Corixoidea; the positions of these taxa are a major problem. Moreover, an unsatisfactory degree of relationships is also evident in the group of naucorids. Unresolved relationships are shown between many naucorid species and the Potamocoridae.
In addition, the superfamily Ochteroidea has been spread out over three polytomous taxa: the Ochteridae, Gelastocorinae, and Nerthrinae.
It can be expected that in future cladistic studies taking into account a wider range of morphological characters will stabilize the positions of most clades that have been recognized so far.
The first step towards achieving such goal can be combining the characters from the present matrix with the morphological matrix used by Hebsgaard et al. [27]. However, it would require further work on a number of significant features so that they would correspond with the list of species that have already been studied.