Host‐associated genomic differentiation in congeneric butterflies: now you see it,now you do not

Ecotypic variation among populations may become associated with widespread genomic differentiation, but theory predicts that this should happen only under particular conditions of gene flow, selection and population size. In closely related species, we might expect the strength of host‐associated genomic differentiation (HAD) to be correlated with the degree of phenotypic differentiation in host‐adaptive traits. Using microsatellite and Amplified Fragment Length Polymorphism (AFLP) markers, and controlling for isolation by distance between populations, we sought HAD in two congeneric species of butterflies with different degrees of host plant specialization. Prior work on Euphydryas editha had shown strong interpopulation differentiation in host‐adapted traits, resulting in incipient reproductive isolation among host‐associated ecotypes. We show here that Euphydryas aurinia had much weaker host‐associated phenotypic differentiation. Contrary to our expectations, we detected HAD in Euphydryas aurinia, but not in E. editha. Even within an E. aurinia population that fed on both hosts, we found weak but significant sympatric HAD that persisted in samples taken 9 years apart. The finding of significantly stronger HAD in the system with less phenotypic differentiation may seem paradoxical. Our findings can be explained by multiple factors, ranging from differences in dispersal or effective population size, to spatial variation in genomic or phenotypic traits and to structure induced by past histories of host‐adapted populations. Other infrequently measured factors, such as differences in recombination rates, may also play a role. Our result adds to recent work as a further caution against assumptions of simple relationships between genomic and adaptive phenotypic differentiation.     

High levels of interspecific gene flow in an endemic cichlid fish adaptive radiation from an extreme lake environment

Studying recent adaptive radiations in isolated insular systems avoids complicating causal events and thus may offer clearer insight into mechanisms generating biological diversity. Here, we investigate evolutionary relationships and genomic differentiation within the recent radiation of Alcolapia cichlid fish that exhibit extensive phenotypic diversification, and which are confined to the extreme soda lakes Magadi and Natron in East Africa. We generated an extensive RAD data set of 96 individuals from multiple sampling sites and found evidence for genetic admixture between species within Lake Natron, with the highest levels of admixture between sympatric populations of the most recently diverged species. Despite considerable environmental separation, populations within Lake Natron do not exhibit isolation by distance, indicating panmixia within the lake, although individuals within lineages clustered by population in phylogenomic analysis. Our results indicate exceptionally low genetic differentiation across the radiation despite considerable phenotypic trophic variation, supporting previous findings from smaller data sets; however, with the increased power of densely sampled SNPs, we identify genomic peaks of differentiation (F_ST outliers) between Alcolapia species. While evidence of ongoing gene flow and interspecies hybridization in certain populations suggests that Alcolapia species are incompletely reproductively isolated, the identification of outlier SNPs under diversifying selection indicates the radiation is undergoing adaptive divergence.     

Diversifying selection drives parallel evolution of gill raker number and body size along the speciation continuum of European whitefish

Adaptive radiation is the evolution of ecological and phenotypical diversity. It arises via ecological opportunity that promotes the exploration of underutilized or novel niches mediating specialization and reproductive isolation. The assumed precondition for rapid local adaptation is diversifying natural selection, but random genetic drift could also be a major driver of this process. We used 27 populations of European whitefish (Coregonus lavaretus) from nine lakes distributed in three neighboring subarctic watercourses in northern Fennoscandia as a model to test the importance of random drift versus diversifying natural selection for parallel evolution of adaptive phenotypic traits. We contrasted variation for two key adaptive phenotypic traits correlated with resource utilization of polymorphic fish; the number of gill rakers and the total length of fish, with the posterior distribution of neutral genetic differentiation from 13 microsatellite loci, to test whether the observed phenotypic divergence could be achieved by random genetic drift alone. Our results show that both traits have been under diversifying selection and that the evolution of these morphs has been driven by isolation through habitat adaptations. We conclude that diversifying selection acting on gill raker number and body size has played a significant role in the ongoing adaptive radiation of European whitefish morphs in this region.     

Patterns of genetic and morphometric differentiation in Melitaea (Mellicta) athalia (Lepidoptera: Nymphalidae)

Patterns of genetic and morphometric differentiation were surveyed in Melitaea (Mellicta) athalia populations of the Carpathian Basin. This species has a wide distribution and exists under a wide variety of ecological conditions. It has two ecotypes in Hungary: with either one or two broods per year. It is of particular interest to reveal the main factors driving differentiation patterns in this species. Samples in our study were obtained from five Hungarian and one Transylvanian (Romanian) regions. Enzyme polymorphism, wing characters and male external genitalia were analysed using traditional morphometric methods. Statistical methods were optimized to compare morphological and genetic data. The results of genetic surveys revealed a clear regional pattern of differentiation in M. athalia. Moreover, the results of principal component analysis, Bayesian clustering and the dendrogram all suggested that the regions can be classified into two groups corresponding to the East or West zones of the Carpathian Basin. In contrast, differentiation between the two ecotypes was less expressed in the genetic variation of M. athalia. Results of the analyses conducted on phenotypic variation also suggested a regional pattern for both sets of morphometric characters (wings and external genitalia). At the same time, neither East–West regional division nor ecotype differentiation was detected in the morphometric studies. In sum, our analyses confirmed that both genetic and phenotypic variations of M. athalia exhibit a regional pattern rather than the differentiation between the two ecotypes. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, ?? , ??–??.     

Population genomic evidence for plant glacial survival in Scandinavia

Quaternary glaciations have played a major role in shaping the genetic diversity and distribution of plant species. Strong palaeoecological and genetic evidence supports a postglacial recolonization of most plant species to northern Europe from southern, eastern and even western glacial refugia. Although highly controversial, the existence of small in situ glacial refugia in northern Europe has recently gained molecular support. We used genomic analyses to examine the phylogeography of a species that is critical in this debate. Carex scirpoidea Michx subsp. scirpoidea is a dioecious, amphi‐Atlantic arctic–alpine sedge that is widely distributed in North America, but absent from most of Eurasia, apart from three extremely disjunct populations in Norway, all well within the limits of the Weichselian ice sheet. Range‐wide population sampling and variation at 5,307 single nucleotide polymorphisms show that the three Norwegian populations comprise unique evolutionary lineages divergent from Greenland with high between‐population divergence. The Norwegian populations have low within‐population genetic diversity consistent with having experienced genetic bottlenecks in glacial refugia, and host private alleles that probably accumulated in long‐term isolated populations. Demographic analyses support a single, pre‐Weichselian colonization into Norway from East Greenland, and subsequent divergence of the three populations in separate refugia. Other refugial areas are identified in North‐east Greenland, Minnesota/Michigan, Colorado and Alaska. Admixed populations in British Columbia and West Greenland indicate postglacial contact. Taken together, evidence from this study strongly indicates in situ glacial survival in Scandinavia.     

Nuclear microsatellites reveal contrasting patterns of genetic structure between western and southeastern European populations of the common ash (Fraxinus excelsior L.)

To determine extant patterns of population genetic structure in common ash and gain insight into postglacial recolonization processes, we applied multilocus-based Bayesian approaches to data from 36 European populations genotyped at five nuclear microsatellite loci. We identified two contrasting patterns in terms of population genetic structure: (1) a large area from the British Isles to Lithuania throughout central Europe constituted effectively a single deme, whereas (2) strong genetic differentiation occurred over short distances in Sweden and southeastern Europe. Concomitant geographical variation was observed in estimates of allelic richness and genetic diversity, which were lowest in populations from southeastern Europe, that is, in regions close to putative ice age refuges, but high in western and central Europe, that is, in more recently recolonized areas. We suggest that in southeastern Europe, restricted postglacial gene flow caused by a rapid expansion of refuge populations in a mountainous topography is responsible for the observed strong genetic structure. In contrast, admixture of previously differentiated gene pools and high gene flow at the onset of postglacial recolonization of western and central Europe would have homogenized the genetic structure and raised the levels of genetic diversity above values in the refuges.     

Genetic Variability and Fluctuating Asymmetry of Morphological Traits in Drosophila melanogasterReared on a Pesticide-Containing Medium

The effects of chlorine-organic insecticide endosulfane (thiodan) on phenotypic and genetic variation in four morphological traits of Drosophila melanogaster(wing length, thorax length, the number of orbital bristles and the number of sternopleural bristles) were examined. In addition, the effect of this pesticide on stability of development measured as fluctuating asymmetry of bilateral traits was estimated. On the medium with endosulfane, phenotypic variation of morphometric traits was significantly higher. No difference in fluctuating asymmetry between the stressed and the control samples was found. The among-line variance of morphometric traits of flies reared on the endosulfane-containing medium was significantly higher as compared to the corresponding variance under control conditions. The efficiency of using fluctuating asymmetry and phenotypic variation of morphometric and meristic traits as indicators of environmental stress in insect populations is discussed.     

Genetic homeostasis and developmental stability in natural populations of bisexual (<Emphasis Type="Italic">Carassius auratus</Emphasis>) and unisexual (<Emphasis Type="Italic">C. gibelio</Emphasis>) goldfishes

One amphymictic diploid Carassius auratus, three apomictic triploids C. gibbelio, and a hybrid triploid Carassius auratus forms of goldfish were compared with respect to a set of characteristics of developmental stability: morphological variation, fluctuating asymmetry, and phenodeviations. It was found that C. gibbelio forms are characterized by lower levels both of morphological variation and fluctuating symmetry and more rare morphological abnormalties. A hybrid form Carassius auratus-gibelio was characterized by the intermediate values of the above-listed features. It is emphasized that phenotypic stabilization of apomictic forms is caused by two factors: its clone structure and actual developmental canalization. It is specified that promoted genetic homeostasis of the C. gibelio form does not give obvious adaptive advantages, as in the basic river basins of East and Central Europe. they are everywhere ousted by an amphimictic C. auratus.     

Biogeography meets conservation: the genetic structure of the endangered lycaenid butterfly Lycaena helle (Denis & Schiffermüller, 1775)

Cold‐adapted species are thought to have had their largest distribution ranges in central Europe during the glacial periods. Postglacial warming caused severe range shifts of such taxa into higher latitudes and altitudes. We selected the boreomontane butterfly Lycaena helle (Denis & Schiffermüller, 1775) as an example to demonstrate the genetic effects of range changes, and to document the recent status of highly fragmented remnant populations. We analysed five polymorphic microsatellite loci in 1059 individuals sampled at 50 different localities scattered over the European distribution area of the species. Genetic differentiation was strong among the mountain ranges of western Europe, but we did not detect similarly distinct genetic groups following a geographical pattern in the more eastern areas. The Fennoscandian populations form a separate genetic group, and provide evidence for a colonization from southern Finland via northern Scandinavia to south‐central Sweden. Species distribution modelling suggests a large extension of the spatial distribution during the last glacial maximum, but highlights strong retractions to a few mountain areas under current conditions. These findings, combined with our genetic data, suggest a more or less continuous distribution of L. helle throughout central Europe at the end of the last ice age. As a consequence of postglacial warming, the species retreated northwards to Fennoscandia and escaped increasing temperatures through altitudinal shifts. Therefore, the species is today restricted to population remnants located at the mountain tops of western Europe, genetically isolated from each other, and evolved into genetically unique entities. Rising temperatures and advancing habitat destruction threaten this wealth of biodiversity. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 101 , 155–168.     

No apparent reduction of gene flow in a hybrid zone between the West and North European karyotypic groups of the common shrew, Sorex araneus

The common shrew, Sorex araneus, exhibits an unusually high level of karyotypic variation. Populations with identical or similar karyotypes are defined as chromosome races, which are, in turn, grouped into larger evolutionary units, karyotypic groups. Using six microsatellite markers, we investigated the genetic structure of a hybrid zone between the Sidensjö and Abisko chromosome races, representatives of two distinct karyotypic groups believed to have been separated during the last glacial maximum, the West European karyotypic group (western group) and the North European karyotypic group (northern group), respectively. Significant F_ST values among populations suggest some weak genetic structure. All hierarchical levels show similar levels of genetic differentiation, equivalent to levels of genetic structure in several intraracial studies of common shrew populations from central Europe. Notably, genetic differentiation was of the same order of magnitude between and within karyotypic groups. Although the genetic differentiation was weak, the correlation between genetic and geographical distance was positive and significant, suggesting that the genetic variation observed between populations is a function of geographical distance rather than racial origin. Hence, considerable chromosomal differences do not seem to prevent extensive gene flow.     

Inferences of evolutionary history of a widely distributed mangrove species,Bruguiera gymnorrhiza,in the Indo‐West Pacific region

Inference of genetic structure and demographic history is fundamental issue in evolutionary biology. We examined the levels and patterns of genetic variation of a widespread mangrove species in the Indo‐West Pacific region, Bruguiera gymnorrhiza, using ten nuclear gene regions. Genetic variation of individual populations covering its distribution range was low, but as the entire species it was comparable to other plant species. Genetic differentiation among the investigated populations was high. They could be divided into two genetic clusters: the West and East clusters of the Malay Peninsula. Our results indicated that these two genetic clusters derived from their ancestral population whose effective size of which was much larger compared to the two extant clusters. The point estimate of speciation time between B. gymnorrhiza and Bruguiera sexangula was two times older than that of divergence time between the two clusters. Migration from the West cluster to the East cluster was much higher than the opposite direction but both estimated migration rates were low. The past Sundaland and/or the present Malay Peninsula are likely to prevent gene flow between the West and East clusters and function as a geographical or land barrier.     

Inbreeding reduces power-law scaling in the distribution of fluctuating asymmetry: an explanation of the basis of developmental instability

The study of fluctuating asymmetry has been controversial because of conflicting results found in much of the primary literature. It has been suggested that the source of this conflict is the fact that the basis of fluctuating asymmetry is poorly understood and that, as a consequence, methodology of fluctuating asymmetry studies may be flawed. A new model for the phenomenological basis of fluctuating asymmetry, that variation in fluctuating asymmetry is in large part due to the random exponential growth of cell populations (geometric Brownian motion) that are terminated randomly around a genetically programmed development time, is presented here. If termination of development has a genetic component, then scaling effects and kurtosis in the distribution of fluctuating asymmetry should increase with genetic redundancy of the population. This model prediction was tested by comparing the distribution of multivariate size and shape fluctuating asymmetry in large samples collected from both wild populations and four moderately inbred lines of Drosophila simulans. It was found that while wild populations were best described by a lognormal distribution with power-law scaled tails, the inbred lines derived from the wild stock were dramatically normalized (half-normal) in three of four cases. As predicted, the scaling exponent of the upper tail of the distribution of fluctuating asymmetry increased with inbreeding while the kurtosis and mean fluctuating asymmetry decreased with inbreeding. The model suggests an additional explanation of leptokurtosis in fluctuating asymmetry. Kurtosis and scaling of the statistical distribution of fluctuating asymmetry in a population is related directly to genetic differences between individuals and these differences affect their ability to buffer the process of development against random perturbations.     

High genetic diversity declines towards the geographic range periphery of Adonis vernalis,a Eurasian dry grassland plant

Genetic diversity is important for species' fitness and evolutionary processes but our knowledge on how it varies across a species' distribution range is limited. The abundant centre hypothesis (ACH) predicts that populations become smaller and more isolated towards the geographic range periphery – a pattern that in turn should be associated with decreasing genetic diversity and increasing genetic differentiation. We tested this hypothesis in Adonis vernalis, a dry grassland plant with an extensive Eurasian distribution. Its life‐history traits and distribution characteristics suggest a low genetic diversity that decreases and a high genetic differentiation that increases towards the range edge. We analysed AFLP fingerprints in 28 populations along a 4698‐km transect from the geographic range core in Russia to the western range periphery in Central and Western Europe. Contrary to our expectation, our analysis revealed high genetic diversity (range of proportion of polymorphic bands = 56–81%, He = 0.168–0.238) and low genetic differentiation across populations (Φ_ST = 0.18). However, in congruence with the genetic predictions of the ACH, genetic diversity decreased and genetic differentiation increased towards the range periphery. Spanish populations were genetically distinct, suggesting a divergent post‐glacial history in this region. The high genetic diversity and low genetic differentiation in the remaining A. vernalis populations is surprising given the species' life‐history traits and points to the possibility that the species has been widely distributed in the studied region or that it has migrated from a diverse source in an East–West direction, in the past.     

Strong genetic impoverishment from the centre of distribution in southern Europe to peripheral Baltic and isolated Scandinavian populations of the pearly heath butterfly

Aim Climatic changes and fluctuations in the past have strongly influenced the distribution of animal and plant species. Such fluctuations are also reflected in the patterns of genetic diversity on both local and global scales. The genetic pattern of the pearly heath butterfly, Coenonympha arcania, was used to evaluate the genetic differentiation of isolated (in north‐western Europe), peripheral (in north‐eastern Europe) and central (in southern Europe) populations in the context of post‐glacial distributional changes of the species. Location Europe (Sweden, Germany, the Baltic states, Italy, Slovenia, Hungary, Romania, Bulgaria). Thus, samples were collected from large parts of the species’ distribution representing the three categories mentioned above. Methods We analysed 18 loci of 569 individuals from 28 populations by allozyme electrophoresis. We used both individual‐based and population‐based analyses, including F‐statistics, various clustering methods and Markov chain Monte Carlo simulations. Results All loci, except Fum, were polymorphic. The mean F_ST for all samples was 0.18. The mean genetic distance among populations was 0.046. Two major genetic lineages were distinguished. Populations from the centre of the distributional range in southern Europe and the northern periphery of the distributional range differed significantly in their level of genetic variability. The central populations of south‐eastern Europe showed high levels of genetic diversity and no differentiation among populations. Main conclusions Most probably the two major genetic lineages evolved during glacial isolation in two disjunct Mediterranean refugia. The lack of genetic differentiation across south‐eastern Europe implies a continuous Würm ice age distribution in this area, thus supporting the functional existence of steppe forests throughout this region. The peripheral‐isolated populations in Sweden seem to have suffered from one or more severe bottlenecks, resulting in substantial genetic impoverishment. The peripheral‐connected eastern Baltic populations, on the other hand, are affected by post‐glacial and possibly recurrent gene flow from more central parts of the distribution.     

Geometric morphometrics and leaf phenotypic plasticity: assessing fluctuating asymmetry and allometry in European white oaks (Quercus)

Because plants are unable to move away from unfavourable habitats and environmental perturbations, leaf phenotypic plasticity facilitates light absorption and gas exchange. Oaks (Quercus spp.) are particularly known for their adaptability and plastic phenotypes, and leaf allometry and developmental instability may represent important mechanisms for their adaptation to environments and evolution. Because of its important role in the adaptation of plant populations to different environments, allometry can be involved in diversifying selection. Developmental instability is related to environmental perturbations and stresses by producing random deviations in structures characterized by bilateral symmetry, such as oak leaves. In addition, developmental instability can also arise from genetic bottlenecks or as a result of hybridization. The splitting of symmetric and asymmetric components of variation and their separate analysis allows the variability in leaf shape traits to be summarized, reducing the variation produced by developmental instability. The geometric morphometric approach is a useful method for the study of leaf asymmetry and allometric patterns. This method provides an important tool for the visualization of shape attributes that characterize species with highly variable leaf phenotypic patterns. In this study, leaf shape and size variability of three white oak species was investigated by means of a two‐dimensional landmark‐based method providing improved knowledge of variance partitioning, species discrimination, fluctuating asymmetry and allometric patterns of variation resulting from the different analyses. © 2015 The Linnean Society of London, Botanical Journal of the Linnean Society, 2015, 179 , 335–348.     

Chloroplast population genetics reveals low levels of genetic variation and conformation to the central–marginal hypothesis in Taxus wallichiana var. mairei,an endangered conifer endemic to China

The central–marginal hypothesis predicts that geographically peripheral populations should exhibit reduced genetic diversity and increased genetic differentiation than central populations due to smaller effective population size and stronger geographical isolation. We evaluated these predictions in the endangered conifer Taxus wallichiana var. mairei. Eight plastid simple sequence repeats (cpSSRs) were used to investigate plastid genetic variation in 22 populations of Taxus wallichiana var. mairei, encompassing nearly its entire distribution range. Low levels of plastid genetic variation and differentiation were detected in the populations, and the findings were attributed to low mutation rates, small population sizes, habitat fragmentation and isolation, and effective pollen or seed dispersal. Hunan and Hubei were identified as major refugia based on the number of private haplotypes and species distribution modeling. Trends in plastid genetic diversity and genetic differentiation from central to peripheral populations supported the predictions of the central–marginal hypothesis. In scenarios wherein the future climate becomes warmer, we predict that some peripheral populations will disappear and southern and southeastern regions will become significantly less habitable. Factors that include the levels of precipitation during the driest month, annual precipitation level, and annual temperature range will be decisive in shaping the future distribution of these populations. This study provides a theoretical basis for the conservation of T. wallichiana var. mairei.     

Global genetic diversity,lineage distribution,and Wolbachia infection of the alfalfa weevil Hypera postica (Coleoptera: Curculionidae)

The alfalfa weevil (Hypera postica) is a well‐known example of a worldwide‐distributed pest with high genetic variation. Based on the mitochondrial genes, the alfalfa weevil clusters into two main mitochondrial lineages. However, there is no clear picture of the global diversity and distribution of these lineages; neither the drivers of its diversification are known. However, it appears likely that historic demographic events including founder effects played a role. In addition, Wolbachia, a widespread intracellular parasite/symbiont, likely played an important role in the evolution of the species. Wolbachia infection so far was only detected in the Western lineage of H. postica with no information on the infecting strain, its frequency, and its consequences on the genetic diversity of the host. We here used a combination of mitochondrial and nuclear sequences of the host and sequence information on Wolbachia to document the distribution of strains and the degree of infection. The Eastern lineage has a higher genetic diversity and is found in the Mediterranean, the Middle East, Eastern Europe, and eastern America, whereas the less diverse Western lineage is found in Central Europe and the western America. Both lineages are infected with the same common strain of Wolbachia belonging to Supergroup B. Based on neutrality tests, selection tests, and the current distribution and diversification of Wolbachia in H. postica, we suggested the Wolbachia infection did not shape genetic diversity of the host. The introduced populations in the United States are generally genetically less diverse, which is in line with founder effects.     

GENETIC VARIATION IN BROMUS TECTORUM (POACEAE): POPULATION DIFFERENTIATION IN ITS NORTH AMERICAN RANGE

Allelic variation in seedlings from 60 North American populations of the alien annual grass Bromus tectorum was determined at 25 loci using starch gel electrophoresis. Populations were collected from four regions; east of the Rocky Mountains, Nevada and California, the Intermountain West, and British Columbia. Compared to other diploid seed plants, genetic variation within these populations of B. tectorum is low: 4.60% of loci are polymorphic per population, with an average of 1.05 alleles per locus and a mean expected heterozygosity of 0.012. Although 2,141 individuals were analyzed, no heterozygous individuals were detected, and consequently, mean observed heterozygosity is 0.000. Extensive deviations from Hardy-Weinberg expectations were observed at every polymorphic locus due to heterozygote deficiencies. The mean genetic identity (Nei's I) between population pairs was 0.980 and indicates a high level of overall genetic similarity among populations. The among-population component of the total gene diversity is high (G_ST ⁼ 0.478), indicating substantial genetic differentiation among populations. These results are consistent with previous reports for highly self-pollinating plants of low genetic variation and substantial genetic differentiation among populations. Despite the lack of genetic variation as measured by enzyme electrophoresis, this weedy grass has become exceedingly abundant in a diverse array of arid environments throughout much of western North America, perhaps due to phenotypic plasticity.     

Geographic patterns in morphometric and genetic variation for coyote populations with emphasis on southeastern coyotes

Prior to 1900, coyotes (Canis latrans) were restricted to the western and central regions of North America, but by the early 2000s, coyotes became ubiquitous throughout the eastern United States. Information regarding morphological and genetic structure of coyote populations in the southeastern United States is limited, and where data exist, they are rarely compared to those from other regions of North America. We assessed geographic patterns in morphology and genetics of coyotes with special consideration of coyotes in the southeastern United States. Mean body mass of coyote populations increased along a west‐to‐east gradient, with southeastern coyotes being intermediate to western and northeastern coyotes. Similarly, principal component analysis of body mass and linear body measurements suggested that southeastern coyotes were intermediate to western and northeastern coyotes in body size but exhibited shorter tails and ears from other populations. Genetic analyses indicated that southeastern coyotes represented a distinct genetic cluster that differentiated strongly from western and northeastern coyotes. We postulate that southeastern coyotes experienced lower immigration from western populations than did northeastern coyotes, and over time, genetically diverged from both western and northeastern populations. Coyotes colonizing eastern North America experienced different selective pressures than did stable populations in the core range, and we offer that the larger body size of eastern coyotes reflects an adaptation that improved dispersal capabilities of individuals in the expanding range.     

Divergent lineages in a semi‐arid mallee species,Eucalyptus behriana,correspond to a major geographic break in southeastern Australia

AimTo infer relationships between populations of the semi‐arid, mallee eucalypt, Eucalyptus behriana, to build hypotheses regarding evolution of major disjunctions in the species'' distribution and to expand understanding of the biogeographical history of southeastern Australia.LocationSoutheastern Australia.Taxon Eucalyptus behriana (Myrtaceae, Angiospermae).MethodsWe developed a large dataset of anonymous genomic loci for 97 samples from 11 populations of E. behriana using double digest restriction site‐associated DNA sequencing (ddRAD‐seq), to determine genetic relationships between the populations. These relationships, along with species distribution models, were used to construct hypotheses regarding environmental processes that have driven fragmentation of the species’ distribution.ResultsGreatest genetic divergence was between populations on either side of the Lower Murray Basin. Populations west of the Basin showed greater genetic divergence between one another than the eastern populations. The most genetically distinct population in the east (Long Forest) was separated from others by the Great Dividing Range. A close relationship was found between the outlying northernmost population (near West Wyalong) and those in the Victorian Goldfields despite a large disjunction between them.ConclusionsPatterns of genetic variation are consistent with a history of vicariant differentiation of disjunct populations. We infer that an early disjunction to develop in the species distribution was that across the Lower Murray Basin, an important biogeographical barrier separating many dry sclerophyll plant taxa in southeastern Australia. Additionally, our results suggest that the western populations fragmented earlier than the eastern ones. Fragmentation, both west and east of the Murray Basin, is likely tied to climatic changes associated with glacial‐interglacial cycles although it remains possible that major geological events including uplift of the Mount Lofty Ranges and basalt flows in the Newer Volcanics Province also played a role.