Smoking changes adaptive immunity with persistent effects.

Individuals differ widely in their immune responses, with age, sex and genetic factors having major roles in this inherent variability1-6. However, the variables that drive such differences in cytokine secretion-a crucial component of the host response to immune challenges-remain poorly defined. Here we investigated 136 variables and identified smoking, cytomegalovirus latent infection and body mass index as major contributors to variability in cytokine response, with effects of comparable magnitudes with age, sex and genetics. We find that smoking influences both innate and adaptive immune responses. Notably, its effect on innate responses is quickly lost after smoking cessation and is specifically associated with plasma levels of CEACAM6, whereas its effect on adaptive responses persists long after individuals quit smoking and is associated with epigenetic memory. This is supported by the association of the past smoking effect on cytokine responses with DNA methylation at specific signal trans-activators and regulators of metabolism. Our findings identify three novel variables associated with cytokine secretion variability and reveal roles for smoking in the short- and long-term regulation of immune responses. These results have potential clinical implications for the risk of developing infections, cancers or autoimmune diseases.

BioConvert: a comprehensive format converter for life sciences.

Bioinformatics is a field known for the numerous standards and formats that have been developed over the years. This plethora of formats, sometimes complementary, and often redundant, poses many challenges to bioinformatics data analysts. They constantly need to find the best tool to convert their data into the suitable format, which is often a complex, technical and time consuming task. Moreover, these small yet important tasks are often difficult to make reproducible. To overcome these difficulties, we initiated BioConvert, a collaborative project to facilitate the conversion of life science data from one format to another. BioConvert aggregates existing software within a single framework and complemented them with original code when needed. It provides a common interface to make the user experience more streamlined instead of having to learn tens of them. Currently, BioConvert supports about 50 formats and 100 direct conversions in areas such as alignment, sequencing, phylogeny, and variant calling. In addition to being useful for end-users, BioConvert can also be utilized by developers as a universal benchmarking framework for evaluating and comparing numerous conversion tools. Additionally, we provide a web server implementing an online user-friendly interface to BioConvert, hence allowing direct use for the community.

Excessive self-grooming, gene dysregulation and imbalance between the striosome and matrix compartments in the striatum of Shank3 mutant mice.

Autism is characterized by atypical social communication and stereotyped behaviors. Mutations in the gene encoding the synaptic scaffolding protein SHANK3 are detected in 1-2% of patients with autism and intellectual disability, but the mechanisms underpinning the symptoms remain largely unknown. Here, we characterized the behavior of Shank3 Δ11/Δ11 mice from 3 to 12 months of age. We observed decreased locomotor activity, increased stereotyped self-grooming and modification of socio-sexual interaction compared to wild-type littermates. We then used RNAseq on four brain regions of the same animals to identify differentially expressed genes (DEGs). DEGs were identified mainly in the striatum and were associated with synaptic transmission (e.g., Grm2, Dlgap1), G-protein-signaling pathways (e.g., Gnal, Prkcg1, and Camk2g), as well as excitation/inhibition balance (e.g., Gad2). Downregulated and upregulated genes were enriched in the gene clusters of medium-sized spiny neurons expressing the dopamine 1 (D1-MSN) and the dopamine 2 receptor (D2-MSN), respectively. Several DEGs (Cnr1, Gnal, Gad2, and Drd4) were reported as striosome markers. By studying the distribution of the glutamate decarboxylase GAD65, encoded by Gad2, we showed that the striosome compartment of Shank3 Δ11/Δ11 mice was enlarged and displayed much higher expression of GAD65 compared to wild-type mice. Altogether, these results indicate altered gene expression in the striatum of Shank3-deficient mice and strongly suggest, for the first time, that the excessive self-grooming of these mice is related to an imbalance in the striatal striosome and matrix compartments.

IRF8 regulates efficacy of therapeutic anti-CD20 monoclonal antibodies.

Anti-CD20 monoclonal antibodies such as Rituximab, Ofatumumab, and Obinutuzumab are widely used to treat lymphomas and autoimmune diseases. They act by depleting B cells, mainly through Fc-dependent effectors functions. Some patients develop resistance to treatment but the underlying mechanisms are poorly understood. Here, we performed a genome-wide CRISPR/Cas9 screen to identify genes regulating the efficacy of anti-CD20 antibodies. We used as a model the killing of RAJI B cells by Rituximab through complement-dependent-cytotoxicity (CDC). As expected, the screen identified MS4A1, encoding CD20, the target of Rituximab. Among other identified genes, the role of Interferon Regulatory Factor 8 (IRF8) was validated in two B-cell lines. IRF8 KO also decreased the efficacy of antibody-dependent cellular cytotoxicity and phagocytosis (ADCC and ADCP) induced by anti-CD20 antibodies. We further show that IRF8 is necessary for efficient CD20 transcription. Levels of IRF8 and CD20 RNA or proteins correlated in normal B cells and in hundreds of malignant B cells. Therefore, IRF8 regulates CD20 expression and controls the depleting capacity of anti-CD20 antibodies. Our results bring novel insights into the pathways underlying resistance to CD20-targeting immunotherapies.

Erythro-myeloid progenitor origin of Hofbauer cells in the early mouse placenta.

Hofbauer cells (HBCs) are tissue macrophages of the placenta thought to be important for fetoplacental vascular development and innate immune protection. The developmental origins of HBCs remain unresolved and could implicate functional diversity of HBCs in placenta development and disease. In this study, we used flow cytometry and paternally inherited reporters to phenotype placenta macrophages and to identify fetal-derived HBCs and placenta-associated maternal macrophages in the mouse. In vivo pulse-labeling traced the ontogeny of HBCs from yolk sac-derived erythro-myeloid progenitors, with a minor contribution from fetal hematopoietic stem cells later on. Single-cell RNA-sequencing revealed transcriptional similarities between placenta macrophages and erythro-myeloid progenitor-derived fetal liver macrophages and microglia. As with other fetal tissue macrophages, HBCs were dependent on the transcription factor Pu.1, the loss-of-function of which in embryos disrupted fetoplacental labyrinth morphology, supporting a role for HBC in labyrinth angiogenesis and/or remodeling. HBC were also sensitive to Pu.1 (Spi1) haploinsufficiency, which caused an initial deficiency in the numbers of macrophages in the early mouse placenta. These results provide groundwork for future investigation into the relationship between HBC ontogeny and function in placenta pathophysiology.

Massively parallel analysis of human 3' UTRs reveals that AU-rich element length and registration predict mRNA destabilization.

AU-rich elements (AREs) are 3' UTR cis-regulatory elements that regulate the stability of mRNAs. Consensus ARE motifs have been determined, but little is known about how differences in 3' UTR sequences that conform to these motifs affect their function. Here, we use functional annotation of sequences from 3' UTRs (fast-UTR), a massively parallel reporter assay (MPRA), to investigate the effects of 41,288 3' UTR sequence fragments from 4653 transcripts on gene expression and mRNA stability in Jurkat and Beas2B cells. Our analyses demonstrate that the length of an ARE and its registration (the first and last nucleotides of the repeating ARE motif) have significant effects on gene expression and stability. Based on this finding, we propose improved ARE classification and concomitant methods to categorize and predict the effect of AREs on gene expression and stability. Finally, to investigate the advantages of our general experimental design we examine other motifs including constitutive decay elements (CDEs), where we show that the length of the CDE stem-loop has a significant impact on steady-state expression and mRNA stability. We conclude that fast-UTR, in conjunction with our analytical approach, can produce improved yet simple sequence-based rules for predicting the activity of human 3' UTRs.

A mouse-specific retrotransposon drives a conserved Cdk2ap1 isoform essential for development.

Retrotransposons mediate gene regulation in important developmental and pathological processes. Here, we characterized the transient retrotransposon induction during preimplantation development of eight mammals. Induced retrotransposons exhibit similar preimplantation profiles across species, conferring gene regulatory activities, particularly through long terminal repeat (LTR) retrotransposon promoters. A mouse-specific MT2B2 retrotransposon promoter generates an N-terminally truncated Cdk2ap1ΔN that peaks in preimplantation embryos and promotes proliferation. In contrast, the canonical Cdk2ap1 peaks in mid-gestation and represses cell proliferation. This MT2B2 promoter, whose deletion abolishes Cdk2ap1ΔN production, reduces cell proliferation and impairs embryo implantation, is developmentally essential. Intriguingly, Cdk2ap1ΔN is evolutionarily conserved in sequence and function yet is driven by different promoters across mammals. The distinct preimplantation Cdk2ap1ΔN expression in each mammalian species correlates with the duration of its preimplantation development. Hence, species-specific transposon promoters can yield evolutionarily conserved, alternative protein isoforms, bestowing them with new functions and species-specific expression to govern essential biological divergence.

A high-risk retinoblastoma subtype with stemness features, dedifferentiated cone states and neuronal/ganglion cell gene expression.

Retinoblastoma is the most frequent intraocular malignancy in children, originating from a maturing cone precursor in the developing retina. Little is known on the molecular basis underlying the biological and clinical behavior of this cancer. Here, using multi-omics data, we demonstrate the existence of two retinoblastoma subtypes. Subtype 1, of earlier onset, includes most of the heritable forms. It harbors few genetic alterations other than the initiating RB1 inactivation and corresponds to differentiated tumors expressing mature cone markers. By contrast, subtype 2 tumors harbor frequent recurrent genetic alterations including MYCN-amplification. They express markers of less differentiated cone together with neuronal/ganglion cell markers with marked inter- and intra-tumor heterogeneity. The cone dedifferentiation in subtype 2 is associated with stemness features including low immune and interferon response, E2F and MYC/MYCN activation and a higher propensity for metastasis. The recognition of these two subtypes, one maintaining a cone-differentiated state, and the other, more aggressive, associated with cone dedifferentiation and expression of neuronal markers, opens up important biological and clinical perspectives for retinoblastomas.

miR-200 deficiency promotes lung cancer metastasis by activating Notch signaling in cancer-associated fibroblasts.

Lung adenocarcinoma, the most prevalent lung cancer subtype, is characterized by its high propensity to metastasize. Despite the importance of metastasis in lung cancer mortality, its underlying cellular and molecular mechanisms remain largely elusive. Here, we identified miR-200 miRNAs as potent suppressors for lung adenocarcinoma metastasis. miR-200 expression is specifically repressed in mouse metastatic lung adenocarcinomas, and miR-200 decrease strongly correlates with poor patient survival. Consistently, deletion of mir-200c/141 in the KrasLSL-G12D/+ ; Trp53flox/flox lung adenocarcinoma mouse model significantly promoted metastasis, generating a desmoplastic tumor stroma highly reminiscent of metastatic human lung cancer. miR-200 deficiency in lung cancer cells promotes the proliferation and activation of adjacent cancer-associated fibroblasts (CAFs), which in turn elevates the metastatic potential of cancer cells. miR-200 regulates the functional interaction between cancer cells and CAFs, at least in part, by targeting Notch ligand Jagged1 and Jagged2 in cancer cells and inducing Notch activation in adjacent CAFs. Hence, the interaction between cancer cells and CAFs constitutes an essential mechanism to promote metastatic potential.

Megakaryocyte production is sustained by direct differentiation from erythromyeloid progenitors in the yolk sac until midgestation.

The extra-embryonic yolk sac contains the first definitive multipotent hematopoietic cells, denominated erythromyeloid progenitors. They originate in situ prior to the emergence of hematopoietic stem cells and give rise to erythroid, monocytes, granulocytes, mast cells and macrophages, the latter in a Myb transcription factor-independent manner. We uncovered here the heterogeneity of yolk sac erythromyeloid progenitors, at the single cell level, and discriminated multipotent from committed progenitors, prior to fetal liver colonization. We identified two temporally distinct megakaryocyte differentiation pathways. The first occurs in the yolk sac, bypasses intermediate bipotent megakaryocyte-erythroid progenitors and, similar to the differentiation of macrophages, is Myb-independent. By contrast, the second originates later, from Myb-dependent bipotent progenitors expressing Csf2rb and colonize the fetal liver, where they give rise to megakaryocytes and to large numbers of erythrocytes. Understanding megakaryocyte development is crucial as they play key functions during vascular development, in particular in separating blood and lymphatic networks.

The antibiotic bedaquiline activates host macrophage innate immune resistance to bacterial infection.

Antibiotics are widely used in the treatment of bacterial infections. Although known for their microbicidal activity, antibiotics may also interfere with the host's immune system. Here, we analyzed the effects of bedaquiline (BDQ), an inhibitor of the mycobacterial ATP synthase, on human macrophages. Genome-wide gene expression analysis revealed that BDQ reprogramed cells into potent bactericidal phagocytes. We found that 579 and 1,495 genes were respectively differentially expressed in naive- and M. tuberculosis-infected macrophages incubated with the drug, with an over-representation of lysosome-associated genes. BDQ treatment triggered a variety of antimicrobial defense mechanisms, including phagosome-lysosome fusion, and autophagy. These effects were associated with activation of transcription factor EB, involved in the transcription of lysosomal genes, resulting in enhanced intracellular killing of different bacterial species that were naturally insensitive to BDQ. Thus, BDQ could be used as a host-directed therapy against a wide range of bacterial infections.

The discovery of antibiotic drugs, which treat diseases caused by bacteria, has been a hugely valuable advance in modern medicine. They work by targeting specific cellular processes in bacteria, ultimately stopping them from multiplying or killing them outright. Antibiotics sometimes also affect their human hosts and can cause side-effects, such as gut problems or skin reactions. Recent evidence suggests that antibiotics also have an impact on the human immune system. This may happen either indirectly, by affecting 'friendly' bacteria normally present in the body, or through direct effects on immune cells. In turn, this could change the effectiveness of drug treatments. For example, if an antibiotic weakens immune cells, the body could have difficulty fighting off the existing infection ­ or become more vulnerable to new ones. However, even though new drugs are being introduced to combat the worldwide rise of antibiotic-resistant bacteria, their effects on immunity are still not well understood. For example, bedaquiline is an antibiotic recently developed to treat tuberculosis infections that are resistant to several drugs. Giraud-Gatineau et al. wanted to determine if bedaquiline altered the human immune response to bacterial infection independently from its direct anti-microbial effects. Macrophages engulf foreign particles like bacteria and break them down using enzymes stored within small internal compartments, or 'lysosomes'. Initial experiments using human macrophages, grown both with and without bedaquiline, showed that the drug did not harm the cells and that they grew normally. A combination of microscope imaging and genetic analysis revealed that exposure to bedaquiline not only increased the number of lysosomes within macrophage cells, but also the activity of genes and proteins that increase lysosomes' ability to break down foreign particles. These results suggested that bedaquiline treatment might make macrophages better at fighting infection, even if the drug itself had no direct effect on bacterial cells. Further studies, where macrophages were first treated with bedaquiline and then exposed to different types of bacteria known to be resistant to the drug, confirmed this hypothesis: in every case, the treated macrophages became efficient bacterial killers. In contrast, older anti-tuberculosis drugs did not have any such potentiating effect on the macrophages. This work sheds new light on our how antibiotic drugs can interact with the cells of the human immune system, and can sometimes even boost our innate defences. Such immune-boosting effects could one day be exploited to make more effective treatments against bacterial infections.

Impact of the gut microbiota on the m6A epitranscriptome of mouse cecum and liver.

The intestinal microbiota modulates host physiology and gene expression via mechanisms that are not fully understood. Here we examine whether host epitranscriptomic marks are affected by the gut microbiota. We use methylated RNA-immunoprecipitation and sequencing (MeRIP-seq) to identify N6-methyladenosine (m6A) modifications in mRNA of mice carrying conventional, modified, or no microbiota. We find that variations in the gut microbiota correlate with m6A modifications in the cecum, and to a lesser extent in the liver, affecting pathways related to metabolism, inflammation and antimicrobial responses. We analyze expression levels of several known writer and eraser enzymes, and find that the methyltransferase Mettl16 is downregulated in absence of a microbiota, and one of its target mRNAs, encoding S-adenosylmethionine synthase Mat2a, is less methylated. We furthermore show that Akkermansia muciniphila and Lactobacillus plantarum affect specific m6A modifications in mono-associated mice. Our results highlight epitranscriptomic modifications as an additional level of interaction between commensal bacteria and their host.

Polygenic Architecture of Human Neuroanatomical Diversity.

We analyzed the genomic architecture of neuroanatomical diversity using magnetic resonance imaging and single nucleotide polymorphism (SNP) data from >26 000 individuals from the UK Biobank project and 5 other projects that had previously participated in the ENIGMA (Enhancing NeuroImaging Genetics through Meta-Analysis) consortium. Our results confirm the polygenic architecture of neuroanatomical diversity, with SNPs capturing from 40% to 54% of regional brain volume variance. Chromosomal length correlated with the amount of phenotypic variance captured, r ~ 0.64 on average, suggesting that at a global scale causal variants are homogeneously distributed across the genome. At a local scale, SNPs within genes (~51%) captured ~1.5 times more genetic variance than the rest, and SNPs with low minor allele frequency (MAF) captured less variance than the rest: the 40% of SNPs with MAF <5% captured

Independent Component Analysis for Unraveling the Complexity of Cancer Omics Datasets.

Independent component analysis (ICA) is a matrix factorization approach where the signals captured by each individual matrix factors are optimized to become as mutually independent as possible. Initially suggested for solving source blind separation problems in various fields, ICA was shown to be successful in analyzing functional magnetic resonance imaging (fMRI) and other types of biomedical data. In the last twenty years, ICA became a part of the standard machine learning toolbox, together with other matrix factorization methods such as principal component analysis (PCA) and non-negative matrix factorization (NMF). Here, we review a number of recent works where ICA was shown to be a useful tool for unraveling the complexity of cancer biology from the analysis of different types of omics data, mainly collected for tumoral samples. Such works highlight the use of ICA in dimensionality reduction, deconvolution, data pre-processing, meta-analysis, and others applied to different data types (transcriptome, methylome, proteome, single-cell data). We particularly focus on the technical aspects of ICA application in omics studies such as using different protocols, determining the optimal number of components, assessing and improving reproducibility of the ICA results, and comparison with other popular matrix factorization techniques. We discuss the emerging ICA applications to the integrative analysis of multi-level omics datasets and introduce a conceptual view on ICA as a tool for defining functional subsystems of a complex biological system and their interactions under various conditions. Our review is accompanied by a Jupyter notebook which illustrates the discussed concepts and provides a practical tool for applying ICA to the analysis of cancer omics datasets.

Broad-Range Papillomavirus Transcriptome as a Biomarker of Papillomavirus-Associated Cervical High-Grade Cytology.

Human papillomaviruses (HPVs) are responsible for >99% of cervical cancers. Molecular diagnostic tests based on the detection of viral DNA or RNA have low positive predictive values for the identification of cancer or precancerous lesions. Triage with the Papanicolaou test lacks sensitivity; and even when combined with molecular detection of high-risk HPV, this results in a significant number of unnecessary colposcopies. We have developed a broad-range detection test of HPV transcripts to take a snapshot of the transcriptome of 16 high-risk or putative high-risk HPVs in cervical lesions (HPVs 16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, 66, 68, 73, and 82). The purpose of this novel molecular assay, named HPV RNA-Seq, is to detect and type HPV-positive samples and to determine a combination of HPV reads at certain specific viral spliced junctions that can better correlate with high-grade cytology, reflecting the presence of precancerous cells. In a proof-of-concept study conducted on 55 patients, starting from cervical smears, we have shown that HPV RNA-Seq can detect papillomaviruses with performances comparable to a widely used HPV reference molecular diagnostic kit; and a combination of the number of sequencing reads at specific early versus late HPV transcripts can be used as a marker of high-grade cytology, with encouraging diagnostic performances as a triage test.

Tracing cellular heterogeneity in pooled genetic screens via multi-level barcoding.

BACKGROUND: While pooled loss- and gain-of-function CRISPR screening approaches have become increasingly popular to systematically investigate mammalian gene function, the large majority of them have thus far not investigated the influence of cellular heterogeneity on screen results. Instead most screens are analyzed by averaging the abundance of perturbed cells from a bulk population of cells. RESULTS: Here we developed multi-level barcoded sgRNA libraries to trace multiple clonal Cas9 cell lines exposed to the same environment. The first level of barcoding allows monitoring growth kinetics and treatment responses of multiplexed clonal cell lines under identical conditions while the second level enables in-sample replication and tracing of sub-clonal lineages of cells expressing the same sgRNA. CONCLUSION: Using our approach, we illustrate how heterogeneity in growth kinetics and treatment response of clonal cell lines impairs the results of pooled genetic screens.

Tri-mannose grafting of chitosan nanocarriers remodels the macrophage response to bacterial infection.

BACKGROUND: Infectious diseases are still a leading cause of death and, with the emergence of drug resistance, pose a great threat to human health. New drugs and strategies are thus urgently needed to improve treatment efficacy and limit drug-associated side effects. Nanotechnology-based drug delivery systems are promising approaches, offering hope in the fight against drug resistant bacteria. However, how nanocarriers influence the response of innate immune cells to bacterial infection is mostly unknown. RESULTS: Here, we used Mycobacterium tuberculosis as a model of bacterial infection to examine the impact of mannose functionalization of chitosan nanocarriers (CS-NCs) on the human macrophage response. Both ungrafted and grafted CS-NCs were similarly internalized by macrophages, via an actin cytoskeleton-dependent process. Although tri-mannose ligands did not modify the capacity of CS-NCs to escape lysosomal degradation, they profoundly remodeled the response of M. tuberculosis-infected macrophages. mRNA sequencing showed nearly 900 genes to be differentially expressed due to tri-mannose grafting. Unexpectedly, the set of modulated genes was enriched for pathways involved in cell metabolism, particularly oxidative phosphorylation and sugar metabolism. CONCLUSIONS: The ability to modulate cell metabolism by grafting ligands at the surface of nanoparticles may thus be a promising strategy to reprogram immune cells and improve the efficacy of encapsulated drugs.

Altered spinogenesis in iPSC-derived cortical neurons from patients with autism carrying de novo SHANK3 mutations.

The synaptic protein SHANK3 encodes a multidomain scaffold protein expressed at the postsynaptic density of neuronal excitatory synapses. We previously identified de novo SHANK3 mutations in patients with autism spectrum disorders (ASD) and showed that SHANK3 represents one of the major genes for ASD. Here, we analyzed the pyramidal cortical neurons derived from induced pluripotent stem cells from four patients with ASD carrying SHANK3 de novo truncating mutations. At 40-45 days after the differentiation of neural stem cells, dendritic spines from pyramidal neurons presented variable morphologies: filopodia, thin, stubby and muschroom, as measured in 3D using GFP labeling and immunofluorescence. As compared to three controls, we observed a significant decrease in SHANK3 mRNA levels (less than 50% of controls) in correlation with a significant reduction in dendritic spine densities and whole spine and spine head volumes. These results, obtained through the analysis of de novo SHANK3 mutations in the patients' genomic background, provide further support for the presence of synaptic abnormalities in a subset of patients with ASD.

Dual gene activation and knockout screen reveals directional dependencies in genetic networks.

Understanding the direction of information flow is essential for characterizing how genetic networks affect phenotypes. However, methods to find genetic interactions largely fail to reveal directional dependencies. We combine two orthogonal Cas9 proteins from Streptococcus pyogenes and Staphylococcus aureus to carry out a dual screen in which one gene is activated while a second gene is deleted in the same cell. We analyze the quantitative effects of activation and knockout to calculate genetic interaction and directionality scores for each gene pair. Based on the results from over 100,000 perturbed gene pairs, we reconstruct a directional dependency network for human K562 leukemia cells and demonstrate how our approach allows the determination of directionality in activating genetic interactions. Our interaction network connects previously uncharacterized genes to well-studied pathways and identifies targets relevant for therapeutic intervention.