Low reproductive isolation and highly variable levels of gene flow reveal limited progress towards speciation between European river and brook lampreys

Ecologically based divergent selection is a factor that could drive reproductive isolation even in the presence of gene flow. Population pairs arrayed along a continuum of divergence provide a good opportunity to address this issue. Here, we used a combination of mating trials, experimental crosses and population genetic analyses to investigate the evolution of reproductive isolation between two closely related species of lampreys with distinct life histories. We used microsatellite markers to genotype over 1000 individuals of the migratory parasitic river lamprey (Lampetra fluviatilis) and freshwater‐resident nonparasitic brook lamprey (Lampetra planeri) distributed in 10 sympatric and parapatric population pairs in France. Mating trials, parentage analyses and artificial fertilizations demonstrated a low level of reproductive isolation between species even though size‐assortative mating may contribute to isolation. Most parapatric population pairs were strongly differentiated due to the joint effects of geographic distance and barriers to migration. In contrast, we found variable levels of gene flow between sympatric populations ranging from panmixia to moderate differentiation, which indicates a gradient of divergence with some population pairs that may correspond to alternative morphs or ecotypes of a single species and others that remain partially isolated. Ecologically based divergent selection may explain these variable levels of divergence among sympatric population pairs, but incomplete genome swamping following secondary contact could have also played a role. Overall, this study illustrates how highly differentiated phenotypes can be maintained despite high levels of gene flow that limit the progress towards speciation.     

Phylogeography: spanning the ecology‐evolution continuum

Synthesis of ecological and evolutionary concepts and tools has led to improved understanding of how diversification, dispersal, community assembly, long‐term coexistence and extinction shape patterns of biological diversity. Phylogeography, with its focus on Quaternary interactions within and between populations, can help elucidate the processes acting between the evolutionary time‐scales on which species arise and the ecological time‐scales on which members of an assemblage interact with each other and their environment. Still, it has yet to be widely incorporated in that synthesis. Here, we highlight three areas where integration of phylogeography with ecological and evolutionary approaches can provide new insights into key questions. First, phylogeography can help clarify the roles of isolation, niche conservatism and environmental stability in generating patterns of alpha‐ and beta‐diversity. Second, phylogeography can help isolate the effects of Quaternary dispersal limitation from other factors driving community assembly and spatial turnover. Third, phylogeography can help identify key processes leading to and resulting from extinction events, including the population dynamics of species range reduction and its effects on the strength and temporal flexibility of networks of species interactions. We conclude with an outlook on the data‐gathering protocols necessary for this collaborative, interdisciplinary research agenda.     

Birds are islands for parasites

Understanding the mechanisms driving the extraordinary diversification of parasites is a major challenge in evolutionary biology. Co-speciation, one proposed mechanism that could contribute to this diversity is hypothesized to result from allopatric co-divergence of host–parasite populations. We found that island populations of the Galápagos hawk (Buteo galapagoensis) and a parasitic feather louse species (Degeeriella regalis) exhibit patterns of co-divergence across variable temporal and spatial scales. Hawks and lice showed nearly identical population genetic structure across the Galápagos Islands. Hawk population genetic structure is explained by isolation by distance among islands. Louse population structure is best explained by hawk population structure, rather than isolation by distance per se, suggesting that lice tightly track the recent population histories of their hosts. Among hawk individuals, louse populations were also highly structured, suggesting that hosts serve as islands for parasites from an evolutionary perspective. Altogether, we found that host and parasite populations may have responded in the same manner to geographical isolation across spatial scales. Allopatric co-divergence is likely one important mechanism driving the diversification of parasites.     

The aggregate site frequency spectrum for comparative population genomic inference

Understanding how assemblages of species responded to past climate change is a central goal of comparative phylogeography and comparative population genomics, an endeavour that has increasing potential to integrate with community ecology. New sequencing technology now provides the potential to perform complex demographic inference at unprecedented resolution across assemblages of nonmodel species. To this end, we introduce the aggregate site frequency spectrum (aSFS), an expansion of the site frequency spectrum to use single nucleotide polymorphism (SNP) data sets collected from multiple, co‐distributed species for assemblage‐level demographic inference. We describe how the aSFS is constructed over an arbitrary number of independent population samples and then demonstrate how the aSFS can differentiate various multispecies demographic histories under a wide range of sampling configurations while allowing effective population sizes and expansion magnitudes to vary independently. We subsequently couple the aSFS with a hierarchical approximate Bayesian computation (hABC) framework to estimate degree of temporal synchronicity in expansion times across taxa, including an empirical demonstration with a data set consisting of five populations of the threespine stickleback (Gasterosteus aculeatus). Corroborating what is generally understood about the recent postglacial origins of these populations, the joint aSFS/hABC analysis strongly suggests that the stickleback data are most consistent with synchronous expansion after the Last Glacial Maximum (posterior probability = 0.99). The aSFS will have general application for multilevel statistical frameworks to test models involving assemblages and/or communities, and as large‐scale SNP data from nonmodel species become routine, the aSFS expands the potential for powerful next‐generation comparative population genomic inference.     

The cophylogeny of populations and cultures: reconstructing the evolution of Iranian tribal craft traditions using trees and jungles

Phylogenetic approaches to culture have shed new light on the role played by population dispersals in the spread and diversification of cultural traditions. However, the fact that cultural inheritance is based on separate mechanisms from genetic inheritance means that socially transmitted traditions have the potential to diverge from population histories. Here, we suggest that associations between these two systems can be reconstructed using techniques developed to study cospeciation between hosts and parasites and related problems in biology. Relationships among the latter are patterned by four main processes: co-divergence, intra-host speciation (duplication), intra-host extinction (sorting) and horizontal transfers. We show that patterns of cultural inheritance are structured by analogous processes, and then demonstrate the applicability of the host-parasite model to culture using empirical data on Iranian tribal populations.     

分子系统地理学及其应用

分子系统地理学是20世纪70年代中期伴随着对线粒体DNA的认识而开始酝酿发展的,本文回顾了该学科自诞生以来的发展简史,阐述了其研究内容,着重从3个方面介绍了该领域在其它相关学科中的研究进展,总结了该学科存在的问题,并提出重点应解决的3个问题：1）将细胞器基因与遗传信息更为丰富的核基因相结合;2）合生理论在非平衡种群中的发展应用;3）提高对分化时间估计的精确性。     

60 million years of co-divergence in the fig-wasp symbiosis

Figs (Ficus; ca 750 species) and fig wasps (Agaoninae) are obligate mutualists: all figs are pollinated by agaonines that feed exclusively on figs. This extraordinary symbiosis is the most extreme example of specialization in a plant-pollinator interaction and has fuelled much speculation about co-divergence. The hypothesis that pollinator specialization led to the parallel diversification of fig and pollinator lineages (co-divergence) has so far not been tested due to the lack of robust and comprehensive phylogenetic hypotheses for both partners. We produced and combined the most comprehensive molecular phylogenetic trees to date with fossil data to generate independent age estimates for fig and pollinator lineages, using both non-parametric rate smoothing and penalized likelihood dating methods. Molecular dating of ten pairs of interacting lineages provides an unparalleled example of plant-insect co-divergence over a geological time frame spanning at least 60 million years.     

Where and when does a ring start and end? Testing the ring-species hypothesis in a species complex of Australian parrots

Speciation, despite ongoing gene flow can be studied directly in nature in ring species that comprise two reproductively isolated populations connected by a chain or ring of intergrading populations. We applied three tiers of spatio-temporal analysis (phylogeny/historical biogeography, phylogeography and landscape/population genetics) to the data from mitochondrial and nuclear genomes of eastern Australian parrots of the Crimson Rosella Platycercus elegans complex to understand the history and present genetic structure of the ring they have long been considered to form. A ring speciation hypothesis does not explain the patterns we have observed in our data (e.g. multiple genetic discontinuities, discordance in genotypic and phenotypic assignments where terminal differentiates meet). However, we cannot reject that a continuous circular distribution has been involved in the group's history or indeed that one was formed through secondary contact at the 'ring's' east and west; however, we reject a simple ring-species hypothesis as traditionally applied, with secondary contact only at its east. We discuss alternative models involving historical allopatry of populations. We suggest that population expansion shown by population genetics parameters in one of these isolates was accompanied by geographical range expansion, secondary contact and hybridization on the eastern and western sides of the ring. Pleistocene landscape and sea-level and habitat changes then established the birds' current distributions and range disjunctions. Populations now show idiosyncratic patterns of selection and drift. We suggest that selection and drift now drive evolution in different populations within what has been considered the ring.     

Estimating Variable Effective Population Sizes from Multiple Genomes: A Sequentially Markov Conditional Sampling Distribution Approach

Throughout history, the population size of modern humans has varied considerably due to changes in environment, culture, and technology. More accurate estimates of population size changes, and when they occurred, should provide a clearer picture of human colonization history and help remove confounding effects from natural selection inference. Demography influences the pattern of genetic variation in a population, and thus genomic data of multiple individuals sampled from one or more present-day populations contain valuable information about the past demographic history. Recently, Li and Durbin developed a coalescent-based hidden Markov model, called the pairwise sequentially Markovian coalescent (PSMC), for a pair of chromosomes (or one diploid individual) to estimate past population sizes. This is an efficient, useful approach, but its accuracy in the very recent past is hampered by the fact that, because of the small sample size, only few coalescence events occur in that period. Multiple genomes from the same population contain more information about the recent past, but are also more computationally challenging to study jointly in a coalescent framework. Here, we present a new coalescent-based method that can efficiently infer population size changes from multiple genomes, providing access to a new store of information about the recent past. Our work generalizes the recently developed sequentially Markov conditional sampling distribution framework, which provides an accurate approximation of the probability of observing a newly sampled haplotype given a set of previously sampled haplotypes. Simulation results demonstrate that we can accurately reconstruct the true population histories, with a significant improvement over the PSMC in the recent past. We apply our method, called diCal, to the genomes of multiple human individuals of European and African ancestry to obtain a detailed population size change history during recent times.     

Genealogy-Based Methods for Inference of Historical Recombination and Gene Flow and Their Application in Saccharomyces cerevisiae

Genetic exchange between isolated populations, or introgression between species, serves as a key source of novel genetic material on which natural selection can act. While detecting historical gene flow from DNA sequence data is of much interest, many existing methods can be limited by requirements for deep population genomic sampling. In this paper, we develop a scalable genealogy-based method to detect candidate signatures of gene flow into a given population when the source of the alleles is unknown. Our method does not require sequenced samples from the source population, provided that the alleles have not reached fixation in the sampled recipient population. The method utilizes recent advances in algorithms for the efficient reconstruction of ancestral recombination graphs, which encode genealogical histories of DNA sequence data at each site, and is capable of detecting the signatures of gene flow whose footprints are of length up to single genes. Further, we employ a theoretical framework based on coalescent theory to test for statistical significance of certain recombination patterns consistent with gene flow from divergent sources. Implementing these methods for application to whole-genome sequences of environmental yeast isolates, we illustrate the power of our approach to highlight loci with unusual recombination histories. By developing innovative theory and methods to analyze signatures of gene flow from population sequence data, our work establishes a foundation for the continued study of introgression and its evolutionary relevance.     

multi‐dice: r package for comparative population genomic inference under hierarchical co‐demographic models of independent single‐population size changes

Population genetic data from multiple taxa can address comparative phylogeographic questions about community‐scale response to environmental shifts, and a useful strategy to this end is to employ hierarchical co‐demographic models that directly test multi‐taxa hypotheses within a single, unified analysis. This approach has been applied to classical phylogeographic data sets such as mitochondrial barcodes as well as reduced‐genome polymorphism data sets that can yield 10,000s of SNPs, produced by emergent technologies such as RAD‐seq and GBS. A strategy for the latter had been accomplished by adapting the site frequency spectrum to a novel summarization of population genomic data across multiple taxa called the aggregate site frequency spectrum (aSFS), which potentially can be deployed under various inferential frameworks including approximate Bayesian computation, random forest and composite likelihood optimization. Here, we introduce the r package multi‐dice , a wrapper program that exploits existing simulation software for flexible execution of hierarchical model‐based inference using the aSFS, which is derived from reduced genome data, as well as mitochondrial data. We validate several novel software features such as applying alternative inferential frameworks, enforcing a minimal threshold of time surrounding co‐demographic pulses and specifying flexible hyperprior distributions. In sum, multi‐dice provides comparative analysis within the familiar R environment while allowing a high degree of user customization, and will thus serve as a tool for comparative phylogeography and population genomics.     

Differences in Quaternary co-divergence reveals community-wide diversification in the mountains of southwest China varied among species

The mountains of southwest China (MSWC) is a biodiversity hotspot with highly complex and unusual terrain. However, with the majority of studies focusing on the biogeographic consequences of massive mountain building, the Quaternary legacy of biodiversity for the MSWC has long been overlooked. Here, we took a statistical comparative phylogeography approach to examine factors that shaped community-wide diversification. With data from 30 vertebrate species, the results reveal spatially concordant genetic structure, and temporally clustered co-divergence events associated with river barriers during severe glacial cycles. This indicates the importance of riverine barriers in the phylogeographic history of the MSWC vertebrate community. We conclude that the repeated glacial cycles are associated with co-divergences that are themselves structured by the heterogeneity of the montane landscape of the MSWC. This orderly process of diversification has profound implications for conservation by highlighting the relative independence of different geographical areas in which some, but not all species in communities have responded similarly to climate change and calls for further comparative phylogeographic investigations to reveal the connection between biological traits and divergence pulses in this biodiversity hotspot.     

The genetic structure of the introduced house sparrow populations in Australia and New Zealand is consistent with historical descriptions of multiple introductions to each country

The house sparrow is one of the most widely introduced vertebrate species around the world, making it an important model species for the study of invasion ecology. Population genetic studies of these invasions provide important insights into colonisation processes and adaptive responses occurring during invasion. Here we use microsatellite data to infer the population structure and invasion history of the introduced house sparrow (Passer domesticus) in Australia and New Zealand. Our results identify stronger population structure within Australia in comparison to New Zealand and patterns are consistent with historical records of multiple introduction sites across both countries. Within the five population clusters identified in Australia, we find declines in genetic diversity as we move away from the reported introduction site within each cluster. This pattern is consistent with sequential founder events. Interestingly, an even stronger decline in genetic diversity is seen across Australia as we move away from the Melbourne introduction site; secondary historical reports suggest this site imported a large number of sparrows and was possible the source of a single range expansion across Australia. However, private allele numbers are highest in the north, away from Melbourne, which could be a result of drift increasing the frequency of rare alleles in areas of smaller population size or due to an independent introduction that seeded or augmented the northern population. This study highlights the difficulties of elucidating population dynamics in introduced species with complex introduction histories and suggests that a combination of historical and genetic data can be useful.     

Glacial history of the North Atlantic marine snail, Littorina saxatilis, inferred from distribution of mitochondrial DNA lineages

The North Atlantic intertidal gastropod, Littorina saxatilis (Olivi, 1792), exhibits extreme morphological variation between and within geographic regions and has become a model for studies of local adaptation; yet a comprehensive analysis of the species' phylogeography is lacking. Here, we examine phylogeographic patterns of the species' populations in the North Atlantic and one remote Mediterranean population using sequence variation in a fragment of the mitochondrial cytochrome b gene (607 bp). We found that, as opposed to many other rocky intertidal species, L. saxatilis has likely had a long and continuous history in the Northwest Atlantic, including survival during the last glacial maximum (LGM), possibly in two refugia. In the Northeast Atlantic, several areas likely harboured refugial populations that recolonized different parts of this region after glacial retreat, resulting in strong population structure. However, the outlying monomorphic Venetian population is likely a recent anthropogenic introduction from northern Europe and not a remnant of an earlier wider distribution in the Mediterranean Sea. Overall, our detailed phylogeography of L. saxatilis adds an important piece to the understanding of Pleistocene history in North Atlantic marine biota as well as being the first study to describe the species' evolutionary history in its natural range. The latter contribution is noteworthy because the snail has recently become an important model species for understanding evolutionary processes of speciation; thus our work provides integral information for such endeavours.     

Cyclostome embryology and early evolutionary history of vertebrates

Modern agnathans include only two groups, the lampreys and thehagfish, that collectively comprise the group Cyclostomata.Although accumulating molecular data support the cyclostomesas a monophyletic group, there remain some unsettled questionsregarding the evolutionary relationships of these animals inthat they differ greatly in anatomical and developmental patternsand in their life histories. In this review, we summarize recentdevelopmental data on the lamprey and discuss some questionsrelated to vertebrate evolutionary development raised by thelimited information available on hagfish embryos. Comparisonof the lamprey and gnathostome developmental patterns suggestssome plesiomorphic traits of vertebrates that would have alreadybeen established in the most recent common ancestor of the vertebrates.Understanding hagfish development will further clarify the,as yet, unrecognized ancestral characters that either the lampreysor hagfishes may have lost. We stress the immediate importanceof hagfish embryology in the determination of the most plausiblescenario for the early history of vertebrate evolution, by addressingquestions about the origins of the neural crest, thyroid, andadenohypophysis as examples.     

Cattle demographic history modelled from autosomal sequence variation

The phylogeography of cattle genetic variants has been extensively described and has informed the history of domestication. However, there remains a dearth of demographic models inferred from such data. Here, we describe sequence diversity at 37 000 bp sampled from 17 genes in cattle from Africa, Europe and India. Clearly distinct population histories are suggested between Bos indicus and Bos taurus, with the former displaying higher diversity statistics. We compare the unfolded site frequency spectra in each to those simulated using a diffusion approximation method and build a best-fitting model of past demography. This implies an earlier, possibly glaciation-induced population bottleneck in B. taurus ancestry with a later, possibly domestication-associated demographic constriction in B. indicus. Strikingly, the modelled indicine history also requires a majority secondary admixture from the South Asian aurochs, indicating a complex, more diffuse domestication process. This perhaps involved multiple domestications and/or introgression from wild oxen to domestic herds; the latter is plausible from archaeological evidence of contemporaneous wild and domestic remains across different regions of South Asia.     

The Ubiquitin-conjugating Enzyme (E2) Ube2w Ubiquitinates the N Terminus of Substrates

Attachment of ubiquitin to substrate is typically thought to occur via formation of an isopeptide bond between the C-terminal glycine residue of ubiquitin and a lysine residue in the substrate. In vitro, Ube2w is nonreactive with free lysine yet readily ubiquitinates substrate. Ube2w also contains novel residues within its active site that are important for its ability to ubiquitinate substrate. To identify the site of modification, we analyzed ubiquitinated substrates by mass spectrometry and found the N-terminal -NH₂ group as the site of conjugation. To confirm N-terminal ubiquitination, we generated lysine-less and N-terminally blocked versions of one substrate, the polyglutamine disease protein ataxin-3, and showed that Ube2w can ubiquitinate a lysine-less, but not N-terminally blocked, ataxin-3. This was confirmed with a second substrate, the neurodegenerative disease protein Tau. Finally, we directly sequenced the N terminus of unmodified and ubiquitinated ataxin-3, demonstrating that Ube2w attaches ubiquitin to the N terminus of its substrates. Together these data demonstrate that Ube2w has novel enzymatic properties that direct ubiquitination of the N terminus of substrates.     

Spatial but not temporal co-divergence of a virus and its mammalian host

Co-divergence between host and parasites suggests that evolutionary processes act across similar spatial and temporal scales. Although there has been considerable work on the extent and correlates of co-divergence of RNA viruses and their mammalian hosts, relatively little is known about the extent to which virus evolution is determined by the phylogeographic history of host species. To test hypotheses related to co-divergence across a variety of spatial and temporal scales, we explored phylogenetic signatures in Andes virus (ANDV) sampled from Chile and its host rodent, Oligoryzomys longicaudatus. ANDV showed strong spatial subdivision, a phylogeographic pattern also recovered in the host using both spatial and genealogical approaches, and despite incomplete lineage sorting. Lineage structure in the virus seemed to be a response to current population dynamics in the host at the spatial scale of ecoregions. However, finer scale analyses revealed contrasting patterns of genetic structure across a latitudinal gradient. As predicted by their higher substitution rates, ANDV showed greater genealogical resolution than the rodent, with topological congruence influenced by the degree of lineage sorting within the host. However, despite these major differences in evolutionary dynamics, the geographic structure of host and virus converged across large spatial scales.     

In vitro and in vivo effects of GABA, muscimol, and bicuculline on lamprey GnRH concentration in the brain of the sea lamprey (Petromyzon marinus)

gamma-Aminobutyric acid (GABA) is a neurotransmitter with a demonstrated neuroregulatory role in reproduction in most representative species of vertebrate classes via the hypothalamus. The role of GABA on the hypothalamus-pituitary axis in lampreys has not been fully elucidated. Recent immunocytochemical and in situ hybridization studies suggest that there may be a neuroregulatory role of GABA on the gonadotropin-releasing hormone (GnRH) system in lampreys. To assess possible GABA-GnRH interactions, the effects of GABA and its analogs on lamprey GnRH in vitro and in vivo were studied in adult female sea lampreys (Petromyzon marinus). In vitro perfusion of GABA and its analogs at increasing concentrations (0.1-100 microM) was performed over a 3-h time course. There was a substantial increase of GnRH-I and GnRH-III following treatment of muscimol at 100 microM. In in vivo studies, GABA or muscimol injected at 200 microg/kg significantly increased lamprey GnRH concentration in the brain 0.5 h after treatment compared to controls in female sea lampreys. No significant change in lamprey GnRH-I or GnRH-III was observed following treatment with bicuculline. These data provide novel physiological data supporting the hypothesis that GABA may influence GnRH in the brain of sea lamprey.