Hybrid constraint programming and metaheuristic methods for large scale optimization problems

This work presents hybrid Constraint Programming (CP) and metaheuristic methods for the solution of Large Scale Optimization Problems; it aims at integrating concepts and mechanisms from the metaheuristic methods to a CP-based tree search environment in order to exploit the advantages of both approaches. The modeling and solution of large scale combinatorial optimization problem is a topic which has arisen the interest of many researcherers in the Operations Research field; combinatorial optimization problems are widely spread in everyday life and the need of solving difficult problems is more and more urgent. Metaheuristic techniques have been developed in the last decades to effectively handle the approximate solution of combinatorial optimization problems; we will examine metaheuristics in detail, focusing on the common aspects of different techniques. Each metaheuristic approach possesses its own peculiarities in designing and guiding the solution process; our work aims at recognizing components which can be extracted from metaheuristic methods and re-used in different contexts. In particular we focus on the possibility of porting metaheuristic elements to constraint programming based environments, as constraint programming is able to deal with feasibility issues of optimization problems in a very effective manner. Moreover, CP offers a general paradigm which allows to easily model any type of problem and solve it with a problem-independent framework, differently from local search and metaheuristic methods which are highly problem specific. In this work we describe the implementation of the Local Branching framework, originally developed for Mixed Integer Programming, in a CP-based environment. Constraint programming specific features are used to ease the search process, still mantaining an absolute generality of the approach. We also propose a search strategy called Sliced Neighborhood Search, SNS, that iteratively explores slices of large neighborhoods of an incumbent solution by performing CP-based tree search and encloses concepts from metaheuristic techniques. SNS can be used as a stand alone search strategy, but it can alternatively be embedded in existing strategies as intensification and diversification mechanism. In particular we show its integration within the CP-based local branching. We provide an extensive experimental evaluation of the proposed approaches on instances of the Asymmetric Traveling Salesman Problem and of the Asymmetric Traveling Salesman Problem with Time Windows. The proposed approaches achieve good results on practical size problem, thus demonstrating the benefit of integrating metaheuristic concepts in CP-based frameworks.

Synthetical engineering of supramolecular properties of large polycyclic aromatic hydrocarbons

An important property for devices is the charge-carrier mobility values for discotic organic materials like hexa-peri-hexabenzocoronenes. A close relation exists between the degree of their columnar self-arrangement of the molecules and their mobilities. Within this first step an induction of a higher order via hydrogen-bonding was considered, which mainly pointed towards the improvement of the intracolumnar stacking of the materials. For the analytics a broad range of methods was used including differential scanning calorimetry (DSC), wide-angle X-ray diffractometry (WAXS), solid-state NMR spectroscopy and scanning tunneling microscopy (STM). Indeed, a specific influence of the hydrogen-bonds could be identified, although in several cases by the cost of a severe reduction of solubility and processability. This effect was dampened by the addition of a long alkyl chain next to the hydrogen-bond exerting functional group, which resulted in an improved columnar arrangement by retention of processability. In contrast to the before mentioned example of inducing a higher intracolumnar order by hydrogen-bonding, the focus was also be set upon larger aromatic systems. The charge-carrier mobility is also in close relation to the size of the aromatic core and larger π-areas are expected to lead to improved mobilities. For photovoltaic applications a high extinction coefficient over a broad range of the spectrum is favorable, which can also be achieved by enlarging the aromatic core component. In addition the stronger π-interactions between the aromatic core components should yield an improved columnar stability and order. However the strengthening of the π-interactions between the aromatic core components led to a reduction of the solubility and the processability due to the stronger aggregation of the molecules. This required the introduction of efficiently solubilizing features in terms of long alkyl chains in the corona of the aromatic entity, in combination of a distortion of the aromatic core moiety by bulky tert-butyl groups. By this approach not only the processing and cleaning of the materials with standard laboratory techniques became possible, but moreover the first structure-rich UV/vis and a resolved 1H-NMR spectra for an aromatic system two times larger than hexa-peri-hexabenzocoronene were recorded. The bulk properties in an extruded fiber as well as on the surface showed a columnar self-assembly including a phase in which a homeotropic alignment on a substrate was observed, which turns the material into an interesting candidate for future applications in electronic devices.

A molecular approach towards tethered bilayer lipid membranes: synthesis and characterization of novel anchor lipids

Tethered bilayer lipid membranes (tBLMs) are a promising model system for the natural cell membrane. They consist of a lipid bilayer that is covalently coupled to a solid support via a spacer group. In this study, we developed a suitable approach to increase the submembrane space in tBLMs. The challenge is to create a membrane with a lower lipid density in order to increase the membrane fluidity, but to avoid defects that might appear due to an increase in the lateral space within the tethered monolayers. Therefore, various synthetic strategies and different monolayer preparation techniques were examined. Synthetical attempts to achieve a large ion reservoir were made in two directions: increasing the spacer length of the tether lipids and increasing the lateral distribution of the lipids in the monolayer. The first resulted in the synthesis of a small library of tether lipids (DPTT, DPHT and DPOT) characterized by 1H and 13C NMR, FD-MS, ATR, DSC and TGA. The synthetic strategy for their preparation includes synthesis of precursor with a double bond anchor that can be easily modified for different substrates (e.g. metal and metaloxide). Here, the double bond was modified into a thiol group suitable for gold surface. Another approach towards the preparation of homogeneous monolayers with decreased two-dimensional packing density was the synthesis of two novel anchor lipids: DPHDL and DDPTT. DPHDL is “self-diluted” tether lipid containing two lipoic anchor moieties. DDPTT has an extended lipophylic part that should lead to the preparation of diluted, leakage free proximal layers that will facilitate the completion of the bilayer. Our tool-box of tether lipids was completed with two fluorescent labeled lipid precursors with respectively one and two phytanyl chains in the hydrophobic region and a dansyl group as a fluorophore. The use of such fluorescently marked lipids is supposed to give additional information for the lipid distribution on the air-water interface. The Langmuir film balance was used to investigate the monolayer properties of four of the synthesized thiolated anchor lipids. The packing density and mixing behaviour were examined. The results have shown that mixing anchor with free lipids can homogeneously dilute the anchor lipid monolayers. Moreover, an increase in the hydrophylicity (PEG chain length) of the anchor lipids leads to a higher packing density. A decrease in the temperature results in a similar trend. However, increasing the number of phytanyl chains per lipid molecule is shown to decrease the packing density. LB-monolayers based on pure and mixed lipids in different ratio and transfer pressure were tested to form tBLMs with diluted inner layers. A combination of the LB-monolayer transfer with the solvent exchange method accomplished successfully the formation of tBLMs based on pure DPOT. Some preliminary investigations of the electrical sealing properties and protein incorporation of self-assembled DPOT and DDPTT-based tBLMs were conducted. The bilayer formation performed by solvent exchange resulted in membranes with high resistances and low capacitances. The appearance of space beneath the membrane is clearly visible in the impedance spectra expressed by a second RC element. The latter brings the conclusion that the longer spacer in DPOT and the bigger lateral space between the DDPTT molecules in the investigated systems essentially influence the electrical parameters of the membrane. Finally, we could show the functional incorporation of the small ion carrier valinomycin in both types of membranes.

From pyrene to large polycyclic aromatic hydrocarbons

Studies of polycyclic aromatic hydrocarbons have shown that the overall size, periphery, and functionalization of PAHs are crucial parameters which significantly alter their electronic structure and chemical reactivity. Therefore, the major direction of this thesis is the synthesis and characterization of extended PAHs: (i) with different functional groups improving their processability, (ii) with different periphery changing their chemical reactivity, (iii) with inclusions of different metal ions, which influence their physical properties. • The cyclodehydrogenation reaction has been proposed for to synthesise polyphenylene ribbons with preplanarized (dibenzo[e,l]pyrene) moieties in the aromatic core with up to 10 nm linear size. The synthetic strategy employed is discussed in Chapter 2 and is based on stoichiometrically controlled DIELS-ALDER cycloaddition. All molecules possessed very good solubility in common organic solvents allowing their characterization by standard analytical techniques. • A new concept was developed to extend PAH’s core. Here the introduction of “zigzag” sites, discussed in Chpater 3 was shown to lower the HOMO-LUMO gap and to result in higher chemical reactivities. This allowed, in Chapters 3, 4 and 5, further functionalization of PAH and enlargement of their aromatic cores up to 224 atoms. Despite the size of these novel molecules, extraordinary solubilities in common organic solvents were obtained due to distortions from planarity of the aromatic cores by bulky tert-butyl groups, which hampered the usually very pronounced aggregation tendency of extended π-systems. All extended PAHs posses the small HOMO-LUMO gap together with good electron affinities making them potential candidates for application in organic FETs. • Another alternative synthetic route has been proposed to obtain extended the metal-PAH complexes. Using the quinoxaline methodology in Chapter 5 three new phenanthroline ligands (up to 60 skeletal atoms) have been synthesized and characterized. Four different (Ru(II), Cu(II) and Pt(II)) complexes were synthesized, allowing to construct a range of large metal complexes by varying the metal as well as the number and nature of ligands.

Heteroatom containing polycyclic aromatic hydrocarbons with positive charge - synthesis and characterization

The synthesis and characterization of various heteroatom containing PAHs with positive charge were investigated in this work: 1. A series of 2-phenyl-benzo[8,9]quinolizino[4,5,6,7-fed]phenanthridinylium (PQP) salts with different alkyl chains and anions were synthesized. The synthesis of the extended derivates of PQP salts with two fused benzene rings, 2-phenyl-naphthacene[1,2]quinolizino[3,4,5,6-def]benzo[i]phenanthridinium (DBPQP) tetrafluoroborate was also developed. The self-assembly behavior of these amphiphilic PAHs was investigated in methanolic solution as well as in the bulk. Various aggregates with different morphologies such as fibers, tubes and vesicals were obtained from their solution. All of these morphology changes could be ascribed to the changes in intermolecular interactions which resulting from the difference in the molecular structures such as aromatic cores, alkyl chains and counterions. 2. The synthetic strategy of oxygen containing positively charged PAHs, benzo[5,6]naphthaceno[1,12,11,10-jklmna]xanthylium (BNAX) salts and its dibenzo derivates, DBNAX salts were developed. With a similar method, sulfur containing benzo[5,6]naphthaceno[1,12,11,10-jklmna]thioxanthylium (BNATX) salts were also synthesized. Various BNAX salts with different alkyl chains could be obtained and their supramolecular behavior were investigated. A discotic liquid crystalline behavior was observed for di- (3-25) and tridodecyl (3-27) substituted BNAX salts and both compounds exhibited large unit cell in their 2D-WAXS patterns which could be attributed to the formation of dimer structures. By drop casting their methanolic solution on silicon wafers, similar nanoscaled fibers from monododecyl substituted BNAX bromide 3-24 and DBNAX bromide 3-35 could be observed. 3. A novel synthetic method toward nitrogen containing 14-phenyl-dibenzo[jk,mn]naphtho[2,1,8-fgh]thebenidinium (DBNT) salts was also developed. In this method, the undehydrogenated precursor of DBNT, dibenzoacridinium salt could be produced directly from the reaction between dibenzoxanthenylium derivates and amine/aniline in reasonable yields. Various DBNT salts with different alkyl and alkylphenyl chains on their nitrogen atom were synthesized in this two-step method. The self-assembly behavior of two alkylated DBNT salts, 4-15a and 4-18b was also studied in this work. Compound 4-15a formed nanoscaled fibers and helical aggregates were obtained from 4-18b in their methanolic solutions. 4. Various ionic complexes were derived by complexing PQP and DBPQP cations with different sulfate/sulfonate group containing anionic surfactants. The ionic complexes resulting from the ionic self-assembly (ISA) method exhibited self-assembly behavior which was controllable by the species and shape of cations and anions. Various aggregates such as nanofibers and spherical aggregates could be produced from their methanolic solution in a defined manner conveniently.

Evolutionary tendencies in African Marantaceae

Die Marantaceae (550 Arten) sind eine weltweit verbreitete Familie von Stauden und Lianen aus dem Unterwuchs tropischer Tieflandregenwälder. Der morphologisch-ökologische Vergleich des basal abzweigenden Sarcophrynium-Astes mit dem in abgeleiteter Position stehenden Marantochloa-Ast, soll beispielhaft evolutionäre Muster in der Familie beleuchten. So wird in der Doktorarbeit zum ersten Mal ein Überblick über die Blütenbiologie und Phylogenie von rund 30 der 40 afrikanischen Marantaceae Arten präsentiert. Die Analysen basieren auf Daten von drei mehrmonatigen Feldaufenthalten in Gabun jeweils zwischen September und Januar. Vier Blütentypen werden beschrieben, die jeweils mit einer spezifischen Bestäubergilde verbunden sind (kleine, mittlere, große Bienen bzw. Vögel). Bestäubungsexperimente belegen, dass 18 Arten selbstkompatibel, aber nur zwei Arten autogam sind, also keine Bestäubungsvermittler benötigen. Der Fruchtansatz ist generell gering (10 -30 %). Die komplexe Synorganisation der Blüte ermöglicht in den Marantaceae einen explosiven Bestäubungsmechanismus. Um dessen ökologische Funktionalität zu verstehen, werden die Blüten von 66 Arten, alle wichtigen Äste der Marantaceae abdeckend, unter einem morphologisch-funktionalen Gesichtspunkt untersucht. Es gibt große Übereinstimmungen zwischen allen untersuchten Arten im Zusammenspiel (Synorganisation) der wichtigsten Bauelemente (Griffel, Kapuzenblatt, Schwielenblatt), die eine präzise Pollenübertragung ermöglichen. Basierend auf Daten von nrDNA (ITS, 5S) und cpDNA (trnL-F) wird für ein nahezu komplettes Artenspektrum die Phylogenie der zwei afrikanischen Äste erstellt. Hierauf werden morphologische und ökologische Merkmale sowie geographischer Verbreitungsmuster nach dem Parsimonieprinzip rekonstruiert, um so deren evolutionäre Bedeutung für die Marantaceae abschätzen zu können. Die Ergebnisse weisen auf die Beteiligung einer Vielzahl verschiedener Artbildungsfaktoren hin.

Large-scale coupled-cluster calculations

Coupled-cluster theory provides one of the most successful concepts in electronic-structure theory. This work covers the parallelization of coupled-cluster energies, gradients, and second derivatives and its application to selected large-scale chemical problems, beside the more practical aspects such as the publication and support of the quantum-chemistry package ACES II MAB and the design and development of a computational environment optimized for coupled-cluster calculations. The main objective of this thesis was to extend the range of applicability of coupled-cluster models to larger molecular systems and their properties and therefore to bring large-scale coupled-cluster calculations into day-to-day routine of computational chemistry. A straightforward strategy for the parallelization of CCSD and CCSD(T) energies, gradients, and second derivatives has been outlined and implemented for closed-shell and open-shell references. Starting from the highly efficient serial implementation of the ACES II MAB computer code an adaptation for affordable workstation clusters has been obtained by parallelizing the most time-consuming steps of the algorithms. Benchmark calculations for systems with up to 1300 basis functions and the presented applications show that the resulting algorithm for energies, gradients and second derivatives at the CCSD and CCSD(T) level of theory exhibits good scaling with the number of processors and substantially extends the range of applicability. Within the framework of the ’High accuracy Extrapolated Ab initio Thermochemistry’ (HEAT) protocols effects of increased basis-set size and higher excitations in the coupled- cluster expansion were investigated. The HEAT scheme was generalized for molecules containing second-row atoms in the case of vinyl chloride. This allowed the different experimental reported values to be discriminated. In the case of the benzene molecule it was shown that even for molecules of this size chemical accuracy can be achieved. Near-quantitative agreement with experiment (about 2 ppm deviation) for the prediction of fluorine-19 nuclear magnetic shielding constants can be achieved by employing the CCSD(T) model together with large basis sets at accurate equilibrium geometries if vibrational averaging and temperature corrections via second-order vibrational perturbation theory are considered. Applying a very similar level of theory for the calculation of the carbon-13 NMR chemical shifts of benzene resulted in quantitative agreement with experimental gas-phase data. The NMR chemical shift study for the bridgehead 1-adamantyl cation at the CCSD(T) level resolved earlier discrepancies of lower-level theoretical treatment. The equilibrium structure of diacetylene has been determined based on the combination of experimental rotational constants of thirteen isotopic species and zero-point vibrational corrections calculated at various quantum-chemical levels. These empirical equilibrium structures agree to within 0.1 pm irrespective of the theoretical level employed. High-level quantum-chemical calculations on the hyperfine structure parameters of the cyanopolyynes were found to be in excellent agreement with experiment. Finally, the theoretically most accurate determination of the molecular equilibrium structure of ferrocene to date is presented.

Life Cycle Assessment of Peach Nectar: a comparative analysis between conventional and bioenergy-from-waste integrated food chains

Modern food systems are characterized by a high energy intensity as well as by the production of large amounts of waste, residuals and food losses. This inefficiency presents major consequences, in terms of GHG emissions, waste disposal, and natural resource depletion. The research hypothesis is that residual biomass material could contribute to the energetic needs of food systems, if recovered as an integrated renewable energy source (RES), leading to a sensitive reduction of the impacts of food systems, primarily in terms of fossil fuel consumption and GHG emissions. In order to assess these effects, a comparative life cycle assessment (LCA) has been conducted to compare two different food systems: a fossil fuel-based system and an integrated system with the use of residual as RES for self-consumption. The food product under analysis has been the peach nectar, from cultivation to end-of-life. The aim of this LCA is twofold. On one hand, it allows an evaluation of the energy inefficiencies related to agro-food waste. On the other hand, it illustrates how the integration of bioenergy into food systems could effectively contribute to reduce this inefficiency. Data about inputs and waste generated has been collected mainly through literature review and databases. Energy balance, GHG emissions (Global Warming Potential) and waste generation have been analyzed in order to identify the relative requirements and contribution of the different segments. An evaluation of the energy “loss” through the different categories of waste allowed to provide details about the consequences associated with its management and/or disposal. Results should provide an insight of the impacts associated with inefficiencies within food systems. The comparison provides a measure of the potential reuse of wasted biomass and the amount of energy recoverable, that could represent a first step for the formulation of specific policies on the integration of bioenergies for self-consumption.

Spatial Disparity and Dynamic Trajectory of Convergence in China: Construction and Application of a Composite Index of Regional Development

China is a large country characterized by remarkable growth and distinct regional diversity. Spatial disparity has always been a hot issue since China has been struggling to follow a balanced growth path but still confronting with unprecedented pressures and challenges. To better understand the inequality level benchmarking spatial distributions of Chinese provinces and municipalities and estimate dynamic trajectory of sustainable development in China, I constructed the Composite Index of Regional Development (CIRD) with five sub pillars/dimensions involving Macroeconomic Index (MEI), Science and Innovation Index (SCI), Environmental Sustainability Index (ESI), Human Capital Index (HCI) and Public Facilities Index (PFI), endeavoring to cover various fields of regional socioeconomic development. Ranking reports on the five sub dimensions and aggregated CIRD were provided in order to better measure the developmental degrees of 31 or 30 Chinese provinces and municipalities over 13 years from 1998 to 2010 as the time interval of three “Five-year Plans”. Further empirical applications of this CIRD focused on clustering and convergence estimation, attempting to fill up the gap in quantifying the developmental levels of regional comprehensive socioeconomics and estimating the dynamic convergence trajectory of regional sustainable development in a long run. Four clusters were benchmarked geographically-oriented in the map on the basis of cluster analysis, and club-convergence was observed in the Chinese provinces and municipalities based on stochastic kernel density estimation.

Genome-wide investigation reveals high evolutionary rates in annual model plants

Background: ;Rates of molecular evolution vary widely among species. While significant deviations from molecular clock have been found in many taxa, effects of life histories on molecular evolution are not fully understood. In plants, annual/perennial life history traits have long been suspected to influence the evolutionary rates at the molecular level. To date, however, the number of genes investigated on this subject is limited and the conclusions are mixed. To evaluate the possible heterogeneity in evolutionary rates between annual and perennial plants at the genomic level, we investigated 85 nuclear housekeeping genes, 10 non-housekeeping families, and 34 chloroplast;genes using the genomic data from model plants including Arabidopsis thaliana and Medicago truncatula for annuals and grape (Vitis vinifera) and popular (Populus trichocarpa) for perennials.;Results: ;According to the cross-comparisons among the four species, 74-82% of the nuclear genes and 71-97% of the chloroplast genes suggested higher rates of molecular evolution in the two annuals than those in the two perennials. The significant heterogeneity in evolutionary rate between annuals and perennials was consistently found both in nonsynonymous sites and synonymous sites. While a linear correlation of evolutionary rates in orthologous genes between species was observed in nonsynonymous sites, the correlation was weak or invisible in synonymous sites. This tendency was clearer in nuclear genes than in chloroplast genes, in which the overall;evolutionary rate was small. The slope of the regression line was consistently lower than unity, further confirming the higher evolutionary rate in annuals at the genomic level.;Conclusions: ;The higher evolutionary rate in annuals than in perennials appears to be a universal phenomenon both in nuclear and chloroplast genomes in the four dicot model plants we investigated. Therefore, such heterogeneity in evolutionary rate should result from factors that have genome-wide influence, most likely those associated with annual/perennial life history. Although we acknowledge current limitations of this kind of study, mainly due to a small sample size available and a distant taxonomic relationship of the model organisms, our results indicate that the genome-wide survey is a promising approach toward further understanding of the;mechanism determining the molecular evolutionary rate at the genomic level.

Targeting siglecs--a novel pharmacological strategy for immuno- and glycotherapy

The immune system must be tightly held in check to avoid bystander tissue damage as well as autoreactivity caused by overwhelming immune reactions. A novel family of immunoregulatory, carbohydrate-binding receptors, the Siglecs (sialic acid binding immunoglobulin-like lectins), has received particular attention in light of their capacity to mediate cell death, anti-proliferative effects and to regulate a variety of cellular activities. Siglec receptors are mainly expressed on leukocytes in a cell type-specific and differentiation-dependent manner. Siglecs might potentially be exploited as targets of novel immune- and glycotherapeutics for cell-directed therapies in autoimmune and allergic diseases, as well as in haematologic malignancies. Here we present novel insights on structural and functional characteristics, expression patterns and evolutionary aspects of Siglecs and their ligands. Pharmacological strategies using Siglec agonistic cross-linking therapeutics, such as monoclonal or engineered antibodies, intravenous immunoglobulin (IVIG), or glycomimetics are discussed. Modulation of immune responses by targeting Siglecs using agonistic or antagonistic therapeutics may have important clinical implications and may pave the way for novel pharmacological avenues for the treatment of autoimmune and allergic diseases or for tumor immunotherapy.

Individual Differences in First and Fourth Year College Women’s Short Term Mating Strategy Preferences and Perceptions

Using survey methodology, a cross sectional study was undertaken to ascertain whether first and fourth year college women have different perceptions and behavior associated with short term mating preferences. It was hypothesized that after incurring significant negative or costly experiences associated with hooking up, fourth year women would prefer men who had qualities associated with a desired long term partner as opposed to characteristics associated with short term mating partners. The results were partially consistent with the hypothesis. Reported preferences in a desired partner and perspective on hooking up differ between first and fourth year groups. No difference was found between frequency and willingness to hookup between the two groups. The findings are explained in terms of evolutionary theory, social exchange theory, and sexual script concepts.

Individual differences in first and fourth year college women’s short term mating strategy preferences and perceptions

Using survey methodology, a cross sectional study was undertaken to ascertain whether first and fourth year college women have different perceptions and behavior associated with short term mating preferences. It was hypothesized that after incurring significant negative or costly experiences associated with hooking up, fourth year women would prefer men who had qualities associated with a desired long term partner as opposed to characteristics associated with short term mating partners. The results were partially consistent with the hypothesis. Reported preferences in a desired partner and perspective on hooking up differ between first and fourth year groups. No difference was found between frequency and willingness to hookup between the two groups. The findings are explained in terms of evolutionary theory, social exchange theory, and sexual script concepts.

Citizen science reveals unexpected continental-scale evolutionary change in a model organism

Organisms provide some of the most sensitive indicators of climate change and evolutionary responses are becoming apparent in species with short generation times. Large datasets on genetic polymorphism that can provide an historical benchmark against which to test for recent evolutionary responses are very rare, but an exception is found in the brown-lipped banded snail (Cepaea nemoralis). This species is sensitive to its thermal environment and exhibits several polymorphisms of shell colour and banding pattern affecting shell albedo in the majority of populations within its native range in Europe. We tested for evolutionary changes in shell albedo that might have been driven by the warming of the climate in Europe over the last half century by compiling an historical dataset for 6,515 native populations of C. nemoralis and comparing this with new data on nearly 3,000 populations. The new data were sampled mainly in 2009 through the Evolution MegaLab, a citizen science project that engaged thousands of volunteers in 15 countries throughout Europe in the biggest such exercise ever undertaken. A known geographic cline in the frequency of the colour phenotype with the highest albedo (yellow) was shown to have persisted and a difference in colour frequency between woodland and more open habitats was confirmed, but there was no general increase in the frequency of yellow shells. This may have been because snails adapted to a warming climate through behavioural thermoregulation. By contrast, we detected an unexpected decrease in the frequency of Unbanded shells and an increase in the Mid-banded morph. Neither of these evolutionary changes appears to be a direct response to climate change, indicating that the influence of other selective agents, possibly related to changing predation pressure and habitat change with effects on micro-climate.

Spatial and temporal variation in population genetic structure of wild Nile tilapia (Oreochromis niloticus) across Africa

Background: Reconstructing the evolutionary history of a species is challenging. It often depends not only on the past biogeographic and climatic events but also the contemporary and ecological factors, such as current connectivity and habitat heterogeneity. In fact, these factors might interact with each other and shape the current species distribution. However, to what extent the current population genetic structure reflects the past and the contemporary factors is largely unknown. Here we investigated spatio-temporal genetic structures of Nile tilapia (Oreochromis niloticus) populations, across their natural distribution in Africa. While its large biogeographic distribution can cause genetic differentiation at the paleo-biogeographic scales, its restricted dispersal capacity might induce a strong genetic structure at micro-geographic scales. Results: Using nine microsatellite loci and 350 samples from ten natural populations, we found the highest genetic differentiation among the three ichthyofaunal provinces and regions (Ethiopian, Nilotic and Sudano-Sahelian) (R(ST) = 0.38 - 0.69). This result suggests the predominant effect of paleo-geographic events at macro-geographic scale. In addition, intermediate divergences were found between rivers and lakes within the regions, presumably reflecting relatively recent interruptions of gene flow between hydrographic basins (R(ST) = 0.24 - 0.32). The lowest differentiations were observed among connected populations within a basin (R(ST) = 0.015 in the Volta basin). Comparison of temporal sample series revealed subtle changes in the gene pools in a few generations (F = 0 - 0.053). The estimated effective population sizes were 23 - 143 and the estimated migration rate was moderate (m similar to 0.094 - 0.097) in the Volta populations. Conclusions: This study revealed clear hierarchical patterns of the population genetic structuring of O. niloticus in Africa. The effects of paleo-geographic and climatic events were predominant at macro-geographic scale, and the significant effect of geographic connectivity was detected at micro-geographic scale. The estimated effective population size, the moderate level of dispersal and the rapid temporal change in genetic composition might reflect a potential effect of life history strategy on population dynamics. This hypothesis deserves further investigation. The dynamic pattern revealed at micro-geographic and temporal scales appears important from a genetic resource management as well as from a biodiversity conservation point of view.