Extensive remodeling of sugar metabolism through gene loss and horizontal gene transfer in a eukaryotic lineage.

BACKGROUND: In yeasts belonging to the subphylum Saccharomycotina, genes encoding components of the main metabolic pathways, like alcoholic fermentation, are usually conserved. However, in fructophilic species belonging to the floral Wickerhamiella and Starmerella genera (W/S clade), alcoholic fermentation was uniquely shaped by events of gene loss and horizontal gene transfer (HGT). RESULTS: Because HGT and gene losses were first identified when only eight W/S-clade genomes were available, we collected publicly available genome data and sequenced the genomes of 36 additional species. A total of 63 genomes, representing most of the species described in the clade, were included in the analyses. Firstly, we inferred the phylogenomic tree of the clade and inspected the genomes for the presence of HGT-derived genes involved in fructophily and alcoholic fermentation. We predicted nine independent HGT events and several instances of secondary loss pertaining to both pathways. To investigate the possible links between gene loss and acquisition events and evolution of sugar metabolism, we conducted phenotypic characterization of 42 W/S-clade species including estimates of sugar consumption rates and fermentation byproduct formation. In some instances, the reconciliation of genotypes and phenotypes yielded unexpected results, such as the discovery of fructophily in the absence of the cornerstone gene (FFZ1) and robust alcoholic fermentation in the absence of the respective canonical pathway. CONCLUSIONS: These observations suggest that reinstatement of alcoholic fermentation in the W/S clade triggered a surge of innovation that goes beyond the utilization of xenologous enzymes, with fructose metabolism playing a key role.

Asynchronous transcription and translation of neurotransmitter-related genes characterize the initial stages of neuronal maturation in Drosophila.

Neuron specification and maturation are essential for proper central nervous system development. However, the precise mechanisms that govern neuronal maturation, essential to shape and maintain neuronal circuitry, remain poorly understood. Here, we analyse early-born secondary neurons in the Drosophila larval brain, revealing that the early maturation of secondary neurons goes through 3 consecutive phases: (1) Immediately after birth, neurons express pan-neuronal markers but do not transcribe terminal differentiation genes; (2) Transcription of terminal differentiation genes, such as neurotransmitter-related genes VGlut, ChAT, or Gad1, starts shortly after neuron birth, but these transcripts are, however, not translated; (3) Translation of neurotransmitter-related genes only begins several hours later in mid-pupa stages in a coordinated manner with animal developmental stage, albeit in an ecdysone-independent manner. These results support a model where temporal regulation of transcription and translation of neurotransmitter-related genes is an important mechanism to coordinate neuron maturation with brain development.

A glimpse at an early stage of microbe domestication revealed in the variable genome of Torulaspora delbrueckii, an emergent industrial yeast.

Microbe domestication has a major applied relevance but is still poorly understood from an evolutionary perspective. The yeast Torulaspora delbrueckii is gaining importance for biotechnology but little is known about its population structure, variation in gene content or possible domestication routes. Here, we show that T. delbrueckii is composed of five major clades. Among the three European clades, a lineage associated with the wild arboreal niche is sister to the two other lineages that are linked to anthropic environments, one to wine fermentations and the other to diverse sources including dairy products and bread dough (Mix-Anthropic clade). Using 64 genomes we assembled the pangenome and the variable genome of T. delbrueckii. A comparison with Saccharomyces cerevisiae indicated that the weight of the variable genome in the pangenome of T. delbrueckii is considerably smaller. An association of gene content and ecology supported the hypothesis that the Mix-Anthropic clade has the most specialized genome and indicated that some of the exclusive genes were implicated in galactose and maltose utilization. More detailed analyses traced the acquisition of a cluster of GAL genes in strains associated with dairy products and the expansion and functional diversification of MAL genes in strains isolated from bread dough. In contrast to S. cerevisiae, domestication in T. delbrueckii is not primarily driven by alcoholic fermentation but rather by adaptation to dairy and bread-production niches. This study expands our views on the processes of microbe domestication and on the trajectories leading to adaptation to anthropic niches.

Phylogenetic analysis and accessory genome diversity reveal insight into the evolutionary history of Streptococcus dysgalactiae.

Streptococcus dysgalactiae (SD) is capable of infecting both humans and animals and causing a wide range of invasive and non-invasive infections. With two subspecies, the taxonomic status of subspecies of SD remains controversial. Subspecies equisimilis (SDSE) is an important human pathogen, while subspecies dysgalactiae (SDSD) has been considered a strictly animal pathogen; however, occasional human infections by this subspecies have been reported in the last few years. Moreover, the differences between the adaptation of SDSD within humans and other animals are still unknown. In this work, we provide a phylogenomic analysis based on the single-copy core genome of 106 isolates from both the subspecies and different infected hosts (animal and human hosts). The accessory genome of this species was also analyzed for screening of genes that could be specifically involved with adaptation to different hosts. Additionally, we searched putatively adaptive traits among prophage regions to infer the importance of transduction in the adaptation of SD to different hosts. Core genome phylogenetic relationships segregate all human SDSE in a single cluster separated from animal SD isolates. The subgroup of bovine SDSD evolved from this later clade and harbors a specialized accessory genome characterized by the presence of specific virulence determinants (e.g., cspZ) and carbohydrate metabolic functions (e.g., fructose operon). Together, our results indicate a host-specific SD and the existence of an SDSD group that causes human-animal cluster infections may be due to opportunistic infections, and that the exact incidence of SDSD human infections may be underestimated due to failures in identification based on the hemolytic patterns. However, more detailed research into the isolation of human SD is needed to assess whether it is a carrier phenomenon or whether the species can be permanently integrated into the human microbiome, making it ready to cause opportunistic infections.

A regeneration-triggered metabolic adaptation is necessary for cell identity transitions and cell cycle re-entry to support blastema formation and bone regeneration.

Regeneration depends on the ability of mature cells at the injury site to respond to injury, generating tissue-specific progenitors that incorporate the blastema and proliferate to reconstitute the original organ architecture. The metabolic microenvironment has been tightly connected to cell function and identity during development and tumorigenesis. Yet, the link between metabolism and cell identity at the mechanistic level in a regenerative context remains unclear. The adult zebrafish caudal fin, and bone cells specifically, have been crucial for the understanding of mature cell contribution to tissue regeneration. Here, we use this model to explore the relevance of glucose metabolism for the cell fate transitions preceding new osteoblast formation and blastema assembly. We show that injury triggers a modulation in the metabolic profile at early stages of regeneration to enhance glycolysis at the expense of mitochondrial oxidation. This metabolic adaptation mediates transcriptional changes that make mature osteoblast amenable to be reprogramed into pre-osteoblasts and induces cell cycle re-entry and progression. Manipulation of the metabolic profile led to severe reduction of the pre-osteoblast pool, diminishing their capacity to generate new osteoblasts, and to a complete abrogation of blastema formation. Overall, our data indicate that metabolic alterations have a powerful instructive role in regulating genetic programs that dictate fate decisions and stimulate proliferation, thereby providing a deeper understanding on the mechanisms regulating blastema formation and bone regeneration.

Phaffia brasiliana sp. nov., a yeast species isolated from soil in a Cerrado-Atlantic Rain Forest ecotone site in Brazil.

During studies of yeasts associated with soil in a Cerrado-Atlantic Rain Forest ecotone site in Brazil, three orange-pigmented yeast strains were isolated from samples collected in Minas Gerais state, Brazil. Molecular analyses combining the 26S rRNA gene (D1/D2 domains) and the internal transcribed spacer (ITS) sequences as well as whole-genome sequence data showed that these strains could not be ascribed to any known species in the basidiomycetous genus Phaffia, and thus they are considered to represent a novel species for which the name Phaffia brasiliana sp. nov. is proposed. The holotype is CBS 16121T and the MycoBank number is MB 839315. The occurrence of P. brasiliana in a tropical region is unique for the genus, since all other species occur in temperate regions. Two factors appear to contribute to the distribution of the novel taxon: first, the region where it was found has relatively moderate temperature ranges and, second, an adaptation to grow or withstand temperatures higher than those of the other species in the genus seems to be in place.

Multiple Pathways to Homothallism in Closely Related Yeast Lineages in the Basidiomycota.

Sexual reproduction in fungi relies on proteins with well-known functions encoded by the mating type (MAT) loci. In the Basidiomycota, MAT loci are often bipartite, with the P/R locus encoding pheromone precursors and pheromone receptors and the HD locus encoding heterodimerizing homeodomain transcription factors (Hd1/Hd2). The interplay between different alleles of these genes within a single species usually generates at least two compatible mating types. However, a minority of species are homothallic, reproducing sexually without an obligate need for a compatible partner. Here, we examine the organization and function of the MAT loci of Cystofilobasidium capitatum, a species in the order Cystofilobasidiales, which is unusually rich in homothallic species. We determined MAT gene content and organization in C. capitatum and found that it resembles a mating type of the closely related heterothallic species Cystofilobasidium ferigula To explain the homothallic sexual reproduction observed in C. capitatum, we examined HD protein interactions in the two Cystofilobasidium species and determined C. capitatum MAT gene expression both in a natural setting and upon heterologous expression in Phaffia rhodozyma, a homothallic species belonging to a clade sister to that of Cystofilobasidium. We conclude that the molecular basis for homothallism in C. capitatum appears to be distinct from that previously established for P. rhodozyma Unlike in the latter species, homothallism in C. capitatum may involve constitutive activation or dispensability of the pheromone receptor and the functional replacement of the usual Hd1/Hd2 heterodimer by an Hd2 homodimer. Overall, our results suggest that homothallism evolved multiple times within the Cystofilobasidiales.IMPORTANCE Sexual reproduction is important for the biology of eukaryotes because it strongly impacts the dynamics of genetic variation. In fungi, although sexual reproduction is usually associated with the fusion between cells belonging to different individuals (heterothallism), sometimes a single individual is capable of completing the sexual cycle alone (homothallism). Homothallic species are unusually common in a fungal lineage named Cystofilobasidiales. Here, we studied the genetic bases of homothallism in one species in this lineage, Cystofilobasidium capitatum, and found it to be different in several aspects from those of another homothallic species, Phaffia rhodozyma, belonging to the genus most closely related to Cystofilobasidium Our results strongly suggest that homothallism evolved independently in Phaffia and Cystofilobasidium, lending support to the idea that transitions between heterothallism and homothallism are not as infrequent as previously thought. Our work also helps to establish the Cystofilobasidiales as a model lineage in which to study these transitions.

The Untapped Australasian Diversity of Astaxanthin-Producing Yeasts with Biotechnological Potential-Phaffia australis sp. nov. and Phaffia tasmanica sp. nov.

Phaffia is an orange-colored basidiomycetous yeast genus of the order Cystofilobasidiales that contains a single species, P. rhodozyma. This species is the only fungus known to produce the economically relevant carotenoid astaxanthin. Although Phaffia was originally found in the Northern hemisphere, its diversity in the southern part of the globe has been shown to be much greater. Here we analyze the genomes of two Australasian lineages that are markedly distinct from P. rhodozyma. The two divergent lineages were investigated within a comprehensive phylogenomic study of representatives of the Cystofilobasidiales that supported the recognition of two novel Phaffia species, for which we propose the names of P. australis sp. nov. and P. tasmanica sp. nov. Comparative genomics and other analyses confirmed that the two new species have the typical Phaffia hallmark-the six genes necessary for the biosynthesis of astaxanthin could be retrieved from the draft genome sequences, and this carotenoid was detected in culture extracts. In addition, the organization of the mating-type (MAT) loci is similar to that of P. rhodozyma, with synteny throughout most regions. Moreover, cases of trans-specific polymorphism involving pheromone receptor genes and pheromone precursor proteins in the three Phaffia species, together with their shared homothallism, provide additional support for their classification in a single genus.

Revisiting the Taxonomic Synonyms and Populations of Saccharomyces cerevisiae-Phylogeny, Phenotypes, Ecology and Domestication.

Saccharomyces cerevisiae-the most emblematic and industrially relevant yeast-has a long list of taxonomical synonyms. Formerly considered as distinct species, some of the synonyms represent variants with important industrial implications, like Saccharomyces boulardii or Saccharomyces diastaticus, but with an unclear status, especially among the fermentation industry, the biotechnology community and biologists not informed on taxonomic matters. Here, we use genomics to investigate a group of 45 reference strains (type strains) of former Saccharomyces species that are currently regarded as conspecific with S. cerevisiae. We show that these variants are distributed across the phylogenetic spectrum of domesticated lineages of S. cerevisiae, with emphasis on the most relevant technological groups, but absent in wild lineages. We analyzed the phylogeny of a representative and well-balanced dataset of S. cerevisiae genomes that deepened our current ecological and biogeographic assessment of wild populations and allowed the distinction, among wild populations, of those associated with low- or high-sugar natural environments. Some wild lineages from China were merged with wild lineages from other regions in Asia and in the New World, thus giving more resolution to the current model of expansion from Asia to the rest of the world. We reassessed several key domestication markers among the different domesticated populations. In some cases, we could trace their origin to wild reservoirs, while in other cases gene inactivation associated with domestication was also found in wild populations, thus suggesting that natural adaptation to sugar-rich environments predated domestication.

The Wickerhamiella/Starmerella clade-A treasure trove for the study of the evolution of yeast metabolism.

The Wickerhamiella and Starmerella genera form a clade (W/S clade) that branches close to Yarrowia lipolytica in the Saccharomycotina species tree. It comprises approximately 90 recognized species and 50 putative new species not formally described yet. The large majority of the members of the W/S clade are ecologically associated with flowers and floricolous insects. Many species exhibit unusual metabolic traits, like fructophily and the production of sophorolipids, which are glycolipids that can be used as environmentally friendly biosurfactants. Genomic data have not only firmly established the W/S clade but have also revealed a tumultuous evolution of metabolism marked by losses and gains of important metabolic pathways, among which alcoholic fermentation. Possibly the most surprising finding brought to light by comparative genomics concerned the large number of genes acquired by some species of the W/S clade from bacteria through horizontal gene transfer, many of which were shown to be functional in their new setting. This was facilitated by the genetic tractability of one species in the clade, Starmerella bombicola, which is used for the industrial production of sophorolipids. We suggest that high-density coverage of genome sequencing in this clade, combined with the possibility to conduct molecular genetics experiments in at least one species, has the potential to set the stage for yet more exciting discoveries concerning the evolution of yeast metabolism.

Identification of novel autoinducer-2 receptors in Clostridia reveals plasticity in the binding site of the LsrB receptor family.

Autoinducer-2 (AI-2) is unique among quorum-sensing signaling molecules, as it is produced and recognized by a wide variety of bacteria and thus facilitates interspecies communication. To date, two classes of AI-2 receptors have been identified: the LuxP-type, present in the Vibrionales, and the LsrB-type, found in a number of phylogenetically distinct bacterial families. Recently, AI-2 was shown to affect the colonization levels of a variety of bacteria in the microbiome of the mouse gut, including members of the genus Clostridium, but no AI-2 receptor had been identified in this genus. Here, we identify a noncanonical, functional LsrB-type receptor in Clostridium saccharobutylicum. This novel LsrB-like receptor is the first one reported with variations in the binding-site amino acid residues that interact with AI-2. The crystal structure of the C. saccharobutylicum receptor determined at 1.35 Å resolution revealed that it binds the same form of AI-2 as the other known LsrB-type receptors, and isothermal titration calorimetry (ITC) assays showed that binding of AI-2 occurs at a submicromolar concentration. Using phylogenetic analysis, we inferred that the newly identified noncanonical LsrB receptor shares a common ancestor with known LsrB receptors and that noncanonical receptors are present in bacteria from different phyla. This led us to identify putative AI-2 receptors in bacterial species in which no receptors were known, as in bacteria belonging to the Spirochaetes and Actinobacteria phyla. Thus, this work represents a significant step toward understanding how AI-2-mediated quorum sensing influences bacterial interactions in complex biological niches.

Genetic Competence Drives Genome Diversity in Bacillus subtilis.

Prokaryote genomes are the result of a dynamic flux of genes, with increases achieved via horizontal gene transfer and reductions occurring through gene loss. The ecological and selective forces that drive this genomic flexibility vary across species. Bacillus subtilis is a naturally competent bacterium that occupies various environments, including plant-associated, soil, and marine niches, and the gut of both invertebrates and vertebrates. Here, we quantify the genomic diversity of B. subtilis and infer the genome dynamics that explain the high genetic and phenotypic diversity observed. Phylogenomic and comparative genomic analyses of 42 B. subtilis genomes uncover a remarkable genome diversity that translates into a core genome of 1,659 genes and an asymptotic pangenome growth rate of 57 new genes per new genome added. This diversity is due to a large proportion of low-frequency genes that are acquired from closely related species. We find no gene-loss bias among wild isolates, which explains why the cloud genome, 43% of the species pangenome, represents only a small proportion of each genome. We show that B. subtilis can acquire xenologous copies of core genes that propagate laterally among strains within a niche. While not excluding the contributions of other mechanisms, our results strongly suggest a process of gene acquisition that is largely driven by competence, where the long-term maintenance of acquired genes depends on local and global fitness effects. This competence-driven genomic diversity provides B. subtilis with its generalist character, enabling it to occupy a wide range of ecological niches and cycle through them.

Integrating Antimicrobial Therapy with Host Immunity to Fight Drug-Resistant Infections: Classical vs. Adaptive Treatment.

Antimicrobial resistance of infectious agents is a growing problem worldwide. To prevent the continuing selection and spread of drug resistance, rational design of antibiotic treatment is needed, and the question of aggressive vs. moderate therapies is currently heatedly debated. Host immunity is an important, but often-overlooked factor in the clearance of drug-resistant infections. In this work, we compare aggressive and moderate antibiotic treatment, accounting for host immunity effects. We use mathematical modelling of within-host infection dynamics to study the interplay between pathogen-dependent host immune responses and antibiotic treatment. We compare classical (fixed dose and duration) and adaptive (coupled to pathogen load) treatment regimes, exploring systematically infection outcomes such as time to clearance, immunopathology, host immunization, and selection of resistant bacteria. Our analysis and simulations uncover effective treatment strategies that promote synergy between the host immune system and the antimicrobial drug in clearing infection. Both in classical and adaptive treatment, we quantify how treatment timing and the strength of the immune response determine the success of moderate therapies. We explain key parameters and dimensions, where an adaptive regime differs from classical treatment, bringing new insight into the ongoing debate of resistance management. Emphasizing the sensitivity of treatment outcomes to the balance between external antibiotic intervention and endogenous natural defenses, our study calls for more empirical attention to host immunity processes.

Rethinking the Niche of Upper-Atmosphere Bacteria: Draft Genome Sequences of Bacillus aryabhattai C765 and Bacillus aerophilus C772, Isolated from Rice Fields.

Here, we report two genome sequences of endospore-forming bacteria isolated from the rice fields of Comporta, Portugal, identified as Bacillus aryabhattai C765 and Bacillus aerophilus C772. Both species were previously identified in air samples from the upper atmosphere, but our findings suggest their presence in a wider range of environmental niches.

Natural genome diversity of AI-2 quorum sensing in Escherichia coli: conserved signal production but labile signal reception.

Quorum sensing (QS) regulates the onset of bacterial social responses in function to cell density having an important impact in virulence. Autoinducer-2 (AI-2) is a signal that has the peculiarity of mediating both intra- and interspecies bacterial QS. We analyzed the diversity of all components of AI-2 QS across 44 complete genomes of Escherichia coli and Shigella strains. We used phylogenetic tools to study its evolution and determined the phenotypes of single-deletion mutants to predict phenotypes of natural strains. Our analysis revealed many likely adaptive polymorphisms both in gene content and in nucleotide sequence. We show that all natural strains possess the signal emitter (the luxS gene), but many lack a functional signal receptor (complete lsr operon) and the ability to regulate extracellular signal concentrations. This result is in striking contrast with the canonical species-specific QS systems where one often finds orphan receptors, without a cognate synthase, but not orphan emitters. Our analysis indicates that selection actively maintains a balanced polymorphism for the presence/absence of a functional lsr operon suggesting diversifying selection on the regulation of signal accumulation and recognition. These results can be explained either by niche-specific adaptation or by selection for a coercive behavior where signal-blind emitters benefit from forcing other individuals in the population to haste in cooperative behaviors.

Mutation accumulation in Tetrahymena.

BACKGROUND: The rate and fitness effects of mutations are key in understanding the evolution of every species. Traditionally, these parameters are estimated in mutation accumulation experiments where replicate lines are propagated in conditions that allow mutations to randomly accumulate without the purging effect of natural selection. These experiments have been performed with many model organisms but we still lack empirical estimates of the rate and effects of mutation in the protists. RESULTS: We performed a mutation accumulation (MA) experiment in Tetrahymena thermophila, a species that can reproduce sexually and asexually in nature, and measured both the mean decline and variance increase in fitness of 20 lines. The results obtained with T. thermophila were compared with T. pyriformis that is an obligate asexual species. We show that MA lines of T. thermophila go to extinction at a rate of 1.25 clonal extinctions per bottleneck. In contrast, populations of T. pyriformis show a much higher resistance to extinction. Variation in gene copy number is likely to be a key factor in explaining these results, and indeed we show that T. pyriformis has a higher mean copy number per cell than T. thermophila. From fitness measurements during the MA experiment, we infer a rate of mutation to copy number variation of 0.0333 per haploid MAC genome of T. thermophila and a mean effect against copy number variation of 0.16. A strong effect of population size in the rate of fitness decline was also found, consistent with the increased power of natural selection. CONCLUSIONS: The rate of clonal extinction measured for T. thermophila is characteristic of a mutational degradation and suggests that this species must undergo sexual reproduction to avoid the deleterious effects detected in the laboratory experiments. We also suggest that an increase in chromosomal copy number associated with the phenotypic assortment of amitotic divisions can provide an alternative mechanism to escape the deleterious effect of random chromosomal copy number variation in species like T. pyriformis that lack the resetting mechanism of sexual reproduction. Our results are relevant to the understanding of cell line longevity and senescence in ciliates.

Identification of functional LsrB-like autoinducer-2 receptors.

Although a variety of bacterial species have been reported to use the interspecies communication signal autoinducer-2 (AI-2) to regulate multiple behaviors, the molecular mechanisms of AI-2 recognition and signal transduction remain poorly understood. To date, two types of AI-2 receptors have been identified: LuxP, present in Vibrio spp., and LsrB, first identified in Salmonella enterica serovar Typhimurium. In S. Typhimurium, LsrB is the ligand binding protein of a transport system that enables the internalization of AI-2. Here, using both sequence analysis and structure prediction, we establish a set of criteria for identifying functional AI-2 receptors. We test our predictions experimentally, assaying key species for their abilities to import AI-2 in vivo, and test their LsrB orthologs for AI-2 binding in vitro. Using these experimental approaches, we were able to identify AI-2 receptors in organisms belonging to phylogenetically distinct families such as the Enterobacteriaceae, Rhizobiaceae, and Bacillaceae. Phylogenetic analysis of LsrB orthologs indicates that this pattern could result from one single origin of the functional LsrB gene in a gammaproteobacterium, suggesting possible posterior independent events of lateral gene transfer to the Alphaproteobacteria and Firmicutes. Finally, we used mutagenesis to show that two AI-2-interacting residues are essential for the AI-2 binding ability. These two residues are conserved in the binding sites of all the functional AI-2 binding proteins but not in the non-AI-2-binding orthologs. Together, these results strongly support our ability to identify functional LsrB-type AI-2 receptors, an important step in investigations of this interspecies signal.

Multilocus phylogeography and phylogenetics using sequence-based markers.

We review recent trends in phylogeography and phylogenetics and argue that these two fields stand to be reunited by the common yardstick provided by sequence and SNP data and by new multilocus methods for phylogenetic analysis. Whereas the modern incarnation of both fields was spawned by PCR approaches applied to mitochondrial DNA in the late 1980s, the two fields diverged during the 1990s largely due to the adoption by phylogeographers of microsatellites, in contrast to the adoption of nuclear sequence data by phylogeneticists. Sequence-based markers possess a number of advantages over microsatellites, even on the recent time scales that are the purview of phylogeography. Using examples primarily from vertebrates, we trace the maturation of nuclear gene phylogeography and phylogenetics and suggest that the abundant instances of gene tree heterogeneity beckon a new generation of phylogenetic methods that focus on estimating species trees as distinct from gene trees. Whole genomes provide a powerful common yardstick on which both phylogeography and phylogenetics can assume their proper place as ends of a continuum.

Contrasting patterns of mitochondrial and microsatellite genetic structure among Western European populations of tawny owls (Strix aluco).

A recent study of mitochondrial phylogeography of tawny owls (Strix aluco) in western Europe suggested that this species survived the Pleistocene glaciations in three allopatric refugia located in Iberia, Italy, and the Balkans, and the latter was likely the predominant source of postglacial colonization of northern Europe. New data from seven microsatellite loci from 184 individual owls distributed among 14 populations were used to assess the genetic congruence between nuclear and mitochondrial DNA (mtDNA) markers. Microsatellites corroborated the major phylogeographical conclusions reached on the basis of the mtDNA sequences, but also showed important differences leading to novel inferences. Microsatellites corroborated the three major refugia and supported the Balkan origin of northern populations. When corrected for differences in effective population size, microsatellites and mtDNA yielded generally congruent overall estimates of population structure (N*ST=0.12 vs. RST=0.16); however, there was substantial heterogeneity in the RST among the seven nuclear loci that was not correlated with heterozygosity. Populations representing the Balkans postglacial expansion interact with populations from the other two refugia forming two clines near the Alps and the Pyrenees. In both cases, the apparent position of the contact zones differed substantially between markers due to the genetic composition of populations sampled in northern Italy and Madrid. Microsatellite data did not corroborate the lower genetic diversity of northern, recently populated regions as was found with mtDNA; this discrepancy was taken as evidence for a recent bottleneck recovery. Finally, this study suggests that congruence among genetic markers should be more likely in cases of range expansion into new areas than when populations interact across contact zones.

The influence of Pleistocene glacial refugia on tawny owl genetic diversity and phylogeography in western Europe.

The glacial refugia hypothesis indicates that during the height of the Pleistocene glaciations the temperate species that are today widespread in western Europe must have survived in small and climatically favourable areas located in the southern peninsulas of Iberia, Italy and Balkans. One such species is the tawny owl, a relatively sedentary, nonmigratory bird presently distributed throughout Europe. It is a tree-nesting species closely associated with deciduous and mixed coniferous woodlands. In this study I used control region mtDNA sequences from 187 individuals distributed among 14 populations to determine whether current genetic patterns in tawny owl populations were consistent with postglacial expansion from peninsular refugia. European, North African and Asian tawny owls were found to represent three distinct lineages, where North Africa is the sister clade to all European owls. Within Europe, I found three well-supported clades that correspond to each of the three allopatric refugia. Expansion patterns indicate that owls from the Balkan refugium repopulated most of northern Europe, while expansion out of Iberia and Italy had only regional effects leading to admixture in France. Estimates of population divergence times between refugia populations are roughly similar, but one order of magnitude smaller between Greece and northern Europe. Based on a wide range of mutation rates and generation times, divergence between refugia appears to date to the Pleistocene.