A multiplexed, single-cell sequencing screen identifies compounds that increase neurogenic reprogramming of murine Muller glia

Retinal degeneration in mammals causes permanent loss of vision, due to an inability to regenerate naturally. Some non-mammalian vertebrates show robust regeneration, via Muller glia (MG). We have recently made significant progress in stimulating adult mouse MG to regenerate functional neurons by transgenic expression of the proneural transcription factor Ascl1. While these results showed that MG can serve as an endogenous source of neuronal replacement, the efficacy of this process is limited. With the goal of improving this in mammals, we designed a small molecule screen using sci-Plex, a method to multiplex up to thousands of single nucleus RNA-seq conditions into a single experiment. We used this technology to screen a library of 92 compounds, identified, and validated two that promote neurogenesis in vivo. Our results demonstrate that high-throughput single-cell molecular profiling can substantially improve the discovery process for molecules and pathways that can stimulate neural regeneration and further demonstrate the potential for this approach to restore vision in patients with retinal disease.

Single-cell sequencing of individual retinal organoids reveals determinants of cell-fate heterogeneity

With a critical need for more complete in vitro models of human development and disease, organoids hold immense potential. Their complex cellular composition makes single-cell sequencing of great utility; however, the limitation of current technologies to a handful of treatment conditions restricts their use in screens or studies of organoid heterogeneity. Here, we apply sci-Plex, a single-cell combinatorial indexing (sci)-based RNA sequencing (RNA-seq) multiplexing method to retinal organoids. We demonstrate that sci-Plex and 10× methods produce highly concordant cell-class compositions and then expand sci-Plex to analyze the cell-class composition of 410 organoids upon modulation of critical developmental pathways. Leveraging individual organoid data, we develop a method to measure organoid heterogeneity, and we identify that activation of Wnt signaling early in retinal organoid cultures increases retinal cell classes up to 6 weeks later. Our data show sci-Plex's potential to dramatically scale up the analysis of treatment conditions on relevant human models.

Single-cell sequencing of individual retinal organoids reveals determinants of cell fate heterogeneity

With a critical need for more complete in vitro models of human development and disease, organoids hold immense potential. Their complex cellular composition makes single-cell sequencing of great utility; however, the limitation of current technologies to a handful of treatment conditions restricts their use in screens or studies of organoid heterogeneity. Here, we apply sci-Plex, a single-cell combinatorial indexing (sci)-based RNA-seq multiplexing method to retinal organoids. We demonstrate that sci-Plex and 10x methods produce highly concordant cell class compositions and then expand sci-Plex to analyze the cell class composition of 410 organoids upon modulation of critical developmental pathways. Leveraging individual organoid data, we develop a method to measure organoid heterogeneity, and we identify that activation of Wnt signaling early in retinal organoid cultures increases retinal cell classes up to six weeks later. Our data show sci-Plex's potential to dramatically scale-up the analysis of treatment conditions on relevant human models.

Long undecoded transcript isoform (LUTI) detection in meiotic budding yeast by direct RNA and transcript leader sequencing

LUTIs (Long Undecoded Transcript Isoforms) are 5'-extended and poorly translated mRNAs that can downregulate transcription from promoters more proximal to a gene's coding sequence (CDS). In this protocol, polyA RNA is extracted from budding yeast cells undergoing highly synchronized meiosis. Using a combination of long-read direct RNA sequencing and transcript leader sequencing (TL-seq), meiosis-specific LUTIs are systematically identified. Following identification, TL-seq is used to quantify the abundance of both LUTI and the more canonical gene-proximal (PROX) transcripts. For complete details on the use and execution of this protocol, please refer to Tresenrider et al. (2021).

Exploring whole-genome duplicate gene retention with complex genetic interaction analysis

Whole-genome duplication has played a central role in the genome evolution of many organisms, including the human genome. Most duplicated genes are eliminated, and factors that influence the retention of persisting duplicates remain poorly understood. We describe a systematic complex genetic interaction analysis with yeast paralogs derived from the whole-genome duplication event. Mapping of digenic interactions for a deletion mutant of each paralog, and of trigenic interactions for the double mutant, provides insight into their roles and a quantitative measure of their functional redundancy. Trigenic interaction analysis distinguishes two classes of paralogs: a more functionally divergent subset and another that retained more functional overlap. Gene feature analysis and modeling suggest that evolutionary trajectories of duplicated genes are dictated by combined functional and structural entanglement factors.

Evidence for an Integrated Gene Repression Mechanism Based on mRNA Isoform Toggling in Human Cells

We recently described an unconventional mode of gene regulation in budding yeast by which transcriptional and translational interference collaborate to down-regulate protein expression. Developmentally timed transcriptional interference inhibited production of a well translated mRNA isoform and resulted in the production of an mRNA isoform containing inhibitory upstream open reading frames (uORFs) that prevented translation of the main ORF. Transcriptional interference and uORF-based translational repression are established mechanisms outside of yeast, but whether this type of integrated regulation was conserved was unknown. Here we find that, indeed, a similar type of regulation occurs at the locus for the human oncogene MDM2 We observe evidence of transcriptional interference between the two MDM2 promoters, which produce a poorly translated distal promoter-derived uORF-containing mRNA isoform and a well-translated proximal promoter-derived transcript. Down-regulation of distal promoter activity markedly up-regulates proximal promoter-driven expression and results in local reduction of histone H3K36 trimethylation. Moreover, we observe that this transcript toggling between the two MDM2 isoforms naturally occurs during human embryonic stem cell differentiation programs.

Transcription of a 5' extended mRNA isoform directs dynamic chromatin changes and interference of a downstream promoter

Cell differentiation programs require dynamic regulation of gene expression. During meiotic prophase in Saccharomyces cerevisiae, expression of the kinetochore complex subunit Ndc80 is downregulated by a 5' extended long undecoded NDC80 transcript isoform. Here we demonstrate a transcriptional interference mechanism that is responsible for inhibiting expression of the coding NDC80 mRNA isoform. Transcription from a distal NDC80 promoter directs Set1-dependent histone H3K4 dimethylation and Set2-dependent H3K36 trimethylation to establish a repressive chromatin state in the downstream canonical NDC80 promoter. As a consequence, NDC80 expression is repressed during meiotic prophase. The transcriptional mechanism described here is rapidly reversible, adaptable to fine-tune gene expression, and relies on Set2 and the Set3 histone deacetylase complex. Thus, expression of a 5' extended mRNA isoform causes transcriptional interference at the downstream promoter. We demonstrate that this is an effective mechanism to promote dynamic changes in gene expression during cell differentiation.

Kinetochore inactivation by expression of a repressive mRNA

Differentiation programs such as meiosis depend on extensive gene regulation to mediate cellular morphogenesis. Meiosis requires transient removal of the outer kinetochore, the complex that connects microtubules to chromosomes. How the meiotic gene expression program temporally restricts kinetochore function is unknown. We discovered that in budding yeast, kinetochore inactivation occurs by reducing the abundance of a limiting subunit, Ndc80. Furthermore, we uncovered an integrated mechanism that acts at the transcriptional and translational level to repress NDC80 expression. Central to this mechanism is the developmentally controlled transcription of an alternate NDC80 mRNA isoform, which itself cannot produce protein due to regulatory upstream ORFs in its extended 5' leader. Instead, transcription of this isoform represses the canonical NDC80 mRNA expression in cis, thereby inhibiting Ndc80 protein synthesis. This model of gene regulation raises the intriguing notion that transcription of an mRNA, despite carrying a canonical coding sequence, can directly cause gene repression.