Plateforme de génomique, IPMC UMR7275
660 Route des Lucioles, SOPHIA ANTIPOLIS, 06560 VALBONNE
tél: 04-93-95-77-77, fax: 04-93-95-77-08

Expertise

La plateforme de génomique fonctionnelle de Nice Sophia Antipolis existe depuis 1999. Initialement orientée vers la conception, la fabrication et l'analyse de puces à ADN, elle a contribué à ouvrir cette nouvelle technologie à une large communauté, mettant à cette occasion en place un système d'information performant (Mediante), capable de gérer de grandes masses de données, et fonctionnant en production depuis plus de 10 ans.

Tout en fournissant encore aujourd'hui un service d'analyse de puces à ADN s'appuyant sur la technologie développée par Agilent, son activité s'est principalement réorientée vers des services de séquencage à haut-débit (Illumina NextSeq500), offrant dans ce contexte de nombreux types d'analyses des acides nucléiques, et une capacité pour analyser de grandes collections d'échantillons, y compris au niveau de la cellule unique. L'activité de routine concerne des applications comme le RNA-seq, le smallRNA-seq, le CHiP-seq, le CLIP-seq, le reséquencage, mais des projets spécifiques peuvent aussi etre mis en place dans des domaines moins standards, comme le séquencage de novo de génomes, ou certains protocoles particuliers : riboSeq, capSeq,... La plateforme se compose de 4 ingénieurs wet lab et de 4 bio-informaticiens.

Equipements

  1. Pré-séquencage : Nanodrop, Bioanalyzer, Qubit, CovarisS2, Ion Chef, NeoPrep, Blue pippin
  2. Analyse Single Cell : 10x Genomics Chromium, Fluidigm C1, Fluidigm Biomark
  3. Séquencage : NextSeq500 Illumina, MinION et PromethION Oxford Nanopore Technology, Chromium 10X Genomics
  4. Puces à ADN : High-Resolution Microarray Scanner Agilent, Station Affymetrix


Les résultats sont stockés automatiquement sur le portail d'informations de la plateforme Mediante. Cela concerne notamment les fichiers .BAM d'alignement, les fichiers .BW de couverture et l'ensemble des fichiers de l'analyse secondaire et des analyses statistiques conduites en partenariat avec le collaborateur. Sur demande l'ensemble des données brutes sont également mises à disposition et une aide est fournit pour la soumission des données vers la base de données publiques GEO (Gene Expression Omnibus).

Related publications




Lebrigand Kevin

  lebrigand@ipmc.cnrs.fr
 04 93 95 77 92
 660 route des lucioles 06560 Valbonne - Sophia-Antipolis

63 publications found

1. Detecting subtle transcriptomic perturbations induced by lncRNAs knock-down in single-cell CRISPRi screening using a new sparse supervised autoencoder neural network., Front Bioinform. 2024 Mar 4;4:1340339. doi: 10.3389/fbinf.2024.1340339. eCollection 2024. (Pubmed: 38501112)
Truchi M, Lacoux C, Gille C, Fassy J, Magnone V, Lopes Goncalves R, Girard-Riboulleau C, Manosalva-Pena I, Gautier-Isola M, Lebrigand K, Barbry P, Spicuglia S, Vassaux G, Rezzonico R, Barlaud M, Mari B

Single-cell CRISPR-based transcriptome screens are potent genetic tools for concomitantly assessing the expression profiles of cells targeted by a set of guides RNA (gRNA), and inferring target gene functions from the observed perturbations. However, due to various limitations, this approach lacks sensitivity in detecting weak perturbations and is essentially reliable when studying master regulators such as transcription factors. To overcome the challenge of detecting subtle gRNA induced transcriptomic perturbations and classifying the most responsive cells, we developed a new supervised autoencoder neural network method. Our Sparse supervised autoencoder (SSAE) neural network provides selection of both relevant features (genes) and actual perturbed cells. We applied this method on an in-house single-cell CRISPR-interference-based (CRISPRi) transcriptome screening (CROP-Seq) focusing on a subset of long non-coding RNAs (lncRNAs) regulated by hypoxia, a condition that promote tumor aggressiveness and drug resistance, in the context of lung adenocarcinoma (LUAD). The CROP-seq library of validated gRNA against a subset of lncRNAs and, as positive controls, HIF1A and HIF2A, the 2 main transcription factors of the hypoxic response, was transduced in A549 LUAD cells cultured in normoxia or exposed to hypoxic conditions during 3, 6 or 24 h. We first validated the SSAE approach on HIF1A and HIF2 by confirming the specific effect of their knock-down during the temporal switch of the hypoxic response. Next, the SSAE method was able to detect stable short hypoxia-dependent transcriptomic signatures induced by the knock-down of some lncRNAs candidates, outperforming previously published machine learning approaches. This proof of concept demonstrates the relevance of the SSAE approach for deciphering weak perturbations in single-cell transcriptomic data readout as part of CRISPR-based screening.


2. Combination of CRISPR-Cas9-RNP and Single-Cell RNAseq to Identify Cell State-Specific FOXJ1 Functions in the Human Airway Epithelium. , Methods Mol Biol. 2024;2725:1-25. doi: 10.1007/978-1-0716-3507-0_1. (Pubmed: 37856015)
Zaragosi LE, Gouleau A, Delin M, Lebrigand K, Arguel MJ, Girard-Riboulleau C, Rios G, Redman E, Plaisant M, Waldmann R, Magnone V, Marcet B, Barbry P, Ponzio G

The study of the airway epithelium in vitro is routinely performed using air-liquid culture (ALI) models from nasal or bronchial basal cells. These 3D experimental models allow to follow the regeneration steps of fully differentiated mucociliary epithelium and to study gene function by performing gene invalidation. Recent progress made with CRISPR-based techniques has overcome the experimental difficulty of this approach, by a direct transfection of ribonucleoprotein complexes combining a mix of synthetic small guide RNAs (sgRNAs) and recombinant Cas9. The approach shows more than 95% efficiency and does not require any selection step. A limitation of this approach is that it generates cell populations that contain heterogeneous deletions, which makes the evaluation of invalidation efficiency difficult. We have successfully used Flongle sequencing (Nanopore) to quantify the number of distinct deletions. We describe the use of CRISPR-Cas9 RNP in combination with single-cell RNA sequencing to functionally characterize the impact of gene invalidation in ALI cultures. The complex ecosystem of the airway epithelium, composed of many cell types, makes single-cell approaches particularly relevant to study cell type, or cell state-specific events. This protocol describes the invalidation of FOXJ1 in ALI cultures through the following steps: (1) Establishment of basal cell cultures from nasal turbinates, (2) CRISPR-Cas9 RNP invalidation of FOXJ1, (3) Quantification of FOXJ1 invalidation efficiency by Nanopore sequencing, (4) Dissociation of ALI cultures and single-cell RNAseq, (5) Analysis of single-cell RNAseq data from FOXJ1-invalidated cells.We confirm here that FOXJ1 invalidation impairs the final differentiation step of multiciliated cells and provides a framework to explore other gene functions.


3. Cholesterol efflux pathways hinder KRAS-driven lung tumor progenitor cell expansion. , Cell Stem Cell. 2023 Jun 1;30(6):800-817.e9. doi: 10.1016/j.stem.2023.05.005. (Pubmed: 37267915)
Guilbaud E, Barouillet T, Ilie M, Borowczyk C, Ivanov S, Sarrazy V, Vaillant N, Ayrault M, Castiglione A, Rignol G, Brest P, Bazioti V, Zaitsev K, Lebrigand K, Dussaud S, Magnone V, Bertolotto C, Marchetti S, Irondelle M, Goldberg I, Huby T, Westerterp M, Gautier EL, Mari B, Barbry P, Hofman P, Yvan-Charvet L

Cholesterol efflux pathways could be exploited in tumor biology to unravel cancer vulnerabilities. A mouse model of lung-tumor-bearing KRASG12D mutation with specific disruption of cholesterol efflux pathways in epithelial progenitor cells promoted tumor growth. Defective cholesterol efflux in epithelial progenitor cells governed their transcriptional landscape to support their expansion and create a pro-tolerogenic tumor microenvironment (TME). Overexpression of the apolipoprotein A-I, to raise HDL levels, protected these mice from tumor development and dire pathologic consequences. Mechanistically, HDL blunted a positive feedback loop between growth factor signaling pathways and cholesterol efflux pathways that cancer cells hijack to expand. Cholesterol removal therapy with cyclodextrin reduced tumor burden in progressing tumor by suppressing the proliferation and expansion of epithelial progenitor cells of tumor origin. Local and systemic perturbations of cholesterol efflux pathways were confirmed in human lung adenocarcinoma (LUAD). Our results position cholesterol removal therapy as a putative metabolic target in lung cancer progenitor cells.


4. The MIR34B/C genomic region contains multiple potential regulators of multiciliogenesis. , FEBS Lett. 2023 Jun;597(12):1623-1637. doi: 10.1002/1873-3468.14630. Epub 2023 May 8. (Pubmed: 37102425)
Cavard A, Redman E, Mercey O, Abelanet S, Plaisant M, Arguel MJ, Magnone V, Ruiz García S, Rios G, Deprez M, Lebrigand K, Ponzio G, Caballero I, Barbry P, Zaragosi LE, Marcet B

The MIR449 genomic locus encompasses several regulators of multiciliated cell (MCC) formation (multiciliogenesis). The miR-449 homologs miR-34b/c represent additional regulators of multiciliogenesis that are transcribed from another locus. Here, we characterized the expression of BTG4, LAYN, and HOATZ, located in the MIR34B/C locus using single-cell RNA-seq and super-resolution microscopy from human, mouse, or pig multiciliogenesis models. BTG4, LAYN, and HOATZ transcripts were expressed in both precursors and mature MCCs. The Layilin/LAYN protein was absent from primary cilia, but it was expressed in apical membrane regions or throughout motile cilia. LAYN silencing altered apical actin cap formation and multiciliogenesis. HOATZ protein was detected in primary cilia or throughout motile cilia. Altogether, our data suggest that the MIR34B/C locus may gather potential actors of multiciliogenesis.


5. The spatial landscape of gene expression isoforms in tissue sections. , Nucleic Acids Res. 2023 May 8;51(8):e47. doi: 10.1093/nar/gkad169. (Pubmed: 36928528)
Lebrigand K, Bergenstråhle J, Thrane K, Mollbrink A, Meletis K, Barbry P, Waldmann R, Lundeberg J

In situ capturing technologies add tissue context to gene expression data, with the potential of providing a greater understanding of complex biological systems. However, splicing variants and full-length sequence heterogeneity cannot be characterized at spatial resolution with current transcriptome profiling methods. To that end, we introduce spatial isoform transcriptomics (SiT), an explorative method for characterizing spatial isoform variation and sequence heterogeneity using long-read sequencing. We show in mouse brain how SiT can be used to profile isoform expression and sequence heterogeneity in different areas of the tissue. SiT reveals regional isoform switching of Plp1 gene between different layers of the olfactory bulb, and the use of external single-cell data allows the nomination of cell types expressing each isoform. Furthermore, SiT identifies differential isoform usage for several major genes implicated in brain function (Snap25, Bin1, Gnas) that are independently validated by in situ sequencing. SiT also provides for the first time an in-depth A-to-I RNA editing map of the adult mouse brain. Data exploration can be performed through an online resource (https://www.isomics.eu), where isoform expression and RNA editing can be visualized in a spatial context.


6. GLI1+ cells are a source of repair-supportive mesenchymal cells (RSMCs) during airway epithelial regeneration , Cell Mol Life Sci. 2022 Nov 5;79(11):581. doi: 10.1007/s00018-022-04599-2. (Pubmed: 36333491)
Chu X, Lingampally A, Moiseenko A, Kheirollahi V, Vazquez-Armendariz AI, Koepke J, Khadim A, Kiliaris G, Shahriari Felordi M, Zabihi M, Shalashova I, Alexopoulos I, Günther S, Lebrigand K, Truchi M, Günther A, Braun T, Mari B, Samakovlis C, Li X, Seeger W, Herold S, Zhang JS, Bellusci S, El Agha E

Repair-supportive mesenchymal cells (RSMCs) have been recently reported in the context of naphthalene (NA)-induced airway injury and regeneration. These cells transiently express smooth muscle actin (Acta2) and are enriched with platelet-derived growth factor receptor alpha (Pdgfra) and fibroblast growth factor 10 (Fgf10) expression. Genetic deletion of Ctnnb1 (gene coding for beta catenin) or Fgf10 in these cells using the Acta2-Cre-ERT2 driver line after injury (defined as NA-Tam condition; Tam refers to tamoxifen) led to impaired repair of the airway epithelium. In this study, we demonstrate that RSMCs are mostly captured using the Acta2-Cre-ERT2 driver when labeling occurs after (NA-Tam condition) rather than before injury (Tam-NA condition), and that their expansion occurs mostly between days 3 and 7 following NA treatment. Previous studies have shown that lineage-traced peribronchial GLI1+ cells are transiently amplified after NA injury. Here, we report that Gli1 expression is enriched in RSMCs. Using lineage tracing with Gli1Cre-ERT2 mice combined with genetic inactivation of Fgf10, we show that GLI1+ cells with Fgf10 deletion fail to amplify around the injured airways, thus resulting in impaired airway epithelial repair. Interestingly, Fgf10 expression is not upregulated in GLI1+ cells following NA treatment, suggesting that epithelial repair is mostly due to the increased number of Fgf10-expressing GLI1+ cells. Co-culture of SCGB1A1+ cells with GLI1+ cells isolated from non-injured or injured lungs showed that GLI1+ cells from these two conditions are similarly capable of supporting bronchiolar organoid (or bronchiolosphere) formation. Single-cell RNA sequencing on sorted lineage-labeled cells showed that the RSMC signature resembles that of alveolar fibroblasts. Altogether, our study provides strong evidence for the involvement of mesenchymal progenitors in airway epithelial regeneration and highlights the critical role played by Fgf10-expressing GLI1+ cells in this context.


7. Coupling live-cell imaging and in situ isolation of the same single cell to profile the transient states of predicted drug-tolerant cells., STAR Protoc. 2022 Aug 18;3(3):101600. doi: 10.1016/j.xpro.2022.101600. eCollection 2022 Sep 16. (Pubmed: 36042886)
Bian B, Paquet A, Arguel MJ, Meyer M, Peyre L, Chalabi A, Péré M, Lebrigand K, Waldmann R, Barbry P, Hofman P, Roux J

Cell response variability is a starting point in cancer drug resistance that has been difficult to analyze because the tolerant cell states are short lived. Here, we present fate-seq, an approach to isolate single cells in their transient states of drug sensitivity or tolerance before profiling. The drug response is predicted in live cells, which are laser-captured by microdissection before any drug-induced change can alter their states. This framework enables the identification of the cell-state signatures causing differential cell decisions upon treatment. For complete details on the use and execution of this protocol, please refer to Meyer et al. (2020).


8. Monitoring SARS-CoV-2 variants alterations in Nice neighborhoods by wastewater nanopore sequencing., Lancet Reg Health Eur. 2021 Aug 17:100202. doi: 10.1016/j.lanepe.2021.100202 (Pubmed: 34423327)
Rios G, Lacoux C, Leclercq V, Diamant A, Lebrigand K, Lazuka A, Soyeux E, Lacroix S, Fassy J, Couesnon A, Thiery R, Mari B, Pradier C, Waldmann R, Barbry P

Background: Wastewater surveillance was proposed as an epidemiological tool to define the prevalence and evolution of the SARS-CoV-2 epidemics. However, most implemented SARS-CoV-2 wastewater surveillance projects were based on qPCR measurement of virus titers and did not address the mutational spectrum of SARS-CoV-2 circulating in the population. Methods: We have implemented a nanopore RNA sequencing monitoring system in the city of Nice (France, 550,000 inhabitants). Between October 2020 and March 2021, we monthly analyzed the SARS-CoV-2 variants in 113 wastewater samples collected in the main wastewater treatment plant and 20 neighborhoods. Findings: We initially detected the lineages predominant in Europe at the end of 2020 (B.1.160, B.1.177, B.1.367, B.1.474, and B.1.221). In January, a localized emergence of a variant (Spike:A522S) of the B.1.1.7 lineage occurred in one neighborhood. It rapidly spread and became dominant all over the city. Other variants of concern (B.1.351, P.1) were also detected in some neighborhoods, but at low frequency. Comparison with individual clinical samples collected during the same week showed that wastewater sequencing correctly identified the same lineages as those found in COVID-19 patients. Interpretation: Wastewater sequencing allowed to document the diversity of SARS-CoV-2 sequences within the different neighborhoods of the city of Nice. Our results illustrate how sequencing of sewage samples can be used to track pathogen sequence diversity in the current pandemics and in future infectious disease outbreaks.


9. The FibromiR miR-214-3p Is Upregulated in Duchenne Muscular Dystrophy and Promotes Differentiation of Human Fibro-Adipogenic Muscle Progenitors. , Cells. 2021 Jul 20;10(7):1832. doi: 10.3390/cells10071832 (Pubmed: 34360002)
Arrighi N, Moratal C, Savary G, Fassy J, Nottet N, Pons N, Clément N, Fellah S, Larrue R, Magnone V, Lebrigand K, Pottier N, Dechesne C, Vassaux G, Dani C, Peraldi P, Mari B

Fibrosis is a deleterious invasion of tissues associated with many pathological conditions, such as Duchenne muscular dystrophy (DMD) for which no cure is at present available for its prevention or its treatment. Fibro-adipogenic progenitors (FAPs) are resident cells in the human skeletal muscle and can differentiate into myofibroblasts, which represent the key cell population responsible for fibrosis. In this study, we delineated the pool of microRNAs (miRNAs) that are specifically modulated by TGFβ1 in FAPs versus myogenic progenitors (MPs) by a global miRNome analysis. A subset of candidates, including several "FibromiRs", was found differentially expressed between FAPs and MPs and was also deregulated in DMD versus healthy biopsies. Among them, the expression of the TGFβ1-induced miR-199a~214 cluster was strongly correlated with the fibrotic score in DMD biopsies. Loss-of-function experiments in FAPs indicated that a miR-214-3p inhibitor efficiently blocked expression of fibrogenic markers in both basal conditions and following TGFβ1 stimulation. We found that FGFR1 is a functional target of miR-214-3p, preventing the signaling of the anti-fibrotic FGF2 pathway during FAP fibrogenesis. Overall, our work demonstrates that the « FibromiR » miR-214-3p is a key activator of FAP fibrogenesis by modulating the FGF2/FGFR1/TGFβ axis, opening new avenues for the treatment of DMD.


10. Single-cell RNA sequencing reveals intratumoral heterogeneity in primary uveal melanomas and identifies HES6 as a driver of the metastatic disease., Cell Death Differ. 2021 Jan 18. doi: 10.1038/s41418-020-00730-7 (Pubmed: 33462406)
Pandiani C, Strub T, Nottet N, Cheli Y, Gambi G, Bille K, Husser C, Dalmasso M, Béranger G, Lassalle S, Magnone V, Pédeutour F, Irondelle M, Maschi C, Nahon-Estève S, Martel A, Caujolle JP, Hofman P, LeBrigand K, Davidson I, Baillif S, Barbry P, Ballotti R, Bertolotto C

Intratumor heterogeneity has been recognized in numerous cancers as a major source of metastatic dissemination. In uveal melanomas, the existence and identity of specific subpopulations, their biological function and their contribution to metastasis remain unknown. Here, in multiscale analyses using single-cell RNA sequencing of six different primary uveal melanomas, we uncover an intratumoral heterogeneity at the genomic and transcriptomic level. We identify distinct transcriptional cell states and diverse tumor-associated populations in a subset of the samples. We also decipher a gene regulatory network underlying an invasive and poor prognosis state driven in part by the transcription factor HES6. HES6 heterogenous expression has been validated by RNAscope assays within primary human uveal melanomas, which further unveils the existence of these cells conveying a dismal prognosis in tumors diagnosed with a favorable outcome using bulk analyses. Depletion of HES6 impairs proliferation, migration and metastatic dissemination in vitro and in vivo using the chick chorioallantoic membrane assay, demonstrating the essential role of HES6 in uveal melanomas. Thus, single-cell analysis offers an unprecedented view of primary uveal melanoma heterogeneity, identifies bona fide biomarkers for metastatic cells in the primary tumor, and reveals targetable modules driving growth and metastasis formation. Significantly, our findings demonstrate that HES6 is a valid target to stop uveal melanoma progression.


11. Identification of a Repair-Supportive Mesenchymal Cell Population during Airway Epithelial Regeneration, Cell Rep. 2020 Dec 22;33(12):108549. doi: 10.1016/j.celrep.2020.108549. (Pubmed: 33357434)
Moiseenko A, Vazquez-Armendariz AI, Kheirollahi V, Chu X, Tata A, Rivetti S, Günther S, Lebrigand K, Herold S, Braun T, Mari B, De Langhe S, Kwapiszewska G, Günther A, Chen C, Seeger W, Tata PR, Zhang JS, Bellusci S, El Agha E

Tissue regeneration requires coordinated and dynamic remodeling of stem and progenitor cells and the surrounding niche. Although the plasticity of epithelial cells has been well explored in many tissues, the dynamic changes occurring in niche cells remain elusive. Here, we show that, during lung repair after naphthalene injury, a population of PDGFRα+ cells emerges in the non-cartilaginous conducting airway niche, which is normally populated by airway smooth muscle cells (ASMCs). This cell population, which we term "repair-supportive mesenchymal cells" (RSMCs), is distinct from conventional ASMCs, which have previously been shown to contribute to epithelial repair. Gene expression analysis on sorted lineage-labeled cells shows that RSMCs express low levels of ASMC markers, but high levels of the pro-regenerative marker Fgf10. Organoid co-cultures demonstrate an enhanced ability for RSMCs in supporting club-cell growth. Our study highlights the dynamics of mesenchymal cells in the airway niche and has implications for chronic airway-injury-associated diseases.


12. Agonist-induced functional analysis and cell sorting associated with single-cell transcriptomics characterizes cell subtypes in normal and pathological brain, Genome Res. 2020 Nov;30(11):1633-1642. doi: 10.1101/gr.262717.120. Epub 2020 Sep 24. (Pubmed: 32973039)
Castagnola S, Cazareth J, Lebrigand K, Jarjat M, Magnone V, Delhaye S, Brau F, Bardoni B, Maurin T

To gain better insight into the dynamic interaction between cells and their environment, we developed the agonist-induced functional analysis and cell sorting (aiFACS) technique, which allows the simultaneous recording and sorting of cells in real-time according to their immediate and individual response to a stimulus. By modulating the aiFACS selection parameters, testing different developmental times, using various stimuli, and multiplying the analysis of readouts, it is possible to analyze cell populations of any normal or pathological tissue. The association of aiFACS with single-cell transcriptomics allows the construction of functional tissue cartography based on specific pharmacological responses of cells. As a proof of concept, we used aiFACS on the dissociated mouse brain, a highly heterogeneous tissue, enriching it in interneurons by stimulation with KCl or with AMPA, an agonist of the glutamate receptors, followed by sorting based on calcium levels. After AMPA stimulus, single-cell transcriptomics of these aiFACS-selected interneurons resulted in a nine-cluster classification. Furthermore, we used aiFACS on interneurons derived from the brain of the Fmr1-KO mouse, a rodent model of fragile X syndrome. We showed that these interneurons manifest a generalized defective response to AMPA compared with wild-type cells, affecting all the analyzed cell clusters at one specific postnatal developmental time.


13. Longevity strategies in response to light in the reef coral Stylophora pistillata, Sci Rep. 2020 Nov 17;10(1):19937. doi: 10.1038/s41598-020-76925-2 (Pubmed: 33203910)
Ottaviani A, Eid R, Zoccola D, Pousse M, Dubal JM, Barajas E, Jamet K, Lebrigand K, Lapébie P, Baudoin C, Giraud-Panis MJ, Rouan A, Beauchef G, Guéré C, Vié K, Barbry P, Tambutté S, Gilson E, Allemand D

Aging is a multifactorial process that results in progressive loss of regenerative capacity and tissue function while simultaneously favoring the development of a large array of age-related diseases. Evidence suggests that the accumulation of senescent cells in tissue promotes both normal and pathological aging. Oxic stress is a key driver of cellular senescence. Because symbiotic long-lived reef corals experience daily hyperoxic and hypoxic transitions, we hypothesized that these long-lived animals have developed specific longevity strategies in response to light. We analyzed transcriptome variation in the reef coral Stylophora pistillata during the day-night cycle and revealed a signature of the FoxO longevity pathway. We confirmed this pathway by immunofluorescence using antibodies against coral FoxO to demonstrate its nuclear translocation. Through qPCR analysis of nycthemeral variations of candidate genes under different light regimens, we found that, among genes that were specifically up- or downregulated upon exposure to light, human orthologs of two "light-up" genes (HEY1 and LONF3) exhibited anti-senescence properties in primary human fibroblasts. Therefore, these genes are interesting candidates for counteracting skin aging. We propose a large screen for other light-up genes and an investigation of the biological response of reef corals to light (e.g., metabolic switching) to elucidate these processes and identify effective interventions for promoting healthy aging in humans.


14. Profiling the Non-genetic Origins of Cancer Drug Resistance with a Single-Cell Functional Genomics Approach Using Predictive Cell Dynamics, Cell Syst. 2020 Oct 21;11(4):367-374.e5.doi: 10.1016/j.cels.2020.08.019. Epub 2020 Sep 23 (Pubmed: 33099406)
Meyer M, Paquet A, Arguel MJ, Peyre L, Gomes-Pereira LC, Lebrigand K, Mograbi B, Brest P, Waldmann R, Barbry P, Hofman P, Roux J

Non-genetic heterogeneity observed in clonal cell populations is an immediate cause of drug resistance that remains challenging to profile because of its transient nature. Here, we coupled three single-cell technologies to link the predicted drug response of a cell to its own genome-wide transcriptomic profile. As a proof of principle, we analyzed the response to tumor-necrosis-factor-related apoptosis-inducing ligand (TRAIL) in HeLa cells to demonstrate that cell dynamics can discriminate the transient transcriptional states at the origin of cell decisions such as sensitivity and resistance. Our same-cell approach, named fate-seq, can reveal the molecular factors regulating the efficacy of a drug in clonal cells, providing therapeutic targets of non-genetic drug resistance otherwise confounded in gene expression noise. A record of this paper's transparent peer review process is included in the Supplemental Information.


15. Identification of oncolytic vaccinia restriction factors in canine high-grade mammary tumor cells using single-cell transcriptomics, PLoS Pathog. 2020 Oct 19;16(10):e1008660. doi: 10.1371/journal.ppat.1008660. eCollection 2020 Oct. (Pubmed: 33075093)
Cambien B, Lebrigand K, Baeri A, Nottet N, Compin C, Lamit A, Ferraris O, Peyrefitte CN, Magnone V, Henriques J, Zaragosi LE, Giorgetti-Peraldi S, Bost F, Gautier-Isola M, Rezzonico R, Barbry P, Barthel R, Mari B, Vassaux G

Mammary carcinoma, including triple-negative breast carcinomas (TNBC) are tumor-types for which human and canine pathologies are closely related at the molecular level. The efficacy of an oncolytic vaccinia virus (VV) was compared in low-passage primary carcinoma cells from TNBC versus non-TNBC. Non-TNBC cells were 28 fold more sensitive to VV than TNBC cells in which VV replication is impaired. Single-cell RNA-seq performed on two different TNBC cell samples, infected or not with VV, highlighted three distinct populations: naïve cells, bystander cells, defined as cells exposed to the virus but not infected and infected cells. The transcriptomes of these three populations showed striking variations in the modulation of pathways regulated by cytokines and growth factors. We hypothesized that the pool of genes expressed in the bystander populations was enriched in antiviral genes. Bioinformatic analysis suggested that the reduced activity of the virus was associated with a higher mesenchymal status of the cells. In addition, we demonstrated experimentally that high expression of one gene, DDIT4, is detrimental to VV production. Considering that DDIT4 is associated with a poor prognosis in various cancers including TNBC, our data highlight DDIT4 as a candidate resistance marker for oncolytic poxvirus therapy. This information could be used to design new generations of oncolytic poxviruses. Beyond the field of gene therapy, this study demonstrates that single-cell transcriptomics can be used to identify cellular factors influencing viral replication.


16. Evidence of early increased sialylation of airway mucins and defective mucociliary clearance in CFTR-deficient piglets, J Cyst Fibros. 2020 Sep 23:S1569-1993(20)30868-7. doi:10.1016/j.jcf.2020.09.009. (Pubmed: 32978064)
Caballero I, Ringot-Destrez B, Si-Tahar M, Barbry P, Guillon A, Lantier I, Berri M, Chevaleyre C, Fleurot I, Barc C, Ramphal R, Pons N, Paquet A, Lebrigand K, Baron C, Bähr A, Klymiuk N, Léonard R, Robbe-Masselot C

Background: Bacterial colonization in cystic fibrosis (CF) lungs has been directly associated to the loss of CFTR function, and/or secondarily linked to repetitive cycles of chronic inflammation/infection. We hypothesized that altered molecular properties of mucins could contribute to this process. Methods: Newborn CFTR+/+ and CFTR-/- were sacrificed before and 6 h after inoculation with luminescent Pseudomonas aeruginosa into the tracheal carina. Tracheal mucosa and the bronchoalveolar lavage (BAL) fluid were collected to determine the level of mucin O-glycosylation, bacteria binding to mucins and the airways transcriptome. Disturbances in mucociliary transport were determined by ex-vivo imaging of luminescent Pseudomonas aeruginosa. Results: We provide evidence of an increased sialylation of CF airway mucins and impaired mucociliary transport that occur before the onset of inflammation. Hypersialylation of mucins was reproduced on tracheal explants from non CF animals treated with GlyH101, an inhibitor of CFTR channel activity, indicating a causal relationship between the absence of CFTR expression and the sialylation of mucins. This increased sialylation was correlated to an increased adherence of P. aeruginosa to mucins. In vivo infection of newborn CF piglets by live luminescent P. aeruginosa demonstrated an impairment of mucociliary transport of this bacterium, with no evidence of pre-existing inflammation. Conclusions: Our results document for the first time in a well-defined CF animal model modifications that affect the O-glycan chains of mucins. These alterations precede infection and inflammation of airway tissues, and provide a favorable context for microbial development in CF lung that hallmarks this disease.


17. High Dimensional Single-Cell Analysis Reveals iNKT Cell Developmental Trajectories and Effector Fate Decision , Cell Rep. 2020 Sep 8;32(10):108116. doi: 10.1016/j.celrep.2020.108116. (Pubmed: 32905761)
Baranek T, Lebrigand K, de Amat Herbozo C, Gonzalez L, Bogard G, Dietrich C, Magnone V, Boisseau C, Jouan Y, Trottein F, Si-Tahar M, Leite-de-Moraes M, Mallevaey T, Paget C

CD1d-restricted invariant Natural Killer T (iNKT) cells represent a unique class of T lymphocytes endowed with potent regulatory and effector immune functions. Although these functions are acquired during thymic ontogeny, the sequence of events that gives rise to discrete effector subsets remains unclear. Using an unbiased single-cell transcriptomic analysis combined with functional assays, we reveal an unappreciated diversity among thymic iNKT cells, especially among iNKT1 cells. Mathematical modeling and biological methods unravel a developmental map whereby iNKT2 cells constitute a transient branching point toward the generation of iNKT1 and iNKT17 cells, which reconciles the two previously proposed models. In addition, we identify the transcription co-factor Four-and-a-half LIM domains protein 2 (FHL2) as a critical cell-intrinsic regulator of iNKT1 specification. Thus, these data illustrate the changing transcriptional network that guides iNKT cell effector fate.


18. High throughput error corrected Nanopore single cell transcriptome sequencing, Nat Commun. 2020 Aug 12;11(1):4025. doi: 10.1038/s41467-020-17800-6. (Pubmed: 32788667)
Lebrigand K, Magnone V, Barbry P, Waldmann R

Droplet-based high throughput single cell sequencing techniques tremendously advanced our insight into cell-to-cell heterogeneity. However, those approaches only allow analysis of one extremity of the transcript after short read sequencing. In consequence, information on splicing and sequence heterogeneity is lost. To overcome this limitation, several approaches that use long-read sequencing were introduced recently. Yet, those techniques are limited by low sequencing depth and/or lacking or inaccurate assignment of unique molecular identifiers (UMIs), which are critical for elimination of PCR bias and artifacts. We introduce ScNaUmi-seq, an approach that combines the high throughput of Oxford Nanopore sequencing with an accurate cell barcode and UMI assignment strategy. UMI guided error correction allows to generate high accuracy full length sequence information with the 10x Genomics single cell isolation system at high sequencing depths. We analyzed transcript isoform diversity in embryonic mouse brain and show that ScNaUmi-seq allows defining splicing and SNVs (RNA editing) at a single cell level.


19. A Single-cell Atlas of the Human Healthy Airways , Am J Respir Crit Care Med. 2020 Jul 29. doi: 10.1164/rccm.201911-2199OC (Pubmed: 32726565)
Deprez M, Zaragosi LE, Truchi M, Becavin C, Ruiz García S, Arguel MJ, Plaisant M, Magnone V, Lebrigand K, Abelanet S, Brau F, Paquet A, Pe'er D, Marquette CH, Leroy S, Barbry P

Rationale: The respiratory tract constitutes an elaborated line of defense that is based on a unique cellular ecosystem. Single-cell profiling methods enable the investigation of cell population distributions and transcriptional changes along the airways. Methods: We have explored the cellular heterogeneity of the human airway epithelium in 10 healthy living volunteers by single-cell RNA profiling. 77,969 cells were collected at 35 distinct locations, from the nose to the 12th division of the airway tree. Results: The resulting atlas is composed of a high percentage of epithelial cells (89.1%), but also immune (6.2%) and stromal (4.7%) cells with distinct cellular proportions in different regions of the airways. It reveals differential gene expression between identical cell types (suprabasal, secretory, and multiciliated cells) from the nose (MUC4, PI3, SIX3) and tracheobronchial (SCGB1A1, TFF3) airways. By contrast, cell-type specific gene expression is stable across all tracheobronchial samples. Our atlas improves the description of ionocytes, pulmonary neuro-endocrine (PNEC) and brush cells, and identifies a related population of NREP-positive cells. We also report the association of KRT13 with dividing cells that are reminiscent of previously described mouse "hillock" cells, and with squamous cells expressing SCEL, SPRR1A/B. Conclusions: Robust characterization of a single-cell cohort in healthy airways establishes a valuable resource for future investigations. The precise description of the continuum existing from the nasal epithelium to successive divisions of the airways and the stable gene expression profile of these regions better defines conditions under which relevant tracheobronchial proxies of human respiratory diseases can be developed.


20. Transcriptomic and Ultrastructural Signatures of K +-Induced Aggregation in Phytophthora parasitica Zoospores, Microorganisms. 2020 Jul 7;8(7):1012. doi: 10.3390/microorganisms8071012. (Pubmed: 32645882)
Bassani I, Rancurel C, Pagnotta S, Orange F, Pons N, Lebrigand K, Panabières F, Counillon L, Noblin X, Galiana E

Most pathogenic oomycetes of the genus Phytophthora spread in water films as flagellated zoospores. Zoospores perceive and produce signals attracting other zoospores, resulting in autoaggregation in vitro or biofilm formation on plant surface. The mechanisms underlying intercellular communication and consequent attraction, adhesion and aggregation are largely unknown. In Phytophthora parasitica, the perception of a K+ gradient induces coordinated motion and aggregation. To define cellular and molecular events associated with oomycete aggregation, we combined transcriptomic and ultrastructural analyses. Results indicate involvement of electroception in K+ sensing. They establish that the transcriptome repertoire required for swimming and aggregation is already fully functional at zoospore release. At the time points analyzed, aggregates are mainly constituted of zoospores. They produce vesicular and fibrillary material discharged at cell-to-cell contacts. Consistently, the signature of transcriptome dynamics during transition to aggregates is an upregulation of genes potentially related to vesicular trafficking. Moreover, transcriptomic and functional analyses show a strong enhancement of carbonic anhydrase activity, indicating that pH homeostasis may contribute to aggregation by acting on both zoospore movement and adhesion. This study poses the molecular and cellular bases of aggregative behavior within oomycetes and expands the current knowledge of ion perception-mediated dissemination of propagules in the rhizosphere.


21. Novel dynamics of human mucociliary differentiation revealed by single-cell RNA sequencing of nasal epithelial cultures., Development. 2019 Sep 26. pii: dev.177428 (Pubmed: 31558434)
Ruiz García S, Deprez M, Lebrigand K, Cavard A, Paquet A, Arguel MJ, Magnone V, Truchi M, Caballero I, Leroy S, Marquette CH, Marcet B, Barbry P, Zaragosi LE

The upper airway epithelium, mainly composed of multiciliated, goblet, club and basal cells, ensures proper mucociliary function and can regenerate upon aggressions. In chronic airway diseases, defective repair leads to tissue remodeling. Delineating key drivers of differentiation dynamics can help understand how normal or pathological regeneration occurs.Using single-cell transcriptomics and lineage inference, we have unraveled trajectories from basal to luminal cells, providing novel markers for specific populations. We report that: (1) a precursor subgroup of multiciliated cells that we have entitled deuterosomal cells, is defined by specific markers, such as DEUP1, FOXN4, YPEL1, HES6 and CDC20B; (2) goblet cells can be precursors of multiciliated cells, thus explaining the presence of hybrid cells that co-express markers of goblet and multiciliated cells; (3) a repertoire of molecules involved in the regeneration process, such as keratins or components of the Notch, Wnt or BMP/TGFβ pathways can be established. Confirmations of our results on fresh human and pig airway samples, and on mouse tracheal cells, extend and confirm our conclusions regarding the molecular and cellular choreography at work during mucociliary epithelial differentiation.


22. The nuclear hypoxia-regulated NLUCAT1 long non-coding RNA contributes to an aggressive phenotype in lung adenocarcinoma through regulation of oxidative stress., Oncogene. 2019 Aug 15. doi: 10.1038/s41388-019-0935-y (Pubmed: 31417181)
Moreno Leon L, Gautier M, Allan R, Ilié M, Nottet N, Pons N, Paquet A, Lebrigand K, Truchi M, Fassy J, Magnone V, Kinnebrew G, Radovich M, Cheok MH, Barbry P, Vassaux G, Marquette CH, Ponzio G, Ivan M, Pottier N, Hofman P, Mari B, Rezzonico R

Lung cancer is the leading cause of cancer death worldwide, with poor prognosis and a high rate of recurrence despite early surgical removal. Hypoxic regions within tumors represent sources of aggressiveness and resistance to therapy. Although long non-coding RNAs (lncRNAs) are increasingly recognized as major gene expression regulators, their regulation and function following hypoxic stress are still largely unexplored. Combining profiling studies on early-stage lung adenocarcinoma (LUAD) biopsies and on A549 LUAD cell lines cultured in normoxic or hypoxic conditions, we identified a subset of lncRNAs that are both correlated with the hypoxic status of tumors and regulated by hypoxia in vitro. We focused on a new transcript, NLUCAT1, which is strongly upregulated by hypoxia in vitro and correlated with hypoxic markers and poor prognosis in LUADs. Full molecular characterization showed that NLUCAT1 is a large nuclear transcript composed of six exons and mainly regulated by NF-κB and NRF2 transcription factors. CRISPR-Cas9-mediated invalidation of NLUCAT1 revealed a decrease in proliferative and invasive properties, an increase in oxidative stress and a higher sensitivity to cisplatin-induced apoptosis. Transcriptome analysis of NLUCAT1-deficient cells showed repressed genes within the antioxidant and/or cisplatin-response networks. We demonstrated that the concomitant knockdown of four of these genes products, GPX2, GLRX, ALDH3A1, and PDK4, significantly increased ROS-dependent caspase activation, thus partially mimicking the consequences of NLUCAT1 inactivation in LUAD cells. Overall, we demonstrate that NLUCAT1 contributes to an aggressive phenotype in early-stage hypoxic tumors, suggesting it may represent a new potential therapeutic target in LUADs.


23. The Long Non-Coding RNA DNM3OS is a Reservoir of FibromiRs with Major Functions in Lung Fibroblast Response to TGF-β and Pulmonary Fibrosis, Am J Respir Crit Care Med. 2019 Apr 9. (Pubmed: 30964696)
Savary G, Dewaeles E, Diazzi S, Buscot M, Nottet N, Fassy J, Courcot E, Henaoui IS, Lemaire J, Martis N, Van der Hauwaert C, Pons N, Magnone V, Leroy S, Hofman V, Plantier L, Lebrigand K, Paquet A, Lino Cardenas CL, Vassaux G, Hofman P, Günther A, Crestani B, Wallaert B, Rezzonico R, Brousseau T, Glowacki F, Bellusci S, Perrais M, Broly F, Barbry P, Marquette CH, Cauffiez C, Mari B, Pottier N

RATIONALE: Given the paucity of effective treatments for Idiopathic Pulmonary Fibrosis (IPF), new insights into the deleterious mechanisms controlling lung fibroblast activation, the key cell type driving the fibrogenic process, are essential to develop new therapeutic strategies. Transforming growth factor β (TGF-β) is the main pro-fibrotic factor, but its inhibition is associated with severe side effects due to its pleiotropic role. OBJECTIVES: We hypothesized that downstream non-coding effectors of TGF-β in fibroblasts may represent new effective therapeutic targets whose modulation may be well-tolerated. METHODS: We investigated the whole non-coding fraction of TGF-β-stimulated lung fibroblast transcriptome to identify new genomic determinants of lung fibroblast differentiation into myofibroblast. Differential expression of the long non-coding RNA DNM3OS and its associated miRNAs was validated in a murine model of pulmonary fibrosis and in IPF tissue samples. Distinct and complementary antisense oligonucleotide-based strategies aiming at interfering with DNM3OS were used to elucidate the role of DNM3OS and its associated miRNAs in IPF pathogenesis. MEASUREMENTS AND MAIN RESULTS: We identified DNM3OS as a fibroblast-specific critical downstream effector of TGF-β-induced lung myofibroblast activation. Mechanistically, DNM3OS regulates this process in trans by giving rise to three distinct profibrotic mature miRNAs (i.e. miR-199a-5p/3p and miR-214-3p), which influence both SMAD and non-SMAD components of TGF-β signaling in a multifaceted way. In vivo, we showed that interfering with DNM3OS function not only prevents lung fibrosis but also improves established pulmonary fibrosis. CONCLUSION: Pharmacological approaches aiming at interfering with DNM3OS may represent new effective therapeutic strategies in IPF.


24. CD4+ T Cells Affect the Thyroid Hormone Transport at the Choroid Plexus in Mice Raised in Enriched Environment, Neuroimmunomodulation. 2019 Jan 31:1-8. doi: 10.1159/000495987 (Pubmed: 30703773)
Zarif H, Paquet A, Lebrigand K, Arguel MJ, Heurteaux C, Glaichenhaus N, Chabry J, Guyon A, Petit-Paitel A

Others and we have shown that T cells have an important role in hippocampal synaptic plasticity, including neurogenesis in the dentate gyrus, spinogenesis, and glutamatergic synaptic function in the CA of the hippocampus. Hippocampus plasticity is particularly involved in the brain effects of the enriched environment (EE), and interestingly CD4+ and CD8+ T cells play essential and differential roles in these effects. However, the precise mechanisms by which they act on the brain remain elusive. OBJECTIVES: We searched for a putative mechanism of action by which CD4+ T cells could influence brain plasticity and hypothesized that they could regulate protein transport at the level of the blood-CSF barrier in the choroid plexus. METHOD: We compared mice housed in EE and deprived of CD4+ T cells using a depleting antibody with a control group injected with the control isotype. We analyzed in the hippocampus the gene expression profiles using the Agilent system, and the expression of target proteins in plasma, CSF, and the choroid plexus using ELISA. RESULTS: We show that CD4+ T cells may influence EE-induced hippocampus plasticity via thyroid hormone signaling by regulating in the choroid plexus the expression of transthyretin, the major transporter of thyroxine (T4) to the brain parenchyma. CONCLUSIONS: Our study highlights the contribution of close interactions between the immune and neuroendocrine systems in brain plasticity and function.


25. Characterization of siRNAs clusters in Arabidopsis thaliana galls induced by the root-knot nematode Meloidogyne incognita, BMC Genomics. 2018 Dec 18;19 (Pubmed: 30563458)
Medina C, da Rocha M, Magliano M, Raptopoulo A, Marteu N, Lebrigand K, Abad P, Favery B, Jaubert-Possamai S

Root-knot nematodes (RKN), genus Meloidogyne, are plant parasitic worms that have the ability to transform root vascular cylinder cells into hypertrophied, multinucleate and metabolically over-active feeding cells. Redifferentiation into feeding cells is the result of a massive transcriptional reprogramming of root cells targeted by RKN. Since RKN are able to induce similar feeding cells in roots of thousands of plant species, these worms are thought to manipulate essential and conserved plant molecular pathways. Small non-coding RNAs of uninfected roots and infected root galls induced by M. incognita from Arabidopsis thaliana were sequenced by high throughput sequencing. SiRNA populations were analysed by using the Shortstack algorithm. We identified siRNA clusters that are differentially expressed in infected roots and evidenced an over-representation of the 23-24 nt siRNAs in infected tissue. This size corresponds to heterochromatic siRNAs (hc-siRNAs) which are known to regulate expression of transposons and genes at the transcriptional level, mainly by inducing DNA methylation. Correlation of siRNA clusters expression profile with transcriptomic data identified several protein coding genes that are candidates to be regulated by siRNAs at the transcriptional level by RNA directed DNA methylation (RdDM) pathway either directly or indirectly via silencing of neighbouring transposable elements.


26. The "one airway, one disease" concept in light of Th2 inflammation., Eur Respir J. 2018 Sep 6. pii: 1800437. doi: 10.1183/13993003.00437-2018. (Pubmed: 30190271)
Giovannini-Chami L, Paquet A, Sanfiorenzo C, Pons N, Cazareth J, Magnone V, Lebrigand K, Chevalier B, Vallauri A, Julia V, Marquette CH, Marcet B, Leroy S, Barbry P

In line with the pathophysiological continuum described between nose and bronchus in allergic respiratory diseases, we assessed whether nasal epithelium could mirror the Th2 status of bronchial epithelium.Nasal and bronchial cells were collected by brushings from patients with allergic rhinitis and asthma (AR, n=12), isolated allergic rhinitis (R, n=14) and healthy controls (C, n=13). Cellular composition was assessed by flow cytometry. Gene expression was analysed by RNA sequencing. Th2, Th17 and interferon signatures were derived from the literature.Infiltration by polymorphonuclear neutrophils in nose excluded 30% of the initial cohort. All bronchial samples from AR group were Th2-high. Nasal samples gene expression profile from the AR group correctly predicted the paired bronchial sample Th2 status in 71% of cases. Nevertheless, nasal cells did not appear as a reliable surrogate of the Th2 response, in particular due to a more robust influence of the interferon response in 14/26 nasal samples. Th2 scores correlated with mast cells counts (p<0.001) and numbers of sensitizations (p=0.006 and 0.002), while Th17 scores correlated with PMN counts (p<0.014).The large variability in nasal cell composition and type of inflammation restricts its use as a surrogate for assessing bronchial Th2 inflammation in AR patients.


27. CD4+ T Cells Have a Permissive Effect on Enriched Environment-Induced Hippocampus Synaptic Plasticity., Front Synaptic Neurosci. 2018 Jun 13;10:14. doi: 10.3389/fnsyn.2018.00014. eCollection 2018. (Pubmed: 29950983)
Zarif H, Hosseiny S, Paquet A, Lebrigand K, Arguel MJ, Cazareth J, Lazzari A, Heurteaux C, Glaichenhaus N, Chabry J, Guyon A, Petit-Paitel A

Living in an enriched environment (EE) benefits health by acting synergistically on various biological systems including the immune and the central nervous systems. The dialog between the brain and the immune cells has recently gained interest and is thought to play a pivotal role in beneficial effects of EE. Recent studies show that T lymphocytes have an important role in hippocampal plasticity, learning, and memory, although the precise mechanisms by which they act on the brain remain elusive. Using a mouse model of EE, we show here that CD4+ T cells are essential for spinogenesis and glutamatergic synaptic function in the CA of the hippocampus. However, CD4+ lymphocytes do not influence EE-induced neurogenesis in the DG of the hippocampus, by contrast to what we previously demonstrated for CD8+ T cells. Importantly, CD4+ T cells located in the choroid plexus have a specific transcriptomic signature as a function of the living environment. Our study highlights the contribution of CD4+ T cells in the brain plasticity and function.


28. Comparative Transcriptome Profiling of Virulent and Attenuated Ehrlichia ruminantium Strains Highlighted Strong Regulation of map1- and Metabolism Related Genes., Front Cell Infect Microbiol. 2018 May 15;8:153. doi: 10.3389/fcimb.2018.00153. eCollection 2018. (Pubmed: 29868509)
Pruneau L, Lebrigand K, Mari B, Lefrançois T, Meyer DF, Vachiery N

The obligate intracellular pathogenic bacterium, Ehrlichia ruminantium, is the causal agent of heartwater, a fatal disease in ruminants transmitted by Amblyomma ticks. So far, three strains have been attenuated by successive passages in mammalian cells. The attenuated strains have improved capacity for growth in vitro, whereas they induced limited clinical signs in vivo and conferred strong protection against homologous challenge. However, the mechanisms of pathogenesis and attenuation remain unknown. In order to improve knowledge of E. ruminantium pathogenesis, we performed a comparative transcriptomic analysis of two distant strains of E. ruminantium, Gardel and Senegal, and their corresponding attenuated strains. Overall, our results showed an upregulation of gene expression encoding for the metabolism pathway in the attenuated strains compared to the virulent strains, which can probably be associated with higher in vitro replicative activity and a better fitness to the host cells. We also observed a significant differential expression of membrane protein-encoding genes between the virulent and attenuated strains. A major downregulation of map1-related genes was observed for the two attenuated strains, whereas upregulation of genes encoding for hypothetical membrane proteins was observed for the four strains. Moreover, CDS_05140, which encodes for a putative porin, displays the highest gene expression in both attenuated strains. For the attenuated strains, the significant downregulation of map1-related gene expression and upregulation of genes encoding other membrane proteins could be important in the implementation of efficient immune responses after vaccination with attenuated vaccines. Moreover, this study revealed an upregulation of gene expression for 8 genes encoding components of Type IV secretion system and 3 potential effectors, mainly in the virulent Gardel strain. Our transcriptomic study, supported by previous proteomic studies, provides and also confirms new information regarding the characterization of genes involved in E. ruminantium virulence and attenuation mechanisms.


29. HITS-CLIP in various brain areas reveals new targets and new modalities of RNA binding by fragile X mental retardation protein, Nucleic Acids Res. 2018 Apr 14. doi: 10.1093/nar/gky267 (Pubmed: 29668986)
Maurin T, Lebrigand K, Castagnola S, Paquet A, Jarjat M, Popa A, Grossi M, Rage F, Bardoni B

Fragile X syndrome (FXS), the most common form of inherited intellectual disability, is due to the functional deficiency of the fragile X mental retardation protein (FMRP), an RNA-binding protein involved in translational regulation of many messenger RNAs, playing key roles in synaptic morphology and plasticity. To date, no effective treatment for FXS is available. We searched for FMRP targets by HITS-CLIP during early development of multiple mouse brain regions (hippocampus, cortex and cerebellum) at a time of brain development when FMRP is most highly expressed and synaptogenesis reaches a peak. We identified the largest dataset of mRNA targets of FMRP available in brain and we defined their cellular origin. We confirmed the G-quadruplex containing structure as an enriched motif in FMRP RNA targets. In addition to four less represented motifs, our study points out that, in the brain, CTGKA is the prominent motif bound by FMRP, which recognizes it when not engaged in Watson-Crick pairing. All of these motifs negatively modulated the expression level of a reporter protein. While the repertoire of FMRP RNA targets in cerebellum is quite divergent, the ones of cortex and hippocampus are vastly overlapping. In these two brain regions, the Phosphodiesterase 2a (Pde2a) mRNA is a prominent target of FMRP, which modulates its translation and intracellular transport. This enzyme regulates the homeostasis of cAMP and cGMP and represents a novel and attractive therapeutic target to treat FXS.


30. Genome evolution across 1,011 Saccharomyces cerevisiae isolates, Nature. 2018 Apr;556(7701):339-344. doi: 10.1038/s41586-018-0030-5. Epub 2018 Apr 11. (Pubmed: 29643504)
Peter J, De Chiara M, Friedrich A, Yue JX, Pflieger D, Bergström A, Sigwalt A, Barre B, Freel K, Llored A, Cruaud C, Labadie K, Aury JM, Istace B, Lebrigand K, Barbry P, Engelen S, Lemainque A, Wincker P, Liti G, Schacherer J

Large-scale population genomic surveys are essential to explore the phenotypic diversity of natural populations. Here we report the whole-genome sequencing and phenotyping of 1,011 Saccharomyces cerevisiae isolates, which together provide an accurate evolutionary picture of the genomic variants that shape the species-wide phenotypic landscape of this yeast. Genomic analyses support a single 'out-of-China' origin for this species, followed by several independent domestication events. Although domesticated isolates exhibit high variation in ploidy, aneuploidy and genome content, genome evolution in wild isolates is mainly driven by the accumulation of single nucleotide polymorphisms. A common feature is the extensive loss of heterozygosity, which represents an essential source of inter-individual variation in this mainly asexual species. Most of the single nucleotide polymorphisms, including experimentally identified functional polymorphisms, are present at very low frequencies. The largest numbers of variants identified by genome-wide association are copy-number changes, which have a greater phenotypic effect than do single nucleotide polymorphisms. This resource will guide future population genomics and genotype-phenotype studies in this classic model system.


31. CD8+ T cells are essential for the effects of enriched environment on hippocampus-dependent behavior, hippocampal neurogenesis and synaptic plasticity., Brain Behav Immun. 2017 Nov 22. pii: S0889-1591(17)30517-2. doi: 10.1016/j.bbi.2017.11.016. (Pubmed: 29175168)
Zarif H, Nicolas S, Guyot M, Hosseiny S, Lazzari A, Canali MM, Cazareth J, Brau F, Golzne V, Dourneau E, Maillaut M, Luci C, Paquet A, Lebrigand K, Arguel MJ, Daoudlarian D, Heurteaux C, Glaichenhaus N, Chabry J, Guyon A, Petit-Paitel A

Enriched environment (EE) induces plasticity changes in the brain. Recently, CD4+ T cells have been shown to be involved in brain plasticity processes. Here, we show that CD8+ T cells are required for EE-induced brain plasticity in mice, as revealed by measurements of hippocampal volume, neurogenesis in the DG of the hippocampus, spinogenesis and glutamatergic synaptic function in the CA of the hippocampus. As a consequence, EE-induced behavioral benefits depend, at least in part, on CD8+ T cells. In addition, we show that spleen CD8+ T cells from mice housed in standard environment (SE) and EE have different properties in terms of 1) TNFα release after in vitro CD3/CD28 or PMA/Iono stimulation 2) in vitro proliferation properties 3) CD8+ CD44+ CD62Llow and CD62Lhi T cells repartition 4) transcriptomic signature as revealed by RNA sequencing. CD8+ T cells purified from the choroid plexus of SE and EE mice also exhibit different transcriptomic profiles as highlighted by single-cell mRNA sequencing. We show that CD8+ T cells are essential mediators of beneficial EE effects on brain plasticity and cognition. Additionally, we propose that EE differentially primes CD8+ T cells leading to behavioral improvement.


32. Characterization of microRNAs from Arabidopsis galls highlights a role for miR159 in the plant response to the root-knot nematode Meloidogyne incognita, New Phytol. 2017 Sep 14. doi: 10.1111/nph.14717. (Pubmed: 28906559)
Medina C, da Rocha M, Magliano M, Ratpopoulo A, Revel B, Marteu N, Magnone V, Lebrigand K, Cabrera J, Barcala M, Silva AC, Millar A, Escobar C, Abad P, Favery B, Jaubert-Possamai S

Root knot nematodes (RKN) are root parasites that induce the genetic reprogramming of vascular cells into giant feeding cells and the development of root galls. MicroRNAs (miRNAs) regulate gene expression during development and plant responses to various stresses. Disruption of post-transcriptional gene silencing in Arabidopsis ago1 or ago2 mutants decrease the infection rate of RKN suggesting a role for this mechanism in the plant-nematode interaction. By sequencing small RNAs from uninfected Arabidopsis roots and from galls 7 and 14 d post infection with Meloidogyne incognita, we identified 24 miRNAs differentially expressed in gall as putative regulators of gall development. Moreover, strong activity within galls was detected for five miRNA promoters. Analyses of nematode development in an Arabidopsis miR159abc mutant had a lower susceptibility to RKN, suggesting a role for the miR159 family in the plant response to M. incognita. Localization of mature miR159 within the giant and surrounding cells suggested a role in giant cell and gall. Finally, overexpression of miR159 in galls at 14 d post inoculation was associated with the repression of the miR159 target MYB33 which expression is restricted to the early stages of infection. Overall, these results implicate the miR159 in plant responses to RKN.


33. Effects of proton versus photon irradiation on (lymph)angiogenic, inflammatory, proliferative and anti-tumor immune responses in head and neck squamous cell carcinoma., Oncogenesis. 2017 Jul 3;6(7):e354. doi: 10.1038/oncsis.2017.56 (Pubmed: 28671677)
Lupu-Plesu M, Claren A, Martial S, N'Diaye PD, Lebrigand K, Pons N, Ambrosetti D, Peyrottes I, Feuillade J, Hérault J, Dufies M, Doyen J, Pagès G

The proximity of organs at risk makes the treatment of head and neck squamous cell carcinoma (HNSCC) challenging by standard radiotherapy. The higher precision in tumor targeting of proton (P) therapy could promote it as the treatment of choice for HNSCC. Besides the physical advantage in dose deposition, few is known about the biological impact of P versus photons (X) in this setting. To investigate the comparative biological effects of P versus X radiation in HNSCC cells, we assessed the relative biological effectiveness (RBE), viability, proliferation and mRNA levels for genes involved in (lymph)angiogenesis, inflammation, proliferation and anti-tumor immunity. These parameters, particularly VEGF-C protein levels and regulations, were documented in freshly irradiated and/or long-term surviving cells receiving low/high-dose, single (SI)/multiple (MI) irradiations with P/X. The RBE was found to be 1.1 Key (lymph)angiogenesis and inflammation genes were downregulated (except for vegf-c) after P and upregulated after X irradiation in MI surviving cells, demonstrating a more favorable profile after P irradiation. Both irradiation types stimulated vegf-c promoter activity in a NF-κB-dependent transcriptional regulation manner, but at a lesser extent after P, as compared to X irradiation, which correlated with mRNA and protein levels. The cells surviving to MI by P or X generated tumors with higher volume, anarchic architecture and increased density of blood vessels. Increased lymphangiogenesis and a transcriptomic analysis in favor of a more aggressive phenotype were observed in tumors generated with X-irradiated cells. Increased detection of lymphatic vessels in relapsed tumors from patients receiving X radiotherapy was consistent with these findings. This study provides new data about the biological advantage of P, as compared to X irradiation. In addition to its physical advantage in dose deposition, P irradiation may help to improve treatment approaches for HNSCC.


34. Characterizing isomiR variants within the microRNA-34/449 family, FEBS Lett. 2017 Mar;591(5):693-705. doi: 10.1002/1873-3468.12595. Epub 2017 Feb 28 (Pubmed: 28192603)
Mercey O, Popa A, Cavard A, Paquet A, Chevalier B, Pons N, Magnone V, Zangari J, Brest P, Zaragosi LE, Ponzio G, Lebrigand K, Barbry P, Marcet B

miR-34/449 microRNAs are conserved regulators of multiciliated cell differentiation. Here, we evidence and characterize expression of two isomiR variant sequences from the miR-34/449 family in human airway epithelial cells. These isomiRs differ from their canonical counterparts miR-34b and miR-449c by one supplemental uridine at their 5'-end, leading to a one-base shift in their seed region. Overexpression of canonical miR-34/449 or 5'-isomiR-34/449 induces distinct gene expression profiles and biological effects. However, some target transcripts and functional activities are shared by both canonical microRNAs and isomiRs. Indeed, both repress important targets that result in cell cycle blockage and Notch pathway inhibition. Our findings suggest that 5'-isomiR-34/449 may represent additional mechanisms by which miR-34/449 family finely controls several pathways to drive multiciliogenesis.


35. A cost effective 5' selective single cell transcriptome profiling approach with improved UMI design, Nucleic Acids Res. 2016 Dec 9. pii: gkw1242. (Pubmed: 27940562)
Arguel MJ, Lebrigand K, Paquet A, Ruiz Garcia S, Zaragosi LE, Barbry P, Waldmann R

Single cell RNA sequencing approaches are instrumental in studies of cell-to-cell variability. 5' selective transcriptome profiling approaches allow simultaneous definition of the transcription start size and have advantages over 3' selective approaches which just provide internal sequences close to the 3' end. The only currently existing 5' selective approach requires costly and labor intensive fragmentation and cell barcoding after cDNA amplification. We developed an optimized 5' selective workflow where all the cell indexing is done prior to fragmentation. With our protocol, cell indexing can be performed in the Fluidigm C1 microfluidic device, resulting in a significant reduction of cost and labor. We also designed optimized unique molecular identifiers that show less sequence bias and vulnerability towards sequencing errors resulting in an improved accuracy of molecule counting. We provide comprehensive experimental workflows for Illumina and Ion Proton sequencers that allow single cell sequencing in a cost range comparable to qPCR assays.


36. miR-200 family controls late steps of postnatal forebrain neurogenesis via Zeb2 inhibition., Sci Rep. 2016 Oct 21;6:35729. doi: 10.1038/srep35729. (Pubmed: 27767083)
Beclin C, Follert P, Stappers E, Barral S, Nathalie C, de Chevigny A, Magnone V, Lebrigand K, Bissels U, Huylebroeck D, Bosio A, Barbry P, Seuntjens E, Cremer H

During neurogenesis, generation, migration and integration of the correct numbers of each neuron sub-type depends on complex molecular interactions in space and time. MicroRNAs represent a key control level allowing the flexibility and stability needed for this process. Insight into the role of this regulatory pathway in the brain is still limited. We performed a sequential experimental approach using postnatal olfactory bulb neurogenesis in mice, starting from global expression analyses to the investigation of functional interactions between defined microRNAs and their targets. Deep sequencing of small RNAs extracted from defined compartments of the postnatal neurogenic system demonstrated that the miR-200 family is specifically induced during late neuronal differentiation stages. Using in vivo strategies we interfered with the entire miR-200 family in loss- and gain-of-function settings, showing a role of miR-200 in neuronal maturation. This function is mediated by targeting the transcription factor Zeb2. Interestingly, so far functional interaction between miR-200 and Zeb2 has been exclusively reported in cancer or cultured stem cells. Our data demonstrate that this regulatory interaction is also active during normal neurogenesis.


37. RiboProfiling: a Bioconductor package for standard Ribo-seq pipeline processing., F1000Res. 2016 Jun 9;5:1309. doi: 10.12688/f1000research.8964.1. eCollection 2016. (Pubmed: 27347386)
Popa A, Lebrigand K, Paquet A, Nottet N, Robbe-Sermesant K, Waldmann R, Barbry P

The ribosome profiling technique (Ribo-seq) allows the selective sequencing of translated RNA regions. Recently, the analysis of genomic sequences associated to Ribo-seq reads has been widely employed to assess their coding potential. These analyses led to the identification of differentially translated transcripts under different experimental conditions, and/or ribosome pausing on codon motifs. In the context of the ever-growing need for tools analyzing Ribo-seq reads, we have developed 'RiboProfiling', a new Bioconductor open-source package. 'RiboProfiling' provides a full pipeline to cover all key steps for the analysis of ribosome footprints. This pipeline has been implemented in a single R workflow. The package takes an alignment (BAM) file as input and performs ribosome footprint quantification at a transcript level. It also identifies footprint accumulation on particular amino acids or multi amino-acids motifs. Report summary graphs and data quantification are generated automatically. The package facilitates quality assessment and quantification of Ribo-seq experiments. Its implementation in Bioconductor enables the modeling and statistical analysis of its output through the vast choice of packages available in R. This article illustrates how to identify codon-motifs accumulating ribosome footprints, based on data from Escherichia coli.


38. Comparative Genomic Analysis of Drechmeria coniospora Reveals Core and Specific Genetic Requirements for Fungal Endoparasitism of Nematodes., PLoS Genet. 2016 May 6;12(5):e1006017. doi: 10.1371/journal.pgen.1006017. eCollection 2016. (Pubmed: 27153332)
Lebrigand K, He le D, Thakur N, Arguel MJ, Polanowska J, Henrissat B, Record E, Magdelenat G, Barbe V, Raffaele S, Barbry P, Ewbank JJ

Drechmeria coniospora is an obligate fungal pathogen that infects nematodes via the adhesion of specialized spores to the host cuticle. D. coniospora is frequently found associated with Caenorhabditis elegans in environmental samples. It is used in the study of the nematode's response to fungal infection. Full understanding of this bi-partite interaction requires knowledge of the pathogen's genome, analysis of its gene expression program and a capacity for genetic engineering. The acquisition of all three is reported here. A phylogenetic analysis placed D. coniospora close to the truffle parasite Tolypocladium ophioglossoides, and Hirsutella minnesotensis, another nematophagous fungus. Ascomycete nematopathogenicity is polyphyletic; D. coniospora represents a branch that has not been molecularly characterized. A detailed in silico functional analysis, comparing D. coniospora to 11 fungal species, revealed genes and gene families potentially involved in virulence and showed it to be a highly specialized pathogen. A targeted comparison with nematophagous fungi highlighted D. coniospora-specific genes and a core set of genes associated with nematode parasitism. A comparative gene expression analysis of samples from fungal spores and mycelia, and infected C. elegans, gave a molecular view of the different stages of the D. coniospora lifecycle. Transformation of D. coniospora allowed targeted gene knock-out and the production of fungus that expresses fluorescent reporter genes. It also permitted the initial characterisation of a potential fungal counter-defensive strategy, involving interference with a host antimicrobial mechanism. This high-quality annotated genome for D. coniospora gives insights into the evolution and virulence of nematode-destroying fungi. Coupled with genetic transformation, it opens the way for molecular dissection of D. coniospora physiology, and will allow both sides of the interaction between D. coniospora and C. elegans, as well as the evolutionary arms race that exists between pathogen and host, to be studied.


39. Sequential activation and distinct functions for distal and proximal modules within the IgH 3' regulatory region., Proc Natl Acad Sci U S A. 2016 Feb 1. pii: 201514090 (Pubmed: 26831080)
Garot A, Marquet M, Saintamand A, Bender S, Le Noir S, Rouaud P, Carrion C1, Oruc Z, Bébin AG, Moreau J, Lebrigand K, Denizot Y, Alt FW, Cogné M, Pinaud E

As a master regulator of functional Ig heavy chain (IgH) expression, the IgH 3' regulatory region (3'RR) controls multiple transcription events at various stages of B-cell ontogeny, from newly formed B cells until the ultimate plasma cell stage. The IgH 3'RR plays a pivotal role in early B-cell receptor expression, germ-line transcription preceding class switch recombination, interactions between targeted switch (S) regions, variable region transcription before somatic hypermutation, and antibody heavy chain production, but the functional ranking of its different elements is still inaccurate, especially that of its evolutionarily conserved quasi-palindromic structure. By comparing relevant previous knockout (KO) mouse models (3'RR KO and hs3b-4 KO) to a novel mutant devoid of the 3'RR quasi-palindromic region (3'PAL KO), we pinpointed common features and differences that specify two distinct regulatory entities acting sequentially during B-cell ontogeny. Independently of exogenous antigens, the 3'RR distal part, including hs4, fine-tuned B-cell receptor expression in newly formed and naïve B-cell subsets. At mature stages, the 3'RR portion including the quasi-palindrome dictated antigen-dependent locus remodeling (global somatic hypermutation and class switch recombination to major isotypes) in activated B cells and antibody production in plasma cells.


40. Pateamine A-sensitive ribosome profiling reveals the scope of translation in mouse embryonic stem cells., BMC Genomics. 2016 Jan 14;17(1):52. doi: 10.1186/s12864-016-2384-0. (Pubmed: 26764022)
Popa A, Lebrigand K, Barbry P, Waldmann R

BACKGROUND: Open reading frames are common in long noncoding RNAs (lncRNAs) and 5'UTRs of protein coding transcripts (uORFs). The question of whether those ORFs are translated was recently addressed by several groups using ribosome profiling. Most of those studies concluded that certain lncRNAs and uORFs are translated, essentially based on computational analysis of ribosome footprints. However, major discrepancies remain on the scope of translation and the translational status of individual ORFs. In consequence, further criteria are required to reliably identify translated ORFs from ribosome profiling data. RESULTS: We examined the effect of the translation inhibitors pateamine A, harringtonine and puromycin on murine ES cell ribosome footprints. We found that pateamine A, a drug that targets eIF4A, allows a far more accurate identification of translated sequences than previously used drugs and computational scoring schemes. Our data show that at least one third but less than two thirds of ES cell lncRNAs are translated. We also identified translated uORFs in hundreds of annotated coding transcripts including key pluripotency transcripts, such as dicer, lin28, trim71, and ctcf. CONCLUSION: Pateamine A inhibition data clearly increase the precision of the detection of translated ORFs in ribosome profiling experiments. Our data show that translation of lncRNAs and uORFs in murine ES cells is rather common although less pervasive than previously suggested. The observation of translated uORFs in several key pluripotency transcripts suggests that translational regulation by uORFs might be part of the network that defines mammalian stem cell identity.


41. MicroRNA target identification: lessons from hypoxamiRs., Antioxid Redox Signal. 2014 Sep 10;21(8):1249-68. doi: 10.1089/ars.2013.5648. Epub 2014 Feb 3. (Pubmed: 24111877)
Bertero T, Robbe-Sermesant K, Le Brigand K, Ponzio G, Pottier N, Rezzonico R, Mazure NM, Barbry P, Mari B

SIGNIFICANCE: MicroRNAs (miRNAs) are small noncoding RNAs that have emerged as key regulators of many physiological and pathological processes, including those relevant to hypoxia such as cancer, neurological dysfunctions, myocardial infarction, and lung diseases. RECENT ADVANCES: During the last 5 years, miRNAs have been shown to play a role in the regulation of the cellular response to hypoxia. The identification of several bona fide targets of these hypoxamiRs has underlined their pleiotropic functions and the complexity of the molecular rules directing miRNA::target transcript pairing. CRITICAL ISSUES: This review outlines the main in silico and experimental approaches used to identify the targetome of hypoxamiRs and presents new recent relevant methodologies for future studies. FUTURE DIRECTIONS: Since hypoxia plays key roles in many pathophysiological conditions, the precise characterization of regulatory hypoxamiRs networks will be instrumental both at a fundamental level and for their future potential therapeutic applications.


42. miR-193b/365a cluster controls progression of epidermal squamous cell carcinoma., Carcinogenesis. 2014 May;35(5):1110-20. doi: 10.1093/carcin/bgt490. Epub 2013 Dec 28. (Pubmed: 24374827)
Gastaldi C, Bertero T, Xu N, Bourget-Ponzio I, Lebrigand K, Fourre S, Popa A, Cardot-Leccia N, Meneguzzi G, Sonkoly E, Pivarcsi A, Mari B, Barbry P, Ponzio G, Rezzonico R

Incidence of cutaneous squamous cell carcinomas (cSCCs) constantly increases in the Caucasian population. Developing preferentially on precancerous lesions such as actinic keratoses due to chronic sunlight exposure, cSCCs result from the malignant transformation of keratinocytes. Although a resection of the primary tumor is usually curative, a subset of aggressive cSCCs shows a high risk of recurrence and metastases. The characterization of the molecular dysfunctions involved in cSCC development should help to identify new relevant targets against these aggressive cSCCs. In that context, we have used small RNA sequencing to identify 100 microRNAs (miRNAs) whose expression was altered during chemically induced mouse skin tumorigenesis. The decreased expression of the miR-193b/365a cluster during tumor progression suggests a tumor suppressor role. Ectopic expression of these miRNAs in tumor cells indeed inhibited their proliferation, clonogenic potential and migration, which were stimulated in normal keratinocytes when these miRNAs were blocked with antisense oligonucleotides. A combination of in silico predictions and transcriptome analyses identified several target genes of interest. We validated KRAS and MAX as direct targets of miR-193b and miR-365a. Repression of these targets using siRNAs mimicked the effects of miR-193b and miR-365a, suggesting that these genes might mediate, at least in part, the tumor-suppressive action of these miRNAs.


43. "Seed-Milarity" confers to hsa-miR-210 and hsa-miR-147b similar functional activity., PLoS One. 2012;7(9):e44919. doi: 10.1371/journal.pone.0044919. Epub 2012 Sep 13. (Pubmed: 23028679)
Bertero T, Grosso S, Robbe-Sermesant K, Lebrigand K, Henaoui IS, Puissegur MP, Fourre S, Zaragosi LE, Mazure NM, Ponzio G, Cardinaud B, Barbry P, Rezzonico R, Mari B

Specificity of interaction between a microRNA (miRNA) and its targets crucially depends on the seed region located in its 5'-end. It is often implicitly considered that two miRNAs sharing the same biological activity should display similarity beyond the strict six nucleotide region that forms the seed, in order to form specific complexes with the same mRNA targets. We have found that expression of hsa-miR-147b and hsa-miR-210, though triggered by different stimuli (i.e. lipopolysaccharides and hypoxia, respectively), induce very similar cellular effects in term of proliferation, migration and apoptosis. Hsa-miR-147b only shares a "minimal" 6-nucleotides seed sequence with hsa-miR-210, but is identical with hsa-miR-147a over 20 nucleotides, except for one base located in the seed region. Phenotypic changes induced after heterologous expression of miR-147a strikingly differ from those induced by miR-147b or miR-210. In particular, miR-147a behaves as a potent inhibitor of cell proliferation and migration. These data fit well with the gene expression profiles observed for miR-147b and miR-210, which are very similar, and the gene expression profile of miR-147a, which is distinct from the two others. Bioinformatics analysis of all human miRNA sequences indicates multiple cases of miRNAs from distinct families exhibiting the same kind of similarity that would need to be further characterized in terms of putative functional redundancy. Besides, it implies that functional impact of some miRNAs can be masked by robust expression of miRNAs belonging to distinct families.


44. Distinct epithelial gene expression phenotypes in childhood respiratory allergy., Eur Respir J. 2012 May;39(5):1197-205. Epub 2011 Oct 17. (Pubmed: 22005912)
Giovannini-Chami L, Marcet B, Moreilhon C, Chevalier B, Illie MI, Lebrigand K, Robbe-Sermesant K, Bourrier T, Michiels JF, Mari B, Crénesse D, Hofman P, de Blic J, Castillo L, Albertini M, Barbry P

Epithelial cell contribution to the natural history of childhood allergic respiratory disease remains poorly understood. Our aims were to identify epithelial pathways that are dysregulated in different phenotypes of respiratory allergy. We established gene expression signatures of nasal brushings from children with dust mite-allergic rhinitis, associated or not associated with controlled or uncontrolled asthma. Supervised learning and unsupervised clustering were used to predict the different subgroups of patients and define altered signalling pathways. These profiles were compared with those of primary cultures of human nasal epithelial cells stimulated with either interleukin (IL)-4, IL-13, interferon (IFN)-α, IFN-β or IFN-γ, or during in vitro differentiation. A supervised method discriminated children with allergic rhinitis from healthy controls (prediction accuracy 91%), based on 61 transcripts, including 21 T-helper cell (Th) type 2-responsive genes. This method was then applied to predict children with controlled or uncontrolled asthma (prediction accuracy 75%), based on 41 transcripts: nine of them, which were down-regulated in uncontrolled asthma, are directly linked to IFN. This group also included GSDML, which is genetically associated with asthma. Our data revealed a Th2-driven epithelial phenotype common to all children with dust mite allergic rhinitis. It highlights the influence of epithelially expressed molecules on the control of asthma, in association with atopy and impaired viral response.


45. B-cell regulator of immunoglobulin heavy chain transcription (Bright)/ARID3a is a direct target of the oncomir microRNA-125b in progenitor B-cells., Leukemia. 2012 Apr 3. doi: 10.1038/leu.2012.95. (Pubmed: 22469780)
Puissegur MP, Eichner R, Quelen C, Coyaud E, Mari B, Lebrigand K, Broccardo C, Nguyen-Khac F, Bousquet M, Brousset P

B-cell acute lymphoblastic leukemia (B-ALL) is often associated with chromosomal translocations leading to the deregulation of proto-oncogenes. MicroRNAs can also be affected by chromosomal alterations and thus contribute to carcinogenesis. The microRNA miR-125b-1 is over-expressed in B-ALL cases with the t(11;14)(q24;q34) translocation, therefore we sought to determine the role of this microRNA in B-cell fate. We used murine pre-BI cells alongside murine and human leukemic B-cell lines to show that miR-125b expression enhances proliferation by targeting Bright/ARID3a, an activator of immunoglobulin heavy-chain transcription. Accordingly, this target gene was down-regulated in B-ALL patients with the t(11;14)(q24;q34) translocation. Repression of Bright/ARID3A blocked differentiation and conferred a survival advantage to Ba/F3 cells under IL3 starvation. In addition, over-expression of miR-125b protected pre-BI and leukemic B-cell lines from apoptosis through blockade of caspase activation via a mechanism that was independent of p53 and BAK1. In summary, miR-125b can act as an oncogene in B-ALL by targeting ARID3a and mediating its repression, thus leading to a blockage in differentiation, increased proliferation and inhibition of apoptosis.Leukemia accepted article preview online, 3 April 2012; doi:10.1038/leu.2012.95.


46. Can the microRNA signature distinguish between thyroid tumors of uncertain malignant potential and other well-differentiated tumors of the thyroid gland?, Endocr Relat Cancer. 2011 Sep 13;18(5):579-94. Print 2011 Oct. (Pubmed: 21778212)
Lassalle S, Hofman V, Ilie M, Bonnetaud C, Puisségur MP, Brest P, Loubatier C, Guevara N, Bordone O, Cardinaud B, Lebrigand K, Rios G, Santini J, Franc B, Mari B, Al Ghuzlan A, Vielh P, Barbry P, Hofman P

The term 'thyroid tumors of uncertain malignant potential' (TT-UMP) was coined by surgical pathologists to define well-differentiated tumors (WDT) showing inconclusive morphological evidence of malignancy or benignity. We have analyzed the expression of microRNA (miRNA) in a training set of 42 WDT of different histological subtypes: seven follicular tumors of UMP (FT-UMP), six WDT-UMP, seven follicular thyroid adenomas (FTA), 11 conventional papillary thyroid carcinomas (C-PTC), five follicular variants of PTC (FV-PTC), and six follicular thyroid carcinomas (FTC), which led to the identification of about 40 deregulated miRNAs. A subset of these altered miRNAs was independently validated by qRT-PCR, which included 18 supplementary TT-UMP (eight WDT-UMP and ten FT-UMP). Supervised clustering techniques were used to predict the first 42 samples. Based on the four possible outcomes (FTA, C-PTC, FV-PTC, and FTC), about 80% of FTA and C-PTC and 50% of FV-PTC and FTC samples were correctly assigned. Analysis of the independent set of 18 WDT-UMP by quantitative RT-PCR for the selection of the six most discriminating miRNAs was unable to separate FT-UMP from WDT-UMP, suggesting that the miRNA signature is insufficient in characterizing these two clinical entities. We conclude that considering FT-UMP and WDT-UMP as distinct and specific clinical entities may improve the diagnosis of WDT of the thyroid gland. In this context, a small set of miRNAs (i.e. miR-7, miR-146a, miR-146b, miR-200b, miR-221, and miR-222) appears to be useful, though not sufficient per se, in distinguishing TT-UMP from other WDT of the thyroid gland.


47. Small RNA sequencing reveals miR-642a-3p as a novel adipocyte-specific microRNA and miR-30 as a key regulator of human adipogenesis., Genome Biol. 2011 Jul 18;12(7):R64. (Pubmed: 21767385)
Zaragosi LE, Wdziekonski B, Brigand KL, Villageois P, Mari B, Waldmann R, Dani C, Barbry P

BACKGROUND: In severe obesity, as well as in normal development, the growth of adipose tissue is the result of an increase in adipocyte size and numbers, which is underlain by the stimulation of adipogenic differentiation of precursor cells. A better knowledge of the pathways that regulate adipogenesis is therefore essential for an improved understanding of adipose tissue expansion. As microRNAs (miRNAs) have a critical role in many differentiation processes, our study aimed to identify the role of miRNA-mediated gene silencing in the regulation of adipogenic differentiation. RESULTS: We used deep sequencing to identify small RNAs that are differentially expressed during adipogenesis of adipose tissue-derived stem cells. This approach revealed the un-annotated miR-642a-3p as a highly adipocyte-specific miRNA. We then focused our study on the miR-30 family, which was also up-regulated during adipogenic differentiation and for which the role in adipogenesis had not yet been elucidated. Inhibition of the miR-30 family blocked adipogenesis, whilst over-expression of miR-30a and miR-30d stimulated this process. We additionally showed that both miR-30a and miR-30d target the transcription factor RUNX2, and stimulate adipogenesis via the modulation of this major regulator of osteogenesis. CONCLUSIONS: Overall, our data suggest that the miR-30 family plays a central role in adipocyte development. Moreover, as adipose tissue-derived stem cells can differentiate into either adipocytes or osteoblasts, the down-regulation of the osteogenesis regulator RUNX2 represents a plausible mechanism by which miR-30 miRNAs may contribute to adipogenic differentiation of adipose tissue-derived stem cells.


48. The human TTAGGG repeat factors 1 and 2 bind to a subset of interstitial telomeric sequences and satellite repeats., Cell Res. 2011 Jul;21(7):1028-38. doi: 10.1038/cr.2011.40. Epub 2011 Mar 22. (Pubmed: 21423270)
Simonet T, Zaragosi LE, Philippe C, Lebrigand K, Schouteden C, Augereau A, Bauwens S, Ye J, Santagostino M, Giulotto E, Magdinier F, Horard B, Barbry P, Waldmann R, Gilson E

The study of the proteins that bind to telomeric DNA in mammals has provided a deep understanding of the mechanisms involved in chromosome-end protection. However, very little is known on the binding of these proteins to nontelomeric DNA sequences. The TTAGGG DNA repeat proteins 1 and 2 (TRF1 and TRF2) bind to mammalian telomeres as part of the shelterin complex and are essential for maintaining chromosome end stability. In this study, we combined chromatin immunoprecipitation with high-throughput sequencing to map at high sensitivity and resolution the human chromosomal sites to which TRF1 and TRF2 bind. While most of the identified sequences correspond to telomeric regions, we showed that these two proteins also bind to extratelomeric sites. The vast majority of these extratelomeric sites contains interstitial telomeric sequences (or ITSs). However, we also identified non-ITS sites, which correspond to centromeric and pericentromeric satellite DNA. Interestingly, the TRF-binding sites are often located in the proximity of genes or within introns. We propose that TRF1 and TRF2 couple the functional state of telomeres to the long-range organization of chromosomes and gene regulation networks by binding to extratelomeric sequences.


49. A synonymous variant in IRGM alters a binding site for miR-196 and causes deregulation of IRGM-dependent xenophagy in Crohn's disease., Nat Genet. 2011 Mar;43(3):242-5. Epub 2011 Jan 30. (Pubmed: 21278745)
Brest P, Lapaquette P, Souidi M, Lebrigand K, Cesaro A, Vouret-Craviari V, Mari B, Barbry P, Mosnier JF, Hébuterne X, Harel-Bellan A, Mograbi B, Darfeuille-Michaud A, Hofman P

Susceptibility to Crohn's disease, a complex inflammatory disease, is influenced by common variants at many loci. The common exonic synonymous SNP (c.313C>T) in IRGM, found in strong linkage disequilibrium with a deletion polymorphism, has been classified as non-causative because of the absence of an alteration in the IRGM protein sequence or splice sites. Here we show that a family of microRNAs (miRNAs), miR-196, is overexpressed in the inflammatory intestinal epithelia of individuals with Crohn's disease and downregulates the IRGM protective variant (c.313C) but not the risk-associated allele (c.313T). Subsequent loss of tight regulation of IRGM expression compromises control of intracellular replication of Crohn's disease-associated adherent invasive Escherichia coli by autophagy. These results suggest that the association of IRGM with Crohn's disease arises from a miRNA-based alteration in IRGM regulation that affects the efficacy of autophagy, thereby implicating a synonymous polymorphism as a likely causal variant.


50. Coxiella burnetii transcriptional analysis reveals serendipity clusters of regulation in intracellular bacteria., PLoS One. 2010 Dec 21;5(12):e15321. (Pubmed: 21203564)
Leroy Q, Lebrigand K, Armougom F, Barbry P, Thiéry R, Raoult D

Coxiella burnetii, the causative agent of the zoonotic disease Q fever, is mainly transmitted to humans through an aerosol route. A spore-like form allows C. burnetii to resist different environmental conditions. Because of this, analysis of the survival strategies used by this bacterium to adapt to new environmental conditions is critical for our understanding of C. burnetii pathogenicity. Here, we report the early transcriptional response of C. burnetii under temperature stresses. Our data show that C. burnetii exhibited minor changes in gene regulation under short exposure to heat or cold shock. While small differences were observed, C. burnetii seemed to respond similarly to cold and heat shock. The expression profiles obtained using microarrays produced in-house were confirmed by quantitative RT-PCR. Under temperature stresses, 190 genes were differentially expressed in at least one condition, with a fold change of up to 4. Globally, the differentially expressed genes in C. burnetii were associated with bacterial division, (p)ppGpp synthesis, wall and membrane biogenesis and, especially, lipopolysaccharide and peptidoglycan synthesis. These findings could be associated with growth arrest and witnessed transformation of the bacteria to a spore-like form. Unexpectedly, clusters of neighboring genes were differentially expressed. These clusters do not belong to operons or genetic networks; they have no evident associated functions and are not under the control of the same promoters. We also found undescribed but comparable clusters of regulation in previously reported transcriptomic analyses of intracellular bacteria, including Rickettsia sp. and Listeria monocytogenes. The transcriptomic patterns of C. burnetii observed under temperature stresses permits the recognition of unpredicted clusters of regulation for which the trigger mechanism remains unidentified but which may be the result of a new mechanism of epigenetic regulation.


51. MiRonTop: mining microRNAs targets across large scale gene expression studies, Bioinformatics. 2010 Oct 19 (Pubmed: 20959382)
Le Brigand K, Robbe-Sermesant K, Mari B, Barbry P

SUMMARY: Current challenges in microRNA (miRNA) research are to improve the identification of in vivo mRNA targets and clarify the complex interplay existing between a specific miRNA and multiple biological networks. MiRonTop is an online java web tool that integrates DNA microarrays or high-throughput sequencing data to identify the potential implication of miRNAs on a specific biological system. It allows a rapid characterization of the most pertinent mRNA targets according to several existing miRNA target prediction approaches. It also provides useful representations of the enrichment scores according to the position of the target site along the 3'- UTR, where the contribution of the sites located in the vicinity of the stop codon and of the polyA tail can be clearly highlighted. It provides different graphs of miRNA enrichment associated with up- or down-regulated transcripts and different summary tables about selections of mRNA targets and their functional annotations by Gene Ontology. AVAILABILITY: http://www.microarray.fr:8080/miRonTop/index


52. miR-210 is overexpressed in late stages of lung cancer and mediates mitochondrial alterations associated with modulation of HIF-1 activity., Cell Death Differ. 2010 Oct 1. (Pubmed: 20885442)
Puissegur MP, Mazure NM, Bertero T, Pradelli L, Grosso S, Robbe-Sermesant K, Maurin T, Lebrigand K, Cardinaud B, Hofman V, Fourre S, Magnone V, Ricci JE, Pouyssegur J, Gounon P, Hofman P, Barbry P, Mari B

Following the identification of a set of hypoxia-regulated microRNAs (miRNAs), recent studies have highlighted the importance of miR-210 and of its transcriptional regulation by the transcription factor hypoxia-inducible factor-1 (HIF-1). We report here that miR-210 is overexpressed at late stages of non-small cell lung cancer. Expression of miR-210 in lung adenocarcinoma A549 cells caused an alteration of cell viability associated with induction of caspase-3/7 activity. miR-210 induced a loss of mitochondrial membrane potential and the apparition of an aberrant mitochondrial phenotype. The expression profiling of cells overexpressing miR-210 revealed a specific signature characterized by enrichment for transcripts related to 'cell death' and 'mitochondrial dysfunction', including several subunits of the electron transport chain (ETC) complexes I and II. The transcript coding for one of these ETC components, SDHD, subunit D of succinate dehydrogenase complex (SDH), was validated as a bona fide miR-210 target. Moreover, SDHD knockdown mimicked miR-210-mediated mitochondrial alterations. Finally, miR-210-dependent targeting of SDHD was able to activate HIF-1, in line with previous studies linking loss-of-function SDH mutations to HIF-1 activation. miR-210 can thus regulate mitochondrial function by targeting key ETC component genes with important consequences on cell metabolism, survival and modulation of HIF-1 activity. These observations help explain contradictory data regarding miR-210 expression and its putative function in solid tumors.


53. miR-34b/miR-34c: a regulator of TCL1 expression in 11q- chronic lymphocytic leukaemia?, Leukemia. 2009 Nov;23(11):2174-7. (Pubmed: 19536169)
Cardinaud B, Moreilhon C, Marcet B, Robbe-Sermesant K, Lebrigand K, Mari B, Eclache V, Cymbalista F, Raynaud S, Barbry P


54. Identification of keratinocyte growth factor as a target of microRNA-155 in lung fibroblasts: implication in epithelial-mesenchymal interactions., PLoS One. 2009 Aug 24;4(8):e6718. (Pubmed: 19701459)
Pottier N, Maurin T, Chevalier B, Puissegur MP, Lebrigand K, Robbe-Sermesant K, Bertero T, Lino Cardenas CL, Courcot E, Rios G, Fourre S, Lo-Guidice JM, Marcet B, Cardinaud B, Barbry P, Mari B

BACKGROUND: Epithelial-mesenchymal interactions are critical in regulating many aspects of vertebrate embryo development, and for the maintenance of homeostatic equilibrium in adult tissues. The interactions between epithelium and mesenchyme are believed to be mediated by paracrine signals such as cytokines and extracellular matrix components secreted from fibroblasts that affect adjacent epithelia. In this study, we sought to identify the repertoire of microRNAs (miRNAs) in normal lung human fibroblasts and their potential regulation by the cytokines TNF-alpha, IL-1beta and TGF-beta. METHODOLOGY/PRINCIPAL FINDINGS: MiR-155 was significantly induced by inflammatory cytokines TNF-alpha and IL-1beta while it was down-regulated by TGF-beta. Ectopic expression of miR-155 in human fibroblasts induced modulation of a large set of genes related to "cell to cell signalling", "cell morphology" and "cellular movement". This was consistent with an induction of caspase-3 activity and with an increase in cell migration in fibroblasts tranfected with miR-155. Using different miRNA bioinformatic target prediction tools, we found a specific enrichment for miR-155 predicted targets among the population of down-regulated transcripts. Among fibroblast-selective targets, one interesting hit was keratinocyte growth factor (KGF, FGF-7), a member of the fibroblast growth factor (FGF) family, which owns two potential binding sites for miR-155 in its 3'-UTR. Luciferase assays experimentally validated that miR-155 can efficiently target KGF 3'-UTR. Site-directed mutagenesis revealed that only one out of the 2 potential sites was truly functional. Functional in vitro assays experimentally validated that miR-155 can efficiently target KGF 3'-UTR. Furthermore, in vivo experiments using a mouse model of lung fibrosis showed that miR-155 expression level was correlated with the degree of lung fibrosis. CONCLUSIONS/SIGNIFICANCE: Our results strongly suggest a physiological function of miR-155 in lung fibroblasts. Altogether, this study implicates this miRNA in the regulation by mesenchymal cells of surrounding lung epithelium, making it a potential key player during tissue injury.


55. Gene expression profiling of imatinib and PD166326-resistant CML cell lines identifies Fyn as a gene associated with resistance to BCR-ABL inhibitors., Mol Cancer Ther. 2009 Jul;8(7):1924-33. Epub 2009 Jun 30. (Pubmed: 19567819)
Grosso S, Puissant A, Dufies M, Colosetti P, Jacquel A, Lebrigand K, Barbry P, Deckert M, Cassuto JP, Mari B, Auberger P

Imatinib is used to treat chronic myelogenous leukemia (CML), but resistance develops in all phases of this disease. The purpose of the present study was to identify the mode of resistance of newly derived imatinib-resistant (IM-R) and PD166326-resistant (PD-R) CML cells. IM-R and PD-R clones exhibited an increase in viability and a decrease in caspase activation in response to various doses of imatinib and PD166326, respectively, as compared with parental K562 cells. Resistance involved neither mutations in BCR-ABL nor increased BCR-ABL, MDR1 or Lyn expression, all known modes of resistance. To gain insight into the resistance mechanisms, we used pangenomic microarrays and identified 281 genes modulated in parental versus IM-R and PD-R cells. The gene signature was similar for IM-R and PD-R cells, accordingly with the cross-sensitivity observed for both inhibitors. These genes were functionally associated with pathways linked to development, cell adhesion, cell growth, and the JAK-STAT cascade. Especially relevant were the increased expression of the tyrosine kinases AXL and Fyn as well as CD44 and HMGA2. Small interfering RNA experiments and pharmacologic approaches identified FYN as a candidate for resistance to imatinib. Our findings provide a comprehensive picture of the transcriptional events associated with imatinib and PD166326 resistance and identify Fyn as a new potential target for therapeutic intervention in CML.


56. Global analysis of DNA methylation and transcription of human repetitive sequences., Epigenetics. 2009 Jul 17;4(5). (Pubmed: 19633427)
Horard B, Eymery A, Fourel G, Vassetzky N, Puechberty J, Roizes G, Lebrigand K, Barbry P, Laugraud A, Gautier C, Simon EB, Devaux F, Magdinier F, Vourc'h C, Gilson E

Half of the human genome consists of repetitive DNA sequences. Recent studies in various organisms highlight the role of chromatin regulation of repetitive DNA in gene regulation as well as in maintainance of chromosomes and genome integrity. Hence, repetitive DNA sequences might be potential "sensors" for chromatin changes associated with pathogenesis. Therefore, we developed a new genomic tool called RepArray. RepArray is a repeat-specific microarray composed of a representative set of human repeated sequences including transposon-derived repeats, simple sequences repeats, tandemly repeated sequences such as centromeres and telomeres. We showed that combined to anti-methylcytosine immunoprecipitation assay, the RepArray can be used to generate repeat-specific methylation maps. Using cell lines impaired chemically or genetically for DNA methyltransferases activities, we were able to distinguish different epigenomes demonstrating that repeats can be used as markers of genome-wide methylation changes. Besides, using a well-documented system model, the thermal stress, we demonstrated that RepArray is also a fast and reliable tool to obtain an overview of overall transcriptional activity on whole repetitive compartment in a given cell type. Thus, the RepArray represents the first valuable tool for systematic and genome-wide analyses of the methylation and transcriptional status of the repetitive counterpart of the human genome.


57. Genetic differences among Staphylococcus aureus isolates from dairy ruminant species: a single-dye DNA microarray approach, Vet Microbiol 133, 105-114 (Pubmed: 18640795)
Vautor E, Magnone V, Rios G, Le Brigand K, Bergonier D, Lina G, Meugnier H, Barbry P, Thiery R, Pepin M

Staphylococcus aureus is recognized worldwide as a major pathogen causing clinical or subclinical intramammary infections in lactating sheep, goats and cows. The present study was carried out to compare 65 S. aureus isolates mainly obtained from nasal carriage and subclinical mastitis in dairy sheep and 43 isolates obtained from subclinical mastitis from 22 goats and 21 cows. A DNA microarray, containing probes against 190 true or putative virulence factors, was used to detect the presence of the virulence genes. Their presence/absence was independently assessed by PCR for the genes of interest. Sheep isolates obtained from the nostrils or the udders did not show any significant tissue specific virulence factor. The dominant pulse-field electrophoresis profile (OV/OV'), associated with spa clonal complex spa-CC 1773, matched mainly with the agr group III and was only found in ovine and caprine isolates. This clone was more specifically characterized by the prevalence of the following virulence genes: lpl4, ssl6, bsaA1, bsaB, bsaP, SAV0812. Moreover, seven virulence-associated genes (lpl1, sel, sec, tst, lukF-PV-like component, lukM, SAV0876) were associated with isolates from small ruminants, while the egc cluster, fhuD1, abiF and SAV2496 with bovine isolates. This genomic study suggests the existence of lineage- and host-specific genes leading to the development of host-specific pathogenic traits of S. aureus isolates.


58. Gene expression profiling of human liver transplants identifies an early transcriptional signature associated with initial poor graft function., Am J Transplant. 2008 Jun;8(6):1221-36. (Pubmed: 18522548)
Defamie V, Cursio R, Le Brigand K, Moreilhon C, Saint-Paul MC, Laurens M, Crenesse D, Cardinaud B, Auberger P, Gugenheim J, Barbry P, Mari B

Liver ischemia-reperfusion injury occurring in orthotopic liver transplantation (OLT) may be responsible for early graft failure. Molecular mechanisms underlying initial poor graft function (IPGF) have been poorly documented in human. The purpose of this study was to identify the major transcriptional alterations occurring in human livers during OLT. Twenty-one RNA extracts derived from liver transplant biopsies taken after graft reperfusion were compared with 7 RNA derived from normal control livers. Three hundred seventy-one genes were significantly modulated and classified in molecular pathways relevant to liver metabolism, inflammatory response, cell proliferation and liver protection. Grafts were then subdivided into two groups based on their peak levels of serum aspartate amino transferase within 72 h after OLT (group 1, non-IPGF: 14 patients; group 2, IPGF: 7 patients). The two corresponding data sets were compared using a supervised prediction method. A new set of genes able to correctly classify 71% of the patients was defined. These genes were functionally associated with oxidative stress, inflammation and inhibition of cell proliferation. This study provides a comprehensive picture of the transcriptional events associated with human OLT and IPGF. We anticipate that such alterations provide a framework for the elucidation of the molecular mechanisms leading to IPGF.


59. Gene expression profiling in human gastric mucosa infected with Helicobacter pylori., Mod Pathol. 2007 Sep;20(9):974-89. Epub 2007 Jul 20. (Pubmed: 17643099)
Hofman VJ, Moreilhon C, Brest PD, Lassalle S, Le Brigand K, Sicard D, Raymond J, Lamarque D, Hebuterne XA, Mari B, Barbry PJ, Hofman PM

Pathogenic mechanisms associated with Helicobacter pylori infection enhance susceptibility of the gastric epithelium to carcinogenic conversion. We have characterized the gene expression profiles of gastric biopsies from 69 French Caucasian patients, of which 43 (62%) were infected with H. pylori. The bacterium was detected in 27 of the 42 antral biopsies examined and in 16 of the 27 fundic biopsies. Infected biopsies were selected for the presence of chronic active gastritis, in absence of metaplasia and dysplasia of the gastric mucosa. Infected antral and fundic biopsies exhibited distinct transcriptional responses. Altered responses were linked with: (1) the extent of polymorphonuclear leukocyte infiltration, (2) bacterial density, and (3) the presence of the virulence factors vacA, babA2, and cagA. Robust modulation of transcripts associated with Toll-like receptors, signal transduction, the immune response, apoptosis, and the cell cycle was consistent with expected responses to Gram-negative bacterial infection. Altered expression of interferon-regulated genes (IFITM1, IRF4, STAT6), indicative of major histocompatibility complex (MHC) II-mediated and Th1-specific responses, as well as altered expression of GATA6, have previously been described in precancerous states. Upregulation of genes abundantly expressed in cancer tissues (UBD, CXCL13, LY96, MAPK8, MMP7, RANKL, CCL18) or in stem cells (IFITM1 and WFDC2) may reveal a molecular switch towards a premalignant state in infected tissues. Tissue microarray analysis of a large number of biopsies, which were either positive or negative for the cag-A virulence factor, when compared to each other and to noninfected controls, confirmed observed gene alterations at the protein level, for eight key transcripts. This study provides 'proof-of-principle' data for identifying molecular mechanisms driving H. pylori-associated carcinogenesis before morphological evidence of changes along the neoplastic progression pathway.Modern Pathology (2007) 20, 974-989; doi:10.1038/modpathol.3800930; published online 20 July 2007.


60. Transcriptional signature of epidermal keratinocytes subjected to in vitro scratch wounding, reveals selective roles for ERK1/2, P38 and PI3K signalling pathways., J Biol Chem. 2007 May 18;282(20):15090-102. (Pubmed: 17363378)
Fitsialos G, Chassot AA, Turchi L, Dayem MA, Lebrigand K, Moreilhon C, Meneguzzi G, Busca R, Mari B, Barbry P, Ponzio G

Covering denuded dermal surfaces after injury requires migration, proliferation and differentiation of skin keratinocytes. To clarify the major traits controlling these intermingled biological events, we surveyed the genomic modifications occurring during the course of a scratch wound closure of cultured human keratinocytes. Using a DNA microarray approach, we report the identification of 161 new markers of epidermal repair. Expression data, combined with functional analysis performed with specific inhibitors of ERK, p38[MAPK] and PI3 kinases, demonstrate that kinase pathways exert very selective functions by precisely controlling the expression of specific genes. Inhibition of the ERK pathway totally blocks the wound closure and inactivates many early transcription factors and EGF-type growth factors. P38[MAPK] inhibition only delays "healing", probably in line with the control of genes involved in the propagation of injury-initiated signalling. In contrast, PI3 kinase inhibition accelerates the scratch closure and potentiates the scratch-dependent stimulation of three genes related to epithelial cell transformation, namely HAS3, HBEGF and Ets1. Our results define in vitro human keratinocyte wound closure as a repair process resulting from a fine balance between positive signals controlled by ERK and p38[MAPK], and negative ones triggered by PI3 kinase. The perturbation of any of these pathways might lead to dysfunction in the healing process, similar to those observed in pathological wounding phenotypes, such as hypertrophic scars or keloids.


61. Mediante: a Web-based Microarray Data manager, Bioinformatics, doi:10.1093/bioinformatics/btm106 (Pubmed: 16855282)
Le Brigand K, Barbry P

Summary : Mediante is a MIAME-compliant microarray data manager that links together annotations and experimental data. Developed as a J2EE three-tier application, Mediante integrates a management system for production of long oligonucleotide microarrays, an experimental data repository suitable for home made or commercial microarrays, and a user interface dedicated to the management of microarrays projects. Several tools allow quality control of hybridizations and submission of validated data to public repositories.
Availability : http://www.microarray.fr.
Supplementary data : http://www.microarray.fr/SP/lebrigand2007/


62. Suppression of microRNA-silencing pathway by HIV-1 during virus replication., Science. 2007 Mar 16;315(5818):1579-82. (Pubmed: 17322031)
Triboulet R, Mari B, Lin YL, Chable-Bessia C, Bennasser Y, Lebrigand K, Cardinaud B, Maurin T, Barbry P, Baillat V, Reynes J, Corbeau P, Jeang KT, Benkirane M

MicroRNAs (miRNAs) are single-stranded noncoding RNAs of 19 to 25 nucleotides that function as gene regulators and as a host cell defense against both RNA and DNA viruses. We provide evidence for a physiological role of the miRNA-silencing machinery in controlling HIV-1 replication. Type III RNAses Dicer and Drosha, responsible for miRNA processing, inhibited virus replication both in peripheral blood mononuclear cells from HIV-1-infected donors and in latently infected cells. In turn, HIV-1 actively suppressed the expression of the polycistronic miRNA cluster miR-17/92. This suppression was found to be required for efficient viral replication and was dependent on the histone acetyltransferase Tat cofactor PCAF. Our results highlight the involvement of the miRNA-silencing pathway in HIV-1 replication and latency.


63. An open-access long oligonucleotide microarray resource for analysis of the human and mouse transcriptomes., Nucleic Acids Res. 2006 Jul 19;34(12):e87 (Pubmed: 17384016)
Le Brigand K, Russell R, Moreilhon C, Rouillard JM, Jost B, Amiot F, Magnone V, Bole-Feysot C, Rostagno P, Virolle V, Defamie V, Dessen P, Williams G, Lyons P, Rios G, Mari B, Gulari E, Kastner P, Gidrol X, Freeman TC, Barbry P

Two collections of oligonucleotides have been designed for preparing pangenomic human and mouse microarrays. A total of 148,993 and 121,703 oligonucleotides were designed against human and mouse transcripts. Quality scores were created in order to select 25,342 human and 24,109 mouse oligonucleotides. They correspond to: (i) a BLAST-specificity score; (ii) the number of expressed sequence tags matching each probe; (iii) the distance to the 3' end of the target mRNA. Scores were also used to compare in silico the two microarrays with commercial microarrays. The sets described here, called RNG/MRC collections, appear at least as specific and sensitive as those from the commercial platforms. The RNG/MRC collections have now been used by an Anglo-French consortium to distribute more than 3500 microarrays to the academic community. Ad hoc identification of tissue-specific transcripts and a approximately 80% correlation with hybridizations performed on Affymetrix GeneChiptrade mark suggest that the RNG/MRC microarrays perform well. This work provides a comprehensive open resource for investigators working on human and mouse transcriptomes, as well as a generic method to generate new microarray collections in other organisms. All information related to these probes, as well as additional information about commercial microarrays have been stored in a freely-accessible database called MEDIANTE.