Phylogeny in Echinocereus (Cactaceae) based on combined morphological and molecular evidence: taxonomic implications

Echinocereus is a morphologically diverse genus that includes 64 species grouped into eight taxonomic sections based on morphological traits. In previous molecular phylogenetic analyses, the relationships amongst Echinocereus species were not entirely revealed and useful characters to recognize clades were not provided. The inclusion of several sources of evidence in a phylogenetic analysis is likely to produce more supported hypotheses. Therefore, we performed a combined phylogenetic analysis with a set of 44 morphological characters and six chloroplast DNA sequences. Topologies from parsimony and Bayesian analyses were mostly congruent. However, the relationships of E. poselgeri were not consistent between analyses. A second Bayesian analysis using a long-branch extraction test resulted in a topology with the morphological position of E. poselgeri congruent with that in parsimony analysis. Parsimony and Bayesian analyses corroborated the monophyly of Echinocereus, which included eight monophyletic groups. The combined phylogeny integrated into different clades those taxa that were not determined in previous analyses and changed the relationships of some recognized clades. The clades did not recover the recent infrageneric classification. In the present study, a new sectional classification for Echinocereus is proposed based on the eight recovered clades, which is supported by a combination of morphological and molecular characters. An identification key for sections in the genus is included.     

Phylogenetic analysis of the Monocotylidae (Monogenea) inferred from 28S rDNA sequences

The current classification of the Monocotylidae (Monogenea) is based on a phylogeny generated from morphological characters. The present study tests the morphological phylogenetic hypothesis using molecular methods. Sequences from domains C2 and D1 and the partial domains C1 and D2 from the 28S rDNA gene for 26 species of monocotylids from six of the seven subfamilies were used. Trees were generated using maximum parsimony, neighbour joining and maximum likelihood algorithms. The maximum parsimony tree, with branches showing less than 70% bootstrap support collapsed, had a topology identical to that obtained using the maximum likelihood analysis. The neighbour joining tree, with branches showing less than 70% support collapsed, differed only in its placement of Heterocotyle capricornensis as the sister group to the Decacotylinae clade. The molecular tree largely supports the subfamilies established using morphological characters. Differences are primarily how the subfamilies are related to each other. The monophyly of the Calicotylinae and Merizocotylinae and their sister group relationship is supported by high bootstrap values in all three methods, but relationships within the Merizocotylinae are unclear. Merizocotyle is paraphyletic and our data suggest that Mycteronastes and Thaumatocotyle, which were synonymized with Merizocotyle after the morphological cladistic analysis, should perhaps be resurrected as valid genera. The monophyly of the Monocotylinae and Decacotylinae is also supported by high bootstrap values. The Decacotylinae, which was considered previously to be the sister group to the Calicotylinae plus Merizocotylinae, is grouped in an unresolved polychotomy with the Monocotylinae and members of the Heterocotylinae. According to our molecular data, the Heterocotylinae is paraphyletic. Molecular data support a sister group relationship between Troglocephalus rhinobatidis and Neoheterocotyle rhinobatidis to the exclusion of the other species of Neoheterocotyle and recognition of Troglocephalus renders Neoheterocotyle paraphyletic. We propose Troglocephalus incertae sedis. An updated classification and full species list of the Monocotylidae is provided.     

Phylogenetic relationships among families of the Scaphopoda (Mollusca)

Phylogenetic relationships among families in the molluscan class Scaphopoda were analysed using morphological characters and cladistic parsimony methods. A maximum parsimony analysis of 34 discrete characters, treated as unordered and equally weighted, from nine ingroup terminal taxa produced a single most parsimonious tree; supplementary analyses of tree length frequency distribution and Bremer support indices indicate a strong phylogenetic signal from the data and moderate to minimally supported clades. The traditional major division of the class, the orders Dentaliida and Gadilida, is supported as both taxa are confirmed as monophyletic clades. Within the Dentaliida, two clades are recognized, the first comprised of the families Dentaliidae and Fustiariidae, the second of the Rhabdidae and Calliodentaliidae; together, these groups comprise a third clade, which has the Gadilinidae as sister. Within the Gadilida, a nested series of relationships is found among [Entalinidae, [Pulsellidae, [Wemersoniellidae, Gadilidae]]]. These results lend cladistic support to earlier hypotheses of shared common ancestry for some families, but are at variance with other previous hypotheses of evolution in the Scaphopoda. Furthermore, analysis of constituent Gadilinidae representatives provide evidence for paraphyly of this family. The relationships supported here provide a working hypothesis that the development of new characters and greater breadth of taxonomic sampling can test, with a suggested primary goal of establishing monophyly at the family level.     

Phylogenetic analysis of the Monocotylidae (Monogenea) inferred from 28S rDNA sequences

The current classification of the Monocotylidae (Monogenea) is based on a phylogeny generated from morphological characters. The present study tests the morphological phylogenetic hypothesis using molecular methods. Sequences from domains C2 and D1 and the partial domains C1 and D2 from the 28S rDNA gene for 26 species of monocotylids from six of the seven subfamilies were used. Trees were generated using maximum parsimony, neighbour joining and maximum likelihood algorithms. The maximum parsimony tree, with branches showing less than 70% bootstrap support collapsed, had a topology identical to that obtained using the maximum likelihood analysis. The neighbour joining tree, with branches showing less than 70% support collapsed, differed only in its placement of Heterocotyle capricornensis as the sister group to the Decacotylinae clade. The molecular tree largely supports the subfamilies established using morphological characters. Differences are primarily how the subfamilies are related to each other. The monophyly of the Calicotylinae and Merizocotylinae and their sister group relationship is supported by high bootstrap values in all three methods, but relationships within the Merizocotylinae are unclear. Merizocotyle is paraphyletic and our data suggest that Mycteronastes and Thaumatocotyle, which were synonymized with Merizocotyle after the morphological cladistic analysis, should perhaps be resurrected as valid genera. The monophyly of the Monocotylinae and Decacotylinae is also supported by high bootstrap values. The Decacotylinae, which was considered previously to be the sister group to the Calicotylinae plus Merizocotylinae, is grouped in an unresolved polychotomy with the Monocotylinae and members of the Heterocotylinae. According to our molecular data, the Heterocotylinae is paraphyletic. Molecular data support a sister group relationship between Troglocephalus rhinobatidis and Neoheterocotyle rhinobatidis to the exclusion of the other species of Neoheterocotyle and recognition of Troglocephalus renders Neoheterocotyle paraphyletic. We propose Troglocephalus incertae sedis. An updated classification and full species list of the Monocotylidae is provided.     

A priori assumptions about characters as a cause of incongruence between molecular and morphological hypotheses of primate interrelationships

When molecules and morphology produce incongruent hypotheses of primate interrelationships, the data are typically viewed as incompatible, and molecular hypotheses are often considered to be better indicators of phylogenetic history. However, it has been demonstrated that the choice of which taxa to include in cladistic analysis as well as assumptions about character weighting, character state transformation order, and outgroup choice all influence hypotheses of relationships and may positively influence tree topology, so that relationships between extant taxa are consistent with those found using molecular data. Thus, the source of incongruence between morphological and molecular trees may lie not in the morphological data themselves but in assumptions surrounding the ways characters evolve and their impact on cladistic analysis. In this study, we investigate the role that assumptions about character polarity and transformation order play in creating incongruence between primate phylogenies based on morphological data and those supported by multiple lines of molecular data. By releasing constraints imposed on published morphological analyses of primates from disparate clades and subjecting those data to parsimony analysis, we test the hypothesis that incongruence between morphology and molecules results from inherent flaws in morphological data. To quantify the difference between incongruent trees, we introduce a new method called branch slide distance (BSD). BSD mitigates many of the limitations attributed to other tree comparison methods, thus allowing for a more accurate measure of topological similarity. We find that releasing a priori constraints on character behavior often produces trees that are consistent with molecular trees. Case studies are presented that illustrate how congruence between molecules and unconstrained morphological data may provide insight into issues of polarity, transformation order, homology, and homoplasy.     

COMBINED AND SEPARATE ANALYSES OF MORPHOLOGICAL AND MOLECULAR DATA IN THE PLANT FAMILY RUBIACEAE

Abstract— The Rubiaceae are one of the largest of the families of angiosperms, with over 10000 species. The tribal and subfamilial classification is provisional due to the lack of phylogenetic hypotheses. The present study of the Rubiaceae is based on 33 genera and three data sets, one morphological and two molecular from chloroplast DNA, restriction sites andrbcL sequences. There is much congruence between the morphological and the molecular data sets, but also conflict. For parsimony reasons, the best phylogenetic hypothesis is a tree based on an analysis of the combined data sets. The so-called “total evidence” criterion for the combined analysis is simply a reiteration of the principle of parsimony. In this particular study, the classification would be almost the same even if based on the separate analyses instead of the combined. Despite the inapplicability of consensus trees or trees from separate analyses for phylogenetic hypotheses and classification, separate analyses may provide important information. It is the best way to reveal conflicts between different data sets. Knowledge of the conflicts can promote further detailed investigation in order to improve understanding of characters and phylogenetic hypotheses. In this study, the tribe Vanguerieae provides such an example; morphological data support a position in the subfamily Cinchonoideae, but DNA and a tree based on the combined data support a position in subfamily Ixoroideae. The tribe's position in the morphological tree is probably due to missing information concerning the correct pollen presentation system. Bootstrap fractions and K. Bremer's branch support values are used to evaluate the stability of particular nodes in the trees. Interestingly these values are not always correlated, e.g. in the morphological tree, the node with the highest branch support value has very low bootstrap fraction. The reasons for these differences are unclear, but large differences are presumably more likely to occur on short branches.     

分支系统学当前的理论和方法概述及华东地区山胡椒属十二个种的分支系统学研究

本文概述了当前分支系统学研究中涉及的主要理论和方法,包括性状的选取、性状状态和极性的确定、数据矩阵的分析计算、结果分支图的处理、分支图可靠性的评价及分支图的应用。本文同时以华东地区樟科山胡椒属Linderal2个种的分支系统学研究为例,讨论了用形态性状进行分支系统学研究中可能遇到的问题,也揭示了一些分支系统学与传统的系统学在应用性状推导进化关系上的不同点。对这12个种的分支系统学研究得出了一些不同于传统系统学方法所推测的山胡椒属内的系统发育关系,如分支系统学研究显示山胡椒组和球果组很近缘。在严格一致性分支图上,杯托组的黑壳楠和江浙山胡椒分别位于最原始和最进化的分支,表明这个组是复系类群。分支图也显示山胡椒组可能是复系类群。     

ANOTHER MONOPHYLY INDEX: REVISITING THE JACKKNIFE

Abstract — Randomization routines have quickly gained wide usage in phylogenetic systematies. Introduced a decade ago, the jackknife has rarely been applied in cladistic methodology. This data resampling technique was re-investigated here as a means to discover the effect that taxon removal may have on the stability of the results obtained from parsimony analyses. This study shows that the removal of even a single taxon in an analysis can cause a solution of relatively few multiple equally parsimonious trees in an inclusive matrix to result in hundreds of equally parsimonious trees with the single removal of a taxon. On the other hand, removal of other taxa can stabilize the results to fewer trees. An index of clade stability, the Jackknife Monophyly Index (JMI) is developed which, like the bootstrap, applies a value to each clade according to its frequency of occurrence in jackknife pseudoreplicates. Unlike the bootstrap and earlier application of the jackknife, alternative suboptimal hypotheses are not forwarded by this method. Only those clades in the most parsimonious tree(s) are given JMI values. The behaviour of this index is investigated both in relation to a hypothetical and a real data set, as well as how it performs in comparison to the bootstrap. The JMI is found to not be influenced by uninformative characters or relative synapomorphy number, unlike the bootstrap.     

A Review of Current Theories and Methods in Cladistics and a Cladistic Study of Twelve Lindera Species in Eastern China

Cladistics has become a widely used method for phylogenetic reconstruction．Because of rapid improvement Of cladistic theories and methodologies,and application of new data,especially,molecular data,it is becoming realistic to reconstruct phylogenies of organisms,and to establish natural classifications based on these phylogenies．This paper reviews some current cladistic theories and methods in a practical way,such as choosing characters,defining character states,polarizing characters,analyzing data matrices, calculating consensus cladograms,choosing among multiple equally most parsimonious cladograms,estimating reliability of cladograms,and applying cladograms to classification, character evolution,and biogeography． Based on 36 morphological characters．a parsimony analysis of 12 species representing six sections in subgenus Lindera and an outgroup species from subgenus lteodaphne of the genus Lindera(Lauraceae)was conducted．The results suggest a close relationship between section Lindera and section Sphaerocarpae,which is different from the previous phylogenetic hypothesis within the genus．In the strict consensus cladogram,two species,L．megaphylla and L.chienii,from section Cupuliformes are in the most primitive and the most advanced clades respectively,indicating that the section is polyphyletic．The cladogram also suggests that section Lindera be a polyphyletic group．     

Phylogenetic relationships of Steinernema Travassos, 1927 (Nematoda: Cephalobina: Steinernematidae) based on nuclear, mitochondrial and morphological data

Entomopathogenic nematodes of the genus Steinernema are lethal parasites of insects that are used as biological control agents of several lepidopteran, dipteran and coleopteran pests. Phylogenetic relationships among 25 Steinernema species were estimated using nucleotide sequences from three genes and 22 morphological characters. Parsimony analysis of 28S (LSU) sequences yielded a well-resolved phylogenetic hypothesis with reliable bootstrap support for 13 clades. Parsimony analysis of mitochondrial DNA sequences (12S rDNA and cox 1 genes) yielded phylogenetic trees with a lower consistency index than for LSU sequences, and with fewer reliably supported clades. Combined phylogenetic analysis of the 3-gene dataset by parsimony and Bayesian methods yielded well-resolved and highly similar trees. Bayesian posterior probabilities were high for most clades; bootstrap (parsimony) support was reliable for approximately half of the internal nodes. Parsimony analysis of the morphological dataset yielded a poorly resolved tree, whereas total evidence analysis (molecular plus morphological data) yielded a phylogenetic hypothesis consistent with, but less resolved than trees inferred from combined molecular data. Parsimony mapping of morphological characters on the 3-gene trees showed that most structural features of steinernematids are highly homoplastic. The distribution of nematode foraging strategies on these trees predicts that S. hermaphroditum, S. diaprepesi and S. longicaudum (US isolate) have cruise forager behaviours.     

A methodological investigation of hominoid craniodental morphology and phylogenetics

The evolutionary relationships of extant great apes and humans have been largely resolved by molecular studies, yet morphology-based phylogenetic analyses continue to provide conflicting results. In order to further investigate this discrepancy we present bootstrap clade support of morphological data based on two quantitative datasets, one dataset consisting of linear measurements of the whole skull from 5 hominoid genera and the second dataset consisting of 3D landmark data from the temporal bone of 5 hominoid genera, including 11 sub-species. Using similar protocols for both datasets, we were able to 1) compare distance-based phylogenetic methods to cladistic parsimony of quantitative data converted into discrete character states, 2) vary outgroup choice to observe its effect on phylogenetic inference, and 3) analyse male and female data separately to observe the effect of sexual dimorphism on phylogenies. Phylogenetic analysis was sensitive to methodological decisions, particularly outgroup selection, where designation of Pongo as an outgroup and removal of Hylobates resulted in greater congruence with the proposed molecular phylogeny. The performance of distance-based methods also justifies their use in phylogenetic analysis of morphological data. It is clear from our analyses that hominoid phylogenetics ought not to be used as an example of conflict between the morphological and molecular, but as an example of how outgroup and methodological choices can affect the outcome of phylogenetic analysis.     

Systematics of Ipomoea subgenus Quamoclit (Convolvulaceae) based on ITS sequence data and a Bayesian phylogenetic analysis

A Bayesian phylogenetic analysis of 36 Ipomoea species using sequence data from the internal transcribed spacer region was compared with classification schemes based on traditional methods and a previously published cpDNA restriction fragment length polymorphism (RFLP) study. These molecular studies support a diversity of groups that were circumscribed on the basis of phenetic principles and agree generally with the results from cpDNA RFLP analyses. The congruence between the phylogenetic hypotheses based on new molecular data and the understanding of relationships developed in earlier studies indicate that these classifications may reflect evolutionary history. Two large clades of species, with one including sections Tricolores, Calonyction, and Pharbitis and the other including sections Mina and Leptocallis, were identified. Furthermore, morphologically distinct groups of Ipomoea species received support from the DNA sequence data. Indices of convergence for the Markov chain Monte Carlo (MCMC) in the Bayesian phylogenetic analysis were evaluated. A limited range of posterior probabilities for each node in the trees from a set of five MCMC samples provides a useful index of convergence. Bayesian node support values were generally higher than bootstrap values from a maximum parsimony analysis. This is consistent with the notion that these measures of support estimate different qualities of the data.     

Phylogeny of sea skaters, Halobates Eschscholtz (Hemiptera, Gerridae), based on mtDNA sequence and morphology

We examined phylogenetic relationships among halobatine water striders (Hemiptera, Gerridae) using molecular and morphological data. The molecular data set was 780 bp DNA sequence data from the 3' half of the mitochondria! gene encoding cytochrome oxidase subunit I from 19 species of sea skaters, Halobates , and one species from each of three related genera, Asclepios annandalei, Austrobates rivularis , and Eurymetra natalensis. The morphological data set was a slightly modified version of a previously published data set. Unweighted parsimony analyses of the molecular data set gave one tree with weak support for most branches. Maximum likelihood analysis of the same data set gave a tree with slightly different topology, but reveiled many of the clades found in parsimony analyses of the morphological data set. Parsimony analyses of the combined molecular + morphology data sets gave a better resolved and better supported tree than did analyses of any single data set. The phytogeny of Halobates presented here allows a more rigorous evaluation of several prior hypotheses about evolutionary processes in marine water striders. In particular, it supports the hypothesis of at least two separate transitions from coastal to oceanic environments.     

Phylogeny of the genera of the parasitic wasps subfamily Doryctinae (Hymenoptera: Braconidae) based on morphological evidence

The evolutionary relationships among most (143 genera) of the currently recognized genera of the braconid wasp subfamily Doryctinae were investigated using maximum parsimony analysis, employing 100 characters from external morphology and four additional, less well‐studied character systems (male genitalia, ovipositor structure, venom apparatus and larval cephalic structure). We investigated the ‘performance’ of characters from external morphology and the other character systems and the effects of abundant missing entries by comparing the data decisiveness, retention and consistency indices of four different character partitions. The results indicate that the performances of the different partitions are not related to the proportions of missing entries, but instead are negatively correlated to their proportion of informative characters, suggesting that the morphological information in this group is subject to high levels of homoplasy. The external morphological partition is significantly incongruent with respect to a data set comprising the other character systems based on the ILD test. Analyses supported neither the monophyly of the large tribes Doryctini and Hecabolini, nor the monophyly of the Spathiini and Heterospilini. Relationships obtained from successive approximation weighting analysis for the complete data differ considerably from the currently accepted tribal and subtribal classifications. The only exceptions were the Ypsistocerini and the Ecphylini, whose recognized members were recovered in single clades. A close relationship between the Binaerini and Holcobraconini, and also Monarea, is consistently supported by venom apparatus and ovipositor structure characters but is not indicated by external morphological data. Low bootstrap values obtained for most of the recovered clades in all analyses do not allow us to propose a meaningful reclassification for the group at this time. A complete list of the recognized genera and their synonymies is given. © 2004 The Linnean Society of London, Zoological Journal of the Linnean Society, 2004, 142 , 369–404.     

Elongation factor-2: a useful gene for arthropod phylogenetics.

Robust resolution of controversial higher-level groupings within Arthropoda requires additional sources of characters. Toward this end, elongation factor-2 sequences (1899 nucleotides) were generated from 17 arthropod taxa (5 chelicerates, 6 crustaceans, 3 hexapods, 3 myriapods) plus an onychophoran and a tardigrade as outgroups. Likelihood and parsimony analyses of nucleotide and amino acid data sets consistently recovered Myriapoda and major chelicerate groups with high bootstrap support. Crustacea + Hexapoda (= Pancrustacea) was recovered with moderate support, whereas the conflicting group Myriapoda + Hexapoda (= Atelocerata) was never recovered and bootstrap values were always <5%. With additional nonarthropod sequences included, one indel supports monophyly of Tardigrada, Onychophora, and Arthropoda relative to molluscan, annelidan, and mammalian outgroups. New and previously published sequences from RNA polymerase II (1038 nucleotides) and elongation factor-1alpha (1092 nucleotides) were analyzed for the same taxa. A comparison of bootstrap values from the three genes analyzed separately revealed widely varying values for some clades, although there was never strong support for conflicting groups. In combined analyses, there was strong bootstrap support for the generally accepted clades Arachnida, Arthropoda, Euchelicerata, Hexapoda, and Pycnogonida, and for Chelicerata, Myriapoda, and Pancrustacea, whose monophyly is more controversial. Recovery of some additional groups was fairly robust to method of analysis but bootstrap values were not high; these included Pancrustacea + Chelicerata, Hexapoda + Cephalocarida + Remipedia, Cephalocarida + Remipedia, and Malaocostraca + Cirripedia. Atelocerata (= Myriapoda + Hexapoda) was never recovered. Elongation factor-2 is now the second protein-encoding, nuclear gene (in addition to RNA polymerase II) to support Pancrustacea over Atelocerata. Atelocerata is widely cited in morphology-based analyses, and the discrepancy between results derived from molecular and morphological data deserves greater attention.     

Phylogenetic relationships of Steinernema Travassos, 1927 (Nematoda: Cephalobina: Steinernematidae) based on nuclear, mitochondrial and morphological data

Entomopathogenic nematodes of the genus Steinernema are lethal parasites of insects that are used as biological control agents of several lepidopteran, dipteran and coleopteran pests. Phylogenetic relationships among 25 Steinernema species were estimated using nucleotide sequences from three genes and 22 morphological characters. Parsimony analysis of 28S (LSU) sequences yielded a well-resolved phylogenetic hypothesis with reliable bootstrap support for 13 clades. Parsimony analysis of mitochondrial DNA sequences (12S rDNA and cox 1 genes) yielded phylogenetic trees with a lower consistency index than for LSU sequences, and with fewer reliably supported clades. Combined phylogenetic analysis of the 3-gene dataset by parsimony and Bayesian methods yielded well-resolved and highly similar trees. Bayesian posterior probabilities were high for most clades; bootstrap (parsimony) support was reliable for approximately half of the internal nodes. Parsimony analysis of the morphological dataset yielded a poorly resolved tree, whereas total evidence analysis (molecular plus morphological data) yielded a phylogenetic hypothesis consistent with, but less resolved than trees inferred from combined molecular data. Parsimony mapping of morphological characters on the 3-gene trees showed that most structural features of steinernematids are highly homoplastic. The distribution of nematode foraging strategies on these trees predicts that S. hermaphroditum, S. diaprepesi and S. longicaudum (US isolate) have cruise forager behaviours.     

Molecular phylogenetics and morphological evolution of St. John’s wort (Hypericum; Hypericaceae)

Phylogenetic hypotheses for the large cosmopolitan genus Hypericum (St. John’s wort) have previously been based on morphology, and molecular studies have thus far included only a few species. In this study, we used 360 sequences of the internal transcribed spacer (ITS) region of nuclear ribosomal DNA (nrDNA) for 206 species representing Hypericum (incl. Triadenum and Thornea) and three other genera of Hypericaceae to generate an explicit phylogenetic hypothesis for the genus using parsimony and model-based methods. The results indicate that the small genus Triadenum is nested in a clade within Hypericum containing most of the New World species. Sister to Hypericum is Thornea from Central America. Within Hypericum, three large clades and two smaller grades were found; these are based on their general morphology, especially characters used previously in taxonomy of the genus. Relative to the most recent classification, around 60% of the sections of Hypericum were monophyletic. We used a Bayesian approach to reconstruct ancestral states of selected morphological characters, which resulted in recognition of characters that support major clades within the genus and a revised interpretation of morphological evolution in Hypericum. The shrubby habit represents the plesiomorphic state from which herbs evolved several times. Arborescent species have radiated convergently in high-elevation habitats in tropical Africa and South America.     

Phylogenetic investigations of Antarctic notothenioid fishes (Perciformes: Notothenioidei) using complete gene sequences of the mitochondrial encoded 16S rRNA

The Notothenioidei dominates the fish fauna of the Antarctic in both biomass and diversity. This clade exhibits adaptations related to metabolic function and freezing avoidance in the subzero Antarctic waters, and is characterized by a high degree of morphological and ecological diversity. Investigating the macroevolutionary processes that may have contributed to the radiation of notothenioid fishes requires a well-resolved phylogenetic hypothesis. To date published molecular and morphological hypotheses of notothenioids are largely congruent, however, there are some areas of significant disagreement regarding higher-level relationships. Also, there are critical areas of the notothenioid phylogeny that are unresolved in both molecular and morphological phylogenetic analyses. Previous molecular phylogenetic analyses of notothenioids using partial mtDNA 12S and 16S rRNA sequence data have resulted in limited phylogenetic resolution and relatively low node support. One particularly controversial result from these analyses is the paraphyly of the Nototheniidae, the most diverse family in the Notothenioidei. It is unclear if the phylogenetic results from the 12S and 16S partial gene sequence dataset are due to limited character sampling, or if they reflect patterns of evolutionary diversification in notothenioids. We sequenced the complete mtDNA 16S rRNA gene for 43 notothenioid species, the largest sampling to-date from all eight taxonomically recognized families. Phylogenetic analyses using both maximum parsimony and maximum likelihood resulted in well-resolved trees with most nodes supported with high bootstrap pseudoreplicate scores and significant Bayesian posterior probabilities. In all analyses the Nototheniidae was monophyletic. Shimodaira–Hasegawa tests were able to reject two hypotheses that resulted from prior morphological analyses. However, despite substantial resolution and node support in the 16S rRNA trees, several phylogenetic hypotheses among closely related species and clades were not rejected. The inability to reject particular hypotheses among species in apical clades is likely due to the lower rate of nucleotide substitution in mtDNA rRNA genes relative to protein coding regions. Nevertheless, with the most extensive notothenioid taxon sampling to date, and the much greater phylogenetic resolution offered by the complete 16S rRNA sequences over the commonly used partial 12S and 16S gene dataset, it would be advantageous for future molecular investigations of notothenioid phylogenetics to utilize at the minimum the complete gene 16S rRNA dataset.