Behavioural vs. molecular sources of conflict between nuclear and mitochondrial DNA: the role of male-biased dispersal in a Holarctic sea duck

Genetic studies of waterfowl (Anatidae) have observed the full spectrum of mitochondrial (mt) DNA population divergence, from apparent panmixia to deep, reciprocally monophyletic lineages. Yet, these studies often found weak or no nuclear (nu) DNA structure, which was often attributed to male‐biased gene flow, a common behaviour within this family. An alternative explanation for this ‘conflict’ is that the smaller effective population size and faster sorting rate of mtDNA relative to nuDNA lead to different signals of population structure. We tested these alternatives by sequencing 12 nuDNA introns for a Holarctic pair of waterfowl subspecies, the European goosander (Mergus merganser merganser) and the North American common merganser (M. m. americanus), which exhibit strong population structure in mtDNA. We inferred effective population sizes, gene flow and divergence times from published mtDNA sequences and simulated expected differentiation for nuDNA based on those histories. Between Europe and North America, nuDNA Ф_ST was 3.4‐fold lower than mtDNA Ф_ST, a result consistent with differences in sorting rates. However, despite geographically structured and monophyletic mtDNA lineages within continents, nuDNA Ф_ST values were generally zero and significantly lower than predicted. This between‐ and within‐continent contrast held when comparing mtDNA and nuDNA among published studies of ducks. Thus, male‐mediated gene flow is a better explanation than slower sorting rates for limited nuDNA differentiation within continents, which is also supported by nonmolecular data. This study illustrates the value of quantitatively testing discrepancies between mtDNA and nuDNA to reject the null hypothesis that conflict simply reflects different sorting rates.     

Population genetic structure, gene flow and sex-biased dispersal in frillneck lizards (Chlamydosaurus kingii)

By using both mitochondrial and nuclear multiloci markers, we explored population genetic structure, gene flow and sex-specific dispersal of frillneck lizards ( Chlamydosaurus kingii ) sampled at three locations, separated by 10 to 50 km, in a homogenous savannah woodland in tropical Australia. Apart from a recombinant lizard, the mitochondrial analyses revealed two nonoverlapping haplotypes/populations, while the nuclear markers showed that the frillneck lizards represented three separate clusters/populations. Due to the small population size of the mtDNA, fixation may occur via founder effects and/or drift. We therefore suggest that either of these two processes, or a combination of the two, are the most likely causes of the discordant results obtained from the mitochondrial and the nuclear markers. In contrast to the nonoverlapping mitochondrial haplotypes, in 12 out of 74 lizards, mixed nuclear genotypes were observed, hence revealing a limited nuclear gene flow. Although gene flow should ultimately result in a blending of the populations, we propose that the distinct nuclear population structure is maintained by frequent fires resulting in local bottlenecks, and concomitant spatial separation of the frillneck lizard populations. Limited mark–recapture data and the difference in distribution of the mitochondrial and nuclear markers suggest that the mixed nuclear genotypes were caused by juvenile male-biased dispersal.     

Ecology and genetics of insecticide resistance inHelicoverpa armigera: Interactions between selection and gene flow

Pyrethroid resistance inHelicoverpa armigera provides a model system in which to study evolution in natural populations. Resistance is seen to evolve as a consequence of selection pressure that varies within and between life-stages and gene flow. Although three different mechanisms are involved, present day fluctuations in phenotype frequency can be explained by variation in only one of these, metabolic resistance, that is inherited as a single, incompletely dominant gene,mfo. Selective mortality of phenotypes occurs in both larvae and adults in the presence of the pyrethroid insecticides. Although most individuals of all three genotypes are killed in young larvae, selection in this age-class contributes significantly to evolution of resistance. While there is some evidence of reduced fitness of resistant pupae during winter diapause, most of the decline in resistance frequencies each spring occurs as a result of immigration of susceptible individuals into insecticide-treated populations.     

Lack of concordance between mtDNA gene flow and population density fluctuations in the bank vole

The genetic structure of bank voles Clethrionomys glareolus was determined from analyses of mitochondrial DNA (mtDNA) sequences, and compared with previous data on geographical synchrony in population density fluctuations. From 31 sample sites evenly spaced out along a 256-km transect in SE Norway a total of 39 distinct mtDNA haplotypes were found. The geographical distribution of the haplotypes was significantly nonrandom, and a cladistic analysis of the evolutionary relationship among haplotypes shows that descendant types were typically limited to a single site, whereas the ancestral types were more widely distributed geographically. This geographical distribution pattern of mtDNA haplotypes strongly indicates that the range and amount of female dispersal is severely restricted and insufficient to account for the previously observed synchrony in population density fluctuations. We conclude that geographical synchrony in this species must be caused by factors that are external to the local population, such as e.g. mobile predators.     

NM+: software implementing parentage-based models for estimating gene dispersal and mating patterns in plants

NM+ is computer software designed for making inferences on plant gene dispersal and mating patterns by modelling parentage probabilities of offspring based on the spatially explicit neighbourhood model. NM+ requires a sample of mapped and genotyped candidate parents and offspring; however, offspring may optionally be assigned to single maternal parents (forming so-called half-sib progeny arrays). Using maximum likelihood, NM+ estimates a number of parameters, including proportions of offspring due to self-fertilization, pollen immigration from outside of a defined study site, parameters of pollen (and/or seed) dispersal kernels (exponential-power, Weibull, geometric or 2Dt) and selection gradients relating covariates (phenotypic traits) with male (and/or female) reproductive success. NM+ allows for missing data both in parents and in offspring. It accounts for null alleles and their frequencies can optionally be considered as estimable parameters. Data files are formatted in a table-like structure so they can be easily prepared in a spreadsheet software. By default NM+ is for studying plant populations, however, it can be used for any organism as long as data requirements and model assumptions are met. NM+ runs under Windows, but it can be launched under Linux using WINE emulator. NM+ can be downloaded free of charge from http://www.genetyka.ukw.edu.pl/index_pliki/software.htm.     

Stable genetic polymorphism in heterogeneous environments: balance between asymmetrical dispersal and selection in the acorn barnacle

Elucidating the processes responsible for maintaining polymorphism at ecologically relevant genes is intimately related to understanding the interplay between selection imposed by habitat heterogeneity and a species' capacity for dispersal in the face of environmental constraints. In this paper, we used a model-based approach to solve equilibria of balanced polymorphism, given values of fitness and larval dispersal among different habitats in the acorn barnacle Semibalanus balanoides from the Gulf of St Lawrence. Our results showed that allele frequencies observed at both MPI* and GPI* loci represented stable equilibria, given empirical estimates of fitness values, and that considerably more larvae dispersed from one region (north) to the other (south) than vice versa. Dispersal conditions were predicted to be similar for the maintenance of polymorphism at both loci. Moreover, the values of asymmetrical dispersal required by the model to reach stable equilibria were compatible with empirical estimates of larval dispersal and oceanic circulation documented in this system. Overall, this study illustrated the usefulness of a modified and computable version of Bulmer's model (1972) in order to test hypotheses of balanced polymorphism resulting from interactions between spatial selection and asymmetrical dispersal.     

Islands and streams: clusters and gene flow in wild barley populations from the Levant

The domestication of plants frequently results in a high level of genetic differentiation between domesticated plants and their wild progenitors. This process is counteracted by gene flow between wild and domesticated plants because they are usually able to inter‐mate and to exchange genes. We investigated the extent of gene flow between wild barley Hordeum spontaneum and cultivated barley Hordeum vulgare, and its effect on population structure in wild barley by analysing a collection of 896 wild barley accessions (Barley1K) from Israel and all available Israeli H. vulgare accessions from the Israeli gene bank. We compared the performance of simple sequence repeats (SSR) and single nucleotide polymorphisms (SNP) marker data genotyped over a core collection in estimating population parameters. Estimates of gene flow rates with SSR markers indicated a high level of introgression from cultivated barley into wild barley. After removing accessions from the wild barley sample that were recently admixed with cultivated barley, the inference of population structure improved significantly. Both SSR and SNP markers showed that the genetic population structure of wild barley in Israel corresponds to the three major ecogeographic regions: the coast, the Mediterranean north and the deserts in the Jordan valley and the South. Gene flow rates were estimated to be higher from north to south than in the opposite direction. As has been observed in other crop species, there is a significant exchange of alleles between the wild species and domesticated varieties that needs to be accounted for in the population genetic analysis of domestication.     

Complete-fosmid and fosmid-end sequences reveal frequent horizontal gene transfers in marine uncultured planktonic archaea

The extent of horizontal gene transfer (HGT) among marine pelagic prokaryotes and the role that HGT may have played in their adaptation to this particular environment remain open questions. This is partly due to the paucity of cultured species and genomic information for many widespread groups of marine bacteria and archaea. Molecular studies have revealed a large diversity and relative abundance of marine planktonic archaea, in particular of Thaumarchaeota (also known as group I Crenarchaeota) and Euryarchaeota of groups II and III, but only one species (the thaumarchaeote Candidatus Nitrosopumilus maritimus) has been isolated in pure culture so far. Therefore, metagenomics remains the most powerful approach to study these environmental groups. To investigate the impact of HGT in marine archaea, we carried out detailed phylogenetic analyses of all open reading frames of 21 archaeal 16S rRNA gene-containing fosmids and, to extend our analysis to other genomic regions, also of fosmid-end sequences of 12 774 fosmids from three different deep-sea locations (South Atlantic and Adriatic Sea at 1000 m depth, and Ionian Sea at 3000 m depth). We found high HGT rates in both marine planktonic Thaumarchaeota and Euryarchaeota, with remarkable converging values estimated from complete-fosmid and fosmid-end sequence analysis (25 and 21% of the genes, respectively). Most HGTs came from bacterial donors (mainly from Proteobacteria, Firmicutes and Chloroflexi) but also from other archaea and eukaryotes. Phylogenetic analyses showed that in most cases HGTs are shared by several representatives of the studied groups, implying that they are ancient and have been conserved over relatively long evolutionary periods. This, together with the functions carried out by these acquired genes (mostly related to energy metabolism and transport of metabolites across membranes), suggests that HGT has played an important role in the adaptation of these archaea to the cold and nutrient-depleted deep marine environment.     

Sex-biased gene flow and colonization in the Formosan lesser horseshoe bat: inference from nuclear and mitochondrial markers

Sex-biased behaviours are expected to play an important role in partitioning genetic variance in animal populations. Comparing genetic structure at markers with different modes of inheritance provides a means of detecting these behaviours and their consequences for population genetic structure. In colonially breeding mammals, the common combination of female philopatry and male vagility can promote contrasting patterns of genetic differentiation between the sexes, both via their effects on recurrent gene flow and on colonization. We examined sex differences in gene flow and structure by comparing maternally inherited mitochondrial DNA (mtDNA) and biparentally inherited autosomal loci in the Formosan lesser horseshoe bat Rhinolophus monoceros . We found that genetic partitioning was higher at mtDNA than autosomal markers in both sexes, indicative of female-biased philopatry and male-biased dispersal. Across Taiwan, isolation-by-distance was detected for all sex/marker combinations but was steeper for mtDNA than for nuclear markers. We suggest that isolation-by-distance shown from mtDNA at large scales is likely to reflect the stepwise founding of new breeding colonies by females during colonization. In contrast, no isolation-by-distance was found at smaller distances of up to 100 km, indicating that gene flow and/or recent shared ancestry homogenises genetic structure among nearby sites. Our results highlight the value of an indirect genetic approach to understanding sex-biased behaviours and their consequences in a little-studied species.     

Linked vs unlinked markers: multilocus microsatellite haplotype-sharing as a tool to estimate gene flow and introgression

We have explored the use of multilocus microsatellite haplotypes to study introgression from cultivated (Malus domestica) into wild apple (Malus sylvestris), and to study gene flow among remnant populations of M. sylvestris. A haplotype consisted of alleles at microsatellite loci along one chromosome. As destruction of haplotypes through recombination occurs much faster than loss of alleles due to genetic drift, the lifespan of a multilocus haplotype is much shorter than that of the underlying alleles. When different populations share the same haplotype, this may indicate recent gene flow between populations. Similarly, haplotypes shared between two species would be a strong signal for introgression. As the expected lifespan of a haplotype depends on the strength of the linkage, the length [in centiMorgans (cM)] of the haplotype shared contains information on the number of generations passed. This application of shared haplotypes is distinct from using haplotype-sharing to detect association between markers and a certain trait. We inferred haplotypes for four to eight microsatellite loci on Linkage Group 10 of apple from genotype data using the program phase, and then identified those haplotypes shared between populations and species. Compared with a Bayesian analysis of unlinked microsatellite loci using the program structure, haplotype-sharing detected a partially different set of putative hybrids. Cultivated haplotypes present in M. sylvestris were short (< 1.5 cM), indicating that introgression had taken place many generations ago, except for two Belgian plants that contained a haplotype of 47.1 cM, indicating recent introgression. In the estimation of gene flow, F(ST) based on unlinked loci indicated small (0.032-0.058) but statistically significant differentiation between some populations only. However, various M. sylvestris haplotypes were shared in nearly all pairwise comparisons of populations, and their length indicated recent gene flow. Hence, all Dutch populations should be considered as one conservation unit. The added value of using sharing of multilocus microsatellite haplotypes as a source of population genetic information is discussed.     

Habitat heterogeneity, predation and gene flow: colour polymorphism in the isopod, Idotea baltica

The colour polymorphic isopod Idotea baltica inhabits the brown alga Fucus vesiculosus which is often colonised by the white epizoite Electra crustulenta (Bryozoa). In an experiment the predation risk for the different colour morphs of I. baltica was highly dependent on background colouration. Morph frequencies and Electra density varied substantially among 10 collecting sites but correlated poorly with each other, suggesting that local selection for cryptic colouration may be counteracted by gene flow. Indeed, estimates based on four polymorphic allozymes suggested rate of gene flow to be high. These results support the hypothesis that locally varying selection for cryptic colouration counteracted by gene flow contributes to the maintenance of colour polymorphism in I. baltica. The visual differences between the microhabitats and the differential microhabitat use between males and females seem to result in different patterns of selection on males and females for cryptic colouration. Also this is likely to play an important role for the polymorphism.     

Recent divergence with gene flow in Tennessee cave salamanders (Plethodontidae: Gyrinophilus) inferred from gene genealogies

Cave organisms occupy a special place in evolutionary biology because convergent morphologies of many species demonstrate repeatability in evolution even as they obscure phylogenetic relationships. The origin of specialized cave-dwelling species also raises the issue of the relative importance of isolation vs. natural selection in speciation. Two alternative hypotheses describe the origin of subterranean species. The 'climate-relict' model proposes allopatric speciation after populations of cold-adapted species become stranded in caves due to climate change. The 'adaptive-shift' model proposes parapatric speciation driven by divergent selection between subterranean and surface habitats. Our study of the Tennessee cave salamander complex shows that the three nominal forms (Gyrinophilus palleucus palleucus, G. p. necturoides, and G. gulolineatus) arose recently and are genealogically nested within the epigean (surface-dwelling) species, G. porphyriticus. Short branch lengths and discordant gene trees were consistent with a complex history involving gene flow between diverging forms. Results of coalescent-based analysis of the distribution of haplotypes among groups reject the allopatric speciation model and support continuous or recurrent genetic exchange during divergence. These results strongly favour the hypothesis that Tennessee cave salamanders originated from spring salamanders via divergence with gene flow.     

Social cohesion among kin,gene flow without dispersal and the evolution of population genetic structure in the killer whale (Orcinus orca)

In social species, breeding system and gregarious behavior are key factors influencing the evolution of large‐scale population genetic structure. The killer whale is a highly social apex predator showing genetic differentiation in sympatry between populations of foraging specialists (ecotypes), and low levels of genetic diversity overall. Our comparative assessments of kinship, parentage and dispersal reveal high levels of kinship within local populations and ongoing male‐mediated gene flow among them, including among ecotypes that are maximally divergent within the mtDNA phylogeny. Dispersal from natal populations was rare, implying that gene flow occurs without dispersal, as a result of reproduction during temporary interactions. Discordance between nuclear and mitochondrial phylogenies was consistent with earlier studies suggesting a stochastic basis for the magnitude of mtDNA differentiation between matrilines. Taken together our results show how the killer whale breeding system, coupled with social, dispersal and foraging behaviour, contributes to the evolution of population genetic structure.     

转基因水稻外源基因向近缘种群漂流和渐渗的研究进展

水稻是我国最重要的粮食作物之一,我国有8亿以上的人口以稻米作为主食。但在水稻生产中,由于病、虫、草害及不良气候等逆境因子的影响,严重制约了水稻的高产、稳产。转基因生物技术的迅速发展,为水稻抗性育种提供了新途径。自20世纪80年代以来,我国全方位地开展了转基因水稻的研发,目前已经培育出大量的抗病、抗虫、抗除草剂和抗逆的转基因水稻品种,这将为提高我国水稻的生产力和确保粮食安全做出重要的贡献。但转基因水稻的基因漂流及其可能带来的生物安全问题备受关注。已有报道证明,外源转基因可以通过异交向非转基因品种和野生近缘种漂流。在不同的试验条件下,抗除草剂基因有0．05％-0．53％逃逸的可能,其向不育系的最大漂移频率可达4．518％。抗虫基因向相邻非转基因水稻的平均漂移频率最高为0．875％。因此,本文对水稻与其近缘野生种的杂交情况,转基因水稻外源基因向非转基因品种、野生近缘种以及野生非近缘种的漂流和渐渗及其潜在的生态环境风险等方面进行了简要分析,并对转基因水稻的发展进行了展望,以期为转基因水稻的安全应用提供参考。     

Population history of Berthelot's pipit: colonization, gene flow and morphological divergence in Macaronesia

The fauna of oceanic islands provide exceptional models with which to examine patterns of dispersal, isolation and diversification, from incipient speciation to species level radiations. Here, we investigate recent differentiation and microevolutionary change in Berthelot's pipit (Anthus berthelotii), an endemic bird species inhabiting three Atlantic archipelagos. Mitochondrial DNA sequence data and microsatellite markers were used to deduce probable colonization pathway, genetic differentiation, and gene flow among the 12 island populations. Phenotypic differentiation was investigated based on eight biologically important morphological traits. We found little mitochondrial DNA variability, with only one and four haplotypes for the control region and cytochrome b, respectively. However, microsatellite data indicated moderate population differentiation (FST=0.069) between the three archipelagos that were identified as genetically distinct units with limited gene flow. Both results, combined with the estimated time of divergence (2.5 millions years ago) from the Anthus campestris (the sister species), suggest that this species has only recently dispersed throughout these islands. The genetic relationships, patterns of allelic richness and exclusive alleles among populations suggest the species originally colonized the Canary Islands and only later spread from there to the Madeiran archipelago and Selvagen Islands. Differentiation has also occurred within archipelagos, although to a lesser degree. Gene flow was observed more among the eastern and central islands of the Canaries than between these and the western islands or the Madeiran Islands. Morphological differences were also more important between than within archipelagos. Concordance between morphological and genetic differentiation provided ambiguous results suggesting that genetic drift alone was not sufficient to explain phenotypic differentiation. The observed genetic and morphological differences may therefore be the result of differing patterns of selection pressures between populations, with Berthelot's pipit undergoing a process of incipient differentiation.     

Pollen and seed dispersal inferred from seedling genotypes: the Bayesian revolution has passed here too

Understanding precisely how plants disperse their seeds and pollen in their neighbourhood is a central question for both ecologists and evolutionary biologists because seed and pollen dispersal governs both the rate of spread of an expanding population and gene flow within and among populations. The concept of a 'dispersal kernel' has become extremely popular in dispersal ecology as a tool that summarizes how dispersal distributes individuals and genes in space and at a given scale. In this issue of Molecular Ecology, the study by Moran & Clark (2011) (M&C in the following) shows how genotypic and spatial data of established seedlings can be analysed in a Bayesian framework to estimate jointly the pollen and seed dispersal kernels and finally derive a parentage analysis from a full-probability approach. This approach applied to red oak shows important dispersal of seeds (138 m on average) and pollen (178 m on average). For seeds, this estimate contrasts with previous results from inverse modelling on seed trap data (9.3 m). This research gathers several methodological advances made in recent years in two research communities and could become a cornerstone for dispersal ecology.     

Microsatellite diversity, population subdivision and gene flow in the Lipizzan horse

Blood samples of 561 Lipizzan horses from subpopulations (studs) of seven European countries representing a large fraction of the breed's population were used to examine the genetic diversity, population subdivision and gene flow in the breed. DNA analysis based on 18 microsatellite loci revealed that genetic diversity (observed heterozygosity = 0.663, gene diversity = 0.675 and the mean number of alleles = 7.056) in the Lipizzan horse is similar to other horse breeds as well as to other domestic animal species. The genetic differentiation between Lipizzan horses from different studs, although moderate, was apparent (pairwise F(ST) coefficients ranged from 0.021 to 0.080). Complementary findings explaining the genetic relationship among studs were revealed by genetic distance and principal component analysis. One genetic cluster consisted of the subpopulations of Austria, Italy and Slovenia, which represent the classical pool of Lipizzan horse breeding. A second cluster was formed by the Croatian, Hungarian and Slovakian subpopulations. The Romanian subpopulation formed a separate unit. The largest genetic differentiation was found between the Romanian and Italian subpopulation. Genetic results are consistent with the known breeding history of the Lipizzan horse. Correct stud assignment was obtained for 80.9% and 92.1% of Lipizzan horses depending on the inclusion or exclusion of migrant horses, respectively. The results of the present study will be useful for the development of breeding strategies, which consider classical horse breeding as well as recent achievements of population and conservation genetics.     

Genetic differentiation in pikeperch (Sander lucioperca): the relative importance of gene flow, drift and common history

To investigate the relative importance of homogenizing factors, such as gene flow, and diversifying factors, such as drift, genetic variation in pikeperch ( Sander lucioperca ) in two Fennoscandian regions (North and South) was analysed with microsatellites. Allelic richness and the degree of differentiation were significantly higher in the North ( F _ST= 0·20) than in the South ( F _ST= 0·064). In northern areas, assignments of genotypes were almost exclusively to the population of origin, but in southern areas, the proportion of correct assignments was significantly lower. Most samples exhibited significant heterozygote deficits, and the level of relatedness was higher than expected from randomness. These combined results suggest that there has been more gene flow between populations in southern areas than in northern areas, where the importance of genetic drift has been greater. Effective population sizes were small ( c. 100) and did not differ between areas. The effect of a common history appears minor, and thus processes such as genetic drift and gene flow have been more influential in shaping the patterns of genetic diversity in this species.