Search

BackgroundAngiosperms employ an astonishing variety of visual and olfactory floral signals that are generally thought to evolve under natural selection. Those morphological and chemical traits can form highly correlated sets of traits. It is not always clear which of these are used by pollinators as primary targets of selection and which would be indirectly selected by being linked to those primary targets. Quantitative genetics tools for predicting multiple traits response to selection have been developed since long and have advanced our understanding of evolution of genetically correlated traits in various biological systems. We use these tools to predict the evolutionary trajectories of floral traits and understand the selection pressures acting on them.ResultsWe used data from an artificial selection and a pollinator (bumblebee, hoverfly) evolution experiment with fast cycling Brassica rapa plants to predict evolutionary changes of 12 floral volatiles and 4 morphological floral traits in response to selection. Using the observed selection gradients and the genetic variance-covariance matrix (G-matrix) of the traits, we showed that the observed responses of most floral traits including volatiles were predicted in the right direction in both artificial- and bumblebee-selection experiment. Genetic covariance had a mix of constraining and facilitating effects on evolutionary responses. We further revealed that G-matrices also evolved in the selection processes.ConclusionsOverall, our integrative study shows that floral signals, especially volatiles, evolve under selection in a mostly predictable way, at least during short term evolution. Evolutionary constraints stemming from genetic covariance affected traits evolutionary trajectories and thus it is important to include genetic covariance for predicting the evolutionary changes of a comprehensive suite of traits. Other processes such as resource limitation and selfing also need to be considered for a better understanding of floral trait evolution.

BackgroundOne hypothesis for the function of sleep is that it serves as a mechanism to conserve energy. Recent studies have suggested that increased sleep can be an adaptive mechanism to improve survival under food deprivation in Drosophila melanogaster. To test the generality of this hypothesis, we compared sleep and its plastic response to starvation in a temperate and tropical population of Drosophila melanogaster.ResultsWe found that flies from the temperate population were more starvation resistant, and hypothesized that they would engage in behaviors that are considered to conserve energy, including increased sleep and reduced movement. Surprisingly, temperate flies slept less and moved more when they were awake compared to tropical flies, both under fed and starved conditions, therefore sleep did not correlate with population-level differences in starvation resistance. In contrast, total sleep and percent change in sleep when starved were strongly positively correlated with starvation resistance within the tropical population, but not within the temperate population. Thus, we observe unexpectedly complex relationships between starvation and sleep that vary both within and across populations. These observations falsify the simple hypothesis of a straightforward relationship between sleep and energy conservation. We also tested the hypothesis that starvation is correlated with metabolic phenotypes by investigating stored lipid and carbohydrate levels, and found that stored metabolites partially contributed towards variation starvation resistance.ConclusionsOur findings demonstrate that the function of sleep under starvation can rapidly evolve on short timescales and raise new questions about the physiological correlates of sleep and the extent to which variation in sleep is shaped by natural selection.

BackgroundEukaryotes acquired the trait of oxygenic photosynthesis through endosymbiosis of the cyanobacterial progenitor of plastid organelles. Despite recent advances in the phylogenomics of Cyanobacteria, the phylogenetic root of plastids remains controversial. Although a single origin of plastids by endosymbiosis is broadly supported, recent phylogenomic studies are contradictory on whether plastids branch early or late within Cyanobacteria. One underlying cause may be poor fit of evolutionary models to complex phylogenomic data.ResultsUsing Posterior Predictive Analysis, we show that recently applied evolutionary models poorly fit three phylogenomic datasets curated from cyanobacteria and plastid genomes because of heterogeneities in both substitution processes across sites and of compositions across lineages. To circumvent these sources of bias, we developed CYANO-MLP, a machine learning algorithm that consistently and accurately phylogenetically classifies ("phyloclassifies") cyanobacterial genomes to their clade of origin based on bioinformatically predicted function-informative features in tRNA gene complements. Classification of cyanobacterial genomes with CYANO-MLP is accurate and robust to deletion of clades, unbalanced sampling, and compositional heterogeneity in input tRNA data. CYANO-MLP consistently classifies plastid genomes into a late-branching cyanobacterial sub-clade containing single-cell, starch-producing, nitrogen-fixing ecotypes, consistent with metabolic and gene transfer data.ConclusionsPhylogenomic data of cyanobacteria and plastids exhibit both site-process heterogeneities and compositional heterogeneities across lineages. These aspects of the data require careful modeling to avoid bias in phylogenomic estimation. Furthermore, we show that amino acid recoding strategies may be insufficient to mitigate bias from compositional heterogeneities. However, the combination of our novel tRNA-specific strategy with machine learning in CYANO-MLP appears robust to these sources of bias with high accuracy in phyloclassification of cyanobacterial genomes. CYANO-MLP consistently classifies plastids as late-branching Cyanobacteria, consistent with independent evidence from signature-based approaches and some previous phylogenetic studies.

BackgroundPhylogenies provide critical information about convergence during adaptive radiation. To test whether there have been multiple origins of a distinctive trophic phenotype in one of the most rapidly radiating groups known, we used ultra-conserved elements (UCEs) to examine the evolutionary affinities of Lake Malawi cichlids lineages exhibiting greatly hypertrophied lips.ResultsThe hypertrophied lip cichlids Cheilochromis euchilus, Eclectochromis ornatus, Placidochromis "Mbenji fatlip", and Placidochromis milomo are all nested within the non-mbuna clade of Malawi cichlids based on both concatenated sequence and single nucleotide polymorphism (SNP) inferred phylogenies. Lichnochromis acuticeps that exhibits slightly hypertrophied lips also appears to have evolutionary affinities to this group. However, Chilotilapia rhoadesii that lacks hypertrophied lips was recovered as nested within the species Cheilochromis euchilus. Species tree reconstructions and analyses of introgression provided largely ambiguous patterns of Malawi cichlid evolution.ConclusionsContrary to mitochondrial DNA phylogenies, bifurcating trees based on our 1024 UCE loci supported close affinities of Lake Malawi lineages with hypertrophied lips. However, incomplete lineage sorting in Malawi tends to render these inferences more tenuous. Phylogenomic analyses will continue to provide powerful inferences about whether phenotypic novelties arose once or multiple times during adaptive radiation.

BackgroundAdvances in genomic technologies have expanded our ability to accurately and exhaustively detect natural genomic variants that can be applied in crop improvement and to increase our knowledge of plant evolution and adaptation. Switchgrass (Panicum virgatum L.), an allotetraploid (2n = 4× = 36) perennial C4 grass (Poaceae family) native to North America and a feedstock crop for cellulosic biofuel production, has a large potential for genetic improvement due to its high genotypic and phenotypic variation. In this study, we analyzed single nucleotide polymorphism (SNP) variation in 372 switchgrass genotypes belonging to 36 accessions for 12 genes putatively involved in biomass production to investigate signatures of selection that could have led to ecotype differentiation and to population adaptation to geographic zones.ResultsA total of 11,682 SNPs were mined from ~ 15 Gb of sequence data, out of which 251 SNPs were retained after filtering. Population structure analysis largely grouped upland accessions into one subpopulation and lowland accessions into two additional subpopulations. The most frequent SNPs were in homozygous state within accessions. Sixty percent of the exonic SNPs were non-synonymous and, of these, 45% led to non-conservative amino acid changes. The non-conservative SNPs were largely in linkage disequilibrium with one haplotype being predominantly present in upland accessions while the other haplotype was commonly present in lowland accessions. Tajima's test of neutrality indicated that PHYB, a gene involved in photoperiod response, was under positive selection in the switchgrass population. PHYB carried a SNP leading to a non-conservative amino acid change in the PAS domain, a region that acts as a sensor for light and oxygen in signal transduction.ConclusionsSeveral non-conservative SNPs in genes potentially involved in plant architecture and adaptation have been identified and led to population structure and genetic differentiation of ecotypes in switchgrass. We suggest here that PHYB is a key gene involved in switchgrass natural selection. Further analyses are needed to determine whether any of the non-conservative SNPs identified play a role in the differential adaptation of upland and lowland switchgrass.

BackgroundAn important feature of eukaryotic evolution is metabolic compartmentalization, in which certain pathways are restricted to the cytosol or specific organelles. Glycolysis in eukaryotes is described as a cytosolic process. The universality of this canon has been challenged by recent genome data that suggest that some glycolytic enzymes made by stramenopiles bear mitochondrial targeting peptides.ResultsMining of oomycete, diatom, and brown algal genomes indicates that stramenopiles encode two forms of enzymes for the second half of glycolysis, one with and the other without mitochondrial targeting peptides. The predicted mitochondrial targeting was confirmed by using fluorescent tags to localize phosphoglycerate kinase, phosphoglycerate mutase, and pyruvate kinase in Phytophthora infestans, the oomycete that causes potato blight. A genome-wide search for other enzymes with atypical mitochondrial locations identified phosphoglycerate dehydrogenase, phosphoserine aminotransferase, and phosphoserine phosphatase, which form a pathway for generating serine from the glycolytic intermediate 3-phosphoglycerate. Fluorescent tags confirmed the delivery of these serine biosynthetic enzymes to P. infestans mitochondria. A cytosolic form of this serine biosynthetic pathway, which occurs in most eukaryotes, is missing from oomycetes and most other stramenopiles. The glycolysis and serine metabolism pathways of oomycetes appear to be mosaics of enzymes with different ancestries. While some of the noncanonical oomycete mitochondrial enzymes have the closest affinity in phylogenetic analyses with proteins from other stramenopiles, others cluster with bacterial, plant, or animal proteins. The genes encoding the mitochondrial phosphoglycerate kinase and serine-forming enzymes are physically linked on oomycete chromosomes, which suggests a shared origin.ConclusionsStramenopile metabolism appears to have been shaped through the acquisition of genes by descent and lateral or endosymbiotic gene transfer, along with the targeting of the proteins to locations that are novel compared to other eukaryotes. Colocalization of the glycolytic and serine biosynthesis enzymes in mitochondria is apparently necessary since they share a common intermediate. The results indicate that descriptions of metabolism in textbooks do not cover the full diversity of eukaryotic biology.

BackgroundSiglecs-11 and -16 are members of the sialic acid recognizing Ig-like lectin family, and expressed in same cells. Siglec-11 functions as an inhibitory receptor, whereas Siglec-16 exhibits activating properties. In humans, SIGLEC11 and SIGLEC16 gene sequences are extremely similar in the region encoding the extracellular domain due to gene conversions. Human SIGLEC11 was converted by the nonfunctional SIGLEC16P allele, and the converted SIGLEC11 allele became fixed in humans, possibly because it provides novel neuroprotective functions in brain microglia. However, the detailed evolutionary history of SIGLEC11 and SIGLEC16 in other primates remains unclear.ResultsWe analyzed SIGLEC11 and SIGLEC16 gene sequences of multiple primate species, and examined glycan binding profiles of these Siglecs. The phylogenetic tree demonstrated that gene conversions between SIGLEC11 and SIGLEC16 occurred in the region including the exon encoding the sialic acid binding domain in every primate examined. Functional assays showed that glycan binding preference is similar between Siglec-11 and Siglec-16 in all analyzed hominid species. Taken together with the fact that Siglec-11 and Siglec-16 are expressed in the same cells, Siglec-11 and Siglec-16 are regarded as paired receptors that have maintained similar ligand binding preferences via gene conversions. Relaxed functional constraints were detected on the SIGLEC11 and SIGLEC16 exons that underwent gene conversions, possibly contributing to the evolutionary acceptance of repeated gene conversions. The frequency of nonfunctional SIGLEC16P alleles is much higher than that of SIGLEC16 alleles in every human population.ConclusionsOur findings indicate that Siglec-11 and Siglec-16 have been maintained as paired receptors by repeated gene conversions under relaxed functional constraints in the primate lineage. The high prevalence of the nonfunctional SIGLEC16P allele and the fixation of the converted SIGLEC11 imply that the loss of Siglec-16 and the gain of Siglec-11 in microglia might have been favored during the evolution of human lineage.

BackgroundLongwing butterflies, Heliconius sp., also called heliconians, are striking examples of diversity and mimicry in butterflies. Heliconians feature strongly colored patterns on their wings, arising from wing scales colored by pigments and/or nanostructures, which serve as an aposematic signal.ResultsHere, we investigate the coloration mechanisms among several species of Heliconius by applying scanning electron microscopy, (micro)spectrophotometry, and imaging scatterometry. We identify seven kinds of colored scales within Heliconius whose coloration is derived from pigments, nanostructures or both. In yellow-, orange- and red-colored wing patches, both cover and ground scales contain wavelength-selective absorbing pigments, 3-OH-kynurenine, xanthommatin and/or dihydroxanthommatin. In blue wing patches, the cover scales are blue either due to interference of light in the thin-film lower lamina (e.g., H. doris) or in the multilayered lamellae in the scale ridges (so-called ridge reflectors, e.g., H. sara and H. erato); the underlying ground scales are black. In the white wing patches, both cover and ground scales are blue due to their thin-film lower lamina, but because they are stacked upon each other and at the wing substrate, a faint bluish to white color results. Lastly, green wing patches (H. doris) have cover scales with blue-reflecting thin films and short-wavelength absorbing 3-OH-kynurenine, together causing a green color.ConclusionsThe pigmentary and structural traits are discussed in relation to their phylogenetic distribution and the evolution of vision in this highly interesting clade of butterflies.

BackgroundThe Spiralia are a large, morphologically diverse group of protostomes (e.g. molluscs, annelids, nemerteans) that share a homologous mode of early development called spiral cleavage. One of the most highly-conserved features of spiralian development is the contribution of the primary quartet cells, 1a-1d, to the anterior region of the embryo (including the brain, eyes, and the anterior ciliary band, called the prototroch). Yet, very few studies have analyzed the ultimate fates of primary quartet sub-lineages, or examined the morphogenetic events that take place in the anterior region of the embryo.ResultsThis study focuses on the caenogastropod slipper snail, Crepidula fornicata, a model for molluscan developmental biology. Through direct lineage tracing of primary quartet daughter cells, and examination of these cells during gastrulation and organogenesis stages, we uncovered behaviors never described before in a spiralian. For the first time, we show that the 1a2-1d2 cells do not contribute to the prototroch (as they do in other species) and are ultimately lost before hatching. During gastrulation and anterior-posterior axial elongation stages, these cells cleavage-arrest and spread dramatically, contributing to a thin provisional epidermis on the dorsal side of the embryo. This spreading is coupled with the displacement of the animal pole, and other pretrochal cells, closer to the ventrally-positioned mouth, and the vegetal pole.ConclusionsThis is the first study to document the behavior and fate of primary quartet sub-lineages among molluscs. We speculate that the function of 1a2-1d2 cells (in addition to two cells derived from 1d12, and the 2b lineage) is to serve as a provisional epithelium that allows for anterior displacement of the other progeny of the primary quartet towards the anterior-ventral side of the embryo. These data support a new and novel mechanism for axial bending, distinct from canonical models in which axial bending is suggested to be driven primarily by differential proliferation of posterior dorsal cells. These data suggest also that examining sub-lineages in other spiralians will reveal greater variation than previously assumed.

BackgroundThe herbivore lifestyle leads to encounters with plant toxins and requires mechanisms to overcome suboptimal nutrient availability in plant tissues. Although the evolution of bacterial endosymbiosis alleviated many of these challenges, the ability to manipulate plant nutrient status has evolved in lineages with and without nutritional symbionts. Whether and how these alternative nutrient acquisition strategies interact or constrain insect evolution is unknown. We studied the transcriptomes of galling and free-living aphidomorphs to characterize how amino acid transporter evolution is influenced by the ability to manipulate plant resource availability.ResultsUsing a comparative approach we found phylloxerids retain nearly all amino acid transporters as other aphidomorphs, despite loss of nutritional endosymbiosis. Free living species show more transporters than galling species within the same genus, family, or infraorder, indicating plant hosts influence the maintenance and evolution of nutrient transport within herbivores. Transcript profiles also show lineage specificity and suggest some genes may facilitate life without endosymbionts or the galling lifestyle.ConclusionsThe transcript abundance profiles we document across fluid feeding herbivores support plant host constraint on insect amino acid transporter evolution. Given amino acid uptake, transport, and catabolism underlie the success of herbivory as a life history strategy, this suggests that plant host nutrient quality, whether constitutive or induced, alters the selective environment surrounding the evolution and maintenance of endosymbiosis.

BackgroundFrom bat wings to whale flippers, limb diversification has been crucial to the evolutionary success of mammals. We performed the first transcriptome-wide study of limb development in multiple species to explore the hypothesis that mammalian limb diversification has proceeded through the differential expression of conserved shared genes, rather than by major changes to limb patterning. Specifically, we investigated the manner in which the expression of shared genes has evolved within and among mammalian species.ResultsWe assembled and compared transcriptomes of bat, mouse, opossum, and pig fore- and hind limbs at the ridge, bud, and paddle stages of development. Results suggest that gene expression patterns exhibit larger variation among species during later than earlier stages of limb development, while within species results are more mixed. Consistent with the former, results also suggest that genes expressed at later developmental stages tend to have a younger evolutionary age than genes expressed at earlier stages. A suite of key limb-patterning genes was identified as being differentially expressed among the homologous limbs of all species. However, only a small subset of shared genes is differentially expressed in the fore- and hind limbs of all examined species. Similarly, a small subset of shared genes is differentially expressed within the fore- and hind limb of a single species and among the forelimbs of different species.ConclusionsTaken together, results of this study do not support the existence of a phylotypic period of limb development ending at chondrogenesis, but do support the hypothesis that the hierarchical nature of development translates into increasing variation among species as development progresses.

BackgroundMolecular markers are revealing a much more diverse and evolutionarily complex picture of marine biodiversity than previously anticipated. Cryptic and/or endemic marine species are continually being found throughout the world oceans, predominantly in inconspicuous tropical groups but also in larger, canopy-forming taxa from well studied temperate regions. Interspecific hybridization has also been found to be prevalent in many marine groups, for instance within dense congeneric assemblages, with introgressive gene-flow being the most common outcome. Here, using a congeneric phylogeographic approach, we investigated two monotypic and geographically complementary sister genera of north-east Pacific intertidal seaweeds (Hesperophycus and Pelvetiopsis), for which preliminary molecular tests revealed unexpected conflicts consistent with unrecognized cryptic diversity and hybridization.ResultsThe three recovered mtDNA clades did not match a priori species delimitations. H. californicus was congruent, whereas widespread P. limitata encompassed two additional narrow-endemic species from California - P. arborescens (here genetically confirmed) and P. hybrida sp. nov. The congruence between the genotypic clusters and the mtDNA clades was absolute. Fixed heterozygosity was apparent in a high proportion of loci in P. limitata and P. hybrida, with genetic analyses showing that the latter was composed of both H. californicus and P. arborescens genomes. All four inferred species could be distinguished based on their general morphology.ConclusionsThis study confirmed additional diversity and reticulation within NE Pacific Hesperophycus/Pelvetiopsis, including the validity of the much endangered, modern climatic relict P. arborescens, and the identification of a new, stable allopolyploid species (P. hybrida) with clearly discernable ancestry (♀ H. californicus x ♂ P. arborescens), morphology, and geographical distribution. Allopolyploid speciation is otherwise completely unknown in brown seaweeds, and its unique occurrence within this genus (P. limitata possibly representing a second example) remains enigmatic. The taxonomic separation of Hesperophycus and Pelvetiopsis is not supported and the genera should be synonymized; we retain only the latter. The transitional coastline between Point Conception and Monterey Bay represented a diversity hotspot for the genus and the likely sites of extraordinary evolutionary events of allopolyploid speciation at sympatric range contact zones. This study pinpoints how much diversity (and evolutionary processes) potentially remains undiscovered even on a conspicuous seaweed genus from the well-studied Californian intertidal shores let alone in other, less studied marine groups and regions/depths.

BackgroundPigeons and doves (Columbiformes) are one of the oldest and most diverse extant lineages of birds. However, the nature and timing of the group's evolutionary radiation remains poorly resolved, despite recent advances in DNA sequencing and assembly and the growing database of pigeon mitochondrial genomes. One challenge has been to generate comparative data from the large number of extinct pigeon lineages, some of which are morphologically unique and therefore difficult to place in a phylogenetic context.ResultsWe used ancient DNA and next generation sequencing approaches to assemble complete mitochondrial genomes for eleven pigeons, including the extinct Ryukyu wood pigeon (Columba jouyi), the thick-billed ground dove (Alopecoenas salamonis), the spotted green pigeon (Caloenas maculata), the Rodrigues solitaire (Pezophaps solitaria), and the dodo (Raphus cucullatus). We used a Bayesian approach to infer the evolutionary relationships among 24 species of living and extinct pigeons and doves.ConclusionsOur analyses indicate that the earliest radiation of the Columbidae crown group most likely occurred during the Oligocene, with continued divergence of major clades into the Miocene, suggesting that diversification within the Columbidae occurred more recently than has been reported previously.

BackgroundFlatfish cranial asymmetry represents one of the most remarkable morphological innovations among vertebrates, and has fueled vigorous debate on the manner and rate at which strikingly divergent phenotypes evolve. A surprising result of many recent molecular phylogenetic studies is the lack of support for flatfish monophyly, where increasingly larger DNA datasets of up to 23 loci have either yielded a weakly supported flatfish clade or indicated the group is polyphyletic. Lack of resolution for flatfish relationships has been attributed to analytical limitations for dealing with processes such as nucleotide non-stationarity and incomplete lineage sorting (ILS). We tackle this phylogenetic problem using a sequence dataset comprising more than 1,000 ultraconserved DNA element (UCE) loci covering 45 carangimorphs, the broader clade containing flatfishes and several other specialized lineages such as remoras, billfishes, and archerfishes.ResultsWe present a phylogeny based on UCE loci that unequivocally supports flatfish monophyly and a single origin of asymmetry. We document similar levels of discordance among UCE loci as in previous, smaller molecular datasets. However, relationships among flatfishes and carangimorphs recovered from multilocus concatenated and species tree analyses of our data are robust to the analytical framework applied and size of data matrix used. By integrating the UCE data with a rich fossil record, we find that the most distinctive carangimorph bodyplans arose rapidly during the Paleogene (66.0-23.03 Ma). Flatfish asymmetry, for example, likely evolved over an interval of no more than 2.97 million years.ConclusionsThe longstanding uncertainty in phylogenetic hypotheses for flatfishes and their carangimorph relatives highlights the limitations of smaller molecular datasets when applied to successive, rapid divergences. Here, we recovered significant support for flatfish monophyly and relationships among carangimorphs through analysis of over 1,000 UCE loci. The resulting time-calibrated phylogeny points to phenotypic divergence early within carangimorph history that broadly matches with the predictions of adaptive models of lineage diversification.

BackgroundSMRT and NCoR are corepressor paralogs that help mediate transcriptional repression by a variety of transcription factors, including the nuclear hormone receptors. The functions of both corepressors are extensively diversified in mice by alternative mRNA splicing, generating a series of protein variants that differ in different tissues and that exert different, even diametrically opposite, biochemical and biological effects from one another.ResultsWe report here that the alternative splicing previously reported for SMRT appears to be a relatively recent evolutionary phenomenon, with only one of these previously identified sites utilized in a teleost fish and a limited additional number of the additional known sites utilized in a bird, reptile, and marsupial. In contrast, extensive SMRT alternative splicing at these sites was detected among the placental mammals. The alternative splicing of NCoR previously identified in mice (and shown to regulate lipid and carbohydrate metabolism) is likely to have arisen separately and after that of SMRT, and includes an example of convergent evolution.ConclusionsWe propose that the functions of both SMRT and NCoR have been diversified by alternative splicing during evolution to allow customization for different purposes in different tissues and different species.

BackgroundDarwin and the architects of the Modern Synthesis found sympatric speciation difficult to explain and suggested it is unlikely to occur. Increasingly, evidence over the past few decades suggest that sympatric speciation can occur under ecological conditions that require at most intraspecific competition for a structured resource. Here we used an individual-based population model with variable foraging strategies to study the evolution of mating behavior among foraging strategy types. Initially, individuals were placed at random on a structureless resource landscape, with subsequent spatial variation induced through foraging activity itself. The fitness of individuals was determined by their biomass at the end of each generational cycle. The model incorporates three diallelic, codominant foraging strategy genes, and one mate-choice or m-trait (i.e. incipient magic trait) gene, where the latter is inactive when random mating is assumed.ResultsUnder non-random mating, the m-trait gene promotes increasing levels of either disassortative or assortative mating when the frequency of m respectively increases or decreases from 0.5. Our evolutionary simulations demonstrate that, under initial random mating conditions, an activated m-trait gene evolves to promote assortative mating because the system, in trying to fit a multipeak adaptive landscape, causes heterozygous individuals to be less fit than homozygous individuals.ConclusionOur results extend our theoretical understanding that sympatric speciation can evolve under nicheless or gradientless resource conditions: i.e. the underlying resource is monomorphic and initially spatially homogeneous. Further the simplicity and generality of our model suggests that sympatric speciation may be more likely than previously thought to occur in mobile, sexually-reproducing organisms.

BackgroundUltraconserved elements (UCEs) have been successfully used in phylogenomics for a variety of taxa, but their power in phylogenetic inference has yet to be extensively compared with that of traditional Sanger sequencing data sets. Moreover, UCE data on invertebrates, including insects, are sparse. We compared the phylogenetic informativeness of 959 UCE loci with a multi-locus data set of ten nuclear markers obtained via Sanger sequencing, testing the ability of these two types of data to resolve and date the evolutionary history of the second most species-rich subfamily of ants in the world, the Formicinae.ResultsPhylogenetic analyses show that UCEs are superior in resolving ancient and shallow relationships in formicine ants, demonstrated by increased node support and a more resolved phylogeny. Phylogenetic informativeness metrics indicate a twofold improvement relative to the 10-gene data matrix generated from the identical set of taxa. We were able to significantly improve formicine classification based on our comprehensive UCE phylogeny. Our divergence age estimations, using both UCE and Sanger data, indicate that crown-group Formicinae are older (104-117 Ma) than previously suggested. Biogeographic analyses infer that the diversification of the subfamily has occurred on all continents with no particular hub of cladogenesis.ConclusionsWe found UCEs to be far superior to the multi-locus data set in estimating formicine relationships. The early history of the clade remains uncertain due to ancient rapid divergence events that are unresolvable even with our genomic-scale data, although this might be largely an effect of several problematic taxa subtended by long branches. Our comparison of divergence ages from both Sanger and UCE data demonstrates the effectiveness of UCEs for dating analyses. This comparative study highlights both the promise and limitations of UCEs for insect phylogenomics, and will prove useful to the growing number of evolutionary biologists considering the transition from Sanger to next-generation sequencing approaches.

BackgroundPhylogeography is an important tool that can be used to reveal cryptic biodiversity and to better understand the processes that promote lineage diversification. We studied the phylogeographic history of the Arboreal Salamander (Aneides lugubris), a wide-ranging species endemic to the California floristic province. We used multi-locus data to reconstruct the evolutionary history of A. lugubris and to discover the geographic location of major genetic breaks within the species. We also used species distribution modeling and comparative phylogeography to better understand the environmental factors that have shaped the genetic history of A. lugubris.ResultsWe found six major mitochondrial clades in A. lugubris. Nuclear loci supported the existence of at least three genetically distinct groups, corresponding to populations north of the San Francisco Bay and in the Sierra Nevada, in the Santa Cruz Mountains, and in the central coast and southern California. All of the genetic breaks in mitochondrial and nuclear loci corresponded to regions where historical barriers to dispersal have been observed in other species. Geologic or water barriers likely were the most important factors restricting gene flow among clades. Climatic unsuitability during glacial maximum may have contributed to the isolation of the mitochondrial clades in the central coast and southern California. A projection of our species distribution model to a future scenario with a moderate amount of climate change suggests that most of the range of A. lugubris will remain climatically suitable, but climatic conditions in the Sierra Nevada and low elevation areas in Southern California are likely to deteriorate.ConclusionsAneides lugubris contains substantial cryptic genetic diversity as a result of historical isolation of populations. At least two (and perhaps three) evolutionarily significant units in A. lugubris merit protection; all six mitochondrial clades should be considered as management units within the species.

BackgroundShort-term experiments have indicated that warmer temperatures can alter fungal biomass production and CO2 respiration, with potential consequences for soil C storage. However, we know little about the capacity of fungi to adapt to warming in ways that may alter C dynamics. Thus, we exposed Neurospora discreta to moderately warm (16 °C) and warm (28 °C) selective temperatures for 1500 mitotic generations, and then examined changes in mycelial growth rate, biomass, spore production, and CO2 respiration. We tested the hypothesis that strains will adapt to its selective temperature. Specifically, we expected that adapted strains would grow faster, and produce more spores per unit biomass (i.e., relative spore production). In contrast, they should generate less CO2 per unit biomass due to higher efficiency in carbon use metabolism (i.e., lower mass specific respiration, MSR).ResultsIndeed, N. discreta adapted to warm temperatures, based on patterns of relative spore production. Adapted strains produced more spores per unit biomass than parental strains in the selective temperature. Contrary to our expectations, this increase in relative spore production was accompanied by an increase in MSR and a reduction in mycelial growth rate and biomass, compared to parental strains.ConclusionsAdaptation of N. discreta to warm temperatures may have elicited a tradeoff between biomass production and relative spore production, possibly because relative spore production required higher MSR rates. Therefore, our results do not support the idea that adaptation to warm temperatures will lead to a more efficient carbon use metabolism. Our data might help improve climate change model simulations and provide more concise predictions of decomposition processes and carbon feedbacks to the atmosphere.

BACKGROUND:In the short-term, organisms acclimate to stress through phenotypic plasticity, but in the longer term they adapt to stress genetically. The mutations that accrue during adaptation may contribute to completely novel phenotypes, or they may instead act to restore the phenotype from a stressed to a pre-stress condition. To better understand the influence of evolution on the diversity and direction of phenotypic change, we used Biolog microarrays to assay 94 phenotypes of 115 Escherichia coli clones that had adapted to high temperature (42.2 °C). We also assayed these same phenotypes in the clones' ancestor under non-stress (37.0 °C) and stress (42.2 °C) conditions. We explored associations between Biolog phenotypes and genotypes, and we also investigated phenotypic differences between clones that have one of two adaptive genetic trajectories: one that is typified by mutations in the RNA polymerase β-subunit (rpoB) and another that is defined by mutations in the rho termination factor. RESULTS:Most (54 %) phenotypic variation was restorative, shifting the phenotype from the acclimated state back toward the unstressed state. Novel phenotypes were more rare, comprising between 5 and 18 % of informative phenotypic variation. Phenotypic variation associated statistically with genetic variation, demonstrating a genetic basis for phenotypic change. Finally, clones with rpoB mutations differed in phenotype from those with rho mutations, largely due to differences in chemical sensitivity. CONCLUSIONS:Our results contribute to previous observations showing that a major component of adaptation in microbial evolution experiments is toward restoration to the unstressed state. In addition, we found that a large deletion strongly affected phenotypic variation. Finally, we demonstrated that the two genetic trajectories leading to thermal adaptation encompass different phenotypes.

BackgroundInsects rely more on chemical signals (semiochemicals) than on any other sensory modality to find, identify, and choose mates. In most insects, pheromone production is typically regulated through biosynthetic pathways, whereas pheromone sensory detection is controlled by the olfactory system. Orchid bees are exceptional in that their semiochemicals are not produced metabolically, but instead male bees collect odoriferous compounds (perfumes) from the environment and store them in specialized hind-leg pockets to subsequently expose during courtship display. Thus, the olfactory sensory system of orchid bees simultaneously controls male perfume traits (sender components) and female preferences (receiver components). This functional linkage increases the opportunities for parallel evolution of male traits and female preferences, particularly in response to genetic changes of chemosensory detection (e.g. Odorant Receptor genes). To identify whether shifts in pheromone composition among related lineages of orchid bees are associated with divergence in chemosensory genes of the olfactory periphery, we searched for patterns of divergent selection across the antennal transcriptomes of two recently diverged sibling species Euglossa dilemma and E. viridissima.ResultsWe identified 3185 orthologous genes including 94 chemosensory loci from five different gene families (Odorant Receptors, Ionotropic Receptors, Gustatory Receptors, Odorant Binding Proteins, and Chemosensory Proteins). Our results revealed that orthologs with signatures of divergent selection between E. dilemma and E. viridissima were significantly enriched for chemosensory genes. Notably, elevated signals of divergent selection were almost exclusively observed among chemosensory receptors (i.e. Odorant Receptors).ConclusionsOur results suggest that rapid changes in the chemosensory gene family occurred among closely related species of orchid bees. These findings are consistent with the hypothesis that strong divergent selection acting on chemosensory receptor genes plays an important role in the evolution and diversification of insect pheromone systems.

BackgroundGene duplication provides raw material for the evolution of functional innovation. We recently developed a phylogenetic method that classifies evolutionary processes driving the retention of duplicate genes by quantifying divergence between their spatial gene expression profiles and that of their single-copy orthologous gene in a closely related sister species.ResultsHere, we apply our classification method to pairs of duplicate genes in eight mammalian genomes, using data from 11 tissues to construct spatial gene expression profiles. We find that young mammalian duplicates are often functionally conserved, and that expression divergence rapidly increases over evolutionary time. Moreover, expression divergence results in increased tissue specificity, with an overrepresentation of expression in male kidney, underrepresentation of expression in female liver, and strong underrepresentation of expression in testis. Thus, duplicate genes acquire a diversity of new tissue-specific functions outside of the testis, possibly contributing to the origin of a multitude of complex phenotypes during mammalian evolution.ConclusionsOur findings reveal that mammalian duplicate genes are initially functionally conserved, and then undergo rapid functional divergence over evolutionary time, acquiring diverse tissue-specific biological roles. These observations are in stark contrast to the much faster expression divergence and acquisition of broad housekeeping roles we previously observed in Drosophila duplicate genes. Due to the smaller effective population sizes of mammals relative to Drosophila, these analyses implicate natural selection in the functional evolution of duplicate genes.

BackgroundThe past decade has witnessed remarkable progress towards resolution of the Tree of Life. However, despite the increased use of genomic scale datasets, some phylogenetic relationships remain difficult to resolve. Here we employ anchored phylogenomics to capture 107 nuclear loci in 29 species of acanthomorph teleost fishes, with 25 of these species sampled from the recently delimited clade Ovalentaria. Previous studies employing multilocus nuclear exon datasets have not been able to resolve the nodes at the base of the Ovalentaria tree with confidence. Here we test whether a phylogenomic approach will provide better support for these nodes, and if not, why this may be.ResultsAfter using a novel method to account for paralogous loci, we estimated phylogenies with maximum likelihood and species tree methods using DNA sequence alignments of over 80,000 base pairs. Several key relationships within Ovalentaria are well resolved, including 1) the sister taxon relationship between Cichlidae and Pholidichthys, 2) a clade containing blennies, grammas, clingfishes, and jawfishes, and 3) monophyly of Atherinomorpha (topminnows, flyingfishes, and silversides). However, many nodes in the phylogeny associated with the early diversification of Ovalentaria are poorly resolved in several analyses. Through the use of rarefaction curves we show that limited phylogenetic resolution among the earliest nodes in the Ovalentaria phylogeny does not appear to be due to a deficiency of data, as average global node support ceases to increase when only 1/3rd of the sampled loci are used in analyses. Instead this lack of resolution may be driven by model misspecification as a Bayesian mixed model analysis of the amino acid dataset provided good support for parts of the base of the Ovalentaria tree.ConclusionsAlthough it does not appear that the limited phylogenetic resolution among the earliest nodes in the Ovalentaria phylogeny is due to a deficiency of data, it may be that both stochastic and systematic error resulting from model misspecification play a role in the poor resolution at the base of the Ovalentaria tree as a Bayesian approach was able to resolve some of the deeper nodes, where the other methods failed.

BackgroundAntarctic notothenioids are an impressive adaptive radiation. While they share recent common ancestry with several species-depauperate lineages that exhibit a relictual distribution in areas peripheral to the Southern Ocean, an understanding of their evolutionary origins and biogeographic history is limited as the sister lineage of notothenioids remains unidentified. The phylogenetic placement of notothenioids among major lineages of perciform fishes, which include sculpins, rockfishes, sticklebacks, eelpouts, scorpionfishes, perches, groupers and soapfishes, remains unresolved. We investigate the phylogenetic position of notothenioids using DNA sequences of 10 protein coding nuclear genes sampled from more than 650 percomorph species. The biogeographic history of notothenioids is reconstructed using a maximum likelihood method that integrates phylogenetic relationships, estimated divergence times, geographic distributions and paleogeographic history.ResultsPercophis brasiliensis is resolved, with strong node support, as the notothenioid sister lineage. The species is endemic to the subtropical and temperate Atlantic coast of southern South America. Biogeographic reconstructions imply the initial diversification of notothenioids involved the western portion of the East Gondwanan Weddellian Province. The geographic disjunctions among the major lineages of notothenioids show biogeographic and temporal correspondence with the fragmentation of East Gondwana.ConclusionsThe phylogenetic resolution of Percophis requires a change in the classification of percomorph fishes and provides evidence for a western Weddellian origin of notothenioids. The biogeographic reconstruction highlights the importance of the geographic and climatic isolation of Antarctica in driving the radiation of cold-adapted notothenioids.

BackgroundThe Mormon Metalmark (Apodemia mormo) species complex occurs as isolated and phenotypically variable colonies in dryland areas across western North America. Lange's Metalmark, A. m. langei, one of the 17 subspecies taxonomically recognized in the complex, is federally listed under the U.S. Endangered Species Act of 1973. Metalmark taxa have traditionally been described based on phenotypic and ecological characteristics, and it is unknown how well this nomenclature reflects their genetic and evolutionary distinctiveness. Genetic variation in six microsatellite loci and mitochondrial cytochrome oxidase subunit I sequence was used to assess the population structure of the A. mormo species complex across 69 localities, and to evaluate A. m. langei's qualifications as an Evolutionarily Significant Unit.ResultsWe discovered substantial genetic divergence within the species complex, especially across the Continental Divide, with population genetic structure corresponding more closely with geographic proximity and local isolation than with taxonomic divisions originally based on wing color and pattern characters. Lange's Metalmark was as genetically divergent as several other locally isolated populations in California, and even the unique phenotype that warranted subspecific and conservation status is reminiscent of the morphological variation found in some other populations.ConclusionsThis study is the first genetic treatment of the A. mormo complex across western North America and potentially provides a foundation for reassessing the taxonomy of the group. Furthermore, these results illustrate the utility of molecular markers to aid in demarcation of biological units below the species level. From a conservation point of view, Apodemia mormo langei's diagnostic taxonomic characteristics may, by themselves, not support its evolutionary significance, which has implications for its formal listing as an Endangered Species.

BackgroundBayesian relaxed-clock dating has significantly influenced our understanding of the timeline of biotic evolution. This approach requires the use of priors on the branching process, yet little is known about their impact on divergence time estimates. We investigated the effect of branching priors using the iconic cycads. We conducted phylogenetic estimations for 237 cycad species using three genes and two calibration strategies incorporating up to six fossil constraints to (i) test the impact of two different branching process priors on age estimates, (ii) assess which branching prior better fits the data, (iii) investigate branching prior impacts on diversification analyses, and (iv) provide insights into the diversification history of cycads.ResultsUsing Bayes factors, we compared divergence time estimates and the inferred dynamics of diversification when using Yule versus birth-death priors. Bayes factors were calculated with marginal likelihood estimated with stepping-stone sampling. We found striking differences in age estimates and diversification dynamics depending on prior choice. Dating with the Yule prior suggested that extant cycad genera diversified in the Paleogene and with two diversification rate shifts. In contrast, dating with the birth-death prior yielded Neogene diversifications, and four rate shifts, one for each of the four richest genera. Nonetheless, dating with the two priors provided similar age estimates for the divergence of cycads from Ginkgo (Carboniferous) and their crown age (Permian). Of these, Bayes factors clearly supported the birth-death prior.ConclusionsThese results suggest the choice of the branching process prior can have a drastic influence on our understanding of evolutionary radiations. Therefore, all dating analyses must involve a model selection process using Bayes factors to select between a Yule or birth-death prior, in particular on ancient clades with a potential pattern of high extinction. We also provide new insights into the history of cycad diversification because we found (i) periods of extinction along the long branches of the genera consistent with fossil data, and (ii) high diversification rates within the Miocene genus radiations.

BackgroundDrosophila melanogaster often shows correlations between latitude and phenotypic or genetic variation on different continents, which suggests local adaptation with respect to a heterogeneous environment. Previous phenotypic analyses of latitudinal clines have investigated mainly physiological, morphological, or life-history traits. Here, we studied latitudinal variation in sleep in D. melanogaster populations from North and Central America. In parallel, we used RNA-seq to identify interpopulation gene expression differences.ResultsWe found that in D. melanogaster the average nighttime sleep bout duration exhibits a latitudinal cline such that sleep bouts of equatorial populations are roughly twice as long as those of temperate populations. Interestingly, this pattern of latitudinal variation is not observed for any daytime measure of activity or sleep. We also found evidence for geographic variation for sunrise anticipation. Our RNA-seq experiment carried out on heads from a low and high latitude population identified a large number of gene expression differences, most of which were time dependent. Differentially expressed genes were enriched in circadian regulated genes and enriched in genes potentially under spatially varying selection.ConclusionOur results are consistent with a mechanistic and selective decoupling of nighttime and daytime activity. Furthermore, the present study suggests that natural selection plays a major role in generating transcriptomic variation associated with circadian behaviors. Finally, we identified genomic variants plausibly causally associated with the observed behavioral and transcriptomic variation.