No evidence for female kin association, indications for extragroup paternity, and sex-biased dispersal patterns in wild western gorillas

Characterizing animal dispersal patterns and the rational behind individuals' transfer choices is a long-standing question of interest in evolutionary biology. In wild western gorillas (Gorilla gorilla), a one-male polygynous species, previous genetic findings suggested that, when dispersing, females might favor groups with female kin to promote cooperation, resulting in higher-than-expected within-group female relatedness. The extent of male dispersal remains unclear with studies showing conflicting results. To investigate male and female dispersal patterns and extragroup paternity, we analyzed long-term field observations, including female spatial proximity data, together with genetic data (10 autosomal microsatellites) on individuals from a unique set of four habituated western gorilla groups, and four additional extragroup males (49 individuals in total). The majority of offspring (25 of 27) were sired by the group male. For two offspring, evidence for extragroup paternity was found. Contrarily to previous findings, adult females were not significantly more related within groups than across groups. Consistently, adult female relatedness within groups did not correlate with their spatial proximity inferred from behavioral data. Adult females were similarly related to adult males from their group than from other groups. Using R ST statistics, we found significant genetic structure and a pattern of isolation by distance, indicating limited dispersal in this species. Comparing relatedness among females and among males revealed that males disperse farer than females, as expected in a polygamous species. Our study on habituated western gorillas shed light on the dispersal dynamics and reproductive behavior of this polygynous species and challenge some of the previous results based on unhabituated groups.

Polygenic Adaptation: Integrating Population Genetics and Gene Regulatory Networks

The adaptation of populations to local environments often relies on the selection of optimal values for polygenic traits. Here, we first summarize the results obtained from different quantitative genetics and population genetics models, about the genetic architecture of polygenic traits and their response to directional selection. We then highlight the contribution of systems biology to the understanding of the molecular bases of polygenic traits and the evolution of gene regulatory networks involved in these traits. Finally, we discuss the need for a unifying framework merging the fields of population genetics, quantitative genetics and systems biology to better understand the molecular bases of polygenic traits adaptation.

Sahelian pastoralism from the perspective of variants associated with lactase persistence

OBJECTIVES

Archeological evidence shows that first nomadic pastoralists came to the African Sahel from northeastern Sahara, where milking is reported by ~7.5 ka. A second wave of pastoralists arrived with the expansion of Arabic tribes in 7th-14th century CE. All Sahelian pastoralists depend on milk production but genetic diversity underlying their lactase persistence (LP) is poorly understood.

MATERIALS AND METHODS

We investigated SNP variants associated with LP in 1,241 individuals from 29 mostly pastoralist populations in the Sahel. Then, we analyzed six SNPs in the neighboring fragment (419 kb) in the Fulani and Tuareg with the -13910*T mutation, reconstructed haplotypes, and calculated expansion age and growth rate of this variant.

RESULTS

Our results reveal a geographic localization of two different LP variants in the Sahel: -13910*T west of Lake Chad (Fulani and Tuareg pastoralists) and -13915*G east of there (mostly Arabic-speaking pastoralists). We show that -13910*T has a more diversified haplotype background among the Fulani than among the Tuareg and that the age estimate for expansion of this variant among the Fulani (~8.5 ka) corresponds to introduction of cattle to the area.

CONCLUSIONS

This is the first study showing that the "Eurasian" LP allele -13910*T is widespread both in northern Europe and in the Sahel; however, it is limited to pastoralists in the Sahel. Since the Fulani haplotype with -13910*T is shared with contemporary Eurasians, its origin could be in a region encompassing the Near East and northeastern Africa in a population ancestral to both Saharan pastoralists and European farmers.

An ABC Method for Whole-Genome Sequence Data: Inferring Paleolithic and Neolithic Human Expansions

Species generally undergo a complex demographic history consisting, in particular, of multiple changes in population size. Genome-wide sequencing data are potentially highly informative for reconstructing this demographic history. A crucial point is to extract the relevant information from these very large data sets. Here, we design an approach for inferring past demographic events from a moderate number of fully sequenced genomes. Our new approach uses Approximate Bayesian Computation, a simulation-based statistical framework that allows 1) identifying the best demographic scenario among several competing scenarios and 2) estimating the best-fitting parameters under the chosen scenario. Approximate Bayesian Computation relies on the computation of summary statistics. Using a cross-validation approach, we show that statistics such as the lengths of haplotypes shared between individuals, or the decay of linkage disequilibrium with distance, can be combined with classical statistics (e.g., heterozygosity and Tajima's D) to accurately infer complex demographic scenarios including bottlenecks and expansion periods. We also demonstrate the importance of simultaneously estimating the genotyping error rate. Applying our method on genome-wide human-sequence databases, we finally show that a model consisting in a bottleneck followed by a Paleolithic and a Neolithic expansion is the most relevant for Eurasian populations.

The Genetic Architecture of Chronic Mountain Sickness in Peru

Chronic mountain sickness (CMS) is a pathological condition resulting from chronic exposure to high-altitude hypoxia. While its prevalence is high in native Andeans (>10%), little is known about the genetic architecture of this disease. Here, we performed the largest genome-wide association study (GWAS) of CMS (166 CMS patients and 146 controls living at 4,380 m in Peru) to detect genetic variants associated with CMS. We highlighted four new candidate loci, including the first CMS-associated variant reaching GWAS statistical significance (rs7304081; P = 4.58 × 10⁻⁹). By looking at differentially expressed genes between CMS patients and controls around these four loci, we suggested AEBP2, CAST, and MCTP2 as candidate CMS causal genes. None of the candidate loci were under strong natural selection, consistent with the observation that CMS affects fitness mainly after the reproductive years. Overall, our results reveal new insights on the genetic architecture of CMS and do not provide evidence that CMS-associated variants are linked to a strong ongoing adaptation to high altitude.

Neutral Theory: From Complex Population History to Natural Selection and Sociocultural Phenomena in Human Populations

Here, we present a synthetic view on how Kimura's Neutral theory has helped us gaining insight on the different evolutionary forces that shape human evolution. We put this perspective in the frame of recent emerging challenges: the use of whole genome data for reconstructing population histories, natural selection on complex polygenic traits, and integrating cultural processes in human evolution.

Cryptic Biological Invasions: a General Model of Hybridization

The dispersal of non-native genes due to hybridization is a form of cryptic invasion with growing concern in evolution and conservation. This includes the spread of transgenic genes and antibiotic resistance. To investigate how genes and phenotypes are transmitted, we developed a general model that, for the first time, considers concurrently: multiple loci, quantitative and qualitative gene expression, assortative mating, dominance/recessivity inheritance and density-dependent demographic effects. Selection acting on alleles or genotypes can also be incorporated. Our results reveal that the conclusions about how hybridization threatens a species can be biased if they are based on single-gene models, while considering two or more genes can correct this bias. We also show that demography can amplify or balance the genetic effects, evidencing the need of jointly incorporating both processes. By implementing our model in a real case, we show that mallard ducks introduced in New Zealand benefit from hybridization to replace native grey-ducks. Total displacement can take a few generations and occurs by interspecific competition and by competition between hybrids and natives, demonstrating how hybridization may facilitate biological invasions. We argue that our general model represents a powerful tool for the study of a wide range of biological and societal questions.

Genetic and linguistic histories in Central Asia inferred using approximate Bayesian computations

Linguistic and genetic data have been widely compared, but the histories underlying these descriptions are rarely jointly inferred. We developed a unique methodological framework for analysing jointly language diversity and genetic polymorphism data, to infer the past history of separation, exchange and admixture events among human populations. This method relies on approximate Bayesian computations that enable the identification of the most probable historical scenario underlying each type of data, and to infer the parameters of these scenarios. For this purpose, we developed a new computer program PopLingSim that simulates the evolution of linguistic diversity, which we coupled with an existing coalescent-based genetic simulation program, to simulate both linguistic and genetic data within a set of populations. Applying this new program to a wide linguistic and genetic dataset of Central Asia, we found several differences between linguistic and genetic histories. In particular, we showed how genetic and linguistic exchanges differed in the past in this area: some cultural exchanges were maintained without genetic exchanges. The methodological framework and the linguistic simulation tool developed here can be used in future work for disentangling complex linguistic and genetic evolutions underlying human biological and cultural histories.

Trait level analysis of multitrait population projection matrices

In most matrix population projection models, individuals are characterized according to, usually, one or two traits such as age, stage, size or location. A broad theory of multitrait population projection matrices (MPPMs) incorporating larger number of traits was long held back by time and space computational complexity issues. As a consequence, no study has yet focused on the influence of the structure of traits describing a life-cycle on population dynamics and life-history evolution. We present here a novel vector-based MPPM building methodology that allows to computationally-efficiently model populations characterized by numerous traits with large distributions, and extend sensitivity analyses for these models. We then present a new method, the trait level analysis consisting in folding an MPPM on any of its traits to create a matrix with alternative trait structure (the number of traits and their characteristics) but similar asymptotic properties. Adding or removing one or several traits to/from the MPPM and analyzing the resulting changes in spectral properties, allows investigating the influence of the trait structure on the evolution of traits. We illustrate this by modeling a 3-trait (age, parity and fecundity) population designed to investigate the implications of parity-fertilitytrade-offs in a context of fecundity heterogeneity in humans. The trait level analysis, comparing models of the same population differing in trait structures, demonstrates that fertility selection gradients differ between cases with or without parity-fertility trade-offs. Moreover it shows that age-specific fertility has seemingly very different evolutionary significance depending on whether heterogeneity is accounted for. This is because trade-offs can vary strongly in strength and even direction depending on the trait structure used to model the population.

Different kinds of genetic markers permit inference of Paleolithic and Neolithic expansions in humans

Recent population genetic studies have provided valuable insights on the demographic history of our species. However, some issues such as the dating of the first demographic expansions in human populations remain puzzling. Indeed, although a few genetic studies argued that the first human expansions were concomitant with the Neolithic transition, many others found signals of expansion events starting during the Palaeolithic. Here we performed a simulation study to show that these contradictory findings may result from the differences in the genetic markers used, especially if two successive expansion events occurred. For a large majority of replicates for each scenario tested, microsatellite data allow only detecting the recent expansion event in that case, whereas sequence data allow only detecting the ancient expansion. Combined with previous real data analyses, our results bring support to the ideas that (i) a first human expansions started during the Palaeolithic period, (ii) a second expansion event occurred later, concomitantly with the Neolithic transition.

Inferring Population Size History from Large Samples of Genome-Wide Molecular Data - An Approximate Bayesian Computation Approach

Inferring the ancestral dynamics of effective population size is a long-standing question in population genetics, which can now be tackled much more accurately thanks to the massive genomic data available in many species. Several promising methods that take advantage of whole-genome sequences have been recently developed in this context. However, they can only be applied to rather small samples, which limits their ability to estimate recent population size history. Besides, they can be very sensitive to sequencing or phasing errors. Here we introduce a new approximate Bayesian computation approach named PopSizeABC that allows estimating the evolution of the effective population size through time, using a large sample of complete genomes. This sample is summarized using the folded allele frequency spectrum and the average zygotic linkage disequilibrium at different bins of physical distance, two classes of statistics that are widely used in population genetics and can be easily computed from unphased and unpolarized SNP data. Our approach provides accurate estimations of past population sizes, from the very first generations before present back to the expected time to the most recent common ancestor of the sample, as shown by simulations under a wide range of demographic scenarios. When applied to samples of 15 or 25 complete genomes in four cattle breeds (Angus, Fleckvieh, Holstein and Jersey), PopSizeABC revealed a series of population declines, related to historical events such as domestication or modern breed creation. We further highlight that our approach is robust to sequencing errors, provided summary statistics are computed from SNPs with common alleles.

Molecular data sampled from extant individuals contains considerable information about their demographic history. In particular, one classical question in population genetics is to reconstruct past population size changes from such data. Relating these changes to various climatic, geological or anthropogenic events allows characterizing the main factors driving genetic diversity and can have major outcomes for conservation. Until recently, mostly very simple histories, including one or two population size changes, could be estimated from genetic data. This has changed with the sequencing of entire genomes in many species, and several methods allow now inferring complex histories consisting of several tens of population size changes. However, analyzing entire genomes, while accounting for recombination, remains a statistical and numerical challenge. These methods, therefore, can only be applied to small samples with a few diploid genomes. We overcome this limitation by using an approximate estimation approach, where observed genomes are summarized using a small number of statistics related to allele frequencies and linkage disequilibrium. In contrast to previous approaches, we show that our method allows us to reconstruct also the most recent part (the last 100 generations) of the population size history. As an illustration, we apply it to large samples of whole-genome sequences in four cattle breeds.

Using genome-wide measures of coancestry to maintain diversity and fitness in endangered and domestic pig populations

Conservation and breeding programs aim at maintaining the most diversity, thereby avoiding deleterious effects of inbreeding while maintaining enough variation from which traits of interest can be selected. Theoretically, the most diversity is maintained using optimal contributions based on many markers to calculate coancestries, but this can decrease fitness by maintaining linked deleterious variants. The heterogeneous patterns of coancestry displayed in pigs make them an excellent model to test these predictions. We propose methods to measure coancestry and fitness from resequencing data and use them in population management. We analyzed the resequencing data of Sus cebifrons, a highly endangered porcine species from the Philippines, and genotype data from the Pietrain domestic breed. By analyzing the demographic history of Sus cebifrons, we inferred two past bottlenecks that resulted in some inbreeding load. In Pietrain, we analyzed signatures of selection possibly associated with commercial traits. We also simulated the management of each population to assess the performance of different optimal contribution methods to maintain diversity, fitness, and selection signatures. Maximum genetic diversity was maintained using marker-by-marker coancestry, and least using genealogical coancestry. Using a measure of coancestry based on shared segments of the genome achieved the best results in terms of diversity and fitness. However, this segment-based management eliminated signatures of selection. We demonstrate that maintaining both diversity and fitness depends on the genomic distribution of deleterious variants, which is shaped by demographic and selection histories. Our findings show the importance of genomic and next-generation sequencing information in the optimal design of breeding or conservation programs.

Statistical inference on genetic data reveals the complex demographic history of human populations in central Asia

The demographic history of modern humans constitutes a combination of expansions, colonizations, contractions, and remigrations. The advent of large scale genetic data combined with statistically refined methods facilitates inference of this complex history. Here we study the demographic history of two genetically admixed ethnic groups in Central Asia, an area characterized by high levels of genetic diversity and a history of recurrent immigration. Using Approximate Bayesian Computation, we infer that the timing of admixture markedly differs between the two groups. Admixture in the traditionally agricultural Tajiks could be dated back to the onset of the Neolithic transition in the region, whereas admixture in Kyrgyz is more recent, and may have involved the westward movement of Turkic peoples. These results are confirmed by a coalescent method that fits an isolation-with-migration model to the genetic data, with both Central Asian groups having received gene flow from the extremities of Eurasia. Interestingly, our analyses also uncover signatures of gene flow from Eastern to Western Eurasia during Paleolithic times. In conclusion, the high genetic diversity currently observed in these two Central Asian peoples most likely reflects the effects of recurrent immigration that likely started before historical times. Conversely, conquests during historical times may have had a relatively limited genetic impact. These results emphasize the need for a better understanding of the genetic consequences of transmission of culture and technological innovations, as well as those of invasions and conquests.

Inference of sex-specific expansion patterns in human populations from Y-chromosome polymorphism

Studying the current distribution of genetic diversity in humans has important implications for our understanding of the history of our species. We analyzed a set of linked STR and SNP loci from the paternally inherited Y chromosome to infer the past demography of 55 African and Eurasian populations, using both the parametric and nonparametric coalescent-based methods implemented in the BEAST application. We inferred expansion events in most sedentary farmer populations, while we found constant effective population sizes for both nomadic hunter-gatherers and seminomadic herders. Our results differed, on several aspects, from previous results on mtDNA and autosomal markers. First, we found more recent expansion patterns in Eurasia than in Africa. This discrepancy, substantially stronger than the ones found with the other kind of markers, may result from a lower effective population size for men, which might have made male-transmitted markers more sensitive to the out-of-Africa bottleneck. Second, we found expansion signals only for sedentary farmers but not for nomadic herders in Central Asia, while these signals were found for both kind of populations in this area when using mtDNA or autosomal markers. Expansion signals in this area may result from spatial expansion processes and may have been erased for the Y chromosome among the herders because of restricted male gene flow.

Microsatellite data show recent demographic expansions in sedentary but not in nomadic human populations in Africa and Eurasia

The transition from hunting and gathering to plant and animal domestication was one of the most important cultural and technological revolutions in human history. According to archeologists and paleoanthropologists, this transition triggered major demographic expansions. However, few genetic studies have found traces of Neolithic expansions in the current repartition of genetic polymorphism, pointing rather toward Paleolithic expansions. Here, we used microsatellite autosomal data to investigate the past demographic history of 87 African and Eurasian human populations with contrasted lifestyles (nomadic hunter-gatherers, semi-nomadic herders and sedentary farmers). Likely due to the combination of a higher mutation rate and the possibility to analyze several loci as independent replicates of the coalescent process, the analysis of microsatellite data allowed us to infer more recent expansions than previous genetic studies, potentially resulting from the Neolithic transition. Despite the variability in their location and environment, we found consistent expansions for all sedentary farmers, while we inferred constant population sizes for all hunter-gatherers and most herders that could result from constraints linked to a nomadic or semi-nomadic lifestyle and/or competition for land between herders and farmers. As an exception, we inferred expansions for Central Asian herders. This might be linked with the arid environment of this area that may have been more favorable to nomadic herders than to sedentary farmers. Alternatively, current Central Asian herders may descent from populations who have first experienced a transition from hunter-gathering to sedentary agropastoralism, and then a second transition to nomadic herding.

Human genetic data reveal contrasting demographic patterns between sedentary and nomadic populations that predate the emergence of farming

Demographic changes are known to leave footprints on genetic polymorphism. Together with the increased availability of large polymorphism data sets, coalescent-based methods allow inferring the past demography of populations from their present-day patterns of genetic diversity. Here, we analyzed both nuclear (20 noncoding regions) and mitochondrial (HVS-I) resequencing data to infer the demographic history of 66 African and Eurasian human populations presenting contrasting lifestyles (nomadic hunter-gatherers, nomadic herders, and sedentary farmers). This allowed us to investigate the relationship between lifestyle and demography and to address the long-standing debate about the chronology of demographic expansions and the Neolithic transition. In Africa, we inferred expansion events for farmers, but constant population sizes or contraction events for hunter-gatherers. In Eurasia, we inferred higher expansion rates for farmers than herders with HVS-I data, except in Central Asia and Korea. Although isolation and admixture processes could have impacted our demographic inferences, these processes alone seem unlikely to explain the contrasted demographic histories inferred in populations with different lifestyles. The small expansion rates or constant population sizes inferred for herders and hunter-gatherers may thus result from constraints linked to nomadism. However, autosomal data revealed contraction events for two sedentary populations in Eurasia, which may be caused by founder effects. Finally, the inferred expansions likely predated the emergence of agriculture and herding. This suggests that human populations could have started to expand in Paleolithic times, and that strong Paleolithic expansions in some populations may have ultimately favored their shift toward agriculture during the Neolithic.

Consecutive five-year analysis of paternal and maternal gene flow and contributions of gametic heterogeneities to overall genetic composition of dispersed seeds of Pinus densiflora (Pinaceae)

PREMISE OF THE STUDY

Genetic variability in monoecious woody plant populations results from the assemblage of individuals issued from asymmetrical male and female reproductive functions, produced during spatially and temporarily heterogeneous reproductive and dispersal events. Here we investigated the dispersal patterns and levels of genetic diversity and differentiation of both paternal and maternal gametes in a natural population of Pinus densiflora at the multiple-year scale as long as five consecutive years. •

METHODS

We analyzed the paternity and maternity for 1576 seeds and 454 candidate adult trees using nuclear DNA polymorphisms of diploid biparental embryos and haploid maternal megagametophytes at eight microsatellite loci. •

KEY RESULTS

Despite the low levels of genetic differentiation among gamete groups, a two-way AMOVA analysis showed that the parental origin (paternal vs. maternal gametes), the year of gamete production and their interaction had significant effects on the genetic composition of the seeds. While maternal gamete groups showed a significant FST value across the 5 years, this was not true for their paternal counterparts. Within the population, we found that the relative reproductive contributions of the paternal vs. the maternal parent differed among adult trees, the maternal contributions showing a larger year-to-year fluctuation. •

CONCLUSIONS

The overall genetic variability of dispersed seeds appeared to result from two sources of heterogeneity: the difference between paternal and maternal patterns of reproduction and gamete dispersal and year-to-year heterogeneity of reproduction of adult trees, especially in their maternal reproduction.

Sociocultural behavior, sex-biased admixture, and effective population sizes in Central African Pygmies and non-Pygmies

Sociocultural phenomena, such as exogamy or phylopatry, can largely determine human sex-specific demography. In Central Africa, diverging patterns of sex-specific genetic variation have been observed between mobile hunter-gatherer Pygmies and sedentary agricultural non-Pygmies. However, their sex-specific demography remains largely unknown. Using population genetics and approximate Bayesian computation approaches, we inferred male and female effective population sizes, sex-specific migration, and admixture rates in 23 Central African Pygmy and non-Pygmy populations, genotyped for autosomal, X-linked, Y-linked, and mitochondrial markers. We found much larger effective population sizes and migration rates among non-Pygmy populations than among Pygmies, in agreement with the recent expansions and migrations of non-Pygmies and, conversely, the isolation and stationary demography of Pygmy groups. We found larger effective sizes and migration rates for males than for females for Pygmies, and vice versa for non-Pygmies. Thus, although most Pygmy populations have patrilocal customs, their sex-specific genetic patterns resemble those of matrilocal populations. In fact, our results are consistent with a lower prevalence of polygyny and patrilocality in Pygmies compared with non-Pygmies and a potential female transmission of reproductive success in Pygmies. Finally, Pygmy populations showed variable admixture levels with the non-Pygmies, with often much larger introgression from male than from female lineages. Social discrimination against Pygmies triggering complex movements of spouses in intermarriages can explain these male-biased admixture patterns in a patrilocal context. We show how gender-related sociocultural phenomena can determine highly variable sex-specific demography among populations, and how population genetic approaches contrasting chromosomal types allow inferring detailed human sex-specific demographic history.

Lactase persistence in central Asia: phenotype, genotype, and evolution

The aim of the present study is to document the evolution of the lactase persistence trait in Central Asia, a geographical area that is thought to have been a region of long-term pastoralism. Several ethnic groups co-exist in this area: Indo-Iranian speakers who are traditionally agriculturist (Tajik) and Turkic speakers who used to be nomadic herders (Kazakh, Karakalpak, Kyrgyz, Turkmen). It was recently demonstrated that horse milking practice existed in the Botai culture of Kazakhstan as early as 5,500 BP ( Outram et al. 2009 ). However, the frequency of the lactase persistence trait and its genetic basis in Central Asian populations remain largely unknown. We propose here the first genotype-phenotype study of lactase persistence in Central Asia based on 183 individuals, as well as the estimation of the time of expansion of the lactase-persistence associated polymorphism. Our results show a remarkable genetic-phenotypic correlation, with the causal polymorphism being the same than in Europe (-13.910C>T, rs4988235). The lactase persistence trait is at low frequency in these populations: between 25% and 32% in the Kazakh population (traditionally herders), according to phenotype used, and between 11% and 30% in the Tajiko-Uzbek population (agriculturalists). The difference in lactase persistence between populations, even if small, is significant when using individuals concordant for both excretion of breath hydrogen and the lactose tolerance blood glucose test phenotypes (P = 0.018, 25% for Kazakh vs. 11% for Tajiko-Uzbeks), and the difference in frequency of the -13.910*T allele is almost significant (P = 0.06, 30% for Kazakhs vs. 19% for Tajiko-Uzbeks). Using the surrounding haplotype, we estimate a date of expansion of the T allele around 6,000-12,000 yrs ago, which is consistent with archaeological records for the emergence of agropastoralism and pastoralism in Central Asia.

Impact of selection on genes involved in regulatory network: a modelling study

Complex phenotypes are often controlled by many interacting genes. One question emerging from such organization is how selection, acting at the phenotypic level, shapes the evolution of genes involved in regulatory networks controlling the phenotypes. We studied this issue through a matrix model of such networks. In a population submitted to selection, we simulated the evolution of a quantitative trait controlled by a set of loci that regulate each other through positive or negative interactions. Investigating several levels of selection intensity on the trait, we studied the evolution of regulation intensity between the genes and the evolution of the genetic diversity of those genes as an indirect measure of the strength of selection acting on them. We show that an increasing intensity of selection on the phenotype leads to an increased level of regulation between the loci. Moreover, we found that the genes responding more strongly to selection within the network were those evolving towards stronger regulatory action on the other genes and/or those that are the less regulated by the other genes. This observation is strongest for an intermediate level of selection. This may explain why several experimental studies have shown evidence of selection on regulatory genes inside gene networks.

The effects of inbreeding, genetic dissimilarity and phenotype on male reproductive success in a dioecious plant

Pollen fate can strongly affect the genetic structure of populations with restricted gene flow and significant inbreeding risk. We established an experimental population of inbred and outbred Silene latifolia plants to evaluate the effects of (i) inbreeding depression, (ii) phenotypic variation and (iii) relatedness between mates on male fitness under natural pollination. Paternity analysis revealed that outbred males sired significantly more offspring than inbred males. Independently of the effects of inbreeding, male fitness depended on several male traits, including a sexually dimorphic (flower number) and a gametophytic trait (in vitro pollen germination rate). In addition, full-sib matings were less frequent than randomly expected. Thus, inbreeding, phenotype and genetic dissimilarity simultaneously affect male fitness in this animal-pollinated plant. While inbreeding depression might threaten population persistence, the deficiency of effective matings between sibs and the higher fitness of outbred males will reduce its occurrence and counter genetic erosion.

A homoallelic FECH mutation in a patient with both erythropoietic protoporphyria and palmar keratoderma

BACKGROUND

Erythropoietic protoporphyria (EPP) is a hereditary disorder caused by the deficiency of ferrochelatase (FECH) in the haem biosynthetic pathway. In the majority of families, EPP is transmitted as a pseudodominant trait. Autosomal recessive form of EPP is found in only about 3% of the families.

OBJECTIVES

In this study, we describe a 6-year-old boy who suffered from both EPP and palmar keratoderma.

METHODS AND RESULTS

A novel homoallelic missense mutation (p.Ser318Tyr) was identified in the FECH gene. In addition, a region of homozygosity of approximately 6.8 Mb was observed in chromosome 18 of the patient by both microsatellite and SNP array. The parents of the patient, both of Palestinian (Jordanian) origin, were heterozygous for the S318Y mutation, although no history of consanguinity was known. Microsatellite genotyping identified a partial haplotype from each parent that corresponds to the region of homozygosity in the patient. Assuming S318Y is a founder mutation, the number of generations separating the two parents from their common ancestor from whom they inherited S318Y was estimated as 21.7 (95% CI 3.42–69.7).

CONCLUSION

EPP was therefore inherited as an autosomal recessive trait in the family. This study confirms the association between palmar keratoderma and autosomal recessive EPP.

In the heartland of Eurasia: the multilocus genetic landscape of Central Asian populations

Located in the Eurasian heartland, Central Asia has played a major role in both the early spread of modern humans out of Africa and the more recent settlements of differentiated populations across Eurasia. A detailed knowledge of the peopling in this vast region would therefore greatly improve our understanding of range expansions, colonizations and recurrent migrations, including the impact of the historical expansion of eastern nomadic groups that occurred in Central Asia. However, despite its presumable importance, little is known about the level and the distribution of genetic variation in this region. We genotyped 26 Indo-Iranian- and Turkic-speaking populations, belonging to six different ethnic groups, at 27 autosomal microsatellite loci. The analysis of genetic variation reveals that Central Asian diversity is mainly shaped by linguistic affiliation, with Turkic-speaking populations forming a cluster more closely related to East-Asian populations and Indo-Iranian speakers forming a cluster closer to Western Eurasians. The scattered position of Uzbeks across Turkic- and Indo-Iranian-speaking populations may reflect their origins from the union of different tribes. We propose that the complex genetic landscape of Central Asian populations results from the movements of eastern, Turkic-speaking groups during historical times, into a long-lasting group of settled populations, which may be represented nowadays by Tajiks and Turkmen. Contrary to what is generally thought, our results suggest that the recurrent expansions of eastern nomadic groups did not result in the complete replacement of local populations, but rather into partial admixture.

Fine-scale spatial genetic structure and gene dispersal in Silene latifolia

Plants are sessile organisms, often characterized by limited dispersal. Seeds and pollen are the critical stages for gene flow. Here we investigate spatial genetic structure, gene dispersal and the relative contribution of pollen vs seed in the movement of genes in a stable metapopulation of the white campion Silene latifolia within its native range. This short-lived perennial plant is dioecious, has gravity-dispersed seeds and moth-mediated pollination. Direct measures of pollen dispersal suggested that large populations receive more pollen than small isolated populations and that most gene flow occurs within tens of meters. However, these studies were performed in the newly colonized range (North America) where the specialist pollinator is absent. In the native range (Europe), gene dispersal could fall on a different spatial scale. We genotyped 258 individuals from large and small (15) subpopulations along a 60 km, elongated metapopulation in Europe using six highly variable microsatellite markers, two X-linked and four autosomal. We found substantial genetic differentiation among subpopulations (global F(ST)=0.11) and a general pattern of isolation by distance over the whole sampled area. Spatial autocorrelation revealed high relatedness among neighboring individuals over hundreds of meters. Estimates of gene dispersal revealed gene flow at the scale of tens of meters (5-30 m), similar to the newly colonized range. Contrary to expectations, estimates of dispersal based on X and autosomal markers showed very similar ranges, suggesting similar levels of pollen and seed dispersal. This may be explained by stochastic events of extensive seed dispersal in this area and limited pollen dispersal.

Genetic traces of east-to-west human expansion waves in Eurasia

In this study, we describe the landscape of human demographic expansions in Eurasia using a large continental Y chromosome and mitochondrial DNA dataset. Variation at these two uniparentally-inherited genetic systems retraces expansions that occurred in the past 60 ky, and shows a clear decrease of expansion ages from east to west Eurasia. To investigate the demographic events at the origin of this westward decrease of expansion ages, the estimated divergence ages between Eurasian populations are compared with the estimated expansion ages within each population. Both markers suggest that the demographic expansion diffused from east to west in Eurasia in a demic way, i.e., through migrations of individuals (and not just through diffusion of new technologies), highlighting the prominent role of eastern regions within Eurasia during Palaeolithic times.

Apparent mutational hotspots and long distance linkage disequilibrium resulting from a bottleneck

Genome wide patterns of nucleotide diversity and recombination reveal considerable variation including hotspots. Some studies suggest that these patterns are primarily dictated by individual locus history related at a broader scale to the population demographic history. Because bottlenecks have occurred in the history of numerous species, we undertook a simulation approach to investigate their impact on the patterns of aggregation of polymorphic sites and linkage disequilibrium (LD). We developed a new index (Polymorphism Aggregation Index) to characterize this aggregation and showed that variation in the density of polymorphic sites results from an interplay between the bottleneck scenario and the recombination rate. Under particular conditions, aggregation is maximized and apparent mutation hotspots resulting in a 50-fold increase in polymorphic sites density can occur. In similar conditions, long distance LD can be detected.

Long distance pollen-mediated gene flow at a landscape level: the weed beet as a case study

Gene flow is a crucial parameter that can affect the organization of genetic diversity in plant species. It has important implications in terms of conservation of genetic resources and of gene exchanges between crop to wild relatives and within crop species complex. In the Beta vulgaris complex, hybridization between crop and wild beets in seed production areas is well documented and the role of the ensuing hybrids, weed beets, as bridges towards wild forms in sugar beet production areas have been shown. Indeed, in contrast to cultivated beets that are bi-annual, weed beets can bolt, flower and reproduce in the same crop season. Nonetheless, the extent of pollen gene dispersal through weedy lineages remains unknown. In this study, the focus is directed towards weed-to-weed gene flow, and we report the results of a pollen-dispersal analysis within an agricultural landscape composed of five sugar beet fields with different levels of infestation by weed beets. Our results, based on paternity analysis of 3240 progenies from 135 maternal plants using 10 microsatellite loci, clearly demonstrate that even if weedy plants are mostly pollinated by individuals from the same field, some mating events occur between weed beets situated several kilometres apart (up to 9.6 km), with rates of interfield-detected paternities ranging from 11.3% to 17.5%. Moreover, we show that pollen flow appears to be more restricted when individuals are aggregated as most mating events occurred only for short-distance classes. The best-fit dispersal curves were fat-tailed geometric functions for populations exhibiting low densities of weed beets and thin-tailed Weibull function for fields with weed beet high densities. Thus, weed beet populations characterized by low density with geographically isolated individuals may be difficult to detect but are likely to act as pollen traps for pollen emitted by close and remote fields. Hence, it appears evident that interfield pollen-mediated gene flow between weed beets is almost unavoidable and could contribute to the diffusion of (trans)genes in the agricultural landscape.

Modelling and estimating pollen movement in oilseed rape (Brassica napus) at the landscape scale using genetic markers

Understanding patterns of pollen movement at the landscape scale is important for establishing management rules following the release of genetically modified (GM) crops. We use here a mating model adapted to cultivated species to estimate dispersal kernels from the genotypes of the progenies of male-sterile plants positioned at different sampling sites within a 10 x 10-km oilseed rape production area. Half of the pollen clouds sampled by the male-sterile plants originated from uncharacterized pollen sources that could consist of both large volunteer and feral populations, and fields within and outside the study area. The geometric dispersal kernel was the most appropriate to predict pollen movement in the study area. It predicted a much larger proportion of long-distance pollination than previously fitted dispersal kernels. This best-fitting mating model underestimated the level of differentiation among pollen clouds but could predict its spatial structure. The estimation method was validated on simulated genotypic data, and proved to provide good estimates of both the shape of the dispersal kernel and the rate and composition of pollen issued from uncharacterized pollen sources. The best dispersal kernel fitted here, the geometric kernel, should now be integrated into models that aim at predicting gene flow at the landscape level, in particular between GM and non-GM crops.

Estimating anisotropic pollen dispersal: a case study in Quercus lobata

The pollen dispersal distribution is an important element of the neighbourhood size of plant populations. Most methods aimed at estimating the dispersal curve assume that pollen dispersal is isotropic, but evidence indicates that this assumption does not hold for many plant species, particularly wind-pollinated species subject to prevailing winds during the pollination season. We propose here a method of detecting anisotropy of pollen dispersal and of gauging its intensity, based on the estimation of the differentiation of maternal pollen clouds (TWOGENER extraction), assuming that pollen dispersal is bivariate and normally distributed. We applied the new method to a case study in Quercus lobata, detecting only a modest level of anisotropy in pollen dispersal in a direction roughly similar to the prevailing wind direction. Finally, we conducted a simulation to explore the conditions under which anisotropy can be detected with this method, and we show that while anisotropy is detectable, in principle, it requires a large volume of data.

From social to genetic structures in central Asia

Pastoral and farmer populations, who have coexisted in Central Asia since the fourth millennium B.C., present not only different lifestyles and means of subsistence but also various types of social organization. Pastoral populations are organized into so-called descent groups (tribes, clans, and lineages) and practice exogamous marriages (a man chooses a bride in a different lineage or clan). In Central Asia, these descent groups are patrilineal: The children are systematically affiliated with the descent groups of the father. By contrast, farmer populations are organized into families (extended or nuclear) and often establish endogamous marriages with cousins. This study aims at better understanding the impact of these differences in lifestyle and social organization on the shaping of genetic diversity. We show that pastoral populations exhibit a substantial loss of Y chromosome diversity in comparison to farmers but that no such a difference is observed at the mitochondrial-DNA level. Our analyses indicate that the dynamics of patrilineal descent groups, which implies different male and female sociodemographic histories, is responsible for these sexually-asymmetric genetic patterns. This molecular signature of the pastoral social organization disappears over a few centuries only after conversion to an agricultural way of life.

Assortative mating and differential male mating success in an ash hybrid zone population

BACKGROUND

The structure and evolution of hybrid zones depend mainly on the relative importance of dispersal and local adaptation, and on the strength of assortative mating. Here, we study the influence of dispersal, temporal isolation, variability in phenotypic traits and parasite attacks on the male mating success of two parental species and hybrids by real-time pollen flow analysis. We focus on a hybrid zone population between the two closely related ash species Fraxinus excelsior L. (common ash) and F. angustifolia Vahl (narrow-leaved ash), which is composed of individuals of the two species and several hybrid types. This population is structured by flowering time: the F. excelsior individuals flower later than the F. angustifolia individuals, and the hybrid types flower in-between. Hybrids are scattered throughout the population, suggesting favorable conditions for their local adaptation. We estimate jointly the best-fitting dispersal kernel, the differences in male fecundity due to variation in phenotypic traits and level of parasite attack, and the strength of assortative mating due to differences in flowering phenology. In addition, we assess the effect of accounting for genotyping error on these estimations.

RESULTS

We detected a very high pollen immigration rate and a fat-tailed dispersal kernel, counter-balanced by slight phenological assortative mating and short-distance pollen dispersal. Early intermediate flowering hybrids, which had the highest male mating success, showed optimal sex allocation and increased selfing rates. We detected asymmetry of gene flow, with early flowering trees participating more as pollen donors than late flowering trees.

CONCLUSION

This study provides striking evidence that long-distance gene flow alone is not sufficient to counter-act the effects of assortative mating and selfing. Phenological assortative mating and short-distance dispersal can create temporal and spatial structuring that appears to maintain this hybrid population. The asymmetry of gene flow, with higher fertility and increased selfing, can potentially confer a selective advantage to early flowering hybrids in the zone. In the event of climate change, hybridization may provide a means for F. angustifolia to further extend its range at the expense of F. excelsior.

Real-time patterns of pollen flow in the wild-service tree, Sorbus torminalis (Rosaceae). III. Mating patterns and the ecological maternal neighborhood

Understanding the role of mother plants as pollen recipients in shaping mating patterns is essential for understanding the evolution of populations and in particular to predict the consequence of habitat fragmentation. Here, we investigated variation in mating patterns due to maternal phenotypic traits, phenological variance, and landscape features in Sorbus torminalis, a hermaphroditic, insect-pollinated and low-density, European temperate forest tree. The diversity and composition of pollen clouds received by maternal trees in S. torminalis were mainly determined by their conspecific neighborhood: isolated individuals sample more diversity through more even paternal contributions, low relatedness among paternal genes, and high rates of long-distance pollen dispersal within their progenies. Maternal phenotypic traits related to pollinator attractiveness also had an effect, but only when competition was strong: in this case, larger mother trees with more flowers sampled more diversity. The floral architecture of S. torminalis, with multiple-seeded fruit, strongly shaped mating patterns, with higher levels of correlated paternity among seeds belonging to the same fruit (30% full sibs) than among seeds belonging to different fruits (14% full sibs). Finally, flowering phenology affected the distribution of diversity among maternal pollen clouds, but the earliest and latest mother trees did not receive less diversity of pollen than the others.

Pollen Dispersal in Spatially Aggregated Populations

We perform a theoretical study of effective pollen dispersal within plant populations exhibiting intraspecific spatial aggregation. We simulate nonuniform distributions of individuals by means of a Poisson cluster process and use an individual-based spatially explicit model of pollen dispersal to assess the effects of different aggregation patterns on the effective pollen pool size (N(ep)) and the axial variance of pollen dispersal (sigma (p)). Results show clear interactions between clumping and both N(ep) and sigma (p), whose precise form and intensity depend on the relative spatial scale of aggregation to pollen dispersal range. If clump size is small relative to dispersal range, clumping results in lower N(ep) and sigma (p) than in randomly distributed populations. Interestingly, by contrast, aggregation may actually enlarge N(ep) and has minimum impact on sigma (p) if clump size is near or above the scale of dispersal. High intraclump to global density ratios enhance the sensitivity of both N(ep) and sigma (p) to clumping, while leptokurtic pollen dispersal generates sharper reductions of both N(ep) and sigma (p) for small clump sizes and stronger increments of N(ep) for larger clump sizes. Overall, our results indicate that isolation-by-distance models in plants should not ignore the effects of intraspecific spatial aggregation on effective dispersal.

A new method of estimating the pollen dispersal curve independently of effective density

We introduce a novel indirect method of estimating the pollen dispersal curve from mother-offspring genotypic data. Unlike an earlier indirect approach (TwoGener), this method is based on a normalized measure of correlated paternity between female pairs whose expectation does not explicitly depend on the unknown effective male population density (d(e)). We investigate the statistical properties of the new method, by comparison with those of TwoGener, considering the sensitivity to reductions of d(e), relative to census density, resulting from unequal male fecundity and asynchronous flowering. Our main results are: (i) it is possible to obtain reliable estimates of the average distance of pollen dispersal, delta, from indirect methods, even under nonuniform male fecundity and variable flowering phenology; (ii) the new method yields more accurate and more precise delta-estimates than TwoGener under a wide range of sampling and flowering scenarios; and (iii) TwoGener can be used to obtain approximate d(e) estimates, if needed for other purposes. Our results also show that accurately estimating the shape of the tail of the pollen dispersal function by means of indirect methods remains a very difficult challenge.

Importance of the life cycle in sympatric host race formation and speciation of pathogens

ABSTRACT Numerous morphological species of pathogenic fungi have been shown to actually encompass several genetically isolated lineages, often specialized on different hosts and, thus, constituting host races or sibling species. In this article, we explore theoretically the importance of some aspects of the life cycle on the conditions of sympatric divergence of host races, particularly in fungal plant pathogens. Because the life cycles classically modeled by theoreticians of sympatric speciation correspond to those of free-living animals, sympatric divergence of host races requires the evolution of active assortative mating or of active host preference if mating takes place on the hosts. With some particular life cycles with restricted dispersal between selection on the host and mating, we show that divergence can occur in sympatry and lead to host race formation, or even speciation, by a mere process of specialization, with strong divergent adaptive selection. Neither active assortative mating nor active habitat choice is required in these cases, and this may explain why the phylo-genetic species concept seems more appropriate than the biological species concept in these organisms.

Deciphering the ancient and complex evolutionary history of human arylamine N-acetyltransferase genes

The human N-acetyltransferase genes NAT1 and NAT2 encode two phase-II enzymes that metabolize various drugs and carcinogens. Functional variability at these genes has been associated with adverse drug reactions and cancer susceptibility. Mutations in NAT2 leading to the so-called slow-acetylation phenotype reach high frequencies worldwide, which questions the significance of altered acetylation in human adaptation. To investigate the role of population history and natural selection in shaping NATs variation, we characterized genetic diversity through the resequencing and genotyping of NAT1, NAT2, and the pseudogene NATP in a collection of 13 different populations with distinct ethnic backgrounds and demographic pasts. This combined study design allowed us to define a detailed map of linkage disequilibrium of the NATs region as well as to perform a number of sequence-based neutrality tests and the long-range haplotype (LRH) test. Our data revealed distinctive patterns of variability for the two genes: the reduced diversity observed at NAT1 is consistent with the action of purifying selection, whereas NAT2 functional variation contributes to high levels of diversity. In addition, the LRH test identified a particular NAT2 haplotype (NAT2*5B) under recent positive selection in western/central Eurasians. This haplotype harbors the mutation 341T-->C and encodes the "slowest-acetylator" NAT2 enzyme, suggesting a general selective advantage for the slow-acetylator phenotype. Interestingly, the NAT2*5B haplotype, which seems to have conferred a selective advantage during the past approximately 6,500 years, exhibits today the strongest association with susceptibility to bladder cancer and adverse drug reactions. On the whole, the patterns observed for NAT2 well illustrate how geographically and temporally fluctuating xenobiotic environments may have influenced not only our genome variability but also our present-day susceptibility to disease.

Pollen flow in the wildservice tree, Sorbus torminalis (L.) Crantz. II. Pollen dispersal and heterogeneity in mating success inferred from parent-offspring analysis

Knowing the extent of gene movements from parents to offspring is essential to understand the potential of a species to adapt rapidly to a changing environment, and to design appropriate conservation strategies. In this study, we develop a nonlinear statistical model to jointly estimate the pollen dispersal kernel and the heterogeneity in fecundity among phenotypically or environmentally defined groups of males. This model uses genotype data from a sample of fruiting plants, a sample of seeds harvested on each of these plants, and all males within a circumscribed area. We apply this model to a scattered, entomophilous woody species, Sorbus torminalis (L.) Crantz, within a natural population covering more than 470 ha. We estimate a high heterogeneity in male fecundity among ecological groups, both due to phenotype (size of the trees and flowering intensity) and landscape factors (stand density within the neighbourhood). We also show that fat-tailed kernels are the most appropriate to depict the important abilities of long-distance pollen dispersal for this species. Finally, our results reveal that the spatial position of a male with respect to females affects as much its mating success as ecological determinants of male fecundity. Our study thus stresses the interest to account for the dispersal kernel when estimating heterogeneity in male fecundity, and reciprocally.

Cultural transmission of fitness: genes take the fast lane

Classical population genetics describes how the fate of an allele is driven by four forces: mutation, migration, selection and drift. However, these are sometimes insufficient to explain how the observed allele frequency changes and, therefore, another factor must be invoked: cultural transmission of fitness (CTF). CTF is the non-genetic transmission of any kind of behaviour that affects reproductive success. There are several clearly documented examples of CTF, and theoretical studies have shown that it affects effective population size, linkage disequilibrium and coalescent times. It is therefore a factor that must be taken into account to explain the structure of genetic diversity. In this article, we will present documented cases of how CTF affects the genetic diversity of populations and yields dramatic changes in allele frequencies.

The genetic or mythical ancestry of descent groups: lessons from the Y chromosome

Traditional societies are often organized into descent groups called "lineages," "clans," and "tribes." Each of these descent groups claims to have a common ancestor, and this ancestry distinguishes the group's members from the rest of the population. To test the hypothesis of common ancestry within these groups, we compared ethnological and genetic data from five Central Asian populations. We show that, although people from the same lineage and clan share generally a recent common ancestor, no such common ancestry is observed at the tribal level. Thus, a tribe might be a conglomerate of clans who subsequently invented a mythical ancestor to strengthen group unity.

Vlax Roma history: what do coalescent-based methods tell us?

Three coalescent-based methods allowed us to infer some aspects of the history of three Bulgarian Gypsies populations belonging to the Vlax linguistic group: the Lom, Rudari and Kalderas. We used several kinds of genetic markers: HV1 sequences of the maternally inherited mitochondrial genome and microsatellites of the paternally inherited Y chromosome and of the biparentally inherited chromosome 8. This allowed us to infer several parameters for men and women: the splitting order of the populations and the ages of the splitting events, the growth rate in each population and the migration rates between populations. Altogether, they enabled us to infer a demographic scenario that could explain the genetic diversity of Vlax Roma: recent splits occurring after the arrival in Europe, asymmetric migration flows especially for males and unequal growth rates. This represents a considerable contribution to the Vlax Roma history in comparison with the inferences from classical population genetics.

Detecting Population Growth, Selection and Inherited Fertility From Haplotypic Data in Humans

The frequency of a rare mutant allele and the level of allelic association between this allele and one or several closely linked markers are frequently measured in genetic epidemiology. Both quantities are related to the time elapsed since the appearance of the mutation in the population and the intrinsic growth rate of the mutation (which may be different from the average population growth rate). Here, we develop a method that uses these two kinds of genetic data to perform a joint estimation of the age of the mutation and the minimum growth rate that is compatible with its present frequency. In absence of demographic data, it provides a useful estimate of population growth rate. When such data are available, contrasts among estimates from several loci allow demographic processes, affecting all loci similarly, to be distinguished from selection, affecting loci differently. Testing these estimates on populations for which data are available for several disorders shows good congruence with demographic data in some cases whereas in others higher growth rates are obtained, which may be the result of selection or hidden demographic processes.

Pollen flow in the wildservice tree,Sorbus torminalis(L.) Crantz. I. Evaluating the paternity analysis procedure in continuous populations

The joint development of polymorphic molecular markers and paternity analysis methods provides new approaches to investigate ongoing patterns of pollen flow in natural plant populations. However, paternity studies are hindered by false paternity assignment and the nondetection of true fathers. To gauge the risk of these two types of errors, we performed a simulation study to investigate the impact on paternity analysis of: (i) the assumed values for the size of the breeding male population (NBMP), and (ii) the rate of scoring error in genotype assessment. Our simulations were based on microsatellite data obtained from a natural population of the entomophilous wild service tree, Sorbus torminalis (L.) Crantz. We show that an accurate estimate of NBMP is required to minimize both types of errors, and we assess the reliability of a technique used to estimate NBMP based on parent-offspring genetic data. We then show that scoring errors in genotype assessment only slightly affect the assessment of paternity relationships, and conclude that it is generally better to neglect the scoring error rate in paternity analyses within a nonisolated population.

Modelling the impact of colonisation on genetic diversity and differentiation of forest trees: interaction of life cycle, pollen flow and seed long-distance dispersal

It was shown previously that the long lifespan and juvenile phase of trees strongly attenuate founder effects during colonisation in a diffusive dispersal model. However, this model yielded too slow a colonisation rate in comparison with palynological data for temperate forest trees. Since rare long-distance dispersal events have been shown to increase considerably colonisation rates in population dynamics models, we investigate here the impact of long-distance dispersal on within-population diversity (H(S)) and among-population differentiation (F(ST)) during the colonisation process. We use a stochastic approach and compare several dispersal strategies, ranging from very rare dispersal events of large amplitude to more frequent events of smaller amplitude. Using a simulation approach, which takes into account tree life-history traits, we show that long-distance dispersal events increase colonisation speed, and yield much larger founder effects in comparison with the diffusive model. The two models that include intermediate- and long-distance dispersal events show stronger deviations from experimental F(ST) values during and at the end of the colonisation process than the model with more frequent events of smaller dispersal variance. Furthermore, the introduction of a high level of pollen flow has a much more limited impact on models that include long-distance dispersal than on a diffusive dispersal model. The relatively high H(S) values that were obtained in all models are discussed according to the assumed mutation rate and effective population size. This study is an example of how observed genetic data can provide additional evidence on the best demographic model for a given species or group of species.

Pollen dispersal of tropical trees (Dinizia excelsa: Fabaceae) by native insects and African honeybees in pristine and fragmented Amazonian rainforest

Tropical rainforest trees typically occur in low population densities and rely on animals for cross-pollination. It is of conservation interest therefore to understand how rainforest fragmentation may alter the pollination and breeding structure of remnant trees. Previous studies of the Amazonian tree Dinizia excelsa (Fabaceae) found African honeybees (Apis mellifera scutellata) as the predominant pollinators of trees in highly disturbed habitats, transporting pollen up to 3.2 km between pasture trees. Here, using microsatellite genotypes of seed arrays, we compare outcrossing rates and pollen dispersal distances of (i) remnant D. excelsa in three large ranches, and (ii) a population in undisturbed forest in which African honeybees were absent. Self-fertilization was more frequent in the disturbed habitats (14%, n = 277 seeds from 12 mothers) than in undisturbed forest (10%, n = 295 seeds from 13 mothers). Pollen dispersal was extensive in all three ranches compared to undisturbed forest, however. Using a twogener analysis, we estimated a mean pollen dispersal distance of 1509 m in Colosso ranch, assuming an exponential dispersal function, and 212 m in undisturbed forest. The low effective density of D. excelsa in undisturbed forest (approximately 0.1 trees/ha) indicates that large areas of rainforest must be preserved to maintain minimum viable populations. Our results also suggest, however, that in highly disturbed habitats Apis mellifera may expand genetic neighbourhood areas, thereby linking fragmented and continuous forest populations.