Large-scale genomic analyses link reproductive aging to hypothalamic signaling, breast cancer susceptibility and BRCA1-mediated DNA repair.

Menopause timing has a substantial impact on infertility and risk of disease, including breast cancer, but the underlying mechanisms are poorly understood. We report a dual strategy in ∼70,000 women to identify common and low-frequency protein-coding variation associated with age at natural menopause (ANM). We identified 44 regions with common variants, including two regions harboring additional rare missense alleles of large effect. We found enrichment of signals in or near genes involved in delayed puberty, highlighting the first molecular links between the onset and end of reproductive lifespan. Pathway analyses identified major association with DNA damage response (DDR) genes, including the first common coding variant in BRCA1 associated with any complex trait. Mendelian randomization analyses supported a causal effect of later ANM on breast cancer risk (∼6% increase in risk per year; P = 3 × 10(-14)), likely mediated by prolonged sex hormone exposure rather than DDR mechanisms.

Hybridization Capture Using RAD Probes (hyRAD), a New Tool for Performing Genomic Analyses on Collection Specimens.

In the recent years, many protocols aimed at reproducibly sequencing reduced-genome subsets in non-model organisms have been published. Among them, RAD-sequencing is one of the most widely used. It relies on digesting DNA with specific restriction enzymes and performing size selection on the resulting fragments. Despite its acknowledged utility, this method is of limited use with degraded DNA samples, such as those isolated from museum specimens, as these samples are less likely to harbor fragments long enough to comprise two restriction sites making possible ligation of the adapter sequences (in the case of double-digest RAD) or performing size selection of the resulting fragments (in the case of single-digest RAD). Here, we address these limitations by presenting a novel method called hybridization RAD (hyRAD). In this approach, biotinylated RAD fragments, covering a random fraction of the genome, are used as baits for capturing homologous fragments from genomic shotgun sequencing libraries. This simple and cost-effective approach allows sequencing of orthologous loci even from highly degraded DNA samples, opening new avenues of research in the field of museum genomics. Not relying on the restriction site presence, it improves among-sample loci coverage. In a trial study, hyRAD allowed us to obtain a large set of orthologous loci from fresh and museum samples from a non-model butterfly species, with a high proportion of single nucleotide polymorphisms present in all eight analyzed specimens, including 58-year-old museum samples. The utility of the method was further validated using 49 museum and fresh samples of a Palearctic grasshopper species for which the spatial genetic structure was previously assessed using mtDNA amplicons. The application of the method is eventually discussed in a wider context. As it does not rely on the restriction site presence, it is therefore not sensitive to among-sample loci polymorphisms in the restriction sites that usually causes loci dropout. This should enable the application of hyRAD to analyses at broader evolutionary scales.

Origin and Evolution of B Chromosomes in the Cichlid Fish Astatotilapia latifasciata Based on Integrated Genomic Analyses

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Biosynthesis of vitamins and cofactors in bacterium-harbouring trypanosomatids depends on the symbiotic association as revealed by genomic analyses

Some non-pathogenic trypanosomatids maintain a mutualistic relationship with a betaproteobacterium of the Alcaligenaceae family. Intensive nutritional exchanges have been reported between the two partners, indicating that these protozoa are excellent biological models to study metabolic co-evolution. We previously sequenced and herein investigate the entire genomes of five trypanosomatids which harbor a symbiotic bacterium (SHTs for Symbiont-Haboring Trypanosomatids) and the respective bacteria (TPEs for Trypanosomatid Proteobacterial Endosymbiont), as well as two trypanosomatids without symbionts (RTs for Regular Trypanosomatids), for the presence of genes of the classical pathways for vitamin biosynthesis. Our data show that genes for the biosynthetic pathways of thiamine, biotin, and nicotinic acid are absent from all trypanosomatid genomes. This is in agreement with the absolute growth requirement for these vitamins in all protozoa of the family. Also absent from the genomes of RTs are the genes for the synthesis of pantothenic acid, folic acid, riboflavin, and vitamin B6. This is also in agreement with the available data showing that RTs are auxotrophic for these essential vitamins. On the other hand, SHTs are autotrophic for such vitamins. Indeed, all the genes of the corresponding biosynthetic pathways were identified, most of them in the symbiont genomes, while a few genes, mostly of eukaryotic origin, were found in the host genomes. The only exceptions to the latter are: the gene coding for the enzyme ketopantoate reductase (EC:1.1.1.169) which is related instead to the Firmicutes bacteria; and two other genes, one involved in the salvage pathway of pantothenic acid and the other in the synthesis of ubiquinone, that are related to Gammaproteobacteria. Their presence in trypanosomatids may result from lateral gene transfer. Taken together, our results reinforce the idea that the low nutritional requirement of SHTs is associated with the presence of the symbiotic bacterium, which contains most genes for vitamin production.

Genomic analyses of flowering transition and inflorescence development in grapevine = Análisis genómicos de la transición floral y del desarrollo de la inflorescencia en la vid

The general objective of this work is to analyze the regulatory processes underlying flowering transition and inflorescence and flower development in grapevine. Most of these crucial developmental events take place within buds growing during two seasons in two consecutive years. During the first season, the shoot apical meristem within the bud differentiates all the basic elements of the shoot including flowering transition in lateral primordia and development of inflorescence primordia. These events practically end with bud dormancy. The second season, buds resume shoot growth associated to flower formation and development. In grapevine, the lateral meristems can give rise either to tendril or inflorescence primordia that are homologous organs. With this purpose, we performed global transcriptome analyses along the bud annual cycle and during inflorescence and tendril development. In addition, we approach the genomic analysis of the MIKC type MADS-box gene family in grapevine to identify all its members and assign them putative biological functions. Regarding buds developmental cycle, the results indicate that the main factors explaining the global gene expression differences were the processes of bud dormancy and active growth as well as stress responses. Non dormant buds exhibited up-regulation in functional categories typical of actively proliferating and growing cells (photosynthesis, cell cycle regulation, chromatin assembly) whereas in dormant ones the main functional categories up-regulated were associated to stress response pathways together with transcripts related to starch catabolism. Major transcriptional changes during the dormancy period were associated to the para/endodormancy, endo/ecodormancy and ecodormancy/bud break transitions. Global transcriptional analyses along tendril and inflorescence development suggested that these two homologous organs share a common transcriptional program related to cell proliferation functions. Both structures showed a progressive decrease in the expression of categories such as cell-cycle, auxin metabolism/signaling, DNA metabolism, chromatin assembly and a cluster of five transcripts belonging to the GROWTH-REGULATING FACTOR (GRF) transcription factor family, that are known to control cell proliferation in other species and determine the size of lateral organs. However, they also showed organ specific transcriptional programs that can be related to their differential organ structure and function. Tendrils showed higher transcription of genes related to photosynthesis, hormone signaling and secondary metabolism than inflorescences, while inflorescences have higher transcriptional activity for genes encoding transcription factors (especially those belonging to the MADS-box gene family). Further analysis along inflorescence development evidenced the relevance of additional functions likely related to processes of flower development such as fatty acid and lipid metabolism, jasmonate signaling and oxylipin biosynthesis. The transcriptional analyses performed highlighted the relevance of several groups of transcriptional regulators in the developmental processes studied. The expression profiles along bud development revealed significant differences for some MADS-box subfamilies in relation to other plant species, like the members of the FLC and SVP subfamilies suggesting new roles for these groups in grapevine. In this way, it was found that VvFLC2 and VvAGL15.1 could participate, together with some members of the SPL-L family, in dormancy regulation, as was shown for some of them in other woody plants. Similarly, the expression patterns of the VvFLC1, VvFUL, VvSOC1.1 (together with VvFT, VvMFT1 and VFL) genes could indicate that they play a role in flowering transition in grapevine, in parallel to their roles in other plant systems. The expression levels of VFL, the grapevine LEAFY homolog, could be crucial to specify the development of inflorescence and flower meristems instead of tendril meristems. MADS-box genes VvAP3.1 and 2, VvPI, VvAG1 and 3, VvSEP1-4, as well as VvBS1 and 2 are likely associated with the events of flower meristems and flower organs differentiation, while VvAP1 and VvFUL-L (together with VvSOC1.1, VvAGL6.2) could be involved on tendril development given their expression patterns. In addition, the biological function ofVvAP1 and VvTFL1A was analyzed using a gene silencing approach in transgenic grapevine plants. Our preliminary results suggested a possible role for both genes in the initiation and differentiation of tendrils. Finally, the genomic analysis of the MADS-box gene family in grapevine revealed differential features regarding number and expression pattern of genes putatively involved in the flowering transition process as compared to those involved in the specification of flower and fruit organ identity. Altogether, the results obtained allow identifying putative candidate genes and pathways regulating grapevine reproductive developmental processes paving the way to future experiments demonstrating specific gene biological functions. RESUMEN El objetivo general de este trabajo es analizar los procesos regulatorios subyacentes a la inducción floral así como al desarrollo de la inflorescencia y la flor en la vid. La mayor parte de estos eventos cruciales tienen lugar en las yemas a lo largo de dos estaciones de crecimiento consecutivas. Durante la primera estación, el meristemo apical contenido en la yema diferencia los elementos básicos del pámpano, lo cual incluye la inducción de la floración en los meristemos laterales y el subsiguiente desarrollo de primordios de inflorescencia. Estos procesos prácticamente cesan con la entrada en dormición de la yema. En la segunda estación, se reanuda el crecimiento del pámpano acompañado por la formación y desarrollo de las flores. En la vid, los meristemos laterales pueden dar lugar a primordios de inflorescencia o de zarcillo que son considerados órganos homólogos. Con este objetivo llevamos a cabo un estudio a nivel del transcriptoma de la yema a lo largo de su ciclo anual, así como a lo largo del desarrollo de la inflorescencia y del zarcillo. Además realizamos un análisis genómico de la familia MADS de factores transcripcionales (concretamente aquellos del tipo MIKC) para identificar todos sus miembros y tratar de asignarles posibles funciones biológicas. En cuanto al ciclo de desarrollo de la yema, los resultados indican que los principales factores que explican las diferencias globales en la expresión génica fueron los procesos de dormición de la yema y el crecimiento activo junto con las respuestas a diversos tipos de estrés. Las yemas no durmientes mostraron un incremento en la expresión de genes contenidos en categorías funcionales típicas de células en proliferación y crecimiento activo (como fotosíntesis, regulación del ciclo celular, ensamblaje de cromatina), mientras que en las yemas durmientes, las principales categorías funcionales activadas estaban asociadas a respuestas a estrés, así como con el catabolismo de almidón. Los mayores cambios observados a nivel de transcriptoma en la yema coincidieron con las transiciones de para/endodormición, endo/ecodormición y ecodormición/brotación. Los análisis transcripcionales globales a lo largo del desarrollo del zarcillo y de la inflorescencia sugirieron que estos dos órganos homólogos comparten un programa transcripcional común, relacionado con funciones de proliferación celular. Ambas estructuras mostraron un descenso progresivo en la expresión de genes pertenecientes a categorías funcionales como regulación del ciclo celular, metabolismo/señalización por auxinas, metabolismo de ADN, ensamblaje de cromatina y un grupo de cinco tránscritos pertenecientes a la familia de factores transcripcionales GROWTH-REGULATING FACTOR (GRF), que han sido asociados con el control de la proliferación celular y en determinar el tamaño de los órganos laterales en otras especies. Sin embargo, también pusieron de manifiesto programas transcripcionales que podrían estar relacionados con la diferente estructura y función de dichos órganos. Los zarcillos mostraron mayor actividad transcripcional de genes relacionados con fotosíntesis, señalización hormonal y metabolismo secundario que las inflorescencias, mientras que éstas presentaron mayor actividad transcripcional de genes codificantes de factores de transcripción (especialmente los pertenecientes a la familia MADS-box). Análisis adicionales a lo largo del desarrollo de la inflorescencia evidenciaron la relevancia de otras funciones posiblemente relacionadas con el desarrollo floral, como el metabolismo de lípidos y ácidos grasos, la señalización mediada por jasmonato y la biosíntesis de oxilipinas. Los análisis transcripcionales llevados a cabo pusieron de manifiesto la relevancia de varios grupos de factores transcripcionales en los procesos estudiados. Los perfiles de expresión estudiados a lo largo del desarrollo de la yema mostraron diferencias significativas en algunas de las subfamilias de genes MADS con respecto a otras especies vegetales, como las observadas en los miembros de las subfamilias FLC y SVP, lo cual sugiere que podrían desempeñar nuevas funciones en la vid. En este sentido, se encontró que los genes VvFLC2 y VvAGL15.1 podrían participar, junto con algunos miembros de la familia SPL-L, en la regulación de la dormición. De un modo similar, los patrones de expresión de los genes VvFLC1, VvFUL, VvSOC1.1 (junto con VvFT, VvMFT1 y VFL) podría indicar que desempeñan un papel en la regulación de la inducción de la floración en la vid, como se ha observado en otros sistemas vegetales. Los niveles de expresión de VFL, el homólogo en vid del gen LEAFY de A. thaliana podrían ser cruciales para la especificación del desarrollo de meristemos de inflorescencia y flor en lugar de meristemos de zarcillo. Los genes VvAP3.1 y 2, VvPI, VvAG1 y 3, VvSEP1-4, así como VvBS1 y 2 parecen estar asociados con los eventos de diferenciación de meristemos y órganos florales, mientras que VvAP1 y VvFUL-L (junto con VvSOC1.1 y VvAGL6.2) podrían estar implicados en el desarrollo del zarcillo dados sus patrones de expresión. Adicionalmente, se analizó la función biológica de los genes VvAP1 y VvTFL1A por medio de una estrategia de silenciamiento génico. Los datos preliminares sugieren un posible papel para ambos genes en la iniciación y diferenciación de los zarcillos. Finalmente, el análisis genómico de la familia MADS en vid evidenció diferencias con respecto a otras especies vegetales en cuanto a número de miembros y patrón de expresión en genes supuestamente implicados en la inducción de la floración, en comparación con aquellos relacionados con la especificación de identidad de órganos florales y desarrollo del fruto. En conjunto, los resultados obtenidos han permitido identificar posibles rutas y genes candidatos a participar en la regulación de los procesos de desarrollo reproductivo de la vid, sentando las bases de futuros experimentos encaminados a conocer la funciones biológicas de genes específicos.

Independent and combined analyses of sequences from all three genomic compartments converge on the root of flowering plant phylogeny

Plant phylogenetic estimates are most likely to be reliable when congruent evidence is obtained independently from the mitochondrial, plastid, and nuclear genomes with all methods of analysis. Here, results are presented from separate and combined genomic analyses of new and previously published data, including six and nine genes (8,911 bp and 12,010 bp, respectively) for different subsets of taxa that suggest Amborella + Nymphaeales (water lilies) are the first-branching angiosperm lineage. Before and after tree-independent noise reduction, most individual genomic compartments and methods of analysis estimated the Amborella + Nymphaeales basal topology with high support. Previous phylogenetic estimates placing Amborella alone as the first extant angiosperm branch may have been misled because of a series of specific problems with paralogy, suboptimal outgroups, long-branch taxa, and method dependence. Ancestral character state reconstructions differ between the two topologies and affect inferences about the features of early angiosperms.

An integrated genomic analysis of human glioblastoma Multiforme

Glioblastoma multiforme ( GBM) is the most common and lethal type of brain cancer. To identify the genetic alterations in GBMs, we sequenced 20,661 protein coding genes, determined the presence of amplifications and deletions using high- density oligonucleotide arrays, and performed gene expression analyses using next- generation sequencing technologies in 22 human tumor samples. This comprehensive analysis led to the discovery of a variety of genes that were not known to be altered in GBMs. Most notably, we found recurrent mutations in the active site of isocitrate dehydrogenase 1 ( IDH1) in 12% of GBM patients. Mutations in IDH1 occurred in a large fraction of young patients and in most patients with secondary GBMs and were associated with an increase in overall survival. These studies demonstrate the value of unbiased genomic analyses in the characterization of human brain cancer and identify a potentially useful genetic alteration for the classification and targeted therapy of GBMs.

I. Identification and characterisation of epigenetically altered tumor supressor genes by proteomic and genomic approches / II. Studies of genes regualted by MeCP2 union to its promoter regions by proteomic and genmomic approches

DNA methylation has an important impact on normal cell physiology, thus any defects in this mechanism may be related to the development of various diseases In this project we are interested in identifying epigeneticaliy modified genes, in general controlled by processes related to the DNA methylation, by means of a new strategy combining protomic and genomic analyses. First, the two Dimensional-Difference Gel Electrophoresis (2-DIGE) protein analyses of extracts obtained from HCT-116 wt and double knockout for DNMT1 and DNMT3b (DKO) cells revealed 34 proteins overexpressed in the condition of DNMTs depletion. From five genes with higher transcript lavels in DKO cells, comparing with HCT-116 wt. oniy AKR1B1, UCHLl and VIM are melhylated in HCT-116. As expected. the DNA methvlation 1s lost in DKO cells. The rneth,vl ation of VIM and UCHLl promoters in some cancer samples has already been repaired, thus further studies has been focused on AKRlBI. AKR1B1 expression due lo DNA methyiaton of promoter region seems to occur specilfically in the colon cancer cell Iines. which was confirmed in the DNA rnethylation status and expression analyses. performed on 32 different cancer cell lines (including colon, breast, lymphoma, leukemia, neuroblastoma, glioma and lung cancer cell Iines) as well as normal colon and normal lymphocytes samples. AKRIBI expression after treatments with DNA demethvlating agent (AZA) was rescued in 5 coloncancer cell lines (including genetic regulation of the candidate gene. The methylation status of the rest of the genes identified in proteomic analysis was checked by methylation specific PCR (MSP) experiment and all appeared to be unmethylated. The similar research has been done also bv means of Mecp2-null mouse model For 14 selected candidate genes the analyses of expression leveis, methylation Status and MeCP2 interaction with promoters are currently being performed.

Genomic analysis for managing small and endangered populations: a case study in Tyrol Grey cattle

Analysis of genomic data is increasingly becoming part of the livestock industry. Therefore, the routine collection of genomic information would be an invaluable resource for effective management of breeding programs in small, endangered populations. The objective of the paper was to demonstrate how genomic data could be used to analyse (1) linkage disequlibrium (LD), LD decay and the effective population size (NeLD); (2) Inbreeding level and effective population size (NeROH) based on runs of homozygosity (ROH); (3) Prediction of genomic breeding values (GEBV) using small within-breed and genomic information from other breeds. The Tyrol Grey population was used as an example, with the goal to highlight the potential of genomic analyses for small breeds. In addition to our own results we discuss additional use of genomics to assess relatedness, admixture proportions, and inheritance of harmful variants. The example data set consisted of 218 Tyrol Grey bull genotypes, which were all available AI bulls in the population. After standard quality control restrictions 34,581 SNPs remained for the analysis. A separate quality control was applied to determine ROH levels based on Illumina GenCall and Illumina GenTrain scores, resulting into 211 bulls and 33,604 SNPs. LD was computed as the squared correlation coefficient between SNPs within a 10 mega base pair (Mb) region. ROHs were derived based on regions covering at least 4, 8, and 16 Mb, suggesting that animals had common ancestors approximately 12, 6, and 3 generations ago, respectively. The corresponding mean inbreeding coefficients (F ROH) were 4.0% for 4 Mb, 2.9% for 8 Mb and 1.6% for 16 Mb runs. With an average generation interval of 5.66 years, estimated NeROH was 125 (NeROH>16 Mb), 186 (NeROH>8 Mb) and 370 (NeROH>4 Mb) indicating strict avoidance of close inbreeding in the population. The LD was used as an alternative method to infer the population history and the Ne. The results show a continuous decrease in NeLD, to 780, 120, and 80 for 100, 10, and 5 generations ago, respectively. Genomic selection was developed for and is working well in large breeds. The same methodology was applied in Tyrol Grey cattle, using different reference populations. Contrary to the expectations, the accuracy of GEBVs with very small within breed reference populations were very high, between 0.13-0.91 and 0.12-0.63, when estimated breeding values and deregressed breeding values were used as pseudo-phenotypes, respectively. Subsequent analyses confirmed the high accuracies being a consequence of low reliabilities of pseudo-phenotypes in the validation set, thus being heavily influenced by parent averages. Multi-breed and across breed reference sets gave inconsistent and lower accuracies. Genomic information may have a crucial role in management of small breeds, even if its primary usage differs from that of large breeds. It allows to assess relatedness between individuals, trends in inbreeding and to take decisions accordingly. These decisions would be based on the real genome architecture, rather than conventional pedigree information, which can be missing or incomplete. We strongly suggest the routine genotyping of all individuals that belong to a small breed in order to facilitate the effective management of endangered livestock populations.

Endosymbiosis in trypanosomatids: the genomic cooperation between bacterium and host in the synthesis of essential amino acids is heavily influenced by multiple horizontal gene transfers

Background Trypanosomatids of the genera Angomonas and Strigomonas live in a mutualistic association characterized by extensive metabolic cooperation with obligate endosymbiotic Betaproteobacteria. However, the role played by the symbiont has been more guessed by indirect means than evidenced. Symbiont-harboring trypanosomatids, in contrast to their counterparts lacking symbionts, exhibit lower nutritional requirements and are autotrophic for essential amino acids. To evidence the symbiont’s contributions to this autotrophy, entire genomes of symbionts and trypanosomatids with and without symbionts were sequenced here. Results Analyses of the essential amino acid pathways revealed that most biosynthetic routes are in the symbiont genome. By contrast, the host trypanosomatid genome contains fewer genes, about half of which originated from different bacterial groups, perhaps only one of which (ornithine cyclodeaminase, EC:4.3.1.12) derived from the symbiont. Nutritional, enzymatic, and genomic data were jointly analyzed to construct an integrated view of essential amino acid metabolism in symbiont-harboring trypanosomatids. This comprehensive analysis showed perfect concordance among all these data, and revealed that the symbiont contains genes for enzymes that complete essential biosynthetic routes for the host amino acid production, thus explaining the low requirement for these elements in symbiont-harboring trypanosomatids. Phylogenetic analyses show that the cooperation between symbionts and their hosts is complemented by multiple horizontal gene transfers, from bacterial lineages to trypanosomatids, that occurred several times in the course of their evolution. Transfers occur preferentially in parts of the pathways that are missing from other eukaryotes. Conclusion We have herein uncovered the genetic and evolutionary bases of essential amino acid biosynthesis in several trypanosomatids with and without endosymbionts, explaining and complementing decades of experimental results. We uncovered the remarkable plasticity in essential amino acid biosynthesis pathway evolution in these protozoans, demonstrating heavy influence of horizontal gene transfer events, from Bacteria to trypanosomatid nuclei, in the evolution of these pathways.

Evolutionary genomics and HIV restriction factors.

PURPOSE OF REVIEW: To provide updated insights into innate antiviral immunity and highlight prototypical evolutionary features of well characterized HIV restriction factors. RECENT FINDINGS: Recently, a new HIV restriction factor, Myxovirus resistance 2, has been discovered and the region/residue responsible for its activity identified using an evolutionary approach. Furthermore, IFI16, an innate immunity protein known to sense several viruses, has been shown to contribute to the defense to HIV-1 by causing cell death upon sensing HIV-1 DNA. SUMMARY: Restriction factors against HIV show characteristic signatures of positive selection. Different patterns of accelerated sequence evolution can distinguish antiviral strategies--offense or defence--as well as the level of specificity of the antiviral properties. Sequence analysis of primate orthologs of restriction factors serves to localize functional domains and sites responsible for antiviral action. We use recent discoveries to illustrate how evolutionary genomic analyses help identify new antiviral genes and their mechanisms of action.

Mitochondrial genome sequence diversity of Indian Plasmodium falciparum isolates

We have analysed the whole mitochondrial (mt) genome sequences (each ~6 kilo nucleotide base pairs in length) of four field isolates of the malaria parasite Plasmodium falciparum collected from different locations in India. Comparative genomic analyses of mt genome sequences revealed three novel India-specific single nucleotide polymorphisms. In general, high mt genome diversity was found in Indian P. falciparum, at a level comparable to African isolates. A population phylogenetic tree placed the presently sequenced Indian P. falciparum with the global isolates, while a previously sequenced Indian isolate was an outlier. Although this preliminary study is limited to a few numbers of isolates, the data have provided fundamental evidence of the mt genome diversity and evolutionary relationships of Indian P. falciparum with that of global isolates.

SlimCodeML: An Optimized Version of CodeML for the Branch-Site Model

CodeML (part of the PAML package) im- plements a maximum likelihood-based approach to de- tect positive selection on a specific branch of a given phylogenetic tree. While CodeML is widely used, it is very compute-intensive. We present SlimCodeML, an optimized version of CodeML for the branch-site model. Our performance analysis shows that SlimCodeML substantially outperforms CodeML (up to 9.38 times faster), especially for large-scale genomic analyses.

VARIABILITY OF PHENOTYPIC, PROTEOMIC AND TRANSCRIPTOMIC EXPRESSION OF STAPHYLOCOCCUS AUREUS SURFACE ADHESINS

Staphylococcus aureus is a highly successful pathogen responsible of a wide variety of diseases, from minor skin infection to life-threatening sepsis or infective endocarditis, as well as food poisoning and toxic shock syndrome. This heterogeneity of infections and the ability of S. aureus to develop antibiotic-resistance to virtually any available drugs reflect its extraordinary capacity to adapt and survive in a great variety of environments. The pathogenesis of S. aureus infection involves a wide range of cell wall-associated adhesins and extracellular toxins that promote host colonization and invasion. In addition, S. aureus is extremely well equipped with regulatory systems that sense environmental conditions and respond by fine tuning the expression of metabolic and virulence determinants. Surface adhesins referred to MSCRAMMs - for Microbial Surface Component Recognizing Adherence Matrix Molecules - mediate binding to the host extracellular matrix or serum components, including fibrinogen, fibronectin, collagen and elastin, and promote tissue colonization and invasion. Major MSCRAMMs include a family of surface-attached proteins covalently bound to the cell wall peptidoglycan via a conserved LPXTG motif. Genomic analyses indicate that S. aureus contain up to 22 LPXTG surface proteins, which could potentially act individually or in synergy to promote infection. In the first part of this study we determined the range of adherence phenotypes to fibrinogen and fibronectin among 30 carriage isolates of S. aureus and compared it to the adherence phenotypes of 30 infective endocarditis and 30 blood culture isolates. Overall there were great variations in in vitro adherence, but no differences were observed between carriage and infection strains. We further determined the relation between in vitro adherence and in vivo infectivity in a rat model of experimental endocarditis, using 4 isolates that displayed either extremely low or high adherence phenotypes. Unexpectedly, no differences were observed between the in vivo infectivity of isolates that were poorly and highly adherent in vitro. We concluded that the natural variability of in vitro adherence to fibrinogen and fibronectin did not correlate with in vivo infectivity, and thus that pathogenic differences between various strains might only be expressed in in vivo conditions, but not in vitro. Therefore, considering the importance of adhesins expression for infection, direct measurement of those adhesins present on the bacterial surface were made by proteomic approach. 5 In the second series of experiments we assessed the physical presence of the LPXTG species at the staphylococcal surface, as measured at various time points during growth in different culture media. S. aureus Newman was grown in either tryptic soy broth (TSB) or in Roswell Park Memorial Institute (RPMI) culture medium, and samples were removed from early exponential growth phase to late stationary phase. Experiments were performed with mutants in the global accessory-gene regulator (agr), surface protein A (Spa) and clumping factor A (ClfA). Peptides of surface proteins were recovered by "trypsin-shaving" of live bacteria, and semi-quantitative proteomic analysis was performed by tandem liquid-chromatography and mass-spectrometry (LC-MS). We also determined in parallel the mRNA expression by microarrays analysis, as well as the phenotypic adherence of the bacteria to fibrinogen in vitro. The surface proteome was highly complex and contained numerous proteins theoretically not belonging to the bacterial envelope, including ribosomal proteins and metabolic enzymes. Sixteen of the 21 known LPXTG species were detected, but were differentially expressed. As expected, 9 known agr-regulated proteins (e.g. including Spa, FnBPA, ClfA, IsdA, IsdB, SasH, SasD, SasG and FmtB) increased up to the late exponential growth phase, and were abrogated in agr-negative mutants. However, only Spa and SasH modified their proteomic and mRNA profiles in parallel in the parent and its agr negative mutant, while all other LPXTG proteins modified their proteomic profiles independently of their mRNA. Moreover, ClfA became highly transcribed and active in in vitro fibrinogen adherence tests during late growth (24h), whereas it remained poorly detected by proteomics. Differential expression was also detected in iron-rich TSB versus iron-poor RPMI. Proteins from the iron-regulated surface determinant (isd) system, including IsdA, IsdB and IsdH were barely expressed in iron-rich TSB, whereas they increased their expression by >10 time in iron-poor RPMI. We conclude that semi-quantitative proteomic analysis of specific protein species is feasible in S. aureus and that proteomic, transcriptomic and adherence phenotypes demonstrated differential profiles in S. aureus. Furthermore, peptide signatures released by trypsin shaving suggested differential protein domain exposures in various environments, which might be relevant for antiadhesins vaccines. A comprehensive understanding of the S. aureus physiology should integrate all these approaches.

Molecular evolution of the Drosophila olfactory system

The olfactory system is an attractive model to study the genetic mechanisms underlying evolution of the nervous system. This sensory system mediates the detection and behavioural responses to an enormous diversity of volatile chemicals in the environment and displays rapid evolution, as species acquire, modify and discard olfactory receptors and circuits to adapt to new olfactory stimuli. Drosophilids provide an attractive model to study these processes. The availability of 12 sequenced genomes of Drosophila species occupying diverse ecological niches provides a rich resource for genomic analyses. Moreover, one of these species, Drosophila melanogaster, is amenable to a powerful combination of genetic and electrophysiological analyses. D. melanogaster has two distinct families of olfactory receptors to detect odours, the well-characterised Odorant Receptors (ORs) and the recently identified lonotropic Receptors (IRs). In my thesis, I have provided new insights into the genetic mechanisms underlying olfactory system evolution through three distinct, but interrelated projects. First, I performed a comparative genomic analysis of the IR repertoire in 12 sequenced Drosophila species, which has revealed that the olfactory IRs are highly conserved across species. By contrast, a large fraction of IRs that are not expressed in the olfactory system - and which may be gustatory receptors - are much more variable in sequence and gene copy number. Second, to identify ligands for IR expressing olfactory sensory neurons, I have performed an electrophysiological screen in D. melanogaster using a panel of over 160 odours. I found that the IRs respond to a number of amines, aldehydes and acids, contrasting with the chemical specificity of the OR repertoire, which is mainly tuned to esters, alcohols and ketones. Finally, the identification of ligands for IRs in this species allowed me to investigate in detail the molecular and functional evolution of a tandem array of IRs, IR75a/IR75b/IR75c, in D. sechellia. This species is endemic to the Seychelles archipelago and highly specialised to breed on the fruits of Morinda citrifolia, which is repulsive and toxic for other Drosophila species. These studies led me to discover that receptor loss, changes in receptor specificity and changes in receptor expression have likely played an important role during the evolution of these IRs in D. sechellia. These changes may explain, in part, the unique chemical ecology of this species. - Le système olfactif est un excellent modèle pour étudier les mécanismes génétiques impliqués dans l'étude de l'évolution du système nerveux. Ce système sensoriel permet la détection de nombreux composés volatils présents dans l'environnement et est à la base des réponses comportementales. Il est propre à chaque espèce et évolue rapidement en modifiant ou en éliminant des récepteurs et leurs circuits olfactifs correspondants pour s'adapter à de nouvelles odeurs. Pour étudier le système olfactif et son évolution, nous avons décidé d'utiliser la drosophile comme modèle. Le séquençage complet de 12 souches de drosophiles habitant différentes niches écologiques permet une analyse génomique conséquente. De plus, l'une de ces espèces Drosophila melanogaster permet la combinaison d'analyses génétiques et électrophysiologiques. En effet, D. melanogaster possède 2 familles distinctes de récepteurs olfactifs qui permettent la détection d'odeurs: les récepteurs olfactifs (ORs) étant les mieux caractérisés et les récepteurs ionotropiques (IRs), plus récemment identifiés. Au cours de ma thèse, j'ai apporté des nouvelles connaissances qui m'ont permis de mieux comprendre les mécanismes génétiques à la base de l'évolution du système olfactif au travers de trois projets différents, mais interdépendants. Premièrement, j'ai réalisé une analyse génomique comparative de l'ensemble des IRs dans les 12 souches de drosophiles séquencées jusqu'à présent. Ceci a montré que les récepteurs olfactifs IRs sont hautement conservés parmi l'ensemble de ces espèces. Au contraire, une grande partie des IRs qui ne sont pas exprimés dans le système olfactif, et qui semblent être des récepteurs gustatifs, sont beaucoup plus variables dans leur séquence et dans le nombre de copie de gènes. Deuxièmement, pour identifier les ligands des récepteurs IRs exprimés par les neurones sensoriels olfactifs, j'ai réalisé une étude électrophysiologique chez D. melanogaster e η testant l'effet de plus de 160 composés chimiques sur les IRs. J'ai trouvé que les IRs répondent à un nombre d'amines, d'aldéhydes et d'acides, contrairement aux récepteurs olfactifs ORs qui eux répondent principalement aux esthers, alcools et cétones. Finalement, l'identification de ligands pour les IRs dans ces espèces m'a permis d'étudier en détail l'évolution fonctionnelle et moléculaire des IR75a/IR75b/IR75c dans D. sechellia. Cette espèce est endémique de l'archipel des Seychelles et se nourrit spécifiquement du fruit Morinda citrifolia qui est répulsif et toxique pour d'autres souches de drosophiles. Ces études m'ont poussé à découvrir que, la perte de IR75a, le changement dans la spécificité de IR75b ainsi que le changement dans l'expression de IR75c ont probablement joué un rôle important dans l'évolution des IRs chez D. sechellia. Ces changements peuvent expliquer, en partie, l'écologie chimique propre à cette espèce. Résumé français large public Le système olfactif permet aux animaux de détecter des milliers de molécules odorantes, les aidant ainsi à trouver de la nourriture, à distinguer si elle est fraîche ou avariée, à trouver des partenaires sexuels, ainsi qu'à éviter les prédateurs. Selon l'environnement et le mode de vie des espèces, le système olfactif doit détecter des odeurs très diverses ; en effet, un moustique qui recherche du sang humain pour se nourrir doit détecter des odeurs bien différentes d'une abeille qui recherche des fleurs. Dans ma thèse, j'ai essayé de comprendre comment les systèmes olfactifs d'une espèce évoluent pour s'adapter aux exigences induites par son environnement. Un très bon modèle pour étudier cela est la drosophile dont les différentes espèces se nichent dans des habitats très divers. Pour ce faire, j'ai étudié les récepteurs olfactifs de différentes espèces de la drosophile. Ces récepteurs sont des protéines qui se lient à des odeurs spécifiques. Lorsqu'ils se lient, ils activent un neurone qui envoie un signal électrique au cerveau. Ce signal est ensuite traité par ce dernier qui indique à la mouche si l'odeur est attractive ou répulsive. J'ai identifié les récepteurs olfactifs de plusieurs espèces de drosophile et étudié s'il y avait des différences entre elles. La plupart des récepteurs sont similaires entre les espèces, cependant dans l'une d'entre elles, certains récepteurs sont différents. Ce fait est particulièrement intéressant car cette espèce de drosophile se nourrit de fruits que les autres espèces n'apprécient pas. Comme nous ne savons pas quels récepteurs se lient à quelles odeurs, j'ai testé un grand nombre de composants odorants. Ceci m'a permis de constater que, effectivement, certains changements produits dans ces récepteurs expliquent pourquoi cette espèce aime particulièrement ces fruits. En outre, mes résultats contribuent à mieux comprendre les changements génétiques qui sont impliqués dans l'évolution du système olfactif.