Influence of pollination and seed dispersal on the genetic population structure of two Commiphora species in Madagascar and South Africa

Pollination and seed dispersal are important ecological processes for the regeneration of plant populations and both vectors for gene exchange between plant populations. For my thesis, I studied the pollination ecology of the South African tree Commiphora harveyi (Burseraceae) and compared it with C. guillauminii from Madagascar. Both species have low visitation rates and a low number of pollinating insect species, resulting in a low fruit set. While their pollination ecology is very similar, they differ in their seed dispersal with a low seed dispersal rate in the Malagasy and a high seed dispersal rate in the South African species. This should be reflected in a stronger genetic differentiation among populations in the Malagasy than in the South African species. My results, based on AFLP markers, contradict these expectations, the overall differentiation was lower in the Malagasy (FST = 0.05) than in the South African species (FST = 0.16). However, at a smaller spatial scale (below 3 km), the Malagasy species was genetically more strongly differentiated than the South African species, which was reflected by the high inter-population variance within the sample site (C. guillauminii: 72.2 - 85.5 %; C. harveyi: 8.4 - 14.5 %). This strong differentiation could arise from limited gene flow, which was confirmed by spatial autocorrelation analyses. The shape of the autocorrelogram suggested that gene exchange between individuals occurred only up to 3 km in the Malagasy species, whereas up to 30 km in the South African species. These results on the genetic structure correspond to the expectations based on seed dispersal data. Thus, seed dispersal seems to be a key factor for the genetic structure in plant populations on a local scale.

Spatial and temporal variability and ecological processes in the epibenthic assemblages of the northern Adriatic Sea

Several coralligenous reefs occur in the soft bottoms of the northern Adriatic continental shelf. Mediterranean coralligenous habitats are characterised by high species diversity and are intrinsically valuable for their biological diversity and for the ecological processes they support. The conservation and management of these habitats require quantifying spatial and temporal variability of their benthic assemblages. This PhD thesis aims to give a relevant contribution to the knowledge of the structure and dynamics of the epibenthic assemblages on the coralligenous subtidal reefs occurring in the northern Adriatic Sea. The epibenthic assemblages showed a spatial variation larger compared to temporal changes, with a temporal persistence of reef-forming organisms. Assemblages spatial heterogeneity has been related to morphological features and geographical location of the reefs, together with variation in the hydrological conditions. Manipulative experiments help to understand the ecological processes structuring the benthic assemblages and maintaining their diversity. In this regards a short and long term experiment on colonization patterns of artificial substrata over a 3-year period has been performed in three reefs, corresponding to the three main types of assemblages detected in the previous study. The first colonisers, largely depending by the different larval supply, played a key role in determining the heterogeneity of the assemblages in the early stage of colonisation. Lateral invasion, from the surrounding assemblages, was the driver in structuring the mature assemblages. These complex colonisation dynamics explained the high heterogeneity of the assemblages dwelling on the northern Adriatic biogenic reefs. The buildup of these coralligenous reefs mainly depends by the bioconstruction-erosion processes that has been analysed through a field experiment. Bioconstruction, largely due to serpulid polychaetes, prevailed on erosion processes and occurred at similar rates in all sites. Similarly, the total energy contents in the benthic communities do not differ among sites, despite being provided by different species. Therefore, we can hypothesise that both bioconstruction processes and energetic storage may be limited by the availability of resources. Finally the major contribution of the zoobenthos compared to the phytobenthos to the total energetic content of assemblages suggests that the energy flow in these benthic habitats is primarily supported by planktonic food web trough the filter feeding invertebrates.

Novel tools for conservation genetics in marine fish: population structure and evolution of the European hake (Merluccius merluccius) inferred by SNP variation and applications to traceability

The research presented in my PhD thesis is part of a wider European project, FishPopTrace, focused on traceability of fish populations and products. My work was aimed at developing and analyzing novel genetic tools for a widely distributed marine fish species, the European hake (Merluccius merluccius), in order to investigate population genetic structure and explore potential applications to traceability scenarios. A total of 395 SNPs (Single Nucleotide Polymorphisms) were discovered from a massive collection of Expressed Sequence Tags, obtained by high-throughput sequencing, and validated on 19 geographic samples from Atlantic and Mediterranean. Genome-scan approaches were applied to identify polymorphisms on genes potentially under divergent selection (outlier SNPs), showing higher genetic differentiation among populations respect to the average observed across loci. Comparative analysis on population structure were carried out on putative neutral and outlier loci at wide (Atlantic and Mediterranean samples) and regional (samples within each basin) spatial scales, to disentangle the effects of demographic and adaptive evolutionary forces on European hake populations genetic structure. Results demonstrated the potential of outlier loci to unveil fine scale genetic structure, possibly identifying locally adapted populations, despite the weak signal showed from putative neutral SNPs. The application of outlier SNPs within the framework of fishery resources management was also explored. A minimum panel of SNP markers showing maximum discriminatory power was selected and applied to a traceability scenario aiming at identifying the basin (and hence the stock) of origin, Atlantic or Mediterranean, of individual fish. This case study illustrates how molecular analytical technologies have operational potential in real-world contexts, and more specifically, potential to support fisheries control and enforcement and fish and fish product traceability.

Effect of mixing, particle interaction, and spatial dimension on the glass transition

In this thesis, the influence of composition changes on the glass transition behavior of binary liquids in two and three spatial dimensions (2D/3D) is studied in the framework of mode-coupling theory (MCT).The well-established MCT equations are generalized to isotropic and homogeneous multicomponent liquids in arbitrary spatial dimensions. Furthermore, a new method is introduced which allows a fast and precise determination of special properties of glass transition lines. The new equations are then applied to the following model systems: binary mixtures of hard disks/spheres in 2D/3D, binary mixtures of dipolar point particles in 2D, and binary mixtures of dipolar hard disks in 2D. Some general features of the glass transition lines are also discussed. The direct comparison of the binary hard disk/sphere models in 2D/3D shows similar qualitative behavior. Particularly, for binary mixtures of hard disks in 2D the same four so-called mixing effects are identified as have been found before by Götze and Voigtmann for binary hard spheres in 3D [Phys. Rev. E 67, 021502 (2003)]. For instance, depending on the size disparity, adding a second component to a one-component liquid may lead to a stabilization of either the liquid or the glassy state. The MCT results for the 2D system are on a qualitative level in agreement with available computer simulation data. Furthermore, the glass transition diagram found for binary hard disks in 2D strongly resembles the corresponding random close packing diagram. Concerning dipolar systems, it is demonstrated that the experimental system of König et al. [Eur. Phys. J. E 18, 287 (2005)] is well described by binary point dipoles in 2D through a comparison between the experimental partial structure factors and those from computer simulations. For such mixtures of point particles it is demonstrated that MCT predicts always a plasticization effect, i.e. a stabilization of the liquid state due to mixing, in contrast to binary hard disks in 2D or binary hard spheres in 3D. It is demonstrated that the predicted plasticization effect is in qualitative agreement with experimental results. Finally, a glass transition diagram for binary mixtures of dipolar hard disks in 2D is calculated. These results demonstrate that at higher packing fractions there is a competition between the mixing effects occurring for binary hard disks in 2D and those for binary point dipoles in 2D.

Structure and dynamics of ionic liquid,electrode interfaces

In the early 20th century, Gouy, Chapman, and Stern developed a theory to describe the capacitance and the spatial ion distribution of diluted electrolytes near an electrode. After a century of research, considerable progress has been made in the understanding of the electrolyte/electrode interface. However, its molecular-scale structure and its variation with an applied potential is still under debate. In particular for room-temperature ionic liquids, a new class of solventless electrolytes, the classical theories for the electrical double layer are not applicable. Recently, molecular dynamics simulations and phenomenological theories have attempted to explain the capacitance of the ionic liquid/electrode interface with the molecular-scale structure and dynamics of the ionic liquid near the electrode. rnHowever, experimental evidence is very limited. rnrnIn the presented study, the ion distribution of an ionic liquid near an electrode and its response to applied potentials was examined with sub-molecular resolution. For this purpose, a new sample chamber was constructed, allowing in situ high energy X-ray reflectivity experiments under potential control, as well as impedance spectroscopy measurements. The combination of structural information and electrochmical data provided a comprehensive picture of the electric double layer in ionic liquids. Oscillatory charge density profiles were found, consisting of alternating anion- and cation-enriched layers at both, cathodic and anodic, potentials. This structure was shown to arise from the same ion-ion correlations dominating the liquid bulk structure that were observed as a distinct X-ray diffraction peak. Therefore, existing physically motivated models were refined and verified by comparison with independent measurements. rnrnThe relaxation dynamics of the interfacial structure upon potential variation were studied by time resolved X-ray reflectivity experiments with sub-millisecond resolution. The observed relaxation times during charging/discharging are consistent with the impedance spectroscopy data revealing three processes of vastly different characteristic time-scales. Initially, the ion transport normal to the interface happens on a millisecond-scale. Another 100-millisecond-scale process is associated with molecular reorientation of electrode-adsorbed cations. Further, a minute-scale relaxation was observed, which is tentatively assigned to lateral ordering within the first layer.

Spatial and functional relationships between air conduits and blood capillaries in the pulmonary gas exchange tissue of adult and developing chickens

The documented data regarding the three-dimensional structure of the air capillaries (ACs), the ultimate sites of gas exchange in the avian lung is contradictory. Further, the mode of gas exchange, described as cross-current has not been clearly elucidated. We studied the temporal and spatial arrangement of the terminal air conduits of the chicken lung and their relationship with the blood capillaries (BCs) in embryos as well as the definitive architecture in adults. Several visualization techniques that included corrosion casting, light microscopy as well as scanning and transmission electron microscopy were used. Two to six infundibulae extend from each atrium and give rise to numerous ACs that spread centrifugally. Majority of the ACs are tubular structures that give off branches, which anastomose with their neighboring cognates. Some ACs have globular shapes and a few are blind-ending tapering tubes. During inauguration, the luminal aspects of the ACs are characterized by numerous microvillus-like microplicae, which are formed during the complex processes of cell attenuation and canalization of the ACs. The parabronchial exchange BCs, initially inaugurated as disorganized meshworks, are reoriented via pillar formation to lie predominantly orthogonal to the long axes of the ACs. The remodeling of the retiform meshworks by intussusceptive angiogenesis essentially accomplishes a cross-current system at the gas exchange interface in the adults, where BCs form ring-like patterns around the ACs, thus establishing a cross-current system. Our findings clarify the mode of gas exchange in the parabronchial mantle and illuminate the basis for the functional efficiency of the avian lung.

A Method for Detecting Population Genetic Structure in Diverse, High Gene-Flow Species

Detecting small amounts of genetic subdivision across geographic space remains a persistent challenge. Often a failure to detect genetic structure is mistaken for evidence of panmixia, when more powerful statistical tests may uncover evidence for subtle geographic differentiation. Such slight subdivision can be demographically and evolutionarily important as well as being critical for management decisions. We introduce here a method, called spatial analysis of shared alleles (SAShA), that detects geographically restricted alleles by comparing the spatial arrangement of allelic co-occurrences with the expectation under panmixia. The approach is allele-based and spatially explicit, eliminating the loss of statistical power that can occur with user-defined populations and statistical averaging within populations. Using simulated data sets generated under a stepping-stone model of gene flow, we show that this method outperforms spatial autocorrelation (SA) and UST under common real-world conditions: at relatively high migration rates when diversity is moderate or high, especially when sampling is poor. We then use this method to show clear differences in the genetic patterns of 2 nearshore Pacific mollusks, Tegula funebralis (5 Chlorostoma funebralis) and Katharina tunicata, whose overall patterns of within-species differentiation are similar according to traditional population genetics analyses. SAShA meaningfully complements UST/FST, SA, and other existing geographic genetic analyses and is especially appropriate for evaluating species with high gene flow and subtle genetic differentiation.

Dynamics of Spatial Game Theory Networks with Novel Modifications

Introduction: Advances in biotechnology have shed light on many biological processes. In biological networks, nodes are used to represent the function of individual entities within a system and have historically been studied in isolation. Network structure adds edges that enable communication between nodes. An emerging fieldis to combine node function and network structure to yield network function. One of the most complex networks known in biology is the neural network within the brain. Modeling neural function will require an understanding of networks, dynamics, andneurophysiology. It is with this work that modeling techniques will be developed to work at this complex intersection. Methods: Spatial game theory was developed by Nowak in the context of modeling evolutionary dynamics, or the way in which species evolve over time. Spatial game theory offers a two dimensional view of analyzingthe state of neighbors and updating based on the surroundings. Our work builds upon this foundation by studying evolutionary game theory networks with respect to neural networks. This novel concept is that neurons may adopt a particular strategy that will allow propagation of information. The strategy may therefore act as the mechanism for gating. Furthermore, the strategy of a neuron, as in a real brain, isimpacted by the strategy of its neighbors. The techniques of spatial game theory already established by Nowak are repeated to explain two basic cases and validate the implementation of code. Two novel modifications are introduced in Chapters 3 and 4 that build on this network and may reflect neural networks. Results: The introduction of two novel modifications, mutation and rewiring, in large parametricstudies resulted in dynamics that had an intermediate amount of nodes firing at any given time. Further, even small mutation rates result in different dynamics more representative of the ideal state hypothesized. Conclusions: In both modificationsto Nowak's model, the results demonstrate the network does not become locked into a particular global state of passing all information or blocking all information. It is hypothesized that normal brain function occurs within this intermediate range and that a number of diseases are the result of moving outside of this range.

How does spatial extent of fMRI datasets affect independent component analysis decomposition?

Spatial independent component analysis (sICA) of functional magnetic resonance imaging (fMRI) time series can generate meaningful activation maps and associated descriptive signals, which are useful to evaluate datasets of the entire brain or selected portions of it. Besides computational implications, variations in the input dataset combined with the multivariate nature of ICA may lead to different spatial or temporal readouts of brain activation phenomena. By reducing and increasing a volume of interest (VOI), we applied sICA to different datasets from real activation experiments with multislice acquisition and single or multiple sensory-motor task-induced blood oxygenation level-dependent (BOLD) signal sources with different spatial and temporal structure. Using receiver operating characteristics (ROC) methodology for accuracy evaluation and multiple regression analysis as benchmark, we compared sICA decompositions of reduced and increased VOI fMRI time-series containing auditory, motor and hemifield visual activation occurring separately or simultaneously in time. Both approaches yielded valid results; however, the results of the increased VOI approach were spatially more accurate compared to the results of the decreased VOI approach. This is consistent with the capability of sICA to take advantage of extended samples of statistical observations and suggests that sICA is more powerful with extended rather than reduced VOI datasets to delineate brain activity.

A BAYESIAN HIERARCHICAL FRAMEWORK FOR SPATIAL MODELING OF fMRI DATA

Functional neuroimaging techniques enable investigations into the neural basis of human cognition, emotions, and behaviors. In practice, applications of functional magnetic resonance imaging (fMRI) have provided novel insights into the neuropathophysiology of major psychiatric,neurological, and substance abuse disorders, as well as into the neural responses to their treatments. Modern activation studies often compare localized task-induced changes in brain activity between experimental groups. One may also extend voxel-level analyses by simultaneously considering the ensemble of voxels constituting an anatomically defined region of interest (ROI) or by considering means or quantiles of the ROI. In this work we present a Bayesian extension of voxel-level analyses that offers several notable benefits. First, it combines whole-brain voxel-by-voxel modeling and ROI analyses within a unified framework. Secondly, an unstructured variance/covariance for regional mean parameters allows for the study of inter-regional functional connectivity, provided enough subjects are available to allow for accurate estimation. Finally, an exchangeable correlation structure within regions allows for the consideration of intra-regional functional connectivity. We perform estimation for our model using Markov Chain Monte Carlo (MCMC) techniques implemented via Gibbs sampling which, despite the high throughput nature of the data, can be executed quickly (less than 30 minutes). We apply our Bayesian hierarchical model to two novel fMRI data sets: one considering inhibitory control in cocaine-dependent men and the second considering verbal memory in subjects at high risk for Alzheimer’s disease. The unifying hierarchical model presented in this manuscript is shown to enhance the interpretation content of these data sets.

BIVARIATE BINOMIAL SPATIAL MODELLING LOA loa PREVALENCE IN TROPICAL AFRICA

We present a state-of-the-art application of smoothing for dependent bivariate binomial spatial data to Loa loa prevalence mapping in West Africa. This application is special because it starts with the non-spatial calibration of survey instruments, continues with the spatial model building and assessment and ends with robust, tested software that will be used by the field scientists of the World Health Organization for online prevalence map updating. From a statistical perspective several important methodological issues were addressed: (a) building spatial models that are complex enough to capture the structure of the data but remain computationally usable; (b)reducing the computational burden in the handling of very large covariate data sets; (c) devising methods for comparing spatial prediction methods for a given exceedance policy threshold.

Endothelial survival factors and spatial completion, but not pericyte coverage of retinal capillaries determine vessel plasticity

Pericyte loss and capillary regression are characteristic for incipient diabetic retinopathy. Pericyte recruitment is involved in vessel maturation, and ligand-receptor systems contributing to pericyte recruitment are survival factors for endothelial cells in pericyte-free in vitro systems. We studied pericyte recruitment in relation to the susceptibility toward hyperoxia-induced vascular remodeling using the pericyte reporter X-LacZ mouse and the mouse model of retinopathy of prematurity (ROP). Pericytes were found in close proximity to vessels, both during formation of the superficial and the deep capillary layers. When exposure of mice to the ROP was delayed by 24 h, i.e., after the deep retinal layer had formed [at postnatal (p) day 8], preretinal neovascularizations were substantially diminished at p18. Mice with a delayed ROP exposure had 50% reduced avascular zones. Formation of the deep capillary layers at p8 was associated with a combined up-regulation of angiopoietin-1 and PDGF-B, while VEGF was almost unchanged during the transition from a susceptible to a resistant capillary network. Inhibition of Tie-2 function either by soluble Tie-2 or by a sulindac analog, an inhibitor of Tie-2 phosphorylation, resensitized retinal vessels to neovascularizations due to a reduction of the deep capillary network. Inhibition of Tie-2 function had no effect on pericyte recruitment. Our data indicate that the final maturation of the retinal vasculature and its resistance to regressive signals such as hyperoxia depend on the completion of the multilayer structure, in particular the deep capillary layers, and are independent of the coverage by pericytes.

Fine-scale genetic structure and dispersal in the common vole (Microtus arvalis)

The genetic structure and demography of local populations is tightly linked to the rate and scale of dispersal. Dispersal parameters are notoriously difficult to determine in the field, and remain often completely unknown for smaller organisms. In this study, we investigate spatial and temporal genetic structure in relation to dispersal patterns among local populations of the probably most abundant European mammals, the common vole (Microtus arvalis). Voles were studied in six natural populations at distances of 0.4-2.5 km in three different seasons (fall, spring, summer) corresponding to different life-history stages. Field observations provided no direct evidence for movements of individuals between populations. The analysis of 10 microsatellite markers revealed a persistent overall genetic structure among populations of 2.9%, 2.5% and 3% FST in the respective season. Pairwise comparisons showed that even the closest populations were significantly differentiated from each other in each season, but there was no evidence for temporal differentiation within populations or isolation by distance among populations. Despite significant genetic structure, assignment analyses identified a relatively high proportion of individuals as being immigrants for the population where they were captured. The immigration rate was not significantly lower for females than for males. We suggest that a generally low and sex-dependent effective dispersal rate as the consequence of only few immigrants reproducing successfully in the new populations together with the social structure within populations may explain the maintenance of genetic differentiation among populations despite migration.

Landscape and forest structure at moose mortality sites on Isle Royale National Park. A LiDAR based assessment

Landscape structure and heterogeneity play a potentially important, but little understood role in predator-prey interactions and behaviourally-mediated habitat selection. For example, habitat complexity may either reduce or enhance the efficiency of a predator's efforts to search, track, capture, kill and consume prey. For prey, structural heterogeneity may affect predator detection, avoidance and defense, escape tactics, and the ability to exploit refuges. This study, investigates whether and how vegetation and topographic structure influence the spatial patterns and distribution of moose (Alces alces) mortality due to predation and malnutrition at the local and landscape levels on Isle Royale National Park. 230 locations where wolves (Canis lupus) killed moose during the winters between 2002 and 2010, and 182 moose starvation death sites for the period 1996-2010, were selected from the extensive Isle Royale Wolf-Moose Project carcass database. A variety of LiDAR-derived metrics were generated and used in an algorithm model (Random Forest) to identify, characterize, and classify three-dimensional variables significant to each of the mortality classes. Furthermore, spatial models to predict and assess the likelihood at the landscape scale of moose mortality were developed. This research found that the patterns of moose mortality by predation and malnutrition across the landscape are non-random, have a high degree of spatial variability, and that both mechanisms operate in contexts of comparable physiographic and vegetation structure. Wolf winter hunting locations on Isle Royale are more likely to be a result of its prey habitat selection, although they seem to prioritize the overall areas with higher moose density in the winter. Furthermore, the findings suggest that the distribution of moose mortality by predation is habitat-specific to moose, and not to wolves. In addition, moose sex, age, and health condition also affect mortality site selection, as revealed by subtle differences between sites in vegetation heights, vegetation density, and topography. Vegetation density in particular appears to differentiate mortality locations for distinct classes of moose. The results also emphasize the significance of fine-scale landscape and habitat features when addressing predator-prey interactions. These finer scale findings would be easily missed if analyses were limited to the broader landscape scale alone.

INVESTIGATING THE IMPACTS OF STREAMBED HABITAT HETEROGENEITY ON ECOSYSTEM STRUCTURE AND PROCESSES USING BASIC AND APPLIED PERSPECTIVES

Stream restoration often focuses on increasing habitat heterogeneity to reverse ecosystem degradation. However, the connection between heterogeneity and ecosystem structure and processes is poorly understood. We looked to investigate this interaction from both applied and basic science perspectives. For the applied study, we examined two culvert replacements designed to mimic natural stream channels, to see if they were better at maintaining ecosystem processes within as well as upstream and downstream of culverts compared to non-replaced culverts. We measured three ecosystem processes (nutrient uptake, hydrologic characteristics, and coarse particulate organic matter retention) and found that stream simulation culvert restoration improved organic matter retention within culverts, and that there were no differences in processes measured upstream and downstream of both restoration designs. Our results suggest that measurements of ecosystem processes are more likely to show a response to restoration if they match the scale of the restoration activity. For the basic science study, we quantified the longitudinal spatial heterogeneity of physical and biofilm characteristics at microhabitat to segment scales on streams with different streambed variability. We found that all physical characteristics and biofilm characteristics were spatially independent at the macro-habitat scale and greater. Together, these studies demonstrate the importance of scale in ecological interactions and the value of incorporating considerations of scale into restoration activities.