Molecular phylogenesis of Mediterranean Octocorals

In order to support the conservation of the Mediterranean octocorals improvements on information regarding their taxonomic units and phylogenetic relationships are strongly needed. In the present thesis work, phylogenetic analyses based on the mitochondrial mtMSH and 16S genes were performed including 15 Mediterranean octocorals species on the 56 recognized to date. Moreover, an extended datasets with Atlanto/Pacific congeners Octocorallia species was implemented to clarify their phylogenetic relationships and estimate the divergence times of the Mediterranean species. Results indicated that: 1) there are similarity and differences among molecular and morphological traits depending on the taxonomical level considered; 2) the molecular phylogeny of the Mediterranean octocorals retrace the previous relationships based on wide octocorals analyses; and 3) the divergence time among Mediterranean and Atlanto/Pacific species varies depending on analysed taxa. At higher taxonomic level, the Mediterranean trees supported the division of the Mediterranean Octocorallia into one major clade (Alcyoniina-Holaxonia) plus two unresolved branch including the single species available of Scleraxonia and Stolonifera respectively. This topology was better supported including the Atlanto/Pacific congeners species. The molecular evidence suggested that Alcyonium palmatum and Corallium rubrum species are the youngest with a divergence time estimated around 4 MYA. Particularly, C. rubrum results were in agreement with the hypothesis that recent orogenesis process of the Mediterranean Sea promoted the allopatric speciation of this specie. Increasing the sample design and implementing the emerging next-generation genomic-sequencing technologies, further studies would be able to improve the understanding of the Mediterranean octocorals phylogenetic relationships and evolution.

Rekonstruktion der eiszeitlichen Verbreitung und Artbildung vier alpiner Primeln durch Artenverbreitungsmodelle und Phylogeographie

Aim: Previous studies revealed that diversification events in the western clade of the alpine Primula sect. Auricula were concentrated in the Quaternary cold periods. This implies that allopatric speciation in isolated glacial refugia was the most common mode of speciation. In the first part of the present dissertation, this hypothesis is further investigated by locating refugial areas of two sister species, Primula marginata & P. latifolia during the last glacial maximum, 21,000 years ago. In the second part, the glacial and postglacial history of P. hirsuta and P. daonensis is investigated. Location: European Alps. Methods: Glacial refugia were located using species distribution models, which are projected to last glacial maximum climate. These refugia are validated with geographic distribution patterns of intra-specific genetic diversity, rarity and variation. Results 1) Speciation: Glacial refugia of the sister taxa Primula marginata and P. latifolia were largely separated, only a small overlapping zone at the southern margin of the former glacier in the Maritime Alps exists. This overlapping zone is too small to indicate sympatric speciation. The largely separated glacial distribution of both species rather confirms our hypothesis of allopatric speciation in isolated glacial refugia. Results 2) Glacial and postglacial history: Surprizingly, the modelled potential refugia of three out of four Primula species are situated within the former ice-shield, except for P. marginata. This indicates that peripheral and central nunataks played an important role for the glacial survival in P. latifolia, P. hirsuta and P. daonensis, while peripheral refugia outside the maximum extend of the glacier were crucial in P. marginata. In P. hirsuta and P. latifolia SDMs allowed to exclude several hypothetical refugial areas that overlap with today’s distribution as potential refugia for the species. In P. marginata, hypothetical refugial areas at the periphery of the former ice-shield that overlap with today’s distribution were confirmed by the models. The results from the SDMs are confirmed by population genetic patterns in three out of four species. P. daonensis represents an exception, where population genetic data contradict the SDMs. Main conclusions: Species distribution models provide species specific scenarios of glacial distribution and postglacial re-colonization, which can be validated using population genetic analyses. This combined approach is useful and helps to understand the complex processes that have lead to the genetic and floristic patterns of biodiversity that is found today in the Alps.

Landscape genetics of Phaedranassa herb. (Amaryllidaceae) in Ecuador

Speciation can be understood as a continuum occurring at different levels, from population to species. The recent molecular revolution in population genetics has opened a pathway towards understanding species evolution. At the same time, speciation patterns can be better explained by incorporating a geographic context, through the use of geographic information systems (GIS). Phaedranassa (Amaryllidaceae) is a genus restricted to one of the world’s most biodiverse hotspots, the Northern Andes. I studied seven Phaedranassa species from Ecuador. Six of these species are endemic to the country. The topographic complexity of the Andes, which creates local microhabitats ranging from moist slopes to dry valleys, might explain the patterns of Phaedranassa species differentiation. With a Bayesian individual assignment approach, I assessed the genetic structure of the genus throughout Ecuador using twelve microsatellite loci. I also used bioclimatic variables and species geographic coordinates under a Maximum Entropy algorithm to generate distribution models of the species. My results show that Phaedranassa species are genetically well-differentiated. Furthermore, with the exception of two species, all Phaedranassa showed non-overlapping distributions. Phaedranassa viridiflora and P. glauciflora were the only species in which the model predicted a broad species distribution, but genetic evidence indicates that these findings are likely an artifact of species delimitation issues. Both genetic differentiation and nonoverlapping geographic distribution suggest that allopatric divergence could be the general model of genetic differentiation. Evidence of sympatric speciation was found in two geographically and genetically distinct groups of P. viridiflora. Additionally, I report the first register of natural hybridization for the genus. The findings of this research show that the genetic differentiation of species in an intricate landscape as the Andes does not necessarily show a unique trend. Although allopatric speciation is the most common form of speciation, I found evidence of sympatric speciation and hybridization. These results show that the processes of speciation in the Andes have followed several pathways. The mixture of these processes contributes to the high biodiversity of the region.

Sampling information, length measurements, DNA alignments, models and phylogenetics trees for Oithona similis s.l.

Traditionally, many small-sized copepod species are considered to be widespread, bipolar or cosmopolitan. However, these large-scale distribution patterns need to be re-examined in view of increasing evidence of cryptic and pseudo-cryptic speciation in pelagic copepods. Here, we present a phylogeographic study of Oithona similis s.l. populations from the Arctic Ocean, the Southern Ocean and its northern boundaries, the North Atlantic and the Mediterrranean Sea. O. similis s.l. is considered as one of the most abundant species in temperate to polar oceans and acts as an important link in the trophic network between the microbial loop and higher trophic levels such as fish larvae. Two gene fragments were analysed: the mitochondrial cytochrome oxidase c subunit I (COI), and the nuclear ribosomal 28S genetic marker. Seven distinct, geographically delimitated, mitochondrial lineages could be identified, with divergences among the lineages ranging from 8 to 24 %, thus representing most likely cryptic or pseudocryptic species within O. similis s.l. Four lineages were identified within or close to the borders of the Southern Ocean, one lineage in the Arctic Ocean and two lineages in the temperate Northern hemisphere. Surprisingly the Arctic lineage was more closely related to lineages from the Southern hemisphere than to the other lineages from the Northern hemisphere, suggesting that geographic proximity is a rather poor predictor of how closely related the clades are on a genetic level. Molecular clock application revealed that the evolutionary history of O. similis s.l. is possibly closely associated with the reorganization of the ocean circulation in the mid Miocene and may be an example of allopatric speciation in the pelagic zone.

Comparison of calling songs in three allopatric populations of Endecous itatibensis (Orthoptera, Phalangopsinae)

Three isolated populations of the cricket Endecous itatibensis Rehn, 1918 had been analyzed to test the hypothesis of divergence in the calling song estridulation and to discuss its implications in the speciation process. The song registers were obtained from specimens in Corumbataí, Piraciacaba and Itatiba cities, state of São Paulo, Brazil. In the three analyzed populations, calling songs are composed by pairs of notes. Specimens from Corumbataí emit composed phrases with 3 to 18 pairs of notes, while those of Piracicaba and Itatiba have 1 to 3 pairs of notes. Inter- and intrapopulational variability in the frequencies of acoustic signals were detected. While recording the calling songs, other stridulations were registered, possibly from species of Orthoptera or Hemiptera. When the tapes were analyzed in the laboratory, it was verified that there was no overlap among the species' calling songs, in such a way that each species had their particular band of frequencies, avoiding interespecific acoustic competition.

Reproductive isolation between sympatric and allopatric Brazilian populations of Lutzomyia longipalpis s.l. (Diptera: Psychodidae)

Lutzomyia longipalpis s.l., the main vector of Leishmania chagasi in Latin America, is a species complex although the exact number of siblings is yet unknown. In Brazil, the siblings differ in male copulatory courtship songs and pheromones that most certainly act as pre-zygotic reproductive barriers. Here we analysed the reproductive isolation between three allopatric and two sympatric populations of Lu. longipalpis s.l. from Brazil. The results indicate a strong copulatory and pre-mating isolation between the three allopatric populations. In addition, the results also indicate a stronger pre-mating isolation between the two sympatric siblings than between the three allopatric ones, suggesting a role for reinforcement in the speciation of the Lu. longipalpis s.l. complex.

Can we continue to neglect genomic variation in introgression rates when inferring the history of speciation? A case study in a Mytilus hybrid zone.

The use of molecular data to reconstruct the history of divergence and gene flow between populations of closely related taxa represents a challenging problem. It has been proposed that the long-standing debate about the geography of speciation can be resolved by comparing the likelihoods of a model of isolation with migration and a model of secondary contact. However, data are commonly only fit to a model of isolation with migration and rarely tested against the secondary contact alternative. Furthermore, most demographic inference methods have neglected variation in introgression rates and assume that the gene flow parameter (Nm) is similar among loci. Here, we show that neglecting this source of variation can give misleading results. We analysed DNA sequences sampled from populations of the marine mussels, Mytilus edulis and M. galloprovincialis, across a well-studied mosaic hybrid zone in Europe and evaluated various scenarios of speciation, with or without variation in introgression rates, using an Approximate Bayesian Computation (ABC) approach. Models with heterogeneous gene flow across loci always outperformed models assuming equal migration rates irrespective of the history of gene flow being considered. By incorporating this heterogeneity, the best-supported scenario was a long period of allopatric isolation during the first three-quarters of the time since divergence followed by secondary contact and introgression during the last quarter. By contrast, constraining migration to be homogeneous failed to discriminate among any of the different models of gene flow tested. Our simulations thus provide statistical support for the secondary contact scenario in the European Mytilus hybrid zone that the standard coalescent approach failed to confirm. Our results demonstrate that genomic variation in introgression rates can have profound impacts on the biological conclusions drawn from inference methods and needs to be incorporated in future studies.

Speciation in the Sorex araneus (Mammalia, Insectivora)

The seven cytologically known species of the Sorex araneus complex each show a different karyotype. The three European species are genetically and morphologically very close, and mostly allopatric, thus having similar ecological requirements. It is suggested that chromosomic changes occasion a stasipatric speciation, the formation of chromosomic races in S. araneus being an illustration of this process. Neither small demes nor high inbreeding are likely in Sorex; thus the Robertsonian changes should spread out by meiotic drive alone, and the hybrids for a change should be almost fully fertile. Chromosomic speciation occurs when two incompatible metacentrics meet together. This model explains the separation between S. araneus and S. coronatus, but geographical isolation must have occured in the case of S. granarius, which keeps a primitive karyotype.

Postmating-prezygotic isolation between two allopatric populations of Drosophila montana: fertilisation success differs under sperm competition

Postmating but prezygotic (PMPZ) interactions are increasingly recognized as a potentially important early-stage barrier in the evolution of reproductive isolation. A recent study described a potential example between populations of the same species: single matings between Drosophila montana populations resulted in differential fertilisation success because of the inability of sperm from one population (Vancouver) to penetrate the eggs of the other population (Colorado). As the natural mating system of D. montana is polyandrous (females remate rapidly), we set up double matings of all possible crosses between the same populations to test whether competitive effects between ejaculates influence this PMPZ isolation. We measured premating isolation in no-choice tests, female fecundity, fertility and egg-to-adult viability after single and double matings as well as second-male paternity success (P-2). Surprisingly, we found no PMPZ reproductive isolation between the two populations under a competitive setting, indicating no difficulty of sperm from Vancouver males to fertilize Colorado eggs after double matings. While there were subtle differences in how P-2 changed over time, suggesting that Vancouver males' sperm are somewhat less competitive in a first-male role within Colorado females, these effects did not translate into differences in overall P-2. Fertilisation success can thus differ dramatically between competitive and noncompetitive conditions, perhaps because the males that mate second produce higher quality ejaculates in response to sperm competition. We suggest that unlike in more divergent species comparisons, where sperm competition typically increases reproductive isolation, ejaculate tailoring can reduce the potential for PMPZ isolation when recently diverged populations interbreed.

Genome-wide evidence for speciation with gene flow in Heliconius butterflies

Most speciation events probably occur gradually, without complete and immediate reproductive isolation, but the full extent of gene flow between diverging species has rarely been characterized on a genome-wide scale. Documenting the extent and timing of admixture between diverging species can clarify the role of geographic isolation in speciation. Here we use new methodology to quantify admixture at different stages of divergence in Heliconius butterflies, based on whole-genome sequences of 31 individuals. Comparisons between sympatric and allopatric populations of H. melpomene, H. cydno, and H. timareta revealed a genome-wide trend of increased shared variation in sympatry, indicative of pervasive interspecific gene flow. Up to 40% of 100-kb genomic windows clustered by geography rather than by species, demonstrating that a very substantial fraction of the genome has been shared between sympatric species. Analyses of genetic variation shared over different time intervals suggested that admixture between these species has continued since early in speciation. Alleles shared between species during recent time intervals displayed higher levels of linkage disequilibrium than those shared over longer time intervals, suggesting that this admixture took place at multiple points during divergence and is probably ongoing. The signal of admixture was significantly reduced around loci controlling divergent wing patterns, as well as throughout the Z chromosome, consistent with strong selection for Müllerian mimicry and with known Z-linked hybrid incompatibility. Overall these results show that species divergence can occur in the face of persistent and genome-wide admixture over long periods of time.

Cyclamen: time, sea and speciation biogeography using a temporally calibrated phylogeny

Aim The Mediterranean region is a species-rich area with a complex geographical history. Geographical barriers have been removed and restored due to sea level changes and local climatic change. Such barriers have been proposed as a plausible mechanism driving the high levels of speciation and endemism in the Mediterranean basin. This raises the fundamental question: is allopatric isolation the mechanism by which speciation occurs? This study explores the potential driving influence of palaeo-geographical events on the speciation of Cyclamen (Myrsinaceae), a group with most species endemic to the Mediterranean region. Cyclamen species have been shown experimentally to have few genetic barriers to hybridization. Location The Mediterranean region, including northern Africa, extending eastwards to the Black Sea coast. Methods A generic level molecular phylogeny of Myrsinaceae and Primulaceae is constructed, using Bayesian approximation, to produce a secondary age estimate for the stem lineage of Cyclamen. This estimate is used to calibrate temporally an infrageneric phylogeny of Cyclamen, built with nrDNA ITS, cpDNA trnL-F and cpDNA rps16 sequences. A biogeographical analysis of Cyclamen is performed using dispersal-vicariance analysis. Results The emergence of the Cyclamen stem lineage is estimated at 30.1-29.2 Ma, and the crown divergence at 12.9-12.2 Ma. The average age of Cyclamen species is 3.7 Myr. Every pair of sister species have mutually exclusive, allopatric distributions relative to each other. This pattern appears typical of divergence events throughout the evolutionary history of the genus. Main conclusions Geographical barriers, such as the varying levels of the Mediterranean Sea, are the most plausible explanation for speciation events throughout the phylogenetic history of Cyclamen. The genus demonstrates distributional patterns congruent with the temporally reticulate palaeogeography of the Mediterranean region.

A comparative structural cytogenetic study in three allopatric populations of Astyanax scabripinnis (Teleostei: Characidae)

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

How does the taxonomic status of allopatric populations influence species richness within African cichlid fish assemblages?

Aim  Current estimates of species richness within rapidly evolving species flocks are often highly dependent on the species status of allopatric populations that differ in phenotypic traits. These traits may be unreliable indicators of biological species status and systematists may have inconsistently assigned species among lineages or locations on the basis of these traits, thus hampering comparative studies of regional species richness and speciation rates. Our aim was to develop a method of generating standardized estimates of regional species richness suitable for comparative analysis, and to use these estimates to examine the extent and consistency of species assignment of allopatric populations within rapidly evolving cichlid fish flocks present in three east African lakes. Location  Lakes Malawi, Victoria and Tanganyika. Methods  Using published taxon co-occurrence data, a novel approach was employed to calculate standardized ‘minimum’ estimates of regional species richness for hard substrate associated complexes of cichlids within each of the lakes. Minimum estimates were based on an explicit assumption that if taxa present on equivalent habitats have disjunct distributions, then they are allopatric forms of the same species. These estimates were compared with current observed ‘high-end’ regional species richness estimates for those complexes to determine the consistency of species assignment of allopatric populations between lineages within a lake. A ‘sympatry’ index was developed to enable comparisons of levels of species assignment of allopatric populations between-lakes to be made. Results  Within each lake, the minimum and high-end estimates for species richness were significantly correlated across complexes, indicating that the complexes that contain more recognized species contain the most genuine biological species. However, comparisons of complexes among lakes revealed considerable differences. For equivalent geographical areas, substantially higher proportions of recognized species were totally allopatric within the studied Lake Malawi and Lake Victoria complexes, than those of Lake Tanganyika. Main Conclusions  Among African lakes, levels of assignment to species status of allopatric populations were found to be distinctly different. It is unclear whether the discrepancies are a consequence of differences between the lake faunas in degrees of phenotypic divergence among allopatric populations, or are simply the result of inconsistent taxonomic practices. In either case, these results have considerable wider relevance for they emphasize that quantitative measures of regional and beta diversity are critically dependent on the species status of allopatric populations, an issue usually neglected in comparative studies of species richness. The technique introduced here can be used to standardize measures of regional diversity of lineages for comparative analyses, potentially enabling more accurate identification of processes influencing rates of speciation.

Can sympatric speciation by disruptive sexual selection explain rapid evolution of cichlid diversity in Lake Victoria?

Rapid speciation can occur on ecological time scales and interfere with ecological processes, resulting in species distribution patterns that are difficult to reconcile with ecological theory. The haplochromine cichlids in East African lakes are an extreme example of rapid speciation. We analyse the causes of their high speciation rates. Various studies have identified disruptive sexual selection acting on colour polymorphisms that might cause sympatric speciation. Using data on geographical distribution, colouration and relatedness from 41 species endemic to Lake Victoria, we test predictions from this hypothesis. Plotting numbers of pairs of closely related species against the amount of distributional overlap between the species reveals a bimodal distribution with modes on allopatric and sympatric. The proportion of sister species pairs that are heteromorphic for the traits under disruptive selection is higher in sympatry than in allopatry. These data support the hypothesis that disruptive sexual selection on colour polymorphisms has caused sympatric speciation and help to explain the rapid radiation of haplochromine species flocks.

Genomics of Rapid Incipient Speciation in Sympatric Threespine Stickleback

Ecological speciation is the process by which reproductively isolated populations emerge as a consequence of divergent natural or ecologically-mediated sexual selection. Most genomic studies of ecological speciation have investigated allopatric populations, making it difficult to infer reproductive isolation. The few studies on sympatric ecotypes have focused on advanced stages of the speciation process after thousands of generations of divergence. As a consequence, we still do not know what genomic signatures of the early onset of ecological speciation look like. Here, we examined genomic differentiation among migratory lake and resident stream ecotypes of threespine stickleback reproducing in sympatry in one stream, and in parapatry in another stream. Importantly, these ecotypes started diverging less than 150 years ago. We obtained 34,756 SNPs with restriction-site associated DNA sequencing and identified genomic islands of differentiation using a Hidden Markov Model approach. Consistent with incipient ecological speciation, we found significant genomic differentiation between ecotypes both in sympatry and parapatry. Of 19 islands of differentiation resisting gene flow in sympatry, all were also differentiated in parapatry and were thus likely driven by divergent selection among habitats. These islands clustered in quantitative trait loci controlling divergent traits among the ecotypes, many of them concentrated in one region with low to intermediate recombination. Our findings suggest that adaptive genomic differentiation at many genetic loci can arise and persist in sympatry at the very early stage of ecotype divergence, and that the genomic architecture of adaptation may facilitate this.