Partitioned Gene-Tree Analyses and Gene-Based Topology Testing Help Resolve Incongruence in a Phylogenomic Study of Host-Specialist Bees (Apidae: Eucerinae)

Incongruence among phylogenetic results has become a common occurrence in analyses of genome-scale data sets. Incongruence originates from uncertainty in underlying evolutionary processes (e.g., incomplete lineage sorting) and from difficulties in determining the best analytical approaches for each situation. To overcome these difficulties, more studies are needed that identify incongruences and demonstrate practical ways to confidently resolve them. Here, we present results of a phylogenomic study based on the analysis 197 taxa and 2,526 ultraconserved element (UCE) loci. We investigate evolutionary relationships of Eucerinae, a diverse subfamily of apid bees (relatives of honey bees and bumble bees) with >1,200 species. We sampled representatives of all tribes within the group and >80% of genera, including two mysterious South American genera, Chilimalopsis and Teratognatha. Initial analysis of the UCE data revealed two conflicting hypotheses for relationships among tribes. To resolve the incongruence, we tested concatenation and species tree approaches and used a variety of additional strategies including locus filtering, partitioned gene-trees searches, and gene-based topological tests. We show that within-locus partitioning improves gene tree and subsequent species-tree estimation, and that this approach, confidently resolves the incongruence observed in our data set. After exploring our proposed analytical strategy on eucerine bees, we validated its efficacy to resolve hard phylogenetic problems by implementing it on a published UCE data set of Adephaga (Insecta: Coleoptera). Our results provide a robust phylogenetic hypothesis for Eucerinae and demonstrate a practical strategy for resolving incongruence in other phylogenomic data sets.     

Ultraconserved elements show utility in phylogenetic inference of Adephaga (Coleoptera) and suggest paraphyly of ‘Hydradephaga’

The beetle suborder Adephaga has been the subject of many phylogenetic reconstructions utilizing a variety of data sources and inference methods. However, no strong consensus has yet emerged on the relationships among major adephagan lineages. Ultraconserved elements (UCEs) have proved useful for inferring difficult or unresolved phylogenies at varying timescales in vertebrates, arachnids and Hymenoptera. Recently, a UCE bait set was developed for Coleoptera using polyphagan genomes and a member of the order Strepsiptera as an outgroup. Here, we examine the utility of UCEs for reconstructing the phylogeny of adephagan families, in the first in vitro application a UCE bait set in Coleoptera. Our final dataset included 305 UCE loci for 18 representatives of all adephagan families except Aspidytidae, and two polyphagan outgroups, with a total concatenated length of 83 547 bp. We inferred trees using maximum likelihood analyses of the concatenated UCE alignment and coalescent species tree methods (astral ii , ASTRID, svdquartets ). Although the coalescent species tree methods had poor resolution and weak support, concatenated analyses produced well‐resolved, highly supported trees. Hydradephaga was recovered as paraphyletic, with Gyrinidae sister to Geadephaga and all other adephagans. Haliplidae was recovered as sister to Dytiscoidea, with Hygrobiidae and Amphizoidae successive sisters to Dytiscidae. Finally, Noteridae was recovered as monophyletic and sister to Meruidae. Given the success of UCE data for resolving phylogenetic relationships within Adephaga, we suggest the potential for further resolution of relationships within Adephaga using UCEs with improved taxon sampling, and by developing Adephaga‐specific probes.     

Advancing mite phylogenomics: Designing ultraconserved elements for Acari phylogeny

Mites (Acari) are one of the most diverse groups of life on Earth; yet, their evolutionary relationships are poorly understood. Also, the resolution of broader arachnid phylogeny has been hindered by an underrepresentation of mite diversity in phylogenomic analyses. To further our understanding of Acari evolution, we design targeted ultraconserved genomic elements (UCEs) probes, intended for resolving the complex relationships between mite lineages and closely related arachnids. We then test our Acari UCE baits in‐silico by constructing a phylogeny using 13 existing Acari genomes, as well as 6 additional taxa from a variety of genomic sources. Our Acari‐specific probe kit improves the recovery of loci within mites over an existing general arachnid UCE probe set. Our initial phylogeny recovers the major mite lineages, yet finds mites to be non‐monophyletic overall, with Opiliones (harvestmen) and Ricinuleidae (hooded tickspiders) rendering Parasitiformes paraphyletic.     

Sequence capture and next‐generation sequencing of ultraconserved elements in a large‐genome salamander

Amidst the rapid advancement in next‐generation sequencing (NGS) technology over the last few years, salamanders have been left behind. Salamanders have enormous genomes—up to 40 times the size of the human genome—and this poses challenges to generating NGS data sets of quality and quantity similar to those of other vertebrates. However, optimization of laboratory protocols is time‐consuming and often cost prohibitive, and continued omission of salamanders from novel phylogeographic research is detrimental to species facing decline. Here, we use a salamander endemic to the southeastern United States, Plethodon serratus, to test the utility of an established protocol for sequence capture of ultraconserved elements (UCEs) in resolving intraspecific phylogeographic relationships and delimiting cryptic species. Without modifying the standard laboratory protocol, we generated a data set consisting of over 600 million reads for 85 P. serratus samples. Species delimitation analyses support recognition of seven species within P. serratus sensu lato, and all phylogenetic relationships among the seven species are fully resolved under a coalescent model. Results also corroborate previous data suggesting nonmonophyly of the Ouachita and Louisiana regions. Our results demonstrate that established UCE protocols can successfully be used in phylogeographic studies of salamander species, providing a powerful tool for future research on evolutionary history of amphibians and other organisms with large genomes.     

Replicated divergence in cichlid radiations mirrors a major vertebrate innovation

Decoupling of the upper jaw bones—jaw kinesis—is a distinctive feature of the ray-finned fishes, but it is not clear how the innovation is related to the extraordinary diversity of feeding behaviours and feeding ecology in this group. We address this issue in a lineage of ray-finned fishes that is well known for its ecological and functional diversity—African rift lake cichlids. We sequenced ultraconserved elements to generate a phylogenomic tree of the Lake Tanganyika and Lake Malawi cichlid radiations. We filmed a diverse array of over 50 cichlid species capturing live prey and quantified the extent of jaw kinesis in the premaxillary and maxillary bones. Our combination of phylogenomic and kinematic data reveals a strong association between biting modes of feeding and reduced jaw kinesis, suggesting that the contrasting demands of biting and suction feeding have strongly influenced cranial evolution in both cichlid radiations.     

Ultraconserved elements anchor thousands of genetic markers spanning multiple evolutionary timescales

Although massively parallel sequencing has facilitated large-scale DNA sequencing, comparisons among distantly related species rely upon small portions of the genome that are easily aligned. Methods are needed to efficiently obtain comparable DNA fragments prior to massively parallel sequencing, particularly for biologists working with non-model organisms. We introduce a new class of molecular marker, anchored by ultraconserved genomic elements (UCEs), that universally enable target enrichment and sequencing of thousands of orthologous loci across species separated by hundreds of millions of years of evolution. Our analyses here focus on use of UCE markers in Amniota because UCEs and phylogenetic relationships are well-known in some amniotes. We perform an in silico experiment to demonstrate that sequence flanking 2030 UCEs contains information sufficient to enable unambiguous recovery of the established primate phylogeny. We extend this experiment by performing an in vitro enrichment of 2386 UCE-anchored loci from nine, non-model avian species. We then use alignments of 854 of these loci to unambiguously recover the established evolutionary relationships within and among three ancient bird lineages. Because many organismal lineages have UCEs, this type of genetic marker and the analytical framework we outline can be applied across the tree of life, potentially reshaping our understanding of phylogeny at many taxonomic levels.     

Faster Protein Splicing with the Nostoc punctiforme DnaE Intein Using Non-native Extein Residues

Inteins are naturally occurring intervening sequences that catalyze a protein splicing reaction resulting in intein excision and concatenation of the flanking polypeptides (exteins) with a native peptide bond. Inteins display a diversity of catalytic mechanisms within a highly conserved fold that is shared with hedgehog autoprocessing proteins. The unusual chemistry of inteins has afforded powerful biotechnology tools for controlling enzyme function upon splicing and allowing peptides of different origins to be coupled in a specific, time-defined manner. The extein sequences immediately flanking the intein affect splicing and can be defined as the intein substrate. Because of the enormous potential complexity of all possible flanking sequences, studying intein substrate specificity has been difficult. Therefore, we developed a genetic selection for splicing-dependent kanamycin resistance with no significant bias when six amino acids that immediately flanked the intein insertion site were randomized. We applied this selection to examine the sequence space of residues flanking the Nostoc punctiforme Npu DnaE intein and found that this intein efficiently splices a much wider range of sequences than previously thought, with little N-extein specificity and only two important C-extein positions. The novel selected extein sequences were sufficient to promote splicing in three unrelated proteins, confirming the generalizable nature of the specificity data and defining new potential insertion sites for any target. Kinetic analysis showed splicing rates with the selected exteins that were as fast or faster than the native extein, refuting past assumptions that the naturally selected flanking extein sequences are optimal for splicing.     

Parallel evolutionary trajectories underlie the origin of giant suspension-feeding whales and bony fishes

Giant suspension feeders such as mysticete whales, basking and whale sharks, and the extinct (indicated by ‘†’) †pachycormiform teleosts are conspicuous members of modern and fossil marine vertebrate faunas. Whether convergent anatomical features common to these clades arose along similar evolutionary pathways has remained unclear because of a lack of information surrounding the origins of all groups of large-bodied suspension feeders apart from baleen whales. New investigation reveals that the enigmatic ray-finned fish †Ohmdenia, from the Lower Jurassic (Toarcian, 183.0–175.6 Ma) Posidonia Shale Lagerstätte, represents the immediate sister group of edentulous †pachycormiforms, the longest lived radiation of large vertebrate suspension feeders. †Ohmdenia bisects the long morphological branch leading to suspension-feeding †pachycormiforms, providing information on the sequence of anatomical transformations preceding this major ecological shift that can be compared to changes associated with the origin of modern mysticetes. Similarities include initial modifications to jaw geometry associated with the reduction of dentition, followed by the loss of teeth. The evolution of largest body sizes within both radiations occurs only after the apparent onset of microphagy. Comparing the fit of contrasting evolutionary models to functionally relevant morphological measurements for whales and †pachycormiform fishes reveals strong support for a common adaptive peak shared by suspension-feeding members of both clades.     

Use of Comparative Genomics to Develop EST-SSRs for Red Drum (Sciaenops ocellatus)

Microsatellites physically linked to expressed sequence tags (EST-SSRs) are an important resource for linkage mapping and comparative genomics, and data mining in publicly available EST databases is a common strategy for EST-SSR discovery. At present, many species lack species-specific EST sequence data needed for the efficient characterization of EST-SSRs. This paper describes the discovery and development of EST-SSRs for red drum (Sciaenops ocellatus), an estuarine-dependent sciaenid species of economic importance in the USA and elsewhere, using a phylogenetically informed, comparative genomics approach to primer design. The approach entailed comparing existing genomic resources from species closely allied phylogenetically to red drum, with resources from more distantly related outgroup species. By taking into account the degree to which flanking regions are conserved across taxa, the efficiency of PCR primer design was increased greatly. The amplification success rate for primers designed for red drum was 100?% when using EST libraries from confamilial species and 92?% when using an EST library from a species in the same suborder. The primers developed also amplified EST-SSRs in a wide range of perciform fishes, suggesting potential use in comparative genomics. This study demonstrates that EST-SSRs can be efficiently developed for an organism when limited species-specific data are available by exploiting genomic resources from well-studied species, even those at extended phylogenetic distances.     

Evolution of gastrulation in the ray-finned (actinopterygian) fishes

Sometime before or during the early Mesozoic era, new lineages of actinopterygian (ray-finned) fishes radically transformed their mode of gastrulation. During this evolutionary transformation, yolky endoderm was a hotspot for ontogenetic change. As holoblastic cleavage patterns were modified into meroblastic cleavage patterns, major changes in cell identity specification occurred within the mesendodermal marginal zone, as well as in the superficial epithelium of the embryo. These cellular identity changes resulted in the appearance of two novel extra-embryonic tissues within the embryos of teleostean fishes: the enveloping layer (EVL) and the yolk syncytial layer (YSL). The generation of these extra-embryonic tissues prompted major morphogenetic changes within the Organizer Region. As these evolutionary changes occurred, the outermost cell layer of the Organizer (the Organizer Epithelium) was apparently retained as a signaling center necessary for the establishment of left-right embryonic asymmetry in the embryo. Conserved and derived features of Organizer morphogenesis and gastrulation within ancient lineages of ray-finned fishes provide important insights into how the genetically encoded cell behaviors of early morphogenesis can be altered during the course of evolution. In particular, a highly divergent form of actinopterygian gastrulation, which is found in the annual fishes of South America, demonstrates that no aspect of vertebrate gastrulation is inherently immutable to evolutionary change.     

Symposium Summary: Evolutionary Patterns in Actinopterygian Fishes

SYNOPSIS. The actinopterygian fishes are an exemplary cladefor the study of structural and functional evolutionary patterns.With over half of all vertebrate species, ray-finned fisheshave diversified into a wide variety of habitats, and considerableprogress has been made over the last fifteen years in understandingthe genealogical relationships of actinopterygians. This symposiumhas contributed to our understanding of phylogenetic patternsin actinopterygians and to knowledge of the major structuraland functional patterns in locomotor, auditory, trophic, andneural systems. A number of key areas for future research havebeen identified. (1) The relationships of "palaeonisciform"fishes, (2) the study of trends in feeding and locomotor systemswithin a phylogenetic context, (3) the identification of primitivepatterns of pharyngeal jaw movement and steady and unsteadylocomotor patterns in actinopterygians, (4) the homologies,identification, and functional significance of neural pathwaysin the telencephalon, and (5) the comparative study of form-functionrelations in the auditory system. The study of teleost fishbiology has proceeded at the expense of data on primitive actinopterygians(e.g., Polypterus, Polyodon, Aapenser, Lepisosteus, Amia) whichare especially important in the analysis of structural and functionalpatterns in ray-finned fishes.     

The nucleotypic effects of cellular DNA content in cartilaginous and ray-finned fishes.

Cytological and organismal characteristics associated with cellular DNA content underpin most adaptionist interpretations of genome size variation. Since fishes are the only group of vertebrate for which relationships between genome size and key cellular parameters are uncertain, the cytological correlates of genome size were examined in this group. The cell and nuclear areas of erythrocytes showed a highly significant positive correlation with each other and with genome size across 22 cartilaginous and 201 ray-finned fishes. Regressions remained significant at all taxonomic levels, as well as among different fish lineages. However, the results revealed that cartilaginous fishes possess higher cytogenomic ratios than ray-finned fishes, as do cold-water fishes relative to their warm-water counterparts. Increases in genome size owing to ploidy shifts were found to influence cell and nucleus size in an immediate and causative manner, an effect that persists in ancient polyploid lineages. These correlations with cytological parameters known to have important influences on organismal phenotypes support an adaptive interpretation for genome size variation in fishes.     

Head size constrains forebrain development and evolution in ray-finned fishes

In ray-finned fishes, which comprise nearly half of all vertebrate species, the telencephalon does not evaginate, as it does in other vertebrates, but instead everts. No detailed explanation for this species difference has ever been offered. Here we propose that telencephalic eversion evolved because ray-finned fish embryos are so small that their telencephalon cannot evaginate but must, instead, squeeze into the space just dorsal to the developing nasal epithelia and rostral to the eyes-morphogenetic movements that amount to eversion. Evidence for this hypothesis derives from cladistic analyses, which show that early ray-finned fishes reduced their adult body size and adopted a novel reproductive strategy, based on the production of myriad minute young. Because body size tends to be inversely proportional to brain:body ratio, this phylogenetic reduction in body size implies that embryonic ray-finned fishes should have proportionately larger brains than embryos of species whose telencephalons evaginate. This prediction was confirmed by comparing serially sectioned heads of representative ray-finned and cartilaginous fish embryos at several stages of development. The brain, excluding its ventricles, occupies 36-46% of the cranial cavity in embryonic ray-finned fishes, but less than 20% in embryonic sharks. Moreover, three-dimensional reconstructions show that in embryonic ray-finned fishes the telencephalon has no room for a full-fledged evagination; instead, it spreads into the spaces just dorsal and caudal to the developing nasal epithelia. These morphogenetic movements, in conjunction with a thinning of the forebrain roof, generate telencephalic eversion.     

Strengths and limitations of molecular sequence comparisons for inferring the phylogeny of the major groups of fishes

Although the phylogenetic relationships of the major groups of fishes have been extensively studied with morphological characters, not all have been convincingly resolved. Analyses of molecular sequences from these groups may provide additional insights into problematical relationships, but are only just beginning to appear. We compare our own results from analyses of 18s ribosomal RNA sequences with those of other studies using globins, parvalbumins, insulin, 28s ribosomal RNA, and portions of two mitochondria1 genes (12S ribosomal RNA and cytochrome b ). Our evaluation of these studies reveals some of the difficulties encountered in reconstructing ancient divergences within the fishes, including unequal rates of evolution (among regions of a molecule as well as among lineages), gene duplication, extinction of lineages, and a possible rapid radiation of gnathostome higher taxa. The importance of evaluating the robustness of particular phylogenetic hypotheses is stressed. Some molecules appear to be inappropriate for investigating higher level divergences within the fishes; others are more promising, but must be examined in more taxa to allow an adequate evaluation of their utility. Convincing support for particular hypotheses of relationship will ultimately require congruence of trees generated from independent molecular data sets.     

Evolution of dimorphisms of the proteasome subunit beta type 8 gene (PSMB8) in basal ray-finned fish

The proteasome subunit beta type 8 (PSMB8) gene encodes a catalytic subunit of immunoproteasome that plays a central role in the processing of antigenic peptides presented by major histocompatibility complex class I molecules. The A- and F-type alleles defined by the 31st amino acid residue determining cleaving specificity have been identified from ray-finned fish, amphibia, and reptiles. These two types show extremely long-term trans-species polymorphism in Polypteriformes, Cypriniformes, and Salmoniformes, suggesting the presence of very ancient lineages termed A and F. To elucidate the evolution of the PSMB8 dimorphism in basal ray-finned fish, we analyzed Pantodon buchholzi (Osteoglossiformes), seven species of Anguilliformes, and Hypomesus nipponensis (Osmeriformes). Both A and F lineage sequences were identified from P. buchholzi and H. nipponensis, confirming that these two lineages have been conserved by basal ray-finned fish. However, both the A- and F-type alleles found in Anguilliformes species belonged to the F lineage irrespective of their types. This apparently suggests that the A lineage was lost in the common ancestor of Anguilliformes, and recovery of the A type within the F lineage occurred in Anguilliformes. The apparent loss of the F lineage and recovery of the F type within the A lineage have already been reported from tetrapods and higher teleosts. However, this is the first report on the reverse situation and reveals the dynamic evolution of the PSMB8 dimorphism.     

Shallow water ray-finned marine fishes follow Bergmann’s rule

Understanding the interspecific variation in body size across macroclimatic gradients has been of paramount importance to naturalists and biogeographers. Bergmann’s rule, which describes a trend of increasing body size polewards, is arguably the best-known ecogeographical rule in terrestrial environments but remains largely unexplored in the marine realm. In this study we tested Bergmann's rule in marine ray-finned fishes (Pisces, Actinopterygii), analyzing the relationship between body size and latitude in 5662 species. To examine possible underlying mechanisms, we adopted a cross-species approach to evaluate the association of body size with four predictors: Sea Surface Temperature, Net Primary Productivity, Salinity, and Human impact. We analyzed the relationships between body size and environmental and anthropogenic variables building mixed linear models, which considered the taxonomic structure in the data. We conducted complementary analyses dividing the data into five latitudinal bands. Actinopterygii showed a clear Bergmannian pattern, with the largest species observed in temperate regions, being the first global analysis on ray-finned fishes showing a pattern consistent with Bergmann’s rule. Sea Surface Temperature and Net Primary Productivity were the best predictors, in accordance with the time to sexual maturity and resource availability hypotheses. Our analyses based on latitudinal bands showed a differential response of body size to the environment, with temperature, salinity and human impact more strongly associated with size variation at cold environments. These results agree with previous studies on Bergmann’s rule for terrestrial ectothermic, freshwater and marine fishes. Our findings suggest that temperature rise in the ocean and growing human impact may have effects on the distribution of body size, thus altering ecosystem functioning. Fundamental differences often assumed to exist between marine and terrestrial systems are not so evidently reflected in the emergence of large-scale body size gradients.     

Opsin gene duplication and divergence in ray-finned fish

Opsin gene sequences were first reported in the 1980s. The goal of that research was to test the hypothesis that human opsins were members of a single gene family and that variation in human color vision was mediated by mutations in these genes. While the new data supported both hypotheses, the greatest contribution of this work was, arguably, that it provided the data necessary for PCR-based surveys in a diversity of other species. Such studies, and recent whole genome sequencing projects, have uncovered exceptionally large opsin gene repertoires in ray-finned fishes (taxon, Actinopterygii). Guppies and zebrafish, for example, have 10 visual opsin genes each. Here we review the duplication and divergence events that have generated these gene collections. Phylogenetic analyses revealed that large opsin gene repertories in fish have been generated by gene duplication and divergence events that span the age of the ray-finned fishes. Data from whole genome sequencing projects and from large-insert clones show that tandem duplication is the primary mode of opsin gene family expansion in fishes. In some instances gene conversion between tandem duplicates has obscured evolutionary relationships among genes and generated unique key-site haplotypes. We mapped amino acid substitutions at so-called key-sites onto phylogenies and this exposed many examples of convergence. We found that dN/dS values were higher on the branches of our trees that followed gene duplication than on branches that followed speciation events, suggesting that duplication relaxes constraints on opsin sequence evolution. Though the focus of the review is opsin sequence evolution, we also note that there are few clear connections between opsin gene repertoires and variation in spectral environment, morphological traits, or life history traits.