Single cell transcriptomics of cerebrospinal fluid cells from patients with recent-onset narcolepsy.

Narcolepsy is a rare cause of hypersomnolence and may be associated or not with cataplexy, i.e. sudden muscle weakness. These forms are designated narcolepsy-type 1 (NT1) and -type 2 (NT2), respectively. Notable characteristics of narcolepsy are that most patients carry the HLA-DQB1*06:02 allele and NT1-patients have strongly decreased levels of hypocretin-1 (synonym orexin-A) in the cerebrospinal fluid (CSF). The pathogenesis of narcolepsy is still not completely understood but the strong HLA-bias and increased frequencies of CD4+ T cells reactive to hypocretin in the peripheral blood suggest autoimmune processes in the hypothalamus. Here we analyzed the transcriptomes of CSF-cells from twelve NT1 and two NT2 patients by single cell RNAseq (scRNAseq). As controls, we used CSF cells from patients with multiple sclerosis, radiologically isolated syndrome, and idiopathic intracranial hypertension. From 27,255 CSF cells, we identified 20 clusters of different cell types and found significant differences in three CD4+ T cell and one monocyte clusters between narcolepsy and multiple sclerosis patients. Over 1000 genes were differentially regulated between patients with NT1 and other diseases. Surprisingly, the most strongly upregulated genes in narcolepsy patients as compared to controls were coding for the genome-encoded MTRNR2L12 and MTRNR2L8 peptides, which are homologous to the mitochondria-encoded HUMANIN peptide that is known playing a role in other neurological diseases including Alzheimer's disease.

Localized osmotic stress activates systemic responses to N limitation in Medicago truncatula-Sinorhizobium symbiotic plants.

In mature symbiotic root nodules, differentiated rhizobia fix atmospheric dinitrogen and provide ammonium to fulfill the plant nitrogen (N) demand. The plant enables this process by providing photosynthates to the nodules. The symbiosis is adjusted to the whole plant N demand thanks to systemic N signaling controlling nodule development. Symbiotic plants under N deficit stimulate nodule expansion and activate nodule senescence under N satiety. Besides, nodules are highly sensitive to drought. Here, we used split-root systems to characterize the systemic responses of symbiotic plants to a localized osmotic stress. We showed that polyéthylène glycol (PEG) application rapidly inhibited the symbiotic dinitrogen fixation activity of nodules locally exposed to the treatment, resulting to the N limitation of the plant supplied exclusively by symbiotic dinitrogen fixation. The localized PEG treatment triggered systemic signaling stimulating nodule development in the distant untreated roots. This response was associated with an enhancement of the sucrose allocation. Our analyses showed that transcriptomic reprogramming associated with PEG and N deficit systemic signaling(s) shared many targets transcripts. Altogether, our study suggests that systemic N signaling is a component of the adaptation of the symbiotic plant to the local variations of its edaphic environment.

The regenerative response of cardiac interstitial cells.

Understanding how certain animals are capable of regenerating their hearts will provide much needed insights into how this process can be induced in humans in order to reverse the damage caused by myocardial infarction. Currently, it is becoming increasingly evident that cardiac interstitial cells play crucial roles during cardiac regeneration. To understand how interstitial cells behave during this process, we performed single-cell RNA sequencing of regenerating zebrafish hearts. Using a combination of immunohistochemistry, chemical inhibition, and novel transgenic animals, we were able to investigate the role of cell type-specific mechanisms during cardiac regeneration. This approach allowed us to identify a number of important regenerative processes within the interstitial cell populations. Here, we provide detailed insight into how interstitial cells behave during cardiac regeneration, which will serve to increase our understanding of how this process could eventually be induced in humans.

The landscape of cancer-associated fibroblasts in colorectal cancer liver metastases.

Rationale: Patients with colorectal cancer die mainly due to liver metastases (CRC-LM). Although the tumor microenvironment (TME) plays an important role in tumor development and therapeutic response, our understanding of the individual TME components, especially cancer-associated fibroblasts (CAFs), remains limited. Methods: We analyzed CRC-LM CAFs and cancer cells by single-cell transcriptomics and used bioinformatics for data analysis and integration with related available single-cell and bulk transcriptomic datasets. We validated key findings by RT-qPCR, western blotting, and immunofluorescence. Results: By single-cell transcriptomic analysis of 4,397 CAFs from six CRC-LM samples, we identified two main CAF populations, contractile CAFs and extracellular matrix (ECM)-remodeling/pro-angiogenic CAFs, and four subpopulations with distinct phenotypes. We found that ECM-remodeling/pro-angiogenic CAFs derive from portal resident fibroblasts. They associate with areas of strong desmoplastic reaction and Wnt signaling in low-proliferating tumor cells engulfed in a stiff extracellular matrix. By integrating public single-cell primary liver tumor data, we propose a model to explain how different liver malignancies recruit CAFs of different origins to this organ. Lastly, we found that LTBP2 plays an important role in modulating collagen biosynthesis, ECM organization, and adhesion pathways. We developed fully human antibodies against LTBP2 that depleted LTBP2+ CAFs in vitro. Conclusion: This study complements recent reports on CRC-LM CAF heterogeneity at the single-cell resolution. The number of sequenced CAFs was more than one order of magnitude larger compared to existing data. LTBP2 targeting by antibodies might create opportunities to deplete ECM-remodeling CAFs in CRC-LMs. This might be combined with other therapies, e.g., anti-angiogenic compounds as already done in CRC. Moreover, we showed that in intrahepatic cholangiocarcinoma, in which ECM-remodeling CAF proportion is similar to that of CRC-LM, several genes expressed by ECM-remodeling CAFs, such as LTBP2, were associated with survival.

Major transcriptomic, epigenetic and metabolic changes underlie the pluripotency continuum in rabbit preimplantation embryos.

Despite the growing interest in the rabbit model for developmental and stem cell biology, the characterization of embryos at the molecular level is still poorly documented. We conducted a transcriptome analysis of rabbit preimplantation embryos from E2.7 (morula stage) to E6.6 (early primitive streak stage) using bulk and single-cell RNA-sequencing. In parallel, we studied oxidative phosphorylation and glycolysis, and analysed active and repressive epigenetic modifications during blastocyst formation and expansion. We generated a transcriptomic, epigenetic and metabolic map of the pluripotency continuum in rabbit preimplantation embryos, and identified novel markers of naive pluripotency that might be instrumental for deriving naive pluripotent stem cell lines. Although the rabbit is evolutionarily closer to mice than to primates, we found that the transcriptome of rabbit epiblast cells shares common features with those of humans and non-human primates.

H3K9me2 genome-wide distribution in the holocentric insect Spodoptera frugiperda (Lepidoptera: Noctuidae).

BACKGROUND: Eukaryotic genomes are packaged by Histone proteins in a structure called chromatin. There are different chromatin types. Euchromatin is typically associated with decondensed, transcriptionally active regions and heterochromatin to more condensed regions of the chromosomes. Methylation of Lysine 9 of Histone H3 (H3K9me) is a conserved biochemical marker of heterochromatin. In many organisms, heterochromatin is usually localized at telomeric as well as pericentromeric regions but can also be found at interstitial chromosomal loci. This distribution may vary in different species depending on their general chromosomal organization. Holocentric species such as Spodoptera frugiperda (Lepidoptera: Noctuidae) possess dispersed centromeres instead of a monocentric one and thus no observable pericentromeric compartment. To identify the localization of heterochromatin in such species we performed ChIP-Seq experiments and analyzed the distribution of the heterochromatin marker H3K9me2 in the Sf9 cell line and whole 4th instar larvae (L4) in relation to RNA-Seq data. RESULTS: In both samples we measured an enrichment of H3K9me2 at the (sub) telomeres, rDNA loci, and satellite DNA sequences, which could represent dispersed centromeric regions. We also observed that density of H3K9me2 is positively correlated with transposable elements and protein-coding genes. But contrary to most model organisms, H3K9me2 density is not correlated with transcriptional repression. CONCLUSION: This is the first genome-wide ChIP-Seq analysis conducted in S. frugiperda for H3K9me2. Compared to model organisms, this mark is found in expected chromosomal compartments such as rDNA and telomeres. However, it is also localized at numerous dispersed regions, instead of the well described large pericentromeric domains, indicating that H3K9me2 might not represent a classical heterochromatin marker in Lepidoptera. (242 words).

The Type III Effectome of the Symbiotic Bradyrhizobium vignae Strain ORS3257.

Many Bradyrhizobium strains are able to establish a Nod factor-independent symbiosis with the leguminous plant Aeschynomene indica by the use of a type III secretion system (T3SS). Recently, an important advance in the understanding of the molecular factors supporting this symbiosis has been achieved by the in silico identification and functional characterization of 27 putative T3SS effectors (T3Es) of Bradyrhizobium vignae ORS3257. In the present study, we experimentally extend this catalog of T3Es by using a multi-omics approach. Transcriptome analysis under non-inducing and inducing conditions in the ORS3257 wild-type strain and the ttsI mutant revealed that the expression of 18 out of the 27 putative effectors previously identified, is under the control of TtsI, the global transcriptional regulator of T3SS and T3Es. Quantitative shotgun proteome analysis of culture supernatant in the wild type and T3SS mutant strains confirmed that 15 of the previously determined candidate T3Es are secreted by the T3SS. Moreover, the combined approaches identified nine additional putative T3Es and one of them was experimentally validated as a novel effector. Our study underscores the power of combined proteome and transcriptome analyses to complement in silico predictions and produce nearly complete effector catalogs. The establishment of the ORS3257 effectome will form the basis for a full appraisal of the symbiotic properties of this strain during its interaction with various host legumes via different processes.

NRG1/ErbB signalling controls the dialogue between macrophages and neural crest-derived cells during zebrafish fin regeneration.

Fish species, such as zebrafish (Danio rerio), can regenerate their appendages after amputation through the formation of a heterogeneous cellular structure named blastema. Here, by combining live imaging of triple transgenic zebrafish embryos and single-cell RNA sequencing we established a detailed cell atlas of the regenerating caudal fin in zebrafish larvae. We confirmed the presence of macrophage subsets that govern zebrafish fin regeneration, and identified a foxd3-positive cell population within the regenerating fin. Genetic depletion of these foxd3-positive neural crest-derived cells (NCdC) showed that they are involved in blastema formation and caudal fin regeneration. Finally, chemical inhibition and transcriptomic analysis demonstrated that these foxd3-positive cells regulate macrophage recruitment and polarization through the NRG1/ErbB pathway. Here, we show the diversity of the cells required for blastema formation, identify a discrete foxd3-positive NCdC population, and reveal the critical function of the NRG1/ErbB pathway in controlling the dialogue between macrophages and NCdC.

Systemic control of nodule formation by plant nitrogen demand requires autoregulation-dependent and independent mechanisms.

In legumes interacting with rhizobia, the formation of symbiotic organs involved in the acquisition of atmospheric nitrogen gas (N2) is dependent on the plant nitrogen (N) demand. We used Medicago truncatula plants cultivated in split-root systems to discriminate between responses to local and systemic N signaling. We evidenced a strong control of nodule formation by systemic N signaling but obtained no clear evidence of a local control by mineral nitrogen. Systemic signaling of the plant N demand controls numerous transcripts involved in root transcriptome reprogramming associated with early rhizobia interaction and nodule formation. SUPER NUMERIC NODULES (SUNN) has an important role in this control, but we found that major systemic N signaling responses remained active in the sunn mutant. Genes involved in the activation of nitrogen fixation are regulated by systemic N signaling in the mutant, explaining why its hypernodulation phenotype is not associated with higher nitrogen fixation of the whole plant. We show that the control of transcriptome reprogramming of nodule formation by systemic N signaling requires other pathway(s) that parallel the SUNN/CLE (CLAVATA3/EMBRYO SURROUNDING REGION-LIKE PEPTIDES) pathway.

Sublethal Exposure Effects of the Neonicotinoid Clothianidin Strongly Modify the Brain Transcriptome and Proteome in the Male Moth Agrotis ipsilon.

Insect pest management relies mainly on neurotoxic insecticides, including neonicotinoids such as clothianidin. The residual accumulation of low concentrations of these insecticides can have positive effects on target pest insects by enhancing various life traits. Because pest insects often rely on sex pheromones for reproduction and olfactory synaptic transmission is cholinergic, neonicotinoid residues could indeed modify chemical communication. We recently showed that treatments with low doses of clothianidin could induce hormetic effects on behavioral and neuronal sex pheromone responses in the male moth, Agrotis ipsilon. In this study, we used high-throughput RNAseq and proteomic analyses from brains of A. ipsilon males that were intoxicated with a low dose of clothianidin to investigate the molecular mechanisms leading to the observed hormetic effect. Our results showed that clothianidin induced significant changes in transcript levels and protein quantity in the brain of treated moths: 1229 genes and 49 proteins were differentially expressed upon clothianidin exposure. In particular, our analyses highlighted a regulation in numerous enzymes as a possible detoxification response to the insecticide and also numerous changes in neuronal processes, which could act as a form of acclimatization to the insecticide-contaminated environment, both leading to enhanced neuronal and behavioral responses to sex pheromone.

Differential immune gene expression associated with contemporary range expansion in two invasive rodents in Senegal.

Biological invasions are major anthropogenic changes associated with threats to biodiversity and health. However, what determines the successful establishment and spread of introduced populations remains unclear. Here, we explore several hypotheses linking invasion success and immune phenotype traits, including those based on the evolution of increased competitive ability concept. We compared gene expression profiles between anciently and recently established populations of two major invading species, the house mouse Mus musculus domesticus and the black rat Rattus rattus, in Senegal (West Africa). Transcriptome analyses identified differential expression between anciently and recently established populations for 364 mouse genes and 83 rat genes. All immune-related genes displaying differential expression along the mouse invasion route were overexpressed at three of the four recently invaded sites studied. Complement activation pathway genes were overrepresented among these genes. By contrast, no particular immunological process was found to be overrepresented among the differentially expressed genes of black rat. Changes in transcriptome profiles were thus observed along invasion routes, but with different specific patterns between the two invasive species. These changes may be driven by increases in infection risks at sites recently invaded by the house mouse, and by stochastic events associated with colonization history for the black rat. These results constitute a first step toward the identification of immune eco-evolutionary processes potentially involved in the invasion success of these two rodent species.

TALC: Transcript-level Aware Long-read Correction.

MOTIVATION: Long-read sequencing technologies are invaluable for determining complex RNA transcript architectures but are error-prone. Numerous 'hybrid correction' algorithms have been developed for genomic data that correct long reads by exploiting the accuracy and depth of short reads sequenced from the same sample. These algorithms are not suited for correcting more complex transcriptome sequencing data. RESULTS: We have created a novel reference-free algorithm called Transcript-level Aware Long-Read Correction (TALC) which models changes in RNA expression and isoform representation in a weighted De Bruijn graph to correct long reads from transcriptome studies. We show that transcript-level aware correction by TALC improves the accuracy of the whole spectrum of downstream RNA-seq applications and is thus necessary for transcriptome analyses that use long read technology. AVAILABILITY AND IMPLEMENTATION: TALC is implemented in C++ and available at https://github.com/lbroseus/TALC. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The role of copy-number variation in the reinforcement of sexual isolation between the two European subspecies of the house mouse.

Reinforcement has the potential to generate strong reproductive isolation through the evolution of barrier traits as a response to selection against maladaptive hybridization, but the genetic changes associated with this process remain largely unexplored. Building upon the increasing evidence for a role of structural variants in adaptation and speciation, we addressed the role of copy-number variation in the reinforcement of sexual isolation evidenced between the two European subspecies of the house mouse. We characterized copy-number divergence between populations of Mus musculus musculus that display assortative mate choice, and those that do not, using whole-genome resequencing data. Updating methods to detect deletions and tandem duplications (collectively: copy-number variants, CNVs) in Pool-Seq data, we developed an analytical pipeline dedicated to identifying genomic regions showing the expected pattern of copy-number displacement under a reinforcement scenario. This strategy allowed us to detect 1824 deletions and seven tandem duplications that showed extreme differences in frequency between behavioural classes across replicate comparisons. A subset of 480 deletions and four tandem duplications were specifically associated with the derived trait of assortative mate choice. These 'Choosiness-associated' CNVs occur in hundreds of genes. Consistent with our hypothesis, such genes included olfactory receptors potentially involved in the olfactory-based assortative mate choice in this system as well as one gene, Sp110, that is known to show patterns of differential expression between behavioural classes in an organ used in mate choice-the vomeronasal organ. These results demonstrate that fine-scale structural changes are common and highly variable within species, despite being under-studied, and may be important targets of reinforcing selection in this system and others. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.

A Transposon Story: From TE Content to TE Dynamic Invasion of Drosophila Genomes Using the Single-Molecule Sequencing Technology from Oxford Nanopore.

Transposable elements (TEs) are the main components of genomes. However, due to their repetitive nature, they are very difficult to study using data obtained with short-read sequencing technologies. Here, we describe an efficient pipeline to accurately recover TE insertion (TEI) sites and sequences from long reads obtained by Oxford Nanopore Technology (ONT) sequencing. With this pipeline, we could precisely describe the landscapes of the most recent TEIs in wild-type strains of Drosophila melanogaster and Drosophila simulans. Their comparison suggests that this subset of TE sequences is more similar than previously thought in these two species. The chromosome assemblies obtained using this pipeline also allowed recovering piRNA cluster sequences, which was impossible using short-read sequencing. Finally, we used our pipeline to analyze ONT sequencing data from a D. melanogaster unstable line in which LTR transposition was derepressed for 73 successive generations. We could rely on single reads to identify new insertions with intact target site duplications. Moreover, the detailed analysis of TEIs in the wild-type strains and the unstable line did not support the trap model claiming that piRNA clusters are hotspots of TE insertions.

Responses of mature symbiotic nodules to the whole-plant systemic nitrogen signaling.

In symbiotic root nodules of legumes, terminally differentiated rhizobia fix atmospheric N2 producing an NH4+ influx that is assimilated by the plant. The plant, in return, provides photosynthates that fuel the symbiotic nitrogen acquisition. Mechanisms responsible for the adjustment of the symbiotic capacity to the plant N demand remain poorly understood. We have investigated the role of systemic signaling of whole-plant N demand on the mature N2-fixing nodules of the model symbiotic association Medicago truncatula/Sinorhizobium using split-root systems. The whole-plant N-satiety signaling rapidly triggers reductions of both N2 fixation and allocation of sugars to the nodule. These responses are associated with the induction of nodule senescence and the activation of plant defenses against microbes, as well as variations in sugars transport and nodule metabolism. The whole-plant N-deficit responses mirror these changes: a rapid increase of sucrose allocation in response to N-deficit is associated with a stimulation of nodule functioning and development resulting in nodule expansion in the long term. Physiological, transcriptomic, and metabolomic data together provide evidence for strong integration of symbiotic nodules into whole-plant nitrogen demand by systemic signaling and suggest roles for sugar allocation and hormones in the signaling mechanisms.

Transcriptome of tambaqui Colossoma macropomum during gonad differentiation: Different molecular signals leading to sex identity.

Tambaqui (Colossoma macropomum) is the major native species in Brazilian aquaculture, and we have shown that females exhibit a higher growth compared to males, opening up the possibility for the production of all-female population. To date, there is no information on the sex determination and differentiation molecular mechanisms of tambaqui. In the present study, transcriptome sequencing of juvenile trunks was performed to understand the molecular network involved in the gonadal sex differentiation. The results showed that before differentiation, components of the Wnt/ß-catenin pathway, fox and fst genes imprint female sex development, whereas antagonistic pathways (gsk3b, wt1 and fgfr2), sox9 and genes for androgen synthesis indicate male differentiation. Hence, in undifferentiated tambaqui, the Wnt/ß-catenin exerts a role on sex differentiation, either upregulated in female-like individuals, or antagonized in male-like individuals.

TAG-RNAi overcomes off-target effects in cancer models.

RNA interference offers therapeutic opportunities for the clinical targeting of otherwise undruggable oncogenes. However RNAi can have off-target effects that considerably increase treatment risks. To manage these side effects and allow an easy subtraction of their activity in healthy tissues, we present here the TAG-RNAi approach where cells that are not designated targets do not have the mRNA tag. Using TAG-RNAi we first established the off-target signatures of three different siRNAs specific to the Cyclin D1 oncogene by RNA-sequencing of cultured cancer cells expressing a FLAG-HA-tagged-Cyclin D1. Then, by symmetrical allografts of tagged-cancer cells and untagged controls on the left and right flanks of model mice, we demonstrate that TAG-RNAi is a reliable approach to study the functional impact of any oncogene without off-target bias. Finally we show, as examples, that mutation-specific TAG-RNAi can be applied to downregulate two oncogenic mutants, KRAS-G12V or BRAF-V600E, while sparing the expression of the wild-type proteins. TAG-RNAi will thus avoid the traditional off-target limitations of RNAi in future experimental approaches.

Seed comparative genomics in three coffee species identify desiccation tolerance mechanisms in intermediate seeds.

In contrast to desiccation-tolerant 'orthodox' seeds, so-called 'intermediate' seeds cannot survive complete drying and are short-lived. All species of the genus Coffea produce intermediate seeds, but they show a considerable variability in seed desiccation tolerance (DT), which may help to decipher the molecular basis of seed DT in plants. We performed a comparative transcriptome analysis of developing seeds in three coffee species with contrasting desiccation tolerance. Seeds of all species shared a major transcriptional switch during late maturation that governs a general slow-down of metabolism. However, numerous key stress-related genes, including those coding for the late embryogenesis abundant protein EM6 and the osmosensitive calcium channel ERD4, were up-regulated during DT acquisition in the two species with high seed DT, C. arabica and C. eugenioides. By contrast, we detected up-regulation of numerous genes involved in the metabolism, transport, and perception of auxin in C. canephora seeds with low DT. Moreover, species with high DT showed a stronger down-regulation of the mitochondrial machinery dedicated to the tricarboxylic acid cycle and oxidative phosphorylation. Accordingly, respiration measurements during seed dehydration demonstrated that intermediate seeds with the highest DT are better prepared to cease respiration and avoid oxidative stresses.

Spodoptera frugiperda transcriptional response to infestation by Steinernema carpocapsae.

Steinernema carpocapsae is an entomopathogenic nematode (EPN) used in biological control of agricultural pest insects. It enters the hemocoel of its host via the intestinal tract and releases its symbiotic bacterium Xenorhabdus nematophila. In order to improve our knowledge about the physiological responses of its different hosts, we examined the transcriptional responses to EPN infestation of the fat body, the hemocytes and the midgut in the lepidopteran pest Spodoptera frugiperda. The tissues poorly respond to the infestation at an early time post-infestation of 8 h with only 5 genes differentially expressed in the fat body of the caterpillars. Strong transcriptional responses are observed at a later time point of 15 h post-infestation in all three tissues. Few genes are differentially expressed in the midgut but tissue-specific panels of induced metalloprotease inhibitors, immune receptors and antimicrobial peptides together with several uncharacterized genes are up-regulated in the fat body and the hemocytes. Among the most up-regulated genes, we identified new potential immune effectors, unique to Lepidoptera, which show homology with bacterial genes of unknown function. Altogether, these results pave the way for further functional studies of the responsive genes' involvement in the interaction with the EPN.

Plant population dynamics on oceanic islands during the Late Quaternary climate changes: genetic evidence from a tree species (Coffea mauritiana) in Reunion Island.

Past climatic fluctuations have played a major role in shaping the current plant biodiversity. Although harbouring an exceptional biota, oceanic islands have received little attention in studies on species demographic history and past vegetation patterns. We investigated the impact of past climatic changes on the effective population size of a tree (Coffea mauritiana) that is endemic to Reunion Island, located in the south-western Indian Ocean (SWIO). Demographic changes were inferred using summary statistics calculated from genomic data. Using ecological niche modelling and the current distribution of genetic diversity, the paleodistribution of the species was also assessed. A reduction in the effective population size of C. mauritiana during the last glaciation maximum was inferred. The distribution of the species was reduced on the western side of the island, due to low rainfall. It appeared that a major reduction in rainfall and a slight temperature decrease prevailed in the SWIO. Our findings indicated that analyses on the current patterns of intraspecific genetic variations can efficiently contribute to past climatic changes characterisation in remote islands. Identifying area with higher resilience in oceanic islands could provide guidance in forest management and conservation faced to the global climate change.