Likelihood Ancestral State Reconstruction Research Articles

Hip joint articular soft tissues of non-dinosaurian Dinosauromorpha and early Dinosauria: evolutionary and biomechanical implications for Saurischia

ABSTRACTDinosauromorphs evolved a wide diversity of hind limb skeletal morphologies, suggesting highly divergent articular soft tissue anatomies. However, poor preservation of articular soft tissues in fossils has hampered any follow-on functional inferences. We reconstruct the hip joint soft tissue anatomy of non-dinosaurian dinosauromorphs and early dinosaurs using osteological correlates derived from extant sauropsids and infer trends in character transitions along the theropod and sauropodomorph lineagues. Femora and pelves of 107 dinosauromorphs and outgroup taxa were digitized using 3D imaging techniques. Key transitions were estimated using maximum likelihood ancestral state reconstruction. The hips of dinosauromorphs possessed wide a disparity of soft tissue morphologies beyond the types and combinations exhibited by extant archosaurs. Early evolution of the dinosauriform hip joint was characterized by the retention of a prominent femoral hyaline cartilage cone in post-neonatal individuals, with the cartilage cone independently reduced within theropods and sauropodomorphs. The femur of Dinosauriformes possessed a fibrocartilage sleeve on the metaphysis, which surrounded a hyaline core. The acetabulum of Dinosauriformes possessed distinct labrum and antitrochanter structures. In sauropodomorphs, hip congruence was maintained by thick hyaline cartilage on the femoral head, whereas theropods relied on acetabular tissues such as ligaments and articular pads. In particular, the craniolaterally ossified hip capsule of non-Avetheropoda neotheropods permitted mostly parasagittal femoral movements. These data indicate that the dinosauromorph hip underwent mosaic evolution within the saurischian lineage and that sauropodomorphs and theropods underwent both convergence and divergence in articular soft tissues, correlated with transitions in body size, locomotor posture, and joint loading.SUPPLEMENTAL DATA—Supplemental materials are available for this article for free at www.tandfonline.com/UJVPCitation for this article: Tsai, H. P., K. M. Middleton, J. R. Hutchinson, and C. M. Holliday. 2018. Hip joint articular soft tissues of non-dinosaurian Dinosauromorpha and early Dinosauria: evolutionary and biomechanical implications for Saurischia. Journal of Vertebrate Paleontology. DOI: 10.1080/02724634.2017.1427593.

Phylogeny of the order Phyllachorales (Ascomycota, Sordariomycetes): among and within order relationships based on five molecular loci.

The order Phyllachorales (Pezizomycotina, Ascomycota) is a group of biotrophic, obligate plant parasitic fungi with a tropical distribution and high host specificity. Traditionally two families are recognised within this order: Phyllachoraceae and Phaeochoraceae, based mostly on morphological and host characteristics. Currently, the position of the order within the class Sordariomycetes is inconclusive, as well as the monophyly of the order, and its internal phylogenetic structure. Here we present a phylogeny of the order Phyllachorales based on sequence data of 29 species with a broad host range resulting from a wide geographical sampling. We inferred Maximum Likelihood and Bayesian phylogenies from data of five DNA regions: nrLSU rDNA, nrSSU rDNA, ITS rDNA, and the protein coding genes RPB2, and TEF1. We found that the order Phyllachorales is monophyletic and related to members of the subclass Sordariomycetidae within Sordariomycetes. Within the order, members of the family Phaeochoraceae form a monophyletic group, and the family Phyllachoraceae is split into two lineages. Maximum Likelihood ancestral state reconstructions indicate that the ancestor of Phyllachorales had a monocotyledonous host plant, immersed perithecia, and a black stroma. Alternative states of these characters evolved multiple times independently within the order. Based on our results we redefine the family Phyllachoraceae and propose the new family Telimenaceae with Telimena erythrinae as type species, resulting in three families in the order. Species of Telimena spp. occur in several monocotyledonous and eudicotyledonous host plants except Poaceae, and generally have enlarged black pseudostroma around the perithecia, a character not present in species of Phyllachoraceae.

Mosasauroid phylogeny under multiple phylogenetic methods provides new insights on the evolution of aquatic adaptations in the group.

Mosasauroids were a successful lineage of squamate reptiles (lizards and snakes) that radiated during the Late Cretaceous (95–66 million years ago). They can be considered one of the few lineages in the evolutionary history of tetrapods to have acquired a fully aquatic lifestyle, similarly to whales, ichthyosaurs and plesiosaurs. Despite a long history of research on this group, their phylogenetic relationships have only been tested so far using traditional (unweighted) maximum parsimony. However, hypotheses of mosasauroid relationships and the recently proposed multiple origins of aquatically adapted pelvic and pedal features in this group can be more thoroughly tested by methods that take into account variation in branch lengths and evolutionary rates. In this study, we present the first mosasauroid phylogenetic analysis performed under different analytical methods, including maximum likelihood, Bayesian inference, and implied weighting maximum parsimony. The results indicate a lack of congruence in the topological position of halisaurines and Dallasaurus. Additionally, the genus Prognathodon is paraphyletic under all hypotheses. Interestingly, a number of traditional mosasauroid clades become weakly supported, or unresolved, under Bayesian analyses. The reduced resolutions in some consensus trees create ambiguities concerning the evolution of fully aquatic pelvic/pedal conditions under many analyses. However, when enough resolution was obtained, reversals of the pelvic/pedal conditions were favoured by parsimony and likelihood ancestral state reconstructions instead of independent origins of aquatic features in mosasauroids. It is concluded that most of the observed discrepancies among the results can be associated with different analytical procedures, but also due to limited postcranial data on halisaurines, yaguarasaurines and Dallasaurus.

Lower Ilium Evolution in Apes and Hominins.

Elucidating the pelvic morphology of the Pan-Homo last common ancestor (LCA) is crucial for understanding ape and human evolution. The pelvis of Ardipithecus ramidus has been the basis of controversial interpretations of the LCA pelvis. In particular, it was proposed that the lower ilium became elongate independently in the orangutan and chimpanzee clades, making these taxa poor analogues for the pelvis of the LCA. This study examines the variation in relative lower ilium height between and within living and fossil hominoid species (and other anthropoids), and models its evolution using available fossil hominoids as calibration points. We find nuanced differences in relative lower ilium height among living hominoids, particularly in regards to gorillas, which do not have elongate lower ilia (because they are likely to represent the plesiomorphic hominoid condition for this trait). We also show that differences in relative lower ilium height among hominoid taxa are not readily explained by differences in size between species. Our maximum likelihood ancestral state reconstructions support inferences that chimpanzees (Pan troglodytes in particular) and orangutans evolved their elongate lower ilia independently. We also find that the predicted lower ilium height of the Pan-Homo LCA is shorter than all great apes except gorillas. This study adds to a growing body of evidence that finds different regions of the body show different evolutionary histories in different hominoids, and underscores that the unique combinations of morphologies of each modern and fossil hominoid species should be considered when reconstructing the mosaic nature of the Pan-Homo LCA. Anat Rec, 300:828-844, 2017. © 2017 Wiley Periodicals, Inc.

Rapid maximum likelihood ancestral state reconstruction of continuous characters: A rerooting-free algorithm.

Ancestral state reconstruction is a method used to study the evolutionary trajectories of quantitative characters on phylogenies. Although efficient methods for univariate ancestral state reconstruction under a Brownian motion model have been described for at least 25 years, to date no generalization has been described to allow more complex evolutionary models, such as multivariate trait evolution, non‐Brownian models, missing data, and within‐species variation. Furthermore, even for simple univariate Brownian motion models, most phylogenetic comparative R packages compute ancestral states via inefficient tree rerooting and full tree traversals at each tree node, making ancestral state reconstruction extremely time‐consuming for large phylogenies. Here, a computationally efficient method for fast maximum likelihood ancestral state reconstruction of continuous characters is described. The algorithm has linear complexity relative to the number of species and outperforms the fastest existing R implementations by several orders of magnitude. The described algorithm is capable of performing ancestral state reconstruction on a 1,000,000‐species phylogeny in fewer than 2 s using a standard laptop, whereas the next fastest R implementation would take several days to complete. The method is generalizable to more complex evolutionary models, such as phylogenetic regression, within‐species variation, non‐Brownian evolutionary models, and multivariate trait evolution. Because this method enables fast repeated computations on phylogenies of virtually any size, implementation of the described algorithm can drastically alleviate the computational burden of many otherwise prohibitively time‐consuming tasks requiring reconstruction of ancestral states, such as phylogenetic imputation of missing data, bootstrapping procedures, Expectation‐Maximization algorithms, and Bayesian estimation. The described ancestral state reconstruction algorithm is implemented in the Rphylopars functions anc.recon and phylopars.

Cytology and molecular phylogenetics of Monoblepharidomycetes provide evidence for multiple independent origins of the hyphal habit in the Fungi

The evolution of filamentous hyphae underlies an astounding diversity of fungal form and function. We studied the cellular structure and evolutionary origins of the filamentous form in the Monoblepharidomycetes (Chytridiomycota), an early-diverging fungal lineage that displays an exceptional range of body types, from crescent-shaped single cells to sprawling hyphae. To do so, we combined light and transmission electron microscopic analyses of hyphal cytoplasm with molecular phylogenetic reconstructions. Hyphae of Monoblepharidomycetes lack a complex aggregation of secretory vesicles at the hyphal apex (i.e. Spitzenkörper), have centrosomes as primary microtubule organizing centers and have stacked Golgi cisternae instead of tubular/fenestrated Golgi equivalents. The cytoplasmic distribution of actin in Monoblepharidomycetes is comparable to the arrangement observed previously in other filamentous fungi. To discern the origins of Monoblepharidomycetes hyphae, we inferred a phylogeny of the fungi based on 18S and 28S ribosomal DNA sequence data with maximum likelihood and Bayesian inference methods. We focused sampling on Monoblepharidomycetes to infer intergeneric relationships within the class and determined 78 new sequences. Analyses showed class Monoblepharidomycetes to be monophyletic and nested within Chytridiomycota. Hyphal Monoblepharidomycetes formed a clade sister to the genera without hyphae, Harpochytrium and Oedogoniomyces. A likelihood ancestral state reconstruction indicated that hyphae arose independently within the Monoblepharidomycetes lineage and in at least two other lineages. Cytological differences among monoblepharidalean and other fungal hyphae are consistent with these convergent origins.

Comparative Morphology and Evolution of Genital Papillae in a Genus of Darters (Percidae:Etheostoma)

A morphological study of female genital papillae of Etheostoma was conducted to examine the evolution of papillae as predicted by molecular phylogenies, and to examine correlations between papilla morphology and spawning behaviors and oviposition. Thirteen characters were used to describe variation in papillae in 128 species. Based on maximum likelihood ancestral state reconstructions, we corroborated several molecular phylogenetic hypotheses with morphology and developed hypotheses on several functional traits/adaptive traits and their associations with spawning behaviors. Papillar synapomorphies supported a close relationship between Allohistium and (Etheostoma + Nanostoma) + Ulocentra, a monophyletic Catonotus, and a sister-group relationship between E. vitreum and Boleosoma. Results suggest that the most recent common ancestor of Etheostoma buried eggs and had a simple tube with a distal, posteriorly oriented genital pore. All egg-burying species, except species of Nothonotus, have tube papillae that ...

Complete forelimb myology of the basal theropod dinosaur Tawa hallae based on a novel robust muscle reconstruction method.

The forelimbs of nonavian theropod dinosaurs have been the subject of considerable study and speculation due to their varied morphology and role in the evolution of flight. Although many studies on the functional morphology of a limb require an understanding of its musculature, comparatively little is known about the forelimb myology of theropods and other bipedal dinosaurs. Previous phylogenetically based myological reconstructions have been limited to the shoulder, restricting their utility in analyses of whole-limb function. The antebrachial and manual musculature in particular have remained largely unstudied due to uncertain muscular homologies in archosaurs. Through analysis of the musculature of extant taxa in a robust statistical framework, this study presents new hypotheses of homology for the distal limb musculature of archosaurs and provides the first complete reconstruction of dinosaurian forelimb musculature, including the antebrachial and intrinsic manual muscles. Data on the forelimb myology of a broad sample of extant birds, crocodylians, lizards, and turtles were analyzed using maximum likelihood ancestral state reconstruction and examined together with the osteology of the early theropod Tawa hallae from the Late Triassic of New Mexico to formulate a complete plesiomorphic myology for the theropod forelimb. Comparisons with previous reconstructions show that the shoulder musculature of basal theropods is more similar to that of basal ornithischians and sauropodomorphs than to that of dromaeosaurids. Greater development of the supracoracoideus and deltoideus musculature in theropods over other bipedal dinosaurs correlates with stronger movements of the forelimb at the shoulder and an emphasis on apprehension of relatively large prey. This emphasis is further supported by the morphology of the antebrachium and the intrinsic manual musculature, which exhibit a high degree of excursion and a robust morphology well-suited for powerful digital flexion. The forelimb myology of Tawa established here helps infer the ancestral conformation of the forelimb musculature and the osteological correlates of major muscle groups in early theropods. These data are critical for investigations addressing questions relating to the evolution of specialized forelimb function across Theropoda.

Water barriers and intra‐island isolation contribute to diversification in the insular Aethopyga sunbirds (Aves: Nectariniidae)

AbstractAimColonization and subsequent isolation across deep‐water barriers is thought to be the primary driver of diversification in insular birds. Shallow‐water barriers and intra‐island isolation are less well‐documented drivers of avian diversification. We examined the relative roles of different geographical barriers in the diversification of Aethopyga sunbirds, a widespread Southeast Asian genus that has its greatest diversity in the Philippine Archipelago.LocationPhilippines, Southeast Asia, Wallacea.MethodsWe reconstructed the phylogenetic relationships among Aethopyga sunbirds with mitochondrial and nuclear DNA sequences. Phylogeny was inferred using concatenated and coalescent frameworks, implemented in maximum likelihood and Bayesian analyses. We used maximum likelihood ancestral state reconstructions to examine the ancestral distribution and colonization history of Aethopyga. To determine whether the diversification at each node occurred within a continent, across a shallow‐water barrier, across a deep‐water barrier or within an island, we used a series of statements based on the phylogeny, current distribution of species and bathymetric reconstructions.ResultsAncestral state reconstructions inferred that the core Aethopyga ancestor was continental, and that the diversity of Aethopyga on oceanic islands is the result of three or four independent colonization events. Dispersal and subsequent isolation across deep‐water barriers was the most common mode of diversification in insular Aethopyga, although intra‐island isolation contributed to diversity, producing a small montane radiation within Mindanao. Analyses inferred only a single unequivocal event of diversification across a shallow‐water barrier. Deep molecular divergences between phenotypically distinct subspecies suggested that the taxonomy of Aethopyga is overly conservative and obscures biogeographical patterns. We recommend elevating five subspecies, all of which are endemic to the Philippines, to full species.Main conclusionsIn addition to corroborating the importance of isolation across deep‐water barriers, these data also underscore a potential role for additional isolating mechanisms in the generation of biodiversity on oceanic islands. A clade of four montane Aethopyga taxa evolved in situ within the Mindanao sky‐island system, an exception to the paradigm that birds do not diversify within the geographical confines of oceanic islands.

ON THE REVERSIBILITY OF MANDIBULAR SYMPHYSEAL FUSION

Experimental and comparative studies suggest that a major determinant of increased ossification of the mandibular symphysis is elevated masticatory stress related to a mechanically challenging diet. However, the morphology of this joint tracks variation in dietary properties in only some mammalian clades. Extant anthropoid primates are a notable exception: synostosis is ubiquitous in this speciose group, despite its great age and diverse array of feeding adaptations. One possible explanation for this pattern is that, once synostosis evolves, reversion to a lesser degree of fusion is unlikely or even constrained. If correct, this has important implications for functional and phylogenetic analyses of the mammalian feeding apparatus. To test this hypothesis, we generated a molecular tree for 76 vespertilionoid and noctilionoid chiropterans using Bayesian phylogenetic analysis and examined character evolution using parsimony and likelihood ancestral-state reconstructions along with the binary state speciation and extinction (BiSSE) model. Results indicate that reversals have occurred within Vespertilionoidea. In contrast, noctilionoids exhibit an anthropoid-like pattern, which suggests that more detailed comparisons of the functional and developmental bases for fusion in these bat clades may provide insight into why fusion is maintained in some lineages but not in others. Potential functional and developmental explanations for the lack of reversal are discussed.

Evolution of anatomical and physiological specialization in the compound eyes of stomatopod crustaceans

Stomatopod crustaceans have complex and diverse visual systems. Among their many unique features are a specialized ommatidial region (the midband) that enables the eye to have multiple overlapping visual fields, as well as sets of spectral filters that are intercalated at two levels between tiers of photoreceptors involved in polychromatic color vision. Although the physiology and visual function of stomatopod eyes have been studied for many years, how these unique visual features originated and diversified is still an open question. In order to investigate how stomatopods have attained the current complexity in visual function, we have combined physiological and morphological information (e.g. number of midband rows, number of filters in the retina, and the spectral properties of filters) with new phylogenetic analyses of relationships among species based on nucleotide sequence data from two nuclear (18S and 28S rDNA) and two mitochondrial [16S and cytochrome oxidase I (COI)] genes. Based on our recovered phylogenetic relationships among species, we propose two new superfamilies within the Stomatopoda: Hemisquilloidea and Pseudosquillodea. Maximum likelihood ancestral state reconstructions indicate that ancestral stomatopod eyes contained six midband rows and four intrarhabdomal filters, illustrating that the visual physiological complexity originated early in stomatopod evolutionary history. While the two distal filters contain conservative sets of filter pigments, the proximal filters show more spectral diversity in filter types, particularly in midband row 2, and are involved in tuning the color vision system to the photic environment. In particular, a set of related gonodactyloid families (Gonodactylidae, Protosquillidae, Takuidae) inhabiting shallow, brightly lit coral reef waters contain the largest diversity of filter pigments, which are spectrally placed relative to the underlying photoreceptors to take advantage of the broad spectrum of light available in the environment.

Molecular systematics of tribe Rubieae (Rubiaceae): Evolution of major clades, development of leaf–like whorls, and biogeography

AbstractRubieae are centered in temperate regions and characterized by whorls of leaf–like structures on their stems. Previous studies that primarily included Old World taxa identified seven major clades with no resolution between and within clades. In this study, a molecular phylogeny of the tribe, based on three chloroplast regions (rpoB–trnC, trnC–psbM, trnL–trnF–ndhJ) from 126 Old and New World taxa, is estimated using parsimony and Bayesian analyses. Seven major clades are strongly supported within the tribe, confirming previous studies. Relationships within and between these seven major clades are also strongly supported. In addition, the position of Callipeltis, a previously unsampled genus, is identified. The resulting phylogeny is used to examine geographic distribution patterns and evolution of leaf–like whorls in the tribe. An Old World origin of the tribe is inferred from parsimony and likelihood ancestral state reconstructions. At least eight subsequent dispersal events into North America occurred from Old World ancestors. From one of these dispersal events, a radiation into North America, followed by subsequent diversification in South America, occurred. Parsimony and likelihood ancestral state reconstructions infer the ancestral whorl morphology of the tribe as composed of six organs. Whorls composed of four organs are derived from whorls with six or more organs. Transitions between four and six or more organs per whorl are common within the tribe, whereas reduction to two leaves at a node is derived and rare.

Evolution of lamina anatomy in the palm family (Arecaceae)

The unique properties of tree building in Arecaceae strongly constrain their architectural lability. Potentially compensating for this limitation, the extensive diversification of leaf anatomical structure within palms involves many characters whose alternate states may confer disparate mechanical or physiological capabilities. In the context of a recent global palm phylogeny, we analyzed the evolution of 10 such lamina anatomical characters and leaf morphology of 161 genera, conducting parsimony and maximum likelihood ancestral state reconstructions, as well as tests of correlated evolution. Lamina morphology evolves independently from anatomy. Although many characters do optimize as synapomorphic for major clades, anatomical evolution is highly homoplasious. Nevertheless, it is not random: analyses indicate the recurrent evolution of different cohorts of correlated character states. Notable are two surface layer (epidermis and hypodermis) types: (1) a parallel-laminated type of rectangular epidermal cells with sinuous anticlinal walls, with fibers present in the hypodermis and (2) a cross-laminated type of hexagonal cells in both layers. Correlated with the cross-laminated type is a remarkable decrease in the volume fraction of fibers, accompanied by changes in the architecture and sheath cell type of the transverse veins. We discuss these and other major patterns of anatomical evolution in relation to their biomechanical and ecophysiological significance.

Molecular phylogeny of the Tylenchina and evolution of the female gonoduct (Nematoda: Rhabditida)

Tylenchina are a morphologically and functionally diverse group of nematode species that range from free-living bacteriovores, over transitory grazing root-hair feeders to highly specialized plant-parasites with complex host associations. We performed phylogenetic analyses of small subunit rDNA sequences from 97 species including an analysis that account for the RNA secondary structure in the models of evolution. The present study confirms the sister relationship of the bacteriovore Cephalobidae with the predominantly plant-parasitic Tylenchomorpha. All analyses appoint the fungal-feeding Aphelenchidae and Aphelenchoididae as being polyphyletic but the morphology based hypothesis of their monophyly could not be significantly rejected. Within the Tylenchomorpha, the families that exclusively parasitize higher plants are joined in a single clade. However, only the monophyletic position of the (super)families Hoplolaimidae and Criconematoidea were supported; Anguinidae, Tylenchidae, Belonolaimidae and Pratylenchidae appeared to be paraphyletic or polyphyletic. Parsimony and likelihood ancestral state reconstruction revealed that burrowing endoparasitism and sedentary endoparasitism each evolved, respectively, at least six and at least three times independently, mostly from migratory ectoparasitic ancestors. Only root-knot nematodes have evolved from burrowing endoparasitic nematodes. Traditional classifications are partially misled by this convergent evolution of feeding type and associated morphology. Contrastingly, mapping attributes of the gonoduct cellular architecture, including newly obtained data of 18 species belonging to the Aphelenchoidea, Criconematoidea, Anguinidae and Panagrolaimidae, revealed a broad congruence of the gonoduct characters and the molecular phylogenetic hypothesis. Yet, the presence of an offset spermatheca and proliferation of uterus cells has evolved multiple times, the latter associated with derived endoparasitic feeding specialization and resulting reproduction mode. Ancestral state reconstruction further revealed that the gonoduct of the morphologically and ecologically dissimilar tylenchid and cephalobid nematodes evolved from a common ancestor.

Gene Duplication Is an Evolutionary Mechanism for Expanding Spectral Diversity in the Long-Wavelength Photopigments of Butterflies

Butterfly long-wavelength (L) photopigments are interesting for comparative studies of adaptive evolution because of the tremendous phenotypic variation that exists in their wavelength of peak absorbance (lambda(max) value). Here we present a comprehensive survey of L photopigment variation by measuring lambda(max) in 12 nymphalid and 1 riodinid species using epi-microspectrophotometry. Together with previous data, we find that L photopigment lambda(max) varies from 510-565 nm in 22 nymphalids, with an even broader 505- to 600-nm range in riodinids. We then surveyed the L opsin genes for which lambda(max) values are available as well as from related taxa and found 2 instances of L opsin gene duplication within nymphalids, in Hermeuptychia hermes and Amathusia phidippus, and 1 instance within riodinids, in the metalmark butterfly Apodemia mormo. Using maximum parsimony and maximum likelihood ancestral state reconstructions to map the evolution of spectral shifts within the L photopigments of nymphalids, we estimate the ancestral pigment had a lambda(max) = 540 nm +/- 10 nm standard error and that blueshifts in wavelength have occurred at least 4 times within the family. We used ancestral state reconstructions to investigate the importance of several amino acid substitutions (Ile17Met, Ala64Ser, Asn70Ser, and Ser137Ala) previously shown to have evolved under positive selection that are correlated with blue spectral shifts. These reconstructions suggest that the Ala64Ser substitution has indeed occurred along the newly identified blueshifted L photopigment lineages. Substitutions at the other 3 sites may also be involved in the functional diversification of L photopigments. Our data strongly suggest that there are limits to the evolution of L photopigment spectral shifts among species with only one L opsin gene and that opsin gene duplication broadens the potential range of lambda(max) values.

A Genomic Fossil Reveals Key Steps in Hemoglobin Loss by the Antarctic Icefishes

Antarctic icefishes are the only vertebrates that do not have hemoglobin and erythrocytes in their blood. These startling phenotypes are associated in several icefish species with deletions of juvenile and adult globin loci, which in red-blooded teleosts are typically composed of tightly linked pairs of alpha- and beta-globin genes. It is unknown if the loss of hemoglobin expression in icefishes was the direct result of such deletions or if other mutational events compromised globin chain synthesis prior to globin gene loss. In this study, we show that 15 of the 16 icefish species have lost the adult beta-globin gene but retain a truncated alpha-globin pseudogene. Surprisingly, a phylogenetically derived icefish species, Neopagetopsis ionah, possesses a complete, but nonfunctional, adult alphabeta-globin complex. This cluster contains 2 distinct beta-globin pseudogenes whose phylogenetic origins span the entire Antarctic notothenioid radiation, consistent with an origin via introgression. Maximum likelihood ancestral state reconstruction supports a scenario of icefish globin gene evolution that involves a single loss of the transcriptionally active adult alphabeta-globin cluster prior to the diversification of the extant species in the clade. Through lineage sorting of ancestral polymorphism, 2 types of alleles became fixed in the clade: 1) the alpha-globin pseudogene of the majority of species and 2) the inactive alphabeta-globin complex of N. ionah. We conclude that the globin pseudogene complex of N. ionah is a "genomic fossil" that reveals key intermediate steps on the pathway to loss of hemoglobin expression by all icefish species.