Recent Intron Gain in Elongation Factor-1α of Colletid Bees (Hymenoptera: Colletidae)

The evolutionary history of the cellophane bee genus Colletes Latreille (Hymenoptera: Colletidae): Molecular phylogeny, biogeography and implications for a global infrageneric classification

Colletes Latreille (Hymenoptera: Colletidae) is a diverse genus with 518 valid species distributed in all biogeographic realms, except Australasia and Antarctica. Here we provide a comprehensive dated phylogeny for Colletes based on Bayesian and maximum likelihood-based analyses of DNA sequence data of six loci: 28S rDNA, cytochrome c oxidase subunit 1, elongation factor-1α copy F2, long-wavelength rhodopsin, RNA polymerase II and wingless. In total, our multilocus matrix consists of 4824 aligned base pairs for 143 species, including 112 Colletes species plus 31 outgroups (one stenotritid and a diverse array of colletids representing all subfamilies). Overall, analyses of each of the six single-locus datasets resulted in poorly resolved consensus trees with conflicting phylogenetic signal. However, our analyses of the multilocus matrix provided strong support for the monophyly of Colletes and show that it can be subdivided into five major clades. The implications of our phylogenetic results for future attempts at infrageneric classification for the Colletes of the world are discussed. We propose species groups for the Neotropical species of Colletes, the only major biogeographic realm for which no species groups have been proposed to date. Our dating analysis indicated that Colletes diverged from its sister taxon, Hemicotelles Toro and Cabezas, in the early Oligocene and that its extant lineages began diversifying only in the late Oligocene. According to our biogeographic reconstruction, Colletes originated in the Neotropics (most likely within South America) and then spread to the Nearctic very early in its evolutionary history. Geodispersal to the Old World occurred soon after colonization of the Northern Hemisphere. Lastly, the historical biogeography of Colletes is analyzed in light of available geological and palaeoenvironmental data.

The roles of barriers, refugia, and chromosomal clines underlying diversification in Atlantic Forest social wasps

Phylogeographic studies have sought to explain the genetic imprints of historical climatic changes and geographic barriers within the Brazilian Atlantic Forest (AF) biota, and consequently two processes of diversification (refugia and barriers) have been proposed. Additionally, there is evidence that eustatic changes influenced the biogeographic history of the AF. Here we evaluate these contrasting diversification processes using two AF social wasp species – the mid-montane Synoeca cyanea and the lowland Synoeca aff. septentrionalis. We analyzed several sources of data including multilocus DNA sequence, climatic niche models and chromosomal features. We find support for idiosyncratic phylogeographic patterns between these wasps, involving different levels of population structure and genetic diversity, contrary suitable climatic conditions during the last glaciation, and contrasting historical movements along the AF. Our data indicate that neotectonics and refugia played distinct roles in shaping the genetic structure of these wasps. However, we argue that eustatic changes influenced the demographic expansion but not population structure in AF biota. Notably, these wasps exhibited chromosomal clines, involving chromosome number and decreasing of GC content, latitudinally oriented along the AF. Together, these results reinforce the need to consider individual organismal histories and indicate that barriers and refugia are significant factors in understanding AF evolution.

Reverse transcriptase and intron number evolution

BACKGROUND

Introns are universal in eukaryotic genomes and play important roles in transcriptional regulation, mRNA export to the cytoplasm, nonsense-mediated decay as both a regulatory and a splicing quality control mechanism, R-loop avoidance, alternative splicing, chromatin structure, and evolution by exon-shuffling.

METHODS

Sixteen complete fungal genomes were used 13 of which were sequenced and annotated by JGI. Ustilago maydis, Cryptococcus neoformans, and Coprinus cinereus (also named Coprinopsis cinerea) were from the Broad Institute. Gene models from JGI-annotated genomes were taken from the GeneCatalog track that contained the best representative gene models. Varying fractions of the GeneCatalog were manually curated by external users. For clarity, we used the JGI unique database identifier.

RESULTS

The last common ancestor of eukaryotes (LECA) has an estimated 6.4 coding exons per gene (EPG) and evolved into the diverse eukaryotic life forms, which is recapitulated by the development of a stem cell. We found a parallel between the simulated reverse transcriptase (RT)-mediated intron loss and the comparative analysis of 16 fungal genomes that spanned a wide range of intron density. Although footprints of RT (RTF) were dynamic, relative intron location (RIL) to the 5'-end of mRNA faithfully traced RT-mediated intron loss and revealed 7.7 EPG for LECA. The mode of exon length distribution was conserved in simulated intron loss, which was exemplified by the shared mode of 75 nt between fungal and Chlamydomonas genomes. The dominant ancient exon length was corroborated by the average exon length of the most intron-rich genes in fungal genomes and consistent with ancient protein modules being ~25 aa. Combined with the conservation of a protein length of 400 aa, the earliest ancestor of eukaryotes could have 16 EPG. During earlier evolution, Ascomycota's ancestor had significantly more 3'-biased RT-mediated intron loss that was followed by dramatic RTF loss. There was a down trend of EPG from more conserved to less conserved genes. Moreover, species-specific genes have higher exon-densities, shorter exons, and longer introns when compared to genes conserved at the phylum level. However, intron length in species-specific genes became shorter than that of genes conserved in all species after genomes experiencing drastic intron loss. The estimated EPG from the most frequent exon length is more than double that from the RIL method.

CONCLUSIONS

This implies significant intron loss during the very early period of eukaryotic evolution. De novo gene-birth contributes to shorter exons, longer introns, and higher exon-density in species-specific genes relative to conserved genes.

Diversification of Ponto-Caspian aquatic fauna: morphology and molecules retrieve congruent evolutionary relationships in Pontogammarus maeoticus (Amphipoda: Pontogammaridae)

The geological history of the Ponto-Caspian region, with alternating cycles of isolation and reconnection among the three main basins (Black and Azov Seas, and the more distant Caspian Sea) as well as between them and the Mediterranean Sea, profoundly affected the diversification of its aquatic fauna, leading to a high degree of endemism. Two alternative hypotheses on the origin of this amazing biodiversity have been proposed, corresponding to phases of allopatric separation of aquatic fauna among sea basins: a Late Miocene origin (10-6 MYA) vs. a more recent Pleistocene ancestry (<2 MYA). Both hypotheses support a vicariant origin of (1) Black+Azov Sea lineages on the one hand, and (2) Caspian Sea lineages on the other. Here, we present a study on the Ponto-Caspian endemic amphipod Pontogammarus maeoticus. We assessed patterns of divergence based on (a) two mitochondrial and one nuclear gene, and (b) a morphometric analysis of 23 morphological traits in 16 populations from South and West Caspian Sea, South Azov Sea and North-West Black Sea. Genetic data indicate a long and independent evolutionary history, dating back from the late Miocene to early Pleistocene (6.6-1.6 MYA), for an unexpected, major split between (i) a Black Sea clade and (ii) a well-supported clade grouping individuals from the Caspian and Azov Seas. Absence of shared haplotypes argues against either recent or human-mediated exchanges between Caspian and Azov Seas. A mismatch distribution analysis supports more stable population demography in the Caspian than in the Black Sea populations. Morphological divergence largely followed patterns of genetic divergence: our analyses grouped samples according to the basin of origin and corroborated the close phylogenetic affinity between Caspian and Azov Sea lineages. Altogether, our results highlight the necessity of careful (group-specific) evaluation of evolutionary trajectories in marine taxa that should certainly not be inferred from the current geographical proximity of sea basins alone.

Convergent intron gains in hymenopteran elongation factor-1α

The eukaryotic translation elongation factor-1α gene (eEF1A) has been used extensively in higher level phylogenetics of insects and other groups, despite being present in two or more copies in several taxa. Orthology assessment has relied heavily on the position of introns, but the basic assumption of low rates of intron loss and absence of convergent intron gains has not been tested thoroughly. Here, we study the evolution of eEF1A based on a broad sample of taxa in the insect order Hymenoptera. The gene is universally present in two copies - F1 and F2 - both of which apparently originated before the emergence of the order. An elevated ratio of non-synonymous versus synonymous substitutions and differences in rates of amino acid replacements between the copies suggest that they evolve independently, and phylogenetic methods clearly cluster the copies separately. The F2 copy appears to be ancient; it is orthologous with the copy known as F1 in Diptera, and is likely present in most insect orders. The hymenopteran F1 copy, which may or may not be unique to this order, apparently originated through retroposition and was originally intron free. During the evolution of the Hymenoptera, it has successively accumulated introns, at least three of which have appeared at the same position as introns in the F2 copy or in eEF1A copies in other insects. The sites of convergent intron gain are characterized by highly conserved nucleotides that strongly resemble specific intron-associated sequence motifs, so-called proto-splice sites. The significant rate of convergent intron gain renders intron-exon structure unreliable as an indicator of orthology in eEF1A, and probably also in other protein-coding genes.

The impact of molecular data on our understanding of bee phylogeny and evolution

Our understanding of bee phylogeny has improved over the past fifteen years as a result of new data, primarily nucleotide sequence data, and new methods, primarily model-based methods of phylogeny reconstruction. Phylogenetic studies based on single or, more commonly, multilocus data sets have helped resolve the placement of bees within the superfamily Apoidea; the relationships among the seven families of bees; and the relationships among bee subfamilies, tribes, genera, and species. In addition, molecular phylogenies have played an important role in inferring evolutionary patterns and processes in bees. Phylogenies have provided the comparative framework for understanding the evolution of host-plant associations and pollen specialization, the evolution of social behavior, and the evolution of parasitism. In this paper, we present an overview of significant discoveries in bee phylogeny based primarily on the application of molecular data. We review the phylogenetic hypotheses family-by-family and then describe how the new phylogenetic insights have altered our understanding of bee biology.

Diversity of phylogenetic information according to the locus and the taxonomic level: an example from a parasitic mesostigmatid mite genus

Molecular markers for cladistic analyses may perform differently according to the taxonomic group considered and the historical level under investigation. Here we evaluate the phylogenetic potential of five different markers for resolving evolutionary relationships within the ectoparasitic genus Dermanyssus at the species level, and their ability to address questions about the evolution of specialization. COI provided 9–18% divergence between species (up to 9% within species), 16S rRNA 10–16% (up to 4% within species), ITS1 and 2 2–9% (up to 1% within species) and Tropomyosin intron n 8–20% (up to 6% within species). EF-1α revealed different non-orthologous copies within individuals of Dermanyssus and Ornithonyssus. Tropomyosin intron n was shown containing consistent phylogenetic signal at the specific level within Dermanyssus and represents a promising marker for future prospects in phylogenetics of Acari. Phylogenetic analyses revealed that the generalist condition is apomorphic and D. gallinae might represent a complex of hybridized lineages. The split into hirsutus-group and gallinae-group in Dermanyssus does not seem to be appropriate based upon these results and D. longipes appears to be composed of two different entities.

Phylogeny of colletid bees (Hymenoptera: Colletidae) inferred from four nuclear genes

Colletidae comprise approximately 2500 species of bees primarily distributed in the southern continents (only two colletid genera are widely distributed: Colletes and Hylaeus). Previously published studies have failed to resolve phylogenetic relationships on a worldwide basis and this has been a major barrier to the progress of research regarding systematics and evolution of colletid bees. For this study, data from four nuclear gene loci: elongation factor-1alpha (F2 copy), opsin, wingless, and 28S rRNA were analyzed for 122 species of colletid bees, representing all subfamilies and tribes currently recognized; 22 species belonging to three other bee families were used as outgroups. Bayesian, maximum likelihood, and parsimony methods were employed to investigate the phylogenetic relationships within Colletidae and resulted in highly congruent and well-resolved trees. The phylogenetic results show that Colletidae are monophyletic and that all traditionally recognized subfamilies (except Paracolletinae) are also strongly supported as monophyletic. Our phylogenetic hypothesis provides a framework within which broad questions related to the taxonomy, biogeography, morphology, evolution, and ecology of colletid bees can be addressed.

Molecular phylogeny of reed beetles (Col., Chrysomelidae, Donaciinae): the signature of ecological specialization and geographical isolation

The Donaciinae consist of approximately 165 species predominantly occurring in the northern hemisphere. We analysed mitochondrial and nuclear DNA (COI, EF-1alpha) of 46 species to investigate their phylogeny and to discuss general topics in the context of insect herbivory (generalists versus specialists, ecological speciation). Phylogenetic reconstructions from various methodical approaches yielded very similar results. Clades corresponding to the traditional tribes/genera were recovered. Within the genus Donacia, species groups with characteristic host plant preference were identified. Estimated divergence times are discussed on the background of geological events. The origin of the Donaciinae is dated to 75-100 million years before present, after which they quickly diversified into the main groups. An initial split of those groups occurred in the Palaeocene. In the Eocene and Oligocene, major lineages specialized on certain host plants, where they radiated in the Miocene. This radiation was enforced by geographic isolation brought about by the final separation of America and Europe, after which there arose continental lineages within three larger species groups. In their evolution based on ecological specialization with a recently superimposed geographic isolation, the Donaciinae follow a pattern of specialists arising from generalists. Host plant shifts show that such a specialization is not necessarily an 'evolutionary dead-end'.

Origins and Evolution of Spliceosomal Introns

Research into the origins of introns is at a critical juncture in the resolution of theories on the evolution of early life (which came first, RNA or DNA?), the identity of LUCA (the last universal common ancestor, was it prokaryotic- or eukaryotic-like?), and the significance of noncoding nucleotide variation. One early notion was that introns would have evolved as a component of an efficient mechanism for the origin of genes. But alternative theories emerged as well. From the debate between the "introns-early" and "introns-late" theories came the proposal that introns arose before the origin of genetically encoded proteins and DNA, and the more recent "introns-first" theory, which postulates the presence of introns at that early evolutionary stage from a reconstruction of the "RNA world." Here we review seminal and recent ideas about intron origins. Recent discoveries about the patterns and causes of intron evolution make this one of the most hotly debated and exciting topics in molecular evolutionary biology today.

Tempo and Mode of Spliceosomal Intron Evolution in Actin of Foraminifera

Spliceosomal introns are present in almost all eukaryotic genes, yet little is known about their origin and turnover in the majority of eukaryotic phyla. There is no agreement whether most introns are ancestral and have been lost in some lineage or have been gained recently. We addressed this question by analyzing the spatial and temporal distribution of introns in actins of foraminifera, a group of testate protists whose exceptionally rich fossil record permits the calibration of molecular phylogenies to date intron origins. We identified 24 introns dispersed along the sequence of two foraminiferan actin paralogues and actin deviating proteins, an unconventional type of fast-evolving actin found in some foraminifera. Comparison of intron positions indicates that 20 of 24 introns are specific to foraminifera. Four introns shared between foraminifera and other eukaryotes were interpreted as parallel gains because they have been found only in single species belonging to phylogenetically distinctive lineages. Moreover, additional recent intron gain due to the transfer between the actin paralogues was observed in two cultured species. Based on a relaxed molecular clock timescale, we conclude that intron gains in actin took place throughout the evolution of foraminifera, with the oldest introns inserted between 550 and 500 million years ago and the youngest ones acquired less than 100 million years ago.

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.

A novel plant major intrinsic protein in Physcomitrella patens most similar to bacterial glycerol channels

A gene encoding a novel fifth type of major intrinsic protein (MIP) in plants has been identified in the moss Physcomitrella patens. Phylogenetic analyses show that this protein, GlpF-like intrinsic protein (GIP1;1), is closely related to a subclass of glycerol transporters in bacteria that in addition to glycerol are highly permeable to water. A likely explanation of the occurrence of this bacterial-like MIP in P. patens is horizontal gene transfer. The expressed P. patens GIP1;1 gene contains five introns and encodes a unique C-loop extension of approximately 110 amino acid residues that has no obvious similarity with any other known protein. Based on alignments and structural comparisons with other MIPs, GIP1;1 is suggested to have retained the permeability for glycerol but not for water. Studies on heterologously expressed GIP1;1 in Xenopus laevis oocytes confirm the predicted substrate specificity. Interestingly, proteins of one of the plant-specific subgroups of MIPs, the NOD26-like intrinsic proteins, are also facilitating the transport of glycerol and have previously been suggested to have evolved from a horizontally transferred bacterial gene. Further studies on localization and searches for GIP1;1 homologs in other plants will clarify the function and significance of this new plant MIP.

Higher-level bee classifications (Hymenoptera, Apoidea, Apidae sensu lato)

A higher-level classification of bees, in which the entire group is treated as a single family - the Apidae - is advocated here. A total of seven subfamilies, 51 tribes and 27 subtribes are recognized. These subfamilies correspond to the families adopted in the traditional classification. Although the proposed changes do not involve any major rearrangement, basically only changing the rank given to the main groups, the new system makes the classification of bees more consistent with that adopted for other major groups of aculeate Hymenoptera. It also departs from the 19th century practice, perpetuated in the traditional classification, of giving family-status to the main groups of bees. A correspondence table associating the taxon names used in the current traditional classification with those of the revised classification is presented. Scrapterini new tribe (type-genus Scrapter Lepeletier & Serville) is proposed to accommodate the southern African genus Scrapter.

Phylogenetic Mapping of Intron Positions: A Case Study of Translation Initiation Factor eIF2γ

Eukaryotic translation initiation factor 2 (eIF2) is a G protein that delivers the methionyl initiator tRNA to the small ribosomal subunit and releases it upon GTP hydrolysis after the recognition of the initiation codon. eIF2 is composed of three subunits, alpha, beta, and gamma. Subunit gamma shows the strongest conservation, and it confers both tRNA and GTP/GDP binding. Using intron positioning and protein sequence alignment, here we show that eIF2gamma is a suitable phylogenetic marker for eukaryotes. We determined or completed the sequences of 13 arthropod eIF2gamma genes. Analyzing the phylogenetic distribution of 52 different intron positions in 55 distantly related eIF2gamma genes, we identified ancient ones and shared derived introns in our data set. Obviously, intron positioning in eIF2gamma is evolutionarily conserved. However, there were episodes of complete and partial intron losses followed by intron gains. We identified 17 clusters of intron positions based on their distribution. The evolution of these clusters appears to be connected with preferred exon length and can be used to estimate the relative timing of intron gain because nearby precursor introns had to be erased from the gene before the new introns could be inserted. Moreover, we identified a putative case of intron sliding that constitutes a synapomorphic character state supporting monophyly of Coleoptera, Lepidoptera, and Diptera excluding Hymenoptera. We also performed tree reconstructions using the eIF2gamma protein sequences and intron positioning as phylogenetic information. Our results support the monophyly of Viridoplantae, Ascomycota, Homobasidiomyceta, and Apicomplexa.