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

Paquet Agnes

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).

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.

Light chain (LC) deposition disease (LCDD) is a rare disorder characterized by glomerular and peritubular amorphous deposits of a monoclonal immunoglobulin LC, leading to nodular glomerulosclerosis and nephrotic syndrome. We developed a transgenic model using site-directed insertion of the variable domain of a pathogenic human LC gene into the mouse immunoglobulin κ locus, ensuring its production by all plasma cells (PCs). High free LC levels were achieved after backcrossing with mice presenting increased PC differentiation and no immunoglobulin heavy chain production. Our mouse model recapitulates the characteristic features of LCDD, including progressive glomerulosclerosis, nephrotic-range proteinuria, and finally kidney failure. The variable domain of the LC bears alone the structural properties involved in its pathogenicity. RNA sequencing conducted on PCs demonstrated that LCDD LC induces endoplasmic reticulum stress, likely accounting for the high efficiency of proteasome inhibitor-based therapy. Accordingly, reduction of circulating pathogenic LC was efficiently achieved and not only preserved renal function but also partially reversed kidney lesions. Finally, transcriptome analysis of presclerotic glomeruli revealed that proliferation and extracellular matrix remodeling represented the first steps of glomerulosclerosis, paving the way for future therapeutic strategies in LCDD and other kidney diseases featuring diffuse glomerulosclerosis, particularly diabetic nephropathy.

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.

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.

Melanoma patients harboring the BRAFV600E mutation are treated with vemurafenib. Almost all of them ultimately acquire resistance, leading to disease progression. Here, we find that a small molecule from a marine sponge, panicein A hydroquinone (PAH), overcomes resistance of BRAFV600E melanoma cells to vemurafenib, leading to tumor elimination in corresponding human xenograft models in mice. We report the synthesis of PAH and demonstrate that this compound inhibits the drug efflux activity of the Hedgehog receptor, Patched. Our SAR study allowed identifying a key pharmacophore responsible for this activity. We showed that Patched is strongly expressed in metastatic samples from a cohort of melanoma patients and is correlated with decreased overall survival. Patched is a multidrug transporter that uses the proton motive force to efflux drugs. This makes its function specific to cancer cells, thereby avoiding toxicity issues that are commonly observed with inhibitors of ABC multidrug transporters. Our data provide strong evidence that PAH is a highly promising lead for the treatment of vemurafenib resistant BRAFV600E melanoma.

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.

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.

Diffuse large B cell lymphoma (DLBCL) is a heterogeneous disease treated with anti-CD20-based immuno-chemotherapy (R-CHOP). We identified that low levels of GAPDH predict a poor response to R-CHOP treatment. Importantly, we demonstrated that GAPDHlow lymphomas use OxPhos metabolism and rely on mTORC1 signaling and glutaminolysis. Consistently, disruptors of OxPhos metabolism (phenformin) or glutaminolysis (L-asparaginase) induce cytotoxic responses in GAPDHlow B cells and improve GAPDHlow B cell-lymphoma-bearing mice survival, while they are low or not efficient on GAPDHhigh B cell lymphomas. Ultimately, we selected four GAPDHlow DLBCL patients, who were refractory to all anti-CD20-based therapies, and targeted DLBCL metabolism using L-asparaginase (K), mTOR inhibitor (T), and metformin (M) (called KTM therapy). Three out of the four patients presented a complete response upon one cycle of KTM. These findings establish that the GAPDH expression level predicts DLBCL patients' response to R-CHOP treatment and their sensitivity to specific metabolic inhibitors.

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.

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.

Multiciliated cells (MCCs) harbor dozens to hundreds of motile cilia, which generate hydrodynamic forces important in animal physiology. In vertebrates, MCC differentiation involves massive centriole production by poorly characterized structures called deuterosomes. Here, single-cell RNA sequencing reveals that human deuterosome stage MCCs are characterized by the expression of many cell cycle-related genes. We further investigated the uncharacterized vertebrate-specific cell division cycle 20B (CDC20B) gene, which hosts microRNA-449abc. We show that CDC20B protein associates to deuterosomes and is required for centriole release and subsequent cilia production in mouse and Xenopus MCCs. CDC20B interacts with PLK1, a kinase known to coordinate centriole disengagement with the protease Separase in mitotic cells. Strikingly, over-expression of Separase rescues centriole disengagement and cilia production in CDC20B-deficient MCCs. This work reveals the shaping of deuterosome-mediated centriole production in vertebrate MCCs, by adaptation of canonical and recently evolved cell cycle-related molecules.

Fragile X syndrome (FXS), the most common form of inherited intellectual disability (ID) and a leading cause of autism, results from the loss of expression of the Fmr1 gene which encodes the RNA-binding protein Fragile X Mental Retardation Protein (FMRP). Among the thousands mRNA targets of FMRP, numerous encode regulators of ion homeostasis. It has also been described that FMRP directly interacts with Ca2+ channels modulating their activity. Collectively these findings suggest that FMRP plays critical roles in Ca2+ homeostasis during nervous system development. We carried out a functional analysis of Ca2+ regulation using a calcium imaging approach in Fmr1-KO cultured neurons and we show that these cells display impaired steady state Ca2+ concentration and an altered entry of Ca2+ after KCl-triggered depolarization. Consistent with these data, we show that the protein product of the Cacna1a gene, the pore-forming subunit of the Cav2.1 channel, is less expressed at the plasma membrane of Fmr1-KO neurons compared to wild-type (WT). Thus, our findings point out the critical role that Cav2.1 plays in the altered Ca2+ flux in Fmr1-KO neurons, impacting Ca2+ homeostasis of these cells. Remarkably, we highlight a new phenotype of cultured Fmr1-KO neurons that can be considered a novel cellular biomarker and is amenable to small molecule screening and identification of new drugs to treat FXS.

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.

Despite constant progress in prognostic risk stratification, children with acute myeloid leukemia (AML) still relapse. Treatment failure and subsequent relapse have been attributed to acute myeloid leukemia-initiating cells (LSC), which harbor stem cell properties and are inherently chemoresistant. Although pediatric and adult AML represent two genetically very distinct diseases, we reasoned that common LSC gene expression programs are shared and consequently, the highly prognostic LSC17 signature score recently developed in adults may also be of clinical interest in childhood AML. Here, we demonstrated prognostic relevance of the LSC17 score in pediatric non-core-binding factor AML using Nanostring technology (ELAM02) and RNA-seq data from the NCI (TARGET-AML). AML were stratified by LSC17 quartile groups (lowest 25%, intermediate 50% and highest 25%) and children with low LSC17 score had significantly better event-free survival (EFS: HR = 3.35 (95%CI = 1.64-6.82), P < 0.001) and overall survival (OS: HR = 3.51 (95%CI = 1.38-8.92), P = 0.008) compared with patients with high LSC17 scores. More importantly, the high LSC17 score was an independent negative EFS and OS prognosticator determined by multivariate Cox model analysis (EFS: HR = 3.42 (95% CI = 1.63-7.16), P = 0.001; OS HR = 3.02 (95%CI = 1.16-7.85), P = 0.026). In conclusion, we have demonstrated the broad applicability of the LSC17 score in the clinical management of AML by extending its prognostic relevance to pediatric AML.

OBJECTIVES: The effect of anti-PD-1/PD-L1 inhibitors on lung adenocarcinomas (LADCs) with KRAS mutations is debatable. We examined the association between specific mutant KRAS proteins and the immune infiltrates with the outcome of patients with LADCs. PATIENTS AND METHODS: In 219 LADCs harboring either wild-type (WT) or mutated KRAS gene, we quantified the density of several immune markers by immunohistochemistry followed by automated digital image analysis. Data were correlated to clinicopathological parameters and outcome of patients. RESULTS: Tumors harboring mutant KRAS-G12 V had a significantly higher PD-L1 expression compared to other tumors (p = 0.044), while mutant KRAS-G12D tumors showed an increase in the density of CD66b+ cells (p = 0.001). High PD-L1 expression in tumor cells was associated to improved overall survival (OS) in KRAS mutant patients (p = 0.012), but not in the WT population (p = 0.385), whereas increased PD-L1 expression in immune cells correlated to poor OS of KRAS-WT patients (p = 0.025), with no difference in patients with KRAS mutations. CONCLUSIONS: KRAS mutational status can affect the immune microenvironment and survival of LADC patients in a heterogeneous way, implying that specific mutant KRAS variants expressed by the tumor should be considered when stratifying patients for immunotherapy.

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.

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.

Although Tacrolimus is an immunosuppressive drug widely used in renal transplantation, its chronic use paradoxically induces nephrotoxic effects, in particular renal fibrosis, which is responsible for chronic allograft dysfunction and represents a major prognostic factor of allograft survival. As molecular pathways and mechanisms involved in Tacrolimus-induced fibrogenic response are poorly elucidated, we assessed whether miRNAs are involved in the nephrotoxic effects mediated by Tacrolimus. Treatment of CD-1 mice with Tacrolimus (1 mg/kg/d for 28 days) resulted in kidney injury and was associated with alteration of a gene expression signature associated with cellular stress, fibrosis and inflammation. Tacrolimus also affected renal miRNA expression, including miRNAs previously involved in fibrotic and inflammatory processes as "fibromirs" such as miR-21-5p, miR-199a-5p and miR-214-3p. In agreement with in vivo data, Renal Proximal Tubular Epithelial cells exposed to Tacrolimus (25 and 50 µM) showed upregulation of miR-21-5p and the concomitant induction of epithelial phenotypic changes, inflammation and oxidative stress. In conclusion, this study suggests for the first time that miRNAs, especially fibromiRs, participate to Tacrolimus-induced nephrotoxic effects. Therefore, targeting miRNAs may be a new therapeutic option to counteract Tacrolimus deleterious effects on kidney.

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.

In therapeutic research, the safety and efficacy of pharmaceutical products are necessarily tested on humans via clinical trials after an extensive and expensive preclinical development period. Methodologies such as computer modeling and clinical trial simulation (CTS) might represent a valuable option to reduce animal and human assays. The relevance of these methods is well recognized in pharmacokinetics and pharmacodynamics from the preclinical phase to postmarketing. However, they are barely used and are poorly regarded for drug approval, despite Food and Drug Administration and European Medicines Agency recommendations. The generalization of CTS could be greatly facilitated by the availability of software for modeling biological systems, by clinical trial studies and hospital databases. Data sharing and data merging raise legal, policy and technical issues that will need to be addressed. Development of future molecules will have to use CTS for faster development and thus enable better patient management. Drug activity modeling coupled with disease modeling, optimal use of medical data and increased computing speed should allow this leap forward. The realization of CTS requires not only bioinformatics tools to allow interconnection and global integration of all clinical data but also a universal legal framework to protect the privacy of every patient. While recognizing that CTS can never replace 'real-life' trials, they should be implemented in future drug development schemes to provide quantitative support for decision-making. This in silico medicine opens the way to the P4 medicine: predictive, preventive, personalized and participatory.

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.

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.

Recent advances in genomic technologies have revolutionized acute myeloid leukemia (AML) understanding by identifying potential novel actionable genomic alterations. Consequently, current risk stratification at diagnosis not only relies on cytogenetics, but also on the inclusion of several of these abnormalities. Despite this progress, AML remains a heterogeneous and complex malignancy with variable response to current therapy. Although copy-number alterations (CNAs) are accepted prognostic markers in cancers, large-scale genomic studies aiming at identifying specific prognostic CNA-based markers in AML are still lacking. Using 367 AML, we identified four recurrent CNA on chromosomes 11 and 21 that predicted outcome even after adjusting for standard prognostic risk factors and potentially delineated two new subclasses of AML with poor prognosis. ERG amplification, the most frequent CNA, was related to cytarabine resistance, a cornerstone drug of AML therapy. These findings were further validated in The Cancer Genome Atlas data. Our results demonstrate that specific CNA are of independent prognostic relevance, and provide new molecular information into the genomic basis of AML and cytarabine response. Finally, these CNA identified two potential novel risk groups of AML, which when confirmed prospectively, may improve the clinical risk stratification and potentially the AML outcome.Leukemia advance online publication, 4 November 2016; doi:10.1038/leu.2016.265.

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.

The SENS-IS test protocol for the in vitro detection of sensitizers is based on a reconstructed human skin model (Episkin) as the test system and on the analysis of the expression of a large panel of genes. Its excellent performance was initially demonstrated with a limited set of test chemicals. Further studies (described here) were organized to confirm these preliminary results and to obtain a detailed statistical analysis of the predictive capacity of the assay. A ring-study was thus organized and performed within three laboratories, using a test set of 19 blind coded chemicals. Data analysis indicated that the assay is robust, easily transferable and offers high predictivity and excellent within- and between-laboratories reproducibility. To further evaluate the predictivity of the test protocol according to Cooper statistics a comprehensive test set of 150 chemicals was then analyzed. Again, data analysis confirmed the excellent capacity of the SENS-IS assay for predicting both hazard and potency characteristics, confirming that this assay should be considered as a serious alternative to the available in vivo sensitization tests.

BACKGROUND: Recent efforts in HIV-1 vaccine design have focused on immunogens that evoke potent neutralizing antibody responses to a broad spectrum of viruses circulating worldwide. However, the development of effective vaccines will depend on the identification and characterization of the neutralizing antibodies and their epitopes. We developed bioinformatics methods to predict epitope networks and antigenic determinants using structural information, as well as corresponding genotypes and phenotypes generated by a highly sensitive and reproducible neutralization assay.282 clonal envelope sequences from a multiclade panel of HIV-1 viruses were tested in viral neutralization assays with an array of broadly neutralizing monoclonal antibodies (mAbs: b12, PG9,16, PGT121 - 128, PGT130 - 131, PGT135 - 137, PGT141 - 145, and PGV04). We correlated IC50 titers with the envelope sequences, and used this information to predict antibody epitope networks. Structural patches were defined as amino acid groups based on solvent-accessibility, radius, atomic depth, and interaction networks within 3D envelope models. We applied a boosted algorithm consisting of multiple machine-learning and statistical models to evaluate these patches as possible antibody epitope regions, evidenced by strong correlations with the neutralization response for each antibody. RESULTS: We identified patch clusters with significant correlation to IC50 titers as sites that impact neutralization sensitivity and therefore are potentially part of the antibody binding sites. Predicted epitope networks were mostly located within the variable loops of the envelope glycoprotein (gp120), particularly in V1/V2. Site-directed mutagenesis experiments involving residues identified as epitope networks across multiple mAbs confirmed association of these residues with loss or gain of neutralization sensitivity. CONCLUSIONS: Computational methods were implemented to rapidly survey protein structures and predict epitope networks associated with response to individual monoclonal antibodies, which resulted in the identification and deeper understanding of immunological hotspots targeted by broadly neutralizing HIV-1 antibodies.

BACKGROUND: Drug resistance testing and co-receptor tropism determination are key components of the management of antiretroviral therapy for HIV-1-infected individuals. The purpose of this study was to examine trends of HIV-1 resistance and viral evolution in the past decade by surveying a large commercial patient testing database. METHODS: Temporal trends of drug resistance, viral fitness and co-receptor usage among samples submitted for routine phenotypic and genotypic resistance testing to protease inhibitors (PIs), nucleoside reverse transcriptase inhibitors (NRTIs) and non-nucleoside reverse transcriptase inhibitors (NNRTIs), as well as for tropism determination were investigated. RESULTS: Within 62,397 resistant viruses reported from 2003 to 2012, we observed a decreasing trend in the prevalence of three-class resistance (from 25% to 9%) driven by decreased resistance to PIs (43% to 21%) and NRTIs (79% to 57%), while observing a slight increase in NNRTI resistance (68% to 75%). The prevalence of CXCR4-mediated entry among tropism testing samples (n=52,945) declined over time from 47% in 2007 to 40% in 2012. A higher proportion of CXCR4-tropic viruses was observed within samples with three-class resistance (50%) compared with the group with no resistance (36%). CONCLUSIONS: Decreased prevalence of three-class resistance and increased prevalence of one-class resistance was observed within samples reported between 2003 and 2012. The fraction of CXCR4-tropic viruses has decreased over time; however, CXCR4 usage was more prevalent among multi-class-resistant samples, which may be due to the more advanced disease stage of treatment-experienced patients. These trends have important implications for clinical practice and future drug discovery and development.