FANTOM5 CAGE profiles of human and mouse reprocessed for GRCh38 and GRCm38 genome assemblies

Mouse promoters are characterised by low occupancy and high turnover of RNA polymerase II

The general transcription machinery and its occupancy at promoters are highly conserved across metazoans. This contrasts with the kinetics of mRNA production that considerably differ between model species such as Drosophila and mouse. The molecular basis for these kinetic differences is currently unknown. Here, we used Single-Molecule Footprinting to measure RNA Polymerase II (Pol II) occupancy, the fraction of DNA molecules bound, at promoters in mouse and Drosophila cell lines. Single-molecule data reveals that Pol II occupancy is on average 3-5 times more frequent at transcriptionally active Drosophila promoters than active mouse promoters. Kinetic modelling of the occupancy states suggests that these differences in Pol II occupancy are determined by the ratio between the transcription initiation and Pol II turnover rates. We used chemical perturbation of transcription initiation to determine Pol II turnover rate in both species. Integration of these data into the model shows that infrequent Pol II occupancy in mouse is explained by the combination of high Pol II turnover and low transcription initiation rates.

Genetically regulated eRNA expression predicts chromatin contact frequency and reveals genetic mechanisms at GWAS loci

The biological functions of extragenic enhancer RNAs and their impact on disease risk remain relatively underexplored. In this work, we develop in silico models of genetically regulated expression of enhancer RNAs across 49 cell and tissue types, characterizing their degree of genetic control. Leveraging the estimated genetically regulated expression for enhancer RNAs and canonical genes in a large-scale DNA biobank (N > 70,000) and high-resolution Hi-C contact data, we train a deep learning-based model of pairwise three-dimensional chromatin contact frequency for enhancer-enhancer and enhancer-gene pairs in cerebellum and whole blood. Notably, the use of genetically regulated expression of enhancer RNAs provides substantial tissue-specific predictive power, supporting a role for these transcripts in modulating spatial chromatin organization. We identify schizophrenia-associated enhancer RNAs independent of GWAS loci using enhancer RNA-based TWAS and determine the causal effects of these enhancer RNAs using Mendelian randomization. Using enhancer RNA-based TWAS, we generate a comprehensive resource of tissue-specific enhancer associations with complex traits in the UK Biobank. Finally, we show that a substantially greater proportion (63%) of GWAS associations colocalize with causal regulatory variation when enhancer RNAs are included.

Identification of a SNAI1 enhancer RNA that drives cancer cell plasticity

Enhancer RNAs (eRNAs) are a pivotal class of enhancer-derived non-coding RNAs that drive gene expression. Here we identify the SNAI1 enhancer RNA (SNAI1e; SCREEM2) as a key activator of SNAI1 expression and a potent enforcer of transforming growth factor-β (TGF-β)/SMAD signaling in cancer cells. SNAI1e depletion impairs TGF-β-induced epithelial-mesenchymal transition (EMT), migration, in vivo extravasation, stemness, and chemotherapy resistance in breast cancer cells. SNAI1e functions as an eRNA to cis-regulate SNAI1 enhancer activity by binding to and strengthening the enrichment of the transcriptional co-activator bromodomain containing protein 4 (BRD4) at the local enhancer. SNAI1e selectively promotes the expression of SNAI1, which encodes the EMT transcription factor SNAI1. Furthermore, we reveal that SNAI1 interacts with and anchors the inhibitory SMAD7 in the nucleus, and thereby prevents TGF-β type I receptor (TβRI) polyubiquitination and proteasomal degradation. Our findings establish SNAI1e as a critical driver of SNAI1 expression and TGF-β-induced cell plasticity.

Hybrid exons evolved by coupling transcription initiation and splicing at the nucleotide level

Exons within transcripts are traditionally classified as first, internal or last exons, each governed by different regulatory mechanisms. We recently described the widespread usage of 'hybrid' exons that serve as terminal or internal exons in different transcripts. Here, we employ an interpretable deep learning pipeline to dissect the sequence features governing the co-regulation of transcription initiation and splicing in hybrid exons. Using ENCODE data from human tissues, we identified 80 000 hybrid first-internal exons. These exons often possess a relaxed chromatin state, allowing transcription initiation within the gene body. Interestingly, transcription start sites of hybrid exons are typically centered at the 3' splice site, suggesting tight coupling between splicing and transcription initiation. We identified two subcategories of hybrid exons: the majority resemble internal exons, maintaining strong 3' splice sites, while a minority show enrichment in promoter elements, resembling first exons. Diving into the evolution of their sequences, we found that human hybrid exons with orthologous first exons in other species usually gained 3' splice sites or whole exons upstream, while those with orthologous internal exons often gained promoter elements. Overall, our findings unveil the intricate regulatory landscape of hybrid exons and reveal stronger connections between transcription initiation and RNA splicing than previously acknowledged.

Arrayed CRISPR libraries for the genome-wide activation, deletion and silencing of human protein-coding genes

Arrayed CRISPR libraries extend the scope of gene-perturbation screens to non-selectable cell phenotypes. However, library generation requires assembling thousands of vectors expressing single-guide RNAs (sgRNAs). Here, by leveraging massively parallel plasmid-cloning methodology, we show that arrayed libraries can be constructed for the genome-wide ablation (19,936 plasmids) of human protein-coding genes and for their activation and epigenetic silencing (22,442 plasmids), with each plasmid encoding an array of four non-overlapping sgRNAs designed to tolerate most human DNA polymorphisms. The quadruple-sgRNA libraries yielded high perturbation efficacies in deletion (75-99%) and silencing (76-92%) experiments and substantial fold changes in activation experiments. Moreover, an arrayed activation screen of 1,634 human transcription factors uncovered 11 novel regulators of the cellular prion protein PrPC, screening with a pooled version of the ablation library led to the identification of 5 novel modifiers of autophagy that otherwise went undetected, and 'post-pooling' individually produced lentiviruses eliminated template-switching artefacts and enhanced the performance of pooled screens for epigenetic silencing. Quadruple-sgRNA arrayed libraries are a powerful and versatile resource for targeted genome-wide perturbations.

m6Am sequesters PCF11 to suppress premature termination and drive neuroblastoma differentiation

N6,2'-O-dimethyladenosine (m6Am) is an abundant mRNA modification that impacts multiple diseases, but its function remains controversial because the m6Am reader is unknown. Using quantitative proteomics, we identified transcriptional terminator premature cleavage factor II (PCF11) as a m6Am-specific reader in human cells. Direct quantification of mature versus nascent RNAs reveals that m6Am does not regulate mRNA stability but promotes nascent transcription. Mechanistically, m6Am functions by sequestering PCF11 away from proximal RNA polymerase II (RNA Pol II). This suppresses PCF11 from dissociating RNA Pol II near transcription start sites, thereby promoting full-length transcription of m6Am-modified RNAs. m6Am's unique relationship with PCF11 means m6Am function is enhanced when PCF11 is reduced, which occurs during all-trans-retinoic-acid (ATRA)-induced neuroblastoma-differentiation therapy. Here, m6Am promotes expression of ATF3, which represses neuroblastoma biomarker MYCN. Depleting m6Am suppresses MYCN repression in ATRA-treated neuroblastoma and maintains their tumor-stem-like properties. Collectively, we characterize m6Am as an anti-terminator RNA modification that suppresses premature termination and modulates neuroblastoma's therapeutic response.

Whole-genome sequencing identifies variants in ANK1, LRRN1, HAS1, and other genes and regulatory regions for stroke in type 1 diabetes

Individuals with type 1 diabetes (T1D) carry a markedly increased risk of stroke, with distinct clinical and neuroimaging characteristics as compared to those without diabetes. Using whole-exome or whole-genome sequencing of 1,051 individuals with T1D, we aimed to find rare and low-frequency genomic variants associated with stroke in T1D. We analysed the genome comprehensively with single-variant analyses, gene aggregate analyses, and aggregate analyses on genomic windows, enhancers and promoters. In addition, we attempted replication in T1D using a genome-wide association study (N = 3,945) and direct genotyping (N = 3,263), and in the general population from the large-scale population-wide FinnGen project and UK Biobank summary statistics. We identified a rare missense variant on SREBF1 exome-wide significantly associated with stroke (rs114001633, p.Pro227Leu, p-value = 7.30 × 10-8), which replicated for hemorrhagic stroke in T1D. Using gene aggregate analysis, we identified exome-wide significant genes: ANK1 and LRRN1 displayed replication evidence in T1D, and LRRN1, HAS1 and UACA in the general population (UK Biobank). Furthermore, we performed sliding-window analyses and identified 14 genome-wide significant windows for stroke on 4q33-34.1, of which two replicated in T1D, and a suggestive genomic window on LINC01500, which replicated in T1D. Finally, we identified a suggestively stroke-associated TRPM2-AS promoter (p-value = 5.78 × 10-6) with borderline significant replication in T1D, which we validated with an in vitro cell-based assay. Due to the rarity of the identified genetic variants, future replication of the genomic regions represented here is required with sequencing of individuals with T1D. Nevertheless, we here report the first genome-wide analysis on stroke in individuals with diabetes.

Lgr6-expressing functional nail stem-like cells differentiated from human-induced pluripotent stem cells

The nail matrix containing stem cell populations produces nails and may contribute to fingertip regeneration. Nails are important tissues that maintain the functions of the hand and foot for handling objects and locomotion. Tumor chemotherapy impairs nail growth and, in many cases, loses them, although not permanently. In this report, we have achieved the successful differentiation of nail stem (NS)-like cells from human-induced pluripotent stem cells (iPSCs) via digit organoids by stepwise stimulation, tracing the molecular processes involved in limb development. Comprehensive mRNA sequencing analysis revealed that the digit organoid global gene expression profile fits human finger development. The NS-like cells expressed Lgr6 mRNA and protein and produced type-I keratin, KRT17, and type-II keratin, KRT81, which are abundant in nails. Furthermore, we succeeded in producing functional Lgr6-reporter human iPSCs. The reporter iPSC-derived Lgr6-positive cells also produced KRT17 and KRT81 proteins in the percutaneously transplanted region. To the best of our knowledge, this is the first report of NS-like cell differentiation from human iPSCs. Our differentiation method and reporter construct enable the discovery of drugs for nail repair and possibly fingertip-regenerative therapy.

OVCH1 Antisense RNA 1 is differentially expressed between non-frail and frail old adults

While some old adults stay healthy and non-frail up to late in life, others experience multimorbidity and frailty often accompanied by a pro-inflammatory state. The underlying molecular mechanisms for those differences are still obscure. Here, we used gene expression analysis to understand the molecular underpinning between non-frail and frail individuals in old age. Twenty-four adults (50% non-frail and 50% frail) from InCHIANTI study were included. Total RNA extracted from whole blood was analyzed by Cap Analysis of Gene Expression (CAGE). CAGE identified transcription start site (TSS) and active enhancer regions. We identified a set of differentially expressed (DE) TSS and enhancer between non-frail and frail and male and female participants. Several DE TSSs were annotated as lncRNA (XIST and TTTY14) and antisense RNAs (ZFX-AS1 and OVCH1 Antisense RNA 1). The promoter region chr6:366,786,54-366,787,97;+ was DE and overlapping the longevity CDKN1A gene. GWAS-LD enrichment analysis identifies overlapping LD-blocks with the DE regions with reported traits in GWAS catalog (isovolumetric relaxation time and urinary tract infection frequency). Furthermore, we used weighted gene co-expression network analysis (WGCNA) to identify changes of gene expression associated with clinical traits and identify key gene modules. We performed functional enrichment analysis of the gene modules with significant trait/module correlation. One gene module is showing a very distinct pattern in hub genes. Glycogen Phosphorylase L (PYGL) was the top ranked hub gene between non-frail and frail. We predicted transcription factor binding sites (TFBS) and motif activity. TF involved in age-related pathways (e.g., FOXO3 and MYC) shows different expression patterns between non-frail and frail participants. Expanding the study of OVCH1 Antisense RNA 1 and PYGL may help understand the mechanisms leading to loss of homeostasis that ultimately causes frailty.

C19ORF84 connects piRNA and DNA methylation machineries to defend the mammalian germ line

In the male mouse germ line, PIWI-interacting RNAs (piRNAs), bound by the PIWI protein MIWI2 (PIWIL4), guide DNA methylation of young active transposons through SPOCD1. However, the underlying mechanisms of SPOCD1-mediated piRNA-directed transposon methylation and whether this pathway functions to protect the human germ line remain unknown. We identified loss-of-function variants in human SPOCD1 that cause defective transposon silencing and male infertility. Through the analysis of these pathogenic alleles, we discovered that the uncharacterized protein C19ORF84 interacts with SPOCD1. DNMT3C, the DNA methyltransferase responsible for transposon methylation, associates with SPOCD1 and C19ORF84 in fetal gonocytes. Furthermore, C19ORF84 is essential for piRNA-directed DNA methylation and male mouse fertility. Finally, C19ORF84 mediates the in vivo association of SPOCD1 with the de novo methylation machinery. In summary, we have discovered a conserved role for the human piRNA pathway in transposon silencing and C19ORF84, an uncharacterized protein essential for orchestrating piRNA-directed DNA methylation.

TE-TSS: an integrated data resource of human and mouse transposable element (TE)-derived transcription start site (TSS)

Transposable elements (TEs) are abundant in the genome and serve as crucial regulatory elements. Some TEs function as epigenetically regulated promoters, and these TE-derived transcription start sites (TSSs) play a crucial role in regulating genes associated with specific functions, such as cancer and embryogenesis. However, the lack of an accessible database that systematically gathers TE-derived TSS data is a current research gap. To address this, we established TE-TSS, an integrated data resource of human and mouse TE-derived TSSs (http://xozhanglab.com/TETSS). TE-TSS has compiled 2681 RNA sequencing datasets, spanning various tissues, cell lines and developmental stages. From these, we identified 5768 human TE-derived TSSs and 2797 mouse TE-derived TSSs, with 47% and 38% being experimentally validated, respectively. TE-TSS enables comprehensive exploration of TSS usage in diverse samples, providing insights into tissue-specific gene expression patterns and transcriptional regulatory elements. Furthermore, TE-TSS compares TE-derived TSS regions across 15 mammalian species, enhancing our understanding of their evolutionary and functional aspects. The establishment of TE-TSS facilitates further investigations into the roles of TEs in shaping the transcriptomic landscape and offers valuable resources for comprehending their involvement in diverse biological processes.

Neuronal membrane proteasome-derived peptides modulate NMDAR-dependent neuronal signaling to promote changes in gene expression

The neuronal membrane proteasome (NMP) degrades intracellular proteins into peptides that are released directly into the extracellular space, whereby they stimulate neurons to promote signaling mechanisms that remain unknown. Here, we demonstrate that neuronal stimulation promotes NMP activity and, subsequently, enhanced production of NMP peptides. We show that these neuronal activity-dependent NMP peptides can rapidly promote N-methyl-D-aspartate receptor (NMDAR)-dependent calcium influx in neurons. This leads to sustained phosphorylation of the well-defined stimulus-induced transcription factor, cyclic AMP response element (CRE)-binding protein (CREB). Downstream of these events, we identified changes to neuronal target genes which included increased expression of immediate early genes (e.g., Fos, Npas4, Egr4) and other genes known to have critical neuroregulatory roles. Further observations led to the discovery that NMP peptide-induced changes in gene expression is dependent on NMDARs and independent of AMPA receptors or voltage-gated sodium channels. These data demonstrate that NMP peptides are endogenous and selective activators of NMDA receptors and act as sufficient and novel stimuli within the context of neuronal activity-dependent signaling. This novel pathway is parallel to classic neuronal activity-dependent programs and points to NMP and its resulting peptides as potential modulators of neuronal function.

Minimum entropy framework identifies a novel class of genomic functional elements and reveals regulatory mechanisms at human disease loci

We introduce CoRE-BED, a framework trained using 19 epigenomic features in 33 major cell and tissue types to predict cell-type-specific regulatory function. CoRE-BED identifies nine functional classes de-novo, capturing both known and new regulatory categories. Notably, we describe a previously undercharacterized class that we term Development Associated Elements (DAEs), which are highly enriched in cell types with elevated regenerative potential and distinguished by the dual presence of either H3K4me2 and H3K9ac (an epigenetic signature associated with kinetochore assembly) or H3K79me3 and H4K20me1 (a signature associated with transcriptional pause release). Unlike bivalent promoters, which represent a transitory state between active and silenced promoters, DAEs transition directly to or from a non-functional state during stem cell differentiation and are proximal to highly expressed genes. CoRE-BED's interpretability facilitates causal inference and functional prioritization. Across 70 complex traits, distal insulators account for the largest mean proportion of SNP heritability (~49%) captured by the GWAS. Collectively, our results demonstrate the value of exploring non-conventional ways of regulatory classification that enrich for trait heritability, to complement existing approaches for cis-regulatory prediction.

Detection of isoforms and genomic alterations by high-throughput full-length single-cell RNA sequencing in ovarian cancer

Understanding the complex background of cancer requires genotype-phenotype information in single-cell resolution. Here, we perform long-read single-cell RNA sequencing (scRNA-seq) on clinical samples from three ovarian cancer patients presenting with omental metastasis and increase the PacBio sequencing depth to 12,000 reads per cell. Our approach captures 152,000 isoforms, of which over 52,000 were not previously reported. Isoform-level analysis accounting for non-coding isoforms reveals 20% overestimation of protein-coding gene expression on average. We also detect cell type-specific isoform and poly-adenylation site usage in tumor and mesothelial cells, and find that mesothelial cells transition into cancer-associated fibroblasts in the metastasis, partly through the TGF-β/miR-29/Collagen axis. Furthermore, we identify gene fusions, including an experimentally validated IGF2BP2::TESPA1 fusion, which is misclassified as high TESPA1 expression in matched short-read data, and call mutations confirmed by targeted NGS cancer gene panel results. With these findings, we envision long-read scRNA-seq to become increasingly relevant in oncology and personalized medicine.

eQTL Catalogue 2023: New datasets, X chromosome QTLs, and improved detection and visualisation of transcript-level QTLs

The eQTL Catalogue is an open database of uniformly processed human molecular quantitative trait loci (QTLs). We are continuously updating the resource to further increase its utility for interpreting genetic associations with complex traits. Over the past two years, we have increased the number of uniformly processed studies from 21 to 31 and added X chromosome QTLs for 19 compatible studies. We have also implemented Leafcutter to directly identify splice-junction usage QTLs in all RNA sequencing datasets. Finally, to improve the interpretability of transcript-level QTLs, we have developed static QTL coverage plots that visualise the association between the genotype and average RNA sequencing read coverage in the region for all 1.7 million fine mapped associations. To illustrate the utility of these updates to the eQTL Catalogue, we performed colocalisation analysis between vitamin D levels in the UK Biobank and all molecular QTLs in the eQTL Catalogue. Although most GWAS loci colocalised both with eQTLs and transcript-level QTLs, we found that visual inspection could sometimes be used to distinguish primary splicing QTLs from those that appear to be secondary consequences of large-effect gene expression QTLs. While these visually confirmed primary splicing QTLs explain just 6/53 of the colocalising signals, they are significantly less pleiotropic than eQTLs and identify a prioritised causal gene in 4/6 cases.

Nanopore sequencing unveils the complexity of the cold-activated murine brown adipose tissue transcriptome

Alternative transcription increases transcriptome complexity by expression of multiple transcripts per gene. Annotation and quantification of transcripts using short-read sequencing is non-trivial. Long-read sequencing aims at overcoming these problems by sequencing full-length transcripts. Activation of brown adipose tissue (BAT) thermogenesis involves major transcriptomic remodeling and positively affects metabolism via increased energy expenditure. We benchmark Oxford Nanopore Technology (ONT) long-read sequencing protocols to Illumina short-read sequencing assessing alignment characteristics, gene and transcript detection and quantification, differential gene and transcript expression, transcriptome reannotation, and differential transcript usage (DTU). We find ONT sequencing is superior to Illumina for transcriptome reassembly, reducing the risk of false-positive events by unambiguously mapping reads to transcripts. We identified novel isoforms of genes undergoing DTU in cold-activated BAT including Cars2, Adtrp, Acsl5, Scp2, Aldoa, and Pde4d, validated by real-time PCR. The reannotated murine BAT transcriptome established here provides a framework for future investigations into the regulation of BAT.

SCN1A: bioinformatically informed revised boundaries for promoter and enhancer regions

Pathogenic variations in the sodium voltage-gated channel alpha subunit 1 (SCN1A) gene are responsible for multiple epilepsy phenotypes, including Dravet syndrome, febrile seizures (FS) and genetic epilepsy with FS plus. Phenotypic heterogeneity is a hallmark of SCN1A-related epilepsies, the causes of which are yet to be clarified. Genetic variation in the non-coding regulatory regions of SCN1A could be one potential causal factor. However, a comprehensive understanding of the SCN1A regulatory landscape is currently lacking. Here, we summarized the current state of knowledge of SCN1A regulation, providing details on its promoter and enhancer regions. We then integrated currently available data on SCN1A promoters by extracting information related to the SCN1A locus from genome-wide repositories and clearly defined the promoter and enhancer regions of SCN1A. Further, we explored the cellular specificity of differential SCN1A promoter usage. We also reviewed and integrated the available human brain-derived enhancer databases and mouse-derived data to provide a comprehensive computationally developed summary of SCN1A brain-active enhancers. By querying genome-wide data repositories, extracting SCN1A-specific data and integrating the different types of independent evidence, we created a comprehensive catalogue that better defines the regulatory landscape of SCN1A, which could be used to explore the role of SCN1A regulatory regions in disease.

Systematic visualisation of molecular QTLs reveals variant mechanisms at GWAS loci

Splicing quantitative trait loci (QTLs) have been implicated as a common mechanism underlying complex trait associations. However, utilising splicing QTLs in target discovery and prioritisation has been challenging due to extensive data normalisation which often renders the direction of the genetic effect as well as its magnitude difficult to interpret. This is further complicated by the fact that strong expression QTLs often manifest as weak splicing QTLs and vice versa, making it difficult to uniquely identify the underlying molecular mechanism at each locus. We find that these ambiguities can be mitigated by visualising the association between the genotype and average RNA sequencing read coverage in the region. Here, we generate these QTL coverage plots for 1.7 million molecular QTL associations in the eQTL Catalogue identified with five quantification methods. We illustrate the utility of these QTL coverage plots by performing colocalisation between vitamin D levels in the UK Biobank and all molecular QTLs in the eQTL Catalogue. We find that while visually confirmed splicing QTLs explain just 6/53 of the colocalising signals, they are significantly less pleiotropic than eQTLs and identify a prioritised causal gene in 4/6 cases. All our association summary statistics and QTL coverage plots are freely available at https://www.ebi.ac.uk/eqtl/.

Pan-cancer analysis identifies tumor-specific antigens derived from transposable elements

Cryptic promoters within transposable elements (TEs) can be transcriptionally reactivated in tumors to create new TE-chimeric transcripts, which can produce immunogenic antigens. We performed a comprehensive screen for these TE exaptation events in 33 TCGA tumor types, 30 GTEx adult tissues and 675 cancer cell lines, and identified 1,068 TE-exapted candidates with the potential to generate shared tumor-specific TE-chimeric antigens (TS-TEAs). Whole-lysate and HLA-pulldown mass spectrometry data confirmed that TS-TEAs are presented on the surface of cancer cells. In addition, we highlight tumor-specific membrane proteins transcribed from TE promoters that constitute aberrant epitopes on the extracellular surface of cancer cells. Altogether, we showcase the high pan-cancer prevalence of TS-TEAs and atypical membrane proteins that could potentially be therapeutically exploited and targeted.

FAVOR: functional annotation of variants online resource and annotator for variation across the human genome

Large biobank-scale whole genome sequencing (WGS) studies are rapidly identifying a multitude of coding and non-coding variants. They provide an unprecedented resource for illuminating the genetic basis of human diseases. Variant functional annotations play a critical role in WGS analysis, result interpretation, and prioritization of disease- or trait-associated causal variants. Existing functional annotation databases have limited scope to perform online queries and functionally annotate the genotype data of large biobank-scale WGS studies. We develop the Functional Annotation of Variants Online Resources (FAVOR) to meet these pressing needs. FAVOR provides a comprehensive multi-faceted variant functional annotation online portal that summarizes and visualizes findings of all possible nine billion single nucleotide variants (SNVs) across the genome. It allows for rapid variant-, gene- and region-level queries of variant functional annotations. FAVOR integrates variant functional information from multiple sources to describe the functional characteristics of variants and facilitates prioritizing plausible causal variants influencing human phenotypes. Furthermore, we provide a scalable annotation tool, FAVORannotator, to functionally annotate large-scale WGS studies and efficiently store the genotype and their variant functional annotation data in a single file using the annotated Genomic Data Structure (aGDS) format, making downstream analysis more convenient. FAVOR and FAVORannotator are available at https://favor.genohub.org.

scDrug: From single-cell RNA-seq to drug response prediction

Single-cell RNA sequencing (scRNA-seq) technology allows massively parallel characterization of thousands of cells at the transcriptome level. scRNA-seq is emerging as an important tool to investigate the cellular components and their interactions in the tumor microenvironment. scRNA-seq is also used to reveal the association between tumor microenvironmental patterns and clinical outcomes and to dissect cell-specific effects of drug treatment in complex tissues. Recent advances in scRNA-seq have driven the discovery of biomarkers in diseases and therapeutic targets. Although methods for prediction of drug response using gene expression of scRNA-seq data have been proposed, an integrated tool from scRNA-seq analysis to drug discovery is required. We present scDrug as a bioinformatics workflow that includes a one-step pipeline to generate cell clustering for scRNA-seq data and two methods to predict drug treatments. The scDrug pipeline consists of three main modules: scRNA-seq analysis for identification of tumor cell subpopulations, functional annotation of cellular subclusters, and prediction of drug responses. scDrug enables the exploration of scRNA-seq data readily and facilitates the drug repurposing process. scDrug is freely available on GitHub at https://github.com/ailabstw/scDrug.

A high-resolution map of human RNA translation

Translated small open reading frames (smORFs) can have important regulatory roles and encode microproteins, yet their genome-wide identification has been challenging. We determined the ribosome locations across six primary human cell types and five tissues and detected 7,767 smORFs with translational profiles matching those of known proteins. The human genome was found to contain highly cell-type- and tissue-specific smORFs and a subset that encodes highly conserved amino acid sequences. Changes in the translational efficiency of upstream-encoded smORFs (uORFs) and the corresponding main ORFs predominantly occur in the same direction. Integration with 456 mass-spectrometry datasets confirms the presence of 603 small peptides at the protein level in humans and provides insights into the subcellular localization of these small proteins. This study provides a comprehensive atlas of high-confidence translated smORFs derived from primary human cells and tissues in order to provide a more complete understanding of the translated human genome.

Metagenomic Next-Generation Sequencing Reveals the Profile of Viral Infections in Kidney Transplant Recipients During the COVID-19 Pandemic

Background

To study the clinical application of metagenomic next-generation sequencing (mNGS) in the detection of viral infections in kidney transplant recipients (KTRs) during the COVID-19 pandemic.

Methods

Using mNGS technology, 50 human fluid samples of KTRs were detected, including 20 bronchoalveolar lavage fluid (BALF) samples, 21 urine samples and 9 blood samples. The detected nucleic acid sequences were compared and analyzed with the existing viral nucleic acid sequences in the database, and the virus infection spectrum of KTRs was drawn.

Results

The viral nucleic acids of 15 types of viruses were detected in 96.00% (48/50) of the samples, of which 11 types of viruses were in BALF (95.00%, 19/20), and the dominant viruses were torque teno virus (TTV) (65.00%; 13/20), cytomegalovirus (CMV) (45.00%; 9/20) and human alphaherpesvirus 1 (25.00%; 5/20). 12 viruses (95.24%, 20/21) were detected in the urine, and the dominant viruses were TTV (52.38%; 11/21), JC polyomavirus (52.38%; 11/21), BK polyomavirus (42.86%; 9/21), CMV (33.33%; 7/21) and human betaherpesvirus 6B (28.57%; 6/21). 7 viruses were detected in the blood (100.00%, 9/9), and the dominant virus was TTV (100.00%; 9/9). Four rare viruses were detected in BALF and urine, including WU polyomavirus, primate bocaparvovirus 1, simian virus 12, and volepox virus. Further analysis showed that TTV infection with high reads indicated a higher risk of acute rejection (P < 0.05).

Conclusions

mNGS detection reveals the rich virus spectrum of infected KTRs, and improves the detection rate of rare viruses. TTV may be a new biomarker for predicting rejection.

The effects of sequencing depth on the assembly of coding and noncoding transcripts in the human genome

Investigating the functions and activities of genes requires proper annotation of the transcribed units. However, transcript assembly efforts have produced a surprisingly large variation in the number of transcripts, and especially so for noncoding transcripts. This heterogeneity in assembled transcript sets might be partially explained by sequencing depth. Here, we used real and simulated short-read sequencing data as well as long-read data to systematically investigate the impact of sequencing depths on the accuracy of assembled transcripts. We assembled and analyzed transcripts from 671 human short-read data sets and four long-read data sets. At the first level, there is a positive correlation between the number of reads and the number of recovered transcripts. However, the effect of the sequencing depth varied based on cell or tissue type, the type of read and the nature and expression levels of the transcripts. The detection of coding transcripts saturated rapidly with both short and long-reads, however, there was no sign of early saturation for noncoding transcripts at any sequencing depth. Increasing long-read sequencing depth specifically benefited transcripts containing transposable elements. Finally, we show how single-cell RNA-seq can be guided by transcripts assembled from bulk long-read samples, and demonstrate that noncoding transcripts are expressed at similar levels to coding transcripts but are expressed in fewer cells. This study highlights the impact of sequencing depth on transcript assembly.

SCN1A overexpression, associated with a genomic region marked by a risk variant for a common epilepsy, raises seizure susceptibility

Mesial temporal lobe epilepsy with hippocampal sclerosis and a history of febrile seizures is associated with common variation at rs7587026, located in the promoter region of SCN1A. We sought to explore possible underlying mechanisms. SCN1A expression was analysed in hippocampal biopsy specimens of individuals with mesial temporal lobe epilepsy with hippocampal sclerosis who underwent surgical treatment, and hippocampal neuronal cell loss was quantitatively assessed using immunohistochemistry. In healthy individuals, hippocampal volume was measured using MRI. Analyses were performed stratified by rs7587026 type. To study the functional consequences of increased SCN1A expression, we generated, using transposon-mediated bacterial artificial chromosome transgenesis, a zebrafish line expressing exogenous scn1a, and performed EEG analysis on larval optic tecta at 4 day post-fertilization. Finally, we used an in vitro promoter analysis to study whether the genetic motif containing rs7587026 influences promoter activity. Hippocampal SCN1A expression differed by rs7587026 genotype (Kruskal-Wallis test P = 0.004). Individuals homozygous for the minor allele showed significantly increased expression compared to those homozygous for the major allele (Dunn's test P = 0.003), and to heterozygotes (Dunn's test P = 0.035). No statistically significant differences in hippocampal neuronal cell loss were observed between the three genotypes. Among 597 healthy participants, individuals homozygous for the minor allele at rs7587026 displayed significantly reduced mean hippocampal volume compared to major allele homozygotes (Cohen's D = - 0.28, P = 0.02), and to heterozygotes (Cohen's D = - 0.36, P = 0.009). Compared to wild type, scn1lab-overexpressing zebrafish larvae exhibited more frequent spontaneous seizures [one-way ANOVA F(4,54) = 6.95 (P < 0.001)]. The number of EEG discharges correlated with the level of scn1lab overexpression [one-way ANOVA F(4,15) = 10.75 (P < 0.001]. Finally, we showed that a 50 bp promoter motif containing rs7587026 exerts a strong regulatory role on SCN1A expression, though we could not directly link this to rs7587026 itself. Our results develop the mechanistic link between rs7587026 and mesial temporal lobe epilepsy with hippocampal sclerosis and a history of febrile seizures. Furthermore, we propose that quantitative precision may be important when increasing SCN1A expression in current strategies aiming to treat seizures in conditions involving SCN1A haploinsufficiency, such as Dravet syndrome.

A framework to score the effects of structural variants in health and disease

While technological advances improved the identification of structural variants (SVs) in the human genome, their interpretation remains challenging. Several methods utilize individual mechanistic principles like the deletion of coding sequence or 3D genome architecture disruptions. However, a comprehensive tool using the broad spectrum of available annotations is missing. Here, we describe CADD-SV, a method to retrieve and integrate a wide set of annotations to predict the effects of SVs. Previously, supervised learning approaches were limited due to a small number and biased set of annotated pathogenic or benign SVs. We overcome this problem by using a surrogate training-objective, the Combined Annotation Dependent Depletion (CADD) of functional variants. We use human and chimpanzee derived SVs as proxy-neutral and contrast them with matched simulated variants as proxy-deleterious, an approach that has proven powerful for short sequence variants. Our tool computes summary statistics over diverse variant annotations and uses random forest models to prioritize deleterious structural variants. The resulting CADD-SV scores correlate with known pathogenic and rare population variants. We further show that we can prioritize somatic cancer variants as well as noncoding variants known to affect gene expression. We provide a website and offline-scoring tool for easy application of CADD-SV.

Preclinical Research in McArdle Disease: A Review of Research Models and Therapeutic Strategies

McArdle disease is an autosomal recessive disorder of muscle glycogen metabolism caused by pathogenic mutations in the PYGM gene, which encodes the skeletal muscle-specific isoform of glycogen phosphorylase. Clinical symptoms are mainly characterized by transient acute “crises” of early fatigue, myalgia and contractures, which can be accompanied by rhabdomyolysis. Owing to the difficulty of performing mechanistic studies in patients that often rely on invasive techniques, preclinical models have been used for decades, thereby contributing to gain insight into the pathophysiology and pathobiology of human diseases. In the present work, we describe the existing in vitro and in vivo preclinical models for McArdle disease and review the insights these models have provided. In addition, despite presenting some differences with the typical patient’s phenotype, these models allow for a deep study of the different features of the disease while representing a necessary preclinical step to assess the efficacy and safety of possible treatments before they are tested in patients.

Integration of high-resolution promoter profiling assays reveals novel, cell type-specific transcription start sites across 115 human cell and tissue types

Accurate transcription start site (TSS) annotations are essential for understanding transcriptional regulation and its role in human disease. Gene collections such as GENCODE contain annotations for tens of thousands of TSSs, but not all of these annotations are experimentally validated nor do they contain information on cell type–specific usage. Therefore, we sought to generate a collection of experimentally validated TSSs by integrating RNA Annotation and Mapping of Promoters for the Analysis of Gene Expression (RAMPAGE) data from 115 cell and tissue types, which resulted in a collection of approximately 50 thousand representative RAMPAGE peaks. These peaks are primarily proximal to GENCODE-annotated TSSs and are concordant with other transcription assays. Because RAMPAGE uses paired-end reads, we were then able to connect peaks to transcripts by analyzing the genomic positions of the 3′ ends of read mates. Using this paired-end information, we classified the vast majority (37 thousand) of our RAMPAGE peaks as verified TSSs, updating TSS annotations for 20% of GENCODE genes. We also found that these updated TSS annotations are supported by epigenomic and other transcriptomic data sets. To show the utility of this RAMPAGE rPeak collection, we intersected it with the NHGRI/EBI genome-wide association study (GWAS) catalog and identified new candidate GWAS genes. Overall, our work shows the importance of integrating experimental data to further refine TSS annotations and provides a valuable resource for the biological community.

Muscle Regeneration and RNA: New Perspectives for Ancient Molecules

The ability of the ribonucleic acid (RNA) to self-replicate, combined with a unique cocktail of chemical properties, suggested the existence of an RNA world at the origin of life. Nowadays, this hypothesis is supported by innovative high-throughput and biochemical approaches, which definitively revealed the essential contribution of RNA-mediated mechanisms to the regulation of fundamental processes of life. With the recent development of SARS-CoV-2 mRNA-based vaccines, the potential of RNA as a therapeutic tool has received public attention. Due to its intrinsic single-stranded nature and the ease with which it is synthesized in vitro, RNA indeed represents the most suitable tool for the development of drugs encompassing every type of human pathology. The maximum effectiveness and biochemical versatility is achieved in the guise of non-coding RNAs (ncRNAs), which are emerging as multifaceted regulators of tissue specification and homeostasis. Here, we report examples of coding and ncRNAs involved in muscle regeneration and discuss their potential as therapeutic tools. Small ncRNAs, such as miRNA and siRNA, have been successfully applied in the treatment of several diseases. The use of longer molecules, such as lncRNA and circRNA, is less advanced. However, based on the peculiar properties discussed below, they represent an innovative pool of RNA biomarkers and possible targets of clinical value.

aniFOUND: analysing the associated proteome and genomic landscape of the repaired nascent non-replicative chromatin

Specific capture of chromatin fractions with distinct and well-defined features has emerged as both challenging and a key strategy towards a comprehensive understanding of genome biology. In this context, we developed aniFOUND (accelerated native isolation of factors on unscheduled nascent DNA), an antibody-free method, which can label, capture, map and characterise nascent chromatin fragments that are synthesized in response to specific cues outside S-phase. We used the 'unscheduled' DNA synthesis (UDS) that takes place during the repair of UV-induced DNA lesions and coupled the captured chromatin to high-throughput analytical technologies. By mass-spectrometry we identified several factors with no previously known role in UVC-DNA damage response (DDR) as well as known DDR proteins. We experimentally validated the repair-dependent recruitment of the chromatin remodeller RSF1 and the cohesin-loader NIPBL at sites of UVC-induced photolesions. Developing aniFOUND-seq, a protocol for mapping UDS activity with high resolution, allowed us to monitor the landscape of UVC repair-synthesis events genome wide. We further resolved repair efficacy of the rather unexplored repeated genome, in particular rDNA and telomeres. In summary, aniFOUND delineates the proteome composition and genomic landscape of chromatin loci with specific features by integrating state-of-the-art 'omics' technologies to promote a comprehensive view of their function.

Intergenic RNA mainly derives from nascent transcripts of known genes

BACKGROUND

Eukaryotic genomes undergo pervasive transcription, leading to the production of many types of stable and unstable RNAs. Transcription is not restricted to regions with annotated gene features but includes almost any genomic context. Currently, the source and function of most RNAs originating from intergenic regions in the human genome remain unclear.

RESULTS

We hypothesize that many intergenic RNAs can be ascribed to the presence of as-yet unannotated genes or the "fuzzy" transcription of known genes that extends beyond the annotated boundaries. To elucidate the contributions of these two sources, we assemble a dataset of more than 2.5 billion publicly available RNA-seq reads across 5 human cell lines and multiple cellular compartments to annotate transcriptional units in the human genome. About 80% of transcripts from unannotated intergenic regions can be attributed to the fuzzy transcription of existing genes; the remaining transcripts originate mainly from putative long non-coding RNA loci that are rarely spliced. We validate the transcriptional activity of these intergenic RNAs using independent measurements, including transcriptional start sites, chromatin signatures, and genomic occupancies of RNA polymerase II in various phosphorylation states. We also analyze the nuclear localization and sensitivities of intergenic transcripts to nucleases to illustrate that they tend to be rapidly degraded either on-chromatin by XRN2 or off-chromatin by the exosome.

CONCLUSIONS

We provide a curated atlas of intergenic RNAs that distinguishes between alternative processing of well-annotated genes from independent transcriptional units based on the combined analysis of chromatin signatures, nuclear RNA localization, and degradation pathways.

Muscle Glycogen Phosphorylase and Its Functional Partners in Health and Disease

Glycogen phosphorylase (PG) is a key enzyme taking part in the first step of glycogenolysis. Muscle glycogen phosphorylase (PYGM) differs from other PG isoforms in expression pattern and biochemical properties. The main role of PYGM is providing sufficient energy for muscle contraction. However, it is expressed in tissues other than muscle, such as the brain, lymphoid tissues, and blood. PYGM is important not only in glycogen metabolism, but also in such diverse processes as the insulin and glucagon signaling pathway, insulin resistance, necroptosis, immune response, and phototransduction. PYGM is implicated in several pathological states, such as muscle glycogen phosphorylase deficiency (McArdle disease), schizophrenia, and cancer. Here we attempt to analyze the available data regarding the protein partners of PYGM to shed light on its possible interactions and functions. We also underline the potential for zebrafish to become a convenient and applicable model to study PYGM functions, especially because of its unique features that can complement data obtained from other approaches.

FANTOM enters 20th year: expansion of transcriptomic atlases and functional annotation of non-coding RNAs

The Functional ANnoTation Of the Mammalian genome (FANTOM) Consortium has continued to provide extensive resources in the pursuit of understanding the transcriptome, and transcriptional regulation, of mammalian genomes for the last 20 years. To share these resources with the research community, the FANTOM web-interfaces and databases are being regularly updated, enhanced and expanded with new data types. In recent years, the FANTOM Consortium's efforts have been mainly focused on creating new non-coding RNA datasets and resources. The existing FANTOM5 human and mouse miRNA atlas was supplemented with rat, dog, and chicken datasets. The sixth (latest) edition of the FANTOM project was launched to assess the function of human long non-coding RNAs (lncRNAs). From its creation until 2020, FANTOM6 has contributed to the research community a large dataset generated from the knock-down of 285 lncRNAs in human dermal fibroblasts; this is followed with extensive expression profiling and cellular phenotyping. Other updates to the FANTOM resource includes the reprocessing of the miRNA and promoter atlases of human, mouse and chicken with the latest reference genome assemblies. To facilitate the use and accessibility of all above resources we further enhanced FANTOM data viewers and web interfaces. The updated FANTOM web resource is publicly available at https://fantom.gsc.riken.jp/.

Epigenetic regulation of Wnt7b expression by the cis-acting long noncoding RNA Lnc-Rewind in muscle stem cells

Skeletal muscle possesses an outstanding capacity to regenerate upon injury due to the adult muscle stem cell (MuSC) activity. This ability requires the proper balance between MuSC expansion and differentiation, which is critical for muscle homeostasis and contributes, if deregulated, to muscle diseases. Here, we functionally characterize a novel chromatin-associated long noncoding RNA (lncRNA), Lnc-Rewind, which is expressed in murine MuSCs and conserved in human. We find that, in mouse, Lnc-Rewind acts as an epigenetic regulator of MuSC proliferation and expansion by influencing the expression of skeletal muscle genes and several components of the WNT (Wingless-INT) signalling pathway. Among them, we identified the nearby Wnt7b gene as a direct Lnc-Rewind target. We show that Lnc-Rewind interacts with the G9a histone lysine methyltransferase and mediates the in cis repression of Wnt7b by H3K9me2 deposition. Overall, these findings provide novel insights into the epigenetic regulation of adult muscle stem cells fate by lncRNAs.

Assessment of transcriptional importance of cell line-specific features based on GTRD and FANTOM5 data

Creating a complete picture of the regulation of transcription seems to be an urgent task of modern biology. Regulation of transcription is a complex process carried out by transcription factors (TFs) and auxiliary proteins. Over the past decade, ChIP-Seq has become the most common experimental technology studying genome-wide interactions between TFs and DNA. We assessed the transcriptional significance of cell line-specific features using regression analysis of ChIP-Seq datasets from the GTRD database and transcriptional start site (TSS) activities from the FANTOM5 expression atlas. For this purpose, we initially generated a large number of features that were defined as the presence or absence of TFs in different promoter regions around TSSs. Using feature selection and regression analysis, we identified sets of the most important TFs that affect expression activity of TSSs in human cell lines such as HepG2, K562 and HEK293. We demonstrated that some TFs can be classified as repressors and activators depending on their location relative to TSS.

Gene expression and in situ protein profiling of candidate SARS-CoV-2 receptors in human airway epithelial cells and lung tissue

In December 2019, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged, causing the coronavirus disease 2019 (COVID-19) pandemic. SARS-CoV, the agent responsible for the 2003 SARS outbreak, utilises angiotensin-converting enzyme 2 (ACE2) and transmembrane serine protease 2 (TMPRSS2) host molecules for viral entry. ACE2 and TMPRSS2 have recently been implicated in SARS-CoV-2 viral infection. Additional host molecules including ADAM17, cathepsin L, CD147 and GRP78 may also function as receptors for SARS-CoV-2.To determine the expression and in situ localisation of candidate SARS-CoV-2 receptors in the respiratory mucosa, we analysed gene expression datasets from airway epithelial cells of 515 healthy subjects, gene promoter activity analysis using the FANTOM5 dataset containing 120 distinct sample types, single cell RNA sequencing (scRNAseq) of 10 healthy subjects, proteomic datasets, immunoblots on multiple airway epithelial cell types, and immunohistochemistry on 98 human lung samples.We demonstrate absent to low ACE2 promoter activity in a variety of lung epithelial cell samples and low ACE2 gene expression in both microarray and scRNAseq datasets of epithelial cell populations. Consistent with gene expression, rare ACE2 protein expression was observed in the airway epithelium and alveoli of human lung, confirmed with proteomics. We present confirmatory evidence for the presence of TMPRSS2, CD147 and GRP78 protein in vitro in airway epithelial cells and confirm broad in situ protein expression of CD147 and GRP78 in the respiratory mucosa.Collectively, our data suggest the presence of a mechanism dynamically regulating ACE2 expression in human lung, perhaps in periods of SARS-CoV-2 infection, and also suggest that alternative receptors for SARS-CoV-2 exist to facilitate initial host cell infection.

Gene expression and in situ protein profiling of candidate SARS-CoV-2 receptors in human airway epithelial cells and lung tissue

In December 2019, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged, causing the coronavirus disease 2019 (COVID-19) pandemic. SARS-CoV, the agent responsible for the 2003 SARS outbreak, utilises angiotensin-converting enzyme 2 (ACE2) and transmembrane serine protease 2 (TMPRSS2) host molecules for viral entry. ACE2 and TMPRSS2 have recently been implicated in SARS-CoV-2 viral infection. Additional host molecules including ADAM17, cathepsin L, CD147 and GRP78 may also function as receptors for SARS-CoV-2.To determine the expression and in situ localisation of candidate SARS-CoV-2 receptors in the respiratory mucosa, we analysed gene expression datasets from airway epithelial cells of 515 healthy subjects, gene promoter activity analysis using the FANTOM5 dataset containing 120 distinct sample types, single cell RNA sequencing (scRNAseq) of 10 healthy subjects, proteomic datasets, immunoblots on multiple airway epithelial cell types, and immunohistochemistry on 98 human lung samples.We demonstrate absent to low ACE2 promoter activity in a variety of lung epithelial cell samples and low ACE2 gene expression in both microarray and scRNAseq datasets of epithelial cell populations. Consistent with gene expression, rare ACE2 protein expression was observed in the airway epithelium and alveoli of human lung, confirmed with proteomics. We present confirmatory evidence for the presence of TMPRSS2, CD147 and GRP78 protein in vitro in airway epithelial cells and confirm broad in situ protein expression of CD147 and GRP78 in the respiratory mucosa.Collectively, our data suggest the presence of a mechanism dynamically regulating ACE2 expression in human lung, perhaps in periods of SARS-CoV-2 infection, and also suggest that alternative receptors for SARS-CoV-2 exist to facilitate initial host cell infection.

Expression profile of immune checkpoint genes and their roles in predicting immunotherapy response

Immune checkpoint genes (ICGs) play critical roles in circumventing self-reactivity and represent a novel target to develop treatments for cancers. However, a comprehensive analysis for the expression profile of ICGs at a pan-cancer level and their correlation with patient response to immune checkpoint blockade (ICB) based therapy is still lacking. In this study, we defined three expression patterns of ICGs using a comprehensive survey of RNA-seq data of tumor and immune cells from the functional annotation of the mammalian genome (FANTOM5) project. The correlation between the expression patterns of ICGs and patients survival and response to ICB therapy was investigated. The expression patterns of ICGs were robust across cancers, and upregulation of ICGs was positively correlated with high lymphocyte infiltration and good prognosis. Furthermore, we built a model (ICGe) to predict the response of patients to ICB therapy using five features of ICG expression. A validation scenario of six independent datasets containing data of 261 patients with CTLA-4 and PD-1 blockade immunotherapies demonstrated that ICGe achieved area under the curves of 0.64-0.82 and showed a robust performance and outperformed other mRNA-based predictors. In conclusion, this work revealed expression patterns of ICGs and underlying correlations between ICGs and response to ICB, which helps to understand the mechanisms of ICGs in ICB signal pathways and other anticancer treatments.

Comparative transcriptomics of primary cells in vertebrates

Gene expression profiles in homologous tissues have been observed to be different between species, which may be due to differences between species in the gene expression program in each cell type, but may also reflect differences in cell type composition of each tissue in different species. Here, we compare expression profiles in matching primary cells in human, mouse, rat, dog, and chicken using Cap Analysis Gene Expression (CAGE) and short RNA (sRNA) sequencing data from FANTOM5. While we find that expression profiles of orthologous genes in different species are highly correlated across cell types, in each cell type many genes were differentially expressed between species. Expression of genes with products involved in transcription, RNA processing, and transcriptional regulation was more likely to be conserved, while expression of genes encoding proteins involved in intercellular communication was more likely to have diverged during evolution. Conservation of expression correlated positively with the evolutionary age of genes, suggesting that divergence in expression levels of genes critical for cell function was restricted during evolution. Motif activity analysis showed that both promoters and enhancers are activated by the same transcription factors in different species. An analysis of expression levels of mature miRNAs and of primary miRNAs identified by CAGE revealed that evolutionary old miRNAs are more likely to have conserved expression patterns than young miRNAs. We conclude that key aspects of the regulatory network are conserved, while differential expression of genes involved in cell-to-cell communication may contribute greatly to phenotypic differences between species.

Bacteria pathogens drive host colonic epithelial cell promoter hypermethylation of tumor suppressor genes in colorectal cancer

BACKGROUND

Altered microbiome composition and aberrant promoter hypermethylation of tumor suppressor genes (TSGs) are two important hallmarks of colorectal cancer (CRC). Here we performed concurrent 16S rRNA gene sequencing and methyl-CpG binding domain-based capture sequencing in 33 tissue biopsies (5 normal colonic mucosa tissues, 4 pairs of adenoma and adenoma-adjacent tissues, and 10 pairs of CRC and CRC-adjacent tissues) to identify significant associations between TSG promoter hypermethylation and CRC-associated bacteria, followed by functional validation of the methylation-associated bacteria.

RESULTS

Fusobacterium nucleatum and Hungatella hathewayi were identified as the top two methylation-regulating bacteria. Targeted analysis on bona fide TSGs revealed that H. hathewayi and Streptococcus spp. significantly correlated with CDX2 and MLH1 promoter hypermethylation, respectively. Mechanistic validation with cell-line and animal models revealed that F. nucleatum and H. hathewayi upregulated DNA methyltransferase. H. hathewayi inoculation also promoted colonic epithelial cell proliferation in germ-free and conventional mice.

CONCLUSION

Our integrative analysis revealed previously unknown epigenetic regulation of TSGs in host cells through inducing DNA methyltransferase by F. nucleatum and H. hathewayi, and established the latter as CRC-promoting bacteria. Video abstract.

Evidence That STK19 Is Not an NRAS-dependent Melanoma Driver

STK19 was proposed to be a cancer driver, and recent work by Yin et al. (2019) in Cell suggested that the frequently recurring STK19 D89N substitution represents a gain-of-function change, allowing increased phosphorylation of NRAS to enhance melanocyte transformation. Here we show that the STK19 gene has been incorrectly annotated, and that the expressed protein is 110 amino acids shorter than indicated by current databases. The "cancer driving" STK19 D89N substitution is thus outside the coding region. We also fail to detect evidence of the mutation affecting STK19 expression; instead, it is a UV signature mutation, found in the promoter of other genes as well. Furthermore, STK19 is exclusively nuclear and chromatin-associated, while no evidence for it being a kinase was found. The data in this Matters Arising article raise fundamental questions about the recently proposed role for STK19 in melanoma progression via a function as an NRAS kinase, suggested by Yin et al. (2019) in Cell. See also the response by Yin et al. (2020), published in this issue.

DPre: computational identification of differentiation bias and genes underlying cell type conversions

SUMMARY

Cells are generally resistant to cell type conversions, but can be converted by the application of growth factors, chemical inhibitors and ectopic expression of genes. However, it remains difficult to accurately identify the destination cell type or differentiation bias when these techniques are used to alter cell type. Consequently, there is demand for computational techniques that can help researchers understand both the cell type and differentiation bias. While advanced tools identifying cell types exist for single cell data and the deconvolution of mixed cell populations, the problem of exploring partially differentiated cells of indeterminate transcriptional identity has not been addressed. To fill this gap, we developed driver-predictor, which relies on scoring per gene transcriptional similarity between RNA-Seq datasets to reveal directional bias of differentiation. By comparing against large cell type transcriptome libraries or a desired target expression profile, the tool enables the user to visualize both the changes in transcriptional identity as well as the genes accounting for the cell type changes. This software will be a powerful tool for researchers to explore in vitro experiments that involve cell type conversions.

AVAILABILITY AND IMPLEMENTATION

Source code is open source under the MIT license and is freely available on https://github.com/LoaloaF/DPre.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A human minisatellite hosts an alternative transcription start site for NPRL3 driving its expression in a repeat number-dependent manner

Minisatellites, also called variable number of tandem repeats (VNTRs), are a class of repetitive elements that may affect gene expression at multiple levels and have been correlated to disease. Their identification and role as expression quantitative trait loci (eQTL) have been limited by their absence in comparative genomic hybridization and single nucleotide polymorphisms arrays. By taking advantage of cap analysis of gene expression (CAGE), we describe a new example of a minisatellite hosting a transcription start site (TSS) which expression is dependent on the repeat number. It is located in the third intron of the gene nitrogen permease regulator like protein 3 (NPRL3). NPRL3 is a component of the GAP activity toward rags 1 protein complex that inhibits mammalian target of rapamycin complex 1 (mTORC1) activity and it is found mutated in familial focal cortical dysplasia and familial focal epilepsy. CAGE tags represent an alternative TSS identifying TAGNPRL3 messenger RNAs (mRNAs). TAGNPRL3 is expressed in red blood cells both at mRNA and protein levels, it interacts with its protein partner NPRL2 and its overexpression inhibits cell proliferation. This study provides an example of a minisatellite that is both a TSS and an eQTL as well as identifies a new VNTR that may modify mTORC1 activity.

TRPV6 as A Target for Cancer Therapy

Two decades ago a class of ion channels, hitherto unsuspected, was discovered. In mammals these Transient Receptor Potential channels (TRPs) have not only expanded in number (to 26 functional channels) but also expanded the view of our interface with the physical and chemical environment. Some are heat and cold sensors while others monitor endogenous and/or exogenous chemical signals. Some TRP channels monitor osmotic potential, and others measure cell movement, stretching, and fluid flow. Many TRP channels are major players in nociception and integration of pain signals. One member of the vanilloid sub-family of channels is TRPV6. This channel is highly selective for divalent cations, particularly calcium, and plays a part in general whole-body calcium homeostasis, capturing calcium in the gut from the diet. TRPV6 can be greatly elevated in a number of cancers deriving from epithelia and considerable study has been made of its role in the cancer phenotype where calcium control is dysfunctional. This review compiles and updates recent published work on TRPV6 as a promising drug target in a number of cancers including those afflicting breast, ovarian, prostate and pancreatic tissues.

Atlas of quantitative single-base-resolution N6-methyl-adenine methylomes

Various methyltransferases and demethylases catalyse methylation and demethylation of N⁶-methyladenosine (m6A) and N⁶,2'-O-dimethyladenosine (m6Am) but precise methylomes uniquely mediated by each methyltransferase/demethylase are still lacking. Here, we develop m6A-Crosslinking-Exonuclease-sequencing (m6ACE-seq) to map transcriptome-wide m6A and m6Am at quantitative single-base-resolution. This allows for the generation of a comprehensive atlas of distinct methylomes uniquely mediated by every individual known methyltransferase or demethylase. Our atlas reveals METTL16 to indirectly impact manifold methylation targets beyond its consensus target motif and highlights the importance of precision in mapping PCIF1-dependent m6Am. Rather than reverse RNA methylation, we find that both ALKBH5 and FTO instead maintain their regulated sites in an unmethylated steady-state. In FTO's absence, anomalous m6Am disrupts snRNA interaction with nuclear export machinery, potentially causing aberrant pre-mRNA splicing events.

Antibody cross-reactivity accounts for widespread appearance of m1A in 5'UTRs

N1-methyladenosine (m1A) was proposed to be a highly prevalent modification in mRNA 5'UTRs based on mapping studies using an m1A-binding antibody. We developed a bioinformatic approach to discover m1A and other modifications in mRNA throughout the transcriptome by analyzing preexisting ultra-deep RNA-Seq data for modification-induced misincorporations. Using this approach, we detected appreciable levels of m1A only in one mRNA: the mitochondrial MT-ND5 transcript. As an alternative approach, we also developed an antibody-based m1A-mapping approach to detect m1A at single-nucleotide resolution, and confirmed that the commonly used m1A antibody maps sites to the transcription-start site in mRNA 5'UTRs. However, further analysis revealed that these were false-positives caused by binding of the antibody to the m7G-cap. A different m1A antibody that lacks cap-binding cross-reactivity does not show enriched binding in 5'UTRs. These results demonstrate that high-stoichiometry m1A sites are exceedingly rare in mRNAs and that previous mappings of m1A to 5'UTRs were the result of antibody cross-reactivity to the 5' cap.

Identification of the m6Am Methyltransferase PCIF1 Reveals the Location and Functions of m6Am in the Transcriptome

mRNAs are regulated by nucleotide modifications that influence their cellular fate. Two of the most abundant modified nucleotides are N⁶-methyladenosine (m⁶A), found within mRNAs, and N⁶,2'-O-dimethyladenosine (m⁶Am), which is found at the first transcribed nucleotide. Distinguishing these modifications in mapping studies has been difficult. Here, we identify and biochemically characterize PCIF1, the methyltransferase that generates m⁶Am. We find that PCIF1 binds and is dependent on the m⁷G cap. By depleting PCIF1, we generated transcriptome-wide maps that distinguish m⁶Am and m⁶A. We find that m⁶A and m⁶Am misannotations arise from mRNA isoforms with alternative transcription start sites (TSSs). These isoforms contain m⁶Am that maps to "internal" sites, increasing the likelihood of misannotation. We find that depleting PCIF1 does not substantially affect mRNA translation but is associated with reduced stability of a subset of m⁶Am-annotated mRNAs. The discovery of PCIF1 and our accurate mapping technique will facilitate future studies to characterize m⁶Am's function.