Long-read sequencing transcriptome quantification with lr-kallisto.

RNA abundance quantification has become routine and affordable thanks to high-throughput "short-read" technologies that provide accurate molecule counts at the gene level. Similarly accurate and affordable quantification of definitive full-length, transcript isoforms has remained a stubborn challenge, despite its obvious biological significance across a wide range of problems. "Long-read" sequencing platforms now produce data-types that can, in principle, drive routine definitive isoform quantification. However some particulars of contemporary long-read datatypes, together with isoform complexity and genetic variation, present bioinformatic challenges. We show here, using ONT data, that fast and accurate quantification of long-read data is possible and that it is improved by exome capture. To perform quantifications we developed lr-kallisto, which adapts the kallisto bulk and single-cell RNA-seq quantification methods for long-read technologies.

Engineered FSHD mutations results in D4Z4 heterochromatin disruption and feedforward DUX4 network activation.

Facioscapulohumeral dystrophy (FSHD) is linked to contraction of D4Z4 repeats on chromosome 4q with SMCHD1 mutations acting as a disease modifier. D4Z4 heterochromatin disruption and abnormal upregulation of the transcription factor DUX4, encoded in the D4Z4 repeat, are the hallmarks of FSHD. However, defining the precise effect of D4Z4 contraction has been difficult because D4Z4 repeats are primate-specific and DUX4 expression is very rare in highly heterogeneous patient myocytes. We generated isogenic mutant cell lines harboring D4Z4 and/or SMCHD1 mutations in a healthy human skeletal myoblast line. We found that the mutations affect D4Z4 heterochromatin differently, and that SMCHD1 mutation or disruption of DNA methylation stabilizes otherwise variegated DUX4 target activation in D4Z4 contraction mutant cells, demonstrating the critical role of modifiers. Our study revealed amplification of the DUX4 signal through downstream targets, H3.X/Y and LEUTX. Our results provide important insights into how rare DUX4 expression leads to FSHD pathogenesis.

The ENCODE4 long-read RNA-seq collection reveals distinct classes of transcript structure diversity.

The majority of mammalian genes encode multiple transcript isoforms that result from differential promoter use, changes in exonic splicing, and alternative 3' end choice. Detecting and quantifying transcript isoforms across tissues, cell types, and species has been extremely challenging because transcripts are much longer than the short reads normally used for RNA-seq. By contrast, long-read RNA-seq (LR-RNA-seq) gives the complete structure of most transcripts. We sequenced 264 LR-RNA-seq PacBio libraries totaling over 1 billion circular consensus reads (CCS) for 81 unique human and mouse samples. We detect at least one full-length transcript from 87.7% of annotated human protein coding genes and a total of 200,000 full-length transcripts, 40% of which have novel exon junction chains. To capture and compute on the three sources of transcript structure diversity, we introduce a gene and transcript annotation framework that uses triplets representing the transcript start site, exon junction chain, and transcript end site of each transcript. Using triplets in a simplex representation demonstrates how promoter selection, splice pattern, and 3' processing are deployed across human tissues, with nearly half of multi-transcript protein coding genes showing a clear bias toward one of the three diversity mechanisms. Evaluated across samples, the predominantly expressed transcript changes for 74% of protein coding genes. In evolution, the human and mouse transcriptomes are globally similar in types of transcript structure diversity, yet among individual orthologous gene pairs, more than half (57.8%) show substantial differences in mechanism of diversification in matching tissues. This initial large-scale survey of human and mouse long-read transcriptomes provides a foundation for further analyses of alternative transcript usage, and is complemented by short-read and microRNA data on the same samples and by epigenome data elsewhere in the ENCODE4 collection.

A Trem2R47H mouse model without cryptic splicing drives age- and disease-dependent tissue damage and synaptic loss in response to plaques.

BACKGROUND: The TREM2 R47H variant is one of the strongest genetic risk factors for late-onset Alzheimer's Disease (AD). Unfortunately, many current Trem2 R47H mouse models are associated with cryptic mRNA splicing of the mutant allele that produces a confounding reduction in protein product. To overcome this issue, we developed the Trem2R47H NSS (Normal Splice Site) mouse model in which the Trem2 allele is expressed at a similar level to the wild-type Trem2 allele without evidence of cryptic splicing products. METHODS: Trem2R47H NSS mice were treated with the demyelinating agent cuprizone, or crossed with the 5xFAD mouse model of amyloidosis, to explore the impact of the TREM2 R47H variant on inflammatory responses to demyelination, plaque development, and the brain's response to plaques. RESULTS: Trem2R47H NSS mice display an appropriate inflammatory response to cuprizone challenge, and do not recapitulate the null allele in terms of impeded inflammatory responses to demyelination. Utilizing the 5xFAD mouse model, we report age- and disease-dependent changes in Trem2R47H NSS mice in response to development of AD-like pathology. At an early (4-month-old) disease stage, hemizygous 5xFAD/homozygous Trem2R47H NSS (5xFAD/Trem2R47H NSS) mice have reduced size and number of microglia that display impaired interaction with plaques compared to microglia in age-matched 5xFAD hemizygous controls. This is associated with a suppressed inflammatory response but increased dystrophic neurites and axonal damage as measured by plasma neurofilament light chain (NfL) level. Homozygosity for Trem2R47H NSS suppressed LTP deficits and loss of presynaptic puncta caused by the 5xFAD transgene array in 4-month-old mice. At a more advanced (12-month-old) disease stage 5xFAD/Trem2R47H NSS mice no longer display impaired plaque-microglia interaction or suppressed inflammatory gene expression, although NfL levels remain elevated, and a unique interferon-related gene expression signature is seen. Twelve-month old Trem2R47H NSS mice also display LTP deficits and postsynaptic loss. CONCLUSIONS: The Trem2R47H NSS mouse is a valuable model that can be used to investigate age-dependent effects of the AD-risk R47H mutation on TREM2 and microglial function including its effects on plaque development, microglial-plaque interaction, production of a unique interferon signature and associated tissue damage.

Mutually Orthogonal Bioconjugation of Vinyl Nucleosides for RNA Metabolic Labeling.

We report a strategy for the orthogonal conjugation of the vinyl nucleosides, 5-vinyluridine (5-VU) and 2-vinyladenosine (2-VA), via selective reactivity with maleimide and tris(2-carboxyethyl)phosphine (TCEP), respectively. The orthogonality was investigated using density functional theory (DFT) and confirmed by reactions with vinyl nucleosides. Further, these chemistries were used to modify RNA for fluorescent cell imaging. These reactions allow for the expanded use of RNA metabolic labeling to study nascent RNA expression within different RNA populations.

A Bump-Hole Strategy for Increased Stringency of Cell-Specific Metabolic Labeling of RNA.

Profiling RNA expression in a cell-specific manner continues to be a grand challenge in biochemical research. Bioorthogonal nucleosides can be utilized to track RNA expression; however, these methods currently have limitations due to background and incorporation of analogs into undesired cells. Herein, we design and demonstrate that uracil phosphoribosyltransferase can be engineered to match 5-vinyluracil for cell-specific metabolic labeling of RNA with exceptional specificity and stringency.

Melanin Concentrating Hormone Signaling Deficits in Schizophrenia: Association With Memory and Social Impairments and Abnormal Sensorimotor Gating.

BACKGROUND: Evidence from anatomical, pharmacological, and genetic studies supports a role for the neuropeptide melanin concentrating hormone system in modulating emotional and cognitive functions. Genome-wide association studies revealed a potential association between the melanin concentrating hormone receptor (MCHR1) gene locus and schizophrenia, and the largest genome-wide association study conducted to date shows a credible genome-wide association. METHODS: We analyzed MCHR1 and pro-melanin concentrating hormone RNA-Seq expression in the prefrontal cortex in schizophrenia patients and healthy controls. Disruptions in the melanin concentrating hormone system were modeled in the mouse brain by germline deletion of MCHR1 and by conditional ablation of melanin concentrating hormone expressing neurons using a Cre-inducible diphtheria toxin system. RESULTS: MCHR1 expression is decreased in the prefrontal cortex of schizophrenia samples (false discovery rate (FDR) P < .05, CommonMind and PsychEncode combined datasets, n = 901) while pro-melanin concentrating hormone is below the detection threshold. MCHR1 expression decreased with aging (P = 6.6E-57) in human dorsolateral prefrontal cortex. The deletion of MCHR1 was found to lead to behavioral abnormalities mimicking schizophrenia-like phenotypes: hyperactivity, increased stereotypic and repetitive behavior, social impairment, impaired sensorimotor gating, and disrupted cognitive functions. Conditional ablation of pro-melanin concentrating hormone neurons increased repetitive behavior and produced a deficit in sensorimotor gating. CONCLUSIONS: Our study indicates that early disruption of the melanin concentrating hormone system interferes with neurodevelopmental processes, which may contribute to the pathogenesis of schizophrenia. Further neurobiological research on the developmental timing and circuits that are affected by melanin concentrating hormone may lead to a therapeutic target for early prevention of schizophrenia.