Combining multiple investigative approaches to unravel functional responses to global change in the understorey of temperate forests.

Plant communities are being exposed to changing environmental conditions all around the globe, leading to alterations in plant diversity, community composition, and ecosystem functioning. For herbaceous understorey communities in temperate forests, responses to global change are postulated to be complex, due to the presence of a tree layer that modulates understorey responses to external pressures such as climate change and changes in atmospheric nitrogen deposition rates. Multiple investigative approaches have been put forward as tools to detect, quantify and predict understorey responses to these global-change drivers, including, among others, distributed resurvey studies and manipulative experiments. These investigative approaches are generally designed and reported upon in isolation, while integration across investigative approaches is rarely considered. In this study, we integrate three investigative approaches (two complementary resurvey approaches and one experimental approach) to investigate how climate warming and changes in nitrogen deposition affect the functional composition of the understorey and how functional responses in the understorey are modulated by canopy disturbance, that is, changes in overstorey canopy openness over time. Our resurvey data reveal that most changes in understorey functional characteristics represent responses to changes in canopy openness with shifts in macroclimate temperature and aerial nitrogen deposition playing secondary roles. Contrary to expectations, we found little evidence that these drivers interact. In addition, experimental findings deviated from the observational findings, suggesting that the forces driving understorey change at the regional scale differ from those driving change at the forest floor (i.e., the experimental treatments). Our study demonstrates that different approaches need to be integrated to acquire a full picture of how understorey communities respond to global change.

Plastic and quantitative genetic divergence mirror environmental gradients among wild, fragmented populations of Impatiens capensis.

PREMISE: Habitat fragmentation generates molecular genetic divergence among isolated populations, but few studies have assessed phenotypic divergence and fitness in populations where the genetic consequences of habitat fragmentation are known. Phenotypic divergence could reflect plasticity, local adaptation, and/or genetic drift. METHODS: We examined patterns and potential drivers of phenotypic divergence among 12 populations of jewelweed (Impatiens capensis) that show strong molecular genetic signals of isolation and drift among fragmented habitats. We measured morphological and reproductive traits in both maternal plants within natural populations and their self-fertilized progeny grown together in a common garden. We also quantified environmental divergence between home sites and the common garden. RESULTS: Populations with less molecular genetic variation expressed less maternal phenotypic variation. Progeny in the common garden converged in phenotypes relative to their wild mothers but retained among-population differences in morphology, survival, and reproduction. Among-population phenotypic variance was 3-10× greater in home sites than in the common garden for 6 of 7 morphological traits measured. Patterns of phenotypic divergence paralleled environmental gradients in ways suggestive of adaptation. Progeny resembled their mothers less as the environmental distance between their home site and the common garden increased. CONCLUSIONS: Despite strong molecular signatures of isolation and drift, phenotypic differences among these Impatiens populations appear to reflect both adaptive quantitative genetic divergence and plasticity. Quantifying the extent of local adaptation and plasticity and how these covary with molecular and phenotypic variation help us predict when populations may lose their adaptive capacity.

Functional traits mediate individualistic species-environment distributions at broad spatial scales while fine-scale species associations remain unpredictable.

PREMISE: Numerous processes influence plant distributions and co-occurrence patterns, including ecological sorting, limiting similarity, and stochastic effects. To discriminate among these processes and determine the spatial scales at which they operate, we investigated how functional traits and phylogenetic relatedness influence the distribution of temperate forest herbs. METHODS: We surveyed understory plant communities across 257 forest stands in Wisconsin and Michigan (USA) and applied Bayesian phylogenetic linear mixed-effects models (PGLMMs) to quantify how functional traits and phylogenetic relatedness influence the environmental distribution of 139 herbaceous plant species along broad edaphic, climatic, and light gradients. These models also allowed us to test how functional and phylogenetic similarity affect species co-occurrence within microsites. RESULTS: Leaf height, specific leaf area, and seed mass all influenced individualistic plant distributions along landscape-scale gradients in soil texture, soil fertility, light availability, and climate. In contrast, phylogenetic relationships did not consistently predict species-environment relationships. Neither functionally similar nor phylogenetically related herbs segregated among microsites within forest stands. CONCLUSIONS: Trait-mediated ecological sorting appears to drive temperate-forest community assembly, generating individualistic plant distributions along regional environmental gradients. This finding links classic studies in plant ecology and prior research in plant physiological ecology to current trait-based approaches in community ecology. However, our results fail to support the common assumption that limiting similarity governs local plant co-occurrences. Strong ecological sorting among forest stands coupled with stochastic fine-scale interactions among species appear to weaken deterministic, niche-based assembly processes at local scales.

Stresses affect inbreeding depression in complex ways: disentangling stress-specific genetic effects from effects of initial size in plants.

The magnitude of inbreeding depression (ID) varies unpredictably among environments. ID often increases in stressful environments suggesting that these expose more deleterious alleles to selection or increase their effects. More simply, ID could increase under conditions that amplify phenotypic variation (CV²), e.g., by accentuating size hierarchies among plants. These mechanisms are difficult to distinguish when stress increases both ID and phenotypic variation. We grew in- and outbred progeny of Mimulus guttatus under six abiotic stress treatments (control, waterlogging, drought, nutrient deficiency, copper addition, and clipping) with and without competition by the grass Poa palustris. ID differed greatly among stress treatments with δ varying from 7% (control) to 61% (waterlogging) but did not consistently increase with stress intensity. Poa competition increased ID under nutrient deficiency but not other stresses. Analyzing effects of initial size on performance of outbred plants suggests that under some conditions (low N, clipping) competition increased ID by amplifying initial size differences. In other cases (e.g., high ID under waterlogging), particular environments amplified the deleterious genetic effects of inbreeding suggesting differential gene expression. Interestingly, conditions that increased the phenotypic variability of inbred progeny regularly increased ID whereas variability among outbred progeny showed no relationship to ID. Our study reconciles the stress- and phenotypic variability hypotheses by demonstrating how specific conditions (rather than stress per se) act to increase ID. Analyzing CV² separately in inbred and outbred progeny while including effects of initial plant size improve our ability to predict how ID and gene expression vary across environments.

Shifts in precipitation and agricultural intensity increase phosphorus concentrations and loads in an agricultural watershed.

Fertilizers and manure applied to cropland to increase yields are often lost via surface erosion, soil leaching, and runoff, increasing nutrient loads in surface and sub-surface waters, degrading water quality, and worsening the 'dead zone' in the Gulf of Mexico. We leverage spatial and temporal variation in agricultural practices and precipitation events to examine how these factors affect stream total phosphorus (TP) concentrations and loads in the Sugar River (Wisconsin), recently listed as impaired. To perform our analysis, we first collected water quality data from 1995 to 2017 from 40 sites along the Sugar River and its tributaries. Starting in 2004, three dairy farms expanded to become concentrated animal feeding operations (CAFOs) in this watershed. We then estimated how time of year, stream position, discharge volume, and proximity to the newly expanded CAFOs affected TP concentrations and loads. Total P concentrations, which ranged from 0.02 to 1.4 mg/L and often exceeded the EPA surface water standard of 0.1 mg/L, increased with increases in stream discharge and proximity to dairy operations, peaking in early spring to mid-summer coincident with extreme precipitation events. Our empirical analysis also shows that TP concentrations downstream from the newly permitted CAFOs increased by 19% relative to upstream concentrations. When examining total daily phosphorus loads (concentration × discharge) from this 780 km2 watershed, we found that loads ranged from 5.88 to 4801 kg. Compared to upstream TP loads, those downstream from the CAFOs increased by 91% after the expansions - over four times that of concentration increases - implying that the rate of downstream phosphorus transfer has increased due to CAFO expansion. Our results argue for standards that focus on loads rather than concentrations and monitoring that includes peak events. As agriculture intensifies and extreme rainfall events become more frequent, it becomes increasingly important to limit soil and TP runoff from manure and fertilizer. Siting CAFOs carefully, limiting their size, and improving farming practices in proximity to CAFOs in spring and early summer could considerably reduce nutrient loads.

Adaptive associations among life history, reproductive traits, environment, and origin in the Wisconsin angiosperm flora.

PREMISE: We tested 25 classic and novel hypotheses regarding trait-origin, trait-trait, and trait-environment relationships to account for flora-wide variation in life history, habit, and especially reproductive traits using a plastid DNA phylogeny of most native (96.6%, or 1494/1547 species) and introduced (87.5%, or 690/789 species) angiosperms in Wisconsin, USA. METHODS: We assembled data on life history, habit, flowering, dispersal, mating system, and occurrence across open/closed/mixed habitats across species in the state phylogeny. We used phylogenetically structured analyses to assess the strength and statistical significance of associations predicted by our models. RESULTS: Introduced species are more likely to be annual herbs, occupy open habitats, have large, visually conspicuous, hermaphroditic flowers, and bear passively dispersed seeds. Among native species, hermaphroditism is associated with larger, more conspicuous flowers; monoecy is associated with small, inconspicuous flowers and passive seed dispersal; and dioecy is associated with small, inconspicuous flowers and fleshy fruits. Larger flowers with more conspicuous colors are more common in open habitats, and in understory species flowering under open (spring) canopies; fleshy fruits are more common in closed habitats. Wind pollination may help favor dioecy in open habitats. CONCLUSIONS: These findings support predictions regarding how breeding systems depend on flower size, flower color, and fruit type, and how those traits depend on habitat. This study is the first to combine flora-wide phylogenies with complete trait databases and phylogenetically structured analyses to provide powerful tests of evolutionary hypotheses about reproductive traits and their variation with geographic source, each other, and environmental conditions.

Drift happens: Molecular genetic diversity and differentiation among populations of jewelweed (Impatiens capensis Meerb.) reflect fragmentation of floodplain forests.

Landscape features often shape patterns of gene flow and genetic differentiation in plant species. Populations that are small and isolated enough also become subject to genetic drift. We examined patterns of gene flow and differentiation among 12 floodplain populations of the selfing annual jewelweed (Impatiens capensis Meerb.) nested within four river systems and two major watersheds in Wisconsin, USA. Floodplain forests and marshes provide a model system for assessing the effects of habitat fragmentation within agricultural/urban landscapes and for testing whether rivers act to genetically connect dispersed populations. We generated a panel of 12,856 single nucleotide polymorphisms and assessed genetic diversity, differentiation, gene flow, and drift. Clustering methods revealed strong population genetic structure with limited admixture and highly differentiated populations (mean multilocus FST  = 0.32, FST ' = 0.33). No signals of isolation by geographic distance or environment emerged, but alleles may flow along rivers given that genetic differentiation increased with river distance. Differentiation also increased in populations with fewer private alleles (R2  = 0.51) and higher local inbreeding (R2  = 0.22). Populations varied greatly in levels of local inbreeding (FIS  = 0.2-0.9) and FIS increased in more isolated populations. These results suggest that genetic drift dominates other forces in structuring these Impatiens populations. In rapidly changing environments, species must migrate or genetically adapt. Habitat fragmentation limits both processes, potentially compromising the ability of species to persist in fragmented landscapes.

Conditional vulnerability of plant diversity to atmospheric nitrogen deposition across the United States.

Atmospheric nitrogen (N) deposition has been shown to decrease plant species richness along regional deposition gradients in Europe and in experimental manipulations. However, the general response of species richness to N deposition across different vegetation types, soil conditions, and climates remains largely unknown even though responses may be contingent on these environmental factors. We assessed the effect of N deposition on herbaceous richness for 15,136 forest, woodland, shrubland, and grassland sites across the continental United States, to address how edaphic and climatic conditions altered vulnerability to this stressor. In our dataset, with N deposition ranging from 1 to 19 kg N⋅ha(-1)⋅y(-1), we found a unimodal relationship; richness increased at low deposition levels and decreased above 8.7 and 13.4 kg N⋅ha(-1)⋅y(-1) in open and closed-canopy vegetation, respectively. N deposition exceeded critical loads for loss of plant species richness in 24% of 15,136 sites examined nationwide. There were negative relationships between species richness and N deposition in 36% of 44 community gradients. Vulnerability to N deposition was consistently higher in more acidic soils whereas the moderating roles of temperature and precipitation varied across scales. We demonstrate here that negative relationships between N deposition and species richness are common, albeit not universal, and that fine-scale processes can moderate vegetation responses to N deposition. Our results highlight the importance of contingent factors when estimating ecosystem vulnerability to N deposition and suggest that N deposition is affecting species richness in forested and nonforested systems across much of the continental United States.

Spatial phylogenetics reveals evolutionary constraints on the assembly of a large regional flora.

PREMISE OF THE STUDY: We used spatial phylogenetics to analyze the assembly of the Wisconsin flora, linking processes of dispersal and niche evolution to spatial patterns of floristic and phylogenetic diversity and testing whether phylogenetic niche conservatism can account for these patterns. METHODS: We used digitized records and a new molecular phylogeny for 93% of vascular plants in Wisconsin to estimate spatial variation in species richness and phylogenetic α and ß diversity in a native flora shaped mainly by postglacial dispersal and response to environmental gradients. We developed distribution models for all species and used these to infer fine-scale variation in potential diversity, phylogenetic distance, and interspecific range overlaps. We identified 11 bioregions based on floristic composition, mapped areas of neo- and paleo-endemism to establish new conservation priorities and predict how community-assembly patterns should shift with climatic change. KEY RESULTS: Spatial phylogenetic turnover most strongly reflects differences in temperature and spatial distance. For all vascular plants, assemblages shift from phylogenetically clustered to overdispersed northward, contrary to most other studies. This pattern is lost for angiosperms alone, illustrating the importance of phylogenetic scale. CONCLUSIONS: Species ranges and assemblage composition appear driven primarily by phylogenetic niche conservatism. Closely related species are ecologically similar and occupy similar territories. The average level and geographic structure of plant phylogenetic diversity within Wisconsin are expected to greatly decline over the next half century, while potential species richness will increase throughout the state. Our methods can be applied to allochthonous communities throughout the world.

Can functional traits account for phylogenetic signal in community composition?

Phylogenetic and functional trait-based analyses inform our understanding of community composition, yet methods for quantifying the overlap in information derived from functional traits and phylogenies remain underdeveloped. Does adding traits to analyses of community composition reduce the phylogenetic signal in the residual variation? If not, then measured functional traits alone may be insufficient to explain community assembly. We propose a general statistical framework to quantify the proportion of phylogenetic pattern in community composition that remains after including measured functional traits. We then illustrate the framework with applications to two empirical data sets. Both data sets showed strong phylogenetic attraction, with related species likely to co-occur in the same communities. In one data set, including traits eliminated all phylogenetic signals in the residual variation of both abundance and presence/absence patterns. In the second data set, including traits reduced phylogenetic signal in residuals by 25% and 98% for abundance and presence/absence data, respectively. Our framework provides an explicit way to estimate how much phylogenetic community pattern remains in the residual variation after including measured functional traits. Knowing that functional traits account for most of the phylogenetic pattern should provide confidence that important traits for phylogenetic community structure have been identified. Conversely, knowing that there is unexplained residual phylogenetic information should spur the search for additional functional traits or other processes underlying community assembly.

Tracking lags in historical plant species' shifts in relation to regional climate change.

Can species shift their distributions fast enough to track changes in climate? We used abundance data from the 1950s and the 2000s in Wisconsin to measure shifts in the distribution and abundance of 78 forest-understory plant species over the last half-century and compare these shifts to changes in climate. We estimated temporal shifts in the geographic distribution of each species using vectors to connect abundance-weighted centroids from the 1950s and 2000s. These shifts in distribution reflect colonization, extirpation, and changes in abundance within sites, separately quantified here. We then applied climate analog analyses to compute vectors representing the climate change that each species experienced. Species shifted mostly to the northwest (mean: 49 ± 29 km) primarily reflecting processes of colonization and changes in local abundance. Analog climates for these species shifted even further to the northwest, however, exceeding species' shifts by an average of 90 ± 40 km. Most species thus failed to match recent rates of climate change. These lags decline in species that have colonized more sites and those with broader site occupancy, larger seed mass, and higher habitat fidelity. Thus, species' traits appear to affect their responses to climate change, but relationships are weak. As climate change accelerates, these lags will likely increase, potentially threatening the persistence of species lacking the capacity to disperse to new sites or locally adapt. However, species with greater lags have not yet declined more in abundance. The extent of these threats will likely depend on how other drivers of ecological change and interactions among species affect their responses to climate change.

Large, connected floodplain forests prone to flooding best sustain plant diversity.

Dams, levees, and water withdrawals disrupt hydrologic regimes and associated floodplain forests. Because these forests are also responding to changes in land use, species invasions, and climate change, the relative effects of these factors are hard to disentangle. Most studies of floodplain forests lack historic data, requiring us to rely on recent data or contemporary spatial relationships to these drivers to infer those causes of vegetation dynamics. Here, we use survey data from the 1950s to reconstruct plant community changes across 40 floodplain forests in Wisconsin. We applied two partial least squares regression (PLS) models to evaluate how current site and landscape scale conditions and changes in these conditions since the 1950s influence contemporary patterns of community diversity and composition. Local site variables were among the most important in explaining current composition metrics and their changes, but historic landscape variables and changes in these were also important. Current local diversity (α) was the highest at sites prone to frequent flooding, even at sites in fragmented landscapes. Sites along sinuous rivers in large watershed areas with more contiguous forest had the highest abundance of wetland indicator plants in the re-survey and had the largest increases in α diversity since the 1950s, despite having the highest presence of exotic species then. These same sites have converged in composition, reflecting increases in wetland indicator plants and common native species. These patterns of increasing α diversity coupled with declines in community distinctiveness are uncommon among long-term studies. Increases in wetland plants may indicate that sites have become wetter with hydrologic changes, but these increases may also reflect improved colonization and establishment processes involving a robust regional pool of generalist wetland taxa. Woody and exotic plants typical of upland forests increased at rarely flooded sites in fragmented and urbanizing landscapes, indicating shifts towards a later-successional conditions and a dampened disturbance regime. This has reduced local species diversity and increased regional distinctness at some sites. As hydrologic connections appear to best maintain native species diversity and composition, even in fragmented landscapes, managers should seek to recreate these whenever feasible.

Combining community resurvey data to advance global change research.

More and more ecologists have started to resurvey communities sampled in earlier decades to determine long-term shifts in community composition and infer the likely drivers of the ecological changes observed. However, to assess the relative importance of, and interactions among, multiple drivers joint analyses of resurvey data from many regions spanning large environmental gradients are needed. In this paper we illustrate how combining resurvey data from multiple regions can increase the likelihood of driver-orthogonality within the design and show that repeatedly surveying across multiple regions provides higher representativeness and comprehensiveness, allowing us to answer more completely a broader range of questions. We provide general guidelines to aid implementation of multi-region resurvey databases. In so doing, we aim to encourage resurvey database development across other community types and biomes to advance global environmental change research.

Microclimate moderates plant responses to macroclimate warming.

Recent global warming is acting across marine, freshwater, and terrestrial ecosystems to favor species adapted to warmer conditions and/or reduce the abundance of cold-adapted organisms (i.e., "thermophilization" of communities). Lack of community responses to increased temperature, however, has also been reported for several taxa and regions, suggesting that "climatic lags" may be frequent. Here we show that microclimatic effects brought about by forest canopy closure can buffer biotic responses to macroclimate warming, thus explaining an apparent climatic lag. Using data from 1,409 vegetation plots in European and North American temperate forests, each surveyed at least twice over an interval of 12-67 y, we document significant thermophilization of ground-layer plant communities. These changes reflect concurrent declines in species adapted to cooler conditions and increases in species adapted to warmer conditions. However, thermophilization, particularly the increase of warm-adapted species, is attenuated in forests whose canopies have become denser, probably reflecting cooler growing-season ground temperatures via increased shading. As standing stocks of trees have increased in many temperate forests in recent decades, local microclimatic effects may commonly be moderating the impacts of macroclimate warming on forest understories. Conversely, increases in harvesting woody biomass--e.g., for bioenergy--may open forest canopies and accelerate thermophilization of temperate forest biodiversity.

Genetic rescue: a safe or risky bet?

Small and isolated populations face threats from genetic drift and inbreeding. To rescue populations from these threats, conservation biologists can augment gene flow into small populations to increase variation and reduce inbreeding depression. Spectacular success stories include greater prairie chickens in Illinois (Westermeier et al. ), adders in Sweden (Madsen et al. ) and panthers in Florida (Johnson et al. ). However, we also know that performing such crosses risks introducing genes that may be poorly adapted to local conditions or genetic backgrounds. A classic example of such 'outbreeding depression' occurred when different subspecies of ibex from Turkey and the Sinai were introduced to assist recovery of an ibex population in Czechoslovakia (Templeton ). Despite being fertile, the hybrids birthed calves too early, causing the whole population to disappear. In the face of uncertainty, conservation biologists have tended to respect genetic identity, shying away from routinely crossing populations. In this issue of Molecular Ecology, Frankham () compiles empirical data from experimental studies to assess the costs and benefits of between-population crosses (Fig. ). Crosses screened to exclude those involving highly divergent populations or distinct habitats show large heterosis with few apparent risks of outbreeding depression. This leads Frankham to advocate for using assisted gene flow more widely. But do the studies analysed in this meta-analysis adequately test for latent outcrossing depression?

Trait--environment relationships remain strong despite 50 years of trait compositional change in temperate forests.

Temperate North American forest communities have changed considerably in response to logging, fragmentation, herbivory, and other global change factors. Significant changes in the structure and composition of seemingly undisturbed Wisconsin forest communities have occurred over the past 50 years, including widespread declines in alpha and beta species diversity. To investigate how shifts in species composition have affected distributions of plant functional traits, we first compiled extensive data on understory plant species traits. We then computed community-weighted trait means and functional diversity metrics for communities in both the 1950s and 2000s. We examined how trait values and diversity varied across environmental gradients and among Wisconsin's four main ecoregions. Trait means and diversity values reflect conspicuous gradients in species composition, soils, and climatic conditions. Over the past 50 years, values of most traits have changed as communities shifted toward species with higher leaf nutrient levels and specific leaf area, particularly in the southern ecoregions. Trait richness and diversity have declined, particularly in historically species- and trait-rich unglaciated southwestern Wisconsin. Reductions in within-site trait diversity may be diminishing the ability of these forest communities to resist or resiliently respond to shifts in environmental conditions. Despite changes in trait and community composition, trait-environment relationships measured directly via fourth-corner analysis remain strong for most plant traits. Nevertheless, accelerating ecological change (including climate change) could outstrip the ability of plant species and traits to match their environment, particularly in more fragmented landscapes.

Long-term nitrogen deposition reduces the diversity of nitrogen-fixing plants.

Biological nitrogen fixation is a fundamental part of ecosystem functioning. Anthropogenic nitrogen deposition and climate change may, however, limit the competitive advantage of nitrogen-fixing plants, leading to reduced relative diversity of nitrogen-fixing plants. Yet, assessments of changes of nitrogen-fixing plant long-term community diversity are rare. Here, we examine temporal trends in the diversity of nitrogen-fixing plants and their relationships with anthropogenic nitrogen deposition while accounting for changes in temperature and aridity. We used forest-floor vegetation resurveys of temperate forests in Europe and the United States spanning multiple decades. Nitrogen-fixer richness declined as nitrogen deposition increased over time but did not respond to changes in climate. Phylogenetic diversity also declined, as distinct lineages of N-fixers were lost between surveys, but the "winners" and "losers" among nitrogen-fixing lineages varied among study sites, suggesting that losses are context dependent. Anthropogenic nitrogen deposition reduces nitrogen-fixing plant diversity in ways that may strongly affect natural nitrogen fixation.

The impacts of an invasive species citizen science training program on participant attitudes, behavior, and science literacy.

Citizen science can make major contributions to informal science education by targeting participants' attitudes and knowledge about science while changing human behavior towards the environment. We examined how training associated with an invasive species citizen science program affected participants in these areas. We found no changes in science literacy or overall attitudes between tests administered just before and after a one-day training program, matching results from other studies. However, we found improvements in science literacy and knowledge using context-specific measures and in self-reported intention to engage in pro-environmental activities. While we noted modest change in knowledge and attitudes, we found comparison and interpretation of these data difficult in the absence of other studies using similar measures. We suggest that alternative survey instruments are needed and should be calibrated appropriately to the pre-existing attitudes, behavior, and levels of knowledge in these relatively sophisticated target groups.

Failure to purge: population and individual inbreeding effects on fitness across generations of wild Impatiens capensis.

Inbreeding exposes deleterious recessive alleles in homozygotes, lowering fitness and generating inbreeding depression (ID). Both purging (via selection) and fixation (via drift) should reduce segregating deleterious mutations and ID in more inbred populations. These theoretical predictions are not well-tested in wild populations, which is concerning given purging/fixation have opposite fitness outcomes. We examined how individual- and population-level inbreeding and genomic heterozygosity affected maternal and progeny fitness within and among 12 wild populations of Impatiens capensis. We quantified maternal fitness in home sites, maternal multilocus heterozygosity (using 12,560 single-nucleotide polymorphisms), and lifetime fitness of selfed and predominantly outcrossed progeny in a common garden. These populations spanned a broad range of individual-level (fi = -0.17-0.98) and population-level inbreeding (FIS = 0.25-0.87). More inbred populations contained fewer polymorphic loci, less fecund mothers, and smaller progeny, suggesting higher fixed loads. However, despite appreciable ID (mean: 8.8 lethal equivalents per gamete), ID did not systematically decline in more inbred population. More heterozygous mothers were more fecund and produced fitter progeny in outcrossed populations, but this pattern unexpectedly reversed in highly inbred populations. These observations suggest that persistent overdominance or some other force acts to forestall purging and fixation in these populations.