EADock : design of a new molecular docking algorithm and some of its applications

3 Summary 3. 1 English The pharmaceutical industry has been facing several challenges during the last years, and the optimization of their drug discovery pipeline is believed to be the only viable solution. High-throughput techniques do participate actively to this optimization, especially when complemented by computational approaches aiming at rationalizing the enormous amount of information that they can produce. In siiico techniques, such as virtual screening or rational drug design, are now routinely used to guide drug discovery. Both heavily rely on the prediction of the molecular interaction (docking) occurring between drug-like molecules and a therapeutically relevant target. Several softwares are available to this end, but despite the very promising picture drawn in most benchmarks, they still hold several hidden weaknesses. As pointed out in several recent reviews, the docking problem is far from being solved, and there is now a need for methods able to identify binding modes with a high accuracy, which is essential to reliably compute the binding free energy of the ligand. This quantity is directly linked to its affinity and can be related to its biological activity. Accurate docking algorithms are thus critical for both the discovery and the rational optimization of new drugs. In this thesis, a new docking software aiming at this goal is presented, EADock. It uses a hybrid evolutionary algorithm with two fitness functions, in combination with a sophisticated management of the diversity. EADock is interfaced with .the CHARMM package for energy calculations and coordinate handling. A validation was carried out on 37 crystallized protein-ligand complexes featuring 11 different proteins. The search space was defined as a sphere of 15 R around the center of mass of the ligand position in the crystal structure, and conversely to other benchmarks, our algorithms was fed with optimized ligand positions up to 10 A root mean square deviation 2MSD) from the crystal structure. This validation illustrates the efficiency of our sampling heuristic, as correct binding modes, defined by a RMSD to the crystal structure lower than 2 A, were identified and ranked first for 68% of the complexes. The success rate increases to 78% when considering the five best-ranked clusters, and 92% when all clusters present in the last generation are taken into account. Most failures in this benchmark could be explained by the presence of crystal contacts in the experimental structure. EADock has been used to understand molecular interactions involved in the regulation of the Na,K ATPase, and in the activation of the nuclear hormone peroxisome proliferatoractivated receptors a (PPARa). It also helped to understand the action of common pollutants (phthalates) on PPARy, and the impact of biotransformations of the anticancer drug Imatinib (Gleevec®) on its binding mode to the Bcr-Abl tyrosine kinase. Finally, a fragment-based rational drug design approach using EADock was developed, and led to the successful design of new peptidic ligands for the a5ß1 integrin, and for the human PPARa. In both cases, the designed peptides presented activities comparable to that of well-established ligands such as the anticancer drug Cilengitide and Wy14,643, respectively. 3.2 French Les récentes difficultés de l'industrie pharmaceutique ne semblent pouvoir se résoudre que par l'optimisation de leur processus de développement de médicaments. Cette dernière implique de plus en plus. de techniques dites "haut-débit", particulièrement efficaces lorsqu'elles sont couplées aux outils informatiques permettant de gérer la masse de données produite. Désormais, les approches in silico telles que le criblage virtuel ou la conception rationnelle de nouvelles molécules sont utilisées couramment. Toutes deux reposent sur la capacité à prédire les détails de l'interaction moléculaire entre une molécule ressemblant à un principe actif (PA) et une protéine cible ayant un intérêt thérapeutique. Les comparatifs de logiciels s'attaquant à cette prédiction sont flatteurs, mais plusieurs problèmes subsistent. La littérature récente tend à remettre en cause leur fiabilité, affirmant l'émergence .d'un besoin pour des approches plus précises du mode d'interaction. Cette précision est essentielle au calcul de l'énergie libre de liaison, qui est directement liée à l'affinité du PA potentiel pour la protéine cible, et indirectement liée à son activité biologique. Une prédiction précise est d'une importance toute particulière pour la découverte et l'optimisation de nouvelles molécules actives. Cette thèse présente un nouveau logiciel, EADock, mettant en avant une telle précision. Cet algorithme évolutionnaire hybride utilise deux pressions de sélections, combinées à une gestion de la diversité sophistiquée. EADock repose sur CHARMM pour les calculs d'énergie et la gestion des coordonnées atomiques. Sa validation a été effectuée sur 37 complexes protéine-ligand cristallisés, incluant 11 protéines différentes. L'espace de recherche a été étendu à une sphère de 151 de rayon autour du centre de masse du ligand cristallisé, et contrairement aux comparatifs habituels, l'algorithme est parti de solutions optimisées présentant un RMSD jusqu'à 10 R par rapport à la structure cristalline. Cette validation a permis de mettre en évidence l'efficacité de notre heuristique de recherche car des modes d'interactions présentant un RMSD inférieur à 2 R par rapport à la structure cristalline ont été classés premier pour 68% des complexes. Lorsque les cinq meilleures solutions sont prises en compte, le taux de succès grimpe à 78%, et 92% lorsque la totalité de la dernière génération est prise en compte. La plupart des erreurs de prédiction sont imputables à la présence de contacts cristallins. Depuis, EADock a été utilisé pour comprendre les mécanismes moléculaires impliqués dans la régulation de la Na,K ATPase et dans l'activation du peroxisome proliferatoractivated receptor a (PPARa). Il a également permis de décrire l'interaction de polluants couramment rencontrés sur PPARy, ainsi que l'influence de la métabolisation de l'Imatinib (PA anticancéreux) sur la fixation à la kinase Bcr-Abl. Une approche basée sur la prédiction des interactions de fragments moléculaires avec protéine cible est également proposée. Elle a permis la découverte de nouveaux ligands peptidiques de PPARa et de l'intégrine a5ß1. Dans les deux cas, l'activité de ces nouveaux peptides est comparable à celles de ligands bien établis, comme le Wy14,643 pour le premier, et le Cilengitide (PA anticancéreux) pour la seconde.

An Automatic Approach for Learning and Tuning Gaussian Interval Type-2 Fuzzy Membership Functions Applied to Lung CAD Classification System

The potential of type-2 fuzzy sets for managing high levels of uncertainty in the subjective knowledge of experts or of numerical information has focused on control and pattern classification systems in recent years. One of the main challenges in designing a type-2 fuzzy logic system is how to estimate the parameters of type-2 fuzzy membership function (T2MF) and the Footprint of Uncertainty (FOU) from imperfect and noisy datasets. This paper presents an automatic approach for learning and tuning Gaussian interval type-2 membership functions (IT2MFs) with application to multi-dimensional pattern classification problems. T2MFs and their FOUs are tuned according to the uncertainties in the training dataset by a combination of genetic algorithm (GA) and crossvalidation techniques. In our GA-based approach, the structure of the chromosome has fewer genes than other GA methods and chromosome initialization is more precise. The proposed approach addresses the application of the interval type-2 fuzzy logic system (IT2FLS) for the problem of nodule classification in a lung Computer Aided Detection (CAD) system. The designed IT2FLS is compared with its type-1 fuzzy logic system (T1FLS) counterpart. The results demonstrate that the IT2FLS outperforms the T1FLS by more than 30% in terms of classification accuracy.

EADock: docking of small molecules into protein active sites with a multiobjective evolutionary optimization.

In recent years, protein-ligand docking has become a powerful tool for drug development. Although several approaches suitable for high throughput screening are available, there is a need for methods able to identify binding modes with high accuracy. This accuracy is essential to reliably compute the binding free energy of the ligand. Such methods are needed when the binding mode of lead compounds is not determined experimentally but is needed for structure-based lead optimization. We present here a new docking software, called EADock, that aims at this goal. It uses an hybrid evolutionary algorithm with two fitness functions, in combination with a sophisticated management of the diversity. EADock is interfaced with the CHARMM package for energy calculations and coordinate handling. A validation was carried out on 37 crystallized protein-ligand complexes featuring 11 different proteins. The search space was defined as a sphere of 15 A around the center of mass of the ligand position in the crystal structure, and on the contrary to other benchmarks, our algorithm was fed with optimized ligand positions up to 10 A root mean square deviation (RMSD) from the crystal structure, excluding the latter. This validation illustrates the efficiency of our sampling strategy, as correct binding modes, defined by a RMSD to the crystal structure lower than 2 A, were identified and ranked first for 68% of the complexes. The success rate increases to 78% when considering the five best ranked clusters, and 92% when all clusters present in the last generation are taken into account. Most failures could be explained by the presence of crystal contacts in the experimental structure. Finally, the ability of EADock to accurately predict binding modes on a real application was illustrated by the successful docking of the RGD cyclic pentapeptide on the alphaVbeta3 integrin, starting far away from the binding pocket.

Effects of sample size on the performance of species distribution models

A wide range of modelling algorithms is used by ecologists, conservation practitioners, and others to predict species ranges from point locality data. Unfortunately, the amount of data available is limited for many taxa and regions, making it essential to quantify the sensitivity of these algorithms to sample size. This is the first study to address this need by rigorously evaluating a broad suite of algorithms with independent presence-absence data from multiple species and regions. We evaluated predictions from 12 algorithms for 46 species (from six different regions of the world) at three sample sizes (100, 30, and 10 records). We used data from natural history collections to run the models, and evaluated the quality of model predictions with area under the receiver operating characteristic curve (AUC). With decreasing sample size, model accuracy decreased and variability increased across species and between models. Novel modelling methods that incorporate both interactions between predictor variables and complex response shapes (i.e. GBM, MARS-INT, BRUTO) performed better than most methods at large sample sizes but not at the smallest sample sizes. Other algorithms were much less sensitive to sample size, including an algorithm based on maximum entropy (MAXENT) that had among the best predictive power across all sample sizes. Relative to other algorithms, a distance metric algorithm (DOMAIN) and a genetic algorithm (OM-GARP) had intermediate performance at the largest sample size and among the best performance at the lowest sample size. No algorithm predicted consistently well with small sample size (n < 30) and this should encourage highly conservative use of predictions based on small sample size and restrict their use to exploratory modelling.

Reciprocal Effects on Neurocognitive and Metabolic Phenotypes in Mouse Models of 16p11.2 Deletion and Duplication Syndromes.

The 16p11.2 600 kb BP4-BP5 deletion and duplication syndromes have been associated with developmental delay; autism spectrum disorders; and reciprocal effects on the body mass index, head circumference and brain volumes. Here, we explored these relationships using novel engineered mouse models carrying a deletion (Del/+) or a duplication (Dup/+) of the Sult1a1-Spn region homologous to the human 16p11.2 BP4-BP5 locus. On a C57BL/6N inbred genetic background, Del/+ mice exhibited reduced weight and impaired adipogenesis, hyperactivity, repetitive behaviors, and recognition memory deficits. In contrast, Dup/+ mice showed largely opposite phenotypes. On a F1 C57BL/6N × C3B hybrid genetic background, we also observed alterations in social interaction in the Del/+ and the Dup/+ animals, with other robust phenotypes affecting recognition memory and weight. To explore the dosage effect of the 16p11.2 genes on metabolism, Del/+ and Dup/+ models were challenged with high fat and high sugar diet, which revealed opposite energy imbalance. Transcriptomic analysis revealed that the majority of the genes located in the Sult1a1-Spn region were sensitive to dosage with a major effect on several pathways associated with neurocognitive and metabolic phenotypes. Whereas the behavioral consequence of the 16p11 region genetic dosage was similar in mice and humans with activity and memory alterations, the metabolic defects were opposite: adult Del/+ mice are lean in comparison to the human obese phenotype and the Dup/+ mice are overweight in comparison to the human underweight phenotype. Together, these data indicate that the dosage imbalance at the 16p11.2 locus perturbs the expression of modifiers outside the CNV that can modulate the penetrance, expressivity and direction of effects in both humans and mice.

Hierarchical analyses of genetic differentiation in a hybrid zone of Sorex araneus (Insectivora : Soricidae)

Microsatellites are used to unravel the fine-scale genetic structure of a hybrid zone between chromosome races Valais and Cordon of the common shrew (Sorex araneus) located in the French Alps. A total of 269 individuals collected between 1992 and 1995 was typed for seven microsatellite loci. A modified version of the classical multiple correspondence analysis is carried out. This analysis clearly shows the dichotomy between the two races. Several approaches are used to study genetic structuring. Gene flow is clearly reduced between these chromosome races and is estimated at one migrant every two generations using X-statistics and one migrant per generation using F-statistics. Hierarchical F- and R-statistics are compared and their efficiency to detect inter- and intraracial patterns of divergence is discussed. Within-race genetic structuring is significant, but remains weak. F-ST displays similar values on both sides of the hybrid zone, although no environmental barriers are found on the Cordon side, whereas the Valais side is divided by several mountain rivers. We introduce the exact G-test to microsatellite data which proved to be a powerful test to detect genetic differentiation within as well as among races. The genetic background of karyotypic hybrids was compared with the genetic background of pure parental forms using a CRT-MCA. Our results indicate that, without knowledge of the karyotypes, we would not have been able to distinguish these hybrids from karyotypically pure samples.

Clinal variation in a brown lemur (Eulemur spp.) hybrid zone: Combining morphological, genetic and climatic data to examine stability.

Studies of hybrid zones can inform our understanding of reproductive isolation and speciation. Two species of brown lemur (Eulemur rufifrons and E. cinereiceps) form an apparently stable hybrid zone in the Andringitra region of south-eastern Madagascar. The aim of this study was to identify factors that contribute to this stability. We sampled animals at 11 sites along a 90-km transect through the hybrid zone and examined variation in 26 microsatellites, the D-loop region of mitochondrial DNA, six pelage and nine morphological traits; we also included samples collected in more distant allopatric sites. Clines in these traits were noncoincident, and there was no increase in either inbreeding coefficients or linkage disequilibrium at the centre of the zone. These results could suggest that the hybrid zone is maintained by weak selection against hybrids, conforming to either the tension zone or geographical selection-gradient model. However, a closer examination of clines in pelage and microsatellites indicates that these clines are not sigmoid or stepped in shape but instead plateau at their centre. Sites within the hybrid zone also occur in a distinct habitat, characterized by greater seasonality in precipitation and lower seasonality in temperature. Together, these findings suggest that the hybrid zone may follow the bounded superiority model, with exogenous selection favouring hybrids within the transitional zone. These findings are noteworthy, as examples supporting the bounded superiority model are rare and may indicate a process of ecologically driven speciation without geographical isolation.

Genetic and karyotypic structure of a hybrid zone between the chromosomal races Cordon and Valais in the common shrew, Sorex araneus

Hybrid zones between genetically differentiated populations provide material to study evolutionary processes. Since the discovery of chromosomal races in Sorex araneus, contact zones have attracted attention of scientists. So far, studies on genetic markers in Sorex hybrid zones are missing. The acrocentric chromosomal race Cordon and the highly metacentric race Valais meet and hybridize at Les Houches in the Western Alps. On a transect through the hybrid zone, 273 shrews were caught at 15 localities over 4 years. Karyotype as well as the nuclear protein loci Alb and Pg were analyzed. F-st and F-is values were calculated by F-statistics. An analysis on pooled samples revealed the genetical differences between the hybridizing races as the only cause of population structuring. Genetical markers show dines with very strong frequency shifts at a mountain torrent, but behave differently through the zone. The performance of the torrent in maintaining the hybrid zone, selection against hybrids, possible assortative mating and linkage of the Valais Pg allele to a diagnostic chromosome arm, are discussed.

A vertebrate reproductive system involving three ploidy levels : hybrid origin of triploids in a contact zone of diploid and tetraploid green toads (Bufo viridis subgroup)

The rise and consequences of polyploidy in vertebrates, whose origin was associated with genome duplications, may be best studied in natural diploid and polyploid populations. In a diploid/tetraploid (2n/4n) geographic contact zone of Palearctic green toads in northern Kyrgyzstan, we examine 4ns and triploids (3n) of unknown genetic composition and origins. Using mitochondrial and nuclear sequence, and nuclear microsatellite markers in 84 individuals, we show that 4n (Bufo pewzowi) are allopolyploids, with a geographically proximate 2n species (B. turanensis) being their maternal ancestor and their paternal ancestor as yet unidentified. Local 3n forms arise through hybridization. Adult 3n mature males (B. turanensis mtDNA) have 2n mothers and 4n fathers, but seem distinguishable by nuclear profiles from partly aneuploid 3n tadpoles (with B. pewzowi mtDNA). These observations suggest multiple pathways to the formation of triploids in the contact zone, involving both reciprocal origins. To explain the phenomena in the system, we favor a hypothesis where 3n males (with B. turanensis mtDNA) backcross with 4n and 2n females. Together with previous studies of a separately evolved, sexually reproducing 3n lineage, these observations reveal complex reproductive interactions among toads of different ploidy levels and multiple pathways to the evolution of polyploid lineages.

Evolution, maintenance and ecological consequences of genetic caste determination in " Pogonomyrmex " harvester ants

Summary Division of labor between reproducers (queens) and helpers (workers) is the main characteristic of social insect societies and at the root of their ecological success. Kin selection models predict that phenotypic differences between queens and workers should result from environmental rather than from genetic differences. However, genetic effects on queen and worker differentiation were found in two populations-of Pogonomyrmex harvester ants. Each of the two populations is composed of two genetically distinct lineages. Queens (which can be of either lineage) generally mate with males of their own and of the alternate lineage and produce two types of female offspring, those fertilized by males of the queens' lineage which develop into queens and those fertilized by males of the alternate lineage which develop into workers. All four lineages were further suggested to be themselves of hybrid origin between-the species P: barbatus and P. rugosus, in which queens and workers do not differ genetically. In a first set of experiments, we tested if female caste determination (the differentiation into queens and workers) in the lineages was genetically hardwired and if it was associated with costs in terms of the ability to optimally allocate resources to the production of queens and workers. To this end we first mated queens of-two lineages to a single male. Queens mated to a male of the alternate lineage successfully raised worker offspring whereas queens mated to a male of their own lineage almost always failed to produce workers. This reveals that pure-lineage individuals have lost the ability to develop into workers. Second, we analyzed offspring produced by naturally mated queens. During the stage of colony founding when only workers are produced, naturally mated queens laid a high proportion of pure-lineage eggs but the large majority of these eggs failed to develop. As a consequence, the number of offspring produced by incipient colonies decreased linearly with the proportion of pure-lineage eggs laid by queens. Moreover, queens of the lineage most commonly represented in a given population produced more pure-lineage eggs, in line with the view that they mate randomly with the two types of males and indiscriminately use their sperm. Altogether these results predict frequency-dependent founding success for pairs of lineages because queens of the more common lineage will produce more pure-lineage eggs and their colonies be less successful during the stage of colony founding. To describe the distribution of populations characterized with genetic caste determination relative to the populations with environmental caste determination we genotyped queens and workers collected during a large survey of -additional populations. Genetic caste determination associated with pairs of interbreeding lineages was frequent and widespread in the studied range and we identified four additional lineages displaying genetic caste determination. Overall, there were thus eight highly differentiated lineages with genetic caste determination. These lineages always co-occurred in the same complementary lineage pairs. Three of the four lineage pairs appeared to have a common origin, while their relationship with the forth could not be resolved. The genetic survey also revealed that, in addition to being genetically isolated from one another, all eight lineages were genetically distinct from P. rugosus and P. barbatus, even when colonies of interbreeding lineages co-occurred with colonies of either putative parent at the same site. This raised the question of the mechanisms involved in the reproductive isolation between the lineages and the parental species and between the two lineages of a lineage pair. At a site where one lineage pair co-occurred with P. rugosus, we identified two pre-zygotic mechanisms (differences in timing for mating flights between P. rugosus and the lineage pair and assortative mating) and one post-zygotic mechanism (high levels of hybrid unviablility) which in combination may largely account for the reproductive isolation between the lineages and their parental species. The mechanisms accounting for the reproductive isolation between the two lineages of a lineage pair varied across lineage pairs. In one lineage pair, inter-lineage individuals exclusively occurred in the sterile worker caste, raising the possibility that inter-lineage eggs have completely lost the ability to develop into queens in this lineage pair and that there is thus no opportunity for gene flow. In each of the three remaining lineage pairs, inter-lineage queens were produced by a minority of colonies. In these lineage pairs, colonies headed by inter-lineage queens failed to grow sufficiently to produce reproductive individuals which may account for the reproductive isolation between co-occurring lineages in three lineage pairs. In conclusion, the results of this thesis show that genetic caste determination is costly but widespread in Pogonomyrmex harvester ants. Reproductive isolation among the lineages and between the lineages and the parental species as well as frequency-dependent founding success for co-occurring lineages may contribute to the persistence of this extraordinary system. Résumé La division du travail entre individus reproducteurs (les reines) et individus non-reproducteurs (ouvrières) représente la caractéristique principale des sociétés d'insectes et est à la base de leur succès écologique. Des modèles de sélection de parentèle prédisent que les différences phénotypiques entre reines et ouvrières devraient provenir d'effets environnementaux plutôt que de différences génétiques. Malgré ce fait, des effets génétiques sur la différentiation entre reines et ouvrières ont été montrés dans deux populations de fourmis moissonneuses du genre Pogonomyrmex. Chacune des deux populations est composée de deux lignées génétiquement distinctes. Les reines de chaque lignée s'accouplent en général avec des mâles de leur propre lignée ainsi qu'avec des mâles de l'autre lignée et produisent deux types d'oeufs, ceux qui sont fécondés par les mâles de leur propre lignée qui se développent en nouvelles reines et ceux qui sont fécondés par les mâles de l'autre lignée qui se développent en ouvrières. Il a été suggéré que les lignées sont elles-mêmes des hybrides entre les deux espèces P. barbatus et P. rugosus. Dans ces deux espèces, les reines et ouvrières ne sont pas génétiquement distinctes. Dans une première série d'expériences, nous avons testé si la détermination de la caste femelle (le développement en reine ou en ouvrière) est génétiquement rigide et si elle est associée à des coûts en terme de capacité à allouer de façon optimale les ressources pour la production de reines et d'ouvrières. Pour cela nous avons accouplé des reines de deux lignées avec un seul mâle. Les reines accouplées avec un mâle de l'autre lignée ont élevé de nouvelles ouvrières avec succès alors que les reines accouplées avec un mâle de leur propre lignée ont presque toujours échoué à produire des ouvrières. Ceci montre que les individus de lignée pure ont perdu la capacité de se développer en ouvrière. Deuxièmement, nous avons analysé la descendance de reines qui se sont accouplées naturellement. Durant le stade de fondation de la colonie, où seules des ouvrières sont élevées, les reines accouplées naturellement ont pondu une grande proportion d'oeufs de lignée pure mais la majorité de ces derniers ne se sont pas développés. En conséquence, le nombre de descendants produits par des colonies fondatrices diminuait linéairement avec la proportion des oeufs de lignée pure pondus par la reine en accord avec l'hypothèse que les reines s'accouplent au hasard avec les deux types de mâles et utilisent leur sperme aléatoirement. Dans l'ensemble; ces résultats prédisent un succès de fondation fréquence-dépendant pour les deux lignées, car les reines de la lignée la plus fréquente produiront .plus d'oeufs de lignée pure et leurs colonies auront moins de succès lors de la fondation de colonies par rapport aux colonies de la lignée la moins fréquente. Pour décrire la distribution des-populations caractérisées par une détermination génétique des castes par rapport aux populations caractérisées par une détermination environnementale des castes, nous avons génotypé des reines et des ouvrières qui ont été collectées lors d'une analyse de populations supplémentaires. La détermination génétique des castes associée à des croisements entre lignées est fréquente et largement répartie dans l'aire étudiée. Nous avons identifié quatre lignées supplémentaires, ayant une détermination génétique des castes, pour un total de huit lignées. Ces huit lignées forment quatre paires de lignées et on ne trouve jamais deux lignées de paires différentes, dans une population. Trois des quatre paires de lignées s'avèrent avoir une origine commune alors que leur relation avec la quatrième paire de lignées n'a pas pu être résolue. L'analyse génétique de populations supplémentaires a également révélé qu'en plus d'être génétiquement isolées les unes des autres, les huit lignées sont génétiquement distinctes de P. rugosus et P. barbatus même si les colonies d'une paire de lignées se trouvent en sympatrie avec l'une ou l'autre des espèces parentales. Ceci relève la question des mécanismes impliqués dans l'isolation reproductive entre les lignées et les espèces parentales ainsi qu'entre les deux lignées d'une paire. En étudiant un site où une paire de lignées se trouve en sympatrie avec P. rugosus, nous avons identifié deux mécanismes pré-zygotiques (des différences dans le timing du vol nuptial entre P. rugosus et les lignées et des accouplements assortis) ainsi qu'un mécanisme post-zygotique (un niveau élevé de non-viabilité des hybrides). En combinaison, ces mécanismes peuvent largement expliquer l'isolement reproductif entre les lignées et leurs espèces parentales. Les mécanismes contribuant à l'isolement reproductif entre les deux lignées d'une paire variaient entre paires de lignées. Dans une paire, les individus de génotype inter-lignée se trouvent uniquement dans la caste stérile des ouvrières, suggérant qu'il n'y a pas d'opportunité pour avoir du flux de gènes entre les deux lignées ce cette paire. Dans chacune des trois autres paires de lignées des nouvelles reines de génotype inter-lignée sont produites par une minorité de colonies. Par contre, les colonies avec une reine mère de génotype inter-lignée ne se développent pas suffisamment pour produire des individus reproducteurs. Ceci peut donc expliquer pourquoi il n'y a pas de flux de gènes entre les deux lignées de trois paires. En conclusion, les résultats de cette thèse montrent que la détermination génétique de la caste est coûteuse mais très répandue chez les fourmis. moissonneuses du genre Pogonomyrmex. L'isolement reproductif des lignées entre elles et avec les espèces parentales, ainsi qu'un succès de fondation fréquence-dépendant contribuent à la persistance de ce système extraordinaire.

Mechanisms of reproductive isolation between an ant species of hybrid origin and one of its parents.

The establishment of new species by hybridization is difficult because it requires the development of reproductive isolation (RI) in sympatry to escape the homogenizing effects of gene flow from the parental species. Here we investigated the role of two pre- and two postzygotic mechanisms of RI in a system comprising two interdependent Pogonomyrmex harvester ant lineages (the H1 and H2 lineages) of hybrid origin and one of their parental species (P. rugosus). Similar to most other ants, P. rugosus is characterized by an environmental system of caste determination with female brood developing either into queens or workers depending on nongenetic factors. By contrast, there is a strong genetic component to caste determination in the H1 and H2 lineages because the developmental fate of female brood depends on the genetic origin of the parents, with interlineage eggs developing into workers and intralineage eggs developing into queens. The study of a mixed mating aggregation revealed strong differences in mating flight timing between P. rugosus and the two lineages as a first mechanism of RI. A second important prezygotic mechanism was assortative mating. Laboratory experiments also provided support for one of the two investigated mechanisms of postzygotic isolation. The majority of offspring produced from the few matings between P. rugosus and the lineages aborted at the egg stage. This hybrid inviability was under maternal influence, with hybrids produced by P. rugosus queens being always inviable whereas a small proportion of H2 lineage queens produced large numbers of adult hybrid offspring. Finally, we found no evidence that genetic caste determination acted as a second postzygotic mechanism reducing gene flow between P. rugosus and the H lineages. The few viable P. rugosus-H hybrids were not preferentially shunted into functionally sterile workers but developed into both workers and queens. Overall, these results reveal that the nearly complete (99.5%) RI between P. rugosus and the two hybrid lineages stems from the combination of two typical prezygotic mechanisms (mating time divergence and assortative mating) and one postzygotic mechanism (hybrid inviability).

Hybrid Multiscale Finite Volume method for two-phase flow in porous media

We present a novel hybrid (or multiphysics) algorithm, which couples pore-scale and Darcy descriptions of two-phase flow in porous media. The flow at the pore-scale is described by the Navier?Stokes equations, and the Volume of Fluid (VOF) method is used to model the evolution of the fluid?fluid interface. An extension of the Multiscale Finite Volume (MsFV) method is employed to construct the Darcy-scale problem. First, a set of local interpolators for pressure and velocity is constructed by solving the Navier?Stokes equations; then, a coarse mass-conservation problem is constructed by averaging the pore-scale velocity over the cells of a coarse grid, which act as control volumes; finally, a conservative pore-scale velocity field is reconstructed and used to advect the fluid?fluid interface. The method relies on the localization assumptions used to compute the interpolators (which are quite straightforward extensions of the standard MsFV) and on the postulate that the coarse-scale fluxes are proportional to the coarse-pressure differences. By numerical simulations of two-phase problems, we demonstrate that these assumptions provide hybrid solutions that are in good agreement with reference pore-scale solutions and are able to model the transition from stable to unstable flow regimes. Our hybrid method can naturally take advantage of several adaptive strategies and allows considering pore-scale fluxes only in some regions, while Darcy fluxes are used in the rest of the domain. Moreover, since the method relies on the assumption that the relationship between coarse-scale fluxes and pressure differences is local, it can be used as a numerical tool to investigate the limits of validity of Darcy's law and to understand the link between pore-scale quantities and their corresponding Darcy-scale variables.

Exogenous selection shapes germination behaviour and seedling traits of populations at different altitudes in a Senecio hybrid zone.

BACKGROUND AND AIMS: The Senecio hybrid zone on Mt Etna, Sicily, is characterized by steep altitudinal clines in quantitative traits and genetic variation. Such clines are thought to be maintained by a combination of 'endogenous' selection arising from genetic incompatibilities and environment-dependent 'exogenous' selection leading to local adaptation. Here, the hypothesis was tested that local adaptation to the altitudinal temperature gradient contributes to maintaining divergence between the parental species, S. chrysanthemifolius and S. aethnensis. METHODS: Intra- and inter-population crosses were performed between five populations from across the hybrid zone and the germination and early seedling growth of the progeny were assessed. KEY RESULTS: Seedlings from higher-altitude populations germinated better under low temperatures (9-13 °C) than those from lower altitude populations. Seedlings from higher-altitude populations had lower survival rates under warm conditions (25/15 °C) than those from lower altitude populations, but also attained greater biomass. There was no altitudinal variation in growth or survival under cold conditions (15/5 °C). Population-level plasticity increased with altitude. Germination, growth and survival of natural hybrids and experimentally generated F(1)s generally exceeded the worse-performing parent. CONCLUSIONS: Limited evidence was found for endogenous selection against hybrids but relatively clear evidence was found for divergence in seed and seedling traits, which is probably adaptive. The combination of low-temperature germination and faster growth in warm conditions might enable high-altitude S. aethnensis to maximize its growth during a shorter growing season, while the slower growth of S. chrysanthemifolius may be an adaptation to drought stress at low altitudes. This study indicates that temperature gradients are likely to be an important environmental factor generating and maintaining adaptive divergence across the Senecio hybrid zone on Mt Etna.

Nosology and classification of genetic skeletal disorders: 2010 revision.

Genetic disorders involving the skeletal system arise through disturbances in the complex processes of skeletal development, growth and homeostasis and remain a diagnostic challenge because of their variety. The Nosology and Classification of Genetic Skeletal Disorders provides an overview of recognized diagnostic entities and groups them by clinical and radiographic features and molecular pathogenesis. The aim is to provide the Genetics, Pediatrics and Radiology community with a list of recognized genetic skeletal disorders that can be of help in the diagnosis of individual cases, in the delineation of novel disorders, and in building bridges between clinicians and scientists interested in skeletal biology. In the 2010 revision, 456 conditions were included and placed in 40 groups defined by molecular, biochemical, and/or radiographic criteria. Of these conditions, 316 were associated with mutations in one or more of 226 different genes, ranging from common, recurrent mutations to "private" found in single families or individuals. Thus, the Nosology is a hybrid between a list of clinically defined disorders, waiting for molecular clarification, and an annotated database documenting the phenotypic spectrum produced by mutations in a given gene. The Nosology should be useful for the diagnosis of patients with genetic skeletal diseases, particularly in view of the information flood expected with the novel sequencing technologies; in the delineation of clinical entities and novel disorders, by providing an overview of established nosologic entities; and for scientists looking for the clinical correlates of genes, proteins and pathways involved in skeletal biology. © 2011 Wiley-Liss, Inc.

Chromosomal rearrangements and gene flow over time in an inter-specific hybrid zone of the Sorex araneus group.

Most hybrid zones have existed for hundreds or thousands of years but have generally been observed for only a short time period. Studies extending over periods long enough to track evolutionary changes in the zones or assess the ultimate outcome of hybridization are scarce. Here, we describe the evolution over time of the level of genetic isolation between two karyotypically different species of shrews (Sorex araneus and Sorex antinorii) at a hybrid zone located in the Swiss Alps. We first evaluated hybrid zone movement by contrasting patterns of gene flow and changes in cline parameters (centre and width) using 24 microsatellite loci, between two periods separated by 10 years apart. Additionally, we tested the role of chromosomal rearrangements on gene flow by analysing microsatellite loci located on both rearranged and common chromosomes to both species. We did not detect any movement of the hybrid zone during the period analysed, suggesting that the zone is a typical tension zone. However, the gene flow was significantly lower among the rearranged than the common chromosomes for the second period, whereas the difference was only marginally significant for the first period. This further supports the role of chromosomal rearrangements on gene flow between these taxa.