Genetic Structure of the Western and Eastern African Sahel/Savannah Belt and the Role of Nomadic Pastoralists as Inferred from the Variation of D-Loop Mitochondrial DNA Sequences

The objective of this study was to provide deeper knowledge of the maternal genetic structure and demographic history of the human populations of the Sahel/Savannah belt, the extensive region lying between the Sahara and tropical rainforests, spanning from the Atlantic Ocean to the Red Sea coast. The study aimed to confirm or disconfirm archaeological and linguistic data indicating that the region's populations underwent diversification as a result of the spread of agropastoral food-producing subsistence lifestyles, over time dividing the region into separate areas of nomadic pastoralism, on the one hand, and sedentary farming, on the other. To perform both descriptive and coalescence analyses from the Sahel/Savannah belt's entire region, including western and eastern rather than just central populations studied previously, we generated a new mitochondrial DNA (mtDNA) data set not only having almost 2,000 samples (875 of which were newly collected) but also encompassing whole mtDNA D-loop segment rather than only the previously studied hypervariable segment 1. While comparing our analyses with previous results from the Lake Chad Basin (central Sahel/Savannah Belt), we found similar intrapopulation diversity measures (i.e., lower values in pastoralists than in farmers). However, the new data set pointed to significant differences in mating strategies between western and eastern pastoralists: our results suggest higher gene flow between the Arabic pastoralists and neighboring farmers in the eastern part than between the Fulani pastoralists and their sedentary neighbors in the western part of the Sahel/Savannah belt. The findings are discussed in light of archaeological and linguistic data, allowing us to postulate that the genetic differentiation of Fulani pastoralists from the common western African agropastoral gene pool occurred at around the same time as the arrival of the Arabic pastoralists to eastern Africa. However, it seems that while the process of divergence of the Fulani pastoralists in the west was accompanied by a loss of Fulani females to other populations, the Arab pastoralists' immigration to the Sahel/Savannah belt conversely resulted in some gain of local females into this Arab population.

Population genetics of Southern Hemisphere tope shark (Galeorhinus galeus): Intercontinental divergence and constrained gene flow at different geographical scales

The tope shark (Galeorhinus galeus Linnaeus, 1758) is a temperate, coastal hound shark found in the Atlantic and Indo-Pacific oceans. In this study, the population structure of Galeorhinus galeus was determined across the entire Southern Hemisphere, where the species is heavily targeted by commercial fisheries, as well as locally, along the South African coastline. Analysis was conducted on a total of 185 samples using 19 microsatellite markers and a 671 bp fragment of the NADH dehydrogenase subunit 2 (ND2) gene. Across the Southern Hemisphere, three geographically distinct clades were recovered, including one from South America (Argentina, Chile), one from Africa (all the South African collections) and an Australia-New Zealand clade. Nuclear data revealed significant population subdivisions (F_ST = 0.192 to 0.376, p<0.05) indicating limited gene flow for tope sharks across ocean basins. Marked population connectivity was however evident across the Indian Ocean based on Bayesian clustering analysis. More locally in South Africa, F-statistics and multivariate analysis supported moderate to high gene flow across the Atlantic/Indian Ocean boundary (F_ST = 0.035 to 0.044, p<0.05), with exception of samples from Struisbaai and Port Elizabeth which differed significantly from the rest. Discriminant and Bayesian clustering analysis indicated admixture in all sampling populations, decreasing from west to east, corroborating possible restriction to gene flow across regional oceanographic barriers. Mitochondrial sequence data recovered seven haplotypes (h = 0.216, π = 0.001) for South Africa, with one major haplotype shared by 87% of the individuals and at least one private haplotype for each sampling location except Port Elizabeth. As with many other coastal shark species with cosmopolitan distribution, this study confirms the lack of both historical dispersal and inter-oceanic gene flow while also implicating contemporary factors such as oceanic currents and thermal fronts to drive local genetic structure of G. galeus on a smaller spatial scale.

Assessment of the potential for gene flow from transgenic maize (Zea mays L.) to eastern gamagrass (Tripsacum dactyloides L.)

Eastern gamagrass (Tripsacum dactyloides L.) belongs to the same tribe of the Poaceae family as maize (Zea mays L.) and grows naturally in the same region where maize is commercially produced in the USA. Although no evidence exists of gene flow from maize to eastern gamagrass in nature, experimental crosses between the two species were produced using specific techniques. As part of environmental risk assessment, the possibility of transgene flow from maize to eastern gamagrass populations in nature was evaluated with the objectives: (1) to assess the seeds of eastern gamagrass populations naturally growing near commercial maize fields for the presence of a transgenic glyphosate-tolerance gene (cp4 epsps) that would indicate cross-pollination between the two species, and (2) to evaluate the possibility of interspecific hybridization between transgenic maize used as male parent and eastern gamagrass used as female parent. A total of 46,643 seeds from 54 eastern gamagrass populations collected in proximity of maize fields in Illinois, USA were planted in a field in 2014 and 2015. Emerged seedlings were treated with glyphosate herbicide and assessed for survival. An additional 48,000 seeds from the same 54 eastern gamagrass populations were tested for the presence of the cp4 epsps transgene markers using TaqMan® PCR method. The results from these trials showed that no seedlings survived the herbicide treatment and no seed indicated presence of the herbicide tolerant cp4 epsps transgene, even though these eastern gamagrass populations were exposed to glyphosate-tolerant maize pollen for years. Furthermore, no interspecific hybrid seeds were produced from 135 hand-pollination attempts involving 1529 eastern gamagrass spikelets exposed to maize pollen. Together, these results indicate that there is no evidence of gene flow from maize to eastern gamagrass in natural habitats. The outcome of this study should be taken in consideration when assessing for environmental risks regarding the consequence of gene flow from transgenic maize to its wild relatives.

No evidence of contemporary interploidy gene flow between the closely related European woodland violets Viola reichenbachiana and V. riviniana (sect. Viola, Violaceae)

Viola reichenbachiana (2n = 4x = 20) and V. riviniana (2n = 8x = 40) are closely related species widely distributed in Europe, often sharing the same habitat throughout their overlapping ranges. It has been suggested in numerous studies that their high intraspecific morphological variability and plasticity might have been further increased by interspecific hybridisation in contact zones, given the sympatry of the species and the incomplete sterility of their hybrid. The aims of this study were to: (i) confirm that V. reichenbachiana and V. riviniana have one 4x genome in common, and (ii) determine the impact of hybridisation and introgression on genetic variation of these two species in selected European populations. For our study, we used 31 Viola populations from four European countries, which were analysed using AFLP and sequencing of a variable plastid intergenic spacer, trnH-psbA. Our analysis revealed that V. reichenbachiana exhibited larger haplotype diversity, having three species-specific haplotypes versus one in V. riviniana. The relationships among haplotypes suggest transfer of common haplotypes into V. riviniana from both V. reichenbachiana and hypothetically the other, now extinct, parental species. AFLP analysis showed low overall genetic diversity of both species, with V. riviniana showing higher among-population diversity. None of the morphologically designated hybrid populations had additive AFLP polymorphisms that would have indicated recent hybridisation. Also, kinship coefficients between both species did not indicate gene flow. V. riviniana showed significant population subdivision and significant isolation by distance, in contrast to V. reichenbachiana. The results indicate lack of gene flow between species, high influence of selfing on genetic variability, as well as probably only localised introgression toward V. riviniana.

The Rice Paradox: Multiple Origins but Single Domestication in Asian Rice

The origin of domesticated Asian rice (Oryza sativa) has been a contentious topic, with conflicting evidence for either single or multiple domestication of this key crop species. We examined the evolutionary history of domesticated rice by analyzing de novo assembled genomes from domesticated rice and its wild progenitors. Our results indicate multiple origins, where each domesticated rice subpopulation (japonica, indica, and aus) arose separately from progenitor O. rufipogon and/or O. nivara. Coalescence-based modeling of demographic parameters estimate that the first domesticated rice population to split off from O. rufipogon was O. sativa ssp. japonica, occurring at ∼13.1-24.1 ka, which is an order of magnitude older then the earliest archeological date of domestication. This date is consistent, however, with the expansion of O. rufipogon populations after the Last Glacial Maximum ∼18 ka and archeological evidence for early wild rice management in China. We also show that there is significant gene flow from japonica to both indica (∼17%) and aus (∼15%), which led to the transfer of domestication alleles from early-domesticated japonica to proto-indica and proto-aus populations. Our results provide support for a model in which different rice subspecies had separate origins, but that de novo domestication occurred only once, in O. sativa ssp. japonica, and introgressive hybridization from early japonica to proto-indica and proto-aus led to domesticated indica and aus rice.

A Working Model of the Deep Relationships of Diverse Modern Human Genetic Lineages Outside of Africa

A major topic of interest in human prehistory is how the large-scale genetic structure of modern populations outside of Africa was established. Demographic models have been developed that capture the relationships among small numbers of populations or within particular geographical regions, but constructing a phylogenetic tree with gene flow events for a wide diversity of non-Africans remains a difficult problem. Here, we report a model that provides a good statistical fit to allele-frequency correlation patterns among East Asians, Australasians, Native Americans, and ancient western and northern Eurasians, together with archaic human groups. The model features a primary eastern/western bifurcation dating to at least 45,000 years ago, with Australasians nested inside the eastern clade, and a parsimonious set of admixture events. While our results still represent a simplified picture, they provide a useful summary of deep Eurasian population history that can serve as a null model for future studies and a baseline for further discoveries.

The relationship between least-cost and resistance distance

Least-cost modelling and circuit theory are common analogs used in ecology and evolution to model gene flow or animal movement across landscapes. Least-cost modelling estimates the least-cost distance, whereas circuit theory estimates resistance distance. The bias added in choosing one method over the other has not been well documented. We designed an experiment to test whether both methods were linearly related. We also tested the sensitivity of these metrics to variation in Euclidean distance, spatial autocorrelation, the number of pixels representing the landscape, and data aggregation. We found that least-cost and resistance distance were not linearly related unless a transformation was applied. Resistance distance was less sensitive to the number of pixels representing a landscape and was also less sensitive than least-cost distance to the Euclidean distance between nodes. Spatial autocorrelation did not affect either method or the relationship between methods. Resistance distance was more sensitive to aggregation in any form compared to least-cost distance. Therefore, the metric used to infer movement or gene flow and the manipulations applied to the data used to calculate these metrics may govern findings.

Rethinking the influence of hydroelectric development on gene flow in a long-lived fish, the Lake Sturgeon Acipenser fulvescens

Many hydroelectric dams have been in place for 50 - >100 years, which for most fish species means that enough generations have passed for fragmentation induced divergence to have accumulated. However, for long-lived species such as Lake Sturgeon, Acipenser fulvescens, it should be possible to discriminate between historical population structuring and contemporary gene flow and improve the broader understanding of anthropogenic influence. On the Winnipeg River, Manitoba, two hypotheses were tested: 1) Measureable quantities of former reservoir dwelling Lake Sturgeon now reside downstream of the Slave Falls Generating Station, and 2) genetically differentiated populations of Lake Sturgeon occur upstream and downstream, a result of historical structuring. Genetic methods based on ten microsatellite markers were employed, and simulations were conducted to provide context. With regards to contemporary upstream to downstream contributions, the inclusion of length-at-age data proved informative. Both pairwise relatedness and Bayesian clustering analysis substantiated that fast-growing outliers, apparently entrained after residing in the upstream reservoir for several years, accounted for ~15% of the Lake Sturgeon 525–750 mm fork length captured downstream. With regards to historical structuring, upstream and downstream populations were found to be differentiated (F_ST = 0.011, and 0.013–0.014 when fast-growing outliers were excluded), and heterozygosity metrics were higher for downstream versus upstream juveniles. Historical asymmetric (downstream) gene flow in the vicinity of the generating station was the most logical explanation for the observed genetic structuring. In this section of the Winnipeg River, construction of a major dam does not appear to have fragmented a previously panmictic Lake Sturgeon population, but alterations to habitat may be influencing upstream to downstream contributions in unexpected ways.

Reduced representation genome sequencing reveals patterns of genetic diversity and selection in apple

Identifying DNA sequence variations is a fundamental step towards deciphering the genetic basis of traits of interest. Here, a total of 20 cultivated and 10 wild apples were genotyped using specific-locus amplified fragment sequencing, and 39,635 single nucleotide polymorphisms with no missing genotypes and evenly distributed along the genome were selected to investigate patterns of genome-wide genetic variations between cultivated and wild apples. Overall, wild apples displayed higher levels of genetic diversity than cultivated apples. Linkage disequilibrium (LD) decays were observed quite rapidly in cultivated and wild apples, with an r² -value below 0.2 at 440 and 280 bp, respectively. Moreover, bidirectional gene flow and different distribution patterns of LD blocks were detected between domesticated and wild apples. Most LD blocks unique to cultivated apples were located within QTL regions controlling fruit quality, thus suggesting that fruit quality had probably undergone selection during apple domestication. The genome of the earliest cultivated apple in China, Nai, was highly similar to that of Malus sieversii, and contained a small portion of genetic material from other wild apple species. This suggested that introgression could have been an important driving force during initial domestication of apple. These findings will facilitate future breeding and genetic dissection of complex traits in apple.

Pollen-mediated gene flow and fine-scale spatial genetic structure in Olea europaea subsp. europaea var. sylvestris

BACKGROUND AND AIMS

Wild olive ( Olea europaea subsp. europaea var. sylvestris ) is important from an economic and ecological point of view. The effects of anthropogenic activities may lead to the genetic erosion of its genetic patrimony, which has high value for breeding programmes. In particular, the consequences of the introgression from cultivated stands are strongly dependent on the extent of gene flow and therefore this work aims at quantitatively describing contemporary gene flow patterns in wild olive natural populations.

METHODS

The studied wild population is located in an undisturbed forest, in southern Spain, considered one of the few extant hotspots of true oleaster diversity. A total of 225 potential father trees and seeds issued from five mother trees were genotyped by eight microsatellite markers. Levels of contemporary pollen flow, in terms of both pollen immigration rates and within-population dynamics, were measured through paternity analyses. Moreover, the extent of fine-scale spatial genetic structure (SGS) was studied to assess the relative importance of seed and pollen dispersal in shaping the spatial distribution of genetic variation.

KEY RESULTS

The results showed that the population under study is characterized by a high genetic diversity, a relatively high pollen immigration rate (0·57), an average within-population pollen dispersal of about 107 m and weak but significant SGS up to 40 m. The population is a mosaic of several intermingled genetic clusters that is likely to be generated by spatially restricted seed dispersal. Moreover, wild oleasters were found to be self-incompatible and preferential mating between some genotypes was revealed.

CONCLUSIONS

Knowledge of the within-population genetic structure and gene flow dynamics will lead to identifying possible strategies aimed at limiting the effect of anthropogenic activities and improving breeding programmes for the conservation of olive tree forest genetic resources.

Clustering of 770,000 genomes reveals post-colonial population structure of North America

Despite strides in characterizing human history from genetic polymorphism data, progress in identifying genetic signatures of recent demography has been limited. Here we identify very recent fine-scale population structure in North America from a network of over 500 million genetic (identity-by-descent, IBD) connections among 770,000 genotyped individuals of US origin. We detect densely connected clusters within the network and annotate these clusters using a database of over 20 million genealogical records. Recent population patterns captured by IBD clustering include immigrants such as Scandinavians and French Canadians; groups with continental admixture such as Puerto Ricans; settlers such as the Amish and Appalachians who experienced geographic or cultural isolation; and broad historical trends, including reduced north-south gene flow. Our results yield a detailed historical portrait of North America after European settlement and support substantial genetic heterogeneity in the United States beyond that uncovered by previous studies.

Limited Genetic Connectivity between Gorgonian Morphotypes along a Depth Gradient

Gorgonian species show a high morphological variability in relation to the environment in which they live. In coastal areas, parameters such as temperature, light, currents, and food availability vary significantly with depth, potentially affecting morphology of the colonies and the structure of the populations, as well as their connectivity patterns. In tropical seas, the existence of connectivity between shallow and deep populations supported the hypothesis that the deep coral reefs could potentially act as (reproductive) refugia fostering re-colonization of shallow areas after mortality events. Moreover, this hypothesis is not so clear accepted in temperate seas. Eunicella singularis is one of the most common gorgonian species in Northwestern Mediterranean Sea, playing an important role as ecosystem engineer by providing biomass and complexity to the coralligenous habitats. It has a wide bathymetric distribution ranging from about 10 m to 100 m. Two depth-related morphotypes have been identified, differing in colony morphology, sclerite size and shape, and occurrence of symbiotic algae, but not in mitochondrial DNA haplotypes. In the present study the genetic structure of E. singularis populations along a horizontal and bathymetric gradient was assessed using microsatellites and ITS1 sequences. Restricted gene flow was found at 30-40 m depth between the two Eunicella morphotypes. Conversely, no genetic structuring has been found among shallow water populations within a spatial scale of ten kilometers. The break in gene flow between shallow and deep populations contributes to explain the morphological variability observed at different depths. Moreover, the limited vertical connectivity hinted that the refugia hypothesis does not apply to E. singularis. Re-colonization of shallow water populations, occasionally affected by mass mortality events, should then be mainly fueled by larvae from other shallow water populations.

Multispecies Outcomes of Sympatric Speciation after Admixture with the Source Population in Two Radiations of Nicaraguan Crater Lake Cichlids

The formation of species in the absence of geographic barriers (i.e. sympatric speciation) remains one of the most controversial topics in evolutionary biology. While theoretical models have shown that this most extreme case of primary divergence-with-gene-flow is possible, only a handful of accepted empirical examples exist. And even for the most convincing examples uncertainties remain; complex histories of isolation and secondary contact can make species falsely appear to have originated by sympatric speciation. This alternative scenario is notoriously difficult to rule out. Midas cichlids inhabiting small and remote crater lakes in Nicaragua are traditionally considered to be one of the best examples of sympatric speciation and lend themselves to test the different evolutionary scenarios that could lead to apparent sympatric speciation since the system is relatively small and the source populations known. Here we reconstruct the evolutionary history of two small-scale radiations of Midas cichlids inhabiting crater lakes Apoyo and Xiloá through a comprehensive genomic data set. We find no signs of differential admixture of any of the sympatric species in the respective radiations. Together with coalescent simulations of different demographic models our results support a scenario of speciation that was initiated in sympatry and does not result from secondary contact of already partly diverged populations. Furthermore, several species seem to have diverged simultaneously, making Midas cichlids an empirical example of multispecies outcomes of sympatric speciation. Importantly, however, the demographic models strongly support an admixture event from the source population into both crater lakes shortly before the onset of the radiations within the lakes. This opens the possibility that the formation of reproductive barriers involved in sympatric speciation was facilitated by genetic variants that evolved in a period of isolation between the initial founding population and the secondary migrants that came from the same source population. Thus, the exact mechanisms by which these species arose might be different from what had been thought before.

Population Structure and Evolution after Speciation of the Hokkaido Salamander (Hynobius retardatus)

The Hokkaido salamander (Hynobius retardatus) is endemic to Hokkaido Island, Japan, and shows intriguing flexible phenotypic plasticity and regional morphological diversity. However, to date, allozymes and partial mitochondria DNA sequences have provided only an outline of its demographic histories and the pattern of its genetic diversification. To understand the finer details of the population structure of this species and its evolution since speciation, we genotyped five regional populations by using 12 recently developed microsatellite polymorphic markers. We found a clear population structure with low gene flow among the five populations, but a close genetic relationship between the Teshio and Kitami populations. Our demographic analysis suggested that Teshio and Erimo had the largest effective population sizes among the five populations. These findings regarding the population structure and demography of H. retardatus improve our understanding of the faunal phylogeography on Hokkaido Island and also provide fundamental genetic information that will be useful for future studies.

Level of tissue differentiation influences the activation of a heat-inducible flower-specific system for genetic containment in poplar (Populus tremula L.)

Differentiation level but not transgene copy number influenced activation of a gene containment system in poplar. Heat treatments promoted CRE gene body methylation. The flower-specific transgene deletion was confirmed. Gene flow between genetic modified trees and their wild relatives is still motive of concern. Therefore, approaches for gene containment are required. In this study, we designed a novel strategy for achieving an inducible and flower-specific transgene removal from poplar trees but still expressing the transgene in the plant body. Hence, pollen carrying transgenes could be used for breeding purposes under controlled conditions in a first phase, and in the second phase genetic modified poplars developing transgene-free pollen grains could be released. This approach is based on the recombination systems CRE/loxP and FLP/frt. Both gene constructs contained a heat-inducible CRE/loxP-based spacer sequence for in vivo assembling of the flower-specific FLP/frt system. This allowed inducible activation of gene containment. The FLP/frt system was under the regulation of a flower-specific promoter, either CGPDHC or PTD. Our results confirmed complete CRE/loxP-based in vivo assembling of the flower-specific transgene excision system after heat treatment in all cells for up to 30 % of regenerants derived from undifferentiated tissue cultures. Degradation of HSP::CRE/loxP spacer after recombination but also persistence as extrachromosomal DNA circles were detected in sub-lines obtained after heat treatments. Furthermore, heat treatment promoted methylation of the CRE gene body. A lower methylation level was detected at CpG sites in transgenic sub-lines showing complete CRE/loxP recombination and persistence of CRE/loxP spacer, compared to sub-lines with incomplete recombination. However, our results suggest that low methylation might be necessary but not sufficient for recombination. The flower-specific FLP/frt-based transgene deletion was confirmed in 6.3 % of flowers.

A Spatial Framework for Understanding Population Structure and Admixture

Geographic patterns of genetic variation within modern populations, produced by complex histories of migration, can be difficult to infer and visually summarize. A general consequence of geographically limited dispersal is that samples from nearby locations tend to be more closely related than samples from distant locations, and so genetic covariance often recapitulates geographic proximity. We use genome-wide polymorphism data to build "geogenetic maps," which, when applied to stationary populations, produces a map of the geographic positions of the populations, but with distances distorted to reflect historical rates of gene flow. In the underlying model, allele frequency covariance is a decreasing function of geogenetic distance, and nonlocal gene flow such as admixture can be identified as anomalously strong covariance over long distances. This admixture is explicitly co-estimated and depicted as arrows, from the source of admixture to the recipient, on the geogenetic map. We demonstrate the utility of this method on a circum-Tibetan sampling of the greenish warbler (Phylloscopus trochiloides), in which we find evidence for gene flow between the adjacent, terminal populations of the ring species. We also analyze a global sampling of human populations, for which we largely recover the geography of the sampling, with support for significant histories of admixture in many samples. This new tool for understanding and visualizing patterns of population structure is implemented in a Bayesian framework in the program SpaceMix.

Home Bodies and Wanderers: Sympatric Lineages of the Deep-Sea Black Coral Leiopathes glaberrima

Colonial corals occur in a wide range of marine benthic habitats from the shallows to the deep ocean, often defining the structure of their local community. The black coral Leiopathes glaberrima is a long-lived foundation species occurring on carbonate outcrops in the Northern Gulf of Mexico (GoM). Multiple color morphs of L. glaberrima grow sympatrically in the region. Morphological, mitochondrial and nuclear ribosomal markers supported the hypothesis that color morphs constituted a single biological species and that colonies, regardless of color, were somewhat genetically differentiated east and west of the Mississippi Canyon. Ten microsatellite loci were used to determine finer-scale population genetic structure and reproductive characteristics. Gene flow was disrupted between and within two nearby (distance = 36.4 km) hardground sites and two sympatric microsatellite lineages, which might constitute cryptic species, were recovered. Lineage one was outbred and found in all sampled locations (N = 5) across 765.6 km in the Northern Gulf of Mexico. Lineage two was inbred, reproducing predominantly by fragmentation, and restricted to sites around Viosca Knoll. In these sites the lineages and the color phenotypes occurred in different microhabitats, and models of maximum entropy suggested that depth and slope influence the distribution of the color phenotypes within the Vioska Knolls. We conclude that L. glaberrima is phenotypically plastic with a mixed reproductive strategy in the Northern GoM. Such strategy might enable this long-lived species to balance local recruitment with occasional long-distance dispersal to colonize new sites in an environment where habitat is limited.

Ancient Humans Influenced the Current Spatial Genetic Structure of Common Walnut Populations in Asia

Common walnut (Juglans regia L) is an economically important species cultivated worldwide for its wood and nuts. It is generally accepted that J. regia survived and grew spontaneously in almost completely isolated stands in its Asian native range after the Last Glacial Maximum. Despite its natural geographic isolation, J. regia evolved over many centuries under the influence of human management and exploitation. We evaluated the hypothesis that the current distribution of natural genetic resources of common walnut in Asia is, at least in part, the product of ancient anthropogenic dispersal, human cultural interactions, and afforestation. Genetic analysis combined with ethno-linguistic and historical data indicated that ancient trade routes such as the Persian Royal Road and Silk Road enabled long-distance dispersal of J. regia from Iran and Trans-Caucasus to Central Asia, and from Western to Eastern China. Ancient commerce also disrupted the local spatial genetic structure of autochthonous walnut populations between Tashkent and Samarkand (Central-Eastern Uzbekistan), where the northern and central routes of the Northern Silk Road converged. A significant association between ancient language phyla and the genetic structure of walnut populations is reported even after adjustment for geographic distances that could have affected both walnut gene flow and human commerce over the centuries. Beyond the economic importance of common walnut, our study delineates an alternative approach for understanding how the genetic resources of long-lived perennial tree species may be affected by the interaction of geography and human history.

Malaria Molecular Epidemiology: Lessons from the International Centers of Excellence for Malaria Research Network

Molecular epidemiology leverages genetic information to study the risk factors that affect the frequency and distribution of malaria cases. This article describes molecular epidemiologic investigations currently being carried out by the International Centers of Excellence for Malaria Research (ICEMR) network in a variety of malaria-endemic settings. First, we discuss various novel approaches to understand malaria incidence and gametocytemia, focusing on Plasmodium falciparum and Plasmodium vivax. Second, we describe and compare different parasite genotyping methods commonly used in malaria epidemiology and population genetics. Finally, we discuss potential applications of molecular epidemiological tools and methods toward malaria control and elimination efforts.

Inferring Population Genetic Structure in Widely and Continuously Distributed Carnivores: The Stone Marten (Martes foina) as a Case Study

The stone marten is a widely distributed mustelid in the Palaearctic region that exhibits variable habitat preferences in different parts of its range. The species is a Holocene immigrant from southwest Asia which, according to fossil remains, followed the expansion of the Neolithic farming cultures into Europe and possibly colonized the Iberian Peninsula during the Early Neolithic (ca. 7,000 years BP). However, the population genetic structure and historical biogeography of this generalist carnivore remains essentially unknown. In this study we have combined mitochondrial DNA (mtDNA) sequencing (621 bp) and microsatellite genotyping (23 polymorphic markers) to infer the population genetic structure of the stone marten within the Iberian Peninsula. The mtDNA data revealed low haplotype and nucleotide diversities and a lack of phylogeographic structure, most likely due to a recent colonization of the Iberian Peninsula by a few mtDNA lineages during the Early Neolithic. The microsatellite data set was analysed with a) spatial and non-spatial Bayesian individual-based clustering (IBC) approaches (STRUCTURE, TESS, BAPS and GENELAND), and b) multivariate methods [discriminant analysis of principal components (DAPC) and spatial principal component analysis (sPCA)]. Additionally, because isolation by distance (IBD) is a common spatial genetic pattern in mobile and continuously distributed species and it may represent a challenge to the performance of the above methods, the microsatellite data set was tested for its presence. Overall, the genetic structure of the stone marten in the Iberian Peninsula was characterized by a NE-SW spatial pattern of IBD, and this may explain the observed disagreement between clustering solutions obtained by the different IBC methods. However, there was significant indication for contemporary genetic structuring, albeit weak, into at least three different subpopulations. The detected subdivision could be attributed to the influence of the rivers Ebro, Tagus and Guadiana, suggesting that main watercourses in the Iberian Peninsula may act as semi-permeable barriers to gene flow in stone martens. To our knowledge, this is the first phylogeographic and population genetic study of the species at a broad regional scale. We also wanted to make the case for the importance and benefits of using and comparing multiple different clustering and multivariate methods in spatial genetic analyses of mobile and continuously distributed species.

A Genome Scan for Genes Underlying Microgeographic-Scale Local Adaptation in a Wild Arabidopsis Species

Adaptive divergence at the microgeographic scale has been generally disregarded because high gene flow is expected to disrupt local adaptation. Yet, growing number of studies reporting adaptive divergence at a small spatial scale highlight the importance of this process in evolutionary biology. To investigate the genetic basis of microgeographic local adaptation, we conducted a genome-wide scan among sets of continuously distributed populations of Arabidopsis halleri subsp. gemmifera that show altitudinal phenotypic divergence despite gene flow. Genomic comparisons were independently conducted in two distinct mountains where similar highland ecotypes are observed, presumably as a result of convergent evolution. Here, we established a de novo reference genome and employed an individual-based resequencing for a total of 56 individuals. Among 527,225 reliable SNP loci, we focused on those showing a unidirectional allele frequency shift across altitudes. Statistical tests on the screened genes showed that our microgeographic population genomic approach successfully retrieve genes with functional annotations that are in line with the known phenotypic and environmental differences between altitudes. Furthermore, comparison between the two distinct mountains enabled us to screen out those genes that are neutral or adaptive only in either mountain, and identify the genes involved in the convergent evolution. Our study demonstrates that the genomic comparison among a set of genetically connected populations, instead of the commonly-performed comparison between two isolated populations, can also offer an effective screening for the genetic basis of local adaptation.

Virulence Structure of Blumeria graminis f. sp. tritici and Its Genetic Diversity by ISSR and SRAP Profiling Analyses

Blumeria graminis f. sp. tritici, which causes wheat powdery mildew, is an obligate biotrophic pathogen that can easily genetically adapt to its host plant. Understanding the virulence structure of and genetic variations in this pathogen is essential for disease control and for breeding resistance to wheat powdery mildew. This study investigated 17 pathogenic populations in Sichuan, China and classified 109 isolates into two distinct groups based on pathogenicity analysis: high virulence (HV, 92 isolates) and low virulence (LV, 17 isolates). Populations from Yibin (Southern region), Xichang (Western region), and Meishan (Middle region) showed lower virulence frequencies than populations from other regions. Many of the previously known resistance genes did not confer resistance in this study. The resistance gene Pm21 displayed an immune response to pathogenic challenge with all populations in Sichuan, and Pm13, Pm5b, Pm2+6, and PmXBD maintained resistance. AMOVA revealed significantly higher levels of variation within populations and lower levels of variation among populations within regions. High levels of gene flow were detected among populations in the four regions. Closely related populations within each region were distinguished by cluster analyses using ISSR and SRAP alleles. Both ISSR and SRAP allele profiling analyses revealed high levels of genetic diversity among pathogenic populations in Sichuan. Although ISSR and SRAP profiling analysis showed similar resolutions, the SRAP alleles appeared to be more informative. We did not detect any significant association between these alleles and the virulence or pathogenicity of the pathogen. Our results suggest that ISSR and SRAP alleles are more efficient for the characterization of small or closely related populations versus distantly related populations.

Refugia persistence of Qinghai-Tibetan plateau by the cold-tolerant bird Tetraogallus tibetanus (Galliformes: Phasianidae)

Most of the temperate species are expected to have moved to lower altitudes during the glacial periods of the Quaternary. Here we tested this hypothesis in a cold-tolerant avian species Tibetan snowcock (Tetraogallus tibetanus) using two segments of mitochondrial gene (a 705bp Cytochrome-b; abbrev. Cyt-b and an 854 bp Control Region; abbrev. CR) and eight microsatellite loci by characterizing population differentiation and gene flow across its range. Combined (Cyt-b + CR) datasets detected several partially lineages with poor support. Microsatellite data, however, identified two distinct lineages congruent with the geographically separated western and central regions of Qinghai-Tibetan Plateau (QTP). The phylogeographic patterns that we observed might be explained by a combination of vicariance events that led to local isolation of T. tibetanus during warm periods and range expansions and population intermixing during cold periods. The results of this study add to our knowledge of population differentiation and connectivity in high altitude mountain ecosystems.

Phylogeography of Rhodiola kirilowii (Crassulaceae): a story of Miocene divergence and quaternary expansion

The evolution and current distribution of the Sino-Tibetan flora have been greatly affected by historical geological events, such as the uplift of the Qinghai-Tibetan Plateau (QTP), and Quaternary climatic oscillations. Rhodiola kirilowii, a perennial herb with its distribution ranging from the southeastern QTP and the Hengduan Mountains (HM) to adjacent northern China and central Asia, provides an excellent model to examine and disentangle the effect of both geological orogeny and climatic oscillation on the evolutionary history of species with such distribution patterns. We here conducted a phylogeographic study using sequences of two chloroplast fragments (trnL-F and trnS-G) and internal transcribed spacers in 29 populations of R. kirilowii. A total of 25 plastid haplotypes and 12 ITS ribotypes were found. Molecular clock estimation revealed deep divergence between the central Asian populations and other populations from the HM and northern China; this split occurred ca. 2.84 million year ago. The majority of populations from the mountains of northern China were dominated by a single haplotype or ribotype, while populations of the HM harbored both high genetic diversity and high haplotype diversity. This distribution pattern indicates that HM was either a diversification center or a refugium for R. kirilowii during the Quaternary climatic oscillations. The present distribution of this species on mountains in northern China may have resulted from a rapid glacial population expansion from the HM. This expansion was confirmed by the mismatch distribution analysis and negative Tajima's D and Fu's FS values, and was dated to ca. 168 thousand years ago. High genetic diversity and population differentiation in both plastid and ITS sequences were revealed; these imply restricted gene flow between populations. A distinct isolation-by-distance pattern was suggested by the Mantel test. Our results show that in old lineages, populations may harbour divergent genetic forms that are sufficient to maintain or even increase overall genetic diversity despite fragmentation and low within-population variation.

Inferred vs realized patterns of gene flow: an analysis of population structure in the Andros Island Rock Iguana

Ecological data, the primary source of information on patterns and rates of migration, can be integrated with genetic data to more accurately describe the realized connectivity between geographically isolated demes. In this paper we implement this approach and discuss its implications for managing populations of the endangered Andros Island Rock Iguana, Cyclura cychlura cychlura. This iguana is endemic to Andros, a highly fragmented landmass of large islands and smaller cays. Field observations suggest that geographically isolated demes were panmictic due to high, inferred rates of gene flow. We expand on these observations using 16 polymorphic microsatellites to investigate the genetic structure and rates of gene flow from 188 Andros Iguanas collected across 23 island sites. Bayesian clustering of specimens assigned individuals to three distinct genotypic clusters. An analysis of molecular variance (AMOVA) indicates that allele frequency differences are responsible for a significant portion of the genetic variance across the three defined clusters (F_st =  0.117, p0.01). These clusters are associated with larger islands and satellite cays isolated by broad water channels with strong currents. These findings imply that broad water channels present greater obstacles to gene flow than was inferred from field observation alone. Additionally, rates of gene flow were indirectly estimated using BAYESASS 3.0. The proportion of individuals originating from within each identified cluster varied from 94.5 to 98.7%, providing further support for local isolation. Our assessment reveals a major disparity between inferred and realized gene flow. We discuss our results in a conservation perspective for species inhabiting highly fragmented landscapes.

Glacial history of the European marine mussels Mytilus, inferred from distribution of mitochondrial DNA lineages

Mussels of the genus Mytilus have been used to assess the circumglacial phylogeography of the intertidal zone. These mussels are representative components of the intertidal zone and have rapidly evolving mitochondrial DNA, suitable for high resolution phylogeographic analyses. In Europe, the three Mytilus species currently share mitochondrial haplotypes, owing to the cases of extensive genetic introgression. Genetic diversity of Mytilus edulis, Mytilus trossulus and Mytilus galloprovincialis was studied using a 900-bp long part of the most variable fragment of the control region from one of their two mitochondrial genomes. To this end, 985 specimens were sampled along the European coasts, at sites ranging from the Black Sea to the White Sea. The relevant DNA fragments were amplified, sequenced and analyzed. Contrary to the earlier findings, our coalescence and nested cladistics results show that only a single M. edulis glacial refugium existed in the Atlantic. Despite that, the species survived the glaciation retaining much of its diversity. Unsurprisingly, M. galloprovincialis survived in the Mediterranean Sea. In a relatively short time period, around the climatic optimum at 10 ky ago, the species underwent rapid expansion coupled with population differentiation. Following the expansion, further contemporary gene flow between populations was limited.

Speciation and introgression between Mimulus nasutus and Mimulus guttatus

Mimulus guttatus and M. nasutus are an evolutionary and ecological model sister species pair differentiated by ecology, mating system, and partial reproductive isolation. Despite extensive research on this system, the history of divergence and differentiation in this sister pair is unclear. We present and analyze a population genomic data set which shows that M. nasutus budded from a central Californian M. guttatus population within the last 200 to 500 thousand years. In this time, the M. nasutus genome has accrued genomic signatures of the transition to predominant selfing, including an elevated proportion of nonsynonymous variants, an accumulation of premature stop codons, and extended levels of linkage disequilibrium. Despite clear biological differentiation, we document genomic signatures of ongoing, bidirectional introgression. We observe a negative relationship between the recombination rate and divergence between M. nasutus and sympatric M. guttatus samples, suggesting that selection acts against M. nasutus ancestry in M. guttatus.

Population genetic structure of a sandstone specialist and a generalist heath species at two levels of sandstone patchiness across the Strait of Gibraltar

Many habitat specialist species are originally composed of small, discontinuous populations because their habitats are naturally fragmented or patchy. They may have suffered the long-term effects of natural patchiness. Mediterranean heathlands, a representative habitat in the Strait of Gibraltar region, are associated with nutrient-poor, acidic sandstone soils. Sandstone soil patches in the African side of the Strait (Tangier) are, in general, smaller and more scattered than in the European side (Algeciras). In this study, we analyze the effect of this sandstone patchiness on the population genetic diversity and structure of two Erica species from these Mediterranean heathlands that differ in their edaphic specificity, E. australis, sandstone specialist, and E. arborea, generalist. Average levels of within-population genetic diversity and gene flow between populations were significantly lower in Tangier (high sandstone patchiness) than in Algeciras (low patchiness) for the sandstone specialist, whereas no differences between both sides of the Strait were detected in the edaphic generalist. Since most endemic species in Mediterranean heathlands of the Strait of Gibraltar are sandstone specialists, these results highlight an increased vulnerability to loss of genetic diversity and local extinction of the heathland endemic flora in the Tangier side of the Strait of Gibraltar.

Admixture and gene flow from Russia in the recovering Northern European brown bear (Ursus arctos)

Large carnivores were persecuted to near extinction during the last centuries, but have now recovered in some countries. It has been proposed earlier that the recovery of the Northern European brown bear is supported by migration from Russia. We tested this hypothesis by obtaining for the first time continuous sampling of the whole Finnish bear population, which is located centrally between the Russian and Scandinavian bear populations. The Finnish population is assumed to experience high gene flow from Russian Karelia. If so, no or a low degree of genetic differentiation between Finnish and Russian bears could be expected. We have genotyped bears extensively from all over Finland using 12 validated microsatellite markers and compared their genetic composition to bears from Russian Karelia, Sweden, and Norway. Our fine masked investigation identified two overlapping genetic clusters structured by isolation-by-distance in Finland (pairwise FST = 0.025). One cluster included Russian bears, and migration analyses showed a high number of migrants from Russia into Finland, providing evidence of eastern gene flow as an important driver during recovery. In comparison, both clusters excluded bears from Sweden and Norway, and we found no migrants from Finland in either country, indicating that eastern gene flow was probably not important for the population recovery in Scandinavia. Our analyses on different spatial scales suggest a continuous bear population in Finland and Russian Karelia, separated from Scandinavia.

The contribution of genetic diversity to subdivide populations living in the silk road of China

There are several indigenous ethnic populations along the silk road in the Northwest of China that display clear differences in culture and social customs, perhaps as a result of geographic isolation and different linguistic traditions. However, extensive trade and other interactions probably facilitated the admixture of different gene pools between these populations over the last two millennia. To further explore the evolutionary relationships of the 13 ethnic populations residing in Northwest China and to reveal the features of population admixture, the 9 most-commonly employed CODIS loci (D3S1358, TH01, D5S818, D13S317, D7S820, CSF1PO, vWA, TPOX, FGA) were selected for genotyping and further analysis. Phylogenetic tree and principal component analysis revealed clear pattern of population differentiation between 4 populations living in Sinkiang Uighur Autonomous Region and other 9 populations dwelled in the upper regions of Silk Road. R matrix regression showed high-level gene flow and population admixture dose exist among these ethic populations in the Northwest region of China. Furthermore, the Mantel test suggests that larger percent of genetic variance (21.58% versus 2.3%) can be explained by geographic isolation than linguistic barriers, which matched with the contribution of geographic factors to other world populations.

Seed dispersal in time can counteract the effect of gene flow between natural populations of Arabidopsis thaliana

Plants may escape unfavorable environments by dispersing to new sites, or by remaining in an ungerminated state at a given site until environmental conditions become favorable. There is limited evidence regarding the occurrence, interplay and relative importance of dispersal processes in time and space in plant populations. Thirty-six natural populations of the annual ruderal species Arabidopsis thaliana were monitored over five consecutive years, sampling both seed bank and above-ground cohorts. We show that immigration rates are considerably higher than previously inferred, averaging 1.7% per population yr(-1). On the other hand, almost one-third of the individuals in a given above-ground cohort result from seeds shed 2 or 3 yr back in time in 10 of the studied populations. Populations that disappeared one year were recolonized by regeneration from the seed bank the subsequent year. Thus, dispersal in both time and space is an important contributor to the structuring of genetic variability in natural populations of A. thaliana, where a high dispersal rate in time may partly counteract the homogenizing effects of spatial seed and pollen dispersal.

Post-glacial expansion and population genetic divergence of mangrove species Avicennia germinans (L.) Stearn and Rhizophora mangle L. along the Mexican coast

Mangrove forests in the Gulf of California, Mexico represent the northernmost populations along the Pacific coast and thus they are likely to be source populations for colonization at higher latitudes as climate becomes more favorable. Today, these populations are relatively small and fragmented and prior research has indicated that they are poor in genetic diversity. Here we set out to investigate whether the low diversity in this region was a result of recent colonization, or fragmentation and genetic drift of once more extensive mangroves due to climatic changes in the recent past. By sampling the two major mangrove species, Rhizophora mangle and Avicennia germinans, along the Pacific and Atlantic coasts of Mexico, we set out to test whether concordant genetic signals could elucidate recent evolution of the ecosystem. Genetic diversity of both mangrove species showed a decreasing trend toward northern latitudes along the Pacific coast. The lowest levels of genetic diversity were found at the range limits around the Gulf of California and the outer Baja California peninsula. Lack of a strong spatial genetic structure in this area and recent northern gene flow in A. germinans suggest recent colonization by this species. On the other hand, lack of a signal of recent northern dispersal in R. mangle, despite the higher dispersal capability of this species, indicates a longer presence of populations, at least in the southern Gulf of California. We suggest that the longer history, together with higher genetic diversity of R. mangle at the range limits, likely provides a gene pool better able to colonize northwards under climate change than A. germinans.

Pre- and post-zygotic reproductive isolation between co-occurring Mussaenda pubescens var. alba and M. shikokiana (Rubiaceae)

Reproductive isolation is a fundamental requirement for speciation and includes several sequential stages. Few studies have determined the relative contributions of pre- and post-zygotic reproductive isolation in plants, especially between relative species with clear differentiation in flower form. To investigate the mechanisms responsible for reproductive isolation in sympatric Mussaenda pubescens var. alba and Mussaenda shikokiana (Rubiaceae) in Guangxi Province, China, we made observations of flowering phenology, patterns of insect visitation, and conducted pollination experiments, including artificial hybridization. The two species had overlapping flowering times and were pollinated by overlapping pollinators; however, their relative importance differed significantly with M. pubescens visited more commonly by bees and M. shikokiana more frequently by butterflies. Using vegetative and floral characters and molecular evidence based on nuclear ribosomal internal and external transcribed spacer regions we detected seven naturally occurring hybrids among a sample of approximately 125 individuals. Hybrids were characterized by morphologies that most closely resembled their maternal parents based on chloroplast evidence. Studies of artificially synthesized and natural hybrids demonstrated that hybrid seed had very low germination rates and naturally occurring hybrids exhibited pollen sterility. Post-zygotic reproductive isolating mechanisms play a primary role in limiting gene exchange between co-occurring species and maintaining species integrity in areas of sympatry.

Adaptive genetic divergence along narrow environmental gradients in four stream insects

A central question linking ecology with evolutionary biology is how environmental heterogeneity can drive adaptive genetic divergence among populations. We examined adaptive divergence of four stream insects from six adjacent catchments in Japan by combining field measures of habitat and resource components with genome scans of non-neutral Amplified Fragment Length Polymorphism (AFLP) loci. Neutral genetic variation was used to measure gene flow and non-neutral genetic variation was used to test for adaptive divergence. We identified the environmental characteristics contributing to divergence by comparing genetic distances at non-neutral loci between sites with Euclidean distances for each of 15 environmental variables. Comparisons were made using partial Mantel tests to control for geographic distance. In all four species, we found strong evidence for non-neutral divergence along environmental gradients at between 6 and 21 loci per species. The relative contribution of these environmental variables to each species' ecological niche was quantified as the specialization index, S, based on ecological data. In each species, the variable most significantly correlated with genetic distance at non-neutral loci was the same variable along which each species was most narrowly distributed (i.e., highest S). These were gradients of elevation (two species), chlorophyll-a, and ammonia-nitrogen. This adaptive divergence occurred in the face of ongoing gene flow (Fst = 0.01-0.04), indicating that selection was strong enough to overcome homogenization at the landscape scale. Our results suggest that adaptive divergence is pronounced, occurs along different environmental gradients for different species, and may consistently occur along the narrowest components of species' niche.

Gene flow and hybridization between numerically imbalanced populations of two duck species on the subantarctic island of South Georgia

Hybridization is common between species of animals, particularly in waterfowl (Anatidae). One factor shown to promote hybridization is restricted mate choice, which can occur when 2 species occur in sympatry but one is rare. According to the Hubbs principle, or "desperation hypothesis," the rarer species is more likely to mate with heterospecifics. We report the second of 2 independent examples of hybridization between 2 species of ducks inhabiting island ecosystems in the Subantarctic and South Atlantic Ocean. Yellow-billed pintails (Anas georgica) and speckled teal (Anas flavirostris) are abundant in continental South America, where they are sympatric and coexist in mixed flocks. But on South Georgia, an isolated island in the Subantarctic, the pintail population of approximately 6000 pairs outnumbers a small breeding population of speckled teal 300∶1. Using 6 genetic loci (mtDNA and 5 nuclear introns) and Bayesian assignment tests coupled with coalescent analyses, we identified hybrid-origin speckled teal alleles in 2 pintails on South Georgia. While it is unclear whether introgression has also occurred into the speckled teal population, our data suggest that this hybridization was not a recent event, but occurred some time ago. We also failed to identify unequivocal evidence of introgression in a much larger sample of pintails and speckled teal from Argentina using a 3-population "Isolation-with-Migration" coalescent analysis. Combined with parallel findings of hybridization between these same 2 duck species in the Falkland Islands, where population ratios are reversed and pintails are outnumbered by speckled teal 1:10, our results provide further support for the desperation hypothesis, which predicts that scarcity in one population and abundance of another will often lead to hybridization. While the South Georgia pintail population appears to be thriving, it's possible that low density of conspecific mates and inverse density dependence (Allee effect) may be one factor limiting the reproductive output of the speckled teal population, and this situation may persist unless speckled teal increase in abundance on South Georgia.

Sequence-based SSR marker development and their application in defining the Introgressions of LA0716 (Solanum pennellii) in the background of cv. M82 (Solanum lycopersicum)

The introgression lines (ILs) from cv. M82 (Solanum lycopersicum) × LA0716 (S. pennellii) have been proven to be exceptionally useful for genetic analysis and gene cloning. The introgressions were originally defined by RFLP markers at their development. The objectives of this study are to develop polymorphic SSR markers, and to re-define the DNA introgression from LA0716 in the ILs. Tomato sequence data was scanned by software to generate SSR markers. In total, 829 SSRs, which could be robustly amplified by PCR, were developed. Among them, 658 SSRs were dinucleotide repeats, 162 were trinucleotide repeats, and nine were tetranucleotide repeats. The 829 SSRs together with 96 published RFLPs were integrated into the physical linkage map of S. lycopersicum. Introgressions of DNA fragments from LA0716 were re-defined among the 75 ILs using the newly developed SSRs. A specific introgression of DNA fragment from LA0716 was identified in 72 ILs as described previously by RFLP, whereas the specific DNA introgression described previously were not detected in the ILs LA4035, LA4059 and LA4091. The physical location of each investigated DNA introgression was finely determined by SSR mapping. Among the 72 ILs, eight ILs showed a shorter and three ILs (IL3-2, IL12-3 and IL12-3-1) revealed a longer DNA introgression than that framed by RFLPs. Furthermore, 54 previously undefined segments were found in 21 ILs, ranging from 1 to 11 DNA introgressions per IL. Generally, the newly developed SSRs provide additional markers for genetic studies of tomatoes, and the fine definition of DNA introgressions from LA0716 would facilitate the use of the ILs for genetic analysis and gene cloning.

Paternal leakage, heteroplasmy, and the evolution of plant mitochondrial genomes

Plant mitochondrial genomes are usually transmitted to the progeny from the maternal parent. However, cases of paternal transmission are known and are perhaps more common than once thought. This review will consider recent evidence, both direct and indirect, of paternal transmission (leakage) of the mitochondrial genome of seed plants, especially in natural populations, and how this can result in offspring that carry a mixture of maternally and paternally derived copies of the genome; a type of heteroplasmy. It will further consider how this heteroplasmy facilitates recombination between genetically distinct partners; a process that can enhance mitochondrial genotypic diversity. This will then form the basis for a discussion of five evolutionary questions that arise from these observations. Questions include how plant mitochondrial genome evolution can be placed on a sexual to asexual continuum, whether cytoplasmic male sterility (CMS) facilitates the evolution of paternal leakage, whether paternal leakage is more likely in populations undergoing admixture, how leakage influences patterns of gene flow, and whether heteroplasmy occurs in natural populations at a frequency greater than predicted by crossing experiments. It is proposed that each of these questions offers fertile ground for future research on a diversity of plant species.

Paternal leakage, heteroplasmy, and the evolution of plant mitochondrial genomes

Plant mitochondrial genomes are usually transmitted to the progeny from the maternal parent. However, cases of paternal transmission are known and are perhaps more common than once thought. This review will consider recent evidence, both direct and indirect, of paternal transmission (leakage) of the mitochondrial genome of seed plants, especially in natural populations, and how this can result in offspring that carry a mixture of maternally and paternally derived copies of the genome; a type of heteroplasmy. It will further consider how this heteroplasmy facilitates recombination between genetically distinct partners; a process that can enhance mitochondrial genotypic diversity. This will then form the basis for a discussion of five evolutionary questions that arise from these observations. Questions include how plant mitochondrial genome evolution can be placed on a sexual to asexual continuum, whether cytoplasmic male sterility (CMS) facilitates the evolution of paternal leakage, whether paternal leakage is more likely in populations undergoing admixture, how leakage influences patterns of gene flow, and whether heteroplasmy occurs in natural populations at a frequency greater than predicted by crossing experiments. It is proposed that each of these questions offers fertile ground for future research on a diversity of plant species.

Complete mitochondrial DNA diversity in Iranians

Due to its pivotal geographical location and proximity to transcontinental migratory routes, Iran has played a key role in subsequent migrations, both prehistoric and historic, between Africa, Asia and Europe. To shed light on the genetic structure of the Iranian population as well as on the expansion patterns and population movements which affected this region, the complete mitochondrial genomes of 352 Iranians were obtained. All Iranian populations studied here exhibit similarly high diversity values comparable to the other groups from the Caucasus, Anatolia and Europe. The results of AMOVA and MDS analyses did not associate any regional and/or linguistic group of populations in the Anatolia/Caucasus and Iran region pointing to close genetic positions of Persians and Qashqais to each other and to Armenians, and Azeris from Iran to Georgians. By reconstructing the complete mtDNA phylogeny of haplogroups R2, N3, U1, U3, U5a1g, U7, H13, HV2, HV12, M5a and C5c we have found a previously unexplored genetic connection between the studied Iranian populations and the Arabian Peninsula, India, Near East and Europe, likely the result of both ancient and recent gene flow. Our results for Persians and Qashqais point to a continuous increase of the population sizes from ∼24 kya to the present, although the phase between 14-24 kya is thought to be hyperarid according to the Gulf Oasis model. Since this would have affected hunter-gatherer ranges and mobility patterns and forced them to increasingly rely on coastal resources, this transition can explain the human expansion across the Persian Gulf region.

Dissimilarity of contemporary and historical gene flow in a wild carrot (Daucus carota) metapopulation under contrasting levels of human disturbance: implications for risk assessment and management of transgene introgression

BACKGROUND AND AIMS

Transgene introgression from crops into wild relatives may increase the resistance of wild plants to herbicides, insects, etc. The chance of transgene introgression depends not only on the rate of hybridization and the establishment of hybrids in local wild populations, but also on the metapopulation dynamics of the wild relative. The aim of the study was to estimate gene flow in a metapopulation for assessing and managing the risks of transgene introgression.

METHODS

Wild carrots (Daucus carota) were sampled from 12 patches in a metapopulation. Eleven microsatellites were used to genotype wild carrots. Genetic structure was estimated based on the FST statistic. Contemporary (over the last several generations) and historical (over many generations) gene flow was estimated with assignment and coalescent methods, respectively.

KEY RESULTS

The genetic structure in the wild carrot metapopulation was moderate (FST = 0·082) and most of the genetic variation resided within patches. A pattern of isolation by distance was detected, suggesting that most of the gene flow occurred between neighbouring patches (≤1 km). The mean contemporary gene flow was 5 times higher than the historical estimate, and the correlation between them was very low. Moreover, the contemporary gene flow in roadsides was twice that in a nature reserve, and the correlation between contemporary and historical estimates was much higher in the nature reserve. Mowing of roadsides may contribute to the increase in contemporary gene flow. Simulations demonstrated that the higher contemporary gene flow could accelerate the process of transgene introgression in the metapopulation.

CONCLUSIONS

Human disturbance such as mowing may alter gene flow patterns in wild populations, affecting the metapopulation dynamics of wild plants and the processes of transgene introgression in the metapopulation. The risk assessment and management of transgene introgression and the control of weeds need to take metapopulation dynamics into consideration.

Genetic diversity and population differentiation of Calanthe tsoongiana, a rare and endemic orchid in China

Calanthe tsoongiana is a rare terrestrial orchid endemic to China, and this species has experienced severe habitat loss and fragmentation. Inter-simple sequence repeat (ISSR) markers were employed to assess the genetic diversity and differentiation of six populations of C. tsoongiana. Based on 124 discernible fragments yielded by eleven selected primers, high genetic diversity was revealed at the species level; however, genetic diversity at the population level was relatively low. High-level genetic differentiation among populations was detected based on analysis of molecular variance (AMOVA), indicating potential limited gene flow. No significant relationship was observed between genetic and geographic distances among the sampled populations. These results suggested that restricted gene flow might be due to habitat fragmentation and reduced population size as a result of human activities. Based on the findings, several conservation strategies were proposed for the preservation of this threatened species.

Genetic structure of Tibeto-Burman populations of Bangladesh: evaluating the gene flow along the sides of Bay-of-Bengal

Human settlement and migrations along sides of Bay-of-Bengal have played a vital role in shaping the genetic landscape of Bangladesh, Eastern India and Southeast Asia. Bangladesh and Northeast India form the vital land bridge between the South and Southeast Asia. To reconstruct the population history of this region and to see whether this diverse region geographically acted as a corridor or barrier for human interaction between South Asia and Southeast Asia, we, for the first time analyzed high resolution uniparental (mtDNA and Y chromosome) and biparental autosomal genetic markers among aboriginal Bangladesh tribes currently speaking Tibeto-Burman language. All the three studied populations; Chakma, Marma and Tripura from Bangladesh showed strikingly high homogeneity among themselves and strong affinities to Northeast Indian Tibeto-Burman groups. However, they show substantially higher molecular diversity than Northeast Indian populations. Unlike Austroasiatic (Munda) speakers of India, we observed equal role of both males and females in shaping the Tibeto-Burman expansion in Southern Asia. Moreover, it is noteworthy that in admixture proportion, TB populations of Bangladesh carry substantially higher mainland Indian ancestry component than Northeast Indian Tibeto-Burmans. Largely similar expansion ages of two major paternal haplogroups (O2a and O3a3c), suggested that they arose before the differentiation of any language group and approximately at the same time. Contrary to the scenario proposed for colonization of Northeast India as male founder effect that occurred within the past 4,000 years, we suggest a significantly deep colonization of this region. Overall, our extensive analysis revealed that the population history of South Asian Tibeto-Burman speakers is more complex than it was suggested before.

Genetic diversity and population structure of rice pathogen Ustilaginoidea virens in China

Rice false smut caused by the fungal pathogen Ustilaginoidea virens is becoming a destructive disease throughout major rice-growing countries. Information about its genetic diversity and population structure is essential for rice breeding and efficient control of the disease. This study compared the genome sequences of two U. virens isolates. Three SNP-rich genomic regions were identified as molecular markers that could be used to analyze the genetic diversity and population structure of U. virens in China. A total of 56 multilocus sequence types (haplotypes) were identified out of 162 representative isolates from 15 provinces covering five major rice-growing areas in China. However, the phylogeny, based on sequences at individual SNP-rich regions, strongly conflicted with each other and there were significant genetic differences between different geographical populations. Gene flow between the different geographical populations and genetic differentiation within each geographical population were also detected. In addition, genetic recombination and genetic isolation resulting from geographic separation was also found.

Structure, gene flow, and recombination among geographic populations of a Russula virescens ally from Southwestern China

Mushrooms that are morphologically indistinguishable from Russula virescens (Schaeff.) Fr. are among the most popular wild edible mushrooms in Yunnan province, southwestern China. However, almost nothing is known about their biology. This study investigated the diversity and population genetics of a R. virescens ally from Yunnan. A total of 210 samples were collected from 13 geographical locations throughout the main distribution range in Yunnan. The patterns of genetic variation within and among these geographic populations were analyzed using sequences from three nuclear and two mitochondrial DNA fragments. Analysis of the ITS sequences revealed that samples from Yunnan showed 3–6% sequence difference from R. virescens samples from North America and Europe and formed a distinct clade. Our multilocus population genetic analyses suggested frequent gene flow among geographic populations of the R. virescens ally in Yunnan. Interestingly, the nuclear and mitochondrial genes exhibited different levels of gene flow and recombination. We discuss the implications of our results for understanding speciation, reproduction and conservation of this important biological resource.

Assessment of genetic diversity in the sorghum reference set using EST-SSR markers

Selection and use of genetically diverse genotypes are key factors in any crop breeding program to develop cultivars with a broad genetic base. Molecular markers play a major role in selecting diverse genotypes. In the present study, a reference set representing a wide range of sorghum genetic diversity was screened with 40 EST-SSR markers to validate both the use of these markers for genetic structure analyses and the population structure of this set. Grouping of accessions is identical in distance-based and model-based clustering methods. Genotypes were grouped primarily based on race within the geographic origins. Accessions derived from the African continent contributed 88.6 % of alleles confirming the African origin of sorghum. In total, 360 alleles were detected in the reference set with an average of 9 alleles per marker. The average PIC value was 0.5230 with a range of 0.1379-0.9483. Sub-race, guinea margaritiferum (Gma) from West Africa formed a separate cluster in close proximity to wild accessions suggesting that the Gma group represents an independent domestication event. Guineas from India and Western Africa formed two distinct clusters. Accessions belongs to the kafir race formed the most homogeneous group as observed in earlier studies. This analysis suggests that the EST-SSR markers used in the present study have greater discriminating power than the genomic SSRs. Genetic variance within the subpopulations was very high (71.7 %) suggesting that the germplasm lines included in the set are more diverse. Thus, this reference set representing the global germplasm is an ideal material for the breeding community, serving as a community resource for trait-specific allele mining as well as genome-wide association mapping.

Genetic structure of Lycorma delicatula (Hemiptera: Fulgoridae) populations in Korea: implication for invasion processes in heterogeneous landscapes

Lycorma delicatula (White) was identified in 2004 as an invasive pest in South Korea, where it causes serious damage to vineyard crops. To investigate the population structure and dispersal pattern of L. delicatula in South Korea, we estimated the population genetic structure and gene flow among nine locations across the country using seven microsatellite markers. Although L. delicatula spread throughout most of its geographical range in South Korea within 5-7 years following invasion, its populations show evidence of genetic structuring across the range with a low but significant global F ST (genetic differentiation across all populations) of 0.0474. Bayesian-based clustering analysis indicates the presence of at least three genetically unique populations in South Korea, including populations in northeastern South Korea, which show a distinct genetic background. However, isolation by distance suggests that populations in South Korea have not yet reached genetic equilibrium. Estimates of the historical rate of gene flow (N e m) indicate that relatively high rates of flow have been maintained among populations within the western region, which may indicate recent range expansion. A population assignment test using the first-generation migrant detection method suggested that long-distance dispersal of L. delicatula may have occurred over large areas of South Korea. More complex dispersal patterns may have occurred during L. delicatula invasion of heterogeneous landscapes in South Korea.

Mitochondrial cox sequences of Nilaparvata lugens and Sogatella furcifera (Hemiptera, Delphacidae): low specificity among Asian planthopper populations

The brown planthoppers (BPH) Nilaparvata lugens (Stål) and the white-backed planthoppers (WBPH) Sogatella furcifera (Horváth) annually migrate from tropical and subtropical regions to temperate regions in Asia, including Japan, Korea and northern China. To elucidate the genetic divergence based on geography of planthoppers and to estimate their migration route on the basis of molecular data, we analysed a part of their mitochondrial genome sequences. Sequences of cytochrome oxidase subunit I (cox1) - transfer RNA for Leu (trnL2) - cox2 were determined for 579 BPH (1,928 bp) and 464 WBPH (1,927 bp) individuals collected from 31 and 25 locations, respectively, in East and Southeast Asia. Thirty and 20 mitochondrial haplotypes were detected for BPH and WBPH, respectively. Single populations of both planthoppers included multiple haplotypes, and many haplotypes were shared in some populations and areas. The most frequently detected haplotypes accounted for approximately 50% of all BPH and WBPH individuals. To evaluate gene flow among planthoppers in different regions in Asia, pairwise fixation index (Fst) values were calculated. For BPH, high Fst values (0.580-0.926) were shown between planthoppers in Papua New Guinea (PNG) and the other areas and moderate Fst values (0.176-0.362) were observed between those in southern Philippines and other areas. For WBPH, the Fst value was the highest between Taiwan and southern Vietnam (0.236), and low among the other areas. AMOVA indicated no genetic structure among eight areas, excluding southern Philippines and PNG, for BPH, and among ten areas for WBPH. These data indicate that both planthoppers do not show much differentiation of local populations and/or have genetically intermixed Asian populations. These data also indicate that it may be difficult to distinguish regional planthopper populations on the basis of differences in mitochondrial sequences.

Gene flow vs. pollution pressure: genetic diversity of Mytilus galloprovincialis in eastern Adriatic

Environmental pollution may modify all the evolutionary processes involved in shaping the genetic patterns of exposed populations. In order to evaluate the pollution impact on the genetic diversity of Mediterranean mussel Mytilus galloprovincialis ten populations inhabiting differently polluted sites along the eastern Adriatic coast, from pristine bays to heavily trafficked harbours, were studied. Pollution pressure was assessed through an integrated study of biological effects and responses across different levels of biological organization. Eight microsatellite markers were analysed to assess genetic diversity of investigated populations. Both the principal component analysis (PCA) of the biomarker data set as well as the biomarker response index (BRI) confirmed substantial pollution pressure at the highly polluted sites, and very low pollution exposure at the three reference sites. Very shallow genetic differentiation was found in respect to maritime distances or pollution status, and this was attributed to a high gene flow among the populations. However, populations inhabiting polluted sites exhibited higher levels of genetic diversity and evolutionary mechanisms underlying this phenomenon are discussed.

When the rule becomes the exception. no evidence of gene flow between two Zerynthia cryptic butterflies suggests the emergence of a new model group

There is increasing evidence that most parapatric cryptic/sister taxa are reproductively compatible across their areas of contact. Consequently, the biological species concept, which assumes absence of interbreeding, is becoming a not so effective criterion in evolutionary ecology. Nevertheless, the few parapatric sister taxa showing complete reproductive barriers represent interesting models to study speciation processes and the evolution of reproductive isolation. In this study, we examined contact populations in northwestern Italy of two butterfly species, Zerynthia polyxena and Z. cassandra, characterized by different genitalic morphotypes. We studied levels of divergence among 21 populations distributed from Sicily to France using three genetic markers (the mitochondrial COI and ND1 genes and the nuclear wingless gene) and genitalic geometric morphometrics. Moreover, we performed species distribution modelling to estimate different climatic requirements of Z. polyxena and Z. cassandra. We projected climatic data into glacial maximum scenarios in order to verify if and to which extent glacial cycles could have contributed to speciation processes. Genetic and morphometric analyses identified two main groups. All specimens showed a concordant pattern of diversification, including those individuals sampled in the contact area. Haplotype distribution and climatic models showed that during glacial maxima both species experienced a strong range contraction and presumably remained separated into different microrefugia in southern France, in the Italian Peninsula and on the islands of Elba and Sicily. Long term separation was probably favoured by reduced dispersal ability and high phylopatry, while genitalic diversification probably favoured interbreeding avoidance. Conversely, the aposematic wing pattern remained almost identical. We compared our results with those obtained in other species and concluded that Z. polyxena and Z. cassandra represent a valuable model in the study of speciation.