Comparative modular analysis of gene expression in vertebrate development

The focus of my PhD research was the concept of modularity. In the last 15 years, modularity has become a classic term in different fields of biology. On the conceptual level, a module is a set of interacting elements that remain mostly independent from the elements outside of the module. I used modular analysis techniques to study gene expression evolution in vertebrates. In particular, I identified ``natural'' modules of gene expression in mouse and human, and I showed that expression of organ-specific and system-specific genes tends to be conserved between such distance vertebrates as mammals and fishes. Also with a modular approach, I studied patterns of developmental constraints on transcriptome evolution. I showed that none of the two commonly accepted models of the evolution of embryonic development (``evo-devo'') are exclusively valid. In particular, I found that the conservation of the sequences of regulatory regions is highest during mid-development of zebrafish, and thus it supports the ``hourglass model''. In contrast, events of gene duplication and new gene introduction are most rare in early development, which supports the ``early conservation model''. In addition to the biological insights on transcriptome evolution, I have also discussed in detail the advantages of modular approaches in large-scale data analysis. Moreover, I re-analyzed several studies (published in high-ranking journals), and showed that their conclusions do not hold out under a detailed analysis. This demonstrates that complex analysis of high-throughput data requires a co-operation between biologists, bioinformaticians, and statisticians.

Characterization of the effect of mammary gland irridiation on experimental brest cancer progression and analysis of the implicated mechanisms

SUMMARY Radiotherapy is commonly and efficiently used to treat solid cancer in the clinic. Experimental evidence however suggests that radiation can promote tumor progression by inducing chronic modifications of the tumor microenvironment. Clinically, these observations are highly relevant to aggressive tumoral lesions relapsing after radiation therapy, a leading cause of patients' death. The investigation and understanding of the biological mechanisms implicated in the malignant progression of post-radiation relapses are therefore of major importance. Here we used a syngeneic (immunocompetent) breast cancer orthotopic xenograft model, to show that local irradiation of the mammary gland promotes the appearance of an invasive and metastatic tumor phenotype. Previous studies in our laboratory revealed that inhibition of tumor-induced angiogenesis and consequent increase in tumor hypoxia promotes metastasis formation through the activation of pro-invasive programs in the tumor cells. Our results extend these observations suggesting that mammary gland irradiation induces the recruitment of CD11b+ cells to both the primary tumor and the lungs at pre-metastatic stages through the hypoxia-dependent induction of Kit-ligand (KITL) expression in primary tumors. Abrogation of KITL expression in tumor cells prevented CD11 b+ cells accumulation in both the primary tumor and lungs and significantly reduced metastases of tumors growing in irradiated mammary gland. Importantly, irradiated mammary gland enhanced tumor-induced mobilization of circulating CD11b+cKit+ myelomonocytic cells through a HIF1- and KITL-dependent process. By cell transfer experiments, mobilized circulating CD11b+cKit+ cells were shown to supply both tumor- and lungs infiltrating CD11b+ cells. Using a blocking antibody against cKit (the KITL receptor), the mobilization of CD11b+cKit+ ceils was prevented as well as lung metastases derived from tumors growing in irradiated mammary gland. Taken together, these results indicate that tumors growing in a pre-irradiated mammary gland partially promote their malignant progression through the distant mobilization of circulating myelomonocytic precursor cells. They identify KITL inhibition and/or cKit receptor neutralization as potentially promising therapeutic approaches for post-radiation relapses. RESUME La radiothérapie est largement utilisée comme traitement de choix de nombreux types de cancers. L'agressivité des récidives tumorales observée en clinique après radiothérapie suggère cependant que le recours à l'irradiation pourrait dans certains cas accélérer la progression tumorale. De récents travaux expérimentaux ont en effet permis d'appuyer cette hypothèse, en montrant notamment l'effet néfaste des modifications chroniques de l'environnement induites par l'irradiation sur la progression tumorale. A l'aide d'un modèle murin syngénique orthotopique de cancer de sein, nous avons pu montrer que l'irradiation locale de la glande mammaire facilite l'invasion et la dissémination métastatique des cellules tumorales en favorisant le recrutement de cellules myéloïdes CD11 b+ vers la tumeur primaire et les poumons à un stade pré-métastatique. Comme mécanisme impliqué dans le recrutement des cellules CD11b+, nous avons pu observer après irradiation locale de la glande mammaire une expression augmentée de Kit-ligand (KITL) dans la tumeur (induite par l'hypoxie) ainsi que la mobilisation de cellules myéloïdes circulantes exprimant le récepteur cKit et précurseurs des cellules CD11b+ infiltrant la tumeur et les poumons. En empêchant la mobilisation par la tumeur de cellules circulantes cKit+ par des approches à la fois génétique et pharmacologique nous avons pu prévenir l'accumulation de cellules myéloïdes CD11 b+ dans la tumeur primaire et les poumons ainsi que la dissémination métastatique induites par' l'irradiation de la glande mammaire. De façon générale, ces résultats montrent que la progression agressive des tumeurs qui se développent dans un environnement irradié repose à la fois sur l'expression tumorale de KITL et la mobilisation de cellules myéloïdes précurseurs cKit*. Ils auront permis d'identifier KITL et/ou cKit comme des cibles thérapeutiques potentielles intéressantes pour le traitement des récidives tumorales après radiothérapie.

2011 lAAF World Championships in Daegu: blood tests for all athletes in the framework of the Athlete Biological Passport.

The 2011 International Association of Athletics Federation (IAAF) World Championships took place in Daegu, Korea. For the first time, all athletes were blood tested prior to the competition in order to give a clear signal to the world athletic community of the wish to enter into the era of the Athlete Biological Passport and fight against doping in their sport. The hematological parameters were measured on site. Thus, a mobile-accredited laboratory for blood testing was created in Daegu. Two serum tubes were collected for clinical chemistry and hormonal analyses in order to build the bases of the endocrine and the androgen (steroid) modules of the Athlete Biological Passport in blood. This paper describes some of the main challenges the project faced with regard to the large number of athletes, competing in different disciplines, and the logistic problems that had to be solved for smart implementation of one of the most complex operations organized in the last decade in the fight against doping.

Fate of Hydrophilic Nanoparticles in Biological Environments

The nanoparticles developed are based on chitosan, a biocompatible and biodegradable polysaccharide. The chitosan nanoparticles are formed in an entirely water-based process by electrostatic interactions with other biocompatible molecules. As a prerequisite to understand the fate of such nanoparticles in cells, comprehensive characterization and stability studies serve to identify quantitatively the impact of the raw material characteristics and preparation conditions on the nanoparticle characteristics. Methods included H-1 NMR spectroscopy, dilution viscometry, particle size analysis and electron microscopy. Cytotoxicity and cell uptake experiments on RAW 264.7 murine macrophages and p23 murine endothelial cells were performed to investigate the correlation with nanoparticle characteristics and effect of surface decoration with alginate. Cytotoxicity was assessed by the MTT survival test; cell uptake was monitored by fluorescent microscopy using labeled polymers.

An evolutionary analysis of the relationship between spite and altruism.

We investigate the selective pressures on a social trait when evolution occurs in a population of constant size. We show that any social trait that is spiteful simultaneously qualifies as altruistic. In other words, any trait that reduces the fitness of less related individuals necessarily increases that of related ones. Our analysis demonstrates that the distinction between "Hamiltonian spite" and "Wilsonian spite" is not justified on the basis of fitness effects. We illustrate this general result with an explicit model for the evolution of a social act that reduces the recipient's survival ("harming trait"). This model shows that the evolution of harming is favoured if local demes are of small size and migration is low (philopatry). Further, deme size and migration rate determine whether harming evolves as a selfish strategy by increasing the fitness of the actor, or as a spiteful/altruistic strategy through its positive effect on the fitness of close kin.

Mining the human phenome using allelic scores that index biological intermediates.

It is common practice in genome-wide association studies (GWAS) to focus on the relationship between disease risk and genetic variants one marker at a time. When relevant genes are identified it is often possible to implicate biological intermediates and pathways likely to be involved in disease aetiology. However, single genetic variants typically explain small amounts of disease risk. Our idea is to construct allelic scores that explain greater proportions of the variance in biological intermediates, and subsequently use these scores to data mine GWAS. To investigate the approach's properties, we indexed three biological intermediates where the results of large GWAS meta-analyses were available: body mass index, C-reactive protein and low density lipoprotein levels. We generated allelic scores in the Avon Longitudinal Study of Parents and Children, and in publicly available data from the first Wellcome Trust Case Control Consortium. We compared the explanatory ability of allelic scores in terms of their capacity to proxy for the intermediate of interest, and the extent to which they associated with disease. We found that allelic scores derived from known variants and allelic scores derived from hundreds of thousands of genetic markers explained significant portions of the variance in biological intermediates of interest, and many of these scores showed expected correlations with disease. Genome-wide allelic scores however tended to lack specificity suggesting that they should be used with caution and perhaps only to proxy biological intermediates for which there are no known individual variants. Power calculations confirm the feasibility of extending our strategy to the analysis of tens of thousands of molecular phenotypes in large genome-wide meta-analyses. We conclude that our method represents a simple way in which potentially tens of thousands of molecular phenotypes could be screened for causal relationships with disease without having to expensively measure these variables in individual disease collections.

Ultrastructural analysis of neuronal plasticity in the adult mouse

ABSTRACT Adult neuronal plasticity is a term that corresponds to a set of biological mechanisms allowing a neuronal circuit to respond and adapt to modifications of the received inputs. Mystacial whiskers of the mouse are the starting point of a major sensory pathway that provides the animal with information from its immediate environment. Through whisking, information is gathered that allows the animal to orientate itself and to recognize objects. This sensory system is crucial for nocturnal behaviour during which vision is not of much use. Sensory information of the whiskers are sent via brainstem and thalamus to the primary somatosensory area (S1) of the cerebral cortex in a strictly topological manner. Cell bodies in the layer N of S 1 are arranged in ring forming structures called barrels. As such, each barrel corresponds to the cortical representation in layer IV of a single whisker follicle. This histological feature allows to identify with uttermost precision the part of the cortex devoted to a given whisker and to study modifications induced by different experimental conditions. The condition used in the studies of my thesis is the passive stimulation of one whisker in the adult mouse for a period of 24 hours. It is performed by glueing a piece of metal on one whisker and placing the awake animal in a cage surrounded by an electromagnetic coil that generates magnetic field burst inducing whisker movement at a given frequency during 24 hours. I analysed the ultrastructure of the barrel corresponding the stimulated whisker using serial sections electron microscopy and computer-based three-dimensional reconstructions; analysis of neighbouring, unstimulated barrels as well as those from unstimulated mice served as control. The following elements were structurally analyzed: the spiny dendrites, the axons of excitatory as well as inhibitory cells, their connections via synapses and the astrocytic processes. The density of synapses and spines is upregulated in a barrel corresponding to a stimulated whisker. This upregulation is absent in the BDNF heterozygote mice, indicating that a certain level of activity-dependent released BDNF is required for synaptogenesis in the adult cerebral cortex. Synpaptogenesis is correlated with a modification of the astrocytes that place themselves in closer vicinity of the excitatory synapses on spines. Biochemical analysis revealed that the astrocytes upregulate the expression of transporters by which they internalise glutamate, the neurotransmitter responsible for the excitatory response of cortical neurons. In the final part of my thesis, I show that synaptogenesis in the stimulated barrel is due to the increase in the size of excitatory axonal boutons that become more frequently multisynaptic, whereas the inhibitory axons do not change their morphology but form more synapses with spines apposed to them. Taken together, my thesis demonstrates that all the cellular elements present in the neuronal tissue of the adult brain contribute to activity-dependent cortical plasticity and form part of a mechanism by which the animal responds to a modified sensory experience. Throughout life, the neuronal circuit keeps the faculty to adapt its function. These adaptations are partially transitory but some aspects remain and could be the structural basis of a memory trace in the cortical circuit. RESUME La plasticité neuronale chez l'adulte désigne un ensemble de mécanismes biologiques qui permettent aux circuits neuronaux de répondre et de s'adapter aux modifications des stimulations reçues. Les vibrisses des souris sont un système crucial fournissant des informations sensorielles au sujet de l'environnement de l'animal. L'information sensorielle collectée par les vibrisses est envoyée via le tronc cérébral et le thalamus à l'aire sensorielle primaire (S 1) du cortex cérébral en respectant strictement la somatotopie. Les corps cellulaires dans la couche IV de S 1 sont organisés en anneaux délimitant des structures nommées tonneaux. Chaque tonneau reçoit l'information d'une seule vibrisse et l'arrangement des tonneaux dans le cortex correspond à l'arrangement des vibrisses sur le museau de la souris. Cette particularité histologique permet de sélectionner avec certitude la partie du cortex dévolue à une vibrisse et de l'étudier dans diverses conditions. Le paradigme expérimental utilisé dans cette thèse est la stimulation passive d'une seule vibrisse durant 24 heures. Pour ce faire, un petit morceau de métal est collé sur une vibrisse et la souris est placée dans une cage entourée d'une bobine électromagnétique générant un champ qui fait vibrer le morceau de métal durant 24 heures. Nous analysons l'ultrastructure du cortex cérébral à l'aide de la microscopie électronique et des coupes sériées permettant la reconstruction tridimensionnelle à l'aide de logiciels informatiques. Nous observons les modifications des structures présentes : les dendrites épineuses, les axones des cellules excitatrices et inhibitrices, leurs connections par des synapses et les astrocytes. Le nombre de synapses et d'épines est augmenté dans un tonneau correspondant à une vibrisse stimulée 24 heures. Basé sur cela, nous montrons dans ces travaux que cette réponse n'est pas observée dans des souris hétérozygotes BDNF+/-. Cette neurotrophine sécrétée en fonction de l'activité neuronale est donc nécessaire pour la synaptogenèse. La synaptogenèse est accompagnée d'une modification des astrocytes qui se rapprochent des synapses excitatrices au niveau des épines dendritiques. Ils expriment également plus de transporteurs chargés d'internaliser le glutamate, le neurotransmetteur responsable de la réponse excitatrice des neurones. Nous montrons aussi que les axones excitateurs deviennent plus larges et forment plus de boutons multi-synaptiques à la suite de la stimulation tandis que les axones inhibiteurs ne changent pas de morphologie mais forment plus de synapses avec des épines apposées à leur membrane. Tous les éléments analysés dans le cerveau adulte ont maintenu la capacité de réagir aux modifications de l'activité neuronale et répondent aux modifications de l'activité permettant une constante adaptation à de nouveaux environnements durant la vie. Les circuits neuronaux gardent la capacité de créer de nouvelles synapses. Ces adaptations peuvent être des réponses transitoires aux stimuli mais peuvent aussi laisser une trace mnésique dans les circuits.

Functional analysis of a missense mutation in the serine protease inhibitor SPINT2 associated with congenital sodium diarrhea.

Membrane-bound serine proteases play important roles in different biological processes. Their regulation by endogenous inhibitors is poorly understood. A Y163C mutation in the SPINT2 gene encoding the serine protease inhibitor Hepatocyte Growth Factor Inhibitor HAI-2 is associated with a congenital sodium diarrhea. The functional consequences of this mutation on HAI-2 activity and its physiological targets are unknown. We established a cellular assay in Xenopus laevis oocytes to study functional interactions between HAI-2 and candidate membrane-bound serine proteases expressed in the gastro-intestinal tract. We found that the wild-type form of HAI-2 is a potent inhibitor of nine gastro-intestinal serine proteases. The Y163C mutation in the second Kunitz domain of HAI-2 resulted in a complete loss of inhibitory activity on two intestinal proteases, prostasin and tmprss13. The effect of the mutation of the homologous Y68C in the first Kunitz domain of HAI-2 is consistent with a differential contribution of the two Kunitz domains of HAI-2 in the inhibition of serine proteases. By contrast to the Tyr to Cys, the Tyr to Ser substitution did not change the inhibitory potency of HAI-2, indicating that the thiol-group of the cysteine rather than the Tyr deletion is responsible for the HAI-2 loss of function. Our functional assay allowed us to identify membrane-bound serine proteases as cellular target for inhibition by HAI-2 wild type and mutants, and to better define the role of the Tyr in the second Kunitz domain in the inhibitory activity of HAI-2.

Fc block treatment, dead cells exclusion, and cell aggregates discrimination concur to prevent phenotypical artifacts in the analysis of subpopulations of tumor-infiltrating CD11b(+) myelomonocytic cells.

It is well established that cancer cells can recruit CD11b(+) myeloid cells to promote tumor angiogenesis and tumor growth. Increasing interest has emerged on the identification of subpopulations of tumor-infiltrating CD11b(+) myeloid cells using flow cytometry techniques. In the literature, however, discrepancies exist on the phenotype of these cells (Coffelt et al., Am J Pathol 2010;176:1564-1576). Since flow cytometry analysis requires particular precautions for accurate sample preparation and trustable data acquisition, analysis, and interpretation, some discrepancies might be due to technical reasons rather than biological grounds. We used the syngenic orthotopic 4T1 mammary tumor model in immunocompetent BALB/c mice to analyze and compare the phenotype of CD11b(+) myeloid cells isolated from peripheral blood and from tumors, using six-color flow cytometry. We report here that the nonspecific antibody binding through Fc receptors, the presence of dead cells and cell doublets in tumor-derived samples concur to generate artifacts in the phenotype of tumor-infiltrating CD11b(+) subpopulations. We show that the heterogeneity of tumor-infiltrating CD11b(+) subpopulations analyzed without particular precautions was greatly reduced upon Fc block treatment, dead cells, and cell doublets exclusion. Phenotyping of tumor-infiltrating CD11b(+) cells was particularly sensitive to these parameters compared to circulating CD11b(+) cells. Taken together, our results identify Fc block treatment, dead cells, and cell doublets exclusion as simple but crucial steps for the proper analysis of tumor-infiltrating CD11b(+) cell populations.

A quantitative analysis of L-glutamate-regulated Na+ dynamics in mouse cortical astrocytes: implications for cellular bioenergetics.

The mode of Na+ entry and the dynamics of intracellular Na+ concentration ([Na+]i) changes consecutive to the application of the neurotransmitter glutamate were investigated in mouse cortical astrocytes in primary culture by video fluorescence microscopy. An elevation of [Na+]i was evoked by glutamate, whose amplitude and initial rate were concentration dependent. The glutamate-evoked Na+ increase was primarily due to Na+-glutamate cotransport, as inhibition of non-NMDA ionotropic receptors by 6-cyano-7-nitroquinoxiline-2,3-dione (CNQX) only weakly diminished the response and D-aspartate, a substrate of the glutamate transporter, produced [Na+]i elevations similar to those evoked by glutamate. Non-NMDA receptor activation could nevertheless be demonstrated by preventing receptor desensitization using cyclothiazide. Thus, in normal conditions non-NMDA receptors do not contribute significantly to the glutamate-evoked Na+ response. The rate of Na+ influx decreased during glutamate application, with kinetics that correlate well with the increase in [Na+]i and which depend on the extracellular concentration of glutamate. A tight coupling between Na+ entry and Na+/K+ ATPase activity was revealed by the massive [Na+]i increase evoked by glutamate when pump activity was inhibited by ouabain. During prolonged glutamate application, [Na+]i remains elevated at a new steady-state where Na+ influx through the transporter matches Na+ extrusion through the Na+/K+ ATPase. A mathematical model of the dynamics of [Na+]i homeostasis is presented which precisely defines the critical role of Na+ influx kinetics in the establishment of the elevated steady state and its consequences on the cellular bioenergetics. Indeed, extracellular glutamate concentrations of 10 microM already markedly increase the energetic demands of the astrocytes.

Quantitative morphometry analysis of the fetal brain using clinical MR imaging

In vivo fetal magnetic resonance imaging provides aunique approach for the study of early human braindevelopment [1]. In utero cerebral morphometry couldpotentially be used as a marker of the cerebralmaturation and help to distinguish between normal andabnormal development in ambiguous situations. However,this quantitative approach is a major challenge becauseof the movement of the fetus inside the amniotic cavity,the poor spatial resolution provided by very fast MRIsequences and the partial volume effect. Extensiveefforts are made to deal with the reconstruction ofhigh-resolution 3D fetal volumes based on severalacquisitions with lower resolution [2,3,4]. Frameworkswere developed for the segmentation of specific regionsof the fetal brain such as posterior fossa, brainstem orgerminal matrix [5,6], or for the entire brain tissue[7,8], applying the Expectation-Maximization MarkovRandom Field (EM-MRF) framework. However, many of theseprevious works focused on the young fetus (i.e. before 24weeks) and use anatomical atlas priors to segment thedifferent tissue or regions. As most of the gyraldevelopment takes place after the 24th week, acomprehensive and clinically meaningful study of thefetal brain should not dismiss the third trimester ofgestation. To cope with the rapidly changing appearanceof the developing brain, some authors proposed a dynamicatlas [8]. To our opinion, this approach however faces arisk of circularity: each brain will be analyzed /deformed using the template of its biological age,potentially biasing the effective developmental delay.Here, we expand our previous work [9] to proposepost-processing pipeline without prior that allow acomprehensive set of morphometric measurement devoted toclinical application. Data set & Methods: Prenatal MRimaging was performed with a 1-T system (GE MedicalSystems, Milwaukee) using single shot fast spin echo(ssFSE) sequences (TR 7000 ms, TE 180 ms, FOV 40 x 40 cm,slice thickness 5.4mm, in plane spatial resolution1.09mm). For each fetus, 6 axial volumes shifted by 1 mmwere acquired under motherâeuro?s sedation (about 1min pervolume). First, each volume is segmentedsemi-automatically using region-growing algorithms toextract fetal brain from surrounding maternal tissues.Inhomogeneity intensity correction [10] and linearintensity normalization are then performed. Brain tissues(CSF, GM and WM) are then segmented based on thelow-resolution volumes as presented in [9]. Ahigh-resolution image with isotropic voxel size of 1.09mm is created as proposed in [2] and using B-splines forthe scattered data interpolation [11]. Basal gangliasegmentation is performed using a levet setimplementation on the high-resolution volume [12]. Theresulting white matter image is then binarized and givenas an input in FreeSurfer software(http://surfer.nmr.mgh.harvard.edu) to providetopologically accurate three-dimensional reconstructionsof the fetal brain according to the local intensitygradient. References: [1] Guibaud, Prenatal Diagnosis29(4) (2009). [2] Rousseau, Acad. Rad. 13(9), 2006. [3]Jiang, IEEE TMI 2007. [4] Warfield IADB, MICCAI 2009. [5]Claude, IEEE Trans. Bio. Eng. 51(4) 2004. [6] Habas,MICCAI 2008. [7] Bertelsen, ISMRM 2009. [8] Habas,Neuroimage 53(2) 2010. [9] Bach Cuadra, IADB, MICCAI2009. [10] Styner, IEEE TMI 19(39 (2000). [11] Lee, IEEETrans. Visual. And Comp. Graph. 3(3), 1997. [12] BachCuadra, ISMRM 2010.

Molecular analysis of a novel gene cluster encoding an insect toxin in plant-associated strains of Pseudomonas fluorescens

Pseudomonas fluorescens CHA0 and the related strain Pf-5 are well-characterized representatives of rhizosphere bacteria that have the capacity to protect crop plants from fungal root diseases, mainly by releasing a variety of exoproducts that are toxic to plant pathogenic fungi. Here, we report that the two plant-beneficial pseudomonads also exhibit potent insecticidal activity. Anti-insect activity is linked to a novel genomic locus encoding a large protein toxin termed Fit (for P. fluorescensinsecticidal toxin) that is related to the insect toxin Mcf (Makes caterpillars floppy) of the entomopathogen Photorhabdus luminescens, a mutualist of insect-invading nematodes. When injected into the haemocoel, even low doses of P. fluorescens CHA0 or Pf-5 killed larvae of the tobacco hornworm Manduca sexta and the greater wax moth Galleria mellonella. In contrast, mutants of CHA0 or Pf-5 with deletions in the Fit toxin gene were significantly less virulent to the larvae. When expressed from an inducible promoter in a non-toxic Escherichia coli host, the Fit toxin gene was sufficient to render the bacterium toxic to both insect hosts. Our findings establish the Fit gene products of P. fluorescens CHA0 and Pf-5 as potent insect toxins that define previously unappreciated anti-insect properties of these plant-colonizing bacteria

Problems of allometric scaling analysis: examples from mammalian reproductive biology.

Biological scaling analyses employing the widely used bivariate allometric model are beset by at least four interacting problems: (1) choice of an appropriate best-fit line with due attention to the influence of outliers; (2) objective recognition of divergent subsets in the data (allometric grades); (3) potential restrictions on statistical independence resulting from phylogenetic inertia; and (4) the need for extreme caution in inferring causation from correlation. A new non-parametric line-fitting technique has been developed that eliminates requirements for normality of distribution, greatly reduces the influence of outliers and permits objective recognition of grade shifts in substantial datasets. This technique is applied in scaling analyses of mammalian gestation periods and of neonatal body mass in primates. These analyses feed into a re-examination, conducted with partial correlation analysis, of the maternal energy hypothesis relating to mammalian brain evolution, which suggests links between body size and brain size in neonates and adults, gestation period and basal metabolic rate. Much has been made of the potential problem of phylogenetic inertia as a confounding factor in scaling analyses. However, this problem may be less severe than suspected earlier because nested analyses of variance conducted on residual variation (rather than on raw values) reveals that there is considerable variance at low taxonomic levels. In fact, limited divergence in body size between closely related species is one of the prime examples of phylogenetic inertia. One common approach to eliminating perceived problems of phylogenetic inertia in allometric analyses has been calculation of 'independent contrast values'. It is demonstrated that the reasoning behind this approach is flawed in several ways. Calculation of contrast values for closely related species of similar body size is, in fact, highly questionable, particularly when there are major deviations from the best-fit line for the scaling relationship under scrutiny.