A timescale for diatom evolution based on four molecular markers: reassessment of ghost lineages and major steps defining diatom evolution.

ABSTRACT. – Phylogenies and molecular clocks of the diatoms have largely been inferred from SSU rDNA sequences. A new phylogeny of diatoms was estimated using four gene markers SSU and LSU rDNA rbcL and psbA (total 4352 bp) with 42 diatom species. The four gene trees analysed with a maximum likelihood (ML) and Baysian (BI) analysis recovered a monophyletic origin of the new diatom classes with high bootstrap support, which has been controversial with single gene markers using single outgroups and alignments that do not take secondary structure of the SSU gene into account. The divergence time of the classes were calculated from a ML tree in the MultliDiv Time program using a Bayesian estimation allowing for simultaneous constraints from the fossil record and varying rates of molecular evolution of different branches in the phylogenetic tree. These divergence times are generally in agreement with those proposed by other clocks using single genes with the exception that the pennates appear much earlier and suggest a longer Cretaceous fossil record that has yet to be sampled. Ghost lineages (i.e. the discrepancy between first appearance (FA) and molecular clock age of origin from an extant taxon) were revealed in the pennate lineage, whereas those ghost lineages in the centric lineages previously reported by others are reviewed and referred to earlier literature.

Host jumps shaped the diversity of extant rust fungi (Pucciniales)

* The aim of this study was to determine the evolutionary time line for rust fungi and date key speciation events using a molecular clock. Evidence is provided that supports a contemporary view for a recent origin of rust fungi, with a common ancestor on a flowering plant. * Divergence times for > 20 genera of rust fungi were studied with Bayesian evolutionary analyses. A relaxed molecular clock was applied to ribosomal and mitochondrial genes, calibrated against estimated divergence times for the hosts of rust fungi, such as Acacia (Fabaceae), angiosperms and the cupressophytes. * Results showed that rust fungi shared a most recent common ancestor with a mean age between 113 and 115 million yr. This dates rust fungi to the Cretaceous period, which is much younger than previous estimations. Host jumps, whether taxonomically large or between host genera in the same family, most probably shaped the diversity of rust genera. Likewise, species diversified by host shifts (through coevolution) or via subsequent host jumps. This is in contrast to strict coevolution with their hosts. * Puccinia psidii was recovered in Sphaerophragmiaceae, a family distinct from Raveneliaceae, which were regarded as confamilial in previous studies.

Comparative phylogeography of alpine invertebrates indicates deep lineage diversification and historical refugia in the Australian Alps

Aim: Comparative phylogeographic analyses of alpine biota from the Northern Hemisphere have linked patterns of genetic diversification to glacial expansion and contraction events in the Pliocene and Pleistocene. Furthermore, the extent of diversification across species groups appears to be associated with vagility. In this study we test whether these patterns apply to a geologically stable system from eastern Australia with comparatively shallow elevational gradients and minimal influence from historical glacial activity. Location: The Australian Alps, Victoria, eastern Australia. Methods: We considered phylogeographic patterns across five alpine invertebrate species based on mitochondrial and nuclear DNA sequence data. Bayesian inference methods were used to estimate species phylogenies and divergence times among lineages. GIS tools were used to map interpopulation genetic divergence and intrapopulation genetic diversity estimates and to visualize spatial patterns across species, providing insights into patterns of endemism and demographic history. Results: Phylogeographic patterns and the timing of lineage diversification were consistent across taxonomic groups. Mountain summits harbour highly differentiated haplogroups, including summits connected by high-elevational plateaus, pointing to diversifications being maintained since the early to mid-Pleistocene. These findings are consistent with previous studies of alpine mammals and reptiles, demonstrating a high degree of endemism in this region, regardless of species vagility. Main conclusions: The fine spatial scales at which deep genetic differentiation among alpine communities was observed in this study are unprecedented. This suggests that glacial periods have had less of an impact on species distributions and genetic diversity than they have in alpine systems in the Northern Hemisphere. Historical gene flow among sky-island populations has been limited despite connecting snowlines during glacial periods, suggesting that factors other than snow cover have influenced patterns of gene flow in this region. These findings emphasize the unique phylogeographic history affecting Victorian alpine biodiversity, and the importance of conserving biodiversity from multiple mountain summits in this region of high endemism.

Concordance in phylogeography and ecological niche modelling identify dispersal corridors for reptiles in arid Australia

Aim: Using the rock-specialist agamid Ctenophorus caudicinctus as a model, we test hypothesized biogeographical dispersal corridors for lizards in the Australian arid zone (across the western sand deserts), and assess how these dispersal routes have shaped phylogeographical structuring. Location: Arid and semi-arid Australia. Methods: We sequenced a c. 1400 bp fragment of mtDNA (ND2) for 134 individuals of C. caudicinctus as well as a subset of each of the mtDNA clades for five nuclear loci (BDNF, BACH1, GAPD, NTF3, and PRLR). We used phylogenetic methods to assess biogeographical patterns within C. caudicinctus, including relaxed molecular clock analyses to estimate divergence times. Ecological niche modelling (Maxent) was employed to estimate the current distribution of suitable climatic envelopes for each lineage. Results: Phylogenetic analyses identified two deeply divergent mtDNA clades within C. caudicinctus - an eastern and western clade - separated by the Western Australian sand deserts. However, divergences pre-date the Pleistocene sand deserts. Phylogenetic analyses of the nuclear DNA data sets generally support major mtDNA clades, suggesting past connections between the western C. c. caudicinctus populations in far eastern Pilbara (EP) and the lineages to the east of the sand deserts. Ecological niche modelling supports the continued suitability of climatic conditions between the Central Ranges and the far EP for C. c. graafi. Main conclusions: Estimates of lineage ages provide evidence of divergence between eastern and western clades during the Miocene with subsequent secondary contact during the Pliocene. Our results suggest that this secondary contact occurred via dispersal between the Central Ranges and the far EP, rather than the more southerly Giles Corridor. These events precede the origins of the western sand deserts and divergence patterns instead appear associated with Miocene and Pliocene climate change.

Propagation of low power low divergence Gaussian fields in unbiased self-defocusing photorefractive media and their interactions

Many optical networks are limited in speed and processing capability due to the necessity for the optical signal to be converted to an electrical signal and back again. In addition, electronically manipulated interconnects in an otherwise optical network lead to overly complicated systems. Optical spatial solitons are optical beams that propagate without spatial divergence. They are capable of phase dependent interactions, and have therefore been extensively researched as suitable all optical interconnects for over 20 years. However, they require additional external components, initially high voltage power sources were required, several years later, high power background illumination had replaced the high voltage. However, these additional components have always remained as the greatest hurdle in realising the applications of the interactions of spatial optical solitons as all optical interconnects. Recently however, self-focusing was observed in an otherwise self-defocusing photorefractive crystal. This observation raises the possibility of the formation of soliton-like fields in unbiased self-defocusing media, without the need for an applied electrical field or background illumination. This thesis will present an examination of the possibility of the formation of soliton-like low divergence fields in unbiased self-defocusing photorefractive media. The optimal incident beam and photorefractive media parameters for the formation of these fields will be presented, together with an analytical and numerical study of the effect of these parameters. In addition, preliminary examination of the interactions of two of these fields will be presented. In order to complete an analytical examination of the field propagating through the photorefractive medium, the spatial profile of the beam after propagation through the medium was determined. For a low power solution, it was found that an incident Gaussian field maintains its Gaussian profile as it propagates. This allowed the beam at all times to be described by an individual complex beam parameter, while also allowing simple analytical solutions to the appropriate wave equation. An analytical model was developed to describe the effect of the photorefractive medium on the Gaussian beam. Using this model, expressions for the required intensity dependent change in both the real and imaginary components of the refractive index were found. Numerical investigation showed that under certain conditions, a low powered Gaussian field could propagate in self-defocusing photorefractive media with divergence of approximately 0.1 % per metre. An investigation into the parameters of a Ce:BaTiO3 crystal showed that the intensity dependent absorption is wavelength dependent, and can in fact transition to intensity dependent transparency. Thus, with careful wavelength selection, the required intensity dependent change in both the real and imaginary components of the refractive index for the formation of a low divergence Gaussian field are physically realisable. A theoretical model incorporating the dependence of the change in real and imaginary components of the refractive index on propagation distance was developed. Analytical and numerical results from this model are congruent with the results from the previous model, showing low divergence fields with divergence less than 0.003 % over the propagation length of the photorefractive medium. In addition, this approach also confirmed the previously mentioned self-focusing effect of the self-defocusing media, and provided an analogy to a negative index GRIN lens with an intensity dependent focal length. Experimental results supported the findings of the numerical analysis. Two low divergence fields were found to possess the ability to interact in a Ce:BaTiO3 crystal in a soliton-like fashion. The strength of these interactions was found to be dependent on the degree of divergence of the individual beams. This research found that low-divergence fields are possible in unbiased self-defocusing photorefractive media, and that soliton-like interactions between two of these fields are possible. However, in order for these types of fields to be used in future all optical interconnects, the manipulation of these interactions, together with the ability for these fields to guide a second beam at a different wavelength, must be investigated.

Substantial genetic divergence between morphologically indistinguishable populations of Fasciola suggests the possibility of cryptic speciation.

The liver flukes, Fasciola hepatica and Fasciola gigantica, are considered to be sister species and between them present a major threat worldwide to livestock production. In this study sequence data have been employed from informative regions of the nuclear and mitochondrial genomes of over 200 morphologically F. hepatica-like or F. gigantica-like flukes from Europe, sub-Saharan Africa and South Asia to assess genetic diversity. Evidence is presented for the existence of four well-separated clades: African gigantica-like flukes, Indian gigantica-like flukes, European hepatica-like flukes and African high-altitude hepatica-like flukes. Application of the Biological Species Concept to trematodes is problematic; however, the degree of separation between these groups was sufficient for them to be considered as distinct species using the four times rule for speciation.

Identifying adaptive genetic divergence among populations from genome scans

The identification of signatures of natural selection in genomic surveys has become an area of intense research, stimulated by the increasing ease with which genetic markers can be typed. Loci identified as subject to selection may be functionally important, and hence (weak) candidates for involvement in disease causation. They can also be useful in determining the adaptive differentiation of populations, and exploring hypotheses about speciation. Adaptive differentiation has traditionally been identified from differences in allele frequencies among different populations, summarised by an estimate of F-ST. Low outliers relative to an appropriate neutral population-genetics model indicate loci subject to balancing selection, whereas high outliers suggest adaptive (directional) selection. However, the problem of identifying statistically significant departures from neutrality is complicated by confounding effects on the distribution of F-ST estimates, and current methods have not yet been tested in large-scale simulation experiments. Here, we simulate data from a structured population at many unlinked, diallelic loci that are predominantly neutral but with some loci subject to adaptive or balancing selection. We develop a hierarchical-Bayesian method, implemented via Markov chain Monte Carlo (MCMC), and assess its performance in distinguishing the loci simulated under selection from the neutral loci. We also compare this performance with that of a frequentist method, based on moment-based estimates of F-ST. We find that both methods can identify loci subject to adaptive selection when the selection coefficient is at least five times the migration rate. Neither method could reliably distinguish loci under balancing selection in our simulations, even when the selection coefficient is twenty times the migration rate.

Coalescent analysis of mtDNA indicates Pleistocene divergence among three species of howler monkey (Alouatta spp.) and population subdivision within the Atlantic Coastal Forest species, A. guariba

We have used coalescent analysis of mtDNA cytochrome b (cyt b) sequences to estimate times of divergence of three species of Alouatta-A. caraya, A. belzebul, and A. guariba-which are in close geographic proximity. A. caraya is inferred to have diverged from the A. guariba/A. belzebul clade approximately 3.83 million years ago (MYA), with the later pair diverging approximately 1.55 MYA. These dates are much more recent than previous dates based on molecular-clock methods. In addition, analyses of new sequences from the Atlantic Coastal Forest species A. guariba indicate the presence of two distinct haplogroups corresponding to northern and southern populations with both haplogroups occurring in sympatry within Sao Paulo state. The time of divergence of these two haplogroups is estimated to be 1.2 MYA and so follows quite closely after the divergence of A. guariba and A. belzebul. These more recent dates point to the importance of Pleistocene environmental events as important factors in the diversification of A. belzebul and A. guariba. We discuss the diversification of the three Alouatta species in the context of recent models of climatic change and with regard to recent molecular phylogeographic analyses of other animal groups distributed in Brazil.

Divergence patterns in Fucus seaweeds (Phaeophyceae) in the northern Baltic Sea and in the Tjongspollen area

To be able to interpret patterns of biodiversity it is important to understand the processes by which new species evolve and how closely related species remain reproductively isolated and ecologically differentiated. Divergence and differentiation can vary during speciation and it can be seen in different stages. Groups of closely related taxa constitute important case studies to understand species and new biodiversity formation. However, it is important to assess the divergence among them at different organismal levels and from an integrative perspective. For this purpose, this study used the brown seaweeds genus Fucus as a model to study speciation, as they constitute a good opportunity to study divergence at different stages. We investigated the divergence patterns in Fucus species from two marginal areas (northern Baltic Sea and the Tjongspollen area), based on phenetic, phylogenetic and biological taxonomical criteria that are respectively characterised by algal morphology, allele frequencies of five microsatellite loci and levels of secondary polyphenolic compounds called phlorotannins. The results from this study showed divergence at morphological and genetic levels to certain extent but complete lack of divergence at biochemical level (i.e. constitutive phlorotannin production) in the Baltic Sea or Norway. Morphological divergence was clearly evident in Tjongspollen (Norway) among putative taxa as they were identified in the field and this divergence corresponds with their neutral genetic divergence. In the Baltic, there are some distinguishable patterns in the morphology of the swedish and finnish individuals according to locality to certain extent but not among putative taxa within localities. Likewise, these morphological patterns have genetic correspondence among localities but not within each locality. At the biochemical level, measured by the phlorotannin contents there were neither evidence of divergence in Norway or the Baltic Sea nor any discernable aggregation pattern among or within localities. Our study have contributed with further understanding of the Baltic Sea Fucus system and its intriguingly rapid and recent divergence as well as of the Tjongspollen area systems where formally undescribed individuals have been observed for the first time; in fact they appear largely differentiated and they may well warrant a new species status. In current times, climate change threatens, peripheral ecosystems, biodiversity, and increased knowledge of processes generating and maintaining biodiversity in those ecosystems seem particularly important and needed.

Intergenomic epistasis causes asynchronous hatch times in whitefish hybrids, but only when parental ecotypes differ

Support for the theory of ecological speciation requires evidence for ecological divergence between species which directly or indirectly causes reproductive isolation. This study investigates effects of ecological vs. genetic disparity of parental species on the presence of endogenous selection (deformation and mortality rates) and potential sources of exogenous selection (growth rates and hatch timing) on hybrids. Hybrid embryonic development is analysed in a common-garden full-sib cross of three species belonging to two different ecotypes within the Coregonus lavaretus species flock in the central Alpine region of Europe. Although hatch timing was similar across the three species, embryonic growth rates and egg sizes differed between ecotypes. This led to a mismatch between embryonic growth rate and egg size in hybrid crosses that reveals epistasis between the maternal and embryonic genomes and transgressive hatch times that were asynchronous with control crosses. A strong constraint of egg size to embryo size at late development was also evident. We argue that this demonstrates potential for coadaptation of a maternal trait (egg size) with offspring growth rate to be an important source of selection against hybridization between ecotypes with different egg sizes. Implications for the measurement and quantification of early life-history traits affected by this additive relationship, such as hatch day and larval size, are also discussed.

Contemporary ecotypic divergence during a recent range expansion was facilitated by adaptive introgression

Although rapid phenotypic evolution during range expansion associated with colonization of contrasting habitats has been documented in several taxa, the evolutionary mechanisms that underlie such phenotypic divergence have less often been investigated. A strong candidate for rapid ecotype formation within an invaded range is the three-spine stickleback in the Lake Geneva region of central Europe. Since its introduction only about 140 years ago, it has undergone a significant expansion of its range and its niche, now forming phenotypically differentiated parapatric ecotypes that occupy either the pelagic zone of the large lake or small inlet streams, respectively. By comparing museum collections from different times with contemporary population samples, we here reconstruct the evolution of parapatric phenotypic divergence through time. Using genetic data from modern samples, we infer the underlying invasion history. We find that parapatric habitat-dependent phenotypic divergence between the lake and stream was already present in the first half of the twentieth century, but the magnitude of differentiation increased through time, particularly in antipredator defence traits. This suggests that divergent selection between the habitats occurred and was stable through much of the time since colonization. Recently, increased phenotypic differentiation in antipredator defence traits likely results from habitat-dependent selection on alleles that arrived through introgression from a distantly related lineage from outside the Lake Geneva region. This illustrates how hybridization can quickly promote phenotypic divergence in a system where adaptation from standing genetic variation was constrained.