Abstract: Development of phylogenetic methods that do not rely on fossils for the study of evolutionary processes through time have revolutionized the field of evolutionary biology and resulted in an unprecedented expansion of our knowledge about the tree of life. These methods have helped to shed light on the macroevolution of many taxonomic groups such as the placentals (Mammalia). However, despite the increase of studies addressing the diversification patterns of organisms, no synthesis has addressed the case of the most diversified mammalian clade: the Rodentia. Here we present a rodent maximum likelihood phylogeny inferred from a molecular supermatrix. It is based on 11 mitochondrial and nuclear genes that covers 1,265 species, i.e., respectively 56% and 81% of the known specific and generic rodent diversity. The inferred topology recovered all Rodentia clades proposed by recent molecular works. A relaxed molecular clock dating approach provided a time framework for speciation events. We found that the Myomorpha clade shows a greater degree of variation in diversification rates than Sciuroidea, Caviomorpha, Castorimorpha and Anomaluromorpha. We identified a number of shifts in diversification rates within the major clades: two in Castorimorpha, three in Ctenohystrica, 6 within the squirrel-related clade and 24 in the Myomorpha clade. The majority of these shifts occurred within the most recent familial rodent radiations: the Cricetidae and Muridae clades. Using the topological imbalances and the time line we discuss the potential role of different diversification factors that might have shaped the rodents radiation. The present glimpse on the diversification pattern of rodents can be used for further comparative meta-analyses. Muroid lineages have a greater degree of variation in their diversification rates than any other rodent group. Different topological signatures suggest distinct diversification processes among rodent lineages. In particular, Muroidea and Sciuroidea display widespread distribution and have undergone evolutionary and adaptive radiation on most of the continents. Our results show that rodents experienced shifts in diversification rate regularly through the Tertiary, but at different periods for each clade. A comparison between the rodent fossil record and our results suggest that extinction led to the loss of diversification signal for most of the Paleogene nodes.

Abstract: Background The large Glycoside Hydrolase family 5 (GH5) groups together a wide range of enzymes acting on β-linked oligo- and polysaccharides, and glycoconjugates from a large spectrum of organisms. The long and complex evolution of this family of enzymes and its broad sequence diversity limits functional prediction. With the objective of improving the differentiation of enzyme specificities in a knowledge-based context, and to obtain new evolutionary insights, we present here a new, robust subfamily classification of family GH5.

Abstract: 45S rDNA sites are the most widely documented chromosomal regions in eukaryotes The analysis of the distribution of these sites along the chromosome in several genera has suggested some bias in their distribution In order to evaluate if these loci are in fact non-randomly distributed and what is the influence of some chromosomal and karyotypic features on the distribution of these sites, a database was built with the position and number of 45S rDNA sites obtained by FISH together with other karyotypic data from 846 plant species In angiosperms the most frequent numbers of sites per diploid karyotype were two and four, suggesting that in spite of the wide dispersion capacity of these sequences the number of rDNA sites tends to be restricted The sites showed a preferential distribution on the short arms, mainly in the terminal regions Curiously, these sites were frequently found on the short arms of acrocentric chromosomes where they usually occupy the whole arm The trend to occupy the terminal region is especially evident in holokinetic chromosomes, where all of them were terminally located In polyploids there is a trend towards reduction in the number of sites per monoploid complement In gymnosperms, however, the distribution of rDNA sites varied strongly among the sampled families The location of 45S rDNA sites do not vary randomly, occurring preferentially on the short arm and in the terminal region of chromosomes in angiosperms The meaning of this preferential location is not known, but some hypotheses are considered and the observed trends are discussed

Abstract: Cave animals converge evolutionarily on a suite of troglomorphic traits, the best known of which are eyelessness and depigmentation. We studied 11 cave and 10 surface populations of Astyanax mexicanus in order to better understand the evolutionary origins of the cave forms, the basic genetic structuring of both cave and surface populations, and the degree to which present day migration among them affects their genetic divergence. To assess the genetic structure within populations and the relationships among them we genotyped individuals at 26 microsatellite loci. We found that surface populations are similar to one another, despite their relatively large geographic separation, whereas the cave populations are better differentiated. The cave populations we studied span the full range of the cave forms in three separate geographic regions and have at least five separate evolutionary origins. Cave populations had lower genetic diversity than surface populations, correlated with their smaller effective population sizes, probably the result of food and space limitations. Some of the cave populations receive migrants from the surface and exchange migrants with one another, especially when geographically close. This admixture results in significant heterozygote deficiencies at numerous loci due to Wahlund effects. Cave populations receiving migrants from the surface contain small numbers of individuals that are intermediate in both phenotype and genotype, affirming at least limited gene flow from the surface. Cave populations of this species are derived from two different surface stocks denoted "old" and "new." The old stock colonized caves at least three times independently while the new stock colonized caves at least twice independently. Thus, the similar cave phenotypes found in these caves are the result of repeated convergences. These phenotypic convergences have occurred in spite of gene flow from surface populations suggesting either strong natural or sexual selection for alleles responsible for the cave phenotype in the cave environment.

Abstract: Proteins of the mammalian PYHIN (IFI200/HIN-200) family are involved in defence against infection through recognition of foreign DNA. The family member absent in melanoma 2 (AIM2) binds cytosolic DNA via its HIN domain and initiates inflammasome formation via its pyrin domain. AIM2 lies within a cluster of related genes, many of which are uncharacterised in mouse. To better understand the evolution, orthology and function of these genes, we have documented the range of PYHIN genes present in representative mammalian species, and undertaken phylogenetic and expression analyses. No PYHIN genes are evident in non-mammals or monotremes, with a single member found in each of three marsupial genomes. Placental mammals show variable family expansions, from one gene in cow to four in human and 14 in mouse. A single HIN domain appears to have evolved in the common ancestor of marsupials and placental mammals, and duplicated to give rise to three distinct forms (HIN-A, -B and -C) in the placental mammal ancestor. Phylogenetic analyses showed that AIM2 HIN-C and pyrin domains clearly diverge from the rest of the family, and it is the only PYHIN protein with orthology across many species. Interestingly, although AIM2 is important in defence against some bacteria and viruses in mice, AIM2 is a pseudogene in cow, sheep, llama, dolphin, dog and elephant. The other 13 mouse genes have arisen by duplication and rearrangement within the lineage, which has allowed some diversification in expression patterns. The role of AIM2 in forming the inflammasome is relatively well understood, but molecular interactions of other PYHIN proteins involved in defence against foreign DNA remain to be defined. The non-AIM2 PYHIN protein sequences are very distinct from AIM2, suggesting they vary in effector mechanism in response to foreign DNA, and may bind different DNA structures. The PYHIN family has highly varied gene composition between mammalian species due to lineage-specific duplication and loss, which probably indicates different adaptations for fighting infectious disease. Non-genomic DNA can indicate infection, or a mutagenic threat. We hypothesise that defence of the genome against endogenous retroelements has been an additional evolutionary driver for PYHIN proteins.

Abstract: Tribe Fabeae comprises about 380 legume species, including some of the most ancient and important crops like lentil, pea, and broad bean. Breeding efforts in legume crops rely on a detailed knowledge of closest wild relatives and geographic origin. Relationships within the tribe, however, are incompletely known and previous molecular results conflicted with the traditional morphology-based classification. Here we analyse the systematics, biogeography, and character evolution in the tribe based on plastid and nuclear DNA sequences. Phylogenetic analyses including c. 70% of the species in the tribe show that the genera Vicia and Lathyrus in their current circumscription are not monophyletic: Pisum and Vavilovia are nested in Lathyrus, the genus Lens is nested in Vicia. A small, well-supported clade including Vicia hirsuta, V. sylvatica, and some Mediterranean endemics, is the sister group to all remaining species in the tribe. Fabeae originated in the East Mediterranean region in the Miocene (23–16 million years ago (Ma)) and spread at least 39 times into Eurasia, seven times to the Americas, twice to tropical Africa and four times to Macaronesia. Broad bean (V. faba) and its sister V. paucijuga originated in Asia and might be sister to V. oroboides. Lentil (Lens culinaris ssp. culinaris) is of Mediterranean origin and together with eight very close relatives forms a clade that is nested in the core Vicia, where it evolved c. 14 Ma. The Pisum clade is nested in Lathyrus in a grade with the Mediterranean L. gloeosperma, L. neurolobus, and L. nissolia. The extinct Azorean endemic V. dennesiana belongs in section Cracca and is nested among Mediterranean species. According to our ancestral character state reconstruction results, ancestors of Fabeae had a basic chromosome number of 2n=14, an annual life form, and evenly hairy, dorsiventrally compressed styles. Fabeae evolved in the Eastern Mediterranean in the middle Miocene and spread from there across Eurasia, into Tropical Africa, and at least seven times to the Americas. The middle-Atlantic islands were colonized four times but apparently did not serve as stepping-stones for Atlantic crossings. Long-distance dispersal events are relatively common in Fabeae (seven per ten million years). Current generic and infrageneric circumscriptions in Fabeae do not reflect monophyletic groups and should be revised. Suggestions for generic level delimitation are offered.

Abstract: Geographic variation in the thermal environment impacts a broad range of biochemical and physiological processes and can be a major selective force leading to local population adaptation. In the intertidal copepod Tigriopus californicus, populations along the coast of California show differences in thermal tolerance that are consistent with adaptation, i.e., southern populations withstand thermal stresses that are lethal to northern populations. To understand the genetic basis of these physiological differences, we use an RNA-seq approach to compare genome-wide patterns of gene expression in two populations known to differ in thermal tolerance. Observed differences in gene expression between the southern (San Diego) and the northern (Santa Cruz) populations included both the number of affected loci as well as the identity of these loci. However, the most pronounced differences concerned the amplitude of up-regulation of genes producing heat shock proteins (Hsps) and genes involved in ubiquitination and proteolysis. Among the hsp genes, orthologous pairs show markedly different thermal responses as the amplitude of hsp response was greatly elevated in the San Diego population, most notably in members of the hsp70 gene family. There was no evidence of accelerated evolution at the sequence level for hsp genes. Among other sets of genes, cuticle genes were up-regulated in SD but down-regulated in SC, and mitochondrial genes were down-regulated in both populations. Marked changes in gene expression were observed in response to acute sub-lethal thermal stress in the copepod T. californicus. Although some qualitative differences were observed between populations, the most pronounced differences involved the magnitude of induction of numerous hsp and ubiquitin genes. These differences in gene expression suggest that evolutionary divergence in the regulatory pathway(s) involved in acute temperature stress may offer at least a partial explanation of population differences in thermal tolerance observed in Tigriopus.

Abstract: Hybridization, genetic mixture of distinct populations, gives rise to myriad recombinant genotypes. Characterizing the genomic composition of hybrids is critical for studies of hybrid zone dynamics, inheritance of traits, and consequences of hybridization for evolution and conservation. Hybrid genomes are often summarized either by an estimate of the proportion of alleles coming from each ancestral population or classification into discrete categories like F1, F2, backcross, or merely “hybrid” vs. “pure”. In most cases, it is not realistic to classify individuals into the restricted set of classes produced in the first two generations of admixture. However, the continuous ancestry index misses an important dimension of the genotype. Joint consideration of ancestry together with interclass heterozygosity (proportion of loci with alleles from both ancestral populations) captures all of the information in the discrete classification without the unrealistic assumption that only two generations of admixture have transpired. I describe a maximum likelihood method for joint estimation of ancestry and interclass heterozygosity. I present two worked examples illustrating the value of the approach for describing variation among hybrid populations and evaluating the validity of the assumption underlying discrete classification. Naively classifying natural hybrids into the standard six line cross categories can be misleading, and false classification can be a serious problem for datasets with few molecular markers. My analysis underscores previous work showing that many (50 or more) ancestry informative markers are needed to avoid erroneous classification. Although classification of hybrids might often be misleading, valuable inferences can be obtained by focusing directly on distributions of ancestry and heterozygosity. Estimating and visualizing the joint distribution of ancestry and interclass heterozygosity is an effective way to compare the genetic structure of hybrid populations and these estimates can be used in classic quantitative genetic methods for assessing additive, dominant, and epistatic genetic effects on hybrid phenotypes and fitness. The methods are implemented in a freely available package “HIest” for the R statistical software ( http://cran.r-project.org/web/packages/HIest/index.html ).

Abstract: Background Many marine meiofaunal species are reported to have wide distributions, which creates a paradox considering their hypothesized low dispersal abilities. Correlated with this paradox is an especially high taxonomic deficit for meiofauna, partly related to a lower taxonomic effort and partly to a high number of putative cryptic species. Molecular-based species delineation and barcoding approaches have been advocated for meiofaunal biodiversity assessments to speed up description processes and uncover cryptic lineages. However, these approaches show sensitivity to sampling coverage (taxonomic and geographic) and the success rate has never been explored on mesopsammic Mollusca.

Abstract: Adaptive radiation is the process by which a single ancestral species diversifies into many descendants adapted to exploit a wide range of habitats. The appearance of ecological opportunities, or the colonisation or adaptation to novel ecological resources, has been documented to promote adaptive radiation in many classic examples. Mutualistic interactions allow species to access resources untapped by competitors, but evidence shows that the effect of mutualism on species diversification can greatly vary among mutualistic systems. Here, we test whether the development of obligate mutualism with sea anemones allowed the clownfishes to radiate adaptively across the Indian and western Pacific oceans reef habitats. We show that clownfishes morphological characters are linked with ecological niches associated with the sea anemones. This pattern is consistent with the ecological speciation hypothesis. Furthermore, the clownfishes show an increase in the rate of species diversification as well as rate of morphological evolution compared to their closest relatives without anemone mutualistic associations. The effect of mutualism on species diversification has only been studied in a limited number of groups. We present a case of adaptive radiation where mutualistic interaction is the likely key innovation, providing new insights into the mechanisms involved in the buildup of biodiversity. Due to a lack of barriers to dispersal, ecological speciation is rare in marine environments. Particular life-history characteristics of clownfishes likely reinforced reproductive isolation between populations, allowing rapid species diversification.

Abstract: The APOBEC3 (A3) genes play a key role in innate antiviral defense in mammals by introducing directed mutations in the DNA. The human genome encodes for seven A3 genes, with multiple splice alternatives. Different A3 proteins display different substrate specificity, but the very basic question on how discerning self from non-self still remains unresolved. Further, the expression of A3 activity/ies shapes the way both viral and host genomes evolve. We present here a detailed temporal analysis of the origin and expansion of the A3 repertoire in mammals. Our data support an evolutionary scenario where the genome of the mammalian ancestor encoded for at least one ancestral A3 gene, and where the genome of the ancestor of placental mammals (and possibly of the ancestor of all mammals) already encoded for an A3Z1-A3Z2-A3Z3 arrangement. Duplication events of the A3 genes have occurred independently in different lineages: humans, cats and horses. In all of them, gene duplication has resulted in changes in enzyme activity and/or substrate specificity, in a paradigmatic example of convergent adaptive evolution at the genomic level. Finally, our results show that evolutionary rates for the three A3Z1, A3Z2 and A3Z3 motifs have significantly decreased in the last 100 Mya. The analysis constitutes a textbook example of the evolution of a gene locus by duplication and sub/neofunctionalization in the context of virus-host arms race. Our results provide a time framework for identifying ancestral and derived genomic arrangements in the APOBEC loci, and to date the expansion of this gene family for different lineages through time, as a response to changes in viral/retroviral/retrotransposon pressure.

Abstract: Rates of molecular evolution are known to vary across taxa and among genes, and this requires rate calibration for each specific dataset based on external information. Calibration is sensitive to evolutionary model parameters, partitioning schemes and clock model. However, the way in which these and other analytical aspects affect both the rates and the resulting clade ages from calibrated phylogenies are not yet well understood. To investigate these aspects we have conducted calibration analyses for the genus Carabus (Coleoptera, Carabidae) on five mitochondrial and four nuclear DNA fragments with 7888 nt total length, testing different clock models and partitioning schemes to select the most suitable using Bayes Factors comparisons. We used these data to investigate the effect of ambiguous character and outgroup inclusion on both the rates of molecular evolution and the TMRCA of Carabus. We found considerable variation in rates of molecular evolution depending on the fragment studied (ranging from 5.02% in cob to 0.26% divergence/My in LSU-A), but also on analytical conditions. Alternative choices of clock model, partitioning scheme, treatment of ambiguous characters, and outgroup inclusion resulted in rate increments ranging from 28% (HUWE1) to 1000% (LSU-B and ITS2) and increments in the TMRCA of Carabus ranging from 8.4% (cox1-A) to 540% (ITS2). Results support an origin of the genus Carabus during the Oligocene in the Eurasian continent followed by a Miocene differentiation that originated all main extant lineages. The combination of several genes is proposed as the best strategy to minimise both the idiosyncratic behaviors of individual markers and the effect of analytical aspects in rate and age estimations. Our results highlight the importance of estimating rates of molecular evolution for each specific dataset, selecting for optimal clock and partitioning models as well as other methodological issues potentially affecting rate estimation.

Abstract: The discovery of giant viruses with genome and physical size comparable to cellular organisms, remnants of protein translation machinery and virus-specific parasites (virophages) have raised intriguing questions about their origin. Evidence advocates for their inclusion into global phylogenomic studies and their consideration as a distinct and ancient form of life. Here we reconstruct phylogenies describing the evolution of proteomes and protein domain structures of cellular organisms and double-stranded DNA viruses with medium-to-very-large proteomes (giant viruses). Trees of proteomes define viruses as a ‘fourth supergroup’ along with superkingdoms Archaea, Bacteria, and Eukarya. Trees of domains indicate they have evolved via massive and primordial reductive evolutionary processes. The distribution of domain structures suggests giant viruses harbor a significant number of protein domains including those with no cellular representation. The genomic and structural diversity embedded in the viral proteomes is comparable to the cellular proteomes of organisms with parasitic lifestyles. Since viral domains are widespread among cellular species, we propose that viruses mediate gene transfer between cells and crucially enhance biodiversity. Results call for a change in the way viruses are perceived. They likely represent a distinct form of life that either predated or coexisted with the last universal common ancestor (LUCA) and constitute a very crucial part of our planet’s biosphere.

Abstract: Fungus farming is an unusual life style in insects that has evolved many times in the wood boring weevils named ‘ambrosia beetles’. Multiple occurrences of this behaviour allow for a detailed comparison of the different origins of fungus farming through time, its directionality, and possible ancestral states. We tested these hypotheses with a phylogeny representing the largest data set to date, nearly 4 kb of nucleotides from COI, EF-1α, CAD, ArgK, 28S, and 200 scolytine taxa. Phylogenetic analyses using Bayesian or parsimony approaches placed the root of Scolytinae close to the tribe Scolytini and Microborus, but otherwise indicated low resolution at older nodes. More recent clades were well resolved, including ten origins of fungus farming. There were no subsequent reversals to bark or phloem feeding in the fungus farming clades. The oldest origin of fungus farming was estimated near 50 Ma, long after the origin of Scolytinae (100-120 Ma). Younger origins included the species rich Xyleborini, dated to 21 Ma. Sister group comparisons and test of independence between traits indicated that neither gregarious larval feeding nor regular inbreeding by sibling mating was strongly correlated with the origin of fungus farming. Origins of fungus farming corresponded mainly with two periods of global warming in the Cenozoic era, which were characterised by broadly distributed tropical forests. Hence, it seems likely that warm climates and expanding tropical angiosperm forests played critical roles in the successful radiation of diverse fungus farming groups. However, further investigation will likely reveal additional biological factors that promote fungus farming.

Abstract: Uncertainty in comparative analyses can come from at least two sources: a) phylogenetic uncertainty in the tree topology or branch lengths, and b) uncertainty due to intraspecific variation in trait values, either due to measurement error or natural individual variation. Most phylogenetic comparative methods do not account for such uncertainties. Not accounting for these sources of uncertainty leads to false perceptions of precision (confidence intervals will be too narrow) and inflated significance in hypothesis testing (e.g. p-values will be too small). Although there is some application-specific software for fitting Bayesian models accounting for phylogenetic error, more general and flexible software is desirable. We developed models to directly incorporate phylogenetic uncertainty into a range of analyses that biologists commonly perform, using a Bayesian framework and Markov Chain Monte Carlo analyses. We demonstrate applications in linear regression, quantification of phylogenetic signal, and measurement error models. Phylogenetic uncertainty was incorporated by applying a prior distribution for the phylogeny, where this distribution consisted of the posterior tree sets from Bayesian phylogenetic tree estimation programs. The models were analysed using simulated data sets, and applied to a real data set on plant traits, from rainforest plant species in Northern Australia. Analyses were performed using the free and open source software OpenBUGS and JAGS. Incorporating phylogenetic uncertainty through an empirical prior distribution of trees leads to more precise estimation of regression model parameters than using a single consensus tree and enables a more realistic estimation of confidence intervals. In addition, models incorporating measurement errors and/or individual variation, in one or both variables, are easily formulated in the Bayesian framework. We show that BUGS is a useful, flexible general purpose tool for phylogenetic comparative analyses, particularly for modelling in the face of phylogenetic uncertainty and accounting for measurement error or individual variation in explanatory variables. Code for all models is provided in the BUGS model description language.

Abstract: The evolutionary history of the biota of North Africa and Arabia is inextricably tied to the complex geological and climatic evolution that gave rise to the prevalent deserts of these areas. Reptiles constitute an exemplary group in the study of the arid environments with numerous well-adapted members, while recent studies using reptiles as models have unveiled interesting biogeographical and diversification patterns. In this study, we include 207 specimens belonging to all 12 recognized species of the genus Stenodactylus. Molecular phylogenies inferred using two mitochondrial (12S rRNA and 16S rRNA) and two nuclear (c-mos and RAG-2) markers are employed to obtain a robust time-calibrated phylogeny, as the base to investigate the inter- and intraspecific relationships and to elucidate the biogeographical history of Stenodactylus, a genus with a large distribution range including the arid and hyper-arid areas of North Africa and Arabia. The phylogenetic analyses of molecular data reveal the existence of three major clades within the genus Stenodactylus, which is supported by previous studies based on morphology. Estimated divergence times between clades and sub-clades are shown to correlate with major geological events of the region, the most important of which is the opening of the Red Sea, while climatic instability in the Miocene is hypothesized to have triggered diversification. High genetic variability is observed in some species, suggesting the existence of some undescribed species. The S. petrii - S. stenurus species complex is in need of a thorough taxonomic revision. New data is presented on the distribution of the sister species S. sthenodactylus and S. mauritanicus. The phylogenetic hypothesis for the genus Stenodactylus presented in this work permits the reconstruction of the biogeographical history of these common desert dwellers and confirms the importance of the opening of the Red Sea and the climatic oscillations of the Miocene as major factors in the diversification of the biota of North Africa and Arabia. Moreover, this study traces the evolution of this widely distributed and highly specialized group, investigates the patterns of its high intraspecific diversity and elucidates its systematics.

Abstract: The critical role of Major Histocompatibility Complex (Mhc) genes in disease resistance and their highly polymorphic nature make them exceptional candidates for studies investigating genetic effects on survival, mate choice and conservation. Species that harbor many Mhc loci and high allelic diversity are particularly intriguing as they are potentially under strong selection and studies of such species provide valuable information as to the mechanisms maintaining Mhc diversity. However comprehensive genotyping of complex multilocus systems has been a major challenge to date with the result that little is known about the consequences of this complexity in terms of fitness effects and disease resistance. In this study, we genotyped the Mhc class I exon 3 of the great tit (Parus major) from two nest-box breeding populations near Oxford, UK that have been monitored for decades. Characterization of Mhc class I exon 3 was adopted and bidirectional sequencing was carried using the 454 sequencing platform. Full analysis of sequences through a stepwise variant validation procedure allowed reliable typing of more than 800 great tits based on 214,357 reads; from duplicates we estimated the repeatability of typing as 0.94. A total of 862 alleles were detected, and the presence of at least 16 functional loci was shown - the highest number characterized in a wild bird species. Finally, the functional alleles were grouped into 17 supertypes based on their antigen binding affinities. We found extreme complexity at the Mhc class I of the great tit both in terms of allelic diversity and gene number. The presence of many functional loci was shown, together with a pseudogene family and putatively non-functional alleles; there was clear evidence that functional alleles were under strong balancing selection. This study is the first step towards an in-depth analysis of this gene complex in this species, which will help understanding how parasite-mediated and sexual selection shape and maintain host genetic variation in nature. We believe that study systems like ours can make important contributions to the field of evolutionary biology and emphasize the necessity of integrating long-term field-based studies with detailed genetic analysis to unravel complex evolutionary processes.

Abstract: Background: The Ampelopsis clade (Ampelopsis and its close allies) of the grape family Vitaceae contains ca. 43 species disjunctly distributed in Asia, Europe, North America, South America, Africa, and Australia, and is a rare example to study both the Northern and the Southern Hemisphere intercontinental disjunctions. We reconstruct the temporal and spatial diversification of the Ampelopsis clade to explore the evolutionary processes that have resulted in their intercontinental disjunctions in six continents. Results: The Bayesian molecular clock dating and the likelihood ancestral area analyses suggest that the Ampelopsis clade most likely originated in North America with its crown group dated at 41.2 Ma (95% HPD 23.4 - 61.0 Ma) in the middle Eocene. Two independent Laurasian migrations into Eurasia are inferred to have occurred in the early Miocene via the North Atlantic land bridges. The ancestor of the Southern Hemisphere lineage migrated from North America to South America in the early Oligocene. The Gondwanan-like pattern of intercontinental disjunction is best explained by two long-distance dispersals: once from South America to Africa estimated at 30.5 Ma (95% HPD 16.9 - 45.9 Ma), and the other from South America to Australia dated to 19.2 Ma (95% HPD 6.7 - 22.3 Ma). Conclusions: The global disjunctions in the Ampelopsis clade are best explained by a diversification model of North American origin, two Laurasian migrations, one migration into South America, and two post-Gondwanan long-distance dispersals. These findings highlight the importance of both vicariance and long distance dispersal in shaping intercontinental disjunctions of flowering plants.

Abstract: Double-stranded (ds) RNA fungal viruses are typically isometric single-shelled particles that are classified into three families, Totiviridae, Partitiviridae and Chrysoviridae, the members of which possess monopartite, bipartite and quadripartite genomes, respectively Recent findings revealed that mycovirus-related dsRNA viruses are more diverse than previously recognized Although an increasing number of viral complete genomic sequences have become available, the evolution of these diverse dsRNA viruses remains to be clarified This is particularly so since there is little evidence for horizontal gene transfer (HGT) among dsRNA viruses In this study, we report the molecular properties of two novel dsRNA mycoviruses that were isolated from a field strain of Sclerotinia sclerotiorum, Sunf-M: one is a large monopartite virus representing a distinct evolutionary lineage of dsRNA viruses; the other is a new member of the family Partitiviridae Comprehensive phylogenetic analysis and genome comparison revealed that there are at least ten monopartite, three bipartite, one tripartite and three quadripartite lineages in the known dsRNA mycoviruses and that the multipartite lineages have possibly evolved from different monopartite dsRNA viruses Moreover, we found that homologs of the S7 Domain, characteristic of members of the genus phytoreovirus in family Reoviridae are widely distributed in diverse dsRNA viral lineages, including chrysoviruses, endornaviruses and some unclassified dsRNA mycoviruses We further provided evidence that multiple HGT events may have occurred among these dsRNA viruses from different families Our study provides an insight into the phylogeny and evolution of mycovirus-related dsRNA viruses and reveals that the occurrence of HGT between different virus species and the development of multipartite genomes during evolution are important macroevolutionary mechanisms in dsRNA viruses

Abstract: Estimates of relationships among Staphylococcus species have been hampered by poor and inconsistent resolution of phylogenies based largely on single gene analyses incorporating only a limited taxon sample. As such, the evolutionary relationships and hierarchical classification schemes among species have not been confidently established. Here, we address these points through analyses of DNA sequence data from multiple loci (16S rRNA gene, dnaJ, rpoB, and tuf gene fragments) using multiple Bayesian and maximum likelihood phylogenetic approaches that incorporate nearly all recognized Staphylococcus taxa. We estimated the phylogeny of fifty-seven Staphylococcus taxa using partitioned-model Bayesian and maximum likelihood analysis, as well as Bayesian gene-tree species-tree methods. Regardless of methodology, we found broad agreement among methods that the current cluster groups require revision, although there was some disagreement among methods in resolution of higher order relationships. Based on our phylogenetic estimates, we propose a refined classification for Staphylococcus with species being classified into 15 cluster groups (based on molecular data) that adhere to six species groups (based on phenotypic properties). Our findings are in general agreement with gene tree-based reports of the staphylococcal phylogeny, although we identify multiple previously unreported relationships among species. Our results support the general importance of such multilocus assessments as a standard in microbial studies to more robustly infer relationships among recognized and newly discovered lineages.

Abstract: Leanchoilia superlata is one of the best known arthropods from the middle Cambrian Burgess Shale of British Columbia Here we re-describe the morphology of L superlata and discuss its possible autecology The re-description follows a standardized scheme, the descriptive matrix approach, designed to provide a template for descriptions of other megacheiran species Our findings differ in several respects from previous interpretations Examples include a more slender body; a possible hypostome; a small specialised second appendage, bringing the number of pairs of head appendages to four; a further sub-division of the great appendage, making it more similar to that of other megacheirans; and a complex joint of the exopod reflecting the arthropod’s swimming capabilities Different aspects of the morphology, for example, the morphology of the great appendage and the presence of a basipod with strong median armature on the biramous appendages indicate that L superlata was an active and agile necto-benthic predator (not a scavenger or deposit feeder as previously interpreted)

Abstract: Bactrocera dorsalis s.s. is a pestiferous tephritid fruit fly distributed from Pakistan to the Pacific, with the Thai/Malay peninsula its southern limit. Sister pest taxa, B. papayae and B. philippinensis, occur in the southeast Asian archipelago and the Philippines, respectively. The relationship among these species is unclear due to their high molecular and morphological similarity. This study analysed population structure of these three species within a southeast Asian biogeographical context to assess potential dispersal patterns and the validity of their current taxonomic status. Geometric morphometric results generated from 15 landmarks for wings of 169 flies revealed significant differences in wing shape between almost all sites following canonical variate analysis. For the combined data set there was a greater isolation-by-distance (IBD) effect under a ‘non-Euclidean’ scenario which used geographical distances within a biogeographical ‘Sundaland context’ (r 2 = 0.772, P < 0.0001) as compared to a ‘Euclidean’ scenario for which direct geographic distances between sample sites was used (r 2 = 0.217, P < 0.01). COI sequence data were obtained for 156 individuals and yielded 83 unique haplotypes with no correlation to current taxonomic designations via a minimum spanning network. beast analysis provided a root age and location of 540kya in northern Thailand, with migration of B. dorsalis s.l. into Malaysia 470kya and Sumatra 270kya. Two migration events into the Philippines are inferred. Sequence data revealed a weak but significant IBD effect under the ‘non-Euclidean’ scenario (r 2 = 0.110, P < 0.05), with no historical migration evident between Taiwan and the Philippines. Results are consistent with those expected at the intra-specific level. Bactrocera dorsalis s.s., B. papayae and B. philippinensis likely represent one species structured around the South China Sea, having migrated from northern Thailand into the southeast Asian archipelago and across into the Philippines. No migration is apparent between the Philippines and Taiwan. This information has implications for quarantine, trade and pest management.

Abstract: Low-complexity regions (LCRs) in proteins are tracts that are highly enriched in one or a few amino acids. Given their high abundance, and their capacity to expand in relatively short periods of time through replication slippage, they can greatly contribute to increase protein sequence space and generate novel protein functions. However, little is known about the global impact of LCRs on protein evolution. We have traced back the evolutionary history of 2,802 LCRs from a large set of homologous protein families from H.sapiens, M.musculus, G.gallus, D.rerio and C.intestinalis. Transcriptional factors and other regulatory functions are overrepresented in proteins containing LCRs. We have found that the gain of novel LCRs is frequently associated with repeat expansion whereas the loss of LCRs is more often due to accumulation of amino acid substitutions as opposed to deletions. This dichotomy results in net protein sequence gain over time. We have detected a significant increase in the rate of accumulation of novel LCRs in the ancestral Amniota and mammalian branches, and a reduction in the chicken branch. Alanine and/or glycine-rich LCRs are overrepresented in recently emerged LCR sets from all branches, suggesting that their expansion is better tolerated than for other LCR types. LCRs enriched in positively charged amino acids show the contrary pattern, indicating an important effect of purifying selection in their maintenance. We have performed the first large-scale study on the evolutionary dynamics of LCRs in protein families. The study has shown that the composition of an LCR is an important determinant of its evolutionary pattern.

Abstract: Hydrocharitaceae is a fully aquatic monocot family, consists of 18 genera with approximately 120 species. The family includes both fresh and marine aquatics and exhibits great diversity in form and habit including annual and perennial life histories; submersed, partially submersed and floating leaf habits and linear to orbicular leaf shapes. The family has a cosmopolitan distribution and is well represented in the Tertiary fossil record in Europe. At present, the historical biogeography of the family is not well understood and the generic relationships remain controversial. In this study we investigated the phylogeny and biogeography of Hydrocharitaceae by integrating fossils and DNA sequences from eight genes. We also conducted ancestral state reconstruction for three morphological characters. Phylogenetic analyses produced a phylogeny with most branches strongly supported by bootstrap values greater than 95 and Bayesian posterior probability values of 1.0. Stratiotes is the first diverging lineage with the remaining genera in two clades, one clade consists of Lagarosiphon, Ottelia, Blyxa, Apalanthe, Elodea and Egeria; and the other consists of Hydrocharis-Limnobium, Thalassia, Enhalus, Halophila, Najas, Hydrilla, Vallisneria, Nechamandra and Maidenia. Biogeographic analyses (DIVA, Mesquite) and divergence time estimates (BEAST) resolved the most recent common ancestor of Hydrocharitaceae as being in Asia during the Late Cretaceous and Palaeocene (54.7-72.6 Ma). Dispersals (including long-distance dispersal and migrations through Tethys seaway and land bridges) probably played major roles in the intercontinental distribution of this family. Ancestral state reconstruction suggested that in Hydrocharitaceae evolution of dioecy is bidirectional, viz., from dioecy to hermaphroditism, and from hermaphroditism to dioecy, and that the aerial-submerged leaf habit and short-linear leaf shape are the ancestral states. Our study has shed light on the previously controversial generic phylogeny of Hydrocharitaceae. The study has resolved the historical biogeography of this family and supported dispersal as the most likely explanation for the intercontinental distribution. We have also provided valuable information for understanding the evolution of breeding system and leaf phenotype in aquatic monocots.

Abstract: Background: The glucosinolate-myrosinase system is an activated chemical defense system found in plants of the Brassicales order. Glucosinolates are stored separately from their hydrolytic enzymes, the myrosinases, in plant tissues. Upon tissue damage, e.g. by herbivory, glucosinolates and myrosinases get mixed and glucosinolates are broken down to an array of biologically active compounds of which isothiocyanates are toxic to a wide range of organisms. Specifier proteins occur in some, but not all glucosinolate-containing plants and promote the formation of biologically active non-isothiocyanate products upon myrosinase-catalyzed glucosinolate breakdown. Results: Based on a phytochemical screening among representatives of the Brassicales order, we selected candidate species for identification of specifier protein cDNAs. We identified ten specifier proteins from a range of species of the Brassicaceae and assigned each of them to one of the three specifier protein types (NSP, nitrile-specifier protein, ESP, epithiospecifier protein, TFP, thiocyanate-forming protein) after heterologous expression in Escherichia coli. Together with nine known specifier proteins and three putative specifier proteins found in databases, we subjected the newly identified specifier proteins to phylogenetic analyses. Specifier proteins formed three major clusters, named AtNSP5-cluster, AtNSP1-cluster, and ESP/TFP cluster. Within the ESP/TFP cluster, specifier proteins grouped according to the Brassicaceae lineage they were identified from. Non-synonymous vs. synonymous substitution rate ratios suggested purifying selection to act on specifier protein genes. Conclusions: Among specifier proteins, NSPs represent the ancestral activity. The data support a monophyletic origin of ESPs from NSPs. The split between NSPs and ESPs/TFPs happened before the radiation of the core Brassicaceae. Future analyses have to show if TFP activity evolved from ESPs at least twice independently in different Brassicaceae lineages as suggested by the phylogeny. The ability to form non-isothiocyanate products by specifier protein activity may provide plants with a selective advantage. The evolution of specifier proteins in the Brassicaceae demonstrates the plasticity of secondary metabolism within an activated plant defense system.

Abstract: Background Although evolutionary models of cooperation build on the intuition that costs of the donor and benefits to the receiver are the most general fundamental parameters, it is largely unknown how they affect the decision of animals to cooperate with an unrelated social partner. Here we test experimentally whether costs to the donor and need of the receiver decide about the amount of help provided by unrelated rats in an iterated prisoner's dilemma game.

Abstract: Bacterial genomes exhibit a remarkable degree of variation in the presence and absence of genes, which probably extends to the level of individual pathways. This variation may be a consequence of the significant evolutionary role played by horizontal gene transfer, but might also be explained by the loss of genes through mutation. A challenge is to understand why there would be variation in gene presence within pathways if they confer a benefit only when complete. Here, we develop a mathematical model to study how variation in pathway content is produced by horizontal transfer, gene loss and partial exposure of a population to a novel environment. We discuss the possibility that variation in gene presence acts as cryptic genetic variation on which selection acts when the appropriate environment occurs. We find that a high level of variation in gene presence can be readily explained by decay of the pathway through mutation when there is no longer exposure to the selective environment, or when selection becomes too weak to maintain the genes. In the context of pathway variation the role of horizontal gene transfer is probably the initial introduction of a complete novel pathway rather than in building up the variation in a genome without the pathway.

Abstract: Arbuscular mycorrhizal fungi (AMF) can form obligate symbioses with the vast majority of land plants, and AMF distribution patterns have received increasing attention from researchers. At the local scale, the distribution of AMF is well documented. Studies at large scales, however, are limited because intensive sampling is difficult. Here, we used ITS rDNA sequence metadata obtained from public databases to study the distribution of AMF at continental and global scales. We also used these sequence metadata to investigate whether host plant is the main factor that affects the distribution of AMF at large scales. We defined 305 ITS virtual taxa (ITS-VTs) among all sequences of the Glomeromycota by using a comprehensive maximum likelihood phylogenetic analysis. Each host taxonomic order averaged about 53% specific ITS-VTs, and approximately 60% of the ITS-VTs were host specific. Those ITS-VTs with wide host range showed wide geographic distribution. Most ITS-VTs occurred in only one type of host functional group. The distributions of most ITS-VTs were limited across ecosystem, across continent, across biogeographical realm, and across climatic zone. Non-metric multidimensional scaling analysis (NMDS) showed that AMF community composition differed among functional groups of hosts, and among ecosystem, continent, biogeographical realm, and climatic zone. The Mantel test showed that AMF community composition was significantly correlated with plant community composition among ecosystem, among continent, among biogeographical realm, and among climatic zone. The structural equation modeling (SEM) showed that the effects of ecosystem, continent, biogeographical realm, and climatic zone were mainly indirect on AMF distribution, but plant had strongly direct effects on AMF. The distribution of AMF as indicated by ITS rDNA sequences showed a pattern of high endemism at large scales. This pattern indicates high specificity of AMF for host at different scales (plant taxonomic order and functional group) and high selectivity from host plants for AMF. The effects of ecosystemic, biogeographical, continental and climatic factors on AMF distribution might be mediated by host plants.

Abstract: The benzoxazinoids 2,4-dihydroxy-1,4-benzoxazin-3-one (DIBOA) and 2,4-dihydroxy-7- methoxy-1,4-benzoxazin-3-one (DIMBOA), are key defense compounds present in major agricultural crops such as maize and wheat. Their biosynthesis involves nine enzymes thought to form a linear pathway leading to the storage of DI(M)BOA as glucoside conjugates. Seven of the genes (Bx1-Bx6 and Bx8) form a cluster at the tip of the short arm of maize chromosome 4 that includes four P450 genes (Bx2-5) belonging to the same CYP71C subfamily. The origin of this cluster is unknown. We show that the pathway appeared following several duplications of the TSA gene (α-subunit of tryptophan synthase) and of a Bx2-like ancestral CYP71C gene and the recruitment of Bx8 before the radiation of Poaceae. The origins of Bx6 and Bx7 remain unclear. We demonstrate that the Bx2-like CYP71C ancestor was not committed to the benzoxazinoid pathway and that after duplications the Bx2-Bx5 genes were under positive selection on a few sites and underwent functional divergence, leading to the current specific biochemical properties of the enzymes. The absence of synteny between available Poaceae genomes involving the Bx gene regions is in contrast with the conserved synteny in the TSA gene region. These results demonstrate that rearrangements following duplications of an IGL/TSA gene and of a CYP71C gene probably resulted in the clustering of the new copies (Bx1 and Bx2) at the tip of a chromosome in an ancestor of grasses. Clustering favored cosegregation and tip chromosomal location favored gene rearrangements that allowed the further recruitment of genes to the pathway. These events, a founding event and elongation events, may have been the key to the subsequent evolution of the benzoxazinoid biosynthetic cluster.