Evolution and phylogenetic utility of alignment gaps within intron sequences of three nuclear genes in bumble bees (Bombus)

Dissecting the ancient rapid radiation of microgastrine wasp genera using additional nuclear genes

Previous estimates of a generic level phylogeny for the ubiquitous parasitoid wasp subfamily Microgastrinae (Hymenoptera) have been problematic due to short internal branches deep in the phylogeny. These short branches might be attributed to a rapid radiation among the taxa, the use of genes that are unsuitable for the levels of divergence being examined, or insufficient quantity of data. We added over 1200 nucleotides from four nuclear genes to a dataset derived from three genes to produce a dataset of over 3000 nucleotides per taxon. While the number of well-supported short branches in the phylogeny increased, we still did not obtain strong bootstrap support for every node. Parametric and nonparametric bootstrap simulations projected that an enormous, and likely unobtainable, amount of data would be required to get bootstrap support greater than 50% for every node. However, a marked increase in the number of well-supported nodes was seen when we conducted a Bayesian analysis of a combined dataset generated from morphological characters added to the seven gene dataset. Our results suggest that, in some cases, combining morphological and genetic characters may be the most practical way to increase support for short branches deep in a phylogeny.

Analysis of family-level relationships in bees (Hymenoptera: Apiformes) using 28S and two previously unexplored nuclear genes: CAD and RNA polymerase II

We analyzed a combined data set of two protein-coding nuclear genes (CAD and RNA polymerase II) and a nuclear ribosomal gene (28S D2-D4 region) for 68 bee species and 11 wasp outgroups. Our taxon sampling included all seven extant bee families, 17 of 20 subfamilies, and diverse tribes. Wasp outgroups included the two families most closely related to bees: Crabronidae and Sphecidae. We analyzed the combined and single gene data sets using parsimony and Bayesian methods, which yielded largely congruent results. Our results provide reasonably strong support for family and subfamily-level relationships among bees. Our data set strongly supports the sister-group relationship of the Colletidae and Stenotritidae, and places Halictidae as sister to this clade combined. Our analyses place the Melittidae and the long-tongued (LT) bee clade (Apidae+Megachilidae) near the base of the tree with Colletidae (and Stenotritidae) in a fairly highly derived position. This topology ("Melittidae-LT basal") was obtained in previous morphological studies under certain methods of character coding. A more widely accepted tree topology that places Colletidae (and/or Stenotritidae) as sister to all other bees ("Colletidae basal") is not supported by our data. The "Melittidae-LT basal" hypothesis may better explain patterns in the bee fossil record as well as historical biogeography of certain bee groups. Our results provide new insights into higher-level bee phylogeny and indicate that CAD, RNA polymerase II, and 28S are useful data sets for resolving Cretaceous-age divergences in bees and other Hymenoptera.

Multiple gene genealogies and AFLPs suggest cryptic speciation and long-distance dispersal in the basidiomycete Serpula himantioides (Boletales)

Serpula himantioides (Boletales, Basidiomycota) produces thin resupinate basidiocarps on dead coniferous wood worldwide and causes damage in buildings as well. In this study, we present evidence for the existence of at least three phylogenetically defined cryptic species (referred to as Sib I-III) within the morphospecies S. himantioides, a conclusion based on analyses of sequence data from four DNA regions and amplified fragment length polymorphisms (AFLPS). A low degree of shared sequence polymorphisms was observed among the three lineages indicating a long-lasting separation. The AFLPs revealed two additional subgroups within Sib III. Results from mating studies were consistent with the molecular data. In Sib III, no correspondence between genetic and geographical distance was observed among isolates worldwide, presumably reflecting recent dispersal events. Our results indicate that at least two of the lineages (Sib II and Sib III) have wide sympatric distributions. A population genetic analysis of Sib III isolates, scoring sequence polymorphisms as codominant SNP markers, indicates that panmictic conditions exist in the Sib III group. This study supports the view that cryptic speciation is a common phenomenon in basidiomycete fungi and that phylogenetic species recognition can be a powerful inference to detect cryptic species. Furthermore, this study shows that AFLP data are a valuable supplement to DNA sequence data in that they may detect a finer level of genetic variation.

Incorporating information from length-mutational events into phylogenetic analysis

With the growing number of phylogenetic studies that use length variable DNA sequences, incorporating information from length-mutational events into phylogenetic analysis is becoming increasingly important. A new method, modified complex indel coding is described that aims at maximizing the phylogenetic information retained from unambiguously aligned sequence regions or regions where the principal relative position of gaps to one another can be safely established. An algorithm is described that allows application of the method to all theoretically possible gap-nucleotide patterns. A platform-independent computer program is introduced that automates the new method as well as several previously published coding schemes. Differences to previously published indel coding approaches as well as to the integration of ambiguously aligned regions into phylogenetic analysis are discussed.

Optimal Intron Analyses in the Trimeresurus Radiation of Asian Pitvipers

Nuclear introns are commonly used as phylogenetic markers, but a number of issues related to alignment strategies, indel treatments, and the incorporation of length-variant heterozygotes (LVHs) are not routinely addressed when generating phylogenetic hypotheses. Topological congruence in relation to an extensive mitochondrial DNA multigene phylogeny (derived from 2,423 bp of 12S, 16S, ND4, and CYTB genes) of the Asian pitviper Trimeresurus radiation was used to compare combinations of "by eye" and edited and unedited ClustalX 1.8 alignments of two nuclear introns. Indels were treated as missing data, fifth character states, and assigned simple and multistate codes. Upon recovery of the optimal alignment and indel treatment strategy, a total evidence approach was used to investigate the phylogenetic utility of the indels and test new generic arrangements within Trimeresurus. Approximately one third of the intron data partitions exhibited LVHs, suggesting that they are common in introns. Furthermore, a simple concatenation approach can facilitate the incorporation of LVHs into phylogenetic analyses to make use of all available data and investigate mechanisms of molecular evolution. Analyses of ClustalX 1.8-assisted alignments were generally more congruent than the "by eye" alignment and the analysis of a simple coded, edited ClustalX 1.8 (gap opening cost 5, gap extension cost 1) alignment revealed the most congruent tree. The total evidence approach supported the new arrangements within Trimeresurus, suggesting that the phylogeny should be considered as a working benchmark in Asian pitviper systematics. Finally, a critical appraisal of the diverse array of indels (56 to 57 per intron, ranging from 1 to 151 bp in length) suggested that they are a combination of Hennigian and homoplasious events unrelated to indel size or location within the intron. [Alignment; indels; intron analysis; length-variant heterozygotes; Trimeresurus.].

Molecular phylogeny of the bumble bee subgenus Pyrobombus (Hymenoptera:Apidae:Bombus) with insights into gene utility for lower-level analysis

Comprising nearly 20% of all bumble bees, the subgenus Pyrobombus is distributed across diverse habitats in the Northern Hemisphere and exhibits considerable morphological and behavioural variation relative to other subgenera. Its size and variation have led to questions concerning its monophyly and intrasubgeneric relationships, but too few known morphological synapomorphies and insufficient taxon sampling have precluded robust answers to these questions. To obtain a robust phylogeny of the group, we obtained DNA sequences for 36 of the 43 species from four genes (mitochondrial 16S rRNA and three nuclear genes: elongation factor – 1α (EF-1α), long wavelength rhodopsin (LW Rh or opsin) and arginine kinase (ArgK)). Both Bayesian and parsimony phylogenies are well resolved and indicate a monophyletic Pyrobombus when assessed against representatives of 20 additional subgenera. The more conserved nuclear genes, especially EF-1α and ArgK, provided good support across all of the taxonomic levels examined, whereas support of the more rapidly evolving mt16S was restricted mostly to close relationships at the tips of the tree. The exon regions of ArgK were the most conserved and may be promising for higher-level phylogenetics. We discuss species relationships within Pyrobombus and its sister-group, Bombus s.s. + Alpinobombus, in relation to previous taxonomic studies.

The evolution of the spindlin gene in birds: sequence analysis of an intron of the spindlin W and Z gene reveals four major divisions of the Psittaciformes

The Psittaciformes (parrots, parakeets) are among the most widely held captive birds. Yet, their evolution and their phylogenetic relationships have been relatively little studied. This paper describes the phylogenetic relationships between a number of Psittaciformes as derived from the sequences of the third intron of the Z-chromosomal and W-chromosomal spindlin genes. The Z-chromosomal sequences of the kakapo (Strigops habroptilus), the kea (Nestor notabilis), and the kaka (Nestor meridionalis) from New Zealand form a cluster which is the sister group to all other Psittaciformes. The results show further that the Z-chromosomal sequences of the other species can be divided into two groups based on the occurrence of a sequence element ACCCT. The group with the insert (A) is mainly from species with an Australasian geographical distribution and includes such species as the Lories (Lorius, etc.), the budgerigar (Melospittacus undulatus), and the rosellas (Platycercus). It also includes the African lovebirds (Agapornidae), which are the only representative of group A outside Australasia. Group B, without the insert, includes the neotropical parrots and parakeets such as the amazons (Amazona, etc.), the macaws (Ara, etc.), and the conures (Aratinga, etc.), the Australian Cacatuini and the African species such as the African grey parrot (Psittacus erithacus) as well as Coracopsis vasa from Madagascar and Psittrichas fulgidus from New Guinea. The W-chromosomal sequence data show that another division of the Psittacidae is found in the replacement of a pyrimidine-rich segment occurring in many non-psittacines as well as the kakapo (S. habroptilus), the kea (N. notabilis), the kaka (N. meridionalis), and the Cacatuini by a microsatellite consisting of a variable number of TATTA monomers in the other Psittaciformes. The results support a Gondwanan origin of the Psittaciformes and the suggestion that paleogeographic events were a major force in psittacine divergence.

Molecular characterization and expression of the UV opsin in bumblebees: three ommatidial subtypes in the retina and a new photoreceptor organ in the lamina

Ultraviolet-sensitive photoreceptors have been shown to be important for a variety of visual tasks performed by bees, such as orientation, color and polarization vision, yet little is known about their spatial distribution in the compound eye or optic lobe. We cloned and sequenced a UV opsin mRNA transcript from Bombus impatiens head-specific cDNA and, using western blot analysis, detected an eye protein band of approximately 41 kDa, corresponding to the predicted molecular mass of the encoded opsin. We then characterized UV opsin expression in the retina, ocelli and brain using immunocytochemistry. In the main retina, we found three different ommatidial types with respect to the number of UV opsin-expressing photoreceptor cells, namely ommatidia containing two, one or no UV opsin-immunoreactive cells. We also observed UV opsin expression in the ocelli. These results indicate that the cloned opsin probably encodes the P350 nm pigment, which was previously characterized by physiological recordings. Surprisingly, in addition to expression in the retina and ocelli, we found opsin expression in different parts of the brain. UV opsin immunoreactivity was detected in the proximal rim of the lamina adjacent to the first optic chiasm, which is where studies in other insects have found expression of proteins involved in the circadian clock, period and cryptochrome. We also found UV opsin immunoreactivity in the core region of the antennal lobe glomeruli and different clusters of perikarya within the protocerebrum, indicating a putative function of these brain regions, together with the lamina organ, in the entrainment of circadian rhythms. In order to test for a possible overlap of clock protein and UV opsin spatial expression, we also examined the expression of the period protein in these regions.

A multi-gene phylogeny of aquiline eagles (Aves: Accipitriformes) reveals extensive paraphyly at the genus level

The phylogeny of the tribe Aquilini (eagles with fully feathered tarsi) was investigated using 4.2 kb of DNA sequence of one mitochondrial (cyt b) and three nuclear loci (RAG-1 coding region, LDH intron 3, and adenylate-kinase intron 5). Phylogenetic signal was highly congruent and complementary between mtDNA and nuclear genes. In addition to single-nucleotide variation, shared deletions in nuclear introns supported one basal and two peripheral clades within the Aquilini. Monophyly of the Aquilini relative to other birds of prey was confirmed. However, all polytypic genera within the tribe, Spizaetus, Aquila, Hieraaetus, turned out to be non-monophyletic. Old World Spizaetus and Stephanoaetus together appear to be the sister group of the rest of the Aquilini. Spizastur melanoleucus and Oroaetus isidori are nested among the New World Spizaetus species and should be merged with that genus. The Old World 'Spizaetus' species should be assigned to the genus Nisaetus (Hodgson, 1836). The sister species of the two spotted eagles (Aquila clanga and Aquila pomarina) is the African Long-crested Eagle (Lophaetus occipitalis). Hieraaetus fasciatus/spilogaster are closest to Aquila verreauxii and should be merged with that genus. Wahlberg's Eagle H. wahlbergi, formerly placed in Aquila, is part of a clade including three small Hieraaetus species (pennatus, ayresii, and morphnoides). The Martial Eagle (Polemaetus bellicosus) is the sister species of the Aquila/Hieraaetus/Lophaetus clade. Basal relationships within this clade remained unresolved. Parsimony reconstruction of the evolution of plumage pattern within Aquilini suggests that: (1) transverse barring of parts of the body plumage was lost in the Palearctic Aquila-Hieraaetus clade, (2) pale underparts in adult plumage evolved three times independently, and (3) dimorphic adult plumage is a derived character of the small-bodied Hieraaetus clade.

Using exon and intron sequences of the gene Mp20 to resolve basal relationships in Cicindela (Coleoptera:Cicindelidae)

The genus Cicindela (Coleoptera: Cicindelidae) is a species-rich cosmopolitan group of tiger beetles useful for comparing clade diversification worldwide. Knowledge about relationships of major groups is important for this analysis but basal nodes in Cicindela have been difficult to resolve with standard mtDNA markers. Here we developed the Mp20 gene, a single-copy nuclear marker coding for a muscle-associated protein in insects, for phylogenetic analysis of basal groups of Cicindela. Nearly full-length sequences were obtained for 51 cicindelids, including major taxonomic groups from all continents. Sequences of Mp20 were between 1.2 and 1.7 kb and spanning three introns. Phylogenetic signal of exon and intron sequences was compared with that from four gene regions of mtDNA (COI, COIII, Cytb, 16S rRNA; 2.4 kb total). Because introns differed in length, sequence alignment was conducted using various procedures of phenetic and parsimony-based character coding of indels to assess their phylogenetic information content, but major nodes were recovered consistently. Mp20 sequences contributed two thirds of the total support of the combined analysis, with most signal from the introns. We found major clades of Cicindela to be geographically largely coincident with continental regions, confined to Australasia, the Holarctic, the Indian subcontinent, Africa, and South and Central America. Clock estimates using various maximum-likelihood (ML) branch length calculations resulted in roughly similar divergence times whether Mp20 exon, introns, or mtDNA were used, and they were not greatly affected by different procedures for coding and optimizing indel characters. Based on existing clock calibrations in Cicindela, basal splits of continental lineages occurred in the mid-Miocene, placing the radiation of basal groups of Cicindela to a period when their open-vegetation habitats expanded globally.

Molecular phylogenetics suggest a North American link between the anthropogenic dry rot fungus Serpula lacrymans and its wild relative S. himantioides

The dry rot fungus Serpula lacrymans causes damages in wooden buildings and constructions in temperate regions worldwide. In this study, the global phylogeography of S. lacrymans and its wild relative S. himantioides has been investigated to clarify genealogical relationships and determine the origin and spread of the building strains. Internal transcribed spacer (ITS) nrDNA and parts of the beta-tubulin (tub) and the translation elongation factor (efa) 1a genes were sequenced, and phylogenetic relationships inferred. Some analyses suggest that S. lacrymans may have originated from an ancient S. himantioides lineage, but most results support that S. lacrymans and S. himantioides are monophyletic sister species. Phylogenetic analysis of the ITS data revealed two subgroups within S. lacrymans corresponding to two earlier described varieties; one group occurring frequently in houses worldwide ('Domesticus'), and one group represented by individuals from forests in Northern California ('Shastensis'). A few collections from nature were included in the Domesticus group as well, among other specimens from two newly discovered localities in Far East Russia and Siberia. In the Domesticus group little sequence variation occurs, suggesting a recent worldwide dispersal, possibly linked to human activity. Phylogenetic analyses indicate that the Domesticus group may have originated from an ancient lineage related closely to the Shastensis group. A remarkable shift in morphology and habitat preferences has occurred during the evolution of the Domesticus lineage, linked to the transition from nature to human-made habitats.

Nonuniform Concerted Evolution and Chloroplast Capture: Heterogeneity of Observed Introgression Patterns in Three Molecular Data Partition Phylogenies of Asian Mitella (Saxifragaceae)

Interspecific hybridization is one of the major factors leading to phylogenetic incongruence among loci, but the knowledge is still limited about the potential of each locus to introgress between species. By directly sequencing three DNA regions: chloroplast DNAs (matK gene and trnL-F noncoding region), the nuclear ribosomal external transcribed spacer (ETS) region, and internal transcribed spacer (ITS) regions, we construct three phylogenetic trees of Asian species of Mitella (Saxifragaceae), a genus of perennials in which natural hybrids are commonly observed. Within this genus, there is a significant topological conflict between chloroplast and nuclear phylogenies and also between the ETS and the ITS, which can be attributed to frequent hybridization within the lineage. Chloroplast DNAs show the most extensive introgression pattern, ITS regions show a moderate pattern, and the ETS region shows no evidence of introgression. Nonuniform concerted evolution best explains the difference in the introgression patterns between the ETS region and ITS regions, as the sequence heterogeneity of the ITS region within an individual genome is estimated to be twice that of an ETS in this lineage. Significant gene conversion patterns between two hybridizing taxa were observed in contiguous arrays of cloned ETS-ITS sequences, further confirming that only ITS regions have introgressed bidirectionally. The relatively slow concerted evolution in the ITS regions probably allows the coexistence of multiple alleles within a genome, whereas the strong concerted evolution in the ETS region rapidly eliminates heterogeneous alleles derived from other species, resulting in species delimitations highly concordant with those based on morphology. This finding indicates that the use of multiple molecular tools has the potential to reveal detailed organismal evolution processes involving interspecific hybridization, as an individual locus varies greatly in its potential to introgress between species.

Single-copy nuclear genes recover cretaceous-age divergences in bees

We analyzed the higher level phylogeny of the bee family Halictidae based on the coding regions of three single-copy nuclear genes (long-wavelength [LW] opsin, wingless, and elongation factor 1-alpha [EF-1 alpha]). Our combined data set consisted of 2,234 aligned nucleotide sites (702 base pairs [bp] for LW opsin, 405 bp for wingless, and 1,127 bp for EF-1 alpha) and 779 parsimony-informative sites. We included 58 species of halictid bees from 33 genera, representing all subfamilies and tribes, and rooted the trees using seven outgroups from other bee families: Colletidae, Andrenidae, Melittidae, and Apidae. We analyzed the separate and combined data sets by a variety of methods, including equal weights parsimony, maximum likelihood, and Bayesian methods. Analysis of the combined data set produced a strong phylogenetic signal with high bootstrap and Bremer support and high posterior probability well into the base of the tree. The phylogeny recovered the monophyly of the Halictidae and of all four subfamilies and both tribes, recovered relationships among the subfamilies and tribes congruent with morphology, and provided robust support for the relationships among the numerous genera in the tribe Halictini, sensu Michener (2000). Using our combined nucleotide data set, several recently described halictid fossils from the Oligocene and Eocene, and recently developed Bayesian methods, we estimated the antiquity of major clades within the family. Our results indicate that each of the four subfamilies arose well before the Cretaceous-Tertiary boundary and suggest that the early radiation of halictid bees involved substantial African-South American interchange roughly coincident with the separation of these two continents in the late Cretaceous. This combination of single-copy nuclear genes is capable of recovering Cretaceous-age divergences in bees with high levels of support. We propose that LW opsin, wingless, and EF-1 alpha(F2 copy) may be useful in resolving relationships among bee families and other Cretaceous-age insect lineages.

Phylogeny of the orchid bees (Hymenoptera: Apinae: Euglossini): DNA and morphology yield equivalent patterns

Orchid bees (Euglossini) are spectacular long-tongued Neotropical bees important in the pollination of Neotropical long-corolla flowers, particularly some orchids. Besides remarkably long tongues, males in particular exhibit other flower-related adaptations, including setal brushes on the foretarsi used for rasping the petals of orchids while collecting aromatic compounds. These compounds are stored in large swollen tibiae and are thought to play an important role in courtship behavior. Euglossini are also unusual in lacking sociality; they are the only tribe among the corbiculate bees that are not eusocial, and two of the genera are cleptoparasitic. Each genus exhibits distinct behavioral traits including nest architecture and host-parasite interactions, yet their evolution is unknown. Despite previous phylogenetic studies of on morphological characters, the relationships among the five euglossine genera remain under debate. We investigate euglossine generic relationships using DNA sequence data from four genes and new morphological characters. The morphological and molecular data yield congruent evolutionary patterns, and combining the data gives a fully resolved and well supported phylogeny of Euglossini.

Phylogeny of a Macaronesian radiation: host-plant use and possible cryptic speciation in Liparthrum bark beetles

The Macaronesian islands are well known for their unique endemic floras of woody plants. Many of these unusual plant groups provide important novel resources for bark and wood boring beetles which breed in dead or moribund parts of their host plants. The bark beetle genus Liparthrum exploits a wide range of unusual host plants and has its highest proportion of species living on the Macaronesian Islands. We used DNA sequences of the mitochondrial Cytochrome Oxidase I gene and the nuclear Elongation Factor 1 alpha gene, and morphological characters, to estimate the phylogenetic relationships among species endemic to these archipelagos, and to trace the evolution of host-plant use. All parsimony and Bayesian analyses of the combined data, and maximum likelihood analyses of the molecular data, showed that species associated with Euphorbia are monophyletic. We also found genetic and subtle morphological evidence for three cases of cryptic speciation in one polyphyletic species associated with different Euphorbia plants, showing that high levels of host specialisation can occur also in insects breeding in older and very dry, dead plant tissues.

How do insect nuclear and mitochondrial gene substitution patterns differ? Insights from Bayesian analyses of combined datasets

We analyzed 12 combined mitochondrial and nuclear gene datasets in seven orders of insects using both equal weights parsimony (to evaluate phylogenetic utility) and Bayesian methods (to investigate substitution patterns). For the Bayesian analyses we used relatively complex models (e.g., general time reversible models with rate variation) that allowed us to quantitatively compare relative rates among genes and codon positions, patterns of rate variation among genes, and substitution patterns within genes. Our analyses indicate that nuclear and mitochondrial genes differ in a number of important ways, some of which are correlated with phylogenetic utility. First and most obviously, nuclear genes generally evolve more slowly than mitochondrial genes (except in one case), making them better markers for deep divergences. Second, nuclear genes showed universally high values of CI and (generally) contribute more to overall tree resolution than mitochondrial genes (as measured by partitioned Bremer support). Third, nuclear genes show more homogeneous patterns of among-site rate variation (higher values of alpha than mitochondrial genes). Finally, nuclear genes show more symmetrical transformation rate matrices than mitochondrial genes. The combination of low values of alpha and highly asymmetrical transformation rate matrices may explain the overall poor performance of mitochondrial genes when compared to nuclear genes in the same analysis. Our analyses indicate that some parameters are highly correlated. For example, A/T bias was positively and significantly associated with relative rate and CI was positively and significantly associated with alpha (the shape of the gamma distribution). These results provide important insights into the substitution patterns that might characterized high quality genes for phylogenetic analysis: high values of alpha, unbiased base composition, and symmetrical transformation rate matrices. We argue that insect molecular systematists should increasingly focus on nuclear rather than mitochondrial gene datasets because nuclear genes do not suffer from the same substitutional biases that characterize mitochondrial genes.