A History of Benthic Research in the NMFS Northeast Fisheries Science Center

The Northeast Fisheries Science Center of NOAA's National Marine Fisheries Service has a long history of research on benthic invertebrates and habitats in support of the management of living marine resources. These studies began in the 1870's under Spencer F. Baird's guidance as part of an effort to characterize the Nation's fisheries and living marine resources and their ecological interactions. This century and a quarter of research has included many benthic invertebrate studies, including community characterizations, shellfish biology and culture, pathology, ecosystem energy budget modeling, habitat evaluations, assessments of human impacts, toxic chemical bioaccumulation in demersal food webs, habitat or endangered species management, benthic autecology, systematics (to define new species and species population boundaries), and other benthic studies. Here we review the scope of past and current studies as a background for strategic research planning and suggest areas for further research to support NOAA's goals of sustainable fisheries management, healthy coastal ecosystems, and protected species populations.

A geo-referenced benthic habitat survey in support of natural resource management: Port Graham Bay, Alaska

The impact of recent changes in climate on the arctic environment and its ecosystems appear to have a dramatic affect on natural populations (National Research Council Committee on the Bering Sea Ecosystem 1996) and pose a serious threat to the continuity of indigenous arctic cultures that are dependent on natural resources for subsistence (Peterson D. L., Johnson 1995). In the northeast Pacific, winter storms have intensified and shifted southward causing fundamental changes in sea surface temperature patterns (Beamish 1993, Francis et al. 1998). Since the mid 1970’s surface waters of the central basin of the Gulf of Alaska (GOA) have warmed and freshened with a consequent increase in stratification and reduced winter entrainment of nutrients (Stabeno et al. 2004). Such physical changes in the structure of the ocean can rapidly affect lower trophic levels and indirectly affect fish and marine mammal populations through impacts on their prey (Benson and Trites 2002). Alaskan natives expect continued and perhaps accelerating changes in resources due to global warming (DFO 2006).and want to develop strategies to cope with their changing environment.

Ontogenetic shifts in natural diet during benthic stages of American lobster (Homarus americanus), off the Magdalen Islands

The natural diet of 506 American lobsters (Homarus americanus) ranging from instar V (4 mm cephalothorax length, CL) to the adult stage (112 mm CL) was determined by stomach content analysis for a site in the Magdalen Islands, Gulf of St. Lawrence, eastern Canada. Cluster and factor analyses determined four size groupings of lobsters based on their diet: <7.5 mm, 7.5 to <22.5 mm, 22.5 to <62.5 mm, and ≥62.5 mm CL. The ontogenetic shift in diet with increasing size of lobsters was especially apparent for the three dominant food items: the contribution of bivalves and animal tissue (flesh) to volume of stomach contents decreased from the smallest lobsters (28% and 39%, respectively) to the largest lobsters (2% and 11%, respectively), whereas the reverse trend was seen for rock crab Cancer irroratus (7% in smallest lobsters to 53% in largest lobsters). Large lobsters also ate larger rock crabs than did small lobsters.

The genetic diversity of Acipenser stellatus populations in the Caspian Sea was studied using microsatellite technique

A total of 361 caudal fin samples were collected from adult A. stellatus specimens caught in the north Caspian Sea, including specimens from Kazakhstan (Ural River), Russia (Volga River), Azerbaijan (Kura River), specimens caught in the south Caspian Sea including specimens from Fishery Zone 1 (from Astara to Anzali), Fishery Zone 2 (from Anzali to Ramsar), Fishery Zone 3 (from Nowshahr to Babolsar), Fishery Zone 4 (from Miyankaleh to Gomishan) as well as from specimens caught in Turkmenistan (all specimens were collected during the sturgeon stock assessment survey). About 2 g of fin tissue was removed from each caudal fin sample, stored in 96% ethyl alcohol and transferred to the genetic laboratory of the International Sturgeon Research Institute. Genomic DNA was extracted using phenol-chloroform method. The quality and quantity of DNA was assessed using 1% Agarose gel electrophoresis and Polymerase Chain Reaction (PCR) was conducted on the target DNA using 15 paired microsatellite primer. PCR products were electrophoresed on polyacrylamide gels (6%) that were stained using silver nitrate. Electrophoretic patterns and DNA bands were analyzed with BioCapt software. Allele count and frequency, genetic diversity, expected heterozygosity and observed heterozygosity allele number, and the effective allele number, genetic similarity and genetic distance, FST and RST were calculated. The Hardy Wienberg Equilibrium based on X2 and Analysis of Molecular Variance (AMOVA) at 10% confidence level was calculated using the Gene Alex software. Dendrogram for genetic distances and identities were calculated using TFPGA program for any level of the hierarchy. It is evident from the results obtained that the 15 paired primers studied, polymorphism was observed in 10 pairs in 12 loci, while one locus did not produce DNA bands. Mean allele number was 13.6. Mean observed and expected heterozygosity was 0.86 and 0.642, respectively. It was also seen that specimens from all regions were not in Hardy Wienberg Equilibrium in most of the loci (P≤0.001). Highest Fst (0.063) was observed when comparing specimens from Fishery Zone 2 and Fishery Zone 4 (Nm=3.7) and lowest FST (0.028) was observed when comparing specimens from the Volga River and those from the Ural River (8.7). Significant differences (P<0.01) were observed between RST recorded in the specimens studied. Highest genetic distance (0.604) and lowest genetic resemblance (0.547) were observed between specimens from Fishery zones 2 and 4. Lowest genetic distance (0.311) and highest genetic resemblance (0.733) was observed between specimens from Turkmenistan and specimens from Fishery zone 1. Based on the genetic dendrogeram tree derived by applying UPGMA algorithm, A. stellatus specimens from Fishery zone 2 or in other words specimens from the Sepidrud River belong to one cluster which divides into two clusters, one of which includes specimens from Fishery zones 1, 3 and 4 and specimens from Turkmenistan while the other cluster includes specimens from Ural, Volga and Kura Rivers. It is thus evident that the main population of this species belongs to the Sepidrud River. Results obtained from the present study show that at least eight different populations of A. stellatus are found in the north and south Caspian Sea, four of which are known populations including the Ural River population, the Volga River population, the Kura River population and the Sepidrud River populations. The four other populations identified belonging to Fishery zones 1, 3, and 4 and to Turkmenistan are most probably late or early spawners of the spring run and autumn run of each of the major rivers mentioned. Specific markers were also identified for each of the populations identified. The Ural River population can be identified using primers Spl-68, 54b and Spl-104, 163 170, 173, the Volga River population can be identified using primers LS-54b and Spl-104, 170, 173 113a and similarly the population from the Kura River can be identified using primers LS-34, 54b and Spl-163, 173 and that from the Sepidrud River can be identified using primers LS-19, 34, 54b and Spl-105, 113b. This study gives evidence of the presence of different populations of this species and calls for serious measures to be taken to protect the genetic stocks of these populations. Considering that the population of A. stellatus in Fishery zone 2 is an independent population of the Sepidrud River in the Gilan Province, the catch of these fishes in the region needs to be controlled and regulated in order to restore the declining stocks of this species.

The genetic structure and phylogenetics of pikeperch (Sander lucioperca) Anzali and Amirkolaye wetlands and perch (Perca fluviatilis) populations in Aras Dam and south-west of the Caspian Sea

The genetic structure of pikeperch (Sander lucioperca) and perch (Perca fluviatilis) populations was studied using microsatellite technique. A total of 207 specimens of adult pikeperch were collected from Aras dam (57 specimens), Anzali wetland (50 specimens), Talesh (50 specimens) and Chaboksar (50 specimens) coasts. Also a total of 158 specimens of adult perch were collected from Anzali (Abkenar (50 specimens)and Hendekhale(48 specimens)) and Amirkolaye(60 specimens) wetlands. About 2 g of each specimen's dorsal fin was removed, stored in 96% ethyl alcohol and transferred to the genetic laboratory of the International Sturgeon Research Institute. Genomic DNA was extracted using ammonium-acetate method. The quality and quantity of DNA was assessed using 1% agarose gel electrophoresis. Polymerase Chain Reaction (PCR) was conducted on the target DNA using 15 pairs of microsatellite primers. PCR products were electrophoresed on poly acryl amide gels (6%) that were stained that were stained using silver nitrate. DNA bands were analyzed with BioCapt software. Allele count and frequency, genetic diversity, expected and observed heterozygosity , allele number and the effective allele number, genetic similarity and genetic distance, Fst, Rst, Hardy Weinberg Equilibrium based on X2 and Analysis of Molecular Variance (AMOVA) at 10% confidence level was calculated using the Gene Alex software. Dendogram for genetic distances and identities were calculated using TFPGA program for any level of hierarchy. The results for P. fluviatilis showed that from 15 pair of primers that were examined 6 polymorphic and 7 monomorphic loci were produced, while 2 loci didn't produce any DNA bands. Mean allele number was 4.1±1.1 and mean observed and expected heterozygosity was 0.56±0.12 and 0.58±0.14 respectively. It was also seen that specimens from all regions were not in Hardy Weinberg Equilibrium in some of loci (P<0.001). Highest Fst (0.095) with Nm=2.37 was observed between Hendekhale and Amirkolaye and the lowest Fst (0.004) with Nm=59.31 was observed between Abkenar and Hendekhale. According to AMOVA Significant difference (P<0.05) was observed between recorded Rst in the studied regions in Anzali and Amirkolaye lagoons. In another words there are two distinct populations of this species in Anzali and Amirkolaye lagoons. The highest genetic distance (0.181) and lowest genetic resemblance (0.834) were observed between specimens from Hendekhale and Amirkolaye and the lowest genetic distance (0.099) and highest genetic 176 resemblance (0.981) were observed between specimens from Abkenar and Hendekhale. Based on the genetic dendogram tree derived by applying UPGMA algorithm, specimens from Anzali and Amirkolaye wetlands have the same ancestor. On the other hand there is no noticeable genetic distance between the specimens of these two regions. Also the results for S. lucioperca showed that from 15 pair of primers that were examined 6 polymorphic and 7 monomorphic loci were produced, while 2 loci didn't produce any DNA bands. Mean allele number was 3.0±0.6 and mean observed and expected heterozygosity was 0.52±0.21 and 0.50±0.14 respectively. It was also seen that specimens from all regions were not in Hardy Weinberg Equilibrium in some of loci (P<0.001). Highest Fst (0.093) with Nm=2.43 was observed between Aras dam and Anzali wetland and the lowest Fst (0.022) with Nm=11.27 was observed between Talesh and Chaboksar coasts. Significant differences (P<0.05) were observed between recorded Rst in the studied regions exept for Talesh and Chaboksar Coasts. In another words there are three distinct populations of this species in Caspian sea, Anzali wetland and Aras dam. Highest genetic distance (0.110) and lowest genetic resemblance (0.896) were observed between specimens from Aras dam and Anzali wetland and the lowest genetic distance (0.034) and highest genetic resemblance (0.966) were observed between specimens from Talesh and Chaboksar coasts. Based on the genetic dendogram tree derived by applying UPGMA algorithm, specimens from Talesh and Chaboksar coasts have the lowest genetic distance. On the other hand the main population of this species belongs to Anzali wetland. Phylogenetic relationship of these two species was inferred using mitochondrial cytochrome b gene sequencing. For this purpose 2 specimens of P. fluviatilis from Anzali wetland, 2 specimens of S. lucioperca from Aras dam and 2 specimens of S. lucioperca from Anzali wetland were sequenced and submitted in Gene Bank. These sequences were aligned with Clustal W. The phylogenic relationships were assessed with Mega 4. The results of evolutionary history studies of these species using Neighbor-Joining and Maximum Parsimony methods showed that the evolutionary origin of pikeperch in Aras Dam and Anzali wetland is common. On the other hand these two species had common ancestor in about 4 million years ago. Also different sequences of any region specimens are supposed as different haplotypes. 177 As a conclusion the results of this study showed that microsatellite and mtDNA sequencing methods respectively are effective in genetic structure and phylogenic studies of P. fluviatilis and S. lucioperca.

Isolation of Y- and X-linked SCAR markers in yellow catfish and application in the production of all-male populations

P>Sex controls have been performed in some farmed fish species because of significant growth differences between females and males. In yellow catfish (Pelteobagrus fulvidraco), adult males are three times larger than female adults. In this study, six Y- and X-linked amplified fragment length polymorphism fragments were screened by sex-genotype pool bulked segregant analysis and individual screening. Interestingly, sequence analysis identified two pairs of allelic genes, Pf33 and Pf62. Furthermore, the cloned flanking sequences revealed several Y- and X-specific polymorphisms, and four Y-linked or X-linked sequence characterized amplified region (SCAR) primer pairs were designed and converted into Y- and X-linked SCAR markers. Consequently, these markers were successfully used to identify genetic sex and YY super-males, and applied to all-male population production. Thus, we developed a novel and simple technique to help commercial production of YY super-males and all-male populations in the yellow catfish.

Genetic differentiation in populations of the cestode Bothriocephalus acheilognathi (Cestoda, Pseudophyllidea) as revealed by eight microsatellite markers

The genetic structure of populations of the fish cestode, Bothriocephalus acheilognathi collected from Bailianhe Reservoir (BLH), Changshou (CSH) and Liangzi (LZH) Lakes was investigated by using 8 microsatellite loci. A total of 108 adult worms were genotyped at each of the 8 loci. For the 3 populations, the mean number of alleles per locus ranged from 2.38 to 5.5, and the mean expected heterozygosity ranged from 0.432 to 0.559. The average polymorphic information content (PIC) was from 0.384 to 0.492. The significant F-is values indicated non-random mating within LZH and BLH populations. On the other hand, when samples were further classified into subpopulations at the level of host fish species, no or little heterozygote deficiency was detected at most loci, showing that cross-fertilization, predominantly, but not exclusively, must have occurred within the subpopulations. Microsatellite markers also revealed an unexpected high level of genetic differentiation, as measured by R-st and N-m values or by deltau(2) genetic distance among subpopulations from different hosts. Factors influencing the population genetic structure and the parasite host specificity are discussed.

The role of the IL-1 type 1 receptor in the life, death and differentiation of adult hippocampal neural precursor cells

Neurogenesis occurs in two distinct regions of the adult brain; the subgranular zone (SGZ) of the dentate gyrus (DG) in the hippocampus, and the subventricular zone (SVZ) lining the lateral ventricles. It is now well-known that adult hippocampal neurogenesis can be modulated by a number of intrinsic and extrinsic factors e.g. local signalling molecules, exercise, environmental enrichment and learning. Moreover, levels of adult hippocampal neurogenesis decrease with age, at least in rodents, and alterations in hippocampal neurogenesis have been reported in animal models and human studies of neuropsychiatric and neurodegenerative conditions. Neuroinflammation is a common pathological feature of these conditions and is also a potent modulator of adult hippocampal neurogenesis. Recently, the orphan nuclear receptor TLX has been identified as an important regulator of adult hippocampal neurogenesis as its expression is necessary to maintain the neural precursor cell (NPC) pool in the adult DG. Likewise, exposure of animals to voluntary exercise has been consistently demonstrated to promote adult hippocampal neurogenesis. Lentivirus (LV)- mediated gene transfer is a useful tool to elucidate gene function and to explore potential therapeutic candidates across an array of conditions as it facilitates sustained gene expression in both dividing and post-mitotic cell populations. Both intrinsic and extrinsic factors are important regulators of adult hippocampal neurogenesis. Examining how these factors are affected by an inflammatory stimulus, and the subsequent effects on adult hippocampal neurogenesis provides important information for the development of novel treatment strategies for neuropsychiatric and neurodegenerative conditions in which adult hippocampal neurogenesis is impaired. The aims of the series of experiments presented in this thesis were to examine the effect of the pro-inflammatory cytokine interleukin-1β (IL-1β) on adult hippocampal NPCs both in vitro and in vivo. In vitro, we have shown that IL-1β reduces proliferation of adult hippocampal NPCs in a dose and time-dependent manner. In addition, we have demonstrated that TLX expression is reduced by IL-1β. Blockade of IL-1β signalling prevented both the IL-1β-induced reduction in cell proliferation and TLX expression. In vivo, we examined the effect of short term and long term exposure to LV-IL-1β in sedentary mice and in mice exposed to voluntary running. We demonstrated that impaired hippocampal neurogenesis is only evident after long term exposure to IL-1β. In mice exposed to voluntary running, hippocampal neurogenesis is significantly increased following short-term but not long-term exposure to running. Moreover, short-term running effectively prevents any IL-1β-induced effects on hippocampal neurogenesis; however, no such effects are seen following long-term exposure to running.

Low demographic variability in wild primate populations: fitness impacts of variation, covariation, and serial correlation in vital rates.

In a stochastic environment, long-term fitness can be influenced by variation, covariation, and serial correlation in vital rates (survival and fertility). Yet no study of an animal population has parsed the contributions of these three aspects of variability to long-term fitness. We do so using a unique database that includes complete life-history information for wild-living individuals of seven primate species that have been the subjects of long-term (22-45 years) behavioral studies. Overall, the estimated levels of vital rate variation had only minor effects on long-term fitness, and the effects of vital rate covariation and serial correlation were even weaker. To explore why, we compared estimated variances of adult survival in primates with values for other vertebrates in the literature and found that adult survival is significantly less variable in primates than it is in the other vertebrates. Finally, we tested the prediction that adult survival, because it more strongly influences fitness in a constant environment, will be less variable than newborn survival, and we found only mixed support for the prediction. Our results suggest that wild primates may be buffered against detrimental fitness effects of environmental stochasticity by their highly developed cognitive abilities, social networks, and broad, flexible diets.

Independent evolutionary origins of landlocked alewife populations and rapid parallel evolution of phenotypic traits.

Alewife, Alosa pseudoharengus, populations occur in two discrete life-history variants, an anadromous form and a landlocked (freshwater resident) form. Landlocked populations display a consistent pattern of life-history divergence from anadromous populations, including earlier age at maturity, smaller adult body size, and reduced fecundity. In Connecticut (USA), dams constructed on coastal streams separate anadromous spawning runs from lake-resident landlocked populations. Here, we used sequence data from the mtDNA control region and allele frequency data from five microsatellite loci to ask whether coastal Connecticut landlocked alewife populations are independently evolved from anadromous populations or whether they share a common freshwater ancestor. We then used microsatellite data to estimate the timing of the divergence between anadromous and landlocked populations. Finally, we examined anadromous and landlocked populations for divergence in foraging morphology and used divergence time estimates to calculate the rate of evolution for foraging traits. Our results indicate that landlocked populations have evolved multiple times independently. Tests of population divergence and estimates of gene flow show that landlocked populations are genetically isolated, whereas anadromous populations exchange genes. These results support a 'phylogenetic raceme' model of landlocked alewife divergence, with anadromous populations forming an ancestral core from which landlocked populations independently diverged. Divergence time estimates suggest that landlocked populations diverged from a common anadromous ancestor no longer than 5000 years ago and perhaps as recently as 300 years ago, depending on the microsatellite mutation rate assumed. Examination of foraging traits reveals landlocked populations to have significantly narrower gapes and smaller gill raker spacings than anadromous populations, suggesting that they are adapted to foraging on smaller prey items. Estimates of evolutionary rates (in haldanes) indicate rapid evolution of foraging traits, possibly in response to changes in available resources.

Imaging Learned Song Representations in Populations of Sensorimotor Neurons Essential to Vocal Communication

Perceiving or producing complex vocalizations such as speech and birdsongs require the coordinated activity of neuronal populations, and these activity patterns can vary over space and time. How learned communication signals are represented by populations of sensorimotor neurons essential to vocal perception and production remains poorly understood. Using a combination of two-photon calcium imaging, intracellular electrophysiological recording and retrograde tracing methods in anesthetized adult male zebra finches (Taeniopygia guttata), I addressed how the bird's own song and its component syllables are represented by the spatiotemporal patterns of activity of two spatially intermingled populations of projection neurons (PNs) in HVC, a sensorimotor area required for song perception and production. These experiments revealed that neighboring PNs can respond at markedly different times to song playback and that different syllables activate spatially intermingled HVC PNs within a small region. Moreover, noise correlation analysis reveals enhanced functional connectivity between PNs that respond most strongly to the same syllable and also provides evidence of a spatial gradient of functional connectivity specific to PNs that project to song motor nucleus (i.e. HVC_RA cells). These findings support a model in which syllabic and temporal features of song are represented by spatially intermingled PNs functionally organized into cell- and syllable-type networks.

Pattern of benthic biomass size spectra from shallow waters in the East China Seas

Benthic biomass size spectra (BSS) and normalized biomass size spectra were constructed, and benthic secondary production was estimated by a size spectrum equation in the shallow waters in the East China Sea, ranging latitudinally from 40A degrees N to 29A degrees N. The BSS patterns were bimodal, two biomass peaks corresponding to meiofauna and macrofauna, respectively, separated by a trough of low biomass at 8-256 mu g individual dry weight which varied in position with median sediment particle size. The BSS also displayed bimodality within meiofauna size ranges, which in most stations was due to the relative proportions of nematodes and other meiofauna taxa. Re-analysis of data from sites in the UK, South Africa, and Antarctic showed a similar bimodality in the adult species body size distribution within the meiofauna size range. Macrofaunal production estimated by the size spectrum equation was very similar to the results of Brey90 empirical equation. However, these production values were much lower than those calculated by Brey01. Different individual dry-to-wet conversion ratios, temperature deviation, and macrofauna taxonomic composition might be responsible for the between-model differences. The macrofaunal P/B ratios calculated by this equation ranged from 0.3 to 3.4 which were in accordance with values from Northern Hemisphere mid-latitudes. Meiofaunal production estimates will need further empirical support.

Models of open populations with space-limited recruitment: extension of theory and application to the barnacle Chthamalus montagui.

1. Barnacles are a good model organism for the study of open populations with space-limited recruitment. These models are applicable to other species with open supply of new individuals and resource limitation. The inclusion of space in models leads to reductions in recruitment with increasing density, and thus predictions of population size and stability are possible. 2. Despite the potential generality of a demographic theory for open space-limited populations, the models currently have a narrow empirical base. In this study, a model for an open population with space-limited recruitment was extended to include size-specific survival and promotions to any size class. The assumptions of this model were tested using data from a pan-European study of the barnacle Chthamalus montagui Southward. Two models were constructed: a 6-month model and a periodic annual model. Predicted equilibria and their stabilities were compared between shores. 3. Tests of model assumptions supported the extension of the theory to include promotions to any size class. Mortality was found to be size-specific and density independent. Studied populations were open, with recruitment proportional to free space. 4. The 6-month model showed a significant interaction between time and location for equilibrium free space. This may have been due to contrasts in the timing of structuring processes (i.e. creating and filling space) between Mediterranean and Atlantic systems. Integration of the 6-month models into a periodic annual model removed the differences in equilibrium-free space between locations. 5. Model predictions show a remarkable similarity between shores at a European scale. Populations were persistent and all solutions were stable. This reflects the apparent absence of density-dependent mortality and a high adult survivorship in C. montagui. As populations are intrinsically stable, observations of fluctuations in density are directly attributable to variations in the environmental forcing of recruitment or mortality

Do non-native invasive fish support elevated lamprey populations?

Summary
1: Managing populations of predators and their prey to achieve conservation or resource management goals is usually technically challenging and frequently socially controversial. This is true even in the simplest ecosystems but can be made much worse when predator–prey relationships are in?uenced by complex interactions, such as biological invasions, population trends or animal movements.
2: Lough Neagh in Northern Ireland is a European stronghold for pollan Coregonus autumnalis, a coregonine ?sh and for river lampreyLampetra ?uviatilis, which feeds parasitically as an adult. Both species are of high conservation importance. Lampreys are known to consume pollan but detailed knowledge of their interactions is scant. While pollan is well known to be a landlocked species in Ireland, the life cycle of normally anadromous river lamprey in Lough Neagh has been unclear. The Lough is also a highly perturbed ecosystem, supporting several invasive, non-native ?sh species that have the potential to in?uence lamprey–pollan interactions.
3: We applied stable isotope techniques to resolve both the movement patterns of lamprey and trophic interactions in this complex community. Recognizing that stable isotope studies are often hampered by high-levels of variability and uncertainty in the systems of interest, we employed novel Bayesian mixing models, which incorporate variability and uncertainty.
4: Stable isotope analyses identi?ed troutSalmo trutta and non-native breamAbramis brama as the main items in lamprey diet. Pollan only represented a major food source for lamprey between May and July.
5: Stable isotope ratios of carbon in tissues from 71 adult lamprey showed no evidence of marine carbon sources, strongly suggesting that Lough Neagh is host to a highly unusual, nonanadromous freshwater population. This ?nding marks out the Lough’s lamprey population as of particular scienti?c interest and enhances the conservation signi?cance of this feature of the Lough.
6: Synthesis and applications.Our Bayesian isotopic mixing models illustrate an unusual pattern of animal movement, enhancing conservation interest in an already threatened population. We have also revealed a complex relationship between lamprey and their food species that is suggestive of hyperpredation, whereby non-native species may sustain high lamprey populations that may in turn be detrimental to native pollan.Long-term conservation of lamprey and pollan in this system is likely to require management intervention, but in light of this exceptional complexity, no simple management options are currently supported. Conservation plans will require better characterization ofpopulation-level interactions and simulation modelling of interventions. More generally, our study demonstrates the importance of considering a full range of possible trophic interactions, particularly in complex ecosystems, and highlights Bayesian isotopic mixing models as powerful tools in resolving trophic relationships.
Key-words: Bayesian, conservation dilemma, Coregonus autumnalis, hyperpredation, Lampetra ?uviatilis, pollan, potamodromous, River lamprey, stable isotope analysis in R, stable isotope