Integrative profiling of Epstein–Barr virus transcriptome using a multiplatform approach

Background Epstein–Barr virus (EBV) is an important human pathogenic gammaherpesvirus with carcinogenic potential. The EBV transcriptome has previously been analyzed using both Illumina-based short read-sequencing and Pacific Biosciences RS II-based long-read sequencing technologies. Since the various sequencing methods have distinct strengths and limitations, the use of multiplatform approaches have proven to be valuable. The aim of this study is to provide a more complete picture on the transcriptomic architecture of EBV. Methods In this work, we apply the Oxford Nanopore Technologies MinION (long-read sequencing) platform for the generation of novel transcriptomic data, and integrate these with other’s data generated by another LRS approach, Pacific BioSciences RSII sequencing and Illumina CAGE-Seq and Poly(A)-Seq approaches. Both amplified and non-amplified cDNA sequencings were applied for the generation of sequencing reads, including both oligo-d(T) and random oligonucleotide-primed reverse transcription. EBV transcripts are identified and annotated using the LoRTIA software suite developed in our laboratory. Results This study detected novel genes embedded into longer host genes containing 5′-truncated in-frame open reading frames, which potentially encode N-terminally truncated proteins. We also detected a number of novel non-coding RNAs and transcript length isoforms encoded by the same genes but differing in their start and/or end sites. This study also reports the discovery of novel splice isoforms, many of which may represent altered coding potential, and of novel replication-origin-associated transcripts. Additionally, novel mono- and multigenic transcripts were identified. An intricate meshwork of transcriptional overlaps was revealed. Conclusions An integrative approach applying multi-technique sequencing technologies is suitable for reliable identification of complex transcriptomes because each techniques has different advantages and limitations, and the they can be used for the validation of the results obtained by a particular approach.

Estradiol Signaling at the Heart of Folliculogenesis: Its Potential Deregulation in Human Ovarian Pathologies

Estradiol (E2) is a major hormone controlling women fertility, in particular folliculogenesis. This steroid, which is locally produced by granulosa cells (GC) within ovarian follicles, controls the development and selection of dominant preovulatory follicles. E2 effects rely on a complex set of nuclear and extra-nuclear signal transduction pathways principally triggered by its nuclear receptors, ERα and ERβ. These transcription factors are differentially expressed within follicles, with ERβ being the predominant ER in GC. Several ERβ splice isoforms have been identified and display specific structural features, which greatly complicates the nature of ERβ-mediated E2 signaling. This review aims at providing a concise overview of the main actions of E2 during follicular growth, maturation, and selection in human. It also describes the current understanding of the various roles of ERβ splice isoforms, especially their influence on cell fate. We finally discuss how E2 signaling deregulation could participate in two ovarian pathogeneses characterized by either a follicular arrest, as in polycystic ovary syndrome, or an excess of GC survival and proliferation, leading to granulosa cell tumors. This review emphasizes the need for further research to better understand the molecular basis of E2 signaling throughout folliculogenesis and to improve the efficiency of ovarian-related disease therapies.

Transcript-Specific Loss-of-Function Variants in VPS16 Are Enriched in Patients With Dystonia

Background and ObjectivesOur objective was to improve rare variant interpretation using statistical measures as well as publicly accessible annotation of expression levels and tissue specificity of different splice isoforms. We describe rare VPS16 variants observed in patients with dystonia and patients without dystonia, elaborate on our interpretation of VPS16 variants affecting different transcripts, and provide detailed clinical description of the movement disorder caused by VPS16 variants.MethodsIn-house exome and genome data sets (n = 11,539) were screened for rare heterozygous missense and putative loss-of-function (pLoF) variants in VPS16. Using pext (proportion expressed across transcripts) values from the Genome Aggregation Database (gnomAD), we differentiated variants affecting weakly and highly expressed exons/transcripts and applied statistical measures to systematically identify disease-associated genetic variation among patients with dystonia (n = 280).ResultsSix different heterozygous pLoFs in VPS16 transcripts were identified in 13 individuals. Three of these pLoFs occurred in 9 individuals with different phenotypes, and 3 pLoFs were identified in 4 unrelated individuals with early-onset dystonia. Although pLoFs were enriched in the dystonia cohort (n = 280; p = 2.04 × 10−4; 4/280 cases vs 9/11,259 controls; Fisher exact test), it was not exome-wide significant. According to the pext values in gnomAD, all 3 pLoFs observed in the patients with dystonia were located in the highly expressed canonical transcript ENST00000380445.3, whereas 2 of 3 pLoFs detected in 8 individuals without dystonia were located in the first exon of the noncanonical transcript ENST00000380443.3 that is weakly expressed across all tissues. Taking these biological implications into account, pLoFs involving the canonical transcript were exome-wide significantly enriched in patients with dystonia (p = 1.67 × 10−6; 4/280 cases vs 1/11,259 controls; Fisher exact test). All VPS16 patients showed mild progressive dystonia with writer's cramp as the presenting symptom between age 7 and 34 years (mean 20 years) that often progressed to generalized dystonia and was even accompanied by hyperkinetic movements and myoclonus in 1 patient.DiscussionOur data provide strong evidence for VPS16 pLoFs to be implicated in dystonia and knowledge on exon resolution expression levels as well as statistical measures proved to be useful for variant interpretation.

Genome annotation with long RNA reads reveals new patterns of gene expression and improves single-cell analyses in an ant brain

Background Functional genomic analyses rely on high-quality genome assemblies and annotations. Highly contiguous genome assemblies have become available for a variety of species, but accurate and complete annotation of gene models, inclusive of alternative splice isoforms and transcription start and termination sites, remains difficult with traditional approaches. Results Here, we utilized full-length isoform sequencing (Iso-Seq), a long-read RNA sequencing technology, to obtain a comprehensive annotation of the transcriptome of the ant Harpegnathos saltator. The improved genome annotations include additional splice isoforms and extended 3′ untranslated regions for more than 4000 genes. Reanalysis of RNA-seq experiments using these annotations revealed several genes with caste-specific differential expression and tissue- or caste-specific splicing patterns that were missed in previous analyses. The extended 3′ untranslated regions afforded great improvements in the analysis of existing single-cell RNA-seq data, resulting in the recovery of the transcriptomes of 18% more cells. The deeper single-cell transcriptomes obtained with these new annotations allowed us to identify additional markers for several cell types in the ant brain, as well as genes differentially expressed across castes in specific cell types. Conclusions Our results demonstrate that Iso-Seq is an efficient and effective approach to improve genome annotations and maximize the amount of information that can be obtained from existing and future genomic datasets in Harpegnathos and other organisms.

Jacob, a Synapto-Nuclear Protein Messenger Linking N-methyl-D-aspartate Receptor Activation to Nuclear Gene Expression

Pyramidal neurons exhibit a complex dendritic tree that is decorated by a huge number of spine synapses receiving excitatory input. Synaptic signals not only act locally but are also conveyed to the nucleus of the postsynaptic neuron to regulate gene expression. This raises the question of how the spatio-temporal integration of synaptic inputs is accomplished at the genomic level and which molecular mechanisms are involved. Protein transport from synapse to nucleus has been shown in several studies and has the potential to encode synaptic signals at the site of origin and decode them in the nucleus. In this review, we summarize the knowledge about the properties of the synapto-nuclear messenger protein Jacob with special emphasis on a putative role in hippocampal neuronal plasticity. We will elaborate on the interactome of Jacob, the signals that control synapto-nuclear trafficking, the mechanisms of transport, and the potential nuclear function. In addition, we will address the organization of the Jacob/NSMF gene, its origin and we will summarize the evidence for the existence of splice isoforms and their expression pattern.

Reconstitution of human atlastin fusion activity reveals autoinhibition by the C terminus

ER network formation depends on membrane fusion by the atlastin (ATL) GTPase. In humans, three paralogs are differentially expressed with divergent N- and C-terminal extensions, but their respective roles remain unknown. This is partly because, unlike Drosophila ATL, the fusion activity of human ATLs has not been reconstituted. Here, we report successful reconstitution of fusion activity by the human ATLs. Unexpectedly, the major splice isoforms of ATL1 and ATL2 are each autoinhibited, albeit to differing degrees. For the more strongly inhibited ATL2, autoinhibition mapped to a C-terminal α-helix is predicted to be continuous with an amphipathic helix required for fusion. Charge reversal of residues in the inhibitory domain strongly activated its fusion activity, and overexpression of this disinhibited version caused ER collapse. Neurons express an ATL2 splice isoform whose sequence differs in the inhibitory domain, and this form showed full fusion activity. These findings reveal autoinhibition and alternate splicing as regulators of atlastin-mediated ER fusion.

A novel TREM2 splice isoform lacking the ligand binding is expressed in brain and shares localization

Background: Genetic variants in TREM2 are strongly associated with Alzheimer's Disease (AD) risk but alternative splicing in TREM2 transcripts has not been comprehensively described. Objective: Recognizing that alternative splice variants can result in reduced gene expression and/or altered function, we sought to fully characterize splice variation in TREM2. Methods: Human blood and anterior cingulate autopsy tissue from 61 donors were used for genotyping and cDNA synthesis followed by both end-point and quantitative PCR to identify and quantify novel TREM2 isoforms. Results: In addition to previously described transcripts lacking exon 3 or exon 4, or retaining part of intron 3, we identified novel isoforms lacking exon 2, along with isoforms lacking multiple exons. Isoforms lacking exon 2 were predominant at approximately 10% of TREM2 mRNA in the brain. Expression of TREM2 and frequency of exon 2 skipping did not differ between AD samples and non-AD controls (p = 0.1268 and p = 0.4909, respectively). Further, these novel splice isoforms were also observed across multiple tissues (brain, liver, lung, kidney, heart, aorta, skeletal muscle) with similar frequency (range 5.3 - 13.0%). Using ectopic expression, we found that the exon 2 skipped isoform D2-TREM2 is translated to protein and localizes similarly to full-length TREM2 protein, and that both D2-TREM2 and FL-TREM2 proteins are primarily retained in the Golgi complex. Conclusion: Since the TREM2 ligand binding domain is encoded by exon 2, and skipping this exon retains reading frame while conserving localization, we hypothesize that D2-TREM2 acts as an inhibitor of TREM2 and that targeting TREM2 splicing may be a novel therapeutic pathway for AD.

High-throughput mutagenesis identifies mutations and RNA binding proteins controlling CD19 splicing and CART-19 therapy resistance

During CART-19 immunotherapy for B-cell acute lymphoblastic leukaemia (B-ALL), many patients relapse due to loss of the cognate CD19 epitope. Since epitope loss can be caused by aberrant CD19 exon 2 processing, we herein investigate the regulatory code that controls CD19 splicing. We combine high-throughput mutagenesis with mathematical modelling to quantitatively disentangle the effects of all mutations in the region comprising CD19 exons 1-3. Thereupon, we identify ~200 single point mutations that alter CD19 splicing and thus could predispose B-ALL patients to CART-19 resistance. Furthermore, we report almost 100 previously unknown splice isoforms that emerge from cryptic splice sites and likely encode non-functional CD19 proteins. We further identify cis-regulatory elements and trans-acting RNA-binding proteins that control CD19 splicing (e.g., PTBP1 and SF3B4) and validate that loss of these factors leads to enhanced CD19 mis-splicing. Our dataset represents a comprehensive resource for potential prognostic factors predicting success of CART-19 therapy.

APPRIS principal isoforms and MANE Select transcripts in clinical variant interpretation

Most coding genes are able to generate multiple alternatively spliced transcripts. Determining which of these transcript variants produces the main protein isoform, and which of a gene’s multiple splice variants are functionally important, is crucial in comparative genomics and essential for clinical variant interpretation.Here we show that the principal isoforms chosen by APPRIS and the MANE Select variants provide the best approximations of the main cellular protein isoforms. Principal isoforms are predicted from conservation and from protein features, and MANE transcripts are chosen from the consensus between teams of expert manual curators. APPRIS principal isoforms coincide in over 94% of coding genes with MANE Select transcripts and the two methods are particularly discriminating when they agree on the main splice variant. Where the two methods agree, the splice variants coincide with the main isoform detected in proteomics experiments in 98.2% of genes with multiple protein isoforms.We also find that almost all ClinVar pathogenic mutations map to MANE Select or APPRIS principal isoforms. Where APPRIS and MANE agree on the main isoform, 99.93% of validated pathogenic variants map to principal rather than alternative exons. MANE Plus Clinical transcripts cover most validated pathogenic mutations in alternative coding exons. TRIFID functional importance scores are particularly useful for distinguishing clinically important alternative isoforms: the highest scoring TRIFID isoforms are more than 300 times more likely to have validated pathogenic mutations.We find that APPRIS, MANE and TRIFID are important for determining the biological relevance of splice isoforms and should be an essential part of clinical variant interpretation.

The Landscapes of Full-Length Transcripts and Splice Isoforms as Well as Transposons Exonization in the Lepidopteran Model System, Bombyx mori

The domesticated silkworm, Bombyx mori, is an important model system for the order Lepidoptera. Currently, based on third-generation sequencing, the chromosome-level genome of Bombyx mori has been released. However, its transcripts were mainly assembled by using short reads of second-generation sequencing and expressed sequence tags which cannot explain the transcript profile accurately. Here, we used PacBio Iso-Seq technology to investigate the transcripts from 45 developmental stages of Bombyx mori. We obtained 25,970 non-redundant high-quality consensus isoforms capturing ∼60% of previous reported RNAs, 15,431 (∼47%) novel transcripts, and identified 7,253 long non-coding RNA (lncRNA) with a large proportion of novel lncRNA (∼56%). In addition, we found that transposable elements (TEs) exonization account for 11,671 (∼45%) transcripts including 5,980 protein-coding transcripts (∼32%) and 5,691 lncRNAs (∼79%). Overall, our results expand the silkworm transcripts and have general implications to understand the interaction between TEs and their host genes. These transcripts resource will promote functional studies of genes and lncRNAs as well as TEs in the silkworm.