Multivariate analysis of regional-scale geochemical data for environmental monitoring

A compositional multivariate approach is used to analyse regional scale soil geochemical data obtained as part of the Tellus Project generated by the Geological Survey Northern Ireland (GSNI). The multi-element total concentration data presented comprise XRF analyses of 6862 rural soil samples collected at 20cm depths on a non-aligned grid at one site per 2 km2. Censored data were imputed using published detection limits. Using these imputed values for 46 elements (including LOI), each soil sample site was assigned to the regional geology map provided by GSNI initially using the dominant lithology for the map polygon. Northern Ireland includes a diversity of geology representing a stratigraphic record from the Mesoproterozoic, up to and including the Palaeogene. However, the advance of ice sheets and their meltwaters over the last 100,000 years has left at least 80% of the bedrock covered by superficial deposits, including glacial till and post-glacial alluvium and peat. The question is to what extent the soil geochemistry reflects the underlying geology or superficial deposits. To address this, the geochemical data were transformed using centered log ratios (clr) to observe the requirements of compositional data analysis and avoid closure issues. Following this, compositional multivariate techniques including compositional Principal Component Analysis (PCA) and minimum/maximum autocorrelation factor (MAF) analysis method were used to determine the influence of underlying geology on the soil geochemistry signature. PCA showed that 72% of the variation was determined by the first four principal components (PC’s) implying “significant” structure in the data. Analysis of variance showed that only 10 PC’s were necessary to classify the soil geochemical data. To consider an improvement over PCA that uses the spatial relationships of the data, a classification based on MAF analysis was undertaken using the first 6 dominant factors. Understanding the relationship between soil geochemistry and superficial deposits is important for environmental monitoring of fragile ecosystems such as peat. To explore whether peat cover could be predicted from the classification, the lithology designation was adapted to include the presence of peat, based on GSNI superficial deposit polygons and linear discriminant analysis (LDA) undertaken. Prediction accuracy for LDA classification improved from 60.98% based on PCA using 10 principal components to 64.73% using MAF based on the 6 most dominant factors. The misclassification of peat may reflect degradation of peat covered areas since the creation of superficial deposit classification. Further work will examine the influence of underlying lithologies on elemental concentrations in peat composition and the effect of this in classification analysis.

Circadian variation with a diurnal bimodal profile on scent-marking behavior in captive common marmosets (callitrhix jacchus)

Scent-marking behavior is associated with different behavioral contexts in callitrichids, including signalizing a territory, location of feeding resources, and social rank. In marmosets and tamarins it is also associated with intersexual communication. Though it appears very important for the daily routine of the individuals, very few researchers have investigated distribution through the 24-h cycle. In a preliminary report, we described a preferential incidence of this behavior 2 h before nocturnal rest in families of common marmosets. We expand the data using 8 family groups (28 subjects), 8 fathers, 6 mothers, 8 nonreproductive adults (4 sons and 4 daughters), and 6 juvenile (3 sons and 3 daughters) offspring that we kept in outdoor cages under natural environmental conditions. We recorded the frequency of anogenital scent marking for each group during the light phase, twice a wk, for 4 consecutive wks, from March 1998 to September 1999. Cosinor test detected 24- and 8-h variations in 89.3% and 85.7% of the subjects, respectively, regardless of sex or reproductive status. The 8-h component is a consequence of the 2 peaks for the behavior, at the beginning and end of the light phase. Daily distribution of scent marking is similar to that others described previously for motor activity in marmosets. The coincident rhythmical patterns for both behaviors seem to be associated with feeding behavior, as described for callitrichids in free-ranging conditions, involving an increase in foraging activities early in the morning and shortly before nocturnal rest

Circadian variation with a diurnal bimodal profile on scent-marking behavior in captive common marmosets (callitrhix jacchus)

Scent-marking behavior is associated with different behavioral contexts in callitrichids, including signalizing a territory, location of feeding resources, and social rank. In marmosets and tamarins it is also associated with intersexual communication. Though it appears very important for the daily routine of the individuals, very few researchers have investigated distribution through the 24-h cycle. In a preliminary report, we described a preferential incidence of this behavior 2 h before nocturnal rest in families of common marmosets. We expand the data using 8 family groups (28 subjects), 8 fathers, 6 mothers, 8 nonreproductive adults (4 sons and 4 daughters), and 6 juvenile (3 sons and 3 daughters) offspring that we kept in outdoor cages under natural environmental conditions. We recorded the frequency of anogenital scent marking for each group during the light phase, twice a wk, for 4 consecutive wks, from March 1998 to September 1999. Cosinor test detected 24- and 8-h variations in 89.3% and 85.7% of the subjects, respectively, regardless of sex or reproductive status. The 8-h component is a consequence of the 2 peaks for the behavior, at the beginning and end of the light phase. Daily distribution of scent marking is similar to that others described previously for motor activity in marmosets. The coincident rhythmical patterns for both behaviors seem to be associated with feeding behavior, as described for callitrichids in free-ranging conditions, involving an increase in foraging activities early in the morning and shortly before nocturnal rest

Variation in composition of macro-benthic invertebrates as an indication of water quality status in three bays in Lake Victoria

Knowledge of how biota can be used to monitor ecosystem health and assess impacts by human alterations such as land use and management measures taken at different spatial scales is critical for improving the ecological quality of aquatic ecosystems. This knowledge in Uganda is very limited or unavailable yet it is needed to better understand the relationship between environmental factors at different spatial scales, assemblage structure and taxon richness of aquatic ecosystems. In this study, benthic invertebrate community patterns were sampled between June 2001 and April 2002 and analysed in relation to water quality and catchment land use patterns from three shallow near-shore bays characterized by three major land uses patterns: urban (Murchison Bay); semi-urban (Fielding Bay); rural (Hannington Bay). Variations in density and guild composition of benthic macro-invertebrates communities were evaluated using GIS techniques along an urban-rural gradient of land use and differences in community composition were related to dissolved oxygen and conductivity variation. Based on numerical abundance and tolerance values, Hilsenhoff's Biotic Index ofthe invertebrates was determined in order to evaluate the relative importance of water quality in the three bays. Murchison Bay supported a relatively taxa-poor invertebrate assemblage mainly comprising stenotopic and eurytopic populations of pollution-tolerant groups such as worms and Chironomus sp. with an overall depression in species diversity. On the contrary, the communities in Fielding and Hannington bays were quite similar and supported distinct and diverse assemblages including pollution-intolerant forms such as Ephemeroptera (mayflies), Odonata (dragonflies). The Hilsenhoff Biotic Index in Murchison Bay was 6.53. (indicating poor water quality) compared to 6.34 for Fielding Bay and 5.78 for Hannington Bay (both indicating fair water quality). The characterization of maximum taxa richness balanced among taxa groups with good representation of intolerant individuals in Hannington Bay relative to Fielding and Murchison bays concludes that the bay is the cleanest in terms of water quality. Contrary, the dominance of few taxa with many tolerant iqdividuals present in Murchison Bay indicates that the bay is degraded in terms of water quality. These result are ofimportance when planning conservation and management measures, implementing large-scale biomonitoring programs, and predicting how human alterations (e.g nutrient loading) affect water ecosystems. Therefore, analysis of water quality in relation to macro-invertebrate community composition patterns as bio-indicators can lead to further understanding of their responses to environmental manipulations and perturbations.

Environmental factors affecting maerl bed structure in Brittany (France)

This study used a large spatial scale approach in order to better quantify the relationships between maerl bed structure and a selection of potentially forcing physical factors. Data on maerl bed structure and morpho-sedimentary characteristics were obtained from recent oceanographic surveys using underwater video recording and grab sampling. Considering the difficulties in carrying out real-time monitoring of highly variable hydrodynamic and physicochemical factors, these were generated by three-dimensional numerical models with high spatial and temporal resolution. The BIOENV procedure indicated that variation in the percentage cover of thalli can best be explained (correlation = 0.76) by a combination of annual mean salinity, annual mean nitrate concentration and annual mean current velocity, while the variation in the proportion of living thalli can best be explained (correlation = 0.47) by a combination of depth and mud content. Linear relationships showed that the percentage cover of maerl thalli was positively correlated with nitrate concentration (R2 = 0.78, P < 0.01) and negatively correlated with salinity (R2 = 0.81, P < 0.01), suggesting a strong effect of estuarine discharge on maerl bed structure, and also negatively correlated with current velocity (R2 = 0.81, P < 0.01). When maerl beds were deeper than 10 m, the proportion of living thalli was always below 30% but when they were shallower than 10 m, it varied between 4 and 100%, and was negatively correlated with mud content (R2 = 0.53, P < 0.01). On the other hand, when mud content was below 10%, the proportion of living thalli showed a negative correlation with depth (R2 = 0.84, P < 0.01). This large spatial scale explanation of maerl bed heterogeneity provides a realistic physical characterization of these ecologically interesting benthic habitats and usable findings for their conservation and management.

Temporal and spatial variation in temperature experienced by macrofauna at main endeavour Hydrothermal vent field

A significant focus of hydrothermal vent ecological studies has been to understand how species cope with various stressors through physiological tolerance and biochemical resistance. Yet, the environmental conditions experienced by vent species have not been well characterized. This objective requires continuous observations over time intervals that can capture environmental variability at scales that are relevant to animals. We used autonomous temperature logger arrays (four roughly parallel linear arrays of 12 loggers spaced every 10–12 cm) to study spatial and temporal variations in the thermal regime experienced by hydrothermal vent macrofauna at a diffuse flow vent. Hourly temperatures were recorded over eight months from 2010 to 2011 at Grotto vent in the Main Endeavour vent field on the Juan de Fuca Ridge, a focus area of the Ocean Networks Canada cabled observatory. The conspicuous animal assemblages in video footage contained Ridgeia piscesae tubeworms, gastropods (primarily Lepetodrilus fucensis), and polychaetes (polynoid scaleworms and the palm worm Paralvinella palmiformis). Two dimensional spatial gradients in temperature were generally stable over the deployment period. The average temperature recorded by all arrays, and in some individual loggers, revealed distinctive fluctuations in temperature that often corresponded with the tidal cycle. We postulate that this may be related to changes in bottom currents or fluctuations in vent discharge. A marked transient temperature increase lasting over a period of days was observed in April 2011. While the distributions and behavior of Juan de Fuca Ridge vent invertebrates may be partially constrained by environmental temperature and temperature tolerance, except for the one transient high-temperature event, observed fluid temperatures were generally similar to the thermal preferences for some species, and typically well below lethal temperatures for all species. Average temperatures of the four arrays ranged from 4.1 to 11.0 °C during the deployment, indicating that on an hourly timescale the temperature conditions in this tubeworm community were fairly moderate and stable. The generality of these findings and behavioural responses of vent organisms to predictable rhythmicity and non-periodic temperature shifts are areas for further investigation

An investigation on benthic macro-invertebrates species diversity as bio-indicators of environmental health in Bahrekan Bay

The purpose of this study was to evaluate benthic macro-invertebrates species diversity as bio-indicators of environmental health in Bahrekan bay (in the Northwest of Persian gulf). Seasonal sediments sampling along 5 transects, 15 stations at 4 replicates (3 replicates for macrobenthos and 1 replicate for sediment analysis) was done from November 2008 to August 2009 by 0.025 m2 Van Veen grab sampler. Physical and chemical parameters of water, grain size analysis, %TOM and Ni and Va concentrations of sediments were assessed through four seasons. Macrobenthic communities after staining and sorting, using stereomicroscope have been identified. Their density in every station and every season calculated. For using of AMBI index, identified macrobenthos according to their sensitivity to stressors and pollutants, categorized into 5 ecological groups and for using of Bentix index categorized into 3 ecological groups. The diversity indices and indicators that showing ecological status were calculated. Also, the differences between physiochemical parameters of sea water, sediments TOM% and grain size, diversity indices in stations and seasons were recorded (P=0.05). The correlation coefficient determined for all parameters. According to the results of grain size analysis, bottom grain size categorized as clay. Highest percent of TOM was belong to autumn (36.39±.075) and lowest was belong to summer (19.01±0.51). Also there was positive correlation (p=0.01) between %TOM and %Clay that showing sediments with lowest size containing highest amounts of organic matters. Ni concentrations in sediments (87.80±21.25)mg/kg showed the amounts over than standards levels but Va concentrations in sediments (53.54±17.60)mg/kg showed the amounts lower than standards level. The highest density of macrobenthos was recorded for summer (8254±485) N/m2 and the lowest density was recorded for spring (3775±172)N/m2. The highest annual density was belong to mollusca (81%) and then polycheates (13%), Others (4%) and crustaceae (2%). The highest diversity was recorded for winter (Simpson index: 0.13±0.01, H':3.47±0.06) and the lowest diversity recorded for autumn (Simpson index: 0.16±0.01, H':3.17±0.06). in all stations, the highest amount of Shanon index was belong to T2S3 station in summer (4.11± 0.32) and the lowest amount was belong to T1S1 station in autumn (2.42± 0.41). The annual mean of Simpson diversity index: (0.15 ±0.04) and Shanon diversity index (3.36±0.03), illustrated that macrobenthos in Bahrekan bay have a good variation. The results of Brilluin and N1 (Number of equally common species) indices confirm the results of Simpson index. For study on the regions that diversity has a little difference between stations, with use of Ni index, the degree of differences could be better ono recognizable. According to the results of AMBI index in all seasons (autumn: 0.46±0.03; summer: 0.22±0.01; annual mean:0.31±0.01) and standards (0.0

Investigation on water and sediment quality variation from coast to off shore and its relation to bent hoses biodiversity: case study Bandar tahery

This paper deals with the qualification of water and sediment particularly those of benthos, as well as their interaction results, through careful laboratory researches within 24 experimental sites around the bank of Bandar Taheri, Persian Gulf water of Iranian borders, under 13th phase south Pars project. Samples were carried to the laboratories and careful experimental tests such as physical chemical, heavy metal, nitrate ammoniac, toc and other biologic tests including various type of benthos count were performed. Data gained through Shanon and Dankan statistical analyses were also studied to determine the water, sediment and pollution rate. Resulted information would classify the area as less polluted area which is rather away from critical environmental zone, another word the area could be liable to change to an undesired one while the density rates of principal metals follows the Fe>Mn>Zn>Cu>Pb>Cd>N pattern.

Dynamique de la variation génétique et épigénétique chez un vertébré kleptogène

La variation phénotypique est essentielle à la persistance des organismes dans le temps ainsi qu’à la colonisation de nouveaux habitats. Les principales sources de variation phénotypique sont la génétique et l'épigénétique. L'épigénétique a été proposé comme un atout important pour les organismes asexués pour compenser le manque de diversité génétique. L'objectif de cette étude est d'évaluer si l’absence de variation génétique est compensée par l'épigénétique en comparant les profils de méthylation d’individus gynogènes et kleptogènes des hybrides de salamandre à points bleus. Les individus échantillonnés s’organisent en cinq groupes génétiquement différenciés, provenant du même haplome paternel A. jeffersonianum. Deux des cinq groupes sont exclusivement gynogènes, pour des raisons écologiques ou génomiques. Les trois autres groupes sont formés d’individus parfois kleptogènes, car ils présentent une variation génétique plus élevée au sein d’un site qu’entre les sites, en plus de porter des allèles très divergents par rapport à la distribution globale des allèles hybrides, trouvés en haute fréquence dans les populations sympatriques de A. laterale. Les patrons épigénétiques sont variables et distincts entre les cinq groupes génétiques. Les groupes gynogènes sont les seuls à présenter un effet environnemental significatif sur leurs patrons épigénétiques, suggérant que ces individus clonaux doivent être en mesure de maximiser leur potentiel de variation épigénétique pour faire face à des variations environnementales.

Method of numerical correction of errors occasioned by delay of records during the monitoring of environmental variables of interest for animal production

As condições de ambiente térmico e aéreo, no interior de instalações para animais, alteram-se durante o dia, devido à influência do ambiente externo. Para que análises estatísticas e geoestatísticas sejam representativas, uma grande quantidade de pontos distribuídos espacialmente na área da instalação deve ser monitorada. Este trabalho propõe que a variação no tempo das variáveis ambientais de interesse para a produção animal, monitoradas no interior de instalações para animais, pode ser modelada com precisão a partir de registros discretos no tempo. O objetivo deste trabalho foi desenvolver um método numérico para corrigir as variações temporais dessas variáveis ambientais, transformando os dados para que tais observações independam do tempo gasto durante a aferição. O método proposto aproximou os valores registrados com retardos de tempo aos esperados no exato momento de interesse, caso os dados fossem medidos simultaneamente neste momento em todos os pontos distribuídos espacialmente. O modelo de correção numérica para variáveis ambientais foi validado para o parâmetro ambiental temperatura do ar, sendo que os valores corrigidos pelo método não diferiram pelo teste Tukey, a 5% de probabilidade dos valores reais registrados por meio de dataloggers.

Genetic diversity and connectivity remain high in Holothuria polii (Delle Chiaje 1823) across a coastal lagoon-open sea environmental gradient

Coastal lagoons represent habitats with widely heterogeneous environmental conditions, particularly as regards salinity and temperature,which fluctuate in both space and time. These characteristics suggest that physical and ecological factors could contribute to the genetic divergence among populations occurring in coastal lagoon and opencoast environments. This study investigates the genetic structure of Holothuria polii at a micro-geographic scale across theMar Menor coastal lagoon and nearbymarine areas, estimating the mitochondrial DNA variation in two gene fragments, cytochrome oxidase I (COI) and 16S rRNA (16S). Dataset of mitochondrial sequences was also used to test the influence of environmental differences between coastal lagoon andmarine waters on population genetic structure. All sampled locations exhibited high levels of haplotype diversity and low values of nucleotide diversity. Both genes showed contrasting signals of genetic differentiation (non-significant differences using COI and slight differences using 16S, which could due to different mutation rates or to differential number of exclusive haplotypes. We detected an excess of recent mutations and exclusive haplotypes, which can be generated as a result of population growth. However, selective processes can be also acting on the gene markers used; highly significant generalized additive models have been obtained considering genetic data from16S gene and independent variables such as temperature and salinity.

GENETIC VARIATION, LOCAL ADAPTATION AND POPULATION STRUCTURE IN NORTH AMERICAN RED OAK SPECIES, QUERCUS RUBRA L. AND Q. ELLIPSOIDALIS E. J. HILL

Forest trees, like oaks, rely on high levels of genetic variation to adapt to varying environmental conditions. Thus, genetic variation and its distribution are important for the long-term survival and adaptability of oak populations. Climate change is projected to lead to increased drought and fire events as well as a northward migration of tree species, including oaks. Additionally, decline in oak regeneration has become increasingly concerning since it may lead to decreased gene flow and increased inbreeding levels. This will in turn lead to lowered levels of genetic diversity, negatively affecting the growth and survival of populations. At the same time, populations at the species’ distribution edge, like those in this study, could possess important stores of genetic diversity and adaptive potential, while also being vulnerable to climatic or anthropogenic changes. A survey of the level and distribution of genetic variation and identification of potentially adaptive genes is needed since adaptive genetic variation is essential for their long-term survival. Oaks possess a remarkable characteristic in that they maintain their species identity and specific environmental adaptations despite their propensity to hybridize. Thus, in the face of interspecific gene flow, some areas of the genome remain differentiated due to selection. This characteristic allows the study of local environmental adaptation through genetic variation analyses. Furthermore, using genic markers with known putative functions makes it possible to link those differentiated markers to potential adaptive traits (e.g., flowering time, drought stress tolerance). Demographic processes like gene flow and genetic drift also play an important role in how genes (including adaptive genes) are maintained or spread. These processes are influenced by disturbances, both natural and anthropogenic. An examination of how genetic variation is geographically distributed can display how these genetic processes and geographical disturbances influence genetic variation patterns. For example, the spatial clustering of closely related trees could promote inbreeding with associated negative effects (inbreeding depression), if gene flow is limited. In turn this can have negative consequences for a species’ ability to adapt to changing environmental conditions. In contrast, interspecific hybridization may also allow the transfer of genes between species that increase their adaptive potential in a changing environment. I have studied the ecologically divergent, interfertile red oaks, Quercus rubra and Q. ellipsoidalis, to identify genes with potential roles in adaptation to abiotic stress through traits such as drought tolerance and flowering time, and to assess the level and distribution of genetic variation. I found evidence for moderate gene flow between the two species and low interspecific genetic differences at most genetic markers (Lind and Gailing 2013). However, the screening of genic markers with potential roles in phenology and drought tolerance led to the identification of a CONSTANS-like (COL) gene, a candidate gene for flowering time and growth. This marker, located in the coding region of the gene, was highly differentiated between the two species in multiple geographical areas, despite interspecific gene flow, and may play a role in reproductive isolation and adaptive divergence between the two species (Lind-Riehl et al. 2014). Since climate change could result in a northward migration of trees species like oaks, this gene could be important in maintaining species identity despite increased contact zones between species (e.g., increased gene flow). Finally I examined differences in spatial genetic structure (SGS) and genetic variation between species and populations subjected to different management strategies and natural disturbances. Diverse management activities combined with various natural disturbances as well as species specific life history traits influenced SGS patterns and inbreeding levels (Lind-Riehl and Gailing submitted).

ANALYSIS OF POPULUS TREMULOIDES CLONAL VARIATION AND DELINEATION IN THE OTTAWA NATIONAL FOREST

The technique of delineating Populus tremuloides (Michx.) clonal colonies based on morphology and phenology has been utilized in many studies and forestry applications since the 1950s. Recently, the availability and robustness of molecular markers has challenged the validity of such approaches for accurate clonal identification. However, genetically sampling an entire stand is largely impractical or impossible. For that reason, it is often necessary to delineate putative genet boundaries for a more selective approach when genetically analyzing a clonal population. Here I re-evaluated the usefulness of phenotypic delineation by: (1) genetically identifying clonal colonies using nuclear microsatellite markers, (2) assessing phenotypic inter- and intraclonal agreement, and (3) determining the accuracy of visible characters to correctly assign ramets to their respective genets. The long-term soil productivity study plot 28 was chosen for analysis and is located in the Ottawa National Forest, MI (46° 37'60.0" N, 89° 12'42.7" W). In total, 32 genets were identified from 181 stems using seven microsatellite markers. The average genet size was 5.5 ramets and six of the largest were selected for phenotypic analyses. Phenotypic analyses included budbreak timing, DBH, bark thickness, bark color or brightness, leaf senescence, leaf serrations, and leaf length ratio. All phenotypic characters, except for DBH, were useful for the analysis of inter- and intraclonal variation and phenotypic delineation. Generally, phenotypic expression was related to genotype with multiple response permutation procedure (MRPP) intraclonal distance values ranging from 0.148 and 0.427 and an observed MRPP delta value=0.221 when the expected delta=0.5. The phenotypic traits, though, overlapped significantly among some clones. When stems were assigned into phenotypic groups, six phenotypic groups were identified with each group containing a dominant genotype or clonal colony. All phenotypic groups contained stems from at least two clonal colonies and no clonal colony was entirely contained within one phenotypic group. These results demonstrate that phenotype varies with genotype and stand clonality can be determined using phenotypic characters, but phenotypic delineation is less precise. I therefore recommend that some genetic identification follow any phenotypic delineation. The amount of genetic identification required for clonal confirmation is likely to vary based on stand and environmental conditions. Further analysis, however, is needed to test these findings in other forest stands and populations.

Life history & environmental effects on telomere dynamics in Atlantic salmon

While much of the study of molecular biology inevitably focuses on the parts of the genome that contain active genes, there are also non-coding regions that nonetheless play an essential role in maintaining genome integrity. One such region are telomeres, which cap the ends of all eukaryotic chromosomes and play an important role in chromosome protection. Telomere loss occurs at each cell division as a result of the ‘end replication problem’ and a relatively short telomere length is indicative of poor biological state. Thus far, the majority of studies on the dynamics and role of telomeres have been biased towards certain taxa. Research to date has mostly focussed on humans, other mammals and birds. There has been far less research on the telomere dynamics of ectotherms. It is important that we do so, especially since ectothermic vertebrates do not seem to down-regulate telomerase expression in the same way as endotherms, suggesting that their telomere dynamics may be less predictable in the later life stages. The main objective of this thesis was therefore to investigate how life history and environmental effects may influence telomere dynamics in Atlantic salmon Salmo salar. I carried out carefully designed experiments, both in the laboratory and in the wild, using a longitudinal approach where possible, in order to address a number of specific questions that are connected to this central theme. In chapter 2, I demonstrate that there can be significant links between parental life history and offspring telomere dynamics. Maternal life history traits, in particular egg size, were most strongly related to offspring telomere length at the embryonic stages. Paternal life history traits, such as early life growth rate, had a greater association with offspring telomere dynamics in the later stages of development. In chapter 3, using a wild Atlantic salmon population, I found that most individuals experienced a reduction in telomere length during the migratory phase of their life cycle; however the relative rate of telomere loss was dependent on sex, with males experiencing a relatively greater loss. Unexpectedly, I also found that juvenile salmon that had the shortest telomeres at the time of outward migration, had the greatest probability of surviving through to the return migration. In chapter 4, again using a wild system involving experimental manipulations of juvenile Atlantic salmon in Scottish streams, I found that telomere length in juvenile fish was influenced by parental traits and by direct environmental effects. Faster-growing fish had shorter telomeres and there was a greater cost (in terms of reduced telomere length) if the growth occurred in a harsher environment. I also found a positive association between offspring telomere length and the growth history of their fathers (but not mothers), represented by the number of years that fathers had spent at sea. Chapter 5 explored the hypotheses that oxidative DNA damage, catalase (CAT) antioxidant activity and cell proliferation rate are underlying mechanisms linking incubation temperature and telomere dynamics in salmon embryos. No evidence was found for any such effects, but telomere lengths in salmon embryos were found to be significantly affected by the temperature of the water in which they were living. There is also evidence that telomere length significantly increases during embryonic development. In summary, this thesis has shown that a complex mix of environmental and parental effects appear to influence telomere dynamics in Atlantic salmon, with parental effects especially evident during early life stages. It also demonstrated that telomeres lengthen through the embryo stages of development before reducing once the fry begin feeding, indicating that the patterns of telomere loss commonly found in endotherms may differ in ectotherms. Reasons for this variation in telomere dynamics are presented in the final Discussion chapter of the thesis.

Obtaining environmental favourability functions from logistic regression

Logistic regression is a statistical tool widely used for predicting species’ potential distributions starting from presence/absence data and a set of independent variables. However, logistic regression equations compute probability values based not only on the values of the predictor variables but also on the relative proportion of presences and absences in the dataset, which does not adequately describe the environmental favourability for or against species presence. A few strategies have been used to circumvent this, but they usually imply an alteration of the original data or the discarding of potentially valuable information. We propose a way to obtain from logistic regression an environmental favourability function whose results are not affected by an uneven proportion of presences and absences. We tested the method on the distribution of virtual species in an imaginary territory. The favourability models yielded similar values regardless of the variation in the presence/absence ratio. We also illustrate with the example of the Pyrenean desman’s (Galemys pyrenaicus) distribution in Spain. The favourability model yielded more realistic potential distribution maps than the logistic regression model. Favourability values can be regarded as the degree of membership of the fuzzy set of sites whose environmental conditions are favourable to the species, which enables applying the rules of fuzzy logic to distribution modelling. They also allow for direct comparisons between models for species with different presence/absence ratios in the study area. This makes themmore useful to estimate the conservation value of areas, to design ecological corridors, or to select appropriate areas for species reintroductions.