Excision of the Staphylococcal Cassette Chromosome mec (SCCmec) and its crucial role in the horizontal transfer of methicillin resistance in staphylococci.

Staphylococcus aureus est un pathogène humain majeur ayant développé des résistances contre la quasi totalité des antibiotiques disponibles, incluant la très importante famille des β- lactamines. La résistance à cette classe d'antibiotiques est conférée par la « Staphylococcal Cassette Chromosome mec » (SCCmec), qui est un élément génétique mobile capable de s'insérer dans le chromosome bactérien et capable d'être transféré horizontalement chez d'autres staphylocoques. Le mécanisme moléculaire impliqué dans ce transfert horizontal demeure largement inconnu. L'une des premières étapes du transfert est l'excision du SCC mec du chromosome bactérien. Cette excision est promue par des enzymes codées par l'élément SCCmec lui- même et appelées de ce fait « Cassette Chromosome Recombinases » (Ccr). L'un des buts de ce travail de thèse a été de comprendre la régulation de l'expression des gènes codant pour les Ccr recombinases. En utilisant des outils moléculaires originaux, nous avons été en mesure de démontrer en premier lieu que les Ccr recombinases étaient exprimées de façon « bistable », c'est à dire qu'uniquement quelques pourcents de cellules dans une population exprimaient ces gènes à un temps donné. Dans un deuxième temps, nous avons également démontré que l'expression de ces gènes était régulée par des facteurs étrangers au SCC mec. L'expression bistable des recombinases est un concept important. Effectivement, cela permet à la majorité des cellules d'une population de conserver l'élément SCC mec, alors que seulement une petite fraction le perd afin de le rendre disponible pour un transfert. Ainsi, alors que l'élément SCC mec continue de se propager avec la multiplication des bactéries Staphylococcus aureus résistant à la méticilline (SARM), il peut être simultanément transmis à des souches susceptibles (Staphylococcus aureus susceptible à la méticilline, SASM), entraînant l'apparition de nouveaux SARM. De façon très intéressante, le fait que cette bistabilité est contrôlée par les bactéries, et non le SCCmec lui-même, montre que la décision de transférer ou non la cassette SCC mec appartient à la bactérie. En conséquence, il doit exister dans la nature des souches qui sont plus ou moins aptes à effectuer ce transfert. En nous appuyant sur ces observations, nous avons montré que l'excision du SCC mec était effectivement régulée de façon très étroite au cours de la division cellulaire, et ne se passait que pendant un temps limité au début de la croissance. Ce résultat est compatible avec une régulation génétique commandée par la densité cellulaire, qui pourrait être dépendante de la production de signaux extracellulaires, du type que l'on rencontre dans le quorum sensing. Les signaux hypothétiques entraînant l'excision du SCC mec restent inconnus à l'heure actuelle. La connaissance de ces signaux pourrait se révéler très importante afin de développer des stratégies pour interférer avec la dissémination de la résistance au β-lactamines. Deux sujets additionnels ont été logiquement investigués au vu de ces premiers résultats. Premièrement, si certaines souches de SARM sont plus ou moins aptes à déclencher l'excision du SCC mec, de même certaines souches de SASM devraient être plus ou moins aptes à acquérir cet élément. Deuxièmement, afin d'étudier ces mécanismes de transfert au niveau épidémiologique, il nous a été nécessaire de développer des outils nous permettant d'explorer le phénomène à une plus large échelle. Concernant le premier point, il a été postulé que certains SASM seraient réfractaires à l'intégration génomique d'un SCC mec en raison de polymorphismes particuliers à proximité du site d'insertion chromosomique (attB). En étudiant plus de 40 isolais de S. aureus, provenant de porteurs sains, nous avons confirmé ce polymorphisme dans l'environnement à'attB. De plus, nous avons pu montrer que ces régions polymorphiques ont évolué parallèlement à des groupes phylogénétiques bien connus. Ainsi, si des telles régions réfractaires à l'intégration de SCC mec existent, celles-ci devraient ségréger dans des complexes clonaux bien définis qui devraient être facilement identifiables au niveau épidémiologique. Concernant le second point, nous avons été capables de construire un système rapporteur de l'excision du SCCmec, en utilisant un plasmide à faible copie. Ce système consistait en un promoteur fort et un gène codant pour une protéine verte fluorescente (GFP) sous le contrôle d'un promoteur fort séparés à l'aide d'un élément SCC artificiel portant trois terminateurs de transcription. Ainsi, la fluorescence ne s'exprime que si l'élément SCC est excisé du plasmide. Ce système a été testé avec succès dans plusieurs types de staphylocoques, et est actuellement évalué dans d'autres souches et conditions stimulant ou inhibant l'excision. De manière générale, cette dissertation représente parcours scientifique à travers plusieurs aspects d'un problème de santé publique majeur en rapport avec la résistance bactérienne aux antibiotiques. Ce travail s'attaque à des problèmes fondamentaux concernant le transfert horizontal de l'élément SCC mec. De plus, il s'intéresse à des aspects plus généraux de cet élément génétique mobile qui pourraient se révéler très importants en terme de mouvement de gènes au sein des staphylocoques, voir d'autres bactéries gram-positives. Finalement ce travail de thèse met en place le fondamentaux requis pour des recherches futures visant à interférer avec le transfert horizontal de la résistance aux β-lactamines. - Staphylococcus aureus is a major human pathogen. Moreover, S. aureus have developed resistance to almost all available antibiotics, including the important family of β-lactam molecules. Intrinsic resistance to β-lactams is conferred by the Staphylococcal Cassette Chromosome mec (SCCmec), which is a mobile genomic island that inserts into the staphylococcal chromosome and can be horizontally transferred into other staphylococci. However, little is known about the molecular mechanisms involved in this horizontal transfer into naïve strains. One of the first steps in SCC mec horizontal transfer is its excision from the chromosome. Excision is mediated by recombinase enzymes that are encoded by SCC mec itself, and named accordingly Ccr recombinases - for Cassette Chromosome recombinases. One goal of this thesis was to understand the regulation these recombinase genes. By using original molecular tools we could demonstrate first that the Ccr recombinases were expressed in a "bistable" manner, i.e. in only few percentages of the bacterial cells at a given time, and second that they were regulated by determinants that were not encoded on the SCC mec element, but elsewhere on the staphylococcal genome. "Bistable" expression Ccr recombinases is an important concept. It allows SCC mec to be excised and thus available for horizontal transfer, while ensuring that only some cells, but not the whole population, loose their valuable SCC mec genes. Thus, while the SCC mec element expands with the multiplication of the MRSA colony, it can simultaneously be transmitted into methicillin-susceptible S. aureus (MSSA), which convert into new MRSA. Most interestingly, the fact that bistability was regulated by the cells, rather than by SCC mec, indicates that it was the choice of the bacteria to trigger or not SCC mec transfer. As a consequence, there must be, in nature, staphylococcal strains that are more or less prone to sustain SCC mec transfer. Following these seminal observations we found that excision was indeed tightly regulated during bacterial division, and occurred only during a limited period of time at the beginning of bacterial growth. This is compatible with cell-density mediated gene regulation, and may depend on the production of extracellular signal molecules that transmit appropriate orders to neighboring cells, such as in quorum sensing. The potential signal triggering SCCmec excision is as yet unknown. However, it could be critical in promoting the horizontal transfer of methicillin resistance, or for the possible development of means to interfere with it. Two additional hypothesis were logically investigated in the view of these first results. First, if some strains of MRSA might be more prone than others to promote SCC mec excision, then some strains of MS SA might be more or less prone to acquire the element as well. Second, to investigate these multiple mechanisms at an epidemiological level, one would need to develop tools amenable to explore S. aureus strains at a larger scale. Regarding the first issue, it was postulated by others that some MSSA might be refractory to SCC mec integration because they had peculiar DNA polymorphisms in the vicinity of the site-specific chromosomal entry point {attB) of SCC mec. By studying >40 S. aureus isolates from healthy carriers, we confirmed the polymorphism of the attB environment. Moreover, we could show that these polymorphic regions co-evolved with well-known phylogenic clonal clusters. Therefore, if SCCwec-refractory attB environments exist, then they would segregate in well- defined S. aureus clonal clusters that would be easy to identify at the epidemiological level. Regarding the second issue, we were able to construct a new excision reporter system in a low copy number S. aureus plasmid. The reporter system consists in a strong promoter driving a green fluorescent protein {gfp) gene, separated by an artificial SCC-like element carrying three transcriptional terminators. Thus, fluorescence is not expressed unless the SCC-like element is excised. The system has been successfully tested in several aureus and non- aureus staphylococci, and is now being applied to more strains and various excision- triggering or inhibiting conditions. Altogether the dissertation is a scientific journey through various aspects of a salient medical problem with regard to antibiotic resistance and public health threat. The research work tackles fundamental issues about the mechanisms of horizontal transfer of the SCC mec element. Moreover, it also addresses more general features of this mobile element, which could be of larger importance with regard to gene trafficking in staphylococci, and maybe other gram-positive bacteria. Finally, the dissertation sets the fundamentals for future work and possible new ways to interfere with the horizontal transfer of methicillin resistance.

An ex-post view of inequality of opportunity in France and its regions

This paper proposes an ex-post measure of inequality of opportunity in France and its regions by assessing the inequality between individuals exerting the same effort. To this end, we define a fair income that fulfils ex-post equality of opportunity requirements. Unfairness is measured by an unfair Gini based on the distance between the actual income and the fair income. Our findings reveal that the measures of ex-post inequality of opportunity largely vary across regions, and that this is due to di_erences in reward schemes and in the impact of the non responsibility factors of income. We find that most regions have actual incomes closer to fair incomes than to average income, excepted Ile de France where the actual income looks poorly related to effort variables. Finally, we find that income inequality and inequality of opportunity are positively correlated among regions.

Genetics of cadmium hyperaccumulation variance in Thlaspi caerulescens natural populations from Switzerland

SUMMARY Heavy metal presence in the environment is a serious concern since some of them can be toxic to plants, animals and humans once accumulated along the food chain. Cadmium (Cd) is one of the most toxic heavy metal. It is naturally present in soils at various levels and its concentration can be increased by human activities. Several plants however have naturally developed strategies allowing them to grow on heavy metal enriched soils. One of them consists in the accumulation and sequestration of heavy metals in the above-ground biomass. Some plants present in addition an extreme strategy by which they accumulate a limited number of heavy metals in their shoots in amounts 100 times superior to those expected for a non-accumulating plant in the same conditions. Understanding the genetic basis of the hyperaccumulation trait - particularly for Cd - remains an important challenge which may lead to biotechnological applications in the soil phytoremediation. In this thesis, Thlaspi caerulescens J. & C. Presl (Brassicaceae) was used as a model plant to study the Cd hyperaccumulation trait, owing to its physiological and genetic characteristics. Twenty-four wild populations were sampled in different regions of Switzerland. They were characterized for environmental and soil parameters as well as intrinsic characteristics of plants (i.e. metal concentrations in shoots). They were as well genetically characterized by AFLPs, plastid DNA polymorphism and genes markers (CAPS and microsatellites) mainly developed in this thesis. Some of the investigated genes were putatively linked to the Cd hyperaccumulation trait. Since the study of the Cd hyperaccumulation in the field is important as it allows the identification of patterns of selection, the present work offered a methodology to define the Cd hyperaccumulation capacity of populations from different habitats permitting thus their comparison in the field. We showed that Cd, Zn, Fe and Cu accumulations were linked and that populations with higher Cd hyperaccumulation capacity had higher shoot and reproductive fitness. Using our genetic data, statistical methods (Beaumont & Nichols's procedure, partial Mantel tests) were applied to identify genomic signatures of natural selection related to the Cd hyperaccumulation capacity. A significant genetic difference between populations related to their Cd hyperaccumulation capacity was revealed based on somè specific markers (AFLP and candidate genes). Polymorphism at the gene encoding IRTl (Iron-transporter also participating to the transport of Zn) was suggested as explaining part of the variation in Cd hyperaccumulation capacity of populations supporting previous physiological investigations. RÉSUMÉ La présence de métaux lourds dans l'environnement est un phénomène préoccupant. En effet, certains métaux lourds - comme le cadmium (Cd) -sont toxiques pour les plantes, les animaux et enfin, accumulés le long de la chaîne alimentaire, pour les hommes. Le Cd est naturellement présent dans le sol et sa concentration peut être accrue par différentes activités humaines. Certaines plantes ont cependant développé des stratégies leur permettant de pousser sur des sols contaminés en métaux lourds. Parmi elles, certaines accumulent et séquestrent les métaux lourds dans leurs parties aériennes. D`autres présentent une stratégie encore plus extrême. Elles accumulent un nombre limité de métaux lourds en quantités 100 fois supérieures à celles attendues pour des espèces non-accumulatrices sous de mêmes conditions. La compréhension des bases génétiques de l'hyperaccumulation -particulièrement celle du Cd - représente un défi important avec des applications concrètes en biotechnologies, tout particulièrement dans le but appliqué de la phytoremediation des sols contaminés. Dans cette thèse, Thlaspi caerulescens J. & C. Presl (Brassicaceae) a été utilisé comme modèle pour l'étude de l'hyperaccumulation du Cd de par ses caractéristiques physiologiques et génétiques. Vingt-quatre populations naturelles ont été échantillonnées en Suisse et pour chacune d'elles les paramètres environnementaux, pédologique et les caractéristiques intrinsèques aux plantes (concentrations en métaux lourds) ont été déterminés. Les populations ont été caractérisées génétiquement par des AFLP, des marqueurs chloroplastiques et des marqueurs de gènes spécifiques, particulièrement ceux potentiellement liés à l'hyperaccumulation du Cd (CAPS et microsatellites). La plupart ont été développés au cours de cette thèse. L'étude de l'hyperaccumulation du Cd en conditions naturelles est importante car elle permet d'identifier la marque, éventuelle de sélection naturelle. Ce travail offre ainsi une méthodologie pour définir et comparer la capacité des populations à hyperaccumuler le Cd dans différents habitats. Nous avons montré que les accumulations du Cd, Zn, Fe et Cu sont liées et que les populations ayant une grande capacité d'hyperaccumuler le Cd ont également une meilleure fitness végétative et reproductive. Des méthodes statistiques (l'approche de Beaumont & Nichols, tests de Martel partiels) ont été utilisées sur les données génétiques pour identifier la signature génomique de la sélection naturelle liée à la capacité d'hyperaccumuler le Cd. Une différenciation génétique des populations liée à leur capacité d'hyperaccumuler le Cd a été mise en évidence sur certains marqueurs spécifiques. En accord avec les études physiologiques connues, le polymorphisme au gène codant IRT1 (un transporteur de Fe impliqué dans le transport du Zn) pourrait expliquer une partie de la variance de la capacité des populations à hyperaccumuler le Cd.

Can we continue to neglect genomic variation in introgression rates when inferring the history of speciation? A case study in a Mytilus hybrid zone.

The use of molecular data to reconstruct the history of divergence and gene flow between populations of closely related taxa represents a challenging problem. It has been proposed that the long-standing debate about the geography of speciation can be resolved by comparing the likelihoods of a model of isolation with migration and a model of secondary contact. However, data are commonly only fit to a model of isolation with migration and rarely tested against the secondary contact alternative. Furthermore, most demographic inference methods have neglected variation in introgression rates and assume that the gene flow parameter (Nm) is similar among loci. Here, we show that neglecting this source of variation can give misleading results. We analysed DNA sequences sampled from populations of the marine mussels, Mytilus edulis and M. galloprovincialis, across a well-studied mosaic hybrid zone in Europe and evaluated various scenarios of speciation, with or without variation in introgression rates, using an Approximate Bayesian Computation (ABC) approach. Models with heterogeneous gene flow across loci always outperformed models assuming equal migration rates irrespective of the history of gene flow being considered. By incorporating this heterogeneity, the best-supported scenario was a long period of allopatric isolation during the first three-quarters of the time since divergence followed by secondary contact and introgression during the last quarter. By contrast, constraining migration to be homogeneous failed to discriminate among any of the different models of gene flow tested. Our simulations thus provide statistical support for the secondary contact scenario in the European Mytilus hybrid zone that the standard coalescent approach failed to confirm. Our results demonstrate that genomic variation in introgression rates can have profound impacts on the biological conclusions drawn from inference methods and needs to be incorporated in future studies.

Novel markers of the human follicle-associated epithelium identified by genomic profiling and microdissection.

BACKGROUND & AIMS: Regulation of gene expression in the follicle-associated epithelium (FAE) over Peyer's patches is largely unknown. CCL20, a chemokine that recruits immature dendritic cells, is one of the few FAE-specific markers described so far. Lymphotoxin beta (LTalpha1beta2) expressed on the membrane of immune cells triggers CCL20 expression in enterocytes. In this study, we measured expression profiles of LTalpha1beta2-treated intestinal epithelial cells and selected CCL20 -coregulated genes to identify new FAE markers. METHODS: Genomic profiles of T84 and Caco-2 cell lines treated with either LTalpha1beta2, flagellin, or tumor necrosis factor alpha were measured using the Affymetrix GeneChip U133A. Clustering analysis was used to select CCL20 -coregulated genes, and laser dissection microscopy and real-time polymerase chain reaction on human biopsy specimens was used to assess the expression of the selected markers. RESULTS: Applying a 2-way analysis of variance, we identified regulated genes upon the different treatments. A subset of genes involved in inflammation and related to the nuclear factor kappaB pathway was coregulated with CCL20 . Among these genes, the antiapoptotic factor TNFAIP3 was highly expressed in the FAE. CCL23 , which was not coregulated in vitro with CCL20 , was also specifically expressed in the FAE. CONCLUSIONS: We have identified 2 novel human FAE specifically expressed genes. Most of the CCL20 -coregulated genes did not show FAE-specific expression, suggesting that other signaling pathways are critical to modulate FAE-specific gene expression.

Untangling the evolution of Rab G proteins: implications of a comprehensive genomic analysis.

BACKGROUND: Membrane-bound organelles are a defining feature of eukaryotic cells, and play a central role in most of their fundamental processes. The Rab G proteins are the single largest family of proteins that participate in the traffic between organelles, with 66 Rabs encoded in the human genome. Rabs direct the organelle-specific recruitment of vesicle tethering factors, motor proteins, and regulators of membrane traffic. Each organelle or vesicle class is typically associated with one or more Rab, with the Rabs present in a particular cell reflecting that cell's complement of organelles and trafficking routes. RESULTS: Through iterative use of hidden Markov models and tree building, we classified Rabs across the eukaryotic kingdom to provide the most comprehensive view of Rab evolution obtained to date. A strikingly large repertoire of at least 20 Rabs appears to have been present in the last eukaryotic common ancestor (LECA), consistent with the 'complexity early' view of eukaryotic evolution. We were able to place these Rabs into six supergroups, giving a deep view into eukaryotic prehistory. CONCLUSIONS: Tracing the fate of the LECA Rabs revealed extensive losses with many extant eukaryotes having fewer Rabs, and none having the full complement. We found that other Rabs have expanded and diversified, including a large expansion at the dawn of metazoans, which could be followed to provide an account of the evolutionary history of all human Rabs. Some Rab changes could be correlated with differences in cellular organization, and the relative lack of variation in other families of membrane-traffic proteins suggests that it is the changes in Rabs that primarily underlies the variation in organelles between species and cell types.

Genomic approaches to the study of HIV‐1 acquisition.

Host genome studies are increasingly available for the study of infectious disease susceptibility. Current technologies include large-scale genotyping, genome-wide screens such as transcriptome and silencing (silencing RNA) studies, and increasingly, the possibility to sequence complete genomes. These approaches are of interest for the study of individuals who remain uninfected despite documented exposure to human immunodeficiency virus type 1. The main limitation remains the ascertainment of exposure and establishing large cohorts of informative individuals. The pattern of enrichment for CCR5 Δ32 homozygosis should serve as the standard for assessing the extent to which a given cohort (of white subjects) includes a large proportion of exposed uninfected individuals.

The genomic impact of 100 million years of social evolution in seven ant species.

Ants (Hymenoptera, Formicidae) represent one of the most successful eusocial taxa in terms of both their geographic distribution and species number. The publication of seven ant genomes within the past year was a quantum leap for socio- and ant genomics. The diversity of social organization in ants makes them excellent model organisms to study the evolution of social systems. Comparing the ant genomes with those of the honeybee, a lineage that evolved eusociality independently from ants, and solitary insects suggests that there are significant differences in key aspects of genome organization between social and solitary insects, as well as among ant species. Altogether, these seven ant genomes open exciting new research avenues and opportunities for understanding the genetic basis and regulation of social species, and adaptive complex systems in general.

Avaluació de l’estructuració genètica humana a la Península Ibèrica i a d’altres regions del sudoest europeu: implicacions poblacionals i epidemiològiques

Els avenços en tècniques de genotipat de polimorfismes genètics a gran escala estan liderant una revolució en el camp de l’epidemiologia genètica i la genètica de poblacions humanes. La informació aportada per aquestes tècniques ha evidenciat l’existència d’estructuracions poblacionals que poden augmentar l’error en els estudis d’associació a escala genòmica (GWAS, genome-wide association studies). Estudis recents han demostrat la presència d’aquestes estructuracions a nivell interregional i intrarregional a Europa. El present projecte ha avaluat el grau d’estructuració genètica en poblacions de la Península Ibèrica i altres regions del sudoest europeu (Itàlia i França) per quantificar l’impacte que aquesta potencial estructuració pot tenir en el disseny d’estudis d’associació GWAS i reconstruir la història demogràfica de les poblacions de la Mediterrània. Per aconseguir aquests objectius, s’han analitzat mostres de DNA de 770 individus de 26 poblacions de la Península Ibèrica, França, Itàlia i d’altres països de la Mediterrània. Aquestes mostres van ser genotipades per 240000 SNPs utilitzant l’array 250K StyI d’Affymetrix en el marc d’aquest projecte o mitjançant altres arrays d’Affymetrix en els projectes internacionals HapMap i POPRES. S’han realitzat anàlisis estadístiques incloent anàlisis de components principals, Fst, identitat per descendència, desequilibri de lligament, barreres genètiques, etc. Aquests resultats han permés construir un marc de referència de la variabilitat en aquesta regió, avaluar el seu impacte en estudis d’associació i proposar mesures per evitar l’increment de qualsevol tipus d’error (tipus I i II) en estudis nacionals i internacionals. A més, també han permés reconstruir la història de les poblacions humanes de la Mediterrània així com analitzar les seves relacions demogràfiques. Donada la duració limitada d’aquesta acció (24 mesos, d’octubre de 2010 a setembre de 2012), els resultats d’aquest projecte es troben actualment en fase de redacció i conduiran a diverses publicacions en revistes internacionals i a la preparació de comunicacions a congressos.

Identification of Leishmania major cysteine proteinases as targets of the immune response in humans.

In this study, we report the identification of two parasite polypeptides recognized by human sera of patients infected with Leishmania major. Isolation and sequencing of the two genes encoding these polypeptides revealed that one of the genes is similar to the L. major cathepsin L-like gene family CPB, whereas the other gene codes for the L. major homologue of the cysteine proteinase a (CPA) of L. mexicana. By restriction enzyme digestion of genomic DNA, we show that the CPB gene is present in multiple copies in contrast to the cysteine proteinase CPA gene which could be unique. Specific antibodies directed against the mature regions of both types expressed in Escherichia coli were used to analyze the expression of these polypeptides in different stages of the parasite's life cycle. Polypeptides of 27 and 40 kDa in size, corresponding to CPA and CPB respectively, were detected at higher level in amastigotes than in stationary phase promastigotes. Purified recombinant CPs were also used to examine the presence of specific antibodies in sera from either recovered or active cases of cutaneous leishmaniasis patients. Unlike sera from healthy uninfected controls, all the sera reacted with recombinant CPA and CPB. This finding indicates that individuals having recovered from cutaneous leishmaniasis or with clinically apparent disease have humoral responses to cysteine proteinases demonstrating the importance of these proteinases as targets of the immune response and also their potential use for serodiagnosis.

Passive modernization? the new human development index and its components in Italy’s regions (1871-2007)

The article presents and discusses estimates of social and economic indicators for Italy’s regions in benchmark years roughly from Unification to the present day: life expectancy, education, GDP per capita at purchasing power parity, and the new Human Development Index (HDI). A broad interpretative hypothesis, based on the distinction between passive and active modernization, is proposed to account for the evolution of regional imbalances over the long-run. In the lack of active modernization, Southern Italy converged thanks to passive modernization, i.e., State intervention: however, this was more effective in life expectancy, less successful in education, expensive and as a whole ineffective in GDP. As a consequence, convergence in the HDI occurred from the late XIX century to the 1970s, but came to a sudden halt in the last decades of the XX century.

Molecular changes underlying mammalian phenotypic innovation

Mammals are characterized by specific phenotypic traits that include lactation, hair, and relatively large brains with unique structures. Individual mammalian lineages have, in turn, evolved characteristic traits that distinguish them from others. These include obvious anatom¬ical differences but also differences related to reproduction, life span, cognitive abilities, be¬havior. and disease susceptibility. However, the molecular basis of the diverse mammalian phenotypes and the selective pressures that shaped their evolution remain largely unknown. In the first part of my thesis, I analyzed the genetic factors associated with the origin of a unique mammalian phenotype lactation and I studied the selective pressures that forged the transition from oviparity to viviparity. Using a comparative genomics approach and evolutionary simulations, I showed that the emergence of lactation, as well as the appear¬ance of the casein gene family, significantly reduced selective pressure on the major egg-yolk proteins (the vitellogenin family). This led to a progressive loss of vitellogenins, which - in oviparous species - act as storage proteins for lipids, amino acids, phosphorous and calcium in the isolated egg. The passage to internal fertilization and placentation in therian mam¬mals rendered vitellogenins completely dispensable, which ended in the loss of the whole gene family in this lineage. As illustrated by the vitellogenin study, changes in gene content are one possible underlying factor for the evolution of mammalian-specific phenotypes. However, more subtle genomic changes, such as mutations in protein-coding sequences, can also greatly affect the phenotypes. In particular, it was proposed that changes at the level of gene reg¬ulation could underlie many (or even most) phenotypic differences between species. In the second part of my thesis, I participated in a major comparative study of mammalian tissue transcriptomes, with the goal of understanding how evolutionary forces affected expression patterns in the past 200 million years of mammalian evolution. I showed that, while com¬parisons of gene expressions are in agreement with the known species phylogeny, the rate of expression evolution varies greatly among lineages. Species with low effective population size, such as monotremes and hominoids, showed significantly accelerated rates of gene expression evolution. The most likely explanation for the high rate of gene expression evolution in these lineages is the accumulation of mildly deleterious mutations in regulatory regions, due to the low efficiency of purifying selection. Thus, our observations are in agreement with the nearly neutral theory of molecular evolution. I also describe substantial differences in evolutionary rates between tissues, with brain being the most constrained (especially in primates) and testis significantly accelerated. The rate of gene expression evolution also varies significantly between chromosomes. In particular, I observed an acceleration of gene expression changes on the X chromosome, probably as a result of adaptive processes associated with the origin of therian sex chromosomes. Lastly, I identified several individual genes as well as co-regulated expression modules that have undergone lineage specific expression changes and likely under¬lie various phenotypic innovations in mammals. The methods developed during my thesis, as well as the comprehensive gene content analyses and transcriptomics datasets made available by our group, will likely prove to be useful for further exploratory analyses of the diverse mammalian phenotypes.

Differential proteoglycan expression in two spinal cord regions after dorsal root injury.

Dorsal root injury leads to reactive gliosis in the spinal cord dorsal root entry zone and dorsal column, two regions that undergo Wallerian degeneration, but have distinct growth-inhibitory properties. This disparity could in part be due to differences in the number of degenerating sensory fibers, differences in glial cell activation, and/or to differential expression of growth-inhibitory molecules such as chondroitin sulfate proteoglycans. Laser capture microdissection of these two spinal cord white matter regions, followed by quantitative analysis of mRNA expression by real-time PCR, revealed that glial marker transcripts were differentially expressed post-injury and that the chondroitin sulfate proteoglycans Brevican and Versican V1 and V2 were preferentially up-regulated in the dorsal root entry zone, but not the dorsal column. These results indicate that reactive gliosis differs between these two regions and that Brevican and Versican are potential key molecules participating in the highly inhibitory properties of the dorsal root entry zone.

An assessment of gene prediction accuracy in large DNA sequences

One of the first useful products from the human genome will be a set of predicted genes. Besides its intrinsic scientific interest, the accuracy and completeness of this data set is of considerable importance for human health and medicine. Though progress has been made on computational gene identification in terms of both methods and accuracy evaluation measures, most of the sequence sets in which the programs are tested are short genomic sequences, and there is concern that these accuracy measures may not extrapolate well to larger, more challenging data sets. Given the absence of experimentally verified large genomic data sets, we constructed a semiartificial test set comprising a number of short single-gene genomic sequences with randomly generated intergenic regions. This test set, which should still present an easier problem than real human genomic sequence, mimics the approximately 200kb long BACs being sequenced. In our experiments with these longer genomic sequences, the accuracy of GENSCAN, one of the most accurate ab initio gene prediction programs, dropped significantly, although its sensitivity remained high. Conversely, the accuracy of similarity-based programs, such as GENEWISE, PROCRUSTES, and BLASTX was not affected significantly by the presence of random intergenic sequence, but depended on the strength of the similarity to the protein homolog. As expected, the accuracy dropped if the models were built using more distant homologs, and we were able to quantitatively estimate this decline. However, the specificities of these techniques are still rather good even when the similarity is weak, which is a desirable characteristic for driving expensive follow-up experiments. Our experiments suggest that though gene prediction will improve with every new protein that is discovered and through improvements in the current set of tools, we still have a long way to go before we can decipher the precise exonic structure of every gene in the human genome using purely computational methodology.