Molecular pedigree reconstruction and estimation of evolutionary parameters in a wild Atlantic salmon river system with incomplete sampling: a power analysis

Background: Pedigree reconstruction using genetic analysis provides a useful means to estimate fundamental population biology parameters relating to population demography, trait heritability and individual fitness when combined with other sources of data. However, there remain limitations to pedigree reconstruction in wild populations, particularly in systems where parent-offspring relationships cannot be directly observed, there is incomplete sampling of individuals, or molecular parentage inference relies on low quality DNA from archived material. While much can still be inferred from incomplete or sparse pedigrees, it is crucial to evaluate the quality and power of available genetic information a priori to testing specific biological hypotheses. Here, we used microsatellite markers to reconstruct a multi-generation pedigree of wild Atlantic salmon (Salmo salar L.) using archived scale samples collected with a total trapping system within a river over a 10 year period. Using a simulation-based approach, we determined the optimal microsatellite marker number for accurate parentage assignment, and evaluated the power of the resulting partial pedigree to investigate important evolutionary and quantitative genetic characteristics of salmon in the system.

Results: We show that at least 20 microsatellites (ave. 12 alleles/locus) are required to maximise parentage assignment and to improve the power to estimate reproductive success and heritability in this study system. We also show that 1.5 fold differences can be detected between groups simulated to have differing reproductive success, and that it is possible to detect moderate heritability values for continuous traits (h(2) similar to 0.40) with more than 80% power when using 28 moderately to highly polymorphic markers.

Conclusion: The methodologies and work flow described provide a robust approach for evaluating archived samples for pedigree-based research, even where only a proportion of the total population is sampled. The results demonstrate the feasibility of pedigree-based studies to address challenging ecological and evolutionary questions in free-living populations, where genealogies can be traced only using molecular tools, and that significant increases in pedigree assignment power can be achieved by using higher numbers of markers.

Complex pattern of genetic structuring in the Atlantic salmon (Salmo salar L.) of the River Foyle system in northwest Ireland: disentangling the evolutionary signal from population stochasticity

Little is known about the microevolutionary processes shaping within river population genetic structure of aquatic organisms characterized by high levels of homing and spawning site fidelity. Using a microsatellite panel, we observed complex and highly significant levels of intrariver population genetic substructure and Isolation-by-Distance, in the Atlantic salmon stock of a large river system. Two evolutionary models have been considered explaining mechanisms promoting genetic substructuring in Atlantic salmon, the member-vagrant and metapopulation models. We show that both models can be simultaneously used to explain patterns and levels of population structuring within the Foyle system. We show that anthropogenic factors have had a large influence on contemporary population structure observed. In an analytical development, we found that the frequently used estimator of genetic differentiation, F-ST, routinely underestimated genetic differentiation by a factor three to four compared to the equivalent statistic Jost's D-est (Jost 2008). These statistics also showed a near-perfect correlation. Despite ongoing discussions regarding the usefulness of "adjusted" F-ST statistics, we argue that these could be useful to identify and quantify qualitative differences between populations, which are important from management and conservation perspectives as an indicator of existence of biologically significant variation among tributary populations or a warning of critical environmental damage.

The replacement of a native freshwater amphipod by an invader: roles for environmental degradation and intraguild predation

We assessed the extent to which an invader, Gammarus pulex (Crustacea: Amphipoda), has replaced a native, Gammarus duebeni celticus, over a 13-year period in a European river system and some of the abiotic and biotic factors that could account for this. Between 1988 and 2001, 56% of mixed-species sites had become invader-only sites, whereas no mixed sites had become native only again. The native dominated areas of higher dissolved oxygen and water quality, with the reciprocal true for the invader. Field transplant experiments revealed that native survivorship was lower in areas where it had been replaced than in areas where the invader does not yet occur. In invader-only areas, native survivorship was lower than that of the invader when kept separately and lowest when both species were kept together. We also observed predation of the native by the invader. Laboratory oxygen manipulation experiments revealed that at 30% saturation, the native's survivorship was two thirds that of the invader. We conclude that decreasing water quality favours replacement of the native by the invader.

Limitations of instantaneous water quality sampling in surface-water catchments: Comparison with near-continuous phosphorus time-series data.

The validity of load estimates from intermittent, instantaneous grab sampling is dependent on adequate spatial coverage by monitoring networks and a sampling frequency that re?ects the variability in the system under study. Catchments with a ?ashy hydrology due to surface runoff pose a particular challenge as intense short duration rainfall events may account for a signi?cant portion of the total diffuse transfer of pollution from soil to water in any hydrological year. This can also be exacerbated by the presence of strong background pollution signals from point sources during low flows. In this paper, a range of sampling methodologies and load estimation techniques are applied to phosphorus data from such a surface water dominated river system, instrumented at three sub-catchments (ranging from 3 to 5 km2 in area) with near-continuous monitoring stations. Systematic and Monte Carlo approaches were applied to simulate grab sampling using multiple strategies and to calculate an estimated load, Le based on established load estimation methods. Comparison with the actual load, Lt, revealed signi?cant average underestimation, of up to 60%, and high variability for all feasible sampling approaches. Further analysis of the time series provides an insight into these observations; revealing peak frequencies and power-law scaling in the distributions of P concentration, discharge and load associated with surface runoff and background transfers. Results indicate that only near-continuous monitoring that re?ects the rapid temporal changes in these river systems is adequate for comparative monitoring and evaluation purposes. While the implications of this analysis may be more tenable to small scale ?ashy systems, this represents an appropriate scale in terms of evaluating catchment mitigation strategies such as agri-environmental policies for managing diffuse P transfers in complex landscapes.

The geoarchaeology of Paleolithic rivers of southwest Britain

In Britain, the majority of Lower and Middle Paleolithic archaeological finds come from river terrace deposits. The impressive “staircase” terrace sequences of southeast England, and research facilitated by aggregate extraction have provided a considerable body of knowledge about the terrace chronology and associated archaeology in that area. Such research has been essential in considering rates of uplift, climatic cycles, archaeological chronologies, and the landscapes in which hominins lived. It has also promoted the view that southeast England was a major hominin route into Britain. By contrast, the terrace deposits of the southwest have been little studied. The Palaeolithic Rivers of South West Britain (PRoSWEB) project employed a range of geoarchaeological methodologies to address similar questions at different scales, focusing on the rivers Exe, Axe, Otter, and the paleo-Doniford, all of which were located south of the maximum Pleistocene glacial limit (marine oxygen isotope stage [MIS] 4–2). Preliminary analysis of the fieldwork results suggests that although the evolution of these catchments is complex, most conform to a standard staircase-type model, with the exception of the Axe, and, to a lesser extent, the paleo-Doniford, which are anomalous. Although the terrace deposits are less extensive than in southeast Britain, differentiation between terraces does exist, and new dates show that some of these terraces are of great antiquity (MIS 10+). The project also reexamined the distribution of artifacts in the region and confirms the distributional bias to the river valleys, and particularly the rivers draining southward to the paleo–Channel River system. This distribution is consistent with a model of periodic occupation of the British peninsula along and up the major river valleys from the paleo–Channel River corridor. These data have a direct impact on our understanding of the paleolandscapes of the southwest region, and therefore our interpretations of the Paleolithic occupation of the edge of the continental landmass.

Habitat-specific adaptation of immune responses of stickleback (Gasterosteus aculeatus) lake and river ecotypes

Freshwater populations of three-spined sticklebacks (Gasterosteus aculeatus) in northern Germany are found as distinct lake and river ecotypes. Adaptation to habitat-specific parasites might influence immune capabilities of stickleback ecotypes. Here, naive laboratory-bred sticklebacks from lake and river populations were exposed reciprocally to parasite environments in a lake and a river habitat. Sticklebacks exposed to lake conditions were infected with higher numbers of parasite species when compared with the river. River sticklebacks in the lake had higher parasite loads than lake sticklebacks in the same habitat. Respiratory burst, granulocyte counts and lymphocyte proliferation of head kidney leucocytes were increased in river sticklebacks exposed to lake when compared with river conditions. Although river sticklebacks exposed to lake conditions showed elevated activation of their immune system, parasites could not be diminished as effectively as by lake sticklebacks in their native habitat. River sticklebacks seem to have reduced their immune-competence potential due to lower parasite diversity in rivers

Technical Note: Assessing a 24/7 solution for monitoring water quality loads in small river catchments

Quantifying nutrient and sediment loads in catchments is dif?cult owing to diffuse controls related to storm hydrology. Coarse sampling and interpolation methods are prone to very high uncertainties due to under-representation of high discharge, short duration events. Additionally, important low-?ow processes such as diurnal signals linked to point source impacts are missed. Here we demonstrate a solution based on a time-integrated approach to sampling with a standard 24 bottle autosampler con?gured to take a sample every 7 h over a week according to a Plynlimon design. This is evaluated with a number of other sampling strategies using a two-year dataset of sub-hourly discharge and phosphorus concentration data. The 24/7 solution is shown to be among the least uncertain in estimating load (inter-quartile range: 96% to 110% of actual load in year 1 and 97% to 104% in year 2) due to the increased frequency raising the probability of sampling storm events and point source signals. The 24/7 solution would appear to be most parsimonious in terms of data coverage and certainty, process signal representation, potential laboratory commitment, technology requirements and the ability to be widely deployed in complex catchments.

UROTENSIN II FROM THE RIVER LAMPREY (LAMPETRA-FLUVIATILIS), THE SEA LAMPREY (PETROMYZON-MARINUS), AND THE PADDLEFISH (POLYODON SPATHULA)

Urotensin II was isolated from extracts of the whole brain of the river lamprey (Lampetra fluviatilis) and the sea lamprey (Petromyzon marinus). The primary structure of the peptide from both species is the same (Asn-Asn-Phe-Ser-Asp-Cys-Phe-Trp-Lys-Tyr-Cys-Val) and this amino acid sequence is identical to that of urotensin II from the dogfish and skate. Consistent with previous morphological studies indicating that the Agnatha lack a caudal neurosecretory system, urotensin II was not detected in an extract of P. marinus spinal cord. The data suggest that the urotensin II may have functioned in the earliest vertebrates as a neurotransmitter/neuromodulator in the central nervous system rather than as a neurohormone of the caudal neurosecretory system. Urotensin II was also isolated from an extract of the spinal cord of a chondrostean fish, the paddlefish (Polyodon spathula). The primary structure of the paddlefish urotensin II (Gly-Ser-Thr-Ser-Glu-Cys-Phe-Trp-Lys-Tyr-Cys-Val) is the same as that of another chondrostean, the sturgeon (Acipenser ruthenus). The study provides further evidence for a widespread distribution of urotensin II in vertebrate species and suggests that the primary structure of the peptide is better conserved in these phylogenetically ancient fish than in teleosts. (C) 1995 Academic Press, Inc.

Modelling the fate of persistent organic pollutants (POPs) in the North Sea system

The environmental fate of selected persistent organic pollutants (POPs) in the North Sea system is modelled with a high resolution Fate and Transport Ocean Model (FANTOM) that uses hydrodynamic model output from the Hamburg Shelf Ocean Model (HAMSOM). Large amounts of POPs enter the North Sea from the surrounding highly populated, industrialised and agricultural countries. Major pathways to the North Sea are atmospheric deposition and river inputs, with additional contributions coming from bottom sediments and adjacent seas. The model domain covers the entire North Sea region, extending northward as far as the Shetland Islands, and includes adjacent basins such as the Skagerrak, Kattegat, and the westernmost part of the Baltic Sea. Model resolution (for both models) is 1.5’ latitude x 2.5’ longitude (approximately 3 km horizontal resolution) with 30 vertical levels. The POP model also has 20 sediment layers. Important model processes controlling the fate of POPs in the North Sea system are discussed. Results focus on Lindane gamma- HCH or gamma-hexachlorocyclohexane) and PCB 153.

Transatlantic distribution of the Alaskan White River Ash

Volcanic ash layers preserved within the geologic record represent precise time markers that correlate disparate depositional environments and enable the investigation of synchronous and/or asynchronous behaviors in Earth system and archaeological sciences. However, it is generally assumed that only exceptionally powerful events, such as supereruptions (≥450 km3 of ejecta as dense-rock equivalent; recurrence interval of ∼105 yr), distribute ash broadly enough to have an impact on human society, or allow us to address geologic, climatic, and cultural questions on an intercontinental scale. Here we use geochemical, age, and morphological evidence to show that the Alaskan White River Ash (eastern lobe; A.D. 833–850) correlates to the “AD860B” ash (A.D. 846–848) found in Greenland and northern Europe. These occurrences represent the distribution of an ash over 7000 km, linking marine, terrestrial, and ice-core records. Our results indicate that tephra from more moderate-size eruptions, with recurrence intervals of ∼100 yr, can have substantially greater distributions than previously thought, with direct implications for volcanic dispersal studies, correlation of widely distributed proxy records, and volcanic hazard assessment.

Just a matter of time: Fungi and roots significantly and rapidly aggregate soil over four decades along the Tagliamento River, NE Italy

Fluvial islands are emergent landforms which form at the interface between the permanently inundated areas of the river channel and the more stable areas of the floodplain as a result of interactions between physical river processes, wood and riparian vegetation. These highly dynamical systems are ideal to study soil structure development in the short to medium term, a process in which soil biota and plants play a substantial role. We investigated soil structure development on islands along a 40 year chronosequence within a 3 km island-braided reach of the Tagliamento River, Northeastern Italy. We used several parameters to capture different aspects of the soil structure, and measured biotic (e.g., fungal and plant root parameters) and abiotic (e.g. organic carbon) factors expected to determine the structure. We estimated models relating soil structure to its determinants, and, in order to confer statistical robustness to our results, we explicitly took into account spatial autocorrelation, which is present due to the space for time substitution inherent in the study of chronosequences and may have confounded results of previous studies. We found that, despite the eroding forces from the hydrological and geomorphological dynamics to which the system is subject, all soil structure variables significantly, and in some case greatly increased with site age. We interpret this as a macroscopic proxy for the major direct and indirect binding effects exerted by root variables and extraradical hyphae of arbuscular mycorrhizal fungi (AMF). Key soil structure parameters such as percentage of water stable aggregates (WSA) can double from the time the island landform is initiated (mean WSA = 30%) to the full 40 years (mean WSA = 64%) covered by our chronosequence. The study demonstrates the fundamental role of soil biota and plant roots in aggregating soils even in a system in which intense short to medium term physical disturbances are common.

14C as a tool to trace terrestrial carbon in a complex lake system: implications for food-web structure and carbon cycling

Globally lakes bury and remineralise significant quantities of terrestrial C, and the associated flux of terrestrial C strongly influences their functioning. Changing deposition chemistry, land use and climate induced impacts on hydrology will affect soil biogeochemistry and terrestrial C export¹ and hence lake ecology with potential feedbacks for regional and global C cycling. C and nitrogen stable isotope analysis (SIA) has identified the terrestrial subsidy of freshwater food webs. The approach relies on different ¹³C fractionation in aquatic and terrestrial primary producers, but also that inorganic C demands of aquatic primary producers are partly met by ¹³C depleted C from respiration of terrestrial C, and ‘old’ C derived from weathering of catchment geology. SIA thus fails to differentiate between the contributions of old and recently fixed terrestrial C. Natural abundance ¹⁴C can be used as an additional biomarker to untangle riverine food webs² where aquatic and terrestrial δ ¹³C overlap, but may also be valuable for examining the age and origin of C in the lake. Primary production in lakes is based on dissolved inorganic C (DIC). DIC in alkaline lakes is partially derived from weathering of carbonaceous bedrock, a proportion of which is¹⁴C-free. The low ¹⁴C activity yields an artificial age offset leading samples to appear hundreds to thousands of years older than their actual age. As such, ¹⁴C can be used to identify the proportion of autochthonous C in the food-web. With terrestrial C inputs likely to increase, the origin and utilisation of ‘fossil’ or ‘recent’ allochthonous C in the food-web can also be determined. Stable isotopes and 14C were measured for biota, particulate organic matter (POM), DIC and dissolved organic carbon (DOC) from Lough Erne, Northern Ireland, a humic alkaline lake. Temporal and spatial variation was evident in DIC, DOC and POM C isotopes with implications for the fluctuation in terrestrial export processes. Ramped pyrolysis of lake surface sediment indicates the burial of two C components. ¹⁴C activity (507 ± 30 BP) of sediment combusted at 400˚C was consistent with algal values and younger than bulk sediment values (1097 ± 30 BP). The sample was subsequently combusted at 850˚C, yielding 14C values (1471 ± 30 BP) older than the bulk sediment age, suggesting that fossil terrestrial carbon is also buried in the sediment. Stable isotopes in the food web indicate that terrestrial organic C is also utilised by lake organisms. High winter δ ¹⁵N values in calanoid zooplankton (δ ¹⁵N = 24%¸) relative to phytoplankton and POM (δ ¹⁵N = 6h and 12h respectively) may reflect several microbial trophic levels between terrestrial C and calanoids. Furthermore winter calanoid 14C ages are consistent with DOC from an inflowing river (75 ± 24 BP), not phytoplankton (367 ± 70 BP). Summer calanoid δ ^13C, δ ¹⁵N and ¹⁴C (345 ± 80 BP) indicate greater reliance on phytoplankton.

¹ Monteith, D.T et al., (2007) Dissolved organic carbon trends resulting from changes in atmospheric deposition chemistry. Nature, 450:537-535

² Caraco, N., et al.,(2010) Millennial-aged organic carbon subsidies to a modern river food web. Ecology,91: 2385-2393.

Storage of sediment-associated nutrients and contaminants in river channel and floodplain systems

Samples of fine-grained channel bed sediment and overbank floodplain deposits were collected along the main channels of the Rivers Aire (and its main tributary, the River Calder) and Swale, in Yorkshire, UK, in order to investigate downstream changes in the storage and deposition of heavy metals (Cr, Cu, Pb, Zn), total P and the sum of selected PCB congeners, and to estimate the total storage of these contaminants within the main channels and floodplains of these river systems. Downstream trends in the contaminant content of the <63 μm fraction of channel bed and floodplain sediment in the study rivers are controlled mainly by the location of the main sources of the contaminants, which varies between rivers. In the Rivers Aire and Calder, the contaminant content of the <63 μm fraction of channel bed and floodplain sediment generally increases in a downstream direction, reflecting the location of the main urban and industrialized areas in the middle and lower parts of the basin. In the River Swale, the concentrations of most of the contaminants examined are approximately constant along the length of the river, due to the relatively unpolluted nature of this river. However, the Pb and Zn content of fine channel bed sediment decreases downstream, due to the location of historic metal mines in the headwaters of this river, and the effect of downstream dilution with uncontaminated sediment. The magnitude and spatial variation of contaminant storage and deposition on channel beds and floodplains are also controlled by the amount of <63 μm sediment stored on the channel bed and deposited on the floodplain during overbank events. Consequently, contaminant deposition and storage are strongly influenced by the surface area of the floodplain and channel bed. Contaminant storage on the channel beds of the study rivers is, therefore, generally greatest in the middle and lower reaches of the rivers, since channel width increases downstream. Comparisons of the estimates of total storage of specific contaminants on the channel beds of the main channel systems of the study rivers with the annual contaminant flux at the catchment outlets indicate that channel storage represents <3% of the outlet flux and is, therefore, of limited importance in regulating that flux. Similar comparisons between the annual deposition flux of specific contaminants to the floodplains of the study rivers and the annual contaminant flux at the catchment outlet, emphasise the potential importance of floodplain deposition as a conveyance loss. In the case of the River Aire the floodplain deposition flux is equivalent to between ca. 2% (PCBs) and 36% (Pb) of the outlet flux. With the exception of PCBs, for which the value is ≅0, the equivalent values for the River Swale range between 18% (P) and 95% (Pb). The study emphasises that knowledge of the fine-grained sediment delivery system operating in a river basin is an essential prerequisite for understanding the transport and storage of sediment-associated contaminants in river systems and that conveyance losses associated with floodplain deposition exert an important control on downstream contaminant fluxes and the fate of such contaminants. © 2003 Elsevier Science Ltd. All rights reserved.