Putting the genes into community genetics

As part of the long-term fusion of evolutionary biology and ecology (Ford, 1964), the field of community genetics has made tremendous progress in describing the impacts of plant genetic variation on community and ecosystem processes. In the "genes-to-ecosystems" framework (Whitham et al., 2003), genetically based traits of plant species have ecological consequences, but previous studies have not identified specific plant genes responsible for community phenotypes. The study by Barker et al. (2019) in this issue of Molecular Ecology uses an impressive common garden experiment of trembling aspen (Figure 1) to test for the genetic basis of tree traits that shape the insect community composition. Using a Genome-Wide Association Study (GWAS), they found that genomic regions associated with phytochemical traits best explain variation in herbivore community composition, and identified specific genes associated with different types of leaf-modifying herbivores and ants. This is one of the first studies to identify candidate genes underlying the heritable plant traits that explain patterns of insect biodiversity.

Identification and characterization of microsatellite loci in the tuliptree, Liriodendron tulipifera (Magnoliaceae)

PREMISE OF THE STUDY

Twenty-three polymorphic microsatellite loci (simple sequence repeats) were identified and characterized for Liriodendron tulipifera (Magnoliaceae), a species native to eastern North America, to investigate its genetic diversity, population structure, and mating system.

METHODS AND RESULTS

Using Illumina HiSeq paired-end reads from genomic DNA, searches for repeat motifs identified approximately 280,000 potentially amplifiable loci. Of 77 loci tested, 51 amplified consistently. When genotyped using 30 to 52 total adult trees from three old-growth populations in Maryland, Virginia, and New Jersey, USA, 23 loci were polymorphic. These loci exhibited four to 13 alleles, and observed and expected heterozygosities ranged from 0.233 to 0.865 and 0.272 to 0.876, respectively.

CONCLUSIONS

The microsatellite marker loci presented here will be valuable in population genetic studies of L. tulipifera because they do not suffer from ascertainment bias and show high polymorphism.

Patterns of Substitution Rate Variation at Many Nuclear Loci in Two Species Trios in the Brassicaceae Partitioned with ANOVA

There are marked variations among loci and among lineages in rates of nucleotide substitution. The generation time hypothesis (GTH) is a neutral explanation for substitution rate heterogeneity that has genomewide application, predicting that species with shorter generation times accumulate DNA sequence substitutions faster than species with longer generation times do since faster genome replication provides more opportunities for mutations to occur and reach fixation by genetic drift. Relatively few studies have rigorously evaluated the GTH in plants, and there are numerous alternative hypotheses for plant substitution rate variation. One major challenge has been finding pairs of closely related plant species with contrasting generation times and appropriate outgroup taxa that all also have DNA sequence data for numerous loci. To test for causes of rate variation, we obtained sequence data for 256 genes for Arabidopsis thaliana, normally reproducing every year, and the biennial Arabidopsis lyrata with three closely related outgroup taxa (Brassica rapa, Capsella grandiflora, and Neslia paniculata) as well as the biennial Brassica oleracea and the annual B. rapa lineage with the outgroup N. paniculata. A sign test indicated that more loci than expected by chance have faster rates of substitution on the branch leading to the annual than to the perennial for one three-species trio but not another. Tajima's 1D and 2D tests, and a likelihood ratio test that incorporated saturation correction, rejected rate homogeneity for up to 26 genes (up to 14 genes when correcting for multiple tests), consistently showing faster rates for the annual lineage in the Arabidopsis species trio. ANOVA showed significant rate heterogeneity between the Arabidopsis and Brassica species trios (about 6 % of rate variation) and among loci (about 26-32 % of rate variation). The lineage-by-locus interaction which would be caused by locus- and lineage-specific natural selection explained about 13 % of substitution rate variation in one ANOVA model using substitution rates from genes partitioned into odd and even codons but was not a significant effect without partitioned genes. Annual/perennial lineage and species trio by annual/perennial lineage each explained about 1 % of substitution rate variation.

Permanent genetic resources added to Molecular Ecology Resources Database 1 August 2010-30 September 2010

This article documents the addition of 229 microsatellite marker loci to the Molecular Ecology Resources Database. Loci were developed for the following species: Acacia auriculiformis × Acacia mangium hybrid, Alabama argillacea, Anoplopoma fimbria, Aplochiton zebra, Brevicoryne brassicae, Bruguiera gymnorhiza, Bucorvus leadbeateri, Delphacodes detecta, Tumidagena minuta, Dictyostelium giganteum, Echinogammarus berilloni, Epimedium sagittatum, Fraxinus excelsior, Labeo chrysophekadion, Oncorhynchus clarki lewisi, Paratrechina longicornis, Phaeocystis antarctica, Pinus roxburghii and Potamilus capax. These loci were cross-tested on the following species: Acacia peregrinalis, Acacia crassicarpa, Bruguiera cylindrica, Delphacodes detecta, Tumidagena minuta, Dictyostelium macrocephalum, Dictyostelium discoideum, Dictyostelium purpureum, Dictyostelium mucoroides, Dictyostelium rosarium, Polysphondylium pallidum, Epimedium brevicornum, Epimedium koreanum, Epimedium pubescens, Epimedium wushanese and Fraxinus angustifolia.

Properties of Weir and Cockerham's Fst estimators and associated bootstrap confidence intervals

Weir and Cockerham introduced single locus and multiloci F(st) estimators for the parameter θ. These estimators are commonly used, but little beyond their bias and variance is known. In this work, we develop formulas that allow us to describe how the underlying value of θ and the genetic diversity of sampled loci affect the distributions of these estimators. We show that in certain settings, these estimators are close to normal, while in others they are far from normal. We use these results to analyze confidence interval construction for θ, showing that the percentile-t bootstrap works well while the BCa bootstrap works poorly. Our results are derived using a novel coalescent based method.

Post-mating reproductive barriers in two unidirectionally hybridizing sunfish (Centrarchidae: Lepomis)

The evolutionary sequence of events in the evolution of reproductive barriers between species is at the core of speciation biology. Where premating barriers fail, post-mating barriers, such as conspecific sperm precedence (CSP), gamete incompatibility (GI) and hybrid inviability (HI) may evolve to prevent the production of (often) costly hybrid offspring with reduced fitness. We tested the role of post-mating mechanisms for the reproductive isolation between two sunfish species [bluegill (BG) Lepomis macrochirus and pumpkinseed (PS) Lepomis gibbosus] and their first-generation hybrids. Performing in vitro sperm competition experiments, we observed asymmetric CSP as main post-mating isolation mechanism when BG and PS sperm were competing for PS eggs, whereas when sperm from both species were competing for BG eggs it was HI. Furthermore, hybrid sperm--although fertile in the absence of competition--were outcompeted by sperm of either parental species. This result may at least partly explain previous observations that natural hybridization in the study system is unidirectional.

Reconsidering the generation time hypothesis based on nuclear ribosomal ITS sequence comparisons in annual and perennial angiosperms

Background

Differences in plant annual/perennial habit are hypothesized to cause a generation time effect on divergence rates. Previous studies that compared rates of divergence for internal transcribed spacer (ITS1 and ITS2) sequences of nuclear ribosomal DNA (nrDNA) in angiosperms have reached contradictory conclusions about whether differences in generation times (or other life history features) are associated with divergence rate heterogeneity. We compared annual/perennial ITS divergence rates using published sequence data, employing sampling criteria to control for possible artifacts that might obscure any actual rate variation caused by annual/perennial differences.

Results

Relative rate tests employing ITS sequences from 16 phylogenetically-independent annual/perennial species pairs rejected rate homogeneity in only a few comparisons, with annuals more frequently exhibiting faster substitution rates. Treating branch length differences categorically (annual faster or perennial faster regardless of magnitude) with a sign test often indicated an excess of annuals with faster substitution rates. Annuals showed an approximately 1.6-fold rate acceleration in nucleotide substitution models for ITS. Relative rates of three nuclear loci and two chloroplast regions for the annual Arabidopsis thaliana compared with two closely related Arabidopsis perennials indicated that divergence was faster for the annual. In contrast, A. thaliana ITS divergence rates were sometimes faster and sometimes slower than the perennial. In simulations, divergence rate differences of at least 3.5-fold were required to reject rate constancy in > 80 % of replicates using a nucleotide substitution model observed for the combination of ITS1 and ITS2. Simulations also showed that categorical treatment of branch length differences detected rate heterogeneity > 80% of the time with a 1.5-fold or greater rate difference.

Conclusion

Although rate homogeneity was not rejected in many comparisons, in cases of significant rate heterogeneity annuals frequently exhibited faster substitution rates. Our results suggest that annual taxa may exhibit a less than 2-fold rate acceleration at ITS. Since the rate difference is small and ITS lacks statistical power to reject rate homogeneity, further studies with greater power will be required to adequately test the hypothesis that annual and perennial plants have heterogeneous substitution rates. Arabidopsis sequence data suggest that relative rate tests based on multiple loci may be able to distinguish a weak acceleration in annual plants. The failure to detect rate heterogeneity with ITS in past studies may be largely a product of low statistical power.

Parallel rate heterogeneity in chloroplast and mitochondrial genomes of Brazil nut trees (Lecythidaceae) is consistent with lineage effects

We investigated whether relative rates of divergence were correlated between the mitochondrial and chloroplast genomes as expected under lineage effects or were genome specific as expected with locus-specific effects. Five mitochondrial noncoding regions (nad1B_C, nad4exon1_2, nad7exon2_3, nad7exon3_4, and rps14-cob) for 21 samples from Lecythidaceae were sequenced. Three chloroplast regions (rpl20-5'rps12, trnS-trnG, and psbA-trnH) were sequenced to expand the taxa in an existing data set. Absolute rates of nucleotide and insertion and deletion (indel) changes were 13 times faster in the chloroplast genome than in the mitochondrial genome. Similar indel length frequency distributions for both organelles suggested that common mechanisms were responsible for generating indels. Molecular clock tests applied to phylogenetic trees estimated from mitochondrial and chloroplast sequences revealed global rate heterogeneity of nucleotide substitution. Maximum likelihood and Tajima's 1D relative rate tests show that Lecythis zabucajo exhibited a rate acceleration for both the mitochondrial and chloroplast sequences. Whereas Eschweilera romeu-cardosoi showed a significant rate slowdown for chloroplast sequences, the mitochondrial sequences for 3 Eschweilera taxa showed evidence for a rate slowdown only when compared with L. zabucajo. Significant rate heterogeneity was also observed for indel changes in the mitochondrial genome but not for the chloroplast. The lack of mitochondrial nucleotide changes for some taxa as well as chloroplast indel homoplasy may have limited the power of relative rate tests to detect rate variation. Relative ratio tests consistently indicated rate proportionality among branch lengths between the mitochondrial and chloroplast phylogenetic trees. The relative ratio tests showed that taxa possessing rate heterogeneity had parallel relative divergence rates in both mitochondrial and chloroplast sequences as expected under lineage effects. A neutral replication-dependent model of rate heterogeneity for both nucleotide and indel changes provides a simple explanation for common patterns of rate heterogeneity across the 2 organelle genomes in Lecythidaceae. The lineage effects observed here were uncoupled from annual/perennial habit because all the species from this study are perennial.

Steady state of homozygosity and for the island model

We examine homozygosity and G(st) for a subdivided population governed by the finite island model. Assuming an infinite allele model and strong mutation we show that the steady state distributions of G(st) and homozygosity have asymptotic expansions in the mutation rate. We use this observation to derive asymptotic expansions for various moments of homozygosity and to derive rigorous formulas for the mean and variance of G(st). We show that G(st) approximately 1/(1+2Nm), similarly to the well known formula of Wright for the infinite island model, and that the variance of G(st) goes to zero as mutation increases.

Hierarchical components of genetic variation at a species boundary: population structure in two sympatric varieties of Lupinus microcarpus (Leguminosae)

Lupinus microcarpus is a self-compatible annual plant that forms a species complex of morphologically variable but indeterminate varieties. In order to examine the hypothesis that varieties of L. microcarpus comprise genetically differentiated and reproductively isolated species, populations of L. microcarpus var. horizontalis and var. densiflorus were sampled from an area of sympatry in central California and genotyped using six microsatellite loci. Bayesian clustering divided the total sample into two groups corresponding to the named varieties with extremely low levels of inferred coancestry. Similarly, maximum likelihood and distance methods for genetic assignment placed individuals in two nonoverlapping groups. Evidence for isolation by distance (IBD) within each variety was found at shorter distance classes, but varieties remained differentiated in sympatry. Furthermore, coalescent estimates of divergence time indicate separation within the past 950-5050 generations, with minimal gene flow after divergence. A four-level hierarchical analysis of molecular variance (amova) found significant levels of genetic differentiation among varieties (theta(P) = 0.292), populations within varieties (theta(S) = 0.449), subpopulations within populations (theta(SS) = 0.623), and individuals within subpopulations (f = 0.421); but the greatest degree of differentiation was at the subpopulation level. Although it is sometimes assumed that the magnitude of genetic differences (e.g. F(ST)) should be greater between species than among populations or subpopulations of the same species, shared ancestral polymorphism may lead to relatively low levels of differentiation at the species level, even as the stochastic effects of genetic drift generate higher levels of differentiation at lower hierarchical levels. These results suggest that L. microcarpus var. horizontalis and var. densiflorus are recently diverged yet reproductively isolated species, with high levels of inbreeding resulting from the combined effects of limited gene flow, demographic bottlenecks, and partial selfing in finite, geographically structured populations.

A finite locus effect diffusion model for the evolution of a quantitative trait

A diffusion model is constructed for the joint distribution of absolute locus effect sizes and allele frequencies for loci contributing to an additive quantitative trait under selection in a haploid, panmictic population. The model is designed to approximate a discrete model exactly in the limit as both population size and the number of loci affecting the trait tend to infinity. For the case when all loci have the same absolute effect size, formal multiple-timescale asymptotics are used to predict the long-time response of the population trait mean to selection. For the case where loci can take on either of two distinct effect sizes, not necessarily with equal probability, numerical solutions of the system indicate that response to selection of a quantitative trait is insensitive to the variability of the distribution of effect sizes when mutation is negligible.

Reconciling nuclear microsatellite and mitochondrial marker estimates of population structure: breeding population structure of Chesapeake Bay striped bass (Morone saxatilis)

Comparative analyses of nuclear and organelle genetic markers may help delineate evolutionarily significant units or management units, although population differentiation estimates from multiple genomes can also conflict. Striped bass (Morone saxatilis) are long-lived, highly migratory anadromous fish recently recovered from a severe decline in population size. Previous studies with protein, nuclear DNA and mitochondrial DNA (mtDNA) markers produced discordant results, and it remains uncertain if the multiple tributaries within Chesapeake Bay constitute distinct management units. Here, 196 young-of-the-year (YOY) striped bass were sampled from Maryland's Choptank, Potomac and Nanticoke Rivers and the north end of Chesapeake Bay in 1999 and from Virginia's Mataponi and Rappahannock Rivers in 2001. A total of 10 microsatellite loci exhibited between two and 27 alleles per locus with observed heterozygosities between 0.255 and 0.893. The 10-locus estimate of R(ST) among the six tributaries was -0.0065 (95% confidence interval -0.0198 to 0.0018). All R(ST) and all but one theta estimates for pairs of populations were not significantly different from zero. Reanalysis of Chesapeake Bay striped bass mtDNA data from two previous studies estimated population differentiation between theta=-0.002 and 0.160, values generally similar to mtDNA population differentiation predicted from microsatellite R(ST) after adjusting for reduced effective population size and uniparental inheritance in organelle genomes. Based on mtDNA differentiation, breeding sex ratios or gene flow may have been slightly male biased in some years. The results reconcile conflicting past studies based on different types of genetic markers, supporting a single Chesapeake Bay management unit encompassing a panmictic striped bass breeding population.

The impact of microsatellite electromorph size homoplasy on multilocus population structure estimates in a tropical tree (Corythophora alta) and an anadromous fish (Morone saxatilis)

Microsatellite allelic states are determined by electrophoretic sizing of polymerase chain reaction fragments to define electromorphs. Numerous studies have documented that identical microsatellite electromorphs are potentially heterogeneous at the DNA sequence level, a phenomenon called electromorph size homoplasy. Few studies have examined the impact of electromorph size homoplasy on estimates of population genetic parameters. We investigated the frequency of microsatellite electromorph size homoplasy for 12 loci in the tropical tree Corythophora alta and 11 loci in the anadromous fish Morone saxatilis by sequencing 14-23 homozygotes per locus sampled from multiple populations for a total of 453 sequences. Sequencing revealed no homoplasy for M. saxatilis loci. Seven C. alta loci exhibited homoplasy, including single and compound repeat motifs both with and without interruptions. Between 12.5 and 42.9% of electromorphs sampled per locus showed size homoplasy. Two methods of correction for homoplasy in C. alta generally produced little or no change in single-locus estimates of RST, except for two loci in which some additional differentiation among populations was revealed. Twelve-locus estimates of RST (including the seven loci corrected for homoplasy) were slightly greater than estimates from uncorrected data, although the 95% confidence intervals overlapped. The frequency of methodological errors such as clerical mistakes or sample mislabelling per genotype scored was estimated at 5.4 and 7.3% for C. alta and M. saxatilis, respectively. Simulations showed that the increase in RST produced by homoplasy correction was only slightly larger than variation in RST estimates expected to be caused by methodological errors.

Patterns and relative rates of nucleotide and insertion/deletion evolution at six chloroplast intergenic regions in new world species of the Lecythidaceae

Insertions and deletions (indels) in chloroplast noncoding regions are common genetic markers to estimate population structure and gene flow, although relatively little is known about indel evolution among recently diverged lineages such as within plant families. Because indel events tend to occur nonrandomly along DNA sequences, recurrent mutations may generate homoplasy for indel haplotypes. This is a potential problem for population studies, because indel haplotypes may be shared among populations after recurrent mutation as well as gene flow. Furthermore, indel haplotypes may differ in fitness and therefore be subject to natural selection detectable as rate heterogeneity among lineages. Such selection could contribute to the spatial patterning of cpDNA haplotypes, greatly complicating the interpretation of cpDNA population structure. This study examined both nucleotide and indel cpDNA variation and divergence at six noncoding regions (psbB-psbH, atpB-rbcL, trnL-trnH, rpl20-5'rps12, trnS-trnG, and trnH-psbA) in 16 individuals from eight species in the Lecythidaceae and a Sapotaceae outgroup. We described patterns of cpDNA changes, assessed the level of indel homoplasy, and tested for rate heterogeneity among lineages and regions. Although regression analysis of branch lengths suggested some degree of indel homoplasy among the most divergent lineages, there was little evidence for indel homoplasy within the Lecythidaceae. Likelihood ratio tests applied to the entire phylogenetic tree revealed a consistent pattern rejecting a molecular clock. Tajima's 1D and 2D tests revealed two taxa with consistent rate heterogeneity, one showing relatively more and one relatively fewer changes than other taxa. In general, nucleotide changes showed more evidence of rate heterogeneity than did indel changes. The rate of evolution was highly variable among the six cpDNA regions examined, with the trnS-trnG and trnH-psbA regions showing as much as 10% and 15% divergence within the Lecythidaceae. Deviations from rate homogeneity in the two taxa were constant across cpDNA regions, consistent with lineage-specific rates of evolution rather than cpDNA region-specific natural selection. There is no evidence that indels are more likely than nucleotide changes to experience homoplasy within the Lecythidaceae. These results support a neutral interpretation of cpDNA indel and nucleotide variation in population studies within species such as Corythophora alta.

Comparing relative rates of pollen and seed gene flow in the island model using nuclear and organelle measures of population structure

We describe a method for comparing nuclear and organelle population differentiation (F(ST)) in seed plants to test the hypothesis that pollen and seed gene flow rates are equal. Wright's infinite island model is used, with arbitrary levels of self-fertilization and biparental organelle inheritance. The comparison can also be applied to gene flow in animals. Since effective population sizes are smaller for organelle genomes than for nuclear genomes and organelles are often uniparentally inherited, organelle F(ST) is expected to be higher at equilibrium than nuclear F(ST) even if pollen and seed gene flow rates are equal. To reject the null hypothesis of equal seed and pollen gene flow rates, nuclear and organelle F(ST)'s must differ significantly from their expected values under this hypothesis. Finite island model simulations indicate that infinite island model expectations are not greatly biased by finite numbers of populations (>/=100 subpopulations). The power to distinguish dissimilar rates of pollen and seed gene flow depends on confidence intervals for fixation index estimates, which shrink as more subpopulations and loci are sampled. Using data from the tropical tree Corythophora alta, we rejected the null hypothesis that seed and pollen gene flow rates are equal but cannot reject the alternative hypothesis that pollen gene flow is 200 times greater than seed gene flow.

Altering fish embryos with aquaporin-3: an essential step toward successful cryopreservation

Fish populations are globally threatened by overharvesting and habitat degradation. The ability to bank fish embryos by cryopreservation could be crucial for preserving species diversity, for aquaculture (allowing circannual fish farming), and for managing fish models used in human biomedical research. However, no nonmammalian embryo has ever been successfully cryopreserved. For fish, low membrane permeability prevents cryoprotectants from entering the yolk to prevent cryodamage. Here, we present evidence of a membrane mechanism hindering cryopreservation of fish and propose a novel solution to this obstacle. Zebrafish (Danio rerio) embryos have rectifying membranes that allow water to leave but not to reenter readily. This feature may be an evolutionary trait that allows freshwater embryos to grow in hypoosmotic environments without osmoregulatory organs. However, this trait may also prevent successful fish embryo cryopreservation because both water and cryoprotectants must move into and out of cells. As a solution, we injected zebrafish embryos with mRNA for the aquaporin-3 water channel protein and demonstrated increased membrane permeability to water and to a cryoprotectant. Modeling indicates that sufficient cryoprotectant enters aquaporin-3-expressing zebrafish embryos to allow cryopreservation.

Cloned microsatellite repeats differ between 4-base restriction endonucleases

Simple sequence repeat (SSR) loci are an important marker type for population genetic studies despite the limitation that development of novel loci requires construction and screening of genomic DNA libraries. The common practice of size fractioning genomic DNA before cloning could lead to differential representation of SSR loci within genomic libraries. In addition, linkage mapping studies have shown that small numbers of SSR markers are not randomly distributed within the genomes from which they are isolated. From attempts to clone five SSR repeat sequences in two wild plant species we show that the numbers and repeat type of potential SSR markers depend on the restriction endonuclease used to sample the genome when constructing DNA libraries. This observation is consistent with unequal sampling of the genome by different restriction enzymes. However, as a group the five SSR repeat sequences are not associated with a given restriction enzyme, suggesting they are not clumped within the genome. Use of multiple restriction enzymes to construct DNA libraries may help ensure that cloned SSR loci are drawn from diverse locations in the genome, helping to meet the assumption of randomly located marker loci required for population genetic inferences.

Universal linker and ligation procedures for construction of genomic DNA libraries enriched for microsatellites

Microsatellite loci are highly informative genetic markers useful for population genetic studies, linkage mapping and parentage determination. Methods to identify novel microsatellite loci commonly use subtractive hybridization to enrich small-insert genomic libraries for repeat sequences. A critical step in enrichment is attachment of an oligonucleotide linker to genomic DNA fragments so that repeat-containing sequences can be recovered by PCR for cloning. Current linkers and ligation methods rely on single restriction enzymes to size-fraction genomic DNA and generate complementary ends. These restriction enzyme/linker combinations are often species-specific, give poor recovery of repeat-enriched DNA and yield library inserts that are not a broad sample of the genome. We have developed a blunt-end linker, named SNX for its restriction sites, that allows the use of combinations of restriction enzymes to digest the majority of genomic DNA into the 200-1000-bp range. SNX is attached to genomic DNA with a simultaneous ligation/restriction reaction that is highly efficient and improves recovery of sequences after subtractive hybridization. SNX can be used for microsatellite enrichment in any species, since ligation is independent of the restriction enzymes used to size-fraction genomic DNA. These methods improve current repeat-enrichment strategies, resulting in representative small-insert libraries with a very high proportion of positive clones.

Relatedness measured by oligonucleotide probe DNA fingerprints and an estimate of the mating system of Sea Lavender (Limonium carolinianum)

Using DNA fingerprint markers within species and populations of wild plants requires information on the relationship between fingerprint similarity and relatedness. We identified a hypervariable marker based on oliog(GATA)4-hybridization of DpnII-cut genomic DNA from Sea Lavender (Limonium carolinianum). Banding patterns were somatically stable and highly variable among unrelated individuals. Band molecular-weight sizing errors (as a percent of band molecular weight) were estimated at 0.44%±0.003 within gels and 0.76%±0.964 between gels. Band sizing errors defined a 99% confidence bin of ±0.95% (1.90% total) of molecular weight. Band-sharing estimates were based on this bin size and on variance estimates that compensate for non-independent comparisons. Band-sharing among nine unrelated individuals (θ) was 0.198±0.O11. Experimental pollinations designed to produce selfed, fulland half-sib progeny groups led to five selfed progeny groups and no outcrossed progeny (mean band-sharing, ovS=0.468±0.074). A linear regression between band-sharing (S) and relatedness (r) assuming 17% inbreeding was r=0.006+0.914*S (R(2)=0.973) and established the maximum amount of inbreeding. ovS(0.392±0.022) estimated from wild pollinated seeds from four maternal families was intermediate to unrelated individuals and experimental selfed progeny, giving evidence for mixed mating in wild plants. More extensive plant pedigrees with known levels of inbreeding will be needed to measure variation in the relationship between S and r among populations and families.