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.

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.

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.

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.

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.

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.

A positive feedback between p53 and miR-34 miRNAs mediates tumor suppression.

As bona fide p53 transcriptional targets, miR-34 microRNAs (miRNAs) exhibit frequent alterations in many human tumor types and elicit multiple p53 downstream effects upon overexpression. Unexpectedly, miR-34 deletion alone fails to impair multiple p53-mediated tumor suppressor effects in mice, possibly due to the considerable redundancy in the p53 pathway. Here, we demonstrate that miR-34a represses HDM4, a potent negative regulator of p53, creating a positive feedback loop acting on p53. In a Kras-induced mouse lung cancer model, miR-34a deficiency alone does not exhibit a strong oncogenic effect. However, miR-34a deficiency strongly promotes tumorigenesis when p53 is haploinsufficient, suggesting that the defective p53-miR-34 feedback loop can enhance oncogenesis in a specific context. The importance of the p53/miR-34/HDM4 feedback loop is further confirmed by an inverse correlation between miR-34 and full-length HDM4 in human lung adenocarcinomas. In addition, human lung adenocarcinomas generate an elevated level of a short HDM4 isoform through alternative polyadenylation. This short HDM4 isoform lacks miR-34-binding sites in the 3' untranslated region (UTR), thereby evading miR-34 regulation to disable the p53-miR-34 positive feedback. Taken together, our results elucidated the intricate cross-talk between p53 and miR-34 miRNAs and revealed an important tumor suppressor effect generated by this positive feedback loop.

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

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.

CRISPR-Cas9-mediated functional dissection of 3'-UTRs.

Many studies using reporter assays have demonstrated that 3' untranslated regions (3'-UTRs) regulate gene expression by controlling mRNA stability and translation. Due to intrinsic limitations of heterologous reporter assays, we sought to develop a gene editing approach to investigate the regulatory activity of 3'-UTRs in their native context. We initially used dual-CRISPR (clustered, regularly interspaced, short palindromic repeats)-Cas9 targeting to delete DNA regions corresponding to nine chemokine 3'-UTRs that destabilized mRNA in a reporter assay. Targeting six chemokine 3'-UTRs increased chemokine mRNA levels as expected. However, targeting CXCL1, CXCL6 and CXCL8 3'-UTRs unexpectedly led to substantial mRNA decreases. Metabolic labeling assays showed that targeting these three 3'-UTRs increased mRNA stability, as predicted by the reporter assay, while also markedly decreasing transcription, demonstrating an unexpected role for 3'-UTR sequences in transcriptional regulation. We further show that CRISPR-Cas9 targeting of specific 3'-UTR elements can be used for modulating gene expression and for highly parallel localization of active 3'-UTR elements in the native context. Our work demonstrates the duality and complexity of 3'-UTR sequences in regulation of gene expression and provides a useful approach for modulating gene expression and for functional annotation of 3'-UTRs in the native context.

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.

The landscape of genomic imprinting across diverse adult human tissues.

Genomic imprinting is an important regulatory mechanism that silences one of the parental copies of a gene. To systematically characterize this phenomenon, we analyze tissue specificity of imprinting from allelic expression data in 1582 primary tissue samples from 178 individuals from the Genotype-Tissue Expression (GTEx) project. We characterize imprinting in 42 genes, including both novel and previously identified genes. Tissue specificity of imprinting is widespread, and gender-specific effects are revealed in a small number of genes in muscle with stronger imprinting in males. IGF2 shows maternal expression in the brain instead of the canonical paternal expression elsewhere. Imprinting appears to have only a subtle impact on tissue-specific expression levels, with genes lacking a systematic expression difference between tissues with imprinted and biallelic expression. In summary, our systematic characterization of imprinting in adult tissues highlights variation in imprinting between genes, individuals, and tissues.

Laser capture microdissection enables transcriptomic analysis of dividing and quiescent liver stages of Plasmodium relapsing species.

Dormant liver stage forms (hypnozoites) of the malaria parasite Plasmodium vivax present major hurdles to control and eradicate infection. Despite major research efforts, the molecular composition of hypnozoites remains ill defined. Here, we applied a combination of state-of-the-art technologies to generate the first transcriptome of hypnozoites. We developed a robust laser dissection microscopy protocol to isolate individual Plasmodium cynomolgi hypnozoites and schizonts from infected monkey hepatocytes and optimized RNA-seq analysis to obtain the first transcriptomes of these stages. Comparative transcriptomic analysis identified 120 transcripts as being differentially expressed in the hypnozoite stage relative to the dividing liver schizont, with 69 and 51 mRNAs being up- or down-regulated, respectively, in the hypnozoites. This lead to the identification of potential markers of commitment to and maintenance of the dormant state of the hypnozoite including three transcriptional regulators of the ApiAP2 family, one of which is unique to P. cynomolgi and P. vivax, and the global translational repressor, eIF2a kinase eIK2, all of which are upregulated in the hypnozoite. Together, this work not only provides a primary experimentally-derived list of molecular markers of hypnozoites but also identifies transcriptional and posttranscriptional regulation of gene expression as potentially being key to establishing and maintaining quiescence.

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.

Gene List significance at-a-glance with GeneValorization.

MOTIVATION: High-throughput technologies provide fundamental informations concerning thousands of genes. Many of the current research laboratories daily use one or more of these technologies and end-up with lists of genes. Assessing the originality of the results obtained includes being aware of the number of publications available concerning individual or multiple genes and accessing information about these publications. Faced with the exponential growth of publications avaliable and number of genes involved in a study, this task is becoming particularly difficult to achieve. RESULTS: We introduce GeneValorization, a web-based tool that gives a clear and handful overview of the bibliography available corresponding to the user input formed by (i) a gene list (expressed by gene names or ids from EntrezGene) and (ii) a context of study (expressed by keywords). From this input, GeneValorization provides a matrix containing the number of publications with co-occurrences of gene names and keywords. Graphics are automatically generated to assess the relative importance of genes within various contexts. Links to publications and other databases offering information on genes and keywords are also available. To illustrate how helpful GeneValorization is, we will consider the gene list of the OncotypeDX prognostic marker test. AVAILABILITY: http://bioguide-project.net/gv CONTACT: cohen@lri.fr SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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