Gene-flow in a mosaic hybrid zone: is local introgression adaptive?

Genome-wide scans of genetic differentiation between hybridizing taxa can identify genome regions with unusual rates of introgression. Regions of high differentiation might represent barriers to gene flow, while regions of low differentiation might indicate adaptive introgression-the spread of selectively beneficial alleles between reproductively isolated genetic backgrounds. Here we conduct a scan for unusual patterns of differentiation in a mosaic hybrid zone between two mussel species, Mytilus edulis and M. galloprovincialis. One outlying locus, mac-1, showed a characteristic footprint of local introgression, with abnormally high frequency of edulis-derived alleles in a patch of M. galloprovincialis enclosed within the mosaic zone, but low frequencies outside of the zone. Further analysis of DNA sequences showed that almost all of the edulis allelic diversity had introgressed into the M. galloprovincialis background in this patch. We then used a variety of approaches to test the hypothesis that there had been adaptive introgression at mac-1. Simulations and model fitting with maximum-likelihood and approximate Bayesian computation approaches suggested that adaptive introgression could generate a "soft sweep," which was qualitatively consistent with our data. Although the migration rate required was high, it was compatible with the functioning of an effective barrier to gene flow as revealed by demographic inferences. As such, adaptive introgression could explain both the reduced intraspecific differentiation around mac-1 and the high diversity of introgressed alleles, although a localized change in barrier strength may also be invoked. Together, our results emphasize the need to account for the complex history of secondary contacts in interpreting outlier loci.

Asymmetric and differential gene introgression at a contact zone between two highly divergent lineages of field voles (Microtus agrestis).

Secondary contact zones have the potential to shed light on the mode and rate at which reproductive isolation accumulates during allopatric speciation. We investigated the population genetics of a contact zone between two highly divergent lineages of field voles (Microtus agrestis) in the Swiss Jura mountains. To shed light on the processes underlying introgression, we used maternally, paternally, and bi-parentally inherited markers. Though the two lineages maintained a strong genetic structure, we found some hybrids and evidence of gene flow. The extent of introgression varied with the mode of inheritance, being highest for mtDNA and absent for the Y chromosome. In addition, introgression was asymmetric, occurring only from the Northern to the Southern lineage. Both patterns seem parsimoniously explained by neutral processes linked to differences in effective sizes and sex-biased dispersal rates. The lineage with lower effective population size was also the more introgressed, and the mode-of-inheritance effect correlated with the male-biased dispersal rate of microtine rodents. We cannot exclude, however, that Haldane's effect contributed to the latter, as we found a marginally significant deficit in males (the heterogametic sex) among hybrids. We propose a possible demographic scenario to account for the patterns documented, and empirical extensions to further investigate this contact zone.

Gene flow and gene flux shape evolutionary patterns of variation in Drosophila subobscura

Gene flow (defined as allele exchange between populations) and gene flux (defined as allele exchange during meiosis in heterokaryotypic females) are important factors decreasing genetic differentiation between populations and inversions. Many chromosomal inversions are under strong selection and their role in recombination reduction enhances the maintenance of their genetic distinctness. Here we analyze levels and patterns of nucleotide diversity, selection and demographic history, using 37 individuals of Drosophila subobscura from Mount Parnes (Greece) and Barcelona (Spain). Our sampling focused on two frequent O-chromosome arrangements that differ by two overlapping inversions (OST and O3+4), which are differentially adapted to the environment as observed by their opposing latitudinal clines in inversion frequencies. The six analyzed genes (Pif1A, Abi, Sqd, Yrt, Atpa and Fmr1) were selected for their location across the O-chromosome and their implication in thermal adaptation. Despite the extensive gene flux detected outside the inverted region, significant genetic differentiation between both arrangements was found inside it. However, high levels of gene flow were detected for all six genes when comparing the same arrangement among populations. These results suggest that the adaptive value of inversions is maintained, regardless of the lack of genetic differentiation within arrangements from different populations, and thus favors the Local Adaptation hypothesis over the Coadapted Genome hypothesis as the basis of the selection acting on inversions in these populations.

Genome-wide evidence for speciation with gene flow in Heliconius butterflies

Most speciation events probably occur gradually, without complete and immediate reproductive isolation, but the full extent of gene flow between diverging species has rarely been characterized on a genome-wide scale. Documenting the extent and timing of admixture between diverging species can clarify the role of geographic isolation in speciation. Here we use new methodology to quantify admixture at different stages of divergence in Heliconius butterflies, based on whole-genome sequences of 31 individuals. Comparisons between sympatric and allopatric populations of H. melpomene, H. cydno, and H. timareta revealed a genome-wide trend of increased shared variation in sympatry, indicative of pervasive interspecific gene flow. Up to 40% of 100-kb genomic windows clustered by geography rather than by species, demonstrating that a very substantial fraction of the genome has been shared between sympatric species. Analyses of genetic variation shared over different time intervals suggested that admixture between these species has continued since early in speciation. Alleles shared between species during recent time intervals displayed higher levels of linkage disequilibrium than those shared over longer time intervals, suggesting that this admixture took place at multiple points during divergence and is probably ongoing. The signal of admixture was significantly reduced around loci controlling divergent wing patterns, as well as throughout the Z chromosome, consistent with strong selection for Müllerian mimicry and with known Z-linked hybrid incompatibility. Overall these results show that species divergence can occur in the face of persistent and genome-wide admixture over long periods of time.

Internal genetic structure and outcrossing rate in a natural population of Araucaria angustifolia (Bert.) O. Kuntze

The internal genetic structure and outcrossing rate of a population of Araucaria angustifolia (Bert.) O. Kuntze were investigated using 16 allozyme loci. Estimates of the mean number of alleles per loci (1.6), percentage of polymorphic loci (43.8%), and expected genetic diversity (0.170) were similar to those obtained for other gymnosperms. The analysis of spatial autocorrelation demonstrated the presence of internal structure in the first distance classes (up to 70 m), suggesting the presence of family structure. The outcrossing rate was high (0.956), as expected for a dioecious species. However, it was different from unity, indicating outcrossings between related individuals and corroborating the presence of internal genetic structure. The results of this study have implications for the methodologies used in conservation collections and for the use or analysis of this forest species. © The American Genetic Association. 2006. All rights reserved.

Study of variability and genetic structure of European populations of Myotis emarginatus and Myotis capaccinii (Chiroptera, Vespertilionidae)

The Geoffroy’s bat Myotis emarginatus is mainly present in southern, south-eastern and central Europe (Červerný, 1999) and is often recorded from northern Spain (Quetglas, 2002; Flaquer et al., 2004). It has demonstrated the species’ preference for forest. Myotis capaccinii, confined to the Mediterranean (Guille´n, 1999), is classified as ‘vulnerable’ on a global scale (Hutson, Mickleburgh & Racey, 2001). In general, the species preferred calm waters bordered by well-developed riparian vegetation and large (> 5 m) inter-bank distances (Biscardi et al. 2007). In this study we present the first results about population genetic structure of these two species of genus Myotis. We used two methods of sampling: invasive and non-invasive techniques. A total of 323 invasive samples and a total of 107 non-invasive samples were collected and analyzed. For Myotis emarginatus we have individuated for the first time a set of 7 microsatellites, which can work on this species, started from a set developed on Myotis myotis (Castella et al. 2000). We developed also a method for analysis of non-invasive samples, that given a good percentage of positive analyzed samples. The results have highlighted for the species Myotis emarginatus the presence on the European territory of two big groups, discovered by using the microsatellites tracers. On this species, 33 haplotypes of Dloop have been identified, some of them are presented only in some colonies. We identified respectively 33 haplotypes of Dloop and 10 of cytB for Myotis emarginatus and 25 of dloop and 15 of cytB for Myotis capaccinii. Myotis emarginatus’ results, both microsatellites and mtDNA, show that there is a strong genetic flow between different colonies across Europe. The results achieved on Myotis capaccinii are very interesting, in this case either for the microsatellites or the mitochondrial DNA sequences, and it has been highlighted a big difference between different colonies.

Speciation with gene-flow in eastern North American red oaks (Quercus section Lobatae)

As foundational species, oaks (Quercus : Fagaceae) support the activities of both humans and wildlife. However, many oaks in North America are declining, a crisis exacerbated by the previous disappearance of other hard mast-producing trees. In addition, the economic demands placed on this drought-tolerant group may intensify if climate change extirpates other, relatively mesophytic species. Genetic tools can help address these management challenges. To this end, we developed a suite of 27 microsatellite markers, of which 22 are derived from expressed sequence tags (ESTs). Many of these markers bear significant homology to known genes and may be able to directly assay functional genetic variation. Markers obtained from enriched microsatellite libraries, on the other hand, are typically located in heterochromatic regions and should reflect demographic processes. Considered jointly, genic and genomic microsatellites can elucidate patterns of gene-flow and natural selection, which are fundamental to both an organism's evolutionary ecology and conservation biology. To this end, we employed the developed markers in an FST-based genome scan to detect the signature of divergent selection among the red oaks (Quercus section Lobatae). Three candidate genes with putative roles in stress responses demonstrated patterns of diversity consistent with adaptation to heterogeneous selective pressures. These genes may be important in both local genetic adaptation within species and divergence among them. Next, we used an isolation-with-migration model to quantify levels of gene-flow among four red oaks species during speciation. Both speciation in allopatry and speciation with gene-flow were found to be major drivers of red oak biodiversity. Loci playing a key role in speciation are also likely to be ecologically important within species

Temporal genetic structure and relatedness in the Tufted Duck Aythya fuligula suggests limited kin association in winter

Conspecific aggregation of waterfowl in winter is a common example of animal flocking behaviour, yet patterns of relatedness and temporal substructure in such social groups remain poorly understood even in common species. A previous study based on mark-recapture data showed that Tufted Ducks Aythya fuligula caught on the same day were re-caught together in subsequent winters more often than expected by chance, suggesting stable assortments of ‘socially familiar’ individuals between wintering periods. The genetic relationships within these social groups were not clear. Based on 191 individuals genotyped at 10 microsatellite markers, we investigated the temporal genetic structure and patterns of relatedness among wintering Tufted Ducks at Lake Sempach, Switzerland, in two consecutive winters. We found no evidence of genetic differentiation between temporal groups within or between winters. The average levels of relatedness in temporal groups were low and not higher than expected in random assortments of individuals. However, Mantel tests performed for each sex separately revealed significant negative correlations between the pairwise relatedness coefficients and the number of days between the capture dates of pairs of wintering Tufted Duck in males and females. This pattern suggests the presence of a small number of co-migrating same-sex sibling pairs in wintering flocks of Tufted Ducks. Our findings provide one of the first genetic analyses of a common duck species outside the breeding season and contribute to the understanding of social interactions in long-distance migratory birds.

The natural killer gene complex genetic locus Chok encodes Ly-49D, a target recognition receptor that activates natural killing

Previously, we established that natural killer (NK) cells from C57BL/6 (B6), but not BALB/c, mice lysed Chinese hamster ovary (CHO) cells, and we mapped the locus that determines this differential CHO-killing capacity to the NK gene complex on chromosome 6. The localization of Chok in the NK gene complex suggested that it may encode either an activating or an inhibitory receptor. Here, results from a lectin-facilitated lysis assay predicted that Chok is an activating B6 NK receptor. Therefore, we immunized BALB/c mice with NK cells from BALB.B6–Cmv1r congenic mice and generated a mAb, designated 4E4, that blocked B6-mediated CHO lysis. mAb 4E4 also redirected lysis of Daudi targets, indicating its reactivity with an activating NK cell receptor. Furthermore, only the 4E4+ B6 NK cell subset mediated CHO killing, and this lysis was abrogated by preincubation with mAb 4E4. Flow cytometric analysis indicated that mAb 4E4 specifically reacts with Ly-49D but not Ly-49A, B, C, E, G, H, or I transfectants. Finally, gene transfer of Ly-49DB6 into BALB/c NK cells conferred cytotoxic capacity against CHO cells, thus establishing that the Ly-49D receptor is sufficient to activate NK cells to lyse this target. Hence, Ly-49D is the Chok gene product and is a mouse NK cell receptor capable of directly triggering natural killing.

Evidence that ultraviolet markings are associated with patterns of molecular gene flow

Recent studies have shown UV vision and markings to be important in vertebrates, particularly birds, where behavioral experiments have demonstrated its potential importance in sexual selection. However, there has been no genetic evidence that UV markings determine patterns of evolution among natural populations. Here we report molecular evidence that UV markings are associated with the pattern of gene flow in the Tenerife lizard (Gallotia galloti). This species has vicariance-induced, approximate east–west lineages in Tenerife closely congruent with the primary lineages of the sympatric gecko species. Against expectations, these molecular phylogeographic lineages (representing geological history) and isolation-by-distance do not appear to influence gene flow. Sexually mature males from populations either side of a latitudinal ecotone have different UV markings and gene flow appears to be linked to this difference in UV markings. It may be that these groups with different UV sexual markings mate assortatively, restricting the gene flow between them. This has implications for debate on the relative importance of vicariance and biotopes in influencing biodiversity, with this evidence supporting the latter.

Measurement of genetic structure within populations using Moran's spatial autocorrelation statistics.

Spatial structure of genetic variation within populations, an important interacting influence on evolutionary and ecological processes, can be analyzed in detail by using spatial autocorrelation statistics. This paper characterizes the statistical properties of spatial autocorrelation statistics in this context and develops estimators of gene dispersal based on data on standing patterns of genetic variation. Large numbers of Monte Carlo simulations and a wide variety of sampling strategies are utilized. The results show that spatial autocorrelation statistics are highly predictable and informative. Thus, strong hypothesis tests for neutral theory can be formulated. Most strikingly, robust estimators of gene dispersal can be obtained with practical sample sizes. Details about optimal sampling strategies are also described.

Gene-flow between populations of cotton bollworm Helicoverpa armigera (Lepidoptera: Noctuidae) is highly variable between years

Both large and small scale migrations of Helicoverpa armigera Hübner in Australia were investigated using AMOVA analysis and genetic assignment tests. Five microsatellite loci were screened across 3142 individuals from 16 localities in eight major cotton and grain growing regions within Australia, over a 38-month period (November 1999 to January 2003). From November 1999 to March 2001 relatively low levels of migration were characterized between growing regions. Substantially higher than average gene-flow rates and limited differentiation between cropping regions characterized the period from April 2001 to March 2002. A reduced migration rate in the year from April 2002 to March 2003 resulted in significant genetic structuring between cropping regions. This differentiation was established within two or three generations. Genetic drift alone is unlikely to drive genetic differentiation over such a small number of generations, unless it is accompanied by extreme bottlenecks and/or selection. Helicoverpa armigera in Australia demonstrated isolation by distance, so immigration into cropping regions is more likely to come from nearby regions than from afar. This effect was most pronounced in years with limited migration. However, there is evidence of long distance dispersal events in periods of high migration (April 2001-March 2002). The implications of highly variable migration patterns for resistance management are considered.

To self, or not to self ... A review of outcrossing and pollen-mediated gene flow in neotropical trees

Despite the typically low population densities and animal-mediated pollination of tropical forest trees, outcrossing and long-distance pollen dispersal are the norm. We reviewed the genetic literature on mating systems and pollen dispersal for neotropical trees to identify the ecological and phylogenetic correlates. The 36 studies surveyed found >90% outcrossed mating for 45 hermaphroditic or monoecious species. Self-fertilization rates varied inversely with population density and showed phylogenetic and geographic trends. The few direct measures of pollen flow (N = 11 studies) suggest that pollen dispersal is widespread among low-density tropical trees, ranging from a mean of 200 m to over 19 km for species pollinated by small insects or bats. Future research needs to examine (1) the effect of inbreeding depression on observed outcrossing rates, (2) pollen dispersal in a wide range of pollination syndromes and ecological classes, (3) and the range of variation of mating system expression at different hierarchical levels, including individual, seasonal, population, ecological, landscape and range wide.

Population dynamics and gene flow of Helicoverpa armigera (Lepidoptera : Noctuidae) on cotton and grain crops in the Murrumbidgee Valley, Australia

The population dynamics of Helicoverpa armigera (Hubner) (Lepidoptera: Noctuidae) in the Murrumbidgee Valley, Australia, has been characterized using five highly variable microsatellite loci. In the 2001-2002 growing season, there were very high levels of migration into the Murrumbidgee Valley with no detectable genetic structuring, consistent with previous analyses on a national scale. By contrast, there was significant genetic structuring over the 2002-2003 growing season, with three distinct genetic types detected. The first type corresponded to the first two generations and was derived from local individuals emerging from diapause and their progeny. The second genetic type corresponded to generation 3 and resulted from substantial immigration into the region. There was another genetic shift in generation 4, which accounts for the third genetic type of the season. This genetic shift occurred despite low levels of immigration. During the third generation of the 2002-2003 growing season, different population dynamics was characterized for H. armigera on maize, Zea mays L., and cotton Gossipium hirsutum L. Populations on cotton tended to cycle independently with very little immigration from outside the region or from maize within the region. Maize acted as a major sink for immigrants from cotton and from outside the region. If resistance were to develop on cotton under these circumstances, susceptible individuals from maize or from other regions would not dilute this resistance. In addition, resistance is likely to be transferred to maize and be perpetuated until diapause, from where it may reemerge next season. If low levels of immigration were to occur on transgenic cotton, this may undermine the effectiveness of refugia, especially noncotton refugia.

A comparative study of life history and gene flow in darters (Pisces: Percidae)

Gene flow, or the exchange of genes between populations, is important because it determines the evolutionary trajectory of a species, including the relative influences of genetic drift and natural selection in the process of population differentiation. Gene flow differs among species because of variation in dispersal capability and abundances across taxa, and historical forces related to geological or lineage history. Both history and ecology influence gene flow in potentially complicated ways, and accounting for their effects remains an important problem in evolutionary biology. This research is a comparative study of gene flow and life-history in a monophyletic group of stream fishes, the darters. As a first step in disentangling historical and ecological effects, I reconstructed the phylogenetic relationships of the study species from nucleotide sequences in the mtDNA control region. I then used this phylogeny and regional glaciation history to infer historical effects on life-history evolution and gene flow in 15 species of darters. Gene flow was estimated indirectly, using information from 20 resolvable and polymorphic allozyme loci. When I accounted for historical effects, comparisons across taxa revealed that gene flow rates were closely associated with differences in clutch sizes and reproductive investment patterns. I hypothesized that differences in larval dispersal among taxa explained this relationship. Results from a field study of larval drift were consistent with this hypothesis. Finally, I asked whether there was an interaction between species' ecology and genetic differentiation across biogeographically distinct regions. Information from allozymes and mtDNA sequences revealed that life history plays an important role in the magnitude of species divergence across biogeographic boundaries. These results suggested an important association between life histories and rates of speciation following an allopatric isolation event. This research, along with other studies from the literature, further illustrates the enormous potential of North American freshwater fishes as a system for studying speciation processes. ^