Abstract

The influence of beaver (Castor canadensis Kuhl) herbivory on Salix reproduction, specifically the stimulation of asexual reproduction via browsed stem fragments, is relatively unknown. This study aimed to determine if beaver herbivory stimulates asexual reproduction of riparian willows and results in mature populations dominated by clones. The survival of seedlings and asexual propagules produced by beaver browse in populations of the riparian willow Salix sitchensis (Sanson in Bongard) were quantified to determine overwinter survival at 6 experimental sites. Salix sitchensis clonal diversity, using five microsatellite markers and the polymerase chain reaction, was calculated to detect if asexual reproduction had been stimulated by beaver herbivory. No sexual propagules survived overwinter in any of our study sites. Numerous asexual propagules were observed and 0–41% survived overwinter. Each sampled individual possessed a unique multilocus genotype, and clonal diversity was 1.0. Beaver herbivory did not create current willow populations dominated by clones. Beaver herbivory and asexual reproduction appeared to have played a minor role in the reproductive strategies of S. sitchensis at our sites in central British Columbia, Canada.

1. Introduction

Whether herbivory positively or negatively influences plant reproduction is a debated topic. Depending on the species, response variable, and the seasonal pattern of the foraging, herbivory can have a negative, [17], positive [1, 5, 8], or neutral [1, 4, 5] influence on plant reproduction. A relatively unexplored aspect of this interaction is the stimulation of asexual reproduction by the production and relocation of vegetative propagules via animal browsing. While some studies have reported the incidence or survival of such translocated propagules [915], it is still unclear if they significantly contribute to plant reproductive strategies.

Riparian habitat in central British Columbia, Canada, presents an ideal system to investigate this relationship. The dominant members of the riparian vegetative community are species of the genus Salix. Willows (Salix spp.) reproduce either sexually through wind dispersed seeds, or asexually through underground runners and severed branches or stems [1012, 1518]. Severed branches may play an important role in Salix colonization [1923], once colonization has occurred, these areas of patchy vegetation can shelter sexual propagules, as well as stabilize and retain sediment. This ecosystem modification creates more habitat conducive to sexual propagule survival [1923]. The opposite sequence of events has been observed for other populations of riparian vegetation. In populations of Populus nigra Linnaeus (Populus is a closely related genus to Salix within the family Salicaceae), vegetative reproduction dominated within sites and sexual reproduction initiated colonization [24]. Other species have shown a mixed reproductive strategy within a site, while sexual reproduction was responsible for colonization [25]. High genetic diversity observed in most species of riparian willows [10, 15, 2629], as well as the abundance of sexual propagules produced by riparian willows [11, 13, 14] suggests that sexual reproduction may play an important role in willow reproduction.

In the riparian habitats of central British Columbia Canada, beavers (Castor canadensis Kuhl) leave evidence of their presence in the form of lodges, dams, ponds, and felled trees. Through selective foraging and the creation of dams, beavers alter forest composition and age, as well as the hydrology, temperature, and chemical properties of water ways [30, 31]. Herbivory by beavers may alter the nutrient and secondary compound concentrations of willows [32], as well as increase basal area, diameter, abundance, and sapling recruitment [33, 34]. Beavers may also influence the reproductive strategies of willows by removing and relocating stems during foraging, creating potential asexual propagules. This mechanism has been hypothesized to contribute to Salix planifolia’s (Pursh) domination of mire habitat [9]. Conversely, other studies have suggested that beaver foraging has limited impact on willow reproduction [1114].

Numerous potential asexual propagules, in the form of branches severed by beavers, are found in the riparian habitats of central British Columbia. It is possible that beaver herbivory may induce willow asexual reproduction and play an important role in willow reproductive strategies. The aim of this study was to examine the impact of beaver browse on willow reproduction with a particular focus on Salix sitchensis (Sanson in Bongard), one of the predominant and most heavily browsed species in riparian habitats in central British Columbia [35, 36]. Specifically, we tested the prediction that the large numbers of potential asexual propagules produced by beaver foraging would lead to mature stands of willows dominated by clonally derived individuals. To do this, we used both field and laboratory analyses to investigate the recent and long-term reproductive history of the site, respectively. In the field we examined the establishment of seedlings and vegetative fragments of S. sitchensis at six riparian sites with beavers present. The pattern of sexual and asexual establishment in the field is a product of the site’s recent reproductive history. Molecular genetics can be used to determine how a population has been reproducing on a longer time scale. Asexual reproduction should lead to a high predominance of clonally derived genotypes (see [15] for an example). We therefore examined clonal diversity using five microsatellite markers in populations of S. sitchensis at six riparian sites.

2. Methods

2.1. Study Sites

Six riparian sites, spread over 8000 km2 of central British Columbia Canada, were chosen for this study (Figure 1; Table 1). No effort was made to discriminate sites based upon beaver abundance. Three of the six sites were established on large rivers: Bowron River and Bowron Bridge on the Bowron River, as well as Penny by the Fraser on the Fraser River. These sites were dominated by sandy and rocky beaches extending 25–50 m from the water to a dense forest composed of Salix, Alnus, Populus, Betula, Pinus, Picea, and Abies species. Present sporadically along the beach were small Salix, Alnus, and Populus individuals. Beavers at these sites did not build dams, but did construct lodges and food caches. Island Park, a floodplain island in the Nechako River, was dominated by a dense grove of tall Salix. Areas not dominated by willows were often covered with grasses and sedges, although substantial expanses of exposed sediment were also observed. Beavers at this site built lodges, several dams along an oxbow, and numerous food caches. Despite the presence of small beaver dams, no substantial beaver pond or beaver-induced flooding was observed at this site. Two sites, Camp Creek and Penny, were located on smaller creeks frequently dammed by beavers. Sporadic clumps of Salix and Alnus were found near the water while a dense mixed forest dominated 50–75 m from the water. The majority of these sites were dominated by shrubs, sedges and grasses, however, small sections of exposed sediment near the creeks were also present. Camp Creek and Penny possessed beaver lodges, dams, food caches, and sizeable beaver ponds. More detailed descriptions of all sites can be found in Gerwing [35, 36].

2.2. Field Survey of the Short-Term Reproductive History

During August 2007 and April–August in 2008, S. sitchensis seedlings were quantified on the portion of the study sites classified as potential germination substrate. Potential substrate was defined as at least partly exposed soil anywhere in the site, not completely covered by grass, sedges, shrubs, or water. The dimensions of this habitat were measured in situ for area calculations. Covered habitats were avoided as Cottrell [9] found no willow reproduction under sedge cover, and Barsoum [11] found the majority of seedlings on fine exposed sediments. The majority of the potential germination substrate was found near the water’s edge.

Beaver-browsed stem fragments were identified by the characteristic diagonal browse pattern left on vegetative propagules. Beaver browse was counted only on potential germination substrate. In the falls of 2007 and 2008, stems browsed and abandoned by beavers, but possessing adventitious roots were denoted as rooted propagules. Those without roots were defined as unrooted propagules. To establish the overwinter survival of beaver-produced asexual propagules within the sites Camp Creek, Bowron Bridge, and Island Park, each category was marked with a characteristic colour of flagging tape in 2007 and located again the following spring. Only S. sitchensis stems were marked and identification was made possible by S. sitchensis’ characteristic leaf and bark morphology [37].

2.3. Laboratory Analysis of Longer Term Reproductive History

Leaf samples for genetic analysis were collected from mature S. sitchensis at all sites except Penny by the Fraser (Table 1). At Camp Creek and Penny, all individuals were sampled and a random subset was selected for analysis. At Bowron River and Bowron Bridge, where the mature populations formed narrow bands of individuals along the rivers, one individual was sampled every 10 m for 160 m and 200 m, respectively. At Island Park ten, , plots separated by 25 m, were established along the center of the island. Along the south side of the island a 200 m transect was also established. All S. sitchensis were sampled within these plots and a random subset of these individuals was analyzed. Samples were collected randomly as our goal was to determine if populations were dominated by clones not to examine fine scale clonal structure.

Leaf tissue was collected and stored at −80°C. Genomic DNA was isolated using QIAGEN DNeasy Plant Mini Kit. Extraction protocol was the unmodified kit instructions. Five microsatellite primers (Karp_SB24, Karp_SB80, Karp_SB85, Karp_SB199, and Karp_SB201) derived for Salix burjatica Nasarow (Barker et al., 2003 [38]), that had previously shown polymorphism in S. sitchensis were used for screening samples in our collection. The polymerase chain reaction (PCR) procedure was 5 ng DNA, 25 ng forward primer, 25 ng reverse primer, 200 μm of each dNTP, 0.5 units of Taq DNA polymerase, 20 μm 10X buffer, and 1.5 mM MgCl2 (reagents from QIAGEN). Total reaction volume was 12.5 μL. Thermocycler conditions were 94°C/2 min; 35 cycles of: 94°C/40 s, 54°C/1 min, and 72°C/2 min, followed by 72°C/20 min [38]. Samples were run on a Beckman-Coulter CEQ800 to determine fragment/allele size. GENECLONE 2.0 [39] was used to determine the number of ramets with identical multilocus genotypes at five microsatellite loci. Clonal diversity was calculated according to Dorken and Eckert as follows [40]: where is the sample size and is the number of distinct genotypes (see also [39]).

3. Results

3.1. Field Survey of the Short-Term Reproductive History

During our study period, populations of S. sitchensis were comprised entirely of established mature individuals. Due to extensive flooding in 2006, 2007, and 2008, no seedlings, of any age, survived over winter in our study sites. Flooding covered all portions of our study sites for several weeks during seed release of S. sitchensis, resulting in no refugia for seedlings or seeds.

Five alive willow translocated stem fragments found at Island Park were not the result of beaver browse. At the remaining sites, ~10 observed fragments were not the result of beaver browse and none of these fragments was alive. Beaver herbivory was the predominant source of severed willow fragments in this study area during this time. Beaver cut willow fragments, both with adventitious roots and without, were observed at all sites and years except at Penny by the Fraser in 2007, and at Bowron River in 2008 (Table 2). In total, 42% of the cut stems produced adventitious roots the same year. Among the stems with roots marked in 2007, 20%, 0%, and 41% of the fragments recovered in 2008 were still displaying roots at Camp Creek, Bowron River, and Island Park, respectively (Table 2). Stems marked in 2007 and relocated in 2008 that were without roots had all died over winter. Marked stems not relocated were assumed to have been washed downstream.

3.2. Laboratory Analysis of Longer Reproductive History

We found 103 alleles for five microsatellite loci among all individuals sampled. All 286 sampled S. sitchensis individuals possessed a unique multilocus genotype and clonal diversity was therefore 1.0. More detailed information on the genetic structure of these populations can be found in Gerwing [35].

4. Discussion

As in our study, low survival of Salicaceae sexual propagules has been documented before [9, 11, 13, 22]. Abundant but variable levels of asexual propagules have also been observed before [11, 12, 14]. On the Drome River in France, 0.056 vegetative propagules per m2 were observed by Barsoum [11], which was within the range of values observed in central British Columbia. Higher survival rates of asexual propagules when compared to sexual propagules have been observed by Cottrell [9], Barsoum [11], and Moggridge and Gurnell [22]. The higher survival of asexual propagules may be due to carbohydrate reserves, preformed root, and shoot primordial [41], or the more resilient nature of branches versus seedlings. Survival of asexual propagules varied from 0–41% and this is most likely a result of variation in topography and the vertical structure of a site. None of the marked asexual propagules survived at Bowron River, which was a bare stretch of beach with few trees, stumps, bends, or other features to retain asexual propagules. It is likely that all of our propagules were washed downstream. Conversely, Island Park had the highest survival rate, and possessed the greatest physical heterogeneity. Island Park contained dense thickets of S. sitchensis, numerous piles of dead trees of various sizes, and abundant topographic variation. The abundance of features able to prevent asexual propagules from being washed downstream is most likely responsible for the higher survival rate of asexual propagules. Twenty percent of the asexual propagules survived at Camp Creek. Camp Creek had ample vertical structure in the form of dense thickets of grasses and sedges, dead trees, as well as numerous bends and elevation changes. While Camp Creek had more features able to retain asexual propagules than Bowron River, it possessed fewer features than Island Park. This is most likely the cause of Camp Creek’s asexual propagule survival being between that of Island Park and Bowron River. While features that are able to retain asexual propagules are most likely responsible for the observed variation in survival, factors such as soil type may also have played an important role. Future studies need to investigate the influence of other factors on the trends observed here.

The findings of the field and genetic analyses are, at first glance, contradictory. The field analysis suggested that the recent reproductive history of these populations was dominated by asexual reproduction, as only asexual propagules survived over winter. However, we did not find any mature individuals with shared genotypes. A large proportion of each population was genotyped, and the absence of shared genotypes suggests that asexual reproduction played a minor role in our populations. Stamati et al. [42] also found that field studies suggested the dominance of asexual reproduction; however, genetic analysis revealed that sexual reproduction dominated in populations of subarctic willows. Stamati et al. [42] postulated this may be due to continuous low levels of unobserved sexual reproduction, or sporadic bursts of sexual reproduction coinciding with optimal conditions for survival. Another possible explanation is that asexual propagules are important in colonization, while sexual propagules are responsible for population growth [1923]. If this was occurring, beavers would be crucial to this process, as nearly 100% of the asexual propagules in our study area were produced by beaver browse. Our study was not designed to differentiate between these hypotheses; however, all are possible and are not necessarily mutually exclusive. Further investigation into these processes is a crucial next step in this system.

Studies examining the influence of beaver browse on riparian plant reproduction have reported a variety of findings. Lin et al. [15] found numerous clonal genotypes in populations of Salix purpurea (Linnaeus). They attributed these clonal genotypes to vegetative reproduction, however, they were unable to differentiate between beaver-produced propagules, mechanically produced propagules, or anthropogenic propagules. Barsoum [11] found that 12% of all asexual propagules and 50% of coppice regrowth in Salix alba (Linnaeus) was due to beaver browse. The 12% of asexual propagules created by beaver foraging observed by Barsoum [11] is contrasted by the nearly 100% of vegetative propagules resulting from beaver foraging in our study. Barsoum [11] also stated that a large proportion of the vegetative propagules of Populus nigra observed on sandbars were a result of beaver placement. Jones et al. [43] observed that 35–40% of beaver-felled willows in Scotland showed regrowth. On the other hand, Rood et al. [12] excavated Populus saplings along the Elk River in British Columbia and found that only 7% of examined saplings grew from beaver browsed asexual propagules. Cooper et al. [14] also examined established vegetative propagules and concluded that beaver browse played no role in willow reproduction. The results of our study are in agreement with the findings of Rood et al. [12] and Cooper et al. [14]. While asexual propagules were abundant and no sexual propagules survived over winter, none of the mature trees sampled for the genetic analysis shared a genotype. The apparent dominance of sexual propagules, despite their poor survival over winter, is likely a result of sporadic peaks of sexual reproduction coinciding with optimal conditions for survival. Regardless, sexual reproduction appeared to dominate in the mature willow stands in this region, and beaver herbivory did not result in stands dominated by clones.

The inter-study disparity over the potential impacts of beaver herbivory on the reproduction of riparian vegetation may be based upon site variability. Differing levels of beaver activity, site topography, weather, and hydrological as well as historic events may have contributed to the observed variation. Different climatic conditions may greatly influence not only propagule survival, but beaver activity as well. Barsoum [11] also investigated a different species of beaver, the European beaver (Castor fiber Linnaeus), while our study looked at C. canadensis. The two species show differences in reproductive rates and dam construction [44, 45], potentially influencing beaver impacts on riparian populations. The finding of Cooper et al. [14], that beaver forage had no impact on willow reproduction, may be due to low levels of beaver activity. Cooper et al. [14] reported that few beaver dams were currently being built in their study area as the population had declined substantial over the past 60 years. The impact of beaver foraging on willow reproduction may have been minimal, as foraging activity would be slight or nonexistent. Differences in soil types may also influence the success of sexual and asexual propagules. Determining the influence of soil type on the trends observed in this study is a crucial next step for this system. Despite the variation within and between studies, beaver herbivory did not result in populations dominated by asexual propagules in our region. Beaver herbivory, in the form of willow asexual propagules, appeared to have minimal influence on the reproduction of S. sitchensis. Understanding the role beaver-produced asexual propagules may play in colonization and site modification is a crucial next step in understanding the influence of beaver herbivory on riparian plants.

More generally, this study adds to the body of literature on plant-herbivore interactions. Our results suggest that beavers played a minor role in willow reproduction in our study sites. However, other studies have postulated that beavers may play a significant role in riparian plant reproduction [9, 11]. This contradiction, combined with the potential role beaver-induced asexual propagules may play in site colonization and modification [1923], further complicates plant-herbivore interactions. As mentioned above, the nature of the plant-herbivore interaction depends on the species studied, the response variables examined, and the seasonal pattern of the foraging [18]. The results of our study, when compared with the literature, add to this list of variables by suggesting that plant-herbivore interactions may vary with species, habitat, successional stage, and method of reproduction. Care must be taken when extrapolating the findings of plant-herbivore interactions as substantial spatial and temporal variations exist within and between systems.

Acknowledgments

The authors thank Annika Klopp, Kathryn Berry, Brock Harpur, Ryan Thorgeirson, and Aydin Maxfield for their help with field work. This study was supported by the Natural Sciences and Engineering Research Council of Canada (NSERC) and Aleza Lake Research Forest (ALRF). They appreciate the helpful and insightful comments on an early draft of this paper by Dr. David Drolet and Dr. Jason Addison.