The study of structured populations--new hope for a difficult and divided science

Population genetic structure and historical demography of the population of forest elephants in Côte d'Ivoire

The population of forest elephant (Loxodonta cyclotis) has continuously declined in Côte d'Ivoire and, the remaining population largely consists of subpopulations that are fragmented and isolated. No data actually exist on the level of genetic diversity and population genetic structure of current forest elephant populations in Côte d'Ivoire. In this sense, determining genetic diversity and the underlying mechanisms of population differentiation is crucial for the initiation of effective conservation management. A total of 158 dung samples of forest elephants were collected at stage 1 of decompositions (dung pile intact, very fresh) in three Classified Forests (CF) (Bossématié, Dassioko and Port-Gauthier) in Côte d'Ivoire. A total of 101 sequences of the mitochondrial DNA control region measuring 600 base pair and 26 haplotypes were obtained. A haplotypic diversity ranging from 0.655 ± 0.050 at Bossématié and 0.859 ± 0.088 at Port Gauthier was obtained. Fifteen (15) out of 26 haplotypes observed were singletons and only the Dassioko and Port Gauthier CFs shared the same haplotypes. The strong genetic connectivity between forest elephant populations of the Dassioko and Port Gauthier CFs is supported by the grouping of these populations into a single cluster by Bayesian analysis. Although populations of L. cyclotis exhibit relatively high genetic diversity, habitat fragmentation could affect the genetic variability of current populations. Urgent measures including the reinforcement/establishment of genetic corridors and the strengthening of protection measures need to be undertaken to save the remaining populations of forest elephants in Côte d'Ivoire.

Genetic signatures of lineage fusion closely resemble population decline

Accurate interpretation of the genetic signatures of past demographic events is crucial for reconstructing evolutionary history. Lineage fusion (complete merging, resulting in a single panmictic population) is a special case of secondary contact that is seldom considered. Here, the circumstances under which lineage fusion can be distinguished from population size constancy, growth, bottleneck, and decline were investigated. Multi-locus haplotype data were simulated under models of lineage fusion with different divergence versus sampling lag times (D:L ratios). These pseudo-observed datasets also differed in their allocation of a fixed amount of sequencing resources (number of sampled alleles, haplotype length, number of loci). Distinguishability of lineage fusion versus each of 10 untrue non-fusion scenarios was quantified based on six summary statistics (neutrality tests). Some datasets were also analyzed using extended Bayesian skyline plots. Results showed that signatures of lineage fusion very closely resemble those of decline-high distinguishability was generally limited to the most favorable scenario (D:L = 9), using the most sensitive summary statistics (F S and Z nS), coupled with the optimal sequencing resource allocation (maximizing number of loci). Also, extended Bayesian skyline plots often erroneously inferred population decline. Awareness of the potential for lineage fusion to carry the hallmarks of population decline is critical.

Demographically explicit scans for barriers to gene flow using gIMble

Identifying regions of the genome that act as barriers to gene flow between recently diverged taxa has remained challenging given the many evolutionary forces that generate variation in genetic diversity and divergence along the genome, and the stochastic nature of this variation. Progress has been impeded by a conceptual and methodological divide between analyses that infer the demographic history of speciation and genome scans aimed at identifying locally maladaptive alleles i.e. genomic barriers to gene flow. Here we implement genomewide IM blockwise likelihood estimation (gIMble), a composite likelihood approach for the quantification of barriers, that bridges this divide. This analytic framework captures background selection and selection against barriers in a model of isolation with migration (IM) as heterogeneity in effective population size (Ne) and effective migration rate (me), respectively. Variation in both effective demographic parameters is estimated in sliding windows via pre-computed likelihood grids. gIMble includes modules for pre-processing/filtering of genomic data and performing parametric bootstraps using coalescent simulations. To demonstrate the new approach, we analyse data from a well-studied pair of sister species of tropical butterflies with a known history of post-divergence gene flow: Heliconius melpomene and H. cydno. Our analyses uncover both large-effect barrier loci (including well-known wing-pattern genes) and a genome-wide signal of a polygenic barrier architecture.

Genetic population structure of Pseudoplatystoma corruscans (Siluriformes: Pimelodidae) and evidence of temporal variation in structure

Surubim (Pseudoplatystoma corruscans, Pimelodidae) are migratory catfish native to the rivers in the La Plata and São Francisco basins. They are piscivores that attain considerable body sizes and are a valuable economic resource. Surubim exhibits extensive migrations during its life cycle that may affect the population structure at vast geographic scales. The authors examined the genetic diversity and population genetic structure of P. corruscans using microsatellite markers from a comprehensive sampling of 260 individuals from the Upper and Lower Paraná River. They identified two well-differentiated genetic clusters corresponding to a natural geographic barrier historically separating Upper and Lower Paraná regions. They also demonstrated temporal variation in population genetic structure at a site in Lower Paraná close to the confluence with the Paraguay River, most likely explained by the influx of migrant fishes at certain times of the year.

The persistent homology of genealogical networks

Genealogical networks (i.e. family trees) are of growing interest, with the largest known data sets now including well over one billion individuals. Interest in family history also supports an 8.5 billion dollar industry whose size is projected to double within 7 years [FutureWise report HC-1137]. Yet little mathematical attention has been paid to the complex network properties of genealogical networks, especially at large scales. The structure of genealogical networks is of particular interest due to the practice of forming unions, e.g. marriages, that are typically well outside one's immediate family. In most other networks, including other social networks, no equivalent restriction exists on the distance at which relationships form. To study the effect this has on genealogical networks we use persistent homology to identify and compare the structure of 101 genealogical and 31 other social networks. Specifically, we introduce the notion of a network's persistence curve, which encodes the network's set of persistence intervals. We find that the persistence curves of genealogical networks have a distinct structure when compared to other social networks. This difference in structure also extends to subnetworks of genealogical and social networks suggesting that, even with incomplete data, persistent homology can be used to meaningfully analyze genealogical networks. Here we also describe how concepts from genealogical networks, such as common ancestor cycles, are represented using persistent homology. We expect that persistent homology tools will become increasingly important in genealogical exploration as popular interest in ancestry research continues to expand.

Multilocus phylogeography, population genetics and niche evolution of Australian brown and black-tailed treecreepers (Aves: Climacteris)

The Carpentarian barrier across north-eastern Australia is a major biogeographic barrier and a generator of biodiversity within the Australian Monsoonal Tropics. Here we present a continent-wide analysis of mitochondrial (control region) and autosomal (14 anonymous loci) sequence and indel variation and niche modelling of brown and black-tailed treecreepers (Climacteris picumnus and Climacteris melanurus), a clade with a classic distribution on either side of the Carpentarian barrier. mtDNA control region sequences exhibited reciprocal monophyly and strong differentiation (Fst = 0.91), and revealed a signature of a recent selective sweep in C. picumnus. A variety of tests support an isolation-with-migration model of divergence, albeit with low levels of gene flow across the Carpentarian barrier and a divergence time between species of ~1.7–2.8 Mya. Palaeoecological niche models show that both range size as measured by available habitat and estimated historical population sizes of both species declined in the past ~600 kyr and that the area of interspecific range overlap was never historically large, perhaps decreasing opportunities for extensive gene flow. The relatively long divergence time and low opportunity for gene flow may have facilitated speciation more so than in other co-distributed bird taxa across the Australian Monsoonal Tropics.

Phylogenetic relationships of sleeper gobies (Eleotridae: Gobiiformes: Gobioidei), with comments on the position of the miniature genus Microphilypnus

Microphilypnus and Leptophilypnion are miniaturized genera within the family Eleotridae. The evolutionary relationships among these taxa are still poorly understood, and molecular analyses are restricted to mitochondrial genes, which have not been conclusive. We compiled both mitochondrial and nuclear genes to study the phylogenetic position of Microphilypnus and the evolutionary history and relationships of eleotrids. We propose that Microphilypnus and Leptophilypnus (a non-miniature genus) are not sister groups as suggested by previous studies, but rather separate lineages that arose in the early Eocene, with Leptophilypnus recovered as a sister group to the other analyzed eleotrids. In fact, Microphilypnus is currently associated with the Neotropical clade Guavina/Dormitator/Gobiomorus. We also identified a well-supported clade that indicated Gobiomorus and Hemieleotris as paraphyletic groups, besides a close relationship among Calumia godeffroyi, Bunaka gyrinoides, Eleotris and Erotelis species. This is the first comprehensive report about the evolutionary relationships in members of the family Eleotridae, including multiloci and multispecies approaches. Therefore, we provided new insights about the phylogenetic position of some taxa absent in previous studies, such as the miniature genus Microphilypnus and a recently described species of Eleotris from South America.

Robust Phylodynamic Analysis of Genetic Sequencing Data from Structured Populations

The multi-type birth–death model with sampling is a phylodynamic model which enables the quantification of past population dynamics in structured populations based on phylogenetic trees. The BEAST 2 package bdmm implements an algorithm for numerically computing the probability density of a phylogenetic tree given the population dynamic parameters under this model. In the initial release of bdmm, analyses were computationally limited to trees consisting of up to approximately 250 genetic samples. We implemented important algorithmic changes to bdmm which dramatically increased the number of genetic samples that could be analyzed and which improved the numerical robustness and efficiency of the calculations. Including more samples led to the improved precision of parameter estimates, particularly for structured models with a high number of inferred parameters. Furthermore, we report on several model extensions to bdmm, inspired by properties common to empirical datasets. We applied this improved algorithm to two partly overlapping datasets of the Influenza A virus HA sequences sampled around the world—one with 500 samples and the other with only 175—for comparison. We report and compare the global migration patterns and seasonal dynamics inferred from each dataset. In this way, we show the information that is gained by analyzing the bigger dataset, which became possible with the presented algorithmic changes to bdmm. In summary, bdmm allows for the robust, faster, and more general phylodynamic inference of larger datasets.

Migration without interbreeding: Evolutionary history of a highly selfing Mediterranean grass inferred from whole genomes

Wild plant populations show extensive genetic subdivision and are far from the ideal of panmixia which permeates population genetic theory. Understanding the spatial and temporal scale of population structure is therefore fundamental for empirical population genetics - and of interest in itself, as it yields insights into the history and biology of a species. In this study we extend the genomic resources for the wild Mediterranean grass Brachypodium distachyon to investigate the scale of population structure and its underlying history at whole-genome resolution. A total of 86 accessions were sampled at local and regional scales in Italy and France, which closes a conspicuous gap in the collection for this model organism. The analysis of 196 accessions, spanning the Mediterranean from Spain to Iraq, suggests that the interplay of high selfing and seed dispersal rates has shaped genetic structure in B. distachyon. At the continental scale, the evolution in B. distachyon is characterized by the independent expansion of three lineages during the Upper Pleistocene. Today, these lineages may occur on the same meadow yet do not interbreed. At the regional scale, dispersal and selfing interact and maintain high genotypic diversity, thus challenging the textbook notion that selfing in finite populations implies reduced diversity. Our study extends the population genomic resources for B. distachyon and suggests that an important use of this wild plant model is to investigate how selfing and dispersal, two processes typically studied separately, interact in colonizing plant species.

The evolution of comparative phylogeography: putting the geography (and more) into comparative population genomics

Comparative population genomics is an ascendant field using genomic comparisons between species to draw inferences about forces regulating genetic variation. Comparative phylogeography, by contrast, focuses on the shared lineage histories of species codistributed geographically and is decidedly organismal in perspective. Comparative phylogeography is approximately 35 years old, and, by some metrics, is showing signs of reduced growth. Here, we contrast the goals and methods of comparative population genomics and comparative phylogeography and argue that comparative phylogeography offers an important perspective on evolutionary history that succeeds in integrating genomics with landscape evolution in ways that complement the suprageographic perspective of comparative population genomics. Focusing primarily on terrestrial vertebrates, we review the history of comparative phylogeography, its milestones and ongoing conceptual innovations, its increasingly global focus, and its status as a bridge between landscape genomics and the process of speciation. We also argue that, as a science with a strong “sense of place,” comparative phylogeography offers abundant “place-based” educational opportunities with its focus on geography and natural history, as well as opportunities for collaboration with local communities and indigenous peoples. Although comparative phylogeography does not yet require whole-genome sequencing for many of its goals, we conclude that it nonetheless plays an important role in grounding our interpretation of genetic variation in the fundamentals of geography and Earth history.

Inferring number of populations and changes in connectivity under the n-island model

Inferring the demographic history of species is one of the greatest challenges in populations genetics. This history is often represented as a history of size changes, ignoring population structure. Alternatively, when structure is assumed, it is defined a priori as a population tree and not inferred. Here we propose a framework based on the IICR (Inverse Instantaneous Coalescence Rate). The IICR can be estimated for a single diploid individual using the PSMC method of Li and Durbin (2011). For an isolated panmictic population, the IICR matches the population size history, and this is how the PSMC outputs are generally interpreted. However, it is increasingly acknowledged that the IICR is a function of the demographic model and sampling scheme with limited connection to population size changes. Our method fits observed IICR curves of diploid individuals with IICR curves obtained under piecewise stationary symmetrical island models. In our models we assume a fixed number of time periods during which gene flow is constant, but gene flow is allowed to change between time periods. We infer the number of islands, their sizes, the periods at which connectivity changes and the corresponding rates of connectivity. Validation with simulated data showed that the method can accurately recover most of the scenario parameters. Our application to a set of five human PSMCs yielded demographic histories that are in agreement with previous studies using similar methods and with recent research suggesting ancient human structure. They are in contrast with the view of human evolution consisting of one ancestral population branching into three large continental and panmictic populations with varying degrees of connectivity and no population structure within each continent.

Geospatial distributions reflect temperatures of linguistic features

Quantifying the speed of linguistic change is challenging because the historical evolution of languages is sparsely documented. Consequently, traditional methods rely on phylogenetic reconstruction. Here, we propose a model-based approach to the problem through the analysis of language change as a stochastic process combining vertical descent, spatial interactions, and mutations in both dimensions. A notion of linguistic temperature emerges naturally from this analysis as a dimensionless measure of the propensity of a linguistic feature to undergo change. We demonstrate how temperatures of linguistic features can be inferred from their present-day geospatial distributions, without recourse to information about their phylogenies. Thus, the evolutionary dynamics of language, operating across thousands of years, leave a measurable geospatial signature. This signature licenses inferences about the historical evolution of languages even in the absence of longitudinal data.

Whole-genome resequencing reveals the pleistocene temporal dynamics of Branchiostoma belcheri and Branchiostoma floridae populations

Global climatic fluctuations governed the ancestral demographic histories of species and contributed to place the current population status into a more extensive ecological and evolutionary context. Genetic variations will leave unambiguous signatures in the patterns of intraspecific genetic variation in extant species since the genome of each individual is an imperfect mosaic of the ancestral genomes. Here, we report the genome sequences of 20 Branchiostoma individuals by whole-genome resequencing strategy. We detected over 140 million genomic variations for each Branchiostoma individual. In particular, we applied the pairwise sequentially Markovian coalescent (PSMC) method to estimate the trajectories of changes in the effective population size (N e) of Branchiostoma population during the Pleistocene. We evaluated the threshold of sequencing depth for proper inference of demographic histories using PSMC was ≥25×. The PSMC results highlight the role of historical global climatic fluctuations in the long-term population dynamics of Branchiostoma. The inferred ancestral N e of the Branchiostoma belcheri populations from Zhanjiang and Xiamen (China) seawaters was different in amplitude before the first (mutation rate = 3 × 10-9) or third glaciation (mutation rate = 9 × 10-9) of the Pleistocene, indicating that the two populations most probably started to evolve in isolation in their respective seas after the first or third glaciation of the Pleistocene. A pronounced population bottleneck coinciding with the last glacial maximum was observed in all Branchiostoma individuals, followed by a population expansion occurred during the late Pleistocene. Species that have experienced long-term declines may be especially vulnerable to recent anthropogenic activities. Recently, the industrial pollution and the exploitation of sea sand have destroyed the harmonious living environment of amphioxus species. In the future, we need to protect the habitat of Branchiostoma and make full use of these detected genetic variations to facilitate the functional study of Branchiostoma for adaptation to local environments.

Phylogeography of the widespread white-eared hummingbird (Hylocharis leucotis): pre-glacial expansion and genetic differentiation of populations separated by the Isthmus of Tehuantepec

The Pleistocene glacial cycles had a strong influence on the demography and genetic structure of many species, particularly on northern-latitude taxa. Here we studied the phylogeography of the white-eared hummingbird (Hylocharis leucotis), a widely distributed species of the highlands of Mexico and Central America. Analysis of mitochondrial DNA (mtDNA) sequences was combined with ecological niche modelling (ENM) to infer the demographic and population differentiation scenarios under present and past conditions. Analyses of 108 samples from 11 geographic locations revealed population structure and genetic differentiation among populations separated by the Isthmus of Tehuantepec (IT) and the Motagua-Polochic-Jocotán (MPJ) fault barriers. ENM predicted a widespread distribution of suitable habitat for H. leucotis since the Last Inter Glacial (LIG), but this habitat noticeably contracted and fragmented at the IT. Models for historical dispersal corridors based on population genetics data and ENM revealed the existence of corridors among populations west of the IT; however, the connectivity of populations across the IT has changed little since the LIG. The shallow geographic structure on either side of the isthmus and a star-like haplotype network, combined with the long-term persistence of populations across time based on genetic data and potential dispersal routes, support a scenario of divergence with migration and subsequent isolation and differentiation in Chiapas and south of the MPJ fault. Our findings corroborate the profound effects of Pleistocene climatic fluctuations on the evolutionary history of montane taxa but challenge the generality of expanded suitable habitat (pine-oak forests) during glacial cycles.

Plant population dynamics on oceanic islands during the Late Quaternary climate changes: genetic evidence from a tree species (Coffea mauritiana) in Reunion Island

Past climatic fluctuations have played a major role in shaping the current plant biodiversity. Although harbouring an exceptional biota, oceanic islands have received little attention in studies on species demographic history and past vegetation patterns. We investigated the impact of past climatic changes on the effective population size of a tree (Coffea mauritiana) that is endemic to Reunion Island, located in the south-western Indian Ocean (SWIO). Demographic changes were inferred using summary statistics calculated from genomic data. Using ecological niche modelling and the current distribution of genetic diversity, the paleodistribution of the species was also assessed. A reduction in the effective population size of C. mauritiana during the last glaciation maximum was inferred. The distribution of the species was reduced on the western side of the island, due to low rainfall. It appeared that a major reduction in rainfall and a slight temperature decrease prevailed in the SWIO. Our findings indicated that analyses on the current patterns of intraspecific genetic variations can efficiently contribute to past climatic changes characterisation in remote islands. Identifying area with higher resilience in oceanic islands could provide guidance in forest management and conservation faced to the global climate change.

The effect of undetected recombination on genealogy sampling and inference under an isolation-with-migration model

Many methods for fitting demographic models to data sets of aligned sequences rely upon an assumption that the data have a branching coalescent history without recombination within regions or loci. To mitigate the effects of the failure of this assumption, a common approach is to filter data and sample regions that pass the four-gamete criterion for recombination, an approach that allows data to run, but that is expected to detect only a minority of recombination events. A series of empirical tests of this approach were conducted using computer simulations with and without recombination for a variety of isolation-with-migration (IM) model for two and three populations. Only the IMa3 program was used, but the general results should apply to related genealogy-sampling-based methods for IM models or subsets of IM models. It was found that the details of sampling intervals that pass a four-gamete filter have a moderate effect, and that schemes that use the longest intervals, or that use overlapping intervals, gave poorer results. A simple approach of using a random nonoverlapping interval returned the smallest difference between results with and without recombination, with the mean difference between parameter estimates usually less than 20% of the true value (usually much less). However, the posterior probability distributions for migration rates were flatter with recombination, suggesting that filtering based on the four-gamete criterion, while necessary for methods like these, leads to reduced resolution on migration. A distinct, alternative approach, of using a finite sites mutation model and not filtering the data, performed quite poorly.

Old divergence and restricted gene flow between torrent duck (Merganetta armata) subspecies in the Central and Southern Andes

AIM

To investigate the structure and rate of gene flow among populations of habitat-specialized species to understand the ecological and evolutionary processes underpinning their population dynamics and historical demography, including speciation and extinction.

LOCATION

Peruvian and Argentine Andes.

TAXON

Two subspecies of torrent duck (Merganetta armata).

METHODS

We sampled 156 individuals in Peru (M. a. leucogenis; Chillón River, n = 57 and Pachachaca River, n = 49) and Argentina (M. a. armata; Arroyo Grande River, n = 33 and Malargüe River, n = 17), and sequenced the mitochondrial DNA (mtDNA) control region to conduct coarse and fine-scale demographic analyses of population structure. Additionally, to test for differences between subspecies, and across genetic markers with distinct inheritance patterns, a subset of individuals (Peru, n = 10 and Argentina, n = 9) was subjected to partial genome resequencing, obtaining 4,027 autosomal and 189 Z-linked double-digest restriction-associated DNA sequences.

RESULTS

Haplotype and nucleotide diversities were higher in Peru than Argentina across all markers. Peruvian and Argentine subspecies showed concordant species-level differences (ΦST mtDNA = 0.82; ΦST autosomal = 0.30; ΦST Z chromosome = 0.45), including no shared mtDNA haplotypes. Demographic parameters estimated for mtDNA using IM and IMa2 analyses, and for autosomal markers using ∂a∂i (isolation-with-migration model), supported an old divergence (mtDNA = 600,000 years before present (ybp), 95% HPD range = 1.2 Mya to 200,000 ybp; and autosomal ∂a∂i = 782,490 ybp), between the two subspecies, characteristic of deeply diverged lineages. The populations were well-differentiated in Argentina but moderately differentiated in Peru, with low unidirectional gene flow in each country.

MAIN CONCLUSIONS

We suggest that the South American Arid Diagonal was preexisting and remains a current phylogeographic barrier between the ranges of the two torrent duck subspecies, and the adult territoriality and breeding site fidelity to the rivers define their population structure.

Amplicon sequencing of 42 nuclear loci supports directional gene flow between South Pacific populations of a hydrothermal vent limpet

In the past few decades, population genetics and phylogeographic studies have improved our knowledge of connectivity and population demography in marine environments. Studies of deep-sea hydrothermal vent populations have identified barriers to gene flow, hybrid zones, and demographic events, such as historical population expansions and contractions. These deep-sea studies, however, used few loci, which limit the amount of information they provided for coalescent analysis and thus our ability to confidently test complex population dynamics scenarios. In this study, we investigated population structure, demographic history, and gene flow directionality among four Western Pacific hydrothermal vent populations of the vent limpet Lepetodrilus aff. schrolli. These vent sites are located in the Manus and Lau back-arc basins, currently of great interest for deep-sea mineral extraction. A total of 42 loci were sequenced from each individual using high-throughput amplicon sequencing. Amplicon sequences were analyzed using both genetic variant clustering methods and evolutionary coalescent approaches. Like most previously investigated vent species in the South Pacific, L. aff. schrolli showed no genetic structure within basins but significant differentiation between basins. We inferred significant directional gene flow from Manus Basin to Lau Basin, with low to no gene flow in the opposite direction. This study is one of the very few marine population studies using >10 loci for coalescent analysis and serves as a guide for future marine population studies.

Phenotypic and genome-wide association with the local environment of Arabidopsis

The environment imposes critical selective forces on all living organisms, and the sessile nature of plants makes them particularly useful for investigating the relationship between genetic variation and environmental adaptation. In the model plant Arabidopsis thaliana, extensive information on phenotypic and genotypic variation is available, but comparable information on environmental variation within the native range of the species is lacking. Here, we compile 204 geoclimatic variables to describe the local environments of Arabidopsis accessions with known collection sites encompassing a wide geo-environmental range, and fully sequenced genomes from the 1001 Genomes Project. We identify candidate adaptive genetic variation associated with these environmental variables, and validate this approach through comparison with previous experimental studies, and by targeted confirmation of a role of the heterotrimeric G-protein γ subunit, AGG3, in cold tolerance, as newly predicted from our environmental genome wide association study (GWAS). To facilitate identification of adaptive variation, we created Arabidopsis CLIMtools : interactive web-based databases of the environment × genome associations and correlations between the local environments and 131 phenotypes compiled from previous experimental GWASs. Our study presents an extensive analysis of the local environments, landscape genomics and phenotypic variation of Arabidopsis, and illustrates how 'in silico GWAS' approaches can inform and complement experimental phenomics studies.

Testing Hypotheses of Diversification in Panamanian Frogs and Freshwater Fishes Using Hierarchical Approximate Bayesian Computation with Model Averaging

Most Neotropical frog and freshwater fish species sampled to date show phylogeographic breaks along the Pacific coast of the Isthmus of Panama, with lineages in Costa Rica and western Panama isolated from central Panama. We examine temporal patterns of diversification of taxa across this ‘western Panama isthmus’ (WPI) break to test hypotheses about the origin of species geographical distributions and genetic structuring in this region. We tested for synchronous diversification of four codistributed frog taxon-pairs and three fish taxon-pairs sharing the WPI break using hierarchical approximate Bayesian computation with model averaging based on mitochondrial DNA sequences. We also estimated lineage divergence times using full-Bayesian models. Several of our results supported synchronous divergences within the frog and freshwater fish assemblages; however, Bayes factor support was equivocal for or against synchronous or asynchronous diversification. Nevertheless, we infer that frog populations were likely isolated by one or multiple Pliocene–Pleistocene events more recently than predicted by previous models, while fish genetic diversity was structured by Pleistocene events. By integrating our results with external information from geology and elevational sea level modeling, we discuss the implications of our findings for understanding the biogeographical scenario of the diversification of Panamanian frogs and fishes. Consistent with the ‘Bermingham/Martin model’ (Molecular Ecology 1998, 7, 499–517), we conclude that the regional fish assemblage was fractured by processes shaping isthmian landscapes during the Pleistocene glaciations, including drainage basin isolation during lowered sea levels.

Unveiling patterns of genetic variation in parasite–host associations: an example with pinworms and Neotropical primates

Patterns of genetic variation among populations can reveal the evolutionary history of species. Pinworm parasites are highly host specific and form strong co-evolutionary associations with their primate hosts. Here, we describe the genetic variation observed in four Trypanoxyuris species infecting different howler and spider monkey subspecies in Central America to determine if historical dispersal processes and speciation in the host could explain the genetic patterns observed in the parasites. Mitochondrial (cox1) and ribosomal (28S) DNA were analysed to assess genetic divergence and phylogenetic history of these parasites. Sequences of the 28S gene were identical within pinworms species regardless of host subspecies. However, phylogenetic analyses, haplotype relationships and genetic divergence with cox1 showed differentiation between pinworm populations according to host subspecies in three of the four Trypanoxyuris species analysed. Haplotype separation between host subspecies was not observed in Trypanoxyuris minutus, nor in Trypanoxyuris atelis from Ateles geoffoyi vellerosus and Ateles geoffoyi yucatanensis. Levels of genetic diversity and divergence in these parasites relate with such estimates reported for their hosts. This study shows how genetic patterns uncovered in parasitic organisms can reflect the host phylogenetic and biogeographic histories.

The IICR and the non-stationary structured coalescent: towards demographic inference with arbitrary changes in population structure

In the last years, a wide range of methods allowing to reconstruct past population size changes from genome-wide data have been developed. At the same time, there has been an increasing recognition that population structure can generate genetic data similar to those produced under models of population size change. Recently, Mazet et al. (Heredity 116:362-371, 2016) showed that, for any model of population structure, it is always possible to find a panmictic model with a particular function of population size changes, having exactly the same distribution of T₂ (the coalescence time for a sample of size two) as that of the structured model. They called this function IICR (Inverse Instantaneous Coalescence Rate) and showed that it does not necessarily correspond to population size changes under non-panmictic models. Besides, most of the methods used to analyse data under models of population structure tend to arbitrarily fix that structure and to minimise or neglect population size changes. Here, we extend the seminal work of Herbots (PhD thesis, University of London, 1994) on the structured coalescent and propose a new framework, the Non-Stationary Structured Coalescent (NSSC) that incorporates demographic events (changes in gene flow and/or deme sizes) to models of nearly any complexity. We show how to compute the IICR under a wide family of stationary and non-stationary models. As an example we address the question of human and Neanderthal evolution and discuss how the NSSC framework allows to interpret genomic data under this new perspective.

Disentangling genetic structure for genetic monitoring of complex populations

Genetic monitoring estimates temporal changes in population parameters from molecular marker information. Most populations are complex in structure and change through time by expanding or contracting their geographic range, becoming fragmented or coalescing, or increasing or decreasing density. Traditional approaches to genetic monitoring rely on quantifying temporal shifts of specific population metrics-heterozygosity, numbers of alleles, effective population size-or measures of geographic differentiation such as FST. However, the accuracy and precision of the results can be heavily influenced by the type of genetic marker used and how closely they adhere to analytical assumptions. Care must be taken to ensure that inferences reflect actual population processes rather than changing molecular techniques or incorrect assumptions of an underlying model of population structure. In many species of conservation concern, true population structure is unknown, or structure might shift over time. In these cases, metrics based on inappropriate assumptions of population structure may not provide quality information regarding the monitored population. Thus, we need an inference model that decouples the complex elements that define population structure from estimation of population parameters of interest and reveals, rather than assumes, fine details of population structure. Encompassing a broad range of possible population structures would enable comparable inferences across biological systems, even in the face of range expansion or contraction, fragmentation, or changes in density. Currently, the best candidate is the spatial Λ-Fleming-Viot (SLFV) model, a spatially explicit individually based coalescent model that allows independent inference of two of the most important elements of population structure: local population density and local dispersal. We support increased use of the SLFV model for genetic monitoring by highlighting its benefits over traditional approaches. We also discuss necessary future directions for model development to support large genomic datasets informing real-world management and conservation issues.

What happened in the South American Gran Chaco? Diversification of the endemic frog genus Lepidobatrachus Budgett, 1899 (Anura: Ceratophryidae)

The Chaco is one the most neglected and least studied regions of the world. This highly-seasonal semiarid biome is an extensive continuous plain without any geographic barrier, and in spite of its high species diversity, the events and processes responsible have never been assessed. Miocene marine introgressions and Pleistocene glaciations have been mentioned as putative drivers of diversification for some groups of vertebrates in adjacent biomes of southern South America. Here we used multilocus data (one mitochondrial and six nuclear loci) from the three species of the endemic frog genus Lepidobatrachus (Lepidobatrachus asper, Lepidobatrachus laevis, and Lepidobatrachus llanensis) to determine if any of the historical events suggested as drivers of vertebrate diversification in southern South America are related to the diversification of the genus and if the Chaco is indeed a biome without barriers. Using fossil calibration in a coalescent framework we estimated that the genus diversified in the second half of the Miocene, coinciding with marine introgressions. Genetic patterns and historical demography suggest an important role of old archs and cratons as refuges during floods. In one species of the genus, L. llanensis, genetic structure reveals some breaks along the landscape, the main one of which corresponds to an area of the central Chaco that may act as a climatic barrier. Additionally, we found differential effects of the main Chacoan rivers on species of Lepidobatrachus that could be related to the time of persistence of populations in the areas influenced by these rivers.

Exact Calculation of the Joint Allele Frequency Spectrum for Isolation with Migration Models

Population genomic datasets collected over the past decade have spurred interest in developing methods that can utilize massive numbers of loci for inference of demographic and selective histories of populations. The allele frequency spectrum (AFS) provides a convenient statistic for such analysis, and, accordingly, much attention has been paid to predicting theoretical expectations of the AFS under a number of different models. However, to date, exact solutions for the joint AFS of two or more populations under models of migration and divergence have not been found. Here, we present a novel Markov chain representation of the coalescent on the state space of the joint AFS that allows for rapid, exact calculation of the joint AFS under isolation with migration (IM) models. In turn, we show how our Markov chain method, in the context of composite likelihood estimation, can be used for accurate inference of parameters of the IM model using SNP data. Lastly, we apply our method to recent whole genome datasets from African Drosophila melanogaster.

A novel genus and cryptic species harboured within the monotypic freshwater crayfish genus Tenuibranchiurus Riek, 1951 (Decapoda: Parastacidae)

Identifying species groups is an important yet difficult task, with there being no single accepted definition as to what constitutes a species, nor a set of criteria by which they should be delineated. Employing the General Lineage Concept somewhat circumvents these issues, as this concept allows multiple concordant lines of evidence to be used as support for species delimitation, where a species is defined as any independently evolving lineage. Genetically diverse groups have previously been identified within the monotypic parastacid genus Tenuibranchiurus Riek, 1951, but no further investigation of this diversity has previously been undertaken. Analysis of two mitochondrial DNA gene regions has previously identified two highly divergent groups within this taxon, representing populations from Queensland (Qld) and New South Wales (NSW), respectively. Additional testing within this study of both mitochondrial and nuclear DNA through species discovery analyses identified genetically diverse groups within these regions, which were further supported by lineage validation methods. The degree of genetic differentiation between Qld and NSW populations supports the recognition of two genera; with Qld retaining the original genus name Tenuibranchiurus, and NSW designated as Gen. nov. until a formal description is completed. Concordance between the species discovery and lineage validation methods supports the presence of six species within Tenuibranchiurus and two within Gen. nov. The recognition of additional species removes the monotypy of the genus, and the methods used can improve species identification within groups of organisms with taxonomic problems and cryptic diversity.

Effects of the population pedigree on genetic signatures of historical demographic events

Genetic variation among loci in the genomes of diploid biparental organisms is the result of mutation and genetic transmission through the genealogy, or population pedigree, of the species. We explore the consequences of this for patterns of variation at unlinked loci for two kinds of demographic events: the occurrence of a very large family or a strong selective sweep that occurred in the recent past. The results indicate that only rather extreme versions of such events can be expected to structure population pedigrees in such a way that unlinked loci will show deviations from the standard predictions of population genetics, which average over population pedigrees. The results also suggest that large samples of individuals and loci increase the chance of picking up signatures of these events, and that very large families may have a unique signature in terms of sample distributions of mutant alleles.

Toward a paradigm shift in comparative phylogeography driven by trait-based hypotheses

For three decades, comparative phylogeography has conceptually and methodologically relied on the concordance criterion for providing insights into the historical/biogeographic processes driving population genetic structure and divergence. Here we discuss how this emphasis, and the corresponding lack of methods for extracting information about biotic/intrinsic contributions to patterns of genetic variation, may bias our general understanding of the factors driving genetic structure. Specifically, this emphasis has promoted a tendency to attribute discordant phylogeographic patterns to the idiosyncracies of history, as well as an adherence to generic null expectations of concordance with reduced predictive power. We advocate that it is time for a paradigm shift in comparative phylogeography, especially given the limited utility of the concordance criterion as genomic data provide ever-increasing levels of resolution. Instead of adhering to the concordance-discordance dichotomy, comparative phylogeography needs to emphasize the contribution of taxon-specific traits that will determine whether concordance is a meaningful criterion for evaluating hypotheses or may predict discordant phylogeographic structure. Through reference to some case studies we illustrate how refined hypotheses based on taxon-specific traits can provide improved predictive frameworks to forecast species responses to climatic change or biogeographic barriers while gaining unique insights about the taxa themselves and their interactions with their environment. We outline a potential avenue toward a synthetic comparative phylogeographic paradigm that includes addressing some important conceptual and methodological challenges related to study design and application of model-based approaches for evaluating support of trait-based hypotheses under the proposed paradigm.

Genetic demographic networks: Mathematical model and applications

Recent improvement in the quality of genetic data obtained from extinct human populations and their ancestors encourages searching for answers to basic questions regarding human population history. The most common and successful are model-based approaches, in which genetic data are compared to the data obtained from the assumed demography model. Using such approach, it is possible to either validate or adjust assumed demography. Model fit to data can be obtained based on reverse-time coalescent simulations or forward-time simulations. In this paper we introduce a computational method based on mathematical equation that allows obtaining joint distributions of pairs of individuals under a specified demography model, each of them characterized by a genetic variant at a chosen locus. The two individuals are randomly sampled from either the same or two different populations. The model assumes three types of demographic events (split, merge and migration). Populations evolve according to the time-continuous Moran model with drift and Markov-process mutation. This latter process is described by the Lyapunov-type equation introduced by O'Brien and generalized in our previous works. Application of this equation constitutes an original contribution. In the result section of the paper we present sample applications of our model to both simulated and literature-based demographies. Among other we include a study of the Slavs-Balts-Finns genetic relationship, in which we model split and migrations between the Balts and Slavs. We also include another example that involves the migration rates between farmers and hunters-gatherers, based on modern and ancient DNA samples. This latter process was previously studied using coalescent simulations. Our results are in general agreement with the previous method, which provides validation of our approach. Although our model is not an alternative to simulation methods in the practical sense, it provides an algorithm to compute pairwise distributions of alleles, in the case of haploid non-recombining loci such as mitochondrial and Y-chromosome loci in humans.

The mutation-drift balance in spatially structured populations

In finite populations the action of neutral mutations is balanced by genetic drift, leading to a stationary distribution of alleles that displays a transition between two different behaviors. For small mutation rates most individuals will carry the same allele at equilibrium, whereas for high mutation rates of the alleles will be randomly distributed with frequencies close to one half for a biallelic gene. For well-mixed haploid populations the mutation threshold is μc=1/2N, where N is the population size. In this paper we study how spatial structure affects this mutation threshold. Specifically, we study the stationary allele distribution for populations placed on regular networks where connected nodes represent potential mating partners. We show that the mutation threshold is sensitive to spatial structure only if the number of potential mates is very small. In this limit, the mutation threshold decreases substantially, increasing the diversity of the population at considerably low mutation rates. Defining kc as the degree of the network for which the mutation threshold drops to half of its value in well-mixed populations we show that kc grows slowly as a function of the population size, following a power law. Our calculations and simulations are based on the Moran model and on a mapping between the Moran model with mutations and the voter model with opinion makers.

PSMC analysis of effective population sizes in molecular ecology and its application to black-and-white Ficedula flycatchers

Climatic fluctuations during the Quaternary period governed the demography of species and contributed to population differentiation and ultimately speciation. Studies of these past processes have previously been hindered by a lack of means and genetic data to model changes in effective population size (Ne ) through time. However, based on diploid genome sequences of high quality, the recently developed pairwise sequentially Markovian coalescent (PSMC) can estimate trajectories of changes in Ne over considerable time periods. We applied this approach to resequencing data from nearly 200 genomes of four species and several populations of the Ficedula species complex of black-and-white flycatchers. Ne curves of Atlas, collared, pied and semicollared flycatcher converged 1-2 million years ago (Ma) at an Ne of ≈ 200 000, likely reflecting the time when all four species last shared a common ancestor. Subsequent separate Ne trajectories are consistent with lineage splitting and speciation. All species showed evidence of population growth up until 100-200 thousand years ago (kya), followed by decline and then start of a new phase of population expansion. However, timing and amplitude of changes in Ne differed among species, and for pied flycatcher, the temporal dynamics of Ne differed between Spanish birds and central/northern European populations. This cautions against extrapolation of demographic inference between lineages and calls for adequate sampling to provide representative pictures of the coalescence process in different species or populations. We also empirically evaluate criteria for proper inference of demographic histories using PSMC and arrive at recommendations of using sequencing data with a mean genome coverage of ≥18X, a per-site filter of ≥10 reads and no more than 25% of missing data.

New Routes to Phylogeography: A Bayesian Structured Coalescent Approximation

Phylogeographic methods aim to infer migration trends and the history of sampled lineages from genetic data. Applications of phylogeography are broad, and in the context of pathogens include the reconstruction of transmission histories and the origin and emergence of outbreaks. Phylogeographic inference based on bottom-up population genetics models is computationally expensive, and as a result faster alternatives based on the evolution of discrete traits have become popular. In this paper, we show that inference of migration rates and root locations based on discrete trait models is extremely unreliable and sensitive to biased sampling. To address this problem, we introduce BASTA (BAyesian STructured coalescent Approximation), a new approach implemented in BEAST2 that combines the accuracy of methods based on the structured coalescent with the computational efficiency required to handle more than just few populations. We illustrate the potentially severe implications of poor model choice for phylogeographic analyses by investigating the zoonotic transmission of Ebola virus. Whereas the structured coalescent analysis correctly infers that successive human Ebola outbreaks have been seeded by a large unsampled non-human reservoir population, the discrete trait analysis implausibly concludes that undetected human-to-human transmission has allowed the virus to persist over the past four decades. As genomics takes on an increasingly prominent role informing the control and prevention of infectious diseases, it will be vital that phylogeographic inference provides robust insights into transmission history.

When studying infectious diseases it is often important to understand how germs spread from location-to-location, person-to-person, or even one part of the body to another. Using phylogeographic methods, it is possible to recover the history of spread of pathogens (or other organisms) by studying their genetic material. Here we reveal that some popular, fast phylogeographic methods are inaccurate, and we introduce a new more reliable method to address the problem. By comparing different phylogeographic methods based on principled population models and fast alternatives, we found that different approaches can give diametrically opposed results, and we offer concrete examples in the context of the ongoing Ebola outbreak in West Africa and the world-wide outbreaks of Avian Influenza Virus and Tomato Yellow Leaf Curl Virus. We found that the most popular phylogeographic method often produces completely inaccurate conclusions. One of the reasons for its popularity has been its computational speed, which has allowed users to analyse large genetic datasets with complex models. More accurate approaches have until now been considerably slower, and therefore we propose a new method called BASTA that achieves good accuracy in a reasonable time. We are relying more and more on genetic sequencing to learn about the origin and spread of infections, and as this role continues to grow, it will be essential to use accurate phylogeographic methods when designing policies to prevent or curb the spread of disease.

Association Analysis of Simple Sequence Repeat (SSR) Markers with Agronomic Traits in Tall Fescue (Festuca arundinacea Schreb.)

Tall fescue is widely used in temperate regions throughout the world as a dominant forage grass as well as a turfgrass, in pastoral and turf industry. However, the utilization of tall fescue was limited because of its leaf roughness, poor regeneration ability and poor stress resistance. New cultivars were desirable in modern pastoral industries exceed the potential of existing cultivars. Therefore, well understanding the agronomic traits and describing germplasms would help to overcome these constraints, and morphological evaluation of tall fescue germplasm is the key component in selecting rational parents for hybridization breeding. However, describing the morphological traits of tall fescue germplasm is costly and time-consuming. Fortunately, biotechnology approaches can supplement conventional breeding efforts for tall fescue improvement. Association mapping, as a powerful approach to identify association between agronomic traits and molecular markers has been widely used for enhancing the utilization, conservation and management of the tall fescue germplasms. Therefore, in the present research, 115 tall fescue accessions from different origins (25 accessions are cultivars; 31 accessions from America; 32 accessions from European; 7 accessions from Africa; 20 accessions from Asia), were evaluated for agronomic traits and genetic diversity with 90 simple sequence repeat (SSR) markers. The panel displayed significant variation in spike count per plant (SCP) and spike weight (SW). However, BCS performed the lowest CV among all the observed agronomic traits. Three subpopulations were identified within the collections but no obvious relative kinship (K) was found. The GLM model was used to describe the association between SSR and agronomic traits. Fifty-one SSR markers associated with agronomic traits were observed. Twelve single-associated markers were associated with PH; six single-associated markers were associated with BCS; eight single-associated markers were associated with SW; five single-associated markers were associated with SC; seven single-associated markers were associated with SCP; three single-associated markers were associated with SL. Especially, we observed that the genetic variation of SW was explained 11.6 % by M37 marker. It is interesting to observe that nine markers (M1, M2, M35, M54 marker was associated with both BCS and SC; M3, M4 markers were associated with BCS, SW, and SC; M19 marker was associated with both pH and PD, M40 marker was associated with both SCP and SW; and M193 marker was associated with both PH and SL) were associated with more than two agronomic traits. Notably, Branch count per spike (BCS) was explained by four markers (M1, M2, M3, and M4) exceeding 10 %. These identified marker alleles associated with agronomic traits could provide important information and markers for molecular-assisted breeding that facilitate the breeding process in tall fescue.

Persistent history of the bird-dispersed arctic-alpine plant Vaccinium vitis-idaea L. (Ericaceae) in Japan

Arctic-alpine plants have expanded and contracted their ranges in response to the Pleistocene climate oscillations. Today, many arctic-alpine plants have vast distributions in the circumarctic region as well as marginal, isolated occurrences in high mountains at lower latitudes. These marginal populations may represent relict, long-standing populations that have persisted for several cycles of cold and warm climate during the Pleistocene, or recent occurrences that either result from southward step-wise migration during the last glacial period or from recent long-distance dispersal. In light of these hypotheses, we investigated the biogeographic history of the marginal Japanese populations of the widespread arctic-alpine plant Vaccinium vitis-idaea (Ericaceae), which is bird-dispersed, potentially over long distances. We sequenced three nuclear loci and one plastid DNA region in 130 individuals from 65 localities covering its entire geographic range, with a focus on its marginal populations in Japan. We found a homogenous genetic pattern across its enormous range based on the loci analysed, in contrast to the geographically structured variation found in a previous study of amplified fragment length polymorphisms in this species. However, we found several unique haplotypes in the Japanese populations, excluding the possibility that these marginal populations result from recent southward migration. Thus, even though V. vitis-idaea is efficiently dispersed via berries, our study suggests that its isolated populations in Japan have persisted during several cycles of cold and warm climate during the Pleistocene.

The effects of purifying selection on patterns of genetic differentiation between Drosophila melanogaster populations

Using the data provided by the Drosophila Population Genomics Project, we investigate factors that affect the genetic differentiation between Rwandan and French populations of D. melanogaster. By examining within-population polymorphisms, we show that sites in long introns (especially those >2000 bp) have significantly lower π (nucleotide diversity) and more low-frequency variants (as measured by Tajima's D, minor allele frequencies, and prevalence of variants that are private to one of the two populations) than short introns, suggesting a positive relationship between intron length and selective constraint. A similar analysis of protein-coding polymorphisms shows that 0-fold (degenerate) sites in more conserved genes are under stronger purifying selection than those in less conserved genes. There is limited evidence that selection on codon bias has an effect on differentiation (as measured by F_ST) at 4-fold (degenerate) sites, and 4-fold sites and sites in 8–30 bp of short introns ⩽65 bp have comparable F_ST values. Consistent with the expected effect of purifying selection, sites in long introns and 0-fold sites in conserved genes are less differentiated than those in short introns and less conserved genes, respectively. Genes in non-crossover regions (for example, the fourth chromosome) have very high F_ST values at both 0-fold and 4-fold degenerate sites, which is probably because of the large reduction in within-population diversity caused by tight linkage between many selected sites. Our analyses also reveal subtle statistical properties of F_ST, which arise when information from multiple single nucleotide polymorphisms is combined and can lead to the masking of important signals of selection.

Phylogeography and biogeography of the lower Central American Neotropics: diversification between two continents and between two seas

Lower Central America (LCA) provides a geologically complex and dynamic, richly biodiverse model for studying the recent assembly and diversification of a Neotropical biota. Here, we review the growing literature of LCA phylogeography studies and their contribution to understanding the origins, assembly, and diversification of the LCA biota against the backdrop of regional geologic and climatic history, and previous biogeographical inquiry. Studies to date reveal that phylogeographical signal within taxa of differing distributions reflects a diversity of patterns and processes rivalling the complexities of LCA landscapes themselves. Even so, phylogeography is providing novel insights into regional diversification (e.g. cryptic lineage divergences), and general evolutionary patterns are emerging. Congruent multi-taxon phylogeographic breaks are found across the Nicaraguan depression, Chorotega volcanic front, western and central Panama, and the Darién isthmus, indicating that a potentially shared history of responses to regional-scale (e.g. geological) processes has shaped the genetic diversity of LCA communities. By contrast, other species show unique demographic histories in response to overriding historical events, including no phylogeographic structure at all. These low-structure or incongruent patterns provide some evidence for a role of local, ecological factors (e.g. long-distance dispersal and gene flow in plants and bats) in shaping LCA communities. Temporally, comparative phylogeographical structuring reflects Pliocene-Pleistocene dispersal and vicariance events consistent with the timeline of emergence of the LCA isthmus and its major physiographic features, e.g. cordilleras. We emphasise the need to improve biogeographic inferences in LCA through in-depth comparative phylogeography projects capitalising on the latest statistical phylogeographical methods. While meeting the challenges of reconstructing the biogeographical history of this complex region, phylogeographers should also take up the critical service to society of applying their work to the conservation of its fascinating biodiversity.

Coalescent analyses show isolation without migration in two closely related tropical orioles: the case of Icterus graduacauda and Icterus chrysater

The Isthmus of Tehuantepec has played an important role in shaping the avian diversity of Mexico, as well as the rest of the Western Hemisphere. It has been both a barrier and a land connector between North and South America for many groups of birds. Furthermore, climatic change over the Pleistocene has resulted in ecological fluctuations that led to periods of connection and isolation of the highlands in this area. Here we studied the divergence of two species of orioles whose distribution in the highlands is separated by the lowlands of the Isthmus of Tehuantepec: Icterus graduacauda (west of the Isthmus) and Icterus chrysater (east of the Isthmus). We sequenced multiple loci (one mitochondrial gene and six nuclear introns) and performed coalescent analyses (Isolation with Migration) to test whether their divergence resulted from prior occupancy of the ancestral area followed by a vicariant event or recent dispersal from one side or the other of this Isthmus. Results strongly indicate a vicariant event roughly 300,000 years ago in the Pleistocene followed by little or no gene flow. Both mitochondrial and nuclear genes show that the Isthmus of Tehuantepec is a strong barrier to gene flow. Thus, these two species appear to not exchange genes despite their recent divergence and the close geographic proximity of their ranges.