Analysis of gene expression patterns in animals

During my PhD, my aim was to provide new tools to increase our capacity to analyse gene expression patterns, and to study on a large-scale basis the evolution of gene expression in animals. Gene expression patterns (when and where a gene is expressed) are a key feature in understanding gene function, notably in development. It appears clear now that the evolution of developmental processes and of phenotypes is shaped both by evolution at the coding sequence level, and at the gene expression level.Studying gene expression evolution in animals, with complex expression patterns over tissues and developmental time, is still challenging. No tools are available to routinely compare expression patterns between different species, with precision, and on a large-scale basis. Studies on gene expression evolution are therefore performed only on small genes datasets, or using imprecise descriptions of expression patterns.The aim of my PhD was thus to develop and use novel bioinformatics resources, to study the evolution of gene expression. To this end, I developed the database Bgee (Base for Gene Expression Evolution). The approach of Bgee is to transform heterogeneous expression data (ESTs, microarrays, and in-situ hybridizations) into present/absent calls, and to annotate them to standard representations of anatomy and development of different species (anatomical ontologies). An extensive mapping between anatomies of species is then developed based on hypothesis of homology. These precise annotations to anatomies, and this extensive mapping between species, are the major assets of Bgee, and have required the involvement of many co-workers over the years. My main personal contribution is the development and the management of both the Bgee database and the web-application.Bgee is now on its ninth release, and includes an important gene expression dataset for 5 species (human, mouse, drosophila, zebrafish, Xenopus), with the most data from mouse, human and zebrafish. Using these three species, I have conducted an analysis of gene expression evolution after duplication in vertebrates.Gene duplication is thought to be a major source of novelty in evolution, and to participate to speciation. It has been suggested that the evolution of gene expression patterns might participate in the retention of duplicate genes. I performed a large-scale comparison of expression patterns of hundreds of duplicated genes to their singleton ortholog in an outgroup, including both small and large-scale duplicates, in three vertebrate species (human, mouse and zebrafish), and using highly accurate descriptions of expression patterns. My results showed unexpectedly high rates of de novo acquisition of expression domains after duplication (neofunctionalization), at least as high or higher than rates of partitioning of expression domains (subfunctionalization). I found differences in the evolution of expression of small- and large-scale duplicates, with small-scale duplicates more prone to neofunctionalization. Duplicates with neofunctionalization seemed to evolve under more relaxed selective pressure on the coding sequence. Finally, even with abundant and precise expression data, the majority fate I recovered was neither neo- nor subfunctionalization of expression domains, suggesting a major role for other mechanisms in duplicate gene retention.

Free-living energy expenditure measured by two independent techniques in pregnant and nonpregnant Gambian women.

Free-living energy expenditure (EE) was assessed in 37 young pregnant Gambian women at the 12th (n = 11, 53.5 +/- 1.7 kg), 24th (n = 14, 54.7 +/- 2.1 kg), and 36th (n = 12, 65.0 +/- 2.6 kg) wk of pregnancy and was compared with nonpregnant nonlactating (NPNL) control women (n = 12, 50.3 +/- 1.6 kg). The following two methods were used to assess EE: 1) the heart rate (HR) method using individual regression lines (HR vs EE) established at different activity levels in a respiration chamber and 2) the doubly labeled water (2H2(18)O) method in a subgroup of 25 pregnant and 7 control women. With the HR method the EE during the agricultural rainy season was found to be 2,408 +/- 87, 2,293 +/- 122, and 2,782 +/- 130 kcal/day at 12, 24, and 36 wk of gestation and were not significantly different from the control group (2,502 +/- 133 kcal/day). These findings were confirmed by the 2H2(18)O measurements, which failed to show any effect of pregnancy on EE. Expressed per unit body weight, the free-living EE was found to be lower (P less than 0.01 with 2H2(18)O method) at 36 wk of gestation than in the NPNL group. It is concluded that, in these Gambian women, energy-sparing mechanisms that contribute to meet the additional energy stress of gestation are operating during pregnancy (e.g., diminished spontaneous physical activity).

Analysis of the mechanisms controlling the activity of flp1, a chromosomal passenger protein in the fission Schizosaccharomyces pombe

ABSTRACT In S. cerevisiae, the protein phosphatase Cdc14pwt is essential far mitotic exit through its contribution to reducing mitotic CDK activity. But Cdc14pwt also acts as a mare general temporal coordinator of mid and late mitotic events by controlling the partitioning of DNA, microtubule stability and cytokinesis. Cdc14pwt orthologs are well conserved from yeasts to humans, and sequence comparison revealed the presence of three domains, A, B and C, of which A and B form the catalytic domain. Cdc14pwt orthologs are regulated (in part) through cell cycle dependent changes in their localization. Some of them are thought to be kept inactive by sequestration in the nucleolus during interphase. This is the case for flp1pwt, the single identified Cdc14pwt ortholog in the fission yeast S. pombe. In early mitosis, flp1pwt leaves the nucleolus and localizes to the kinetochores, the contractile ring and the mitotic spindle, suggesting that it has multiple substrates and regulates many mitotic processes. flp1D cells show a high chromosome loss rate and septation defects, suggesting a role for flp1wt in the fidelity of chromosome transmission and cytokinesis. The aim of this study is to characterize the mechanisms underlying flp1pwt functions and the control of its activity. A structure-function analysis has revealed that the presence of both A and B domains is required for biological function and for proper flp1pwt mitotic localization. In contrast, the C domain of flp1pwt is responsible for its proper nucleolar localization in G2/interphase. My data suggest that dephosphorylation of substrates by flp1pwt is not necessary for any changes in localization of flp1pwt except that at the medial ring. In that particular case, the catalytic activity of flp1pwt is required for efficient localization, therefore revealing an additional level of regulation. All the functions of flp1pwt assayed to date require its catalytic activity, emphasizing the importance of further identification of its substrates. As described for other orthologs, the capability of selfinteraction and phosphorylation status might help to control flp1pwt activity. My data suggest that flp1pwt forms oligomers in vivo and that phosphorylation is not essential far localization changes of the protein. In addition, the hypophosphorylated form of flp1pwt might be specifically involved in the promotion of cytokinesis. The results of this study suggest that multiple modes of regulation including localization, selfassociation and phosphorylation allow a fine-tuning regulation of flp1pwt phosphatase activity, and more generally that of Cdc14pwt family of phosphatases. RESUME Chez la levure S. cerevisiae, la protéine phosphatase Cdc14pwt est essentielle pour la sortie de mitose du fait de sa contribution dans la réduction d'activité des CDK mitotiques. Comme elle contrôle également le partage de l'ADN, la stabilité des microtubules et la cytokinèse, Cdc14pwt est en fait considérée comme un coordinateur temporel général des évènements de milieu et de fin de mitose. Les orthologues de Cdc14pwt sont bien conservés, des levures jusqu'à l'espèce humaine. Des comparaisons de séquence ont révélé la présence de trois domaines A, B et C, les deux premiers constituant le domaine catalytique. Ils sont régulés (en partie) via des changements dans leur localisation, eux-mêmes dépendants du cycle cellulaire. Plusieurs de ces orthologues sont supposés inactivés par séquestration dans le nucléole en interphase, ce qui est le cas de flp1pwt le seul orthologue de Cdc14pwt identifié chez la levure fissipare S, pombe. En début de mitose, flp1pwt quitte le nucléole et localise au niveau des kinetochores, de l'anneau contractile d'actine et du fuseau mitotique, ce qui laisse supposer de multiples substrats et fonctions. Comme les cellules délétées pour le gène flp1wt présentent un taux élevé de perte de chromosome et des défauts de septation, flp1pwt semble jouer un rôle dans la fidélité de la transmission du matériel génétique et la cytokinèse. Le but de cette étude est de caractériser les mécanismes impliqués dans les fonctions assurées par flp1pwt d'une part, et dans le contrôle de son activité d'autre part. Une analyse structure-fonction a révélé que la présence simultanée des deux domaines A et B est requise pour la fonction biologique de flp1pwt et sa localisation correcte pendant la mitose. Par contre, le domaine C de flp1pwt confère une localisation nucléolaire adéquate en G2/interphase. Mes données suggèrent que la déphosphorylation de substrats par flp1pwt est dispensable pour sa localisation correcte excepté celle à l'anneau médian, qui requiert dans ce cas, l'activité catalytique de flp1pwt, révélant ainsi un niveau de régulation supplémentaire. Toutes les fonctions de flp1 pwt testées jusqu'à présent nécessitent également son activité catalytique, ce qui accentue l'importance de l'identification future de ses substrats. Comme cela a déjà été décrit pour d'autres orthologues, la capacité d'auto-intéraction et le niveau de phosphorylation pourraient contrôler l'activité de flp1pwt. En effet, mes données suggèrent que flp1pwt forme des oligomères in vivo et que la phosphorylation n'est pas essentielle pour les changements de localisation observés pour la protéine. De plus, la forme hypophosphorylée de flp1pwt pourrait être spécifiquement impliquée dans la promotion de la cytokinèse. De multiples modes de régulation incluant la localisation, l'auto-association et la phosphorylation semblent permettre un contrôle fin et subtil de l'activité de la phosphatase flp1pwt, et plus généralement celle des protéines de la famille de Cdc14pwt.

Trace element partitioning in HP-LT metamorphic assemblages during subduction-related metamorphism, Ile de Groix, France: a detailed LA-ICPMS study

Devolatilization reactions and subsequent transfer of fluid from subducted oceanic crust into the overlying mantle wedge are important processes, which are responsible for the specific geochemical characteristics of subduction-related metamorphic rocks, as well as those of arc magmatism. To better understand the geochemical fingerprint induced by fluid mobilization during dehydration and rehydration processes related to subduction zone metamorphism, the trace element and rare earth element (REE) distribution patterns in HP-LT metamorphic assemblages in eclogite-, blueschist- and greenschist-facies rocks of the Ile de Groix were obtained by laser ablation inductively coupled plasma mass spectrometry (LA-ICPMS) analysis. This study focuses on 10 massive basic rocks representing former hydrothermally altered mid-ocean ridge basalts (MORB), four banded basic rocks of volcano-sedimentary origin and one micaschist. The main hosts for incompatible trace elements are epidote (REE, Th, U, Pb, Sr), garnet [Y, heavy REE (HREE)], phengite (Cs, Rb, Ba, B), titanite [Ti, Nb, Ta, REE; HREE > LREE (light REE)], rutile (Ti, Nb, Ta) and apatite (REE, Sr). The trace element contents of omphacite, amphibole, albite and chlorite are low. The incompatible trace element contents of minerals are controlled by the stable metamorphic mineral assemblage and directly related to the appearance, disappearance and reappearance of minerals, especially epidote, garnet, titanite, rutile and phengite, during subduction zone metamorphism. Epidote is a key mineral in the trace element exchange process because of its large stability field, ranging from lower greenschist- to blueschist- and eclogite-facies conditions. Different generations of epidote are generally observed and related to the coexisting phases at different stages of the metamorphic cycle (e.g. lawsonite, garnet, titanite). Epidote thus controls most of the REE budget during the changing P-T conditions along the prograde and retrograde path. Phengite also plays an important role in determining the large ion lithophile element (LILE) budget, as it is stable to high P-T conditions. The breakdown of phengite causes the release of LILE during retrogression. A comparison of trace element abundances in whole-rocks and minerals shows that the HP-LT metamorphic rocks largely retain the geochemical characteristics of their basic, volcano-sedimentary and pelitic protoliths, including a hydrothermal alteration overprint before the subduction process. A large part of the incompatible trace elements remained trapped in the rocks and was recycled within the various metamorphic assemblages stable under changing metamorphic conditions during the subduction process, indicating that devolatilization reactions in massive basic rocks do not necessarily imply significant simultaneous trace element and REE release.

Prediction of resting energy expenditure from fat-free mass and fat mass.

On the basis of literature values, the relationship between fat-free mass (FFM), fat mass (FM), and resting energy expenditure [REE (kJ/24 h)] was determined for 213 adults (86 males, 127 females). The objectives were to develop a mathematical model to predict REE based on body composition and to evaluate the contribution of FFM and FM to REE. The following regression equations were derived: 1) REE = 1265 + (93.3 x FFM) (r2 = 0.727, P < 0.001); 2) REE = 1114 + (90.4 x FFM) + (13.2 x FM) (R2 = 0.743, P < 0.001); and 3) REE = (108 x FFM) + (16.9 x FM) (R2 = 0.986, P < 0.001). FM explained only a small part of the variation remaining after FFM was accounted for. The models that include both FFM and FM are useful in examination of the changes in REE that occur with a change in both the FFM and FM. To account for more of the variability in REE, FFM will have to be divided into organ mass and skeletal muscle mass in future analyses.

The use of uniaxial accelerometry for the assessment of physical-activity-related energy expenditure: a validation study against whole-body indirect calorimetry.

Assessing the total energy expenditure (TEE) and the levels of physical activity in free-living conditions with non-invasive techniques remains a challenge. The purpose of the present study was to investigate the accuracy of a new uniaxial accelerometer for assessing TEE and physical-activity-related energy expenditure (PAEE) over a 24 h period in a respiratory chamber, and to establish activity levels based on the accelerometry ranges corresponding to the operationally defined metabolic equivalent (MET) categories. In study 1, measurement of the 24 h energy expenditure of seventy-nine Japanese subjects (40 (SD 12) years old) was performed in a large respiratory chamber. During the measurements, the subjects wore a uniaxial accelerometer (Lifecorder; Suzuken Co. Ltd, Nagoya, Japan) on their belt. Two moderate walking exercises of 30 min each were performed on a horizontal treadmill. In study 2, ten male subjects walked at six different speeds and ran at three different speeds on a treadmill for 4 min, with the same accelerometer. O2 consumption was measured during the last minute of each stage and was expressed in MET. The measured TEE was 8447 (SD 1337) kJ/d. The accelerometer significantly underestimated TEE and PAEE (91.9 (SD 5.4) and 92.7 (SD 17.8) % chamber value respectively); however, there was a significant correlation between the two values (r 0.928 and 0.564 respectively; P<0.001). There was a strong correlation between the activity levels and the measured MET while walking (r(2) 0.93; P<0.001). Although TEE and PAEE were systematically underestimated during the 24 h period, the accelerometer assessed energy expenditure well during both the exercise period and the non-structured activities. Individual calibration factors may help to improve the accuracy of TEE estimation, but the average calibration factor for the group is probably sufficient for epidemiological research. This method is also important for assessing the diurnal profile of physical activity.

Complete nucleotide sequence of Birmingham IncP alpha plasmids. Compilation and comparative analysis.

The IncP alpha promiscuous plasmid (R18, R68, RK2, RP1 and RP4) comprises 60,099 bp of nucleotide sequence, encoding at least 74 genes. About 40 kb of the genome, designated the IncP core and including all essential replication and transfer functions, can be aligned with equivalent sequences in the IncP beta plasmid R751. The compiled IncP alpha sequence revealed several previously unidentified reading frames that are potential genes. IncP alpha plasmids carry genetic information very efficiently: the coding sequences of the genes are closely packed but rarely overlap, and occupy almost 86% of the genome's nucleotide sequence. All of the 74 genes should be expressed, although there is as yet experimental evidence for expression of only 60 of them. Six examples of tandem-in-frame initiation sites specifying two gene products each are known. Two overlapping gene arrangements occupy different reading frames of the same region. Intergenic regions include most of the 25 promoters; transcripts are usually polycistronic. Translation of most of the open reading frames seems to be initiated independently, each from its own ribosomal binding and initiation site, although, a few cases of coupled translation have been reported. The most frequently used initiation codon is AUG but translation for a few open reading frames begins at GUG or UUG. The most common stop-codon is UGA followed by UAA and then UAG. Regulatory circuits are complex and largely dependent on two components of the central control operon. KorA and KorB are transcriptional repressors controlling at least seven operons. KorA and KorB act synergistically in several cases by recognizing and binding to conserved nucleotide sequences. Twelve KorB binding sites were found around the IncP alpha sequence and these are conserved in R751 (IncP beta) with respect to both sequence and location. Replication of IncP alpha plasmids requires oriV and the plasmid-encoded initiator protein TrfA in combination with the host-encoded replication machinery. Conjugative plasmid transfer depends on two separate regions occupying about half of the genome. The primary segregational stability system designated Par/Mrs consists of a putative site-specific recombinase, a possible partitioning apparatus and a post-segregational lethality mechanism, all encoded in two divergent operons. Proteins related to the products of F sop and P1 par partitioning genes are separately encoded in the central control operon.

Discordances in the application of different criteria for quantification of paediatric abdominal obesity: an analysis of two Swiss studies.

Several definitions of paediatric abdominal obesity have been proposed but it is unclear whether they lead to similar results. We assessed the prevalence of abdominal obesity using five different waist circumference-based definitions and their agreement with total body fat (TBF) and abdominal fat (AF). Data from 190 girls and 162 boys (Ballabeina), and from 134 girls and 113 boys (Kinder-Sportstudie, KISS) aged 5-11 years were used. TBF was assessed by bioimpedance (Ballabeina) or dual energy X-ray absorption (KISS). On the basis of the definition used, the prevalence of abdominal obesity varied between 3.1 and 49.4% in boys, and 4.7 and 55.5% in girls (Ballabeina), and between 1.8 and 36.3% in boys and 4.5 and 37.3% in girls (KISS). Among children considered as abdominally obese by at least one definition, 32.0 (Ballabeina) and 44.7% (KISS) were considered as such by at least two (out of five possible) definitions. Using excess TBF or AF as reference, the areas under the receiver operating curve varied between 0.577 and 0.762 (Ballabeina), and 0.583 and 0.818 (KISS). We conclude that current definitions of abdominal obesity in children lead to wide prevalence estimates and should not be used until a standard definition can be proposed.

Use of phycoerythrin and allophycocyanin for fluorescence resonance energy transfer analyzed by flow cytometry: advantages and limitations.

BACKGROUND: This study validates the use of phycoerythrin (PE) and allophycocyanin (APC) for fluorescence energy transfer (FRET) analyzed by flow cytometry. METHODS: FRET was detected when a pair of antibody conjugates directed against two noncompetitive epitopes on the same CD8alpha chain was used. FRET was also detected between antibody conjugate pairs specific for the two chains of the heterodimeric alpha (4)beta(1) integrin. Similarly, the association of T-cell receptor (TCR) with a soluble antigen ligand was detected by FRET when anti-TCR antibody and MHC class I/peptide complexes (&lt;&lt;tetramers&gt;&gt;) were used. RESULTS: FRET efficiency was always less than 10%, probably because of steric effects associated with the size and structure of PE and APC. Some suggestions are given to take into account this and other effects (e.g., donor and acceptor concentrations) for a better interpretation of FRET results obtained with this pair of fluorochromes. CONCLUSIONS: We conclude that FRET assays can be carried out easily with commercially available antibodies and flow cytometers to study arrays of multimolecular complexes.

Energy gap in the aetiology of body weight gain and obesity: a challenging concept with a complex evaluation and pitfalls.

The concept of energy gap(s) is useful for understanding the consequence of a small daily, weekly, or monthly positive energy balance and the inconspicuous shift in weight gain ultimately leading to overweight and obesity. Energy gap is a dynamic concept: an initial positive energy gap incurred via an increase in energy intake (or a decrease in physical activity) is not constant, may fade out with time if the initial conditions are maintained, and depends on the 'efficiency' with which the readjustment of the energy imbalance gap occurs with time. The metabolic response to an energy imbalance gap and the magnitude of the energy gap(s) can be estimated by at least two methods, i.e. i) assessment by longitudinal overfeeding studies, imposing (by design) an initial positive energy imbalance gap; ii) retrospective assessment based on epidemiological surveys, whereby the accumulated endogenous energy storage per unit of time is calculated from the change in body weight and body composition. In order to illustrate the difficulty of accurately assessing an energy gap we have used, as an illustrative example, a recent epidemiological study which tracked changes in total energy intake (estimated by gross food availability) and body weight over 3 decades in the US, combined with total energy expenditure prediction from body weight using doubly labelled water data. At the population level, the study attempted to assess the cause of the energy gap purported to be entirely due to increased food intake. Based on an estimate of change in energy intake judged to be more reliable (i.e. in the same study population) and together with calculations of simple energetic indices, our analysis suggests that conclusions about the fundamental causes of obesity development in a population (excess intake vs. low physical activity or both) is clouded by a high level of uncertainty.

Major mite allergen Der f 1 concentration is reduced in buildings with improved energy performance.

BACKGROUND: Environmental conditions play a crucial role in mite growth, and optimal environmental control is key in the prevention of airway inflammation in chronic allergic rhinoconjunctivitis or asthma. OBJECTIVE: To evaluate the relationship between building energy performance and indoor mite allergen concentration in a cross-sectional study. METHODS: Major allergen concentration (Der f 1, Der p 1, mite group 2, Fel d 1 and Bla g 2) was determined by quantitative dot blot analysis from mattress and carpet dust samples in five buildings designed for low energy use (LEB) and in six control buildings (CB). Inhabitants had received 4 weeks prior to mite measurement a personal validated questionnaire related to the perceived state of health and comfort of living. RESULTS: Cumulative mite allergen concentration (with Der f 1 as the major contributor) was significantly lower in LEB as compared with CB both in mattresses and in carpets. In contrast, the two categories of buildings did not differ in Bla g 2 and Fel d 1 concentration, in the amount of dust and airborne mould collected. Whereas temperature was higher in LEB, relative humidity was significantly lower than in CB. Perceived overall comfort was better in LEB. CONCLUSIONS: Major mite allergen Der f 1 preferentially accumulates in buildings not specifically designed for low energy use, reaching levels at risk for sensitization. We hypothesize that controlled mechanical ventilation present in all audited LEB may favour lower air humidity and hence lower mite growth and allergen concentration, while preserving optimal perceived comfort.

Energy deficit and length of hospital stay can be reduced by a two-step quality improvement of nutrition therapy: the intensive care unit dietitian can make the difference.

OBJECTIVE: Critically ill patients are at high risk of malnutrition. Insufficient nutritional support still remains a widespread problem despite guidelines. The aim of this study was to measure the clinical impact of a two-step interdisciplinary quality nutrition program. DESIGN: Prospective interventional study over three periods (A, baseline; B and C, intervention periods). SETTING: Mixed intensive care unit within a university hospital. PATIENTS: Five hundred seventy-two patients (age 59 ± 17 yrs) requiring >72 hrs of intensive care unit treatment. INTERVENTION: Two-step quality program: 1) bottom-up implementation of feeding guideline; and 2) additional presence of an intensive care unit dietitian. The nutrition protocol was based on the European guidelines. MEASUREMENTS AND MAIN RESULTS: Anthropometric data, intensive care unit severity scores, energy delivery, and cumulated energy balance (daily, day 7, and discharge), feeding route (enteral, parenteral, combined, none-oral), length of intensive care unit and hospital stay, and mortality were collected. Altogether 5800 intensive care unit days were analyzed. Patients in period A were healthier with lower Simplified Acute Physiologic Scale and proportion of "rapidly fatal" McCabe scores. Energy delivery and balance increased gradually: impact was particularly marked on cumulated energy deficit on day 7 which improved from -5870 kcal to -3950 kcal (p < .001). Feeding technique changed significantly with progressive increase of days with nutrition therapy (A: 59% days, B: 69%, C: 71%, p < .001), use of enteral nutrition increased from A to B (stable in C), and days on combined and parenteral nutrition increased progressively. Oral energy intakes were low (mean: 385 kcal*day, 6 kcal*kg*day ). Hospital mortality increased with severity of condition in periods B and C. CONCLUSION: A bottom-up protocol improved nutritional support. The presence of the intensive care unit dietitian provided significant additional progression, which were related to early introduction and route of feeding, and which achieved overall better early energy balance.

Multi-scale systems analysis of plant development : beyond cellular signaling and tissue morphodynamics

Les plantes sont essentielles pour les sociétés humaines. Notre alimentation quotidienne, les matériaux de constructions et les sources énergétiques dérivent de la biomasse végétale. En revanche, la compréhension des multiples aspects développementaux des plantes est encore peu exploitée et représente un sujet de recherche majeur pour la science. L'émergence des technologies à haut débit pour le séquençage de génome à grande échelle ou l'imagerie de haute résolution permet à présent de produire des quantités énormes d'information. L'analyse informatique est une façon d'intégrer ces données et de réduire la complexité apparente vers une échelle d'abstraction appropriée, dont la finalité est de fournir des perspectives de recherches ciblées. Ceci représente la raison première de cette thèse. En d'autres termes, nous appliquons des méthodes descriptives et prédictives combinées à des simulations numériques afin d'apporter des solutions originales à des problèmes relatifs à la morphogénèse à l'échelle de la cellule et de l'organe. Nous nous sommes fixés parmi les objectifs principaux de cette thèse d'élucider de quelle manière l'interaction croisée des phytohormones auxine et brassinosteroïdes (BRs) détermine la croissance de la cellule dans la racine du méristème apical d'Arabidopsis thaliana, l'organisme modèle de référence pour les études moléculaires en plantes. Pour reconstruire le réseau de signalement cellulaire, nous avons extrait de la littérature les informations pertinentes concernant les relations entre les protéines impliquées dans la transduction des signaux hormonaux. Le réseau a ensuite été modélisé en utilisant un formalisme logique et qualitatif pour pallier l'absence de données quantitatives. Tout d'abord, Les résultats ont permis de confirmer que l'auxine et les BRs agissent en synergie pour contrôler la croissance de la cellule, puis, d'expliquer des observations phénotypiques paradoxales et au final, de mettre à jour une interaction clef entre deux protéines dans la maintenance du méristème de la racine. Une étude ultérieure chez la plante modèle Brachypodium dystachion (Brachypo- dium) a révélé l'ajustement du réseau d'interaction croisée entre auxine et éthylène par rapport à Arabidopsis. Chez ce dernier, interférer avec la biosynthèse de l'auxine mène à la formation d'une racine courte. Néanmoins, nous avons isolé chez Brachypodium un mutant hypomorphique dans la biosynthèse de l'auxine qui affiche une racine plus longue. Nous avons alors conduit une analyse morphométrique qui a confirmé que des cellules plus anisotropique (plus fines et longues) sont à l'origine de ce phénotype racinaire. Des analyses plus approfondies ont démontré que la différence phénotypique entre Brachypodium et Arabidopsis s'explique par une inversion de la fonction régulatrice dans la relation entre le réseau de signalisation par l'éthylène et la biosynthèse de l'auxine. L'analyse morphométrique utilisée dans l'étude précédente exploite le pipeline de traitement d'image de notre méthode d'histologie quantitative. Pendant la croissance secondaire, la symétrie bilatérale de l'hypocotyle est remplacée par une symétrie radiale et une organisation concentrique des tissus constitutifs. Ces tissus sont initialement composés d'une douzaine de cellules mais peuvent aisément atteindre des dizaines de milliers dans les derniers stades du développement. Cette échelle dépasse largement le seuil d'investigation par les moyens dits 'traditionnels' comme l'imagerie directe de tissus en profondeur. L'étude de ce système pendant cette phase de développement ne peut se faire qu'en réalisant des coupes fines de l'organe, ce qui empêche une compréhension des phénomènes cellulaires dynamiques sous-jacents. Nous y avons remédié en proposant une stratégie originale nommée, histologie quantitative. De fait, nous avons extrait l'information contenue dans des images de très haute résolution de sections transverses d'hypocotyles en utilisant un pipeline d'analyse et de segmentation d'image à grande échelle. Nous l'avons ensuite combiné avec un algorithme de reconnaissance automatique des cellules. Cet outil nous a permis de réaliser une description quantitative de la progression de la croissance secondaire révélant des schémas développementales non-apparents avec une inspection visuelle classique. La formation de pôle de phloèmes en structure répétée et espacée entre eux d'une longueur constante illustre les bénéfices de notre approche. Par ailleurs, l'exploitation approfondie de ces résultats a montré un changement de croissance anisotropique des cellules du cambium et du phloème qui semble en phase avec l'expansion du xylème. Combinant des outils génétiques et de la modélisation biomécanique, nous avons démontré que seule la croissance plus rapide des tissus internes peut produire une réorientation de l'axe de croissance anisotropique des tissus périphériques. Cette prédiction a été confirmée par le calcul du ratio des taux de croissance du xylème et du phloème au cours de développement secondaire ; des ratios élevés sont effectivement observés et concomitant à l'établissement progressif et tangentiel du cambium. Ces résultats suggèrent un mécanisme d'auto-organisation établi par un gradient de division méristématique qui génèrent une distribution de contraintes mécaniques. Ceci réoriente la croissance anisotropique des tissus périphériques pour supporter la croissance secondaire. - Plants are essential for human society, because our daily food, construction materials and sustainable energy are derived from plant biomass. Yet, despite this importance, the multiple developmental aspects of plants are still poorly understood and represent a major challenge for science. With the emergence of high throughput devices for genome sequencing and high-resolution imaging, data has never been so easy to collect, generating huge amounts of information. Computational analysis is one way to integrate those data and to decrease the apparent complexity towards an appropriate scale of abstraction with the aim to eventually provide new answers and direct further research perspectives. This is the motivation behind this thesis work, i.e. the application of descriptive and predictive analytics combined with computational modeling to answer problems that revolve around morphogenesis at the subcellular and organ scale. One of the goals of this thesis is to elucidate how the auxin-brassinosteroid phytohormone interaction determines the cell growth in the root apical meristem of Arabidopsis thaliana (Arabidopsis), the plant model of reference for molecular studies. The pertinent information about signaling protein relationships was obtained through the literature to reconstruct the entire hormonal crosstalk. Due to a lack of quantitative information, we employed a qualitative modeling formalism. This work permitted to confirm the synergistic effect of the hormonal crosstalk on cell elongation, to explain some of our paradoxical mutant phenotypes and to predict a novel interaction between the BREVIS RADIX (BRX) protein and the transcription factor MONOPTEROS (MP),which turned out to be critical for the maintenance of the root meristem. On the same subcellular scale, another study in the monocot model Brachypodium dystachion (Brachypodium) revealed an alternative wiring of auxin-ethylene crosstalk as compared to Arabidopsis. In the latter, increasing interference with auxin biosynthesis results in progressively shorter roots. By contrast, a hypomorphic Brachypodium mutant isolated in this study in an enzyme of the auxin biosynthesis pathway displayed a dramatically longer seminal root. Our morphometric analysis confirmed that more anisotropic cells (thinner and longer) are principally responsible for the mutant root phenotype. Further characterization pointed towards an inverted regulatory logic in the relation between ethylene signaling and auxin biosynthesis in Brachypodium as compared to Arabidopsis, which explains the phenotypic discrepancy. Finally, the morphometric analysis of hypocotyl secondary growth that we applied in this study was performed with the image-processing pipeline of our quantitative histology method. During its secondary growth, the hypocotyl reorganizes its primary bilateral symmetry to a radial symmetry of highly specialized tissues comprising several thousand cells, starting with a few dozens. However, such a scale only permits observations in thin cross-sections, severely hampering a comprehensive analysis of the morphodynamics involved. Our quantitative histology strategy overcomes this limitation. We acquired hypocotyl cross-sections from tiled high-resolution images and extracted their information content using custom high-throughput image processing and segmentation. Coupled with an automated cell type recognition algorithm, it allows precise quantitative characterization of vascular development and reveals developmental patterns that were not evident from visual inspection, for example the steady interspace distance of the phloem poles. Further analyses indicated a change in growth anisotropy of cambial and phloem cells, which appeared in phase with the expansion of xylem. Combining genetic tools and computational modeling, we showed that the reorientation of growth anisotropy axis of peripheral tissue layers only occurs when the growth rate of central tissue is higher than the peripheral one. This was confirmed by the calculation of the ratio of the growth rate xylem to phloem throughout secondary growth. High ratios are indeed observed and concomitant with the homogenization of cambium anisotropy. These results suggest a self-organization mechanism, promoted by a gradient of division in the cambium that generates a pattern of mechanical constraints. This, in turn, reorients the growth anisotropy of peripheral tissues to sustain the secondary growth.

Trophic resource partitioning and competition between the two sibling bat species Myotis myotis and Myotis blythii

1. Niche theory predicts that the stable coexistence of species within a guild should be associated, if resources are limited, with a mechanism of resource partitioning. Using extensive data on diets, the present study attempts: (i) to test the hypothesis that, in sympatry, the interspecific overlap between the trophic niches of the sibling bat species Myotis myotis and M. blythii-which coexist intimately in their roosts-is effectively lower than the two intraspecific overlaps; (ii) to assess the role played by interspecific competition in resource partitioning through the study of trophic niche displacement between several sympatric and allopatric populations. 2. Diets were determined by the analysis of faecal samples collected in the field from individual bats captured in various geographical areas. Trophic niche overlaps were calculated monthly for all possible intraspecific and interspecific pairs of individuals from sympatric populations. Niche breadth was estimated from: (i) every faecal sample; (ii) all the faecal samples collected per month in a given population (geographical area). 3. In every population, the bulk of the diets of M. myotis and M. blythii consisted of, respectively, terrestrial (e.g. carabid beetles) and grass-dwelling (mostly bush crickets) prey. All intraspecific trophic niche overlaps were significantly greater than the interspecific one, except in Switzerland in May when both species exploited mass concentrations of cockchafers, a non-limiting food source. This clearcut partitioning of resources may allow the stable, intimate coexistence observed under sympatric conditions. 4. Relative proportions of ground-and grass-dwelling prey, as well as niche breadths (either individual or population), did not differ significantly between sympatry and allopatry, showing that, under allopatric conditions, niche expansion does not take place. This suggests that active interspecific competition is not the underlying mechanism responsible for the niche partitioning which is currently observed between M. myotis and M. blythii.