Evolutionary footprint of coevolving positions in genes.

MOTIVATION: The analysis of molecular coevolution provides information on the potential functional and structural implication of positions along DNA sequences, and several methods are available to identify coevolving positions using probabilistic or combinatorial approaches. The specific nucleotide or amino acid profile associated with the coevolution process is, however, not estimated, but only known profiles, such as the Watson-Crick constraint, are usually considered a priori in current measures of coevolution. RESULTS: Here, we propose a new probabilistic model, Coev, to identify coevolving positions and their associated profile in DNA sequences while incorporating the underlying phylogenetic relationships. The process of coevolution is modeled by a 16 × 16 instantaneous rate matrix that includes rates of transition as well as a profile of coevolution. We used simulated, empirical and illustrative data to evaluate our model and to compare it with a model of 'independent' evolution using Akaike Information Criterion. We showed that the Coev model is able to discriminate between coevolving and non-coevolving positions and provides better specificity and specificity than other available approaches. We further demonstrate that the identification of the profile of coevolution can shed new light on the process of dependent substitution during lineage evolution.

Multilocus phylogeography of the common lizard Zootoca vivipara at the Ibero-Pyrenean suture zone reveals lowland barriers and high-elevation introgression.

BACKGROUND: The geographic distribution of evolutionary lineages and the patterns of gene flow upon secondary contact provide insight into the process of divergence and speciation. We explore the evolutionary history of the common lizard Zootoca vivipara (= Lacerta vivipara) in the Iberian Peninsula and test the role of the Pyrenees and the Cantabrian Mountains in restricting gene flow and driving lineage isolation and divergence. We also assess patterns of introgression among lineages upon secondary contact, and test for the role of high-elevation trans-mountain colonisations in explaining spatial patterns of genetic diversity. We use mtDNA sequence data and genome-wide AFLP loci to reconstruct phylogenetic relationships among lineages, and measure genetic structure RESULTS: The main genetic split in mtDNA corresponds generally to the French and Spanish sides of the Pyrenees as previously reported, in contrast to genome-wide AFLP data, which show a major division between NW Spain and the rest. Both types of markers support the existence of four distinct and geographically congruent genetic groups, which are consistent with major topographic barriers. Both datasets reveal the presence of three independent contact zones between lineages in the Pyrenean region, one in the Basque lowlands, one in the low-elevation mountains of the western Pyrenees, and one in the French side of the central Pyrenees. The latter shows genetic evidence of a recent, high-altitude trans-Pyrenean incursion from Spain into France. CONCLUSIONS: The distribution and age of major lineages is consistent with a Pleistocene origin and a role for both the Pyrenees and the Cantabrian Mountains in driving isolation and differentiation of Z. vivipara lineages at large geographic scales. However, mountain ranges are not always effective barriers to dispersal, and have not prevented a recent high-elevation trans-Pyrenean incursion that has led to asymmetrical introgression among divergent lineages. Cytonuclear discordance in patterns of genetic structure and introgression at contact zones suggests selection may be involved at various scales. Suture zones are important areas for the study of lineage formation and speciation, and our results show that biogeographic barriers can yield markedly different phylogeographic patterns in different vertebrate and invertebrate taxa.

Comparative transcriptomics of early dipteran development

BACKGROUND: Modern sequencing technologies have massively increased the amount of data available for comparative genomics. Whole-transcriptome shotgun sequencing (RNA-seq) provides a powerful basis for comparative studies. In particular, this approach holds great promise for emerging model species in fields such as evolutionary developmental biology (evo-devo). RESULTS: We have sequenced early embryonic transcriptomes of two non-drosophilid dipteran species: the moth midge Clogmia albipunctata, and the scuttle fly Megaselia abdita. Our analysis includes a third, published, transcriptome for the hoverfly Episyrphus balteatus. These emerging models for comparative developmental studies close an important phylogenetic gap between Drosophila melanogaster and other insect model systems. In this paper, we provide a comparative analysis of early embryonic transcriptomes across species, and use our data for a phylogenomic re-evaluation of dipteran phylogenetic relationships. CONCLUSIONS: We show how comparative transcriptomics can be used to create useful resources for evo-devo, and to investigate phylogenetic relationships. Our results demonstrate that de novo assembly of short (Illumina) reads yields high-quality, high-coverage transcriptomic data sets. We use these data to investigate deep dipteran phylogenetic relationships. Our results, based on a concatenation of 160 orthologous genes, provide support for the traditional view of Clogmia being the sister group of Brachycera (Megaselia, Episyrphus, Drosophila), rather than that of Culicomorpha (which includes mosquitoes and blackflies).

The odd couple: contrasting phylogeographic patterns in two sympatric sibling species of woodlouse-­‐hunter spiders in the Canary Islands

Comparative phylogeography seeks for commonalities in the spatial demographic history of sympatric organisms to characterize the mechanisms that shaped such patterns. The unveiling of incongruent phylogeographic patterns in co-occurring species, on the other hand, may hint to overlooked differences in their life histories or microhabitat preferences. The woodlouse-hunter spiders of the genus Dysdera have undergone a major diversi cation on the Canary Islands. The species pair Dysdera alegranzaensis and Dysdera nesiotes are endemic to the island of Lanzarote and nearby islets, where they co-occur at most of their known localities. The two species stand in sharp contrast to other sympatric endemic Dysdera in showing no evidence of somatic (non-genitalic) differentiation. Phylogenetic and population genetic analyses of mitochondrial cox1 sequences from an exhaustive sample of D. alegranzaensis and D. nesiotes specimens, and additional mitochondrial (16S, L1, nad1) and nuclear genes (28S, H3) were analysed to reveal their phylogeographic patterns and clarify their phylogenetic relationships. Relaxed molecular clock models using ve calibration points were further used to estimate divergence times between species and populations. Striking differences in phylogeography and population structure between the two species were observed. Dysdera nesiotes displayed a metapopulation-like structure, while D. alegranzaensis was characterized by a weaker geographical structure but greater genetic divergences among its main haplotype lineages, suggesting more complex population dynamics. Our study con rms that co-distributed sibling species may exhibit contrasting phylogeographic patterns in the absence of somatic differentiation. Further ecological studies, however, will be necessary to clarify whether the contrasting phylogeographies may hint at an overlooked niche partitioning between the two species. In addition, further comparisons with available phylogeographic data of other eastern Canarian Dysdera endemics con rm the key role of lava ows in structuring local populations in oceanic islands and identify localities that acted as refugia during volcanic eruptions

Evaluation of Xylella fastidiosa genetic diversity by fAFPL markers

The first phytopathogenic bacterium with its DNA entirely sequenced is being detected and isolated from different host plants in several geographic regions. Although it causes diseases in cultures of economic importance, such as citrus, coffee, and grapevine little is known about the genetic relationships among different strains. Actually, all strains are grouped as a single species, Xylella fastidiosa, despite colonizing different hosts, developing symptoms, and different physiological and microbiological observed conditions. The existence of genetic diversity among X. fastidiosa strains was detected by different methodological techniques, since cultural to molecular methods. However, little is know about the phylogenetic relationships developed by Brazilian strains obtained from coffee and citrus plants. In order to evaluate it, fAFLP markers were used to verify genetic diversity and phylogenetic relationships developed by Brazilian and strange strains. fAFLP is an efficient technique, with high reproducibility that is currently used for bacterial typing and classification. The obtained results showed that Brazilian strains present genetic diversity and that the strains from this study were grouped distinctly according host and geographical origin like citrus-coffee, temecula-grapevine-mulberry and plum-elm.

Metal dealing at the origin of the Chordata phylum: The metallothionein system and metal overload response in Amphioxus.

Non-vertebrate chordates, specifically amphioxus, are considered of the utmost interest for gaining insight into the evolutionary trends, i.e. differentiation and specialization, of gene/protein systems. In this work, MTs (metallothioneins), the most important metal binding proteins, are characterized for the first time in the cephalochordate subphylum at both gene and protein level, together with the main features defining the amphioxus response to cadmium and copper overload. Two MT genes (BfMT1 and BfMT2) have been identified in a contiguous region of the genome, as well as several ARE (antioxidant response element) and MRE (metal response element) located upstream the transcribed region. Their corresponding cDNAs exhibit identical sequence in the two lancelet species (B. floridae and B. lanceolatum), BfMT2 cDNA resulting from an alternative splicing event. BfMT1 is a polyvalent metal binding peptide that coordinates any of the studied metal ions (Zn, Cd or Cu) rendering complexes stable enough to last in physiological environments, which is fully concordant with the constitutive expression of its gene, and therefore, with a metal homeostasis housekeeping role. On the contrary, BfMT2 exhibits a clear ability to coordinate Cd(II) ions, while it is absolutely unable to fold into stable Cu (I) complexes, even as mixed species. This identifies it as an essential detoxification agent, which is consequently only induced in emergency situations. The cephalochordate MTs are not directly related to vertebrate MTs, neither by gene structure, protein similarity nor metal-binding behavior of the encoded peptides. The closest relative is the echinoderm MT, which confirm proposed phylogenetic relationships between these two groups. The current findings support the existence in most organisms of two types of MTs as for their metal binding preferences, devoted to different biological functions: multivalent MTs for housekeeping roles, and specialized MTs that evolve either as Cd-thioneins or Cu-thioneins, according to the ecophysiological needs of each kind of organisms.

Ultrastructural characters of the spermatozoon of the digenean Hypocreadium caputvadum Kacem et al., 2011 (Lepocreadioidea: Lepocreadiidae), an intestinal parasite of Balistes capriscus in Tunisia

The ultrastructural organization of the spermatozoon of the digenean Hypocreadium caputvadum (Lepocreadioidea: Lepocreadiidae) is described. Live digeneans were collected from Balistes capriscus (Teleostei: Balistidae) from the Gulf of Gabès, Tunisia (Eastern Mediterranean Sea). The mature spermatozoon of H. caputvadum shows several ultrastructural characters such as two axonemes of different lengths exhibiting the classical 9 +"1" trepaxonematan pattern, a nucleus, two mitochondria, granules of glycogen, external ornamentation of the plasma membrane and two bundles of parallel cortical microtubules. Moreover, in the anterior extremity, the second axoneme is partly surrounded by a discontinuous and submembranous layer of electron-dense material. Our study provides new data on the spermatozoon of H. caputvadum in order to improve the understanding of phylogenetic relationships in the Digenea, particularly in the superfamily Lepocreadioidea. In this context, the electron-dense material surrounding one of the axonemes in the anterior spermatozoon extremity constitutes the unique distinguishing ultrastructural character of lepocreadioideans, and it is present in spermatozoa of lepocreadiids, aephnidiogenids and gyliauchenids.

Ultrastructural characters of the spermatozoon of the digenean Hypocreadium caputvadum Kacem et al., 2011 (Lepocreadioidea: Lepocreadiidae), an intestinal parasite of Balistes capriscus in Tunisia

The ultrastructural organization of the spermatozoon of the digenean Hypocreadium caputvadum (Lepocreadioidea: Lepocreadiidae) is described. Live digeneans were collected from Balistes capriscus (Teleostei: Balistidae) from the Gulf of Gabès, Tunisia (Eastern Mediterranean Sea). The mature spermatozoon of H. caputvadum shows several ultrastructural characters such as two axonemes of different lengths exhibiting the classical 9 +"1" trepaxonematan pattern, a nucleus, two mitochondria, granules of glycogen, external ornamentation of the plasma membrane and two bundles of parallel cortical microtubules. Moreover, in the anterior extremity, the second axoneme is partly surrounded by a discontinuous and submembranous layer of electron-dense material. Our study provides new data on the spermatozoon of H. caputvadum in order to improve the understanding of phylogenetic relationships in the Digenea, particularly in the superfamily Lepocreadioidea. In this context, the electron-dense material surrounding one of the axonemes in the anterior spermatozoon extremity constitutes the unique distinguishing ultrastructural character of lepocreadioideans, and it is present in spermatozoa of lepocreadiids, aephnidiogenids and gyliauchenids.

The odd couple: contrasting phylogeographic patterns in two sympatric sibling species of woodlouse-­‐hunter spiders in the Canary Islands

Comparative phylogeography seeks for commonalities in the spatial demographic history of sympatric organisms to characterize the mechanisms that shaped such patterns. The unveiling of incongruent phylogeographic patterns in co-occurring species, on the other hand, may hint to overlooked differences in their life histories or microhabitat preferences. The woodlouse-hunter spiders of the genus Dysdera have undergone a major diversi cation on the Canary Islands. The species pair Dysdera alegranzaensis and Dysdera nesiotes are endemic to the island of Lanzarote and nearby islets, where they co-occur at most of their known localities. The two species stand in sharp contrast to other sympatric endemic Dysdera in showing no evidence of somatic (non-genitalic) differentiation. Phylogenetic and population genetic analyses of mitochondrial cox1 sequences from an exhaustive sample of D. alegranzaensis and D. nesiotes specimens, and additional mitochondrial (16S, L1, nad1) and nuclear genes (28S, H3) were analysed to reveal their phylogeographic patterns and clarify their phylogenetic relationships. Relaxed molecular clock models using ve calibration points were further used to estimate divergence times between species and populations. Striking differences in phylogeography and population structure between the two species were observed. Dysdera nesiotes displayed a metapopulation-like structure, while D. alegranzaensis was characterized by a weaker geographical structure but greater genetic divergences among its main haplotype lineages, suggesting more complex population dynamics. Our study con rms that co-distributed sibling species may exhibit contrasting phylogeographic patterns in the absence of somatic differentiation. Further ecological studies, however, will be necessary to clarify whether the contrasting phylogeographies may hint at an overlooked niche partitioning between the two species. In addition, further comparisons with available phylogeographic data of other eastern Canarian Dysdera endemics con rm the key role of lava ows in structuring local populations in oceanic islands and identify localities that acted as refugia during volcanic eruptions

Climate niche conservatism does not explain restricted distribution patterns in Tynanthus (Bignonieae, Bignoniaceae)

Studies on niche evolution allow us to establish how species niches have changed over time as well as to identify how long-term evolutionary processes have led to present-day species distributions. Here, we investigate the patterns of climatic niche evolution in Tynanthus (Bignonieae, Bignoniaceae), a genus comprising narrowly distributed species. We test the hypothesis that niche conservatism has played an important role in the diversification history of this group of Neotropical lianas. For that, we perform univariate and multivariate comparisons between species' climatic niches and associated environmental data with information on species' phylogenetic relationships. We encountered considerable divergence in niches among species, indicating that niche conservatism in climatic variables has does not seem to havenot played a key role in the diversification of the genus. Our results are used as a basis to discuss patterns of ecological niche evolution in the group and to suggest novel approaches for future analyses.

Planarian (Platyhelminthes, Tricladida) Diversity and Molecular Markers: A New View of an Old Group

Planarians are a group of free-living platyhelminths (triclads) best-known largely due to long-standing regeneration and pattern formation research. However, the group"s diversity and evolutionary history has been mostly overlooked. A few taxonomists have focused on certain groups, resulting in the description of many species and the establishment of higher-level groups within the Tricladida. However, the scarcity of morphological features precludes inference of phylogenetic relationships among these taxa. The incorporation of molecular markers to study their diversity and phylogenetic relationships has facilitated disentangling many conundrums related to planarians and even allowed their use as phylogeographic model organisms. Here, we present some case examples ranging from delimiting species in an integrative style, and barcoding them, to analysing their evolutionary history on a lower scale to infer processes affecting biodiversity origin, or on a higher scale to understand the genus level or even higher relationships. In many cases, these studies have allowed proposing better classifications and resulted in taxonomical changes. We also explain shortcomings resulting in a lack of resolution or power to apply the most up-to-date data analyses. Next-generation sequencing methodologies may help improve this situation and accelerate their use as model organisms.

Typhlonesticus gocmeni sp. n., a new cave-dwelling blind spider species from the Aegean region of Turkey (Araneae, Nesticidae)

A new species of the troglobitic spider genus Typhlonesticus is described from specimens found in Keloğlan Cave (Denizli Province, Dodurgalar Town), Turkey. Typhlonesticus gocmeni sp. n. is described on the basis of both sexes; and its phylogenetic relationships with closely related European genera and species are discussed based on morphological and molecular data (the cox1, rrnL and H3 genes). Three new combinations are proposed: Typhlonesticus idriacus (Roewer, 1931), comb. n., Typhlonesticus morisii (Brignoli, 1975) comb. n. and Typhlonesticus obcaecatus (Simon, 1907), comb. n. all ex Nesticus.

Combining conservative and variable markers to inferthe evolutionary history of Prunus subgen. Amygdalus s.l.under domestication

The genus Prunus L. is large and economically important. However, phylogenetic relationships within Prunus at low taxonomic level, particularly in the subgenus Amygdalus L. s.l., remain poorly investigated. This paper attempts to document the evolutionary history of Amygdalus s.l. and establishes a temporal framework, by assembling molecular data from conservative and variable molecular markers. The nuclear s6pdh gene in combination with the plastid trnSG spacer are analyzed with bayesian and maximum likelihood methods. Since previous phylogenetic analysis with these markers lacked resolution, we additionally analyzed 13 nuclear SSR loci with the δµ2 distance, followed by an unweighted pair group method using arithmetic averages algorithm. Our phylogenetic analysis with both sequence and SSR loci confirms the split between sections Amygdalus and Persica, comprising almonds and peaches, respectively. This result is in agreement with biogeographic data showing that each of the two sections is naturally distributed on each side of the Central Asian Massif chain. Using coalescent based estimations, divergence times between the two sections strongly varied when considering sequence data only or combined with SSR. The sequence-only based estimate (5 million years ago) was congruent with the Central Asian Massif orogeny and subsequent climate change. Given the low level of differentiation within the two sections using both marker types, the utility of combining microsatellites and data sequences to address phylogenetic relationships at low taxonomic level within Amygdalus is discussed. The recent evolutionary histories of almond and peach are discussed in view of the domestication processes that arose in these two phenotypically-diverging gene pools: almonds and peaches were domesticated from the Amygdalus s.s. and Persica sections, respectively. Such economically important crops may serve as good model to study divergent domestication process in close genetic pool.

Phylogenomics Controlling for Base Compositional Bias Reveals a Single Origin of Eusociality in Corbiculate Bees.

As increasingly large molecular data sets are collected for phylogenomics, the conflicting phylogenetic signal among gene trees poses challenges to resolve some difficult nodes of the Tree of Life. Among these nodes, the phylogenetic position of the honey bees (Apini) within the corbiculate bee group remains controversial, despite its considerable importance for understanding the emergence and maintenance of eusociality. Here, we show that this controversy stems in part from pervasive phylogenetic conflicts among GC-rich gene trees. GC-rich genes typically have a high nucleotidic heterogeneity among species, which can induce topological conflicts among gene trees. When retaining only the most GC-homogeneous genes or using a nonhomogeneous model of sequence evolution, our analyses reveal a monophyletic group of the three lineages with a eusocial lifestyle (honey bees, bumble bees, and stingless bees). These phylogenetic relationships strongly suggest a single origin of eusociality in the corbiculate bees, with no reversal to solitary living in this group. To accurately reconstruct other important evolutionary steps across the Tree of Life, we suggest removing GC-rich and GC-heterogeneous genes from large phylogenomic data sets. Interpreted as a consequence of genome-wide variations in recombination rates, this GC effect can affect all taxa featuring GC-biased gene conversion, which is common in eukaryotes.

Genetic diversity of begomovirus infecting tomato and associated weeds in Southeastern Brazil

The genetic diversity of begomovirus isolates from tomato (Lycopersicon esculentum) fields in the Southeastern region of Brazil was analyzed by direct sequencing of PCR fragments amplified by using universal oligonucleotides for the begomovirus DNA-A, and subsequent computer-aided phylogenetic analysis. Samples of tomato plants and associated weeds showing typical symptoms of virus infection were collected at seven locations in the states of Minas Gerais, Espírito Santo and Rio de Janeiro. A total of 137 out of 369 samples were infected with a begomovirus based on PCR analysis. Phylogenetic analysis indicated a high degree of genetic diversity among begomoviruses infecting tomatoes in the sampled area. One species (Tomato chlorotic mottle virus, TCMV) occurs predominantly in Minas Gerais, whereas in Rio de Janeiro and Espírito Santo a distinct species, not yet fully characterized, predominates. Phylogenetic analysis further indicates the presence of an additional four possible new species. This high degree of genetic diversity suggests a recent transfer of indigenous begomovirus from wild hosts into tomatoes. The close phylogenetic relationship verified between begomovirus infecting tomato and associated weeds favors this hypothesis.