The exon quantification pipeline (EQP): a comprehensive approach to the quantification of gene, exon and junction expression from RNA-seq data

Photoaffinity enabled transcriptome-wide identification of splice modulating small molecule-RNA binding events in native cells

Splice modulating small molecules have been developed to promote the U1 snRNP to engage with pre-mRNAs with strong and altered sequence preference. Transcriptomic profiling of bulk RNA from compound treated cells enables detection of RNAs impacted; however, it is difficult to delineate whether transcriptional changes are a consequence of direct compound treatment or trans-acting effects. To identify RNA targets that bind directly with splice modulating compounds, we deployed a photoaffinity labeling (PAL)-based Chem-CLIP approach. Through this workflow, we identify the telomerase lncRNA (TERC) as a previously unknown target of this class of clinically relevant small molecules. Using cellular ΔSHAPE-MaP, we orthogonally validate and further define the compound binding site as likely to be the conserved CR4/5 domain. Additionally, a thorough analysis of the PAL-based Chem-CLIP data reveals that considering competed RNAs, irrespective of magnitude of enrichment, adds a rich dimension of hit calling.

Sulfinyl Aziridines as Stereoselective Covalent Destabilizing Degraders of the Oncogenic Transcription Factor MYC

While MYC is a significant oncogenic transcription factor driver of cancer, directly targeting MYC has remained challenging due to its intrinsic disorder and poorly defined structure, deeming it "undruggable." Whether transient pockets formed within intrinsically disordered and unstructured regions of proteins can be selectively targeted with small molecules remains an outstanding challenge. Here, we developed a bespoke stereochemically-paired spirocyclic oxindole aziridine covalent library and screened this library for degradation of MYC. Through this screen, we identified a hit covalent ligand KL2-236, bearing a unique sulfinyl aziridine warhead, that engaged MYC in vitro as pure MYC/MAX protein complex and in situ in cancer cells to destabilize MYC, inhibit MYC transcriptional activity and degrade MYC in a proteasome-dependent manner through targeting intrinsically disordered C203 and D205 residues. Notably, this reactivity was most pronounced for specific stereoisomers of KL2-236 with a diastereomer KL4-019 that was largely inactive. Mutagenesis of both C203 and D205 completely attenuated KL2-236-mediated MYC degradation. We have also optimized our initial KL2-236 hit compound to generate a more durable MYC degrader KL4-219A in cancer cells. Our results reveal a novel ligandable site within MYC and indicate that certain intrinsically disordered regions within high-value protein targets, such as MYC, can be interrogated by isomerically unique chiral small molecules, leading to destabilization and degradation.

Association of Genetically Predicted Skipping of COL4A4 Exon 27 with Hematuria and Albuminuria

Key Points

Using transcriptome-wide association studies, we identified an association between splicing out of exon 27 of COL4A4 and hematuria. We confirmed the presence of COL4A4 exon 27 splicing in an independent cohort. Functional assays revealed that the COL4A4 transcript with exon 27 spliced out affects collagen IV trimer assembly and secretion.

Background

Hematuria is an established sign of glomerular disease and can be associated with kidney failure, but there has been limited scientific study of this trait.

Methods

Here, we combined genetic data from the UK Biobank with predicted gene expression and splicing from Genotype Tissue Expression kidney cortex samples (n =65) in a transcriptome-wide association study to identify additional potential biological mechanisms influencing hematuria.

Results

The transcriptome-wide association study using kidney cortex identified significant associations for five genes in expression and three significant splicing events. Notably, we identified an association between the skipping of COL4A4 exon 27, which is genetically predicted by intronic rs11898094 (minor allele frequency 13%), and hematuria. Association between this variant was also found with urinary albumin excretion. We found independent evidence supporting the same variant predicting this skipping event in glomeruli-derived mRNA transcriptomics data (n =245) from the Nephrotic Syndrome Study Network. The functional significance of loss of exon 27 was demonstrated using the split NanoLuc-based α 3α 4α 5(IV) heterotrimer assay, in which type IV collagen heterotrimer formation was quantified by luminescence. The causal splicing variant for this skipping event is yet to be identified.

Conclusions

In summary, by integrating multiple data types, we identified a potential splicing event associated with hematuria and albuminuria.

IL-26 Potentiates Type 2 Skin Inflammation in the Presence of IL-1β

Atopic dermatitis (AD) is a debilitating inflammatory skin disorder. Biologics targeting the IL-4/IL-13 axis are effective in AD, but there is still a large proportion of patients who do not respond to IL-4R blockade. Further exploration of potentially pathogenic T-cell-derived cytokines in AD may lead to new effective treatments. This study aimed to investigate the downstream effects of IL-26 on skin in the context of type 2 skin inflammation. We found that IL-26 alone exhibited limited inflammatory activity in the skin. However, in the presence of IL-1β, IL-26 potentiated the secretion of TSLP, CXCL1, and CCL20 from human epidermis through Jak/signal transducer and activator of transcription signaling. Moreover, in an in vivo AD-like skin inflammation model, IL-26 exacerbated skin pathology and locally increased type 2 cytokines, most notably of IL13 in skin T helper cells. Neutralization of IL-1β abrogated IL-26-mediated effects, indicating that the presence of IL-1β is required for full IL-26 downstream action in vivo. These findings suggest that the presence of IL-1β enables IL-26 to be a key amplifier of inflammation in the skin. As such, IL-26 may contribute to the development and pathogenesis of inflammatory skin disorders such as AD.

The transcription factor ZNF469 regulates collagen production in liver fibrosis

Non-alcoholic fatty liver disease (NAFLD) - characterized by excess accumulation of fat in the liver - now affects one third of the world's population. As NAFLD progresses, extracellular matrix components including collagen accumulate in the liver causing tissue fibrosis, a major determinant of disease severity and mortality. To identify transcriptional regulators of fibrosis, we computationally inferred the activity of transcription factors (TFs) relevant to fibrosis by profiling the matched transcriptomes and epigenomes of 108 human liver biopsies from a deeply-characterized cohort of patients spanning the full histopathologic spectrum of NAFLD. CRISPR-based genetic knockout of the top 100 TFs identified ZNF469 as a regulator of collagen expression in primary human hepatic stellate cells (HSCs). Gain- and loss-of-function studies established that ZNF469 regulates collagen genes and genes involved in matrix homeostasis through direct binding to gene bodies and regulatory elements. By integrating multiomic large-scale profiling of human biopsies with extensive experimental validation we demonstrate that ZNF469 is a transcriptional regulator of collagen in HSCs. Overall, these data nominate ZNF469 as a previously unrecognized determinant of NAFLD-associated liver fibrosis.

Molecular characterization of chronic cutaneous wounds reveals subregion- and wound type-specific differential gene expression

A limited understanding of the pathology underlying chronic wounds has hindered the development of effective diagnostic markers and pharmaceutical interventions. This study aimed to elucidate the molecular composition of various common chronic ulcer types to facilitate drug discovery strategies. We conducted a comprehensive analysis of leg ulcers (LUs), encompassing venous and arterial ulcers, foot ulcers (FUs), pressure ulcers (PUs), and compared them with surgical wound healing complications (WHCs). To explore the pathophysiological mechanisms and identify similarities or differences within wounds, we dissected wounds into distinct subregions, including the wound bed, border, and peri-wound areas, and compared them against intact skin. By correlating histopathology, RNA sequencing (RNA-Seq), and immunohistochemistry (IHC), we identified unique genes, pathways, and cell type abundance patterns in each wound type and subregion. These correlations aim to aid clinicians in selecting targeted treatment options and informing the design of future preclinical and clinical studies in wound healing. Notably, specific genes, such as PITX1 and UPP1, exhibited exclusive upregulation in LUs and FUs, potentially offering significant benefits to specialists in limb preservation and clinical treatment decisions. In contrast, comparisons between different wound subregions, regardless of wound type, revealed distinct expression profiles. The pleiotropic chemokine-like ligand GPR15L (C10orf99) and transmembrane serine proteases TMPRSS11A/D were significantly upregulated in wound border subregions. Interestingly, WHCs exhibited a nearly identical transcriptome to PUs, indicating clinical relevance. Histological examination revealed blood vessel occlusions with impaired angiogenesis in chronic wounds, alongside elevated expression of genes and immunoreactive markers related to blood vessel and lymphatic epithelial cells in wound bed subregions. Additionally, inflammatory and epithelial markers indicated heightened inflammatory responses in wound bed and border subregions and reduced wound bed epithelialization. In summary, chronic wounds from diverse anatomical sites share common aspects of wound pathophysiology but also exhibit distinct molecular differences. These unique molecular characteristics present promising opportunities for drug discovery and treatment, particularly for patients suffering from chronic wounds. The identified diagnostic markers hold the potential to enhance preclinical and clinical trials in the field of wound healing.

Molecular mechanisms of atrazine toxicity on H19-7 hippocampal neurons revealed by integrated miRNA and mRNA "omics"

Atrazine (ATR) is a widely applied herbicide in Asia and South America with slow natural degradation and documented deleterious effects on human and animal health, including hippocampal toxicity. However, relatively little is known about the molecular mechanisms responsible for ATR-induced hippocampal damage. Screening for differentially expressed mRNAs and microRNAs (miRNAs), and construction of potential miRNA-mRNA regulatory networks can reveal such mechanisms, so we analyzed the mRNA and miRNA expression profiles of rat hippocampus-derived H19-7 cells in response to ATR (500 μM) and conducted Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes enrichment (KEGG) analyses. Integration of miRNA sequencing (miRNA-seq) and mRNA sequencing (mRNA-seq) results identified 114 differentially expressed miRNAs (DEMIs, 40 upregulated and 74 downregulated), and 510 differentially expressed mRNAs (DEMs, 177 upregulated and 333 downregulated) targeted by these DEMIs. The top 10 hub mRNAs (Fos, Prkcb, Ncf1, Vcam1, Atf3, Pak3, Pak1, Cacna1s, Junb, and Ccl2) and 19 related miRNAs (rno-miR-194-5p, rno-miR-24-3p, rno-miR-3074, rno-miR-1949, rno-miR-218a-1-3p, rno-miR-1843a-5p, rno-miR-1843b-5p, rno-miR-296-3p, rno-miR-320-3p, rno-miR-219a-1-3p, rno-miR-122-5p, rno-miR-1839-5p, rno-miR-1843a-3p, rno-miR-215, rno-miR-3583-3p, rno-miR-194-3p, rno-miR-128-1-5p, rno-miR-1956-5p, and rno-miR-466b-2-3p) were validated by quantitative real-time PCR. GO analysis indicated that these DEMs were enriched in genes associated with synaptic plasticity and antioxidant capacity, while KEGG analysis suggested that enriched DEMs were involved in calcium signaling, axon guidance, MAPK signaling, and glial carcinogenesis. The miRNA-mRNA regulatory network identified here may provide potential biomarkers and novel strategies for the treatment of hippocampal neurotoxicity induced by ATR.

Identification of marker genes to monitor residual iPSCs in iPSC-derived products

BACKGROUND

Engineered tissues and cell therapies based on human induced pluripotent stem cells (iPSCs) represent a promising approach for novel medicines. However, iPSC-derived cells and tissues may contain residual undifferentiated iPSCs that could lead to teratoma formation after implantation into patients. As a consequence, highly sensitive and specific methods for detecting residual undifferentiated iPSCs are indispensable for safety evaluations of iPSC-based therapies. The present study provides an approach for identifying potential marker genes for iPSC impurities in iPSC-derived cells using RNA sequencing data from iPSCs and various differentiated cell types.

METHODS

Identifying iPSC marker genes for each cell type individually provided a larger and more specific set of potential marker genes than considering all cell types in the analysis. Thus, the authors focused on identifying markers for iPSC impurities in iPSC-derived cardiomyocytes (iCMs) and validated the selected genes by reverse transcription quantitative polymerase chain reaction. The sensitivity of the candidate genes was determined by spiking different amounts of iPSCs into iCMs and their performance was compared with the previously suggested marker lin-28 homolog A (LIN28A).

RESULTS

Embryonic stem cell-related gene (ESRG), long intergenic non-protein coding RNA 678 (LINC00678), CaM kinase-like vesicle-associated (CAMKV), indoleamine 2,3-dioxygenase 1 (IDO1), chondromodulin (CNMD), LINE1-type transposase domain containing 1 (L1DT1), LIN28A, lymphocyte-specific protein tyrosine kinase (LCK), vertebrae development-associated (VRTN) and zinc finger and SCAN domain containing 10 (ZSCAN10) detected contaminant iPSCs among iCMs with a limit of detection that ranged from 0.001% to 0.1% depending on the gene and iCM batch used.

CONCLUSIONS

Using the example of iCMs, the authors provide a strategy for identifying a set of highly specific and sensitive markers that can be used for quality assessment of iPSC-derived products.

Transcriptional profiling of hepatocytes infected with the replicative form of the malaria parasite Plasmodium cynomolgi

BACKGROUND

The zoonotic simian parasite Plasmodium cynomolgi develops into replicating schizonts and dormant hypnozoites during the infection of hepatocytes and is used as a model organism to study relapsing malaria. The transcriptional profiling of P. cynomolgi liver stages was previously reported and revealed many important biological features of the parasite but left out the host response to malaria infection.

METHODS

Previously published RNA sequencing data were used to quantify the expression of host genes in rhesus macaque hepatocytes infected with P. cynomolgi in comparison to either cells from uninfected samples or uninfected bystander cells.

RESULTS

Although the dataset could not be used to resolve the transcriptional profile of hypnozoite-infected hepatocytes, it provided a snapshot of the host response to liver stage schizonts at 9-10 day post-infection and identified specific host pathways that are modulated during the exo-erythrocytic stage of P. cynomolgi.

CONCLUSIONS

This study constitutes a valuable resource characterizing the hepatocyte response to P. cynomolgi infection and provides a framework to build on future research that aims at understanding hepatocyte-parasite interactions during relapsing malaria infection.

An orally available, brain penetrant, small molecule lowers huntingtin levels by enhancing pseudoexon inclusion

Huntington's Disease (HD) is a progressive neurodegenerative disorder caused by CAG trinucleotide repeat expansions in exon 1 of the huntingtin (HTT) gene. The mutant HTT (mHTT) protein causes neuronal dysfunction, causing progressive motor, cognitive and behavioral abnormalities. Current treatments for HD only alleviate symptoms, but cerebral spinal fluid (CSF) or central nervous system (CNS) delivery of antisense oligonucleotides (ASOs) or virus vectors expressing RNA-induced silencing (RNAi) moieties designed to induce mHTT mRNA lowering have progressed to clinical trials. Here, we present an alternative disease modifying therapy the orally available, brain penetrant small molecule branaplam. By promoting inclusion of a pseudoexon in the primary transcript, branaplam lowers mHTT protein levels in HD patient cells, in an HD mouse model and in blood samples from Spinal Muscular Atrophy (SMA) Type I patients dosed orally for SMA (NCT02268552). Our work paves the way for evaluating branaplam's utility as an  HD therapy, leveraging small molecule splicing modulators to reduce expression of dominant disease genes by driving pseudoexon inclusion.

STING regulates peripheral nerve regeneration and colony stimulating factor 1 receptor (CSF1R) processing in microglia

Inflammatory responses are crucial for regeneration following peripheral nerve injury (PNI). PNI triggers inflammatory responses at the site of injury. The DNA-sensing receptor cyclic GMP-AMP synthase (cGAS) and its downstream effector stimulator of interferon genes (STING) sense foreign and self-DNA and trigger type I interferon (IFN) immune responses. We demonstrate here that following PNI, the cGAS/STING pathway is upregulated in the sciatic nerve of naive rats and dysregulated in old rats. In a nerve crush mouse model where STING is knocked out, myelin content in sciatic nerve is increased resulting in accelerated functional axon recovery. STING KO mice have lower macrophage number in sciatic nerve and decreased microglia activation in spinal cord 1 week post injury. STING activation regulated processing of colony stimulating factor 1 receptor (CSF1R) and microglia survival in vitro. Taking together, these data highlight a previously unrecognized role of STING in the regulation of nerve regeneration.

•

The cGAS/STING pathway is upregulated in sciatic nerve post nerve injury and in aging

•

STING ablation increases myelin content and accelerates functional axon recovery

•

STING KO mice reduces macrophage number in sciatic nerve and microglia activation post injury

Single-cell and bulk transcriptomics of the liver reveals potential targets of NASH with fibrosis

Fibrosis is characterized by the excessive production of collagen and other extracellular matrix (ECM) components and represents a leading cause of morbidity and mortality worldwide. Previous studies of nonalcoholic steatohepatitis (NASH) with fibrosis were largely restricted to bulk transcriptome profiles. Thus, our understanding of this disease is limited by an incomplete characterization of liver cell types in general and hepatic stellate cells (HSCs) in particular, given that activated HSCs are the major hepatic fibrogenic cell population. To help fill this gap, we profiled 17,810 non-parenchymal cells derived from six healthy human livers. In conjunction with public single-cell data of fibrotic/cirrhotic human livers, these profiles enable the identification of potential intercellular signaling axes (e.g., ITGAV-LAMC1, TNFRSF11B-VWF and NOTCH2-DLL4) and master regulators (e.g., RUNX1 and CREB3L1) responsible for the activation of HSCs during fibrogenesis. Bulk RNA-seq data of NASH patient livers and rodent models for liver fibrosis of diverse etiologies allowed us to evaluate the translatability of candidate therapeutic targets for NASH-related fibrosis. We identified 61 liver fibrosis-associated genes (e.g., AEBP1, PRRX1 and LARP6) that may serve as a repertoire of translatable drug target candidates. Consistent with the above regulon results, gene regulatory network analysis allowed the identification of CREB3L1 as a master regulator of many of the 61 genes. Together, this study highlights potential cell-cell interactions and master regulators that underlie HSC activation and reveals genes that may represent prospective hallmark signatures for liver fibrosis.

ZNRF3 and RNF43 cooperate to safeguard metabolic liver zonation and hepatocyte proliferation

AXIN2 and LGR5 mark intestinal stem cells (ISCs) that require WNT/β-Catenin signaling for constant homeostatic proliferation. In contrast, AXIN2/LGR5+ pericentral hepatocytes show low proliferation rates despite a WNT/β-Catenin activity gradient required for metabolic liver zonation. The mechanisms restricting proliferation in AXIN2+ hepatocytes and metabolic gene expression in AXIN2+ ISCs remained elusive. We now show that restricted chromatin accessibility in ISCs prevents the expression of β-Catenin-regulated metabolic enzymes, whereas fine-tuning of WNT/β-Catenin activity by ZNRF3 and RNF43 restricts proliferation in chromatin-permissive AXIN2+ hepatocytes, while preserving metabolic function. ZNRF3 deletion promotes hepatocyte proliferation, which in turn becomes limited by RNF43 upregulation. Concomitant deletion of RNF43 in ZNRF3 mutant mice results in metabolic reprogramming of periportal hepatocytes and induces clonal expansion in a subset of hepatocytes, ultimately promoting liver tumors. Together, ZNRF3 and RNF43 cooperate to safeguard liver homeostasis by spatially and temporally restricting WNT/β-Catenin activity, balancing metabolic function and hepatocyte proliferation.

CYP27A1-dependent anti-melanoma activity of limonoid natural products targets mitochondrial metabolism

Three limonoid natural products with selective anti-proliferative activity against BRAF(V600E) and NRAS(Q61K)-mutation-dependent melanoma cell lines were identified. Differential transcriptome analysis revealed dependency of compound activity on expression of the mitochondrial cytochrome P450 oxidase CYP27A1, a transcriptional target of melanogenesis-associated transcription factor (MITF). We determined that CYP27A1 activity is necessary for the generation of a reactive metabolite that proceeds to inhibit cellular proliferation. A genome-wide small interfering RNA screen in combination with chemical proteomics experiments revealed gene-drug functional epistasis, suggesting that these compounds target mitochondrial biogenesis and inhibit tumor bioenergetics through a covalent mechanism. Our work suggests a strategy for melanoma-specific targeting by exploiting the expression of MITF target gene CYP27A1 and inhibiting mitochondrial oxidative phosphorylation in BRAF mutant melanomas.

AXIN2+ Pericentral Hepatocytes Have Limited Contributions to Liver Homeostasis and Regeneration

The existence of specialized liver stem cell populations, including AXIN2⁺ pericentral hepatocytes, that safeguard homeostasis and repair has been controversial. Here, using AXIN2 lineage tracing in BAC-transgenic mice, we confirm the regenerative potential of intestinal stem cells (ISCs) but find limited roles for pericentral hepatocytes in liver parenchyma homeostasis. Liver regrowth following partial hepatectomy is enabled by proliferation of hepatocytes throughout the liver, rather than by a pericentral population. Periportal hepatocyte injury triggers local repair as well as auxiliary proliferation in all liver zones. DTA-mediated ablation of AXIN2⁺ pericentral hepatocytes transiently disrupts this zone, which is reestablished by conversion of pericentral vein-juxtaposed glutamine synthetase (GS)^- hepatocytes into GS⁺ hepatocytes and by compensatory proliferation of hepatocytes across liver zones. These findings show hepatocytes throughout the liver can upregulate AXIN2 and LGR5 after injury and contribute to liver regeneration on demand, without zonal dominance by a putative pericentral stem cell population.

UTS2B Defines a Novel Enteroendocrine Cell Population and Regulates GLP-1 Secretion Through SSTR5 in Male Mice

The gut-pancreas axis plays a key role in the regulation of glucose homeostasis and may be therapeutically exploited to treat not only type 2 diabetes but also hypoglycemia and hyperinsulinemia. We identify a novel enteroendocrine cell type expressing the peptide hormone urotensin 2B (UTS2B). UTS2B inhibits glucagon-like peptide-1 (GLP-1) secretion in mouse intestinal crypts and organoids, not by signaling through its cognate receptor UTS2R but through the activation of the somatostatin receptor (SSTR) 5. Circulating UTS2B concentrations in mice are physiologically regulated during starvation, further linking this peptide hormone to metabolism. Furthermore, administration of UTS2B to starved mice demonstrates that it is capable of regulating blood glucose and plasma concentrations of GLP-1 and insulin in vivo. Altogether, our results identify a novel cellular source of UTS2B in the gut, which acts in a paracrine manner to regulate GLP-1 secretion through SSTR5. These findings uncover a fine-tuning mechanism mediated by a ligand-receptor pair in the regulation of gut hormone secretion, which can potentially be exploited to correct metabolic unbalance caused by overactivation of the gut-pancreas axis.

TORC1 inhibition enhances immune function and reduces infections in the elderly

Inhibition of the mechanistic target of rapamycin (mTOR) protein kinase extends life span and ameliorates aging-related pathologies including declining immune function in model organisms. The objective of this phase 2a randomized, placebo-controlled clinical trial was to determine whether low-dose mTOR inhibitor therapy enhanced immune function and decreased infection rates in 264 elderly subjects given the study drugs for 6 weeks. A low-dose combination of a catalytic (BEZ235) plus an allosteric (RAD001) mTOR inhibitor that selectively inhibits target of rapamycin complex 1 (TORC1) downstream of mTOR was safe and was associated with a significant (P = 0.001) decrease in the rate of infections reported by elderly subjects for a year after study drug initiation. In addition, we observed an up-regulation of antiviral gene expression and an improvement in the response to influenza vaccination in this treatment group. Thus, selective TORC1 inhibition has the potential to improve immune function and reduce infections in the elderly.

Age-Related Gene Expression Signature in Rats Demonstrate Early, Late, and Linear Transcriptional Changes from Multiple Tissues

To understand the changes in gene expression that occur as a result of age, which might create a permissive or causal environment for age-related diseases, we produce a multi-time point age-related gene expression signature (AGES) from liver, kidney, skeletal muscle, and hippocampus of rats, comparing 6-, 9-, 12-, 18-, 21-, 24-, and 27-month-old animals. We focus on genes that changed in one direction throughout the lifespan of the animal, either early in life (early logistic changes), at mid-age (mid-logistic), late in life (late-logistic), or linearly, throughout the lifespan of the animal. The pathways perturbed because of chronological age demonstrate organ-specific and more-global effects of aging and point to mechanisms that could potentially be counter-regulated pharmacologically to treat age-associated diseases. A small number of genes are regulated by aging in the same manner in every tissue, suggesting they may be more-universal markers of aging.

CellSIUS provides sensitive and specific detection of rare cell populations from complex single-cell RNA-seq data

We develop CellSIUS (Cell Subtype Identification from Upregulated gene Sets) to fill a methodology gap for rare cell population identification for scRNA-seq data. CellSIUS outperforms existing algorithms for specificity and selectivity for rare cell types and their transcriptomic signature identification in synthetic and complex biological data. Characterization of a human pluripotent cell differentiation protocol recapitulating deep-layer corticogenesis using CellSIUS reveals unrecognized complexity in human stem cell-derived cellular populations. CellSIUS enables identification of novel rare cell populations and their signature genes providing the means to study those populations in vitro in light of their role in health and disease.

YAP, but Not RSPO-LGR4/5, Signaling in Biliary Epithelial Cells Promotes a Ductular Reaction in Response to Liver Injury

Biliary epithelial cells (BECs) form bile ducts in the liver and are facultative liver stem cells that establish a ductular reaction (DR) to support liver regeneration following injury. Liver damage induces periportal LGR5+ putative liver stem cells that can form BEC-like organoids, suggesting that RSPO-LGR4/5-mediated WNT/β-catenin activity is important for a DR. We addressed the roles of this and other signaling pathways in a DR by performing a focused CRISPR-based loss-of-function screen in BEC-like organoids, followed by in vivo validation and single-cell RNA sequencing. We found that BECs lack and do not require LGR4/5-mediated WNT/β-catenin signaling during a DR, whereas YAP and mTORC1 signaling are required for this process. Upregulation of AXIN2 and LGR5 is required in hepatocytes to enable their regenerative capacity in response to injury. Together, these data highlight heterogeneity within the BEC pool, delineate signaling pathways involved in a DR, and clarify the identity and roles of injury-induced periportal LGR5+ cells.

Transcriptomic analysis reveals reduced transcriptional activity in the malaria parasite Plasmodium cynomolgi during progression into dormancy

Relapses of Plasmodium dormant liver hypnozoites compromise malaria eradication efforts. New radical cure drugs are urgently needed, yet the vast gap in knowledge of hypnozoite biology impedes drug discovery. We previously unraveled the transcriptome of 6 to 7 day-old P. cynomolgi liver stages, highlighting pathways associated with hypnozoite dormancy (Voorberg-van der Wel et al., 2017). We now extend these findings by transcriptome profiling of 9 to 10 day-old liver stage parasites, thus revealing for the first time the maturation of the dormant stage over time. Although progression of dormancy leads to a 10-fold decrease in transcription and expression of only 840 genes, including genes associated with housekeeping functions, we show that pathways involved in quiescence, energy metabolism and maintenance of genome integrity remain the prevalent pathways active in mature hypnozoites.

TRRAP is a central regulator of human multiciliated cell formation

Multiciliated cells (MCCs) function to promote directional fluid flow across epithelial tissues. Wang et al. show that TRRAP, a component of multiple histone acetyltransferase complexes, is required for airway MCC formation and regulates a network of genes involved in MCC differentiation and function.

The multiciliated cell (MCC) is an evolutionarily conserved cell type, which in vertebrates functions to promote directional fluid flow across epithelial tissues. In the conducting airway, MCCs are generated by basal stem/progenitor cells and act in concert with secretory cells to perform mucociliary clearance to expel pathogens from the lung. Studies in multiple systems, including Xenopus laevis epidermis, murine trachea, and zebrafish kidney, have uncovered a transcriptional network that regulates multiple steps of multiciliogenesis, ultimately leading to an MCC with hundreds of motile cilia extended from their apical surface, which beat in a coordinated fashion. Here, we used a pool-based short hairpin RNA screening approach and identified TRRAP, an essential component of multiple histone acetyltransferase complexes, as a central regulator of MCC formation. Using a combination of immunofluorescence, signaling pathway modulation, and genomic approaches, we show that (a) TRRAP acts downstream of the Notch2-mediated basal progenitor cell fate decision and upstream of Multicilin to control MCC differentiation; and (b) TRRAP binds to the promoters and regulates the expression of a network of genes involved in MCC differentiation and function, including several genes associated with human ciliopathies.

Genome-wide CRISPR screen for PARKIN regulators reveals transcriptional repression as a determinant of mitophagy

In mitophagy, damaged mitochondria are targeted for disposal by the autophagy machinery. PARKIN promotes signaling of mitochondrial damage to the autophagy machinery for engagement, and PARKIN mutations cause Parkinson’s disease, possibly because damaged mitochondria accumulate in neurons. Because regulation of PARKIN abundance and the impact on signaling are poorly understood, we performed a genetic screen to identify PARKIN abundance regulators. Both positive and negative regulators were identified and will help us to further understand mitophagy and Parkinson’s disease. We show that some of the identified genes negatively regulate PARKIN gene expression, which impacts signaling of mitochondrial damage in mitophagy. This link between transcriptional repression and mitophagy is also apparent in neurons in culture, bearing implications for disease.

PARKIN, an E3 ligase mutated in familial Parkinson’s disease, promotes mitophagy by ubiquitinating mitochondrial proteins for efficient engagement of the autophagy machinery. Specifically, PARKIN-synthesized ubiquitin chains represent targets for the PINK1 kinase generating phosphoS65-ubiquitin (pUb), which constitutes the mitophagy signal. Physiological regulation of PARKIN abundance, however, and the impact on pUb accumulation are poorly understood. Using cells designed to discover physiological regulators of PARKIN abundance, we performed a pooled genome-wide CRISPR/Cas9 knockout screen. Testing identified genes individually resulted in a list of 53 positive and negative regulators. A transcriptional repressor network including THAP11 was identified and negatively regulates endogenous PARKIN abundance. RNAseq analysis revealed the PARKIN-encoding locus as a prime THAP11 target, and THAP11 CRISPR knockout in multiple cell types enhanced pUb accumulation. Thus, our work demonstrates the critical role of PARKIN abundance, identifies regulating genes, and reveals a link between transcriptional repression and mitophagy, which is also apparent in human induced pluripotent stem cell-derived neurons, a disease-relevant cell type.

A comparative transcriptomic analysis of replicating and dormant liver stages of the relapsing malaria parasite Plasmodium cynomolgi

Plasmodium liver hypnozoites, which cause disease relapse, are widely considered to be the last barrier towards malaria eradication. The biology of this quiescent form of the parasite is poorly understood which hinders drug discovery. We report a comparative transcriptomic dataset of replicating liver schizonts and dormant hypnozoites of the relapsing parasite Plasmodium cynomolgi. Hypnozoites express only 34% of Plasmodium physiological pathways, while 91% are expressed in replicating schizonts. Few known malaria drug targets are expressed in quiescent parasites, but pathways involved in microbial dormancy, maintenance of genome integrity and ATP homeostasis were robustly expressed. Several transcripts encoding heavy metal transporters were expressed in hypnozoites and the copper chelator neocuproine was cidal to all liver stage parasites. This transcriptomic dataset is a valuable resource for the discovery of vaccines and effective treatments to combat vivax malaria.

Profiling RNA-Seq at multiple resolutions markedly increases the number of causal eQTLs in autoimmune disease

Genome-wide association studies have identified hundreds of risk loci for autoimmune disease, yet only a minority (~25%) share genetic effects with changes to gene expression (eQTLs) in immune cells. RNA-Seq based quantification at whole-gene resolution, where abundance is estimated by culminating expression of all transcripts or exons of the same gene, is likely to account for this observed lack of colocalisation as subtle isoform switches and expression variation in independent exons can be concealed. We performed integrative cis-eQTL analysis using association statistics from twenty autoimmune diseases (560 independent loci) and RNA-Seq data from 373 individuals of the Geuvadis cohort profiled at gene-, isoform-, exon-, junction-, and intron-level resolution in lymphoblastoid cell lines. After stringently testing for a shared causal variant using both the Joint Likelihood Mapping and Regulatory Trait Concordance frameworks, we found that gene-level quantification significantly underestimated the number of causal cis-eQTLs. Only 5.0-5.3% of loci were found to share a causal cis-eQTL at gene-level compared to 12.9-18.4% at exon-level and 9.6-10.5% at junction-level. More than a fifth of autoimmune loci shared an underlying causal variant in a single cell type by combining all five quantification types; a marked increase over current estimates of steady-state causal cis-eQTLs. Causal cis-eQTLs detected at different quantification types localised to discrete epigenetic annotations. We applied a linear mixed-effects model to distinguish cis-eQTLs modulating all expression elements of a gene from those where the signal is only evident in a subset of elements. Exon-level analysis detected disease-associated cis-eQTLs that subtly altered transcription globally across the target gene. We dissected in detail the genetic associations of systemic lupus erythematosus and functionally annotated the candidate genes. Many of the known and novel genes were concealed at gene-level (e.g. IKZF2, TYK2, LYST). Our findings are provided as a web resource.

A comprehensive assessment of RNA-seq protocols for degraded and low-quantity samples

BACKGROUND

RNA-sequencing (RNA-seq) has emerged as one of the most sensitive tool for gene expression analysis. Among the library preparation methods available, the standard poly(A) + enrichment provides a comprehensive, detailed, and accurate view of polyadenylated RNAs. However, on samples of suboptimal quality ribosomal RNA depletion and exon capture methods have recently been reported as better alternatives.

METHODS

We compared for the first time three commercial Illumina library preparation kits (TruSeq Stranded mRNA, TruSeq Ribo-Zero rRNA Removal, and TruSeq RNA Access) as representatives of these three different approaches using well-established human reference RNA samples from the MAQC/SEQC consortium on a wide range of input amounts (from 100 ng down to 1 ng) and degradation levels (intact, degraded, and highly degraded).

RESULTS

We assessed the accuracy of the generated expression values by comparison to gold standard TaqMan qPCR measurements and gained unprecedented insight into the limits of applicability in terms of input quantity and sample quality of each protocol. We found that each protocol generates highly reproducible results (R ² > 0.92) on intact RNA samples down to input amounts of 10 ng. For degraded RNA samples, Ribo-Zero showed clear performance advantages over the other two protocols as it generated more accurate and better reproducible gene expression results even at very low input amounts such as 1 ng and 2 ng. For highly degraded RNA samples, RNA Access performed best generating reliable data down to 5 ng input.

CONCLUSIONS

We found that the ribosomal RNA depletion protocol from Illumina works very well at amounts far below recommendation and over a good range of intact and degraded material. We also infer that the exome-capture protocol (RNA Access, Illumina) performs better than other methods on highly degraded and low amount samples.

Identification of oncogenic driver mutations by genome-wide CRISPR-Cas9 dropout screening

BACKGROUND

Genome-wide CRISPR-Cas9 dropout screens can identify genes whose knockout affects cell viability. Recent CRISPR screens detected thousands of essential genes required for cellular survival and key cellular processes; however discovering novel lineage-specific genetic dependencies from the many hits still remains a challenge.

RESULTS

To assess whether CRISPR-Cas9 dropout screens can help identify cancer dependencies, we screened two human cancer cell lines carrying known and distinct oncogenic mutations using a genome-wide sgRNA library. We found that the gRNA targeting the driver mutation EGFR was one of the highest-ranking candidates in the EGFR-mutant HCC-827 lung adenocarcinoma cell line. Likewise, sgRNAs for NRAS and MAP2K1 (MEK1), a downstream kinase of mutant NRAS, were identified among the top hits in the NRAS-mutant neuroblastoma cell line CHP-212. Depletion of these genes targeted by the sgRNAs strongly correlated with the sensitivity to specific kinase inhibitors of the EGFR or RAS pathway in cell viability assays. In addition, we describe other dependencies such as TBK1 in HCC-827 cells and TRIB2 in CHP-212 cells which merit further investigation.

CONCLUSIONS

We show that genome-wide CRISPR dropout screens are suitable for the identification of oncogenic drivers and other essential genes.