Systems biology of plants : from intraspecific genome variation to computational morphodynamics

Recently, the introduction of second generation sequencing and further advance-ments in confocal microscopy have enabled system-level studies for the functional characterization of genes. The degree of complexity intrinsic to these approaches needs the development of bioinformatics methodologies and computational models for extracting meaningful biological knowledge from the enormous amount of experi¬mental data which is continuously generated. This PhD thesis presents several novel bioinformatics methods and computational models to address specific biological questions in Plant Biology by using the plant Arabidopsis thaliana as a model system. First, a spatio-temporal qualitative analysis of quantitative transcript and protein profiles is applied to show the role of the BREVIS RADIX (BRX) protein in the auxin- cytokinin crosstalk for root meristem growth. Core of this PhD work is the functional characterization of the interplay between the BRX protein and the plant hormone auxin in the root meristem by using a computational model based on experimental evidence. Hyphotesis generated by the modelled to the discovery of a differential endocytosis pattern in the root meristem that splits the auxin transcriptional response via the plasma membrane to nucleus partitioning of BRX. This positional information system creates an auxin transcriptional pattern that deviates from the canonical auxin response and is necessary to sustain the expression of a subset of BRX-dependent auxin-responsive genes to drive root meristem growth. In the second part of this PhD thesis, we characterized the genome-wide impact of large scale deletions on four divergent Arabidopsis natural strains, through the integration of Ultra-High Throughput Sequencing data with data from genomic hybridizations on tiling arrays. Analysis of the identified deletions revealed a considerable portion of protein coding genes affected and supported a history of genomic rearrangements shaped by evolution. In the last part of the thesis, we showed that VIP3 gene in Arabidopsis has an evo-lutionary conserved role in the 3' to 5' mRNA degradation machinery, by applying a novel approach for the analysis of mRNA-Seq data from random-primed mRNA. Altogether, this PhD research contains major advancements in the study of natural genomic variation in plants and in the application of computational morphodynamics models for the functional characterization of biological pathways essential for the plant. - Récemment, l'introduction du séquençage de seconde génération et les avancées dans la microscopie confocale ont permis des études à l'échelle des différents systèmes cellulaires pour la caractérisation fonctionnelle de gènes. Le degrés de complexité intrinsèque à ces approches ont requis le développement de méthodologies bioinformatiques et de modèles mathématiques afin d'extraire de la masse de données expérimentale générée, des information biologiques significatives. Ce doctorat présente à la fois des méthodes bioinformatiques originales et des modèles mathématiques pour répondre à certaines questions spécifiques de Biologie Végétale en utilisant la plante Arabidopsis thaliana comme modèle. Premièrement, une analyse qualitative spatio-temporelle de profiles quantitatifs de transcripts et de protéines est utilisée pour montrer le rôle de la protéine BREVIS RADIX (BRX) dans le dialogue entre l'auxine et les cytokinines, des phytohormones, dans la croissance du méristème racinaire. Le noyau de ce travail de thèse est la caractérisation fonctionnelle de l'interaction entre la protéine BRX et la phytohormone auxine dans le méristème de la racine en utilisant des modèles informatiques basés sur des preuves expérimentales. Les hypothèses produites par le modèle ont mené à la découverte d'un schéma différentiel d'endocytose dans le méristème racinaire qui divise la réponse transcriptionnelle à l'auxine par le partitionnement de BRX de la membrane plasmique au noyau de la cellule. Cette information positionnelle crée une réponse transcriptionnelle à l'auxine qui dévie de la réponse canonique à l'auxine et est nécessaire pour soutenir l'expression d'un sous ensemble de gènes répondant à l'auxine et dépendant de BRX pour conduire la croissance du méristème. Dans la seconde partie de cette thèse de doctorat, nous avons caractérisé l'impact sur l'ensemble du génome des délétions à grande échelle sur quatre souches divergentes naturelles d'Arabidopsis, à travers l'intégration du séquençage à ultra-haut-débit avec l'hybridation génomique sur puces ADN. L'analyse des délétions identifiées a révélé qu'une proportion considérable de gènes codant était affectée, supportant l'idée d'un historique de réarrangement génomique modelé durant l'évolution. Dans la dernière partie de cette thèse, nous avons montré que le gène VÏP3 dans Arabidopsis a conservé un rôle évolutif dans la machinerie de dégradation des ARNm dans le sens 3' à 5', en appliquant une nouvelle approche pour l'analyse des données de séquençage d'ARNm issue de transcripts amplifiés aléatoirement. Dans son ensemble, cette recherche de doctorat contient des avancées majeures dans l'étude des variations génomiques naturelles des plantes et dans l'application de modèles morphodynamiques informatiques pour la caractérisation de réseaux biologiques essentiels à la plante. - Le développement des plantes est écrit dans leurs codes génétiques. Pour comprendre comment les plantes sont capables de s'adapter aux changements environnementaux, il est essentiel d'étudier comment leurs gènes gouvernent leur formation. Plus nous essayons de comprendre le fonctionnement d'une plante, plus nous réalisons la complexité des mécanismes biologiques, à tel point que l'utilisation d'outils et de modèles mathématiques devient indispensable. Dans ce travail, avec l'utilisation de la plante modèle Arabidopsis thalicinci nous avons résolu des problèmes biologiques spécifiques à travers le développement et l'application de méthodes informatiques concrètes. Dans un premier temps, nous avons investigué comment le gène BREVIS RADIX (BRX) régule le développement de la racine en contrôlant la réponse à deux hormones : l'auxine et la cytokinine. Nous avons employé une analyse statistique sur des mesures quantitatives de transcripts et de produits de gènes afin de démontrer que BRX joue un rôle antagonisant dans le dialogue entre ces deux hormones. Lorsque ce-dialogue moléculaire est perturbé, la racine primaire voit sa longueur dramatiquement réduite. Pour comprendre comment BRX répond à l'auxine, nous avons développé un modèle informatique basé sur des résultats expérimentaux. Les simulations successives ont mené à la découverte d'un signal positionnel qui contrôle la réponse de la racine à l'auxine par la régulation du mouvement intracellulaire de BRX. Dans la seconde partie de cette thèse, nous avons analysé le génome entier de quatre souches naturelles d'Arabidopsis et nous avons trouvé qu'une grande partie de leurs gènes étaient manquant par rapport à la souche de référence. Ce résultat indique que l'historique des modifications génomiques conduites par l'évolution détermine une disponibilité différentielle des gènes fonctionnels dans ces plantes. Dans la dernière partie de ce travail, nous avons analysé les données du transcriptome de la plante où le gène VIP3 était non fonctionnel. Ceci nous a permis de découvrir le rôle double de VIP3 dans la régulation de l'initiation de la transcription et dans la dégradation des transcripts. Ce rôle double n'avait jusqu'alors été démontrée que chez l'homme. Ce travail de doctorat supporte le développement et l'application de méthodologies informatiques comme outils inestimables pour résoudre la complexité des problèmes biologiques dans la recherche végétale. L'intégration de la biologie végétale et l'informatique est devenue de plus en plus importante pour l'avancée de nos connaissances sur le fonctionnement et le développement des plantes.

Textural classification of land cover using support vector machines: an empirical comparison with parametric, non parametric and hybrid classifiers in the Bolivian Amazon

Land cover classification is a key research field in remote sensing and land change science as thematic maps derived from remotely sensed data have become the basis for analyzing many socio-ecological issues. However, land cover classification remains a difficult task and it is especially challenging in heterogeneous tropical landscapes where nonetheless such maps are of great importance. The present study aims to establish an efficient classification approach to accurately map all broad land cover classes in a large, heterogeneous tropical area of Bolivia, as a basis for further studies (e.g., land cover-land use change). Specifically, we compare the performance of parametric (maximum likelihood), non-parametric (k-nearest neighbour and four different support vector machines - SVM), and hybrid classifiers, using both hard and soft (fuzzy) accuracy assessments. In addition, we test whether the inclusion of a textural index (homogeneity) in the classifications improves their performance. We classified Landsat imagery for two dates corresponding to dry and wet seasons and found that non-parametric, and particularly SVM classifiers, outperformed both parametric and hybrid classifiers. We also found that the use of the homogeneity index along with reflectance bands significantly increased the overall accuracy of all the classifications, but particularly of SVM algorithms. We observed that improvements in producer’s and user’s accuracies through the inclusion of the homogeneity index were different depending on land cover classes. Earlygrowth/degraded forests, pastures, grasslands and savanna were the classes most improved, especially with the SVM radial basis function and SVM sigmoid classifiers, though with both classifiers all land cover classes were mapped with producer’s and user’s accuracies of around 90%. Our approach seems very well suited to accurately map land cover in tropical regions, thus having the potential to contribute to conservation initiatives, climate change mitigation schemes such as REDD+, and rural development policies.

The Body Action and Posture Coding System (BAP): Development and Reliability

Several methods are available for coding body movement in nonverbal behavior research, but there is no consensus on a reliable coding system that can be used for the study of emotion expression. Adopting an integrative approach, we developed a new method, the Body Action and Posture (BAP) coding system, for the time-aligned micro description of body movement on an anatomical level (different articulations of body parts), a form level (direction and orientation of movement), and a functional level (communicative and self-regulatory functions). We applied the system to a new corpus of acted emotion portrayals, examined its comprehensiveness and demonstrated intercoder reliability at three levels: a) occurrence, b) temporal precision and c) segmentation. We discuss issues for further validation and propose some research applications.

Vertebrate conserved non coding DNA regions have a high persistence length and a short persistence time.

BACKGROUND: The comparison of complete genomes has revealed surprisingly large numbers of conserved non-protein-coding (CNC) DNA regions. However, the biological function of CNC remains elusive. CNC differ in two aspects from conserved protein-coding regions. They are not conserved across phylum boundaries, and they do not contain readily detectable sub-domains. Here we characterize the persistence length and time of CNC and conserved protein-coding regions in the vertebrate and insect lineages. RESULTS: The persistence length is the length of a genome region over which a certain level of sequence identity is consistently maintained. The persistence time is the evolutionary period during which a conserved region evolves under the same selective constraints.Our main findings are: (i) Insect genomes contain 1.60 times less conserved information than vertebrates; (ii) Vertebrate CNC have a higher persistence length than conserved coding regions or insect CNC; (iii) CNC have shorter persistence times as compared to conserved coding regions in both lineages. CONCLUSION: Higher persistence length of vertebrate CNC indicates that the conserved information in vertebrates and insects is organized in functional elements of different lengths. These findings might be related to the higher morphological complexity of vertebrates and give clues about the structure of active CNC elements.Shorter persistence time might explain the previously puzzling observations of highly conserved CNC within each phylum, and of a lack of conservation between phyla. It suggests that CNC divergence might be a key factor in vertebrate evolution. Further evolutionary studies will help to relate individual CNC to specific developmental processes.

An mRNA region of the canine distemper virus fusion protein gene lacking AUG codons can promote protein expression.

Canine distemper virus (CDV) produces a glycosylated type I fusion protein (F) with an internal hydrophobic signal sequence beginning around 115 residues downstream of the first AUG used for translation initiation. Cleavage of the signal sequence yields the F0 molecule, which is cleaved into the F1 and F2 subunits. Surprisingly, when all in-frame AUGs located in the first third of the F gene were mutated a protein of the same molecular size as the F0 molecule was still expressed from both the Onderstepoort (OP) and A75/17-CDV F genes. We designated this protein, which is initiated from a non-AUG codon protein Fx. Site-directed mutagenesis allowed to identify codon 85, a GCC codon coding for alanine, as the most likely position from which translation initiation of Fx occurs in OP-CDV. Deletion analysis demonstrated that at least 60 nucleotides upstream of the GCC codon are required for efficient Fx translation. This sequence is GC-rich, suggesting extensive folding. Secondary structure may therefore be important for translation initiation at codon 85.

Characterization of the murine high Km glucose transporter GLUT2 gene and its transcriptional regulation by glucose in a differentiated insulin-secreting cell line.

In pancreatic beta-cells, the high Km glucose transporter GLUT2 catalyzes the first step in glucose-induced insulin secretion by glucose uptake. Expression of the transporter has been reported to be modulated by glucose either at the protein or mRNA levels. In this study we used the differentiated insulinoma cell line INS-1 which expresses high levels of GLUT2 and show that the expression of GLUT2 is regulated by glucose at the transcriptional level. By run-on transcription assays we showed that glucose induced GLUT2 gene transcription 3-4-fold in INS-1 cells which was paralleled by a 1.7-2.3-fold increase in cytoplasmic GLUT2 mRNA levels. To determine whether glucose regulatory sequences were present in the promoter region of GLUT2, we cloned and characterized a 1.4-kilobase region of mouse genomic DNA located 5' of the translation initiation site. By RNase protection assays and primer extension, we determined that multiple transcription initiation sites were present at positions -55, -64, and -115 from the first coding ATG and which were identified in liver, intestine, kidney, and beta-cells mRNAs. Plasmids were constructed with the mouse promoter region linked to the reporter gene chloramphenicol acetyltransferase (CAT), and transiently and stably transfected in the INS-1 cells. Glucose induced a concentration-dependent increase in CAT activity which reached a maximum of 3.6-fold at 20 mM glucose. Similar CAT constructs made of the human GLUT2 promoter region and the CAT gene displayed the same glucose-dependent increase in transcriptional activity when transfected into INS-1 cells. Comparison of the mouse and human promoter regions revealed sequence identity restricted to a few stretches of sequences which suggests that the glucose responsive element(s) may be conserved in these common sequences.

Amniotic fluid insulin-like growth factor binding protein 3 concentration as early indicator of fetal growth restriction.

OBJECTIVE: Insulin-like growth factor-I (IGF-I) is an important regulator of fetal growth and its bioavailability depends on insulin-like growth factor binding proteins (IGFBPs). Genes coding for IGF-I and IGFBP3 are polymorphic. We hypothesized that either amniotic fluid protein concentration at the beginning of the second trimester or genotype of one of these two genes could be predictive of abnormal fetal growth. STUDY DESIGN: Amniotic fluid samples (14-18 weeks of pregnancy) from 123 patients with appropriate for gestational age (AGA) fetuses, 39 patients with small for gestational age (SGA) fetuses and 34 patients with large for gestational age (LGA) were analyzed. Protein concentrations were evaluated by ELISA and gene polymorphisms by PCR. RESULTS: Amniotic fluid IGFBP3 concentrations were significantly higher in SGA compared to AGA group (P=0.030), and this was even more significant when adjusted to gestational age at the time of amniocentesis and other covariates (ANCOVA analysis: P=0.009). Genotypic distribution of IGF-I variable number of tandem repeats (VNTR) polymorphism was significantly different in SGA compared to AGA group (P=0.029). 19CA/20CA genotype frequency was threefold decreased in SGA compared to AGA group and the risk of SGA occurrence of this genotype was decreased accordingly: OR=0.289, 95%CI=0.1-0.9, P=0.032. Genotype distribution of IGFBP3(A-202C) polymorphism was similar in all three groups. CONCLUSIONS: High IGFBP3 concentrations in amniotic fluid at the beginning of the second trimester are associated with increased risks of SGA while 19CA/20CA genotype at IGF-I VNTR polymorphism is associated with reduced risks of SGA. Neither IGFBP3 concentrations, nor IGF-I/IGFBP3 polymorphisms are associated with modified risks of LGA.

Cloning and characterization of SmZF1, a gene encoding a Schistosoma mansoni zinc finger protein

The zinc finger motifs (Cys2His2) are found in several proteins playing a role in the regulation of transcripton. SmZF1, a Schistosoma mansoni gene encoding a zinc finger protein was initially isolated from an adult worm cDNA library, as a partial cDNA. The full sequence of the gene was obtained by subcloning and sequencing cDNA and genomic fragments. The collated gene sequence is 2181 nt and the complete cDNA sequence is 705 bp containing the full open reading frame of the gene. Analysis of the genome sequence revealed the presence of three introns interrupting the coding region. The open reading frame theoretically encodes a protein of 164 amino acids, with a calculated molecular mass of 18,667Da. The predicted protein contains three zinc finger motifs, usually present in transcription regulatory proteins. PCR amplification with specific primers for the gene allowed for the detection of the target in egg, cercariae, schistosomulum and adult worm cDNA libraries indicating the expression of the mRNA in these life cycle stages of S. mansoni. This pattern of expression suggests the gene plays a role in vital functions of different life cycle stages of the parasite. Future research will be directed to elucidate the functional role of SmZF1.

A Heminested Polymerase Chain Reaction for the Detection of Brazilian Rabies Isolates from Vampire Bats and Herbivores

A heminested-PCR (hn-PCR) using primers to the nucleoprotein-coding gene in a nested set was evaluated in the detection of Brazilian strains of rabies virus (RV). A representative number of RV nucleoprotein sequences belonging to genotype 1 were aligned. Based on such alignment, primers were directed to highly conserved regions. All 42 clinical samples positive by both fluorescent antibody and mouse inoculation tests were also positive by the hn-PCR. Brain tissue that had been left to decompose, obtained from an experimentally inoculated mouse was tested by hn-PCR and yielded positive results. In conclusion, primers designed here were capable of amplifying Brazilian RV isolates obtained from a rural epidemiological cycle.

The 245 kb amplified chromosome of Leishmania (V.) braziliensis contains a biopterin transporter gene

Leishmania (V.) braziliensis M2903 presents a small linear and stable 245 kb chromosome originating from a genomic amplification. Similar amplifications present in other species of Leishmania contain a gene coding for a biopterin transporter. Since Leishmania is auxotrophic for this metabolite, this amplification could result from the need to better capture biotpterin from growth media under specific circumstances. In this paper we show that this gene is also present in L. (V.) braziliensis small chromosome, which shares sequences with other genomic amplifications already described.

Differential requirement for CD4 help in the development of an antigen-specific CD8+ T cell response depending on the route of immunization.

Previous studies in our laboratory have shown that DBA/2 mice injected i.p. with syngeneic P815 tumor cells transfected with the HLA-CW3 gene (P815-CW3) showed a dramatic expansion of activated CD8+CD62L- T cells expressing exclusively the Vbeta10 segment. We have used this model to study the regulatory mechanisms involved in the development of the CW3-specific CD8+ response, with respect to different routes of immunization. Whereas both intradermal (i.d.) and i.p. immunization of DBA/2 mice with P815-CW3 cells led to a strong expansion of CD8+CD62L-Vbeta10+ cells, only the i.d. route allowed this expansion after immunization with P815 cells transfected with a minigene coding for the antigenic epitope CW3 170-179 (P815 miniCW3). Furthermore, depletion of CD4+ T cells in vivo completely abolished the specific response of CD8+CD62L-Vbeta10+ cells and prevented the rejection of P815-CW3 tumor cells injected i.p., whereas it did not affect CD8S+CD62L-Vbeta10+ cell expansion after i.d. immunization with either P815-CW3 or P815 miniCW3. Finally, the CW3-specific CD8+ memory response was identical whether or not CD4+ T cells were depleted during the primary response. Collectively, these results suggest that the CD8+ T cell response to P815-CW3 tumor cells injected i.p. is strictly dependent upon recognition of a helper epitope by CD4+ T cells, whereas no such requirement is observed for i.d. injection.

cAMP prevents the glucose-mediated stimulation of GLUT2 gene transcription in hepatocytes.

Glucose homoeostasis necessitates the presence in the liver of the high Km glucose transporter GLUT2. In hepatocytes, we and others have demonstrated that glucose stimulates GLUT2 gene expression in vivo and in vitro. This effect is transcriptionally regulated and requires glucose metabolism within the hepatocytes. In this report, we further characterized the cis-elements of the murine GLUT2 promoter, which confers glucose responsiveness on a reporter gene coding the chloramphenicol acetyl transferase (CAT) gene. 5'-Deletions of the murine GLUT2 promoter linked to the CAT reporter gene were transfected into a GLUT2 expressing hepatoma cell line (mhAT3F) and into primary cultured rat hepatocytes, and subsequently incubated at low and high glucose concentrations. Glucose stimulates gene transcription in a manner similar to that observed for the endogenous GLUT2 mRNA in both cell types; the -1308 to -212 bp region of the promoter contains the glucose-responsive elements. Furthermore, the -1308 to -338 bp region of the promoter contains repressor elements when tested in an heterologous thymidine kinase promoter. The glucose-induced GLUT2 mRNA accumulation was decreased by dibutyryl-cAMP both in mhAT3F cells and in primary hepatocytes. A putative cAMP-responsive element (CRE) is localized at the -1074/-1068 bp region of the promoter. The inhibitory effect of cAMP on GLUT2 gene expression was observed in hepatocytes transfected with constructs containing this CRE (-1308/+49 bp fragment), as well as with constructs not containing the consensus CRE (-312/+49 bp fragment). This suggests that the inhibitory effect of cAMP is not mediated by the putative binding site located in the repressor fragment of the GLUT2 promoter. Taken together, these data demonstrate that the elements conferring glucose and cAMP responsiveness on the GLUT2 gene are located within the -312/+49 region of the promoter.

Identification of small RNAs in Pseudomonas species

Summary: Bacterial small RNAs (sRNAs) are transcripts most of which have regulatory functions. Sequence and secondary structure elements enable numerous sRNAs to interact with mRNAs or with regulatory proteins resulting in diverse regulatory effects on virulence, iron storage, organization of cell envelope proteins or stress response. sRNAs having high affinity for RsmA-like RNA-binding proteins are important for posttranscriptional regulation in various Gram-negative bacteria. In Pseudomonas spp., the GacS/GacA two component system positively controls the production of such sRNAs. They titrate RsmA-like proteins and thus overcome translational repression due to these proteins. As a consequence, secondary metabolites can be produced that are implicated in the biocontrol capacity of P. fluorescens or in the virulence of P. aeruginosa. A genome-wide search carried out in P. aeruginosa PAO1 and in closely related Pseudomonas spp. resulted in the identification of 15 genes coding for sRNAs. Eight of these are novel, the remaining seven have previously been observed. Among them, the 1698 sRNA gene was expressed under GacA control, whereas the transcription of 1887 sRNA gene was transcribed under the control of the anaerobic regulator Anr in an oxygen-limited environment. Overexpression of 1698 sRNA in P. fluorescens strain CHAO did not affect the expression of the GacA-regulated hcnA gene (first gene of the operon coding for HCN synthase), indicating that 1698 sRNA is probably not part of the secondary metabolite regulation pathway. The expression of 1698 sRNA was positively regulated by RpoS in both P. aeruginosa PAO 1 and P. ,fluorescens CHAO and appeared to be modulated temporarily by oxidative stress conditions. However, the effect of 1698 sRNA on oxidative stress survival has not yet been established. Hfq protein interacted with 1698 sRNA in vitro and improved 1698 sRNA expression in vivo in P. aeruginosa. In P. fluorescens, GacA and Hfq were both required for expression of rpoS and GacA showed a positively control on the hfq expression; therefore, at least in this organism, GacA control of 1698 sRNA expression may act indirectly via Hfq and RpoS. Different methods were employed to find abase-pairing target for 1698 sRNA. In a proteomic analysis carried out in P. aeruginosa, positive regulation by 1698 sRNA was observed for Soda, the iron-associated superoxide dismutase, an enzyme involved in oxidative stress resistance. A sequence complementary with 1698 sRNA was predicted to be located in the 5' leader of soda mRNA. However, base-pairing between soda mRNA and 1698 sRNA remains to be proven. In conclusion, this work has revealed eight novel sRNAs and novel functions of two sRNAs in Pseudomonas spp. Résumé Les petits ARNs non-codants (sRNAs) produits par les bactéries sont des transcrits ayant pour la plupart des activités régulatrices importantes. Leurs séquences nucléotidiques ainsi que leurs structures secondaires permettent aux sRNAs d'interagir soit avec des RNA messagers (mRNAs), de sorte à modifier l'expression des protéines pour lesquelles ils codent, soit avec des protéines régulatrices liant des rnRNAs, ce qui a pour effet de modifier l'expression de ces mRNAs. Des sRNAs sont impliqués dans diverses voies de régulation, telles que celles qui régissent la virulence, le stockage du fer, l'organisation des protéines de l'enveloppe bactérienne ou la réponse au stress. Chez les Pseudomonas spp., le système à deux composantes GacS/GacA contrôle la production de métabolites secondaires. Ceux-ci sont engagés dans l'établissement du biocontrôle, chez P. fluorescens, ou. de la virulence, chez P. aeruginosa. La régulation génique dirigée par le système GacS/GacA fait intervenir les sRNAs du type RsmZ, capables de contrecarrer l'action au niveau traductionnel exercée par les protéines régulatrices du type RsmA. Une recherche au niveau du génome a été menée chez P. aeruginosa PAO1 de même que chez des espèces qui lui sont étroitement apparentées, débouchant sur la mise en évidence de 15 gènes codant pour des sRNAs. Parmi ceux-ci, huit ont été découverts pour la première fois et sept confirment des travaux publiés. L'expression du gène du sRNAs 1698 s'avère être régulée par GacA, vraisemblablement de manière indirecte. La transcription du gène du sRNA 1887 montre une dépendance envers Anr, régulateur de l'anaérobiose, et envers une carence en oxygène. La surexpression du sRNA 1698 chez P. fluorescens CHAO n'affecte pas l'expression de hcnA, un gène du régulon GacA, laissant supposer que le sRNA n'intervient pas dans la régulation des métabolites secondaires. Chez P. aeruginosa PAOI et chez P. fluorescens CHAO, RpoS, le facteur sigma du stress, est nécessaire à l'expression du sRNA 1698, et la concentration de ce dernier est modulée par des conditions de stress oxydatif. Toutefois, un effet du sRNA 1698 quant à la survie suite au stress oxydatif n'a pas été établi. Par ailleurs, l'interaction entre le sRNA 1698 et Hfq, la protéine chaperone de RNAs, in vitro ainsi qu'un rôle positif de Hfq pour l'expression du sRNA 1698 in vivo ont été démontrés chez P. aeruginosa. L'induction de l'expression par GacA de rpoS et de hfq a été confirmée chez P. fluorescens CHAO, suggérant que la régulation par GacA du sRNA 1698 pourrait se faire par l'intermédiaire de RpoS et Hfq. Diverses méthodes ont été employées pour identifier un transcrit qui puisse être apparié par le sRNA 1698. Une analyse de protéome chez P. aeruginosa montre que l'expression de Soda, la superoxyde dismutase associée au fer, est positivement régulée par le sRNA 1698. Soda est une enzyme impliquée dans la résistance au stress oxydatif. Une séquence de complémentarité avec le sRNA 1698 a bien été prédite sur le leader 5' du mRNA de soda. Cependant, l'appariement entre le sRNA et son transcrit cible est encore à prouver. En conclusion, ce travail a dévoilé huit nouveaux sRNAs et de nouvelles fonctions pour deux sRNAs chez les Pseudomonas.

Communally breeding Bechstein's bats have a stable social system that is independent from the postglacial history and location of the populations.

Investigating macro-geographical genetic structures of animal populations is crucial to reconstruct population histories and to identify significant units for conservation. This approach may also provide information about the intraspecific flexibility of social systems. We investigated the history and current structure of a large number of populations in the communally breeding Bechstein's bat (Myotis bechsteinii). Our aim was to understand which factors shape the species' social system over a large ecological and geographical range. Using sequence data from one coding and one noncoding mitochondrial DNA region, we identified the Balkan Peninsula as the main and probably only glacial refugium of the species in Europe. Sequence data also suggest the presence of a cryptic taxon in the Caucasus and Anatolia. In a second step, we used seven autosomal and two mitochondrial microsatellite loci to compare population structures inside and outside of the Balkan glacial refugium. Central European and Balkan populations both were more strongly differentiated for mitochondrial DNA than for nuclear DNA, had higher genetic diversities and lower levels of relatedness at swarming (mating) sites than in maternity (breeding) colonies, and showed more differentiation between colonies than between swarming sites. All these suggest that populations are shaped by strong female philopatry, male dispersal, and outbreeding throughout their European range. We conclude that Bechstein's bats have a stable social system that is independent from the postglacial history and location of the populations. Our findings have implications for the understanding of the benefits of sociality in female Bechstein's bats and for the conservation of this endangered species.

BlastR--fast and accurate database searches for non-coding RNAs.

We present and validate BlastR, a method for efficiently and accurately searching non-coding RNAs. Our approach relies on the comparison of di-nucleotides using BlosumR, a new log-odd substitution matrix. In order to use BlosumR for comparison, we recoded RNA sequences into protein-like sequences. We then showed that BlosumR can be used along with the BlastP algorithm in order to search non-coding RNA sequences. Using Rfam as a gold standard, we benchmarked this approach and show BlastR to be more sensitive than BlastN. We also show that BlastR is both faster and more sensitive than BlastP used with a single nucleotide log-odd substitution matrix. BlastR, when used in combination with WU-BlastP, is about 5% more accurate than WU-BlastN and about 50 times slower. The approach shown here is equally effective when combined with the NCBI-Blast package. The software is an open source freeware available from www.tcoffee.org/blastr.html.