Patterns Of Nucleotide Diversity Research Articles

Evolutionary patterns of plastome uncover diploid-polyploid maternal relationships in Triticeae

To investigate the diploid-polyploid relationships and the role of maternal progenitors in establishment of polyploid richness in Triticeae, 35 polyploids representing almost all genomic constitutions together with 48 diploid taxa representing 20 basic genomes in the tribe were analyzed. Phylogenomic reconstruction, genetic distance matrix, and nucleotide diversity patterns of plastome sequences indicated that (1) The maternal donor of the annual polyploid species with the U- and D-genome are related to extant Ae. umbellulata and Ae. tauschii, respectively. The maternal donor to the annual polyploid species with the S-, G-, and B-genome originated from the species of Sitopsis section of the genus Aegilops. The annual species with the Xe-containing polyploids were donated by Eremopyrum as the female parent; (2) Pseudoroegneria and Psathyrostachys were the maternal donor of perennial species with the St- and Ns-containing polyploids, respectively; (3) The Lophopyrum, Thinopyrum and Dasypyrum genomes contributed cytoplasm genome to Pseudoroegneria species as a result of incomplete lineage sorting and/or chloroplast captures, and these lineages were genetically transmitted to the St-containing polyploid species via polyploidization; (4) There is a reticulate relationship among the St-containing polyploid species. It can be suggested that genetic heterogeneity might associate with the richness of the polyploids in Triticeae.

Processing and Analyzing Multiple Genomes Alignments with MafFilter.

As the number of available genome sequences from both closely related species and individuals within species increased, theoretical and methodological convergences between the fields of phylogenomics and population genomics emerged. Population genomics typically focuses on the analysis of variants, while phylogenomics heavily relies on genome alignments. However, these are playing an increasingly important role in studies at the population level. Multiple genome alignments of individuals are used when structural variation is of primary interest and when genome architecture permits to assemble de novo genome sequences. Here I describe MafFilter, a command-line-driven program allowing to process genome alignments in the Multiple Alignment Format (MAF). Using concrete examples based on publicly available datasets, I demonstrate how MafFilter can be used to develop efficient and reproducible pipelines with quality assurance for downstream analyses. I further show how MafFilter can be used to perform both basic and advanced population genomic analyses in order to infer the patterns of nucleotide diversity along genomes.

Genetic variation in the Chinese soft-shell turtles (Pelodiscus spp.) revealed by sequences of mitochondrial Cytb gene

Genetic diversity is an important component of biodiversity and investigating and protecting the wild genetic diversity of species has always been an important research topic in conservation biology. Due to habitat destruction and over catching, wild Chinese soft-shell turtles (Pelodiscus spp.) have been severely damaged, resulting in the species being listed on the IUCN Red List of Threatened Species in 2000. However, only few studies with contradictory results on the genetic diversity of Pelodiscus turtles have been reported. To investigate the genetic diversity and population structure of Pelodiscus turtles, 123 specimens were collected from five localities in East Asia, and their genetic variation was analyzed on the basis of a 922-bp partial sequence of the mitochondrial cytb gene. Forty-nine polymorphic sites were detected, revealing 56 haplotypes. A pattern of high haplotype diversity (Hd = 0.994) and nucleotide diversity (π = 0.01655) was found in the examined range. A maximum likelihood phylogenetic tree and median-joining network analyses indicated that the Pelodiscus turtles divided in four populations throughout East Asia, and the Pelodiscus turtles probably originate from the Yangtze River and was introduced to Yellow River, Taiwan, and Japan through Xijiang River. These results served a helpful resource for conservation of Pelodiscus turtles.

Heterogeneous patterns of genetic variation at nuclear genes and quantitative traits in a Scots pine provenance trial

We studied the genetic variation in a set of nuclear genes analyzed from 16 populations of Scots pine derived from a 50-year-old provenance trial in Poland. At the same set of loci, the pattern of genetic variation was compared to several reference populations from a latitudinal gradient in Northern and Central Europe. Similar levels of nucleotide diversity were observed between the defined groups of Polish populations representing three climatic regions (π<sub>total</sub> = 0.0040–0.0051) in comparison with the reference samples (π<sub>total</sub> = 0.0054–0.0058). Polish populations showed minor but heterogeneous patterns of genetic variation between regional groups (<em>F</em><sub>ST</sub> up to 6%), which were caused by differentiation at specific loci. When outlier loci were excluded from between group comparisons, there were no differences between the Polish populations. Loci related to glycosyltransferase and laccase were identified as outliers, and were correlated with phenotypic differentiation using mixed-linear models. Moreover, these genes were also found as being potentially under selection across the Scots pine distribution range as the patterns of nucleotide variation correlated with latitude and altitude of the maternal stands. The provenance trial measurements have characterized a set of growth and developmental traits over 50 years and forms a suitable experimental system for detailed genetic studies.

Nucleotide diversity patterns at the DREB1 transcriptional factor gene in the genome donor species of wheat (Triticum aestivum L).

Bread wheat (AABBDD) originated from the diploid progenitor Triticum urartu (AA), a relative of Aegilops speltoides (BB), and Ae. tauschii (DD). The DREB1 transcriptional factor plays key regulatory role in low-temperature tolerance. The modern breeding strategies resulted in serious decrease of the agricultural biodiversity, which led to a loss of elite genes underlying abiotic stress tolerance in crops. However, knowledge of this gene’s natural diversity is largely unknown in the genome donor species of wheat. We characterized the dehydration response element binding protein 1 (DREB1) gene-diversity pattern in Ae. speltoides, Ae. tauschii, T. monococcum and T. urartu. The highest nucleotide diversity value was detected in Ae. speltoides, followed by Ae. tauschii and T. monococcum. The lowest nucleotide diversity value was observed in T. urartu. Nucleotide diversity and haplotype data might suggest no reduction of nucleotide diversity during T. monococcum domestication. Alignment of the 68 DREB1 sequences found a large-size (70 bp) insertion/deletion in the accession PI486264 of Ae. speltoides, which was different from the copy of sequences from other accessions of Ae. speltoides, suggesting a likely existence of two different ancestral Ae. speltoides forms. Implication of sequences variation of Ae. speltoides on origination of B genome in wheat was discussed.

Identification and analysis of tRNA genes provide new insights into oil biosynthesis in tung tree (Vernicia fordii Hemsl.)

Tung tree (Vernicia fordii Hemsl.) is a potential commercial source of biodiesel which is widely planted in China and many other countries. Investigation of the regulatory mechanism of lipid biosynthesis is an essential prerequisite to improve tung tree varieties by genetic engineering. In this study, the transfer RNA (tRNA) genes were identified in four genome-sequenced species (tung tree, rubber tree, castor bean and physic nut) in Euphorbiaceae. Gene clusters were observed for the four species tRNAs. Though tRNAs were highly conserved, contrasting patterns of nucleotide diversity were observed among the tRNA family. A total of 86 differentially expressed tRNA-derived RNA fragments (tRFs) were identified of which the target genes were enriched in lipid biosynthesis by Gene Ontology (GO) enrichment analysis. Interestingly, 11 lipid genes previously reported were identified and found to be regulated by 17 tRFs. A tRF-mRNA regulatory network was constructed, which provides a new insight into the genetic mechanism of lipid biosynthesis in tung tree.

Understanding the Factors That Shape Patterns of Nucleotide Diversity in the House Mouse Genome.

A major goal of population genetics has been to determine the extent by which selection at linked sites influences patterns of neutral nucleotide diversity in the genome. Multiple lines of evidence suggest that diversity is influenced by both positive and negative selection. For example, in many species there are troughs in diversity surrounding functional genomic elements, consistent with the action of either background selection (BGS) or selective sweeps. In this study, we investigated the causes of the diversity troughs that are observed in the wild house mouse genome. Using the unfolded site frequency spectrum, we estimated the strength and frequencies of deleterious and advantageous mutations occurring in different functional elements in the genome. We then used these estimates to parameterize forward-in-time simulations of chromosomes, using realistic distributions of functional elements and recombination rate variation in order to determine whether selection at linked sites can explain the observed patterns of nucleotide diversity. The simulations suggest that BGS alone cannot explain the dips in diversity around either exons or conserved noncoding elements. A combination of BGS and selective sweeps produces deeper dips in diversity than BGS alone, but the inferred parameters of selection cannot fully explain the patterns observed in the genome. Our results provide evidence of sweeps shaping patterns of nucleotide diversity across the mouse genome and also suggest that infrequent, strongly advantageous mutations play an important role in this. The limitations of using the unfolded site frequency spectrum for inferring the frequency and effects of advantageous mutations are discussed.

Effects of Geological and Environmental Events on the Diversity and Genetic Divergence of Four Closely Related Pines: Pinus koraiensis, P. armandii, P. griffithii, and P. pumila.

The effects of mountain uplift and environmental oscillations on nucleotide variability and species divergence remain largely unknown in East Asia. In this study, based on multiple nuclear DNA markers, we investigated the levels and patterns of nucleotide diversity and interspecific divergence in four closely related pines in China, i.e., Pinus koraiensis, P. armandii, P. griffithii, and P. pumila. The four pine taxa shared low levels of nucleotide polymorphisms at the species level. P. pumila had the highest silent nucleotide diversity (πsil = 0.00661) whereas P. griffithii had the lowest (πsil = 0.00175), while the levels of genetic polymorphism in P. armandii (πsil = 0.00508) and P. koraiensis (πsil = 0.00652) were intermediate between the other two species. Population genetic structure analysis showed that variations primarily existed within populations of the four pine species, presumably due to habitat fragmentation or the island-like distributions of Pinus species. Population divergence (FST) analysis showed that the genetic divergence between P. griffithii and P. koraiensis was much greater than that between P. koraiensis and the other two pines species. Isolation-with-migration analysis suggested that asymmetric gene flow had occurred between any two pairs of pine species. Phylogenetic analyses indicated that the four allied species split into two groups about 1.37 million years ago, where P. armandii and P. pumila were closer and clustered as sister species, whereas P. koraiensis and P. griffithii were clustered on another branch. Our results and those obtained in previous studies suggest that mountain uplift and geological climate oscillations may have led to the patterns of genetic divergence and nucleotide variations in these four pine species.

Modeling SNP array ascertainment with Approximate Bayesian Computation for demographic inference

Single nucleotide polymorphisms (SNPs) in commercial arrays have often been discovered in a small number of samples from selected populations. This ascertainment skews patterns of nucleotide diversity and affects population genetic inferences. We propose a demographic inference pipeline that explicitly models the SNP discovery protocol in an Approximate Bayesian Computation (ABC) framework. We simulated genomic regions according to a demographic model incorporating parameters for the divergence of three well-characterized HapMap populations and recreated the SNP distribution of a commercial array by varying the number of haploid samples and the allele frequency cut-off in the given regions. We then calculated summary statistics obtained from both the ascertained and genomic data and inferred ascertainment and demographic parameters. We implemented our pipeline to study the admixture process that gave rise to the present-day Mexican population. Our estimate of the time of admixture is closer to the historical dates than those in previous works which did not consider ascertainment bias. Although the use of whole genome sequences for demographic inference is becoming the norm, there are still underrepresented areas of the world from where only SNP array data are available. Our inference framework is applicable to those cases and will help with the demographic inference.

Comparative Genomic Analysis of 130 Bacteriophages Infecting Bacteria in the Genus Pseudomonas

Bacteria of the genus Pseudomonas are genetically diverse and ubiquitous in the environment. Like other bacteria, those of the genus Pseudomonas are susceptible to bacteriophages which can significantly affect their host in many ways, ranging from cell lysis to major changes in morphology and virulence. Insights into phage genomes, evolution, and functional relationships with their hosts have the potential to contribute to a broader understanding of Pseudomonas biology, and the development of novel phage therapy strategies. Here we provide a broad-based comparative and evolutionary analysis of 130 complete Pseudomonas phage genome sequences available in online databases. We discovered extensive variation in genome size (ranging from 3 to 316 kb), G + C percentage (ranging from 37 to 66%), and overall gene content (ranging from 81–96% of genome space). Based on overall nucleotide similarity and the numbers of shared gene products, 100 out of 130 genome sequences were grouped into 12 different clusters; 30 were characterized as singletons, which do not have close relationships with other phage genomes. For 5/12 clusters, constituent phage members originated from two or more different Pseudomonas host species, suggesting that phage in these clusters can traverse bacterial species boundaries. An analysis of CRISPR spacers in Pseudomonas bacterial genome sequences supported this finding. Substantial diversity was revealed in analyses of phage gene families; out of 4,462 total families, the largest had only 39 members and there were 2,992 families with only one member. An evolutionary analysis of 72 phage gene families, based on patterns of nucleotide diversity at non-synonymous and synonymous sites, revealed strong and consistent signals for purifying selection. Our study revealed highly diverse and dynamic Pseudomonas phage genomes, and evidence for a dominant role of purifying selection in shaping the evolution of genes encoded in them.

Structural patterns of selection and diversity for Plasmodium vivax antigens DBP and AMA1

BackgroundPlasmodium vivax is a significant contributor to the global malaria burden, and a vaccine targeting vivax malaria is urgently needed. An understanding of the targets of functional immune responses during the course of natural infection will aid in the development of a vaccine. Antibodies play a key role in this process, with responses against particular epitopes leading to immune selection pressure on these epitopes. A number of techniques exist to estimate levels of immune selection pressure on particular epitopes, with a sliding window analysis often used to determine particular regions likely to be under immune pressure. However, such analysis neglects protein three-dimensional structural information. With this in mind, a newly developed tool, BioStructMap, was applied to two key antigens from Plasmodium vivax: PvAMA1 and PvDBP Region II. This tool incorporates structural information into tests of selection pressure.ResultsSequences from a number of populations were analysed, examining spatially-derived nucleotide diversity and Tajima’s D over protein structures for PvAMA1 and PvDBP. Structural patterns of nucleotide diversity were similar across all populations examined, with Domain I of PvAMA1 having the highest nucleotide diversity and displaying significant signatures of immune selection pressure (Tajima’s D > 0). Nucleotide diversity for PvDBP was highest bordering the dimerization and DARC-binding interface, although there was less evidence of immune selection pressure on PvDBP compared with PvAMA1. This study supports previous work that has identified Domain I as the main target of immune-mediated selection pressure for PvAMA1, and also supports studies that have identified functional epitopes within PvDBP Region II.ConclusionsThe BioStructMap tool was applied to leading vaccine candidates from P. vivax, to examine structural patterns of selection and diversity across a number of geographic populations. There were striking similarities in structural patterns of diversity across multiple populations. Furthermore, whilst regions of high diversity tended to surround conserved binding interfaces, a number of protein regions with very low diversity were also identified, and these may be useful targets for further vaccine development, given previous evidence of functional antibody responses against these regions.

Population genetic structure and natural selection of Plasmodium falciparum apical membrane antigen-1 in Myanmar isolates

BackgroundPlasmodium falciparum apical membrane antigen-1 (PfAMA-1) is one of leading blood stage malaria vaccine candidates. However, genetic variation and antigenic diversity identified in global PfAMA-1 are major hurdles in the development of an effective vaccine based on this antigen. In this study, genetic structure and the effect of natural selection of PfAMA-1 among Myanmar P. falciparum isolates were analysed.MethodsBlood samples were collected from 58 Myanmar patients with falciparum malaria. Full-length PfAMA-1 gene was amplified by polymerase chain reaction and cloned into a TA cloning vector. PfAMA-1 sequence of each isolate was sequenced. Polymorphic characteristics and effect of natural selection were analysed with using DNASTAR, MEGA4, and DnaSP programs. Polymorphic nature and natural selection in 459 global PfAMA-1 were also analysed.ResultsThirty-seven different haplotypes of PfAMA-1 were identified in 58 Myanmar P. falciparum isolates. Most amino acid changes identified in Myanmar PfAMA-1 were found in domains I and III. Overall patterns of amino acid changes in Myanmar PfAMA-1 were similar to those in global PfAMA-1. However, frequencies of amino acid changes differed by country. Novel amino acid changes in Myanmar PfAMA-1 were also identified. Evidences for natural selection and recombination event were observed in global PfAMA-1. Among 51 commonly identified amino acid changes in global PfAMA-1 sequences, 43 were found in predicted RBC-binding sites, B-cell epitopes, or IUR regions.ConclusionsMyanmar PfAMA-1 showed similar patterns of nucleotide diversity and amino acid polymorphisms compared to those of global PfAMA-1. Balancing natural selection and intragenic recombination across PfAMA-1 are likely to play major roles in generating genetic diversity in global PfAMA-1. Most common amino acid changes in global PfAMA-1 were located in predicted B-cell epitopes where high levels of nucleotide diversity and balancing natural selection were found. These results highlight the strong selective pressure of host immunity on the PfAMA-1 gene. These results have significant implications in understanding the nature of Myanmar PfAMA-1 along with global PfAMA-1. They also provide useful information for the development of effective malaria vaccine based on this antigen.

Patterns of Genome-Wide Nucleotide Diversity in the Gynodioecious Plant Thymus vulgaris Are Compatible with Recent Sweeps of Cytoplasmic Genes.

Gynodioecy is a sexual dimorphism where females coexist with hermaphrodite individuals. In most cases, this dimorphism involves the interaction of cytoplasmic male sterility (CMS) genes and nuclear restorer genes. Two scenarios can account for how these interactions maintain gynodioecy. Either CMS genes recurrently enter populations at low frequency via mutation or migration and go to fixation unimpeded (successive sweeps), or CMS genes maintain polymorphism over evolutionary time through interactions with a nuclear restorer allele (balanced polymorphism). To distinguish between these scenarios, we used transcriptome sequencing in gynodioecious Thymus vulgaris and surveyed genome-wide diversity in 18 naturally occurring individuals sampled from populations at a local geographic scale. We contrast the amount and patterns of nucleotide diversity in the nuclear and cytoplasmic genome, and find ample diversity at the nuclear level (π = 0.019 at synonymous sites) but reduced genetic diversity and an excess of rare polymorphisms in the cytoplasmic genome relative to the nuclear genome. Our finding is incompatible with the maintenance of gynodioecy via scenarios invoking long-term balancing selection, and instead suggests the recent fixation of CMS lineages in the populations studied.

High rate of adaptive evolution in two widespread European pines.

Comparing related organisms with differing ecological requirements and evolutionary histories can shed light on the mechanisms and drivers underlying genetic adaptation. Here, by examining a common set of hundreds of loci, we compare patterns of nucleotide diversity and molecular adaptation of two European conifers (Scots pine and maritime pine) living in contrasted environments and characterized by distinct population genetic structure (low and clinal in Scots pine, high and ecotypic in maritime pine) and demographic histories. We found higher nucleotide diversity in Scots pine than in maritime pine, whereas rates of new adaptive substitutions (ωa ), as estimated from the distribution of fitness effects, were similar across species and among the highest found in plants. Sample size and population genetic structure did not appear to have resulted in significant bias in estimates of ωa . Moreover, population contraction-expansion dynamics for each species did not affect differentially the rate of adaptive substitution in these two pines. Several methodological and biological factors may underlie the unusually high rate of adaptive evolution of Scots pine and maritime pine. By providing two new case studies with contrasting evolutionary histories, we contribute to disentangling the multiple factors potentially affecting adaptive evolution in natural plant populations.

Single-Molecule Sequencing Reveals the Chromosome-Scale Genomic Architecture of the Nematode Model Organism Pristionchus pacificus

The nematode Pristionchus pacificus is an established model for integrative evolutionary biology and comparative studies with Caenorhabditis elegans. While an existing genome draft facilitated the identification of several genes controlling various developmental processes, its high degree of fragmentation complicated virtually all genomic analyses. Here, we present a de novo genome assembly from single-molecule, long-read sequencing data consisting of 135 P.pacificus contigs. When combined with a genetic linkage map, 99% of the assembly could be ordered and oriented into six chromosomes. This allowed us to robustly characterize chromosomal patterns of gene density, repeat content, nucleotide diversity, linkage disequilibrium, and macrosynteny in P.pacificus. Despite widespread conservation of synteny between P.pacificus and C.elegans, we identified one major translocation from an autosome to the sex chromosome in the lineage leading to C.elegans. This highlights the potential of the chromosome-scale assembly for future genomic studies of P.pacificus.

Signatures of soft sweeps across the Dt1 locus underlying determinate growth habit in soya bean [Glycine max (L.) Merr.

Determinate growth habit is an agronomically important trait associated with domestication in soya bean. Previous studies have demonstrated that the emergence of determinacy is correlated with artificial selection on four nonsynonymous mutations in the Dt1 gene. To better understand the signatures of the soft sweeps across the Dt1 locus and track the origins of the determinate alleles, we examined patterns of nucleotide variation in Dt1 and the surrounding genomic region of approximately 800kb. Four local, asymmetrical hard sweeps on four determinate alleles, sized approximately 660, 120, 220 and 150kb, were identified, which constitute the soft sweeps for the adaptation. These variable-sized sweeps substantially reflected the strength and timing of selection and indicated that the selection on the alleles had been completed rapidly within half a century. Statistics of EHH, iHS, H12 and H2/H1 based on haplotype data had the power to detect the soft sweeps, revealing distinct signatures of extensive long-range LD and haplotype homozygosity, and multiple frequent adaptive haplotypes. A haplotype network constructed for Dt1 and a phylogenetic tree based on its extended haplotype block implied independent sources of the adaptive alleles through de novo mutations or rare standing variation in quick succession during the selective phase, strongly supporting multiple origins of the determinacy. We propose that the adaptation of soya bean determinacy is guided by a model of soft sweeps and that this model might be indispensable during crop domestication or evolution.

No obvious genetic erosion, but evident relict status at the westernmost range edge of the Pontic-Pannonian steppe plant Linum flavum L. (Linaceae) in Central Europe.

We investigate patterns of genetic variation along an east–west transect of Central European populations of Linum flavum and interpret the Quaternary history of its peripheral populations, especially those at the westernmost isolated range edge, discussing their migrations and possible relict status. We defined our peripheral transect across three study regions from Central Hungary, eastern Austria to southwestern Germany. Using AFLP fingerprinting and cpDNA sequence variation (rpL16 intron, atpI‐H), we analyzed 267 and 95 individuals, respectively, representing each study region by four populations. Hierarchical AMOVA (AFLPs) indicated significant variation among study regions (12% of total variance) and moderate differentiation between populations (10%). Population differentiation was high at the westernmost range edge (11.5%, Germany), but also in the east (13.4%, Hungary), compared to the Austrian study region (8.6%). Correspondingly, AFLP diversity was highest in the center of the study transect in eastern Austria. CpDNA haplotypes support a pattern of regional structuring with the strongest separation of the westernmost range edge, and some haplotype sharing among Austrian and Hungarian individuals. Equilibrating nucleotide versus haplotype diversity patterns, the highly diverse populations at the Pannonian range edge (Austria) indicate long‐term persistence, while Central Pannonian populations are obviously effected by recent bottlenecks. Intermediate nucleotide, but high haplotype diversity within the westernmost exclave (Swabian Alb), is indicative of a founder bottleneck during its pre‐LGM or early postglacial migration history, followed by sufficient time to accumulate cpDNA variation. The not obviously reduced genetic diversity and distinctiveness of L. flavum at the westernmost range edge suggest a long‐term persistence (relict status) of populations in this region, where the species has survived probably even the Würm glaciation in extra‐Mediterranean refugia. This genetic relict variation represents an important part of the overall genetic diversity found in the western periphery of this steppe plant and highlights the high conservation priority of respective gene pools.

Domestication rewired gene expression and nucleotide diversity patterns in tomato

Plant domestication has led to considerable phenotypic modifications from wild species to modern varieties. However, although changes in key traits have been well documented, less is known about the underlying molecular mechanisms, such as the reduction of molecular diversity or global gene co-expression patterns. In this study, we used a combination of gene expression and population genetics in wild and crop tomato to decipher the footprints of domestication. We found a set of 1729 differentially expressed genes (DEG) between the two genetic groups, belonging to 17 clusters of co-expressed DEG, suggesting that domestication affected not only individual genes but also regulatory networks. Five co-expression clusters were enriched in functional terms involving carbohydrate metabolism or epigenetic regulation of gene expression. We detected differences in nucleotide diversity between the crop and wild groups specific to DEG. Our study provides an extensive profiling of the rewiring of gene co-expression induced by the domestication syndrome in one of the main crop species.

Genomic islands of divergence or opportunities for introgression?

In animals, introgression between species is often perceived as the breakdown of reproductive isolating mechanisms, but gene flow between incipient species can also represent a source for potentially beneficial alleles. Recently, genome-wide datasets have revealed clusters of differentiated loci ('genomic islands of divergence') that are thought to play a role in reproductive isolation and therefore have reduced gene flow. We use simulations to further examine the evolutionary forces that shape and maintain genomic islands of divergence between two subspecies of the migratory songbird, Swainson's thrush (Catharus ustulatus), which have come into secondary contact since the last glacial maximum. We find that, contrary to expectation, gene flow is high within islands and is highly asymmetric. In addition, patterns of nucleotide diversity at highly differentiated loci suggest selection was more frequent in a single ecotype. We propose a mechanism whereby beneficial alleles spread via selective sweeps following a post-glacial demographic expansion in one subspecies and move preferentially across the hybrid zone. We find no evidence that genomic islands are the result of divergent selection or reproductive isolation, rather our results suggest that differentiated loci both within and outside islands could provide opportunities for adaptive introgression across porous species boundaries.

Molecular evolution of glutathione peroxidase (Gpx) family genes in desert poplar (Populus euphratica Oliv.)

Populus euphratica Oliv. is a poplar species that is distributed mainly in deserts, making it an interesting model in which to investigate molecular mechanisms underlying different stress responses. Here, we used molecular population genetic methods to detect potential selection in candidate genes belonging to the P. euphratica glutathione (GSH) peroxidase (Gpx) family, which are associated with an enzymatic mechanism that combats oxidative damage caused by reactive oxygen species (ROS) in plant cells; earlier studies have shown that Gpx proteins play important roles in coping with increased ROS levels during biotic and abiotic stresses in plants. We analyzed the nucleotide diversity and divergence patterns of five loci encoding Gpx genes, and 16 reference loci used as controls, in order to detect departures from the neutral expectation. Gpx1 has an excess of mid-frequency alleles; high intraspecific nucleotide diversity, distributed in the upper tail of the simulated neutral model; and extensive LD, showing strong evidence of balancing selection/local adaptation. The Gpx3.2 gene exhibits very low nucleotide diversity and divergence, suggesting that it has evolved under strong purifying selection. We failed to detect any evidence for natural selection at the other loci (Gpx2, Gpx4, and Gpx5) compared with the reference loci. The results show that nucleotide diversity and/or divergence differ greatly between members of the Gpx gene family, resulting from differential selective pressure acting on genes, and that adaptive evolution influenced the distribution of P. euphratica in desert regions.