A dual-reference study design for understanding and improving AAV genome size analysis.

Recombinant adeno-associated virus (rAAV) is the leading platform of gene delivery for its long-lasting gene transformation and low immunogenicity. Characterization of the integrity and purity of the rAAV genome is critical to ensure clinical potency and safety. However, current rAAV genome characterization methods that can provide size assessment are either time-consuming or not easily accessible to general labs. Additionally, there is a lack of right reference standard for analyzing long single-stranded DNA (ssDNA) fragments. Here, we have developed an ssDNA assay on a microfluidic capillary electrophoresis platform using ssDNA reference standard. This assay provides size calling for ssDNA fragment, a detection sensitivity at ∼89 pg/µL (3 × 1010  GC/mL AAV) for 5.1 kb ssDNA fragment, and a turnaround time at ∼100 s per sample with a high throughput sample analyzing capability. Moreover, we have observed that the annealing of AAV ssDNA subsequent to its release from the capsid might introduce an additional double-stranded DNA (dsDNA) peak. This phenomenon is dependent on the sample processing workflow. To avoid the risk of mischaracterization, we recommend the use of dual-reference standards in combination with other orthogonal methods to have a comprehensive understanding of the rAAV genome size and integrity.

Endoscopic modified total laminoplasty for symptomatic lumbar spinal stenosis.

Context/objective: At present, there is no consensus on the most effective surgical method for treating symptomatic lumbar spinal stenosis (LSS). Total laminectomy, which is frequently used at this time, destroys the posterior midline structure, causing many postoperative complications. We have designed a new surgical approach instead of total laminectomy. In this paper, we aimed to describe the surgical method of endoscopic modified total laminectomy for lumbar spinal stenosis as well as to explore its early efficacy.Participants: Patients with symptomatic LSS who underwent endoscopic modified total laminoplasty between August 2016 and August 2017 were eligible for our study.Outcome measures: Before surgery and one year after surgery, we measured lower limb pain and back pain by visual analog scale (VAS), disability via Oswestry Disability Index (ODI), and severity of back pain according to the Japanese Orthopedic Association Score for Back Pain (JOA), while any complications were also assessed.Results: Endoscopic modified total laminoplasty was performed on 22 LSS patients, including eight males and 14 females(mean age = 59.3 ± 9.6 years). We found statistically significant differences before and one year after surgery for VAS lower limb pain and back pain, ODI and JOA scores(P < 0.001). Complications included intraoperative dural tears(n = 1),and weak fusion between the lamina and the vertebral body (n = 1).Conclusion: Endoscopic modified total laminectomy is a promising surgical approach which reduces patient suffering and improves patient quality of life.

[Impact of JAK2V617F Mutation Burden on Clinial Presentation and Survival in ET Patients].

OBJECTIVE: To better define the effect of JAK2V617F mutant allele burden on clinical presentation of patients with essential thrombo cythamia (ET), especially thrombosis. METHODS: Two ml of heparin anti-coagulated bone marrow was collected from 229 ET cases, who were diagnosed and treated in the First People's Hospital of Yunnan Province during 2013.10 to 2016.12. and then the mononuclear cells were separated by Red Blood Cell Lysis Buffer, genomic DNA was extracted from mononuclear cells by using a commercial DNA isolation kit and amplified by allele specific polymerase chain reaction (PCR). According to the size of molecular weight, the amplified products were separated by electrophoresis on a 2% agarose gel to screen the JAK2V617F mutation, then the JAK2V617F mutation burden was detected by real-time polymerase chain reaction (RT-PCR) in 120 patients with JAK2V617F mutation. Meanwhile, these samples were sequenced in order to verify the accuracy of the PCR screewing. RESULTS: ET patients with thrombotic events had significantly higher JAK2V617F allele burden than those without thrombosis (23.2% vs 14.2%) ( P<0.05). Meanwhile, ET patients showed increased JAK2V617F allele burden in the group with higher leukocytosis (WBC > 10×109/L) (P<0.001) and hemoglobin (> 150 g/L) (P<0.05). JAK2V617F mutation burden in 17 patients with splenomegaly was higher than that in 45 patients without splenomegaly (28.1% vs 11.8%) (P<0.05). but the JAK2V617F mutation burden was regatively correlated with platelet count (P<0.05). On the other hand, no correlation was found between JAK2V617F mutation burden and sex (P > 0.05). Univariate analysis showed that the JAK2V617F allele burden did not affect survival. Multivariable analysis showed that prognostic variable including WBC counts, hemoglobin level, age, sex, and splenomegaly not affected survival, (P > 0.05). CONCLUSION: The clinical presentations of ET patients, such as WBC counts, hemoglobin level and splenomegaly, are influenced by the JAK2V617F mutation burden. ET patients with thrombotic events has significantly higher JAK2V617F allele burden than those in ET palients without thrombosis.JAK2V617F mutation burden has no relations with sex and age..

Concerted effects of heterogeneous nuclear ribonucleoprotein C1/C2 to control vitamin D-directed gene transcription and RNA splicing in human bone cells.

Traditionally recognized as an RNA splicing regulator, heterogeneous nuclear ribonucleoprotein C1/C2 (hnRNPC1/C2) can also bind to double-stranded DNA and function in trans as a vitamin D response element (VDRE)-binding protein. As such, hnRNPC1/C2 may couple transcription induced by the active form of vitamin D, 1,25-dihydroxyvitamin D (1,25(OH)2D) with subsequent RNA splicing. In MG63 osteoblastic cells, increased expression of the 1,25(OH)2D target gene CYP24A1 involved immunoprecipitation of hnRNPC1/C2 with CYP24A1 chromatin and RNA. Knockdown of hnRNPC1/C2 suppressed expression of CYP24A1, but also increased expression of an exon 10-skipped CYP24A1 splice variant; in a minigene model the latter was attenuated by a functional VDRE in the CYP24A1 promoter. In genome-wide analyses, knockdown of hnRNPC1/C2 resulted in 3500 differentially expressed genes and 2232 differentially spliced genes, with significant commonality between groups. 1,25(OH)2D induced 324 differentially expressed genes, with 187 also observed following hnRNPC1/C2 knockdown, and a further 168 unique to hnRNPC1/C2 knockdown. However, 1,25(OH)2D induced only 10 differentially spliced genes, with no overlap with differentially expressed genes. These data indicate that hnRNPC1/C2 binds to both DNA and RNA and influences both gene expression and RNA splicing, but these actions do not appear to be linked through 1,25(OH)2D-mediated induction of transcription.

m(6)A-LAIC-seq reveals the census and complexity of the m(6)A epitranscriptome.

N(6)-Methyladenosine (m(6)A) is a widespread, reversible chemical modification of RNA molecules, implicated in many aspects of RNA metabolism. Little quantitative information exists as to either how many transcript copies of particular genes are m(6)A modified ('m(6)A levels') or the relationship of m(6)A modification(s) to alternative RNA isoforms. To deconvolute the m(6)A epitranscriptome, we developed m(6)A-level and isoform-characterization sequencing (m(6)A-LAIC-seq). We found that cells exhibit a broad range of nonstoichiometric m(6)A levels with cell-type specificity. At the level of isoform characterization, we discovered widespread differences in the use of tandem alternative polyadenylation (APA) sites by methylated and nonmethylated transcript isoforms of individual genes. Strikingly, there is a strong bias for methylated transcripts to be coupled with proximal APA sites, resulting in shortened 3' untranslated regions, while nonmethylated transcript isoforms tend to use distal APA sites. m(6)A-LAIC-seq yields a new perspective on transcriptome complexity and links APA usage to m(6)A modifications.

The contribution of Alu exons to the human proteome.

BACKGROUND: Alu elements are major contributors to lineage-specific new exons in primate and human genomes. Recent studies indicate that some Alu exons have high transcript inclusion levels or tissue-specific splicing profiles, and may play important regulatory roles in modulating mRNA degradation or translational efficiency. However, the contribution of Alu exons to the human proteome remains unclear and controversial. The prevailing view is that exons derived from young repetitive elements, such as Alu elements, are restricted to regulatory functions and have not had adequate evolutionary time to be incorporated into stable, functional proteins. RESULTS: We adopt a proteotranscriptomics approach to systematically assess the contribution of Alu exons to the human proteome. Using RNA sequencing, ribosome profiling, and proteomics data from human tissues and cell lines, we provide evidence for the translational activities of Alu exons and the presence of Alu exon derived peptides in human proteins. These Alu exon peptides represent species-specific protein differences between primates and other mammals, and in certain instances between humans and closely related primates. In the case of the RNA editing enzyme ADARB1, which contains an Alu exon peptide in its catalytic domain, RNA sequencing analyses of A-to-I editing demonstrate that both the Alu exon skipping and inclusion isoforms encode active enzymes. The Alu exon derived peptide may fine tune the overall editing activity and, in limited cases, the site selectivity of ADARB1 protein products. CONCLUSIONS: Our data indicate that Alu elements have contributed to the acquisition of novel protein sequences during primate and human evolution.

Transcriptome-wide landscape of pre-mRNA alternative splicing associated with metastatic colonization.

UNLABELLED: Metastatic colonization is an ominous feature of cancer progression. Recent studies have established the importance of pre-mRNA alternative splicing (AS) in cancer biology. However, little is known about the transcriptome-wide landscape of AS associated with metastatic colonization. Both in vitro and in vivo models of metastatic colonization were utilized to study AS regulation associated with cancer metastasis. Transcriptome profiling of prostate cancer cells and derivatives crossing in vitro or in vivo barriers of metastasis revealed splicing factors with significant gene expression changes associated with metastatic colonization. These include splicing factors known to be differentially regulated in epithelial-mesenchymal transition (ESRP1, ESRP2, and RBFOX2), a cellular process critical for cancer metastasis, as well as novel findings (NOVA1 and MBNL3). Finally, RNA-seq indicated a large network of AS events regulated by multiple splicing factors with altered gene expression or protein activity. These AS events are enriched for pathways important for cell motility and signaling, and affect key regulators of the invasive phenotype such as CD44 and GRHL1. IMPLICATIONS: Transcriptome-wide remodeling of AS is an integral regulatory process underlying metastatic colonization, and AS events affect the metastatic behavior of cancer cells.

rMATS: robust and flexible detection of differential alternative splicing from replicate RNA-Seq data.

Ultra-deep RNA sequencing (RNA-Seq) has become a powerful approach for genome-wide analysis of pre-mRNA alternative splicing. We previously developed multivariate analysis of transcript splicing (MATS), a statistical method for detecting differential alternative splicing between two RNA-Seq samples. Here we describe a new statistical model and computer program, replicate MATS (rMATS), designed for detection of differential alternative splicing from replicate RNA-Seq data. rMATS uses a hierarchical model to simultaneously account for sampling uncertainty in individual replicates and variability among replicates. In addition to the analysis of unpaired replicates, rMATS also includes a model specifically designed for paired replicates between sample groups. The hypothesis-testing framework of rMATS is flexible and can assess the statistical significance over any user-defined magnitude of splicing change. The performance of rMATS is evaluated by the analysis of simulated and real RNA-Seq data. rMATS outperformed two existing methods for replicate RNA-Seq data in all simulation settings, and RT-PCR yielded a high validation rate (94%) in an RNA-Seq dataset of prostate cancer cell lines. Our data also provide guiding principles for designing RNA-Seq studies of alternative splicing. We demonstrate that it is essential to incorporate biological replicates in the study design. Of note, pooling RNAs or merging RNA-Seq data from multiple replicates is not an effective approach to account for variability, and the result is particularly sensitive to outliers. The rMATS source code is freely available at rnaseq-mats.sourceforge.net/. As the popularity of RNA-Seq continues to grow, we expect rMATS will be useful for studies of alternative splicing in diverse RNA-Seq projects.

MATS: a Bayesian framework for flexible detection of differential alternative splicing from RNA-Seq data.

Ultra-deep RNA sequencing has become a powerful approach for genome-wide analysis of pre-mRNA alternative splicing. We develop MATS (multivariate analysis of transcript splicing), a bayesian statistical framework for flexible hypothesis testing of differential alternative splicing patterns on RNA-Seq data. MATS uses a multivariate uniform prior to model the between-sample correlation in exon splicing patterns, and a Markov chain Monte Carlo (MCMC) method coupled with a simulation-based adaptive sampling procedure to calculate the P-value and false discovery rate (FDR) of differential alternative splicing. Importantly, the MATS approach is applicable to almost any type of null hypotheses of interest, providing the flexibility to identify differential alternative splicing events that match a given user-defined pattern. We evaluated the performance of MATS using simulated and real RNA-Seq data sets. In the RNA-Seq analysis of alternative splicing events regulated by the epithelial-specific splicing factor ESRP1, we obtained a high RT-PCR validation rate of 86% for differential exon skipping events with a MATS FDR of <10%. Additionally, over the full list of RT-PCR tested exons, the MATS FDR estimates matched well with the experimental validation rate. Our results demonstrate that MATS is an effective and flexible approach for detecting differential alternative splicing from RNA-Seq data.

Context-dependent robustness to 5' splice site polymorphisms in human populations.

There has been growing evidence for extensive diversity of alternative splicing in human populations. Genetic variants within the 5' splice site can cause splicing differences among human individuals and constitute an important class of human disease mutations. In this study, we explored whether natural variations of splicing could reveal important signals of 5' splice site recognition. In seven lymphoblastoid cell lines of Asian, European and African ancestry, we identified 1174 single nucleotide polymorphisms (SNPs) within the consensus 5' splice site. We selected 129 SNPs predicted to significantly alter the splice site activity, and quantitatively examined their splicing impact in the seven individuals. Surprisingly, outside of the essential GT dinucleotide position, only ∼14% of the tested SNPs altered splicing. Bioinformatic and minigene analyses identified signals that could modify the impact of 5' splice site polymorphisms, most notably a strong 3' splice site and the presence of intronic motifs downstream of the 5' splice site. Strikingly, we found that the poly-G run, a known intronic splicing enhancer, was the most significantly enriched motif downstream of exons unaffected by 5' splice site SNPs. In TRIM62, the upstream 3' splice site and downstream intronic poly-G runs functioned redundantly to protect an exon from its 5' splice site polymorphism. Collectively, our study reveals widespread context-dependent robustness to 5' splice site polymorphisms in human transcriptomes. Consequently, certain exons are more susceptible to 5' splice site mutations. Additionally, our work demonstrates that genetic diversity of alternative splicing can provide significant insights into the splicing code of mammalian cells.