Consequences of introgression and gene flow on the genetic structure and diversity of Lima bean (Phaseolus lunatus L.) in its Mesoamerican diversity area.

We evaluated the role of gene flow and wild-crop introgression on the structure and genetic diversity of Lima bean (Phaseolus lunatus) in the Yucatan Peninsula, an important Mesoamerican diversity area for this crop, using a genotyping-by-sequencing approach (15,168 SNP markers) and two scales. At the local scale, STRUCTURE and NGSEP analyses showed predominantly crop-to-wild introgression, but also evidence of a bidirectional gene flow in the two wild-weedy-crop complexes studied (Itzinté and Dzitnup). The ABBA-BABA tests showed a higher introgression in Itzinté (the older complex) than in Dzitnup (the younger one); at the allelic level, the wild-crop introgression in Itzinté was similar in both directions, in Dzitnup it was higher from crop-to-wild; and at the chromosomal level, introgression in Itzinté was from wild-to-crop, whereas in Dzitnup it occured in the opposite direction. Also, we found H E values slightly higher in the domesticated accessions than in the wild ones, in both complexes (Itzinté: wild = 0.31, domesticated = 0.34; Dzinup: wild = 0.27, domesticated = 0.36), but %P and π estimators were higher in the wild accessions than in the domesticated ones. At a regional scale, STRUCTURE and MIGRATE showed a low gene flow, predominantly from crop-to-wild; and STRUCTURE, Neighbor-Joining and PCoA analyses indicated the existence of two wild groups and one domesticated group, with a marked genetic structure based in the existence of domesticated MI and wild MII gene pools. Also, at the regional scale, we found a higher genetic diversity in the wild accessions than in the domesticated ones, in all estimators used (e.g., H E = 0.27 and H E = 0.17, respectively). Our results indicate that gene flow and introgression are playing an important role at the local scale, but its consequences on the structure and genetic diversity of the Lima bean are not clearly reflected at the regional scale, where diversity patterns between wild and domesticated populations could be reflecting historical events.

High phylogeographic and genetic diversity of Tidestromia lanuginosa supports full-glacial refugia for arid-adapted plants in southern and central Coahuila, Mexico.

PREMISE: Recent phylogeographic work suggests the existence of latitudinal gradients in genetic diversity in northern Mexican plants, but very few studies have examined plants of the Chihuahuan Desert. Tidestromia lanuginosa is a morphologically variable annual species whose distribution includes the Chihuahuan Desert Region. Here we undertook phylogeographic analyses of chloroplast loci in this species to test whether genetic diversity and differentiation of Mexican populations of T. lanuginosa change along a latitudinal gradient and whether diversity is higher in Coahuila, consistent with ideas of lower plant community turnover during the Pleistocene. METHODS: Haplotype network, maximum likelihood tree, and Bayesian phylogenetic haplotype were reconstructed, and genetic diversity was assessed among 26 populations. Barrier analysis was used to explore barriers to gene flow. RESULTS: Four major population groups were identified, corresponding with physiographic provinces in Mexico. Each population group displayed high levels of genetic structure, haplotype, and nucleotide diversity. Diversity was highest in southern populations across the species as a whole and among the Chihuahuan Desert populations. CONCLUSIONS: Tidestromia lanuginosa provides an important example of high phylogeographic and genetic diversity in plants of northern Mexico. Barriers to gene flow among the major population groups have most likely been due to a combination of orographic, climatic, and edaphic variables. The high genetic diversity of T. lanuginosa in southern and central Coahuila is consistent with the hypothesis of full-glacial refugia for arid-adapted plants in this area, and highlights the importance of this region as a center of diversity for the Chihuahuan Desert flora.

Spatial detection of outlier loci with Moran eigenvector maps.

The spatial signature of microevolutionary processes structuring genetic variation may play an important role in the detection of loci under selection. However, the spatial location of samples has not yet been used to quantify this. Here, we present a new two-step method of spatial outlier detection at the individual and deme levels using the power spectrum of Moran eigenvector maps (MEM). The MEM power spectrum quantifies how the variation in a variable, such as the frequency of an allele at a SNP locus, is distributed across a range of spatial scales defined by MEM spatial eigenvectors. The first step (Moran spectral outlier detection: MSOD) uses genetic and spatial information to identify outlier loci by their unusual power spectrum. The second step uses Moran spectral randomization (MSR) to test the association between outlier loci and environmental predictors, accounting for spatial autocorrelation. Using simulated data from two published papers, we tested this two-step method in different scenarios of landscape configuration, selection strength, dispersal capacity and sampling design. Under scenarios that included spatial structure, MSOD alone was sufficient to detect outlier loci at the individual and deme levels without the need for incorporating environmental predictors. Follow-up with MSR generally reduced (already low) false-positive rates, though in some cases led to a reduction in power. The results were surprisingly robust to differences in sample size and sampling design. Our method represents a new tool for detecting potential loci under selection with individual-based and population-based sampling by leveraging spatial information that has hitherto been neglected.

Life history and past demography maintain genetic structure, outcrossing rate, contemporary pollen gene flow of an understory herb in a highly fragmented rainforest.

INTRODUCTION: Theory predicts that habitat fragmentation, by reducing population size and increasing isolation among remnant populations, can alter their genetic diversity and structure. A cascade of effects is expected: genetic drift and inbreeding after a population bottleneck, changes in biotic interactions that may affect, as in the case of plants, pollen dynamics, mating system, reproductive success. The detection of the effects of contemporary habitat fragmentation on the genetic structure of populations are conditioned by the magnitude of change, given the few number of generations since the onset of fragmentation, especially for long-lived organisms. However, the present-day genetic structure of populations may bear the signature of past demography events. Here, we examine the effects of rainforest fragmentation on the genetic diversity, population structure, mating system (outcrossing rate), indirect gene flow and contemporary pollen dynamics in the understory herb Aphelandra aurantiaca. Also, we assessed its present-day genetic structure under different past demographic scenarios. METHODS: Twelve populations of A. aurantiaca were sampled in large (4), medium (3), and small (5) forest fragments in the lowland tropical rainforest at Los Tuxtlas region. Variation at 11 microsatellite loci was assessed in 28-30 reproductive plants per population. In two medium- and two large-size fragments we estimated the density of reproductive plants, and the mating system by analyzing the progeny of different mother plants per population. RESULTS: Despite prevailing habitat fragmentation, populations of A. aurantiaca possess high genetic variation (He  = 0.61), weak genetic structure (Rst  = 0.037), and slight inbreeding in small fragments. Effective population sizes (Ne ) were large, but slightly lower in small fragments. Migrants derive mostly from large and medium size fragments. Gene dispersal is highly restricted but long distance gene dispersal events were detected. Aphelandra aurantiaca shows a mixed mating system (tm  = 0.81) and the outcrossing rate have not been affected by habitat fragmentation. A strong pollen pool structure was detected due to few effective pollen donors (Nep ) and low distance pollen movement, pointing that most plants received pollen from close neighbors. Past demographic fluctuations may have affected the present population genetic structure as Bayesian coalescent analysis revealed the signature of past population expansion, possibly during warmer conditions after the last glacial maximum. DISCUSSION: Habitat fragmentation has not increased genetic differentiation or reduced genetic diversity of A. aurantiaca despite dozens of generations since the onset of fragmentation in the region of Los Tuxtlas. Instead, past population expansion is compatible with the lack of observed genetic structure. The predicted negative effects of rainforest fragmentation on genetic diversity and population structure of A. aurantiaca seem to have been buffered owing to its large effective populations and long-distance dispersal events. In particular, its mixed-mating system, mostly of outcrossing, suggests high efficiency of pollinators promoting connectivity and reducing inbreeding. However, some results point that the effects of fragmentation are underway, as two small fragments showed higher membership probabilities to their population of origin, suggesting genetic isolation. Our findings underscore the importance of fragment size to maintain genetic connectivity across the landscape.

Genetic diversity and structure of wild populations of Carica papaya in Northern Mesoamerica inferred by nuclear microsatellites and chloroplast markers.

BACKGROUND AND AIMS: Few studies have evaluated the genetic structure and evolutionary history of wild varieties of important crop species. The wild papaya (Carica papaya) is a key element of early successional tropical and sub-tropical forests in Mexico, and constitutes the genetic reservoir for evolutionary potential of the species. In this study we aimed to determine how diverse and structured is the genetic variability of wild populations of C. papaya in Northern Mesoamerica. Moreover, we assessed if genetic structure and evolutionary history coincide with hypothetized (1) pre-Pleistocene events (Isthmus of Tehuantepec sinking), (2) Pleistocene refugia or (3) recent patterns. METHODS: We used six nuclear and two chloroplast (cp) DNA markers to assess the genetic diversity and phylogeographical structure of 19 wild populations of C. papaya in its natural distribution in Northern Mesoamerica. KEY RESULTS: We found high genetic diversity (Ho = 0·681 for nuclear markers, and h = 0·701 for cpDNA markers) and gene flow between populations of C. papaya (migration r up to 420 km). A lack of phylogeographical structure was found with the cpDNA markers (NST < GST), whereas a recent population structure was inferred with the nuclear markers. Evidence indicates that pre-Pleistocene events or refugia did not play an important role in the genetic structuring of wild papaya. CONCLUSIONS: Because of its life history characteristics and lack of an ancient phylogeographical structure found with the cpDNA markers, we suggest that C. papaya was dispersed throughout the lowland rain forests of Mexico (along the coastal plains and foothills of Sierras). This scenario supports the hypothesis that tropical forests in Northern Mesoamerica did not experience important climate fluctuations during the Pleistocene, and that the life history of C. papaya could have promoted long-distance dispersal and rapid colonization of lowland rainforests. Moreover, the results obtained with the nuclear markers suggest recent human disturbances. The fragmentation of tropical habitats in Northern Mesoamerica appears to be the main driver of genetic structuring, and the major threat to the dispersion and survival of the species in the wild.

Habitat fragmentation threatens wild populations of Carica papaya (Caricaceae) in a lowland rainforest.

• Premise of the study: Wild populations of domesticated species constitute a genetic reservoir and are fundamental to the evolutionary potential of species. Wild papaya (Carica papaya) is a rare, short-lived, gap-colonizing, dioecious tree that persists in the forest by continuous dispersal. Theoretically, these life-history characteristics render wild papaya highly susceptible to habitat fragmentation, with anticipated negative effects on its gene pool. Further, species dioecy may cause founder effects to generate local biases in sex ratio, decreasing effective population size.• Methods: We contrasted the genetic diversity and structure of C. papaya between wild populations from rainforest fragments and continuous forest at Los Tuxtlas, Mexico. We evaluated recent migration rates among populations as well as landscape resistance to gene flow. Finally, we calculated the sex ratio of the populations in both habitats.• Key results: Populations of wild papaya in rainforest fragments showed lower genetic diversity and higher population differentiation than populations in continuous rainforest. Estimates of recent migration rates showed a higher percentage of migrants moving from the continuous forest to the forest fragments than in the opposite direction. Agricultural land and cattle pasture were found to be the most resistant matrices to gene flow. Finally, biased sex ratios were seen to affect the effective population size in both habitats.• Conclusions: The mating system, rarity, and short life cycle of C. papaya are exacerbating the effects of rainforest fragmentation on its genetic diversity, threatening the persistence of its natural populations in the proposed place of origin as well as its genetic reservoir.