Computational analysis of super-resolved in situ sequencing data reveals genes modified by immune-tumor contact events.

Cancer cells can manipulate immune cells and escape from the immune system response. Quantifying the molecular changes that occur when an immune cell touches a tumor cell can increase our understanding of the underlying mechanisms. Recently, it became possible to perform such measurements in situ-for example, using expansion sequencing, which enabled in situ sequencing of genes with super-resolution. We systematically examined whether individual immune cells from specific cell types express genes differently when in physical proximity to individual tumor cells. First, we demonstrated that a dense mapping of genes in situ can be used for the segmentation of cell bodies in 3D, thus improving our ability to detect likely touching cells. Next, we used three different computational approaches to detect the molecular changes that are triggered by proximity: differential expression analysis, tree-based machine learning classifiers, and matrix factorization analysis. This systematic analysis revealed tens of genes, in specific cell types, whose expression separates immune cells that are proximal to tumor cells from those that are not proximal, with a significant overlap between the different detection methods. Remarkably, an order of magnitude more genes are triggered by proximity to tumor cells in CD8 T cells compared to CD4 T cells, in line with the ability of CD8 T cells to directly bind major histocompatibility complex (MHC) class I on tumor cells. Thus, in situ sequencing of an individual biopsy can be used to detect genes likely involved in immune-tumor cell-cell interactions. The data used in this manuscript and the code of the InSituSeg, machine learning, cNMF, and Moran's I methods are publicly available at doi:10.5281/zenodo.7497981.

MolOptimizer: A Molecular Optimization Toolkit for Fragment-Based Drug Design.

MolOptimizer is a user-friendly computational toolkit designed to streamline the hit-to-lead optimization process in drug discovery. MolOptimizer extracts features and trains machine learning models using a user-provided, labeled, and small-molecule dataset to accurately predict the binding values of new small molecules that share similar scaffolds with the target in focus. Hosted on the Azure web-based server, MolOptimizer emerges as a vital resource, accelerating the discovery and development of novel drug candidates with improved binding properties.

Behavioral Neuroscience in the Era of Genomics: Tools and Lessons for Analyzing High-Dimensional Datasets.

Behavioral neuroscience underwent a technology-driven revolution with the emergence of machine-vision and machine-learning technologies. These technological advances facilitated the generation of high-resolution, high-throughput capture and analysis of complex behaviors. Therefore, behavioral neuroscience is becoming a data-rich field. While behavioral researchers use advanced computational tools to analyze the resulting datasets, the search for robust and standardized analysis tools is still ongoing. At the same time, the field of genomics exploded with a plethora of technologies which enabled the generation of massive datasets. This growth of genomics data drove the emergence of powerful computational approaches to analyze these data. Here, we discuss the composition of a large behavioral dataset, and the differences and similarities between behavioral and genomics data. We then give examples of genomics-related tools that might be of use for behavioral analysis and discuss concepts that might emerge when considering the two fields together.

Urbanization comprehensively impairs biological rhythms in coral holobionts.

Coral reefs are in global decline due to climate change and anthropogenic influences (Hughes et al., Conservation Biology, 27: 261-269, 2013). Near coastal cities or other densely populated areas, coral reefs face a range of additional challenges. While considerable progress has been made in understanding coral responses to acute individual stressors (Dominoni et al., Nature Ecology & Evolution, 4: 502-511, 2020), the impacts of chronic exposure to varying combinations of sensory pollutants are largely unknown. To investigate the impacts of urban proximity on corals, we conducted a year-long in-natura study-incorporating sampling at diel, monthly, and seasonal time points-in which we compared corals from an urban area to corals from a proximal non-urban area. Here we reveal that despite appearing relatively healthy, natural biorhythms and environmental sensory systems were extensively disturbed in corals from the urban environment. Transcriptomic data indicated poor symbiont performance, disturbance to gametogenic cycles, and loss or shifted seasonality of vital biological processes. Altered seasonality patterns were also observed in the microbiomes of the urban coral population, signifying the impact of urbanization on the holobiont, rather than the coral host alone. These results should raise alarm regarding the largely unknown long-term impacts of sensory pollution on the resilience and survival of coral reefs close to coastal communities.

Barcoded oligonucleotides ligated on RNA amplified for multiplexed and parallel in situ analyses.

We present barcoded oligonucleotides ligated on RNA amplified for multiplexed and parallel insitu analyses (BOLORAMIS), a reverse transcription-free method for spatially-resolved, targeted, in situ RNA identification of single or multiple targets. BOLORAMIS was demonstrated on a range of cell types and human cerebral organoids. Singleplex experiments to detect coding and non-coding RNAs in human iPSCs showed a stem-cell signature pattern. Specificity of BOLORAMIS was found to be 92% as illustrated by a clear distinction between human and mouse housekeeping genes in a co-culture system, as well as by recapitulation of subcellular localization of lncRNA MALAT1. Sensitivity of BOLORAMIS was quantified by comparing with single molecule FISH experiments and found to be 11%, 12% and 35% for GAPDH, TFRC and POLR2A, respectively. To demonstrate BOLORAMIS for multiplexed gene analysis, we targeted 96 mRNAs within a co-culture of iNGN neurons and HMC3 human microglial cells. We used fluorescence in situ sequencing to detect error-robust 8-base barcodes associated with each of these genes. We then used this data to uncover the spatial relationship among cells and transcripts by performing single-cell clustering and gene-gene proximity analyses. We anticipate the BOLORAMIS technology for in situ RNA detection to find applications in basic and translational research.

Expansion sequencing: Spatially precise in situ transcriptomics in intact biological systems.

Methods for highly multiplexed RNA imaging are limited in spatial resolution and thus in their ability to localize transcripts to nanoscale and subcellular compartments. We adapt expansion microscopy, which physically expands biological specimens, for long-read untargeted and targeted in situ RNA sequencing. We applied untargeted expansion sequencing (ExSeq) to the mouse brain, which yielded the readout of thousands of genes, including splice variants. Targeted ExSeq yielded nanoscale-resolution maps of RNAs throughout dendrites and spines in the neurons of the mouse hippocampus, revealing patterns across multiple cell types, layer-specific cell types across the mouse visual cortex, and the organization and position-dependent states of tumor and immune cells in a human metastatic breast cancer biopsy. Thus, ExSeq enables highly multiplexed mapping of RNAs from nanoscale to system scale.

Resource: A multi-species multi-timepoint transcriptome database and webpage for the pineal gland and retina.

The website and database https://snengs.nichd.nih.gov provides RNA sequencing data from multi-species analysis of the pineal glands from zebrafish (Danio rerio), chicken (White Leghorn), rat (Rattus novegicus), mouse (Mus musculus), rhesus macaque (Macaca mulatta), and human (Homo sapiens); in most cases, retinal data are also included along with results of the analysis of a mixture of RNA from tissues. Studies cover day and night conditions; in addition, a time series over multiple hours, a developmental time series and pharmacological experiments on rats are included. The data have been uniformly re-processed using the latest methods and assemblies to allow for comparisons between experiments and to reduce processing differences. The website presents search functionality, graphical representations, Excel tables, and track hubs of all data for detailed visualization in the UCSC Genome Browser. As more data are collected from investigators and improved genomes become available in the future, the website will be updated. This database is in the public domain and elements can be reproduced by citing the URL and this report. This effort makes the results of 21st century transcriptome profiling widely available in a user-friendly format that is expected to broadly influence pineal research.

Expansion microscopy: enabling single cell analysis in intact biological systems.

There is a need for single cell analysis methods that enable the identification and localization of different kinds of biomolecules throughout cells and intact tissues, thereby allowing characterization and classification of individual cells and their relationships to each other within intact systems. Expansion microscopy (ExM) is a technology that physically magnifies tissues in an isotropic way, thereby achieving super-resolution microscopy on diffraction-limited microscopes, enabling rapid image acquisition and large field of view. As a result, ExM is well-positioned to integrate molecular content and cellular morphology, with the spatial precision sufficient to resolve individual biological building blocks, and the scale and accessibility required to deploy over extended 3-D objects like tissues and organs.

ADAR2/miR-589-3p axis controls glioblastoma cell migration/invasion.

Recent studies have reported the emerging role of microRNAs (miRNAs) in human cancers. We systematically characterized miRNA expression and editing in the human brain, which displays the highest number of A-to-I RNA editing sites among human tissues, and in de novo glioblastoma brain cancer. We identified 299 miRNAs altered in their expression and 24 miRNAs differently edited in human brain compared to glioblastoma tissues. We focused on the editing site within the miR-589-3p seed. MiR-589-3p is a unique miRNA almost fully edited (∼100%) in normal brain and with a consistent editing decrease in glioblastoma. The edited version of miR-589-3p inhibits glioblastoma cell proliferation, migration and invasion, while the unedited version boosts cell proliferation and motility/invasion, thus being a potential cancer-promoting factor. We demonstrated that the editing of this miRNA is mediated by ADAR2, and retargets miR-589-3p from the tumor-suppressor PCDH9 to ADAM12, which codes for the metalloproteinase 12 promoting glioblastoma invasion. Overall, our study dissects the role of a unique brain-specific editing site within miR-589-3p, with important anticancer features, and highlights the importance of RNA editing as an essential player not only for diversifying the genomic message but also for correcting not-tolerable/critical genomic coding sites.

A-to-I RNA Editing in the Earliest-Diverging Eumetazoan Phyla.

The highly conserved ADAR enzymes, found in all multicellular metazoans, catalyze the editing of mRNA transcripts by the deamination of adenosines to inosines. This type of editing has two general outcomes: site specific editing, which frequently leads to recoding, and clustered editing, which is usually found in transcribed genomic repeats. Here, for the first time, we looked for both editing of isolated sites and clustered, non-specific sites in a basal metazoan, the coral Acropora millepora during spawning event, in order to reveal its editing pattern. We found that the coral editome resembles the mammalian one: it contains more than 500,000 sites, virtually all of which are clustered in non-coding regions that are enriched for predicted dsRNA structures. RNA editing levels were increased during spawning and increased further still in newly released gametes. This may suggest that editing plays a role in introducing variability in coral gametes.

Trade-off between Transcriptome Plasticity and Genome Evolution in Cephalopods.

RNA editing, a post-transcriptional process, allows the diversification of proteomes beyond the genomic blueprint; however it is infrequently used among animals for this purpose. Recent reports suggesting increased levels of RNA editing in squids thus raise the question of the nature and effects of these events. We here show that RNA editing is particularly common in behaviorally sophisticated coleoid cephalopods, with tens of thousands of evolutionarily conserved sites. Editing is enriched in the nervous system, affecting molecules pertinent for excitability and neuronal morphology. The genomic sequence flanking editing sites is highly conserved, suggesting that the process confers a selective advantage. Due to the large number of sites, the surrounding conservation greatly reduces the number of mutations and genomic polymorphisms in protein-coding regions. This trade-off between genome evolution and transcriptome plasticity highlights the importance of RNA recoding as a strategy for diversifying proteins, particularly those associated with neural function. PAPERCLIP.

Genetically Blocking the Zebrafish Pineal Clock Affects Circadian Behavior.

The master circadian clock in fish has been considered to reside in the pineal gland. This dogma is challenged, however, by the finding that most zebrafish tissues contain molecular clocks that are directly reset by light. To further examine the role of the pineal gland oscillator in the zebrafish circadian system, we generated a transgenic line in which the molecular clock is selectively blocked in the melatonin-producing cells of the pineal gland by a dominant-negative strategy. As a result, clock-controlled rhythms of melatonin production in the adult pineal gland were disrupted. Moreover, transcriptome analysis revealed that the circadian expression pattern of the majority of clock-controlled genes in the adult pineal gland is abolished. Importantly, circadian rhythms of behavior in zebrafish larvae were affected: rhythms of place preference under constant darkness were eliminated, and rhythms of locomotor activity under constant dark and constant dim light conditions were markedly attenuated. On the other hand, global peripheral molecular oscillators, as measured in whole larvae, were unaffected in this model. In conclusion, characterization of this novel transgenic model provides evidence that the molecular clock in the melatonin-producing cells of the pineal gland plays a key role, possibly as part of a multiple pacemaker system, in modulating circadian rhythms of behavior.

Nanoscale imaging of RNA with expansion microscopy.

The ability to image RNA identity and location with nanoscale precision in intact tissues is of great interest for defining cell types and states in normal and pathological biological settings. Here, we present a strategy for expansion microscopy of RNA. We developed a small-molecule linker that enables RNA to be covalently attached to a swellable polyelectrolyte gel synthesized throughout a biological specimen. Then, postexpansion, fluorescent in situ hybridization (FISH) imaging of RNA can be performed with high yield and specificity as well as single-molecule precision in both cultured cells and intact brain tissue. Expansion FISH (ExFISH) separates RNAs and supports amplification of single-molecule signals (i.e., via hybridization chain reaction) as well as multiplexed RNA FISH readout. ExFISH thus enables super-resolution imaging of RNA structure and location with diffraction-limited microscopes in thick specimens, such as intact brain tissue and other tissues of importance to biology and medicine.

Knockdown of unc119c results in visual impairment and early-onset retinal dystrophy in zebrafish.

PURPOSE: UNC119 proteins are involved in G protein trafficking in mouse retinal photoreceptors and Caenorhabditis elegans olfactory neurons. An Unc119 null allele is associated with cone-rod dystrophy in mouse, but the mechanism leading to disease is not understood. We studied the role of Unc119 paralogs and Arl3l2 in zebrafish vision and retinal organization resulting from unc119c and arl3l2 knockdown. METHODS: Zebrafish unc119c was amplified by PCR from retina and pineal gland cDNA. Its expression pattern in the eye and pineal gland was determined by whole-mount in-situ hybridization. unc119c and arl3l2 were knocked down using morpholino-modified oligonucleotides (MO). Their visual function was assessed with a quantitative optomotor assay on 6 days post-fertilization larvae. Retinal morphology was analyzed using immunohistochemistry with anti-cone arrestin (zpr-1) and anti-cone transducin-α (GNAT2) antibodies. RESULTS: The zebrafish genome contains four genes encoding unc119 paralogs located on different chromosomes. The exon/intron arrangements of these genes are identical. Three Unc119 paralogs are expressed in the zebrafish retina, termed Unc119a-c. Based on sequence similarity, Unc119a and Unc119b are orthologs of mammalian UNC119a and UNC119b, respectively. A third, Unc119c, is unique and not present in mammals. Whole mount in-situ hybridization revealed that unc119a and unc119b RNA are ubiquitously expressed in the CNS, and unc119c is specifically expressed in photoreceptive tissues (pineal gland and retina). A Unc119 interactant, Arl3l2 also localizes to the pineal gland and the retina. As measured by the optomotor response, unc119c and arl3l2 knockdown resulted in significantly lower vision compared to wild-type zebrafish larvae and control morpholino (MO). Immunohistological analysis with anti-cone transducin and anti-cone arrestin (zpr-1) indicates that knockdown of unc119c leads to photoreceptor degeneration mostly affecting cones. CONCLUSIONS: Our results suggest that Unc119c is the only Unc119 paralog that is highly specific to the retina in zebrafish. Unc119c and Arl3l2 proteins are important for the function of cones.

Signaling cascades and the importance of moonlight in coral broadcast mass spawning.

Many reef-building corals participate in a mass-spawning event that occurs yearly on the Great Barrier Reef. This coral reproductive event is one of earth's most prominent examples of synchronised behavior, and coral reproductive success is vital to the persistence of coral reef ecosystems. Although several environmental cues have been implicated in the timing of mass spawning, the specific sensory cues that function together with endogenous clock mechanisms to ensure accurate timing of gamete release are largely unknown. Here, we show that moonlight is an important external stimulus for mass spawning synchrony and describe the potential mechanisms underlying the ability of corals to detect environmental triggers for the signaling cascades that ultimately result in gamete release. Our study increases the understanding of reproductive chronobiology in corals and strongly supports the hypothesis that coral gamete release is achieved by a complex array of potential neurohormones and light-sensing molecules.

Hypocretin neuron-specific transcriptome profiling identifies the sleep modulator Kcnh4a.

Sleep has been conserved throughout evolution; however, the molecular and neuronal mechanisms of sleep are largely unknown. The hypothalamic hypocretin/orexin (Hcrt) neurons regulate sleep\wake states, feeding, stress, and reward. To elucidate the mechanism that enables these various functions and to identify sleep regulators, we combined fluorescence cell sorting and RNA-seq in hcrt:EGFP zebrafish. Dozens of Hcrt-neuron-specific transcripts were identified and comprehensive high-resolution imaging revealed gene-specific localization in all or subsets of Hcrt neurons. Clusters of Hcrt-neuron-specific genes are predicted to be regulated by shared transcription factors. These findings show that Hcrt neurons are heterogeneous and that integrative molecular mechanisms orchestrate their diverse functions. The voltage-gated potassium channel Kcnh4a, which is expressed in all Hcrt neurons, was silenced by the CRISPR-mediated gene inactivation system. The mutant kcnh4a (kcnh4a(-/-)) larvae showed reduced sleep time and consolidation, specifically during the night, suggesting that Kcnh4a regulates sleep.

Profiling molecular and behavioral circadian rhythms in the non-symbiotic sea anemone Nematostella vectensis.

Endogenous circadian clocks are poorly understood within early-diverging animal lineages. We have characterized circadian behavioral patterns and identified potential components of the circadian clock in the starlet sea anemone, Nematostella vectensis: a model cnidarian which lacks algal symbionts. Using automatic video tracking we showed that Nematostella exhibits rhythmic circadian locomotor activity, which is persistent in constant dark, shifted or disrupted by external dark/light cues and maintained the same rate at two different temperatures. This activity was inhibited by a casein kinase 1δ/ε inhibitor, suggesting a role for CK1 homologue(s) in Nematostella clock. Using high-throughput sequencing we profiled Nematostella transcriptomes over 48 hours under a light-dark cycle. We identified 180 Nematostella diurnally-oscillated transcripts and compared them with previously established databases of adult and larvae of the symbiotic coral Acropora millepora, revealing both shared homologues and unique rhythmic genes. Taken together, this study further establishes Nematostella as a non-symbiotic model organism to study circadian rhythms and increases our understanding about the fundamental elements of circadian regulation and their evolution within the Metazoa.

DREAM: a webserver for the identification of editing sites in mature miRNAs using deep sequencing data.

DREAM: detecting RNA editing associated with microRNAs, is a webserver for the identification of mature microRNA editing events using deep sequencing data. Raw microRNA sequencing reads can be provided as input, the reads are aligned against the genome and custom scripts process the data, search for potential editing sites and assess the statistical significance of the findings. The output is a text file with the location and the statistical description of all the putative editing sites detected.

The majority of transcripts in the squid nervous system are extensively recoded by A-to-I RNA editing.

RNA editing by adenosine deamination alters genetic information from the genomic blueprint. When it recodes mRNAs, it gives organisms the option to express diverse, functionally distinct, protein isoforms. All eumetazoans, from cnidarians to humans, express RNA editing enzymes. However, transcriptome-wide screens have only uncovered about 25 transcripts harboring conserved recoding RNA editing sites in mammals and several hundred recoding sites in Drosophila. These studies on few established models have led to the general assumption that recoding by RNA editing is extremely rare. Here we employ a novel bioinformatic approach with extensive validation to show that the squid Doryteuthis pealeii recodes proteins by RNA editing to an unprecedented extent. We identify 57,108 recoding sites in the nervous system, affecting the majority of the proteins studied. Recoding is tissue-dependent, and enriched in genes with neuronal and cytoskeletal functions, suggesting it plays an important role in brain physiology.