What causes mating system shifts in plants? Arabidopsis lyrata as a case study.

The genetic breakdown of self-incompatibility (SI) and subsequent mating system shifts to inbreeding has intrigued evolutionary geneticists for decades. Most of our knowledge is derived from interspecific comparisons between inbreeding species and their outcrossing relatives, where inferences may be confounded by secondary mutations that arose after the initial loss of SI. Here, we study an intraspecific breakdown of SI and its consequences in North American Arabidopsis lyrata to test whether: (1) particular S-locus haplotypes are associated with the loss of SI and/or the shift to inbreeding; (2) a population bottleneck may have played a role in driving the transition to inbreeding; and (3) the mutation(s) underlying the loss of SI are likely to have occurred at the S-locus. Combining multiple approaches for genotyping, we found that outcrossing populations on average harbour 5 to 9 S-locus receptor kinase (SRK) alleles, but only two, S1 and S19, are shared by most inbreeding populations. Self-compatibility (SC) behaved genetically as a recessive trait, as expected from a loss-of-function mutation. Bulked segregant analysis in SC × SI F2 individuals using deep sequencing confirmed that all SC plants were S1 homozygotes but not all S1 homozygotes were SC. This was also revealed in population surveys, where only a few S1 homozygotes were SC. Together with crossing data, this suggests that there is a recessive factor that causes SC that is physically unlinked to the S-locus. Overall, our results emphasise the value of combining classical genetics with advanced sequencing approaches to resolve long outstanding questions in evolutionary biology.

Local adaptation with high gene flow: temperature parameters drive adaptation to altitude in the common frog (Rana temporaria).

Both environmental and genetic influences can result in phenotypic variation. Quantifying the relative contributions of local adaptation and phenotypic plasticity to phenotypes is key to understanding the effect of environmental variation on populations. Identifying the selective pressures that drive divergence is an important, but often lacking, next step. High gene flow between high- and low-altitude common frog (Rana temporaria) breeding sites has previously been demonstrated in Scotland. The aim of this study was to assess whether local adaptation occurs in the face of high gene flow and to identify potential environmental selection pressures that drive adaptation. Phenotypic variation in larval traits was quantified in R. temporaria from paired high- and low-altitude sites using three common temperature treatments. Local adaptation was assessed using Q(ST)-F(ST) analyses, and quantitative phenotypic divergence was related to environmental parameters using Mantel tests. Although evidence of local adaptation was found for all traits measured, only variation in larval period and growth rate was consistent with adaptation to altitude. Moreover, this was only evident in the three mountains with the highest high-altitude sites. This variation was correlated with mean summer and winter temperatures, suggesting that temperature parameters are potentially strong selective pressures maintaining local adaptation, despite high gene flow.

Using genetic variation to infer associations with climate in the common frog, Rana temporaria.

Recent and historical species' associations with climate can be inferred using molecular markers. This knowledge of population and species-level responses to climatic variables can then be used to predict the potential consequences of ongoing climate change. The aim of this study was to predict responses of Rana temporaria to environmental change in Scotland by inferring historical and contemporary patterns of gene flow in relation to current variation in local thermal conditions. We first inferred colonization patterns within Europe following the last glacial maximum by combining new and previously published mitochondrial DNA sequences. We found that sequences from our Scottish samples were identical to (92%), or clustered with, the common haplotype previously identified from Western Europe. This clade showed very low mitochondrial variation, which did not allow inference of historical colonization routes but did allow interpretation of patterns of current fine-scale population structure without consideration of confounding historical variation. Second, we assessed fine-scale microsatellite-based patterns of genetic variation in relation to current altitudinal temperature gradients. No population structure was found within altitudinal gradients (average FST=0.02), despite a mean annual temperature difference of 4.5 °C between low- and high-altitude sites. Levels of genetic diversity were considerable and did not vary between sites. The panmictic population structure observed, even along temperature gradients, is a potentially positive sign for R. temporaria persistence in Scotland in the face of a changing climate. This study demonstrates that within taxonomic groups, thought to be at high risk from environmental change, levels of vulnerability can vary, even within species.

Inbreeding depression in self-incompatible North-American Arabidopsis lyrata: disentangling genomic and S-locus-specific genetic load.

Newly formed selfing lineages may express recessive genetic load and suffer inbreeding depression. This can have a genome-wide genetic basis, or be due to loci linked to genes under balancing selection. Understanding the genetic architecture of inbreeding depression is important in the context of the maintenance of self-incompatibility and understanding the evolutionary dynamics of S-alleles. We addressed this using North-American subspecies of Arabidopsis lyrata. This species is normally self-incompatible and outcrossing, but some populations have undergone a transition to selfing. The goals of this study were to: (1) quantify the strength of inbreeding depression in North-American populations of A. lyrata; and (2) disentangle the relative contribution of S-linked genetic load compared with overall inbreeding depression. We enforced selfing in self-incompatible plants with known S-locus genotype by treatment with CO(2), and compared the performance of selfed vs outcrossed progeny. We found significant inbreeding depression for germination rate (δ=0.33), survival rate to 4 weeks (δ=0.45) and early growth (δ=0.07), but not for flowering rate. For two out of four S-alleles in our design, we detected significant S-linked load reflected by an under-representation of S-locus homozygotes in selfed progeny. The presence or absence of S-linked load could not be explained by the dominance level of S-alleles. Instead, the random nature of the mutation process may explain differences in the recessive deleterious load among lineages.

The effect of mating system on growth of Arabidopsis lyrata in response to inoculation with the biotrophic parasite Albugo candida.

The effects of variation in host reproductive systems on response to pathogens are not well understood. We inoculated individuals from outcrossing and inbreeding populations of North American Arabidopsis lyrata with Albugo candida (white blister rust) to test the effect of mating system and heterozygosity on disease response. We observed three host infection phenotypes, classified as fully resistant, partially resistant and fully susceptible. Overall, inbreeding populations had more susceptible and fewer partially resistant individuals than outcrossing populations, but the highest proportion of resistant individuals was found in two of the inbreeding populations. Mating system did not affect relative growth rate of inoculated plants, but there were strong effects of population and infection phenotype. We conclude that mating system per se does not determine the resistance of natural A. lyrata populations to infection by Albugo, but that the increased variability in responses among inbreeding populations may be due to reduced effective population size.

Sporophytic self-incompatibility genes and mating system variation in Arabis alpina.

BACKGROUND AND AIMS: Sporophytic self-incompatibility (SI) prevents inbreeding in many members of the Brassicaceae, and has been well documented in a variety of high-profile species. Arabis alpina is currently being developed as a model system for studying the ecological genetics of arctic-alpine environments, and is the focus of numerous studies on population structure and alpine phylogeography. Although it is highly inbreeding throughout most of its range, populations in central Italy have been identified that show inbreeding coefficients (F(IS)) more typical of self-incompatible relatives. The purpose of this study was to establish whether this variation is due to a functioning SI system. METHODS: Outcrossing rate estimates were calculated based on 16 allozyme loci and self-compatibility assessed based on controlled pollinations for six Italian populations that have previously been shown to vary in F(IS) values. Putative SRK alleles (the gene controlling the female component of SI in other Brassicaceae) amplified from A. alpina were compared with those published for other species. Linkage of putative SRK alleles and SI phenotypes was assessed using a diallel cross. KEY RESULTS: Functional avoidance of inbreeding is demonstrated in three populations of A. alpina, corresponding with previous F(IS) values. The presence is described of 15 putative SRK-like alleles, which show high sequence identity to known alleles from Brassica and Arabidopsis and the high levels of synonymous and nonsynonymous variation typical of genes under balancing selection. Also, orthologues of two other members of the S-receptor kinase gene family, Aly8 (ARK3) and Aly9 (AtS1) are identified. Further to this, co-segregation between some of the putative S-alleles and compatibility phenotypes was demonstrated using a full-sibling cross design. CONCLUSIONS: The results strongly suggest that, as with other species in the Brassicaceae, A. alpina has a sporophytic SI system but shows variation in the strength of SI within and between populations.

Bayesian paternity analysis and mating patterns in a parasitic nematode, Trichostrongylus tenuis.

Mating behaviour is a fundamental aspect of the evolutionary ecology of sexually reproducing species, but one that has been under-researched in parasitic nematodes. We analysed mating behaviour in the parasitic nematode Trichostrongylus tenuis by performing a paternity analysis in a population from a single red grouse host. Paternity of the 150 larval offspring of 25 mothers (sampled from one of the two host caeca) was assigned among 294 candidate fathers (sampled from both caeca). Each candidate father's probability of paternity of each offspring was estimated from 10-locus microsatellite genotypes. Seventy-six (51%) offspring were assigned a father with a probability of >0.8, and the estimated number of unsampled males was 136 (95% credible interval (CI) 77-219). The probability of a male from one caecum fathering an offspring in the other caecum was estimated as 0.024 (95% CI 0.003-0.077), indicating that the junction of the caeca is a strong barrier to dispersal. Levels of promiscuity (defined as the probability of two of an adult's offspring sharing only one parent) were high for both sexes. Variance in male reproductive success was moderately high, possibly because of a combination of random mating and high variance in post-copulatory reproductive success. These results provide the first data on individual mating behaviour among parasitic nematodes.

Multiple losses of self-incompatibility in North-American Arabidopsis lyrata?: phylogeographic context and population genetic consequences.

Arabidopsis lyrata is mostly outcrossing due to a sporophytic self-incompatibility (SI) system but around the Great Lakes of North America some populations have experienced a loss of SI. We researched the loss of SI in a phylogeographic context. We used cpDNA and microsatellite markers to test if populations of North-American A. lyrata around the Great Lakes have experienced different (recent) histories, and linked this with individually established selfing phenotype and population level realized outcrossing rates calculated based on variation in progeny arrays at multi-locus microsatellite markers. We found three chloroplast haplotypes, in two of which the loss of self-incompatibility had occurred independently. Shifts to high rates of inbreeding were most apparent in one of these lineages but individuals showing loss of SI occurred in all three. Self-compatible individuals usually showed a reduction of observed heterozygosity (H(O)) compared to outcrossing individuals. In the lineage that included the populations with the highest levels of inbreeding, this reduction was more substantial. This may indicate that the loss of SI in this lineage did not occur as recently as in the other lineage. The geographic distribution of the haplotypes suggested that there had been at least two independent colonization routes to the north of the Great Lakes following the last glaciation. This is consistent with postglacial migration patterns that have been suggested for other organisms with limited dispersal, such as reptiles and amphibians.

Identification and characterization of a polymorphic receptor kinase gene linked to the self-incompatibility locus of Arabidopsis lyrata.

We study the segregation of variants of a putative self-incompatibility gene in Arabidopsis lyrata. This gene encodes a sequence that is homologous to the protein encoded by the SRK gene involved in self-incompatibility in Brassica species. We show by diallel pollinations of plants in several full-sib families that seven different sequences of the gene in A. lyrata are linked to different S-alleles, and segregation analysis in further sibships shows that four other sequences behave as allelic to these. The family data on incompatibility provide evidence for dominance classes among the S-alleles, as expected for a sporophytic SI system. We observe no division into pollen-dominant and pollen-recessive classes of alleles as has been found in Brassica, but our alleles fall into at least three dominance classes in both pollen and stigma expression. The diversity among sequences of the A. lyrata putative S-alleles is greater than among the published Brassica SRK sequences, and, unlike Brassica, the alleles do not cluster into groups with similar dominance.

Ploidy evolution in the yeast Saccharomyces cerevisiae: a test of the nutrient limitation hypothesis.

The nutrient limitation hypothesis provides a nongenetic explanation for the evolution of life cycles that retain both haploid and diploid phases: differences in nutrient requirements and uptake allow haploids to override the potential genetic advantages provided by diploidy under certain nutrient limiting conditions. The relative fitness of an isogenic series of haploid, diploid and tetraploid yeast cells (Saccharomyces cerevisiae), which were also equivalent at the mating type locus, was measured. Fitness was measured both by growth rate against a common competitor and by intrinsic growth rate in isolated cultures, under four environmental conditions: (1) rich medium (YPD) at the preferred growth temperature (30 °C); (2) nutrient poor medium (MM) at 30 °C; (3) YPD at a nonpreferred temperature (37 °C); and (4) MM at 37 °C. In contrast to the predictions of the nutrient limitation hypothesis, haploids grew significantly faster than diploids under nutrient rich conditions, but there were no apparent differences between them when fitness was determined by relative competitive ability. In addition, temperature affected the relative growth of haploids and diploids, with haploids growing proportionately faster at higher temperatures. Tetraploids performed very poorly under all conditions compared. Cell geometric parameters were not consistent predictors of fitness under the conditions measured.

Validation of swabs as a non-destructive and relatively non-invasive DNA sampling method in fish.

Non-destructive methods of collecting DNA from small fish species can be problematic, as fin clips can potentially affect behaviour or survivorship in the wild. Swabbing body mucus may provide a less invasive method of DNA collection. However, risk of contamination from other individuals in high density groups could give erroneous genotyping results. We compared multilocus microsatellite genotypes from the same individuals when collected at low and high density and compared this with fin clips. We found no differences between these categories, with a genotyping error rate of 0.42%, validating the use of body mucus swabbing for DNA collection in fish.

Masking and purging mutations following EMS treatment in haploid, diploid and tetraploid yeast (Saccharomyces cerevisiae).

The yeast, Saccharomyces cerevisiae, was used as a model to investigate theories of ploidy evolution. Mutagenesis experiments using the alkylating agent EMS (ethane methyl sulphonate) were conducted to assess the relative importance that masking of deleterious mutations has on response to and recovery from DNA damage. In particular, we tested whether cells with higher ploidy levels have relatively higher fitnesses after mutagenesis, whether the advantages of masking are more pronounced in tetraploids than in diploids, and whether purging of mutations allows more rapid recovery of haploid cells than cells with higher ploidy levels. Separate experiments were performed on asexually propagating stationary phase cells using (1) prototrophic haploid (MAT alpha) and diploid (MATa/alpha) strains and (2) isogenic haploid, diploid and tetraploid strains lacking a functional mating type locus. In both sets of experiments, haploids showed a more pronounced decrease in apparent growth rate than diploids, but both haploids and diploids appeared to recover very rapidly. Tetraploids did not show increased benefits of masking compared with diploids but volume measurements and FACScan analyses on the auxotrophic strains indicated that all treated tetraploid strains decreased in ploidy level and that some of the treated haploid lines increased in ploidy level. Results from these experiments confirm that while masking deleterious mutations provides an immediate advantage to higher ploidy levels in the presence of mutagens, selection is extremely efficient at removing induced mutations, leading growth rates to increase rapidly over time at all ploidy levels. Furthermore, ploidy level is itself a mutable trait in the presence of EMS, with both haploids and tetraploids often evolving towards diploidy (the ancestral state of S. cerevisiae) during the course of the experiment.

Hybridization between tetraploid and diploid species of treefrogs (Genus hyla).

An indirect analysis of female meiotic mechanisms in the tetraploid (4n = 48) grey treefrog, Hyla versicolor, was performed by examining artificially produced hybrids. First generation hybrids between a H. versicolor female from Canada (4n = 48) and a H. arborea male from France (2n = 24) were all triploid and appeared to contain two sets of chromosomes from H. versicolor and one set from H. arborea. Males and females were produced in equal numbers but testes in general were more completely developed than ovaries. Electrophoretic analysis of selected allozyme loci suggested that gene products from the two parents were not equally expressed. Backcross hybrids were triploid, tetraploid, and pentaploid but did not appear to contain recognizable chromosomes from the H. arborea grandparent. Allozymes from these hybrids indicated that only H. versicolor alleles were expressed, as none of the distinctive H. arborea alleles present in the triploid male parent were present in the offspring. It was concluded that preferential pairing of chromosomes and gene regulatory biases may help to explain factors that relate to the ability of tetraploids to hybridize with even distantly related taxa and may be involved in the rediploidization process that usually follows polyploidization.

Estimating the number, frequency, and dominance of S-alleles in a natural population of Arabidopsis lyrata(Brassicaceae) with sporophytic control of self-incompatibility.

In homomorphic plant self-incompatibility (SI) systems, large numbers of alleles may be maintained at a single Mendelian locus. Most estimators of the number of alleles present in natural populations are designed for gametophytic self-incompatibility systems (GSI) in which the recognition phenotype of the pollen is determined by its own haploid genotype. In sporophytic systems (SSI), the recognition phenotype of the pollen is determined by the diploid genotype of its parent, and dominance differs among alleles. We describe research aimed at estimates of S-allele numbers in a natural population of Arabidopsis lyrata (Brassicaceae), whose SSI system has recently been described. Using a combination of pollination studies and PCR-based identification of alleles at a locus equivalent to the Brassica SRK gene, we identified and sequenced 11 putative alleles in a sample of 20 individuals from different maternal seed sets. The pollination results indicate that we have not amplified all alleles that must be present. Extensive partial incompatibility, nonreciprocal compatibility differences, and evidence of weakened expression of SI in some genotypes, prevent us from determining the exact number of missing alleles based only on cross-pollination data. Although we show that none of the theoretical models currently proposed is completely appropriate for estimating the number of alleles in this system, we estimate that there are between 13 and 16 different S-alleles in our sample, probably between 16 and 25 alleles in the population, and discuss the relative frequency of alleles in relation to dominance.

Phylogenetic analysis of reptilian hemoglobins: trees, rates, and divergences.

Phylogenetic relationships among reptiles were examined using previously published and newly determined hemoglobin sequences. Trees reconstructed from these sequences using maximum-parsimony, neighbor-joining, and maximum-likelihood algorithms were compared with a phylogenetic tree of Amniota, which was assembled on the basis of published morphological data. All analyses differentiated alpha chains into alphaA and alphaD types, which are present in all reptiles except crocodiles, where only alphaA chains are expressed. The occurrence of the alphaD chain in squamates (lizards and snakes only in this study) appears to be a general characteristic of these species. Lizards and snakes also express two types of beta chains (betaI and betaII), while only one type of beta chain is present in birds and crocodiles. Reconstructed hemoglobin trees for both alpha and beta sequences did not yield the monophyletic Archosauria (i.e., crocodilians + birds) and Lepidosauria (i.e., Sphenodon + squamates) groups defined by the morphology tree. This discrepancy, as well as some other poorly resolved nodes, might be due to substantial heterogeneity in evolutionary rates among single hemoglobin lineages. Estimation of branch lengths based on uncorrected amino acid substitutions and on distances corrected for multiple substitutions (PAM distances) revealed that relative rates for squamate alphaA and alphaD chains and crocodilian beta chains are at least twice as high as those of the rest of the chains considered. In contrast to these rate inequalities between reptilian orders, little variation was found within squamates, which allowed determination of absolute evolutionary rates for this subset of hemoglobins. Rate estimates for hemoglobins of lizards and snakes yielded 1.7 (alphaA) and 3.3 (beta) million years/PAM when calibrated with published divergence time vs. PAM distance correlates for several speciation events within snakes and for the squamate left and right arrow sphenodontid split. This suggests that hemoglobin chains of squamate reptiles evolved approximately 3.5 (alphaA) or approximately 1.7 times (beta) faster than their mammalian equivalents. These data also were used to obtain a first estimate of some intrasquamate divergence times.

Inheritance and dominance of self-incompatibility alleles in polyploid Arabidopsis lyrata.

Natural populations of diploid Arabidopsis lyrata exhibit the sporophytic type of self-incompatibility system characteristic of Brassicaceae, in which complicated dominance interactions among alleles in the diploid parent determine self-recognition phenotypes of both pollen and stigma. The purpose of this study was to investigate how polyploidy affects this already complex system. One tetraploid population (Arabidopsis lyrata ssp kawasakiana from Japan) showed complete self-compatibility and produced viable selfed progeny for at least three generations subsequent to field collection. In contrast, individuals from a second tetraploid population (A. lyrata ssp petraea from Austria) were strongly self-incompatible (SI). Segregation of SI genotypes in this population followed Mendelian patterns based on a tetrasomic model of inheritance, with two to four alleles per individual, independent segregation of alleles, and little evidence of dosage effects of alleles found in multiple copies. Similar to results from diploids, anomalous compatibility patterns involving particular combinations of individuals occurred at a low frequency in the tetraploids, suggesting altered dominance in certain genetic backgrounds that could be due to the influence of a modifier locus. Overall, dominance relationships among S-alleles in self-incompatible tetraploid families were remarkably similar to those in related diploids, suggesting that this very important and complicated locus has not undergone extensive modification subsequent to polyploidization.