Bootstrap and permutation tests of independence for point processes

Motivated by a neuroscience question about synchrony detection in spike train analysis, we deal with the independence testing problem for point processes. We introduce non-parametric test statistics, which are rescaled general $U$-statistics, whose corresponding critical values are constructed from bootstrap and randomization/permutation approaches, making as few assumptions as possible on the underlying distribution of the point processes. We derive general consistency results for the bootstrap and for the permutation w.r.t. to Wasserstein's metric, which induce weak convergence as well as convergence of second order moments. The obtained bootstrap or permutation independence tests are thus proved to be asymptotically of the prescribed size, and to be consistent against any reasonable alternative. A simulation study is performed to illustrate the derived theoretical results, and to compare the performance of our new tests with existing ones in the neuroscientific literature

Genetic Structuration, Demography and Evolutionary History of Mycobacterium tuberculosis LAM9 Sublineage in the Americas as Two Distinct Subpopulations Revealed by Bayesian Analyses.

International audienceTuberculosis (TB) remains broadly present in the Americas despite intense global efforts for its control and elimination. Starting from a large dataset comprising spoligotyping (n = 21183 isolates) and 12-loci MIRU-VNTRs data (n = 4022 isolates) from a total of 31 countries of the Americas (data extracted from the SITVIT2 database), this study aimed to get an overview of lineages circulating in the Americas. A total of 17119 (80.8%) strains belonged to the Euro-American lineage 4, among which the most predominant genotypic family belonged to the Latin American and Mediterranean (LAM) lineage (n = 6386, 30.1% of strains). By combining classical phylogenetic analyses and Bayesian approaches, this study revealed for the first time a clear genetic structuration of LAM9 sublineage into two subpopulations named LAM9C1 and LAM9C2, with distinct genetic characteristics. LAM9C1 was predominant in Chile, Colombia and USA, while LAM9C2 was predominant in Brazil, Dominican Republic, Guadeloupe and French Guiana. Globally, LAM9C2 was characterized by higher allelic richness as compared to LAM9C1 isolates. Moreover, LAM9C2 sublineage appeared to expand close to twenty times more than LAM9C1 and showed older traces of expansion. Interestingly, a significant proportion of LAM9C2 isolates presented typical signature of ancestral LAM-RDRio MIRU-VNTR type (224226153321). Further studies based on Whole Genome Sequencing of LAM strains will provide the needed resolution to decipher the biogeographical structure and evolutionary history of this successful family

Surrogate data methods based on a shuffling of the trials for synchrony detection: the centering issue

International audienceWe investigate several distribution-free dependence detection procedures, all based on a shuffling of the trials, from a statistical point of view. The mathematical justification of such procedures lies in the bootstrap principle and its approximation properties. In particular, we show that such a shuffling has mainly to be done on centered quantities-that is, quantities with zero mean under independence-to construct correct p-values, meaning that the corresponding tests control their false positive (FP) rate. Thanks to this study, we introduce a method, named permutation UE, which consists of a multiple testing procedure based on permutation of experimental trials and delayed coincidence count. Each involved single test of this procedure achieves the prescribed level, so that the corresponding multiple testing procedure controls the false discovery rate (FDR), and this with as few assumptions as possible on the underneath distribution, except independence and identical distribution across trials. The mathematical meaning of this assumption is discussed, and it is in particular argued that it does not mean what is commonly referred in neuroscience to as cross-trials stationarity. Some simulations show, moreover, that permutation UE outperforms the trial-shuffling of Pipa and Grün ( 2003 ) and the MTGAUE method of Tuleau-Malot et al. ( 2014 ) in terms of single levels and FDR, for a comparable amount of false negatives. Application to real data is also provided

Sweet and Sour Ehrlichia: Glycoproteomics and Phosphoproteomics Reveal New Players in Ehrlichia ruminantium Physiology and Pathogenesis

Ehrlichia ruminantium; N-glycoproteins; O-GlcNAcylated proteinsEhrlichia ruminantium; N-glicoproteïnes; Proteïnes O-GlcNAciladesEhrlichia ruminantium; N-glicoproteínas; Proteínas O-GlcNAciladasUnraveling which proteins and post-translational modifications (PTMs) affect bacterial pathogenesis and physiology in diverse environments is a tough challenge. Herein, we used mass spectrometry-based assays to study protein phosphorylation and glycosylation in Ehrlichia ruminantium Gardel virulent (ERGvir) and attenuated (ERGatt) variants and, how they can modulate Ehrlichia biological processes. The characterization of the S/T/Y phosphoproteome revealed that both strains share the same set of phosphoproteins (n = 58), 36% being overexpressed in ERGvir. The percentage of tyrosine phosphorylation is high (23%) and 66% of the identified peptides are multi-phosphorylated. Glycoproteomics revealed a high percentage of glycoproteins (67% in ERGvir) with a subset of glycoproteins being specific to ERGvir (n = 64/371) and ERGatt (n = 36/343). These glycoproteins are involved in key biological processes such as protein, amino-acid and purine biosynthesis, translation, virulence, DNA repair, and replication. Label-free quantitative analysis revealed over-expression in 31 proteins in ERGvir and 8 in ERGatt. While further PNGase digestion confidently localized 2 and 5 N-glycoproteins in ERGvir and ERGatt, respectively, western blotting suggests that many glycoproteins are O-GlcNAcylated. Twenty-three proteins were detected in both the phospho- and glycoproteome, for the two variants. This work represents the first comprehensive assessment of PTMs on Ehrlichia biology, rising interesting questions regarding ER-host interactions. Phosphoproteome characterization demonstrates an increased versatility of ER phosphoproteins to participate in different mechanisms. The high number of glycoproteins and the lack of glycosyltransferases-coding genes highlight ER dependence on the host and/or vector cellular machinery for its own protein glycosylation. Moreover, these glycoproteins could be crucial to interact and respond to changes in ER environment. PTMs crosstalk between of O-GlcNAcylation and phosphorylation could be used as a major cellular signaling mechanism in ER. As little is known about the Ehrlichia proteins/proteome and its signaling biology, the results presented herein provide a useful resource for further hypothesis-driven exploration of Ehrlichia protein regulation by phosphorylation and glycosylation events. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium with the data set identifier PXD012589

Sweet and Sour Ehrlichia: Glycoproteomics and Phosphoproteomics Reveal New Players in Ehrlichia ruminantium Physiology and Pathogenesis

Unraveling which proteins and post-translational modifications (PTMs) affect bacterial pathogenesis and physiology in diverse environments is a tough challenge. Herein, we used mass spectrometry-based assays to study protein phosphorylation and glycosylation in Ehrlichia ruminantium Gardel virulent (ERGvir) and attenuated (ERGatt) variants and, how they can modulate Ehrlichia biological processes. The characterization of the S/T/Y phosphoproteome revealed that both strains share the same set of phosphoproteins (n = 58), 36% being overexpressed in ERGvir. The percentage of tyrosine phosphorylation is high (23%) and 66% of the identified peptides are multi-phosphorylated. Glycoproteomics revealed a high percentage of glycoproteins (67% in ERGvir) with a subset of glycoproteins being specific to ERGvir (n = 64/371) and ERGatt (n = 36/343). These glycoproteins are involved in key biological processes such as protein, amino-acid and purine biosynthesis, translation, virulence, DNA repair, and replication. Label-free quantitative analysis revealed over-expression in 31 proteins in ERGvir and 8 in ERGatt. While further PNGase digestion confidently localized 2 and 5 N-glycoproteins in ERGvir and ERGatt, respectively, western blotting suggests that many glycoproteins are O-GlcNAcylated. Twenty-three proteins were detected in both the phospho- and glycoproteome, for the two variants. This work represents the first comprehensive assessment of PTMs on Ehrlichia biology, rising interesting questions regarding ER–host interactions. Phosphoproteome characterization demonstrates an increased versatility of ER phosphoproteins to participate in different mechanisms. The high number of glycoproteins and the lack of glycosyltransferases-coding genes highlight ER dependence on the host and/or vector cellular machinery for its own protein glycosylation. Moreover, these glycoproteins could be crucial to interact and respond to changes in ER environment. PTMs crosstalk between of O-GlcNAcylation and phosphorylation could be used as a major cellular signaling mechanism in ER. As little is known about the Ehrlichia proteins/proteome and its signaling biology, the results presented herein provide a useful resource for further hypothesis-driven exploration of Ehrlichia protein regulation by phosphorylation and glycosylation events. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium with the data set identifier PXD012589

Two tales: Worldwide distribution of Central Asian (CAS) versus ancestral East-African Indian (EAI) lineages of Mycobacterium tuberculosis underlines a remarkable cleavage for phylogeographical, epidemiological and demographical characteristics.

The East African Indian (EAI) and Central Asian (CAS) lineages of Mycobacterium tuberculosis complex (MTBC) mainly infect tuberculosis (TB) patients in the eastern hemisphere which contains many of the 22 high TB burden countries including China and India. We investigated if phylogeographical, epidemiological and demographical characteristics for these 2 lineages differed in SITVIT2 database. Genotyping results and associated data (age, sex, HIV serology, drug resistance) on EAI and CAS lineages (n = 10,974 strains) were extracted. Phylogenetic and Bayesian, and other statistical analyses were used to compare isolates. The male/female sex ratio was 907/433 (2.09) for the EAI group vs. 881/544 (1.62) for CAS (p-value<0.002). The proportion of younger patients aged 0-20 yrs. with CAS lineage was significantly higher than for EAI lineage (18.07% vs. 10.85%, p-value<0.0001). The proportion of multidrug resistant and extensively drug resistant TB among CAS group (30.63% and 1.03%, respectively) was significantly higher than in the EAI group (12.14% and 0.29%, respectively; p-value<0.0001). Lastly, the proportion of HIV+ patients was 20.34% among the EAI group vs. 3.46% in the CAS group (p-value<0.0001). This remarkable split observed between various parameters for these 2 lineages was further corroborated by their geographic distribution profile (EAI being predominantly found in Eastern-Coast of Africa, South-India and Southeast Asia, while CAS was predominantly found in Afghanistan, Pakistan, North India, Nepal, Middle-east, Libya, Sudan, Ethiopia, Kenya and Tanzania). Some geo-specificities were highlighted. This study demonstrated a remarkable cleavage for aforementioned characteristics of EAI and CAS lineages, showing a North-South divide along the tropic of cancer in Eastern hemisphere-mainly in Asia, and partly prolonged along the horn of Africa. Such studies would be helpful to better comprehend prevailing TB epidemic in context of its historical spread and evolutionary features, and provide clues to better treatment and patient-care in countries and regions concerned by these lineages

Plan d'échantillonnage national de la surveillance sanitaire microbiologique des zones de production de coquillages (REMI)

Le plan d’échantillonnage national du REMI est un document d’application de la procédure nationale de la surveillance REMI. Il définit, pour chaque zone classée, les lieux de prélèvement identifiés par un nom et un code d’identification, les espèces de coquillages et les fréquences à respecter et rappelle le classement en vigueur

Phylogenetic analyses of Mycobacterium ulcerans from French Guiana using combination of core and accessory genomes

International audienceIn South America, French Guiana is the only Buruli ulcer endemic area, even though the incidence remains low over the years, with a total of 245 cases reported between 1969 and 2013. In a previous study, we reported evidence for uncommon heterogeneity in M. ulcerans (usually characterized by a very low level of genetic variability) and we identified two main contrasted sublineages. Nevertheless, this study was based on classical molecular typing approaches which can sometimes lead to discrepancies in Mycobacterium classification. To overcome these issues, and to move forward in the exploration of M. ulcerans genetic variability, we performed cgMLST (Core Genome Multi Locus Sequence Typing) using whole genome sequences of 25 M. ulcerans strains originated from French Guiana (n=5) and from other countries in the world were Buruli is endemic (n=20). Maximum likelihood phylogeny revealed a total of 5 lineages clearly delineated. As expected lineage L3 previously described by Doig et al. (2012) comprised all isolates from Africa and from Papua New Guinea and, lineage L2 includes Asian isolates from China and Japan. Concerning the lineage L1, Doig et al. initially proposed a clusterization of the M. ulcerans strains Mu_1G897 isolated in French Guiana with MPM isolates (Mycolactone Producing Mycobacteria: M. liflandii, M. pseudoshotsii, M. shinshuense), known to be pathogenic to ectotherms fishes and frogs. In our study, as we included 5 additional M. ulcerans samples from French Guiana that clearly clustered together, we were able to reveal for the first time an independent phylogenetic group called lineage L1.2, while MPM isolates were grouped into L1.1. Finally, the Mexican isolate and the strain Mu_FG102 recently identified in French Guiana constitute the lineage L4; this lineage was also characterized for the first time in this work using WGS