QuickMIRSeq: a pipeline for quick and accurate quantification of both known miRNAs and isomiRs by jointly processing multiple samples from microRNA sequencing.

BACKGROUND: Genome-wide miRNA expression data can be used to study miRNA dysregulation comprehensively. Although many open-source tools for microRNA (miRNA)-seq data analyses are available, challenges remain in accurate miRNA quantification from large-scale miRNA-seq dataset. We implemented a pipeline called QuickMIRSeq for accurate quantification of known miRNAs and miRNA isoforms (isomiRs) from multiple samples simultaneously. RESULTS: QuickMIRSeq considers the unique nature of miRNAs and combines many important features into its implementation. First, it takes advantage of high redundancy of miRNA reads and introduces joint mapping of multiple samples to reduce computational time. Second, it incorporates the strand information in the alignment step for more accurate quantification. Third, reads potentially arising from background noise are filtered out to improve the reliability of miRNA detection. Fourth, sequences aligned to miRNAs with mismatches are remapped to a reference genome to further reduce false positives. Finally, QuickMIRSeq generates a rich set of QC metrics and publication-ready plots. CONCLUSIONS: The rich visualization features implemented allow end users to interactively explore the results and gain more insights into miRNA-seq data analyses. The high degree of automation and interactivity in QuickMIRSeq leads to a substantial reduction in the time and effort required for miRNA-seq data analysis.

Modeling the clinical phenotype of BTK inhibition in the mature murine immune system.

Inhibitors of Bruton's tyrosine kinase (BTK) possess much promise for the treatment of oncologic and autoimmune indications. However, our current knowledge of the role of BTK in immune competence has been gathered in the context of genetic inactivation of btk in both mice and man. Using the novel BTK inhibitor PF-303, we model the clinical phenotype of BTK inhibition by systematically examining the impact of PF-303 on the mature immune system in mice. We implicate BTK in tonic BCR signaling, demonstrate dependence of the T3 B cell subset and IgM surface expression on BTK activity, and find that B1 cells survive and function independently of BTK. Although BTK inhibition does not impact humoral memory survival, Ag-driven clonal expansion of memory B cells and Ab-secreting cell generation are inhibited. These data define the role of BTK in the mature immune system and mechanistically predict the clinical phenotype of chronic BTK inhibition.

Comparison of stranded and non-stranded RNA-seq transcriptome profiling and investigation of gene overlap.

BACKGROUND: While RNA-sequencing (RNA-seq) is becoming a powerful technology in transcriptome profiling, one significant shortcoming of the first-generation RNA-seq protocol is that it does not retain the strand specificity of origin for each transcript. Without strand information it is difficult and sometimes impossible to accurately quantify gene expression levels for genes with overlapping genomic loci that are transcribed from opposite strands. It has recently become possible to retain the strand information by modifying the RNA-seq protocol, known as strand-specific or stranded RNA-seq. Here, we evaluated the advantages of stranded RNA-seq in transcriptome profiling of whole blood RNA samples compared with non-stranded RNA-seq, and investigated the influence of gene overlaps on gene expression profiling results based on practical RNA-seq datasets and also from a theoretical perspective. RESULTS: Our results demonstrated a substantial impact of stranded RNA-seq on transcriptome profiling and gene expression measurements. As many as 1751 genes in Gencode Release 19 were identified to be differentially expressed when comparing stranded and non-stranded RNA-seq whole blood samples. Antisense and pseudogenes were significantly enriched in differential expression analyses. Because stranded RNA-seq retains strand information of a read, we can resolve read ambiguity in overlapping genes transcribed from opposite strands, which provides a more accurate quantification of gene expression levels compared with traditional non-stranded RNA-seq. In the human genome, it is not uncommon to find genomic loci where both strands encode distinct genes. Among the over 57,800 annotated genes in Gencode release 19, there are an estimated 19 % (about 11,000) of overlapping genes transcribed from the opposite strands. Based on our whole blood mRNA-seq datasets, the fraction of overlapping nucleotide bases on the same and opposite strands were estimated at 2.94 % and 3.1 %, respectively. The corresponding theoretical estimations are 3 % and 3.6 %, well in agreement with our own findings. CONCLUSIONS: Stranded RNA-seq provides a more accurate estimate of transcript expression compared with non-stranded RNA-seq, and is therefore the recommended RNA-seq approach for future mRNA-seq studies.

An enhancer-AAV toolbox to target and manipulate distinct interneuron subtypes.

In recent years, we and others have identified a number of enhancers that, when incorporated into rAAV vectors, can restrict the transgene expression to particular neuronal populations. Yet, viral tools to access and manipulate fine neuronal subtypes are still limited. Here, we performed systematic analysis of single cell genomic data to identify enhancer candidates for each of the cortical interneuron subtypes. We established a set of enhancer-AAV tools that are highly specific for distinct cortical interneuron populations and striatal cholinergic neurons. These enhancers, when used in the context of different effectors, can target (fluorescent proteins), observe activity (GCaMP) and manipulate (opto- or chemo-genetics) specific neuronal subtypes. We also validated our enhancer-AAV tools across species. Thus, we provide the field with a powerful set of tools to study neural circuits and functions and to develop precise and targeted therapy.

Systemic Copper Disorders Influence the Olfactory Function in Adult Rats: Roles of Altered Adult Neurogenesis and Neurochemical Imbalance.

Disrupted systemic copper (Cu) homeostasis underlies neurodegenerative diseases with early symptoms including olfactory dysfunction. This study investigated the impact of Cu dyshomeostasis on olfactory function, adult neurogenesis, and neurochemical balance. Models of Cu deficiency (CuD) and Cu overload (CuO) were established by feeding adult rats with Cu-restricted diets plus ip. injection of a Cu chelator (ammonium tetrathiomolybdate) and excess Cu, respectively. CuD reduced Cu levels in the olfactory bulb (OB), subventricular zone (SVZ), rostral migratory stream (RMS), and striatum, while CuO increased Cu levels in these areas. The buried pellet test revealed both CuD and CuO prolonged the latency to uncover food. CuD increased neural proliferation and stem cells in the SVZ and newly differentiated neurons in the OB, whereas CuO caused opposite alterations, suggesting a "switch"-type function of Cu in regulating adult neurogenesis. CuO increased GABA in the OB, while both CuD and CuO reduced DOPAC, HVA, 5-HT and the DA turnover rate in olfactory-associated brain regions. Altered mRNA expression of Cu transport and storage proteins in tested brain areas were observed under both conditions. Together, results support an association between systemic Cu dyshomeostasis and olfactory dysfunction. Specifically, altered adult neurogenesis along the SVZ-RMS-OB pathway and neurochemical imbalance could be the factors that may contribute to olfactory dysfunction.

Mucosal Gene Expression in Pediatric and Adult Patients With Ulcerative Colitis Permits Modeling of Ideal Biopsy Collection Strategy for Transcriptomic Analysis.

Background: Transcriptional profiling has been performed on biopsies from ulcerative colitis patients. Limitations in prior studies include the variability introduced by inflammation, anatomic site of biopsy, extent of disease, and medications. We sought to more globally understand the variability of gene expression from patients with ulcerative colitis to advance our understanding of its pathogenesis and to guide clinical study design. Methods: We performed transcriptional profiling on 13 subjects, including pediatric and adult patients from 2 hospital sites. For each patient, we collected 6 biopsies from macroscopically inflamed tissue and 4 biopsies from macroscopically healthy-appearing tissue. Isolated RNA was used for microarray gene expression analysis utilizing Affymetrix Human Primeview microarrays. Ingenuity pathway analysis was used to assess over-representation of gene ontology and biological pathways. RNAseq was also performed, and differential analysis was assessed to compare affected vs unaffected samples. Finally, we modeled the minimum number of biopsies required to reliably detect gene expression across different subject numbers. Results: Transcriptional profiles co-clustered independently of the hospital collection site, patient age, sex, and colonic location, which parallels prior gene expression findings. A small set of genes not previously described was identified. Our modeling analysis reveals the number of biopsies and patients per cohort to yield reliable results in clinical studies. Conclusions: Key findings include concordance, including some expansion, of previously published gene expression studies and similarity among different age groups. We also established a reliable statistical model for biopsy collection for future clinical studies.

In vitro and in vivo effects of 2,4 diaminoquinazoline inhibitors of the decapping scavenger enzyme DcpS: Context-specific modulation of SMN transcript levels.

C5-substituted 2,4-diaminoquinazoline inhibitors of the decapping scavenger enzyme DcpS (DAQ-DcpSi) have been developed for the treatment of spinal muscular atrophy (SMA), which is caused by genetic deficiency in the Survival Motor Neuron (SMN) protein. These compounds are claimed to act as SMN2 transcriptional activators but data underlying that claim are equivocal. In addition it is unclear whether the claimed effects on SMN2 are a direct consequence of DcpS inhibitor or might be a consequence of lysosomotropism, which is known to be neuroprotective. DAQ-DcpSi effects were characterized in cells in vitro utilizing DcpS knockdown and 7-methyl analogues as probes for DcpS vs non-DcpS-mediated effects. We also performed analysis of Smn transcript levels, RNA-Seq analysis of the transcriptome and SMN protein in order to identify affected pathways underlying the therapeutic effect, and studied lysosomotropic and non-lysosomotropic DAQ-DCpSi effects in 2B/- SMA mice. Treatment of cells caused modest and transient SMN2 mRNA increases with either no change or a decrease in SMNΔ7 and no change in SMN1 transcripts or SMN protein. RNA-Seq analysis of DAQ-DcpSi-treated N2a cells revealed significant changes in expression (both up and down) of approximately 2,000 genes across a broad range of pathways. Treatment of 2B/- SMA mice with both lysomotropic and non-lysosomotropic DAQ-DcpSi compounds had similar effects on disease phenotype indicating that the therapeutic mechanism of action is not a consequence of lysosomotropism. In striking contrast to the findings in vitro, Smn transcripts were robustly changed in tissues but there was no increase in SMN protein levels in spinal cord. We conclude that DAQ-DcpSi have reproducible benefit in SMA mice and a broad spectrum of biological effects in vitro and in vivo, but these are complex, context specific, and not the result of simple SMN2 transcriptional activation.

The ligand-dependent interaction of mineralocorticoid receptor with coactivator and corepressor peptides suggests multiple activation mechanisms.

We investigated the coregulator (coactivator and corepressor) interactions with the mineralocorticoid receptor (MR) that lead to activation and inhibition of the receptor in the presence of agonist and/or antagonist. Our results indicate that MR ligand binding domain (LBD) interacts strongly with only a few specific coactivator peptides in the presence of the agonist aldosterone and that these interactions are blocked by the antagonist eplerenone. We also discovered that cortisol, the preferred physiological ligand for the glucocorticoid receptor in humans, is a partial MR agonist/antagonist, providing a possible molecular explanation of the tissue-selective effects of glucocorticoids on MR. However, when we examined the coactivator and corepressor peptide interactions in the presence of cortisol, we found that MR bound with cortisol or aldosterone interacted with the same set of peptides. Thus, the partial agonism shown by cortisol is unlikely to be the result of differential interaction with known coactivators and corepressors. On the other hand, we have identified coactivator binding groove mutations that are critical for cortisol activation but not for aldosterone activation, suggesting that the two steroids induce different MR LBD conformations. In addition, we also show that cortisol becomes full agonist when S810L mutation is introduced in the LBD of MR. Interestingly, MR antagonists, such as eplerenone and progesterone, become partial agonist/antagonist of S810L but are still able to recruit LXXLL peptides to the mutant receptor. Together, these findings suggest a model to explain the MR activation by various ligands.

The antagonists but not partial agonists of glucocorticoid receptor ligands show substantial side effect dissociation.

A synthetic glucocorticoid receptor (GR) ligand with the efficacy of a glucocorticoid, but without the accompanying side effects, would meet an unmet medical need for the treatment of inflammatory diseases. It was hypothesized that a GR ligand that shifted helix 12 in a manner distinct from an agonist and an antagonist would confer a distinct GR conformation, resulting in differential gene expression and, ultimately, dissociation of antiinflammatory activity from side effects. A structural feature expected to interfere with helix 12 was incorporated into a nonsteroidal, tricyclic scaffold to create novel, high-affinity, and selective GR ligands that manifested a dual function in cellular assays, partial but robust agonist activity for inflammatory cytokine inhibition, and full antagonist activity for reporter gene activation. In contrast, analogs not likely to hinder helix 12 exhibited partial agonist activity for reporter gene activation. The requirement of full antagonist activity for substantial side effect dissociation was demonstrated in primary human preadipocytes, hepatocytes, and osteoblasts in which effects on adipogenesis, key genes involved in gluconeogenesis, and genes important for bone formation were examined, respectively. The dissociated GR ligands, despite lacking significant reporter gene activation, weakly recruit a limited number of coactivators such as peroxisomal proliferator-activated receptor-γ coactivator 1α. Transcriptional activation was sensitive to both peroxisomal proliferator-activated receptor-γ coactivator 1α and GR levels, providing a basis for cell-selective modulation of gene expression. The antiinflammatory activity of the dissociated ligands was further demonstrated in mouse models of inflammation. Together these results suggest that these ligands are promising candidates with robust antiinflammatory activity and likely dissociation against glucocorticoid-induced side effects.

A number of marketed dihydropyridine calcium channel blockers have mineralocorticoid receptor antagonist activity.

Calcium channel blockers are widely used antihypertensives. Mineralocorticoid receptor antagonists are also used to treat hypertension and heart failure. We report here that a number of widely used dihydropyridine class calcium channel blockers are able to inhibit aldosterone-induced activation of mineralocorticoid receptor. These dihydropyridines varied in the extent of their effect on mineralocorticoid receptor, with nimodipine and felodipine the most potent and amlodipine the least. In contrast, both diltiazem and verapamil, nondihydropyridine calcium channel blockers, had no effect on mineralocorticoid receptor. These dihydropyridines compete with aldosterone for binding and block aldosterone-induced coactivator recruitment to mineralocorticoid receptor. The mineralocorticoid receptor S810L mutant, which is activated by steroidal mineralocorticoid receptor antagonist such as eplerenone, is inhibited by these drugs. Furthermore, nimodipine decreased aldosterone-induced expression of the mineralocorticoid receptor target gene epithelial sodium channel gamma subunit in adrenalectomized rats, demonstrating that dihydropyridine calcium channel blockers can function as mineralocorticoid receptor antagonists in vivo. Molecular modeling indicates that dihydropyridines dock into the ligand binding domain of mineralocorticoid receptor in a consensus pose that partially overlaps with steroidal mineralocorticoid receptor antagonists. Together, our data suggest that, in addition to their calcium channel blocking activity, a number of dihydropyridine calcium channel blockers also have mineralocorticoid receptor antagonist activity at high doses, a finding which may thus prove useful for the design of novel antihypertensive drugs in the future.

Optimization and validation of a quantitative capillary electrophoresis sodium dodecyl sulfate method for quality control and stability monitoring of monoclonal antibodies.

In previous work, a capillary electrophoresis sodium dodecyl sulfate (CE-SDS) method using precolumn labeling and laser-induced fluorescence (LIF) detection was developed at Genentech Inc. as part of the control system for the quality control release of a recombinant monoclonal antibody (rMAb) (Hunt, G.; Nashabeh, W. Anal. Chem. 1999, 71, 2390-2397.). In the current work, a generic and quantitative CE-SDS assay with LIF detection of rMAbs with improved accuracy and precision is described. The implementation of an alkylating step with iodoacetamide and optimization of the incubation temperature and time, in the presence of SDS, greatly decrease any thermally induced fragmentation of nonreduced labeled rMAb samples. In addition, a quantitative study of the effects of sample buffer pH on rMAb fragmentation is also presented. Furthermore, the performance of alternative CE-SDS polymer solutions and instrumentation for quantitative analysis of rMAbs is shown in this article. The validation of this method, under the guidelines of the International Committee on Harmonization (ICH), demonstrates that the assay quantitatively determines the consistency of rMAb manufacture as it relates to size heterogeneity and product purity.