Effect of balancing selection on spatial genetic structure within populations: theoretical investigations on the self-incompatibility locus and empirical studies in Arabidopsis halleri

Zygosity-based sex determination in a butterfly drives hypervariability of Masculinizer

Nature has devised many ways of producing males and females. Here, we report on a previously undescribed mechanism for Lepidoptera that functions without a female-specific gene. The number of alleles or allele heterozygosity in a single Z-linked gene (BaMasc) is the primary sex-determining switch in Bicyclus anynana butterflies. Embryos carrying a single BaMasc allele develop into WZ (or Z0) females, those carrying two distinct alleles develop into ZZ males, while (ZZ) homozygotes initiate female development, have mismatched dosage compensation, and die as embryos. Consequently, selection against homozygotes has favored the evolution of spectacular allelic diversity: 205 different coding sequences of BaMasc were detected in a sample of 246 females. The structural similarity of a hypervariable region (HVR) in BaMasc to the HVR in Apis mellifera csd suggests molecular convergence between deeply diverged insect lineages. Our discovery of this primary switch highlights the fascinating diversity of sex-determining mechanisms and underlying evolutionary drivers.

Origin and persistence of polymorphism in loci targeted by disassortative preference: a general model

The emergence and persistence of polymorphism within populations generally requires specific regimes of natural or sexual selection. Here, we develop a unified theoretical framework to explore how polymorphism at targeted loci can be generated and maintained by either disassortative mating choice or balancing selection due to, for example, heterozygote advantage. To this aim, we model the dynamics of alleles at a single locus A in a population of haploid individuals, where reproductive success depends on the combination of alleles carried by the parents at locus A. Our theoretical study of the model confirms that the conditions for the persistence of a given level of allelic polymorphism depend on the relative reproductive advantages among pairs of individuals. Interestingly, equilibria with unbalanced allelic frequencies were shown to emerge from successive introduction of mutants. We then investigate the role of the function linking allelic divergence to reproductive advantage on the evolutionary fate of alleles within the population. Our results highlight the significance of the shape of this function for both the number of alleles maintained and their level of genetic divergence. Large number of alleles are maintained with substantial replacement of alleles, when disassortative advantage slowly increases with allelic differentiation . In contrast, few highly differentiated alleles are predicted to be maintained when genetic differentiation has a strong effect on disassortative advantage. These opposite effects predicted by our model explain how disassortative mate choice may lead to various levels of allelic differentiation and polymorphism, and shed light on the effect of mate preferences on the persistence of balanced and unbalanced polymorphism in natural population.

Natural Populations of Astrocaryum aculeatum Meyer in Amazonia: Genetic Diversity and Conservation

Astrocaryum aculeatum, a palm tree incipiently domesticated from upland ecosystems in the Brazilian Amazon, is especially adapted to anthropized areas. The pulp of the fruit, obtained by extractivism, is consumed fresh by the Amazonian population. The objective of the study is to evaluate the diversity and genetic structure of the natural populations of A. aculeatum, exploited by extractive farmers in Amazonas, Brazil, seeking to suggest conservation and management strategies for this species. A total of 218 plants were sampled in 15 populations in 14 municipalities in the state of Amazonas, evaluated by 12 microsatellite loci. A total of 101 alleles were observed. The means of the observed heterozygosities (H_O = 0.6390) were higher than expected (H_E = 0.557), with high levels of heterozygotes in the populations. The fixation index in the loci and populations was negative. The F_ST (0.07) and AMOVA showed moderate population structure. Bayesian analysis indicated the grouping k = 4 as the most adequate. There is a high genetic diversity in populations, with a moderate genetic structure due to possible historical events, which could be related to the process of subpopulation formation, possibly presenting three historical moments: before and after the beginning of deforestation and today. The conservation and management policies of this species must be carried out at a watershed level.

Infection phenotypes of a coevolving parasite are highly diverse, structured, and specific

Understanding how diversity is maintained in natural populations is a major goal of evolutionary biology. In coevolving hosts and parasites, negative frequency-dependent selection is one mechanism predicted to maintain genetic variation. While much is known about host diversity, parasite diversity remains understudied in coevolutionary research. Here, we survey natural diversity in a bacterial parasite by characterizing infection phenotypes for over 50 isolates in relation to 12 genotypes of their host, Daphnia magna. We find striking phenotypic variation among parasite isolates, and we discover the parasite can infect its host through at least five different attachment sites. Variation in attachment success at each site is explained to varying degrees by host and parasite genotypes. A spatial correlation analysis showed that infectivity of different isolates does not correlate with geographic distance, meaning isolates from widespread populations are equally able to infect the host. Overall, our results reveal that infection phenotypes of this parasite are highly diverse. Our results are consistent with the prediction that under Red Queen coevolutionary dynamics both the host and the parasite should show high genetic diversity for traits of functional importance in their interactions.

Genome-Wide Association Analysis Identifies a Genetic Basis of Infectivity in a Model Bacterial Pathogen

Knowledge of the genetic architecture of pathogen infectivity and host resistance is essential for a mechanistic understanding of coevolutionary processes, yet the genetic basis of these interacting traits remains unknown for most host-pathogen systems. We used a comparative genomic approach to explore the genetic basis of infectivity in Pasteuria ramosa, a Gram-positive bacterial pathogen of planktonic crustaceans that has been established as a model for studies of Red Queen host-pathogen coevolution. We sequenced the genomes of a geographically, phenotypically, and genetically diverse collection of P. ramosa strains and performed a genome-wide association study to identify genetic correlates of infection phenotype. We found multiple polymorphisms within a single gene, Pcl7, that correlate perfectly with one common and widespread infection phenotype. We then confirmed this perfect association via Sanger sequencing in a large and diverse sample set of P. ramosa clones. Pcl7 codes for a collagen-like protein, a class of adhesion proteins known or suspected to be involved in the infection mechanisms of a number of important bacterial pathogens. Consistent with expectations under Red Queen coevolution, sequence variation of Pcl7 shows evidence of balancing selection, including extraordinarily high diversity and absence of geographic structure. Based on structural homology with a collagen-like protein of Bacillus anthracis, we propose a hypothesis for the structure of Pcl7 and the physical location of the phenotype-associated polymorphisms. Our results offer strong evidence for a gene governing infectivity and provide a molecular basis for further study of Red Queen dynamics in this model host-pathogen system.

Genetic Structure in Populations of Euterpe precatoria Mart. in the Brazilian Amazon

Euterpe precatoria is a palm tree belonging to the Arecaceae family, occurring in Western and Central Brazilian Amazonia. Its fruit, which is very appreciated in the Amazon region, produces pulp that is consumed in fresh form. Its production is carried out almost exclusively by extractive farmers. In order to establish adequate strategies to sustain this genetic resource, we need knowledge about the diversity and genetic structure in natural populations. This study aimed to evaluate the influence of geographic distance on genetic structure in the main extractive populations of E. precatoria in the Brazilian Amazon. Leaves from 377 plants were collected in 19 populations located in 16 municipalities in the State of Amazonas and three in the State of Rondônia. Twelve microsatellite loci were used to genotype the plants. The diversity and genetic structure among populations were estimated. The average number of alleles per locus was 5.97. The observed heterozygosity means (HO) were higher than expected (HE) at the population level (HO = 0.72, HE = 0.66) and fixation index (f = -0.100) was negative. The FST value (0.1820) and the AMOVA results (Φ = 0.1796) showed population structure. The populations were clustered into three groups (K = 3) in the Bayesian analysis. The Discriminant Analysis of Principal Components (DAPC) confirmed eight clusters, with the populations close to those identified by the Bayesian analysis. The geographic differentiation was confirmed by the groupings obtained in the Structure analysis and the DACP function. Information related to phenotypic, genetic and environmental characterization of populations is important to guide conservation and management strategies and the formulation of public species management policies in Amazonia.

Host-parasite co-evolution and its genomic signature

Studies in diverse biological systems have indicated that host-parasite co-evolution is responsible for the extraordinary genetic diversity seen in some genomic regions, such as major histocompatibility (MHC) genes in jawed vertebrates and resistance genes in plants. This diversity is believed to evolve under balancing selection on hosts by parasites. However, the mechanisms that link the genomic signatures in these regions to the underlying co-evolutionary process are only slowly emerging. We still lack a clear picture of the co-evolutionary concepts and of the genetic basis of the co-evolving phenotypic traits in the interacting antagonists. Emerging genomic tools that provide new options for identifying underlying genes will contribute to a fuller understanding of the co-evolutionary process.

Paternity analysis, pollen flow, and spatial genetic structure of a natural population of Euterpe precatoria in the Brazilian Amazon

Euterpe precatoria, known as açaí do Amazonas, is a regionally important palm of the Amazon rainforest for the fruit production through extractive agriculture. Little information is available with regard to genetic diversity, gene flow, and spatial genetic structure (SGS) of açaí populations, which are essential for the use, management, and conservation of genetic resources of the species. This research aimed to assess the genetic diversity, inbreeding level, SGS, and gene flow in four ontogenetic stages of a natural E. precatoria population in the Brazilian Amazon, based on 18 microsatellite loci. The study was carried out in a natural population dispersed in an area of about 10 ha. Leaf tissues of 248 plants were mapped and sampled and classified into four ontogenetic stages: reproductive (59), immature (70), young (60), and seedling (59). Genetic diversity indices were high for all ontogenetic stages. The fixation index (F) for all ontogenetic stages was not significantly different from zero, indicating the absence of inbreeding. A significant SGS was found for all ontogenetic stages (68-110 m), indicating seed dispersal over short distances. Paternity analysis detected pollen immigration of 39.1%, a selfing rate of 4.2%, and a mean pollen dispersal distance within the population of 531 m. The results indicate substantial allele input in the population via pollen immigration, contributing to the maintenance of the genetic diversity of the population. However, within a population, the renewal with new progenies selected from seed plants spaced at least 110 m apart is important to avoid collecting seeds from related plants.

Balancing selection and introgression of newt immune-response genes

The importance of interspecific introgression as a source of adaptive variation is increasingly recognized. Theory predicts that beneficial genetic variants cross species boundaries easily even when interspecific hybridization is rare and gene flow is strongly constrained throughout the genome. However, it remains unclear whether certain classes of genes are particularly prone to adaptive introgression. Genes affected by balancing selection (BS) may constitute such a class, because forms of BS that favour novel, initially rare alleles, should facilitate introgression. We tested this hypothesis in hybridizing newts by comparing 13 genes with signatures of BS, in particular an excess of common non-synonymous polymorphisms, to the genomic background (154 genes). Parapatric hybridizing taxa were less differentiated in BS candidate genes than more closely related allopatric lineages, while the opposite was observed in the control genes. Coalescent and forward simulations that explored neutral and BS scenarios under isolation and migration showed that processes other than differential gene flow are unlikely to account for this pattern. We conclude that BS, probably involving a form of novel allele advantage, promotes introgression. This mechanism may be a source of adaptively relevant variation in hybridizing species over prolonged periods.

Genetic architecture and balancing selection: the life and death of differentiated variants

Balancing selection describes any form of natural selection, which results in the persistence of multiple variants of a trait at intermediate frequencies within populations. By offering up a snapshot of multiple co-occurring functional variants and their interactions, systems under balancing selection can reveal the evolutionary mechanisms favouring the emergence and persistence of adaptive variation in natural populations. We here focus on the mechanisms by which several functional variants for a given trait can arise, a process typically requiring multiple epistatic mutations. We highlight how balancing selection can favour specific features in the genetic architecture and review the evolutionary and molecular mechanisms shaping this architecture. First, balancing selection affects the number of loci underlying differentiated traits and their respective effects. Control by one or few loci favours the persistence of differentiated functional variants by limiting intergenic recombination, or its impact, and may sometimes lead to the evolution of supergenes. Chromosomal rearrangements, particularly inversions, preventing adaptive combinations from being dissociated are increasingly being noted as features of such systems. Similarly, due to the frequency of heterozygotes maintained by balancing selection, dominance may be a key property of adaptive variants. High heterozygosity and limited recombination also influence associated genetic load, as linked recessive deleterious mutations may be sheltered. The capture of deleterious elements in a locus under balancing selection may reinforce polymorphism by further promoting heterozygotes. Finally, according to recent genomewide scans, balanced polymorphism might be more pervasive than generally thought. We stress the need for both functional and ecological studies to characterize the evolutionary mechanisms operating in these systems.

Population genetic structure of cotton pink bollworm, Pectinophora gossypiella (Saunders) (Lepidoptera: Gelechiidae) using mitochondrial cytochrome oxidase I (COI) gene sequences from India

Pink bollworm (PBW), Pectinophora gossypiella is one of the most destructive pest's globally inflicting huge economic losses in cotton even during later stages of crop growth. In the present investigation, the population genetic structure, distribution, and genetic diversity of P. gossypiella in cotton growing zones of India using partial mitochondrial DNA cytochrome oxidase-I (COI) gene was addressed. The overall haplotype (Hd), number of nucleotide differences (K), and nucleotide diversity (π) were 0.3028, 0.327, and 0.00047, respectively which suggest that entire population exhibited low level of genetic diversity. Zone-wise clustering of population revealed that central zone recorded low level of Hd (0.2730) as compared to north (0.3619) and south (0.3028) zones. The most common haplotype (H1) reported in all 19 locations could be proposed as ancestral/original haplotype. This haplotype with one mutational step formed star-like phylogeny connected with 11 other haplotypes. The phylogenetic relationship studies revealed that most haplotypes of populations are closely related to each other. Haplotype 5 was exclusively present in Dharwad (South zone) shared with populations of Hanumangarh and Bathinda (North zone). The result indicated that there is no isolation by distance effect among the Indian populations of PBW. The present study reports a low genetic diversity among PBW populations of India and H1, as ancestral haplotype from which other haplotypes have evolved suggests that the migration and dispersal over long distance and invasiveness are major factors.

The evolution of selfing from outcrossing ancestors in Brassicaceae: what have we learned from variation at the S-locus?

Evolutionary transitions between mating systems have occurred repetitively and independently in flowering plants. One of the most spectacular advances of the recent empirical literature in the field was the discovery of the underlying genetic machinery, which provides the opportunity to retrospectively document the scenario of the outcrossing to selfing transitions in a phylogenetic perspective. In this review, we explore the literature describing patterns of polymorphism and molecular evolution of the locus controlling self-incompatibility (S-locus) in selfing species of the Brassicaceae family in order to document the transition from outcrossing to selfing, a retrospective approach that we describe as the 'mating system genes approach'. The data point to strikingly contrasted scenarios of transition from outcrossing to selfing. We also perform original analyses of the fully sequenced genomes of four species showing self-compatibility, to compare the orthologous S-locus region with that of functional S-locus haplotypes. Phylogenetic analyses suggest that all species we investigated evolved independently towards loss of self-incompatibility, and in most cases almost intact sequences of either of the two S-locus genes suggest that these transitions occurred relatively recently. The S-locus region in Aethionema arabicum, representing the most basal lineage of Brassicaceae, showed unusual patterns so that our analysis could not determine whether self-incompatibility was lost secondarily, or evolved in the core Brassicaceae after the split with this basal lineage. Although the approach we detail can only be used when mating system genes have been identified in a clade, we suggest that its integration with phylogenetic and population genetic approaches should help determine the main routes of this predominant mating system shift in plants.

Pollen limitation and reduced reproductive success are associated with local genetic effects in Prunus virginiana, a widely distributed self-incompatible shrub

BACKGROUND AND AIMS

A vast quantity of empirical evidence suggests that insufficient quantity or quality of pollen may lead to a reduction in fruit set, in particular for self-incompatible species. This study uses an integrative approach that combines field research with marker gene analysis to understand the factors affecting reproductive success in a widely distributed self-incompatible species, Prunus virginiana (Rosaceae).

METHODS

Twelve patches of P. virginiana distributed within three populations that differed in degree of disturbance were examined. Two of the sites were small (7-35 km(2)) remnants of forest in an intensively used agricultural landscape, while the third was continuous (350 km(2)) and less disturbed. Field studies (natural and hand cross-pollinations) were combined with marker gene analyses (microsatellites and S-locus) in order to explore potential factors affecting pollen delivery and consequently reproductive success at landscape (between populations) and fine scales (within populations).

KEY RESULTS

Reductions in reproductive output were found in the two fragments compared with the continuous population, and suggest that pollen is an important factor limiting fruit production. Genetic analyses carried out in one of the fragments and in the continuous site suggest that even though S-allele diversity is high in both populations, the fragment exhibits an increase in biparental inbreeding and correlated paternity. The increase in biparental inbreeding in the fragment is potentially attributable to variation in the density of individuals and/or the spatial distribution of genotypes among populations, both of which could alter mating dynamics.

CONCLUSIONS

By using a novel integrative approach, this study shows that even though P. virginiana is a widespread species, fragmented populations can experience significant reductions in fruit set and pollen limitation in the field. Deatiled examination of one fragmented population suggests that these linitations may be explained by an increase in biparental inbreeding, correlated paternity and fine-scale genetic structure. The consistency of the field and fine-scale genetic analyses, and the consistency of the results within patches and across years, suggest that these are important processes driving pollen limitation in the fragment.

The population genetics of sporophytic self-incompatibility in three hybridizing senecio (asteraceae) species with contrasting population histories

Hybridization generates evolutionary novelty and spreads adaptive variation. By promoting outcrossing, plant self-incompatibility (SI) systems also favor interspecific hybridization because the S locus is under strong negative frequency-dependent balancing selection. This study investigates the SI mating systems of three hybridizing Senecio species with contrasting population histories. Senecio aethnensis and S. chrysanthemifolius native to Sicily, form a hybrid zone at intermediate altitudes on Mount Etna, and their neo-homoploid hybrid species, S. squalidus, has colonized disturbed urban habitats in the UK during the last 150 years. We show that all three species express sporophytic SI (SSI), where pollen incompatibility is controlled by the diploid parental genome, and that SSI is inherited and functions normally in hybrids. Large-scale crossing studies of wild sampled populations allowed direct comparison of SSI between species and found that the main impacts of colonization in S. squalidus compared to Sicilian Senecio was a reduced number of S alleles, increased S allele frequencies, and increased interpopulation S allele sharing. In general, many S alleles were shared between species and the S locus showed reduced intra- and interspecific population genetic structure compared to molecular genetic markers, indicative of enhanced effective gene flow due to balancing selection.

Impact of negative frequency-dependent selection on mating pattern and genetic structure: a comparative analysis of the S-locus and nuclear SSR loci in Prunus lannesiana var. speciosa

Mating processes of local demes and spatial genetic structure of island populations at the self-incompatibility (S-) locus under negative frequency-dependent selection (NFDS) were evaluated in Prunus lannesiana var. speciosa in comparison with nuclear simple sequence repeat (SSR) loci that seemed to be evolutionarily neutral. Our observations of local mating patterns indicated that male-female pair fecundity was influenced by not only self-incompatibility, but also various factors, such as kinship, pollen production and flowering synchrony. In spite of the mating bias caused by these factors, the NFDS effect on changes in allele frequencies from potential mates to mating pollen was detected at the S-locus but not at the SSR loci, although the changes from adult to juvenile cohorts were not apparent at any loci. Genetic differentiation and isolation-by-distance over various spatial scales were smaller at the S-locus than at the SSR loci, as expected under the NFDS. Allele-sharing distributions among the populations also had a unimodal pattern at the S-locus, indicating the NFDS effect except for alleles unique to individual populations probably due to isolation among islands, although this pattern was not exhibited by the SSR loci. Our results suggest that the NFDS at the S-locus has an impact on both the mating patterns and the genetic structure in the P. lannesiana populations studied.

Microevolution of S-allele frequencies in wild cherry populations: respective impacts of negative frequency dependent selection and genetic drift

Negative frequency dependent selection (NFDS) is supposed to be the main force controlling allele evolution at the gametophytic self-incompatibility locus (S-locus) in strictly outcrossing species. Genetic drift also influences S-allele evolution. In perennial sessile organisms, evolution of allelic frequencies over two generations is mainly shaped by individual fecundities and spatial processes. Using wild cherry populations between two successive generations, we tested whether S-alleles evolved following NFDS qualitative and quantitative predictions. We showed that allelic variation was negatively correlated with parental allelic frequency as expected under NFDS. However, NFDS predictions in finite population failed to predict more than half S-allele quantitative evolution. We developed a spatially explicit mating model that included the S-locus. We studied the effects of self-incompatibility and local drift within populations due to pollen dispersal in spatially distributed individuals, and variation in male fecundity on male mating success and allelic frequency evolution. Male mating success was negatively related to male allelic frequency as expected under NFDS. Spatial genetic structure combined with self-incompatibility resulted in higher effective pollen dispersal. Limited pollen dispersal in structured distributions of individuals and genotypes and unequal pollen production significantly contributed to S-allele frequency evolution by creating local drift effects strong enough to counteract the NFDS effect on some alleles.

Evolution of the S-locus region in Arabidopsis relatives

The S locus, a single polymorphic locus, is responsible for self-incompatibility (SI) in the Brassicaceae family and many related plant families. Despite its importance, our knowledge of S-locus evolution is largely restricted to the causal genes encoding the S-locus receptor kinase (SRK) receptor and S-locus cysteine-rich protein (SCR) ligand of the SI system. Here, we present high-quality sequences of the genomic region of six S-locus haplotypes: Arabidopsis (Arabidopsis thaliana; one haplotype), Arabidopsis lyrata (four haplotypes), and Capsella rubella (one haplotype). We compared these with reference S-locus haplotypes of the self-compatible Arabidopsis and its SI congener A. lyrata. We subsequently reconstructed the likely genomic organization of the S locus in the most recent common ancestor of Arabidopsis and Capsella. As previously reported, the two SI-determining genes, SCR and SRK, showed a pattern of coevolution. In addition, consistent with previous studies, we found that duplication, gene conversion, and positive selection have been important factors in the evolution of these two genes and appear to contribute to the generation of new recognition specificities. Intriguingly, the inactive pseudo-S-locus haplotype in the self-compatible species C. rubella is likely to be an old S-locus haplotype that only very recently became fixed when C. rubella split off from its SI ancestor, Capsella grandiflora.